1、精選優質文檔-傾情為你奉上 定積分的近似計算方法 摘要 本文主要討論了一元函數常見的數值積分方法,例如插值型求積公式、龍貝格求積公式、高斯求積公式等近似計算方法,在用這些方法計算定積分時,會產生一些誤差,為了減少誤差, 可以利用復化求積公式、復化高斯公式等.本文圍繞這些方法,系統介紹它們的計算公式以及截斷誤差,并用例題分析它們產生誤差的大小、計算量等.關鍵詞 插值型積分 龍貝格積分 高斯積分 誤差分析 近似計算 1引言 在計算定積分的值時,常常根據微積分學基本定理求出的一個原函數,再用牛頓-萊布尼茨公式求的積分,.但在實際應用中,這種方法只限于解決一小部分定積分的求值問題.當函數沒有具體表達式

2、,只是一些實驗測得數據形成的表格或圖形或者是無法用初等函數表示,例如,等等,這就需要我們用一些近似方法求的積分值.與數值積分一樣,把積分區間細分,在每個小區間上,找到簡單函數來近似代替,且的值容易求的.這樣就把計算復雜的轉化為求簡單的積分值.因此,定積分的近似計算實質上是就是被積函數的近似計算問題.2常見數值方法 2.1牛頓-科茨數值方法牛頓-科茨求積公式是求積節點等距離分布的插值型求積公式.利用插值多項式來構造數值積分公式是最常用、最基本的方法,具體做法是： 給定區間上一組節點,以及節點處函數,作的次拉格朗日多項式 , 其中 ,將插值公式專心-專注-專業 . 其中 ,依賴于變量, 上式積分得

3、 若記 . (1) , (2)則有 (3)稱式(3)為插值求型公式,其中. 與無關,叫求積系數, 為求積節點, 為求積公式余項,其中求積系數由(1)決定. 2.1.1梯形求積公式1梯形公式當插值節點分別選取區間端點時,由式(3)分別求出求積系數 , .從而的求積公式 . (4)稱求積公式(4)為梯形求積公式,簡稱梯形公式. 2梯形公式截斷誤差: . (5) 3梯形求積公式的代數精度:1當時,式(5)中 .精確成立.2.1.2 辛普森求積公式 1辛普森求積公式 當選取節點為時,由式(1)求下列求積系數 , . . 從而求積公式 . （6）稱式(6)為拋物線積分公式或辛普森積分公式.2拋物線求積公

4、式誤差估計 定理1.若在上有四階連續導數,則拋物線公式(6)的余項為: . (7)3拋物線公式的代數精度為3. 易驗證,當時,式(6)精確成立,而當時,式(6)不能精確成立. 2.1.3 牛頓-科茨公式 1牛頓-科茨公式 在等距離節點下,其中. .作為變量替換,那么由求積公式(1),得系數: (8) 則 (9) 于是差值求積公式為: (10)稱公式(10)為牛頓-科茨求積公式,其中稱為科茨系數.顯然,科茨系數與被積函數及積分區間無關,它指依賴于,且為多項式積分.因此,只要給出,就能看出,并寫出相應地牛頓-科茨公式.2牛頓-科茨公式的截斷誤差與代數精度. 當與情況分析牛頓-科茨公式的截斷誤差為

5、牛頓-科茨公式的截斷誤差還可以寫成 為奇數） (11)其中,且不依賴于,對為任何并不超過次多項式,均有,因而,即精確成立,也就是說,牛頓-科茨公式的代數精度至少為,牛頓-科茨公式在為偶數時,至少具有次代數精度,在為奇數情況時,至少具有次代數精度.2.1.4復化梯形求積公式將區間等分,節點為 (步長),)在每個小區間上采用梯形公式(4)得 (12) 稱式(12)為復化梯形公式.復化梯形公式余項為 (13)2.1.5復化辛普森求積公式在每個小區間上,辛普森公式(6)得 （14） 記 (15)式中,為的中點,即.式(15)稱為復化辛普森公式,其余項為 , 故 (16) 為復化辛普森的截斷誤差. 2.

6、1.6復化科茨求積公式 將區間等分, ,為正整數,在每個子區間上用科茨求積公式得到復化求積公式: (17)其中 , 其截斷誤差為. 2.1.7 變步長復化求積方法復化求積公式雖然計算簡單,也達到了提高精度的目的,但為了滿足精度要求必須顧及誤差,利用誤差公式往往很困難,因為誤差表達式中含有未知函數的導數,而估計各階導數的最大值不太容易.我們可以采取把積分的區間細分的辦法,在計算積分時將步長逐步折半,利用前后兩次結果進行誤差估計,如此繼續,直到相鄰兩次結果相差不大,取最小的步長算出的結果為積分值,這種方法稱為變步長積分法.以復化梯形公式為例,把區間分成等分,設復化梯形公式的近似值為,原積分值為,由

7、復化梯形公式誤差公式(14)知: 再把區間分成等分,得近似值,則 假定在上變化不大,既有.由上式得 . 于是 (18)式(18)表明若用作為的近似值,其截斷誤差約為 (19) 2.2 龍貝格求積公式龍貝格積分法的基本思想是采用復化梯形求積方法不斷折半步長過程中,在積分結果中加入時候誤差估計值進行補償,使積分計算的收斂性加速,就可以加工出精度較高的積分結果.由式(19), 的誤差大致為,因此,可用這個誤差值作為的一種補償,加到上,則可得到積分準確值,比的更好近似值. (20)式(20)左端時 記 恰好為上應用辛普生公式(16)的結果.在每個小區間應用辛普生公式: 代入式(20)的左端得 從而復化

8、辛普森公式與復化梯形公式公式有以下關系式 (21)類似也可以推證,在辛普森序列基礎上,利用以下關系式 (22)可以造出收斂速度更快的科茨序列將此推行下去,在科茨序列基礎上,通過 (23)構造出收斂速度比科茨序列更快的龍貝格序列.以上這種通過逐步構造龍貝格序列的積分近似值法就稱為龍貝格積分法.2.3高斯求積公式 由定理知,插值型求積公式的代數精度與求積節點的個數有關,具有個節點的插值型求積公式至少具有次代數精度.不僅如此,代數精度與節點的選取有關,在構造牛頓-科茨求積公式時,為了簡化處理過程,限定用等分節點作為求積節點,這樣做,雖然公式確實得到簡化,但同時也限制了公式的代數精度.設積分本段討論如

9、下求積公式 (24) 對任意積分區間,通過變 可以轉換到區間上,這時 此時,求積公式寫為 若一組節點使插值型求積公式(24)具有次代數精度,則稱此組節點為高斯點,并稱相應求積公式(24)為高斯求積公式. 2.3.1 高斯求積公式的余項 其中 ,且不依賴于. 2.3.2 復化高斯求積公式 復化高斯求積公式的基本思想是:將積分區間分成個等長小區間,然后在低階()高斯求積公式算出近似值,最后將他們相加的積分的近似值,即 (25)其中,與可由書中表中查出.3 應用3.1插值型積分的應用例1 用牛頓-科茨公式()計算積分. 解 時 例2 利用復化梯形求積公式計算積分 解 設,分點個數為1,2,4,5時,

10、求出相應積分, 列表如下:=1的計算結果見表1-1所列10.50.00.51.00.80.45=2的表格如下20.250.000.250.501.000.0.800. =4時計算結果如下表40.1250.000.1250.250.3750.501.000.0.0.0.800.= 5時計算結果如下50.10.00.10.20.30.40.51.00.0.0.917430.0.80. 例3 利用復化求積公式,問積分區間為多少等分才能得證有5位有效數字？ 解 由式（14）知 有,當時,在,所以 由于的準確值具有一位整數,所以要使近似值具有5位有效數字,必須滿足 取對數有 .即將區間19等分可滿足給定

11、的精度要求. 例4 利用復化拋物線求積公式計算 .解 設,取=1,2, 3時,公式當=1,2,3時結果如下表所示當=1時10.251.00.0.800.當=2時20.1251.00.0.0.0.800.當=3時30.83331.00.0.0.0.90.852070.80.4636 例5 用復化梯形公式,辛普森公式和科茨公式計算積分的近似值. 解按精度要求確定分多少等分,即確定步長,要使,只需令,則 所以只要取=4即可, 當時,在每個子區間上用式(25）,或(14),或(17),結果3.2 龍貝格積分公式應用 例6 用龍貝格算法計算積分的近似值,要求誤差小于. 解 步驟如下: 得 計算由此得 .

12、 (3)算出從而 (4)計算從而得到: , (5)再計算 從而得到: , , 所以 3.3高斯求積公式的應用例7 用兩點復化高斯求積公式計算要求允許誤差 解 在本算法中取時,其中 =2時, =, =3時, =. 3.4 幾種方法的比較分析例8 計算積分，精確到0.001.(1)利用矩形公式計算, 因為對于,有(如果1<<2),所以按照公式 . 0<<. 如果取=10,則我們公式的余項的余數得,我們還必須加進由于在計算函數值實行四舍五入所產生的誤差的界限相差于0.16,為了這個目的只要計算的值到四位小數精確到0.00005就夠了.我們有 和6.9284 (2) 按照梯形公

13、式作同樣的計算,在這種情況下,作公式 在這兒也試一試取=10,雖然此時僅可以證,縱坐標是 和 (3) 用辛普森公式做同樣的計算 作公式 并且=5時有.實行計算到五位數字,精確到0. . 由此可見，用辛普森公式計算得到的值誤差最小，計算量相對一般；而用矩形公式計算得到的值誤差較大，計算量也比較大；用梯形公式計算的值誤差比用矩形公式得到的值要誤差小，計算量也是如此.所以我們計算定積分時用辛普森公式往往得到的值誤差小，而對沒有要求誤差大小的，則可以選擇辛普森或者是梯形公式，因為這兩種方法計算量相對較小. 結 束 語 本文只討論了一些一維數值積分方法及其它們的應用,誤差分析等有關內容.其中最常用的方法

14、是插值型積分以及復化方法、龍貝格積分方法和高斯積分方法,并討論了相關求積方法的代數精度和誤差分析,并給出了一些例題,分析各種方法的近似值,得出誤差分析最小的近似方法.由于篇幅有限,對于高維數值積分方法本文便不再討論. 參考文獻 1 華東師范大學數學系，數學分析（第一版）M，北京:高等教育出版社,2001.2 李慶陽,關治,白峰杉，數值計算原理（第二版）M，北京: 清華大學出版社, 2008.3 肖筱南，現代數值計算方法（第一版）M，北京: 北京大學出版社, 1999.4 菲赫金格爾茨，微積分學教程（第三版）M，北京: 高等教育出版社, 2005. 5 裴禮文，數學分析中的典型問題與方法（第一版

15、）M ，北京: 北京大學出版社,2004.6 李桂成，計算方法（第三版）M，北京: 高等教育出版社,2010.7 ,.Calculating Skillfully the Curve Integral and Surface Integral Type 2 by Symmetry, ,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 met