1000次元程度の行列の実固有値を数値的にすべて求めたいのですが、
いい方法はあるでしょうか?
行列に特に対称性は仮定しません。
通常はQR法+ハウスホルダー変換を使うと思います。
しかし、以下のページ
http://www.ma.is.saga-u.ac.jp/minamoto/book/book …
のコードで実行して見ましたが、虚数解を含むような場合には
計算が終了しないようです。
No.1ベストアンサー
- 回答日時:
にいろいろ、解法が書いてあるみたいです。
ただ、質問者さんの引用されてあるプログラムは
C言語による
「数値計算入門」
~ 解法・アルゴリズム・プログラム ~
http://www.saiensu.co.jp/?page=book_details&ISBN …
のものようですが、これは、対称行列でない場合にも、適用出できるのでしょうか?
もともとの原理は、対称行列の場合は固有値は実数で、固有ベクトルが基底になる(つまり直交行列で対角化可能)
ですから、一番大きな固有値でわった行列をかけることで他の固有ベクトルは0に収束するというのが出発点のようですが。
一般の行列は、対角化できるとは限らない、複素行列を使わないといけない。とかあるので、固有値は共役複素数をつかって実数化する云々は、前述の論文には書いてありますが。
この回答への補足
ありがとうございます。
早速読んでみます。
>対称行列でない場合にも、適用出できるのでしょうか?
固有値が全部実数の場合には
事前にハウスホルダー変換することにより
適用できるようです。
念の為、
http://keisan.casio.jp/has10/SpecExec.cgi?path=0 …
と比較して見ましたが、特に問題ありませんでした。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・人生のプチ美学を教えてください!!
- ・10秒目をつむったら…
- ・あなたの習慣について教えてください!!
- ・牛、豚、鶏、どれか一つ食べられなくなるとしたら?
- ・【大喜利】【投稿~9/18】 おとぎ話『桃太郎』の知られざるエピソード
- ・街中で見かけて「グッときた人」の思い出
- ・「一気に最後まで読んだ」本、教えて下さい!
- ・幼稚園時代「何組」でしたか?
- ・激凹みから立ち直る方法
- ・1つだけ過去を変えられるとしたら?
- ・【あるあるbot連動企画】あるあるbotに投稿したけど採用されなかったあるある募集
- ・【あるあるbot連動企画】フォロワー20万人のアカウントであなたのあるあるを披露してみませんか?
- ・映画のエンドロール観る派?観ない派?
- ・海外旅行から帰ってきたら、まず何を食べる?
- ・誕生日にもらった意外なもの
- ・天使と悪魔選手権
- ・ちょっと先の未来クイズ第2問
- ・【大喜利】【投稿~9/7】 ロボットの住む世界で流行ってる罰ゲームとは?
- ・推しミネラルウォーターはありますか?
- ・都道府県穴埋めゲーム
- ・この人頭いいなと思ったエピソード
- ・準・究極の選択
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
行列 線形代数 "diag"って何...
-
正値ってなんですか
-
行列の読み方
-
行列の性質に、二つの行列A,Bが...
-
行列の三角化について
-
次元とランク
-
1、)A×Aの転置行列=Aの転置行...
-
一般的な行列の逆行列に関する質問
-
反対称行列の行列式
-
対角成分が等しい対称行列の正...
-
可換で対角化可能な2つの行列...
-
【線形代数】 基本変形によって...
-
3x3行列の固有値と固有ベクト...
-
固有値、固有ベクトルおよび対...
-
合成変換とケーリーハミルトン
-
Mn(R)の次元
-
2の6乗の答えと計算方法
-
電卓の使い方 乗数はどうした...
-
e^(x^2)の積分に関して
-
積分で1/x^2 はどうなるのでし...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
おすすめ情報