正规方程组与最小二乘解

设 A 是 m\times n 实矩阵,b 是 m 维向量。最小二乘问题是

\min_x \|Ax-b\|_2^2

对应的正规方程组(也叫法方程)是

A^T A x = A^T b.

其中 A^T 是 A 的转置。


1. 为什么正规方程组一定有解

我们知道,一般的线性方程组 Ax=b 不一定有解,这时候我们就去找它的最小二乘解,最小二乘解对应的就是下面这个正规方程组的解
A^T A x = A^T b

验证这件事,就需要说明r(A^T A) = r(A^T b) 换言之,A^T b 一定落在矩阵 A^T A 的列空间中。

证明: 若 z\in N(A^T A),即 A^T A z=0,则

0=z^T A^T A z=(Az)^T(Az)=\|Az\|^2,

所以 Az=0。于是

z^T A^T b=(Az)^T b=0.

这说明 A^T b 与 N(A^T A) 正交,因此

A^T b\in N(A^T A)^\perp=\operatorname{Col}(A^T A).

所以存在 x,使得

A^T A x=A^T b.

也就是说,正规方程组一定有解。

不过要注意:有解不一定唯一。
如果 A 的列向量线性无关,即 A 列满秩,那么 A^T A 可逆,解唯一:

x=(A^T A)^{-1}A^T b.

如果 A 列不满秩,则 A^T A 奇异,正规方程有无穷多解,但它们给出的 Ax 是同一个向量。


2. 正规方程组的解七号是最小二乘解

对于最小二乘解,我们要最小化残差平方和

f(x)=\|Ax-b\|_2^2.

将这个表达式展开得到:

f(x)=(Ax-b)^T(Ax-b) = x^T A^T A x-2x^T A^T b+b^T b.

对 x 求梯度:

\nabla f(x)=2A^T A x-2A^T b.

令梯度为零:

2A^T A x-2A^T b=0,

得到

A^T A x=A^T b.

这就是正规方程组。

由于 A^T A 是半正定矩阵,f(x) 是凸二次函数,所以满足梯度为零的点就是全局最小值点。


我们看下面的例子:

用数据 (0,1),(1,2),(2,2) 拟合一条直线

y=c+dt.

设计矩阵和向量为

A= \begin{bmatrix} 1&0\\ 1&1\\ 1&2 \end{bmatrix}, \qquad b= \begin{bmatrix} 1\\ 2\\ 2 \end{bmatrix}, \qquad x= \begin{bmatrix} c\\ d \end{bmatrix}.

计算:

A^T A= \begin{bmatrix} 1&1&1\\ 0&1&2 \end{bmatrix} \begin{bmatrix} 1&0\\ 1&1\\ 1&2 \end{bmatrix} = \begin{bmatrix} 3&3\\ 3&5 \end{bmatrix},

A^T b= \begin{bmatrix} 1&1&1\\ 0&1&2 \end{bmatrix} \begin{bmatrix} 1\\ 2\\ 2 \end{bmatrix} = \begin{bmatrix} 5\\ 6 \end{bmatrix}.

正规方程组为

\begin{bmatrix} 3&3\\ 3&5 \end{bmatrix} \begin{bmatrix} c\\ d \end{bmatrix} = \begin{bmatrix} 5\\ 6 \end{bmatrix}.

即

3c+3d=5,

3c+5d=6.

两式相减得

2d=1,\qquad d=\frac12.

代回:

3c+3\cdot \frac12=5 \Rightarrow 3c+\frac32=5 \Rightarrow 3c=\frac72 \Rightarrow c=\frac76.

所以最小二乘直线为

y=\frac76+\frac12 t.

预测值:

t=0:\ y=\frac76,

t=1:\ y=\frac76+\frac12=\frac53,

t=2:\ y=\frac76+1=\frac{13}{6}.

残差:

r=b-Ax= \begin{bmatrix} 1-\frac76\\ 2-\frac53\\ 2-\frac{13}{6} \end{bmatrix} = \begin{bmatrix} -\frac16\\ \frac13\\ -\frac16 \end{bmatrix}.

验证残差与列空间正交:

A^T r= \begin{bmatrix} 1&1&1\\ 0&1&2 \end{bmatrix} \begin{bmatrix} -\frac16\\ \frac13\\ -\frac16 \end{bmatrix} = \begin{bmatrix} 0\\ 0 \end{bmatrix}.

因此正规方程给出的 c=\frac76,d=\frac12 确实使残差平方和最小。

残差平方和为

\left(-\frac16\right)^2+\left(\frac13\right)^2+\left(-\frac16\right)^2 = \frac1{36}+\frac4{36}+\frac1{36} = \frac16.


总结:

  • 正规方程组是 A^T A x=A^T b。
  • 它一定有解,因为 A^T b 属于 A^T A 的列空间。
  • 它的解使 \|Ax-b\|^2 最小,因为残差与 A 的列空间正交。
  • 若 A 列满秩,解唯一;否则解不唯一,但 Ax 是唯一的最小二乘投影。
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

相关阅读更多精彩内容

友情链接更多精彩内容