插值型數(shù)值求積公式

1、1 插值型數(shù)值求積公式 教學(xué)目的 1. 會求插值型數(shù)值求積公式及gauss型數(shù)值求積公式并會討論它們的代數(shù)精度； 2. 理解復(fù)化梯形數(shù)值求積公式及復(fù)化simpson數(shù)值求積公式和余項的推導(dǎo)的基礎(chǔ)上掌握它們； 3. 理解數(shù)值微分公式推導(dǎo)的基礎(chǔ)上掌握一階、二階數(shù)值微分公式及余項； 4. 了解外推原理。教學(xué)重點及難點 重點是插值型數(shù)值求積公式及gauss型數(shù)值求積公式的求解及它們代數(shù)精度的討論；難點是gauss型數(shù)值求積公式節(jié)點的求解方法的推導(dǎo)及求解方法。教學(xué)時數(shù) 12學(xué)時 教學(xué)過程 11一般求積公式及其代數(shù)精度設(shè)是上的權(quán)函數(shù)，是上具有一定光滑度的函數(shù)。用數(shù)值方逑下積分的最一般方法是用在節(jié)點上函數(shù)值

2、的某種線性組合來近似其中是獨立于函數(shù)的常數(shù)，稱為積分系數(shù)，而節(jié)點稱為求積節(jié)點。我們也可將（12）寫成帶余項的形式（12）和（13）都稱之為數(shù)值求積公式或機械求積公式。更一般些的求積公式還可以包含函數(shù)在某些點的低階導(dǎo)數(shù)值。在（1.3）中余項也稱為求積公式的截斷誤差。一個很自然的想法是數(shù)值求積公式要對低次多項式精確成立這就導(dǎo)出了求積公式數(shù)精度的概念。定義1 若求積公式（1.2）對任意不高于次的代數(shù)多項式都精確成立，而對不能精確成立，則稱該求積公式具有次代數(shù)精度。一個求積公式的代數(shù)精度越高，就會對越多的代數(shù)多項式精確成立。例1 確定求積公式的代數(shù)精度。解 。從而該求積公式的代數(shù)精度為。對給定節(jié)點，如

3、何選擇求積系數(shù)使求積公式代數(shù)精度盡可能高，對此可用插值型求積公式來實現(xiàn)。1 2 插值型求積公式對給定求積節(jié)點構(gòu)造求積公式的一種簡單方法是利用插值多項式的準許確積分來作為數(shù)值積分值。設(shè)是關(guān)于的lagrange插值多項式其中為lagrange基函數(shù)。取 其中。定義2 對給定互異求積節(jié)點，若求積系數(shù)是由（1.4）給出的，則稱該求積公式是插值型的。定理1 數(shù)值求積公式（1.2）或（1.3）是插值型的當且僅當它的代數(shù)精度。證明 假設(shè)求積公式（1.2）是插值型的，則上面我們假設(shè)了。從而當為次數(shù)的代數(shù)多項時必精確成立，故有。假設(shè)。注意到多項式的次數(shù)為，對=數(shù)值求積精確成立，從而即其求積系數(shù)由（14）給出。推

4、論1 對給定求積節(jié)點，代精度最高的求積公式是插值型求積公式。例2 求插值型求積公式并確定其代數(shù)精度。解 。從而求積公式為且。對從而。若我們利用hermite插值多項式的準確積分作為數(shù)值積分值，我們可以類似地建立帶有函數(shù)在某些節(jié)點導(dǎo)數(shù)值的插值型求積分式。推論2 若是插值型求積公式，則有余項公式其中13 newton-cotes 求積公式在a,b上的插值型求積公式應(yīng)用最方便、最廣泛，稱之為newton-cotes求積公式。設(shè)令則求積系數(shù)為其中因此，newton-cotes公式為 其中 由（16）給出。求職系數(shù)獨平于區(qū)間a,b 稱之為cotes系數(shù)。cotes系數(shù)可以用（16）計算或查（見表4-1）

