《數(shù)值計算方法》課件7數(shù)值積分與數(shù)值微分_第1頁
《數(shù)值計算方法》課件7數(shù)值積分與數(shù)值微分_第2頁
《數(shù)值計算方法》課件7數(shù)值積分與數(shù)值微分_第3頁
《數(shù)值計算方法》課件7數(shù)值積分與數(shù)值微分_第4頁
《數(shù)值計算方法》課件7數(shù)值積分與數(shù)值微分_第5頁
已閱讀5頁,還剩17頁未讀, 繼續(xù)免費閱讀

下載本文檔

版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認領(lǐng)

文檔簡介

7.1

數(shù)值求積公式與代數(shù)精度

數(shù)值求積的必要性由微積分理論可知,只要被積函數(shù)在區(qū)間[a,b]連續(xù),就可以使用牛頓—萊布尼茲(Newton-Leibniz)公式計算定積分。然而,在許多實際問題中,這種解析方法是無能為力的,或者是非常麻煩的。體現(xiàn)在:(1)找不到原函數(shù)F(x),如第一個例子中的第二類橢圓積分等,有些被積函數(shù)是以表格形式給出,沒有有限的解析表達式。(2)雖然可以找到原函數(shù)F(x),但它比被積函數(shù)復雜得多,而且有時難以給出最后的數(shù)值結(jié)果。(3)除一些特殊的無窮積分外,通常很難求無窮積分的值。7.1

數(shù)值求積公式與代數(shù)精度

數(shù)值求積基本思想求積區(qū)間離散化借助函數(shù)插值求定積分得到近似求積公式及余項7.1.1

數(shù)值求積的基本思想數(shù)值求積公式積分余項

基于積分中值定理構(gòu)造的求積公式左、右矩形公式中點公式梯形公式辛卜生公式7.1

數(shù)值求積公式與代數(shù)精度7.1.2

求積公式的代數(shù)精度定義7-1

m次代數(shù)精度定理7-1例7-1

確定求積公式的代數(shù)精度例7-2

確定求積公式的代數(shù)精度例

確定下面2個求積公式的代數(shù)精度7.1.3

插值型求積公式定理7-2p1487.1

數(shù)值求積公式與代數(shù)精度求積系數(shù)求積余項容易得到7.1.4

求積公式的收斂性與穩(wěn)定性定義7-2

求積公式收斂定義7-3

求積公式穩(wěn)定定理7-3

若Ak>0,則求積公式是穩(wěn)定的。7.2.1

牛頓—柯特斯求積公式與求積系數(shù)7.2

牛頓—柯特斯求積公式不難推得n=1,2,4時的梯形、辛卜生和柯特斯求積系數(shù)。表7-1牛頓--柯特斯系數(shù)表。觀察柯特斯系數(shù)的對稱性、和為1,n>7時有負值,影響公式的穩(wěn)定性。7.2.2

偶數(shù)階牛頓—柯特斯公式的代數(shù)精度7.2牛頓—柯特斯求積公式具有3次代數(shù)精度辛卜生公式被廣泛使用。

7.2.3

低階牛頓—柯特斯公式的余項7.2牛頓—柯特斯求積公式

梯形求積公式的余項7.2.3

低階牛頓—柯特斯公式的余項7.2牛頓—柯特斯求積公式

辛卜生求積公式的余項同理可求出柯特斯求積公式的余項7.2.4

復化求積公式及其余項7.2牛頓—柯特斯求積公式問題:牛頓—柯特斯公式無論從收斂性、還是穩(wěn)定性上都不能得到保證解決方法:將積分區(qū)間分成n等分,使每等分的寬度盡可能的小每個小區(qū)間使用低階的牛頓—柯特斯公式1.復化梯形公式及余項1次代數(shù)精度7.2.4

復化求積公式及其余項7.2牛頓—柯特斯求積公式2.復化辛卜生公式及余項3次代數(shù)精度7.2.4

復化求積公式及其余項7.2牛頓—柯特斯求積公式3.復化柯特斯公式及余項5次代數(shù)精度復化辛卜生算法流程圖7-1例7-3例7-4

