Notes LT1
Notes LT1
Definition and properties of Laplace Transform, piecewise continuous functions, the Laplace Transform method of solving initial value problems The method of Laplace transforms is a system that relies on algebra (rather than calculus-based methods) to solve linear differential equations. While it might seem to be a somewhat cumbersome method at times, it is a very powerful tool that enables us to readily deal with linear differential equations with discontinuous forcing functions.
Definition: Let f (t) be defined for t 0. The Laplace transform of f (t), denoted by F(s) or L{f (t)}, is an integral transform given by
st L {f (t)} = F ( s ) = e f (t ) dt . 0
Provided that this (improper) integral exists, i.e. that the integral is convergent. For functions continuous on [0, ), the above transformation is one-to-one. That is, different continuous functions will have different transforms.
C-1 - 1
1 , s
s > 0.
1 st e s
0
L{f (t)} =
e st f (t ) dt = e st dt =
0
1 1 1 0 e0 = ( 1) = = F ( s ) . s s s
1 , s2
s > 0.
1 , sa
s > a.
1 ( as )t e as
0
L{f (t)} =
e st e at dt = e ( a s )t dt =
0
1 1 1 0 e0 = ( 1) = = F ( s) . as as sa
C-1 - 2
Definition: A function f (t) is called piecewise continuous if it only has finitely many (or none whatsoever a continuous function is considered to be piecewise continuous!) discontinuities on any interval [a, b], and that both one-sided limits exist as t approaches each of those discontinuity from within the interval. The last part of the definition means that f could have removable and/or jump discontinuities only; it cannot have any infinity discontinuity.
Theorem: Suppose that 1. f is piecewise continuous on the interval 0 t A for any A > 0. 2. f (t) K e at when t M, for any real constant a, and some positive constants K and M. (This means that f is of exponential order, i.e. its rate of growth is no faster than that of exponential functions.) Then the Laplace transform, F(s) = L{f (t)}, exists for s > a.
Note: The above theorem gives a sufficient condition for the existence of Laplace transforms. It is not a necessary condition. A function does not need to satisfy the two conditions in order to have a Laplace transform. Examples of such functions that nevertheless have Laplace transforms are logarithmic functions and the unit impulse function.
C-1 - 3
1. L {0} = 0 2. L {f (t) g(t)} = L {f (t)} L {g(t)} 3. L {c f (t)} = c L {f (t)}, for any constant c.
Properties 2 and 3 together means that the Laplace transform is linear.
or, equivalently
L {t f (t)} = F (s)
Example: L {t } = (L {t}) =
d 1 2 2 = = 3 ds s 2 s3 s
at
C-1 - 4
Exercises C-1.1: 1 5 Use the (integral transformation) definition of the Laplace transform to find the Laplace transform of each function below. 2 6t 1. t 2. t e 3. cos 3t
it
4. e sin 2t
5.* e , where i and are constants, i = 1 . 6 8 Each function F(s) below is defined by a definite integral. Without integrating, find an explicit expression for each F(s). [Hint: each expression is the Laplace transform of a certain function. Use your knowledge of Laplace Transformation, or with the help of a table of common Laplace transforms to find the answer.] 6. 8.
e ( s +7 ) t dt 4e st sin 6t dt
7.
t 2 e ( s 3) t dt
Answers C-1.1: 1. 3.
2 s3
s s +9 s +i 2 5. 2 2 s + s + 2
2
1 ( s 6) 2 2 4. 2 s + 2s + 5
2.
Note: Since the Eulers formula says that e = cos t + i sin t, therefore, L{e it} = L{cos t + i sin t}. That is, the real part of its Laplace transform corresponds to that of cos t, the imaginary part corresponds to that of sin t. (Check it for yourself!) 6.
1 s+7
it
7.
2 ( s 3) 3
8.
24 s + 36
2
C-1 - 5
L {f (t)} = s2 L {f (t)} s f (0) f (0), L {f (t)} = s3 L {f (t)} s2 f (0) s f (0) f (0), L {f (n)(t)} = s n L {f (t)} s n 1 f (0) s n 2 f (0)
s2 f
(n 3)
: :
(0) s f
(n 2)
(0) f
(n 1)
(0).
This is an extremely useful aspect of the Laplace transform: that it changes differentiation with respect to t into multiplication by s (and, as seen a little earlier, differentiation with respect to s into multiplication by t, on the other hand). Equally importantly, it says that the Laplace transform, when applied to a differential equation, would change derivatives into algebraic expressions in terms of s and (the transform of) the dependent variable itself. Thus, it can transform a differential equation into an algebraic equation. We are now ready to see how the Laplace transform can be used to solve differentiation equations.
C-1 - 6
1. Take the Laplace transforms of both sides of an equation. 2. Simplify algebraically the result to solve for L{y} = Y(s) in terms of s. 3. Find the inverse transform of Y(s). (Or, rather, find a function y(t) whose Laplace transform matches the expression of Y(s).) This inverse transform, y(t), is the solution of the given differential equation.
The nice thing is that the same 3-step procedure works whether or not the differential equation is homogeneous or nonhomogeneous. The first two steps in the procedure are rather mechanical. The last step is the heart of the process, and it will take some practice. Lets get started.
C-1 - 7
C-1 - 8
Example:
y 6y + 5y = 0,
y(0) = 1, y(0) = 3
s 9 L{y} = 2 s 6s + 5
[Step 3] Find the inverse transform y(t) Use partial fractions to simplify,
s9 a b = + L{y} = s 2 6 s + 5 s 1 s 5
s 9 a ( s 5) b( s 1) = + s 2 6 s + 5 ( s 1)( s 5) ( s 5)( s 1)
s 9 = a ( s 5) + b( s 1) = ( a + b) s + ( 5a b)
C-1 - 9
1=a+b 9 = 5a b
Hence,
a=2 b = 1
2 1
L{y} = s 2 6 s + 5 = s 1 s 5 .
The last expression corresponds to the Laplace transform of 2e t e 5t. Therefore, it must be that
s9
y(t) = 2e t e 5t.
Many of the observant students no doubt have noticed an interesting aspect (out of many) of the method of Laplace transform: that it finds the particular solution of an initial value problem directly, without solving for the general solution first. Indeed, it usually takes more effort to find the general solution of an equation than it takes to find a particular solution! The Laplace Transform method can be used to solve linear differential equations of any order, rather than just second order equations as in the previous example. The method will also solve a nonhomogeneous linear differential equation directly, using the exact same three basic steps, without having to separately solve for the complementary and particular solutions. These points are illustrated in the next two examples.
C-1 - 10
Example:
y + 2y = 4t e2t,
y(0) = 3.
}=
4 ( s + 2) 2
4 ( s + 2) 2
(s + 2) L{y} + 3 =
4 ( s + 2) 2
(s + 2) L{y} =
4 3 ( s + 2) 2
4 3 4 3( s + 2) 2 3s 2 12 s 8 L{y} = ( s + 2) 3 s + 2 = ( s + 2) 3 = ( s + 2) 3
3s 2 12 s 8 a b c = + + L{y} = ( s + 2) 3 ( s + 2) 3 ( s + 2) 2 s + 2 .
C-1 - 11
3s 2 12 s 8 a b( s + 2) c ( s + 2) 2 = + + ( s + 2) 3 ( s + 2) 3 ( s + 2) 3 ( s + 2) 3 a + bs + 2b + cs 2 + 4cs + 4c cs 2 + (b + 4c ) s + ( a + 2b + 4c ) = = ( s + 2) 3 ( s + 2) 3
3 = c
12 = b + 4c 8 = a + 2b + 4c
a=4 b=0 c = 3
3s 2 12 s 8 4 3 = L{y} = ( s + 2) 3 ( s + 2) 3 s + 2 .
This expression corresponds to the Laplace transform of 2t2 e2 t 3e 2t. Therefore,
n! Note: L {t e } = ( s a ) n +1
n at
C-1 - 12
Example:
y 3 y + 2 y = e 3t ,
y(0) = 1, y(0) = 0
(s 3 s + 2) L{y} s + 3 = 1 / (s 3)
2 1 ( s 3 ) +1 2 = (s 3 s + 2) L{y} = s 3 + s 3 s 3
Therefore, y (t ) =
C-1 - 13
For the next example, we will need the following Laplace transforms:
s>0
s>0
s>a
s>a
Note: The values of a and b in the last two expressions denominators can be determined without using the method of completing the squares. Any irreducible quadratic polynomial s2 + Bs + C can always be written in the required from of (s a)2 + b2 by using the quadratic formula to find (necessarily complex-valued) s. The value a is the real part of s, and the value b is just the absolute value of the imaginary part of s. That is, if s = i, then a = and b = .
C-1 - 14
Example:
y 2y + 2y = cos(t),
y(0) = 1, y(0) = 0
s 3 2s 2 + 2s 2 1 s2 4s 6 + = L{y} = ( s 2 + 1)( s 2 2 s + 2) 5 s 2 + 1 s 2 2 s + 2
1 s 2 4( s 1) 2 = 2 2 + 2 2 5 s + 1 s + 1 s 2s + 2 s 2s + 2
which corresponds to
y (t ) =
]
C-1 - 15
2s 5 s 2 + 4s + 8
F (s) =
2s 5 2( s + 2) 9 = ( s + 2) 2 + 2 2 ( s + 2) 2 + 2 2 ( s + 2) 2 + 2 2
9 2t e sin( 2t ) 2
2t Answer: f (t ) = 2e cos( 2t )
s+4 F ( s ) = (ii) ( s 2) 3
2t 2 2t Answer: f (t ) = te + 3t e
C-1 - 16
Appendix A
Some Additional Properties of Laplace Transforms I. Suppose f (t) is discontinuous at t = 0, then the Laplace transform of its derivative becomes
L {f (t)} =
1 1 e sT
e st f (t ) dt .
Comment: This property is convenient when finding the Laplace transform of a discontinuous periodic function, whose discontinuities (necessarily there are infinitely many, due to the periodic nature of the function) would make the usual approach of integrating from 0 to unwieldy. III. Let c > 0 be a constant, the time-scaling property of Laplace transform states that
L {f (c t)} = F ( ).
IV. We have known that L {t f (t)} = F (s). Taking the inverse transforms on both sides yields t f(t) = L 1{F (s)}. Therefore, f (t) =
1 1 L {F (s)}. t
1 c
s c
C-1 - 17
That is, if we know how to invert the function F (s), then we also know how to find the inverse of its anti-derivative F(s). Formally,
f (t ) L{ }= t
F (u ) du .
L { 0 f (u ) du } =
F (s) . s
Comment: Therefore, dividing a function by its independent variable has the effect of anti-differentiation with respect to the other independent variable for Laplace transforms. These two properties (IV and V) are the counter parts of the multiplication-corresponds-to-differentiation properties seen earlier. VI. After learning the fact that L {f (t) g(t)} L{f (t)} L{g(t)}, one might have wondered whether there is an operation of two functions f and g whose result has a Laplace transform equal to the product of the individual transforms of f and g. Such operation does exist:
L { 0 f ( ) g (t ) d } = L{f (t)}L{g(t)}.
The integral on the left is called convolution, usually denoted by f * g (the asterisk is the convolution operator, not a multiplication sign!). Hence,
L{f * g} = L{g * f }.
C-1 - 18
Exercises C-1.2: 1 8 Find the Laplace transform of each function below. 1. f (t) = t3 t2 + 5t + 2 2. f (t) = 2cos 6t + 5sin 6t 3. f (t) = 3et 4e2t + 2e4t 4. f (t) = 7t + 6et 2et 10 5. (a) f (t) = sin 2t sin 3t 6. (a) f (t) = cos2 5t 7. (a) f (t) = t e at cos bt, 8. f (t) = t3 cos 2t 9 14 Find the inverse Laplace transform of each function below. 1 4s + 2 9. F ( s ) = 2 10. F ( s ) = ( s 2)( s 4)( s 8) s + 6 s + 34 11. F ( s ) =
3s 7 4s 2 + 1
(b) f (t) = sin 2t cos 2t (b) f (t) = t cos2 5t (b) f (t) = t e at sin bt
12. F ( s ) =
s2 + s 6 s 3 + 2s 2 + s
1 13. F ( s ) = 4 s 81
s3 14. F ( s ) = 4 s 16
15 27 Use the method of Laplace transforms to solve each IVP. 15. y + 10y = t2, y(0) = 0 16. y + 2y = t et, 17. y + 2y = t e t , 18. y + 4y 5y = 0, 19. y 2y + y = 2t 3,
2008 Zachary S Tseng
20. y + 8y + 25y = 13e2t, 21. y + 6y + 34y = 0, 22. y + 4y = 8cos 2t 8e2t, 23. y 4y + 4y = 8t2 16t + 4, 24. y 4y 5y = 3t3, 25. y + 3 y + 3 y + y = 0, 26. y + 4 y 5 y = 0, 27. y y + 4y 4y = 26e3t,
28. Prove the time-scaling property (property III, Appendix A). 1 vu Hint: Let u = c t, and v = s / c, then show L {f (c t)} = 0 e f (u ) du . c 29. Use property IV of Appendix A to verify that
1 t
Answers C-1.2:
2 s 3 + 5s 2 2 s + 6 1. F ( s ) = s4 s 2 12 s 4 F ( s ) = 3. ( s 1)( s 2)( s + 4) s s 5. (a) F ( s ) = 2 2 2( s + 1) 2( s + 25)
2. F ( s ) =
2 s + 30 s 2 + 36 6 s 3 + 15 s 2 + 10 s 7 F ( s ) = 4. s4 s2 2 (b) F ( s ) = 2 s + 16
C-1 - 20
s 1 1 6. (a) F ( s ) = + 2 , 2 s s + 100
7. (a) F ( s ) = 8. F ( s ) =
( s a) 2 b 2 , (( s a ) 2 + b 2 ) 2
(b) F ( s ) =
2b( s a ) (( s a ) 2 + b 2 ) 2
6 s 4 144 s 2 + 96 ( s 2 + 4) 4 9. f (t) = 4e3t cos 5t 2e3t sin 5t 1 2 t 1 4 t 1 8t 10. f (t ) = e e + e 12 8 24 3 t 7 t 11. f (t ) = cos sin 4 2 2 2 t t 12. f (t) = 7e + 6t e 6 1 3t 1 3t 1 e e sin 3t 13. f (t ) = 108 108 54 1 2t 1 2t 1 14. f (t ) = e + e cos 2t 4 4 2 2 t t 1 1 10t e 15. y = + 10 50 500 500 16. y = t et et + 3e2t 17. y = t e t + 6 18. y = 4et + e5t 19. y = 4et + 5t et + 2t + 1 20. y = 2 e4t cos 3t + 4 e4t sin 3t + e2t 21. y = e 3t cos 5t + 2e 3t sin 5t 22. y = cos 2t sin 2t + 2t sin 2t e2t 23. y = t e2t + 2t2 24. y = e5t + 2et + 3t3 25. y = 7e t + 2t2 e t 26. y = 5 2e t + e 5t 27. y = 4et + e3t + cos 2t + 2sin 2t
C-1 - 21