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

236. Etudiez la nature d’une conique (1/2)

Considérez l’équation suivante :

x^2+\sqrt{3}xy+x-2=0.

En effectuant une rotation du repère, vous allez supprimer le terme en $xy$ où les variables sont croisées.

Construisez la formule de changement de repère

19/03/2022 - Rotation dun repere dun angle theta

Partant d’une base $(\vv{i}, \vv{j})$ orthonormée directe, que vous faites tourner d’un angle $\theta$ dans le sens direct, vous obtenez une nouvelle base orthonormée $(\vv{u}, \vv{v}).$

D’après les coordonnées indiquées ci-dessus, vous obtenez :

\left\{
\begin{align*}
\vv{u} &= \cos \theta \vv{i} + \sin \theta \vv{j}\\
\vv{v} &= -\sin \theta \vv{i} + \cos \theta \vv{j}.
\end{align*}\right.

Soit maintenant $M$ un point du plan. Notez $(x,y)$ les coordonnées de $M$ dans le repère $(O,\vv{i},\vv{j})$ et $(X,Y)$ celles de $M$ dans le repère $(O,\vv{u},\vv{v}).$

Les calculs fournissent successivement :

\begin{align*}
\vv{OM} &= X \vv{u} + Y \vv{v}\\
&=X( \cos \theta \vv{i} + \sin \theta \vv{j})+Y(-\sin \theta \vv{i} + \cos \theta \vv{j})\\
&=(X\cos \theta - Y\sin \theta)\vv{i}+(X\sin \theta+Y\cos \theta)\vv{j}.
\end{align*}

Le point $M$ a donc pour coordonnées $(X\cos \theta – Y\sin \theta,X\sin \theta+Y\cos \theta)$ dans le repère $(O,\vv{i},\vv{j}).$

L’unicité des coordonnées dans un repère fournit le changement suivant :

\boxed{\begin{align*}
x&=X\cos \theta - Y\sin \theta \\
y&=X\sin \theta+Y\cos \theta.
\end{align*}}

Note. Si vous faites tourner la base $(\vv{u},\vv{v})$ d’un angle de rotation égal à $-\theta$ vous vous attendez à retomber sur la base de départ $(\vv{i},\vv{j}).$ Compte tenu des relations $\cos(-\theta) = \cos (\theta)$ et $\sin(-\theta) = -\sin (\theta)$ vous vous attendez à trouver :

\begin{align*}
X&=x\cos \theta + y\sin \theta \\
Y&=-x\sin \theta+y\cos \theta.
\end{align*}

Ce résultat peut être trouvé directement grâce à l’identité fondamentale $\cos^2\theta + \sin^2\theta = 1.$

En effet :

\begin{align*}
x\cos \theta + y\sin\theta &= (X\cos \theta- Y\sin\theta)\cos \theta + (X\sin\theta+Y\cos\theta)\sin\theta\\
&= X\cos^2 \theta+ X\sin^2\theta\\
&=X(\cos^2 \theta+ \sin^2\theta)\\
&=X.
\end{align*}

De même :

\begin{align*}
-x\sin \theta + y\cos\theta &= (-X\cos\theta +Y\sin\theta)\sin\theta + (X\sin \theta+Y\cos \theta)\cos\theta\\
&= Y\sin^2\theta + Y\cos^2\theta\\
&=Y(\sin^2\theta+\cos^2\theta)\\
&=Y.
\end{align*}

Construisez une équation équivalente à celle de départ

Soit $\theta$ un réel tel que $\forall k\in\Z, \theta \neq \frac{\pi}{2}+k\pi.$

Pour tout point $M$ de coordonnées $(x,y):$

\begin{align*}
x^2+\sqrt{3}xy+x-2=0 &\Longleftrightarrow (X\cos \theta - Y\sin \theta)^2+ \sqrt{3}(X\cos \theta - Y\sin \theta)(X\sin \theta+Y\cos \theta)+X\cos \theta - Y\sin \theta - 2 = 0\\
&\Longleftrightarrow X^2\cos^2\theta+Y^2\sin^2\theta-2XY\sin\theta \cos\theta+ \sqrt{3}(X^2\sin\theta \cos\theta-Y^2\sin\theta\cos\theta + XY(\cos^2\theta - \sin^2\theta))+X\cos \theta - Y\sin \theta - 2 = 0\\
&\Longleftrightarrow (\cos^2\theta+\sqrt{3}\sin\theta\cos\theta)X^2+(\sin^2\theta-\sqrt{3}\sin\theta\cos\theta)Y^2+(\sqrt{3}(\cos^2\theta - \sin^2\theta)-2\sin\theta\cos\theta)XY+X\cos \theta - Y\sin \theta - 2 = 0.
\end{align*}

Vous allez maintenant poser $u = \tan\theta$ qui est bien défini compte tenu de la définition de $\theta.$

Les identités trigonométriques sont permettre de simplifier les calculs. Dans un souci de complétude, elles sont détaillées ici. Tout d’abord :

\begin{align*}
\cos^2\theta &=\frac{\cos^2\theta}{1}\\
&=\frac{\cos^2\theta}{\sin^2\theta + \cos^2\theta}\\
&=\frac{1}{u^2 + 1} 
\end{align*} 
\begin{align*}
\sin^2\theta &=1 - \cos^2\theta\\
&=1-\frac{1}{1+u^2}\\
&=\frac{u^2}{1+u^2}.
\end{align*} 

Vous déduisez :

\begin{align*}
\cos^2\theta - \sin^2\theta &=\frac{1-u^2}{1+u^2}.
\end{align*}

D’autre part :

\begin{align*}
\sin\theta\cos\theta&=\frac{\sin\theta}{\cos\theta}\cos^2\theta\\
&= u\cos^2\theta\\
&=\frac{u}{1+u^2}.
\end{align*}

Vous déduisez que :

\begin{align*}
x^2+\sqrt{3}xy+x-2=0 
&\Longleftrightarrow \left(\frac{1}{1+u^2}+\frac{\sqrt{3}u}{1+u^2}\right)X^2+\left(\frac{u^2}{1+u^2}-\frac{\sqrt{3}u}{1+u^2}\right)Y^2+\left(\sqrt{3}\frac{1-u^2}{1+u^2}-2\frac{u}{1+u^2}\right)XY+X\cos \theta - Y\sin \theta - 2 = 0\\
&\Longleftrightarrow (1+\sqrt{3}u)X^2+(u^2-\sqrt{3}u)Y^2+(\sqrt{3}(1-u^2)-2u)XY+X\cos \theta (1+u^2)- Y\sin \theta (1+u^2)- 2(1+u^2) = 0.
\end{align*}

Pour éliminer le terme croisé $XY$ vous allez chercher un réel $u$ tel que :

\begin{align*}
\sqrt{3}(1-u^2)-2u &=0\\
3(1-u^2)-2\sqrt{3}u &=0\\
3-3u^2-2\sqrt{3}u &=0\\
3u^2+2\sqrt{3}u -3&=0\\
(\sqrt{3}u+1)^2-4&=0\\
(\sqrt{3}u+1)^2-2^2&=0\\
(\sqrt{3}u-1)(\sqrt{3}u+3)&=0\\
(3u-\sqrt{3})(3u+3\sqrt{3})&=0\\
(3u-\sqrt{3})(u+\sqrt{3})&=0\\
\end{align*}

Ainsi vous choisissez $\theta = \arctan \frac{\sqrt{3}}{3}.$ Vous obtenez immédiatement $u = \tan \theta = \frac{\sqrt{3}}{3}.$

Ainsi vous avez bien $\sqrt{3}(1-u^2)-2u =0$ et le terme en $XY$ a disparu.

Note. Vous aurez remarqué que dans ce cas précis, $\theta = \frac{\pi}{6}.$ Notez qu’il est possible de poursuivre sans calculer $\theta$ explicitement, c’est ce qui va être effectué.

Eliminez le terme croisé

De part la positivité du nombre $\frac{\sqrt{3}}{3}$ vous déduisez que $\theta$ il appartient à $[0,\pi/2[$ et par conséquent $\sin\theta$ et $\cos\theta$ sont positifs.

Ainsi:

\begin{align*}
\cos\theta &= \frac{1}{\sqrt{1+u^2}}\\
 &= \frac{1}{\sqrt{1+\frac{1}{3}}}\\
&= \frac{1}{\sqrt{\frac{4}{3}}}\\
&= \frac{\sqrt{3}}{2}.
\end{align*}

Du coup:

\begin{align*}
\sin\theta &=u\cos\theta\\
 &=\frac{\sqrt{3}}{3}\times \frac{\sqrt{3}}{2}\\
&= \frac{1}{2}.
\end{align*}

Vous calculez les coefficients respectifs:

\begin{align*}
1+\sqrt{3}u &=1+\sqrt{3} \times \frac{\sqrt{3}}{3}\\
&=2.
\end{align*}
\begin{align*}
u^2-\sqrt{3}u &= \frac{1}{3}- 1\\
&=\frac{-2}{3}.
\end{align*}
\begin{align*}
1+u^2 &= 1+\frac{1}{3}\times 1\\
&=\frac{4}{3.}
\end{align*}
\begin{align*}
x^2+\sqrt{3}xy+x-2=0 
&\Longleftrightarrow 2X^2-\frac{2}{3}Y^2+\frac{\sqrt{3}}{2}\times \frac{4}{3}X- \frac{1}{2}\times \frac{4}{3}Y-\frac{8}{3} = 0\\
&\Longleftrightarrow 2X^2-\frac{2}{3}Y^2+\frac{2\sqrt{3}}{3}X- \frac{2}{3}Y-\frac{8}{3} = 0\\
&\Longleftrightarrow 6X^2-2Y^2+2\sqrt{3}X- 2Y-8 = 0\\
&\Longleftrightarrow 3X^2-Y^2+\sqrt{3}X- Y-4 = 0.
\end{align*}

Concluez

Via une rotation du repère vous avez obtenu:

\boxed{x^2+\sqrt{3}xy+x-2=0   \Longleftrightarrow 3X^2-Y^2+\sqrt{3}X- Y-4 = 0.}

Notez l’apparition du terme $Y^2.$

Prolongement

Malgré la disparition du terme $xy$ il vous faut encore effectuer autre chose, comme une opération afin d’éliminer les termes en $X$ et en $Y.$ Cela va vous permettre de déduire les éléments de la conique étudiée.

Pour en savoir davantage, reportez au contenu rédigé dans l'article 237.

235. Théorème de Weierstrass par les polynômes de Hermite

Soit $f: [-1,1]\to \R$ une fonction continue. Considérez un réel strictement positif $\varepsilon$ fixé.

Vous allez démontrer qu’il existe un polynôme $P\in\R[X]$ tel que :

\forall x\in[-1,1], \vert f(x)-P(x)\vert \leq \varepsilon.

Utilisez la continuité uniforme de la fonction $f$

Via le théorème de Heine qui a été démontré dans l'article 192, la fonction $f$ est uniformément continue sur l’intervalle $[-1,1].$ Ainsi, il existe un réel $\delta$ strictement positif tel que :

\forall (x,y)\in[-1,1]^2, \vert x-y \vert \leq \delta \implies\vert f(x)-f(y)\vert \leq \frac{\varepsilon}{2}.

La continuité de $f$ implique aussi le caractère majoré de la fonction $\vert f\vert:$

\exists M>0, \forall x\in[-1,1], \vert f(x) \vert \leq M.

Soit $n$ un entier naturel tel que $n >2+\frac{8M}{\varepsilon \delta^2}$, afin d’obtenir $\frac{4M}{n\delta^2} < \frac{\varepsilon}{2}$ et pour avoir au moins 2 points d’interpolation.

Réalisez une interpolation de la fonction $f$ aux points de Tchebychev

Pour tout $k\in\llbracket 1, n\rrbracket$ vous posez :

x_k = \cos\frac{(2k-1)\pi}{2n}.

Vous réalisez une interpolation de la fonction $f$ aux points de Tchebychev par le polynôme de Hermite qui y possède des dérivées toutes nulles en ces points.

Pour tout $k\in\llbracket 1,n \rrbracket$ vous définissez le polynôme de Hermite $H_k$ qui vaut $1$ en $x_k$ et qui possède les propriétés décrites et démontrées dans l'article 232.

Pour rappel, le polynôme $H_k$ est défini par l’expression suivante, à l’aide d’un polynôme $Q_k:$

\begin{align*}
Q_k(X)&=\prod_{\substack{1\leq i \leq n \\i\neq k}} (X-x_i)^2\\
H_k(X) &= \frac{(Q_k(x_k) - Q'(x_k)(X-x_k))Q_k(X)}{Q_k(x_k)^2}.
\end{align*}

Le polynôme $P$ suivant est un polynôme interpolateur de la fonction $f$ au points $x_1,\dots,x_n:$

P(X) = \sum_{k=1}^n f(x_k)H_k(X).

Majorez l’écart obtenu

Soit $x$ un nombre réel appartenant à l’intervalle $[-1,1].$ Utilisez d’abord la définition du polynôme $P:$

\lvert f(x)-P(x)\vert \leq \left\vert f(x) - \sum_{k=1}^n f(x_k)H_k(x) \right\vert.

Ensuite, vous utilisez le fait que la somme des $H_k$ est égale à $1$, $\sum_{k=1}^n H_k(X) = 1$, comme cela a été vu au sein du contenu écrit dans l'article 230. Vous obtenez alors en combinaison avec l’inégalité triangulaire :

\begin{align*}
\lvert f(x)-P(x)\vert &\leq \left\vert f(x) \sum_{k=1}^n H_k(x) - \sum_{k=1}^n f(x_k)H_k(x) \right\vert\\
&\leq \left\vert  \sum_{k=1}^n (f(x)-f(x_k))H_k(x)\right\vert\\
&\leq   \sum_{k=1}^n \left\vert (f(x)-f(x_k))H_k(x)\right\vert.
\end{align*}

Dans le cas d’interpolation avec les points de Tchebychev, il a été établi au sein du contenu écrit dans l'article 232 que $\forall k\in\llbracket 1, n\rrbracket, H_k(x)\in[0,1]$, d’où :

\begin{align*}
\lvert f(x)-P(x)\vert &\leq   \sum_{k=1}^n \left\vert f(x)-f(x_k)\right\vert H_k(x)\\
&\leq   \sum_{\substack{1\leq k\leq n\\\vert x-x_k\vert \leq \delta}} \left\vert f(x)-f(x_k)\right\vert H_k(x)+\sum_{\substack{1\leq k\leq n\\\vert x-x_k\vert >\delta}} \left\vert f(x)-f(x_k)\right\vert H_k(x)\\
&\leq   \sum_{\substack{1\leq k\leq n\\\vert x-x_k\vert \leq \delta}} \frac{\varepsilon}{2} H_k(x)+\sum_{\substack{1\leq k\leq n\\\vert x-x_k\vert >\delta}} (\vert f(x)\vert+\vert f(x_k)\vert) H_k(x)\\
&\leq   \frac{\varepsilon}{2}\sum_{\substack{1\leq k\leq n\\\vert x-x_k\vert \leq \delta}}  H_k(x)+\sum_{\substack{1\leq k\leq n\\\vert x-x_k\vert >\delta}} 2M H_k(x)\\
&\leq   \frac{\varepsilon}{2}\sum_{k=1}^n  H_k(x)+\sum_{\substack{1\leq k\leq n\\\vert x-x_k\vert >\delta}} 2M H_k(x)\\
&\leq   \frac{\varepsilon}{2}\sum_{k=1}^n  H_k(x)+2M \sum_{\substack{1\leq k\leq n\\\vert x-x_k\vert >\delta}}  H_k(x)\\
&\leq   \frac{\varepsilon}{2}+2M \sum_{\substack{1\leq k\leq n\\\vert x-x_k\vert >\delta}}  H_k(x).
\end{align*}

Or, dans l'article 232 il a été vu que :

\forall k\in\llbracket 1,n\rrbracket, x_k\neq x\implies\lvert H_k(x)\vert \leq \frac{2}{n^2(x-x_k)^2}.

Vous déduisez finalement :

\begin{align*}
\lvert f(x)-P(x)\vert  &\leq   \frac{\varepsilon}{2}+2M \sum_{\substack{1\leq k\leq n\\\vert x-x_k\vert >\delta}}  \frac{2}{n^2(x-x_k)^2}\\
&\leq   \frac{\varepsilon}{2}+2M \sum_{\substack{1\leq k\leq n\\\vert x-x_k\vert >\delta}}  \frac{2}{n^2\delta^2}\\
&\leq   \frac{\varepsilon}{2}+2M \sum_{k=1}^n  \frac{2}{n^2\delta^2}\\
&\leq   \frac{\varepsilon}{2}+\frac{4M}{n\delta^2}\\
&\leq   \frac{\varepsilon}{2}+\frac{\varepsilon}{2}\\
&\leq \varepsilon.
\end{align*}

Concluez

Il a été démontré dans cet article que toute fonction continue sur l’intervalle $[-1,1]$ est limite uniforme de polynômes.

Plus précisément, pour toute fonction continue $f$ sur l’intervalle $[-1,1]:$

\boxed{ \forall \varepsilon>0, \exists P\in\R[X], \forall x\in[-1,1], \vert f(x)-P(x)\vert \leq \varepsilon.}

Visualisations graphiques

La fonction $x\mapsto \frac{1}{1+25x^2}$ est représentée ci-dessous avec plusieurs approximations uniformes basées sur les polynômes de Hermite calculés précédemment.

234. Interpolation de Hermite par différences divisées

Dans cet article vous supposerez que les différences divisées généralisées permettent de calculer un polynôme interpolateur de Hermite dans le cas où le nombre de conditions est égal à $3.$

Rappelez-vous que les différences divisées généralisées d’une fonction $f$, dans le cas où les paramètres sont identiques, fournissent : $f[x] = f(x)$, $f[x,x] = f'(x)$ et $f[x,x,x] = \frac{f »(x)}{2}.$ Dans le cas où $x$ serait répété $n$ fois, $n$ désignant un nombre entier supérieur ou égal à $2$, vous auriez $f[x,\dots,x] = \frac{f^{(n)}(x)}{n!}.$

Ce rapprochement est à mettre en rapport avec la formule de Taylor.

Cela étant posé, le but de cet article est de construire un polynôme $P$ vérifiant les conditions suivantes :

\left\{\begin{align*}
P(1) &= 2\\
P'(1) &= 7\\
P''(1)&=1\\
P(2)&=3.
\end{align*}\right.

Scindez les conditions en deux parties

Supposez que vous avez déjà à votre disposition un polynôme $R$ tel que :

\left\{\begin{align*}
R(1) &= 2\\
R'(1) &= 7\\
R''(1)&=1.
\end{align*}\right.

Supposez de plus que vous avez déjà à votre disposition un polynôme $S$ tel que :

\left\{\begin{align*}
S(1) &= 2\\
S'(1) &= 7\\
S(2)&=3.
\end{align*}\right.

Alors, dans l’esprit du contenu se trouvant dans l'article 233, il suffit de poser :

P(X)=\frac{R(X)(X-2)-S(X)(X-1)}{1-2}.

Le polynôme $P$ va alors convenir.

Vérifiez les premières conditions sans dérivation :

\begin{align*}
P(1)&=\frac{R(1)(1-2)-S(1)(1-1)}{1-2}\\
&=\frac{R(1)(1-2)}{1-2}\\
&=R(1)\\
&=2.
\end{align*}

Ensuite :

\begin{align*}
P(2)&=\frac{R(2)(2-2)-S(2)(2-1)}{1-2}\\
&=\frac{-S(2)(2-1)}{1-2}\\
&=S(2)\\
&=3.
\end{align*}

Pour vérifier la condition avec la dérivée première, vous dérivez une fois puis vous substituez :

\begin{align*}
P'(X) &= \frac{R'(X)(X-2)+R(X)-S'(X)(X-1)-S(X)}{1-2}\\
P'(1) &= \frac{R'(1)(1-2)+R(1)-S(1)}{1-2}\\
 &= \frac{R'(1)(1-2)}{1-2}\\
 &= R'(1)\\
 &= 7.
\end{align*}

Pour la condition avec la dérivée seconde, vous utilisez la formule de Leibniz sur les deux produits.

\begin{align*}
P''(X) &=\frac{R''(X)(X-2)+2R'(X)-S''(X)(X-1)-2S'(X)}{1-2}\\
P''(1) &=\frac{R''(1)(1-2)+2R'(1)-2S'(1)}{1-2}\\
 &=\frac{R''(1)(1-2)}{1-2}\\
 &=R''(1)\\
&=1.
\end{align*}

Reconstruisez le polynôme $P$

Il s’agit de comprendre pourquoi le polynôme $P$ se calcule lui aussi avec des différences divisées.

Le polynôme $R$ s’obtient directement par la formule de Taylor :

\begin{align*}
R(X) = R(1)+R'(1)(X-1)+\frac{R''(1)}{2}(X-1)^2.
\end{align*}

Vous déduisez immédiatement que :

\begin{align*}
R(X) = P[1]+P[1,1](X-1)+P[1,1,1](X-1)^2.
\end{align*}

Admettant la possibilité de construire les polynômes de Hermite par différence divisée pour 3 conditions, le polynôme $S$ est égal à :

\begin{align*}
S(X) = S[1]+S[1,1](X-1)+S[1,1,2](X-1)^2.
\end{align*}

Le tableau des différences divisées de $S$ est alors le suivant, compte tenu des égalités $S(1)=2, S(2)=3$ et $S'(1)=7.$

\begin{array}{c|cc}
1 & 2 & 7 & \times \\
1 & 2 & \times \\
2 & 3
\end{array}

Vous le complétez au fur et à mesure :

\begin{array}{c|cc}
1 & 2 & 7 & \times \\
1 & 2 & \frac{3-2}{2-1}\\
2 & 3
\end{array}
\begin{array}{c|cc}
1 & 2 & 7 & \times \\
1 & 2 & 1\\
2 & 3
\end{array}
\begin{array}{c|cc}
1 & 2 & 7 & \frac{1-7}{2-1}\\
1 & 2 & 1\\
2 & 3
\end{array}
\begin{array}{c|cc}
1 & 2 & 7 & -6\\
1 & 2 & 1\\
2 & 3
\end{array}

Vous déduisez de ces calculs que $S(X) = 2 + 7(X-1)-6(X-1)^2.$

De l’expression de $P$ suivante :

P(X)=\frac{R(X)(X-2)-S(X)(X-1)}{1-2}.

Vous déduisez que le coefficient dominant de $P$ n’est autre que :

\frac{P[1,1,1] - S[1,1,2]}{1-2}

Le tableau des différences divisées de $P$ fournit :

\begin{array}{c|ccc}
1 & 2 & 7 & 1/2 & P[1,1,1,2]\\
1 & 2 & 7 & -6 &  \\
1 & 2  & 1 & \\
2 & 3
\end{array}

Il s’agit du tableau des différences divisées de $S$ avec une ligne en plus, si bien que :

\begin{align*}
\frac{P[1,1,1] - S[1,1,2]}{1-2} &= \frac{ S[1,1,2] - P[1,1,1]}{2-1}\\
&= \frac{-6-1/2}{2-1}\\
&=\frac{-13}{2}\\
&=P[1,1,1,2].
\end{align*}

Maintenant que vous savez que le coefficient dominant de $P$ est $P[1,1,1,2] = \frac{-13}{2}$ il reste à trouver une expression plus agréable pour $P$ que celle donnée avec $R$ et $S.$ Les polynômes $P$ et $R$ sont égaux en $1$, ainsi que leurs dérivées première et seconde.

Posez alors :

\begin{align*}
Q(X) &= P(X)- P[1]- P[1,1](X-1)- P[1,1,1](X-1)^2\\
&= P(X)-R(X).
\end{align*}

Comme $Q(1) = Q'(1) = Q »(1) = 0$ le nombre $1$ est racine triple de $Q$ donc $Q$ est factorisable par $(X-1)^3.$

Comme $R\in\R_2[X]$, l’expression $P(X)=\frac{R(X)(X-2)-S(X)(X-1)}{1-2}$ montre que $P\in\R_3[X].$ Pour des raisons de degré, $Q$ est donc égal à une constante $C$ multipliée par $(X-1)^3.$ Comme $R\in\R_2[X]$ le coefficient de $X^3$ est le même dans $Q$ et dans $P$. Donc la constante $C$ est nécessairement le coefficient dominant de $P.$

De cette analyse, vous déduisez :

P(X)=P[1]+P[1,1](X-1)+P[1,1,1](X-1)^2+P[1,1,1,2](X-1)^3.

Concrètement, le polynôme $P$ cherché est :

P(X)=2+7(X-1)+\frac{1}{2}(X-1)^2-\frac{13}{2}(X-1)^3.

En développant, il est égal à :

\boxed{P(X) = -\frac{13 X^3}{2}+20 X^2-\frac{27 X}{2}+2.}

La représentation graphique de ce polynôme est présentée ci-dessous.

17/03/2022 - Polynome de hermite avec 4 conditions 64d324db00af295894b0a855e51c181b71fe36f3

Prolongement

A vous de trouver l’expression développée d’un polynôme de Hermite $H$ tel que

\left\{\begin{align*}
H(1) &= 1\\
H'(1) &= 0\\
H(2) &= 3\\
H'(2) &= 0\\
H(3) &=2\\
H'(3)&=0.
\end{align*}\right.

Ci-dessous, vous trouverez la représentation graphique d’un tel polynôme.

17/03/2022 - Polynome de hermite avec 6 conditions aef3160f30b48f69c384972a10b9e12bb6455325

233. Construction des formules de différences divisées

Dans cet article, $x_1, x_2$ et $x_3$ sont trois nombres réels deux à deux distincts.

Vous vous fixez trois réels quelconques notés $y_1, y_2$ et $y_3$, qui peuvent être identiques ou non.

Vous cherchez à construire un polynôme $P$ tel que $P(x_1)=y_1$, $P(x_2)=y_2$ et $P(x_3) = y_3.$

Vous allez y parvenir en suivant les sections suivantes :

  • La première étape va consister à construire un polynôme $R$ tel que $R(x_1)=y_1$ et $R(x_2)=y_2;$
  • La deuxième étape consiste à construire un polynôme $S$ tel que $S(x_2)=y_2$ et $S(x_3)=y_3;$
  • La troisième étape va donner le polynôme $P$ à partir des polynômes $R$ et $S$ en leur ajoutant une racine ;
  • La dernière étape permet de trouver une expression plus agréable et plus systématique.

Etape 1 : construisez un polynôme $R$ tel que $R(x_1)=y_1$ et $R(x_2)=y_2$

Analyse

En cherchant ce polynôme sous la forme $R(X)=aX+b$, vous êtes amené à résoudre le système :

\left\{\begin{align*}
ax_1+b&=y_1\\
ax_2+b&=y_2.
\end{align*}\right.

Par soustraction, il vient $ax_1-ax_2 = y_1-y_2$ d’où $a(x_1-x_2)=y_1-y_2$ et comme $x_1\neq x_2$ il est possible de diviser par $x_1-x_2$ ce qui fournit $a = \frac{y_1-y_2}{x_1-x_2}.$

Ce nombre est si important qu’il sera noté dans la suite $P[x_1,x_2].$ Ainsi :

 \boxed{P[x_1,x_2] = \frac{y_1-y_2}{x_1-x_2}.}

$P[x_1,x_2] = \frac{y_1-y_2}{x_1-x_2}.$

Alors il vient $b = y_1 – ax_1 = y_1-P[x_1,x_2] x_1.$

Ainsi un bon candidat pour $R(X)$ est $R(X) = aX+b = P[x_1,x_2] X + y_1 – P[x_1,x_2] x_1 = y_1 + P[x_1,x_2](X-x_1).$

Par souci de cohérence avec le résultat précédent vous notez :

 \boxed{P[x_1] = y_1.}

Synthèse

Posez :

R(X) = P[x_1]+P[x_1,x_2](X-x_1).

En substituant $x_1$, il vient $R(x_1) = P[x_1] = y_1.$

En substituant $x_2$, vous obtenez :

\begin{align*} 
R(x_2)  &= y_1 + \frac{y_1-y_2}{x_1-x_2}\times (x_2-x_1)\\
&= y_1-(y_1-y_2)\\
&=y_2.
\end{align*}

Par conséquent le polynôme $R$ convient bien au problème posé.

Etape 2 : construisez un polynôme $S$ tel que $S(x_2)=y_2$ et $S(x_3)=y_3$

En suivant exactement la même démarche que précédemment, vous posez $P[x_2] = y_2$ et $P[x_2,x_3] = \frac{y_2-y_3}{x_2-x_3}.$

Ainsi, le polynôme $S$ suivant va convenir :

S(X) = P[x_2]+P[x_2,x_3](X-x_2).

Etape 3 : construisez le polynôme $P$ recherché

Vous avez un polynôme $R$ qui vérifie :

\begin{align*}
R(x_1)&=y_1\\
R(x_2)&=y_2.
\end{align*}

Vous avez un polynôme $S$ qui vérifie :

\begin{align*}
S(x_2)&=y_2\\
S(x_3)&=y_3.
\end{align*}

Posez alors $\mathscr{R}(X) = R(X)(X-x_3)$ et $\mathscr{S}(X) = S(X)(X-x_1).$ Du coup :

\begin{align*}
\mathscr{R}(x_1)&=y_1(x_1-x_3)\\
\mathscr{R}(x_2)&=y_2(x_2-x_3)\\
\mathscr{R}(x_3)&=0\\
\mathscr{S}(x_1)&=0\\
\mathscr{S}(x_2)&=y_2(x_2-x_1)\\
\mathscr{S}(x_3)&=y_3(x_3-x_1)\\
\end{align*}

Par soustraction, vous déduisez :

\begin{align*}
\mathscr{(R-S)}(x_1)&=y_1(x_1-x_3)\\
\mathscr{(R-S)}(x_2)&=y_2(x_2-x_3)-y_2(x_2-x_1)\\
\mathscr{(R-S)}(x_3)&=-y_3(x_3-x_1).
\end{align*}

Vous n’obtenez que des multiples de $x_1-x_3$ :

\begin{align*}
\mathscr{(R-S)}(x_1)&=y_1(x_1-x_3)\\
\mathscr{(R-S)}(x_2)&=y_2(x_1-x_3)\\
\mathscr{(R-S)}(x_3)&=y_3(x_1-x_3).
\end{align*}

Comme $x_1\neq x_3$, en divisant, vous obtenez le polynôme cherché :

P(X) = \frac{R(X)(X-x_3)-S(X)(X-x_1)}{x_1-x_3}.

Le coefficient dominant du polynôme $P$ est alors égal à $\frac{P[x_1,x_2]-P[x_2,x_3]}{x_1-x_3}.$

De part l’importance de ce résultat, appelé différence divisée, vous posez :

\boxed{P[x_1,x_2,x_3] = \frac{P[x_1,x_2]-P[x_2,x_3]}{x_1-x_3}.}

Etape 4 : déterminez une expression plus agréable du polynôme $P$

Comme $R$ et $S$ sont de degré $1$, le polynôme $P$ appartient à $\R_2[X].$

Par soustraction, le polynôme $Q = P-R$ appartient aussi à $\R_2[X].$

Or $Q(x_1) = P(x_1)-R(x_1) = y_1-y_1 = 0$ et $Q(x_2) = P(x_2)-R(x_2) = y_2-y_2 = 0$ donc $Q$ admet deux racines distinctes $x_1$ et $x_2$, donc il est factorisable par $(X-x_1)(X-x_2).$ Il existe donc un polynôme $T$ tel que $Q(X)=(X-x_1)(X-x_2)T(X).$ Pour une question de degré, $T$ est constant, sinon $Q$ n’appartiendrait pas à $R_2[X].$

Notez $C = T(X).$ En développant $C(X-x_1)(X-x_2)$ vous déduisez que $C$ est le coefficient dominant de $Q.$

Or il a été vu que ce dernier est $P[x_1,x_2,x_3].$

Donc $Q(X) = P[x_1,x_2,x_3](X-x_1)(X-x_2).$

Comme $P = R+Q$ vous déduisez :

\boxed{P(X) = P[x_1]+P[x_1,x_2](X-x_1)+P[x_1,x_2,x_3](X-x_1)(X-x_2).}

Cette expression est attribuée à Newton.

Prolongement

Soit $n$ un entier naturel supérieur ou égal à $2.$ Considérez $n$ réels deux à deux distincts notés $x_1,\dots,x_n$. Soient $n$ réels $y_1,\dots y_n$ quelconques. Pourriez-vous généraliser en justifiant non seulement de l’existence mais d’un procédé de calcul permettant de déterminer l’expression d’un polynôme $P$ de degré $n+1$ vérifiant $\forall j\in\llbracket 1, n\rrbracket, P(x_j)=y_j$ ?

231. Factorisez les polynômes de Tchebychev

Dans le cadre du prolongement du contenu écrit dans l'article 227 vous avez établi qu’il existe une suite de polynômes notée $(T_n)_{n\in\N}$ vérifiant les propriétés suivantes :

T_0(X)=1\\
T_1(X)=X\\
\forall n\in\N, T_{n+2}(X)=2XT_{n+1}(X)-T_n(X).

Pour tout $n\in\N^{*}$ le polynôme $T_n$ est de degré $n$ et son coefficient dominant est $2^{n-1}.$

De plus, pour tout réel $x$ et pour tout entier naturel $n$, $T_n(\cos x) = \cos nx.$

Cette dernière propriété va vous permettre de factoriser complètement les polynômes $T_n$ où $n\geq 1.$

Calculez les racines des polynômes de Tchebychev

Soit $n$ un entier naturel fixé et non nul.

Vous cherchez un nombre $y\in\R$ tel que $T_n(y)=0.$

Compte tenu de la propriété précitée, il semble naturel de chercher $y$ sous la forme $y=\cos x.$

Dès lors, $y$ sera racine de $T_n$ dès que $\cos nx$ sera nul.

Or, vous avez la série d’équivalences suivante :

\begin{align*}
\forall x\in\R, \cos nx = 0 &\Longleftrightarrow \exists k\in\Z, nx = \frac{\pi}{2}+k\pi\\
&\Longleftrightarrow \exists k\in\Z, x = \frac{\pi}{2n}+\frac{k\pi}{n}\\
&\Longleftrightarrow \exists k\in\Z, x = \frac{(2k+1)\pi}{2n}\\
&\Longleftrightarrow \exists k\in\Z, x = \frac{(2k-1)\pi}{2n}.
\end{align*}

Maintenant, pour tout $k\in\llbracket 1, n\rrbracket$ posez $x_k = \frac{(2k-1)\pi}{2n}.$

Alors la suite $(x_k)_{1\leq k \leq n}$ est strictement croissante.

Soit $k\in\llbracket 1, n\rrbracket.$ Comme $\frac{\pi}{2n} \leq x_k \leq \frac{(2n-1)\pi}{2n}$ vous déduisez que $\forall k\in \llbracket 1, n\rrbracket x_k\in]0,\pi[.$

Pour tout $k\in\llbracket 1, n\rrbracket$ posez $y_k = \cos x_k = \cos\left(\frac{(2k-1)\pi}{2n}\right).$

Comme la fonction cosinus est strictement monotone sur l’intervalle $]0,\pi[$ vous déduisez que les $n$ réels $y_1, \dots, y_n$ sont deux à deux distincts. Ils appartiennent tous à l’intervalle ouvert $]-1, 1[$ et ils vérifient ceci: $\forall k\in\llbracket 1, n\rrbracket, T_n(y_k) = 0.$

Finalisez la factorisation des polynômes de Tchebychev

Soit $n$ un entier naturel fixé et non nul.

Le polynôme $T_n$ étant de degré $n$ et possédant $n$ racines deux à deux distinctes $y_1, \dots, y_n$ vous déduisez l’existence d’une constante $C\in\R $ telle que:

 T_n(X) = C \prod_{i=1}^n (X-y_i).

En développant, il apparaît que $C$ est le coefficient dominant de $T_n$ et donc $C = 2^{n-1}$ compte tenu du contenu écrit dans l'article 227 .

En définitive vous avez obtenu la factorisation complète :

\boxed{\forall n\in\N^{*}, T_n(X) = 2^{n-1}\prod_{k=1}^n\left(X-\cos \frac{(2k-1)\pi}{2n}\right).}

De ce qui précède, vous avez montré que les polynômes de Tchebychev sont scindés, à racines simples appartenant toutes à l’intervalle ouvert $]-1,1[.$

Visualisez les racines

12/03/2022 - Polynome de tchebychev numero 5 et ses racines 03c2505a4050c4a540f679440baa629927982a91
Tracé du polynôme de Tchebychev $T_5(X) = 16 X^5-20 X^3+5 X$ et de ses 5 racines
12/03/2022 - Polynome de tchebychev numero 10 et ses racines bccdc3edf69801607cfd01b42f6f3a5fe452e7a4
Tracé du polynôme de Tchebychev $T_{10}(X) = 512 X^{10}-1280 X^8+1120 X^6-400 X^4+50 X^2-1$ et de ses 10 racines
12/03/2022 - Polynome de tchebychev numero 16 et ses racines 7a05435b670895df91cfb334f94db526178a7495
Tracé du polynôme de Tchebychev $T_{16}(X) = 32768 X^{16}-131072 X^{14}+212992 X^{12}-180224 X^{10}+84480 X^8-21504 X^6+2688 X^4-128 X^2+1$ et de ses 16 racines

Les racines de $T_n$ ne sont pas réparties de façon uniforme. Plus $n$ augmente, plus le nombre de racines proches des bords de l’intervalle $[-1,1]$ est élevé.

230. La somme des polynômes de Hermite est égale à 1

Dans cet article, $n$ désigne un nombre entier supérieur ou égal à $2.$

Vous considérez $n$ réels deux à deux distincts que vous notez $x_1,\dots x_n.$

Pour tout $k\in\llbracket 1, n\rrbracket$ vous définissez les polynômes $Q_k$ et $H_k$ comme suit :

\begin{align*}
Q_k(X) &= \prod_{\substack{1\leq i \leq n\\ i\neq k}} (X-x_i)^2\\
H_k(X) &= \frac{(Q_k(x_k)-Q'(x_k)(X-x_k))Q_k(X)}{Q_k(x_k)^2}.
\end{align*}

Pour tout $k\in\llbracket 1, n\rrbracket$, le polynôme $H_k$ est appelé $k$ième polynôme de Hermite.
Il vérifie les propriétés suivantes, d’après le contenu que vous trouverez dans l'article 229 :

H_k(x_k) = 1\\
\forall i\in\llbracket 1, n\rrbracket, i\neq k \implies H_k(x_i) = 0\\
\forall i\in\llbracket 1, n\rrbracket, H'_k(x_i) = 0.

Pour démontrer que $\sum_{k=1}^n H_k(X)=1$ vous considérez le polynôme suivant :

P(X) = \sum_{k=1}^n H_k(X)-1.

Soit $m\in\llbracket 1, n \rrbracket.$ Notez tout d’abord que :

\begin{align*}
P(x_m) &= \sum_{k=1}^n H_k(x_m)-1\\
&=H_m(x_m) + \sum_{k\neq m} H_k(x_m)-1\\
&=H_m(x_m) + \sum_{k\neq m} 0-1\\
&=H_m(x_m) -1\\
&=1 -1\\
&=0.
\end{align*}

Ensuite :

\begin{align*}
P'(x_m) &= \sum_{k=1}^n H'_k(x_m)\\
&=\sum_{k=1}^n 0\\
&=0.
\end{align*}

Ainsi $\boxed{\forall m\in\llbracket 1,n\rrbracket, P(x_m)=P'(x_m)=0.}$

Note. Pour tout $m\in\llbracket 1,n\rrbracket$ le réel $x_m$ est racine double de $P.$ Ainsi $P$ admet $n$ racines doubles deux à deux distinctes. Cela va entraîner la factorisation de $P$ par le polynôme $\prod_{i=1}^n (X-x_i)^2.$ Par souci de complétude, la preuve de cette factorisation va être détaillée en utilisant une récurrence limitée.

Démontrez que $P$ est factorisable par $\prod_{i=1}^n (X-x_i)^2$

Pour tout $m\in\llbracket 1,n\rrbracket$ notez $\mathscr{P}(m)$ la propriété : « $\exists R\in\R[X], P(X) = R(X) \prod_{i=1}^m (X-x_i)^2$. »

Initialisation. Pour $m=1.$ Comme $x_1$ est racine de $P$, il existe un polynôme $S\in\R[X]$ tel que $P(X) = (X-x_1)S(X).$

En dérivant, il vient :

\begin{align*}
P'(X)&=S(X)+(X-x_1)S'(X)\\
P'(x_1)&=S(x_1)\\
0 &= S(x_1).
\end{align*}

Comme $x_1$ est racine de $S$, il existe un polynôme $R\in\R[X]$ tel que $S(X) = (X-x_1)R(X).$

Du coup, $P(X)=(X-x_1)^2R(X)$ ce qui prouve $\mathscr{P}(1).$

Hérédité. Soit $m\in\llbracket 1, n-1\rrbracket$ tel que $\mathscr{P}(m)$ soit vérifiée.

Il existe un polynôme $S\in\R[X]$ tel que $P(X) = S(X) \prod_{i=1}^m (X-x_i)^2.$ Alors :

\begin{align*}
P(x_{m+1}) &=S(x_{m+1}) \prod_{i=1}^m (x_{m+1}-x_i)^2\\
0 &=S(x_{m+1}) \prod_{i=1}^m (x_{m+1}-x_i)^2\\
0 &= S(x_{m+1}).
\end{align*}

Le polynôme $S$ est donc factorisable par $X-x_{m+1}.$

Il existe un polynôme $T\in\R[X]$ tel que $S(X)=(X-x_{m+1})T(X).$

Alors $P(X) = T(X) (X-x_{m+1})\prod_{i=1}^m (X-x_i)^2.$ En dérivant :

\begin{align*}
P'(X) &= T'(X) (X-x_{m+1})\prod_{i=1}^m (X-x_i)^2 + T(X)\prod_{i=1}^m (X-x_i)^2+T(X) (X-x_{m+1})\left[\sum_{i=1}^m 2(X-x_i)\prod_{k\neq i}(X-x_k)^2 \right] \\
P'(x_{m+1}) &=T(x_{m+1})\prod_{i=1}^m (x_{m+1}-x_i)^2\\
0 &=T(x_{m+1})\prod_{i=1}^m (x_{m+1}-x_i)^2\\
0 &=T(x_{m+1}).
\end{align*}

Il existe donc un polynôme $R\in\R[X]$ tel que $T(X) = (X-x_{m+1})R(X).$ Du coup $S(X)=(X-x_{m+1})^2 R(X).$

Donc $P(X) = R(X) (X-x_{m+1})^2 \prod_{i=1}^m (X-x_i)^2 = R(X)\prod_{i=1}^{m+1} (X-x_i)^2.$

Ainsi $\mathscr{P}(m+1)$ est vérifiée.

Par récurrence, vous avez montré que $\mathscr{P}(n)$ est vérifiée.

Il existe donc un polynôme $R\in\R[X]$ tel que $P(X) = R(X) \prod_{i=1}^n (X-x_i)^2.$

Déduisez-en que $P$ est nul

Soit $k\in\llbracket 1,n \rrbracket.$ Le polynôme $Q_k$ appartient à $\R_{2n-2}[X]$ et par produit, $H_k$ appartient à $\R_{2n-1}[X].$ Le polynôme constant $1$ appartient à $\R_{2n-1}[X]$ donc par somme et par différence $P\in\R_{2n-1}[X].$

Comme le degré du polynôme $\prod_{i=1}^n (X-x_i)^2$ est égal à $2n$, si $R$ n’était pas nul, par produit, le degré de $P$ serait supérieur ou égal à $2n.$ Cela contredirait le fait que $P\in\R_{2n-1}[X].$

Donc $R = 0$ et par suite le résultat suivant est établi :

\boxed{\sum_{k=1}^n H_k(X) = 1.}

229. Polynômes de Hermite avec des dérivées nulles en certains points (2/2)

Soit $n$ un entier supérieur ou égal à $2$ fixé dans tout cet article.

Soit $(x_i)_{1\leq i \leq n}$ une famille de $n$ réels deux à deux distincts.

L’objectif de cet article est de vous démontrer qu’il est possible de construire petit à petit un polynôme, dit de Hermite et noté $H$, tel que :

\left\{\begin{align*}
H(x_1)&=1 \\
\forall i\in\llbracket 2, n\rrbracket, H(x_i)&=0\\
\forall i\in\llbracket 1, n\rrbracket, H'(x_i)&=0.
\end{align*}\right.

Construire directement un tel polynôme n’est a priori pas chose aisée, parce qu’il vous faut rendre la dérivée nulle en $x_1$, alors que la valeur prise en $x_1$ par le polynôme n’est pas nulle. Cela est possible si vous utilisez les formules de Taylor et l’identité remarquable $a^2-b^2$ mais c’est une autre histoire.

Pour parvenir à ce but, vous allez d’abord résoudre une autre question légèrement différente, qui sera très utile pour conclure. Il semble dans cette approche pertinent de partir d’un polynôme dont on connaît bien les racines.

Résolvez d’abord un problème différent

Dans cette section, vous allez construire un polynôme $S$ tel que :

\left\{\begin{align*}
S'(x_1)&=1 \\
\forall i\in\llbracket 1, n\rrbracket, S(x_i)&=0\\
\forall i\in\llbracket 2, n\rrbracket, S'(x_i)&=0.
\end{align*}\right.

Il s’agit de trouver un polynôme $S$ qui admet chaque $x_i$, quand $i$ varie de $2$ à $n$, comme racine double.

D’autre part, $x_1$ est racine simple du polynôme $S$ puisque $S(x_1)=0$ alors que $S'(x_1)\neq 0.$

Il semble commode de considérer le polynôme suivant :

T(X)=(X-x_1)\prod_{i=2}^n(X-x_i)^2.

Le calcul du polynôme dérivé $T’$ fournit :

T'(X)=\prod_{i=2}^n(X-x_i)^2 + \sum_{i=2}^n 2(X-x_1)(X-x_i)\prod_{k\neq i} (X-x_k)^2.

L’évaluation en $x_1$ fournit :

T'(x_1) = \prod_{i=2}^n (x_1-x_i)^2.

Pour tout $i\in\llbracket 2, n\rrbracket, x_1\neq x_i$ donc $T'(x_i)\neq 0.$

En définitive, un polynôme satisfaisant est trouvé dès que vous posez :

\boxed{S(X) = \frac{(X-x_1)\displaystyle\prod_{i=2}^n(X-x_i)^2}{ \displaystyle\prod_{i=2}^n (x_1-x_i)^2}.}

Revenez au problème initial

Pour tout $i\in\llbracket 2, n\rrbracket$ le nombre $x_i$ est racine double du polynôme $H.$

Vous allez donc poser $\boxed{Q(X) = \prod_{i=2}^n(X-x_i)^2}$ ce qui construit un polynôme $Q$ tel que :

\left\{\begin{align*}
\forall i\in\llbracket 2, n\rrbracket, Q(x_i)&=0\\
\forall i\in\llbracket 2, n\rrbracket, Q'(x_i)&=0.
\end{align*}\right.

Vous allez maintenant poser $R(X) = Q(X)-Q'(x_1)S(X).$ Alors $R'(X) = Q'(X)-Q'(x_1)S'(X).$

Pour tout $i\in\llbracket 2, n\rrbracket, R(x_i) = Q(x_i)-Q'(x_1)S(x_i) = 0-Q'(x_1)\times 0 = 0.$

D’autre part, $R'(x_1) = Q'(x_1)-Q'(x_1)S'(x_1)=Q'(x_1)-Q'(x_1)\times 1 = 0.$

Le polynôme $R$ satisfait les conditions suivantes :

\left\{\begin{align*}
\forall i\in\llbracket 2, n\rrbracket, R(x_i)&=0\\
\forall i\in\llbracket 1, n\rrbracket, R'(x_i)&=0.
\end{align*}\right.

Maintenant, $R(x_1) = Q(x_1)-Q'(x_1)S(x_1) = Q(x_1)-Q'(x_1)\times 0 = Q(x_1)=\prod_{i=2}^n (x_1-x_i)^2$ qui est non nul.

Le polynôme de Hermite cherché est :

H(X)=\frac{Q(X)-Q'(x_1)S(X)}{Q(x_1)}.

Avec les notations prises sur le polynôme $Q:$

\begin{align*}
S(X) &= \frac{(X-x_1)\prod_{i=2}^n(X-x_i)^2}{ \prod_{i=2}^n (x_1-x_i)^2}\\
 &= \frac{(X-x_1)Q(X)}{Q(x_1)}.
\end{align*}

Du coup vous déduisez :

\begin{align*}
H(X)&=\frac{Q(X)-Q'(x_1)S(X)}{Q(x_1)} \\
H(X)&=\frac{Q(X)-Q'(x_1)\frac{(X-x_1)Q(X)}{Q(x_1)}}{Q(x_1)} \\
H(X) &=\frac{Q(X)Q(x_1)-Q'(x_1)(X-x_1)Q(X)}{Q(x_1)^2} \\
\end{align*}

et finalement :

\boxed{H(X) = \frac{\left[Q(x_1)-Q'(x_1)(X-x_1)\right]Q(X)}{Q(x_1)^2}.}

Vérifiez le résultat construit

En dérivant le polynôme $H$ il vient :

\begin{align*}
H'(X) &= \frac{-Q'(x_1)Q(X)+\left[Q(x_1)-Q'(x_1)(X-x_1)\right]Q'(X)}{Q(x_1)^2}\\
H'(x_1) &=\frac{-Q'(x_1)Q(x_1)+Q(x_1)Q'(x_1)}{Q(x_1)^2}\\
H'(x_1) &=0.
\end{align*}

D’autre part, $H(x_1) = \frac{Q(x_1)^2}{Q(x_1)^2} = 1.$

Pour tout $i\in\llbracket 2, n\rrbracket$ vous observez que $x_i$ est racine double de $Q$, donc de $H$ et donc $H(x_i)=H'(x_i)=0.$

Illustrez ce résultat en visualisant trois polynômes de Hermite

Choisissez $x_1 = 1,$ $x_2=2$ et $x_3=3$ comme points d’interpolation.

Vous posez $Q_1(X) = (X-2)^2(X-3)^2$ et associez le polynôme de Hermite $H_1(X) = \frac{(Q_1(1)-Q’_1(1)(X-1))Q_1(X)}{Q_1(1)^2}.$ Ce polynôme vérifie $H_1(1)=1$ puis $H_1(2)=H_1(3)=0.$ Mieux qu’un polynôme de Lagrange, il vérifie $H’_1(1)=H’_1(2)=H’_1(3)=0$, c’est un polynôme interpolateur avec des dérivées nulles aux points d’interpolation.

11/03/2022 - Polynome de hermite 1 sur 3
Représentation graphique du polynôme de Hermite $H_1$ défini par $H_1(X) = -18 + 57 X – \frac{127}{2} X^2 + \frac{131}{4} X^3 – 8 X^4 + \frac{3}{4} X^5$

Vous posez $Q_2(X) = (X-1)^2(X-3)^2$ et associez le polynôme de Hermite $H_2(X) = \frac{(Q_2(2)-Q’_2(2)(X-2))Q_2(X)}{Q_2(2)^2}.$
Ce polynôme vérifie $H_2(2)=1$ puis $H_2(1)=H_2(3)=0$ avec $H’_2(1)=H’_2(2)=H’_2(3)=0.$

11/03/2022 - Polynome de hermite 2 sur 3
Représentation graphique du polynôme de Hermite $H_2$ défini par $H_2(X) = 9 – 24 X + 22 X^2 – 8 X^3 + X^4$

Vous posez $Q_3(X) = (X-1)^2(X-2)^2$ et associez le polynôme de Hermite $H_3(X) = \frac{(Q_3(3)-Q’_3(3)(X-3))Q_3(X)}{Q_3(3)^2}.$
Ce polynôme vérifie $H_3(3)=1$ puis $H_3(1)=H_3(2)=0$ avec $H’_3(1)=H’_3(2)=H’_3(3)=0.$

11/03/2022 - Polynome de hermite 3 sur 3
Représentation graphique du polynôme de Hermite $H_3$ défini par $H_3(X) = 10 – 33 X + \frac{83}{2} X^2 – \frac{99}{4} X^3 + 7 X^4 – \frac{3}{4} X^5$

Prolongement

Vous pouvez constater l’égalité suivante : $H_1(X)+H_2(X)+H_3(X) =1$ à partir des trois polynômes de Hermite explicités ci-dessus.

Or, il se trouve qu’un tel résultat est général. Il servira à démontrer que toute fonction continue sur un segment est limite uniforme de polynômes.

Pour en savoir davantage, n’hésitez pas à consulter le contenu rédigé dans l'article 230.

228. Polynômes de Hermite avec des dérivées nulles en certains points (1/2)

ll s’agit d’une amélioration des polynômes de Lagrange qui ont été utilisés par Fejér au début du XXème siècle.

Pour fixer les idées, fixez trois réels $x_1$, $x_2$ et $x_3$ deux à deux distincts.

Traitez le cas avec trois points

Vous cherchez un polynôme $H$ tel que :

\left\{\begin{align*}
H(x_1) &= 0\\
H(x_2) &= 0\\
H(x_3) &= 0\\
H'(x_1) &= 1\\
H'(x_2) &= 0\\
H'(x_3) &= 0.
\end{align*}\right.

Les réels $x_2$ et $x_3$ sont des racines doubles de $H$ alors que $x_1$ est racine simple.

Il convient de poser $P(X) = (X-x_1)(X-x_2)^2(X-x_3)^2.$

En dérivant il vient :

P'(X) = (X-x_2)^2(X-x_3)^2+2(X-x_1)(X-x_2)(X-x_3)^2+2(X-x_1)(X-x_2)^2(X-x_3).

En substituant $x_1$ il vient :

P'(x_1) = (x_1-x_2)^2(x_1-x_3)^2.

Comme $x_1-x_2\neq 0$ et $x_1-x_3\neq 0$ il suffit de poser :

H(X) = \frac{(X-x_1)(X-x_2)^2(X-x_3)^2}{(x_1-x_2)^2(x_1-x_3)^2}.

Le polynôme $H$ convient.

10/03/2022 - Polynome de hermite avec 3 racines dont 2 doubles
Représentation graphique du polynôme $H$ lorsque $x_1=1,$ $x_2=2$ et $x_3=3.$

Traitez un autre cas avec trois points

Vous cherchez un polynôme $K$ tel que :

\left\{\begin{align*}
K(x_1) &= 1\\
K(x_2) &= 0\\
K(x_3) &= 0\\
K'(x_1) &= 0\\
K'(x_2) &= 0\\
K'(x_3) &= 0.
\end{align*}\right.

Vous constatez que $x_2$ et $x_3$ sont des racines doubles du polynôme $K$.

Il semble judicieux de travailler avec le polynôme $Q$ défini par $Q(X)=(X-x_2)^2(X-x_3)^2.$

Ce polynôme satisfait quatre conditions :

\left\{\begin{align*}
Q(x_2) &= 0\\
Q(x_3) &= 0\\
Q'(x_2) &= 0\\
Q'(x_3) &= 0.
\end{align*}\right.

Vous allez maintenant poser $R(X) = Q(X)-Q'(x_1)H(X)$ où $H$ est le polynôme de Hermite défini dans le paragraphe précédent. En dérivant, $R'(X) = Q'(X)-Q'(x_1)H'(X)$ puis :

\left\{\begin{align*}
R(x_2) &= Q(x_2)-Q'(x_1)H(x_2) = 0-Q'(x_1)\times 0 = 0\\
R(x_3) &= Q(x_3)-Q'(x_1)H(x_3) = 0 - Q'(x_1)\times 0 = 0\\
R'(x_2) &= Q'(x_2)-Q'(x_1)H'(x_2) = 0 - Q'(x_1)\times 0 = 0\\
R'(x_3) &= Q'(x_3)-Q'(x_1)H'(x_3) = 0 - Q'(x_1)\times 0 = 0\\
R'(x_1) &= Q'(x_1)-Q'(x_1)H'(x_1) =  Q'(x_1)-Q'(x_1)\times 1 = 0.
\end{align*}\right.

D’autre part, $R(x_1) = Q(x_1)-Q'(x_1)H(x_1) = Q(x_1)-Q'(x_1)\times 0 = Q(x_1) = (x_1-x_2)^2(x_1-x_3)^2$ qui est non nul.

En posant $K(X) = \frac{Q(X)-Q'(x_1)H(X)}{(x_1-x_2)^2(x_1-x_3)^2}$ vous constatez que le polynôme $K$ convient.

10/03/2022 - Polynome de hermite avec 2 racines dont 2 doubles
Courbe du polynôme $K$ lorsque $x_1 = 1,$ $x_2=2$ et $x_3=3.$

Positivité et négativité

Les représentations graphiques précédentes peuvent amener à penser que, si $m$ est le plus petit des $x_i$ et si $M$ le plus grand des $x_i$, alors les polynômes de Hermite seraient positifs sur l’intervalle $[m,M].$

Cette conjecture s’avère erronée.

Par exemple, le polynôme $H$ défini par $H(X) = \frac{(X-9)^2 (X-2)^2 (X+2)^2 (17 X-74)}{364500}$ est un polynôme de Hermite associé aux quatre points d’interpolation suivants : $x_1 = -2, x_2=2, x_3=7$ et $x_4=9.$ La présence du polynôme du premier degré $17X-74$ dans la factorisation de $H$ montre un changement de signe à l’intérieur de l’intervalle $[-2,9].$

La représentation graphique du polynôme $H$ ci-dessous permet de faire état d’un tel revirement.

12/03/2022 - Changement de signe dun polynome de hermite d95b20aef4f3c49f5c2745410d16a3bc3d445412
La courbe du polynôme $H$ est tantôt au-dessous de l’axe des abscisses, tantôt au-dessus

Prolongement

Le cas général sera traité dans l'article 229. De tels polynômes sont utiles dans la démonstration de l’approximation uniforme d’une fonction continue par un polynôme sur un segment. Les polynômes de Lagrange ne suffisent pas à cause du phénomène explicité par Runge.

227. Polynômes de Tchebychev, relation de récurrence, degré et coefficient dominant

Dans cet article vous cherchez à exprimer, pour tout entier naturel $n$ et pour tout réel $x$, $\cos nx$ en fonction de $\cos x.$

Pour y parvenir vous utiliserez les formules d’addition et de soustraction de la fonction cosinus :

\begin{align*}
\forall (a,b)\in\R^2, \cos(a+b) &=\cos a\cos b-\sin a\sin b\\
\forall (a,b)\in\R^2, \cos(a-b) &=\cos a\cos b+\sin a\sin b.
\end{align*}

Exprimez $\cos 2x$ en fonction de $\cos x$

Soit $x$ un nombre réel.

D’après la formule d’addition :

\cos 2x=\cos(x+x) = \cos x \cos x - \sin x \sin x.

D’après la formule de soustraction :

1 = \cos 0=\cos(x-x) = \cos x \cos x + \sin x \sin x.

Par somme des deux égalités ci-dessus, il vient :

1+\cos 2x = 2\cos^2 x\\
\cos 2x = 2\cos^2 x-1.

Notez $\boxed{P_2(X) = 2X^2-1.}$ Alors $\forall x\in\R, \cos(2x) = P_2(\cos x).$

En notant $\boxed{P_1(X) =X}$ vous avez aussi $\forall x\in\R, \cos(1x) = P_1(\cos x).$

De même, en notant $\boxed{P_0(X)=1}$ vous avez encore $\forall x\in\R, 1 = \cos(0x) = P_0(\cos x).$

Construisez une preuve par récurrence

Vous allez utiliser une preuve à deux crans.

Pour tout entier naturel $n$, notez $\mathscr{P}(n)$ la propriété : « $\exists P_n\in\R[X], \forall x\in\R, \cos nx = P_n(\cos x)$ ».

Pour $n=0$, vous posez $P_0(X)=1$ et vous constatez que $\forall x\in\R, 1 = \cos(0x) = P_0(\cos x).$ Donc $\mathscr{P}(0)$ est vérifiée.

Pour $n=1$, vous posez $P_1(X) =X$ vous avez aussi $\forall x\in\R, \cos(1x) = P_1(\cos x).$ Donc $\mathscr{P}(1)$ est vérifiée.

Soit $n$ un entier naturel. Supposez $\mathscr{P}(n)$ et $\mathscr{P}(n+1).$

Il existe un polynôme $P_n\in\R[X]$ tel que $\forall x\in\R, \cos nx = P_n(\cos x).$

Il existe un polynôme $P_{n+1}\in\R[X]$ tel que $\forall x\in\R, \cos ((n+1)x) = P_{n+1}(\cos x).$

Soit $x\in\R.$ Alors :

\begin{align*}
\cos((n+2)x) &= \cos((n+1)x+x) = \cos ((n+1)x) \cos x-\sin ((n+1)x) \sin x\\
\cos nx &= \cos((n+1)x-x) =\cos ((n+1)x) \cos x + \sin ((n+1)x) \sin x.
\end{align*}

Par somme, il vient :

\begin{align*}
\cos((n+2)x)+\cos nx&=2\cos ((n+1)x) \cos x \\
\cos ((n+2)x) &= 2(\cos x) P_{n+1}(\cos x)-P_n(\cos x)
\end{align*}

Posez $P_{n+2}(X) = 2X P_{n+1}(X)-P_n(X).$

Vous déduisez $\cos ((n+2)x) = P_{n+2}(\cos x).$

Ainsi, pour tout entier naturel $n$, la propriété $\mathscr{P}(n)$ est vraie.

Définissez la suite des polynômes de Tchebychev

Vous définissez la suite de polynômes de $\Z[X]$ dite de Tchebychev, en posant $P_0 = 1,$ $P_1 = X$ puis en utilisant la relation de récurrence ci-dessous :

\boxed{\forall n\in\N, P_{n+2} = 2XP_{n+1}-P_n.}

D’après la récurrence qui a été effectuée ci-dessus, cette suite de polynômes vérifie l’égalité suivante :

\boxed{\forall n\in\N, \forall x\in\R, \cos(nx) = P_n(\cos x).}

Déterminez le coefficient dominant de $P_n$

Les premiers calculs montrent que :

\begin{align*}
P_0(X) &= 1\\
P_1(X) &= X\\
P_2(X) &= 2X^2-1\\
P_3(X) &= 2X(2X^2-1)-X = 4X^3-3X\\
P_4(X) &=2X(4X^3-3X)-(2X^2-1) = 8X^4-8X^2+1.
\end{align*}

Il semble que pour tout entier naturel $n$ supérieur ou égal à $1,$ le degré du polynôme $P_n$ soit égal à $n$ et que son coefficient dominant soit égal à $2^{n-1}.$

Vous allez de nouveau démontrer cette propriété par récurrence à deux crans.

Pour tout entier naturel $n$ non nul, notez $\mathscr{P}(n)$ la propriété : « le polynôme $P_n$ est de degré $n$ et son coefficient dominant est égal à $2^{n-1}$ ».

Pour $n=1$, $P_1 = X = 2^0 X$ donc la propriété $\mathscr{P}(1)$ est vérifiée.

Pour $n=2$, $P_2 = 2X^2-1 = 2^1 X^2-1$ donc la propriété $\mathscr{P}(2)$ est vérifiée.

Soit maintenant $n$ un entier supérieur ou égal à $1.$ Supposez $\mathscr{P}(n)$ et $\mathscr{P}(n+1).$

Par hypothèse de récurrence, il existe $(a_0,\dots,a_{n})\in\R^{n+1}, P_n(X) = \sum_{k=0}^{n} a_k X^k.$

De même, il existe $(b_0,\dots,b_{n})\in\R^{n+1}, P_{n+1}(X) = 2^n X^{n+1} + \sum_{k=0}^n b_kX^k.$

Alors :

\begin{align*}
P_{n+2} &= 2XP_{n+1}-P_n\\
&=2X\left(2^n X^{n+1} + \sum_{k=0}^n b_kX^k\right)  - \sum_{k=0}^{n} a_k X^k\\
&= 2^{n+1}X^{n+2}+ \sum_{k=0}^n 2b_kX^{k+1}- \sum_{k=0}^{n} a_k X^k\\
&= 2^{n+1}X^{n+2}+ \sum_{k=1}^{n+1} 2b_{k-1}X^{k}- \sum_{k=0}^{n} a_k X^k\\
&= 2^{n+1}X^{n+2}+ 2b_nX^{n+1}+\sum_{k=1}^{n} 2b_{k-1}X^{k}- \sum_{k=1}^{n} a_k X^k - a_0\\
&= 2^{n+1}X^{n+2}+ 2b_nX^{n+1}+\sum_{k=1}^{n} (2b_{k-1}-a_k)X^{k} - a_0.
\end{align*}

Le polynôme $P_{n+2}$ est de degré $n+2$ et admet $2^{n+1}$ pour coefficient dominant.

Par récurrence, vous avez démontré que, pour tout entier naturel $n$ supérieur ou égal à $1,$ le degré du polynôme $P_n$ est égal à $n$ et son coefficient dominant est égal à $2^{n-1}.$

226. Calculez une intégrale avec la méthode de Simpson (3/3)

Considérez une fonction $f$ définie et quatre fois dérivable sur un intervalle $[a,b]$ où $a$ et $b$ sont des réels tels que $a<b.$

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

Dans cet article vous allez déterminer une majoration de l’erreur commise par la méthode de Simpson, d’abord avec l’intervalle $[a,b]$ en entier, puis lorsque celui-ci est subdivisé en $n$ parties de même amplitude.

Pour plus de commodité, notez $m = \frac{a+b}{2}.$

Il a vu, au sein du contenu écrit dans l'article 224 que l’on peut considérer que :

\begin{align*}
\int_a^bf(t)\dt\approx(b-a)\left(\frac{1}{6}f(a)+\frac{2}{3}f\left(m\right)+\frac{1}{6}f(b)\right).
\end{align*}

Les calculs faits dans l’article précité ont montré qu’il existe un polynôme $P\in\R_2[X]$ tel que :

\left\{\begin{align*}
\int_a^b P(t)\dt&=(b-a)\left(\frac{1}{6}f(a)+\frac{2}{3}f\left(m\right)+\frac{1}{6}f(b)\right)\\
P(a)&=f(a)\\
P(b)&=f(b)\\
P\left(m\right)&=f\left(m\right).
\end{align*}\right.

Dans cet article, vous allez estimer l’erreur en utilisant le théorème de Rolle. Mais avant, pour faire apparaître la dérivée quatrième de $f$, il convient d’utiliser un polynôme du quatrième degré, et non pas du troisième, comme cela avait été le cas dans le contenu écrit dans l'article 225.

Une intégrale nulle

Vous allez établir que le polynôme $(X-a)(X-m)(X-b)$ a une intégrale nulle sur l’intervalle $[a,b].$

Vous effectuez le changement de variable $u=t-m$ ce qui conduit à :

t-a = u+m-a=u+\frac{b-a}{2}\\
t-b = u+m-b = u-\frac{b-a}{2}.
\begin{align*}
\int_a^b (t-a)(t-m)(t-b)\dt &= \int_{-(b-a)/2}^{(b-a)/2} \left(u^2-\frac{(b-a)^2}{4}\right)u\du\\
&= \int_{-(b-a)/2}^{(b-a)/2}u^3\du-\frac{(b-a)^2}{4} \int_{-(b-a)/2}^{(b-a)/2}u\du\\
&=0-\frac{(b-a)^2}{4}\times 0\\
&=0.
\end{align*}

Construisez un réel $\lambda$ bien choisi

Vous allez construire une fonction dont la dérivée s’annule en $m.$

Pour comprendre ceci, notez $R(X)$ le polynôme défini par :

R(X)=(X-a)(X-m)(X-b).

Le polynôme dérivé de $R$ est le polynôme suivant :

R'(X)=(X-m)(X-b)+(X-a)(X-b)+(X-a)(X-m).

De cette égalité vous déduisez $R'(m) = (m-a)(m-b)$ qui est non nul.

Soit alors $\lambda$ le réel défini par l’égalité suivante :

\lambda=\frac{f'(m)-P'(m)}{R'(m)}.

Utilisez le théorème de Rolle

Considérez la fonction $g$, définie sur l’intervalle $[a,b]$, par $g(x) = f(x)-P(x)-\lambda R(x).$

Comme $\forall x\in[a,b], g'(x) = f'(x)-P'(x)-\lambda R'(x)$ il vient $g'(m)=0.$

D’autre part :

\left\{\begin{align*}
g(a) &= f(a)-P(a) = 0\\
g(b) &= f(b)-P(b) = 0\\
g(m) &= f(m)-P(m) = 0.
\end{align*}\right.

Soit maintenant un réel $t\in[a,b]$ fixé tel que $t\notin\{a,b,m\}.$

Vous considérez le réel $A$ défini par :

A = \frac{f(t)-P(t)-\lambda R(t)}{(t-a)(t-m)^2(t-b)}.

Notez alors $h$ la fonction définie sur $[a,b]$ par :

\begin{align*}
\forall x\in[a,b], h(x)&=f(x)-P(x)-\lambda R(x) - A(x-a)(x-m)^2(x-b)\\
&=g(x)-A(x-a)(x-m)^2(x-b).
\end{align*}

Notez que $h(a) = h(b) = h(m) = h(t) = 0.$

La fonction $h$ admet au moins 4 zéros dans l’intervalle $[a,b].$

Par conséquent, l’application du théorème de Rolle à la fonction $h$ permet d’affirmer que la fonction $h’$ possède au moins 3 zéros sur l’intervalle ouvert $]a,b[$, aucun d’entre eux n’appartenant à l’ensemble $\{a,m,b,t\}.$

Or, $m$ est racine double du polynôme $(X-a)(X-m)^2(X-b)$ donc $m$ est racine de son polynôme dérivé.

Notez $S(X) = (X-a)(X-m)^2(X-b)$ vous avez $S'(m) = 0.$

Comme $\forall x\in [a,b], h(x) = g(x)-AS(x)$ vous déduisez $\forall x\in [a,b], h'(x) = g'(x)-AS'(x)$ et donc $h'(m)=g'(m)-AS'(m).$ Comme $g'(m)=S'(m) = 0$ vous déduisez $h'(m)=0.$

Vous déduisez donc que la fonction $h’$ admet au moins 4 zéros dans l’intervalle $]a,b[.$

En appliquant le théorème de Rolle à $h’$, vous déduisez que la fonction $h »$ admet au moins 3 zéros dans l’intervalle $]a,b[.$

En appliquant le théorème de Rolle à $h »$, vous déduisez que la fonction $h »’$ admet au moins 2 zéros dans l’intervalle $]a,b[.$

En appliquant le théorème de Rolle à $h »’$, vous déduisez que la fonction $h » »$ admet au moins un zéro dans l’intervalle $]a,b[.$ Il existe un réel $c\in]a,b[$ tel que $h » »(c)=0.$

Or, $\forall x\in[a,b], h » »(x) = f » »(x)-P » »(x)-\lambda R » »(x)-24A.$

Comme $P\in\R_2[X]$ il vient $P » »(X) = 0.$ Comme $R\in\R_3[X]$ il vient $R » »(X)=0.$

Par suite il vient $h » »(c)=f » »(c)-24A$ d’où $A = \frac{f » »(c)}{24}.$

Ainsi $h(t)=0$ fournit $f(t) = P(t)+\lambda R(t) + \frac{f » »(c)}{24}(t-a)(t-m)^2(t-b).$

Maintenant soit $t\in\{a,b,m\}.$ Comme $f(t)-P(t)-\lambda R(t) = 0$ et comme $(t-a)(t-m)^2(t-b)=0$ vous posez $c=m$ et vous avez encore $f(t) = P(t)+\lambda R(t) + \frac{f » »(c)}{24}(t-a)(t-m)^2(t-b).$

En définitive :

\boxed{\begin{align*}
\forall t\in[a,b], \exists c\in]a,b[, f(t)&=P(t)+\lambda R(t)+\frac{f''''(c)}{24}S(t)\\
&=P(t)+\lambda (t-a)(t-m)(t-b)+\frac{f''''(c)}{24}(t-a)(t-m)^2(t-b).
\end{align*}}

Vous déduisez la majoration suivante:

\forall t\in[a,b], \left\vert f(t)-P(t)-\lambda(t-a)(t-m)(t-b)\right\vert \leq \frac{M_4}{24}\vert t-a\vert (t-m)^2\vert t-b\vert.

Procédez à la majoration effective de l’erreur

Compte tenu du fait que

\int_a^b (t-a)(t-m)(t-b)\dt = 0

vous déduisez successivement que :

\begin{align*}
\left\vert \int_a^b f(t)\dt - (b-a)\left(\frac{1}{6}f(a)+\frac{2}{3}f\left(m\right)+\frac{1}{6}f(b)\right) \right\vert &\leq \left\vert \int_a^b f(t)\dt - \int_a^b P(t)\dt \right\vert\\
&\leq \left\vert \int_a^b f(t)\dt - \int_a^b P(t)\dt - \lambda \int_a^b (t-a)(t-m)(t-b)\dt \right\vert\\
&\leq \left\vert  \int_a^b \left[f(t)- P(t) - \lambda  (t-a)(t-m)(t-b)\right]\dt \right\vert\\
&\leq  \int_a^b \left\vert f(t)- P(t) - \lambda  (t-a)(t-m)(t-b)\right\vert\dt \\
&\leq \frac{M_4}{24} \int_a^b \vert t-a\vert (t-m)^2\vert t-b\vert\dt.
\end{align*} 

Vous effectuez le changement de variable $u=t-m$ pour avoir à nouveau $t-a = u+\frac{b-a}{2}$ et $t-b = u-\frac{b-a}{2}.$

\begin{align*}
\int_a^b \vert t-a\vert (t-m)^2\vert t-b\vert\dt &= \int_{-(b-a)/2}^{(b-a)/2}   \left\vert u+\frac{b-a}{2}\right\vert u^2\left\vert u-\frac{b-a}{2}\right\vert \du\\
&= \int_{-(b-a)/2}^{(b-a)/2}   \left\vert u^2-\frac{(b-a)^2}{4}\right \vert u^2 \du\\
&=2 \int_{0}^{(b-a)/2}   \left\vert u^2-\frac{(b-a)^2}{4}\right \vert u^2 \du\\
&=2 \int_{0}^{(b-a)/2}   \left( \frac{(b-a)^2}{4}-u^2\right) u^2 \du\\
&= \frac{(b-a)^2}{4}  \int_{-(b-a)/2}^{(b-a)/2}   u^2\du - \int_{-(b-a)/2}^{(b-a)/2}   u^4  \du \\
&=  \frac{(b-a)^2}{2}  \int_{0}^{(b-a)/2}   u^2\du - 2\int_{0}^{(b-a)/2}   u^4  \du \\
&= \frac{(b-a)^2}{2} \frac{((b-a)/2)^3}{3} - \frac{2((b-a)/2)^5}{5} \\
&= \frac{(b-a)^5}{2^4\times 3} - \frac{(b-a)^5}{5\times 2^4}  \\
&= \frac{5(b-a)^5}{2^4\times 3\times 5} - \frac{3(b-a)^5}{3\times5\times 2^4}  \\
&= \frac{2(b-a)^5}{2^4\times 3\times 5} \\
&= \frac{(b-a)^5}{2^3\times 3\times 5} \\
&= \frac{(b-a)^5}{2^2\times 3\times 10} \\
&= \frac{(b-a)^5}{120}.
\end{align*}

Mis bout à bout vous obtenez :

\begin{align*}
\left\vert \int_a^b f(t)\dt - (b-a)\left(\frac{1}{6}f(a)+\frac{2}{3}f\left(m\right)+\frac{1}{6}f(b)\right) \right\vert 
&\leq \frac{M_4}{24} \int_a^b \vert t-a\vert (t-m)^2\vert t-b\vert\dt \\
&\leq \frac{M_4}{24} \frac{(b-a)^5}{120}\\
&\leq  \frac{M_4(b-a)^5}{2880}.
\end{align*} 

Ce résultat constitue la majoration recherchée autour de cet article.

Pour toute fonction $f$ quatre fois dérivable sur l’intervalle $[a,b]:$

\boxed{\left\vert \int_a^b f(t)\dt - (b-a)\left(\frac{1}{6}f(a)+\frac{2}{3}f\left(m\right)+\frac{1}{6}f(b)\right) \right\vert \leq  \frac{M_4(b-a)^5}{2880}.} 

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}.$

Autrement dit, $\forall k\in\llbracket 0, n-1\rrbracket, t_k = \frac{a_k+a_{k+1}}{2}.$

\begin{align*}
\left\vert\int_a^b f(t)\dt-\frac{b-a}{n}\sum_{k=0}^{n-1} \left(\frac{1}{6}f(a_k)+\frac{2}{3}f(t_k)+\frac{1}{6}f(a_{k+1})\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} \left(\frac{1}{6}f(a_k)+\frac{2}{3}f(t_k)+\frac{1}{6}f(a_{k+1})\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} \left(\frac{1}{6}f(a_k)+\frac{2}{3}f(t_k)+\frac{1}{6}f(a_{k+1}) \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} \left(\frac{1}{6}f(a_k)+\frac{2}{3}f(t_k)+\frac{1}{6}f(a_{k+1}) \right) \right\vert \\
&\leq  \sum_{k=0}^{n-1}\frac{M_4((b-a)/n)^5}{2880}\\
&\leq  \sum_{k=0}^{n-1}\frac{M_4(b-a)^5}{2880n^5}\\
&\leq  \frac{M_4(b-a)^5}{2880n^4}.
\end{align*}

Dans le cas d’un intervalle $[a,b]$ partagé en $n$ parties égales, la méthode de Simpson amène à une estimation de l’intégrale avec une majoration qui est en $O\left(\frac{1}{n^4}\right).$

Vous en déduisez la formule composite de Simpson :

\boxed{
\left\vert\int_a^b f(t)\dt-\frac{b-a}{n}\sum_{k=0}^{n-1} \left(\frac{1}{6}f(a_k)+\frac{2}{3}f(t_k)+\frac{1}{6}f(a_{k+1})\right)\right\vert \leq  \frac{M_4(b-a)^5}{2880n^4}.}

Prolongement

Pourriez-vous utiliser ce résultat pour déterminer les cinq premières décimales exactes du nombre $\pi$ ? Quelle fonction à intégrer utiliseriez-vous ? Sur quel intervalle ?