ATLAS Ouvrir dans l'atlas →

ExerciceMaths expertes, chapitre 6 — Applications des matrices · Pour aller plus loin

Exercice 25 Interpolation de Lagrange

DémonstrationCalcul

Soient y0y_{0}, y1y_{1} et y2y_{2} trois réels. On cherche une fonction polynôme PP de degré au plus 22, P(x)=ax2+bx+cP(x)=ax^{2}+bx+c, telle que P(−1)=y0P(-1)=y_{0}, P(0)=y1P(0)=y_{1} et P(1)=y2P(1)=y_{2}.

1. Montrer que ce problème s'écrit VX=YVX=Y, avec X=(abc)X=\begin{pmatrix} a \\ b \\ c \end{pmatrix} et Y=(y0y1y2)Y=\begin{pmatrix} y_{0} \\ y_{1} \\ y_{2} \end{pmatrix}, où VV est une matrice que l'on précisera.

2. Vérifier, sans calculatrice, que la matrice W=(12−112−12012010)W=\begin{pmatrix} \frac{1}{2} & -1 & \frac{1}{2} \\[1mm] -\frac{1}{2} & 0 & \frac{1}{2} \\[1mm] 0 & 1 & 0 \end{pmatrix} est l'inverse de VV. En déduire que le problème admet une unique solution, et exprimer aa, bb et cc en fonction de y0y_{0}, y1y_{1} et y2y_{2}.

3. On pose L0(x)=x(x−1)2L_{0}(x)=\dfrac{x(x-1)}{2}, L1(x)=1−x2L_{1}(x)=1-x^{2} et L2(x)=x(x+1)2L_{2}(x)=\dfrac{x(x+1)}{2}. Calculer Lj(−1)L_{j}(-1), Lj(0)L_{j}(0) et Lj(1)L_{j}(1) pour j∈{0;1;2}j\in\{0;1;2\}, puis montrer que P=y0L0+y1L1+y2L2P=y_{0}L_{0}+y_{1}L_{1}+y_{2}L_{2}. Quel lien y a-t-il entre les polynômes LjL_{j} et la matrice WW ?

4. Exprimer P(2)P(2) en fonction de y0y_{0}, y1y_{1} et y2y_{2}. Contrôler ce résultat avec y0=1y_{0}=1, y1=4y_{1}=4 et y2=9y_{2}=9.

Source : inspiré de poly « Mathématiques : du lycée aux CPGE scientifiques » (LLG / H4), n° 442 (transposition matricielle)

Corrigé

1. P(−1)=a−b+cP(-1)=a-b+c, P(0)=cP(0)=c et P(1)=a+b+cP(1)=a+b+c, d'où : V=(1−11001111)\boxed{V=\begin{pmatrix} 1 & -1 & 1 \\ 0 & 0 & 1 \\ 1 & 1 & 1 \end{pmatrix}}

2. On calcule VWVW ligne par ligne : VW=(12+12+0−1+0+112−12+00+0+00+0+10+0+012−12+0−1+0+112+12+0)=I3VW=\begin{pmatrix} \frac{1}{2}+\frac{1}{2}+0 & -1+0+1 & \frac{1}{2}-\frac{1}{2}+0 \\[1mm] 0+0+0 & 0+0+1 & 0+0+0 \\[1mm] \frac{1}{2}-\frac{1}{2}+0 & -1+0+1 & \frac{1}{2}+\frac{1}{2}+0 \end{pmatrix}=I_{3} Une matrice carrée WW telle que VW=I3VW=I_{3} est l'inverse de VV (chapitre 3), donc VV est inversible. D'après la Proposition 1, le problème admet une unique solution, X=WYX=WY : a=y0−2y1+y22,b=y2−y02,c=y1\boxed{a=\frac{y_{0}-2y_{1}+y_{2}}{2}, \qquad b=\frac{y_{2}-y_{0}}{2}, \qquad c=y_{1}}

3. On obtient le tableau suivant :

x=−1x=0x=1L0(x)100L1(x)010L2(x)001\begin{array}{|c||c|c|c|} \hline & x=-1 & x=0 & x=1 \\ \hline L_{0}(x) & 1 & 0 & 0 \\ \hline L_{1}(x) & 0 & 1 & 0 \\ \hline L_{2}(x) & 0 & 0 & 1 \\ \hline \end{array}

Le polynôme Q=y0L0+y1L1+y2L2Q=y_{0}L_{0}+y_{1}L_{1}+y_{2}L_{2} est de degré au plus 22 et vérifie Q(−1)=y0Q(-1)=y_{0}, Q(0)=y1Q(0)=y_{1}, Q(1)=y2Q(1)=y_{2}. Par unicité (question 2), Q=P\boxed{Q=P}. Chaque LjL_{j} vaut 11 en un point et 00 aux deux autres, ce qui rend la formule immédiate.

Développons : L0(x)=12x2−12xL_{0}(x)=\frac{1}{2}x^{2}-\frac{1}{2}x, L1(x)=−x2+1L_{1}(x)=-x^{2}+1, L2(x)=12x2+12xL_{2}(x)=\frac{1}{2}x^{2}+\frac{1}{2}x. Les coefficients (a ;b ;c)(a\,;b\,;c) de L0L_{0}, L1L_{1}, L2L_{2} sont exactement les colonnes de WW. C'est logique : L0L_{0} est le polynôme obtenu pour Y=(100) ⁣⊤Y=\begin{pmatrix} 1 & 0 & 0 \end{pmatrix}^{\!\top}, et WYWY est alors la première colonne de WW.

4. L0(2)=2×12=1L_{0}(2)=\frac{2\times 1}{2}=1, L1(2)=1−4=−3L_{1}(2)=1-4=-3 et L2(2)=2×32=3L_{2}(2)=\frac{2\times 3}{2}=3, donc : P(2)=y0−3y1+3y2\boxed{P(2)=y_{0}-3y_{1}+3y_{2}} Avec y0=1y_{0}=1, y1=4y_{1}=4, y2=9y_{2}=9, on trouve P(2)=1−12+27=16P(2)=1-12+27=16. C'est cohérent : ces valeurs sont celles de P(x)=(x+2)2P(x)=(x+2)^{2} en −1-1, 00, 11, et (2+2)2=16(2+2)^{2}=16. La formule permet de « prolonger » une suite de trois valeurs par un polynôme de degré 22 : on retrouve que la suite des carrés 11, 44, 99 se poursuit par 1616.

Dans l'édifice

Thèmes

matrice inverse, interpolation

Notions

inversible

Ouvrir dans l'atlas interactif →