Votre navigateur n'accepte pas le Javascript.La navigation sur ce site risque de ne pas fonctionner correctement.

223. Calculez une intégrale avec la méthode du point milieu (2/2)

Soient $a$ et $b$ deux réels tels que $a<b$ et $f$ une fonction deux fois dérivable sur l’intervalle $[a,b].$

Notez $M_2 =\mathrm{Sup}\{\vert f^{\pprime}(x)\vert, x\in[a,b]\} \in \R_{+}\cup \{+\infty\}.$

Le but de cet article est de donner une majoration explicite de l’erreur commise en remplaçant le calcul de l’intégrale $\int_{a}^b f(t)\dt$ par $(b-a)f(m)$ où $m = \frac{a+b}{2}$ désigne le milieu de l’intervalle $[a,b]$ puis en subdivisant de façon régulière l’intervalle $[a,b]$ de façon uniforme.

Utilisez la tangente à la courbe de $f$ au point d’abscisse $m$

La tangente au point d’abscisse $m$ est la représentation graphique de la fonction $g$ définie par :

g(x)=f'(m)(x-m)+f(m).

En calculant l’intégrale de $g$ sur l’intervalle $[a,b]$ vous obtenez :

\begin{align*}
\int_a^b g(t)\dt&=f'(m)\int_a^b (t-m)\dt+(b-a)f(m)\\
&=f'(m)\int_{-(b-a)/2}^{(b-a)/2}u\,\mathrm{d}u+(b-a)f(m)\\
&=f'(m)\times 0+(b-a)f(m)\\
&=(b-a)f(m).
\end{align*}

Remarque : dans le cas d’une fonction positive, il exprime le fait que l’aire du rectangle obtenu avec le point milieu est exactement l’aire obtenue avec le trapèze formé par l’axe des abscisses, la tangente susmentionnée et les droites verticales d’équation $x=a$ et $x=b.$

Ce calcul montre que :

\begin{align*}
\left\vert\int_a^b f(t)\dt-(b-a)f(m)\right\vert &\leq \left\vert\int_a^b f(t)\dt-\int_a^b g(t)\dt\right\vert\\
&\leq \left\vert\int_a^b (f(t)- g(t))\dt\right\vert\\
&\leq \int_a^b \left\vert f(t)- g(t)\right\vert\dt.
\end{align*}

Majorez l’erreur commise entre la fonction $f$ et la fonction $g$

