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

Ch7 內積空間

內積

定義

VV 為 vector space over FF
,:V×VF\langle \cdot,\cdot\rangle:V\times V\to F 為 function,且對所有 u,v,wV\vec{u},\vec{v},\vec{w}\in Vc,dFc,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=cu,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}
  • v0\vec{v}\ne \vec{0},則 v,v>0\langle \vec{v},\vec{v}\rangle>0

,\langle \cdot,\cdot\rangleVV 上之 inner product

Note
  • 左線性: cu+dv,w=cu,w+dv,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=yHxin Fn×1\langle \vec{x},\vec{y}\rangle=\vec{y}^{H}\vec{x}\quad \text{in }F^{n\times 1}
    • matrix form: x,y=xyHin 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=1nj=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_nx1,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(x1)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\rangleVV 上的 inner product。

  • v,0=0=0,v\langle \vec{v},\vec{0}\rangle=0=\langle \vec{0},\vec{v}\rangle
  • cu,v=cu,v\langle c\vec{u},\vec{v}\rangle=c\langle \vec{u},\vec{v}\rangle
  • iciui,v=iciui,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,icivi=icˉiu,vi\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,v0\langle \vec{v},\vec{v}\rangle\ge 0,且 v,v=0v=0\langle \vec{v},\vec{v}\rangle=0\Longleftrightarrow\vec{v}=\vec{0}
  • u,v=0\langle \vec{u},\vec{v}\rangle=0u,vVu=0\vec{u},\vec{v}\in V\Longleftrightarrow\vec{u}=\vec{0}v=0\vec{v}=\vec{0}
  • β={v1,,vn}\beta=\{\vec{v}_1,\ldots,\vec{v}_n\}VV 的 basis,且 u,vi=0\langle \vec{u},\vec{v}_i\rangle=0i=1,,ni=1,\ldots,n,則 u=0\vec{u}=\vec{0}
定義

,\langle \cdot,\cdot\rangleVV 上的 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,記作 uv\vec{u}\perp\vec{v}
  • S={v1,,vk}S=\{\vec{v}_1,\ldots,\vec{v}_k\},若 vivj\vec{v}_i\perp\vec{v}_jij\forall i\ne j,稱 SSorthogonal set,允許 0\vec{0}
  • SSorthogonal setvi=1\lVert\vec{v}_i\rVert=1i=1,,ki=1,\ldots,k,稱 SSorthonormal set,不允許 0\vec{0}
  • u\vec{u}v\vec{v} 的 distance 為 d(u,v)=uvd(\vec{u},\vec{v})=\lVert\vec{u}-\vec{v}\rVert
Note

S={v1,,vk}S=\{\vec{v}_1,\ldots,\vec{v}_k\},則

S 為 orthonormal setvi,vj=δij={1,i=j,0,ijS\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,cosx,cos2x,}S=\left\{\frac{1}{\sqrt{2}},\cos x,\cos 2x,\ldots\right\}orthonormal set

Note

FnF^n 的 standard inner product 下:

x=x,x=(x12++xn2)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。

也可定義:

  • x1=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|\}
  • xp=(x1p++xnp)1/p\lVert\vec{x}\rVert_p=\left(|x_1|^p+\cdots+|x_n|^p\right)^{1/p}

norm 需滿足:

  • x0\lVert\vec{x}\rVert\ge 0,且 x=0x=0\lVert\vec{x}\rVert=0\Longleftrightarrow\vec{x}=\vec{0}
  • cx=cx\lVert c\vec{x}\rVert=|c|\lVert\vec{x}\rVert
  • x+yx+y\lVert\vec{x}+\vec{y}\rVert\le\lVert\vec{x}\rVert+\lVert\vec{y}\rVert

定理

Pythagorean theorem:
,\langle\cdot,\cdot\rangleVV 上的 inner product。若 uv\vec{u}\perp\vec{v},則 u+v2=u2+v2. \lVert\vec{u}+\vec{v}\rVert^2=\lVert\vec{u}\rVert^2+\lVert\vec{v}\rVert^2.

Note

{v1,,vk}\{\vec{v}_1,\ldots,\vec{v}_k\}orthogonal set,則

c1v1++ckvk2=c1v12++ckvk2=c12v12++ck2vk2\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\rangleVV 上的 inner product,則對所有 u,vV\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+vu+v.\lVert\vec{u}+\vec{v}\rVert\le\lVert\vec{u}\rVert+\lVert\vec{v}\rVert.

證明
  • v=0\vec{v}=\vec{0},結論顯然。若 v0\vec{v}\ne\vec{0},取 α=u,vv,v.\alpha=\frac{\langle\vec{u},\vec{v}\rangle}{\langle\vec{v},\vec{v}\rangle}.
    0uαv20\le\lVert\vec{u}-\alpha\vec{v}\rVert^2,可得 u,v2u2v2.|\langle\vec{u},\vec{v}\rangle|^2\le\lVert\vec{u}\rVert^2\lVert\vec{v}\rVert^2.

  • 展開後套用 Cauchy-Schwarz inequality:
    u+v2=u2+2Reu,v+v2u2+2uv+v2=(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)2dx)1/2(abg(x)2dx)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,v0\vec{u},\vec{v}\ne\vec{0},由 Cauchy-Schwarz inequality 可得

