486_极值问题.docVIP

实验二 极值问题 实验的目的 1、了解多元函数极值的求法； 2、了解多元函数条件极值的求法； 3、进一步掌握利用Matlab进行多元函数偏导数的求解与应用。 实验的基本理论与方法 1、多元函数极值的充分条件： 设在点的某邻域内具有二阶连续偏导数，又，，令，，，则在点处 1）AC-B20时具有极值，当A0时有极大值，当A0时有极小值； 2）AC-B20时没有极值； 3）AC-B2=0时不能确定是否有极值。 2、条件极值与拉格朗日(Lagrange)乘数法： 求函数在约束条件，下的条件极值的方法是，构造，求解方程组： 即可得到可能的极值点。 以上的方法可以推广到更多元的情形。 实验使用的函数与命令 函数求导指令diff； 2、方程求解指令sovle； 3、显示文本指令 disp； 4、创建二维等高线指令contour： 1）contour(Z)，参数Z为一个矩阵，表示相对于XY平面的高度，Z最小为2行2列的矩阵； 2）contour(Z, n)，根据矩阵绘制n组等高线； 3）contour(Z, v)，根据矢量v绘制指定等高线； 4）contour(X, Y, Z)或contour(X, Y, Z, n)或contour(X, Y, Z, v)，其中矢量X,Y,分别表示两个坐标范围。如果它们为矩阵，必须与矩阵Z大小相同，此时的Z为一般用函数surf创建的面。 实验指导 例1 解函数的极值。 解：首先求解偏导数。Matlab的M文件程序及结果： syms x y; f=x^3-y^3+3*x^2+3*y^2-9*x; diff(f,x) diff(f,y) ans =3*x^2+6*x-9 ans =-3*y^2+6*y 其次求解驻点坐标： [x,y]=solve(3*x^2+6*x-9=0,-3*y^2+6*y=0,x,y) 得到四个驻点为 P ( 1,0),Q(-3,0),R(1,2),S(-3,2) 第三步，求借二阶偏导数，并输出结果 A=diff(f,x,2) B=diff(diff(f,x),y) C=diff(f,y,2) A =6*x+6 B =0 C =-6*y+6 第四步，分别判别P, Q, R, S四点是否为极值，建立M文件，自动判断P, Q, R, S四点的极值情况： xx=[1 -3 1 -3]; %驻点横坐标 yy=[0 0 2 2]; %驻点纵坐标 for i=1:4 D=(6*xx(i)+6)*(-6*yy(i)+6) if D0 if (6*xx(i)+6)0 x=xx(i) y=yy(i) disp(为极大值点；) disp(极大值为) fmax=x^3-y^3+3*x^2+3*y^2-9*x elseif (6*xx(i)+6)0 x=xx(i) y=yy(i) disp(为极小值点；) disp(极小值为) fmin=x^3-y^3+3*x^2+3*y^2-9*x end end if D0 x=xx(i) y=yy(i) disp(该点不是极值点；) end if D==0 x=xx(i) y=yy(i) disp(无法确定！) end end 运行输出结果为 x =1 y =0 为极小值点； 极小值为 fmin = -5 x = -3 y =0 该点不是极值点； x =1 y =2 该点不是极值点； x = -3 y =2 为极大值点； 极大值为 fmax =31 下面绘出函数图形观察极值点和鞍点的情形，在函数曲面图3.1左中，观察不到细节。而右图的等高线图中有两个极值点(1, 0 )，(-3, 2 )，又因为极值点有等高线环绕，而(-3, 0 )，(1, 2)周围没有等高线，故不是极值点，是鞍点。 x=-5:0.1:5; y=-1:0.1:3; [X,Y]=meshgrid(x,y); Z=X.^3-Y.^3+3*X.^2+3*Y.^2-9*X; subplot(2,1,1); mesh(X,Y,Z); title(函数曲面图) subplot(2,1,2); contour(X,Y,Z,200); xlabel(x); ylabel(y); title(等高线图) 例2 抛物面被平面截成一椭圆，求原点到该椭圆的最长与最短距离。 解：先绘出图形 x