線形代数I/行列の階数 のバックアップの現在との差分(No.2)

更新


  • 追加された行はこの色です。
  • 削除された行はこの色です。
[[線形代数I]]

#contents
&katex();

* 階段行列 [#h48fd710]
培風館「教養の線形代数(五訂版)」に沿って行っている授業の授業ノート(の一部)です。

* 2.3 変形定理(掃き出しの定理)・階段行列 [#h48fd710]

[[線形代数I/連立一次方程式]] でやったように、
行列に対してガウスの消去法を適用することで、
例えば以下のような形に変形できる。
任意の行列 &math(A); に対してガウスの消去法を適用することで、
例えば以下のような形に変形できる。(*は任意の数値)

&math(\begin{bmatrix}0&1&*&0&*&*&0&*\\0&0&0&1&*&*&0&*\\0&0&0&0&0&0&1&*\\0&0&0&0&0&0&0&0\\0&0&0&0&0&0&0&0\\\end{bmatrix});
&math(A\sim\begin{bmatrix}0&1&*&0&*&*&0&*\\0&0&0&1&*&*&0&*\\0&0&0&0&0&0&1&*\\0&0&0&0&0&0&0&0\\0&0&0&0&0&0&0&0\\\end{bmatrix});

このように「これ以上掃き出せない」形の行列を「階段行列」と呼ぶ。
行に対する基本変形では「これ以上掃き出せない」このような行列を「階段行列」と呼ぶ。

一般形は、

&math(\left[\begin{array}{c@{\,}c@{\,}ccc@{\,}c@{\,}ccc@{\,}c@{\,}cccccc@{\,}c@{\,}c}0&\cdots&0&1&*&\cdots&*&0&*&\cdots&*&0&*&{\ \ \ \ \ \ \ \ \ \ \ }&0&*&\cdots&*\\\vdots&&\vdots&0&0&\cdots&0&1&*&\cdots&*&0&*&&\vdots&\vdots&&\vdots\\\vdots&&\vdots&\vdots&\vdots&&\vdots&0&0&\cdots&0&1&*&\cdots&\vdots&\vdots&&\vdots\\\vdots&&\vdots&\vdots&\vdots&&\vdots&\vdots&\vdots&&\vdots&0&0&&\vdots&\vdots&&\vdots\\\vdots&&\vdots&\vdots&\vdots&&\vdots&\vdots&\vdots&&\vdots&\vdots&\vdots&&\vdots&\vdots&&\vdots\\\vdots&&\vdots&\vdots&\vdots&&\vdots&\vdots&\vdots&&\vdots&\vdots&\vdots&&0&\vdots&&\vdots\\\vdots&&\vdots&\vdots&\vdots&&\vdots&\vdots&\vdots&&\vdots&\vdots&\vdots&&1&*&\cdots&*\\\vdots&&\vdots&\vdots&\vdots&&\vdots&\vdots&\vdots&&\vdots&\vdots&\vdots&\cdots&0&0&\cdots&0\\\vdots&&\vdots&\vdots&\vdots&&\vdots&\vdots&\vdots&&\vdots&\vdots&\vdots&&\vdots&\vdots&&\vdots\\0&\cdots&0&0&0&\cdots&0&0&0&\cdots&0&0&0&&0&0&\cdots&0\\\end{array}\right]);

である。

** 階段行列の形 [#hb2834aa]

- 掃き出せた列
-- 1が一つあって、後は0
-- 左の列よりも1つ段を上がる
- 掃き出せなかった列
-- 左の列より段は上がらない

** 行に対する基本変形 [#h8a08c91]

ガウスの消去法で用いる

+ ある行を定数倍(&math(\neq 0);)する
+ ある行を定数倍して別の行に加える
+ ある行を別の式を入れ替える

の3つの操作(同値変形)を、「行に対する基本変形」と呼ぶ。

** 行列の階段化 [#xe6a7a74]

任意の行列は「行に対する基本変形」により「階段行列」にできる~
→ これを「行列の階段化」と呼ぼう

* 2.2 基本変形を表す行列(基本行列) [#n979851f]

行に対する3つの基本変形を、行列のかけ算で表すことができる。

** $i$ 行目を $c$ 倍する ( $c \ne 0$ ) [#g727c5f5]

&math(P_i(c) = \left[\begin{array}{ccccccc}1&&&&&&\\&\ddots&&&&0&\\&&1&&&&\\&&&c&&&\\&&&&1&&\\&0&&&&\ddots&\\&&&&&&1\end{array}\right]);

は、単位行列の ''i 行目を &math(c); 倍した''行列~
単位行列の &math((i,i)); 要素が &math(c); になっている。

対象の行列に左から &math(P_i(c)); を掛けることで基本変形が行われる。

例:&math(P_3(2)A=\begin{bmatrix}1&0&0\\0&1&0\\0&0&2\end{bmatrix}\begin{bmatrix}a&b\\c&d\\e&f\end{bmatrix}=\begin{bmatrix}a&b\\c&d\\2e&2f\end{bmatrix});

ここで &math(P_i(c)^{-1}=P_i(1/c)); であるから、 &math(P_i(c)); は正則である。(&math(c\ne 0);)

例:&math(P_3(2)P_3(1/2)=\begin{bmatrix}1&0&0\\0&1&0\\0&0&2\end{bmatrix}\begin{bmatrix}1&0&0\\0&1&0\\0&0&1/2\end{bmatrix}=\begin{bmatrix}1&0&0\\0&1&0\\0&0&1\end{bmatrix}); 逆から掛けても同様。

** $i$ 行目に $j$ 行目の $c$ 倍を加える [#y2269bc7]

&math(P_{ij}(c) = \left[\begin{array}{ccccc}1&&&&0\\&\ddots&&&\\&&\ddots&c&\\&&&\ddots&\\0&&&&1\end{array}\right]);

は、単位行列の ''&math(i); 行目に &math(j); 行目の &math(c); 倍を加えた''行列~
単位行列の&math((i,j));要素が &math(c); になっている。

対象の行列に左から &math(P_{ij}(c)); を掛けることで基本変形が行われる。

例:&math(P_{13}(2)A=\begin{bmatrix}1&0&2\\0&1&0\\0&0&1\end{bmatrix}\begin{bmatrix}a&b\\c&d\\e&f\end{bmatrix}=\begin{bmatrix}a+2e&b+2f\\c&d\\e&f\end{bmatrix});

ここで &math(P_{ij}(c)^{-1}=P_{ij}(-c)); であるから、 &math(P_{ij}(c)); は正則である。

例:&math(P_{13}(2)P_{13}(-2)=\begin{bmatrix}1&0&2\\0&1&0\\0&0&1\end{bmatrix}\begin{bmatrix}1&0&-2\\0&1&0\\0&0&1\end{bmatrix}=\begin{bmatrix}1&0&0\\0&1&0\\0&0&1\end{bmatrix}); 逆から掛けても同様。

** $i$ 行目と $j$ 行目を入れ替える [#la920c65]

&math(P_{ij} = \left[\begin{array}{ccccccc}1&&&&&&\\&\ddots&&&&0&\\&&0&&1&&\\&&&\ddots&&&\\&&1&&0&&\\&0&&&&\ddots&\\&&&&&&1\end{array}\right]);

は、単位行列の ''&math(i); 行目と &math(j); 行目を入れ替えた''行列~
単位行列の&math((i,i), (j,j));要素が 0 になり、&math((i,j), (j,i)); 要素が 1 になっている。

対象の行列に左から &math(P_{ij}); を掛けることで基本変形が行われる。

例:&math(P_{23}A=\begin{bmatrix}1&0&0\\0&0&1\\0&1&0\end{bmatrix}\begin{bmatrix}a&b\\c&d\\e&f\end{bmatrix}=\begin{bmatrix}a&b\\e&f\\c&d\end{bmatrix});

ここで &math(P_{ij}^{-1}=P_{ij}); であるから、 &math(P_{ij}); は正則である。

例:&math(P_{23}P_{23}=\begin{bmatrix}1&0&0\\0&0&1\\0&1&0\end{bmatrix}\begin{bmatrix}1&0&0\\0&0&1\\0&1&0\end{bmatrix}=\begin{bmatrix}1&0&0\\0&1&0\\0&0&1\end{bmatrix});

** 基本行列・左基本変形 [#n45a17bb]

このように、行に対する基本変形は、「ある正則行列」を左から掛けることと同一視できる。

そこで、「行に対する基本変形」を「左基本変形」とも呼ぶ。~
(後に出てくる「列に対する基本変形」を「右基本変形」と呼ぶ)

「ある正則行列」は上で見た &math(P_i(c));、&math(P_{ij}(c));、&math(P_{ij}); のいずれかであり、
これらの行列は「基本行列」と呼ばれる。

・「基本行列」は正則である

* 2.3 行列で表した掃出しの定理 (定理 2.2 変形定理) [#d59a814a]

・任意の行列 &math(A); に、左から基本行列を複数掛けることで、階段行列に変形できる。

&math(\underbrace{P_kP_{k-1}\cdots P_2P_1}_PA=(階段行列));

&math(P_kP_{k-1}\cdots P_2P_1=P); と置けば、&math(P); は正則なので(∵ &math(P^{-1}=P_1^{-1}P_2^{-1}\cdots P_k^{-1});)

「任意の行列 &math(A); に対して、&math(PA); が階段行列になるような正則行列 &math(P); が存在する」

と言うこともできる。

この定理は非常に強力であり、以下に見るようにこれを用いて数々の重要な定理を導ける。

* 行列の階段化 [#e84747fb]

行列 &math(A); に左から正則行列 &math(P); を掛けて &math(PA); を階段行列にすることを &math(A); を「階段化する」と言う。

* A が正則なら、階段化すると単位行列になる(定理2.6の一部) [#t0a8974e]

単位行列は階段行列の一種であることに注意せよ。

さて、正則行列とは逆行列を持つ正方行列であった。

&math(A); を階段化して &math(PA=X); となったとしよう。

&math(A); も &math(P); も正則なので &math(X); も正則である。(∵&math(X^{-1}=A^{-1}P^{-1});)

正則な階段行列は単位行列しかない。

なぜなら、正則な行列は必ず正方行列であり、また正方行列では
- 全ての列が掃き出されれば単位行列となる
- 掃き出されない列が1つでもあると、最後の行がゼロになる~
→ ゼロ行ベクトルを含む行列は正則行列になり得ない

したがって、&math(A); が正則なら、階段化すると &math(PA=I); となる。

** (補題) ゼロ行ベクトルを含む行列は正則行列になり得ない [#xfa76cc6]

&math(A); の &math(k); 行目がゼロだとすると、
任意の &math(B); について &math(AB); の &math(k); 行目もゼロになってしまう。

&math(
AB=
\begin{matrix}
\\
\\
\scriptstyle k\\
\\
\end{matrix}
\underbrace{
\begin{bmatrix}
\cdots&\cdots&\cdots\\
\cdots&\cdots&\cdots\\
0&\cdots&0\\
\cdots&\cdots&\cdots
\end{bmatrix}}_{\displaystyle A}
B=
\begin{matrix}
\\
\\
\scriptstyle k\\
\\
\end{matrix}\underbrace{
\begin{bmatrix}\cdots&\cdots&\cdots\\\cdots&\cdots&\cdots\\0&\cdots&0\\\cdots&\cdots&\cdots\end{bmatrix}}_{\displaystyle AB}
);

したがって、任意の &math(B); について &math(AB\ne I); で、
&math(A); は逆行列を持たない。

** 逆行列の求め方 (再) [#w803fe41]

この定理から、&math(A); が正則であるとき、&math(A); を階段化して
&math(PA=I); と書いた時の &math(P); が逆行列(の候補)となる。

これを利用して、&math(A); と &math(I); を並べた行列を階段化すれば、

&math(\begin{bmatrix}A&I\end{bmatrix}\sim \begin{bmatrix}PA&P\end{bmatrix}\sim \begin{bmatrix}I&A^{-1}\end{bmatrix});

となり、左側が単位行列になった時、右側に逆行列(の候補)&math(P); が得られる。

次項にて、&math(P); が本当に逆行列であることを言うために &math(AP=I); を確かめよう

* A, B が正方行列の時、AB = I なら BA = I である (2章問6) [#za525b69]

∵ &math(A); を階段化して &math(PA=X); となったとする。

&math(AB=I);

に左から &math(P); を掛ければ、

&math(PAB&=(PA)B=XB\\&=P(AB)=PI=P);~
より

&math(XB=P);

&math(X); は階段行列でかつ正方行列なので、もし単位行列でなければ最終行がゼロになる。

&math(X); がゼロ行ベクトルを含む場合、&math(XB); の同じ行がゼロとなり、
&math(P); が正則であるとの条件と矛盾する。

したがって、&math(X); は単位行列である &math(X=I);

すると、&math(XB=IB=B=P); となり、

&math(BA=I); を得る。

つまり、&math(AB=I); であれば、&math(B=A^{-1}); である。

(&math(A^{-1}); の定義は、&math(AA^{-1}=I); と &math(A^{-1}A=I); 
が両方成り立たつことであった。)

- この定理により、&math(PA=I); となる &math(P); を求めれば、それが逆行列であることを言える

* A が正則ならば、A を基本行列の積で表すことができる (定理2.7) [#m1920bc7]

&math(A); を階段化すれば &math(PA=I);

すなわち &math(P=P_k\cdots P_1=A^{-1}); であるから、、

&math(A=P_1^{-1}\cdots P_k^{-1}); と表せる。

基本行列の逆行列は基本行列であり、これは &math(A); 
を基本行列の積として表したことになる。

* 行列の階数 [#o8ed4035]

行列 &math(A); を正則行列 &math(P); を用いて

&math(PA = ); (階段行列)

の形に階段化した際の、階段行列の段数(ゼロでない行数)を「&math(A); の階数」と呼び、

&math({\rm rank}A);

と表す。

-「階数」は一意に決まるか???

** ある行列 A を階段化した形は一意に決定される (定理2.11) [#r407ed63]

> 一般に、ある行列を階段化する手順は複数あるものの、~
最終的な階段行列の形は手順によらず同じ形になる。

この定理は重要であるが、証明は時間の都合で割愛する。~
教科書では 定理2.11 として扱われており、その証明は面倒ではあるものの初等的である。

この定理から導かれる結論として、

- 階数は一意に決まる

** 階数の性質 [#oaf41683]

- &math({\rm rank}A); は「掃き出しの行えた数」に等しいから、
-- &math({\rm rank}A\le ); ( A の列数 )
-- &math({\rm rank}A\le ); ( A の行数 )

- &math(n);次の正方行列 &math(A); に対して、以下の3つは同値な条件
-- &math(A); は正則
-- &math(PA=I_n);
-- &math({\rm rank}A=n);

** 正則行列を掛けても rank は変化しない [#fd0869f7]

&math(B); が正則なら &math({\rm rank}(BA)={\rm rank}A);

∵&math(PA=X); を階段行列とすると、~
&math((PB^{-1})BA=PA=X); となり、&math(X); は &math(PB^{-1}); により &math(BA); を階段化した階段行列でもある。

** 転置により rank は変化しない [#beb5d621]

この定理を証明するために、教科書では行列の標準形の概念を導入しているが、ここでは割愛する。

- &math({\rm rank}\,{}^t\!\!A={\rm rank}\,A); (定理2.13)

この定理から、&math(B);、&math(C); が正則なとき、
- &math({\rm rank}\,BAC={\rm rank}\,A);

つまり、行列 &math(A); に左右どちらから正則行列を掛けても rank は変化しないことが導かれる。

∵ &math(\rank BAC=\rank AC=\rank \transpose (AC)=\rank (\transpose C \transpose A)=\rank \transpose A=\rank A);~

** 行列の積により rank が減ることはあっても増えることはない [#f3032696]

定理2.13
- &math({\rm rank}\,BA\le{\rm rank}\,A);
- &math({\rm rank}\,BA\le{\rm rank}\,B);

これも重要な定理であるが、証明は割愛して結果のみを示す。

この定理から2つの重要な定理を導いておく。

** 逆行列を持つのは正方行列に限る [#zdedb319]

&math(A); の逆行列 &math(A^{-1}); とは、
&math(AA^{-1}=I);、&math(A^{-1}A=I); を満たす行列として定義された。

&math(A); を &math(m\times n);行列、~
&math(B); を &math(n\times m);行列とすれば、~
&math(m\ne n); であっても、~
&math(m);次正方行列 &math(AB); と~
&math(n);次正方行列 &math(BA); はどちらも定義される。~
しかし、この両方が単位行列になることはありえない。

∵ &math(m<n); のとき、~
&math({\rm rank}\,A\le m);~
&math({\rm rank}\,B\le m);~
より、~
&math({\rm rank}\,BA\le m); ~
である。

一方、&math(n);次単位行列 &math(I_n); の rank は &math(n); だから、
&math(BA=I_n); となることはありえない。

注)&math(m<n); のとき、&math(AB=I_m); となることはあり得る。

例:&math(\begin{bmatrix}1&0&0\\0&1&0\end{bmatrix}\begin{bmatrix}1&0\\0&1\\0&0\end{bmatrix}=\begin{bmatrix}1&0\\0&1\end{bmatrix});

** 正方行列の積が正則であれば、元の行列も正則である [#s001bd94]

&math(n); 次正方行列 &math(A,B); があり、&math(A); が正則でないとする。

すなわち &math({\rm rank}\,A<n);

すると、&math({\rm rank}\,AB\le{\rm rank}\,A<n); となり、&math(AB); も正則でない。

したがって、正方行列 &math(A_1, A_2, \cdots, A_k); の積
&math(A_1 A_2 \cdots A_k); が正則であるなら、元の全ての行列が正則である。

すなわち、正方行列 &math(A_1, A_2, \cdots, A_k); について、次の2つの条件は同値である。
- &math(A_1, A_2, \cdots, A_k); がすべて正則
- $k$ 個の行列の積 &math(A_1 A_2 \cdots A_k); が正則

これらと、次の2つの条件とを比べてみると良い。
- &math(a_1, a_2, \cdots, a_k); がすべてゼロでない
- $k$ 個の数の積 &math(a_1 a_2 \cdots a_k); がゼロでない

* 行列の階数と連立一次方程式の解 [#r53ca169]

&math(\begin{bmatrix}A&\bm b\end{bmatrix}\sim\begin{bmatrix}PA&P\bm b\end{bmatrix});

で、右辺が階段行列になったとする。~
このとき、&math(PA); も階段行列化されていることに注意。

- &math({\rm rank}\,A<{\rm rank}\,\begin{bmatrix}A&\bm b\end{bmatrix}); ならば
最後の行まで掃き出せてしまい、解なしとなる
- &math({\rm rank}\,A={\rm rank}\,\begin{bmatrix}A&\bm b\end{bmatrix}); ならば
(A の列数)−(rank A) だけパラメータが入った一般解を得る

&math(\bm x=\bm x_0+k_1\bm x_1+k_2\bm x_2+k_3\bm x_3+\cdots);

** 解の自由度 [#rcc4e5a5]

一般解中のパラメータの数を「解の自由度」と呼ぶ。

&math(\text{自由度\,} &= \text{(A の列数)}−(\text{rank\,} A)\\
&=\text{(A の列数)}−\text{(掃き出しの行えた回数)}\\
&=\text{(掃き出しの行えなかった列数)}\\
);

ガウスの消去法では掃き出しの行えなかった列をパラメータに置くことで解の一般形が得られたことを思い出そう。

このように、解の自由度は &math(\bm b); によらず、&math(A); によってのみ決まる。

** 斉次方程式と非斉次方程式 [#pcfa889a]

斉次の「斉」は、一斉とか校歌斉唱とかの斉で、そろっているという意味。

斉次方程式とは、変数の次数がそろっているということで、
今の場合1次方程式なので、0次項がないということ。

- &math(A\bm x=\bm o); は斉次方程式
- &math(A\bm x=\bm b); は非斉次方程式 &math((\bm b\ne \bm o));

*** 斉次方程式の解 [#d82c1b4e]

斉次連立一次方程式 &math(A\bm x=\bm o); は「自明な解」&math(\bm x=\bm o); を持つ。

&math({\rm rank}\,A); が &math(A); の列数と等しければ、&math(A); は正則なので左から &math(A^{-1}); を掛けることで &math(\bm x=\bm o); が得られ、解はこれ以外に存在しない。

一方、&math({\rm rank}\,A); が &math(A); の列数より小さければ、自由度が1以上となり、自明な解以外にも解を持つ。

*** 非斉次方程式の解 [#q447dcd7]

斉次方程式 &math(A\bm x=\bm o); の一般解を &math(\bm x_0); とし、~
非斉次方程式 &math(A\bm x=\bm b); の一つの特殊解を &math(\bm x^*); とすると、~
&math(\bm x=\bm x^*+\bm x_0); は非斉次方程式 &math(A\bm x=\bm b); の一般解となる。

∵&math(A\bm x=\bm b=A\bm x^*); ならば、~
&math(A(\bm x-\bm x^*)=\bm o); だから、~
&math(\bm x=\bm x^*+\bm x_0);

すなわち、非斉次方程式の特殊解に斉次方程式の一般解を加えることで非斉次方程式の一般解を作れる。

* コメント [#y1cc64e2]

#article_kcaptcha
**行列の階数について [#v7882d2a]
>[[Kobe]] (&timetag(2021-10-01T08:29:31+09:00, 2021-10-01 (金) 17:29:31);)~
~
「n,lを与えられた自然数とするとき,n次行列Aが正則である必要十分条件~
はrank(AB)=rank(B)が任意の(n,l)型行列Bに対して成りたつことを示せ。」という問題に対して、十分条件を示すにはどのようにアプローチしていけばいいですか?~

//
- rank は行列の積によって減ることはない、つまり rank(AB) <= min(rank(A), rank(B)) であること、および、任意の行列 C に正則行列 P を掛けても rank が変化しないこと rank(PC)=rank(C)、は非常に重要な定理でほぼすべての教科書に載っているのではないかと思います。ただ、教科書によって証明の順番(どの定理をどの定理よりも先に証明するか)は異なりますので、この問題への解答としてはお使いの教科書の流れに沿ったアプローチが求められるのではないかと思います。 -- [[武内(管理人)]] &new{2021-10-05 (火) 12:07:03};

#comment_kcaptcha

**(補題) ゼロ行ベクトルを含む行列は正則行列になり得ない †の部分について [#h77bdbb4]
>[[rum]] (&timetag(2018-07-05T15:26:02+09:00, 2018-07-06 (金) 00:26:02);)~
~
(補題) ゼロ行ベクトルを含む行列は正則行列になり得ない †の部分についてなのですが、AB=()、B=()となっていますが、A=()、AB=()ではないでしょうか?私の理解が及んでおらず間違っていたらすみません。~

//
- わかりにくい表記ですみません。ここでは2つの式の間にカンマはなく、AB=[ ]B=[ ] のつもりでした。分かりやすくなるよう改訂しましたので参照して下さい。 -- [[武内(管理人)]] &new{2018-07-10 (火) 13:34:24};

#comment_kcaptcha

**定理2.11はPの一意性 [#b41d38a8]
>[[Betsumiya]] (&timetag(2015-08-06T05:36:47+09:00, 2015-08-06 (木) 14:36:47);)~
~
定理2.11はPの一意性までは主張していないと思います。~
P.S. いつも参考にさせて頂いています。~

//
- コメントありがとうございます。おっしゃる通りですね。簡単な反例として、&math(A=O); はそのまま階段行列ですが、この場合任意の正則行列を &math(P); として採用できます。本文の誤りを修正いたしました。 -- [[武内(管理人)]] &new{2015-08-06 (木) 15:33:10};

#comment_kcaptcha


Counter: 93757 (from 2010/06/03), today: 13, yesterday: 0