4 - Time Domain Analysis
4 - Time Domain Analysis
Slide 2
4.1 Time Response of Continuous Systems
• In order to analyse 𝑐 𝑡 and 𝑐 𝑡 , let us define
𝑁(𝑠)
𝐺 𝑠 =
𝐷(𝑠)
• Poles of 𝐺 𝑠
• the values of the Laplace transform variable, s, that cause the transfer function to become
infinite or
• any roots of the denominator of the transfer function that are common to roots of the
numerator.
• Zeros of 𝐺 𝑠
• the values of the Laplace transform variable, s, that cause the transfer function to become
zero, or
• any roots of the numerator of the transfer function that are common to roots of the
denominator.
Slide 3
4.1 Time Response of Continuous Systems
• Time response of a system is defined in two parts
• Transient State / Natural Response
• Steady State/ Forced Response
𝑐 𝑡 =𝑐 𝑡 +𝑐 𝑡
• In order to analyse 𝑐 𝑡 and 𝑐 𝑡 , we need to also define the following deterministic test
signals
• Step function Input
• Ramp function Input
• Parabolic function Input
Slide 4
4.1 Time Response of Continuous Systems
Slide 5
Motivating Example
Slide 6
4.2 Time Response of LTI Systems
A) Time Response of first order systems
• Let 1st order dynamic system described as follows:
𝑑𝑐
= 𝑎𝑐 𝑡 + 𝑏𝑢 𝑡 , 𝑎 < 0
𝑑𝑡
or
𝑑𝑐 1
= − 𝑐 𝑡 + 𝑏𝑢 𝑡 , 𝜏 > 0
𝑑𝑡 𝜏
• We can define:
𝑏𝜏 𝐾
𝐺 𝑠 = ≡
𝜏𝑠 + 1 𝜏𝑠 + 1
Slide 7
A) Time Response of first order systems
• The step response of this system is as follows:
𝐾𝑀 𝐾𝑀 𝐾𝑀
𝐶 𝑠 = = −
𝑠(𝜏𝑠 + 1) 𝑠 1
𝑠+𝜏
𝐾𝑀
𝑐 𝑡 = 𝐾𝑀(1 − 𝑒 / )
• Let us examine the significance of 𝜏
0.632𝐾𝑀
𝑐 𝑡 = 𝐾𝑀 1 − 𝑒 = 0.632𝐾𝑀
• 𝜏 defines the speed of the system response
• 𝜏: Time Constant
• The slope of the system response at t=0
𝑑𝑐 𝑡 𝐾𝑀
= 𝑒
𝑑𝑡 𝜏
𝑑𝑐 0 𝐾𝑀 𝐾𝑀
= 𝑒 =
𝑑𝑡 𝜏 𝜏
Slide 8
A) Time Response of first order systems
• The step response of this system is as follows:
𝐾𝑀 𝐾𝑀 𝐾𝑀
𝐶 𝑠 = = −
𝑠(𝜏𝑠 + 1) 𝑠 1
𝑠+𝜏
𝐾𝑀
𝑐 𝑡 = 𝐾𝑀(1 − 𝑒 / )
• Let us examine the significance of 𝜏
0.632𝐾𝑀
𝑐 𝑡 = 𝐾𝑀 1 − 𝑒 = 0.632𝐾𝑀
• 𝜏 defines the speed of the system response
• 𝜏: Time Constant
• The exponential terms can be neglected after
t = 4𝜏. Thus, the settling time of the system is:
𝑇 ≅ 4𝜏
Slide 9
B) Time Response of second order systems
• Motivating Example
Slide 10
B) Time Response of second order systems
• Motivating Example
Slide 12
B) Time Response of second order systems
• Let 2nd order dynamic system described as follows:
𝐾𝑤
𝐺 𝑠 =
𝑠 + 2𝜁𝑤 𝑠 + 𝑤
• K: Gain of the system
• 𝜁: Damping ratio of a system
• 𝑤 : Natural Frequency
• The time response of unit step input is
𝐾𝑤
𝐶 𝑠 =
𝑠(𝑠 + 2𝜁𝑤 𝑠 + 𝑤 )
After inverse Laplace transformation:
𝜁
𝑐 𝑡 = 𝐾[1 − 𝑒 (cos 𝑤 𝑡 + sin(𝑤 𝑡) ]
1−𝜁
𝑤 =𝑤 1−𝜁
• 𝑤 : Damped natural frequency
Slide 13
B) Time Response of second order systems
• The time response of unit step input is
𝐾𝑤
𝐶 𝑠 =
𝑠(𝑠 + 2𝜁𝑤 𝑠 + 𝑤 )
After inverse Laplace transformation:
𝜁
𝑐 𝑡 = 𝐾[1 − 𝑒 (cos 𝑤 𝑡 + sin(𝑤 𝑡) ]
1−𝜁
𝑤 =𝑤 1−𝜁
• 𝑤 : Damped natural frequency
Or in a compact form:
𝑒
𝑐 𝑡 = 𝐾[1 − sin(𝑤 𝑡 − Φ) ]
1−𝜁
Φ = atan(𝜁/ 1 − 𝜁 )
Slide 14
Pole-Zero Map ( )
𝐾𝑤
𝐺 𝑠 =
𝑠 + 2𝜁𝑤 𝑠 + 𝑤
Poles:
𝑠 , = −𝜁𝑤 ± 𝑗𝑤 1−𝜁
Damping ratio:
𝜁 = cos(𝜃)
Slide 15
Connecting the pole-zero map to responses
Slide 16
Time Domain Specifications
Slide 17
Time Domain Specifications
• Evaluation of Tp: Tp is found by differentiating c(t) and finding the first zero crossing after
t=0.
Slide 18
Time Domain Specifications
• Evaluation of % OS: Maximum overshoot is
Slide 19
Time Domain Specifications
For a step unit and 𝐾 = 1, 𝑐 = 1 and thus
Slide 20
Time Domain Specifications
• Evaluation of 𝑇𝒓 : Rise time is the time required for the system
to rise from 10% to 90% of its final value
0.8 + 2.5𝜁
𝑇 ≅
𝑤
Slide 21
Time Domain Specifications
• Evaluation of 𝑇 : we must find the time for which c(t) reaches
and stays within ±2% of the steady-state value cfinal.
Slide 22
Time Domain Specifications
• The settling time is the time it takes for the amplitude of
decaying sinusoid to reach 0.02.
𝒆 𝜻𝝎𝒏 𝒕
𝑐(𝑡) = 1 − cos 𝜔 1 − 𝜁 𝑡 − 𝜙
𝟏 − 𝜻𝟐
• Thus,
𝑒
= 0.02
1−𝜁
−ln 0.02 1 − 𝜁
𝑇 =
𝜁𝜔
• You can verify that the numerator varies from 3.91 to 4.74 as
varies from 0 to 0.9. Let us agree on an approximation for
the settling time:
4
𝑇 =
𝜁𝜔
Slide 23
Connecting the pole-zero map to response specifications
𝐾𝑤
𝐺 𝑠 =
𝑠 + 2𝜁𝑤 𝑠 + 𝑤
Poles:
𝑠 , = −𝜁𝑤 ± 𝑗𝑤 1−𝜁
Response specifications:
Slide 24
Connecting the pole-zero map to response specifications
Lines of constant
• Peak time,𝑇𝑝
• Settling time,𝑇𝑠 ,
• Percent overshoot, %𝑂𝑆
𝑇 <𝑇
𝑇 <𝑇
%𝑂𝑆 < %𝑂𝑆
Slide 25
Connecting the pole-zero map to response specifications
Step responses of second-order underdamped
systems as poles move:
a. with constant real part;
b. with constant imaginary part;
c. with constant damping ratio
Slide 26
Connecting the pole-zero map to response specifications: Example-1
Given the pole plot shown in Figure, find , 𝑤 ,
𝑇 , % 𝑂𝑆 and 𝑇
Slide 27
C) Time Response of high-order/complex LTI systems
• Let us handle the following system
𝐾
𝐺 𝑠 = ,0 ≤ 𝜁 < 1
𝑠(𝜏 𝑠 + 1)(𝜏 𝑠 + 2𝜁𝜏 𝑠 + 1)
• For any input, the system response will include
• A constant term from the 𝑠 factor
• An 𝑒 / term from (𝜏 𝑠 + 1)
• The following two terms from (𝜏 𝑠 + 2𝜁𝜏 𝑠 + 1)
/
𝑡 1−𝜉
𝑒 sin
𝜏
/
𝑡 1−𝜉
𝑒 cos
𝜏
• Not easy to comment on the system response!
• It would be really nice if we could extract the dominant dynamics of the system to
comment on the system response
𝜏 𝑠+1
𝐺(𝑠) = 𝐾
𝜏 𝑠+1 𝜏 𝑠+1
𝐾 𝜏 +𝐾 𝜏
𝐾 =𝐾 +𝐾 ; 𝜏 =
𝐾 +𝐾
• However, they are also quite useful in stabilizing closed loop control systems!
Slide 44
C) Time Response of high-order/complex LTI systems: extra zero in the LHP
Slide 45
C) Time Response of high-order/complex LTI systems: extra zero in the RHP
Slide 46
C) Time Response of high-order/complex LTI systems: Time Delay
• Time Delay: delay between the commanded response and the start of the output response
𝐺 𝑠 =𝑒 𝐺 𝑠
• This is NOT a rational transfer function!
• We know that it easier to work with polynomials
• Extracting poles and zeros
• Root-Locus, Routh-Hurwitz
𝐿 𝐿𝑠
𝑒 / 1− 𝑠+
≈ 2 12
𝑒 / 𝐿 𝐿𝑠
1+ 𝑠+
2 12
Slide 48
C) Time Response of high-order/complex LTI systems: Time Delays
• Pole-Zero Maps
• 1st Pade approximation
𝐿
1− 𝑠
𝑒 ≈ 2
𝐿 −2 2
1 + 2𝑠 𝐿 𝐿
Slide 49
4.3 Steady State Responses and Errors
• The steady state error is a measure of the accuracy of the
control system
• In general, the steady state errors of LTI control systems
depend on the “type of the input” and “type of the system”.
Slide 50
4.3 Steady State Responses and Errors
• To find 𝐸(𝑠), the error between the input, 𝑅(𝑠), and the output, 𝐶(𝑠), we write
Slide 51
4.3 Steady State Responses and Errors
• To find 𝐸(𝑠), the error between the input, 𝑅(𝑠), and the output, 𝐶(𝑠), we write
Slide 52
4.3 Steady State Responses and Errors
Note that 𝑠𝐸(𝑠) must have not any poles on 𝑗𝑤 axis or in the RHP of the s-domain (i.e. Stability).
Slide 53
4.3 Steady State Responses and Errors
Slide 54
A) Steady State Errors for Step Function Inputs
• Step Input with magnitude “1”, R(s)=1/s
Slide 55
B) Steady State Errors for Ramp Function Inputs
• Ramp Input with magnitude “1”, R(s)=1/s2
Slide 56
C) Steady State Errors for Parabolic Function Inputs
• Parabolic Input with magnitude “1”, R(s)=1/s3
Slide 57
Example
Find steady-state errors for inputs of 5 u(t), 5t u(t) , and 5t2 u(t) to the system shown in Figure. The
function u(t) is the unit step.
Slide 58
The Type of the Control System
• Since steady-state errors are dependent upon the number of integration in the forward path, we
give a name to this system attribute.
• We define system type to be value of 𝑛 in the denominator, or equivalently, the number of pure
integration in the forward path.
Slide 59
The Type of the Control System
• We define system type to be value of 𝑛 in the denominator, or equivalently, the number of pure
integration in the forward path.
Slide 60
Example: Steady State Specifications
Given the control system in Figure, find the value of K so that there is 10% error in the steady-state.
Since the system is Type 1, the error stated in the problem must apply to a ramp input; only a ramp
yields a finite error in a Type 1 system. Thus,
Slide 61
D) Steady State Errors for Disturbances
• Feedback control systems are used to compensate for disturbances or unwanted inputs that
enter a system
Slide 62
D) Steady State Errors for Disturbances
• eR(), is the steady-state error due to R(s), which have already obtained.
• eD(), is the steady-state error due to the disturbance.
Slide 63
D) Steady State Errors for Disturbances
• Let us explore the conditions on eD() that must exist to reduce to error due to the disturbance.
• We must make an assumption on D(s) now
• Step disturbance, D(s) = 1 / s.
• Then,
Slide 64
Example:
Find the steady-state error due to a step reference signal and disturbance for the system
Slide 65
E) Steady State Errors for Nonunity Feedback Systems
• A general feedback system, showing the input transducer, G1(s), controller and plant, G2(s), and
feedback, H1(s), is shown in figure
• To employ the presented approach, we must forming an equivalent unity feedback system from
a general nonunity feedback system
Slide 66
E) Steady State Errors for Nonunity Feedback Systems
• To employ the presented approach, we must forming an equivalent unity feedback system from
a general nonunity feedback system
Slide 67
Example:
Find the system type, the appropriate error constant associated with the system type, and steady-
state error for a unit step input.
Thus, the system is Type 0, since there are no pure integrations. The appropriate static error
constant is then Kp, whose value is
Slide 68