Matlab快速入门教程.pptVIP

左图中当α=1时，L的曲面（网格线），其最小值（红色箭头所指）出现在x=(1,0)处。右图中，不断更改α的值，可以看到L的曲面开口逐渐收缩，最小值也逐渐下沉。这个现象从公式上不难理解，由于L=f(x)+α*g(x)，在g(x)0的区域，α越大，L比f小得越多（极值点下沉），而在在g(x)0的区域，α越大，L比f大得越多（开口收缩）。我们可以找到每一个L的极值点，将其连线，即得到右图中的红线，它实际就是θd(α)。【第二步】我们最大化θd(α)，此时可以将θd(α)看成是一个一元函数求极值(无条件)的问题，且α0。用拉格朗日乘数法即可求解。设D=θd(α)，则D对α偏导并令其为0有：求得α=1，另外一个解α=-3不满足条件。从而可知在x=(1,0)处可以取到极大值。clearall;%清除所有内存变量symsx1x2alphaalpha1betafx=x1.^2+x2.^2;%准备f(x)的表达式gx=(x1-2).^2+x2.^2-1;%准备g(x)的表达式Lx=fx+alpha*gx;%准备L(x)的表达式%%Lx求偏导数，得到三个方程dLx=[diff(Lx,x1);diff(Lx,x2);]%%solve函数解方程组[x1T,x2T]=solve(dLx,