考研相關文章參考資料為 wjungle 大神提供的筆記

Ch7 內積空間

內積

定義

設 VV 為 vector space over FF。
⟨⋅,⋅⟩:V×V→F\langle \cdot,\cdot\rangle:V\times V\to F 為 function,且對所有 u⃗,v⃗,w⃗∈V\vec{u},\vec{v},\vec{w}\in V,c,d∈Fc,d\in F,滿足:

  • ⟨u⃗+v⃗,w⃗⟩=⟨u⃗,w⃗⟩+⟨v⃗,w⃗⟩\langle \vec{u}+\vec{v},\vec{w}\rangle=\langle \vec{u},\vec{w}\rangle+\langle \vec{v},\vec{w}\rangle。
  • ⟨cu⃗,v⃗⟩=c⟨u⃗,v⃗⟩\langle c\vec{u},\vec{v}\rangle=c\langle \vec{u},\vec{v}\rangle。
  • ⟨u⃗,v⃗⟩=⟨v⃗,u⃗⟩‾\langle \vec{u},\vec{v}\rangle=\overline{\langle \vec{v},\vec{u}\rangle}。
  • 若 v⃗≠0⃗\vec{v}\ne \vec{0},則 ⟨v⃗,v⃗⟩>0\langle \vec{v},\vec{v}\rangle>0。

稱 ⟨⋅,⋅⟩\langle \cdot,\cdot\rangle 為 VV 上之 inner product。

Note
  • 左線性: ⟨cu⃗+dv⃗,w⃗⟩=c⟨u⃗,w⃗⟩+d⟨v⃗,w⃗⟩\langle c\vec{u}+d\vec{v},\vec{w}\rangle=c\langle \vec{u},\vec{w}\rangle+d\langle \vec{v},\vec{w}\rangle
  • 右共軛線性: ⟨u⃗,cv⃗+dw⃗⟩=⟨cv⃗+dw⃗,u⃗⟩‾=cˉ⟨u⃗,v⃗⟩+dˉ⟨u⃗,w⃗⟩\begin{aligned}\langle \vec{u},c\vec{v}+d\vec{w}\rangle=\overline{\langle c\vec{v}+d\vec{w},\vec{u}\rangle}=\bar{c}\langle \vec{u},\vec{v}\rangle+\bar{d}\langle \vec{u},\vec{w}\rangle\end{aligned}

常見內積

  • V=FnV=F^n (standard inner product): ⟨x⃗,y⃗⟩=x1yˉ1+x2yˉ2+⋯+xnyˉn\langle \vec{x},\vec{y}\rangle=x_1\bar{y}_1+x_2\bar{y}_2+\cdots+x_n\bar{y}_n

    • matrix form: ⟨x⃗,y⃗⟩=y⃗Hx⃗in Fn×1\langle \vec{x},\vec{y}\rangle=\vec{y}^{H}\vec{x}\quad \text{in }F^{n\times 1}
    • matrix form: ⟨x⃗,y⃗⟩=x⃗y⃗Hin F1×n\langle \vec{x},\vec{y}\rangle=\vec{x}\vec{y}^{H}\quad \text{in }F^{1\times n}
  • V=Fn×nV=F^{n\times n}: ⟨A,B⟩=∑i=1n∑j=1naijbˉij\langle A,B\rangle=\sum_{i=1}^{n}\sum_{j=1}^{n}a_{ij}\bar{b}_{ij}

    • matrix form: ⟨A,B⟩=tr⁡(ABH)\langle A,B\rangle=\operatorname{tr}(AB^H)
  • V=C[a,b]V=C[a,b]: ⟨f,g⟩=∫abf(x)gˉ(x) dx\langle f,g\rangle=\int_a^b f(x)\bar{g}(x)\,dx

  • V=PnV=P_n,x1,x2,…,xn+1x_1,x_2,\ldots,x_{n+1} 相異: ⟨P,Q⟩=∑i=1n+1P(xi)Q(xi)‾\langle P,Q\rangle=\sum_{i=1}^{n+1}P(x_i)\overline{Q(x_i)}

    • 若取點數不夠,可能不成 inner product。例如 V=P2V=P_2,只取 x0=0,x1=1x_0=0,x_1=1,令 P=x(x−1)≠0P=x(x-1)\ne 0,則 ⟨P,P⟩=P(0)P(0)+P(1)P(1)=0\langle P,P\rangle=P(0)P(0)+P(1)P(1)=0 不滿足正定性。

性質

