正交基令內積和坐標計算更容易,但題目給出的基底通常只是一組任意的線性無關向量。
Gram-Schmidt 按次序逐步扣除已由較早方向解釋的成分,將它變成正交向量組。輸入次序
很重要:次序改變可能令輸出不同,但輸出仍然張成同一個子空間。
投影想法
先考慮一個非零向量 v 和一個向量 w。要從 w 中扣除 v 的一個倍數,並讓餘量垂直於 v,係數 c 必須滿足
0=⟨w−cv,v⟩=⟨w,v⟩−c∥v∥2,c=∥v∥2⟨w,v⟩.
由於 v=0,分母為正,係數亦唯一。對多個兩兩正交的方向,扣除其中一個方向的分量不會改變與其餘方向的內積。
假設 v1,…,vk 已經是非零而且兩兩正交的向量,並令 w∈Rm。w
沿 vi 的分量係數是
∥vi∥2⟨w,vi⟩.
扣除所有舊有正交方向後留下的部分是
r=w−i=1∑k∥vi∥2⟨w,vi⟩vi.
分母合法,是因為正交集不含零向量。當 k=0 時,這個和是空和,r=w。這對零子空間
很有用,但不會要求我們把零向量正規化。
定理
正交餘量定理
設 S={v1,…,vk} 是 Rm 中的正交集,而 w∈Rm。定義
r=w−j=1∑k∥vj∥2⟨w,vj⟩vj.則對每個 i 都有 ⟨r,vi⟩=0,即 r 與 S 中每個向量正交。
固定一個 vi 展開內積:
⟨r,vi⟩=⟨w,vi⟩−j=1∑k∥vj∥2⟨w,vj⟩⟨vj,vi⟩.
當 j=i 時的各項均為零,剩下
⟨w,vi⟩−⟨w,vi⟩=0。若 r=0,只表示 w 已在舊有張成空間
中;對任意 w 而言這種情況可能發生。Gram-Schmidt 要求輸入線性無關,正是為了避免
每一步得到零餘量。
投影與張成空間判別
對正交列表 S={v1,…,vk},同一公式亦給出成員資格判別。定義
r(w)=w−i=1∑k∥vi∥2⟨w,vi⟩vi.
正交餘量定理給出 r(w)⊥spanS,而且
r(w)=0⟺w∈spanS.
正向由投影公式立即得到。反過來,若 w=∑iαivi,與 vi 取內積,
因交叉項為零而得到 αi=⟨w,vi⟩/∥vi∥2;代回便有 r(w)=0。
所以非零餘量正好記錄 w 在目前張成空間以外的部分。
證明
投影與正交餘量
令 U=spanS,並設
p(w)=i=1∑k∥vi∥2⟨w,vi⟩vi,r=w−p(w).則 p(w)∈U;正交餘量定理給出 r 與每個生成向量 vi 正交。由線性性,
r⊥U,所以 w=p+r。若 S 為空,和為 0、U={0},同一結論就是
w=0+w。
若還有 w=p′+r′,其中 p′∈U 且 r′⊥U,令
z=p−p′=r′−r。則 z∈U 且 z⊥U,因此
⟨z,z⟩=0,從而 z=0;所以 p=p′ 且 r=r′。分解為 U 中部分
與正交餘量是唯一的。
這個分解精確說明「刪去舊方向」的含義。p(w) 是由已有正交生成向量重構的部分,而 r
是它們的任何線性組合都不能解釋的部分。若 w∈U,成員資格判別給出 r=0;若
w∈/U,則 r=0,它提供一個垂直於 U 內所有方向的新方向。因此在
Gram-Schmidt 歸納中,獨立輸入不能消失。唯一性亦說明符合這些條件的正交部分不會有兩個
不同答案:z∈U∩U⊥ 時,正定內積迫使 ⟨z,z⟩=0,於是 z=0。
Gram-Schmidt 過程
定理
Gram-Schmidt 正交化過程
設 k≥1,而 w1,…,wk 是 Rm 中一組線性無關
向量。定義
v1=w1,並對 ℓ=2,…,k 定義
vℓ=wℓ−j=1∑ℓ−1∥vj∥2⟨wℓ,vj⟩vj.則每個 vℓ 都非零,v1,…,vk 兩兩正交,而且對每個 ℓ,
span{w1,…,wℓ}=span{v1,…,vℓ}.
空輸入列表另行處理:它的張成空間是零子空間 {0},輸出是空的正交列表。
這情況沒有要定義的 v1,亦沒有向量需要正規化。
每一步必須使用已經構造的正交向量 vj,不能把原來的 wj 當作互相正交。否則
分別扣除各個 wj 的投影,並不能扣掉正確的成分。
證明
歸納中的三個不變量
我們同時證明第 ℓ 步的三件事:(i) 新向量非零;(ii) 已構造向量兩兩正交;
(iii) 兩個部分張成空間相等。基例是 v1=w1=0,兩個張成空間都是同一條直線。
假設這三件事已對 ℓ−1 成立。先前的 vj 非零,所以 vℓ 定義式中的分母均
合法。正交餘量定理給出 vℓ⊥vj(j<ℓ),新的正交性因而成立。
定義式說明 vℓ 是 wℓ 與先前 vj 的線性組合,因此新的 v-張成空間包含
於新的 w-張成空間。整理定義式:
wℓ=vℓ+j=1∑ℓ−1∥vj∥2⟨wℓ,vj⟩vj,便得到反向包含,兩個部分張成空間相等。最後,若 vℓ=0,上式會把 wℓ 放進
span{v1,…,vℓ−1}=span{w1,…,wℓ−1},
與輸入線性無關矛盾。因此新向量非零,三個部分在不預先假設結論的情況下同時完成。
因此正交性不只是若干個點積為零;正交集還必須不含零向量。這個非零檢查令後面的正規化
有合法分母。
在三維空間中扣除已有方向
例題
R3 中的三個線性無關向量
取有序向量組
w1=101,w2=111,w3=012.以它們為欄的矩陣行列式是 2,所以輸入線性無關。先有
v1=w1=101,∥v1∥2=2.第二步的係數為
∥v1∥2⟨w2,v1⟩=22=1,v2=w2−v1=010.第三步必須同時處理兩個已構造的方向:
∥v1∥2⟨w3,v1⟩=22=1,∥v2∥2⟨w3,v2⟩=11=1,所以
v3=w3−v1−v2=−101.三個平方範數是 2,1,2,而且
⟨v1,v2⟩=⟨v1,v3⟩=⟨v2,v3⟩=0.張成空間保持也可直接看出:
w2=v1+v2,w3=v1+v2+v3.因此輸出與輸入張成完全相同的 R3。
第三步的算術可在接受答案前逐項核對:⟨w3,v1⟩=0⋅1+1⋅0+2⋅1=2,
而 ⟨w3,v2⟩=0⋅0+1⋅1+2⋅0=1。所以扣除
(1,0,1)T+(0,1,0)T 後確實得到 (−1,0,1)T,負的第一個座標不是猜出的符號。
反例模式
為甚麼不能用原來的 w 向量代替 v 向量
若錯誤地用原向量計算第三步,會得到
w3−∥w1∥2⟨w3,w1⟩w1−∥w2∥2⟨w3,w2⟩w2=w3−w1−w2=−200.這個向量與 w1 的點積是 −2,所以並不正交。修正方法是對 v1,v2 扣除投影,
所得 v3=(−1,0,1)T 才正確。兩組向量不同,但每一步的部分張成空間相同。
準確保留投影係數中的分數
這裏的投影係數是分數。每次相減時保留精確分數,再核對正交性與張成空間相等。
例題
三維空間中的分數投影係數
從
w1=(1,1,0)T,w2=(1,0,1)T,w3=(0,1,1)T開始。
以這些向量為欄的行列式是 −2,所以它們線性無關。這裏 ∥v2∥2=3/2,而
⟨w3,v2⟩=0(1/2)+1(−1/2)+1(1)=1/2;第三步的第二個系數是
(1/2)/(3/2)=1/3。所得 v3 的平方範數為 4/3,直接計算亦給出
v1⋅v2=v1⋅v3=v2⋅v3=0。
正規化與邊界情況
定理
由正交基得到標準正交基
若 v1,…,vk 是正交基,定義
ui=∥vi∥vi.每個 ui 的範數都是 1,縮放不改變兩兩正交性,所以得到同一子空間的標準正交基。
例題
正規化第一個 R3 輸出
對第一個 R3 例子,
∥v1∥=2,∥v2∥=1,∥v3∥=2.因此
u1=21101,u2=010,u3=21−101.這裏沿用第一個輸出,而不是靜態表格的另一組輸入。
因為範數為正,ui 與 vi 方向相同;乘回 ∥vi∥ 便得到原向量,而不是它的相反向量。
單位長度可直接核對:∥u1∥2=(1/2)∥v1∥2=1、∥u2∥2=1、
∥u3∥2=(1/2)∥v3∥2=1。而且 vi=∥vi∥ui,把 ui 乘回範數便得到
vi;因此正規化保持每一個張成空間等式的兩個方向。
例題
由 R4 約束得到標準正交基
令
V=N[1111]={x∈R4:x1+x2+x3+x4=0}.寫成 x=(−a−b−c,a,b,c)T,可見 V 的一組基為
w1=−1100,w2=−1010,w3=−1001.這些向量線性無關:若 a1w1+a2w2+a3w3=0,第 2、3、4 個座標分別迫使
a1=a2=a3=0。約束方程亦直接說明它們都在 V 內。
Gram-Schmidt 給出
v1=w1,v2=w2−21v1=−1/2−1/210,這裏 ⟨w2,v1⟩=1、∥v1∥2=2,所以係數為 1/2。另外
⟨w3,v1⟩=1,而 ⟨w3,v2⟩=1/2、∥v2∥2=3/2,
所以係數為 (1/2)/(3/2)=1/3。於是
v3=w3−21v1−31v2=−1/3−1/3−1/31.平方範數為 2,3/2,4/3,正規化後得到
u1=21−1100,u2=61−1−120,u3=121−1−1−13.例如 ∥v2∥2=1/4+1/4+1=3/2,∥v3∥2=1/9+1/9+1/9+1=4/3;
顯示的縮放因子確實令三個範數都是 1。每個向量的坐標和為零,每一對的點積為零。又因
vi=∥vi∥ui,把 ui 乘回範數即可得到 vi,所以正規化保持每一個張成空間等式
的兩個方向。由於零空間維數是 3,這三個向量構成 V 的標準正交基。
R4 計算亦可雙向核對張成空間。由三個定義式,
w2=v2+21v1,w3=v3+21v1+31v2.反過來,每個 vi 都由相應的 wi 與較早向量構成;這是對角系數全為 1 的三角變換,
所以沒有丟失任何方向。正規化後的向量仍在 V 內:u1 的座標和為
(−1+1)/2=0,u2 的座標和為 (−1−1+2)/6=0,u3 的座標和為
(−1−1−1+3)/12=0。它們的平方範數分別是 2/2、(1+1+4)/6 及
(1+1+1+9)/12,全部等於 1;這些座標核對獨立確認了約束與正規化。
定理
每個有限維子空間都有標準正交基
Rm 的每個子空間都有有限基。正維數時,Gram-Schmidt 產生非零正交基,再正規化
便得到標準正交基。零子空間是邊界情況:它的空集基已是正交基和標準正交基,不需要對
任何向量除以範數。
常見錯誤
相關輸入會產生零餘量
若 w1=(1,0)T、w2=(2,0)T,則
v2=w2−∥v1∥2⟨w2,v1⟩v1=(2,0)T−2(1,0)T=0.下一步正規化便沒有定義。應先按線性無關性檢查並刪去冗餘向量;Gram-Schmidt 的輸入應是
線性無關的有序向量組。
快速檢查
思考檢查
Gram-Schmidt 的第一個向量 v1 是甚麼?
解答 · 答案
v1=w1;第一個輸入方向原樣保留。
思考檢查
Gram-Schmidt 除了正交性還保持甚麼性質?
解答 · 答案
它保持每個部分張成空間:
span{w1,…,wℓ}=span{v1,…,vℓ}。
解答 · 答案
對每個非零 vi 除以範數:ui=vi/∥vi∥。
思考檢查
對 w1=(1,2,0) 與 w2=(2,0,1) 做前兩步 Gram-Schmidt。
解答 · 引導解答
v1=(1,2,0)T,投影係數為
⟨w2,v1⟩/∥v1∥2=2/5,因此
v2=(2,0,1)T−52(1,2,0)T=(8/5,−4/5,1)T.核對 v1⋅v2=8/5−8/5=0。而且 w2=v2+52v1,
所以新舊兩組向量張成同一個平面。
思考檢查
為甚麼獨立輸入列表出現零餘量就表示有問題?
解答 · 引導解答
餘量為零當且僅當目前向量在先前向量的張成空間中,這與線性無關性矛盾。
思考檢查
若正交向量 v1,v2 的範數分別是 3 和 4,對應的單位向量是甚麼?
解答 · 引導解答
它們是 u1=v1/3 與 u2=v2/4;兩者範數為 1,仍然正交。
相關筆記
先閱讀
9.2 正交集與標準正交基,
因為 Gram-Schmidt 建立在正交係數公式上。約束空間例子也使用
6.2 子空間與
6.5 基底與維數的零空間與基方法。
下一節
9.4 Cauchy-Schwarz 與三角不等式
會發展控制這些內積計算的估計。