Лабораторная работа - Безусловная многомерная оптимизация второго порядка (Вариант 10)

Подождите немного. Документ загружается.

f2[2][2]=dfx2x2(x[k].x1,x[k].x2);

N2=N2+3;

m2:

d=f2[1][1]*f2[2][2]-pow(f2[1][2],2);

d1=f1_1*f2[2][2]-f1_2*f2[1][2];

d2=f2[1][1]*f1_2-f2[2][1]*f1_1;

p1=d1/d;

p2=d2/d;

//opredel A,B

y1=f(x[k].x1,x[k].x2);

i=0;

m5:

y2=f(x[k].x1-((i+1)*da)*p1,x[k].x2-((i+1)*da)*p2);

N0++;

if (y2<y1) {i++;y1=y2;goto m5;}

else {B=(1+i)*da;

if (i>0) A=(i-1)*da;

else A=0;

}

//cout<<"A="<<A<<" B="<<B<<endl;

d0=eps/100;

{a=(A+B)/2.0;

a1=a-d0;

a2=a+d0;

f1=f(x[k].x1-a1*p1,x[k].x2-a1*p2);

f3=f(x[k].x1-a2*p1,x[k].x2-a2*p2);

N0=N0+2;

if (f1<f3) B=a2;

else A=a1;

}

while (B-A>2*eps);

amin=(A+B)/2.0;

x[k+1].x1=x[k].x1-amin*p1;

x[k+1].x2=x[k].x2-amin*p2;

fout<<x[k+1].x1<<"; "<<x[k+1].x2<<endl;

cout<<(k+1)<<") "<<x[k+1].x1<<" "<<x[k+1].x2<<endl;

f1_1=dfx1(x[k+1].x1,x[k+1].x2);

f1_2=dfx2(x[k+1].x1,x[k+1].x2);

N1=N1+2;

F=f1_1*f1_1+f1_2*f1_2;

dlina=sqrt(F);

k++;

if (dlina>eps) goto m2;

x_=x[k];

y_=f(x_.x1,x_.x2);

N=N1+N2+N0;

fout.close();

cout<<"Result:\n";

cout<<"x_(x1,x2) = x_( "<<x_.x1<<","<<x_.x2<<" )"<<endl;

cout<<"y_ = "<<y_<<endl;

cout<<"Kol-vo icpitanii: N0 = "<<N0<<endl;

cout<<"Kol-vo icpitanii: N1 = "<<N1<<endl;

cout<<"Kol-vo icpitanii: N2 = "<<N2<<endl;

cout<<"Kol-vo icpitanii: N = "<<N<<endl;

cout<<"k="<<k<<endl;

}

void s_modif_new()

