Digitalcontrolsystems Problem
Digitalcontrolsystems Problem
1. Obtain the unit impulse response of the open-loop sampled data system shown in
Zero-order 1/[s(s+1)]
hold
Fig.1
Solution:
From Fig.
G(s)=(1-exp(-st))/[s2(s+1)] =(1-exp(-st))[1/s2 –1/s +1/(s+1)]
The response is
C(0)=0
C(T)= .3678
C(2T)= .7675
C(3T)= .9145 etc.
2. Characteristics of a digital controller (in the difference equations form) to give a specified
response of the system to a unit step function.Fig. 2 shows a digital control system.
Fig. 2
(a) Show that the z transform, D(z), of the digital controller is given by
(s+2)/[s(s+1)].
Determine the characteristics of a digital controller (in the difference equations form) such that
the response of the system to a unit step function will be c (t) = 5(1-e-2t). The sampling period is
1 s.
Use the following table of Laplace and z transforms.
u (t)-1/s-z/(z-1)
t-1/s2-zT/(z-1)2
e-at-1/(s+a)-z/(z-e-aT)
Solution:
For the plant and zero order hold G(s) = G1(s)G2(s) = (1- e-Ts) (s+2)/[s2(s+1)]
G1(s) = 1-e-Ts.
=(1.368z-0.104)/(z2-1.368z+0.368)
For
C(s) =5[(1/s)-(1/(s+2))],
D(z) =[z2-1.368z+.368]/[1.368z-.104].[4.32/(z-4.45)]
= [4.32z2-5.91z+1.59]/[1.368z2-6.202z+.46]
m(k) =4.53m(k-1)-.34m(k-2)+3.16e(k)-4.32e(k-1)+1.16e(k-2)
(a) Derive a state variable equation for a system governed by the second –order difference
equation
(b) Fig.4 shows a position servo operated by an input via a sampler and a zero order hold with
the position c and velocity dc/dt represented by x1and x2 respectively. If the states at the
sampling instants alone are required ,determine the A and B matrices in the following
equation
(c)
u(t)
u*(t)Zero-order hold
c(t) = x1
Process 2/(s2+3s+2) x1(kT)
T dc(t)/dt = x2x2(kT)
T
Fig.4
Solution:
(a)
X2(k+1) =-2x1(k)-3x2(k0+u(k)
(b)
(T) =L-1(sI-A)-1= s -1 -1
2 s+3
Setting T =1,
(T) = 0.1998
0.2326
4. To obtain the state equation of a discrete data system described by its difference equation
(a) Write a short note on the existence and uniqueness of the solution of state equations.
where u (k) is the reference input and c (k) is the output. Show that the state equation of the
system is
x (k+1) = 0 1 x (k) + 1 u (k)
-3 -5 -3
Assume zero initial conditions and describe the technique used to arrive at the above equation.
Solution:
(b)
Taking z-transform of the difference equation, and assuming zero initial conditions yields
=0
Substitution of the state variables into the original difference eqn. gives
(a) Fig.5 shows the block diagram of a control system with feed forward and
feedback paths. Apply Mason’s rule to determine the transfer function, C(s)/R(s).
G2
+
R(s) + + C(s)
G1
-
Fig.5
0.97
y (k) = 0 0 0
6. In the above, xi (k) is the number of students in the ith year of the four-year degree
program in academic year k, u (k) represents the input (that is, the enrollment in the first-
year class in year k) and y (k) is the output (that is, the number of students graduating in
year k).In July 1992, the college had 188 students in the first year class, 126, in the
second year class, 103 in the third year class, and 72 in the fourth year class.
The table below gives the history of admissions to the first year class in the subsequent years.
Year-1993-1994-1995-1996-1997-1998-1999
Admissions-275-320-350-400-450-450-430
Use this model to determine the number of students graduating in May of each year from 1994 to
2002.
Solution:
(a) There are two forward paths from R(s) to C(s), with transmittance G1 and G2. There is on
loop with transmittance –G1H, and this is touching both the forward paths. Hence,
C(s)/R(s) = (G1+G2)/(1+G1H)
x (1992)=
188
126
103
72
(a) Derive a state variable equation for a system governed by the second-order difference
equation:
8. Draw the block diagram of the system, showing the A, B, C, D matrices on it.
C= [1 0 ] D=0
Solution:
A
Double arrow
(b)
As u (k) = 0, we have
x (2) = -4
13
x (3)= 13
40
x (4) = -40
121 and y = (0, 1, -4, 13, -40 etc)
Q = K (P1-P2) ½
Q P1 P2
Fig.6
Fig.6 shows the variables, and K is a constant. Determine a linear approximation to the fluid
flow equation.
(b) A slack loop is maintained between two sets of drive rolls, Fig.7, in a certain computer tape
deck.
1(t) 2(t)
O O
O O
L (t)
Fig.7
All rolls have diameter D. Under normal operating conditions drive speeds are 1(t)= 2(t)= ’
and slack loop length is L (t) = L’. Starting at time t0, there are small independent variations in
the drive speeds so that 1(t)= ’ + 1(t) and 2(t)= ’ + 2(t). Derive an expression for the
resulting change in the loop length l (t).
(c) Obtain a discrete model of the tape-drive system described in Part (b) when a digital
controller is used,
Q=K/(2 P01/2) P
L = L’ +l
kT+T
l (kT+T) = l (KT) +(D/2) [ 1(t)- 2(t)] dt
kT
- (t-kT)/
1(t)- 2(t) =[ 1(kT) des- 2(kT) des](1-e ) Therefore,
(a) How would you modify the Routh-Hurwitz criterion to study the stability of a closed-loop
discrete-time system?
(b) The characteristic polynomial of a closed-loop discrete-time system is given
by z3 + 5z2 + 3z +2 =0.
Solution:
(b)
f(z)= z3+5z2+3z+2=0
Construct Routh table. From Routh table, there is one sign change in the first column.
11. Stability from the characteristic polynomial corresponding to the pulse transfer
function of a sampled data system
transformation z = (r+1)/(r+2)
Solution:
(b)
((r+1)/(r-1))3+5,94((r+1)/(r-1))2+7.7(r+1)/(r-1)-.368=0
or
f( r)=14.27r3+2.3 r2-11.47r+3.13=0
Routh’s array
14.27 -11.47
r3
2.3 3.13
r2
-31.1
r1 3.13
r0
Two sign changes. Two zeroes of he characteristic polynomial f(z) will remain outside the unit
circle. System is UNSTABLE.
Fig.9
Fig.10
Solution:
(b)
A= 0 1
0 -1
1 t=T
= Laplace-1 s -1
0 s+1 t=T
= 1/(s+1) -1/(s(s+1))
0 1/s t=T
= 1 1-e-T = 1 .632
0 e-T 0 .368
T
(T)= ( )Bd
0
T
1-e-
(T)= 1 e- 0 d
0 0 1
T + e-T-1
= -e-T+1
T=1
= .368
.632
Solution:
i. Continuous
ii. Discrete
iii. Continuous
(b) A discrete –time system can be simulated on a digital computer in the same way as a
continuous-time system on an analogue computer. The only difference is that integrators are
replaced by delays. This implies that the blocks containing s-1 are replaced by blocks
containing z-1. Use this approach to draw the simulation diagram for a system represented
by the transfer function
Solution:
4
+
y
u z-1 z-1 z-1
+
1.2
-.45
.05
Problem 9.14: To obtain the equivalent discrete time equation from the state space
description of a continuous time plant
(a) Describe a numerical method for the solution of continuous time state equations.
(b)The state space description for a continuous-time plant is
A= 0 1 B= 0
0 -1 1
Obtain the equivalent discrete-time equation. Take the sampling instant, T = 1s.
Solution:
Setting t0= kT and t = (k+1)T and assuming u(t) = u(kT) for kT less than/equal to t less
than(k+1)T,
we get
(T) = exp(A ) Bd
F(z) = 10z/(z-1)(z-2)
(c) If G*(s) = f(nT) e-nTs, show that the closed form solution
is
G*(s) = 1/(1-e T(s+p))
Solution:
(a) g(t)=e-pt.
Then for t = T
G* = -npT e-nTs
n=0 e
= -nT (p+s)
n=0 e
n -T (p
= n=0 X where X= e
+ s) n=0
X n =1/(1-X)
(c ) X(z)/z = 10 (z-1)(z-2)
or, x(0)=0
x(T)= 10
x(2T)= 30
x(3T)= 70
x(4T)= 150
TOP
(a) The discrete-time system of Fig.12 has the open-loop transfer function:
G(s) = 10
s (s+1)
R + T=1 G(s) C
Fig.12
(b) Obtain the pulse transfer function of the closed-loop system shown in Fig.13.
R +T G(s) C
-
H(s)
T T
Fig. 13
Solution:
(a)
exp(-1)=.368
therefore, z2 +4.952z +.368 =0
z = 0.076, -4.876
Magnitude of one root 1
Therefore, unstable.
(b)
E
R +T G(s) C
B-
H(s)
T T