跳到正文
格致开物MATHWIKI

半正定矩阵

半正定矩阵(positive semidefinite matrix)在实数情形通常指满足 x𝖳Ax≥0对每个 x∈ℝn 的实对称矩阵 A,记为 A⪰0。复数情形要求 A=A∗,并以 x∗Ax≥0 定义。若对每个非零 x 不等式都严格成立,就称正定。半正定允许非零方向的二次型为零。[1]

例如 A=(1111) 给出 x𝖳Ax=(x1+x2)2≥0,所以半正定;但取 x=(1,−1)𝖳 得到零,因此不正定。这里检查的是二次型对所有向量的值,而不是要求矩阵每个元素都非负。

特征值与平方范数的等价刻画

对实对称矩阵,以下条件等价:二次型处处非负;全部特征值非负;存在某个实矩阵 B 使 A=B𝖳B。

由谱定理可写 A=QΛQ𝖳,其中 Q 正交,Λ=diag⁡(λ1,…,λn)。令 y=Q𝖳x,得到 x𝖳Ax=∑i=1nλiyi2. 若所有特征值非负,和式就非负;若某个 λj<0,取 x=Qej 就得到负值,违反半正定性。这证明前两者等价。

当 λi≥0 时,令 B=Λ1/2Q𝖳,直接相乘有 B𝖳B=A。反向若有这种因子,则 x𝖳Ax=‖Bx‖22≥0。这些步骤也给出正半定平方根 A1/2=QΛ1/2Q𝖳。

原例特征方向为 (1,1) 与 (1,−1),特征值分别为二和零。沿前一方向有正的二次型,沿后一方向完全平坦。

平面中x一加x二为零的直线是零方向,两侧平行等值线分别对应平方值一和四。
二次型只取决于两个坐标的和;沿零方向移动不会改变它。

图中的每条等值线都平行于核空间方向。固定 x1+x2=c 后,二次型为 c2;因此两条标记为一的线对应 c=±1,标记为四的对应 c=±2。图没有把这些无界直线误画成正定矩阵的封闭椭圆。

主子式检验与边界反例

实对称矩阵半正定,当且仅当所有主子式非负。主子矩阵要求同时选取相同索引的行与列;只检查左上角逐阶的顺序主子式不足以判定半正定。例如 diag⁡(0,−1) 的左上角一阶子式与全行列式均为零,但在 x=(0,1) 处二次型为负。[2]

二阶对称矩阵 (accb) 的全部主子式恰为 a,b,ab−c2,所以条件为三者均非负。变式 (2−2−22) 虽有负元素,但二次型是 2(x1−x2)2,仍半正定。相反,(1221) 所有元素都正,却在 (1,−1) 处取负二。

仅检查特征值非负也要保留对称条件。矩阵 (1401) 的两个特征值都是一,但在 x=(1,−1) 处二次型为负二;它本来就不对称,不能使用上面的谱定理刻画。若只关心一般实矩阵的二次型,真正起作用的是对称部分 (A+A𝖳)/2。

半正定锥与矩阵比较

若 A,B⪰0 且 α,β≥0,则对每个 x 有 x𝖳(αA+βB)x≥0,故半正定矩阵构成凸锥。定义 A⪰B 表示 A−B⪰0,称为半正定序。它比较的是所有方向上的二次型,不是逐元素比较。

例如 A=(2−1−12) 满足 A⪰I,因为 A−I 的二次型为 (x1−x2)2。虽然 A 的非对角元小于 I 的对应元素,矩阵序关系仍成立。

对任意实矩阵 C,合同变换 C𝖳AC⪰0,因为其二次型是 (Cx)𝖳A(Cx)。若 C 可逆,反向也成立。普通乘积则不同:两个半正定矩阵的乘积一般不对称,因而不一定属于半正定锥;额外的可交换条件会改变这一结论。

半正定性也自然出现在协方差矩阵中。若随机向量 X 有有限二阶矩,协方差矩阵 Σ 满足 u𝖳Σu=Var⁡(u𝖳X)≥0。零方向表示某个线性组合无随机波动;“半正定”保留了这种完全线性依赖,而正定性将它排除。

参考文献

  1. ↑ Stephen Boyd、Lieven Vandenberghe,Convex Optimization,附录A.5及2.2.5节。
  2. ↑ Nathan Barczi,14.102 Math for Economists: Lecture Notes, 10/7/2004,MIT,定理176。