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

Необходимо изменить программу в C++ - C++

Восстановить пароль Регистрация
 
Рейтинг: Рейтинг темы: голосов - 18, средняя оценка - 4.94
sl1m_dogg
0 / 0 / 0
Регистрация: 20.01.2011
Сообщений: 18
20.01.2011, 19:09     Необходимо изменить программу в C++ #1
Здравствуйте. Собственно есть программа, которая минимизирует модель методом сопряженных направлений Пауэлла.
Код:
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
#include <iostream.h> //Метод сопряженных направлений
#include <math.h>
#include <conio.h>
 
double f(double  x1,double x2) {
    return 4*(x1-5)*(x1-5)+(x2-6)*(x2-6);}
int main(){
    double E=0.1, min, var[200], t;
    int i, j, n=2, k=0;
    int s=0, I[100];
 
    double d[3][3]={{0,1}, {1,0}, {0,1}};    cout <<d[0][0]<<" "<<d[0][1]<<" "<<
                            d[1][0]<<" "<<d[1][1]<<" "<<
                            d[2][0]<<" "<<d[2][1]<<endl;
    double X[20][20], Y[20][40], Xz[1][2];
    I[s]=0;
    cout <<"X0:"<<endl;           //задаем начальную точку Хо
    for(i=0;i<1;i++)
      for(int j=0;j<2;j++)
    cin >>X[i][j];           cout <<"X("<<X[0][0]<<";"<<X[0][1]<<")\n";
      for(i=0;i<1;i++)             //Y0 = Xo
      for(j=0;j<2;j++)
    Y[i][j]=X[i][j];         cout <<"Y("<<Y[0][0]<<";"<<Y[0][1]<<")\n";
//.2.........................
TWO:  cout <<" __2__"<<endl;
      min=f(Y[I[s]][0],Y[I[s]][1]);
      i=0;
      for(t=-50;t<50;t++,i++)    //вывод различных значений ф-ции
    {
      var[i]=f( Y[I[s]][0]+t*d[I[s]][0], Y[I[s]][1]+t*d[I[s]][1]);
      cout <<"var ["<<i<<"] = "<<var[i]<<endl;
    }
      for(i=0;i<100;i++)      //вывод наименьшего значения ф-ции
      if(var[i] < min)
        min = var[i];
        cout <<"min var[] = "<<min<<endl;
      for(int index=0;index<100;index++)  //вывод t
        { 
          if( min==var[index] )
          t=index-50;
        }
      cout <<"T = "<<t<<endl;
      Y[I[s]+1][0] = Y[I[s]][0] + t*d[I[s]][0];
      Y[I[s]+1][1] = Y[I[s]][1] + t*d[I[s]][1];  
                                           cout <<"Y["<<I[s]+1<<"] = ("<<Y[I[s]+1][0]<<";";
                                           cout <<Y[I[s]+1][1]<<")"<<endl;
//.3.................
THREE:  cout <<" __3__"<<endl;  
        if(I[s] < n-1)
        {
            cout <<"I[s] < n-1"<<endl;    
            I[s]=I[s]+1;
            cout <<"I[s] = "<<I[s]<<endl;
            goto TWO;    
        }   //*/
        if(I[s] == n-1)
        {   
            cout <<"I[s] == n-1"<<endl;
            if( (Y[n][0]==Y[0][0]) && (Y[n][1]==Y[0][1]) ){
               Xz[0][0]=Y[n][0];        cout <<"Xz("<<Xz[0][0]<<";";
               Xz[0][1]=Y[n][1];        cout <<Xz[0][1]<<")"<<endl;
            }
            else {
            cout <<"Yn != Y0 "<<endl;
            I[s]=I[s]+1;
             cout <<"I[s] = "<<I[s]<<endl;
            goto TWO; }  
        }
        if(I[s] == n)
        {
            cout <<"I[s]  == n"<<endl;    
            if( (Y[n+1][0]==Y[1][0]) && (Y[n+1][1]==Y[1][1]) ){
               Xz[0][0]=Y[n+1][0];        cout <<"Xz("<<Xz[0][0]<<";";
               Xz[0][1]=Y[n+1][1];        cout <<Xz[0][1]<<")"<<endl;
            }
            if( (Y[n+1][0]!=Y[1][0]) || (Y[n+1][1]!=Y[1][1]) ){
            goto Fourth; 
            }
        }
//.4.........................        
Fourth:  cout <<" __4__"<<endl;  
         X[k+1][0]=Y[n+1][0];    cout <<"X["<<(k+1)<<"] = ("<<X[k+1][0];
         X[k+1][1]=Y[n+1][1];    cout <<";"<<X[k+1][1]<<")"<<endl; 
//...........a...............
         if( pow( ( pow((X[k+1][0]-X[k][0]),2)+pow((X[k+1][1]-X[k][1]),2) ),1/2.0)<E ){
         cout <<"X(k+1) - X(k) = "<<pow( ( pow((X[k+1][0]-X[k+1][0]),2)+
                                    pow((X[k+1][1]-X[k+1][1]),2) ),1/2.0);
             Xz[0][0]=X[k+1][0];
             Xz[0][1]=X[k+1][1];    
         }
//..................b.....................
         if( pow( ( pow((X[k+1][0]-X[k][0]),2)+pow((X[k+1][1]-X[k][1]),2) ),1/2.0)>E ){
         d[0][0]=d[n][0]=Y[n+1][0]-Y[1][0];    cout <<"d0 = ("<<d[0][0];
             d[0][1]=d[n][1]=Y[n+1][1]-Y[1][1];    cout <<";"<<d[0][1]<<")"<<endl;
             
             d[I[s]][0]=d[I[s]+1][0];
             d[I[s]][1]=d[I[s]+1][1];
             if( (d[1][0]*d[2][1]-d[1][1]*d[2][0]) != 0){
                cout <<"rang(D1,...,Dn) = n"<<endl;
                Y[0][0]=X[k+1][0];  cout <<"Y0 = ("<<Y[0][0]<<";";
                Y[0][1]=X[k+1][1];  cout <<Y[0][1]<<")\n";
                k=k+1;
                I[s]=0;
                goto TWO;                                    }  
             if( (d[1][0]*d[2][1]-d[1][1]*d[2][0]) == n){
                cout <<"rang(D1,...,Dn) < n"<<endl;
                Y[0][0]=X[k+1][0];
                Y[0][1]=X[k+1][1];
                k=k+1; I[s]=0;
goto TWO;                                      
             } 
         }   
 cout <<"\nXz("<<Xz[0][0]<<";"<<Xz[0][1]<<")"<<endl;    
     getch();}
