Lagrange 乘子法求最优解

clc
clear
syms x y z r1 r2 w
f=x^+y^+z^+w^;
g1=*x-y+z-w-;
g2=x+y-z+w-;
h=f-r1*g1 -r2*g2;

hx=diff(h,x);
hy=diff(h,y);
hz=diff(h,z);
hw=diff(h,w);
hr1=diff(h,r1);
hr2=diff(h,r2);

r=solve([hx==,hy==,hz==,hw==,hr1==,hr2==],[x,y,z,w,r1,r2]);

arr_x=double(r.x)
arr_y=double(r.y)
arr_z=double(r.z)
arr_w=double(r.w)

arr_fv=[];
for i= : length(arr_x)
xv=arr_x(i);
yv=arr_y(i);
zv=arr_z(i);
wv=arr_w(i);
fv=subs(f,[x,y,z,w],[xv,yv,zv,wv]);
arr_fv(i)= double(fv);
g1v=subs(g1,[x,y,z,w],[xv,yv,zv,wv])
g2v=subs(g2,[x,y,z,w],[xv,yv,zv,wv])
end
arr_fv
disp('after sort:')
sort(arr_fv)