Control Systems - Lecture Note
Control Systems - Lecture Note
Engineering
Eeng 3182
x1 (t ) y1 (t ) x2 (t ) y2 (t )
x1 (t ) + x2 (t ) y1 (t ) + y2 (t )
Example:
x(t ) y (t )
y (t ) = x(t ) x(t − 1)
let x(t ) = x1 (t )
y1 (t ) = x1 (t ) x1 (t − 1)
let x(t ) = ax1 (t )
y (t ) = ax1 (t )ax1 (t − 1) = a x1 (t ) x1 (t − 1) = a y1 (t )
2 2
y (t ) ≠ ay1 (t )
Non linear system
37 Bahirdar University, Institute of Technology_IOT_SCEE 4/8/2014
Time Invariant:
A system is said to be time invariant if a time
delay or time advance of the input signal leads
to an identical time shift in the output signal.
x(t ) y (t )
Time invariant
system
x(t − t0 ) y (t − t0 )
t0 t0
38 Bahirdar University, Institute of Technology_IOT_SCEE 4/8/2014
Example: y (t )
x(t ) y (t ) =
x(t )
R (t )
x1 (t )
y1 (t ) =
R (t )
x2 (t ) = x1 (t − t0 )
x2 (t ) x1 (t − t0 )
⇒ y2 (t ) = =
R(t ) R(t )
x1 (t − t0 )
but y1 (t ) =
R (t − t0 )
y1 (t − t0 ) ≠ y2 (t ), for t0 ≠ 0
U (s ) 1 Y (s )
LCs 2 + RCs + 1
Block diagram
42 Bahirdar University, Institute of Technology_IOT_SCEE 4/8/2014
Impulse Response:
Output signals from the system due to Input signals and
initial conditions.
1. The point of view of Mathematic:
Homogenous solution yh (t ) + Particular solution y p (t )
Where;
Laplace transform
of basic signals
Solution
Taking the Laplace transform of the differential at
zero initial conditions
R(s) C(s)
Block diagram
53 Bahirdar University, Institute of Technology_IOT_SCEE 4/8/2014
Using transfer function, we can analyze the system
response.
Taking the above example, calculate the System Response
C(t), if the input be unit step r(t)=U(t), which is step
response.
Laplace of r(t), R(s)=1/s, then
Expanding
Solution:
Using KVL around the loop
Substituting current:
For Resistor:
For Inductor:
Solving V/I
Solution
V 2( s ) −1
V 1( s ) RCs
V 2( s )
−RCs
V 1( s )
Torques on J1
due only to the
motion of J2
Final free-body
diagram for J1
Where;
Substituting ϴ2 with ϴ1
After simplification
Which gives
representation in
the form of part b
78 Bahirdar University, Institute of Technology_IOT_SCEE 4/8/2014
Example: Find the transfer function for the following
system taking the angular displacement of the second
inertia as an output when we apply an input torque on the
first inertia:
Voltage-Speed: e b = K bω m
where Kt: torque constant, Kb: velocity constant For an
ideal motor.
Combing previous equations results in the following
mathematical model:
di a
La + Ra i a + K bω m = ea
dt
J ω& m + Dωm - K t ia = 0
Taking Laplace transform of the system’s differential
equations with zero initial conditions gives:
83 Bahirdar University, Institute of Technology_IOT_SCEE 4/8/2014
... Ra La
(La s + Ra )I a ( s ) + K b Ω m ( s ) = Ea ( s ) B
ia
(Js + D )Ω m (s) - K t I a ( s ) = 0
Kt
ea
ω
Liquid Tank
Flow rate = Rate of change of volume
Laminar Flow:
(Linear)
Turbulent Flow:
(Non-linear)
q ICs + RCs + 1 2
i
Pickoff Point
Closed Loop
Transfer Function
Solution:
s3 +1
Cascade/Series: S² S4 + S2 + S
S3 +1
Parallel: s3 +1
s3 +1
s
Feedback: S3 + S +1
2S 3 + S + 2
108 Bahirdar University, Institute of Technology_IOT_SCEE 4/8/2014
R(s) − +1 C (s)
Cascade/Series
−
3
S 1
S S +
Parallel
S
R(s ) +1 1 C (s )
−
3
S
S 3
+ S +1 S
Feedback
S
R (s ) S 3 +1 C (s)
− S 4 + S 2 + s
Cascade/series
S
R(s) C (s )
S3 +1
2S 3 + S + 2
Feedback
Hence; s3 +1
G (s) = 3
2S + S + 2
NB:
All variables, summing points and takeoff points are also
represented by nodes.
The negative & positive sign on the summing junction is
assigned on the system transfer function.
112 Bahirdar University, Institute of Technology_IOT_SCEE 4/8/2014
Block diagrams Signal flow graphs
a. Block diagram
L3= G7 H7
Non-touching loops are
L1 & L2, L1 & L3, L2 &L3
Signal flow graph
Self loop: L2H3
Forward Path: direct path from R(s) to C(s)
Forward path gain 1: G1G2G3G4
Forward path gain 2: G5G6G7G8
118 Bahirdar University, Institute of Technology_IOT_SCEE 4/8/2014
Example 2:
Where;
k= number of forward paths
Tk= the kth forward-path gain
Calculating
Where;
Thus;
Settling Time, Ts
Settling time is defined as the time for the response to
reach, and stay within, 2% of its final value.
Letting c(t) = 0.98 and solving for time, t, we find the
settling time to be
and,
C(s) ω n2
= 2
R( s ) s + 2ζω n s + ω n2
ωn un-damped natural frequency of the second order system, which
is the frequency of oscillation of the system without damping.
ζ damping ratio of the second order system, which is a measure of
the degree of resistance to change in the system output.
− ω nζ ± ω n ζ 2 − 1
1
T =
ζω n Real Part Imaginary Part
NB: The damping ratio and the natural frequency of a system can be
related with the pole location of the system : which means if the
information about these two parameter is known we can calculate
the poles of the system.
Solution:
0.5 0.38 0.28 0.17 0.08 2.5
0.64
2
2
C( s ) ω n2 4 0.8
1.5
= 2 = 2 1
R( s ) s + 2ζω n s + ω n
1
2
s + 2s + 4 0.94 0.5
0.94 0.5
− ωnζ + ωn ζ 2 − 1 = −1 + j1.73 -1
0.8
1
1.5
− ωnζ − ωn ζ 2 − 1 = −1 − j1.73
-2
2
0.64
0.5 0.38 0.28 0.17 0.08 2.5
-3
-2 -1.5 -1 -0.5 30
149 Bahirdar University, Institute of Technology_IOT_SCEE 4/8/2014
Second order response as a function of damping ratio
1 s + 2ζω n
C( s ) = −
s (s + ζω n )2 + ω n2 1 − ζ 2 ( )
1 s + 2ζω n
C( s ) = −
s (s + ζω n )2 + ω d2
Where ωd=ωn√(1-ζ), is the frequency of transient oscillations and is called
damped natural frequency.
152 Bahirdar University, Institute of Technology_IOT_SCEE 4/8/2014
1 s + ζω n ζω n
C( s ) = − −
s (s + ζω n ) + ω d (s + ζω n )2 + ω d2
2 2
1 s + ζω n ζω n
C( s ) = − −
s (s + ζω n ) + ω d (s + ζω n )2 + ω d2
2 2
ζ
ω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 −ζωn t sin ω d t
1−ζ 2
153 Bahirdar University, Institute of Technology_IOT_SCEE 4/8/2014
−ζωn t ζ
c(t ) = 1 − e cos ω d t − e −ζωnt sin ω d t
1−ζ 2
ζ
c(t ) = 1 − e −ζωn t cos ω d t + sin ω d t
1 − ζ 2
where;
1.6
1.2 1.2
1.4
1 1
1.2
1 0.8 0.8
0.8 0.6
0.6
0.6
0.4 0.4
0.4
0.2 0.2
0.2
0 0 0
0 2 4 6 8 10 0 2 4 6 8 10 0 2 4 6 8 10
1.2
0.8
0.6
0.4
0.2
0
0 1 2 3 4 5 6 7 8 9 10
We can see the relationship between the value of ϛ and the type of
time response obtained: The lower the value of ϛ, the more
oscillatory the response.
The natural frequency is a time-axis scale factor and does not
affect the nature of the response other than to scale it in time.
157 Bahirdar University, Institute of Technology_IOT_SCEE 4/8/2014
Time response for underdamped response for different
natural frequency.
1.4
1.2
0.8
wn=0.5
0.6 wn=1
wn=1.5
wn=2
0.4 wn=2.5
0.2
0
0 1 2 3 4 5 6 7 8 9 10
158 Bahirdar University, Institute of Technology_IOT_SCEE 4/8/2014
Other parameters associated with the
underdamped response
Rise time, Tr: The time required for the waveform to go from
0% or 10% of the final value to 90% or 100% of the final value.
Peak time, Tp: The time required to reach the first the
maximum, peak.
Percent overshoot, %OS: The amount that the
waveform overshoots the steady-state, or final value at the
peak time, expressed as a percentage of the steady-state
value.
Settling time, Ts: The time required for the transient's
damped oscillations to reach and stay within ±2% ( ± 5%)
of the steady-state value.
159 Bahirdar University, Institute of Technology_IOT_SCEE 4/8/2014
For 0< ζ <1 and ωn > 0, the 2nd order system’s response due to a unit step
Therefore:
Hence;
162 Bahirdar University, Institute of Technology_IOT_SCEE 4/8/2014
The primary peak occurs at n=1;
How ϛ calculated
For a unit step input, from %OS ?
Then;
163 Bahirdar University, Institute of Technology_IOT_SCEE 4/8/2014
Settling time, Ts
The time for which C(t) in reaches and stays within
±2% of the steady-state value, Cfinal.
Using this definition, the settling time is the time it takes
for the amplitude of the decaying sinusoid reaches 0.02.
= 2.838%
= 0.533secs
The normalized rise time (Trωn) for 0.75 damping ratio is
around 2.3 seconds.
Tr=2.3/ωn=0.23secs
166 Bahirdar University, Institute of Technology_IOT_SCEE 4/8/2014
Relationship b/n %OS, Tp and Ts to the location of
the poles of the system
The pole plot for a general, underdamped second-order
system, From the Pythagorean theorem
that the radial distance from the
origin to the pole is the natural
frequency, ωn, and the cosine of
the angle formed b/n this line &
the real axis be the damping ratio.
σd
ς = cosθ =
ωn
Where;
ωd is the imaginary part of the pole and is
called the damped frequency of oscillation,
σd is the magnitude of the real part of the pole
and is the exponential damping frequency.
168 Bahirdar University, Institute of Technology_IOT_SCEE 4/8/2014
From the equations, Tp is inversely proportional to the
imaginary part of the pole. Since horizontal lines on the s-
plane are lines of constant imaginary value, they are also
lines of constant peak time.
The settling time is inversely proportional to the real
part of the pole.
Since vertical lines on the s-plane are lines of constant
real value, they are also lines of constant settling time.
Finally, since ϛ = cosѳ, radial lines are lines of constant
ϛ.
Since percent overshoot is only a function of ϛ, radial
lines are thus lines of constant percent overshoot, %OS.
169 Bahirdar University, Institute of Technology_IOT_SCEE 4/8/2014
Lines of constant peak time, Tp, settling time, Ts, and
percent overshoot, %OS.
• Stepresponses of second-
order underdamped systems
as poles move with constant
damping ratio.
171 Bahirdar University, Institute of Technology_IOT_SCEE 4/8/2014
Example:
For the system shown in the Figure, determine the values
of gain K and velocity-feedback constant Kh so that the
maximum overshoot in the unit-step response is 0.2 and
the peak time is 1 sec. With these values of K and Kh,
obtain the rise time and settling time. Assume that J=1 kg-
m2 and B=1 N-m/rad/sec.
ω n = 3.53
3.53 = K
π
tr =
ωn 1 − ζ 2
3.532 = K
t r = 0.65s
K = 12.
12.5 4
ts =
(1 + KK h ) ζω n
ζ =
2 K
Then;
Solution:
Solution:
1 5
R(s) = s and T(s)= s 2 + 7 s + 10