HIGHER-ORDER LINEAR ORDINARY DIFFERENTIAL EQUATIONS IV: Laplace Transform Method. David Levermore Department of Mathematics University of Maryland

Size: px
Start display at page:

Download "HIGHER-ORDER LINEAR ORDINARY DIFFERENTIAL EQUATIONS IV: Laplace Transform Method. David Levermore Department of Mathematics University of Maryland"

Transcription

1 HIGHER-ORDER LINEAR ORDINARY DIFFERENTIAL EQUATIONS IV: Laplace Transform Method David Levermore Department of Mathematics University of Maryland 6 April Because the presentation of this material in lecture will differ from that in the book, I felt that notes that closely follow the lecture presentation might be appreciated. 8. Laplace Transform Method 8.. Definition of the Transform 8.. Properties of the Transform Existence and Differentiability Transform of Derivatives Application to Initial-Value Problems 8.6. Piecewise Defined Forcing Inverse Transform Computing Green Functions Convolutions (not covered 8.. Green Functions and Natural Fundamental Sets (not covered 3

2 8. Laplace Transform Method The Laplace transform will allow us to transform an initial-value problem for a linear ordinary differential equation with constant coefficients into a linear algebaric equation that can be easily solved. The solution of an initial-value problem can then be obtained from the solution of the algebaric equation by taking its so-called inverse Laplace transform. 8.. Definition of the Transform. The Laplace transform of a function f(t defined over t is another function L[f](s that is formally defined by L[f](s = e st f(t dt. (8. You should recall from calculus that the above definite integral is improper because its upper endpoint is. The proper definition of the Laplace transform is therefore L[f](s = lim T e st f(t dt, (8. provided that the definite integrals over [, T] appearing in the above limit are proper. The Laplace transform L[f](s is defined only at those s for which the limit in (8. exists. Example. Use definition (8. to compute L[e at ](s for any real a. Solution. From (8. you see that for any s a one has L[e at ](s = lim e st e at dt = lim e (a st dt T T = lim T while for s = a one has e (a st a s T t= [ = lim T s a e(a st s a L[e at ](s = lim e (s at dt = lim T T Therefore L[e at ](s is only defined for s > a with ] for s > a, = s a for s < a, dt = lim T T =. L[e at ](s = s a for s > a.

3 3 Example. Use definition (8. to compute L[t e at ](s for any real a. Solution. From (8. you see that for any s a one has L[t e at ](s = lim T = lim T [ = lim T ( t t e st e at dt = lim T T e (a st a s (a s (s a t e (a st dt t= ( T s a + (s a ] e (a st = (s a for s > a, for s < a, while for s = a one has L[t e at ](s = lim t e (s at dt = lim t dt = lim T T T T =. Therefore L[t e at ](s is only defined for s > a with L[t e at ](s = (s a for s > a. Example. Use definition (8. to compute L[e ibt ](s for any real b. Solution. For b you see from (8. that for any real s we have L[e ibt ](s = lim e st e ibt dt = lim T T [ ] = lim T s ib e (s ibt s ib e (s ibt dt = lim T for s >, = s ib undefined for s. ( e (s ibt T s ib The case b = is identical to our first example with a =. In every case L[e ibt ](s is only defined for s > with L[e ibt ](s = for s >. s ib 8.. Properties of the Transform. If we always had to return to the definition of the Laplace transform everytime we wanted to apply it, it would not be easy to use. Rather, we will use the definition to compute the Laplace transform for a few basic functions and to establish some general properties that will allow us to build formulas for more complicated functions. The most important such property is the fact that the Laplace transform L is a linear operator. t=

4 4 Theorem. If L[f](s and L[g](s exist for some s then so does L[f + g](s and L[cf](s for every constant c with L[f + g](s = L[f](s + L[g](s, L[cf](s = cl[f](s. (8.3 Proof. This follows directly from definition (8. and the facts that definite integrals and limits depend linearly on their arguments. Specifically, one sees that L[f + g](s = lim T = lim T L[cf](s = lim T e st( f(t + g(t dt e st f(t dt + lim T e st cf(t dt = c lim T e st g(t dt = L[f](s + L[g](s, e st f(t dt = cl[f](s. Example. Compute L[cos(bt](s and L[sin(bt](s for any real b. Solution. This can be done by using the Euler identity e ibt = cos(bt + i sin(bt and the linearity of L. Then L[cos(bt](s + il[sin(bt](s = L[e ibt ](s = s ib = s + ib s + b for s >. By equating the real and imaginary parts above, we see that s L[cos(bt](s = s + b for s >, b L[sin(bt](s = s + b for s >. Another general property of the Laplace transform is that it turns multiplication by an exponential in t into a translation of s. Theorem. If L[f](s exists for every s > α and a is any real number then L[e at f(t](s exists for every s > α + a with L[e at f(t](s = L[f](s a. Proof. This follows directly from definition (8.. Specifically, one sees that L[e at f](s = lim e st e at f(t dt = lim e (s at f(t dt = L[f](s a. T T

5 5 Examples. From our previous examples and the above theorem we see that L[e (a+ibt ](s = s a ib for s > a, L[e at s a cos(bt](s = (s a + b for s > a, L[e at sin(bt](s = b (s a + b for s > a. The Laplace transform also turns a translation of t into multiplication by an exponential in s. Notice that L[f](s only depends on the values of f(t over [.. Therefore before we translate f we multiply it by the unit step function u(t defined by { for t, u(t = for t <. (8.4 Because the functions uf and f agree over [,, it is clear that L[uf](s = L[f](s. We now consider the Laplace transform of the translation u(t cf(t c for every c >. Theorem. If L[f](s exists for every s > α and c > then L[u(t cf(t c](s exists for every s > α with L[u(t cf(t c](s = e cs L[f](s for s > α. Proof. For every T > c one has Therefore e st u(t cf(t c dt = c e st f(t c dt = e cs e s(t c f(t c dt c = e cs e st f(t dt. L[u(t cf(t c](s = lim e st u(t cf(t c dt T = e cs lim T c c e st f(t dt = e cs L[f](s for s > α.

6 Existence and Differentiablity of the Transform. In each of the above examples the definite integrals over [, T] that appear in the limit (8. were proper. Indeed, we were able to evaluate the definite integrals analytically and determine the limit (8. for every real s. More generally, from calculus we know that a definite integral over [, T] is proper whenever its integrand is: bounded over [, T], continuous at all but a finite number of points in [, T]. Such an integrand is said to be piecewise continuous over [, T]. Because for e st is continuous (and therefore bounded function of t over every [, T] for each real s, the definite integrals over [, T] that appear in the limit (8. will be proper whenever f(t is piecewise continuous over every [, T]. Example. The function { for t < π, f(t = cos(t for t π, is piecewise continuous over every [, T] because it is clearly bounded over [, and its only discontinuity is at the point t = π. Example. The so-called sawtooth function f(t = t k for k t < k + where k =,,, 3,, is piecewise continuous over every [, T] because it is clearly bounded over [, and has discontinuities at the points t =,, 3,, only a finite number of which lie in each [, T]. If we assume that f(t is piecewise continuous over every [, T], we still have to give a condition under which the limit (8. will exist for certain s. Such a condition is provided by the following definition. Definition. A function f(t defined over [, is said to be of exponential order α as t provided that for every σ > α there exist K σ and T σ such that f(t K σ e σt for every t T σ. (8.5 Roughly speaking, a function is of exponential order α as t if its absolute value does not grow faster than e σt as t for every σ > α. Example. The function e at is of exponential order a as t because (8.5 holds with K σ = and T σ = for every σ > a. Example. The function cos(bt is of exponential order as t because (8.5 holds with K σ = and T σ = for every σ >.

7 Example. For every p > the function t p is of exponential order as t. Indeed, for every σ > the function e σt t p takes on its maximum over [, at t = p/σ, whereby ( p p e σt t p for every t. eσ It follows that (8.5 holds with K σ = ( p eσ p and T σ = for every σ >. One can show that if functions f and g are of exponential orders α and β respectively as t then the function f + g is of exponential order max{α, β} as t, while the function fg is of exponential order α + β as t. Example. For every real a the function e at + e at is of exponential order a as t. This is because the functions e at and e at are exponential orders a and a respectively as t, and because a = max{a, a}. Example. For every p > and every real a and b the function t p e at cos(bt is of exponential order a as t. This is because the functions t p, e at, and cos(bt are of exponential orders, a, and respectively as t. The fact you should know about the existence of the Laplace transform for certain s is the following. Theorem. Let f(t be piecewise continuous over every [, T], of exponential order α as t. Then for every positive integer k the function t k f(t has these same properties. The function F(s = L[f](s is defined for every s > α. Moreover, F(s is infinitely differentiable over s > α with L[t k f(t](s = ( k dk F(s for s > α. (8.6 dsk Proof. Formula (8.6 can be derived formally by differentiating the integrands: F(s = d ds F(s = d ds F(s =. d k F(s = ( k dsk e st f(t dt, t e st f(t dt, t e st f(t dt, t k e st f(t dt. 7

8 8 A correct proof would require some justification of taking the derivatives inside the above imporper integrals. We will not go into those details here. We will however give an easier proof of the fact that F(s is defined for s > α. The proof uses the direct comparison test for the convergence of improper integrals. That test implies that if g(t and G(t are piecewise continuous over every [, T] such that g(t G(t for every t then G(t dt converges = g(t dt converges. Let s > α and apply this test to g(t = e st f(t. Pick σ so that α < σ < s. Because f(t is of exponential order α as t and σ > α there exist K σ and T σ such that (8.5 holds. Because g(t = e st f(t is bounded over [, T σ ] there exists B σ such that g(t B σ over [, T σ ]. It thereby follows that { g(t = e st Bσ for t < T σ f(t G(t K σ e (σ st for t T σ. Because s > σ for this G(t you can show that G(t dt = lim T G(t dt converges. It follows that the limit in (8. converges, whereby F(s = L[f](s is defined at s. Example. Because for every real a and b we have L[e (a+ibt ](s = s a ib for s > a, it follows from the above theorem that for every nonnegative integer k L[t k e (a+ibt ](s = ( k dk ds k s a ib = k! (s a ib k+ for s > a. This formula implies that for every real a and b and every nonnegative integer k L[t k ](s = k! for s >, s k+ L[t k e at k! ](s = (s a k+ for s > a. ( L[t k e at k! cos(bt](s = Re for s > a, (s a ib ( k+ L[t k e at k! sin(bt](s = Im for s > a. (s a ib k+

9 8.4. Transform of Dervatives. The previous result shows that the Laplace transform turns a multiplication by t into a derivative with respect to s. The next result shows the Laplace transform turns a derivative with respect to t into a multiplication by s. Theorem. Let f(t be continuous over [, such that f(t is of exponential order α as t, f (t is piecewise continuous over every [, T]. Then L[f ](s is defined for every s > α with L[f ](s = s L[f](s f(. 9 Proof. Let s > α. By definition (8., an integration by parts, the fact that f(t is of exponential order α as t, and the fact that L[f](s exists, one sees that [ T L[f ](s = lim e st f (t dt = lim e st f(t T + s T T t= = lim T e st f(t f( + s lim T e st f(t dt e st f(t dt = f( + s L[f](s. If f(t is sufficiently differentiable then the previous result can be applied repeatedly. For example, if f(t is twice differentiable then L[f ](s = s L[f ](s f ( = s ( s L[f](s f( f ( If f(t is thrice differentiable then = s L[f](s s f( f (. L[f ](s = s L[f ](s f ( = s ( s L[f](s s f( f ( f ( = s 3 L[f](s s f( s f ( f (. Proceeding in this way you can use induction to prove the following. Theorem. Let f(t be n-times differentiable over [, such that f(t, f (t,, f (n (t are of exponential order α as t, f (n (t is piecewise continuous over every interval [, T]. Then L[f (n ](s is defined for every s > α with L[f (n ](s = s n L[f](s s n f( s n f ( s f (n ( f (n (. (8.8 ]

10 This means that if you know that a function y(t is n-times differentiable and that it and its first n derivatives are of exponential order as t then if Y (s = L[y](s one has L[y ](s = s Y (s y(, L[y ](s = s Y (s s y( y (, L[y ](s = s 3 Y (s s y( s y ( y (,. L[y (n ](s = s n Y (s s n y( s n y ( s y (n ( y (n (. ( Application to Initial-Value Problems. Because the Laplace transform turns derivatives with respect to t into multiplications by s, it transforms initial-value problems into algebraic problems. Suppose that y(t is the solution of the initial-value problem y (n + a y (n + + a n y + a n y = f(t, y( = y, y ( = y, y (n ( = y n. It can be shown that if f(t is piecewise continuous over [, and is of exponential order as t then y(t is n-times differentiable and that it and its first n derivatives are of exponential order as t. You can thereby use (8.9 to find Y (s = L[y](s in terms of the inital data y, y,, y n, and the Laplace transform of the forcing, F(s = L[f](s. Indeed, the fact L is a linear operator implies that the Laplace transform of the initial-value problem is If we use (8.9 then we find L[y (n ] + a L[y (n ] + + a n L[y ] + a n L[y] = L[f]. where p(s is the characteristic polynomial p(sy (s = q(s + F(s, p(s = s n + a s n + + a n s + a n. and q(s is the polynomial given in terms of the initial data by q(s = ( s n + a s n + + a n s + a n y + ( s n + a s n a n 3 s + a n y. + ( s + a s + a yn 3 + ( s + a yn + y n. You therefore find that Y (s = q(s + F(s p(s. (8.

11 Later we will see how to determine y(t from Y (s, but first we will illustrate how to compute Y (s. The hardest part of doing this is often computing F(s = L[f](s in (8.. Usually f(t can be expressed as a combination of the basic forms whose Laplace transform we have already computed. Some of these basic forms are L[t n ](s = n! s n+ L[cos(bt](s = L[sin(bt](s = s s + b b s + b L[e at j(t](s = J(s a for s >, for s >, for s >, where J(s = L[j(t](s, L[t n j(t](s = ( n J (n (s L[u(t cj(t c](s = e cs J(s L[e at t n ](s = L[e at cos(bt](s = L[e at sin(bt](s = n! (s a n+ s a (s a + b b (s a + b where J(s = L[j(t](s, where J(s = L[j(t](s and u(t is the unit step function, for s > a, for s > a, for s > a. (8. One can use these to build up a much longer table of basic forms such as the one given in the book. However, the above table contains all the forms you really need. In fact, you can argue the first three entries are redundent because they follow by setting a = in the last three. Alternatively, you can argue that the last three entries are redundent because they follow immediately from the first three and the fourth. On exams you will be given a table that includes at least the last six entries above, so there is no need to memorize this table. However, you should learn how to use it efficiently. Example. Find the Laplace transform Y (s of the solution y(t of the initial-value problem y y = e 5t, y( = 3. Solution. By setting a = 5 and n = in the seventh entry of table (8. you see that L[e 5t ](s = /(s 5. The Laplace transform of the initial-value problem therefore is L[y ](s L[y](s = L[e 5t ](s = s 5,

12 where we see from (8.9 that It follows that L[y](s = Y (s, L[y ](s = s Y (s y( = s Y (s 3. (s Y (s 3 = s 5, = Y (s = (s (s s. Example. Find the Laplace transform Y (s of the solution y(t of the initial-value problem y y 8y =, y( = 3, y ( = 7. Solution. Here there is no forcing. The Laplace transform of the initial-value problem is where we see from (8.9 that It follows that L[y ](s L[y ](s 8L[y](s =, L[y](s = Y (s, L[y ](s = s Y (s y( = s Y (s 3, L[y ](s = s Y (s s y( y ( = s Y (s 3s 7. (s s 8Y (s 3s =, = Y (s = 3s + s s 8. Example. Find the Laplace transform Y (s of the solution y(t of the initial-value problem y + 4y = sin(3t, y( = y ( =. Solution. By setting b = 3 in the third entry of table (8. you see that L[sin(3t](s = 3/(s + 9. The Laplace transform of the initial-value problem therefore is where we see from (8.9 that L[y ](s + 4L[y](s = L[sin(3t](s = 3 s + 3 = 3 s + 9, L[y](s = Y (s, L[y ](s = s Y (s sy( y ( = s Y (s. It follows that (s + 4Y (s = 3 s + 9, = Y (s = 3 (s + 4(s + 9.

13 8.6. Piecewise-Defined Forcing. The Laplace transform method can be to solve initialvalue problems of the form y (n + a y (n + + a n y + a n y = f(t, y( = y, y ( = y, y (n ( = y n, where the forcing f(t is piecewise-defined over [, by a list of cases given in the form f (t for t < c, f (t for c t < c, f(t =.. f m (t for c m t < c m, f m (t for c m t <, where = c < c < c < < c m <. We assume that for each k =,,, m the function f k is continuous and bounded over [c k, c k+, while the function f m is continuous over [c m, and is of exponential order as t. Here we show how to compute the Laplace transform F(s = L[f](s for such a function. There are three steps. The first step is to express f(t in terms of translations of the unit step u(t. How this is done becomes clear once you see that for every c < d one has { for c t < d, u(t c u(t d = otherwise. In other words, the function u(t c u(t d is a switch that turns on at t = c and turns off at t = d. So for any given function g(t one has ( u(t c u(t d g(t = { g(t for c t < d, This observation allows you to express f(t as otherwise. f(t = ( u(t u(t c f (t + ( u(t c u(t c f (t + + ( u(t c m u(t c m f m (t + u(t c m f m (t. By grouping terms above that involve the same u(t c k, you bring f(t into the form f(t = f (t + u(t c h (t + + u(t c m h m (t, (8. where h k (t = f k (t f k (t for k =,,, m. This is the form you want. You can get to it either by carrying out the grouping indicated above or by memorizing (8.. It helps to recall that each term u(t c k h k (t appearing in (8. simply changes the forcing from f k (t to f k (t at time t = c k because h k (t = f k (t f k (t. 3

14 4 The idea of the second step is to bring (8. into a form that allows you to use the sixth entry in table (8.. That entry states that L[u(t cj(t c](s = e cs L[j](s. You must therefore recast (8. into the form f(t = f (t + u(t c j (t c + + u(t c m j m (t c m. Each function j k (t is obtained from the h k (t appearing in (8. by j k (t = h k (t + c k for k =,,, m. Indeed, this formula implies j k (t c k = h k (t, which is what appears in (8.. Once you have found all the j k (t then the final step is to compute L[f ](s and each L[j k ](s, and use the fact that to compute L[f](s as L[u(t c k j k (t c k ](s = e c ks L[j k ](s for k =,,, m, L[f](s = L[f ](s + e c s L[j ](s + + e c ms L[j m ](s. Often you will have to use identities to express f (t and each j k (t in forms that allows you to compute their Laplace transforms from table (8.. Example. Find the Laplace transform Y (s of the solution y(t of the initial-value problem y + 4y = f(t, y( = 7, y ( = 5, where t for t <, f(t = t for t < 4, 4 for 4 t. Solution. The Laplace transform of the initial-value problem is where F(s = L[f](s and L[y](s = Y (s, L[y ](s + 4L[y](s = F(s, L[y ](s = sy (s y( = sy (s 7, L[y ](s = s Y (s sy( y ( = s Y (s 7s 5. The Laplace transform of the initial-value problem thereby becomes (s + 4Y (s 7s 5 = F(s, = Y (s = ( 7s F(s. s + 4

15 All that remains to be done is to compute F(s. The first step is to express f(t in terms of unit step functions as f(t = ( u(t u(t t + ( u(t u(t 4 t + u(t 4 4 The second step is to write = t + u(t (t t + u(t 4 (4 t. f(t = t + u(t j (t + u(t 4 j (t 4, 5 where j (t = (t + (t + = t + 4 t 4t 4 = t t, j (t = 4 (t + 4 = t 4. Here we obtained j (t by replacing t with t+ in the factor (t t and j (t by replacing t with t + 4 in the factor (4 t. Finally, the above form for f(t allows you to use the sixth entry of table (8. to compute F(s = L[f](s as It follows that F(s = L[t ](s + L[u(t j (t ](s + L[u(t 4 j (t 4](s = L[t ](s e s L[t + t](s e 4s L[t + 4](s = ( s 3 e s s 3 + ( s e 4s s + 4 s = ( e s s 3 ( e s + e 4s s 4 e 4s s. Y (s = 7s + 5 s ( e s s 3 (s + 4 ( e s + e 4s s (s + 4 e 4s 4 s (s Inverse Transform. The process of determining y(t from Y (s is called taking the inverse Laplace transform. It is important to know that this process has a unique result. Indeed, we will use the following theorem. Theorem. Let f(t and g(t be two functions over [, and α a real number such that f(t and g(t are of exponential order α as t, f(t and g(t are piecewise continuous over every [, T], L[f](s = L[g](s for every s > α. Then f(t = g(t for every t in [,.

16 6 The proof of this result requires tools from complex variables that are beyond the scope of this course. Fortunately, you do not need to know how to prove this result to use it! Its usefulness stems from the fact that solutions y(t to the initial-value problems we are considering lie within the class of functions considered above namely, they are functions that are of exponential order as t and that are piecewise continuous over every [, T]. In fact, they are continuous and piecewise differentiable over every [, T]. This means that if we succeed in finding a function y(t within this class such that L[y](s = Y (s then it will be the unique solution of the initial-value problem that we seek. Because the above result states there is a unique f(t that is of exponential order as t and is piecewise continuous over every [, T] such that L[f](s = F(s, we introduce the notation f(t = L [F](t. The operator L denotes the inverse Laplace transform. Because it undoes the Laplace transform L, it inherits many properties from L. For example, it is linear. You can also easily read-off from the first and last three entries in table (8. of basic forms that [ ] n! L s n+ (t = t n, [ ] L s s + b (t = cos(bt, [ ] L b s + b (t = sin(bt, It is also clear from the sixth entry of table (8. that [ ] L n! (s a n+ (t = e at t n, [ ] L s a (s a + b (t = e at cos(bt, [ ] L b (s a + b (t = e at sin(bt. (8.3 L [e cs J(s](t = u(t c j(t c, where j(t = L [J](t. (8.4 For us, the process of computing y(t = L [Y ](t for a given Y (s will be one of expressing Y (s as a sum of terms that will allow us to read off y(t from the basic forms above. To illustrate this process, we will compute y(t = L [Y ](t for the Y (s found in the examples given in the previous section, thereby completing our solution of the initial-value problems. Example. Find y(t = L [Y ](t for Y (s = (s (s s. Solution. By the partial fraction identity (s (s 5 = 3 s s,

17 7 you can express Y (s as Y (s = 3 s s. The top right entry of table (8.3 with a = 5 and a = then yields [ ] [ ] y(t = L [Y (s](t = 3 L (t + 8 s 5 3 L (t = s 3 e5t et. Example. Find y(t = L [Y ](t for Solution. By the partial fraction identity you can express Y (s as Y (s = 3s + s s s + s s 8 = 3s + (s 4(s + = 6 s s +, Y (s = 3 6 s s +. The top right entry of table (8.3 with a = 4 and a = then yields [ ] [ ] y(t = L [Y (s](t = 3 6 L (t + 5 s 4 6 L (t = 3 s + 6 e4t e t. Example. Find y(t = L [Y ](t for Y (s = Solution. By the partial fraction identity 3 (s + 4(s + 9. you can express Y (s as 3 3 (z + 4(z + 9 = 5 (z (z + 9, Y (s = s s + 9 = 3 s s + 3. The bottom left entry of table (8.3 with b = and b = 3 then yields [ [ ] y(t = L [Y (s](t = 3 ](t L 5 3 s + L (t = 3 s + 3 sin(t 5 sin(3t.

18 8 Example. Find y(t = L [Y ](t for Y (s = 7s + 5 s ( e s s 3 (s + 4 ( e s + e 4s s (s + 4 e 4s 4 s (s + 4. Solution. You first derive the partial fraction identities 7s + 5 s + 4 = 7s s s + 4, s 3 (s + 4 = s 8 3 s + 8 s s + 4, s (s + 4 = s s + 4, 4 s (s + 4 = s s s + 4. The top left of these is straightforward. The top right identity only involves s, so is simply the identity z(z + 4 = z z + 4, evaluated at z = s. The bottom right identity is simply s times the top right one. Finally, the bottom left identity is obtained by first dividing the top right one by s and then employing the bottom right one divided by 8 to the last term. These partial fraction identites allow you to express Y (s as Y (s = 7s s s ( e s ( s 8 3 s + 8 s s + 4 ( e s + e 4s( s = 7 s ( e s + e 4s( s s + 4 s + + ( e s s 4 s + The formulas in the first column of table (8.3 show that ( e 4s s s s + 4 ( 4 s 3 8 s + s 8 s + ( e 4s s s. s + L [7 s s ] (t = 7 cos(t + 5 s + cos(t, [ L 4 s 3 8 s + ] s (t = 8 s + 4 t cos(t, [ L s ] (t = 4 s + t 4 sin(t, [ ] L s s (t = cos(t. s +

19 9 By combining these facts with formula (8.4, it follows that y(t = 7 cos(t + 5 cos(t + ( 4 t cos(t u(t ( 4 (t cos((t u(t ( (t 4 sin((t u(t 4 ( (t 4 4 sin((t 4 u(t 4 ( cos((t Computing Green Functions. The Laplace transform can be used to efficiently compute Green functions for differential operators with constant coefficients. Recall that given the n th -order differential operator L with constant coefficients given by L = D n + a D n + + a n D + a n, the Green function g(t associated with L is the solution of the initial-value problem g (n + a g (n + + a n g + a n g =, g( =, g ( =, g (n ( =, g (n ( =. The Laplace transform of this initial-value problem is L [ g (n] (s + a L [ g (n ] (s + + a n L[g ](s + L[g](s =, where if G(s = L[g](s then L[g ](s = s G(s g( = s G(s, L [ g ] (s = s G(s s g( g ( = s G(s,. L [ g (n ] (s = s n G(s s n g( s n 3 g ( g (n ( = s n G(s, L [ g (n] (s = s n G(s s n g( s n g ( g (n ( = s n G(s. We thereby see that G(s satisfies p(sg(s =, = G(s = p(s, where p(s is the characteristic polynomial of L, which is given by p(s = s n + a s n + + a n s + a n. In other words, the Laplace transform of the Green function of L is the reciprocal of the characteristic polynomial of L.

20 The problem of computing a Green function is thereby reduced to the problem of finding an inverse Laplace transform. This can often be done quickly. Example. Find the Green function g(t for the operator L = D + 6D + 3. Solution. Because p(s = s + 6s + 3 = (s + 3 +, the bottom right entry of table (8.3 with a = 3 and b = shows that the Green function is given by [ ] g(t = L = [ ] p(s L = (s e 3t sin(t. Example. Find the Green function g(t for the operator L = D + D 5. Solution. Because p(s = s +s 5 = (s 3(s+5, we use the partial fraction identity (s 3(s + 5 = 8 s 3 8 s + 5. The top right entry of table (8.3 with a = 3 and with a = 5 shows that the Green function is given by [ ] g(t = L = [ ] p(s 8 L [ ] s 3 8 L = e3t e 5t. s Example. Find the Green function g(t for the operator L = D 4 + 3D Solution. Because p(s = s 4 + 3s + 36 = (s + 4(s + 9 only depends on s, we can use the partial fraction identity (z + 4(z + 9 = 5 z z + 9 at z = s. The bottom left entry of table (8.3 with b = and with b = 3 shows that the Green function is given by [ ] g(t = L = [ ] p(s L [ ] 3 s + 5 L = sin(t sin(3t. s Convolutions. Let f(t and g(t be any two functions defined over the interval [,. Their convolution is a third function (f g(t that is defined by the formula (f g(t = t f(t τg(τ dτ, (8.5 whenever the above integral makes sense for every t. In particular, the convolution of f and g will be defined whenever both f and g are piecewise continuous over every [, T].

21 The convolution can be thought of a some kind of product between two functions. It is easily checked that this so-called convolution product satisfies some of the properties of ordinary multipulcation. For example, for any functions f, g, amd h that are piecewise continuous over every [, T] we have g f = f g h (f + g = h f + h g h (g f = (h g f commutative law, distributive law, associative law. The commutative law is proved by introducing τ = t τ as a new variable of integration, whereby one sees that (g f(t = t g(t τf(τ dτ = t g(τ f(t τ dτ = (f g(t. Verification of the distributive and associative laws is left to you. The convolution differs from ordinary multipulcation in some respects too. For example, it is not generaly true that f = f or that f f. Indeed, one sees that ( (t = t dτ = t, and that t t t (sin sin(t = sin(t τ sin(τ dτ = sin(t cos(τ sin(τ dτ + cos(t sin(τ dτ = sin(t3 + t cos(t sin(t cos(t for every t >. The main result of this section is that the Laplace transform of a convolution of two functions is the ordinary product of their Laplace transforms. In other words, the Laplace transform maps convolutions to multiplication. Convolution Theorem. Let f(t and g(t be piecewise continuous over every [, T] of exponential order α as t. Then L[f g](s is defined for every s > α with L[f g](s = F(sG(s, where F(s = L[f](s and G(s = L[g](s. (8.6 Proof. For every T > definition (8.5 of convolution implies that e st (f g(t dt = t e st f(t τg(τ dτ dt = t e st f(t τg(τ dτ dt.

22 We now exchange the order of the definite integrals over τ and t on the right-hand side. As you recall from Calculus III, this should be done carefully because the upper endpoint of the inner integral depends on the variable of integration t of the outer integral. When viewed in the (τ, t-plane, the domain over which the double integral is being taken is the triangle given by τ t T. In general, when the order of definite integrals is exchanged over this domain we have t dτ dt = τ dt dτ, where denotes any appropriate integrand. We thereby obtain e st (f g(t dt = τ e st f(t τg(τ dt dτ. We now factor e st as e st = e s(t τ e sτ, and group the factor e s(t τ with f(t τ and the factor e sτ with g(τ, whereby e st (f g(t dt = = τ e s(t τ f(t τ e sτ g(τ dt dτ e sτ g(τ τ e s(t τ f(t τ dt dτ. We then make the change of variable t = t τ in the inner definite integral to obtain e st (f g(t dt = e sτ g(τ τ e st f(t dt dτ. Upon formally letting T above, definition (8. of the Laplace transform shows that the inner integral converges to F(s, which is independent of τ. The double integral thereby converges to G(sF(s, yielding (8.6. Remark. Because the upper endpoint of the inner integral depends on the variable of integration τ of the outer integral, properly passing to the limit above requires greater care than we took here. The techniques one needs are taught in Advanced Calculus courses. The argument given above suits our purposes because it illuminates why (8.6 holds. The convolution theorem can be used to help evaluate inverse Laplace transforms. For example, suppose that you know for a given F(s and G(s that f(t = L [F](t and g(t = L [G](t. Then (8.6 implies that L [ F(sG(s ] (t = (f g(t. (8.7 You can use this fact to express inverse Laplace transforms as convolutions. You may still have to evaluate the convolution integral, but some of you might find that easier than using partial fraction identities to express F(sG(s in basic forms.

23 3 Example. Compute y(t = L [Y ](t for Y (s = s (s + 4. Solution. Because you know from table (8.3 that [ ] [ ] L s = t, L s + = sin(t, it follows from (8.7 and an integration by parts that [ ] [ ] y(t = L s (s = L t + 4 s s + = (t τ sin(τ dτ = (τ t cos(τ t t cos(τ dτ = t sin(t. 4 This is the same result we got on page 8 using a partial fraction identity. 8.. Green Functions and Natural Fundamental Sets. The convolution theorem gives us other ways to understand Green functions. We have used the Green function to construct a particular solution of the nonhomogeneous equation Ly = p(d = f(t by the formula y P (t = t g(t τf(τ dτ. Notice that the right-hand side above is exactly (g f(t. Taking the Laplace transform of this formula, the Convolution Theorem then yields L[y P ](s = L[g f](s = G(sF(s = F(s p(s, where F(s = L[f](s. But this agrees with formula (8.. Indeed, because y P (t given by the above formula satisfies the initial conditions y P ( =, y P ( =, y(n P ( =, y (n P ( =, it follows that the polynomial q(s appearing in (8. vanishes. Formula (8. can generally be recast as [ ] q(s y(t = y H (t + y P (t, where y H (t = L (t, p(s [ ] F(s y P (t = L (t. p(s

24 4 This is the decomposition of y(t into the solution y H (t of the associated homogeneous equation whose initial data agree with y(t and the particular solution y P (t whose initial data are zero. Because G(s = /p(s, q(s = ( s n + a s n + + a n 3 s + a n s + a n y + ( s n + a s n a n 3 s + a n y. + ( s + a s + a yn 3 + ( s + a yn + y n, and because D k g(t = L [s k G(s](t = L [s k /p(s](t for every k =,,, n, the function y H (t = L [q(s/p(s](t can be expressed in terms of g(t as ( y H (t = y D n + a D n + + a n 3 D + a n D + a n g(t ( + y D n + a D n a n 3 D + a n g(t. + y n 3 ( D + a D + a g(t + y n ( D + a g(t + yn g(t. The natural fundamental set of solutions associated with the homogeneous equation Ly = is therefore given in terms of the Green function g by N (t = (D n + a D n + + a n 3 D + a n D + a n g(t, N (t = (D n + a D n a n 3 D + a n g(t,. N n 3 (t = (D + a D + a g(t, N n (t = ( D + a g(t, (8.8 N n (t = g(t. The solution of the general initial-value problem Ly = f(t, y( = y, y ( = y, y (n ( = y n, y (n ( = y n, can then be expressed as y(t = y N (t + y N (t + + y n N n (t + y n N n (t + (N n f(t, where N (t, N (t,, N n (t is the natural fundamental set of solutions to the associated homogeneous equation, which is given in terms of the Green function g(t by (8.8. Recall that W[N, N,, N n ]( =, which implies W[N, N,, N n ](t = e a t by the Abel formula.

HIGHER-ORDER LINEAR ORDINARY DIFFERENTIAL EQUATIONS IV: Laplace Transform Method David Levermore Department of Mathematics University of Maryland

HIGHER-ORDER LINEAR ORDINARY DIFFERENTIAL EQUATIONS IV: Laplace Transform Method David Levermore Department of Mathematics University of Maryland HIGHER-ORDER LINEAR ORDINARY DIFFERENTIAL EQUATIONS IV: Laplace Transform Method David Levermore Department of Mathematics University of Maryland 9 December Because the presentation of this material in

More information

(f g)(t) = Example 4.5.1: Find f g the convolution of the functions f(t) = e t and g(t) = sin(t). Solution: The definition of convolution is,

(f g)(t) = Example 4.5.1: Find f g the convolution of the functions f(t) = e t and g(t) = sin(t). Solution: The definition of convolution is, .5. Convolutions and Solutions Solutions of initial value problems for linear nonhomogeneous differential equations can be decomposed in a nice way. The part of the solution coming from the initial data

More information

The Laplace transform

The Laplace transform The Laplace transform Samy Tindel Purdue University Differential equations - MA 266 Taken from Elementary differential equations by Boyce and DiPrima Samy T. Laplace transform Differential equations 1

More information

Solutions of the Sample Problems for the Third In-Class Exam Math 246, Fall 2017, Professor David Levermore

Solutions of the Sample Problems for the Third In-Class Exam Math 246, Fall 2017, Professor David Levermore Solutions of the Sample Problems for the Third In-Class Exam Math 6 Fall 07 Professor David Levermore Compute the Laplace transform of ft t e t ut from its definition Solution The definition of the Laplace

More information

HIGHER-ORDER LINEAR ORDINARY DIFFERENTIAL EQUATIONS IV: Variable Coefficient Nonhomogeneous Case

HIGHER-ORDER LINEAR ORDINARY DIFFERENTIAL EQUATIONS IV: Variable Coefficient Nonhomogeneous Case HIGHER-ORDER LINEAR ORDINARY DIFFERENTIAL EQUATIONS IV: Variable Coefficient Nonhomogeneous Case David Levermore Department of Mathematics University of Maryland 15 March 2009 Because the presentation

More information

Section 7.4: Inverse Laplace Transform

Section 7.4: Inverse Laplace Transform Section 74: Inverse Laplace Transform A natural question to ask about any function is whether it has an inverse function We now ask this question about the Laplace transform: given a function F (s), will

More information

HIGHER-ORDER LINEAR ORDINARY DIFFERENTIAL EQUATIONS II: Nonhomogeneous Equations. David Levermore Department of Mathematics University of Maryland

HIGHER-ORDER LINEAR ORDINARY DIFFERENTIAL EQUATIONS II: Nonhomogeneous Equations. David Levermore Department of Mathematics University of Maryland HIGHER-ORDER LINEAR ORDINARY DIFFERENTIAL EQUATIONS II: Nonhomogeneous Equations David Levermore Department of Mathematics University of Maryland 14 March 2012 Because the presentation of this material

More information

The Laplace Transform

The Laplace Transform C H A P T E R 6 The Laplace Transform Many practical engineering problems involve mechanical or electrical systems acted on by discontinuous or impulsive forcing terms. For such problems the methods described

More information

Third In-Class Exam Solutions Math 246, Professor David Levermore Thursday, 3 December 2009 (1) [6] Given that 2 is an eigenvalue of the matrix

Third In-Class Exam Solutions Math 246, Professor David Levermore Thursday, 3 December 2009 (1) [6] Given that 2 is an eigenvalue of the matrix Third In-Class Exam Solutions Math 26, Professor David Levermore Thursday, December 2009 ) [6] Given that 2 is an eigenvalue of the matrix A 2, 0 find all the eigenvectors of A associated with 2. Solution.

More information

Computing inverse Laplace Transforms.

Computing inverse Laplace Transforms. Review Exam 3. Sections 4.-4.5 in Lecture Notes. 60 minutes. 7 problems. 70 grade attempts. (0 attempts per problem. No partial grading. (Exceptions allowed, ask you TA. Integration table included. Complete

More information

HIGHER-ORDER LINEAR ORDINARY DIFFERENTIAL EQUATIONS. David Levermore Department of Mathematics University of Maryland.

HIGHER-ORDER LINEAR ORDINARY DIFFERENTIAL EQUATIONS. David Levermore Department of Mathematics University of Maryland. HIGHER-ORDER LINEAR ORDINARY DIFFERENTIAL EQUATIONS David Levermore Department of Mathematics University of Maryland 28 March 2008 The following is a review of some of the material that we covered on higher-order

More information

Laplace Transform. Chapter 4

Laplace Transform. Chapter 4 Chapter 4 Laplace Transform It s time to stop guessing solutions and find a systematic way of finding solutions to non homogeneous linear ODEs. We define the Laplace transform of a function f in the following

More information

Math 2C03 - Differential Equations. Slides shown in class - Winter Laplace Transforms. March 4, 5, 9, 11, 12, 16,

Math 2C03 - Differential Equations. Slides shown in class - Winter Laplace Transforms. March 4, 5, 9, 11, 12, 16, Math 2C03 - Differential Equations Slides shown in class - Winter 2015 Laplace Transforms March 4, 5, 9, 11, 12, 16, 18... 2015 Laplace Transform used to solve linear ODEs and systems of linear ODEs with

More information

+ + LAPLACE TRANSFORM. Differentiation & Integration of Transforms; Convolution; Partial Fraction Formulas; Systems of DEs; Periodic Functions.

+ + LAPLACE TRANSFORM. Differentiation & Integration of Transforms; Convolution; Partial Fraction Formulas; Systems of DEs; Periodic Functions. COLOR LAYER red LAPLACE TRANSFORM Differentiation & Integration of Transforms; Convolution; Partial Fraction Formulas; Systems of DEs; Periodic Functions. + Differentiation of Transforms. F (s) e st f(t)

More information

Math 353 Lecture Notes Week 6 Laplace Transform: Fundamentals

Math 353 Lecture Notes Week 6 Laplace Transform: Fundamentals Math 353 Lecture Notes Week 6 Laplace Transform: Fundamentals J. Wong (Fall 217) October 7, 217 What did we cover this week? Introduction to the Laplace transform Basic theory Domain and range of L Key

More information

Math Exam 3 Solutions

Math Exam 3 Solutions Math 6 - Exam 3 Solutions Thursday, July 3rd, 0 Recast the following higher-order differential equations into first order systems If the equation is linear, be sure to give the coefficient matrix At and

More information

(an improper integral)

(an improper integral) Chapter 7 Laplace Transforms 7.1 Introduction: A Mixing Problem 7.2 Definition of the Laplace Transform Def 7.1. Let f(t) be a function on [, ). The Laplace transform of f is the function F (s) defined

More information

f(t)e st dt. (4.1) Note that the integral defining the Laplace transform converges for s s 0 provided f(t) Ke s 0t for some constant K.

f(t)e st dt. (4.1) Note that the integral defining the Laplace transform converges for s s 0 provided f(t) Ke s 0t for some constant K. 4 Laplace transforms 4. Definition and basic properties The Laplace transform is a useful tool for solving differential equations, in particular initial value problems. It also provides an example of integral

More information

37. f(t) sin 2t cos 2t 38. f(t) cos 2 t. 39. f(t) sin(4t 5) 40.

37. f(t) sin 2t cos 2t 38. f(t) cos 2 t. 39. f(t) sin(4t 5) 40. 28 CHAPTER 7 THE LAPLACE TRANSFORM EXERCISES 7 In Problems 8 use Definition 7 to find {f(t)} 2 3 4 5 6 7 8 9 f (t),, f (t) 4,, f (t) t,, f (t) 2t,, f (t) sin t,, f (t), cos t, t t t 2 t 2 t t t t t t t

More information

e st f (t) dt = e st tf(t) dt = L {t f(t)} s

e st f (t) dt = e st tf(t) dt = L {t f(t)} s Additional operational properties How to find the Laplace transform of a function f (t) that is multiplied by a monomial t n, the transform of a special type of integral, and the transform of a periodic

More information

APPM 2360: Midterm exam 3 April 19, 2017

APPM 2360: Midterm exam 3 April 19, 2017 APPM 36: Midterm exam 3 April 19, 17 On the front of your Bluebook write: (1) your name, () your instructor s name, (3) your lecture section number and (4) a grading table. Text books, class notes, cell

More information

First-Order Ordinary Differntial Equations: Classification and Linear Equations. David Levermore Department of Mathematics University of Maryland

First-Order Ordinary Differntial Equations: Classification and Linear Equations. David Levermore Department of Mathematics University of Maryland First-Order Ordinary Differntial Equations: Classification and Linear Equations David Levermore Department of Mathematics University of Maryland 1 February 2009 These notes cover some of the material that

More information

UCLA Math 135, Winter 2015 Ordinary Differential Equations

UCLA Math 135, Winter 2015 Ordinary Differential Equations UCLA Math 135, Winter 2015 Ordinary Differential Equations C. David Levermore Department of Mathematics University of Maryland January 5, 2015 Contents 1.1. Normal Forms and Solutions 1.2. Initial-Value

More information

Name: Solutions Exam 3

Name: Solutions Exam 3 Instructions. Answer each of the questions on your own paper. Put your name on each page of your paper. Be sure to show your work so that partial credit can be adequately assessed. Credit will not be given

More information

20. The pole diagram and the Laplace transform

20. The pole diagram and the Laplace transform 95 0. The pole diagram and the Laplace transform When working with the Laplace transform, it is best to think of the variable s in F (s) as ranging over the complex numbers. In the first section below

More information

The Laplace Transform (Sect. 4.1). The Laplace Transform (Sect. 4.1).

The Laplace Transform (Sect. 4.1). The Laplace Transform (Sect. 4.1). The Laplace Transform (Sect. 4.1). s of Laplace Transforms. The Laplace Transform (Sect. 4.1). s of Laplace Transforms. The definition of the Laplace Transform. Definition The function F : D F R is the

More information

ORDINARY DIFFERENTIAL EQUATION: Introduction and First-Order Equations. David Levermore Department of Mathematics University of Maryland

ORDINARY DIFFERENTIAL EQUATION: Introduction and First-Order Equations. David Levermore Department of Mathematics University of Maryland ORDINARY DIFFERENTIAL EQUATION: Introduction and First-Order Equations David Levermore Department of Mathematics University of Maryland 7 September 2009 Because the presentation of this material in class

More information

ORDINARY DIFFERENTIAL EQUATIONS: Introduction and First-Order Equations. David Levermore Department of Mathematics University of Maryland

ORDINARY DIFFERENTIAL EQUATIONS: Introduction and First-Order Equations. David Levermore Department of Mathematics University of Maryland ORDINARY DIFFERENTIAL EQUATIONS: Introduction and First-Order Equations David Levermore Department of Mathematics University of Maryland 1 February 2011 Because the presentation of this material in class

More information

Math 308 Exam II Practice Problems

Math 308 Exam II Practice Problems Math 38 Exam II Practice Problems This review should not be used as your sole source for preparation for the exam. You should also re-work all examples given in lecture and all suggested homework problems..

More information

The Laplace Transform and the IVP (Sect. 6.2).

The Laplace Transform and the IVP (Sect. 6.2). The Laplace Transform and the IVP (Sect..2). Solving differential equations using L ]. Homogeneous IVP. First, second, higher order equations. Non-homogeneous IVP. Recall: Partial fraction decompositions.

More information

Laplace Transforms and use in Automatic Control

Laplace Transforms and use in Automatic Control Laplace Transforms and use in Automatic Control P.S. Gandhi Mechanical Engineering IIT Bombay Acknowledgements: P.Santosh Krishna, SYSCON Recap Fourier series Fourier transform: aperiodic Convolution integral

More information

FIRST-ORDER ORDINARY DIFFERENTIAL EQUATIONS I: Introduction and Analytic Methods David Levermore Department of Mathematics University of Maryland

FIRST-ORDER ORDINARY DIFFERENTIAL EQUATIONS I: Introduction and Analytic Methods David Levermore Department of Mathematics University of Maryland FIRST-ORDER ORDINARY DIFFERENTIAL EQUATIONS I: Introduction and Analytic Methods David Levermore Department of Mathematics University of Maryland 4 September 2012 Because the presentation of this material

More information

27. The pole diagram and the Laplace transform

27. The pole diagram and the Laplace transform 124 27. The pole diagram and the Laplace transform When working with the Laplace transform, it is best to think of the variable s in F (s) as ranging over the complex numbers. In the first section below

More information

Ordinary Differential Equations. Session 7

Ordinary Differential Equations. Session 7 Ordinary Differential Equations. Session 7 Dr. Marco A Roque Sol 11/16/2018 Laplace Transform Among the tools that are very useful for solving linear differential equations are integral transforms. An

More information

Math 3313: Differential Equations Laplace transforms

Math 3313: Differential Equations Laplace transforms Math 3313: Differential Equations Laplace transforms Thomas W. Carr Department of Mathematics Southern Methodist University Dallas, TX Outline Introduction Inverse Laplace transform Solving ODEs with Laplace

More information

Chapter 31. The Laplace Transform The Laplace Transform. The Laplace transform of the function f(t) is defined. e st f(t) dt, L[f(t)] =

Chapter 31. The Laplace Transform The Laplace Transform. The Laplace transform of the function f(t) is defined. e st f(t) dt, L[f(t)] = Chapter 3 The Laplace Transform 3. The Laplace Transform The Laplace transform of the function f(t) is defined L[f(t)] = e st f(t) dt, for all values of s for which the integral exists. The Laplace transform

More information

Ordinary Differential Equations

Ordinary Differential Equations Ordinary Differential Equations for Engineers and Scientists Gregg Waterman Oregon Institute of Technology c 2017 Gregg Waterman This work is licensed under the Creative Commons Attribution 4.0 International

More information

Laplace Transform Theory - 1

Laplace Transform Theory - 1 Laplace Transform Theory - 1 Existence of Laplace Transforms Before continuing our use of Laplace transforms for solving DEs, it is worth digressing through a quick investigation of which functions actually

More information

= AΦ, Φ(0) = I, = k! tk A k = I + ta t2 A t3 A t4 A 4 +, (3) e ta 1

= AΦ, Φ(0) = I, = k! tk A k = I + ta t2 A t3 A t4 A 4 +, (3) e ta 1 Matrix Exponentials Math 26, Fall 200, Professor David Levermore We now consider the homogeneous constant coefficient, vector-valued initial-value problem dx (1) = Ax, x(t I) = x I, where A is a constant

More information

Unit 2: Modeling in the Frequency Domain Part 2: The Laplace Transform. The Laplace Transform. The need for Laplace

Unit 2: Modeling in the Frequency Domain Part 2: The Laplace Transform. The Laplace Transform. The need for Laplace Unit : Modeling in the Frequency Domain Part : Engineering 81: Control Systems I Faculty of Engineering & Applied Science Memorial University of Newfoundland January 1, 010 1 Pair Table Unit, Part : Unit,

More information

Math Shifting theorems

Math Shifting theorems Math 37 - Shifting theorems Erik Kjær Pedersen November 29, 2005 Let us recall the Dirac delta function. It is a function δ(t) which is 0 everywhere but at t = 0 it is so large that b a (δ(t)dt = when

More information

Lecture 23 Using Laplace Transform to Solve Differential Equations

Lecture 23 Using Laplace Transform to Solve Differential Equations Lecture 23 Using Laplace Transform to Solve Differential Equations /02/20 Laplace transform and Its properties. Definition: L{f (s) 6 F(s) 6 0 e st f(t)dt. () Whether L{f (s) or F(s) is used depends on

More information

First-Order Ordinary Differntial Equations II: Autonomous Case. David Levermore Department of Mathematics University of Maryland.

First-Order Ordinary Differntial Equations II: Autonomous Case. David Levermore Department of Mathematics University of Maryland. First-Order Ordinary Differntial Equations II: Autonomous Case David Levermore Department of Mathematics University of Maryland 25 February 2009 These notes cover some of the material that we covered in

More information

Name: Solutions Final Exam

Name: Solutions Final Exam Instructions. Answer each of the questions on your own paper. Put your name on each page of your paper. Be sure to show your work so that partial credit can be adequately assessed. Credit will not be given

More information

The Laplace Transform

The Laplace Transform The Laplace Transform Introduction There are two common approaches to the developing and understanding the Laplace transform It can be viewed as a generalization of the CTFT to include some signals with

More information

ECEN 420 LINEAR CONTROL SYSTEMS. Lecture 2 Laplace Transform I 1/52

ECEN 420 LINEAR CONTROL SYSTEMS. Lecture 2 Laplace Transform I 1/52 1/52 ECEN 420 LINEAR CONTROL SYSTEMS Lecture 2 Laplace Transform I Linear Time Invariant Systems A general LTI system may be described by the linear constant coefficient differential equation: a n d n

More information

HIGHER-ORDER LINEAR ORDINARY DIFFERENTIAL EQUATIONS II: Nonhomogeneous Case and Vibrations

HIGHER-ORDER LINEAR ORDINARY DIFFERENTIAL EQUATIONS II: Nonhomogeneous Case and Vibrations HIGHER-ORDER LINEAR ORDINARY DIFFERENTIAL EQUATIONS II: Nonhomogeneous Case and Vibrations David Levermore Department of Mathematics University of Maryland 22 October 28 Because the presentation of this

More information

MA 201, Mathematics III, July-November 2018, Laplace Transform (Contd.)

MA 201, Mathematics III, July-November 2018, Laplace Transform (Contd.) MA 201, Mathematics III, July-November 2018, Laplace Transform (Contd.) Lecture 19 Lecture 19 MA 201, PDE (2018) 1 / 24 Application of Laplace transform in solving ODEs ODEs with constant coefficients

More information

First and Second Order Differential Equations Lecture 4

First and Second Order Differential Equations Lecture 4 First and Second Order Differential Equations Lecture 4 Dibyajyoti Deb 4.1. Outline of Lecture The Existence and the Uniqueness Theorem Homogeneous Equations with Constant Coefficients 4.2. The Existence

More information

Generalized sources (Sect. 6.5). The Dirac delta generalized function. Definition Consider the sequence of functions for n 1, Remarks:

Generalized sources (Sect. 6.5). The Dirac delta generalized function. Definition Consider the sequence of functions for n 1, Remarks: Generalized sources (Sect. 6.5). The Dirac delta generalized function. Definition Consider the sequence of functions for n, d n, t < δ n (t) = n, t 3 d3 d n, t > n. d t The Dirac delta generalized function

More information

Second In-Class Exam Solutions Math 246, Professor David Levermore Thursday, 31 March 2011

Second In-Class Exam Solutions Math 246, Professor David Levermore Thursday, 31 March 2011 Second In-Class Exam Solutions Math 246, Professor David Levermore Thursday, 31 March 211 (1) [6] Give the interval of definition for the solution of the initial-value problem d 4 y dt 4 + 7 1 t 2 dy dt

More information

e st t u(t 2) dt = lim t dt = T 2 2 e st = T e st lim + e st

e st t u(t 2) dt = lim t dt = T 2 2 e st = T e st lim + e st Math 46, Profeor David Levermore Anwer to Quetion for Dicuion Friday, 7 October 7 Firt Set of Quetion ( Ue the definition of the Laplace tranform to compute Lf]( for the function f(t = u(t t, where u i

More information

The Laplace Transform. Background: Improper Integrals

The Laplace Transform. Background: Improper Integrals The Laplace Transform Background: Improper Integrals Recall: Definite Integral: a, b real numbers, a b; f continuous on [a, b] b a f(x) dx 1 Improper integrals: Type I Infinite interval of integration

More information

ENGIN 211, Engineering Math. Laplace Transforms

ENGIN 211, Engineering Math. Laplace Transforms ENGIN 211, Engineering Math Laplace Transforms 1 Why Laplace Transform? Laplace transform converts a function in the time domain to its frequency domain. It is a powerful, systematic method in solving

More information

9.5 The Transfer Function

9.5 The Transfer Function Lecture Notes on Control Systems/D. Ghose/2012 0 9.5 The Transfer Function Consider the n-th order linear, time-invariant dynamical system. dy a 0 y + a 1 dt + a d 2 y 2 dt + + a d n y 2 n dt b du 0u +

More information

Math 256: Applied Differential Equations: Final Review

Math 256: Applied Differential Equations: Final Review Math 256: Applied Differential Equations: Final Review Chapter 1: Introduction, Sec 1.1, 1.2, 1.3 (a) Differential Equation, Mathematical Model (b) Direction (Slope) Field, Equilibrium Solution (c) Rate

More information

Honors Differential Equations

Honors Differential Equations MIT OpenCourseWare http://ocw.mit.edu 8.034 Honors Differential Equations Spring 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. LECTURE 20. TRANSFORM

More information

Linear algebra and differential equations (Math 54): Lecture 20

Linear algebra and differential equations (Math 54): Lecture 20 Linear algebra and differential equations (Math 54): Lecture 20 Vivek Shende April 7, 2016 Hello and welcome to class! Last time We started discussing differential equations. We found a complete set of

More information

Differential Equations Class Notes

Differential Equations Class Notes Differential Equations Class Notes Dan Wysocki Spring 213 Contents 1 Introduction 2 2 Classification of Differential Equations 6 2.1 Linear vs. Non-Linear.................................. 7 2.2 Seperable

More information

Linear Filters. L[e iωt ] = 2π ĥ(ω)eiωt. Proof: Let L[e iωt ] = ẽ ω (t). Because L is time-invariant, we have that. L[e iω(t a) ] = ẽ ω (t a).

Linear Filters. L[e iωt ] = 2π ĥ(ω)eiωt. Proof: Let L[e iωt ] = ẽ ω (t). Because L is time-invariant, we have that. L[e iω(t a) ] = ẽ ω (t a). Linear Filters 1. Convolutions and filters. A filter is a black box that takes an input signal, processes it, and then returns an output signal that in some way modifies the input. For example, if the

More information

CHAPTER 8 Laplace Transforms

CHAPTER 8 Laplace Transforms CHAPTER Laplace Transforms IN THIS CHAPTER we study the method of Laplace transforms, which illustrates one of the basic problem solving techniques in mathematics: transform a difficult problem into an

More information

4.3 Linear, Homogeneous Equations with Constant Coefficients. Jiwen He

4.3 Linear, Homogeneous Equations with Constant Coefficients. Jiwen He 4.3 Exercises Math 3331 Differential Equations 4.3 Linear, Homogeneous Equations with Constant Coefficients Jiwen He Department of Mathematics, University of Houston jiwenhe@math.uh.edu math.uh.edu/ jiwenhe/math3331

More information

FIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS III: Autonomous Planar Systems David Levermore Department of Mathematics University of Maryland

FIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS III: Autonomous Planar Systems David Levermore Department of Mathematics University of Maryland FIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS III: Autonomous Planar Systems David Levermore Department of Mathematics University of Maryland 4 May 2012 Because the presentation of this material

More information

The Laplace Transform. Background: Improper Integrals

The Laplace Transform. Background: Improper Integrals The Laplace Transform Background: Improper Integrals Recall: Definite Integral: a, b real numbers, a b; f continuous on [a, b] b a f(x) dx 1 Improper integrals: Type I Infinite interval of integration

More information

MATH 251 Examination II April 4, 2016 FORM A. Name: Student Number: Section:

MATH 251 Examination II April 4, 2016 FORM A. Name: Student Number: Section: MATH 251 Examination II April 4, 2016 FORM A Name: Student Number: Section: This exam has 12 questions for a total of 100 points. In order to obtain full credit for partial credit problems, all work must

More information

MA 201, Mathematics III, July-November 2016, Laplace Transform

MA 201, Mathematics III, July-November 2016, Laplace Transform MA 21, Mathematics III, July-November 216, Laplace Transform Lecture 18 Lecture 18 MA 21, PDE (216) 1 / 21 Laplace Transform Let F : [, ) R. If F(t) satisfies the following conditions: F(t) is piecewise

More information

2.161 Signal Processing: Continuous and Discrete Fall 2008

2.161 Signal Processing: Continuous and Discrete Fall 2008 MIT OpenCourseWare http://ocw.mit.edu 2.6 Signal Processing: Continuous and Discrete Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. MASSACHUSETTS

More information

10 Transfer Matrix Models

10 Transfer Matrix Models MIT EECS 6.241 (FALL 26) LECTURE NOTES BY A. MEGRETSKI 1 Transfer Matrix Models So far, transfer matrices were introduced for finite order state space LTI models, in which case they serve as an important

More information

Math 216 Second Midterm 16 November, 2017

Math 216 Second Midterm 16 November, 2017 Math 216 Second Midterm 16 November, 2017 This sample exam is provided to serve as one component of your studying for this exam in this course. Please note that it is not guaranteed to cover the material

More information

Chemical Engineering 436 Laplace Transforms (1)

Chemical Engineering 436 Laplace Transforms (1) Chemical Engineering 436 Laplace Transforms () Why Laplace Transforms?? ) Converts differential equations to algebraic equations- facilitates combination of multiple components in a system to get the total

More information

MATH 251 Examination II April 3, 2017 FORM A. Name: Student Number: Section:

MATH 251 Examination II April 3, 2017 FORM A. Name: Student Number: Section: MATH 251 Examination II April 3, 2017 FORM A Name: Student Number: Section: This exam has 12 questions for a total of 100 points. In order to obtain full credit for partial credit problems, all work must

More information

MTH 464: Computational Linear Algebra

MTH 464: Computational Linear Algebra MTH 464: Computational Linear Algebra Lecture Outlines Exam 4 Material Prof. M. Beauregard Department of Mathematics & Statistics Stephen F. Austin State University April 15, 2018 Linear Algebra (MTH 464)

More information

The Laplace Transform

The Laplace Transform The Laplace Transform Generalizing the Fourier Transform The CTFT expresses a time-domain signal as a linear combination of complex sinusoids of the form e jωt. In the generalization of the CTFT to the

More information

Higher Order Linear Equations

Higher Order Linear Equations C H A P T E R 4 Higher Order Linear Equations 4.1 1. The differential equation is in standard form. Its coefficients, as well as the function g(t) = t, are continuous everywhere. Hence solutions are valid

More information

f (t) K(t, u) d t. f (t) K 1 (t, u) d u. Integral Transform Inverse Fourier Transform

f (t) K(t, u) d t. f (t) K 1 (t, u) d u. Integral Transform Inverse Fourier Transform Integral Transforms Massoud Malek An integral transform maps an equation from its original domain into another domain, where it might be manipulated and solved much more easily than in the original domain.

More information

Math Fall Linear Filters

Math Fall Linear Filters Math 658-6 Fall 212 Linear Filters 1. Convolutions and filters. A filter is a black box that takes an input signal, processes it, and then returns an output signal that in some way modifies the input.

More information

Ch 6.2: Solution of Initial Value Problems

Ch 6.2: Solution of Initial Value Problems Ch 6.2: Solution of Initial Value Problems! The Laplace transform is named for the French mathematician Laplace, who studied this transform in 1782.! The techniques described in this chapter were developed

More information

MAT Linear Algebra Collection of sample exams

MAT Linear Algebra Collection of sample exams MAT 342 - Linear Algebra Collection of sample exams A-x. (0 pts Give the precise definition of the row echelon form. 2. ( 0 pts After performing row reductions on the augmented matrix for a certain system

More information

3.5 Undetermined Coefficients

3.5 Undetermined Coefficients 3.5. UNDETERMINED COEFFICIENTS 153 11. t 2 y + ty + 4y = 0, y(1) = 3, y (1) = 4 12. t 2 y 4ty + 6y = 0, y(0) = 1, y (0) = 1 3.5 Undetermined Coefficients In this section and the next we consider the nonhomogeneous

More information

Lecture 7: Laplace Transform and Its Applications Dr.-Ing. Sudchai Boonto

Lecture 7: Laplace Transform and Its Applications Dr.-Ing. Sudchai Boonto Dr-Ing Sudchai Boonto Department of Control System and Instrumentation Engineering King Mongkut s Unniversity of Technology Thonburi Thailand Outline Motivation The Laplace Transform The Laplace Transform

More information

Math 20D: Form B Final Exam Dec.11 (3:00pm-5:50pm), Show all of your work. No credit will be given for unsupported answers.

Math 20D: Form B Final Exam Dec.11 (3:00pm-5:50pm), Show all of your work. No credit will be given for unsupported answers. Turn off and put away your cell phone. No electronic devices during the exam. No books or other assistance during the exam. Show all of your work. No credit will be given for unsupported answers. Write

More information

8.3 Partial Fraction Decomposition

8.3 Partial Fraction Decomposition 8.3 partial fraction decomposition 575 8.3 Partial Fraction Decomposition Rational functions (polynomials divided by polynomials) and their integrals play important roles in mathematics and applications,

More information

Linear Systems Theory

Linear Systems Theory ME 3253 Linear Systems Theory Review Class Overview and Introduction 1. How to build dynamic system model for physical system? 2. How to analyze the dynamic system? -- Time domain -- Frequency domain (Laplace

More information

The formulas for derivatives are particularly useful because they reduce ODEs to algebraic expressions. Consider the following ODE d 2 dx + p d

The formulas for derivatives are particularly useful because they reduce ODEs to algebraic expressions. Consider the following ODE d 2 dx + p d Solving ODEs using Fourier Transforms The formulas for derivatives are particularly useful because they reduce ODEs to algebraic expressions. Consider the following ODE d 2 dx + p d 2 dx + q f (x) R(x)

More information

Vectors in Function Spaces

Vectors in Function Spaces Jim Lambers MAT 66 Spring Semester 15-16 Lecture 18 Notes These notes correspond to Section 6.3 in the text. Vectors in Function Spaces We begin with some necessary terminology. A vector space V, also

More information

FIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS II: Homogeneous Linear Systems with Constant Coefficients

FIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS II: Homogeneous Linear Systems with Constant Coefficients FIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS II: Homogeneous Linear Systems with Constant Coefficients David Levermore Department of Mathematics University of Maryland 28 January 2012 Because

More information

CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION. Professor Dae Ryook Yang

CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION. Professor Dae Ryook Yang CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dept. of Chemical and Biological Engineering 5-1 Road Map of the Lecture V Laplace Transform and Transfer

More information

Math 307 Lecture 19. Laplace Transforms of Discontinuous Functions. W.R. Casper. Department of Mathematics University of Washington.

Math 307 Lecture 19. Laplace Transforms of Discontinuous Functions. W.R. Casper. Department of Mathematics University of Washington. Math 307 Lecture 19 Laplace Transforms of Discontinuous Functions W.R. Casper Department of Mathematics University of Washington November 26, 2014 Today! Last time: Step Functions This time: Laplace Transforms

More information

Matemáticas II: Segundo del Grado en Ingeniería Aeroespacial

Matemáticas II: Segundo del Grado en Ingeniería Aeroespacial Matemáticas II: Segundo del Grado en Ingeniería Aeroespacial Sergio Blanes, Dolors Roselló Ecuaciones diferenciales y transformadas de Laplace con aplicaciones L.M. Sánchez, M.P. Legua Ref.: 211-798 Capítulo

More information

Infinite series, improper integrals, and Taylor series

Infinite series, improper integrals, and Taylor series Chapter 2 Infinite series, improper integrals, and Taylor series 2. Introduction to series In studying calculus, we have explored a variety of functions. Among the most basic are polynomials, i.e. functions

More information

Chapter 6 The Laplace Transform

Chapter 6 The Laplace Transform Ordinary Differential Equations (Math 2302) 2017-2016 Chapter 6 The Laplace Transform Many practical engineering problems involve mechanical or electrical systems acted on by discontinuous or impulsive

More information

2tdt 1 y = t2 + C y = which implies C = 1 and the solution is y = 1

2tdt 1 y = t2 + C y = which implies C = 1 and the solution is y = 1 Lectures - Week 11 General First Order ODEs & Numerical Methods for IVPs In general, nonlinear problems are much more difficult to solve than linear ones. Unfortunately many phenomena exhibit nonlinear

More information

( ) f (k) = FT (R(x)) = R(k)

( ) f (k) = FT (R(x)) = R(k) Solving ODEs using Fourier Transforms The formulas for derivatives are particularly useful because they reduce ODEs to algebraic expressions. Consider the following ODE d 2 dx + p d 2 dx + q f (x) = R(x)

More information

MA 266 Review Topics - Exam # 2 (updated)

MA 266 Review Topics - Exam # 2 (updated) MA 66 Reiew Topics - Exam # updated Spring First Order Differential Equations Separable, st Order Linear, Homogeneous, Exact Second Order Linear Homogeneous with Equations Constant Coefficients The differential

More information

Outline. Control systems. Lecture-4 Stability. V. Sankaranarayanan. V. Sankaranarayanan Control system

Outline. Control systems. Lecture-4 Stability. V. Sankaranarayanan. V. Sankaranarayanan Control system Outline Control systems Lecture-4 Stability V. Sankaranarayanan Outline Outline 1 Outline Outline 1 2 Concept of Stability Zero State Response: The zero-state response is due to the input only; all the

More information

Laplace Transforms Chapter 3

Laplace Transforms Chapter 3 Laplace Transforms Important analytical method for solving linear ordinary differential equations. - Application to nonlinear ODEs? Must linearize first. Laplace transforms play a key role in important

More information

EXAM 2 MARCH 17, 2004

EXAM 2 MARCH 17, 2004 8.034 EXAM MARCH 7, 004 Name: Problem : /30 Problem : /0 Problem 3: /5 Problem 4: /5 Total: /00 Instructions: Please write your name at the top of every page of the exam. The exam is closed book, closed

More information

Modes and Roots ... mx + bx + kx = 0. (2)

Modes and Roots ... mx + bx + kx = 0. (2) A solution of the form x(t) = ce rt to the homogeneous constant coefficient linear equation a x (n) + (n 1). n a n 1 x + + a 1 x + a 0 x = 0 (1) is called a modal solution and ce rt is called a mode of

More information

REVIEW NOTES FOR MATH 266

REVIEW NOTES FOR MATH 266 REVIEW NOTES FOR MATH 266 MELVIN LEOK 1.1: Some Basic Mathematical Models; Direction Fields 1. You should be able to match direction fields to differential equations. (see, for example, Problems 15-20).

More information

Chapter 7. Homogeneous equations with constant coefficients

Chapter 7. Homogeneous equations with constant coefficients Chapter 7. Homogeneous equations with constant coefficients It has already been remarked that we can write down a formula for the general solution of any linear second differential equation y + a(t)y +

More information