可換な行列の定義・具体例・性質・計算問題について

はるか
はるか
可換な行列って何?
ふゅか
ふゅか
そうね。行列の積が交換できるやつね!例えば、逆行列と行列みたいな感じの!

1. 可換な行列とは

行列 A A B B が可換であるとは、行列の積が順序に依存せず同じ結果となることを意味します。以下の条件が満たされるとき、行列 A A B B は可換であると言います。

AB=BA AB = BA

1.1. 可換な行列の具体例

1.1.1. スカラー行列と任意の行列

スカラー行列 kI kI I I は単位行列)と任意の行列 A A は常に可換です。すなわち、次の式が成り立ちます。 kIA=AkI kI \cdot A = A \cdot kI これはスカラーの性質によるものです。

1.1.2. 対角行列同士

対角行列同士は常に可換です。例えば、 A=(a00b),B=(c00d) A = \begin{pmatrix} a & 0 \\ 0 & b \end{pmatrix}, \quad B = \begin{pmatrix} c & 0 \\ 0 & d \end{pmatrix} の場合、行列の積 AB AB BA BA は次のようになります。 AB=(ac00bd),BA=(ac00bd) AB = \begin{pmatrix} ac & 0 \\ 0 & bd \end{pmatrix}, \quad BA = \begin{pmatrix} ac & 0 \\ 0 & bd \end{pmatrix} したがって、AB=BA AB = BA となり、可換です。

1.1.3. 回転行列とその逆行列

2Dの回転行列 R(θ) R(\theta) とその逆行列 R(θ) R(-\theta) は可換です。例えば、 R(θ)=(cosθsinθsinθcosθ),R(θ)=(cosθsinθsinθcosθ) R(\theta) = \begin{pmatrix} \cos\theta & -\sin\theta \\ \sin\theta & \cos\theta \end{pmatrix}, \quad R(-\theta) = \begin{pmatrix} \cos\theta & \sin\theta \\ -\sin\theta & \cos\theta \end{pmatrix} これらの積は順序に依存せず、単位行列になります。 R(θ)R(θ)=R(θ)R(θ)=I R(\theta)R(-\theta) = R(-\theta)R(\theta) = I

はるか
はるか
なるほど。正則行列でもあるのか。

2. 可換な行列の性質

2.1. 二項定理が成り立つ

行列 A A B B が可換である場合、つまり AB=BA AB = BA であるとき、次の二項定理のような等式が成り立つ。

(A+B)n=k=0n(nk)AnkBk (A+B)^n = \sum_{k=0}^{n} \binom{n}{k} A^{n-k} B^k

ふゅか
ふゅか
数学的帰納法で示してみよう!

n=1 n = 1 のとき、 (A+B)1=A+B (A+B)^1 = A + B です。一方で、右辺を n=1 n = 1 で評価すると、 k=01(1k)A1kBk=(10)A10B0+(11)A11B1=A+B \sum_{k=0}^{1} \binom{1}{k} A^{1-k} B^k = \binom{1}{0} A^{1-0} B^0 + \binom{1}{1} A^{1-1} B^1 = A + B これにより、n=1のときは成り立ちます。

次に、n=m n = m のときに命題が成り立つと仮定します。

仮定により、

(A+B)m=k=0m(mk)AmkBk (A+B)^m = \sum_{k=0}^{m} \binom{m}{k} A^{m-k} B^k が成り立つとします。このとき、(A+B)m+1 (A+B)^{m+1} を次のように書けます。

(A+B)m+1=(A+B)(A+B)m (A+B)^{m+1} = (A+B)(A+B)^m

仮定を代入すると、

(A+B)m+1=(A+B)k=0m(mk)AmkBk (A+B)^{m+1} = (A+B)\sum_{k=0}^{m} \binom{m}{k} A^{m-k} B^k

右辺を分配法則を使って展開すると、

