鐵嶺市

[拼音]:shuzhi jifen

[外文]:numerical integration

用被積函式的有限個抽樣值的離散和或加權平均值近似地代替定積分的值。在求函式ƒ(x)的定積分

時,常常無法用初等函式表示原函式

,因此能按牛頓-萊布尼茨公式

(1)

計算積分值的定積分是不多的。另外,當ƒ(x)是列表函式時,也不能使用式(1)計算它的積分值。上述事實說明,必須研究近似估算積分的數值積分方法。歷史上,阿基米德、I.牛頓、L.尤拉、C.F.高斯、∏.Л.切比雪夫等人都對此有過貢獻。

數值積分公式

一般是形如

(2)

的近似公式,又稱求積公式,xj和Aj(i=0,1,…,m)分別稱為求積結點和求積係數,通常xj∈[α,b];式(2)右端稱為求積和;兩端之差

稱為求積餘項或求積誤差;區間[α,b]可以是有限的或無限的。構造求積公式的問題就是確定xj和Aj使得 E(ƒ)在某種意義下儘可能地小。

代數精度

若式(2)對ƒ(x)=xk(k=0,1,…,d)精確成立,亦即E(ƒ)=0,而當ƒ(x)=x

時(2)不再是精確等式,則說求積公式(2)的代數精度是d。根據K.外爾斯特拉斯的多項式逼近定理,就一般的連續函式ƒ而言,d越大E(ƒ)越小,因此可以用代數精度的高低說明求積公式的優劣。

插值型求積公式

通過插值途徑構成的求積公式。用ƒ(x)的以x0,x1,…,xm為結點的插值多項式

近似替代ƒ(x)後,經過積分可以得到形如(2)的插值型求積公式,其中求積係數

。 (3)

特別,若所有的xj都屬於[α,b],則稱它為內插型求積公式。這是一類最基本的求積公式。由於m+1個結點的插值型求積公式的代數精度至少是m,所以具有一定代數精度的求積公式總是存在的。

牛頓-科茨公式

等距結點情形下的權函式為1的內插型求積公式。設[α,b]為有限區間,ω(x)呏1。取

,Aj由式(3)確定,則求積公式

(4)

稱為[α,b]上的m+1點牛頓-科茨公式,它的代數精度至少是m。當m=1時,式(4)變成

此式右端等於以ƒ(α)和ƒ(b)為底,以b-α為高的梯形的面積值,故通稱為梯形公式,它的代數精度是1。若ƒ″(x)在[α,b]上連續,則通過積分插值餘項,可知它的求積誤差為

當m=2時,式(4)變成

這是辛普森公式,由於求積結點選得恰當,它的代數精度是3。當ƒ(4)(x)在[α,b]上連續時,它的求積誤差為

當m≥10,牛頓-科茨公式中的求積係數總有一些是負的。這樣的公式在計算上會帶來較大的誤差,一般不被採用。

由上述兩個求積公式的誤差表示式看出,積分割槽間越小,求積誤差就越小。因此為了提高求積精度,可使用復化求積公式。若用分點

將[α,b]n等分,然後對每個子區間[xj,xj+1]應用梯形公式,並對i=0,1,…, n-1求和,即得復化梯形公式

若用分點

將[α,b]2n等分,然後對子區間[x2j,x2j+2]應用辛普森公式,並對i=0,1,…,n-1求和,即得復化辛普森公式

逐次分半演算法和龍貝格公式

遞推關係和逐次分半演算法是數值方法的重要技巧,可用以節省計算時間和計算機的儲存量。龍貝格求積方法正是利用逐次分半演算法和遞推關係構成的一種在現代計算機上十分有效的數值積分法。

下面以梯形公式為例說明逐次分半演算法。在整個區間[α,b]上應用梯形公式算出積分近似值T1;將[α,b]二等分,應用n=2的復化梯形公式算出T2;再將每個小區間二等分(即將[α,b]四等分),應用n=4的復化梯形公式算出T4,如此進行,可得T1,T2,T4,…。在計算T2n時可利用已算出的Tn值:

式中

為復化中矩形公式,這樣,只需要計算ƒ(x)的n個新值即可從Tn得到T2n。顯然,逐次分半演算法充分地利用了前次的計算結果。

比較復化公式S2n、T2n和Tn發現, 適當地組合T2n與Tn可得到代數精度為3的辛普森公式,即有

同樣,適當組合S4n與S2n可得到代數精度為5的求積公式

如此可以引出一系列新公式(遞推關係):

此處,T

呏Tn。上式的代數精度是2k+1。通常稱上式為逐次分半加速公式或龍貝格公式。實際計算可按表1

所示進行:當對角線上相鄰兩個近似值

之差的絕對值小於允許誤差時,計算即可停止,並取

為積分近似值。

高斯型公式

一類具有最高的代數精度的內插型求積公式(表2

)。求積公式(2)含有2(m+1)個自由引數(xj和Aj),恰當選擇這些引數,能使公式(2)的代數精度達到2m+1。高斯求積理論中的一個基本定理斷言:只要把結點x0,x1,…,xm取為區間[α,b]上關於權函式 ω(x)的m+1次正交多項式的零點,內插型求積公式(2)即達到最高代數精度2m+1。這裡[α,b]可以是有限或無限區間,ω(x)為取正值的權函式。

許多有關數值積分的論著都列舉出各種高斯型公式的結點和係數的數值。可以證明:對每個連續函式,當結點個數趨於無窮時,高斯型公式所給出的近似值序列收斂到相應積分的精確值,而牛頓-科茨公式則不具有這種性質。

高維數值積分的主要方法有蒙特卡羅法、代數方法和數論方法。