Digital and Non-Linear Control: Time Domain Analysis
Digital and Non-Linear Control: Time Domain Analysis
A
A t0
(t )
0 t0
0 t
characteristic of actual
A
input signal.
A t0 t
u( t ) 0
0 t0
At t0
r(t )
r(t)
0 t0
ramp signal with slope A
At 2 p(t)
t0
p(t ) 2
0 t0
parabolic signal with slope A
p(t)
• Impulse
A t0
(t )
0 t0
L{ (t )} ( s) A
• Step
A t0
u( t )
0 t0
A
L{u(t )} U ( s )
S
Laplace Transform of Test Signals
• Ramp At t0
r(t )
0 t0
A
L{ r(t )} R( s )
s2
• Parabolic At 2
t0
p(t ) 2
0 t0
A
L{ p (t )} P( s ) 3
S
Time Response of Control Systems
• Time response of a dynamic system response to an input
expressed as a function of time.
System
Step Input
5
• The response of the
0
0 2 4 6 8 10 12 14 16 18 20
Time (sec)
Time Response of Control Systems
• Transient response depend upon the system poles only and not
on the type of input.
K
R(s ) C(s )
1
Ts 1
t
0
R( s ) ( s ) 1
K
C( s )
Ts 1
Impulse Response of 1st Order System
K
C( s )
Ts 1
• Re-arrange following equation as
K /T
C( s )
s 1/ T
1
c(t)
0.5
0
0 2 4 6 8 10
Time
Step Response of 1st Order System
• Consider the following 1st order system
K
R(s ) C(s )
Ts 1
1
R( s ) U ( s )
s
K
C( s )
sTs 1
• In order to find out the inverse Laplace of the above equation, we
need to break it into partial fraction expansion (page 867 in the
Textbook) K KT
C( s )
s Ts 1
Step Response of 1st Order System
1 T
C( s ) K
s Ts 1
• Taking Inverse Laplace of above equation
c(t ) K u(t ) et / T
c(t ) K 1 e
• Where u(t)=1
t / T
c(t ) K 1 e1 0.632K
Step Response of 1st Order System
• If K=10 and T=1.5s then
c(t ) K 1 e t / T
K*(1-exp(-t/T))
11
10
9 Step Response
8
steady state output 10
7 D.C Gain K
63% Input 1
6
c(t)
2
Unit Step Input
1
0
0 1 2 3 4 5 6 7 8 9 10
Time
Step Response of 1st order System
• System takes five time constants to reach its
final value.
Step Response of 1st Order System
• If K=10 and T=1, 3, 5, 7
c(t ) K 1 e t / T
K*(1-exp(-t/T))
11
10
T=1s
9
8 T=3s
7
T=5s
6
c(t)
5 T=7s
4
3
2
1
0
0 5 10 15
Time
Step Response of 1st Order System
• If K=1, 3, 5, 10 and T=1
c(t ) K 1 e t / T
K*(1-exp(-t/T))
11
10
K=10
9
8
7
6
K=5
c(t)
5
4
K=3
3
2
K=1
1
0
0 5 10 15
Time
Relation Between Step and impulse
response
• The step response of the first order system is
c(t ) K 1 e t / T K Ket / T
• v(t)=0.5v0 t = 0.69RC
– i.e., delay = 0.69RC (50% delay)
v(t)=0.1v0 t = 0.1RC
v(t)=0.9v0 t = 2.3RC
– i.e., rise time = 2.2RC (if defined as time from 10% to 90% of Vdd)
• For simplicity, industry uses
TD = RC (= Elmore delay)
Elmore Delay
1. 50%-50%
point delay
Delay 2. Delay=0.69
RC
Example 1
• Impulse response of a 1st order system is given below.
c(t ) 3e 0.5t
• Find out
– Time constant T
– D.C Gain K
– Transfer Function
– Step Response
Example 1
• The Laplace Transform of Impulse response of a
system is actually the transfer function of the system.
• Therefore taking Laplace Transform of the impulse
response given by following equation.
c(t ) 3e 0.5t
3 3
C( s ) 1 (s)
S 0.5 S 0.5
C( s ) C( s ) 3
( s ) R( s ) S 0.5
C( s ) 6
R( s ) 2 S 1
Example 1
• Impulse response of a 1st order system is given below.
c(t ) 3e 0.5t
• Find out
– Time constant T=2
– D.C Gain K=6
– Transfer Function C ( s ) 6
R( s ) 2 S 1
– Step Response
Example 1
• For step response integrate impulse response
c(t ) 3e 0.5t
0.5t
c( t )dt 3 e dt
cs (t ) 6e 0.5t C
0 6e 0.50 C
C6
cs (t ) 6 6e 0.5t
Example 1
• If initial conditions are not known then partial fraction
expansion is a better choice
C( s ) 6
R( s ) 2 S 1
1
since R( s ) is a step input , R( s )
s
6
C( s )
s2S 1
6 A B
s2S 1 s 2s 1
6 6 6
s2S 1 s s 0.5
c(t ) 6 6e 0.5t
Ramp Response of 1st Order System
• Consider the following 1st order system
K
R(s ) C(s )
Ts 1
1
R( s )
s2
K
C( s )
s 2 Ts 1
• The ramp response is given as
c(t ) K t T Tet / T
Parabolic Response of 1st Order System
• Consider the following 1st order system
K
R(s ) C(s )
Ts 1
1 K
R( s ) Therefore, C( s )
s 3
s 3 Ts 1
Practical Determination of Transfer
Function of 1st Order Systems
• Often it is not possible or practical to obtain a system's
transfer function analytically.
• With a step input, we can measure the time constant and the
steady-state value, from which the transfer function can be
calculated.
Practical Determination of Transfer
Function of 1st Order Systems
• If we can identify T and K empirically we can obtain the
transfer function of the system.
C( s ) K
R( s ) Ts 1
Practical Determination of Transfer Function
of 1st Order Systems
• For example, assume the unit
step response given in figure. K=0.72
• From the response, we can
measure the time constant, that
is, the time for the amplitude to
reach 63% of its final value.
• Since the final value is about T=0.13s
Unit Step
Step Response
t
td
First Order System With Delays
Step Response
C (s) 10 2 s
e
R( s ) 3s 1 10
C (s)
10
e 2 s K 10
s (3s 1)
8
L1[e s F ( s )] f (t )u (t )
10 10 2 s
L1[( )e ]
s s 1/ 3 6
Amplitude
t d 2s
T 3s
0
0 5 10 15
Time (sec)
Second Order System
• We have already discussed the affect of location of poles and zeros on
the transient response of 1st order systems.
C( s ) n2
2
R( s ) s 2 n s n2
C( s ) 4
2
R( s ) s 2s 4
C( s ) n2
2
R( s ) s 2 n s n2
n2 4 n 2 2 n s 2s
n 1
s 2 2 n s n2 s 2 2s 4
0.5 46
Introduction
C( s ) n2
2
R( s ) s 2 n s n2
n n 2 1
n n 2 1
47
Introduction
n n 2 1
n n 2 1
• According the value of , a second-order system can be set into
one of the four categories (page 169 in the textbook):
1. Overdamped - when the system has two real distinct poles ( >1).
jω
δ
-c -b -a
48
Introduction
n n 2 1
n n 2 1
• According the value of , a second-order system can be set into
one of the four categories (page 169 in the textbook):
2. Underdamped - when the system has two complex conjugate poles (0 < <1)
jω
δ
-c -b -a
49
Introduction
n n 2 1
n n 2 1
• According the value of , a second-order system can be set into
one of the four categories (page 169 in the textbook):
δ
-c -b -a
50
Introduction
n n 2 1
n n 2 1
• According the value of , a second-order system can be set into
one of the four categories (page 169 in the textbook):
4. Critically damped - when the system has two real but equal poles ( = 1).
jω
δ
-c -b -a
51
Underdamped System
For 0< <1 and ωn > 0, the 2nd order system’s response due to a
unit step input is as follows.
Important timing characteristics: delay time, rise time, peak
time, maximum overshoot, and settling time.
52
Delay Time
• The delay (td) time is the time required for the response to
reach half the final value the very first time.
53
Rise Time
• The rise time is the time required for the response to rise from 10%
to 90%, 5% to 95%, or 0% to 100% of its final value.
• For underdamped second order systems, the 0% to 100% rise time is
normally used. For overdamped systems, the 10% to 90% rise time is
commonly used.
Peak Time
• The peak time is the time required for the response to reach
the first peak of the overshoot.
55
55
Maximum Overshoot
The maximum overshoot is the maximum peak value of the
response curve measured from unity. If the final steady-state
value of the response differs from unity, then it is common to
use the maximum percent overshoot. It is defined by
56
Settling Time
• The settling time is the time required for the response curve
to reach and stay within a range about the final value of size
specified by absolute percentage of the final value (usually 2%
or 5%).
57
Step Response of underdamped System
C( s ) n2 Step Response n2
2 C( s )
R( s ) s 2 n s n2
s s 2 2 n s n2
• The partial fraction expansion of above equation is given as
1 s 2 n
C( s ) 2
s s 2 n s n2
n2 1 2
1 s 2 n
C( s ) 2
s 2 n 2 s s 2 n s 2 n2 n2 2 n2
1 s 2 n
C( s )
s s n 2 n2 1 2 58
Step Response of underdamped System
1 s 2 n
C( s )
s s n 2 n2 1 2
• Above equation can be written as
1 s 2 n
C( s )
s s n 2 d2
• Where d n 1 2 , is the frequency of transient oscillations
and is called damped natural frequency.
n 1 2
1 s n 1 2
C( s )
s s n d
2 2
s n 2 d2
1 s n d
C( s )
s s n d
2 2
1 2 s 2 2
n d
n t
c(t ) 1 e cos d t e nt sin d t
1 2
60
Step Response of underdamped System
n t
c(t ) 1 e cos d t e nt sin d t
1 2
c(t ) 1 e nt cos d t sin d t
1 2
• When 0
d n 1 2
n
c(t ) 1 cos n t
61
Step Response of underdamped System
c(t ) 1 e nt cos d t sin d t
1 2
if 0.1 and n 3
1.8
1.6
1.4
1.2
0.8
0.6
0.4
0.2
62
0
0 2 4 6 8 10
Step Response of underdamped System
c(t ) 1 e nt cos d t sin d t
1 2
if 0.5 and n 3
1.4
1.2
0.8
0.6
0.4
0.2
63
0
0 2 4 6 8 10
Step Response of underdamped System
c(t ) 1 e nt cos d t sin d t
1 2
if 0.9 and n 3
1.4
1.2
0.8
0.6
0.4
0.2
0 64
0 2 4 6 8 10
Step Response of underdamped System
c(t ) 1 e nt cos d t sin d t
1 2
65
S-Plane (Underdamped System)
Since 𝜔2 𝜁 2 − 𝜔2 𝜁 2 − 1 = 𝜔2 , the distance
n n 2 1 from the pole to the origin is 𝜔 and 𝜁 = 𝑐𝑜𝑠𝛽
n n 2 1
66
Analytical Solution
• Page 171 in the textbook c(t ) 1 e nt
cos d t
1 2
sin d t
𝜋−𝛽
• Rise time: set c(t)=1, we have 𝑡𝑟 = d n 1 2
𝜔𝑑
𝑑𝑐(𝑡) 𝜋
• Peak time: set = 0, we have 𝑡𝑝 =
𝑑𝑡 𝜔𝑑
• Maximum overshoot: M𝑝 = 𝑐 𝑡𝑝 − 1 =
𝑒 −(𝜁𝜔/𝜔𝑑 )𝜋 (for unity output)
• Settling time: the time for the outputs always
4
within 2% of the final value is approximately
𝜁𝜔
Empirical Solution Using MATLAB
• Page 242 in the textbook
Steady State Error
• If the output of a control system at steady state does not
exactly match with the input, the system is said to have
steady state error
• The transfer function between the error signal E(s) and the
input signal R(s) is E ( s ) 1
R( s ) 1 G( s )
• The final-value theorem provides a convenient way to find
the steady-state performance of a stable system.
• Since E(s) is
• The steady state error is
Static Error Constants
• The static error constants are figures of merit of
control systems. The higher the constants, the
smaller the steady-state error.
• In a given system, the output may be the position,
velocity, pressure, temperature, or the like.
• Therefore, in what follows, we shall call the output
“position,” the rate of change of the output
“velocity,” and so on.
• This means that in a temperature control system
“position” represents the output temperature,
“velocity” represents the rate of change of the
output temperature, and so on.
Static Position Error Constant (Kp)
• The steady-state error of the system for a unit-step input is
100 ( s 2)( s 5)
R(S) C(S)
2
s ( s 8)( s 12 )
-
Example 2
100 ( s 2)( s 5)
G( s )
s 2 ( s 8)( s 12 )
K p lim G( s )
s 0 K v lim sG( s )
s 0
100 ( s 2 )( s 5)
K p lim 2 100 s( s 2 )( s 5)
s 0 s ( s 8)( s 12 ) K v lim 2
s 0 s ( s 8)( s 12 )
Kp
Kv
0
0
0.09