Statevariable PDF
Statevariable PDF
WEBSITE takes you to start page after you have read this chapter.
Start page has links to other chapters.
Problem: 7.1 State equations of a mechanical system
Write the state equations of the system, choosing displacement and velocity of the mass as
state variables, x1 and x2, respectively.
Y (t)
K M
F (t)
Fig.1
(b) Write the difference equation obeyed by a sum of money in a bank account at a yearly
interest rate of R percent, compounded every T years. Denote the amount immediately
following the kth interest period as A (kT). Find the differential equation which results as
T 0.
Solution
fM+ fB + fs = F
Let x1 = y, x2 = dy/dt
dx1/dt = x2
equation
TOP
Problem: 7.2 Optimal control system based on quadratic performance index
dx/dt = Ax + Bu
where x = x1, A = 0 1 B= 0
x2 0 0 1
U =-Ktx = -k1x1-k2x2.
Determine the constants k1and k2 so that the following performance index is minimized:
J = xtx dt
0
The initial condition is x(0) = 2
0
Solution
dx/dt =Ax-BKtx = 0 1 x1
-k1 -k2 x2 (1)
A-BKt = 0 1
-k1 -k2
Therefore A –BKt = 0 1
-4 -k2
P= (A-BKt )t P +P (A-BK) = -1
Or,
0 -4 p11 p12 + p11 p12 0 1 = -1 0
1 -k2 p12 p22 p12 p22 -4 - k2 0 -1
Therefore,
= ( 5/k2 ) + ( k2/4 )
TOP
Problem : 7.3 Transfer function of the system from the state equations
(a)
(ii) Is the obtained state variable form in Part (i) a unique description of the differential
equation?
dx/dt = Ax + Bu
y = Cx
where x = x1 , A= 0 1 0
x2 0 0 1
x3 0 -4 -5
B= 1 , C= 1 0 0
-1
5
Solution
(a)
(ii) No
(b) (i)
G(s) = C (sI-A)-1 B
(sI –A) = s -1 0
0 s -1
0 4 s+5
TOP
Problem : 7.4 State equations of the system from the transfer function
dx/dt =Ax + Bu
y =Cx + Du
where
A = -1 1 B= 1 0 1
0 -2 0 1 1
C= 1 2
1 0
1 1
In addition, D is a null matrix of proper order. Find the transfer function matrix,
[Y(s)/U(s)].
(b) The overall transfer function of a single input single output control system is given by
Y(s)/U(s) = ( s2+4s+4)/(s3+5s2+4s)
i Find the A,B and C matrices in the state variable description of the system
Solution
(b)
X1(s)
1/s 1
-1
1/s X3(s)
1/3
-4
dx1/dt = u
dx2/dt = x2+u
dx3/dt = - 4x3 +u
and y = x1-x2/3+x3/3
Problem : 7.5 Damping ratio of the system by minimizing the quadratic performance
index
1
u + x c
s (s + )
-
Fig.2
Determine the damping ratio 0 so that, for unit-step input u (t), the following performance
index is minimized:
J = xT (t) Q x (t) dt
0
where
x1 x 1 0
x = = , Q = ( 0)
x2 dx /dt 0
Solution
d2c/dt2 + dc/dt + c = u
Noting that
d2x/dt2 + dx/dt +x =0 (t 0)
dx2 -1 - x2 1
dt
x1 (0) 1
x2 (0) 0
where A = 0 1
-1 -
At P + PA = -Q
or,
-2p12 = -1
2p12-2 p22=-
0.5 1+
2
TOP
Problem : 7.6 Discrete time state equation from the continuous time state equation
dx1 1 1 x1 1
dt
= + u
dx2 0 1 x2 1
dt
where u (t) is a unit step. The initial conditions are xT (0) = [0 1].
dx1 0 1 x1 0
dt
= + u
dx2 0 -1 x2 1
dt
Obtain the equivalent discrete–time equation, suitable for a time step, T , of 1 sec.
Solution
(a)
x1 et + tet -1
= 2 Answer
x2 et - 0.5
(b) The solution of the equation is
t
X (t) = exp (A (t-t0)) x (t0) + exp (A (t- )) Bu ( ) d
To
Setting t0 = kT and t = (k+1) T and assuming u (t) = U(kT) for kT t (k+1)T, we get
=
x (k+1) = x1 et +tet -1
=2 Answer
x2 et - 0.5
For given A,
1 1- e-t
exp (At) = Setting T= 1 sec. we get
-t
0 e
(T) =
1 0.632
exp (At) = Setting T= 1 sec. we get
0 0.368
(T)= exp (A ) Bd
t
(1-exp (- ) d
0
t
= exp (- ) -d
0
0.368
=
0.632
Therefore,
TOP
Problem : 7.7 Equivalent electrical network and state equations for a mechanical
system
The suspension system for an automobile is designed to increase the comfort of the
passengers by absorbing the vibrations due to the terrain of the road. Fig.3 shows
the model of such a system. Here M represents the masses, B represents the
dashpot damping coefficients, and K represents the spring constants respectively.
M2
M1
u2 (t) = force due to winds
Passenger area
x2
B2 K2
Lower body
x1
B1 K1
Wheels
Road
Solution
K2
u1 B1 K1 M1 M2 u2
(b)
We need four state variables, x1, dx1/dt, x2, and dx2/dt. The state eqns. are:
dx/dt=Ax+Bu
where A= 0 1 0 0
-(K1+K2)/M1 -(B1+B2)/M K2/M1 B2/M1
0 0 0 1
K2/M2 B2/M2 -K2/M2 -B2/M2
B= 0 0 dx/dt = x1 u= u1
1/M1 0 dx1/dt u2
0 0 x2
0 1/M2 dx2/dt
Problem : 7.8 Controllability, observability, and stability from the state equations
0 0 -8 2
dx/dt = 1 0 -14 x + 1 u
0 1 -7 0
c = (0 0 1) x
Solution
U= 2 0 -8 Rank U = 2, Uncontrollable
1 2 -14
0 1 -5
Stable
TOP
Problem : 7.9 State equations of the system from the state diagram
(a) Fig.4 shows a feedback control system. The open –loop transfer function is
G(s) =
2(s+1)(s+3)
s(s+2)(s+4)
What is the closed-loop transfer function of the system?
G(s)
U(s) + C(s)
Fig.4
(b) Fig.5 shows the state diagram of the system of Part (a), The state variables are used to
provide the output signal.
b
U(s) 1 1/s 1/s 1/s C(s)
x3 x2 x1 d g
e
f
Fig.5
Write the state variable equations of the system in the form:
c (t) = D x (t)
Solution
T(s) =(2s2+8s+6)/(s3+8s2+16s+6)
Dx/dt = 0 1 0 0
0 0 1 + 0 u(t)
-6 -16 -8 1
TOP
Problem : 7.10 State diagram of the system from the state equations
where x1 (t), x2 (t) are the system states and r (t) is the reference input. Determine the system
states for t 0 when the input r (t) =1 for t 0; that is,
r (t) =u (t). The initial state of the system is
x1 (0) = 1
x2 (0) = 1
(b) Draw the state diagram for the system described in Part (a) by its state equation and initial
condition. Assign the state variables as the outputs of the integrators and the reference input and
the initial conditions as the system inputs.
Solution
(a)
A= 0 1 B= 0 x (0) = 1
-2 -3 1 1
sI-A = s -1
2 s+3
x (t) =
r 1 dx2/dt dx1/dt = x2 x1
1/s 1/s
TOP
Problem : 7.11 Analogue simulation of the system from the transfer function
(a) A robot-arm drive system for one joint can be represented by the differential equation:
where v (t) =velocity, y (t) =position, and i (t) is the control-motor current. Write the second-
Solution
(a) Let x1 = y, x2= dy/dt=v, u= i. Then
dx1/dt = x2
The state equations lead to the following simple analogue computer simulation.
+u integrator
x3 integrator
x2 integrator
x1= y
7 14 8
+ +
+ +
TOP
Problem : 7.12 Damping ratio of the system by minimizing the quadratic performance
index
Determine the value of the damping ratio 0 so that when the system shown in Fig.6 is
subjected to a unit step input, the following performance index is minimized:
where
x = x1 = e Q= 1 0
x2 de/dt , 0 2
Assume the system to be at rest initially.
Fig.6
Solution
d2c/dt2 +3 dc/dt+c=r
A= 0 1
-1 -3
Note that, since e= r-c, and applied input is a step function, the initial conditions for the state
variables are:
x1 (0) = 1
x2 (0) = 0
J = xT (0) P x (0)
P is determined from
ATP+PA= -Q
(Elements of symmetric P matrix are p11, p12, and p22)
Solving the above eqn. we get
p11 = (1/2 ) +3 /2
p12 = ½
p22 = 1/2
Problem : 7.13 Transfer function of the system from the state equations
(a) An integrator whose transfer function is 1/s has a negative feedback path with an adjustable
gain K.
(i) Demonstrate the manner in which the gain K affects the closed-loop transfer function and the
response of the system, when it is varied from zero to plus and minus infinity.
(ii) Draw the locus of the closed loop pole in the complex s- plane.
(b) Obtain the transfer function Y(s)/U(s)of the system shown in Fig.7 where x1 and x2 are the
state variables.
-5
+
+ x1
2 dt
+
u + y
3
+
-1
+
dt 2
5
+
+
-1 x2
Fig.7
Solution
(a)
The closed loop transfer function is T(s) =1/(s+K). Pole is at s = -K. Addition of the feedback
path has caused the transfer function pole to shift from the origin to the left in the s-plane as K
is increased from some positive value to . While increasing K, the system has been changed from
an integrator to a first-order low –pass filter whose half-power bandwidth is K rad/sec. For
negative K , system will be unstable. Locus of closed loop poles in the complex plane as K is
varied is sown below.
K= 40 K= 30 K= 20 K= 10 K= 0
(b)
dx2/dt = 3x1-x2 + 5u
y = x1+ 2x2
-1
= [1 2] s+5 1 2 = (12s+59)/(s+2) (s +4)
-3 s+1 5
TOP
dx1/dt 0 1 0 x1 0
dx2/dt = 0 0 1 x2 + 0 u
dx3/dt -a3 -a2 -a1 x3 1
dx2/dt = x1-2x2 +u
y1 = 3x1+x2
y2 = x1+2x2
y3 = x2
Find the time solution of the output variables when the input u =1 for t 0 and 0 for t 0.
Solution
(a)
dx1/dt = x2
dx2/dt = x3
dx3/dt =- a3x1- a2x2- a1x3 + u.
These first order equations are then written in the matrix-vector form.
(b)
Y3=(2/3) +a + (b/3)
TOP
(a) Describe a system involving feedback that arises in the following items:
i. Automobiles
ii. Economics
Attempt to identify the plant or process, the input (both reference and disturbance), and the
feedback path(s).
(b) A multi-input /multi-output plant with a state-variable feedback controller is shown in Fig.8.
W U X Y
B 1/s C
Controller
Fig.8
0 1 0 0
A= 0 0 1 B= 0
0 -30 -11 1
C= 1 0 0
Find the feedback matrix K such that the poles are at s= -2,-5, -6.
Solution
(a)
TOP
A hoist may be thought of as a motor lifting a mass in the presence of gravity, as shown in Fig.9
Kb
R L Jdrum
KT
I JM D
V
Load
M
g
Fig.9
Write the appropriate equations and develop a state variable model. Kb is the back emf constant
and KT is the torque constant. JM and Jdrum are the moments of inertia of the motor and drum
respectively.
Solution
dx/dt= 0 1 0 x1 - 0 MgD/2J + 0 V
0 0 KT/J x2 1 0
0 -Kb/L -R/L x3 0 1/L
TOP
Problem : 7.17 Design of the observer for given observer poles
( a) State the conditions for complete observability in terms of transfer functions or transfer
matrices.
(b)For a unity feedback closed loop system having an open loop transfer function
2
G(s) = 100/(s(s +0.1s +10)),
write the state equations and design an observer for observer poles at - 50, -10-j10, -10+ j10.
Solution
dx1/dt = x2
dx2/dt = 1x2+x3
K0 K1 K2
.1
100
s (s2 +.1s +10)
PLANT
s3 +70s2+1200s +10000=0
TOP
A d.c shunt motor, Fig.10, drives an inertia load. The shaft connecting motor to loads is not
rigid and must be considered as spring. The applied voltage is V.
R=armature resistance
L=armature inductance
Kb=back emf constant
KT= Torque constant
Ks=spring constant
JM=rotor moment of inertia
JL=load moment of inertia
I =armature current
M=rotor angle
L=load angle, Write the state equations of the control system using M, L, I, d M/dt, d L/dt as
the state variables.
I
V R L JM Ks JL
M L
FIG. 10
Solution
V= IR +L dI/dt + Kb d M/dt
KTI = JM d M/dt + Ks ( M- L)
0= JL d L/dt - Ks ( M- L)
x1 = L
x2 = d L/dt
x3 = M
x4 = d M/dt
x5 =I
0 1 0 0 0
0 0 - Ks/ JL Ks/ JL 0
A= 0 0 0 1 0
Ks/ JM 0 Ks/ JM 0 KT/ JM
0 0 0 Kb/L R/L
BT= ( 0 0 0 0 Kb/L )
TOP
Problem : 7.19 Transfer function of the system from the state diagram
(a) Write down the state and output equations of the system shown in Fig.11.
(b) Obtain its transfer function, C(s)/R(s)
- 10
+ dt
4 +
+
r (t) 6 + c (t)
-2
+
+ +
10 4
+
dt
-2
Fig.11
Solution
:
dX1/dt=-10x1-2x2+44
dX2/dt=-6x1-2x2+104
y=x1+4 x2
A= -10, -2
6, -2 B= 4
10 C= 1,4
TOP
A = 1 1 , b= 0 ,
2 0 1
Ct = [1 0]
The plant input, u and the system input, w are related by
U = w + Kx .
.
w u + dx/dt x y
b 1/s
Ct
+ +
Fig.12
s = -1,-2.
Solution
A= 1 1 B= 0
2 0 1
C= 1 0
Anew = A + BK = 1 1
2+K1 K2
det (sI – Anew) = 0 . This gives s2-(K2+1) s-2K1=0. Compare with desired equation and obtain
K1=0, K2=-4
TOP
Problem : 7.21 Transfer function of the system from the state diagram
dx/dt = Ax+Bu
y = Cx+Du
1/8
+
+ x1 +
y1 (t)
u1(t) + +
-3/8 1/8 2 4
+ +
x2
+ + y2 (t)
u2(t) +
5/8
A= 1/8 -3/8 B= 1 0
1/8 5/8 0 1
C= 1 4
2 1
D = a null matrix of proper order , find the transfer function matrix relating the input vector
with the output vector.
Solution
= a/e, b/e
c/e, d/e
where a= 8s-1,b=32s-1,c=16s-9,d=8s-7,e=8s2-6s+1
TOP
5 0 2
A = 0 2 -3
2 1 2
(b) Find the characteristic polynomial and the eigenvalues of matrix B.
0 1 0
B= 0 0 1
-4 -6 -4
(c) Find the state-transition matrix of the system described by the following state
equations:
dx1/dt = x2
dx2/dt = -2x1-3x2
Solution
(a)
Inverse =(1/27) 7 2 -4
-6 6 15
-4 -5 10
TOP
(a) State space analysis is not applicable to linear time invariant systems
(b) Conventional root-locus and frequency response methods are applicable to the design
of optimal and adaptive control systems that are mostly time varying and/or nonlinear.
( c) State space analysis is a powerful method for time-varying systems, nonlinear
systems, and multiple- input-multiple output systems
(d)Frequency response methods are a powerful method for time-varying systems,
nonlinear systems, and multiple- input-multiple output systems
Ans: ( c)
Fig.6
where ξ is the damping ratio. The system is initially at rest and is subjected to
a unit –step input. The performance index of the system
t │e(t)│dt
0
is
(a) 2 ξ
(b) 4 ξ 2 -1 +2ξ
(c) 4 ξ 2 -1
(d) ξ
Ans: ( c )
dx1/dt =x1
dx2/dt =-x1 +u
y = x1
Ans: ( b )
(iv ) For the third order plant which of the following may be a set of state variables?
Ans: ( a )
(v) Which one of the following minimum information is necessary to formulate the
problem of control systems optimization?
(a) System state equation, output equation and the control vector only
(b) System state equation, output equation, control vector, and the performance index only
(c) System state equation, output equation, control vector, Constraints of the problem, the
performance index, and system parameters
Ans: ( c )
(vi) Which one of the following characteristics are the attributes of an ideal relay?
(a)
OUTPUT
+
INPUT
-
(b)
OUTPUT
+
INPUT
-
(c)
OUTPUT
+
INPUT
Ans: (a)
(vii) In a series RLC circuit containing a voltage source v (t), the current i and the voltage
vc across the capacitor are chosen as the state variables x1 and x2 respectively. Show that
the state equations are:
(a) dx1 R 1 x1 1
dt = L L + L v (t)
dx2 1 0 x2 0
dt C
(b) dx1 -R -1 x1 1
dt = L L + L v (t)
dx2 1 0 x2 0
dt C
(c) dx1 -R -1 x1 1
dt = L L + C v (t)
dx2 1 0 x2 0
dt L
Ans: ( b )
(viii) In the RLC circuit of the above question, part (v), show that the characteristic
equation of the circuit is
(a) RCs2+LCs +1 =0
(c) LCs2+RCs +1 =0
Ans: ( c )
(ix )
dx1/dt -1 0 x1 2 u
dx2/dt = 0 -2 x2 + 5
dx1/dt -1 0 x1 2 u
dx2/dt = 0 -2 x2 + 0
(c) The following system is not completely observable:
dx1/dt 1 1 x1 0 u
dx2/dt = -2 -1 x2 + 1
y = (1 0) x1
x2
Ans: ( b )
( x ) With the knowledge of the state variable description of a plant, the transfer function of the plant
Ans: (b )
(a) 0 1
-5 -4
(b) 0 1
5 4
(c) -5 -4
0 1
Ans: ( a)
( xiii ) In a control system
Ans: (c )
( xiv ) The state transition matrix (t) possesses the following properties:
Ans: ( a )
(a) The frequency domain techniques are applicable to non-linear and time-varying systems.
(b) The frequency domain techniques are particularly useful for multivariable control systems
(c) Time- domain methods cannot be utilized to design an optimum control system
(d) Time- domain techniques can be utilized for non-linear, time-varying, and multivariable
systems.
Ans:( d )
dx/dt=Ax + Bu
y=Ctx
Which one of the following describes the transfer function, G(s)= y(s)/u(s), for this system?
(a) can be uniquely determined from the state variable equations and the state
variable equations can also be uniquely determined from the transfer function.
(b) for a given set of state variable equations is not unique, but the state variable
equations can be uniquely determined from the transfer function.
(c) can be uniquely determined from the state variable equations, but the state
variable equations for a given transfer function are not unique.
Ans: (a)
dx1/dt =-3x1-x2+u
dx2/dt =2x1
y = x1-x2
is controllable
dx1/dt = -3x1-x2+u
dx2/dt = 2x1
y = x1-x2
is observable.
dx1/dt=-2x1+3x2+u
dx2/dt=x1
y=x1+u
is unobservable.
Ans: (c)
(xix) For the state-transition equation,
t
x (t) = (t-t0) x (t0)+ (t- ) Br ( ) d ,
t0
the Laplace transformed equation is
Ans: (b)
( xx ) For a third-order plant, which of the following may be a set of state variables?
Ans: (a)
(a) dx1/dt= x2
dx2/dt =-x1+u
(b) dx1/dt=-2x1+x2+u
dx2/dt =-x2
(c) dx1/dt=-0.5x1
dx2/dt=-2x2+u
(d) dx1/dt=-x1+2u
dx2/dt=-2x2
Ans: (a)
Ans: (b)
(a) can be uniquely determined from the state variable equations and the state variable
equations can also be uniquely determined from the transfer function.
(b) for a given set of state variable equations is not unique, but the state variable equations
can be uniquely determined from the transfer function.
(c) can be uniquely determined from the state variable equations, but the state variable
equations for a given transfer function are not unique.
(d) cannot be uniquely determined from the state variable equations, and the state
variable equations cannot be uniquely determined from the transfer function.
Ans:( c)
dx1/dt =x1
dx2/dt = -2x1-3x2
is
(a) 2a -b
2b -a
( xxv ) The number of rabbits is x1 , and if left alone would grow indefinitely until the food
supply was exhausted.
dx1/dt = Kx1
where x2 is the number of foxes. Now, if the foxes must have rabbits to exist, we have
(a) h K, ab Kh
(b) h K, ab Kh
(c) h K, ab Kh
Ans:a
(xxvi ) In the problem of rabbits and foxes given in Part (xv), if K=1, h=3, and a=b=2, then
Ans :b
(xxvii ) In a series RLC circuit, if the current in the circuit and the voltage across the capacitor
are selected as the state variables, x1 and x2, respectively, then the state space description of the
circuit is given by
(a)
(b)
(c)
Ans: c
dx1/dt =-3x1-x2+u
dx2/dt =2x1
y=x1-x2
is controllable.
dx1/dt = -2x1+3x2+u
dx2/dt = x1+ u
y = x1
is uncontrollable.
is observable.
dx1/dt =-2x1+3x2+u
dx2/dt =x1
y = x1+u
is unobservable.
Ans: (d)
(xxix) Which of the systems with the following state space descriptions is completely
controllable?
(a)
A= -1 0 B= 2
0 -2 5
(b)
A= -1 0 B= 2
0 -2 0
( c)
A= 0 1 B= 1
2.5 -1.5 1
Ans:a
(a)
A= -1 0 Ct = 1
0 -2 3
(b)
A= -1 0 Ct = 0
0 -2 1
(c)
A= 1 1 Ct = 1
-2 -1 0
B= 0
1
Ans: b
TOP