跳到主要内容

大阪大学 基礎工学研究科 数理科学 (システム創成専攻) 2019年8月実施 数理科学 II [2]

Author​

祭音Myyura (co-authored with GPT 5.6 SOL)

Description​

m≥nm\ge n は自然数、b∈Rmb\in\mathbb R^m とする。AA は rank⁡A=n\operatorname{rank}A=n を満たす m×nm\times n 実行列で、W={Ax∣x∈Rn}W=\{Ax\mid x\in\mathbb R^n\} とおく。内積は標準内積とする。

(1) y∈Rny\in\mathbb R^n とする。b−Ay⊥Wb-Ay\perp W と ATb=ATAyA^Tb=A^TAy が同値であることを示せ。

(2) x∈Rnx\in\mathbb R^n とする。xTATAx=0x^TA^TAx=0 ならば x=0x=0 を示せ。

(3) ATAA^TA が正則であることを示せ。

(4) bb の WW への直交射影は A(ATA)−1ATbA(A^TA)^{-1}A^Tb であることを示せ。

Kai​

(1)​

b−Ay⊥W  ⟺  xTAT(b−Ay)=0 (∀x∈Rn)  ⟺  ATb=ATAy.b-Ay\perp W\iff x^TA^T(b-Ay)=0\ (\forall x\in\mathbb R^n) \iff A^Tb=A^TAy.

(2)​

xTATAx=∥Ax∥2x^TA^TAx=\|Ax\|^2 なので Ax=0Ax=0。列が一次独立であるから x=0x=0。

(3)​

ATAx=0A^TAx=0 なら (2) より x=0x=0。したがって核は {0}\{0\} であり正則である。

(4)​

y=(ATA)−1ATby=(A^TA)^{-1}A^Tb とおくと u=Ay∈Wu=Ay\in W で、(1) より v=b−u⊥Wv=b-u\perp W。よって

proj⁡Wb=A(ATA)−1ATb.\boxed{\operatorname{proj}_Wb=A(A^TA)^{-1}A^Tb}.