试验三多元函数微积分2.docVIP

附录Ⅰ 大学数学实验指导书 项目三 多元函数微积分 实验1 多元函数微分学（基础实验） 实验目的 掌握利用Mathematica计算多元函数偏导数和全微分的方法, 掌握计算二元 函数极值和条件极值的方法. 理解和掌握曲面的切平面的作法. 通过作图和观察, 理解二元 函数的性质、方向导数、梯度和等高线的概念. 求多元函数的偏导数与全微分 例1.1 (教材 例1.1) 设求 输入 Clear[z]; z=Sin[x*y]+Cos[x*y]^2; D[z,x] D[z,y] D[z,{x,2}] D[z,x,y] 则输出所求结果. 例1.2 设求和全微分dz. 输入 Clear[z];z=(1+x*y)^y; D[z,x] D[z,y] 则有输出 再输入 Dt[z] 则得到输出 例1.3 (教材 例1.2) 设其中a是常数, 求dz. 输入 Clear[z,a];z=(a+x*y)^y; wf=Dt[z,Constants-{a}]//Simplify 则输出结果: (a+xy)-1+y(y2Dt[x,Constants-{a}]+ Dt[y,Constants-{a}](xy+(a+xy)Log[a+xy])) 其中Dt[x,Constants-{a}]就是dx, Dt[y,Constants-{a}]就是dy. 可以用代换命令“/.”把它们 换掉. 输入 wf/.{Dt[x,Constants-{a}]-dx,Dt[y,Constants-{a}]-dy} 输出为 (a+xy)-1+y(dxy2+dy(xy+(a+xy)Log[a+xy])) 例1.4 (教材 例1.3) 设,求 输入 eq1=D[x==E^u+u*Sin[v],x,NonConstants-{u,v}] (*第一个方程两边对x求导数, 把u,v看成x,y的函数*) eq2=D[y==E^u-u*Cos[v],x,NonConstants-{u,v}] (*第二个方程两边对x求导数, 把u,v看成x,y的函数*) Solve[{eq1,eq2},{D[u,x,NonConstants-{u,v}], D[v,x,NonConstants-{u,v}]}]//Simplify (*解求导以后由eq1,eq2组成的方程组*) 则输出 其中D[u,x,NonConstants-{u,v}]表示u对x的偏导数, 而D[v,x,NonCosnstants-{u,v}]表示v 对x的偏导数. 类似地可求得u,v对y的偏导数. 微分学的几何应用 例1.5 求出曲面在点(1,1)处的切平面、法线方程, 并画出图形. 解(1) 画出曲面的图形. 曲面的参数方程为 输入命令 Clear[f]; f[x_,y_]=2x^2+y^2;p1=Plot3D[f[x,y], {x,-2,2},{y,-2,2}]; g1=ParametricPlot3D[{r*Sin[u]/Sqrt[2.],r*Cos[u],r^2}, {u,0,2*Pi},{r,0,2}] 则输出相应图形(图1.2). 图1.2 (2) 画出切平面的图形. 输入命令 a=D[f[x,y],x]/.{x-1,y-1}; b=D[f[x,y],y]/.{x-1,y-1}; p[x_,y_]=f[1,1]+a(x-1)+b(y-1); g2=Plot3D[p[x,y],{x,-2,2},{y,-2,2}]; 则输出切平面方程为及相应图形(图1.3). 图1.3 (3) 画出法线的图形. 输入命令 ly[x_]=1+b(x-1)/a;lz[x_]=f[1,1]-(x-1)/a; g3=ParametricPlot3D[{x,ly[x],lz[x]},{x,-2,2}]; Show[p1,g2,g3,AspectRatio-Automatic, ViewPoint-{-2.530,-1.025,2.000}]; 则输出相应图形(图1.4). 图1.4 例1.6 (教材 例1.4) 求曲面在点处的切平面方程, 并把曲面和它的切平面作在同一图形里. 输入 Clear[k,z]; k[x_,y_]=4/(x^2+y^2+1); (*定义函数k(x,y)*) kx=D[k[x,y],x]/.{x-1/4,y-1/2}; (*求函数k(x,y)对x的偏导数, 并代入在指定点的值*) ky=D[k[x,y],y]/.{x-1/4,y-1/2}; (*求函数k(x,y)对y的偏导数, 并代入在指定的值*) z=kx*(x-1/4)+ky*(y-1/2)+k[1/4,1/2]; (*定义在指定点的切平面函数*) 再输入 qm=Plot3D[k[x,y],{x,-2,2},{y,-2,2},PlotRange-{0,4}, Bo