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.
Étant donnés 
points
d'abscisses distinctes, les différences divisées sont définies de la manière suivante :
![{\displaystyle [y_{\nu }]=y_{\nu }\qquad (\nu =0,\ldots ,n)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ebf0e53e70626ac3be6b262d58968a6f434514ad)
![{\displaystyle [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).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d4e2f628167b5bd3e106aa3e616d47e91fd05744)
Pour toute fonction 
telle que 
, on note parfois ![{\displaystyle f[x_{0},\dots ,x_{n}]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/683c749df05d7cf790d0a90a9b90a8e4006d8b5d)
la différence divisée ![{\displaystyle [y_{0},\dots ,y_{n}]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e2bbd5be83762398c1f21265e7d625288cd3a555)
.
D'après le théorème d'interpolation de Newton, la différence divisée associée à points est égale au coefficient de degré 
du polynôme d'interpolation de Lagrange de ces points. Autrement dit :
![{\displaystyle [y_{0},\dots ,y_{n}]=\sum _{j=0}^{n}{\frac {y_{j}}{\prod _{0\leq i\leq n,\,i\neq j}(x_{j}-x_{i})}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/52c3554da2961eaebd2eedf1e0579a0b2c4d4896)
.
Cette égalité a des conséquences remarquables :
- invariance par permutation des indices :
![{\displaystyle f[x_{\sigma (0)},\dots ,x_{\sigma (n)}]=f[x_{0},\dots ,x_{n}]\quad (\sigma \in S_{\{0,\dots ,n\}})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e3bd200e79024c5df90f51d5175499b65c41c37d)
;
- linéarité :
![{\displaystyle (af+bg)[x_{0},\dots ,x_{n}]=a\,f[x_{0},\dots ,x_{n}]+b\,g[x_{0},\dots ,x_{n}]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/046c53fa8acdc64644111522565975ecacd7e658)
;
- règle de Leibniz :
![{\displaystyle (fg)[x_{0},\dots ,x_{n}]=\sum _{j=0}^{n}f[x_{0},\dots ,x_{j}]g[x_{j},\dots ,x_{n}]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/70510e11397c7183ca704543841400bf22366e3a)
;
- théorème de la moyenne : pour n ≥ 1 et
, si est de classe Cn–1 sur ![{\displaystyle [x_{0},x_{n}]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9edac6e40dff00287ca6ec6d62184b9217f610b5)
et possède une dérivée n-ième sur ![{\displaystyle ]x_{0},x_{n}[}](https://wikimedia.org/api/rest_v1/media/math/render/svg/56bab931a3856fbc023d483c1c6d0ca86c8533ea)
, il existe ![{\displaystyle c\in ]x_{0},x_{n}[}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2682ddf6b6cd8d4e6de4fa978d97655a851c7563)
tel que![{\displaystyle f[x_{0},\dots ,x_{n}]={\frac {f^{(n)}(c)}{n!}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/103e0df6254ad08c55e660905661dd4e0e9b6f48)
.
Les premières itérations donnent :
- Ordre 0 :
![{\displaystyle [y_{0}]=y_{0}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/be6b71d95d874915cb8659f9bd97988482733218)
- Ordre 1 :
![{\displaystyle [y_{0},y_{1}]={\frac {y_{1}-y_{0}}{x_{1}-x_{0}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b5cebb2b013c27e5c8c3943356fe2faee854b01e)
- Ordre 2 :
![{\displaystyle [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}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/17976ad41f73c89d4a4bd40611dba48569a10baa)
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 :
![{\displaystyle {\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}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d60ef92b00038923edcedaecccbadd25373b07c1)
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
![{\displaystyle \Delta _{h}^{0}[f]=f\quad {\text{et}}\quad \triangle _{h}^{n+1}[f](x)=\triangle _{h}^{n}[f](x+h)-\triangle _{h}^{n}[f](x)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9f6dc9595a3ea50aa00c9079d016bbfae5b4fd05)
,
par la relation (immédiate par récurrence) :
![{\displaystyle f[x,x+h,\ldots ,x+nh]={\frac {1}{n!h^{n}}}\triangle _{h}^{n}[f](x)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2f9af1fefb8f5a5dabe45cebad952ed656e705bc)
.
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.
Interpolation polynomiale (sic) de type Newton et différences divisées
