释义 |
Newton-Cotes FormulasThe Newton-Cotes formulas are an extremely useful and straightforward family of Numerical Integration techniques.
To integrate a function over some interval , divide it into equal parts such that and . Then find Polynomials which approximate the tabulated function, and integrate them toapproximate the Area under the curve. To find the fitting Polynomials, use LagrangeInterpolating Polynomials. The resulting formulas are called Newton-Cotes formulas, orQuadrature Formulas.
Newton-Cotes formulas may be ``closed'' if the interval is included in the fit, ``open'' if the points are used, or a variation of these two. If the formula uses points (closed or open), theCoefficients of terms sum to .
If the function is given explicitly instead of simply being tabulated at the values , the best numericalmethod of integration is called Gaussian Quadrature. By picking the intervals at which to sample the function,this procedure produces more accurate approximations (but is significantly more complicated to implement).
The 2-point closed Newton-Cotes formula is called the Trapezoidal Rule because it approximates the area under acurve by a Trapezoid with horizontal base and sloped top (connecting the endpoints and ). If the firstpoint is , then the other endpoint will be located at
and the Lagrange Interpolating Polynomial through the points and is
Integrating over the interval (i.e., finding the area of the trapezoid) then gives
This is the trapezoidal rule, with the final term giving the amount of error (which, since , is no worsethan the maximum value of in this range).
The 3-point rule is known as Simpson's Rule. The Abscissas are
and the Lagrange Interpolating Polynomial isIntegrating and simplifying gives
 | (7) |
The 4-point closed rule is Simpson's 3/8 Rule,
 | (8) |
The 5-point closed rule is Bode's Rule,
 | (9) |
(Abramowitz and Stegun 1972, p. 886). Higher order rules include the 6-point
 | (10) |
7-point | |  | (11) | 8-point | |  | | | (12) | 9-point | |  | (13) | 10-pointand 11-pointrules.
Closed ``extended'' rules use multiple copies of lower order closed rules to build up higher order rules. Byappropriately tailoring this process, rules with particularly nice properties can be constructed. For tabulated points, usingthe Trapezoidal Rule times and adding the results gives Using a series of refinements on the extended Trapezoidal Rule gives the method known as Romberg Integration.A 3-point extended rule for Odd isApplying Simpson's 3/8 Rule, then Simpson's Rule (3-point) twice, and adding givesTaking the next Simpson's 3/8 step then gives
 | (19) |
Combining with the previous result giveswhere terms up to have now been completely determined. Continuing gives | |  | (21) | Now average with the 3-point result
 | (22) |
to obtain | |  | (23) | Note that all the middle terms now have unity Coefficients.Similarly, combining a 4-point with the (2+4)-point rule gives
 | (24) |
Other Newton-Cotes rules occasionally encountered include Durand's Rule
 | (25) |
(Beyer 1987), Hardy's Rule | |  | (26) | and Weddle's Rule
 | (27) |
(Beyer 1987).
The open Newton-Cotes rules use points outside the integration interval, yielding the 1-point
 | (28) |
2-point3-point
 | (30) |
4-point
 | (31) |
5-point
 | (32) |
6-point | |  | (33) | and 7-point | |  | (34) | rules.
A 2-point open extended formula is  | |  | | | (35) |
Single interval extrapolative rules estimate the integral in an interval based on the points around it. Anexample of such a rule is
 | (36) |
 | (37) |
 | (38) |
 | (39) |
See also Bode's Rule, Difference Equation, Durand's Rule, Finite Difference, GaussianQuadrature, Hardy's Rule, Lagrange Interpolating Polynomial, Numerical Integration, Simpson'sRule, Simpson's 3/8 Rule, Trapezoidal Rule, Weddle's Rule References
Abramowitz, M. and Stegun, C. A. (Eds.). ``Integration.'' §25.4 in Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, 9th printing. New York: Dover, pp. 885-887, 1972.Beyer, W. H. (Ed.) CRC Standard Mathematical Tables, 28th ed. Boca Raton, FL: CRC Press, p. 127, 1987. Hildebrand, F. B. Introduction to Numerical Analysis. New York: McGraw-Hill, pp. 160-161, 1956. Press, W. H.; Flannery, B. P.; Teukolsky, S. A.; and Vetterling, W. T. ``Classical Formulas for Equally Spaced Abscissas.'' §4.1 in Numerical Recipes in FORTRAN: The Art of Scientific Computing, 2nd ed. Cambridge, England: Cambridge University Press, pp. 124-130, 1992. |