|
|
| Line 1: |
Line 1: |
|
| |
|
| [[File:Riemann sum convergence.png|right|thumb|300px|Four of the Riemann summation [[#Methods|methods]] for approximating the area under curves. <span style="color:#0081cd">'''Right''' </span> and <span style="color:#fec200">'''left'''</span> methods make the approximation using the right and left endpoints of each subinterval, respectively. <span style="color:#009246">'''Maximum'''</span> and <span style="color:#bc1e47">'''minimum'''</span> methods make the approximation using the largest and smallest endpoint values of each subinterval, respectively. The values of the sums converge as the subintervals halve from top-left to bottom-right.]]
| |
|
| |
|
| In [[mathematics]], a '''Riemann sum''' is an [[approximation]] of the [[area]] of a region, often the region underneath a curve. It is named after German mathematician [[Bernhard Riemann]].
| | Hi there. The author's name is Eusebio but rather he never really liked that name. In his professional life he is also a people manager. He's always loved living inside Guam and he owns everything that he prerequisites there. To drive is one of those things he loves a large [http://Www.answers.com/topic/percentage percentage] of. He's been working around his website for some people time now. Check understand it out here: http://prometeu.net<br><br>Here is my homepage :: [http://prometeu.net clash of clans hack ipad download] |
| | |
| The sum is calculated by dividing the region up into shapes ([[rectangle]]s or [[trapezoid]]s) that together form a region that is similar to the region being measured, then calculating the area for each of these shapes, and finally adding all of these small areas together. This approach can be used to find a numerical approximation for a [[definite integral]] even if the [[fundamental theorem of calculus]] does not make it easy to find a [[Closed-form_expression|closed-form solution]]. | |
| | |
| Because the region filled by the small shapes is usually not exactly the same shape as the region being measured, the Riemann sum will differ from the area being measured. This error can be reduced by dividing up the region more finely, using smaller and smaller shapes. As the shapes get smaller and smaller, the sum approaches the [[Riemann integral]].
| |
| | |
| ==Definition==
| |
| Let ''f'' : ''D'' → '''R''' be a function defined on a subset, ''D'', of the real line, '''R'''. Let ''I'' = [''a'', ''b''] be a closed interval contained in ''D'', and let
| |
| :<math>P= \left \{[x_0,x_1),[x_1,x_2),\dots,[x_{n-1},x_{n}] \right \},</math>
| |
| be a partition of ''I'', where
| |
| :<math>a=x_0<x_1<x_2<\cdots<x_n=b.</math>
| |
| The '''Riemann sum''' of ''f'' over ''I'' with partition ''P'' is defined as
| |
| :<math>S = \sum_{i=1}^{n} f(x_i^*)(x_{i}-x_{i-1}), \quad x_{i-1}\le x_i^* \le x_i.</math>
| |
| The choice of <math>x_i^*</math> in the interval <math>[x_{i-1},x_i]</math> is arbitrary.
| |
| | |
| '''Example:''' Specific choices of <math>x_i^*</math> give us different types of Riemann sums:
| |
| * If <math>x_i^*=x_{i-1}</math> for all ''i'', then ''S'' is called a '''left Riemann sum'''.
| |
| * If <math>x_i^*=x_i</math> for all ''i'', then ''S'' is called a '''right Riemann sum'''.
| |
| * If <math>x_i^*=\tfrac{1}{2}(x_i+x_{i-1})</math> for all ''i'', then ''S'' is called a '''middle Riemann sum'''.
| |
| * The average of the left and right Riemann sum is the '''trapezoidal sum'''.
| |
| * If it is given that
| |
| ::<math>S = \sum_{i=1}^{n} v_i(x_{i}-x_{i-1}),</math>
| |
| :where <math>v_i</math> is the [[supremum]] of ''f'' over <math>[x_{i-1},x_i]</math>, then ''S'' is defined to be an '''upper Riemann sum'''.
| |
| * Similarly, if <math>v_i</math> is the [[infimum]] of ''f'' over <math>[x_{i-1},x_i]</math>, then ''S'' is a '''lower Riemann sum'''.
| |
| | |
| {{multiple image
| |
| | align = right
| |
| | direction = vertical
| |
| | header = Riemann sum methods of ''x''<sup>3</sup> over [0,2] using 4 subdivisions
| |
| | width = 200
| |
| | image1 = LeftRiemann2.svg
| |
| | caption1 = Left
| |
| | image2 = RightRiemann2.svg
| |
| | caption2 = Right
| |
| | image3 = MidRiemann2.svg
| |
| | caption3 = Middle
| |
| | image4 = TrapRiemann2.svg
| |
| | caption4 = Trapezoidal
| |
| | image5 = Simpsons_method_illustration.png
| |
| | caption5 = Simpson's Method
| |
| }}
| |
| Any Riemann sum on a given partition (that is, for any choice of <math>x_i^*</math> between <math>x_{i-1}</math> and <math>x_i</math>) is contained between the lower and the upper Riemann sums. A function is defined to be ''[[Riemann integral|Riemann integrable]]'' if the lower and upper Riemann sums get ever closer as the partition gets finer and finer. This fact can also be used for [[numerical integration]].
| |
| | |
| ==Methods==
| |
| The four methods of Riemann summation are usually best approached with partitions of equal size. The interval [''a'', ''b''] is therefore divided into ''n'' subintervals, each of length
| |
| | |
| :<math>\Delta x = \frac{b-a}{n}.</math>
| |
| | |
| The points in the partition will then be
| |
| | |
| : <math>a, a + \Delta x, a + 2 \Delta x, \ldots, a + (n-2) \Delta x, a + (n-1) \Delta x, b.</math>
| |
| | |
| ===Left Riemann Sum===
| |
| For the left Riemann sum, approximating the function by its value at the left-end point gives multiple rectangles with base Δ''x'' and height ''f''(''a'' + ''i''Δ''x''). Doing this for ''i'' = 0, 1, ..., ''n'' − 1, and adding up the resulting areas gives
| |
| | |
| :<math>\Delta x \left[f(a) + f(a + \Delta x) + f(a + 2 \Delta x)+\cdots+f(b - \Delta x)\right].</math>
| |
| | |
| The left Riemann sum amounts to an overestimation if ''f'' is [[monotonically decreasing]] on this interval, and an underestimation if it is [[monotonically increasing]].
| |
| | |
| ===Right Riemann Sum===
| |
| ''f'' is here approximated by the value at the right endpoint. This gives multiple rectangles with base Δ''x'' and height ''f''(''a'' + ''i''Δ''x''). Doing this for ''i'' = 1, ..., ''n'', and adding up the resulting areas produces
| |
| :<math>\Delta x \left[ f( a + \Delta x ) + f(a + 2 \Delta x)+\cdots+f(b) \right].</math> | |
| | |
| The right Riemann sum amounts to an underestimation if ''f'' is [[monotonically decreasing]], and an overestimation if it is [[monotonically increasing]].
| |
| | |
| ===Middle sum===
| |
| Approximating ''f'' at the midpoint of intervals gives ''f''(''a'' + ''Q''/2) for the first interval, for the next one ''f''(''a'' + 3''Q''/2), and so on until ''f''(''b'' − ''Q''/2). Summing up the areas gives
| |
| | |
| :<math>Q\left[f(a + \tfrac{Q}{2}) + f(a + \tfrac{3Q}{2})+\cdots+f(b-\tfrac{Q}{2})\right].</math>
| |
| | |
| The error of this formula will be
| |
| | |
| :<math>\left \vert \int_{a}^{b} f(x) \, dx - A_\mathrm{mid} \right \vert \le \frac{M_2(b-a)^3}{24n^2}, </math>
| |
| | |
| where <math>M_2</math> is the maximum value of the [[absolute value]] of <math>f^{\prime\prime}(x)</math> on the interval.
| |
| | |
| ===Trapezoidal Rule===
| |
| In this case, the values of the function ''f'' on an interval are approximated by the average of the values at the left and right endpoints. In the same manner as above, a simple calculation using the area formula
| |
| :<math>A=\tfrac{1}{2}h(b_1+b_2)</math>
| |
| for a [[trapezoid|trapezium]] with parallel sides ''b''<sub>1</sub>, ''b''<sub>2</sub> and height ''h'' produces
| |
| | |
| :<math>\tfrac{1}{2}Q\left[f(a) + 2f(a+Q) + 2f(a+2Q) + 2f(a+3Q)+\cdots+f(b)\right].</math>
| |
| | |
| The error of this formula will be
| |
| | |
| :<math>\left \vert \int_{a}^{b} f(x) \, dx - A_\mathrm{trap} \right \vert \le \frac{M_2(b-a)^3}{12n^2},</math> | |
| | |
| where <math>M_2</math> is the maximum value of the absolute value of <math>f^{\prime\prime}(x).</math>
| |
| | |
| The approximation obtained with the trapezoid rule for a function is the same as the average of the left hand and right hand sums of that function.
| |
| | |
| ==Example==
| |
| | |
| [[Image:Area-under-curve-for-x-squared.svg|thumb|200px|right|A visual representation of the area under the curve ''y'' = ''x''<sup>2</sup> for the interval from 0 to 2]]
| |
| [[File:Riemann sum error.svg|thumb|The value of the Riemann sum under the curve ''y'' = ''x''<sup>2</sup> from 0 to 2. As the number of rectangles increases, it approaches the exact area of 8/3.]]
| |
| Taking an example, the area under the curve of ''y'' = ''x''<sup>2</sup> between 0 and 2 can be procedurally computed using Riemann's method.
| |
| | |
| The interval [0, 2] is firstly divided into ''n'' subintervals, each of which is given a width of <math>\tfrac{2}{n}</math>; these are the widths of the Riemann rectangles (hereafter "boxes"). Because the right Riemann sum is to be used, the sequence of ''x'' coordinates for the boxes will be <math>x_1, x_2, \ldots, x_n</math>. Therefore, the sequence of the heights of the boxes will be <math>x_1^2, x_2^2, \ldots, x_n^2</math>. It is an important fact that <math>x_i = \tfrac{2i}{n}</math>, and <math>x_n = 2</math>.
| |
| | |
| The area of each box will be <math>\tfrac{2}{n} \times x_i^2</math> and therefore the ''n''th right Riemann sum will be:
| |
| :<math>\begin{align}
| |
| S &= \frac{2}{n} \times \left(\frac{2}{n}\right)^2 + \cdots + \frac{2}{n} \times \left(\frac{2i}{n}\right)^2 + \cdots + \frac{2}{n} \times \left(\frac{2n}{n}\right)^2 \\
| |
| &= \frac{8}{n^3} \left(1 + \cdots + i^2 + \cdots + n^2\right)\\
| |
| &= \frac{8}{n^3} \left(\frac{n(n+1)(2n+1)}{6}\right)\\
| |
| &= \frac{8}{n^3} \left(\frac{2n^3+3n^2+n}{6}\right)\\
| |
| &= \frac{8}{3} + \frac{4}{n} + \frac{4}{3n^2}
| |
| \end{align}</math>
| |
| | |
| If the limit is viewed as ''n'' → ∞, it can be concluded that the approximation approaches the actual value of the area under the curve as the number of boxes increases. Hence:
| |
| :<math>\lim_{n \to \infty} S = \lim_{n \to \infty}\left(\frac{8}{3} + \frac{4}{n} + \frac{4}{3n^2}\right) = \frac{8}{3}</math>
| |
| | |
| This method agrees with the definite integral as calculated in more mechanical ways:
| |
| | |
| :<math>\int_0^2 x^2\, dx = \frac{8}{3}</math>
| |
| | |
| == Animations ==
| |
| <gallery>
| |
| Image:Riemann sum (leftbox).gif|Left sum
| |
| Image:Riemann sum (rightbox).gif|Right sum
| |
| Image:Riemann sum (middlebox).gif|Middle sum
| |
| Image:Riemann sum (y=x^2).gif|With <math>y=x^2</math>
| |
| </gallery>
| |
| | |
| ==See also==
| |
| * [[Riemann integral]]
| |
| * [[Riemann–Stieltjes integral]]
| |
| * [[Lebesgue integral]]
| |
| * [[Simpson's rule]]
| |
| * [[Euler method]] and [[midpoint method]], related methods for solving differential equations
| |
| | |
| ==References==
| |
| *{{citation| last1=Thomas|first1=George B. Jr.|last2=Finney|first2=Ross L.| title=Calculus and Analytic Geometry (9th ed.)| publisher=Addison Wesley| year=1996| isbn=0-201-53174-7}}
| |
| | |
| ==External links==
| |
| *[http://www.vias.org/simulations/simusoft_riemannsum.html A simulation showing the convergence of Riemann sums]
| |
| | |
| [[Category:Integral calculus]]
| |