56
15

Quadratic functional estimation in inverse problems

Abstract

We consider in this paper a Gaussian sequence model of observations YiY_i, i1i\geq 1 having mean (or signal) θi\theta_i and variance σi\sigma_i which is growing polynomially like iγi^\gamma, γ>0\gamma >0. This model describes a large panel of inverse problems. We estimate the quadratic functional of the unknown signal i1θi2\sum_{i\geq 1}\theta_i^2 when the signal belongs to ellipsoids of both finite smoothness functions (polynomial weights iαi^\alpha, α>0\alpha>0) and infinite smoothness (exponential weights eβire^{\beta i^r}, β>0\beta >0, 0<r20<r \leq 2). We propose a Pinsker type projection estimator in each case and study its quadratic risk. When the signal is sufficiently smoother than the difficulty of the inverse problem (α>γ+1/4\alpha>\gamma+1/4 or in the case of exponential weights), we obtain the parametric rate and the efficiency constant associated to it. Moreover, we give upper bounds of the second order term in the risk and conjecture that they are asymptotically sharp minimax. When the signal is finitely smooth with αγ+1/4\alpha \leq \gamma +1/4, we compute non parametric upper bounds of the risk of and we presume also that the constant is asymptotically sharp.

View on arXiv
Comments on this paper