[论文解读] New method of computing prolate spheroidal wavefunctions and bandlimited extrapolation: A general approach to bandlimited Fredholm kernels

Vishal Vaibhav|arXiv (Cornell University)|Apr 8, 2018

Numerical methods in inverse problems被引用 2

一句话总结

本文提出了一种新颖的退化核方法，利用球贝塞尔函数计算抛物椭球函数（PSWFs）并求解带限弗雷德霍姆积分方程。通过使用球贝塞尔基函数对核进行离散化，该方法能够高效求解相关特征值问题及带限函数的反问题，且在实轴上对一般带限、平方可积核具有高精度。

ABSTRACT

The paper deals with numerical solution of the Fredholm integral equation associated with the classical problem of extrapolating bandlimited functions known on $(-1,1)$ to the entire real line. The approach presented can be characterized as the degenerate kernel method using the spherical Bessel functions as basis functions. This discretization also facilitates the solution of the associated eigenvalue problem whose eigenfunctions are the prolate spheroidal wave functions of order zero, thus yielding a new method of computing these functions on the entire real line. These ideas are then extended to Fredholm integral equations whose kernel belongs to a class of bandlimited functions that are square integrable. The proposed discretization scheme is used to solve the associated eigenvalue problem as well as the inverse problem that arises in the estimation of object function from its image function.

研究动机与目标

• 开发一种在实轴上对零阶抛物椭球函数进行数值稳定计算的方法。
• 求解在带限函数外推中超出区间 $(-1,1)$ 的弗雷德霍姆积分方程。
• 将该方法推广至一般带限、平方可积核，以求解特征值问题与反问题。
• 通过使用球贝塞尔函数进行离散化，建立系统化框架以处理带限弗雷德霍姆核。

提出的方法

• 该方法通过将带限核展开为球贝塞尔函数作为基函数，采用退化核近似。
• 离散化将连续的弗雷德霍姆积分方程转化为矩阵特征值问题，从而实现抛物椭球函数的数值求解。
• 所得矩阵系统的特征函数近似于经典零阶PSWFs，使其可在整个实轴上计算。
• 通过应用相同的球贝塞尔基函数展开，该方法可推广至任意带限、平方可积核。
• 反问题——即从带限核下的图像估计原函数——通过相同的离散化框架求解。
• 该方法通过利用球贝塞尔函数在相关函数空间中的正交性与完备性，确保了高精度。

实验结果

研究问题

• RQ1 如何高效且准确地在实轴上而非仅在 $[-1,1]$ 上计算零阶抛物椭球函数？
• RQ2 对带限弗雷德霍姆核进行离散化的最优方式是什么，以保持其谱特性并实现相关特征值问题的数值求解？
• RQ3 能否将使用球贝塞尔函数的退化核方法推广至任意带限、平方可积核？
• RQ4 与现有数值技术相比，该方法在PSWF计算和带限外推中的精度与稳定性如何？
• RQ5 该框架在信号处理中的反问题（如带限核去卷积）中可应用到何种程度？

主要发现

• 所提方法可在整个实轴上准确计算零阶抛物椭球函数，克服了传统方法仅限于 $[-1,1]$ 的限制。
• 使用球贝塞尔函数作为基函数可导致条件良好的矩阵特征值问题，从而增强数值稳定性。
• 该方法成功求解了从其带限图像重建原函数的反问题，对一般带限核表现出鲁棒性。
• 离散化方案保持了核的带限特性，并在基维数增加时收敛于真实解。
• 该框架可扩展至任意平方可积、带限核，为信号处理与数学物理中的弗雷德霍姆方程提供通用工具。