No.1ベストアンサー
- 回答日時:
n×n行列Aの固有値問題。
困り度3ということですから、中途半端ながらres付けてしまいます。
ヤコビ法は大きいnには手間が大きくて辛いでしょうね。
QR法が便利だけれど、一手間掛けまして、原点を移動する。つまり狙いの固有値を絶対値≒0になるように座標変換しておいてからQR法を使うのが旨いやり方だと思います。つまり予想される値μを使って、Aの代わりにA-μI の固有値を求めて、得られた答えにμを加えればよい。
反復計算で固有値と固有ベクトルを求める段階では、固有値の絶対値同士がどのぐらい散らばっているかによって計算の手間が全然違う。絶対値の順に固有値を並べた時、隣り合う固有値の比が1に近いと手間が掛かります。固有値が一個ずつ出てくるのを使って問題を減次するのも、nが大きいときには結構な手間です。なかなか一概には論じきれないな。
また固有値と固有ベクトルの近似値μ,vを改良するという仕上げは欠かせません。つまり
(A-μI)x = v
を解いてxを求め、改良された固有ベクトルの近似値xと、固有値の近似値(Ax)[j]/x[j]を求める。これを繰り返します。(これは最初の近似が良ければ凄く速く収束します)。QR法はこの計算にも向いてます。
スパースな行列(つまり零要素が殆ど、という場合)だと、また話が違うようです。
いずれにせよ、誤差の累積が怖いので、計算の精度(有効桁数)を変えて2度計算し、両者の結果が一致することを確認すべきです。もし食い違うようなら、有効桁数が足りない訳です。
No.2
- 回答日時:
大きい固有値だけに限って求められれば十分、という応用が沢山あります。
また、制御工学では不安定な固有値成分だけ取り出したい。そういう場合どうするのか、ちょっと調べてみたら「ランチョス法に、等角写像を組み合わせる。」
というテクニックがあるみたいです。
等角写像
y = (x+c) / (x-c) (x,yは複素数、cは適当な正の定数)
によって、不安定なやつを単位円の外に、安定な奴を単位円の中に写像しておいて、ランチョス法(絶対値が大きい固有値から順に出てきます)を適用するらしい。
まだ調査中ですが、お急ぎのようなのでとりあえず。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・人生のプチ美学を教えてください!!
- ・10秒目をつむったら…
- ・あなたの習慣について教えてください!!
- ・牛、豚、鶏、どれか一つ食べられなくなるとしたら?
- ・【大喜利】【投稿~9/18】 おとぎ話『桃太郎』の知られざるエピソード
- ・街中で見かけて「グッときた人」の思い出
- ・「一気に最後まで読んだ」本、教えて下さい!
- ・幼稚園時代「何組」でしたか?
- ・激凹みから立ち直る方法
- ・1つだけ過去を変えられるとしたら?
- ・【あるあるbot連動企画】あるあるbotに投稿したけど採用されなかったあるある募集
- ・【あるあるbot連動企画】フォロワー20万人のアカウントであなたのあるあるを披露してみませんか?
- ・映画のエンドロール観る派?観ない派?
- ・海外旅行から帰ってきたら、まず何を食べる?
- ・誕生日にもらった意外なもの
- ・天使と悪魔選手権
- ・ちょっと先の未来クイズ第2問
- ・【大喜利】【投稿~9/7】 ロボットの住む世界で流行ってる罰ゲームとは?
- ・推しミネラルウォーターはありますか?
- ・都道府県穴埋めゲーム
- ・この人頭いいなと思ったエピソード
- ・準・究極の選択
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
固有値の値について
-
固有値の求める順番?
-
行列の2分の1乗の計算の仕方
-
A・B=B・AならばAの固有...
-
固有値が複素数のときの固有ベ...
-
行列A^2からAを求める方法について
-
縮退行列のときのジョルダン形...
-
文字(定数)入りの行列の固有値
-
固有値の問題に関して
-
近似固有値と固有値の違い
-
極値の判定でヘッシアンの値が...
-
固有ベクトルが複数の場合
-
行列A^2の固有値は、Aの各固有...
-
この行列の固有値が全て実数と...
-
固有値が出ません…。
-
線形代数12
-
行列の固有値
-
数学について質問です。 nを正...
-
パソコンで行列はどう書けばいい?
-
sinθ―√3cosθ=a(θ+α)の形にした...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
おすすめ情報