ATLAS Ouvrir dans l'atlas →

ExerciceMaths Sup – MPSI, chapitre 18 — Intégration · Sommes de Riemann et calcul approché

Exercice 45

ProblèmeDémonstration

Inégalités de Hermite-Hadamard et méthode du point milieu. Soit f:[a ,b]→ℝf : [a\,,b] \to \mathbb{R} et m=a+b2m = \frac{a + b}{2}.

1. On suppose ff continue et convexe. Montrer que f(m)⩽1b−a∫abf⩽f(a)+f(b)2.f(m) \leqslant \frac{1}{b - a}\int_{a}^{b} f \leqslant \frac{f(a) + f(b)}{2}. Pour l'inégalité de gauche, on pourra comparer f(t)f(t) et f(a+b−t)f(a + b - t) à f(m)f(m).

2. On suppose ff de classe 𝒞2\mathcal{C}^{2}, et l'on note M2=max⁡[a,b]∣f′′∣M_{2} = \max\limits_{[a, b]} |f''|. Montrer que ∣∫abf−(b−a)f(m)∣⩽M2 (b−a)324.\Big|\int_{a}^{b} f - (b - a)f(m)\Big| \leqslant M_{2}\,\frac{(b - a)^{3}}{24}.

3. Pour n⩾1n \geqslant 1, on note xk=a+k b−anx_{k} = a + k\,\frac{b - a}{n}. En déduire que ∣∫abf−b−an∑k=0n−1f(xk+xk+12)∣⩽M2 (b−a)324n2.\Big|\int_{a}^{b} f - \frac{b - a}{n}\sum_{k=0}^{n-1} f\Big(\frac{x_{k} + x_{k+1}}{2}\Big)\Big| \leqslant M_{2}\,\frac{(b - a)^{3}}{24n^{2}}. C'est la majoration de l'erreur de la méthode du point milieu, admise dans le cours.

4. Montrer que, si ff est continue et convexe, la méthode du point milieu donne une valeur approchée par défaut de ∫abf\int_{a}^{b} f, et celle des trapèzes une valeur approchée par excès.

5. Pour f(x)=41+x2f(x) = \frac{4}{1 + x^{2}} sur [0 ,1][0\,,1], calculer la valeur approchée de π\pi donnée par la méthode du point milieu avec n=4n = 4, et la comparer à la valeur T4(f)≈3,1312T_{4}(f) \approx 3{,}1312 donnée par la méthode des trapèzes dans le cours.

Corrigé

1. Pour l'inégalité de droite, on utilise la position de la courbe sous ses cordes (chapitre 9, proposition 9). Pour t∈[a ,b]t \in [a\,,b], la corde joignant les points d'abscisses aa et bb donne f(t)⩽b−tb−a f(a)+t−ab−a f(b).f(t) \leqslant \frac{b - t}{b - a}\,f(a) + \frac{t - a}{b - a}\,f(b). On intègre sur [a ,b][a\,,b], avec ∫ab(b−t) dt=∫ab(t−a) dt=(b−a)22\int_{a}^{b} (b - t)\,\mathrm{d} t = \int_{a}^{b} (t - a)\,\mathrm{d} t = \frac{(b - a)^{2}}{2} : ∫abf⩽b−a2(f(a)+f(b)).\int_{a}^{b} f \leqslant \frac{b - a}{2}\big(f(a) + f(b)\big). Pour l'inégalité de gauche, on remarque que mm est le milieu de tt et de a+b−ta + b - t. La convexité, avec les poids 12\frac{1}{2} et 12\frac{1}{2}, donne f(m)⩽f(t)+f(a+b−t)2pour tout t∈[a ,b].f(m) \leqslant \frac{f(t) + f(a + b - t)}{2} \qquad \text{pour tout } t \in [a\,,b]. On intègre sur [a ,b][a\,,b]. Le changement de variable affine u=a+b−tu = a + b - t, qui échange aa et bb, montre que ∫abf(a+b−t) dt=∫abf(u) du\int_{a}^{b} f(a + b - t)\,\mathrm{d} t = \int_{a}^{b} f(u)\,\mathrm{d} u : on obtient (b−a)f(m)⩽∫abf(b - a)f(m) \leqslant \int_{a}^{b} f. f(m)⩽1b−a∫abf⩽f(a)+f(b)2\boxed{f(m) \leqslant \frac{1}{b - a}\int_{a}^{b} f \leqslant \frac{f(a) + f(b)}{2}} Géométriquement, l'aire sous la courbe est comprise entre celle du rectangle de hauteur f(m)f(m), qui est aussi celle du trapèze limité par la tangente en mm quand ff est dérivable, et celle du trapèze limité par la corde.

2. On applique l'inégalité de Taylor-Lagrange à l'ordre 11 (théorème 8) en mm : pour tout t∈[a ,b]t \in [a\,,b], ∣f(t)−f(m)−f′(m)(t−m)∣⩽M2 (t−m)22.\big|f(t) - f(m) - f'(m)(t - m)\big| \leqslant M_{2}\,\frac{(t - m)^{2}}{2}. On intègre sur [a ,b][a\,,b]. Le terme du premier ordre disparaît, car ∫ab(t−m) dt=0\int_{a}^{b} (t - m)\,\mathrm{d} t = 0 par symétrie : c'est tout l'intérêt du point milieu. Par ailleurs, en notant ℓ=b−a2\ell = \frac{b - a}{2}, ∫ab(t−m)22 dt=[(t−m)36]ab=2ℓ36=(b−a)324.\int_{a}^{b} \frac{(t - m)^{2}}{2}\,\mathrm{d} t = \Big[\frac{(t - m)^{3}}{6}\Big]_{a}^{b} = \frac{2\ell^{3}}{6} = \frac{(b - a)^{3}}{24}. Par l'inégalité triangulaire (théorème 3 (iv)), ∣∫abf−(b−a)f(m)∣⩽M2 (b−a)324\boxed{\Big|\int_{a}^{b} f - (b - a)f(m)\Big| \leqslant M_{2}\,\frac{(b - a)^{3}}{24}}

