尤拉公式推導及實現程式碼

來源:酷知科普網 1.34W

尤拉方法的基本思想是在小區間上用數值微分的前差公式代替方程左端的導數,右端函式f(x,y(x))中x取中的某一點。

操作方法

(01)用數值前差表示左端

尤拉公式推導及實現程式碼

(02)從微分方程開始推導

尤拉公式推導及實現程式碼 第2張

(03)由(1)(2)可以得到

尤拉公式推導及實現程式碼 第3張

(04)推出向前尤拉公式

尤拉公式推導及實現程式碼 第4張

(05)以此類推得到向後尤拉公式和中心尤拉公式

尤拉公式推導及實現程式碼 第5張

(06)向前尤拉公式MATLAB程式碼Function [x,y]=uler(fun,ts,y0,h)%fun是常微分函式,ts為解區間y0為初值,h為步長X=ts(1):h:ts(2)Y(1)=y(0);For n=1:length(x);Y(n+1)=y(n)+h*fun(x(n),y(n));end

特別提示

向後尤拉法使用起來比較麻煩,可是穩定性比向前尤拉好

熱門標籤