Lab 3 Laplace Transform - v3
Lab 3 Laplace Transform - v3
LAPLACE
LAPLACE TRANSFORM
Name:
__________________________________
ID #:
__________________________________
Section #:
__________________________________
Due Date:
__________________________________
Instructor
__________________________________
LAPLACE TRANSFORM
OBJECTIVES
Figure 1
Laplace Transform
G (s ) =
num
den
s +5
3
s + 6s 2 + 11s + 6
(1)
Solution:
In Eq. (1) num and den denotes respectively, the numerator and denominator of
G (s ) .
G (s ) = 0 . For
num = s + 5 = 0 s = 5 . Therefore s = 5
Zeros:
values
Poles:
values of
of
G (s ) = 0
or
for
which
for which
G (s )
G( s) = 0
is the zero of G (s ) .
instance,
den = s 3 + 6 s 2 + 11s + 6 = 0 ).
>> den=[1
6
>> roots(den)
ans=
-3.000
-2.000
-1.000
11
6];
G (s ) =
where the residues
methods:
G (s )
num
s +5
A
B
C
=
=
+
+
den (s + 3)(s + 2 )(s + 1) (s + 3) (s + 2 ) (s + 1)
A, B
and
(2)
B=
Laplace Transform
(s + 3)(s + 5)
(s + 3)(s + 2)(s + 1)
(s + 2)(s + 5)
(s + 3)(s + 2)(s + 1)
( 3 + 5) = (2) = 1
( 3 + 2 )( 3 + 1) ( 1)( 2)
( 2 + 5) = (3) = 3
( 2 + 3)( 2 + 1) (1)( 1)
s =3
s = 2
C=
Substituting
A, B
and
(s + 1)(s + 5)
(s + 3)(s + 2)(s + 1)
=
s =1
( 1 + 5) = (4) = 2
( 1 + 3)( 1 + 2) (2)(1)
G( s) =
Therefore, the expression of
1
3
2
+
( s + 3 ) ( s + 2 ) ( s + 1)
(3)
g (t ) is given by
g (t ) =
G ( s ) = e3t 3e2t + 2e t
(4)
g (t ) =
where ri , pi and
k (s)
L [G (s )] =
A, B
and
. Equation
r
r1
r
+ 2 + 3 + k (s )
s p1 s p2 s p3
(5)
are the residues, poles and direct terms of the partial fraction
num ( s )
and
den ( s ) .
Notice that
ri , pi
and
k (s)
1.0000
-3.0000
2.0000
p=
-3.0000
-2.0000
-1.0000
k=[]
Notice that k = [ ] means that
k (s )
G (s ) =
Laplace Transform
1
3
2
+
s + 3 s + 2 s +1
G (s ) =
num
1
= 2
den s (s + 2 )2
(1)
Solution:
Zeros:
Poles:
None
2
s =0
(s + 2)
s = 0, 0
2
s = 0 is a pole of multiplicity 2
= 0 s = 2, 2 s = 2
is a pole of multiplicity 2
G(s)
can be written in
G (s ) =
num
1
A A2
B1
B2
= 2
= 1+ 2 +
+
2
den s (s + 2 )
s
s
(s + 2 ) (s + 2 ) 2
(2)
A2 =
s2
s 2 (s + 2 )
2
s =0
( 0 + 2)
1
4
0 (s + 2 ) 2 (s + 2 )
d
s2
d 1
1
=
A1 = 2
=
=
2
4
2
ds s (s + 2 )
ds (s + 2 )
4
(s + 2 )
s =0
s =0
s =0
B2 =
(s + 2 )
s (s + 2 )
2
( 2 )
1
4
s =2
d (s + 2 )
B2 =
ds s 2 (s + 2 ) 2
G (s ) =
The expression of
d 1
0 s 2 2s
1
= 2
=
=
2
2
ds s s =2
(s ) s =2 4
s =2
A i and B i into Eq. (2) gives
num
1
1/ 4 1/ 4
1/ 4
1/ 4
= 2
=
+ 2 +
+
2
den s (s + 2 )
s
s
(s + 2 ) (s + 2 ) 2
(3)
g (t ) is given by
1 1
1
1
g (t ) = + t + e 2t + te 2t
4 4
4
4
(4)
MATLAB command
resi2(num,den,pole,n,k)
num = [1];
a = conv([1 0], [1
b = conv([1 2], [1
den = conv(a, b);
0]);
2]);
A1 = resi2(num,den,0,2,1)
>> A2 = resi2(num,den,0,2,2)
>> B1 = resi2(num,den,-2,2,1)
>> B2 = resi2(num,den,-2,2,2)
Laplace Transform
% the expression of s2
% the expression of (s+2)2
% the expression of s2 (s+2)2
A1
gives A2
gives B1
gives B2
gives
= 0.2500,
= 0.2500,
= 0.2500,
= 0.2500,
5
Substituting the above values into Eq. (2), one can get
G( s) =
num
1
1/ 4
1/ 4
1/ 4 1/ 4
=
+ 2 +
+
= 2
2
(s + 2) (s + 2)2
den
s
s
s (s + 2)
sinh(t ) =
et e t
2
and
cosh(t ) =
et + e t
2
1 2t 1 2t 1
e t e t + e t e t e t 1
g (t ) =
t
= t + te e
2
2
2 4
4
4
4
which is similar to Eq. (4).
G(s) =
where
Solution:
The expression of
Laplace Transform
Bs + C
s + 2as + a 2 + 2
2
(1)
G ( s ) can be written as
6
G (s) =
Bs + C
Bs + C
=
2
2
( s + 2as + a ) + ( s + a )2 + 2
2
(2)
It is clear that the denominator is the sum of two complete squares. The roots of
the denominator (poles) are:
s = a + j
1
_
s2 = a + j
2
(s + a) + 2 = 0
(3)
where
s1
s2
G (s ) =
G ( s ) is written as
Bs + C
2
B ( s + a ) + C Ba
2
(s + a ) +
(s + a ) + 2
B (s + a )
C Ba
=
+
2
( s + a ) + 2 ( s + a )2 + 2
=
B (s + a )
(s + a )
+ 2
C Ba
2
( s + a ) + 2
g (t ) = Be at cos (t ) +
Laplace Transform
(4)
g (t ) is given by
(C Ba ) e at sin t
( )
(5)
written as
G (s ) =
where
and
Bs + C
( s + a + j )( s + a j )
s + a + j
(6)
s + a j
(s + a + j)(Bs + C )
B ( a j) + C
(C Ba )
B
=
=
= + j
2
(s + a + j)(s + a j) s = a j ( a j + a j) 2
(7)
(s + a j)(Bs + C )
B( a + j) + C
(C Ba )
B
=
=
= j
2
(s + a + j)(s + a j) s = a + j ( a + j + a + j) 2
(8)
and
Notice that
and
gives
and
(C Ba ) B j (C Ba )
B
+ j
2 + 2
2
G(s) = 2
s + a + j
s + a j
(9)
B
(C Ba ) e ( a j ) t + B j (C Ba ) e ( a + j ) t
g (t ) = + j
2
2
2
2
Arranging terms, the above expression can be written in the form
g (t ) =
(C Ba ) e at e jt e jt
B at jt jt
e
e
+
e
2
2
B
g (t ) = e at
2
e jt + e jt
(C Ba ) e at
jj
2
2
e jt e jt
2j
(10)
Laplace Transform
cos t =
the expression of
e jt + e jt
e jt e jt
, sin t =
2
2j
(11)
g (t ) is written as:
g (t ) = Be at cos t +
(C Ba ) e at sin t
(12)
y + 3 y + 2 y = f (t ),
where
y (0) = 0, y (0) = 0
f (t )
(1)
shown in Fig. 3.
Solution:
t
1,
f (t ) =
0,
0 t <1
t 1
Figure 3 Graph of
f (t ) .
st
Alternatively,
st
=
0
(2)
f (t )
u ( t 1)
t
0
Figure 4
Laplace Transform
e s 1 1
= (1 e s )
s
s
t
1
Therefore,
f (t ) = u (t ) u (t 1)
(3)
F (s) =
1 es 1
= 1 e s
s
s
s
Step 2: Remember:
Zero Initial Conditions
L
L
L
y (0) = 0, y (0) = 0
L
L
L
{ y} = s 2 Y ( s)
{ y} = sY ( s)
{ y} = Y ( s)
(s
+ 3s + 2 ) Y ( s ) =
1
1 e s )
(
s
Y (s) =
1
1 e s ) = H ( s ) (1 e s )
(
2
s ( s + 3s + 2 )
where
H ( s) =
1/ 2
1
1/ 2
= +
+
=
+
s
s +1 s + 2
s ( s + 3s + 2 ) s s + 1 s + 2
2
Hence
h(t ) =
1
1
e t + e 2t
2
2
L [ h (t )u (t ) ] = e s H (s ) L
For this example
Laplace Transform
=1
e s H (s ) = h (t )u (t ) , 0
e s H ( s )
is given by
10
0 t <1
t >1
y (t ) = L
1 (t 1) 1 2(t 1)
0 t <1
e
+ e
2
[Y (s )] = f (t ) f (t 1)u (t 1) = 2
K e t K e 2t
t >1
2
1
where K 1 = e 1 and K 2 = (e 2 1) / 2 .
Laplace Transform
11
PART 2: ASSIGNMENT
Consider the following transfer function (Instructors are required to fill in the
blanks such that the poles are real and distinct)
G(s) =
num ( s )
den ( s )
s +
s + s2 + s +
3
G(s)
G(s) .
in partial fractions
num ( s )
A
B
C
+
+
den ( s ) s + p1 s + p2 s + p3
Where pi are the poles of G ( s ) and A , B and C are its residues.
By hand calculations, find A , B and C .
G( s) =
3.
4. Check your finding by using the MATLAB residue command (See help
residue) to find A , B and C .
5. Using the results of part (3) above, find
g (t ) = L
[G ( s )] .
g (t ) = L
Laplace Transform
12
Laplace Transform
13
Laplace Transform
14
Consider the following transfer function (Instructors are required to fill in the
blanks such that the poles are real and distinct)
G(s) =
num ( s )
den ( s )
G(s)
s +
s 2 ( s + )
G(s) .
in partial fractions
G ( s ) by hand calculations.
4. Check your finding by using the MATLAB residue command resi2 (See
help resi2) to find the partial fraction expansion of G ( s ) .
5. Using the results of part (3) above, find
g (t ) = L
[G ( s )] .
6. Use the symbolic MATLAB command ilaplace (See help ilaplace), find
g (t ) = L
Laplace Transform
15
Laplace Transform
16
Laplace Transform
17
Laplace Transform
18
Consider the following transfer function (Instructors are required to fill in the
blanks such that the poles are complex conjugates)
G( s) =
num ( s )
den ( s )
s +
s2 + s +
g (t ) = L
[G ( s )] .
1
L [G ( s )]
g (t ) = L
Laplace Transform
19
Laplace Transform
20
Laplace Transform
21
Consider the following transfer function (Instructors are required to fill in the
blanks). Use the Initial-Value and Final-Value-Theorems to find g ( 0 ) and
g ( ) where
num ( s )
s2 + s +
G(s) =
=
den ( s ) s 3 + s 2 + s +
Laplace Transform
22
f ( t ) ).
f (t )
y + y + y = f (t ) ,
y ( 0 ) = , y ( 0 ) =
where
t
0
Laplace Transform
23
..
Laplace Transform
24