ここで行き詰まった? AI で解決しよう
Algebra · real student question

A を m 行 n 列の行列とする。ABA = A を満たす n 行 m 列の行列 B が存在することを示せ。

Question

AA を m×nm\times n 行列とする。

ABA=AABA=A

を満たす n×mn\times m 行列 BB が存在することを示せ。

Step-by-step solution

  1. 標準形に帰着させる。 r=rank⁡(A)r=\operatorname{rank}(A) とおく。どの行列もその階数標準形と同値である。すなわち正則な PP(m×mm\times m)と QQ(n×nn\times n)が存在して

    PAQ=D:=(Ir000)m×nPAQ=D:=\begin{pmatrix}I_{r}&0\\ 0&0\end{pmatrix}_{m\times n}

    となる。PP と QQ は、AA をこのブロック形まで掃き出す行基本変形と列基本変形を記録したものである。

  2. まず DD について問題を解く。 n×mn\times m 行列

    D−:=(Ir000)n×mD^{-}:=\begin{pmatrix}I_{r}&0\\ 0&0\end{pmatrix}_{n\times m}

    を定める。同じブロックの形で寸法だけ入れ替えたものであり、これによって積の型がそろう。ブロックごとの掛け算から DD−=(Ir000)m×mDD^{-}=\begin{pmatrix}I_{r}&0\\0&0\end{pmatrix}_{m\times m} となり、したがって

    DD−D=DDD^{-}D=D

    が成り立つ。この冪等な射影を DD に左から掛けても、唯一の非零ブロックはそのまま残るからである。

  3. 解を AA に戻す。 PAQ=DPAQ=D から A=P−1DQ−1A=P^{-1}DQ^{-1} が得られる。ここで

    B:=QD−PB:=QD^{-}P

    と定めると、これは (n×n)(n×m)(m×m)=n×m(n\times n)(n\times m)(m\times m)=n\times m であり、まさに求める型である。B=QD−PB=QD^{-}P であって QDPQDP ではないことに注意する。寸法を合わせるには、型を入れ替えたブロック行列が不可欠である。

  4. 代入して ABA=AABA=A を確かめる。 両方の表示を代入し、逆行列の対を消していく。

    ABA=(P−1DQ−1)(QD−P)(P−1DQ−1)=P−1D D−D Q−1=P−1DQ−1=AABA=\left(P^{-1}DQ^{-1}\right)\left(QD^{-}P\right)\left(P^{-1}DQ^{-1}\right)=P^{-1}D\,D^{-}D\,Q^{-1}=P^{-1}DQ^{-1}=A

    ここで Q−1Q=InQ^{-1}Q=I_{n}、PP−1=ImPP^{-1}=I_{m}、およびステップ 2 の等式 DD−D=DDD^{-}D=D を用いた。これで証明が完成する。

  5. この主張が言っていること・言っていないことを確認する。 BB は AA の一般化逆行列(あるいは {1}\{1\} 逆行列)と呼ばれる。一般に一意ではない。P,QP,Q の選び方を変えれば QD−PQD^{-}P の形の BB はどれも条件を満たし、r<min⁡(m,n)r<\min(m,n) のときは自由なブロックのぶんだけ無数に存在する。これが A−1A^{-1} と一致するのは、AA が正方かつ正則なとき(r=m=nr=m=n)だけである。

  6. 数値による確認。 ムーア・ペンローズの擬似逆行列は、そのような BB の 1 つの具体例である。4×34\times 3(階数 2)、3×53\times 5(階数 3)、5×55\times 5(階数 1)の形のランダムな行列について A A+A−AA\,A^{+}A-A を計算すると、成分の大きさは 10−1510^{-15} 以下になる。定理が予言するとおり、これは機械的な零である。

Answer

B=QD−P, ただし PAQ=D=(Ir000)m×n かつ D−=(Ir000)n×mB=QD^{-}P,\ \text{ただし } PAQ=D=\begin{pmatrix}I_{r}&0\\ 0&0\end{pmatrix}_{m\times n}\ \text{かつ}\ D^{-}=\begin{pmatrix}I_{r}&0\\ 0&0\end{pmatrix}_{n\times m}

Need to solve a different problem like this? Open the solver →