Speciale Functies – Orthogonale polynomen – Legendre polynomen

De Legendre polynomen zijn orthogonaal op het interval \((-1,1)\) met betrekking tot de gewichtsfunctie \(w(x)=1\). Deze kunnen worden gedefinieerd door middel van hun Rodrigues formule:

\[P_n(x)=\frac{(-1)^n}{2^n\,n!}\,\frac{1}{w(x)}\,D^n\left[w(x)\,(1-x^2)^n\right] =\frac{(-1)^n}{2^n\,n!}\,D^n\left[(1-x^2)^n\right],\quad n=0,1,2,\ldots.\tag1\]

Dit is een speciaal geval \(\alpha=\beta=0\) van de Jacobi polynomen:

\[P_n(x)=P_n^{(0,0)}(x)={}_2F_1\left(\genfrac{}{}{0pt}{}{-n,n+1}{1}\,;\,\frac{1-x}{2}\right),\quad n=0,1,2,\ldots.\tag2\]

Verder geldt

\[P_n(-x)=(-1)^nP_n(x),\quad P_n(1)=1\quad\text{en}\quad P_n(-1)=(-1)^n,\quad n=0,1,2,\ldots.\]

De kopcoëfficiënt van het polynoom \(P_n(x)\) is gelijk aan

\[k_n=\frac{(-n)_n\,(n+1)_n}{(1)_n\,n!}\,\frac{(-1)^n}{2^n}=\frac{(2n)!}{2^n\,(n!)^2},\quad n=0,1,2,\ldots.\]

De orthogonaliteitsrelatie is

\[\int_{-1}^1P_m(x)P_n(x)\,dx=\frac{2}{2n+1}\,\delta_{mn},\quad m,n\in\{0,1,2,\ldots\}.\tag3\]

Dit kan worden aangetoond met behulp van de Rodrigues formule (1) en partiële integratie

\[\int_{-1}^1P_m(x)P_n(x)\,dx=\frac{(-1)^n}{2^n\,n!}\,\int_{-1}^1P_m(x)D^n\left[(1-x^2)^n\right]\,dx =\frac{1}{2^n\,n!}\,\int_{-1}^1D^nP_m(x)\,(1-x^2)^n\,dx,\]

hetgeen nul is voor \(m < n\). Voor \(m=n\) geldt

\[\int_{-1}^1D^nP_n(x)\,(1-x^2)^n\,dx=k_n\,n!\,\int_{-1}^1(1-x^2)^n\,dx=\frac{(2n)!}{2^n\,n!}\,\int_{-1}^1(1-x^2)^n\,dx.\]

Uiteindelijk volgt dan met de substitutie \(1-x=2t\)

\begin{align*} \int_{-1}^1(1-x^2)^n\,dx&=\int_{-1}^1(1-x)^n(1+x)^n\,dx=\int_0^1(2t)^n(2-2t)^n\,2\,dx=2^{2n+1}\,B(n+1,n+1)\\[2.5mm] &=2^{2n+1}\,\frac{\Gamma(n+1)\,\Gamma(n+1)}{\Gamma(2n+2)}=\frac{2^{2n+1}\,(n!)^2}{(2n+1)!},\quad n=0,1,2,\ldots. \end{align*}

Dus geldt

\[\int_{-1}^1\left\{P_n(x)\right\}^2\,dx=\frac{(2n)!}{2^{2n}\,(n!)^2}\,\frac{2^{2n+1}\,(n!)^2}{(2n+1)!}=\frac{2}{2n+1},\quad n=0,1,2,\ldots.\]

Dit bewijst de orthogonaliteitsrelatie (3).

Om een genererende functie voor de Legendre polynomen te vinden, gebruiken we de hypergeometrische representatie (2) en vinden

\begin{align*} \sum_{n=0}^{\infty}P_n(x)t^n&=\sum_{n=0}^{\infty}\sum_{k=0}^n\frac{(-n)_k\,(n+1)_k}{(1)_k\,k!}\left(\frac{1-x}{2}\right)^kt^n =\sum_{k=0}^{\infty}\sum_{n=k}^{\infty}\frac{(-n)_k\,(n+1)_k}{k!\,k!}\left(\frac{1-x}{2}\right)^kt^n\\[2.5mm] &=\sum_{k=0}^{\infty}\sum_{n=0}^{\infty}\frac{(-n-k)_k\,(n+k+1)_k}{k!\,k!}\left(\frac{1-x}{2}\right)^kt^{n+k} =\sum_{k=0}^{\infty}\frac{(2k)!}{k!\,k!}\left(\frac{x-1}{2}\right)^kt^k\sum_{n=0}^{\infty}\frac{(2k+1)_n}{n!}\,t^n\\[2.5mm] &=\sum_{k=0}^{\infty}\frac{(2k)!}{k!\,k!}\left(\frac{x-1}{2}\right)^kt^k(1-t)^{-2k-1} =\sum_{k=0}^{\infty}\frac{(1/2)_k}{k!}\,\left[2(x-1)t\right]^k\,(1-t)^{-2k-1}\\[2.5mm] &=(1-t)^{-1}\left[1-\frac{2(x-1)t}{(1-t)^2}\right]^{-1/2}=\left[(1-t)^2-2(x-1)t\right]^{-1/2}=(1-2xt+t^2)^{-1/2} =\frac{1}{\sqrt{1-2xt+t^2}}. \end{align*}

Dit bewijst de genererende functie

\[\frac{1}{\sqrt{1-2xt+t^2}}=\sum_{n=0}^{\infty}P_n(x)t^n.\tag4\]

Definieer nu \(F(x,t)=(1-2xt+t^2)^{-1/2}\), dan geldt

\[\frac{\partial}{\partial t}F(x,t)=-\frac{1}{2}(1-2xt+t^2)^{-3/2}\,(-2x+2t)=\frac{x-t}{(1-2xt+t^2)^{3/2}}.\]

Hieruit volgt dat

\[(1-2xt+t^2)\frac{\partial}{\partial t}F(x,t)=(x-t)F(x,t).\]

Met behulp van (4) volgt dan

\[(1-2xt+t^2)\sum_{n=1}^{\infty}nP_n(x)t^{n-1}=(x-t)\sum_{n=0}^{\infty}P_n(x)t^n.\]

Dit kan ook worden geschreven als

\[\sum_{n=1}^{\infty}nP_n(x)t^{n-1}-2x\sum_{n=1}^{\infty}nP_n(x)t^n+\sum_{n=1}^{\infty}nP_n(x)t^{n+1} =x\sum_{n=0}^{\infty}P_n(x)t^n-\sum_{n=0}^{\infty}P_n(x)t^{n+1}\]

oftewel

\[\sum_{n=0}^{\infty}(n+1)P_{n+1}(x)t^n-x\sum_{n=0}^{\infty}(2n+1)P_n(x)t^n+\sum_{n=0}^{\infty}(n+1)P_n(x)t^{n+1}=0.\]

Dit leidt tot \(P_1(x)=xP_0(x)\) en de drieterms recurrente betrekking

\[(n+1)P_{n+1}(x)-(2n+1)xP_n(x)+nP_{n-1}(x)=0,\quad n=1,2,3,\ldots.\]

De Legendre polynomen \(P_2(x)\), \(P_3(x)\) en \(P_4(x)\).

---

Last modified on 22 mei 2021