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

@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;
        }
    }

Новые блоги и статьи Все статьи Все блоги /
PhpStorm 2025.3: WSL Terminal всегда стартует в ~ and_y87 14.12.2025 PhpStorm 2025. 3: WSL Terminal всегда стартует в ~ (home), игнорируя директорию проекта Симптом: После обновления до PhpStorm 2025. 3 встроенный терминал WSL открывается в домашней директории. . .	Как объединить две одинаковые БД Access с разными данными VikBal 11.12.2025 Помогите пожалуйста !! Как объединить 2 одинаковые БД Access с разными данными.	Новый ноутбук volvo 07.12.2025 Всем привет. По скидке в "черную пятницу" взял себе новый ноутбук Lenovo ThinkBook 16 G7 на Амазоне: Ryzen 5 7533HS 64 Gb DDR5 1Tb NVMe 16" Full HD Display Win11 Pro	Музыка, написанная Искусственным Интеллектом volvo 04.12.2025 Всем привет. Некоторое время назад меня заинтересовало, что уже умеет ИИ в плане написания музыки для песен, и, собственно, исполнения этих самых песен. Стихов у нас много, уже вышли 4 книги, еще 3. . .	От async/await к виртуальным потокам в Python IndentationError 23.11.2025 Армин Ронахер поставил под сомнение async/ await. Создатель Flask заявляет: цветные функции - провал, виртуальные потоки - решение. Не threading-динозавры, а новое поколение лёгких потоков. Откат?. . .
Поиск "дружественных имён" СОМ портов Argus19 22.11.2025 Поиск "дружественных имён" СОМ портов На странице: https:/ / norseev. ru/ 2018/ 01/ 04/ comportlist_windows/ нашёл схожую тему. Там приведён код на С++, который показывает только имена СОМ портов, типа,. . .	Сколько Государство потратило денег на меня, обеспечивая инсулином. Programma_Boinc 20.11.2025 Сколько Государство потратило денег на меня, обеспечивая инсулином. Вот решила сделать интересный приблизительный подсчет, сколько государство потратило на меня денег на покупку инсулинов. . . .	Ломающие изменения в C#.NStar Alpha Etyuhibosecyu 20.11.2025 Уже можно не только тестировать, но и пользоваться C#. NStar - писать оконные приложения, содержащие надписи, кнопки, текстовые поля и даже изображения, например, моя игра "Три в ряд" написана на этом. . .	Мысли в слух kumehtar 18.11.2025 Кстати, совсем недавно имел разговор на тему медитаций с людьми. И обнаружил, что они вообще не понимают что такое медитация и зачем она нужна. Самые базовые вещи. Для них это - когда просто люди. . .	Создание Single Page Application на фреймах krapotkin 16.11.2025 Статья исключительно для начинающих. Подходы оригинальностью не блещут. В век Веб все очень привыкли к дизайну Single-Page-Application . Быстренько разберем подход "на фреймах". Мы делаем одну. . .