設 ⟨⋅,⋅⟩\langle \cdot,\cdot\rangle 為 VV 上的 inner product。

  • ⟨v⃗,0⃗⟩=0=⟨0⃗,v⃗⟩\langle \vec{v},\vec{0}\rangle=0=\langle \vec{0},\vec{v}\rangle。
  • ⟨cu⃗,v⃗⟩=c⟨u⃗,v⃗⟩\langle c\vec{u},\vec{v}\rangle=c\langle \vec{u},\vec{v}\rangle。
  • ⟨∑iciu⃗i,v⃗⟩=∑ici⟨u⃗i,v⃗⟩\left\langle \sum_i c_i\vec{u}_i,\vec{v}\right\rangle=\sum_i c_i\langle \vec{u}_i,\vec{v}\rangle。
  • ⟨u⃗,cv⃗⟩=cˉ⟨u⃗,v⃗⟩\langle \vec{u},c\vec{v}\rangle=\bar{c}\langle \vec{u},\vec{v}\rangle。
  • ⟨u⃗,∑iciv⃗i⟩=∑icˉi⟨u⃗,v⃗i⟩\left\langle \vec{u},\sum_i c_i\vec{v}_i\right\rangle=\sum_i\bar{c}_i\langle \vec{u},\vec{v}_i\rangle。
  • ⟨v⃗,v⃗⟩≥0\langle \vec{v},\vec{v}\rangle\ge 0,且 ⟨v⃗,v⃗⟩=0⟺v⃗=0⃗\langle \vec{v},\vec{v}\rangle=0\Longleftrightarrow\vec{v}=\vec{0}。
  • ⟨u⃗,v⃗⟩=0\langle \vec{u},\vec{v}\rangle=0,u⃗,v⃗∈V⟺u⃗=0⃗\vec{u},\vec{v}\in V\Longleftrightarrow\vec{u}=\vec{0} 或 v⃗=0⃗\vec{v}=\vec{0}。
  • 若 β={v⃗1,…,v⃗n}\beta=\{\vec{v}_1,\ldots,\vec{v}_n\} 為 VV 的 basis,且 ⟨u⃗,v⃗i⟩=0\langle \vec{u},\vec{v}_i\rangle=0,i=1,…,ni=1,\ldots,n,則 u⃗=0⃗\vec{u}=\vec{0}。
定義

設 ⟨⋅,⋅⟩\langle \cdot,\cdot\rangle 為 VV 上的 inner product。

  • ∥v⃗∥=⟨v⃗,v⃗⟩\lVert\vec{v}\rVert=\sqrt{\langle \vec{v},\vec{v}\rangle} 稱為 v⃗\vec{v} 的 norm。若 ∥v⃗∥=1\lVert\vec{v}\rVert=1,稱 v⃗\vec{v} 為 unit vector。
  • 若 ⟨u⃗,v⃗⟩=0\langle \vec{u},\vec{v}\rangle=0,稱 u⃗,v⃗\vec{u},\vec{v} 為 orthogonal,記作 u⃗⊥v⃗\vec{u}\perp\vec{v}。
  • S={v⃗1,…,v⃗k}S=\{\vec{v}_1,\ldots,\vec{v}_k\},若 v⃗i⊥v⃗j\vec{v}_i\perp\vec{v}_j,∀i≠j\forall i\ne j,稱 SS 為 orthogonal set,允許 0⃗\vec{0}。
  • 若 SS 為 orthogonal set 且 ∥v⃗i∥=1\lVert\vec{v}_i\rVert=1,i=1,…,ki=1,\ldots,k,稱 SS 為 orthonormal set,不允許 0⃗\vec{0}。
  • u⃗\vec{u} 與 v⃗\vec{v} 的 distance 為 d(u⃗,v⃗)=∥u⃗−v⃗∥d(\vec{u},\vec{v})=\lVert\vec{u}-\vec{v}\rVert。
Note

若 S={v⃗1,…,v⃗k}S=\{\vec{v}_1,\ldots,\vec{v}_k\},則

S 為 orthonormal set⟺⟨v⃗i,v⃗j⟩=δij={1,i=j,0,i≠jS\text{ 為 orthonormal set} \Longleftrightarrow \langle \vec{v}_i,\vec{v}_j\rangle=\delta_{ij}= \begin{cases} 1,& i=j,\\ 0,& i\ne j \end{cases}


Note

在 C[−π,π]C[-\pi,\pi] 上定義 ⟨f,g⟩=1π∫−ππf(x)g(x)‾ dx.\langle f,g\rangle=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)\overline{g(x)}\,dx.

則 S={12,cos⁡x,cos⁡2x,…}S=\left\{\frac{1}{\sqrt{2}},\cos x,\cos 2x,\ldots\right\} 為 orthonormal set。

Note

在 FnF^n 的 standard inner product 下:

