数学实验教程_实验7(多元函数微分学).docVIP

实验7 多元函数微分学 实验目的 1．加深理解偏导数的定义及几何意义 2．学会多元函数偏导数的计算方法 3．学会绘制曲面的切平面与法线 4．学会求多元函数的极值 实验准备 1．复习偏导数的定义以及几何意义 2．复习偏导数的计算方法，如链式法则，隐函数求导等 3．复习空间曲面的切平面以及法线的求法 4．复习求解无条件极值问题的方法 5．复习求解条件极值问题的拉格朗日乘子法 实验内容 1．利用几何图形演示偏导数的定义及几何意义 2．多元函数偏导数和全微分的计算 3．求曲面的切平面与法线方程并绘制图形 4．求多元函数的梯度，并绘制梯度场 5．求多元函数的极值 6．偏导数的应用 软件命令 表7-1 Matlab多元函数微分学命令 函数名称 调用格式 说 明 syms Syms 变量名1，变量名2,… 定义符号变量 sym sym(x,…) 定义符号变量 diff diff(f,x,n) 求函数f对x的n阶导数 plot3 plot3(x,y,z,可选项s) 绘制空间参数曲线 mesh mesh(x,y,z) 用两组相交的平行平面上的网状线方式来表示曲面 surf surf(x,y,z) 用网状线与补片填充色彩的方式来表示曲面 fminsearch fminsearch(fun,x0,options) 求多变量函数极小值 jacobian jacobian(f,v) 求符号函数的Jacobian矩阵 实验示例 【例7.1】偏导数的定义及几何意义 利用函数说明偏导数的几何意义。 【程序】：参见Exm07Demo01.m 【例7.2】偏导数的计算 （1）已知，求； （2）已知，且具有连续的二阶偏导数，求； （3）已知，求全微分。 【程序】： clear clc syms x y; % 第一小题 z=x^4+y^4-cos(2*x+3*y); zx=diff(z,x); zy=diff(z,y); zxx=diff(z,x,2); zyy=diff(z,y,2); zxy=diff(diff(z,x),y); zyx=diff(diff(z,y),x); zxyx=diff(zxy,x); % 第二小题 u=sym(f(y/x,x^2*y));% 定义复合抽象函数 ux=diff(u,x); uy=diff(u,y); uxx=diff(ux,x); uxy=diff(ux,y); uyx=diff(uy,x); uyy=diff(uy,y); % 第三小题 clear z; syms z dx dy dz; u=asin(z/(x^2+y^2)^(1/2)); ux=diff(u,x); uy=diff(u,y); uz=diff(u,z); du=[ux uy uz]*[dx dy dz]; 【例7.3】曲面的切平面与法线 求出旋转抛物面在处的切平面、法线方程并画出它们的图形。 【原理】： 一般曲面在点处的切平面方程为： 法线方程为： 【步骤】： 【Step1】：计算在指定点（1，1）处的偏导数值、切平面以及法线方程数据 clear;clc; syms x y; z=x^2+y^2; a=-2;b=2;c=-2;d=2; n=40; [X,Y]=meshgrid(linspace(a,b,n),linspace(-c,d,n)); % 曲面数据 Z=X.^2+Y.^2; % 切平面数据 x0=1;y0=1;z0=x0^2+y0^2; DFx=diff(z,x); DFy=diff(z,y); Fx0=subs(DFx,x,x0); Fy0=subs(DFy,y,y0); Z1=Fx0*(X-x0)+Fy0*(Y-y0)+z0; % 法线数据 t=-0.5:0.1:0.5; Xt=x0+Fx0*t; Yt=y0+Fy0*t; Zt=z0-t; 【Step2】：绘制旋转抛物面、切平面和法线 % 绘制图形 hold on grid on box on axis equal axis([Xmin Xmax Ymin Ymax Zmin Zmax]) surf(X,Y,Z); shading interp surf(X,Y,Z1,FaceColor,y); plot3(Xt,Yt,Zt,r,LineWidth,2); hold off 【输出】： 图7-3 切平面与法线  【例7.4】梯度、梯度线与等值线的实际应用 某景区计划堆积一座人工高地，其表面形状可近似地用曲面 表示（单位：米），在高地表面上的点（-24，-2.5，10.1）处设置人工喷泉口，让泉水自由地从高地上流下，直到水流到达水平高度为0.3米的地面某处。为了防止水流冲刷地表，施工时需要在水流路线上采取一定措施。为此，请帮助施工人员绘制出人工高