1u,vuv1-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,vuv\cos\theta=\frac{\langle\vec{u},\vec{v}\rangle}{\lVert\vec{u}\rVert\lVert\vec{v}\rVert}

θ\thetau\vec{u}v\vec{v} 之間的角度。


定理

S={v1,,vk}S=\{\vec{v}_1,\ldots,\vec{v}_k\}nonzero orthogonal set,且 vspan(S)\vec{v}\in\operatorname{span}(S)。若 v=c1v1++ckvk\vec{v}=c_1\vec{v}_1+\cdots+c_k\vec{v}_k

ci=v,vivi,vi,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_iFourier coefficient

SSorthonormal set,則 ci=v,vic_i=\langle\vec{v},\vec{v}_i\rangle

定理

S={v1,,vk}S=\{\vec{v}_1,\ldots,\vec{v}_k\}nonzero orthogonal set,則 SS 為 LI。

Gram-Schmidt 正交化

定理

Gram-Schmidt orthogonalization Process:
給定 {v1,,vn}\{\vec{v}_1,\ldots,\vec{v}_n\} 為 LI。令

u1=v1,ui=vij=1i1vi,ujuj,ujuj,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}

{u1,,un}\{\vec{u}_1,\ldots,\vec{u}_n\}nonzero orthogonal set,且 span{u1,,un}=span{v1,,vn}\operatorname{span}\{\vec{u}_1,\ldots,\vec{u}_n\}=\operatorname{span}\{\vec{v}_1,\ldots,\vec{v}_n\}

Note

u0\vec{u}\ne\vec{0},令 z=uu\vec{z}=\frac{\vec{u}}{\lVert\vec{u}\rVert}

稱為 normalized vector,且 z=1\lVert\vec{z}\rVert=1
因此 {u1u1,,unun}\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:

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

其中 QFm×nQ\in F^{m\times n} 的 columns 為 orthonormalRFn×nR\in F^{n\times n} 為 upper triangular 且可逆。

正交投影

定義

WVW\subseteq V,且 vV\vec{v}\in V。若存在 v0W\vec{v}_0\in W,使得 vv0,w=0,wW\langle\vec{v}-\vec{v}_0,\vec{w}\rangle=0,\quad\forall\vec{w}\in W

vv0W\vec{v}-\vec{v}_0\perp W,稱 v0\vec{v}_0v\vec{v} on WWorthogonal projection vector

定理

WVW\subseteq VvV\vec{v}\in V,且 {v1,,vk}\{\vec{v}_1,\ldots,\vec{v}_k\}WWorthogonal basis。令

v0=i=1kv,vivi,vivi\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

v0\vec{v}_0v\vec{v}WW 上唯一的 orthogonal projection vector,記作 projWv\operatorname{proj}_W\vec{v}

W=span{v}W=\operatorname{span}\{\vec{v}\} 時:

projWu=u,vv,vv\operatorname{proj}_W\vec{u} =\frac{\langle\vec{u},\vec{v}\rangle}{\langle\vec{v},\vec{v}\rangle}\vec{v}

Note
  • vWprojWv=v\vec{v}\in W\Longleftrightarrow\operatorname{proj}_W\vec{v}=\vec{v}
  • v=0projWv=0\vec{v}=\vec{0}\Longrightarrow\operatorname{proj}_W\vec{v}=\vec{0}
  • vWprojWv=0\vec{v}\perp W\Longleftrightarrow\operatorname{proj}_W\vec{v}=\vec{0}

定義

WVW\subseteq V。定義 P:VV,P(v)=projWvP:V\to V,\quad P(\vec{v})=\operatorname{proj}_W\vec{v}

PPVVWW 上之 orthogonal projection operator

Note

P:VVP:V\to VVVWW 上之 orthogonal projection operator

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

P:VVP:V\to VVVWW 上之 orthogonal projection operator
則對所有 vV\vec{v}\in VwW\vec{w}\in W: vP(v)vw\lVert\vec{v}-P(\vec{v})\rVert\le\lVert\vec{v}-\vec{w}\rVert

證明

因為 P(v)wWP(\vec{v})-\vec{w}\in W,且 vP(v)W\vec{v}-P(\vec{v})\perp W,所以

vw2=(vP(v))+(P(v)w)2=vP(v)2+P(v)w2vP(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


定理

AFm×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)

定理

AFm×nA\in F^{m\times n},則 A 的 columns 為 LIAHA 為可逆A\text{ 的 columns 為 LI}\Longleftrightarrow A^HA\text{ 為可逆}

證明

