惩罚函数地内点法

《惩罚函数地内点法》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实用标准文案2013-2014(1)专业课程实践论文内点法精彩文档实用标准文案一、算法理论内点法总是从可行域的内点出发，并保持在可行域内进行搜索，因此这种方法适用于只有不等式约束条件的问题内点法据图计算步骤：1．给定初，允许误差，初始参数缩小系数；2．以为初始点，求解问题MinS.t.3．若则停，得近似解；否则令回2.二、算法框图求满足精彩文档实用标准文案从出发，求改变约束极值方法输出结果停三、算法程序clc精彩文档实用标准文案m=zeros(1,50);a=zeros(1,50);b=zeros(1,50);f0=

2、zeros(1,50);symsx1x2e;m(1)=1;c=0.2;a(1)=2;b(1)=-3;f=x1^2+x2^2-e*(1/(2*x1+x2-2)+1/(1-x1));f0(1)=15;fx1=diff(f,'x1');fx2=diff(f,'x2');fx1x1=diff(fx1,'x1');fx1x2=diff(fx1,'x2');fx2x1=diff(fx2,'x1');fx2x2=diff(fx2,'x2');fork=1:100x1=a(k);x2=b(k);e=m(k);forn=1:100f1

3、=subs(fx1);f2=subs(fx2);f11=subs(fx1x1);f12=subs(fx1x2);f21=subs(fx2x1);f22=subs(fx2x2);if(double(sqrt(f1^2+f2^2))<=0.002)a(k+1)=double(x1);b(k+1)=double(x2);f0(k+1)=double(subs(f));break;elseX=[x1x2]'-inv([f11f12;f21f22])*[f1f2]';x1=X(1,1);x2=X(2,1);endendif(d

4、ouble(sqrt((a(k+1)-a(k))^2+(b(k+1)-b(k))^2))<=0.001)&&(double(abs((f0(k+1)-f0(k))/f0(k)))<=0.001)a(k+1)b(k+1)kf0(k+1)break;elsem(k+1)=c*m(k);end精彩文档实用标准文案end四、算法实现精彩文档实用标准文案例1.利用内点法求解s.t解：改变算法中f=x1^2+x2^2-e*(1/(2*x1+x2-2)+1/(1-x1));回车完成结果复制粘贴代码，回车出现结果例2.利用内点法求解

5、精彩文档实用标准文案解：改变算法中f=x1^2-x2^2-e*(1/(x1+x2-3)+1/(-x1));回车完成结果例3.利用内点法求解精彩文档实用标准文案解：改变算法中f=2*x1^2-x2^2-e*(1/(x1+4*x2)+1/(x2-x1));回车完成结果例4.利用内点法求解精彩文档实用标准文案解：改变算法中f=3*x1^2+x2^2-e*(1/(2*x1+x2)+1/(x2-x1));回车完成结果精彩文档