∥x⃗∥=⟨x⃗,x⃗⟩=(∣x1∣2+⋯+∣xn∣2)1/2\lVert\vec{x}\rVert=\sqrt{\langle\vec{x},\vec{x}\rangle}=\left(|x_1|^2+\cdots+|x_n|^2\right)^{1/2} 稱為 22-norm。

也可定義:

  • ∥x⃗∥1=∣x1∣+⋯+∣xn∣\lVert\vec{x}\rVert_1=|x_1|+\cdots+|x_n|
  • ∥x⃗∥∞=max⁡{∣x1∣,…,∣xn∣}\lVert\vec{x}\rVert_\infty=\max\{|x_1|,\ldots,|x_n|\}
  • ∥x⃗∥p=(∣x1∣p+⋯+∣xn∣p)1/p\lVert\vec{x}\rVert_p=\left(|x_1|^p+\cdots+|x_n|^p\right)^{1/p}

norm 需滿足:

  • ∥x⃗∥≥0\lVert\vec{x}\rVert\ge 0,且 ∥x⃗∥=0⟺x⃗=0⃗\lVert\vec{x}\rVert=0\Longleftrightarrow\vec{x}=\vec{0}。
  • ∥cx⃗∥=∣c∣∥x⃗∥\lVert c\vec{x}\rVert=|c|\lVert\vec{x}\rVert。
  • ∥x⃗+y⃗∥≤∥x⃗∥+∥y⃗∥\lVert\vec{x}+\vec{y}\rVert\le\lVert\vec{x}\rVert+\lVert\vec{y}\rVert。

定理

Pythagorean theorem:
設 ⟨⋅,⋅⟩\langle\cdot,\cdot\rangle 為 VV 上的 inner product。若 u⃗⊥v⃗\vec{u}\perp\vec{v},則 ∥u⃗+v⃗∥2=∥u⃗∥2+∥v⃗∥2. \lVert\vec{u}+\vec{v}\rVert^2=\lVert\vec{u}\rVert^2+\lVert\vec{v}\rVert^2.

Note

若 {v⃗1,…,v⃗k}\{\vec{v}_1,\ldots,\vec{v}_k\} 為 orthogonal set,則

∥c1v⃗1+⋯+ckv⃗k∥2=∥c1v⃗1∥2+⋯+∥ckv⃗k∥2=∣c1∣2∥v⃗1∥2+⋯+∣ck∣2∥v⃗k∥2\left\lVert c_1\vec{v}_1+\cdots+c_k\vec{v}_k\right\rVert^2 =\lVert c_1\vec{v}_1\rVert^2+\cdots+\lVert c_k\vec{v}_k\rVert^2 =|c_1|^2\lVert\vec{v}_1\rVert^2+\cdots+|c_k|^2\lVert\vec{v}_k\rVert^2

定理

設 ⟨⋅,⋅⟩\langle\cdot,\cdot\rangle 為 VV 上的 inner product,則對所有 u⃗,v⃗∈V\vec{u},\vec{v}\in V:

  • Cauchy-Schwarz inequality: $ |\langle\vec{u},\vec{v}\rangle|\le\lVert\vec{u}\rVert\lVert\vec{v}\rVert.$

  • Triangle inequality: ∥u⃗+v⃗∥≤∥u⃗∥+∥v⃗∥.\lVert\vec{u}+\vec{v}\rVert\le\lVert\vec{u}\rVert+\lVert\vec{v}\rVert.

證明
  • 若 v⃗=0⃗\vec{v}=\vec{0},結論顯然。若 v⃗≠0⃗\vec{v}\ne\vec{0},取 α=⟨u⃗,v⃗⟩⟨v⃗,v⃗⟩.\alpha=\frac{\langle\vec{u},\vec{v}\rangle}{\langle\vec{v},\vec{v}\rangle}.
    由 0≤∥u⃗−αv⃗∥20\le\lVert\vec{u}-\alpha\vec{v}\rVert^2,可得 ∣⟨u⃗,v⃗⟩∣2≤∥u⃗∥2∥v⃗∥2.|\langle\vec{u},\vec{v}\rangle|^2\le\lVert\vec{u}\rVert^2\lVert\vec{v}\rVert^2.

  • 展開後套用 Cauchy-Schwarz inequality:
    ∥u⃗+v⃗∥2=∥u⃗∥2+2Re⁡⟨u⃗,v⃗⟩+∥v⃗∥2≤∥u⃗∥2+2∥u⃗∥∥v⃗∥+∥v⃗∥2=(∥u⃗∥+∥v⃗∥)2.\lVert\vec{u}+\vec{v}\rVert^2=\lVert\vec{u}\rVert^2+2\operatorname{Re}\langle\vec{u},\vec{v}\rangle+\lVert\vec{v}\rVert^2\le\lVert\vec{u}\rVert^2+2\lVert\vec{u}\rVert\lVert\vec{v}\rVert+\lVert\vec{v}\rVert^2=(\lVert\vec{u}\rVert+\lVert\vec{v}\rVert)^2.

