Différences divisées

{{#invoke:Bandeau|ébauche}} En mathématiques, les différences divisées correspondent à une discrétisation des dérivées successives d'une fonction. Ce sont des quantités définies et calculées de manière récursive en généralisant la formule du taux d'accroissement. Elles sont utilisées en particulier en interpolation newtonienne.

Définition

Étant donnés <math>n+1</math> points

<math>(x_0, y_0),\ldots,(x_n, y_n),</math>

d'abscisses distinctes, les différences divisées sont définies de la manière suivante :

<math>[y_{\nu}] = y_{\nu}\qquad(\nu=0,\ldots,n)</math>

<math>[y_{\nu},\ldots,y_{\nu+j}] = \frac{[y_{\nu+1},\ldots y_{\nu+j}] - [y_{\nu},\ldots y_{\nu+j-1}]}{x_{\nu+j}-x_{\nu}}\qquad(j=1,\ldots,n,\qquad \nu = 0,\ldots,n-j).</math>

Pour toute fonction <math>f</math> telle que <math>y_i=f(x_i)\qquad(i=0,\ldots,n)</math>, on note parfois <math>f[x_0,\dots,x_n]</math> la différence divisée <math>[y_0,\dots,y_n]</math>.

Propriétés

D'après le théorème d'interpolation de Newton, la différence divisée associée à <math>n+1</math> points est égale au coefficient de degré <math>n</math> du polynôme d'interpolation de Lagrange de ces points. Autrement dit :

<math>[y_0,\dots,y_n]=\sum_{j=0}^n\frac{y_j}{\prod_{0\le i\le n,\,i\neq j}(x_j-x_i)}</math>.

Cette égalité a des conséquences remarquables :

• invariance par permutation des indices : <math>f[x_{\sigma(0)},\dots,x_{\sigma(n)}]=f[x_0,\dots,x_n]\quad(\sigma\in S_{\{0,\dots,n\}})</math> ;
• linéarité : <math>(af+bg)[x_0,\dots,x_n]=a\,f[x_0,\dots,x_n]+b\,g[x_0,\dots,x_n]</math> ;
• règle de Leibniz : <math>(fg)[x_0,\dots,x_n]=\sum_{j=0}^nf[x_0,\dots,x_j]g[x_j,\dots,x_n]</math> ;
• théorème de la moyenne : pour n ≥ 1 et <math>x_0<\dots<x_n</math>, si <math>f</math> est de classe C^n–1 sur <math>[x_0,x_n]</math> et possède une dérivée n-ième sur <math>]x_0,x_n[</math>, il existe <math>c\in]x_0,x_n[</math> tel queModèle:Retrait

Exemples

Les premières itérations donnent :

Ordre 0 : <math>[y_0] = y_0</math>

Ordre 1 : <math>[y_0,y_1] = \frac{y_1-y_0}{x_1-x_0}</math>

Ordre 2 : <math>[y_0,y_1,y_2] = \frac{\frac{y_2-y_1}{x_2-x_1}-\frac{y_1-y_0}{x_1-x_0}}{x_2-x_0}</math>

Pour expliciter le processus récursif, les différences divisées peuvent être calculées en les disposant de la manière suivante dans un tableau :

<math>

\begin{matrix} x_0 & y_0 = [y_0] & & & \\

       &       & [y_0,y_1] &               & \\

x_1 & y_1 = [y_1] & & [y_0,y_1,y_2] & \\

       &       & [y_1,y_2] &               & [y_0,y_1,y_2,y_3]\\

x_2 & y_2 = [y_2] & & [y_1,y_2,y_3] & \\

       &       & [y_2,y_3] &               & \\

x_3 & y_3 = [y_3] & & & \\ \end{matrix} </math>

Dans le cas où les abscisses sont en progression arithmétique, les différences divisées sont reliées aux différences finies, définies par

<math>\Delta^0_h[f]=f\quad\text{et}\quad\triangle^{n+1}_h[f](x)=\triangle^n_h[f](x+h)-\triangle^n_h[f](x)</math>,

par la relation (immédiate par récurrence) :

<math>f[x,x+h,\ldots,x+nh]=\frac1{n!h^n}\triangle^n_h[f](x)</math>.

Application

Les différences divisées interviennent dans la formulation du théorème d'interpolation de Newton, qui donne une expression particulière du polynôme d'interpolation de Lagrange, permettant par exemple de démontrer que toute fonction polynomiale est égale à sa série de Newton.

Voir aussi

Lien externe

Interpolation polynomiale (sic) de type Newton et différences divisées