Control Systems Chapter 07 - FME - GIKI
Control Systems Chapter 07 - FME - GIKI
• Next, we will determine the dynamic response of linear systems using analytical techniques (i.e., the
total solution of linear ODEs)
– In particular, we will investigate the analytical system response to “standard” inputs such as the step, ramp,
and impulse
• We will focus on solving “simple” systems (1st- and 2nd-order ODEs) “by hand” using analytic methods
• System designers should understand the basic parameters governing 1st- and 2nd-order
system response
– How does the “resistance” parameter effect time to reach the steady-state response?
• Often the response of a complex, nonlinear model “looks like” a standard 1 st- or 2nd-order
response, and therefore the complex system can be approximated and replaced by a
reduced-order model
– System designers can trade-off model complexity vs. model accuracy, as well as obtain a physical,
intuitive feel for dynamic system response
Complete solution
a2 y a1 y a0 y 0 y P (t )
Related to input u(t)
y H (t )
Complete solution
7.2 Analytical Solutions to
Linear Differential Equations
• As an example, consider the 2nd-order I/O equation (ODE)
2 y 12 y 16 y 3u (t )
• The total solution y(t) is comprised of two parts: the homogeneous (or natural) solution
yH(t) and the particular (or forced) solution yP(t)
y (t ) y H (t ) y P (t ) Total solution
• Homogeneous solution: the solution to the homogeneous ODE, where the input u(t) is
zero. Therefore, the homogeneous solution is the natural response, and is strictly due to
the “natural dynamics”.
Homogeneous Solution
• Homogeneous solution: the solution to the homogeneous ODE, where the input u(t) is
zero
• Therefore, solve the unforced ODE with zero input:
Complete solution
2 y 12 y 16 y 0
y (t ) y H (t ) y P (t )
• Assume an exponential solution form for yH(t)
y H (t ) ce rt
• Therefore, y H cre rt and y H cr 2 e rt
2 rt rt
2cr e 12cre 16ce 0 rt
2r 2
12r 16ce rt 0
Because cert cannot be zero for all time t, the bracket term must = 0
Characteristic Equation
• The Characteristic Equation is the bracket term set to zero:
2r 2 12r 16 0
• The characteristic roots (r1 = -2 and r2 = -4 in this case) give us the complimentary
solution:
The constants c1 and c2 are determined from the two initial conditions (ICs) after the
particular response has been determined
4 y 8 y 6
2
y 8 y 6 y 0 For Homogenous solution with zero input.
yP a sin 4t b cos 4t
Need to find the constant a and b
2( 16a sin 4t 16b cos 4t ) 8(4a cos 4t 4b sin 4t ) 6(a sin 4t b cos 4t ) 10 sin 4t
Solving an LTI System: Example 7.2
( 32b 32a 6b) cos 4t ( 32a 32b 6a) sin 4t 0 cos 4t 10 sin 4t
32b 32a 6b 0
Roots (a 0.1529, b 0.1882)
32a 32b 6a 10
yP 0.153sin 4t 0.188cos 4t
t 3t y(0) = 2
y(t ) c1e c2 e 0.15sin 4t 0.188cos 4t Complete solution
y (0) 1
y( 0 ) c1e 0 c2 e 3( 0) 0.15sin 4(0) 0.188cos 4(0) Initial conditions
2 c1 c2 0.188 A
• Steady-state response: part of the total response that remains as time t infinity
Transient Steady-state
response response
• The values of s that make the denominator polynomial a(s) equal zero
are called the poles of the transfer function.
• Example 7.3: 3-way valve with applied force f(t) 0.04 y 16 y 7000 y f (t )
1
Transfer function: G (s)
0.04s 2 16s 7000
r1 200 j 367.42
Poles of TF G(s) are the same as
r2 200 j 367.42 the characteristic roots
DC Gain and the Steady State
• An extremely useful analysis technique is the use of the DC gain to compute a system’s
steady-state response to a constant input. The term “DC gain” is derived from circuit
analysis (DC implies a constant, non-oscillating input)
• Definition: the system DC gain is the steady-state gain to a constant input, for the case
when the output has a constant value at steady state.
fSS = 16 N
• Therefore, at steady state, the solenoid dynamics can be replaced by a constant gain = 8,
and therefore the steady-state EM force is f ss 2 8 16N
DC Gain: Example 7.4 (2)
•The DC gain of the valve TF is simply 1/ 7000 1.42 10 4
12 1
Of course, the DC gain of the entire system =
5 7000
• We can easily find the natural response by writing the characteristic equation for the
standard 1st-order model:
For the case with ZERO input, the constant c is simply the initial condition (IC)
for the dynamic variable, c = y(0) = y0
First-Order System: Time Constant
• Repeat the natural response (for zero input)
yH(t) = y0e-t/
y y 0
Decayed to 2% e 4 0.018
y y bu (t ) y y Ab
y (t ) y H (t ) y P (t ) Total response
yH = ce-t/ yP Ab yss
– Sketch an exponential response from the initial condition y0 to the steady-state value yss
– The total response will “decay down” if y0 > yss , or show an “exponential rise” if yss > y0
1st-Order System Step Response
Example 7.5
• The electrical input is voltage ein(t), and the output is force f (in N). An appropriate first-
order model is
the time constant is = 0.002 sec, and the settling time tS = 4 = 0.008 sec = 8 ms
• The step response (from Simulink) is below. The response is clearly exponential and
the settling time is indeed 8 ms, and the steady-state response is 16 N
2 y 10 y 3u (t )
y a1 y a0 y b0u (t )
• The natural or free response (zero input) will depend on the roots of the characteristic
equation
a1 a12 4a 0 Recall: j 1
r
2
1) real and distinct; 2) real and repeated (radicand = 0);
3) complex conjugates (radicand < 0); 4) purely complex (coefficient a1 = 0 and a0 > 0)
Case1:
2nd-Order Root Locations Case 1 real and distinct
Example: y 4 y 3 y u (t )
Therefore, the natural response (no input) is composed of two exponential functions, with time
constants 1 sec and 1/3 sec
r1 < 0 , r2 < 0
Case1:
2nd-Order Root Locations Case 1 real and distinct
Example: y 2 y 3 y u (t )
Characteristic Equation r 2 2r 3 0
y H (t ) c1e t c2 e 3t
Characteristic roots r1 = 1, r2 = -3
unstable
Therefore, the natural response (no input) is composed of two exponential functions; but et
diverges to infinity (unstable)
Case1:
2nd-Order Roots: Case 1 real and distinct
r1 > 0 , r2 < 0
Case1:
2nd-Order Root Locations Case 1 real and distinct
• Both roots are real, one root is zero, one root is negative
Example: y 3 y u (t )
Therefore, the natural response (no input) is composed of a constant and one exponential function
(marginally stable)
Case1:
2nd-Order Root Locations Case 1 real and distinct
r1 = 0 , r2 < 0
Case2
2 -Order Root Locations Case 2
nd
real and repeated
• Two real repeated roots (both are negative for this example)
Example: y 6 y 9 y u (t )
Characteristic Equation r 2 6r 9 0
Characteristic roots r1 = -3, r2 = -3
y H (t ) c1te 3t c2 e 3t
Case2
real and repeated
2nd-Order Root Locations Case 2
r1 = r2 < 0
Case2
real and repeated
2nd-Order Root Locations Case 2
r1 = r2 > 0
Case2
2nd-Order Root Locations Case 2 real and repeated
r1 = r2 = 0
Case 3:
2nd-Order Root Locations Case 3 complex conjugates
Example: y 4 y 29 y u (t )
Out put
Characteristic Equation r 2 4r 29 0
Characteristic roots r1, 2 2 j 5
Therefore, the natural response (no input) is a stable, damped sinusoidal function, where the real
part is the exp decay, and imaginary part is the frequency
2nd-Order Root Locations Case 3 Case 3:
complex conjugates
Example: y 4 y 29 y u (t )
Characteristic Equation r 2 4r 29 0
Characteristic roots r1, 2 2 j 5
y H (t ) c1e 2 j 5 t c2 e 2 j 5 t e 2t c1e j 5t c2 e j 5t
y H (t ) e 2t c3 cos 5t c4 sin 5t
unstable
2nd-Order Root Locations Case 3:
complex conjugates
Case 3 (con’t)
Complex conjugate roots with POSITIVE real part
Example: y 9 y u (t )
Characteristic Equation r 2 9 0
Characteristic roots r1, 2 j 3
y H (t ) c1e j 3t c2 e j 3t
r1, 2 j
Damping Ratio and
Undamped Natural Frequency
• We can characterize the transient response of a second-order system by two basic
parameters
– Damping ratio
– Undamped natural frequency, n
a1 a12 4a 0
Root locations: r
2
Critically Damped System
a1 a12 4a 0
r
2
If a12 < 4a0 then the roots If a12 > 4a0 then the roots
If a1 = 4a0 then the roots
2
• Therefore, the condition a12 = 4a0 is the critical transition point between a damped oscillating
response and a damped exponential response
a1
“Damping ratio”
2 a0
• If 0 < < 1, the system is underdamped (oscillations are present in the natural or
transient response)
• If > 1, the system is overdamped (no oscillations in the natural or transient response)
2nd-Order System: Standard Form
• Like the 1st-order system, we can define a “standard form” for the 2 nd-order system:
• For the undamped ()/underdamped (0<) case, the (complex) roots are
r n jn 1 2
y H (t ) ce nt cosd t
Underdamped 2nd-Order System
• The natural (free) response is (repeated)
• Note that the real part of the complex roots determines the exponential decay rate of the transient
response
• Note that the imaginary part of the complex roots determines the frequency of oscillation of the
transient response
“Damped frequency”
d n 1 2
frequency of oscillation of
the transient response
Underdamped System: Root Locations
• Root locations in the “complex plane”, where +x is the Real-number axis, and +y is the Imaginary-number axis
Imaginary
n , n 1 2
r n jn 1 2 X
2nd-order underdamped Complex Plane
root locations b c
θ
Right triangle: a2 + b2 = c2 a 0
Real
a = n Real part
0 Real
X
X
Circles:
All roots have
the same undamped X Complex roots ALWAYS appear
natural frequency ωn, increasing as conjugate pairs symmetric
radius increases ωn about the Real axis
Step Response of an Underdamped 2nd Order System
• Performance criteria for the transient response of an underdamped second-order system are
often based on the step response
y t 1 e cos t sin t
t Unit Step response of an underdamped
system
n ; n 1 2
y t 1 e cos t sin t
t
/ d
ymax 1 e cos sin Put n , d n 1 2
ymax 1 e / 1 2
d
Function of ζ only
Step Response of an Underdamped 2nd Order System
/ 1 2
M os e
ymax yss 1 M os
Step Response of an Underdamped 2nd Order System
2 1 2
N cycles Function of ζ only
Step Response of an Underdamped 2nd Order System
The concepts and formulation discussed so far are for underdamped system only
Therefore, first determine if the 2nd-order system is underdamped (check roots or ζ ) before applying these
transient-response equations!
2nd order system
k
G ( s)
Ymax s 2 2n s n2
Peak time tp
n 1 2
Yss
4
Settling time ts
n
Maximum overshot
ymax yss
M os
yss
tr t p ts / 1 2
Maximum overshot M os e
2nd-Order Underdamped System:
Example 7.8
• Sketch the step response
/ 1 2
Max overshoot: M OS e = 0.498 (48.9% overshoot)
4
Settling time: tS = 1.0 sec
n
2
Period: Tperiod = 0.3574 sec Ncycles = 2.8
n 1 2
2nd-Order Underdamped System:
Example 7.8 (3)
48.9% overshoot
Peak response max = 0.0573 rad
Peak time tp = 0.179 sec
x Ax
• Following our previous methods for solving ODEs, let’s assume the solutions are
exponential functions of time:
x1 (t ) c1e t
x2 (t ) c2e t
or, in compact
x(t ) ce t
notation:
t
xn (t ) cn e
SSR and Eigenvalues (2)
• The time-derivative of the assumed solution is
x (t ) ce t
• Sub this expression into the state equation
x Ax ce Acet t
• Move all terms to the left-hand side and factor out cet
I Ace t
0
I A detI A 0
• The n values of that satisfy the determinant equation are called the
eigenvalues of matrix A
• The characteristic roots are computed from the characteristic equation, which can be
derived from the system’s I/O equation, transfer function, or state-space
representation
– The location of the roots in the complex plane determines the system’s response speed and
damping characteristics