Hamburger Moment Problem
   HOME

TheInfoList



OR:

In
mathematics Mathematics is an area of knowledge that includes the topics of numbers, formulas and related structures, shapes and the spaces in which they are contained, and quantities and their changes. These topics are represented in modern mathematics ...
, the Hamburger
moment problem In mathematics, a moment problem arises as the result of trying to invert the mapping that takes a measure ''μ'' to the sequences of moments :m_n = \int_^\infty x^n \,d\mu(x)\,. More generally, one may consider :m_n = \int_^\infty M_n(x) ...
, named after
Hans Ludwig Hamburger Hans Ludwig Hamburger (5 August 1889, Berlin – 14 August 1956, Cologne) was a German mathematician. He was a professor at universities in Berlin, Cologne and Ankara.. Biography Hans was the elder son of Karl Hamburger and Margarethe Levy. He was ...
, is formulated as follows: given a
sequence In mathematics, a sequence is an enumerated collection of objects in which repetitions are allowed and order matters. Like a set, it contains members (also called ''elements'', or ''terms''). The number of elements (possibly infinite) is calle ...
(''m''0, ''m''1, ''m''2, ...), does there exist a positive
Borel measure In mathematics, specifically in measure theory, a Borel measure on a topological space is a measure that is defined on all open sets (and thus on all Borel sets). Some authors require additional restrictions on the measure, as described below. F ...
''μ'' (for instance, the measure determined by the
cumulative distribution function In probability theory and statistics, the cumulative distribution function (CDF) of a real-valued random variable X, or just distribution function of X, evaluated at x, is the probability that X will take a value less than or equal to x. Ev ...
of a
random variable A random variable (also called random quantity, aleatory variable, or stochastic variable) is a mathematical formalization of a quantity or object which depends on random events. It is a mapping or a function from possible outcomes (e.g., the po ...
) on the
real line In elementary mathematics, a number line is a picture of a graduated straight line (geometry), line that serves as visual representation of the real numbers. Every point of a number line is assumed to correspond to a real number, and every real ...
such that :m_n = \int_^\infty x^n\,d \mu(x) \text In other words, an affirmative answer to the problem means that (''m''0, ''m''1, ''m''2, ...) is the sequence of moments of some positive Borel measure ''μ''. The Stieltjes moment problem, Vorobyev moment problem, and the
Hausdorff moment problem In mathematics, the Hausdorff moment problem, named after Felix Hausdorff, asks for necessary and sufficient conditions that a given sequence be the sequence of moments :m_n = \int_0^1 x^n\,d\mu(x) of some Borel measure supported on the clo ...
are similar but replace the real line by ,+\infty) (Stieltjes and Vorobyev; but Vorobyev formulates the problem in the terms of matrix theory), or a bounded interval (Hausdorff).


Characterization

