Lecture04 A First Analysis of Feedback
Lecture04 A First Analysis of Feedback
1
Open-loop Control System
𝑌𝑜𝑙 = 𝐺𝐷𝑜𝑙 𝑅 + 𝐺𝑊
𝐸𝑜𝑙 = 𝑅 − 𝑌𝑜𝑙
= 𝑅 − 𝐺𝐷𝑜𝑙 𝑅 + 𝐺𝑊
= 1 − 𝐺𝐷𝑜𝑙 𝑅 − 𝐺𝑊
= 1 − 𝑇𝑜𝑙 𝑅 − 𝐺𝑊
open-loop transfer function
2
Feedback Control System
𝑇𝑐𝑙 : closed-loop
transfer function
𝐺𝐷𝑐𝑙 𝐺 𝐺𝐷𝑐𝑙
𝐸𝑐𝑙 = 𝑅 − 𝑌𝑐𝑙 = 𝑅 − 𝑅+ 𝑊− 𝑉
1 + 𝐺𝐷𝑐𝑙 1 + 𝐺𝐷𝑐𝑙 1 + 𝐺𝐷𝑐𝑙
1 𝐺 𝐺𝐷𝑐𝑙
= 𝑅− 𝑊+ 𝑉
1 + 𝐺𝐷𝑐𝑙 1 + 𝐺𝐷𝑐𝑙 1 + 𝐺𝐷𝑐𝑙
3
Stability
4
Tracking
• The tracking problem is to cause the output to follow the reference input as
closely as possible.
• For the open-loop system, if the plant is stable and has neither poles nor zeros
in the RHP, the controller can be selected to cancel the transfer function of the
plant and substitute the desired transfer function.
– The controller function must be proper, i.e. no more zeros than poles.
– Cannot request an unrealistically fast design which might cause large control efforts
to saturate the system.
– Non-perfect pole-zero cancelation might cause stability issues.
• Exercise (for unity feedback system)
1 𝑐2 𝑠 2 + 𝑐1 𝑠 + 𝑐0
𝐺 𝑠 = 2 , 𝐷𝑐𝑙 𝑠 =
𝑠 + 3𝑠 + 9 𝑠 𝑠 + 𝑑1
Desired characteristic equation: 𝑠 + 6 𝑠 + 3 𝑠 2 + 3𝑠 + 9 = 0
5
Regulation
• The problem of regulation is to keep the error small when the reference is at
most a constant set point and disturbances are present.
• For the open-loop system, the controller has no influence at all on the system
response to either of the disturbance. Thus useless for regulation.
• For the closed-loop system, we find a conflict between 𝑤 and 𝑣 in the search for
a good controller.
1 𝐺 𝐺𝐷𝑐𝑙
𝐸𝑐𝑙 = 𝑅− 𝑊+ 𝑉
1 + 𝐺𝐷𝑐𝑙 1 + 𝐺𝐷𝑐𝑙 1 + 𝐺𝐷𝑐𝑙
The resolution of the dilemma is to observe that each of these terms is a
function of frequency so one of them can be large for some frequencies and
small for others.
• Exercise for a constant bias 𝑤
If 𝐷𝑐𝑙 has a pole at 𝑠 = 0, then the error due to this bias will be zero.
If 𝐺 has a pole at 𝑠 = 0, it does not help with a disturbance bias.
6
Sensitivity
Open-loop system
Closed-loop system
𝐹 can be designed to
improve the tracking
accuracy.
𝑇 1 𝐺𝐷𝑐𝑙 𝐻
𝑇
𝑆𝐹 𝑐𝑙 = 1, 𝑆𝐺 𝑐𝑙 = , 𝑆𝐻𝑇𝑐𝑙 =
1 + 𝐺𝐷𝑐𝑙 𝐻 1 + 𝐺𝐷𝑐𝑙 𝐻
8
Control of Steady-State Error
• Previously
1 𝐺 𝐷𝑐𝑙 𝐺
𝐸𝑐𝑙 = 𝑅− 𝑊+ 𝑉
1 + 𝐷𝑐𝑙 𝐺 1 + 𝐷𝑐𝑙 𝐺 1 + 𝐷𝑐𝑙 𝐺
1 𝐷𝑐𝑙 𝐺
Sensitivity: 𝑆 = Complementary Sensitivity: 𝑇 = 1 − 𝑆 =
1 + 𝐷𝑐𝑙 𝐺 1 + 𝐷𝑐𝑙 𝐺
1
𝐸𝑐𝑙 = 𝑆𝑅 − 𝑆𝐺𝑊 + 𝑇𝑉 𝐸𝑐𝑙 = 𝑅 If 𝑊 = 𝑉 = 0 and 𝐿 = 𝐷𝑐𝑙 𝐺
1+𝐿
loop transfer function
• Consider polynomial inputs, 𝑅 = 1 𝑠 𝑘+1
– 𝑘 = 0, step input or position input
– 𝑘 = 1, ramp input or velocity input
– 𝑘 = 2, parabola input or acceleration input
1 1 1 1 𝑠𝑛 1
𝑒𝑠𝑠 = lim𝑠𝐸 𝑠 = lim𝑠 = lim = lim
𝑠→0 𝑠→0 1 + 𝐿 𝑠 𝑠 𝑘+1 𝑠→0 𝐿𝑜 𝑠 𝑠 𝑘 𝑠→0 𝑠 𝑛 + 𝐾𝑛 𝑠 𝑘
1 + 𝑠𝑛
System type: 𝑛 = 0 type 0
𝑛 = 1 type 1
𝑛 = 2 type 2
9
Errors as a Function of System Type
𝐾𝑝 = lim𝐿 𝑠 , 𝑛 = 0
𝑠→0
𝐾𝑣 = lim𝑠𝐿 𝑠 , 𝑛 = 1
𝑠→0
𝐾𝑎 = lim𝑠 2 𝐿 𝑠 , 𝑛 = 2
𝑠→0
10
System Type Examples
11
System Type Examples
1 + 𝐻 − 1 𝐷𝐺 1 1 𝑠 𝜏𝑠 + 1 + 1 + 𝑘𝑡 𝑠 − 1 𝑘𝑝
𝑒𝑠𝑠 = lim𝑠 = lim
𝑠→0 1 + 𝐻𝐷𝐺 𝑠 𝑘+1 𝑠→0 𝑠 𝑘 𝑠 𝜏𝑠 + 1 + 1 + 𝑘𝑡 𝑠 𝑘𝑝
= 0, 𝑘=0
1 + 𝑘𝑡 𝑘𝑝 𝑘𝑝
= , 𝑘=1 Type 1, 𝐾𝑣 =
𝑘𝑝 1 + 𝑘𝑡 𝑘𝑝
12
System Type with Respect to Disturbances
• The transfer function from the disturbance input 𝑊(𝑠) to the error 𝐸(𝑠) is
formulated as
𝐸 𝑠 −𝑌 𝑠
= = 𝑇𝑤 𝑠
𝑊 𝑠 𝑊 𝑠
If the disturbance-to-error transfer function can be written as
1 1
𝑇𝑤 𝑠 = 𝑠 𝑛 𝑇𝑜,𝑤 𝑠 , 𝑇𝑜,𝑤 0 = , 𝑊 𝑠 = 𝑘+1
𝐾𝑛,𝑤 𝑠
1 𝑠𝑛
𝑒𝑠𝑠 = lim𝑠𝑇𝑤 𝑠 = lim𝑇𝑜,𝑤 𝑠 𝑘
𝑠→0 𝑠 𝑘+1 𝑠→0 𝑠
13
System Type Example w.r.t. Disturbance
𝐴𝑘𝐼
Type 1, 𝐾𝑛,𝑤 =−
𝐵
14
PID Control
15
Proportional Control (P)
𝑠 2 + 2𝜁𝜔𝑛 𝑠 + 𝜔𝑛2 = 0
• The designer can control the constant term in this equation, which determines
the natural frequency, but cannot control the damping of the equation.
• The system is type 0 and if 𝑘𝑃 is made large to get adequately small steady-
state error. However, the damping may be too low for satisfactory transient
response with p control alone.
16
Integral Control (I)
17
History of System Error
18
Derivative Control (D)
19
Proportional Plus Integral Control (PI)
𝐴𝑘𝐼 𝐴𝑘 + 1
𝜔𝑛 = , 𝜁= 𝑃
𝜏 2𝜏𝜔𝑛
• If the plant is second order, the characteristic equation with PI control is
𝑘𝐼 𝐴
1 + 𝑘𝑃 + =0
𝑠 𝑠 2 + 𝑎1 𝑠 + 𝑎2
𝑠 3 + 𝑎1 𝑠 2 + (𝑎2 + 𝐴𝑘𝑃 )𝑠 + 𝐴𝐾𝐼 = 0
PI control can only be used to set two coefficients. 20
PID Control
𝑘𝐼 𝐴
1 + 𝑘𝑃 + + 𝑘𝐷 𝑠 2 =0
𝑠 𝑠 + 𝑎1 𝑠 + 𝑎2
21
PID Control of Motor Speed
• Ex. 4.5
Step disturbance input Step reference input
• Ex. 4.6
𝐷 𝑠 = 𝑘𝑃
−𝐵
𝑇𝑤 𝑠 = = 𝑠 0 𝑇𝑜,𝑤
𝑠 𝜏𝑠 + 1 + 𝐴𝑘𝑃 ℎ
𝐴𝑘𝑃 ℎ
Type 0, 𝐾0,𝑤 = −
𝐵
𝑘𝐼 −ℎ
If 𝐷 𝑠 = 𝑘𝑃 +
𝑠
−𝐵𝑠 1𝑇
𝑇𝑤 𝑠 = = 𝑠 𝑜,𝑤
𝑠 2 𝜏𝑠 + 1 + 𝑘𝑃 𝑠 + 𝑘𝐼 𝐴ℎ
𝐴𝑘𝐼 ℎ
Type 1, 𝐾𝑛,𝑤 = −
𝐵
23
Satellite Attitude Control
• Ex. 4.7
If 𝐷 𝑠 = 𝑘𝑃 + 𝑘𝐷 𝑠
1 0
𝑇𝑤 𝑠 = = 𝑠 𝑇𝑜,𝑤
𝐽𝑠 2 + 𝑘𝐷 𝑠 + 𝑘𝑃
Type 0, 𝐾0,𝑤 = 𝑘𝑃
𝑘𝐼
If 𝐷 𝑠 = 𝑘𝑃 + + 𝑘𝐷 𝑠
𝑠
𝑠 1
𝑇𝑤 𝑠 = 3 = 𝑠 𝑇𝑜,𝑤
𝐽𝑠 + 𝑘𝐷 𝑠 2 + 𝑘𝑃 𝑠 + 𝑘𝐼
Type 1, 𝐾𝑛,𝑤 = 𝑘𝐼
24
Ziegler-Nichols Tuning of PID Regulators
25
Quarter Decay Ratio
• Ex. 4.8
1
𝑅≅
90
𝐿 ≅ 13
After tuning
1
P control: 𝑘𝑃 = = 6.92
𝑅𝐿
0.9
PI control: 𝑘𝑃 = = 6.22
𝑅𝐿
𝐿 13
𝑇𝐼 = = = 43.3
0.3 0.3
26
Ziegler-Nichols Tuning of PID Regulators
Impulse response
• Ultimate Sensitivity Method
Neutral stable
27
Ultimate Sensitivity Method
• Ex. 4.9
28
Feedforward Control
𝐴
𝐺 𝑠 = 2
𝑠 + 𝑎1 𝑠 + 𝑎2
Possible
problems?
29
Integrator Anti-Windup
with anti-windup
• Ex. 9.9
31
Digital Control
• Block Diagram
Digital controller
• Ex.
𝑈 𝑧 23𝑧 − 21 5 23 21
𝐷𝑑 𝑧 = = 𝑢 𝑘+1 = 𝑢 𝑘 + 𝑒 𝑘+1 − 𝑒 𝑘
𝐸 𝑧 7𝑧 − 5 7 7 7
Difference equation
• Ex.
𝑈 𝑧 1.105𝑧 − 0.895
𝐷𝑑 𝑧 = = 1.4
𝐸 𝑧 𝑧−1
𝑢 𝑘 + 1 = 𝑢 𝑘 + 1.4 1.105𝑒 𝑘 + 1 − 0.895𝑒 𝑘
32
Continuous vs. Discrete Controllers
33
Homework #6
• Problem 4.6 Part (b) has been modified to compute the system 𝐾𝑣 in terms of
parameters 𝐾 and 𝑘𝑡 .
• Problem 4.29 Part (a) has been modified with 𝑒𝑠𝑠 ≤ 0.01 rad/sec.
• Problem 4.30
• Problem 4.32 The expression of 𝑣𝑎 has been modified using non-inverting
𝑡
amplifier, i.e. 𝑣𝑎 = 𝑘𝑃 𝑒 + 𝑘𝐼 0
𝑒𝑑𝑡.
34