Lecture - 4 2
Lecture - 4 2
FACULTY OF ENGINEERING
ELECTRICAL & ELECTRONICS ENGINEERING DEPT.
Instructor:
Assist. Prof. Dr. Mohammed ALKRUNZ
Electrical & Electronics Engineering Department
Roadmap,,,
Modeling:
Mathematical representation for the physical system to be analyzed in the next steps.
As it is given in the previous lectures: Transfer Function or State Space.
In time response analysis, transfer function modeling will be used.
Analysis:
Evaluate transient and steady state responses to see if they meet the performance
requirements (Satisfactory or not).
To be studied starting from this lecture.
Design:
Design controller that satisfies all the given transient and steady state responses.
Dr. M. Alkrunz
March 29, 2023 2
EEE362: Control Systems
Introduction,,,
𝑦 𝑡 = 𝑦𝑡 𝑡 + 𝑦𝑠𝑠 (𝑡)
Dr. M. Alkrunz
March 29, 2023 3
EEE362: Control Systems
Introduction,,,
Step Response
1.5
1
Amplitude
0.5
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4
Time (seconds)
As it is clear from the figure that the transient time is the time from zero (reset
state) till the response is settled down, while the steady state time is the remain
part after the transient died out.
Dr. M. Alkrunz
March 29, 2023 4
EEE362: Control Systems
Time Response – Input and Output,,,
Step Response
1.5
Input Output
System 1
Amplitude
𝒓(𝒕) 𝐺(𝑠) 𝒚(𝒕) 0.5
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4
Time (seconds)
We would like to analyze a system property by applying a test input 𝑟(𝑡) and
observing the time response 𝑦(𝑡).
If a step input is applied to the system as an input signal, then the output is
called “Step Response”.
Poles of a transfer function are the roots of the characteristic polynomial in the
denominator Q(s).
Zeros of a transfer function are the roots of the characteristic polynomial in the
numerator P(s).
Example:
𝒑𝒐𝒍𝒚𝒏𝒐𝒎𝒊𝒂𝒍 𝒐𝒇 (𝒛𝒆𝒓𝒐𝒔)
+6𝑠+5 𝑠2 (𝑠 + 1)(𝑠 + 5)
𝐺 𝑠 = 3 =
𝑠 + 9 𝑠 2 + 26 𝑠 + 24 (𝑠 + 2)(𝑠 + 3)(𝑠 + 4)
𝒑𝒐𝒍𝒚𝒏𝒐𝒎𝒊𝒂𝒍 𝒐𝒇 (𝒑𝒐𝒍𝒆𝒔)
Poles are: 𝑠 = −2, 𝑠 = −3, and 𝑠 = −4
Zeros are: 𝑠 = −1 and 𝑠 = −5
Dr. M. Alkrunz
March 29, 2023 6
EEE362: Control Systems
Time Response - s-Plane,,,
𝒊𝒎𝒈(𝒔)
LHP RHP
𝑠 ≡ 𝑟𝑒𝑎𝑙 𝑝𝑎𝑟𝑡 + 𝑖𝑚𝑎𝑔𝑖𝑛𝑎𝑟𝑦 𝑝𝑎𝑟𝑡
𝑠 = 𝜎 + 𝑗𝑤
𝑹𝒆(𝒔)
𝑳𝒆𝒇𝒕 𝑯𝒂𝒍𝒇 𝑷𝒍𝒂𝒏𝒆 𝑹𝒊𝒈𝒉𝒕 𝑯𝒂𝒍𝒇 𝑷𝒍𝒂𝒏𝒆
(LHP) (RHP)
𝑗𝑤
𝒔 − 𝒑𝒍𝒂𝒏𝒆 𝒑𝒐𝒍𝒆𝒔
In the s-plane: × 𝐨 𝐨 × 𝜎
Poles are indicated as “x”
Zeros are indicated as “o” 𝒛𝒆𝒓𝒐𝒔
1. First Order Systems [First order means that we have just one pole]
2. Second Order Systems [Second order means that we have two poles]
𝒌𝒅𝒄
𝒌𝒅𝒄 𝝉 𝒌𝒅𝒄 𝒂
𝑮 𝒔 = = =
𝝉𝒔+𝟏 𝒔+𝟏 𝝉 𝒔+𝒂
where:
𝑘𝑑𝑐 is the dc gain
𝜏 is the time constant
𝑎=1 𝜏
Dr. M. Alkrunz
March 29, 2023 8
EEE362: Control Systems
Time Response – System Orders,,,
𝒌𝒅𝒄
𝒌𝒅𝒄 𝝉 𝒌𝒅𝒄 𝒂
𝑮 𝒔 = = =
𝝉𝒔+𝟏 𝒔+𝟏 𝝉 𝒔+𝒂
4
Settling Time “𝑻𝒔 ”: 𝑇𝑠 = or 4𝜏 𝝉
𝑎
2.2
Rise Time “𝑻𝒓 ”: 𝑇𝑟 = or 2.2𝜏
𝑎
Example:
50
𝐺 𝑠 =
𝑠 + 50
Find the dc gain “𝑘𝑑𝑐 ”,time constant “𝜏”, settling time “𝑇𝑠 ”, and rise time “𝑇𝑟 ”.
Solution:
50 𝑘𝑑𝑐 𝑎
The system is a first order system, then: 𝐺 𝑠 = =
𝑠 + 50 𝑠 + 𝑎
𝑘𝑑𝑐 𝑎 = 50, since 𝑎 = 50 then 𝑘𝑑𝑐 = 1
1 1
𝜏= = = 0.02 𝑠𝑒𝑐
𝑎 50
𝑇𝑠 = 4 𝜏 = 4 × 0.02 = 0.08 𝑠𝑒𝑐
𝑇𝑟 = 2.2 𝜏 = 2.2 × 0.02 = 0.044 𝑠𝑒𝑐
Dr. M. Alkrunz
March 29, 2023 10
EEE362: Control Systems
Time Response – System Orders,,,
𝒌𝒅𝒄 𝝎𝒏 𝟐
𝑮 𝒔 = 𝟐
𝒔 + 𝟐𝜻𝝎𝒏 𝒔 + 𝝎𝒏 𝟐
where:
𝜎
𝜔𝑛 is the natural frequency (rad/sec)
𝜁 is the damping ratio
𝜔𝑑 is the damping frequency (rad/sec)
Dr. M. Alkrunz
March 29, 2023 11
EEE362: Control Systems
Time Response – System Orders,,,
The above figure illustrates the transient and steady state responses analysis.
Let us investigate that in details.
Dr. M. Alkrunz
March 29, 2023 12
EEE362: Control Systems
Time Response – System Orders,,,
Dr. M. Alkrunz
March 29, 2023 13
EEE362: Control Systems
Time Response – System Orders,,,
Peak Overshoot“𝑶𝑺%”:
−𝜋𝜁
𝑦𝑝 −𝑦𝑠𝑠
𝑂𝑆% = 𝑒 1−𝜁 2 × 100% Or 𝑂𝑆% = × 100%
𝑦𝑠𝑠
Dr. M. Alkrunz
March 29, 2023 14
EEE362: Control Systems
Time Response – System Orders,,,
Example:
36
𝐺 𝑠 =
𝑠 2 + 4.2 𝑠 + 36
Find 𝜁 and 𝜔𝑛 .
Solution:
36 𝑘𝑑𝑐 𝜔𝑛 2
𝐺 𝑠 = 2 =
𝑠 + 4.2 𝑠 + 36 𝑠 2 + 2𝜁𝜔𝑛 𝑠 + 𝜔𝑛 2
𝜔𝑛 2 = 36 ≫≫≫ 𝜔𝑛 = 36 = 𝟔 𝒓𝒂𝒅/𝒔𝒆𝒄
4.2
2𝜁𝜔𝑛 = 4.2 ≫≫≫ 𝜁 = = 𝟎. 𝟑𝟓
2𝜔𝑛
Dr. M. Alkrunz
March 29, 2023 15
EEE362: Control Systems
Time Response – Types of Step Response,,,
Here, the types of step response for the second order system:
Dr. M. Alkrunz
March 29, 2023 16
EEE362: Control Systems
Time Response – Types of Step Response,,,
Poles of the transfer function are real and repeated (two real poles at same
position).
𝑠1,2 = −𝜎 = −𝜁𝜔𝑛
Dr. M. Alkrunz
March 29, 2023 17
EEE362: Control Systems
Time Response – Types of Step Response,,,
Poles of the transfer function are real and distinct (two real poles at different
positions).
𝑠1 = 𝜎1 and 𝑠2 = 𝜎2
𝜎1
𝜎2
4. UndampeStep Response:
𝜁=0
Since there is no real part, high and periodic oscillation will be produced.
1. Overshoot has direct relationship with the damping ratio “𝜁”. Therefore, it is
important to note that damping ratio “𝜁” changes between 0 and 1. Thus, the
overshoot is changing from 0% (𝜁 = 1) to 100% (𝜁 = 0). See the figure,
𝜻 Step response
0 < 𝜁 < 1 Underdamped
𝜁=1 Critically damped
𝜁=0 Undamped
𝜁>1 Overdamped
Dr. M. Alkrunz
March 29, 2023 20
EEE362: Control Systems
Time Response – Types of Step Response,,,
Some Important Comments
2. The system is stable if all poles are located in the left half plane (LHP), and if
any pole of them is located in the right half plane (RHP), directly the system
will be unstable. Namely, the system needs at least one pole in the (RHP) to be
unstable while it needs all the poles to be located in (LHP) to be stable.
𝒊𝒎𝒂𝒈𝒊𝒏𝒂𝒓𝒚
𝒂𝒙𝒊𝒔
𝒔 − 𝒑𝒍𝒂𝒏𝒆
Dr. M. Alkrunz
March 29, 2023 21
EEE362: Control Systems
Time Response,,,
Another possible way to find 𝝎𝒏 and 𝜻 from the pole location in s-plane
= 𝜎 2 + 𝜔𝑑 2 = 𝜁 2 𝜔𝑛 2 + 𝜔𝑛 2 1 − 𝜁 2 = 𝜔𝑛 2 = 𝜔𝑛
×
𝜎 𝜁𝜔𝑛
𝝎𝒏 cos 𝛽 = = =𝜁
𝜔𝑛 𝜔𝑛
𝜁 = cos 𝛽
Dr. M. Alkrunz
March 29, 2023 22
EEE362: Control Systems
Open Loop and Closed Loop Systems,,,
𝑌(𝑠) 𝑮(𝒔)
𝐺 𝑠 =
𝑋(𝑠) 𝑿(𝒔) 𝒀(𝒔)
𝑿(𝒔) + 𝒀(𝒔)
𝑌(𝑠) 𝐺(𝑠) 𝑮(𝒔)
𝑇 𝑠 = =
𝑋(𝑠) 1 + 𝐺 𝑠 𝐻(𝑠) −
𝐇(𝒔)
Dr. M. Alkrunz
March 29, 2023 23
EEE362: Control Systems
Examples,,,
Example 𝑮(𝒔)
𝑿(𝒔) 𝒀(𝒔)
𝒀(𝒔) 𝟏𝟖
𝑮 𝒔 = =
𝑿(𝒔) 𝒔𝟐 + 𝟒 𝒔 + 𝟗
Solution:
(a)
18 𝑘𝑑𝑐 𝜔𝑛 2
𝐺 𝑠 = 2 =
𝑠 + 4 𝑠 + 9 𝑠 2 + 2𝜁𝜔𝑛 𝑠 + 𝜔𝑛 2
Dr. M. Alkrunz
March 29, 2023 24
EEE362: Control Systems
Examples,,,
Solution (Continue):
18 𝑘𝑑𝑐 𝜔𝑛 2
𝐺 𝑠 = 2 =
𝑠 + 4 𝑠 + 9 𝑠 2 + 2𝜁𝜔𝑛 𝑠 + 𝜔𝑛 2
𝜔𝑛 2 = 9 ≫≫ 𝜔𝑛 = 9 = 𝟑 𝒓𝒂𝒅/𝒔𝒆𝒄
4 4
2𝜁𝜔𝑛 = 4 ≫≫ 𝜁 = = = 𝟎. 𝟔𝟕
2𝜔𝑛 2 3
4 4
𝑡𝑠 = = = 2 𝑠𝑒𝑐
𝜁𝜔𝑛 (0.67)(3)
−𝜋𝜁
𝑂𝑆% = exp × 100% = 6 %
1− 𝜁2
Dr. M. Alkrunz
March 29, 2023 25
EEE362: Control Systems
Examples,,,
Solution (Continue):
(b)
Since 𝜁 = 0.67 which means that 0 < 𝜁 < 1, then the system is underdamped.
If you check the poles of the system, you will find them complex poles with
real and imaginary part which mean again the response is underdamped.
(c)
3 18 (3)(18)
𝑦𝑠𝑠 = lim 𝑠 𝑌(𝑠) = lim 𝑠 𝑋 𝑠 𝐺 𝑠 = lim 𝑠 × × 2 = =𝟔
𝑠→0 𝑠→0 𝑠→0 𝑠 𝑠 +4𝑠+9 9
Dr. M. Alkrunz
March 29, 2023 26
EEE362: Control Systems
Examples,,,
Example
Solution:
𝜔𝑑 = 7 and 𝜎 = 3
𝜔𝑛 = 𝜎 2 + 𝜔𝑑 2 = 32 + 72 = 7.616 rad/sec
𝜎 3
𝜁 = cos 𝜃 = = = 0.394
𝜔𝑛 7.616
4 4
𝑡𝑠 = = = 1.33 𝑠𝑒𝑐
𝜁𝜔𝑛 0.349 × 7.616
−𝜋𝜁
𝑂𝑆% = 𝑒𝑥𝑝 × 100% = 26 %
1− 𝜁2
Dr. M. Alkrunz
March 29, 2023 27
EEE362: Control Systems
Examples,,,
Example
+ + 500 1
𝑹(𝒔) 𝑘 𝐺 𝑠 = 𝒀(𝒔)
𝑠+5 20 𝑠
− −
𝑘𝑡
Dr. M. Alkrunz
March 29, 2023 28
EEE362: Control Systems
Examples,,,
Solution (Continue):
𝑨(𝒔)
+ + 500 1
𝑹(𝒔) 𝑘 𝐺 𝑠 = 𝒀(𝒔)
𝑠+5 20 𝑠
− −
𝑘𝑡
500
𝐺(𝑠) 𝑠+5 500
𝐴 𝑠 = = =
1 + 𝑘𝑡 𝐺(𝑠) 1 + 𝑘 500 𝑠 + 500𝑘𝑡 + 5
𝑡𝑠+5
Dr. M. Alkrunz
March 29, 2023 29
EEE362: Control Systems
Examples,,,
Solution (Continue):
+ 1 𝒀(𝒔) +
𝑹(𝒔) 𝑘 𝐴(𝑠) 𝐵(𝑠)
20 𝑠
− −
1 500𝑘 25𝑘
𝐵 𝑠 =𝑘𝐴 𝑠 = = 2
20 𝑠 20𝑠(𝑠 + 500𝑘𝑡 + 5 ) 𝑠 + 500𝑘𝑡 + 5 𝑠
Dr. M. Alkrunz
March 29, 2023 30
EEE362: Control Systems
Examples,,,
Solution (Continue):
𝑹(𝒔) 25𝑘 𝒀(𝒔)
𝑠 2 + 500𝑘𝑡 + 5 𝑠 + 25𝑘
3.147 1 − 𝜁 2 = 𝜋𝜁
4.3 −𝜋𝜁 2
= exp 3.147 1 − 𝜁2 = 𝜋 2 𝜁2
100 1 − 𝜁2
𝜁 2 𝜋 2 + 9.9 = 9.9
4.3 −𝜋𝜁
ln =
100 1 − 𝜁2 9.9
𝜁= = 𝟎. 𝟕𝟎𝟕
𝜋 2 + 9.9
−𝜋𝜁
−3.147 =
1 − 𝜁2
Dr. M. Alkrunz
March 29, 2023 31
EEE362: Control Systems
Examples,,,
Solution (Continue):
4
𝑡𝑠 = = 0.628
𝜁𝜔𝑛
4 4
𝜔𝑛 = = = 𝟗 𝑟𝑎𝑑/𝑠𝑒𝑐
0.628 × 𝜁 0.628 × 0.707
Dr. M. Alkrunz
March 29, 2023 32
EEE362: Control Systems
Examples,,,
Solution (Continue):
25𝑘 𝑘𝑑𝑐 𝜔𝑛 2
=
𝑠 2 + 500𝑘𝑡 + 5 𝑠 + 25𝑘 𝑠 2 + 2𝜁𝜔𝑛 𝑠 + 𝜔𝑛 2
25𝑘 = 𝜔𝑛 2
𝜔𝑛 2 9 2
𝑘= = = 𝟑. 𝟐𝟓
25 25
500𝑘𝑡 + 5 = 2𝜁𝜔𝑛
2𝜁𝜔𝑛 − 5 2 0.707 9 − 5
𝑘𝑡 = = = 𝟎. 𝟎𝟏𝟓𝟒
500 500
Dr. M. Alkrunz
March 29, 2023 33
EEE362: Control Systems
*Questions
Dr. M. Alkrunz
March 29, 2023 34
EEE362: Control Systems
*Thanks
Dr. M. Alkrunz
March 29, 2023 35
EEE362: Control Systems