(A+B)m+1=k=0m(mk)Amk+1Bk+k=0m(mk)AmkBk+1 (A+B)^{m+1} = \sum_{k=0}^{m} \binom{m}{k} A^{m-k+1} B^k + \sum_{k=0}^{m} \binom{m}{k} A^{m-k} B^{k+1}

ここで、式を一つの総和にまとめるために、各項の k k を調整します。まず、左の項は k k のままで、右の項は k k k+1 k+1 に置き換えます。

(A+B)m+1=k=0m(mk)A(m+1)kBk+k=1m+1(mk1)Am(k1)Bk (A+B)^{m+1} = \sum_{k=0}^{m} \binom{m}{k} A^{(m+1)-k} B^k + \sum_{k=1}^{m+1} \binom{m}{k-1} A^{m-(k-1)} B^k

これをさらに整理すると、

(A+B)m+1=Am+1+k=1m((mk)+(mk1))A(m+1)kBk+Bm+1 (A+B)^{m+1} = A^{m+1} + \sum_{k=1}^{m} \left( \binom{m}{k} + \binom{m}{k-1} \right) A^{(m+1)-k} B^k + B^{m+1}

ここで、二項係数の性質 (mk)+(mk1)=(m+1k)\binom{m}{k} + \binom{m}{k-1} = \binom{m+1}{k} を使うと、

(A+B)m+1=Am+1+k=1m(m+1k)A(m+1)kBk+Bm+1 (A+B)^{m+1} = A^{m+1} + \sum_{k=1}^{m} \binom{m+1}{k} A^{(m+1)-k} B^k + B^{m+1}

これをまとめると、

(A+B)m+1=k=0m+1(m+1k)A(m+1)kBk (A+B)^{m+1} = \sum_{k=0}^{m+1} \binom{m+1}{k} A^{(m+1)-k} B^k

したがって、全ての自然数 n n に対して、

(A+B)n=k=0n(nk)AnkBk (A+B)^n = \sum_{k=0}^{n} \binom{n}{k} A^{n-k} B^k

が成り立つことが数学的帰納法により証明されました。

2.2. 指数関数

行列の指数関数 eA e^A eB e^B は、行列 A A B B が可換であれば次の性質を持ちます。

eA+B=eAeB=eBeA e^{A+B} = e^A e^B = e^B e^A

ただし、行列 A A の指数関数 eA e^A は、次の無限級数で定義されます。

eA=n=0Ann! e^A = \sum_{n=0}^{\infty} \frac{A^n}{n!}

行列 B B についても同様です。

行列 A A B B が可換である、すなわち AB=BA AB = BA が成立するとき、先ほど示したように二項定理が行列にも適用できます。次のように表せます。

(A+B)n=m=0nn!m!(nm)!AmBnm (A + B)^n = \sum_{m=0}^{n} \frac{n!}{m!(n-m)!} A^m B^{n-m}

eA+B e^{A+B} の定義に従って展開すると、

eA+B=n=0(A+B)nn! e^{A+B} = \sum_{n=0}^{\infty} \frac{(A+B)^n}{n!}

これに先ほどの二項定理を適用すると、

eA+B=n=01n!m=0nn!m!(nm)!AmBnm e^{A+B} = \sum_{n=0}^{\infty} \frac{1}{n!} \sum_{m=0}^{n} \frac{n!}{m!(n-m)!} A^m B^{n-m}

となります。この式を整理すると、

eA+B=n=0m=0nAmm!Bnm(nm)! e^{A+B} = \sum_{n=0}^{\infty} \sum_{m=0}^{n} \frac{A^m}{m!} \frac{B^{n-m}}{(n-m)!}

eA+B=(m=0Amm!)(n=0Bnn!) e^{A+B} = \left(\sum_{m=0}^{\infty} \frac{A^m}{m!}\right) \left(\sum_{n=0}^{\infty} \frac{B^n}{n!}\right)

したがって、

