数値解析 第12回 (6) ヤコビ法の収束性
反復行列
Def.12 $A$ の対角部分を
$\dps{
D=\mat{cc}{
\begin{array}{cc}
a_{11} & \\
& a_{22} \\
\end{array}
&
\Huge{0}
\\[-0.5em]
\Huge{0}
&
\begin{array}[b]{cc}
\ddots & \\
& a_{nn} \\
\end{array}
}
}$
とし、
$\dps{H=(h_{ij})=-D^{-1}(A-D)}$
とおく。$H$ をヤコビ法の「反復行列」と呼ぶ。
ちょっと見、複雑そうですが、
$\dps{
h_{ij}=
\left\{
\begin{array}{cl}
\dps{-\left(\frac{a_{ij}}{a_{ii}}\right)} & \quad \mbox{if }\ i \neq j \\
0 & \quad \mbox{otherwise} \\
\end{array}
\right.
}$
です。なぜ反復行列と呼ぶかと言うと、
Th.13 ヤコビ法の $k$ ステップ目の誤差ベクトルを
$\dps{\eee^{(k)}=\xxx^{(k)} - }$ ( 真の $\xxx$ )
と置くと
$\eee^{(k+1)} = H \eee^{(k)}$.
従って
$\eee^{(k)} = H^k \eee^{(0)}$.
証明 ヤコビ法のアイデア は、$D$ を用いて
$D\xxx^{(k+1)} = \bbb - (A-D)\xxx^{(k)}$
と書くことができます。真の解 $\xxx$ は
$D\xxx = \bbb - (A-D)\xxx$
を満たすので、辺々引くと
$D(\xxx^{(k+1)}-\xxx) = - (A-D)(\xxx^{(k)}-\xxx)$
$\therefore\quad D\eee^{(k+1)} = - (A-D)\eee^{(k)}$
$\therefore\quad \eee^{(k+1)} = - D^{-1}(A-D)\eee^{(k)} = H\eee^{(k)}$. (証明終)
ヤコビ法が収束する十分条件
Th.14 $A=(a_{ij})$ が次の条件を満たせばヤコビ法は必ず収束する:
$\forall\, i$ について
$\dps{
\sum_{j \neq i} |\,a_{ij}\,| \lt |\,a_{ii}\,|}$ $\cdots\cdots$ $(\sharp)$
Ex. $A=\mat{rrr}{5 & 1 & 2 \\ 2 & 7 & -3 \\ 1 & -1 & 4 }$ のとき
- $i=1$ : $|\,1\,|+|\,2\,| = 3 \lt |\,5\,|$
- $i=2$ : $|\,2\,|+|\,-3\,| = 5 \lt |\,7\,|$
- $i=3$ : $|\,1\,|+|\,-1\,| = 2 \lt |\,4\,|$
となって $(\sharp)$ が成り立ちますのでヤコビ法は収束します。
証明 $\require{color}$
$\dps{
||\,H\,||_{\infty}
\mathop{=}^{\textcolor{red}{\mbox{Th.9}}}
\max_{1 \leqq i \leqq n} \left( \sum_{j=1}^n |\,h_{ij}\,| \right)
\mathop{=}^{\textcolor{red}{\mbox{Def.12}}}
\max_{1 \leqq i \leqq n} \left( \frac{1}{|\,a_{ii}\,|} \sum_{j \neq i} |\,a_{ij}\,| \right)
\ \mathop{\lt}^{\textcolor{red}{\Large{(\sharp)}}}\ 1.
}$
すると
Th.13 より
$\dps{
\left|\,\eee^{(k)}\,\right|_{\infty}
= \left|\,H^k \eee^{(0)}\,\right|_{\infty} \leqq \left(||\,H\,||_{\infty} \right)^k \left|\,\eee^{(0)}\,\right|
\longrightarrow 0 \qquad ( k \rightarrow \infty )
}$
よって
Cor.11 より
$\eee^{(k)}
\longrightarrow \ooo \qquad ( k \rightarrow \infty )$
$\therefore$
$\xxx^{(k)}
\longrightarrow \xxx \qquad ( k \rightarrow \infty )$
(証明終)
$\require{color}$
注と例
Rem.15 $(\sharp)$ が成り立てばガウス・ザイデル法も収束する。
証明は下記のプリント pp.4-5 を見てください。
Rem.16 $(\sharp)$ はあくまで十分条件であって、
$(\sharp)$ が成り立たなくてヤコビ法やガウス・ザイデル法が収束することも少なくありません。