[论文解读] Decay bounds for non-Hermitian matrix functions

Stefano Pozza, Valeria Simoncini|arXiv (Cornell University)|May 5, 2016

Matrix Theory and Algorithms被引用 7

一句话总结

该论文利用切比雪夫多项式级数，推导出非厄米带状矩阵函数条目的一般先验衰减界，将已知的厄米特情形结果扩展至非厄米情形。即使在矩阵不可对角化的情况下，该研究也建立了严格的衰减估计，并将其应用于改进不精确Krylov近似，数值实验验证了边界的紧致性。

ABSTRACT

The derivation of a-priori decay bounds for the entries of functions of banded matrices is of interest in a variety of applications. While decay bounds for functions of Hermitian banded matrices have been known for some time, the non-Hermitian case is an especially challenging setting. By using Faber polynomial series we explore the bounds obtainable by extending results for Hermitian matrices to banded non-Hermitian (not necessarily diagonalizable) matrices. Several special cases are treated, together with an application to the inexact Krylov approximation of matrix function evaluations. Numerical experiments illustrate the quality of all new bounds.

研究动机与目标

• 将矩阵函数衰减界从厄米特矩阵推广至非厄米带状矩阵，此类矩阵此前尚无相关衰减界结果。
• 利用切比雪夫多项式级数建立理论框架，分析非厄米矩阵函数的条目衰减行为，包括不可对角化矩阵的情形。
• 提供可实际应用的衰减估计，以提升不精确Krylov子空间方法在矩阵函数计算中的效率与精度。
• 通过数值实验验证理论边界，展示其紧致性及在实际场景中的适用性。

提出的方法

• 该方法采用切比雪夫多项式展开来表示矩阵函数，利用函数的解析性质和矩阵的谱性质。
• 通过在复平面上进行多项式逼近，将厄米特矩阵理论中的技术扩展至非厄米情形。
• 通过分析数值域和谱凝聚效应，考虑非正规性和不可对角化性对切比雪夫级数的影响。
• 利用柯西积分公式和切比雪夫多项式系数的界，推导出条目级的衰减率。
• 结合矩阵的带状结构，定位衰减行为并推导条目特定的边界。
• 通过具有不同谱性质的测试矩阵上的矩阵函数计算，对理论边界进行数值验证。

实验结果

研究问题

• RQ1 能否为非厄米带状矩阵的矩阵函数条目严格推导出衰减界，即使矩阵不可对角化？
• RQ2 切比雪夫多项式级数如何实现将厄米特矩阵衰减结果推广至非厄米情形？
• RQ3 非正规性和谱凝聚效应对矩阵函数条目衰减行为有何影响？
• RQ4 在实际中，所推导的边界有多紧致？与数值实验中观测到的衰减行为相比如何？
• RQ5 这些边界在多大程度上能提升不精确Krylov子空间近似在矩阵函数计算中的收敛性和效率？

主要发现

• 该论文成功利用切比雪夫多项式级数，为非厄米带状矩阵函数的条目推导出条目级衰减界，将先前仅适用于厄米特情形的结果进行了扩展。
• 即使在矩阵不可对角化的情况下，该边界依然有效，解决了现有理论中的关键局限。
• 数值实验表明，所推导的边界具有紧致性，并能准确预测各类矩阵类型和函数的衰减行为。
• 边界在估计不精确Krylov近似中矩阵函数的稀疏模式和截断误差方面表现出有效性。
• 理论框架表明，衰减行为由数值域的几何形状和函数的解析区域共同决定。
• 该方法为量化矩阵函数条目的局域化提供了系统性方法，使大规模问题中的高效矩阵函数计算成为可能。