The Hamburger moment problem is solvable (that is, (''m''''n'') is a sequence of moments) if and only if the corresponding Hankel kernel on the nonnegative integers : A = \left(\begin m_0 & m_1 & m_2 & \cdots \\ m_1 & m_2 & m_3 & \cdots \\ m_2 & m_3 & m_4 & \cdots \\ \vdots & \vdots & \vdots & \ddots \end\right) is positive definite kernel, positive definite, i.e., : \sum_m_c_j\overline\ge0 for every arbitrary sequence (''c''''j'')''j'' ≥ 0 of
complex number In mathematics, a complex number is an element of a number system that extends the real numbers with a specific element denoted , called the imaginary unit and satisfying the equation i^= -1; every complex number can be expressed in the form ...
s that are finitary (i.e. ''cj'' = 0 except for finitely many values of ''j''). For the "only if" part of the claims simply note that : \sum_m_c_j \overline = \int_^\infty \left, \sum_ c_j x^j\^2\,d \mu(x) which is non-negative if \mu is non-negative. We sketch an argument for the converse. Let Z+ be the nonnegative integers and ''F''0(Z+) denote the family of complex valued sequences with finitary support. The positive Hankel kernel ''A'' induces a (possibly degenerate)
sesquilinear In mathematics, a sesquilinear form is a generalization of a bilinear form that, in turn, is a generalization of the concept of the dot product of Euclidean space. A bilinear form is linear in each of its arguments, but a sesquilinear form allows ...
product on the family of complex-valued sequences with finite support. This in turn gives a
Hilbert space In mathematics, Hilbert spaces (named after David Hilbert) allow generalizing the methods of linear algebra and calculus from (finite-dimensional) Euclidean vector spaces to spaces that may be infinite-dimensional. Hilbert spaces arise natural ...
:(\mathcal, \langle\; , \; \rangle) whose typical element is an
equivalence class In mathematics, when the elements of some set S have a notion of equivalence (formalized as an equivalence relation), then one may naturally split the set S into equivalence classes. These equivalence classes are constructed so that elements a ...
denoted by  'f'' Let ''en'' be the element in ''F''0(Z+) defined by ''en''(''m'') = ''δ''''nm''. One notices that :\langle _ _m\rangle = A_ = m_ = \langle _n _rangle. Therefore, the "shift" operator ''T'' on \mathcal, with ''T'' 'en''nbsp;=  'e''''n'' + 1 is
symmetric Symmetry (from grc, συμμετρία "agreement in dimensions, due proportion, arrangement") in everyday language refers to a sense of harmonious and beautiful proportion and balance. In mathematics, "symmetry" has a more precise definiti ...
. On the other hand, the desired expression :m_n = \int_^\infty x^n\,d \mu(x) suggests that ''μ'' is the
spectral measure In mathematics, the spectral theory of ordinary differential equations is the part of spectral theory concerned with the determination of the spectrum and eigenfunction expansion associated with a linear ordinary differential equation. In his disse ...
of a
self-adjoint operator In mathematics, a self-adjoint operator on an infinite-dimensional complex vector space ''V'' with inner product \langle\cdot,\cdot\rangle (equivalently, a Hermitian operator in the finite-dimensional case) is a linear map ''A'' (from ''V'' to itse ...
. (More precisely stated, ''μ'' is the spectral measure for an operator \overline defined below and the vector ). If we can find a "function model" such that the symmetric operator ''T'' is multiplication by ''x'', then the spectral resolution of a self-adjoint extension of ''T'' proves the claim. A function model is given by the natural isomorphism from ''F''0(Z+) to the family of
polynomial In mathematics, a polynomial is an expression consisting of indeterminates (also called variables) and coefficients, that involves only the operations of addition, subtraction, multiplication, and positive-integer powers of variables. An exa ...
s, in one single real variable and complex coefficients: for ''n'' ≥ 0, identify ''en'' with ''xn''. In the model, the operator ''T'' is multiplication by ''x'' and a densely defined symmetric operator. It can be shown that ''T'' always has self-adjoint extensions. Let \overline be one of them and ''μ'' be its spectral measure. So :\langle \overline^n \rangle = \int x^n d \mu(x). On the other hand, : \langle \overline^n \rangle = \langle T^n _0 _0\rangle = m_n. For an alternative proof of the existence that only uses Stieltjes integrals, see also, in particular theorem 3.2.


Uniqueness of solutions

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 In linear algebra, a Hankel matrix (or catalecticant matrix), named after Hermann Hankel, is a square matrix in which each ascending skew-diagonal from left to right is constant, e.g.: \qquad\begin a & b & c & d & e \\ b & c & d & e & f \\ c & d & ...
:\Delta_n = \left begin m_0 & m_1 & m_2 & \cdots & m_ \\ m_1 & m_2 & m_3 & \cdots & m_ \\ m_2 & m_3 & m_4 & \cdots & m_ \\ \vdots & \vdots & \vdots & \ddots & \vdots \\ m_ & m_ & m_ & \cdots & m_ \end\right Positivity of ''A'' means that for each ''n'', det(Δ''n'') ≥ 0. If det(Δ''n'') = 0, for some ''n'', then :(\mathcal, \langle \; , \; \rangle) is finite-dimensional and ''T'' is self-adjoint. So in this case the solution to the Hamburger moment problem is unique and ''μ'', being the spectral measure of ''T'', has finite support. More generally, the solution is unique if there are constants ''C'' and ''D'' such that for all ''n'', , ''m''''n'', ≤ ''CD''''n''''n''! . This follows from the more general
Carleman's condition In mathematics, particularly, in analysis, Carleman's condition gives a sufficient condition for the determinacy of the moment problem. That is, if a measure \mu satisfies Carleman's condition, there is no other measure \nu having the same moment ...
. There are examples where the solution is not unique; see e.g.


Further results

One can see that the Hamburger moment problem is intimately related to
orthogonal polynomials In mathematics, an orthogonal polynomial sequence is a family of polynomials such that any two different polynomials in the sequence are orthogonality, orthogonal to each other under some inner product. The most widely used orthogonal polynomial ...
on the real line. The Gram–Schmidt procedure gives a basis of orthogonal polynomials in which the operator: \overline has a tridiagonal ''Jacobi matrix representation''. This in turn leads to a ''tridiagonal model'' of positive Hankel kernels. An explicit calculation of the
Cayley transform In mathematics, the Cayley transform, named after Arthur Cayley, is any of a cluster of related things. As originally described by , the Cayley transform is a mapping between skew-symmetric matrices and special orthogonal matrices. The transform ...
of ''T'' shows the connection with what is called the '' Nevanlinna class'' of analytic functions on the left half plane. Passing to the non-commutative setting, this motivates ''Krein'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 inverse
Laplace transform In mathematics, the Laplace transform, named after its discoverer Pierre-Simon Laplace (), is an integral transform In mathematics, an integral transform maps a function from its original function space into another function space via integra ...
to the moment generating function :m(t) = \sum_m_n\frac, provided that this function converges.


References

* * * . {{DEFAULTSORT:Hamburger Moment Problem Probability problems Measure theory Functional analysis Moment (mathematics) Mathematical problems