无约束最优化方法可变单纯形法(simple)nelder-mead.docxVIP

v1.0可编辑可修改 v1.0可编辑可修改 PAGE PAGE # 无约束最优化方法可变单纯形法(simplex ) Nelder-Mead 可爱的馒头 本程序是用C + +编写的，从编写的算例来看，应该是没有问题的。所采用的原理和步骤是参考华南理工 大学岀版社蒋金山等编写的最优化计算方法第8章第三节可变单纯形法。欢迎各位批评指正。 #include #include #include int i,j; double d[3][100]={{0,1,0,0},{0,0,1,0},{0,0,0,1}},f[100];//d[][] 为单纯形的顶点，本算例中未知数 个数为3，则顶点个数为4 double g,h,l,q,s=1,t=2,u=,v=,y=0;//s 为反射系数，t为扩展系数，u为压缩系数，v为允许误差 int o,F,r,D,e,lj=0,N=4;//N 为顶点的个数，o为最大值点的位置，F为最小值点的位置，r为次大值点的 位置 void function1(int e)// 求解函数 f[e] { f[e]=(d[0] [e]-3)*(d[0] [e]-3)+2*(d[1] [ e]+2)*(d[1] [ e]+2)+(d[2 ][e]-4)*(d[2 ][ e]-4);// 函数为 f=(x1-3)A2+2(x2+2)A2+(x3-4)A2, 求其最小值 } void function2() { while((++lj)100)// 最大迭代次数 { for(i=0,g=f[i];iN-1;i++)〃 求最大值点 { if(gf[i+1]) { g=f[i+1];o=i+1; } else if(i==0) o=i; } for(i=0,h=f[i];iN-1;i++)〃 求最小值点 { if(hf[i+1]) { h=f[i+1];F=i+1; } else if(i==0) F=i; } for(i=0,l=f[i];iN-1;i++)〃 求次大值点 { if(i==oi==0) { l=f[i+1]; r=i+1; continue; } if(lf[i+1](i+1)!=o) { l=f[i+1];r=i+1; } else if(i==0) r=i; } for(i=0;iN-1;i++)〃 求除最大值点，其余点的形心 { d[i][N]=0; for(j=0;jN;j++) { if(j!=o) { d【i]【N]=d【i]【N]+d【i][j]; } } d[i][N]=d[i][N]/(N-1);〃 平均值 } functionl(N); for(i=0;iN-1;i++) { d[i][N+1]=d[i][N]+(d[i][N]-d[i][o])*s;〃 反射 } function1(N+1); if(f[N+1]vf[F]) { for(i=0;iN-1;i++) { d[i][N+2]=d[i][N]+(d[i][N+1]-d[i][N])*t;〃 扩展 } function1(N+2); if(f[N+2]vf[N+1]) { for(i=0;iN-1;i++) { d【i]【o]=d【i]【N+2]; } f[o]=f[N+2]; y=0; for(i=0;iN;i++) { y=y+(f[i]-f[N])*(f[i]-f[N]); } y=sqrt(y/(N*); if(yv) { for(i=0;iN-1;i++) { printf(d[%d][%d]=%f,i,F,d[i][F]); printf(\n); } printf(f[%d]=%f,F,f[F]); break; } else { function2(); break; } } else { for(i=0;iN-1;i++) { d【i]【o]=d【i]【N+1]; } f[o]=f[N+1]; y=0; for(i=0;iN;i++) { y=y+(f[i]-f[N])*(f[i]-f[N]); } y=sqrt(y/(N*); v1.0可编辑可修改 v1.0可编辑可修改 PAGE PAGE # v1.0可编辑可修改 v1.0可编辑可修改 { PAGE { PAGE # if(yw) { for(i=0;iN-1;i++) { printf(d[%d][%d]=%f,i,F,d[i][F]); printf(\n); } printf(f[%d]=%f,F,f[F]); break; } else { function2(); break; } } } if(f[N+1]=f[F ] f[N+1]v=f[r]) { for(i=0;iN-1;i++) d【i]【o]=d【i]【