- Différence finie
-
En mathématique, une différence finie est une expression de la forme f(x + b) − f(x + a). (où f est une fonction numérique) ; la même expression divisée par b − a s'appelle un taux d'accroissement (ou taux de variation), et il est possible, plus généralement, de définir de même des différences divisées. L'approximation des dérivées par des différences finies joue un rôle central dans les méthodes des différences finies utilisées pour la résolution numérique des équations différentielles, tout particulièrement pour les problèmes de conditions aux limites.
Les opérateurs mettant en jeu des différences finies, appelés opérateurs de différence, font l'objet d'une branche de l'analyse mathématique possédant de nombreux points communs avec la théorie des opérateurs différentiels ; certaines de ces analogies ont été systématisées et exploitées par le calcul ombral.
Sommaire
Opérateurs de différence finie
On considère en général uniquement les différences « en avant », « en arrière » et « centrées ».
L'opérateur de différence avant, Δh, est défini par
 = f(x + h) - f(x), \](7/c57d07fb3b28bc4d17061cb79540f558.png)
où h est le plus souvent une constante, qu'on appelle le pas.
De même, l'opérateur de différence arrière,
, est défini parEnfin, l'opérateur de différence centrée, δh, est défini par
Relation avec les dérivées
Le nombre dérivé d'une fonction f en un point x est défini par la limite
Si h a une valeur fixée (non nulle), le terme de droite est
Ainsi, la différence avant divisée par h approxime le nombre dérivé quand h est petit. L'erreur de cette approximation peut être déterminée à l'aide de la formule de Taylor. Supposant que f est continument dérivable,
La même formule est valable pour la différence arrière :
Cependant, la différence centrale donne une approximation plus précise ; son erreur est proportionnelle au carré du pas (si f est de classe
) :Différences d'ordre supérieur
De manière analogue, on peut obtenir des approximations de dérivées d'ordre supérieur. Ainsi, appliquant la formule de différence centrale précédente à la dérivée de f' en x, et utilisant les différences centrales pour f'(x + h / 2) et f'(x − h / 2), on obtient une approximation (par différences centrales) de la dérivée seconde de f:
Plus généralement, les différences avant, arrière et centrale d'ordre n sont respectivement données par
On verra (dans le paragraphe consacré au calcul symbolique des différences) que la notation
correspond au n-ème itéré de Δh, c'est-à-dire que
(n fois).Il faut remarquer que dans le cas des différences centrales, pour tout n impair, h est multiplié par des non-entiers; ce qui peut être un problème dans des applications où cela change l'intervalle de discrétisation ; dans ce cas, une solution peut être de prendre la moyenne de δn[f](x − h / 2) et de δn[f](x + h / 2).
La relation entre ces différences itérées et les dérivées d'ordre correspondant est très simple :
Les différences d'ordre supérieur peuvent être également utilisées pour obtenir de meilleures approximations. Comme on l'a vu, les différences du premier ordre approximent les dérivées (premières) à un terme d'ordre h près. Mais la combinaison
est une approximation de f'(x) à un terme d'ordre h2 près ; ceci peut se montrer en développant l'expression précédente par la formule de Taylor, ou en utilisant le calcul symbolique des différences finies, qui sera exposé plus bas.
Ces différences itérées possèdent des propriétés formelles utiles :
- Pour tous k et n entiers positifs
Il est enfin possible d'utiliser les différences itérées pour approximer des dérivées de tout ordre, et même d'employer un nombre de termes dissymétrique autour de la valeur centrale, ce qui est utile lorsqu'on effectue des calculs sur une grille, pour compenser les effets de bord ; les détails sont analysés dans cette série de notes (en).
Méthodes des différences finies
Article détaillé : méthode des différences finies.En analyse numérique, une importante application des différences finies est la résolution numérique des équations différentielles et des équations aux dérivées partielles : l'idée est de remplacer les dérivées apparaissant dans l'équation par des différences finies qui les approximent. Les diverses méthodes qui en résultent sont appelées méthodes des différences finies.
De nombreuses applications de ces méthodes se rencontrent en théorie du calcul numérique, ainsi que dans des disciplines scientifiques très variées, telles que la mécanique des fluides, et ses applications à la météorologie ou à l'aérodynamique.
Le calcul symbolique des différences finies
L'opérateur de différence avant envoie la fonction f vers Δh[f]. Cet opérateur vérifie Δh = Th − I, où Th est l'opérateur de translation de pas h, défini par Th[f](x) = f(x + h), et I est l'identité (I[f] = f, et I = T0). Δh est un opérateur linéaire, et il satisfait à la règle du produit (on dit parfois que c'est une dérivation; voir aussi à ce sujet l'article corps différentiel). Les opérateurs de différence arrière et centrale obéissent à des règles similaires.
On a
![\Delta^n_h = [T_h-I]^n](9/b39cbedb0e2eaf9401098cac66733cf3.png)
, ![\nabla^n_h = [I-T_{-h}]^n](d/9ed8321dfd1947dea1f133295c8578b8.png)
, et ![\delta^n_h = [T_{h/2}-T_{-h/2}]^n](c/44c38139f209bc33511d3d7e628f729f.png)
, expressions qui peuvent se développer par la formule du binôme (appliquée à l'algèbre des opérateurs), puisque Th et Tk commutent. On vérifie alors que les formules données plus haut pour les différences d'ordre supérieur coïncident avec ces développements, ce qui justifie les notations choisies.Appliquant le théorème de Taylor (par rapport à h), on obtient
où D désigne l'opérateur de dérivation usuel, envoyant f vers sa dérivée f'. Inversant (formellement) l'exponentielle, on peut donc penser que
Cette formule est correcte quand on l'applique à une fonction f polynômiale, car alors la série se réduit à une somme finie. En revanche, même si f est analytique, il n'est pas garanti que la série converge : ce pourrait n'être qu'un développement asymptotique. Elle peut cependant être utilisée pour obtenir des approximations plus précises de la dérivée de f. Ainsi, les deux premiers termes de la série donnent l'approximation du second ordre à f'(x) mentionnée à la fin de la section sur les différences d'ordre supérieur.
On a des formules analogues pour les différences arrière et centrale :
Toutes ces formules se généralisent encore à celles du calcul ombral et à d'autres formules de l'analyse combinatoire ; ainsi, l'opérateur inverse de l'opérateur de différence avant est celui de somme indéfinie[1], ce qui a permis à Lagrange une "démonstration" abstraite de la formule d'Euler-Maclaurin, en l'interprétant comme l'inverse symbolique de la formule de Taylor[2].
Généralisations
Une différence finie généralisée est d'habitude définie par
où
est le vecteur de coefficients qui lui est associé. Une différence infinie en est une généralisation supplémentaire, où la somme finie précédente est remplacée par une série (infinie). Une autre sorte de généralisation consiste à faire dépendre les coefficients μk du point x : μk = μk(x), on parle alors de différences finies pondérées ; on peut aussi faire dépendre le pas h du point x : h = h(x). Toutes ces généralisations sont utiles pour construire différents modules de continuité[3].Différences finies à plusieurs variables
Il est possible de calculer des différences finies dans le cas de fonctions à plusieurs variables ; elles sont l'analogue discret des dérivées partielles. On obtient ainsi les approximations suivantes :
Notes
- Voir pour l'instant l'article indefinite sum(en).
- Joseph-Louis Lagrange, Œuvres, vol 3, p.441 et suivantes
- Voir pour l'instant l'article modulus of continuity (en).
Références
- William F. Ames, Numerical Methods for Partial Differential Equations, Section 1.6. Academic Press, New York, 1977. ISBN 0-12-056760-1.
- Francis B. Hildebrand, Finite-Difference Equations and Simulations, Section 2.2, Prentice-Hall, Englewood Cliffs, New Jersey, 1968.
- Boole, George, A Treatise On The Calculus of Finite Differences, 2nd ed., Macmillan and Company, 1872. [See also: Dover edition 1960].
- Robert D. Richtmyer et K. W. Morton, Difference Methods for Initial Value Problems, 2nd ed., Wiley, New York, 1967.
Voir aussi
Liens externes
(les liens suivants renvoient pour l'instant vers des sites anglophones)
- Finite Difference Method
- Finite Difference Method for Boundary Value Problems
- Table of useful finite difference formula generated using Mathematica
- Finite Calculus: A Tutorial for Solving Nasty Sums
- (en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Finite difference » (voir la liste des auteurs)
 = f(x) - f(x-h). \](2/c8202ccc4ee962447c3a7463a47ecbff.png)
 = f(x+\tfrac12h)-f(x-\tfrac12h). \](0/03003621f1d1c21db25c422bac0c580e.png)
}{h}.](7/757400f2745c7cc6ba4f4e43baf4cabf.png)
}{h} - f'(x) = O(h) \quad (h \to 0).](5/475d8465a4ea17726430ae3cb71ca792.png)
}{h} - f'(x) = O(h).](5/415fabcde94e8a24b84ab476254fe50d.png)
}{h} - f'(x) = O(h^{2}) . \!](4/7a45867fa8cb69b0122039189f0a581c.png)
}{h^2} = \frac{f(x+h) - 2 f(x) + f(x-h)}{h^{2}} .](5/1953a53dd7b94b8f9a1a3dff7d3d322b.png)
 =
