A~をハウスホルダ変換によって得られるAと相似な行列とする。A~の固有ベクトルXkは
以下の計算を繰り返せば近似的に求められる。
 Xk = (μI - A~)-1Xk
ただし、
  μはXkに対応する固有値μkの任意の近似値である。

説明すると長くなるので省略するが、計算は以下のように3つのステップに分けられる。

 ステップ1：(これはA~を左下三角行列と右上三角行列の積に変換するステップである)
 for(i=1; k < n; k++)
       {
       mk = βk / αk
       αk+1 = αk+1 - βk×mk
       }
 

 ステップ2： 
 for(i=2; i < n; i++)
       Xi = Xi - mi-1×Xi-1;
 Xn = Xn/αn;
 
 ステップ3：
 for(i=n-1; i > 0; i--)
       Xi = (Xi-βi×Xi+1)/αi;

 

戻る