數(shù)值分析--第4章數(shù)值積分與數(shù)值微分(文件)
【正文】 誤差，要求的高階導(dǎo)數(shù)，由于所以有于是由復(fù)化梯形誤差公式(416)得由復(fù)化辛普森誤差公式(419)得例3 若用復(fù)化求積分公式計算積分的近似值，要求計算結(jié)果有四位有效數(shù)字，應(yīng)取多大？解 因為當(dāng)時，有于是要求計算結(jié)果有四位有效數(shù)字，即要求誤差不超過。故應(yīng)取。注意到每個子區(qū)間經(jīng)過二分只增加了一個分點，用復(fù)化梯形公式求得該子區(qū)間上的積分值為注意，這里代表二分前的步長，將每個子區(qū)間上的積分值相加得即 (420)這表明，將步長由縮小為時，等于的一半再加新增加節(jié)點處的函數(shù)值乘以當(dāng)前步長。令由式(423)，以作為的近似值，其誤差至少為，因此收斂于的速度比快。進(jìn)一步分析，我們有如下歐拉—麥克勞林（EulerMaclaurin）公式： 設(shè)，則有其中系數(shù)與無關(guān)。Romberg公式的計算過程見下表42表4201234(1) 輸入(2) 置(3) 計算對(4) ，輸出停機(jī)；否則，返回(3)。試想如果對節(jié)點不加限制，并適當(dāng)選擇求積系數(shù)，可能會提高求積公式的精度。另一方面，對式(1)，不管如何選擇與，最高精度不可能超過。設(shè)，則，因此，如果是高斯點，則式(1)對于精確成立，即有故式(429)成立。定理表明在上關(guān)于權(quán)的正交多項式系中的次多項式的零點就是求積公式(427)的高斯點。設(shè)在節(jié)點上的次Hermite插值多項式為，即由Hermite余項公式有定理6 高斯求積公式的求積系數(shù)全是正的。 圖1 差商示意圖現(xiàn)在考察式(434)計算近似導(dǎo)數(shù)所產(chǎn)生的截斷誤差，首先分別將在處作Talor展開，有16。定理7 設(shè)，高斯求積公式是收斂的，即6 數(shù)值微分在微積分學(xué)里，求函數(shù)的導(dǎo)數(shù)一般來講是容易辦到的，但若所給函數(shù)由表格給出，則就不那么容易了，這種對列表函數(shù)求導(dǎo)數(shù)通常稱為數(shù)值微分。有了求積節(jié)點后，可如下確定求積系數(shù)其中。對于，用除，記商為，余式為，即，其中，由式(429)可得 (430)由于所給求積公式(427) 是插值型的，它對于是精確成立的，即再注意到，知，從而由(430)有可見求積公式(427)對一切次數(shù)不超過的多項式精確成立，因此為高斯點。定義4 如果求積分公式(427)具有次代數(shù)精度，則稱這組節(jié)點為Gauss點，相應(yīng)公式(427)稱為帶權(quán)的高斯求積公式。設(shè)有個互異節(jié)點的機(jī)械求積分公式 (227)具有次代數(shù)精度，那么有取，式(1)精確成立，即 (228)式(2)構(gòu)成階的非線性方程組，且具有個未知數(shù)，所以當(dāng)給定后
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1