🔴  Lives #BAC2024

À partir du 12 mai, révise le bac avec nous sur YouTube tous les soirs à 19h30 ! Découvrir la chaîne →

Intégration <FONT color="red"><b>(il faut avoir étudier au préalable le chapitre Primitives)</b></FONT>

Eléments fondamentaux - Exercice 1

40 min
60
Illustrons, dans un cas très simple, donc très pédagogique, les différentes approches numériques du calcul approchée d'une intégrale.
On note par ff la fonction suivante :
f:x[0;1]exf : x \in [0 \,;\, 1] \longmapsto e^x
On pose alors :
I=01f(x)dxI = \int_0^1 f(x) \, dx
Question 1

Calculer la valeur de II par la méthode usuelle des primitives.

Correction
On sait que la fonction exponentielle est définie sur R\mathbb{R} tout entier, et qu'elle y est infiniment dérivable, et de fait elle y est continue.
En outre, on sait que :
xR,(ex)=ex\forall x \in \mathbb{R}, \,\, (e^x)' = e^x
Donc :
I=01f(x)dx=[ex]01=e1e0I = \int_0^1 f(x) \, dx = \big[ e^x \big]_0^1 = e^1 - e^0
Finalement :
I=01f(x)dx=e1u.a.I = \int_0^1 f(x) \, dx = e - 1 \,\, u.a.
De manière approchée, on a :
I=01f(x)dx1,71828u.a.I = \int_0^1 f(x) \, dx \simeq 1,71828 \,\, u.a.
Graphiquement, cela nous donne :
Question 2

Calculer la valeur de II par une somme de RiemannRiemann.

Correction
On sait que :
limn+(1ni=1n1f(a+i×ban))=1baabf(x)dx\lim_{n \longrightarrow + \infty} \left( \dfrac{1}{n} \sum_{i = 1}^{n-1} f\left( a + i \times \dfrac{b-a}{n} \right) \right) = \dfrac{1}{b-a} \int_a^b f(x) \, dx
Dans notre cas, ceci nous donne :
limn+(1ni=1n1e(0+i×10n))=11001f(x)dx\lim_{n \longrightarrow + \infty} \left( \dfrac{1}{n} \sum_{i = 1}^{n-1} e^{\left( 0 + i \times \frac{1-0}{n} \right)} \right)= \dfrac{1}{1-0} \int_0^1 f(x) \, dx
Soit :
I=limn+(1ni=1n1ein)I = \lim_{n \longrightarrow + \infty} \left( \dfrac{1}{n} \sum_{i = 1}^{n-1} e^{\frac{i}{n} } \right)
Et on a :
i=1n1ein=i=1n1(e1n)i=(e1n)1+(e1n)2++(e1n)n1\sum_{i = 1}^{n-1} e^{\frac{i}{n} } = \sum_{i = 1}^{n-1} \left( e^{\frac{1}{n} } \right)^i = \left( e^{\frac{1}{n}} \right)^1 + \left( e^{\frac{1}{n}} \right)^2 + \cdots + \left( e^{\frac{1}{n}} \right)^{n-1}
Il s'agit de la somme des n1n-1 premiers termes d'une suite géométrique de premier terme e1ne^{\frac{1}{n}} et de raison e1ne^{\frac{1}{n}}. On a alors :
i=1n1ein=e1n1(e1n)n11e1n=e1n1en1n1e1n=e1n1e11n1e1n=e1n1e11n1e1n=e1n1ee1n1e1n\sum_{i = 1}^{n-1} e^{\frac{i}{n} } = e^{\frac{1}{n}} \dfrac{1 - \left( e^{\frac{1}{n}} \right)^{n-1}}{1 - e^{\frac{1}{n}}} = e^{\frac{1}{n}} \dfrac{1 - e^{\frac{n-1}{n}} }{1 - e^{\frac{1}{n}}} = e^{\frac{1}{n}} \dfrac{1 - e^{1 - \frac{1}{n}} }{1 - e^{\frac{1}{n}}} = e^{\frac{1}{n}} \dfrac{1 - e^{1 - \frac{1}{n}} }{1 - e^{\frac{1}{n}}} = e^{\frac{1}{n}} \dfrac{1 - ee^{- \frac{1}{n}} }{1 -e^{\frac{1}{n}} }
Soit :
i=1n1ein=e1nee1ne1n1e1n=e1nee1n1n1e1n=e1nee01e1n=e1ne1e1n=ee1ne1n1\sum_{i = 1}^{n-1} e^{\frac{i}{n} } = \dfrac{e^{\frac{1}{n}} - ee^{\frac{1}{n}}e^{- \frac{1}{n}} }{1 - e^{\frac{1}{n}} } = \dfrac{e^{\frac{1}{n}} - ee^{\frac{1}{n}- \frac{1}{n}} }{1 - e^{\frac{1}{n}} } = \dfrac{e^{\frac{1}{n}} - ee^{0} }{1 - e^{\frac{1}{n}} } = \dfrac{e^{\frac{1}{n}} - e }{1 - e^{\frac{1}{n}} } = \dfrac{e - e^{\frac{1}{n}} }{e^{\frac{1}{n}} -1}
Ceci nous permet d'écrire que :
limn+(1ni=1n1ein)=limn+1nee1ne1n1\lim_{n \longrightarrow + \infty} \left( \dfrac{1}{n} \sum_{i = 1}^{n-1} e^{\frac{i}{n} } \right) = \lim_{n \longrightarrow + \infty} \dfrac{1}{n} \dfrac{e - e^{\frac{1}{n}} }{e^{\frac{1}{n}} -1}
Posons X=1nX = \dfrac{1}{n}. Si n+n \longrightarrow + \infty alors X0X \longrightarrow 0. Donc :
limn+(1ni=1n1ein)=limn+1nee1ne1n1=limX0XeeXeX1\lim_{n \longrightarrow + \infty} \left( \dfrac{1}{n} \sum_{i = 1}^{n-1} e^{\frac{i}{n} } \right) = \lim_{n \longrightarrow + \infty} \dfrac{1}{n}\dfrac{e - e^{\frac{1}{n}} }{e^{\frac{1}{n}} -1} = \lim_{X \longrightarrow 0} X\dfrac{e - e^X}{e^X -1}
Or, lorsque X0X \longrightarrow 0 on a l'équivalence suivante :
eX01+Xe^X \underset{0}{\sim} 1 + X
Donc :
limX0XeeXeX1=limX0Xe(1+X)1+X1=limX0Xe1XX=limX0(e1X)=e1\lim_{X \longrightarrow 0} X\dfrac{e - e^X}{e^X -1} = \lim_{X \longrightarrow 0} X\dfrac{e - (1 + X)}{1 + X - 1} = \lim_{X \longrightarrow 0} X\dfrac{e - 1 - X}{ X } = \lim_{X \longrightarrow 0}(e - 1 - X) = e - 1
Finalement, on a bien :
I=e1u.a.I = e - 1 \,\, u.a.
Question 3