Note

在 C[a,b]C[a,b] 上,若 ⟨f,g⟩=∫abf(x)g(x) dx\langle f,g\rangle=\int_a^b f(x)g(x)\,dx
則 Cauchy-Schwarz inequality 為

∣∫abf(x)g(x) dx∣≤(∫abf(x)2 dx)1/2(∫abg(x)2 dx)1/2\left|\int_a^b f(x)g(x)\,dx\right| \le\left(\int_a^b f(x)^2\,dx\right)^{1/2} \left(\int_a^b g(x)^2\,dx\right)^{1/2}

Note

若 u⃗,v⃗≠0⃗\vec{u},\vec{v}\ne\vec{0},由 Cauchy-Schwarz inequality 可得

−1≤⟨u⃗,v⃗⟩∥u⃗∥∥v⃗∥≤1-1\le\frac{\langle\vec{u},\vec{v}\rangle}{\lVert\vec{u}\rVert\lVert\vec{v}\rVert}\le1

故存在 θ∈[0,π]\theta\in[0,\pi],使得

cos⁡θ=⟨u⃗,v⃗⟩∥u⃗∥∥v⃗∥\cos\theta=\frac{\langle\vec{u},\vec{v}\rangle}{\lVert\vec{u}\rVert\lVert\vec{v}\rVert}

稱 θ\theta 為 u⃗\vec{u} 與 v⃗\vec{v} 之間的角度。


定理

設 S={v⃗1,…,v⃗k}S=\{\vec{v}_1,\ldots,\vec{v}_k\} 為 nonzero orthogonal set,且 v⃗∈span⁡(S)\vec{v}\in\operatorname{span}(S)。若 v⃗=c1v⃗1+⋯+ckv⃗k\vec{v}=c_1\vec{v}_1+\cdots+c_k\vec{v}_k 則

ci=⟨v⃗,v⃗i⟩⟨v⃗i,v⃗i⟩,i=1,…,kc_i=\frac{\langle\vec{v},\vec{v}_i\rangle}{\langle\vec{v}_i,\vec{v}_i\rangle},\quad i=1,\ldots,k

稱 cic_i 為 Fourier coefficient。

若 SS 為 orthonormal set,則 ci=⟨v⃗,v⃗i⟩c_i=\langle\vec{v},\vec{v}_i\rangle

定理

設 S={v⃗1,…,v⃗k}S=\{\vec{v}_1,\ldots,\vec{v}_k\} 為 nonzero orthogonal set,則 SS 為 LI。

Gram-Schmidt 正交化

定理

Gram-Schmidt orthogonalization Process:
給定 {v⃗1,…,v⃗n}\{\vec{v}_1,\ldots,\vec{v}_n\} 為 LI。令

u⃗1=v⃗1,u⃗i=v⃗i−∑j=1i−1⟨v⃗i,u⃗j⟩⟨u⃗j,u⃗j⟩u⃗j,i=2,…,n\begin{aligned} \vec{u}_1&=\vec{v}_1,\\ \vec{u}_i&=\vec{v}_i-\sum_{j=1}^{i-1} \frac{\langle\vec{v}_i,\vec{u}_j\rangle} {\langle\vec{u}_j,\vec{u}_j\rangle}\vec{u}_j,\quad i=2,\ldots,n \end{aligned}

則 {u⃗1,…,u⃗n}\{\vec{u}_1,\ldots,\vec{u}_n\} 為 nonzero orthogonal set,且 span⁡{u⃗1,…,u⃗n}=span⁡{v⃗1,…,v⃗n}\operatorname{span}\{\vec{u}_1,\ldots,\vec{u}_n\}=\operatorname{span}\{\vec{v}_1,\ldots,\vec{v}_n\}

Note

若 u⃗≠0⃗\vec{u}\ne\vec{0},令 z⃗=u⃗∥u⃗∥\vec{z}=\frac{\vec{u}}{\lVert\vec{u}\rVert}

稱為 normalized vector,且 ∥z⃗∥=1\lVert\vec{z}\rVert=1。
因此 {u⃗1∥u⃗1∥,…,u⃗n∥u⃗n∥}\left\{\frac{\vec{u}_1}{\lVert\vec{u}_1\rVert},\ldots,\frac{\vec{u}_n}{\lVert\vec{u}_n\rVert}\right\} 為 orthonormal set。


