CH.6 Laplace Transform

Transcription

1 CH.6 Laplace Transform Where does the Laplace transform come from? How to solve this mistery that where the Laplace transform come from? The starting point is thinking about power series. The power series looks like n= f n x n = f + f x + f 2 x 2 + f 3 x 3 + This series is the summation of x s with increasing positive integer powers. assigned numbers at n =,, 2,. If we write it using a computer notation. The coefficients here are n= f[nx n = f[ + f[x + f[2x 2 + f[3x 3 + Next, we would like to investigate the relation of F(x) and the coefficients. coefficients f[n = [f[, f[, f[2,, how to calculate F(x)? For example, if f[n is a unit vector, (e.g. f[n = [,,, ) Specifically, for different n= x n = x except for this is the wrong answer! It is not true for every value of x. It is only true for the series to converge. This answer is true, if and only if its suport set is restricted to x <. As such n= x n = x, x < If f[n = n! results in f[n = [,, 2!, 3!,, the power series looks like n= n! xn = e x This identity here is valid for any arbitrary value x I am not satisfied with this discrete representation. I want more! I want a continuous analogy. The discrete variable n =,, 2, are replaced by a continuous variable t. And, in the continuous case, summation is approximated by integral. + f(t)x t dt

2 In general, we do not not like the variable x appearing at the base,making the integral very difficult to compute. How to fix this? + f(t)x t dt = + f(t) ( e ln x) t dt One problem solved, but here comes another two. Given any value x, can I calculate the integral? For any x, since t goes to infinity, the integral is quite unlikely to converge. Then, we have to restrict that x <. Another problem is if < x <, since t can value from to, a negative number with fractional power may give you a complex number, which is very undesired. For example 2 = = j Hence, we have to restrict that < x <. Look at the blue expression again. If x (, ), we have ln x (, ). The integral is taken over all the positive numbers, but the power is negative??? How confusing! Then, we apply variable replacement s = ln x (, + ) As such, everything is all set. Then, what happened to the integral? F(s) = This is the unilateral Laplace transform of f(t). + f(t)e st dt 2 Calculate Laplace Transform Laplace transform is defined over the whole complex plane (e.g. the variable in Laplace transform s = σ+jω is a complex number). Fourier transform is a special case, when σ = if and only if the ROC of Laplace transform contains the imagine axis. In your text book, the author stated Laplace transform of a signal may exist when the Fourier transform of the same signal does not. This can be explained by the ROC of Laplace transform does not include the imagine axis. As such, we can calculate Laplace transform following similar steps in carrying out the Fourier transform. In addition, we need to discuss the region of convergence of corresponding Laplace transform. Example.Calculate the Laplace transform of the following signals (a). x(t) = cos(ω t)u(t) (b). x(t) = t e 2τ u(τ)dτ (c). x(t) = e at u( t + b), where a is a positive real constant and b is a real constant Solution (a). By definition, the Laplace transform of x(t) is equivalent to solve the following integral L [cos(ω t)u(t) = cos(ω t)u(t)e st dt = cos(ω t)e st dt = ( e jω t + e ) jω t e st dt s = σ + jω 2 [ = + 2 jω σ jω ej(ω ω)t e σt + jω σ jω ej( ω ω)t e σt + 2

3 I will use the first part as an illustration to discuss the ROC when you compute Laplace transform by definition. + jω σ jω ej(ω ω)t e σt = jω σ jω ej(ω ω)+ e σ+ e j(ω ω)t is a complex number with magnitude one. As t grows bigger, this complex number runs along unit ring anti-clock-wise. I do not care the exact value of it when t approaches infinity. I just know, the limit, e j(ω ω), is still a complex number on this unit ring. And it is not infinity. e σt is another story. The sketch or the limit at t = depends on σ as illustrated in the figure bellow If σ <, σ >. As t goes to infinity, the limit is also infinity, which means the integral does not converge. The limit exists only when σ >, which is equivalent to Re[s >. Correspondingly, the limit e σ is when Re[s >. Restrict the real part of s, the integral can be carried out as L [cos(ω t)u(t) = [ 2 jω σ jω + jω σ jω = [ 2 jω + s jω s s =, Re[s > s 2 + ω 2 (b). x(t) = t f (τ)dτ, where f (t) = e 2t f 2 (t) and f 2 (t) = u(t). Apply Time domain integration, we have Apply Laplace domain shifting, we obtain Laplace transform of f 2 (t) is specified as L x (s) = s L f (s), R x = R f Re[s > L f (s) = L f2 (s + 2), R f = R f2 2 L f2 (s) = L [u(t) = s, R f 2 = Re[s > 3

4 Now, we are ready to characterize the Laplace transform of x(t) The ROC can be determined by L x (s) = s L f (s) = s L f 2 (s + 2) = s(s + 2) R x = R f Re[s > = [R f2 2 Re[s > = Re[s > 2 Re[s > = Re[s > The final answer is L x (s) = s(s + 2), Re[s > (c). x(t) = e at f (t), where f (t) = f 2 ( t), f 2 (t) = f 3 (t + b) and f 3 (t) = u(t) Apply Laplace domain shifting and linearity, we have L x = L f (s a), R x = R f + a Apply Time domain scaling, we have L f (s) = L f2 ( s), R f = R f2 Apply Time domain shifting, we have L f2 (s) = e bs L f3 (s), R f2 = R f3 The Laplace transform of f 3 (t) can be calculated as L f3 (s) = s, R f 3 = Re[s > By substituting each equation to the previous one, we obtain The ROC of x(t) is determined as L f2 (s) = e bs s, L f = e bs s, L x = e b(s a) s a R f = Re[s <, as such, R x = Re[s < a The final result is L x (s) = e b(s a) s a, Re[s < a As I stated, calculation Laplace transform is similar to calculating Fourier transform, except that we need to determine the ROC of corresponding Laplace transform. 3 Inverse Laplace Transform Calculate inverse Laplace transform by the definition is torturing as the contour integration is not easily to solve. However, with partial fraction expansion, we can construct a lot of Laplace transform pairs. Example.2 Find all possible inverse Laplace transform of H(s) = 7s s 2 4

5 Solution With partial fractional expansion where the coefficients are H(s) = 7s s 2 = 7s (s )(s + ) = A s + A 2 s + A = (s )H(s) = 7s = 3 s + s = s= A 2 = (s + )H(s) = 7s = 4 s s+= s= We obtain this expansion as H(s) = 3 s + 4 s + We can easily figure out that H(s) has two poles s = and s =. The possible ROC should be Re[s <, Re[s > or < Re[s <. From the Laplace transform table, we note that ( ) { e at L u(t), Re[s > a = s + a e at () u( t), Re[s < a For Re[s <, We note that ( ) h(t) = 3L s = [ 3e t 4e t u( t), Re[s < ( ) + 4L = 3[ e t u( t) + 4[ e t u( t) s + I will solve this by baby-steps. Hopefully, you can figure out how to proceed for the rest of the two cases. ( ) L, Compare with Eq. (), a =, a = ; s Accordingly, The ROC is Re[s<, of course Re[s< Re[s< a ( ) Then, L = e t u( t) Re[s < s ( ) L, Compare with Eq. (), a =, a = ; s + The ROC is Re[s< Re[s< a ( ) Then, L = e t u( t) Re[s < s + Substitute these two inverse Laplace transform back to the expression of h(t), we got the result in Eq. (2). For Re[s >, We note that ( ) ( ) h(t) = 3L + 4L = 3[e t u(t) + 4[e t u(t) s s + = [3e t + 4e t u(t), Re[s > For < Re[s < ( ) h(t) = 3L s ( ) + 4L = 3[ e t u( t) + 4[e t u(t) s + = 3e t u( t) + 4e t u(t), < Re[s < This example tells us, ROC means a lot!!!the mathematical expression combined with its ROC defines a Laplace transform pair. 5 (2)

6 4 LTI system analysis I will skip the differential equation part. We can solve it following exactly identical steps using Fourier analysis. 4. bilateral Laplace transform for all zeros initial condition In your coming study, courses like ELEC 36, control theory, focuses on characterizing system input-output relation with Laplace transform. Let us see how to proceed Example.3 Consider the LTI system with input x(t), output y(t), and system function H(s), as shown in the figure below. Suppose that H and H 2 are causal LTI systems with system functions H (s) and H 2 (s), respectively, given by H (s) = s, H 2(s) = A (a). Find an expression for H(s) in terms of H (s) and H 2 (s). (b). Determine for what values of A the system is BIBO stable. SolutionFrom the diagram, we have the following equations y(t) = v(t) h (t), v(t) = x(t) + y(t) h 2 (t) Take Laplace transform on both sides L y (s) = L v (s)l h (s), L v (s) = L x (s) + L y (s)l h2 (s) 6

7 Replace L v (s) by L x (s) + L y (s)l h2 (s), we obtain [ L y (s) = L x (s) + L y (s)l h2 (s) L h (s) We obtain the system function as = L x (s)l h (s) + L y (s)l h2 (s)l h (s) [ L h2 (s)l h (s) L y (s) = L x (s)l h (s) H(s) = Y (s) X(s) = H (s) H (s)h 2 (s) = s = A s s (A + ) Since the two sub-system H and H 2 are causal. we have the ROC as Re[s > A +. (b). System is stable if and only if the ROC includes imagine axis. To guarantee Re[s > A + include the imagine axis, A + <. We have A < 4.2 unilateral Laplace transform for non-zero initial conditions Example.4 Suppose that we have a incrementally-linear TI system with input x(t) and output y(t) characterized by the differential equation y (t) + 7y (t) + 2y(t) = x(t) If y( ) =, y ( ) = and x(t) = u(t), find y(t) Solution Take unilateral Laplace transform on both sides U L [y (t) + 7U L [y (t) + 2U L [y(t) = U L [x(t) From the unilateral Laplace transform table, we have U L [y (t) = s 2 Y (s) sy( ) y ( ) U L [y (t) = sy (s) y( ) U L [y(t) = Y (s) U L [x(t) = U L [u(t) = s I have to stress, you can check out everything in your text book except fot U L [y (t). Denote z(t) = y (t), the unilateral Laplace transform of z(t) is Apply the time domain differentiation, we have Z[s = U L [y (t) = sy (s) y( ) U L [y (t) = U L [z (t) = sz(s) z( ) = s[sy (s) y( ) y ( ) = s 2 Y (s) sy( ) y ( ) Substituting all the unilateral Laplace transform in differential equation and after several carrying out manipulations, the Laplace transform of output y(t) is given by Y (s) = s2 7s + s(s + 3)(s + 4) 7

8 The partial fractional expansion of Y (s) is Y (s) = A s + A 2 s A 3 s + 4 Solve the coefficients by A = sy (s) = s2 7s + = (s + 3)(s + 4) 2 s= s= A 2 = (s + 3)Y (s) = s2 7s + = 3 s(s + 4) 3 s= 3 s= 3 A 3 = (s + 4)Y (s) = s2 7s + = 3 s(s + 3) 4 s= 4 s= 4 The Laplace transform can be calculated as Y (s) = ( ) 2 s 3 3 ( ) + 3 ( ) s s + 4 y(t) can be calculated by taking inverse unilateral Laplace transform y(t) = [ 2 U L 3 [ s 3 U L s + 3 = 3 u(t) 2 3 e 3t u(t) e 4t u(t) U L [ s + 4 8