Laplace Transform Methods - University of North Florida

1 CHAPTER 1. Laplace Transform Methods Laplace Transform is a method frequently employed by engineers. By applying the Laplace Transform , one can change an ordinary dif- ferential equation into an algebraic equation, as algebraic equation is generally easier to deal with. Another advantage of Laplace Transform is in dealing the external force is either impulsive , (the force lasts a very shot time period such as the bat hits a baseball) or the force is on and off for some regular or irregular period of time. 1. The Laplace Transform If f (t) is defined over interval [0, ), the Laplace Transform of f , denoted as f (s), is Z . b L(f ) = f (s) = e st f (t) dt 0. Our first theorem states when Laplace Transform can be performed, Theorem If f (t) is (piecewise) continuous and there are pos- itive numbers M, a such that |f (t)| M eat for all t c Then fb(s) is defined for all s > c The next result shows that Laplace Transform is unique in the sense that different continuous functions will have different Laplace trans- form.]

2 Theorem If fb(s) = gb(s) for all s > c, then f (t) = g(t) at all t where both are continuous. Notice if f (t) and g(t) are piecewise continuous (continuous except at finite points where left and right limits exists,) their Laplace trans- forms can be same for all s > c even if they are different at the isolated discontinuous point. Since solutions of ordinary equations must be continuous, so this is of no important concern. 1. 2 1. Laplace Transform Methods . Due the uniqueness, we can define the inverse Laplace Transform L 1 as L 1 (fb)(t) = f (t). Theorem If both fb(s) and gb(s) exist for all s > c, then af (t)+. bg(t) has Laplace Transform for all constant a and b and \. af + bg(s) = afb(s) + bb g (s), for all s > c So to find Laplace Transform of summation, we just need to find Laplace Transform of each term.

3 The next result shows that Laplace Transform changes derivative into scalar multiplication, it is this property enable Laplace Transform to change ODE into algebraic equation. Theorem Suppose that the function f (t) is continuous and piecewise smooth (f 0 (t) is piecewise continuous) for all t 0 and there are constants M, a such that |f (t)| M eat for t T, Then fb0 (s) is defined for all s > a and fb0 (s) = sfb(s) f (0). Now if nth derivative f n (t) is piecewise smooth, then fd (n) (s) = sn fb(s) sn 1 f (0) sn 2 f 0 (0) f (n 1) (0). For example, When n = 2. fd (2) (s) = s2 fb(s) sf (0) f 0 (0). When n = 3. fd (3) (s) = s3 fb(s) s2 f (0) sf 0 (0) f 00 (0). When n = 4. fd (4) (s) = s4 fb(s) s3 f (0) s2 f 0 (0) sf 00 (0) f (3) (0). The following table list most commonly used functions, and . 0 if t<a ua (t) = u(t a) =.

4 1 if t a is the Heaviside function. 1. THE Laplace Transform 3. f = L 1 (fb) fb f = L 1 (fb) fb 1 1. 1 s (s > 0) t s2. (s > 0). n! (a+1). tn (n 0) sn+1. (s > 0) ta (a > 1) sa+1. (s > 0). 1 k eat s a (s > 0) sin(kt) s2 +k2. (s > 0). s k cos(kt) s2 +k2. (s > 0) sinh(kt) s2 k2. (s > |k|). s eas cosh(kt) s2 k2. (s > |k|) u(t a) s (s > 0). b(s) if x00 +2x0 3x = et , x(0) = 1, x0 (0) = 1. Example Find x Solution Apply Laplace Transform on both side of the equation, L(x00 + 2x0 3x)(s) = L(et )(s). \. Using the linear property af + bg(s) = afb(s) + bb g (s) we have L(x00 + 2x0 3x)(s) = xb00 (s) + 2xb0 (s) 3b x(s). Together with xb00 (s) = s2 x b(s) sx(0) x0 (0) and xb0 (s) = sb x(s) x(0). 0. We have, due to x(0) = 1, x (0) = 1, L(x00 + 2x0 3x)(s) = xb00 (s) + 2xb0 (s) 3b x(s). 2. = [s x b(s) sx(0) x0 (0)] + 2[sb x(s) x(0)] 3b x(s).

5 2. = (s + 2s 3)b x(s) s 1. So we have an algebraic equation for x b(s), 1. x(s) s 1 = ebt (s) =. (s2 + 2s 3)b s 1. Solve this equation for x b(s), x(s) s 1 = ebt (s) =. (s2 + 2s 3)b 1. s 1. 1. (s2 + 2s 3)b x(s) = s 1 +s+1. 2. (s + 2s 3)bx(s) = s 1 + (s+1)(s 1). 1. s 1. x(s) = 1+(s+1)(s 1). (s2 + 2s 3)b s 1. 2 1+(s2 1). (s + 2s 3)bx(s) = s 1. 2 s2. (s + 2s 3)bx(s) = s 1. 4 1. Laplace Transform Methods . we get s2. x b(s) = (s 1)(s2 +2s 3). a 2. Using Mathcad Mathcad can help us in find both Laplace Transform and inverse Laplace Transform . To use Mathcad to find Laplace Transform , we first enter the expres- sion of the function, then press [Shift][Ctrl][.], in the place holder type the key word Laplace followed by comma(,) and the variable name. For example, to find the Laplace of f (t) = t2 sin(at), you firs enter the expression t2 sin(at) by typing, t^2*sin(a*t), then press [Shift][Ctrl][.]

