Abstract
Finite rate of innovation (FRI) is a powerful reconstruction framework enabling the recovery of sparse Dirac streams from uniform low-pass filtered samples. An extension of this framework, called generalised FRI (genFRI), has been recently proposed for handling cases with arbitrary linear measurement models. In this context, signal reconstruction amounts to solving a joint constrained optimisation problem, yielding estimates of both the Fourier series coefficients of the Dirac stream and its so-called annihilating filter, involved in the regularisation term. This optimisation problem is however highly non convex and non linear in the data. Moreover, the proposed numerical solver is computationally intensive and without convergence guarantee. In this work, we propose an implicit formulation of the genFRI problem. To this end, we leverage a novel regularisation term which does not depend explicitly on the unknown annihilating filter yet enforces sufficient structure in the solution for stable recovery. The resulting optimisation problem is still non convex, but simpler since linear in the data and with less unknowns. We solve it by means of a provably convergent proximal gradient descent (PGD) method. Since the proximal step does not admit a simple closed-form expression, we propose an inexact PGD method, coined as Cadzow plug-and-play gradient descent (CPGD). The latter approximates the proximal steps by means of Cadzow denoising, a well-known denoising algorithm in FRI. We provide local fixed-point convergence guarantees for CPGD. Through extensive numerical simulations, we demonstrate the superiority of CPGD against the state-of-the-art in the case of non uniform time samples.
Highlights
S AMPLING theorems lie at the foundation of modern digital signal processing as they permit the convenient navigation between the analogue and digital worlds [1], [2]
This fact was brought to the attention of the signal processing community in [4], where the authors introduced the finite rate of innovation (FRI) framework
We propose an implicit version of the generalised finite-rateof-innovation problem for the recovery of the Fourier series coefficients of sparse Dirac streams with arbitrary linear sensing
Summary
S AMPLING theorems lie at the foundation of modern digital signal processing as they permit the convenient navigation between the analogue and digital worlds [1], [2]. The gain in accuracy comes at the price of significantly higher computational cost, the Douglas-Rachford splitting method requiring many more iterations to converge than Cadzow algorithm In addition to their somewhat heuristic nature, neither Cadzow denoising nor its upgrade can handle more general types of input measurements as considered in the generalised FRI (genFRI) framework introduced by Pan et al in [23]. We propose an implicit formulation of the genFRI problem in which only the Fourier coefficients to be annihilated are recovered This formulation does not rely explicitly on the unknown annihilating filter but rather leverages a structured low-rank regularisation constraint based on a “Toeplitzification” linear operator, guaranteeing non-trivial solutions to the annihilating equation. All experiments and simulations are fully reproducible using the benchmarking routines provided in our GitHub repository [30]
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.