Реализация метода Хука-Дживса. Поиск ошибки в коде

@Doc_ff · Регистрация: 26.10.2015

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

Попытался реализовать простой способ нахождения минимума с помощью метода Хука - Дживса. Но не как не могу понять почему он работает не правильно. Просьба как говорится, грязью не кидаться, пошуршал по поиску, но слишком всё как то закрученно. Хочется, нечто проще реализовать. Вот пример моего кода, может кто сможет мне помочь найти ошибку.

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
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
# include <iostream>
# include <cmath>
# include <vector>
 
using namespace std;
using standard_type = double ;
 
vector <standard_type> F (const vector <standard_type>& v)
{
    return { v[0] * v[0] + v[1] * v[1] } ;
}
 
vector <standard_type> get_initial_condition ()
{
    return { -5, -5 } ;
}
 
// vector <standard_type> F (const vector <standard_type>& v)
// {
//     return { 8 * pow ( v[0], 2 ) + 4 * ( v[0] * v[1] ) + 5 * pow ( v[1], 2 ) } ;
// }
 
// vector <standard_type> get_initial_condition ()
// {
//     return { -4, -4 } ;
// }
 
void print_debug_info (const vector <standard_type>& v, string name_vector)
{
    standard_type size_vector = v.size() ;
    cout << " \t print debug info: " << endl;
    cout << name_vector << " = { ";
    for (int i = 0; i < size_vector - 1; i++)
    {
        cout << v[i] << " ;" ;
    }
    cout << v[size_vector - 1] << " } ;" << endl;
    cout << " \t print debug info ends. " << endl;  
}
 
void print_vector (const vector<standard_type>& v, const string name_vector)
{
    cout << name_vector << " = { ";
    for (int i = 0; i < v.size() - 1; i++)
    {
        cout << v[i] << " ;" ;
    }
    cout << v.back() << " } ;" << endl;
}
 
vector<standard_type> search_by_pattern(vector<standard_type> previous_point, vector<standard_type> found_point);
 
vector <standard_type> HJ_method (standard_type step = 1, standard_type eps = 0.1)
{
    standard_type number_of_iterations = 0 ;
    standard_type number_of_variables = get_initial_condition().size();
    vector <standard_type> x_basis = get_initial_condition() ;
    vector <standard_type> F_result_basis = F(x_basis) ;
    vector <standard_type> new_point (number_of_variables);
    vector <standard_type> F_result_new_point(F_result_basis.size()) ;
 
    do
    {
        number_of_iterations++;
        F_result_basis = F(x_basis);
 
        for (size_t i = 0; i < number_of_variables; i++)
        {
            new_point.clear();
            F_result_new_point.clear();
            new_point = x_basis ;
            new_point[i] = x_basis[i] + step;
            print_debug_info(x_basis, "x_basis");
            print_debug_info(new_point, "new_point");
            F_result_new_point = F(new_point);
            print_debug_info(F_result_new_point, "F_result_new_point");
            if (F_result_basis > F_result_new_point)
            {
                cout << "F_result_basis > F_result_new_point => " << F_result_basis[0] << " > " << F_result_new_point[0] << endl;
                x_basis = search_by_pattern(x_basis, new_point);
                print_debug_info(x_basis, " NEW x_basis ");
                cout << "number_of_iterations = " << number_of_iterations << endl;
                break;
            }
 
            new_point.clear();
            F_result_new_point.clear();
            new_point = x_basis ;
            new_point[i] = x_basis[i] - step;
            print_debug_info(x_basis, "x_basis");
            print_debug_info(new_point, "new_point");
            F_result_new_point = F(new_point);
            print_debug_info(F_result_new_point, "F_result_new_point");
            if (F_result_basis > F_result_new_point)
            {
                cout << "F_result_basis > F_result_new_point => " << F_result_basis[0] << " > " << F_result_new_point[0] << endl;
                x_basis = search_by_pattern(x_basis, new_point);
                print_debug_info(x_basis, " NEW x_basis ");
                cout << "number_of_iterations = " << number_of_iterations << endl;
                break;
            }
        }
 
        step = step / 2 ;
        cout << "step_reduction. NEW step = " << step << endl;
 
 
        
    } while (step > eps);
 
    cout << "step > eps => " << step << " > " << eps << endl;
    cout << "number_of_iterations = " << number_of_iterations << endl;
 
    return x_basis ;
}
 