6 ], and entering keyword Laplace followed by comma(,) and t, you will get 8as2 2a t2 sin(at) Laplace , t 2 2 3. 2. (s + a ) (s + a2 )2. If you want a simplified result, you type comma(,) after entering the variable name and keyword simplify, so the following input t^2*sin(a*t)[Shift][Ctrl][.] Laplace ,t,si mplify will produce the result 2a(a2 3s2 ). t2 sin(at) Laplace , t, simplif y . (s2 + a2 )3. Similarly, to use Mathcad to find inverse Laplace Transform , we first enter the expression, then press [Shift][Ctrl][.], in the place holder type the key word invlaplace followed by comma(,) and the variable name. 2. For example, to find the inverse Laplace of (s28as +a2 )3. , you firs enter the 8as2. expression (s2 +a2 )3. as, 8a*s^2/(s^2+a^2)^3, then press [Shift][Ctrl][.], and entering keyword invlaplace followed by comma(,) and s, and for simplifying the result, by another comma(,).

7 And keyword simplify, we get . 8as2 1 2 2. invlaplace, s, simplif y 2 sin(at) ta cos(at)+t a sin(at). (s2 + a2 )3 a 2. USING Mathcad 5. When using Mathcad together with Laplace Transform to solve an ODE. an x(n) + an 1 x(n 1) + + a1 x0 + a0 x = f (t). we follow these steps, Step One: Apply Laplace to both sides of equation. Using Mathcad to find Laplace Transform of f (t). Step Two: Using the linear property and fd (n) (s) = sn fb(s) sn 1 f (0) sn 2 f 0 (0) f (n 1) (0). to find an algebraic equation for x b(s), Step Three: Solve the equation obtained in Step Two for x b(s) and using Mathcad to find inverse Transform which will be the solution x(t). Example Find general solution to x00 + 2x0 + 3x = t2 et . Solution Since we want to find general solution, we set x(0) = a and x0 (0) = b Step One Apply Laplace Transform to both sides of the equa- tion and find Laplace Transform for t2 et.

8 L(x00 + 2x0 + 3x)(s) = L(t2 e2 )(s). Type, t^2*e^t[Shift][Ctrl][.] Laplace ,t,simplif y we find 2. L(t2 e2 )(s) =. (s 1)3. Step Two: Using linear property to find an equation for x b(s). L(x00 + 2x0 + 3x)(s) = (xb00 + 2xb0 + 3b x(s). = s2 x b(s) sa b + 2(sb x(s) a) + 3b x(s). 2. = (s + 2s + 3)b x(s) sa b 2a The equation for x b(s) is 2. (s2 + 2s + 3)b x(s) sa b 2a = (s 1)3. Hence, sa+b+2a 2. x b(s) = (s2 +2s+3). + (s 1)3 (s2 +2s+3). Step Three: Using Mathcad to find inverse Laplace trans- form and x(t), we enter, s*a+b+a/(s^2+2s+3)+2/(s-1)^3(s^2+2s+3)[S hift][Ctrl][.] Laplace ,s 6 1. Laplace Transform Methods . and get . 2.. x(t) = ae t sin( 2t) + cos( 2t) + b 22 e t sin( 2t). 2.. +et 3 t + 9 + 9 cos( 2t) + 182 e t sin( 2t). 1 2 2. a Solving System of equations. We can use Laplace trans- form method to solve system of differential equations.)

9 The procedure is the same as solving a higher order ODE . But, after applying Laplace Transform to each equation, we get a system of linear equations whose unknowns are the Laplace Transform of the unknown functions. The following example shows how we can use Laplace method with Mathcad to solve system of differential equations. Example Find solution to x0 = 3x 4y + t2. y 0 = x + 5y + et x(0) = 1, y(0) = 1. Solution Step One: Apply Laplace Transform to both sides of the equa- tions and use Mathcad to find Laplace Transform of t2 and et , L(x0 ) = L(3x 4y + t2 ). L(y 0 ) = L( x + 5y + et ). and L(t2 )(s) = s63 and L(et )(s) = s 1. 1. , Step Two: Apply linear property to get a system of equations for xb(s) and yb(s), due to L(x0 ) = sbx(s) x(0) and L(y 0 ) =. sb y (s) y(0), sb x(s) 1 = 3b x(s) 4by (s) + s63. 1. sb y + 1 = bx(s) + 5b y (s) + s 1.

10 From which we get (s 3)bx(s) + 4by (s) = 1 + s63. 1. x b(s) + (s 5)b y = 1 + s 1.. b x(s) s 3 4. Define xb (s) = ,A= , and b(s) =. yb(s) 1 s 5.. 1 + s63. 1 , we have, in matrix form, 1 + s 1. Ab x(s) = b(s). 3. THE PIECEWISE DEFINED FUNCTIONS 7. Step Three: Solve the algebraic equation and apply the in- verse Laplace Transform , using Mathcad we have " #. (s 5)(s3 +6) 2 s 1 4 s 1. b (s) = A-1 b(s) = 2. x s3. (s 3)(2 s) 3. s 8s + 11 s 1. s s+6. 3.. x(t). x(t) =. y(t) . 15 2 174 1062 4t 3742. 2334.. 11 t 121 t 1331 + e 1331. cosh( 5t) 6655 sinh( 5t) et . = . 3 2.. 11 t 121 t 1331 e4t 2662 cosh( 5t) + 1331 sinh( 5t) 12 et 48 318 695 8143. a 3. The Piecewise defined Functions When solving an ODE an x(n) + an 1 x(n 1) + + a1 x0 + a0 x = f (t). with multiple defined function f (t), we need to change f (t) into linear ua (t).