定理

QR-decomposition:

若 A∈Fm×nA\in F^{m\times n} 的 columns 為 LI,則 A=QRA=QR

其中 Q∈Fm×nQ\in F^{m\times n} 的 columns 為 orthonormal,R∈Fn×nR\in F^{n\times n} 為 upper triangular 且可逆。

正交投影

定義

設 W⊆VW\subseteq V,且 v⃗∈V\vec{v}\in V。若存在 v⃗0∈W\vec{v}_0\in W,使得 ⟨v⃗−v⃗0,w⃗⟩=0,∀w⃗∈W\langle\vec{v}-\vec{v}_0,\vec{w}\rangle=0,\quad\forall\vec{w}\in W

即 v⃗−v⃗0⊥W\vec{v}-\vec{v}_0\perp W,稱 v⃗0\vec{v}_0 為 v⃗\vec{v} on WW 的 orthogonal projection vector。

定理

設 W⊆VW\subseteq V,v⃗∈V\vec{v}\in V,且 {v⃗1,…,v⃗k}\{\vec{v}_1,\ldots,\vec{v}_k\} 為 WW 的 orthogonal basis。令

v⃗0=∑i=1k⟨v⃗,v⃗i⟩⟨v⃗i,v⃗i⟩v⃗i\vec{v}_0=\sum_{i=1}^k \frac{\langle\vec{v},\vec{v}_i\rangle} {\langle\vec{v}_i,\vec{v}_i\rangle}\vec{v}_i

則 v⃗0\vec{v}_0 為 v⃗\vec{v} 在 WW 上唯一的 orthogonal projection vector,記作 proj⁡Wv⃗\operatorname{proj}_W\vec{v}。

當 W=span⁡{v⃗}W=\operatorname{span}\{\vec{v}\} 時:

proj⁡Wu⃗=⟨u⃗,v⃗⟩⟨v⃗,v⃗⟩v⃗\operatorname{proj}_W\vec{u} =\frac{\langle\vec{u},\vec{v}\rangle}{\langle\vec{v},\vec{v}\rangle}\vec{v}

Note
  • v⃗∈W⟺proj⁡Wv⃗=v⃗\vec{v}\in W\Longleftrightarrow\operatorname{proj}_W\vec{v}=\vec{v}。
  • v⃗=0⃗⟹proj⁡Wv⃗=0⃗\vec{v}=\vec{0}\Longrightarrow\operatorname{proj}_W\vec{v}=\vec{0}。
  • v⃗⊥W⟺proj⁡Wv⃗=0⃗\vec{v}\perp W\Longleftrightarrow\operatorname{proj}_W\vec{v}=\vec{0}。

定義

設 W⊆VW\subseteq V。定義 P:V→V,P(v⃗)=proj⁡Wv⃗P:V\to V,\quad P(\vec{v})=\operatorname{proj}_W\vec{v}

稱 PP 為 VV 在 WW 上之 orthogonal projection operator。

Note

設 P:V→VP:V\to V 為 VV 在 WW 上之 orthogonal projection operator。

  • PP is linear。
  • Im⁡(P)=W\operatorname{Im}(P)=W。
  • 對所有 v⃗∈V\vec{v}\in V,P(v⃗)∈WP(\vec{v})\in W,故 P(P(v⃗))=P(v⃗)P(P(\vec{v}))=P(\vec{v})
    因此 P2=PP^2=P (idempotent)。
定理

設 P:V→VP:V\to V 為 VV 在 WW 上之 orthogonal projection operator。
則對所有 v⃗∈V\vec{v}\in V 與 w⃗∈W\vec{w}\in W: ∥v⃗−P(v⃗)∥≤∥v⃗−w⃗∥\lVert\vec{v}-P(\vec{v})\rVert\le\lVert\vec{v}-\vec{w}\rVert

證明

因為 P(v⃗)−w⃗∈WP(\vec{v})-\vec{w}\in W,且 v⃗−P(v⃗)⊥W\vec{v}-P(\vec{v})\perp W,所以

∥v⃗−w⃗∥2=∥(v⃗−P(v⃗))+(P(v⃗)−w⃗)∥2=∥v⃗−P(v⃗)∥2+∥P(v⃗)−w⃗∥2≥∥v⃗−P(v⃗)∥2\begin{aligned} \lVert\vec{v}-\vec{w}\rVert^2 &=\lVert(\vec{v}-P(\vec{v}))+(P(\vec{v})-\vec{w})\rVert^2\\ &=\lVert\vec{v}-P(\vec{v})\rVert^2+\lVert P(\vec{v})-\vec{w}\rVert^2\\ &\ge\lVert\vec{v}-P(\vec{v})\rVert^2 \end{aligned}

Note

Bessel's inequality: ∥P(v⃗)∥≤∥v⃗∥\lVert P(\vec{v})\rVert\le\lVert\vec{v}\rVert


定理

設 A∈Fm×nA\in F^{m\times n},則

  • ker⁡(A)=ker⁡(AHA)\ker(A)=\ker(A^HA)。
  • rank⁡(A)=rank⁡(AHA)\operatorname{rank}(A)=\operatorname{rank}(A^HA)。
Note

rank⁡(A)=rank⁡(AT)=rank⁡(AAT)=rank⁡(ATA)\operatorname{rank}(A)=\operatorname{rank}(A^T)=\operatorname{rank}(AA^T)=\operatorname{rank}(A^TA)

但一般而言,rank⁡(AB)≠rank⁡(BA)\operatorname{rank}(AB)\ne\operatorname{rank}(BA)。

定理

設 A∈Fm×nA\in F^{m\times n},則 A 的 columns 為 LI⟺AHA 為可逆A\text{ 的 columns 為 LI}\Longleftrightarrow A^HA\text{ 為可逆}

證明

A 的 columns 為 LI⟺ker⁡(A)={0⃗}⟺ker⁡(AHA)={0⃗}⟺AHA 為可逆A\text{ 的 columns 為 LI}\Longleftrightarrow\ker(A)=\{\vec{0}\}\Longleftrightarrow\ker(A^HA)=\{\vec{0}\}\Longleftrightarrow A^HA\text{ 為可逆}

Note

設 A∈Fm×nA\in F^{m\times n}:

AA: 行獨立 AA: 列獨立
AA 列生成 AA 行生成
AA 具左反 AA 具右反
rank⁡(A)=n\operatorname{rank}(A)=n rank⁡(A)=m\operatorname{rank}(A)=m
AA nonsingular AA nonsingular
ker⁡(A)={0⃗}\ker(A)=\{\vec{0}\} Lker⁡(A)={0⃗}\operatorname{Lker}(A)=\{\vec{0}\}
AHAA^HA 為可逆 AAHAA^H 為可逆

定理

矩陣版正交投影公式:
設 A∈Fm×nA\in F^{m\times n} 的 columns 為 LI,W=R(A)W=R(A),且 b⃗∈Fm×1\vec{b}\in F^{m\times1}。

則 proj⁡Wb⃗=A(AHA)−1AHb⃗\operatorname{proj}_W\vec{b}=A(A^HA)^{-1}A^H\vec{b}

定理

設 A∈Fm×nA\in F^{m\times n},W=R(A)W=R(A),且 b⃗∈Fm×1\vec{b}\in F^{m\times1}。

則 ∥Ax⃗−b⃗∥ 最小⟺Ax⃗=proj⁡Wb⃗⟺AHAx⃗=AHb⃗\lVert A\vec{x}-\vec{b}\rVert\text{ 最小}\Longleftrightarrow A\vec{x}=\operatorname{proj}_W\vec{b}\Longleftrightarrow A^HA\vec{x}=A^H\vec{b}

Note
  • 如何求 proj⁡R(A)b⃗\operatorname{proj}_{R(A)}\vec{b}:

    • Step 1: solve AHAx⃗=AHb⃗A^HA\vec{x}=A^H\vec{b},得 x⃗0\vec{x}_0
    • Step 2: proj⁡R(A)b⃗=Ax⃗0\operatorname{proj}_{R(A)}\vec{b}=A\vec{x}_0
  • 稱 AHAx⃗=AHb⃗A^HA\vec{x}=A^H\vec{b} 為 normal equation。

  • normal equation 一定有解。

  • 當 AA 的 columns 為 LI 時,AHAA^HA 可逆
    故 x⃗0=(AHA)−1AHb⃗\vec{x}_0=(A^HA)^{-1}A^H\vec{b}
    proj⁡R(A)b⃗=Ax⃗0=A(AHA)−1AHb⃗\operatorname{proj}_{R(A)}\vec{b}=A\vec{x}_0=A(A^HA)^{-1}A^H\vec{b}

  • 當 AA 的 columns 不為 LI 時,AHAx⃗=AHb⃗A^HA\vec{x}=A^H\vec{b} 有無限多解,但 Ax⃗0A\vec{x}_0 唯一。

  • 若 x⃗\vec{x} 使得 ∥Ax⃗−b⃗∥\lVert A\vec{x}-\vec{b}\rVert 最小,稱 x⃗\vec{x} 為 Ax⃗=b⃗A\vec{x}=\vec{b} 的 least squares solution。


