13/02/21 12:16:43.83 .net
皆遊んでいる様なので質問を・・・
固有値、固有ベクトルの計算をする時の固定条件の入れ方を教えて頂けないでしょうか？
ググっているのですがそれらしい文献を見つけられずに困っています。
詳しく書けと言われそうなので少し書いておきますが長いので面倒なら読み飛ばしてください
固有値、固有ベクトルの計算プログラムを作っているのですが、今までは固有ベクトルがぐちゃぐちゃだったのを
逆反復法で計算することでそれっぽい(と言ってもまだ怪しい)ものが出るようになったのですがoctaveで
プログラムしているのでとても遅いです。そこで比較的規模が大きくても良さそうなサブスペース法で組むことに
したのですがどうも小行列を作るとそれがコレスキー分解出来ないのです境界条件を設定していなく、剛体モード
を含んでいるのが原因ではないかなと思っています。その他にヤコビ法で小行列の固有値、固有ベクトルを求めると
サブスペース法の繰り返し計算が収束しない様で困っています。エラーが出ない行列を手作業で与えてやると
ビルトインのeigの結果と比較して符号が違う固有ベクトルがでていました。ハウスホルダー＋QR法＋逆反復法で
出した結果はeigと比べてパターンは違いますが符号が違う固有ベクトルが出ていますがサブスペース法は収束します
ここらもエスパーレスしてくれると児童買春・児童ポルノ法違反有難いです