Soit $t\in[a,b].$ Si $t=m$ comme $f(m) = g(m)$ vous déduisez que $f(m)-g(m)=\frac{f^{\pprime}(m)}{2}(m-m)^2.$ Donc il existe $c\in]a,b[$ tel que $f(t)-g(t)=\frac{f^{\pprime}(c)}{2}(t-m)^2.$

Note : cette égalité est aussi appelée « égalité de Taylor-Lagrange. »

Soit $t\in[a,b]\setminus \{m\}.$

L’idée est de montrer que cette relation reste valable.

Pour générer une dérivée seconde, vous pouvez chercher à appliquer deux fois le théorème de Rolle.

Comme $g$ est un polynôme de degré $1$, il sera nul en dérivant deux fois. Il va falloir donc poser $h(x) = f(x)-g(x)-AP(x)$ où $A$ est une constante choisie telle que $h(t)=0.$ Vous souhaitez aussi avoir $h(m)=0$, ce qui donne $f(m)-g(m)-AP(m)=0$ ce qui se produit dès que $P(m)=0.$ Or vous souhaitez prendre $P$ de degré $2$ qui sera dérivé deux fois. Un bon candidat pour $P$ consiste à prendre $(x-m)^2.$

En définitive : posez $h(x) = f(x)-g(x)-A(x-m)^2$ où $A$ est défini par $A=\frac{f(t)-g(t)}{(t-m)^2}.$ Vous avez déjà $h(t)=0.$

D’autre part, $h(m) = f(m)-g(m) = 0.$ Comme $t\neq m$ vous déduisez l’existence d’un réel $d$ strictement compris entre $t$ et $m$ tel que $h^{\prime}(d) = 0.$

D’autre part :

\begin{align*}
h^{\prime}(x) &= f^{\prime}(x)-g^{\prime}(x)-2A(x-m)\\
&= f^{\prime}(x)-f^{\prime}(m)-2A(x-m).
\end{align*}

Si bien que $h^{\prime}(m) = 0.$ Or, $d\neq m$ donc l’application du théorème de Rolle à la fonction $h^{\prime}$ fournit l’existence d’un réel $c$ strictement compris entre $m$ et $d$ tel que $h^{\pprime}(c)=0.$

Or :

h^{\pprime}(x) = f^{\pprime}(x)-2A.

Et par suite $A = \frac{f^{\pprime}(c)}{2}.$ La relation $h(t)=0$ fournit :

\begin{align*}
f(t)-g(t)-\frac{f^{\pprime}(c)}{2}(t-m)^2&=0\\
f(t)-g(t)&=\frac{f^{\pprime}(c)}{2}(t-m)^2.
\end{align*}

Ainsi vous avez montré que :

\forall t\in[a,b], \exists c\in]a,b[, f(t)-g(t)=\frac{f^{\pprime}(c)}{2}(t-m)^2.

Vous déduisez la majoration:

\forall t\in[a,b], \vert f(t)-g(t)\vert \leq \frac{M_2}{2}(t-m)^2.

Etablissez le contrôle souhaité

Vous déduisez les majorations suivantes :

\begin{align*}
\left\vert\int_a^b f(t)\dt-(b-a)f(m)\right\vert 
&\leq \int_a^b \left\vert f(t)- g(t)\right\vert\dt\\
&\leq \frac{M_2}{2}\int_a^b (t-m)^2\dt\\
&\leq \frac{M_2}{2}\int_{-(b-a)/2}^{(b-a)/2} u^2\,\mathrm{d}u\\
&\leq M_2\int_{0}^{(b-a)/2} u^2\,\mathrm{d}u\\
&\leq M_2\times \frac{((b-a)/2)^3}{3}\\
&\leq \frac{M_2(b-a)^3}{24}.
\end{align*}

Subdivisez l’intervalle $[a,b]$ et concluez

Soit $n\in\NN.$ Une subdivision uniforme de l’intervalle $[a,b]$ conduit à poser $\forall k\in\llbracket 0,n\rrbracket, a_k = a+k\frac{b-a}{n}$ et $\forall k\in\llbracket 0, n-1\rrbracket, t_k=a+\left(k+\frac{1}{2}\right) \frac{b-a}{n}.$

Du coup :

\begin{align*}
\left\vert\int_a^b f(t)\dt-\frac{b-a}{n}\sum_{k=0}^{n-1} f(t_k)\right\vert &\leq \left\vert \sum_{k=0}^{n-1}\int_{a_k}^{a_{k+1}} f(t)\dt-\sum_{k=0}^{n-1}\frac{b-a}{n} f(t_k)\right\vert \\
&\leq \left\vert \sum_{k=0}^{n-1} \left[\int_{a_k}^{a_{k+1}} f(t)\dt-\frac{b-a}{n} f(t_k)\right]\right\vert \\
&\leq  \sum_{k=0}^{n-1} \left\vert\int_{a_k}^{a_{k+1}} f(t)\dt-\frac{b-a}{n} f(t_k)\right\vert \\
&\leq  \sum_{k=0}^{n-1} \left\vert\int_{a_k}^{a_{k+1}} f(t)\dt-(a_{k+1}-a_k) f(t_k)\right\vert \\
&\leq  \sum_{k=0}^{n-1} \frac{M_2 ((b-a)/n)^3}{24} \\
&\leq  \sum_{k=0}^{n-1} \frac{M_2 (b-a)^3}{24n^3} \\
&\leq  \frac{M_2 (b-a)^3}{24n^2}.
\end{align*}

En définitive, vous avez obtenu :

\boxed{\forall n\in\NN, \left\vert\int_a^b f(t)\dt-\frac{b-a}{n}\sum_{k=0}^{n-1} f(t_k)\right\vert \leq\frac{M_2 (b-a)^3}{24n^2}.}

La méthode du point milieu fournit le remarquable résultat : une approximation au deuxième ordre du calcul d’une intégrale.

222. Calculez une intégrale avec la méthode du point milieu (1/2)

Soient $a$ et $b$ deux réels tels que $a<b$ et $f$ une fonction dérivable sur l’intervalle $[a,b].$

Le but de cet article est de donner une majoration explicite de l’erreur commise en remplaçant le calcul de l’intégrale $\int_{a}^b f(t)\dt$ par $(b-a)f(m)$ où $m = \frac{a+b}{2}$ désigne le milieu de l’intervalle $[a,b]$ puis en subdivisant de façon régulière l’intervalle $[a,b]$ de façon uniforme.

Dans la suite, vous noterez $M_1 =\mathrm{Sup}\{\vert f^{\prime}(x)\vert, x\in[a,b]\} \in \R_{+}\cup \{+\infty\}.$

Trouvez une première majoration

L’idée de départ consiste à majorer la quantité :

\left\vert\int_a^b f(t)\dt - (b-a)f(m)\right\vert.

Vous y parvenez en remarquant que $b-a$ est l’intégrale $\int_a^b \dt.$

Du coup :

\begin{align*}
\left\vert\int_a^b f(t)\dt - (b-a)f(m)\right\vert &= \left\vert\int_a^b f(t)\dt - f(m)\int_a^b \dt\right\vert \\
&= \left\vert\int_a^b f(t)\dt - \int_a^b f(m)\dt\right\vert \\
&= \left\vert\int_a^b (f(t)- f(m))\dt\right\vert \\
&= \int_a^b \left\vert f(t)- f(m) \right\vert\dt \\
\end{align*}

ll vous reste à majorer la quantité $\left\vert f(t)- f(m) \right\vert$, lorsque $t\in[a,b].$

Fixez $t\in[a,b]\setminus \{m\}$ et considérez la fonction $g$ définie sur $[a,b]$ par $g(x) = f(x)-f(m)-A(x-m)$ où $A$ est le réel défini par $\frac{f(t)-f(m)}{t-m}$, de sorte que $g(t) =0.$ Comme $g(m)=0$ et que $t\neq m$, l’application du théorème de Rolle à la fonction $g$ fournit l’existence d’un réel $c$ strictement compris entre $t$ et $m$ tel que $g'(c)=0.$ Du coup, $f^{\prime}(c)-A = 0.$ La nullité de $g(t)$ fournit l’égalité $f(t)-f(m) = f'(c)(t-m).$

Si $t=m$, l’égalité $f(t)-f(m) = f'(m)(t-m)$ est satisfaite. Vous avez donc montré :

\boxed{\forall t\in[a,b], \exists c\in]a,b[,f(t)-f(m) = f'(c)(t-m).}

Vous déduisez donc:

\boxed{\forall t\in[a,b], \vert f(t)-f(m) \vert\leq M_1\vert t-m\vert.}

Note : cette inégalité est aussi appelée « inégalité des accroissements finis. »

Vous poursuivez les majorations :

\begin{align*}
\left\vert\int_a^b f(t)\dt - (b-a)f(m)\right\vert &\leq M_1\int_a^b \vert t-m\vert\dt\\
&\leq M_1\int_{-(b-a)/2}^{(b-a)/2} \vert u\vert\,\mathrm{d}u\\
&\leq 2M_1\int_{0}^{(b-a)/2} \vert u\vert\,\mathrm{d}u\\
&\leq 2M_1\int_{0}^{(b-a)/2}  u\,\mathrm{d}u\\
&\leq 2M_1\times \frac{((b-a)/2)^2}{2}\\
&\leq\frac{M_1(b-a)^2}{4}.
\end{align*}

Soit $n\in\NN.$ Subdivisez l’intervalle $[a,b]$ en $n$ parties égales, en posant $a_k = a+k\frac{b-a}{n}$, lorsque $k$ parcourt l’intervalle $\llbracket 0, n\rrbracket.$

Pour plus de lisibilité, notez, pour tout $k\in\llbracket 0,n-1\rrbracket, t_k = a+\left(k+\frac{1}{2}\right)\frac{b-a}{n}$ les points milieux.

Soit maintenant $k\in\llbracket 0,n-1\rrbracket.$ D’après ce qui précède :

\begin{align*}
\left\vert\int_{a_k}^{a_{k+1}} f(t)\dt - \frac{b-a}{n}f\left(t_k\right)\right\vert 
&\leq\frac{M_1(b-a)^2}{4n^2}.
\end{align*}

Finalement :

\begin{align*}
\left\vert\int_{a}^{b} f(t)\dt -\frac{b-a}{n}\sum_{k=0}^{n-1} f\left(t_k\right)\right\vert 
&\leq\left\vert\sum_{k=0}^{n-1}\int_{a_k}^{a_{k+1}}f(t)\dt-\sum_{k=0}^{n-1}\frac{b-a}{n}f\left(t_k\right)\right\vert\\
&\leq\left\vert\sum_{k=0}^{n-1}\left[\int_{a_k}^{a_{k+1}}f(t)\dt-\frac{b-a}{n}f\left(t_k\right)\right]\right\vert\\
&\leq\sum_{k=0}^{n-1}\left\vert\int_{a_k}^{a_{k+1}}f(t)\dt-\frac{b-a}{n}f\left(t_k\right)\right\rvert\\
&\leq \sum_{k=0}^{n-1}\frac{M_1(b-a)^2}{4n^2}\\
&\leq\frac{M_1(b-a)^2}{4n}.
\end{align*}

Pour résumer :

\boxed{\left\vert\int_{a}^{b} f(t)\dt -\frac{b-a}{n}\sum_{k=0}^{n-1} f\left(\frac{a_k+a_{k+1}}{2}\right)\right\vert 
\leq\frac{M_1(b-a)^2}{4n}.}

A première vue, vous obtenez une majoration de la forme $O\left(\frac{1}{n}\right)$ qui est d’ordre $1.$

Cependant, vous allez voir qu’il est possible d’affiner cette majoration pour obtenir une majoration de l’erreur en $O\left(\frac{1}{n^2}\right)$ dès que vous supposez que la fonction $f$ est deux fois dérivable, et en considérant $M_2 =\mathrm{Sup}\{\vert f^{\pprime}(x)\vert, x\in[a,b]\} \in \R_{+}\cup \{+\infty\}.$

Poursuivez votre lecture en allant dans l'article 223.

221. Calculez les premiers polynômes de Legendre

La formule de Rodrigues, vue dans l'article 214 permet de définir la suite des polynômes de Legendre suivante :

\forall n\in \N, \mathscr{L}_n(X) = \frac{1}{2^n n!}\frac{\mathrm{d}^n}{\mathrm{d}X^n}\left[(X^2-1)^n\right].

Vous déduisez immédiatement que :

\begin{align*}
\mathscr{L}_0(X) &= 1\\
\mathscr{L}_1(X) &= X.
\end{align*}

Utilisez une relation de récurrence pour calculer le début des autres polynômes

Il a été établi dans l'article 220 la relation suivante :

\forall n\in\NN,   (n+1)\mathscr{L}_{n+1}(X)- (2n+1)X\mathscr{L}_{n}(X)+n\mathscr{L}_{n-1}(X)=0.

Pour $n=1$, cela s’écrit :

\begin{align*}
2\mathscr{L}_{2}(X)- 3X\mathscr{L}_{1}(X)+\mathscr{L}_{0}(X)&=0\\
2\mathscr{L}_{2}(X)- 3X^2+1&=0\\
2\mathscr{L}_{2}(X)&=3X^2-1\\
\mathscr{L}_{2}(X)&=\frac{3X^2-1}{2}.
\end{align*}

Pour $n=2$, cela s’écrit :

\begin{align*}
3\mathscr{L}_{3}(X)- 5X\mathscr{L}_{2}(X)+2\mathscr{L}_{1}(X)&=0\\
3\mathscr{L}_{3}(X)- 5X\times \frac{3X^2-1}{2}+2X&=0\\
3\mathscr{L}_{3}(X)-  \frac{15X^3-5X}{2}+2X&=0\\
6\mathscr{L}_{3}(X)- (15X^3-5X)+4X&=0\\
6\mathscr{L}_{3}(X)&=15X^3-5X-4X\\
6\mathscr{L}_{3}(X)&=15X^3-9X\\
2\mathscr{L}_{3}(X)&=5X^3-3X\\
\mathscr{L}_{3}(X)&=\frac{5X^3-3X}{2}.
\end{align*}

Pour $n=3$, cela s’écrit :

\begin{align*}
4\mathscr{L}_{4}(X)- 7X\mathscr{L}_{3}(X)+3\mathscr{L}_{2}(X)&=0\\
4\mathscr{L}_{4}(X)- 7X\times \frac{5X^3-3X}{2} +3\times \frac{3X^2-1}{2}&=0\\
4\mathscr{L}_{4}(X)+ \frac{-35X^4+21X^2}{2} +\frac{9X^2-3}{2}&=0\\
8\mathscr{L}_{4}(X)-35X^4+21X^2 +9X^2-3&=0\\
8\mathscr{L}_{4}(X)-35X^4+30X^2 -3&=0\\
8\mathscr{L}_{4}(X)&=35X^4-30X^2 +3\\
\mathscr{L}_{4}(X)&=\frac{35X^4-30X^2 +3}{8}.\\
\end{align*}

Pour $n=4$, cela s’écrit :

\begin{align*}
5\mathscr{L}_{5}(X)- 9X\mathscr{L}_{4}(X)+4\mathscr{L}_{3}(X)&=0\\
5\mathscr{L}_{5}(X)- 9X\times \frac{35X^4-30X^2 +3}{8}+4\times \frac{5X^3-3X}{2}&=0\\
5\mathscr{L}_{5}(X)+ \frac{-315X^5+270X^3 -27X}{8}+2\times (5X^3-3X)&=0\\
40\mathscr{L}_{5}(X)-315X^5+270X^3 -27X+16\times (5X^3-3X)&=0\\
40\mathscr{L}_{5}(X)-315X^5+270X^3 -27X+80X^3-48X&=0\\
40\mathscr{L}_{5}(X)-315X^5+350X^3 -75X&=0\\
8\mathscr{L}_{5}(X)-63X^5+70X^3 -15X&=0\\
8\mathscr{L}_{5}(X)&=63X^5-70X^3 +15X\\
\mathscr{L}_{5}(X)&=\frac{63X^5-70X^3 +15X}{8}.
\end{align*}

220. Déterminez une relation de récurrence entre les polynômes de Legendre

Dans le cadre du contenu écrit dans l'article 214, les polynômes de Legendre peuvent être définis par la formule de Rodrigues :

\forall n\in \N, \mathscr{L}_n(X) = \frac{1}{2^n n!}\frac{\mathrm{d}^n}{\mathrm{d}X^n}\left[(X^2-1)^n\right].

Il a été démontré, via les paragraphes contenus dans l'article 216, que ces polynômes sont orthogonaux relativement au produit scalaire défini par :

\forall (P,Q)\in (\R[X])^2, \langle P, Q \rangle = \int_{-1}^1 P(t)Q(t)\,\mathrm{d}t.

D’autre part, il a été établi dans l'article 218 que les normes vérifient :

\forall n\in\N, \Vert  \mathscr{L}_n\Vert^2=\frac{2}{2n+1}.

Déterminez le coefficient dominant

Soit $n\in\N.$ Le coefficient dominant du polynôme $(X^2-1)^n$ est $1$ il est attaché au monôme $X^{2n}.$ En dérivant $n$ fois, vous déduisez que le coefficient dominant du polynôme $\frac{\mathrm{d}^n}{\mathrm{d}X^n}\left[(X^2-1)^n\right]$ est $\frac{(2n)!}{n!}$ attaché au monôme $X^n.$ Multipliant par $\frac{1}{2^n n!}$, vous en tirez que le coefficient dominant de $\mathscr{L_n}$ est :

\frac{1}{2^n n!}\times \frac{(2n)!}{n!} = \boxed{\frac{(2n)!}{2^n(n!)^2}.}

Décomposez le polynôme $X\mathscr{L}_{n}(X)$ sur les polynômes de Legendre précédents

Soit $n\in\N.$ La famille $(\mathscr{L}_{0},\dots, \mathscr{L}_{n+1})$ est étagée de degré en degré, par conséquent c’est une base de l’espace des polynômes réels dont les degrés sont inférieurs ou égaux à $n+1.$

Or, le polynôme $X\mathscr{L}_{n}(X)$ est de degré $n+1.$

Il existe donc $(a_0,\dots,a_{n+1})\in\R^{n+2}$ tel que :

X\mathscr{L}_{n}(X) = \sum_{k=0}^{n+1} a_k \mathscr{L}_{k}(X).

Le coefficient dominant du polynôme $\mathscr{L}_{n}(X)$ est $\frac{(2n)!}{2^n(n!)^2}$ il en est donc de même du polynôme $X\mathscr{L}_{n}(X).$

Le coefficient dominant du polynôme $\mathscr{L}_{n+1}(X)$ est $\frac{(2n+2)!}{2\times 2^n (n+1)^2 (n!)^2}.$

Par identification du coefficient dominant, il vient :

\begin{align*}
\frac{(2n)!}{2^n(n!)^2} &= a_{n+1}\times \frac{(2n+2)!}{2\times 2^n (n+1)^2 (n!)^2}\\
\frac{(2n)!}{2^n} &= a_{n+1}\times \frac{(2n+2)!}{2\times 2^n (n+1)^2}\\
(2n)! &= a_{n+1}\times \frac{(2n+2)!}{2 (n+1)^2}\\
1 &= a_{n+1}\times \frac{(2n+2)(2n+1)}{2 (n+1)^2}\\
1 &= a_{n+1}\times \frac{2(n+1)(2n+1)}{2 (n+1)^2}\\
1 &= a_{n+1}\times \frac{2(2n+1)}{2 (n+1)}\\
1 &= a_{n+1}\times \frac{2n+1}{n+1}\\
\frac{n+1}{2n+1} &=a_{n+1}.
\end{align*}

De ce qui précède, vous avez :

\boxed{\forall n\in\N, \exists (a_0,\dots,a_n)\in\R^{n+1},X\mathscr{L}_{n}(X) =\frac{n+1}{2n+1}\mathscr{L}_{n+1}(X)+ \sum_{k=0}^{n} a_k \mathscr{L}_{k}(X).}

L’orthogonalité des polynômes de Legendre va avoir pour conséquence étonnante le fait que les coefficients $a_k$, pour $k$ décrivant l’intervalle d’entiers $\llbracket 0, n-2\rrbracket$ sont tous nuls.

Démontrez que $\forall n\geq 2, \forall k\in\llbracket 0, n-2\rrbracket, a_k = 0$

Soit $n\in\N$ tel que $n\geq 2$ et soit $k\in \llbracket 0, n-2\rrbracket.$ Posez $a_{n+1} = \frac{n+1}{2n+1}.$ Il existe $(a_0,\dots,a_n)\in\R^{n+1}$ tel que :

\begin{aligned}
\langle X\mathscr{L}_{n}(X) , \mathscr{L}_{k}\rangle &= \left\langle \sum_{\ell=0}^{n+1} a_{\ell} \mathscr{L}_{\ell}, \mathscr{L}_{k}\right\rangle\\
&=a_{\ell} \sum_{\ell=0}^{n+1} \langle\mathscr{L}_{\ell}, \mathscr{L}_{k}\rangle\\
&=a_{k}\Vert \mathscr{L}_{k} \Vert^2\\
&=\frac{2a_{k}}{2k+1}.
\end{aligned}

Or :

\begin{aligned}
\langle X\mathscr{L}_{n}(X) , \mathscr{L}_{k}\rangle &= \int_{-1}^{1} t\mathscr{L}_{n}(t)\mathscr{L}_{k}(t)\dt\\
 &= \int_{-1}^{1} \mathscr{L}_{n}(t)(t\mathscr{L}_{k}(t))\dt\\
&=\langle \mathscr{L}_{n}, X\mathscr{L}_{k}(X)\rangle.
\end{aligned}

Le polynôme $X\mathscr{L}_{k}(X)$ est de degré $k+1$, donc il existe $(b_0\dots, b_{k+1})\in\R^{k+2}$ tel que :

X\mathscr{L}_{k}(X) = \sum_{\ell=0}^{k+1} b_{\ell} \mathscr{L}_{\ell}(X).

Utilisant la linéarité du produit scalaire :

\begin{align*}
\langle X\mathscr{L}_{n}(X) , \mathscr{L}_{k}\rangle &=\langle \mathscr{L}_{n}, X\mathscr{L}_{k}(X)\rangle \\
&= \left\langle\mathscr{L}_{n}, \sum_{\ell=0}^{k+1} b_{\ell} \mathscr{L}_{\ell}(X)\right\rangle\\
&= b_{\ell} \sum_{\ell=0}^{k+1} \langle\mathscr{L}_{n}, \mathscr{L}_{\ell}\rangle
\end{align*}

Or, pour tout $\ell\in\llbracket 0, k+1\rrbracket$ vous avez $\ell\in\llbracket 0, n-1\rrbracket$ et donc $\forall \ell\in\llbracket 0, k+1\rrbracket, n\neq \ell.$

Par orthogonalité des polynômes de Legendre, il vient $\langle X\mathscr{L}_{n}(X) , \mathscr{L}_{k}\rangle = 0.$ Du coup, $\frac{2a_{k}}{2k+1}=0$ et $a_k = 0.$

Formez la récurrence à deux crans

D’après ce qui précède, vous avez obtenu :

\boxed{\forall n\in\NN, \exists (a_{n-1}, a_n)\in\R^2,X\mathscr{L}_{n}(X) = \frac{n+1}{2n+1}\mathscr{L}_{n+1}(X)+a_{n}\mathscr{L}_{n}(X)+a_{n-1}\mathscr{L}_{n-1}(X).}

Montrez que $a_n = 0$

Soit $n$ un entier supérieur ou égal à $1.$ Posez $a_{n+1} = \frac{n+1}{2n+1}$ et utilisez l’existence de $(a_{n-1},a_n)\in\R^2$ dans la relation précédente et le produit scalaire à droite avec $\mathscr{L}_{n}$ :

\begin{align*}
\langle X\mathscr{L}_{n}(X), \mathscr{L}_{n}\rangle &= \langle a_{n+1}\mathscr{L}_{n+1}(X)+a_{n}\mathscr{L}_{n}(X)+a_{n-1}\mathscr{L}_{n-1}(X), \mathscr{L}_{n} \rangle\\
&=a_{n+1}\langle \mathscr{L}_{n+1}, \mathscr{L}_{n}\rangle+a_n\Vert \mathscr{L}_{n}\Vert^2+a_{n-1}\langle \mathscr{L}_{n-1}, \mathscr{L}_{n}\rangle\\
&=a_n\Vert \mathscr{L}_{n}\Vert^2\\
&=\frac{2a_n}{2n+1}.
\end{align*}

Utilisant la définition du produit scalaire :

\begin{align*}
\langle X\mathscr{L}_{n}(X), \mathscr{L}_{n}\rangle &= \int_{-1}^{1} t(\mathscr{L}_{n}(t))^2\dt.
\end{align*}

Or, la fonction $t\mapsto t(\mathscr{L}_{n}(t))^2$ possède une parité.

Plus précisément, si vous posez $Q(X) = (X^2-1)^n$, vous avez une fonction paire : $Q(X)=Q(-X).$ Du coup, en dérivant, pour tout $k\in\N, Q^{(k)}(X) = (-1)^k Q(-X)$ et par suite $\mathscr{L}_{n}(X) = (-1)^n \mathscr{L}_{n}(-X)$ d’où $(\mathscr{L}_{n}(X))^2 = (\mathscr{L}_{n}(-X))^2.$

Vous déduisez que la fonction $t\mapsto t(\mathscr{L}_{n}(t))^2$ est polynômiale et impaire.

Du coup, $\int_{-1}^{1} t(\mathscr{L}_{n}(t))^2\dt = 0$ puis $\frac{2a_n}{2n+1}=0$ et ainsi, $a_n = 0.$

Finalement :

\boxed{\forall n\in\NN, \exists a_{n-1}\in\R, X\mathscr{L}_{n}(X) = \frac{n+1}{2n+1}\mathscr{L}_{n+1}(X)+a_{n-1}\mathscr{L}_{n-1}(X).}

Calculez $a_{n-1}$

Soit $n$ un entier supérieur ou égal à $1.$ Posez encore $a_{n+1} = \frac{n+1}{2n+1}$ et utilisez l’existence de $a_{n-1}\in\R$ dans la relation précédente et le produit scalaire à droite avec $\mathscr{L}_{n-1}$ pour obtenir :

\begin{align*}
\langle X\mathscr{L}_{n}(X), \mathscr{L}_{n-1}\rangle &= \langle a_{n+1}\mathscr{L}_{n+1}+a_{n-1}\mathscr{L}_{n-1}, \mathscr{L}_{n-1} \rangle\\
&= a_{n+1}\langle \mathscr{L}_{n+1},\mathscr{L}_{n-1}\rangle +a_{n-1}\Vert\mathscr{L}_{n-1}\Vert^2\\
&= a_{n-1}\Vert\mathscr{L}_{n-1}\Vert^2\\
&=\frac{2a_{n-1}}{2n-1}.
\end{align*}

D’autre part :

\begin{align*}
\langle X\mathscr{L}_{n}(X), \mathscr{L}_{n-1}\rangle &= \int_{-1}^{1}t\mathscr{L}_{n}(t)\mathscr{L}_{n-1}(t)\dt\\
&=\int_{-1}^{1}\mathscr{L}_{n}(t)(t\mathscr{L}_{n-1}(t))\dt\\
&=\langle \mathscr{L}_{n}, X\mathscr{L}_{n-1}(X)\rangle.
\end{align*}

Premier cas : $n\geq 2.$ Il existe un couple $(u,v)\in\R^2$ tel que :

X\mathscr{L}_{n-1}(X) = \frac{n}{2n-1}\mathscr{L}_{n}(X)+u\mathscr{L}_{n-1}(X)+v\mathscr{L}_{n-2}(X).

Utilisant le produit scalaire à gauche avec $\mathscr{L}_{n}$, il vient :

\begin{align*}
\langle \mathscr{L}_{n}, X\mathscr{L}_{n-1}(X)\rangle &= \left\langle \mathscr{L}_{n}, \frac{n}{2n-1}\mathscr{L}_{n}(X)+u\mathscr{L}_{n-1}(X)+v\mathscr{L}_{n-2}(X) \right\rangle\\
&=\frac{n}{2n-1}\Vert\mathscr{L}_{n}\Vert^2+u\langle\mathscr{L}_{n}, \mathscr{L}_{n-1}\rangle+v\langle \mathscr{L}_{n}
,\mathscr{L}_{n-2}\rangle\\
&=\frac{n}{2n-1}\Vert\mathscr{L}_{n}\Vert^2\\
&=\frac{n}{2n-1}\times \frac{2}{2n+1}.
\end{align*}

Second cas : $n=1.$ Le polynôme $X\mathscr{L}_{n-1}(X)$ est égal à $X = \mathscr{L}_{1}(X).$

\begin{align*}
\langle \mathscr{L}_{n}, X\mathscr{L}_{n-1}(X)\rangle &= \Vert  \mathscr{L}_{1} \Vert^2\\
&=\frac{2}{3}\\
&=\frac{n}{2n-1}\times \frac{2}{2n+1}.
\end{align*}

Vous déduisez dans tous les cas que :

\begin{align*}
\frac{2a_{n-1}}{2n-1} &= \frac{2n}{(2n-1)(2n+1)}\\
2a_{n-1} &= \frac{2n}{2n+1}\\
a_{n-1} &=  \frac{n}{2n+1}.
\end{align*}

Concluez

D’après ce qui précède, vous avez :

\forall n\in\NN,  X\mathscr{L}_{n}(X) = \frac{n+1}{2n+1}\mathscr{L}_{n+1}(X)+\frac{n}{2n+1}\mathscr{L}_{n-1}(X).

Et finalement, vous obtenez la relation de récurrence suivante liant les polynômes de Legendre entre eux :

\boxed{\forall n\in\NN,   (n+1)\mathscr{L}_{n+1}(X)- (2n+1)X\mathscr{L}_{n}(X)+n\mathscr{L}_{n-1}(X)=0.}

219. Evaluez une intégrale qui intervient dans les polynômes de Legendre

Pour tout entier naturel $n$, vous posez $I_n = \int_{-1}^{1} (t^2-1)^n\,\mathrm{d}t.$

Cette intégrale apparaît dans le calcul de la norme des polynômes de Legendre. Pour en savoir davantage, reportez-vous au contenu écrit dans l'article 218.

Calculez $I_0$

\begin{aligned}
I_0 &= \int_{-1}^{1} (t^2-1)^0\,\mathrm{d}t \\
&= \int_{-1}^{1} \mathrm{d}t\\
&=2.
\end{aligned}

Vous avez $\boxed{I_0 = 2.}$

Trouvez une relation entre $I_{n+1}$ et $I_n$

Soit $n\in\N.$ Partez de $I_{n+1}$ et effectuez une intégration par parties.

\begin{aligned}
I_{n+1} &= \int_{-1}^{1} (t^2-1)^{n+1}\,\mathrm{d}t \\
&= \int_{-1}^{1} 1 \times (t^2-1)^{n+1}\,\mathrm{d}t \\
&=\left[t \times (t^2-1)^{n+1}\right]_{-1}^{1}  - \int_{-1}^{1} t  (n+1)\times (t^2-1)^{n} \times 2t \,\mathrm{d}t \\
&=  -2(n+1)  \int_{-1}^{1} t^2  (t^2-1)^{n}  \,\mathrm{d}t \\
&=  -2(n+1)  \int_{-1}^{1} (t^2-1+1)  (t^2-1)^{n}  \,\mathrm{d}t \\
&=  -2(n+1)  \left[\int_{-1}^{1} (t^2-1)  (t^2-1)^{n}\,\mathrm{d}t +\int_{-1}^{1}  (t^2-1)^{n}\,\mathrm{d}t \right]  \\
&=  -2(n+1)  \left[\int_{-1}^{1}   (t^2-1)^{n+1}\,\mathrm{d}t + I_n \right]  \\
&=  -2(n+1)  (I_{n+1} + I_n).
\end{aligned}

Donc :

$(1 + 2(n+1))I_{n+1} = -(2n+2)I_n.$

Finalement, $\boxed{\forall n\in\N, I_{n+1} = -\frac{2n+2}{2n+3}I_n.}$

Déduisez-en la valeur de $I_n$ avec une récurrence

L’expression précédente suggère la conduction d’une récurrence. Par souci de commodité, il sera utile de faire appel à la notion de factorielle double.

Pour tout entier naturel $n\geq 2$, notez $(2n)!!$ le produit de tous les entiers pairs compris entre $2$ et $2n$, avec la convention $0!!=1$ et $2!!=2.$

De même, pour tout entier naturel $n\geq 1$ notez $(2n+1)!!$ le produit de tous les entiers impairs compris entre $1$ et $2n+1$, avec la convention $1!!=1.$

Pour tout entier naturel $n$, notez $\mathscr{P}(n)$ le prédicat :
« $I_n = 2 (-1)^n \frac{(2n)!!}{(2n+1)!!}$. »

Initialisation. Pour $n=0$ :

\begin{aligned}
2 (-1)^n \frac{(2n)!!}{(2n+1)!!} &= 2 (-1)^0 \frac{0!!}{(1)!!} \\
&= 2\\
&= I_0.
\end{aligned}

Donc $\mathscr{P}(0)$ est vérifié.

Hérédité. Soit $n$ un entier naturel. Supposez $\mathscr{P}(n).$

\begin{aligned}
I_{n+1} &= -\frac{2n+2}{2n+3}I_n \\
&=-2(-1)^n\frac{2n+2}{2n+3} \times \frac{(2n)!!}{(2n+1)!!}\\
&=2(-1)^{n+1} \times \frac{(2n+2)!!}{(2n+3)!!}.
\end{aligned}

Donc $\mathscr{P}(n+1)$ est vérifié.

Ainsi, il est établi que $\forall n\in\N, I_n = 2 (-1)^n \frac{(2n)!!}{(2n+1)!!}.$

Ecrivez $I_n$ avec des factorielles

Soit $n\in\N.$ Vous allez utiliser le fait que le produit des pairs avec les impairs donne le produit de tous les entiers. Autrement dit, $(2n+1)! = (2n+1)!!(2n)!!.$ Du coup :

\begin{aligned}
\frac{(2n)!!}{(2n+1)!!} &= \frac{(2n)!! (2n)!!}{(2n+1)!!(2n)!!} \\
&= \frac{((2n)!!)^2}{(2n+1)!}\\
&= \frac{(2^n n!)^2}{(2n+1)!}\\
&= \frac{2^{2n} (n!)^2}{(2n+1)!}.
\end{aligned}

Vous déduisez $\boxed{\forall n\in\N, \int_{-1}^{1} (t^2-1)^n\,\mathrm{d}t = (-1)^{n} \frac{2^{2n+1} (n!)^2}{(2n+1)!}.}$

218. Calculez la norme d’un polynôme de Legendre

Dans le cadre du contenu se trouvant dans l'article 214, les polynômes de Legendre sont définis, en utilisant la formule de Rodrigues, par :

$\forall n\in \N, \mathscr{L}_n(X) = \frac{1}{2^n n!}\frac{\mathrm{d}^n}{\mathrm{d}X^n}\left[(X^2-1)^n\right].$

L’espace vectoriel $\R[X]$ est muni du produit scalaire suivant :

$\forall (P,Q)\in(\R[X])^2, \langle P, Q \rangle = \int_{-1}^{1} P(t)Q(t)\,\mathrm{d}t.$

La norme associée est définie par :

$\forall P\in\R[X], \Vert P \Vert = \sqrt{\langle P, P \rangle}.$

L’objectif de cet article est de trouver un moyen qui permet de calculer, quand $n\in\N$, le nombre $\Vert \mathscr{L}_n \Vert^2.$

Dans la suite de ce paragraphe, vous fixez $n\in\N.$

Utilisez des notations adaptées

Posez $P(X) = (X^2-1)^n$ et $c = \frac{1}{2^n n!}$, de sorte que $\mathscr{L}_n(X) = cP^{(n)}(X).$

Le calcul de la norme au carré fournit :

\begin{aligned}
\Vert \mathscr{L}_n \Vert^2 &= \langle \mathscr{L}_n, \mathscr{L}_n \rangle \\
&= c^2 \langle P^{(n)}, P^{(n)} \rangle.
\end{aligned}

Dérivée d’ordre $2n$

Utilisant la formule du binôme, $P(X) = \sum_{k=0}^n \binom{n}{k}(-1)^{n-k}X^{2k}.$
Isolant un terme des autres et posant $A(X) = \sum_{k=0}^{n-1} \binom{n}{k}(-1)^{n-k}X^{2k}$ il vient $P(X) = X^{2n} + A(X).$

En dérivant $2n$ fois, le polynôme $A$ qui a un degré strictement inférieur à $2n$, est nul.
Il vient donc $\boxed{P^{(2n)}(X) = (2n)!.}$

Evaluation en $1$ et en $-1$ de certaines dérivées

Notez que $P(X) = (X+1)^n(X-1)^n$ ce qui montre que $1$ est racine de $P$ avec une multiplicité égale à $n.$

Ainsi, de part l’étude des racines multiples effectuée dans l'article 217, $\forall k\in \llbracket 0,n-1 \rrbracket , P^{(k)}(1)=0.$

D’autre part, $P(-X) = P(X).$

En dérivant, il vient $P^{(n)}(X) = (-1)^n P(-X)$ donc en définitive :

$\boxed{\forall k\in \llbracket 0,n-1 \rrbracket, P^{(k)}(-1)=P^{(k)}(1) = 0.}$

Utilisez une récurrence limitée

En utilisant des intégrations par parties successives, vous allez intégrer le $P^{(n)}$ de gauche et dériver $P^{(n)}$. En intégrant $n$ fois le polynôme $P^{(n)}$, vous allez arriver sur $P.$ En dérivant $n$ fois le polynôme $P^{(n)}$ vous allez tomber sur $P^{(2n)}$ qui est un polynôme constant, puisque $P$ est de degré $2n.$ Cela conduit à la récurrence limitée suivante.

Pour tout entier $k$ compris entre $0$ et $n$, vous notez $\mathscr{P}(k)$ le prédicat :
« $\Vert \mathscr{L}_n \Vert^2 = c^2 (-1)^k\langle P^{(n-k)}, P^{(n+k)} \rangle.$ »

Initialisation. Quand $k = 0$ :

\begin{aligned}
c^2 (-1)^k\langle P^{(n-k)}, P^{(n+k)} \rangle &= c^2 (-1)^0\langle P^{(n-0)}, P^{(n+0)} \rangle \\
&= c^2 \langle P^{(n)}, P^{(n)} \rangle \\
&= \Vert \mathscr{L}_n \Vert^2.
\end{aligned}

Hérédité. Soit $k$ un entier compris entre $0$ et $n-1.$ Supposez $\mathscr{P}(k).$ Commencez par une intégration par parties :

\begin{aligned}
\langle P^{(n-k)}, P^{(n+k)} \rangle &= \int_{-1}^{1} P^{(n-k)}(t) P^{(n+k)}(t)\,\mathrm{d}t \\
&= \left[P^{(n-k-1)}(t) P^{(n+k)}(t) \right]_{-1}^{1}- \int_{-1}^{1} P^{(n-k-1)}(t) P^{(n+k+1)}(t)\,\mathrm{d}t.
\end{aligned}

Notez que, comme $0\leq k\leq n-1$, le nombre $n-k-1$ est compris entre $0$ et $n-1$, donc $P^{(n-k-1)}(1) = P^{(n-k-1)}(-1) = 0.$

\begin{aligned}
\langle P^{(n-k)}, P^{(n+k)} \rangle &= – \int_{-1}^{1} P^{(n-k-1)}(t) P^{(n+k+1)}(t)\,\mathrm{d}t \\
&= – \langle P^{(n-k-1)}, P^{(n+k+1)} \rangle\\
&= – \langle P^{(n-(k+1))}, P^{(n+(k+1))}\rangle.
\end{aligned}

Du coup :

\begin{aligned}
\Vert \mathscr{L}_n \Vert^2 &= c^2 (-1)^k\langle P^{(n-k)}, P^{(n+k)} \rangle \\
&= – c^2 (-1)^k \langle P^{(n-(k+1))}, P^{(n+(k+1))}\rangle\\
&= c^2 (-1)^{k+1}\langle P^{(n-(k+1))}, P^{(n+(k+1))}\rangle.
\end{aligned}

Le prédicat $\mathscr{P}(k+1)$ est vérifié.

D’après ce qui précède, le prédicat $\mathscr{P}(n)$ est vérifié. Donc :

\begin{aligned}
\Vert \mathscr{L}_n \Vert^2 &= c^2 (-1)^n \langle P, P^{(2n)} \rangle \\
&= (2n)! c^2 (-1)^n \langle P, 1 \rangle \\
&= (2n)! c^2 (-1)^n \int_{-1}^{1} (t^2-1)^n\,\mathrm{d}t.
\end{aligned}

L’évaluation de l’intégrale $\int_{-1}^{1} (t^2-1)^n\,\mathrm{d}t$ est effectuée dans l'article 219. Vous avez $\int_{-1}^{1} (t^2-1)^n\,\mathrm{d}t = (-1)^{n} \frac{2^{2n+1} (n!)^2}{(2n+1)!}.$

\begin{aligned}
\Vert \mathscr{L}_n \Vert^2 &= (2n)! c^2 (-1)^n (-1)^{n} \frac{2^{2n+1} (n!)^2}{(2n+1)!} \\
&= (2n)! \times \frac{1}{2^{2n} (n!)^2} \times \frac{2^{2n+1} (n!)^2}{(2n+1)!}\\
&= (2n)! \times \frac{1}{2^{2n} } \times \frac{2^{2n+1} }{(2n+1)!}\\
&= \frac{ (2n)!}{(2n+1)!} \times \frac{2^{2n+1}}{2^{2n} } \\
&= \frac{2}{2n+1}.
\end{aligned}

Ainsi, $\boxed{\forall n\in\N, \Vert \mathscr{L}_n \Vert^2 = \frac{2}{2n+1}.}$

217. Caractérisez les racines multiples d’un polynôme avec ses dérivées successives

Par définition, un polynôme $P\in\R[X]$ admet une racine $a\in\R$ de multiplicité $n\in\N^{*}$, si et seulement si, il existe un polynôme $Q\in\R[X]$ tel que $P(X)=(X-a)^nQ(X)$ et $Q(a)\neq 0.$

Démontrez que les $n$ dérivées successives du polynôme $P$ sont nulles et que $P^{(n)}(a)$ est non nul

Soit $n$ un entier naturel non nul et $P$ un polynôme de $\R[X]$ admettant une racine $a\in\R$ de multiplicité $n.$

Alors il existe un polynôme $Q\in\R[X]$ tel que $P(X)=(X-a)^nQ(X)$ et $Q(a)\neq 0.$

Par commodité, il sera utile de noter $T(X) = (X-a)^n$ et de remarquer que, pour tout entier naturel $k\in[0,n-1], T^{(k)}(X) = \frac{n!}{(n-k)!}(X-a)^{n-k}.$

Partez de $P = TQ.$

Soit maintenant $k\in[0,n-1]\in \N$ et appliquez la formule de Leibniz.

\begin{aligned}
P^{(k)}(X) &= \sum_{i=0}^k \binom{k}{i} T^{(i)}(X) Q^{(k-i)}(X).
\end{aligned}

Du coup :

\begin{aligned}
P^{(k)}(a) &= \sum_{i=0}^k \binom{k}{i} T^{(i)}(a) Q^{(k-i)}(a).
\end{aligned}

Pour tout $i\in[0,k]\cap \N, i\in[0,n-1]$ donc $T^{(i)}(a) = 0.$

Il vient $P^{(k)}(a) = 0.$

Vous avez donc montré que $\boxed{\forall k\in[0,n-1]\in\N, P^{(k)}(a) = 0.}$

De plus,

\begin{aligned}
P^{(n)}(a) &= \sum_{i=0}^n \binom{n}{i} T^{(i)}(a) Q^{(n-i)}(a) \\
&= \sum_{i=0}^{n-1} \binom{n}{i} T^{(i)}(a) Q^{(n-i)}(a) + \binom{n}{n} T^{(n)}(a) Q^{(0)}(a) \\
&= \sum_{i=0}^{n-1}0 \times Q^{(n-i)}(a) + n! Q(a) \\
&= n! Q(a).
\end{aligned}

Et par suite $\boxed{P^{(n)}(a) \neq 0.}$

Etablissez la réciproque

Etablissez d’abord un lemme

Pour tout entier naturel $n\geq 1$, notez $\mathscr{P}(n)$ le prédicat :

« $\forall a\in\R, \forall P\in\R[X], (\forall k\in[0,n-1]\cap \N, P^{(k)}(a) = 0) \implies (\exists Q\in\R[X], P(X) = (X-a)^nQ(X).)$ »

Initialisation. Supposez $n=1$, soit $a\in\R$ et soit $P\in\R[X]$ tel que $P(a)=0.$

Le nombre $a$ étant racine de $P$, il existe $Q\in\R[X]$ tel que $P(X)=(X-a)Q(X).$

Le prédicat $\mathscr{P}(1)$ est vérifié.

Hérédité. Soit $n$ un entier naturel tel que $n\geq 1.$ Supposez que $\mathscr{P}(n)$ est vérifié.

Soient $a\in\R$ et $P\in\R[X]$ tels que $\forall k\in[0,n]\cap \N, P^{(k)}(a) = 0.$

Il convient de remarquer d’abord que $\forall k\in[0,n-1]\cap \N, P^{(k)}(a) = 0.$

D’après l’hypothèse de récurrence, il existe un polynôme $Q\in\R[X], P(X) = (X-a)^{n}Q(X).$

En utilisant la formule de Leibniz, et en évaluant en $a$, il vient $P^{(n)}(a) = n! Q(a).$ Comme $P^{(n)}(a) = 0$, vous avez $Q(a) = 0.$ Donc il existe un polynôme $R\in\R[X], Q(X) = (X-a)R(X).$

Mis bout à bout, il vient $P(X) = (X-a)^{n+1}R(X)$ et le prédicat $\mathscr{P}(n+1)$ est vérifié.

Par récurrence, vous avez montré que $\forall n\in\N^{*}, \forall a\in\R, \forall P\in\R[X], (\forall k\in[0,n-1]\cap \N, P^{(k)}(a) = 0) \implies (\exists Q\in\R[X], P(X) = (X-a)^nQ(X).)$

Concluez

Soit $n$ un entier naturel non nul, soit $a\in\R$ et soit $P\in\R[X].$

Supposez que $\forall k\in[0,n-1]\cap \N, P^{(k)}(a) = 0$ et que $P^{(n)}(a)\neq 0.$

D’après le lemme établi ci-dessus, il existe $Q\in\R[X]$ tel que $P(X) = (X-a)^nQ(X).$

Or, toujours d’après la formule de Leibniz, et après évaluation en $a$, il vient $P^{(n)}(a) = n! Q(a).$

Comme $P^{(n)}(a) \neq 0$, vous avez $Q(a)\neq 0$ et la réciproque est démontrée.

Enoncez le résultat final

D’après ce qui a été établi, vous avez l’équivalence suivante : pour tout polynôme $P\in\R[X]$, pour tout réel $a$ et pour tout nombre entier naturel $n$ non nul, le nombre $a$ est racine du polynôme $P$ avec une multiplicité égale à $n$, si et seulement si, $P^{(n)}(a)\neq 0$ et pour tout entier naturel $k$ inférieur ou égal à $n-1, P^{(k)}(a) = 0.$

Prolongement

Le résultat de cet article subsiste-t-il si on remplace l’ensemble $\R$ par l’ensemble $\C$ ou par un corps $\K$ commutatif ?

216. Orthogonalité des polynômes de Legendre

Dans le cadre de la suite du contenu écrit dans l'article 214, vous définissez les polynômes de Legendre par :

$\forall n\in \N, \mathscr{L}_n(X) = \frac{1}{2^n n!}\frac{\mathrm{d}^n}{\mathrm{d}X^n}\left[(X^2-1)^n\right].$

Dans cette étude, l’orthogonalité est relative au produit scalaire suivant, défini sur l’ensemble des polynômes réels $\R[X]$ :

$\forall (P,Q)\in (\R[X])^2, \langle P, Q \rangle = \int_{-1}^1 P(t)Q(t)\,\mathrm{d}t.$

Démontrez l’orthogonalité

Soient $n$ et $m$ deux entiers naturels tels que $n>m.$ Il s’agit de comprendre pourquoi $\langle \mathscr{L}_n, \mathscr{L}_m \rangle = 0.$

Par définition du produit scalaire :

\begin{aligned}
\langle \mathscr{L}_n, \mathscr{L}_m \rangle &= \int_{-1}^1 \mathscr{L}_n(t) \mathscr{L}_m(t)\,\mathrm{d}t\\
&= \frac{1}{2^n n!} \int_{-1}^1 \frac{\mathrm{d}^n}{\mathrm{d}t^n}\left[(t^2-1)^n\right] \mathscr{L}_m(t)\,\mathrm{d}t.
\end{aligned}

Posez $I = \int_{-1}^1 \frac{\mathrm{d}^n}{\mathrm{d}t^n}\left[(t^2-1)^n\right] \mathscr{L}_m(t)\,\mathrm{d}t$ et remarquez que $ \langle \mathscr{L}_n, \mathscr{L}_m \rangle = \frac{I}{2^n n!}.$

Afin d’évaluer l’intégrale $I$, il sera utile de poser $Q(t) = (t^2-1)^n$ et d’intégrer par parties.

La clé du raisonnement est de constater que $Q(t) = (t-1)^n(t+1)^n$, ce qui prouve que le polynôme $Q$ est factorisable par $(t-1)^n$ et par $(t+1)^n$, d’où $\forall k\in[0,n-1]\cap\N, Q^{(k)}(1) = Q^{(k)}(-1)=0.$

Le calcul suivant va inciter à conduire une récurrence limitée, l’idée étant de faire remonter le nombre de dérivations imposées au polynôme $Q$, dérivant par la suite le polynôme $\mathscr{L}_m.$

\begin{aligned}
I &= \int_{-1}^1 \frac{\mathrm{d}^n}{\mathrm{d}t^n}\left[(t^2-1)^n\right] \mathscr{L}_m(t)\,\mathrm{d}t\\
&= \int_{-1}^1 Q^{(n)}(t) \mathscr{L}_m(t)\,\mathrm{d}t\\
&= \langle Q^{(n)}, \mathscr{L}_m \rangle\\
&=\left[ Q^{(n-1)}(t) \mathscr{L}_m(t)\right]_{-1}^1 – \int_{-1}^1 Q^{(n-1)}(t) (\mathscr{L}_m)'(t)\,\mathrm{d}t\\
&= – \int_{-1}^1 Q^{(n-1)}(t) (\mathscr{L}_m)'(t)\,\mathrm{d}t\\
&= – \langle Q^{(n-1)}, \mathscr{L}_m’ \rangle.
\end{aligned}

Vous allez maintenant considérer, pour tout entier $k\in[0,n]$ le prédicat suivant $P(k)$ : « $I = (-1)^k \langle Q^{(n-k)}, \mathscr{L}_m^{(k)} \rangle.$ »

Pour $k=0$ :

\begin{aligned}
(-1)^0 \langle Q^{(n-0)}, (\mathscr{L}_m)^{(0)} \rangle &= \langle Q^{(n)}, \mathscr{L}_m \rangle \\
&=I.
\end{aligned}

Donc $P(0)$ est vérifié.

Hérédité. Soit $k$ un entier appartenant à l’intervalle $[0,n-1]$. Supposez $P(k).$

\begin{aligned}
I &= (-1)^k \langle Q^{(n-k)}, \mathscr{L}_m^{(k)} \rangle\\
&= (-1)^k \int_{-1}^1 Q^{(n-k)}(t) (\mathscr{L}_m)^{(k)}(t)\,\mathrm{d}t\\
&= (-1)^k \left[ Q^{(n-k-1)}(t) (\mathscr{L}_m(t))^{(k)}\right]_{-1}^1 – (-1)^k \int_{-1}^1 Q^{(n-k-1)}(t) (\mathscr{L}_m)^{(k+1)}(t)\,\mathrm{d}t.
\end{aligned}

Comme $0\leq n-k-1 \leq n-1$, vous avez $Q^{(n-k-1)}(1)=Q^{(n-k-1)}(-1) = 0.$

D’où :

\begin{aligned}
I &= (-1)^{k+1} \int_{-1}^1 Q^{(n-(k+1))}(t) (\mathscr{L}_m)^{(k+1)}(t)\,\mathrm{d}t\\
&= (-1)^{k+1} \langle Q^{(n-(k+1))}, \mathscr{L}_m^{(k+1)} \rangle.
\end{aligned}

Le prédicat $P(k+1)$ est vérifié.

Ainsi, par récurrence limitée, $P(n)$ est vérifié.

Concluez

De ce qui précède :

\begin{aligned}
I &=(-1)^n \langle Q^{(0)}, \mathscr{L}_m^{(n)} \rangle\\
&= (-1)^n \langle Q, \mathscr{L}_m^{(n)} \rangle.
\end{aligned}

Comme le polynôme $\mathscr{L}_m$ est degré $m$ et que $n>m$, en dérivant $n$ fois ce polynôme, vous obtenez $\mathscr{L}_m^{(n)}=0.$

Par suite, $I = (-1)^n \langle Q, 0 \rangle = 0$ et donc $\boxed{\langle \mathscr{L}_n, \mathscr{L}_m \rangle = \frac{I}{2^n n!} = 0.}$

Prolongement

En utilisant l’orthogonalité des polynômes de Legendre, pourriez-vous démontrer que, pour tout entier naturel $n$ non nul, le polynôme $\mathscr{L}_n$ est scindé à racines simples ?

215. Les polynômes de Legendre sont scindés à racines simples

Dans cet article vous vous basez sur la définition des polynômes de Legendre utilisant la formule de Rodrigues vue dans l'article 214.

Pour tout entier naturel $n$, vous posez $L_n(X) = \frac{1}{2^n n!}\frac{\mathrm{d}^n}{\mathrm{d}X^n}\left[(X^2-1)^n\right].$

Il s’agit de comprendre pourquoi, pour tout entier naturel $n\geq 1$, le polynôme $L_n$, de degré $n$, admet précisément $n$ racines réelles deux à deux distinctes, appartenant toutes à l’intervalle ouvert $]-1,1[.$

Utilisez une récurrence limitée

Soit $n\geq 1$ un entier naturel fixé. Vous posez $Q(X) = (X^2-1)^n = (X-1)^n(X+1)^n.$

Pour tout entier $k$ compris entre $0$ et $n$, vous notez $P(k)$ le prédicat : « le polynôme $Q^{(k)}$ admet $k$ racines appartenant à l’intervalle $]-1,1[$. »

Représentation graphique de deux polynômes dérivés
Illustration quand $n=2$, $Q’$ possède une racine dans $]-1,1[$ et $Q^{\prime\prime}$ possède deux racines dans $]-1,1[.$

Pour $k=0$, $Q^{(k)} = Q.$ Or les racines de $Q$ sont les mêmes que celles du polynôme $X^2-1=(X-1)(X+1).$ Les deux racines étant $1$ et $-1$, le polynôme $Q$ n’admet pas de racine appartenant à l’intervalle $]-1,1[$ donc $P(0)$ est vérifié.

Soit $k$ un entier naturel compris entre $0$ et $n-1$. Vous supposez que $P(k)$ est vérifié.

Alors le polynôme $Q^{(k)}$ admet $k$ racines appartenant à l’intervalle $]-1,1[.$ Notez-les $x_1, \dots, x_k$ où $-1<x_1<\dots < x_k<1.$

Cependant, le polynôme $Q$ est factorisable par $(X-1)^n$. Comme $k$ est inférieur à $n-1$, vous avez $Q^{(k)}(1)=0.$

De même, le polynôme $Q$ est factorisable par $(X+1)^n$. Comme $k$ est inférieur à $n-1$, vous avez $Q^{(k)}(-1)=0.$

Le polynôme $Q^{(k)}$ admet donc au moins $k+2$ racines, deux à deux distinctes, qui sont $x_0 = -1$, $x_1, \dots, x_k$ et $x_{k+1}= 1.$

Pour tout $\ell$ compris entre $1$ et $k+1$, la fonction $Q^{(k)}$ est une fonction continue sur $[x_{\ell-1}, x_{\ell}]$ et dérivable sur $]x_{\ell-1}, x_{\ell}[$ puisque c’est un polynôme, avec $Q^{(k)}(x_{\ell-1}) = Q^{(k)}(x_{\ell}) = 0.$

L’application du théorème de Rolle fournit l’existence d’un réel $y_{\ell}\in]x_{\ell-1}, x_{\ell}[$ tel que $\left(Q^{(k)}\right)’ = Q^{(k+1)}(y_{\ell}) = 0.$

Vous avez les inégalités suivantes : $x_0 < y_1 < x_1 < \cdots < x_k < y_{k+1} < x_{k+1}$ ce qui prouve que les $y_{\ell}$ sont deux à deux distincts quand $\ell$ parcourt tous les entiers compris entre $1$ et $k+1.$

De ce qui précède, vous déduisez que le polynôme $Q^{(k+1)}$ admet $k+1$ racines appartenant à l’intervalle $]-1,1[$, deux à deux distinctes.

Concluez

La récurrence limitée effectuée ci-dessus montre, en particulier, que $P(n)$ est vérifié.

Du coup, le polynôme $Q^{(n)}$ possède $n$ racines deux à deux distinctes qui appartiennent toutes à l’intervalle $]-1,1[.$

Après multiplication par $\frac{1}{2^n n!}$, vous déduisez que le polynôme de Legendre $L_n$ admet bien $n$ racines deux à deux distinctes qui appartiennent toutes à l’intervalle $]-1,1[.$ Comme $L_n$ est de degré $n$, $L_n$ est bien scindé à racines simples.

Prolongement

Si vous utilisez le fait que les polynômes de Legendre sont orthogonaux pour le produit scalaire défini par $\forall (P,Q)\in\R[X]^2, \langle P, Q \rangle = \int_{-1}^1 P(t)Q(t)\,\mathrm{d}t$, pourriez-vous démontrer qu’ils sont scindés ?

214. Normalisation et construction des polynômes de Legendre à partir de la formule de Rodrigues

Suivant la formule attribuée à Rodrigues, pour tout entier naturel $n$ vous définissez le polynôme $P_n$, de degré $n$, comme étant la dérivée $n$-ième du polynôme $(X^2-1)^n$ de degré $2n$.

$\forall n\in \N, P_n(X) = \frac{\mathrm{d}^n}{\mathrm{d}X^n}\left[(X^2-1)^n\right].$

Les premiers calculs montrent que :

\begin{aligned}
P_0 &= 1\\
P_1 &= 2X\\
P_2 &= 12X^2-4.
\end{aligned}

Comme $P_1(1) = 2$, et que $P_2(1) = 8$, ces polynômes ne sont pas normalisés.

Par contre, $L_1 = \frac{1}{2}P_1$ et $L_2=\frac{1}{8}P_2$ le sont, puisque $L_1(1) = L_2(1) = 1.$

Les polynômes de Legendre sont ceux construits à partir des polynômes $P_n$ quand ils sont normalisés.

Déterminez la valeur de $P_n(1)$

Soit $n$ un entier naturel.

Accéder à $P_n(1)$ n’est pas chose aisée car le calcul explicite de $P_n$, par la formule du binôme sur $(X^2-1)^n$, aboutit à un sigma peu attractif, qu’il faudra dériver $n$ fois… ce qui aboutit à une impasse.

Une autre approche possible consiste à utiliser la formule de Taylor.

Posez $Q(X) = (X^2-1)^n = (X-1)^n(X+1)^n$, de sorte que $P_n(X) = Q^{(n)}(X).$

Alors d’une part :

\begin{aligned}
Q(X) &= Q(1)+Q'(1)(X-1)+\dots+\frac{Q^{(2n)}(1)}{(2n)!}(X-1)^{2n}\\
&=\sum_{i=0}^{2n} \frac{Q^{(i)}(1)}{i!}(X-1)^{i}\\
\end{aligned}

D’autre part, par la formule du binôme :

\begin{aligned}
(X+1)^n &= ((X-1)+2)^n\\
&=\sum_{k=0}^n \binom{n}{k}2^{n-k}(X-1)^k
\end{aligned}

D’où, après multiplication par $(X-1)^n$ :

\begin{aligned}
Q(X) &= (X-1)^n(X+1)^n \\
&=\sum_{k=0}^n \binom{n}{k}2^{n-k}(X-1)^{n+k}
\end{aligned}

L’écriture de $Q(X)$, scindée en deux, fournit des renseignements intéressants.

\begin{aligned}
Q(X) &=\sum_{i=0}^{2n} \frac{Q^{(i)}(1)}{i!}(X-1)^{i}\\
&= \sum_{i=0}^{n-1} \frac{Q^{(i)}(1)}{i!}(X-1)^{i} + \sum_{i=n}^{2n} \frac{Q^{(i)}(1)}{i!}(X-1)^{i}\\
&= \sum_{i=0}^{n-1} \frac{Q^{(i)}(1)}{i!}(X-1)^{i} + \sum_{k=0}^{n} \frac{Q^{(k+n)}(1)}{(k+n)!}(X-1)^{k+n}\\
\end{aligned}

Vous avez obtenu deux décompositions de $Q(X)$ dans la base $(1,(X-1), \dots, (X-1)^{2n})$ du $\R$-espace vectoriel formé par les polynômes dont le degré est inférieur ou égal à $2n.$

Cette famille est bien une base car les polynômes qui la forment est étagée en degré et le nombre de polynômes de cette famille est égal à $2n+1$, qui est la dimension de l’espace vectoriel précité.

Par unicité de la décomposition d’un vecteur dans une base, vous déduisez :

Pour tout entier $i$ compris entre $0$ et $n-1$, $Q^{(i)}(1) = 0$ et $\frac{Q^{(n)}(1)}{n!} = 2^n.$

Vous déduisez que $\boxed{\forall n\in \N, P_n(1) = 2^n n!.}$

Concluez

Pour tout entier naturel $n$ vous posez $\boxed{L_n(X) = \frac{1}{2^n n!}\frac{\mathrm{d}^n}{\mathrm{d}X^n}\left[(X^2-1)^n\right]}$ appelé $n$-ième polynôme de Legendre.

De cet article il résulte que, pour tout entier naturel $n$, $L_n$ est un polynôme de degré $n$ et il est normalisé, autrement dit, $L_n(1)=1.$

Prolongement

Pourriez-vous démontrer que la famille des polynômes de Legendre est orthogonale vis-à-vis du produit scalaire défini par $\forall (P,Q)\in\R[X], \langle P, Q \rangle = \int_{-1}^1 P(t)Q(t)\,\mathrm{d}t$ ?

Autrement dit, quels que soient les entiers naturels $n$ et $m$ tels que $n\leq m$, pourriez-vous démontrer que $\langle L_n, L_m \rangle = 0$ ?