实验七重积分及其MATLAB计算

实验七 重积分及其MATLAB计算 一、实验目的 用MATLAB实现二重积分的计算，学习高等数学中有关三重积分的定义和计算方法，掌握利用MATLAB软件进行三重积分的计算。 二、相关知识 在高等数学中我们还知道，二重积分是化为二次定积分运算来实现的，因此我们利用MATLAB的函数int( )，还可以完成二重积分的计算。 例1、计算二重积分 ，其中D为由曲线 和 所围成的区域。 解：区域可用不等式表示为： 其所围的区域如图7-1所示。 所以，计算该积分的 MATLAB程序为： 二、相关知识 clear syms x y f=x*x+y; int(int(f,y,x*x,sqrt(x)),x,0,1) 计算结果为：33/140 二、相关知识 例2 计算 ，其中D是由圆周 所围成的闭区域。 解：积分区域如图2-2所示这个 积分需要化为极坐标来进行积分 化为极坐标后，积分为： 所以，可以用如下的程序 来完成本题的计算： 二、相关知识 clear syms r s int(int(r*r,r,0,cos(s)),s,-pi/2,pi/2) 计算结果为：4/9。 二、相关知识 定积分及二重积分作为和式极限的概念，可以很自然地推广到三重积分。 设 是空间有界闭区域 上的有界函数，将 任意分成 n 个小闭区域 其中 表示第i个小闭区域，也表示它的体积。在每个 上任取一点 ，作乘积 ，并作和 。如果当各小闭区域直径中的最大值 趋于零时这和的极限总存在，则称此极限为函数 在闭区域 上的三重积分。记作 ，即 其中 叫做体积元素。 二、相关知识 计算三重积分的基本方法是将三重积分化为三次积分来计算。下面按利用不同的坐标来分别讨论将三重积分化为三次积分的方法。 1. 利用直角坐标计算三重积分 积分区域为： 则三重积分的计算公式为： 先对z、次对y、最后对x的三次积分 二、相关知识 如果平行于x轴或y轴且穿过闭区域Ω内部的直线与Ω的边界曲面S相交不多于两点，也可把闭区域Ω投影到yOz面上或xOz面上，这样便可把三重积分化为按其他顺序的三次积分。如果平行于坐标轴且穿过闭区域Ω内部的直线与边界面S的交点多于两个，也可像处理二重积分那样，把Ω分成若干部分，使Ω上的三重积分化为各部分闭区域上的三重积分的和。 二、相关知识 例3 计算三重积分 ，其中Ω为三个坐标面及平面 所围成的闭区域。 解：作闭区域Ω如图7.2所示。 将Ω投影到面xOy上，得投影区域 为三角形闭区域OAB。直线OA、 OB及AB的方程依次为 y=0、x=0 及 x+2y=1 ，所以 二、相关知识 在 内任取一点 ，过此点作平行于z轴的直线，该直线通过平面z=0穿入Ω内，然后通过平面z=1-x-2y穿出Ω外。 于是，由公式（7.2）得 此时，编写MATLAB程序如下： syms x y z I=int(int(int(x,z,0,1-x-2*y),y,0,(1-x)/2),x,0,1) 得到结果为 I=1/48 二、相关知识 例4 计算三重积分 ，其中Ω是由椭球面 所围成的空间闭区域。 解：空间闭区域Ω可表示为 其投影区域 为： ， 因此可表示为 故三重积分 其中： ， 二、相关知识 此时，编写MATLAB程序如下： syms a b c x y z d=c*sqrt(1-x*x/(a*a)-y*y/(b*b)) e=b*sqrt(1-x*x/(a*a)) A=int(int(int(z*z,z,-d,d),y,-e,e),x,-a,a) 运行结果为： A= 4/15*c^3*pi/(1/b^2)^(1/2)*a 二、相关知识 2. 利用柱面坐标计算三重积分 空间曲面的方程可以用多种形成的坐标来表达，其中一种表达方式称为柱面坐标，某些情况下，利用柱面坐标可使计算简化。 设 为空间内一点，并设点M在xOy面上的投影P的极坐标为 ，则这样的三个数 就叫做点