矩阵对数
矩阵对数(matrix logarithm)是矩阵指数的逆问题:给定方阵 ,寻找满足 的矩阵 。它通常不唯一,而且“有复对数”“有实对数”“存在所选分支的主对数”是不同命题。逐元素取标量对数一般不解决这个问题。
例如 ,。令 由于 与 交换,且幂零级数截断,。右上角是二分之一,而不是对数一所得的零。
可逆性与分支
矩阵指数总可逆,因为 。所以奇异矩阵不可能有对数。反过来,每个可逆复矩阵都有复对数:在其约旦分解中,写单块 ,其中 , 幂零。选择一个标量值 使 ,并置 在商代数 中,标量形式幂级数的恒等式 截断后依然成立,故 。分别处理各块并相似变换回去,即构造出一个复对数。
如果 是对数,则 也是(整数 ),已经说明不唯一。重复特征值还可能允许更多不同结构的对数。
主对数与正定情形
若 的谱不与非正实轴 相交,则存在唯一的对数,其特征值的虚部都严格位于 ,称为主对数。这里必须排除零特征值;主分支的条件比“矩阵可逆”更强。[1]
对于 Hermitian 正定矩阵 ,所有 ,主对数为 ,仍 Hermitian。它是唯一的 Hermitian 对数:任何 Hermitian 可酉对角化,实指数在实轴上一一对应,故其谱与特征子空间上的作用都由 确定。
实对数并非总存在
一阶实矩阵 没有实对数,因为实指数恒正。二阶的 却有:令 ,则 。同一个负特征值在不同维数及块结构下可有不同结果,不能把标量结论原样推广。
对接近单位矩阵的情形,若某次乘性范数满足 ,主对数可按绝对收敛级数计算 直接截断的误差可用尾项范数和控制;当矩阵不满足这一范围时,不能只因每个元素较小就宣称级数收敛。
变式取 ,虽然 不满足上述充分条件,仍有 ,所以有限式 完全正确。范数小于一是保证一般无限级数的方便条件,不是所有对数存在的必要条件。
参考文献
Sheldon Axler,Linear Algebra Done Right,第四版,§8C(约旦结构);Higham矩阵对数作者说明;本文例子与图解独立构造。
- ↑ Nicholas J. Higham,What Is the Matrix Logarithm?,关于主对数、实对数与幂级数的说明。