Riemann Sums. Outline. James K. Peterson. September 15, Riemann Sums. Riemann Sums In MatLab

Transcription

1 Riemann Sums James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University September 15, 2013 Outline Riemann Sums Riemann Sums In MatLab

2 Abstract This lecture introduces Riemann sums and Riemann integrals. We talked earlier about the special limiting process that leads to a Riemann integral and we have gone over derivatives and antiderivatives of a wide class of functions. But we need desperately a connection between these ideas. We are going to find the easy way to calculate a Riemann integral. No limits at all! Here is the rule b f (t) dt = F (b) F (a) a where F is any antiderivative of f! All we have to do is guess one and we have had a lot of practice doing just that in the last chapter. Of course, this is not magic and it follows from careful reasoning. So get ready for a nice long read. Remember the punch line: this new tool helps us and understanding how a tool is made helps us master it. We must learn learn to make them do our bidding. So here we go with another step on your road to becoming craftsman.

3 Let s go back and be formal for a bit. To study the integration of a function f, there are two intellectually separate ideas: the primitive or antiderivative and the Riemann integral. They seem quite different, don t they? We have 1. the idea of a Primitive or Antiderivative: f. This is any function F which is differentiable and satisfies F (t) = f (t) at all points in the domain of f. Normally, the domain of f is a finite interval of the form [a, b], although it could also be an infinite interval like all of R or [1, ) and so on. Note that an antiderivative does not require any understanding of the process of Riemann integration at all only what differentiation is! 2. The idea of the Riemann integral of a function. Riemann integration is far more complicated to setup than the process of guessing a primitive or antiderivative. To define a Riemann integral properly, first, we start with a bounded function f on a finite interval [a, b]. This kind of function f need not be continuous! Then select a finite number of points from the interval [a, b], {t 0, t 1,,..., t n 1, t n}. We don t know how many points there are, so a different selection from the interval would possibly gives us more or less points. But for convenience, we will just call the last point t n and the first point t 0. These points are not arbitrary t 0 is always a, t n is always b and they are ordered like this: t 0 = a < t 1 < t 2 <... < t n 1 < t n = b The collection of points from the interval [a, b] is called a Partition of [a, b] and is denoted by some letter here we will use the letter P.

4 So if we say P is a partition of [a, b], we know it will have n + 1 points in it, they will be labeled from t 0 to t n and they will be ordered left to right with strict inequalities. But, we will not know what value the positive integer n actually is. The simplest Partition P is the two point partition {a, b}. Note these things also: 1. Each partition of n + 1 points determines n subintervals of [a, b] 2. The lengths of these subintervals always adds up to the length of [a, b] itself, b a. 3. These subintervals can be represented as {[t0, t1], [t1, t2],..., [tn 1, tn]} or more abstractly as [ti, ti+1] where the index i ranges from 0 to n The length of each subinterval is ti+1 ti for the indices i in the range 0 to n The largest subinterval length is called the norm of the partition and we denote it by the symbol P. Now from each subinterval [t i, t i+1] determined by the Partition P, select any point you want and call it s i. This will give us the points s 0 from [t 0, t 1], s 1 from [t 1, t 2] and so on up to the last point, s n 1 from [t n 1, t n]. At each of these points, we can evaluate the function f to get the value f (s j). Call these points an Evaluation Set for the partition P. Let s denote such an evaluation set by the letter E. If the function f was nice enough to be positive always and continuous, then the product f (s i) (t i+1 t i) can be interpreted as the area of a rectangle.

5 Then, if we add up all these rectangle areas we get a sum which is useful enough to be given a special name: the Riemann Sum for the function f associated with the Partition P and our choice of evaluation set E = {s 0,..., s n 1}. This sum is represented by the symbol S(f, P, E) where the things inside the parenthesis are there to remind us that this sum depends on our choice of the function f, the partition P and the evaluation set E. We have S(f, P, E) = n 1 f (s i) (t i+1 t i) i=0 Definition The Riemann sum for the bounded function f, the partition P and the evaluation set E = {s 0,..., s n 1} from P{t 0, t 1,,..., t n 1, t n} is defined by S(f, P, E) = n 1 f (s i) (t i+1 t i) i=0 It is pretty misleading to write the Riemann sum this way as it can make us think that the n is always the same when in fact it can change value each time we select a different P. So many of us write the definition this way instead S(f, P, E) = i P f (si) (ti+1 ti) and we just remember that the choice of P will determine the size of n.

