Analyse numérique et calcul scientifique/Interpolation polynomiale

Une page de Wikiversité.


Interpolation polynomiale
Nuvola apps edu mathematics-p.svg
Chapitre 1
Leçon : Analyse numérique et calcul scientifique
Retour au Sommaire
Chap. suiv. : Splines
Icon falscher Titel.svg

En raison de limitations techniques, la typographie souhaitable du titre, « Analyse numérique et calcul scientifique : Interpolation polynomiale
Analyse numérique et calcul scientifique/Interpolation polynomiale
 », n'a pu être restituée correctement ci-dessus.

Sommaire

[modifier] Interpolation polynomiale

[modifier] Introduction

Dans ce chapitre, le but est d'interpoler un ensemble de points (x_i,y_i)_{i\in[0,n]} par une fonction polynomiale P. C'est-à-dire trouver les coefficients \{a_i\}_{i\in[0,n+1]} définissant P telle que P(x)=\sum_{i=1}^{n+1} a_{i}x^i et \forall i\in[0,n]\; P(x_i)=y_i.

On pourra aussi interpoler une fonction f en un ensemble de points \{x_i\}_{i\in[0,n]}, c'est-à-dire trouver P telle que \forall i\in[0,n]\; P(x_i)=f(x_i)


[modifier] Matrice de Vandermonde

On peut exprimer sous la forme d'une matrice :

\begin{bmatrix}
x_0^{n+1} & x_0^{n} & x_0^{n-1} & \ldots & x_0 & 1 \\
x_1^{n+1} & x_1^{n} & x_1^{n-1} & \ldots & x_1 & 1 \\
\vdots & \vdots & \vdots & & \vdots & \vdots \\
x_n^{n+1} & x_n^{n} & x_n^{n-1} & \ldots & x_n & 1 
\end{bmatrix}
\begin{bmatrix}
a_{n+1} \\
a_{n} \\
\vdots \\
a_0
\end{bmatrix}
=
\begin{bmatrix}
y_0 \\
y_1 \\
\vdots \\
y_n
\end{bmatrix}

appelée matrice de Vandermonde.
Son déterminant vaut \prod_{1\le i<j\le n} (x_j-x_i).

Le système admet une solution unique si le déterminant de Vandermonde est non nul. Ce qui prouve que pour faire passer un polynôme unique par n+1 points distincts celui-ci doit être au plus de degré n.

[modifier] Interpolation Lagrangienne

Soient les n+1 points (x_i,y_i)_{i\in[0,n]} à interpoler par un polynôme P de degré n.

Soient les n+1 polynômes l_j(x)\;\; j\in [0,n]

l_{j}(x):=\prod_{i=0,i\neq j}^n\frac{x-x_i}{x_j-x_i}

Les principales propriétés de ces polynômes sont:

  • l_{j}(x_{i})=\delta_{i,j}=\begin{cases} 
1 & \mbox{si } i=j  \\ 
0 & \mbox{si } i \neq j \end{cases}
  • lj est de degré n pour tout j

On définit le polynôme d'interpolation de Lagrange :

P(x) = \sum_{j=0}^{n} y_j l_j(x)

tel que : P(x_j) = \sum_{j=0}^{n} y_{j} l_{j}(x_j)=y_j