有约束优化问题

《有约束优化问题》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、任课教师签字：课程作业学年学期：2017——2018学年第二学期课程名称：优化理论作业名称：作业四学生姓名：学号：提交时间:一、问题重述形如下式的寻优问题称为无约束最优化问题，这类问题的最优解称为约束最优解。minf(x)s・t・gj(x)SO,i=l,.・.phy(x)=O,j=l,...q〃>0,g〉0约束优化问题的最优性条件是，在满足灯饰和不等式约束条件下，其目标函数值最小的点所必须满足的条件。约束优化设计问题求解方式有两种，间接法和直接法。直接法是在满足不等式约束的可行设计区域内直接搜索问题的最优解和最小值，常用的方法有随机方向法、复合形法。间接法是将优化问题转

2、化为一系列无约束优化问题来求解，常'用的方法有内惩罚函数法、外惩罚函数法以及混合惩罚函数法。本次作业以为例，介绍无约束最优化问题的寻优方法。二、算法原理复合形法是求解约束非线性寻优问题的一种重要的直接方法。复合形法的核心在于可行域内构造的不断逼近最优点的复合形。每次迭代,计算各顶点的目标函数值，找到目标函数值最大的顶点(称最坏点)，然后按相应的原则求出目标函数可行的下降点，以此代替最坏点，构成新的复合形。复合形每改变一次，各个顶点就向最优点移动一步，直至满足终止条件，找到最优点。复合形法的顶点数K通常取/2+1

3、者确定或者随机产生的，但一定要保证在可行域内。如果随机产生的初始点没有在可行域内，可以通过以下步骤将其调入可行域内。(1)计算在可行域内点的初始点集中心X(s);(2)将可行域外的点向X®靠拢，每次前进间距的一半，直至进入可行域内。复合行法的终止条件可以有以下几种形式，满足终止条件后，可将最后复合形的好点及其函数值作为最优解输出。(1)各顶点与好点函数值之差的均方根小于误差限；(2)各顶点与好点的函数之差平和小于误差限；(1)各顶点与好点函数值差的绝对值之和小于误差限。另外，需要注意的是复合形法只能求解仅含不等式约束的问题。三、算法流程复合形算法的流程图:图1复合形法算

4、法流程四、实验验证用Matlab自带的fmincon函数求得的约束最优解作为标准值，以此检验复合形算法的准确度。用复合形法求解下列问题的约束最优解min/(x)=60-10X]-4兀2-x}x2Mgi(x)=-X]<0tg2(x)=-x2+l<0

5、%%%%部分参数初始值area=L3;error=0.0000000000001;%误差限k二4;%复合形顶点数dimensions^;%搜索空间的维度a=[-10,-10];%初始点1b二[10,10];%初始点2c二zeros(1,2);x_in=[];%%%%%%%%%%%%%产生初始随机点，保证至少有一个点在可行域内whi1eisempty(x_in)x_in二[];x_out=[];fori=l:k%产生初始随机点c(l)=a(l)+rand*(b(l)-a(l));c(2)=a(2)+rand*(b(2)-a(2));ifInequal(c)>=0xin二[

6、xin;c];elsex_out二[x_out;c];endendend%%%%%%%%%%%%%%将不在Rj行域内的随机点调入町行域内while^isempty(x_out)xinmean=mean(xin,1);x_out(1,:)=x_in_mean+0・5*(x_out(1,:)-x_in_mean);ifInequal(xout(1,:))>0x_in二[x_in;x_out(1,:)];x_out(1,:)=[];endend都在可行域内的随机点x_zero二x_in;%%%%%%%%%%%%迭代复合形算法F=zeros(1,k);whi1e(1)fori=

7、l:kF(i)=Fmin(x_zero(i,:));end[F_L,x_L_i]=min(F);%好点[FH,xHi]=max(F);%坏点J=(F-F_L)*(F-F_D，;%各顶点与好点的函数值之差的平方和ifJ