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

Методом Нелдера-Мида. Ошибка в коде. - C++

Войти
Регистрация
Восстановить пароль
Другие темы раздела
C++ функция вычисления моды ряда http://www.cyberforum.ru/cpp-beginners/thread411017.html
Здравствуйте. есть задание - написать функцию нахождение моды ряда((наиболее часто встречаемое число) параметры функции - ряд и количество элементов в ряде) накидал код, но как дальше его реализовать именно в функции не понимаю. #include <iostream.h> #include <conio.h> int moda(int n,int **mas); //prototip void main(){ } int moda(int n,int **mas){ int x,i,j,b;
C++ очериди C++ Смоделировать очередь на базе статического массива, в которой реализован такой механизм изъятия элемента из очереди: начало очереди всегда находится в первом элементе массива; при изъятии одного элемента из очереди, все остальные элементы перемещаются на одну позицию ближе к началу массива. После обработки всей заданной входной последовательности найти среднее арифметическое всех элементов... http://www.cyberforum.ru/cpp-beginners/thread411012.html
поиск связных компонент графа C++
Нужна програмная реализация алгоритма поиска связных компонент графа написала вот что-то такое: #pragma hdrstop #pragma argsused #include <iostream> #include <conio.h> #include <time.h>
Перенаправление потока вывода C++
Здравствуйте! Задача такова - нужно перенаправить поток вывода с cout на файловый в одном из конструкторов класса. Итог - необработанное исключение 0xC0000005: Нарушение прав доступа при чтении "0xccccccd0". Вот сами исходники: 1. logger.h #ifndef LOGGER_H #define LOGGER_H #include <iostream> #include <fstream> class Logger {
C++ Динамический массив http://www.cyberforum.ru/cpp-beginners/thread410993.html
Разработать класс содержащий: - компоненты данные - методы: а) конструктор по умолчанию; б) конструктор с параметрами; в) конструктор копирования; г) деструктор; д) методы класса для работы с данными. Данные класса должны иметь атрибут private. Задача 1. Создать 2 объекта разработанного класса. Класс – динамический вектор (одномерный массив). Найти в каждом объекте MAX значение элемента...
C++ Модуль в массиве Добрый вечер. Подскажите мне пожалуйста вот с этой задачей: В одномерном массиве, состоящем из n вещественных элементов, вычислить: • номер минимального элемента; • сумму элементов массива, расположенных между первым и вторым отрицательными элементами. • Преобразовать таким образом, чтобы сначала располагались все элементы, модуль которых не превышает 1, а потом – все остальные. ... подробнее

Показать сообщение отдельно
YouDoItWrong
46 / 46 / 7
Регистрация: 29.10.2011
Сообщений: 154

Методом Нелдера-Мида. Ошибка в коде. - C++

18.12.2011, 21:50. Просмотров 1630. Ответов 1
Метки (Все метки)

Ест схема (см. вложения) метода Нелдера-Мида. По ней я написал код, но в нем есть ошибка логическая и/или математическая. Пожалуйста, помогите найти, я что-то не могу ее найти . Должно получиться (x,y,z) ~ (2,-6,-7). Заранее спасибо.
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
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
#include <vector>
#include <numeric>
#include <algorithm>
#include <iostream>
#include <cmath>
#define N 3
#define T 1.0
 
using namespace std;
 
struct point {
    double x,y,z;
    point( double cX, double cY, double cZ ) {
        x = cX; y = cY; z = cZ;
    }
    point point::operator +( point &other ) {
        return point( x + other.x, y + other.y, z + other.z );
    }
    point point::operator *( double &other ) {
        return point( x*other, y*other, z*other);
    }
    point point::operator -( point &other ) {
        return point( x - other.x, y - other.y, z - other.z );
    }
    void point::operator =( point &other ) {
        x = other.x;
        y = other.y;
        z = other.z;
    }
    void point::operator +=( point &other ) {
        x += other.x;
        y += other.y;
        z += other.z;
    }
    void Show() {
        cout << x << " " << y << " " << z << endl;
    }
    double Function() {
        return 3 * pow( x - 2, 2 ) + 3 * pow( y + 6, 2 ) + 4 * pow( z + 7, 2 );
    }
 
};
 
double Test(vector<double> vec, double C) {
    double tmp = double();
    for( int i = 0; i < vec.size(); i ++) {
        tmp += pow(vec[i] - C, 2.0);
    }
    return sqrt( ( 1.0 / ((double)N + 1.0) ) * tmp );
}
 
int NotHorL(int I) {
    for( int i = 0; i <= N; i++) {
        if( i != I ) {
            return i;
        }
    }
    return -1;
}
 
int main() {
    vector<point> vectorOfPoints, vectorOfC;
    vector<double> vectorOfFunctions, tmp;
    point R( 0, 0, 0 ), C( 0, 0, 0 ), S( 0, 0, 0 ), E( 0, 0, 0 ), O( 0, 0, 0 );
    double dFh, dFl, temp, eps, alpha = 1, betha = 0.5, gama = 3.0;
    int iL, iH;
 
    
    for(int j = 0; j <=N; j++) {
        for(int k = 0; k < N; k++) {
            if(j == 0) {
                tmp.push_back(0);
                continue;
            }
            if(j == k+1) {
                tmp.push_back( (T / (N * sqrt(2.0))) * ( sqrt( (double)N + 1.0 - 1 + (double)N ) ) );
                continue;
            }
            tmp.push_back( (T / (N * sqrt(2.0))) * ( sqrt( (double)N + 1.0 - 1 ) ) );
        }
    }
 
    for(int i = 0; i < tmp.size(); i++) {
        vectorOfPoints.push_back(point( tmp[i], tmp[i + 1], tmp[i + 2] ) );
        i += 2;
    }
    tmp.clear();
    //
    cout << "eps: ";
    cin >> eps;
    while(true) {
        for(int i = 0; i < vectorOfPoints.size(); i++) {
            vectorOfFunctions.push_back(vectorOfPoints[i].Function());
        }
        
        dFh = *max_element(vectorOfFunctions.begin(), vectorOfFunctions.end());
        dFl = *min_element(vectorOfFunctions.begin(), vectorOfFunctions.end());
        iL = (min_element(vectorOfFunctions.begin(), vectorOfFunctions.end()) - vectorOfFunctions.begin());
        iH = (max_element(vectorOfFunctions.begin(), vectorOfFunctions.end()) - vectorOfFunctions.begin());
 
        for(int i = 0; i < vectorOfPoints.size(); i++) {
            O += vectorOfPoints[i];
        }
 
        temp =  (double)( 1.0 / (double)N );
        C = ( O * temp) - vectorOfPoints[iH];
 
        cout << Test(vectorOfFunctions, C.Function()) << endl;
        if(Test(vectorOfFunctions, C.Function()) < eps) {
            break;
        }
 
        O = C - vectorOfPoints[iH];
        R = C + ( O * alpha );
 
        if(R.Function() < dFl) { //Ветвь 1
            O = R - C;
            E = C + ( O * gama );
 
            if( E.Function() < R.Function() ) {
                vectorOfPoints[iH] = E;
            } else {
                vectorOfPoints[iH] = R;
            }
        } else { // Ветвь 2
            if( R.Function() > vectorOfFunctions[NotHorL(iH)] ) {
                if( R.Function() < vectorOfPoints[iH].Function() ) {
                    vectorOfPoints[iH] = R;
                }
 
                O = vectorOfPoints[iH] - C;
                S = C + ( O * betha );
 
                if( S.Function() < vectorOfFunctions[iH] ) {
                    vectorOfPoints[iH] = S;
                } else {
                    O = vectorOfPoints[NotHorL(iL)] - vectorOfPoints[iL]; temp = 0.5;
                    vectorOfPoints[NotHorL(iL)] = vectorOfPoints[iL] + ( O * temp );
                }
            } else {
                vectorOfPoints[iH] = R;
            }
        }
        vectorOfFunctions.clear();
        O = point( 0, 0, 0 );
        C.Show();
    }
    cout << "Final result: ";
    C.Show();
    return 0;
}
Миниатюры
Методом Нелдера-Мида. Ошибка в коде.  
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
 
Текущее время: 01:13. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru