階数因数分解

提供: testwiki
ナビゲーションに移動 検索に移動

階数因数分解(かいすういんすうぶんかい、テンプレート:Lang-en-short)あるいは階数分解(rank decomposition)とは、数学線型代数学の分野において、階数r のある与えられた m×n 行列 A のある m×r 行列 Cr×n 行列 F の積としての表示 A=CF のことを言う。

全ての有限次元行列には階数因数分解が存在するA を、列階数r であるような m×n 行列とする。すなわち、A には r 個の線型独立な列が含まれる。あるいは同じ意味であるが、A列空間次元r である。c1,c2,,cr を、A の列空間の任意の基底とし、それらを列ベクトルとして m×r 行列 C=[c1:c2::cr] を構成する。したがって、A の全ての列ベクトルは、C の列の線型結合である。正確に言うと、A=[a1:a2::an] を第 j 列が aj であるような m×n 行列とすれば、

aj=f1jc1+f2jc2++frjcr,

となる。ただし fij は、基底 c1,c2,,cr に関する aj のスカラー係数である。このことは、fij(i,j)-成分とする行列 F によって A=CF が得られることを意味する。

rank(A) = rank(AT)

階数因数分解から直ちに従う帰結として、A の階数はその転置行列 AT の階数と等しい、というものがある。すると A の列は AT の行であることから、A列階数行階数は等しいことが分かる。

証明:これが真であることを示すために、はじめに行列の「階数」とはその「列階数」を意味するものであると定義しておく。A=CF より、AT=FTCT が従う。テンプレート:仮リンクの定義から、この等式は AT の各列が FT の列の線型結合であることを意味する。したがって、AT の列空間は FT の列空間に含まれるものであることが分かり、したがって rank(AT) ≤ rank(FT) が成立する。今 FTn×r 行列であるので、FT には r 個の列が存在し、したがって rank(AT) ≤ r = rank(A) が成立する。これより rank(AT) ≤ rank(A) が示された。続いて、その逆の不等式が成立することを示すために、AT に対して上述の結果を適用する。(AT)T = A なので、rank(A) = rank((AT)T) ≤ rank(AT) と書くことが出来る。このことから rank(A) ≤ rank(AT) が示される。したがって、rank(AT) ≤ rank(A) かつ rank(A) ≤ rank(AT) であることから、rank(A) = rank(AT) が示された。

行階段形からの階数因数分解

実際、特定の階数因数分解を次の手順で構成することが出来る:A行既約階段形 B は計算することで得られる。このとき、上述の行列 CA から全ての非ピボット列を除くことで得られ、FB から全てのゼロ行を除くことで得られる。

行列

A=[1314273915311208][1020011000010000]=B

を考える。B は既約階段形である。このとき、C は、A の唯一つの非ピボット列である第三列を除くことで得られ、F は最後のゼロ行を除くことで得られる。すなわち、

C=[134279151128],F=[102001100001]

が得られる。次の関係式が、直ちに確かめられる:

A=[1314273915311208]=[134279151128][102001100001]=CF.

証明

P を、ブロック区分けの形式で AP=(C,D) が成立するような n×n 置換行列とする。ただし C はその列が Ar 個のピボット列であるものとする。D の全ての列は C の列の線形結合であり、したがって D=CG が成立するような行列 G が存在する。ただし G の列は、それら各線形結合の係数を含むものである。すると、r×r 単位行列 Ir によって AP=(C,CG)=C(Ir,G) と書くことが出来る。続いて、(Ir,G)=FP を証明する。

AP を、基本行列の積であるような行列 E を左から掛けることによって、行既約階段形へと変換する。すなわち、EAP=BP=EC(Ir,G) が得られる。ただし EC=[Ir0] である。すると、BP=[IrG00] と書くことが出来、したがって (Ir,G)=FP ということが分かる。これはすなわち、A に対して行ったものと同様の置換を列に対して行うことで得られる、行既約階段形に含まれる非ゼロの r 個の行である。したがって、AP=CFP が得られ、P が可逆であることから A=CF が得られる。以上で証明は完成された。

参考文献

テンプレート:Refbegin

テンプレート:Refend