数值分析第四章 数值积分与数值微分ppt课件

1、第四章数值积分数值积分与数值微分与数值微分1 1 引引 言言一、数值积分的必要性一、数值积分的必要性本章主要讨论如下方式的一元函数积分本章主要讨论如下方式的一元函数积分badxxffI)()( )( )( )( )baI ff x dxF bF a在微积分里，按在微积分里，按Newton-Leibniz公式求定积分公式求定积分 f x F x F x F x f x实践问题实践问题例如函数例如函数:2,1,ln1,sin,cos,sin322xexxxxxx2 sinf xx这个问题就是要求由函数这个问题就是要求由函数0 x 48x dxxdxxfL48024802)(cos1)(1 由微积分

2、学我们知道由微积分学我们知道,所求的弧长可表示为所求的弧长可表示为:Whats the Original function?!Its so complex that we can not get it.3222xx)322ln(21693216332412222xxxxxx f x f xx1423454.5688.5原来经过原函数来原来经过原函数来计算积分有它的局计算积分有它的局限性。那限性。那怎样办呢？怎样办呢？呵呵呵呵这就需求积这就需求积分的数值方法来协分的数值方法来协助啦。助啦。二、数值积分的根本思想二、数值积分的根本思想1、定积分的几何意义、定积分的几何意义badxxffI)()(a

3、bxyo f x2、数值积分的实际根据、数值积分的实际根据 f x根据积分中值定理根据积分中值定理, 对于延续函数对于延续函数 ,在在 内存在一点内存在一点 ,使得使得, a b)()()()(fabdxxffIba称称 为区间为区间 的平均高度的平均高度. f, a b ?f3、求积公式的构造、求积公式的构造 假设简单项选择取区间端点或中点的函数值作为平均高假设简单项选择取区间端点或中点的函数值作为平均高度，那么可得一点求积公式如下：度，那么可得一点求积公式如下：左矩形公式：左矩形公式：中矩形公式：中矩形公式：右矩形公式：右矩形公式： Iff aba 2abIffba Iff bbaxyOa

4、b f x f a左矩形公式：左矩形公式： Iff abaxyOab f x2abf2ab中矩形公式：中矩形公式： 2abIffbaxyOab f x f b右矩形公式：右矩形公式： Iff bba 假设取假设取 两点，并令两点，并令 ，那么可得，那么可得梯形公式两点求积公式梯形公式两点求积公式, a b 2f af bf 2f af bIfbaxyOab f x f a f b那么可得那么可得Simpson公式公式(三点求积公式三点求积公式), ,2aba b c 46f af cf bf 假设取三点，假设取三点， 并令并令 46f af cf bIfba 普通地普通地 ，取区间，取区间 内

5、内 个点个点, a b1n ,0,1,2,.,ixin ,0,1,.,if xin f处的高度处的高度经过加权平均的方法近似地得出平均高度经过加权平均的方法近似地得出平均高度这类求积方法称为机械求积这类求积方法称为机械求积:)()()(0ibaniixfabdxxf 或写成或写成: :数值积分公式数值积分公式求积系数求积系数 求积节点求积节点 )()(0kbankkxfAdxxf记记0( )()nnkkkIfA f x0( )( )( )( )(),nbnkkakR fI fIff x dxA f x称为数值称为数值求积公式求积公式称为求积公称为求积公式余项式余项(误误差差).三、求积公式的代

6、数精度三、求积公式的代数精度1、问题的提出、问题的提出构造或确定一个求积公式，要讨论处理的问题有构造或确定一个求积公式，要讨论处理的问题有:(i) 确定求积系数确定求积系数 和求积节点和求积节点 kAkx；(iii)求积公式的误差估计和收敛性分析求积公式的误差估计和收敛性分析.(ii) 断定求积公式精度的衡量规范；断定求积公式精度的衡量规范； 称求积公式称求积公式 具有具有m次代数精度次代数精度,假假设它满足如下两个条件设它满足如下两个条件:2、定义、定义0( )()nnkkkIfA f x(i) 对一切次数对一切次数m次的多项式次的多项式 ,有有)(xPm0)()()(mnmmPIPIPR(

7、ii)存在存在m+1次多项式次多项式 ,使得使得)(1xPm0)()()(111mnmmPIPIPR上述定义中的条件上述定义中的条件(i),(ii)等价于等价于:1( )()0miiR x( )()()()0,(0)kkkniR xI xIxkm2 2 插值型求积公式插值型求积公式一、定义一、定义在积分区间在积分区间 上，上， ,a b取取 个节点个节点1n,0,1,2,.,ix in作作 的的 次代数插值多项式拉格朗日插值公式次代数插值多项式拉格朗日插值公式: f xnnjjjnxfxlxL0)()()(那么有那么有)()()(xRxLxfnn其中，其中，)()!1()()(1)1(xwnf

8、xRnnn为插值余项。为插值余项。于是有：于是有： bajnjbajbanbanbadxxRxfdxxldxxRdxxLdxxf)()()()()()(0取取 babaknkkdxxlxfdxxf)()()(0Ak0()()nbikaikiikxxAdxxx由由 节点节点 决议，决议，与与 无关。无关。 f x称为插值称为插值型求积公型求积公式式二、截断误差与代数精度二、截断误差与代数精度1、截断误差、截断误差0(1)0( )()( )( )()()(1)!nbbkknaaknnbxkakR ff x dxA f xf xLxdxfxxdxn2、代数精度、代数精度0( )nkkkA f x(

9、)bkkaAl x dx推论推论 求积系数求积系数 满足满足:0nkkAba 形如形如 的求积公式至少有的求积公式至少有 n 次代数精度次代数精度 该公式为插值型即：该公式为插值型即： 定理定理kA3 Newton-Cotes3 Newton-Cotes公式公式一、一、Cotes系数系数取节点为等距分布：取节点为等距分布：,0,1,.,ib axa ih hinn 由此构造的插值型求积公式称为由此构造的插值型求积公式称为Newton-Cotes公式公式,此此时求积系数：时求积系数：0()()nxjixj iijxxAdxxx令令htax 00()()( 1)()()!()!n innijijt

10、j hbah dttj dtij hn i niCotes系数系数( ) nkC二、二、Newton-Cotes公式公式1、定义：、定义：记记dtktnjnjCnnjkkjnnj0, 0)()()!( !) 1(那那么么njCabAnjj, 2 , 1 , 0,)()(求积公式变为求积公式变为( )0( )()()nbnjjajf x dxbaCf x称上式为称上式为n阶闭型阶闭型Newton-Cotes求积公式。求积公式。dtktnjnjCnnjkkjnnj0, 0)()()!( !) 1(留意留意:由式由式确定的确定的Cotes系数只与系数只与 和和 有关有关,jn 与与 和积分区间和积分

11、区间 f x, a b无关，无关， 且满足且满足: ()021nnjjC 1nnknkCC2、截断误差、截断误差Newton-Cotes公式的误差为公式的误差为:),(,)()()!1()()!1()()(00)1(21)1(badtjtfnhdxxwnffRnnjnnnban与与x x有关有关3、代数精度、代数精度作为插值型求积公式，作为插值型求积公式，具有具有 次代数精度，次代数精度，n阶阶Newton-Cotes公式至少公式至少n而实践的代数精度能否可以进一步而实践的代数精度能否可以进一步提高呢？提高呢？定理定理当阶数当阶数 为偶数时为偶数时,nNewton-Cotes公式至少具有公式至

12、少具有次代数精度。次代数精度。1n 证明证明:只需验证当只需验证当 为偶数时为偶数时,Newton-Cotes公式对公式对的余项为零。的余项为零。n 1nf xx由于由于 ,所以所以 1nf xx 11 !nfxn即得即得 nnjndtjthfR002)()(引进变换引进变换 ,由于由于 为偶数为偶数,故故 为整数为整数,2ntu2nn于是有于是有 2202)2()(nnnjndujnuhfR据此可断定据此可断定 ,由于上述被积函数是个奇函数由于上述被积函数是个奇函数. 0R f4、数值稳定性、数值稳定性如今讨论舍入误差对计算结果产生的影响如今讨论舍入误差对计算结果产生的影响.设用公式设用公式

13、 njjnjnxfCabfI0)()()()(近似计算积分近似计算积分badxxffI)()(时时, 其中计算函数值其中计算函数值 有误差有误差那么在那么在 的计算中的计算中,由由 引起的误差为引起的误差为jf x0,1,2,.jjn没有误差没有误差, 中间计算过程中的舍入误差也不思索中间计算过程中的舍入误差也不思索, njC计算计算( )nIfj，而，而njjnjnjjjnjnjjnjnCabxfCabxfCabe0)(0)(0)()()()()()(假设假设 都是正数都是正数,并设并设 njC0max|jj n 那么那么有有)(|)(|0)(abCabenjnjn故故 是有界的是有界的,n

14、ejba7n n njC即由即由 引起的误差遭到控制引起的误差遭到控制,的的 倍倍,不超越不超越保证了数值计算的稳定性。保证了数值计算的稳定性。将出现负数将出现负数,而当而当 时时,njnjC0)(|将随将随 增大增大,因此不能保证数值稳定性因此不能保证数值稳定性.故高阶公式不宜采用故高阶公式不宜采用,有适用价值的仅仅是几种低阶的求有适用价值的仅仅是几种低阶的求积公式积公式.三、几种常用的低阶求积公式三、几种常用的低阶求积公式(1)(1)0111,22CCn = 1:( ) ( )( )2bab af x dxf af b梯形公式梯形公式() ()()2!bxafR fxa x b dx/*

15、令令 x = a+th, h = ba, 用中用中值定理值定理 */31( ), , ,121bah fa bh 代数精度代数精度 = 1n = 2:(2)(2)(2)012121,636CCC2( ) ( ) 4 ()( )6ba bab af x dxf aff bSimpson 公式公式代数精度代数精度 = 35(4)1 ( ) ,( , ) ,902baR fh fa bh n = 4: Cotes 公公式式 7(6)8( )945R fh f 01234( )7 ()32 ( ) 12 ()32 ()7 ()90babaf x dxf xf xf xf xf x代数精度代数精度 =

16、5,4kbaxakh h这里这里四、复化求积公式四、复化求积公式 高次插值有高次插值有Runge 景象，怎样办？景象，怎样办？可采用分段低次插值来处理可采用分段低次插值来处理高阶高阶Newton-Cotes公式会出现公式会出现数值不稳定。数值不稳定。而低阶而低阶Newton-Cotes公式公式有时又不能满足精度要求，怎样办？有时又不能满足精度要求，怎样办？可将积分区间可将积分区间 分成假设干小分成假设干小区间，在每个小区间上用区间，在每个小区间上用低阶求积公式计算，然后求和。低阶求积公式计算，然后求和。, a b 复化梯形公式：复化梯形公式：,(0,., )ibahxaihinn在每个在每个

17、上用梯形公式：上用梯形公式：,1iixx111( ) ( )(),0,.,12iixiiiixxxf x dxf xf xin11( )2()( )2nkihf af xf b110( ) ( )()2nbiiaihf x dxf xf x= Tn1321002() ()()1212()( ),( , )12niniiifhhR ffbanhba fa b /*中值定理中值定理*/ 复化梯形公式积分法 复化复化 Simpson 公式：公式：),., 0(,nkhkaxnabhk)()(4)(6)(1211kkkxxxfxfxfhdxxfkkkx21 kx1 kx44444 )()(2)(4)(

18、6)(1010121nknkkkbabfxfxfafhdxxf= Sn)(2180)4(4fhabfR 复化Simpson公式积分法 复化复化 Cotes公式：公式：,(0,., )kbahxak hknn)(7)(32)(12)(32)(790)(14321411kkkkkxxxfxfxfxfxfhdxxfkk101)(7)(32)(12)(32)(790)(432141nkkkkkkbaxfxfxfxfxfhdxxf= Cn),(, )(4945)(2)6(6bafhabfR 收敛速度与误差估计：收敛速度与误差估计：定义：定义：假设一个积分公式的误差满足假设一个积分公式的误差满足 ，0li

19、mphR fCh 且且 ，那么称该公式是，那么称该公式是 p 阶收敛的。阶收敛的。0C 246() ,() ,()nnnTO hSO hCO h例例:利用数据表利用数据表kxkf x01/83/81/25/83/47/811/422.265492.460002.876403.200003.506853.764703.938464计算积分计算积分1*2041Idxx解：解：这个问题有明显的答案这个问题有明显的答案*104arctg |3.1415926.Ix取取n = 8用复化梯形公用复化梯形公式式 1872432852212832412812)0(21818fffffffffT= 3.9884

20、94取取n=4 用辛卜生公式用辛卜生公式 1874432854212834412814)0(61414fffffffffS= 3.141592502运算量根运算量根本一样本一样复化梯形公式的误差估计复化梯形公式的误差估计给定精度给定精度 ，如何取，如何取 ?n例如：要求例如：要求 ，如何判别，如何判别 n = ？ |nTI1、误差先验估计式、误差先验估计式)()(12)()(2fabhfTfIfRn 2max( )a x bMfx 记记那么那么 222)(12)()(12Mabhfabh)()(122 fabhfR ?21()12nkkhfh 22( )( )( )1212bahhfx dxf

21、 bf a 上例中假设要求上例中假设要求 ，那么，那么6| 10nI T226| |(1)(0)|10126nhhR fff0.00244949h 即：取即：取 n = 409通常采取将区间不断对分的方法，即取通常采取将区间不断对分的方法，即取 n = 2k上例中上例中2k 409 k = 9 时，时，T512 = 3.14159202S4 = 3.141592502留意到区间再次对分时留意到区间再次对分时2211 ( )( ) 1224nnhRff bf aR f214nnITIT)(141)(31222nnnnnTTTTTI可用来判别迭代可用来判别迭代能否停顿。能否停顿。2、误差后验估计式

22、、误差后验估计式复化复化Simpson公式的误差估计公式的误差估计1、误差先验估计式、误差先验估计式2、误差后验估计式、误差后验估计式)(2180)4(4fhabSIfRn4(3)(3)1 ( )( )1802nhR fISfbfa )(141)(1512222nnnnnSSSSSI复化复化Cotes公式的误差估计公式的误差估计1、误差先验估计式、误差先验估计式),(, )(4945)(2)6(6bafhabCIfRn2、误差后验估计式、误差后验估计式6)5()5(4)()(9452hafbfCIfRn)(141)(6312322nnnnnCCCCCI四、龙贝格积分四、龙贝格积分例例:计算计算

23、21410 xdx知对于知对于 = 106 须将区间对分须将区间对分 9 次，得到次，得到 T512 = 3.14159202调查调查412 nnTITI由由 来计算来计算 I 效果能否好些？效果能否好些？224414 133nnnnTTITT844133TT= 3.141592502= S4普通有：普通有：nnnSTT 1442nnnCSS 144222nnnRCC 144323Romberg求积求积公式公式 Romberg 算法：算法： ? ? ? T1 =)0(0T T8 =)3(0T T4 =)2(0T T2 =)1(0T S1 =)0(1T R1 =)0(3T S2 =)1(1T C

24、1 =)0(2T C2 =)1(2T S4 =)2(1T 理查德森外推法理查德森外推法利用低阶公式产生高精度的结果。利用低阶公式产生高精度的结果。由由Taylor展开得到：展开得到： i 与与 h 无关无关现将现将 对分，得：对分，得：h0h 0Th 230123.ThIhhh230123.2222hhhhTI设对于某一设对于某一 ， 有公式有公式 近似计算某一未知值近似计算某一未知值 。I如何将公式精度由如何将公式精度由 提高到提高到 ？.432112)()(23322020 hhIhTTh 即：即：230021122( )( )( ).21hTT hT hIhh34212( ).T hIh

25、h211222( )( )21hTT h1212( ).mmmThIhh1122( )( )21mhmmmTTh O h2O h计算步骤：计算步骤：1取 ，计算0hba00( )( )2hTf af b2对对k = 1, 2, 计算以下各步计算以下各步12( )(1)00001121222kkkkkihiTTfah3对对n = 0, 1, 2, k = n 1, n 2, (1)( )( )11441nkkknnnnTTT4收敛控制收敛控制(0)(0)1kkTT(0)(0)1(0)kkkTTT假设假设或或那么输出积分值那么输出积分值 ，否那么转，否那么转3 3。 )0(kTNewton-Cot

26、es公式采用公式采用等距节点作为求积节点代等距节点作为求积节点代数精度至多可到达数精度至多可到达 。 为偶数为偶数1nn那么，在节点个数一定的情那么，在节点个数一定的情况下，能否可以在况下，能否可以在 上自上自在选择节点的位置，使求积在选择节点的位置，使求积公式的精度提得更高公式的精度提得更高 ？, a b例例 ：求形如求形如100111( )()( )f x dxA f xA f x的两点求积公式。的两点求积公式。 1用梯形公式即以用梯形公式即以x0 = -1，x1 = 1为节点的插值型为节点的插值型 求积公式立刻可得求积公式立刻可得 。11( )( 1)(1)f x dxff只具有一次代数

27、只具有一次代数准确度！准确度！2假设对求积公式中的四个待定系数假设对求积公式中的四个待定系数A0, A1, x0, x1适中选取，适中选取，使求积公式对使求积公式对f (x) = 1f (x) = 1，x x，x2x2，x3x3都准确成立，那么都准确成立，那么0101,A A xx需满足如下方程组：需满足如下方程组：0122001 12233001 13344001 1121314AAbaA xAxbaA xAxbaA xAxbaxyo f x11AB0 x1x五、高斯型积分五、高斯型积分0( ) ( )()nbkkakx f x dxA f x构造具有构造具有2n+1次代数精度的求积公式次代

28、数精度的求积公式将节点将节点 以及系数以及系数 都作为待定系数。都作为待定系数。0.nxx0.nAA令令 代入可求解，代入可求解， 2211, ,.,nfxx xx得到的公式得到的公式具有具有 次代数精度。次代数精度。21n节点称为节点称为Gauss 点点此公式称为此公式称为Gauss 型求型求积公式积公式例：求例：求 的的 2 点点 Gauss 公式。公式。dxxfx)(10 解：设解：设 ，应有，应有 3 次代数精度。次代数精度。 101100)()()(xfAxfAdxxfx代入代入 f (x) = 1, x, x2, x3 31130092211200721100521032xAxAx

29、AxAxAxAAA2776. 03891. 02899. 08212. 01010 AAxx不是线性方程组，不是线性方程组，不易求解。不易求解。定理：定理： x0 xn 为为 Gauss 点点 与恣意次数与恣意次数不大于不大于n 的多项式的多项式 P(x) 带权正交。带权正交。nkkxxxw0)()(证明：证明： “ x0 xn 为为 Gauss 点点, 那么公式那么公式 至少有至少有 2n+1 次代数精度。次代数精度。bankkkxfAdxxfx0)()()(对恣意次数不大于对恣意次数不大于n 的多项式的多项式 Pm(x)， Pm(x) w(x)的次的次数不大于数不大于2n+1，那么代入公式

30、应准确成立：，那么代入公式应准确成立：nkkkmkbamxwxPAdxxwxPx0)()()()()(= 00 求求 Gauss 点点 求求w(x)不大于不大于 的多项式的多项式 准确成立，即证明：准确成立，即证明：“要证明要证明 为为 Gauss 点，点， 即要证公式对恣意次数即要证公式对恣意次数设设)()()()(xrxqxwxPmbababamdxxrxdxxqxwxdxxPx)()()()()()()(0nkkkxrA0)(nkkmkxPA0)( nkkmkbamxPAdxxPx0)()()(0,.,nxx21n mPx 正交多项式族正交多项式族 0, 1, , n, 有性质：恣有性质

31、：恣意次数不大于意次数不大于n 的多项式的多项式 P(x) 必与必与n+1 正交。正交。假设取假设取 w(x) 为其中的为其中的n+1，那么，那么n+1的根就是的根就是 Gauss 点。点。53 a0)(10 dxaxx0),(10 j jj j 1021102100)(53(0),(0)(0),(dxcbxxxxdxcbxxxj jj jj jj j215910 cb即：即：22105( )921xxxjStep 1：构造正交多项式：构造正交多项式2设设cbxxxaxxx 2210)(,)(, 1)(j jj jj j再解上例：再解上例： 101100)()()(xfAxfAdxxfxStep 2：求：求2 = 0 的的 2 个根，即为个根，即为 Gauss 点点 x0 ，x1221/20)9/10(9/1021;0 xStep 3：代入：代入 f (x) = 1, x 以求解以求解 A0 ，A1解线性方程组，解线性方程组，简单。简单。结果与前一方法一样：结果与前一方法一样：2776. 0,3891. 0,2899. 0,8212. 01010 AAx