terminate.m

来自「Lagrange乘子法 用于解约束最优化问题」· M 代码 · 共 17 行

function y=terminate(x)
global lumda sigma l m;
c=confun(x);
sum1=0;
for i=1:l
    sum1=sum1+c(i)^2;
end
sum2=0;
for i=(l+1):m 
    if c(i)>=lumda(i)/sigma;
        temp=c(i)^2;
    else
        temp=(lumda(i)/sigma)^2;
        sum2=sum2+temp;
    end
end
y=sqrt(sum1+sum2);