Она минимизирует функцию 4((x1-5)^2)+(x2-6)^2. В моем случае нужно минимизировать функцию (9*х1*x1)+(35*х2*x2)+(32*х1*х2)–(88*х1)–(176*х2)+242, пытался заменить функцию в строке
C++
1
return 4*(x1-5)*(x1-5)+(x2-6)*(x2-6);}
но выдает или неверный результат, или кучу ошибок.
Помогите кто сможет.
P.S.: результат должен быть x1=2.207 x2=1.507
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
killboss
 Аватар для killboss
57 / 57 / 2
Регистрация: 18.12.2010
Сообщений: 127
20.01.2011, 20:20     Необходимо изменить программу в C++ #2
Цитата Сообщение от sl1m_dogg Посмотреть сообщение
if( pow( ( pow((X[k+1][0]-X[k][0]),2)+pow((X[k+1][1]-X[k][1]),2) ),1/2.0)<E )
Разделите это выражение на составляющие. Должно быть:
C++
1
pow(x1, x2)
т.е. внутри pow() никаких выражений не должно быть.
sl1m_dogg
0 / 0 / 0
Регистрация: 20.01.2011
Сообщений: 18
20.01.2011, 20:33  [ТС]     Необходимо изменить программу в C++ #3
А не могли бы вы продемонстрировать как должен выглядеть этот if? А то я что то не особо понял
killboss
 Аватар для killboss
