Analyse · Somme de Riemann

Somme de Riemann

Une somme de Riemann additionne f(xᵢ*)·Δx sur les morceaux de [a, b] : les aires algébriques de rectangles dont les hauteurs sont des valeurs échantillonnées de f. À mesure que les morceaux s’affinent, les sommes d’une fonction intégrable tendent vers l’intégrale définie.

01

La somme de Riemann dans cette visualisation

Une somme de Riemann approche l’intégrale de ff sur [a,b][a, b]. Découpez l’intervalle en sous-intervalles [xi−1,xi][x_{i-1}, x_i], choisissez un point xi∗x_i^* dans chacun, puis additionnez f(xi∗) Δxif(x_i^*)\,\Delta x_i. Cette page utilise nn sous-intervalles égaux, donc chaque morceau a la même largeur :

Somme de Riemann
Sn=∑i=1nf(xi∗) Δx,Δx=b−an.S_n = \sum_{i=1}^{n} f(x_i^*)\,\Delta x, \qquad \Delta x = \frac{b-a}{n}.

La définition générale autorise des sous-intervalles de largeurs différentes et un point quelconque dans chacun, puis prend la limite lorsque la largeur maximale, λ=max⁡Δxi\lambda = \max \Delta x_i, tend vers 0 ; avec des sous-intervalles égaux comme ici, λ→0\lambda \to 0 équivaut à n→∞n \to \infty.

Chaque rectangle jaune repose sur un sous-intervalle et atteint f(xi∗)f(x_i^*), le point jaune de la courbe ; SnS_n additionne leurs aires algébriques : un rectangle sous l’axe des x compte donc négativement, et un rectangle peut dépasser la courbe ou rester en dessous. La méthode détermine la position du point choisi : extrémité gauche xi−1x_{i-1}, extrémité droite xix_i, ou milieu. La somme des trapèzes relie plutôt f(xi−1)f(x_{i-1}) et f(xi)f(x_i) par une corde et égale la moyenne des sommes à gauche et à droite :

Somme des trapèzes
Tn=∑i=1nf(xi−1)+f(xi)2 Δx=Ln+Rn2.T_n = \sum_{i=1}^{n} \frac{f(x_{i-1}) + f(x_i)}{2}\,\Delta x = \frac{L_n + R_n}{2}.

n → ∞ double nn étape par étape jusqu’à n=256n = 256, chaque morceau se divisant donc en deux (depuis un nn qui n’est pas une puissance de deux, la première étape atteint la puissance de deux suivante). Pour une fonction intégrable, les sommes tendent vers l’intégrale définie, affichée en bleu :

Intégrale définie
∫abf(x) dx=lim⁡n→∞Sn=F(b)−F(a).\int_a^b f(x)\,dx = \lim_{n \to \infty} S_n = F(b) - F(a).

La page calcule l’intégrale exactement à partir d’une primitive FF, jamais d’une somme avec un grand nn ; l’erreur est SnS_n moins cette valeur.

02

À observer

  • L’aire est algébrique. Lorsque f<0f < 0, un morceau pend sous l’axe des x ; il est hachuré et compte négativement. Sur [0,2π][0, 2\pi], les deux arches de sin⁡x\sin x se compensent, et l’intégrale vaut 0.
  • Pour une fonction croissante, la somme à gauche est trop petite et celle à droite trop grande ; pour une fonction décroissante, c’est l’inverse. L’intégrale est comprise entre les deux.
  • La rapidité de diminution de l’erreur dépend de la méthode. Pour une fonction régulière et nn grand, les erreurs à gauche et à droite se comportent comme C/nC/n, et celles du milieu et du trapèze comme C/n2C/n^2 : chaque doublement finit donc par les diviser par deux ou par quatre. Au début, le facteur peut différer (1.6 pour les premières sommes à gauche de x2x^2 sur [0,1][0, 1]), et lorsque le terme dominant s’annule, l’erreur diminue plus vite : pour x2x^2 sur [−1,1][-1, 1], les erreurs à gauche et à droite sont déjà en 1/n21/n^2. Le graphique d’erreur dans la barre de valeurs a deux échelles logarithmiques ; le facteur apparaît donc dans la pente de chaque ligne : son logarithme en base 2, pente 1 pour un facteur 2 et pente 2 pour un facteur 4. Gauche et droite descendent doucement, milieu et trapèze environ deux fois plus raide, sauf lorsque le terme dominant s’annule.
  • Les erreurs du milieu et du trapèze ont généralement des signes opposés. Pour f(x)=x2f(x) = x^2, l’erreur au milieu vaut exactement −12-\tfrac12 de l’erreur des trapèzes, pour tout nn.
  • La définition de l’intégrale autorise des morceaux inégaux et un point quelconque dans chacun ; les morceaux égaux avec une règle fixe sont le cas particulier dessiné ici. La limite doit être la même pour tous ces choix lorsque le plus grand morceau tend vers une largeur nulle.
  • Au sens strict, la somme des trapèzes n’est pas une somme de Riemann (pas un seul point choisi par morceau), mais étant la moyenne des sommes à gauche et à droite, elle a la même limite. Les méthodes des rectangles et des trapèzes en intégration numérique sont ces sommes.
03

Commandes

  • Déplacez les lignes pointillées x=ax = a et x=bx = b, ou leurs marques sur le bord inférieur, tournez les deux roues sous Intervalle, ou saisissez les extrémités : π, 2π, π/2 et e sont reconnus. Les touches sous les roues les placent sur des valeurs exactes.
  • Les touches de méthode choisissent la position du point d’échantillonnage ; sous chacune figurent sa somme et son erreur pour le nn actuel, permettant de comparer les quatre d’un coup d’œil. Le graphique entre la somme et l’intégrale montre la valeur absolue de l’erreur selon nn pour les quatre méthodes.
  • Le bouton règle nn à 1, 2, 4 … 256, un cran par doublement ; ses flèches passent au cran suivant. Toute autre valeur de nn jusqu’à 10 000 peut être saisie. n → ∞ double nn étape par étape jusqu’à 256 (d’abord jusqu’à la puissance de deux suivante) et y reste ; un appui à 256 relance depuis 1. Le Mode présentation ouvre une vue en plein écran pour la projection (Espace lance n → ∞, ← → divisent n par deux et le doublent, H masque les réponses).
04

Voir aussi