Experiment 8: To Compute The Steady-State Error of Transfer Functions Using MATLAB
Experiment 8: To Compute The Steady-State Error of Transfer Functions Using MATLAB
E=R ( s )−C ( s )
If we call the output at steady state then the steady state error will be given by
c steady state=K e steady state
And
c steady state
K=
esteady state
This shows that by increasing the gain the steady state error reduces and output increases.
Consider the figure below where G(s) is a system given the error signal is given by
Lab Manual of ‘FEEDBACK CONTROL SYSTEMS’
E=R ( s )−C ( s )
C ( s ) =E ( s ) G ( s )
Substituting it in above mentioned equation. We get the following error equation
R( s)
E ( s )=
1+G( s)
By applying the final value theorem we can find the final value without using its inverse
transform.
Put the values and you will get the desired error for any input.
In-Lab Exercise
1.1: Find the steady state error for the following given system as transfer function excited by step
input? Use MATLAB to evaluate it.
5
G ( s )= 2
s + 7 s+ 10
MATLAB:
G = tf(5 , [1 7 10]);
Kp = dcgain(G)
Ess = 1 / (1+Kp)
Kp = 0.5000
Ess = 0.6667
Lab Manual of ‘FEEDBACK CONTROL SYSTEMS’
120(s+2)
1.2: For the given system G(s), find the steady state error for step input when G ( s )=
(s+3)( s+ 4)
.verify using MATLAB
MATLAB:
G = tf(120*([1 2]) , [1 7 12]);
Kp = dcgain(G)
Ess = 1 / (1+Kp)
Kp = 20
Ess = 0.0476
Lab Manual of ‘FEEDBACK CONTROL SYSTEMS’
1.3: Find the esteady state for 2 u(t), 5u(t) and 7u(t). verify using MATLAB
For 2u(t)
MATLAB:
G = tf(5 , [1 7 10]);
Kp = dcgain(G)
Ess = 2 / (1+Kp)
Kp = 0.5000
Ess = 1.3333
For 5u(t)
Lab Manual of ‘FEEDBACK CONTROL SYSTEMS’
MATLAB:
G = tf(5 , [1 7 10]);
Kp = dcgain(G)
Ess = 5 / (1+Kp)
Kp = 0.5
Ess = 3.3
For 7u(t)
MATLAB:
G = tf(5 , [1 7 10]);
Kp = dcgain(G)
Ess = 7 / (1+Kp)
Kp = 0.5
Ess = 4.66
Lab Manual of ‘FEEDBACK CONTROL SYSTEMS’
1.4: What happens to the steady state error when the system is excited by ramp and parabola
input? Give facts mathematically and verify using MATLAB. Give transfer function of q1.1.
For Ramp
s = tf('s');
G = tf(5 , [1 7 10]);
Kv = dcgain(s*G)
Ess = 1 / Kv
Kv = 0
Ess = Infinity
Lab Manual of ‘FEEDBACK CONTROL SYSTEMS’
For Parabola
s = tf('s');
G = tf(5 , [1 7 10]);
Ka = dcgain((s^2)*G)
Ess = 1 / Ka
Ka = 0
Ess=infinity
To comprehend
the theoretical
Correctly With minor Partially Poorly
basis for Doesn’t
1 analyzing the
answers all the mistakes answers answers the answers the
answer
questions all the questions questions questions
given control
systems