- 締切済み
行列について教えてください。
taka41の回答
- taka41
- ベストアンサー率27% (3/11)
基本的な線形代数の本で、 『線形代数とその応用』産業図書 G.ストラング著 ISBN 4-7828-0502-0 がわかりやすかったと思います。 ぼくも線型ニガテなんで。。。
関連するQ&A
- 行列の固有ベクトルの解法
現在行列の固有値と固有ベクトルをもとめるプログラムを作成しています。 手順としては、入力行列をハウスホルダー法により三重対角行列に変換し、その後ＱＲ法で対角化を行い固有値を求めます。 固有ベクトルはＬＵ分解を使用して固有値ごとに求めていこうと考えました。 現状固有値を求めるプログラムは作成できました（そして正しく求められていることも確認しました）。そして行列のLU分解を行うプログラムまで作成できたのですが、LU分解後の行列から固有ベクトルを求める方法がわかりません。 詳しく説明します Ax = λx を (A - Iλ)x = 0 として、この(A - Iλ)をLU分解しました。 すると式は LUx = 0 となり 最終的に Ux = 0 をとく問題になります。 ここで行列Uは上三角行列なので、１次の連立方程式を解くように、行列Uの右下の要素を使って計算を始めていくのですが、自分がなにか勘違いをしているのだと思うのですがこの方法で計算すると固有ベクトルが全て０になってしまいます。 行列Ｕ x ０ | 2 3 4 5 | |x1| ＝ |0| | 0 4 2 9 | |x2| ＝ |0| | 0 0 7 5 | |x3| ＝ |0| | 0 0 0 8 | |x4| ＝ |0| このような図式になり、固有ベクトルであるxを求めていくのですが、x4から順にもとめても０にしかならないんです。 下記のサイトを参考に学んでいたんですが、この部分が分からずにいます。 http://hooktail.org/computer/index.php?KL%C5%B8%B3%AB2 どこを勘違いしているんでしょうか？ アドバイスをお願いします。
- ベストアンサー
- 数学・算数
- 不完全ＬＵ分解前処理つき双共役勾配法についておしえてください。
連立方程式を解くために不完全ＬＵ分解前処理つき双共役勾配法 について勉強しています。 前処理の際に、行列Ａを不完全ＬＵ分解しその逆行列(LU)^(-1)というのを使用します。LU分解まではできたのですが、この逆行列は普通にＬＵ分解＋直接法という形でもとめるのでしょうか。だとしたら、直接法をつかっていてあまり高速化が期待できない様な気がしました。 不完全コレスキー分解つき共役勾配法(ICCG)のときは、不完全コレスキー分解後、間接的にAの逆行列をもとめて使用する方法がありましたのでなにかいい方法があるのかと思い質問しました。 はじめてのプログラミングで見当違いなことをいっているかもしれませんがよろしくおねがいします。
- ベストアンサー
- 数学・算数
- すみません。行列式のLU分解について質問なんですが。
すみません。行列式のLU分解について質問なんですが。 ある行列式AをLU分解するまではいいのですが、その後、どうやって固有値を求めたらいいのでしょうか？ QR分解で求める方法の方が効率が良いのはわかっているのですが、気になるので、教えてもらえたら幸いです。 あと、あまり数学が得意でないので、噛み砕いて説明していただけるとほんと助かります。
- ベストアンサー
- 数学・算数
- 行列演算： 固有ベクトルの解法
現在、対称行列の固有値、固有ベクトルを求めるプログラムを作成し、つい最近完成しました。 しかし、とても使い物にならないプログラムになってしまいました。 理由はとても遅いのです。 解法の手順として、まず固有値を求めてから固有ベクトルを求めるようと考え、入力の対称行列をHouseHolder法により三重対角行列に変換し、それをQR法により対角化してまず固有値を求めました。 固有値を求めることができたので、次に固有ベクトルを求めます。手順として、固有値ごとに入力対称行列の対角成分から固有値を減算した行列をＬＵ分解し、連立一次方程式を解くように固有ベクトルを求めていきます。 この一連の手順で、対称行列の固有値、固有ベクトルを求めることができたのですが、とても時間がかかってしまいます。 ただし、対称行列の固有値を求めるまでの時間はとても高速です。 ５００×５００の行列の固有値、固有ベクトルを求めるのに３０分はかかってしまいますが、その中で固有値を求める時間は２秒しかかかりません。 つまり今固有値がわかっている状態で、固有ベクトルを高速に求めたいと考えています。 なにか高速に固有ベクトルを求める方法（アルゴリズム）はあるでしょうか？
- ベストアンサー
- 数学・算数
- コレスキー分解
現在、数学で固有値問題などを勉強しています。様々な参考書を借りて勉強しているのですが具体的な計算例がほとんど載っていないため理解できません。コレスキー分解というものでつまずいていますので質問します。対称行列においてCC'=Aのように上三角行列の転置行列である下三角行列C'とCの積に分解することをコレスキー分解というそうです。 例 4 2 2 A= 2 5 3 という行列をコレスキー分解すると 2 3 11 2 0 0 2 1 1 C'= 1 2 0 C= 0 2 1 となるのはどのような計算をしている 1 1 3 0 0 3 のですか。 CかC'のどちらかが求められればもう一方が求まるぐらいしか分かりません。方針でも構いませんので教えてください、お願いします。
- ベストアンサー
- 数学・算数
- （非）対称行列の実固有値をすべて求めたい
1000次元程度の行列の実固有値を数値的にすべて求めたいのですが、 いい方法はあるでしょうか？ 行列に特に対称性は仮定しません。 通常はQR法+ハウスホルダー変換を使うと思います。 しかし、以下のページ http://www.ma.is.saga-u.ac.jp/minamoto/book/book8/program/chap9/program94.c のコードで実行して見ましたが、虚数解を含むような場合には 計算が終了しないようです。
- ベストアンサー
- 数学・算数
- 行列の積の固有値
「nxn正則行列A,B,X=ABに関して，AとBの固有値と固有ベクトルが分かっているときに，Xの固有値と固有ベクトルを求めよ．」 という問題が解けません．どなたかお分かりになる方いらっしゃいますでしょうか？ 補足 難しさは変わらないと思うのですが，Aは上三角行列，Bは下三角行列，X,A,Bは確率行列（全ての行を足すと(1,1,...,1)となる）という条件がありました． ですので，X,A,Bはそれぞれ固有値１を持つというのは分かります．この条件を使わないほうがいいのですが，使った場合でもありがたいです． なお，反復法などを用いて数値計算的に解きたいのではなくて，解析的にビシッと解きたいので，よろしくお願いします． ヤフー知恵袋のほうでも質問させていただいたのですが，有効な解が得られれませんでした． 教科書や，問題集の練習問題などではなく私が個人的に解きたい問題です．
- 締切済み
- 数学・算数