数値計算 差分法

 

目次

数値微分

微分はコンピュータでは正確な計算で行えないため、差分法で近似して計算する。

1階微分

前進差分:$\displaystyle\frac{df}{dx}\simeq \displaystyle\frac{f(x+h)-f(x)}{h}$

後退差分:$\displaystyle\frac{df}{dx}\simeq \displaystyle\frac{f(x)-f(x-h)}{h}$

中心差分:$\displaystyle\frac{df}{dx}\simeq \displaystyle\frac{f(x+h)-f(x-h)}{2h}$

 

テイラー展開すると、前進差分(後退差分)は次のように1次のオーダーの誤差になる。

$$\displaystyle\frac{df}{dx}\simeq \displaystyle\frac{f(x+h)-f(x)}{h}\simeq f^{\prime}(x)+\displaystyle\frac{h}{2}f^{\prime\prime}(x)+\cdots$$

中心差分は2次のオーダーの誤差になる。

$$\displaystyle\frac{df}{dx}\simeq \displaystyle\frac{f(x+h)-f(x-h)}{2h}\simeq f^{\prime}(x)+\displaystyle\frac{h^2}{3!}f^{\prime\prime\prime}(x)+\cdots$$

 

2階微分

次の式のように、近似計算する。

$$\displaystyle\frac{d^2 f}{dx^2}\simeq \displaystyle\frac{f(x+h)+f(x-h)-2f(x)}{2h}$$

 

 

 

タイトルとURLをコピーしました