57 / 57 / 2
Регистрация: 18.12.2010
Сообщений: 127
20.01.2011, 20:42     Необходимо изменить программу в C++ #4
Жди минутку...
sl1m_dogg
0 / 0 / 0
Регистрация: 20.01.2011
Сообщений: 18
20.01.2011, 20:43  [ТС]     Необходимо изменить программу в C++ #5
Может я не совсем правильно описал. Программа правильно работает для первой функции (которая в коде прописана), она выводит правильное значение т. (5;6). Но когда я в return прописываю значение второй функции программа выводит неверное значение. Первоначальную точку вводил (0;0).
killboss
 Аватар для killboss
57 / 57 / 2
Регистрация: 18.12.2010
Сообщений: 127
20.01.2011, 20:46     Необходимо изменить программу в C++ #6
Цитата Сообщение от killboss Посмотреть сообщение
if( pow( ( pow((X[k+1][0]-X[k][0]),2)+pow((X[k+1][1]-X[k][1]),2) ),1/2.0)<E )
Вот так:
C++
1
2
3
4
float n1 = X[k+1][0]-X[k][0];
float n2 = X[k+1][1]-X[k][1];
float n3 = pow(n1,2)+pow(n2,2);
if(pow(n3,1/2) < E)
И так сделать по всему алгоритму!!!
sl1m_dogg
0 / 0 / 0
Регистрация: 20.01.2011
Сообщений: 18
20.01.2011, 20:54  [ТС]     Необходимо изменить программу в C++ #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
#include <iostream.h> //Метод сопряженных направлений
#include <math.h>
#include <conio.h>
 
double f(double  x1,double x2) {
    return (9*х1*x1)+(35*х2*x2)+(32*х1*х2)(88*х1)(176*х2)+242;}
