Linear Control System ] 2016-EE-113

LAB 09
OBJECT: Determine the step response of 2nd order system by creating a model in Simulink and by
generating a code in matlab.
TASK 01:-

Find the step response of following 2nd order open loop system by creating model in Simulink.
A) G1(s)= 8/S^2+16S+8
SIMULATION:
RESULT:

B) G2(s)=4/s^2+4
SIMULATION:

RESULT:

G(S)=1/s^2+2s+1
SIMULATION:

RESULT:

C) G(S)= 16/s^2+16s+16
SIMULATION:

RESULT:

SIR SYED UNIVERSITY OF ENGINEERING AND TECHNOLOGY Page 1

Linear Control System ] 2016-EE-113

TASK 02:-

Find the step response of following 2nd order closed loop system by creating model in Simulink.
A) G(s)=8/S^2+16S+8 , H(s)=1
SIMULATION:

RESULT:

G(s)=4/s^2+4, H(s)=0.5/5
SIMULATION:
RESULT:

G(s)=1/s^2+2s+1 , H(s)=2
SIMULATION:
RESULT:

G(s)=16/s^2+16s+16 , H(s)=1/s+2
SIMULATION:

SIR SYED UNIVERSITY OF ENGINEERING AND TECHNOLOGY Page 2

Linear Control System ] 2016-EE-113

RESULT:

TASK 03:-
A) For the given open loop system specifies the transient response and steady state response in
the step response graph
G(s)= 8/s^2+s+6
SOURCE CODE:
s=tf('s')
g=(8)/(s*s+s+6)
step(g)
gtext('Transient Area')
gtend('Steady State')
RESULT:

(b)for the given closed loop system specifies the transient response and steady state response in
the step response graph
G(s)=8/s^2+s+6 , H(s)=2
SOURCE CODE:
s=tf('s')
g=(8)/(s*s+s+6)
h=2
y=feedback(g,h)
step(y)
gtext('Transient Area')
gtend('Steady State')
RESULT:

TASK 4:-
A) Plot the step response of standard 2nd order system G(s)=Wn^2/s^2+2EWnS+Wn^2 With
different values of E= 0,0.3,0.5,1,2,2.5 & VALUE of Wn = 4 without using the plot command and
show the command in one window of figure.
SOURCE CODE:

t=[0.1 0.5 0.8 2]

for i=1:4; t1=0:0.001:5;
n1=[2.5]

SIR SYED UNIVERSITY OF ENGINEERING AND TECHNOLOGY Page 3

Linear Control System ] 2016-EE-113

d1=[t(i) 1]
g=tf(n1,d1)
subplot(4,1,i)
step(g,t1)
gtext('transient area')
gtext('steady state')
end
RESULT:

B) Plot the step response of standard 2nd order system G(s)=Wn^2/s^2+2EWnS+Wn^2 With
different values of Wn=1,2,3,4,5,6 & VALUE of E= 0.6 without using the plot command and
show the command in one window of figure.
SOURCE CODE:
k=[0.1 0.5 1 2]
for i=1:4; t1=0:0.001:5;
n1=[k(i)]
d1=[0.5 1]
g=tf(n1,d1)
subplot(4,1,i)
step(g,t1)
gtext('transient area')
gtext('steady state')
end
RESULT:

TASK NO 05:-
A) Generate a matlab code to plot the step response of the open loop 2nd order system by using
the plot command
G(s)= 9/S^2+2.5s+9
SOURCE CODE:
s=tf('s')
g=(9)/(s^2+2.5*s+9)
step(g)
gtext('Transient Area')
gtext('Steady State')
y=stepinfo(g)
plot(y)
RESULT:

SIR SYED UNIVERSITY OF ENGINEERING AND TECHNOLOGY Page 4

Linear Control System ] 2016-EE-113

B) G(s)= 9/s^2+2.5s+9, H(s)=1/s

SOURCE CODE:
s=tf('s')
g=(9)/(s^2+2.5*s+9)
h=1/s
y=feedback(g,h)
step(y)
gtext('Transient Area')
gtext('Steady State')
z=stepinfo(y)
plot(z)
RESULT:

LAB ASSIGNMENT
TASK 01:-

Find the step response of following 2nd order open loop system by creating model in Simulink.
A) G(s) =3/s^2+2s+3
SIMULATION:

RESULT:

G(s)=5/s^2+5
SIMULATION:

RESULT:

SIR SYED UNIVERSITY OF ENGINEERING AND TECHNOLOGY Page 5

Linear Control System ] 2016-EE-113

G(s)= 4/s^2+ 4s+4

SIMULATION:

RESULT:

G(s)=36/s^2+24s+36
SIMULATION:

RESULT:

TASK 02:-

Find the step response of following 2nd order closed loop system by creating model in Simulink.
A) G(s)=3/s^2+2s+ 3 H(s)=1
SIMULATION:

RESULT:

G(s)=5/s^2+5 , H(s)=2
SIMULATION:

SIR SYED UNIVERSITY OF ENGINEERING AND TECHNOLOGY Page 6

Linear Control System ] 2016-EE-113

RESULT:

G(s)= 4/s^2+ 4s+ 4 , H(s)= 3

SIMULATION:

RESULT:

G(s)= 36/s^2+ 24s+ 36, H(s)= 4

SIMULATION:

RESULT:

TASK 03:-
A) For the given open loop system specifies the transient response and steady state response in
the step response graph
G(s)=2/s^2+s+3
SOURCE CODE:
s=tf('s')
g=2/(1*s^2+1*s+3)
step(g)
gtext('Transient area')
gtext('steady state')
RESULT:

SIR SYED UNIVERSITY OF ENGINEERING AND TECHNOLOGY Page 7

Linear Control System ] 2016-EE-113

For the given closed loop system specifies the transient response and steady state response in
the step response graph
G(s)= 2/S^2+s+3 ,H(s)=3
SOURCE CODE:
s=tf('s')
g=2/(1*s^2+s+3)
h=3
y=feedback(g,h)
step(g)
gtext('Transient area')
gtext('steady state')
RESULT:
Step Response
1

0.9 Transient area

System: g
Peak amplitude: 0.925
0.8 Overshoot (%): 38.8
At time (seconds): 1.89 System: g
0.7 Settling time (seconds): 7.87
steady state

0.6 System: g
Final value: 0.667
Amplitude

System: g
0.5 Rise time (seconds): 0.759

0.4

0.3

0.2

0.1

0
0 2 4 6 8 10 12
Time (seconds)

TASK 04:-
A) Plot the step response of standard 2nd order system G(s)=Wn^2/s^2+2EWnS+Wn^2 With
different values of E= 0,0.2,0.4,1,2,2.8 & VALUE of Wn = 3 without using the plot command and
show the command in one window of figure
SOURCE CODE:

t=[0.2 0.4 1 3]
for i=1:4; t1=0:0.001:5;
n1=[1.5]
d1=[t(i) 1]
g=tf(n1,d1)
subplot(4,1,i)
step(g,t1)
gtext('transient area')
gtext('steady state')
end
RESULT:

SIR SYED UNIVERSITY OF ENGINEERING AND TECHNOLOGY Page 8

Linear Control System ] 2016-EE-113

B) Plot the step response of standard 2nd order system G(s)=Wn^2/s^2+2EWnS+Wn^2 With
different values of Wn=2,3,4,5,6,7 & VALUE of E= 0.7 without using the plot command and
show the command in one window of figure.
SOURCE CODE:
k=[0.2 0.4 0.8 3]
for i=1:4; t1=0:0.001:5;
n1=[k(i)]
d1=[0.8 1]
g=tf(n1,d1)
subplot(4,1,i)
step(g,t1)
gtext('transient area')
gtext('steady state')
end
RESULT:

TASK NO 05:-
A) Generate a matlab code to plot the step response of the open loop 2nd order system by using
the plot command
G(s)= 6/s^2+3.5s+6
SOURCE CODE:
s=tf('s')
g=(6)/(s^2+3.5*s+6)
step(g)
gtext('Transient Area')
gtext('Steady State')
y=stepinfo(g)
plot(y)
RESULT:

SIR SYED UNIVERSITY OF ENGINEERING AND TECHNOLOGY Page 9

Linear Control System ] 2016-EE-113

(b) Generate a matalb code to plot the step respone of the given closed loop 2nd order system by
using plot command.
G(s)= 6/s^2+3.5s+6 H(s)=5
SOURCE CODE:
s=tf('s')
g=(6)/(s^2+3.5*s+6)
h=5
y=feedback(g,h)
step(y)
gtext('Transient Area')
gtext('Steady State')
z=stepinfo(y)
plot(z)
RESULT:

SIR SYED UNIVERSITY OF ENGINEERING AND TECHNOLOGY Page 10