Calculer, avec une précision meilleure que le millième, la valeur approchée de II par la méthode numérique des rectangles.

Correction
On sait que la fonction exponentielle est définie sur R\mathbb{R} tout entier, et qu'elle y est infiniment dérivable, et de fait elle y est continue.
En outre, on sait que :
xR,(ex)=ex\forall x \in \mathbb{R}, \,\, (e^x)' = e^x
Donc f(x)=exf'(x) = e^x. Ce qui implique que ff' est croissante sur l'intervalle d'intégration [0;1][0\,;\,1], et de fait :
supx[0;1]f(x)=e1=e\underset{x \in [0\,;\,1]}{\sup}\left| f'(x) \right| = e^1 = e
De plus, on sait que :
IRnsupx[a;b]f(x)×(ba)22n\big| I - \mathcal{R}_n \big| \leqslant \underset{x \in [a\,;\,b]}{\sup}\left| f'(x) \right| \times\dfrac{(b-a)^2}{2n}
Donc, pour notre situation :
supx[0;1]f(x)×(10)22n=e2n\underset{x \in [0\,;\,1]}{\sup}\left| f'(x) \right| \times\dfrac{(1-0)^2}{2n} = \dfrac{e}{2n}
Comme on désire la valeur approchée de II, meilleure que le millième, on a donc :
e2n<110001000e2<n500e<n\dfrac{e}{2n} < \dfrac{1}{1000} \,\,\, \Longleftrightarrow \,\,\, \dfrac{1000e}{2} < n \,\,\, \Longleftrightarrow \,\,\, 500e < n
Ce qui nous permet de choisir n=1360n = 1360 car nNn \in \mathbb{N}^\star. Et on a alors :
h=101360=11360h = \dfrac{1-0}{1360} = \dfrac{1}{1360}
On a alors :
Rn=hk=0n1f(0+kh)=11360k=013601f(k1360)=11360k=01359f(k1360)=11360k=01359e(k1360)=11360k=01359(e11360)k\mathcal{R}_n = h \sum_{k=0}^{n-1} f(0 + kh) = \dfrac{1}{1360} \sum_{k=0}^{1360-1} f\left (\dfrac{k}{1360} \right) = \dfrac{1}{1360} \sum_{k=0}^{1359} f\left (\dfrac{k}{1360} \right) = \dfrac{1}{1360} \sum_{k=0}^{1359} e^{\left (\frac{k}{1360} \right)} = \dfrac{1}{1360} \sum_{k=0}^{1359} \left( e^{ \frac{1}{1360}} \right)^k
Soit :
Rn=11360((e11360)0+(e11360)1++(e11360)1359)\mathcal{R}_n = \dfrac{1}{1360} \left( \left( e^{ \frac{1}{1360}} \right)^0 + \left( e^{ \frac{1}{1360}} \right)^1 + \cdots + \left( e^{ \frac{1}{1360}} \right)^{1359} \right)
On constate la présence de la somme des 13601360 premiers termes d'une suite géométrique de premier terme (e11360)0=1\left( e^{ \frac{1}{1360}} \right)^0 = 1 et de raison e11360e^{\frac{1}{1360}}. On a alors :
Rn=11360(1×1(e11360)13601e11360)=11360(1×1e136013601e11360)=11360(1×1e11e11360)=11360×1e1e11360\mathcal{R}_n = \dfrac{1}{1360} \left( 1 \times \dfrac{1 - \left( e^{ \frac{1}{1360}} \right)^{1360} }{ 1 - e^{\frac{1}{1360}} } \right) = \dfrac{1}{1360} \left( 1 \times \dfrac{1 - e^{ \frac{1360}{1360}} }{ 1 - e^{\frac{1}{1360}} } \right) = \dfrac{1}{1360} \left( 1 \times \dfrac{1 - e^{1} }{ 1 - e^{\frac{1}{1360}} } \right) = \dfrac{1}{1360} \times \dfrac{1 - e }{ 1 - e^{\frac{1}{1360}} }
Soit :
Rn=11360×e1e113601u.a.\mathcal{R}_n = \dfrac{1}{1360} \times \dfrac{e - 1}{ e^{\frac{1}{1360}} - 1 } \,\, u.a.
De manière approchée, on a :
Rn1,71765u.a.\mathcal{R}_n \simeq 1,71765 \,\, u.a.
Question 4

