ヤコビ法のソースを表示

{{出典の明記|date=2022年3月}}
{{distinguish|ヤコビ法 (固有値問題)}}
'''ヤコビ法'''（ヤコビほう）とは<math>n</math>元の[[連立一次方程式]]<math>A\vec{x}=\vec{b}</math>を[[反復法_(数値計算)|反復法]]で解く手法の1つである。[[ドイツ]]の[[数学者]][[カール・グスタフ・ヤコブ・ヤコビ]]の名前にちなむ。

<math>n</math>次[[正方行列]]<math>A</math>は、上[[三角行列]]<math>U</math>、下三角行列<math>L</math>、[[対角行列]]を<math>D</math>とすると、''<math>A=L+D+U</math>''と書ける。このようにすると、まず以下のような変形ができる。

<center><math>
\begin{array}{ccc}
(L+D+U) \vec{x} &=& \vec{b} \\
D\vec{x} &=& \vec{b}-(L+U)\vec{x} \\
\end{array}
</math></center>

この式を満たす''<math>\ x</math>''を求める。初期値<math>\vec{x}^{(0)}</math>に対して、
<math>\ k</math>回目の反復で得られた<math>x_1</math>の値を<math>x_1^{(k)}</math>と書くと、
以下のような反復法の[[漸化式]]ができる。
<center><math>
D\vec{x}^{(k+1)} = \vec{b}-(L+U)\vec{x}^{(k)}
</math></center>
この式は以下のように変形できる。
<center><math>
\vec{x}^{(k+1)} = D^{-1} \{ \vec{b}-(L+U)\vec{x}^{(k)} \}
</math></center>
もし、解が収束した場合、その場合は<math>x_1^{(k+1)}</math>と<math>x_1^{(k)}</math>は共通の値<math>x_1^{(*)}</math>を持つことになる。このとき、
<center><math>
\vec{x}^{(*)} = D^{-1} \{ \vec{b}-(L+U)\vec{x}^{(*)} \}
</math></center>
となり、変形していくと元の連立方程式の形に戻る。
したがって、ヤコビ法で解が収束した場合、その解は連立方程式の解となる。
また、その収束の[[十分条件]]は、係数行列の対角要素の絶対値が非対角要素の絶対値よりも相対的に大きい場合、すなわち対角優位な行列である場合に収束する。これは[[ガウス＝ザイデル法]]も同様である。

ヤコビ法の式は[[幾何ベクトル|ベクトル]]<math>\vec{x}</math>の各成分ごとに次のような式で書くことができ、[[数値解析]]ではこの式が用いられる。
<center><math>
x_{i}^{(k+1)} = \frac{1}{a_{ii}} \left( b_{i} - \sum_{j=1}^{i-1} a_{ij}x_{j}^{(k)} - \sum_{j=i+1}^{n} a_{ij}x_{j}^{(k)} \right) = \frac{1}{a_{ii}} \left( b_{i} - \sum_{j \neq i}^{n} a_{ij}x_{j}^{(k)} \right)
</math></center>

ガウス＝ザイデル法とヤコビ法を加速する方法としては[[SOR法]]が知られている。

== 具体例 ==
3元の連立一次方程式、すなわち、

<math>\left( \begin{array}{ccc} a_{11} & a_{12} & a_{13} \\ a_{21} & a_{22} & a_{23} \\ a_{31} & a_{32} & a_{33} \end{array} \right) \left( \begin{array}{c} x_1 \\ x_2 \\ x_3 \end{array} \right) = \left( \begin{array}{c} b_1 \\ b_2 \\ b_3 \end{array} \right)</math>

を解くことを考える。<math>k</math>回目の反復で得られた<math>x_1</math>の値を<math>x_1^{(k)}</math>と書く。
初期値<math>\vec{x}^{(0)}</math>は、適当な値、例えばゼロベクトルでもかまわない。

<math>x_1^{(k+1)} = (b_1 - a_{12} x_2^{(k)} - a_{13} x_3^{(k)}) / a_{11}</math>

<math>x_2^{(k+1)} = (b_2 - a_{21} x_1^{(k)} - a_{23} x_3^{(k)}) / a_{22}</math>

