首先,我們可以知道一階龍格庫塔法的表示式上述式子可以改寫成如下式子
這裡我們可以看成k是在步長為h情況下,函式的增長量,f(x,y)為函式y的斜率。這個式子其實比較好理解的。
然後我們來看二階龍格庫塔的推導過程:
二階情況下一般表示式中就只有k1和k2了,二階龍格庫塔我們可以看成是一階表示式的乙個變形,由原來的乙個k變成了k1,k2,用k1,k2來表示k,c1和c2相當於k1和k2所佔的權重一樣,這樣理解是不是好一點。那k1和k2分別表示什麼??其實k1表示xi點處的增長量,k2表示的是xi+a2*h處的增長量。在這基礎上,因為c1、c2、a2和
這樣二階龍格庫塔的迭代式子就可以表達出來了,寫完收工。
龍格 庫塔(Runge Kutta)法解微分方程
龍格 庫塔 runge kutta 方法是一種在工程上應用廣泛的高精度單步演算法。由於此演算法精度高,採取措施對誤差進行抑制,所以其實現原理也較複雜。該演算法是構建在數學支援的基礎之上的。對於一階精度的尤拉公式有 yi 1 yi h k1 k1 f xi,yi 當用點xi處的斜率近似值k1與右端點x...
4階經典龍格庫塔公式求解微分方程
本次求解微分方程,用4階龍格庫塔方法,該方法的代數精度高,比尤拉法,三階龍格庫塔都高,本次的實現很粗糙,由於本人沒有熟練掌握函式指標的方式,對於不同微分方程,並不能做到乙個普適的輸入版本,所以本演算法只能對本次的微分方程適用,等掌握了函式指標的用法,再來重寫,這次的演算法很簡陋,也沒有做一些越界的異...
二階偏微分方程組 龍格庫塔法 四階龍格庫塔法
建議閱讀原文 預備知識中點法解常微分方程 組 龍格庫塔法是一類數值解微分方程的演算法,其中較常見的是四階龍格庫塔法 這裡不進行推導,僅僅給出公式如下 其中 由以上兩式,不難把該演算法拓展到方程組的情況 對於 元微分方程組 我們可以把該式記為向量函式的形式 現在我們僅需要把式 1 和 式 2 中的所有...