ch9 二重积分-(4月13日-15日)
第一节 二重积分的概念与性质
三、计算题
1．利用二重积分性质估计积分()2
22d d D
I x
y x y =++??的值，其中1x y +≤.
2. 根据二重积分的性质，比较
2()d D
x y σ+??与3
()d D
x y σ+??的大小，其中D 由圆周22(2)(1)2x y -+-=围成.
第二节 二重积分的计算（一）
一、填空题
1. 设积分区域D 由x 轴，直线1x =及曲线2
y x =所围成，其X 型区域可用不等式表示为 ，Y 型区域可用不等式表示为 . 2. 设D 为直线y x =及抛物线24y x =所围成的闭区域，化二重积分(,)d d D
f x y x y ??为直角
坐标下的二次积分为 . 3．设220
d (,)d x
x
I x f x y y =
??
，改变积分次序，则I = .
二、选择题
1．设D 是以(0,0)O ，(1,0)A ，(1,2)B ，(0,1)C 为顶点的梯形所围成的闭区域，则
(,)d D
f x y σ??化成二次积分是 （ ）
A.110
1
d (,)d x
x f x y y +??
； B.120
1
d (,)d x
x f x y y +??
；
C.
1
121
1
1
d (,)d d (,)d y y f x y x y f x y x -+?
???
； D.11
2
1
1
d (,)d d (,)d y y f x y x y f x y x -+????
.
2. 设
10
()d f x x a =?
，{}(,)01,01D x y x y =≤≤≤≤，则()()d d D
f x f y x y =??（ ）
A. a ；
B. 2
a ； C.
2
12
a ； D. 0. 三、计算题
1. 计算d d D
xy x y ??
，其中D 是由直线1y x =-与抛物线2
26y x =+所围成的闭区域.
2. 计算??+D
y x d e σ23．其中2||,2||:≤≤y x D ．
第二节 二重积分的计算法(二)
一、填空题
1. 设积分区域D 由x 轴,y =
D 用极坐标型的不等式表示为 .
2. 设D 是{}
22
(,)14x y x y ≤+≤,则积分
(,)d d D
f x y x y ??表示为极坐标形式的二次积分
为 ________.
3. 化
20
d d x
x f
y ?
?
为极坐标形式的二次积分为 .
二、单项选择题
1. 设区域{}
222
(,),0,0D x y x y a a y =+≤>≥，则
()2
2d d =D
x
y x y +?? .
A ．
30
d d a π
θρρ?
? B ．20
d d a
πθρρ??
C.
3
20
2
d d a
π
πθρρ-??
D.
220
2
d d a
π
πθρρ-??
2. 二次积分cos 20
d (cos ,sin )d f πθ
θρθρθρρ??
可以写成 .
A.
1
d (,)d y f x y x ? B.
1
d (,)d y f x y x ?
C.
11
d (,)d x f x y y ??
D.
1
d (,)d x f x y y ?
3. 设积分区域D ：2
2
4x y +≤，则22
e
d d =x y D
x y +?? .
A.
4(e 1)2
π
- B. 42(e 1)π- C. 4(e 1)π- D. 4e π
三、计算题
1. 计算d D
x y ,其中{}22(,)0,2D x y y x x y x =≤≤+≤.
2．??++D
d y x
σ)1ln(22
．4:22≤+y x D ，0≥x ，0≥y ．.
3．交换下列累次积分的积分顺序： ⑴?
?
--a
a
x a dy y x f dx 2
20
),(． ⑵??
??
-+31
30
10
20
),(),(y
y
dx y x f dy dx y x f dy ．
4. 选择适当的坐标计算??+D
d y x σ)(22．其中1||||:≤+y x D ．
5．计算2
2|4|d d D
x
y x y +-??，其中D 是区域229x y +≤.
二重积分的概念
第一节 二重积分的概念与性质 一、内容要点 1、引例 例1曲顶柱体的体积 例2平面薄片的质量 通过两个实际意义不同的例子，引出所求量可归结为同一形式的和式的极限，进而一般地抽象出二重积分的定义。 2、二重积分的概念：注意讲清楚定义中两个“任意性”及和式极限中各符号的意义。 3、二重积分的性质1-6，注意将其与定积分性质加以比较。 例3关于估值定理的应用 例4关于中值定理的应用 4、二重积分的几何意义——曲顶柱体的体积。 二、教学要求和注意点 理解二重积分，了解重积分的性质，了解二重积分的中值定理。 第二节 二重积分的计算法 一、内容要点 利用直角坐标计算二重积分 1、从几何入手，利用计算“平行截面面积为已知的立体的体积”方法，将二重分化为二次积分： ①若D 为X —型区域：{}b x a x y x y x ≤≤≤≤),()(),(21?? 则 ????=D x x b a dy y x f dx d y x f )()(21),(),(??σ ②若D 为Y —型区域：{}d y c y x y y x ≤≤≤≤),()(),(21?? 则 ????=D y y d c dx y x f dy d y x f )()(21),(),(??σ ③若D 既非X —型，又非Y —型区域，则将D 划分为若干子区域，使每一个子区域为X —型或Y —型。 2、介绍“对称性”在二重积分计算中的应用。 例1化二重积分为二次积分并求值，通过例子说明确定积分限的方法。 例2更换积分次序并计算，通过该例说明选择积分次序的重要性。
例3关于利用对称性计算二重积分的例子。 例4被积函数为绝对值函数、符号函数，取最大值或最小值等函数的例子。 利用极坐标计算二重积分 1、介绍极坐标下二重积分的换元公式。 2、何时选用极坐标进行计算，一般说来，当积分域D 的边界曲线用极坐标方程表示比较简单或被积函数用极坐标表示比较简单，可考虑用积坐标计算。 3、确定积分上下限的办法。 例1将直角坐标系下的二次积分化为极坐标系下的二次积分 例2利用二重积分计算概率积分 dx e x 2 0-+∞? 例3将极坐标系下的二次积分化为直角坐标系下的二次积分 例4利用极坐标计算二重积分 二、教学要求和注意点 1、掌握二重积分（直角坐标、极坐标）的计算方法 2、将重积分化为累次积分计算时，积分限的确定要保持每个单积分的下限小于上限，因此在交换二次积分次序时应注意符号问题。 3、在二重积分的计算时应尽量利用区域和被积函数的对称性以简化计算。 第四节 三重积分 一、内容要点 1、三重积分的概念，存在性及性质 2、三重积分在直角坐标系下的计算 ①先单积分后二重积分 ②先二重积分后单积分 3、更换积分次序 例1将三重积分化为三次积分 例2更换积分次序 例3先二重积分后单积分 4、柱面坐标系下三重积分的计算。 5、何时选用柱面坐标——当Ω是柱形，锥形或旋转体且在坐标面上的投影是圆域或其部分，或者被积函数含有式子)(22y x +?等时，常用柱面坐标计算。 6、球面坐标系下三重积分的计算。 7、何时选用球面坐标——当Ω是球体或其部分，或被积函数含有式子)(222z y x ++?
用递推公式计算定积分(matlab版)
用递推公式计算定积分 实验目的： 1．充分理解不稳定的计算方法会造成误差的积累，在计算过程中会导致误差的迅速增加，从而使结果产生较大的误差。 2．在选择数值计算公式来进行近似计算时，应学会选用那些在计算过程中不会导致误差迅速增长的计算公式。 3．理解不稳定的计算公式造成误差积累的来源及具体过程； 4．掌握简单的matlab语言进行数值计算的方法。 实验题目： 对n=0,1,2,…,20，计算定积分： 实验原理： 由于y(n)= = – 在计算时有两种迭代方法，如下： 方法一： y(n)=– 5*y(n-1),n=1,2,3, (20) 取y(0)= = ln6-ln5 ≈ 0.182322 方法二： 利用递推公式：y(n-1)=-*y(n),n=20,19, (1) 而且，由 = * ≤≤* =
可取：y(20)≈*()≈0.008730. 实验容： 对算法一,程序代码如下： function [y,n]=funa() syms k n t; t=0.182322; n=0; y=zeros(1,20); y(1)=t; for k=2:20 y(k)=1/k-5*y(k-1); n=n+1; end y(1:6) y(7:11) 对算法二，程序代码如下： %计算定积分； %n--表示迭代次数； %y用来存储结果； function [y,n]=f(); syms k y_20;
y=zeros(21,1); n=1; y_20=(1/105+1/126)/2; y(21)=y_20; for k=21:-1:2 y(k-1)=1/(5*(k-1))-y(k)/5; n=n+1; end 实验结果： 由于计算过程中，前11个数字太小，后9个数字比较大，造成前面几个数字只显示0.0000的现象，所以先输出前6个，再输出7—11个，这样就能全部显示出来了。 算法一结果： [y,n]=funa %先显示一y(1)—y(6) ans = 0.1823 -0.4116 2.3914 -11.7069 58.7346
数值积分算法与MATLAB实现陈悦5133201讲解
东北大学秦皇岛分校 数值计算课程设计报告 数值积分算法及MATLAB实现 学院数学与统计学院 专业信息与计算科学 学号5133201 姓名陈悦 指导教师姜玉山张建波 成绩 教师评语： 指导教师签字： 2015年07月14日
1 绪论 数值分析是计算数学的一个主要部分，计算数学是数学科学的一个分支，它研究用计算机求解各种数学问题的数值检索方其理论与软件的实现.而数值分析主要研究数值计算. 现科学技术的发展与进步提出了越来越多的复杂的数值计算问题，这些问题的圆满解决已远人工手算所能胜任，必须依靠电子计算机快速准确的数据处理能力.这种用计算机处理数值问题的方法，成为科学计算.今天，科学计算的应用范围非常广泛，天气预报、工程设计、流体计算、经济规划和预测以及国防尖端的一些科研项目，如核武器的研制、导弹和火箭的发射等，始终是科学计算最为活跃的领域. 1.1 数值积分介绍 数值积分是数值分析的重要环节，实际问题当中常常需要计算积分，有些数值方法，如微分方程和积分方程的求解，也都和积分计算相联系. 求某函数的定积分时，在多数情况下，被积函数的原函数很难用初等函数表达出来，因此能够借助微积分学的牛顿-莱布尼兹公式计算定积分的机会是不多的.另外，许多实际问题中的被积函数往往是列表函数或其他形式的非连续函数，对这类函数的定积分，也不能用不定积分方法求解.由于以上原因，数值积分的理论与方法一直是计算数学研究的基本课题.对微积分学做出杰出贡献的数学大师，如I.牛顿、L.欧拉、C.F.高斯、拉格朗日等人都在数值积分这个领域作出了各自的贡献，并奠定了这个分支的理论基础. 构造数值积分公式最通常的方法是用积分区间上的n 次插值多项式代替被积函数，由此导出的求积公式称为插值型求积公式.特别在节点分布等距的情形称为牛顿-科特斯公式，例如梯形公式(Trapezoidal Approximations)与抛物线公式(Approximations Using Parabolas)就是最基本的近似公式.但它们的精度较差.龙贝格算法是在区间逐次分半过程中，对梯形公式的近似值进行加权平均获得准确程度较高的积分近似值的一种方法，它具有公式简练、计算结果准确、使用方便、稳定性好等优点，因此在等距情形宜采用龙贝格求积公式(Rhomberg Integration).当用不等距节点进行计算时，常用高斯型求积公式计算，它在节点数目相同情况下，准确程度较高，稳定性好，而且还可以计算无穷积分.数值积分还是微分方程数值解法的重要依据.许多重要公式都可以用数值积分方程导出.现探讨数值积分算法以及运用MATLAB软件的具体实现
matlab求定积分之实例说明
一、符号积分 符号积分由函数int来实现。该函数的一般调用格式为： int(s)：没有指定积分变量和积分阶数时，系统按findsym函数指示的默认变量对被积函数或符号表达式s求不定积分； int(s,v)：以v为自变量，对被积函数或符号表达式s求不定积分； int(s,v,a,b)：求定积分运算。a,b分别表示定积分的下限和上限。该函数求被积函数在区间[a,b]上的定积分。a和b可以是两个具体的数，也可以是一个符号表达式，还可以是无穷(inf)。当函数f关于变量x在闭区间[a,b]上可积时，函数返回一个定积分结果。当a,b中有一个是inf时，函数返回一个广义积分。当a,b中有一个符号表达式时，函数返回一个符号函数。 例： 求函数x^2+y^2+z^2的三重积分。内积分上下限都是函数，对z积分下限是sqrt(x*y)，积分上限是x^2*y；对y积分下限是sqrt(x)，积分上限是x^2；对x的积分下限1，上限是2，求解如下： >>syms x y z %定义符号变量 >>F2=int(int(int(x^2+y^2+z^2,z,sqrt(x*y),x^2*y),y,sqrt(x),x^2),x,1,2) %注意定积分的书写格式 F2 = 1610027357/6563700-6072064/348075*2^(1/2)+14912/4641*2^(1/4)+64/225*2 ^(3/4) %给出有理数解 >>VF2=vpa(F2) %给出默认精度的数值解 VF2 = 224.92153573331143159790710032805 二、数值积分 1.数值积分基本原理 求解定积分的数值方法多种多样，如简单的梯形法、辛普生(Simpson)?法、牛顿－柯特斯(Newton-Cotes)法等都是经常采用的方法。它们的基本思想都是将整个积分区间[a,b]分成n个子区间[xi,xi+1]，i=1,2,…,n，其中x1=a，xn+1=b。这样求定积分问题就分解为求和问题。 2.数值积分的实现方法 基于变步长辛普生法，MATLAB给出了quad函数来求定积分。该函数的调用格式为： [I,n]=quad('fname',a,b,tol,trace) 基于变步长、牛顿－柯特斯(Newton-Cotes)法，MATLAB给出了quadl函数来求定积分。该函数的调用格式为： [I,n]=quadl('fname',a,b,tol,trace) 其中fname是被积函数名。a和b分别是定积分的下限和上限。tol用来控制积分精度，缺省时取tol=0.001。trace控制是否展现积分过程，若取非0则展现积分过程，取0则不展现，缺省时取trace=0。返回参数I即定积分值，n为被积函数的调用次数。 例： 求函数'exp(-x*x)的定积分，积分下限为0，积分上限为1。 >>fun=inline('exp(-x.*x)','x'); %用内联函数定义被积函数fname
二重积分学习总结
高等数学论文 《二重积分学习总结》 ：徐琛豪 班级：安全工程02班 学号：1201050221 完成时间：2013年6月2日
二重积分 【本章学习目标】 ⒈理解二重积分的概念与性质，了解二重积分的几何意义以及二重积分与定积分之间的联系，会用性质比较二重积分的大小，估计二重积分的取值围。 ⒉领会将二重积分化为二次积分时如何确定积分次序和积分限，如何改换二次积分的积分次序，并且如何根据被积函数和积分区域的特征选择坐标系。熟练掌握直角坐标系和极坐标系下重积分的计算方法。 ⒊掌握曲顶柱体体积的求法，会求由曲面围成的空间区域的体积。 1 二重积分的概念与性质 1．二重积分定义 为了更好地理解二重积分的定义，必须首先引入二重积分的两个“原型”，一个是几何的“原型”－曲顶柱体的体积如何计算，另一个是物理的“原型”—平面薄片的质量如何求。从这两个“原型”出发，对所抽象出来的二重积分的定义就易于理解了。 在二重积分的定义中，必须要特别注意其中的两个“任意”，一是将区域D 成n 个小区域12,,,n σσσ???的分法要任意，二是在每个 小区域i σ?上的点(,)i i i ξησ∈?的取法也要任意。有了这两个“任意”， 如果所对应的积分和当各小区域的直径中的最大值0λ→时总有同一个极限，才能称二元函数(,)f x y 在区域D 上的二重积分存在。 2．明确二重积分的几何意义。 (1) 若在D 上(,)f x y ≥0，则(,)d D f x y σ??表示以区域D 为底，以 (,)f x y 为曲顶的曲顶柱体的体积。特别地，当(,)f x y ＝1时，(,)d D f x y σ ??表示平面区域D 的面积。
用MATLAB算多元函数积分
用MATLAB 计算多元函数的积分 三重积分的计算最终是化成累次积分来完成的，因此只要能正确的得出各累次积分的积分限，便可在MA TLAB 中通过多次使用int 命令来求得计算结果。但三重积分的积分域Ω是一个三维空间区域，当其形状较复杂时，要确定各累次积分的积分限会遇到一定困难，此时，可以借助MATLAB 的三维绘图命令，先在屏幕上绘出Ω的三维立体图，然后执行命令 rotate3d on ↙ 便可拖动鼠标使Ω的图形在屏幕上作任意的三维旋转，并且可用下述命令将Ω的图形向三个坐标平面进行投影： view(0,0),向XOZ 平面投影； view(90,0),向YOZ 平面投影； view(0,90),向XOY 平面投影. 综合运用上述方法，一般应能正确得出各累次积分的积分限。 例11.6.1计算zdv Ω ???，其中Ω是由圆锥曲面222z x y =+与平面z=1围成的闭区域 解 首先用MA TLAB 来绘制Ω的三维图形，画圆锥曲面的命令可以是： syms x y z ↙ z=sqrt(x^2+y^2); ↙ ezsurf(z,[-1.5,1.5]) ↙ 画第二个曲面之前，为保持先画的图形不会被清除，需要执行命令 hold on ↙ 然后用下述命令就可以将平面z=1与圆锥面的图形画在一个图形窗口内： [x1,y1]=meshgrid(-1.5:1/4:1.5); ↙ z1=ones(size(x1)); ↙ surf(x1,y1,z1) ↙ 于是得到Ω的三维图形如图:
由该图很容易将原三重积分化成累次积分： 111zdv dy -Ω=???? 于是可用下述命令求解此三重积分： clear all ↙ syms x y z ↙ f=z; ↙ f1=int(f,z.,sqrt(x^2+ y^2),1); ↙ f2=int(f1,x,-sqrt(1- y^2), sqrt(1- y^2)); ↙ int(f2,y,-1,1) ↙ ans= 1/4*pi 计算结果为4 π 对于第一类曲线积分和第一类曲面积分，其计算都归结为求解特定形式的定积分和二重积分，因此可完全类似的使用int 命令进行计算，并可用diff 命令求解中间所需的各偏导数。 例11.6.2用MATLAB 求解教材例11.3.1 解 求解过程如下 syms a b t ↙ x=a*cos(t); ↙ y=a*sin(t); ↙ z=b*t; ↙ f=x^2 +y^2+z^2; ↙ xt=diff(x,t); ↙ yt=diff(y,t); ↙ zt=diff(z,t); ↙ int(f*sqrt(xt^2 +yt^2+zt^2),t,0,2*pi) ↙ ans= 2/3*( a^2 +b^2)^1/2*a^2*pi+8/3*( a^2 +b^2)^1/2*b^2*pi^3 对此结果可用factor 命令进行合并化简： factor （ans ） ans= 2/3*( a^2 +b^2)^1/2*pi*(3* a^2 +4*b^2*pi^2) 例11.6.3用MATLAB 求解教材例11.4.1 解 求解过程如下 syms x y z1 z2↙ f= x^2 +y^2; ↙ z1=sqrt(x^2 +y^2); ↙ z2=1; ↙ z1x=diff(z1,x); ↙ z1y=diff(z1,y); ↙ z2x=diff(z2,x); ↙ z2y=diff(z2,y); ↙
二重积分的概念及性质
二重积分的概念及性质 前面我们已经知道了，定积分与曲边梯形的面积有关。下面我们通过曲顶柱体的体积来引出二重积分的概念，在此我们不作详述，请大家参考有关书籍。 二重积分的定义 设z=f(x,y)为有界闭区域(σ)上的有界函数： (1)把区域(σ)任意划分成n个子域(△σk)(k=1,2,3,…,n）,其面积记作△σk(k=1,2,3,…,n）； (2)在每一个子域(△σk)上任取一点，作乘积； (3)把所有这些乘积相加,即作出和数 (4)记子域的最大直径d.如果不论子域怎样划分以及怎样选取,上述和数当n→＋∞且d→0时的极限存在,那末称此极限为函数f(x,y)在区域(σ)上的二重积分.记作: 即：= 其中x与y称为积分变量，函数f(x,y)称为被积函数,f(x,y)dσ称为被积表达式,(σ)称为积分区域. 关于二重积分的问题 对于二重积分的定义,我们并没有f(x,y)≥0的限.容易看出,当f(x,y)≥0时,二重积分在几何上就是以z=f(x,y)为曲顶，以(σ)为底且母线平行于z轴的曲顶柱体的体积。 上述就是二重积分的几何意义。
如果被积函数f(x,y)在积分区域(σ)上连续，那末二重积分必定存在。 二重积分的性质 (1).被积函数中的常数因子可以提到二重积分符号外面去. (2).有限个函数代数和的二重积分等于各函数二重积分的代数和. (3).如果把积分区域(σ)分成两个子域(σ1)与(σ2),即(σ)=(σ1)+(σ2),那末: (4).如果在(σ)上有f(x,y)≤g(x,y),那末: ≤ (5).设f(x,y)在闭域(σ)上连续,则在(σ)上至少存在一点(ξ,η),使 其中σ是区域(σ)的面积. 二重积分的计算法 直角坐标系中的计算方法 这里我们采取的方法是累次积分法。也就是先把x看成常量，对y进行积分，然后在对x进行积分，或者是先把y看成常量，对x进行积分，然后在对y进行积分。为此我们有积分公式，如下：
二重积分的计算方法
重庆三峡学院数学分析课程论文 二重积分的计算方法 院系数学与统计学院 专业数学与应用数学（师范） 姓名 年级 2010级 学号 指导教师刘学飞 2014年5月
二重积分的计算方法 (重庆三峡学院数学与统计学院10级数本1班) 摘 要 ：本文总结出了求二重积分的几种方法，比如用定义、公式、定理、性质求极限. 关键词 ：函数极限；计算方法；洛必达法则; 四则运算 引言 二重积分的概念和计算是多元函数微积分学的重要部分,在几何、物理、力学等方面有着重 要的应用.重积分是由一元函数积分推广而来的,但与一元函数相比,计算重积分的难度除了与被 积函数有关外,还与积分区域的特点有关,计算重积分的主要思想方法是化重积分为累次积分.求 二重积分的方法很多且非常灵活，本文归纳了二重积分计算的一些常见方法和技巧. 1. 预备知识 1.1二重积分的定义 设(),f x y 是定义在可求面积的有界区域D 上的函数. J 是一个确定的数,若对任给的正数 ε,总存在某个正数δ,使对于D 的任意分割T ,当它的细度T δ<时,属于T 的所有积分和都有 ()1 ,n i i i i f J ξησ ε=?-<∑, 则称(),f x y 在D 上可积,数J 称为函数(),f x y 在D 上的二重积分,记作(),D J f x y d σ= ??, 其中(),f x y 称为二重积分的被积函数, ,x y 称为积分变量, D 称为积分区域. 1.2二重积分的若干性质 1.21若(),f x y 在区域D 上可积, k 为常数,则(),kf x y 在D 上也可积,且 (),D kf x y d σ??(),D k f x y d σ=??. 1.22 若(),f x y ,(),g x y 在D 上都可积,则()(),,f x y g x y ±在D 上也可积,且 ()()[,,]D f x y g x y d σ±??()(),,D D f x y d g x y d σσ=±????.
详解Matlab求积分的各种方法
详解Matlab求积分的各种方法 一、符号积分由函数int来实现。 该函数的一般调用格式为： int(s)： 没有指定积分变量和积分阶数时，系统按findsym函数指示的默认变量对被积函数或符号表达式s求不定积分；int(s,v)： 以v为自变量，对被积函数或符号表达式s求不定积分；int(s,v,a,b)： 求定积分运算。 a,b分别表示定积分的下限和上限。 该函数求被积函数在区间[a,b]上的定积分。 a和b可以是两个具体的数，也可以是一个符号表达式，还可以是无穷(inf)。 当函数f关于变量x在闭区间[a,b]上可积时，函数返回一个定积分结果。 当a,b中有一个是inf时，函数返回一个广义积分。 当a,b中有一个符号表达式时，函数返回一个符号函数。 例： 求函数x^2+y^2+z^2的三重积分。 内积分上下限都是函数，对z积分下限是sqrt(x*y)，积分上限是x^2*y；对y积分下限是sqrt(x)，积分上限是x^2；对x的积分下限1，上限是2，求解如下： >>syms x y z %定义符号变 量>>F2=int(int(int(x^2+y^2+z^2,z,sqrt(x*y),x^2*y),y,sqrt(x),x^2),x,1,2) %注意定积分的书写格式F2 =57/-
/348075*2^(1/2)+14912/4641*2^(1/4)+64/225*2^(3/4) %给出有理数 解>>VF2=vpa(F2) %给出默认精度的数值解VF2 = 224.9 232805二、数值积分 1.数值积分基本原理求解定积分的数值方法多种多样，如简单的梯形法、辛普生(Simpson)?法、牛顿－柯特斯(Newton-Cotes)法等都是经常采用的方法。 它们的基本思想都是将整个积分区间[a,b]分成n个子区间[xi,xi+1]， i=1,2,…,n，其中x1=a，xn+1=b。 这样求定积分问题就分解为求和问题。 2.数值积分的实现方法基于变步长辛普生法，MATLAB给出了quad函数来求定积分。 该函数的调用格式为： [I,n]=quad('fname',a,b,tol,trace)基于变步长、牛顿－柯特斯(Newton-Cotes)法，MATLAB给出了quadl函数来求定积分。 该函数的调用格式为： [I,n]=quadl('fname',a,b,tol,trace)其中fname是被积函数名。 a和b分别是定积分的下限和上限。 tol用来控制积分精度，缺省时取tol= 0.0 01。 trace控制是否展现积分过程，若取非0则展现积分过程，取0则不展现，缺省时取trace= 0。
数值积分的matlab实现
实验10 数值积分 实验目的： 1．了解数值积分的基本原理； 2．熟练掌握数值积分的MATLAB 实现； 3．会用数值积分方法解决一些实际问题。 实验内容： 积分是数学中的一个基本概念，在实际问题中也有很广泛的应用。同微分一样，在《微积分》中，它也是通过极限定义的，由于实际问题中遇到的函数一般都以列表形式给出，所以常常不能用来直接进行积分。此外有些函数虽然有解析式，但其原函数不是初等函数，所以仍然得不到积分的精确值，如不定积分?1 0 d sin x x x 。这时我们一般考虑用数值方法计算其 近似值，称为数值积分。 10.1 数值微分简介 设函数()y f x =在* x 可导，则其导数为 h x f h x f x f h ) ()(lim )(**0* -+='→ （10.1） 如果函数()y f x =以列表形式给出（见表10-1），则其精确值无法求得，但可由下式求得其近似值 h x f h x f x f ) ()()(*** -+≈' （10.2） 表 10-1 一般的，步长h 越小，所得结果越精确。（10.2）式右端项的分子称为函数()y f x =在 *x 的差分，分母称为自变量在*x 的差分，所以右端项又称为差商。数值微分即用差商近似 代替微商。常用的差商公式为： 000()() ()2f x h f x h f x h +--'≈ （10.3） h y y y x f 243)(2 100-+-≈ ' （10.4）
h y y y x f n n n n 234)(12+-≈ '-- （10.5） 其误差均为2 ()O h ，称为统称三点公式。 10.2 数值微分的MATLAB 实现 MATLAB 提供了一个指令求解一阶向前差分，其使用格式为： dx=diff(x) 其中x 是n 维数组，dx 为1n -维数组[]21321,, ,n x x x x x x ---,这样基于两点的数值导 数可通过指令diff(x)/h 实现。对于三点公式，读者可参考例1的M 函数文件diff3.m 。 例1 用三点公式计算()y f x =在=x 1.0,1.2,1.4处的导数值，()f x 的值由下表给 解：建立三点公式的M 函数文件diff3.m 如下： function f=diff3(x,y) n=length(x);h=x(2)-x(1); f(1)=(-3*y(1)+4*y(2)-y(3))/(2*h); for j=2:n-1 f(j)=(y(j+1)-y(j-1))/(2*h); end f(n)=(y(n-2)-4*y(n-1)+3*y(n))/(2*h); 在MATLAB 指令窗中输入指令： x=[1.0,1.1,1.2,1.3,1.4];y=[0.2500,0.2268,0.2066,0.1890,0.1736];diff3(x,y) 运行得各点的导数值为：-0.2470，-0.2170，-0.1890，-0.1650，-0.0014。所以()y f x =在=x 1.0,1.2,1.4处的导数值分别为-0.2470，-0.1890和-0.0014。 对于高阶导数，MATLAB 提供了几个指令借助于样条函数进行求导，详细使用步骤如下： step1：对给定数据点（x,y ），利用指令pp=spline(x,y)，获得三次样条函数数据pp ，供后面ppval 等指令使用。其中，pp 是一个分段多项式所对应的行向量，它包含此多项式的阶数、段数、节点的横坐标值和各段多项式的系数。 step2：对于上面所求的数据向量pp ，利用指令[breaks,coefs,m,n]=unmkpp(pp)进行处理，生成几个有序的分段多项式pp 。 step3：对各个分段多项式pp 的系数，利用函数ppval 生成其相应导数分段多项式的系数，再利用指令mkpp 生成相应的导数分段多项式 step4：将待求点xx 代入此导数多项式，即得样条导数值。 上述过程可建立M 函数文件ppd.m 实现如下： function dy=ppd(pp) [breaks,coefs,m]=unmkpp(pp);
二重积分的计算方法
第二节 二重积分的计算法 教学目的：熟练掌握二重积分的计算方法 教学重点：利用直角坐标和极坐标计算二重积分 教学难点：化二重积分为二次积分的定限问题 教学内容： 利用二重积分的定义来计算二重积分显然是不实际的,二重积分的计算是通过两个定积分的计算(即二次积分)来实现的. 一、利用直角坐标计算二重积分 我们用几何观点来讨论二重积分的计算问题. 讨论中,我们假定 ； 假定积分区域可用不等式 表示, 其中, 在上连续. 据二重积分的几何意义可知,的值等于以为底,以曲面为顶的曲顶柱体的体积. 在区间上任意取定一个点,作平行于面的平面,这平面截曲顶柱体所得截面是一个以区间为底,曲线为曲边的曲边梯形,其面积为
一般地,过区间上任一点且平行于面的平面截曲顶柱体所得截面的面积为 利用计算平行截面面积为已知的立体之体积的方法,该曲顶柱体的体积为 从而有 (1) 上述积分叫做先对Y,后对X的二次积分,即先把看作常数,只看作的函数,对 计算从到的定积分,然后把所得的结果( 它是的函数 )再对从到计算定积分. 这个先对, 后对的二次积分也常记作 在上述讨论中,假定了，利用二重积分的几何意义,导出了二重积分的计算公式(1).但实际上,公式(1)并不受此条件限制,对一般的(在上连续),公式(1)总是成立的. 例如:计算 解: 类似地,如果积分区域可以用下述不等式 表示,且函数,在上连续,在上连续,则 (2)
显然,(2)式是先对,后对的二次积分. 二重积分化二次积分时应注意的问题 1、积分区域的形状 前面所画的两类积分区域的形状具有一个共同点: 对于I型(或II型)区域, 用平行于轴(轴 )的直线穿过区域内部,直线与区域的边界相交不多于两点. 如果积分区域不满足这一条件时,可对区域进行剖分,化归为I型(或II型)区域的并集. 2、积分限的确定 二重积分化二次积分, 确定两个定积分的限是关键.这里,我们介绍配置二 次积分限的方法 -- 几何法.画出积分区域的图形(假设的图形如下 ) 在上任取一点,过作平行于轴的直线,该直线穿过区域,与区域的边界有两个交 点与,这里的、就是将,看作常数而对积分时的下限和上限； 又因是在区间上任意取的,所以再将看作变量而对积分时,积分的下限为、上限为 . 例1计算,其中是由轴,轴和抛物线在第一象限内所围成的区域.
MATLAB计算积分
函数的积分和椭圆的周长 1．正弦函数的积分 [问题]求正弦函数从0到π的积分 y = sin x 当x = 0时，积分为0，画出积分的函数曲线。 [数学模型] 定积分的结果为 ππ00 sin d cos 2S x x x ==-=? 不定积分的结果为 sin d cos I x x x C ==-+? 其中C 是积分常量，由初始条件决定。当x = 0时，积分为I = 0，必有C = 1。结果为 I = -cos x + 1 [算法]根据积分的基本概念，将积分区域分为多份，用矩形法求曲线下的近似面积表示积分的近似值 1()n i i S f x x ==?∑ 矩形法的函数是sum(f)。 用梯形法求曲线下的近似面积表示积分的近似值 1 101[()()]2 n i i i S f x f x x -+==+?∑ 梯形法的函数是trapz(f)。 用数值积分的函数是quad 和quadl ，常用使用格式是 S = quad(f,a,b) 其中，f 表示被积函数，a 表示积分的下限，b 表示积分的下限。 用符号的函数是int ，常用使用格式是 S = int(f,a,b) [程序]zqy4_1.m 如下。 %正弦函数的积分 clear %清除变量 x=linspace(0,pi); %自变量向量 dx=x(2); %间隔 y=sin(x); %被积函数 s1=sum(y)*dx %矩形法积分 s2=trapz(y)*dx %梯形法积分 f=inline('sin(x)'); %被积的内线函数 s3=quad(f,0,pi) %数值定积分
s4=int('sin(x)',0,pi) %符号积分 sc1=cumsum(y)*dx; %矩形法累积积分(精度稍差) sc2=cumtrapz(y)*dx; %梯形法累积积分 figure %创建图形窗口 plot(x,-cos(x)+1,x,sc1,'.',x,sc2,'o') %画解析式和矩阵法以及梯形法积分曲线 s=int('sin(x)') %符号积分 sc3=subs(s,'x',x); %替换数值求符号积分的值 C=-sc3(1) %求积分常数 hold on %保持图像 plot(x,sc3+C,'c*') %画符号法积分曲线 grid on %加网格 fs=16; %字体大小 xlabel('\itx','FontSize',fs) %横坐标 ylabel('\intsin\itx\rmd\itx','FontSize',fs)%纵坐标 title('正弦函数的积分','FontSize',fs) %标题 legend('解析解','矩形法','梯形法','符号法')%图例 zqy4.1图 zqy4.2图 2.三角函数和指数的积分 [问题]求如下函数的积分 y = e ax sin bx 其中a = 0.5，b = 2。积分下限为0。画出积分的函数曲线。 [数学模型] 设 11e sin d sin de {e sin e cos d }ax ax ax ax I bx x bx bx b bx x a a == =-??? 11{e sin cos de }{e sin [e cos e sin d ]}ax ax ax ax ax b b bx bx bx bx b bx x a a a a =-=-+?? 因此不定积分为 221e (sin cos )ax I a bx b bx C a b =-++ 当x = 0时，I 应该为零，所以 22b C a b =+
MATLAB实验三-定积分的近似计算
实验三定积分的近似计算 一、问题背景与实验目的 利用牛顿—莱布尼兹公式虽然可以精确地计算定积分的值，但它仅适用于被积函数的原函数能用初等函数表达出来的情形．如果这点办不到或者不容易办到，这就有必要考虑近似计算的方法．在定积分的很多应用问题中，被积函数甚至没有解析表达式，可能只是一条实验记录曲线，或者是一组离散的采样值，这时只能应用近似方法去计算相应的定积分． 本实验将主要研究定积分的三种近似计算算法：矩形法、梯形法、抛物线法．对于定积分的近似数值计算，Matlab有专门函数可用． 二、相关函数（命令）及简介 1．sum(a)：求数组a的和． 2．format long：长格式，即屏幕显示15位有效数字． （注：由于本实验要比较近似解法和精确求解间的误差，需要更高的精度）．3．double()：若输入的是字符则转化为相应的ASCII码；若输入的是整型数值则转化为相应的实型数值． 4．quad()：抛物线法求数值积分． 格式： quad(fun,a,b) ，注意此处的fun是函数，并且为数值形式的，所以使用*、/、^等运算时要在其前加上小数点，即 .*、./、.^等．例：Q = quad('1./(x.^3-2*x-5)',0,2); 5．trapz()：梯形法求数值积分． 格式：trapz(x,y) 其中x为带有步长的积分区间；y为数值形式的运算（相当于上面介绍的函数fun） 例：计算 0sin()d x x π ? x=0:pi/100:pi;y=sin(x); trapz(x,y) 6．dblquad()：抛物线法求二重数值积分． 格式：dblquad(fun,xmin,xmax,ymin,ymax)，fun可以用inline定义，也可以通过某个函数文件的句柄传递． 例1：Q1 = dblquad(inline('y*sin(x)'), pi, 2*pi, 0, pi) 顺便计算下面的Q2，通过计算，比较Q1 与Q2结果（或加上手工验算），找出积分变量x、y的上下限的函数代入方法． Q2 = dblquad(inline('y*sin(x)'), 0, pi, pi, 2*pi)例2：Q3 = dblquad(@integrnd, pi, 2*pi, 0, pi) 这时必须存在一个函数文件integrnd.m：
二重积分的概念与性质教案
7.1二重积分的基本概念（教案） 主讲人：孙杰华 教学目的：理解二重积分的概念、性质 教学重难点：二重积分的概念、二重积分的几何意义. 教学方法：讲授为主 教学内容： 一、二重积分的概念 1．曲顶柱体的体积 设有一空间立体Ω,它的底是xoy 面上的有界区域D ,它的侧面是以D 的边界曲线为准线,而母线平行于z 轴的柱面,它的顶是曲面(.)z f x y =，称这种立体为曲顶柱体. 与求曲边梯形的面积的方法类似，我们可以这样来求曲顶柱体的体积V : (1)用任意一组曲线网将区域D 分成n 个小区域1σ?，2σ?， ，n σ?，以这些小区 域的边界曲线为准线,作母线平行于z 轴的柱面,这些柱面将原来的曲顶柱体Ω分划成n 个小曲顶柱体1?Ω，2?Ω， ，n ?Ω. （假设i σ?所对应的小曲顶柱体为i ?Ω,这里i σ?既代表第i 个小区域,又表示它的面积值, i ?Ω既代表第i 个小曲顶柱体,又代表它的体积值.)，从而1 n i i V ==?Ω∑． 图7.1 (2)由于(,)f x y 连续,对于同一个小区域来说,函数值的变化不大．因此,可以将小曲顶柱体近似地看作小平顶柱体,于是
(,),((,))i i i i i i i f ξησξησ?Ω≈??∈?. (3)整个曲顶柱体的体积近似值为 1 (,)n i i i i V f ξησ=≈?∑. (4)为得到的精确值,只需让这个小区域越来越小,即让每个小区域向某点收缩．为此,我们引入区域直径的概念: 一个闭区域的直径是指区域上任意两点距离的最大者. 所谓让区域向一点收缩性地变小,意指让区域的直径趋向于零. 设n 个小区域直径中的最大者为λ,则 1 lim (,),(,)n i i i i i i i V f λξησξησ→==??∈?∑. 2．二重积分的定义 设(),f x y 是闭区域D 上的有界函数, 将区域D 分成个小区域 12,,,,n σσσ??? 其中,i σ?既表示第i 个小区域,也表示它的面积, i λ表示它的直径. 1max{}(,)i i i i i n λλξησ≤≤=?∈?， 作乘积(,)(1,2 ,)i i i f i n ξησ?=， 作和式 1 (,)n i i i i f ξησ =?∑， 若极限()0 1 lim ,n i i i i f λξησ →=?∑存在,则称此极限值为函数(),f x y 在区域D 上的二重积分,记 作 (),D f x y d σ??.即 (),D f x y d σ=??()0 1 lim ,n i i i i f λξησ →=?∑． 其中:(),f x y 称之为被积函数,(),f x y d σ称之为被积表达式,d σ称之为面积元素, ,x y 称之为积分变量,D 称之为积分区域. V n
[整理]Matlab积分.
一.数值积分的实现方法 1．变步长辛普生法 基于变步长辛普生法，MA TLAB给出了quad函数来求定积分。该函数的调用格式为：[I,n]=quad('fname',a,b,tol,trace) 其中fname是被积函数名。a和b分别是定积分的下限和上限。tol用来控制积分精度，缺省时取tol=0.001。trace控制是否展现积分过程，若取非0则展现积分过程，取0则不展现，缺省时取trace=0。返回参数I即定积分值，n为被积函数的调用次数。 例8-1 求定积分。 (1) 建立被积函数文件fesin.m。 function f=fesin(x) f=exp(-0.5*x).*sin(x+pi/6); (2) 调用数值积分函数quad求定积分。 [S,n]=quad('fesin',0,3*pi) S = 0.9008 n = 77 2．牛顿－柯特斯法 基于牛顿－柯特斯法，MA TLAB给出了quad8函数来求定积分。该函数的调用格式为：[I,n]=quad8('fname',a,b,tol,trace) 其中参数的含义和quad函数相似，只是tol的缺省值取10-6。?该函数可以更精确地求出定积分的值，且一般情况下函数调用的步数明显小于quad函数，从而保证能以更高的效率求出所需的定积分值。 (1) 被积函数文件fx.m。 function f=fx(x) f=x.*sin(x)./(1+cos(x).*cos(x)); (2) 调用函数quad8求定积分。 I=quad8('fx',0,pi) I = 2.4674 分别用quad函数和quad8函数求定积分的近似值，并在相同的积分精度下，比较函数的调用次数。 调用函数quad求定积分： format long; fx=inline('exp(-x)'); [I,n]=quad(fx,1,2.5,1e-10) I = 0.28579444254766 n = 65 调用函数quad8求定积分： format long; fx=inline('exp(-x)'); [I,n]=quad8(fx,1,2.5,1e-10) I = 0.28579444254754 n = 33
二重积分的概念及计算法(一)
习题9-1，9-2 二重积分的概念及计算法（一） 1．填空题： （1）由二重积分的几何意义得 ∫∫≤+=??122221y x d y x σ . （2）根据二重积分的性质，比较下列积分的大小： ① ，其中是三角形区域，三顶点为（1,0），（1,1），（2,0），则 ∫∫+=D d y x I σ)ln(1∫∫ +=D d y x I σ22)][ln(D 1I 2I . ②，，其中是由∫∫++=D d y x I σ21)1(∫∫ ++=D d y x I σ32)1(D x 轴与直线围成的区域，则 1,0?==+x y x 1I 2I . （3）化二重积分为两种不同次序下的二次积分，其中是直线D 2,==x x y 及双曲线)0(1f x x y =所围成的闭区域，= ∫∫d y x f σ),(D = （4）①交换积分次序： ∫∫??=22221),(x x x dy y x f dx ②交换积分次序： ∫∫∫∫?=+y y dx y x f dy dx y x f dy 20313010),(),( 2．利用二重积分的性质，估计积分的值： ∫∫++=D d y x I σ)94(22，其中是圆形闭区域：. D 422≤+y x 3．计算下列二重积分： （1）∫∫+= D d x x y I σ2)1(cos ，其中是顶点分别为（0,0）,(1,0),(1,2)和（0,1）的梯形闭区域. D （2），其中是由∫∫+=D y x d e I σD 1≤+y x 所确定的闭区域. 4．计算二次积分∫∫101dx e dy y x y . 5．交换积分次序，证明： ∫∫∫???=a y a x a m x a m dx x f e x a dx x f e dy 000)()()()()(. 6．设平面薄片所占的闭区域是由直线D x y y x ==+,2和x 轴所围成，它的面密度
二重积分计算中的积分限的确定
二重积分计算中积分限的确定 摘要:二重积分计算中积分限的确定对于初学者是一个重点更是一个难点.本文旨在介绍一种二重积分计算中确定积分限的简单易行的方法. 关键词:二重积分累次积分积分限积分次序 引言：高等数学学习过程中，二重积分计算是个难点。原因在于将二重积分化为累次积分时，对于积分限的确定学生难以掌握。本人结合自己的教学过程和自己的学习体会总结出一个口诀，发现在教学过程中效果不错可以很好的帮助学生解决这一难题。 1．高等数学中计算二重积分的方法 在高等数学课本中,在直角坐标系下计算二重积分的步骤为:]1[。 (1)画出积分区域 (2)确定积分区域是否为X-型或Y-型区域,如既不是X-型也不是Y-型区域,则要将 积分区域化成几个X-型和Y-型区域,并用不等式组表示每个X-型和Y-型区域. (3)用公式化二重积分为累次积分. (4)计算累次积分的值. 在教学的过程中我发现学生对于此种方法掌握的很不好,尤其是在第二步中,确定积分区域从而确定累次积分的积分限是一个薄弱环节.下面就本人在教学中的体会谈谈在这方面的一点心得. 2．教学过程中总结的方法 本人的心得可用下面的口诀概括:后积先定限,限内画条线,先交下限取,后交上限见.下面简单解释一下该口诀,然后以具体的例题加以说明.在将二重积分转化为累次积分的时候对于两个积分变量必然会有个先后顺序,这就要求对后积分的那个变量我们要根据积分区域确定其上下限(所谓确定是指根据积分区域图将其上下限定为常数).确定了这个变量的上下限以后,我们在其上下限内画一条和上下限平行的直线,该直线沿着坐标轴的正方向画过来,这样该直线如果和积分区域总是有两个交点,先交的即为另一个积分变量的积分下限,后交的即为其积分上限. 3．例题解析 例1 计算?? D xydxdy,其中D是由直线x y y x= = =,1 ,2所围成的区域. 解:作出积分区域D的图形 x 页脚内容1