用拉格朗日乘子法求解最优化程序

用拉格朗日乘子法求解以下约束最优化问题：

解：先构造函数

 分别对L中的x和u求偏导数，则有

解上述方程组，得

拉格朗日乘子法的M文件：

syms u

f=x1.^2+x2.^2-x1.*x2-10.*x1-4.*x2+60;

h=x1+x2-6;

L=f+u.*h;

dldx1=diff(L,'x1');

dldx2=diff(L,'x2');

dldu=diff(L,'u');

s1='2*x1-x2-10+u';

s2='2*x2-4+u';

s3='x1+x2-6';

[u,x1,x2]=solve(s1,s2,s3);

f=x1.^2+x2.^2-x1.*x2-10.*x1-4.*x2+60;

h=x1+x2-6;

L=f+u.*h;

x1

x2

u

f

h

L

运行结果为：

x1 =24/5

x2 = 6/5

u =8/5

f = 8/25

h = 0

L =8/25