Formule de Newton-Cotes
En analyse numérique, les formules de Newton-Cotes, du nom d'Isaac Newton et de Roger Cotes, servent au calcul numérique d'une intégrale sur un intervalle réel Modèle:Math, ceci à l’aide d’une interpolation polynomiale de la fonction en des points répartis uniformément.
Méthodologie
La fonction Modèle:Math est évaluée en des points équidistants Modèle:Math, pour Modèle:Math et Modèle:Math. La formule de degré Modèle:Math est définie ainsi :
où les Modèle:Math sont appelés les coefficients de quadrature. Ils se déduisent d'une base de polynômes de Lagrange et sont indépendants de la fonction Modèle:Math.
Plus précisément, si Modèle:Math est l'interpolation lagrangienne aux points Modèle:Math et <math>l_i(X) = \prod_{j=0, j\neq i}^n\frac{X-x_j}{x_i-x_j}</math>, alors :
&= \int_a^b \sum_{i=0}^n f(x_i)\, l_i(x)~{\rm d}x\\ &=\sum_{i=0}^n \int_a^b f(x_i) l_i(x)~{\rm d}x\\
&=\sum_{i=0}^n f(x_i) \underbrace{\int_a^b l_i(x)~{\rm d}x}_{w_i}.\end{align}</math>Ainsi ; Modèle:Retrait Le changement de variable <math>y=\frac{x-a}{\Delta}</math> conduit à l'expression<ref name="Wolfram" />: Modèle:Retrait{i!(n-i)!}\int_0^n\prod_{k=0,k\ne i}^n(y-k)~{\rm d}y.</math>}}
Application pour n = 1
En calculant l'expression précédente lorsque Modèle:Math et Modèle:Math, on obtient
- <math>
\begin{align} w_0 &= (b-a) \frac{(-1)^{1-0}}{0!\, (1-0)!} \int_0^1 \prod_{k=0,k \ne 0}^1 (y - k)~{\rm d}y \\
&= -(b-a) \int_0^1(y-1)~{\rm d}y \\
&= -(b-a) \left[ \frac{(y-1)^2}{2} \right]^1_0 \\
&= \frac{b-a}2.
\end{align} </math>
On obtient de la même manière <math>w_1 = \frac{b-a}2</math>. On a ainsi retrouvé les coefficients de quadrature de la méthode des trapèzes.
Premières formules de Newton-Cotes
Soit un intervalle Modèle:Math séparé en n intervalles de longueur Modèle:Math. On note Modèle:Math et ξ un élément indéterminé de Modèle:Math. Les formules relatives aux premiers degrés sont résumées dans le tableau suivant :
| Degré | Nom commun | Formule | Terme d'erreur |
| 1 | Méthode des trapèzes | <math> \frac{b-a}{2} (f_0 + f_1) </math> | <math>-\frac{(b-a)^3}{12}\,f^{(2)}(\xi)</math> |
| 2 | Méthode de Simpson 1/3 | <math> \frac{b-a}{6} (f_0 + 4 f_1 + f_2) </math> | <math>-\frac{(b-a)^5}{2880}\,f^{(4)}(\xi)</math> |
| 3 | Méthode de Simpson 3/8 | <math> \frac{b-a}{8} (f_0 + 3 f_1 + 3 f_2 + f_3) </math> | <math>-\frac{(b-a)^5}{6480}\,f^{(4)}(\xi)</math> |
| 4 | Méthode de Boole-Villarceau | <math> \frac{b-a}{90} (7 f_0 + 32 f_1 + 12 f_2 + 32 f_3 + 7 f_4) </math> | <math>-\frac{(b-a)^7}{1935360}\,f^{(6)}(\xi)</math> |
| 6 | Méthode de Weddle-Hardy | <math> \frac{b-a}{840} (41 f_0 + 216 f_1 + 27 f_2 + 272 f_3 + 27 f_4 + 216 f_5 + 41 f_6)</math> | <math>-\frac{(b-a)^9}{1567641600}\,f^{(8)}(\xi)</math> |
Les formules relatives aux degrés supérieurs sont donnés dans le tableau suivant :
| Degré | Nombre de points | Formule | Terme d'erreur |
| 7 | Méthode à 8 points<ref name="Wolfram">Weisstein, Eric W. "Newton-Cotes Formulas." From MathWorld--A Wolfram Web Resource</ref> | <math> \frac{b-a}{17 280} (751 (f_0+f_7) + 3577 (f_1+f_6) + 1323 (f_2+f_5) + 2989 (f_3 + f_4)) </math> | <math>-\frac{8183}{518 400}\frac{(b-a)^9}{7^9}\,f^{(8)}(\xi)</math> |
| 8 | Méthode à 9 points<ref name="Wolfram" /> | <math> \frac{b-a}{28 350} (989 (f_0+f_8) + 5888 (f_1+f_7) -928 (f_2+f_6)+ 10496 (f_3+f_5) - 4540 f_4 </math> | <math>-\frac{2 368}{467 775} \frac{(b-a)^{11}}{8^{11}}\,f^{(10)}(\xi)</math> |
| 9 | Méthode à 10 points<ref name="Wolfram" /> | <math> \frac{b-a}{89 600} ( 2857 (f_0+f_9) + 15 741 (f_1+f_8)+1080 (f_2+f_7)+ 19 344( f_3+f_6)+5778(f_4+f_5))</math> | <math>-\frac{519}{394 240} \frac{(b-a)^{11}}{9^{10}}\,f^{(10)}(\xi)</math> |
| 10 | Méthode à 11 points<ref name="Wolfram" /> | <math> \frac{b-a}{598 752} (16067 (f_0+f_{10}) + 106300 (f_1+f_9) - 48525 (f_2+f_8) + 272400 (f_3+f_7) -260550 (f_4+f_6) + 427368 f_5)</math> | <math>-\frac{1 346 350}{326 918 592} \frac{(b-a)^{13}}{{10}^{13}}\,f^{(12)}(\xi)</math> |
Ordre de la méthode
L'ordre d'une formule de quadrature est définie comme le plus grand entier Modèle:Math pour lequel la valeur calculée par la formule vaut exactement l'intégrale recherchée pour tout polynôme de degré inférieur ou égal à Modèle:Math.
L'ordre de la formule de Newton-Cotes de degré Modèle:Math est supérieur ou égal à Modèle:Math, car on a alors Modèle:Math pour tout Modèle:Math polynôme de degré inférieur ou égal à Modèle:Math.
On peut en fait montrer le résultat suivant<ref name=Demailly>Modèle:Ouvrage.</ref>: Modèle:Énoncé
L'ordre donne une indication de l'efficacité d'une formule de quadrature. Les formules de Newton-Cotes sont donc généralement utilisées pour des degrés pairs.
Convergence
Bien qu'une formule de Newton-Cotes puisse être établie pour n'importe quel degré, l'utilisation de degrés supérieurs peut causer des erreurs d'arrondi<ref name=Demailly/>, et la convergence n’est pas assurée lorsque le degré augmente à cause du phénomène de Runge. Pour cette raison, il est généralement préférable de se restreindre aux premiers degrés, et d'utiliser des formules composites pour améliorer la précision de la formule de quadrature. Toutefois, la méthode de Newton-Cotes d'ordre 8 est employée dans le livre Computer Methods for Mathematical Computations, de Forsythe, Malcolm et Moler, qui a joui d'un succès certain dans les années 70 et 80. Elle y apparaît sous la forme d'une méthode adaptative : QUANC8<ref>Code source de QUANC8</ref>.
Références
Liens externes
- Formules de Newton-Cotes sur Math-Linux.com
- Modèle:Mathworld
