我们从一个很小的例子开始。
设
$$ A= \begin{bmatrix} 1 & 2\\ 2 & 4\\ 3 & 6 \end{bmatrix} $$这个矩阵有两列:
$$ a_1= \begin{bmatrix} 1\\2\\3 \end{bmatrix}, \qquad a_2= \begin{bmatrix} 2\\4\\6 \end{bmatrix} $$显然
$$ a_2=2a_1 $$也就是说,第二列并没有提供新的方向。虽然 $A$ 有两列,但它真正独立的列方向只有一个。
这就是“不满列秩”的最简单例子。
这件事看似只是“有一列是另一列的倍数”,但它会引出一串重要结论:
$$ A\text{ 不满列秩} \Longleftrightarrow \exists z\neq 0,\ Az=0 \Longleftrightarrow A^TA\text{ 奇异} \Longleftrightarrow A^TA\text{ 不可逆} \Longleftrightarrow A^TA\text{ 有零特征值} \Longleftrightarrow A^TA\text{ 半正定但不正定} $$这些结论正是理解普通最小二乘、Tikhonov 正则化、Moore-Penrose 伪逆的基础。
向量线性组合
设
$$ A\in\mathbb{R}^{m\times n} $$也就是说,$A$ 有 $m$ 行、$n$ 列。把它按列写成:
$$ A= \begin{bmatrix} | & | & & |\\ a_1 & a_2 & \cdots & a_n\\ | & | & & | \end{bmatrix} $$对于任意
$$ x= \begin{bmatrix} x_1\\x_2\\ \vdots\\x_n \end{bmatrix} $$矩阵乘法
$$ Ax $$本质上是列向量的线性组合:
$$ Ax=x_1a_1+x_2a_2+\cdots+x_na_n $$所以矩阵 $A$ 可以理解为一个线性变换:
$$ A:\mathbb{R}^n\to\mathbb{R}^m $$它把参数空间中的 $x$,映射到输出空间中的 $Ax$。
矩阵的基本空间
理解矩阵性质,必须先理解四个空间。
列空间
列空间记为
$$ \mathcal{C}(A) $$定义为
$$ \mathcal{C}(A)=\{Ax:x\in\mathbb{R}^n\} $$也就是 $A$ 的所有列向量的线性组合构成的空间。
换句话说:
$$ \mathcal{C}(A)=\operatorname{span}\{a_1,a_2,\dots,a_n\} $$它告诉我们:所有可能的输出 $Ax$ 落在哪个空间里。
在最小二乘问题
$$ Ax\approx b $$中,真正能被 $Ax$ 表示出来的向量,只能落在 $\mathcal{C}(A)$ 中。
零空间
零空间记为
$$ \mathcal{N}(A) $$定义为
$$ \mathcal{N}(A)=\{x:Ax=0\} $$它告诉我们:哪些输入方向会被 $A$ 压成零。
如果存在非零向量
$$ z\neq 0 $$使得
$$ Az=0 $$那么说明 $A$ 在方向 $z$ 上丢失了信息。
因为对于任意 $x$,都有
$$ A(x+z)=Ax+Az=Ax $$也就是说,$x$ 和 $x+z$ 经过 $A$ 之后完全一样。
这就是“不唯一”的根源。
行空间
行空间是 $A$ 的行向量张成的空间。它也等于 $A^T$ 的列空间:
$$ \mathcal{R}(A)=\mathcal{C}(A^T) $$它位于输入空间 $\mathbb{R}^n$ 中。
行空间与零空间互相正交:
$$ \mathcal{R}(A)\perp \mathcal{N}(A) $$更准确地说:
$$ \mathbb{R}^n=\mathcal{R}(A)\oplus \mathcal{N}(A) $$这里的 $\oplus$ 表示直和分解。
这句话很重要:任意一个输入向量 $x$,都可以拆成两部分:
$$ x=x_r+x_n $$其中
$$ x_r\in\mathcal{R}(A),\qquad x_n\in\mathcal{N}(A) $$而
$$ Ax=A x_r+A x_n=A x_r $$所以零空间分量 $x_n$ 对输出没有任何影响。
左零空间
左零空间是
$$ \mathcal{N}(A^T) $$定义为
$$ \mathcal{N}(A^T)=\{y:A^Ty=0\} $$它位于输出空间 $\mathbb{R}^m$ 中。
列空间和左零空间互相正交:
$$ \mathcal{C}(A)\perp \mathcal{N}(A^T) $$并且
$$ \mathbb{R}^m=\mathcal{C}(A)\oplus \mathcal{N}(A^T) $$普通最小二乘中的正交投影,就是在这个分解下发生的。
什么是秩
矩阵的秩记为
$$ \operatorname{rank}(A) $$它可以理解为 $A$ 真正保留下来的独立方向数。
等价地说:
$$ \operatorname{rank}(A)=\dim \mathcal{C}(A) $$也就是列空间的维数。
虽然矩阵 $A$ 有 $n$ 列,但这些列不一定都独立。
如果
$$ \operatorname{rank}(A)=n $$说明 $A$ 的 $n$ 个列向量全部线性无关,此时称 $A$ 列满秩。
如果
$$ \operatorname{rank}(A) < n $$说明 $A$ 的列向量之间存在冗余,此时称 $A$ 不满列秩。
注意:列满秩要求
$$ \operatorname{rank}(A)=\text{列数} $$而不是要求 $A$ 是方阵。
例如
$$ A\in\mathbb{R}^{100\times 3} $$只要三列线性无关,它就是列满秩。
列满秩和零空间的关系
核心结论:
$$ A\text{ 列满秩} \Longleftrightarrow \mathcal{N}(A)=\{0\} $$也就是说:
$$ \operatorname{rank}(A)=n \Longleftrightarrow Az=0\text{ 只有零解} $$反过来:
$$ A\text{ 不满列秩} \Longleftrightarrow \exists z\neq 0,\ Az=0 $$证明很直接。
如果 $A$ 不满列秩,那么列向量线性相关,所以存在一组不全为零的系数 $z_1,\dots,z_n$,使得
$$ z_1a_1+\cdots+z_na_n=0 $$写成矩阵形式就是
$$ Az=0 $$其中
$$ z= \begin{bmatrix} z_1\\ \vdots\\ z_n \end{bmatrix} \neq 0 $$所以不满列秩推出存在非零零空间向量。
反过来,如果存在非零 $z$,使得
$$ Az=0 $$那么
$$ z_1a_1+\cdots+z_na_n=0 $$这就是列向量之间存在非平凡线性关系,因此列向量线性相关,所以
$$ \operatorname{rank}(A) < n $$于是二者等价。
秩-零化度定理
对任意
$$ A\in\mathbb{R}^{m\times n} $$有
$$ \operatorname{rank}(A)+\dim\mathcal{N}(A)=n $$这叫秩-零化度定理。
它的含义是:
输入空间一共有 $n$ 个维度,其中一部分被 $A$ 有效映射出去,另一部分被 $A$ 压成零。
如果 $A$ 列满秩:
$$ \operatorname{rank}(A)=n $$那么
$$ \dim\mathcal{N}(A)=0 $$所以
$$ \mathcal{N}(A)=\{0\} $$如果 $A$ 不满列秩:
$$ \operatorname{rank}(A) < n $$那么
$$ \dim\mathcal{N}(A)>0 $$所以存在非零零空间方向。
奇异性和可逆性
严格来说,“奇异”和“可逆”只对方阵讨论。
设
$$ M\in\mathbb{R}^{n\times n} $$如果存在矩阵 $M^{-1}$,使得
$$ M^{-1}M=MM^{-1}=I $$那么 $M$ 可逆。
如果 $M$ 不可逆,则称 $M$ 奇异。
所以:
$$ M\text{ 奇异} \Longleftrightarrow M\text{ 不可逆} $$对于方阵 $M$,下面这些条件等价:
$$ M\text{ 可逆} $$$$ \Longleftrightarrow \mathcal{N}(M)=\{0\} $$$$ \Longleftrightarrow Mx=0\text{ 只有零解} $$$$ \Longleftrightarrow \operatorname{rank}(M)=n $$$$ \Longleftrightarrow \det(M)\neq 0 $$$$ \Longleftrightarrow 0\text{ 不是 }M\text{ 的特征值} $$反过来,下面这些条件也等价:
$$ M\text{ 奇异} $$$$ \Longleftrightarrow M\text{ 不可逆} $$$$ \Longleftrightarrow \exists z\neq 0,\ Mz=0 $$$$ \Longleftrightarrow \operatorname{rank}(M) < n $$$$ \Longleftrightarrow \det(M)=0 $$$$ \Longleftrightarrow 0\text{ 是 }M\text{ 的特征值} $$为什么有零特征值就不可逆
特征值定义为:
$$ Mv=\lambda v $$其中
$$ v\neq 0 $$如果
$$ \lambda=0 $$那么
$$ Mv=0 $$也就是说,存在非零向量 $v$,被 $M$ 压成零。
如果 $M$ 可逆,则两边左乘 $M^{-1}$:
$$ M^{-1}Mv=M^{-1}0 $$得到
$$ v=0 $$这和 $v\neq 0$ 矛盾。
所以:
$$ 0\text{ 是特征值} \Longrightarrow M\text{ 不可逆} $$反过来,如果 $M$ 不可逆,那么存在非零 $v$,使得
$$ Mv=0 $$这正是
$$ Mv=0\cdot v $$所以 $0$ 是特征值。
因此:
$$ M\text{ 不可逆} \Longleftrightarrow 0\text{ 是 }M\text{ 的特征值} $$为什么 $A^TA$ 特别重要
普通最小二乘中,我们要求
$$ \min_x\|Ax-b\|^2 $$目标函数为
$$ J(x)=\|Ax-b\|^2 $$求导得到法方程:
$$ A^TAx=A^Tb $$所以 $A^TA$ 是普通最小二乘中的核心矩阵。
另一方面,$A^TA$ 还有明确的几何意义。
如果
$$ A= \begin{bmatrix} | & | & & |\\ a_1 & a_2 & \cdots & a_n\\ | & | & & | \end{bmatrix} $$那么
$$ A^TA $$的第 $i,j$ 个元素是
$$ a_i^Ta_j $$也就是第 $i$ 列和第 $j$ 列的内积。
因此 $A^TA$ 记录了 $A$ 的列向量之间的相关关系。
它也叫 Gram 矩阵。
为什么 $A^TA$ 一定半正定
对任意向量 $z$,都有
$$ z^TA^TAz=(Az)^T(Az)=\|Az\|^2 $$而范数平方一定非负,所以
$$ z^TA^TAz\ge 0 $$因此:
$$ A^TA\text{ 一定半正定} $$这个结论不需要 $A$ 满秩,也不需要 $A$ 是方阵。
只要 $A$ 是实矩阵,$A^TA$ 就是对称半正定矩阵。
正定与半正定
设 $M\in\mathbb{R}^{n\times n}$ 是实对称矩阵。
如果对任意 $z$,都有
$$ z^TMz\ge 0 $$则称 $M$ 是半正定矩阵,记作
$$ M\succeq 0 $$如果对任意非零 $z\neq 0$,都有
$$ z^TMz>0 $$则称 $M$ 是正定矩阵,记作
$$ M\succ 0 $$区别在于:
$$ \text{半正定允许某些非零方向上等于 }0 $$$$ \text{正定要求所有非零方向上都严格大于 }0 $$因此:
$$ M\text{ 正定} \Longrightarrow M\text{ 半正定} $$但反过来不一定成立。
例如
$$ M= \begin{bmatrix} 1 & 0\\ 0 & 0 \end{bmatrix} $$对任意
$$ z= \begin{bmatrix} z_1\\z_2 \end{bmatrix} $$有
$$ z^TMz=z_1^2\ge 0 $$所以它半正定。
但取
$$ z= \begin{bmatrix} 0\\1 \end{bmatrix} \neq 0 $$有
$$ z^TMz=0 $$所以它不是正定。
因此:
$$ \text{正定是比半正定更强的条件} $$半正定、正定和特征值
对实对称矩阵 $M$,有以下等价关系:
$$ M\succeq 0 \Longleftrightarrow M\text{ 的所有特征值 }\lambda_i\ge 0 $$$$ M\succ 0 \Longleftrightarrow M\text{ 的所有特征值 }\lambda_i>0 $$这两个结论要求 $M$ 是实对称矩阵,或者复数情况下是 Hermitian 矩阵。
简要证明如下。
因为实对称矩阵可以正交对角化:
$$ M=Q\Lambda Q^T $$其中
$$ Q^TQ=I $$$$ \Lambda=\operatorname{diag}(\lambda_1,\dots,\lambda_n) $$于是
$$ z^TMz=z^TQ\Lambda Q^Tz $$令
$$ y=Q^Tz $$由于 $Q$ 是正交矩阵,$z\neq 0$ 等价于 $y\neq 0$。
于是
$$ z^TMz=y^T\Lambda y $$也就是
$$ z^TMz=\sum_{i=1}^n \lambda_i y_i^2 $$如果所有
$$ \lambda_i\ge 0 $$则显然
$$ z^TMz\ge 0 $$所以 $M$ 半正定。
如果所有
$$ \lambda_i>0 $$并且 $z\neq 0$,则 $y\neq 0$,至少有一个 $y_i\neq 0$,因此
$$ \sum_{i=1}^n \lambda_i y_i^2>0 $$所以 $M$ 正定。
反过来,如果 $M$ 半正定,那么取 $z$ 为某个特征向量 $v_i$,有
$$ v_i^TMv_i=v_i^T\lambda_i v_i=\lambda_i\|v_i\|^2\ge 0 $$由于
$$ \|v_i\|^2>0 $$所以
$$ \lambda_i\ge 0 $$正定情形同理可得
$$ \lambda_i>0 $$为什么正定矩阵一定可逆
设 $M$ 是正定矩阵。
正定意味着:
$$ \forall z\neq 0,\quad z^TMz>0 $$现在反设 $M$ 不可逆。
如果 $M$ 不可逆,那么存在非零向量 $z\neq 0$,使得
$$ Mz=0 $$于是
$$ z^TMz=z^T0=0 $$这和正定性要求的
$$ z^TMz>0 $$矛盾。
所以 $M$ 不可能不可逆。
因此:
$$ M\succ 0 \Longrightarrow M\text{ 可逆} $$从特征值角度看也一样。
正定矩阵的所有特征值都满足
$$ \lambda_i>0 $$所以没有零特征值。
没有零特征值,就可逆。
因此:
$$ M\succ 0 \Longrightarrow 0\text{ 不是特征值} \Longrightarrow M\text{ 可逆} $$$A$ 满列秩与 $A^TA$ 正定等价
这是最重要的桥梁之一。
结论:
$$ A\text{ 满列秩} \Longleftrightarrow A^TA\text{ 正定} $$证明如下。
先证明
$$ A\text{ 满列秩} \Longrightarrow A^TA\text{ 正定} $$若 $A$ 满列秩,则
$$ \mathcal{N}(A)=\{0\} $$所以对任意非零 $z\neq 0$,都有
$$ Az\neq 0 $$因此
$$ z^TA^TAz=\|Az\|^2>0 $$所以
$$ A^TA\text{ 正定} $$再证明反方向:
$$ A^TA\text{ 正定} \Longrightarrow A\text{ 满列秩} $$若 $A^TA$ 正定,则对任意非零 $z\neq 0$,有
$$ z^TA^TAz>0 $$而
$$ z^TA^TAz=\|Az\|^2 $$所以
$$ \|Az\|^2>0 $$即
$$ Az\neq 0 $$这说明不存在非零 $z$,使得
$$ Az=0 $$因此
$$ \mathcal{N}(A)=\{0\} $$所以 $A$ 满列秩。
于是二者等价。
$A$ 不满列秩与 $A^TA$ 奇异等价
结论:
$$ A\text{ 不满列秩} \Longleftrightarrow A^TA\text{ 奇异} $$证明如下。
先证明:
$$ A\text{ 不满列秩} \Longrightarrow A^TA\text{ 奇异} $$如果 $A$ 不满列秩,则存在
$$ z\neq 0 $$使得
$$ Az=0 $$于是
$$ A^TAz=A^T0=0 $$所以存在非零向量 $z$,使得
$$ A^TAz=0 $$因此
$$ A^TA $$不可逆,也就是奇异。
再证明反方向:
$$ A^TA\text{ 奇异} \Longrightarrow A\text{ 不满列秩} $$如果 $A^TA$ 奇异,则存在
$$ z\neq 0 $$使得
$$ A^TAz=0 $$左乘 $z^T$,得到
$$ z^TA^TAz=0 $$也就是
$$ \|Az\|^2=0 $$所以
$$ Az=0 $$于是存在非零 $z$,使得
$$ Az=0 $$因此
$$ A\text{ 不满列秩} $$所以二者等价。
相关性质的充分、必要和条件范围
很多矩阵结论必须注意条件范围。
例如:
$$ A\text{ 满列秩} \Longleftrightarrow A^TA\text{ 可逆} $$这个结论对任意
$$ A\in\mathbb{R}^{m\times n} $$都成立。
但是:
$$ A\text{ 可逆} $$这个说法只在 $A$ 是方阵时才有意义。
如果 $A$ 是 $m\times n$ 的长方矩阵,那么不能说 $A$ 可逆或奇异。只能说它列满秩、行满秩、满秩、零空间是否平凡。
又如:
$$ M\succ 0 \Longrightarrow M\text{ 可逆} $$但
$$ M\text{ 可逆} \not\Longrightarrow M\succ 0 $$例如
$$ M= \begin{bmatrix} 1 & 0\\ 0 & -1 \end{bmatrix} $$它可逆,因为特征值为 $1$ 和 $-1$,都不为零。
但它不是正定,因为取
$$ z= \begin{bmatrix} 0\\1 \end{bmatrix} $$有
$$ z^TMz=-1 < 0 $$所以可逆不等于正定。
再如:
$$ M\succ 0 \Longrightarrow M\succeq 0 $$但
$$ M\succeq 0 \not\Longrightarrow M\succ 0 $$因为半正定允许零特征值,而正定不允许零特征值。
Moore-Penrose 伪逆
普通逆矩阵只对可逆方阵存在。
但在很多问题中,矩阵可能不是方阵,或者虽然是方阵但不可逆。此时仍然希望有一个“类似逆矩阵”的对象。
这就是 Moore-Penrose 伪逆,记为
$$ A^+ $$它对任意矩阵都存在。
如果 $A$ 是可逆方阵,那么
$$ A^+=A^{-1} $$如果 $A$ 不可逆或者不是方阵,$A^+$ 仍然可以定义。
最清楚的定义来自奇异值分解。
设
$$ A=U\Sigma V^T $$其中 $\Sigma$ 是奇异值矩阵:
$$ \Sigma= \operatorname{diag}(\sigma_1,\dots,\sigma_r,0,\dots) $$正奇异值满足
$$ \sigma_i>0 $$伪逆定义为
$$ A^+=V\Sigma^+U^T $$其中 $\Sigma^+$ 的做法是:
$$ \sigma_i\mapsto \frac{1}{\sigma_i} $$而零奇异值仍然保持为零。
也就是说,普通逆矩阵要求所有方向都能反解;伪逆只在 $A$ 没有丢失的方向上做反解,对已经被压成零的方向不强行反解。
伪逆和零空间的关系
如果 $A$ 不满列秩,则存在非零
$$ z\in\mathcal{N}(A) $$所以
$$ A(x+z)=Ax $$这说明很多不同的 $x$ 对应同一个输出。
此时如果我们想解
$$ Ax=b $$或者最小化
$$ \|Ax-b\|^2 $$解可能不唯一。
Moore-Penrose 伪逆给出的解
$$ x=A^+b $$具有一个重要性质:
它是所有最小二乘解中范数最小的那个。
如果有多个 $x$ 都能给出同样的最小残差,那么 $A^+b$ 会选择长度最短的那个。
从空间上说,它会去掉零空间中的任意漂移分量。
因为零空间分量不会改变 $Ax$,但会增大 $\|x\|$。所以最小范数解不会包含零空间分量。
普通最小二乘的解性质
普通最小二乘问题为:
$$ \min_x\|Ax-b\|^2 $$几何上,它是在列空间 $\mathcal{C}(A)$ 中寻找一个向量
$$ A\hat{x} $$使得它最接近 $b$。
也就是说,
$$ A\hat{x} $$是 $b$ 在列空间 $\mathcal{C}(A)$ 上的正交投影。
因此残差
$$ r=b-A\hat{x} $$必须垂直于列空间:
$$ r\perp \mathcal{C}(A) $$这等价于
$$ A^Tr=0 $$代入
$$ r=b-A\hat{x} $$得到
$$ A^T(b-A\hat{x})=0 $$也就是
$$ A^TA\hat{x}=A^Tb $$这就是法方程。
普通最小二乘在不同秩条件下的结果
$A$ 列满秩
如果
$$ \operatorname{rank}(A)=n $$那么
$$ A^TA\text{ 正定} $$所以
$$ A^TA\text{ 可逆} $$于是最小二乘解唯一:
$$ \hat{x}=(A^TA)^{-1}A^Tb $$这时不存在非零 $z$,使得
$$ Az=0 $$所以不存在“加一个零空间方向但输出不变”的自由度。
$A$ 不满列秩
如果
$$ \operatorname{rank}(A) < n $$那么存在
$$ z\neq 0 $$使得
$$ Az=0 $$因此如果 $\hat{x}$ 是一个最小二乘解,那么
$$ \hat{x}+z $$也是一个最小二乘解,因为
$$ A(\hat{x}+z)=A\hat{x} $$这时解不唯一。
同时
$$ A^TA $$奇异,所以不能写
$$ (A^TA)^{-1}A^Tb $$但可以写
$$ \hat{x}=A^+b $$其中 $A^+b$ 是最小范数最小二乘解。
更一般地,所有最小二乘解可表示为
$$ x=A^+b+(I-A^+A)w $$其中 $w$ 是任意向量。
这里
$$ (I-A^+A)w $$属于 $A$ 的零空间部分。它不会改变 $Ax$,只是在不同的等价解之间移动。
一个具体例子串起来
回到开头的矩阵:
$$ A= \begin{bmatrix} 1 & 2\\ 2 & 4\\ 3 & 6 \end{bmatrix} $$它的第二列是第一列的两倍,所以
$$ \operatorname{rank}(A)=1 < 2 $$因此 $A$ 不满列秩。
我们可以找到非零向量
$$ z= \begin{bmatrix} 2\\ -1 \end{bmatrix} $$使得
$$ Az=0 $$因为
$$ 2a_1-a_2=0 $$这说明 $A$ 的零空间非平凡:
$$ \mathcal{N}(A)\neq \{0\} $$于是 $A^TA$ 必然奇异。
实际计算:
$$ A^TA= \begin{bmatrix} 1 & 2 & 3\ 2 & 4 & 6 \end{bmatrix} \begin{bmatrix} 1 & 2\ 2 & 4\ 3 & 6 \end{bmatrix}
\begin{bmatrix} 14 & 28\ 28 & 56 \end{bmatrix} $$
第二列仍然是第一列的两倍:
$$ \begin{bmatrix} 28\56 \end{bmatrix}
2 \begin{bmatrix} 14\28 \end{bmatrix} $$
所以
$$ A^TA $$也不满秩,因而不可逆。
同时
$$ z^TA^TAz=\|Az\|^2=0 $$其中 $z\neq 0$,所以 $A^TA$ 不是正定。
但对任意 $u$,仍然有
$$ u^TA^TAu=\|Au\|^2\ge 0 $$所以 $A^TA$ 是半正定。
因此它的特征值满足:
$$ \lambda_i\ge 0 $$但至少有一个特征值等于零。
实际上,
$$ A^TA= \begin{bmatrix} 14 & 28\\ 28 & 56 \end{bmatrix} $$的特征值是
$$ 70,\quad 0 $$这正好对应:
$$ A^TA\text{ 半正定但不正定} $$$$ A^TA\text{ 有零特征值} $$$$ A^TA\text{ 奇异} $$$$ A^TA\text{ 不可逆} $$$$ A\text{ 不满列秩} $$加上 Tikhonov 正则化之后
虽然本文主要讲矩阵性质,但可以顺手看一下 Tikhonov 正则化为什么有效。
普通最小二乘的法方程是:
$$ A^TAx=A^Tb $$如果 $A$ 不满列秩,则 $A^TA$ 不可逆。
Tikhonov 正则化改为:
$$ \min_x\|Ax-b\|^2+\lambda\|x\|^2 $$求导得到:
$$ (A^TA+\lambda I)x=A^Tb $$其中
$$ \lambda>0 $$关键是:
$$ A^TA+\lambda I $$一定正定。
因为对任意非零 $z\neq 0$:
$$ zT(ATA+\lambda I)z
zTATAz+\lambda z^Tz $$
也就是
$$ zT(ATA+\lambda I)z
|Az|^2+\lambda|z|^2 $$
其中
$$ \|Az\|^2\ge 0 $$而
$$ \lambda\|z\|^2>0 $$所以
$$ \|Az\|^2+\lambda\|z\|^2>0 $$因此
$$ A^TA+\lambda I $$正定。
正定矩阵一定可逆,所以
$$ (A^TA+\lambda I)^{-1} $$一定存在。
这就是 Tikhonov 正则化能把病态或不唯一问题变成稳定唯一问题的原因。
各性质关系图
下面这张图可以作为总览。
flowchart TD
A["A ∈ R^{m×n}"] --> B["列向量 a1,...,an"]
B --> C["列空间 C(A)"]
B --> D["秩 rank(A)=dim C(A)"]
D --> E{"rank(A)=n ?"}
E -->|是| F["A 列满秩"]
F --> G["N(A)={0}"]
G --> H["不存在 z≠0 使 Az=0"]
H --> I["A^T A 正定"]
I --> J["A^T A 所有特征值 > 0"]
J --> K["A^T A 非奇异"]
K --> L["A^T A 可逆"]
L --> M["OLS 解唯一: x=(A^T A)^(-1)A^T b"]
E -->|否| N["A 不满列秩"]
N --> O["N(A) 非平凡"]
O --> P["存在 z≠0 使 Az=0"]
P --> Q["A^T A z=0"]
Q --> R["A^T A 有零特征值"]
R --> S["A^T A 奇异"]
S --> T["A^T A 不可逆"]
T --> U["OLS 不能用普通逆公式"]
U --> V["用 Moore-Penrose 伪逆 A^+b"]
方阵性质关系图
对于方阵
$$ M\in\mathbb{R}^{n\times n} $$有下面的等价关系:
flowchart LR
A["M 可逆"] <--> B["M 非奇异"]
B <--> C["rank(M)=n"]
C <--> D["N(M)={0}"]
D <--> E["Mx=0 只有零解"]
E <--> F["det(M)≠0"]
F <--> G["0 不是 M 的特征值"]
对应地:
flowchart LR
A["M 不可逆"] <--> B["M 奇异"]
B <--> C["rank(M)<n"]
C <--> D["N(M) 非平凡"]
D <--> E["存在 x≠0 使 Mx=0"]
E <--> F["det(M)=0"]
F <--> G["0 是 M 的特征值"]
正定、半正定、特征值关系图
对实对称矩阵 $M$:
flowchart TD
A["M 正定 M ≻ 0"] --> B["∀z≠0, z^T M z > 0"]
B --> C["所有特征值 λ_i > 0"]
C --> D["无零特征值"]
D --> E["M 可逆"]
F["M 半正定 M ⪰ 0"] --> G["∀z, z^T M z ≥ 0"]
G --> H["所有特征值 λ_i ≥ 0"]
A --> F
C --> H
注意这里有几个非等价关系:
$$ M\succ 0 \Longrightarrow M\succeq 0 $$但
$$ M\succeq 0 \not\Longrightarrow M\succ 0 $$$$ M\succ 0 \Longrightarrow M\text{ 可逆} $$但
$$ M\text{ 可逆}\not\Longrightarrow M\succ 0 $$$A$、$A^TA$、最小二乘之间的关系表
| 条件 | 等价和推出关系 | 结果 |
|---|---|---|
| $A$ 列满秩 | $\operatorname{rank}(A)=n$ | 列向量线性无关 |
| $A$ 列满秩 | $\mathcal{N}(A)=\{0\}$ | 没有非零输入被压成零 |
| $A$ 列满秩 | $A^TA\succ 0$ | $A^TA$ 正定 |
| $A$ 列满秩 | $A^TA$ 可逆 | OLS 普通逆公式成立 |
| $A$ 不满列秩 | $\operatorname{rank}(A) < n$ | 列向量线性相关 |
| $A$ 不满列秩 | $\exists z\neq 0,\ Az=0$ | 存在零空间自由度 |
| $A$ 不满列秩 | $A^TA$ 奇异 | 不能直接求 $(A^TA)^{-1}$ |
| $A$ 不满列秩 | $A^TA\succeq 0$ 但非正定 | 有零特征值 |
| $A^TA$ 奇异 | $A^TA$ 不可逆 | OLS 解可能不唯一 |
| OLS 解不唯一 | 存在零空间自由度 | 可用 $A^+b$ 选最小范数解 |
总结
矩阵 $A$ 的列向量决定了它能表示哪些输出。列向量独立,说明每个参数方向都有独立作用;列向量相关,说明存在某些参数方向对输出没有影响。
如果
$$ A\text{ 不满列秩} $$那么存在非零方向
$$ z\neq 0 $$使得
$$ Az=0 $$这个方向就是零空间方向。沿着这个方向移动,输出 $Ax$ 不变。
因此在最小二乘中,多个不同的 $x$ 可能产生同一个 $Ax$,解就可能不唯一。
同时,由于
$$ z^TA^TAz=\|Az\|^2 $$所以 $A^TA$ 一定半正定。
如果 $A$ 满列秩,则非零 $z$ 不可能被 $A$ 压成零,所以
$$ z^TA^TAz>0 $$于是 $A^TA$ 正定、可逆,普通最小二乘解唯一。
如果 $A$ 不满列秩,则存在非零 $z$,使得
$$ Az=0 $$于是
$$ z^TA^TAz=0 $$所以 $A^TA$ 只是半正定,不是正定;它有零特征值,因此奇异、不可逆。
这时普通逆公式失效,但 Moore-Penrose 伪逆仍然可以给出最小范数最小二乘解。
而 Tikhonov 正则化通过加入
$$ \lambda I $$把
$$ A^TA $$变成
$$ A^TA+\lambda I $$从而把半正定提升为正定,保证可逆性和稳定性。
总结:
$$ \text{列相关} \Rightarrow \text{零空间非平凡} \Rightarrow \text{信息丢失} \Rightarrow A^TA\text{ 奇异} \Rightarrow \text{普通逆失效} \Rightarrow \text{伪逆或正则化} $$