A 的 columns 為 LIker(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

AFm×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 為可逆

定理

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

projWb=A(AHA)1AHb\operatorname{proj}_W\vec{b}=A(A^HA)^{-1}A^H\vec{b}

定理

AFm×nA\in F^{m\times n}W=R(A)W=R(A),且 bFm×1\vec{b}\in F^{m\times1}

Axb 最小Ax=projWbAHAx=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
  • 如何求 projR(A)b\operatorname{proj}_{R(A)}\vec{b}:

    • Step 1: solve AHAx=AHbA^HA\vec{x}=A^H\vec{b},得 x0\vec{x}_0
    • Step 2: projR(A)b=Ax0\operatorname{proj}_{R(A)}\vec{b}=A\vec{x}_0
  • AHAx=AHbA^HA\vec{x}=A^H\vec{b}normal equation

  • normal equation 一定有解。

  • AA 的 columns 為 LI 時,AHAA^HA 可逆,故
    $\vec{x}_0=(A^HA)^{-1}A^H\vec{b} $
    projR(A)b=Ax0=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=AHbA^HA\vec{x}=A^H\vec{b} 有無限多解,但 Ax0A\vec{x}_0 唯一。

  • x\vec{x} 使得 Axb\lVert A\vec{x}-\vec{b}\rVert 最小,稱 x\vec{x}Ax=bA\vec{x}=\vec{b}least squares solution


Note

PFm×mP\in F^{m\times m} 滿足 Px=projWx,xP\vec{x}=\operatorname{proj}_W\vec{x}, \forall\vec{x}

PP 為投影在 WW 上之 orthogonal projection matrix

Note

AFm×nA\in F^{m\times n} 的 columns 為 LI
projR(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^TR(A)R(A) 上之 orthogonal projection matrix

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

AFm×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+bxiyi)2\sum_{i=1}^n(a+bx_i-y_i)^2 為最小。

定理

A=QRFm×nA=QR\in F^{m\times n} 的 columns 為 LI,其中 QQ 的 columns 為 orthonormalRR 為 upper triangular 且可逆。則 Axb 最小Rx=QTb\lVert A\vec{x}-\vec{b}\rVert\text{ 最小} \Longleftrightarrow R\vec{x}=Q^T\vec{b}

x=R1QTb\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,則 projWb=A(ATA)1ATb=AATb\operatorname{proj}_W\vec{b}=A(A^TA)^{-1}A^T\vec{b}=AA^T\vec{b}

Note

bprojWb\lVert\vec{b}-\operatorname{proj}_W\vec{b}\rVert 稱為 least square error

正交補空間

定義

SVS\subseteq V,定義
S={vVv,s=0, sS}S^\perp=\{\vec{v}\in V\mid \langle\vec{v},\vec{s}\rangle=0,\ \forall\vec{s}\in S\}
稱為 SSorthogonal 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-axisS=y-axisS=\text{x-axis} \Longrightarrow S^\perp=\text{y-axis}
Note
  • SVS^\perp\subseteq V
  • SSS\subseteq S^{\perp\perp},其中 S=(S)S^{\perp\perp}=(S^\perp)^\perp
  • SS{0}S\cap S^\perp\subseteq\{\vec{0}\}

Note
  • W={(x,y,z)R32xy+3z=0}W=\{(x,y,z)\in\mathbb{R}^3\mid 2x-y+3z=0\}
    WWhyperplane,且 W=span{(2,1,3)}W^\perp=\operatorname{span}\{(2,-1,3)\}
  • W={(x,y,z)R3|xy+z=02xy+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)R3xy+2z=3}W=\{(x,y,z)\in\mathbb{R}^3\mid x-y+2z=3\},則其法向量為 (1,1,2)(1,-1,2)。此 WW 不是 subspace,不可寫成 span
定理

WVW\subseteq VP:VVP:V\to Vorthogonal 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)=WWV=\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=WW=W^{\perp\perp}
  • 對所有 vV\vec{v}\in Vv=projWv+(vprojWv)\vec{v}=\operatorname{proj}_W\vec{v}+(\vec{v}-\operatorname{proj}_W\vec{v})
    其中 projWvW\operatorname{proj}_W\vec{v}\in WvprojWvW\vec{v}-\operatorname{proj}_W\vec{v}\in W^\perp
  • v=projWv+projWv\vec{v}=\operatorname{proj}_W\vec{v}+\operatorname{proj}_{W^\perp}\vec{v}
  • Q(v)=vP(v)Q(\vec{v})=\vec{v}-P(\vec{v}),則 Q=IPQ=I-P

定理

AFm×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,則 projR(A)v=A(ATA)1ATv\operatorname{proj}_{R(A)}\vec{v}=A(A^TA)^{-1}A^T\vec{v}
  • AA 的 rows 為 LI,則 projR(AT)v=AT(AAT)1Av\operatorname{proj}_{R(A^T)}\vec{v}=A^T(AA^T)^{-1}A\vec{v}
  1. AA 的 rows 為 LI,則 projN(A)v=vAT(AAT)1Av\operatorname{proj}_{N(A)}\vec{v}=\vec{v}-A^T(AA^T)^{-1}A\vec{v}
  2. AA 的 columns 為 LI,則 projN(AT)v=vA(ATA)1ATv\operatorname{proj}_{N(A^T)}\vec{v}=\vec{v}-A(A^TA)^{-1}A^T\vec{v}

Note

Ax=bA\vec{x}=\vec{b}:

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

Ax=bA\vec{x}=\vec{b} 有解,則:

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