Модифицированный метод наилучшей пробы

@xSimplex · Регистрация: 24.07.2013

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

Доброго времени суток.
Вопросы по моей реализации:
1) Есть ли хоть доля правды в том, что я написал.(правильно ли я понял математику

)
2) Прошу помочь вытащить необходимый мне \xi.
P.S. когда смогу вытащить необходимый \xi, тогда буду менять аргументы фу-ии

.

Собственно, то что я написал

.

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
100
#include <iostream>
#include <conio.h>
#include <iomanip>
 
using namespace std;
 
double f(double x1, double x2){
    return 8 * (x1 - 1)*(x1 - 1) + 2 * (x2 - 4)*(x2 - 4);
}
 
//ограничения
double g1(double x1, double x2){
    return x1 + 2 * x2 - 12;
}
double g2(double x1, double x2){
    return -x1 + 2;
}
double g3(double x1, double x2){
    return -x2 - 1;
}
 
//штрафы в виде gi(x1,x2)*gi(x1,x2) , i=1,2,3
double sg1(double x1, double x2){
    return (g1(x1, x2) > 0) ? g1(x1, x2)*g1(x1, x2) : 0;
}
 
double sg2(double x1, double x2){
    return (g2(x1, x2) > 0) ? g2(x1, x2)*g2(x1, x2) : 0;
}
 
double sg3(double x1, double x2){
    return (g3(x1, x2) > 0) ? g3(x1, x2)*g3(x1, x2) : 0;
}
 
//вспомогательная функция
double F(double r, double x1, double x2){
    return f(x1, x2) + r*(sg1(x1, x2) + sg2(x1, x2) + sg3(x1, x2));
}
 
 
 
int main()
{
    setlocale(LC_ALL, "Russian");
 
    double x1 = .3, x2 = .1, r = 1.;
 
///////////////////////////////////////////////////////////////////////////////////////////////////////////////
/*                            Модифицированный метод наилучшей пробы                                         */
                                                                                                             
                                                                                                                            
    ///Исходные данные////                                                                                     
    const int n = 10;   //                                                                                   
    double h = .1;      //                                                                                   
    //////////////////////                                                                                   
                                                                                                             
    double psi[n], psiw[n], /*args[2],*/ y[n];                                                                   
    double min_psi[n];                                                                                           
                                                                                                             
    for (int k = 0; k < 5; k++)
    {
        cout << "Итерация " << k << "\n\n";
 
        for (int i = 0; i < n; i++)                                                                              
        {                                                                                                        
            psi[i] = ((rand() % 100 / 100. > .5) ? (1) : (-1))*(rand() % 100) / 100.;                            
            cout << "psi  "<< setw(6) << psi[i] << endl;                                                                     
        }                                                                                                        
        
        cout << endl;                                                                                            
        
        for (int j = 0; j < n; j++)                                                                              
        {                                                                                                        
            if (F(r, x1 + h*psi[j], x2 + h*psi[j]) < F(r, x1, x2))                                           
                y[j] = psi[j] * h;                                                                               
            else                                                                                                 
            {                                                                                                    
                y[j] = 0;                                                                                           
            }                                                                                                    
            
            if (F(r, x1 + y[j], x2 + y[j]) < F(r, x1, x2))                                                       
                min_psi[j] = y[j] * (1 / h);                                                                         
            else                                                                                                 
                min_psi[j] = 0;                                                                                  
            
            cout << F(r, x1, x2) << "  ||  " << setw(7) << F(r, x1 + y[j], x2 + y[j])
                << "  ||  " << min_psi[j] << endl;
            
        }
 
        //x1 += min_psi*h;
        //x2 += min_psi*h;
    
        //cout << "\n" << min_psi << "\n\n";    
        cout << "-----------------------------------" << endl;
/////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
    }
    _getch();
    return 0;
}

Добавлено через 16 часов 13 минут
Я тут, это доработал. Вроде все работает (мб даже правильно).
Проверьте математику пожалуйста.
Если математика верна, то тогда работает правильно.

