#include "iostream.h"
#include "matrix.h"
double f(vector& x);
void g(vector& x,vector& y);//专为LMF方法所用
void main(void)
{
vector x(5);
double minf;
int i;
cout.precision(15);
x.Init(10.0,20.0,30.0,40.0,50.0);
minf=x.FR(f);
cout<<"FR方法计算结果:"<<endl;
for(i=1;i<=5;i++)
cout<<"x["<<i<<"]="<<x[i]<<endl;
cout<<"minf = "<<f(x)<<endl;
cout<<endl;
x.Init(10.0,20.0,30.0,40.0,50.0);
minf=x.DFP(f);
cout<<"DFP方法计算结果:"<<endl;
for(i=1;i<=5;i++)
cout<<"x["<<i<<"]="<<x[i]<<endl;
cout<<"minf = "<<f(x)<<endl;
cout<<endl;
x.Init(10.0,20.0,30.0,40.0,50.0);
minf=x.BFGS(f);
cout<<"BFGS方法计算结果:"<<endl;
for(i=1;i<=5;i++)
cout<<"x["<<i<<"]="<<x[i]<<endl;
cout<<"minf = "<<f(x)<<endl;
cout<<endl;
x.Init(10.0,20.0,30.0,40.0,50.0);
minf=x.TR(f);
cout<<"TR方法计算结果:"<<endl;
for(i=1;i<=5;i++)
cout<<"x["<<i<<"]="<<x[i]<<endl;
cout<<"minf = "<<f(x)<<endl;
cout<<endl;
x.Init(10.0,20.0,30.0,40.0,50.0);
minf=x.Newton(f);
cout<<"Newton方法计算结果:"<<endl;
for(i=1;i<=5;i++)
cout<<"x["<<i<<"]="<<x[i]<<endl;
cout<<"minf = "<<f(x)<<endl;
cout<<endl;
x.Init(10.0,20.0,30.0,40.0,50.0);
minf=x.Powell(f);
cout<<"Powell方法计算结果:"<<endl;
for(i=1;i<=5;i++)
cout<<"x["<<i<<"]="<<x[i]<<endl;
cout<<"minf = "<<f(x)<<endl;
cout<<endl;
x.Init(10.0,20.0,30.0,40.0,50.0);
minf=x.LMF(g,5);
cout<<"LMF方法计算结果:"<<endl;
for(i=1;i<=5;i++)
cout<<"x["<<i<<"]="<<x[i]<<endl;
cout<<"minf = "<<f(x)<<endl;
}
double f(vector& x)
{
double v=0;
for(int i=1;i<=4;i++)
v+=(x[i+1]-x[i]*x[i])*(x[i+1]-x[i]*x[i]);
v*=100;
v+=(1-x[1])*(1-x[1]);
return v;
}
void g(vector& x,vector& y)
{
y[1]=10*(x[2]-x[1]*x[1]);
y[2]=10*(x[3]-x[2]*x[2]);
y[3]=10*(x[4]-x[3]*x[3]);
y[4]=10*(x[5]-x[4]*x[4]);
y[5]=1-x[1];
}
数值与非数值分析VC++类库
5星 · 超过95%的资源 需积分: 31 31 浏览量
2015-11-02
19:36:16
上传
评论 3
收藏 769KB RAR 举报 
数值与非数值分析VC++类库 （280个子文件） 
5.33.cpp 2KB 5.35.cpp 1KB 7.09.cpp 1KB 6.14.cpp 1KB 7.15.cpp 1KB 5.20.cpp 1KB 5.21.cpp 1KB 5.08.cpp 1002B 4.22.cpp 986B 7.19.cpp 960B 5.19.cpp 955B 5.18.cpp 954B 5.22.cpp 943B 5.17.cpp 848B 5.28.cpp 842B 7.61.cpp 818B 5.25.cpp 813B 7.64.cpp 812B 5.34.cpp 811B 7.59.cpp 798B 7.03.cpp 793B 7.48.cpp 786B 5.06.cpp 734B 7.11.cpp 730B 4.12.cpp 724B 6.05.cpp 701B 6.01.cpp 700B 7.20.cpp 692B 4.09.cpp 689B 7.49.cpp 686B 4.10.cpp 685B 7.01.cpp 685B 6.04.cpp 663B 7.18.cpp 642B 7.14.cpp 641B 7.51.cpp 619B 4.13.cpp 616B 5.04.cpp 612B 7.10.cpp 599B 5.32.cpp 597B 7.67.cpp 586B 4.19.cpp 578B 6.11.cpp 574B 6.03.cpp 571B 6.02.cpp 565B 7.37.cpp 550B 4.21.cpp 545B 7.25.cpp 539B 5.29.cpp 537B 4.20.cpp 535B 4.16.cpp 533B 6.06.cpp 520B 5.07.cpp 517B 5.02.cpp 510B 7.02.cpp 510B 4.11.cpp 508B 5.27.cpp 507B 5.23.cpp 504B 5.03.cpp 486B 7.57.cpp 483B 5.05.cpp 478B 1.04.cpp 472B 4.15.cpp 468B 7.28.cpp 458B 5.26.cpp 458B 7.50.cpp 456B 7.35.cpp 456B 7.43.cpp 455B 6.12.cpp 453B 2.53.cpp 435B 2.04.cpp 415B 2.39.cpp 413B 4.23.cpp 412B 6.08.cpp 410B 6.07.cpp 402B 5.15.cpp 401B 4.04.cpp 397B 6.13.cpp 390B 5.31.cpp 389B 7.46.cpp 386B 5.30.cpp 385B 4.07.cpp 384B 1.01.cpp 384B 5.13.cpp 382B 7.63.cpp 373B 4.17.cpp 372B 5.14.cpp 367B 2.10.cpp 363B 4.01.cpp 347B 6.09.cpp 347B 7.44.cpp 338B 5.09.cpp 332B 7.47.cpp 331B 2.01.cpp 328B 7.06.cpp 325B 7.53.cpp 325B 5.16.cpp 324B 4.18.cpp 319B 3.11.cpp 311B 7.62.cpp 307B共 280 条
- 1
- 2
- 3
- 1
- 2
前往页