[论文解读] Numerical stability bounds for algebraic systems of Prony type and their accurate solution by decimation

Dmitry Batenkov|arXiv (Cornell University)|Sep 10, 2014

Image and Signal Denoising Methods参考文献 45被引用 3

一句话总结

本文为高阶Prony系统提供了逐分量的渐近稳定性界，确立了任何方法的绝对分辨率极限。提出了一种降采样技术——通过等差数列选择方程——以准确求解过定Prony系统，实现接近最优的超分辨率，且不损失数值精度。

ABSTRACT

The Prony system of equations and its higher-order (confluent) generalizations appear prominently in many theoretical and applied problems. For instance, in signal processing these systems arise in recovery of sums of Diracs from a finite number of their Fourier measurements. The accurate and robust numerical solution of Prony type systems is considered to be a challenging problem, in particular reconstructing closely spaced nonlinear parameters (``nodes'', e.g. the support of the Diracs) in the presence of perturbed data. Our first contribution is providing component-wise asymptotic estimates for the numerical condition of the high-order Prony system, when the number of equations can in general be greater than the number of unknowns. These results provide, in particular, an absolute resolution limit for any method whatsoever. Our second contribution is proposing a technique for the overdetermined Prony problem with closely spaced nodes by ``decimation'', i.e. taking subsets of the equations with indices belonging to arithmetic progressions, and subsequently solving the resulting square systems. We show that solution of a decimated system is as accurate as the solution to the full overdetermined problem. Thus, decimation provides a tool to achieve near-optimal super-resolution.

研究动机与目标

• 推导高阶Prony系统的逐分量渐近条件估计，尤其针对方程数多于未知数的过定情形。
• 确立任何求解Prony型系统的数值方法的绝对分辨率极限。
• 开发一种鲁棒的数值技术，用于求解在噪声数据下节点紧密排列的过定Prony系统。
• 证明通过等差数列选择方程的降采样方法，其解的精度与求解完整过定系统相当。

提出的方法

• 利用扰动分析，推导高阶Prony系统数值条件数的逐分量渐近界。
• 提出一种降采样策略，基于索引的等差数列选择方程子集，形成方阵子系统。
• 使用标准数值求解器求解所得方阵Prony系统，利用降采样系统的结构以保持精度。
• 证明在节点分离满足弱假设的条件下，降采样子系统的解与完整过定系统的解精度相当。
• 分析降采样子系统的稳定性与分辨率极限，关联条件数与节点聚集情况。
• 使用重合Prony系统建模具有紧密排列支撑的Dirac函数之和，与信号处理中的超分辨率密切相关。

实验结果

研究问题

• RQ1 求解具有紧密排列节点的高阶Prony系统时，其基本数值稳定性极限是什么？
• RQ2 过定Prony系统中的一组方程能否产生与求解完整系统相当的解？
• RQ3 通过等差数列选择的降采样如何影响Prony系统解的分辨率与精度？
• RQ4 任何求解Prony型系统的理论分辨率极限是什么？其如何从条件数估计中推导得出？
• RQ5 在何种条件下，降采样能保持完整过定系统解的精度？

主要发现

• 本文建立了高阶Prony系统条件数的逐分量渐近界，为任何数值方法提供了理论分辨率极限。
• 通过等差数列选择方程的降采样技术，其解的精度与完整过定系统解相当。
• 该方法通过有效缓解紧密排列节点系统中的病态性，实现了接近最优的超分辨率。
• 即使原始系统因节点聚集而过定且病态，降采样子系统的精度仍能保持。
• 从条件数估计中推导出的理论分辨率极限适用于所有求解Prony型系统的算法。
• 降采样方法能够从噪声傅里叶测量中鲁棒且准确地恢复非线性参数（如Dirac支撑）