P.S. Желательно побыстрее

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
int main()
{
    setlocale(LC_ALL, "Russian");
 
    double x1 = .3, x2 = .1, r = 1.;
 
///////////////////////////////////////////////////////////////////////////////////////////////////////////////
/*                            Модифицированный метод наилучшей пробы                                         */
                                                                                                             
    for (int l = 0; l < 5; l++)
    {
    ///Исходные данные////                                                                                     
    const int n = 10;   //                                                                                   
    double h = .1;      //                                                                                   
    //////////////////////                                                                                   
                                                                                                             
    double psi[n], y[n];                                                                                                                                                                 
                                                                                                             
 
        cout << "----------------------" << endl;
        cout << "    Итерация " << l << "\t     |" <<"\n";
        cout << "----------------------" << endl;
 
        for (int i = 0; i < n; i++)                                                                              
        {                                                                                                        
            psi[i] = ((rand() % 100 / 100. > .5) ? (1) : (-1))*(rand() % 100) / 100.;                            
            cout << "psi  "<< setw(6) << psi[i] << endl;                                                                     
        }                                                                                                        
            
        cout << "------------------" << endl;
        
        for (int j = 0; j < n; j++)                                                                              
        {                                                                                                        
            if (F(r, x1 + h*psi[j], x2 + h*psi[j]) < F(r, x1, x2))
                y[j] = psi[j] * h;
            else
                y[j] = 0;
        
            
            cout << F(r, x1, x2) << "  ||  " << setw(7) << F(r, x1 + y[j], x2 + y[j])
                << endl;
            
        }
 
        double min_xi = 0, min_fun = F(r, x1, x2);
 
        cout << "------------------";
 
        for (int k = 0; k < n; k++)
        {
            if (min_fun > F(r, x1 + y[k], x2 + y[k]))
            { 
                min_xi = psi[k];
                min_fun = F(r, x1 + y[k], x2 + y[k]);
            }
        }
 
        cout << "\n  " << min_xi << "  ||  " << setw(7) << F(r, x1 + h*min_xi, x2 + h*min_xi) << endl;
    }
 
    _getch();
    return 0;
}

@xSimplex · 10.09.2015, 22:51 **[ТС]**

Снято с вопроса

Рабочий код:

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
#include <iomanip>
#include <math.h>
#include <iostream>
#include <iomanip>
 
using namespace std;
 
double get_a(double x1, double x2)
{
    double a = 0;
    
    if ((x1 + x2) > 3) a = 1;
    if ( x2 < -1) a = 1;
    if ( x1 < -4 ) a = 1;
 
    return a;
}
 
//впомогательная функция
double F(double x1, double x2)
{
    double a = get_a(x1, x2);
 
    return 3*(x1 - 1)*(x1 - 1) + 9 * (x2 - 6)*(x2 - 6) +
        a*((x1 + x2 - 3)*(x1 + x2 - 3) + (x1 + 4)*(x1 + 4) + (x2 + 1)*(x2 + 1));;
}
 
 
// Модифицированный метод наилучшей пробы
double MMP2(double x1, double x2)
{
    const int n = 10;
    double h = .1, xi[n], chek = F(x1, x2), min_xi;
 
    for (int i = 0; i < n; i++)
    {
        xi[i] = ((rand() % 100 / 100. > .5) ? (1) : (-1))*(rand() % 100) / 100.;
 
        if (F(x1 + h*xi[i], x2 + h*xi[i]) < chek)
        {
            min_xi = xi[i]*h;
            chek = F(x1 + h*xi[i], x2 + h*xi[i]);
        }
        else
        {
            min_xi = 0;
            chek = F(x1 + h*xi[i], x2 + h*xi[i]);
        }
 
        if (F(x1 + h*xi[i], x2 + h*xi[i]) >= chek)
        {
            h *= rand() % 100 / 100. + .5;
        }
 
    }
 
    return min_xi;
}

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
#include <iostream>
#include <conio.h>
 
#include "sup.h"
 
using namespace std;
 
int main()
{
    setlocale(LC_ALL, "Russian");
 
    double x1 = -3., x2 = 4., xi;
 
    for (int j = 0; j < 148; j++)
    {
        xi = MMP2(x1, x2);
 
        x1 += xi;
        x2 += xi;
 
            cout << "-----------------------------------------" << endl;
            cout << "|| " << setw(9) << x1 << " || " << setw(9) << x2 << " || " << setw(9) << F(x1, x2) << " || " << endl;
 
    }
    
 
    _getch();
    return 0;
 
}