\sum_{i = 0}^{n} (-1)^i \binom{n}{i} f(x + (n - i) h),](5/a25c0e912e46484437a236128eda063d.png)
 =
\sum_{i = 0}^{n} (-1)^i \binom{n}{i} f(x - ih),](6/90621614d2f6c04fe52246fc9e8ca7e6.png)
 =
\sum_{i = 0}^{n} (-1)^i \binom{n}{i} f\left(x + \left(\frac{n}{2} - i\right) h\right).](2/682eb520a31cb9c590269e244010b0a9.png)
}{h^n}+O(h)](1/2d120533be03a29b422296ee3487ae36.png)
}{h^n}+O(h)](0/630870d63a734696aac48668362453a4.png)
}{h^n} + O(h^2).](3/393e90ca5911b6e9458ecf95592a4a68.png)
 - \frac12 \Delta_h^2[f](x)}{h} = - \frac{f(x+2h)-4f(x+h)+3f(x)}{2h}](7/fc70d2c77bc379db24104f3ec369d28d.png)
.](2/bc2474a67b708e11398aa2b9052f7524.png)
 \Delta^{n-k}_h[g](x+kh).](b/49b5e98e5b8baad1aa8e8d50f44fef59.png)
 = \sum_{k=0}^N \mu_k f(x+kh),](8/708de88797ca34cc08637cae669d4868.png)
Wikimedia Foundation. 2010.