vector<standard_type> search_by_pattern(vector<standard_type> previous_point, vector<standard_type> found_point)
{
    vector<standard_type> v(previous_point.size());
 
    for (int i = 0; i < v.size(); i++)
    {
        v[i] = previous_point[i] + 2 * (found_point[i] - previous_point[i]);
    }
 
    cout << " --------- " << endl;
    print_debug_info(v, "search_by_pattern () -> found_point");
    cout << " --------- " << endl;
 
    return v ;
}

@БедолагаЖека · 22.10.2019, 10:06

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

using standard_type = double ;

а так можно?
typedef вроде надо...

Добавлено через 2 минуты

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

vector <standard_type> F (const vector <standard_type>& v)
vector <standard_type> get_initial_condition ()

Эти вещи тоже вроде токо через обертку надо делать

Алексей1153 · 22.10.2019, 10:12

Не по теме:

Сообщение от БедолагаЖека

а так можно

можно по новым стандартам

Сообщение от БедолагаЖека

Эти вещи тоже вроде токо через обертку надо делать

вектор на стеке мало места занимает, так что можно. А оптимизатор вообще может сам ссылку на внешний объект потом передать в релизе

@БедолагаЖека · 22.10.2019, 10:41

Сообщение от Алексей1153

вектор на стеке мало места занимает, так что можно. А оптимизатор вообще может сам ссылку на внешний объект потом передать в релизе

Да это я вообще гоню... показалось переопределение методов...

Новые блоги и статьи Все статьи Все блоги /
Программный контроль заполнения реквизита табличной части документа Maks 02.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: реализовать контроль заполнения реквизита "ПричинаСписания". . .	wmic не является внутренней или внешней командой Maks 02.04.2026 Решение: DISM / Online / Add-Capability / CapabilityName:WMIC~~~~ Отсюда: https:/ / winitpro. ru/ index. php/ 2025/ 02/ 14/ komanda-wmic-ne-naydena/	Программная установка даты и запрет ее изменения Maks 02.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: при создании документов установить период списания автоматически. . .	Вывод данных в справочнике через динамический список Maks 01.04.2026 Реализация из решения ниже выполнена на примере нетипового справочника "Спецтехника" разработанного в конфигурации КА2. Задача: вывести данные из ТЧ нетипового документа. . .
Программное заполнения текстового поля в реквизите формы документа Maks 01.04.2026 Алгоритм из решения ниже реализован на нетиповом документе "ВыдачаОборудованияНаСпецтехнику" разработанного в конфигурации КА2, в дополнении к предыдущему решению. На форме документа создается. . .	К слову об оптимизации kumehtar 01.04.2026 Вспоминаю начало 2000-х, университет, когда я писал на Delphi. Тогда среди программистов на форумах активно обсуждали аккуратную работу с памятью: нужно было следить за переменными, вовремя. . .	Идея фильтра интернета (сервер = слой+фильтр). Hrethgir 31.03.2026 Суть идеи заключается в том, чтобы запустить свой сервер, о чём я если честно мечтал давно и давно приобрёл книгу как это сделать. Но не было причин его запускать. Очумелые учёные напечатали на. . .	Модель здравосоХранения 6. ESG-повестка и устойчивое развитие; углублённый анализ кадрового бренда anaschu 31.03.2026 В прикрепленном документе раздумья о том, как можно поменять модель в будущем