第四章 數(shù)值積分

1、 計(jì)算方法4 數(shù)值積分0 0、基本概念：求積公式、代數(shù)精度的概念、基本概念：求積公式、代數(shù)精度的概念1、各種求積公式的代數(shù)精度、各種求積公式的代數(shù)精度2、牛頓柯特斯公式、復(fù)、牛頓柯特斯公式、復(fù)合合求積公式、龍貝格公式求積公式、龍貝格公式、高斯型求積公式、高斯型求積公式 主要知識(shí)點(diǎn)主要知識(shí)點(diǎn)2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第2頁(yè)，共45頁(yè)第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法4.1 引言 在一元函數(shù)的積分學(xué)中,我們已經(jīng)熟知,若函數(shù)f(x)在區(qū)間a, b 上連續(xù)且其原?函數(shù)為F(x) ,則可用牛頓萊布尼茲公式 近似計(jì)算( )( )( )(51)baIfx dxF a

2、F b 但是在許多實(shí)際問(wèn)題經(jīng)常遇到下列情況：（1）原函數(shù)存在但不能用初等函數(shù)表示； （2）原函數(shù)可以用初等函數(shù)表示，但結(jié)構(gòu)復(fù)雜：（3）被積函數(shù)沒(méi)有表達(dá)式，僅僅是一張函數(shù)表或圖形。 dxedxxxx10102sin和2222222139( )23( )2323ln( 223)41616 2f xxxF xxxxxx xx Newton-Leibniz公式所不能或很難解決的積分問(wèn)題, 這時(shí)需要用數(shù)值解法來(lái)建立積分的近似計(jì)算方法。., 1 , 0),(,(nixfxii 4.1.1、必要性2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第3頁(yè)，共48頁(yè)2022-3-12第第4 4章章 數(shù)

3、值積分?jǐn)?shù)值積分 計(jì)算方法4.1.2 數(shù)值積分的基本思想 積分值 定義為： 其幾何意義（如圖所示） 由x=a,x=b,y=0以及y=f(x)這四條邊所圍成的曲邊梯形面積。 難于計(jì)算因?yàn)樗幸粭l曲邊y=f(x)。badxxfI)( y=f(x) a b 0110111:( )lim( ),max,nnbkkiiixakiiiii naxxxbf x dxx fxxxxx xx ab曲邊梯形的面積( )yf x 2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第4頁(yè)，共48頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法 建立數(shù)值積分公式的途徑比較多, 其中最常用的有兩種：

4、 （1）由積分中值定理可知，對(duì)于連續(xù)函數(shù)f(x)，在積分區(qū)間a,b內(nèi)存在一點(diǎn)，使得即所求的曲邊梯形的面積恰好等于底為(b-a)，高為 的矩形面積。但是點(diǎn)的具體位置一般是未知的, 因而 的值也是未知的, 稱 為f(x) 在區(qū)間a,b上的平均高度。那么只要對(duì)平均高度 提供一種算法，相應(yīng)地就獲得一種數(shù)值求積方法bafabdxxfba,)()()()(f)(f)(f)(fy=f(x)yab(a+b)/2ab(a+b)/2y=f(x)by(a+b)/2abxaby=f(x)aby=f(x)(1) 1( )()( )( )2baf x dxb af af b梯形公式(2)( )() ()2baabf x

5、dxba f中距公式1( )()( ) 4 ()( )62baa bf x dxb af aff b（）公式2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第5頁(yè)，共48頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法 中矩形公式把a(bǔ),b 的中點(diǎn)處函數(shù)值 作為平均高度f(wàn)()的近似值而獲得的一種數(shù)值積分方法。(矩形近似，還有左距、右距公式) )()(21bfaf)2(baf 梯形公式把f(a), f(b)的加權(quán)平均值 作為平均高度f(wàn)()的近似值而獲得的一種數(shù)值積分方法。(矩形近似) Simpson公式是以函數(shù)f(x)在a, b, (a+b)/2三點(diǎn)的函數(shù)值f(a), f

6、(b), 的加權(quán)平均值 作為平均高度f(wàn)()的近似值而獲得的一種數(shù)值積分方法。（拋物線近似）1()4()( )62abfaffb)2(baf （2）先用某個(gè)簡(jiǎn)單函數(shù) 近似逼近f(x), 用 代替原被積函數(shù)f(x)，即 以此構(gòu)造數(shù)值算法。 從數(shù)值計(jì)算的角度考慮,函數(shù) 應(yīng)對(duì)f(x)有充分的逼近程度,并且容易計(jì)算其積分,因此將 選取為插值多項(xiàng)式, 這樣f(x)的積分就可以用其插值多項(xiàng)式的積分來(lái)近似代替。)(x)(xbabadxxdxxf)()()(x)(x2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第6頁(yè)，共48頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法 上述公式的

7、共同形式是利用上述公式的共同形式是利用f(x)在在a,b 的若干個(gè)點(diǎn)的若干個(gè)點(diǎn)xk上的值進(jìn)行加上的值進(jìn)行加權(quán)平均，其一般形式為：權(quán)平均，其一般形式為：0( )()()nnkkkIfA f x ( )baf x dx 011nnaxxxxb 求積節(jié)點(diǎn)xk求積系數(shù)0 1 , ,kAkn 僅與求積節(jié)點(diǎn)有關(guān)求積公式的截?cái)嗾`差或余項(xiàng)：0()( )()nbnkkakRff x dxA f x (#)(#)形如（形如（*）式的公式稱為）式的公式稱為機(jī)械求積公式機(jī)械求積公式。4.1.3、余項(xiàng)2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第7頁(yè)，共48頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)

8、值積分 計(jì)算方法由定積分定義0001100( )lim( )(1).(2)( )(3)( )nbiiaxiniiiiiiinnniiiiif x dxfxaxxxbsfxxxxSsfx 分割近似求和思路整理2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第8頁(yè)，共48頁(yè)第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法baniiixnxinidxxfxfSxx)()(limlimmax)4(0001求極限001100( )()().()()()()( )bnnaniiiniiiiibafx dxA fxA fxA fxR fA fxR fAA fxfxfx dx由此想到機(jī)械求積公式其中權(quán)系

9、數(shù)，是加權(quán)和，也是的近似值。2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第9頁(yè)，共48頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法數(shù)值積分的思想圖形描述：分割、近似、求和ab取右端點(diǎn)矩形近似復(fù)化型求積公式ab取左端點(diǎn)矩形近似( )y f x -2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第10頁(yè)，共48頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法衡量某種方法好壞的標(biāo)準(zhǔn)： （1）代數(shù)精度 （2）數(shù)值穩(wěn)定性（3）收斂性 對(duì)多少次多項(xiàng)式該方法無(wú)誤差，即計(jì)算值是準(zhǔn)確的。 或者

10、說(shuō)成靈敏度如何，也就是看舍入誤差對(duì)計(jì)算結(jié)果影響的大小。比如病態(tài)方程組，當(dāng)系數(shù)矩陣中的元素有微小變化時(shí)，引起方程組無(wú)解。這實(shí)際上是由舍入誤差或者說(shuō)成舍入誤差的傳遞造成的。 即是截?cái)嗾`差的大小。4.1.4、代數(shù)精度2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第12頁(yè)，共48頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法代數(shù)精度的判別方法求積公式的代數(shù)精度（/*Algebraic Precision */）1Def 如果求積公式對(duì)一切不高于m次的多項(xiàng)式都恒成立，而對(duì)于某個(gè)m+1次多項(xiàng)式不能精確成立，則稱該求積公式具有m次代數(shù)精度。0()()nnkkkIfA f x 0

