A procedure which takes a nonorthogonal set of Linearly Independent functions and constructs an Orthogonal
Basis over an arbitrary interval with respect to an arbitrary Weighting Function . Given an original set
of linearly independent functions , let denote the orthogonalized (but not normalized) functions and
the orthonormalized functions.
Take
|
(3) |
where we require
|
(4) |
By definition,
|
(5) |
so
|
(6) |
The first orthogonalized function is therefore
|
(7) |
and the corresponding normalized function is
|
(8) |
By mathematical induction, it follows that
|
(9) |
where
|
(10) |
and
|
(11) |
If the functions are normalized to instead of 1, then
|
(12) |
|
(13) |
|
(14) |
Orthogonal Polynomials are especially easy to generate using Gram-Schmidt Orthonormalization. Use the notation
|
(15) |
where is a Weighting Function, and
define the first few Polynomials,
As defined, and are Orthogonal Polynomials, as can be seen from
Now use the Recurrence Relation
|
(19) |
to construct all higher order Polynomials.
To verify that this procedure does indeed produce Orthogonal Polynomials, examine
since
. Therefore, all the Polynomials are orthogonal. Many common
Orthogonal Polynomials of mathematical physics can be generated in this manner. However, the process is
numerically unstable (Golub and van Loan 1989).
See also Gram Determinant, Gram's Inequality, Orthogonal Polynomials
References
Arfken, G. ``Gram-Schmidt Orthogonalization.'' §9.3 in Mathematical Methods for Physicists, 3rd ed.
Orlando, FL: Academic Press, pp. 516-520, 1985.
Golub, G. H. and van Loan, C. F. Matrix Computations, 3rd ed. Baltimore, MD: Johns Hopkins, 1989.
© 1996-9 Eric W. Weisstein
1999-05-25