Intégration numérique

Calcul numérique d'une intégrale

En analyse numérique, il existe toute une famille d'algorithmes permettant d'approcher la valeur numérique d'une intégrale. Toutes consistent à approcher l'intégrale $I = \int f(x) \, dx$ par une formule dite de quadrature, du type $I(f) = \sum_{i=0}^p \omega_i f(x_i)$ . Le choix de p, des pondérations $ω i$ et des nœuds $x i$ dépendent de la méthode employée. Il conviendra aussi de s'intéresser à la précision des formules utilisées.

Une première indication de l'efficacité d'une méthode est donnée par son ordre. Par définition, une méthode de quadrature est d'ordre n quand elle donne la valeur exacte de l'intégrale pour tout polynôme de degré inférieur ou égal à n, et un résultat faux pour au moins un polynôme de degré n+1.

Sommaire

1 Méthode de calcul d'intégrale à une dimension
- 1.1 Formules simples
- 1.2 Formules composites
2 Autres méthodes de quadrature numérique
- 2.1 Méthode de calcul d'intégrale à plusieurs dimensions
- 2.2 Méthode de calcul d'intégrale de forme particulière
3 Voir aussi

Méthode de calcul d'intégrale à une dimension

Ces méthodes utilisent l'interpolation des fonctions à intégrer. Généralement, les fonctions sont interpolées par des polynômes dont on connait facilement la primitive.

Formules simples

Formules du rectangle et du point milieu

C'est la méthode la plus simple qui consiste à interpoler la fonction $f$ à intégrer par une fonction constante (polynôme de degré 0). Soit $ξ$ le point d'interpolation; la formule devient alors:

$I(f) = (b-a) f(\xi)\,$

Le choix du point a de l'importance pour la détermination du terme d'erreur $E(f) = I - I(f)\,$ :

Si $\xi = a\,$ ou $\xi = b\,$ , l'erreur est $E(f) = \frac{(b-a)^2}{2} f'(\eta), \quad \eta \in [a,b]$ . C'est la méthode du rectangle;
Si $\xi = (a+b)/2\,$ , alors l'erreur devient $E(f) = \frac{(b-a)^3}{24} f''(\eta), \quad \eta \in [a,b]\,$ . Il s'agit de la méthode du point milieu.

Ainsi, le choix du point milieu améliore l'ordre de la méthode. La méthode du rectangle est exacte (c'est-à-dire $E (f) = 0$ ) pour les fonctions constantes et celle du point milieu est exacte pour les polynômes de degré inférieur ou égal à 1. Ceci s'explique par le fait que pour l'intégration de $x$ , la méthode du point milieu donne lieu à deux erreurs d'évaluation, égales en valeur absolue et opposées en signe.

Formule du trapèze

Si on interpole f par un polynôme de degré un (fonction affine), on a besoin de deux points d'interpolation, à savoir (a, f(a)) et (b, f(b)). L'intégrale est alors approchée par l'aire du polynôme interpolateur, en l'occurrence un trapèze. Ceci justifie le nom de méthode des trapèzes:

$I(f) = (b-a) \, \frac{f(a) + f(b)}{2}$

L'erreur commise est

$E(f) = - \frac{(b-a)^3}{12} f''(\eta), \quad \eta \in [a,b]$

L'erreur s'annule pour tout polynôme de degré inférieur ou égal à un. Selon ce critère, la méthode des trapèzes est donc moins performante que celle du point milieu, étant donné que les degrés d'exactitude sont les mêmes et qu'il y a une évaluation de plus avec la méthode des trapèzes.

Formule de Simpson

La fonction f est maintenant remplacée par une parabole, qui nécessite trois points d'interpolation. Les extrémités a, b, et leur milieu m sont choisis. La méthode de Simpson consiste alors à remplacer l'intégrale par

$I(f) = \frac{(b-a)}{6} \left[ f(a) + 4 f(m) + f(b) \right]$

L'erreur est

$E(f) = - \frac{(b-a)^5}{2880} f^{(4)}(\eta), \quad \eta \in [a,b]$

Le degré d'exactitude est de 3 pour cette méthode, pour 3 évaluations de f.

Le tableau suivant résume les performances de chaque méthode

Nom de la Méthode	Degré du polynome	Nombre de points	Degré d’exactitude	Degré d’erreur	Coefficient d’erreur
Rectangle	0	1	0	2	1/2
Point Milieu	0	1	1	3	1/24
Trapèze	1	2	1	3	-1/12
Simpson	2	3	3	5	-1/2880

Les formules de Newton-Cotes permettent de généraliser ces résultats sur des intervalles constants, où la fonction f est interpolée par des polynômes de degré de plus en plus élevé. Pour des questions de stabilité numérique, il est préferable de limiter le degré du polynôme d'interpolation en subdivisant l'intervalle en sous-intervalles, pour lesquels une interpolation linéaire est suffisante.

Formules composites

Pour chacune des méthodes précédentes, le terme d'erreur dépend de b-a. Si cette amplitude est trop élevée, on peut réduire simplement l'erreur en découpant l'intervalle [a, b] en n sous-intervalles, sur lesquels on calculera la valeur approchée de l'intégrale. On parle alors de formule composite. La valeur sur l'intervalle [a, b] sera la somme de la valeur sur chaque sous-intervalle.

Pour la méthode du point milieu, la formule devient

