惩罚函数的内点法.docVIP

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