數(shù)值分析牛頓柯特斯求積公式

1、數(shù)值分析數(shù)值分析數(shù)值分析數(shù)值分析第七章第七章 數(shù)值積分與數(shù)值微分?jǐn)?shù)值積分與數(shù)值微分第一節(jié)第一節(jié)等距節(jié)點(diǎn)的等距節(jié)點(diǎn)的Newton-Cotes求積公式求積公式第二節(jié)第二節(jié)復(fù)化求積公式復(fù)化求積公式第三節(jié)第三節(jié)外推算法外推算法第四節(jié)第四節(jié)Gauss型求積公式型求積公式第五節(jié)第五節(jié)數(shù)值微分?jǐn)?shù)值微分?jǐn)?shù)值分析數(shù)值分析數(shù)值分析數(shù)值分析 引引 言言( )( )( )( )( )baNF xf xf x dxF bFewtonLeibn tzai 其其中中為為的的原原函函數(shù)數(shù)公公式式2020,)txedxt 例例如如，對(duì)對(duì)概概率率積積分分 由于被積函數(shù)的原函數(shù)由于被積函數(shù)的原函數(shù)F(x)不可能找到，牛頓不可能找到

2、，牛頓-萊布尼茲公式也就無(wú)能為力了。萊布尼茲公式也就無(wú)能為力了。數(shù)值分析數(shù)值分析數(shù)值分析數(shù)值分析 0 , ()( )()iinbiiaia bxf xf x dxA f x 所所謂謂，從從近近似似計(jì)計(jì)算算的的角角度度看看，就就是是在在區(qū)區(qū)間間上上適適當(dāng)當(dāng)?shù)氐剡x選取取若若干干個(gè)個(gè)點(diǎn)點(diǎn) ，然然后后用用這這些些節(jié)節(jié)點(diǎn)點(diǎn)上上的的函函數(shù)數(shù)值值的的加加權(quán)權(quán)平平均均方方法法獲獲得得定定積積分分的的近近似似值值，即即數(shù)數(shù)值值積積分分( )( )( )( )bbaaxf xf x dxx dx 從從數(shù)數(shù)值值逼逼近近的的觀觀點(diǎn)點(diǎn)看看, ,所所謂謂數(shù)數(shù)值值積積分分，就就是是用用一一個(gè)個(gè)具具有有一一定定精精度度的的簡(jiǎn)

