MATLAB程序 惩罚函数法 内点法

%主调函数

%fun6_2_1.m

global r

r0=1;

c=0.1;

km=10;

x0=[2;4];

for k=1:km;

r=r0*c^(k-1);

[x,fv]=fminsearch(@fun,x0);

fprintf(1,'μú%3.0f′?μü′ú\n',k) fmin=fv

x0=x

end

%功能函数

function f=fun(x)

global r

f=x(1)^2+2*x(2)^2+r/(x(1)+x(2)-1);