Modern Control Systems Overview
Modern Control Systems Overview
net/publication/393003170
CITATIONS READS
0 135
1 author:
SEE PROFILE
All content following this page was uploaded by Daniel Abebe Beyene on 07 August 2025.
This lecture note is a valuable resource for anyone seeking to understand the State Space
Representation of Linear Time-Invariant (LTI) Systems. It covers a wide range of topics,
including Mathematical Formulation of Control System Problems, Mathematical Model-
ing and State Space representation of Dynamic Systems, and State Space Representation
From Transfer Function.
The lecture note also provides a comprehensive review of Eigenvalues and Eigenvectors,
Solution of State Space Equations, Similarity Transformation, Phase Plane, and Stability
Analysis of LTI systems. These topics are crucial for understanding the behavior of LTI
systems.
In addition, the lecture note covers Controllability and Observability of LTI Systems,
Modality Tests, Kalman Decomposition, and Weaker Conditions for LTI Systems. It
also discusses Controllability and Observability in the Sense of Transfer Function. These
topics are essential for designing control systems that can achieve desired performance.
The lecture notes also include Full-state Feedback Controller Design for LTI Systems, Pole
Placement State Feedback Controller Design, Desired State Tracking, Servo-mechanism
Controller Design, and Determining Poles of Second Order LTI Systems From Time Do-
main Specification. These topics are crucial for designing controllers that can achieve
desired performance.
Furthermore, the lecture note covers State Estimators/Observers Design for LTI systems,
Full Order State Observer Design, Reduced Order State Observer Design, Continuous
Time Kalman Filtering, Closed Loop State Feedback Controller Design with Full State
Observer, and Closed Loop State Feedback Controller Design with Reduced Order State
Observer. These topics are essential for designing state estimators/observers that can
estimate the system’s state accurately.
The lecture note also provides an Introduction to Optimal Control, Performance Indices,
Elements of Calculus of Variation, Solution to Optimal Linear Quadratic Regulator Prob-
lem, Optimal Steady State Reference Tracking, and Optimal State Estimation. These
topics are crucial for designing optimal control systems that can achieve desired perfor-
mance.
In addition, the lecture notes cover Introduction to MIMO Systems and Nonlinear Sys-
tems, Introduction to Multiple Input Multiple Output Systems, and Introduction to
Nonlinear Systems. These topics are essential for understanding the behavior of complex
systems.
i
Finally, the lecture note includes Analog Simulation, Analog Computers and Simulation,
Solving Differential Equations with an Analog Computer, Physical Realization of Analog
Computer, Programming Differential Equation Using Analog Computer, Magnitude and
Time Scaling, Analog Simulation of Transfer Functions, and State Space representation
of Systems. These topics are crucial for understanding how analog computers can be
used to simulate complex systems.
Overall, this lecture note serves as an essential guide for anyone interested in understand-
ing the State Space Representation of linear time-invariant (LTI) Systems. It provides
a comprehensive overview of various topics related to LTI systems and is an excellent
resource for students and professionals alike.
In these lectures, the author has taken into consideration students who may not have
a strong background in state space representation and analysis for different systems.
Therefore, the author has deliberately avoided introducing basic concepts in state space
representation and state feedback control for linear systems whenever possible. While
some may view this as a drawback, the author hopes that users will not be too critical.
It is important to note that these are only introductory lectures and examples on the
topic. Due to time constraints and the availability of other books on the subject, the
author has chosen to focus on the most critical aspects and applications of the topic.
This decision was made to make the discussion more interesting from both an analytic
point of view.
The author hopes that this brief introduction will encourage interested students to ex-
plore more about this fascinating mathematical tool for analyzing and designing state
feedback controllers for different systems using state space representation. It is impor-
tant to note that there are many other books available on this topic for further reading,
such as ”Modern Control Engineering” by Katsuhiko Ogata [1] and ”Modern Control
Systems” by Richard C. Dorf and Robert H. Bishop [2].
ii
Acronym’s (Abbreviation)
Abbreviation Description
DC Direct Current
BIBO Bounded Input Bounded Output
FBD Free Body Diagram
2DOF 2 Degree of Freedom
IPC Inverted Pendulum on Cart
LQR Linear Quadratic Regulator
LTI Linear Time Invariant
MATLAB Matrix Laboratory
MIMO Multiple Inputs and Multiple Outputs
SIMO Single Input and Multiple Outputs
iii
Contents
Preface i
iv
1.15.1 Stability analysis of LTI systems . . . . . . . . . . . . . . . . . . 52
1.15.2 Phase plane analysis of LTI systems . . . . . . . . . . . . . . . . . 53
1.16 Examples . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
1.17 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72
v
5.5 Optimal State Estimation . . . . . . . . . . . . . . . . . . . . . . . . . . 153
5.6 Examples . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 155
5.7 Exercises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 164
vi
A.7 Euler-Lagrangian Equation for systems with Generalized Forces . . . . . vii
A.8 Example:-1 Inverted Pendulum on a Cart System Dynamics . . . . . . . viii
A.8.1 System Description . . . . . . . . . . . . . . . . . . . . . . . . . . viii
A.8.2 Dynamic Model of the System . . . . . . . . . . . . . . . . . . . . viii
A.9 Example:-2 Ball and Beam System . . . . . . . . . . . . . . . . . . . . . xi
A.9.1 System Description . . . . . . . . . . . . . . . . . . . . . . . . . . xi
A.9.2 Dynamic Model of the System . . . . . . . . . . . . . . . . . . . . xi
vii
Chapter 1
• Empirical Models: These models are constructed based on observed data rather
than derived from first principles. They often employ statistical methods to fit
data, making them effective for predicting system behavior within the range of
observed conditions.
• Hybrid Models: These models incorporate both continuous and discrete elements,
making them suitable for systems that exhibit both behaviors.
1
Daniel Abebe Modern Control Systems Lecture Note
on the internal state of the system and employs a set of first-order differential equations.
While you have previously covered analysis and controller design using transfer function
models in classical control systems, this course will specifically emphasize analysis and
controller design within the state space framework.
A system is a thing having parts that can be considered as a single entity. The system
must be defined clearly with a particular set of relationships that completely characterize
the conceivable states of the system.
The state of the system describes the condition(status) of the system at any time t ≥ 0
and provides information which (together with the knowledge of the equation describing
the system) enables us to calculate/analyze the future behavior of the system under
certain control input.
In practice, it is not possible to get all state variables; instead, we get controlled variables
(outputs) which depend on someway on measured states.
Equation 1.3.1 is called the state equation of a dynamic system. The function F (x(t), u(t), t)
can be linear or nonlinear, it can be separable or not(i.e it may or may not written as
sum of input mapping function and state mapping function explicitly as F (x(t), u(t), t) =
f (x(t), t) + g(u(t), t)).
Since all states of the system may not be accessible directly, we define the output (acces-
sible states) of a dynamic system as
Equation 1.3.2 is called the output equation of a dynamic system. y(t) is an output
(accessible states/controlled states) of a dynamic system, and h(x(t), u(t), t) is an output
mapping function. The function h(x(t), u(t), t) defines the relation between the output
Chapter 1 2
Daniel Abebe Modern Control Systems Lecture Note
of the dynamic system, the state of the system, and the control input.
Thus the control system problem is finding admissible (inputs that we can apply to the
dynamic system) control input u(t) that derive the evolution of the dynamic system in
the desired way or toward the desired state, by using the relation between the state
dynamics and control input as well as the relation between the controlled states (output)
of the system and the control input.
In this course, we only deal with linear time-invariant (LTI) dynamic systems. LTI
dynamic systems have linear time-invariant state mapping and input mapping functions.
This implies that the evolution of LTI dynamic systems is linear concerning the state of
the system and control input. Which simplifies the state equation 1.3.1
for LTI dynamic systems. Matrices C, D are proportionality constant matrices. Matrix
C is called an output matrix, and matrix D is called an input transmission matrix. Their
size is defined based on the number of states, outputs, and inputs a dynamic system
has. For a dynamic system with n states, p number of outputs, and m number of control
inputs, matrix C is a p × n matrix, while matrix D is a p × m matrix.
Any LTI dynamic system can be represented by a state equation and output equation,
which is called the state-space representation of LTI dynamic systems.
For simplicity of expressions and convenience in the rest of the document, any function
of time f (t) is written as f , unless expressed explicitly. With this, the state space
representation of LTI dynamic systems is written as:
ẋ =Ax + Bu (1.3.7)
y =Cx + Du (1.3.8)
Chapter 1 3
Daniel Abebe Modern Control Systems Lecture Note
D
u x + y
+
B ʃ C +
+
A
Figure 1.1: State Space Representation of a Dynamic System in Block Diagram
vR = iR R (1.4.1)
• Kirchhoff ’s Current Law (KCL): States that the total current entering a junc-
tion equals the total current leaving the junction.
• Kirchhoff ’s Voltage Law (KVL): States that the sum of the electrical potential
differences (voltage) around any closed network is zero..
diL
vL = L (1.4.2)
dt
where L is the inductance of an inductor.
dvC
iC = C (1.4.3)
dt
where C is the capacitance of a capacitor.
Chapter 1 4
Daniel Abebe Modern Control Systems Lecture Note
1.4.2 Example
RLC Circuit
Consider the series RLC circuit in Figure 1.2, if the concern of the designer is the voltage
across the capacitor, From Kirchhoff’s current law
R L
vs C vo
d
vs = iR + L i + vc
dt
d R 1 1
i = − i − vc + vs (1.4.4)
dt L L L
d
C vc = i
dt
d 1
vc = i (1.4.5)
dt C
The state space representation will be
d R 1 1
dt i − L − L i L
d = 1 + vs (1.4.6)
vc 0 vc 0
dt C
h i i
y= 0 1 (1.4.7)
vc
Chapter 1 5
Daniel Abebe Modern Control Systems Lecture Note
Φ = Bm · A (1.5.1)
where Bm is measured in teslas (T) and A in square meters (m2 ). The direction of the
magnetic flux is determined by the orientation of the magnetic field lines.
Fmmf = N · i (1.5.2)
where N is the number of turns in the coil and I is the current flowing through the coil,
measured in amperes (A).
Fmmf = Rm · Φ (1.5.4)
1.5.4 Example
Hinged Armature Type Relay
The hinged armature type relay in Figure 1.3 operates as an electromagnetic circuit
consisting of a coil (electromagnet), an armature (movable part), and a core. The strength
of this magnetic field is proportional to the current flowing through the coil (I) and the
number of turns in the coil (N ):
N ·i
Bm = µ ·
l
where Bm is the magnetic flux density, µ is the permeability of the core material, and l
is the length of the magnetic path.
Chapter 1 6
Daniel Abebe Modern Control Systems Lecture Note
The force (Fem ) exerted on the armature can be calculated using the magnetic flux and
the area of the armature:
B2 · A
Fem =
2µ0
where A is the area of the armature and µ0 is the permeability of free space.
The electromagnetic force is balanced by the spring force. When voltage is applied to
the coil, it generates a magnetic field that attracts the pivoted armature, which is held
away from the core by a spring. As the armature moves towards the core, it closes the
relay contacts, allowing current to flow through the connected circuit and activating the
device. When the current is removed, the magnetic field collapses, and the armature
returns to its original position by the spring, reopening the contacts and interrupting the
current flow.
Chapter 1 7
Daniel Abebe Modern Control Systems Lecture Note
F = −kx (1.6.1)
τ = −kθ (1.6.2)
where F and τ are the restoring force and restoring torque, respectively, k is the spring
constant (a measure of stiffness), and x and θ are the linear and angular displacements,
respectively. This relationship indicates that the force exerted by the spring opposes the
x
k
k
F
(a) (b)
If both ends are free to move or rotate, the equations describing the motion are:
F = −k(x1 − x2 ) (1.6.3)
τ = −k(θ1 − θ2 ) (1.6.4)
These equations account for the relative displacements between two points, which is
1
1
k
x1 k x2
F F
(a) (b)
important in systems where springs may be subjected to forces from multiple directions.
Chapter 1 8
Daniel Abebe Modern Control Systems Lecture Note
Fd = −dẋ (1.6.5)
τd = −dθ̇ (1.6.6)
where Fd and τd are the damping force and damping torque, respectively, b is the damping
coefficient, and ẋ and θ̇ are the linear and angular velocities. This relationship indicates
that the damping force opposes the motion, providing a stabilizing effect in dynamic
systems. Dampers are essential for controlling vibrations and ensuring system stability.
If both ends are free to move or rotate, the equations describing the motion are:
x1 x2 1 2
b b
(a) (b)
These equations illustrate how dampers respond to relative velocities, which is crucial in
applications such as automotive shock absorbers, where they mitigate the effects of road
irregularities on vehicle stability and comfort.
Mechanical Transformers
Mechanical transformers, such as levers and gears, are devices that change the direction
or magnitude of forces and motion. These components are fundamental in mechanical
systems, allowing for the efficient transfer and modification of energy.
• Levers: A lever, as shown in Figure 1.7, amplifies an input force to lift a load. The
relationship between the input force (F1 ), output force (F2 ), input distance (l1 ),
and output distance (l2 ) can be described by the law of the lever:
F1 l1 = F2 l2 (1.6.9)
This equation indicates that the moments about the fulcrum are equal, allowing
for force amplification depending on the distances from the fulcrum. The longer the
input arm (l1 ), the smaller the input force needed to lift a given load. This principle
is utilized in various applications, such as crowbars and seesaws, demonstrating the
mechanical advantage achieved through lever systems.
Chapter 1 9
Daniel Abebe Modern Control Systems Lecture Note
l1 F1
F2 l2
• Gears: Gears, as shown in Figure 1.8, transmit torque and change the speed of
rotation. The relationship between the input and output angular velocities (ω1 , ω2 )
and torques (τ1 , τ2 ) can be expressed as:
ω1 r2
= and τ1 ω1 = τ2 ω2 (1.6.10)
ω2 r1
where r1 and r2 are the number of teeth on the input and output gears (gear radius
ratio), respectively.
r1
2 2
1 1
r2
This relationship shows how gears can change speed and torque, allowing for me-
chanical advantage in various applications. For instance, in a bicycle, the gear
system enables different pedaling speeds and forces, optimizing performance based
on terrain.
Chapter 1 10
Daniel Abebe Modern Control Systems Lecture Note
1.6.3 Examples
Two Degree of Freedom Quarter Car Suspension System
A vehicle suspension system is a complex vibration system having multiple degrees of
freedom. The purpose of the suspension system is to isolate the vehicle body from the
road imperfections. Suspension systems serve a dual purpose
• Contributing to the vehicle’s road holding/handling and braking for good active
safety and driving pleasure, and
• Keeping vehicle occupants comfortable and reasonably well isolated from road noise,
bumps, and vibrations.
• Semi-active suspension systems include devices such as air springs and switchable
shock absorbers, various self-leveling solutions, as well as systems like hydro pneu-
matic, hydro elastic, and hydro-gas suspensions or
Consider a fully active suspension system shown in Figure 1.9, the free body diagram of
the system is shown in Figure 1.10.
y
mc
y1
k2 d
mw
y2
k1
z
x
Figure 1.9: 2-DOF Quarter Car Suspension Model
Using Newton’s second law of motion, the force on the body of the car is
ΣFy = mc ÿ1
fk2 + fd + u = mc ÿ1
k2 (y2 − y1 ) + d(ẏ2 − ẏ1 ) + u = mc ÿ1
k2 k2 d d 1
ÿ1 = − y1 + y2 − ẏ1 + ẏ2 + u (1.6.11)
mc mc mc mc mc
Chapter 1 11
Daniel Abebe Modern Control Systems Lecture Note
fk2 u fd
mc mw
y1 y2
fk2 u fd fk1
ΣFy = mw ÿ2
− fk2 − fd − u + fk1 = mw ÿ2
− k2 (y2 − y1 ) − d(ẏ2 − ẏ1 ) − u + k1 (z − y2 ) = mw ÿ2
k2 k1 + k2 d d k1 1
ÿ2 = y1 − y2 + ẏ1 − ẏ2 + z− u (1.6.12)
mw mw mw mw mw mw
Let y1 = x1 , ẏ1 = x2 , y2 = x3 and ẏ2 = x4 , the equations 1.6.11 and 1.6.12 simplify
ẋ1 =x2
k2 k2 d d 1
ẋ2 = − x1 + x3 − x2 + x4 + u
mc mc mc mc mc
ẋ3 =x4
k2 k1 + k2 d d k1 1
ẋ4 = x1 − x3 + x2 − x4 + z− u (1.6.13)
mw mw mw mw mw mw
The state space equation of the system is then
0 1 0 0 0
ẋ
1 x1 0
k d k d 1
− 2 −
2
0
ẋ2 m mc mc mc x 2
mc u +
= c +
z
(1.6.14)
ẋ3 0 0 0 1
x3 0 0
k2 d k1 + k2 d
1 k1
ẋ4 − − x4 −
mw mw mw mw mw mw
If the output h is the position of the body of the car concerning the road x1 , the output
equation is
x
1
h i
x2
h= 1 0 0 0
(1.6.15)
x3
x4
The goal of the designer is to design a controller that keeps the car body in a predefined
range, irrespective of the roughness of the road.
Chapter 1 12
Daniel Abebe Modern Control Systems Lecture Note
The relationship between the angular speeds and torques of the turbine and generator is
essential for assessing the efficiency and overall performance of the system.
In a typical configuration, the turbine’s angular speed (θ̇1 ) is associated with the aero-
dynamic torque (τaero ) generated by the wind acting on the blades. This torque can be
expressed as a function of wind speed, blade characteristics, and the turbine’s design
parameters. The turbine is generally coupled to the generator through a gearbox, which
adjusts the speed and torque to enhance generator performance.
The generator, in turn, possesses its angular speed (θ̇2 ) and produces a generator torque
r2
(τg ). The gearbox ratio (N = ) is critical in relating these parameters, where r1 and
r1
r2 are the radii of the turbine and generator gears, respectively.
1
wind
aero r1
b
2
g
fc
r2
The mathematical model of the system can be divided into two main components: the
turbine system and the generator system. The mathematical model for the turbine system
is given by:
where fc represents the contact force between the gears, b1 is a damping coefficient, and
J1 is the moment of inertia of the turbine system.
Similarly, the mathematical model for the generator system is expressed as:
Chapter 1 13
Daniel Abebe Modern Control Systems Lecture Note
To solve for fc from equation 1.6.17, we can substitute fc into equation 1.6.16:
1
τaero − b1 θ̇1 − τg + b2 θ̇2 + J2 θ̈2 = J1 θ̈1 (1.6.18)
N
Using the relationship between the angular speeds, we have θ̇1 = N θ̇2 , where N is the
gearbox ratio. The final expression for τaero in terms of θ̇2 is:
b2 1 J2
τaero = b1 N + θ̇2 + τg + + J1 N θ̈2
N N N
b2 + b1 N 2
N 1
θ̈2 = τaero − τg − θ̇2 (1.6.19)
J2 + J1 N 2 J2 + J1 N 2 J2 + J1 N 2
b2 + b1 N 2
N 1
Let γ1 = , γ2 = , γ3 = and the states of system
J2 + J1 N 2 J2 + J1 N 2 J2 + J1 N 2
be x1 = θ2 and x2 = θ̇2 the position and velocity of the mass of the spool, then the sate
space representation will be
ẋ1 = x2
ẋ2 = −γ3 x2 + γ1 τaero − γ2 τg (1.6.20)
ẋ1 0 1 x1 0 0 τaero
= + (1.6.21)
ẋ2 0 −γ3 x2 γ1 −γ2 τg
h i x1
y= 0 1 (1.6.22)
x2
Chapter 1 14
Daniel Abebe Modern Control Systems Lecture Note
lsinθ lsinθ
y Jα
lcosθ
θ
θ
lcosθ
P N
mpg
N
P
x
u
x ff
The free-body diagram of the IPC is shown in Figure 1.12b. If we denote the coefficient
of friction between the cart and the rail by b, the normal force exerted by the cart on the
pendulum by P , and the horizontal force exerted by the cart on the pendulum by N .
Using Newton’s second law of motion on a pendulum free body diagram, we have,
ΣFpx = mp ẍp
N = mp ẍp = mp (ẍ + θ̈l cos θ − θ̇2 l sin θ) (1.6.25)
ΣFpy = mp ÿp
P = mp g − mp θ̈l sin θ − mp lθ̇2 cosθ (1.6.26)
Using Newton’s second law of motion forces exerted on a cart are
ΣFcx = mc ẍ
u = mc ẍ + bẋ + N (1.6.27)
ΣFcy = mc ÿ = 0 (1.6.28)
Substituting for N in equation 1.6.27 we have,
u = (mc + mp )ẍ + bẋ + mp lθ̈ cos θ − mp lθ̇2 sin θ (1.6.29)
Assuming the friction between the pendulum and the cart is negligible, the momentum
of the pendulum rod about the pivot is described by,
Jα = J θ̈ = P l sin θ − N l cos θ (1.6.30)
Combining equation 1.6.25 , equation 1.6.26 and equation 1.6.30 we have
(J + mp l2 )θ̈ − mp lg sin θ = −mp lẍ cos θ (1.6.31)
Solving for ẍ and θ̈ in-terms of ẋ, x, θ̇, θ, we get the two governing dynamic equations for
IPC ,
(mp l2 + I)mp lθ̇2 sin θ − (mp l)2 g sin θ cos θ − b(mp l2 + J)ẋ + (mp l2 + I)u
ẍ = (1.6.32)
(mp + mc )(mp l2 + I) − (mp l cos θ)2
−(mp l)2 θ̇2 g sin θ cos θ + bẋmp l cos θ + (mp + mc )mp lg sin θ − mp l cos θu
θ̈ = (1.6.33)
(mp + mc )(mp l2 + I) − (mp l cos θ)2
Chapter 1 15
Daniel Abebe Modern Control Systems Lecture Note
Using small angle approximation about θ = 0 deg, we know that for small angle cos θ ≈ 1,
sin θ ≈ θ and θ̇ ≈ 0. Replacing cos θ = 1, sin θ = θ and θ̇ = 0 in equation 1.6.32 and
equation 1.6.33 we have
mp l 2 + J b(mp l2 + J) (mp l)2 g
ẍ = u− ẋ − θ (1.6.34)
q q q
bmp l mp lg(mp + mc ) mp l
θ̈ = ẋ + θ̇ − u (1.6.35)
q q q
where q = (mc + mp )(mp l2 + J) − (mp l)2 .
To represent equation 1.6.34 and equation 1.6.35 in state space let x1 = x, x2 = ẋ, x3 = θ
and x4 = θ̇. The state space representation of IPC is;
0 1 0 0 0
ẋ x1
1 2
(mp l)2 g mp l2 + J
0 − b(mp l + J) − 0 x
ẋ2 q q 2 q
= + u
ẋ3 0 0 0 1 x3 0
ẋ4 bmp l mp lg(mp + mc ) mp l
0 x4
0 −
q q q
(1.6.36)
x
1
y1 1 0 0 0 x2
y= =
(1.6.37)
y2 0 0 1 0 x3
x4
Planar Drone
The planar quadrotor is a simplified model of a quadrotor drone that operates within
a two-dimensional plane, making it a valuable tool in control theory and robotics for
studying dynamics, control strategies, and navigation algorithms. This model features a
frame with four rotors positioned at the corners, each generating lift and capable of in-
dependent control. The quadrotor possesses two degrees of freedom in position, allowing
movement along the x and y axes, and one degree of freedom in orientation, enabling
rotation about the vertical axis (yaw).
Control inputs consist of thrust produced by each rotor, with differential thrust utilized
to manage rotation and directional movement. The dynamics of the planar quadrotor
can be described by equations derived from Newton’s laws or Lagrangian mechanics (see
appendix A), incorporating forces from thrust and moments from rotor speeds. These
equations govern its position, velocity, orientation, and angular velocity, accounting for
factors such as mass, acceleration due to gravity, and moment of inertia.
The planar quadrotor finds applications in various fields, including aerial robotics re-
search, flight dynamics simulation, and the development of control algorithms. Its sim-
plicity allows for easier analysis and the creation of control strategies that can be adapted
Chapter 1 16
Daniel Abebe Modern Control Systems Lecture Note
for more complex three-dimensional systems, making it an essential model for understand-
ing the principles of flight control in aerial vehicles.
The equations of motion are almost trivial, since it is only a single rigid body, and cer-
tainly fit into our standard manipulator equations:
mẍ = −mg sin θ (1.6.38)
mÿ = (u1 + u2 ) − mg cos θ (1.6.39)
I θ̈ = r(u1 − u2 ) (1.6.40)
To express the given equations in state-space representation let x1 = x, x2 = ẋ, x3 =
y, x4 = ẏ, x5 = θ, x6 = θ̇ and let F = u1 + u2 the lift force and τ = u1 − u2 the torque,
that cause the drone to pitch.
ẋ1 = x2
ẋ2 = −g sin x5
ẋ3 = x4
F
ẋ4 = − g cos x5
m
ẋ5 = x6
r
ẋ6 = τ (1.6.41)
I
Linearizing the given state-space representation using the small angle approximation, we
assume that the angle θ or x5 is small enough such that: cos x5 ≈ 1 and sin x5 ≈ x5 .
Using these approximations, we can rewrite the equations as follows:
ẋ1 = x2
ẋ2 = −gx5
ẋ3 = x4
F
ẋ4 = −g
m
ẋ5 = x6
r
ẋ6 = τ (1.6.42)
I
Chapter 1 17
Daniel Abebe Modern Control Systems Lecture Note
I,m
x
Let I1 , m, and I2 denote the inertia of the beam, the mass of the ball, and the inertia
of the ball, respectively. The position and velocity of the center of mass of the ball are
Chapter 1 18
Daniel Abebe Modern Control Systems Lecture Note
2
1 1 ẋ
Krot = I1 θ̇2 + I2 (1.6.48)
2 2 r
The total kinetic energy of the system is
1 2 2 1 I2
K = Ktrans + Krot = (I1 + mx )θ̇ + m 1 + ẋ2 (1.6.49)
2 2 mr2
I2
Let λ = 1 + is a constant.
mr2
The potential energy of the system is
V (x, θ) = mgx sin θ (1.6.50)
The Lagrangian of the system is
L = K(x, θ) − V (x, θ)
1 2 2 1 I2
= (I1 + mx )θ̇ + m 1 + ẋ2 − mgx sin θ (1.6.51)
2 2 mr2
The only variable that can be actuated is the angle of the beam; the position of the ball
is an unactuated variable. Thus the Euler-Lagrangian equation (see appendix A) of the
systems is
d ∂L ∂L
− = Qi (1.6.52)
dt ∂ q̇ ∂q
where Qi represents the generalized forces or torques, including friction and external
inputs. Assuming that the pivot friction is negligible, the only externalforce
applied to
τ
the system is the torque at the pivot joint τ . Therefore, we have Qi = . By solving
0
the Euler-Lagrange equation, we obtain:
(I1 + mx2 )θ̈ + 2mxẋθ̇ + mgx cos θ = τ (1.6.53)
λẍ − xθ̇2 + g sin θ = 0 (1.6.54)
With state variables x1 = x, x2 = ẋ, x3 = θ, x4 = θ̇, the dynamics of the ball and beam
system can be described by the following set of nonlinear differential equations:
ẋ1 = x2
mx24 x1 − g sin x3
ẋ2 =
λ
ẋ3 = x4
τ − 2mx1 x2 x3 − mgx1 cos x3
ẋ4 = (1.6.55)
I + mx21
Chapter 1 19
Daniel Abebe Modern Control Systems Lecture Note
Using small angle approximation(x3 = 0) and the ball is at the center of the beam
(x1 = 0), which means (sin x3 ≈ x3 and cos x3 ≈ 1), the linearized state-space equations
become:
ẋ1 = x2
−gx3
ẋ2 =
λ
ẋ3 = x4
τ − mgx1
ẋ4 = (1.6.56)
I
The state-space representation of the system is:
0 1 0 0 0
−g
0 0 0 0
ẋ = λ x + τ (1.6.57)
0
0 0 0 1
−mg 1
0 0 0
I I
The control objective is to regulate the position of the ball; thus, the output must be x1 ,
which results in the output equation
h i
y= 1 0 0 0 x (1.6.58)
The Acrobot
The Acrobot is representative of the primary challenge in underactuated robots. In
order to swing up and balance the entire system, the controller must reason about and
exploit the state-dependent coupling between the actuated degree of freedom and the
unactuated degree of freedom. It is also an important system because, as we will see, it
closely resembles one of the simplest models of a walking robot.
Figure 1.15 illustrates the model parameters used in our analysis. θ1 is the shoulder
joint angle, θ2 is the elbow (relative) joint angle, and we will use q = [θ1 , θ2 ]⊤ , x = [q, q̇]⊤ .
Chapter 1 20
Daniel Abebe Modern Control Systems Lecture Note
The zero configuration is with both links pointed directly down. The moments of iner-
tia, I1 , I2 , are taken about the pivots. The task is to stabilize the unstable fixed point
x = [π, 0, 0, 0]⊤ .
We will derive the equations of motion for the Acrobot using the method of Lagrange.
The locations of the center of mass of each link, pc1 , pc2 , are given by the kinematics:
lc1 s1 l1 s1 + lc2 s1+2
pc1 = , pc2 = , (1.6.59)
−lc1 c1 −l1 c1 − lc2 c1+2
where s1 is shorthand for sin(θ1 ), c1+2 is shorthand for cos(θ1 + θ2 ), etc. The energy
is given by:
1
T = T1 + T2 , T1 = I1 q̇12 (1.6.60)
2
1 1
T2 = (m2 l12 + I2 + 2m2 l1 lc2 c2 )q̇12 + I2 q̇22 + (I2 + m2 l1 lc2 c2 )q̇1 q̇2 (1.6.61)
2 2
Entering these quantities into the Lagrangian yields the equations of motion:
(I1 + I2 + m2 l12 + 2m2 l1 lc2 c2 )q̈1 + (I2 + m2 l1 lc2 c2 )q̈2 − 2m2 l1 lc2 s2 q̇1 q̇2 (1.6.63)
Chapter 1 21
Daniel Abebe Modern Control Systems Lecture Note
The nonlinear equations of motion for the Acrobot are linearized around the equilib-
rium point q = [π, 0]⊤ and q̇ = [0, 0]⊤ . The linearization is performed using small-angle
approximations. Assume small deviations from the equilibrium point:
θ1 = π + ∆θ1 , θ2 = 0 + ∆θ2 , (1.6.71)
where ∆θ1 and ∆θ2 are small angles. For small angles, we use the following approxima-
tions:
sin(θ1 ) ≈ −∆θ1 , cos(θ1 ) ≈ −1, (1.6.72)
Chapter 1 22
Daniel Abebe Modern Control Systems Lecture Note
Where:
m1 glc1 + m2 gl1 − m2 glc2
a31 = ,
M11
m2 glc2 −m2 glc2 m2 glc2
a32 = , a41 = , a42 = ,
M11 M22 M22
1 1
b31 = , b32 = 0, b41 = 0, b42 = .
M11 M22
1.7.1 Example
Solenoid Actuator-Valve System
The 3/2 solenoid actuator valve, depicted in Figure 1.16a, plays a crucial role in regulating
fluid flow through three ports: one inlet and two outlets.
Spool Valve
- Fs
Fem m
Plunger
Fd
Drain
Chapter 1 23
Daniel Abebe Modern Control Systems Lecture Note
over actuators and cylinders in various pneumatic and hydraulic applications, establish-
ing the solenoid actuator valve as a vital component in automation and fluid management
systems.
Since the spool valve operates within a fluid medium, which introduces viscous friction
that influences its dynamic behavior. The free body diagram of the system, illustrated
in Figure 1.16b, elucidates the forces acting on the valve and provides insights into its
dynamics, particularly regarding the impact of fluid resistance on its operation.
Applying Newton’s second law of motion, the force acting on the spool valve can be
expressed as:
ΣF = Fem − Fs − Fd = mẍ
N i − kx − dẋ = mẍ (1.7.1)
where m is mass of the spool, x is displacement of the spool from its equilibrium position,
N is number of turns in the solenoid, i is the current flowing through the solenoid, k is
spring constant and d is damping coefficient.
Let the states of the system be x1 = x and x2 = ẋ, the position and velocity of the mass
of the spool. The state space representation will be
ẋ1 = x2
k d N
ẋ2 = − x1 − x2 + i (1.7.2)
m m m
ẋ 0 1 x 0
1 = 1 + i (1.7.3)
k d N
ẋ2 − − x2
m m m
h i x1
y= 1 0 (1.7.4)
x2
DC Motor System
Consider an armature voltage-controlled DC motor with armature resistance Ra and
armature inductance La used to drive a mechanical load with load torque τL having
rotational inertia J through a coupling gear with DC motor shaft gear to load side gear
rm
ratio at a speed ωL . Let the mechanical friction coefficient between the rotor and the
rL
load be b, as shown by the equivalent circuit in Figure 1.17. The control input to the DC
motor is armature voltage Va , and the output is the speed of the load.
The desired goal is to find Va that can derive the load at the required speed. To determine
Va , we have to know the relation between the speed of the DC motor ωm and armature
voltage Va , which, of course, needs the mathematical model of the DC motor.
From the electrical characteristics of the DC motor, the torque produced by the motor τm
is proportional to armature current ia , while the back emf Eb of the motor is proportional
Chapter 1 24
Daniel Abebe Modern Control Systems Lecture Note
Ra ia La
+ τm , ωm rm
Va M
- τL , ωL
Load
rL
τm =kt ia (1.7.5)
Eb =kb ωm (1.7.6)
where kt and kb are torque and back emf proportionality constants, respectively.
Using Kirchhoff’s current law
d
Va = La ia + Ra ia + Eb (1.7.7)
dt
and from the mechanical analysis of the load side
d
τm = J ωL + bωL + τL (1.7.8)
dt
Replacing for Eb and τm in equation 1.7.7 and equation 1.7.8 respectively and rearranging;
d Ra kb 1
ia = − ia − ωm + Va (1.7.9)
dt La La La
d kt b 1
ωL = ia − ωL − τL (1.7.10)
dt J J J
From the relation between gear radius ratio and angular speed of the gears, we know that
rL
ωm = ωL (1.7.11)
rm
replacing for ωm in equation 1.7.9, we have
d Ra kb rL 1
ia = − ia − ωL + Va (1.7.12)
dt La La rm La
d kt b 1
ωL = ia − ωL − τL (1.7.13)
dt J J J
If we consider load torque τL as an uncontrolled (disturbance) input, the state space
representation of DC motor dynamics is
Ra kb rL 1
i̇a − − i a 0
= kLa La rm
b + La Va + 1 τL (1.7.14)
t −
ω̇ − ωL 0
J J J
Chapter 1 25
Daniel Abebe Modern Control Systems Lecture Note
h i ia
y= 0 1 (1.7.15)
ωL
As we can see from equation 1.7.14 and equation 1.7.15, the dc motor dynamics has two
states, namely ia and ωL .
If the position of the DC motor is required for control purposes, like in servo motors, we
add one state θ̇ = ωm . In this case, the state space representation becomes
R kb rL
−
a
− 0 1
i̇a La i
La rm a La 0
k b 1
t
ω̇L = − 0 ωL + V + − τL (1.7.16)
0
a
J J J
rL
θ̇L 0 0 θL 0 0
rm
i
h i a
y = 0 0 1 ωL (1.7.17)
θL
+ ia
Ra La
+
− +
Rf ωm
Vs Va M Load
Ri
-
-
(a) Armature Voltage Controlled DC Mo- (b) Remote Controlled
tor With Pre-amplifier Projector Screen
control the position of a projector screen. The pre-amplifier boosts the input armature
voltage supplied to the DC motor, enhancing its performance and responsiveness. The
Chapter 1 26
Daniel Abebe Modern Control Systems Lecture Note
position of the leading edge of the projector screen is proportional to the angular position
of the coupling gear. Mathematically, this relationship is expressed as:
y = θrL (1.7.18)
where: y is the position of the leading edge of the screen, theta is the angular position
of the coupling gear, rL is the radius of the coupling gear.
The state-space model describes the dynamics of the system, including how the armature
current, angular velocity, and screen position evolve in response to the input voltage and
load torque. The output equation specifies that the position of the leading edge of the
screen y(t) is the measured output of the system.
R kb rL
−
a
− 0 Av
i̇ La La rm ia La 0
a
kt b 1
ω̇L = − + V + − τL (1.7.19)
0 ωL
0 a J
J J
rL
ẏ 0 rL 0 y 0 0
rm
i
h i a
y(t) = 0 0 1 ωL (1.7.20)
y
• Continuity Equation: This principle states that mass cannot be created or de-
stroyed in a closed system. For incompressible fluids, the continuity equation can
be expressed as:
∇·v =0 (1.8.1)
where v is the fluid velocity vector. This equation ensures that the mass flow rate
remains constant throughout a streamline.
Chapter 1 27
Daniel Abebe Modern Control Systems Lecture Note
• Bernoulli’s Equation: This equation relates the pressure, velocity, and elevation
in a moving fluid and is derived from the conservation of energy. For an incom-
pressible, non-viscous fluid, it can be expressed as:
1
P + ρv 2 + ρgh = K (1.8.2)
2
where P is the pressure, ρ is the fluid density, v is the flow velocity, g is the
acceleration due to gravity, and h is the height above a reference level.
where ρ is the fluid density, p is the pressure, µ is the dynamic viscosity, and f
represents body forces (e.g., gravity).
• Law of Conservation of Mass: This law states that the mass of a closed system
must remain constant over time. In fluid dynamics, this is often represented by the
continuity equation mentioned above.
Fluid Inertia
Fluid inertia refers to the tendency of a fluid to resist changes in its state of motion. It
is quantified by the mass of the fluid and can be represented in terms of fluid pressure
(P ) and volumetric flow rate (Q̇):
∆P
Lf = (1.8.4)
∆Q̇
where Lf is the fluid inertia. This property is crucial in analyzing transient flows and
dynamic responses in fluid systems.
Fluid Capacitance
Fluid capacitance describes the ability of a fluid system to store volume in response to
pressure changes. It is defined as the ratio of the change in volume (∆V ) to the change
in pressure (∆P ):
∆V
C= (1.8.5)
∆P
where C is the fluid capacitance. This property is important in systems like hydraulic
accumulators, where fluid storage plays a significant role in performance.
Chapter 1 28
Daniel Abebe Modern Control Systems Lecture Note
Laminar Flow
Laminar flow occurs when a fluid flows in parallel layers, with minimal disruption between
them. This type of flow is characterized by smooth, orderly motion and is typically
observed when the source pressure and outlet pressure are the same, or when the diameter
of the pipe is larger than the length of the pipe. The laminar flow resistance is proportional
to the ratio of pressure difference to volumetric flow rate:
∆P
RL = (1.8.7)
Q
Turbulent Flow
Turbulent flow is characterized by chaotic, irregular fluid motion. This type of flow
involves mixing and eddies, leading to increased energy dissipation and resistance. The
transition from laminar to turbulent flow is influenced by factors such as flow velocity,
fluid viscosity, and pipe diameter. The Turbulent flow resistance is proportional to the
ratio of the pressure difference to the square of the volumetric flow rate:
∆P
RT = (1.8.8)
Q2
1.8.4 Examples
Hydraulic Piston and Cylinder System
Hydraulic systems are essential in various industrial applications, facilitating the efficient
transfer of power through the controlled movement of fluids. Central to many hydraulic
systems is the hydraulic piston and cylinder arrangement, which converts hydraulic energy
into mechanical work. This system operates on the principles of fluid mechanics and
adheres to fundamental laws, including Pascal’s principle, which asserts that pressure
Chapter 1 29
Daniel Abebe Modern Control Systems Lecture Note
Qin CV
Using the law of conservation of mass, we can express the rate of accumulation of mass
as the difference between the mass inflow rate and the mass outflow rate. In this case,
since there is no mass leaving the system (no outflow), the dynamic equation governing
the motion of the piston can be formulated as:
ṁ = ρQin (1.8.9)
Here, ṁ represents the rate of change of mass, ρ is the fluid density, and Qin is the
volumetric flow rate into the system.
Next, we can substitute the mass m with the product of density and volume, m = ρV .
This leads to the following relationship:
In this equation, V̇ denotes the rate of change of volume, and V is the volume of the fluid
within the cylinder.
To further refine our model, we can replace the density ρ using the relationship derived
from the bulk modulus of elasticity, which is given by:
ρo dP
ρ= (1.8.11)
β dt
dP
where ρo is the reference density, β is the bulk modulus, and represents the rate of
dt
change of pressure. Additionally, we can express the rate of change of volume V̇ in terms
of the piston head area A and the displacement of the piston head x:
V̇ = Aẋ (1.8.12)
Chapter 1 30
Daniel Abebe Modern Control Systems Lecture Note
In this final equation, Ṗ represents the rate of change of pressure, A is the cross-sectional
area of the piston, and ẋ is the velocity of the piston head. This equation captures
the dynamic relationship between pressure changes, fluid density, and the movement of
the piston, allowing us to analyze the behavior of the hydraulic system under varying
conditions.
By integrating these principles, we can effectively model the hydraulic piston and cylinder
system, facilitating predictions of its performance in practical applications.
k Fatm
Fs
Qin Fh FL
FL Fv
The free body diagram of the system is illustrated in Figure 1.20b. This diagram helps
visualize the forces acting on the actuator, which include the hydraulic force, atmospheric
pressure force, viscous force, spring force, and the load force.
The mathematical model describing the motion of the actuator is derived using New-
ton’s first law of motion, which states that the sum of forces acting on an object is equal
to the mass of the object multiplied by its acceleration. This can be expressed as:
Where Fh is the hydraulic force acting on the piston, Fatm is the force due to atmospheric
pressure, Fv is the viscous force, Fs is the spring force, FL is the load force, m is the mass
of the piston, and ẍ is the acceleration of the piston.
By replacing Fh with P A (where P is the pressure and A is the piston head area), and
substituting Fatm = Patm A for the atmospheric pressure force, the equation becomes:
Chapter 1 31
Daniel Abebe Modern Control Systems Lecture Note
Selecting the states of the system x1 = x, x2 = ẋ, the sate space representation will be
ẋ1 = x2
k b PA 1 Patm A
ẋ2 = − x1 − x2 + − FL − (1.8.16)
m m m m m
ẋ1 0 1 x1 0 0 P 0
= + − Patm (1.8.17)
k b A 1 A
ẋ2 − − x2 − FL
m m m m m
h i x1
y= 1 0 (1.8.18)
x2
To find the pressure required to hold the piston at its nominal position, we can set the
velocity ẋ = 0 and assume x = x0 (the equilibrium position). This simplifies our equation
to:
FL k
P = Patm + + x0 (1.8.19)
A A
This equation indicates that the pressure P must balance the load force per unit area,
plus the atmospheric pressure acting on the piston. Understanding this relationship
is essential for designing hydraulic systems that effectively manage loads and maintain
desired positions.
Note:- TZB (through the boundary of the system), IZB (In the boundary of the system)
To formulate the dynamic equation, we have to decide the control boundary and which
balance we have to use. Once the boundary is set, if the interest is,
• Concentration→component balance,
Chapter 1 32
Daniel Abebe Modern Control Systems Lecture Note
With concentrations c̄A = nA /V , cA,in = FA,in /qin , cA,out = FA,out /qout and the generation
of component as average reaction rate times the volume GA = r̄A V , the equation will be
d(c̄A V )
= cA,in qin − cA,out qout + r̄A V [molA/s] (1.9.8)
dt
1
Volumetric flow rate is the product of flow rate and cross-sectional area perpendicular to the flow
rate
Chapter 1 33
Daniel Abebe Modern Control Systems Lecture Note
where H[J] is enthalpy of the system (H = U(internal energy) + pV ), Hin − Hout is the
sum of internal energy in the streams plus the flow work that the streams perform on the
system as they are ”pushed” in or out of the system. Q [J/s] is supplied heat (through
the system’s wall) and Ws [J/s] is supplied useful mechanical work (usually shaft work,
for example, from a compressor, pump, or turbine).
dV dp
The term −(pex − p) +V occurs if the volume and pressure in the boundary of
dt dt
control changes, which for may process it is negligible compared to the other parameters.
H
Using specific enthalpy ℏ = , we can combine mass balance and energy balance equa-
m
tions, since
dH dℏ dm
=m +ℏ (1.9.10)
dt dt dt
dm
by replacing for = win − wout
dt
dℏ dV dp
m = win (ℏin − ℏ) − wout (ℏout − ℏ) + Q + Ws − (pex − p) +V (1.9.11)
dt dt dt
1.9.2 Examples
Stirred-Tank Heating Process: Constant Holdup process
The liquid inlet stream consists of a single component with a mass flow rate wi and an
inlet temperature Ti . The tank contents are agitated and heated using an electrical heater
that provides a heating rate,Q. The modeling assumptions
• Perfect mixing; thus, the exit temperature T is also the temperature of the tank
contents.
• The inlet and outlet flow rates are equal; thus, wi = w and the liquid holdup V is
constant.
• The density ρ and heat capacity C of the liquid are assumed to be constant. Thus,
their temperature dependence is neglected.
Chapter 1 34
Daniel Abebe Modern Control Systems Lecture Note
Ti
wi
T
w
The enthalpy (total energy of the system) is the sum of internal energy and the product of
volume and pressure ℏ = Ūint + pV , but for low pressure system, it can be approximated
to ℏ ≈ Ūint , internal energy per unit mass.
The system enthalphy can be considered as a function of only temperature and pressure,
ℏ(T, p), hence
dℏ ∂ℏ dT ∂ℏ dp
= f (p) + f (T ) (1.9.12)
dt ∂T dt ∂p dt
∂ℏ ∂ℏ
For the fluid and solids ≈ 0, and f (p) = Cp (heat capacity at constant pressure),
∂p ∂T
For the total volume we have
dH dℏ dT
=m = ρV Cp (1.9.13)
dt dt dt
Solving for H by integrating we have H − Href = ρV Cp (T − Tref ), in specific enthalpy,
ℏ−ℏref = Cp (T −Tref ), for inlet stream ℏi −href = Cp (Ti −Tref ) and ℏ−href = Cp (T −Tref )
outlet stream.
Chapter 1 35
Daniel Abebe Modern Control Systems Lecture Note
We know that ℏ − ℏref = Cp (T − Tref ), without loss of generality, we can set ℏref = 0,
dℏ dT
hence ℏ = Cp (T − Tref ), replacing this values in equation 1.9.11 and = Cp ,
dt dt
dV dT
ρCp (T − Tref ) + ρV Cp = wi Cp (Ti − Tref ) − wCp (T − Tref ) + Q (1.9.18)
dt dt
dV
replacing for ρ with wi − w from equation 1.9.15 and rearranging
dt
1 1 1
Ṫ = wi Ti − T+ Q (1.9.19)
ρV Cp ρV Cp ρV Cp
1 1
V̇ = wi − w (1.9.20)
ρ ρ
Chapter 1 36
Daniel Abebe Modern Control Systems Lecture Note
temperature, not the liquid contents. Temperature gradients in the heating element are
neglected, assuming it has a uniform temperature (Te), which represents the average tem-
perature for the heating. Based on this new assumption, and the previous assumptions,
the unsteady-state energy balances for the tank and the heating element can be written
as
dT
mCp = wCp (Ti − T ) + he Ae (Te − T ) (1.9.21)
dt
dTe
me Ce = Q − he Ae (Te − T ) (1.9.22)
dt
Where m = ρV and me Ce is the product of the mass of metal in the heating element and
its specific heat. The term he Ae is the product of the heat transfer coefficient and area
available for heat transfer.
Note that mC and me Ce are the thermal capacitance of the tank contents and the heating
element, respectively. Q is an input variable, the thermal equivalent of the instantaneous
electrical power dissipation in the heating element.
To simplify the analysis, some common assumptions must be made. These assumptions
include:
• It is assumed that the contents of the reactor are perfectly mixed, meaning that
the concentration of reactants and products is uniform throughout the reactor at
any given time.
• The reactor operates under steady-state conditions, where the operating parameters
such as flow rates, temperatures, and concentrations do not change with time.
• The volume of the reactor remains constant during the reaction, and there are no
significant volume changes due to the reaction.
• The reactor is assumed to behave ideally, with no heat losses, no mass transfer
limitations, and no non-ideal mixing effects.
r = kCA (1.9.23)
Chapter 1 37
Daniel Abebe Modern Control Systems Lecture Note
Pure A
q, c Ai , Ti
Mixture of A and B
q, c A , T
Tc
where r is the reaction rate of A per unit volume, k is the reaction constant, and CA is
the concentration of A. For single-phase reactions, the rate constant is typically a strong
function of reaction temperature given by the Arrhenius relation,
−E
k = ko e RT (1.9.24)
where ko is the frequency factor and E is the activation energy, usually determined by
fitting experimental data, and R is the gas constant.
The inlet stream consists of pure component A with molar concentration, CAi . A cooling
coil is used to maintain the reaction mixture at the desired operating temperature by
removing heat that is released in the exothermic reaction.
Chapter 1 38
Daniel Abebe Modern Control Systems Lecture Note
Note:- both flow rates are denoted by the symbol q in the Figure.
For the stated assumptions, the unsteady-state component balances for species A (in
molar concentration units) are
dCA
V = q(CAi − CA ) − V kCA (1.9.26)
dt
This balance is a special case of the general component balance.
For the unsteady-state energy balance for the CSTR, we need an additional five assump-
tions:
• The thermal capacitance of the coolant and the cooling coil wall is negligible com-
pared to the thermal capacitance of the liquid in the tank.
• All of the coolant is at a uniform temperature, Tc . (That is, the increase in coolant
temperature as the coolant passes through the coil is neglected.)
• The rate of heat transfer from the reactor contents to the coolant is given by
Q = U A(Tc − T )
where U is the overall heat transfer coefficient and A is the heat transfer area. Both
of these model parameters are assumed to be constant.
• The enthalpy change associated with the mixing of the feed and the liquid in the
tank is negligible compared with the enthalpy change for the chemical reaction.
• Shaft work and heat losses to the ambient can be neglected.
With the use of the new assumptions and inclusion of energy balance, the dynamic model
of CSTR is given by
dT
V ρC = wC(Ti − T ) + (−∆HR )V kCA + U A(Tc − T ) (1.9.27)
dt
This model is nonlinear as a result of the many product terms and the exponential
temperature dependence of k; hence, it should be solved by numerical integration. The
model will be more complex if
Chapter 1 39
Daniel Abebe Modern Control Systems Lecture Note
So we have
q1in = γ1 k1 V1 q2in = γ2 k2 V2
q3in = (1 − γ2 )k2 V2 q4in = (1 − γ1 )k1 V1
From fluid dynamics, the rate of mass accumulation in each tank is equal to the difference
between mass inflow and mass outflow.
ṁ = ρ(qin − qout ) (1.9.28)
where qin is volumetric inflow rate and qout is volumetric outflow rate. The volumetric
outflow rate is determined by considering each tank as an isolated system. Using the
Bernoulli equation
1 1
P1 + gh1 + v12 = P2 + gh2 + v22 (1.9.29)
2 2
2
Volumetric flow rate q is equal to the product of the flow rate of the fluid (ν) times the area (a)
perpendicular to the flow. (q = ν × a)
Chapter 1 40
Daniel Abebe Modern Control Systems Lecture Note
If we take the top as point 2 and the bottom of the tank as point 1 (h1 = 0), both of
them are exposed to atmospheric pressure, so P1 = P2 . Similarly, since there is no inflow
from the top, v2 = 0. The fluid flows out of each tank due to gravity only.
with these assumption’s equation 1.9.29 simplifies to
p
v1 = 2gh2 (1.9.30)
Let the outlet have a cross-sectional area a, the outflow rate will be
p
qout = av = 2gh (1.9.31)
m m
from basic physics ρ = = , rate of accumulation of mass ṁ is
V Ah
dh
ṁ = ρA (1.9.32)
dt
Combining equation 1.9.28 and equation 1.9.32, the dynamic equation for each tank will
be
√ √
a3 2gh3 a1 2gh1
γ1 k1
− 0
ḣ
1
√A1 √A1 A1
a 2gh a 2gh γ2 k2
4 4 2 2
− 0
ḣ2 A A2
A2
V
1
= 2 √ +
(1.9.34)
a3 2gh3
(1 − γ2 )k2
ḣ3 −
0 V
2
A3
A
√ 3
ḣ4 a4 2gh4 (1 − γ1 )k1
− 0
A4 A4
If the levels of the liquid in tank 1 and tank 2 are required to be controlled, the output
equation will be
h
1
y1 1 0 0 0 h2
y= =
(1.9.35)
y2 0 1 0 0 h3
h4
Chapter 1 41
Daniel Abebe Modern Control Systems Lecture Note
Chapter 1 42
Daniel Abebe Modern Control Systems Lecture Note
Y (s) = bn U (s) + s−1 (bn−1 U (s) − an−1 Y (s) + s−1 (bn−2 U (s) − an−2 Y (s) + . . .
+s−1 (b1 U (s) − a1 Y (s) + s−1 (b0 U (s) − a0 Y (s)))))
(1.10.10)
Let
Solving for Y (s) from equation 1.10.12 and replace for Y (s) in equation 1.10.11 and
rearranging we have
Chapter 1 43
Daniel Abebe Modern Control Systems Lecture Note
Chapter 1 44
Daniel Abebe Modern Control Systems Lecture Note
Let
1
X1 (s) = U (s)
s + p1
1
X2 (s) = U (s)
(s + p1 )2
which can be written as
1
X2 (s) = X1 (s)
s + p1
in similar way
1 1
X3 (s) = X2 (s), . . . , Xr (s) = X(r−1) (s)
s + p1 s + p1
Since the rest poles are distinct, the derivation follows the way we drive for diagonal
normal form. In the time domain, the states will be
ẋ1 = −p1 x1 + u
ẋ2 = −p1 x2 + x1
..
.
ẋr = −p1 xr + xr−1
ẋr+1 = −p2 xr+1 + u
..
.
ẋn = −pn−r xn−r + u (1.10.21)
Chapter 1 45
Daniel Abebe Modern Control Systems Lecture Note
If the initial condition is zero(x0 = 0, which is the case for all transfer functions), the
transfer function from input to output is
Y(s)
= C(sI − A)−1 B + D
U(s)
Cadj(sI − A)B + det(sI − A)D
= (1.11.3)
det(sI − A)
Av = λv (1.12.1)
(A − λI)v = 0 (1.12.2)
Chapter 1 46
Daniel Abebe Modern Control Systems Lecture Note
Left eigenvectors describe invariant directions under the row transformation induced by
A (as opposed to column transformation for right eigenvectors).
Chapter 1 47
Daniel Abebe Modern Control Systems Lecture Note
The right eigenvectors define mode shapes, which represent specific combinations of state
variables that evolve exponentially according to their corresponding eigenvalues. Physi-
cally, these manifest as observable vibration patterns.
Left eigenvectors determine the controllability properties of the system through their inner
products ν i B; when this product is zero, the corresponding mode becomes uncontrollable.
This eigen-structure underpins key control techniques, including pole placement (eigen-
value assignment), model reduction (mode truncation, approximation of the system with
dominant poles), and fault detection (since right eigenvectors project the states to the
observable subspace).
Chapter 1 48
Daniel Abebe Modern Control Systems Lecture Note
The CeAt x(0) inZ equation 1.13.3 is called the zero-input response of the dynamic system,
t
while the term (CeA(t−τ ) B + D)u(τ )dτ is called the zero initial condition response.
t=0
The frequency domain representation of equation 1.13.3 is equation 1.11.2.
The term eAt denoted as Φ(t) is called state transition matrix as it transfer the state of a
dynamic system from any initial condition x(t0 ) to a final state x(t) whether the dynamic
system is driven by an input or not. A few properties of the state transition matrix are
stated below (It is a property due to the exponential nature of the function).
• Φ(t)n =Φ(nt)
i=n
At
X (At)i
Φ(t) = e = (1.13.5)
i=0
i!
Chapter 1 49
Daniel Abebe Modern Control Systems Lecture Note
f (A) can be represented by a finite sum of powers of the matrix not exceeding order
n − 1.
Where i = 1, 2, . . . n.
If the matrix A r multiplicity of eigenvalue, the polynomial with the multiplicity
eigenvalues is differentiated r times and used as an additional equation, to have an
equal number of equations as the number of variables.
y = CTξ + Du
y = C̄ξ + D̄u (1.14.2)
Chapter 1 50
Daniel Abebe Modern Control Systems Lecture Note
dynamic systems are complex in some coordinate systems, but simpler in other coordinate
systems.
Similarity transformation only changes the coordinates in which the dynamic system is
expressed. It does not change
det(λI − Ā) = 0
det(λI − T−1 AT) = 0
det(T−1 λT − T−1 AT) = 0
det(T−1 )det(λI − A)det(T) = 0
det(λI − A) = 0
Note Since the state of the system has changed solution to the state equation is different.
But we can derive a solution for the original state equation from the solution of the new
state equation by x(t) = Tξ(t).
Chapter 1 51
Daniel Abebe Modern Control Systems Lecture Note
Internal Stability
Internal stability deals with the stability of the states of the dynamic system. It analyzes
the behavior of the system state under no external inputs(zero-input response of the
states). Consider a dynamic system ẋ = Ax evolving from an initial condition x0 ;
the equilibrium point xeq (t) satisfies the state equation of the dynamic system. This
equilibrium point is
• Laypunov stable if and only if given ϵ > 0 there exist σ > 0 such that ||x0 || < σ
implies that ||x(t)|| < ϵ for all time t > 0. This means for every initial condition in
a sphere of radius σ, the final state of the system should be in a sphere of radius
ϵ. LTI system is Laypunov stable if and only if trace3 (A) < 0, det(A) ≥ 0 and
rank(A) = rank(A2 ) (i.e. All eigenvalues of system matrix A should have zero or
negative real parts).
• Exponentially stable if it is stable and there exist σ > 0, k > 0, λ > 0 such that
||x0 || < σ implies ||x(t)|| < ke−λt ||x0 || for all t > 0. If σ can be any real number
(anywhere), then the system is globally exponentially stable.
To determine whether the LTI system is Lyapunov stable or not, Let us take a quadratic
function V = x⊤ Px, where P are symmetric4 positive definite5 matrix, then
V̇ = x⊤ Pẋ + ẋ⊤ Px
= x⊤ P(Ax) + (Ax)⊤ Px
= x⊤ (PA + A⊤ P)x (1.15.1)
3
trace(A) is the sum of all diagonal elements of A
4
A matrix P is called symmetric if and only if P = P⊤
5
A matrix P is positive definite if and only if all eigenvalues of P are positive or leading principal
minors of the are positive
Chapter 1 52
Daniel Abebe Modern Control Systems Lecture Note
PA + A⊤ P + R = 0 (1.15.2)
BIBO/External Stability
It is a stability analysis of the response of the dynamic system to bounded input.
Definition:- A dynamic system is bounded-input input bounded-output, and stable if
and only if for zero initial condition, if bounded input always invokes bounded output.
Consider an LTI dynamic system with zero initial condition response
Z t
y(t) = C eA(t−τ ) Bu(τ )dτ + Du(t) (1.15.3)
0
Which implies that forZ an LTI system to be BIBO stable the weighting matrix pattern
t
must be bounded (ı.e CeA(t−τ ) Bdτ < ∞).
0
For the LTI system, asymptotic stability always implies bounded-input input bounded-
output stability. While minimal bounded input and bounded output imply asymptotic
stability.
ẋ1 = f1 (x1 , x2 )
ẋ2 = f1 (x1 , x2 ) (1.15.5)
where x1 and x2 are states of the system and f1 and f2 are linear (for LTI systems) and
nonlinear for nonlinear autonomous systems.
Since there is no control input, the phase plane analysis is used to understand the prop-
erties of the system under various initial conditions. For linear time invariant systems
equation 1.15.5 simplifies to
Geometrically, the state space of this system is a plane having x1 , x2 as coordinates. This
plane is called a phase plane. The solution of equation 1.15.5 or equation 1.15.6 with
time varying from zero to infinity can be represented as a curve in the phase plane. Such
Chapter 1 53
Daniel Abebe Modern Control Systems Lecture Note
a curve is called a phase plane trajectory. A family of phase plane trajectories is called
a phase portrait of a system.
Equation 1.15.6 can be written as
ẍ1 =a11 ẋ1 + a12 ẋ2 = a11 ẋ1 + a12 (a21 x1 + a22 x2 )
ẋ1 − a11 x1
ẍ1 =a11 ẋ1 + a12 (a21 x1 + a22 )
a12
(a11 a12 + a12 a22 ) (a12 a21 − a12 a22 a11 )
ẍ1 = ẋ1 + x1
a12 a12
ẍ1 +ξ1 ẋ1 + ξ2 x1 = 0 (1.15.7)
(a11 a12 + a12 a22 ) (a12 a21 − a12 a22 a11 )
where ξ1 = and ξ2 = .
a12 a12
The solution of equation 1.15.7 has the form
x(t) = k1 eλ1 t + k2 eλ2 t f or λ1 ̸= λ2 (1.15.8)
λ1 t λ2 t
x(t) = k1 e + tk2 e f or λ1 = λ2 (1.15.9)
A point in the phase plane where ẋ = 0 is called a singular point or equilibrium point of
the system. For a Linear system, we have only one singular/equilibrium point, namely
the origin. The trajectories in the vicinity of this singularity point can display quite
different characteristics, depending on the values of λ1 and λ2 . The following cases can
occur
• λ1 ,λ2 are both real and have the same sign (+ or -)
• λ1 ,λ2 are both real and have opposite sign
• λ1 ,λ2 are complex conjugates with non-zero real parts
• λ1 ,λ2 are complex conjugates with real parts equal to 0
The characteristics of the system under these cases are described below
1. The first condition corresponds to a node
(a) λ1 < 0,λ2 < 0 singularity point is called stable node. In this, every phase
portrait goes to the singularity point as shown in Figure 1.24, which means
the system is asymptotically stable.
(b) λ1 > 0,λ2 > 0 singularity point is called unstable node. In this, every phase
portrait goes away from the singularity point as shown in Figure 1.25, which
means the system is unstable.
Chapter 1 54
Daniel Abebe Modern Control Systems Lecture Note
(a) Re(λ1 < 0,λ2 < 0): stable focus, where trajectories circulate the point and
finally tend to the point as shown in Figure 1.27. The circulation shows that
there is an oscillation in the system response, and the convergence shows that
these oscillations will die out and the system becomes stable over time.
(b) Re(λ1 > 0,λ2 > 0): unstable focus, where trajectories circulate the point and
diverge to infinity as shown in Figure 1.28. The circulation shows that there
is an oscillation in the system response, and the divergence shows that these
oscillations grow large and the system becomes unstable over time.
4. The last condition corresponds to a center point. All trajectories are ellipses, and
the singularity point is the center of these ellipses as shown in Figure 1.29.
The stability characteristics of linear systems are uniquely determined by the nature of
their singularity points.
Chapter 1 55
Daniel Abebe Modern Control Systems Lecture Note
1.16 Examples
1. Find the controllable canonical form, observable canonical form, and diagonal nor-
mal form or Jordan canonical form of the following transfer function.
Y (s) 2s3 + 5s2 + 3s + 1
= 4
U (s) 4s + 2s2 + 3s + 6
0.5s3 + 1.25s2 + 0.75s + 0.25
=
s4 + 0.5s2 + 0.75s + 1.5
Solution:
ẋ 0 1 0 0 x1 0
1
ẋ2 0 0 1 0 x2 0
= + u
ẋ3 0 0 0 1 x3 0
ẋ4 −1.5 −0.75 −0.5 0 x4 1
x
1
h i
x2
y = 0.25 0.75 1.25 0.5
x3
x4
(b) Observable canonical form can be determined from controllable canonical form
by Ao = A⊤ ⊤ ⊤
c , Bo = Cc , Co = Bc and Do = Dc .
⊤
ẋ 0 0 0 −1.5 x 0.25
1 1
ẋ2 1 0 0 −0.75 x2 0.75
= + u
ẋ3 0 1 0 −0.5 x3 1.25
ẋ4 0 0 1 0 x4 0.5
x
1
h i
x2
y= 0 0 0 1
x3
x4
(c) Since the poles of the system are distinct, p1 = 0.71 + i,p2 = 0.71 − i, p3 =
−0.71 + 0.70i and p4 = −0.71 − 0.70i. We have a diagonal normal form. The
Chapter 1 56
Daniel Abebe Modern Control Systems Lecture Note
2. Find the state transition matrix for the following dynamic system
0 1
ẋ = x
−2 −3
Solution
−1
s 0 0 1
Φ(t) = L−1 ((sI − A)−1 ) = L−1 −
0 s −2 −3
−1
s −1
= L−1
2 s+3
s −1
adj
2 s+3
−1
=L
det s −1
2 s+3
Chapter 1 57
Daniel Abebe Modern Control Systems Lecture Note
s+3 1
(s + 2)(s + 1) (s + 2)(s + 1)
= L−1 −2 s
(s + 2)(s + 1) (s + 2)(s + 1)
−1 2 −1 1
+ 2 + s + 1 s + 2 + s + 1
= L−1 s−2 −2 2 −1
+ +
s + 2 s + 1 s +2 s + 1
2e−t − e−2t e−t − e−2t
Φ(t) =
2e−2t − 2e−t 2e−2t − e−t
2
i=2
X (At) i 1 0 0 −1 0 −1
Φ(t) = = + t + t2
i=0
i! 0 1 −2 −3 −2 −3
1 0 0 −1 −2 −3 2
= + t + 1 t
0 1 −2 −3 2 6 7
2 3 2
1−t −t − t
2
= 7
−2t + 3t 1 − 3t + t2
2
2
i=2
at
X (at)i
From Taylor expansion of e = ,
i=0
i!
1 4
e−t = 1 − t + t2 , e−2t = 1 − 2t + t2
2 2
−t −2t −3 2
2
1 − t = 2e − e , −t + t = e−t − e−2t
2
−2t −t 7
2
− 2t + 3t = 2e − 2e , 1 − 3t + t2 = 2e−2t − e−t
2
So,
−3 2
2
1−t −t + t
−t −2t −t −2t
2 2e − e e − e
Φ(t) =
=
7 2 −2t −t −2t −t
−2t + 3t2 1 − 3t + t 2e − 2e 2e − e
2
(c) Method 3:-Using Cayley Hamilton Theorem
Determining the eigenvalues
det (λI − A) = (λ + 2)(λ + 1) = 0
Which result in eigenvalues of matrix A, λ1 = −2 and λ2 = −1. Hence the
simultaneous equations from eλi t = α0 + α1 λi are
e−2t = α0 − 2α1
e−t = α0 − α1
Chapter 1 58
Daniel Abebe Modern Control Systems Lecture Note
Φ(t) = α0 I + α1 A
1 0 0 −1
= (2e−t − e−2t ) + (e−t − e−2t )
0 1 −2 −3
−t −2t −t −2t
2e − e e −e
Φ(t) =
−2t −t −2t −t
2e − 2e 2e − e
Hence the modality M matrix, it’s inverse M−1 and eΛt are
−2t
1 1 −1 1 e 0
M= , M−1 = , eΛt = ,
2 1 2 −1 0 e−t
Chapter 1 59
Daniel Abebe Modern Control Systems Lecture Note
Finally
−2t
−1 1 e 0 1 1
Φ(t) = M−1 eΛt M =
−t
2 −1 0 e 2 1
2e−t − e−2t e−t − e−2t
=
−2t −t −2t −t
2e − 2e 2e − e
Solution
2
s − s + 3 −1 2−s 0
h i
0 1 1 −2s − 2 s2 − 1 7 − 5s 1
2
−2 s − 1 s − 3s + 2 0
= +2
s3 − 2s2 + 4s − 5
s2 + s − 2
= +2
s3 − 2s2 + 4s − 5
Chapter 1 60
Daniel Abebe Modern Control Systems Lecture Note
(λ1 I − A)ν 1 = 0
−1 0 1 −1 ν
− 11 = 0
0 −1 0 −1 ν12
−2 1 ν
11 = 0
0 0 ν12
Which results in
Hence ν 1
ν11 1
ν1 = = m
ν12 2
(λ2 I − A)ν 2 = 0
1 0 1 −1 ν
− 21 = 0
0 1 0 −1 ν22
0 1 ν
21 = 0
0 2 ν22
Which results in
ν22 = 0, 2ν22 = 0
5. Determine the eigenvalues, right eigenvector, and left eigenvectors of the following
matrix
1 2
A=
3 0
Chapter 1 61
Daniel Abebe Modern Control Systems Lecture Note
Solution
λ 0
1 2 λ − 1 −2
det( − = det = 0
0 λ 3 0 −3 λ
λ(λ − 1) + 5 = 0
λ2 − λ + 6 = 0 ⇒ (λ − 3)(λ + 2) = 0
λ1 = 3 λ2 = −2
the eigenvectors associated with λ1 = 3, ν 1 is
(λI − A) ν 1 = 0
3 0 1 2 ν11 2 −2 ν11
− = 0 ⇒ = 0
0 3 3 0 ν12 −3 3 ν12
which results in
2ν11 − 2ν12 = 0, −3ν11 + 3ν12 = 0
ν11 1
ν11 = ν12 = 1 ν 1 ⇒ = ν11
ν11 1
−2 0 1 2 ν11 −3 −2 ν21
− = 0 ⇒ = 0
0 −2 3 0 ν12 −3 −2 ν22
which results in
− 3ν21 − 2ν22 = 0, −3ν21 − 2ν22 = 0
−2
−2 ν22 −2
ν21 = ν22 ν2 ⇒ 3 = ν22
3 ν22 3 3
6. For the following matrix determine the state transition matrix Φ(t), eigenvectors
−6 2
A=
−6 1
Chapter 1 62
Daniel Abebe Modern Control Systems Lecture Note
Solution
λ 0 −6 2 λ + 6 −2
det( − = det = 0
0 λ −6 1 6 λ−1
(λ + 6)(λ − 1) + 12 = 0
λ2 5λ + 6 = 0 ⇒ (λ + 3)(λ + 2) = 0
λ1 = −3 λ2 = −2
(λI − A) ν 1 = 0
ν11 −6 2 ν
−3 = 11
ν12 −6 1 ν12
which results in
(λI − A) ν 2 = 0
[λI − A] ν 1 = 0
ν11 −6 2 ν
−2 = 21
ν12 −6 1 ν22
which results in
Chapter 1 63
Daniel Abebe Modern Control Systems Lecture Note
Solution:
The characteristic equation of the matrix is
λ+1 0 0
det 0 0 = 0
λ+1
−1 0 λ+1
(λ + 1)(λ + 1)(λ + 1) = 0
⇒ λ1 = λ2 = λ3 = −1
When any matrix has a multiplicity of eigenvalues, the method we don’t rely on is
using the modality matrix since the determinant becomes zero. So let us use the
Chapter 1 64
Daniel Abebe Modern Control Systems Lecture Note
8. For the following dynamic system, find eigenvalues, zero-input response, and zero-
state response to step input for t > 0.
−1 −3 5 h i
ẋ = x + u, y = 3 4 x
0 −2 6
h i⊤
with initial condition x(0) = 4 5 .
Solution:
The characteristic equation of the system is,
Chapter 1 65
Daniel Abebe Modern Control Systems Lecture Note
Solution:
using Lyapunov stability test: use R = I
AP + PA⊤ = −R
3 −2 p11 p12 p p12 3 1 −1 0
+ 11 =
1 −1 p21 p22 p21 p22 −2 −1 0 −1
Chapter 1 66
Daniel Abebe Modern Control Systems Lecture Note
−13
P is symmetric but not positive definite because det(P) = < 0. Hence, the A
16
is not stable. Checking the eigenvalues
√ √
||λI − A|| = 0 ⇒ (λ − 3)(λ + 1) + 2 = 0 ⇒ λ1 = 1 − 2, λ2 = 1 + 2
both eigenvalues has strictly positive real parts, so A is unstable.
10. For dc motor model derive on 1.7.1, check stability with Ra = 0.25, J = 0.0025, La =
0.003, kb = 0.0067, kt = kb .
Solution:
using MATLAB command [ν, D] = eig(A), give us the eigenvaluesλ1 = −83.25, λ2 =
−8.08 and corresponding eigenvectors of the dc-motor dynamics
−1 0.03
ν1 = ν2 =
0.034 −1
Since all eigenvalues have strictly negative real parts, the dynamic system is stable.
11. Check stability of IPC dynamic system derived on 1.6.3, with parameters mp =
0.5, mc = 0.2, b = 0.1, l = 0.3, j = 0.5mp l2 , g = 9.81.
Solution:
using MATLAB command [ν, D] = eig(A), give us the eigenvalues λ1 = 0, λ2 =
−0.14 λ3 = −5.21, λ4 = 5.17 and corresponding eigenvectors of the IPC dynamic
system
1 −0.99 −0.01661 0.0158
0 0.1414 0.0862 0.0816
ν1 = ν2 =
ν3 =
ν4 =
−0.1879
0 0.0021 0.1891
0 −0.0003 0.9782 0.9784
Since not all eigenvalues doesn’t have strictly negative real parts, the dynamic
system is unstable.
12. Consider the following simple mass measuring system, with system parameters k =
20N/m, ms = 1kg, and d = 9N s/m shown in Figure C.4b. (Mid-Exam 2021)
ms x Scale
k d
The system shows the mass of the object to be measured on the scale. The reading
on the scale is proportional to the mass of the object(y = αmo ). Where α is a
proportionality constant (for this problem, use α = 1).
Chapter 1 67
Daniel Abebe Modern Control Systems Lecture Note
Solution:
ms x
fk fd
u − kx − dẋ = ms ẍ
−k d u
ẍ = x− ẋ +
ms ms m
Let x1 = x, x2 = ẋ = ẋ1 , then the state space representation of the system is:
ẋ 0 1 x1 0
1 = −k −d + 1 u
ẋ2 x2
ms ms m
u
The output is y = αmo , but mo = (since u is the weight of the measured
g
αu
mass), which means y = . The scale shows the result after the transient
g
response dies out (under steady-state conditions). At steady-state conditions,
Chapter 1 68
Daniel Abebe Modern Control Systems Lecture Note
the only forces on the system are the spring force and the weight of the system.
Hence,
u = mg = kx
αkx αkx1
y= ⇒ y=
g g
x1
αk
y= 0
g x2
Chapter 1 69
Daniel Abebe Modern Control Systems Lecture Note
s+9 1
(s + 5)(s + 4) (s + 5)(s + 4)
= L−1 −20 s
(s + 5)(s + 4) (s + 5)(s + 4)
20 5 −1 1
+ +
= L−1 s−20
+ 5 s + 4 s + 5 s + 4
5 5 −4
+ +
s+5 s+4 s + 5 s + 4
−4e−5t + 5e−4t −e−5t + e−4t
=
−20e−5t + 20e−4t −5t
5e − 4e −4t
Chapter 1 70
Daniel Abebe Modern Control Systems Lecture Note
(g) The maximum mass of an object that the system can handle is
(h) If we want to increase the maximum mass the system can measure, we have to
replace the spring with another spring having a higher spring constant, because
the spring force is the only force that supports the weight of the object at a
steady state. But this will affect the sensitivity of the system, which has to be
calibrated again. In this case, we have to replace the spring with a constant
kn = 5k = 100N/m.
The eigenvalues of the new system are
λ 0 0 1
||λI − A|| = 0 ⇒ − =0
0 λ −100 −9
λ(λ + 9) + 100 = 0
λ1 = −4.5 + j8.93 λ2 = −4.5 − j8.93
The system is stable, but there will be an oscillation since the poles are complex
conjugates. To solve this problem, we have to use a damper with a higher
damping coefficient. If we take a damper with a damping coefficient d =
20(the minimum value for this case). The eigenvalues of the system will be
λ1 = λ2 = −10. This damper will remove the oscillation.
Chapter 1 71
Daniel Abebe Modern Control Systems Lecture Note
1.17 Exercises
1. Find the characteristic polynomial, eigenvalues, and the corresponding eigenvectors
for each of the following matrices.
1 4 2 2 0
(a) A =
2 3 (e) F = 1 2 1
2 −2 1 2 1
(b) C =
2 2 1 1 1
(f) H =
2 2 1 1 1
(c) D =
−1 −1 1 1 1
0 −1 cos θ sin θ
(d) E = (g) I =
ab a + b − sin θ cos θ
3. Find the general solution to the following initialized systems with step input.
Chapter 1 72
Daniel Abebe Modern Control Systems Lecture Note
(a)
−1 −1 1 h i 0.5
ẋ = x + u, y = 0 1 x, x0 =
2 −4 0 0.5
(b)
2 −1 0 h i 0
ẋ = x + u, y = 1 1 x, x0 =
−1 2 1 0.5
(c)
0 1 1 h i 1
ẋ = x + u, y = 1 0 x, x0 =
1 0 1 1
(d)
1 0 −1 1 0
h i
ẋ = 1 2 1 x + 0 u, y = 0 0 1 x, x0 = 1
2 2 3 0 0
z̈ + ωz = 0
ẋ1 = x1 + 4x2 + u
ẋ2 = 3x1 + 2x2
Find the state transition matrix and the solution of the state equation
is asymptotically stable.
Chapter 1 73
Daniel Abebe Modern Control Systems Lecture Note
is asymptotically stable.
(a)
ẋ1 = −x1 + x2
ẋ2 = (x1 + x2 ) sin k − 3x2
(b)
ẋ1 = −kx1 + x2
ẋ2 = −ax1 − bx2
11. Consider an RLC Circuit in Figure B.3, with system parameters L = 0.1H, C =
0.1F and R = 1Ω. If the designer aims to regulate the voltage across the capacitor.
R L
vs C vo
Chapter 1 74
Daniel Abebe Modern Control Systems Lecture Note
(f) Determine the zero state response of the system under step u = e−2t .
(g) What is the maximum value of voltage the capacitor can be charged to?
(h) If the output is the voltage across the inductor, determine the state space
representation of the system .
(i) Compare the result on 11h with the result on 11a and state their difference.
(j) Determine the zero state response of the state space representation derived in
part 11a and part 11h for u = cos t.
Chapter 1 75
Chapter 2
Definition:-A linear time invariant system with state space representation on equation
1.3.6 and equation 1.3.7 is completely controllable if for any t0 any initial state x(t0 ) = x0
can be transferred to a final state x(tf ) at tf using control input u (i.e. if there exist a
control signal u that transfer a state of a dynamic system from any initial condition x0
to a final state x(tf )).
76
Daniel Abebe Modern Control Systems Lecture Note
Z t1
Letting α(τ )u(τ )dτ = βk .
0
we have
k=n−1
X
x0 = − Ak Bβk (2.1.4)
k=0
β
1
i
β2
.. .. .. ..
h
2 n−1
x0 = − B . AB . A B . · · · . A B .
(2.1.5)
..
βk
If the system is completely controllable, given any initial state x0 , equation 2.1.4 must
be satisfied, which requires the matrix
. . . .
h i
C o = B .. AB .. A2 B .. · · · .. An−1 B (2.1.6)
to have full rank. This imply that the vectors B, AB, A2 B, · · · , An−1 B should be linearly
independent. The matrix in equation 2.1.6 is called the controllability matrix. If the
controllability matrix has full rank, then the system defined by (A, B) is completely
controllable.
Or if a matrix Q is right Lyapunov and the solution for Q of equation 2.1.7
Chapter 2 77
Daniel Abebe Modern Control Systems Lecture Note
Definition:- A linear time invariant system with state space representation on equation
1.3.6 and equation 1.3.7 is completely observable if for any t0 and any initial state x0 ,
there exist finite time t1 > t0 such that the knowledge of u(·) and y(·) for t ∈ [t0 , t1 ]
suffices to determine x0 properly.
From the solution to the output equation
Z t1
At
y(t) = Ce x0 + CeA(t−τ ) Bu(τ )dτ + Du(t) (2.2.1)
0
Knowing u(t) and y(t) will give information about x0 , if the dynamic system defined by
(A, C) is completely observable. Let the input be zero to derive the required condition.
This means the system response is only the response to an initial condition. So,
y(t) = CeAt x0 (2.2.2)
k=n−1
X
At
Replacing e = α(t)Ak
k=0
k=n−1
X
y(t) = C α(t)Ak x0
k=0
y(t) = α0 (t)Cx0 + α0 (t)Cx0 + α0 (t)CA2 x0 + · · · + α0 (t)CAn−1 x0
C
CA
h i
2
y(t) = α0 (t) α1 (t) · · · αn−1 (t) CA x0
(2.2.3)
..
.
n−1
CA
to have full rank. This implies that the co-vectors1 C, CA, CA2 . . . CAn−1 should be
linearly independent. The matrix in equation 2.2.4 is called the observability matrix. If
1
A co-vector (or dual vector ) is an element of the dual space V ∗ of a vector space V . It is a linear
functional that maps vectors to scalars.
Chapter 2 78
Daniel Abebe Modern Control Systems Lecture Note
the observability matrix has full rank, then the system defined by (A, C) is completely
observable.
Or if a matrix P is left Lyapunov and the solution for P of the equation 2.2.5
A⊤ P + QA + C⊤ C = 0 (2.2.5)
is positive definite then the system defined by (A, C) is completely observable. P is de-
termined from co-variance of y(t) which is called observability grammian matrix, defined
as
Z ∞ Z ∞
⊤ ⊤
Q= y (t)y(t)dt = (eA t C⊤ CeAt )dt (2.2.6)
0 0
Consider a linear time invariant system on equation 1.3.6 and equation 1.3.7, if there
is a transformation V that converts a matrix A to a diagonal matrix Λ. Where
V−1 ΛV = A (2.3.1)
The system is completely controllable if and only if no left (row) eigenvectors of A are
orthogonal to B. If ν is left eigenvector of matrix A corresponding to any eigenvalue λ,
for the system to be completely controllable
ν · B ̸= 0 (2.3.3)
The system is completely observable if and only if no right (column) eigenvectors of A are
orthogonal to C. If q is a right eigenvector of matrix A corresponding to an eigenvalue
λ, for the system to be completely observable
C · q ̸= 0 (2.3.4)
These tests are also called PHB (Popov Belvitch Hautus) tests.
Chapter 2 79
Daniel Abebe Modern Control Systems Lecture Note
u + xc y
Bc ʃ Cc
+
+ Ac
+
A12
Ac
+ xc y
+
0 ʃ Cc
Figure 2.1: Block Diagram Representation of Controllable and Uncontrollable Modes of a System
Similarly, if a system is not completely observable, the system can be decomposed into
its observable states and unobservable states using a similarity transformation. Con-
sider LTI system having r observable states xo , and n − r unobservable states xō , with a
transformation matrix T, the system can be transformed into
ẋo Ao 0 x B
= o + o u (2.4.3)
ẋō A21 Aō xō Bō
h i xo
y = Co 0 + Du (2.4.4)
xō
The Block diagram representation of observable and unobservable modes decomposition
is shown in Figure 2.2.
Chapter 2 80
Daniel Abebe Modern Control Systems Lecture Note
u + xo y
Bo ʃ Co
+
Ao
A21
+ +
Ao
+ + xo
Bo ʃ 0
Figure 2.2: Block Diagram Representation of Observable and Unobservable Modes of a System
• The observable subspace is defined by the right eigenvectors qi for which Cqi ̸=
0. These modes can be detected through the output matrix C.
• The unobservable subspace comprises right eigenvectors qi where Cqi = 0,
indicating that these states cannot be inferred from output measurements.
Chapter 2 81
Daniel Abebe Modern Control Systems Lecture Note
If the unobservable mode of the system is stable, or conversely if the unstable mode is
observable, then the system is said to be detectable.
Stabilizability and detectability are called weak conditions since it doesn’t impose tight
conditions like controllability and observability on system dynamics. If the system state
is uncontrollable and stable, there is no need to put additional restrictions since the rest
of the states are controllable. Similarly, as long as the unstable states are observable,
there is no need to add a restriction on unobservable stable states.
Stability of uncontrollable states can be deduced from the eigenvalues of the uncontrol-
lable subsystem matrix (Ac̄ , while the stability of unobservable dynamics is deduced
from the eigenvalues of the unobservable subsystem matrix Aō ). If the matrices have
eigenvalues with a negative real part, then the dynamics are stable.
1. Column Space (or Range): The column space, denoted Col(A), consists of all
possible linear combinations of the columns of A. It lives in Rm and represents all
possible outputs the matrix can produce when multiplied by a vector. Its dimension
is equal to the rank of A, which tells us how many linearly independent columns
A has. Essentially, it shows us where the matrix can ”reach” in the output space.
2. Null Space (or Kernel): The null space, denoted Null(A), contains all vectors
x in Rn that get mapped to the zero vector when multiplied by A. It shows us
which inputs completely disappear when transformed by A. The dimension of
the null space is called the nullity and equals n − rank(A). This tells us how much
information gets ”lost” when the matrix transforms vectors. The null space consists
of all eigenvectors for λ = 0 (plus the zero vector) for a square matrix.
Ax = 0.
Chapter 2 82
Daniel Abebe Modern Control Systems Lecture Note
3. Row Space: The row space is the span of all row vectors of A and lives in Rn .
Surprisingly, it always has the same dimension as the column space - both equal to
the rank of A. This space represents all possible linear combinations of the matrix’s
rows, which corresponds to all possible ”measurements” or ”constraints” that the
matrix’s rows can represent.
4. Left Null Space: This less obvious space, denoted Null(A⊤ ), consists of all vectors
y in Rm that get mapped to zero when multiplied by A⊤ (the transpose of A).
These vectors are orthogonal to all columns of A, and the space has dimension
m − rank(A). It represents all possible ”constraints” or ”tests” that vectors in the
output space must satisfy to be in the column space of A.
A⊤ y = 0.
Example
For matrix:
1 0 1
B= (2.4.5)
0 1 1
x1 + x3 = 0
x2 + x3 = 0
Chapter 2 83
Daniel Abebe Modern Control Systems Lecture Note
Note: In linear algebra, when a matrix is transformed into its reduced row echelon
form (RREF), the variables corresponding to pivot columns are called basic vari-
ables, while those corresponding to non-pivot columns are free variables. The basic
variables can be expressed uniquely in terms of the free variables, which represent
the independent degrees of freedom in the solution space.
−1 −1
Null(B) = x3 −1 | x3 ∈ R → span −1 (2.4.8)
1
1
• Column Space The column space Col(B) is the span of the column vectors of B.
The columns are:
1 0 1
, , (2.4.9)
0 1 1
Since the first two columns are linearly independent and the third is their sum, the
column space is:
1 0
Col(B) = span , = R2
(2.4.10)
0 1
• Row Space The row space Row(B) is the span of the row vectors of B.
The rows are:
h i h i
1 0 1 , 0 1 1 (2.4.11)
• Left Null Space The left null space Null(B⊤ ) is the set of solutions to B⊤ y = 0.
Compute B⊤ :
1 0
⊤
B = 0 1 (2.4.13)
1 1
Solve:
1 0 0
y
1
= 0 (2.4.14)
0 1
y2
1 1 0
Chapter 2 84
Daniel Abebe Modern Control Systems Lecture Note
This gives:
y1 = 0
y2 = 0
y1 + y2 = 0
ẋ = Ax + Bu (2.4.16)
y = Cx + Du (2.4.17)
1. Controllability: The controllability matrix Co is defined as:
h i
C o = B AB A2 B · · · An−1 B (2.4.18)
Controllability Decomposition
The concept of controllability fundamentally concerns whether a system can be driven
from any initial state x(0) to any desired final state x(tf ) in finite time through appro-
priate control inputs u(t). This property is characterized by the reachable subspace of
the system by control inputs u(t), which corresponds precisely to the column space of
the controllability matrix C o .
The state space decomposition into controllable and uncontrollable subspaces can be ob-
tained through a similarity transformation matrix Tc constructed from the basis vectors
spanning Col(C o ) (the controllable subspace) and the basis vectors spanning Null(C ⊤
o)
(the uncontrollable subspace, being the orthogonal complement).
h i
⊤
Tc = basis(Col(C o )) basis(Null(C o )) (2.4.20)
Chapter 2 85
Daniel Abebe Modern Control Systems Lecture Note
Observability Decomposition
The concept of observability fundamentally addresses whether a system’s initial state
x(0) can be uniquely determined from knowledge of the input u(t) and output y(t) over
a finite time interval [0, tf ]. This property characterizes which state information is visible
through the output measurements, defining the observable subspace of the system.
The state space decomposition into observable and unobservable components can be
obtained through a similarity transformation matrix To constructed from the basis vec-
tors spanning Col(O ⊤v ) (the observable subspace, equivalent to the row space of O v ) and
the basis vectors spanning Null(O v ) (the unobservable subspace).
h i
To = basis(Col(O ⊤ v )) basis(Null(O v )) (2.4.21)
For a system to be completely controllable, adj(sI − A)B and det(sI − A) shouldn’t have
common factors (i.e, there shouldn’t be a pole-zero cancellation in a transfer function
from input to state of the system).
Consider a transfer function from U (s) to Y (s)
For a system to be completely observable, Cadj(sI − A)B and det(sI − A) shouldn’t have
common factors (i.e, there shouldn’t be a pole-zero cancellation in the transfer function
from input to output of the system).
2.6 Examples
1. Consider the following dynamic system, check whether the system is completely
controllable and completely observable.
−1 −3 5 h i
ẋ = x + u, y = 3 4 x
0 −2 6
Solution:
the controllability matrix C o is
h 5 −23
i
C o = B AB =
6 −12
Chapter 2 86
Daniel Abebe Modern Control Systems Lecture Note
5 −23
α1 + α2 = 0 ⇒ α1 = α2 = 0
6 −12
They are linearly independent, hence the controllability matrix has full rank, which
means the system is completely controllable.
3 4
α1 + α2 = 0 ⇒ α1 = α2 = 0
−3 −17
They are linearly independent, hence the observability matrix has full rank, which
means the system is completely observable.
2. Check controllability and observability of the following dynamic system.
−2 −2 −2 1
h i
ẋ = 1 x + u y = 1 1 0 x
0 0 0
1 1 −1 0
The controllability matrix is,
1 −2 0
C o = 0 1 −2
0 1 −2
The rank rank(C o ) = 2, hence the dynamic system has an uncontrollable state.
The observability matrix is,
1 1 0
O v = −1 −2 −2
−2 0 4
having rank rank(O v ) = 2, hence the dynamic system have unobservable state.
The transformation matrix Tc that decomposes the system into controllable and
uncontrollable subspaces is
h i
Tc = basis(Col(C o )) basis(Null(C ⊤
o ))
Chapter 2 87
Daniel Abebe Modern Control Systems Lecture Note
The basis of (Col(C o )) is any linearly independent columns of C o , taking the first
two columns and the basis vector that span the null space of C ⊤
o , the transformation
matrix Tc becomes,
1 −2 0
Tc = 0 1 −1
0 1 1
To solve for T−1c , transform the augmented matrix [Tc | I] into [I | F ] using the
Gauss-Jordan elimination method, where I is the identity matrix of the same size
as Tc . The resulting matrix F will be the inverse of Tc . Using this method, the
inverse of Tc is
2 2 2
T−1 = 0.5
c 0 1 1
0 −1 1
With this transformation matrix, the decomposed system system matrix Ā, input
matrix B̄ and C̄ is,
..
0 −4 . −2 1
..
1 −2 −1
. 0
Āc = T−1 AT =
,
B̄c = T−1 B =
. . . . . . . . . . . . . . .
..
0 0 . −1 0
.
h i
C̄c = CT = 1 −1 .. −1
Chapter 2 88
Daniel Abebe Modern Control Systems Lecture Note
from the decomposition, we see that the uncontrollable state is stable (with eigen-
value −1), hence the system is stabilizable.
The transformation matrix To that decomposes the system into controllable and
uncontrollable subspaces is
h i
To = basis(Col(O ⊤
v )) basis(Null(O v ))
With this transformation matrix, the decomposed system system matrix Ā, input
Chapter 2 89
Daniel Abebe Modern Control Systems Lecture Note
The transformation matrix Tc that decomposes the system into controllable and
uncontrollable subspaces is
h i
Tc = basis(Col(C o )) basis(Null(C ⊤
o ))
Chapter 2 90
Daniel Abebe Modern Control Systems Lecture Note
The basis of (Col(C o )) is any linearly independent columns of C o , taking the first
two column and the basis vector that span the null space of C ⊤
o , the transformation
matrix Tc and T−1
c becomes,
1 1 0 0 −1 0
Tc = 1 0 0 T−1 = −1 1
c 0
0 0 1 0 0 −1
With this transformation matrix, the decomposed system system matrix Ā, input
matrix B̄ and C̄ is,
..
0 1 . 1 1
..
. −1
1 0 0
Āc = T−1 AT =
, B̄c = T−1 B =
. . . . . . . . . . . . . . .
..
0 0 . −2 0
.
h i
C̄c = CT = 1 1 .. 0
From the decomposition, we see that the uncontrollable state is stable (with eigen-
value −2), hence the system is stabilizable.
The transformation matrix To that decomposes the system into controllable and
uncontrollable subspaces is
h i
To = basis(Col(O ⊤
v )) basis(Null(O v ))
With this transformation matrix, the decomposed system system matrix Ā, input
Chapter 2 91
Daniel Abebe Modern Control Systems Lecture Note
from the decomposition we see that the unobservable states are stable(with eigen-
values −1 and −2 ), hence the system is detectable.
4. Check the detectability and stabilizability of the dynamic system
−2 0 −3 −1
ẋ = 0 0 0 x + 0 u
3 5
0 −2
2 2
21
y= 1 0 x
2
The controllability matrix of the system is
−1 8 3.5
Co = 0
0 0
−2 −6.5 −4.25
The rank of the controllability matrix is 2, hence the system has one uncontrollable
dynamics.
The transformation matrix Tc that decomposes the system into controllable and
uncontrollable subspaces is
h i
Tc = basis(Col(C o )) basis(Null(C ⊤
o ))
Chapter 2 92
Daniel Abebe Modern Control Systems Lecture Note
The basis of (Col(C o )) is any linearly independent columns of C o , taking the first
two column and the basis vector that span the null space of C ⊤
o , the transformation
matrix Tc and T−1
c becomes,
−1 8 0 6.5 0 8
−1 −1
Tc = 0 T = −2
0 0 c 0 1
22.5
−2 −6.5 1 0 −22.5 0
With this transformation matrix, the decomposed system system matrix Ā, input
matrix B̄ and C̄ is,
..
0 −11.25 . 0 −22.5
..
−22.5 −11.25 .
−1 −1 0 −1 −1 0
Āc = T AT = , B̄c = T B =
22.5
... ...
. . . . . . 22.5
...
..
0 0 . 0 0
.
h i
C̄c = CT = −1 8 .. 10.5
From the decomposition, we see that the uncontrollable state is marginally (with
eigenvalue 0), hence the system is stabilizable.
which has rank 2, hence the dynamic system has one unobservable mode.
Chapter 2 93
Daniel Abebe Modern Control Systems Lecture Note
The transformation matrix To that decomposes the system into controllable and
uncontrollable subspaces is
h i
To = basis(Col(O ⊤
v )) basis(Null(O v ))
With this tranformation matrix, the decomposed system system matrix Ā, input
matrix B̄ and C̄ is,
..
−56 −77 . 0 2
..
−5
−1 119 140 . 0
Āo = T−1 AT = , B̄o = T−1 B = −1
21
. . . . . . . . . . . . 21
. . .
.
91 112 .. 42 −4
.
h i
C̄o = CT = 5 2 .. 0
From the decomposition, we see that the unobservable state is stable(with eigen-
value −2), hence the system is detectable.
6. Given the following LTI system, determine the value of α for which the system is
not c.c.
2 α−1 1
ẋ = x + u
0 2 α
Solution:
The controllability matrix of the system is
1 2 + α(α − 1)
Co =
α 2α
Chapter 2 94
Daniel Abebe Modern Control Systems Lecture Note
The determinant of the controllability matrix must be non-zero if the system is c.c.
So, to find the value of α that makes the system not c.c.
det(C o ) = 0 ⇒ 2α − α3 + α2 − 2α = 0
⇒ α2 (1 − α) = 0 ⇒ α = 0, or α = 1
The same will hold for observability of the system det(O v ) ̸= 0 if the system is c.o.
2.7 Exercises
1. Verify that the following control system is completely controllable
0 1 0 0
ẋ = x + u
0 0 1 −1
(b)
2 α−3 1 1
ẋ = x + u
2
0 2 0 α −2
Chapter 2 95
Daniel Abebe Modern Control Systems Lecture Note
(b)
0 0 0 h i
ẋ = x+ u y= 1 1 x
0 −1 1
(c)
4 3 5
−2 1 h i
ẋ = 1 −2 −3 x + u y = 0 1 1 x
−1 0
2 1 8
(d)
0 1 5
0 h i
ẋ = 0 1 x + u y = 1 0 0 x
0
1
−18 −15 −2
Chapter 2 96
Chapter 3
We know that we can affect the performance and behavior of a dynamic system by
changing the location of the closed-loop poles of the system. We can cancel the unstable
poles by multiplying the system transfer function by another transfer function having a
zero at the location of the unstable pole of the system. This other system is called a
controller or compensator.
Consider a completely controllable system on equation 1.3.6 and equation 1.3.7, suppose
we apply a linear feedback control law, which means that the control law is a linear
combination of state variables,
u = −Kf x (3.0.1)
where Kf is m × n feedback matrix. The closed-loop system under this control law
becomes,
ẋ = Ax + B(−Kf x)
ẋ = (A − BKf )x (3.0.2)
The block diagram representation of the system with a full state feedback controller is
shown in Figure 3.1.
ud + u + x y
- B ʃ C
+
Kf
Figure 3.1: Block Diagram Representation of System with Full State Feedback Controller
97
Daniel Abebe Modern Control Systems Lecture Note
Pole shifting theorem:- for any controllable system, it is possible to obtain arbi-
trary eigenvalues for matrix A − BKf , using suitable feedback matrix Kf , subjected to
the only obvious constraint that complex eigenvalues must appear in pair.
The state feedback controller design can be done either by selecting the location of the
poles and determine the feedback gain Kf that locate the closed-loop poles of the system
at the desired location, or determining the gain Kf that determine the optimal location
of the closed-loop poles and that minimize u at the same time.
The pole placement or eigenvalue assignment problem can be define as: let λ1 , λ2 , . . . , λn
be the eigenvalues of the matrix A of the open loop LTI system and λ̃1 , λ̃2 , . . . , λ̃n be
desired eigenvalues of the matrix A − BKf of the closed loop system. Where all complex
eigenvalues exist in complex conjugate pairs. Let p(s) and p̂(s) be respective character-
istics polynomial of matrix A and A − BKf respectively.
n
Y
p(s) = (s − λi ) = ∥sI − A∥ = sn + a1 sn−1 + . . . + an−1 s + an (3.1.1)
i=1
Yn
p̂(s) = (s − λ̃i ) = ∥sI − A + BKf ∥ = sn + â1 sn−1 + . . . + ân−1 s + ân (3.1.2)
i=1
Kf is determined if and only if the open loop system is controllable (i.e., the controllability
matrix C o has full rank).
The feedback matrix Kf can be determined by one of the following methods.
Kf = (W⊤ C ⊤ −1
o ) (â − a) (3.1.3)
Where
1 a1 a2 a3 . . . an−1 â1 a
1
0 1 a1 a2 . . . an−2 â2 a2
W = 0 0 1 a1 . . . an−3 , â = â3 , a = a3
.. .. .. .. .. .. .. ..
. . . . . . . .
0 0 0 0 ... 1 ân an
Chapter 3 98
Daniel Abebe Modern Control Systems Lecture Note
2. The phase canonical form formula: Consider the case where A is in controllable
canonical form Ac , bc .
0 1 0 ... 0 0
0 0 1 ... 0 0
.. .. .. .. .. .
Ac = . . . . . , bc = ..
0 0 0 ... 1 0
−a0 −a1 −a2 . . . −an1 1
Which means:
Kf = (W⊤ C ⊤ −1 ˜
o ) (â − a) = I = (â − a)
Kf = (â − a) (3.1.4)
3. Ackerman’s formula:
⊤ −1
K⊤
f = e C o p̂(A) (3.1.5)
Where
e⊤ = [0, 0, . . . , 1],
and p̂(A) is obtained from polynomial on equation 3.1.1 by replacing s with matrix
A, which means
An LTI system defined by (A, B) is stabilizable if there exists a matrix Kf such that
(A − BKf ) has a strictly negative real part.
Chapter 3 99
Daniel Abebe Modern Control Systems Lecture Note
ud + u + x y
- B ʃ C
+
Kf
Figure 3.2: Block Diagram Representation of Full State Feedback with Reference State
Chapter 3 100
Daniel Abebe Modern Control Systems Lecture Note
r + + u + x y
ʃ k1 - B ʃ C
- +
Kf
A − BKf Bk1 A 0 B h i
⇒ − Kf −k1 (3.3.6)
−C 0 −C 0 0
which means
A 0 B
, (3.3.7)
−C 0 0
To find the value of k1 , we can analyze the steady state response of the system. At steady
state, the output of the system should follow the reference input y = r, which means
that ξ̇ = 0. This implies that the output of the integrator is a constant value equal to
r. Therefore, the states of the system should also follow the desired states of the system
xd , given by xd = k1 r.
However, at this point, the output of the system is given by y = r = Cxd + Du, where
u = xd − Kf xd . Thus, we can write:
Chapter 3 101
Daniel Abebe Modern Control Systems Lecture Note
ωn2
G(s) = (3.4.1)
s2 + 2ξωn s + ωn2
Since any second-order system can be written in this form, we can determine the desired
poles of the system from the time domain specification of the system.
Note:- The goal of the designer is to find Kf such that a closed loop time domain
specification of the given second order system meets the requirement.
Chapter 3 102
Daniel Abebe Modern Control Systems Lecture Note
3.5 Examples
1. Given a continuous-time system
1 2 0
ẋ = + u
3 1 1
h i
y= 1 3 x
Solution
Check for controllability
h 0 2
i
C o = B AB =
1 1
Chapter 3 103
Daniel Abebe Modern Control Systems Lecture Note
Solution
We have to cancel zero at z = −2 by a pole at p = −2, cancel the pole at p = −1
by a zero at z = −1, cancel the pole at p = 2 by a zero at z = 2, and add one pole
at p = −2 to get G(s) from H(s).
(s + 1)(s − 2) (s + 1)(s − 2) (s − 1)(s + 2)
G(s) = H(s) =
(s + 2)2 (s + 2)2 (s + 1)(s − 2)(s + 3)
(s − 1)(s + 2)
=
(s + 2)2 (s + 3)
Since the degree of the numerator is less than the degree of the denominator, it is
possible to have such a kind of compensator. Hence it is possible to get G(s) from
H(s).
The state-space representation of H(s) is
−2 5 6 1
ẋ = 1 0 0 x + 0 u
0 1 0 0
h i
y = 1 1 −2 x
We need a state feedback control law u = −Kf x that converts the system H(s) to
G(s), hence the desired location of poles of the closed-loop system is at p = −2, p =
−2, p = −3. Using the Base-Gura method,
det(sI − A + BKf ) = (s + 2)2 (s + 3)
s 0 0 −2 5 6 1
h i
det 0 s 0 − 1 0 0 + 0 k1 k2 k3 = s3 + 7s2 + 16s + 12
0 0 s 0 1 0 0
s + 2 + k1 −5 + k2 −6 + k3
3 2
det −1 0 = s + 7s + 16s + 12
s
0 −1 s
s3 + (2 + k1 )s2 + (−5 + 6)s + (−6 + k3 ) = s3 + 7s2 + 16s + 12
Chapter 3 104
Daniel Abebe Modern Control Systems Lecture Note
2 + k1 = 7 ⇒ k1 = 5, −5 + 6 = 16 ⇒ k2 = 21, −6 + k3 = 12 ⇒ k3 = 18
h i h i
Kf = k1 k2 k3 = 5 21 18
The closed loop dynamics of the system with full state feedback gain matrix Kf
are
ẋ = (A − BKf )x
ẋ −2 5 6 5 21 8 x
1 1
ẋ2 = 1 0 0 − 0 0 0 x2
ẋ3 0 1 0 0 0 0 x3
ẋ −7 −17 −2 x1
1
ẋ2 = 1
0 0 x2
ẋ3 0 1 0 x3
3. For IPC dynamic system on 1.6.3 with parameters mc = 0.5, mp = 0.2, b = 0.1, l =
0.3, I = 0.5mp l2 , g = 9.81, Find a feedback gain matrix that place the poles of the
closed loop system at p1 = −1, p2 = −2, p3 = −1 − i, p4 = −1 + i.
Solution
To check for complete controllability of IPC dynamic system, use C o = ctrb(A, B),
and rank(C o ) commands on MATLAB to check controllability and the rank of
controllability matrix, respectively. With the given parameters
0 1.00 0 0 0
0 −0.17 2.31
0 1.75
ẋ =
x +
u
0 0 0 1.00 0
0 −0.39 26.93 0 3.92
h i
y= 1 0 1 0 x
Chapter 3 105
Daniel Abebe Modern Control Systems Lecture Note
The closed loop dynamics of the system with full state feedback gain matrix Kf
are
ẋ = (A − BKf )x
ẋ1 0 1.00 0 0 0 0 0 0 x
1
ẋ2 0 −0.17 2.31 0 −0.18 −0.64 16.7 2.46 x2
= −
ẋ3 0 0 0 1.00 0 0 0 0 x3
ẋ4 0 −0.39 26.93 0 −0.41 −1.41 37.11 5.46 x4
ẋ 0 1 0 0 x
1 1
ẋ2 0.18 0.46 −14.39 −2.46 x2
=
ẋ3 0 0 0 1 x3
ẋ4 0.41 1.02 −10.18 −5.46 x4
Figure 3.4: Closed Loop Response of IPC for Impulse Input with Kf
4. For the following system, determine a feedback gain matrix that places the closed-
loop poles of the system at p = −2 ± 2i.
1 0 0 1
ẋ = 1 −1 1 x + 1 u
0 0 −2 0
h i
y= 1 0 0 x
Solution
As shown in chapter two example 3, the system is not completely controllable. So
we have to decompose the system into controllable and uncontrollable modes as
shown in chapter two example 3. So we can design a state feedback controller for
Chapter 3 106
Daniel Abebe Modern Control Systems Lecture Note
(s − 1 + k1 )(s + 1 + k2 ) − k2(−1 + k1 ) = s2 + 4s + 8
s2 + (k1 + k2 )s − 1 + k1 = s2 + 4s + 8
k1 − 1 = 8 ⇒ k1 = 9 k1 + k2 = 4 ⇒ k2 = −5
h i
Kf = 9 −5
What happens if we want to design a full state feedback controller for as a whole?
Let us assume the third pole is at p = α, then
which results in
Hence, the system requires the stability of the uncontrollable mode only. k3 does
not affect the system. If the uncontrollable mode is stable, the feedback gain corre-
sponding to the uncontrollable mode does not affect the system closed closed-loop
dynamics. This example proves this scenario.
The system only requires the stability of uncontrollable mode; no need to design
feedback gain Kf for the uncontrollable mode, even if the uncontrollable mode is
not stable. We cannot impose any effect on the uncontrollable mode (that is why
Chapter 3 107
Daniel Abebe Modern Control Systems Lecture Note
The closed loop dynamics of the system with full state feedback gain matrix Kf
are
ẋ = (A − BKf )x
x1
ẋ 1 0 0 9 −5 k3
1 x2
ẋ2 = 1 −1 1 − 9 −5 k3
. . .
ẋ3 0 0 −2 0 0 0
x3
.
−8 5 .. −k3 x
ẋ1 1
.
−8 4 .. 1 − k x
3
2
ẋ2 =
. . . . . . ...
. . . . . .
ẋ3
..
0 0 . −2 x3
5. Determine a full state feedback matrix for the following second-order system such
that the closed loop system has settling time ts = 4sec and peak time tp = 1.813sec.
1 0 1
ẋ = x + u
1 1 0
h i
y= 1 0 x
Solution
From given time domain specifications
4 4
ξωn = = = 1
ts 4
π 3.14
ωd = = = 1.732
tp 1.813
p
ωn 1 − ξ 2 = 1.732 ⇒ 1 − ξ 2 = 3ξ 2
⇒ ξ = 0.5 ⇒ ωn = 2
Thus the closed loop √ system has characteristics
√ equation of s2 + 2ξωn + ωn2 =
s2 + 2s + 4 = (s + 1 + 3i)(s + 1 − 3i). The state feedback controller gain must
locate the closed-loop system poles at s = −1 ± i. We determine Kf as follows
∥sI − A + BKf ∥ = s2 + 2s + 4
s−1 0 k1 k2
+ = s2 + 2s + 4
−1 s − 1 0 0
s − 1 + k1 k2
= s2 + 2s + 4
−1 s−1
s2 + (−2 + k1 )s + 1 − k1 + k2 = s2 + 2s + 4
− 2 + k1 = 2 ⇒ k1 = 4 1 − k1 + k2 = 4 ⇒ k2 = 7
Chapter 3 108
Daniel Abebe Modern Control Systems Lecture Note
h i
Kf = 4 7
The closed loop dynamics of the system with full state feedback gain matrix Kf
are
ẋ = (A − BKf )x
ẋ 1 0 4 7 x
1 = − 1
ẋ2 1 1 0 0 x2
ẋ −3 −7 x1
1 =
ẋ2 1 1 x2
3.6 Exercises
1. Consider a system
−1 −2 1
ẋ = x + u
2 −4 3
find a feedback gain matrix Kf such that the closed loop system has eigenvalues at
−4 and −5.
2. For LTI system
1 0 −1 1
ẋ = 1 2 1 x + 0 u
2 2 3 1
find state feedback gain matrix Kf such that the eigenvalues of closed loop system
at −1, −1 ± 2i.
3. For the dynamic system
0 1 0 1 0
ẋ = 0 0 1 x + 0 1 u
6 −11 6 1 1
Find the state feedback gain matrix Kf such that the eigenvalues of the closed loop
system are at 1,2, and 3, then check for stability of the system with the control law.
4. Determine whether the system described by
−1 0 3 1
ẋ = 0 −3 0 x + 1 u
1 0 −3 −1
Chapter 3 109
Daniel Abebe Modern Control Systems Lecture Note
is c.c. and show that under linear feedback of the form u = αx1 + βx3 , the closed
loop system has two fixed eigenvalues, one of which is equal to −3. Determine the
second fixed eigenvalues and also values of α and β such that the third closed loop
eigenvalue is −4.
5. Determine the range of values of k, the linear dynamic system
0 1 0
ẋ = 0 1 x
0
−k −1 −2
h i⊤
is asymptotically stable. If k = −1 and control term 0 1 0 u is added, find
linear feedback gain Kf that makes the closed loop system to have poles at −1.
6. For the control system described by
ẋ1 = x1 − x3 + u
ẋ2 = x1 + 2x2 + x3
ẋ3 = 2x1 + 2x2 + 3x3 + u
find a suitable feedback matrix Kf such that the closed loop system has eigenvalues
−1 and −1 ± 2i.
7. For the control system described by
ẋ1 = x2 + u1
ẋ2 = x3 + u2
ẋ3 = 6x1 − 11x2 + 6x3 + u1 + u2
find a suitable feedback matrix Kf such that the closed loop system has eigenvalues
−1 and −2 ± i.
8. For the following system, find a state feedback controller gain matrix so that the
closed-loop system has maximum percentage overshoot mp = 14% and rise time
tr = 1.25sec.
0 −1 0
ẋ = x + u
−3 −2 1
Chapter 3 110
Chapter 4
To design a full-state feedback control for a given system, all the states of the system must
be available as output (must be measured directly). But in many real-world problems, all
states of the system are not available as output. If the system is completely observable,
the states of the system are related to the output of the system in some way. So we can
recover the unavailable states from the output of the system. The method of recovering
these states from the output of the system is called the observer or state estimator design
problem. If all states are required to be observed or estimated, the problem is a full state
observer problem; if we can measure some states directly and the others are to be observed
or estimated, then the problem becomes a reduced-order state observer/estimator.
If the system is partially observable, we can design an observer for the observable mode
of the system. Since we cannot say anything about the relation between the output of
the system and the unobservable mode of the system, we cannot design an observer for
the unobservable mode of the system.
ẋ = Ax + Bu (4.1.1)
y = Cx (4.1.2)
x̃˙ = ẋ − x̂˙
x̃˙ = Ax + Bu − Âx̂ − B̂u − Lc y
x̃˙ = Ax + Bu − Âx̂ − B̂u − Lc Cx (4.1.4)
111
Daniel Abebe Modern Control Systems Lecture Note
Thus, the design problem is finding a suitable Lc such that the estimation error dynamics
go to zero as t → ∞, or determining a suitable Lc that places the closed loop eigenvalues
of estimation error dynamics x̃˙ = (A − Lc C)x at the desired location.
Replacing  with A − Lc C in the full state observer equation on 4.1.3 we have
Lc = K⊤ (4.1.9)
If the desired poles of x̃˙ = (A − Lc C)x and ẋ = (A − BKf )x are the same. For
a system expressed in observable canonical form Lc equal to the transpose of the full
state feedback matrix K⊤ f , that places the closed-loop eigenvalues of the corresponding
controllable canonical form representation of the same system at the desired location.
An LTI system defined by (A, C) is detectable, if there exists a matrix Lc such that
(A − Lc C) has a strictly negative real part.
Chapter 4 112
Daniel Abebe Modern Control Systems Lecture Note
The block diagram representation of a full-order state observer is shown in Figure 4.1.
u + x +
B + ʃ C
y
+
A
+
Lc
--
x̂
+ ʃ C
ŷ
A
Figure 4.1: Block Diagram Representation of Full Order State Observer
The states to be estimated are related to output so, the output of the state estimator
is η̂ 2 = Lc y + ω, where ω is the state of reduced order dynamics. In terms of η 1 and
η 2 , the dynamic system can be decomposed into directly measured state dynamics and
directly unmeasured state dynamics as follows,
η̇ 1 = A11 η 1 + A12 η 2 + B1 u
η̇ 2 = A21 η 1 + A22 η 2 + B2 u
y = η1 (4.2.1)
Since η 1 states are directly measured, the estimated state is equal to the output, which
means
η̂ 1 = η 1 = y (4.2.2)
Chapter 4 113
Daniel Abebe Modern Control Systems Lecture Note
Let ω̇ = Jω + Ny + Mu, which is called minimum order state observer dynamics. The
estimation error dynamics become
η̃˙ 2 = η̇ 2 − η̂˙ 2
= (A21 η 1 + A22 η 2 + B2 u) − (Lc ẏ + ω̇)
= (A21 η 1 + A22 η 2 + B2 u) − (Lc A11 η 1 + Lc A12 η 2 + Lc Bu + Jω + Nη 1 + Mu)
(4.2.6)
replacing for ω = η 2 − Lc η 1 − η̃ 2 , then we have
η̃˙ 2 =(A21 η 1 + A22 η 2 + B2 u) − (Lc A11 η 1 + Lc A12 η 12 + Lc B1 u)
− (Jη 2 − JLc η 1 − Jη̃ 2 + Nη 1 + Mu)
η̃˙ 2 =(A12 − Lc A + JLc − N)η 1 + (A22 − Lc A12 − J)η 2
+ (B2 − Lc B1 − M)u + Jη̃ 2 (4.2.7)
For estimation error dynamics to go to zero independent of x and u, the following
conditions should be full-filled
A12 − Lc A11 + JLc − N = 0 ⇒ N = A21 − Lc A11 + JLc
B2 − Lc B1 − M = 0 ⇒ M = B2 − Lc B1
A22 − J − Lc A12 = 0 ⇒ J = A22 − Lc A12
If these conditions are fulfilled, the estimation error dynamics become,
η̃˙ 2 = Jη̃ 2
η̃˙ 2 = (A22 − Lc A12 )η̃ 2 (4.2.8)
As we can see, Lc is the state feedback gain matrix that locates the estimation error
dynamics closed-loop poles at the desired location. The reduced order observer state
equation will be,
ω̇ = Jω + Nη 1 + Mu
ω̇ = (A22 − Lc A12 )ω + (A21 + JLc − Lc A11 )η 1 + (B2 − Lc B1 )u (4.2.9)
If J has eigenvalues with a negative real part, the error dynamics go to zero asymptot-
ically. So we determine Lc such that the closed-loop error dynamics A22 − Lc A12 has
poles at the desired location (i.e., the poles of the reduced-order observer location). In
determining Lc , A22 and A12 take a role of A and C is full order observer respectively.
Chapter 4 114
Daniel Abebe Modern Control Systems Lecture Note
Finally, the combined system state equation (state dynamics) with reduced-order
observer dynamic equation (state equation) is given as
η̇ A11 A12 0 η B1
1 1
η̇ 2 = η 2 + u
A21 A22 0 B2
ω̇ [A12 + JLc − Lc A11 ] 0 [A22 − Lc A12 ] ω [B2 − Lc B1 ]
(4.2.10)
The block diagram representation of the reduced-order state observer is shown in Figure
4.2.
u + 1 y
B1 ʃ
+
+ A11
+
A12
A21
+ +
A22
+
x2
+
B2 ʃ
N
Lc
ω ̂ 2
M + ʃ +
Chapter 4 115
Daniel Abebe Modern Control Systems Lecture Note
Chapter 4 116
Daniel Abebe Modern Control Systems Lecture Note
ẋ = Ax + Bu
y = Cx (4.4.1)
u = −Kf x̂ (4.4.3)
ẋ = Ax − BKf x̂
ẋ = (A − BKf )x + BKf x̃ (4.4.4)
and
This shows that Kf and Lc can be designed independently of each other. If (A, B) is
controllable, Kf can be chosen in a way that the poles of the closed-loop system are at the
desired location, and if (A, C) is observable, Lc can be chosen so that the error dynamics
go to zero rapidly. This property, where the two design problems (observer design problem
and state feedback design problem of a system) can be handled independently, is called
the separation principle. The dynamics of the closed-loop system of full state feedback
with full order state observer is then,
ẋ A − BKf BKf x
=
x̃˙ 0 A − Lc C x̃
Chapter 4 117
Daniel Abebe Modern Control Systems Lecture Note
The block diagram representation of full state feedback with full order state observer
is shown in Figure 4.3.
+ u + x
xud
-
B + ʃ C
y
A
+
Lc
-
+ ʃ
x̂ C
ŷ
A
Kf
Figure 4.3: Block Diagram Representation of Full State Feedback Controller with Full Order
State Observer
u = −K(S1 η 1 + S2 η̂ 2 )
= −K(S2 η 1 + S2 (η 2 − η̃ 2 )) (4.5.1)
and
Chapter 4 118
Daniel Abebe Modern Control Systems Lecture Note
The closed loop dynamics (closed loop state equation) of the system with reduced order
observer (reduced order observer state equation) is
η̇ A − B1 KS1 A12 − B1 KS2 0 B1 KS2 η
1 11 1
η̇ 2 A21 − B2 KS1 A22 − B2 KS2 0
B2 KS2 η 2
= (4.5.5)
˙
η̃ 1 0 0 0 0 η̃ 1
˙η̃ 2 0 0 0 A22 − Lc A12 η̃ 2
The block diagram representation of full state feedback with reduced order state observer
is shown in Figure 4.4.
uxd + u + 1 y
B1 ʃ
- +
+ A11
+
A12
+
Kf S1
+ A21
+ +
A22
+ + x2
B2 ʃ
Lc
̂2
M + ʃ +
S2
Figure 4.4: Block Diagram Representation of Full State Feedback Controller with Reduced Order
State Observer
Chapter 4 119
Daniel Abebe Modern Control Systems Lecture Note
4.6 Examples
1. For a dynamic system
0 1 0
ẋ = x + u
−1 −2 1
h i
y= 0 1 x
Solution
Check for observability of the system
C 0 1
Ov = =
CA −1 −2
s 0 0 1 l1 h i
det − + 0 1 = s2 + 0.75s + 0.125
0 s −1 −2 l2
s −1 + l1
det = s2 + 0.75s + 0.125
1 s + 2 + l2
s2 + (2 + l2 )s + (1 − l1 ) = s2 + 0.75s + 0.125
The state observer gain matrix that puts the poles of the state observer at the
desired location is then,
Chapter 4 120
Daniel Abebe Modern Control Systems Lecture Note
Determine the gain Lc that result in estimation error dynamics with poles at p1 =
−1, p2 = −2, p3 = −3 and closed loop dynamics of the system with full state
feedback using full state observer, if the feedback gain Kf required to place the
closed loop poles of the system at s = −2, s = −2 ± 2i.
Solution
Check for observability of the system
C 1 0 0
O v = CA = 0 1 0
2
CA 0 0 1
Chapter 4 121
Daniel Abebe Modern Control Systems Lecture Note
s 0 0 0 1 0 l
1 h i
3 2
det 0 s 0 − 0 1 + l2 1 0 0 = s + 6s + 11s + 6
0
0 0 s −4 −4 −1 l3
s + l1 −1 0
3 2
det l2 −1 = s + 6s + 11s + 6
s
l3 + 4 4 s+1
The state observer gain matrix that puts the poles of the state observer at the
desired location is then,
1 + l1 = 6 ⇒ l1 = 5,
l1 + l2 + 4 = 11 ⇒ l2 = 2
4l1 + l2 + l3 + 4 = 6 ⇒ l3 = −20
l 5
1
Lc = l2 = 2
l3 −20
Chapter 4 122
Daniel Abebe Modern Control Systems Lecture Note
The closed-loop dynamics with a full-order state observer based on state feedback
control are
ẋ A − BK BKf x
=
x̃˙ 0 A − Lc C x̃
Chapter 4 123
Daniel Abebe Modern Control Systems Lecture Note
Design a reduce order state observer for the system that has poles at p1 = −5, p2 =
−5 and closed loop dynamics of the system with full state feedback using reduced
order state observer, if the feedback gain Kf required to place the closed loop poles
of the system at s = −2, s = −2 ± 2i.
Solution
Decomposing the system into a state that can be measured directly from the output
and that can’t be measured directly from the output,
..
0 . 1 0 1
. . . . . . . . . . . . . . .
A=
..
, B =
0 . 0 1 0
..
−4 . −4 −1 0
.
h i
C = 1 .. 0 0
We have
h i
A11 = 0, A12 = 1 0 B1 = 1 C1 = 1
0 0 1 0 h i
A21 = A22 = B2 = C2 = 0 0
−4 −4 −1 0
As we can see from the decomposition of the system, x2 and x3 can’t be measured
directly, while x1 can be measured directly. So the observer is designed for x2 and
x3 .
Chapter 4 124
Daniel Abebe Modern Control Systems Lecture Note
s2 + (1 + l1 )s + (l1 + l2 + 4) = s2 + 10s + 25
The state observer gain matrix that puts the poles of the state observer at the
desired location is then,
1 + l1 = 10 ⇒ L1 = 9, l1 + l2 + 4 = 25 ⇒ l2 = 12
l1 9
Lc = =
l2 12
The matrices
−9 1
J = A22 − Lc A12 =
−16 −1
−69
N = A21 + JLc + Lc A11 =
−160
0
M = B2 − Lc B1 =
10
Chapter 4 125
Daniel Abebe Modern Control Systems Lecture Note
The combined open-loop system dynamics of the system and reduced-order observer
dynamics are given as,
..
ẋ 0 1 0 . 0 0 x 0
1 1
..
ẋ2 0 0 1 . 0 0 x2 0
..
ẋ3 −4 −4 −1 .
0 0 x3 1
= + u
. . . . . . . . . . . . . . . . . . . . . . . . . . .
..
ω̇1 −69 0 −9 1 ω1 0
0 .
..
ω̇2 −160 0 0 . −16 −1 ω2 1
h h i i
A12 − B1 Kf S2 = 1 0 − 12 12 5 S2
⊤
1 1 1
which imply S2 is 3 × 2 matrix, let S2 = we have,
1 1 1
1 1
h i h i h i h i h i
A12 − B1 KS2 = 1 0 − 12 12 5 1 1 = 1 0 − 29 29 = −28 −29
1 1
1 1
h i h i
B1 KS2 = 12 12 5 1 1 = 29 29
1 1
1
0 0 h i 0
A21 − B2 KS1 = − 12 12 5 1 =
−4 0 −4
1
Chapter 4 126
Daniel Abebe Modern Control Systems Lecture Note
1 1
0 1 0 h i 0 1
A22 − B2 KS2 = − 12 12 5 1 1 =
−4 −1 0 −4 −1
1 1
1 1
0 h i 0 0
B2 KS2 = 12 12 5 1 1 =
0 0 0
1 1
0 1 9 h i −9 1
A22 − Lc A12 = − 1 0 =
−4 −1 12 −16 −1
Resulting in the closed-loop system dynamics with reduced state observer after
h i⊤
replacing η1 = x1 (directly measured states)and η 2 = x2 x3 (estimated states)
we have,
.. .. ..
−29 . −28 −29 . 0 . 29 29
ẋ1 . . . . . . . . . . . . . . . . . . . . . . . . . . . x1
.. .. ..
ẋ2 0 . 0 1 . 0 . 0 0 x2
.. .. ..
ẋ3 −4 . −4 −1 . 0 . 0 0 x3
. . . = . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
.. .. ..
˙
x̃1 0 . 0 0 . 0 . 0 0 x̃1
˙
x̃2 . . . . . . . . . . . . . . . . . . . . . . . . . . . x̃2
.. .. ..
˙x̃3
0 . 0 0 . 0 . −9 1 x̃3
.. .. ..
0 . 0 0 . 0 . −16 −1
Design a reduce order state observer for the system that has poles at p1 = −5, p2 =
−5 and closed loop dynamics of the system with full state feedback using reduced
Chapter 4 127
Daniel Abebe Modern Control Systems Lecture Note
order state observer, if the feedback gain Kf required to place the closed loop poles
of the system at s = −2, s = −2 ± 2i.
Solution
Check the controllability of the system
0 0 1
0 1 −1
1 −1 −3
..
0 . 0 1 0
..
η̇ . . . . . . . . . . η1 . . .
1 = + u
.
1 .. 0
η̇ 2 0 η2 0
..
−4 . −4 −1 1
i η1
.
h
y = 1 .. 0 0
η2
Chapter 4 128
Daniel Abebe Modern Control Systems Lecture Note
h i
A11 = 0 A12 = 0 1 B1 = 0 C 1 = 1
1 0 0 0 h i
A21 = A22 = B2 = C2 = 0 0
−4 −4 −1 1
we have
The value of the reduced observer matrix gain Lc that places the poles of the
estimation reduced-order state observer error dynamics at the desired location is,
∥sI − J∥ = s2 + 10s + 25
∥sI − A22 + Lc A12 ∥ = s2 + 10s + 25
s 0 0 0 l h i
− + 1 0 1 = s2 + 10s + 25
0 s −4 −1 l2
s l1
= s2 + 10s + 25
4 s + 1 + l2
s2 + (l2 + 1)s − 4l1 = s2 + 10s + 25
−6.25
l2 + 1 = 10 ⇒ l2 = 9 − 4l1 = 25 ⇒ l1 = −6.25 Lc =
9
then
0 6.25
J = A22 − Lc A12 =
−4 −10
57.25
N = A21 − Lc A11 + JLc =
−69
0
M = B2 − Lc B1 =
1
Chapter 4 129
Daniel Abebe Modern Control Systems Lecture Note
The open-loop system dynamics combined with minimum state observer dynamics
are
..
ẋ 0 1 0 . 0 0 x 0
1 1
..
ẋ2 0 0 1 . 0 0 x2 1
..
ẋ3 −4 −4 −1
. 0 0 x3 0
= + u
. . . . . . . . . . . . . . . . . . . . . . . . . . .
.
0 57.25 ..
ω̇1 0 0 6.25 ω1 0
..
ω̇2 0 0 −69 . −4 −10 ω2 1
h i
Since matrix B is not changed Kf = 12 12 5 , the closed loop dynamics of the
system with reduced order observer are determined as,
h i
A11 − B1 Kf S1 = 0 − 12 12 5 S1
h i⊤
which imply S1 is 3 × 1 matrix, let S1 = 1 1 1 , with this
1
h i
A11 − B1 Kf S1 = 0 − [0] 12 12 5 1 = 0
1
h i h i
A12 − B1 Kf S2 = 0 1 − [0] 12 12 5 S2
⊤
1 1 1
which imply S2 is 3 × 2 matrix, let S2 = we have,
1 1 1
1 1
h i h i h i
A12 − B1 Kf S2 = 0 1 − [0] 12 12 5 1 1 = 0 1
1 1
1 1
h i h i
B1 Kf S2 = [0] 12 12 5 1 1 = 0 0
1 1
Chapter 4 130
Daniel Abebe Modern Control Systems Lecture Note
1
1 0 h i 1
A21 − B2 Kf S1 = − 12 12 5 1 =
−4 1 −33
1
1 1
0 0 0 h i 0 0
A22 − B2 Kf S2 = − 12 12 5 1 1 =
−4 −1 1 −33 −30
1 1
1 1
0 h i 0 0
B2 Kf S2 = 12 12 5 1 1 =
1 29 29
1 1
0 0 −6.25 h i 0 6.25
A22 − Lc A12 = − 0 1 =
−4 −1 9 −4 −10
Finally, the closed-loop system dynamics with reduced state observer after replacing
h i⊤
η1 = x2 (directly measured states)and η 2 = x1 x3 (unmeasured states) is,
.. .. ..
0 . 0 1 . 0 . 0 0
ẋ2 . . . . . . . . . . . . . . . . . . . . . . . . . . . x2
.. .. ..
ẋ1 1 . 0 0 . 0 . 0 0 x1
.. .. ..
ẋ3 −33 . −33 −30 . 0 . 29 29 x3
. . . = . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
.. .. ..
˙
x̃2 0 . 0 0 . 0 . 0 0 x̃2
˙
x̃1 . . . . . . . . . . . . . . . . . . . . . . . . . . . x̃1
.. .. ..
˙x̃3
0 . 0 0 . 0 .
0 6.25 x̃3
.. .. ..
0 . 0 0 . 0 . −4 −10
Chapter 4 131
Daniel Abebe Modern Control Systems Lecture Note
Design a reduce order state observer for the system that has poles at p1 = −5, p2 =
−5 and closed loop dynamics of the system with full state feedback using reduced
order state observer, if the feedback gain Kf required to place the closed loop poles
of the system at s = −2, s = −2 ± 2i.
Solution
Check the controllability of the system
0 0 1
0 1 −1
1 −1 −3
Chapter 4 132
Daniel Abebe Modern Control Systems Lecture Note
h i
A11 = −1 A12 = −4 −4 B1 = 1 C1 = 1
0 1 0 0 h i
A21 = A22 = B2 = C2 = 0 0
1 0 0 0
So we have
The value of the reduced observer matrix gain Lc that place the poles of estimation
reduced order state observer error dynamics at the desired location is
∥sI − J∥ = s2 + 10s + 25
∥sI − A22 + Lc A12 ∥ = s2 + 10s + 25
s 0 1 0 l h i
− + 1 −4 −4 = s2 + 10s + 25
0 s 0 0 l2
s + 1 − 4l1 −4l1
= s2 + 10s + 25
−4l2 s − 4l2
s2 + (1 − 4l1 − 4l2 )s − 4l2 = s2 + 10s + 25
4
1 − 4l1 − 4l2 = 10 ⇒ l2 = 4 − 4l2 = 25 ⇒ l2 = −6.25 Lc =
−6.25
then
17 16
J = A22 − Lc A12 =
−25 −25
−28
N = A21 − Lc A11 + JLc =
49
−4
M = B2 − Lc B1 =
6.25
Chapter 4 133
Daniel Abebe Modern Control Systems Lecture Note
The open-loop system dynamics combined with minimum state observer dynamics
are
..
ẋ1 0 1 0 . 0 0 x 0
1
..
ẋ2 0 0 1 . 0 0 x2 1
..
ẋ3 −4 −4 −1
. 0 0 x3 0
= + u
. . . . . . . . . . . . . . . . . . . . . . . . . . .
.
0 −28 ..
ω̇1 0 17 16 ω1 0
..
ω̇2 0 0 49 . −25 −25 ω2 1
h i
Since matrix B is not changed Kf = 12 12 5 , the closed loop dynamics of the
system with reduced order observer are
h i⊤
by letting S1 = 1 1 1 , with this
h i
B1 KS2 = 29 29
0
A21 − B2 KS1 =
−1
1 0
A22 − B2 KS2 =
0 0
0 0
B2 KS2 =
0 0
Chapter 4 134
Daniel Abebe Modern Control Systems Lecture Note
17 16
A22 − Lc A12 =
−25 −25
h i⊤
Replacing η1 = x3 (directly measured states)and η 2 = x1 x2 (unmeasured
states),, the closed-d-loop system with reduced state observer is,
.. .. ..
−30 . −33 −33 . 0 . 29 29
ẋ3 . . . . . . . . . . . . . . . . . . . . . . . . . . . x3
.. .. ..
ẋ1 0 . 1 0 . 0 . 0 0 x1
.. .. ..
ẋ2 −1 . 0 0 . 0 . 0 0 x2
. . . = . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
.. .. ..
˙
x̃3 0 . 0 0 . 0 . 0 0 x̃3
˙
x̃1 . . . . . . . . . . . . . . . . . . . . . . . . . . . x̃1
.. .. ..
˙x̃2
0 . 0 0 . 0 . 17
16 x̃2
.. .. ..
0 . 0 0 . 0 . −25 −25
Solution
The observability matrix of the system is
1 0 0
0 1 0
0 0 1
which has full rank, implying that the system is completely observable. But from
the system output matrix, only one state can be directly measured. The two states
are not available. Hence, the reduced order observer has order 2. If the observer is
required to give a pole at p = −2 ± 2i, then we have
..
0 . 1 0 0
. . . . . . . . . . . . . . .
ẋ =
..
x + u
0 . 0 1 1
..
−6 . −11 −6 −3
.
h i
y = 1 .. 0 0
Chapter 4 135
Daniel Abebe Modern Control Systems Lecture Note
Implying that
0 1 0
h i
A22 = A11 = 0 A12 = 1 0
A21 =
−11 −6 −6
1 h i
B1 = 0 B2 = C1 = 1 C2 = 0 0
−3
Finally
0 13 13
N = A21 + JLc + Lc A11 = + =
−6 4 −2
1
M = B2 − Lc B1 =
−3
Chapter 4 136
Daniel Abebe Modern Control Systems Lecture Note
The open-loop dynamics of the system, combined with reduced-order state observer
dynamics, are
..
ẋ1 0 1 0 . 0 0 x 0
1
..
ẋ2 0 0 1 . 0 0 x2 1
..
ẋ3 −6 −11 −6 . 0 x3 −3
0
= + u
. . . . . . . . . . . . . . . . . . . . . . . . . . .
..
−2 1 ω1 1
ω̇1 13 0 0 .
..
ω̇2 −2 0 0 . −20 −6 ω2 −3
Solution
As shown in chapter 2 on example 2, the system has one unobservable mode. So
we cannot design an observer for the unobservable mode of the system, but we can
design for the observable modes of the system. As we can see from the output
matrix of the decomposed system, one state can be directly measured, while one
cannot. Thus, we design a state observer of order 1. The dynamics look of the
observable mode of the system is
..
2 . 1 1
ẋ = . . . . . . . . . x + . . . u
..
−4 . 0 1
.
h i
y = 1 .. 0 x
From this A11 = 2, A12 = 1, A21 = −4 A22 = 0. If the observer is required to have
a pole at p = −1, then
sI − A22 + lA12 = s + 1
s − 0 + l = s + 1 ⇒ l = 1 ⇒ Lc = l = 1
J = A22 − Ac A12 = 0 − 1 = −1
N = A21 + JLc + Lc A11 = −4 − 1 + 2 = −3
M = B2 − Lc B1 = 1 − 1 = 0
ω̇ = −ω − 3x1
Chapter 4 137
Daniel Abebe Modern Control Systems Lecture Note
The open-loop system state equation, combined with the reduced-order state ob-
server equation, is
..
ẋ1 2 1 . 0 x 1
1
..
ẋ2 −4 0
. 0 x2 1
= + u
. . . . . . . . . . . . . . . . . . . . .
..
ω̇ −3 0 . 1 ω 0
4.7 Exercises
1. Consider a system
0.16 2.16 −1 h i
ẋ = x+ u y= 1 1 x
−0.16 −1.16 1
2. Design a full state observer for IPC dynamics on 1.6.3 that has poles at −1, −1, −1±
i. Use any valid parameter.
3. Design a full state observer for DC motor dynamics having poles at −1 ± 2i. Use
any valid parameters.
4. Design a reduced-order state observer for the dynamics of the DC motor on 1.7.1,
where the speed is directly measured and the armature current ia is required to be
estimated. The pole of the observer must be at −1, use any valid parameters.
5. Design a reduced state observer for IPC dynamics on 1.6.3, where the position of the
cart and position of the pendulum are directly measured and the angular velocity
of the pendulum and linear velocity of the cart are not. The poles of the observer
are at−1 ± 2i.
6. Design a Kalman-Bucy filter for IPC dynamics and DC motor dynamics. Use any
valid parameters.
Chapter 4 138
Daniel Abebe Modern Control Systems Lecture Note
Chapter 4 139
Chapter 5
Optimal control deals with the most celebrated problem of modern control theory, ”deter-
mination of best control strategy, selecting optimum control strategy, selecting optimum
control vector u(t), that minimizes cost function or performance index”. Consider an
LTI system with initial condition x0 :
ẋ = Ax + Bu
y = Cx + Du (5.0.1)
The objective of optimal control is to determine the control vector u that forces the
behavior of the dynamic system to the required final state by maximizing the required
cost function and at the same time fulfilling the physical constraints.
The cost function has a general representation of:
Z tf
t=tf
J = θ[x(t), t] |t=t0 + φ(x(t), u(t), t)dt (5.0.2)
t0
t=t
The first term θ[x(t), t] |t=tf0 is called terminal error optimization problem (terminal cost
Z tf
function) and the second term φ(x(t), u(t), t)dt is called optimal cost function over
t0
the entire trajectory (the function to be optimized over the entire trajectory).
140
Daniel Abebe Modern Control Systems Lecture Note
Chapter 5 141
Daniel Abebe Modern Control Systems Lecture Note
∆f depends both on t and ∆t, hence it is written as ∆f (t, ∆t). In similar way for a
functional J(x(t), t) the incremental of a functional is defined as:
δx(t) is called variation of the function x(t). Since the increment of a functional depends
on x(t) and δx(t) it is written as J((x(t), δx(t)).
1 ∂ 2J 2
∗ ∂J
∆J = J(x (t)) + δx + 2
δx + . . . − J(x∗ (t)) (5.2.5)
∂x ∗
x (t) 2! ∂x ∗
x (t)
∂J 1 ∂ 2J 2
δJ = δx is called the first variation, and δ 2 J = δx is called the second variation
∂x ∗
2! ∂x2
of the functional about a function x (t).
Chapter 5 142
Daniel Abebe Modern Control Systems Lecture Note
∂J
∆J ≈ δJ = δx (5.2.6)
∂x
Similar conditions hold for functional to have optimum value, the necessary condition
is that its first variation must vanish (∆J = 0) at x∗ (t) and the sufficient condition is
that ∆2 J < 0 for local maximum or ∆2 J > 0 for local minimum, for all x(t), such that
(x(t) − x∗ (t) < ϵ, for ϵ > 0).
A functional may be a function of one or more functions, like a function of several variables
V (x, y, z). In a similar manner, we have a functional which can be J(x(t), y(t), z(t)). The
only difference is in function; the variables x, y, and z are independent variables, while
in functional, they are a function of another independent variable.
Chapter 5 143
Daniel Abebe Modern Control Systems Lecture Note
∂φ
By letting u = and dv = d(δx),
∂ ẋ
Z tf Z tf
∂φ ∂φ tf d ∂φ
d (δx) = δx |t0 − δx dt
t0 ∂ ẋ ∂ ẋ t0 dt ∂ ẋ
Z tf Z tf
∂φ ∂φ tf d ∂φ
d (δx) = δx |t0 − δx dt (5.2.14)
t0 ∂ ẋ ∂ ẋ t0 dt ∂ ẋ
Replacing equation 5.2.14 in equation 5.2.13 for the second part and rearranging we have,
Z tf
∂φ d ∂φ ∂φ t
− δx dt + δx |tf0 = 0 (5.2.15)
t0 ∂x dt ∂ ẋ ∂ ẋ
Chapter 5 144
Daniel Abebe Modern Control Systems Lecture Note
Equation A.4.5 is called the Euler-Lagrange equation, while equation A.4.6 is called the
∂φ
boundary condition equation. For the boundary condition to be zero, either must be
∂ ẋ
zero or δx must be zero at the boundary. This leads to four cases:
1. Case-I: x is fixed at t0 and tf (fixed end trajectory).
Here, x(t0 ) = C1 and x(tf ) = C2 , where C1 and C2 are constants. Since δx = 0 at
∂φ
both t0 and tf , there are no restrictions on .
∂ ẋ
2. Case-II: x is fixed at t0 and free at tf .
∂φ
Here, x(t0 ) = C1 , and = 0 at tf because δx is unknown at tf . No restrictions
∂ ẋ
∂φ
apply to at t0 .
∂ ẋ
3. Case-III: x is fixed at tf and free at t0 .
∂φ
Here, x(tf ) = C2 , and = 0 at t0 because δx is unknown at t0 . No restrictions
∂ ẋ
∂φ
apply to at tf .
∂ ẋ
4. Case-IV: x is free at t0 and tf .
∂φ
Here, = 0 at both t0 and tf because δx is unknown at both times.
∂ ẋ
The sufficient condition for the functional to have optimum value is derived using δ 2 J
(second variation),
Z tf 2
∂ 2φ ∂ 2φ
2 1 ∂ φ 2 2
δ J= (δx) + (δ ẋ) + 2 δ ẋδx dt (5.2.18)
t0 2 ∂x2 ∂ ẋ ∂ ẋ∂x
∂ 2φ
Solving the last term using integration by part u = δ ẋδx and dv = δ ẋ and rearrang-
∂x
ing,
tf
∂ 2φ ∂ 2φ ∂ 2φ
Z
2 1 d 2 2
δ J= 2
− (δx) + (δ ẋ) dt (5.2.19)
2 t0 ∂x dt ∂ ẋ∂x ∂ ẋ2
which means
∂ 2φ ∂ 2φ ∂ 2φ
d
2
− ≥0 and ≥0
∂x dt ∂x∂ ẋ ∂ ẋ2
for local minimum and
∂ 2φ ∂ 2φ ∂ 2φ
d
2
− ≤ 0 and ≤0
∂x dt ∂x∂ ẋ ∂ ẋ2
for local maximum.
Chapter 5 145
Daniel Abebe Modern Control Systems Lecture Note
subjected to
ẋ = Ax + Bu (5.2.21)
Let us define a general function that includes the constraint and the function to be
optimized as
Chapter 5 146
Daniel Abebe Modern Control Systems Lecture Note
Matrix Q is a weighing matrix for state error, and R is a weighing matrix for input.
The Choice depended on the optimization problem and design specifications. If state
error is required to be minimized matrix Q should be high, and if input to the system is
required to be minimized more than the state tracking error matrix R should be high.
In state feedback controller design, penalizing input is the same as penalizing state, since
the control signal is a linear combination of the state of the system.
The possible starting of guessing the values of the matrices can be one of the following
• Simplest choice: Q = I, R = ρI ⇒ L = x⊤ · x + ρu⊤ u then vary ρ to get something
that has good response.
• Diagonal weights
q 0 0 0
1
0 q2 0 0
Q=
0 0 q3 0
0 0 0 q4
1
Choose each qi = , where x̃i is the tolerable error of each state of the system.
x̃2i
• Output weighting
Let z = Hx be the output you want to keep small. Assume (A, H) observable. Use
Q = x⊤ x, R = ρI ⇒ trade off z⊤ · z vs ρu2 .
Where H = x⊤ Qx + u⊤ Ru + λ⊤ and λ(t) = P(t)x(t)
• Trial and error (on weights)
Most of the time Q = C⊤ C and the weight of the matrix is adjusted based on the per-
formance requirement of the system.
Our objective is to determine a control law u(t) that optimizes the optimal linear quadratic
regulator cost function subjected to ẋ = Ax + Bu. To do so, we use Euler-Lagrangian
equations:
∂φ d ∂φ ⊤
− = 0 ⇒ x⊤ Q + x⊤ Q⊤ + λ⊤ A − λ̇ = 0 (5.3.2)
∂x dt ∂ ẋ
∂φ d ∂φ
− = 0 ⇒ u⊤ Q + u⊤ Q⊤ + λ⊤ B = 0 (5.3.3)
∂u dt ∂ u̇
∂φ d ∂φ
− = 0 ⇒ ẋ = Ax + Bu (5.3.4)
∂λ⊤ dt ∂ λ̇⊤
Chapter 5 147
Daniel Abebe Modern Control Systems Lecture Note
⊤
Let λ = Px, where P is positive definite symmetric matrix, So λ⊤ = x⊤ P and λ̇ =
Q + Q⊤ R + R⊤
ẋ⊤ P = (Ax + Bu)⊤ P = (x⊤ A⊤ + u⊤ B⊤ )P, replacing for Q = ,R= ,
2 2
⊤
λ⊤ and λ̇ , the Euler Lagrangian equation become:
x⊤ Q + x⊤ PA + (x⊤ A⊤ + u⊤ B⊤ )P = 0
x⊤ (Q + pA + A⊤ P) + u⊤ B⊤ P = 0 (5.3.5)
⊤ ⊤
u R + x PB = 0 (5.3.6)
ẋ = Ax + Bu (5.3.7)
If the control law is design based on state feedback u = −Kf x, optimizing u(t) is sames
as determining the optimal state feedback gain matrix Kf , hence with state feedback
control law, equation 5.3.5 becomes
x⊤ (Q + PA + A⊤ P)x⊤ K⊤ ⊤
fB P = 0
x⊤ (Q + PA + A⊤ P − K⊤ B⊤ P) = 0
Q + PA + A⊤ P − K⊤ ⊤
fB P = 0 (5.3.8)
− x⊤ K⊤ ⊤
f R + x PB = 0
x⊤ (−K⊤
f R + PB) = 0
K = R−1 B⊤ P (5.3.9)
The control law u(t) that optimizes the functional (cost function) is,
K is the optimal feedback gain matrix that optimizes u(t), which in turn optimizes our
optimal linear quadratic regulator cost function, subjected to ẋ = Ax + Bu.
Replacing Kf in equation 5.3.8 results in a Riccati equation in which we solve for the
matrix P.
Q + PA + A⊤ P − PBR−1 B⊤ P = 0 (5.3.11)
The other method of finding solution for linear quadratic optimal regulator without using
d ⊤
Euler-Lagrangian equations achieved is by letting x⊤ (Q+K⊤ f RKf )x = − (x Px) after
dt
replacing for u in equation 5.1.7 and final state xf = 0 (same control objective bring the
system to zero final state). which results in
x⊤ (Q + K⊤ ⊤ ⊤
f RKf )x = −ẋ Px − x Pẋ
x⊤ (Q + K⊤ ⊤ ⊤ ⊤ ⊤ ⊤
f RKf )x = −x (A − Kf B )Px − x P(A − BKf )x
A⊤ P − K⊤ ⊤ ⊤
f B P + PA − PBKf = −(Q + Kf RKf ) (5.3.12)
Chapter 5 148
Daniel Abebe Modern Control Systems Lecture Note
Finding u that minimize J is same as finding Kf that minimize J, which in-turn minimizes
equation 5.3.11. Derivative of equation 5.3.11 with respect to Kf to get Kf that minimize
the equation result in:
Kf = R−1 B⊤ P (5.3.13)
So the control signal that minimizes the optimal quadratic regulator function is
ẋ = Ax + Bu ⇒ ẋ = (A − BKf )x
ud + u + x y
- B ʃ C
+
Kf
Figure 5.1: Block Diagram Representation of LTI System with Optimal State Feedback Controller
Chapter 5 149
Daniel Abebe Modern Control Systems Lecture Note
ẋ = Ax + Bu
y = Cx + Du (5.4.1)
At this point, the rate of change of the state variables will be zero, and the output is
expected to be equal to the reference values. Hence for a pair (xd ,ud ) equation 5.4.1 will
be
Axd + Bxd = 0
r = Cxd + Dud (5.4.2)
(5.4.3)
−A B −x 0
d = (5.4.4)
−C D ud r
• Case-I: When the number of inputs m, is strictly smaller than the number of
controlled outputs p, we have an under-actuated system. In this case, the system of
equations 5.4.4 generally does not have a solution because it presents more equations
than unknowns.
• Case-II: When the number of inputs m is equal to the number of controlled outputs
p, equation 5.4.4 always has a solution as long as the Rosenbrock’s 1 system matrix
sI − A B
P (s) = (5.4.6)
−C D
1
Rosenbrock’s matrix is obtained by regarding the controller output as the only output of the system.
Chapter 5 150
Daniel Abebe Modern Control Systems Lecture Note
• Case-III: When the number of inputs m is strictly larger than the number of
controlled outputs p, we have an over-actuated system and equation 5.4.4 generally
has multiple solutions.
When P(0) is full row-rank, i.e., when it has n + p linearly independent rows, the
(n + p) × (n + p) matrix P(0)P(0)⊤ is non-singular and one solution to equation
5.4.4 is given by
x 0
d = P(0)⊤ (P(0)P(0)⊤ )−1 (5.4.7)
ud r
Also in this case, s = 0 should not be an invariant zero of the system because
otherwise P(0) cannot be full rank.
The optimal set-point problem can be reduced to that of optimal regulation by considering
an auxiliary system with state x̃ = x − xd , ũ = u − ud whose dynamics are:
We can then regard equation 5.4.5 and equation 5.4.8 as an optimal regulation problem
for which the optimal solution is given by:
ũ = −Kf x̃ (5.4.9)
Going back to the original input and state variables u and x, we conclude that the optimal
control for the set-point problem defined by equation 5.4.1 and equation 5.4.5 is given by
u = −Kf (x − xd ) + ud (5.4.10)
xd = Fr (5.4.11)
ud = Nr (5.4.12)
The values of F and N are determined using equation 5.4.2, by replacing for xe = Fr and
ud = Nr:
Solving for F from equation 5.4.13 and replacing in equation 5.4.14, the values for matrices
N and F is:
Chapter 5 151
Daniel Abebe Modern Control Systems Lecture Note
ud
N
r xd + u + x y
F + -
Kf + B ʃ C
+
u = −Kf x + Ko r (5.4.17)
as time goes to infinity (t → ∞) the state x, goes to the desired state xd and the input
u, tends toward desired input ud , having dynamic equation:
since ẋd is zero at t → ∞, then the dynamic equation 5.4.19, can be written as
Form equation 5.4.24, by replacing for ud and xd ,the value of matrix Ko is:
In cases where the states of the dynamic system are not directly accessible, a state
observer can be utilized to estimate the state of the system, thereby allowing the use of
the estimated state for control purposes.
Chapter 5 152
Daniel Abebe Modern Control Systems Lecture Note
r u + x y
Ko
+ B ʃ C
- +
Kf
The LQR approach is a widely used method for designing state feedback controllers. It
involves minimizing a quadratic cost function that captures the system’s performance and
control effort. The optimal state feedback controller gain matrix is obtained by solving
the associated algebraic Riccati equation. Once the controller gain matrix is obtained, it
can be used to stabilize the system and achieve the desired performance.
On the other hand, the problem of designing an optimal state estimator gain matrix in-
volves minimizing a similar quadratic cost function that captures the difference between
the actual state and the estimated state. The optimal estimator gain matrix is obtained
by solving another algebraic Riccati equation. Once the estimator gain matrix is ob-
tained, it can be used to estimate the system’s state accurately.
The duality property between these two problems arises because the same cost function
is used in both cases, but with different weighting matrices. In other words, the solution
to one problem can be used to solve the other problem by simply swapping the weighting
matrices.
ẋ = Ax + Bu
y = Cx + Du (5.5.1)
where x is the state vector, u is the control input, y is the output, and A, B, C, and D
are matrices that describe the system’s dynamics.
Chapter 5 153
Daniel Abebe Modern Control Systems Lecture Note
To design an optimal state feedback controller gain matrix using LQR, we first define
the cost function as follows:
Z tf
J= ((x̂ − x)⊤ Q(x̂ − x) + u⊤ Ru) dt (5.5.2)
t0
where Q and R are weighting matrices that capture the relative importance of the state
and control effort, respectively. The optimal controller gain matrix Kf is obtained by
solving the following algebraic Riccati equation:
Q + PA + A⊤ P − PBR−1 B⊤ P = 0 (5.5.3)
To design an optimal state estimator gain matrix using LQR, we define a similar cost
function as follows:
Z tf
J= ((x̂ − x)⊤ Q(x̂ − x) + (y⊤ Ry y)) dt (5.5.4)
t0
where Q and Ry are weighting matrices that capture the relative importance of the
output and state estimation error, respectively. The optimal estimator gain matrix Lc is
obtained by solving the following algebraic Riccati equation:
Q + PA + A⊤ P − PC⊤ R−1
y CP = 0 (5.5.5)
The duality property between the two problems arises from the fact that they share
the same mathematical structure. Specifically, the state feedback controller gain matrix
design problem can be formulated as a quadratic optimization problem using the linear
quadratic regulator (LQR) technique. Similarly, the state estimator gain matrix design
problem can be formulated as a quadratic optimization problem using the linear quadratic
estimator (LQE) technique.
The LQR and LQE techniques are based on minimizing a cost function that measures
the deviation of the system’s states from a desired trajectory and the deviation of the
estimated states from the actual states, respectively. The cost function is formulated as
a quadratic function of the system’s states and inputs, and its coefficients are determined
by the system’s dynamics and performance specifications.
The duality property between the two problems implies that the optimal gain matrices
obtained from the LQR and LQE techniques are related to each other. Specifically, the
optimal state feedback controller gain matrix is equal to the transpose of the optimal
state estimator gain matrix, up to a scaling factor. This means that if we solve the state
feedback controller gain matrix design problem using LQR, we can obtain the optimal
state estimator gain matrix by transposing the solution and scaling it appropriately.
Chapter 5 154
Daniel Abebe Modern Control Systems Lecture Note
5.6 Examples
1. Determine the increment of the following functional
Z tf
J= (2x2 (t) + 1)dt (5.6.1)
t0
Solution
Solution
∂J
δJ = δx(t)
Z∂xtf
∂
= (4x(t)δx(t) + 2(δx(t))2 )δx(t) dt
t ∂x
Z 0tf
= (4δx(t) + 4δx(t))δx(t) dt
t0
Z tf
= 8(δx(t))2 dt
t0
π
with boundary condition x(0) = 0, x( ) = 1
2
Solution:
From the cost function we know that φ(x, ẋ) = x2 − ẋ2 , from Euler Lagrangian
equations
∂φ d ∂φ
− =0
∂x dt ∂ ẋ
Chapter 5 155
Daniel Abebe Modern Control Systems Lecture Note
ẍ + x = 0
λ2 + λ = 0 ⇒ λ1 = −i, λ2 = i
x = C1 eλ1 t + C2 eλ2 t
= C1 e−it + C2 eit
= (C1 + C2 ) cos t + i(C2 − C1 ) sin t
Solving the two equations simultaneously results in C1 = 0.5 and C2 = −0.5. Hence,
the value of x that minimizes the cost function is
x(t) = sin(t)
subjected to
0 1 0
ẋ = x + u
0 0 1
with initial and final states x(0) = [1 1]⊤ and x(1) = [0 0]⊤ .
Chapter 5 156
Daniel Abebe Modern Control Systems Lecture Note
Solution:
1
From constrained optimization we know that φ(x, ẋ, u) = u2 + λ⊤ (Ax + Bu − ẋ).
2
The control law we determine should minimize
Z 1
1 2
J= u + λ1 (x2 − ẋ1 ) + λ2 (u − ẋ2 ) dt
0 2
Using the Euler-Lagrange equations
∂φ d ∂φ
− = 0 ⇒ 0 + λ̇1 = 0 ⇒ λ̇1 = 0
∂x1 dt ∂ ẋ1
∂φ d ∂φ
− = 0 ⇒ x2 − ẋ1 = 0 ⇒ ẋ1 = x2
∂λ1 dt ∂ λ̇1
result in the state equation.
∂φ d ∂φ
− = 0 ⇒ λ1 + λ̇2 = 0 ⇒ λ̇2 = −λ1
∂x2 dt ∂ ẋ2
∂φ d ∂φ
− = 0 ⇒ u − ẋ2 = 0 ⇒ ẋ2 = u
∂λ2 dt ∂ λ̇2
result in the state equation.
∂φ d ∂φ
− = 0 ⇒ u − λ2 = 0 ⇒ λ2 = −u
∂u dt ∂ u̇
λ̇1 = 0 ⇒ λ1 = C1 ,
λ̇2 = −λ1 ⇒ λ2 = −C1 t + C2
u(t) = −λ2 ⇒ u(t) = C1 t + C2
solving this equations C1 = 18 and C2 = 10. Hence, the control law that optimizes
the constrained cost function is
u(t) = 18t − 10
Chapter 5 157
Daniel Abebe Modern Control Systems Lecture Note
5. For the following system, determine an optimal full state feedback gain matrix Kf .
−1 0 0
ẋ = x + u
0 −2 1
h i
y= 1 0 x
Solution
From Riccati equation
Q + PA + A⊤ P − PBR−1 B⊤ P = 0
Let
p11 p12 1 0
P= Q = C⊤
1 C1 =
R = 0.5;
p21 p22 0 0
1 0 p p −1 0 −1 0 p p
+ 11 12 + 11 12
0 0 p21 p22 0 −2 0 −2 p21 p22
p11 p12 0 h i p p
− 2 0 1 11 12 = 0
p21 p22 1 p21 p22
1 0 −p −2p12 −p11 −p12 p p p p
+ 11 + − 2 12 21 12 22 = 0
0 0 −p21 −2p22 −2p21 −2p22 p22 p21 p22 p22
1 − 2p11 − 2p12 p21 −3p12 − 2p12 p22 0 0
=
2
−3p12 − 2p12 p22 −4p22 − 2p22 0 0
Solving for P
− 4p22 − 2p222 = 0 ⇒ p22 = 0 or p22 = −2
− 3p12 − 2p12 p22 = 0
1 − 2p11 − 2p212 = 0
Solving for p12 and p11
(a) Case I p22 = 0
− 3p12 − 2p12 p22 = 0 ⇒ p12 (3 + 2p22 ) = 0 ⇒ p12 = 0
1 − 2p11 − 2p212 = 0 ⇒ 1 − 2p11 = 0 ⇒ p11 = 0.5
Chapter 5 158
Daniel Abebe Modern Control Systems Lecture Note
So we have
0.5 0 0.5 0
P1 = P2 =
0 0 0 −0.8
The control law becomes u = −Kf x = 0, with this control law the closed loop
system becomes
ẋ = Ax + Bu ⇒ ẋ = (A − BKf )x
∥λI − (A − BKf )∥ = 0
λ 0 −1 0
− =0
0 λ 0 −2
λ+1 0
=0
0 λ+2
(λ + 1)(λ + 2) = 0 ⇒ λ1 = −1 λ2 = −2
If the error x̃1 = 0.01 and x̃2 = 0.1 is tolerable for the system, then Q will be
1000 0
Q=
0 100
u = −Kf x ⇒ u = −12.2829x2 , with this control law the eigenvalues of closed loop
Chapter 5 159
Daniel Abebe Modern Control Systems Lecture Note
system ẋ = (A − BKf ) x is
∥λI − (A − BKf )∥ = 0
λ 0 −1 0 0 0
− + =0
0 λ 0 −2 0 12.2829
λ+1 0
=0
0 λ + 14.2829
(λ + 1)(λ + 14.2829) = 0 ⇒ λ1 = −1 λ2 = −14.2829
If
2 0
Q=
0 2
and R = 20 the value of P and Kf is
500 0 h i
P= K= 0 1
0 20
u = −Kf x ⇒ u = −x2 , with this control law the eigenvalues of closed loop system
ẋ = (A − BKf ) x is
∥λI − (A − BKf )∥ = 0
λ 0 −1 0 0 0
− + =0
0 λ 0 −2 0 1
λ+1 0
=0
0 λ+3
(λ + 1)(λ + 3) = 0 ⇒ λ1 = −1 λ2 = −3
6. For the following system, determine an optimal full state feedback gain matrix Kf .
1 1 0
ẋ = x + u
0 2 1
h i
y= 1 0 x
Solution
From Riccati equation
Q + PA + A⊤ P − PBR−1 B⊤ P = 0
Let
p11 p12 1 0
P= Q = C⊤
1 C1 =
R = 1;
p21 p22 0 0
Chapter 5 160
Daniel Abebe Modern Control Systems Lecture Note
1 0 p p 1 1 1 0 p p
+ 11 12 + 11 12
0 0 p21 p22 0 2 1 2 p21 p22
p11 p12 0 h i p p
− 0 1 11 12 = 0
p21 p22 1 p21 p22
1 0 p p + 2p12 p11 p12 p p p p
+ 11 11 + − 12 21 12 22 = 0
0 0 p21 p21 + 2p22 p11 + 2p21 p21 + 2p22 p22 p21 p22 p22
1 + 2p11 − p12 p21 p11 + 3p12 − p12 p22 0 0
=
2
p11 + 3p12 − p12 p22 2(p21 + 2p22 ) − p22 0 0
Solving for P results in two values
−0.3702 −0.5095 19.4317 6.3138
P1 = P2 =
−0.5095 3.7265 6.3138 6.0777
1.4259 −1.9626 −0.4875 0.1584
P3 = P4 =
−1.9626 2.2735 0.1584 −0.0777
Since P has to be positive definite (all eigenvalues of P must be greater than or
equal to zero). Checking for the positive definiteness
−0.4326 3.5653
eig(P1 ) = eig(P2 ) =
3.7890 21.9442
−0.1582 −0.5415
eig(P3 ) = eig(P4 ) =
3.8575 −0.0236
19.4317 6.3138
Hence P = P2 = with this value of P, the value of feedback
6.3138 6.0777
gain matrix Kf is
h i
Kf = R−1 B⊤ P = 6.3138 6.0777
The control law becomes u = −Kf x = −6.3138x1 − 6.0777x2 , with this control law
the eigenvalues of closed loop system ẋ = (A − BKf ) x is
∥λI − (A − BKf )∥ = 0
λ 0 1 1 0 0
− + =0
0 λ 0 2 6.3138 6.0777
λ−1 −1
=0
6.3138 λ + 4.0777
(λ − 1)(λ + 4.0777) + 6.3138 = 0 ⇒ λ1 = −2.9990 λ2 = −0.0787
Chapter 5 161
Daniel Abebe Modern Control Systems Lecture Note
7. For the following system, determine an optimal full state feedback gain matrix K.
1 0 0 1
ẋ = 1 −1 1 x + 1 u
0 0 −2 0
h i
y= 1 0 0 x
Solution
As shown in chapter two example 3, the system is not completely controllable. So,
we have to decompose the system into controllable and uncontrollable modes as
shown in chapter two example 3. So we can design a state feedback controller for
controllable mode (A11 , B1 ). From Riccati equation
Q + PA + A⊤ P − PBR−1 B⊤ P = 0
Let
p11 p12 1 1
P= Q = C⊤
1 C1 =
R = 2;
p21 p22 1 1
1 1 p p 2 −4 2 1 p p
+ 11 12 + 11 12
1 1 p21 p22 1 0 −4 0 p21 p22
p11 p12 1 h i p p
− 2 1 0 11 12 = 0
p21 p22 0 p21 p22
p211
(p11 p12 )
− + 4p11 + 2p12 + 1 2p12 − 4p11 + p22 − + 1
2 2 =0
(p11 p12 ) p212
2p12 − 4p11 + p22 − +1 − − 8p12 + 1
2 2
Chapter 5 162
Daniel Abebe Modern Control Systems Lecture Note
Since P has to be a positive definite matrix, the valid answer is (check using all
eigenvalues of P (eig(P)))
8.3007 0.1240
P2 =
0.1240 32.4696
Chapter 5 163
Daniel Abebe Modern Control Systems Lecture Note
u = −Kf x ⇒ u = −7.2915x1 − 2x2 , with this control law the eigenvalues of closed
loop system ẋ = (A − BKf ) x is
∥λI − (A − BKf )∥ = 0
λ 0 2 −4 7.2915 2
− + =0
0 λ 1 0 0 0
λ + 5.2915 6
=0
−1 λ
λ2 + 5.2915λ + 6 = 0 ⇒ λ1 = −3.6457 λ2 = −1.6458
5.7 Exercises
1. A system described by
ẋ = −2x + u
ẋ = −2x + 2u
ẋ1 = x2
ẋ2 = −x2 + u
Chapter 5 164
Daniel Abebe Modern Control Systems Lecture Note
ẋ1 = −x1 + u
ẋ2 = x1
5. Use Riccati equation to determine the feedback control law for a system
ẋ = −x + u
which minimize
Z 1
1
J= (3x2 + u2 ) dt
2 0
6. Find an optimal feedback gain Kf for the following systems with the given initial
conditions
(a)
0 1 0 1
ẋ = x + u, x0 =
−4 0 1 1
that minimize
Z ∞
1
J= (x⊤ Qx + u⊤ Ru) dt
2 0
(b)
0 1 0 1
ẋ = x + u, x0 =
2 −1 1 1
that minimize
Z ∞
1
J= (x⊤ Qx + Ru2 ) dt
2 0
(c)
−1 0 1 0
ẋ = x + u, x0 =
0 −4 1 1
that minimize
Z ∞
1
J= (x⊤ Qx + u⊤ Ru) dt
2 0
Chapter 5 165
Daniel Abebe Modern Control Systems Lecture Note
(d)
0 1 0 1
ẋ = x + u, x0 =
−6 −8 1 1
that minimize
Z ∞
1
J= (x⊤ Qx + u⊤ Ru) dt
2 0
(e)
0 1 0 1
ẋ = x + u, x0 =
−6 0 1 1
that minimize
Z ∞
1
J= (x⊤ Qx + u⊤ Ru) dt
2 0
7. Design an optimal state feedback matrix Kf that minimize the cost function
1 ∞ ⊤
Z
J= (x Qx + u⊤ Ru) dt
2 0
for
using MATLAB.
Chapter 5 166
Chapter 6
167
Daniel Abebe Modern Control Systems Lecture Note
ẋ3 = x4 = 0
Hence, the equilibrium point for the IPC dynamic system is (x1 , 0, nπ, 0). We can’t draw
any conclusion about x1 since it vanishes during the equilibrium point calculation.
Since the tanks are coupled, the equilibrium point of each tank cannot be set individually,
for instance, to set the tank 4 to a height h4e , a voltage V1e is required which in turn
affect the level of tank 1 and tank 2.
p 1 γ1 k1 V1e + (1 − γ2 )k2 V2e
γ1 k1 V1e + (1 − γ2 )k2 V2e = a1 2gh1e =⇒ h1e =
2g a1
p 1 γ2 k2 V2e + (1 − γ1 )k1 V1e
γ2 k2 V2e + (1 − γ1 )k1 V1e = a2 2gh2e =⇒ h2e =
2g a2
ẋ = f (x) (6.1.6)
x = xe + ∆x (6.1.7)
Chapter 6 168
Daniel Abebe Modern Control Systems Lecture Note
The Taylor series expansion of the function about this point xe is given as
′ 1 ′′
ẋe + ∆ẋ = f (x + ∆x) = f (xe ) + f (xe )∆x + f xe (∆x)2 + HOT (6.1.8)
2!
Neglecting higher-order terms, the function can be approximated by a linear function as
′
ẋe + ∆ẋ = f (xe ) + f (xe )∆x (6.1.9)
Similarly for multi-variable function f (x, u) the approximation of the function about
operating point xe and ue given as
∂f ∂f
f (x, u) = f (xe , ue ) + (x − xe ) + (u − ue ) (6.1.12)
∂x ∂u
If the operating point is the equilibrium point ẋ = f (xe , ue ) = 0, which means
∂f ∂f
ẋ = f (x, u) = x+ u (6.1.13)
∂x ∂u
∂f ∂f
Which implies that A = and B = evaluated at the equilibrium point (xe , ue ).
∂x ∂u
We have a nonlinear system with n number of states dynamics ẋi = fi (xi , uj ) i =
1, 2, . . . n, j = 1, 2, 3, . . . , m
ẋ f (x , x , . . . , xn , u1 , u2 , . . . , um )
1 1 1 2
ẋ2 f2 (x1 , x2 , . . . , xn , u1 , u2 , . . . , um )
= (6.1.14)
.. ..
. .
ẋn f3 (x1 , x2 , . . . , xn , u1 , u2 , . . . , um )
Chapter 6 169
Daniel Abebe Modern Control Systems Lecture Note
While
∂f1 ∂f1 ∂f1
∂u1 ...
∂u2 ∂
∂f ∂f2 ∂f2
2
B= ... (6.1.18)
∂u1 ∂u2 ∂um
. .. ..
..
. .
∂fn ∂fn ∂fn
...
∂u1 ∂ ∂ |(xeq ,ueq )
In similar way the nonlinear system has p outputs yi = hi (x, u),i = 1, 2, . . . p the approx-
imation about equilibrium point will be
∂h1 ∂h1 ∂h1 ∂h1 ∂h1 ∂h1
y x1 x2 . . . xn u1 u2 . . . um
1 ∂x1 ∂x2 ∂xn ∂h1 ∂u2 ∂um
∂h ∂h2 ∂h2 ∂h2 ∂h2 ∂h2
2
y2 = x1 x2 . . . xn u1 u2 . . . um (6.1.19)
∂x1 ∂x2 ∂xn ∂u1 ∂u2 ∂um
.. . .. .. .. .. ..
. ..
. . . . .
∂hp ∂h ∂h ∂h ∂h ∂h
yp p p p p p
x1 x2 . . . xn u1 u2 . . . um
∂x1 ∂x2 ∂xn ∂u1 ∂u2 ∂um |(xeq ,ueq )
Chapter 6 170
Daniel Abebe Modern Control Systems Lecture Note
Chapter 6 171
Daniel Abebe Modern Control Systems Lecture Note
Chapter 6 172
Daniel Abebe Modern Control Systems Lecture Note
Chapter 6 173
Daniel Abebe Modern Control Systems Lecture Note
Thus the linearized model of IPC about equilibrium point (xe , ue ) = (x1 , 0, 0, 0, 0) is given
as
0 1 0 0 0
ẋ x1
1 2 2 2
0 − b(mp l + J) (mp l) g mp l + J
ẋ2 − 0 x
= q q 2
+ q
u
ẋ3 0 0 0 1 x3 0
ẋ4 bmp l mp lg(mp + mc ) mp l
0 x4
0 −
q q q
(6.1.40)
x
1
y1 1 0 0 0 x2 0
y= = + u
(6.1.41)
y2 0 0 1 0 x3 0
x4
The matrix A ∈ R4×4 is obtained by taking partial derivatives of each equation concerning
the states:
∂fi
Aij =
∂xj xe ,ue
p p
∂ ḣ1 a1 g a1 g/2h1e ∂ ḣ1 a3 g a3 g/2h3e
=− ·p =− , = ·p =
∂h1 A1 2gh1e A1 ∂h3 A1 2gh3e A1
p p
∂ ḣ2 a2 g a2 g/2h2e ∂ ḣ2 a4 g a4 g/2h4e
=− ·p =− , = ·p =
∂h2 A2 2gh2e A2 ∂h4 A2 2gh4e A2
p p
∂ ḣ3 a3 g a3 g/2h3e ∂ ḣ4 a4 g a4 g/2h4e
=− ·p =− , =− ·p =−
∂h3 A3 2gh3e A3 ∂h4 A4 2gh4e A4
Chapter 6 174
Daniel Abebe Modern Control Systems Lecture Note
Let
p p p
a1
g/2h1e a2
g/2h2e a3 g/2h3e
α1 = , α2 = , α3 =
A1 A2 A
p p p 1
a4 g/2h4e a3 g/2h3e a4 g/2h4e
α4 = , α5 = , α6 =
A2 A3 A4
With this he system matrix A is then,
−α1 0 α3 0
0 −α2
0 α4
A= (6.1.42)
0 −α5
0 0
0 0 0 −α6
The matrix B ∈ R4×2 is obtained by taking partial derivatives concerning the inputs:
∂fi
Bij =
∂uj xe ,ue
γ1 k1
0
A1
γ2 k2
0
B=
A2
(6.1.43)
(1 − γ2 )k2
0
A3
(1 − γ1 )k1
0
A4
with system outputs y1 = h1 (x, u) = h1 and y2 = h2 (x, u) = h2 , the output matrix
becomes
∂hi 1 0 0 0
Cij = = (6.1.44)
∂xj xe ,ue
0 1 0 0
Chapter 6 175
Daniel Abebe Modern Control Systems Lecture Note
Yi (s)
= Gij (s) (6.2.1)
Ui (s)
To determine the largest multiplicity pole, determine the largest multiplicity pole among
1×1 minors of the transfer function, then the largest multiplicity pole among 2×2 minors
of the transfer function so on. Stop at minors of size k × k if k is the first size for which
the previous multiplicity of the pole is greater than the current minor’s pole multiplicity.
The largest of all the multiplicities of the minors is the multiplicity of the pole.
Determining zeros of the MIMO system is a tricky problem. To start with, we know
that zeros of the system prohibit the control signal from affecting the dynamics of the
system. Hence, from this knowledge, a MIMO system has zero at location z0 if the z0
causes a reduction in rank of the MIMO system transfer function (i.e., if the z0 causes
zero rows in the transfer function of the MIMO system). Otherwise, the zero of one
transfer function will not be considered as the zero of the MIMO system.
Consider IPC dynamics derived in section 1.6.3, the transfer function of IPC will be
Y2 (s) mp l
= (6.2.4)
U (s) (mp + mc )mp lg − qs
Chapter 6 176
Daniel Abebe Modern Control Systems Lecture Note
As we can see, the system has only one input; hence, it is a single-input multiple-output
system (SIMO). The poles of the system are the poles of equation 6.2.3 and equation
6.2.4. The zero of the system is zero of equation 6.2.3, which will result in one zero row
in the system transfer function on equation 6.2.5.
(mp l)3 g mp l 2 + I
G11 (s) = − +
qs((mp + mc )mp lg − qs) sq
3 2
−(mp l) g + (mp l + I)((mp + mc )mp lg − qs)
=
sq((mp + mc )mp lg − qs)
For the linearized model of the quadruple tank level control system in section 6.1.2, if
height h1 and h2 are required to be controlled, we have two outputs and two [Link]
transfer function matrix G(s) relates inputs to outputs in the Laplace domain:
Chapter 6 177
Daniel Abebe Modern Control Systems Lecture Note
γ1 k1 α3 (1 − γ2 )k2
A1 (s + α1 ) A3 (s + α1 )(s + α5 )
G(s) = C(sI − A)−1 B = (6.2.10)
α4 (1 − γ1 )k1 γ2 k2
A4 (s + α2 )(s + α6 ) A2 (s + α2 )
The transmission zeros are the values of s that make the transfer function matrix lose
rank. We find them by solving:
γ1 γ2 k1 k2 α3 α4 (1 − γ1 )(1 − γ2 )k1 k2
det G(s) = − (6.2.11)
A1 A2 (s + α1 )(s + α2 ) A3 A4 (s + α1 )(s + α2 )(s + α5 )(s + α6 )
Chapter 6 178
Daniel Abebe Modern Control Systems Lecture Note
solving for the transmission zeros and poles of the system by setting det G(s),
1.8200 × 10−06
0.0069
||G(s)|| = −
s + 0.028 (s + 0.028)(s + 0.020)
2
0.006889(s + 0.04s + 0.0001358)
=
(s + 0.028)2 (s + 0.020)2
Chapter 6 179
Daniel Abebe Modern Control Systems Lecture Note
6.3 Examples
1. Find the poles and zeros of the following system
ẋ1 = x2 + u1
ẋ2 = x3 + u2
ẋ3 = 6x1 − 11x2 + 6x3 + u1 + u2
y1 = x1
y2 = x3
Solution
The state space representation of the system is
ẋ
1
0 1 0 x1 1 0
ẋ2 u
1
= 0 0 1 x2 + 0 1
ẋ3 u2
6 −11 6 x3 1 1
ẋ4
x1
y1 1 0 0 0 0 u
= x2 +
1
y2 0 0 1 0 0 u2
x3
The transfer function of the system is
Cadj(sI − A)B
G(s) =
det(sI − A)
s2 − 6s + 11 s + 6 1 1 0
1 0 0 2
s − 6s s 0 1
6
0 0 1
2
6s s s 1 1
=
s3 − 6s2+ 11s − 6
s2 − 6s + 12 s + 7
s2 + 6s s + s2
=
s3 − 6s2 + 11s − 6
The only number that can cause a zero row in the transfer function matrix is
s = 0. Hence s = 0 is the zero of the system. The poles of the system are zeros of
s3 − 6s2 + 11s − 6 = (s + 1)(s + 2)(s + 3) = 0, hence s = 1, s = 2 and s = 3 are the
poles of the system.
Chapter 6 180
Daniel Abebe Modern Control Systems Lecture Note
Solution
No number causes zero rows or zero columns in the transfer function matrix, hence
the system has no zeros. The poles of the system are
s s2 + 3 1 s+3
• poles of first minors, , , 2 and
s + 1 (s + 1)(s + 4) s + 1 s(s + 5)
• poles of the second minors
s2 + 3
s s+3 1
−
s+1 s(s + 5) s2 + 1 (s + 1)(s + 4)
Thus, the poles of the system are s = 0, s = −1, s = −4, s = −5, and s = ±i.
3. Find the equilibrium point/s and linearized model of the system at the equilibrium
point/s.
Solution
At equilibrium point/s the the state dynamics is zero ẋ1 = 0, ẋ2 = 0, ueq = 0, thus
syms x1 x2 u
eqn1 = 7 ∗ x1 + 2 ∗ sin(x2 ) − x42 + u3 ;
eqn2 = exp(x1 ) − 3 ∗ x2 − 1 + 5 ∗ x21 + sin(u)
[x1 , x2 , u] = solve(eqn1 == 0, eqn1 == 0, [x1 , x2 , u])
the only equilibrium point is (x1eq , x2eq , ueq = 0) = (0, 0, 0). The linearized model
Chapter 6 181
Daniel Abebe Modern Control Systems Lecture Note
about (0, 0, 0) is
∂ ∂
∂
f (x , x , u) f (x , x , u)
ẋ
1 ∂x1
1 1 2
∂x2
1 1 2
x1 ∂u f1 (x1 , x2 , u)
= + u
∂ ∂ ∂
ẋ2 f2 (x1 , x2 , u) f2 (x1 , x2 , u) x2 f2 (x1 , x2 , u)
∂x1 ∂x2 |xeq ∂u |ueq=0
7 2 cos x2 − 4x32 x1 3u2
= + u
x1
e + 10x1 −3 x2 cos u
|(x1 =0,x2 =0) |ueq=0
7 1 x 0
= 1 + u
1 −3 x2 1
x
∂ ∂ 1
y= h(x1 , x2 , u) h(x1 , x2 , u)
∂x1 ∂x2 |(x1 =0,x2 =0) x2
h i x
= ex1 ex2 1
|(x1 =0,x2 =0) x
2
h i x1
= 1 1
x2
4. Find the equilibrium point/s of the following system and the linearized model of
the system at the point/s.
Solution
At equilibrium point/s the the state dynamics is zero ẋ1 = 0, ẋ2 = 0, thus
x1 = 0, x2 = 0, x21 − x22 = 0
Which means all of the points where x21 x22 = 0 ⇒ x1 = x2 , or x1 = −x2 are
equilibrium points of the system. This set of points is called a limit cycle of the
system.
Chapter 6 182
Daniel Abebe Modern Control Systems Lecture Note
or at (x1 = −x2 ) is
∂ ∂
ẋ f1 (x1 , x2 ) f1 (x1 , x2 )
1 ∂x1 ∂x2
x
1
=
∂ ∂ x2
ẋ2 f2 (x1 , x2 ) f2 (x1 , x2 )
∂x1 ∂x2 |x
eq
3 + 2x1 2x2 x
= 1
4x1 x2 4x2 x1 x2
|(x1 =−x2 )
3 − 2x2 −2x2 x
= 1
−4x22 −4x22 x2
Chapter 6 183
Daniel Abebe Modern Control Systems Lecture Note
6.4 Exercises
1. For the following nonlinear system, find the equilibrium point and the linear model
at the equilibrium point.
(a) ẋ1 = −x1 + x2
(f) ẋ1 = −3x2
ẋ2 = −x1 − x2 − x21
ẋ2 = x1 − α(x32 − x2 )
(b) ẋ1 = −3x2
(g) ẋ1 = −x1 + x2 − x3 + αx31
ẋ2 = x1 − x2 − 2x32
ẋ2 = −x1 − x2
(c) ẋ1 = −x31 − 3x1 + x2 ẋ3 = x1 − x2 + βx3 + u
ẋ2 = −2x1
(h) ẋ1 = −x1 + x2
(d) ẋ1 = x21 − x22 ẋ2 = (x1 + x2 ) sin x1 − 3x2
ẋ2 = −2x1 x2
(i) z̈ + ż + z 3 + u = 0
(e) ẋ1 = −x1 + x2
(j) ẍ − µ(1 − x2 )ẋ + x = 0
ẋ2 = −x1 − x2 − x21
2. For DC motor dynamics on 1.7.1, find the transfer function from input voltage to
output speed and input torque to output speed.
3. For the following systems, find a transfer function from each input to outputs, zeros,
and poles of the system.
(a)
0 1 0 0 0
1 0 0
ẋ = 0 x + 0 1 u, y= x
0 1
0 1 0
−4 −4 −1 1 1
(b)
0 1 0 0 0
1 0 −2
ẋ = 0 1 x + 0 1 u, y= x
0
0 2 1
−4 −4 −1 1 0
(c)
0 1 0 0
−3 1 0
ẋ = 0 1 x + 1 u, y= x
0
4 0 1
−6 −11 −6 −3
(d)
1 1 0 0 1 1 0 −0.5
ẋ = 0 1 0 x + 1 0 u, y = −1 1 0 x
0 1 1 0 1 0 1 1
Chapter 6 184
Chapter 7
Analog Simulation
Analog computers are particularly well-suited for solving differential equations, which
are equations that describe how a quantity changes over time. This is because analog
computers can continuously measure and adjust the voltage in the circuit to represent
the changing quantity.
Analog computers were widely used in the 20th century to solve problems in engineering,
physics, and mathematics. However, they have been largely replaced by digital comput-
ers, which are more accurate and versatile. However, analog computers are still used in
some specialized applications, such as aircraft flight control systems.
• The mathematical model of the problem is converted into an analog block diagram.
The block diagram shows how the different parts of the problem are related to each
other.
• The analog block diagram is used to create an electrical circuit. The electrical
components are connected together in a way that represents the mathematical re-
lationships in the block diagram.
• The electrical circuit is powered up, and the voltages at different points in the
circuit are measured. The voltages represent the solutions to the problem.
• The voltages are converted back to the original problem variables. This gives the
final solution to the problem.
Analog computers have a number of advantages over digital computers for certain types
of problems. For example, analog computers can solve differential equations in real time.
This makes them ideal for applications such as aircraft flight control systems, where the
computer needs to be able to respond to changes in the environment very quickly.
185
Daniel Abebe Modern Control Systems Lecture Note
However, analog computers also have a number of disadvantages. For example, analog
computers are less accurate than digital computers. They are also more difficult to
program and maintain.
x2
+ x x
+
x1 - y x1 x2 x3
x3
Summation Multiplication
x0
x x x
1
t
x x0 x( )d
Integration 0 Sign Inversion
Consider an unforced first-order system (zero input signal), with initial condition x0
ẋ = x x0 = 1 (7.2.1)
Integrating both sides
Z t
x = x0 + x(τ ) dτ (7.2.2)
0
Chapter 7 186
Daniel Abebe Modern Control Systems Lecture Note
x0
t
x x0 x ( )d
0
As we can understand from Figure 7.2, it is not necessary to know the input of the
integrator to solve equation 7.2.1. The only requirement is that the input to the integrator
must always be equal to the output. Hence, the concept of feeding the unknown output
back to the integrator is the solution to the differential equation in an analog computer.
ẍ + ẋ + x = 0 x0 = 1 ẋ0 = 0 (7.2.3)
ẋ1 x2
ẋ2 = −x1 − x2 (7.2.4)
x2 (0) x1 (0)
x2 x1
+
Chapter 7 187
Daniel Abebe Modern Control Systems Lecture Note
R1 Rf
v1
R2
v2 −
.. vo
. +
Rn
vn
2. Integration
C
R1
v1
R2
v2 −
.. vo
. +
Rn
vn
Chapter 7 188
Daniel Abebe Modern Control Systems Lecture Note
C1
Rf
v1
v2 −
vo
C2 +
..
.
vn
Cn
3. Differentiation
For the differentiator circuit, the voltage transfer function is,
dv1 dv2 dvn
vo = −Rf C1 − Rf C2 − . . . − Rf Cn
Zdt t dt dt
1 1 1
vo = vo (0) − v1 (τ ) + v2 (τ ) + . . . + vn (τ ) dτ (7.3.5)
0 CR1 CR2 CRn
4. Sign inversion
R1 Rf 1 R2 Rf 2
v1
− −
vo
+ +
For the sign inversion circuit, the voltage transfer function is,
Rf 2 Rf 1
vo = v1 (7.3.7)
R1 R2
Rf 2 Rf 1
If = 1, we have vo = v1 .
R1 R2
Chapter 7 189
Daniel Abebe Modern Control Systems Lecture Note
The circuit for an nth order integrator consists of n operational amplifiers (op-amps)
connected in series with capacitors and resistors. The output of each op-amp is connected
to the input of the next op-amp, and the input signal is applied to the first op-amp. The
output of the last op-amp is the integrated signal.
The circuit for an nth order differentiator is similar, but with capacitors and resistors
connected in parallel instead of in series. The input signal is applied to the first capacitor,
and the output of each op-amp is connected to the input of the next op-amp. The output
of the last op-amp is the differentiated signal.
By combining these circuits, we can solve differential equations by converting them into
an analog block diagram, creating an electrical circuit based on the block diagram, and
measuring the voltages at different points in the circuit to obtain the solution to the
problem.
dn dn−1 d
n
x + n−1
x + ... + x + x = 0 (7.4.1)
dt dt dt
with n initial conditions.
dn dn
If x is known, then all the rest variables can be determined by integrating x, or
dtn dtn
dn−1 d dn
if n−1 x, . . ., x and x are known, we can determine n x. This implies that we are
dt dt dt
only expected to know the input and output of the integrator. The analog circuit that
dn
integrates n x is shown in Figure 7.8.
dt
dn x dn−1 x dn−2 x
− x
dtn dtn−1 dtn−2 . . .
Chapter 7 190
Daniel Abebe Modern Control Systems Lecture Note
dn−1 d dn
The analog circuit that can sum n−1 x, . . ., x and x to get n x the is shown in
dt dt dt
Figure 7.8.
R1 Rf
dn−1 x
dtn−1 R2
dn−2 x
−
dtn−2 dn x
.. dtn
. +
Rn
x
Combining the two analog circuits will give us an analog circuit used to solve nth order
ordinary differential equations shown in Figure 7.10.
d nx d n1 x dx
dt n dt n1 dt x
Chapter 7 191
Daniel Abebe Modern Control Systems Lecture Note
For instance, if we were to solve a problem involving rocket speed on an analog computer,
it would be impossible to relate 1 volt on the computer to 1cm/sec in the real problem.
However, we could relate 1 volt on the computer to 1km/sec in the problem. This illus-
trates the need to scale the problem so that it will fit the computer.
Similarly, in the same problem, an acceleration term might have a maximum value of
1cm/sec2 . Using the basic unit of miles for length as determined for scaling the velocity
would give a voltage of approximately 0.002 volts, which corresponds to the acceleration
term. However, this signal is in the noise level of the computer and would yield erroneous
results. Thus, it is necessary to scale up some parts of a problem and scale down other
parts.
Time scaling is also necessary in certain situations. For example, when studying a phys-
ical phenomenon that evolves much faster or slower than is convenient for the computer
operator, time scaling can be used. If we wish to study planetary motion, waiting for
a complete period of revolution about the sun would be impractical. Instead, an actual
time of one year may be represented by 10 seconds of computer time.
Equipment constraints and limitations require careful consideration when solving prob-
lems on analog computers. Scaling techniques must be employed to ensure that the
problem fits within the computer’s limitations. Time scaling can also be used to study
physical phenomena that evolve at different rates than is convenient for the computer
operator.
Chapter 7 192
Daniel Abebe Modern Control Systems Lecture Note
Scaling is performed based on the maximum voltage supplied to the analog computer and
the maximum magnitude that appears in the solution of the differential equation. The
input and output of the analog computer are scaled down using the voltage divider rule,
while operational amplifiers are used to scale them up. Similarly, the initial conditions
can be scaled up or down.
vin
vin vo
vo vin
We must scale the problem so that we can represent 500 centimeters, ymax , as an analog
voltage. We also know
V x = Lx x Vy = Ly y
Chapter 7 193
Daniel Abebe Modern Control Systems Lecture Note
Vx Lx x Vy Ly y
kGLx
Which means a = . So the gain of the potentiometer must have a setting
Ly
aLy
k=
Lx
In this example
10 10 aLy 100 ∗ 0.1
Lx = =2 Ly = = 0.1 ⇒ k = = =5
5 100 Lx 2
Generally
Level Out COEF
k= (7.5.1)
Level In gain
Where
• Level Out is the level of the variable immediately following the pot to be set,
Chapter 7 194
Daniel Abebe Modern Control Systems Lecture Note
τc = βt (7.5.2)
1
Time scaling can be achieved by multiplying the gains of every integrator by .
β
Therefore, we can adjust the previous pot setting formula to take care of magnitude and
time scaling. We simply use a single pot for each input to the integrators such that:
Level Out COEF 1
k= (7.5.3)
Level In gain β
Since there is no time variation for initial conditions and non-integrating operational
amplifiers, time scaling will not be applied to them.
1. Convert the transfer function into differential equations using direct or indirect
programming or standard programming methods.
4. Draw integrators needed and put a pot between each pair of integrators.
5. Determine the scaling factors based on the maximum variable occurring in the
differential equation, as described in step 3.
To develop a scaled analog simulation diagram of state space representation, use the
following procedures.
4. Draw all integrators needed and put a pot between each pair of integrators.
5. Determine the scaling factors based on the maximum variable occurring in the
differential equation, as described in step 3.
Chapter 7 195
Daniel Abebe Modern Control Systems Lecture Note
and also
We note that the initial states of the integrators can be rewritten in terms of the initial
states of y and ẏ, and upon the input and its derivative. Note that when we take the
Laplace transform of a signal (like u(t)), we assume the signal to be zero and all its
derivatives to be zero. With is assumptions equations 7.6.2 and 7.6.3 simplifies to
Let a11 = a0 −a2 b0 , a12 = a1 −a2 b1 , a21 = −b0 (a1 −a2 b1 ) and a22 = a0 −a2 b0 −b1 (a1 −a2 b1 ),
thus we obtain
α22 y(0− ) − α12 ẏ(0− )
x(0− ) = (7.6.9)
α11 α22 − α12 α21
−α12 y(0− ) − α22 ẏ(0− )
ẋ(0− ) = (7.6.10)
α11 α22 − α12 α21
Chapter 7 196
Daniel Abebe Modern Control Systems Lecture Note
7.7 Examples
1. Determine the circuit to produce the output voltage given by
Solution:
The analog circuit diagram that must be patched on the analog computer
50k 50k
v1
5k
v2 vo
50k 50k
v3
5k
v4 - 50k
Solution:
The analog circuit diagram must be patched on the analog computer
v1 50k
.350
5k
50k
v2 5k
5k
.524 - vo
-
10v 5k
.260
Chapter 7 197
Daniel Abebe Modern Control Systems Lecture Note
.350
v1 1
.524 vo
v2 10
-10v 1
3. For the following second-order linear differential equation and initial conditions,
determine the analog circuit diagram.
ẍ + 0.5ẋ + x = 4 ẋ0 = 1 and x0 = 0 (7.7.1)
Solution:
The analog circuit diagram that must be patched on the analog computer
0.4 x0 x0
10v
x x x
0.5
1
Chapter 7 198
Daniel Abebe Modern Control Systems Lecture Note
Y (s) X (s)
X (s) U (s)
This implies that the analog circuit diagram of the two differential equations can
be designed same way and combined as follows
The analog circuit diagram of the transfer function from the state to the output is
a2
x – 1 y
a1 1
–
a0
Figure 7.18: Analog Circuit Diagram of Transfer Function From State to Input
The analog circuit diagram of the transfer function from the input to the state is
u x x x
– –
b1
b2
1
Figure 7.19: Analog Circuit Diagram of Transfer Function from Input to State
Combining the two analog circuit diagrams, the transfer function from input to
output is,
Chapter 7 199
Daniel Abebe Modern Control Systems Lecture Note
a2
u y
x x x
– – 1
b1 a0
1
b2
a1
Figure 7.20: Analog Circuit Diagram of Transfer Function from Input to Output
5. Simulate the following system and magnitude scale so that no overload occurs and
the maximum amplifier range is used
Solution:
10V 10V
5 4
10V G1
6
x x
G2
– 1 1 – 2
G3
0.8
2
0.25
10
3 5 -1
(b) Step 2: This step is not necessary since the maximum values are already
specified in the problem statement.
(c) Step 3: Prepare a scaling table
(d) Step 4: Write the scaled equations for each amplifier.
1/5 1
k1 = , if we select G1 = 10, then k1 = 0.5. Similarly, computing for
1/25 G1
Chapter 7 200
Daniel Abebe Modern Control Systems Lecture Note
x x
5
25
10V 10V
5 0.6 4 0.8
0.2
10V 10
6
x 10 x
10
– 1 1 – 2
10
0.5
0.8
2
0.25
1
3 5 -1
0.5
Chapter 7 201
Daniel Abebe Modern Control Systems Lecture Note
7.8 Exercises
1. For the following system, find an analog circuit diagram
Chapter 202
Appendix A
where:
• Is a scalar quantity, making it easier to work with than vector quantities like force.
i
A.1.3 Connection to Newton’s Second Law
Consider a particle with mass m moving in a generalized coordinate system. The kinetic
1
energy of the particle is T = mq̇2 , and its potential energy is V = V (q). The Lagrangian
2
of the system is:
1
L = mq̇2 − V (q) (A.1.3)
2
Applying the Euler-Lagrange equation:
∂V d
− − (mq̇) = 0 =⇒ F = mq̈ (A.1.4)
∂q dt
This recovers Newton’s second law.
ii
A.2.3 Differential of a Function
For a function f (t), the Taylor series expansion about a point t∗ is given as:
1 ∂ 2f
∗ ∂f
∆f = f (t ) + ∆t + ∆t2 + · · · − f (t∗ ) (A.2.3)
∂t t∗ 2! ∂t2 t∗
• Necessary Condition: The first derivative must vanish at the point (i.e., df = 0
at t∗ ).
• Sufficient Condition:
1 ∂ 2J 2
∗ ∂J
∆J = J(x (t)) + δx + 2
δx + . . . − J(x∗ (t))
∂x x∗ (t) 2! ∂x x∗ (t)
∂J 1 ∂ 2J 2
Here, δJ = δx is called the first variation, and δ 2 J = δx is called the second
∂x ∗
2! ∂x2
variation of the functional about a function x (t).
∂J
∆J = δJ = δx (A.3.2)
∂x
iii
A.4 Optimization of an Unconstrained Functional
For a function to have an optimum value:
• Sufficient Condition:
Using integration by parts for the second term, we derive the Euler-Lagrange equation:
∂φ d ∂φ
− =0 (A.4.5)
∂x dt ∂ ẋ
iv
Boundary Conditions
∂φ
For the boundary condition to be zero, either or δx must be zero at the boundary.
∂ ẋ
This leads to four cases:
1. Case-I: x is fixed at t0 and tf (fixed end trajectory). Here, x(t0 ) = C1 and x(tf ) =
C2 , where C1 and C2 are constants. Since δx = 0 at both t0 and tf , there are no
∂φ
restrictions on .
∂ ẋ
∂φ
2. Case-II: x is fixed at t0 and free at tf . Here, x(t0 ) = C1 , and = 0 at tf because
∂ ẋ
∂φ
δx is unknown at tf . No restrictions apply to at t0 .
∂ ẋ
∂φ
3. Case-III: x is fixed at tf and free at t0 . Here, x(tf ) = C2 , and = 0 at t0
∂ ẋ
∂φ
because δx is unknown at t0 . No restrictions apply to at tf .
∂ ẋ
∂φ
4. Case-IV: x is free at t0 and tf . Here, = 0 at both t0 and tf because δx is
∂ ẋ
unknown at both times.
tf
∂ 2φ ∂ 2φ ∂ 2φ
Z
2 1 d 2 2
δ J= 2
− (δx) + (δ ẋ) dt (A.4.7)
2 t0 ∂x dt ∂ ẋ∂x ∂ ẋ2
∂ 2φ ∂ 2φ ∂ 2φ
d
− ≥ 0 and ≥0 (A.4.8)
∂x2 dt ∂x∂ ẋ ∂ ẋ2
∂ 2φ ∂ 2φ ∂ 2φ
d
2
− ≤ 0 and ≤0 (A.4.9)
∂x dt ∂x∂ ẋ ∂ ẋ2
v
A.5 Typical Functional in Control Systems
In control systems, a typical functional has the form J(x, ẋ, u, u̇, t):
Z tf
J(x, ẋ, u, u̇, t) = φ(x, ẋ, u, u̇, t) dt (A.5.1)
t0
Z tf
∆J = (φ(x + δx, ẋ + δ ẋ, u + δu, u̇ + δ u̇, t + δt) − φ(x, ẋ, u, u̇, t)) dt (A.5.2)
t0
Z tf
∆J = (φ(x + δx, ẋ + δ ẋ, u + δu, u̇ + δ u̇) − φ(x, ẋ, u, u̇)) dt (A.5.3)
t0
Using Taylor Series expansions, the first variation of the functional is:
Z tf
∂φ ∂φ ∂φ ∂φ
∆J = δx + δ ẋ + δu + δ u̇ + hot dt = 0 (A.5.4)
t0 ∂x ∂ ẋ ∂u ∂ u̇
The first variation of the functional is:
Z tf Z tf
∂φ ∂φ ∂φ ∂φ
δJ = δx + δ ẋ dt + δu + δ u̇ dt = 0 (A.5.5)
t0 ∂x ∂ ẋ t0 ∂u ∂ u̇
Z tf
Υ(x, ẋ, u)dt (A.6.1)
t0
ẋ = Ax + Bu (A.6.2)
vi
A.6.2 General Function Definition
To incorporate the constraints into the optimization problem, we define a general function
that includes both the cost function and the constraints:
Z tf
Υ(x, ẋ, u) + λ⊤ (Ax + Bu − ẋ) dt
J=
t
Z 0tf
= φ(x, ẋ, u)dt (A.6.4)
t0
∂φ d ∂φ
− =0 (A.6.5)
∂x dt ∂ ẋ
∂φ d ∂φ
− =0 (A.6.6)
∂u dt ∂ u̇
∂φ d ∂φ
− =0 (A.6.7)
∂λ⊤ dt ∂ λ̇⊤
where Fi are the generalized forces, including friction and external input.
vii
A.8 Example:-1 Inverted Pendulum on a Cart System Dynam-
ics
The objective of the inverted pendulum on a cart control system is to maintain the
pendulum in an upright position while simultaneously controlling the movement of the
cart along a track. This system is a classic example in control theory and is used to study
stability and control strategies.
lsinθ
y
θ
lcosθ
x
viii
Kinetic and Potential Energy
The total kinetic energy of the system is:
1 1 1
Ke = mc ẋ2c + mp (ẋ2p + ẏp2 ) + I q̇22 (A.8.4)
2 2 2
Substituting the velocities:
1 1
Ke = (mc + mp )q̇12 + mp lq̇1 q̇2 cos q2 + (mp l2 + I)q̇22 (A.8.5)
2 2
V = mp gl cos q2 (A.8.6)
Lagrangian Formulation
The Lagrangian is defined as:
L = Ke − V (A.8.7)
Substituting Ke and V :
1 1
L = (mc + mp )q̇12 + mp lq̇1 q̇2 cos q2 + (mp l2 + I)q̇22 − mp gl cos q2 (A.8.8)
2 2
where Fi are the generalized forces, including friction and external input.
Generalized Forces
• For q1 (cart motion):
F1 = u − bc q̇1 (A.8.10)
where u is the external force applied to the cart, and bc q̇1 is the friction force.
ix
Equations of Motion
The equations of motion for the system are derived as follows:
State-Space Representation
Let:
Linearization
Using the small angle approximation about θ = 0◦ , we assume cos θ ≈ 1, sin θ ≈ θ, and
θ̇ ≈ 0. Substituting these approximations, the linear model is:
0 1 0 0 0
ẋ x1
1 2 2 2
b c (m p l + I) (m p l) g b p m p l m p l + I
0
ẋ2 − −
x2
= λ λ λ + λ
u
ẋ3 0 0 0 1 x3 0
ẋ4
bc mp l mp lg(mp + mc ) bp (mp + mc ) x mp l
0 − 4 −
λ λ λ λ
(A.8.17)
x
1
y1 1 0 0 0 x2
y= =
(A.8.18)
y2 0 0 1 0 x3
x4
x
A.9 Example:-2 Ball and Beam System
The ball and beam system is a classic control theory problem that involves balancing a
ball on a beam, which can tilt to control the ball’s position.
I,m
x
xi
Total Kinetic and Potential Energy
• Total kinetic energy:
1 2 2 1 I2
K = Ktrans + Krot = (I1 + mx )θ̇ + m 1 + ẋ2 (A.9.5)
2 2 mr2
I2
• Let λ = 1 + be a constant.
mr2
• Potential energy:
Euler-Lagrange Equations
The Euler-Lagrange equation of the system is:
d ∂L ∂L
− = Fi (A.9.8)
dt ∂ q̇ ∂q
Where Fi are the generalized forces or torques, including friction and external input. In
this case, assuming the pivot friction is negligible, the only external force applied to the
system is the torque at the pivot joint τ .
State-Space Representation
With state variables x1 = x, x2 = ẋ, x3 = θ, and x4 = θ̇, the dynamics of the system are
described by:
ẋ1 = x2
mx24 x1 − g sin x3
ẋ2 =
λ
ẋ3 = x4
τ − 2mx1 x2 x3 − mgx1 cos x3
ẋ4 =
I + mx21
xii
Linearized Model State-Space Representation
Using small-angle approximations (sin x3 ≈ x3 and cos x3 ≈ 1) and assuming the ball is
at the center (x1 = 0), the linearized state-space equations become:
ẋ1 = x2
−gx3
ẋ2 =
λ
ẋ3 = x4
τ − mgx1
ẋ4 =
I
The state-space representation of the linearized system is:
0 1 0 0 0
−g
0 0 0 0
ẋ =
λ x + τ (A.9.11)
0
0 0 0 1
−mg 1
0 0 0
I I
The control objective is to regulate the position of the ball; thus, the output must be
x1 . The output equation is:
h i
y= 1 0 0 0 x (A.9.12)
xiii
Appendix B
Gdc ( s ) Gr ( s )
U(s) Y(s)
xiv
B.1.2 Test 2 (2024 G.C.)
Part I:-Conceptual Question(3 points)
1. Define stabilizability and detectability for a dynamic system and the requirements
that the system dynamics should have to be detectable and stabilizable. (3 points)
1. Show that the solution of the state equation will be changed with a similarity
transformation. (2.5 points)
2. Derive an expression for the solution of the state space equation.(2.5 points)
Consider the RLC Circuit, with system parameters L = 0.1H, C = 0.1F , and R = 1Ω.
R L
vs C vo
xv
1. If the designer aims to regulate the voltage across the capacitor.
2. If the output is the voltage across the inductor, determine the output state space
representation of the system. (2 points)
3. Compare the result on 2 with the result on 1a and state their difference. (0.5 points)
u = mo g
ms
x Scale
k d
If the output, which is shown on the scale, is the mass of the object to be measured
kx
(mo = ). (use g = 10m/s2 )
g
(a) Derive a state space representation of the system.(2.5 points)
(b) Find the eigenvalues and eigenvectors of the system.(1.5 points)
(c) Determine whether the system is stable or not.(0.5 points)
(d) Determine the state transition matrix of the system.(2.5 points)
xvi
(e) Determine zero input response of the system for initial condition x(0) =
h i⊤
1 0 .(2.5 points)
(f) Determine the zero state response of the system under a step input.(2.5 points)
2. If the maximum compression of the spring is 10cm, determine the maximum mass
the system can handle.(0.5 points)
3. If you are the designer and you want to increase the maximum mass the system
can handle 5 times,
1. Consider the RLC Circuit, with system parameters L = 0.1H, C = 0.1F and
R = 1Ω. If the output of the system is the voltage across the capacitor.
R L
vs C vo
xvii
(e) If C is very large, check for system observability and controllability and state
your conclusion? (2 points)
2. For the following completely controllable and completely observable system, (10
points)
0 1 0 0
ẋ = 0 1 x + 0 u
0
−5 −10 −5 1
h i
y= 0 0 1 x (B.2.1)
(a) Determine the feedback gain Kf so that the control law is u = −Kf x, assuming
the desired poles of the closed loop system are at s = −2±j4 and s = −2.(3.5
points)
(b) Design an full order observer gain Ke or Lc such that the eigenvalues of the
observer are located at s = −5 ± j4 and s = −10. (3.5 points)
(c) Draw the block diagram implementation of the full state feedback with full
state observer. (3 points)
(a) Determine the reduced observer gain matrix Ke or Lc assume that desired
observer poles to be at s = −5, s = −5. (4 points)
(b) Determine the reduced-order state observer dynamics and reduced-order state
observer output equations(open loop system dynamics).(3 points)
(c) Draw the overall block diagram for the above system (for question 3).(2
points)
xviii
(b) Design the linear quadratic Regulator (LQR) feedback controller Kf . (2
points)
(c) The control signal u. (1 points)
(d) The closed loop eigenvalues of the system. (2 points)
(e) Draw the block diagram for the LQR controller. (2 points)
xix
5. Design a reduced state observer based on the poles p1 , p2 identified in the previous
part for the system in equation B.2.6. (7 points)
(a) Determine the reduced state observer gain matrix Lc , that locate the pole of
the reduced order observer at p = p1 + p2 .
(b) Determine the output equation of the observer.
(c) Determine the reduced-order state observer dynamic equation.
(d) Determine the state feedback control law u with full state observer.
(e) Draw a block diagram representation of the full state feedback system con-
troller with reduced state observer.
ωn2
G2order (s) = (B.2.7)
s2 + 2ξωn s + ωn2
xx
B.2.3 Final Exam 3 (2020 G.C.)
Part I:-Conceptual Questions(10 points)
1. Consider the following simple mass spring damper system (11 points)
m u
b
x
2. For the following completely controllable and completely observable system, (10
points)
0 1 0 0
ẋ = 0 1 x + 0 u
0
−5 −10 −5 1
h i
y= 0 1 0 x (B.2.8)
(a) Determine the feedback gain Kf so that the control law is u = −Kf x, as-
suming the desired pole of the closed loop system is at s = −2 ± j4 and
s = −2.(3.5 points)
(b) Design an full order observer gain Ke or Lc such that the eigenvalues of the
observer are located at s = −5 ± j4 and s = −10. (3.5 points)
xxi
(c) Draw the block diagram implementation of the full state feedback with full
state observer. (3 points)
(a) Determine the reduced observer gain matrix Ke or Lc assume that desired
observer poles to be at s = −5, s = −5. (4 points)
(b) Determine the reduced-order state observer dynamics and reduced-order state
observer output equations(open loop system dynamics).(3 points)
(c) Draw the overall block diagram for the above system (for question 3).(2
points)
3. What is optimal control? List and give an example of the most popular optimal
control cost functions (performance indexes). (6 points)
xxii
Part II:- Calculation (34 points)
ẋ1 = x2
ẋ2 = x2 + u
y = x1 (B.2.13)
(a) Determine the feedback gain Kf so that the control law is u = −Kf x, assuming
the desired poles of the closed loop system are at s = −4 and s = −1.(3.5
points)
(b) Design an full order observer gain Lc such that the eigenvalues of the observer
are located at s = −4 and s = −1. (3.5 points)
(c) Draw the block diagram implementation of the system with full state feedback
using estimated states. (3 points)
(a) Determine the reduced observer gain matrix Lc assume that desired observer
pole is at s = −5. (4 points)
(b) Determine the reduced-order state observer dynamics and reduced-order state
observer output equation(open loop system dynamics).(3 points)
(c) Draw the overall block diagram for the system with reduced state observer.
(1 points)
4. For the system shown on equation (B.2.13) (10 points) If the designer wants to
optimize the performance index J,
1 ∞ ⊤
Z
J= (x Qx + u⊤ Ru) dt
2 0
where
1 0
Q= , R=1
0 1
xxiii
Appendix C
The 3/2 solenoid actuator valve, depicted in Figure C.1, plays a crucial role in regulating
fluid flow through three ports: one inlet and two outlets. When an electrical current is
Spool Valve
Air gap sping
Seated Position Coil Supply Drain
Plunger
Drain
applied to the solenoid, it generates a magnetic field that actuates the internal mechanism
of the valve. This movement facilitates fluid flow from the inlet to one outlet while
simultaneously blocking the other. Upon deactivation of the current, a spring mechanism
restores the valve to its original position, thereby redirecting the flow to the second
outlet. This efficient operation enables precise control over actuators and cylinders in
various pneumatic and hydraulic applications, establishing the solenoid actuator valve as
a vital component in automation and fluid management systems.
Since the spool valve operates within a fluid medium, which introduces viscous friction
that influences its dynamic behavior. The free body diagram of the system, illustrated
in Figure C.2, elucidates the forces acting on the valve and provides insights into its
dynamics, particularly regarding the impact of fluid resistance on its operation.
xxiv
R
x
i k
Fem
vs L iL m
d
(a) Electrical Subsystem (b) Mechanical Subsystem
(a) Derive the transfer function of the system from the input (source voltage Vs (s))
to the output (Spool Valve position X(s)) in zero-pole-gain form (zpk form).
(3 points)
(b) Derive the transfer function, poles, and zeros of the system. (3 points)
(c) Determine the DC gain of the system. Which parameter affects the DC gain
of the system? What conclusions can you draw from your result? (2 points)
(d) Approximate the system to second order using dominant poles approximation.
(Note:- The DC gain of the original system transfer function (1a) and reduced
order (approximated) system should be equal.) (5 points)
(e) Which system parameter affects the steady-state error (ess ) for a step input?
(2 points)
(a) Does it affect the steady-state error of the system for a step input? Compare
your results with (1d) and draw a conclusion. (3 points)
(b) Does it affect the DC gain of the system? Compare your results with (1c) and
draw a conclusion. (2 points)
xxv
C.1.2 Test 2 (2025 G.C)
A shock absorber reduces unwanted oscillations in mechanical systems, enhancing sta-
bility and performance. It works by absorbing and dissipating energy from vibrations
through a combination of hydraulic fluid and mechanical components. When a force
is applied, the shock absorber compresses, allowing the fluid to flow through valves,
which dampens the motion. The stiffness and damping coefficients can vary based on
application, design, and specific requirements, enabling tailored performance for different
environments, such as vehicles or industrial machinery.
The 3D model and Mechanical System Model of a typical shock absorber are shown
in Figures C.3a and C.3b, respectively.
moving piston
fluid
SECTION A-A x
spring
d
k m u
A
(a) Derive the transfer function, poles, and zeros of the system. (3 points)
(b) Determine the DC gain of the system. Which parameter affects the DC gain
of the system? What conclusions can you draw from your result? (2 points)
(c) Determine the conditions for the system to be undamped, under-damped,
critically damped, and overdamped. How does each scenario influence the
system’s poles and its response to inputs?(5 points)
(d) Which system parameter affects the steady-state error (ess ) for a step input?
(2 points)
xxvi
(c) What would be the effect on peak time (tp )? (2 points)
(a) Does it affect the steady-state error of the system for a step input? Compare
your results with (1d) and draw a conclusion. (3 points)
(b) Does it affect the DC gain of the system? Compare your results with (1b) and
draw a conclusion. (2 points)
xxvii
C.2 Control Systems Exams Final
C.2.1 Final Exam 1 (2025 G.C.)
Solenoid Actuator Valve System
The operational parameters of the solenoid actuator valve system in Figure C.1 are given
in Table C.2, including the parameters and voltage gain of the non-inverting operational
amplifier configuration.
Parameter Value
Stiffness Coefficient (K) 100N/m
Number of Turn (N ) 400
Winding Resistance (R) 1Ω
Winding Inductance (L) 0.5H
Spool valve mass (m) 1kg
Damping Friction Coefficient (b) 25N s/m
Electromagnetic Force Fem = N iL
1. For the solenoid actuator valve system shown in Figure C.1; (12 points)
(a) Derive the transfer function of the system from the input (source voltage Vs (s))
to the output (Spool Valve position X(s)) in zero-pole-gain form (zpk form).
(8 points)
(b) Draw a block diagram representation of the system.
(c) Approximate the system to second order using dominant poles approxima-
tion.(4 points)
(a) Determine the time constants for the electrical and mechanical subsystems.
What conclusions can you draw from your results about part 1a? (5 points)
(b) Determine the DC gain of the system. (2 points)
(c) Determine the damping factor (ξ) and natural frequency (ωn ) of the system.
(5 points)
(d) Is the system Stable? Why? (3 points)
(e) Determine the steady-state error (ess ) of the system for a step input. (2 points)
(f) Determine the rise time, maximum percentage overshoot, and steady-state
error of the system with unity feedback for a step input.(5 points)
xxviii
(g) For the closed-loop frequency response, determine the resonant frequency, the
magnitude of the system at the resonant frequency, and the phase angle at the
resonant frequency. What implications does your answer have for the closed-
loop time-domain analysis of the system? Relate your answer to part 2f. (6
points)
(h) Determine the response of the closed loop system for r(t) = 2 cos(4t) input.
(4 points)
(i) Determine the gain crossover frequency (ωgc ), phase crossover frequency (ωpc ),
gain margin (GM ), and phase margin (P M ) of the open loop transfer function
of the system and interpret the implications of the results. (7 points)
3. Design a lead compensator to ensure that the steady-state error for a step input is
less than 0.02, while achieving a desired phase margin of at least 45◦ at 20rad/s
frequency for the system you derived in part 1b. (20 points)
+ ia
Ra La
+
− +
Rf ωm
Vs Va M
Ri
-
-
(a) Armature Voltage Controlled DC Mo-
tor With Pre-amplifier (b) DC-Fan
In this configuration, the weight of the fan is considered negligible (Zero load torque),
which contributes to optimal efficiency. The motor is directly coupled to the fan blades,
ensuring smooth and responsive operation. The pre-amplifier is instrumental in stepping
up the input voltage, allowing for refined control over the motor’s speed and improving
the overall effectiveness of the system.
The operational parameters of the DC motor are given in Table C.2, including the pa-
rameters and voltage gain of the non-inverting operational amplifier configuration.
xxix
Table C.2: DC Motor Parameters
Parameter Value
Back EMF constant (Kb ) 0.001 V-s/rad
Torque constant (Kt ) 0.001 V-s/rad
Winding Resistance (Ra ) 1Ω
Winding Inductance (La ) 0.5 H
Motor Inertia (Jm ) 0.001kg/m2
Windage friction (b) 0.01 N-s
Input resistance (Ri ) 1 kΩ
Feedback Resistance (Rf ) 470 kΩ
Rf
Non-Inverting Op-amp voltage gain (Av ) 1+
Ri
1. For the DC fan speed control system shown in Figure C.4; (30 points)
(a) Derive the transfer function of the system from the input (source voltage Vs (s))
to the output (motor speed ωm (s)) in zero-pole-gain form (zpk form). (8
points)
(b) Draw a block diagram representation of the system.
(c) Determine the time constants for the electrical and mechanical subsystems.
What conclusions can you draw from your results about part 1a? (5 points)
(d) Determine the DC gain of the system. (2 points)
(e) Determine the damping factor (ξ) and natural frequency (ωn ) of the system.
(5 points)
(f) Is the system Stable? Why? (3 points)
(g) Determine the steady-state error (ess ) of the system for a step input. (2 points)
(h) Determine the settling time, maximum percentage overshoot, and steady-state
error of the system with unity feedback for a step input.(5 points)
(a) For the closed-loop frequency response, determine the resonant frequency, the
magnitude of the system at the resonant frequency, and the phase angle at the
resonant frequency. What implications does your answer have for the closed-
loop time-domain analysis of the system? Relate your answer to part 2d. (4
points)
(b) Determine the response of the closed loop system for r(t) = 2 cos(4t) input.
(4 points)
(c) Determine the gain crossover frequency (ωgc ), phase crossover frequency (ωpc ),
gain margin (GM ), and phase margin (P M ) of the open loop transfer function
of the system and interpret the implications of the results. (8 points)
xxx
3. Design a lead compensator to ensure that the steady-state error for a step input is
less than 0.02, while achieving a desired phase margin of at least 45◦ . (20 points)
Consider a projector screen positioning system in Figure 1.18, using the operational
parameters of the DC motor detailed in Table C.2, which includes specifications and the
voltage gain of the non-inverting operational amplifier configuration.
1. For the Projector Screen Positioning System shown in Figure C.4; (30 points)
(a) Derive the transfer function of the system from the input (source voltage Vs (s))
to the output (motor speed ωm (s)) in zero-pole-gain form (zpk form). (8
points)
(b) Approximate the system to second order using dominant poles approximation.
(4 points)
(c) Determine the DC gain of the system (1a) and the reduced order system (1b).
(2 points)
(d) Determine the steady-state error (ess ) of both systems for a step input. (2
points)
(a) Determine the time constants for the electrical and mechanical subsystems.
What conclusions can you draw from your results about part 1b? (5 points)
(b) Determine the damping factor (ξ) and natural frequency (ωn ) of the reduced
order system. (5 points)
(c) Is the reduced order system Stable? Why? (3 points)
(d) Determine the settling time, maximum percentage overshoot, and steady-state
error of the reduced-order system with unity feedback for a step input.(5
points)
(e) For the closed-loop frequency response, determine the resonant frequency, the
magnitude of the system at the resonant frequency, and the phase angle at the
resonant frequency. What implications does your answer have for the closed-
loop time-domain analysis of the reduced-order system? Relate your answer
to part 2d. (4 points)
(f) Determine the gain crossover frequency (ωgc ), phase crossover frequency (ωpc ),
gain margin (GM ), and phase margin (P M ) of the open loop transfer function
of the reduced order system and interpret the implication of the results. (8
points)
3. Design a lead compensator to ensure that the steady-state error for a step input
is less than 0.02, while achieving a desired phase margin of at least 45◦ for the
reduced-order system in (1b). (20 points)
xxxi
Appendix D
xxxii
(d) If the system is not stable the design a controller using the direct Raggazzini
method with one pole of the closed loop transfer function at z = 0.5. (2 marks)
3. Given a sampling period T = 0.5 sec. determine the pulse transfer function for,
R(S) C(S)
Figure D.1: a
R(S) C(S)
Figure D.2: b
xxxiii
D.1.2 Final Exam 2
1. Batter making system has two inputs called water and flour. The system mixes
the two inputs and results in one output called batter. The discrete-time equation
governing the relation between the batter and flour in the discrete time domain is
as follows. y(k + 2) + 5y(k + 1) + 6y(k) = u1 (k + 1) + 2u1 (k)
(a) Derive a controllable canonical form state space representation for the system
(5 points)
(b) Derive diagonal canonical form state space representation for the system (5
points)
(c) Determine discrete time state space transition matrix. (7 points)
(d) For constant flour input 1kg/min determines the solution for the states vector
in the discrete time domain and the output equation in the discrete time
domain. (7 points)
2. For system in question 1, if the desired location of the poles is at z1 = −0.5 and
z2 = 0.4
3. Derive an expression for diagonal controllable canonical form, assuming that the
system has distinct poles p1 , p2 , ... pn . (5 points)
4. Derive an expression for the state transition matrix using the z-transform method.
( 4 points)
xxxiv
2. For the system derived in 1a,
(a) Determine the state feedback gain Kf that locates the poles of the discrete-
time system at s1 = −3 and s2 = −0.5. (5 points)
(b) Determine the predictive estimator gain Lp that locates the poles of the esti-
mator at s1 = −3 and s2 = −0.5. (5 points)
(c) Determine the current estimator gain Lc that locates the poles of the estimator
at s1 = −3 and s2 = −0.5. (6 points)
(d) Compare the results from 2b and 2c and derive a conclusion. (3 points)
3. If a PID controller is used to control the plant derived in 1a. Determine the
discrete-time equivalent of the controller. (5 points)
xxxv
"Strive for perfection in everything you do. Take the best that exists
and make it better. When it does not exist, design it."
-Sir Henry Royce, English engineer and car designer.
Remember always! We are yesterday, and you are tomorrow; work hard before you end
up being yesterday!
Good luck and be happy being an engineer!
I am very happy to address or clarify any confusing concepts or anything that feels
wrong. If you feel so, contact me using email: abebe.daniel7@[Link]
Daniel Abebe Beyene (Science Ambassador, MSc in Control Engineering)
Lecturer at AASTU, Electro-mechanical Engineering Department
Thank you very much!!!
xxxvi
References
xxxvii
Poles of a reduced-order state observer are chosen based on desired response characteristics such as speed and damping. The observer's characteristic equation should match the desired polynomial that specifies pole locations. These serve to optimize the observer's convergence speed and stability, ensuring the observer states estimate the actual states effectively under closed-loop dynamics .
The Riccati equation is solved to find a matrix P that defines the optimal state feedback gain for an LQR. This gain minimizes a quadratic cost function that balances state error and control effort. The Riccati equation links system dynamics with performance weighting matrices, yielding a feedback scheme that refines system control towards optimal performance .
In optimal control, the state weighting matrix Q and input weighting matrix R in the cost function are chosen based on design priorities. A higher weight in Q prioritizes reducing state errors, whereas a higher weight in R emphasizes minimizing control effort. The precise values depend on the system’s performance specifications and desired trade-offs between state deviation and control energy .
The controllability matrix, Co, assesses if all state variables of a system can be controlled by manipulating inputs. If the rank of Co is equal to the number of state variables, the system is fully controllable. Similarly, the observability matrix determines if the state variables can be observed through the outputs. A system is completely observable if the observability matrix has full rank equivalent to the number of states. These matrices help identify uncontrollable or unobservable subsystems and influence the design of controllers and observers .
The duality property asserts that the mathematical structure of linear quadratic regulator (LQR) and linear quadratic estimator (LQE) problems is equivalent, allowing the optimal estimator gain matrix Lc to be derived by transposing and appropriately scaling the LQR's feedback gain matrix. This relationship simplifies the design process by providing a direct connection between control and estimation solutions .
A shock absorber system shifts from under-damped to critically damped when its damping ratio increases from less than one to exactly one. At under-damping, poles have a significant imaginary component, leading to oscillations. At critical damping, poles are real and repeated, minimizing oscillations and achieving the fastest convergence without overshoot .
Similarity transformations utilize reversible linear transformations to reorganize system state-space matrices without altering the inherent transfer function or dynamic characteristics. These transformations simplify system analysis by indicating equivalences such as diagonalization or transformation to canonical forms, fostering easier application of control theory .
A dynamic system is stabilizable if the uncontrollable states do not destabilize the system, implying that their eigenvalues have negative real parts. These eigenvalues indicate how system dynamics evolve over time and are crucial for determining if the system can be made stable through partial control of the controllable states .
The transformation matrix facilitates the decoupling of a system into controllable and uncontrollable subspaces by applying changes of basis that isolate these subspaces. With the matrix, parts of the system’s dynamic matrices are reordered and transformed into canonical forms, simplifying the analysis and design of controllers and observers for each subspace separately .
Unity negative feedback reduces the steady-state error by adjusting the system's output closer to the reference input, improving accuracy and response precision. It modifies the loop gain, altering the error sensitivity while enhancing system stability .