記事ソース/ジョルダン細胞のn乗
をテンプレートにして作成
査読
rst2hooktail
進行表
執筆中
かぎマニュ
物理のかぎプロジェクト
トップ
最近の更新
ヘルプ
開始行:
#rst2hooktail_source
=========================================================...
ジョルダン細胞のn乗
=========================================================...
ジョルダンの標準形で有名なジョルダン細胞のn乗を求めます。
ジョルダン細胞
=====================
ジョルダン細胞とは、次のk次正方行列のことを言います。
<tex>
J_k =
\begin{pmatrix}
\lambda & 1 & 0 & \cdots & 0 \\
0 & \lambda & 1 & \cdots & 0 \\
0 & 0 & \lambda & \cdots & 0 \\
\vdots & \vdots & \vdots & \ddots & \vdots \\
0 & 0 & 0 & \cdots & \lambda
\end{pmatrix} \tag{##}
</tex>
このn乗を求めてみましょう。
注目する性質は、対角行列(単位行列の定数倍) $\Lambda$
<tex>
\Lambda =
\begin{pmatrix}
\lambda & 0 & 0 & \cdots & 0 \\
0 & \lambda & 0 & \cdots & 0 \\
0 & 0 & \lambda & \cdots & 0 \\
\vdots & \vdots & \vdots & \ddots & \vdots \\
0 & 0 & 0 & \cdots & \lambda
\end{pmatrix} \tag{##}
</tex>
のどんな行列とも可換な性質と、
べきゼロ行列 $N$
<tex>
N =
\begin{pmatrix}
0 & 1 & 0 & \cdots & 0 \\
0 & 0 & 1 & \cdots & 0 \\
0 & 0 & 0 & \cdots & 0 \\
\vdots & \vdots & \vdots & \ddots & \vdots \\
0 & 0 & 0 & \cdots & 0
\end{pmatrix} \tag{##}
</tex>
の持つ、何乗かするとゼロになる性質です。
ためしに四次のべきゼロ行列のべき乗を求めてみましょう。
<tex>
N =
\begin{pmatrix}
0 & 1 & 0 & 0 \\
0 & 0 & 1 & 0 \\
0 & 0 & 0 & 1 \\
0 & 0 & 0 & 0
\end{pmatrix} \tag{##}
</tex>
<tex>
N^2 =
\begin{pmatrix}
0 & 0 & 1 & 0 \\
0 & 0 & 0 & 1 \\
0 & 0 & 0 & 0 \\
0 & 0 & 0 & 0
\end{pmatrix} \tag{##}
</tex>
<tex>
N^3 =
\begin{pmatrix}
0 & 0 & 0 & 1 \\
0 & 0 & 0 & 0 \\
0 & 0 & 0 & 0 \\
0 & 0 & 0 & 0
\end{pmatrix} \tag{##}
</tex>
<tex>
N^4 =
\begin{pmatrix}
0 & 0 & 0 & 0 \\
0 & 0 & 0 & 0 \\
0 & 0 & 0 & 0 \\
0 & 0 & 0 & 0
\end{pmatrix} \tag{##}
</tex>
とこの様に、べき乗すると1のなすラインが上がっていきます。
ここで求めたいのは、 $J_k$ の $n$ 乗、
<tex>
(J_k)^n = (\Lambda + N)^n \tag{##}
</tex>
です。二項定理を用います。
<tex>
(J_k)^n &= (\Lambda + N)^n \\
&= _n C_0 \Lambda^n + _n C_1 \Lambda^{n-1} N^{1} + _n C_2...
</tex>
ここで $N$ のべき数を昇順にならべました。あるところからは...
簡単な例
================
<tex>
J_4 =
\begin{pmatrix}
\lambda & 1 & 0 & 0 \\
0 & \lambda & 1 & 0 \\
0 & 0 & \lambda & 1 \\
0 & 0 & 0 & \lambda
\end{pmatrix} \tag{##}
</tex>
<tex>
(J_4)^2 =
\begin{pmatrix}
\lambda^2 & 2 \lambda & 1 & 0 \\
0 & \lambda^2 & 2 \lambda & 1 \\
0 & 0 & \lambda^2 & 2 \lambda \\
0 & 0 & 0 & \lambda^2
\end{pmatrix} \tag{##}
</tex>
<tex>
(J_4)^3 =
\begin{pmatrix}
\lambda^3 & 3 \lambda^2 & 3 \lambda & 1 \\
0 & \lambda^3 & 3 \lambda^2 & 3 \lambda \\
0 & 0 & \lambda^3 & 3 \lambda^2 \\
0 & 0 & 0 & \lambda^3
\end{pmatrix} \tag{##}
</tex>
<tex>
(J_4)^4 =
\begin{pmatrix}
\lambda^4 & 4 \lambda^3 & 6 \lambda^2 & 4 \lambda \\
0 & \lambda^4 & 4 \lambda^3 & 6 \lambda^2 \\
0 & 0 & \lambda^4 & 4 \lambda^3 \\
0 & 0 & 0 & \lambda^4
\end{pmatrix} \tag{##}
</tex>
とこの様に簡単にべき乗が求まります。
行列の指数関数
=======================
行列の指数関数がジョルダン細胞の場合にも、求まったので書...
<tex>
\exp (t J_k) &\equiv \sum_{n=0}^\infty \dfrac{(t J_k)^n}{...
&= \sum_{n=0}^\infty \dfrac{1}{n!} \sum_{i=0}^n \dfrac{n!...
&= \sum_{i=0}^\infty N^i \dfrac{t^i}{i!} \sum_{n=i}^\inft...
&= \sum_{i=0}^\infty N^i \dfrac{t^i}{i!} \exp (t \lambda)...
&= \sum_{i=0}^{k-1} N^i \dfrac{t^i}{i!} \exp (t \lambda) ...
</tex>
よって、例えば、四次なら、
<tex>
\exp (t J_4)
=
\begin{pmatrix}
\exp (t \lambda) & \dfrac{t}{1!} \exp(t \lambda) & \dfra...
0 & \exp(t \lambda) & \dfrac{t}{1!} \exp(t \lambda) & \...
0 & 0 & \exp(t \lambda) & \dfrac{t}{1!} \exp(t \lambda)...
0 & 0 & 0 & \exp(t \lambda)
\end{pmatrix} \tag{##}
</tex>
となります。以上でこの話は終わりです。 続々ベクトルの回転...
ここまで読んだなら、その応用をぜひ知ってください。 ジョル...
今日はここまで、お疲れ様でした。
.. _続々ベクトルの回転: http://hooktail.sub.jp/vectoranal...
.. _ジョルダン標準形の指数関数の応用: http://hooktail.sub...
@@author:クロメル@@
@@accept:2013-03-11@@
@@category:物理数学@@
@@id:jordanCalcu@@
終了行:
#rst2hooktail_source
=========================================================...
ジョルダン細胞のn乗
=========================================================...
ジョルダンの標準形で有名なジョルダン細胞のn乗を求めます。
ジョルダン細胞
=====================
ジョルダン細胞とは、次のk次正方行列のことを言います。
<tex>
J_k =
\begin{pmatrix}
\lambda & 1 & 0 & \cdots & 0 \\
0 & \lambda & 1 & \cdots & 0 \\
0 & 0 & \lambda & \cdots & 0 \\
\vdots & \vdots & \vdots & \ddots & \vdots \\
0 & 0 & 0 & \cdots & \lambda
\end{pmatrix} \tag{##}
</tex>
このn乗を求めてみましょう。
注目する性質は、対角行列(単位行列の定数倍) $\Lambda$
<tex>
\Lambda =
\begin{pmatrix}
\lambda & 0 & 0 & \cdots & 0 \\
0 & \lambda & 0 & \cdots & 0 \\
0 & 0 & \lambda & \cdots & 0 \\
\vdots & \vdots & \vdots & \ddots & \vdots \\
0 & 0 & 0 & \cdots & \lambda
\end{pmatrix} \tag{##}
</tex>
のどんな行列とも可換な性質と、
べきゼロ行列 $N$
<tex>
N =
\begin{pmatrix}
0 & 1 & 0 & \cdots & 0 \\
0 & 0 & 1 & \cdots & 0 \\
0 & 0 & 0 & \cdots & 0 \\
\vdots & \vdots & \vdots & \ddots & \vdots \\
0 & 0 & 0 & \cdots & 0
\end{pmatrix} \tag{##}
</tex>
の持つ、何乗かするとゼロになる性質です。
ためしに四次のべきゼロ行列のべき乗を求めてみましょう。
<tex>
N =
\begin{pmatrix}
0 & 1 & 0 & 0 \\
0 & 0 & 1 & 0 \\
0 & 0 & 0 & 1 \\
0 & 0 & 0 & 0
\end{pmatrix} \tag{##}
</tex>
<tex>
N^2 =
\begin{pmatrix}
0 & 0 & 1 & 0 \\
0 & 0 & 0 & 1 \\
0 & 0 & 0 & 0 \\
0 & 0 & 0 & 0
\end{pmatrix} \tag{##}
</tex>
<tex>
N^3 =
\begin{pmatrix}
0 & 0 & 0 & 1 \\
0 & 0 & 0 & 0 \\
0 & 0 & 0 & 0 \\
0 & 0 & 0 & 0
\end{pmatrix} \tag{##}
</tex>
<tex>
N^4 =
\begin{pmatrix}
0 & 0 & 0 & 0 \\
0 & 0 & 0 & 0 \\
0 & 0 & 0 & 0 \\
0 & 0 & 0 & 0
\end{pmatrix} \tag{##}
</tex>
とこの様に、べき乗すると1のなすラインが上がっていきます。
ここで求めたいのは、 $J_k$ の $n$ 乗、
<tex>
(J_k)^n = (\Lambda + N)^n \tag{##}
</tex>
です。二項定理を用います。
<tex>
(J_k)^n &= (\Lambda + N)^n \\
&= _n C_0 \Lambda^n + _n C_1 \Lambda^{n-1} N^{1} + _n C_2...
</tex>
ここで $N$ のべき数を昇順にならべました。あるところからは...
簡単な例
================
<tex>
J_4 =
\begin{pmatrix}
\lambda & 1 & 0 & 0 \\
0 & \lambda & 1 & 0 \\
0 & 0 & \lambda & 1 \\
0 & 0 & 0 & \lambda
\end{pmatrix} \tag{##}
</tex>
<tex>
(J_4)^2 =
\begin{pmatrix}
\lambda^2 & 2 \lambda & 1 & 0 \\
0 & \lambda^2 & 2 \lambda & 1 \\
0 & 0 & \lambda^2 & 2 \lambda \\
0 & 0 & 0 & \lambda^2
\end{pmatrix} \tag{##}
</tex>
<tex>
(J_4)^3 =
\begin{pmatrix}
\lambda^3 & 3 \lambda^2 & 3 \lambda & 1 \\
0 & \lambda^3 & 3 \lambda^2 & 3 \lambda \\
0 & 0 & \lambda^3 & 3 \lambda^2 \\
0 & 0 & 0 & \lambda^3
\end{pmatrix} \tag{##}
</tex>
<tex>
(J_4)^4 =
\begin{pmatrix}
\lambda^4 & 4 \lambda^3 & 6 \lambda^2 & 4 \lambda \\
0 & \lambda^4 & 4 \lambda^3 & 6 \lambda^2 \\
0 & 0 & \lambda^4 & 4 \lambda^3 \\
0 & 0 & 0 & \lambda^4
\end{pmatrix} \tag{##}
</tex>
とこの様に簡単にべき乗が求まります。
行列の指数関数
=======================
行列の指数関数がジョルダン細胞の場合にも、求まったので書...
<tex>
\exp (t J_k) &\equiv \sum_{n=0}^\infty \dfrac{(t J_k)^n}{...
&= \sum_{n=0}^\infty \dfrac{1}{n!} \sum_{i=0}^n \dfrac{n!...
&= \sum_{i=0}^\infty N^i \dfrac{t^i}{i!} \sum_{n=i}^\inft...
&= \sum_{i=0}^\infty N^i \dfrac{t^i}{i!} \exp (t \lambda)...
&= \sum_{i=0}^{k-1} N^i \dfrac{t^i}{i!} \exp (t \lambda) ...
</tex>
よって、例えば、四次なら、
<tex>
\exp (t J_4)
=
\begin{pmatrix}
\exp (t \lambda) & \dfrac{t}{1!} \exp(t \lambda) & \dfra...
0 & \exp(t \lambda) & \dfrac{t}{1!} \exp(t \lambda) & \...
0 & 0 & \exp(t \lambda) & \dfrac{t}{1!} \exp(t \lambda)...
0 & 0 & 0 & \exp(t \lambda)
\end{pmatrix} \tag{##}
</tex>
となります。以上でこの話は終わりです。 続々ベクトルの回転...
ここまで読んだなら、その応用をぜひ知ってください。 ジョル...
今日はここまで、お疲れ様でした。
.. _続々ベクトルの回転: http://hooktail.sub.jp/vectoranal...
.. _ジョルダン標準形の指数関数の応用: http://hooktail.sub...
@@author:クロメル@@
@@accept:2013-03-11@@
@@category:物理数学@@
@@id:jordanCalcu@@
ページ名:
Modified by
物理のかぎプロジェクト
PukiWiki 1.4.6
Copyright © 2001-2005
PukiWiki Developers Team
. License is
GPL
.
Based on "PukiWiki" 1.3 by
yu-ji
Powered by PHP 5.3.29 HTML convert time to 0.002 sec.