3. On applique la question 2 sur chaque segment [xk ,xk+1][x_{k}\,,x_{k+1}], de longueur h=b−anh = \frac{b - a}{n} et de milieu mk=xk+xk+12m_{k} = \frac{x_{k} + x_{k+1}}{2}, le réel M2M_{2} majorant encore ∣f′′∣|f''| sur ce segment. Par la relation de Chasles et l'inégalité triangulaire, ∣∫abf−h∑k=0n−1f(mk)∣⩽∑k=0n−1∣∫xkxk+1f−hf(mk)∣⩽n M2h324,\Big|\int_{a}^{b} f - h\sum_{k=0}^{n-1} f(m_{k})\Big| \leqslant \sum_{k=0}^{n-1} \Big|\int_{x_{k}}^{x_{k+1}} f - hf(m_{k})\Big| \leqslant n\,\frac{M_{2}h^{3}}{24},

∣∫abf−b−an∑k=0n−1f(xk+xk+12)∣⩽M2 (b−a)324n2\boxed{\Big|\int_{a}^{b} f - \frac{b - a}{n}\sum_{k=0}^{n-1} f\Big(\frac{x_{k} + x_{k+1}}{2}\Big)\Big| \leqslant M_{2}\,\frac{(b - a)^{3}}{24n^{2}}} C'est la majoration admise dans la deuxième remarque qui suit la proposition 14.

4. On suppose ff convexe et, comme à la question 1, continue. On note Mn(f)M_{n}(f) la somme du point milieu. L'inégalité de gauche de la question 1, sur chaque [xk ,xk+1][x_{k}\,,x_{k+1}], donne hf(mk)⩽∫xkxk+1fhf(m_{k}) \leqslant \int_{x_{k}}^{x_{k+1}} f ; celle de droite donne ∫xkxk+1f⩽h f(xk)+f(xk+1)2\int_{x_{k}}^{x_{k+1}} f \leqslant h\,\frac{f(x_{k}) + f(x_{k+1})}{2}. En sommant, f convexe  ⟹  Mn(f)⩽∫abf⩽Tn(f)\boxed{f \text{ convexe} \implies M_{n}(f) \leqslant \int_{a}^{b} f \leqslant T_{n}(f)} La méthode du point milieu donne une valeur par défaut, celle des trapèzes une valeur par excès. On retrouve, sans supposer ff de classe 𝒞2\mathcal{C}^{2}, la première remarque qui suit la proposition 14.

5. Les milieux sont 18\frac{1}{8}, 38\frac{3}{8}, 58\frac{5}{8} et 78\frac{7}{8}, et f(j8)=25664+j2f\big(\frac{j}{8}\big) = \frac{256}{64 + j^{2}}, d'où M4(f)=14∑k=03f(2k+18)=6465+6473+6489+64113,M4(f)≈3,146 801M_{4}(f) = \frac{1}{4}\sum_{k=0}^{3} f\Big(\frac{2k + 1}{8}\Big) = \frac{64}{65} + \frac{64}{73} + \frac{64}{89} + \frac{64}{113}, \qquad \boxed{M_{4}(f) \approx 3{,}146\,801}

M4(f)−π≈+5,2×10−3,T4(f)−π≈3,131 176−π≈−1,04×10−2.M_{4}(f) - \pi \approx +5{,}2 \times 10^{-3}, \qquad T_{4}(f) - \pi \approx 3{,}131\,176 - \pi \approx -1{,}04 \times 10^{-2}. L'erreur du point milieu est environ la moitié de celle des trapèzes, et de signe opposé, ce qu'annonçaient les constantes 124\frac{1}{24} et 112\frac{1}{12}. Avec M2=max⁡[0,1]∣f′′∣=8M_{2} = \max_{[0, 1]} |f''| = 8, la question 3 garantit une erreur au plus 824×16≈2,1×10−2\frac{8}{24 \times 16} \approx 2{,}1 \times 10^{-2}, quatre fois plus que l'erreur réelle.

Piège. La question 4 ne s'applique pas ici : f′′(x)=8(3x2−1)(1+x2)3f''(x) = \frac{8(3x^{2} - 1)}{(1 + x^{2})^{3}} est négative sur [0 ,13]\big[0\,,\frac{1}{\sqrt{3}}\big], et ff n'est pas convexe sur [0 ,1][0\,,1]. Les signes observés, point milieu par excès et trapèzes par défaut, sont ceux d'une fonction concave, ce que ff est sur la plus grande partie de l'intervalle.

À retenir. Les deux erreurs étant à peu près dans le rapport −12-\frac{1}{2}, la combinaison 2M4(f)+T4(f)3\frac{2M_{4}(f) + T_{4}(f)}{3} les compense : elle vaut 3,141 592 503{,}141\,592\,50, à 1,5×10−71{,}5 \times 10^{-7} de π\pi, avec seulement neuf évaluations de ff. C'est la méthode de Simpson, que l'on rencontrera en informatique.

Dans l'édifice

Thèmes

sommes de riemann, intégrales, calcul approché, suites d'intégrales

Notions

somme de Riemann, fonction convexe, dérivées successives

Ouvrir dans l'atlas interactif →