北京交通大學數(shù)值分析 6微分方程數(shù)值解習題課_第1頁
北京交通大學數(shù)值分析 6微分方程數(shù)值解習題課_第2頁
北京交通大學數(shù)值分析 6微分方程數(shù)值解習題課_第3頁
北京交通大學數(shù)值分析 6微分方程數(shù)值解習題課_第4頁
北京交通大學數(shù)值分析 6微分方程數(shù)值解習題課_第5頁
已閱讀5頁,還剩4頁未讀, 繼續(xù)免費閱讀

下載本文檔

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

文檔簡介

微分方程初值問題數(shù)值解習題課一、應用向前歐拉法和改進歐拉法求由如下積分所確定的函數(shù)y在點x=0.5,1.0,1.5的近似值。解:該積分問題等價于常微分方程初值問題其中h=0.5。其向前歐拉格式為改進歐拉格式為將兩種計算格式所得結(jié)果列于下表向前歐拉法改進歐拉法000010.50.50.4447021.00.889400.7313731.51.073340.84969二、應用4階4步阿達姆斯顯格式求解初值問題取步長h=0.1.解:4步顯式法必須有4個起步值,已知,其他3個用4階龍格庫塔方法求出。本題的信息有:步長h=0.1;結(jié)點;經(jīng)典的4階龍格庫塔公式為算得,,4階4步阿達姆斯顯格式由此算出三、用Euler方法求問步長應該如何選取,才能保證算法的穩(wěn)定性?解:本題本題的絕對穩(wěn)定域為得,故步長應滿足求梯形方法的絕對穩(wěn)定域。證明:將Euler公式用于試驗方程,得到整理設計算時有舍入誤差,則有據(jù)穩(wěn)定性定義,要想,只須因此方法絕對穩(wěn)定域為復平面的整個左半平面(?),是A-穩(wěn)定的。五、對初值問題證明:用梯形公式求得的數(shù)值解為并證明當步長時,收斂于該初值問題的精確解證明:由梯形公式,有整理,得由此遞推公式和初值條件,有,則有在區(qū)間上有,步長,由前面結(jié)果有由x的任意性,得所證。六。常微分方程初值問題的單步法為試求其局部截斷誤差主項并回答它是幾階精度的?解該單步公式的局部截斷誤差是故局部截斷誤差主項是,方法是一階的。七、對于微分方程,已知在等距結(jié)點處的y的值為,h為步長。試建立求的線性多步顯格式與與隱格式。解:取積分區(qū)間,對兩端積分:對右端作的二次插值并積分得到線性4步顯格式若對右端在兩點上作線性插值并積分,有由此產(chǎn)生隱格式八、證明線性多步法存在的一個值,使方法是4階的。解

溫馨提示

  • 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

提交評論