線形代数

一次独立

いずれも 0\bold{0} でないベクトル a1,a2,…,an\bold{a_1,a_2,\dots,a_n} について c1a1+c2a2+⋯+cnan=0c_1\bold{a_1}+c_2\bold{a_2}+\cdots+c_n\bold{a_n}=\bold{0} が成り立つのが c1=c2=⋯=cn=0c_1=c_2=\cdots=c_n=0 だけの場合、a1,a2,…,an\bold{a_1,a_2,\dots,a_n} を一次独立 (linearly independent) という。 また、一次独立でないものを一次従属 (linearly dependent) という。

定理

あるベクトル b\bold{b} が一次独立なベクトル a1,a2,…,an\bold{a_1,a_2,\dots,a_n} の線型結合で表されるとき、その表し方は1通りである。

階数

任意の行列 AA は行基本変形を繰り返すことによって階段行列にすることができる。

行列 AA の階段行列の中の少なくとも1つは 00 でない成分をもつ行の個数を行列 AA の階数 (rank) といい、以下 rank(A)\text{rank}(A) で表す。

nn 元連立方程式の係数行列を AA、定数項ベクトルを b\bold{b}、係数行列と定数項ベクトルを連結した拡大係数行列を AbAb とおくと、以下の関係が成り立つ。

