大学数学实验指导书.docVIP

附录Ⅰ 大学数学实验指导书 项目三 多元函数微积分 实验1 多元函数微积分（基础实验） 实验目的 掌握利用Mathematica计算多元函数偏导数和全微分的方法, 掌握计算二元 函数极值和条件极值的方法. 通过作图和观察, 理解二元函数的性质. 掌握利用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.4) 求的极值. 输入 Clear[f]; f[x_,y_]=x^3-y^3+3x^2+3y^2-9x; fx=D[f[x,y],x] fy=D[f[x,y],y] critpts=Solve[{fx==0,fy==0}] 则分别输出所求偏导数和驻点: {{x--3,y-0},{x--3,y-2},{x-1,y-0},{x-1,y-2}} 再输入求二阶偏导数和定义判别式的命令 fxx=D[f[x,y],{x,2}]; fyy=D[f[x,y],{y,2}]; fxy=D[f[x,y],x,y]; disc=fxx*fyy-fxy^2 输出为判别式函数的形式: (6+6x)(6-6y) 再输入 data={x,y,fxx,disc,f[x,y]}/.critpts; TableForm[data,TableHeadings-{None,{ x , y , fxx , disc , f }}] 最后我们得到了四个驻点处的判别式与的值并以表格形式列出. X y fxx disc f -3 0 -12 -72 27 -3 2 -12 72 31 1 0 12 72 -5 1 2 12 -72 -1 易见,当时判别式disc=72, 函数有极大值31; 当时判别式disc=72, 函数有极小值-5; 当和时, 判别式disc=-72, 函数在这些点没有极值. 注:在项目一的实验4中,我们曾用命令FindMinimum来求一元函数的极值, 实际上,也可 以用它求多元函数的极值, 不过输入的初值要在极值点的附近. 对本例,可以输入以下命令 FindMinimum[f[x,y],{x,-1},{y,1}] 则输出 {-5.,{x-1.,y--2.36603×10-8}} 从中看到在的附近函数有极小值-5, 但y的精度不够好. *例1.6 求函数在条件下的极值. 输入 Clear[f,g,la]; f[x_,y_]=x^2+y^2; g[x_,y_]=x^2+y^2+x+y-1; la[x_,y_,r_]=f[x,y