舒尔分解
舒尔分解(Schur decomposition)把任意复方阵写成 ,其中 为酉矩阵, 为上三角矩阵。它不要求矩阵可对角化; 的对角元是特征值,非对角项保留特征方向之间未能消去的耦合。
对于数值计算,酉矩阵换基保持欧氏长度,因而不会仅因坐标变换本身放大误差。这是舒尔分解相对于可能高度病态的特征向量基的一项优势,但不保证原矩阵的所有特征值问题都条件良好。
复数域上的存在性证明
对阶数归纳。一阶矩阵已上三角。设 ,由代数基本定理选特征值 及单位特征向量 ,将其扩充成标准正交基,组成 。首列 使 对 阶矩阵 使用归纳假设,写 。再令 ,得到所需上三角分解。每一步都选正交补足,不要求先找到 个特征向量。
由于相似变换保留特征多项式,且上三角行列式是对角元之积,。这证明对角元按代数重数给出全部特征值。
一个不可对角化的算例
取 。其特征多项式为 ,单位特征向量可取 ,再取正交的 。 可以逐列核对:,而 。后式解释了右上角的二; 并非另一个特征向量。 的核只有一维,所以不存在对角化,却完全不影响这个舒尔分解。
正规矩阵时进一步变为对角
若 正规,酉相似后 仍正规。比较 与 左上角,得 ,迫使第一行非对角元全零;在余下块中重复, 即为对角矩阵。反过来,酉相似于对角矩阵必正规。所以舒尔分解说明了复谱定理的边界:一般矩阵只能保证上三角,正规性才保证进一步消去全部耦合。
实方阵不一定存在实上三角的正交舒尔分解,因为实数域可能没有特征值。实舒尔形式改用一阶与二阶实对角块的准上三角矩阵;例如旋转九十度矩阵的非实共轭特征值保留在一个二阶块中。本文的存在性证明针对复数域,不把它当成实上三角化的证明。
结果的解释与变式
主例 的对角元都为二,但 。只保留对角元会遗漏第二项,也就遗漏瞬时耦合对幂和指数的影响。
变式把 改成 ,用同样两列组成的 ,得到 。这次矩阵实对称,三角形中的非对角项确实可以全零。两例的区别由结构定理解释,而非某次计算偶然找到更漂亮的坐标。
参考文献
Sheldon Axler,Linear Algebra Done Right,第四版,§6B(Schur定理)、7B(谱定理);本文例子与图解独立构造。