[论文解读] A fast alternating projection method for complex frequency estimation

Фредрик Андерссон, Marcus Carlsson+1 more|arXiv (Cornell University)|Jul 8, 2011

Direction-of-Arrival Estimation Techniques参考文献 26被引用 5

一句话总结

该论文提出了一种用于通过复指数和估计信号中复频率的快速交替投影方法。通过在低秩对称矩阵与汉克尔矩阵之间迭代投影——利用快速傅里叶变换（FFT）和复对称Lanczos算法——该方法实现了快速收敛和精确的频率估计，在数值实验中优于MUSIC和ESPRIT方法，并在一般条件下具备理论收敛保证。

ABSTRACT

The problem of approximating a sampled function using sums of a fixed number of complex exponentials is considered. We use alternating projections between fixed rank matrices and Hankel matrices to obtain such an approximation. Convergence, convergence rates and error estimates for this technique are proven, and fast algorithms are developed. We compare the numerical results obtain with the MUSIC and ESPRIT methods.

研究动机与目标

• 开发一种快速且精确的方法，用于使用固定数量的复指数和来近似采样信号。
• 建立在复对称秩-k矩阵流形与汉克尔矩阵流形之间交替投影的收敛性与误差估计。
• 设计计算高效的投影步骤算法，利用快速变换与迭代求解器。
• 在准确性和计算速度方面，将所提方法与MUSIC和ESPRIT等成熟技术进行比较。
• 通过在一般设定下验证横截性条件，为收敛性提供理论依据。

提出的方法

• 该方法在复对称秩-k矩阵集合与汉克尔矩阵集合之间交替进行投影，以将信号近似为k个复指数的和。
• 利用由k个复指数构成的信号会产生一个秩为k的汉克尔矩阵的事实，从而可采用低秩逼近方法。
• 通过快速傅里叶变换（FFT）实现汉克尔矩阵-向量乘法与低秩投影的快速计算。
• 采用定制化的复对称版本Lanczos算法，以高效地实现对低秩矩阵的投影。
• 利用文献[2]中的结果，建立理论收敛性，验证在弱横截性条件下，交替投影方案在一般情况下收敛。
• 该方法被应用于加权函数空间，包括高斯型权重，以实现在非均匀采样设置下的逼近。

实验结果

研究问题

• RQ1 在低秩对称矩阵与汉克尔矩阵之间交替投影，是否能收敛到由k个复指数构成的信号的有意义近似？
• RQ2 在复频率估计背景下，该交替投影方案的收敛速率与误差界是什么？
• RQ3 如何利用FFT与迭代求解器加速投影步骤，使该方法适用于大规模问题？
• RQ4 所提方法在准确性和速度方面与MUSIC和ESPRIT等成熟方法相比如何？
• RQ5 在何种一般条件下，交替投影方案能收敛到接近真实信号的解？

主要发现

• 交替投影方法收敛至低秩对称矩阵流形与汉克尔矩阵流形交集中的解，且在一般横截性条件下已证明收敛。
• 通过利用基于FFT的乘法运算与复对称Lanczos算法进行低秩投影，该方法实现了快速收敛，显著降低了计算成本。
• 数值比较显示，该方法在噪声信号与有限长度信号下的精度优于MUSIC和ESPRIT。
• 推导出理论误差估计，表明近似误差有界，且依赖于初始信号与秩k。
• 该方法适用于加权函数空间，包括具有高斯型权重的空间，从而可在非均匀采样条件下实现逼近。
• 证明了参数化映射的导数的秩为2kn - k²，确认了流形结构，并支持使用横截性论证进行收敛性分析。