3、簡(jiǎn)單單函函數(shù)數(shù)代代替替被被積積函函數(shù)數(shù)，而而求求出出定定積積分分的的近近似似值值，即即( )( )( ),( )( )nnbbnaaxpxpxf xf x dxpx dx 插插值值型型求求積積公公式式，取?。?）= =得得即即：用用插插值值多多項(xiàng)項(xiàng)式式數(shù)值分析數(shù)值分析數(shù)值分析數(shù)值分析下面推導(dǎo)插值型求積公式下面推導(dǎo)插值型求積公式設(shè)設(shè)x0 ,x1 ,xna,b,pn(x)是是f(x)的的n次次Lagrange插值多項(xiàng)式插值多項(xiàng)式0( )() ( )nniiipxf x l x 則有則有(1)1101( ( )( )( )( )(1)!( )()()(),( )nnnnnfxf xpxwxnwxxx

4、xxxxaxb (1)1( ( )( )( )( )(1)!nbbbnnaaafxf x dxpx dxwx dxn (1)101() ( )( ( )( )(1)!nbbniinaaif x lx dxfxwx dxn (1)101()( ()()(1)!nbniinaiA f xfxwx dxn 數(shù)值分析數(shù)值分析數(shù)值分析數(shù)值分析插值型求積公式插值型求積公式00( )()( )()(1)nnbiiiiaiif x dxA f xR fA f x 其中其中( )0,1,(2)biiaAl x dxin 截?cái)嗾`差或余項(xiàng)為截?cái)嗾`差或余項(xiàng)為(1)11( )( ( )( )(3)(1)!bnnaR f

5、fx wx dxn li(x)為為L(zhǎng)agrange插值基函數(shù)。插值基函數(shù)。(1)101()()( ()()(1)!nbbniinaaif x dxA f xfxwx dxn 數(shù)值分析數(shù)值分析數(shù)值分析數(shù)值分析Ai (i=0,1,n)稱為稱為求積系數(shù)求積系數(shù)，xi (i=0,1,n)稱為稱為求積節(jié)點(diǎn)求積節(jié)點(diǎn)。0( )()nbiiaif x dxA f x 數(shù)數(shù)值值求求積積公公式式的的一一般般形形式式數(shù)值分析數(shù)值分析數(shù)值分析數(shù)值分析一、一、 牛頓牛頓柯特斯求積公式的導(dǎo)出柯特斯求積公式的導(dǎo)出將積分區(qū)間將積分區(qū)間a,b n等分，節(jié)點(diǎn)等分，節(jié)點(diǎn)xi為為xi=a+ih, i=0,1,2,n其中其中h=(b

6、a)/n。有。有第一節(jié)第一節(jié) 等距節(jié)點(diǎn)的牛頓等距節(jié)點(diǎn)的牛頓柯特斯求積公式柯特斯求積公式當(dāng)求積節(jié)點(diǎn)等距分布時(shí)，插值型求積公式稱為當(dāng)求積節(jié)點(diǎn)等距分布時(shí)，插值型求積公式稱為牛頓牛頓柯特斯柯特斯(Newton-Cotes)求積公式。求積公式。0( )()(4)nbiiaif x dxA f x 其中其中( )biiaAl x dx 數(shù)值分析數(shù)值分析數(shù)值分析數(shù)值分析0n()00()(),0,1,bbnjiijijaajinnijjixxAlx dxdxxxtjhdtba Cinij 000011(5)0,1,ninnnn(n)ijjjijitj()Cdt(tj)dtniji!(ni)!nin Ci(n)

7、 稱為柯特斯系數(shù)稱為柯特斯系數(shù)。( )00( )()()()(6)innbniiiaiif x dxA f xbaCf x 于是于是牛頓牛頓柯特斯求積公式為柯特斯求積公式為引進(jìn)變換引進(jìn)變換x=a+th , 0tnxj=a+jh, j=0,1,2,n數(shù)值分析數(shù)值分析數(shù)值分析數(shù)值分析二、兩種特殊的數(shù)值求積公式二、兩種特殊的數(shù)值求積公式: :（1）梯形公式）梯形公式(n=1)x0 =a, x1=b, h= b- a, c0(1)=c1(1) =1/2I=( )( )2babaf(x)dxf af bT 梯形公式的幾何意義梯形公式的幾何意義是用四邊梯形是用四邊梯形x0 ABx1的的面積代替曲邊梯形的面

8、積。面積代替曲邊梯形的面積。xy0ABy=P1(x)y=f(x)f0f1x0=ax1=b圖圖1數(shù)值分析數(shù)值分析數(shù)值分析數(shù)值分析（2）辛卜生公式）辛卜生公式(n=2)辛卜生公式又稱為拋物線公式辛卜生公式又稱為拋物線公式。I=( )4 ()( )62( )4 ()( )32bababaabf(x)dx(f aff b )ShabIf(x)dx(f aff b )S 或或 x0 =a, x1=a+h, x2=b, h= (b-a)/2 C0(2)=1/6,C1(2)=4/6,C2(2)=1/6數(shù)值分析數(shù)值分析數(shù)值分析數(shù)值分析 辛卜生公式的幾何意義是用拋物線辛卜生公式的幾何意義是用拋物線y=P2(x)

9、圍成的圍成的曲邊梯形面積代替由曲邊梯形面積代替由y=f(x)圍成的曲邊梯形面積圖圍成的曲邊梯形面積圖2。)()2(4)(6)(bfbafafabdxxfbaxyx0 x2x1y=P2(x)y=f(x)0圖圖2數(shù)值分析數(shù)值分析數(shù)值分析數(shù)值分析例例：用梯形公式與用梯形公式與辛卜生公式辛卜生公式求求321xIedx 的近似值。的近似值。解：解：辛卜生公式辛卜生公式3123222212(4)0.7665755056xIedxeee I=0.766801031322212()0.8296608192xIedxee 梯形公式梯形公式數(shù)值分析數(shù)值分析數(shù)值分析數(shù)值分析nc0c1c2c3c4c5c6c7c812

10、345678三、牛頓三、牛頓柯特斯系數(shù)柯特斯系數(shù)數(shù)值分析數(shù)值分析數(shù)值分析數(shù)值分析例例n=3為為3/8辛卜生公式辛卜生公式300123( )(33)8xxbaf x dxffff x0 =a, x1=a+h, x2=a+2h, x3=b , h= (b-a)/3n=4為為Cotes公式公式x0 =a, x1=a+h, x2=a+2h, x3=a+3h, x4=b , h= (b-a)/4 430012473212327)90 xxbaf(x)dx(fffff 數(shù)值分析數(shù)值分析數(shù)值分析數(shù)值分析例：例：用用Newton-CotesNewton-Cotes公式計(jì)算公式計(jì)算 解：解：當(dāng)當(dāng)n n取不同值時(shí)

11、，計(jì)算結(jié)果如下所示。取不同值時(shí)，計(jì)算結(jié)果如下所示。 I I準(zhǔn)準(zhǔn)=0.9460831=0.946083110sinxIdxx n 近似結(jié)果近似結(jié)果1 0.92703542 0.94613593 0.94611094 0.94608305 0.9460830數(shù)值分析數(shù)值分析數(shù)值分析數(shù)值分析四、代數(shù)精度四、代數(shù)精度0()()nbiiaifx dxA fx 定義定義1：若求積公式若求積公式 對(duì)一對(duì)一切不高于切不高于m次的多項(xiàng)式次的多項(xiàng)式p(x)都等號(hào)成立，即都等號(hào)成立，即R(p (x)=0; ;而對(duì)于某個(gè)而對(duì)于某個(gè)m+1次多項(xiàng)式等號(hào)不成立，則稱此公式的次多項(xiàng)式等號(hào)不成立，則稱此公式的代數(shù)精度為代數(shù)精度

12、為m.代數(shù)精度代數(shù)精度求法求法 從從(x)=1,x,x2,x3依次驗(yàn)證求積公依次驗(yàn)證求積公式是否成立，若第一個(gè)不成立的等式是式是否成立，若第一個(gè)不成立的等式是xm, ,則其代數(shù)則其代數(shù)精度是精度是m-1. .代數(shù)精度越高，數(shù)值求積公式越精確代數(shù)精度越高，數(shù)值求積公式越精確定義定義2：若求積公式若求積公式 對(duì)對(duì)(x)=1,x,x2,x3xm, 都等號(hào)成立，即都等號(hào)成立，即R(xi)=0; ;而對(duì)于而對(duì)于xm+1 等號(hào)不成立，則稱此公式等號(hào)不成立，則稱此公式 的代數(shù)精度為的代數(shù)精度為m. .0()()nbiiaifx dxA fx 數(shù)值分析數(shù)值分析數(shù)值分析數(shù)值分析定義定義1定義定義221證證義義明

13、明：定定定定義義01110( ),0,1,( ),nbkkkiiainbmmmiiaif xxx dxA xkmf xxxdxA x 已已知知：對(duì)對(duì)有有對(duì)對(duì)有有000( )( )mkkkmmbbbkkkkaaakkf xxf x dxx dxx dx 對(duì)對(duì)，有有00000()mnnmnkkkiiikiiikiikiA xAxA f x 1110( ),nbmmmiiaif xxxdxA x 對(duì)對(duì)有有數(shù)值分析數(shù)值分析數(shù)值分析數(shù)值分析例例1：證明下面數(shù)值求積公式證明下面數(shù)值求積公式具有具有1 1次代數(shù)精度次代數(shù)精度. .101( )( (0)(1)2f x dxff 所以求積公式具有所以求積公式具

14、有1次次代數(shù)精度。代數(shù)精度。10( )11=1( )( (0)(1)12f xf x dxff 取取，左左解解：右右10( )111=( )( (0)(1)222f xxf x dxff 取取，左左右右210( )111=( )( (0)(1)322f xxf x dxff 取取，左左右右數(shù)值分析數(shù)值分析數(shù)值分析數(shù)值分析例例2：設(shè)有設(shè)有成立，確定成立，確定 A0、 A1 、 A2，使上述數(shù)值求積公式的代數(shù)使上述數(shù)值求積公式的代數(shù)精度盡可能高，精度盡可能高，并求代數(shù)精度并求代數(shù)精度。解：解：分別取分別取 (x)=1，x，x2，則有，則有 A0 +A1 + A2=2 -A0 + A2=0 A0 +

15、 A2=2/3解得解得A0 =1/3，A1 =4/3， A2=1/3；111( )( ( 1)4 (0)(1)3f x dxfff 則則取取 (x)=x3，左，左= =右右=0=0； (x)=x4，左左= =-11x4dx=2/5 =2/5 右右=2/3=2/3所以具有所以具有3 3次代數(shù)精度。次代數(shù)精度。10121()( 1)(0)(1)f x dxA fA fA f 數(shù)值分析數(shù)值分析數(shù)值分析數(shù)值分析Newton-Cotes公式的代數(shù)精度公式的代數(shù)精度(1)11(1)!()( )( )bnnnaR ffx dx 其其中中0()()nbjjajf(x)dxA f xR f 因因 為為證證明明：

16、nj 0( )()bijaf x dxA f x 其中其中 n+1(x)= (x-x0)(x-x1).(x-xn-1)(x-xn)即求積公式即求積公式至少具有至少具有n次代次代數(shù)精度。數(shù)精度。定理定理1：由由n+1個(gè)個(gè)互異節(jié)點(diǎn)互異節(jié)點(diǎn)x0 、x1 、x n構(gòu)造的插值構(gòu)造的插值型求積公式的代數(shù)精度至少為型求積公式的代數(shù)精度至少為n。這里系數(shù)這里系數(shù)Aj只依賴于求積節(jié)點(diǎn)與積分區(qū)間只依賴于求積節(jié)點(diǎn)與積分區(qū)間,與與f(x)無(wú)關(guān)。無(wú)關(guān)。顯然當(dāng)顯然當(dāng)f(x)是任何一個(gè)不超過(guò)是任何一個(gè)不超過(guò)n次的多項(xiàng)式時(shí)次的多項(xiàng)式時(shí),余項(xiàng)余項(xiàng)(1)11(1)!( )( )( )0bnnnaR ffx dx 數(shù)值分析數(shù)值分析

17、數(shù)值分析數(shù)值分析 由于由于Newton-Cotes公式是其特殊情形公式是其特殊情形( (等距節(jié)點(diǎn)等距節(jié)點(diǎn)),),它的代數(shù)精度至少是它的代數(shù)精度至少是n,n,還可以證明還可以證明當(dāng)當(dāng)n n 為偶數(shù)時(shí)為偶數(shù)時(shí)Newton-CotesNewton-Cotes公式的代數(shù)精度至少是公式的代數(shù)精度至少是n+1.n+1. 定理定理2：當(dāng)當(dāng)n為偶數(shù)時(shí)為偶數(shù)時(shí)，由由n+1個(gè)等距節(jié)點(diǎn)個(gè)等距節(jié)點(diǎn)x0 、x1 、x n構(gòu)造的牛頓構(gòu)造的牛頓-柯特斯求積公式的代數(shù)精度至少為柯特斯求積公式的代數(shù)精度至少為n+1。數(shù)值分析數(shù)值分析數(shù)值分析數(shù)值分析五、五、( ) , ( ) , ( , )( ) ( )( )( )bbaaf

18、xa bg xa ba bf x g x dxfg x dx （第第二二積積分分中中值值定定理理）如如果果函函數(shù)數(shù)在在上上連連續(xù)續(xù)，函函數(shù)數(shù)在在上上可可積積且且不不變變號(hào)號(hào)，則則存存在在使使引引理理：3( ), ()( )( ( )( )2()( )12bTaf xa bbaRff x dxf af bbaf 設(shè)設(shè)在在 上上有有二二階階連連續(xù)續(xù)導(dǎo)導(dǎo)數(shù)數(shù)，則則梯梯形形求求積積公公式式的的截截?cái)鄶嗾`誤差差定定為為理理3 3：數(shù)值分析數(shù)值分析數(shù)值分析數(shù)值分析帶誤差項(xiàng)的梯形公式是帶誤差項(xiàng)的梯形公式是3( )( )( )212bababaf(x)dxf af bf （）1,(3)( )()()(), ,

19、 2bTanfRfxaxb dxa b 證證明明: :由由截截?cái)鄶嗾`誤差差公公式式有有3( )()()()()( )212bTafbaRfxaxb dxf 證證畢畢( ) , ()()0 , , fxa bxaxbxa ba b 由由于于是是依依賴賴于于 的的函函數(shù)數(shù)且且在在上上連連續(xù)續(xù)，又又，由由引引理理知知，在在區(qū)區(qū)間間上上存存在在一一點(diǎn)點(diǎn) 使使得得數(shù)值分析數(shù)值分析數(shù)值分析數(shù)值分析證：證：已知辛卜生求積公式的代數(shù)精度為已知辛卜生求積公式的代數(shù)精度為3，因此考，因此考慮構(gòu)造一個(gè)三次插值多項(xiàng)式慮構(gòu)造一個(gè)三次插值多項(xiàng)式p3(x)滿足下列條件滿足下列條件根據(jù)插值余項(xiàng)定理得：根據(jù)插值余項(xiàng)定理得：3(

20、 )( )p af a 3( )( )p bf b 322()()ababpf 322()()ababpf (4)3( )24!2( )( )()() ()fa bf xp xx a xx bab 5(4)( ), ()( )( ( )4 ()( )62()2880bSaf xa bbaabRff x dxf aff bbaf 設(shè)設(shè)在在 上上有有4 4階階連連續(xù)續(xù)導(dǎo)導(dǎo)數(shù)數(shù)，則則辛辛卜卜生生求求積積公公式式的的截截?cái)鄶嗾`誤差差為為定定理理4 4：數(shù)值分析數(shù)值分析數(shù)值分析數(shù)值分析得到截?cái)嗾`差得到截?cái)嗾`差3333( )( ) 4()( )( ) 4 ()( )6262bab aabb aabp x

21、dxp app bf aff b (4)21( )( )()() ()4!2baabR ffxa xxb dx 3( )px因因?yàn)闉槭鞘鞘鞘侨未味喽囗?xiàng)項(xiàng)式式，所所以以(4)2321( )( )( )()() ()4!bbba baaaf x dxp x dxfxaxxb dx 兩邊求定積分得兩邊求定積分得 (4)( ),fa b 假設(shè)在區(qū)間上連續(xù)，假設(shè)在區(qū)間上連續(xù)，數(shù)值分析數(shù)值分析數(shù)值分析數(shù)值分析 2(4)2(4)5(4)1()( )()4!21()4!22880babaabR ffxaxxb dxabfxaxxb dxbafab 因此辛卜生求積公式的截?cái)嗾`差為因此辛卜生求積公式的截?cái)嗾`差

22、為 5(4)S(),2880baRffab 2,0,a bxa bxaxxba b 而且當(dāng)時(shí)而且當(dāng)時(shí)由引理知，在上總存在一點(diǎn) 使得由引理知，在上總存在一點(diǎn) 使得數(shù)值分析數(shù)值分析數(shù)值分析數(shù)值分析帶誤差項(xiàng)的辛卜生公式是帶誤差項(xiàng)的辛卜生公式是5(4)()( )4 ()( )( )622880babaabbaf(x)dx(f aff b )fab 023(2)2( )0012(1)( )001,(). , ,( )( )()()(1)()(2)! , ,( )( )()()(1)()12(1)!nnnnnnbnjjajnnnnnbnjjajxa xb hbanfCa bnhff x dxbaftttn

23、 dtCxnfCa bnhff x dxbaft ttn dtCxn 設(shè)設(shè)若若是是偶偶數(shù)數(shù), ,則則若若（ ）（ ）是是奇奇數(shù)數(shù), ,則則數(shù)值分析數(shù)值分析數(shù)值分析數(shù)值分析( )001nnnjjjjAbaC ，例例證證明明：1,( )1,nf x 證證取?。? )00( )()()()innbniiiaiif x dxA f xbaCf x 由由(1)11( )( ( )( )(1)!bnnaR ffx wx dxn 及及( )0R f 知知00( )( )00( )()()()()iinnbiiiaiinnnniiibaf x dxA f xAbaCf xbaC 所所以以( )001nnnjjjjAbaC ，數(shù)值分析數(shù)值分析數(shù)值分析數(shù)值分析 六、六、 初步看來(lái)初步看來(lái)似乎似乎n n值越大，代數(shù)精度越高。是不是值越大，代數(shù)精度越高。是不是 n n 越大越好呢？答案是否定的?？疾煸酱笤胶媚兀看鸢甘欠穸ǖ?。考察Newton-Cotes公式的數(shù)值穩(wěn)定性，即討論舍入誤差對(duì)計(jì)算結(jié)果的公式的數(shù)值穩(wěn)定性，即討論舍入誤差對(duì)計(jì)算結(jié)果的影響。影響。( )00( )()()()innbniiinaiiIf x dxA f xbaCf xI ()kf x n nk