Solutions Manual to Accompany
Homework
For Linear Control
System Engineering
Liu Sijiu
HIT 651
Copyright 2006
P1-2 Redesign the turbine speed control system discussed in Sample
Problem 1.1, but replace the fly-ball governor with the tachometer
shown in Fig.1.2. A tachometer consists basically of a small DC
motor operated in reverse as a generator, the shaft being rotated
continuously, producing a DC voltage proportional to shaft speed.
Solution
2
P1-4. Shown in Fig.P1.4 is a water-level control system comprising a tank, inlet pipe, slide valve, and float.
Details of the operation of the flow valve are
also shown in the figure. Draw a block diagram
of the feedback control system and identify the
main elements, writing down the mathematical
transfer functions where appropriate.
Solution
P1-5 A method of producing a displacement
proportional to an input displacement but with a much
larger output force is shown in Fig.P1.5. The input
displacement x causes the movement of the spool valve
to produce a differential flow to the hydraulic actuator.
Draw a block diagram of the system, label the principal
parts of the control loop, and identify as many of the
transfer functions as possible. State any assumptions
made in the analysis.
Solution
3
P2-1 A system unknown transfer function is shown in Fig.P2.1. If a unit impulse applied at the input produces
at the output a signal described by the time function
, determine the unknown transfer function.
Solution
)t(c
t3e2
−
t3e2
−
=>
)t(c
=
=
=
C
R
2
3s
+
P2-2. Find the solution of the differential equation
+
2
xd
2
dt
)t(z
When
+
dx
dt
t2e
−=
x8
=
dz
dt
+
z3
and all other initial conditions are zero.
Solution
)s(Z
=
1
2s
+
⋅
3s
+
8s
++
2
s
=
1
2s
+
⎛
⎜⎜
⎝
−
Cs
+
2
2
8s
++
⎞
⋅⎟⎟
⎠
s
C
1
=
⎛
⎜
⎜
⎜
⎜
⎝
1
2s
+
−
11
2
+
s
−
1
2
)
31
4
s(
+
⋅
1
10
⎞
⎟
⎟
⎟
⎟
⎠
)t(y
=
e1.0
−
t2
−
e1.0
−
1
2
t
⎛
⎜
⎜
⎝
cos
31
2
t
−
23
31
sin
31
2
t
⎞
⎟
⎟
⎠
P2-3 For the system shown in Fig.P2.3., determine the relationship between voltage and current, express this
relationship in the form of a transfer function and determine the current as a function of time when the voltage
is a step change from zero to 10V.
Solution
1
sLR
+
+
=
/1
Cs
Cs
+
RCs
,
+
1
2
LCs
)s(Y
=
10
s
2
LCs
Cs
+
RCs
=
+
1
−
10
10
26
s
−
10
⋅
10
+
6
−
3
1s
+
=
2
s
+
10
3
s
10
6
+
10
)t(y
=
10
−
2
2
3
−
500
t
e
sin(
3
2
1000
)t
=
.0
01154
⋅
e
−
500
t
sin(
866
)t
P2-8 For the system shown in Fig2.8 determine the closed –loop
transfer function C/R.
Solution
C
R
=
G
HG1
+
1
1
⎛
1
⎜⎜
+⋅
⎝
G
G
2
1
⎞
=⎟⎟
⎠
+
GG
2
1
HG1
+
1
P2-9 For the single input system shown in Fig2.9, find the
transfer function of output to input C/R.
Solution
GG
1
2
2
3
HGG1
+
C
R
=
HGG1
+
23
2
+
GGG
21
+
2
3
HGG
121
=>
32
HGG1
+
2
GGG
21
3
HGGGHGG
121
33
21
+
4
P3.2 Determine the output of the open-loop system G(s)=a/(1+sT) to the input r(t)=t. Sketch
both input and output as function of time, and determine the steady-state error between the
input and output. Compare the result with that given by Fig.3.7.
Solution:
)s(R)s(c
=
⋅
e
−=−=
c
r
t
Tt(a
−+
a
)sT1(
+
T/t−
Te
As
∞→t
,
ess
−−=
at
t
aT
=
a
)sT1(
+⋅
)t(c
+−⋅=
Tt(a
Te
T/t−
)
for
for
1a
=
1a
≠
2
=
s
)
T
⎧
⎨
∞
⎩
P3.5 An open-loop first-order system is characterized b the transfer function
)S(G
=
1
s
τ+
1
, where the time
s5=τ
. Calculate the steady-state error when the system input is r(t)=1+6t. Confirm the result by
constant is
using the final-value theorem.
Solution: By the superposition theorem, the system output c(t) could be considered as a sum of a step
response and a ramp response. That is
−
5/t
e1(
−−
−
5/t
)
=
30
−
e5
−
5/t
, yields
(e
=∞
)
30
r
c
)t(e
=−=
e5t(6)1t6(
+−⋅−+
)
By the final value theorem for Laplace’s transform,
⎛
1)
−⋅
⎜
⎝
)s(G1)s(R)s(C)s(R)s(E
[
−⋅
1
s
−
=
+
=
1
1s5
+
⎞
=⎟
⎠
5)s6(
⋅+
)1s5(s
+
⋅
]
=
6(
2
s
5)s6(s
⋅+
⋅
)1s5(s
+
⋅
⋅
0
=
=
=
30
)t(e
Lim
s
∞→
)S(EsLim
s
→
Lim
t
∞→
P3.7 One definition of the bandwidth of a system is the frequency range over which the amplitude of the output
signal is greater than 70% of the input signal amplitude when a system is subjected to a harmonic input. Find a
relationship between the bandwidth and time constant of a first-order system. What is the phase angle at the
bandwidth frequency?
1
=
.0
707
≈
or
1
=τω+
22
b
2
2
2
Solution:
1
s
τ+
1
=
=ω
b
/1
τ
i.e.
1
22
τω+
b
tan
−=ω∠
)j(G
−
1
(
−=ωτ
)
tan
−
1
)1(
45
°−=
, phase lag is
°45
at the bandwidth frequency.
P3.8 Figure P3.8 shows the experimentally obtained voltage output of an unknown system subjected to a step
input of +10V. Determine the transfer function of the system and locate its pole on the complex plane.
Solution: System appears to be first order. So suppose the transfer function is as follow:
C)s(G
R
/t τ−−⋅
e1(K10
K
s
τ+
, then
)t(c
)
=
=
=
1
From final value theory
Lim
t
∞→
)t(c
=
)s(CsLim
s
→
0
⋅
=
Now consider time constant,
⋅
sLim
0
s
→
t
⋅=τ
10
s
K
1s
+τ
=
V5.2K10
=
⋅
, hence K=0.25
1
1
K10
)t(c
⎞
⎟
⎠
−
⎛ −
1ln
⎜
⎝
=1.4983
Matlab: K=0.25, c=[1.22, 1.84, 2.16, 2.33, 2.41], t=1:5, tao=t./(-log(ones(1,5)-C/10/K))*ones(5,1)/5
5
P4.1 Figure P4.1 shows a closed-loop feedback system with a
second-order plant. Determine the damped natural frequency and
damping ratio of the closed-loop response.
2
)s(C
ω
n
)s(R
2
ζω+
ζω
=
10
s6
+
16.3
2
s
ω+
n
n
2/6
=
G
GH
+
rad/s,
=ω
Hence
=>
s
=
=ζ
10
/3
=
+
=
=
3
1
s
2
2
n
n
ω=ω
d
n
1
=ζ−
10
⋅
9.01
−
=
1
rad/s
10
2
10
=
3.0
10
P4.4 Calculate the required value of gain K shown in follow such that the
closed-loop response of the system to a step input is limited to no more than
10% overshoot. Plot the closed-loop poles on the complex plane for the same
value of gain.
)s(C
)s(R
G
GH
+
K
sKs2
2
ω
n
ζω+
12/2
=
ζω
=
=
+
+
=
=
2
1
s
s
2
2
n
,
2
ω+
n
n
n
=
1
ζ
=ω
6.0=ζ
By textbook fig.4.17, 10% overshoot is responding to
5
3
4j1
3
The closed-loop poles will be
ω=ω
=ς−
±−=
25
9
25
9
=>
67.1
78.2
5
3
s 2
s2
K
,
,
=
=
=
=
+
+
1
0
s
n
d
2
6.01
−
,
K2
n =ω
2
=
4
3
=
33.1
or
s
ςω−=
n
±−=ω±
33.1j1
j
d
the close-loop poles on
P4.5 A unity-feedback control system has the forward-path transfer
function G(s)=K/s(s+10). Find the closed –loop transfer function, and
develop expression for the damping ratio and damped natural frequency in
terms of K. Plot
the complex plane for
K=0,10,25,50,100. For each value of K calculate the corresponding damping
ratio and damped natural frequency. What conclusions can you draw from
the plot?
)s(C
)s(R
K
sKs10
2
ω
n
ζω+
ω+
=
+
+
=
2
s
s
2
2
2
n
n
10
5
K
Hence
n =ω
K
and
2
ζω
n =
or
=ζ
Closed loop poles located at
s
ζω−=
1j
±
2
n
n
±−=ως−
25Kj5
−
P4.7 Prove that for an under damped second-order system subject to a step input, the percentage overshoot
above the steady-state output (as shown in textbook Fig.P4.7) is a function only of the damping ratio ζ.
By step response
1)t(c
−=
At the peak time,
ζω−
e
n
ζω−
1
t
n
1
2
ζ−
ζω−
t
n
e
sin(
)
φ+ω
t
d
sin(
ω+φ+ω
)
t
ζω−
t
n
e
d
cos(
=φ+ω
)
0
d
Tpd
t
d
π
ω
d
(tg
ω=φ+ω
)
t
d
/
ςω
n
d
φ=
=>
ω
=
,0
2,
ππ
...
=>
pT
=
ζω−
n
e
π
ω
d
PO
=
1
2
1
ζ−
sin(
ω
d
π
ω
d
=φ+
)
1
2
1
ζ−
e
ζπ−
2
1
ζ−
ζπ−
2
1
ζ−
sin
=φ
e
6
P5.2 A second-order system is given by the transfer function
)s(C
)s(R
=
10
9s2
+
+
=
10
9
2
2
9
9s2
+
s
+
s
Determine:
a. damping ratio
b. The damped and undamped natural frequency.
c. Maximum peak modulus Mp
d. The frequency at which Mp occurs.
e. The bandwidth, defined as the frequency range over which the modulus does not fall more
then 3db below the low-frequency (DC) value.
f. The steady-state output for a unit step input.
g. The location of the closed-loop poles on the complex plane.
h. The rise time.
i. The 2% settling time
Solution:
s/
9/10
ω=ω
3/1=ζ
n =ω
ζω−
19
=−
and
rad3
83.2
22
rad
ζω
s/
,
=
=
1
2
n
2
n
d
)
(
=
77.1
,
ω=ω
p
n
21
=ζ−
2
7
=2.65
n =
10
82
=
, so
5
4
2
=
M
p
=
2
⋅ς
1
ζ−
2
, where r=ω/ωn
22
r
x 2
−
1
3
=
++
x)9/4(1x2
=
29.4
3
07.77
+
=
2
=>
x 2
−
01x)9/14(
=−
,
=ω⋅=ω
r
29.4
b
n
From
.0M
=
707
=
r1(
−
1
22
)
Squaring, and setting
81
x=
±
7
9
49
+
9
x =
7
=
f. R(s)=1/s,
e
ss
=
)t(c
lim
t
∞→
=
g. closed-loop poles for
s2
4
ζ+
2r
±
,
, gives
130
9
r
1s
lim
⋅
s
0
s
→
, 9s2
+
+
⋅
h. Rise time
T
r
=
1
ω
d
i. 2% settling time
T
s
⎛
⎜
⎜⎜
⎝
=
1
1
−
tan
−π
ζ−
ζ
4
ζω
n
=
sec4
10
9
s
s
2
2
=
10
9s2
+
+
22j1
±−=
⎞
⎟
⎟⎟
⎠
1
22
=
⋅
91.1
rad
=
s68.0
P5.3 Consider each of the following closed-loop transfer functions. By considering the location of the poles on
the complex plane, sketch the unit step response, explaining the results obtained:
Solution:
a.
e
ss
b.
)s(C
)s(R
=
20
s)(2s(
+
+
)t(c
=
=
lim
t
∞→
)s(C
)s(R
=
lim
0
s
→
6
+
)3s)(2s)(1s(
+
+
)10
1s
⋅
s
⋅
20
s)(2s(
+
+
)10
=
1
Dominant time constant=1s,
e
ss
=
lim
s
0
→
1s
⋅
s
⋅
6
+
)3s)(2s)(1s(
+
+
=
1
7
c.
Poles:
=
)s(C
)s(R
s
−=
=ω
n
from Fig4.7,
1s
⋅
s
=
e
ss
⋅
lim
0
s
→
)s(C
)s(R
=
d.
Solving
2
s
9
3s
++
5.0j5.0
±
,3
=ζ
PO
=
9
3s
++
10
+
+
5s4
s
2
2
11
732.1/5.0
,%35
c
Δ
=
3
+
)5.0s)(5s4
=+
0
gives
=ζ
sin
=θ
=
54.0
=
89.0
s(
s 2
+
2
5
=
=
289.0
05.4
s
±−=
2
j
No overshoot, however, real axis pole will dominate with a time constant at 2s.ess=10/2.5=4
e.
)s(C
)s(R
=
Solving
gives
1
2
+
=+
s)(5s(
+
s2
5s2
+
2j1
s
±−=
)5s2
+
0
=ζ
sin
=θ
1
5
=
45.0
The complex pair will dominate the real axis pole,
ess=1/25=0.04
=
f.
)s(C
)s(R
Solving
Solving
100
2
s)(13
2
s(
+
+
s6
+
+
s2
13
s6
+
s2
40
s12
+
pair closer
+
)40
s12
+
j23
s
±−=
s
j26
±−=
to im
aginary axis will
,gives
, gives
Complex
dominate those further away.
ζ
dom
=
3
13
=
83.0
High damping, no or little overshoot, Steady state output=100/(13*40)=0.192
P5.6 Determine the magnitude of the first overshoot of an under damped, second-order system subjected to a
unit impulse input. Plot a graph of this magnitude as a function of damping ratio ζ.
Solution:
)s(C
)s(R
=
2
ω
n
ζω+
2
2
s
ω+
2
n
n
If R(s)=1 then
)s(C
=
ω
ζω+
n
2
n
)
s(
2
ω+
2
d
)t(c
let
2
n
=
ω
ω
d
)t(dc
dt
ζω−
t
n
e
sin
ω
t
d
where
ω=ω
d
n
2
1 ζ−
= ,
0
ζω−
e
n
ζω−
t
n
sin
ω+ω
t
d
ζω−
t
n
e
d
cos
=ω
t
d
0
hence
=pT
π
φ
at this time ,output
M
=
ςω−
n
φ
ω
d
ω
ω
2
n
d
e
sin
ω
d
φ
ω
d
This is a functi
on of
nω and ζ ,so plot M/
nω
−
ω=
n
⎛
⎜
exp
⎜
⎜
⎝
against
ζ
.
1
−
tg
ς
2
1
ζ−
2
ζ
1
−
ς
⎞
⎟
⎟
⎟
⎠
8