Practical-0
Simulate open loop systems with step input for system identification of plant model various plant
transfer function.
Consider system 1/(s+1)^8 :
Take approximation models of step response are two parameter, three parameter, four
parameter models.
Step response of all model given in figure
Step Response
1
System: g1
0.9 Time (seconds): 20
Amplitude: 0.999
0.8
System: g1
Time (seconds): 8.63
0.7 Amplitude: 0.631
0.6
0.5
magnitude
System: g1
0.4
Time (seconds): 6.19
Amplitude: 0.283
0.3
0.2
0.1
0
0 2 4 6 8 10 12 14 16 18 20
FIGURE 1 time (seconds)
MATLAB program for system 1/(s+1)^8 given below:
clc;
clear;
syms s
num=[1]
den=poly([-1 -1 -1 -1 -1 -1 -1 -1])
g1=tf(num,den);
printsys(num,den);
[n2,d2]=feedback(num,den,1,1);
g2=tf(n2,d2);
printsys(n2,d2);
step(g1);
xlabel('time');
ylabel('magnitude');
s=stepinfo(g1)
Analysis :
Two parameter Model: K/(Tars+1)
Here Tar is average residence time approximated as time is taken by system to reach
steady state value of step.
K is the gain of system which is equal to steady state value of step.
Here from the graph we get K=1, Tar = 20 sec.
Put value in the model we find : 1/(20s+1) model.
Three parameter Model: K/(Ts+1) *e-LS
Here T is time constant approximated as time t1 and time t2 is taken by system to reach
28.3% and 63.2% steady state value of step respectively.
Than T = 3/2 (t2-t1) and L = t2-T
L is lag or transportation delay of system
K is the gain of system which is equal to steady state value of step.
Here from the graph we get K=1, t1 = 6.19 sec, t2 = 8.64 sec.
From equation of T & L we get T = 3.67 sec, L = 4.965 sec
Put value in the model we find : 1/(3.67s+1) *e-4.965 S model.
Three parameter Model: K/(t1s+1) (t2s+1) *e-LS
Here T is time constant approximated as time t1 and time t2 is taken by system to reach
28.3% and 63.2% steady state value of step respectively.
L is lag or transportation delay of system
K is the gain of system which is equal to steady state value of step.
Here from the graph we get K=1, t1 = 6.19 sec, t2 = 8.64 sec.
Put value in the model we find : 1/(6.19s+1) (8.64s+1) *e-4.965 S model.