eA+B=eAeB e^{A+B} = e^A e^B

となります。

同様の議論を eB+A e^{B+A} についても行うと、eA+B=eB+A=eBeA e^{A+B} = e^{B+A} = e^B e^A も成立するため、

eA+B=eAeB=eBeA e^{A+B} = e^A e^B = e^B e^A

が示されました。

2.3. 対角行列同士が可換

A=diag(a1,a2,,an)A = \text{diag}(a_1, a_2, \dots, a_n)B=diag(b1,b2,,bn)B = \text{diag}(b_1, b_2, \dots, b_n)と置いたとき、A,BA,Bは可換である。

ここで、対角行列 AABB の要素は次のようになっています。

A=(a1000a2000an) A = \begin{pmatrix} a_1 & 0 & \cdots & 0 \\ 0 & a_2 & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & a_n \end{pmatrix} B=(b1000b2000bn) B = \begin{pmatrix} b_1 & 0 & \cdots & 0 \\ 0 & b_2 & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & b_n \end{pmatrix}

まず、行列 ABAB を計算します。

AB=(a1000a2000an)(b1000b2000bn) AB = \begin{pmatrix} a_1 & 0 & \cdots & 0 \\ 0 & a_2 & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & a_n \end{pmatrix} \begin{pmatrix} b_1 & 0 & \cdots & 0 \\ 0 & b_2 & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & b_n \end{pmatrix}

行列の積を行うと、対角成分はそれぞれの成分の積となります。

AB=(a1b1000a2b2000anbn) AB = \begin{pmatrix} a_1b_1 & 0 & \cdots & 0 \\ 0 & a_2b_2 & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & a_nb_n \end{pmatrix}

次に、行列 BABA を計算します。

BA=(b1000b2000bn)(a1000a2000an) BA = \begin{pmatrix} b_1 & 0 & \cdots & 0 \\ 0 & b_2 & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & b_n \end{pmatrix} \begin{pmatrix} a_1 & 0 & \cdots & 0 \\ 0 & a_2 & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & a_n \end{pmatrix}

同様に計算すると、

BA=(b1a1000b2a2000bnan) BA = \begin{pmatrix} b_1a_1 & 0 & \cdots & 0 \\ 0 & b_2a_2 & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & b_n a_n \end{pmatrix}

行列の積は成分ごとの積であり、aibi=biaia_i b_i = b_i a_i が成り立つため、 AB=BAAB = BA です。したがって、対角行列同士は可換であることが示されました。

3. 可換な行列の例題

ふゅか
ふゅか
可換な行列の性質を具体的に体感してみよう!

3.1. 例題1(二項定理が成り立つことを確認)

行列AA,BBが可換であるとき(A+B)m(A+B)^mm=24m=2\cdots 4のとき、二項定理が成り立つことを確認しなさい。

m=2 m = 2 の場合
(A+B)2=(A+B)(A+B) (A + B)^2 = (A + B)(A + B)

これを分配法則に従って展開すると、

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

=A2+AB+BA+B2 =A^2 + AB + BA + B^2

ここで AB=BA AB = BA なので、

(A+B)2=A2+2AB+B2 (A + B)^2 = A^2 + 2AB + B^2

m=3 m = 3 の場合

(A+B)3=(A+B)(A+B)(A+B) (A + B)^3 = (A + B)(A + B)(A + B)

=(A2+2AB+B2)(A+B) =(A^2 + 2AB + B^2)(A + B)
これを分配法則に従って展開すると、

=A2(A+B)+2AB(A+B)+B2(A+B) =A^2(A + B) + 2AB(A + B) + B^2(A + B)

さらに展開すると、

=A3+A2B+2ABA+2AB2+B2A+B3 =A^3 + A^2B + 2ABA + 2AB^2 + B^2A + B^3

ここで AB=BA AB = BA なので、整理して、

