歐拉法

首先先來介紹歐拉法，如果對(duì)于一條曲線，用$\Delta$$t$來將線段劃分成許多小段，那么我們就可以近似認(rèn)為每一段的斜率都是常數(shù)，且對(duì)于第$i$段的斜率表示為
$$\frac{d{y}_i}{dt}=\frac{y_{i+1}?y_i}{\Delta t}$$
以最為常見的一階微分方程為例
$$\frac{dy}{dt}=a_{0}y+f(t)$$
我們?cè)倮脷W拉法求解的時(shí)候需要第$i$段的斜率，因此利用上式，帶入$(t_i,y_i)$求得$\frac{d{y}_i}{dt}$，綜合可得
$$\frac{y_{i+1}?y_i}{\Delta t}=a_0{y}_i+f(t_i)$$
再變一下型
$$y_{i+1}=(a_0\Delta t+1)y_i+f(t_i)\Delta t$$
這就是歐拉法，其誤差的大小取決于你$\Delta$$t$選擇的大小。

改進(jìn)的歐拉法

首先計(jì)算出第$i$個(gè)點(diǎn)的斜率
$$\frac{d{y}_i}{dt}=a_0{y}_i+f(t_i)$$
然后再利用
$$y_{i+1}=\frac{d{y}_i}{dt}\Delta t+y_i$$
求出第$i+1$個(gè)點(diǎn)處的斜率
$$\frac{d{y}_{i+1}}{dt}=a_0{y}_{i+1}+f(t_i)$$
這樣我們?cè)谟?jì)算這一段時(shí)就會(huì)有兩個(gè)斜率可以用，所以就取這兩個(gè)斜率的平均值來作為這個(gè)線段的斜率
$$y_{i+1}=y_i+\frac{\Delta t}{2}(\frac{d{y}_{i+1}}{dt}+\frac{d{y}_i}{dt})$$

四階龍格-庫塔積分法

為了方便這里設(shè)$s_i=f(y_i,t_i)=a_0{y}_i+f(t_i)$，那么先計(jì)算
$$s_i=f(y_i,t_i)$$
再計(jì)算
$$s_{i+1}=f(y_i+\frac{\Delta t}{2}{s}_i,t_i+\frac{\Delta t}{2})??沿著s_i這個(gè)斜率向前走\frac{\Delta t}{2}處的斜率$$
再計(jì)算
$$s_{i+2}=f(y_i+\frac{\Delta t}{2}{s}_{i+1},t_i+\frac{\Delta t}{2})$$
$$s_{i+3}=f(y_i+\Delta t{s}_{i+2},t_i+\Delta t)$$
得到這4個(gè)斜率后再用下式計(jì)算
$$y_{i+1}=y_i+\frac{\Delta t}{6}(s_i+2s_{i+1}+2s_{i+2}+s_{i+3})$$
這種方法比前兩種算法都來得精準(zhǔn)，但是增加了計(jì)算量，以及計(jì)算程序變得復(fù)雜。

總結(jié)

以上是生活随笔為你收集整理的解微分方程数值解法（理论部分）的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。