第五章约束优化方法2讲课文档.pptVIP

3.4复合形法的迭代步骤(1)构造初始复合形；(2)计算各顶点的函数值F(X(j))，j=1,2,….,K。选出好点X(L)和坏点X(H)；(3)计算坏点外的其余各顶点的中心点X(0)。第28页，共72页。(4)计算映射点X(R)检查X(R)是否在可行域内。若X(R)为非可行点，将映射系数减半后再按上式改变映射点，直到X(R)进入可行域内为止。(5)构造新的复合形计算映射点的函数值F(X(R))，并与坏点的函数值F(X(H))比较，可能存在两种情况：1）映射点优于坏点F(X(R))F(X(H))在此情况，用X(R)代替X(`H)，构成新的复合形。第29页，共72页。若经过多次的映射系数减半，仍不能使映射点由于坏点，则说明该映射方向不利，此时，应改变映射方向，取对次坏点的映射。再转回本步骤的开始处，直到构成新的复合形。2）映射点次于坏点F(X(R))F(X(H))这种情况由于映射点过远引起的，减半映射系数，若有F(X(R))F(X(H))，这又转化为第一种情况。第30页，共72页。3.5判断终止条件1)各顶点与好点函数值之差的均方根值小于误差限，即2）各顶点与好点的函数值之差的平方和小于误差限，即3）各顶点与好点函数值差的绝对值之和小于误差限，即如果不满足终止迭代条件，则返回步骤2继续进行下一次迭代；否则，可将最后复合形的好点X(L)及其函数值F(X(L))作为最优解输出。第31页，共72页。比较复合形各顶点的函数值，找出好点XL，坏点XHXH=XRα=0.5α找出次坏点XSH，XH=XSH满足终止条件？X*=XL,F*=F(XL)结束3.6流程图是否给定K，δ，α，ε，ai,bii=1,2,…n产生初始复合形顶点Xj,j=1,2,…,K计算复合形各顶点的函数值F(Xj),j=1,2,…,K是是是否否否XR∈DFRF(XH)第32页，共72页。3.7复合形法的Matlab程序程序清单function[xo,fo,go]=opt_complex(f,g_cons,x0,xl,xu,TolX,TolFun,MaxIter)N=length(x0);M=size(g_cons);M=length(M(:,1));k1=0;k=N+1;%单纯形顶点个数gx=ones(M,1);whilemax(gx)0x0=xl+rand(N,1).*xu;gx=feval(g_cons,x0);end第33页，共72页。[x1,fx]=gen_complex(x0,k,f,g_cons);flag1=1;flag2=1;flag3=1;k1=0fxx1fprintf(此处暂停，请按下任意键继续\n)pausewhilek1MaxIterflag1=1;flag2=1;flag3=1;k1=k1+1[fx,I]=sort(fx);fori=1:kx2(:,i)=x1(:,I(i));endx1=x2;fmax1=fx(k);imax1=I(k);fmin=fx(1);imin=I(1);fmax2=fx(k-1);imax2=I(k-1);%计算形心xc=zeros(N,1);fori=1:kxc=xc+x1(:,i);endxc=xc-x1(:,imax1);xc=xc/(k-1);gxc=feval(g_cons,xc);alpha=1.31;%反射xr=xc+alpha*(xc-x1(:,imax1));gxr=feval(g_cons,xr)ifmax(gxr)0fxr=feval(f,xr);iffxrfmax1fprintf(反射成功\n)fmax1,fxrfmax1=fxr;fx(imax1)=fxr;x1(:,imax1)=xr;flag1=-1;else%反射失败flgg1=1;end第34页，共72页。else%反射失败flag1=1;endgama=0.7;ifflag1==-1fprintf(延伸\n)xe=xr+gama*(xr-xc);gxe=feval(g_cons,xe)ifmax(gxe)0fxe=feval(f,xe);iffxefmax1fprintf(延伸成功\n)fxe,fmax1fx(imax1)=fxe;fmax1=fxe