Добавил:
Studfiles2
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:Различные алгоритмы минимизации, реализованные в С++ / all / FIBDAV2
.CPP#include <conio.h>
#include <iostream.h>
#include <math.h>
struct Vector
{
double x1,x2;
Vector()
{
x1=0;
x2=0;
}
Vector(double ax1,double ax2)
{
x1=ax1;
x2=ax2;
}
const Vector operator * (const double &vec) const
{
return Vector(x1*vec,x2*vec);
}
const void operator = (const Vector &vec)
{
x1=vec.x1;
x2=vec.x2;
}
const double operator * (const Vector &vec) const
{
return (x1*vec.x1+x2*vec.x2);
}
const Vector operator + (const Vector &vec) const
{
return Vector(x1+vec.x1,x2+vec.x2);
}
const Vector operator - (const Vector &vec) const
{
return Vector(x1-vec.x1,x2-vec.x2);
}
};
double function(Vector vec)
{
return ((vec.x1-2)*(vec.x1-2)*(vec.x1-2)*(vec.x1-2)+(vec.x1-2*vec.x2)*(vec.x1-2*vec.x2));
}
double f_dx1(Vector vec)
{
return (4*(vec.x1-2)*(vec.x1-2)*(vec.x1-2)+2*(vec.x1-2*vec.x2));
}
double f_dx2(Vector vec)
{
return (-4*(vec.x1-2*vec.x2));
}
double proizv(Vector vec,Vector p)
{
Vector grad;
grad=Vector(f_dx1(vec),f_dx2(vec));
return (grad*p);
}
void Swann4(Vector& x1,double& a,double& b,Vector& p)
{
double pr;
Vector x2;
pr=proizv(x1,p);
if(pr>0) p=Vector(0,0)-p;
b=fabs(function(x1)/pr-1);
if(b>1) b=1;
x2=x1+p*b;
while(proizv(x1,p)*proizv(x2,p)>0)
{
a=b;
b*=2;
x1=x2;
x2=x1+p*b;
}
}
void Davidon(double& a,double& b,Vector& x1,Vector& x2,Vector p)
{
double e=0.0001,w,z,ar,k;
int n=1;
Vector r;
z=proizv(x1,p)+proizv(x2,p)+3*(function(x1)-function(x2))/(b-a);
w=sqrt(z*z-proizv(x1,p)*proizv(x2,p));
ar=b*(w-proizv(x1,p)+z)/(2*w-proizv(x1,p)+proizv(x2,p));
r=x1+p*ar;
k=proizv(r,p);
while(fabs(k)>e)
{
if(k<0)
{
a=ar;
x1=r;
}
else
{
b=ar;
x2=r;
}
z=proizv(x1,p)+proizv(x2,p)+3*(function(x1)-function(x2))/(b-a);
w=sqrt(z*z-proizv(x1,p)*proizv(x2,p));
ar=b*(w-proizv(x1,p)+z)/(2*w-proizv(x1,p)+proizv(x2,p));
r=x1+p*ar;
k=proizv(r,p);
n++;
}
cout<<r.x1<<endl<<r.x2;
}
void main()
{
clrscr();
Vector p,x1,x2;
double a=0,b=0;
x1=Vector(0,3);
p=Vector(1,0);
Swann4(x1,a,b,p);
Davidon(a,b,x1,x2,p);
getch();
}
Соседние файлы в папке all