11、()()nnkkkIfA f x 41 .Th 求積公式具有次m代數(shù)精度的充要條件是 為 時(shí)求積公式精確成立，而 為 時(shí)求積公式不能成為等式。( )f x231mxxxx、 、( )f x1mx 可以驗(yàn)證：可以驗(yàn)證：梯形公式、中矩形公式的代數(shù)精確度均為梯形公式、中矩形公式的代數(shù)精確度均為1 1，SimpsonSimpson公式代公式代數(shù)精確度為數(shù)精確度為3 3。 2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第13頁(yè)，共48頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法代入 f(x)0 = 1： baabdx1112ba =代入 f(x)1 = x ：=代入 f(

12、x)2 = x2 ：222babaxdx 2baab 3323babax dx 222baab 代數(shù)精度代數(shù)精度 = 1。如梯形公式：2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第14頁(yè)，共48頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法 定理定理2 2 對(duì)任意給定的對(duì)任意給定的n+1個(gè)互異節(jié)點(diǎn)：個(gè)互異節(jié)點(diǎn)： bxxxan 10 總存在總存在n+1個(gè)相應(yīng)的求積系數(shù)個(gè)相應(yīng)的求積系數(shù) ，使求積公式（，使求積公式（*）至少具有）至少具有n次次代數(shù)精度代數(shù)精度. kA證明：證明：分別將分別將 代入代入( (* *) )，并令其精確，并令其精確成立，于是得到關(guān)于成立，于

