- ベストアンサー
行列のベキ乗
最小多項式が重解をもったり、固有値が虚数になる場合などの対角化不可能な整数の行列（たとえば２×２）のベキ乗の一般項を、手計算で求める方法はあるのでしょうか？ 数値によっては、ケイリーハミルトンなどで次数を下げたときに、等差数列などに帰着できれば可能かと思いますが、、、。
- sunasearch
- お礼率62% (82/132)
- 数学・算数
- 回答数5
- ありがとう数5
- みんなの回答 （5）
- 専門家の回答
質問者が選んだベストアンサー
いずれにせよn=2など特別の場合でなければ相当困難だと思います。困難というのはがんばれば不可能ではない、というぐらいの意味です。たとえジョルダン標準形を求めたとしても、ジョルダン行列のk乗はもはやジョルダン行列ではありませんので、n次正方行列の場合は、ベキ零行列をn乗以下した行列と対角行列のべき乗が複雑に入り組んだものが出てきます。この議論はある種の行列の収束の議論の役には立つこともありますが、直接n乗計算するのにはあまり向いていないような気もします。対角成分だけは綺麗ですけれど。 ケーリーハミルトンで次数下げというのはよい方法です。nが具体的に与えられた（あまり大きくない整数なら）役に立つと思います。 n=2（もっと一般にも無理ではないですが、複雑なので割愛）の場合は標準形のn乗は比較的簡単に計算できます。（以下少し行列がずれて見にくいかも） [a_n c_n] [0 b_n] が行列Aのn乗だとしましょう。Aは [a 1] [0 b] であるとします。A^{n+1}は [a・a_n a_n+b・c_n] [0. b・b_n ] となりますから(1,2)成分を計算すればよいですが、それは漸化式 c_{n+1}=a_n+b・c_n をとくことに帰着されます。a_n=a^nに注意すれば、この両辺をa^{n+1}で割ってやって c_{n+1}/a^{n+1}=(b/a)・(c_n/a^n)+1/a となります。新たにd_n=c_n/a^nとおいてやればこれは線型二項間漸化式 d_{n+1}=(b/a)d_n+1/a のことです。初項d_1=c_1/a^1=1/aとわかっていますから、これは計算可能です。あとは標準化した操作を逆にしてn乗の公式を得ます。 行列のランクが上がっていくと連立漸化式を解くことになりそうですが、うまく線形化できたとして、そして一つの数列に関して解けばn項間漸化式のようなものがでてきます。そもそもn次方程式を代数的に解くのがほとんど困難な場合しかない以上、一般には解くのが不可能な場合がほとんどなのではないかと個人的には思います。2×2じゃないときはコンピュータに任せるのが賢明なのではないかと。
その他の回答 (4)
- keyguy
- ベストアンサー率28% (135/469)
任意の正方行列Aは適当な正則行列Uとジョルダン標準形Jによって A=U・J・U^(-1) となります A^n=U・J^n・U^(-1) です J1,J2,…を適当なジョルダン細胞とすると J=diag(J1,J2,…) です そして J^n=diag(J1^n,J2^n,…) なおジョルダン細胞のn乗は2項定理の係数がでてきれいにもとまります [a100] [0a10] [00a1] [000a] を2乗,3乗,4乗,… していくと規則性が見えます
お礼
回答ありがとうございます。 一般的なｎ乗の形について、理解できたような気がします。 手で計算するのは難しそうですね。
- onakyuu
- ベストアンサー率45% (36/80)
>対角化不可能な場合は、ジョルダン化の後にｎ乗を >計算することができるようになるということでしょ >うか? そのとおりです。 ジョルダン標準形について調べてみるとよいでしょう。
- keyguy
- ベストアンサー率28% (135/469)
おっしゃるとおり要素が整数の行列： [0.5 1.0] [0.0 0.5] は要素が整数でない行列ですがこれを含まないと言うことならば私の回答はなかったことにしてください
- keyguy
- ベストアンサー率28% (135/469)
質問の意味が難解なので回答者がいません たとえば整数の行列とは要素が整数の行列のことを言うのですがそうでない意味で使っているような気がします 対角化できなくてもジョルダン化はできます ジョルダン細胞のｎ乗はきれいな式になるでしょう
お礼
回答ありがとうございます。 整数の行列とは、おっしゃるとおり要素が整数の行列です。 ジョルダン化は知らないのですが、対角化不可能な場合は、ジョルダン化の後にｎ乗を計算することができるようになるということでしょうか？
関連するQ&A
- 行列のベキ乗（固有値が虚数の時）
先に、 http://oshiete1.goo.ne.jp/kotaeru.php3?q=1137580 で質問させていただいた者です。 ２x２の行列で、固有値が実数のときには、ジョルダン化を用いてベキ乗の計算ができるようになりました。 が、固有値が虚数のときにも同様の方法で解けるのでしょうか？ 実部と虚部に分けて計算したり、虚数平面を用いたりすることで計算できるのでしょうか？
- ベストアンサー
- 数学・算数
- 線形代数 ケーリー・ハミルトンの定理 詳しい参考書
※線形代数の参考書についての質問はかなり出ていますので，既出でしたらもすいません．（探し出すことができませんでした．） 現在，制御の研究をしているものですが，その中で，ケーリー・ハミルトンの定理を使って行列計算をしなければなりません． この定理は，A^nの次数を次数n-1以下の行列で表現するために用います． 今，A行列の固有値が重解を持つため，不可能がと思っていましたが，重解をもつ場合の方法をネットで調べると，同様の方法でできるという内容がすぐに出てきました．（http://www.cfv21.com/math/hamcayley.htm） できれば，参考書や論文などで，確認したいので，ケーリー・ハミルトンの定理について詳しく記されているものがあれば紹介いただけたらと思います．よろしくお願いします．
- 締切済み
- 数学・算数
- 行列固有値問題
Aは、3×3行列で、3つの固有値のうち2つが同じ(1組が重解)で、もう一つが異なる解、つまり固有値λ1、λ2、λ3で λ1=λ2 λ3≠λ1 の場合、 Aが対称行列ではないもの具体例を示して下さい。また、その具体例の行列を対角化する行列Pも示して下さい。 この時、求める最小多項式は重解はないものとします。 つまり、(A-λ1E)(A-λ3E)=0 をみたし、 対角化した行列は、λ1=λ2、λ1≠λ3で [λ1 0 0] [0 λ2 0] [0 0 λ3] になります。 このようなAでなおかつ対称行列でないものをあげて欲しいのですが、存在しますか? 対称行列だったら、いくつか列があったのですが、そうでない具体例が知りたいのです。
- 締切済み
- 数学・算数
- 行列・対角化可能の条件は？
行列で対角化可能の時の条件を教えて下さい。 問題で固有値、固有ベクトル、対角化可能の場合は対角化する正則行列を求めよ、とあります。 ３×３行列で固有値が３つ、全て異なる場合は対角化可能。 固有値が１つ（３重解）の場合は対角化不可。 では、固有値が２つの場合は対角化可能と不可の場合がありますが、これはどのようにして見分けるのでしょうか？ 例えば －３ －２ －２ B＝[ ２ １ ２ ] ２ ２ １ の時、固有値は１、－１（重解）ですが対角化可能です。なぜでしょうか？宜しくお願いします。
- ベストアンサー
- 数学・算数
- 行列の対角化
┌1 －2 －2┐ A＝│1 2 2│ └(-2) 2 1┘ という行列なのですが、対角化できるのでしょうか？ 何度も何度も解きなおしてるんですけど対角化できません。 Aの固有方程式の解で重解になっているものがないので対角化は・・可能ですよね？ 固有値として-1、±√7が求まるのですが、±√7に対する固有空間を考えるとどうしても固有ベクトルとして成分がすべて0の（3,1）行列しか出てこなく、対角化行列が ┌0 0 0┐ P＝│1 0 0│ └(-1) 0 0┘ といったような行列になってしまうのですが、この場合P^(-1)が存在しないためP^(-1)*A*Pは存在しない事になり、Aは対角化不可能ということになってしまいますよね？？ 多分どこか間違った理解をしているところがあると思います。 どなたかご教授お願いできないでしょうか？
- 締切済み
- 数学・算数
- ケーリーハミルトンの定理を用いる問題
A=[ 2 1] [-7 -3] とする。このAに対する固有多項式をgA(t)とおく。(Aは、2×2の行列を表している) ケーリーハミルトンの定理を用いて、f(A)を計算せよ。 f(t)=t^20 固有多項式は、gA(t)={t^2}+t+1なので、ケーリーハミルトンの定理から、単位行列をEとおくと、gA(A)={A^2}+A+E=０,A^3=Eであることを用いると、答えには書いてあるのですが、A^3＝Eということがどうして出てくるのかが、分かりません。 どなたかよろしくお願いします。
- ベストアンサー
- 数学・算数
- 行列の問題です、よろしくお願いします。
A=(a b ; c d) （←2次の正方行列をこのように表すとします) (1)行列Aが固有値λ1、λ2 (λ1≠0、λ2≠0)を持つとするとき、ケーリー・ハミルトンの式を用いて、 tr(A)=λ1+λ2、 det(A)=λ1・λ2 となることを示せ。 (λ1の"1"などはλの添え字だとします) (2)上記の条件の下で、(λ1-λ2)・A＾n=((λ1)＾n-(λ2)＾n)・A-((λ1)＾n・λ2-λ1・(λ2)＾n)・E が成り立つことを示せ。ただしnは正の整数とする。 ( "＾n" はn乗を、"E"は単位行列を表しています) という問題がよくわかりません。 (1)は、僕なりの解としては、ケーリー・ハミルトンを用いなければ A-λE=(a-λ b ; c d-λ) det(A-λE)=(a-λ)(d-λ)-bc =λ＾2-(a+d)λ+ad-bc=0 この方程式の2解はλ1、λ2なので、解と係数の関係より λ1+λ2=a+d=tr(B) λ1・λ2=ad-bc=det(B) としましたが、ケーリー・ハミルトンを用いるとどのようになるのでしょうか？ (2)は、全然方針が思い浮かびません…どのように解くのでしょうか？よろしくお願いします。
- 締切済み
- 物理学
お礼
詳細な回答と、具体的なご説明、 ありがとうございます。 ２ｘ２の場合にはできるということですね。 ジョルダン化を勉強して試してみたいと思います。