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;
戻る