練習問題メモ 12(行列の指数関数)

Contents

練習問題メモ 12(行列の指数関数)#

https://www2.itc.kansai-u.ac.jp/~afujioka/2012-2016/2012/la1/120627la1.pdf

行列指数関数

exp(A)=n=01n!An=I+A+A22!+

定理:exp(A+B)=exp(A)exp(B)

2つのn次正方行列A,Bについて、両者が可換ならば

exp(A+B)=exp(A)exp(B)

定理:exp(A)1=exp(A)

exp(A)exp(A)=exp(AA)=exp(O)=I

よって

exp(A)1=exp(A)

定理:exp(A)T=exp(AT)

k=1,2,について

(AT)k=(Ak)T

であるため

exp(A)T=exp(AT)

問1#

正方行列 (λ10λ) の指数関数を求めよ。

A=(λ10λ)

とする。

数学的帰納法により、

An=(λnnλn10λn)

であるため

証明

n=1のとき $A1=(λ11×λ00λ1)n = 2$ A^2 = \begin{pmatrix} \lambda & 1\ 0 & \lambda \end{pmatrix} \begin{pmatrix} \lambda & 1\ 0 & \lambda \end{pmatrix}

\begin{pmatrix} \lambda^2 & 2 \lambda\ 0 & \lambda^2 \end{pmatrix} $n=kAk=Ak1A=(λk1k1λk20λk1)(λ10λ)=(λkkλk10λk)n = k+1Ak+1=AkA=(λkkλk10λk)(λ10λ)=(λk+1(k+1)λk0λk+1)An=(λnnλn10λn)$

Aの指数関数eA

eA=n=01n!An=I+A+A22!+=(1001)+(λ10λ)+12!(λ22λ0λ2)+

(1,1)要素と(2,2)要素は

A11=A22=1+λ+λ22!+λ33!+=n=0xnn!=eλ

である。

(1,2)要素もまた

A12=0+1+2λ2!+3λ23!+=1+λ+λ22!+=n=0xnn!=eλ

よって

eA=(eλeλ0eλ)
行列指数関数

n次正方行列Xの指数関数eXは以下のように定義される

eX=n=01n!Xn

例:

A=(3004)

のような対角行列なら

Ak=(3k004k)

となるので、

eA=I+A+A22!+=(1001)+(3004)+12!(320042)+=(e300e4)
別解

https://www2.itc.kansai-u.ac.jp/~afujioka/2012-2016/2012/la1/120627la1.pdf

2つのn次正方行列A,Bについて、両者が可換ならば

exp(A+B)=exp(A)exp(B)

という定理を使う方法。

N=(0100)

とおいて、

(λ10λ)=λI+N

として、

exp(λ10λ)=exp(λI+N)=exp(λI)exp(N)(λIN)=(eλI)(I+N)(N2=O)=(eλeλ0eλ)

とする。

問2#

A が交代行列ならば、exp(A)は直交行列であることを示せ。

  • 交代行列:AT=A

  • 直交行列:AT=A1

Aが交代行列、すなわちAT=Aであれば

exp(A)T=exp(AT)((AB)T=BTAT(An)T=(AT)n)=exp(A)(A)=exp(A)1(exp(A)exp(A)=exp(AA)=Iexp(A)=exp(A)1)
(eA)T=0とはならないのか?
(eA)T=(n=01n!An)T=n=01n!(An)T=n=01n!(A)n=n=01n!(1)nAn=n{xN|x}1n!Ann{xN|x}1n!An

問3#

次の文章のを埋めよ。

A,Bn 次の正方行列とする。トレースについて、次の 13 が成り立つことを 示す。

  1. tr(AT)=trA

  2. tr(AB)=tr(BA)

  3. B が正則行列のとき、 tr(B1AB)=trA

  4. AT(i,i) 成分は A1 成分に一致するから、トレースの定義より、 1 が成り立つ。

  5. A,B(i,j) 成分をそれぞれ aij,bij とおくと、

tr(AB)=i=1nj=1n2=j=1ni=1nbjiaij=tr(3)

よって、2が成り立つ。

  1. B が正則行列ならば、 B の逆行列 B1 が存在し、

tr(B1AB)=tr(B1(AB))=tr(4)=tr(5)

よって、3が成り立つ。

  • 1:(i,i)

  • 2:aijbji

  • 3:BA

  • 4:(AB)B1

  • 5:A