(A+B)3=A3+3A2B+3AB2+B3 (A + B)^3 = A^3 + 3A^2B + 3AB^2 + B^3

m=4 m = 4 の場合

(A+B)4=(A+B)(A+B)(A+B)(A+B) (A + B)^4 = (A + B)(A + B)(A + B)(A + B)

=(A3+3A2B+3AB2+B3)(A+B) =(A^3 + 3A^2B + 3AB^2 + B^3)(A + B)

これを展開すると、

=A4+A3B+3A2BA+3A2B2+3ABA2+3AB3+B3A+B4 =A^4 + A^3B + 3A^2BA + 3A^2B^2 + 3ABA^2 + 3AB^3 + B^3A + B^4

ここで AB=BA AB = BA を用いて整理すると、

(A+B)4=A4+4A3B+6A2B2+4AB3+B4 (A + B)^4 = A^4 + 4A^3B + 6A^2B^2 + 4AB^3 + B^4

3.2. 例題2(可換な行列に対する指数関数の性質)

次のような対角行列があるとします。

A=(1002),B=(3004) A = \begin{pmatrix} 1 & 0 \\ 0 & 2 \end{pmatrix}, \quad B = \begin{pmatrix} 3 & 0 \\ 0 & 4 \end{pmatrix}

eAe^A eBe^BeA+Be^{A+B}を計算しなさい。

まず、行列 A A B B の積を計算します。

AB=(1002)(3004)=(3008) AB = \begin{pmatrix} 1 & 0 \\ 0 & 2 \end{pmatrix} \begin{pmatrix} 3 & 0 \\ 0 & 4 \end{pmatrix} = \begin{pmatrix} 3 & 0 \\ 0 & 8 \end{pmatrix}

次に、行列 B B A A の積を計算します。

BA=(3004)(1002)=(3008) BA = \begin{pmatrix} 3 & 0 \\ 0 & 4 \end{pmatrix} \begin{pmatrix} 1 & 0 \\ 0 & 2 \end{pmatrix} = \begin{pmatrix} 3 & 0 \\ 0 & 8 \end{pmatrix}

AB=BA AB = BA であることが確認できます。したがって、これらの行列 A A B B は可換であると言えます。

ふゅか
ふゅか
次に、行列の指数関数 eA e^A eB e^B を計算してみましょう。

行列の指数関数は次のように定義されます。

eA=n=0Ann! e^A = \sum_{n=0}^{\infty} \frac{A^n}{n!}

対角行列に対する指数関数は、各対角成分に対する指数関数を計算することで簡単に求められます。したがって、

eA=(e100e2),eB=(e300e4) e^A = \begin{pmatrix} e^1 & 0 \\ 0 & e^2 \end{pmatrix}, \quad e^B = \begin{pmatrix} e^3 & 0 \\ 0 & e^4 \end{pmatrix}

次に、行列 A A B B の和 A+B A + B を求め、その指数関数を計算します。

A+B=(4006) A + B = \begin{pmatrix} 4 & 0 \\ 0 & 6 \end{pmatrix}

これにより、

eA+B=(e400e6) e^{A+B} = \begin{pmatrix} e^4 & 0 \\ 0 & e^6 \end{pmatrix}

次に、行列の積 eAeB e^A e^B を計算します。

eAeB=(e100e2)(e300e4)=(e1e300e2e4)=(e400e6) e^A e^B = \begin{pmatrix} e^1 & 0 \\ 0 & e^2 \end{pmatrix} \begin{pmatrix} e^3 & 0 \\ 0 & e^4 \end{pmatrix} = \begin{pmatrix} e^1 \cdot e^3 & 0 \\ 0 & e^2 \cdot e^4 \end{pmatrix} = \begin{pmatrix} e^4 & 0 \\ 0 & e^6 \end{pmatrix}

この結果から、eA+B=eAeB e^{A+B} = e^A e^B が成り立つことが確認できました。

PR