變步長的龍格庫塔法.ppt_第1頁
變步長的龍格庫塔法.ppt_第2頁
變步長的龍格庫塔法.ppt_第3頁
變步長的龍格庫塔法.ppt_第4頁
變步長的龍格庫塔法.ppt_第5頁
已閱讀5頁，還剩11頁未讀， 繼續(xù)免費閱讀

下載本文檔

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

文檔簡介

變步長的龍格 庫塔方法 公式 7 4 6變步長的龍格 庫塔法在微分方程的數(shù)值解中 選擇適當(dāng)?shù)牟介L是非常重要的 單從每一步看 步長越小 截斷誤差就越小 但隨著步長的縮小 在一定的求解區(qū)間內(nèi)所要完成的步數(shù)就增加了 這樣會引起計算量的增大 并且會引起舍入誤差的大量積累與傳播 因此微分方程數(shù)值解法也有選擇步長的問題 以經(jīng)典的四階龍格 庫塔法 7 20 為例 從節(jié)點xi出發(fā) 先以h為步長求出一個近似值 記為 由于局部截斷誤差為 故有 當(dāng)h值不大時 式中的系數(shù)c可近似地看作為常數(shù) 然后將步長折半 即以為步長 從節(jié)點xi出發(fā) 跨兩步到節(jié)點xi 1 再求得一個近似值 每跨一步的截斷誤差是 因此有 這樣 由此可得 這表明以作為的近似值 其誤差可用步長折半前后兩次計算結(jié)果的偏差 來判斷所選步長是否適當(dāng) 當(dāng)要求的數(shù)值精度為 時 1 如果 反復(fù)將步長折半進(jìn)行計算 直至 為止 并以上一次步長的計算結(jié)果作為 這種通過步長加倍或折半來處理步長的方法稱為變步長法 表面上看 為了選擇步長 每一步都要反復(fù)判斷 增加了計算工作量 但在方程的解y x 變化劇烈的情況下 總的計算工作量得到減少 結(jié)果還是合算的 其中 i i 1 m i i 2 m 和 ij i 2 m j 1 i 1 均為待定系數(shù) 確定這些系數(shù)的步驟與前面相似 2Runge KuttaMethod 高階Runge KuttaMethod Gill公式 4階經(jīng)典龍格 庫塔公式的一種改進(jìn) 2Runge KuttaMethod 最常用為四級4階經(jīng)典龍格 庫塔法 ClassicalRunge KuttaMethod 2Runge KuttaMethod 由于龍格 庫塔法的導(dǎo)出基于泰勒展開 故精度主要受解函數(shù)的光滑性影響 對于光滑性不太好的解 最好采用低階算法而將步長h取小 2Runge KuttaMethod 變步長的Runge KuttaMethod Q 由局部截斷誤差可以看出 步長h越小 局部截斷誤差越小 但步長減小 在一定求解范圍 區(qū)間 內(nèi)要完成的步數(shù)就增加了 步數(shù)增加會引起計算量增大 導(dǎo)致舍入誤差積累 因此要選取適當(dāng)?shù)牟介L 選擇步長時要考慮兩個問題

溫馨提示

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

評論

0/150

提交評論