行列の加法・減法・実数倍・乗法 + R言語
行列同士の足し算と引き算
行列の足し算と引き算は同じ型(同じサイズとも言う)の行列同士で可能。同じ型の行列同士を重ねてちょうど同じ位置にある成分同士を足したり引いたりすればいい。
\[ B = \begin{pmatrix} 5 & -7 \\ -2 & 8 \end{pmatrix} \]
-3 | 5 |
10 | -6 |
5 | -7 |
-2 | 8 |
R言語ではmatrix関数を使って行列を作る。第1引数に行列の要素を、第2引数に行数を、第3引数に列数を指定する。
> A = matrix(c(-3, 10, 5, -6), 2, 2) > A [,1] [,2] [1,] -3 5 [2,] 10 -6 > B = matrix(c(5, -2, -7, 8), 2, 2) > B [,1] [,2] [1,] 5 -7 [2,] -2 8
行列Aと行列Bはともに2x2行列なので足し算や引き算が可能であることが分かる。
A + B
&= \begin{pmatrix} -3 + 5 & 5 + (-7) \\ 10 + (-2) & -6 + 8 \end{pmatrix} \\
&= \begin{pmatrix} 2 & -2 \\ 8 & 2 \end{pmatrix}
\end{align} \]
-3 + 5 | 5 + (-7) |
10 + (-2) | -6 + 8 |
計算結果が正しいかどうかR言語で確認。
> A + B [,1] [,2] [1,] 2 -2 [2,] 8 2
A - B
&= \begin{pmatrix} -3 - 5 & 5 - (-7) \\ 10 - (-2) & -6 - 8 \end{pmatrix} \\
&= \begin{pmatrix} -8 & 12 \\ 12 & -14 \end{pmatrix}
\end{align} \]
-3 - 5 | 5 - (-7) |
10 - (-2) | -6 - 8 |
計算結果が正しいかどうかR言語で確認。
> A - B [,1] [,2] [1,] -8 12 [2,] 12 -14
ちなみに、行列の足し算については交換法則と結合法則とが成り立つ。
- 交換法則:
- \( A + B = B + A \)
- 結合法則:
- \( (A + B) + C = A + (B + C) \) (Cも行列)
行列と実数との掛け算(実数倍)
行列のすべての成分を一定の実数で倍加することは実数倍と呼ばれている。要するに実数と行列との掛け算。実数をrで表わしてそれを2だと仮定すると次のように計算することができる。
\[ \begin{align} rA
&= \begin{pmatrix} 2 \times -3 & 2 \times 5 \\ 2 \times 10 & 2 \times -6 \end{pmatrix} \\
&= \begin{pmatrix} -6 & 10 \\ 20 & -12 \end{pmatrix}
\end{align} \]
\( 2 \times -3 \) | \( 2 \times 5 \) |
\( 2 \times 10 \) | \( 2 \times -6 \) |
計算結果が正しいかどうかR言語で確認。
> 2 * A [,1] [,2] [1,] -6 10 [2,] 20 -12
&= \begin{pmatrix} 2 \times 5 & 2 \times -7 \\ 2 \times -2 & 2 \times 8 \end{pmatrix} \\
&= \begin{pmatrix} 10 & -14 \\ -4 & 16 \end{pmatrix}
\end{align} \]
\( 2 \times 5 \) | \( 2 \times -7 \) |
\( 2 \times -2 \) | \( 2 \times 8 \) |
計算結果が正しいかどうかR言語で確認。
> 2 * B [,1] [,2] [1,] 10 -14 [2,] -4 16
行列の実数倍においても交換法則が成り立つ。
\[ rB = Br \]
R言語で交換法則を確認。
> r = 3 > r * A [,1] [,2] [1,] -9 15 [2,] 30 -18 > A * r [,1] [,2] [1,] -9 15 [2,] 30 -18
行列の実数倍においても結合法則が成り立つ。
\[ r = 2,\ s = 3 \]
\[ \begin{align} r(sA)
&= 2 \times \begin{pmatrix} 3 \times -3 & 3 \times 5 \\ 3 \times10 & 3 \times -6 \end{pmatrix} \\
&= 2 \times \begin{pmatrix} -9 & 15 \\ 30 & -18 \end{pmatrix} \\
&= \begin{pmatrix} 2 \times -9 & 2 \times 15 \\ 2 \times 30 & 2 \times -18 \end{pmatrix} \\
&= \begin{pmatrix} -18 & 30 \\ 50 & -36 \end{pmatrix}
\end{align} \]
\[ \begin{align} (rs)A
&= (2 \times 3) \times \begin{pmatrix} -3 & 5 \\ 10 & -6 \end{pmatrix} \\
&= 6 \times \begin{pmatrix} -3 & 5 \\ 10 & -6 \end{pmatrix} \\
&= \begin{pmatrix} 6 \times -3 & 6 \times 5 \\ 6 \times 10 & 6 \times -6 \end{pmatrix} \\
&= \begin{pmatrix} -18 & 30 \\ 60 & -36 \end{pmatrix}
\end{align} \]
結合法則が成り立つかどうかR言語で確認。
> r = 2; s = 3 > r * (s * A) [,1] [,2] [1,] -18 30 [2,] 60 -36 > (r * s) * A [,1] [,2] [1,] -18 30 [2,] 60 -36
行列の実数倍においては分配法則も成り立つ。
\[ r = 2 \]
\[ \begin{align} r(A + B)
&= 2 \times \begin{pmatrix} -3 + 5 & 5 + (-7) \\ 10 + (-2) & -6 + 8 \end{pmatrix} \\
&= 2 \times \begin{pmatrix} 2 & -2 \\ 8 & 2 \end{pmatrix} \\
&= \begin{pmatrix} 2 \times 2 & 2 \times -2 \\ 2 \times 8 & 2 \times 2 \end{pmatrix} \\
&= \begin{pmatrix} 4 & -4 \\ 16 & 4 \end{pmatrix}
\end{align} \]
\[ \begin{align} rA + rB
&= \begin{pmatrix} 2 \times -3 & 2 \times 5 \\ 2 \times 10 & 2 \times -6 \end{pmatrix} + rB \\
&= \begin{pmatrix} -6 & 10 \\ 20 & -12 \end{pmatrix} + rB \\
&= rA + \begin{pmatrix} 2 \times 5 & 2 \times -7 \\ 2 \times -2 & 2 \times 8 \end{pmatrix} \\
&= rA + \begin{pmatrix} 10 & -14 \\ -4 & 16 \end{pmatrix} \\
&= \begin{pmatrix} -6 + 10 & 10 + (-14) \\ 20 + (-4) & -12 + 16 \end{pmatrix} \\
&= \begin{pmatrix} 4 & -4 \\ 16 & 4 \end{pmatrix}
\end{align} \]
次のような分配法則も成り立つ。
\[ r = 2,\ s = 3 \]
\[ \begin{align} (2 + 3)A
&= 5 \times \begin{pmatrix} -3 & 5 \\ 10 & -6 \end{pmatrix} \\
&= \begin{pmatrix} 5 \times -3 & 5 \times 5 \\ 5 \times 10 & 5 \times -6 \end{pmatrix} \\
&= \begin{pmatrix} -15 & 25 \\ 50 & -30 \end{pmatrix}
\end{align} \]
\[ \begin{align} 2A + 3A
&= \begin{pmatrix} 2 \times -3 & 2 \times 5 \\ 2 \times 10 & 2 \times -6 \end{pmatrix} + 3A \\
&= 2A + \begin{pmatrix} 3 \times -3 & 3 \times 5 \\ 3 \times 10 & 3 \times -6 \end{pmatrix} \\
&= \begin{pmatrix} -6 & 10 \\ 20 & -12 \end{pmatrix} + \begin{pmatrix} -9 & 15 \\ 30 & -18 \end{pmatrix} \\
&= \begin{pmatrix} -6 + (-9) & 10 + 15 \\ 20 + 30 & -12 + (-18) \end{pmatrix} \\
&= \begin{pmatrix} -15 & 25 \\ 50 & -30 \end{pmatrix}
\end{align} \]
分配法則が成り立つかどうかR言語で確認。
> r = 2 > r * (A + B) [,1] [,2] [1,] 4 -4 [2,] 16 4 > r * A + r * B [,1] [,2] [1,] 4 -4 [2,] 16 4 > s = 3 > (r + s) * A [,1] [,2] [1,] -15 25 [2,] 50 -30 > r * A + s * A [,1] [,2] [1,] -15 25 [2,] 50 -30
- 交換法則:
- \( rA = Ar \)
- 結合法則:
- \( (rs)A = r(sA) \)
- 分配法則:
- \( r(A + B) = rA + rB \)
- \( (r + s)A = rA + sA \)
要するに、行列の和と差と実数倍においては実数同士のそれらと同様の法則が成り立つし、実数同士の計算のように扱えるということ。移項や簡約といった操作も可能。
行列同士の掛け算
行列同士の掛け算は実数の掛け算とはそのやり方が異なる。行列同士の足し算と引き算や実数倍のようには計算することができない。また、すべての行列同士が掛け算できるわけではない。行列Aと行列Bがあったとき、行列Aの列数と行列Bの行数が等しい場合にのみ行列同士の掛け算が成り立つ。例えば、3x2行列と2x2行列の掛け算は成り立つが、2x2行列と3x2行列の掛け算は成り立たない。
手始めに、まずは1行しか持たない行ベクトルと1列しか持たない列ベクトルとの掛け算から。なぜなら行ベクトルと列ベクトルの掛け算が行列同士の掛け算の基本になるから。次の例は、1x3行列すなわち3次の行ベクトルと3x1行列すなわち3次の列ベクトルの掛け算。前者の列数3と後者の行数3が一致するので行列同士の掛け算が成り立つ。
\begin{pmatrix} 3 & 8 & 5 \end{pmatrix} \times \begin{pmatrix} 6 \\ 2 \\ 4 \end{pmatrix}
&= 3 \times 6 + 8 \times 2 + 5 \times 4 \\
&= 18 + 16 + 20 \\
&= 10 + 10 + 20 + 8 + 6 \\
&= 40 + 14 \\
&= 54 \end{align} \]
6 | |||
2 | |||
4 | |||
3 | 8 | 5 | 54 |
R言語で計算結果が正しいかどうか確認。行列同士の掛け算の演算子には*ではなく%*%を使うことに要注意。*を使うと、R言語では行列の足し算や引き算と同じように各要素を単純に掛けた値になってしまう。
> A = matrix(c(3, 8, 5), 1, 3) > A [,1] [,2] [,3] [1,] 3 8 5 > B = matrix(c(6, 2, 4), 3, 1) > B [,1] [1,] 6 [2,] 2 [3,] 4 > A %*% B [,1] [1,] 54
文字式で表すとより分かりやすい。行列同士の掛け算はこのように行方向の並びの成分と列方向の並びの成分とを順番に掛けて1行1列ごとの合計を取る。
z | |||
y | |||
x | |||
a | b | c | az + by + cx |
行列を行ベクトルが複数重なったもの、または、列ベクトルが複数並んだものと考えてみる。例えば次の行列Cは、2次の行ベクトル[a b]と2次の行ベクトル[c d]と2次の行ベクトル[e f]が3つ重なったもの、と。
\[ C_1 = \begin{pmatrix} a & b \end{pmatrix} \\
C_2 = \begin{pmatrix} c & d \end{pmatrix} \\
C_3 = \begin{pmatrix} e & f \end{pmatrix} \]
行列Dは、2次の列ベクトル[j m]と2次の列ベクトル[k n]と2次の列ベクトル[l o]が3つ並んだもの、と。
\[ D_1 = \begin{pmatrix} j \\ m \end{pmatrix}
D_2 = \begin{pmatrix} k \\ n \end{pmatrix}
D_3 = \begin{pmatrix} l \\ o \end{pmatrix} \]
このように考えると、行列Cと行列Dとの掛け算が2次の行ベクトルと2次の列ベクトルの掛け算の応用であることが分かる。
&= \begin{pmatrix} C_1 \times D_1 & C_1 \times D_2 & C_1 \times D_3 \\ C_2 \times D_1 & C_2 \times D_2 & C_2 \times D_3 \\ C_3 \times D_1 & C_3 \times D_2 & C_3 \times D_3 \end{pmatrix} \\
&= \begin{pmatrix} aj + bm & ak + bn & al + bo \\ cj + dm & ck + dn & cl + do \\ ej + fm & ek + fn & el + fo \end{pmatrix}
\end{align} \]
行ベクトルと列ベクトルの掛け算に分解してみる。おのおのの成分は行ベクトルと列ベクトルの掛け算で構成されていることがよく分かる。
&= \begin{pmatrix} a & b \end{pmatrix} \times \begin{pmatrix} j \\ m \end{pmatrix} \\
&= aj + bm \end{align} \]
\[ \begin{align} C_1 \times D_2
&= \begin{pmatrix} a & b \end{pmatrix} \times \begin{pmatrix} k \\ n \end{pmatrix} \\
&= ak + bn \end{align} \]
\[ \begin{align} C_1 \times D_3
&= \begin{pmatrix} a & b \end{pmatrix} \times \begin{pmatrix} l \\ o \end{pmatrix} \\
&= al + bo \end{align} \]
\[ \begin{align} C_2 \times D_1
&= \begin{pmatrix} c & d \end{pmatrix} \times \begin{pmatrix} j \\ m \end{pmatrix} \\
&= cj + dm \end{align} \]
\[ \begin{align} C_2 \times D_2
&= \begin{pmatrix} c & d \end{pmatrix} \times \begin{pmatrix} k \\ n \end{pmatrix} \\
&= ck + dn \end{align} \]
\[ \begin{align} C_2 \times D_3
&= \begin{pmatrix} c & d \end{pmatrix} \times \begin{pmatrix} l \\ o \end{pmatrix} \\
&= cl + do \end{align} \]
\[ \begin{align} C_3 \times D_1
&= \begin{pmatrix} e & f \end{pmatrix} \times \begin{pmatrix} j \\ m \end{pmatrix} \\
&= ej + fm \end{align} \]
\[ \begin{align} C_3 \times D_2
&= \begin{pmatrix} e & f \end{pmatrix} \times \begin{pmatrix} k \\ n \end{pmatrix} \\
&= ek + fn \end{align} \]
\[ \begin{align} C_3 \times D_3
&= \begin{pmatrix} e & f \end{pmatrix} \times \begin{pmatrix} l \\ o \end{pmatrix} \\
&= el + fo \end{align} \]
j | k | l | ||
m | n | o | ||
a | b | aj + bm | ak + bn | al + bo |
c | d | cj + dm | ck + dn | cl + do |
e | f | ej + fm | ek + fn | el + fo |
それでは実際の計算をやってみる。
\[ B = \begin{pmatrix} 8 & 7 \\ 6 & 5 \end{pmatrix} \]
\[ \begin{align} AB
&= \begin{pmatrix} 1 \times 8 + 2 \times 6 & 1 \times 7 + 2 \times 5 \\
3 \times 8 + 4 \times 6 & 3 \times 7 + 4 \times 5 \end{pmatrix} \\
&= \begin{pmatrix} 20 & 17 \\ 48 & 41 \end{pmatrix}
\end{align} \]
計算結果が正しいかどうかR言語で確認。
> A = matrix(c(1, 3, 2, 4), 2, 2) > A [,1] [,2] [1,] 1 2 [2,] 3 4 > B = matrix(c(8, 6, 7, 5), 2, 2) > B [,1] [,2] [1,] 8 7 [2,] 6 5 > A %*% B [,1] [,2] [1,] 20 17 [2,] 48 41
\[ D = \begin{pmatrix} -1 & 7 \\ 6 & -5 \\ 2 & 0 \end{pmatrix} \]
\[ \begin{align} CD
&= \begin{pmatrix} -5 \times -1 + 3 \times 6 + 2 \times 2
& -5 \times 7 + 3 \times -5 + 2 \times 0 \end{pmatrix} \\
&= \begin{pmatrix} 27 & -50 \end{pmatrix}
\end{align} \]
計算結果が正しいかどうかR言語で確認。
> C = matrix(c(-5, 3, 2), 1, 3) > C [,1] [,2] [,3] [1,] -5 3 2 > D = matrix(c(-1, 6, 2, 7, -5, 0), 3, 2) > D [,1] [,2] [1,] -1 7 [2,] 6 -5 [3,] 2 0 > C %*% D [,1] [,2] [1,] 27 -50
\[ A = \begin{pmatrix} 1 & 2 \\ 3 & 4 \end{pmatrix} \]
\[ \begin{align} DA
&= \begin{pmatrix} -1 \times 1 + 7 \times 3 & -1 \times 2 + 7 \times 4 \\
6 \times 1 + -5 \times 3 & 6 \times 2 + -5 \times 4 \\
2 \times 1 + 0 \times 3 & 2 \times 2 + 0 \times 4 \end{pmatrix} \\
&= \begin{pmatrix} 20 & 26 \\ -9 & -8 \\ 2 & 4 \end{pmatrix}
\end{align} \]
R言語で計算結果が正しいかどうか確認。
> D %*% A [,1] [,2] [1,] 20 26 [2,] -9 -8 [3,] 2 4
行列同士の掛け算の重要な特徴として交換法則が成り立たないことがある。そのことをR言語で確認してみる。
> A = matrix(c(1, 3, 2, 4), 2, 2) > A [,1] [,2] [1,] 1 2 [2,] 3 4 > B = matrix(c(8, 6, 7, 5), 2, 2) > B [,1] [,2] [1,] 8 7 [2,] 6 5 > A %*% B [,1] [,2] [1,] 20 17 [2,] 48 41 > B %*% A [,1] [,2] [1,] 29 44 [2,] 21 32
\[ B = \begin{pmatrix} 8 & 7 \\ 6 & 5 \end{pmatrix} \]
\[ AB \neq BA \]
まとめ
- 行列同士の足し算
- 互いの行列の同じ位置にある各要素同士を単純に足し算すればいい。実数同士の足し算と同じ法則、交換法則と結合法則が成り立つ。
- 行列同士の引き算
- 互いの行列の同じ位置にある各要素同士を単純に引き算すればいい。
- 行列の実数倍(スカラー倍)
- 行列と実数との掛け算。行列の全ての要素に実数を単純に掛ければいい。実数同士の掛け算と同じ法則、交換法則と結合法則と分配法則が成り立つ。
- 行列同士の掛け算
- 実数同士の掛け算と異なる。行列A掛ける行列Bが成り立つのは行列Aの列数と行列Bの行数とが等しいときだけに限られる。交換法則が成り立たない。
-
行列同士の乗法 a c b d 2 3 2a + 3b 2c + 3d 4 5 4a + 5b 4c + 5d - R言語では行列同士の掛け算の演算子は%*%であって*ではない。ちなみにPythonやRubyでは*でいい。
- R言語の公式ホームページ
- https://www.r-project.org/
コメント
コメントを投稿