13、是得到關(guān)于 的線性代數(shù)方程組：的線性代數(shù)方程組： ), 2 , 1 , 0()(nixxfi kA, 1 , 0),(11011 nkiibaiikkniabidxxxA（1 1） 其系數(shù)行列式為范德蒙行列式其系數(shù)行列式為范德蒙行列式 njiijnnnnnnnxxxxxxxxxxxxxxV021022221202100)(1111 所以（所以（1 1）有惟一解）有惟一解 。顯然，以此組。顯然，以此組 構(gòu)成的（構(gòu)成的（* *）對(duì)所有次）對(duì)所有次數(shù)數(shù) 的多項(xiàng)式都是精確的的多項(xiàng)式都是精確的. . kAkAn 2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第15頁(yè)，共48頁(yè)2022-3-

14、12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法16利用下面方法求 nkkA0 badxxf)( bandxxL)()()()(0 xlxfxLknkkn nkjjjkjknkxxxxxl0, 2 , 1 , 0,)(（3） nkdxxlAbakk, 2 , 1 , 0,)( （4）定義2 若（*）中的 滿足（3），（4），則稱求積公式（*）為插值型的. kA沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第16頁(yè)，共48頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法例1111( )( 1)4(0)(1).3f x dxfff確定求積公式的代數(shù)精度分析：由等價(jià)定義求代數(shù)精度，只對(duì)最簡(jiǎn)單的函數(shù)xm來(lái)

15、驗(yàn)證。1)1(1111 kdxxIkkk解：0) 1041(31 ;32) 101 (31) 1 ()0(4) 1(31),2()(22Ifffkxxf 時(shí)時(shí)當(dāng)當(dāng) ;0) 101(31) 1 () 0(4) 1(31),3()(33Ifffkxxf 時(shí)時(shí)當(dāng)當(dāng) 。時(shí)時(shí)當(dāng)當(dāng)445232) 101 (31) 1 ()0(4) 1(31),4()(Ifffkxxf 所以該求積公式的代數(shù)精度m=3。),0(1)( kxf時(shí)時(shí)當(dāng)當(dāng) ) 1 ()0(4) 1(31fff 2) 1141 (31 ) 1 () 0(4) 1(31fff ),1()( kxxf時(shí)時(shí)當(dāng)當(dāng) 為為偶偶數(shù)數(shù)為為奇奇數(shù)數(shù)kkk,12, 0

16、;0I ;1I 2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第17頁(yè)，共48頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法例2: 試確定一個(gè)具有3次代數(shù)精度的公式301230( )(0)(1)(2)(3).f x dxA fA fA fA f 解：因?yàn)楫?dāng)30981(0,1,2,3), 3, , 9,24iix dx 時(shí)時(shí)分分別別為為和和故要公式具有3次代數(shù)精度，則必須有0123123123123=3,923=,249=9,81827=.4AAAAAAAAAAAAA 解此方程組得：01343993=, =, =, =.8888AAAA故所求公式為：303993(

17、)(0)(1)(2)(3).8888f x dxffff 2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第18頁(yè)，共48頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法例3 設(shè)積分區(qū)間a, b為0, 2，取 時(shí), 分別用梯形和辛卜生公式 xexxxxxf, 1)(43220)2()0()(ffdxxf20)2() 1 (4)0(31)(fffdxxf計(jì)算其積分結(jié)果并與準(zhǔn)確值進(jìn)行比較解:梯形公式和辛卜生的計(jì)算結(jié)果與準(zhǔn)確值比較如下表所示 f(x) 1 x x2 x3 x4 ex 準(zhǔn)確值準(zhǔn)確值 2 2 2.67 4 6.40 6.389 梯形公式計(jì)算值梯形公式計(jì)算值 2

18、2 4 8 16 8.389 辛卜生公式計(jì)算值辛卜生公式計(jì)算值 2 2 2.67 4 6.67 6.421 從表中可以看出,當(dāng)f(x)是 時(shí),辛卜生公式比梯形公式更精確。 一般說(shuō)來(lái)，代數(shù)精度越高，求積公式越精確。梯形公式和中矩形公式具有1次代數(shù)精度，辛卜生公式有3次代數(shù)精度。下面驗(yàn)證梯形公式。432,xxx2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第19頁(yè)，共48頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法babfafabdxxf)()(2)(取取f(x)f(x)=1時(shí)，時(shí)， abababdxba) 11 (2,1兩端相等兩端相等 取取f(x)=xf(x)=

19、x時(shí)時(shí), , )(21)(2),(212222abbaababxdxba取取f(x)=xf(x)=x2 2 時(shí)時(shí), , baabbabaababdxx)(21)(2),(312222332兩端不相等兩端不相等 所以梯形公式只有所以梯形公式只有1 1次代數(shù)精度。次代數(shù)精度。 兩端相等兩端相等 2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第20頁(yè)，共48頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法例4 試確定一個(gè)至少具有2次代數(shù)精度的公式 40)3()1()0()(CfBfAfdxxf解: 要使公式具有2次代數(shù)精度,則對(duì)f(x)=1,x,x2,求積公式準(zhǔn)確成立，即

20、得如下方程組。 3649834CBCBCBA920,34,94CBA解之得， ) 3(20) 1 (12)0(491)(40fffdxxf所求公式為： 2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第21頁(yè)，共48頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法4.2.0 插值型求積公式 由插值理論可知，任一函數(shù) 給定一組節(jié)點(diǎn) 后，可用一n次多項(xiàng)式 ,對(duì)其插值，即( )f x01naxxxb( )nP x( )( )( )nnfxPxRx因此()()()bbbnnaaafx dxPx dxRx dx當(dāng) 為拉格朗日插值多項(xiàng)式時(shí)，即nP1()()nnkkkPlxfx則(

21、1)1( )( )( )()( )(1)!nnbbbkkaaakffx dxlx fxdxx dxn1( )()nbkknaklx dxfxRf 1()nkknkAfxRf4.2 牛頓 柯特斯(NewtonCotes) 公式其中011011()()()()( )()()()()bbkknkkaakkkkkknxxxxxxxxAl x dxdxxxxxxxxx(1)( ) ( )(1)!bnaf nRfx dxn通常稱公式為插值型求積公式。2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第22頁(yè)，共48頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法為簡(jiǎn)便起見，將積分區(qū)

22、間a,b 劃分為n等分,即 步長(zhǎng)為：求積節(jié)點(diǎn)為：計(jì)算系數(shù)Ak： 由 , 有：nabh), 1 ,0(nkkhaxkhikxxik)( nknnkkkkkkhknkxxxxxxxx)!( !) 1()()()(110作變量代換 ，當(dāng) 時(shí),有,于是可得 thaxkbax,nt, 00000( )(1)(1)(1)(1)()!()!(1)()()!()!nbbikkaaikiiknknnnnknniikxxAlx dxdxxxt ttktktn h hdtknkhbatidtnknk4.2.1 Cotes系數(shù)系數(shù)2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第23頁(yè)，共48頁(yè)2022-

23、3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法dtitknnkCnnkiiknk00)()!(!)1( ( k=0,1,n ) 代入插值求積公式有 nkkkbaxfCabxxf0)()(d)(稱為牛頓-柯特斯求積公式,Ck稱為柯特斯系數(shù)引進(jìn)記號(hào)：kkCabA)( ( k=0,1,n ) 則2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第24頁(yè)，共48頁(yè)容易驗(yàn)證 10nkkC bakkkkdxxlAAabC)(1 0001( )1( )111nnbkkakknbkakbaClx dxbalx dxbadxba2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法顯然, C

24、k是不依賴于積分區(qū)間a,b以及被積函數(shù)f(x)的常數(shù),只要給出n,就可以算出柯特斯系數(shù)。例如：當(dāng)n=1時(shí) 1011002121) 1(! 1! 011tdtCdttC2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第25頁(yè)，共48頁(yè)當(dāng)n=2時(shí) 220012100220(1)1(1)(2 )20 ! 2 !6(1)2(2 )2 1! 1!3(1)1(1)22 ! 0 !6Cttd tCt td tCt td t表-1給出了n從18的柯特斯系數(shù)。 當(dāng)n = 8時(shí)，從表中可以看出出現(xiàn)了負(fù)系數(shù)，從而影響穩(wěn)定性和收斂性，因此實(shí)用的只是低階公式。 2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值

25、積分 計(jì)算方法以此類推得Cotes系數(shù)表:2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第26頁(yè)，共48頁(yè)nkkkbaxfCabxxf0)()(d)(第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法4.2.2 常用的幾個(gè)積分公式 梯形公式(n=1)1101312( ) 11 ()( )( )( ( )( )222() ( )( , )12bTaTCCf x dxT fRfbaT fbaf af bf af bbaRffa b 。且 Simpson公式(n=2)(2)(2)(2)0125(4)141,666( ) ( ( )4 ()( )62() ( )( , )2880bSaSCCCf

26、 x dxS fRfbaabS ff aff bbaRffa b 且2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第27頁(yè)，共48頁(yè)第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法 Newton公式(n=3)(3)(3)(3)(3)01231331,8888( ) ( ( )3 ()3 (2 )( )83bNaCCCCf x dxN fRfbaN ff af ahf ahf bbah，且 其中。 Cotes公式(n=4)(4)(4)(4)(4)(4)0123473212327,9090909090( ) (7 ( )32 ()12 (2 )90 32 (2 )7 ( )4bCaCCCC

27、Cf x dxC fRfbaC ff af ahf ahf ahf bbah，且 其中。2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第28頁(yè)，共48頁(yè)第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法4.2.3 Newton-Cotes公式截?cái)嗾`差及代數(shù)精度公式截?cái)嗾`差及代數(shù)精度3(4)55(4)(4)( )72( ) , (-)-( )12( ) , ,1()()( )( )9022880( ) , ,8()9454SCfxa bbaRfffxa bSimpsonbabaRffffxa bCotesbaRff 定理 若在上連續(xù)，則梯形公式的余項(xiàng)（截?cái)嗾`差）為設(shè)在上連續(xù) 則公式的余項(xiàng)（

28、截?cái)嗾`差）為設(shè)在上連續(xù) 則公式的余項(xiàng)（截?cái)嗾`差）為7(6)(6)3(2)02(1)0()( )( )483840,( )()(1)(2).()(2)!2 ( )(1)(2).()(1)!, , nnnnnnnbafn NewtonCoteshfnttttn dtnnRfhftttn dtnnbaha bn 一般地對(duì)任意公式的截?cái)嗾`差為為偶數(shù)）為奇即數(shù)）其中，幾個(gè)常用的求積公式的截?cái)嗾`差2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第29頁(yè)，共48頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法幾個(gè)常用的求積公式的代數(shù)精度1 T 公式的代數(shù)精度2222233322(

29、)11( )()22 ( ( )( )()( )22( ) 11 ( )()33 ( ( )( )()( )22bbbaaababbbaaabaf xxf x dxxdxxbababaT ff af babf x dxf xxf x dxx dxxbababaT ff af babf x dxT當(dāng)時(shí) 當(dāng)時(shí)所以 公式具有一次的代數(shù)精度ab1次代數(shù)精度用梯形面積近似( )yf x -2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第30頁(yè)，共48頁(yè)第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法2 S公式的代數(shù)精度22222233222331( )( )()21 ( ( )4 ()( )=(

30、4)()626221( )( )()31 ( ( )4 ()( ) =(4()=()62623( )bbaabbaaf xxf x dxxdxbabababaabS ff aff babbaf xxf x dxx dxbabababaabS ff aff babbaf xx當(dāng)時(shí) 當(dāng)時(shí) 當(dāng)3244333441( )()41 ( ( )4 ()( )(4()()62624( ) bbaabaf x dxx dxbabababaabS ff aff babbaf x dxS f時(shí) 即精確成立因此S-公式具有三次代數(shù)精度。同理可得N-公式具有三次代數(shù)精度，C-公式具有五次代數(shù)精度。ab3次代數(shù)精度用拋

31、物形面積近似( )y f x -ab5次代數(shù)精度( )y f x -近似等于曲邊梯形的面積2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第31頁(yè)，共48頁(yè)第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法例4-2-1212111 ln20.693147181 110.75;2 211 111(4)0.6944;36 1221111(1 33)0.69375458233IdxxNewtonLeibnizIdxxISimpsonINewtonI計(jì)算。解：由公式得 由梯形公式 （） 由公式 由公式 0.693175CotesI 由公式得 2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教

32、研室第32頁(yè)，共48頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法例例4-2-2 分別用梯形公式、辛卜生公式和柯特斯公式計(jì)算定積分分別用梯形公式、辛卜生公式和柯特斯公式計(jì)算定積分 的近似值的近似值 (計(jì)算結(jié)果取計(jì)算結(jié)果取5位有效數(shù)字位有效數(shù)字) 15 . 0dxx(1) 用梯形公式計(jì)算用梯形公式計(jì)算 4267767. 0 170711. 025. 0)1 () 5 . 0(25 . 01d15 . 0ffxx(2) 用辛卜生公式用辛卜生公式 10.510.5d0.54(0.51) / 2161 0.707 1140.866 0310.4309340312x x(3) 用柯特斯公

33、式計(jì)算，系數(shù)為用柯特斯公式計(jì)算，系數(shù)為 ,10.51 0.5d70.5 320.625 120.75 320.875 71901 =4.94975 25.29822 10.39223 29.93326 7 0.43096407180 x x 積分的準(zhǔn)確值為積分的準(zhǔn)確值為 43096441. 032d15 . 02315 . 0 xxx三個(gè)求積公式的精度逐漸提高三個(gè)求積公式的精度逐漸提高2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第33頁(yè)，共48頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法例例4-2-3 用辛卜生公式和柯特斯公式計(jì)算定積分用辛卜生公式和柯特斯公式

34、計(jì)算定積分的近似值的近似值,并估計(jì)其誤差并估計(jì)其誤差(計(jì)算結(jié)果取計(jì)算結(jié)果取5位小數(shù)位小數(shù))3321(275)dxxxx解解: (1) 辛卜生公式辛卜生公式 322036225941613)(24)(6 bfbafafabS由于由于 由辛卜生公式余項(xiàng)由辛卜生公式余項(xiàng) 572)(23xxxxf0)()4(xfbafabfR,),(2880)()()4(5知其誤差為知其誤差為 0)(fR(2) 柯特斯公式柯特斯公式 知其誤差為知其誤差為 0)(fR3 17 (1) 32 (1.5) 12 (2) 32 (2.5) 7 (3)901351252 7 3212 9 327 92045883Cfffff

35、2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第34頁(yè)，共48頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法4.3 復(fù)化求積公式/*Compound Quadrature Formula */思想將積分區(qū)間 分成若干個(gè)小區(qū)間，然后在每個(gè)小區(qū)間上采用低階的NewtonCotes公式,a b一、復(fù)化梯形公式：/*Compound Trapezoidal Formula */將積分區(qū)間 n等分： , a b分點(diǎn),kbaxakh hn 在區(qū)間 上采用梯形公式1011 , ,kkx xkn ()( )baI ff x dx 110( )kknxxkf x dx 1102()

36、()()nkknkhf xf xRf 2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第35頁(yè)，共84頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法復(fù)化梯形公式1122( )( )()( )nnkkhTff af xf b 復(fù)化梯形公式的幾何意義小梯形面積之和近似( )y f x -2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第36頁(yè)，共84頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法復(fù)化梯形公式的余項(xiàng)31012()()()nnnkkhRfITff 設(shè)2( ) , f xCa b 101min( )()( )nka x ba x b

37、kmfxfmax fxMn 由介值定理 , a b 101( )()nkkffn 212( )( )( )nnbaRfITfh f 余項(xiàng)估計(jì)式2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第37頁(yè)，共84頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法復(fù)化梯形公式的收斂性1122( )( )()( )nnkkhTff af xf b 10112()()nnkkkkbabaf xf xnn ( )()baf x dx n 01( )lim()nbkkakf x dxfx 其中定積分與區(qū)間分法和 的取法無(wú)關(guān)k 設(shè)( ) , f xC a b 2008年10月6日12時(shí)7

38、分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第38頁(yè)，共84頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法二、復(fù)化Simpson公式： /*Compound Simpon Formula */分點(diǎn),kbaxakh hn 在區(qū)間 上采用Simpson公式10 11, ,kkxxkn ( )( )baI ff x dx 110( )kknxxkf x dx 1110246()()()( )nkknkkhf xf xf xRf 其中122()()kkhf xf x 將積分區(qū)間 n等分： , a b2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第39頁(yè)，共84頁(yè)2022-3-12第第

39、4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法復(fù)化Simpson公式111012426( )( )()()( )nnnkkkkhSff af xf xf b 復(fù)化Simpson公式的幾何意義小拋物面積之和近似( )y f x -2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第40頁(yè)，共84頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法復(fù)化Simpson公式的余項(xiàng)4140180 2( )( )( )()nnnkkhhRfISff 設(shè)4()(),fxCab 144401()()()m in()()()nkaxbaxbkmfxfm ax fxMn 由介值定理 , a b 144

40、01( )( )( )()nkkffn 441802( )( )( )( )nnbahRfISff 余項(xiàng)估計(jì)式2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第41頁(yè)，共84頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法復(fù)化Simpson公式的收斂性111012426( )( )()()( )nnnkkkkhSff af xf xf b ( )()baf x dx n () ,fxC a b 類似地可以得到復(fù)化Cotes公式111100427321290()()()()nnnkkkkhCff afxfx 1130143214()()( )nnkkkkfxfxf b

41、2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第42頁(yè)，共84頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法例2：分別利用復(fù)化梯形公式、復(fù)化Simpson公式計(jì)算積分 的近似值，要求按復(fù)化Simpson公式計(jì)算時(shí)誤差不超過(guò) 。10sin xIdxx 60 510. 解：首先來(lái)確定步長(zhǎng)1bahnn 444418021802( )( )( )nb a hb a hR ffM 復(fù)化Simpson公式的余項(xiàng)：44()max()axbMfx 其中2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第43頁(yè)，共84頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分

42、計(jì)算方法4M本題 的求法：sin( )xf xx 10costxdt 11002( )sincos()fxttxdtttxdt 11220022( )coscos()fxttxdtttxdt 由歸納法知102( )( )cos()kkkfxttxdt 1100121( )( )cos()kkkkfxttxdtt dtk 415M 2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第44頁(yè)，共84頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法4441111180 2900 2()nRfMnn 60 5 10. 4n 解不等式得將區(qū)間 8等分，分別采用復(fù)化Simpson、

43、梯形公式0 1, 0 1/8 1/4 3/8 10.9973980.9896880.976727 1/2 5/8 6/8 7/8 10.9588510.9361560.9088580.8771930.841471ix()if x2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第45頁(yè)，共84頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法復(fù)化梯形公式(n=8)復(fù)化Simpson公式(n=4)81113022 8848153712848( )( )( )( )( ) ( )( )( )( )( )T ffffffffff 18h 0 945692. 411357046

44、4888811321424( )( )( )( )( )( ) ( )( )( )( )Sffffffffff 0 9460832. 14h 0946083070367.2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第46頁(yè)，共84頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法復(fù)化梯形公式1122()( )()( )nnkkhTff af xf b 復(fù)化Simpson公式111012426( )( )()()( )nnnkkkkhSff af xf xf b 2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第47頁(yè)，共84頁(yè)2022-3-12第第4 4

45、章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法（1）使用復(fù)化梯形公式、Simpson公式，首先要確定步長(zhǎng) ；h（2）而步長(zhǎng)要根據(jù)余項(xiàng)確定，這就涉及到高階導(dǎo)數(shù)的估計(jì)；（3）高階導(dǎo)數(shù)的估計(jì)一般比較困難，且估計(jì)值往往偏大；（4）計(jì)算機(jī)上實(shí)現(xiàn)起來(lái)不方便，通常采用“事后估計(jì)法”。三、積分步長(zhǎng)的自動(dòng)選?。鹤⒁馐马?xiàng)：基本思想：將積分區(qū)間逐次分半終止法則：前后兩次近似值的誤差小于已知精度2nnII 2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第48頁(yè)，共84頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法具體過(guò)程（以復(fù)化梯形公式為例）1、首先將區(qū)間 n等分： , a bbahn 1122( )

46、()( )nnkkhTf af xf b 2、再將區(qū)間 2n等分，即步長(zhǎng)減半： , a b12hh 11121102222( )()()( )nnnkkkkhTf af xf xf b 1102122()nnkkbaTf xn 122()()kkhf xf x 2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第49頁(yè)，共84頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法上述條件滿足，程序終止；否則，繼續(xù)分半計(jì)算。3、終止條件：由復(fù)化梯形公式的余項(xiàng)知2112()()nba baITfn 222122()()nba baITfn 24nnITIT ( )fx 變化不大時(shí)

47、由此得到近似關(guān)系式2214 1()nnnITTT 誤差控制條件214 1()nnTT 2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第50頁(yè)，共84頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法收斂速度慢對(duì)于復(fù)化Simpson公式、Cotes公式可以類似得到222141()nnnISSS 223141()nnnICCC 不足2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第51頁(yè)，共84頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法加速收斂2214 1()nnnITTT 222141()nnnISSS 223141()nnnICCC 應(yīng)

48、用步長(zhǎng)逐次減半得到的復(fù)化梯形值、復(fù)化Simpson值、復(fù)化Cotes值與精確值的比較nS nC nR 2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第52頁(yè)，共84頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法4 Romberg積分法/*Romberg Integration Method */ Romberg積分思想由上節(jié)分析知，用復(fù)化梯形公式計(jì)算積分值I2nT的誤差大約為：213()nnTT 令2213()nnnI TTT 243nnTT 由復(fù)化梯形公式知12102122()nnnkkbaTTf xn 2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研

49、室第53頁(yè)，共84頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法12102412333()()nnnnkkTTbaTf xn 1122( )()( )nnkkbaTf af xf bn 11102246( )()( )()nnkkkkbaf af xf bf xn nS 111012426( )( )()()( )nnnkkkkhSff af xf xf b 2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第54頁(yè)，共84頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法梯形加速公式：2244134 14 1nnnnnTTSTT 利用復(fù)化梯形公式前后兩次

50、積分近似值 和 ，按2nTnT照上式作出的線性組合得到了具有更高精度的積分值。上述公式說(shuō)明：Romberg積分公式正是由此產(chǎn)生2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第55頁(yè)，共84頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法Romberg 值序列Simpson加速公式：222441nnSS Cotes加速公式：323441nnnCCR 類似于梯形加速公式的處理方法，得到：222141()nnnISSS 222441nnnSSC 223141()nnnICCC 323441nnCC 2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第56頁(yè)，共8

51、4頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法通過(guò)上述3個(gè)積分值序列求積分近似值的方法，稱之為Romberg積分法。4個(gè)積分值序列： 2kT梯形值序列Simpson值序列Romberg值序列Cotes值序列 2kS 2kC 2kR122244 1kkkTTS 122222441kkkSSC 132232441kkkCCR 2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第57頁(yè)，共84頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法Romberg積分法的一般公式11111441,jm jmjm jjTTT 其中1121,()mmTTm 2222,()

52、mmTSm 3323,()mmTCm 4424,()mmTRm 2 3 4, , ;jmj2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第58頁(yè)，共84頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法Romberg積分表1 1,T2 1,T3 1,T4 1,T5 1,T2 2,T3 2,T4 2,T5 2,T3 3,T4 3,T5 3,T4 4,T5 4,T21111222,()babaTTf 12102122()nnnkkbaTTf xn 2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第59頁(yè)，共84頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積

53、分 計(jì)算方法例3：利用Romberg 積分法式計(jì)算積分 ,要求精確到小數(shù)點(diǎn)后面7位。1 5011.Idxx 解：11( )f xx 111 501 51 052,. ( )( . ).Tff 根據(jù)Romberg 積分法計(jì)算得2 11111 50 750 9535714292,.( .).TTf 2 1112 240 9214285713,.TTT 312 110 750 751 1250 9259835752,. ( .)( .).TTff 3 12 13 240 9167876243,.TTT 3 22 23 3160 9164782283,.TTT 具體結(jié)果見下表2008年10月6日12時(shí)

54、7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第60頁(yè)，共84頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法105.1,kT2,kT3,kT4,kTk132450953571429.0921428571.0925983575.0918741799.0916905342.0916787624.0916327874.0916293190.0916478228.0916297224.0916290077.0916294351.0916290776.5 40 916290762,.T 2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第61頁(yè)，共84頁(yè)2022-3-12第第4 4章章 數(shù)

55、值積分?jǐn)?shù)值積分 計(jì)算方法高精度的求積公式 不失一般性，設(shè)不失一般性，設(shè) ，考慮下列求積公式，考慮下列求積公式 我們將會(huì)看到，適當(dāng)?shù)倪x取求積節(jié)點(diǎn)我們將會(huì)看到，適當(dāng)?shù)倪x取求積節(jié)點(diǎn) 可以使上述求積公式具有可以使上述求積公式具有 次代數(shù)精度，這種高精度的求次代數(shù)精度，這種高精度的求積公式稱為高斯（積公式稱為高斯（Gauss）公式，高斯公式的求積節(jié)點(diǎn)稱為高）公式，高斯公式的求積節(jié)點(diǎn)稱為高斯點(diǎn)。斯點(diǎn)。1,1ab 111nkkkf x dxA f x1,2,kxkn21n高斯點(diǎn)的基本特性 盡管高斯點(diǎn)的確定原則上可以化為代數(shù)問(wèn)題，但是由于所歸結(jié)盡管高斯點(diǎn)的確定原則上可以化為代數(shù)問(wèn)題，但是由于所歸結(jié)的方程組是非

56、線性的，而它的求解存在實(shí)質(zhì)性的困難，所以我的方程組是非線性的，而它的求解存在實(shí)質(zhì)性的困難，所以我們要從研究高斯點(diǎn)的基本特性著手解決高斯公式的構(gòu)造問(wèn)題。們要從研究高斯點(diǎn)的基本特性著手解決高斯公式的構(gòu)造問(wèn)題。 設(shè)設(shè) 是求積公式中的高斯點(diǎn)，令是求積公式中的高斯點(diǎn)，令 則有如下結(jié)論：則有如下結(jié)論： 定理定理 節(jié)點(diǎn)節(jié)點(diǎn) 是高斯點(diǎn)的充分必要條件是高斯點(diǎn)的充分必要條件是多項(xiàng)式是多項(xiàng)式 與一切次數(shù)與一切次數(shù) 的多項(xiàng)式的多項(xiàng)式 正交，即正交，即成立成立(1,2, )kx kn 12nxxxxxxx(1,2, )kx kn x1n P x 110f x P x dx 2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院

57、飛機(jī)設(shè)計(jì)教研室第62頁(yè)，共84頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法63勒讓德多項(xiàng)式 以高斯點(diǎn)以高斯點(diǎn) 為零點(diǎn)的為零點(diǎn)的 次多項(xiàng)式次多項(xiàng)式稱為勒讓德稱為勒讓德(Legendre)多項(xiàng)式。多項(xiàng)式。 一般的，勒讓德多項(xiàng)式可以依據(jù)一般的，勒讓德多項(xiàng)式可以依據(jù)來(lái)求得。來(lái)求得。(1,2, )kx knn 12nnPxxxxxxx 2!12!nnnnndPxxndx2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第63頁(yè)，共84頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法64n 牛頓柯特斯型求積公式是封閉型的（區(qū)間a,b的兩端點(diǎn)a, b均是求積節(jié)點(diǎn)）而

58、且要求求積節(jié)點(diǎn)是等距的,受此限制，牛頓柯特斯型求積公式的代數(shù)精確度只能是n(n為奇數(shù))或n+1(n為偶數(shù))。而如果對(duì)求積節(jié)點(diǎn)也適當(dāng)?shù)倪x取,即在求積公式中不僅Ak而且xk也加以選取,這就可以增加自由度，從而可提高求積公式的代數(shù)精確度。 0 ( )() 22,(0,1, ),(0,1, )nbkkakkf x dxA f xnAknknkk求積公式含有個(gè)待定參數(shù)x適當(dāng)選擇這些參數(shù)使其具有2n+1次代數(shù)精度。這類求積公式稱為高斯公式。x是高斯點(diǎn)。高斯公式高斯公式2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第64頁(yè)，共84頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法0

59、0 ( )()(0,1,)()( )0)()(nbkkakknkkbafx dxA fxxknxxPnxx dxx 定 理 ： 插 值 型 求 積 公 式其 節(jié) 點(diǎn)是 高 斯 點(diǎn) 的 充 分 必 要 條件 是 以 這 些 點(diǎn) 為 零 點(diǎn) 的 多 項(xiàng) 式 與 任 意 次 數(shù) 不 超 過(guò) 的 多 項(xiàng) 式 P(x)均 正 交 ： 高斯公式高斯公式(續(xù)1）2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第65頁(yè)，共84頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法0( )( ) ( )21 ( ) ( )() ()( )0(0,1, ), ( ) ( ) 0nbkkkakkbaP xP xxP xx dxA P xxxknP xx dxk必要性證明： 設(shè)是次數(shù)不超過(guò)n的多項(xiàng)式則次數(shù)不超過(guò) n 。若x 是高斯點(diǎn)，則有 又因故有 21 f(x)=P(x)(x)+Q(x ) ,n充分性證明：對(duì)一次數(shù)不超過(guò)的多項(xiàng)數(shù)f(x),用 （x)除f(x)，則有高斯公式高斯公式(續(xù)2）2008年10月6日12時(shí)7分沈陽(yáng)航空工業(yè)學(xué)院飛機(jī)設(shè)計(jì)教研室第66頁(yè)，共84頁(yè)2022-3-12第第4 4章章 數(shù)值積分?jǐn)?shù)值積分 計(jì)算方法000()(