Добавлено через 9 минут
Только это, кому надо сами доведите до ума как прерывать, цикл при достижении минимума, а то я это еще не дописал.
Если могете скинте, а то я ща больной сижу голова совсем не варит

Мой вариант

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
    const int M = 500;
    double x1[M], x2[M], xi, x1c, x2c;
 
    x1[0] = -3.;
    x2[0] = 4.;
 
    for (int j = 0; j <= M; j++)
    {
        xi =  MMP2(x1[j], x2[j]);
 
        x1[j] += xi;
        x2[j] += xi;
 
        x1c = x1[j];
        x2c = x2[j];
 
        if (F(x1c, x2c) < F(x1[j], x2[j]))
        {
        cout << "-----------------------------------------" << endl;
        cout << "|| " << setw(9) << x1 << " || " << setw(9) << x2 << " || " << setw(9) << F(x1[j], x2[j]) << " || " << endl; 
        }
        else
        {
            cout << "Минимум найден на итерации: " << j - 1 << "\n" <<"Значение фу - ии в минимуме :" << F(x1[j], x2[j]) 
                << "\n"<< "Значение минимумов :"<< x1 << x2 << endl;
            break;
        }
    }

Новые блоги и статьи Все статьи Все блоги /
SDL3 для Web (WebAssembly): Работа со звуком через SDL3_mixer 8Observer8 08.02.2026 Содержание блога Пошагово создадим проект для загрузки звукового файла и воспроизведения звука с помощью библиотеки SDL3_mixer. Звук будет воспроизводиться по клику мышки по холсту на Desktop и по. . .	SDL3 для Web (WebAssembly): Основы отладки веб-приложений на SDL3 по USB и Wi-Fi, запущенных в браузере мобильных устройств 8Observer8 07.02.2026 Содержание блога Браузер Chrome имеет средства для отладки мобильных веб-приложений по USB. В этой пошаговой инструкции ограничимся работой с консолью. Вывод в консоль - это часть процесса. . .	SDL3 для Web (WebAssembly): Обработчик клика мыши в браузере ПК и касания экрана в браузере на мобильном устройстве 8Observer8 02.02.2026 Содержание блога Для начала пошагово создадим рабочий пример для подготовки к экспериментам в браузере ПК и в браузере мобильного устройства. Потом напишем обработчик клика мыши и обработчик. . .	Философия технологии iceja 01.02.2026 На мой взгляд у человека в технических проектах остается роль генерального директора. Все остальное нейронки делают уже лучше человека. Они не могут нести предпринимательские риски, не могут. . .
SDL3 для Web (WebAssembly): Вывод текста со шрифтом TTF с помощью SDL3_ttf 8Observer8 01.02.2026 Содержание блога В этой пошаговой инструкции создадим с нуля веб-приложение, которое выводит текст в окне браузера. Запустим на Android на локальном сервере. Загрузим Release на бесплатный. . .	SDL3 для Web (WebAssembly): Сборка C/C++ проекта из консоли 8Observer8 30.01.2026 Содержание блога Если вы откроете примеры для начинающих на официальном репозитории SDL3 в папке: examples, то вы увидите, что все примеры используют следующие четыре обязательные функции, а. . .	SDL3 для Web (WebAssembly): Установка Emscripten SDK (emsdk) и CMake для сборки C и C++ приложений в Wasm 8Observer8 30.01.2026 Содержание блога Для того чтобы скачать Emscripten SDK (emsdk) необходимо сначало скачать и уставить Git: Install for Windows. Следуйте стандартной процедуре установки Git через установщик. . . .	SDL3 для Android: Подключение Box2D v3, физика и отрисовка коллайдеров 8Observer8 29.01.2026 Содержание блога Box2D - это библиотека для 2D физики для анимаций и игр. С её помощью можно определять были ли коллизии между конкретными объектами. Версия v3 была полностью переписана на Си, в. . .