{

const int L=1500;

znach_x x[L];

znach_x x_;

long double eps,y_,f1_1,f1_2,d,d1,d2,p1,p2,dlina,F;

long double f2[2][2];

int N1,N2,k,N,m=9;

ofstream fout;

cout<<"Vvedite bazovuyu tochku:"<<endl;

cout<<"x0(x1,x2) ";cout<<endl;

cout<<"x1 = ";

cin>>x[0].x1;

cout<<"x2 = ";

cin>>x[0].x2;

eps = 0.0004;

fout.open("result7.txt");

if(fout.fail()){ cout<<"result7.txt "<<"Mistake!";

getch(); exit(1);}

f1_1=dfx1(x[0].x1,x[0].x2);

f1_2=dfx2(x[0].x1,x[0].x2);

/*cout<<"x[k]:"<<endl<<x[0].x1<<endl<<x[0].x2<<endl;

cout<<"f1_1="<<f1_1<<endl;

cout<<"f1_2="<<f1_2<<endl;

k=0;N1=0;N2=0;

m1:

f2[1][1]=dfx1x1(x[k].x1,x[k].x2);

f2[1][2]=dfx1x2(x[k].x1,x[k].x2);

f2[2][1]=f2[1][2];

f2[2][2]=dfx2x2(x[k].x1,x[k].x2);

/*cout<<"f211="<<f2[1][1]<<endl;

cout<<"f212="<<f2[1][2]<<endl;

cout<<"f221="<<f2[2][1]<<endl;

cout<<"f222="<<f2[2][2]<<endl;

N2=N2+3;

m2:

d=f2[1][1]*f2[2][2]-pow(f2[1][2],2);

d1=f1_1*f2[2][2]-f1_2*f2[1][2];

d2=f2[1][1]*f1_2-f2[2][1]*f1_1;

/*cout<<"x[k]:"<<endl<<x[k].x1<<endl<<x[k].x2<<endl;

cout<<"d="<<d<<endl;

cout<<"d1="<<d1<<endl;

cout<<"d2="<<d2<<endl;

p1=d1/d;

p2=d2/d;

x[k+1].x1=x[k].x1-p1;

x[k+1].x2=x[k].x2-p2;

fout<<x[k+1].x1<<"; "<<x[k+1].x2<<endl;

cout<<(k+1)<<") "<<x[k+1].x1<<endl<<x[k+1].x2<<endl;

f1_1=dfx1(x[k+1].x1,x[k+1].x2);

f1_2=dfx2(x[k+1].x1,x[k+1].x2);

N1=N1+2;

F=f1_1*f1_1+f1_2*f1_2;

dlina=sqrt(F);

//cout<<"F="<<F<<endl;

k++;

if (dlina>eps)

{if (k%m==0)goto m1;

else goto m2; }

x_=x[k];

y_=f(x_.x1,x_.x2);

N=N1+N2;

fout.close();

cout<<"Result:\n";

cout<<"x_(x1,x2) = x_( "<<x_.x1<<","<<x_.x2<<" )"<<endl;

cout<<"y_ = "<<y_<<endl;

cout<<"Kol-vo icpitanii: N1 = "<<N1<<endl;

cout<<"Kol-vo icpitanii: N2 = "<<N2<<endl;

cout<<"Kol-vo icpitanii: N = "<<N<<endl;

cout<<"k="<<k<<endl;

}

void s_modif_new_rafs_drobl()

{

const int L=1500;

znach_x x[L];

znach_x x_;

double eps,y_,f1_1,f1_2,d,d1,d2,p1,p2,dlina,F,delta,a,y,fp;

double f2[2][2];

int N1,N2,k,N,N0,m=2;

ofstream fout;

cout<<"Vvedite bazovuyu tochku:"<<endl;

cout<<"x0(x1,x2) ";cout<<endl;

cout<<"x1 = ";

cin>>x[0].x1;

cout<<"x2 = ";

cin>>x[0].x2;

eps=0.0004;

cout<<"Vvedite delta':\n";

cout<<"delta = "; cin>>delta;

fout.open("result8.txt");

if(fout.fail()){ cout<<"result8.txt "<<"Mistake!";

getch(); exit(1);}

f1_1=dfx1(x[0].x1,x[0].x2);

f1_2=dfx2(x[0].x1,x[0].x2);

k=0;N1=0;N2=0;a=1; N0=0;

m2:

f2[1][1]=dfx1x1(x[k].x1,x[k].x2);

f2[1][2]=dfx1x2(x[k].x1,x[k].x2);

f2[2][1]=f2[1][2];

f2[2][2]=dfx2x2(x[k].x1,x[k].x2);

N2=N2+3;

m6:

d=f2[1][1]*f2[2][2]-pow(f2[1][2],2);

d1=f1_1*f2[2][2]-f1_2*f2[1][2];

d2=f2[1][1]*f1_2-f2[2][1]*f1_1;

p1=d1/d;

p2=d2/d;

y=f(x[k].x1,x[k].x2);

N0++;

fp=f1_1*p1+f1_2*p2; //scal.

m1:

if ((f(x[k].x1-a*p1,x[k].x2-a*p2)-y)>(-delta*a*fp)) {a=0.5*a;

N0++;goto m1;}

x[k+1].x1=x[k].x1-a*p1;

x[k+1].x2=x[k].x2-a*p2;

fout<<x[k+1].x1<<"; "<<x[k+1].x2<<endl;

cout<<(k+1)<<") "<<x[k+1].x1<<" "<<x[k+1].x2<<endl;

f1_1=dfx1(x[k+1].x1,x[k+1].x2);

f1_2=dfx2(x[k+1].x1,x[k+1].x2);

N1=N1+2;

F=f1_1*f1_1+f1_2*f1_2;

dlina=sqrt(F);

k++;

if (dlina>eps) { if (k%m==0) goto m2;

else goto m6; }

x_=x[k];

y_=f(x_.x1,x_.x2);

N=N1+N2+N0;

fout.close();

cout<<"Result:\n";

cout<<"x_(x1,x2) = x_( "<<x_.x1<<","<<x_.x2<<" )"<<endl;

cout<<"y_ = "<<y_<<endl;

cout<<"Kol-vo icpitanii: N0 = "<<N0<<endl;

cout<<"Kol-vo icpitanii: N1 = "<<N1<<endl;

cout<<"Kol-vo icpitanii: N2 = "<<N2<<endl;

cout<<"Kol-vo icpitanii: N = "<<N<<endl;

cout<<"k="<<k<<endl;

}

void s_modif_new_rafs_optim()

{

const int L=1500;

znach_x x[L];

znach_x x_;

double eps,y_,f1_1,f1_2,d,d1,d2,p1,p2,dlina,F,y1,y2,A,B;

double a1,a2,amin,da,d0,a,f1,f3;

double f2[2][2];

int i,N1,N2,k,N,N0,m=2;

ofstream fout;

cout<<"Vvedite bazovuyu tochku:"<<endl;

cout<<"x0(x1,x2) ";cout<<endl;

cout<<"x1 = ";

cin>>x[0].x1;

cout<<"x2 = ";

cin>>x[0].x2;

cout<<"Vvedite da:\n";

cout<<"da = "; cin>>da;

eps=0.0004;

fout.open("result9.txt");

if(fout.fail()){ cout<<"result9.txt "<<"Mistake!";

getch(); exit(1);}

f1_1=dfx1(x[0].x1,x[0].x2);

f1_2=dfx2(x[0].x1,x[0].x2);

k=0;N1=0;N2=0;N0=0;

m2:

f2[1][1]=dfx1x1(x[k].x1,x[k].x2);

f2[1][2]=dfx1x2(x[k].x1,x[k].x2);

f2[2][1]=f2[1][2];

f2[2][2]=dfx2x2(x[k].x1,x[k].x2);

N2=N2+3;

m6:

d=f2[1][1]*f2[2][2]-pow(f2[1][2],2);

d1=f1_1*f2[2][2]-f1_2*f2[1][2];

d2=f2[1][1]*f1_2-f2[2][1]*f1_1;

p1=d1/d;

p2=d2/d;

//opredel A,B

y1=f(x[k].x1,x[k].x2);

i=0;

m5:

y2=f(x[k].x1-((i+1)*da)*p1,x[k].x2-((i+1)*da)*p2);

N0++;

if (y2<y1) {i++;y1=y2;goto m5;}

else {B=(1+i)*da;

if (i>0) A=(i-1)*da;

else A=0;

}

//cout<<"A="<<A<<" B="<<B<<endl;

d0=eps/100;

{a=(A+B)/2.0;

a1=a-d0;

a2=a+d0;

f1=f(x[k].x1-a1*p1,x[k].x2-a1*p2);

f3=f(x[k].x1-a2*p1,x[k].x2-a2*p2);

N0=N0+2;

if (f1<f3) B=a2;

else A=a1;

}

while (B-A>2*eps);

amin=(A+B)/2.0;

x[k+1].x1=x[k].x1-amin*p1;

x[k+1].x2=x[k].x2-amin*p2;

fout<<x[k+1].x1<<"; "<<x[k+1].x2<<endl;

cout<<(k+1)<<") "<<x[k+1].x1<<" "<<x[k+1].x2<<endl;

f1_1=dfx1(x[k+1].x1,x[k+1].x2);

f1_2=dfx2(x[k+1].x1,x[k+1].x2);

N1=N1+2;

F=f1_1*f1_1+f1_2*f1_2;

dlina=sqrt(F);

k++;

if (dlina>eps) if (k%m==0) goto m2;

else goto m6;

x_=x[k];

y_=f(x_.x1,x_.x2);

N=N1+N2+N0;

fout.close();

cout<<"Result:\n";

cout<<"x_(x1,x2) = x_( "<<x_.x1<<","<<x_.x2<<" )"<<endl;

cout<<"y_ = "<<y_<<endl;

cout<<"Kol-vo icpitanii: N0 = "<<N0<<endl;

cout<<"Kol-vo icpitanii: N1 = "<<N1<<endl;

cout<<"Kol-vo icpitanii: N2 = "<<N2<<endl;

cout<<"Kol-vo icpitanii: N = "<<N<<endl;

cout<<"k="<<k<<endl;

}

void main()

{int j;

while (1)

{clrscr();

cout<<endl;

cout<<"1.Method Newton\n";

cout<<"2.Metod Newton-Rafson s drobleniem\n";

cout<<"3.Metod Newton-Rafson optim\n";

cout<<"4.1 modification Newton\n";

cout<<"5.1 modification Newton-Rafson s drobleniem\n";

cout<<"6.1 modification Newton-Rafson optim\n";

cout<<"7.2 modification Newton\n";

cout<<"8.2 modification Newton-Rafson s drobleniem\n";

cout<<"9.2 modification Newton-Rafson optim\n";

cout<<"0.Exit\n";

cout<<endl;

cout<<"Your choice 0-9: ";

cin>>j;

switch(j)

{

case 1:newton();getch(); break;

case 2:newton_rafs_drobl();getch(); break;

case 3:newton_rafs_optim();getch(); break;

case 4:f_modif_new();getch(); break;

case 5:f_modif_new_rafs_drobl();getch(); break;

case 6:f_modif_new_rafs_optim();getch(); break;

case 7:s_modif_new();getch(); break;

case 8:s_modif_new_rafs_drobl();getch(); break;

case 9:s_modif_new_rafs_optim();getch(); break;

case 0: cout<<"Over"; getch();exit(0);

default:{cout<<"Net takogo punkta!";exit(1);

}

getch();}

}

Сравнение результатов

Название

метода

xmin ymin Общее кол-во

испытаний

Кол-во итераций

x1 x2

1. Поиск по

образцу

-2,712158 -0,097412 -54,722806 99 28

2. Деформирмир.

симплекс

-2,713434 -0,101445 -54,722682 50 24

3.Регул. симплекс -2,712158 -0,097656 -54,722806 50 22

4. Градиентный с

дроблением шага

-2,71213 -0,097677 -54,722806 86 27

5. Наискорейший

спуск

-2,717127 -0,097659 -54,722806 71 22

6. Гаусс-Зейдель -2,71213 -0,0976367 -54,7228 278 3

7.Эвристический -2,71213 -0,0976367 -54,7228 506 109

Название метода

xmin ymin Общее кол-во

испытаний

Кол-во итераций

x1 x2

1. Ньютона -2,71213 -0,0976373 -54,7228 385 77

2. Ньютона(1я)

нач. приближ. (-2,7;-0,09)

-2,71213 -0,0976364 -54,7228 25 11

3. Ньютона(2я)

нач. приближ. (-2,7;-0,09)

m=9

-2,71214 -0,0976304 -54,7228 26 10

4. Ньютона-Рафсона

с дробл.

-2.71214 -0,0976427 -54,7228 224 37

5. Ньютона-Рафсона

с дробл.(1я)

-2,71213 -0,0976771 -54,7228 80 24

6. Ньютона-Рафсона

с дробл.(2я)

-2,71212 -0,097647 -54,7228 1253 277

7. Ньютона-Рафсона

с оптим.

-2,71213 -0,0976377 -54,7228 88 3

8. Ньютона-Рафсона

с оптим.(1я)

-2,71214 -0,0976128 -54,7228 98 5

9. Ньютона-Рафсона

с оптим.(2я)

-2,71213 -0,0976369 -54,7228 76 3

Вывод

На основе полученных результатов, приведённых в таблице, можно сделать вывод,

что наилучшим методом многомерной безусловной оптимизации с точки зрения

количества итераций является метод Гаусса-Зейделя, а лучший метод с точки зрения

количества экспериментов – это метод деформируемого симплекса.

Наилучшим методом многомерной безусловной оптимизации второго порядка с

точки зрения числа испытаний является первая модификация метода Ньютона, а с

точки зрения количества итераций - метод Ньютона-Рафсона с оптимальным шагом

и его вторая модификация.