int main(){
    double E=0.0001, min, var[200], t;
    int i, j, n=2, k=0;
    int s=0, I[100];
 
    double d[3][3]={{0,1}, {1,0}, {0,1}};
    cout <<d[0][0]<<" "<<d[0][1]<<" "<<d[1][0]<<" "<<d[1][1]<<" "<<d[2][0]<<" "<<d[2][1]<<endl;
    double X[100][100], Y[100][200], Xz[1][2];
 
    I[s]=0;
    cout <<"X0:"<<endl;           //задаем начальную точку Хо
    for(i=0;i<1;i++)
      for(j=0;j<2;j++)
    cin >>X[i][j];           cout <<"X("<<X[0][0]<<";"<<X[0][1]<<")\n";
 
    for(i=0;i<1;i++)             //Y0 = Xo
      for(j=0;j<2;j++)
    Y[i][j]=X[i][j];         cout <<"Y("<<Y[0][0]<<";"<<Y[0][1]<<")\n";
//.2.........................
TWO:  cout <<" __2__"<<endl;
      min=f(Y[I[s]][0],Y[I[s]][1]);
 
      i=0;
      for(t=-50;t<50;t++,i++)    //вывод различных значений ф-ции
    {
      var[i]=f( Y[I[s]][0]+t*d[I[s]][0], Y[I[s]][1]+t*d[I[s]][1]);
      cout <<"var ["<<i<<"] = "<<var[i]<<endl;
    }
      for(i=0;i<100;i++)      //вывод наименьшего значения ф-ции
      if(var[i] < min)
        min = var[i];
        cout <<"min var[] = "<<min<<endl;
      for(int index=0;index<100;index++)  //вывод t
    {
      if( min==var[index] )
      t=index-50;
    }
 
      cout <<"T = "<<t<<endl;
 
      Y[I[s]+1][0] = Y[I[s]][0] + t*d[I[s]][0];
      Y[I[s]+1][1] = Y[I[s]][1] + t*d[I[s]][1];
                       cout <<"Y["<<I[s]+1<<"] = ("<<Y[I[s]+1][0]<<";";
                       cout <<Y[I[s]+1][1]<<")"<<endl;
//.3............................
THREE:  cout <<" __3__"<<endl;
    if(I[s] < n-1)
    {
        cout <<"I[s] < n-1"<<endl;
        I[s]=I[s]+1;
        cout <<"I[s] = "<<I[s]<<endl;
        goto TWO;
    }   //*/
    if(I[s] == n-1)
    {
        cout <<"I[s] == n-1"<<endl;
        if( (Y[n][0]==Y[0][0]) && (Y[n][1]==Y[0][1]) ){
           Xz[0][0]=Y[n][0];        cout <<"Xz("<<Xz[0][0]<<";";
           Xz[0][1]=Y[n][1];        cout <<Xz[0][1]<<")"<<endl;
        }
        else {
        cout <<"Yn != Y0 "<<endl;
        I[s]=I[s]+1;
         cout <<"I[s] = "<<I[s]<<endl;
        goto TWO; }
    }
    if(I[s] == n)
    {
        cout <<"I[s]  == n"<<endl;
        if( (Y[n+1][0]==Y[1][0]) && (Y[n+1][1]==Y[1][1]) ){
           Xz[0][0]=Y[n+1][0];        cout <<"Xz("<<Xz[0][0]<<";";
           Xz[0][1]=Y[n+1][1];        cout <<Xz[0][1]<<")"<<endl;
        }
        if( (Y[n+1][0]!=Y[1][0]) || (Y[n+1][1]!=Y[1][1]) ){
        goto Fourth;
        }
    }
//.4.........................
Fourth:  cout <<" __4__"<<endl;
     X[k+1][0]=Y[n+1][0];    cout <<"X["<<(k+1)<<"] = ("<<X[k+1][0];
     X[k+1][1]=Y[n+1][1];    cout <<";"<<X[k+1][1]<<")"<<endl;
//...........a...............
double n1 = X[k+1][0]-X[k][0];
double n2 = X[k+1][1]-X[k][1];
double n3 = pow(n1,2)+pow(n2,2);
if(pow(n3,1/2) < E)
{
     cout <<"X(k+1) - X(k) = "<<pow( ( pow((X[k+1][0]-X[k+1][0]),2)+
                    pow((X[k+1][1]-X[k+1][1]),2) ),1/2.0);
 
         Xz[0][0]=X[k+1][0];
         Xz[0][1]=X[k+1][1];
     }
//..................b.....................
 n1 = X[k+1][0]-X[k][0];
 n2 = X[k+1][1]-X[k][1];
 n3 = pow(n1,2)+pow(n2,2);
if(pow(n3,1/2) < E)
{
 
         d[0][0]=d[n][0]=Y[n+1][0]-Y[1][0];    cout <<"d0 = ("<<d[0][0];
         d[0][1]=d[n][1]=Y[n+1][1]-Y[1][1];    cout <<";"<<d[0][1]<<")"<<endl;
 
         d[I[s]][0]=d[I[s]+1][0];
         d[I[s]][1]=d[I[s]+1][1];
         if( (d[1][0]*d[2][1]-d[1][1]*d[2][0]) != 0){
        cout <<"rang(D1,...,Dn) = n"<<endl;
        Y[0][0]=X[k+1][0];  cout <<"Y0 = ("<<Y[0][0]<<";";
        Y[0][1]=X[k+1][1];  cout <<Y[0][1]<<")\n";
        k=k+1;
        I[s]=0;
        goto TWO;
         }
         if( (d[1][0]*d[2][1]-d[1][1]*d[2][0]) == n){
        cout <<"rang(D1,...,Dn) < n"<<endl;
        Y[0][0]=X[k+1][0];
        Y[0][1]=X[k+1][1];
        k=k+1;
        I[s]=0;
        goto TWO;
         }
     }
     cout <<"\nXz("<<Xz[0][0]<<";"<<Xz[0][1]<<")"<<endl;
     getch();}
как бы вот так переписал. не знаю правильно ли. но все тоже остается, 9 ошибок для второй функции.
killboss
 Аватар для killboss
57 / 57 / 2
Регистрация: 18.12.2010
Сообщений: 127
20.01.2011, 20:56     Необходимо изменить программу в C++ #8
На что ругается компилятор? Скинь сообщение...
sl1m_dogg
0 / 0 / 0
Регистрация: 20.01.2011
Сообщений: 18
20.01.2011, 21:00  [ТС]     Необходимо изменить программу в C++ #9
в общем суть видимо не в функции pow. для первой модели оно все отлично считало и с этим кодом:
C++
1
2
3
 if( pow( ( pow((X[k+1][0]-X[k][0]),2)+pow((X[k+1][1]-X[k][1]),2) ),1/2.0)<E ){
     cout <<"X(k+1) - X(k) = "<<pow( ( pow((X[k+1][0]-X[k+1][0]),2)+
                    pow((X[k+1][1]-X[k+1][1]),2) ),1/2.0);
.
Но когда в return прописываю вместо:
C++
1
 return 4*(x1-5)*(x1-5)+(x2-6)*(x2-6);}
