用被積函數的有限個抽樣值的離散和或加權平均值近似地代替定積分的值。在求函數f(x)的定積分
時,常常無法用初等函數表示原函數 ,因此能按牛頓-萊佈尼尼茨公式 (1) 計算積分值的定積分是不多的。另外,當 f( x)是列表函數時,也不能使用式(1)計算它的積分值。上述事實說明,必須研究近似估算積分的數值積分方法。歷史上, 阿基米德、 I.牛頓、 L.歐拉、 C.F.高斯、 ∏.Л.切比雪夫等人都對此有過貢獻。數值積分公式 一般是形如
(2) 的近似公式,又稱求積公式, x j和 A j( i=0,1,…, m)分別稱為求積結點和求積系數,通常 x j∈[ α, b];式(2)右端稱為求積和;兩端之差 稱為求積餘項或求積誤差;區間[ α, b]可以是有限的或無限的。 構造求積公式的問題就是確定 x j和 A j使得 E( f)在某種意義下盡可能地小。代數精度 若式(2)對f(x)=xk(k=0,1,…,d)精確成立,亦即E(f)=0,而當f(x)=x
時(2)不再是精確等式,則說求積公式(2)的代數精度是 d。根據K.外爾斯特拉斯的多項式逼近定理,就一般的連續函數 f而言, d越大 E( f)越小,因此可以用代數精度的高低說明求積公式的優劣。插值型求積公式 通過插值途徑構成的求積公式。用f(x)的以x0,x1,…,xm為結點的插值多項式
, , , 近似替代 f( x)後,經過積分可以得到形如(2)的插值型求積公式,其中求積系數 。 (3) 特別,若所有的 x j都屬於[ α, b],則稱它為內插型求積公式。這是一類最基本的求積公式。由於 m+1個結點的插值型求積公式的代數精度至少是 m,所以具有一定代數精度的求積公式總是存在的。牛頓-科茨公式 等距結點情形下的權函數為1的內插型求積公式。設[α,b]為有限區間,ω(x)≡1。取
, A j由式(3)確定,則求積公式 (4) 稱為[ α, b]上的 m+1點牛頓-科茨公式,它的代數精度至少是 m。當 m=1時,式(4)變成 , 此式右端等於以 f( α)和 f( b)為底,以 b- α為高的梯形的面積值,故通稱為梯形公式,它的代數精度是1。若 f″( x)在[ α, b]上連續,則通過積分插值餘項,可知它的求積誤差為 。 當 m=2時,式(4)變成 , 這是辛普森公式,由於求積結點選得恰當,它的代數精度是3。當 f (4)( x)在[ α, b]上連續時,它的求積誤差為 。 當 m≥10,牛頓-科茨公式中的求積系數總有一些是負的。這樣的公式在計算上會帶來較大的誤差,一般不被采用。由上述兩個求積公式的誤差表達式看出,積分區間越小,求積誤差就越小。因此為瞭提高求積精度,可使用復化求積公式。若用分點
將[ α, b] n等分,然後對每個子區間[ x j, x j +1]應用梯形公式,並對 i=0,1,…, n-1求和,即得復化梯形公式 。 若用分點 將[ α, b]2 n等分,然後對子區間[ x 2 j, x 2 j +2]應用辛普森公式,並對 i=0,1,…, n-1求和,即得復化辛普森公式逐次分半算法和龍貝格公式 遞推關系和逐次分半算法是數值方法的重要技巧,可用以節省計算時間和計算機的存儲量。龍貝格求積方法正是利用逐次分半算法和遞推關系構成的一種在現代計算機上十分有效的數值積分法。
下面以梯形公式為例說明逐次分半算法。在整個區間[α,b]上應用梯形公式算出積分近似值T1;將[α,b]二等分,應用n=2的復化梯形公式算出T2;再將每個小區間二等分(即將[α,b]四等分),應用n=4的復化梯形公式算出T4,如此進行,可得T1,T2,T4,…。在計算T2n時可利用已算出的Tn值:
, 式中 為復化中矩形公式,這樣,隻需要計算 f( x)的 n個新值即可從 T n得到 T 2 n。顯然,逐次分半算法充分地利用瞭前次的計算結果。比較復化公式S2n、T2n和Tn發現,適當地組合T2n與Tn可得到代數精度為3的辛普森公式,即有
表1 龍貝格表 同樣,適當組合 S 4 n與 S 2 n可得到代數精度為5的求積公式 表1 龍貝格表 如此可以引出一系列新公式(遞推關系): 表1 龍貝格表 此處, T 表1 龍貝格表 ≡ T n。上式的代數精度是 2 k+1。通常稱上式為逐次分半加速公式或龍貝格公式。實際計算可按表 1 表1 龍貝格表 所示進行:當對角線上相鄰兩個近似值 和 之差的絕對值小於允許誤差時,計算即可停止,並取 為積分近似值。高斯型公式 一類具有最高的代數精度的內插型求積公式(表2
表2 幾種常用的高斯型求積公式 )。求積公式(2)含有2( m+1)個自由參數( x j和 A j),恰當選擇這些參數,能使公式(2)的代數精度達到2 m+1。高斯求積理論中的一個基本定理斷言:隻要把結點 x 0, x 1,…, x m取為區間[ α, b]上關於權函數 ω( x)的 m+1次正交多項式的零點,內插型求積公式(2)即達到最高代數精度 2 m+1。這裡[ α, b]可以是有限或無限區間, ω( x)為取正值的權函數。許多有關數值積分的論著都列舉出各種高斯型公式的結點和系數的數值。可以證明:對每個連續函數,當結點個數趨於無窮時,高斯型公式所給出的近似值序列收斂到相應積分的精確值,而牛頓-科茨公式則不具有這種性質。
高維數值積分的主要方法有蒙特卡羅法、代數方法和數論方法。