<math>x_3^{(k+1)} = (b_3 - a_{31} x_1^{(k)} - a_{32} x_2^{(k)}) / a_{33}</math>

という反復を繰り返していく。
ヤコビ法は、直列計算ではガウス＝ザイデル法よりも遅いが、ガウス＝ザイデル法と異なり各式が他の式に依存せず並列性があるため[[並列コンピューティング|並列計算]]でも用いられる。

== 固有値問題 ==
実対称行列の固有値および固有ベクトルを求める繰り返し計算手法においても'''ヤコビ法'''と呼ばれる解法がある（紛らわしさを避けるためにはヤコビ対角化法という）。
<math>\ n</math>次の実[[対称行列]]<math>\ A</math>について次のように <math>\ G(p, q, \theta)</math> による相似変換、すなわち[[ギブンス回転]]を実行することにより、非対角要素 <math>\ a_{ij}(i \ne j)</math> の最大値 <math>\ a_{pq}</math> が0となるようにする。

:<math>\ B=G^TAG</math>

これによって行列 <math>\ B</math> の各要素は次のようになる。但し、 <math>\ i,j \ne p,q </math> である。

:<math>\begin{cases} b_{pp} = a_{pp}\cos^2\theta + a_{qq}\sin^2\theta - 2a_{pq}\sin\theta\cos\theta ,\\ 
                     b_{qq} = a_{pp}\sin^2\theta + a_{qq}\cos^2\theta + 2a_{pq}\sin\theta\cos\theta ,\\
                     b_{pq} = b_{qp} = \frac{1}{2}(a_{pp}-a_{qq})\sin2\theta + a_{pq}\cos2\theta ,\\
                     b_{ij} = a_{ij} ,\\
                     b_{pj} = a_{pj}\cos\theta - a_{qj}\sin\theta ,\\
                     b_{qj} = a_{pj}\sin\theta + a_{qj}\cos\theta ,\\
                     b_{ip} = a_{ip}\cos\theta - a_{iq}\sin\theta ,\\
                     b_{iq} = a_{ip}\sin\theta + a_{iq}\cos\theta
\end{cases}</math>

ここで、<math>\ a_{pq} \ne 0</math> のとき <math>\ b_{pq}=0</math> となる <math>\ \theta</math> は上式より

:<math> \tan 2\theta = \frac{-2a_{pq}}{(a_{pp}-a_{qq})}</math>

から求められることがわかる。
ギブンス回転をすべての非対角要素がほぼ0になるまで繰り返せば、実対称行列 <math>A\quad</math> が対角化された形となるから、その対角要素が <math>A\quad</math> の固有値となる。また、 <math>A\quad</math> がk回変換された行列を <math>A_{k}\quad</math> 、k回目のギブンス回転を表す直交行列を <math>G_{k}\quad</math> と表せば、

:<math> A_k = G_k^T A_{k-1} G_k = U_k^TAU_k</math>
:ここに、 <math> U_k = G_1 G_2 G_3 \cdots G_{k-1} G_k</math>

となる。 <math>A_{k}\quad</math> のすべての非対角要素がほぼ0となったとき、 <math>U_{k}\quad</math> は固有ベクトルを並べた行列となっている。
なお、ギブンス回転の繰り返し過程において、一度は0になった要素がその後の変換により0でなくなることもあるが、変換の繰り返しによって非対角項は0に近づいてゆく。

なお上記のように、ヤコビの対角化法は実対称（あるいは複素エルミート）の場合が最も良く知られておりその場合だけしか適用できないものと考えられるかもしれない。しかし非対称な行列に対するヤコビ法もあって研究もされプログラムも公開されていたがQR法が登場した後では行列が非対称な場合にはほとんどQR法だけが使われているため今日では非対称な場合についてはほとんど認知されていないようである（複素エルミートの場合についてもその具体的な実装に言及している文献は稀であり，もっぱら実対称行列の場合だけがよく知られている）。

== 関連項目 ==
*[[数値解析]]
*[[ガウス＝ザイデル法]]

{{Linear algebra}}

{{デフォルトソート:やこひほう}}
[[Category:数値線形代数]]
[[Category:緩和法 (反復法)]]
[[Category:数学に関する記事]]
[[Category:数学のエポニム]]