Interpolation Lagrange
Interpolation Lagrange
Interpolation Lagrange
1) Isomorphisme fondamental
Alors lapplication linaire u : Kn 1 [X] ! K n P 7 ! (P (a1 ); :::; P (an )) est bijective (isomorphisme linaire).
Autrement dit, un polynme de degr n 1 est entirement dni par (P (a1 ); :::; P (an )) 2 K n :
dem : u est injective : Si P 2 Ker u, alors les ai sont racines de P , et comme deg P n 1, alors P = 0:
Remarque : Si on prend u : K[X] ! K n P 7 ! (P (a1 ); :::; P (an )), le noyau Ker u est le sev des multiples de
Lorsque deux aj sont gaux, alors V (a1 ; a2 ; :::; an ) contient deux lignes gales donc nest pas inversible.
On en dduit que V (a1 ; a2 ; :::; an ) est inversible ssi les aj sont deux deux distincts .
Remarque : La matrice de Van der Monde est donc la matrice du systme linaire dont les inconnues sont les
coe cients de P = + n 1,
j 0 1X + ::: + nX et dont les quations sont les P (ai ) = yi pour tout i 2 [[1; n]].
Q
Remarque : On peut dmontrer par ailleurs que det V (a1 ; a2 ; :::; an ) = i<j (aj ai ):
On retrouve ainsi que la matrice V (a1 ; a2 ; :::; an ) est inversible ssi les aj sont deux deux distincts.
a) Polynmes de Lagrange
Q X aj
On considre les polynmes de Lagrange 8i 2 [[1; n]]; Li (X) = j6=i 2 Kn 1 [X].
ai aj
On a 8(i; j) 2 [[1; n]]2 , Li (aj ) = ij (cest--dire 1 si i = j et 0 sinon).
Ainsi, u(Li ) = Ei = (0; 0; :::; 0; 1; 0:; :::; 0) le i-ime vecteur de la base canonique de K n :
b) Isomorphisme rciproque
On cherche expliciter lunique polynme P 2 Kn 1 [X] tel que u(P ) = Y , cest--dire 8j, P (ai ) = yi :
On applique le principe de superposition :
P
En crivant Y = ni=1 yi Ei , on se ramne dterminer les solutions des quations u(Li ) = Ei :
est
(X b)(X c) (X a)(X c) (X a)(X b)
P = + +
(a b)(a c) (b a)(b c) (c a)(c b)
B, cest--dire les coe cients du polynme Lj . Par exemple, les coe cients situs sur la dernire ligne sont les
Q
coe cients dominants des Lj , cest--dire les 1 i n et i6=j aj 1 ai :
Tout polynme P 2 Kn 1 [X] est lunique polynme valant P (ai ) en ai , pour tout i 2 [[1; n]].
Pn
On a 8P 2 Kn 1 [X], P (X) = i=1 P (ai )Li (X)
1) Interpolation de Lagrange
a) Proprit fondamentale
Pour tout (y1 ; :::; yn ) 2 K n , il existe un unique P 2 Kn 1 [X] tel que 8i 2 [[1; n]], P (ai ) = yi :
Preuve : On montre que u : Kn 1 [X] ! K n P 7 ! (P (a1 ); :::; P (an )) est linaire bijective.
Le polynme dinterpolation de f en les aj est lunique polynme P 2 Kn 1 [X] tel que 8i 2 [[1; n]], P (ai ) = f (ai ):
n
X Y X aj
P (X) = f (ai )Li (X), avec Li (X) =
ai aj
i=1 j6=i
(x b) (x a)
Exemple : La fonction a ne interpolant f en a et b est L(x) = f (a) + f (b) :
(a b) (b a)
Preuve : On a M = BQ + R, avec deg R n 1 et R(ai ) = M (ai ) pour tout i 2 f1; 2; :::; ng.
f (b) f (a)
L(x) = + (x a), avec = f (a) et =
b a
Exemple : Pour dterminer le polynme P de degr 2 vriant P (a) = f (a), P (b) = f (b) et P (c) = f (c), il est
2) Interpolation de Taylor
a) Proprit fondamentale
Prop : Soient a 2 K et n 2 N.
Pour tout (y0 ; y1 ; :::; yn 1) 2 K n , il existe un unique P 2 Kn 1 [X] tel que 8k 2 [[0; n 1]], P (k) (a) = yk :
Preuve : On montre que u : Kn 1 [X] ! K n P 7 ! (P (a); P 0 (a); :::; P (n 1) (a)) est linaire bijective.
(X a)k (j)
Le polynme Rk (X) = vrie Rk (a) = jk :
k!
Pn 1
On en dduit que lunique P 2 Kn 1 [X] tel que 8k 2 [[0; n 1]], P (k) (a) = yk est P (X) = k=0 yk Rk (X):
Soit f : K ! R.
Le polynme de Taylor de f en a lordre n est lunique P 2 Kn 1 [X] tel que 8k 2 [[0; n 1]], P (k) (a) = f (k) (a):
Pn (k) (a) (X a)k
On a alors P (X) = k=0 f .
k!
Remarque : En particulier, un polynme de degr n 1 est son propre polynme de Taylor.
Pn 1 (k) (X a)k
On obtient ainsi la formule de Taylor algbrique 8P 2 Kn 1 [X], P (X) = k=0 P (a) :
k!
3) Interpolation de Hermite
Lapplication u est injective : Si u(P ) = 0, alors les aj sont des racines au moins doubles de P .
Donc P admet au moins 2n racines comptes avec multiplicit, donc par degr, P = 0:
Comme dim K2n 1 [X] = 2n = dim K 2n , alors u est bijective (car injective).
Prenons le cas n = 2, avec P 2 K3 [X] dni par (P (a); P (b); P 0 (a); P 0 (b)) = (y1 ; y2 ; y3 ; y4 ):
Par le principe de superposition, on se ramne traiter les cas o (y1 ; y2 ; y3 ; y4 ) est lun des quatre vecteurs de la
base canonique. Par exemple, en prenant P = (X a)(X b)2 , on a P (a) = P (b) = P 0 (b) = 0. On explicite de
1
sorte que P 0 (a) = 0, do = (a b)2
:
4) Lien entre interpolation dun polynme et reste de ce polynme par une division euclidienne
bn an ban abn
Les valeurs obtenues pour a = b sont dailleurs les limites de et lorsque b tend vers a.
b a b a