вот это код:
C++
1
 return (9*х1*x1)+(35*х2*x2)+(32*х1*х2)(88*х1)(176*х2)+242;}
то программа выдает неверную точку минимума (правильная точка уже известна).
killboss
 Аватар для killboss
57 / 57 / 2
Регистрация: 18.12.2010
Сообщений: 127
20.01.2011, 21:06     Необходимо изменить программу в C++ #10
Нет разбей pow() на части и пересмотри ещё раз математику.
sl1m_dogg
0 / 0 / 0
Регистрация: 20.01.2011
Сообщений: 18
20.01.2011, 22:49  [ТС]     Необходимо изменить программу в C++ #11
Переписал четвертую часть в таком виде:
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
//.4.........................
Fourth:  cout <<" __4__"<<endl;
     X[k+1][0]=Y[n+1][0];    cout <<"X["<<(k+1)<<"] = ("<<X[k+1][0];
     X[k+1][1]=Y[n+1][1];    cout <<";"<<X[k+1][1]<<")"<<endl;
//...........a...............
double n1 = X[k+1][0]-X[k][0];
double n2 = X[k+1][1]-X[k][1];
double n3 = pow(n1,2)+pow(n2,2);
if(pow(n3,1/2) < E)
{
     cout <<"X(k+1) - X(k) = "<<n3;
 
         Xz[0][0]=X[k+1][0];
         Xz[0][1]=X[k+1][1];
     }
//..................b.....................
 n1 = X[k+1][0]-X[k][0];
 n2 = X[k+1][1]-X[k][1];
 n3 = pow(n1,2)+pow(n2,2);
if(pow(n3,1/2) < E)
{
 
             d[0][0]=d[n][0]=Y[n+1][0]-Y[1][0];    cout <<"d0 = ("<<d[0][0];
             d[0][1]=d[n][1]=Y[n+1][1]-Y[1][1];    cout <<";"<<d[0][1]<<")"<<endl;
 
             d[I[s]][0]=d[I[s]+1][0];
             d[I[s]][1]=d[I[s]+1][1];
             if( (d[1][0]*d[2][1]-d[1][1]*d[2][0]) != 0){
                cout <<"rang(D1,...,Dn) = n"<<endl;
                Y[0][0]=X[k+1][0];  cout <<"Y0 = ("<<Y[0][0]<<";";
                Y[0][1]=X[k+1][1];  cout <<Y[0][1]<<")\n";
                k=k+1;
                I[s]=0;
                goto TWO;
             }
             if( (d[1][0]*d[2][1]-d[1][1]*d[2][0]) == n){
                cout <<"rang(D1,...,Dn) < n"<<endl;
                Y[0][0]=X[k+1][0];
                Y[0][1]=X[k+1][1];
                k=k+1;
                I[s]=0;
                goto TWO;
             }
         }
     cout <<"\nXz("<<Xz[0][0]<<";"<<Xz[0][1]<<")"<<endl;
     getch();}
минимальное значение функции не меняется (не верное)

Добавлено через 1 час 28 минут
killboss, можешь посмотреть что происходить с переменной t. По моему она считается только на первой итерации, а на всех остальных постоянно равна 0.
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
20.01.2011, 23:02     Необходимо изменить программу в C++
Еще ссылки по теме:

необходимо изменить программу C++
C++ Необходимо написать программу в C++
C++ Работа с файлом. Необходимо изменить программу

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

Или воспользуйтесь поиском по форуму:
killboss
 Аватар для killboss
57 / 57 / 2
Регистрация: 18.12.2010
Сообщений: 127
20.01.2011, 23:02     Необходимо изменить программу в C++ #12
Использование GOTO - плохой стиль программирования. Достаточно трудно понять алгоритм, в котором использунтся оператор безусловного перехода!
Yandex
Объявления
20.01.2011, 23:02     Необходимо изменить программу в C++
Ответ Создать тему
Опции темы

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