5、給出。 n=1,2 的newton-cotes求積是常用公式。n=1的公式稱為梯形公式，其幾何意義是用直邊梯形的面積來近似曲邊梯形面積（圖4-1）。即 表4-1 （1.8）的newton-cotes公式稱為simpson公式: （1.9）simpson公式的幾何意義是用以插值拋物線為曲邊的曲邊梯形面積來近似為曲邊的曲邊梯形面積（如圖42），因此simpson求積公式也稱為拋物線公式。newtoncotes公式分別為simpson法則（公式）和cotes公式。1.4 newtoncotes求積公式的余項定理若，則梯形公式(1.8)的余項為 (1.10)證明由插值型求積公式的余項得利用在上不變的號

6、，由積分中值定理得定理若，則simpson公式(1.9)的余項為 (1.11)證明由例知simpson公式的代數(shù)的精度為。令為的三次hermite插值多項式，滿足插值條件：對多項式，simpson公式精確成立，即：從而利用上小于等于零，由積分中值定理給出可以證明，對一般的，只要充分光滑，newtoncotes公式的余項為 （n 為奇數(shù)）（n為偶數(shù)） （12）其中。例3 用、2、3、4、5 相應(yīng)的newtoncotes公式計算積分解 、2、3、4、5相應(yīng)newtoncotes公式所得積分近似值見表4-2表4-2n積分近似值10.920735420.946145930.946110940.9460

7、83050.9460830積分的準確值是0.9460830。容易發(fā)現(xiàn)的結(jié)果比有顯著改進，但相比較沒有實質(zhì)性的進展。對充分光滑的被積函數(shù)，為了既保證精度又節(jié)約時間，應(yīng)盡量選用n是偶數(shù)的情形。1.5 newtoncotes公式的數(shù)值穩(wěn)定性和收斂性求積分式（1.2）的數(shù)值穩(wěn)定性是指的誤差對數(shù)值積分結(jié)果的影響。若影響很大，就稱該數(shù)值求積公式不穩(wěn)定。設(shè)的近似值。由近似值所得數(shù)值積分值為其誤差為。在的前提下，最大可達一般求積公式對準確成立。因此有對newton-cotes公式來說， 從而當時，是數(shù)值穩(wěn)定的。當時，有正有負，而且有從而高階newton-cotes公式是數(shù)值不穩(wěn)定的。我們可以證明，存在上的連續(xù)

8、函數(shù)，對newton-cotes公式來說，不成立。即newton-cotes公式當時，對連續(xù)函數(shù)的數(shù)值積分不能保證收斂?；谏鲜龇€(wěn)定性、收斂性原因，在職實際計算中，很少采用高階newton-cotes求積人以式，而是采用gauss型求積公式或復(fù)化求積公式來提高數(shù)值積分的精度。2 gauss型求積公式21最高代數(shù)精度求積公式由推論1知，插值型求積公式的代數(shù)精度完全由求積節(jié)點的分布所決定。節(jié)點數(shù)目固定后，節(jié)點分布不同，所達到的確良代數(shù)精度也不同。例4 求節(jié)點使插值型求積公式 （2.1）具有盡可能高的代數(shù)精度。解 首先有由于是插值型的，其代數(shù)精度。令，有，及故只要有，就有。進一步取，有就有。上述方程

9、的解為，對應(yīng)的求積公式為 對于。因此二個節(jié)點的求積公式，代數(shù)精度最高為。 對于任意求積節(jié)點，任意求積系數(shù)，求積公式的代數(shù)精度必小于。這是因為對于有 而 是次代數(shù)多項式，從而。在例4中，這是最高能達到的代數(shù)精度了。下面我們利用正交多項式的根來構(gòu)造代數(shù)精度能達到最高的求積公式。引理1 若是上關(guān)于權(quán)函數(shù)的次正交多項式的根，則插值型求積公式具有代數(shù)精度。證明 設(shè)為任一次數(shù)的代數(shù)多項式，則有其是和為次數(shù)的多項式。于是 其中表示與在上帶權(quán)的內(nèi)積，由于是次正交多項式，次數(shù)小于等于，它們的內(nèi)積為0，而次數(shù)不高于。對于插值型求積公式（2.2）有從而對所有次數(shù)的代數(shù)多項式成立.定義3 個節(jié)點的求積公式(2.2)稱