rank(A)=rank(Ab)∧{rank(A)=n  ⟺  解は1つrank(A)<n  ⟺  解は不定であり、n−rank(A) 個の任意定数で表せるrank(A)≠rank(Ab)  ⟺  解なし (不能)\begin{aligned} \text{rank}(A)&=\text{rank}(A\bold{b})\land\begin{cases} \text{rank}(A)=n&\iff\text{解は1つ}\\ \text{rank}(A)\lt n&\iff\text{解は不定であり、}n-\text{rank}(A)\,\text{個の任意定数で表せる} \end{cases}\\ \text{rank}(A)&\ne\text{rank}(A\bold{b})\iff\text{解なし (不能)} \end{aligned}
補足

特に b=0\bold{b}=\bold{0} のとき、必ず rank(A)=rank(Ab)\text{rank}(A)=\text{rank}(A\bold{b}) であるため、連立方程式は必ず解をもつ。

行列式

行列式 (determinant) とは正方行列に対して定義される量であり、以下行列 AA の行列式を ∣A∣|A| で表す。

性質

  • 転置不変性

    ∣A∣=∣At∣|A|=|A^t|
  • 交代性

    行 or 列を入れ替えると、行列式の値は −1-1 倍になる。

  • 多重線形性

    ある行 or 列の定数倍を他の行 or 列に加えても行列式の値は変化しない。

求め方: サラスの方法

サラスの方法は、2次あるいは3次正方行列において、左下から右下へ向かう方向に ++, 右上から左下に向う方向に −- の符号を付けて積をとり、それらの和をとることで行列式を求める。

例

2次正方行列 A=(a11a12a21a22)A=\left(\begin{matrix} a_{11}&a_{12}\\ a_{21}&a_{22} \end{matrix}\right) の行列式は

∣A∣=a11a22−a12a21|A|=a_{11}a_{22}-a_{12}a_{21}

で計算できる。

3次正方行列 A=(a11a12a13a21a22a23a31a32a33)A=\left(\begin{matrix} a_{11}&a_{12}&a_{13}\\ a_{21}&a_{22}&a_{23}\\ a_{31}&a_{32}&a_{33} \end{matrix}\right) の行列式は

∣A∣=a11a22a33+a12a23a31+a13a21a32−a31a22a13−a32a23a11−a33a21a12|A|=a_{11}a_{22}a_{33}+a_{12}a_{23}a_{31}+a_{13}a_{21}a_{32}-a_{31}a_{22}a_{13}-a_{32}a_{23}a_{11}-a_{33}a_{21}a_{12}

で計算できる。

ただし、4次以上の正方行列に対してサラスの方法は使えない。

求め方: 余因子展開

行列 AA から ii 行 jj 列を取り除いた小行列式を Mi,jM_{i,j} と表すとき

a~i,j=(−1)i+j∣Mi,j∣\tilde{a}_{i,j}=(-1)^{i+j}|M_{i,j}|

を AA の**(i,j)(i,j) 余因子 (cofactor)** といい、a~i,j\tilde{a}_{i,j} で表す。

行列 AA を nn次正方行列とすると、AA の行列式は次のように与えられる。

第 ii 列に沿った余因子展開

∣A∣=ai,1(−1)i+1∣Mi,1∣+ai,2(−1)i+2∣Mi,2∣+⋯+ai,n(−1)i+n∣Mi,n∣=∑j′=1nai,j′(−1)i+j′∣Mi,j′∣=∑j′=1nai,j′a~i,j′\begin{aligned} |A|&=a_{i,1}(-1)^{i+1}|M_{i,1}|+a_{i,2}(-1)^{i+2}|M_{i,2}|+\cdots+a_{i,n}(-1)^{i+n}|M_{i,n}|\\ &=\sum_{j'=1}^n{a_{i,j'}(-1)^{i+j'}|M_{i,j'}|}\\ &=\sum_{j'=1}^n{a_{i,j'}\tilde{a}_{i,j'}} \end{aligned}

第 jj 行に沿った余因子展開

∣A∣=a1,j(−1)1+j∣M1,j∣+a2,j(−1)2,j∣M2,j∣+⋯+an,j(−1)n+j∣Mn,j∣=∑i′=1nai′,j(−1)i′+j∣Mi′+j∣=∑i′=1nai′,ja~i′,j\begin{aligned} |A|&=a_{1,j}(-1)^{1+j}|M_{1,j}|+a_{2,j}(-1)^{2,j}|M_{2,j}|+\cdots+a_{n,j}(-1)^{n+j}|M_{n,j}|\\ &=\sum_{i'=1}^n{a_{i',j}(-1)^{i'+j}|M_{i'+j}|}\\ &=\sum_{i'=1}^n{a_{i',j}\tilde{a}_{i',j}} \end{aligned}

例題

  • ∣102−1−11212∣\begin{vmatrix}1&0&2\\-1&-1&1\\2&1&2\end{vmatrix} を第1行に沿って余因子展開せよ
  • ∣147258369∣\begin{vmatrix}1&4&7\\2&5&8\\3&6&9\end{vmatrix} を第1列に沿って余因子展開せよ
  • ∣−210210−12−4231201−1∣\begin{vmatrix}-2&1&0&2\\1&0&-1&2\\-4&2&3&1\\2&0&1&-1\end{vmatrix} を第2列に沿って余因子展開せよ

逆行列

正方行列 AA に対して AX=XA=EAX=XA=E となる正方行列 XX が存在するとき、行列 AA は「正則である (regular)」という。 このとき、XX を AA の逆行列 (inverse matrix) といい、以下 A−1A^{-1} で表す。

nn次正方行列 AA が正則であるとき、その逆行列 A−1A^{-1} は

A−1=1∣A∣A~A^{-1}=\frac1{|A|}\tilde{A}

と表せる。

ここで、A~\tilde{A} を AA の**余因子行列 (cofactor matrix)**といい、以下のように定義される。

A~=(a~11a~12⋯a~1na~21a~22⋯a~2n⋮⋮⋱⋮a~n1a~n2⋯a~nn)t=(a~11a~21⋯a~n1a~12a~22⋯a~n2⋮⋮⋱⋮a~1na~2n⋯a~nn)\tilde{A}=\left(\begin{array}{ccccc} \tilde{a}_{11}&\tilde{a}_{12}&\cdots&\tilde{a}_{1n}\\ \tilde{a}_{21}&\tilde{a}_{22}&\cdots&\tilde{a}_{2n}\\ \vdots&\vdots&\ddots&\vdots\\ \tilde{a}_{n1}&\tilde{a}_{n2}&\cdots&\tilde{a}_{nn}\\ \end{array}\right)^t=\left(\begin{array}{ccccc} \tilde{a}_{11}&\tilde{a}_{21}&\cdots&\tilde{a}_{n1}\\ \tilde{a}_{12}&\tilde{a}_{22}&\cdots&\tilde{a}_{n2}\\ \vdots&\vdots&\ddots&\vdots\\ \tilde{a}_{1n}&\tilde{a}_{2n}&\cdots&\tilde{a}_{nn}\\ \end{array}\right)

性質

n次正方行列が正則  ⟺  rank(A)=n  ⟺  ∣A∣≠0n\text{次正方行列が正則}\iff\text{rank}(A)=n\iff|A|\ne0

例題

  • A=(abcd)A=\left(\begin{matrix}a&b\\c&d\end{matrix}\right) の逆行列を求めよ

固有ベクトルと固有値

ある nn次正方行列 AA に対し Ax=λxA\bold{x}=\lambda\bold{x} を満たす nn 次元列ベクトル x  (x≠0)\bold{x}\;(\bold{x}\ne\bold{0}) が存在するとする。 このとき、λ\lambda を AA の固有値 (eigenvalue) といい、x\bold{x} を λ\lambda に対する固有ベクトル (eigenvector) という。

また、AA の相異なる固有ベクトル λ1,λ2,…,λk\lambda_1,\lambda_2,\dots,\lambda_k は一次独立である。

補足

固有値と固有ベクトルの意味の考察

Ax=λxA\bold{x}=\lambda\bold{x}

ベクトル x\bold{x} を行列 AA によって変換したときに、変換後のベクトル AxA\bold{x} が変換前のベクトルのスカラー倍になる (変換前後でベクトルの向きが変わらない) とき、変換前のベクトルが固有ベクトル x\bold{x} でそのときの倍率が固有値 λ\lambda となる。

求め方

Ax=λx(A−λE)x=0\begin{aligned} A\bold{x}&=\lambda\bold{x}\\ (A-\lambda E)\bold{x}&=\bold{0} \end{aligned}

ここで、上記の式が x=0\bold{x}=\bold{0} 以外の解をもつことから

∣A−λE∣=0\begin{aligned} |A-\lambda E|=0 \end{aligned}

が導け、これを固有方程式という。 固有方程式は λ\lambda についての nn 元連立方程式である。

例題

  • A=(2341)A=\left(\begin{matrix}2&3\\4&1\end{matrix}\right) の固有値・固有ベクトルを求めよ
  • A=(211121112)A=\left(\begin{matrix}2&1&1\\1&2&1\\1&1&2\end{matrix}\right) の固有値・固有ベクトルを求めよ

対角化

nn次正方行列 AA に対し、適切な正則行列 PP が存在して

P−1AP=(λ10⋯00λ2⋯0⋮⋮⋱⋮00⋯λn)P^{-1}AP=\left(\begin{array}{cccc} \lambda_1&0&\cdots&0\\ 0&\lambda_2&\cdots&0\\ \vdots&\vdots&\ddots&\vdots\\ 0&0&\cdots&\lambda_n \end{array}\right)

のような対角行列にできる時、行列 AA は対角化可能であるといい、このときの行列 PP を変換行列という。

ここで λ0,λ1 ⋯ ,λn\lambda_0,\lambda_1\,\cdots,\lambda_n は行列 AA の固有値である。

定理

nn次正方行列 AA の nn 個の一次独立な固有ベクトルを x1,x2,⋯ ,xn\bold{x_1},\bold{x_2},\cdots,\bold{x_n} とする。 それらを並べた行列 (x1,x2,⋯ ,xn)(\bold{x_1},\bold{x_2},\cdots,\bold{x_n}) を PP とすると、行列 AA は次のように対角化できる。

P−1AP=(λ10⋯00λ2⋯0⋮⋮⋱⋮00⋯λn)P^{-1}AP=\left(\begin{array}{cccc} \lambda_1&0&\cdots&0\\ 0&\lambda_2&\cdots&0\\ \vdots&\vdots&\ddots&\vdots\\ 0&0&\cdots&\lambda_n \end{array}\right)

対角化できない例

A=(−3−11−1)A=\left(\begin{matrix} -3&-1\\ 1&-1 \end{matrix}\right)

行列 AA の固有方程式を解くと

∣−3−λ−11−1−λ∣=(−3−λ)(−1−λ)+1=λ2+4λ+4=(λ+2)2=0∴λ=−2(−1−111)(xy)=(00)  ⟺  x+y=0\begin{aligned} \begin{vmatrix} -3-\lambda&-1\\ 1&-1-\lambda \end{vmatrix}&=(-3-\lambda)(-1-\lambda)+1\\ &=\lambda^2+4\lambda+4\\ &=(\lambda+2)^2=0\\ \therefore\lambda&=-2\\ \left(\begin{matrix} -1&-1\\ 1&1 \end{matrix}\right)\left(\begin{matrix} x\\ y \end{matrix}\right)&=\left(\begin{matrix} 0\\ 0 \end{matrix}\right)\iff{x+y=0} \end{aligned}

ここで x=sx=s とおくと y=−sy=-s

∴(xy)=s(1−1)\therefore{\left(\begin{matrix} x\\ y \end{matrix}\right)=s\left(\begin{matrix} 1\\ -1 \end{matrix}\right)}

これ以外に独立な固有ベクトルがとれないので対角化不可能。

応用例: 対角化を使った A^n の計算

(P−1AP)k=(P−1AP)(P−1AP)⋯(P−1AP)=P−1AkPAk=P(P−1AP)kP−1=P(λ1k0⋯00λ2k⋯0⋮⋮⋱⋮00⋯λnk)P−1\begin{aligned} (P^{-1}AP)^k&=(P^{-1}AP)(P^{-1}AP)\cdots(P^{-1}AP)\\ &=P^{-1}A^kP\\ A^k&=P(P^{-1}AP)^kP^{-1}\\ &=P\left(\begin{matrix} \lambda_1^k&0&\cdots&0\\ 0&\lambda_2^k&\cdots&0\\ \vdots&\vdots&\ddots&\vdots\\ 0&0&\cdots&\lambda_n^k \end{matrix}\right)P^{-1} \end{aligned}

と求められる。

例題

  • 行列 A=(−2152)A=\left(\begin{matrix}-2&1\\5&2\end{matrix}\right) を変換行列 PP を用いて対角化せよ
  • 行列 A=(−224−232−214)A=\left(\begin{matrix}-2&2&4\\-2&3&2\\-2&1&4\end{matrix}\right) を変換行列 PP を用いて対角化せよ

参考文献

前の記事

微分方程式

次の記事

高校数学