Note

設 P∈Fm×mP\in F^{m\times m} 滿足 Px⃗=proj⁡Wx⃗,∀x⃗P\vec{x}=\operatorname{proj}_W\vec{x}, \forall\vec{x}

稱 PP 為投影在 WW 上之 orthogonal projection matrix。

Note

若 A∈Fm×nA\in F^{m\times n} 的 columns 為 LI
則 proj⁡R(A)b⃗=A(ATA)−1ATb⃗\operatorname{proj}_{R(A)}\vec{b}=A(A^TA)^{-1}A^T\vec{b}

因此 P=A(ATA)−1ATP=A(A^TA)^{-1}A^T 為 R(A)R(A) 上之 orthogonal projection matrix。

定理
  • 若 A∈Fm×nA\in F^{m\times n} 的 columns 為 LI,且 P=A(ATA)−1ATP=A(A^TA)^{-1}A^T
    則 P2=PP^2=P 且 PT=PP^T=P。
  • 若 P∈Fm×mP\in F^{m\times m},P2=PP^2=P 且 PT=PP^T=P,則 PP 為 W=R(P)W=R(P) 上之 orthogonal projection matrix。
Note

設 A∈Fm×nA\in F^{m\times n} 的 columns 為 LI
且 P=A(ATA)−1ATP=A(A^TA)^{-1}A^T,則:

  • R(P)=R(A)R(P)=R(A)。
  • rank⁡(P)=rank⁡(A)=n\operatorname{rank}(P)=\operatorname{rank}(A)=n。

最小平方直線

迴歸線

給定 nn 個 points (x1,y1),…,(xn,yn)(x_1,y_1),\ldots,(x_n,y_n),其中 x1,…,xnx_1,\ldots,x_n 不全相同。求一條直線 y=a+bx y=a+bx

使得各點到直線的高度平方和 ∑i=1n(a+bxi−yi)2\sum_{i=1}^n(a+bx_i-y_i)^2 為最小。

定理

設 A=QR∈Fm×nA=QR\in F^{m\times n} 的 columns 為 LI,其中 QQ 的 columns 為 orthonormal,RR 為 upper triangular 且可逆。則 ∥Ax⃗−b⃗∥ 最小⟺Rx⃗=QTb⃗\lVert A\vec{x}-\vec{b}\rVert\text{ 最小} \Longleftrightarrow R\vec{x}=Q^T\vec{b}

故 x⃗=R−1QTb⃗\vec{x}=R^{-1}Q^T\vec{b}

  • QTQ=InQ^TQ=I_n: columns 為 orthonormal,自己與自己內積才為 1,否則為 0
Note

若 W=R(A)W=R(A),且 AA 的 columns 為 orthonormal,則 proj⁡Wb⃗=A(ATA)−1ATb⃗=AATb⃗\operatorname{proj}_W\vec{b}=A(A^TA)^{-1}A^T\vec{b}=AA^T\vec{b}

Note

∥b⃗−proj⁡Wb⃗∥\lVert\vec{b}-\operatorname{proj}_W\vec{b}\rVert 稱為 least square error。

正交補空間

定義

設 S⊆VS\subseteq V,定義
S⊥={v⃗∈V∣⟨v⃗,s⃗⟩=0, ∀s⃗∈S}S^\perp=\{\vec{v}\in V\mid \langle\vec{v},\vec{s}\rangle=0,\ \forall\vec{s}\in S\}
稱為 SS 之 orthogonal complement。

若 V=R2V=\mathbb{R}^2,則

  • S={(1,0)}⟹S⊥=y-axisS=\{(1,0)\} \Longrightarrow S^\perp=\text{y-axis}
  • S={(0,1)}⟹S⊥=x-axisS=\{(0,1)\} \Longrightarrow S^\perp=\text{x-axis}
  • S=x-axis⟹S⊥=y-axisS=\text{x-axis} \Longrightarrow S^\perp=\text{y-axis}
Note
  • S⊥⊆VS^\perp\subseteq V
  • S⊆S⊥⊥S\subseteq S^{\perp\perp},其中 S⊥⊥=(S⊥)⊥S^{\perp\perp}=(S^\perp)^\perp
  • S∩S⊥⊆{0⃗}S\cap S^\perp\subseteq\{\vec{0}\}

