惩罚函数内点法c语言,分享：惩罚函数法（内点法、外点法）求解约束优化问题最优值...

方法一：外点牛顿法：

clc

m=zeros(1,50);a=zeros(1,50);b=zeros(1,50);f0=zeros(1,50);%a b为最优点坐标，f0为最优点函数值，f1 f2最优点梯度。

syms x1 x2 e; %e为罚因子。

m(1)=1;c=10;a(1)=0;b(1)=0; %c为递增系数。赋初值。

f=x1^2+x2^2+e*(1-x1)^2;f0(1)=1;

fx1=diff(f,'x1');fx2=diff(f,'x2');fx1x1=diff(fx1,'x1');fx1x2=diff(fx1,'x2');fx2x1=diff(fx2,'x1');fx2x2=diff(fx2,'x2');%求偏导、海森元素。

for k=1:100 %外点法e迭代循环.

x1=a(k);x2=b(k);e=m(k);

for n=1:100 %梯度法求最优值。

f1=subs(fx1); %求解梯度值和海森矩阵

f2=subs(fx2);

f11=subs(fx1x1);

f12=subs(fx1x2);

f21=subs(fx2x1);

f22=subs(fx2x2);

if(double