We consider the multi-dimensional Hausdorff moment problem over the unit cube: to reconstruct an unknown function from the (inaccurately) given values of the integrals of the unknown function multiplied by all power-products of the independent variables. We describe a regularization scheme using orthogonalization by the tensor product of (shifted) Legendre polynomials and "approximation" of the unknown function by a finite sum, the dimension of the space of approximation playing the role of the regularization parameter. For the case of square integrability of the unknown function we present an estimate of the regularization error that implies convergence if the data error tends to zero.
Read full abstract