2-數(shù)值積分.ppt

1 第2章數(shù)值積分 1 2 2 1引言 利用牛頓 萊布尼茲 Newton Leibniz 公式 2 1 解決函數(shù)在上的積分問題在理論和應(yīng)用上都有重大的意義 然而 在實際問題中 往往會遇到一些困難 有些形式上較簡單的函數(shù) 其原函數(shù)不易求出或不能用初等函數(shù)表示成有限形式 有些被積函數(shù)的原函數(shù)過于復(fù)雜 而有些函數(shù)的函數(shù)值是由實驗 觀測等方法得出 并沒有給出具體的解析表達式 這些情形說明公式 5 1 在應(yīng)用上是有局限性的 因此研究定積分的數(shù)值計算問題就顯得十分必要 本章主要介紹一些常用的數(shù)值積分方法 包括梯形積分法 辛卜生積分法 變步長積分法 牛頓 柯特斯積分法 高斯積分法 龍貝格積分法 2 3 2 2梯形積分法 2 2 1梯形積分法的基本思想梯形積分法的基本思想 在積分區(qū)間上 根據(jù)給定的插值條件和 構(gòu)造一個一次二項式 并以的積分值近似地代替 從幾何角度而言 是以梯形面積近似地代替曲邊梯形的面積 3 圖2 1 4 2 2 2梯形求積公式依據(jù)梯形積分法的基本思想 將區(qū)間分成個相等的小區(qū)間 則每個小區(qū)間的長度為 對每個小區(qū)間均實施如下的梯形求積 將這些小梯形的求積值加起來 可以得到如下梯形求積公式 4 5 2 2 3實現(xiàn)梯形積分法的基本步驟 1 輸入?yún)^(qū)間的端點值以及分割數(shù) 2 將區(qū)間等分成個小區(qū)間 每一個小區(qū)間的長度 3 計算每一個等分點的函數(shù)值 4 計算 5 輸出的值 6 結(jié)束 5 6 例2 1使用梯形求積公式求下列定積分的值 6 7 辛卜生積分法的基本思想 在積分區(qū)間上 根據(jù)給定的插值條件 和 構(gòu)造一個二次插值求積多項式 并以的積分值近似地代替 從幾何角度而言 是用過三點的拋物線面積近似地代替積分的曲邊面積 2 3辛卜生 Simpson 積分法2 3 1辛卜生積分法的基本思想 7 圖2 3 8 2 3 2辛卜生求積公式依據(jù)辛卜生積分法的基本思想 將區(qū)間分成 必須是偶數(shù) 個相等的小區(qū)間 則每個小區(qū)間的長度為 在小區(qū)間均實施如下的辛卜生求積 將這些求積值加起來 可以得到如下辛卜生求積公式 其中 為奇數(shù)項的函數(shù)值之和 為偶數(shù)項的函數(shù)值之和 8 9 2 3 3實現(xiàn)辛卜生積分法的基本步驟 1 輸入?yún)^(qū)間的端點的值以及分割數(shù) 2 將區(qū)間等分成個小區(qū)間 每一個小區(qū)間的長度 3 計算每一個等分點的函數(shù)值 4 計算 計算奇數(shù)項的函數(shù)值之和 計算偶數(shù)項的函數(shù)值之和 5 計算 6 輸出的值 7 結(jié)束 9 10 例2 2使用辛卜生求積公式求下列定積分的值 10 11 2 4變步長求積分法 2 4 1變步長求積分法的基本思想變步長求積分法是以梯形公式為基礎(chǔ) 逐步改變步長 以達到預(yù)先所要求的精度 變步長求積分法主要有變步長梯形求積分法和變步長辛卜生求積分法 本節(jié)我們將介紹這兩種方法 11 12 2 4 2變步長梯形求積分法變步長梯形求積分法的基本過程 1 利用梯形公式 將積分區(qū)間等分 即其中 2 將每一個求積小區(qū)間再二等分一次 即由原來的等分變成等分 則有 12 13 其中 為再二等分一次后新增加的結(jié)點 它們都是原來各小區(qū)間的中點 由上式可以看出 在對每一個小區(qū)間在二等分后 在積分值的第一項中只包含再二等分之前的各結(jié)點上的函數(shù)值 并且第一項的值正好是再二等分之前積分值的一半 顯然 這一項中所包含的函數(shù)值就不必計算了 再二等分后需要計算的函數(shù)都包含在第二項中 它們都是二等分后出現(xiàn)的新的結(jié)點 因此有 3 若 二等分后的積分值就是最后的結(jié)果 否則保存當前的等分數(shù) 積分值與步長 即轉(zhuǎn)到第 2 步繼續(xù)做二等分處理 13 14 2 4 3實現(xiàn)變步長梯形積分法的基本步驟 1 輸入?yún)^(qū)間的端點的值以及容許誤差 2 計算區(qū)間的長度 3 首先將區(qū)間進行一等分 并設(shè) 4 再將每一個求積小區(qū)間 由原來的等分變成等分 5 計算 6 成立 則繼續(xù) 否則 轉(zhuǎn) 4 7 輸出的值 8 結(jié)束 14 15 2 4 4變步長辛卜生求積分法變步長辛卜生求積分法的基本過程 1 利用梯形公式 將積分區(qū)間一等分 2 將其中的每一個求積小區(qū)間再二等分一次 3 根據(jù)上面兩式和 可以推導(dǎo)出如下的變步長辛卜生求積公式 進一步得到再二次等分一次后的變步長辛卜生求積公式為 4 若 二等分后的積分值就是最后的結(jié)果 否則保存當前的變步長梯形積分值 等分數(shù) 積分值與步長 轉(zhuǎn)到第 2 步繼續(xù)做二等分處理 15 16 2 4 5實現(xiàn)變步長辛卜生積分法的基本步驟 16 17 17 18 例2 4使用變步長辛卜生求積分法求下列定積分的值 18 19 2 5牛頓 柯特斯 Newton Cotes 積分法 2 5 1牛頓 柯特斯積分法的基本思想牛頓 柯特斯積分法的基本思想 用高次的插值求積多項式去逼近被積函數(shù) 以獲得高精度的積分值 事實上 梯形積分是當時的牛頓 柯特斯積分 辛卜生積分是當時的牛頓 柯特斯積分 它們都是牛頓 柯特斯積分的特例 19 20 2 5 2牛頓 柯特斯求積公式下面給出三到五階牛頓 柯特斯求積公式 20 21 實現(xiàn)三階牛頓 柯特斯求積公式的基本步驟如下 21 22 例2 5使用牛頓 柯特斯求積公式求下列定積分的值 22 23 2 6龍貝格 Romberg 積分法 2 6 1龍貝格積分法的基本思想前面講述的各種求積方法是插值求積的思想 而龍貝格積分法的基本思想是 使用一個諸如梯形求積法等代數(shù)精度較低的求積公式 相繼以步長和求得定積分的兩個近似結(jié)果 然后再做它們適當?shù)木€性組合 就可以得到一個代數(shù)精度更高的公式 23 24 根據(jù)復(fù)化梯形公式的余項表達式可知 2 6龍貝格 Romberg 積分法 25 即 依此類推 26 這樣直接用計算結(jié)果來估計誤差的方法通常稱作誤差的事后估計法 27 同理由復(fù)化辛普森公式的余項 可得 由復(fù)化Cotes公式的余項 得 28 例 根據(jù)如下函數(shù)值表 利用復(fù)化梯形公式計算積分的近似值 要求誤差不超過 原積分的精確值為 29 然后將區(qū)間二等分 利用梯形公式的遞推公式求出 遞推公式 進一步二分積分區(qū)間 類似可求出 如此不斷二分并利用遞推公式 可得下表中的結(jié)果 解 先在整個區(qū)間上用梯形公式 30 k表示二分次數(shù) 區(qū)間數(shù) 由表中可以看出 對分8次和對分7次之間的差 因而是滿足精度要求的解 31 收斂速度慢 對于復(fù)化辛蒲生公式 柯特斯公式可以類似得到 32 加速收斂 應(yīng)用步長逐次減半得到的復(fù)化梯形值 復(fù)化辛蒲生值 復(fù)化柯特斯值與精確值的比較 33 例將以上三個加速公式用于求 從表中可以看出三次加速求得R1 0 9460831每位數(shù)字都是有效數(shù)字 34 上述用若干個積分近似值算出更精確的積分近似值的方法 稱之為外推法 4個積分值序列 梯形值序列 辛蒲生值序列 龍貝格值序列 柯特斯值序列 35 外推法的計算步驟 36 例 利用龍貝格積分法式計算積分要求精確到小數(shù)點后面7位 解 根據(jù)龍貝格積分法計算得 37 具體結(jié)果見下表 38 精確值為0 91629073187415506518352721176801 39 2 6 2實現(xiàn)龍貝格積分法的基本步驟 1 輸入?yún)^(qū)間的端點的值 最大迭代次數(shù)以及容許誤差 2 計算區(qū)間的長度 3 用梯形積分法計算積分近似值 4 對計算對計算 如果 則退出循環(huán) 5 如果 則繼續(xù) 否則輸出無解信息 轉(zhuǎn) 7 6 輸出的值 7 結(jié)束 39 40 本章小結(jié) 本章首先簡要介紹了數(shù)值積分在實際應(yīng)用中的重要性 并對在本章中介紹的各種數(shù)值積分的基本思想作了詳細的說明 作為數(shù)值積分的方法 介紹了梯形積分法 辛卜生積分法 變步長積分法 牛頓 柯特斯積分法 高斯積分法 龍貝格積分法 梯形積分法和辛卜生積分法是最基本的求積方法 其不足之處是當積分區(qū)間的長度較大時 直接使用這兩種方法 則計算結(jié)果的精度難以滿足要求 因此 為了提高精度 引進了不斷二等分每個子區(qū)間的變步長梯形積分法和變步長辛卜生積分法 牛頓 柯特斯積分法是一種可以獲得高精