Inmathematics, theHamburgermoment problem, named afterHans Ludwig Hamburger, is formulated as follows: given asequence(m0,m1,m2, ...), does there exist a positiveBorel measureμ (for instance, the measure determined by thecumulative distribution function of arandom variable) on thereal line such that
In other words, an affirmative answer to the problem means that(m0,m1,m2, ...) is the sequence ofmoments of some positive Borel measureμ.
TheStieltjes moment problem,Vorobyev moment problem, and theHausdorff moment problem are similar but replace the real line by (Stieltjes and Vorobyev; but Vorobyev formulates the problem in the terms of matrix theory), or abounded interval (Hausdorff).
The Hamburger moment problem is solvable (that is,(mn) is a sequence ofmoments) if and only if the correspondingHankel kernel on the nonnegativeintegers
ispositive definite, i.e.,
for every arbitrary sequence(cj)j ≥ 0 ofcomplex numbers that are finitary (i.e.,cj = 0 except for finitely many values ofj).
For the "only if" part of the claims simply note that
which is non-negative if is non-negative.
We sketch an argument for the converse. LetZ+ be the nonnegative integers andF0(Z+) denote the family of complex valued sequences with finitary support. The positive Hankel kernelA induces a (possibly degenerate)sesquilinear product on the family of complex-valued sequences with finite support. This in turn gives aHilbert space
whose typical element is anequivalence class denoted by[f].
Leten be the element inF0(Z+) defined byen(m) =δnm. One notices that
Therefore, theshift operatorT on, withT[en] = [en + 1], issymmetric.
On the other hand, the desired expression
suggests thatμ is thespectral measure of aself-adjoint operator. (More precisely stated,μ is the spectral measure for an operator defined below and the vector [1], (Reed & Simon 1975, p. 145)). If we can find a "function model" such that the symmetric operatorT ismultiplication byx, then the spectral resolution of aself-adjoint extension ofT proves the claim.
A function model is given by the natural isomorphism fromF0(Z+) to the family ofpolynomials, in one single real variable and complex coefficients: forn ≥ 0, identifyen withxn. In the model, the operatorT is multiplication byx and a densely defined symmetric operator. It can be shown thatT always has self-adjoint extensions. Let be one of them andμ be its spectral measure. So
On the other hand,
For an alternative proof of the existence that only usesStieltjes integrals, see also,[1] in particular theorem 3.2.
The solutions form a convex set, so the problem has either infinitely many solutions or a unique solution.
Consider the(n + 1) × (n + 1)Hankel matrix
Positivity ofA means that, for eachn,det(Δn) ≥ 0. Ifdet(Δn) = 0, for somen, then
is finite-dimensional andT is self-adjoint. So in this case the solution to the Hamburger moment problem is unique andμ, being the spectral measure ofT, has finite support.
More generally, the solution is unique if there are constantsC andD such that, for alln,|mn| ≤CDnn! (Reed & Simon 1975, p. 205). This follows from the more generalCarleman's condition.
There are examples where the solution is not unique; see e.g.[2]
The Hamburger moment problem is intimately related toorthogonal polynomials on the real line. That is, assume is the moment sequence of some positive measure on. Then for any polynomial it holds that such that the Hankel matrix is positive semidefinite. This is a necessary condition for a sequence to be a moment sequence and a sufficient condition for the existence of a positive measure.[3]
TheGram–Schmidt procedure gives a basis of orthogonal polynomials in which the operator: has a tridiagonalJacobi matrix representation. This in turn leads to atridiagonal model of positive Hankel kernels.
An explicit calculation of theCayley transform ofT shows the connection with what is called theNevanlinna class of analytic functions on the left half plane. Passing to the non-commutative setting, this motivatesKrein's formula which parametrizes the extensions of partial isometries.
The cumulative distribution function and the probability density function can often be found by applying the inverseLaplace transform to the moment generating function
provided that this function converges.
{{citation}}:ISBN / Date incompatibility (help).