Algebraic description of the finite Stieltjes moment problem

  • ID: 2801, RIV: 10393271
  • ISSN: 0024-3795, ISBN: not specified
  • source: Linear Algebra and Its Applications
  • keywords: Stieltjes moment problem; Continued fractions; Hankel matrices; Jacobi matrices; Cholesky factorization of moment matrices
  • authors: Stefano Pozza, Zdeněk Strakoš
  • authors from KNM: Strakoš Zdeněk

Abstract

The Stieltjes problem of moments seeks for a nondecreasing positive distribution function mu(lambda) on the semi-axis [0, +infinity) so that its moments match a given infinite sequence of positive real numbers m(0), m(l), . . . . In his seminal paper Investigations on continued fractions published in 1894 Stieltjes gave a complete solution including the conditions for the existence and uniqueness in relation to his main goal, the convergence theory of continued fractions. One can also reformulate the Stieltjes problem of moments as looking for a sequence of positive distribution functions mu((1))(lambda), mu((2))(lambda), . . . , where the nth distribution function has n points of increase and, m(0), m(1), . . . , m(2n-1 )represent its (first) 2n moments, i.e., as the sequence of the finite Stieltjes moment problems. This view can be linked to iterative solution of (large) linear algebraic systems. Providing that m(0), m(1), . . . , are moments of some linear, self-adjoint and coercive operator A on a Hilbert space with respect to a given vector f , the finite Stieltjes moment problems determine the iterations of the conjugate gradient method applied for solving Au = f, and vice versa. Here the existence and uniqueness is guaranteed by the properties of the operator A (reformulation for finite sequences, matrices and finite vectors is obvious). This fundamental link raises a question on how the solution of the finite Stieltjes moment problem can be described purely algebraically. This has motivated the presented exposition built upon ideas published previously by several authors. Since the description uses matrices of moments, it is not intended for numerical computations. (C) 2018 Elsevier Inc. All rights reserved.