Department of Electrical Engineering
Control System – I Laboratory (EE593)
Simulation of Step Response of a Second Order System with Unity
Feedback and Calculation of Control System Specification like Time
constant, % Peak Overshoot, Settling Time etc. from the Response
Objective:
1. To determine the step response second order system using MATLAB
2. Calculation of transient response specifications from the step response
Minimum Software Requirement: MATLAB
Pre-Lab:
1. Find the pole zero mapping and transient response specification of the following systems
9 9 9 9
𝐺" 𝑠 = & ,𝐺 𝑠 = & ,𝐺 𝑠 = & ,𝐺 𝑠 = &
𝑠 + 6𝑠 + 9 & 𝑠 + 2𝑠 + 9 + 𝑠 + 9𝑠 + 9 " 𝑠 +9
2. Poles of the second order system are given as
𝑎 − 1 ± 𝑗2.8 𝑏 − 2 ± 𝑗2.8 𝑐 − 0.5 ± 𝑗2.8
3. Poles of the second order system are given as
𝑎 − 1 ± 𝑗2.8 𝑏 − 1 ± 𝑗4 𝑐 − 1 ± 𝑗
4. Damping factors of 1st three transfer function of 1. Are
𝑎 𝜉 = 0.3 (𝑏) 𝜉 = 0.5 𝑐 𝜉 = 0.1
Lab:
1. Find the pz map and step response of pre-lab1 using MATLAB to find transient
response specifications
2. Find the pz map and step response of pre-lab2 using MATLAB to find transient
response specifications
3. Find the pz map and step response of pre-lab3 using MATLAB to find transient
response specifications
4. Find the pz map and step response of pre-lab4 using MATLAB to find transient
response specifications
Post-Lab:
1. Record the poles and transient response specification of lab1 in tabular form
2. Record the poles and transient response specification of lab2 in tabular form
3. Record the poles and transient response specification of lab3 in tabular form
4. Record the poles and transient response specification of lab4 in tabular form
Conclusion:
1
Post-Lab:
1. Record the poles and transient response specification of lab1 in tabular form
Rise Peak Settling
SL Pole-Zero Maximum
System Time Time Time Comment
No. (sec) Overshoot
(sec) (sec) (sec)
1.
2.
3.
4.
2. Record the poles and transient response specification of lab2 in tabular form
Rise Peak Settling
SL Pole-Zero Maximum
System Time Time Time Comment
No. (sec) Overshoot
(sec) (sec) (sec)
1.
2.
3.
3. Record the poles and transient response specification of lab3 in tabular form
Rise Peak Settling
SL Pole-Zero Maximum
System Time Time Time Comment
No. (sec) Overshoot
(sec) (sec) (sec)
1.
2.
3.
4. Record the poles and transient response specification of lab4 in tabular form
Rise Peak Settling
SL Pole-Zero Maximum
System Time Time Time Comment
No. (sec) Overshoot
(sec) (sec) (sec)
1.
2.
3.
Conclusion:
MATLAB Input Command Code:
s=tf('s')
G1 = 9/(s^2+6*s+9)
G2 = 9/(s^2+2*s+9)
G3 = 9/(s^2+9*s+9)
G4 = 9/(s^2+9)
G_2a = 9/(s^2+2*s+9)
G_2b = 11.84/(s^2+4*s+11.84)
G_2c = 8.09/(s^2+s+8.09)
G_3a = 9/(s^2+2*s+9)
G_3b = 17/(s^2+2*s+17)
G_3c = 2/(s^2+2*s+2)
G_4a = 9/(s^2+1.8*s+9)
G_4b = 9/(s^2+3*s+9)
G_4c = 9/(s^2+0.6*s+9)
MATLAB Command Window Output:
>> EXP2C
s=
Continuous-time transfer function.
G1 =
9
-------------
s^2 + 6 s + 9
Continuous-time transfer function.
G2 =
9
-------------
s^2 + 2 s + 9
Continuous-time transfer function.
G3 =
9
-------------
s^2 + 9 s + 9
Continuous-time transfer function.
G4 =
9
-------
s^2 + 9
Continuous-time transfer function.
G_2a =
9
-------------
s^2 + 2 s + 9
Continuous-time transfer function.
G_2b =
11.84
-----------------
s^2 + 4 s + 11.84
Continuous-time transfer function.
G_2c =
8.09
--------------
s^2 + s + 8.09
Continuous-time transfer function.
G_3a =
9
-------------
s^2 + 2 s + 9
Continuous-time transfer function.
G_3b =
17
--------------
s^2 + 2 s + 17
Continuous-time transfer function.
G_3c =
2
-------------
s^2 + 2 s + 2
Continuous-time transfer function.
G_4a =
9
---------------
s^2 + 1.8 s + 9
Continuous-time transfer function.
G_4b =
9
-------------
s^2 + 3 s + 9
Continuous-time transfer function.
G_4c =
9
---------------
s^2 + 0.6 s + 9
Continuous-time transfer function.
>> ltiview