6 Example Let s look at an example of all this. Here we see the graph of a typical function which is always positive on some finite interval [a, b]. (a, f (a)) (b, f (b)) A generic curve f on the interval [a, b] which is always positive. Note the area under this curve is the shaded region. a b Figure: The Area Under The Curve f Example Continued Next, let s set the interval to be [1, 6] and compute the Riemann Sum for a particular choice of Partition P and evaluation set E. The partition (gray) is P = {1.0, 1.5, 2.6, 3.8, 4.3, 5.6, 6.0}. Hence, we have subinterval lengths of t 1 t 0 = 0.5, t 2 t 1 = 1.1, t 3 t 2 = 1.2, t 4 t 3 = 0.5, t 5 t 4 = 1.3 and t 6 t 5 = 0.4, giving P = 1.3. For the evaluation set (red) E = {1.1, 1.8, 3.0, 4.1, 5.3, 5.8} the Riemann sum is S(f, P, E) = 5 f (s i) (t i+1 t i) i=0 = f (1.1) f (1.8) f (3.0) f (4.1) f (5.3) f (5.8) 0.4

7 Example Continued (1, f (1)) (6, f (6)) The partition points are in gray and the evaluation points are in red. 1 6 Figure: A Simple Riemann Sum Example Continued We can also interpret the Riemann sum as an approximation to the area under the curve. The partition (gray) is P = {1.0, 1.5, 2.6, 3.8, 4.3, 5.6, 6.0}. For the evaluation set (red) E = {1.1, 1.8, 3.0, 4.1, 5.3, 5.8} (1, f (1)) (6, f (6)) 1 6 Figure: The Riemann Sum As An Approximate Area

8 Example Example We let f (t) = t 2 + 6t 8 on the interval [2, 4] with P = {2, 2.5, 3.0, 3.7, 4.0} and E = {2.2, 2.8, 3.3, 3.8}. Solution The partition determines subinterval lengths of t 1 t 0 = 0.5, t 2 t 1 = 0.5, t 3 t 2 = 0.7, and t 4 t 3 = 0.3, giving P = 0.7. For E and P, we have the Riemann sum S(f, P, E) = 3 i=0 f (si) (ti+1 ti. Thus S(f, P, E) = f (2.2) f (2.8) f (3.3) f (3.8) 0.3 = Example Example Let f (t) = 3t 2 on the interval [ 1, 2] with P = { 1, 0.3, 0.6, 1.2, 2.0} and E = { 0.7, 0.2, 0.9, 1.6}. Find the Riemann sum. Solution The partition determines subinterval lengths of t 1 t 0 = 0.7, t 2 t 1 = 0.9, t 3 t 2 = 0.6, and t 4 t 3 = 0.8, giving P = 0.9. For the evaluation set E the Riemann sum is S(f, P, E) = f ( 0.7) f (0.2) f (0.9) f (1.6) 0.8 = 8.739

9 Homework 24 For the given function f, partition P and evaluation set E, do the following. 1. Find S(f, P, E) for the partition P and evaluation set E. 2. Find P. 3. Sketch a graph of this Riemann sum as an approximation to the area under the curve f. Do a nice graph with appropriate use of color. Homework 24 Continued 24.1 Let f (t) = t on the interval [1, 3] with P = {1, 1.5, 2.0, 2.5, 3.0} and E = {1.2, 1.8, 2.3, 2.8} Let f (t) = t on the interval [1, 3] with P = {1, 1.6, 2.3, 2.8, 3.0} and E = {1.2, 1.9, 2.5, 2.85} Let f (t) = 3t 2 + 2t on the interval [1, 2] with P = {1, 1.2, 1.5, 1.8, 2.0} and E = {1.1, 1.3, 1.7, 1.9} 24.4 Let f (t) = 3t 2 + t on the interval [1, 4] with P = {1, 1.2, 1.5, 2.8, 4.0} and E = {1.1, 1.3, 2.3, 3.2}

10 Let s start with Riemann sums in Matlab. As you have seen, doing these by hand is tedious. Let s look at how we might do them using MatLab. Here is a typical MatLab session to do this. Let s calculate the Riemann sum for the function f (x) = x 2 on the interval [1, 3] using the partition P = {1, 1.5, 2.1, 2.8, 3.0} and evaluation set E = {1.2, 1.7, 2.5, 2.9}. First, set up our function. MatLab allows us to define a function inside the MatLab environment as follows >> f x ) ( x. ˆ 2 ) ; This defines the function f (x) = x 2. If we had wanted to define g(x) = 2x 2 + 3, we would have used >> g x ) ( 2 x. ˆ 2 + 3) ; To square X, we would write X.^2 to square each component creating a new vector with each entry squared. >> X. ˆ 2 ans = The way we set up the function f (x.^2); makes use of this. The variable X may or may not be a vector. So we write x.^2 so that if x is a vector, multiplication is done component wise and if not, it is just the squaring of a number. So for our function, to find f for all the values in X, we just type >> f (X) ans = 1 4 9

