
In
numerical analysis
Numerical analysis is the study of algorithms that use numerical approximation (as opposed to symbolic computation, symbolic manipulations) for the problems of mathematical analysis (as distinguished from discrete mathematics). It is the study of ...
, Gauss–Hermite quadrature is a form of
Gaussian quadrature
In numerical analysis, an -point Gaussian quadrature rule, named after Carl Friedrich Gauss, is a quadrature rule constructed to yield an exact result for polynomials of degree or less by a suitable choice of the nodes and weights for .
Th ...
for approximating the value of integrals of the following kind:
:
In this case
:
where ''n'' is the number of sample points used. The ''x''
''i'' are the roots of the physicists' version of the
Hermite polynomial
In mathematics, the Hermite polynomials are a classical orthogonal polynomial sequence.
The polynomials arise in:
* signal processing as Hermitian wavelets for wavelet transform analysis
* probability, such as the Edgeworth series, as well a ...
''H''
''n''(''x'') (''i'' = 1,2,...,''n''), and the associated weights ''w''
''i'' are given by
[ Abramowitz, M & Stegun, I A, ''Handbook of Mathematical Functions'', 10th printing with corrections (1972), Dover, . Equation 25.4.46.]
:
Example with change of variable
Consider a function ''h(y)'', where the variable ''y'' is
Normally distributed
In probability theory and statistics, a normal distribution or Gaussian distribution is a type of continuous probability distribution for a real number, real-valued random variable. The general form of its probability density function is
f(x ...
:
. The
expectation of ''h'' corresponds to the following integral:
As this does not exactly correspond to the Hermite polynomial, we need to change variables:
Coupled with the
integration by substitution
In calculus, integration by substitution, also known as ''u''-substitution, reverse chain rule or change of variables, is a method for evaluating integrals and antiderivatives. It is the counterpart to the chain rule for differentiation, and c ...
, we obtain:
leading to:
As an illustration, in the simplest non-trivial case, with
, we have
and
, so the estimate reduces to:
– i.e. the average of the function's values one standard deviation below and above the mean.
References
*
*
*
*
External links
* For tables of Gauss-Hermite abscissae and weights up to order ''n'' = 32 see http://www.efunda.com/math/num_integration/findgausshermite.cfm.
Generalized Gauss–Hermite quadrature free software
Free software, libre software, libreware sometimes known as freedom-respecting software is computer software distributed open-source license, under terms that allow users to run the software for any purpose as well as to study, change, distribut ...
in C++, Fortran, and Matlab
{{DEFAULTSORT:Gauss-Hermite quadrature
Numerical integration
Estimation methods