Lecture Slides - 05 System Modeling Using Block Diagrams &. 2nd Order Systems
Lecture Slides - 05 System Modeling Using Block Diagrams &. 2nd Order Systems
decompose
large system ) * smaller blocks (subsystems)
compose
System Modeling Diagrams
decompose
large system ) * smaller blocks (subsystems)
compose
decompose
large system ) * smaller blocks (subsystems)
compose
We will take smaller blocks from some given library and play
with them to create/build more complicated systems.
All-Integrator Diagrams
Two warnings:
I We can (and will) work either with u, y (time domain) or
with U, Y (s-domain) — will often go back and forth
I When working with block diagrams, we typically ignore
initial conditions.
All-Integrator Diagrams
Two warnings:
I We can (and will) work either with u, y (time domain) or
with U, Y (s-domain) — will often go back and forth
I When working with block diagrams, we typically ignore
initial conditions.
ÿ = u () s2 Y = U
-
Example 1
ÿ = u () s2 Y = U
This is obvious:
ẏ
u 1/s 1/s y
or
sY
U 1/s 1/s Y
Example 2
(building on Example 1)
ÿ + a1 ẏ + a0 y = u () s2 Y + a1 sY + a0 Y = U
U (s)
or Y (s) = 2
s + a1 s + a0
Example 2
(building on Example 1)
ÿ + a1 ẏ + a0 y = u () s2 Y + a1 sY + a0 Y = U
U (s)
or Y (s) = 2
s + a1 s + a0
Always solve for the highest derivative:
ÿ = a1 ẏ a y+u
| {z 0 }
=v
Example 2
(building on Example 1)
ÿ + a1 ẏ + a0 y = u () s2 Y + a1 sY + a0 Y = U
U (s)
or Y (s) = 2
s + a1 s + a0
Always solve for the highest derivative:
ÿ = a1 ẏ a y+u
| {z 0 }
=v
+ v ẏ
u 1/s 1/s y
a1
a0
Example 2
(building on Example 1)
ÿ + a1 ẏ + a0 y = u () s2 Y + a1 sY + a0 Y = U
U (s)
or Y (s) = 2
s + a1 s + a0
Always solve for the highest derivative:
ÿ = a1 ẏ a y+u
| {z 0 }
=v
+ V 1/s sY
U 1/s Y
a1
a0
Example 3
1 X
U s2 + a1 s + a0
b1 s + b 0 Y
1 X
U s2 + a1 s + a0
b1 s + b 0 Y
1 X
U s2 + a1 s + a0
b1 s + b 0 Y
Example 3, continued
1
Step 1: decompose H(s) = · (b1 s + b0 )
s2 + a1 s + a0
1 X
U s2 + a1 s + a0
b1 s + b 0 Y
1 X
U s2 + a1 s + a0
b1 s + b 0 Y
2
+s X sX
U 1/s 1/s X
a1
a0
Example 3, continued
Step 3: now we notice that
b1
2 +
+s X sX X +
U 1/s 1/s b0 Y
a1
a0
Example 3, continued
b1 s + b0
All-integrator diagram for H(s) =
s2 + a1 s + a0
b1
2 +
+s X sX X +
U 1/s 1/s b0 Y
a1
a0
Example 3, continued
b1 s + b0
All-integrator diagram for H(s) =
s2 + a1 s + a0
b1
2 +
+s X sX X +
U 1/s 1/s b0 Y
a1
a0
b1
2 +
+s X sX X +
U 1/s 1/s b0 Y
a1
a0
State-space model:
Example 3, continued
b1
2 +
+s X sX X +
U 1/s 1/s b0 Y
a1
a0
State-space model:
s2 X = U a1 sX a0 X
Example 3, continued
b1
2 +
+s X sX X +
U 1/s 1/s b0 Y
a1
a0
State-space model:
s2 X = U a1 sX a0 X
ẍ = a1 ẋ a0 x + u
Example 3, continued
b1
2 +
+s X sX X +
U 1/s 1/s b0 Y
a1
a0
State-space model:
s2 X = U a1 sX a0 X Y = b1 sX + b0 X
ẍ = a1 ẋ a0 x + u
Example 3, continued
b1
2 +
+s X sX X +
U 1/s 1/s b0 Y
a1
a0
State-space model:
s2 X = U a1 sX a0 X Y = b1 sX + b0 X
ẍ = a1 ẋ a0 x + u y = b1 ẋ + b0 x
Example 3, continued
State-space model:
ẍ = a1 ẋ a0 x + u y = b1 ẋ + b0 x
Example 3, continued
State-space model:
ẍ = a1 ẋ a0 x + u y = b1 ẋ + b0 x
x1 = x, x2 = ẋ
Example 3, continued
State-space model:
ẍ = a1 ẋ a0 x + u y = b1 ẋ + b0 x
x1 = x, x2 = ẋ
✓ ◆ ✓ ◆✓ ◆ ✓ ◆
ẋ1 0 1 x1 0
= + u
ẋ2 a0 a1 x2 1
Example 3, continued
State-space model:
ẍ = a1 ẋ a0 x + u y = b1 ẋ + b0 x
x1 = x, x2 = ẋ
✓ ◆ ✓ ◆✓ ◆ ✓ ◆ ✓ ◆
ẋ1 0 1 x1 0 x1
= + u y = b0 b1
ẋ2 a0 a1 x2 1 x2
Example 3, continued
State-space model:
ẍ = a1 ẋ a0 x + u y = b1 ẋ + b0 x
x1 = x, x2 = ẋ
✓ ◆ ✓ ◆✓ ◆ ✓ ◆ ✓ ◆
ẋ1 0 1 x1 0 x1
= + u y = b0 b1
ẋ2 a0 a1 x2 1 x2
State-space model:
ẍ = a1 ẋ a0 x + u y = b1 ẋ + b0 x
x1 = x, x2 = ẋ
✓ ◆ ✓ ◆✓ ◆ ✓ ◆ ✓ ◆
ẋ1
ẋ2
=
0
a0
- -
1
a1
x1
x2
+
0
1
u
-
0
y = b0 b1
x1
x2 -
I
Easily generalizes to dimension > 1 Ines
The reason behind the name will be made clear later in the
semester
Example 3, wrap-up
b1 s + b0
All-integrator diagram for H(s) =
s 2 + a1 s + a0
b1
2 +
+s X sX X +
U 1/s 1/s b0 Y
a1
a0
State-space model:
✓ ◆ ✓ ◆✓ ◆ ✓ ◆ ✓ ◆
ẋ1 0 1 x1 0 x1
= + u y = b0 b1
ẋ2 a0 a1 x2 1 x2
Example 3, wrap-up
b1 s + b0
All-integrator diagram for H(s) =
s 2 + a1 s + a0
b1
2 +
+s X sX X +
U 1/s 1/s b0 Y
a1
a0
State-space model:
✓ ◆ ✓ ◆✓ ◆ ✓ ◆ ✓ ◆
ẋ1 0 1 x1 0 x1
= + u y = b0 b1
ẋ2 a0 a1 x2 1 x2
U G 1 G2 Y
Y
= G1 G2 (for SISO systems, the order of G1 and G2
U
does not matter)
Basic System Interconnections: Series & Parallel
Series connection
(G is common
U G1 G2 Y
notation for t.f.’s)
U G 1 G2 Y
Y
= G1 G2 (for SISO systems, the order of G1 and G2
U
does not matter)
Parallel connection
Basic System Interconnections: Series & Parallel
Series connection
(G is common
U G1 G2 Y
notation for t.f.’s)
U G 1 G2 Y
Y
= G1 G2 (for SISO systems, the order of G1 and G2
U
does not matter)
Parallel connection
G1
+
U Y
+
G2
Basic System Interconnections: Series & Parallel
Series connection
(G is common
U G1 G2 Y
notation for t.f.’s)
U G 1 G2 Y
Y
= G1 G2 (for SISO systems, the order of G1 and G2
U
does not matter)
Parallel connection
G1
+
U Y
+
G2
Y U G1 + G2 Y
= G1 + G2
U
Basic System Interconnections: Negative Feedback
+ U
R G1 Y
Find the transfer function from R
W
(reference) to Y
G2
Basic System Interconnections: Negative Feedback
+ U
R G1 Y
Find the transfer function from R
W
(reference) to Y
G2
U =R W
Basic System Interconnections: Negative Feedback
+ U
R G1 Y
Find the transfer function from R
W
(reference) to Y
G2
U =R W
Y = G1 U
Basic System Interconnections: Negative Feedback
+ U
R G1 Y
Find the transfer function from R
W
(reference) to Y
G2
U =R W
Y = G1 U
= G1 (R W)
Basic System Interconnections: Negative Feedback
+ U
R G1 Y
Find the transfer function from R
W
(reference) to Y
G2
U =R W
Y = G1 U
= G1 (R W)
= G1 R G1 G2 Y
Basic System Interconnections: Negative Feedback
+ U
R G1 Y
Find the transfer function from R
W
(reference) to Y
G2
G1
=) Y = R
1 + G1 G2
U =R W
Y = G1 U
= G1 (R W) G1
= G1 R G1 G2 Y U Y
1 + G1 G2
Basic System Interconnections: Negative Feedback
+ U
R G1 Y G1
W =) Y = R
1 + G1 G2
G2
Basic System Interconnections: Negative Feedback
+ U
R G1 Y G1
W =) Y = R
1 + G1 G2
G2
forward gain
1 + loop gain
Basic System Interconnections: Negative Feedback
+ U
R G1 Y G1
W =) Y = R
1 + G1 G2
G2
forward gain
1 + loop gain
forward gain
Let’s practice with deriving transfer functions:
1 + loop gain
Unity Feedback
+ E U
R G2 G1 Y
forward gain
Let’s practice with deriving transfer functions:
1 + loop gain
I Reference R to output Y :
Unity Feedback
+ E U
R G2 G1 Y
forward gain
Let’s practice with deriving transfer functions:
1 + loop gain
I Reference R to output Y :
Y G1 G 2
=
R 1 + G1 G 2
Unity Feedback
+ E U
R G2 G1 Y
forward gain
Let’s practice with deriving transfer functions:
1 + loop gain
I Reference R to output Y :
Y G1 G 2
=
R 1 + G1 G 2
forward gain
Let’s practice with deriving transfer functions:
1 + loop gain
I Reference R to output Y :
Y G1 G 2
=
R 1 + G1 G 2
forward gain
Let’s practice with deriving transfer functions:
1 + loop gain
I Reference R to output Y :
Y G1 G 2
=
R 1 + G1 G 2
I Error E to output Y :
Unity Feedback
+ E U
R G2 G1 Y
forward gain
Let’s practice with deriving transfer functions:
1 + loop gain
I Reference R to output Y :
Y G1 G 2
=
R 1 + G1 G 2
I Error E to output Y :
Y
= G1 G 2 (no feedback path)
E
Block Diagram Reduction
examples. ↳ Norman S .
Nise
General strategy:
I Name all the variables in the diagram
I Write down as many relationships between these variables
as you can
I Learn to recognize series, parallel, and feedback
interconnections
I Replace them by their equivalents
I Repeat
Prototype 2nd-Order System
For now, we will only look at second-order systems, but this will
be sufficient to develop some nontrivial intuition (dominant
poles).
Prototype 2nd-Order System
For now, we will only look at second-order systems, but this will
be sufficient to develop some nontrivial intuition (dominant
poles).
!n2
H(s) =
s2 + 2⇣!n s + !n2
Prototype 2nd-Order System
!n2
H(s) =
s2 + 2⇣!n s + !n2
Comments:
Prototype 2nd-Order System
!n2
H(s) =
s2 + 2⇣!n s + !n2
Comments:
I ⇣ > 0, !n > 0 are arbitrary parameters
Prototype 2nd-Order System
!n2
H(s) =
s2 + 2⇣!n s + !n2
Comments:
I ⇣ > 0, !n > 0 are arbitrary parameters
I the denominator is a general 2nd-degree monic polynomial,
just written in a weird way ↓ variable
polynomial
in
single
ehighest
active
Prototype 2nd-Order System
!n2
H(s) =
s2 + 2⇣!n s + !n2
Comments:
I ⇣ > 0, !n > 0 are arbitrary parameters
I the denominator is a general 2nd-degree monic polynomial,
just written in a weird way
I H(s) is normalized to have DC gain = 1 (provided DC gain
exists)
Prototype 2nd-Order System
!n2
H(s) =
s2 + 2⇣!n s + !n2
Prototype 2nd-Order System
!n2
H(s) =
s2 + 2⇣!n s + !n2
By the quadratic formula, the poles are:
p
s= ⇣!n ± !n ⇣ 2 1
⇣ p ⌘
= !n ⇣ ± ⇣ 2 1
Prototype 2nd-Order System
!n2
H(s) =
s2 + 2⇣!n s + !n2
By the quadratic formula, the poles are:
p
s= ⇣!n ± !n ⇣ 2 1
⇣ p ⌘
= !n ⇣ ± ⇣ 2 1
p
s= ⇣!n ± !n ⇣ 2 1
⇣ p ⌘
= !n ⇣ ± ⇣ 2 1
p
s= ⇣!n ± !n ⇣ 2 1
⇣ p ⌘
= !n ⇣ ± ⇣ 2 1
p
s= ⇣!n ± !n ⇣ 2 1
⇣ p ⌘
= !n ⇣ ± ⇣ 2 1
p
s= ⇣!n ± !n ⇣ 2 1
⇣ p ⌘
= !n ⇣ ± ⇣ 2 1
p
s= ⇣!n ± !n ⇣ 2 1
⇣ p ⌘
= !n ⇣ ± ⇣ 2 1
p
s= ⇣!n ± !n ⇣ 2 1
⇣ p ⌘
= !n ⇣ ± ⇣ 2 1
p
s= ⇣!n ± !n ⇣ 2 1
⇣ p ⌘
= !n ⇣ ± ⇣ 2 1
s= ± j!d
p
where = ⇣!n , !d = !n 1 ⇣2
Prototype 2nd-Order System
!n2
H(s) = , ⇣<1
s2 + 2⇣!n s + !n2
Prototype 2nd-Order System
!n2
H(s) = , ⇣<1
s2 + 2⇣!n s + !n2
The poles are
p
s= ⇣!n ± j!n 1 ⇣2 = ± j!d
Prototype 2nd-Order System
!n2
H(s) = , ⇣<1
s2 + 2⇣!n s + !n2
The poles are
p
s= ⇣!n ± j!n 1 ⇣2 = ± j!d
Im
p
!d = !n 1 ⇣2 Note that
!n
2
' + !d2 = ⇣ 2 !n2 + !n2 ⇣ 2 !n2
= ⇣!n Re
0 = !n2
⇣!n
cos ' = =⇣
!n
2nd-Order Response
Let’s compute the system’s impulse and step response:
2nd-Order Response
Let’s compute the system’s impulse and step response:
!n2 !n2
H(s) = =
s2 + 2⇣!n s + !n2 (s + )2 + !d2
2nd-Order Response
Let’s compute the system’s impulse and step response:
!n2 !n2
H(s) = =
s2 + 2⇣!n s + !n2 (s + )2 + !d2
I Impulse response:
⇢
1 1 !n2
h(t) = L {H(s)} = L
(s + )2 + !d2
2nd-Order Response
Let’s compute the system’s impulse and step response:
!n2 !n2
H(s) = =
s2 + 2⇣!n s + !n2 (s + )2 + !d2
I Impulse response:
⇢
1 1 (!n2 /!d )!d
h(t) = L {H(s)} = L
(s + )2 + !d2
2nd-Order Response
Let’s compute the system’s impulse and step response:
!n2 !n2
H(s) = =
s2 + 2⇣!n s + !n2 (s + )2 + !d2
I Impulse response:
⇢
1 1 (!n2 /!d )!d
h(t) = L {H(s)} = L
(s + )2 + !d2
!n2 t
= e sin(!d t)
!d
2nd-Order Response
Let’s compute the system’s impulse and step response:
!n2 !n2
H(s) = =
s2 + 2⇣!n s + !n2 (s + )2 + !d2
I Impulse response:
⇢
1 1 (!n2 /!d )!d
h(t) = L {H(s)} = L
(s + )2 + !d2
!n2 t
= e sin(!d t) (table, # 20)
!d
2nd-Order Response
Let’s compute the system’s impulse and step response:
!n2 !n2
H(s) = =
s2 + 2⇣!n s + !n2 (s + )2 + !d2
I Impulse response:
⇢
1 1 (!n2 /!d )!d
h(t) = L {H(s)} = L
(s + )2 + !d2
!n2 t
= e sin(!d t) (table, # 20)
!d
↳
Laplace table)
I Step response:
⇢ ⇢
1 H(s) 1 !n2
L =L
s s[(s + )2 + !d2 ]
2nd-Order Response
Let’s compute the system’s impulse and step response:
!n2 !n2
H(s) = =
s2 + 2⇣!n s + !n2 (s + )2 + !d2
I Impulse response:
⇢
1 1 (!n2 /!d )!d
h(t) = L {H(s)} = L
(s + )2 + !d2
!n2 t
= e sin(!d t) (table, # 20)
!d
↳ Laplaces
I Step response:
⇢ ⇢ 2
1 H(s) 1 + !d2
L =L
s s[(s + )2 + !d2 ]
2nd-Order Response
Let’s compute the system’s impulse and step response:
!n2 !n2
H(s) = =
s2 + 2⇣!n s + !n2 (s + )2 + !d2
I Impulse response:
⇢
1 1 (!n2 /!d )!d
h(t) = L {H(s)} = L
(s + )2 + !d2
!n2 t
= e sin(!d t) (table, # 20)
!d
↳ Laplace)
I Step response:
⇢ ⇢ 2
H(s) + !d2
L 1 =L 1
s s[(s + )2 + !d2 ]
✓ ◆
t
=1 e cos(!d t) + sin(!d t)
!d
2nd-Order Response
Let’s compute the system’s impulse and step response:
!n2 !n2
H(s) = =
s2 + 2⇣!n s + !n2 (s + )2 + !d2
I Impulse response:
⇢
1 1 (!n2 /!d )!d
h(t) = L {H(s)} = L
(s + )2 + !d2
!n2 t
= e sin(!d t) (table, # 20)
!d
I Step response:
⇢ ⇢ 2
H(s) + !d2
L 1 =L 1
s s[(s + )2 + !d2 ]
✓ ◆
t
=1 e cos(!d t) + sin(!d t) (table, #21)
!d
2nd-Order Step Response
!n2 !n2
H(s) = =
s2 + 2⇣!n s + !n2 (s + )2 + !d2
✓ ◆
t
u(t) = 1(t) ! y(t) = 1 e cos(!d t) + sin(!d t)
!d
p
where = ⇣!n and !d = !n 1 ⇣ 2 (damped frequency)
2nd-Order Step Response
!n2 !n2
H(s) = =
s2 + 2⇣!n s + !n2 (s + )2 + !d2
✓ ◆
t
u(t) = 1(t) ! y(t) = 1 e cos(!d t) + sin(!d t)
!d
p
where = ⇣!n and !d = !n 1 ⇣ 2 (damped frequency)
yHtL
The parameter ⇣ is called
1.5 the damping ratio
I ⇣ > 1: system is
1.0 overdamped
I ⇣ < 1: system is
z=0.1
0.5 underdamped
z=0.9
I ⇣ = 0: no damping
z=1
t (!d = !n )
2 4 6 8 10 12 14
2nd-Order Step Response
!n2 !n2
H(s) = =
s2 + 2⇣!n s + !n2 (s + )2 + !d2
✓ ◆
t
u(t) = 1(t) ! y(t) = 1 e cos(!d t) + sin(!d t)
!d
p
where = ⇣!n and !d = !n 1 ⇣ 2 (damped frequency)
yHtL
The parameter ⇣ is called
1.5 the damping ratio
I ⇣ > 1: system is
1.0 overdamped
I ⇣ < 1: system is
z=0.1
0.5 underdamped
z=0.9
I ⇣ = 0: no damping
z=1
t (!d = !n )
2 4 6 8 10 12 14