10、為gauss型求積公式,若其代數(shù)精度達,即達最高.并稱其節(jié)點gauss點.2.2gauss點與正交多項式的聯(lián)系利用正交多項式零點作插值型求積公式,可使其代數(shù)精度達到最高.下面我們給出gauss點與正交多項式零點的聯(lián)系.定理4 求積公式(2.2)是gauss型的,當且僅當gauss點是上關(guān)于權(quán)的次正交多項式的根.證明 充分性即引理1的結(jié)論.下證必要性.置.任取次數(shù)的多項式有用內(nèi)積術(shù)語來描述,即對一切次數(shù)不高于的代數(shù)多項式成立,從而是上關(guān)于權(quán)的次正交多項式. gauss點是次正交多項式的根.2.3gauss求積公式的余項定理5 若,則gauss求積公式(2.2)的余項為 (2.3)證明 取的her

11、mite插值多項式，滿足插值條件由得得利用由積分中值定理即得式（23）。24 gauss求積公式的數(shù)值穩(wěn)定性和收斂性設(shè)為lagrange基函數(shù)。為次代數(shù)理化多項式。于是由知gauss型求積公式是數(shù)值穩(wěn)定的。設(shè)上關(guān)于權(quán)的次正交多項式根為，對應(yīng)的gauss求積公式為引理2 對于有限閉區(qū)間上的任何連續(xù)函數(shù)有證明 上的連續(xù)函數(shù)可以用代數(shù)多基式一致逼近。對任意給定的存在某個多項式，有當時，有從而 上面應(yīng)用了 及 由的任意性得（21）。證畢。定理6 gauss型求積公式是數(shù)值穩(wěn)定的；且對（有限閉區(qū)間上的）連續(xù)函數(shù)，gauss求積的數(shù)值隨節(jié)點數(shù)目的增加而收斂到準確積公值。gauss型求積公式有很多優(yōu)點，但對

12、一般的權(quán)函數(shù)，gauss節(jié)點不容易求。gauss求積系數(shù)多為無理數(shù)，因此不如newton- cotes 求積公式的等距節(jié)點和cotes系數(shù)。當函數(shù)賦值計算量大或計算的積公多，這時gauss型求積公式常被優(yōu)先選取。25 幾個常用的gauss型求積公式常用gauss型求積公式有g(shù)auss-laguerre 求積公式和gauss-laguerre 求積公式等。gauss-laguerre 求積公式：-1，1上關(guān)于權(quán)的gauss型求積公式對應(yīng)的gauss點和求積系數(shù)列在表4-3中表4-3102521364 對于一般區(qū)間上帶權(quán)的gauss權(quán)型求積公式，可通過變量變換，由gauss-legendre求積公

13、式得到：例5 用二點、三點gauss型求積公式計算 解 令用二節(jié)點、三節(jié)點計算結(jié)果列在表4-4中。表4-4與newton-cotes公式相比較，近似值要精確得多。gauss-chebyshev求積公式： gauss-laguerre求積公式：表4.5 1116.28994508290.010389256520.58578643763.4142356240.85355339060.146446609440.32254768961.74576110124.53662029699.39507091230.60315410430.35741869240.03888790850.000539294730

14、.41577455682.29428036030.71109300990.2785177336gauss-hermite 求積公式：上關(guān)于權(quán)函數(shù)的gauss型求積公式。對應(yīng)的gauss點和求積系數(shù)列在表4-6中。表4-6 101.772453850952.02018287050.958572464600.019953242060.39361932320.945308720520.7071678120.886226925531.224744871400.29540897521.181635900662.35060497371.33584907400.43607741190.00453000990

15、724629595241.65068012390.52464762330.081312835450.80491409003 復(fù)化數(shù)值求積公式3.1 復(fù)化數(shù)值求積法無論用newtoncotes求積公式或causs型求積公式，提高數(shù)值積分精度的一個途徑是增加求積節(jié)點數(shù)目。當增大時，newtoncotes公式的數(shù)值穩(wěn)定性變差，也不能保證能提高精度而causs型救只公式的causs點、求積系數(shù)通常是無理數(shù)，查找、計算都不方便。當?shù)馁x值不太復(fù)雜時，提高數(shù)值積分精度的另一個途徑是利用復(fù)化求積公式。復(fù)化求積公式的幫派則是把求積區(qū)間進行等距細分：在每個小區(qū)間上用相同的“基本”求積公式

