OnderwijsSpeciale Functies – Orthogonale polynomen – Jacobi polynomen
De Jacobi polynomen zijn orthogonaal op het interval \((-1,1)\) met betrekking tot de beta verdeling \(w(x)=(1-x)^{\alpha}(1+x)^{\beta}\). Deze kunnen worden gedefinieerd door middel van hun Rodrigues formule:
\[P_n^{(\alpha,\beta)}(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!}\,(1-x)^{-\alpha}(1+x)^{-\beta}\,D^n\left[(1-x)^{n+\alpha}(1+x)^{n+\beta}\right],\quad n=0,1,2,\ldots.\tag1\]Met behulp van de regel van Leibniz volgt
\begin{align*} D^n\left[(1-x)^{n+\alpha}(1+x)^{n+\beta}\right]&=\sum_{k=0}^n\binom{n}{k}D^k(1-x)^{n+\alpha}D^{n-k}(1+x)^{n+\beta}\\[2.5mm] &=\sum_{k=0}^n\binom{n}{k}(-1)^k(n+\alpha)(n+\alpha-1)\cdots(n+\alpha-k+1)(1-x)^{n+\alpha-k}\\[2.5mm] &{}\hspace{50mm}{}\times(n+\beta)(n+\beta-1)\cdots(\beta+k+1)(1+x)^{\beta+k}\\[2.5mm] &=n!\,\sum_{k=0}^n(-1)^k\binom{n+\alpha}{k}\binom{n+\beta}{n-k}(1-x)^{n+\alpha-k}(1+x)^{\beta+k},\quad n=0,1,2,\ldots. \end{align*}Hieruit volgt dat
\[P_n^{(\alpha,\beta)}(x)=\frac{(-1)^n}{2^n}\,\sum_{k=0}^n(-1)^k\binom{n+\alpha}{k}\binom{n+\beta}{n-k}(1-x)^{n-k}(1+x)^k,\quad n=0,1,2,\ldots.\]This shows that \(P_n^{(\alpha,\beta)}(x)\) is a polynomial of degree \(n\). Note that we have the symmetry
\[P_n^{(\alpha,\beta)}(-x)=(-1)^nP_n^{(\beta,\alpha)}(x),\quad n=0,1,2,\ldots\tag2\]en
\[P_n^{(\alpha,\beta)}(1)=\binom{n+\alpha}{n}\quad\text{en}\quad P_n^{(\alpha,\beta)}(-1)=(-1)^n\binom{n+\beta}{n},\quad n=0,1,2,\ldots.\]Om een hypergeometrische representatie te vinden schrijven we voor \(x\neq 1\)
\[P_n^{(\alpha,\beta)}(x)=\left(\frac{x-1}{2}\right)^n\sum_{k=0}^n\binom{n+\alpha}{n} \binom{n+\beta}{n-k}\left(\frac{x+1}{x-1}\right)^k,\quad n=0,1,2,\ldots.\]Nu geldt voor \(x\neq 1\) \[\left(\frac{x+1}{x-1}\right)^k=\left(1+\frac{2}{x-1}\right)^k =\sum_{i=0}^k\binom{k}{i}\left(\frac{2}{x-1}\right)^i,\quad k=0,1,2,\ldots.\]
Nu volgt door verwisseling van de volgorde van sommatie voor \(x\neq 1\)
\begin{align*} P_n^{(\alpha,\beta)}(x)&=\left(\frac{x-1}{2}\right)^n\sum_{i=0}^n\sum_{k=i}^n \binom{n+\alpha}{k}\binom{n+\beta}{n-k}\binom{k}{i}\left(\frac{2}{x-1}\right)^i\\[2.5mm] &=\left(\frac{x-1}{2}\right)^n\sum_{i=0}^n\sum_{k=0}^{n-i} \binom{n+\alpha}{i+k}\binom{n+\beta}{n-i-k}\binom{i+k}{i}\left(\frac{2}{x-1}\right)^i,\quad n=0,1,2,\ldots. \end{align*}Nu verwisselen we de volgorde in de eerste som en vinden voor \(x\neq 1\) en \(n=0,1,2,\ldots\)
\begin{align*} P_n^{(\alpha,\beta)}(x)&=\left(\frac{x-1}{2}\right)^n\sum_{i=0}^n\sum_{k=0}^n \binom{n+\alpha}{n-i+k}\binom{n+\beta}{i-k}\binom{n-i+k}{n-i}\left(\frac{2}{x-1}\right)^{n-i}\\ &=\sum_{i=0}^n\sum_{k=0}^n\binom{n+\alpha}{n-i+k}\binom{n+\beta}{i-k} \binom{n-i+k}{n-i}\left(\frac{x-1}{2}\right)^i\\ &=\sum_{i=0}^n\sum_{k=0}^n\frac{\Gamma(n+\alpha+1)}{(n-i+k)!\,\Gamma(i-k+\alpha+1)}\frac{\Gamma(n+\beta+1)}{(i-k)!\,\Gamma(n-i+k+\beta+1)}\, \frac{(n-i+k)!}{(n-i)!\,k!}\left(\frac{x-1}{2}\right)^i\\ &=\frac{\Gamma(n+\alpha+1)\Gamma(n+\beta+1)}{n!}\sum_{i=0}^n \frac{(-n)_i}{\Gamma(i+\alpha+1)\,i!\,\Gamma(n-i+\beta+1)}\left(\frac{1-x}{2}\right)^i\sum_{k=0}^n\frac{(-i)_k(-i-\alpha-1)_k}{(n-i+\beta+1)_k\,k!}. \end{align*}Omdat \(i\in\{0,1,2,\ldots,n\}\) volgt uit de Chu-Vandermonde sommatieformule
\[\sum_{k=0}^n\frac{(-i)_k(-i-\alpha-1)_k}{(n-i+\beta+1)_k\,k!} ={}_2F_1\left(\genfrac{}{}{0pt}{}{-i,-i-\alpha-1}{n-i+\beta+1}\,;\,1\right)=\frac{(n+\alpha+\beta+1)_i}{(n-i+\beta+1)_i}.\]Dus volgt met \(\Gamma(n-i+\beta+1)\,(n-i+\beta+1)_i=\Gamma(n+\beta+1)\) dat
\begin{align*} P_n^{(\alpha,\beta)}(x)&=\frac{\Gamma(n+\alpha+1)}{n!}\sum_{i=0}^n \frac{(-n)_i(n+\alpha+\beta+1)_i}{\Gamma(i+\alpha+1)\,i!}\left(\frac{1-x}{2}\right)^i\\ &=\frac{\Gamma(n+\alpha+1)}{n!\,\Gamma(\alpha+1)}\sum_{i=0}^n \frac{(-n)_i(n+\alpha+\beta+1)_i}{(\alpha+1)_i\,i!}\left(\frac{1-x}{2}\right)^i,\quad n=0,1,2,\ldots. \end{align*}Dit bewijst de hypergeometrische representatie
\[P_n^{(\alpha,\beta)}(x)=\binom{n+\alpha}{n}\,{}_2F_1\left(\genfrac{}{}{0pt}{}{-n,n+\alpha+\beta+1}{\alpha+1}\,;\,\frac{1-x}{2}\right),\quad n=0,1,2,\ldots.\]Merk op dat dit ook geldt voor \(x=1\). Vanwege de symmetrie (2) geldt ook
\[P_n^{(\alpha,\beta)}(x)=(-1)^n\binom{n+\beta}{n}\,{}_2F_1\left(\genfrac{}{}{0pt}{}{-n,n+\alpha+\beta+1}{\beta+1}\,;\,\frac{1+x}{2}\right),\quad n=0,1,2,\ldots.\]Merk op dat uit de hypergeometrische representatie volgt dat
\begin{align*} \frac{d}{dx}P_n^{(\alpha,\beta)}(x)&=\binom{n+\alpha}{n}\,\frac{(-n)(n+\alpha+\beta+1)}{\alpha+1}\cdot\left(-\frac{1}{2}\right) {}_2F_1\left(\genfrac{}{}{0pt}{}{-n+1,n+\alpha+\beta+2}{\alpha+2}\,;\,\frac{1-x}{2}\right)\\ &=\frac{n+\alpha+\beta+1}{2}\,\binom{n+\alpha}{n-1}\,{}_2F_1\left(\genfrac{}{}{0pt}{}{-n+1,n+\alpha+\beta+2}{\alpha+2}\,;\,\frac{1-x}{2}\right)\\ &=\frac{n+\alpha+\beta+1}{2}\,P_{n-1}^{(\alpha+1,\beta+1)}(x),\quad n=1,2,3,\ldots. \end{align*}Een ander gevolg van de hypergeometrische representatie is dat de kopcoëfficiënt van het polynoom \(P_n^{(\alpha,\beta)}(x)\) gelijk is aan
\[k_n=\binom{n+\alpha}{n}\,\frac{(-n)_n(n+\alpha+\beta+1)_n}{(\alpha+1)_n\,n!}\,\frac{(-1)^n}{2^n} =\frac{(n+\alpha+\beta+1)_n}{2^n\,n!},\quad n=0,1,2,\ldots.\]Nu kan worden aangetoond dat de Jacobi polynomen voldoen aan de orthogonaliteitsrelatie
\[\int_{-1}^1(1-x)^{\alpha}(1+x)^{\beta}P_m^{(\alpha,\beta)}(x)P_n^{(\alpha,\beta)}(x)\,dx =\frac{2^{\alpha+\beta+1}\,\Gamma(n+\alpha+1)\,\Gamma(n+\beta+1)} {(2n+\alpha+\beta+1)\,\Gamma(n+\alpha+\beta+1)\,n!}\,\delta_{mn}\]voor \(\alpha>-1\), \(\beta>-1\) en \(m,n\in\{0,1,2,\ldots\}\). Dit kan worden aangetoond met behulp van de definitie (1) en partiële integratie. De waarde van de integraal in het geval \(m=n\) kan worden berekend met behulp van de kopcoëffici&eum;nt en door vervolgens de integraal te schrijven in termen van een beta integraal:
\begin{align*} \int_{-1}^1(1-x)^{\alpha}(1+x)^{\beta}\left\{P_n^{(\alpha,\beta)}(x)\right\}^2\,dx &=\frac{(-1)^n}{2^n\,n!}\int_{-1}^1P_n^{(\alpha,\beta)}(x)D^n\left[(1-x)^{n+\alpha}(1+x)^{n+\beta}\right]\,dx\\ &=\frac{1}{2^n\,n!}\int_{-1}^1D^nP_n^{(\alpha,\beta)}(x)(1-x)^{n+\alpha}(1+x)^{n+\beta}\,dx\\ &=\frac{(n+\alpha+\beta+1)_n}{2^{2n}\,n!}\int_{-1}^1(1-x)^{n+\alpha}(1+x)^{n+\beta}\,dx\\ &=\frac{\Gamma(2n+\alpha+\beta+1)}{\Gamma(n+\alpha+\beta+1)\,2^{2n}\,n!}\int_{-1}^1(1-x)^{n+\alpha}(1+x)^{n+\beta}\,dx,\quad n=0,1,2,\ldots \end{align*}en met behulp van de substitutie \(1-x=2t\)
\begin{align*} \int_{-1}^1(1-x)^{n+\alpha}(1+x)^{n+\beta}\,dx &=\int_0^1(2t)^{n+\alpha}(2-2t)^{n+\beta}\,2\,dt=2^{2n+\alpha+\beta+1}\int_0^1t^{n+\alpha}(1-t)^{n+\beta}\,dt\\[2.5mm] &=2^{2n+\alpha+\beta+1}\,B(n+\alpha+1,n+\beta+1) =2^{2n+\alpha+\beta+1}\,\frac{\Gamma(n+\alpha+1)\Gamma(n+\beta+1)}{\Gamma(2n+\alpha+\beta+2)}\\ &=2^{2n+\alpha+\beta+1}\,\frac{\Gamma(n+\alpha+1)\Gamma(n+\beta+1)}{(2n+\alpha+\beta+1)\,\Gamma(2n+\alpha+\beta+1)},\quad n=0,1,2,\ldots. \end{align*}De Jacobi polynomen \(P_n^{(\alpha,\beta)}(x)\) voldoen aan de tweede orde lineaire differentiaalvergelijking
\[(1-x^2)y''(x)+\left[\beta-\alpha-(\alpha+\beta+2)x\right]y'(x)+n(n+\alpha+\beta+1)y(x)=0,\quad n=0,1,2,\ldots.\]Een genererende functie voor de Jacobi polynomen wordt gegeven door
\[\frac{2^{\alpha+\beta}}{R(1+R-t)^{\alpha}(1+R+t)^{\beta}} =\sum_{n=0}^{\infty}P_n^{(\alpha,\beta)}(x)t^n,\quad R=\sqrt{1-2xt+t^2}.\]Last modified on 22 mei 2021
