1 / 22

Lecture 8

Lecture 8. Integration. Indefinite Integrals Indefinite Integrals of a function are functions that differ from each other by a constant. Definite Integrals Definite Integrals are numbers. Why Numerical Integration?.

missy
Télécharger la présentation

Lecture 8

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Lecture 8

  2. Integration Indefinite Integrals Indefinite Integrals of a function are functionsthat differ from each other by a constant. Definite Integrals Definite Integrals are numbers.

  3. Why Numerical Integration? • Very often, the function f(x) to differentiate or the integrand to integrate is too complex to derive exact analytical solutions. • In most cases in engineering, the function f(x) is only available in a tabulated form with values known only at discrete points. Numerical Solution

  4. Numerical Integration The general form of numerical integration of a function f (x) over some interval [a, b] is a weighted sum of the function values at a finite number (n) of sample points (nodes), referred to as ‘quadrature’:

  5. One interpretation of the definite integral is Integral = area under the curve f(x) a b

  6. Newton-Cotes Integration • Common numerical integration scheme • Based on the strategy of replacing a complicated function or tabulated data with some approximating function that is easy to integrate Pn(x) is an nth order polynomial

  7. Trapezoidal Rule • Corresponds to the case where the polynomial is a first order F(b) h F(a)

  8. From the trapezoidal rule we can obtain for the total area of (n-1) intervals where there are nequally spaced base points.

  9. Error Estimate in the trapezoidal rule Error It can be obtained by integrating the interpolation error we defined in previous chapter for Lagrange polynomial as

  10. Example

  11. 0 1

  12. 1 2 1 1/2 0 1 1 1 2 2 2 1/2 3/4 0 1/4 1

  13. Remark 1: in this example instead of re-computation of some function values when h is changed to h/2 we observe that

  14. Simpson’s Rules Simpson’s 1/3 rule can be obtained by passing a parabolic interpolant through three adjacent nodes. The area under the parabola is

  15. To obtain the total area of (n-1) even intervals we apply the following general Simpson’s 1/3 rule 1 Note: f(x1), f(xn) 4 f(x2), f(x4), f(x6),.. 2 f(x3), f(x5), f(x7),..

  16. Remark 2:Simpson’s 1/3 rule requires the number of intervals to be even. If this condition is not satisfied, we can integrate over the first (or last) three intervals with Simpson’s 3/8rule which can be obtained by passing a cubic interpolant through four adjacent nodes, and defined by The error in the Simpson’s rule is

  17. Simpson’s 1/3 rule Simpson’s 3/8 rule • Because the number of panels is odd, we compute the integral over the first three intervals by Simpson’s 3/8 rule, and use the 1/3 rule for the last two intervals:

  18. Summary • Newton Cotes formulae for Numerical integration. • Trapezoidal Rule • Simpson’s Rules. • Romberg Integration • Double Integrals To be continued in Lecture 9

More Related