Orthogonale polynomen

In de wiskunde is een stelsel orthogonale polynomen een rij polynomen van toenemende graad die onderling orthogonaal zijn met betrekking tot een of ander inproduct. Veel gebruikte en bekende stelsls zijn de hermite-polynomen, de laguerre-polynomen, de legendre-polynomen, de jacobi-polynomen en de chebyshev-polynomen. Orthogonale polynomen treden op als oplossingen van speciale differentiaalvergelijkingen en vinden toepassing in numerieke benaderingen van integralen.

Definitie

Een stelsel polynomen $P_0,P_1,\dots ,P_n,\dots$, waarin $P_n$ een polynoom van de graad $n$ is, heet orthogonaal op het interval $[a,b]$ met betrekking tot de gewichtsfunctie $w(x)$, als voor $n\ne m$ geldt:

$\underset{a}{\overset{b}{\int }}w(x)P_n(x)P_m(x)\mathrm{d}x=0$

Als de polynomen genormeerd zijn, d.w.z.

$\Vert P{\Vert }^2=\underset{a}{\overset{b}{\int }}w(x)(P_n(x){)}^2\mathrm{d}x=1$

dan heet het stelsel orthonormaal:

$\underset{a}{\overset{b}{\int }}w(x)P_n(x)P_m(x)\mathrm{d}x={\delta }_{nm}$

met ${\delta }_{nm}$ de kroneckerdelta, dus 1 als $n=m$ en 0 $n\ne m$.

Constructie

Een stelsel orthogonale polynomen kan geconstrueerd worden door orthogonalisering van de rij eenvormen $1,x,x^2,\dots ,x^n,\dots$ met behulp van de Gram-Schmidtmethode.

De coëfficiënten van de polynoom $P_{n+1}$ uit een orthogonaal stelsel volgen, op een schaalfactor na, ook uit de eis dat $P_{n+1}$ orthogonaal moet zijn met de voorgaande polynomen $P_0,P_1,\dots ,P_n$.

Voorbeeld

Een orthogonaal stelsel voor het interval $[-1,1]$ en gewichtsfunctie $w(x)=1$ wordt bepaald door:

$P_0(x)=1$

$P_1(x)=x+a=x$, want

$\underset{-1}{\overset{1}{\int }}{P}_0(x)P_1(x)\mathrm{d}x=\underset{-1}{\overset{1}{\int }}(x+a)\mathrm{d}x=2a=0$, dus $a=0$

$P_2(x)=x^2+bx+a=x^2-\frac{1}{3}$, want

$\underset{-1}{\overset{1}{\int }}{P}_0(x)P_2(x)\mathrm{d}x=\underset{-1}{\overset{1}{\int }}(x^2+bx+a)\mathrm{d}x=\frac{2}{3}+2a=0$, dus $a=-\frac{1}{3}$.

$\underset{-1}{\overset{1}{\int }}{P}_1(x)P_2(x)\mathrm{d}x=\underset{-1}{\overset{1}{\int }}x(x^2+bx+a)\mathrm{d}x=\frac{2}{3}b=0$, dus $b=0$.

Enzovoort; steeds is 1 als coëfficiënt van de hoogste macht gekozen.

De berekeningen kunnen sterk vereenvoudigd worden door de constatering dat

• vanwege de orthogonaliteit met $P_0(x)=1$, voor alle overige polynomen geldt:

$\underset{-1}{\overset{1}{\int }}{P}_n(x)\mathrm{d}x=0$

• vanwege de orthogonaliteit met $P_0(x)=1$ en $P_1(x)=x$, alle polynomen slecht uit alleen even machten of alleen oneven machten van $x$ bestaan.

Zo krijgt men:

$P_3(x)=x^3+ax=x^3-\frac{3}{5}x$, want

$\underset{-1}{\overset{1}{\int }}{P}_1(x)P_3(x)\mathrm{d}x=\underset{-1}{\overset{1}{\int }}(x^4+a{x}^2)\mathrm{d}x=\frac{2}{5}+\frac{2}{3}a=0$, dus $a=-\frac{3}{5}$.

Met de Gram-Schmidtmethode, met het inproduct

$⟨P_i,P_j⟩=\underset{-1}{\overset{1}{\int }}{P}_i(x)P_j(x)\mathrm{d}x$

krijgt men:

$P_0(x)=1$

$P_1(x)=x$, want $⟨P_0,x⟩=0$

$P_2(x)=x^2-\frac{⟨P_0,x^2⟩}{⟨P_0,P_0⟩}{P}_0(x)=x^2-\frac{\frac{2}{3}}{2}\cdot 1=x^2-\frac{1}{3}$

$P_3(x)=x^3-\frac{⟨1,x^3⟩}{2}\cdot 1-\frac{⟨x,x^3⟩}{\frac{2}{3}}x-\frac{⟨x^2-\frac{1}{3},x^3⟩}{⟨x^2-\frac{1}{3},x^2-\frac{1}{3}⟩}(x^2-\frac{1}{3})=x^3-\frac{3}{5}x$,

want $⟨x,x^3⟩=\frac{2}{5}$,$⟨1,x^3⟩=0$ en $⟨x^2-\frac{1}{3},x^3⟩=0$

Voorbeelden van orthogonale stelsels

Tabel

Integratiegrenzen		gewichtsfunctie	polynomen
$a$	$b$	$w(x)$	$R_n(x)$
$-1$	$1$	$1$	legendre-polynoom
$a$	$b$	$(x-a{)}^p(b-x{)}^q$	jacobi-polynoom
$-1$	$1$	$\frac{1}{\sqrt{1-x^2}}$	chebyshev-polynoom eerste soort
$-1$	$1$	$\sqrt{1-x^2}$	chebyshev-polynoom tweede soort
$-\mathrm{\infty }$	$\mathrm{\infty }$	$\exp (-x^2)$	hermite-polynoom
$0$	$\mathrm{\infty }$	$\exp (-x)$	laguerre-polynoom
$0$	$\mathrm{\infty }$	$x^p\exp (-x)$	geassocieerd laguerre-polynoom

De coëfficiënten van de polynomen en van hun afgeleiden zijn evenals hun nulpunten in een tabel te vinden.