形式陈述
设 v 1 , … , v k 是内积空间中的线性无关向量。先令
w 1 = v 1 , u 1 = w 1 ‖ w 1 ‖ , 再对 j = 2 , … , k 递归定义
w j = v j − ∑ i < j ⟨ u i , v j ⟩ u i , u j = w j ‖ w j ‖ . 这里采用本页第一变量共轭线性、第二变量线性的约定,因此 ⟨ u i , v j ⟩ u i 正是 v j 到方向 u i 的正交投影 公理库 正交投影 Orthogonal projection 把向量映到子空间上最近点并使误差与子空间正交的线性算子。 。所得 ( u 1 , … , u k ) 是正交规范组,并且对每个 j 都有
span ( u 1 , … , u j ) = span ( v 1 , … , v j ) . 对 i < j ,将残差与 u i 配对,正交性给出 ⟨ u i , w j ⟩ = ⟨ u i , v j ⟩ − ⟨ u i , v j ⟩ = 0 。若 w j = 0 ,则 v j 已由前面的 u i 、也就是前面的 v i 张成,与线性无关矛盾。因此每一步都得到非零正交残差,可以归一化。
直觉
每一步都把新向量拆成两部分:一部分已经能由旧方向解释,另一部分是旧子空间从未见过的新方向。把前者逐项减去便得到残差 w j ;它与所有旧的 u i 正交,归一化后成为新的单位坐标轴。这个过程只删除重复方向,没有删除新信息,所以每一步的张成空间都保持不变。
图片加载失败 Gram–Schmidt 的投影消去
例子与边界
由 v 1 = ( 1 , 1 ) 、v 2 = ( 1 , 0 ) 出发,先得到
w 1 = ( 1 , 1 ) , u 1 = ( 1 , 1 ) 2 . 第二步为
w 2 = v 2 − ⟨ u 1 , v 2 ⟩ u 1 = ( 1 , 0 ) − 1 2 ( 1 , 1 ) 2 = ( 1 2 , − 1 2 ) , 所以 u 2 = ( 1 , − 1 ) / 2 。若输入中含有依赖向量,则相应残差为零:该向量已在先前的张成空间内。跳过零残差、继续处理后续向量,仍可为全部输入的张成空间提取正交规范基。
上述递推是经典 Gram–Schmidt:先用原向量 v j 算出全部投影系数,再一起减去。浮点运算中,相近向量相减会放大相对误差,使残差偏离精确的正交方向。改进 Gram–Schmidt 则令 r = v j ,每次用当前残差计算 ⟨ u i , r ⟩ 并立即更新 r ;这种逐次消去通常更能保持正交性。稠密矩阵的数值分解还常用Householder 反射 公理库 Householder 反射 Householder reflection · Householder transformation 用一个向量紧凑表示酉反射,把整段向量化到单一坐标方向,并作为稳定 QR 的基本原语。 构造QR 分解 公理库 QR 分解 QR factorization · QR decomposition · Economy-size QR 把长方矩阵分解为正交列与上三角因子,并区分经济型表示、数值算法和秩亏边界。 。
推论与应用
该过程为线性无关 公理库 线性无关 Linear independence 只有全零系数能产生零向量的向量组。 组的张成空间构造正交规范基 公理库 正交规范基 Orthonormal basis 由单位长度且两两正交的向量组成的基。 。把 k 个输入向量排成 A ∈ F m × k 的列,输出排成 Q ∈ F m × k 的列,就有 A = Q R 。这里 R ∈ F k × k 为上三角矩阵:第 j 列的上方条目是投影系数 ⟨ u i , v j ⟩ ,对角条目是残差长度 ‖ w j ‖ ,下方条目为零。
在满列秩的最小二乘 公理库 最小二乘与正规方程 Least squares · Normal equations 将目标向量正交投影到矩阵列空间,并以残差正交条件导出正规方程。 问题中,把 b 投影到 Q 的列空间后,求解上三角方程 R x = Q ∗ b 即可得到最优系数。这条QR 计算路线 公理库 用 QR 与 SVD 求最小二乘 Least squares via QR · Least squares via SVD · Numerical least squares 以 QR 作为满列秩最小二乘的默认计算路线,并用 SVD 处理秩亏、欠定和最小范数解。 直接使用正交坐标,保留了投影的几何结构。
对 1 , x , x 2 , … 按积分内积依次正交化,则得到正交多项式:第 j 步减去低次多项式方向,保留一个与它们全部正交的新多项式。
参考资料