Mean Squared Error Minimization for Inverse Moment Problems

Henrion, Didier; Lasserre, Jean; Mevissen, Martin
August 2014
Applied Mathematics & Optimization;Aug2014, Vol. 70 Issue 1, p83
Academic Journal
We consider the problem of approximating the unknown density $$u\in L^2(\Omega ,\lambda )$$ of a measure $$\mu $$ on $$\Omega \subset \mathbb {R}^n$$ , absolutely continuous with respect to some given reference measure $$\lambda $$ , only from the knowledge of finitely many moments of $$\mu $$ . Given $$d\in \mathbb {N}$$ and moments of order $$d$$ , we provide a polynomial $$p_d$$ which minimizes the mean square error $$\int (u-p)^2d\lambda $$ over all polynomials $$p$$ of degree at most $$d$$ . If there is no additional requirement, $$p_d$$ is obtained as solution of a linear system. In addition, if $$p_d$$ is expressed in the basis of polynomials that are orthonormal with respect to $$\lambda $$ , its vector of coefficients is just the vector of given moments and no computation is needed. Moreover $$p_d\rightarrow u$$ in $$L^2(\Omega ,\lambda )$$ as $$d\rightarrow \infty $$ . In general nonnegativity of $$p_d$$ is not guaranteed even though $$u$$ is nonnegative. However, with this additional nonnegativity requirement one obtains analogous results but computing $$p_d\ge 0$$ that minimizes $$\int (u-p)^2d\lambda $$ now requires solving an appropriate semidefinite program. We have tested the approach on some applications arising from the reconstruction of geometrical objects and the approximation of solutions of nonlinear differential equations. In all cases our results are significantly better than those obtained with the maximum entropy technique for estimating $$u$$ .


Related Articles

  • BIOLUMINESCENCE TOMOGRAPHY: BIOMEDICAL BACKGROUND, MATHEMATICAL THEORY, AND NUMERICAL APPROXIMATION. Weimin Han; Ge Wang // Journal of Computational Mathematics;May2008, Vol. 26 Issue 3, p324 

    Over the last couple of years molecular imaging has been rapidly developed to study physiological and pathological processes in vivo at the cellular and molecular levels. Among molecular imaging modalities, optical imaging stands out for its unique advantages, especially performance and...

  • Normalization of mean squared differences to measure agreement for continuous data. Almehrizi, Rashid // Statistical Methods in Medical Research;Oct2016, Vol. 25 Issue 5, p1955 

    Agreement among observations on two variables for reliability or validation purposes is usually assessed by the evaluation of the mean squared differences (MSD). Many transformations of MSD have been proposed to interpret and make statistical inferences about the agreement between the two...

  • A SPECTRAL REGULARIZATION METHOD FOR A HEAT EQUATION BACKWARD IN TIME ON THE PLANE. Nguyen Huy Tuan; Pham Hoang Quan // ROMAI Journal;Jun2011, Vol. 7 Issue 1, p171 

    For a two-dimensional heat conduction problem on the plane, we consider its initial boundary value problem and the related inverse problem of determining the initial temperature distribution from transient temperature measurements. The stability for this inverse problem and the error analysis...

  • ON RECONSTRUCTING AN UNKNOWN COORDINATE OF A NONLINEAR SYSTEM OF DIFFERENTIAL EQUATIONS. Blizorukova, Marina; Kuklin, Alexander; Maksimov, Vyacheslav // Opuscula Mathematica;2014, Vol. 34 Issue 2, p257 

    The paper discusses a method of auxiliary controlled models and the application of this method to solving problems of dynamical reconstruction of an unknown coordinate in a nonlinear system of differential equations. The solving algorithm, which is stable with respect to informational noises and...

  • ON A SINGULAR NONLINEAR NEUMANN PROBLEM. Chabrowski, Jan // Opuscula Mathematica;2014, Vol. 34 Issue 2, p271 

    We investigate the solvability of the Neumann problem involving two critical exponents: Sobolev and Hardy-Sobolev. We establish the existence of a solution in three cases: (i) 2 < p + 1 < 2* (s), (ii) p + 1 = 2* (s) and (iii) 2* (s) < p + 1 ≤ 2*, where 2* (s) = 2(N - s)/N-2, 0 < s < 2, and...

  • Real-time inverse hysteresis compensation of piezoelectric actuators with a modified Prandtl-Ishlinskii model. Gu, Guo-Ying; Yang, Mei-Ju; Zhu, Li-Min // Review of Scientific Instruments;Jun2012, Vol. 83 Issue 6, p065106 

    This paper presents a novel real-time inverse hysteresis compensation method for piezoelectric actuators exhibiting asymmetric hysteresis effect. The proposed method directly utilizes a modified Prandtl-Ishlinskii hysteresis model to characterize the inverse hysteresis effect of piezoelectric...

  • Semi-Iterative Method for Computing the Generalized Inverse AT,S(2). Liu, Xiaoji; Jiang, Caijing // Abstract & Applied Analysis;1/29/2015, Vol. 2015, p1 

    The main aim of this paper is to compute the generalized inverse AT,S(2) over Banach spaces by using semi-iterative method and to present the error bounds of the semi-iterative method for approximating AT,S(2).

  • Application of the RBF Method to the Estimation of Temperature on the External Surface in Laminar Pipe Flow. Shiqin Lyu; Can Wu; Sufang Zhang // Mathematical Problems in Engineering;2013, p1 

    The inverse heat conduction problem on the heat transfer characteristics of cooled/heated laminar flows through finite length thick-walled circular tubes is studied, using temperature measurements taken at several different locations within the fluid in this paper. The method of radial basis...

  • An Estimate of the Misclassification Error with Hinge and Square Loss. Sahoo, J. K.; Singh, Arindama // International Journal of Computational & Applied Mathematics;2010, Vol. 5 Issue 1, p89 

    In this paper we study how Tikhonov regularization can be effectively used in the context of learning theory, especially in support vector machine classification problems. We review formulations of classification problems in the context of ill-posed inverse problem and regularization. The main...


Read the Article


Sorry, but this item is not currently available from your library.

Try another library?
Sign out of this library

Other Topics