跳到正文
格致开物MATHWIKI

最小范数解

最小范数解是在给定方程的全部解中,使指定范数最小的解。在线性代数中,通常指相容线性方程 Ax=b 的最小欧几里得范数解,即求 min⁡{‖x‖2:Ax=b},‖x‖2=x𝖳x. 必须说明使用哪种范数、方程是否要求精确成立。方程不相容时,“残差最小的解中再取范数最小者”是最小范数最小二乘解,属于另一层选择,详见伪逆。[1]

考虑一个教学约束 x1+2x2=5。它有无穷多解,例如 (5,0)、(1,2)、(−1,3)。若两个坐标使用相同尺度,最小范数要求从这条解直线上挑出离原点最近的点。

一条直线上的完整计算

令 x2=t,则 x1=5−2t。长度平方为 ‖x‖22=(5−2t)2+t2=5(t−2)2+5. 平方项在 t=2 时唯一为零,所以最小范数解为 x∗=(1,2)𝖳,最小范数是 5。把它代回约束,1+2⋅2=5,确实是精确解。

解直线x一加二倍x二等于五,原点到一二的线段与直线垂直,另一解三一位于标出最短距离的圆外。
所有点都满足同一约束;从原点出发的垂线选出最短的那个解。

直线方向可取 z=(−2,1)𝖳,而 x∗𝖳z=−2+2=0。任一解都可写成 x∗+sz,正交性给出 ‖x∗+sz‖22=‖x∗‖22+s2‖z‖22=5+5s2. 这不是图形测量,而是全体解上的精确比较。

核空间的正交分解

定理 若有限维实线性方程 Ax=b 相容,则它有唯一的最小欧几里得范数解;该解是全部解中唯一垂直于 ker⁡A 的向量。

证明先取任一解 x0。两个解之差在核空间中,故全体解为 x0+ker⁡A。把 x0 正交分解为 x∗+z0,其中 z0∈ker⁡A、x∗⊥ker⁡A。因为 Az0=0,Ax∗=Ax0=b,所以 x∗ 仍是解。

任意其他解可写为 x∗+z,z∈ker⁡A。由勾股等式 ‖x∗+z‖22=‖x∗‖22+‖z‖22≥‖x∗‖22, 等号只在 z=0 成立,所以最小者存在且唯一。证明真正使用的是有限维内积空间的正交分解;如果更换范数,这个正交判别不能照搬。

满行秩时的显式公式

设 A∈ℝm×n 满行秩,因而 m≤n。矩阵 AA𝖳 正定:对非零 y∈ℝm, y𝖳AA𝖳y=‖A𝖳y‖22>0, 其中严格大于零来自 A𝖳 的列线性无关。因此 AA𝖳 可逆,且 x∗=A𝖳(AA𝖳)−1b. 左乘 A 得 Ax∗=b;若 Az=0,则 x∗𝖳z=b𝖳(AA𝖳)−1Az=0,所以它满足刚证明的最小性判别。这里用的是行秩条件,而不是要求 A 为可逆方阵。

原例中 A=(1  2),AA𝖳=5,故 x∗=(1,2)𝖳⋅5/5=(1,2)𝖳。变式将右端换成十,解线平移为 x1+2x2=10,最短解同比变为 (2,4),范数变为 25。

若行不独立,例如 A=(1224),则 AA𝖳 奇异,不能使用上述逆矩阵公式;但相容右端 b=(5,10)𝖳 仍有唯一最小范数解 (1,2)。定理未失效,失效的是满行秩下的特殊表达式。一般情况可通过奇异值分解或伪逆处理。

范数反映了坐标的尺度

最小欧几里得范数隐含两个坐标在同一尺度下比较。若第一坐标每单位的代价更高,可以改为最小化 4x12+x22。对同一约束令 x1=5−2t,目标为 4(5−2t)2+t2=17(t−4017)2+10017. 于是解变为 (5/17,40/17),而不是 (1,2)。这不是前一次计算错误,而是长度的定义改变了。

更一般地,对正定权重矩阵 W,满行秩约束的最小 W-范数解为 W−1A𝖳(AW−1A𝖳)−1b;它可以通过变量替换 u=W1/2x 转回欧几里得范数推得。若使用 ℓ1 或 ℓ∞ 范数,唯一性也可能消失。例如约束 x1+x2=1 下,所有非负解的 ℓ1 范数都为一,整条线段都是最小者。

参考文献

  1. ↑ Stephen Boyd,Least-norm solutions of underdetermined equations,Stanford EE263 讲义。