logo

線形代数における射影 📂行列代数

線形代数における射影

定義

正方行列 PCm×mP \in \mathbb{C}^{m \times m}P2=PP^2 = P であれば、射影作用素projectorという。

説明

代数学の用語では、冪等元idempotentと表現し、同様に a2=aa^2 = a のような元を指す。

PP が射影であるなら、(IP)2=I2P+P2=I2P+P=(IP)(I-P)^2 = I - 2P + P^2 = I - 2P + P = (I-P) であり、従って (IP)(I-P) も射影であることがわかる。

このような射影作用素 (IP)(I - P)PP補射影作用素complementary Projectorと呼ぶ。

射影を幾何学的に考えると、空間図形に光を当ててその影を得ることである。例えば、f(x,y,z):=(x,y,0)f(x,y,z) := (x,y,0) のような関数は zz 軸の方向に光を投げて、xyxy 平面にできる影を表す。その影を再度射影しても結果は同じで、その意味で P2=PP^2 = P が射影作用素である定義は理にかなっていると言える。

性質

射影 PCm×mP \in \mathbb{C}^{m \times m} とその補射影 IPI-P は次の性質を満たす。

(a) C(IP)=N(P)\mathcal{C} (I-P) = \mathcal{N} (P)

(b) N(1P)=C(P)\mathcal{N} (1-P) = \mathcal{C} (P)

(c) N(1P)N(P)={0}\mathcal{N} (1-P) \cap \mathcal{N} (P) = \left\{ 0 \right\}

(d) C(P)N(P)={0}\mathcal{C} (P) \cap \mathcal{N} (P) = \left\{ 0 \right\}

(e) C(P)N(P)=Cm\mathcal{C} (P) \oplus \mathcal{N} (P) = \mathbb{C}^{m}

証明

(a)(b)

C(IP)\mathcal{C} ( I - P)N(P)\mathcal{N}(P) が互いに含まれることを示せばよい。任意のベクトル vCm\mathbf{v} \in \mathbb{C}^{m} に対して、

(IP)v=vPv (I - P) \mathbf{v} = \mathbf{v} - P \mathbf{v}

もし vN(P)\mathbf{v} \in \mathcal{N} (P) ならば、Pv=0P \mathbf{v} = \mathbb{0} となり、従って (IP)v=v(I - P) \mathbf{v} = \mathbf{v}、つまり vC(IP)\mathbf{v} \in \mathcal{C} (I - P) であり、

N(P)C(IP) \mathcal{N} (P) \subset \mathcal{C} (I - P)

wC(IP)\mathbf{w} \in \mathcal{C} (I-P) とすると、w=(IP)v\mathbf{w} = (I - P) \mathbf{v} を満たす v\mathbf{v}C(IP)\mathcal{C} (I-P) に存在する。w=(IP)v\mathbf{w} = (I - P) \mathbf{v} に射影 PP を適用すると、

Pw=PvP2v=PvPv=0 P \mathbf{w} = P \mathbf{v} - P^2 \mathbf{v} = P \mathbf{v} - P \mathbf{v} = \mathbb{0}

つまり wN(P)\mathbf{w} \in \mathcal{N} (P) であり、従って、

C(IP)N(P) \mathcal{C} (I - P) \subset \mathcal{N} (P)

これにより (1) が証明され、P=I(IP)P = I - (I- P) であるから、射影 PP の補射影 (IP)(I - P) について考えると、すぐに (2) が証明される。

(c)(d)

vN(IP)N(P)\mathbf{v} \in \mathcal{N} (I - P) \cap \mathcal{N} (P) が零ベクトルではないと仮定する。しかし、vN(IP)\mathbf{v} \in \mathcal{N} (I - P) であるから、(IP)v=0(I-P) \mathbf{v} = \mathbb{0} であり、vN(P)\mathbf{v} \in \mathcal{N} (P) であるから、

Pv=0 P \mathbf{v} = \mathbb{0}

両辺を足せば (IP)v+Pv=v=0(I-P) \mathbf{v} + P \mathbf{v} = \mathbf{v} = \mathbb{0} となり、これは仮定に矛盾する。

これにより (3) が証明され、12 によってすぐに 4 も証明される。

(e)

直和の定義に従い、存在性、排他性、一意性を示せばよい。

  • (i) 存在性

    (a) により N(P)=C(IP)\mathcal{N} (P) = \mathcal{C} (I - P) であり、任意のベクトル sCm\mathbf{s} \in \mathbb{C}^{m} に対して、PsC(P)P \mathbf{s} \in \mathcal{C}(P) そして (IP)sC(IP)(I - P)\mathbf{s} \in \mathcal{C} (I - P) である。

一方、Ps+(IP)s=Ps+sPs=sCmP \mathbf{s} + (I - P) \mathbf{s} = P \mathbf{s} + \mathbf{s} - P \mathbf{s} = s \in \mathbb{C}^{m} であるから、s\mathbf{s} は常に C(P)\mathcal{C}(P)C(IP)\mathcal{C} (I - P) の和として表すことができる。

  • (ii) 排他性

    既に (d) で証明されている。

  • (iii) 一意性

    上記の (ii) により、sCms \in \mathbb{C}^{m} に対して、

    s=c1+n1=c2+n2 \mathbf{s} = \mathbf{c}_{1} + \mathbf{n}_{1} = \mathbf{c}_{2} + \mathbf{n}_{2}

    を満たす c1,c2C(P)\mathbf{c}_{1} , \mathbf{c}_{2} \in \mathcal{C}(P)n1,n2N(P)\mathbf{n}_{1}, \mathbf{n}_{2} \in \mathcal{N}(P) が存在する。

ここで c1c2\mathbf{c}_{1} \ne \mathbf{c}_{2} n1n2\mathbf{n}_{1} \ne \mathbf{n}_{2} と仮定する。

c1+n1=c2+n2\mathbf{c}_{1} + \mathbf{n}_{1} = \mathbf{c}_{2} + \mathbf{n}_{2} の両辺に PP を乗じると、

Pc1+Pn1=Pc2+Pn2 P\mathbf{c}_{1} + P\mathbf{n}_{1} = P\mathbf{c}_{2} + P\mathbf{n}_{2}

一方、n1,n2N(P)\mathbf{n}_{1}, \mathbf{n}_{2} \in \mathcal{N} (P) であるから、

Pc1=Pc2 P\mathbf{c}_{1} = P\mathbf{c}_{2}

つまり P(c1c2)=0P( \mathbf{c}_{1} - \mathbf{c}_{2} ) = \mathbb{0}。零空間の定義から、(c1c2)N(P)( \mathbf{c}_{1} - \mathbf{c}_{2}) \in \mathcal{N}(P) であり、ベクトル空間の性質から、(c1c2)C(P)( \mathbf{c}_{1} - \mathbf{c}_{2}) \in \mathcal{C}(P) だが、(排他性) により、c1c2=0\mathbf{c}_{1} - \mathbf{c}_{2} = \mathbb{0} でなければならない。

これは c1c2\mathbf{c}_{1} \ne \mathbf{c}_{2} という仮定と矛盾するし、同様に n1=n2\mathbf{n}_{1} = \mathbf{n}_{2} も証明できる。

参考