Lecture Slides 04 Dynamic Response - Transient + Steady State
Lecture Slides 04 Dynamic Response - Transient + Steady State
u h y
u h y
(e
-
-
0, t<0
Z 1
- - )
-
1 1 1
st st
L {1(t)} = e dt = e = (pole at s = 0)
0 s 0 s
= (pole at s = j)
j s
Z 1
jt jt st
L {e }= e e dt
0
Laplace Transforms Revisited
Example: f (t) = cos t
⇢
1 jt 1 jt
L {cos t} = L e + e (Euler’s formula)
2 2
1 1
= L {ejt } + L {e jt } (linearity)
2 2
Z 1 Z 1
1 1
jt jt st
L {e } = e e dt = e(j s)t
dt = e(j s)t
0 0 j s 0
1
= (pole at s = j)
j s
Z 1 Z 1
jt jt st (j+s)t
L {e }= e e dt = e dt
0 0
Laplace Transforms Revisited
Example: f (t) = cos t
⇢
1 jt 1 jt
L {cos t} = L e + e (Euler’s formula)
2 2
1 1
= L {ejt } + L {e jt } (linearity)
2 2
Z 1 Z 1
1 1
jt jt st
L {e } = e e dt = e(j s)t
dt = e(j s)t
0 0 j s 0
1
= (pole at s = j)
j s
Z 1 Z 1
1 1
jt jt st (j+s)t (j+s)t
L {e }= e e dt = e dt = e
0 0 j+s 0
Laplace Transforms Revisited
Example: f (t) = cos t
⇢
1 jt 1 jt
L {cos t} = L e + e (Euler’s formula)
2 2
1 1
= L {ejt } + L {e jt } (linearity)
2 2
Z 1 Z 1
1 1
jt jt st
L {e } = e e dt = e(j s)t
dt = e(j s)t
0 0 j s 0
1
= (pole at s = j)
j s
Z 1 Z 1
1 1
jt jt st (j+s)t (j+s)t
L {e }= e e dt = e dt = e
0 0 j+s 0
1
=
j+s
Laplace Transforms Revisited
Example: f (t) = cos t
⇢
1 jt 1 jt
L {cos t} = L e + e (Euler’s formula)
2 2
1 1
= L {ejt } + L {e jt } (linearity)
2 2
Z 1 Z 1
1 1
jt jt st
L {e } = e e dt = e(j s)t
dt = e(j s)t
0 0 j s 0
1
= (pole at s = j)
j s
Z 1 Z 1
1 1
jt jt st (j+s)t (j+s)t
L {e }= e e dt = e dt = e
0 0 j+s 0
1
= (pole at s = j)
j+s
Laplace Transforms Revisited
Example: f (t) = cos t
⇢
1 jt 1 jt
L {cos t} = L e + e (Euler’s formula)
2 2
1 1
= L {ejt } + L {e jt } (linearity)
2 2
Z 1 Z 1
1 1
jt jt st
L {e } = e e dt = e(j s)t
dt = e(j s)t
0 0 j s 0
1
= (pole at s = j)
j s
Z 1 Z 1
1 1
jt jt st (j+s)t (j+s)t
L {e }= e e dt = e dt = e
0 0 j+s 0
1
= (pole at s = j)
j+s
— in both cases, require Re(s) > 0, i.e., s must lie in the right
half-plane (RHP)
Laplace Transforms Revisited
We already know
1
H(s) = (from earlier example)
s+1
s
U (s) = 2 (just proved)
s +1
Laplace Transforms Revisited
Convolution: L {f ? g} = L {f }L {g}
(useful because Y (s) = H(s)U (s))
We already know
1
H(s) = (from earlier example)
s+1
s
U (s) = 2 (just proved)
s +1
s
=) Y (s) = H(s)U (s) =
(s + 1)(s2 + 1)
1
y(t) = L {Y }
Laplace Transforms Revisited
Convolution: L {f ? g} = L {f }L {g}
(useful because Y (s) = H(s)U (s))
We already know
1
H(s) = (from earlier example)
s+1
s
U (s) = 2 (just proved)
s +1
s
=) Y (s) = H(s)U (s) =
(s + 1)(s2 + 1)
1
y(t) = L {Y }
s (s + 1)(as + b)
(s + 1)Y (s) = =a+
s2 + 1 (s2 + 1)
Method of Partial Fractions
Problem: compute L 1 {Y (s)}, where
s
Y (s) =
(s + 1)(s2 + 1)
s (s + 1)(as + b)
(s + 1)Y (s) = =a+
s2 + 1 (s2 + 1)
!
-
at
bj For - I
=
- tart
!
a
Magitude
W
⑳
=>
part
7
-
M(WS /
= a 5
--
)
=N =-G(w) Fant (a
=
-
M(1 )! =>
I()
= tar"(-1)
=Y
*
-
:
-
Transient and Steady-State Response
Consider the system ẏ = y+u y(0) = 0
1 t 1
u(t) = cos t ! y(t) = e + p cos(t ⇡/4)
| 2{z } | 2 {z }
transient steady-state
response response
Similarly:
L {f 00 (t)} = L {(f 0 (t))0 } = sL {f 0 (t)} f 0 (0)
= s2 F (s) sf (0) f 0 (0)
Example
Y (s)
Let’s compute the transfer function: H(s) =
U (s)
Example
Y (s)
Let’s compute the transfer function: H(s) =
U (s)
— take Laplace transform of both sides (zero I.C.’s):
Example
Y (s)
Let’s compute the transfer function: H(s) =
U (s)
— take Laplace transform of both sides (zero I.C.’s):
Y (s)
Let’s compute the transfer function: H(s) =
U (s)
— take Laplace transform of both sides (zero I.C.’s):
Y (s) 1
s2 Y (s) + 3sY (s) + 2Y (s) = U (s) H(s) = = 2
U (s) s + 3s + 2
Example (continued)
ÿ + 3ẏ + 2y = u, y(0) = ↵, ẏ(0) =
Example (continued)
ÿ + 3ẏ + 2y = u, y(0) = ↵, ẏ(0) =
Compute the step response, i.e., response to u(t) = 1(t)
Example (continued)
ÿ + 3ẏ + 2y = u, y(0) = ↵, ẏ(0) =
Compute the step response, i.e., response to u(t) = 1(t)
1
↵s + (3↵ + ) + s
Y (s) =
s2 + 3s + 2
Example (continued)
ÿ + 3ẏ + 2y = u, y(0) = ↵, ẏ(0) =
Compute the step response, i.e., response to u(t) = 1(t)
1
↵s + (3↵ + ) + s ↵s2 + (3↵ + )s + 1
Y (s) = =
s2 + 3s + 2 s(s + 1)(s + 2)
Example (continued)
ÿ + 3ẏ + 2y = u, y(0) = ↵, ẏ(0) =
Compute the step response, i.e., response to u(t) = 1(t)
1
↵s + (3↵ + ) + s ↵s2 + (3↵ + )s + 1
Y (s) = =
s2 + 3s + 2 s(s + 1)(s + 2)
1
Note: if ↵ = = 0, then Y (s) = = H(s)U (s)
s(s + 1)(s + 2)
Example (continued)
Compute the step response of
↵s2 + (3↵ + )s + 1 1
Y (s) = y(t) = L {Y (s)}
s(s + 1)(s + 2)
Example (continued)
Compute the step response of
↵s2 + (3↵ + )s + 1 1
Y (s) = y(t) = L {Y (s)}
s(s + 1)(s + 2)
Use the method of partial fractions:
↵s2 + (3↵ + )s + 1 a b c
= + +
s(s + 1)(s + 2) s s+1 s+2
Example (continued)
Compute the step response of
↵s2 + (3↵ + )s + 1 1
Y (s) = y(t) = L {Y (s)}
s(s + 1)(s + 2)
Use the method of partial fractions:
↵s2 + (3↵ + )s + 1 a b c
= + +
s(s + 1)(s + 2) s s+1 s+2
— this gives a = 1/2, b = 2↵ + 1, c = ↵ + 1/2
Example (continued)
Compute the step response of
↵s2 + (3↵ + )s + 1 1
Y (s) = y(t) = L {Y (s)}
s(s + 1)(s + 2)
Use the method of partial fractions:
↵s2 + (3↵ + )s + 1 a b c
= + +
s(s + 1)(s + 2) s s+1 s+2
— this gives a = 1/2, b = 2↵ + 1, c = ↵ + 1/2
1 1 ↵ + 1/2
Y (s) = + (2↵ + 1) +
2s s+1 s+2
Example (continued)
Compute the step response of
↵s2 + (3↵ + )s + 1 1
Y (s) = y(t) = L {Y (s)}
s(s + 1)(s + 2)
Use the method of partial fractions:
↵s2 + (3↵ + )s + 1 a b c
= + +
s(s + 1)(s + 2) s s+1 s+2
— this gives a = 1/2, b = 2↵ + 1, c = ↵ + 1/2
1 1 ↵ + 1/2
Y (s) = + (2↵ + 1) +
2s s+1 s+2
1 1 t 2t
y(t) = L {Y (s)} = 1(t) + (2↵ + 1)e + (1/2 ↵ )e
2
Example (continued)
The step response of
is given by
1 t 2t
y(t) = 1(t) + (2↵ + 1)e + (1/2 ↵ )e
2
Example (continued)
The step response of
is given by
1 t 2t
y(t) = 1(t) + (2↵ + 1)e + (1/2 ↵ )e
2
What are the transient and the steady-state terms?
Example (continued)
The step response of
is given by
1 t 2t
y(t) = 1(t) + (2↵ + 1)e + (1/2 ↵ )e
2
What are the transient and the steady-state terms?
I The transient terms are e t , e 2t
(decay to zero at exponential
rates 1 and 2)
Example (continued)
The step response of
is given by
1 t 2t
y(t) = 1(t) + (2↵ + 1)e + (1/2 ↵ )e
2
What are the transient and the steady-state terms?
I The transient terms are e t , e 2t
(decay to zero at exponential
rates 1 and 2)
1
Note the poles of H(s) = at s = 1 and s = 2
(s + 1)(s + 2)
Example (continued)
The step response of
is given by
1 t 2t
y(t) = 1(t) + (2↵ + 1)e + (1/2 ↵ )e
2
What are the transient and the steady-state terms?
I The transient terms are e t , e 2t
(decay to zero at exponential
rates 1 and 2)
1
Note the poles of H(s) = at s = 1 and s = 2
(s + 1)(s + 2)
— these are stable poles (both lie in LHP)
Example (continued)
The step response of
is given by
1 t 2t
y(t) = 1(t) + (2↵ + 1)e + (1/2 ↵ )e
2
What are the transient and the steady-state terms?
I The transient terms are e t , e 2t
(decay to zero at exponential
rates 1 and 2)
1
Note the poles of H(s) = at s = 1 and s = 2
(s + 1)(s + 2)
— these are stable poles (both lie in LHP)
I the steady-state part is 12 1(t)
Example (continued)
The step response of
is given by
1 t 2t
y(t) = 1(t) + (2↵ + 1)e + (1/2 ↵ )e
2
What are the transient and the steady-state terms?
I The transient terms are e t , e 2t
(decay to zero at exponential
rates 1 and 2)
1
Note the poles of H(s) = at s = 1 and s = 2
(s + 1)(s + 2)
— these are stable poles (both lie in LHP)
I the steady-state part is 12 1(t) — converges to steady-state value
of 1/2
DC Gain
u h y
u h y
u h y
u h y
1 H(s)
u(t) = 1(t) U (s) = =) Y (s) =
s s
Steady-State Value
u h y
1 H(s)
u(t) = 1(t) U (s) = =) Y (s) =
s s
— can we compute y(1) from Y (s)?
Steady-State Value
u h y
1 H(s)
u(t) = 1(t) U (s) = =) Y (s) =
s s
— can we compute y(1) from Y (s)?
Let’s look at some examples:
Steady-State Value
u h y
1 H(s)
u(t) = 1(t) U (s) = =) Y (s) =
s s
— can we compute y(1) from Y (s)?
Let’s look at some examples:
1
I Y (s) = , a>0 (pole at s = a < 0)
s+a
Steady-State Value
u h y
1 H(s)
u(t) = 1(t) U (s) = =) Y (s) =
s s
— can we compute y(1) from Y (s)?
Let’s look at some examples:
1
I Y (s) = , a > 0 (pole at s = a < 0)
s+a
at
y(t) = e =) y(1) = 0
Steady-State Value
u h y
1 H(s)
u(t) = 1(t) U (s) = =) Y (s) =
s s
— can we compute y(1) from Y (s)?
Let’s look at some examples:
1
I Y (s) = , a > 0 (pole at s = a < 0)
s+a
at
y(t) = e =) y(1) = 0
1
I Y (s) = , a < 0 (pole at s = a > 0)
s+a
Steady-State Value
u h y
1 H(s)
u(t) = 1(t) U (s) = =) Y (s) =
s s
— can we compute y(1) from Y (s)?
Let’s look at some examples:
1
I Y (s) = , a > 0 (pole at s = a < 0)
s+a
at
y(t) = e =) y(1) = 0
1
I Y (s) = , a < 0 (pole at s = a > 0)
s+a
at
y(t) = e =) y(1) = 1
Steady-State Value
u h y
1 H(s)
u(t) = 1(t) U (s) = =) Y (s) =
s s
— can we compute y(1) from Y (s)?
Let’s look at some examples:
1
I Y (s) = , a > 0 (pole at s = a < 0)
s+a
at
y(t) = e =) y(1) = 0
1
I Y (s) = , a < 0 (pole at s = a > 0)
s+a
at
y(t) = e =) y(1) = 1
1
I Y (s) = , ! 2 R (poles at s = ±j!, purely imaginary)
s2 + ! 2
Steady-State Value
u h y
1 H(s)
u(t) = 1(t) U (s) = =) Y (s) =
s s
— can we compute y(1) from Y (s)?
Let’s look at some examples:
1
I Y (s) = , a > 0 (pole at s = a < 0)
s+a
at
y(t) = e =) y(1) = 0
1
I Y (s) = , a < 0 (pole at s = a > 0)
s+a
at
y(t) = e =) y(1) = 1
1
I Y (s) = , ! 2 R (poles at s = ±j!, purely imaginary)
s2 + ! 2
y(t) = sin(!t) =) y(1) does not exist
Steady-State Value
u h y
1 H(s)
u(t) = 1(t) U (s) = =) Y (s) =
s s
— can we compute y(1) from Y (s)?
#
Let’s look at some examples:
1
I Y (s) = , a > 0 (pole at s = a < 0)
s+a
at
y(t) = e =) y(1) = 0
1
I Y (s) = , a < 0 (pole at s = a > 0)
s+a
at
y(t) = e =) y(1) = 1
1
I Y (s) = , ! 2 R (poles at s = ±j!, purely imaginary)
s2 + ! 2
y(t) = sin(!t) =) y(1) does not exist
c
I Y (s) = (pole at the origin, s = 0)
s
~
ooh
S P
Fr
.
.
3
-
-
-
- - /
M 2
-
I -
De
Ko
Steady-State Value
u h y
1 H(s)
u(t) = 1(t) U (s) = =) Y (s) =
s s
— can we compute y(1) from Y (s)?
Let’s look at some examples:
1
I Y (s) = , a > 0 (pole at s = a < 0)
s+a
at
y(t) = e =) y(1) = 0
1
I Y (s) = , a < 0 (pole at s = a > 0)
s+a
at
y(t) = e =) y(1) = 1
1
I Y (s) = , ! 2 R (poles at s = ±j!, purely imaginary)
s2 + ! 2
y(t) = sin(!t) =) y(1) does not exist
c
I Y (s) = (pole at the origin, s = 0)
s
↳ y(t) = c1(t) =) y(1) = c
The Final Value Theorem
We can now deduce the Final Value Theorem (FVT):
If all poles of sY (s) are strictly stable or lie in the open left
half-plane (OLHP), i.e., have Re(s) < 0, then
15
0
y(1) = lim sY (s).
s!0 <
2
y(x) -
lim Y(+)
t-1 D
The Final Value Theorem
We can now deduce the Final Value Theorem (FVT):
If all poles of sY (s) are strictly stable or lie in the open left
half-plane (OLHP), i.e., have Re(s) < 0, then
I Y (s) =
c
sY (s) = c
s
poles at infinity, so y(1) = c –I FVT gives correct answer
Back to DC Gain
u h y
Back to DC Gain
u h y
H(s)
Step response: Y (s) =
s
Back to DC Gain
u h y
H(s)
Step response: Y (s) =
s
— if all poles of sY (s) = H(s) are strictly stable, then
y(1) = lim H(s)
s!0
by the FVT.
Back to DC Gain
u h y
H(s)
Step response: Y (s) =
s
— if all poles of sY (s) = H(s) are strictly stable, then
y(1) = lim H(s)
s!0
by the FVT.
Example: compute DC gain of the system with transfer function
s2 + 5s + 3
H(s) =
s3 + 4s + 2s + 5
Back to DC Gain
u h y
H(s)
Step response: Y (s) =
s
— if all poles of sY (s) = H(s) are strictly stable, then
y(1) = lim H(s)
s!0
by the FVT.
Example: compute DC gain of the system with transfer function
s2 + 5s + 3
H(s) =
s3 + 4s + 2s + 5
All poles of H(s) are strictly stable (we will see this later using
the Routh–Hurwitz criterion), so
3
y(1) = H(s) = .
s=0 5