Time and covariance smoothing for restoration of bivariate signals

In many applications and physical phenomena, bivariate signals are polarized, i.e. they trace an elliptical trajectory over time when viewed in the 2D planes of their two components. The smooth evolution of this elliptical trajectory, called polarization ellipse, is highly informative to solve ill-posed inverse problems involving bivariate signals where the signal is collected through indirect, noisy or incomplete measurements. This work proposes a novel formulation and an efficient algorithm for reconstructing bivariate signals with polarization regularization. The proposed formulation leverages the compact representation of polarization through the instantaneous covariance matrices. To address the resulting quartic optimization problem, we propose a well-suited parameter splitting strategy which leads to an efficient iterative algorithm (alternating direction method of multipliers (ADMM)) with convex subproblems at each iteration. The performance of the proposed method is illustrated on numerical synthetic data experiments.