Note
  • 若 W={(x,y,z)∈R3∣2x−y+3z=0}W=\{(x,y,z)\in\mathbb{R}^3\mid 2x-y+3z=0\}
    則 WW 為 hyperplane,且 W⊥=span⁡{(2,−1,3)}W^\perp=\operatorname{span}\{(2,-1,3)\}
  • 若 W={(x,y,z)∈R3|x−y+z=02x−y+3z=0}W=\left\{(x,y,z)\in\mathbb{R}^3\middle|\begin{aligned}x-y+z&=0\\2x-y+3z&=0\end{aligned}\right\}
    則 W⊥=span⁡{(1,−1,1),(2,−1,3)}W^\perp=\operatorname{span}\{(1,-1,1),(2,-1,3)\}
  • 若 W={(x,y,z)∈R3∣x−y+2z=3}W=\{(x,y,z)\in\mathbb{R}^3\mid x-y+2z=3\},則其法向量為 (1,−1,2)(1,-1,2)。此 WW 不是 subspace,不可寫成 span
定理

設 W⊆VW\subseteq V,P:V→VP:V\to V 為 orthogonal projection on WW
則 Im⁡(P)=W,ker⁡(P)=W⊥\operatorname{Im}(P)=W, \ker(P)=W^\perp

Note
  • P2=PP^2=P,且 V=Im⁡(P)⊕ker⁡(P)=W⊕W⊥V=\operatorname{Im}(P)\oplus\ker(P)=W\oplus W^\perp
  • dim⁡(V)=dim⁡(W)+dim⁡(W⊥)\dim(V)=\dim(W)+\dim(W^\perp)
  • dim⁡(V)=dim⁡(W⊥)+dim⁡(W⊥⊥)\dim(V)=\dim(W^\perp)+\dim(W^{\perp\perp}),故 W=W⊥⊥W=W^{\perp\perp}
  • 對所有 v⃗∈V\vec{v}\in V,v⃗=proj⁡Wv⃗+(v⃗−proj⁡Wv⃗)\vec{v}=\operatorname{proj}_W\vec{v}+(\vec{v}-\operatorname{proj}_W\vec{v})
    其中 proj⁡Wv⃗∈W\operatorname{proj}_W\vec{v}\in W,v⃗−proj⁡Wv⃗∈W⊥\vec{v}-\operatorname{proj}_W\vec{v}\in W^\perp
  • v⃗=proj⁡Wv⃗+proj⁡W⊥v⃗\vec{v}=\operatorname{proj}_W\vec{v}+\operatorname{proj}_{W^\perp}\vec{v}
  • 若 Q(v⃗)=v⃗−P(v⃗)Q(\vec{v})=\vec{v}-P(\vec{v}),則 Q=I−PQ=I-P

定理

設 A∈Fm×nA\in F^{m\times n},則

  • R(A)⊥=N(AT)R(A)^\perp=N(A^T)
  • R(AT)⊥=N(A)R(A^T)^\perp=N(A)
  • N(A)⊥=R(AT)N(A)^\perp=R(A^T)
  • N(AT)⊥=R(A)N(A^T)^\perp=R(A)
Note
  • 若 AA 的 columns 為 LI,則 proj⁡R(A)v⃗=A(ATA)−1ATv⃗\operatorname{proj}_{R(A)}\vec{v}=A(A^TA)^{-1}A^T\vec{v}
  • 若 AA 的 rows 為 LI,則 proj⁡R(AT)v⃗=AT(AAT)−1Av⃗\operatorname{proj}_{R(A^T)}\vec{v}=A^T(AA^T)^{-1}A\vec{v}
  1. 若 AA 的 rows 為 LI,則 proj⁡N(A)v⃗=v⃗−AT(AAT)−1Av⃗\operatorname{proj}_{N(A)}\vec{v}=\vec{v}-A^T(AA^T)^{-1}A\vec{v}
  2. 若 AA 的 columns 為 LI,則 proj⁡N(AT)v⃗=v⃗−A(ATA)−1ATv⃗\operatorname{proj}_{N(A^T)}\vec{v}=\vec{v}-A(A^TA)^{-1}A^T\vec{v}

Note

對 Ax⃗=b⃗A\vec{x}=\vec{b}:

  • 若唯一解,解之
  • 若無解,求 least squares solution 當近似解
  • 若無限多解,求其中使 ∥x⃗∥2\lVert\vec{x}\rVert_2 最小者,稱為 minimal solution
定理

若 Ax⃗=b⃗A\vec{x}=\vec{b} 有解,則:

  • 存在唯一 s⃗∈R(AT)\vec{s}\in R(A^T),為 Ax⃗=b⃗A\vec{x}=\vec{b} 之 minimal solution
  • 若 u⃗\vec{u} 為 AATu⃗=b⃗AA^T\vec{u}=\vec{b} 之解,則 s⃗=ATu⃗\vec{s}=A^T\vec{u}