Ddisna Dopgtr Dopmtr Dorgtr Dormtr Dpteqr Dsbtrd Dsptrd Dstebz Dstedc Dstein Dstemr Dsteqr Dsterf Dsytrd

Dsteqr() Sub Dsteqr ( Compz As  String, N As  Long, D() As  Double, E() As  Double, Z() As  Double, Info As  Long  ) 対称三重対角行列の固有値・固有ベクトル (QL/QR法) 目的 本ルーチンは対称行列のすべての固有値, および, 必要により固有ベクトルを Implicit QL または QR法により求める. Dsytrd, Dsptrd または Dsbtrd を使用して三重対角形に変換した場合は, 元の対称行列または対称帯行列の固有ベクトルを求めることもできる. 引数 [in] Compz = "N": 固有値のみ求める. = "V": 元の対称行列の固有値・固有ベクトルを求める. 元の行列から三重対角形への変換を行った直交行列を Z() に入れておかなければならない. = "I": 三重対角行列の固有値・固有ベクトルを求める. Z() は単位行列に初期化される. [in] N 行列の行および列数. (N >= 0) (N = 0 の場合, 処理を行わずに戻る) [in,out] D() 配列 D(LD - 1) (LD >= N) [in] 三重対角行列の対角要素. [out] Info = 0 であれば固有値が昇順に入る. [in,out] E() 配列 E(LE - 1) (LE >= N - 1) [in] 三重対角行列の N - 1 個の副対角要素. [out] E() は破壊される. [in,out] Z() 配列 Z(LZ1 - 1, LZ2 - 1) (LZ1 >= N, LZ2 >= N) [in] Compz = "V" の場合, 三重対角形への変換を行った直交行列を Z() に入れる. [out] Info = 0 であれば, Compz = "V" の場合, 元の対称行列の正規直交化した固有ベクトルが Z() に入る. Compz = "I" の場合, 対称三重対角行列の正規直交化した固有ベクトルが Z() に入る. Compz = "N" の場合, Z() は参照されない. [out] Info = 0: 正常終了. = -1: パラメータ Compz の誤り. (Compz <> "N", "V" および "I") = -2: パラメータ N の誤り. (N < 0) = -3: パラメータ D() の誤り. = -4: パラメータ E() の誤り. = -5: パラメータ Z() の誤り. = i > 0: 30*N 回以内の反復ですべての固有値を見つけることはできなかった. E() の i 個の要素が 0 に収束しなかった. 終了時, D() と E() には元の行列に直交相似な対称三重対角行列の要素が入る. 出典 LAPACK 使用例 対称行列 A の固有値・固有ベクトルを求める. ただし, ( 2.20 -0.11 -0.32 ) A = ( -0.11 2.93 0.81 ) ( -0.32 0.81 2.37 ) とする. Dsytrdで対称三重対角形に変換したのち, DorgtrとDsteqrを適用する. Sub Ex_Dsytrd_Dsteqr() Const N = 3 Dim A(N - 1, N - 1) As Double, Info As Long Dim D(N - 1) As Double, E(N - 2) As Double, Tau(N - 2) As Double A(0, 0) = 2.2 A(1, 0) = -0.11: A(1, 1) = 2.93 A(2, 0) = -0.32: A(2, 1) = 0.81: A(2, 2) = 2.37 Call Dsytrd("L", N, A(), D(), E(), Tau(), Info) If Info <> 0 Then Debug.Print "Error in Dsytrd: Info =", Info Exit Sub End If Call Dorgtr("L", N, A(), Tau(), Info) If Info <> 0 Then Debug.Print "Error in Dorgtr: Info =", Info Exit Sub End If Call Dsteqr("V", N, D(), E(), A(), Info) If Info <> 0 Then Debug.Print "Error in Dsterf: Info =", Info Exit Sub End If Debug.Print "Eigenvalues =", D(0), D(1), D(2) Debug.Print "Eigenvectors =" Debug.Print A(0, 0), A(0, 1), A(0, 2) Debug.Print A(1, 0), A(1, 1), A(1, 2) Debug.Print A(2, 0), A(2, 1), A(2, 2) End Sub 実行結果 Eigenvalues = 1.70705954911046 2.22943643244226 3.56350401844728 Eigenvectors = 0.399322077213382 0.894521385341403 0.200931256445799 -0.481026444340548 0.390994588677271 -0.78468897753835 0.780484105216166 -0.216690384632057 -0.586421212707156