Ce document pédagogique propose une série d'exercices corrigés en Analyse Numérique, destinée aux étudiants universitaires. Il explore les techniques fondamentales d'évaluation numérique des intégrales.
Le contenu couvre des notions essentielles de l'intégration numérique, incluant :
- L'application des méthodes des trapèzes, de Simpson et du point milieu.
- Le calcul et l'estimation de l'erreur d'approximation.
- La détermination du nombre de subdivisions nécessaires pour atteindre une précision requise.
- L'impact de la convexité des fonctions sur l'exactitude des méthodes.
Comprendre les Méthodes d'Intégration Numérique : Trapèzes, Simpson et Point Milieu
L'intégration numérique est une branche essentielle de l'analyse numérique, permettant d'estimer la valeur d'intégrales définies lorsque le calcul analytique est difficile ou impossible. Ce guide explore les principes et applications des méthodes des trapèzes, de Simpson et du point milieu, ainsi que l'estimation des erreurs associées.
1. La Méthode des Trapèzes
La méthode des trapèzes approxime l'aire sous une courbe en divisant l'intervalle d'intégration en sous-intervalles et en remplaçant la fonction par un segment de droite sur chacun d'eux, formant ainsi des trapèzes. La formule composite des trapèzes pour une intégrale ∫ab f(x) dx avec un pas h est :
I ≈ h/2 × [f(x0) + 2f(x1) + ... + 2f(xn-1) + f(xn)]
Exemple : Intégrale I = ∫12 (1/x) dx
Considérons l'intégrale I = ∫12 (1/x) dx. La valeur exacte de cette intégrale est ln(2) ≈ 0.693.
Observation : Surestimation par la méthode des trapèzes
Lorsqu'on évalue I par la méthode des trapèzes, la valeur numérique obtenue est supérieure à ln(2). Cette observation s'explique visuellement : la courbe de la fonction y = 1/x est convexe sur l'intervalle [1, 2]. Cela signifie que les segments de droite qui forment le haut des trapèzes sont situés au-dessus de la courbe. Par conséquent, la somme des aires des trapèzes sera toujours supérieure à l'aire exacte sous la courbe.
Généralisation de la Surestimation
Cette surestimation est systématique et demeure vraie quel que soit le pas h choisi, tant que la fonction est convexe sur l'intervalle d'intégration. Une fonction est convexe si sa dérivée seconde est positive (f''(x) > 0). Pour toute fonction f(x) dont la dérivée seconde est positive sur [a,b], la méthode des trapèzes surestimera toujours l'intégrale.
2. La Méthode de Simpson
La méthode de Simpson est une méthode d'intégration numérique qui offre généralement une meilleure précision que celle des trapèzes. Au lieu d'approximer la fonction par des segments de droite, elle utilise des polynômes de degré 2 (paraboles) sur des paires de sous-intervalles. Le nombre de subdivisions 'n' doit être pair.
La formule composite de Simpson pour une intégrale ∫ab f(x) dx avec un pas h est :
I ≈ h/3 × [f(x0) + 4f(x1) + 2f(x2) + 4f(x3) + ... + 2f(xn-2) + 4f(xn-1) + f(xn)]
Détermination du pas pour une erreur donnée
Pour évaluer une intégrale avec la méthode de Simpson et garantir une erreur inférieure à une certaine tolérance (par exemple, 10-4 ou 10-5), le choix du pas h (et donc du nombre de subdivisions n) est crucial. Le pas h est déterminé en fonction de la borne supérieure de la dérivée quatrième de la fonction sur l'intervalle d'intégration. Plus la précision souhaitée est élevée, plus le nombre de subdivisions 'n' doit être grand.
3. La Méthode du Point Milieu
La méthode du point milieu (ou méthode des rectangles au point milieu) approxime l'aire sous la courbe en utilisant des rectangles dont la hauteur est la valeur de la fonction au milieu de chaque sous-intervalle.
La formule composite du point milieu pour une intégrale ∫ab f(x) dx avec un pas h est :
I ≈ h × [f(x1/2) + f(x3/2) + ... + f(xn-1/2)]
Calcul du nombre de subdivisions pour une précision donnée
Pour évaluer des intégrales telles que F(1) où F(x) = ∫0x te-t dt avec une précision de 10-5 en utilisant la méthode du point milieu, il est nécessaire de calculer le nombre de subdivisions. Ce calcul dépend de la borne supérieure de la dérivée seconde de la fonction à intégrer sur l'intervalle.
Formules Composites et Ordre des Erreurs
La précision d'une méthode d'intégration numérique est caractérisée par l'ordre de son erreur. Les formules composites sont plus efficaces pour des intervalles plus grands en divisant le problème en sous-intervalles plus petits.
- Méthode du Point Milieu (Composite) : L'erreur est de l'ordre de O(h²) ou O(1/n²).
- Méthode des Trapèzes (Composite) : L'erreur est également de l'ordre de O(h²) ou O(1/n²).
- Méthode de Simpson (Composite) : L'erreur est de l'ordre de O(h4) ou O(1/n4).
Où h est le pas de subdivision et n est le nombre de subdivisions. La méthode de Simpson converge donc plus rapidement vers la valeur exacte pour des valeurs de h suffisamment petites.
Application : Évaluation de ∫01 t cos(t) dt avec une erreur inférieure à 10-3
Pour évaluer une intégrale comme ∫01 t cos(t) dt avec une précision de 10-3, le nombre de subdivisions requis dépendra de la méthode choisie :
- Méthode des Trapèzes : Nécessitera un nombre de subdivisions calculé en utilisant la borne supérieure de la dérivée seconde de f(t) = t cos(t) sur [0, 1].
- Méthode de Simpson : Nécessitera un nombre de subdivisions, généralement inférieur à celui des trapèzes pour la même précision, en utilisant la borne supérieure de la dérivée quatrième de f(t) = t cos(t) sur [0, 1].
FAQ - Questions Fréquentes sur l'Intégration Numérique
- Quelle est la différence fondamentale entre la méthode des trapèzes et celle de Simpson ?
- La méthode des trapèzes approxime la fonction par des segments de droite sur chaque sous-intervalle. La méthode de Simpson utilise des paraboles (polynômes de degré 2) pour une meilleure approximation de la courbe, ce qui lui confère généralement une plus grande précision pour un même nombre de subdivisions.
- Quand la méthode des trapèzes surestime-t-elle la valeur exacte d'une intégrale ?
- La méthode des trapèzes surestime la valeur exacte d'une intégrale lorsque la fonction est convexe sur l'intervalle d'intégration, c'est-à-dire quand sa dérivée seconde est positive (f''(x) > 0). Dans ce cas, les segments de droite reliant les points de la courbe sont toujours au-dessus de la courbe réelle.
- Pourquoi la méthode de Simpson est-elle souvent préférée pour des calculs de haute précision ?
- La méthode de Simpson possède un ordre de convergence plus élevé (O(h4)) comparée aux méthodes des trapèzes et du point milieu (O(h²)). Cela signifie qu'elle atteint une précision donnée avec un nombre de subdivisions (n) bien plus faible, réduisant ainsi le temps de calcul nécessaire.