Форум программистов, компьютерный форум CyberForum.ru

Метод Хука Дживса - C++

Восстановить пароль Регистрация
 
Рейтинг: Рейтинг темы: голосов - 19, средняя оценка - 4.84
Demacr_s
0 / 0 / 0
Регистрация: 13.11.2012
Сообщений: 21
28.12.2012, 23:32     Метод Хука Дживса #1
Здравствуйте!Пишу я значит программу по нахождению минимума методом Хука Дживса.Есть книга с примером данного метода на бейсике.Но не работает она как надо.У меня уже мозг кипит.Можете помочь?
Книжка не прикрепляется(называется банди методы оптимизации)(страница 39 в книжке)
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
#include <iostream>
#include <math.h>
using namespace std;
 
double myFun(double *x);
void xy(double *x,double *y);
 
void main()
{
    setlocale(LC_ALL,"Russian");
 
    int n,ps,bs;
 
    cout << "Введите количество координат: ";
    cin >> n;
    //n++;    // One cell for value of function 
    
    double *x,*y,*b,*p;
    x=(double*)calloc(n,sizeof(int));
    y=(double*)calloc(n,sizeof(int));
    b=(double*)calloc(n,sizeof(int));
    p=(double*)calloc(n,sizeof(int));
 
    int i,j;
    cout << "Введите начальные точки: ";
    for(i=0; i<n; i++)
        cin >> x[i];
 
    double h,k,fi,fz,fb;
    cout << "Введите длинну шага: ";
    cin >> h;
    k=h;
 
    for(i=0; i<n; i++){
        y[i]=x[i];
        p[i]=x[i];
        b[i]=x[i];
    }
 
    fi=myFun(x);
 
    cout << "Начальное значение функции " << fi << endl;
    for(i=0;i<n;i++)
        cout << x[i] << endl;
 
    ps=0;
    bs=0;
    j=0;
    fb=fi;
    label_name:
    //Точка возврата
    x[j]=y[j]+k;
 
    fz=myFun(x);
 
    if (fz>=fi){
        x[j]=y[j]-k;
        fz=myFun(x);
        if(fz>=fi)
            x[j]=y[j];
    }
    else
        y[j]=x[j];
    fi=myFun(x);
    
    cout << "Исследующий поиск" << endl;
    for(i=0; i<n; i++)
        cout << x[i] << endl;
 
    if (j!=n)
        j=j++;//Возврат
    cout << j << "  " << n << endl;
    goto label_name;
        if (fi>=fb-(1e-08)){
        if (ps==0&&bs==1)
                for(i=0; i<n; i++){
                    p[i]=b[i];
                    y[i]=b[i];
                    x[i]=b[i];
                }
    bs=1;
    ps=0;
    fi=myFun(x);
    fb=myFun(x);
    cout << "Замена базисной точки" << endl;
    for(i=0; i<n; i++)
        cout << x[i] << endl;
    j++;
    goto label_name;//Возврат
    }
    
    k=k/10;
    cout << "Уменьшить длину шага" << endl;
    if (k>=(1e-08)){
        j=1;
        goto label_name;//Возврат
        for(i=0;i<n;i++){
            p[i]=2*y[i]-b[i];
            b[i]=y[i];
            x[i]=p[i];
            y[i]=x[i];
        }
        fb=fi;
        ps=1;
        bs=0;
        fi=myFun(x);
        cout << "Поиск по образцу" << fi << endl;
        for(i=0;i<n;i++)
            cout << x[i] << endl;
        j=1;
        goto label_name;//Возврат
        cout << "Минимум функции найден" << endl;
        for (i=0;i<n;i++)
            cout << "X=" << p[i] << endl;
        cout << "Минимум функции = " << fb << endl;
    }
    system("pause");
    free(x);
    free(b);
    free(p);
}
 
 
double myFun(double *x)
{
    double result=pow((x[0]-2),2.)+pow((x[1]-5),2.)+pow((x[2]+2),4.) ;
    return result;
}
Миниатюры
Метод Хука Дживса  
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
28.12.2012, 23:32     Метод Хука Дживса
Посмотрите здесь:

метод деления отрезка пополам и метод итерации C++
Метод Эйлера, и Метод Лагранжа, в долгу не останусь C++
Метод деления отрезка пополам для решения нелинейных уравнений (метод дихотомии) C++
C++ Метод медиан из трех элементов VS улучшенный быстрый метод сортировки(метод Бентли-Макилроя)
C++ Производный класс: метод возведения в произвольную степень, и метод для вычисления логарифма числа
Класс vector (поля: координаты, 2 конструктора, метод нахождения длины вектора и метод вывода координат на экран) C++
C++ СЛАУ. Метод обратной матрицы, метод Гаусса, метод Крамера, метод Зейделя
Нахождения корней уравнения: метод половинного деления (бисекции) или метод хорд C++

Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Ответ Создать тему
Опции темы

Текущее время: 13:19. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru