Root Locus
Root Locus
ROOT LOCUS
H(s)
C(s) K ⋅G(s)
=
R(s) 1 + K ⋅ G(s)⋅ H(s)
k = 0,1,2,
C2
Example:
K
K ⋅ G(s) ⋅ H(s) =
s ⋅ (s + 1)
1 + K ⋅ G(s) ⋅H(s) = 0 ⇒ s2 + s + K = 0
1 1
s 1,2 = − ± ⋅ 1−4⋅K
2 2
Angle condition:
K
∠ = − ∠s − ∠s +1=− (180−θ1 ) −θ2 =±180,
s ⋅ (s +1)
s-plane
K=1
θ2 θ1
K=0 K=0
-1
K=1/4
C3
K ⋅( s + z1 )
K ⋅ G(s) ⋅ H(s) =
(s + p1 ) ⋅(s + p2 )⋅ (s + p3 )⋅ (s + p4 )
K ⋅ B1
K ⋅ G(s) ⋅ H(s) = = 1
A1 ⋅ A2 ⋅ A3 ⋅ A4
for k = 0, 1, 2, …
C4
B( s)
K H(s) ⋅ G(s) = K
A(s)
I
m
Asymptote
γ Real
σa
±180° ⋅(2k + 1)
γ = , k = 0, 1, 2, ...
n −m
n m
∑ (− p ) − ∑ (− z )
i=1
i
i=1
i
σa=
n−m
Subtract from 180° the sum of all angles from all other
zeros and poles of the open-loop system to the complex
pole (or zero) with appropriate signs.
Consider:
H(s)
B( s)
∑b ⋅ s
i =0
i
m −i
K ⋅ H ( s ) ⋅ G ( s) = K ⋅ =K⋅ n
A( s )
∑α
i =0
i ⋅ s n −i
ϕ4
s2 ϕ1
ϕ2 s1
ϕ3
ϕ5
ϕ =180 , ϕ = 0 = ϕ , ϕ +ϕ 5 = 180⋅2
s1 { ϕ1+ϕ −ϕ +ϕ2 +ϕ = 33⋅180 4
1 2 3 4 5
Magnitude condition:
m
∏ (s + z i )
B(s) 1 i=1
= = n
A(s) K
∏ (s + pi )
i= 1
K=0 open loop poles
K=∞ m open loop zeros
m-n branches approach infinity
B(s)
K⋅ = −1
A(s)
K ⋅ B(s) K
lim = lim s n− m = −1
s →∞ A(s) s→ ∞
s n−m = − K for s → ∞
Using the angle condition:
(n − m) ⋅ ∠s = ± 180, ⋅ (2 ⋅ k + 1)
leading to
±180, ⋅ (2 ⋅ k + 1)
∠s= γ =
n− m
m n n − m −1
− ∑ zi − ∑ pi ⋅ s
n− m
s + ... = − K
i=1 i =1
sn −m = −K
m n n m
− ∑ z + ∑ p ∑− p − ∑− z
i =1
i
i =1
i
i =1
i
i =1
i
this implies σ a =− =
n−m n−m
'
⇒ for K = −
A ( s)
, f(s) has multiple roots
B ' (s )
A(s)
K =− from f (s) = 0
B(s)
This gives:
Extended Rule 6:
Consider
180,
θy = ±
y
θ3
θ2 = 90 ,
θ3 = 180, − (θ1 + θ 2 − ϕ 1 )
ϕ1 θ1
θ2
From
n
A( s) + K ⋅ B( s)= ∏ (s + r )
i=1
i
we have
n m n
∏ (s + p ) + K ⋅∏ (s + r )=∏ (s + r )
i =1
i
i =1
i
i =1
i
with
n
∏ (s + p )
m
A( s)=
i =1
i
and
B (s)= ∏
i =1
(s + zi )
the following:
Positional servo.
Closed-loop poles:
s = − 0.1 ± j ⋅ 0.995
5
Open-loop of system I: GI (s) = s ⋅ ( 5 ⋅ s + 1)
5 ⋅ (1 + 0.8s)
G(s) =
Open-loop of systems II and III: s ⋅ ( 5 ⋅ s + 1)
C18
jω
j4
X X
0 σ
-2 -1
-j1
a) System I
jω
Closed-loop zero
j1
X X
0 σ
-2 -1
-j1
b) System II
C19
jω
Open-loop zero
j1
X X
0 σ
-2 -1
-j1
c) System III
Closed-loop zeros:
System I: none
System II: 1+0.8s=0
System III: none
Observations:
R(s) C(s)
+ K G(s)
-
O unstable
X X X
X
O unstable
stable
x(t) y(t)
G(s)
X(s) Y(s)
ωX
x(t) = X sin(ω t) X(s) =
s +ω
2 2
ωX a a bi
Y (s) = G(s) ⋅ X(s) = G(s) ⋅
s +ω
2 2 = +
s + jω s − jω
+ ∑s+s
i i
+ ∑ bi e
− jω t jω t −s i t
y(t) = ae + ae
i
stable system ⇔ Re(−si ) < 0 for all i
− jω t
for t→∞ ⇒ y(t ) = ae + a e jω t
ωX X G(− j ω)
a = G(s)⋅ 2 ⋅( s+ j ω ) = −
s +ω
2
( s= − j ω )
2j
ωX X G( j ω )
a = G(s)⋅ 2 2 (s− j ω) =
s +ω (s = j ω )
2j
C23
Im(G( j ω ))
G( jω ) = G ( j ω ) ⋅ e ϕ
j
ϕ = tan −1
Re( G( j ω ))
and
−jϕ
G(− j ω ) = G ( j ω ) ⋅ e
ω +ϕ − ω +ϕ
e j( t ) − e j( t )
y(t) = X ⋅ G(jω) ⋅ = Y ⋅ sin(ωt + ϕ)
2j
Y( j ω )
G( j ω ) =
X( jω )
Y (j ω )
G( j ω ) =
X( j ω ) Magnitude response
Y( j ω )
ϕ = ∠ (G(j ω ) ) = ∠
X( j ω ) Phase response
C24
K(s + z)
G(s) =
s(s + p)
K ⋅ jω + z
G( j ω ) =
j ω ⋅ jω + p
∠ G( j ω ) = ϕ − θ1 − θ 2
• Bode Diagrams
• Polar Plots (Nyquist Plots)
• Log-Magnitude-Versus-Phase Plots (Nichols Plots)
C25
Bode Diagrams
• Magnitude response G( j ω )
20 log G( j ω ) in dB
• Gain K
±1
• Integral or derivative factors ( jω)
±1
• First-order factors (1 + jωT )
±1
2
jω jω
1 + 2 j ω + ω
• Quadratic factors n n
1. Gain Factor K
Phase is zero
C26
±1
2. Integral or derivative factors ( jω)
• ( jω )−
1
1
20 log = − 20 log ω
jω
magnitude: straight line with slope –20 dB/decade
phase: -90o
• ( jω )
20 log j ω = 20log ω
magnitude: straight line with slope 20 dB/decade
phase: +90o
C27
±1
3. First order factors (1 + jω T)
−1
• (1 + jω T)
Magnitude:
1
20 log = − 20 log 1 + ω 2T 2 dB
1 + j ωT
−1
for ω << T ⇒ 0 dB magnitude
⇒ − 20 log(ω T ) dB magnitude
−1
for ω >> T
Phase:
,
for ω = 0 ϕ = 0
1 −1 T
for ω = ⇒ − tan =1 ϕ = − 45,
T T
,
for ω = ∞ ϕ = − 90
C28
+1
• (1 + jω T)
1
Using 20 log 1+ jω T = − 20log
1 + jω T
1
∠ (1+ jω T ) = tan (ω T ) = − ∠
−1
1 + jω T
C29
4. Quadratic Factors
1
G ( jω ) = 2 0 <ς <1
ω jω
1 + 2ς j +
ω
n n ω
Magnitude:
2
ω 2 ω
2
20 log G( j ω ) = − 20 log 1 − 2 − 2 j
ω n ω n
for
ω << ω n ⇒ 0 dB
for
ω2 ω
ω >> ω n ⇒ − 20 log 2 = − 40log dB
ω n ωn
Phase:
2ς ω
ωn
ϕ = tan −1 ∠G( jω ) = − tan −1 2
ω
1 −
ω n
Resonant Frequency:
ωT =ωn 1 − 2ς 2
Example:
10( jω + 3)
G( jω ) =
(
( jω ) ⋅ ( jω + 2 )⋅ ( j ω ) + jω + 2
2
)
( jω )
7.5 ⋅ + 1
3
G ( jω ) =
jω ( jω ) jω
2
jω + 1 + + 1
2 2 2
C31
C32
Minimum phase systems have all poles and zeros in the left half s-plane
and were discussed before
Consider
Consider
1 1 1
G1 (s) = G2 (s) = G3 (s) =
1 + Ts 1 − Ts Ts − 1
then
∠ G2 ( j ω ) = − ∠ G1 ( jω )
∠ G3 ( j ω ) = − 180 o − ∠ G1 ( jω )
Type 0
lim G( jω) = K p
ω→ 0
C35
Type 1
Kv = lim jω G( jω )
ω→ 0
Kv
G( jω ) = ω <<
jω for
Kv
20log K v = 20log for ω = 1
jω
Type 2
Acceleration Error Coefficient Ka
K a = lim ( jω )2 G ( jω) ≠ 0
ω →0
( K p = Kv = ∞ )
Ka
G( jω ) = ω << 1
( jω ) 2 for
Ka
20log K a = 20log for ω = 1
( jω ) 2
G( jω ) = G( jω ) ∗ ∠ G( jω )
= Re[ G( jω )] + Im[ G( jω )]
1 1 1
G ( jω ) = = −j = ∠ − 90,
jω ω ω
,
G( jω ) = jω = ω ∠ 90
C38
±1
First order factors (1 + jω T)
1
G( jω ) = = X + jY
1 + jω T
1 −ω T
X = Y =
1+ ω 2 T 2 ,
1 +ω 2 T2
2 2 2
It can be show that (X − 0.5) + Y = (0.5)
⇒ Polar plot is a circle with
• Center (1/2, 0) and
• Radius 0.5.
C39
Quadratic Factors
1
G( jω ) =
2
ω jω 1 >ζ > 0
1 2
+ ζ j +
ωn ωn
lim G( jω ) = 1 ∠ 0
,
ω→0
+1
ω ω
2
1 + 2ς j
+
j
ωn ωn
lim G( jω ) = 1 ∠0,
ω →0
,
lim G( jω ) = ∞ ∠ 180
ω→ ∞
C41
Re [G( jω )] for ω → 0
Example:
• Bode Diagram
• Polar Plot
• Log-Magnitude-Versus-Phase Plots
C45
R(s) + C(s)
G(s)
-
H(s)
Closed-loop: G(s)
1+ G(s)⋅ H(s )
Preview
Mathematical Background
• Mapping theorem
• Nyquist path
Z=N+P
Mapping Theorem
A(s)
F(s) =
Q(s)
VSODQH ) V SODQH
FRQWRXU
)V
[
R
[
R
[
R
[
]HURV $ V
SROHV 4 V
C48
VSODQH ) V SODQH
M ω
5 ∞
M ω V + V
R=∞ 1+G(jω)H(jω)
-1 0 Re
0 1
G(jω)H(j ω)
F(s)=1+G(s)H(s)
-jω
Consider
R(s) + C(s)
G(s)
-
H(s)
Z=N+P
(...)
G(s)⋅ H (s) =
λ
s (...)
Gain margin in dB: Kg in dB= −20logG jω
1
R(s) + ωn C(s)
- s (s + 2 ζω n )
C (s ) ωn 2 C( jω )
=
R (s ) s 2 + 2ζω n s + ωn 2 = M( ω ) ⋅ e ( ω )
jα
R( jω )
Phase margin:
ωn 2
G ( jω ) =
jω ( jω + 2ζω n ) becomes unity for ω1 = ω n 1 + 4ζ 4 − 2ζ 2
γ = tan −1
2ζω n = tan −1
2ζ
and ω1
1 + ζ 4 − 2ζ 2
Æ γ depends only on ζ
C59
Performance specifications in the frequency
domain:
0dB Mr
-3dB
ωr ωc
ωc:
ωc : Cutoff Frequency
0 ≤ ω ≤ ωc : Bandwidth
R(s) + C(s)
G(s)
-
→
G(s) OA
=
1+ G(s) →
PA
C61
Closed-Loop Frequency Response:
C(jω )
= M ⋅e jα
G( jω )
=
R(jω ) 1+ G( jω )
G( jω ) = X + jY
X + jY
M= = const
1 + X + jY
2
M2 M2
X + + Y2 =
M2 −1 2 2
M −1
G( jω ) = X + jY Æ ∠e jα = ∠
X + jY
1+ X + jY
= const
2 2 2
X + 1 + Y + 1 = 1 + 1 , N = tan α
2 2N 4 2N
C62
Figure:
(a) G(jω) locus superimposed on a family of M circles;
(b) G(jω) locus superimposed om a family of N circles;
(c) Closed-loop frequency-response curves
C64
Experimental Determination of Transfer Function
Æ
• Minimum phase: ω ∞ phase = -90 (n - m)
(n-m) difference in the order of denominator and
numerator.