Найти точки экстремума sin(x) на заданном интервале

@KrotZaDroT · Регистрация: 18.04.2014

Студворк — интернет-сервис помощи студентам

переделал программу из нахождения экстремумов параболы в sin(x). выдает ошибку "Run-Time Check Failure #3 - The variable 'b' is being used without being initialized." Помогите доделать)

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
#include <iostream>
#include <conio.h>
#include <cmath>
#include <locale.h>
using namespace std;
 
double f ( double x)
{
    
        return sin (x);
}
double extremum ( double a, double b)
{
    float x; 
    double eps=0.001;
    //int k;
    //if (A>=0) k=1; else k=-1;
    while (abs(b-a)>eps)
    {
        x=(a+b)/2;
        if (f(x-eps)<f(x+eps))
            b=x;
        else 
            a=x;
    }
    return (a+b)/2;
}
 
 
void main() 
{
    setlocale(LC_ALL,"Rus");
    float a,b,y,z,x;
    cout<<"Введите промежуток на котором искать экстремум[a,b]:"<<endl;
    cout<<"a: "; cin>>y;
    cout<<"b: "; cin>>z;
    x=extremum(a,b);
    cout<<"\nextr x:"<<x<<endl;
    cout<<"f(extr x): "<<f(x)<<endl;
    getch();
}

@SuurKissat · 21.01.2015, 17:10

1. Ошибка означает, что переменная b используется, хотя не инициализирована.
2. Вижу в коде строки

C++
1
2
3
4
5
    float a,b,y,z,x;
    cout<<"Введите промежуток на котором искать экстремум[a,b]:"<<endl;
    cout<<"a: "; cin>>y;
    cout<<"b: "; cin>>z;
    x=extremum(a,b);

Вы почему-то предлагаете ввести a и b, а код как написан - вводит это в x и y. a и b остались неинициализированными, и понятно следующая строка x=extremum(a,b); дает ошибку.

@ilbm · 21.01.2015, 17:21

то есть так

Сообщение от SuurKissat

x=extremum(y,z);

@KrotZaDroT · 21.01.2015, 17:23 **[ТС]**

Спасибо )сам не заметил)

Добавлено через 1 минуту
SuurKissat,
это мы наши минимум, а для максимума надо в условии поменять знак, но как все это записать, что бы не запутаться?

C++
1
 if (f(x-eps)<f(x+eps))

ввести переменную x2 как максимум?

@The535 · 21.01.2015, 20:14

Почитайте о методе золотого сечения. Если найду, скину программу, я как-то ее писал.

Добавлено через 11 минут
Нашел таки программу.
В программе учитывается шаг сетки. Программа работает для любой функции. Вы задаете любую точку, программа ищет самый близкий экстремум.

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
#include <stdio.h>
#include <math.h>
 
double f (double x);
double otrezok (double a);
double minimum (double a, double b);
double waq (double a);
 
int main (void) {
    double a, b, x;
    printf("Vvedite A: ");
    scanf("%lg", &a);
    x = otrezok(a);
    printf("\nX = %lg\n", x);
    return 0;
}
 
double otrezok (double a) {
    double b, c, t, d, k=(1 + sqrt(5.0))/2;
    int i=0;
    if(a > 1)    d = 1e6 * waq(a);
    else         d = k;
    b = a + d;
    if(f(b) > f(a)) {
        d = -d;
        c = a;
        a = b;
        b = c;
    }
    while(f(b) < f(a)) {
        t = a;
        a = b;
        b += d;
        d *= k;
        i++;
    }
 
    c = minimum(b, t);
    return c;
}
 
double minimum (double a, double b) {
    double c, d, fa, fb, fc, fd, eps = 2*1e-15, x, k = (3.0 - sqrt(5.0))/2.0;
    if(b<a){c=b; b=a; a=c;}
    fa = f(a);
    fb = f(b);
    c = a + k * (b - a);
    d = b - k * (b - a);
    fc = f(c);
    fd = f(d);
    while((fabs(b - a) > eps) && (fabs(b - a) > eps * fabs((a + b)/2)) && (d > c)) {
        if(fd < fc) {
            a = c;
            fa = fc;
            c = d;
            d = b -  (b - a) * k;
            if (c>d) {
                c = a + k * (b - a);
                d = b - k * (b - a);
                fc = f(c);
                fd = f(d);
            }
            else { fc = fd;        fd = f(d);}
 
        }
        else
        {
            b = d;
            fb = fd;
            d = c;
            c = a +  (b - a) * k;
            if (c>d) {
                c = a + k * (b - a);
                d = b - k * (b - a);
                fc = f(c);
                fd = f(d);
            }
            else {  fd = fc; fc = f(c); }
            if (c>d) return 0;
        }
    }
 
    return (a + b)/2;
}
 
double f (double x) {
    return sin(x);
}
 
double waq (double a) {
    double b, c, d;
    b = a;
    while(a + b > a) {
        c = a + b;
        d = c - a;
        b /= 2;
    }
    return d;
}

П.С. Где вы видели параболу sin(x)?)) Может синусоида?)

Добавлено через 6 минут
Функция waq - считает шаг сетки, т.е. машинное эпсилон. Разница между двумя соседними машинно представимыми числами.

Новые блоги и статьи Все статьи Все блоги /
Символьное дифференцирование igorrr37 13.02.2026 / * Программа принимает математическое выражение в виде строки и выдаёт его производную в виде строки и вычисляет значение производной при заданном х Логарифм записывается как: (x-2)log(x^2+2) -. . .	Камера Toupcam IUA500KMA Eddy_Em 12.02.2026 Т. к. у всяких "хикроботов" слишком уж мелкий пиксель, для подсмотра в ESPriF они вообще плохо годятся: уже 14 величину можно рассмотреть еле-еле лишь на экспозициях под 3 секунды (а то и больше),. . .	И ясному Солнцу zbw 12.02.2026 И ясному Солнцу, и светлой Луне. В мире покоя нет и люди не могут жить в тишине. А жить им немного лет.	«Знание-Сила» zbw 12.02.2026 «Знание-Сила» «Время-Деньги» «Деньги -Пуля»
SDL3 для Web (WebAssembly): Подключение Box2D v3, физика и отрисовка коллайдеров 8Observer8 12.02.2026 Содержание блога Box2D - это библиотека для 2D физики для анимаций и игр. С её помощью можно определять были ли коллизии между конкретными объектами и вызывать обработчики событий столкновения. . . .	SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL_LoadPNG (без SDL3_image) 8Observer8 11.02.2026 Содержание блога Библиотека SDL3 содержит встроенные инструменты для базовой работы с изображениями - без использования библиотеки SDL3_image. Пошагово создадим проект для загрузки изображения. . .	SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL3_image 8Observer8 10.02.2026 Содержание блога Библиотека SDL3_image содержит инструменты для расширенной работы с изображениями. Пошагово создадим проект для загрузки изображения формата PNG с альфа-каналом (с прозрачным. . .	Установка Qt-версии Lazarus IDE в Debian Trixie Xfce volvo 10.02.2026 В общем, достали меня глюки IDE Лазаруса, собранной с использованием набора виджетов Gtk2 (конкретно: если набирать текст в редакторе и вызвать подсказку через Ctrl+Space, то после закрытия окошка. . .