# orthogonal polynomials

## 1 Orthogonal Polynomials

Polynomials of degree $n$ are analytic functions^{} that can be written in the form

$${p}_{n}(x)={a}_{0}+{a}_{1}x+{a}_{2}{x}^{2}+\mathrm{\cdots}+{a}_{n}{x}^{n}$$ |

They can be differentiated and integrated for any value of $x$, and are fully determined by the $n+1$ coefficients ${a}_{0}\mathrm{\dots}{a}_{n}$ . For this simplicity they are frequently used to approximate more complicated or unknown functions^{}. In approximations, the necessary degree $n$ of the polynomial is not normally defined by criteria other than the quality of the approximation.

Using polynomials as defined above tends to lead into numerical difficulties when determining the ${a}_{i}$, even for small values of $n$. It is therefore customary to stabilize results numerically by using orthogonal polynomials over an interval $[a,b]$, defined with respect to a positive weighting function $W(x)>0$ by

$${\int}_{a}^{b}{p}_{n}(x){p}_{m}(x)W(x)\mathit{d}x=0\text{for}n\ne m$$ |

Orthogonal polynomials are obtained in the following way: define the scalar product.

$$(f,g)={\int}_{a}^{b}f(x)g(x)W(x)\mathit{d}x$$ |

between the functions $f$ and $g$, where $W(x)$ is a weight factor. Starting with the polynomials ${p}_{0}(x)=1$, ${p}_{1}(x)=x$, ${p}_{2}(x)={x}^{2}$, etc., from the Gram-Schmidt decomposition one obtains a sequence of orthogonal polynomials ${q}_{0}(x),{q}_{1}(x),\mathrm{\dots}$, such that $({q}_{m},{q}_{n})={N}_{n}{\delta}_{mn}$. The normalization factors ${N}_{n}$ are arbitrary. When all ${N}_{i}$ are equal to one, the polynomials are called orthonormal.

Some important orthogonal polynomials are:

$a$ | $b$ | $W(x)$ | name |
---|---|---|---|

-1 | 1 | 1 | Legendre polynomials^{} |

-1 | 1 | ${(1-{x}^{2})}^{-1/2}$ | Chebyshev polynomials |

$-\mathrm{\infty}$ | $\mathrm{\infty}$ | ${e}^{-{x}^{2}}$ | Hermite polynomials^{} |

Orthogonal polynomials of successive orders can be expressed by a recurrence relation

$${p}_{n}=({A}_{n}+{B}_{n}x){p}_{n-1}+{C}_{n}{p}_{n-2}$$ |

This relation can be used to compute a finite series

$${a}_{0}{p}_{0}+{a}_{1}{p}_{1}+\mathrm{\cdots}+{a}_{n}{p}_{n}$$ |

with arbitrary coefficients ${a}_{i}$, without computing explicitly every polynomial ${p}_{j}$ (Horner’s Rule).

Chebyshev polynomials ${T}_{n}(x)$ are also orthogonal with respect to discrete values ${x}_{i}$:

$$ |

where the ${x}_{i}$ depend on $M$.

## References

- 1 Originally from The Data Analysis Briefbook (http://rkb.home.cern.ch/rkb/titleA.htmlhttp://rkb.home.cern.ch/rkb/titleA.html)
- 2 M. Abramowitz and I.A. Stegun (Eds.), Handbook of Mathematical Functions, National Bureau of Standards, Dover, New York, 1974.
- 3 W.H. Press, S.A. Teukolsky, W.T. Vetterling, and B.P. Flannery, Numerical Recipes in C, Second edition, Cambridge University Press, 1995. (The same book exists for the Fortran language). There is also an Internet version which you can work from.

Title | orthogonal polynomials |
---|---|

Canonical name | OrthogonalPolynomials |

Date of creation | 2013-03-22 12:06:20 |

Last modified on | 2013-03-22 12:06:20 |

Owner | akrowne (2) |

Last modified by | akrowne (2) |

Numerical id | 11 |

Author | akrowne (2) |

Entry type | Definition |

Classification | msc 42C05 |

Classification | msc 33D45 |

Related topic | OrthogonalMatrices |

Related topic | Orthonormal |