下一节:特征问题的具体化 上一级:厄米特征问题 上一节:特征分解

良态与病态

矩阵 A 的特征值总是处于良态, 这意味着在范数上将 A 改变最多 \epsilon 时, 任何特征值的变化也不会超过 \epsilon。 我们参考第4.8节中的技术定义。

这对于大多数目的来说已经足够了,除非用户对小特征值的前几位有效数字感兴趣, 即特征值的大小不超过 \epsilon。 例如,计算 \lambda_i = 10^{-5} 的误差范围在正负 \epsilon = 10^{-4} 内意味着 计算出的 \lambda_i 的前几位有效数字可能都不正确。 参见[114,118]关于小特征值敏感性及其前几位有效数字何时能准确计算的讨论。

另一方面,特征向量和特征空间可能是病态的。 例如,做如下改变

A_0 = \begin{bmatrix} 2 & 0 & 0 \\ 0 & 1 & \epsilon \\ 0 & \epsilon & 1 \end{bmatrix} \quad \to \quad A_1 = \begin{bmatrix} 2 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1+\epsilon \end{bmatrix}
会使对应于接近1的两个特征值的两个特征向量旋转\pi/4,无论\epsilon有多小。因此它们对微小变化非常敏感。 特征向量的条件数取决于其特征值与最近其他特征值之间的间隔:间隔越小,特征向量越敏感。 在这个例子中,接近1的两个特征值非常接近,因此它们的间隔很小,对应的特征向量也很敏感。 但它们所张成的二维不变子空间对A的变化非常不敏感(因为它们的特征值都接近1,与下一个最近的特征值2相距甚远)。 因此,当对应于一组接近特征值的特征向量过于病态时,用户可能希望计算它们所张成的不变子空间的基,而不是单个特征向量。



下一节:特征问题的具体化 上一级:厄米特征问题 上一节:特征分解
Susan Blackford 2000-11-20