Notes
Notes
Unit I
State feedback controller design
Pole placement approach
In the conventional approach to the design of a single input, single-output control system,
we design a controller (compensator) such that the dominant closed-loop poles have a desired
damping ratio and a desired undamped natural frequency n. In this approach, the order of the
system may be raised by 1 or 2 unless pole–zero cancellation takes place. Note that in this approach
we assume the effects on the responses of nondominant closed-loop poles to be negligible. Different
from specifying only dominant closed-loop poles (the conventional design approach), the present
pole-placement approach specifies all closed-loop poles.
The pole-placement method is somewhat similar to the root-locus method in that we place
closed-loop poles at desired locations. The basic difference is that in the root-locus design we place
only the dominant closed-loop poles at the desired locations, while in the pole-placement design
we place all closed-loop poles at desired locations.
There is a requirement on the part of the system for the closed-loop poles to be placed at
arbitrarily chosen locations. The requirement is
1. The system must be completely state controllable.
2. All the state variables are measurable and available for feedback.
𝑥̇ = 𝐴𝑥 + 𝐵𝑢
𝑦 = 𝐶𝑥 + 𝐷𝑢 (1)
This means that the control signal u is determined by an instantaneous state. Such a scheme
is called state feedback. The 1×n matrix K is called the state feedback gain matrix. We assume that
all state variables are available for feedback. A block diagram for this system is shown in figure.
𝑥̇ = (𝐴 − 𝐵𝐾) 𝑥(𝑡)
The solution of this equation is given by
𝑥(𝑡) = 𝑒 (𝐴−𝐵𝐾)𝑡 𝑥(0)
where, 𝑥(0) is the initial state caused by external disturbances. The stability and transient response
characteristics are determined by the eigenvalues of matrix(𝐴 − 𝐵𝐾). If matrix K is chosen
properly, the matrix (𝐴 − 𝐵𝐾) can be made an asymptotically stable matrix, and for all 𝑥(0) ≠ 0x,
it is possible to make 𝑥(𝑡) approach 0 as t approaches infinity. The eigenvalues of matrix (𝐴 − 𝐵𝐾)
are called the regulator poles. If these regulator poles are placed in the left-half s plane, then 𝑥(𝑡)
approaches 0 as t approaches infinity. The problem of placing the regulator poles (closed-loop
poles) at the desired location is called a pole-placement problem.
𝑥̇ = 𝐴𝑥 + 𝐵𝑢
The system is completely state controllable and the desired closed loop poles are at 𝑠 = 𝜇1 , 𝑠 =
𝜇2 , 𝑠 = 𝜇3 , … … , 𝑠 = 𝜇𝑛 . Use of the state feedback control
𝑢 = −𝐾𝑥
modifies the system equation to
𝑥̇ = (𝐴 − 𝐵𝐾)𝑥 (1)
Let us define 𝐴̃ = 𝐴 − 𝐵𝐾. The desired characteristic equation is
|𝑠𝐼 − 𝐴 + 𝐵𝐾| = |𝑠𝐼 − 𝐴̃| = (𝑠 − 𝜇1 ) (𝑠 − 𝜇2 ) (𝑠 − 𝜇3 ) . .. (𝑠 − 𝜇𝑛 )
= 𝑠 𝑛 + 𝛼1 𝑠 𝑛−1 + 𝛼2 𝑠 𝑛−2 + ⋯ . . +𝛼𝑛−1 𝑠 + 𝛼𝑛 = 0
Since Cayley-Hamilton theorem states that 𝐴̃ satisfies its own characteristic equation.
We shall utilize Equation (2) to derive Ackermann’s formula. To simplify the derivation, we
consider the case where n=3. (For any other positive integer n, the following derivation can be
easily extended.)
Multiplying the preceding equations in order by 𝛼3 , 𝛼2 , 𝛼1 , and 𝛼0 respectively, and adding the
results, we obtain,
𝛼3 𝐼 + 𝛼2 𝐴̃ + 𝛼1 𝐴̃2 + 𝐴̃3
= 𝛼3 𝐼 + 𝛼2 (𝐴 − 𝐵𝐾) + 𝛼1 (𝐴2 − 𝐴𝐵𝐾 − 𝐵𝐾𝐴̃) + 𝐴3 − 𝐴2 𝐵𝐾 − 𝐴𝐵𝐾𝐴̃ − 𝐵𝐾𝐴̃2
[Link]’s College of Engineering 2 VI Semester ICE
IC6601 Advanced Control System Department of ICE 2018-19
= 𝛼3 𝐼 + 𝛼2 𝐴 + 𝛼1 𝐴2 + 𝐴3 − 𝛼2 𝐵𝐾 − 𝛼1 𝐴𝐵𝐾 − 𝛼1 𝐵𝐾𝐴̃ − 𝐴2 𝐵𝐾 − 𝐴𝐵𝐾𝐴̃ − 𝐵𝐾𝐴̃2 (3)
𝛼2 𝐾 + 𝛼1 𝐾𝐴̃ + 𝐾𝐴̃2
[𝐵 𝐴𝐵 𝐴 𝐵 ] [ 2
𝛼1 𝐾 + 𝐾𝐴̃ ] (4)
𝐾
Since the system is completely state controllable, the inverse of the controllability matrix
[𝐵 𝐴𝐵 𝐴2 𝐵 ]
exists. Premultiplying both sides of equation (4) by the inverse of the controllability matrix, we
obtain
𝛼2 𝐾 + 𝛼1 𝐾𝐴̃ + 𝐾𝐴̃2
[𝐵 𝐴𝐵 𝐴2 𝐵 ∅(𝐴) = [ ]−1 𝛼1 𝐾 + 𝐾𝐴̃ ]
𝐾
Premultiplying both sides of this last equation by [0 0 1] , we obtain
𝛼2 𝐾 + 𝛼1 𝐾𝐴̃ + 𝐾𝐴̃2
[0 0 1][𝐵 𝐴𝐵 𝐴2 𝐵 ]−1 ∅(𝐴) = [0 0 1] [ 𝛼1 𝐾 + 𝐾𝐴̃ ]=𝐾
𝐾
Which can be written as
𝐾 = [0 0. . . 0 1 ] [𝐵 𝐴𝐵 … … … … . . 𝐴𝑛−1 𝐵 ]−1 ∅(𝐴) (5)
This equation gives the required state feedback gain matrix K. this is known as Ackermann’s
formula for the determination of the state feedback gain matrix K
Step 1: Check the controllability of the given system. If it is controllable, then go to step2
Step 3: Determine the transformation matrix ′𝑇′ which transforms system state equation into
controllable canonical form.
[Link]’s College of Engineering 3 VI Semester ICE
IC6601 Advanced Control System Department of ICE 2018-19
𝑎𝑛−1 𝑎𝑛−2 ,…… 𝑎1 1
𝑎𝑛−2 𝑎𝑛−3 ,…… 1 0
𝑇 = 𝑀 × 𝑊 = [𝐵 𝐴𝐵 𝐴2 𝐵 𝑛−1 ]
,……….. 𝐴 𝐵 ; ; ,……. ; ;
𝑎1 1 ,……. 0 0
[ 1 0 ,……. 0 0]
Step 4: Using desired closed loop poles, write the desired polynomial. i.e.
(𝑆 − 𝜇1 )(𝑆 − 𝜇2 )(𝑆 − 𝜇3 ), … … … … … … (𝑆 − 𝜇𝑛 )
Step 1: Check the controllability of the given system. If it is controllable, then go to step2
Note: If the system is of low order (𝑛 ≤ 3), direct substitution of matrix K into the desired
characteristics polynomial may be simpler.
Step 3: Substitute this 𝐾 matrix into the desired characteristics polynomial |𝑆𝐼 − 𝐴 + 𝐵𝐾| and
equate it to (𝑠 − 𝜇1 )(𝑠 − 𝜇2 )(𝑠 − 𝜇3 ).
Step 4: By equating the coefficients of like power of ′𝑆′ on both sides, it is possible to determine
the values of 𝐾1 , 𝐾2 , 𝐾3
Step 1: Check the controllability of the given system. If it is controllable, then go to step2
Step 2: Determine the desired characteristic equation from the specified closed loop poles,
𝜇1 , 𝜇2 , 𝜇3 , . . , 𝜇𝑛
i.e., (𝑠 − 𝜇1 ) (𝑠 − 𝜇2 ) (𝑠 − 𝜇3 ). . . . . . . . (𝑠 − 𝜇𝑛 ) = 𝑠 𝑛 + 𝛼1 𝑠 𝑛−1 +, … … . . , +𝛼𝑛−1 𝑠1 + 𝛼𝑛
and determine the values of 𝛼1 , 𝛼2 , … … . , 𝛼𝑛 .
Step 3: Determine the matrix ∅(𝐴) using the coefficient of desired characteristics polynomial.
𝑥̇ = 𝐴𝑥 + 𝐵𝑢
𝑦 = 𝐶𝑥 (1)
Where, 𝐾 = [𝑘1 𝑘2 , , , … . . 𝑘𝑛 ]
Assume that the reference input is applied at t= 0. Then, for 𝑡 > 0, the system dynamics can
be described by equation (1) and (2), or
𝑥̇ = 𝐴𝑥 + 𝐵𝑢 = (𝐴 − 𝐵𝐾)𝑥 + 𝐵𝑘1 𝑟 (3)
Now we can design a type 1 servo system such that the closed loop poles are located at desired
positions. The designed system will be an asymptotically stable system, 𝑦(∞) will approach
the constant value 𝑟, and 𝑢(∞) will approach 0. At steadystate,
𝑒̇ = (𝐴 − 𝐵𝐾)𝑒 (6)
Equation (6) describes error dynamics. The design of the type 1 servo system here is converted
to the design of an asymptotically stable regulator system such that 𝑒(𝑡) approaches zero, given
any initial conditions 𝑒(0). If the system defined by equation (1) is completely state
controllable, then, by specifying the desired eigenvalues 𝜇1 , 𝜇2 , , … . . 𝜇𝑛 for the matrix 𝐴 − 𝐵𝐾,
matrix 𝐾 can be determined by the pole placement technique. The steadystate values of 𝑥(𝑡)
and 𝑢(𝑡) can be found as follows. At steady state 𝑡 = ∞, from equation (2)
Since the desired eigenvalues of 𝐴 − 𝐵𝐾 are all in the left half s-plane, the inverse of matrix
𝐴 − 𝐵𝐾 exists. Consequently, 𝑥(∞) can be determined as
𝑢(∞) = −𝐾𝑥(∞) + 𝑘1 𝑟 = 0
If the plant has no integrator (type 0 plant), the basic principle of the design of a type 1 servo
system is to insert an integrator in the feedforward path between the error comparator and the plant,
as shown in Figure 1. (The block diagram of Figure 1 is a basic form of the type 1 servo system
where the plant has no integrator.)
𝑥̇ = 𝐴𝑥 + 𝐵𝑢 (1)
𝑦 = 𝐶𝑥 (2)
𝑢 = −𝐾𝑥 + 𝑘1 (3)
= 𝑟 − 𝑦 = 𝑟 − 𝐶𝑥 (4)
Define
𝑥(𝑡) − 𝑥(∞) = 𝑥𝑒 (𝑡)
(𝑡) − (∞) = 𝑒 (𝑡)
𝑢(𝑡) − 𝑢(∞) = 𝑢𝑒 (𝑡)
Then, equation (7) can be written as
𝑥𝑒̇ (𝑡) 𝐴 0 𝑥𝑒 (𝑡) 𝐵
[̇ ]=[ ][ ] + [ ] 𝑢𝑒 (𝑡) (8)
𝑒 (t) −𝐶 0 𝑒 (𝑡) 0
Where, 𝑢𝑒 (𝑡) = −𝐾𝑥𝑒 (𝑡) + 𝑘𝑟 𝑒 (𝑡) (9)
Define a new (𝑛 + 1) th order error vector 𝑒(𝑡) by
𝑥𝑒 (𝑡)
𝑒(𝑡) = [ (𝑡)] = (𝑛 + 1) vector
𝑒
There are times when this will not be necessary, when we will need observation of only the
unmeasurable state variables, but not of those that are directly measurable as well. For example,
since the output variables are observable and they are linearly related to the state variables, we need
not observe all state variables, but observe only 𝑛 − 𝑚 state variables, where 𝑛 is the dimension of
the state vector and 𝑚 is the dimension of the output vector.
An observer that estimates fewer than 𝑛 state variables, where 𝑛 is the dimension of the
state vector, is called a reduced-order state observer or simply, a reduced-order observer. If the
order of the reduced-order state observer is the minimum possible, the observer is called a
minimum-order state observer or minimum-order observer.
State Observer Design: A state observer estimates the state variables based on the measurements
of the output and control variables. Here the concept of observability plays an important role. State
observers can be designed if and only if the observability condition is satisfied.
In many practical cases, the observed state vector is used in the state feedback to generate
the desired control vector. Consider the plant defined by
𝑥̇ = 𝐴𝑥 + 𝐵𝑢 (1)
𝑦 = 𝐶𝑥 (2)
The observer is a subsystem to reconstruct the state vector of the plant. The mathematical
model of the observer is basically the same as that of the plant, except that we include an additional
term that includes the estimation error to compensate for inaccuracies in matrices A and B and the
lack of the initial error. The estimation error or observation error is the difference between the
measured output and the estimated output. The initial error is the difference between the initial state
and the initial estimated state. Thus, we define the mathematical model of the observer to be
𝑥̃ = 𝐴𝑥̃ + 𝐵𝑢 + 𝐾𝑒 (𝑦 − 𝐶𝑥̃)
= (𝐴 − 𝐾𝑒 𝐶)𝑥̃ + 𝐵𝑢 + 𝐾𝑒 𝑦 (3)
where 𝑥̃ is the estimated state and 𝐶 𝑥̃ is the estimated output. The inputs to the observer
are the output y and the control input u. Matrix 𝐾𝑒 , which is called the observer gain matrix, is a
weighting matrix to the correction term involving the difference between the measured output y
and the estimated output 𝐶 𝑥̃. This term continuously corrects the model output and improves the
performance of the observer. Figure shows the block diagram of the system and the full-order state
observer.
Full-Order State Observer. The order of the state observer that will be discussed here is the same
as that of the plant. Assume that the plant is defined by Equations
𝑥̇ = 𝐴𝑥 + 𝐵𝑢 (1)
𝑦 = 𝐶𝑥 (2)
𝑥̃ = 𝐴𝑥̃ + 𝐵𝑢 + 𝐾𝑒 (𝑦 − 𝐶𝑥̃)
= (𝐴 − 𝐾𝑒 𝐶)𝑥̃ + 𝐵𝑢 + 𝐾𝑒 𝑦 (3)
To obtain the observer error equation, let us subtract Equation (3) from Equation (1)
(4)
Define the difference between x and 𝑥̃ as the error vector e, or
𝑒 = 𝑥 − 𝑥̃
Then Equation (4) becomes
𝑒̇ = (𝐴 − 𝐾𝑒 𝐶)𝑒 (5)
From Equation (5), the dynamic behavior of the error vector is determined by the
eigenvalues of matrix 𝐴 − 𝐾𝑒 𝐶. If matrix 𝐴 − 𝐾𝑒 𝐶 is a stable matrix, the error vector will converge
to zero for any initial error vector e(0). That is, will converge to x(t) regardless of the values of x(0)
and If the eigenvalues of matrix 𝐴 − 𝐾𝑒 𝐶 are chosen in such a way that the dynamic behavior of
the error vector is asymptotically stable and is adequately fast, then any error vector will tend to
zero (the origin) with an adequate speed.
If the plant is completely observable, then it can be proved that it is possible to choose
matrix 𝐾𝑒 such that 𝐴 − 𝐾𝑒 𝐶 has arbitrarily desired eigenvalues. That is, the observer gain matrix
𝐾𝑒 can be determined to yield the desired matrix 𝐴 − 𝐾𝑒 𝐶.
Duality Property
The design of the full-order observer becomes that of determining an appropriate 𝐾𝑒 such
that 𝐴 − 𝐾𝑒 𝐶 has desired eigenvalues. Thus, the problem here becomes the same as the pole-
placement problem. In fact, the two problems are mathematically the same. This property is called
duality.
In designing the full-order state observer, we may solve the dual problem, that is, solve the
pole-placement problem for the dual system.
𝑧̇ = 𝐴∗ 𝑧 + 𝐶 ∗ 𝑣
𝑛 = 𝐵∗𝑧
If the dual system is completely state controllable, then the state feedback gain matrix 𝐾
can be determined such that matrix 𝐴∗ − 𝐶 ∗ 𝐾 will yield a set of the desired eigenvalues. If
𝜇1 , 𝜇2 , 𝜇3 … . . 𝜇𝑛 are the desired eigenvalues of the state observer matrix, then by taking the same
𝜇𝑖′ 𝑠 as the desired eigenvalues of the state-feedback gain matrix of the dual system, we obtain
Noting that the eigenvalues of 𝐴∗ − 𝐶 ∗ 𝐾 and those of 𝐴 − 𝐾 ∗ 𝐶 are the same, we have,
Comparing the characteristic polynomial |𝑠𝐼 − (𝐴 − 𝐾 ∗ 𝐶)| and the characteristic polynomial for
the observer system |𝑠𝐼 − (𝐴 − 𝐾𝑒 𝐶)|, we find that 𝐾𝑒 and 𝐾 ∗ are related by
Thus, using the matrix 𝐾 determined by the pole-placement approach in the dual system, the
observer gain matrix 𝐾𝑒 for the original system can be determined by using the relationship
𝐾𝑒 = 𝐾 ∗
Procedure for calculating Observer gain
Method I : Using Transformation Matrix T
Step 1: Check the observability of the given system. If it is observable, then go to step2
𝐶
𝐶𝐴
| 𝐶𝐴2 |
Observability matrix |𝑶𝑴| = ≠𝟎
| ; |
𝐶𝐴𝑛−2
𝐶𝐴𝑛−1
Step 2: Find the characteristic polynomial
𝑄 = (𝑊𝑁 ∗ )−1
and
Step 4: Using desired closed loop poles, write the desired polynomial. i.e.
(𝒔 − 𝝁𝟏 )(𝒔 − 𝝁𝟐 )(𝒔 − 𝝁𝟑 ), … … … … … … (𝒔 − 𝝁𝒏 )
Step 1: Check the observability of the given system. If it is observable, then go to step2
𝐶
𝐶𝐴
| 𝐶𝐴2 |
Observability matrix |𝑶𝑴| = ≠𝟎
| ; |
𝐶𝐴𝑛−2
𝐶𝐴𝑛−1
𝒌𝒆𝟏
Step 2: 𝒏 = 𝟐, then 𝑲𝒆 = [ ]
𝒌𝒆𝟐
Step 3: Substitute this 𝑲𝒆 matrix into the desired characteristics polynomial |𝒔𝑰 − (𝑨 − 𝑲𝒆 𝑪)|
and equate it to (𝒔 − 𝝁𝟏 )(𝒔 − 𝝁𝟐 )(𝒔 − 𝝁𝟑 ).
Step 4: By equating the coefficients of like power of ′𝑺′ on both sides, it is possible to determine
the values of 𝑲𝒆𝟏 , 𝑲𝒆𝟐
Step 1: Check the observability of the given system. If it is observable, then go to step2
Step 3: Determine the matrix ∅(𝑨) using the coefficient of desired characteristics polynomial.
∅(𝑨) = 𝑨𝒏 + 𝜶𝟏 𝑨𝒏−𝟏 + 𝜶𝟐 𝑨𝒏−𝟐 +, … … . . , +𝜶𝒏−𝟏 𝑨𝟏 + 𝜶𝒏 𝑰
Step 4: Calculate the state feedback gain matrix, ′𝑲𝒆 ′ using the Ackermann’s formula
𝐶 −𝟏 𝟎
𝐶𝐴 𝟎
𝑲𝒆 = [∅(𝑨)] 𝐶𝐴2 ;
; ;
[𝐶𝐴 ]𝑛 [ 𝟏]
𝑥̇ = 𝐴𝑥 + 𝐵𝑢
𝑦 = 𝐶𝑥 (1)
(2)
Where, 𝐴𝑎𝑎 = scalar
𝐴𝑎𝑏 = 1 × (𝑛 − 1) matrix
𝐴𝑏𝑎 = (𝑛 − 1) × 1 matrix
𝐴𝑏𝑏 = (𝑛 − 1) × (𝑛 − 1) matrix
𝐵𝑎 = scalar
𝐵𝑏 = (𝑛 − 1) × 1 matrix
From (2) the equation for the measure portion of the state becomes
[Link]’s College of Engineering 13 VI Semester ICE
IC6601 Advanced Control System Department of ICE 2018-19
𝑥̇ 𝑎 = 𝐴𝑎𝑎 𝑥𝑎 + 𝐴𝑎𝑏 𝑋𝑏 + 𝐵𝑎 𝑢
The term on the left side of (3) can be measured. It act as output equation. In designing minimum
order observer, the left side of (3) to be known quantities. Thus it relates the measurable quantites
and unmeasurable quantities of the state. From (1) the equation for the unmeasured portion of the
state becomes
𝐴𝑏𝑎 𝑥𝑎 and 𝐵𝑏 𝑢 are known quantities and it describes the dynamics of unmeasured portion of the
state.
Design of Minimum order observer
The state equation for the full order observer is
𝑥̇ = 𝐴𝑥 + 𝐵𝑢
The state equation for the minimum order observer is
𝑥̇ 𝑏 = 𝐴𝑏𝑎 𝑥𝑎 + 𝐴𝑏𝑏 𝑥𝑏 + 𝐵𝑏 𝑢
The output equation for the full-order observer is
𝑦 = 𝐶𝑥
The output equation for the minimum-order observer is
𝑥̇ 𝑎 − 𝐴𝑎𝑎 𝑥𝑎 − 𝐵𝑎 𝑢 = 𝐴𝑎𝑏 𝑋𝑏
Where the state observer gain matrix 𝐾𝑒 is an (𝑛 − 1) × 1 matrix. To estimate 𝑥̂𝑏 we need the
derivative of 𝑥𝑎 . This presents a difficulty, because differentiation amplifies noise, if 𝑥𝑎 (= 𝑦)is
noisy, the use of 𝑥̇ 𝑎 is unacceptable.
Table 1. List of necessary substitutions for writing the observer equation for the minimum-order
observer
Define
Solved problem
[Link]’s College of Engineering 15 VI Semester ICE
IC6601 Advanced Control System Department of ICE 2018-19
1. Consider the regulator system shown in following figure. The plant is given by
𝒙̇ 𝟏 𝟎 𝟏 𝟎 𝒙𝟏 𝟎
[𝒙̇ 𝟐 ] = [ 𝟎 𝟎 𝟏 ] [𝒙𝟐 ] + [𝟎] 𝒖(𝒕);
𝒙̇ 𝟑 −𝟏 −𝟓 −𝟔 𝒙𝟑 𝟏
The system uses the state feedback control 𝒖 = −𝑲𝒙. The desired eigenvalues are 𝝁𝟏 = −𝟐 +
𝒋𝟒, 𝝁𝟏 = −𝟐 − 𝒋𝟒, 𝝁𝟑 = −𝟏𝟎. Determine the state feedback gain matrix K
Solution
0 1 0 0
Given 𝐴 = [ 0 0 1 ]; 𝐵 = [0]; 𝜇1 = −2 + 𝑗4, 𝜇2 = −2 − 𝑗4, 𝜇3 = −10; 𝑂𝑟𝑑𝑒𝑟, 𝑛 = 3
−1 −5 −6 1
Method I : Using Transformation Matrix T
Step 1: Check the controllability of the given system. If it is controllable, then go to step2
Controllability matrix 𝐶𝑀 = |𝐵 𝐴𝐵 𝐴2 𝐵 | ≠ 0
0
𝐵 = [0] ;
1
0 1 0 0 0
𝐴×𝐵 = [ 0 0 1 ] × [0] = [ 1 ]
−1 −5 −6 1 −6
0 1 0 0 1
𝐴2 𝐵 = 𝐴 × 𝐴𝐵 = [ 0 0 1 ] × [ 1 ] = [−6]
−1 −5 −6 −6 31
0 0 1
1 −6 0 −6 0 1
|𝐶𝑀| = |𝐵 𝐴𝐵 𝐴 𝐵 | = [0 1 −6] = 0 |
2 | − 0| | + 1| | = −1 ≠ 0
−6 31 1 31 1 −6
1 −6 31
Hence the system is completely state controllable.
𝑠 0 0 0 1 0 𝑠 −1 0
|𝑆𝐼 − 𝐴| = |[0 𝑠 0] − [ 0 0 1 ]| = |0 𝑠 −1 |
0 0 𝑠 −1 −5 −6 1 5 𝑠+6
𝑠 −1 0 −1 0 𝑠
= 𝑠| | − (−1) | | + 0| | = 𝑠(𝑠 2 + 6𝑠 + 5) + 1
5 𝑠+6 1 𝑠+6 1 5
= 𝑠 3 + 6𝑠 2 + 5𝑠 + 1 = 0
Step 4: Using desired closed loop poles, write the desired polynomial. i.e.
(𝒔 − 𝝁𝟏 )(𝒔 − 𝝁𝟐 )(𝒔 − 𝝁𝟑 ), … … … … … … (𝒔 − 𝝁𝒏 )
𝐾 = [𝛼3 − 𝑎3 𝛼2 − 𝑎 2 𝛼1 − 𝑎1 ]𝑇 −1
1 0 0
Since the system is in controllable canonical form 𝑇 = 𝐼 = [0 1 0]
0 0 1
1 0 0
𝐾 = [200 − 1 ]
60 − 5 14 − 6 [0 1 0]
0 0 1
Hence, 𝐾 = [199 55 8]
𝑥1
𝑢 = −𝐾𝑥 = −[𝐾1 𝐾2 𝐾3 ] [𝑥2 ]
𝑥3
𝑥1
= −[199 55 8] [𝑥2 ]
𝑥3
𝑢 = −199𝑥1 − 55𝑥2 − 8𝑥3
Step 1: Check the controllability of the given system. If it is controllable, then go to step2
Controllability matrix 𝐶𝑀 = |𝐵 𝐴𝐵 𝐴2 𝐵 | ≠ 0
0
𝐵 = [0] ;
1
0 1 0 0 0
𝐴×𝐵 = [ 0 0 1 ] × [0] = [ 1 ]
−1 −5 −6 1 −6
0 1 0 0 1
𝐴2 𝐵 = 𝐴 × 𝐴𝐵 = [ 0 0 1 ] × [ 1 ] = [ −6]
−1 −5 −6 −6 31
0 0 1
1 −6 0 −6 0 1
|𝐶𝑀| = |𝐵 𝐴𝐵 𝐴2 𝐵 | = [0 1 −6] = 0 | | − 0| | + 1| | = −1 ≠ 0
−6 31 1 31 1 −6
1 −6 31
Hence the system is completely state controllable.
Step 3: Substitute this 𝑲 matrix into the desired characteristics polynomial |𝒔𝑰 − 𝑨 + 𝑩𝑲|
and equate it to (𝒔 − 𝝁𝟏 )(𝒔 − 𝝁𝟐 )(𝒔 − 𝝁𝟑 ).
𝑠 0 0 0 1 0 0
[𝑠𝐼 − 𝐴 + 𝐵𝐾] = [0 𝑠 0] − [ 0 0 1 ] + [0] [𝑘1 𝑘2 𝑘3 ]
0 0 𝑠 −1 −5 −6 1
𝑠 −1 0 0 0 0 𝑠 −1 0
= [0 𝑠 −1 ] + [ 0 0 0]=[ 0 𝑠 −1 ]
1 5 𝑠+6 𝑘1 𝑘2 𝑘3 1 + 𝑘1 5 + 𝑘2 𝑠 + 6 + 𝑘2
𝑠 −1 0
| 0 𝑠 −1 | = 𝑠 3 + (6 + 𝑘3 )𝑠 2 + (5 + 𝑘2 )𝑠 + 1 + 𝑘1
1 + 𝑘1 5 + 𝑘2 𝑠 + 6 + 𝑘2
Step 4: By equating the coefficients of like power of ′𝒔′ on both sides, it is possible to determine
the values of 𝑲𝟏 , 𝑲𝟐 , 𝑲𝟑
6 + 𝑘3 = 14; 𝑘3 = 14 − 6 = 8
[Link]’s College of Engineering 18 VI Semester ICE
IC6601 Advanced Control System Department of ICE 2018-19
5 + 𝑘2 = 60; 𝑘2 = 60 − 5 = 55
Hence 𝐾 = [𝑘1 𝑘2
𝑘3 ] = [199 55 8]
𝑥1
𝑢 = −𝐾𝑥 = −[𝐾1 𝐾2 𝐾3 ] [𝑥2 ]
𝑥3
𝑥1
= −[199 55 8] [𝑥2 ]
𝑥3
𝑢 = −199𝑥1 − 55𝑥2 − 8𝑥3
Step 1: Check the controllability of the given system. If it is controllable, then go to step 2
Controllability matrix 𝐶𝑀 = |𝐵 𝐴𝐵 𝐴2 𝐵 | ≠ 0
0
𝐵 = [0] ;
1
0 1 0 0 0
𝐴×𝐵 = [ 0 0 1 ] × [0] = [ 1 ]
−1 −5 −6 1 −6
0 1 0 0 1
𝐴2 𝐵 = 𝐴 × 𝐴𝐵 = [ 0 0 1 ] × [ 1 ] = [ −6]
−1 −5 −6 −6 31
0 0 1
1 −6 0 −6 0 1
|𝐶𝑀| = |𝐵 𝐴𝐵 𝐴2 𝐵 | = [0 1 −6] = 0 | | − 0| | + 1| | = −1 ≠ 0
−6 31 1 31 1 −6
1 −6 31
Hence the system is completely state controllable.
Step 2: Determine the desired characteristic equation from the specified closed loop poles,
𝝁𝟏 , 𝝁𝟐 , 𝝁𝟑 , . . , 𝝁𝒏
Step 3: Determine the matrix ∅(𝑨) using the coefficient of desired characteristics polynomial.
0 1 0
𝐴=[ 0 0 1]
−1 −5 −6
0 1 0 0 1 0 0 0 1
𝐴2 = 𝐴 × 𝐴 = [ 0 0 1 ] × [ 0 0 1 ] = [ −1 −5 −6]
−1 −5 −6 −1 −5 −6 6 29 31
0 0 1 0 1 0 −1 −5 −6
𝐴3 = 𝐴2 × 𝐴 = [−1 −5 −6] × [ 0 0 1 ]=[ 6 29 31 ]
6 29 31 −1 −5 −6 −31 −149 −157
∅(𝑨) = 𝑨𝟑 + 𝜶𝟏 𝑨𝟐 + 𝜶𝟐 𝑨 + 𝜶𝟑 𝑰
0 1 0 3 0 1 0 2 0 1 0 1 0 0
=[ 0 0 1 ] + 14 [ 0 0 1 ] + 60 [ 0 0 1 ] + 200 [0 1 0]
−1 −5 −6 −1 −5 −6 −1 −5 −6 0 0 1
−1 −5 −6 0 0 1 0 1 0 1 0 0
=[ 6 29 31 ] + 14 [ −1 −5 −6] + 60 [ 0 0 1 ] + 200 [ 0 1 0]
−31 −149 −157 6 29 31 −1 −5 −6 0 0 1
199 55 8
= [ −8 159 7 ]
−7 −43 117
Step 4: Calculate the state feedback gain matrix, ′𝑲′ using the Ackermann’s formula
𝑲 = [𝟎 𝟎 𝟎 , … … . 𝟏] [𝑩 𝑨𝑩 𝑨𝟐 𝑩 , … … . . 𝑨𝒏−𝟏 𝑩]−𝟏 [∅(𝑨)]
𝑲 = [𝟎 𝟎 𝟏] [𝑪𝑴]−𝟏 [∅(𝑨)]
0 0 1 −1 𝐴𝑑𝑗(𝐶𝑀)
𝐶𝑀−1 = [0 1 −6] =
|𝐶𝑀|
1 −6 31
0 0 1
1 −6 0 −6 0 1
|𝐶𝑀| = |𝐵 𝐴𝐵 2 |
𝐴 𝐵 = |0 1 −6| = 0 | | − 0| | +1| | = −1
−6 31 1 31 1 −6
1 −6 31
1 −6 0 −6 0 1 𝑇
+| | −| | +| |
−6 31 1 31 1 −6 −5 −6 −1 𝑇 −5 −6 −1
0 1 0 1 0 0
𝐴𝑑𝑗(𝐶𝑀) = − | | +| | −| | = [−6 −1 0 ] = [−6 −1 0 ]
−6 31 1 31 1 −6
0 1 0 1 0 0 −1 0 0 −1 0 0
[ + |1 −6
| −|
0 −6
| +|
0 1
|]
𝐴𝑑𝑗(𝐶𝑀) −5 −6 −1 5 6 1
𝐶𝑀−1 = = −1 [−6 −1 0 ] = [6 1 0]
|𝐶𝑀|
−1 0 0 1 0 0
5 6 1 199 55 8
𝐾 = [0 0 1] × [ 6 1 0] × [ −8 159 7 ]
1 0 0 −7 −43 117
199 55 8
𝐾 = [1 0 0] × [ −8 159 7 ] = [199 55 8]
−7 −43 117
[Link]’s College of Engineering 20 VI Semester ICE
IC6601 Advanced Control System Department of ICE 2018-19
𝐾 = [𝑘1 𝑘2 𝑘3 ] = [199 55 8]
𝑥1
= −𝐾𝑥 = −[𝐾1 𝐾2 𝐾3 ] [𝑥2 ]
𝑥3
𝑥1
= −[199 55 8] [𝑥2 ]
𝑥3
𝑢 = −199𝑥1 − 55𝑥2 − 8𝑥3
𝟎 𝟎 −𝟔 𝟏
With 𝑨 = [𝟏 𝟎 −𝟏𝟏] ; 𝑩 = [𝟎] ; 𝑪 = [𝟎 𝟎 𝟏]
𝟎 𝟏 −𝟔 𝟎
Compute 𝑲 so that the control law 𝒖 = −𝑲𝒙, places the control loop poles at −𝟐 ±
𝐣 𝟑. 𝟒𝟔𝟒, −𝟑. Give the state variable model of the closed loop system
Solution
0 0 −6 1
Given 𝐴 = [1 0 −11]; 𝐵 = [0]; 𝜇1 = −2 + 𝑗3.464, 𝜇2 = −2 − 𝑗3.464, 𝜇3 = −3; 𝑛 = 3
0 1 −6 0
Step 1: Check the controllability of the given system. If it is controllable, then go to step 2
Controllability matrix 𝐶𝑀 = |𝐵 𝐴𝐵 𝐴2 𝐵 | ≠ 0
1
𝐵 = [0] ;
0
0 0 −6 1 0
𝐴 × 𝐵 = [1 0 −11] × [0] = [1]
0 1 −6 0 0
0 0 −6 0 0
2
𝐴 𝐵 = 𝐴 × 𝐴𝐵 = [1 0 −11] × [1] = [0]
0 1 −6 0 1
1 0 0
1 0 0 0 0 1
|𝐶𝑀| = |𝐵 𝐴𝐵 𝐴 𝐵2| = |0 1 0| = 1 | | − 0| |+ 0| |=1≠0
0 1 0 1 0 0
0 0 1
Hence the system is completely state controllable.
Step 2: Determine the desired characteristic equation from the specified closed loop poles,
𝝁𝟏 , 𝝁𝟐 , 𝝁𝟑 , . . , 𝝁𝒏
i.e., (𝒔 − 𝝁𝟏 ) (𝒔 − 𝝁𝟐 ) (𝒔 − 𝝁𝟑 ). . . . . . . . (𝒔 − 𝝁𝒏 ) = 𝒔𝒏 + 𝜶𝟏 𝒔𝒏−𝟏 +, … … . . , +𝜶𝒏−𝟏 𝒔𝟏 + 𝜶𝒏
and determine the values of 𝜶𝟏 , 𝜶𝟐 , … … . , 𝜶𝒏 .
Step 3: Determine the matrix ∅(𝑨) using the coefficient of desired characteristics polynomial.
∅(𝑨) = 𝑨𝒏 + 𝜶𝟏 𝑨𝒏−𝟏 + 𝜶𝟐 𝑨𝒏−𝟐 +, … … . . , +𝜶𝒏−𝟏 𝑨𝟏 + 𝜶𝒏 𝑰
0 0 −6
𝐴 = [1 0 −11]
0 1 −6
0 0 −6 0 0 −6 0 −6 36
𝐴2 = 𝐴 × 𝐴 = [1 0 −11] × [1 0 −11] = [0 −11 60]
0 1 −6 0 1 −6 1 −6 25
0 −6 36 0 0 −6 −6 36 −150
3 2
𝐴 = 𝐴 × 𝐴 = [0 −11 60] × [1 0 −11] = [−11 60 −239]
1 −6 25 0 1 −6 −6 25 −90
∅(𝑨) = 𝑨𝟑 + 𝜶𝟏 𝑨𝟐 + 𝜶𝟐 𝑨 + 𝜶𝟑 𝑰
0 0 −6 3 0 0 −6 2 0 0 −6 1 0 0
=[1 0 −11] + 7 [1 0 −11] + 28 [1 0 −11] + 48 [0 1 0]
0 1 −6 0 1 −6 0 1 −6 0 0 1
−6 36 −150 0 −6 36 0 0 −6 1 0 0
= [−11 60 −239] + 7 [0 −11 60] + 28 [1 0 −11] + 48 [0 1 0]
−6 25 −90 1 −6 25 0 1 −6 0 0 1
42 −6 −66
= [17 31 −127]
1 11 −35
Step 4: Calculate the state feedback gain matrix, ′𝑲′ using the Ackermann’s formula
𝑲 = [𝟎 𝟎 𝟎 , … … . 𝟏] [𝑩 𝑨𝑩 𝑨𝟐 𝑩 , … … . . 𝑨𝒏−𝟏 𝑩]−𝟏 [∅(𝑨)]
𝑲 = [𝟎 𝟎 𝟏] [𝑪𝑴]−𝟏 [∅(𝑨)]
1 0 0 −1 𝐴𝑑𝑗(𝐶𝑀)
𝐶𝑀−1 = [0 1 0] =
|𝐶𝑀|
0 0 1
1 0 0
1 0 0 0 0 1
|𝐶𝑀| = |𝐵 𝐴𝐵 2
𝐴 𝐵 | = |0 1 0| = 1 | | − 0| |+ 0| |=1
0 1 0 1 0 0
0 0 1
1 0 0 0 0 1 𝑇
+| | −| | +| |
0 1 0 1 0 0 1 0 0𝑇 1 0 0
0 0 1 0 1 0
𝐴𝑑𝑗(𝐶𝑀) = − | | +| | −| | = [0 1 0] = [0 1 0]
0 1 0 1 0 0
0 0 1 0 1 0 0 0 1 0 0 1
[+ |1 0
| −|
0 0
| +|
0 1
|]
𝐴𝑑𝑗(𝐶𝑀) 1 0 0 1 0 0
𝐶𝑀−1 = = 1 [0 1 0 ] = [ 0 1 0]
|𝐶𝑀|
0 0 1 0 0 1
1 0 0 42 −6 −66
𝐾 = [0 0 1] × [0 1 0] × [17 31 −127]
0 0 1 1 11 −35
[Link]’s College of Engineering 22 VI Semester ICE
IC6601 Advanced Control System Department of ICE 2018-19
42 −6 −66
𝐾 = [0 0 1] × [17 31 −127] = [1 11 −35]
1 11 −35
𝐾 = [𝑘1 𝑘2 𝑘3 ] = [1 11 −35]
𝑥1
𝑢 = −𝐾𝑥 = −[𝐾1 𝐾2 𝐾3 ] [𝑥2 ]
𝑥3
𝑥1
= −[1 11 −35] [𝑥2 ]
𝑥3
𝑢 = −1𝑥1 − 11𝑥2 + 35𝑥3
𝐱̇ 𝟏 𝟏 𝟐 𝟏 𝐱𝟏 𝟎
[𝐱̇ 𝟐 ] = [𝟎 𝟏 𝟑] [𝐱 𝟐 ] + [𝟎] 𝐮(𝐭);
𝐱̇ 𝟑 𝟏 𝟏 𝟏 𝐱𝟑 𝟏
Determine the state feedback gain for each state variable to place the poles at −𝟏 + 𝐣, −𝟏 −
𝐣, −𝟑
Solution
1 2 1 0
Given 𝐴 = [0 1 3]; 𝐵 = [0]; 𝜇1 = −1 + 𝑗, 𝜇2 = −1 − 𝑗, 𝜇3 = −3; 𝑜𝑟𝑑𝑒𝑟, 𝑛 = 3
1 1 1 1
Step 1: Check the controllability of the given system. If it is controllable, then go to step 2
Controllability matrix 𝐶𝑀 = |𝐵 𝐴𝐵 𝐴2 𝐵 | ≠ 0
0
𝐵 = [0] ;
1
1 2 1 0 1
𝐴 × 𝐵 = [0 1 3] × [0] = [3]
1 1 1 1 1
1 2 1 1 8
𝐴2 𝐵 = 𝐴 × 𝐴𝐵 = [0 1 3] × [ 3] = [ 6]
1 1 1 1 5
0 1 8
3 6 0 6 0 3
|𝐶𝑀| = |𝐵 𝐴𝐵 𝐴2 𝐵 | = |0 3 6| = 0 | | − 1| |+ 8| | = −18 ≠ 0
1 5 1 5 1 1
1 1 5
[Link]’s College of Engineering 23 VI Semester ICE
IC6601 Advanced Control System Department of ICE 2018-19
Hence the system is completely state controllable.
Step 2: Determine the desired characteristic equation from the specified closed loop poles,
𝝁𝟏 , 𝝁𝟐 , 𝝁𝟑 , . . , 𝝁𝒏
i.e., (𝒔 − 𝝁𝟏 ) (𝒔 − 𝝁𝟐 ) (𝒔 − 𝝁𝟑 ). . . . . . . . (𝒔 − 𝝁𝒏 ) = 𝒔𝒏 + 𝜶𝟏 𝒔𝒏−𝟏 +, … … . . , +𝜶𝒏−𝟏 𝒔𝟏 + 𝜶𝒏
and determine the values of 𝜶𝟏 , 𝜶𝟐 , … … . , 𝜶𝒏 .
(𝑠 + 1 − 𝑗) (𝑠 + 1 + 𝑗) (𝑠 + 3) = (𝑠 2 + 2𝑠 + 2) (𝑠 + 3) = 0
= 𝑠 3 + 5𝑠 2 + 8𝑠 + 6 = 0
= 𝑠 3 + 𝛼1 𝑠 2 + 𝛼2 𝑠 + 𝛼3 = 0
Hence, 𝛼1 = 5; 𝛼2 = 8; 𝛼3 = 6
Step 3: Determine the matrix ∅(𝑨) using the coefficient of desired characteristics polynomial.
∅(𝑨) = 𝑨𝒏 + 𝜶𝟏 𝑨𝒏−𝟏 + 𝜶𝟐 𝑨𝒏−𝟐 +, … … . . , +𝜶𝒏−𝟏 𝑨𝟏 + 𝜶𝒏 𝑰
1 2 1
𝐴 = [0 1 3]
1 1 1
1 2 1 1 2 1 2 5 8
2
𝐴 = 𝐴 × 𝐴 = [0 1 3 ] × [0 1 3] = [3 4 6]
1 1 1 1 1 1 2 4 5
2 5 8 1 2 1 10 17 25
𝐴3 = 𝐴2 × 𝐴 = [3 4 6] × [0 1 3] = [ 9 16 21]
2 4 5 1 1 1 7 13 19
∅(𝑨) = 𝑨𝟑 + 𝜶𝟏 𝑨𝟐 + 𝜶𝟐 𝑨 + 𝜶𝟑 𝑰
1 2 13 1 2 12 1 2 1 1 0 0
=[0 1 3] + 5 [0 1 3] + 8 [0 1 3] + 6 [0 1 0]
1 1 1 1 1 1 1 1 1 0 0 1
10 17 25 2 5 8 1 2 1 1 0 0 34 58 73
=[9 16 21] + 5 [3 4 6] + 8 [0 1 3] + 6 [0 1 0] = [24 50 75]
7 13 19 2 4 5 1 1 1 0 0 1 25 41 58
Step 4: Calculate the state feedback gain matrix, ′𝑲′ using the Ackermann’s formula
𝑲 = [𝟎 𝟎 𝟎 , … … . 𝟏] [𝑩 𝑨𝑩 𝑨𝟐 𝑩 , … … . . 𝑨𝒏−𝟏 𝑩]−𝟏 [∅(𝑨)]
𝑲 = [𝟎 𝟎 𝟏] [𝑪𝑴]−𝟏 [∅(𝑨)]
0 1 8 −1 𝐴𝑑𝑗(𝐶𝑀)
𝐶𝑀−1 = [0 3 6] =
|𝐶𝑀|
1 1 5
0 1 8
3 6 0 6 0 3
|𝐶𝑀| = |𝐵 𝐴𝐵 2
𝐴 𝐵 | = [ 0 3 6] = 0 | | − 1| | + 8| | = −18
1 5 1 5 1 1
1 1 5
−1
𝐴𝑑𝑗(𝐶𝑀) 1 9 3 −18 −0.5 −0.1667 1
𝐶𝑀 == = − [ 6 −8 0 ] = [−0.3333 0.4444 0]
|𝐶𝑀| 18
−3 2 0 0.1667 −0.0556 0
−0.5 −0.1667 1 34 58 73
𝐾 = [0 0 1] × [−0.3333 0.4444 0] × [24 50 75]
0.1667 −0.0556 0 25 41 58
34 58 73
𝐾 = [0.1667 −0.0556 0] × [24 50 75] = [4.3334 6.8886 7.9991]
25 41 58
𝐾 = [𝑘1 𝑘2 𝑘3 ] = [4.3334 6.8886 7.9991]
−𝟎. 𝟏𝟒 𝟎. 𝟑𝟑 −𝟎. 𝟑𝟑 𝟎
𝑿̇ = [ 𝟎. 𝟏 −𝟎. 𝟐𝟖 𝟎 ] 𝑿 + [ 𝟎 ] 𝑼; 𝒀 = [𝟐 𝟎 𝟎] 𝑿
𝟎 𝟏. 𝟕 −𝟎. 𝟕𝟕 −𝟎. 𝟎𝟐𝟓
Design a state feedback controller which will give closed loop poles at −𝟎. 𝟔𝟕, −𝟎. 𝟔𝟕 & −
𝟎. 𝟔𝟕
Solution
Step 1: Check the controllability of the given system. If it is controllable, then go to step 2
Controllability matrix 𝐶𝑀 = |𝐵 𝐴𝐵 𝐴2 𝐵 | ≠ 0
0
𝐵=[ 0 ]
−0.025
−0.14 0.33 −0.33 0 0.0083
𝐴 × 𝐵 = [ 0.1 −0.28 0 ] ×[ 0 ] =[ 0 ]
0 1.7 −0.77 −0.025 0.0193
−0.14 0.33 −0.33 0.0083 −0.0075
𝐴2 𝐵 = 𝐴 × 𝐴𝐵 = [ 0.1 −0.28 0 ] × [ 0 ] = [ 0.0008 ]
0 1.7 −0.77 0.0193 −0.0148
0 0.0083 −0.0075
|𝐶𝑀| = |𝐵 𝐴𝐵 𝐴 𝐵| = | 0
2 0 0.0008 |
−0.025 0.0193 −0.0148
0 0.0008 0 0.0008 0 0
= 0| | − 0.0083 | | + −0.0075 | |
0.0193 −0.0148 −0.025 −0.0148 −0.025 0.0193
= −1.66 × 10−7 ≠ 0
Step 2: Determine the desired characteristic equation from the specified closed loop poles,
𝝁𝟏 , 𝝁𝟐 , 𝝁𝟑 , . . , 𝝁𝒏
i.e., (𝒔 − 𝝁𝟏 ) (𝒔 − 𝝁𝟐 ) (𝒔 − 𝝁𝟑 ). . . . . . . . (𝒔 − 𝝁𝒏 ) = 𝒔𝒏 + 𝜶𝟏 𝒔𝒏−𝟏 +, … … . . , +𝜶𝒏−𝟏 𝒔𝟏 + 𝜶𝒏
and determine the values of 𝜶𝟏 , 𝜶𝟐 , … … . , 𝜶𝒏 .
The desired characteristic equation is
(𝑠 + 0.67) (𝑠 + 0.67) (𝑠 + 0.67) = 𝑠 3 + 2.01𝑠 2 + 1.3467𝑠 + 0.3008 = 0
= 𝑠 3 + 𝛼1 𝑠 2 + 𝛼2 𝑠 + 𝛼3 = 0
Step 3: Determine the matrix ∅(𝑨) using the coefficient of desired characteristics polynomial.
∅(𝑨) = 𝑨𝒏 + 𝜶𝟏 𝑨𝒏−𝟏 + 𝜶𝟐 𝑨𝒏−𝟐 +, … … . . , +𝜶𝒏−𝟏 𝑨𝟏 + 𝜶𝒏 𝑰
∅(𝑨) = 𝑨𝟑 + 𝜶𝟏 𝑨𝟐 + 𝜶𝟐 𝑨 + 𝜶𝟑 𝑰
𝑲 = [𝟎 𝟎 𝟏] [𝑪𝑴]−𝟏 [∅(𝑨)]
Solution
1 1 −2 1
Given, 𝐴 = [0 1 1 ] ; 𝐵 = [0] ; 𝜇1 = −2, 𝜇2 = −1 + 𝑗, 𝜇3 = −1 − 𝑗; 𝑛 = 3
0 0 1 1
Step 1: Check the controllability of the given system. If it is controllable, then go to step 2
Controllability matrix 𝐶𝑀 = |𝐵 𝐴𝐵 𝐴2 𝐵 | ≠ 0
1
𝐵 = [0]
1
1 1 −2 1 −1
𝐴 × 𝐵 = [0 1 1 ] × [0 ] = [ 1 ]
0 0 1 1 1
1 1 −2 −1 −2
2
𝐴 𝐵 = 𝐴 × 𝐴𝐵 = [0 1 1 ] ×[ 1 ]=[ 2 ]
0 0 1 1 1
1 −1 −2
|𝐶𝑀| = |𝐵 | 1 2 0 2 0 1
𝐴𝐵 2
𝐴 𝐵 = |0 1 2 | = 1| | − (−1) | | − 2| | = −1 ≠ 0
1 1 1 1 1 1
1 1 1
Hence the system is completely state controllable.
Step 2: Determine the desired characteristic equation from the specified closed loop poles,
𝝁𝟏 , 𝝁𝟐 , 𝝁𝟑 , . . , 𝝁𝒏
i.e., (𝒔 − 𝝁𝟏 ) (𝒔 − 𝝁𝟐 ) (𝒔 − 𝝁𝟑 ). . . . . . . . (𝒔 − 𝝁𝒏 ) = 𝒔𝒏 + 𝜶𝟏 𝒔𝒏−𝟏 +, … … . . , +𝜶𝒏−𝟏 𝒔𝟏 + 𝜶𝒏
and determine the values of 𝜶𝟏 , 𝜶𝟐 , … … . , 𝜶𝒏 .
The desired characteristic equation is
(𝑠 + 2) (𝑠 + 1 + 𝑗) (𝑠 + 1 − 𝑗) = 0
= 𝑠 3 + 4𝑠 2 + 6𝑠 + 4 = 0
= 𝑠 3 + 𝛼1 𝑠 2 + 𝛼2 𝑠 + 𝛼3 = 0
Hence, 𝛼1 = 4; 𝛼2 = 6; 𝛼3 = 4
Step 3: Determine the matrix ∅(𝑨) using the coefficient of desired characteristics polynomial.
∅(𝑨) = 𝑨𝒏 + 𝜶𝟏 𝑨𝒏−𝟏 + 𝜶𝟐 𝑨𝒏−𝟐 +, … … . . , +𝜶𝒏−𝟏 𝑨𝟏 + 𝜶𝒏 𝑰
∅(𝑨) = 𝑨𝟑 + 𝜶𝟏 𝑨𝟐 + 𝜶𝟐 𝑨 + 𝜶𝟑 𝑰
1 1 −2
𝐴 = [0 1 1]
0 0 1
1 1 −2 1 1 −2 1 2 −3
𝐴2 = 𝐴 × 𝐴 = [0 1 1 ] × [0 1 1 ] = [0 1 2 ]
0 0 1 0 0 1 0 0 1
1 2 −3 1 1 −2 1 3 −3
𝐴3 = 𝐴2 × 𝐴 = [0 1 2 ] × [0 1 1 ] = [0 1 3]
0 0 1 0 0 1 0 0 1
∅(𝑨) = 𝑨𝟑 + 𝜶𝟏 𝑨𝟐 + 𝜶𝟐 𝑨 + 𝜶𝟑 𝑰
1 1 −2 3 1 1 −2 2 1 1 −2 1 0 0
= [0 1 1 ] + 4 [0 1 1 ] + 6 [0 1 1 ] + 4 [0 1 0]
0 0 1 0 0 1 0 0 1 0 0 1
1 3 −3 1 2 −3 1 1 −2 1 0 0 15 17 −27
=[0 1 3 ] + 4 [0 1 2 ] + 6 [0 1 1 ] + 4 [0 1 0] = [ 0 15 17 ]
0 0 1 0 0 1 0 0 1 0 0 1 0 0 15
Step 4: Calculate the state feedback gain matrix, ′𝑲′ using the Ackermann’s formula
𝑲 = [𝟎 𝟎 𝟎 , … … . 𝟏] [𝑩 𝑨𝑩 𝑨𝟐 𝑩 , … … . . 𝑨𝒏−𝟏 𝑩]−𝟏 [∅(𝑨)]
𝑲 = [𝟎 𝟎 𝟏] [𝑪𝑴]−𝟏 [∅(𝑨)]
1 −1 −2 −1 𝐴𝑑𝑗(𝐶𝑀)
𝐶𝑀−1 = [0 1 2] =
|𝐶𝑀|
1 1 1
1 −1 −2
|𝐶𝑀| = |0 1 1 2 0 2 0 1
2 | = 1| | − (−1) | | − 2| | = −1
1 1 1 1 1 1
1 1 1
1 2 0 2 0 1 𝑇
+| | −| | +| |
1 1 1 1 1 1 −1 2 −1 𝑇 −1 −1 0
−1 −2 1 −2 1 −1
𝐴𝑑𝑗(𝐶𝑀) = − | | +| | −| | = [−1 3 −2] = [ 2 3 −2]
1 1 1 1 1 1
−1 −2 1 −2 1 −1 0 −2 1 −1 −2 1
[+ | 1 2
| −|
0 2
| +|
0 1
|]
𝐴𝑑𝑗(𝐶𝑀) −1 −1 0 1 1 0
𝐶𝑀−1 == = −1 [ 2 3 −2 ] = [ −2 −3 2 ]
|𝐶𝑀|
−1 −2 1 1 2 −1
1 1 0 15 17 −27
𝐾 = [0 0 1] × [−2 −3 2 ] × [ 0 15 17 ]
1 2 −1 0 0 15
15 17 −27
𝐾 = [1 2 ]
−1 × [ 0 15 17 ] = [12 47 −8]
0 0 15
𝐾 = [𝑘1 𝑘2 𝑘3 ] = [12 47 −8]
−𝟏 𝟎 𝟎 𝟎
̇𝑿 = [ 𝟎 −𝟐 𝟎 ] 𝑿 + [𝟎] 𝑼; 𝒀 = [𝟏 𝟎 𝟎]𝑿
𝟎 𝟎 −𝟑 𝟏
Solution
−1 0 0 0
Given, 𝐴 = [ 0 −2 0 ] ; 𝐵 = [0] ; 𝜇1 = −6, 𝜇2 = −1 + 2𝑗, 𝜇3 = −1 − 2𝑗; 𝑛 = 3
0 0 −3 1
Step 1: Check the controllability of the given system. If it is controllable, then go to step 2
Controllability matrix 𝐶𝑀 = |𝐵 𝐴𝐵 𝐴2 𝐵 | ≠ 0
0
𝐵 = [0]
1
−1 0 0 0 0
𝐴 × 𝐵 = [ 0 −2 0 ] × [0] = [ 0 ]
0 0 −3 1 −3
−1 0 0 0 0
2
𝐴 𝐵 = 𝐴 × 𝐴𝐵 = [ 0 −2 0 ] × [ 0 ] = [0]
0 0 −3 −3 9
0 0 0
0 0 0 0 0 0
|𝐶𝑀| = |𝐵 𝐴𝐵 2
𝐴 𝐵 | = |0 0 0| = 0 | | − 0| | + 0| |=0
−3 9 1 9 1 −3
1 −3 9
Hence the system is NOT state controllable. Design of state feedback controller is not possible.
𝐱̇ 𝟐 𝟏 𝐱𝟏 𝟏 𝐱𝟏
7. Given [ 𝟏 ] = [ ] [𝐱 ] + [ ] 𝐮; 𝐲 = [𝟏 𝟏] [𝐱 ]. Find the state feedback gain K, so
𝐱̇ 𝟐 −𝟏 𝟏 𝟐 𝟐 𝟐
that the state feedback system has -1 and -2 as its eigenvalues.
Solution
2 1 1
Given, 𝐴 = [ ] ; 𝐵 = [ ] ; 𝜇1 = −1, 𝜇2 = −2
−1 1 2
Step 1: Check the controllability of the given system. If it is controllable, then go to step 2
Controllability matrix 𝐶𝑀 = |𝐵 𝐴𝐵 | ≠ 0
1
𝐵=[ ]
2
2 1 1 4
𝐴×𝐵 = [ ] ×[ ]= [ ]
−1 1 2 1
|𝐶𝑀| = |𝐵 1 4
𝐴𝐵 | = [ ] = −7 ≠ 0
2 1
Hence the system is completely state controllable.
[Link]’s College of Engineering 30 VI Semester ICE
IC6601 Advanced Control System Department of ICE 2018-19
Step 2: Determine the desired characteristic equation from the specified closed loop poles,
𝝁𝟏 , 𝝁𝟐 , 𝝁𝟑 , . . , 𝝁𝒏
i.e., (𝒔 − 𝝁𝟏 ) (𝒔 − 𝝁𝟐 ) (𝒔 − 𝝁𝟑 ). . . . . . . . (𝒔 − 𝝁𝒏 ) = 𝒔𝒏 + 𝜶𝟏 𝒔𝒏−𝟏 +, … … . . , +𝜶𝒏−𝟏 𝒔𝟏 + 𝜶𝒏
and determine the values of 𝜶𝟏 , 𝜶𝟐 , … … . , 𝜶𝒏 .
The desired characteristic equation is
(𝑠 + 1) (𝑠 + 2) = 0
= 𝑠 2 + 3𝑠 + 2 = 0
= 𝑠 2 + 𝛼1 𝑠 + 𝛼2 = 0
Hence, 𝛼1 = 3; 𝛼2 = 2
Step 3: Determine the matrix ∅(𝑨) using the coefficient of desired characteristics polynomial.
∅(𝑨) = 𝑨𝒏 + 𝜶𝟏 𝑨𝒏−𝟏 + 𝜶𝟐 𝑨𝒏−𝟐 +, … … . . , +𝜶𝒏−𝟏 𝑨𝟏 + 𝜶𝒏 𝑰
∅(𝑨) = 𝑨𝟐 + 𝜶𝟏 𝑨 + 𝜶𝟐 𝑰
2 1
𝐴=[ ]
−1 1
2 1 2 1 3 3
𝐴2 = 𝐴 × 𝐴 = [ ]×[ ]=[ ]
−1 1 −1 1 −3 0
∅(𝑨) = 𝑨𝟐 + 𝜶𝟏 𝑨 + 𝜶𝟐 𝑰
2 12 2 1 1 0
=[ ] + 3[ ] + 2[ ]
−1 1 −1 1 0 1
3 3 2 1 1 0 11 6
=[ ] + 3[ ]+ 2[ ]=[ ]
−3 0 −1 1 0 1 −6 5
Step 4: Calculate the state feedback gain matrix, ′𝑲′ using the Ackermann’s formula
𝑲 = [𝟎 𝟎 𝟎 , … … . 𝟏] [𝑩 𝑨𝑩 𝑨𝟐 𝑩 , … … . . 𝑨𝒏−𝟏 𝑩]−𝟏 [∅(𝑨)]
𝑲 = [𝟎 𝟏] [𝑪𝑴]−𝟏 [∅(𝑨)]
1 4 −1 𝐴𝑑𝑗(𝐶𝑀)
𝐶𝑀−1 = [ ] =
2 1 |𝐶𝑀|
|𝐶𝑀| = |𝐵 1 4
𝐴𝐵 | = | | = −7
2 1
1 4 1 −4
𝐴𝑑𝑗(𝐶𝑀) = 𝑎𝑑𝑗 [ ]=[ ]
2 1 −2 1
𝐴𝑑𝑗(𝐶𝑀) 1 1 −4 −0.1429 0.5714
𝐶𝑀−1 = =− [ ]=[ ]
|𝐶𝑀| 7 −2 1 0.2857 −0.1429
−0.1429 0.5714 11 6
𝐾 = [0 1] × [ ]×[ ]
0.2857 −0.1429 −6 5
11 6
𝐾 = [0.2857 −0.1429] × [ ] = [4 1]
−6 5
𝐾 = [𝑘1 𝑘2 ] = [4 1]
|𝑺𝑰 − 𝑨| = | 𝟎 𝟏
| = 𝑺𝟐 − 𝟐𝟎. 𝟔 = 𝟎
𝟐𝟎. 𝟔 𝒔
Since the characteristic roots are 𝑺 = ±𝟒. 𝟓𝟑𝟗, the system is unstable. By using the state
feedback control 𝒖 = −𝑲𝒙, it is desired to have closed loop poles at 𝑺 = −𝟏. 𝟖 ± 𝒋𝟐. 𝟒 ( that is
the eigenvalues of 𝑨 − 𝑩𝑲 to be 𝝁𝟏 = −𝟏. 𝟖 + 𝒋𝟐. 𝟒 and 𝝁𝟐 = −𝟏. 𝟖 − 𝒋𝟐. 𝟒. determine the
state feedback gain matrix K
Solution
𝟎 𝟏 𝟏
Given, 𝐴 = [ ] ; 𝐵 = [ ] ; 𝝁𝟏 = −𝟏. 𝟖 + 𝒋𝟐. 𝟒 ; 𝝁𝟐 = −𝟏. 𝟖 − 𝒋𝟐. 𝟒
𝟐𝟎. 𝟔 𝟎 𝟎
Step 1: Check the controllability of the given system. If it is controllable, then go to step 2
Controllability matrix 𝐶𝑀 = |𝐵 𝐴𝐵 | ≠ 0
1
𝐵=[ ]
0
0 1 1 0
𝐴×𝐵 = [ ] ×[ ]=[ ]
20.6 0 0 20.6
1 0
|𝐶𝑀| = |𝐵 𝐴𝐵 | = [ ] = 20.6 ≠ 0
0 20.6
Hence the system is completely state controllable.
Step 2: Determine the desired characteristic equation from the specified closed loop poles,
𝝁𝟏 , 𝝁𝟐 , 𝝁𝟑 , . . , 𝝁𝒏
i.e., (𝒔 − 𝝁𝟏 ) (𝒔 − 𝝁𝟐 ) (𝒔 − 𝝁𝟑 ). . . . . . . . (𝒔 − 𝝁𝒏 ) = 𝒔𝒏 + 𝜶𝟏 𝒔𝒏−𝟏 +, … … . . , +𝜶𝒏−𝟏 𝒔𝟏 + 𝜶𝒏
and determine the values of 𝜶𝟏 , 𝜶𝟐 , … … . , 𝜶𝒏 .
The desired characteristic equation is
(𝑠 + 1.8 + 𝑗2.4) (𝑠 + 1.8 − 𝑗2.4) = 0
= 𝑠 2 + 3.6𝑠 + 9 = 0
= 𝑠 2 + 𝛼1 𝑠 + 𝛼2 = 0
Hence, 𝛼1 = 3.6; 𝛼2 = 9
Step 3: Determine the matrix ∅(𝑨) using the coefficient of desired characteristics polynomial.
∅(𝑨) = 𝑨𝒏 + 𝜶𝟏 𝑨𝒏−𝟏 + 𝜶𝟐 𝑨𝒏−𝟐 +, … … . . , +𝜶𝒏−𝟏 𝑨𝟏 + 𝜶𝒏 𝑰
∅(𝑨) = 𝑨𝟐 + 𝜶𝟏 𝑨 + 𝜶𝟐 𝑰
0 1
𝐴=[ ]
20.6 0
0 1 0 1 20.6 0
𝐴2 = 𝐴 × 𝐴 = [ ]×[ ]=[ ]
20.6 0 20.6 0 0 20.6
∅(𝑨) = 𝑨𝟐 + 𝜶𝟏 𝑨 + 𝜶𝟐 𝑰
0 12 0 1 1 0
=[ ] + 3.6 [ ] + 9[ ]
20.6 0 20.6 0 0 1
[Link]’s College of Engineering 32 VI Semester ICE
IC6601 Advanced Control System Department of ICE 2018-19
20.6 0 0 1 1 0 29.6 3.6
=[ ] + 3.6 [ ] + 9[ ]=[ ]
0 20.6 20.6 0 0 1 74.16 29.6
Step 4: Calculate the state feedback gain matrix, ′𝑲′ using the Ackermann’s formula
𝑲 = [𝟎 𝟎 𝟎 , … … . 𝟏] [𝑩 𝑨𝑩 𝑨𝟐 𝑩 , … … . . 𝑨𝒏−𝟏 𝑩]−𝟏 [∅(𝑨)]
𝑲 = [𝟎 𝟏] [𝑪𝑴]−𝟏 [∅(𝑨)]
1 0 −1 𝐴𝑑𝑗(𝐶𝑀)
𝐶𝑀−1 = [ ] =
0 20.6 |𝐶𝑀|
1 0
|𝐶𝑀| = |𝐵 𝐴𝐵 | = | | = 20.6
0 20.6
1 0 20.6 0
𝐴𝑑𝑗(𝐶𝑀) = 𝑎𝑑𝑗 [ ]=[ ]
0 20.6 0 1
𝐴𝑑𝑗(𝐶𝑀) 1 20.6 0 1 0
𝐶𝑀−1 = = [ ]=[ ]
|𝐶𝑀| 20.6 0 1 0 0.0485
1 0 29.6 3.6
𝐾 = [0 1] × [ ]×[ ]
0 0.0485 74.16 29.6
29.6 3.6
𝐾 = [0 0.0485] × [ ] = [3.6 1.4369]
74.16 29.6
𝐾 = [𝐾1 𝐾2 ] = [3.6 1.4369]
𝑥1
𝑢 = −𝐾𝑥 = −[𝐾1 𝐾2 ] [𝑥 ]
2
𝑥1
= −[3.6 1.4369] [𝑥 ]
2
𝑢 = −3.6𝑥1 − 1.4369𝑥2
𝒙̇ 𝟏 −𝟏 𝟎 𝟎 𝒙𝟏 𝟏𝟎
[𝒙̇ 𝟐 ] = = [ 𝟏 −𝟐 𝟎 ] [𝒙𝟐 ] + [ 𝟏 ] 𝑼
𝒙̇ 𝟑 𝟐 𝟏 −𝟑 𝒙𝟑 𝟎
Design a state feedback controller which will give closed loop poles at −𝟏 + 𝒋𝟐, −𝟏 − 𝒋𝟐, 𝟐 (May
2017)
Step 1: Check the controllability of the given system. If it is controllable, then go to step2
Note: If the system is of low order (𝑛 ≤ 3), direct substitution of matrix K into the desired
characteristics polynomial may be simpler.
Controllability matrix 𝐶𝑀 = |𝐵 𝐴𝐵 𝐴2 𝐵| ≠ 0
10
𝐵 = [ 1 ];
0
−1 0 0 10 −10
𝐴 × 𝐵 = [ 1 −2 0 ] × [ 1 ] = [ 8 ]
2 1 −3 1 21
−1 0 0 −10 10
𝐴2 𝐵 = 𝐴 × 𝐴𝐵 = [ 1 −2 0 ] × [ 8 ] = [−26]
2 1 −3 21 −75
10 −10 10
|𝐶𝑀| = |𝐵 𝐴𝐵 2
𝐴 𝐵 | = [ 1 8 −26]
0 21 −75
8 −26 1 −26 1 8
= 10 | | − (−10) | | + 10 | | = −1080 ≠ 0
21 −75 0 −75 0 21
Step 2: Determine the desired characteristic equation from the specified closed loop poles, 𝜇1 , 𝜇2 , 𝜇3 , . . , 𝜇𝑛
= 𝑠 3 + 𝛼1 𝑠 2 + 𝛼2 𝑠 + 𝛼3 = 0
Hence, 𝛼1 = 0; 𝛼2 = 1; 𝛼3 = −10
Step 3: Determine the matrix ∅(𝐴) using the coefficient of desired characteristics polynomial.
−1 0 0
𝐴 = [ 1 −2 0 ]
2 1 −3
−1 0 0 −1 0 0 1 0 0
𝐴2 = 𝐴 × 𝐴 = [ 1 −2 0 ] × [ 1 −2 0 ] = [−3 4 0]
2 1 −3 2 1 −3 −7 −5 9
1 0 0 −1 0 0 −1 0 0
𝐴3 = 𝐴2 × 𝐴 = [−3 4 0] × [ 1 −2 0 ] = [ 7 −8 0 ]
−7 −5 9 2 1 −3 20 −19 −27
∅(𝑨) = 𝑨𝟑 + 𝜶𝟏 𝑨𝟐 + 𝜶𝟐 𝑨 + 𝜶𝟑 𝑰
−1 0 0 3 −1 0 0 2 −1 0 0 1 0 0
=[ 1 −2 0 ] + 0 [ 1 −2 0 ] + 1 [ 1 −2 0 ] − 10 [0 1 0]
2 1 −3 2 1 −3 2 1 −3 0 0 1
−1 0 0 1 0 0 −1 0 0 1 0 0
=[ 7 −8 0 ] + 0 [−3 4 0] + 1 [ 1 −2 0 ] − 10 [0 1 0]
20 −19 −27 −7 −5 9 2 1 −3 0 0 1
−12 0 0
=[ 8 −20 0 ]
20 19 −27
Step 4: Calculate the state feedback gain matrix, ′𝐾′ using the Ackermann’s formula
𝐾 = [0 0 0 , … … . 1] [𝐵 𝐴𝐵 𝐴2 𝐵 , … … . . 𝐴𝑛−1 𝐵]−1 [∅(𝐴)]
𝐾 = [0 0 1] [𝐶𝑀]−1 [∅(𝐴)]
[Link]’s College of Engineering 34 VI Semester ICE
IC6601 Advanced Control System Department of ICE 2018-19
10 −10 10 −1 𝐴𝑑𝑗(𝐶𝑀)
−1
𝐶𝑀 =[1 8 −26] =
|𝐶𝑀|
0 21 −75
8 −26 1 −26 1 8
|𝐶𝑀| = |𝐵 𝐴2 𝐵| == 10 |21 −75| − (−10) |0 −75| + 10 |0 21| = −1080
𝐴𝐵
𝑇
8 −26 1 −26 1 8
+| | −| | +| |
21 −75 0 −75 0 21 −54 75 21 𝑇
−10 10 10 10 10 −10
𝐴𝑑𝑗(𝐶𝑀) = − | | +| | −| | = [−540 −750 −210]
21 −75 0 −75 0 21
−10 10 10 10 10 −10 180 270 90
[+ | 8 −26
| −|
1 −26
| +|
0 21
|]
−54 −540 180
= [ 75 −750 270]
21 −210 90
−54 −540 180 0.05 0.5 −0.1667
𝐴𝑑𝑗(𝐶𝑀) 1
-𝐶𝑀−1 = |𝐶𝑀|
= − 1080 [ 75 −750 270] = [−0.0694 0.6944 −0.25 ]
21 −210 90 −0.0194 0.1944 −0.0833
0.05 0.5 −0.1667 −12 0 0
𝐾 = [0 0 1] × [−0.0694 0.6944 −0.25 ] × [ 8 −20 0 ]
−0.0194 0.1944 −0.0833 20 19 −27
−12 0 0
𝐾 = [−0.0194 0.1944 −0.0833] × [ 8 −20 0 ] = [−0.0444 −5.5556 3.3333]
20 19 −27
𝐾 = [𝑘1 𝑘2
𝑘3 ] = [−0.0444 −5.5556 3.3333]
𝑥1 𝑥1
𝑥 𝑥
𝑢 = −𝐾𝑥 = −[𝐾1 𝐾2 𝐾3 ] [ 2 ] = −[−0.0444 −5.5556 3.3333] [ 2 ]
𝑥3 𝑥3
𝑢 = 0.0444𝑥1 + 5.5556𝑥2 − 3.3333𝑥3
𝟎 𝟏 𝟎 𝟏
𝑿̇ = [ ] 𝑿 + [ ] 𝒖𝟏 + [ ] 𝒖 𝟐
−𝟔 −𝟓 𝟏 𝟏
Design a state feedback controller to place the closed loop poles at −𝟏𝟎 and − 𝟏𝟏 for the three
following cases. (i) Using input one only (ii) Using input two only Case (i): Using input one
only.
Solution
0 1 0
Given, 𝐴 = [ ] ; 𝐵 = [ ] ; 𝜇1 = −1.8 + 𝑗2.4 ; 𝜇2 = −1.8 − 𝑗2.4
−6 −5 1
Step 1: Check the controllability of the given system. If it is controllable, then go to step 2
Controllability matrix 𝐶𝑀 = |𝐵 𝐴𝐵 | ≠ 0
0
𝐵=[ ]
1
0 1 0 1
𝐴×𝐵 = [ ] ×[ ]=[ ]
−6 −5 1 −5
0 1
|𝐶𝑀| = |𝐵 𝐴𝐵 | = [ ] = −1 ≠ 0
1 −5
[Link]’s College of Engineering 35 VI Semester ICE
IC6601 Advanced Control System Department of ICE 2018-19
Hence the system is completely state controllable.
Step 2: Determine the desired characteristic equation from the specified closed loop poles,
𝝁𝟏 , 𝝁𝟐 , 𝝁𝟑 , . . , 𝝁𝒏
i.e., (𝒔 − 𝝁𝟏 ) (𝒔 − 𝝁𝟐 ) (𝒔 − 𝝁𝟑 ). . . . . . . . (𝒔 − 𝝁𝒏 ) = 𝒔𝒏 + 𝜶𝟏 𝒔𝒏−𝟏 +, … … . . , +𝜶𝒏−𝟏 𝒔𝟏 + 𝜶𝒏
and determine the values of 𝜶𝟏 , 𝜶𝟐 , … … . , 𝜶𝒏 .
The desired characteristic equation is
(𝑠 + 10) (𝑠 + 11) = 0
= 𝑠 2 + 21𝑠 + 110 = 0
= 𝑠 2 + 𝛼1 𝑠 + 𝛼2 = 0
Step 3: Determine the matrix ∅(𝑨) using the coefficient of desired characteristics polynomial.
∅(𝑨) = 𝑨𝒏 + 𝜶𝟏 𝑨𝒏−𝟏 + 𝜶𝟐 𝑨𝒏−𝟐 +, … … . . , +𝜶𝒏−𝟏 𝑨𝟏 + 𝜶𝒏 𝑰
∅(𝑨) = 𝑨𝟐 + 𝜶𝟏 𝑨 + 𝜶𝟐 𝑰
0 1
𝐴=[ ]
−6 −5
0 1 0 1 −6 −5
𝐴2 = 𝐴 × 𝐴 = [ ]×[ ]=[ ]
−6 −5 −6 −5 30 19
∅(𝑨) = 𝑨𝟐 + 𝜶𝟏 𝑨 + 𝜶𝟐 𝑰
0 1 2 0 1 1 0
=[ ] + 21 [ ] + 110 [ ]
−6 −5 −6 −5 0 1
−6 −5 0 1 1 0 104 16
=[ ] + 21 [ ] + 110 [ ]=[ ]
30 19 −6 −5 0 1 −96 24
Step 4: Calculate the state feedback gain matrix, ′𝑲′ using the Ackermann’s formula
𝑲 = [𝟎 𝟎 𝟎 , … … . 𝟏] [𝑩 𝑨𝑩 𝑨𝟐 𝑩 , … … . . 𝑨𝒏−𝟏 𝑩]−𝟏 [∅(𝑨)]
𝑲 = [𝟎 𝟏] [𝑪𝑴]−𝟏 [∅(𝑨)]
0 1 −1 𝐴𝑑𝑗(𝐶𝑀)
𝐶𝑀−1 = [ ] =
1 −5 |𝐶𝑀|
0 1
|𝐶𝑀| = |𝐵 𝐴𝐵 | = [ ] = −1
1 −5
0 1 −5 −1
𝐴𝑑𝑗(𝐶𝑀) = 𝑎𝑑𝑗 [ ]=[ ]
1 −5 −1 0
𝐴𝑑𝑗(𝐶𝑀) −5 −1 5 1
𝐶𝑀−1 = = −1 [ ]=[ ]
|𝐶𝑀| −1 0 1 0
5 1 104 16
𝐾 = [0 1] × [ ]×[ ]
1 0 −96 24
104 16
𝐾 = [1 0] × [ ] = [104 16]
−96 24
𝐾 = [𝑘1 𝑘2 ] = [104 16]
Solution
0 1 1
Given, 𝐴 = [ ] ; 𝐵 = [ ] ; 𝜇1 = −1.8 + 𝑗2.4 ; 𝜇2 = −1.8 − 𝑗2.4
−6 −5 1
Step 1: Check the controllability of the given system. If it is controllable, then go to step 2
Controllability matrix 𝐶𝑀 = |𝐵 𝐴𝐵 | ≠ 0
1
𝐵=[ ]
1
0 1 1 1
𝐴×𝐵 = [ ] ×[ ]=[ ]
−6 −5 1 −11
|𝐶𝑀| = |𝐵 1 1
𝐴𝐵 | = [ ] = −12 ≠ 0
1 −11
Hence the system is completely state controllable.
Step 2: Determine the desired characteristic equation from the specified closed loop poles,
𝝁𝟏 , 𝝁𝟐 , 𝝁𝟑 , . . , 𝝁𝒏
i.e., (𝒔 − 𝝁𝟏 ) (𝒔 − 𝝁𝟐 ) (𝒔 − 𝝁𝟑 ). . . . . . . . (𝒔 − 𝝁𝒏 ) = 𝒔𝒏 + 𝜶𝟏 𝒔𝒏−𝟏 +, … … . . , +𝜶𝒏−𝟏 𝒔𝟏 + 𝜶𝒏
and determine the values of 𝜶𝟏 , 𝜶𝟐 , … … . , 𝜶𝒏 .
The desired characteristic equation is
(𝑠 + 10) (𝑠 + 11) = 0
= 𝑠 2 + 21𝑠 + 110 = 0
= 𝑠 2 + 𝛼1 𝑠 + 𝛼2 = 0
Step 3: Determine the matrix ∅(𝑨) using the coefficient of desired characteristics polynomial.
∅(𝑨) = 𝑨𝒏 + 𝜶𝟏 𝑨𝒏−𝟏 + 𝜶𝟐 𝑨𝒏−𝟐 +, … … . . , +𝜶𝒏−𝟏 𝑨𝟏 + 𝜶𝒏 𝑰
∅(𝑨) = 𝑨𝟐 + 𝜶𝟏 𝑨 + 𝜶𝟐 𝑰
0 1
𝐴=[ ]
−6 −5
0 1 0 1 −6 −5
𝐴2 = 𝐴 × 𝐴 = [ ]×[ ]=[ ]
−6 −5 −6 −5 30 19
∅(𝑨) = 𝑨𝟐 + 𝜶𝟏 𝑨 + 𝜶𝟐 𝑰
0 1 2 0 1 1 0
=[ ] + 21 [ ] + 110 [ ]
−6 −5 −6 −5 0 1
−6 −5 0 1 1 0 104 16
=[ ] + 21 [ ] + 110 [ ]=[ ]
30 19 −6 −5 0 1 −96 24
Step 4: Calculate the state feedback gain matrix, ′𝑲′ using the Ackermann’s formula
𝑲 = [𝟎 𝟎 𝟎 , … … . 𝟏] [𝑩 𝑨𝑩 𝑨𝟐 𝑩 , … … . . 𝑨𝒏−𝟏 𝑩]−𝟏 [∅(𝑨)]
𝐾 = [0 1] [𝐶𝑀]−1 [∅(𝐴)]
1 1 −1 𝐴𝑑𝑗(𝐶𝑀)
𝐶𝑀−1 = [ ] =
1 −11 |𝐶𝑀|
|𝐶𝑀| = |𝐵 1 1
𝐴𝐵 | = [ ] = −12
1 −11
1 1 −11 −1
𝐴𝑑𝑗(𝐶𝑀) = 𝑎𝑑𝑗 [ ]=[ ]
1 −11 −1 1
𝐴𝑑𝑗(𝐶𝑀) 1 −11 −1 0.9167 0.0833
𝐶𝑀−1 = =− [ ]=[ ]
|𝐶𝑀| 12 −1 1 0.0833 −0.0833
Solution
0 1 0 1
𝑋̇ = [ ] 𝑋 + [ ] 𝑢1 + [ ] 𝑢2
−6 −5 1 1
𝑋̇ 0 1 𝑋1 0 1
[ 1] = [ ] [ ] + [ ] 𝑢1 + [ ] 𝑢2
𝑋̇2 −6 −5 𝑋2 1 1
Ẋ1 = 0 × X1 + 1 × X 2 + 0 × u1 + 1 × u2
Ẋ 2 = −6 × X1 − 5 × X 2 + 1 × u1 + 1 × u2
𝑋̇ 0 1 𝑋1 0 1 𝑢1
[ 1] = [ ][ ] + [ ][ ]
̇
𝑋2 −6 −5 𝑋2 1 1 𝑢2
𝟎 𝟏 𝟎 𝟏
Given, 𝐴 = [ ]; 𝑩 = [ ] ; 𝝁𝟏 = −𝟏𝟎 ; 𝝁𝟐 = −𝟏𝟏
−𝟔 −𝟓 𝟏 𝟏
Use Direct substitution method for multi input systems (Ackerman’s formula cannot be used
Step 1: Check the controllability of the given system. If it is controllable, then go to step 2
0 1
𝐵=[ ]
1 1
0 1 0 1 1 1
𝐴×𝐵 =[ ] ×[ ]=[ ]
−6 −5 1 1 −5 −11
0 1 1 1
|𝐶𝑀| = |𝐵 𝐴𝐵| = | |≠0
1 1 −5 −11
𝐾 𝐾2
𝑲=[ 1 ]
𝐾3 𝐾4
Step 3: Substitute this 𝑲 matrix into the desired characteristics polynomial |𝒔𝑰 − 𝑨 + 𝑩𝑲| and equate
it to (𝒔 − 𝝁𝟏 )(𝒔 − 𝝁𝟐 ).
[𝑠𝐼 − 𝐴 + 𝐵𝐾] = [ 𝑠 0 0 1 0 1 𝐾1 𝐾2
]−[ ]+[ ][ ]
0 𝑠 −6 −5 1 1 𝐾3 𝐾4
𝑠 −1 𝐾3 𝐾4 𝑠 + 𝐾3 −1 + 𝐾4
=[ ]+[ ]=[ ]
6 𝑠+5 𝐾1 + 𝐾3 𝐾2 + 𝐾4 6 + 𝐾1 + 𝐾3 𝑠 + 5 + 𝐾2 + 𝐾4
𝑠 + 𝐾3 −1 + 𝐾4
| | = (𝑠 + 𝐾3 ) (𝑠 + 5 + 𝐾2 + 𝐾4 ) − (6 + 𝐾1 + 𝐾3 )(−1 + 𝐾4 )
6 + 𝐾1 + 𝐾3 𝑠 + 5 + 𝐾2 + 𝐾4
Step 4: By equating the coefficients of like power of ′𝒔′ on both sides, it is possible to determine the
values of 𝑲𝟏 , 𝑲𝟐 , 𝑲𝟑 , 𝑲𝟒 , 𝑲𝟓 , 𝑲𝟔
5 + 𝐾2 + 𝐾3 + 𝐾4 = 21 (1)
K K2 −16 5
K=[ 1 ]=[ ]
K3 K4 10 1
Solution
1 0 0 1 0
Given 𝐴 = [0 2 0 ] ; 𝐵 = [ 0 1] ; 𝜇1 = −1; 𝜇2 = −2; 𝜇3 = −3
0 0 3 1 1
Step 1: Check the controllability of the given system. If it is controllable, then go to step2
Controllability matrix 𝐶𝑀 = |𝐵 𝐴𝐵 𝐴2 𝐵| ≠ 0
1 0
𝐵 = [ 0 1] ;
1 1
1 0 0 1 0 1 0
𝐴 × 𝐵 = [0 2 0 ] × [ 0 1] = [ 0 2]
0 0 3 1 1 3 3
1 0 0 1 0 1 0
𝐴2 𝐵 = 𝐴 × 𝐴𝐵 = [0 2 0] × [0 2] = [0 4]
0 0 3 3 3 9 9
1 0 1 0 1 0
|𝐶𝑀| = |𝐵 𝐴𝐵 2
𝐴 𝐵 | = |0 1 0 2 0 4| ≠ 0
1 1 3 3 9 9
Hence the system is completely state controllable.
𝐾1 𝐾2 𝐾3
𝑲=[ ]
𝐾4 𝐾5 𝐾6
Step 3: Substitute this 𝑲 matrix into the desired characteristics polynomial |𝒔𝑰 − 𝑨 + 𝑩𝑲| and equate
it to (𝒔 − 𝝁𝟏 )(𝒔 − 𝝁𝟐 )(𝒔 − 𝝁𝟑 ).
𝑠 0 0 1 0 0 1 0 𝐾 𝐾2 𝐾3
[𝑠𝐼 − 𝐴 + 𝐵𝐾] = [0 𝑠 0] − [0 2 0] + [0 1] [𝐾1 𝐾5 𝐾6
]
4
0 0 𝑠 0 0 3 1 1
𝑠−1 0 0 𝑘1 𝑘2 𝑘3
=[ 0 𝑠−2 0 ] + [ 𝑘4 𝑘5 𝑘6 ]
0 0 𝑠−3 𝑘1 + 𝑘4 𝑘2 + 𝑘5 𝑘3 + 𝑘6
Step 4: By equating the coefficients of like power of ′𝒔′ on both sides, it is possible to determine the
values of 𝑲𝟏 , 𝑲𝟐 , 𝑲𝟑
(𝑠 − 1) + 𝑘1 𝑘2 𝑘3
| 𝑘4 (𝑠 − 2) + 𝑘5 𝑘6 | = (𝑠 + 1)(𝑠 + 2)(𝑠 + 3)
𝑘1 + 𝑘4 𝑘2 + 𝑘5 (𝑠 − 3) + 𝑘3 + 𝑘6
𝐱̇ 𝟏 𝟎 𝟐𝟎. 𝟔 𝐱 𝟏 𝟎 𝐱𝟏
12. Consider the system [ ]=[ ] [𝐱 ] + [ ] 𝐮(𝐭); 𝐲 = [𝟎 𝟏] [𝐱 ]. Use observer
𝐱̇ 𝟐 𝟏 𝟎 𝟐 𝟏 𝟐
based approach to design state feedback control such that 𝐮 = −𝐊𝐱̌. Design a full order
state observer assume that the desired eigenvalues of the observer matrix are 𝛃𝟏 =
−𝟏𝟎 𝐚𝐧𝐝 𝛃𝟐 = −𝟏𝟎
Step 1: Check the observability of the given system. If it is observable, then go to step2
𝐶
𝐶𝐴
| 𝐶𝐴2 |
Observability matrix |𝑶𝑴| = ≠𝟎
| ; |
𝐶𝐴𝑛−2
𝐶𝐴𝑛−1
𝐶
Observability matrix |𝑂𝑀| = | | ≠ 0
𝐶𝐴
𝐶 = [0 1];
0 20.6
𝐶 × 𝐴 = [0 1] × [ ] = [1 0]
1 0
|𝑂𝑀| = | 𝐶 | = [0 1] = −1 ≠ 0
𝐶𝐴 1 0
Hence the system is completely state observable.
𝑎1 = 0; 𝑎2 = −20.6
The given system is already in observable canonical form. Hence, the transformation
matrix 𝑄 = (𝑊𝑁 ∗ )−1 is 𝐼
Step 4: Using desired closed loop poles, write the desired polynomial. i.e.
(𝒔 − 𝝁𝟏 )(𝒔 − 𝝁𝟐 )(𝒔 − 𝝁𝟑 ), … … … … … … (𝒔 − 𝝁𝒏 )
𝛼1 = 20; 𝛼2 = 100
1 0
Since the system is in observable canonical form (𝑊𝑁 ∗ )−1 = 𝐼 = [ ]
0 1
1 0 100 − (−20) 120.6
𝐾𝑒 = [ ][ ]=[ ]
0 1 20 − 0 20
𝑘𝑒1 120.6
Hence, 𝐾𝑒 = [ ]=[ ]
𝑘𝑒2 20
Step 1: Check the observability of the given system. If it is observable, then go to step2
𝐶
𝐶𝐴
| 𝐶𝐴2 |
Observability matrix |𝑶𝑴| = ≠𝟎
| ; |
𝐶𝐴𝑛−2
𝐶𝐴𝑛−1
𝐶
Observability matrix |𝑂𝑀| = | | ≠ 0
𝐶𝐴
[Link]’s College of Engineering 42 VI Semester ICE
IC6601 Advanced Control System Department of ICE 2018-19
𝐶 = [0 1];
0 20.6
𝐶 × 𝐴 = [0 1] × [ ] = [1 0]
1 0
|𝑂𝑀| = | 𝐶 | = [0 1] = −1 ≠ 0
𝐶𝐴 1 0
Hence the system is completely state observable.
𝒌𝒆𝟏
Step 2: 𝒏 = 𝟐, then 𝑲𝒆 = [ ]
𝒌𝒆𝟐
Step 3: Substitute this 𝑲𝒆 matrix into the desired characteristics polynomial |𝒔𝑰 − (𝑨 − 𝑲𝒆 𝑪)|
and equate it to (𝒔 − 𝝁𝟏 )(𝒔 − 𝝁𝟐 )(𝒔 − 𝝁𝟑 ).
𝑘 𝑠 −20.6 + 𝑘𝑒1
|𝑠𝐼 − (𝐴 − 𝐾𝑒 𝐶)| = |[ 𝑠 0] − [0 20.6] + [ 𝑒1 ] [0 1]| = [ ]
0 𝑠 1 0 𝑘𝑒2 −1 𝑠 + 𝑘𝑒2
Step 4: By equating the coefficients of like power of ′𝑺′ on both sides, it is possible to determine
the values of 𝑲𝒆𝟏 , 𝑲𝒆𝟐
𝑘𝑒2 = 20;
𝑘𝑒1 = 120.6
𝑘𝑒1 120.6
∴ 𝑘𝑒 = [ ]=[ ]
𝑘𝑒2 20
Step 1: Check the observability of the given system. If it is observable, then go to step2
𝐶
𝐶𝐴
| 𝐶𝐴2 |
Observability matrix |𝑶𝑴| = ≠𝟎
| ; |
𝐶𝐴𝑛−2
𝐶𝐴𝑛−1
𝐶
Observability matrix |𝑂𝑀| = | | ≠ 0
𝐶𝐴
𝐶 = [0 1];
[Link]’s College of Engineering 43 VI Semester ICE
IC6601 Advanced Control System Department of ICE 2018-19
0 20.6
𝐶 × 𝐴 = [0 1] × [ ] = [1 0]
1 0
|𝑂𝑀| = | 𝐶 | = [0 1] = −1 ≠ 0
𝐶𝐴 1 0
Hence the system is completely state observable.
Step 2: Determine the desired characteristic equation from the specified closed loop poles,
𝝁𝟏 , 𝝁𝟐 , 𝝁𝟑 , . . , 𝝁𝒏
i.e., (𝒔 − 𝝁𝟏 ) (𝒔 − 𝝁𝟐 ) (𝒔 − 𝝁𝟑 ). . . . . . . . (𝒔 − 𝝁𝒏 ) = 𝒔𝒏 + 𝜶𝟏 𝒔𝒏−𝟏 +, … … . . , +𝜶𝒏−𝟏 𝒔𝟏 +
𝜶𝒏
and determine the values of 𝜶𝟏 , 𝜶𝟐 , … … . , 𝜶𝒏 .
(𝑠 − 𝜇1 ) (𝑠 − 𝜇2 ) (𝑠 − 𝜇3 ) = (𝑠 + 10)(𝑠 + 10) = 𝑠 2 + 20𝑠 + 100
Step 3: Determine the matrix ∅(𝑨) using the coefficient of desired characteristics polynomial.
∅(𝑨) = 𝑨𝒏 + 𝜶𝟏 𝑨𝒏−𝟏 + 𝜶𝟐 𝑨𝒏−𝟐 +, … … . . , +𝜶𝒏−𝟏 𝑨𝟏 + 𝜶𝒏 𝑰
0 20.6
𝐴=[ ]
1 0
0 20.6 0 20.6 20.6 0
𝐴2 = 𝐴 × 𝐴 = [ ]×[ ]=[ ]
1 0 1 0 0 20.6
20.6 0 0 20.6 0 424.36
𝐴3 = 𝐴2 × 𝐴 = [ ]×[ ]=[ ]
0 20.6 1 0 20.6 0
∅(𝑨) = 𝑨𝟐 + 𝜶𝟏 𝑨 + 𝜶𝟐 𝑰
0 20.6 2 20.6 0 1 0
=[ ] + 20 [ ] + 100 [ ]
1 0 0 20.6 0 1
0 424.36 20.6 0 1 0
=[ ] + 20 [ ] + 100 [ ]
20.6 0 0 20.6 0 1
120.6 412
=[ ]
20 120.6
Step 4: Calculate the state feedback gain matrix, ′𝑲𝒆 ′ using the Ackermann’s formula
𝐶 −𝟏 𝟎
𝐶𝐴 𝟎
𝑲𝒆 = [∅(𝑨)] 𝐶𝐴2 ;
; ;
[𝐶𝐴𝑛 ] [𝟏]
𝟎
𝑲𝒆 = [∅(𝑨)] [𝑶𝑴]−𝟏 [ ]
𝟏
𝑘𝑒1 120.6 412 0 1 0 120.6
𝑲𝒆 = [ ]=[ ]×[ ]×[ ] =[ ]
𝑘𝑒2 20 120.6 1 0 1 20
𝟏 𝟐 𝟎 𝟐
𝒙̇ = [𝟑 −𝟏 𝟏] 𝒙 + [𝟐] 𝒖; 𝒚 = [𝟎 𝟎 𝟏] 𝒙
𝟎 𝟐 𝟎 𝟏
Design a state observer so that eigen values of state matrix (𝑨 − 𝑮𝑪) are −𝟒, −𝟑 ± 𝑱𝟏
Solution
1 2 0 2
Given 𝐴 = [3 −1 1] ; 𝐵 = [2] ; 𝐶 = [0 0 1]; 𝜇1 = −4, 𝜇2 = −3 − 𝑗1, 𝜇3 = −3 −
0 2 0 1
𝑗1; 𝑛 = 3
Step 1: Check the Observability of the given system. If it is observable, then go to step 2
𝐶
Observability matrix |𝑂𝑀| = | 𝐶𝐴 | ≠ 0
𝐶𝐴2
𝐶 = [0 0 1]
1 2 0
𝐶𝐴 = 𝐶 × 𝐴 = [0 0 1] × [3 −1 1] = [0 2 0]
0 2 0
1 2 0
𝐶𝐴2 = 𝐶𝐴 × 𝐴 = [0 2 0] × [3 −1 1] = [6 −2 2]
0 2 0
𝐶 0 0 1
|𝑂𝑀| = | 𝐶𝐴 | = |0 2 0| = 0 | 2 0| − 0 |0 0 0 2
| + 1| | = −12 ≠ 0
−2 2 6 2 6 −2
𝐶𝐴2 6 −2 2
Hence the system is completely state Observable.
Step 2: Determine the desired characteristic equation from the specified closed loop poles,
𝝁𝟏 , 𝝁𝟐 , 𝝁𝟑 , . . , 𝝁𝒏
i.e., (𝒔 − 𝝁𝟏 ) (𝒔 − 𝝁𝟐 ) (𝒔 − 𝝁𝟑 ). . . . . . . . (𝒔 − 𝝁𝒏 ) = 𝒔𝒏 + 𝜶𝟏 𝒔𝒏−𝟏 +, … … . . , +𝜶𝒏−𝟏 𝒔𝟏 + 𝜶𝒏
and determine the values of 𝜶𝟏 , 𝜶𝟐 , … … . , 𝜶𝒏 .
Step 3: Determine the matrix ∅(𝑨) using the coefficient of desired characteristics polynomial.
1 2 0
𝐴 = [3 −1 1]
0 2 0
1 2 0 1 2 0 7 0 2
𝐴2 = 𝐴 × 𝐴 = [3 −1 1] × [3 −1 1] = [0 9 −1]
0 2 0 0 2 0 6 −2 2
7 0 2 1 2 0 7 18 0
𝐴3 = 𝐴2 × 𝐴 = [0 9 −1] × [3 −1 1] = [27 −11 9]
6 −2 2 0 2 0 0 18 −2
∅(𝑨) = 𝑨𝟑 + 𝜶𝟏 𝑨𝟐 + 𝜶𝟐 𝑨 + 𝜶𝟑 𝑰
1 2 0 3 1 2 02 1 2 0 1 0 0
=[3 −1 1] + 10 [3 −1 1] + 34 [3 −1 1] + 40 [0 1 0]
0 2 0 0 2 0 0 2 0 0 0 1
7 18 0 7 0 2 1 2 0 1 0 0
= [27 −11 9 ] + 10 [0 9 −1] + 34 [3 −1 1] + 40 [0 1 0]
0 18 −2 6 −2 2 0 2 0 0 0 1
151 86 20
= [129 85 33]
60 66 58
Step 4: Calculate the state feedback gain matrix, ′𝑲′ using the Ackermann’s formula
𝐶 −𝟏 𝟎
𝐶𝐴 𝟎
𝑲= [∅(𝑨)] 𝐶𝐴2 ;
; ;
[𝐶𝐴𝑛 ] [𝟏]
𝟎
[∅(𝑨)] [𝑶𝑴] −𝟏
𝑲= [𝟎 ]
𝟏
0 0 1 −1 𝐴𝑑𝑗(𝑂𝑀)
[𝑂𝑀]−1 = [0 2 0] =
|𝑂𝑀|
6 −2 2
𝐶 0 0 1
|𝑂𝑀| = | 𝐶𝐴 | = |0 2 0| = 0 | 2 0| − 0 |0 0| + 1 |0 2 | = −12
−2 2 6 2 6 −2
𝐶𝐴2 6 −2 2
2 0 0 0 0 2 𝑇
+| | −| | +| |
−2 2 6 2 6 −2 4 0 −12 𝑇 4 2 −2
0 1 0 1 0 0
𝐴𝑑𝑗(𝑂𝑀) = − | | +| | −| | = [ 2 −6 0 ] =[ 0 −6 0 ]
−2 2 6 2 6 −2
0 1 0 1 0 0 −2 0 0 −12 0 0
[ + |2 0
| −|
0 0
| +|
0 2
|]
𝟎 𝟏 𝟎 𝟎
̇𝑿 = [ 𝟎 𝟎 𝟏 ] 𝑿 + [𝟎] 𝑼; 𝒀 = [𝟏 𝟎 𝟎] 𝑿
−𝟔 −𝟏𝟏 −𝟔 𝟏
Design a full order observer with eigenvalues −𝟓, −𝟐 ± 𝟑. 𝟒𝟔𝟒𝒋
Solution
𝟎 𝟏 𝟎 𝟎
Given 𝐴 = [ 𝟎 𝟎 𝟏 ] ; 𝐵 = [ 𝟎 ] ; 𝐶 = [𝟏 𝟎 𝟎]; 𝜇1 = −5, 𝜇2 = −2 − 3.464𝑗, 𝜇3 =
−𝟔 −𝟏𝟏 −𝟔 𝟏
−2 − 3.464𝑗; 𝑛 = 3
Step 1: Check the Observability of the given system. If it is observable, then go to step 2
𝐶
Observability matrix |𝑂𝑀| = | 𝐶𝐴 | ≠ 0
𝐶𝐴2
𝐶 = [𝟏 𝟎 𝟎]
𝟎 𝟏 𝟎
𝐶𝐴 = 𝐶 × 𝐴 = [𝟏 𝟎 𝟎] × [ 𝟎 𝟎 𝟏 ] = [0 1 0]
−𝟔 −𝟏𝟏 −𝟔
𝟎 𝟏 𝟎
2
𝐶𝐴 = 𝐶𝐴 × 𝐴 = [0 1 0] × [ 𝟎 𝟎 𝟏 ] = [0 0 1]
−𝟔 −𝟏𝟏 −𝟔
𝐶 1 0 0
1 0
|𝑂𝑀| = | 𝐶𝐴 | = |0 1 0| = 1 | |=1≠0
0 1
𝐶𝐴2 0 0 1
Hence the system is completely state Observable.
Step 2: Determine the desired characteristic equation from the specified closed loop poles,
𝝁𝟏 , 𝝁𝟐 , 𝝁 𝟑 , . . , 𝝁𝒏
i.e., (𝒔 − 𝝁𝟏 ) (𝒔 − 𝝁𝟐 ) (𝒔 − 𝝁𝟑 ). . . . . . . . (𝒔 − 𝝁𝒏 ) = 𝒔𝒏 + 𝜶𝟏 𝒔𝒏−𝟏 +, … … . . , +𝜶𝒏−𝟏 𝒔𝟏 + 𝜶𝒏
and determine the values of 𝜶𝟏 , 𝜶𝟐 , … … . , 𝜶𝒏 .
Hence, 𝛼1 = 9; 𝛼2 = 36; 𝛼3 = 80
[Link]’s College of Engineering 47 VI Semester ICE
IC6601 Advanced Control System Department of ICE 2018-19
Step 3: Determine the matrix ∅(𝑨) using the coefficient of desired characteristics polynomial.
∅(𝑨) = 𝑨𝒏 + 𝜶𝟏 𝑨𝒏−𝟏 + 𝜶𝟐 𝑨𝒏−𝟐 +, … … . . , +𝜶𝒏−𝟏 𝑨𝟏 + 𝜶𝒏 𝑰
𝟎 𝟏 𝟎
𝐴=[ 𝟎 𝟎 𝟏]
−𝟔 −𝟏𝟏 −𝟔
𝟎 𝟏 𝟎 𝟎 𝟏 𝟎 0 0 1
2
𝐴 =𝐴×𝐴=[ 𝟎 𝟎 𝟏 ]×[ 𝟎 𝟎 𝟏 ] = [−6 −11 −6]
−𝟔 −𝟏𝟏 −𝟔 −𝟔 −𝟏𝟏 −𝟔 36 60 25
0 0 1 𝟎 𝟏 𝟎 −6 −11 −6
𝐴3 = 𝐴2 × 𝐴 = [−6 −11 −6] × [ 𝟎 𝟎 𝟏 ] = [ 36 60 25 ]
36 60 25 −𝟔 −𝟏𝟏 −𝟔 −150 −239 −90
∅(𝑨) = 𝑨𝟑 + 𝜶𝟏 𝑨𝟐 + 𝜶𝟐 𝑨 + 𝜶𝟑 𝑰
𝟎 𝟏 𝟎 3 𝟎 𝟏 𝟎 2 𝟎 𝟏 𝟎 1 0 0
=[ 𝟎 𝟎 𝟏 ] +9[ 𝟎 𝟎 𝟏 ] + 36 [ 𝟎 𝟎 𝟏 ] + 80 [0 1 0]
−𝟔 −𝟏𝟏 −𝟔 −𝟔 −𝟏𝟏 −𝟔 −𝟔 −𝟏𝟏 −𝟔 0 0 1
−6 −11 −6 0 0 1 𝟎 𝟏 𝟎 1 0 0
= [ 36 60 25 ] + 9 [−6 −11 −6] + 36 [ 𝟎 𝟎 𝟏 ] + 80 [ 0 1 0]
−150 −239 −90 36 60 25 −𝟔 −𝟏𝟏 −𝟔 0 0 1
74 25 3
= [−18 41 7]
−42 −95 −1
Step 4: Calculate the state feedback gain matrix, ′𝑲′ using the Ackermann’s formula
𝐶 −𝟏 𝟎
𝐶𝐴 𝟎
𝑲 = [∅(𝑨)] 𝐶𝐴2 ;
; ;
[𝐶𝐴 ]𝑛 [ 𝟏]
𝟎
𝑲 = [∅(𝑨)] [𝑶𝑴]−𝟏 [𝟎]
𝟏
1 0 0 −1 𝐴𝑑𝑗(𝑂𝑀)
[𝑂𝑀]−1 = [0 1 0] =
|𝑂𝑀|
0 0 1
𝐶 1 0 0
1 0
|𝑂𝑀| = | 𝐶𝐴 | = |0 1 0| = 1 | |=1
2 0 1
𝐶𝐴 0 0 1
1 0 0
𝐴𝑑𝑗(𝑂𝑀) = [0 1 0]
0 0 1
𝐴𝑑𝑗(𝑂𝑀) 1 1 0 0 1 0 0
𝑂𝑀−1 = = [0 1 0 ] = [0 1 0 ]
|𝑂𝑀| 1
0 0 1 0 0 1
74 25 3 1 0 0 0
𝐾𝑒 = [−18 41 7 ] × [0 1 0 ] × [0 ]
−42 −95 −1 0 0 1 1
74 25 3 0 3
𝐾𝑒 = [−18 41 7 ] × [0] = [ 7 ]
−42 −95 −1 1 −1
𝑘𝑒1 3
𝐾𝑒 = [𝑘𝑒2 ] = [ 7 ]
𝑘𝑒3 −1
Nonlinear Systems
Nonlinear systems are systems which does not obey the principle of superposition. Principle of
superposition implies that if a system has responses 𝑦1 (𝑡) and 𝑦2 (𝑡) to any two inputs 𝑥1 (𝑡) and
𝑥2 (𝑡) respectively then the system response to the linear combination of these inputs 𝛼1 (𝑡)𝑥1 (𝑡) +
𝛼2 (𝑡)𝑥2 (𝑡) is given by the linear combination of the individual outputs. i.e., 𝛼1 (𝑡)𝑦1 (𝑡) +
𝛼2 (𝑡)𝑦2 (𝑡) where 𝛼1 and 𝛼2 constants.
1. The response of nonlinear system to a particular test signal is no guide to their behaviour to
their inputs, since the principle of superposition does not hold good for nonlinear systems.
2. Nonlinear system response may be highly sensitive to input amplitude. The stability study
of nonlinear systems requires the information about the type and amplitude of the
anticipated inputs, initial conditions, etc., in addition to the usual requirement of the
mathematical model.
5. The output of a nonlinear system will have harmonics and sub-harmonics when excited by
sinusoidal signals.
6. Nonlinear systems will exhibit phenomena like frequency entrainment and asynchronous
quenching.
The differential equation governing the dynamics of the system may be written as
𝑀𝑥̈ + 𝐵𝑥̇ + 𝐾𝑥 + 𝐾 ′ 𝑥 3 = 0
The parameters M, B and K are positive constants. The parameter K' may be positive or
negative. If K' is positive, the spring is called hard spring and if K' is negative the spring is called
soft spring. The equation is nonlinear differential equation and it is called Duffing’s equation.
When the system has non zero initial conditions, the free response is damped oscillatory. The
frequency of free oscillations depends on the amplitude of oscillations. When K' < 0 (soft spring)
the frequency decreases with decreasing amplitude. When K' > 0 (hard spring) the frequency
increases with decreasing amplitude. When K' = 0 (corresponding to linear system) the frequency
remains unchanged as the amplitude of free oscillation decreases. The frequency-amplitude
dependence characteristic of nonlinear mechanical system is shown in figure.
Jump resonance
In the frequency response of nonlinear systems, the amplitude of the response (output) may
jump from one point to another for increasing or decreasing values of frequency, . This
phenomenon is called jump resonance and it can be observed in the frequency response of the
system shown in figure when it is subjected to sinusoidal input. Let the mechanical system shown
in figure 1 be subjected to an input of type 𝐴 𝐶𝑜𝑠 𝜔𝑡. Now the differential equation governing
the mechanical system is
Let X be the amplitude of the response or output of the system. In frequency response studies,
the amplitude, A of the input is held constant, while its is varied and the amplitude X of the
output is observed. The frequency response curve is plotted between X and . The frequency
response curves of the mechanical system of figure 1are shown in figure 2.
In the frequency response curve shown in figure, as the frequency is increased, the
amplitude X increases, until point-2 is reached. A further increase in frequency will cause a jump
from point-2 to point-3. This phenomenon is called jump resonance. As the frequency is
increased further, the amplitude X follows the curve from point-3 towards point-4. When the
frequency is reduced starting from a high value corresponding to point-4, the amplitude X slowly
increases through point-3, until point-5 is reached. A further decrease in will cause another
jump from point-5 to point-6. This phenomenon is called jump resonance. After this jump, the
amplitude X decreases with and follows the curve from point-6 towards point-1. For jump
resonance to take place, it is necessary that the damping term be small and the amplitude of the
forcing function be large enough to drive the system into a region of appreciably nonlinear
operation.
Subharmonic oscillations
When a nonlinear system is excited by a sinusoidal signal, the response or output will have
steadystate oscillation whose frequency is an integral submultiple of the forcing frequency.
These oscillations are called subharmonic oscillations. The generation of subharmonic
oscillations depends on the system parameters and initial conditions. It also depends on
amplitude and frequency of the forcing functions.
Limit Cycles
The response of nonlinear system may exhibit oscillations with fixed amplitude and
frequency. These oscillations are called limit cycles. Consider a mechanical system with a
nonlinear damping and described by the equation.
where M, B and K are positive constants. The equation is called the van der pol equation. For small
values of x the damping will be negative which implies the stored energy in the damper is fed to
the system. For large values of x the damping is positive which implies that it absorbs energy from
the system. Thus, it can be expected that such, a system may exhibit a sustained oscillation. Since
the system explained above is not a forced system, this oscillation is called a self-excited oscillation
or zero input limit cycle.
Asynchronous quenching
In a nonlinear system that exhibits a limit cycle of frequency l it is possible to quench (stop
or eliminate) the limit cycle oscillation by forcing the system of a frequency q where q and ,
are not related to each other. The phenomenon is called signal stabilization or asynchronous
quenching
The intentional nonlinearities are those which are deliberately inserted in the system to
modify system characteristics. The most common example of this type of nonlinearity is a relay.
Saturation: In this type of nonlinearity the output is proportional to input for a limited range
of input signals. When the input exceeds this range, the output tends to become nearly constant
as shown in figure.
All devices when driven by sufficiently large signals, exhibit the phenomenon of saturation due
to limitations of their physical capabilities. Saturation in the output of electronic, rotating and
flow (hydraulic and pneumatic) amplifiers, speed and torque saturation in electric and hydraulic
motors, saturation in the output of sensors for measuring position, velocity, temperature, etc., are
the well known examples.
Dead Zone : The dead zone is the region in which the output is zero for a given input. Many
physical devices do not respond to small signals, i.e., if the input amplitude is less than some small
value, there will be no output. The region in which the output is zero is called dead zone. When the
input is increased beyond this dead zone value, the output will be linear. Figure shows the Dead zone
nonlinearity and combination of dead zone a nd saturation nonlinearity.
Friction: Friction exists in any system when there is relative motion between contacting surfaces.
The different types of friction are viscous friction, coulomb friction and stiction. The viscous
friction is linear in nature and the frictional force is directly proportional to relative velocity of the
sliding surfaces. The coulomb friction and stiction are nonlinear frictions. The coulomb friction
offers a constant retarding force only when the motion is initiated. Due to interlocking of surface
irregularities, more force is required to move an object from rest than to maintain it in motion.
Hence the force of stiction is always greater than that of coulomb friction.
In actual practice, the stiction force gradually decreases with velocity and changes over to coulomb
friction at reasonably low velocities as shown in figure. The composite characteristics of various
frictions are shown in figure.
METHOD OF LINEARISATION
Local linearization has been extensively used while analyzing nonlinear phenomena in many
engineering situations. Sometimes, this may also be called perturbation analysis. Consider a time
invariant nonlinear system
𝑥̇ = 𝑓(𝑥, 𝑢) (1)
And ℎ(𝑥𝑒 , 𝑢0 , 𝛿𝑥, 𝛿𝑢) contains the remainder higher order terms. Assume that as 𝛿𝑥 and 𝛿𝑢 go to
zero, the term ℎ(𝑥𝑒 , 𝑢0 , 𝛿𝑥, 𝛿𝑢) tends to zero at a faster rate that the first order terms on the right
hand side of equation (4). Specifically, we assume that;
Then, in a sufficiently small region about 𝑥𝑒 and for sufficiently small values of |𝛿𝑢|, the behaviour
of the perturbed system can be approximated by the locally linearized equation;
𝜕𝑓 𝜕𝑓
𝛿𝑥̇ (𝑡) = 𝜕𝑥 (𝑥𝑒 , 𝑢0 )𝛿𝑥 + 𝜕𝑢 (𝑥𝑒 , 𝑢0 )𝛿𝑢 (6)
This being a set of linear and time invariant system of equations makes it much simpler to study
from an analytical point of view. All the methods available for linear systems like analysis based
on eigen value or roots of the characteristic equation are applicable to determine the stability of this
linearized system.
The phase plane method of analysis is a graphical method for the analysis of linear and
nonlinear systems. The analysis is carried by constructing phase trajectories. It gives an idea about
the transient behaviour and stability of the system. The phase plane analysis is usually restricted
to second order systems excited by step or ramp inputs. This analysis technique can be extended
to a higher order system if it is approximated as a second order system. The dynamics of control
systems can be represented by differential equations. A second order linear system can be
represented by the differential equation
𝑑2 𝑥 𝑑𝑥
+ 2𝜔𝑛 𝑑𝑡 + 𝜔𝑛2 𝑥 = 0 (1)
𝑑𝑡 2
[Link]’s College of Engineering 54 VI Semester ICE
IC6601 Advanced Control System Department of ICE 2018-19
where,
𝑥 = One of the system variable (e.g, displacement in mechanical system, current in
electrical system, etc.,)
= Damping ratio
n = Natural frequency of oscillation.
The state of the second order system represented by equ(1) can be described by choosing
two state variables. In state space modelling using phase variables we choose one of the system
variable and its derivatives as state variables. Let 𝑥1 and 𝑥2 be the state variables of the second
order system.
𝑑𝑥
Here 𝑥1 = 𝑥 and 𝑥2 = (2)
𝑑𝑡
On substituting the state variables in equation (1)
𝑥̇ 2 + 2𝜔𝑛 𝑥2 + 𝜔𝑛2 𝑥1 = 0 (3)
The state equations of the system are obtained from equation (2) and (3). The state equations
are,
𝑥̇ 1 = 𝑥2 (4)
𝑥̇ 2 = −𝜔𝑛 𝑥1 − 2𝜔𝑛 𝑥2
2
(5)
For linear systems the state equations are a set of first order linear differential equations
and solutions of state equations can be easily obtained by integration. But for nonlinear systems,
the state equations are a set of first-order nonlinear differential equations and solving the nonlinear
differential equations will not be an easy task. Hence, for nonlinear systems the phase plane
method of analysis will be a useful tool.
Phase plane: The coordinate plane with state variables 𝑥1 and 𝑥2 as two axis is called the phase
plane.
Phase Trajectory: The curve describing the state point (𝑥1 , 𝑥2 ) in the phase plane with time as
running parameter is called phase trajectory.
Phase Portrait: A trajectory can be constructed in the phase plane for each set of initial conditions.
Hence a family of trajectories can be constructed for a system in a phase plane and such a family
of trajectories is called the ‘Phase Portrait’.
Singular point: A point in phase plane at which the derivatives of all state variables are zero is
called a singular point. It is also called equilibrium point. At the singular point, the slope of the
trajectory is indeterminant. If the system is placed at a singular or equilibrium point, it will
continue to lie there if left undisturbed (i.e., the derivatives of all the phase variables are zero and
so the system state remains unchanged).
The singular points are classified as Nodal point, Saddle point, Focus point and Centre
or Vortex point depending on the eigen values of the system matrix.
Table 1. Phase portrait of systems with various types of singular points.
[Link] Eigen values of system Type of singular Phase portrait of the system
matrix point with singular point at origin
Complex-conjugate
4 with negative real part Stable focus
Complex conjugate
5 with positive real part Unstable focus
For linear time invariant systems, the concept of stability can be defined as follows,
1. When the input is zero, the system is stable for arbitrary initial conditions if the resulting
trajectory tends towards the equilibrium state.
2. When the system is excited by a bounded input, the system is stable if the system output is
bounded.
In nonlinear systems the concept of stability is not clear-cut. There are many types of
stability definitions in the literature. The linear autonomous system has only one equilibrium state.
[Link]’s College of Engineering 56 VI Semester ICE
IC6601 Advanced Control System Department of ICE 2018-19
The behaviour of linear system about the equilibrium state completely determines the qualitative
behaviour in the entire state-plane.
In nonlinear systems there may be multiple equilibrium state. The behaviour of nonlinear
system about the equilibrium point may be different for small deviations and large deviations about
the equilibrium point. In nonlinear systems with multiple equilibrium states, the system trajectories
may move away from one equilibrium point and tend to other as time progresses. Hence in
nonlinear systems, stability is discussed relative to the equilibrium state and the general stability of
system cannot be defined.
Consider an autonomous system described by the state equation, x = f(x). Let us assume
that the system has one equilibrium point and the origin of phase plane is the equilibrium point. For
this system, the following definitions of stability are proposed.
1. The autonomous system defined by equation 𝑥 = 𝑓(𝑥) is stable at the origin, if for every
initial state 𝑥(𝑡0 ) which is sufficiently close to origin, 𝑥(t) remains near the origin for all t.
2. The autonomous system defined by equation 𝑥 = 𝑓(𝑥) is asymptotically stable if 𝑥(t)
approaches the origin as 𝑡 → .
3. The autonomous system defined by equation 𝑥 = 𝑓(𝑥) is asymptotically stable in the large
if it is asymptotically stable for every initial state regardless of how near or for it is from the
origin.
The limit cycles are oscillations of fixed amplitude and period. The existence of limit cycle
in nonlinear systems can be predicted from closed trajectories in the phase portrait. In linear
systems, when oscillations occur, the resulting trajectories will be closed curves as shown in figure
2. The amplitude of the oscillations is not fixed. It changes with the size of the initial conditions.
Slight changes in system parameters will destroy the oscillations
The limit cycles are oscillations of the response of nonlinear systems with fixed amplitude and
frequency. If these oscillations or limit cycles exists when there is no input then they are called
zero input limit cycles.
In nonlinear systems, there can be limit cycles (oscillations) that are independent of the size
of initial conditions as shown in figure 3. These limit cycles are usually less sensitive to system
parameter variations. Limit cycles of fixed amplitude and period can be sustained over a finite range
of system parameters. The limit cycle is stable if the paths in its neighbourhood converge towards
the limit cycle as shown in figure 3. The limit cycle is unstable if the paths in the neighbourhood
of a limit cycle diverge away from it as shown in figure 4.
Figure 1. Phase portrait showing limit Figure 2. Unstable limit Figure 3. Stable limit cycle
cycle behavior in linear system cycle
The equation (4) defines the slope of phase trajectory at every point in the phase- plane,
except at singular points. At singular points, the slope of the phase trajectory is indeterminate. The
phase trajectory can be constructed using the slope equation (4) either analytically or graphically.
In analytical method of construction, the equation (4) is integrated and the resulting equation is
used to construct phase trajectories for a given set of initial conditions. In many cases, it is not
possible to perform integration of the slope equation (4). Hence a number-of graphical methods
have been developed for construction of phase trajectories and the two popular methods are isocline
method and delta method.
The analytical method is used if the differential equations describing the system can be
approximated by piecewise linear differential equations. (i.e., the equations are linearized for small
regions). The slope equation 𝑑𝑥2 ⁄𝑑𝑥1 is formed from the state equations. Then the slope equation
is splitted into sections of linear equations. Each section of linear slope equation is directly
integrated to get the solution of state equations. Each section of solution equation is used to
construct a section of phase trajectories for various sets of initial conditions. The solution of slope
equation will be a function of 𝑥1 and 𝑥2 . For a given set of initial conditions, assume different
values of 𝑥1 , and calculate 𝑥2 for each value of 𝑥1 , using the solution of slope equation. The
values of 𝑥1 and 𝑥2 are tabulated. In an ordinary graph sheet, take 𝑥1 and 𝑥2 is, choose appropriate
scales and mark the phase points (𝑥1 and 𝑥2 ). Join all the points by a smooth curve. For each set
of initial conditions one phase trajectory can be constructed using the procedure described above:
Eventhough, time domain solutions are obtained by direct integration, the construction of
phase-trajectories (or phase plane analysis) will help in investigating system behaviour and the
design of system parameters to achieve a desired response. Also the existence of limit cycles is
sharply brought into focus by the phase-portrait.
[Link]’s College of Engineering 58 VI Semester ICE
IC6601 Advanced Control System Department of ICE 2018-19
Example 1
Consider a system with an ideal relay as shown in figure. Determine the singular point.
Construct phase trajectories, corresponding to initial conditions, (𝒊) 𝒄(𝟎) = 𝟐, 𝒄̇ (𝟎) =
𝟏 𝒂𝒏𝒅 (𝒊𝒊) 𝒄(𝟎) = 𝟐, 𝒄̇ (𝟎) = 𝟏. 𝟓. Take 𝒓 = 𝟐 volts amd M = 𝟏. 𝟐 volts.
SOLUTION
𝑑𝑥 𝑓 (𝑥 ,𝑥 )
Slope 𝑆 = 𝑑𝑥2 = 𝑓2 (𝑥1 ,𝑥2) (1)
1 1 1 2
Equation (2) defines the locus of all such points in phase plane at which the slope of the
phase trajectory is S1. A locus passing through the points of same slope in phase plane is called
isocline. The slope of a phase trajectory at the crossing point of an isocline will be the slope of
corresponding isocline. A typical plot of isoclines for various values of slope, S is shown in figure.
Using these isoclines the phase trajectories can be constructed as explained below.
The phase trajectory starts at a point corresponding to initial conditions. (For each set of
initial conditions one phase trajectory can be constructed).
Let, S1, S2, S3 etc, be the slopes associated with isoclines 1,2,3, etc.,
In figure, let point A on isocline-1be the point corresponding to a set of initial conditions.
The phase-trajectory will leave the point A at a slope S1. When the trajectory reaches the isocline-
2, the slope changes to S2. Draw two lines from point A one at a slope of S 1 (i.e., at angle of 𝛼1 =
𝑡𝑎𝑛−1 (𝑆1 ) and the other at a slope of S2, (i.e., at angle of 𝛼1 = 𝑡𝑎𝑛−1 (𝑆1 )). Let these two lines
meet the isocline- 2 at p and q. Now we can say that the trajectory would cross the iscoline-2 at a
point midway between p and q. Mark the point B on the isocline- 2 approximately midway between
p and q. The constructional procedure is now repeated at B to find the crossing point C on the
isocline-3. By similar procedures the crossing points on the isolines are determined.
A smooth curve drawn through the crossing points gives the phase-trajectory starting at
point A. The accuracy of the trajectory is closely related to the spacing of the isoclines. The phase
trajectory will be more accurate if large numbers of isoclines are used, which are very close to each
other. It should be noted that using a set of isoclines, any number of trajectories can be constructed.
𝒆̈ + 𝟐𝝎𝒏 𝒆̇ + 𝝎𝟐𝒏 𝒆 = 𝟎
𝒓𝒂𝒅
Where = 𝟎. 𝟏𝟓, 𝝎𝒏 = 𝟏 𝒔𝒆𝒄 , 𝒆(𝟎) = 𝟏. 𝟓 and 𝒆̇ (𝟎) = 𝟎. Determine the singular point.
Construct the phase trajectory using the method of isoclines.
Solution
Let 𝑥1 and 𝑥2 be the state variables of the system and they are related to the system variables, 𝑒
𝑥1 = 𝑒 (1)
𝑥2 = 𝑒̇ (2)
On differentiating (1)
𝑥̇ 1 = 𝑒̇ = 𝑥2 (3)
On differentiating (2)
𝑥̇ 2 = 𝑒̈ (4)
Given that
𝑒̈ + 2 𝜔𝑛 𝑒̇ + 𝜔2𝑛 𝑒 = 0 (5)
[Link]’s College of Engineering 63 VI Semester ICE
IC6601 Advanced Control System Department of ICE 2018-19
Substitute equation (1), (2) and (4) in (5)
𝑥̇ 2 + 2 𝜔𝑛 𝑥2 + 𝜔2𝑛 𝑥1 = 0
𝑥̇ 2 = −2 𝜔𝑛 𝑥2 − 𝜔2𝑛 𝑥1 (6)
The state equation of the system are given by equation (3) and (6)
𝑥̇ 1 = 𝑥2
𝑥̇ 2 = −2 𝜔𝑛 𝑥2 − 𝜔2𝑛 𝑥1
0 = −2 𝜔𝑛 𝑥2 − 𝜔𝑛2 𝑥1
Therefore the coordinates of singular point are (0,0) and so the origin is the singular point.
𝑑𝑥 ⁄𝑑𝑡 𝑥̇
𝑆 = 𝑑𝑥2⁄𝑑𝑡 = 𝑥̇ 2 (7)
1 1
On substituting for 𝑥1̇ and 𝑥̇ 2 from equations (3) and (6) in equation (7), we get
From equation (9) we can conclude that the isoclines are straight lines. For each value of S
we can draw one isocline. Using equation (9), the coordinates (𝑥1 , 𝑥2 ) in the phase plane for various
slopes can be calculated. Since there are three variables. Let us assume two variables and calculate
the third variable.
For each value of S, choose two values of 𝑥1 and calculate 𝑥2 using equation (9_. The value
of 𝑆, 𝑥1 and 𝑥2 are tabulated. The slope angle, 𝛼 is calculated for each value of S, using the
expression, 𝛼 = tan−1(𝑆), and tabulated in table.
S -2.0 -1.0 -0.5 0 0.5 1.0 2.0
𝛼 -63 -45 -27 0 27 45 63
𝑥1 𝑥2 𝑥1 𝑥2 𝑥1 𝑥2 𝑥1 𝑥2 𝑥1 𝑥2 𝑥1 𝑥2 𝑥1 𝑥2
1.0 0.6 1.0 1.4 0.25 1.25 0.25 -0.8 1.0 -1.25 1.0 -0.77 1.0 -0.43
2.0 1.2 1.5 2.1 0.5 2.5 0.75 -2.5 2.0 -2.5 2.0 -1.54 2.0 -0.86
Result
1. The singular point lies at origin
2. From figure it is observed that the phase trajectory spiral towards the origin, hence the type
of singular point is stable focus.
[Link]’s College of Engineering 64 VI Semester ICE
IC6601 Advanced Control System Department of ICE 2018-19
Unit III
DESCRIBING FUNCTION FOR COMMON NON-LINEARITIES
SATURATION
It is the most commonly encountered nonlinearity in control system. It is often associated with
amplifiers and actuators.
In transistor amplifiers, the output varies linearly with the input,
only for a small amplitude limits.
When the input amplitude gets out of the linear range of the
amplifier, the output changes very little and stays close to its
maximum value.
Most of the actuators display saturation characteristics.
For eg., the output torque of a servo motor cannot increase
infinitely and tends to saturate due to the properties of magnetic
material. Similarly, valve controlled hydraulic actuators are saturated by the maximum flow
rate.
DEAD ZONE
The dead-zone is the region in which the output is zero for a
given input. When the input is increased beyond this dead-
zone value, the output will be linear. The region where the
applied voltage has no effect is called the deadband e.g.
sensors, friction in all components
Dead zone nonlinearity may occur in sensors, amplifiers and
actuators.
In a de motor, we assume that any voltage applied to armature winding will cause armature to
rotate if the field current is maintained constant.
In reality, due to static friction at the motor shaft, rotation will occur only if the torque provided
by the motor is sufficiently large. This corresponds to dead zone.
This phenomenon occurs also in valve controlled pneumatic and hydraulic actuators.
BACKLASH
FRICTION
Friction exists in a system, when there is relative motion between contacting surfaces.
The different types of friction are viscous friction, coulomb friction and stiction. The viscous
fiction is linear in nature and the frictional force is directly proportional to relative velocity of
sliding surfaces.
The coulomb friction is a drag force which opposes motion (For eg: drag due to rubbing contact
between brushes and the commutator).
Due to interlocking of surface irregularities, more force is required to move an object from rest
than to maintain it in motion. Hence, force of stiction is always greater than that of Coulomb
friction.
Describing Function (DF) method provides a linear approximation to the nonlinear element based
on the assumption that the input to nonlinear element is a sinusoidal input of known constant amplitude.
The fundamental harmonic element's output is compared with input sinusoidal to determine steady state
amplitude and phase relation. This relation is the describing function for the nonlinear element. Thus, the
method can also be viewed as 'harmonic linearization' of a non- linear element. DF method is based on
'Fourier series. DF method predicts whether limit cycle oscillation exist or not and gives numerical
estimates of amplitude and frequency of oscillation, when limit cycles are predicted.
To discuss the basic concept of DF analysis, consider the block diagram of a nonlinear system
shown in figure below, where blocks 𝐺1 (𝑠)) and 𝐺2 (𝑠) represent linear elements while the block N
represents nonlinear element.
Let 𝑥 = 𝑋 sin 𝜔𝑡 be a sinusoidal input given to the nonlinear system. The output ‘y' of the nonlinear
element will be a non-sinusoidal periodic function, which may be expressed in terms of Fourier series.
If the nonlinearity is symmetrical the average value of y is zero and hence the output y is given by
In the absence of external input (i.e., when 𝑟 = 0) the output y of the nonlinearity N is feedback to
input through 𝐺1 (𝑠) and 𝐺2 (𝑠) in tandem. If 𝐺1 (𝑠) 𝐺2 (𝑠) has low pass characteristics, then all the
harmonics of y are filtered, so that the input 𝑥 to the nonlinear element N is mainly contributed by
fundamental component of 𝑦 and hence 𝑥 remains sinusoidal. Under such conditions the harmonics of the
output are neglected and the fundamental component of 𝑦 alone considered for the purpose of analysis.
𝑦 = 𝑦1 = 𝐴1 sin 𝜔𝑡 + 𝐵1 cos 𝜔𝑡 = 𝑌1 ∅1 = 𝑌1 sin(𝜔𝑡 + ∅)
𝐵
Where, Y1 = √A21 + 𝐵12 and ∅1 = tan−1 𝐴1
1
∅1 =Phase shift of the fundamental harmonic component of the output with respect to the input.
When the input 𝑥 to the nonlinearity is sinusoidal, the describing function of the nonlinearity is
defined as,
𝑌1
𝐾𝑁 (𝑋, 𝜔) = ∅1
𝑋
The nonlinear element N in the system can be replaced by the following describing function.
[Link]’s College of Engineering 68 VI Semester ICE
IC6601 Advanced Control System Department of ICE 2018-19
If the nonlinearity is replaced by the describing function then all linear theory frequency domain
techniques can be used for the analysis of the system. The describing functions are used only for stability
analysis and it is not directly applied to the optimization of the system design. The describing function is a
frequency domain approach and no general correlation is possible between time and frequency response.
Saturation Nonlinearity
The input-output relation is linear for 𝑥 = 0 to 𝑆. When the input 𝑥 > 𝑆, the output reaches a saturated value of 𝐾𝑆.
The input is a sinusoidal,
𝑥 = 𝑋 sin 𝜔𝑡 (1)
Where 𝑋 is the maximum value of input.
When 𝜔𝑡 = 𝛽, 𝑥 = 𝑆
Hence equation (1) can be written as,
𝑆 = 𝑋 sin 𝛽 (2)
𝑆 −1 𝑆
sin 𝛽 = 𝑜𝑟 𝛽 = sin ( ) (3)
𝑋 𝑋
The output 𝑦 of the nonlinearity can be divided into three regions in a period of 𝜋. The output equation for the three
regions are
𝐾𝑥 ; 0 ≤ 𝜔𝑡 ≤ 𝛽
𝑦 = { 𝐾𝑆 ; 𝛽 ≤ 𝜔𝑡 ≤ (𝜋 − 𝛽) (4)
𝐾𝑥 ; (𝜋 − 𝛽) ≤ 𝜔𝑡 ≤ 𝜋
2𝐾 2
−1 𝑆 𝑆
√1 − ( 𝑆 )
∴ 𝐾𝑁 (𝑋, 𝜔) = 𝜋
[sin (𝑋
) + (𝑋
)] 𝑋
𝑓𝑜𝑟 𝑋 > 𝑆 (17)
Backlash Nonlinearity
The output equation has no symmetry so that its Fourier series contains both sine and cosine terms. The average
value of the output is zero and hence 𝐴0 = 0
2 𝜋
𝐴1 = 𝜋 ∫0 𝑦 cos 𝜔𝑡 𝑑(𝜔𝑡)
2 𝜋+𝛽 −4𝑁
𝐴1 = 𝜋 ∫𝛽 𝑁 cos 𝜔𝑡 𝑑(𝜔𝑡) = 𝜋
sin 𝛽 (3)
2 𝜋
𝐵1 = 𝜋 ∫0 𝑦 sin 𝜔𝑡 𝑑(𝜔𝑡)
2 𝜋+𝛽
= 𝜋 ∫𝛽 𝑁 sin 𝜃 𝑑𝜃 (4)
Let 𝑌1 = Amplitude of the fundamental harmonic component of the output
∅1 = Phase shift of the fundamental harmonic component of the output with respect to the I input.
𝑌
The describing function is given by, 𝐾𝑁 (𝑋, 𝜔) = 𝑋1 ∅1
𝐵
Where 𝑌1 = √𝐴12 + 𝐵12 and ∅1 = tan−1 (𝐴1 )
1
The output y has half wave and quarter wave symmetries
4𝑁 4𝑁 −𝑗𝛽
(cos 𝛽 − 𝑗 𝑠𝑖𝑛𝛽) = 𝑒
𝜋𝑀 𝜋𝑀
4𝑁
𝜋𝑀
−𝛽 ; 𝑀 ≥𝑎 (5)
Describing functions of nonlinear elements can be used for stability analysis of nonlinear control
systems. Also it is used to predict the sustained oscillations or limit cycles in the output of the system.
Consider a unity feedback system shown, in which the nonlinearity is represented by its describing function
𝐾𝑁 (𝑋, 𝜔) or 𝐾𝑁 . Let 𝐶(𝑗𝜔)/𝑅(𝑗𝜔) be the closed loop sinusoidal transfer function of the system
𝐶(𝑗𝜔) 𝐾𝑁 G(𝑗𝜔)
=
𝑅(𝑗𝜔) 1 + 𝐾𝑁 (𝑗𝜔)
The characteristic equation of the system is obtained by equating the denominator to zero.
1 + 𝐾𝑁 𝐺(𝑗𝜔) = 0
The Nyquist stability criteria can also be extended to the stability analysis of nonlinear systems. According
to the Nyquist stability criterion the system will exhibit systained oscillations or limit cycles when
𝐾𝑁 𝐺(𝑗𝜔) = −1
This implies that the sustained oscillations or limit cycles will occur if 𝐾𝑁 𝐺(𝑗𝜔) locus pass through the
critical point, −1 + 𝑗0 in the complex plane. Above equation can be written as
1
𝐺(𝑗𝜔) = −
𝐾𝑁
This implies that the critical point, −1 + 𝑗0 becomes the critical locus which is the locus of −1/𝐾𝑁 . Hence
the intersection point of 𝐺(𝑗𝜔) locus and −1/𝐾𝑁 locus will give the amplitude and frequency of limit cycle.
To determine the stability of the system due to nonlinearity sketch the −1/𝐾𝑁 locus and 𝐺(𝑗𝜔) locus (polar
plot) in complex plane. From the sketches the following conclusions can be obtained.
1. If the −1⁄𝐾 locus is not enclosed by the 𝐺(𝑗𝜔) locus then the system is stable or there is no limit
𝑁
cycle at steady state.
2. If the −1⁄𝐾 locus is enclosed by the 𝐺(𝑗𝜔) locus then the system is unstable.
𝑁
3. If the −1⁄𝐾 locus and the 𝐺(𝑗𝜔) locus intersect, then the system output may exhibit a sustained
𝑁
oscillation or a limit cycle. The amplitude of the limit cycle is, given by the value of −1⁄𝐾 locus
𝑁
at the intersection point. The frequency of the limit cycle is given by the frequency of G(jω)
corresponding to the intersection point.
CONCEPT OF ENCLOSURE
In a complex plane the −1/𝐾𝑁 locus is said to be enclosed by𝐺(𝑗𝜔) locus if it lies in the region to the
right of an observer travelling through 𝐺(𝑗𝜔) locus in the direction of increasing .
[Link]’s College of Engineering 86 VI Semester ICE
IC6601 Advanced Control System Department of ICE 2018-19
In a complex plane −1/𝐾𝑁 locus is not enclosed by 𝐺(𝑗𝜔) if it lies in the region to the left of an
observer travelling through 𝐺(𝑗𝜔) locus in the direction of increasing .
If the −1/𝐾𝑁 l ocus and 𝐺(𝑗𝜔) locus intersect as shown in figure, then for an observer
travelling through 𝐺(𝑗𝜔) locus in the direction of increasing , the region on the Right is unstable
region and the region on the left is stable region .
Problem : A servo system used for positioning a load has backlash characteristics as shown in figure
1. The block diagram of the system is shown in figure 2. The magnitude and phase of the describing
function of backlash nonlinearity for various values of b/X are listed in table-1, where X = Maximum
value of input sinusoidal signal to the nonlinearity.
Show that the system is stable if K=1. Also show that limit cycle exists when K=2. Investigate the
stability of these limit cycles and determine their frequency and b/X.
Solution
The magnitude of 𝐺(𝑗𝜔) in 𝑑𝑏 and phase of 𝐺(𝑗𝜔) are calculated for various values of 𝜔 and
tabulated. The Nichols pot of 𝐺(𝑗𝜔) is sketched in an ordinary graph sheet.
When K=2 the magnitude of G(jω) increases by an amount 20 log2=6db. The phase of G(jω) is not
altered. The increase in magnitude is independent of frequency. Hence, G(jω) locus when K=2 is obtained
by shifting the locus of G(jω) when K=1 , by 6db upwards.
The magnitude and phase of the describing function of backlash, KN is listed in the problem in table for
various values of b/X. Using these values of KN and KN the values of 1/KN in db and 1/KN are
calculated for various values of b/X are listed in table. Using these values the locus of -1/Kn is sketched.
STABILITY ANALYSIS
Case (i) when K = 1
From the Nichols plots it is observed that when K = 1, 𝐺(𝑗𝜔) locus does not enclose −1/𝐾𝑛 locus. Hence
the system is stable.
From the Nichols plots it is observed that when K = 2; 𝐺(𝑗𝜔) locus intersects −1/𝐾𝑛 locus at two points.
At one intersection point unstable limit cycle exists and at another intersection point stable limit cycle
exists.
RESULT
• The unstable limit cycle exists when 𝑏⁄𝑋 = 0.316 and the frequency of oscillation is 0.36 rad/sec.
• The stable limit cycle exists when 𝑏⁄𝑋 = 0.464 and the frequency of oscillation is 1.07 rad/sec.
For the nonlinear system given below, investigate the possibility of a limit cycle. If limit cycle exists, find
its magnitude and frequency. (May 2017)
Solution
10
Given 𝑁 = 1, 𝐺(𝑠) = 𝑠(𝑠+1)(𝑠+2)
4𝑁
The describing function of the nonlinear element is 𝐾𝑁 = 𝜋𝑀
1 𝜋𝑀 𝜋𝑀 𝜋𝑀
− =− =− =−
𝐾𝑁 4𝑁 4×1 4
𝟏
Locus of −
𝑲𝑵
1
When 𝑀 → 0, − 𝐾 → 0
𝑁
1
When M → ∞, − K → ∞
N
The intersection point P can be determined by solving the equation. At the intersection point,
∴ |G(jω)| = 1 and G(jω) = −180°
Let, 𝜔𝑙1 = frequency when G(jω) = −1
At 𝜔 = 𝜔𝑙1 ,
G(jω) = −90° − tan−1 𝜔𝑙1 − tan−1 0.5𝜔𝑙1 = −180°
− tan−1 𝜔𝑙1 − tan−1 0.5𝜔𝑙1 = −180° + 90°
tan−1 𝜔𝑙1 + tan−1 0.5𝜔𝑙1 = 90°
On taking 𝑡𝑎𝑛 on either side,
𝑡𝑎𝑛(tan−1 𝜔𝑙1 + tan−1 0.5𝜔𝑙1 ) = tan 90°
tan(tan−1 𝜔𝑙1 ) + tan(tan−1 0.5𝜔𝑙1 )
= tan 90°
1 − tan(tan−1 𝜔𝑙1 ) × tan(tan−1 0.5𝜔𝑙1 )
𝜔𝑙1 + 0.5𝜔𝑙1 1
=∞=
1 − 𝜔𝑙1 × 0.5𝜔𝑙1 0
2 2 2
1
1 − 𝜔𝑙1 × 0.5𝜔𝑙1 = 0; 1 − 0.5𝜔𝑙1 = 0; 0.5𝜔𝑙1 = 1; 𝜔𝑙1 = ; 𝜔𝑙1 = √2 𝑟𝑎𝑑/𝑠𝑒𝑐
0.5
The frequency of limit cycle = √2 𝑟𝑎𝑑/𝑠𝑒𝑐
At 𝜔 = 𝜔𝑙1 , |G(jω)| = 1
5
|G(jω)| = 2
ω √1+𝜔 √1+0.25𝜔2
5
|G(jω)| =
ω √1+𝜔2 √1+0.25𝜔2
5 5 5
= 2 ×√3×1.5
= 3 = 1.66
2 2 √
√2 √1+(√2) √1+0.25(√2)
1
The 𝐺(𝑗𝜔) locus intersects the − K at (−1.66,0). Limit cycle exists at the intersection point
N
𝜋𝑀
Amplitude of limit cycle= − 4
= −1.66
Hence 𝑀 = 2.1146. the limit cycle is −2.1146 sin √2 𝑡
Unit IV
Introduction to Optimal Control
INTRODUCTION
The sources available on the planet earth are being consumed at an alarming rate and in order to
preserve whatever is left for the future generations, it is essential that consumption has to be reduced to
the minimum possible. The situation applies to almost all resources, energy, fuel, food and so on. Any
system that is designed in these years should confirm to stringent conditions on energy, efficiency and
quality to avoid unnecessary wastage of resources. For this, it is essential that design is made according
to strict constraints rather than trial and error approach as followed in the early development of control
systems. The requirements on the system formulated as specifications should be satisfied and all the
physical and mathematical constraints met simultaneously. This requires formulation of the design problem
as maximization or minimization of a performance measure subject to the set of constraints. Such a design
is said to be an optimal design.
In the classical control using transfer function and frequency domain approach some attempts were
made to arrive at an optimal design. The performance measures were based on the integral of the error
and its adaptations. Some of the performance measures used in classical control are the Integral of Error
(IE), Integral of the Absolute Error (IAE), Integral of the Squared Error (ISE) and Int egral of Time
Absolute Error (ITAE) criteria are discussed below.
∞
J1 = ∫ |e(t)|dt
0
Positive and negative errors are given equal weightage in this measure.
3. Integral of the Squared Error (ISE)
∞
J1 = ∫ e(t)2 dt
0
Penalty for small errors is significantly reduced as t he square of small quantities will be still smaller in
relation to the square of large errors which becomes relatively larger when squared.
4. Integral of the Time Absolute Error (ITAE)
∞
J1 = ∫ t e(t)dt
0
∞
J1 = ∫ t e(t)2 dt
0
Among these the 1SE criterion has been reasonably well studied due to this being a quadratic function and
the solution is mathematically elegant and tractable.
15. b. For the system of the following figure, compute the value of 𝐾 that minimizes 𝐼𝑆𝐸 for
the unit step input. (May 2013)
Problem 1: Consider the feedback system shown. Compute the value of K that minimizes ISE for the
unit step input.
Solution
For a continuous time system, the state equation in general can be written as
ẋ = f(x(t), u(t), t), t ∈ [t 0 , t i ]
and for a discrete time system
x(k + 1) = f(x(k), u(k), k), k ∈ [k 0 , k j ]
The problem is to control the plant so as to achieve certain specified objectives while satisfying the
operational constraints. The steps involved in the solution of the optimal control problem are:
1. Compute the performance index 𝐽 as a function of the free parameters 𝑘1 , 𝑘2 , … .. , 𝑘𝑛 of the system
with fixed configuration:
𝐽 = 𝐽(𝑘1 , 𝑘2 , . … .. , 𝑘𝑛 )
2. Determine the solution set 𝑘𝑖 of the equations
𝜕𝐽
= 0; 𝑖 = 1, 2, … … … , 𝑛
𝜕𝑘𝑖
This gives the necessary conditions for J to be minimum. From the solution set of these equations, find
the subset that satisfies the sufficient conditions which require that the Hessian matrix is positive
definite.
𝜕 2𝐽 𝜕 2𝐽 𝜕 2𝐽
,,,,,,
𝜕𝑘12 𝜕𝑘1 𝜕𝑘2 𝜕𝑘1 𝜕𝑘𝑛
𝜕 2𝐽 𝜕 2𝐽 𝜕 2𝐽
𝐻 = 𝜕𝑘2 𝜕𝑘1 ,,,,,,
𝜕𝑘22 𝜕𝑘2 𝜕𝑘𝑛
; ; ;
,,,,,,
𝜕 2𝐽 𝜕 2𝐽 𝜕 2𝐽
,,,,,,
[𝜕𝑘𝑛 𝜕𝑘1 𝜕𝑘𝑛 𝜕𝑘2 𝜕𝑘𝑛2 ]
𝜕2 𝐽 𝜕2 𝐽
Since = 𝜕𝑘
𝜕𝑘𝑖 𝜕𝑘𝑗 𝑗 𝜕𝑘𝑖
The design of an optimal controller is based on the following factors relating to the plant and to
the nature of its connection with the controller.
[Link]’s College of Engineering 99 VI Semester ICE
IC6601 Advanced Control System Department of ICE 2018-19
1. Characteristics of the plant.
2. The requirements to be satisfied by the plant.
3. The nature of information supplied to the controller about the plant.
Constraints can be laid on the states 𝑥(𝑡) also during the interval of operation of the optimal control.
A state trajectory that satisfies the state variable constraints during the time interval 𝑡 ∈ [𝑡0 , 𝑡1 ] will be
called the admissible trajectory. The set of admissible state trajectory is denoted by 𝑋. Then 𝑥(𝑡) is
admissible if 𝑥(𝑡) ∈ 𝑋 for all 𝑡 ∈ [𝑡0 , 𝑡1 ].
2. Requirements on the Plant
The designer translates the requirements or specifications to be satisfied by the plant into a
mathematical performance criterion or index and then tries to find a solution that optimises the performance
measure. In certain cases, the performance measure can be chosen without difficulty, but in some other
cases it may have to be selected with difficulty with lot of subjectivity involved. In some cases several
performance measures will have to be attempted to find a satisfactory and suitable solution.
14 a (i) Explain briefly about the minimum-time control of linear time-invariant systems. (8) May 2016
The time required for the solution of the problem is to be minimized similar to reaching particular
destination from a starting point. The starting point is specified as initial state 𝑥(𝑡0 ) = 𝑥0 . The final state
is required to be on as specific region 𝑆 of the state space. 𝑆 may be called a target set. If the final state
is fixed, the target is a straight line. The objective is to transfer the system from the initial stale 𝑥(𝑡0 ) =
𝑥0 to the specified target in the minimum possible time. The performance measure thus becomes:
𝑡1
𝐽 = 𝑡1 − 𝑡0 = ∫ 𝑑𝑡
𝑡0
[Link]’s College of Engineering 100 VI Semester ICE
IC6601 Advanced Control System Department of ICE 2018-19
where 𝑡0 is the first instant of time when 𝑥(𝑡) is in 𝑆 of interest.
The interception of an attacking aircraft by a missile is a typical example. The attacking aircraft
has to be destroyed before it can do any damage to human beings or machinery. This type of problem is
always associated with inequality constraints on the magnitude of the input-output vector elements. In the
absence of any constraints, the change can be made in zero time which is impossible as no physical
device can give infinite energy to a device.
2.1.2 Minimum Energy Problem
The objective is to transfer the system from the initial state 𝑥(𝑡0 ) = 𝑥0 to the specified target with
minimum use of energy during the process. The performance measure thus becomes:
𝑡1
𝐽 = ∫ 𝑢2 (𝑡)𝑑𝑡
𝑡0
If there are several inputs, i.e. 𝑢(𝑡) is a vector, then the performance measure is modified as:
𝑡1
𝐽 = ∫ [𝑢𝑇 (𝑡) 𝑢(𝑡)]𝑑𝑡
𝑡0
If it is required to give different weightage or penalty for different components of the input, we may
introduce a weighting matrix 𝑅 into the performance measure as:
𝑡1
𝐽 = ∫ [𝑢𝑇 (𝑡) 𝑅 𝑢(𝑡)]𝑑𝑡
𝑡0
Here 𝑅 is a diagonal matrix with suitable values on the diagonal to give proper penalty for the different
components of 𝑢.
2.1.3 Minimum Fuel Problem
The quantity of fuel consumed in rockets and propulsion of space crafts is to be optimised to the
maximum extent as the mass of fuel is a significant part of the whole vehicle and any reduction of the fuel
consumed can result in a better utilization of the payload to carry additional equipment for scientific
investigations. The rate of consumption of fuel in jet engines can be taken as directly proportional to the
thrust (input to the space vehicle):
𝑡1
𝐽 = ∫ |𝑢(𝑡)|𝑑𝑡
𝑡0
Allowing the various jets present to different consumption ratio, it can be modified to:
𝑡1
𝐽 = ∫ [𝑘1 |𝑢1 (𝑡)| + 𝑘2 |𝑢2 (𝑡)| + 𝑘3 |𝑢3 (𝑡)|+, … … . . 𝑘𝑝 |𝑢𝑝 (𝑡)|]𝑑𝑡
𝑡0
Here the objective of the optimal control problem is to transfer a system from an initial state 𝑥(𝑡0 ) =
𝑥 0 to the desired state 𝑥(𝑡1 ) = 𝑥1 (may be the equilibrium point) with minimal integral squared error.
[Link]’s College of Engineering 101 VI Semester ICE
IC6601 Advanced Control System Department of ICE 2018-19
1 1)
Relative to the state 𝑥 the quantity (𝑥(𝑡) − 𝑥 can be viewed as instantaneous error. If we transform the
system co-ordinates such that 𝑥1 can be the origin, 𝑥(𝑡) itself will be the error.
𝑡1 𝑛 𝑡1
𝐽 = ∫ [∑(𝑥𝑖 (𝑡)) ] 𝑑𝑡 = ∫ [𝑥 𝑇 (𝑡)𝑥(𝑡)] 𝑑𝑡
2
𝑡0 𝑖=1 𝑡0
where 𝑄 is a real symmetric positive definite matrix, a diagonal matrix in the simplest case. Further, if it
is also desired to minimize the deviation of the final state, the performance index can be modified as:
𝑡1
𝐽 = 𝑥 𝑇 (𝑡1 ) 𝐻 𝑥(𝑡1 ) + ∫ [𝑥 𝑇 (𝑡) 𝑄(𝑥) 𝑥(𝑡)] 𝑑𝑡
𝑡0
Where the first term is introduced to minimize the deviation from the final state and 𝐻 is again a positive
definite or positive semi definite matrix. A more realistic performance measure may also put some penalty
on the energy consumed as well.
𝑡1
𝐽 = 𝑥 𝑇 (𝑡1 ) 𝐻 𝑥(𝑡1 ) + ∫ [𝑥 𝑇 (𝑡) 𝑄(𝑥) 𝑥(𝑡) + 𝑢𝑇 𝑅 𝑢] 𝑑𝑡
𝑡0
In output regulator, the concern is to make the output 𝑦(𝑡) to be as close as possible to the desired
output. The performance measure can be modified as:
1 𝑡1
𝐽 = 𝑦 𝑇 (𝑡1) 𝐻 𝑦(𝑡1 ) + ∫ [𝑦 𝑇 (𝑡) 𝑄(𝑥) 𝑦(𝑡) + 𝑢𝑇 𝑅 𝑢] 𝑑𝑡
2 𝑡0
2.1.6 Tracking Problem
Here the objective is to maintain the system state 𝑥(𝑡) as close as possible to a desired state or
reference state 𝑥(𝑡) in t he interval [𝑡0 , 𝑡1 ]. i.e.
1 𝑇 1 𝑡1
𝐽= 𝑒 (𝑡1 ) 𝐻 𝑒(𝑡1 ) + ∫ [𝑒(𝑡) 𝑄(𝑥) 𝑒(𝑡) + 𝑢𝑇 𝑅 𝑢] 𝑑𝑡
2 2 𝑡0
Where 𝑒(𝑡) = 𝑥(𝑡) − 𝑟(𝑡)
2.2 Choice of Performance Measure
The choice of the performance index or measure is dictated by the physical requirements which
have to be translated to a mathematical function. In some cases very complex functions will be required
to meet all the specifications.
If the optimal control is determined as a function of time for a specified initial state,
𝑢∗ (𝑡) = 𝑓(𝑥(𝑡0 ), 𝑡)
[Link]’s College of Engineering 102 VI Semester ICE
IC6601 Advanced Control System Department of ICE 2018-19
then the optimal control is said to be of the open loop form. This requires no additional information and in
the absence of meaningful disturbances and errors, open loop control can be successful. If a functional
relationship of the form:
𝑢∗ (𝑡) = 𝑓(𝑥(𝑡), 𝑡)
can be found for the optimal control at time 𝑡, then the optimal control is said to be closed loop and
𝑓() is called an optimal control law. Here the controller gets information on the actual state 𝑥(𝑡) of the
process and will be more effective. The effect of disturbances and errors can be significantly reduced in
closed loop control.
Definition l: The optimal control problem can thus be formally defined as: Find an admissible control 𝑢∗ (𝑡)
which causes the system 𝑥̇ (𝑡) = 𝑓(𝑥(𝑡), 𝑢(𝑡), 𝑡) 𝑑𝑡 to follow an admissible trajectory 𝑥 ∗ (𝑡) and minimizes
the performance measure:
𝑡1
𝐽 = ℎ(𝑥(𝑡1 ), 𝑡1 ) + ∫ 𝑔(𝑥(𝑡), 𝑢(𝑡), 𝑡) 𝑑𝑡
𝑡0
Several approaches are available for the solution of the optimal control problem based on the nature of
the performance measure and constraints. Some of the commonly used methods are:
1. Method of calculus of variations.
2. Maximum principle due to Pontryagin.
3. Dynamic Programming approach of Bellman,
4. Lyapunov's method.
5. Linear Programming approach
6. Integer programming approach etc.
15. b. Derive the Ricatti’s equation as a solution of steady state optimal control problem. (Nov
2015)
Quadratic performance measures have been extensively used in the design of optimal control
systems because of its mathematical tractability and the ability to accommodate realistic functions such
as energy function and squared errors or deviations. When squared errors are used as performance
measures, large errors are given a much higher penalty than small errors as the square of large quantities
are very large and squares of small quantities are smaller.
Consider the system described by
𝑥̇ = 𝐴𝑥 + 𝐵𝑢 (1)
where x is the n × 1 state vector, 𝑢 is the 𝑝 × 1 input vector, A and B are respectively, 𝑛 × 𝑛 and 𝑛 × 𝑝
real constant matrices. Find the control law,
𝑢 = −𝐾𝑥 (2)
[Link]’s College of Engineering 103 VI Semester ICE
IC6601 Advanced Control System Department of ICE 2018-19
Where 𝐾 is 𝑝 × 𝑛 real constant unconstrained gain matrix, that minimizes the following performance index
subject to the initial conditions 𝑥(0) ≜ 𝑥 0 :
1 ∞
𝐽 = 2 ∫0 [𝑥 𝑇 𝑄 𝑥 + 𝑢𝑇 𝑅 𝑢] 𝑑𝑡 (3)
where, 𝑄 is a positive definite, real, symmetric, constant matrix and 𝑅 is the 𝑝 × 𝑝 positive definite, real,
symmetric, constant matrix.
There are several ways to solve this optimal control problem and here we use Laypunov approach.
Substitute equation (2) in (1)
𝑥̇ = 𝐴𝑥 − 𝐵𝐾𝑥 = (𝐴 − 𝐵𝐾) 𝑥 (4)
Since, (𝐴, 𝐵) pair is completely controllable, there exists a feedback matrix 𝐾 such that (𝐴 − 𝐵𝐾) is a
stable matrix. The controllability of the given plant, thus, ensures the existence of a 𝐾 that minimizes 𝐽.
The time derivative of Laypunov function is
1
𝑉̇ (𝑥) = − 2 𝑥 𝑇 (𝑄 + 𝐾 𝑇 𝑅𝐾)𝑥 (5)
Since feedback matrix 𝐾 is unconstrained, the optimum value of 𝐽 is independent of initial conditions. The
optimal 𝐾𝑖𝑗 ’s are obtained from the equations,
𝑃
= 0, 𝑓𝑜𝑟 𝑎𝑙𝑙 𝑖, 𝑗 (9)
𝑘𝑖𝑗
1 ∞
• Performance index 𝐽 = 2 ∫0 (𝑥 𝑇 𝑄𝑥 + 𝑢𝑇 𝑅𝑢)𝑑𝑡
14. a (ii) Explain the properties and use of LQR. (8) (May 2016)
Robustness: The LQR achieves infinite gain margin: 𝑘𝑔 = ∞, implying that the loci of (PC) (scalar case)
or (det(1 + 𝑃𝐶) − 1) (MIMO case) approach the origin along the imaginary axis. The LQR also
guarantees phase margin 𝛾 = 60 degrees. This is in good agreement with the practical guidelines for control
system design.
Output Variables: In many cases, it is not the states 𝑥 which are to be minimized, but the output variables
𝑦. In this case, we set the weighting matrix 𝑄 = 𝐶 𝑇 𝑄 ′ 𝐶, since 𝑦 = 𝐶𝑥, and the auxiliary matrix 𝑄 ′ weights
the plant output.
Behavior of Closed-Loop Poles: Expensive Control. When 𝑅 ≫ 𝐶 𝑇 𝑄 ′ 𝐶R the cost function is dominated
by the control effort 𝑢, and so the controller minimizes the control action itself. In the case of a completely
stable plant, the gain will indeed go to zero, so that the closed-loop poles approach the open-loop plant
poles in a manner consistent with the scalar root locus.
The optimal control must always stabilize the closed-loop system, however, so there should be some
account made for unstable plant poles. The expensive control solution puts stable closed-loop poles at the
mirror images of the unstable plant poles.
The Butterworth pattern refers to an arc in the stable left-half plane. The angular separation of n closed-
loop poles on the arc is constant, and equal to 180°/𝑛. An angle 90°/𝑛 separates the most lightly-damped
poles from the imaginary axis.
𝒙̇ 𝟎 𝟏 𝒙𝟏 𝟎
Problem 2. Consider a double integrator plant having state model [ 𝟏 ] = [ ] [𝒙 ] + [ ] 𝒖; 𝑪 =
𝒙̇ 𝟐 𝟎 𝟎 𝟐 𝟏
[𝟏 𝟎] . It is desired to design an optimal state feedback controller of the form 𝒖(𝒕) = −𝑲𝒙(𝒕) such
∞ 𝟏 𝟎
that the performance measure 𝑱 = ∫𝟎 [𝒙𝑻 𝑸𝒙 + 𝒖𝑻 𝒖]𝒅𝒕 where 𝑸 = [ ] , 𝑹 = 𝟏 is minimized.
𝟎 𝟎
SOLUTION
0 1 0 1 0 𝑝1 𝑝2
Given 𝐴 = [ ];𝐵 = [ ]; 𝑄 = [ ] ; 𝑅 = 1; 𝑃 = [𝑝 𝑝3 ]
0 0 1 0 0 2
0 0 𝑝1 𝑝2 𝑝1 𝑝2 0 1 𝑝1 𝑝2 0 𝑝1 𝑝2 1 0 0 0
[
1
][
0 𝑝2 𝑝3 ] + [𝑝2 𝑝3 ] [0 0
] − [𝑝
2 𝑝3 ] [1] [1] [0 1] [𝑝2 𝑝3 ] + [0 0] = [0 0
]
0 0 0 0 𝑝1 𝑝2 1 0 𝑝1 𝑝2 1 0 0 0
[ ]+[ ] − [𝑝 𝑝3 ] [0 0] [𝑝2 𝑝3 ] + [0 0] = [0 0]
𝑝1 𝑝2 0 𝑝2 2
0 0 𝑝 2 𝑝1 𝑝2 1 0 0 0
[ ]− [ 1 2 ] + [0 ]=[ ]
𝑝1 2𝑝2 𝑝1 𝑝2 𝑝2 0 0 0
Expanding
1 − 𝑝22 = 0 𝑝2 = 1
𝑝1 − 𝑝2 𝑝3 = 0
2𝑝2 − 𝑝32 = 0
definite.
𝐾 = [1] [0 √2 1 ] = [
1] [ 1 √2]
1 √2
Which is independent of any initial state under the given performance index. Under state feedback, the
system matrix becomes:
0 1
𝐴 − 𝐵𝐾 = [ ]
−1 −√2
The characteristic equation of the closed loop system is :
2 + √2 + 1 = 0
√2
Which has the roots at = (−1 ± 𝑗) with damping ratio of 0.707 which gives a good performance.
2
∞
15 a. (ii) Derive the control law which minimizes the performance index 𝐽 = ∫0 (𝑥12 + 𝑢2 ) 𝑑𝑡
for the system as follows. (Nov 2014)
𝑥̇ 1 0 1 𝑥1 0
[ ]=[ ][ ] + [ ]𝑢
𝑥̇ 2 0 0 𝑥2 1
1
(𝑠 + )
𝐾∅(𝑠)𝐵 = 𝐾(𝑠𝐼 − 𝐴)−1 𝐵 = √2 √2
𝑠2
∞
Problem 3. Obtain the control law which minimize the performance index 𝑱 = ∫𝟎 (𝒙𝟐𝟏 + 𝒖𝟐 )
𝟎 𝟏 𝟎
for the system 𝑿̇ = [ ]𝑿 + [ ]𝒖
𝟎 𝟎 𝟏
SOLUTION
0 1 0
Given, 𝐴 = [ ] ;𝐵 = [ ]
0 0 1
Modify the equation to
1 ∞ 𝑇
𝐽= ∫ [𝑥 𝑄 𝑥 + 𝑢𝑇 𝑅 𝑢] 𝑑𝑡
2 0
1 ∞ 2
𝐽 = ∫ [2𝑥1 + 2𝑢2 ] 𝑑𝑡
2 0
On comparing the above equations,
2 0
𝑄=[ ] ; 𝑅 = [2]
0 0
The reduced matrix Riccati equation is
𝐴𝑇 𝑃 + 𝑃𝐴 − 𝑃𝐵𝑅 −1 𝐵 𝑇 𝑃 + 𝑄 = 0
0 1 𝑇 𝑝11 𝑝12 𝑝11 𝑝12 0 1 𝑝11 𝑝12 0 1 𝑝11 𝑝12 2 0 0 0
[
0
] [
0 𝑝12 𝑝22 ] + [𝑝12 𝑝22 ] [0 0] − [𝑝12 𝑝22 ] [1] [2] [0 1] [𝑝
12 𝑝22 ] + [0 0
]=[
0 0
]
𝑝12
0
[
0 𝑝11
]−[ 2 ] [𝑝11 𝑝12 2 0 0 0
𝑝11 2𝑝12 𝑝22 𝑝12 𝑝22 ] + [0 0
]=[
0 0
]
0
2
2
𝑝12 𝑝12 𝑝22
0 𝑝11 0+ 0+
2 2 2 0 0 0
[ ]− 2 +[ ]=[ ]
𝑝11 2𝑝12 𝑝12 𝑝22 𝑝22 0 0 0 0
[ 0 + 0+
2 2 ]
Expanding the above,
2
𝑝12
0− +2=0 (1)
2
𝑝12 𝑝22
𝑝11 − +0=0 (2)
2
𝑝12 𝑝22
𝑝11 − +0=0 (3)
2
2
𝑝22
2𝑝12 − +0=0 (4)
2
From (1),
2
𝑝12 2
− + 2 = 0; 𝑝12 = 4; 𝑝12 = 2
2
From (4)
2
𝑝22
2𝑝12 − =0
2
2
𝑝22 2 2
= 2𝑝12 ; 𝑝22 = 4𝑝12 ; 𝑝22 = 8; 𝑝12 = 2√2
2
From (2)
𝑝12 𝑝22
𝑝11 − =0
2
𝑝12 𝑝22 2 × 2√2
𝑝11 = ; = 2√2
2 2
𝑝11 𝑝12
] = [2√2 2 ]
𝑃 = [𝑝
12 𝑝22 2 2√2
1 2√2 2 𝑥 (𝑡)
= − [2] [0 1] [ ][ 1 ]
2 2√2 𝑥2 (𝑡)
2√2 2 𝑥 (𝑡)
= − [0 1⁄2] [ ][ 1 ]
2 2√2 𝑥2 (𝑡)
𝑥 (𝑡)
=−[1 √2] [ 1 ]
𝑥2 (𝑡)
𝑢∗ = −𝑥1 (𝑡) − √2𝑥2 (𝑡)
∴ 𝐾 = [1 √2]
where, 𝑄 is a positive definite, real, symmetric, constant matrix and 𝑅 is the 𝑝 × 𝑝 positive definite, real,
symmetric, constant matrix.
With the control law, the closed loop system becomes
𝑥̇ = 𝐴𝑥 − 𝐵𝐾𝑥 = (𝐴 − 𝐵𝐾) 𝑥 (4)
Since, (𝐴, 𝐵) pair is completely controllable, there exists a feedback matrix 𝐾 such that (𝐴 − 𝐵𝐾) is a
stable matrix. The controllability of the given plant, thus, ensures the existence of a 𝐾 that minimizes 𝐽.
The optimization problem is to determine the values of free parameters of the matrix 𝐾 so as to minimize
the performance index.
Substituting the control vector 𝑢 from equation (2) in the performance index 𝐽 of equation (3)
1 ∞
𝐽 = 2 ∫0 [𝑥 𝑇 𝑄 𝑥 + 𝑥 𝑇 𝐾 𝑇 𝑅 𝐾𝑥] 𝑑𝑡
1 ∞
𝐽 = 2 ∫0 𝑥 𝑇 (𝑄 + 𝐾 𝑇 𝑅 𝐾)𝑥 𝑑𝑡 (5)
The value of the performance index for the system trajectory starting at 𝑥(0) is 𝑉(𝑥(0))
[Link]’s College of Engineering 109 VI Semester ICE
IC6601 Advanced Control System Department of ICE 2018-19
The time derivative of Laypunov function is
1 ∞
𝑉̇ (𝑥) = 2 𝑥 𝑇 (𝑄 + 𝐾 𝑇 𝑅𝐾)𝑥|
𝑡
1 𝑇 1
= 𝑥 (∞)[𝑄 + 𝐾 𝑇 𝑅 𝐾] 𝑥(∞) − 𝑥 𝑇 (𝑡)[𝑄 + 𝐾 𝑇 𝑅 𝐾] 𝑥(𝑡)
2 2
The matrix (𝐴 − 𝐵𝐾) is stable. From equation (4)
𝑥(∞) → 0
1
Therefore, 𝑉̇ (𝑥) = − 2 𝑥 𝑇 (𝑄 + 𝐾 𝑇 𝑅𝐾)𝑥 (6)
Since 𝑉̇ (𝑥) is a quadratic in 𝑥 and the plant equation in linear, let us assume that 𝑉(𝑥) is also given by
quadratic form:
1
𝑉(𝑥) = 2 𝑥 𝑇 𝑃 𝑥 (7)
This equation is of the form of Lyapunov equation. Equation (9) has 𝑛2 nonlinear algebraic equation. Since
𝑛(𝑛+1)
𝑛 × 𝑛 matrix 𝑃 is symmetric, we need to solve only equations for the elements 𝑃𝑖𝑗 as functions of
2
𝑉(𝑥(0)) is the value of performance index for the system trajectory starting at 𝑥(𝑡), from equation (7)
1
𝐽 = 2 𝑥 𝑇 (0) 𝑃 𝑥(0) (10)
A suboptimal control law may be obtained by minimizing 𝐽 with respect to the available elements 𝑘𝑖𝑗 of 𝐾,
[𝑥 𝑇 (0)𝑃 𝑥(0)]
=0
𝑘𝑖𝑗
If for the suboptimal solution thus obtained, the matrix (A − BK) is stable, then the minimization of J as
the procedure described above gives the correct result.
Additional Problems
Problem 4. Consider the second order system shown in figure. Determine the value of damping ratio
> 𝟎 which minimizes the integral square error
∞
𝑱 = ∫ 𝒆𝟐 (𝒕)𝒅𝒕
𝟎
For the initial conditions 𝒚(𝟎) = 𝟏, 𝒚̇ (𝟎) = 𝟎. What is the minimum value of the performance
index?
Solution
The differential equation governing the system
1⁄
𝑌(𝑠) 𝐺(𝑠) 𝑠(𝑠 + 2) 1
= = = 2
𝑅(𝑠) 1 + 𝐺(𝑠)𝐻(𝑠) 1 + 1⁄ 𝑠 + 2𝑠 + 1
𝑠(𝑠 + 2)
[Link]’s College of Engineering 111 VI Semester ICE
IC6601 Advanced Control System Department of ICE 2018-19
𝑅(𝑠) = 0 ; ∴ 𝑌(𝑠)[𝑠 + 2𝑠 + 1] = 0
2
1 1 1 1
𝐽= [2 + 1] [ ] = +
2 2 0 4
To find the minimum 𝐽, differentiate and w.r.t and equate to 0
𝜕𝐽 1 𝜕 2𝐽
= 1 − 2 = 0; = 0.5; 2 > 0; 𝐽𝑚𝑖𝑛 = 1
𝜕 4 𝜕
Problem 5. A unity feedback system has the forward path transfer function
𝑲
𝑮(𝒔) =
𝒔(𝒔 + 𝜶)
The reference input is a unit step function, and the system is at rest initially.
(i) Taking K as constant, determine the value of 𝜶 which minimizes the integral square error.
(ii) Taking 𝜶 as constant, determine the value of K which minimizes the integral square error.
Solution
The closed loop transfer function,
𝑌(𝑠) 𝐾
= 2
𝑅(𝑠) 𝑠 + 𝛼𝑠 + 𝐾
𝑌(𝑠)[𝑠 2 + 𝛼𝑠 + 𝐾] = 𝐾 𝑅(𝑠)
Take inverse Laplace Transform,
𝑦̈ + 𝛼𝑦̇ + 𝐾 𝑦 = 𝐾 𝑟 (1)
𝑒 = 𝑟 − 𝑦; 𝑒̇ = 𝑦̇ ; 𝑒̈ = 𝑦̈ (2)
Substitute (2) in (1), 𝑒̈ + 𝛼𝑒̇ + 𝐾 𝑒 = 0 (3)
Assume the state variable as, 𝑥1 = 𝑒; 𝑥2 = 𝑒̇
Equation (3) becomes, 𝑥̇ 2 + 𝛼𝑥2 + 𝐾 𝑥1 = 0
𝑥̇ 2 = −𝛼𝑥2 − 𝐾 𝑥1 and 𝑥̇ 1 = 𝑥2
The state equation,
𝑥̇ 0 1 𝑥1
𝑥̇ = [ 1 ] = [ ][ ];
𝑥̇ 2 −𝐾 −𝛼 𝑥2
0 1 1
𝐴=[ ] 𝑎𝑛𝑑 𝑥(0) = [ ]
−𝐾 −𝛼 0
Modify the equation to
[Link]’s College of Engineering 113 VI Semester ICE
IC6601 Advanced Control System Department of ICE 2018-19
∞
1
𝐽= ∫ [𝑥 𝑇 𝑄 𝑥 + 𝑢𝑇 𝑅 𝑢] 𝑑𝑡
2 0
The performance index,
∞ 𝟏 ∞
𝑱 = ∫𝟎 𝒆𝟐 (𝒕)𝒅𝒕 = 𝟐 ∫𝟎 𝟐𝒙𝟐𝟏 (𝒕)𝒅𝒕
2 0
From the above equation, 𝑄 = [ ] ; 𝑅 = 0; 𝐾 = 0
0 0
The Lyapunov equation, [𝑨 − 𝑩𝑲]𝑻 𝑷 + 𝑷[𝑨 − 𝑩𝑲] + 𝑸 + 𝑲𝑻 𝑹𝑲 = 𝟎
Since, 𝐾 = 0, the Lyapunov equation becomes,
∴ 𝑨𝑻 𝑷 + 𝑷𝑨 + 𝑸 = 𝟎
1 𝛼 1 1 1 𝛼 1
𝐽= [ + ] [ ]= +
2 𝐾 𝛼 𝐾 0 2𝐾 2𝛼
(ii) Taking 𝜶 as constant, determine the value of K which minimizes the integral square error.
𝜕𝐽 𝛼 1
= (− 2 ) = 0; 𝐾 → ∞
𝜕𝐾 2 𝐾
SOLUTION
Modify the equation to
1 ∞ 𝑇
𝐽 = ∫ [𝑥 𝑄 𝑥 + 𝑢𝑇 𝑅 𝑢] 𝑑𝑡
2 0
1 ∞ T 1 ∞
𝐽= ∫ x x 𝑑𝑡 = ∫ (𝑥12 + 𝑥22 ) 𝑑𝑡
2 0 2 0
On comparing the above equations,
1 0
𝑄=[ ] ; 𝑅 = [0]
0 1
𝑥1
u(t) = −x1 − kx2 = −[1 𝐾 ] [𝑥 ] = −𝐾𝑥
2
0 1 0 1
Given 𝐴 = [ ] ; 𝐵 = [ ] ; 𝑥(0) = [ ] 𝐾 = [1 𝐾 ]
0 0 1 1
The Lyapunov equation, [𝑨 − 𝑩𝑲]𝑻 𝑷 + 𝑷[𝑨 − 𝑩𝑲] + 𝑸 + 𝑲𝑻 𝑹𝑲
Since, 𝑅 = 0, the Lyapunov equation becomes,
[𝑨 − 𝑩𝑲]𝑻 𝑷 + 𝑷[𝑨 − 𝑩𝑲] + 𝑸 = 𝟎
[𝐴 − 𝐵𝐾] = [0 1] − [0] [1 0 1 0 0 0 1
𝐾] = [ ]−[ ]=[ ]
0 0 1 0 0 1 𝐾 −1 −𝐾
𝐾2 + 2 1
𝑝 𝑝12 2𝐾 2]
∴ 𝑃 = [𝑝11 𝑝22 ] = [ 1
12 1
2 𝐾
𝟏
(ii) 𝑱 = 𝟐 𝒙𝑻 (𝟎)𝑷𝒙(𝟎)
𝐾2 +2 1
1 2𝐾 2 1
=2 [1 1] [ 1 1
] [ ]
1
2 𝐾
1 𝐾2 + 2 1 1 1 1 𝐾 2 + 2𝐾 + 4
𝐽= [ + ] [1] =
2 2𝐾 + 2 2 𝐾 4𝐾
𝜕𝐽 𝜕 𝐾 2 + 2𝐾 + 4 3
= ( ) = 0; 𝐾 = 2; 𝐽𝑀𝑖𝑛 =
𝜕𝐾 𝜕𝐾 4𝐾 2
To check the stability: eigen values should be negative. i. e. |I − (A − BK)| = 0
Solving this equation results,
2 + 2 + 1 = 0; 1 = −1; 2 = −1
Hence the system is stable.
To find the sensitivity
For 𝐾 = 1.5, 𝐽 = 1.54
For 𝐾 = 2, 𝐽 = 1.5
[Link]’s College of Engineering 116 VI Semester ICE
IC6601 Advanced Control System Department of ICE 2018-19
∆𝑱⁄ 𝟎.𝟎𝟒⁄
𝑶𝒑𝒕 𝑱 𝟏.𝟒
Sensitivity 𝑺𝑲 = ∆𝑲⁄ = 𝟎.𝟓⁄ = 𝟎. 𝟏𝟎𝟕
𝑲 𝟐
SOLUTION
Modify the equation to
1 ∞ 𝑇
𝐽= ∫ [𝑥 𝑄 𝑥 + 𝑢𝑇 𝑅 𝑢] 𝑑𝑡
2 0
1 ∞ 2
𝐽 = ∫ (𝑥1 + 𝑥22 ) 𝑑𝑡
2 0
On comparing the above equations,
1 0
𝑄=[ ] ; 𝑅 = [0]
0 1
𝑥1
u(t) = −𝑘1 x1 − k 2 x2 = −[𝑘1 𝑘2 ] [𝑥 ] = −𝐾𝑥
2
0 1 0 1
𝐴=[ ] ; 𝐵 = [ ] ; 𝑥(0) = [ ] ; 𝐾 = [𝑘1 𝑘2 ]
0 0 1 0
0 1 0 0 1 0 0 0 1
[𝐴 − 𝐵𝐾] = [ ] − [ ] [𝑘1 𝑘2 ] = [ ]−[ ]=[ ]
0 0 1 0 0 𝑘1 𝑘2 −𝑘1 −𝑘2
0 1
ẋ = [ ] x or 𝑥̇ 1 = 𝑥2 ; 𝑥̇ 2 = −𝑘1 𝑥1 − 𝑘2 𝑥2 ; 𝑥̇ 2 = 𝑥̈ 1
−𝑘1 −𝑘2
Since the undamped natural frequency is specified as 2 𝑟𝑎𝑑/𝑠𝑒𝑐, we obtain 𝑘1 = 4 (compare with
standard second order system. i.e. 𝑥̈ 1 + 2𝜔𝑛 𝑥̇ 1 + 𝜔𝑛2 𝑥1 = 0
0 1
Therefore, [𝐴 − 𝐵𝐾] = [ ]
−4 −𝑘2
The Lyapunov equation, [𝑨 − 𝑩𝑲]𝑻 𝑷 + 𝑷[𝑨 − 𝑩𝑲] + 𝑸 + 𝑲𝑻 𝑹𝑲
Since, 𝑅 = 0, the Lyapunov equation becomes,
[𝑨 − 𝑩𝑲]𝑻 𝑷 + 𝑷[𝑨 − 𝑩𝑲] + 𝑸 = 𝟎
5 𝑘2 1
𝑝 𝑝12
+
2𝑘2 8 8
∴ 𝑃 = [𝑝11 𝑝22 ] =
12 1 5
[ 8 8𝑘2 ]
𝟏 𝑻
𝑱= 𝒙 (𝟎)𝑷𝒙(𝟎)
𝟐
5 𝑘2 1
+
1 2𝑘2 8 8 1
=2 [1 0] [ 1 5 ][ ]
0
8 8𝑘2
1 5 𝑘2 1 1 1 5 𝑘2
𝐽= [ + ] [ ]= ( + )
2 2𝑘2 8 8 0 2 2𝑘2 8
𝜕𝐽 𝜕 1 5 𝑘2 𝜕2 𝐽
= 𝜕𝑘 (2 (2𝑘 + )) = 0; 𝐾 = √20; >0
𝜕𝑘2 2 2 8 𝜕𝑘22
√5
For this value of 𝑘2 , 𝐽𝑀𝑖𝑛 = 4
15. a. Consider the system described by the state model. (May 2013)
𝟎 𝟏 𝟎
𝒙̇ = [ ] 𝒙 + [ ] 𝒖 ; 𝒚 = [𝟏 𝟎]𝒙
𝟎 −𝟐 𝟐𝟎
Find the optimal control law that minimizes
𝟏 ∞
𝑱 = 𝟐 ∫𝟎 (𝒚𝟐 + 𝒖𝟐 ) 𝒅𝒕
SOLUTION
For output regulator problem the performance index is given by,
1 ∞
J = ∫ (𝑦 𝑇 𝑄 𝑦 + 𝑢𝑇 𝑅 𝑢)dt
2 0
Modify the equation to
1 ∞ 𝑇
𝐽 = ∫ [𝑥 𝑄 𝑥 + 𝑢𝑇 𝑅 𝑢] 𝑑𝑡
2 0
Substituting 𝑦 = 𝐶𝑥, the performance index becomes,
1 ∞ 𝑇 𝑇
J = ∫ (𝑥 𝐶 𝑄 𝐶 𝑥 + 𝑢𝑇 𝑅 𝑢)dt
2 0
𝑄 is replace by 𝐶 𝑇 𝑄 𝐶
1 ∞ 𝑇 1
J = ∫ (𝑥 [ ] [1 0]𝑥 + 𝑢𝑇 𝑅 𝑢) dt
2 0 0
1 ∞ 1 0
J = ∫ (𝑥 𝑇 [ ] 𝑥 + 𝑢𝑇 𝑅 𝑢) dt
2 0 0 0
On comparing the above equations,
1 0
𝑄=[ ] ; 𝑅 = [1]
0 0
0 1 0
Given, 𝐴 = [ ] ;𝐵 = [ ]
0 −2 20
The reduced matrix Riccati equation is
𝐴𝑇 𝑃 + 𝑃𝐴 − 𝑃𝐵𝑅 −1 𝐵 𝑇 𝑃 + 𝑄 = 0
2
0 0 0 𝑝11 − 2𝑝12 400𝑝12 400𝑝12 𝑝22 1 0 0 0
[ ]+[ ]− [ ]+[ ]=[ ]
𝑝11 − 2𝑝12 𝑝12 − 2𝑝22 0 𝑝12 − 2𝑝22 400𝑝12 𝑝22 2
400𝑝22 0 0 0 0
Expanding the above,
2
0 − 400𝑝12 +1=0 (1)
2
𝑝12 − 2𝑝22 + 𝑝12 − 2𝑝22 − 400𝑝22 =0 (3)
From (1),
2 2 2
1 1
0 − 400𝑝12 + 1 = 0; 400𝑝12 = 1; 𝑝12 = ; 𝑝12 = = 0.05
400 20
From (3)
2
𝑝12 − 2𝑝22 + 𝑝12 − 2𝑝22 − 400𝑝22 =0
2
2𝑝12 − 4𝑝22 − 400𝑝22 =0
2
400𝑝22 + 4𝑝22 − 2𝑝12 = 0
2
400𝑝22 + 4𝑝22 − 0.1 = 0
= 0.1 + 20 × 0.0116
∴ 𝐾 = [1 0.232]
𝑝12
−𝑝12 −𝑝22 0
[𝑝
−𝑝12 𝑝11 − 2𝑝12 2 ] [𝑝11 𝑝12 2 0 0 0
11 − 2𝑝12 𝑝12 − 2𝑝22 ] + [−𝑝22 𝑝12 − 2𝑝22
]− [ 𝑝22 𝑝12 𝑝22 ] + [0 0
]=[
0 0
]
0
2
2
𝑝12 𝑝12 𝑝22
−𝑝12 −𝑝22 −𝑝12 𝑝11 − 2𝑝12 2 2 2 0 0 0
[𝑝 𝑝12 − 2𝑝22 ] + [−𝑝22 ]− +[ ]=[ ]
11 − 2𝑝12 𝑝12 − 2𝑝22 𝑝12 𝑝22 2
𝑝22 0 0 0 0
[ 2 2 ]
Expanding the above,
2
𝑝12
−𝑝12 − 𝑝12 − +2=0 (1)
2
2
𝑝22
𝑝12 − 2𝑝22 + 𝑝12 − 2𝑝22 − +0=0 (3)
2
2
𝑝12
From (1), −𝑝12 − 𝑝12 − +2=0
2
2
𝑝12
+ 2𝑝12 − 2 = 0
2
Solving the above equation yields, 𝑝12 = −4.828 𝑎𝑛𝑑 0.8284. Since P is a positive matrix, choose
positive value for 𝑝12
𝑖. 𝑒., 𝑝12 = 0.8284
2
𝑝22
From (3), 𝑝12 − 2𝑝22 + 𝑝12 − 2𝑝22 − +0=0
2
2
𝑝22
2𝑝12 − 4𝑝22 − =0
2
2
𝑝22
+ 4𝑝22 − 1.6568 = 0
2
Solving the above equation yields, 𝑝22 = −8.3947 𝑎𝑛𝑑 0.3947. Since P is a positive matrix, choose
positive value for 𝑝22
𝑖. 𝑒., 𝑝22 = 0.3947
From (2)
𝑝12 𝑝22
−𝑝22 + 𝑝11 − 2𝑝12 +
+0=0
2
𝑝12 𝑝22 0.8284 × 0.3947
𝑝11 =− + 2𝑝12 + 𝑝22 ; − + (2 × 0.8284) + 0.3947 = 2.2151
2 2
𝑝11 𝑝12 2.2151 0.8284
𝑃 = [𝑝 𝑝 ]= [ ]
12 22 0.8284 0.3947
15.a. Determine the optimal control law for the system described by (May 2014)
𝟎 𝟏 𝟎
𝒙̇ (𝒕) = [ ] 𝒙(𝒕) + [ ] 𝒖(𝒕),
−𝟐 −𝟑 𝟏
such that the following performance index is minimized.
∞
𝑱 = ∫𝟎 (𝒙𝑻 𝒙 + 𝒖𝟐 ) 𝒅𝒕.
Problem 10. Determine the optimal control law for the system described by
𝟎 𝟏 𝟎
𝐱̇ (𝐭) = [ ] 𝐱(𝐭) + [ ] 𝐮(𝐭)
−𝟐 −𝟑 𝟏
such that the following performance index is minimized.
∞
𝐉 = ∫ (𝐱 𝐓 𝐱 + 𝒖𝟐 )𝐝𝐭
𝟎
SOLUTION
Modify the equation to
1 ∞ 𝑇
𝐽= ∫ [𝑥 𝑄 𝑥 + 𝑢𝑇 𝑅 𝑢] 𝑑𝑡
2 0
1 ∞ 𝑇 1 ∞ 2
𝐽 = ∫ [2𝑥 𝑥 + 2𝑢 ] 𝑑𝑡 = ∫ [2𝑥1 + 2𝑥22 + 2𝑢2 ] 𝑑𝑡
2
2 0 2 0
On comparing the above equations,
2 0
𝑄=[ ] ; 𝑅 = [2]
0 2
0 1 0
Given 𝐴 = [ ] ;𝐵 = [ ] ;
−2 −3 1
The reduced matrix Riccati equation is
𝐴𝑇 𝑃 + 𝑃𝐴 − 𝑃𝐵𝑅 −1 𝐵 𝑇 𝑃 + 𝑄 = 0
0 −2 p11 p12 p11 p12 0 1 p11 p12 0 1 𝑝11 𝑝12 2 0 0 0
[
1
][
−3 p12 p22 ] + [p12 p22 ] [−2 −3] − [p12 p22 ] [1] [2] [0 1] [𝑝12 𝑝22 ] + [0 2] = [0 0
]
𝑝12
0
[
−2𝑝12 −2𝑝22
]+[
−2𝑝12 𝑝11 − 3𝑝12
]− [ 2 ] [𝑝11 𝑝12 ] + [2 0] = [0 0]
𝑝11 − 3𝑝12 𝑝12 − 3𝑝22 −2𝑝22 𝑝12 − 3𝑝22 𝑝 𝑝12 𝑝22 0 2 0 0
0 22
2
2
𝑝12
−2𝑝12 − 2𝑝12 − +2=0 (1)
2
2
𝑝22
𝑝12 − 3𝑝22 + 𝑝12 − 3𝑝22 − +2=0 (3)
2
2
𝑝12
From (1), −2𝑝12 − 2𝑝12 − +2=0
2
2
𝑝12
+ 4𝑝12 − 2 = 0
2
Solving the above equation yields, 𝑝12 = −8.4721 𝑎𝑛𝑑 0.4721. Since P is a positive matrix, choose
positive value for 𝑝12
𝑖. 𝑒., 𝑝12 = 0.4721
2
𝑝22
From (3), 𝑝12 − 3𝑝22 + 𝑝12 − 3𝑝22 − +2=0
2
2
𝑝22
2𝑝12 − 6𝑝22 − +2=0
2
2
𝑝22
+ 6𝑝22 − (2 × 0.4721) + 2 = 0
2
2
𝑝22
+ 6𝑝22 − 2.9442 = 0
2
Solving the above equation yields, 𝑝22 = −12.4721 𝑎𝑛𝑑 0.4721. Since P is a positive matrix, choose
positive value for 𝑝22
𝑖. 𝑒., 𝑝22 = 0.4721
From (2)
𝑝12 𝑝22
−2𝑝22 + 𝑝11 − 3𝑝12 ++0=0
2
𝑝12 𝑝22 0.4721 × 0.4721
𝑝11 =− + 3𝑝12 + 2𝑝22 ; − + (3 × 0.4721) + (2 × 0.4721) = 2.4721
2 2
𝑝11 𝑝12 2.4721 0.4721
𝑃 = [𝑝 𝑝 ]= [ ]
12 22 0.4721 0.4721
The optimal control law 𝑢∗ = −𝐾𝑥(𝑡) = −𝑅 −1 𝐵𝑇 𝑃𝑥(𝑡)
1 2.4721 0.4721 𝑥1 (𝑡)
=− [2] [0 1] [ ][ ]
0.4721 0.4721 𝑥2 (𝑡)
2.2151 0.8284 𝑥1 (𝑡)
=− [0 1⁄2] [ ][ ]
0.8284 0.3947 𝑥2 (𝑡)
𝑥 (𝑡)
=− [0.2361 0.2361] [ 1 ]
𝑥2 (𝑡)
𝑢∗ = −0.2361𝑥1 (𝑡) − 0.2361𝑥2 (𝑡)
∴ 𝐾 = [0.2361 0.2361]
Estimation deals with finding an approximate value (called an estimate) of a quantity from
observations or measurements which contain the information on the quantity to be estimated. All
measurements are prone to errors and therefore an estimate can be only as close as possible to the correct
value depending on the methods used for the estimation. Estimation is to be done generally for the purpose
of implementing a control of a physical system.
Figure shows the general problem of estimation. The physical system is subjected to two types
of inputs, a control input which can be easily manipulated and a disturbance input which accounts for the
presence of internal or external phenomena which cannot be easily determined. The disturbance may be
inherent in the system or may be due to the environment, such as noise in electronic circuits, signal
interference due to radiation, turbulence of aircraft etc., which occur in an unpredictable or random
manner. The system variables which are the outputs of the system are then measured by a measurement
process or system. The measuring system not being the perfect one introduces its own errors which may
again be unpredictable or random. Some of the errors in the measurement may be systematic and can be
corrected for, but many estimation problems deal with the measurement errors also as random.
The estimation problem here is to correctly calculate the true values of the system variables in the
presence of the disturbances and the errors introduced due to error in measurement. Often the time history
of the measurements is available and the estimates are to be obtained from these. Usually, a performance
measure is defined to assess the quality of the approximation or the estimate and the estimate is derived
that will maximize or minimize this performance measure. In such a case, the estimate is called an optimal
estimate. What is being done in an optimal estimation problem is to determine an algorithm by which the
estimate is obtained from the measurements based on the knowledge of the dynamics of the physical system
and whatever information one has on the disturbance process and the measurement process or-the errors
caused by these. In an associated control problem, the estimates obtained are used to control the system in
a desired manner.
1. Development of the models: This involves the specification of the models for (a) the physical
1. Communication systems: To extract a message from received signal. The transmitted or received
signal will contain the message to be extracted and will be contaminated by unwanted disturbance
and measurement errors in different stages of the Transmission and reception.
2. Navigation: In a typical navigation problem the position and velocity of a vehicle have to be
correctly estimated from the available measurement on the two. Aircraft, space crafts, surface ships
and submarines use these for their movement.
3. Post experimental data analysis: Here the recoded data from an experiment is analyzed in detail
to assess the success. For example tracking and telemetry data available from the launching of a
spacecraft provides valuable information for subsequent missions.
4. Process control: Successful operation of large chemical processes require regulation and control for the
purpose of maintaining the efficiency of operation, quality of products and attainment of other specified
goals. Specific examples are machine tool control, aircraft and space craft flight control and multistage
chemical processes.
Kalman Filter
Since the time of its introduction, the Kalman filter has been the subject of extensive research and
application, particularly in the area of autonomous or assisted navigation. This is likely due in large part to
advances in digital computing that made the use of the filter practical, but also to the relative simplicity
and robust nature of the filter itself. Rarely do the conditions necessary for optimality actually exist, and
yet the filter apparently works well for many applications in spite of this situation.
Linear Quadratic Gaussian (LQG) control is a modern state space technique for designing optimal
dynamic regulators. It enables you to trade off regulation performance and control effort, and to take into
account process and measurement noise. Like pole placement, LQG design requires a state-space model of
the plant. The LQG controller design methodology is based on the Kalman filter who in 1960 published
his famous paper describing a recursive solution to the discrete-data linear filtering problem.
The LQG controller is simply the combination of a Kalman filter i.e. a linear-quadratic estimator
(LQE) with a linear-quadratic regulator (LQR). When we use the combination of an optimal estimator
(Kalman filter) and an optimal regulator (LQR) to design the controller, the compensator is called Linear
Quadratic Gaussian (LQG).
Drawbacks of LQG
LQG robustness is not guaranteed. It is thus it is important to look at the stability radius and spectral
values sets of the closed loop design. This is related to the interaction of the transients of the true state, the
controller action and the observer.