16、計算出的近似值。并取當權(quán)函數(shù)時，不易構(gòu)造復(fù)化求積公式。下面講座一些常用復(fù)化求積公式。3.2 復(fù)化梯形公式記在上采用梯形公式得即復(fù)化梯形求積公式為 （3.1）設(shè)由得定理7 若，則復(fù)化梯形公式的余項為 （3.2）及漸近估計式 （3.3）當區(qū)間細分節(jié)點加密一倍時，得 其中為復(fù)化中矩求積公式.3.3復(fù)化simpson公式在每個小區(qū)間上采用simpson公式,得復(fù)化simpson求積公式 (3.4)利用復(fù)化梯形公式和復(fù)化中矩公式有 其中 對應(yīng)于復(fù)化simpson公式有如下余項定理. 定理8 當時,復(fù)化simpson公式的余項有表達式 (3.5)及漸近估計式(3.6)類似地我們可以建立復(fù)化cotes公式,

17、復(fù)化causs-legendre求積公式等,同時給出相應(yīng)的余項估計.例7 利用9點函數(shù)值,用復(fù)化梯形公式和復(fù)化simpson公式計算解 經(jīng)計算得t8=0.9456909,s4=0.9460832,c2=0.9460829三種方法所用函數(shù)值個數(shù)一樣多,與積分準確值0.9460831相比較,復(fù)化simpson公式的結(jié)果與復(fù)化梯形公式的結(jié)果相比,復(fù)化simpson公式的結(jié)果要準確得多.故在實際使用中,復(fù)化simpson公式應(yīng)用較廣泛.3.4復(fù)化求積公式的收斂階對上的任何連續(xù)函數(shù),.都有但對代數(shù)多項式因此復(fù)化求積公式不能用代數(shù)精度來決定其優(yōu)劣.對復(fù)化求積公式我們用收斂階來刻劃其收斂性.定義4 設(shè)是將等

18、分,用某一基本求積公式生成的復(fù)化求積公式,我們稱該復(fù)化求積公式具有收斂階，若對充分光滑的被積函數(shù)，有 （3.7）其中獨立于,依賴于.根據(jù)定義,復(fù)化梯形公式的收斂階是2(當時收斂階大于2);復(fù)化simpson公式的收斂階是4(當時大于4);復(fù)化節(jié)點causs-legendre求積公式的收斂階為.收斂階越高,當區(qū)間劃分加密時,積分近似值就越精確. 4 外推方法在用復(fù)化梯形求積公式時,記為,其中.利用定積分定義有在數(shù)值計算中,經(jīng)常會遇到類似情況:精確值是所要求的,但不能用有限計算量算出來,而對某些,卻可以很方便地計算出來.如何從已知的推出的近似值,為些介紹數(shù)值計算中的重要方法;外推方法.4.1外推原

19、理定理9 若逼近有下述余項展開 (4.1)其中,設(shè)為相近的互異正數(shù),則可用 (4.2)來近似，其中滿足 （4.3）而且. 對于的情況,我們有：定理10 若逼近的余項能寫成漸近形式 （4.4）及是獨立于的常數(shù)，則由 （4.5）定義的序列隨增大以更快的速度收斂于： （4.6）其中 （4.7）定理10也稱richardson外推法。42復(fù)化梯形公式余項的漸近展開利用euler-maclaurin求和公式，可以證明復(fù)化梯形公式的余項具有漸近展開。定理11 若則有 (4.8) 其中為數(shù)，而 （4.9）若記 從而可用richardson外推法提高精度。43 romberg算法在復(fù)化梯形公式中，選取記注意到

20、得 （4.10）計算順序如表4-7所示。表4-7 0 1 2 3 0123romberg算法：1輸入外推次數(shù)（一般取為3），控制精度；2置計算，??；3計算；4對進行外推計算 ；5若，輸出數(shù)值積分值，停機；6置，7轉(zhuǎn)3。在romberg算法中，第一列對應(yīng)于復(fù)化梯形序列，第二列對應(yīng)于復(fù)化simpson序列，第三列對應(yīng)于cotes序列，第四列稱為romberg序列。在實際使用中常常只計算到第4段列（即?。?，更高的列較少用。romberg算法中止準則，一般取同列或同行相鄰兩高值的的誤差絕以值小于事先給定的精度要求。romberg算法是數(shù)值穩(wěn)定的，且對任意連續(xù)函數(shù)，都能保證數(shù)值積分收斂到準確值。romb

