[论文解读] Decay bounds for functions of matrices with banded or Kronecker structure

Michele Benzi, Valeria Simoncini|arXiv (Cornell University)|Jan 29, 2015

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

一句话总结

本文提出了针对带状或克罗内克和结构的厄米特矩阵的矩阵函数条目之基于积分的衰减界。这些边界显著优于以往的估计，能够捕捉非单调性和超指数衰减，尤其适用于反平方根和矩阵指数等函数。

ABSTRACT

We present decay bounds for a broad class of Hermitian matrix functions where the matrix argument is banded or a Kronecker sum of banded matrices. Besides being significantly tighter than previous estimates, the new bounds closely capture the actual (non-monotonic) decay behavior of the entries of functions of matrices with Kronecker sum structure. We also discuss extensions to more general sparse matrices.

研究动机与目标

• 解决以往衰减边界过于保守且无法捕捉矩阵函数条目实际非单调衰减行为的局限性。
• 为带状和稀疏厄米特矩阵的函数发展更精确的衰减估计，特别是具有克罗内克和结构的函数。
• 提供能准确反映整函数在克罗内克和设置下超指数衰减的边界，以及实际中观察到的振荡衰减模式。
• 提出一种基于积分表示的框架，可对特定函数进行数值计算，提升实际适用性。

提出的方法

• 利用矩阵函数的积分表示推导衰减边界，特别针对严格完全单调函数和马尔可夫函数。
• 利用拉普拉斯–斯蒂尔杰斯和柯西–斯蒂尔杰斯积分形式表达 $ f(oldsymbol{/mathcal{A}}) $ 的条目，其中 $ oldsymbol{/mathcal{A}} $ 是带状矩阵的克罗内克和。
• 通过舒尔积和张量结构，对克罗内克和情形应用矩阵指数条目的边界，将问题分解为一维分量。
• 利用已知的 $ ext{exp}(- au M) $ 条目衰减估计（$ M $ 为带状矩阵），通过双重积分界约束 $ f(oldsymbol{/mathcal{A}}) $ 的条目。
• 通过将多克罗内克积之和递归分解为两分量和，推广该方法至多克罗内克积之和的情形。
• 对所得积分进行数值计算，获得可用于实际应用的、对非对角线矩阵条目的紧致上界。

实验结果

研究问题

• RQ1 如何改进克罗内克和结构矩阵的矩阵函数衰减边界，以捕捉非单调衰减行为？
• RQ2 基于积分的表示能否为 $ f(x) = x^{-1/2} $ 和 $ f(x) = e^{-x} $ 等函数提供比现有指数边界更紧致的衰减估计？
• RQ3 在克罗内克和设置下，新边界在多大程度上能反映整函数的超指数衰减？
• RQ4 如何实际应用边界积分形式，以最小化紧致性损失，应用于特定矩阵函数？

主要发现

• 所提出的基于积分的衰减边界显著优于以往估计，尤其适用于具有非单调衰减模式的克罗内克和矩阵函数。
• 对于 $ f(x) = x^{-1/2} $，边界成功捕捉了数值实验中观察到的振荡衰减行为，如图 7 所示。
• 边界能有效反映整函数（如矩阵指数）在克罗内克和矩阵上的超指数衰减。
• 该方法允许对积分边界进行数值计算，从而在不依赖渐近或单调近似的情况下，实现对非对角线矩阵条目的实用估计。
• 通过将多克罗内克积之和递归分解为两分量和，可实现对多克罗内克积之和的扩展，同时保持边界的结构。
• 若能获得矩阵指数和预解式条目上的良好逐元素边界，该框架可推广至非厄米特矩阵，尽管此工作留待未来研究。