行列がn×nの場合の計算は、どうなるのでしょうか。計算手順としては、まず、 非対角要素で絶対値が最大の要素を探します。例えば、[A]の非対角でaklが絶対値で一番大きかったとします。 この場合の[P]行列は、以下になります。
上の[P]を用いて[P]-1[A][P]=[A']を計算すると、以下が得られます。
上式に最後の式をa'kl=0と置くことにより、回転角度θが得られます。式で 書くと以下のようになります。
また、 固有ベクトルの初期の行列として以下のように設定しておきます。
BACK | NEXT |
---|
Menu | LU Decompo | Stiff | 3D Solid | 3D Fluid | Eigen&Lanczos | Sound Eigen | Solid Eigen | Solid Axisym |