11 Now let s setup the partition with the command >> P = [ 1 ; 1. 5 ; 2. 1 ; 2. 8 ; 3. 0 ] P = The command diff in MatLab is applied to a vector to create the differences we have called the x i s. >> dx = d i f f (P) dx = Next, we set up the evaluation set E. >> E = [ 1. 2 ; 1. 7 ; 2. 5 ; 2. 9 ] E = Find f (E); a new vector with the values f (s i) s. Use f (E). dx to create the new vector with components f (s i) x i. >> g = f ( E ). dx g =

12 Finally, we add all these components together to get the Riemann sum. In MatLab, we add up the entries of a vector g with the command sum(g). >> RS = sum ( g ) RS = Without the comments, the MatLab session is not too long. >> f x ) ( x. ˆ 2 ) ; >> P = [ 1 ; 1. 5 ; 2. 1 ; 2. 8 ; 3. 0 ] ; >> dx = d i f f (P) ; >> E = [ 1. 2 ; 1. 7 ; 2. 5 ; 2. 9 ] ; >> g = f ( E ). dx ; >> RS = sum ( g ) ; Homework 25 For the given function f, partition P and evaluation set E, do the following: use Matlab to find S(f, P, E) for the partition P and evaluation set E. Create a new word document for this homework. Do the document in single space. Do matlab fragments in bold font. The document starts with your name, MTHSC 106-Bio, Section Number, Date and Homework number.

13 Homework 25 Continued Then answer the problems like this: 1. State Problem 1. insert into your doc the matlab commands you use to solve the problem. Do this in bold. before each line of matlab add explanatory comments so we can check to see you know what you re doing. 2. State Problem 2. same stuff Homework 25 Example Jim Peterson MTHSC 106-B, Section Number today s date and HW Number, Problem 1: Let f (t) = sin(5t) on the interval [1, 3] with P = {1, 1.5, 2.0, 2.5, 3.0} and E = {1.2, 1.8, 2.3, 2.8}. Find S(f, P, E), % add e x p l a n a t i o n h e r e >> f x ) s i n ( 5 x ) ; % add e x p l a n a t i o n h e r e >> P = [ 1 ; 1. 5 ; 2. 0 ; 2. 5 ; 3. 0 ] ; % add e x p l a n a t i o n h e r e >> E = [ 1. 2 ; 1. 8 ; 2. 3 ; 2. 8 ] ; % add e x p l a n a t i o n h e r e >> dx = d i f f (P) ; % add e x p l a n a t i o n h e r e >> g = f ( E ). dx ; % add e x p l a n a t i o n h e r e >> RS = sum ( g ) RS =

14 Homework 25 Continued 25.1 Let f (t) = t on the interval [1, 3] with P = {1, 1.5, 2.0, 2.5, 3.0} and E = {1.2, 1.8, 2.3, 2.8} Let f (t) = t on the interval [1, 3] with P = {1, 1.6, 2.3, 2.8, 3.0} and E = {1.2, 1.9, 2.5, 2.85} Let f (t) = 3t 2 + 2t on the interval [1, 2] with P = {1, 1.2, 1.5, 1.8, 2.0} and E = {1.1, 1.3, 1.7, 1.9} 25.4 Let f (t) = 3t 2 + t on the interval [1, 4] with P = {1, 1.2, 1.5, 2.8, 4.0} and E = {1.1, 1.3, 2.3, 3.2}