ℝ² 中点在通过原点的直线上的正交投影和 ℝⁿ 的情形非常类似.对给定向量 y 和 ℝⁿ 中子空间 W ,存在属于 W 的向量 haty 满足:(1) W 中有唯一向量 haty ,使得 y - haty 与 W 正交,(2) haty 是 W 中唯一最接近 y 的向量,见图 6-17. haty 的这两个性质提供了本章介绍性实例中提到的求线性方程组的最小二乘解的方法,完整的内容在6.5节中学习.
为准备第一个定理,我们注意到,当 y 表示成 Rⁿ 空间中向量 u₁, …, uₙ 的线性组合时,y 的和中的各项可分为两部分,使得 y 可写成
mathbf y = mathbf z _ 1 + mathbf z _ 2其中, z₁ 是其中一些 uᵢ 的线性组合, z₂ 是其余 uᵢ 的线性组合。当 {u₁,…,uₙ} 是正交基时,这个思路特别有用。回忆在6.1节, W^(perp) 表示所有与W正交的向量集合。
假若 {u₁,…,u₅} 是 R⁵ 中的正交基,且
boldsymbol y = c _ 1 boldsymbol u _ 1 + … + c _ 5 boldsymbol u _ 5考虑子空间 W = operatornameSpan{u₁, u₂} ,并将 y 写成 W 中向量 z₁ 与 W^(perp) 中向量 z₂ 的和.
解
其中
r l & ±b z _ 1 = c _ 1 ±b u _ 1 + c _ 2 ±b u _ 2 text 属于 Span {±b u _ 1, ±b u _ 2 } & ±b z _ 2 = c _ 3 ±b u _ 3 + c _ 4 ±b u _ 4 + c _ 5 ±b u _ 5 text 属于 Span {±b u _ 3, ±b u _ 4, ±b u _ 5 }为证明 z₂ 属于 W^(perp) ,只需证明 z₂ 与以 {u₁, u₂} 为基的空间 W 中的向量正交(见 6.1 节)。利用内积的性质计算
r l mathbf z _ 2 · mathbf u _ 1 & = (c _ 3 mathbf u _ 3 + c _ 4 mathbf u _ 4 + c _ 5 mathbf u _ 5) · mathbf u _ 1 & = c _ 3 mathbf u _ 3 · mathbf u _ 1 + c _ 4 mathbf u _ 4 · mathbf u _ 1 + c _ 5 mathbf u _ 5 · mathbf u _ 1 & = 0上式中利用了 u₁ 与 u₃ , u₄ 和 u₅ 的正交性。同样的计算可以证明 z₂ · u₂ = 0 ,从而说明 z₂ 属于 W^(perp) 。下面的定理表明,在没有 Rⁿ 中的正交基时,对例1 中的分解 y = z₁ + z₂ 同样可以计算,只需 W 中有一个正交基即可。
若 W 是 Rⁿ 的一个子空间,那么 Rⁿ 中每一个向量 y 可以唯一表示为
mathbf y = hat mathbf y + z (1)其中 haty 属于 W 而 z 属于 W^(perp).实际上,如果 {u₁,…,uₚ} 是 W 的任意正交基,那么
hat mathbf y = frac mathbf y · mathbf u _ 1mathbf u _ 1 · mathbf u _ 1 mathbf u _ 1 + … + frac mathbf y · mathbf u _ pmathbf u _ p · mathbf u _ p mathbf u _ p (2)且 z = y - haty
(1)式中的 haty 称为 y 在 W 上的正交投影,常记作 proj_(w)y 。见图 6-18。当 W 是一维子空间时, haty 的公式和 6.2 节中的公式一致。
证明
若 {u₁,…,uₚ} 是 W 的正交基,且 haty 用(2)式定义。 ^(circledR) 由于 haty 是基 u₁,…,uₚ 的线性组合,所以 haty 属于 W。令 z=y-haty ,由于 u₁ 正交于 u₂,…,uₚ ,因此从(2)式可以得出
mathbf z · mathbf u _ 1 = (mathbf y - hat mathbf y) · mathbf u _ 1 = mathbf y · mathbf u _ 1 - (frac mathbf y · mathbf u _ 1mathbf u _ 1 · mathbf u _ 1) mathbf u _ 1 · mathbf u _ 1 - 0 - … - 0 = mathbf y · mathbf u _ 1 - mathbf y · mathbf u _ 1 = 0从而 z 与 u₁ 正交.类似地,z 与空间 W 的每个基向量 uⱼ 正交,因此 z 与 W 中任何向量正交,即 z 属于 W^(perp) .
为证明分解(1)是唯一的,设 y 也可以写成 y = haty₁ + z₁ ,其中 haty₁ 属于 W 而 z₁ 属于 W^(perp) ,那么 haty + z = haty₁ + z₁ (由于两边的 y 相同),从而得出
hat boldsymbol y - hat boldsymbol y _ 1 = boldsymbol z _ 1 - boldsymbol z这个等式表明向量 v = haty - haty₁ 属于 W,同时又属于 W^(perp) (由于 z₁ 和 z 属于 W^(perp) 且 W^(perp) 是子空间)。因此 v · v = 0 ,这说明 v = 0,即证明 haty = haty₁ 且 z₁ = z 。
分解式(1)的唯一性表明,正交投影 haty 仅依赖于 W,而不依赖于(2)中使用的特殊基.
假设 ±bu₁ = 2 5 -1 , ±bu₂ = -2 1 1 和 ±by = 1 2 3.注意到 {±bu₁, ±bu₂} 是 W = operatornameSpan{±bu₁, ±bu₂} 的正交基,将 ±by 写成属于 W 的向量与正交于 W 的向量之和.
解
y 在 W 上的正交投影是
r l hat mathbf y & = frac mathbf y · mathbf u _ 1mathbf u _ 1 · mathbf u _ 1 mathbf u _ 1 + frac mathbf y · mathbf u _ 2mathbf u _ 2 · mathbf u _ 2 mathbf u _ 2 & = frac 93 0 [ l 2 5 - 1 ] + frac 36 [ l - 2 1 1 ] = frac 93 0 [ l 2 5 - 1 ] + frac 1 53 0 [ l - 2 1 1 ] = [ l - 2 / 5 2 1 / 5 ]且
mathbf y - hat mathbf y = [ c 1 2 3 ] - [ c - 2 / 5 2 1 / 5 ] = [ c 7 / 5 0 1 4 / 5 ]保证 y - haty 属于 W^(perp) 。为检验计算,一个好的方法是验证 y - haty 正交于 u₁ 和 u₂ ,因而正交于 W 。所期望的 y 的分解式是
mathbf y = [ l 1 2 3 ] = [ c - 2 / 5 2 1 / 5 ] + [ c 7 / 5 0 1 4 / 5 ]正交投影的几何解释
当 W 是一维子空间时,公式(2)的 proj_(W)y 仅包含一项。这样,当 dim W > 1 时,(2)式中的每一项都是自身,即向量 y 在由 W 的一个基 u 所生成的一维子空间上的正交投影。图 6-19 表示当 W 是 R³ 中由 u₁ 和 u₂ 所生成的子空间的情形。 haty₁ 和 haty₂ 分别表示 y 在由 u₁ 和 u₂ 所生成的直线上的投影,向量 y 在 W 上的正交投影 haty 在两个相互正交的一维子空间上各自投影之和。图 6 对应于 6.2 节图 6-13 中的向量 y 在两个相互正交的一维子空间上各自投影之和。图 6-19 中的向量 haty 对应于 6.2 节图 6-13 中的向量 y,其原因是现在 haty 属于 W。
正交投影的性质
Section titled “正交投影的性质”如果 {u₁,…,uₚ} 是 W 的正交基,且若 y 正好属于 W,那么 proj_(w)y 的公式和 6.2 节中定理5 里 y 的表达式完全一致。这种情形下, proj_(w)y=y 。
如果 y 属于 W = Span{u₁, …, uₚ} ,那么 proj_(w) y = y
这个结论可以从下列定理得出.
假设 W 是 Rⁿ 的一个子空间,y 是 Rⁿ 中的任意向量, haty 是 y 在 W 上的正交投影,那么 haty 是 W 中最接近 y 的点,也就是
| mathbf y - hat mathbf y | < | mathbf y - mathbf v | (3)对所有属于W又异于 haty 的v成立.
定理9 中的向量 haty 称为 W 中元素对 y 的最佳逼近.在后面章节中,我们会研究这样的问题,即对给定的元素 y,可以被某个给定子空间 W 中的向量 v 代替或 “逼近”. 用 |y-v| 表示的从 y 到 v 的距离可以认为是用 v 代替 的 “误差”. 定理9 说明误差在 的 “误差”. 定理 9 说明误差在 v=haty 处取得最小值.
方程(3)给出了一个新的证明,即 haty 的计算不依赖于特定正交基。如果 W 的一个不同的正交基被用于构造 y 的正交投影,那么这个投影就是 W 中离 y 最近的点,记为 haty 。
证明
取 W 中的 v (不同于 haty ),见图 6-20,则 haty-v 是 W 中的向量。由正交分解定理, y-haty 正交于 W,特别地, y-haty 正交于 haty-v 。由于
boldsymbol y - boldsymbol v = (boldsymbol y - hat boldsymbol y) + (hat boldsymbol y - boldsymbol v)故由勾股定理可知
| mathbf y - mathbf v | ^ 2 = | mathbf y - hat mathbf y | ^ 2 + | hat mathbf y - mathbf v | ^ 2(见图 6-20中右边的直角三角形.每个边的长度均已标出.)因为 haty -v≠ 0 ,所以 | haty -v| ² >0 ,从而不等式(3)立刻推出.
如果 ±bu₁ = 2 5 -1 , ±bu₂ = -2 1 1 , ±by = 1 2 3 和 W = operatornameSpan{±bu₁, ±bu₂} ,类似例2,则 W 中离 ±by 最近的点是
hat mathbf y = frac mathbf y · mathbf u _ 1mathbf u _ 1 · mathbf u _ 1 mathbf u _ 1 + frac mathbf y · mathbf u _ 2mathbf u _ 2 · mathbf u _ 2 mathbf u _ 2 = [ c - 2 / 5 2 1 / 5 ]ℝⁿ 中的一个点 y 到一个子空间 W 的距离定义为从 y 到空间 W 中最近点的距离,求出 y 到 W = operatornameSpan{±bu₁, ±bu₂} 的距离,其中
mathbf y = [ l - 1 - 5 1 0 ], mathbf u _ 1 = [ l 5 - 2 1 ], mathbf u _ 2 = [ l 1 2 - 1 ]解
由最佳逼近定理,从 y 到 W 的距离是 |y-haty| ,其中 haty=operatornameproj_(W)y 。由于 {u₁,u₂} 是 W 的正交基,我们得到
hat mathbf y = frac 1 53 0 mathbf u _ 1 + frac - 2 16 mathbf u _ 2 = frac 12 [ c 5 - 2 1 ] - frac 72 [ c 1 2 - 1 ] = [ c - 1 - 8 4 ]mathbf y - hat mathbf y = [ c - 1 - 5 1 0 ] - [ c - 1 - 8 4 ] = [ c 0 3 6 ]| mathbf y - hat mathbf y | ^ 2 = 3 ^ 2 + 6 ^ 2 = 4 5从 y 到 W 的距离是 √(45) = 3√(5) .
本节最后的定理说明,关于 proj_(w)y 的公式(2)当 W 的基是单位正交集时如何被简化.
如果 {u₁,…,uₚ} 是 Rⁿ 中子空间 W 的单位正交基,那么
operatornameproj _ W boldsymbol y = (boldsymbol y · boldsymbol u _ 1) boldsymbol u _ 1 + (boldsymbol y · boldsymbol u _ 2) boldsymbol u _ 2 + … + (boldsymbol y · boldsymbol u _ p) boldsymbol u _ p (4)如果 U = [u₁ u₂ … uₚ] ,则
proj _ w ±b y = U U ^ T ±b y, text 对所有 ±b y ∈ mathbb R ^ n text 成立 (5)证明
公式(4)可由定理8中(2)式立刻推出,而且,(4)式说明 proj_(w)y 是 U 的列的线性组合,对应权分别为 y · u₁, y · u₂, …, y · uₚ 。这些权同样可写成 u₁^(T)y, u₂^(T)y, …, uₚ^(T)y 。表明它们是
U^(T)y 中的元素,从而证明了(5)式.
若 U 是列单位正交的 n × p 矩阵,W 是 U 的列空间,那么
U^(T)U±bx = Iₚ±bx = ±bx ,对所有属于 ℝ^p 的 ±bx 成立 定理6
UU^(T)±b y = proj_(W)±by ,对所有属于 ℝⁿ 的 ±by 成立 定理10
如果 U 是 n × n (方)阵,且列单位正交,那么 U 是正交矩阵,列空间 W 是全部 Rⁿ 空间,且 UU^(T)y = Iy = y ,对所有 y ∈ Rⁿ 成立.
尽管公式(4)在理论上十分重要,实际上,它常包含数字的平方根运算(在 uᵢ 的元素中),用手工运算时我们推荐使用公式(2).
-
令 ±bu₁ = -7 1 4 , ±bu₂ = -1 1 -2 , ±b y = -9 1 6 ,且 W = Span{±bu₁,±bu₂} ,利用 ±bu₁ 和 ±bu₂ 正交计算projw y.
-
设 W 是 Rⁿ 的子空间,x 和 y 是 Rⁿ 中的向量,且 z = x + y 。如果 u 是 x 在 W 上的投影,v 是 y 在 W 上的投影。证明 u + v 是 z 在 W 上的投影。
在习题 1 和习题 2 中,我们假设 {u₁,…,u₄} 是 R⁴ 中的正交基.
1. boldsymbol u _ 1 = [ c 0 1 - 4 - 1 ], boldsymbol u _ 2 = [ c 3 5 1 1 ], boldsymbol u _ 3 = [ c 1 0 1 - 4 ], boldsymbol u _ 4 = [ c 5 - 3 - 1 1 ],x=10-820.将 dotx 写成两个向量之和,一个属于
Span {u₁,u₂,u₃} ,另一个属于Span {u₄} .
2. boldsymbol u _ 1 = [ l 1 2 1 1 ], boldsymbol u _ 2 = [ c - 2 1 - 1 1 ], boldsymbol u _ 3 = [ c 1 1 - 2 - 1 ], boldsymbol u _ 4 = [ c - 1 1 1 - 2 ]v=4 5 -3 3.将 v 写成两个向量之和,一个属于
boldsymbol y = [ c - 1 4 3 ], boldsymbol u _ 1 = [ c 1 1 0 ], boldsymbol u _ 2 = [ c - 1 1 0 ]Span {u₁} ,另一个属于 Span {u₂, u₃, u₄} .
并找出 y 在 operatornameSpan{u₁, u₂} 上的正交投影.
在习题 3~6 中,验证 {u₁, u₂} 是一个正交集,
boldsymbol y = [ c 6 3 - 2 ], boldsymbol u _ 1 = [ c 3 4 0 ], boldsymbol u _ 2 = [ c - 4 3 0 ]mathbf y = [ c - 1 2 6 ], mathbf u _ 1 = [ c 3 - 1 2 ], mathbf u _ 2 = [ c 1 - 1 - 2 ]mathbf y = [ l 6 4 1 ], mathbf u _ 1 = [ l - 4 - 1 1 ], mathbf u _ 2 = [ l 0 1 1 ]在习题 7~10 中,设 W 是由 u 生成的子空间,将 y 写成 W 中的向量与正交于 W 的向量之和.
mathbf y = [ l 1 3 5 ], mathbf u _ 1 = [ c 1 3 - 2 ], mathbf u _ 2 = [ l 5 1 4 ]在习题 11 和习题 12 中,在 v₁ 和 v₂ 所生成的子空间 W 中,找出距离 y 最近的点.
1 1. mathbf y = [ l 3 1 5 1 ], mathbf v _ 1 = [ c 3 1 - 1 1 ], mathbf v _ 2 = [ l 1 - 1 1 - 1 ]1 2. boldsymbol y = [ r 3 - 1 1 1 3 ], boldsymbol v _ 1 = [ r 1 - 2 - 1 2 ], boldsymbol v _ 2 = [ r - 4 1 0 3 ]在习题 13 和习题 14 中,在形如 c₁v₁ + c₂v₂ 的向量中找出最接近 z 的向量.
mathbf z = [ c 3 - 7 2 3 ], mathbf v _ 1 = [ c 2 - 1 - 3 1 ], mathbf v _ 2 = [ c 1 1 0 - 1 ]1 4. z = [ c 2 4 0 - 1 ], v _ 1 = [ c 2 0 - 1 - 3 ], v _ 2 = [ c 5 - 2 4 2 ]- 假设 y=5 -9 5 , u₁=-3 -5 1 , u₂=-3 2 1 ,找出向
-
若 y, v₁ 和 v₂ 如习题 12 中所示,找出 y 到 R⁴ 中由 v₁ 和 v₂ 所生成的子空间的距离.
-
设 y=4 8 1 , u₁=2/3 1/3 2/3 , u₂=-2/3 2/3 1/3 ,且 W = Span{u₁, u₂} . a. 若 U = [u₁ u₂] ,计算 U^(T)U 和 UU^(T) . b. 计算 proj_(w)y 和 (UU^(T))y .
-
设 y=7 9, u₁=1/√(10) -3/√(10) 和 W=operatornameSpan{±bu₁} . a. 若 U 是 2 × 1 矩阵,其列向量为 u₁ ,计算 U^(T)U 和 UU^(T) . b. 计算 proj_(w)y 和 (UU^(T))y .
-
假设 u₁=1 1 -2,u₂=5 -1 2 和 u₃=0 0 1 ,注意 u₁ 和 u₂ 正交,但 u₃ 和 u₁ 或 u₂ 不正交,可以证明 u₃ 不属于 u₁ 和 u₂ 所生成的子空间。利用这个事实,构造 R³ 中与 u₁ 和 u₂ 正交的非零向量 v。
-
设 u₁ 和 u₂ 如习题 19 所示,且 u₄=0 1 0 ,可以证明 u₄ 不属于 u₁ 和 u₂ 所生成的子空间 W.利用这个结论构造 R³ 中与 u₁ 和 u₂ 正交的非零向量 v. 在习题 21~22 中,所有向量和子空间均属于 Rⁿ ,判断下列结论的正误,验证每个结论.
-
a.如果 z 与 ±bu₁ 和 ±bu₂ 正交,且 W = operatornameSpan{±bu₁, ±bu₂} ,那么 z 一定属于 W^perp .
b.对任一向量 y 和任一子空间 W,向量 y - proj_(w)y 正交于 W.
c. y 在子空间 W 上的正交投影 haty 的计算有时可能依赖于 W 正交基的选取.
d.如果 y 属于子空间 W,那么 y 在 W 上的正交投影是 y 本身.
e.如果 n × p 矩阵 U 的列是单位正交的,那么 UU^(T)y 是 y 在 U 的列空间上的正交投影.
- a.如果 W 是 Rⁿ 的子空间,且 v 属于 W 和 W^(perp) ,那么 v 一定是零向量.
b.在正交分解定理中,公式(2)中 haty 的每一项在 W 的子空间上的正交投影是它本身.
c.如果 y = z₁ + z₂ ,其中 z₁ 属于子空间 W 且 z₂ 属于子空间 W^(perp) ,那么 z₁ 一定是 y 在 W 上的正交投影.
d.子空间 W 中的元素对向量 y 的最佳逼近是向量 y - proj_(w)y .
e.如果一个 n × p 矩阵 U 有单位正交列,那么
UU^(T)±b x = ±bx 对所有 ±b x∈ ℝⁿ 成立.
-
若 A 是 m × n 矩阵,证明 ℝⁿ 中任一向量 x 可以写成 x = p + u 的形式,其中 p 属于Row A 且 u 属于Nul A 。同样可以证明:如果方程 Ax = b 是相容的,那么存在唯一向量 p 属于Row A ,使得 Ap = b 。
-
假若 W 是 Rⁿ 的一个子空间且 {w₁,…,wₚ} 是 W 的一个正交基, {v₁,…,v_(q)} 是 W^(perp) 的正交基.
a.说明为什么 {±bw₁, …, ±bwₚ, ±bv₁, …, ±bv_q} 是正交集.
b. 说明为什么(a)中的集合可以生成 ℝⁿ .
c. 证明 dim W + dim W^perp = n .
-
[M] 若 U 是6.2节中习题36中的 8 × 4 矩阵,在Col U 中找出最接近 y = (1,1,1,1,1,1,1,1,1) 的点,写出你解决该问题所使用的键击内容或命令.
-
[M] 假设 U 是习题 25 中的矩阵,求出从 ±bb = (1,1,1,1,-1,-1,-1,-1) 到 operatornameCol U 的距离.
- 计算
在这种情形下,y恰好是 u₁ 和 u₂ 的线性组合,从而y属于W,即W中最接近y的点是y本身.
- 利用定理10,设 U 是一个矩阵,它的列向量构成 W 的一组标准正交基.则