“用复 Schur 分解写 $A=QTQ^$、$B=USU^$。置 $Y=Q^ XU$、$D=Q^ CU$,得到”
形式陈述
对任意复方阵
对实矩阵
其中
若
浮点算法以 Schur 形而非 Jordan 形作为一般全谱目标。正交相似变换不放大
量级的后向稳定结果,常数依矩阵规模与实现而定。这里稳定的是“邻近矩阵具有所算 Schur 形”这一关系;当特征值聚集时,具体 Schur 向量和块排序仍可能对扰动敏感。
直觉
对角化试图找到一组彼此独立的特征轴,一般矩阵却可能没有足够特征向量,或者特征向量基极度歪斜。Schur 分解退一步:只要求一组正交坐标把矩阵变成上三角。对角线交付全部特征值,上三角非对角部分保留方向之间尚未解耦的作用,而坐标变换始终保持长度与夹角。
这一步退让换来稳定性。Jordan 形要判断某些量究竟是精确零还是极小非零,微小扰动就可能改变块尺寸;Schur 形不试图做这项不连续判决,而是把全部耦合诚实保留在上三角部分。
例子与边界
剪切矩阵
本身已经是复 Schur 形,可以取
实旋转
没有实特征向量,实 Schur 形只能保留这个
近缺陷族
在
推论与应用
Jordan 标准形精确编码广义特征向量链,适合代数分类与矩阵函数公式;Schur 分解不提供 Jordan 块尺寸,却是一般浮点矩阵的稳定计算目标。二者回答不同问题,不能用“Schur 近似 Jordan”混成一条主线。
QR 特征值算法先把矩阵化为 Hessenberg 形,再通过隐式移位与 deflation 收敛到 Schur 形。这实现了本页分解的计算接口:返回
Sylvester 方程的谱不交判据说明何时可以逐块求解,Bartels–Stewart 算法给出两侧 Schur 变换、右端变换及块回代的完整执行。Schur–Parlett 算法则先求对角块上的函数,再用 Sylvester 方程恢复跨块耦合;接近特征值需要一起处理,不能直接除以很小的对角差。
参考资料
- Lloyd N. Trefethen and David Bau III, Numerical Linear Algebra, SIAM, 1997, Lectures 24 and 25.
- Roger A. Horn and Charles R. Johnson, Matrix Analysis, 2nd ed., Cambridge University Press, 2013, Ch. 2.
- LAPACK Users’ Guide, 3rd ed., SIAM, 1999, Schur Factorization.