復化求積公式收斂的階7.3龍貝格求積公式7.3.1

變步長求積公式7.3龍貝格求積公式例7-5P1567.3.2

龍貝格求積公式7.3龍貝格求積公式

修正梯形公式

修正辛卜生公式不難驗證,修正梯形公式正好就是辛卜生公式Sn,可以設(shè)想對辛卜生公式修正得到柯特斯公式,視為對梯形公式的二次修正。7.3.2

龍貝格求積公式7.3龍貝格求積公式

龍貝格(Romberg)公式

龍貝格算法示意圖對柯特斯公式進行修正,得到精度更好的求積公式,但它不是n=8時的牛頓—柯特斯公式,由龍貝格首先提出,稱為龍貝格公式。

龍貝格算法流程框圖例7-6P1587.4.1

高斯求積公式與高斯點7.4高斯求積公式定理7-5插值型求積公式能達到的最高代數(shù)精度為2n+1次。定義7-4

具有2n+1次代數(shù)精度的求積公式稱為高斯求積公式,其求積節(jié)點稱為高斯點。為構(gòu)造具有最高代數(shù)精度的高斯型求積公式,直接用待定系數(shù)法求解n+1個系數(shù)Ak和n+1個求積節(jié)點xk通常是不可能的。通常的方法是:(1)先確定求積區(qū)間上的高斯點xk;(2)在求n+1個求積系數(shù)Ak(待定系數(shù)法和公式法均可)定理7-6高斯點的充要條件。p1607.4.2

高斯求積公式的構(gòu)造構(gòu)造被積區(qū)間[a,b]上的n+1次帶權(quán)正交多項式gn+1(x)求gn+1(x)的n+1個零點,作為高斯點xk用待定系數(shù)法或?qū)窭嗜栈瘮?shù)求積分來確定求積系數(shù)Ak例7-7P1617.4.3

高斯—勒讓德求積公式7.4高斯求積公式定理7-7[-1,1]上的n+1次勒讓德多項式與任一不超過n次的多項式正交。(權(quán)函數(shù)為1)當求積節(jié)點xk為n+1次勒讓德正交多項式的零點時,求積公式稱為高斯—勒讓德求積公式。例7-8高斯—勒讓德一點、兩點、三點公式當積分區(qū)間為[a,b]時,需引進變換則t在區(qū)間[-1,1]上。例7-9P1647.4.4

高斯—切比雪夫求積公式7.4高斯求積公式其中,求積區(qū)間為[-1,1],權(quán)函數(shù)、高斯點分別為7.4.5

高斯—埃爾米特求積公式其中,求積區(qū)間為[-∞,+∞],高斯點為n+1次埃爾米特正交多項式的零點,權(quán)函數(shù)、求積系數(shù)分別為例7-10P1657.4.6

高斯求積公式的余項與穩(wěn)定性7.4高斯求積公式定理7-8帶權(quán)高斯公式的余項

勒讓德公式的余項

切比雪夫公式的余項

埃爾米特公式的余項一、二、三點勒讓德公式的余項定理7-9高斯求積公式是穩(wěn)定的7.5.1

基于Taylor展開的微分公式7.5數(shù)值微分Taylor展開式向前差商公式向后差商公式中心差商公式二階差商公式7.5.1基于Taylor展開的微分公式7.5數(shù)值微分

方法誤差

舍入誤差

整體誤差由截斷誤差項可以看出,前兩個公式的精度為O(h),后兩個公式的精度為O(h2)。如果僅從截斷誤差上考慮,h的冪次越高、h越小,計算精度就越高。從穩(wěn)定性的角度講,還應考慮舍入誤差。

h越小,f(x0+h)、f(x0-h)和f(x0)的值越接近,它們相減后的有效數(shù)字損失越嚴重。在實際計算時,步長h不宜過大,也不宜過小,應綜合考慮截斷誤差和數(shù)值穩(wěn)定性這兩個重要因素。中心差商公式整體誤

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責。
  • 6. 下載文件中如有侵權(quán)或不適當內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論