实验2 多元函数积分学基础实验.docVIP

项目三 多元函数微积分 实验2 多元函数积分学（基础实验） 实验目的 掌握用Mathematica计算二重积分与三重积分的方法; 深入理解曲线积分、曲面积分的 概念和计算方法. 提高应用重积分和曲线、曲面积分解决各种问题的能力. 基本命令 1. 计算重积分的命令lntegrate和NIntegrate 例如,计算, 输入 Integrate[x*y^2,{x,0,1},{y,0,x}] 则输出 又如,计算的近似值, 输入 NIntegrate[Sin[x*y^2],{x,0,1},{y,0,1}] 则输出 0.160839 注: Integrate命令先对后边的变量积分. 计算三重积分时,命令Integrate的使用格式与计算二重积分时类似. 由此可见, 利用 Mathematica计算重积分表示成,的函数. 例如,在直角坐标系中方 程是一旋转抛物面, 在柱坐标系中它的方程为. 因此,输入 CylindricalPlot3D[r^2,{r,0,2},{t,0,2Pi}] 则在柱坐标系中作出了该旋转抛物面的图形. 3. 球面坐标系中作三维图形命令SphericalPlot3D 使用命令SphericalPlot3D, 首先要调出作图软件包. 输入 Graphics`ParametricPlot3D` 执行成功后便可继续下面的工作. 命令SphericalPlot3D的基本格式为 SphericalPlot3D[r[, { 其中r[是曲面的球面坐标方程, 使用时一定要把球面坐标中的表示成、的函数. 例如,在球面坐标系中作出球面 输入 Sphericalplot3D[2,{u,0,pi},|v,0,2,pi|,plotpoints-40] 则在球面坐标系中作出了该球面的图形. 4. 向量的内积 用“.”表示两个向量的内积. 例如,输入 vecl={al,bl,cl} vec2={a2,b2,c2} 则定义了两个三维向量, 再输入 vec1. vec2 则得到它们的内积 a1a2+b1b2+c1c2 实验举例 计算重积分 例2.1 (教材 例2.1) 计算 其中为由 所围成的有界区域.先作出区域的草图, 易直接确定积分限,且应先对积分, 因此, 输入 Integrate[x*y^2,{y,1,2},{x,2-y,Sqrt[y]}] 则输出所求二重积分的计算结果 例2.2 (教材 例2.2) 计算 其中为 如果用直角坐标计算, 输入 Clear[f,r]; f[x,y]=Exp[-(x^2+y^2)]; Integrate[f[x,y],{x,-1,1},{y,-Sqrt[1-x^2],Sqrt[1-x^2]}] 则输出为 其中Erf是误差函数. 显然积分遇到了困难. 如果改用极坐标来计算, 也可用手工确定积分限. 输入 Integrate[(f[x,y]/.{x-r*Cos[t],y-r*Sin[t]})*r,{t,0,2 Pi},{r,0,1}] 则输出所求二重积分的计算结果 如果输入 NIntegrate[(f[x,y]/.{x-r*Cos[t],y-r*Sin[t]})*r,{t,0,2 Pi},{r,0,1}] 则输出积分的近似值 1.98587 例2.3 (教材 例2.3) 计算, 其中由曲面与围成. 先作出区域的图形. 输入 g1=ParametricPlot3D[{Sqrt[2]*Sin[fi]*Cos[th], Sqrt[2]*Sin[fi]*Sin[th], Sqrt[2]*Cos[fi]},{fi,0,Pi/4},{th,0,2Pi}] g2=ParametricPlot3D[{z*Cos[t],z*Sin[t],z},{z,0,1},{t,0,2Pi}] Show[g1,g2,ViewPoint-{1.3,-2.4,1.0}] 则分别输出三个图形（图2.1(a), (b), (c)）. (a) (b) 图2.1 考察上述图形, 可用手工确定积分限. 如果用直角坐标计算, 输入 g[x_,y_,z_]=x^2+y^2+z; Integrate[g[x,y,z],{x,-1,1},{y