Three Dimensional Finite Element Method
eigenvaluelanczos
-10

行列が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