21、erg算法程序簡單，當函數(shù)值不太復(fù)雜時，romberg算法是常用的實用方法。5 自適應(yīng)求積方法51自適應(yīng)計算問題若要計算要求誤差不超過，我們可以用newton-cotes公式，gauss型求積公式，復(fù)化求積公式,romberg算法等來實現(xiàn)。當充分光滑時，利用余項公式可以確定或區(qū)間等公數(shù)。這兒有些不足之外，首先高階導(dǎo)數(shù)不易估計，即使給出了估計，估計式也把誤差放大到一個誤差限；其次上述求積公方法全是把被積函數(shù)在整個區(qū)間上作整體處理的。函數(shù)在上性質(zhì)可能差異很大。例如在不等長分劃下，若每個小區(qū)間上用simpson公式求積，有一個明顯的啟示是在取值小的區(qū)域，步長可以取大些，而在取值大的區(qū)域，步長應(yīng)當取小

22、些。本節(jié)介紹的自適算法，就是根據(jù)在不同子區(qū)間的性質(zhì)作不同的處理，使在較少的函數(shù)計算前提下達到誤差要求。在上面我們對每個小區(qū)間采用simpson公式為基本求積公式，也可用其它的基本求積公式來實現(xiàn)。52自適應(yīng)算法將要計算的積公記為同理，上積分記為 （5.2）在上某一基本求積公式記為 ， （5.3）設(shè)基本求積公式的代數(shù)精度為，則有 （5.4）其中是獨立于和具的常數(shù)。將等分成兩個小區(qū)間，每個小區(qū)間，每個區(qū)間上用上述同一基本求積公式，得。其余項為（5.5）綜合（5.4），（5.5）得 （5.6）記，設(shè)預(yù)先規(guī)定誤差。階代數(shù)精度基本求積公式的自適應(yīng)算法從開始計算，分別計算出和，若，則取的近似值為，計算結(jié)束。

23、8 數(shù)值微分列表函數(shù)的數(shù)值微分多取插值函數(shù)微分；非列表函數(shù)的數(shù)值微分多取差分近似。數(shù)值微分的數(shù)值穩(wěn)定性差，經(jīng)常利用外推法來提高精度。當同時計算等距節(jié)點上的導(dǎo)數(shù)，隱式方法具有較高的精度。81插值函數(shù)法設(shè)是的lagrange插值多項式或hermite插值多項式，在插值區(qū)間內(nèi)，一種很直觀的數(shù)值微分公式是取 （8.1）1二點公式設(shè)于是 (8.2)其中。利用廣義的均差微公公式（證明從略）我們有當時，有特別地 （8.3）2.三點公式設(shè)充分光滑，利用關(guān)于三點的lagrange插值多項式,我們可建立數(shù)值微分公式 （8.4） （8.5）當時：（8.6）及 (8.7)從(8.3),(8.7)可以看出二點公式中和三

24、點公式中具有較高的精度,這二個節(jié)點(指)處在插值節(jié)點的“中心”位置。利用函數(shù)的三次樣條插值函數(shù)，我們也可以建立數(shù)值微分公式 （8.8）若記.當時,具體微分公式為由于與可能不一樣,但相差不大.對于數(shù)值微分公式（8.8）的余項有：定理12 若是的一型或二型邊值插值三次樣條函數(shù)，則成立 (8.9)其中82差分算子近似微分算子法求解微分方程的一種主要方法是差分法。差分法就是利用函數(shù)的差分逼近函數(shù)的微分或偏微分。利用向前差分算子，向后差分算子或中心差分算子來近似微分算子，可以很方便地建立數(shù)值微分公式。例如對函數(shù)的一階微分有 （8.10）例如對將所有微分用中心差分代替，可建立數(shù)值微分公式 = 特別地當時， （8.11）若用逼近也可建立公式由（83）和（87）可以看出，利用中心差