%%%ex1208.m 拉各朗日乘子法求最優(yōu)化解
x=zeros(1,2);
%用syms表示出轉(zhuǎn)化后的無約束函數(shù)
syms x y lama
f=x+y+lama*(x^2+y^2-2);
%分別求函數(shù)關(guān)于x、y、lama的偏導(dǎo)
dx=diff(f,x);
dy=diff(f,y);
dlama=diff(f,lama);
%令偏導(dǎo)為零求解x、y
xx=solve(dx,x); %將x表示為lama函數(shù)
yy=solve(dy,y);  %將y表示為lama函數(shù)
ff=subs(dlama,{x,y},{xx,yy}); %代入dlama得關(guān)于lama的一元函數(shù)
lamao=solve(ff); %求解得lamao
xo=subs(xx,lama,lamao) %求得取極值處的xo
yo=subs(yy,lama,lamao)  %取極值處的yo
fo=subs(f,{x,y,lama},{xo,yo,lamao}) %極值點函數(shù)值