下一节:GEMV形式的Lanczos方法 上一级:隐式重启Lanczos 上一节:隐式重启

位移的选择

选择应用于QR步骤的位移 \{\mu_j\} 有许多方法。几乎任何显式的多项式重启方案都可以通过这种隐式机制来应用。选择精确位移取得了相当大的成功。这种选择是通过将 T_m 的特征值分成两个不相交的集合,即 k 个“想要的”和 p 个“不想要的”特征值,并使用这 p 个不想要的特征值作为位移。通过这种选择,p 次位移应用后,T_k^+ 的特征值谱就包含了这 k 个想要的特征值。随着收敛的进行,T_k 的次对角线元素趋向于零,而最期望的特征值近似值作为 T_m 的前 k \times k 块的对角线上的特征值出现。基向量 V_k 趋向于成为 A 的正交特征向量。

“想要的集合”的指定示例包括:

代数上最小的 k 个特征值,
代数上最大的 k 个特征值,
模最大的 k 个特征值,
模最小的 k 个特征值。
在选择这些选项时,请记住,谱两端的特征值首先收敛,因此任何避免这些特征值的选择都可能导致收敛缓慢。当寻求内部特征值时,建议使用SI(位移-逆)。

其他有趣的位移选择 \{\mu_j\} 包括切比雪夫多项式的根[383]、调和Ritz值[331,337,349,411]、Leja多项式的根[23]、最小二乘多项式的根[384]以及精炼位移[244]。特别是,Leja多项式和调和Ritz值已被用于估计 A 的内部特征值。

另一种解释源于以下事实:每一次位移循环都会隐式地将一个 A 的多项式(次数为 p)应用于初始向量:

v_1 \leftarrow \psi(A) v_1, \quad \text{其中} \quad \psi(\lambda) = \prod_{j=1}^p (\lambda - \mu_j).
这个多项式的根就是在QR过程中使用的位移,这些位移可以选择来增强初始向量在期望特征值对应的特征向量方向上的分量,并抑制不希望方向上的分量。当然,这是可取的,因为它迫使初始向量进入与期望特征值相关的不变子空间。这反过来又使得 r_k 变小,从而实现收敛。详细信息可以在[419]中找到。

还有一种隐式重启的方法,即厚重启,由Wu和Simon[463]描述。在那里,计算了 T_m 的特征值分解(4.20),并保留对应于想要的特征值的部分作为箭头矩阵,其余部分则被丢弃。如果选择了相同的想要的和不想要的特征值,厚重启在数学上等价于精确位移IRLM。



下一节:GEMV形式的Lanczos方法 上一级:隐式重启Lanczos 上一节:隐式重启
Susan Blackford 2000-11-20