基于MATLAB的乘子法

function [xmulambdaoutput]=multphr（funhfgfdfundhfdgfx0）
maxk=500;
sigma=2.0;
eta=2.0;theta=0.8;
k=0;ink=0;
epsilon=1e-5;
x=x0;
he=feval（hfx）;
gi=feval（gfx）;
n=length（x）;
l=length（he）;
m=length（gi）;
mu=0.1*ones（11）;
lambda=0.1*ones（m1）;
btak=10;
btaold=10;
while（btak>epsilon & k    [xvik]=bfgs（‘mpsi‘‘dmpsi‘x0funhfgfdfundhfdgfmulambdasigma）;
    ink=ink+ik;
    he=feval（hfx）;
    gi=feval（gfx）;
    btak=0.0;
    for（i=1:1）btak=btak+he（i）^2;
    end
    for（i=1:m）
        temp=min（gi（i）lambda（i）/sigma）;
        btak=btak+temp^2;
    end
    btak=sqrt（btak）;
    if btak>epsilon
        if（k>=2&btak>theta*btaold）
            sigma=eta*sigma;
        end
        for（i=1:1）mu（i）=mu（i）-sigma*he（i）;
        end
        for（i=1