$I(f) = \frac{(b-a)}{n} \sum_{k=0}^{n-1} f( m_k )$

où $m k$ est le milieu du k-ième sous-intervalle. Puisque les n sous-intervalles sont identiques, ils sont de la forme [a + k h, a + (k+1) h], avec h = (b-a)/n et k = 0, 1, 2, ..., n-1. Ceci entraîne finalement que $m k = a + k h + h / 2$ . Le terme d'erreur s'écrit $E(h)= h^2 \frac{(b-a)}{24} f''(\eta),\quad \eta\in [a,b]$ . La formule composite a un ordre 1, comme précédemment. L'agrégation a fait baisser d'une puissance le terme en (a-b).

Pour la méthode des trapèzes, la formule composite est

$I(f) = \frac{(b-a)}{n} \left( {f(a) + f(b) \over 2} + \sum_{k=1}^{n-1} f \left( a + k h \right) \right)$

Le terme d'erreur s'écrit $E(f)=-h^2 \frac{(b-a)}{12} f''(\eta),\quad \eta\in [a,b]$ .

La formule composite de Simpson prend la forme

$I(f) = \frac{h}{6} \left[f(a)+f(b) + 2 \sum_{i=1}^{n-1} f(x_i) + 4 \sum_{i=0}^{n-1} f(x_i + h/2) \right]$

et l'erreur devient

$E(f) = -h^4 \frac{(b-a)}{2880} f^{(4)}(\eta),\quad \eta\in [a,b]$

Autres méthodes de quadrature numérique

Méthode de calcul d'intégrale à plusieurs dimensions

Méthode de Monte-Carlo

Méthode de calcul d'intégrale de forme particulière

Méthode de Laplace pour les intégrales du type $\int_a^b\! e^{M f(x)}\, \mbox{d}x\,$ ;
Méthode du point col pour les intégrales du type $I(\lambda) = \int_\mathcal{C} f(z) e^{\lambda g(z)} \, \mbox{d}z\,$ .

Voir aussi

Corrections:cas d'exactitude

Intégrale
- Table d'intégrales
- Calcul intégral
- Calcul numérique d'une intégrale
Primitive
- Table de primitives
Règles de Bioche

Portail des mathématiques

Ce document provient de « Calcul num%C3%A9rique d%27une int%C3%A9grale ».

Catégories : Algorithme numérique | Intégration numérique

Wikimedia Foundation. 2010.

Contenu soumis à la licence CC-BY-SA. Source : Article Intégration numérique de Wikipédia en français (auteurs)

Игры ⚽ Поможем написать реферат

Regardez d'autres dictionnaires:

intégration numérique — skaitmeninis integravimas statusas T sritis fizika atitikmenys: angl. numerical integration vok. numerische Integration, f rus. численное интегрирование, n pranc. intégration numérique, f … Fizikos terminų žodynas
Intégration (mathématiques) — Pour les articles homonymes, voir intégration. interprétée comme l’aire sous la courbe de f … Wikipédia en Français
numérique — [ nymerik ] adj. • 1616; du lat. numerus « nombre » 1 ♦ Math. Qui est représenté par un nombre, se fait avec des nombres. Valeur numérique. Calcul numérique. 2 ♦ Qui concerne les nombres arithmétiques. Table numérique : table de correspondance… … Encyclopédie Universelle
NUMÉRIQUE (ANALYSE) — Les problèmes et les méthodes numériques ne délimitent pas un secteur spécifique des mathématiques; ils interviennent en effet non seulement dans les domaines traditionnels (analyse classique et équations fonctionnelles), mais aussi en algèbre,… … Encyclopédie Universelle
Integration par changement de variable — Intégration par changement de variable Pour les articles homonymes, voir Intégration et Changement de variable. En mathématiques, le changement de variable est un procédé qui consiste à remplacer une variable ou même une fonction par une autre… … Wikipédia en Français
Intégration Par Changement De Variable — Pour les articles homonymes, voir Intégration et Changement de variable. En mathématiques, le changement de variable est un procédé qui consiste à remplacer une variable ou même une fonction par une autre fonction de celle ci ou d un autre… … Wikipédia en Français
Intégration par changement de variables — Intégration par changement de variable Pour les articles homonymes, voir Intégration et Changement de variable. En mathématiques, le changement de variable est un procédé qui consiste à remplacer une variable ou même une fonction par une autre… … Wikipédia en Français
Intégration de Verlet — L intégration de Verlet est un schéma d intégration qui permet de calculer la trajectoire de particules en simulation dynamique moléculaire. Cette méthode offre une meilleure stabilité que la plus simple méthode d Euler, de même que d importantes … Wikipédia en Français
Intégration par changement de variable — Pour les articles homonymes, voir Intégration et Changement de variable. En mathématiques, le changement de variable est un procédé qui consiste à remplacer une variable ou même une fonction par une autre fonction de celle ci ou d un autre… … Wikipédia en Français
Integration d'une forme differentielle — Forme différentielle Pour les articles homonymes, voir Forme. En géométrie différentielle, une forme différentielle est la donnée d un champ d applications multilinéaires alternées sur les espaces tangents d une variété différentielle possédant… … Wikipédia en Français

Dictionnaires et Encyclopédies sur 'Academic'

Intégration numérique

Calcul numérique d'une intégrale

Sommaire