数値解析 第5回 (2) １階微分
前進差分・後退差分
まずは１階微分の近似式から。
アイデア 微分の定義は
$\dps{y'(x)=\lim_{t \rightarrow 0} \frac{y(x+t)-y(x)}{t}}$
なので、刻み幅 $h \gt 0$ が十分小さければ、$t=h$ として
$y'(x)$ ≒ $\dps{\frac{y(x+h)-y(x)}{h}}$
であろうし、また $t=-h$ として
$y'(x)$ ≒ $\dps{\frac{y(x-h)-y(x)}{-h}=\frac{y(x)-y(x-h)}{h}}$
であろうと思われます。
$y'(x)$ の近似式として、
Def.1
$\dps{\frac{\mathit{\Delta} y}{\mathit{\Delta} x}=\frac{y(x+h)-y(x)}{h}}$ $\cdots\cdots$ $(1)$
を前進差分 ( 前方差分 )、
$\dps{\frac{\nabla y}{\nabla x}=\frac{y(x)-y(x-h)}{h}}$ $\cdots\cdots$ $(2)$
を後退差分 ( 後方差分 )と呼ぶ。
$h$ を大き目にして誇張して描けば図のようになります：
前進差分・後退差分の誤差評価
Prop.2
$\dps{\frac{\mathit{\Delta} y}{\mathit{\Delta} x}}$, $\dps{\frac{\nabla y}{\nabla x}}$ の近似誤差は、いずれもおおよそ
$\dps{\left|\,\frac{1}{2}y''(x)h\,\right|}$
であり、刻み幅 $h$ の 1 乗オーダー $O(h)$ になる。
証明 $x$ での $y(x)$ のテイラー展開から
$\dps{y(x+h)=y(x)+y'(x)h+\frac{1}{2}y''(x)h^2+O(h^3)}$
従って
$\dps{\frac{\mathit{\Delta} y}{\mathit{\Delta} x}
=\frac{y(x+h)-y(x)}{h}=y'(x)+\frac{1}{2}y''(x)h +O(h^2)}$ $\cdots\cdots$ $(3)$
∴ $\dps{\left|\,\frac{\mathit{\Delta} y}{\mathit{\Delta} x}-y'(x)\,\right|
=\left|\,\frac{1}{2}y''(x)h +O(h^2)\,\right|}$ ≒ $\left|\,\frac{1}{2}y''(x)h\,\right|$
$h$ の代わりに $-h$ を入れた式
$\dps{y(x-h)=y(x)-y'(x)h+\frac{1}{2}y''(x)h^2+O(h^3)}$
からは
$\dps{\frac{\nabla y}{\nabla x}
=\frac{y(x)-y(x-h)}{h}=y'(x)-\frac{1}{2}y''(x)h +O(h^2)}$ $\cdots\cdots$ $(4)$
となって同様です。(証明終)
中心差分
前進差分・後退差分はいまいち精度が良くありません。そこで
次のアイデア $(3)$ と $(4)$ の平均を取ると $y''$ の項が消えます：
$\dps{
\frac{1}{2}\left(
\frac{\mathit{\Delta} y}{\mathit{\Delta} x}
+\frac{\nabla y}{\nabla x}
\right)
=\frac{y(x+h)-y(x-h)}{2h}=y'(x)+O(h^2)}$ $\cdots\cdots$ $(5)$
名前を付けて
Def.3
$\dps{\frac{\delta y}{\delta x}=\frac{y(x+h)-y(x-h)}{2h}}$ を中心差分と呼ぶ。
中心差分は $h$ が多少大きくてもかなり $y'(x)$ に近いことが図からもわかります：
そして、テイラー展開を $h^3$ の項まで書けば誤差評価ができます。
Prop.4
$\dps{\frac{\delta y}{\delta x}}$ の近似誤差はおおよそ
$\dps{\left|\,\frac{1}{6}y'''(x)h^2\,\right|}$
であり、刻み幅 $h$ の 2 乗オーダー $O(h^2)$ になる。