Calculer, avec une précision meilleure que le millième, la valeur approchée de II par la méthode numérique des trapèzes.

Correction
On sait que la fonction exponentielle est définie sur R\mathbb{R} tout entier, et qu'elle y est infiniment dérivable, et de fait elle y est continue.
En outre, on sait que :
xR,(ex)=ex\forall x \in \mathbb{R}, \,\, (e^x)'' = e^x
Donc f(x)=exf''(x) = e^x. Ce qui implique que ff'' est croissante sur l'intervalle d'intégration [0;1][0\,;\,1], et de fait :
supx[0;1]f(x)=e1=e\underset{x \in [0\,;\,1]}{\sup}\left| f''(x) \right| = e^1 = e
De plus, on sait que :
ITnsupx[a;b]f(x)×(ba)312n2\big| I - \mathcal{T}_n \big| \leqslant \underset{x \in [a\,;\,b]}{\sup}\left| f''(x) \right| \times\dfrac{(b-a)^3}{12n^2}
Donc, pour notre situation :
supx[0;1]f(x)×(10)312n2=e12n2\underset{x \in [0\,;\,1]}{\sup}\left| f''(x) \right| \times\dfrac{(1-0)^3}{12n^2} = \dfrac{e}{12n^2}
Comme on désire la valeur approchée de II, meilleure que le millième, on a donc :
e12n2<110001000e12<n2250e12<n2\dfrac{e}{12n^2} < \dfrac{1}{1000} \,\,\, \Longleftrightarrow \,\,\, \dfrac{1000e}{12} < n^2 \,\,\, \Longleftrightarrow \,\,\, \dfrac{250e}{12} < n^2
Soit n2>226,52n^2 > 226,52. Ce qui nous permet de choisir n=16n = 16 car nNn \in \mathbb{N}^\star. Et on a alors :
h=1016=116h = \dfrac{1-0}{16} = \dfrac{1}{16}
On a alors :
Tn=116(e0+e12+k=1161e(0+k16))=116(1+e2+k=115ek16)=116(1+e2+k=115(e116)k)\mathcal{T}_n = \dfrac{1}{16} \left( \dfrac{e^0 + e^1}{2} + \sum_{k=1}^{16-1} e^{\left(0 + \frac{k}{16} \right)} \right) = \dfrac{1}{16} \left( \dfrac{1 + e}{2} + \sum_{k=1}^{15} e^{\frac{k}{16}} \right) = \dfrac{1}{16} \left( \dfrac{1 + e}{2} + \sum_{k=1}^{15} \left( e^{\frac{1}{16}} \right)^k\right)
Avec :
k=115(e116)k=(e116)1+(e116)1++(e116)15\sum_{k=1}^{15} \left( e^{\frac{1}{16}} \right)^k = \left( e^{\frac{1}{16}} \right)^1 + \left( e^{\frac{1}{16}} \right)^1 + \cdots + \left( e^{\frac{1}{16}} \right)^{15}
Il s'agit de la somme des 1515 premiers termes d'une suite géométrique de premier terme e116e^{ \frac{1}{16}} et de raison e116e^{\frac{1}{16}}. On a alors :
k=115(e116)k=e116×1(e116)151e116=e116×1e15161e116=e116×e15161e1161=e116e1516e116e1161=e1616e116e1161\sum_{k=1}^{15} \left( e^{\frac{1}{16}} \right)^k = e^{\frac{1}{16}} \times \dfrac{1 - \left( e^{ \frac{1}{16}} \right)^{15} }{ 1 - e^{\frac{1}{16}} } = e^{\frac{1}{16}} \times \dfrac{1 - e^{ \frac{15}{16}} }{ 1 - e^{\frac{1}{16}} } = e^{\frac{1}{16}} \times \dfrac{e^{ \frac{15}{16}} - 1}{ e^{\frac{1}{16}} - 1} = \dfrac{e^{\frac{1}{16}}e^{ \frac{15}{16}} - e^{\frac{1}{16}}}{ e^{\frac{1}{16}} - 1} = \dfrac{e^{\frac{16}{16}} - e^{\frac{1}{16}}}{ e^{\frac{1}{16}} - 1}
Donc :
k=115(e116)k=e1e116e1161=ee116e1161\sum_{k=1}^{15} \left( e^{\frac{1}{16}} \right)^k = \dfrac{e^{1} - e^{\frac{1}{16}}}{ e^{\frac{1}{16}} - 1} = \dfrac{e - e^{\frac{1}{16}}}{ e^{\frac{1}{16}} - 1}
Ainsi :
Tn=116(1+e2+ee116e1161)u.a.\mathcal{T}_n = \dfrac{1}{16} \left( \dfrac{1 + e}{2} + \dfrac{e - e^{\frac{1}{16}} }{ e^{\frac{1}{16}} - 1} \right) \,\, u.a.
De manière approchée :
Tn1,71884u.a.\mathcal{T}_n \simeq 1,71884 \,\, u.a.
Question 5

