數(shù)值積分的計(jì)算方法論文

1、數(shù)值積分本文應(yīng)用插值積分法和逼近論的思想，簡(jiǎn)單重述了推導(dǎo) 本文應(yīng)用插值積分法和逼近論的思想，簡(jiǎn)單重述了推導(dǎo) Newton-Cotes 公式 公式和 Gauss-Legendre 求積公式的過(guò)程，以及這兩個(gè)公式的系數(shù)、精度等問(wèn)題。并 求積公式的過(guò)程，以及這兩個(gè)公式的系數(shù)、精度等問(wèn)題。并以這兩種數(shù)值積分的求解方法為基礎(chǔ)，應(yīng)用 以這兩種數(shù)值積分的求解方法為基礎(chǔ)，應(yīng)用 quad、guass 函數(shù)編寫(xiě)具體 函數(shù)編寫(xiě)具體 Matlab程序，通過(guò)計(jì)算

2、機(jī)軟件計(jì)算出所給題目的近似數(shù)值積分。對(duì)二者所得的結(jié)果進(jìn) 程序，通過(guò)計(jì)算機(jī)軟件計(jì)算出所給題目的近似數(shù)值積分。對(duì)二者所得的結(jié)果進(jìn)行比較，從而研究了用 行比較，從而研究了用 Newton-Cotes 和 Gauss-Legendre 公式求積分的方法和 公式求積分的方法和二者的精確度問(wèn)題。得知，這兩種求積公式所得的結(jié)果在精度上的確存在差 二者的精確度問(wèn)題。得知，這兩種求積公式所得的結(jié)果在精度上的確存在差異，結(jié)合理論部分更加充分地說(shuō)明了， 異，

3、結(jié)合理論部分更加充分地說(shuō)明了，n 相同時(shí) 相同時(shí) Gauss-Legendre 公式比 公式比Newton-Cotes 公式具有更高的代數(shù)精度，但當(dāng)代數(shù)精度相同時(shí)，二者計(jì)算的結(jié) 公式具有更高的代數(shù)精度，但當(dāng)代數(shù)精度相同時(shí)，二者計(jì)算的結(jié)果仍存在細(xì)微的差異。 果仍存在細(xì)微的差異。理論依據(jù)逼近論——構(gòu)造一個(gè)簡(jiǎn)單函數(shù) p(x)近似表示 f(x),然后對(duì) p(x)求積分得到 f(x)的積分的近似值?；诓逯翟?，推導(dǎo)出數(shù)值積分的基本公式。

4、7;1 插值求積公式為了用數(shù)值方法求 ，對(duì)被積函數(shù) f(x)在給定的 n+1 個(gè)節(jié)ba I(f)= f(x)dx ?點(diǎn)上作 Lagrange 插值，用插值函數(shù) Pn(x)代替 f(x)，就可用 I（Pn(x)）構(gòu)造求積公式，近似地計(jì)算定積分 I(f(x))。§2Newton—Cotes 公式§2.1Newton—Cotes 公式的推導(dǎo)當(dāng)§1.1 插值求積公式的插值節(jié)點(diǎn)為等距節(jié)點(diǎn)時(shí)，就得到 Newton—Co

5、tes公式。將區(qū)間[a,b]n 等分， ，n+1 個(gè)節(jié)點(diǎn)為b a h n? ?xk=a+kh (k=0,1,…,n)在節(jié)點(diǎn)上對(duì) f(x)的 Lagrange 插值多項(xiàng)式是：0 0( ) ( ) ( )n njn kk j k jj kx x p x f x x x ? ? ?? ? ? ? ?用 Pn(x)代替 f(x)構(gòu)造求積公式：（1-3）( 1)0( ) ( ) ( ) ( 1)!n n bk a kf R f x x dx n?

6、 ??? ? ? ? ?討論舍入誤差對(duì)計(jì)算結(jié)果產(chǎn)生的影響，設(shè)（1-2）式近似計(jì)算 ( )ba f x dx ?其中計(jì)算函數(shù)值 f(xn)有誤差值 （k=0,1,2, …,n） 。在（1-2）式中令 εn? 設(shè)計(jì)算 無(wú)誤差，舍入誤差也忽略，則，由（1-2）式計(jì) f(x) ≡ 1,pn(x) = 1n ∑k = 0Cnk Cnk算時(shí) 引式的誤差為 εn( ) ( ) ( ) ( )0 00 0( )[ ( ) ( ( ) ) ( )( ..

7、. )n nn n n nn k k k k n n nk ke b a C f x C f x b a C C ? ? ?? ?? ? ? ? ? ? ? ? ? ? ?如果 皆為正，并設(shè) ，則 ，故 有界，即 Cnk ε = max0 ≤ k ≤ n|εk| |en| ≤ ε(b ? a)n ∑k = 0|Cnk| = ε(b ? a) en引起的誤差受控制，不超過(guò) 倍。保證了數(shù)值計(jì)算的穩(wěn)定性。 εk ( ) b a ? ?但當(dāng) n

8、 8 時(shí)， 將出現(xiàn)負(fù)數(shù)，這時(shí)，數(shù)值計(jì)算的穩(wěn)定性不能保證，所以節(jié) ? Cnk點(diǎn)超過(guò) 8 時(shí) Newton—Cotes 公式不能用。當(dāng) n 為偶數(shù)時(shí)，Newton—Cotes 積分公式具有 n+1 次代數(shù)精度。§2.3 經(jīng)典 Newton—Cotes 公式當(dāng) n=4,5 點(diǎn)公式稱(chēng)為經(jīng)典 Newton—Cotes 公式0 1 2 3 4( ) ( )0 0(7 ( ) 32 ( ) 12 ( ) 32 ( ) 7 ( )) 90(

9、) ( ) ( ( ) 1, ( ) 1 1n nn nk k n k k n kk kb a C f x f x f x f x f xy f x I b a C y R f x p x C? ?? ? ? ? ? ?? ? ? ? ? ? ? ? ?其中(k=0,1,…,4)，它具有 5 次代數(shù)精度。 xk = a + kb ? a4§3 Gauss-Legendre 求積公式在積分區(qū)間[a,b]內(nèi)對(duì)積分節(jié)點(diǎn)不作限制，不取