跳到正文
格致开物MATHWIKI

矩阵平方根

矩阵平方根(matrix square root)是满足 X2=A 的方阵 X。矩阵乘法中的交叉项使这一问题不同于对每个元素开平方;解可能有多个,也可能不存在。对于 Hermitian 半正定矩阵,则有唯一的 Hermitian 半正定平方根,通常记作 A1/2。

取实对称矩阵 A=(5445)。候选 S=(2112) 满足 S2=A。它的元素并不是 A 各元素的标量平方根;正确性须通过完整的矩阵乘法检查。

谱构造与完整算例

Hermitian 半正定矩阵有谱分解 A=Udiag⁡(λ1,…,λn)U∗,其中 λj≥0。定义 S=Udiag⁡(λ1,…,λn)U∗. 平方时中间的 U∗U 消去,立刻得到 S2=A,并且谱非负,故 S 半正定。

主例的单位特征方向是 u+=(1,1)𝖳/2 和 u−=(1,−1)𝖳/2,特征值分别为九和一。因此 S=3u+u+𝖳+u−u−𝖳=12(3333)+12(1−1−11)=(2112). 对向量 u+,一次作用 S 放大三倍,两次便是九倍,恰与 A 相同。这给出了“平方根”在线性变换中的直接含义。

半正定根的唯一性

若 R 也 Hermitian 半正定且 R2=A,则 RA=R3=AR,所以 R 保持 A 的每个特征空间。在 A 的 λ 特征空间上,限制算子仍 Hermitian 半正定,可用谱定理对角化,其每个特征值 r 都满足 r2=λ,r≥0,故只能为 λ。于是它在整个该空间上等于 λI,与 S 一致。这同时覆盖零特征值。

唯一性须保留“半正定”限制。主例还存在平方根 −S,以及 T=(1221);后者的特征值为三和负一,平方仍给九和一,所以 T2=A,却不是半正定根。

奇异矩阵可能没有根

令 N=(0100)。假设 X2=N,则 XN=X3=NX。将一般二阶矩阵代入交换关系可得 X=(ab0a)。其平方为 X2=(a22ab0a2). 若它等于 N,对角元要求 a=0,右上角却因此变零,与一矛盾。故这个实矩阵甚至没有复矩阵平方根。不能把“复数都有平方根”推广成“复方阵都有平方根”。

若矩阵的谱避开非正实轴,可定义唯一谱位于右半平面的主平方根;这是比任意解更具体的分支选择。对奇异半正定矩阵,前述唯一半正定根仍存在,不要把两种定理的不同条件混为一谈。[1]

变式取 A=diag⁡(0,16),半正定根为 diag⁡(0,4)。它虽不可逆,却满足唯一性定理;问题不是“奇异就没有根”,而是奇异矩阵的结构需要逐一判断。

参考文献

Sheldon Axler,Linear Algebra Done Right,第四版,§7C、8C(正算子平方根与一般平方根);本文例子与图解独立构造。

  1. ↑ Nicholas J. Higham,What Is a Matrix Square Root?,关于主平方根及不存在、非唯一平方根的实例。