「行列（線形）の収束について」

Question

行列P(n×ｎ行列)は、全ての固有値の絶対値が１より小さければ、Pの累乗は零行列に収束する。
（実際この文では行列Pの成分は全てマイナスなしの、０より大きく１より小さい値としている）

　ｒ→∞　　⇒　　P^r（Pのｒ乗）→０（零行列）

とある文にさらっと書いてあったのですが、これをどう説明（証明）すればいいか困っています。たしかに固有値が全て１より小さい（分数）本などを調べてもわからずで、、アドバイスでもいいのでいただけると幸いです。ジョルダンの標準形を使うといいのでは？という意見もありましたが、はたしてどうなのでしょうか？

ringohatimitu · Accepted Answer

作用素環の一般論（今の場合は行列なので有限次元ヒルベルト空間上の有界線形作用素となります）から従うスペクトル半径の公式：
max{|λ|：λはPのスペクトル}=lim(P^nの作用素ノルム)^{1/n} (n→∞)

今の場合は行列なのでスペクトルは固有値全体と一致し、固有値の絶対値がすべて１より小さいので上の公式における極限は１より小さくなります。したがってP^nの作用素ノルムｎ→∞のときに０に収束することが明らかでこれよりP^nが行列として０に収束することが言えます。　
ちなみに上の公式は比較的簡単に得られるもので関数解析や作用素論あるいは作用素環論の大抵の入門書に載ってると思うので証明など参照してみてください。

toranekosan · Answer

＃５の回答は間違いです。あほな事書いたなぁ…^^;

toranekosan · Answer

対角化が出来なくても、上の対角成分付近に０でない要素を寄せ集める行列Uが存在し、

100 010 001
010 001 000
001 000 000

というような、行列で展開できるかと思います。

それが大体N個できるわけで、
E(1)(これは対角成分のみ）
E(2)(これは対角成分は０でその一こ上の成分だけが１)
という風にして
ある行列Uの元に
UPU(-1)=Σa(i)E(i)＝A
と成りますね。当然a(i)は一より小さいです。
E(i)*E(j)=(i=1あるいはj=1ならば、E(j)あるいはE(i)それいがいは０行列）であるので、これを利用して
A^(r)の展開を考えてやれば…いいような気がしないでも
ないです。ただちゃんとやってないのでどうなのでしょう。こういうアプローチもありかと思います。
かなりいい加減です…。

ringohatimitu · Answer

すいません、下に書いた回答においては暗に行列Pが対角化可能であることを仮定してしまっていました。一般の場合にこのような議論がうまくいくか分かり次第また回答したいと思います。混乱させてしまって申し訳ありません。

ringohatimitu · Answer

少し行列計算から離れてみますと：
R^nにおける任意のベクトルｘはker(P)と各固有空間の直和成分の和で書けるのでそれをｘ＝Σa(j)e(j)と表します、ここでe(j)は核空間と各固有値α(j)に対応する固有空間からとった直交基底、すなわちPe(j)=α(j)e(j)を満たしているとします。
このとき|x|≦1に対して
|<P^r(x),x>|≦Σ|a(j)|^2 |α(j)|^r≦Σ|α(j)|^r
となります。ここで和はｊに関して1からnまでで固定されているのでこれによりPの作用素ノルムはr→∞のとき０に収束します。作用素ノルムとmaxノルム（成分の絶対値の最大値）は行列に対して同値なのでr→∞のとき各成分が一斉に０に収束することが示せます。

adinat · Answer

対角化できる場合は♯1さんの方法が分かりよいですが、一般に対角化できるとは限らないので、その場合はジョルダン標準形を使うのが定石です。ジョルダン標準形の対角成分はすべて0に収束するのは明らかですが、非対角成分に1を含みますので、r乗したときの数字は少しややこしいですが、二項展開を用いて適切なorder評価をしてやることによっていくらでも小さくなることが示せます。正直少し面倒くさいですが、ジョルダン標準形の応用として誰でも一度はやる問題だとも思いますので、じっくり計算されてみては、と思います。

toranekosan · Answer

まずは対角化します。
UPU^(-1)=A
Aは対角行列です。

A^nは、対角成分以外は０です。
nが大きくなると、対角成分の値（つまり固有値）は
１より小さいので０に収束します。

A^r(r→∞)=O
A^r=UPU^(-1)*UPU^(-1)…UPU^(-1)=UP^(r)U^(-1)
なので、
P^r=U^(-1)A^(r)U
となり、
P^r(r→∞)=U^(-)A^r(r→∞)U=U0U=0

となるのではないでしょうか・・・

「行列（線形）の収束について」

作用素環の一般論（今の場合は行列なので有限次元ヒルベルト空間上の有界線形作用素となります）から従うスペクトル半径の公式：

＃５の回答は間違いです。

対角化が出来なくても、上の対角成分付近に０でない要素を寄せ集める行列Uが存在し、

すいません、下に書いた回答においては暗に行列Pが対角化可能であることを仮定してしまっていました。

少し行列計算から離れてみますと：

対角化できる場合は♯1さんの方法が分かりよいですが、一般に対角化できるとは限らないので、その場合はジョルダン標準形を使うのが定石です。

まずは対角化します。

関連するカテゴリからQ&Aを探す

デイリーランキングこのカテゴリの人気デイリーQ&Aランキング

マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング