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

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
Рейтинг: Рейтинг темы: голосов - 18, средняя оценка - 4.94
sl1m_dogg
0 / 0 / 0
Регистрация: 20.01.2011
Сообщений: 18
#1

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

20.01.2011, 19:09. Просмотров 2685. Ответов 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
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
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
20.01.2011, 19:09
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Необходимо изменить программу в C++ (C++):

необходимо изменить программу - C++
#include &lt;iostream&gt; int main() { const int N = 5; using namespace std; double n; int count = N; float...

Работа с файлом. Необходимо изменить программу - C++
Мне преподаватель сказал, что программа сделана правильно. Теперь необходимо ее переделать так, чтобы после создания файла выполнялась...

Что необходимо изменить в программе? - C++
Подскажите пожалуйста что необходимо сделать чтобы высчитывался стаж, как связать два класса? #include &quot;stdafx.h&quot; #include &quot;cstdlib&quot; ...

Необходимо изменить код программы. - C++
Есть программа для измерения напряжения до 25.5В и тока до 2.55А , а мне необходимо измерять напряжение и ток до 50В и 5А соответственно....

Необходимо изменить размеры окна консоли - C++
Приветствую. Необходимо изменить размеры окна консоли, допустим, на 50 * 50 средствами WinAPI, поэтому вот такой код system...

Решение интеграла (необходимо немного изменить код) - C++
есть код, как я думал правильный (он компилируется без ошибок)! Но в итоге выдает решение интеграла равно 0, что я думаю не есть правильно,...

11
killboss
60 / 60 / 2
Регистрация: 18.12.2010
Сообщений: 127
20.01.2011, 20:20 #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() никаких выражений не должно быть.
1
sl1m_dogg
0 / 0 / 0
Регистрация: 20.01.2011
Сообщений: 18
20.01.2011, 20:33  [ТС] #3
А не могли бы вы продемонстрировать как должен выглядеть этот if? А то я что то не особо понял
0
killboss
60 / 60 / 2
Регистрация: 18.12.2010
Сообщений: 127
20.01.2011, 20:42 #4
Жди минутку...
0
sl1m_dogg
0 / 0 / 0
Регистрация: 20.01.2011
Сообщений: 18
20.01.2011, 20:43  [ТС] #5
Может я не совсем правильно описал. Программа правильно работает для первой функции (которая в коде прописана), она выводит правильное значение т. (5;6). Но когда я в return прописываю значение второй функции программа выводит неверное значение. Первоначальную точку вводил (0;0).
0
killboss
60 / 60 / 2
Регистрация: 18.12.2010
Сообщений: 127
20.01.2011, 20:46 #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)
И так сделать по всему алгоритму!!!
1
sl1m_dogg
0 / 0 / 0
Регистрация: 20.01.2011
Сообщений: 18
20.01.2011, 20:54  [ТС] #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 ошибок для второй функции.
0
killboss
60 / 60 / 2
Регистрация: 18.12.2010
Сообщений: 127
20.01.2011, 20:56 #8
На что ругается компилятор? Скинь сообщение...
1
sl1m_dogg
0 / 0 / 0
Регистрация: 20.01.2011
Сообщений: 18
20.01.2011, 21:00  [ТС] #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;}
то программа выдает неверную точку минимума (правильная точка уже известна).
0
killboss
60 / 60 / 2
Регистрация: 18.12.2010
Сообщений: 127
20.01.2011, 21:06 #10
Нет разбей pow() на части и пересмотри ещё раз математику.
1
sl1m_dogg
0 / 0 / 0
Регистрация: 20.01.2011
Сообщений: 18
20.01.2011, 22:49  [ТС] #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.
0
killboss
60 / 60 / 2
Регистрация: 18.12.2010
Сообщений: 127
20.01.2011, 23:02 #12
Использование GOTO - плохой стиль программирования. Достаточно трудно понять алгоритм, в котором использунтся оператор безусловного перехода!
1
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
20.01.2011, 23:02
Привет! Вот еще темы с ответами:

Работа со стеком. Если массив нужен float, а не int, что необходимо изменить - C++
#include &lt;iostream&gt; #include &lt;stdlib.h&gt; using namespace std; // ****************СТЕК НА ОСНОВЕ МАССИВА*************** const int N =...

Необходимо изменить строку так, чтобы длина строки стала равна заданной длине - C++
Дана строка текста, в котором нет начальных и конечных пробелов. Необходимо изменить ее так, чтобы длина строки стала равна заданной длине...

Необходимо разобрать программу - C++
Добрый день . встретились такие операторы как&quot;ch+=(char) 32&quot; и &quot;ch-=(char) 32&quot; Объясните принцип работы этих операторов. В моей книге об...

Необходимо подправить программу - C++
я написал программу: #include &lt;iostream&gt; #include &quot;liquid.h&quot; #include &quot;SpNapitki.h&quot; using namespace std; void...


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

Или воспользуйтесь поиском по форуму:
12
Yandex
Объявления
20.01.2011, 23:02
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru