[论文解读] New Recurrence Relationships between Orthogonal Polynomials which Lead to New Lanczos-type Algorithms

Muhammad Farooq, Abdellah Salhi|arXiv (Cornell University)|Mar 3, 2014

Matrix Theory and Algorithms参考文献 33被引用 6

一句话总结

本文提出了一组关于形式正交多项式（FOPs）的新递推关系，以开发求解线性系统的新型Lanczos型算法。通过突破传统次数差约束的限制，该研究识别出此前未被探索的递推结构——具体为 $A_{12}, A_{13}, A_{16}, A_{19}, B_{13}, B_{15}, B_{16}$——这些结构在数学上有效且适用于算法实现，从而可构建新型Lanczos方法，有望提升稳定性和收敛性。

ABSTRACT

Lanczos methods for solving $ extit{A} extbf{x}= extbf{b}$ consist in constructing a sequence of vectors $( extbf{x}_k), k=1,...$ such that $ extbf{r}_{k}= extbf{b}- extit{A} extbf{x}_{k}= extit{P}_{k}( extit{A}) extbf{r}_{0}$,, where $ extit{P}_{k}$ is the orthogonal polynomial of degree at most $k$ with respect to the linear functional $c$ defined as $c(ξ^i)=( extbf{y}, extit{A}^i extbf{r}_{0})$. Let $ extit{P}^{(1)}_{k}$ be the regular monic polynomial of degree $k$ belonging to the family of formal orthogonal polynomials (FOP) with respect to $c^{(1)}$ defined as $c^{(1)}(ξ^{i})=c(ξ^{i+1})$. All Lanczos-type algorithms are characterized by the choice of one or two recurrence relationships, one for $ extit{P}_{k}$ and one for $ extit{P}^{(1)}_{k}$. We shall study some new recurrence relations involving $ extit{P}_{k}$ and $ extit{P}^{(1)}_{k}$ and their possible combination to obtain new Lanczos-type algorithms. We will show that some recurrence relations exist, but cannot be used to derive Lanczos-type algorithms, while others do not exist at all.

研究动机与目标

• 探索形式正交多项式（FOPs）$P_k$ 与 $P_k^{(1)}$ 之间超越标准次数差约束的新递推关系。
• 确定这些新递推关系中哪些在数学上有效且可用于构建Lanczos型算法。
• 识别并排除那些虽满足正交性和归一化条件，但不适合用于计算残差向量 $\mathbf{r}_k$ 和迭代向量 $\mathbf{x}_k$ 的关系。
• 将新递推关系分类为：存在但不可用、根本不存在，以及可用于算法设计的可行关系。
• 提供一种系统化框架，通过有效 $A_i/B_j$ 递推对的组合生成新型Lanczos型算法。

提出的方法

• 将形式正交多项式（FOP）框架扩展至分析多项式两侧次数差达三的递推关系，而非传统的一至两个次数差。
• 利用正交性和归一化条件，推导并评估 $P_k(x)$ 的递推关系 $A_i$ 和 $P_k^{(1)}(x)$ 的 $B_j$。
• 应用线性泛函 $c^{(1)}(\xi^i) = c(\xi^{i+1})$ 定义 $P_k^{(1)}$ 的内积结构，从而实现递归构造。
• 通过系数匹配与正交性约束，推导高阶递推中递推系数（如 $C_k, D_k, F_k, G_k$）的方程组。
• 采用行列式分析（$\Delta'_k$）评估系数系统的可解性，确保非奇异以实现有效递推构造。
• 提出如 $A_{13}/B_{13}$、$A_{16}/B_{16}$ 和 $A_{19}/B_{15}$ 等算法组合，从有效递推对生成新型Lanczos型求解器。

实验结果

研究问题

• RQ1 在标准次数差约束之外，哪些关于形式正交多项式 $P_k$ 与 $P_k^{(1)}$ 的新递推关系存在？
• RQ2 这些新识别出的递推关系中，哪些在数学上有效但不适合用于Lanczos型算法实现？
• RQ3 哪些递推关系因违反正交性或归一化条件而根本不存在？
• RQ4 能否系统性地组合可行的递推对（$A_i/B_j$）以生成新型Lanczos型算法？
• RQ5 在FOP-based Lanczos方法中，确保高阶递推系数系统可解性与稳定性的条件是什么？

主要发现

• 关系 $A_{11}, A_{17}, B_{11}, B_{12}$ 因正交性与归一化条件的根本不一致而不存在。
• 关系 $A_{14}, A_{15}, A_{18}, B_{14}$ 虽在数学上成立，但无法用于Lanczos型算法，因其无法生成有效的 $\mathbf{r}_k$ 和 $\mathbf{x}_k$ 向量。
• 关系 $A_{12}, A_{13}, A_{16}, A_{19}, B_{13}, B_{15}, B_{16}$ 在数学上有效且适用于算法实现。
• 关系 $A_{12}$ 具有自洽性，仅凭其自身即可生成一种新型Lanczos型算法，如文献[10:Salhi]所示。
• 可行的算法组合包括 $A_{13}/B_{13}$、$A_{13}/B_{15}$、$A_{13}/B_{16}$、$A_{16}/B_{13}$、$A_{16}/B_{15}$、$A_{16}/B_{16}$、$A_{19}/B_{13}$、$A_{19}/B_{15}$、$A_{19}/B_{16}$，所有组合均可生成新型Lanczos型求解器。
• 递推系数系统的可解性由行列式 $\Delta'_k$ 非零决定，确保 $D_k, F_k, G_k$ 系数有唯一解。