3 Time Response Stability and Steady State Error
3 Time Response Stability and Steady State Error
Output response = forced response (e.g. constant) + natural response (e.g. exponential)
roots of the denominator (characteristic
• Poles of a Transfer Function (TF): The values of s that cause polynomial) of the transfer function
• Zeros of a TF: the values of s that cause TF = 0. roots of the numerator of the transfer function
where and
System showing input and output
Input function
generates
Input poles Forced response
( pole at the origin generated a step function at the output)
generates
System poles Natural response
Transfer function
Problem:
Given the following system, write the output, c(t), in general terms.
Specify the forced and natural parts of the solution.
Solution:
Forced Natural
response response
Taking inverse Laplace transform,
When
1. Time constant : Time it takes for the step response to rise to 63% of its final value. 1
𝑇 𝑐=
𝑎
• we can call the parameter (system pole) the exponential frequency (The reciprocal of the time constant)
• is related to the speed at which the system responds to a step input.
2. Rise Time : Rise time is defined as the time for the response to go from 0.1 to 0.9 of its final value.
found by solving for the difference in time at c() = 0.9 and c() = 0.1
𝑡𝑖𝑚𝑒 𝑓𝑜𝑟 𝐶(𝑡)=0.1
ln ( 0.1 ) 2.31
𝐶 ( 𝑡 2 )= 0.9=1 − 𝑒
− 𝑎𝑡 2
𝑡 2=− = 2.13 0.11 2.2 2.2
𝑎 𝑎 𝑇 𝑟 =𝑡 2 −𝑡 1 = − = 𝑇 =
−𝑎 𝑡1 ln ( 0.9) 0.11 𝑎 𝑎 𝑎 Rise time: 𝑟
𝑎
𝐶 ( 𝑡 1 ) = 0.1=1 − 𝑒 𝑡 1 =− = 𝑡𝑖𝑚𝑒 𝑓𝑜𝑟 𝐶(𝑡)=0.9
𝑎 𝑎
3. Settling time : The time for the response to reach, and stay within, 2% of its final value.
Lettingand solving for time, , we find the settling time to be
ln ( 0. 98 ) 4 4
𝐶 ( 𝑇 𝑠 ) =0.98=1 − 𝑒
− 𝑎 𝑇𝑠 𝑇 𝑠 =− = 𝑇 𝑠=
𝑎
𝑎 𝑎
1
𝑎= =7.7
0.13
Natural frequency
( 𝛏=
exponential decay
𝑛𝑎𝑡𝑢𝑟𝑎𝑙 𝑓𝑟𝑒𝑞𝑢𝑒𝑛𝑐𝑦 )
: (Natural Frequency ) the frequency of oscillation of the system without damping.
(D) dimensionless measure describing how oscillations in a system decay.
exponential decay , real part of the pole
Poles
𝑏=𝝎 𝒏 2
𝑎=2 ξ 𝜔 𝑛
Canonical form
CEN455: Dr. Nassim Ammour 9
(two finite poles and no zeros)
Second-Order System2
The sign of the discriminant of the denominator polynomial depends on the damping ratio , three cases.
𝑎 √ 𝑎 2
−4𝑏 𝑎=2 ξ 𝜔 𝑛
System poles 𝑠 =− ∓
1, 2
2 2 2
𝑏=𝝎 𝒏
Case1: Overdamped system: () Two real poles
𝑎=9 𝑎𝑛𝑑 𝑏=9 the general case
(two finite poles and no zeros)
1 0.171 1 .171
𝐶 ( 𝑠 )= + −
𝑠 𝑠 +7.854 𝑠 +1.146
− 7.854 𝑡 − 1.146𝑡
𝑐 (𝑡 )=1+ 0.171 𝑒 −0 1 .171 𝑒
Overdamped system
𝑎=2𝑎𝑛𝑑𝑏=9
𝑎=6𝑎𝑛𝑑𝑏=9
𝑎=0𝑎𝑛𝑑𝑏=9 two system poles on the imaginary axis 𝑖𝑛𝑝𝑢𝑡 𝑝𝑜𝑙𝑒 𝑠𝑦𝑠𝑡𝑒𝑚𝑝𝑜𝑙𝑒𝑠
𝛏 >1
Under-damped responses
Two complex
Two poles at
0 < 𝛏 <1
Un-damped responses
Two imaginary
Two poles at 𝛏 =0
Example
For the system find the value of and report the kind of
response expected.
2
𝜔𝑛 12
We have 𝐺 ( 𝑠 ) = 2 2
= 2
𝑠 + 2 𝛏 𝜔𝑛 𝑠 + 𝜔𝑛 𝑠 + 8 𝑠 +12
𝜔 2𝑛 =12 𝜔 𝑛 =√ 12
and 2 𝛏 𝜔𝑛 =8 𝛏=
8
=
2
>1
2 √ 12 √ 3
System is over-damped.
Where,
CEN455: Dr. Nassim Ammour Second-order underdamped responses for damping ratio values 15
Underdamped Second-Order Systems
Specifications
• Other parameters associated with the underdamped response are rise time, peak time, percent overshoot, and settling time.
𝜎 𝑑 =𝛏 𝜔𝑛 𝑛 𝑑 𝑛 𝑛
Damping Ratio
damped frequency of oscillation,
𝜎 𝑑 𝛏 𝜔𝑛
cos ( 𝜃 )= = =𝛏 Natural frequency
𝜔𝑛 𝜔𝑛
inversely proportional to the
imaginary part of the pole.
Solution:
Damping ration,
2 2
Natural frequency, 𝜔 𝑛 =𝜔 𝑑 + 𝜎 𝑑
2
𝜔 𝑛 =√ 𝜔 +𝜎 =√ 7 +3 =7.616
2
𝑑
2
𝑑
2 2
Peak time,
• We need to approximate that system to a second-order system that has just two dominant complex poles.
21
CEN455: Dr. Nassim Ammour
Evaluating Pole-Zero Cancellation
Effect of a zero on the system: A system with a zero consists of the derivative of the original response and the scaled version of the original
response. If the zero is very large, the Laplace transform of the response is approximately the scaled version of the original response. As the
zero becomes smaller, the derivative term contributes more to the response and has a greater effect. ( 𝑠+ 𝑎 ) 𝐶 ( 𝑠 )=𝑠𝐶 ( 𝑠 ) +𝑎𝐶 (𝑠)
derivative scaled
pole-zero cancellation response response
Problem: For any function for which pole-zero cancellation is valid, find the approximate response.
Solution:
Stable system: If natural response approaches zero as time approaches infinity (LTI System).
Marginally stable system: If natural response neither decays nor grows but remains constant or
oscillates as time approaches infinity.
Stable systems have closed-loop transfer functions with poles only in the left-half plane.
How many poles are in left / right plane or in jw axis, not where
Routh Table Generation:
Denominator:
10
Example: T (s)
s 5 2 s 4 3s 3 6 s 2 5s 3
Example:
10
T (s)
s 5 2 s 4 3s 3 6 s 2 5s 3
Reverse coefficients:
D ( s ) 3s 5 5 s 4 6 s 3 3 s 2 2 s 1
Example: 10
T (s)
P( s) s 6s 8
4 2
s 5 7 s 4 6 s 3 42 s 2 8s 56
dP( s )
4 s 3 12s 0 Derivative of the polynomial of the row above the zeros row
ds
entire row
consists of zeros
Stable system.
CEN455: Dr. Nassim Ammour 28
Pole Distribution via Routh Table
PROBLEM:
with Row of Zeros
Tell how many poles are in the right half-plane, in the
left half-plane, and on the jw-axis.
20
T (s)
s s 12 s 22 s 39 s 4 59 s 3 48s 2 38s 20
8 7 6 5
P( s ) s 4 3s 2 2 :Even polynomial
Taking the derivative
interpretation
No sign change
PROBLEM: Find the range of gain, K, for the system that will cause the system to be stable, unstable, and
marginally stable. Assume K > 0.
Variable gain K
1. If K < 1386, then stable system.
2. If K > 1386, then two sign changes; two right-half plane
poles and one left-half plane pole. Unstable system.
3. If K = 1386, an entire row of zeros j poles.
K
T (s) replacing K=1386
s 3 18s 2 77 s K
dP( s )
P( s ) 18s 2 1386 36 s 0
ds
No sign change
can be positive,
zero, or negative
Row of
zeros
P( s ) s 2 10 (2)
PROBLEM: find out how many poles are in the left half-plane, in the right half-plane, and on the jw-axis.
0 3 1 10
x 2 8 1 x 0 u
10 5 2 0
y 1 0 0x
SOLUTION:
s 0 0 0 3 1 s 3 1
det( sI A) s 3 6 s 2 7 s 52
( sI A) 0 s 0 2 8 1 2 s 8 1
0 0 s 10 5 2 10 5 s 2
Using this polynomial, form the Routh table
step input
ramp input
Steady-state error
CEN455: Dr. Nassim Ammour 33
Evaluating Steady-State Errors1
s 2 7s 5
Applying final value theorem [ f () lim sF ( s )] E (s)
s0 s ( s 2 7 s 10)
T(s) is stable, hence E(s) is stable.
e() lim e(t ) lim sE ( s ) lim sR ( s )[1 T ( s )] Applying final value theorem,
t s 0 s 0
e( ) 1 / 2
CEN455: Dr. Nassim Ammour 34
Evaluating Steady-State Errors2
Problem: Find steady-state errors for inputs 5u(t), 5tu(t), and 5t2u(t) to the above system.
Solution: Let, the system is stable.
5 5 5
5u (t ) : e()
1 lim G ( s ) 1 20 21
s 0
5 5
5tu (t ) : e()
lim sG ( s ) 0
s 0
10 10
5t 2u (t ) : e( )
lim s 2G ( s ) 0
s 0
One integration, s1
5 5 1
5tu (t ) : e() No integration will make it infinity.
lim sG ( s ) 100 20 one integration makes it constant.
s 0
K p lim G ( s )
s 0
500 2 5 6
K v lim sG ( s ) 31.25
s 0 8 10 12
Position constant, Kp:
K a lim s 2G ( s ) 0
s 0
K p lim G ( s )
s 0
1
For step input, e() 0
1 K p
Velocity constant, Kv:
1 1
For ramp input, e() 0.032
K v lim sG ( s ) K v 31.25
s 0
1
Acceleration constant, Ka: For parabolic input, e()
Ka
K a lim s 2G ( s )
s 0
CEN455: Dr. Nassim Ammour 38
System Type
39
CEN455: Dr. Nassim Ammour
System Type
Relationships between input, system type, static error constants, and steady-state errors
Problem:
Find the steady-state error component
due to a step disturbance
PROBLEM:
Find the system type, error constant, and the steady-state error for a unit step input.
G ( s) 100( s 5)
Ge ( s ) 3
1 G ( s ) H ( s ) G ( s ) s 15s 2 50s 400