Polynomials

Orthogonal sets

Two polynomials are orthogonal relative to a given functional if $⟨ λ | p q ⟩ = 0$ . It is often convenient to create sets where there is one polynomial for every possible degree and polynomials of different degree are orthogonal. The further requirement that the coefficient of highest degree is 1, creates a set of “monic” polynomials. If the coefficient ring is an algebraic field, one can divide any polynomial by the value of its highest degree coefficient to give a monic polynomial. For convenience, we define $λ_{n} = ⟨ λ | x^{n} ⟩$ .

So let’s start. The monic definition requires that the initial zero-degree polynomial $p_{0} (x) = 1$ . For $n > 0$ , $⟨ λ | p_{0} p_{n} ⟩ = ⟨ λ | p_{n} ⟩ = 0$ .

For $p_{1} (x) = a + x$ , this determines $a$ : $⟨ λ | a + x ⟩ = λ_{0} a + λ_{1} = 0 ⟹ a = - λ_{1} / λ_{0}$

For this to make sense, we need $λ_{0} \neq 0$ . We could redefine the functional so that $⟨ λ | 1 ⟩ = 1$ .

We have $p_{n} (x) = x^{n} + \sum_{k = 0}^{n - 1} a_{k} p_{k} (x)$

If $m > n$ , $⟨ λ | p_{m} x^{n} ⟩ = ⟨ λ | p_{m} (p_{n} - \sum_{k = 0}^{n - 1} a_{k} p_{k}) ⟩ = 0,$ by orthogonality.

Since the leading terms of the $p_{n}$ have been chosen to be $x^{n}$ , we also have $⟨ λ | p_{n} x^{n} ⟩ = ⟨ λ | p_{n} x^{m} p_{n - m} ⟩$ and in particular, $⟨ λ | p_{n} x^{n} ⟩ = ⟨ λ | p_{n}^{2} ⟩$

We can investigate $p_{n + 1}$ in stages. Now $p_{n + 1} - x p_{n}$ is a polynomial of degree n. If $k < n - 1$ , $⟨ λ | p_{k} p_{n + 1} ⟩ = 0$ and $⟨ λ | p_{k} x p_{n} ⟩$ is a sum of terms $⟨ λ | x^{k} p_{n} ⟩ = 0$ with $k < n$ .

One of the main uses of orthogonality is to simplify the expression of a general polynomial in terms of a given monic set ${p_{n}}$ . In this, we can certainly express, $\deg (q) = M$ as: $q = \sum_{k = 0}^{M} a_{k} p_{k}$

Orthogonality then gives $⟨ λ | p_{k} q ⟩ = a_{k} ⟨ λ | p_{k}^{2} ⟩$ . So $a_{k} = ⟨ λ | p_{k} q ⟩ / ⟨ λ | p_{k}^{2} ⟩$ is well defined and unique, so long as $⟨ λ | p_{k}^{2} ⟩ \neq 0$

Applying this reasoning to $p_{n + 1} - x p_{n}$ , we already have that there are no terms $a_{k} p_{k}$ for $k < n - 1$ . Hence $p_{n + 1} = (x + a_{n}) p_{n} + b_{n} p_{n - 1}$ . This is a “three-term recurrence relation”.

Assuming we already know $λ, p_{n}, p_{n - 1}$ , we can apply orthogonality to determine $a_{n}, b_{n}$ : $0 = ⟨ λ | p_{n + 1} p_{n - 1} ⟩ = ⟨ λ | x p_{n} p_{n - 1} ⟩ + b_{n} ⟨ λ | p_{n - 1}^{2} ⟩ ⟹ b_{n} = - ⟨ λ | x p_{n} p_{n - 1} ⟩ / ⟨ λ | p_{n - 1}^{2} ⟩ = - ⟨ λ | p_{n}^{2} ⟩ / ⟨ λ | p_{n - 1}^{2} ⟩$ $0 = ⟨ λ | p_{n + 1} p_{n} ⟩ = ⟨ λ | x p_{n}^{2} ⟩ + a_{n} ⟨ λ | p_{n}^{2} ⟩ ⟹ a_{n} = - ⟨ λ | x p_{n}^{2} ⟩ / ⟨ λ | p_{n}^{2} ⟩$