Calculer, avec une précision meilleure que le millième, la valeur approchée de II par la méthode numérique de SimpsonSimpson.

Correction
On sait que la fonction exponentielle est définie sur R\mathbb{R} tout entier, et qu'elle y est infiniment dérivable, et de fait elle y est continue.
En outre, on sait que :
xR,(ex)=ex\forall x \in \mathbb{R}, \,\, (e^x)'''' = e^x
Donc f(x)=exf''''(x) = e^x. Ce qui implique que ff'''' est croissante sur l'intervalle d'intégration [0;1][0\,;\,1], et de fait :
supx[0;1]f(x)=e1=e\underset{x \in [0\,;\,1]}{\sup}\left| f''''(x) \right| = e^1 = e
De plus, on sait que :
IS2nsupx[a;b]f(x)×(ba)5180(2n)4\big| I - \mathcal{S}_{2n} \big| \leqslant \underset{x \in [a\,;\,b]}{\sup}\left| f''''(x) \right| \times\dfrac{(b-a)^5}{180(2n)^4}
Donc, pour notre situation :
supx[0;1]f(x)×(10)5180(2n)4=e180×24×n4=e180×16×n4=e2880n4\underset{x \in [0\,;\,1]}{\sup}\left| f''''(x) \right| \times\dfrac{(1-0)^5}{180(2n)^4} = \dfrac{e}{180 \times 2^4 \times n^4} = \dfrac{e}{180 \times 16 \times n^4} = \dfrac{e}{2880 n^4}
Comme on désire la valeur approchée de II, meilleure que le millième, on a donc :
e2880n4<110001000e2880<n425e72<n4\dfrac{e}{2880 n^4} < \dfrac{1}{1000} \,\,\, \Longleftrightarrow \,\,\, \dfrac{1000e}{2880} < n^4 \,\,\, \Longleftrightarrow \,\,\, \dfrac{25e}{72} < n^4
Soit n4>0,94n^4 > 0,94, d'où n>0,98n > 0,98. Ce qui nous permet de choisir n=1n = 1 car nNn \in \mathbb{N}^\star. Et on a alors :
h=102×1=12h = \dfrac{1-0}{2 \times 1} = \dfrac{1}{2}
De plus, on sait que, de manière générale, on a :
S2n=h3(f(a)+f(b)+4k=0n1f(a+(2k+1)h)+2k=1n1f(a+(2k)h))\mathcal{S}_{2n} = \dfrac{h}{3} \left( f(a) + f(b) + 4 \sum_{k=0}^{n-1} f\big(a + (2k+1)h \big) + 2 \sum_{k=1}^{n-1} f\big(a + (2k)h \big) \right)
Dans notre cas à nous (n=1)(n=1), on obtient donc :
S2n=123(e0+e1+4k=011e(0+(2k+1)12)+2k=111e(0+(2k)12))=16(1+e+4k=00e12(2k+1)+2k=10ek)\mathcal{S}_{2n} = \dfrac{\dfrac{1}{2}}{3} \left( e^0 + e^1 + 4 \sum_{k=0}^{1-1} e^{\big(0 + (2k+1)\frac{1}{2} \big)} + 2 \sum_{k=1}^{1-1} e^{\big(0 + (2k)\frac{1}{2} \big)} \right) = \dfrac{1}{6} \left( 1 + e + 4 \sum_{k=0}^{0} e^{\frac{1}{2}(2k+1)} + 2 \sum_{k=1}^{0} e^{k} \right)
On constate que la dernière somme k=10ek\sum_{k=1}^{0} e^{k} n'existe pas. De fait, on obtient :
S2n=16(1+e+4k=00e12(2k+1))\mathcal{S}_{2n} = \dfrac{1}{6} \left( 1 + e + 4 \sum_{k=0}^{0} e^{\frac{1}{2}(2k+1)} \right)
Puis, on constate également que :
k=00e12(2k+1)=e12(2×0+1)=e12(0+1)=e12=e\sum_{k=0}^{0} e^{\frac{1}{2}(2k+1)} = e^{\frac{1}{2}(2\times0+1)} = e^{\frac{1}{2}(0+1)} = e^{\frac{1}{2}} = \sqrt{e}
Ainsi, on trouve que :
S2n=16(1+e+4e)u.a.\mathcal{S}_{2n} = \dfrac{1}{6} \left( 1 + e + 4 \sqrt{e} \right) \,\, u.a.
De manière approchée, on obtient :
S2n1,71886u.a.\mathcal{S}_{2n} \simeq 1,71886 \,\, u.a.
Question 6

Quelle(s) conclusion(s) pouvez-vous faire à partir de cet exemple simple ?

Correction
On constate que :
- la méthode de SimpsonSimpson est celle qui est la plus rapide à mettre oeuvre ;
- la méthode des rectangles est celle qui est la plus intuitive mais la moins rapide à mettre oeuvre ;
- La méthode par les sommes de RiemannRiemann s'avère plus technique ;
- toutes ces approches sont parfaitement cohérentes entre-elles.

Il existe d'autres méthode, citons par exemple :
- la méthode des tangentes ;
- la méthode de BooleVillarceauBoole-Villarceau ;
- la méthode de WeddleHardyWeddle-Hardy ;
- la méthode de NewtonCotesNewton-Cotes ;
- la méthode de RombergRomberg.