A half-quadratic block-coordinate descent method for spectral estimation - CEA - Commissariat à l’énergie atomique et aux énergies alternatives Access content directly
Journal Articles Signal Processing Year : 2002

A half-quadratic block-coordinate descent method for spectral estimation


In short-time spectral estimation, Sacchi et al. (IEEE Trans. Signal Process. 46(1) (1998) 31) and Ciuciu et al. (IEEE Trans. Signal Process. 49 (2001) 2202) derived new nonlinear spectral estimators de4ned as minimizers ofpenalized criteria. The 4rst contributors have introduced separable penalizations for line spectra (LS) recovering, whereas the latter have proposed circular Gibbs–Markov functions for smooth spectra (SS) restoration, and combined both contributions for estimation of “mixed” spectra (MS), i.e., frequency peaks superimposed on a homogeneous background (Ciuciu et al., 2001). Sacchi et al. resorted to the iteratively reweighted least squares (IRLS) algorithm for the minimization stage. Here, we show that IRLS is a block-coordinate descent (BCD) method performing the minimization of a half-quadratic (HQ) energy. The latter, derived from the Geman and Reynolds construction, has the same minimizer as the initial criterion but depends on more variables. After proving that such a construction is not available for Gibbs–Markov penalizations, we extend the pioneering work ofGeman and Yang (IEEE Trans. Image Process. 4(7) (1995) 932) that leads to a suitable HQ energy for any kind ofpenalization encountered in Ciuciu et al. (2001). The BCD algorithm used for minimizing such HQ criteria is actually an original residual steepest descent (RSD) procedure (IEEE Trans. Acoust. Speech Signal Process. ASSP-33(1) (1985) 174) and thus converges in any convex case. A comparison between RSD, IRLS when available, and a pseudo-conjugate gradient algorithm is addressed in any case.
Fichier principal
Vignette du fichier
SP4339_ciuciu.pdf (473.58 Ko) Télécharger le fichier
Origin : Files produced by the author(s)

Dates and versions

cea-00333746 , version 1 (24-10-2008)


  • HAL Id : cea-00333746 , version 1


Philippe Ciuciu, J. Idier. A half-quadratic block-coordinate descent method for spectral estimation. Signal Processing, 2002, 82, pp.941-959. ⟨cea-00333746⟩
166 View
270 Download


Gmail Facebook Twitter LinkedIn More