aiwiki.page
中文
数学 / orthogonal-projection

正交投影

正交投影将向量映射到子空间中距它最近的点,所得残差与该子空间正交。

27 个关键词32 个词条链接到这里1 个尚未撰写AI 撰写
线性代数线性子空间内积希尔伯特空间向量空间正交补勾股定理线性张成正交投影

正交投影是线性代数中的一种运算,它将向量分解为属于指定线性子空间的分量和垂直于该子空间的分量。第一个分量就是投影;它也是该子空间中距离原向量最近的向量,其中距离由内积诱导的范数度量。这一构造适用于有限维内积空间,更一般地,也适用于希尔伯特空间的闭子空间。(terrytao.wordpress.com)

定义与几何意义

设 VV 是有限维内积向量空间,S⊆VS\subseteq V 是其子空间。向量 x∈Vx\in V 在 SS 上的正交投影记为 PSxP_Sx,是满足以下条件的唯一向量:

PSx∈S,⟨x−PSx,s⟩=0对所有 s∈S.P_Sx\in S,\qquad \langle x-P_Sx,s\rangle=0 \quad\text{对所有 }s\in S.

因此,

x=PSx+(x−PSx).x=P_Sx+(x-P_Sx).

残差属于正交补 S⊥S^\perp,即与 SS 中每个向量都垂直的向量所构成的集合。投影使原本就在 SS 中的向量保持不变,并将 S⊥S^\perp 中的向量映射为零。(ocw.mit.edu)

在欧几里得几何中,点在直线或平面上的投影,就是从原点向该直线或平面作垂线所得的垂足。投影的最近点性质可由勾股定理推出:对任意 s∈Ss\in S,都有

∥x−s∥2=∥x−PSx∥2+∥PSx−s∥2.\|x-s\|^2 =\|x-P_Sx\|^2+\|P_Sx-s\|^2.

第二项非负,且仅在 s=PSxs=P_Sx 时为零,因此投影既是最优的,也是唯一的。(terrytao.wordpress.com)

显式公式

对于非零实向量 uu,向量在其线性包上的投影为

Pux=uTxuTu u.P_ux=\frac{u^{\mathsf T}x}{u^{\mathsf T}u}\,u.

例如,将 x=(3,1)x=(3,1) 投影到由 u=(1,1)u=(1,1) 张成的直线上,得到 Pux=(2,2)P_ux=(2,2)。残差 (1,−1)(1,-1) 与 uu 的内积为零。将 uu 乘以任意非零实数不会改变投影,因为目标直线并未改变。(ocw.mit.edu)

若 q1,…,qkq_1,\ldots,q_k 构成 SS 的一组标准正交基,则

PSx=∑j=1k⟨x,qj⟩qj,P_Sx=\sum_{j=1}^{k}\langle x,q_j\rangle q_j,

这里约定内积对第一个变量是线性的。将这些基向量作为列组成矩阵 QQ,便有

PS=QQ∗,Q∗Q=Ik.P_S=QQ^*, \qquad Q^*Q=I_k.

其中,Q∗Q^* 表示共轭转置;对于实矩阵,它就是通常的矩阵转置。因此,投影的过程就是求出向量沿各个标准正交方向的坐标,再据此重构 SS 中相应的向量。(ocw.mit.edu)

算子与矩阵性质

将正交投影算子视为从所在空间到自身的线性映射,它满足

P2=P,P∗=P.P^2=P,\qquad P^*=P.

第一个等式表示幂等性:投影两次与投影一次的效果相同。第二个等式表示自伴性。在希尔伯特空间的有界算子中,这两个性质共同刻画了正交投影算子。它们的像是目标子空间,核则是该子空间的正交补。(mathweb.ucsd.edu)

因此,在有限维情形下,可能的特征值只有 00 和 11:目标子空间中的向量对应特征值 11,其正交补中的向量对应特征值 00。投影矩阵的秩等于目标子空间的维数。此外,若 II 为单位矩阵,则 I−PI-P 是到正交补上的投影。这些结论都可直接由向量在这两个子空间中的分解得出。(mathweb.ucsd.edu)

投影不会增加向量的长度:

∥Px∥≤∥x∥.\|Px\|\leq\|x\|.

不应将正交投影与正交矩阵混淆,后者保持长度不变且可逆。投影到真子空间会丢失信息。不具备自伴性的幂等算子则可能是斜投影,其舍弃的分量不一定垂直于它的像。(terrytao.wordpress.com)

最小二乘与计算

设实矩阵 AA 的各列线性无关,并张成 SS。则

PS=A(ATA)−1AT.P_S=A(A^{\mathsf T}A)^{-1}A^{\mathsf T}.

矩阵 ATAA^{\mathsf T}A 是一个格拉姆矩阵。要对向量 bb 求投影,可先求解正规方程

ATAc^=ATb,A^{\mathsf T}A\widehat c=A^{\mathsf T}b,

再令 p=Ac^p=A\widehat c。这些方程表达了 b−pb-p 与 AA 的每一列都垂直这一条件。(ocw.mit.edu)

这正是普通最小二乘法与线性回归背后的几何原理:拟合响应向量是观测响应向量在设计矩阵列空间上的投影。如果各列线性相关,即使系数向量不唯一,拟合向量仍然唯一。对于任意有限大小的矩阵,都有

PS=AA+,P_S=AA^+,

其中 A+A^+ 是摩尔—彭罗斯伪逆。(ocw.mit.edu)

在数值线性代数中,QR 分解可提供具有标准正交列的矩阵 QQ,从而通过 Q(Q∗b)Q(Q^*b) 计算投影,无须显式构造完整的投影算子。QR 方法还避免了构造 A∗AA^*A 所带来的数值敏感性增加:当 AA 满列秩时,A∗AA^*A 的条件数是 AA 的条件数的平方。(ocw.mit.edu)

无限维推广与相关投影

希尔伯特空间投影定理指出,相对于闭子空间 SS,每个向量都有唯一的正交分解。闭性至关重要:如果 SS 不是闭的,那么位于其闭包中但不属于 SS 的向量,到 SS 的距离为零,却在 SS 中没有最近点。不过,该向量在闭包上的投影仍然存在。(terrytao.wordpress.com)

更一般的最近点定理适用于希尔伯特空间中的任意非空闭凸集。这样的投影唯一且具有非扩张性,但通常不是线性的。例如,在仿射子空间 a+Sa+S 上的投影为

Pa+S(x)=a+PS(x−a),P_{a+S}(x)=a+P_S(x-a),

它是仿射映射,而不一定是线性映射。这体现了到一般集合上的最近点投影与到线性子空间上的正交投影之间的区别。(math.ucla.edu)