高等数学-第七版-课件-10-6 定积分的近似计算
高等数学 定积分
第五章 定积分第一节 定积分的概念第二节 定积分的性质和中值定理第三节 微积分基本公式第四节 定积分的换元法第五节 定积分的分部积分法第六节 定积分的近似计算第七节 广义积分问题的提出定积分的定义 几何意义定积分存在定理第一节 定积分的概念abxyo?=A 曲边梯形由连续曲线实例1 (求曲边梯形的面积))(x f y =)0)((≥x f 、x 轴与两条直线a x =、b x =所围成.一、问题的提出)(x f y =ab xyoab x yo用矩形面积近似取代曲边梯形面积显然,小矩形越多,矩形总面积越接近曲边梯形面积.(四个小矩形)(九个小矩形)观察下列演示过程,注意当分割加细时,矩形面积和与曲边梯形面积的关系.观察下列演示过程,注意当分割加细时,矩形面积和与曲边梯形面积的关系.观察下列演示过程,注意当分割加细时,矩形面积和与曲边梯形面积的关系.观察下列演示过程,注意当分割加细时,矩形面积和与曲边梯形面积的关系.观察下列演示过程,注意当分割加细时,矩形面积和与曲边梯形面积的关系.曲边梯形如图所示,,],[1210b x x x x x a b a n n =<<<<<=- 个分点,内插入若干在区间a bxyoi ξi x 1x 1-i x 1-n x ;],[],[11---=∆i i i i i x x x x x n b a 长度为,个小区间分成把区间形面积,曲边梯形面积用小矩上任取一点在每个小区间i i i x x ξ-],[1ii i x f A ∆ξ≈)(:))(],[(1近似为高为底,以i i i f x x ξ-(1)分割(2)近似ini i x f A ∆≈∑=)(1ξ曲边梯形面积的近似值为ini i x f A ∆=∑=→)(lim 10ξλ时,趋近于零即小区间的最大长度当分割无限加细)0(},,max{,21→∆∆∆=λλn x x x 曲边梯形面积为(3)求和(4)取极限实例2 (求变速直线运动的路程)设某物体作直线运动,已知速度)(t v v =是时间间隔],[21T T 上t 的一个连续函数,且0)(≥t v ,求物体在这段时间内所经过的路程.思路:把整段时间分割成若干小段,每小段上速度看作不变,求出各小段的路程再相加,便得到路程的近似值,最后通过对时间的无限细分过程求得路程的精确值.(1)分割212101T t t t t t T n n =<<<<<=- 1--=∆i i i t t t ii i t v s ∆≈∆)(τ部分路程值某时刻的速度(3)求和ii ni t v s ∆≈∑=)(1τ(4)取极限},,,max{21n t t t ∆∆∆= λini i t v s ∆=∑=→)(lim 10τλ路程的精确值(2)近似设函数)(x f 在],[b a 上有界,记},,,max{21n x x x ∆∆∆= λ,如果不论对],[b a 在],[b a 中任意插入若干个分点bx xx x x a nn =<<<<<=-121把区间],[b a 分成n 个小区间,各小区间的长度依次为1--=∆i i i x x x ,),2,1( =i ,在各小区间上任取一点i ξ(i i x ∆∈ξ),作乘积i i x f ∆)(ξ ),2,1( =i 并作和i i ni x f S∆=∑=)(1ξ,二、定积分的定义定义怎样的分法,⎰==ba I dx x f )(ii ni x f ∆∑=→)(lim 10ξλ被积函数被积表达式积分变量积分区间],[b a 也不论在小区间],[1i i x x -上点i ξ怎样的取法,只要当0→λ时,和S 总趋于确定的极限I ,我们称这个极限I 为函数)(x f 在区间],[b a 上的定积分,记为积分上限积分下限积分和几点说明:(1) 定积分是一个数值,它仅与被积函数及积分区间有关,⎰b a dx x f )(⎰=b a dt t f )(⎰=ba duu f )(而与积分变量的字母无关.)( ,)()( 2⎰⎰⎰=-=aaabbadx x f dx x f dx x f 规定:)(.],[)(],[)( 3的取法无关的分法及的和式的极限与所表示上可积,则在区间若)(i bab a dx x f b a x f ξ⎰,0)(≥x f ⎰=ba Adx x f )(曲边梯形的面积,0)(≤x f ⎰-=ba Adx x f )(曲边梯形的面积的负值a b xyo)(x f y =AxyoabA -)(x f y =三、定积分的几何意义1A 2A 3A 4A 4321)(A A A A dx x f ba ⎰=-+-,],[)(变号时在区间b a x f 三、定积分的几何意义.)(是面积的代数和⎰badx x f几何意义:积取负号.轴下方的面在轴上方的面积取正号;在数和.之间的各部分面积的代直线的图形及两条轴、函数它是介于x x b x a x x f x ==,)(++--当函数)(x f 在区间],[b a 上连续时,定理1定理2 设函数)(x f 在区间],[b a 上有界,且只有有限个间断点,则)(x f 在四、定积分的存在定理区间],[b a 上可积.例1 利用定义计算定积分.12dx x ⎰解将]1,0[n 等分,分点为nix i =,(n i ,,2,1 =)小区间],[1i i x x -的长度nx i 1=∆,(n i ,,2,1 =)取i i x =ξ,(n i ,,2,1 =)i i n i x f ∆∑=)(1ξi i ni x ∆=∑=21ξ,12i ni ix x ∆=∑=.,102的选取无关及法故和式极限与区间的分可积因为i dx x ξ⎰n n i ni 121⋅⎪⎭⎫ ⎝⎛=∑=∑==n i i n 12316)12)(1(13++⋅=n n n n ,121161⎪⎭⎫ ⎝⎛+⎪⎭⎫ ⎝⎛+=n n ∞→⇒→n 0λdx x ⎰102i i ni x ∆=∑=→210lim ξλ⎪⎭⎫ ⎝⎛+⎪⎭⎫ ⎝⎛+=∞→n n n 121161lim .31= 几何上是曲线y=x 2,直线x=1及x 轴围成的曲边三角形面积.例2 利用定义计算定积分.121dx x⎰解在]2,1[中插入分点 12,,,-n q q q ,典型小区间为],[1ii q q -,(n i ,,2,1 =)小区间的长度)1(11-=-=∆--q qq q x i i i i ,取1-=i i qξ,(n i ,,2,1 =)i i ni x f ∆∑=)(1ξi ni ix ∆=∑=11ξ)1(1111-=-=-∑q q q i ni i ∑=-=ni q 1)1()1(-=q n 取2=nq即nq 12=),12(1-=n n )12(lim 1-+∞→xx x x xx 112lim1-=+∞→,2ln =)12(lim 1-∴∞→nn n ,2ln =dx x ⎰211i ni ix ∆=∑=→101lim ξλ)12(lim 1-=∞→n n n .2ln =i i ni x f ∆∑=)(1ξ原式⎥⎦⎤⎢⎣⎡π+π-++π+π=∞→n n n n n n n nsin )1(sin 2sin sin 1lim π=∑=∞→n i n n i n 1sin 1lim n n i ni n π⋅⎪⎭⎫ ⎝⎛ππ=∑=∞→1sin lim 1.sin 10⎰ππ=xdx ix ∆i ξ例3:将下列和式极限表示成定积分.⎥⎦⎤⎢⎣⎡-+++∞→n n n n n n πππ)(sin sin sin lim121 :五、小结1.定积分的实质:特殊和式的极限.2.定积分的思想和方法:分割化整为零求和积零为整取极限精确值——定积分求近似以直(不变)代曲(变)取极限Z .思考n n n n f n f n f ⎪⎭⎫ ⎝⎛⎪⎭⎫ ⎝⎛⋅⎪⎭⎫ ⎝⎛∞→ 21lim 试证.1)(ln ⎰=dxx f e 2:将和式极限,表示成定积分.⎥⎦⎤⎢⎣⎡-++-+-∞→2222241241141lim n n n n n 证明n n n n f n f n f ⎪⎭⎫⎝⎛⎪⎭⎫ ⎝⎛⋅⎪⎭⎫⎝⎛∞→ 21lim ⎪⎭⎫⎝⎛⎪⎭⎫ ⎝⎛⎪⎭⎫ ⎝⎛⋅⎪⎭⎫ ⎝⎛∞→=n n n n f n f n f e21lim ln n n n n f n f n f ⎪⎭⎫⎝⎛⎪⎭⎫ ⎝⎛⋅⎪⎭⎫ ⎝⎛∞→ 21lim 试证.1)(ln ⎰=dx x f e 利用对数的性质得⎪⎭⎫⎝⎛∑==∞→n i f n ni n e1ln 1lim n n i f ni n e1ln lim 1⋅⎪⎭⎫ ⎝⎛∑==∞→ 指数上可理解为:)(ln x f 在]1,0[区间上的一个积分和.分割是将]1,0[n 等分分点为nix i =,(n i ,,2,1 =)⎪⎭⎫⎝⎛⎪⎭⎫ ⎝⎛⎪⎭⎫ ⎝⎛⋅⎪⎭⎫ ⎝⎛∞→=n n n n f n f n f e21ln lim 极限运算与对数运算换序得nn i f n i n 1ln lim 1⋅⎪⎭⎫ ⎝⎛∑=∞→⎰=10)(ln dx x f 故nn n n f n f n f ⎪⎭⎫ ⎝⎛⎪⎭⎫ ⎝⎛⋅⎪⎭⎫ ⎝⎛∞→ 21lim.10)(ln ⎰=dxx f e 因为)(x f 在区间]1,0[上连续,且0)(>x f 所以)(ln x f 在]1,0[上有意义且可积 ,2:将和式极限,表示成定积分.⎥⎦⎤⎢⎣⎡-++-+-∞→2222241241141lim n n n n n ⎰∑-=-=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-++-+-=⎥⎦⎤⎢⎣⎡-++-+-=∞→∞→∞→1021222222222411)(41lim )(41)2(41)1(411lim 41241141lim dxx n ni n n n n n n n n n n i n n n 解第二节 定积分的性质、中值定理1.定积分性质2.中值定理对定积分的补充规定:(1)当b a =时,0)(=⎰ba dx x f ;(2)当b a >时,⎰⎰-=abb adx x f dx x f )()(.说明 在下面的性质中,假定定积分都存在,且不考虑积分上下限的大小.一、定积分性质和中值定理证⎰±ba dxx g x f )]()([i i i ni x g f ∆±=∑=→)]()([lim 10ξξλi i ni x f ∆=∑=→)(lim 10ξλii ni x g ∆±∑=→)(lim 10ξλ⎰=ba dx x f )(.)(⎰±ba dx x g ⎰±b a dx x g x f )]()([⎰=b a dx x f )(⎰±ba dx x g )(.(此性质可以推广到有限多个函数作和的情况)性质1⎰⎰=ba b a dx x f k dx x kf )()( (k 为常数).证⎰ba dx x kf )(ii ni x kf ∆=∑=→)(lim 10ξλi i n i x f k ∆=∑=→)(lim 1ξλii ni x f k ∆=∑=→)(lim 10ξλ.)(⎰=ba dx x f k 性质2⎰ba dx x f )(⎰⎰+=bcca dx x f dx x f )()(.补充:不论 的相对位置如何, 上式总成立.c b a ,,例 若,c b a <<⎰c a dx x f )(⎰⎰+=cb b a dx x f dx x f )()(⎰b a dx x f )(⎰⎰-=cb c a dxx f dx x f )()(.)()(⎰⎰+=bc ca dx x f dx x f (定积分对于积分区间具有可加性)假设bc a <<性质3dx b a ⋅⎰1dx ba⎰=a b -=.则0)(≥⎰dx x f ba. )(b a <证,0)(≥x f ,0)(≥ξ∴i f ),,2,1(n i =,0≥∆i x ,0)(1≥∆ξ∴∑=i i ni x f },,,max{21n x x x ∆∆∆= λi i ni x f ∆∴∑=→)(lim 1ξλ.0)(⎰≥=ba dx x f 性质4性质5如果在区间],[b a 上0)(≥x f ,例1 比较积分值dx e x⎰-20和dx x ⎰-20的大小.解令,)(x e x f x -=]0,2[-∈x ,0)(>x f ,0)(02>-∴⎰-dx x exdx ex⎰-∴2,02dx x ⎰->于是dx e x ⎰-2.20dx x ⎰-<性质5的推论:证),()(x g x f ≤ ,0)()(≥-∴x f x g ,0)]()([≥-∴⎰dx x f x g ba ,0)()(≥-⎰⎰ba ba dx x f dx x g 于是 dx x f ba ⎰)( dx x g ba ⎰≤)(.则dx x f ba ⎰)( dx x g ba ⎰≤)(. )(b a <如果在区间],[b a 上)()(x g x f ≤,(1)dx x f b a ⎰)(dx x f ba⎰≤)(.)(b a <证,)()()(x f x f x f ≤≤- ,)()()(dx x f dx x f dx x f ba ba ba ⎰⎰⎰≤≤-∴即dx x f ba ⎰)(dx x f ba⎰≤)(.说明: 可积性是显然的.|)(x f |在区间],[b a 上的性质5的推论:(2)设M 及m 分别是函数证,)(M x f m ≤≤ ,)(⎰⎰⎰≤≤∴ba ba b a Mdx dx x f dx m ).()()(a b M dx x f a b m ba -≤≤-⎰(此性质可用于估计积分值的大致范围)则 )()()(a b M dx x f a b m ba -≤≤-⎰.)(x f 在区间],[b a 上的最大值及最小值,性质6例2 估计积分dx x⎰π+03sin 31值的范围.解,sin 31)(3xx f +=],,0[π∈∀x ,1sin 03≤≤x ,31sin 31413≤+≤x ,31sin 31410030dx dx x dx ⎰⎰⎰πππ≤+≤.3sin 31403π≤+≤π∴⎰πdx x例3 估计积分dx xx⎰ππ24sin 值的范围.解,sin )(xx x f =2sin cos )(x x x x x f -='2)tan (cos x x x x -=⎥⎦⎤⎢⎣⎡∈2,4ππx ,0<)(x f 在]2,4[ππ上单调下降,,22)4(π=π=f M ,2)2(π=π=f m ,442π=π-π=-a b ,422sin 4224π⋅π≤≤π⋅π∴⎰ππdx x x .22sin 2124≤≤∴⎰ππdx x x 如果函数)(x f 在闭区间],[b a 上连续,上的平均值在],[)()(1b a x f dxx f a b ba⎰-则在积分区间],[b a 上至少存在一个点 ξ,使dx x f b a ⎰)())((a b f -=ξ. )(b a ≤≤ξ性质7(定积分中值定理)积分中值公式证Mdx x f a b m ba≤-≤∴⎰)(1)()()(a b M dx x f a b m ba -≤≤-⎰ 由闭区间上连续函数的介值定理知在区间],[b a 上至少存在一个点 ξ,)(1)(⎰-=ξbadx x f a b f dx x f ba ⎰)())((ab f -=ξ.)(b a ≤≤ξ即在区间],[b a 上至少存在一个点ξ,1. 积分中值公式的几何解释:xyoa b ξ)(ξf 使得以区间],[b a 为以曲线)(x f y =底边,为曲边的曲边梯形的面积等于同一底边而高为)(ξf 的一个矩形的面积。
高等数学第七版上册总复习PPT
第一页,共50页。
(一)函数
根本初等函数 复合函数
函数的定义
反函数
反函数与直接 函数之间关系
函数
的性质
奇偶性 单调性 有界性 周期性
第二页,共50页。
(二)数列极限
第三页,共50页。
第四页,共50页。
第五页,共50页。
(三)函数极限
第六页,共50页。
第七页,共50页。
振荡型
第十七页,共50页。
第十八页,共50页。
第十九页,共50页。
第二十页,共50页。
第二章 导数与微分
第二十一页,共50页。
第二十二页,共50页。
第二十三页,共50页。
第二十四页,共50页。
第二十五页,共50页。
第二十六页,共50页。
第二十七页,共50页。
第三章 中值定理和导数的应用
第四十八页,共50页。
第四十九页,共50页。
谢谢各位的聆听
第五十页,共50页。
第八页,共50页。
第九页,共50页。
第十页,共50页。
第十一页,共50页。
第十二页,共50页。
第十三页,共50页。
第十四页,共50页。
(三)连续与间断
第十五页,共50页。
第十六页,共50页。
第
y
一
可去型
类
间
断
点
o x0
x
y
第 二 类 间 断o 点
x0
x
无穷型
y
跳跃型
o
x0
x
y
o
x
第四十页,共50页。
四 、有理函数与可化为有理函数的积分
〔1〕有理函数的积分
《高数定积分》课件
05
广义积分及其收敛性判别法
广义积分的概念及分类
广义积分的定义
广义积分是相对于正常积分而言的一种特殊积分,其积分区间可能包含无穷大或者无界 函数。
广义积分的分类
根据被积函数和积分区间的不同,广义积分可分为无穷限广分的收敛性判别法
比较判别法
通过比较被积函数与已知收敛或发散的函数,来判断广义积分的收敛性。
换元法求解定积分
01
换元法的基本思想
通过变量代换简化定积分的计算 。
02
常见的换元方法
03
换元法的注意事项
三角函数代换、倒代换、根式代 换等。
代换后需调整积分上下限,并验 证代换的可行性。
分部积分法求解定积分
分部积分法的基本思想
将复杂函数拆分为简单函数 进行积分。
常见的分部积分公式
幂函数与三角函数、幂函数 与指数函数、幂函数与对数 函数等。
06
定积分在经济学等领域的应用
由边际函数求原经济函数
边际函数与定积分的关系
边际函数描述的是经济量变化的瞬时速率,而定积分则可用于求取原经济函数,即总量 函数。
求原经济函数的步骤
首先确定边际函数的表达式,然后根据定积分的定义,对边际函数进行积分,得到原经 济函数的表达式。
示例
已知某产品的边际收益函数为MR(q),通过对其进行定积分,可以得到总收益函数 TR(q)。
曲线的长度、图形的面积等。
THANKS
感谢观看
原函数与不定积分概念
原函数定义
原函数是指一个函数的导数等于给定函数的函数。根据微积分基本定理,不定积分就是求原函数的过 程。
不定积分性质
不定积分具有线性性质、常数倍性质和积分区间可加性。这些性质在求解复杂函数的定积分时非常有 用。
最新同济大学高等数学第七版上册定积分精品课件
则有
b
a f ( x)dx f [(t)] (t)dt
b
a f ( x)dx f [(t)] (t)dt
证 因为 f ( x) 在 [a, b] 上连续,故原函数存在,设 F( x) 是 f (x) 的一个原函数,则有
f [(t)](t)dt
f [(t)]d(t)
F[(t)] F[( )] F[( )]
2
三、小结
1、使用定积分的换元法时要注意积分限的对 应。
2、不引入新的变量记号,积分限不变;引入 新的变量记号,积分限跟着变。
3、定积分分部积分公式的用法与不定积分分 部积分公式的用法类似。
作业
P254 1 (4) , (10) , (16) ,(24) ; 3 ; 6; 7 (4), (9), (10)
2 2arctant 1 2 .
0
2
例7
设
f
(x)
12xx, 1 x
,
x0 x0,
2
求 f ( x 1)dx . 0
解 令 x1 t,
原式
1
f (t)dt
1
f ( x)dx
1
1
1
2xdx
0 1 x dx
0
1 1 x
x2 1
0
(1
2 ) dx
0
1
1 x
1 1 2 ln(1 x) 0 2 ln 2 . 1
T f ( x)dx .
a
0
aT
证 a f ( x)dx
0
T
aT
a f ( x)dx 0 f ( x)dx T f ( x)dx ,
aT
f ( x)dx
xT t
定积分近似计算方法
定积分的近似计算方法摘要 本文主要讨论了一元函数常见的数值积分方法,例如插值型求积公式、龙贝格求积公式、高斯求积公式等近似计算方法,在用这些方法计算定积分时,会产生一些误差,为了减少误差, 可以利用复化求积公式、复化高斯公式等.本文围绕这些方法,系统介绍它们的计算公式以及截断误差,并用例题分析它们产生误差的大小、计算量等.关键词 插值型积分 龙贝格积分 高斯积分 误差分析 近似计算1引言在计算定积分的值()b aI f x dx =⎰时,常常根据微积分学基本定理求出)(x f 的一个原函数)(x F ,再用牛顿-莱布尼茨公式求的积分,()()()baI f x dx F b F a ==-⎰.但在实际应用中,这种方法只限于解决一小部分定积分的求值问题.当函数没有具体表达式,只是一些实验测得数据形成的表格或图形或者是()F x 无法用初等函数表示,例如,2bx ae dx ⎰,2sin ba x dx ⎰等等,这就需要我们用一些近似方法求的积分值.与数值积分一样,把积分区间细分,在每个小区间上,找到简单函数)(x ϕ来近似代替()f x ,且()bax dx ϕ⎰的值容易求的.这样就把计算复杂的()baf x dx ⎰转化为求简单的积分值()bax dx ϕ⎰.因此,定积分的近似计算实质上是就是被积函数的近似计算问题.2常见数值方法 2.1牛顿-科茨数值方法牛顿-科茨求积公式是求积节点等距离分布的插值型求积公式.利用插值多项式来构造数值积分公式是最常用、最基本的方法,具体做法是:给定区间[,]a b 上一组节点01...n a x x x b =<<<=,以及节点处函数()(0,1,2,i f x i n =,作()f x 的n 次拉格朗日多项式()()()nn i i i x f x l x ϕ==∑,其中 011011()()()()()()()()()i i n i i i i i i i n x x L x x x x L x x l x x x L x x x x L x x -+-+----=----,将插值公式(1)1()()()()(1)!n n n f f x x x n ξϕω++=++. 其中 1012()()()()()n n x xx x xx x L x x ω+=----,[,]a b ξ∈,依赖于变量x , 上式积分得(1)1()()()()(1)!n bb bn n aa af f x dx x dx x dx n ξϕω++=++⎰⎰⎰(1)(1)0()()()()(1)!n nb biiin aai f f x l x dx x dx n ξω++==++∑⎰⎰(1)(1)0()()()()(1)!n nb bi i n aai f f x b l x dx x dxn ξω++==++∑⎰⎰若记 (),(0,1,2,bi ia A l x dx i ==⎰….. )n (1)(1)1()[]()(1)!n bn af R f x dxn ξω++=+⎰, (2)则有()()[]nbi i ai f x dx A f x R f ==+∑⎰(3)称式(3)为插值求型公式,其中(0,1,2,i A i =…. )n 与()f x 无关,叫求积系数, i x 为求积节点,[]R f 为求积公式余项,其中求积系数由(1)决定.2.1.1梯形求积公式1梯形公式当插值节点01,x x 分别选取区间端点,a b 时,由式(3)分别求出求积系数10012bb aa x x xb b aA dx dx x x a b ---===--⎰⎰,01102bb aa x x x ab a A dx dx x x b a ---===--⎰⎰.从而的求积公式()[()()]2bab af x dx f a f b -≈+⎰. (4) 称求积公式(4)为梯形求积公式,简称梯形公式.2梯形公式截断误差: 3*()[](),12b a R f f ξ-''=- *[,]a b ξ∈. (5) 3梯形求积公式的代数精度:1 当()1f x =时,式(5)中 1(1)2bab adx b a x b a -=-=+=-⎰. 精确成立.2.1.2 辛普森求积公式1辛普森求积公式当选取节点为012,,2a bx a x x b +===时,由式(1)求下列求积系数 1200102()()()()2()()6()()2b b a a a b x x b x x x x b a A dx dx a b x x x x a a b +-----===+----⎰⎰,0211002()()()()2()()()3()()22bb aa x x x x x a xb b a A dx dx a b a b x x x x a b -----===++----⎰⎰.0122021()()()()2()()6()()22b b a a a bx a x x x x x b a A dx dx a b a b x x x x a b +-----===++----⎰⎰ .从而求积公式()[()4()()]62bab a a bf x dx f a f f b -+≈++⎰. (6)称式(6)为抛物线积分公式或辛普森积分公式.2抛物线求积公式误差估计定理1.若()f x 在[,]a b 上有四阶连续导数,则抛物线公式(6)的余项为:5(4)**()[](),[,]2880b a R f f a b ξξ--=∈. (7) 3抛物线公式的代数精度为3.易验证,当23()1,,,f x x x x =时,式(6)精确成立,而当4()f x x =时,式(6)不能精确成立.2.1.3 牛顿-科茨公式1牛顿-科茨公式在等距离节点i x a ih =+下,其中(0,1,2b ah i n-==…. )n .作为变量替换x a th =+,那么由求积公式(1),得系数:10(1)(1)(1)()!(1)(1)!ni n t t t i t i t n A h dt i n ---+---==--⎰10(1)(1)...(1)(1)...()(0,1,2,...)!(1)!n nb a t t t i t i t n dt i n n i n -----+---=-⎰ (8)则 ()()n i iA b a C =- (9)于是差值求积公式为:()0()()()[]nbn i i ai f x dx b a C f x R f ==-+∑⎰(10)称公式(10)为牛顿-科茨求积公式,其中()n iC 称为科茨系数.显然,科茨系数与被积函数()f x 及积分区间[,]a b 无关,它指依赖于n ,且为多项式积分.因此,只要给出n ,就能看出i A ,并写出相应地牛顿-科茨公式.2牛顿-科茨公式的截断误差与代数精度.当1n =与2n =情况分析牛顿-科茨公式的截断误差为(1)()[]()()()(1)!n b b bn aaaf R f f x dx x dx x dxn ξϕω+=-=+⎰⎰⎰牛顿-科茨公式的截断误差还可以写成(2)*1()[]()((2)!n bn a f R f x dx n n ξω++=+⎰为偶数)(1)*1()[]()(1)!n bn af R f x dx n ξω++=+⎰ (n 为奇数) (11) 其中*[,]a b ξ∈,且不依赖于x ,101()()()...()n n x x x x x x x ω+=---,对()f x 为任何并不超过n 次多项式,均有(1)()0n fx +≡,因而[]0R f ≡,即0()()nbi i ai f x dx A f x ==∑⎰精确成立,也就是说,牛顿-科茨公式的代数精度至少为n ,牛顿-科茨公式在n 为偶数时,至少具有1n +次代数精度,在n 为奇数情况时,至少具有n 次代数精度.2.1.4复化梯形求积公式将区间[,]a b 等分,节点为i x a ih =+ (步长b ah n-=),0,1,2...,i n =)在每个小区间1[,]i i x x -上采用梯形公式(4)得11111()()[(()()]2ii nnbx i i i i ax i i x x f x dx f x dx f x f x ---==-=≈+=∑∑⎰⎰11[()()]2ni i i hf x f x +=+=∑11[()2()()]2n i n i hf a f x f b T -=++=∑ (12)称式(12)为复化梯形公式. 复化梯形公式余项为()2()()()12i n b a R f h f η-''=-(13) 2.1.5复化辛普森求积公式在每个小区间],[1+i i x x 上,辛普森公式(6)得11102()[()4()()]6n bi i ai i hf x dx f x f x f x -++==++∑⎰(14)111012[()4()2((6)]6n n i i i i hf a f x f x f --+===+++∑∑记 )]()(2)(4)([6111021b f x f x f a f hS n i i n i i n +++=∑∑-=-=+ (15)式中,21+i x为],[1+i i x x 的中点,即h x x i i 2121+=+.式(15)称为复化辛普森公式,其余项为∑-=-=-=10)4(4)()2(180)()(n i i n n f h h S f I f R η, 1(,).i i i x x η+∈ 故 ),(),()2(180)(R )4(4b a f h a b f n ∈--=ηη (16) 为复化辛普森的截断误差. 2.1.6复化科茨求积公式将区间[,]a b n 等分, 4n m =,m 为正整数,在每个子区间444[,]k k x x -上用科茨求积公式得到复化求积公式:412()[7()7()32()45mbk ak hf x dx f a f b f x -≈++∑⎰14241411112()32()14()mmm k k k N k k k f xf x f x C ---===+++=∑∑∑ (17)其中 4b a b ah n m--==, k x a kh =+ 其截断误差为6(6)2()[,](),()945n b a R f C h f a b ηη-=-<. 2.1.7 变步长复化求积方法复化求积公式虽然计算简单,也达到了提高精度的目的,但为了满足精度要求必须顾及误差,利用误差公式往往很困难,因为误差表达式中含有未知函数的导数,而估计各阶导数的最大值不太容易.我们可以采取把积分的区间[,]a b 细分的办法,在计算积分时将步长逐步折半,利用前后两次结果进行误差估计,如此继续,直到相邻两次结果相差不大,取最小的步长算出的结果为积分值,这种方法称为变步长积分法.以复化梯形公式为例,把区间[,]a b 分成n 等分,设复化梯形公式的近似值为n T ,原积分值为I ,由复化梯形公式误差公式(14)知:2"11()()()n b a b a I T f a b N N ηη--=-<<再把区间[,]a b 分成2n 等分,得近似值2n T ,则2222()()()122k b a b a I T f a b nηη--''=-<< 假定()f x ''在[,]a b 上变化不大,既有12()()f f ηη''''≈. 由上式得 .24kkI T I T -≈-于是 222211()()341n n n n n n I T T T T T T ≈+-=+-- (18) 式(18)表明若用2n T 作为I 的近似值,其截断误差约为2()3n n T T - (19)2.2 龙贝格求积公式龙贝格积分法的基本思想是采用复化梯形求积方法不断折半步长过程中,在积分结果中加入时候误差估计值进行补偿,使积分计算的收敛性加速,就可以加工出,,,...n n n S C R 精度较高的积分结果.由式(19), 2n T 的误差大致为23n nT T -,因此,可用这个误差值作为2n T 的一种补偿,加到2n T 上,则可得到积分准确值I ,比2n T 的更好近似值~T .222141()333n n n n nT T T T T T =+-=-2221(2)21n n T T =-- (20)式(20)左端1n =时 记122121141()333S T T T T T =+-=- 112()()332a b T b a f +=+- [()4()()]62b a a b f a f f b -+=++恰好为[,]a b 上应用辛普生公式(16)的结果.在每个小区间应用辛普生公式:11[()2()()]2n n k k hT f a f x f b -==++∑121()112[()2()()2()]4n n n k k k k hT f a f x f b f x --===+++∑∑代入式(20)的左端得11111[()2()()2()32n nk k k k h f a f x f b f x -==+++--∑∑ 11[()2()()]2n k k h f a f x f b -++∑11111[()4()2()()]62n n k k k k f a f x f x f b -===+-++∑∑nS =从而复化辛普森公式与复化梯形公式公式有以下关系式2441n nn T T S -=- (21)类似也可以推证,在辛普森序列基础上,利用以下关系式22242161151541n n n n n S S C S S -=-=- (22)可以造出收敛速度更快的科茨序列12,...,...n C C C 将此推行下去,在科茨序列基础上,通过243431n nn C C R -=- (23)构造出收敛速度比科茨序列更快的龙贝格序列12,,......n R R R .以上这种通过逐步构造龙贝格序列的积分近似值法就称为龙贝格积分法.2.3高斯求积公式由定理()()()baf x F b F a =-⎰知,插值型求积公式的代数精度与求积节点的个数有关,具有1n +个节点的插值型求积公式至少具有n 次代数精度.不仅如此,代数精度与节点的选取有关,在构造牛顿-科茨求积公式时,为了简化处理过程,限定用等分节点作为求积节点,这样做,虽然公式确实得到简化,但同时也限制了公式的代数精度. 设积分,1,1=-=b a 本段讨论如下求积公式11()()ni i i f x A f x -==∑⎰(24)对任意积分区间[,]a b ,通过变 22ba t ab x ++-= 可以转换到区间]1,1[-上,这时11()()222bab a b a a bf x dx f t dt ---+=+⎰⎰ 此时,求积公式写为0()()222n bii ai b a a b b af x dx A f t =-+-=+∑⎰若一组节点]1,1[.....,10-∈n x x x 使插值型求积公式(24)具有21n +次代数精度,则称此组节点为高斯点,并称相应求积公式(24)为高斯求积公式.2.3.1 高斯求积公式的余项(2)2()[]()()()(22)!n nbb k k aa k f R f f x dx A f x x dx n ηω+==-=+∑⎰⎰ 其中 01()()()...(),[,]n x x x x x x x ab ωη=---∈,且不依赖于x .2.3.2 复化高斯求积公式复化高斯求积公式的基本思想是:将积分区间[,]a b 分成n个等长小区间1[,](1,...)i i t t i m -=,然后在低阶(2n =)高斯求积公式算出近似值,最后将他们相加的积分()baf t dt ⎰的近似值m G ,即11111111()()[]222ii mmbt i i i i i i at i i t t t t t t f t dt f t dt dt -----==-+-==+∑∑⎰⎰⎰1111[()]222m i h ha i h x dx-==+-+∑⎰101[()]222m n j j mi j h hA f a i h x G ==≈+-+≈∑∑ (25)其中mab h -=,j A 与(0,1,2,...,)j t j n =可由书中表中查出. 3 应用3.1插值型积分的应用例1 用牛顿-科茨公式(1,2,4n =)计算积分12211I x =+⎰. 解 1n =时2210112[]0.4512101()2I -≈+=++2n =时22211112[4]0.463725116101()1()42I -≈++=+++4n =时2222111112[7321232]0.46363311390101()1()1()848I =++++≈++++例2 利用复化梯形求积公式计算积分 12211I dx x =+⎰解 设211)(xx f +=,分点个数为n =1,2,4,5时,求出相应积分n T , 111[(()())],21,2(),.n n i i i i i T f a f b f h b a h n n f x f x a ih ih -=⎧=++⎪⎪-⎪==⎨⎪=⎪⎪=+=⎩∑列表如下:n =1的计算结果见表1-1所列 n h0x 1x 0f1f1T10.50.00.51.0 0.8 0.45n =2的表格如下 n hx1x2xf1f2f2T20.250.00 0.25 0.50 1.00 0.941765 0.80 0.460294n =4时计算结果如下表 n h 0x1x2x3x4x40.1250.00 0.125 0.25 0.375 0.50f1f2f3f4f4T1.00 0.9846154 0.9411765 0.876712 0.80 0.462813n = 5时计算结果如下 n hx1x2x3x4x5x50.10.0 0.1 0.2 0.3 0.40.5f1f2f3f4f5f5T1.0 0.990099 0.9615385 0.91743 0.862069 0.80.463114例3 利用复化求积公式120x e dx ⎰,问积分区间为多少等分才能得证有5位有效数字?解 由式(14)知322()[],()()1212n b a b a R f h f n f n n--''''=-=- 有1(),(),2x x f x e f x e b a ''==-=,当]21,0[∈x 时,在12|()|f x e ''≤,所以122|[]|96n eR f n≤ 由于120x e dx ⎰的准确值具有一位整数,所以要使近似值具有5位有效数字,n 必须满足4242211048,102196⨯≥⨯≤-e n n e 或 取对数有 19=n .即将区间]21,0[19等分可满足给定的精度要求.例4 利用复化抛物线求积公式计算 120211I dx x =+⎰. 解 设11)(2+=x x f ,取m =1,2, 3时,公式()⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎨⎧++=+=====-=+++=+---=-=+∑∑.)12(,2),(),(),(,,242[31221212221111,1222h i a x ih a x x f f x f f b f f a f f m a b n f f f f h S i i i i i i b a m i m i i b a m当m =1,2,3时结果如下表所示 当m =1时m h(0.0)f )25.0(f )5.0(f2S1 0.25 1.0 0.9411765 0.80 0.463725当m =2时mh(0.0)f(0.125)f (0.025)f (0.35)f )5.0(f4S20.125 1.0 0.9846154 0.9411765 0.8767123 0.80 0.463653当m =3时mh(0.0)f(0.08333)f (0.16667)f (0.35)f(0.33333)f (0.14166667)f )5.0(f4S30.83331.00.99310340.9729730.9411760.90.852070.80.4636例5 用复化梯形公式,辛普森公式和科茨公式计算积分10sin xdx x ⎰的近似值.解按精度要求确定]1,0[分多少等分,即确定步长,要使6441021)1(28801|],[|-⨯≤≤M m S f R n ,只需.4642880102M m ⨯≥令10sin ()cos xf x txdt x==⎰,则1()0sin ()()(cos )k kk k k d xd fx tx dt dx x dx==⎰ 1cos().2k t tx kdt π=+⎰dt ktx t x f k k |)2cos(|max )(|max 10)(π+≤⎰11.1k t d t t≤=+⎰)10(≤≤x (4)1max |()| 5.f x ≤所以只要,9.13831288010264=⨯⨯≥-m 取m =4即可, 当4n =时,在每个子区间上用式(25),或(14),或(17),结果.9460829.0,9460833.0,9456911.0888===C S T3.2 龙贝格积分公式应用例6 用龙贝格算法计算积分1241I dx x=+⎰的近似值,要求误差小于510-. 解 .3,0,14)(2==+=b a x x f 步骤如下:2)1(,4)0()1(==f f 得.3)]1()0([211=+=f f T )2(计算,1.3)]21([21,516)21(12=+==f T T f 由此得301333334121=-=T T S . (3)算出),(43),41(f f 从而,3013118)]43()41([412124=++=f f T T,14157.334242=-=T T S .30142121516121=-=S S C(4)计算),87(),85(),83(),81(f f f f 从而得到:13899.3)]87()85()83()81([812148=++++=f f f f T T ,,14159.334482=-=T T S ,14059.31516242=-=S S C.1458.36364121=-=C C R (5)再计算),1615(),1613(),1611(),169(),167(),165(),163(),161(f f f f f f f f 从而得到: 14094.316=T30141598=S ,,14159.3,14159.324==R C 51210||-≤-R R , 所以12043.14159.1dx x ≈+⎰3.3高斯求积公式的应用例7 用两点复化高斯求积公式计算10,x I e dx =⎰要求允许误差.106-=ε解 在本算法中取21=+n 时,,110==A A 其中;,)(mab h e x f x-== =++--=∑=)22(2201j jj b a x a b f A a b G.87189637800.1][21)32121()32121(=++-eem =2时, h =21, ]4121)21([4120202j i j j x i f A G +⨯-=∑∑==.57182571650.1)(41341333413341333413=+++=++--eeee m =3时, h =31. .37182769352.1]631)21([6130203=+⨯-=∑∑==j i j j x i f A G.101027.71||||56323--<⨯≈+-G G G3.4 几种方法的比较分析例8 计算积分211ln 2dx x =⎰,精确到0.001.(1)利用矩形公式计算, 因为对于x x f 1)(=,有320()2f x x''<=<(如果1<x <2),所以按照公式0)2(S =+-dx ba xb a . 0<n R <2112n . 如果取n =10,则我们公式的余项的余数得31010.84101200R -<<⨯,我们还必须加进由于在计算函数值实行四舍五入所产生的误差的界限相差于0.16⨯310-,为了这个目的只要计算1x的值到四位小数精确到0.00005就够了.我们有1232527292132152172192 1.051.151.251.351.551.651.751.851.95x x x x x x x x x =========5128.05405.05714.06061.06897.07407.08.08696.09524.02192172152132927252321=========y y y y y y y y y和6.928469284.0109284.6= (2) 按照梯形公式作同样的计算,在这种情况下,作公式 210,||6n n R R n<<在这儿也试一试取n =10,虽然此时仅可以证3107.16001||-⨯<<n R ,纵坐标是9.18.17.16.15.14.13.12.11.1987654321=========x x x x x x x x x 5263.05556.05882.06250.06667.07143.07692.08333.09091.0987654321=========y y y y y y y y y和1877.669377.01877.621500101=+)( (3) 用辛普森公式做同样的计算作公式 .0))(()2(180)()4(45<≤≤⨯--=n n R b a f n a b R ξξ 并且n =5时有55104.1||-⨯<R .实行计算到五位数字,精确到0.0000058.16.14.12.14321====x x x x 45636.555556.062500.071429.083333.04321和====y y y y 9.17.15.13.11.12927252321=====x x x x x83820.1352632.058824.066667.076923.090909.029********和=====y y y y y.20.150==x x 50000.150000.060000.150和==y y6931525.083820.345636.550000.1301=++)(. 由此可见,用辛普森公式计算得到的值误差最小,计算量相对一般;而用矩形公式计算得到的值误差较大,计算量也比较大;用梯形公式计算的值误差比用矩形公式得到的值要误差小,计算量也是如此.所以我们计算定积分时用辛普森公式往往得到的值误差小,而对没有要求误差大小的,则可以选择辛普森或者是梯形公式,因为这两种方法计算量相对较小.结 束 语本文只讨论了一些一维数值积分方法及其它们的应用,误差分析等有关内容.其中最常用的方法是插值型积分以及复化方法、龙贝格积分方法和高斯积分方法,并讨论了相关求积方法的代数精度和误差分析,并给出了一些例题,分析各种方法的近似值,得出误差分析最小的近似方法.由于篇幅有限,对于高维数值积分方法本文便不再讨论.参考文献[1] 华东师范大学数学系,数学分析(第一版)[M],北京:高等教育出版社,2001. [2] 李庆阳,关治,白峰杉,数值计算原理(第二版)[M],北京: 清华大学出版社, 2008. [3] 肖筱南,现代数值计算方法(第一版)[M],北京: 北京大学出版社, 1999.[4] 菲赫金格尔茨,微积分学教程(第三版)[M],北京: 高等教育出版社, 2005. [5] 裴礼文,数学分析中的典型问题与方法(第一版)[M] ,北京: 北京大学出版社,2004. [6] 李桂成,计算方法(第三版)[M],北京: 高等教育出版社,2010.[7] Yin Y uezhu ,Yang Zhonglian.Calculating Skillfully the Curve Integral and Surface Integral Type 2 bySymmetry, SCIENCE & TECHNOLOGY INFORMATION ,2008(30)The Approximate Numerical Method of the Definite IntegralAbstract This paper mainly discusses common numerical methods of unary function, such as approximate calculation method of interpolation integral, Lebesgue integral and Gauss integration. With these methods in calculating the integral, it will produce some error. In order to reduce the error, we can use after the formula for product and after the Gauss formula. This paper focus on these methods introducing formula of introduction and truncation errors .In addition they can provide examples to analysis size of the error and computation.Keywords interpolation integral Lebesgue integral Gauss integral error analysis approximate computation。
《高数》定积分-课件
性质 6 设M和m分别是函 f(x)数 在[a,b] 上的最大值和则 最小值,
b
m(ba)a f(x)dxM(ba)
又称为定积分理 的估值定
性质( 7 积分中值定理) 设函数f (x)在[a,b]上连续,则
(a,b)内至少存在一 ,点使得
b
a f (x)dx f ()(ba)
例题1 利用定积分的性质,比较下列积分大小
y
yf(x)
x x ax0 xi1 i i
xnb
(2)、近似 在每个小区[x间i1,xi]上任取
一点i,则小曲边梯形的面 Ai积可用以
f (i)为高,以xi 为底的小矩形的面积
f (i)xi 来近似代替,即
Ai f (i)xi (i 1,2,,n)
(3)、求和 把n个小矩形的面积加,
便得曲边梯形面 A的积近似值,即
2、定积分是一种特定 式的 极和 限,它 值仅与被积函 f (x数 )及积分区[a间 ,b]有关 而与积分变量用什 母么 表字 示无关,即
b
b
b
b, 定a当 了 b时,
b
a
a f(x)dxb f(x)dx
4 、 当 ab时,bf规 (x)d定 x0 a
4
(
l
nx)2dx
3
3
例题2 估计下列各积分的值 5π
1) π 4 (1sin2x)dx 解:4在区间[π,5π]上,函数 f (x) 1 sin 2 x 44 之最大值和最小值分别 为
M f (π) 112 2, m f (π) 1 2
积分区间 b a 5πππ 44
5π
π
4 π
定积分的几何意义
对于[a区 , b]上 间的连f(续 x),函 其数 定
第七讲-定积分的近似计算
quad 举例
例:用 quad 计算定积分:
dx 0 1 x 2
1
解:
>> quad('1./(1+x.^2)',0,1)
>> quad('1./(1+x.^2)',0,1,10e-10) 函数表达式一定要用 单引号 括起来! 涉及的运算一定要用 数组运算!
dblquad
i 1
n
通常我们取
x1 x2 xn
h ba n
点 i [ xi 1, xi ] 可以任意选取,常见的取法有: 左端点 xi 1 ,右端点
xi 和中点 ( xi 1 xi ) / 2 。
中点法
左点法
右点法
左点法、右点法和中点法
步长
xi h (b a) / n xi a ih, i 1,2, n
抛物线法
设过以上三点的抛物线方程为: y = x2 + x + = p1(x)
则在区间 [x0, x2] 上,有
x2
x0
f ( x)dx p1 ( x)dx x ( x2 x )dx
x2
0
x2
x0
x x x
3 2
x2
x2 x0 (y0 4y1 y2 ) 6 ba (y0 4 y1 y2 ) 6n
i 1 n
x2 i 2
f ( x )dx
ba ( y2i 2 4 y2i 1 y2i ) i 1 6n
抛物线法
整理后可得:
baຫໍສະໝຸດ b a f ( x)dx [ y0 y2n 4( y1 y3 y2n1 ) 6n 2( y2 y4 y2n2 )]
高中数学(人教版)定积分的近似计算课件
yi1 2
yi
Δxi
,i
1, 2,L
, n.
于是,整个曲边梯形面积的近似值为
b a
f ( x)dx
n i 1
yi
yi 2
1
Δxi
,
即
b a
f ( x)dx
b a ( y0 n2
y1 L
yn1
yn ) , 2
以上近似式称为定积分的梯形法公式.
二、抛物线法
由梯形法求定积分的近似值, 当 y f ( x)为凸曲
P0 , P1,L , P2n ( Pi ( xi , yi ), i 0,1,L , 2n).
现把区间 [ x0 , x2 ] 上的曲线 y f ( x) 用通过三点
P0( x0 , y0 ), P1( x1, y1 ), P2( x2 , y2 )
的抛物线 p1( x) 1x2 1x 1 来近似替代, 便有
x2 f ( x)dx x0
x2 x0
p1( x)dx
x2 x0
(1 x2
1x
1)dx
1
3
( x23
x03 )
1
2
( x22
x03 )
1( x2
x0 )
x2
6
x0 [(1 x02
1 x0
1 ) (1 x22
1 x2
1)
1( x0 x2 )2 21( x0 x2 ) 4 1]
0.5
(1) 用矩形法公式
1 dx
0 1 x2
1 10
(
y0
y1 L
y9 ) 0.8099
(或
1 10
(
y1
y2
L
《定积分的近似计算》课件
03
梯形法与矩形法
梯形法的基本思想
01
梯形法是一种基于几何直观的 定积分近似计算方法。
02
基本思想是将积分区间[a, b]分 成若干个小区间,每个小区间 的长度记为Δx。
03
在每个小区间上,取一个梯形 作为该小区间的近似面积,然 后将所有梯形面积相加,得到 定积分的近似值。
矩形法的基本思想
矩形法也是一种基于几何直 观的定积分近似计算方法。
在积分区间上变化的,与被积函数有关。
02
牛顿-莱布尼兹公式
牛顿-莱布尼兹公式的推导
1 2
背景
牛顿和莱布尼兹时代,微积分的基础尚未完善, 许多概念需要进一步明确。
推导过程
基于无穷小量的思想,通过连续函数的性质和极 限理论,逐步推导出牛顿-莱布尼兹公式。
3
关键步骤
利用定积分的定义,将积分转化为求和的形式, 再利用极限的性质,得到积分的近似值。
、[b, c]和[c, d],有 ∫abf(x)dx+∫bcdf(x)dx=∫adf(x)dx;积分 的可加性是指对于任意两个区间[a, b]和[b, c],有∫abf(x)dx=∫abf(x)dx+∫bcf(x)dx。
定积分的几何意义
总结词
定积分的几何意义是函数曲线与x轴所夹的面积。
详细描述
定积分的几何意义可以通过微元法来理解。微元法是将积分区间分成许多小的区间,每 个小区间上取一个点,做函数在这些点的值的矩形,这个矩形的面积就是微元。所有微 元的面积的和就是定积分的值。这个值也就是函数曲线与x轴所夹的面积。这个面积是
THANKS
感谢观看
梯形法和矩形法的区别 在于对每个小区间上的 近似形状选择不同,梯 形法选择梯形作为近似 形状,而矩形法选择矩 形作为近似形状。
高等数学-第七版-课件-高等数学课件介绍
示例一:导数概念
(1)变速直线运动的速度 (2)平面曲线的切线
匀速运动: v s
t
物 理 问 题 变速运动: v(t0) ?
f (t0)
f (t0 t)
s f (t)
t0
t t0 t
t
s f (t0 t) f (t0 )
s v t
v(t0)
lim
t 0
R
,
2
k)
z
例5 求曲线 x t , y t 2 , z t3在点(1,1,1)处
的切线方程和法平面方程.
o
x
y
切线方程
x x0
(t0 )
y y0
(t0 )
z z0
(t0 )
法平面方程 (t0 )(x x0 ) (t0 ) ( y y0 ) (t0 )(z z0 ) 0
示例二:直线、平面的相互关系
本课件是为教师课堂教学而设计的,不是供学生学习的教案.
设计时,避免让课件“说话”,造成课件与讲授的冲突,而是
给教师讲授留出足够的空间.
这里为不此妨,啰采嗦取几了句许:多方法,比如:将要讲授的道理变成各种 流现 授程在课图某时、些就框课给图件人、常“表常念格把课、要件动讲”画的的;大感课段觉件原。中话其仅放实出在,现课如一件果个里真简。是明这这的样样论,的断在课, 教件师,再那围么绕听这众个多论半断会展不开由讲自解主等地等自.己“念课件”,而不再听 讲。老师的讲课反而影响了听众的“念”。不仅如此,由于 老师另不外知,听随众时念注到意了课哪件里的,播只放顾与自讲己解翻的屏同,步倒.是更加阻碍了 听众。这会导致不折不扣的“冲突”。因此,作者认为: “不让课件说话”是设计课件的一个重要原则
(整理)定积分的近似计算.
实验二定积分的近似计算一、问题背景与实验目的利用牛顿—莱布尼兹公式虽然可以精确地计算定积分的值,但它仅适用于被积函数的原函数能用初等函数表达出来的情形.如果这点办不到或者不容易办到,这就有必要考虑近似计算的方法.在定积分的很多应用问题中,被积函数甚至没有解析表达式,可能只是一条实验记录曲线,或者是一组离散的采样值,这时只能应用近似方法去计算相应的定积分.本实验将主要研究定积分的三种近似计算算法:矩形法、梯形法、抛物线法.对于定积分的近似数值计算,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)例:计算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:function z = integrnd(x, y)z = y*sin(x);7.fprintf(文件地址,格式,写入的变量):把数据写入指定文件.例:x = 0:.1:1;y = [x; exp(x)];fid = fopen('exp.txt','w'); %打开文件fprintf(fid,'%6.2f %12.8f\n',y); %写入fclose(fid) %关闭文件8.syms 变量1 变量2 …:定义变量为符号.9.sym('表达式'):将表达式定义为符号.解释:Matlab中的符号运算事实上是借用了Maple的软件包,所以当在Matlab中要对符号进行运算时,必须先把要用到的变量定义为符号.10.int(f,v,a,b):求f关于v积分,积分区间由a到b.11.subs(f,'x',a):将 a 的值赋给符号表达式 f 中的 x,并计算出值.若简单地使用subs(f),则将f的所有符号变量用可能的数值代入,并计算出值。
高等数学第七版上册总复习PPT课件.ppt
基本初等函数 复合函数
函数的定义
反函数
反函数与直接 函数之间关系
函数 的性质
奇偶性 单调性 有界性 周期性
1
(二)数列极限
2
3
4
(三)函数极限
5
6
7
8
9
10
11
12
13
(三)连续与间断
14
15
第y 一
可去型
类
间
断
点
o x0
x
y
第 二 类 间 断o 点
x0
x
无穷型
x) x)
a0 xn b0 x m
a1 x n1 b1 x m1
an1 x an bm1 x bm
其中m 、n 都是非负整数;a0 , a1 ,, an 及 b0 , b1 ,, bm 都是实数,并且a0 0 ,b0 0 .
真分式化为部分分式之和的待定系数法
40
(2) 三角函数有理式的积分
y 跳跃型
o
x0
x
y
o
x
振荡型
16
17
18
19
第二章 导数与微分
20
21
22
23
24
25
26
第三章 中值定理和导数的应用
27
28
29
30
31
32
33
34
35
第四章 不定积分
36
37
38
39
四 、有理函数与可化为有理函数的积分
(1)有理函数的积分
P( Q(
讨论类型: R( x, n ax b) R( x, n ax b ) cx e
求定积分近似值的常用方
求定积分近似值的常用方
积分是数学中用于解释面积的概念,也是求解不定积分的基本工具。
它分成定积分和不定积分,定积分是一种求积分的方法,它的结果是一个常数,是一种完美的积分方法,可以用来计算更复杂的积分。
但定积分有时不容易计算出最终结果,这时可以采用不定积分来近似求解定积分。
求定积分近似值的常用方法共有四种:梯形法、抛物线法、Simpson 法和 Gauss 法。
梯形法是一种定积分近似值计算的最简单最常用的方法。
它是把积分积分区间划分为若干小区间,每个小区间内都保持函数的连续性,把区间[a,b] 分为 n 个小区间,求出每个小区间的积分函数值之和,即可得到定积分的近似值。
抛物线法是采用抛物线去近似求取定积分的值,它也是一种比较简单的方法,其基本思想是把计算的积分区间分成多个小的子区间,将每个小子区间内的函数曲线用抛物线近似代替,并计算抛物线的积分,最后求得积分和,即可得到定积分的近似值。
Simpson 法是由 Simpson 提出的,是一种使用多项式对函数进行拟合近似求解定积分的方法,其基本思想是将定积分积分区间划分为多个子区
间,使用抛物线来近似求取子区间的积分,最后求得抛物线积分的和,即可得到定积分的近似值。
最后,Gauss 法是一种基于多项式拟合近似求定积分的方法,其基本思想是将定积分积分区间划分为多个子区间,使用多项式来近似求取子区间的积分,最后求得多项式积分的和,从而计算出定积分的近似值。
Gauss 法在计算精确度方面要远超过抛物线法和多项式法。
总之,求定积分近似值的常用方法有梯形法、抛物线法、Simpson 法和 Gauss 法。
这四种方法各有优劣,在使用这些方法时应根据不同的需要来选择合适的方法。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
数学分析 第十章 定积分的应用
高等教育出版社
1
(x x )
1
x0
2 3 ( x2 x0 ) 1 ( x 2 x0 )
*§6 定积分的近似计算
同样地在[ x2 i 2, x2 i ]上用 pi ( x ) i x 2 i x i 替代
二、抛物线法
由梯形法求定积分的近似值, 当 y f ( x ) 为凸曲 线时数值偏大, 为凹曲线时数值偏小. 用抛物线法可克服上述缺点. 将积分区间 [a, b] 作 2n 等分,分点为:
a x0 x1 x2 n ba b, Δxi . 2n
相应的被积函数的值记为 y0 , y1 ,, y2 n ( yi f ( xi ), i 0,1,,2n), 曲线 y f ( x ) 上相应的点记为 P0 , P1 ,, P2 n ( Pi ( xi , yi ), i 0,1,,2n).
b
a
f ( x )dx f ( xi 1 )Δxi .
i 1
在几何意义上,这是用一系列小矩形来近似小曲边
梯形面积的结果, 所以把这个近似计算法称为矩形 法. 矩形法的精度较差,通常使用下面着重介绍的 两种方法.
数学分析 第十章 定积分的应用
高等教育出版社 后退 前进 目录 退出
*§6 定积分的近似计算
0.6
0.7352941
0.7
0.6711409
0.8
0.6097561
0.9
0.5524862
1
0.5
数学分析 第十章 定积分的应用
高等教育出版社
*§6 定积分的近似计算
(1) 用矩形法公式
dx 1 0 1 x 2 10 ( y0 y1 y9 ) 0.8099 1 (或 ( y1 y2 y10 ) 0.7600). 10
曲线 y f ( x ) 上相应的点记为 将曲线上每一段 Pi 1 Pi 用 Pi 1 Pi 替代, 这使每个小
曲边梯形换成了梯形, 其面积为
数学分析 第十章 定积分的应用
高等教育出版社
*§6 定积分的近似计算
y i 1 y i Δxi , i 1, 2, , n. 2
于是,整个曲边梯形面积的近似值为
n
即
b
a
f ( x )dx
i 1
y i y i 1 Δxi , 2
b
a
yn b a y0 f ( x )dx ( y1 yn 1 ) , n 2 2
以上近似式称为定积分的梯形法公式.
数学分析 第十章 定积分的应用
高等教育出版社
*§6 定积分的近似计算
数学分析 第十章 定积分的应用
*§6 定积分的近似计算
利用牛顿 - 莱布尼茨公式虽然可以精确计算 定积分的值,但它仅适合被积函数的原函数能够求 出的情形.我们这里介绍定积分的近似计算方法.
*§6 定积分的近似计算
根据定积分的定义,
或
b
a
f ( x )dx f ( xi )Δxi ,
i 1 n
数学分析 第十章 定积分的应用
高等教育出版社
*§6 定积分的近似计算
现把区间 [ x0 , x2 ] 上的曲线 y f ( x ) 用通过三点 P0 ( x0 , y0 ), P1 ( x1 , y1 ), P2 ( x2 , y2 )
2 p ( x ) x 1 x 1 来近似替代, 便有 的抛物线 1 1
x
x2
0
f ( x )dx
3 2 3 0
x2 x0
p1 ( x )dx (1 x 2 1 x 1 )dx
x2
3 2 x 2 x0 2 2 [(1 x0 1 x0 1 ) (1 x2 1 x2 1 ) 6 1 ( x0 x2 )2 21 ( x0 x2 ) 4 1 ]
1
2( y2 y4 y8 )) 0.7853982.
与精确值
dx π 0 1 x 2 4 0.78539816
1
相比较, 矩形法只有一位有效数字是准确的; 梯形法有三位有效数字是准确的; 而抛物线法有六位有效数字是准确的.
数学分析 第十章 定积分的应用
高等教育出版社
一、梯形法
将积分区间 [a, b] 作 n 等分,分点为
ba a x0 x1 xn b, Δxi . n 相应的被积函数值记为
y0 , y1 ,, yn P0 , P1 ,, Pn ( yi f ( xi ), i 0,1,, n), ( Pi ( xi , yi ), i 0,1,, n).
1
(2) 用梯形法
y10 dx 1 y0 0 1 x 2 10 ( 2 y1 y9 2 ) 0.7850.
1
(3) 用抛物线法
数学分析 第十章 定积分的应用
高等教育出版社
*§6 定积分的近似计算
dx 1 0 1 x 2 30 ( y0 y10 4 y1 y3 y9
曲线 y f ( x ), 将得到
x
x2 i
2i2
f ( x )dx
x2 i x2 i 2
ba pi ( x )dx y2 i 2 4 y2 i 1 y2 i . 6n
n x2 i x2 i 2
最后得到
b
a
f ( x )dx
i 1
f ( x )dx
数学分析 第十章 定积分的应用
高等教育出版社
*§6 定积分的近似计算
dx 例 计算 0 2 的近似值. 1 x
1
解 将区间 [0, 1] 十等分,各分点上被积函数的值列 表如下:
xi yi xi yi
0 1 0.1 0.9900990 0.2 0.9615385 0.3 0.9174312 0.4 0.8620690 0.5 0.8000000
ba y0 y2 n 4( y1 y3 y ) 即 f ( x )dx 2 n 1 a 6n 2( y2 y4 y2 n2 ) . 这就是抛物线公式,亦称为辛普森公式.
b
ba n ( y2 i 2 4 y2 i 1 y2 i ). 6n i 1