前処理行列
線型代数、数値解析において、行列Aの前処理行列Pとは、P−1AがAより小さな条件数を持つ行列を指す。 前処理は、大規模疎行列を係数とする連立一次方程式
[math] Ax = b[/math]
を解くために反復法を用いる場合に有効である。これは、ほとんどの反復法で行列の条件数が増大するに従って収束率が低下するためである。具体的には、元の方程式を解く代わりに、左前処理を適用した方程式
[math] P^{-1}Ax = P^{-1}b,\,[/math]
すなわち
[math]c=P^{-1}b, \qquad (P^{-1}A)x=c[/math]
を解くか、もしくは右前処理を適用した方程式
[math] AP^{-1}Px = b,\,[/math]
すなわち
[math](AP^{-1})y=b, \qquad x=P^{-1}y[/math]
を解く。これらは、前処理行列Pが正則なら元の方程式と同値である。
これらの前処理の目的は、前処理を施した行列
[math]P^{-1}A[/math]
もしくは
[math]AP^{-1}[/math]
の条件数を小さくすることにある。
通常、Pの選択に関してはトレードオフがある。P-1は反復法の各ステップで適用する必要があるため、コストを抑えるためには計算しやすいものでなければならない。最も効率のよい前処理は
[math] P=I[/math]
もしくは
[math] P^{-1}=I[/math]
であるが、これは元の方程式と同じで前処理行列は何もしない。一方、
[math] P=A[/math]
すなわち
[math] \quad P^{-1}A = AP^{-1} = I[/math]
とすると条件数は最適な1となり、1回の反復で収束するが、
[math] P^{-1}=A^{-1}[/math]
の計算は元の方程式の求解と同程度に難しい。
そこで行列Pをこれらの中間から選び、P-1ができるだけ簡単に計算でき、かつ最小の反復回数となるように取る。
上の議論で、前処理行列
[math] P^{-1}A[/math]
もしくは
[math] AP^{-1}[/math]
は明に計算されないことに注意されたい。すなわち反復法では、与えられたベクトルに対する前処理の適用P-1だけが必要である。
また、Aが対称な場合、前処理の効果は、[math]P^{-1}A[/math]の固有値を互いに近づけることに相当する [1]。
Contents
例
以下の前処理では、Aを下三角部分L、対角部分D、上三角部分Uに分割する。
ヤコビ前処理
ヤコビ前処理は前処理の最も単純な形態の一つで、前処理行列
[math] P = D[/math]
は係数行列の対角要素のみからなる。つまり
[math]P_{ij} = A_{ii}\delta_{ij} = \begin{cases}A_{ii} & i=j \\ 0 &\mbox{otherwise}\end{cases}[/math]
である。
SOR
逐次過緩和（SOR）前処理では、Pを
[math]P=\left(\frac{D}{\omega}+L\right)\frac{\omega}{2-\omega}D^{-1}\left(\frac{D}{\omega}+U\right)[/math]
となるように取る。ωは[math]0 \lt \omega \lt 2[/math]を満たす緩和パラメータである。
Aが対称な場合を対称逐次過緩和前処理もしくはSSOR前処理という。
関連項目
外部リンク
- Preconditioned Conjugate Gradient – math-linux.com
- An Introduction to the Conjugate Gradient Method Without the Agonizing Pain by Jonathan Richard Shewchuk.
- www.preconditioners.com