[论文解读] The asymptotics of Wilkinson's shift iteration

Ricardo S. Leite, Nicolau C. Saldanha+1 more|ArXiv.org|Dec 27, 2004

Advanced Topics in Algebra参考文献 15被引用 5

一句话总结

本文分析了雅可比矩阵上威尔金森位移QR迭代的收敛速率，证明当谱为无三元等差数列（AP-free）时，收敛为三阶。对于某些具有等差数列谱的3×3矩阵，收敛严格为二阶，且收敛至不动点的初始条件集合构成一个符号序列的康托集，其动力学由符号动力学和双对角坐标控制。

ABSTRACT

We study the rate of convergence of Wilkinson's shift iteration acting on Jacobi matrices with simple spectrum. We show that for AP-free spectra (i.e., simple spectra containing no arithmetic progression with 3 terms), convergence is cubic. In order 3, there exists a tridiagonal symmetric matrix P_0 which is the limit of a sequence of a Wilkinson iteration, with the additional property that all iterations converging to P_0 are strictly quadratic. Among tridiagonal matrices near P_0, the set X of initial conditions with convergence to P_0 is rather thin: it is a union of disjoint arcs X_s meeting at P_0, where s ranges over the Cantor set of sign sequences s: N -> {1,-1}. Wilkinson's step takes X_s to X_{s'}, where s' is the left shift of s. Among tridiagonal matrices conjugate to P_0, initial conditions near P_0 but not in X converge at a cubic rate.

研究动机与目标

• 确定威尔金森位移QR迭代在实对称三对角矩阵（具有单谱）上的渐近收敛速率。
• 解决长期存在的猜想：对于通用矩阵，收敛是否为三阶，特别是当谱中不存在等差数列时。
• 刻画收敛至不动点的初始条件结构，尤其在3×3情形下收敛可能为二阶时。
• 应用符号动力学与双对角坐标分析迭代在奇点和不连续点附近的性质。
• 描述收敛至特定不动点的初始条件集合的分形结构，表明其为由符号序列索引的康托集。

提出的方法

• 使用双对角坐标（特征值与次对角参数）将全局收敛问题转化为光滑坐标图中的局部分析。
• 将威尔金森位移定义为最下角2×2子矩阵中距离(n,n)元素最近的特征值，以确保QR步保持对称三对角。
• 在不动点附近对迭代映射进行泰勒展开，推导三阶收敛估计，尤其当位移函数在最后一个次对角参数中为奇函数时。
• 分析在集合Y（特征值等距）和Z（位移等于特征值）处的不连续性，表明在特定条件下迭代可光滑延拓过Z。
• 使用符号动力学建模迭代对收敛行为符号序列的作用，将每个初始条件映射至康托符号序列集中的一个序列。
• 建立收敛至不动点p₀的初始条件集合X为由符号序列s索引的互不相交弧Xₛ的并集，且映射W将Xₛ映射至Xₛ′，其中s′为s的左移。

实验结果

研究问题

• RQ1 威尔金森位移QR迭代对所有具有单谱的雅可比矩阵是否都呈三阶收敛？
• RQ2 在3×3情形下，收敛至不动点的初始条件集合结构如何，特别是当收敛为二阶时？
• RQ3 当位移函数在某些点出现不连续或奇点时，迭代行为如何？
• RQ4 是否可使用符号动力学描述收敛动力学？收敛初始条件集合的拓扑性质为何？
• RQ5 收敛速率与谱结构之间存在何种关系，特别是特征值中是否存在等差数列？

主要发现

• 对于无三元等差数列谱（AP-free）的情形，威尔金森位移QR迭代的收敛为三阶：|(T_{k+1})_{n,n-1}| ≤ C|(T_k)_{n,n-1}|^3，其中C > 0为某常数。
• 在3×3情形下若谱中存在等差数列，则存在一个不动点p₀，使得收敛至p₀为严格二阶，且收敛至p₀的初始条件集合为符号序列的康托集。
• 收敛至p₀的初始条件集合X为由符号序列s ∈ {1, -1}^ℕ索引的互不相交弧Xₛ的并集，在自然拓扑下构成康托集。
• 迭代映射W将Xₛ映射至Xₛ′，其中s′为s的左移，从而在收敛集合上建立符号动力学结构。
• 对于靠近p₀但不在X中的初始矩阵，收敛严格为三阶，表明二阶收敛被限制在薄的、分形状集合中。
• 集合X的豪斯多夫维数在p₀附近至少为1，数值计算表明其为光滑曲线Xₛ的并集，由利普希茨函数参数化。