07实验七 多元函数微分学.pptVIP

因此有两个极值可疑点。再输入: x=3^(1/2)/2-1/2; y=3^(1/2)/2-1/2; f1=eval(g) x=-3^(1/2)/2-1/2; y=-3^(1/2)/2-1/2; f2=eval(g) 得到输出: f1=0.2679 f2=3.7321 即得到两个可能是条件极值的函数值 但是否真的取到条件极值呢?可利用等高线作图来判断。 输入: [x,y]=meshgrid(-2:0.1:2, -2:0.1:2); z=x.^2+y.^2; contour(x,y,z,30) hold on ezplot(x^2+y^2+x+y-1) 输出如图7.4所示。 图7.4 从图7.4可以看到，在极值可疑点处，函数z=g(x,y)的等高线与曲线h(x,y)=0相切。函数z=g(x,y)的等高线是一系列同心圆，由里向外，函数值在增大， 在 的附近观察，可以得出z=g(x,y)取条件极大的结论。 在 的附近观察，可以得出z=g(x,y)取条件极小的结论。 MATLAB 高等数学实验 实验七 多元函数微分学 实验目的 掌握用MATLAB计算多元函数偏导数和全微分的方法，并掌握计算二元函数极值和条件极值的方法。理解和掌握曲面的切平面的作法。通过作图和观察，理解方向导数、梯度和等高线的概念。 7.1 学习MATLAB命令 7.1.1 求偏导数命令命令 diff既可以用于求一元函数的导数，也可用于求多元函数的偏导数。用于求偏导数时，可根据需要分别采用如下几种形式。 若求f(x,y,z)对x的偏导数，输入diff(f(x,y,z),x) 若求f(x,y,z)对y的偏导数，输入diff(f(x,y,z),y) 若求f(x,y,z)对x的二阶偏导数，输入diff(diff(f(x,y,z),x),x) 或者diff(f(x,y,z),x,2) 若求f(x,y,z)对x,y的混合偏导数，输入diff(diff(f(x,y,z),x),y) 其余类推。 7.1.2 在xoy平面上作二元函数 z=f(x,y)等高线的命令 contour的命令格式类似于mesh和surf这两个命令。 例如输入: [x,y]=meshgrid(-2:0.1:2,-2:0.1:2); z=x.^2-y.^2; contour(x,y,z,20) %参数20是等高线的数量 便作出了函数z=x^2-y^2的等高线(见图7.1)。 图7.1 7.1.3 解符号形式的代数方程组的命令 solve命令用于求解符号形式的代数方程组，其格式如下。 s=solve(eq1,eq2,…,eqN,var1,var2,…,varN) 它对方程组eq1，eq2，…，eqN中指定的N个变量var1，var2，…，varN求解。s返回解的结构，其内容通过阅读其域得到(输入s.var1，s.var2，…，等等)。当系统求不出解析解时，会自动求原点附近的一个近似解。 例如，输入: s=solve(‘x^2+x*y+y=3’,‘x^2-4*x+3=0’) %或者 s=solve(x^2+x*y+y-3,x^2-4*x+3) 输出为: s = x: [2x1 sym] y: [2x1 sym] 输入: s.x,s.y %结构的具体内容 输出为: ans= 1 3 ans= 1 -3/2 solve有另外一种输出形式。输入: [x,y]=solve(‘x^2+x*y+y=3’,‘x^2-4*x+3=0’) %或者 [x,y]=solve(x^2+x*y+y-3,x^2-4*x+3) 输出为: x= 1 3 y= 1 -3/2 7.2 实验内容 7.2.1 求多元函数的偏导数与全微分 【例1】设 ，求 输入: syms x y z=sin(x*y)+(cos(x*y))^2 diff(z,x) diff(z,y) diff(z,x,2) diff(diff(z,x),y) 便依次得到函数表达式及所求的四个偏导数结果: z=sin(x*y)+(cos(x*y))^2 ans=y*cos(x*y)-2*y*cos(x*y)*sin(x*y) ans=x*cos(x*y)-2*x*cos(x*y)*sin(x*y) ans=-2*y^2*cos(x*y)^2+2*y^2*sin(x*y)^2- y^2*sin(x*y) ans=-2*x*y