Systems Theory in Automation Explained
Systems Theory in Automation Explained
The present notes introduce the study of stationary linear dynamic systems.
Classes of models that can be assumed to represent, possibly in
approximate way, the behaviors of various phenomena and natural processes
and artificial. The proposed analysis methods are preparatory to the techniques of
design of automatic systems and control.
Salvatore MONACO
professor of Systems Theory
Department of Computer, Automation and Management Engineering 'Antonio
Ruberti
University of Rome 'La Sapienza'
via Ariosto 25 - 00185 Rome
[Link]: [Link]@[Link]
2. System Theory in Automation
.. in Automation
Planning
supervision
intervention
Figure 2.1
A decision-making and intervention process that does not consider this flow
the inverse of information would imply a perfect correspondence between
planning and execution of the task which is not possible in the presence of
uncertainties and disturbances. In reality, in fact, the knowledge of the way of
the functioning of the process is not free from uncertainties (model uncertainties) and
Interaction with the environment can give rise to unforeseen disturbances.
It is therefore necessary to foresee that during normal operation
the decision-making and intervention process should be able to adapt, with elaborations
ad hoc that use measures of the deviation between expected behavior and
the effective feedback that has been mentioned.
This is represented in the diagram of figure (1.1) where inside each
stage, and in the connection between them, the interventions (u(t)) are activated not only
from the information of the descending flow, v(t) and z(t), but also from information
in feedback, m(t) (feedback). To better understand the three-tier structure
the characteristic of an automatic system, think of a robotic system for
the exploration of a planet. Imagine a mobile platform that needs to
carry out excursions, collect soil samples, capture images, carry out
elementary analysis on site and transmit the results... In accordance with the intuition
we can imagine the completion of an exploration mission like the
result of a complex procedure composed of the three indicated phases.
6 2. System Theory in Automation
An initial planning will have to define the path that the robot will have to follow.
seguire e il sequenziamento delle operazioni; le elaborazioni necessarie saranno
activated by external inputs that represent the desired objectives (number of
samples to be taken, experiments, ..), and will take into account the characteristics
of the soil detected through images or other types of data. This corresponds
to the first phase of the proposed plan. The path and sequencing have been planned.
During the operation, the devices of the robotic system must ensure
to the organization of the sequencing of activities, to activate and coordinate them
the execution; what corresponds to the second and third phase of the scheme with-
proposed settlement. The close interaction in the execution of the
verses phases and the need to reactivate processes based on the results that
Obstacles and unforeseen barriers are detected along the way.
It may be necessary at a high level to redefine the trajectory or to replan.
that the execution of scheduled soil survey experiments, as well as
the malfunction of an actuator or the modification of the characteristics of
soil makes it necessary to modify the execution interventions at the level
lower according to automatic modes. In this context, and with reference-
related to the architecture indicated, in the teachings of Automatic Controls
The Foundations of Automation specifically study the methods that
they are at the core of the execution level project. The methodological approach
of the project, and many of the concepts used are common to the three phases, that is
gives the study characteristics of adequate generality.
C P m(t)
Figure 2.2
What basic examples come to mind regarding the problem of maintaining a preset?
temperature profile in an oven or refrigerator, compensating for the effect of
8 2. System Theory in Automation
Since the models of the system and the various components are an approximation
of physical reality and do not describe with absolute precision the links between the vari-
abilities, the control subsystem calculated from this schema must
garantire il mantenimento delle modalita` di comportamento richieste entro
predefine variability margins of component models. This pro-
the robustness of the control system is indicated as necessary
would ensure the fulfillment of the project specifications of the system
physical.
of any measurable disturbances and, on the other hand, of the reference quantities
v, which represent the external commands.
The theory of systems and control encompasses a wide range of
project methodologies with reference to predefined families of models as-
it's used to represent the plant or process. If the plant is described by a...
Linear systems of differential equations are well documented.
synthesis methodologies of controllers. In many cases, it is possible to refer to
to these methodologies even if the process is not linear, using its
approximate representation.
The main analysis methods for this class of systems are the subject of
queste note preliminari allo studio dei metodi dell’automatica.
An implicit representation with the state of a system of dimension
finished, adopting the vector notation, is of the type
∆x(t) = f(x(t), u(t), z(t)) m(t) = w(x(t), z(t)) y(t) = h(x(t), z(t)) (1)
in which the time can be assumed continuous (coinciding with the set
of real numbers) or discrete (coinciding with the set of relative numbers); in
d(x(t)
correspondence will be the derivative operator dt or the operator of
unitary advance x(t+!). The behavior of the system with respect to time t is
therefore described by five types of variables: u(.) is the vector of the variables of
control, they represent the channels through which it is possible to modify
the behavior of the process; z(.) that of the disturbances, magnitudes that influence
they characterize the process behavior, they are sometimes measurable, but they cannot
be modified; m(.) represents the measured quantities; y(.) the outputs of
process, these are a subset of the measured variables and represent the
quantities on which control is to be exercised; x(.), finally, are the variables
of the state and represent, with their evolution over time, the behavior
of the system in its most complete form. The function f, the generating function,
h, the output transformation, ew, the transformation of the measures, have a
structure corresponding to classes of processes and depend on parameters that
we will support notes in an initial phase. If, for example, such functions are linear
11
In this chapter, we want to show, with the help of simple examples, that
the behavior of numerous physical systems, phenomena and processes, from different
disciplinary sectors can be described by systems of differential equations,
continuous-time systems, or prime difference systems, discrete-time systems,
type
The interaction with the environment, which is realized through the inputs, u(t),
and the outputs, y(t), give life to a cause-and-effect relationship that is expressed completely
using a set of auxiliary variables, the state variables.
In the hypothesis of linearity of the functions: Rn xRp→Rn edh:Rn xRp→
Rq on their respective domains of definition, the representation with the state (1)
take on the form
x
Figure 3.1
deriving and taking into account the balance of forces, one obtains
ẋ(t)
1 = x2 (t)
ẋ(t)
2 = -k x1 (t) + -b x2 (t) + u(t)
M M M
y(t) = x 1 (t)
0 1 0
A= ! -k -b " B= ! 1 " C= ( 1 0 )
M M M
16 3. Representations with the state of finite-dimensional linear stationary systems
u R2 y
C
- -
Figure 3.2
u(t) − vc (t)
u(t) = R1 I1 (t) +vc (t)=⇒i1 (t) =
R1
from the sum of the currents at the node
to
u(t) = Ri(t) + ec (t) +L
dt
dovec (t) = kc ω(t). The electrical and mechanical transduction occurs for generating
action of a couplem , of intensity proportional to the intensity of the current
i(t)
Cm (t) = km i(t)
the dynamics of the mechanical part is described by
dω(t)
J +Fω(t) = Cm (t)
dt
It is ultimately obtained
R kc 1
i̇(t) = -i(t) - ω(t) + u(t)
L L L
km F
=)t(ω̇ i(t) - ω(t)
J J
We ultimately obtain a stationary linear system of size
two with statox(t) = (i(t) ω(t) ′ .
R Kc 1
A= ! − L − L
B= ! L C= (O1 )
km -F " 0"
J J
3.2. Equivalent representations: continuous-time systems 19
u I Cm ω
EE km
m
M
kc
c
Figure 3.3
A
=
+
B)xutẋ(
y(t) =Cx(t) +Du(t)
with a linear transformation of state variables through a matrix T
non-singular constant
z=Tx |T|≠0
one obtains the representation
ż()=
tTẋ()=
tTAx()+ tTAT−1 z(t) + TBu(t)
tTBu()=
there is still a system of the same type with dynamic matrix TAT−1 , ma-
triple of inputs TB, and matrix of outputs CT−1 The student will not have
difficulty in applying the previous considerations to some of the previous examples
teeth.
An important observation that helps in the calculation of a transformation of
The state variable is the following: if you remember that x is a n-tuple that represents
a vector, siav, with respect to a reference frame, that is
3.3. Equivalent models of different physical systems 21
0
1 0
⎛ 1⎞
⎛ 0⎞ ⎛ ⎞...
x = Σxi eI= . x1+ ⎜ 0⎟ x2+...+ xn
⎜ ⎟. ⎜ ⎟.. ⎜ 0⎟
⎝.⎠
⎜ 0⎟ ⎜ ⎟ ⎜ 1⎟
⎝.⎠
⎜ 0⎟ ⎝ ⎠
people
x=T−1 z
0
⎛ .. ⎞
.
⎜ 0⎟
⎜ ⎟
(T−1 )i=T−1 · ⎜⎜ 1 ⎟ i−ma
⎟
⎜ 0⎟
⎜ .. ⎟⎟
⎜ .⎟
⎜
⎜ 0 ⎟⎠
⎝
(T−1 )i , the i-th column of (T −1 and the representation of the i-th vet-
tower of the new base (the one in which the vector is represented by z) with respect to
to the old base (the one in which the vector is represented by dax). This means
to perform a linear transformation of the state variable of the type
z=T means, from an operational point of view, choosing a new basis.
∼ ∼
of the state space that is represented by1 , ..., e n in the coordinates ten;
that is
∼ ∼
T −1= (e 1 , ..., e n )
+
b)ut)=
(axt()tẋ( z=ex x=lnz
ż()=
t ex (ax(t) + bu(t)) = az(t)lnz(t) + bz(t)u(t)
22 3. Representations with the state of finite-dimensional linear stationary systems
and clarify what has been stated, namely that there are non-linear representations
of linear systems.
1 dT(t) q(t)
q(t) = (T2 (t)−T1 (t)) =
R dt C
trovano una corrispondenza immediata nelle equivalenti elettriche
1 dv(t) i(t)
i(t) = (v2 (t)−v1 (t)) =
R dt C
where it is assumed that they represent resistance and thermal capacity
electric power, expressed in degree/calorie and calorie per degree, ohm farad,
respectively.
With the established correspondences (and0 (t) correspond to aT0 (t)) it's not difficult
check that the following equations define a different representation
circuit branch valid for the thermal system
3.4. Linear approximations of nonlinear systems: continuous time systems 23
1 1 1 1 1
v̇(t)
1 =− ( + )v1 (t) + v2 (t) + e0 (t)
C1R1 R2 C1 R3 R1 C1
1 1 1 1 1
v̇(t)
2 = v1 (t)− ( + )v2 (t) e0 (t)
C2 R3 C2R2 R3 C2 R2
df
y = f(x) = f(xe ) + (x−xe ) +....
dx ))xe
)
)
df
the postoxa=x−x e ,ya=y−y e ,m= dx a linear relationship is obtained
next that approximates the curve around adxe
ya=mxa .
u)x(=
,fẋ f ( x e , ue ) = 0
obtaining:
24 3. Representations with the state of finite-dimensional stationary linear systems
∂f∂f x−xe
f(x, u) = f(xe , ue ) + ! ∂x∂u " ) ! +...
)xe ,ue u−ue "
)
)
∂f ∂f
=f(xe , ue ) + (x−x ) + (u−ue ) +...
∂x ))xe ,ue ∂u ))xe ,ue
e
) )
and place ) )
z = x - xe , v = u - ue
∂f ∂f
A= , B=
∂x ))xe ,ue ∂x ))xe ,ue
) )
if so: ) )
ż=
()tAz+
()tBv+
()t.
In the same way for the exit
∂h∂h x−xe
h(x, u) = h(xe , ue ) + ! ∂x∂u " ) ! +...
)xe ,ue u−ue "
)
)
∂h ∂h
=hand+ (x−xe ) + (u−ue ) +...
∂x ))xe ,ue ∂u ))xe ,ue
) )
the place ) )
∂h ∂h
C= , D=
∂x ))xe ,ue ∂x ))xe ,ue
) )
)
the linear approximation is obtained )
This is the linearized model that approximates around the pair (xe , ue ) he
∂f ∂f they are the derivatives of the vector
behavior of the nonlinear system. ∂x e ∂u
compared to the state vector and the output vector, respectively.
If for the first, called the Jacobian matrix,
∂f 1 ∂f 1
... ∂x
∂f ∂x 1
= ⎛ ⎞
n
J(x) = ..
∂x ∂f n ∂f n
⎝ ∂x 1 ... ∂x n ⎠
3.4. Linear approximations of nonlinear systems: continuous-time systems 25
Pendulum
Si consideri un pendolo di massam, sospeso ad un’asta rigida di peso
negligible in length on which it acts according to the tangent to the motion a
external force. A mathematical model is to be written that represents
feel the movement on the horizontal, originating from the equilibrium point, of the
projection along the vertical of the mass of the pendulum. The balance of forces
along the tangent to the motion allows us to write
Where represents the external force acting, and it is an efficient coefficient of attraction.
˙ concerning the
dynamic. Setting x1 (t) = θ(t) ex2 (t) = θ(t) if it is,
equations of state,
ẋ(t)
1 = x2 (t)
k
-g thank you1 (t)−x 1
ẋ(t)
2 = 2 (t) + u(t)
l m ml
that is a nonlinear representation of the type
26 3. Representations with the state of linear, stationary systems of finite dimension
0
+)xt(=
f)tẋ( ! 1 " u(t)
ml
in which the output is equal to
ẋ(t)
1 = x2 (t)
g k 1
ẋ(t)
2 = x1 (t)−x2 (t) + u(t)
l m ml
y(t) = lx 1 (t)
for others:
ẋ(t)
1 = x2 (t)
g k 1
dx2ot(t) = -x1 (t)−x2 (t) + u(t)
l m ml
y(t) = lx 1 (t)
The following model describes the dynamics of two interacting species, prey and
predator, it was introduced in 1926 by the Italian mathematician Vito Volterra and
takes its name. This model is able to interpret the most salient aspect.
of such a phenomenon: it concerns the presence of equilibrium situations that, if
perturbate, they see the emergence of oscillation phenomena (alternation of development
among the species).
It is assumed, in the formulation of the model, that:
The prey grows, in the absence of predation, according to the so-called equation
logistics
3.5. Further examples of continuous-time systems* 27
ẋ(t)
2 = -cx2 (t)
the predation rate is proportional to the product of the individuals of
two species
With these assumptions, the equations that describe the evolution of species
can be easily deduced from the previous equations. It is obtained:
ẋ(t)
1 = ax1 (t)−bx1 (t)x2 (t)
dx2ot(t) = -cx2 (t) +dx1 (t)x2 (t)
y1 (t) = x1 (t) y2 (t) = x2 (t)
So once again there is a nonlinear model of the same type
from the previous one. It is easy to verify that in this case there are two points of
equilibrium: one trivial, corresponding to the absence of species, the other
c
x2 (-c+dx1 = 0 → x1e =
d
ad−kc
x1 (a−kx1−bx2 0 = x2e =
bd
Around such a point of equilibrium, the evolution is approximated by a sis-
linear theme in free evolution with
kc bc
A= ! − d − d
ad−kc
b 0 "
28 3. Representations with the state of finite-dimensional stationary linear systems
iC (t) = il (t)−iT(t)
The current flow in the diode depends on the voltage across its terminals.
according to a nonlinear relationship, IT =h(vTreported in the graph of Figure
6. From the previous equalities, with obvious steps, we obtain
R L
iT
vC
C
Figure 3.4
dvc (t)
C =il (t)−h(vC (t)
dt
ofl (t)
L =−vC(t)−Ril (t) + u(t)
dt
With positions similar to the previous ones
h(x1 )x2
0 =− +
C C
x1 R 1
0 =− x2+ U
L − L L
Solving with respect to x2one obtains the equality
U 1
x1=h(x1 )
R − R
which for fixed values of UedRammette has one or three solutions as it turns out
evident from Figure 6
iR
e1
e2
e3
vR
Figure 3.5
The circuit under examination from an experimental point of view shows a behavior
so-called bi-stable actuator on the two equilibrium states 1ede3in response to
sufficiently large impulsive stresses.
How it can be verified based on the theory outlined in chapter 7,
such behavior depends on the properties of the equilibrium states resulting
e1 ede3establish2unstable.
30 3. Representations with the state of finite-dimensional linear stationary systems
Reserving the right to resume the study with one of the investigation techniques
what we will study, the trend of effective behavior can be verified
cat through numerical simulation. To this end, one can assume the
following parameter values: u = 1.2V, R = 1.5Kohm, C = 2pF, L = 5µH
scaling the time to nanoseconds and measuring the currents in
milliamps, the circuit equations become
ẋ(t)
2 = 0.2(-x1-1.5x2+ 1.2)
Perh(x1 it can be assumed
3 5
h(x1 ) = 17.76x1-103.79x1 2+ 229.62x1-226.31x1 4+ 83.72x1
feky
k(1−a2 y 2 )y |ay|<1
which models a decrease in elastic force as the displacement increases
corresponding to a failure of the elastic characteristics; or
a relationship of the type
k(1 + a2 y 2 )y
external forces the mass remains stationary until when|f e |µs mg. In motion the
friction force takes the value of µk mg, µk dynamic friction coefficient.
the conclusion is reached
ẋ=x
1 2
ẋ=x
1 2
ẋ=x
1 2
ẋ=
2
-k x1 (t) + -b x2 (t) + u(t) +µk g
M M M
3.5. Further examples of continuous-time systems* 33
u1
r u2
ϑ
Figure 3.6
˙2 β
r¨(t) = r(t)θ(t)− +u1 (t)
r(t) 2
ṙ()tθ()tu2 (t)˙
θ̈(t) = -2 +
r(t) r(t)
In the absence of radial thrust1and the ring road2the solutions are ellipses,
hyperbolas or parabolas. The simplest orbit is a circle
r(t) = c 1 θ(t) = c 2
1
Conr(0) = r0 r˙(0) = 0, θ(0) = θ0 ,θ(0) =ω0˙eω0= ( )
β 2
r30
The nominal orbit is thus described by
˙
Postox1=r,x2=r˙,x3=θ,x4the following representation is obtained
station with the state
ẋ(t)
1 = x2 (t)
β
ẋ(t)
2 = x1 (t)x42(t)− +u1 (t)
x12(t)
ẋ(t)
3 = x4 (t)
x2 (t)x4 t u2 (t)
ẋ(t)
4 = -2 +
x1 (t) x1 (t)
The linearization of the dynamics, that is, the description of the deviations
compared to the nominal orbit
r0
0 0 ⎛ ⎞
ur (t) = ! " x(t)
r =
0 ω0 t+θ0
⎜ ⎟
⎝ ω0 ⎠
from the following stationary linear representation
0 1 0 0 0 0
⎛ 3ω20 0 0 2r0 ω0 ⎞ ⎛ 1 0⎞
ż()=
t z(t) + u(t)
0 0 0 1 0 0
⎜ ⎟ ⎜ ⎟
⎝ 0 -2ωr
0
0
0 0 ⎠ ⎝ 0r0 ⎠
h0
Figure 3.7
3.5. Further examples of continuous-time systems* 35
x2
⎛ ve u(t) ⎞
=)tẋ(-g+ x3
⎝ u(t) ⎠
.
In the particular case, very frequent in practice, where one assumes a
mass variation, u0 constant, that which corresponds to assuming that
m(t) = m0+u0 t
it is obtained
m(t)v˙(t) = -m(t)g + ve u0
and, postox1=h,x2=v
0 1 0
=)tẋ( ! x(t) + !
0" "
ve u0
0 - (m0 +u0 t)
y(t) = ( 1 0 )x(t)x(0) = 0
0 1 0 0 0
ż()=
t ⎛ 0 0−(m+you teav)e2uz(t)+
0 ⎞ ⎛ ve ⎞ v(t)z(0) = x(0) - ⎛ 0 ⎞
0 0 (m0 +u0 t)
⎝ 0 0 0 ⎠ ⎝ 1 ⎠ ⎝ m0 ⎠
iR
Rs
ω Ls vs
Figure 3.8
One wants to calculate the mathematical model that describes the dynamics of
motore elettrico in corrente continua di figura. Si assume che l’ingresso sia la
excitation circuit voltage, u(t) = vs (t). For the excitation circuit
the following equation holds for the loop that expresses the balance of tensions
ofs (t)
Vs (t) = Rs Is (t) +Ls
dt
The generated flow can be expressed by
Φ(t) = Ls is (t)
For the armature circuit, the following equation holds that expresses the equilibrium
of the tensions:
ofr (t)
Vr (t)−ec (t) = Rr ir (t) +Lr
dt
in which the electromotive force is equal to
ec (t) = kc Φ(t)ω(t).
Finally, regarding the mechanical part, there is
3.5. Further examples of continuous-time systems* 37
dω(t)
J +Fω(t) = Cm (t)
dt
in which the drive pair, Cm , and the date of
Rs 1
ẋ(t)
1 = -x1 (t) + u(t)
Ls Ls
Rr Vr KL s
dx2ot(t) = -x2 (t) + − x1 (t)x3 (t)
Lr Lr Lr
F KL s
dx3ot(t) = -x3 (t) + x1 (t)x2 (t)
J J
Figure 3.9
Figure 3.10
writes
Mx¨(t) = f(x(t), i(t)) − Mg
which indicates the acceleration due to gravity.
Around axleethe linear approximation perf(x, i)
∂f ∂f
f(x, i) = f(xe , ie ) + (xe , ie )(x−xe ) + (xe , Ie )(i−ie )
∂x ∂i
The derivatives with respect to ax,kx , and adi,ki they are calculated in the following [Link] the
slope at constant current equal to adieat the point of x-axis e ; in the
the case in question is equal to 14 N/m. kicorresponds to the variation in strength,
inxe with respect to the current; in the case at hand peri1= 700mA the force is equal
a 12210−3Don't worry3= 500mA The force is equal to 4210−3N. So it has therefore
80−3
k i= = 0.4N/A
700−500
By substituting the numerical values, the following differential equation is obtained,
which governs the motion around the equilibrium (ζ=x−x eeu=δi)
¨ = 14ζ(t) + 0.4δi(t)
(3210−3)ζ(t) ζ̈(t) = 1667ζ(t) + 47.6u(t)
v tn=ω0
vn =
v0 t
If it is observed that the effect of scaling with respect to time means that
ζn (t)
(ω20 ζ0 d2 ) = (1667ζ0)ζn (t) + (47.6u0 )u(t)
dtn2
that is
d2 ζn (t) 1667 47.6u0
2 = 2
ζn + u(t)
dtn ω0 ω20 ζ0
Choosingω20= 1667 (ω0 approximately equal to 40) the coefficient of ζ(t) is even
1
to one that corresponds to taking a timeline equal to 40 seconds.
If moreover the displacements are measured in cm and the currents in Amperes,
that is chosen0= 0.01 eu0= 1, the equation is obtained
d2 ζn (t)
=z(t) + 2.86u(t)
dtn 2
Set in scale
Regarding the general case of a system of differential equations
first-order linear in the state variable x ∈ R nand entered∈Rp ,
the scaling operation leads to the following system
ẋ(t)
n = 1/ω0 S−1 −1
x ASx+ 1 over ω0 Sx BS u un
in which the real number ω0represents the variation of the temporal scale and Sx
edSuthey are square diagonal matrices, (n×n) and (p×p), of the coefficients of
scalature associated with transformation
3.6. Discrete time systems 41
xn=S−1
x x un=S−1
u u.
If the time variable is discrete and under the assumption of linearity of the function
generator of the product space (X×U) the representation takes on the
shape.
D
u (t ) x (t+ 1) x(t ) +
1
B + z C + and
( )t
+
A
Figure 3.11
figure 2.13; in that case the discrete model approximately describes the
behavior of the continuous-time system championship.
H S
TC
u y
Clock
S TD
Figure 3.12
A
=
+
B)xutẋ(
. x(t0 =x0
y(t) = Cx(t)
t t0=kT
x(t) = eA(t−t ) x00+ / eA(t−τ ) Bu(τ)dτ
t0 t = (k + 1)T
(k+1)T
x(k+1)T=e ATx(kT)+ / eA(k+1)T−τ
+ dτBu(kT)
, (k+1)T - τ = ξ
+ , kT
T
x(k + 1)T = eATx(kT) + / eAξ dξBu(k)
+ , 0
ove
T
AD=eAT BD= / eAξ dξB CD=C
0
Some considerations on the calculation of the discretized system starting from the
representation with the state. In practical cases it is common to calculate AD , BD
by numerical means starting from A, BeTed using the expression of the
series development dieAT .
T2 2 Tn n
eATI AT A +. . .+ A+
2 n!
It is clear that stopping the calculation at an insufficient point
the mistake that is made is small.
A systematic procedure for performing the approximation is the one here.
as indicated below. It is defined
TA TA TA TA TA T2 A2
ψn=I+ 0 I+ ... I+ 1
. . . =I+ + +. . .
2 3 ! n−1 ! n "" 2 3!
T∥A∥
in cuin`e s c e lt o in m o d o c h e s ia s o d d is fa t t a la s e g u e n t e d is e g u an g lianza
<
It is easy to verify that:
TA T2 A2
AD=I+TA(I+ + +. . .)
2 3!
and therefore with good approximation
ADI+TA·ψ
=
∼ n
A
=ẋ+
B
xu
↓
x(k+ 1) = x(k) + TAx(k) + TBu(k)
3.6. Discrete time systems 45
It is also common to assume that the increase in output is proportional to the investment.
secondor, growth factor
P(t+ 1) = 1 + r(1−m)P(t)−rG(t)
* -
A simple model of microeconomics
The price dynamics under equilibrium conditions between demand and supply.
In the case that there is only one good indicating with unit price, condla
demand and supply, it is usual to assume that
d(p) = d0 −ap
s(p) = s0 +bp, a, b > 0
s(t+ 1) = s0+bp(t)
d(t + 1) = d0-ap(t+ 1)
Taking the equilibrium means assuming that s(t + 1) = d(t + 1) and therefore
b d0−s0
p(t + 1) = -p(t) + .
a a
There exists an equilibrium for the price dynamics given by
p(t + 1) = p(t) = d0 -0 .
a+b
The simple model introduced can be used to predict the trend
of the price in conditions of disturbance compared to equilibrium. How easy it is
it will check if there will be a return to the state of equilibrium (stability) if < a,
a removal (instability) seb > a. It is interesting to observe that the two
The coefficients a and b represent, respectively, the propensity to consume.
(of the buyer) and the propensity (ambition) of profit of the entrepreneur.
The economic interpretation of the condition of stability of the equilibrium before
espressa leads to the following conclusion: excessive ambition for profit
the entrepreneur creates market instability.
3.6. Discrete-time systems 47
∂f∂f x−xe
f(x, u) = f(xe , ue ) + )
! partial" derivative
! " +...
)xe ,ue u−u
of ex with respect to u
)
)
∂f ∂f
=f(xe , ue ) + (x−xe ) + (u−ue ) +...
∂x ))xe ,ue ∂u ))xe ,ue
) )
) )
∂h∂h x−xe
h(x, u) = h(xe , ue ) + ! ∂x∂u " ) ! " +...
)xe ,ue u−ue
)
)
and postoz=x−xe ,v=u−u e ,ya =y−h(x e )
∂f ∂f ∂h ∂h
A= B= C= x ,u D=
∂x ))xe ,ue ∂u ))xe ,ue ∂x | e e ∂u ))xe ,ue
) ) )
) ) )
yes, infinite:
ya=Cx(t) + Du(t)
48 3. Representations with the state of finite-dimensional stationary linear systems
A
=
+
Bt)(xutẋ() =x0∈Rn (3.1.a)
0x
y(t) = Cx(t) + Du(t) (3.1.b)
The implicit model introduced immediately suggests a scheme.
physical means through which to realize the system using devices that perform the
50 4. The Linear Time-Invariant Representations in the Time Domain
D
.
u (t ) x(t ) x(t ) +
B + +
! C + and
( )t
A
Figure 4.1
We will now verify the evolutions over time dix(t) ey(t), associated with
a fixed initial state0at the time0and a defined entrance date0from now on,
are described by a model, explicit model, which has the structure
following
t
x(t) = Φ(t−t0 )x(t0 ) + ! H(t−τ)u(τ)dτ(3.2.a)
t0
t
y(t) = Ψ(t−t 0 )x(t0 ) + ! W(t−τ)u(τ)dτ(3.2.b)
t0
where the matrices (Φ, H, Ψ, and W) take on the calculated expressions below.
To this end, let's define the following matrix function,
inf
t Ak t k
I at A2+...= =eAt(3.3)
2 "0 k!
By differentiating the series term by term, it is easy to verify that one obtains
4.1. Continuous-time systems: implicit and explicit representations 51
d At t2
e =A(I+At+A 2+· · ·) =AeAt
dt 2
And it follows immediately that
x(t) = eAtx0
e s o lu z io n e d i
˙ = Ax(t)
x(t) x(0) = x0
This solution is, as is known, unique and globally defined in the linear case.
studio. It is therefore proven that Φ in (3.2.a) takes the expression
Φ(t) = eAt
Place, moreover,
W ( t ) = C Φ ( t ) B + D δ ( t ) = C e AtB+Dδ(t)
t
˙ = AΦ(t−t0 )x(t0 ) +
x(t) AH(t−τ)u(τ)dτ+H(0)u(t)
!
t0
Highlighting on the left between the first two addends and recalling that, in
Based on the definition of H(t), B := H(0), it follows that
˙ = Ax(t) + Bu(t)
x(t)
t
y(t) = Ce A(t−t ) x(t
0
0) + ! (ItAtB+Dδ(t−τ))u(τ)dτ
t0
t t
=C(eA(t−t) x(t00 ) + ! eAtBu(τ)dτ) + ! Dδ(t−τ)u(τ)dτ=Cx(t) +Du(t)
t0 t0
where the property of the unit impulse δ has been taken into account: the integral of
a function multiplied by δ, within the interval of integration, is even
to the value of the function at the moment when the impulse is centered (zero value
of the topic).
It is advisable to note from now on that the reverse process, which in practice-
ICA is presented as the reconstruction of a representation with the state.
(A, B, C, D) starting from a model of forced behavior, it is also,
a significant problem for applications. In fact, as already observed, to a
representation with the state corresponds to the existence of a realization scheme
simulation. The reconstruction of an implicit representation in
starting from a model of forced behavior, it is the problem of the
realization in that, once the differential representation is obtained,
the referenced schema allows for both a physical implementation using the components in-
indicate both the implementation through a computing system for simulation
of the system.
A=T−1 DT
and therefore that
−1
DTt t2 tk
eA t=eT =T−1 T+tT−1 DT+ (T−1 DT) 2+. ..+ (T −1 DT) k+. . .=
2 k!
54 4. The Stationary Linear Representations in the time domain
tk
=T−1 T+tT −1 DT...+T−k D k T k= . . .=T−1 eDt T
k!
α ω
# $t cos(ωt) sin(ωt)
e −α ω =eαt # $
-senωtcosωt
In the presence of real eigenvalues and pairs of conjugate complex numbers, we have
so pearsDtthe expression
4.1. Continuous-time systems: implicit and explicit representations 55
eλ 1t
⎛ .. ⎞
.
⎜ e λ m 1t ⎟
⎜ ⎟
⎜ eα 1t cosine of omega1 t eα tsinω1 t
1
⎟
⎜ ⎟
⎜ 1
−eα tsinω1 t eα t cosω1 t
1
⎟
⎜ .. ⎟
⎜ . ⎟
⎜ ⎟
⎜ ⎟
⎜ eα mt cosine
2 of omegam2t eαmtsinω
2
m 2t ⎟
⎜ ⎟
⎝ ⎠
m m 2
t 2
−eα sinωmt eα t2cosine of omegam2t
and developing the calculations
m1 m2
eAt = eλ it ui v′i+ αt
j
′ ′ bj av ′ -ubv ′ )
"j=1 e cosω
+ j t(uav +ubj v a)j + sinωj j t(u bj aj
"i=1
j j
,
An example
Consider the simple mechanical system composed of mass-spring and
damping addressed in paragraph (2.2). Starting from the dynamic matrix
0 1
A= # k b $
− M − M
√
λ1/2 = -b b2-4kM
2M
and for the corresponding eigenvectors, we obtain
√ √
λ1= -b + b2-4kM -b-b2-4kM
2M 2M →u1= # 2k $
√ √
λ2= −b b2-4kM -b+b2-4kM
2M − 2M →u2= # 2k $
Post
∆=b2-4kM
56 4. The Linear Time-Invariant Representations in the Time Domain
The following three different cases corresponding to parameter values are considered.
for which the value of di∆` is negative, positive, and zero.
If Δ < 0, then it has
√
ua= -b ub= # −∆
# 2k $ 0 $
and with
√
T −1 = # -b-∆ $
2k 0
it is obtained
1 √∆
TAT −1 = -b
√
2M # -∆-b $
and therefore
b cos(ωt) sin(ωt)
eAt=e− t T −12M # $T
-senωtcosωt
−bt
√∆ √∆
=e(cos
2M t(ua v′a+ub vb′) +sen ′ ′
2M 2M t(ua vb−ub va))
If ∆>0, the calculation of the eigenvectors immediately leads to:
√ √∆
T −1 = # -b-Δ-b+ $
2k 2k
1 √∆ 0
TAT −1 = -b+ √∆
2M # 0 -b-
$
tk−1
I i
eΛ tyouΛ t· · · · · · (k−1)! e Λi t
⎛ .. ⎞
.. ..
0 . . .
⎜ .. ··· .. ⎟
eB i t = ⎜ .. .. ⎟
⎜ . ··· . . . ⎟
⎜ .. .. ⎟
⎜ . . youiΛ t ⎟
⎜ ··· ··· ⎟
⎜ 0 0 eΛ it ⎟
⎝ ··· ··· ⎠
oveΛ=λ
i iin the case of real eigenvalue and
αiωi
Λi= # $
−ωiαI
for a pair of conjugate complexes.
For example, consider a four-dimensional matrix with two
pairs of complex eigenvalues coinciding with an even geometric multiplicity
due.
λ1/2=αjω
u1=ua1+jub1(A−α−jω)u1= 0
u2=ua2+jub2(A−α−jω)u2=u1
⇓
α ω 1 0
⎛ −ω α 0 1⎞
TAT −1 = =ΛCj
0 0 α ω
⎜ ⎟
⎝ 0 0 −ω α⎠
cosωt sinωt tcosωt tsinωt
C ⎛ -sin(ωt)cos(ωt) - tsin(ωt)cos(ωt) ⎞
e Λj t
=eαt
0 0 cosωt senωt
⎜ ⎟
⎝ 0 0 -senωtcosωt ⎠
To conclude, let us consider again the previous example in the case where
-b
2M it has multiplicity two. It can
∆ = 0, we are in the case where the eigenvalue is
in this simple case, verify the existence of the following Jordan form
58 4. The Linear Stationary Representations in the Time Domain
1 0
T −1 = #
−-
k
M 1$
k
−- 1
TAT −1 = ⎛ ⎞
M
k
0
⎝ −- M ⎠
The calculation of the exponential can be completed without difficulty.
It is elementary the case in which the matrix Asi reduces to a scalar. We have
A=λ
the free evolution of the state takes place according to the expression
> 0
e t
=0
1
<0
t
Figure 4.2
4.1. Continuous time systems: implicit and explicit representations 59
To characterize the speed at which the system evolves, both growing and decreasing
crescendo, a parameter, τ, the time constant is introduced
1
τ=−
λ
With respect to this parameter
t
x(t) = e− x(0)
τ
and it is clear that it represents, in the case where the state decreases, the time
1
necessary for the state to become of the
e initial value.
In the simplest case, where the size of the system is the matrix A.
it has real and distinct eigenvalues, as seen in the previous paragraph, we have
eλ 1t v′1 n
⎛ .. ⎞ ⎛ .. ⎞ =
eAt=T−1 eΛt T= (u1... un) . . eλ it ui v′i
"
⎠ ⎜⎝ v′n ⎟⎠
λn t
⎝ ee i=1
ci eλ ti
> 0 i 0 = < 0
i x0 x0 x0 i
u ui ui
i
Figure 4.3
1
τi=−
λi
"
in which the time constants play the previously emphasized role.
62 4. The Stationary Linear Representations in the Time Domain
How much is the value in the case of real eigenvalues. Suppose now that there are
also complex eigenvalues, real and complex conjugates, conn=µ+2ν,
µ ν
xl (t) = ci eλ ti ui+ mk eα tksin(ω)
# k t+ϕk )uka+ cos (ωkt+ϕk )ukb $
"i=1 "k=1
having placed
cka
cka=v′ka x0 sinφk=
2 +c2 mk
mk=cka ckb
- kb ckb=v′kb x0 cos φk =
mk
It has been shown that in the case of distinct, real, and paired eigenvalues.
of complex conjugates,
k > 0 k
< 0 k= 0
x0 x0
uk b x0
uk a
Figure 4.4
k0 k= 0 k
< 0
t t t
Figure 4.5
The pulse of the temporal law is ω.k the imaginary coefficient of the eigenvalue
complex.
In conclusion, in the general case of distinct eigenvalues, the free evolution
in the state it turns out to be a linear combination of natural periodic modes,
associated with real eigenvalues, natural pseudoperiodic modes, associated
all pairs of complex eigenvalues. We will say in this case that they have modes.
64 4. The Linear Time-Invariant Representations in the Time Domain
ξ=sinθ= −α
√ α2+ω2
sin =
ωn
Figure 4.6
α=−ξω n e ω=ωn1 - ξ2
.
4.1. Continuous-time systems: implicit and explicit representations 65
In the event that the dynamic matrix is regular, that is, the eigenvalues,
currently not distinct, they have geometric multiplicity, mi , unitary, and possi-
calculate corresponding to each real eigenvalue (pair of eigenvalues
conjugate complexes) as many eigenvectors (eigenspaces of dimension two) as
and the multiplicity of the algebraic eigenvalue of this self-value (of the pair of eigenvalues).
In this situation, the matrix A is still a diagonalizable matrix, as
it can be easily verified using the same arguments as in the previous case. For
natural modes correspond to expressions analogous to the previous ones being able to
coinciding the laws of motion of different natural modes. The natural mode is, in
this case, characterized by the associated self-space.
In conclusion, if the dynamic matrix A has coincident eigenvalues, but
regularly correspond to each eigenvalue with algebraic multiplicity
greater than one, equally simple natural ways, characterized by
from different self-spaces, but with coinciding temporal laws.
We will now examine the case where there are eigenvalues with geometric multiplicity.
greater than one. For this purpose, the procedure is recalled below,
selection of a basis with respect to which a given linear operator takes
the Jordan form.
As is known, such a basis corresponds to a decomposition of [Link]
disjoint subspaces (autospaces) each associated with an eigenvalue and of
dimension equal to its algebraic multiplicity (of dimension equal to double
66 4. The Linear Stationary Representations in the Time Domain
onixtalerm
It is worth remembering that each subspace ofrf tcerdiasU
i tiinad
each associated with a chain of generalized eigenvectors. Furthermore, each
of these last ones, generated by a first-order eigenvector, contains eigenspaces
of dimension one, two,..., up to the maximum order of the generalized eigenvector
because it belongs to that chain.
The two extreme situations that occur in the presence of an eigenvalue
the algebraic multiplicity greater than one is that in which the multiplicity
geometric and unitary, that is, what corresponds to the existence of many eigenspaces of
dimension one in a number equal to the algebraic multiplicity, and it is the case
Already treated previously (Aregular), or the geometric multiplicity
it corresponds to the existence of a single chain
of generalized eigenvectors and the existence of a sequence of eigenspaces.
eralizzati of increasing size up to the algebraic multiplicity. In cases
intermediates have more chains of generalized eigenvectors, each defines
a sequence of autospaces of increasing size. The number of chains,
as is known, it is defined as the index, I, associated with the fixed eigenvalue
or, defines the number of Jordan blocks that make up the most
simple of the operator.
4.2. Discrete-time systems: transition to the explicit model 67
The free evolution is the sum of evolutions associated with it.i Each
It defines a natural multiple way of order, according to the order.
with whichiand the self-regenerating autovettore. The natural mode multiplies.
the expression
tk−1
eAtci=eAt−λ t eλ t ci i=eλi t t(A−λ#i I)ei λ t +. . .+ i
(A−λi I) k−1 eλ tci i $
(k−1)
Note that in this case the temporal law of the mode is of an exponential type
a polynomial coefficient in the time variable.
Each eigenvalue is associated with multiple natural modes, one for each
with autoplace. Therefore, there is an even number of natural ways equal to the index.
The order of the natural mode is assumed to be equal to the dimension of the eigenspace and
it corresponds to the order of the polynomial that characterizes the temporal law.
There will therefore be a set of natural ways in ascending order up to the
algebraic multiplicity. The number of ways, as mentioned, is equal to the index.
(to the number of Jordan blocks) and the size is that of the corresponding
Jordan block. In this case, we will say that we are in the presence of nat-
multiple distinct urals.
In the two limiting cases of geometric multiplicity equal to algebraic multiplicity
or equal to one we will say that there is only one natural multiple more
number equal to the algebraic multiplicity) natural multiple ways with laws
coinciding temporals.
We will now verify that the evolutions over time dix(t) ey(t) associated with
a fixed initial state0 at the time0 and a defined entrance date0in
then, they are described by the following explicit model
Φ(t) = At
W ( t ) = C Φ ( t − 1 ) B = C A (t−1)B t >0W(0) =D
Even in this case, the observations made for the systems can be repeated.
at a continuous tempo:
the responses in the state and outgoing result are the sum of two
evolutions: free evolution and forced evolution; forced evolution has
the form of a convolution summation;
HeW, impulse response matrices in the state and output, respect-
again, they have columns that can be interpreted as forced responses
at impulsive entrances;
For forced responses, the principle of superposition of effects applies;
The forced behavior of entry and exit is described by W, response
impulsive (matrix of impulsive responses), a model of behavior
forcibly
In the terminology of Systems Theory, the model (1.a) - (1.b) is
implicit or differential said, the corresponding model (2.a) - (2.b) is said
explicit;
In the calculation of the explicit model, power plays a central role.
of the dynamic matrix.
4.2. Discrete-time systems: transition to the explicit model 69
A=T−1 DT ⇒A=
t (T−1 DT) t. . .=T−1 D t T.
Being block diagonal, the power is equal to the diagonal of the powers.
of the individual blocks. There are no issues with the power of the diagonal part
resulting equal to the diagonal of the powers of the elements on the diagonal.
It is only necessary to understand how to calculate the power of the block (2×2)
associated with a pair of complex conjugate eigenvalues, given
λ=α+jω=σe jθ =σ(cosθ+jsenθ)
70 4. The Stationary Linear Representations in the Time Domain
⇒α=σcosθ ω=σsinθ
it results
α ω cosθsenθ
# $ =σ # $
−ω α -sin(θ)cos(θ)
it's not difficult to calculate
t
cosθsenθ cosθt sinθt
# $ = # $
-sinθcosθ -sinθcosθ
So for the generic block (2×2) associated with a complex pair
He/she is tired.
cosθt sinθt
Λt=σt # $
−sin(θ)cos(θ)
It therefore results
λt1
⎛ .. ⎞
.
⎜ λtm1 ⎟
⎜ ⎟
⎜
t ⎜ σt1cosθ1 tσt1sinθ1 t ⎟
D= ⎜ ⎟
⎟
⎜ −σt1sinθ1 tσt1cosθ1 t ⎟
⎜ .. ⎟
⎜ . ⎟
⎜ ⎟
⎜ σtmcosθ t
m2tσmsinθm2t
⎟
⎜ ⎟
⎜ ⎟
2 2
⎝ −σtmsinθ t 2
2 mtσcosθ
m m 2t 2 ⎠
and developing the calculations
m1 m2
At= λitui v′i + t ′ ′ bj av ′ −ubv ′ )
"j=1 σcosθ
j + j t(uav +ubvaj) + sinθj t(u bj aj
"i=1
j j j j
,
The student will not have difficulty revisiting, with reference to the calculation of
power of a matrix, the considerations made earlier in the case of
eigenvalues with geometric multiplicity greater than one.
For the general case, the same considerations already made for the calculation apply.
color of the exponential that in the specific case gives for the generic block, of
dimensions,
4.2. Discrete-time systems: transition to the explicit model 71
t[k−1]
ΛItt[1]Λi t−1
· · · · · · (k−1)!Λi
t−k+1
⎛ .. ⎞
.. ..
0 . . ··· .
⎜ .. .. ⎟
B=i t ⎜ .. .. ⎟
⎜ . ··· . . . ⎟
⎜ .. .. ⎟
⎜ . . t Λi t−1
[1] ⎟
⎜ ··· ··· ⎟
⎜ 0 ΛI t ⎟
⎝ ··· ··· 0 ⎠
oveΛ=λ
i iin the case of real eigenvalue and
cosθiscentI
ΛI=σi # $
-sinθicosθI
for a pair of conjugate complexes and, furthermore, t[k] denotes the polynomial
factorial of a defined order as
t[k] t(t−1)...(t−k+1)
love
λℓ =|λ ℓ |ejθ
ℓ
, 0<thetaℓ <π
Im [ ]
1Re [ ]
Figure 4.7
In the case of eigenvalues with geometric multiplicity greater than one, with
Similar considerations to those made in the continuous time case are calculated:
t
t h t−h
" # h$ (A−λi I) ci λi
h=0
4.3. Further remarks * 73
t t(t−1)...(t−h+1) t(h)
(# $ ) = :=
h h! h!
Therefore, the laws of motion are of the type power with coefficients that are polynomials.
factorials.
D λCI T
λ=e
i
n
C
AD=eAT = eλ iTui v′I
"i=1
in fact
C
eλi T
→λiD
The concept of natural mode now introduced lends itself to carrying out a first
investigation into what was mentioned earlier regarding the representativeness of
model consisting only of forced response (impulse response matrix).
4.4. Natural modes in the state and outgoing 75
Let us first consider the case in which the natural ways are distinct; that is,
which corresponds to the coincidence between the algebraic and geometric multiplicities of
an eigenvalue. In this case, in fact, the temporal laws associated with the modes are
distinct.
It is evident that while changing dix0in the free evolution of
the system will be able to show all the temporal laws of natural ways, such as
it is said that free evolution is characterized by all natural ways, for
As for the forced evolution in the state, that is no longer true. In fact,
the forced response is characterized by all and only the 'temporal' laws that
they appear ineAtB.
Such modes are said to be excitability (and it is implied 'with input impulses')
as an impulsive input is applied, it turns out
t
x(t) = ! eA(t−τ ) Bδ(τ)dτ=e AtB
0
that is, the response coincides with the free evolution that moves from the initial state
x0=B, in the case of single-input systems, from an initial state equal
to a fixed linear combination of the columns of B, in the case of systems with
more entries. The mode associated with the eigenvalue λi , and therefore excitable if
influenza the evolution in the state. More precisely, we will say it is excitable if the
your temporal law appears in H(t).
The excitability condition of the generic way naturally in the case of
distinct eigenvalues are
v′iB ≠ 0
In fact, from the expression
n
eAtB= eλ it ui v′iB
"i=1
76 4. The Stationary Linear Representations in the Time Domain
Cu inot equal to 0
what is evident from the expression
n
This At= eλ it Cu i v′i
"i=1
One can only wonder what are the temporal laws that characterize
W(t); and therefore the forced response at the output. It is clear that only the natural modes,
and therefore the relative temporal laws, which are simultaneously excitable and
Observables may appear in the forced response.
It obviously concerns all and only the natural ways for which it results
v′iB̸= 0 Cu i≠ 0
The previous considerations actually apply in the case where the matrix of-
she is regular with the caution of distinguishing between natural ways that have
the same temporal law. In this case, it will not be the presence (inΦ, H,Ψ,
of the law of motion corresponding to the natural way to certify it
property
1 1
A= # -1 $ B= # $ C= ( 0 1 )
0 -1 0
There are two natural ways, one of first order and one of second order; that one
first order and temporary law− it can be seen, it is not observable, while
the second order is not excitable, but observable. Therefore, there is no...
simultaneously excitable and observable, which is reflected in the fact that
the impulsive response is identically zero.
e−tte−t 1
(01) # 0 =0
e−t $ # $0
4.5. Appendix 3
(A−λI)v = 0
The previous equation has a solution if λ is zero from the characteristic polynomial.
isto, d(λ), associated with the matrix A
µr
d(λ) = det(λI - A) = (s - λ1 )µ . . .(s−λ1 r )
mr
m(λ) = (λ−λ1 )m. . .(λ−λ
1
r)
80 4. The Stationary Linear Representations in the time domain
Almost diagonal form with real coefficients, D, (diagonal for the eigenvalues
real and block diagonal of size two for the pairs of co-complexes
A transformation to put this into canonical form exists.
the matrix is said to be regular if all the eigenvalues
they have unitary geometric multiplicity. How could one show the con-
the definition of unitary geometric multiplicity implies that in correspondence to
each eigenvalue λiof algebraic multiplicityi It is possible to calculate proprioµI
first order auto-vectors, that is, independent vectors solutions of
(A−λi I)u I= 0
The transformation we are looking for is the one that corresponds to as-
to take as a new basis of the space that associated with the eigenvectors of the first
[Link] that, as can be demonstrated, are independent of each other.
v′1
⎛ .. ⎞
.
⎜ vm′ 1 ⎟⎟
⎜
D=TAT−1 = ⎜⎜ v′1a ⎟⎟ ·A·(u 1... umu1au1bone
. . . u m a um 2 2 )
⎜ v′1b ⎟⎟
b
⎜
⎜ .. ′ ⎟
⎜ .vm2 a ⎟⎟
⎜
⎜ ′
vm ⎟
⎝ 2b ⎠
λ1
⎛ ..
. ⎞
⎜ λ m1 ⎟
⎜ ⎟
⎜ α1ω1 ⎟
= ⎜ ⎟
⎜ −ω1α1 ⎟
⎜ ⎟
⎜ .. ⎟
⎜ . ⎟
⎜ ⎟
⎜ α m2 ωm2 ⎟
⎜ ⎟
⎝ −ωmαm 2 2 ⎠
where do you go′I Ii′the rows of the matrix and the columns of T are indicated−1 ,
respectively.
The previous representation constitutes the almost diagonal form at co-
real eigenvalues of a matrix. The verification is soon done starting from the following
equality obtained from the previous by multiplying both members by
T −1
A(u1. . . u mu1au1b1 . . . u m a um b ) = 2 2
= (u1. . . u mu1au1b1 . . . u m a um b )· 2 2
λ1
⎛ ..
. ⎞
⎜ λ m1 ⎟
⎜ ⎟
⎜ α1ω1 ⎟
· ⎜⎜ ⎟
⎟
⎜ −ω1α1 ⎟
⎜ .. ⎟
⎜ . ⎟
⎜ ⎟
⎜ α m2 ωm2 ⎟
⎜ ⎟
⎝ −ωmα2m 2 ⎠
82 4. The Linear Time-Invariant Representations in the Time Domain
The equality among the primaries1columns of products in the first and second member
the assured fax to c h e λiIithey are eigenvalue and corresponding eigenvector
e, therefore, (A−λi I)u i= 0 that is equivalent to writing Aui=λi ui(beingAuIthe
i-ma column of the product to the first member eλi uilaima column of the product
a second member). The equality between the remaining columns is ensured by
property on the real and imaginary part to which the eigenvectors satisfy
complex. If, in fact, λ = (α + jω) is the eigenvalue associated with the eigenvector
u=ua+jub , yes it has
(A - (α + jω)I)(ua+jub ) = 0
what is equivalent to the two equalities between real and imaginary part
(A−αI)ua=−ωub (A−αI)ub=ωua
In the general case of a non-regular matrix, the presence of eigenvalues with multiplicity
geometric greater than one, the simplest canonical form it takes on
the operator is the so-called Jordan blocks, J.
In that case, by appropriately choosing the new coordinates in the au-
subspaces associated with eigenvalues, eigenspaces that can be further decomposed into sub-
associated spaces to "chains of generalized eigenvectors," a form is obtained
of the type:
J1 0
J= ⎛ .. ⎞
.
⎝ 0 Jm ⎠
4.5. Appendix 3 83
ΛI I
⎛ .. .. ⎞
. . αiωi
Bi= ⎜ .. ⎟ , Λi= # $
⎜ . I ⎟ −ωiαi
⎜ .. ⎟
⎜ . ΛI ⎟⎠
⎝
for a pair of complex conjugate eigenvalues. The maximum dimension
the blocksBiit coincides with the maximum order of the generalized eigenvectors
associated with the eigenvalue λi .
It is remembered that a generalized eigenvector of order k satisfies the equations
j (j−1)
(A−λi I) i u=u
i i
which is defined for Re[s] > α,f being αf a real number associated
adf(t) is called the convergence abscissa. This operation is linear and holds true.
important properties that are the basis of the calculation of transformations of classes of
functions. Among these, the following, known as the comet theorem of derivation, connects
the transform of the derivative of a function to its transform. Deriving
both members of the previous identity, the derivative's transform takes on
the following expression
d
L " dt f ( t ) #= s F ( s ) - f ( 0 ) .
These simple calls allow for the calculation of the Laplace transform of
both members of the equations that describe a linear representation,
stationary, finite size, regular. It is obtained, indicating with the letters
capital letters the transforms of the corresponding time functions,
t t
x(t) = eAtx0+ ! eA(t−τ ) Bu(τ)dτ=φ(t)x0 + ! H(t−τ)u(τ)dτ
0 0
t t
y(t) = Ce Atx0+ ! This A(t−τ )Bu(τ)dτ + Du(t) = ψ(t)x0+ ! W(t−τ)u(τ)dτ
0 0
88 5. The stationary linear representations in the complex domain
e
t
L "! w(t−τ)u(τ)dτ=W(s)U(s)
#
0
What has been presented so far provides an alternative method for calculating the
answers or, what is equivalent, for the transition to the representation ex-
explicit. Once an input is assigned, and its Laplace transform is calculated, the
calculation of the evolution at output starting from an initial state0behave
the calculation of (sI−A) −1 , simple multiplications and an operation of anti-
transformation. The antitransformation operation is particularly sem-
splice in the presence of inputs that admit rational transforms; in that case
the free and enforced response are themselves rational functions that can
to be anti-transformed without difficulty once the corresponding ones are calculated
expansions in partial fractions.
Consider, as an example, a system that has a transfer function.
1
W(s) =
s+ 1
and calculate the forced response to the unit step input, u(t) = δ−1 (t).
1
After all, how easy it is to calculate L[δ1 (t)] = , if
s has
-1 1 1 1
yf (t) = L [ ] =L−1 [− ] = (1−e−t )δ−1 (t)
s(s + 1) s s+ 1
It is observed that in the calculation it is good practice in the operation of anti-transfer
definition, to emphasize that the function obtained is defined from zero in
then multiply on the right by the function δ−1 (t).
The following delves into the aspects of computation in the domain com-
the junction and the connections are established with the analysis conducted in the chapter
precedent.
Let it be assumed that the matrix A, (n×n), has eigenvalues (λ1 , . . . ,λ r ) with
algebraic multiplicity (µ1 , . . . , µ r ) (µ1+. . .+µ r=n). Since the calculation of the
transition matrix in the complex variable reduces to the calculation of the inverse
by (sI−A), one obtains
T
(sI−A)a E(s) E(s)
(sI−A)−1 = = &1 ' mr
= & '
|sI−A| (s−λ1 )m. . .(s−λr) m(s)
90 5. The stationary linear representations in the complex domain
mr
m(s) = (s−λ1 )m ...(s−λ
1
r)
least common multiple of the polynomials in the denominator of the elements of the
matrix (sI−A) −1 Thsi poylnomail si ceald hte mnm
i al [Link] tI
it deals with a polynomial that is a factor of the characteristic polynomial
µr
d(s) = (s−λ1 )µ ...(s−λ
1
r)
since it has the same roots as the latter, but each with a
multiplicityi geometric multiplicity, which is less than or equal to
the algebraic (µi An indirect demonstration of the fact that all self-
values appear between the zeros of the minimum polynomial is observed that if so
if it weren't for a fixed eigenvalue then the corresponding exponential law
would not appear in the matrix's exponential; what we know not to be
possible.
In the particular case of unitary geometric multiplicity, there is
E(s)
(sI−A)−1 = & '
(s−λ1 . . .(s−λr )
µ ν
RI Rk R∗
k
= + ) + * µ + 2ν = r
(i=1 s−λi (k=1 s−λk s−λ∗k
Rk=Rka+jRkb λk=αk+jωk
it is rewritten
µ ν
RI 2Rka (s−αk )−2Rkb ωk
(sI−A)−1 = + ) *
(i=1 s−λI (
k=1
(s−αk )2+ω2k
5.2. The transition matrix in the complex domain 91
Precision in the following, the laws of motion of natural modes are more complex.
in fact, the expansion in partial fractions is composed not only of addends with
first degree polynomials in the denominator (simple fractions), but also from
fractions with a denominator of degree greater than one (compound fractions) that are
the transforms of exponential functions to polynomial coefficients int. Is this
the case that arises when the matrix is regular, that is, it does not admit
the diagonal form, but that of Jordan. In this case the size of the
The Jordan block of the largest size coincides with the geometric multiplicity.
To verify with a simple example how it can happen that in the expression
the polynomial in the denominator does not coincide with the characteristic one and
if degree < n, consider
0 1 s+1 0
A= + -1 , (sI−A)−1 = + =
0 -1 (s + 1) 2 0 s+ 1 ,
1
(s+1) 0 1
= + 1 , = (s+ 1) I
0 (s+1)
It is a matrix of size two with the only eigenvalue -1 that pre-
algebraic multiplicity two and geometric multiplicity one. The geometric multiplicity
unitary accounts for the presence of a simple exponential law in the
exponential of this matrix.
This does not happen in this other simple case.
1 1 s+1 1
A= + -1 , (sI−A)−1 = =
0 -1 (s+ 1)2 + 0 s+ 1 ,
1 1
(s+1) (s+1)2 1 (s + 1) 1
= + , = (s+ 1)2 +
0 1
(s+1)
0 (s + 1) ,
in which the geometric multiplicity of the eigenvalue is equal to the algebraic one and
ofm
egnechertsesxim
rhetat ofphispornshet ,esacshitnI [Link]
una legge esponenziale non semplice nell’elemento fuori dalla diagonale.
92 5. The stationary linear representations in the complex domain
λ1·µ1
E(s)
A mi= 1⇔(sI−A)−1 = & '
- (s−λ1 ...(s−λr )
-
λr·µr
r
R1 R2 Rr Ri
= + +. . .+ = =L(eAt)
s−λ1 s−λ2 s−λr (i=1 s−λI
If it has indeed
r
L−1(sI−A)−1 = Ri eλ ti
& ' (i=1
µi
Ri= uij v′ij=Ui Vi′
(j=1
doveij(j= 1, . . . , µI they are the right eigenvectors, all of the first order in
this case, associated with the eigenvalue λi . The matrixRi (n×n), has rank µI
being the vectorsijindependent.
A basic example:
1
A= + -1 , (sI−A) = +
s+ 1 -1
0 -2 0 s+2 ,
1 s+2 1 R1 R2
(sI−A)−1 = = +
(s + 1)(s + 2) + 0 s+ 1 , s+ 1 s + 2
conR1eR2residue matrices (2×2) of unit rank (distinct eigenvalues)
obtained by expansion into simple fractions
1 1 0-1
R1 R2 + 0 0, +0 1 ,
φ(s) = + = +
s+ 1 s + 2 s+ 1 s+2
1 1 0-1
φ(t) = + e−t + + e−2t
0 0, 0 1 ,
µ ν
RI 2Rka (s−αk )−2Rkb ωk
(sI−A)−1 = +
(i=1 s−λI (k=1 (s−αk )2+ω2k
Recalling the expression of the transition matrix int, it is not difficult to verify
make sure that
′ ′
vka vkb
Rk=Rka+jRkb= (uka+jukb ) -j 2
&2 '
and, developing the calculation, we have
1 1
Rka=(uka v′ka+ukb v′kb ) Rkb(ukb v2′ka-uka v′kb ).
2
Multiplying by the initial state:
Ri x0=ci ui
= mk(s−α
& k
)senϕk +ωk cos phikuka+mk(s−αk )cosϕk −ωk senϕkukb
' & '
(s−αk )2+ω2k
One has taken into account the definitions of the parameters. In fact, by antitransforming
according to Laplace, aperiodic and pseudoperiodic modes are obtained
µ ν 2Rka s 2Rkb
RI τi ωnkωnk( +ζk )− ωnk 1−ζk2
(sI−A)−1 = + 2zk s2
.
(i=1 1 + τi s (k=1 1+ s+
ωnk ω2nk
5.4. The transfer function 95
1 dm−k
i
−1 m
RI= lim + (sI−A) (s−λI ) i ,
s→λi i
(mi-k)!dsm−k
By identity with the calculated expressions, one obtains the link between the
residualIand the vectorsicomponents ten0in the autospacesi ) It is obtained
ci=Ri1 x0 (A−λi I)c i=RI2 x0 (A−λi I) m−1 cim i−1 =RI amxI0
Based on what has been observed, the transfer function can be calculated
side, in a so-called black box modeling approach in which it is assumed
of not having other information about the system other than the input measurements and
corresponding outputs, starting from the knowledge of the transform of a
any forced input-output pair. In fact, for a system with a single input
and an exit one has
Yf(s)
W(s) =
U(s)
and if the input is the unit impulse δ(t), since L(δ(t)) = 1, the transform
the forced exit is precisely the transfer function.
For systems with more inputs and more outputs, the characterization of components
the matrix of transfer functions allows to define the experts-
minds that it is necessary to conduct in order to calculate the W. From the expression of the
matrix of transfer functions
w11 (s)... w1p (s)
W(s) =⎛ ... ... ... ⎞
⎝ wq1 (s). . . wqp s ⎠
5.4. The transfer function 97
it is deduced
yi
yi=wi1 u1+. . .+wipup⇒wij= uj --uk =0k ≠ j
-
-
In reality, the identification procedure outlined now consists of the calculation-
create the transfer function from input-output experiments,
corresponds, in the time domain, to calculating the average impulse response
before deconvolution operations. Such an identification procedure and
not easy to apply in practice due to the intrinsic complexity of the measurements in
dynamic regime and the related elaborations, particularly in pre-
without errors and measurement uncertainties. We will see in the following of the chapter how
the values that the transfer function of a linear system takes
They can be directly connected to a part of the answer in
time, the steady-state response. This allows us to enable
a simple and robust experimental procedure regarding uncertainties
of measure for the identification of the transfer function.
Before delving into the details of these aspects that are at the core of a
characterization of the system behavior in relation to its behavior
In frequency, it is important to observe that the transfer function model
This can also be obtained directly through the so-called approach
a transparent box in which it is assumed to have a model of the compo-
focus on the system and its connections are known. Procedures of this type,
based on graph manipulation, are available for modeling of
mechanical and electrical systems. Among these, the most well-known in engineering is the
method of complex impedances characteristic of circuit theory. With
reference to this method is indeed known that associated with each component
a complex impedance is constructed the graph, which expresses the connections of
circuit, appropriate manipulations define a systematic procedure of
calculation of the complex impedance of the circuit, which is nothing more than the function
of transfer.
Since the model in the domain of the complex variable can, on par with the
representation with the state, being the starting point in the study of a
assigned physical system and it is necessary to specify its representativeness.
A simple example that illustrates to what extent the function
98 5. The stationary linear representations in the complex domain
the poles of the transfer function (at the denominator in the general case)
the poles of the transfer function matrix). In light of what
Observed, the poles are a subset of the eigenvalues of the system.
The zeros of the polynomial in the numerator in the case p=q=1 (the common zeros
determinants of polynomial matrices, of size equal to the minimum trap
equivalent, which are obtained by considering all the combinations of columns among the ep
available seq < p, which are obtained by considering all the combinations of dip
lines between the possible sep < q), are called zeros of the transfer function.
matrix of transfer function transmission zeros.
The meaning of the poles is known. The poles are the eigenvalues that appear in the
transfer function: the corresponding laws of motion, aperiodic or pseudo-
periodic, convergent, constant or divergent, characterize the behavior of
forced amendment.
The zeros of the transfer function also allow for an interpretation.
in terms of behavior input output. Siazuno zero (z1ez2a couple
of complex conjugate zeros) and siau(t) an input that admits transformed of
Laplace
1 1
( )
(s−z) (s−z1 ) (s−z2 )
Now calculate the forced response and observe that, following the cancellation-
the relationship between the transformed input and zero (the pair of zeros)
the trend over time of such a response does not contain the terms that correspond
both in the development in partial fractions at the poles that belong to the transform
of the entrance. In other words, that type of entrance generates a forced response.
which boils down to a free evolution. As we will have the opportunity to better c
subsequently taking the presence of one or more zeros translates into the ability of the
system for filtering the effect of a class of input functions. This characteristic
characteristic, which can be thought of as a lack of reactivity of the system
with respect to certain stresses, it is even more evident (and in some ways
surprising) when one thinks of such a filtering action due to the
presence of zeros in the positive real part. In this case, in fact, the filter action
occurs in the presence of an entrance whose width tends to infinity as it increases
of time.
100 5. The stationary linear representations in the complex domain
It will return to these aspects in the next chapter where, with reference
AI systems with multiple inputs and outputs, the study will be resumed based on the
following definition of zero based on state representation. The
The transmission frequencies of the system coincide with the values that cancel the
determinant of the so-called system matrix. These are the values
such
s∗ I−A−B
det ) *=0
C 0
As we will have the opportunity to show, such a definition is equivalent to the prayer.
tooth.
Ratio of polynomials
This denomination is taken from the representation already introduced in the para-
graph 5.1.1
B0 +B1 s+. . .+Bm sm
W(s) =
a0 +a1 s+. . .+sN
oveBi , (q, p), are appropriate coefficient matrices.
As an example, consider
s(s + 2) (s2-1)
1 s−1 + s(s + 1) 0 ,
s+1 s2 +2s
W(s) =+ 1 , =
s+2 0 s(s + 1)(s + 2)
0-1 2 0 1 1 2
+ 0 0 , + + 1 0 , s+ + 1 0,
s
W(s) =
s3+ 3s2+ 2s
Bode's form
Let us initially consider the following factorization.
real coefficients of a polynomial:
=h(s−λ
5I i ) (s2-2α 2 2
5 k s+αk+ωk)
k
=h(−λ
5II ) (1 5 +τi s) (s2-2α5k s+α2k+ω2k)
i k
2zk s2
p(s) = h(-λi ) (ω2nk ) (1 +τi s) ) 1+ s+ 2 *
5i 5k 5i 5k ωnk ωnk
αk2+ω2k=ω2nk
ωnk ζk=−αk
Starting from the previous expression, one obtains
2zk′ s2
h′ (1 + τi′ s) ) 1+ ′
ωnk s+ ′2
ωnk *
4i 4k
W(s) =
2zk s2
h (1 + τi s) ) 1+ ωnk s+ ω2nk
*
4i 4k
ν′ µ′
2z k′ s2 ′
(1+τi′ s) ) 1+ ′
ω nk s+ 2
ω nk
*
4
i=1 4
k=1
W(s) = k
ν µ
2zk s2
sr (1+τi s) ) 1+ ωnk s+ ω 2nk *
4
i=1 4
k=1
The study in the complex domain allows us to understand how the characteristics
characteristics of the response of a dynamic system are related to the behavior
treatment regarding a particular class of solicitations: the solicitations
periodic. These aspects are the basis of the so-called frequency analysis
of a dynamic system; an approach that engineers use in
study of dynamic systems. The starting point is represented by the charact-
characterization of the physical meaning of the transfer function evaluated on
points of the imaginary axis and rests on the concept of steady-state response
none.
The steady-state response to a given input is defined as
that time function at which, regardless of the initial state,
the response tends to increase over time.
If you remember that the output response is composed of free evolution and
forced, it is understood that in order to achieve independence from the initial state
it is necessary that you read the code that accompanies inCeAttend to zero at
growing of time. In other words, the condition of independence from the state
Initially, it is reduced to requiring that the observable modes are associated with eigenvalues.
a negative real part. In fact, under this condition, in
t
y(t) = Ce A(t−t) x00+ ! W(t−τ)u(τ)dτ
t0
5.5. The steady-state response 103
the generic free evolution for a fixed point0tends to zero as the increases
time. In reality, for a physical system, the existence of the steady state does not
can do without an additional property that makes correct function possible
system zoning: it concerns the limitation of all evolutions
To understand this, as can be grasped in light of the...
God led in natural ways and as we will have the opportunity to specify in the
successive phases of our study, having already imposed that the eigenvalues asso-
If observable modes have a negative real part, it is necessary that the others
eigenvalues have strictly negative real part if with geometric multiplicity
metric greater than one, real part less than or equal to zero if with multiplicity
unitary geometry. In other words, that property which is known is necessary.
as internal stability of the system. We will therefore assume internal stability
In addition, to ensure the aforementioned independence from the initial state as-
we will assume that the eigenvalues associated with the observable modes have a real part
strictly negative regardless of their multiplicity.
Limited in this way to the class of systems to which we will refer, those
which have eigenvalues with negative real part and possibly zero those do not
observable with unit geometric multiplicity, we will define the steady-state response
permanent like that function obtained by taking the limit as it approaches.0how'sitgoing
to the response of the system.
More precisely, we define
t
yr (t) = limit ! W(t−τ)u(τ)dτ
t0 →−∞ t0
-u(τ)dτ- ! W ( ξ ) u ( T a+t−ξ)dξ
0
from which follows
∞
t
yr (t) = lim ! W(t−τ)eȷ ωτ dτ
t0 →−∞ t0
and put (t−τ) = ξ, developing the calculation remembering that due to the regularity of the
the functions involved the limit shifts to the boundary of integration, it is obtained
∞
yr (t) = eȷ ωt
! W ( ξ ) e −ȷ ωξ dξ=eȷ ωt W ( s )
0 -s=ȷ ω
-
The last equality is guaranteed by the fact that being the abscissa of convergence
the Laplace transform is negative being W(t) a linear combination
of exponentials with negative real part. By exploiting linearity and observing that
W(ȷ ω) in the quality of the transform of a real function verifies
W ( - j ω ) = M ( ω ) e −ȷ φ(ω)
106 5. The stationary linear representations in the complex domain
ove
M(ω) = |W(ȷ ω)|, φ(ω) = W(ȷ
̸ ω)
it is obtained
It will proceed now, assuming that the system meets the aforementioned condition of
existence, to the calculation of the steady-state response to so-called inputs
canonical; that is of the type
tk
u(t) = δ−1 (t)
k!
108 5. The stationary linear representations in the complex domain
t
rk
yr (t) = lim ! W(t−r)dr= t−r=ξ
t0 →−∞ t0 k!
∞
(t−ξ)k
= ! W(ξ) dξ
0 k!
∞ k
t k−i i i k
= ! W(ξ)· ξ -1 + i , dξ
0 (i=0 k!
k ∞
1 k k−i i
= t -1 W ( ξ ) ξ i dξ
(i=0 k! + i , !
0
k
1
= -1i tk−i Mi
(i=0 i!(k−i)!
Sia therefore obtains a polynomial of degree k whose coefficients are the parameters.
Mi
∞
M i= ! W ( ξ ) ξ i dξ
0
determinants of the impulsive response. It should be noted that such coefficients are
I define it in terms of the assumptions about the system in the impulse response com-
They seem to be temporal laws of modes associated with eigenvalues with negative real parts;
what ensures the boundedness of the integral.
The moments of the impulsive response are connected by a simple re-
Relation to the coefficients of the Maclaurin series expansion of W(s). It follows that
in fact:
i
dW(s)
Mi-1I
dsi --s=0
-
Recalling the development of Mc Laurin W(s) = - c si , in which
i
6
i
1d(W(s)
c i= -
i!dsI -s=0
-
-
5.5. The steady-state response 109
it is obtained
-1i Mi=i!ci
⇓
MI-1i i!cI
=L(W−1⋆u′ )
110 5. The stationary linear representations in the complex domain
where the prime indicates the derivative with respect to time and
W(s)
W−1 (s) =
s
So the forced response, equal to the convolution of the impulse response with W.
siva), pu`o also be obtained as the convolution of W−1(forced response)
to the input at unit step) with the derivative of the input if, as occurs
In the case of the inputs being considered here, the derivative at zero is zero.
These topics can be generalized to higher order derivatives.
or the inputs and to the responses to canonical inputs of higher order if so
Note that the Laplace transform of the canonical input of order k:
tk
u(t) = δ−1 (t)
k!
vale
tk 1
L + k! δ−1 (t) =, sk+1
In fact
1 1
Lδ−1 (t) = Ltδ−1 (t) =
& ' s & ' s2
in fact−1 (t)`e the derivative of ditδ−1 (t) and therefore
1 1
Lδ−1 (t) = s· s2 −0 = s
& '
Generalizing
t (k−1) 1 tk 1
L (k−1)! δ−1 (t) = sk ⇒L &k!
δ−1 (t) =
s( k+ 1)
& ' '
Returning to the study of the steady-state response, consider the response
forced at the entrance of ordinek,
From this, it is evident that, having the polypsinegative real part, as time increases
the response tends to take on an even trend
-1 A0+. . .+Ak sk
yr (t) = L + ,
sk+1
W ( s ) k+1
A0= s - W(0) = k gain
sk+1 -s=0
-
-
d W ( s ) k+1 dW
A1= s s=0 =
+
ds s k+1 , | ds --s=0
-
-
|
1dk
Ak= W ( s )-
kids k -s=0
-
-
↓
tk tk−1
yr (t) = A0+A1 +. . .+Ak
k! (k−1)!
Ultimately, the previous expression is obtained where with AiI am here
indicate the coefficientsiof the series expansion of Mc Laurin. To the calculation of such
the coefficients reduce the calculation of the steady-state response.
W ( s ) A0 A1 An
y(s) = = + +. . .+
s s s−p1 s−pn
where doveiThe poles of the transfer function have been indicated. Again a
turn, sincei )<0 it is understood that the anti-transformed of the response
forced is a function that tends to the constant value A as time increases0
=k, profit. In fact:
n
yf (t) = A0 δ−1 (t) + A1 eptδi −1 (t)
(i=1
limf(t) = limsF(s)
t→∞ s→0
it is obtained immediately
moreover
y ′ (0) = lims sy(s)−y(0) = limsW(s) =
s→∞ & ' s approaches infinity
=bn−1let m = n - 1 = 0 if m < n−1
5.5. The steady-state response 113
..
.
y (n−m) (0) = bm
ẋ=x
1 2
7 k b u(t)
̇ -x1−
2x
= x2+
M M M
y=x1
-1
s -1 0
W(s) = ( 1 0 +
) k
M s+ Mb , +
M
1 ,
b
1 s+ m 1 0
= (10) +
s s + b
M + k
M −
k
M s ,+ 1
M
,
& '
1
M
⇒W(s) = s2+ Ms+ b k
M
The gain, the steady-state displacement under the action of a unit force, is
to Paris
1
A0=
k
while the behavior around zero is characterized by
1
y(0) = 0 y ′ (0) = 0 y ′′(0) =
M
climbing time, ts time needed to reach for the first time the
steady state value (in the oscillating case); time to go from 0.1 to 0.9 of
final value in the non-oscillating case.
λ1m1≤µ1
.. .. ..
. . .
λrmr≤µr
The calculation of motion laws in the general case involves the anti-transformation.
the second member in which RIthey have expressions
1 dm −k
i
−1 m
Rik= lim + (sI−A) (s−λI )
I
,
s→λi I
(mI-k)!dsm−k
5.7. Discrete-time systems: analysis in the complex variable domain 115
It is obtained
r mi k
1
φ(t) = Rik L−1 ) *
(i=1i=1
( s−λi
Because
tk−1 1
L) eλ it * =
(k−1)! (s−λI )k
yes, ultimately
r mi
tk−1
φ(t) = RI eλ it
(i=1i=1
( (k−1)!
The temporary laws are thus of the type
which is defined for |z|>ρ f,ρf radius of convergence associated with the fun-
zionef.
116 5. The stationary linear representations in the complex domain
A Zt = (zI−A)−1 z
$ %
it is also implicitly demonstrated that the Z-transform of a con-
evolution is equal to the product of the transforms
discreet
For stationary discrete-time linear representations, it is possible to
develop similar considerations to those made using the Z-transform.
Assigned a function f(t) defined in Z+ ∪0 is defined Z transformed
dif(t) the function of complex variable
∞
Z[f(t)] = F(z) = f ( t ) z −t
(t=0
which is defined for |z|>ρ f,ρfradius of convergence associated with the function
zionef.
The transformation induced by such an operation is linear and holds true in-between
essential properties that are the basis for the calculation of transforms of classes of
5.7. Discrete-time systems: analysis in the domain of complex variables 117
These simple calls allow for the calculation of the zeta transform of
both the members of the equations that describe a linear representation,
stationary, at a finite size. It is obtained, indicating with capital letters
the transforms of the corresponding functions
x(t) = At x0+Σt−1
0 A
(t−τ −1) Bu(τ)
y(t) = A t x0+Σt−1
0 CA
(t−τ −1) Bu(τ) + Du(t)
ZAt(zI−A)−1 z
$ %
118 5. The stationary linear representations in the complex domain
e
ZΣ0t w(t−τ)u(τ) := (A(t−1) ∗Bu(t)) = W(z)U(z)
$ %
These identities express two properties of the Laplace transform.
which explain the interest in the use in the study of linear dynamic systems.
transformed power functions, including constant functions
periodic and quasi-periodic) are proper rational functions; the transform
of a convolution integral, and this is the type of forced response that weighs
the values of the input at τ compared to those of the impulse response at −τ are equal to
product of transforms, convolution theorem.
In the last equality, the well-known property has been taken into account.
rema of the translation to the right
F(z)
Z[f(t−1)] = z
entries that admit rational transforms; in this case the expressions from an-
Titrasformare are rational functions that can be antitransformed without
It's difficult once the corresponding expansions in partial fractions are calculated.
The same considerations made for continuous time systems apply where
please keep in mind that
Z(A t ) = (zI−A)−1 z
Let's begin our analysis with the free response in the state:
xl (t) = At x0
what it implies, in
Xl (z) = (zI - A)−1 zx0
1 0
A= + ,
-1 0.5
in the calculation, it is advisable to refer to
aT
z−1 0 z−0.5 0
φ(z) + 1 z−0.5 , + z−1 ,
(zI - A)−1 = = -1
z (z−1)(z−0.5) (z−1)(z−0.5)
and do the expansion into partial fractions; in fact, multiplying by will have
expressions of the known type, z , of which the inverse can be calculated. In the example
z−λ
Hello studio,
1 0 0 0
+ , + +2 1 ,
φ(z) = −2 0 +
(z−1) (z-0.5)
120 5. The stationary linear representations in the complex domain
1 0 0 0
z+ , z+
-2 0 +2 1 ,
φ(t) = Z ) + * (R1+R20.5t )δ−1 (t)
(z−1) (z-0.5)
hay
1
C= ( 1 0 ) B= + 1,
and it is obtained through simple calculations
1
W(z) =
z−1
While the system is characterized by two natural modes and therefore two laws of
motto, only one is excitable and observable and that is the one associated with the way
constant with eigenvalue 1. Ultimately, we obtain a transfer function
ment characterized by a single pole.
The numerical integrator performs a sum of the values of the input signal.
just like the continuous integrator, characterized by a transfer function
1
imentoW(s) = , performs
s the integral
1 )=
ut()tẋ(
W(s) = ⇒ 8 W(t) = 1
s y(t) = x(t)
5.7. Discrete-time systems: analysis in the complex variable domain 121
t
y(t) = ! u(τ)dr
0
One observes the correspondence between the operations of discrete integration and with-
continue: the transfer function of a discrete integrator (summator)
at a pole equal to 1, the transfer function inside an integrator has a
polo in 0.
The following example shows that the same applies for discrete-time systems.
the considerations already made for those with continuous time regarding the features
filtering characteristics of a given system with respect to inputs whose transform
admits poles coinciding with some of the zeros of the transfer function.
Consider the system with transfer function
z+2
W(z) =
(z−0.5)(z+1)
yf(z) 1 z+2
=
z z - a(z - 0.5)(z + 1)
It remains evident that part of the response has the same time trend.
rale dell’ingresso. Ci`o non accade solo se si verifica una cancellazione con uno
zero of the transfer function. If it were, for example
z
u(t) = (−2)t⇒U(z) = z+2
it would be, coming out
yf (t) = r1 0.5t+r2 -1t
it deals with those parts of the entries whose transformation has one or more poles
coinciding with one or more zeros of the transfer function.
As a further example, the forced response to the unit step will be calculated.
the initial answer,
z
yf(z) = W(z)U(z) U(z) = u(t) = δ−1 (t)
z−1
yf(z) W(z) K r1 rn
= = + +. . .+
z z−1 z−1z−p 1 z−pn
N
|pi |<1 W=
D
n
yf(t) = Kδ−1 (t) + rI pit yr (t) = Kδ−1 (t)
(i=1
Similar behavior to continuous time: as time increases, if the poles
dellaW(z) all have modulus <1, the answer settles around
a constant trend,K,
W(z)
K= (z−1)|z=1=W(1)
z−1
detailed gain of the discrete-time system, that is, the value to which it tends
response to the growth of time when the transfer function has all the
police with module<1.
Perfect parallel therefore with continuous time systems where the gain is
W(0) and we find the most cited correspondence between 0 and 1 at the end.
Let's proceed with our analysis by noting that assigned
The value at zero can be calculated using the initial value theorem.
(in appendix)
f(0) = lim F(z)
|z| → ∞
5.7. Discrete-Time Systems: Analysis in the Complex Variable Domain 123
that exist for definite classes of input functions. Let's consider classes of
special entries
u(t) = sin(θt) = ejθt-−jθt
2j
it is obtained
t
yr (t) = W ( t − τ ) e jθτ
(
τ=−∞
5.7. Discrete time systems: analysis in the complex variable domain 125
and place
t-τ=ξ
∞ ∞
yr (t) = W ( ξ ) e jθ(t−ξ) =ejθt W ( ξ ) e −jθξ
(ξ=0 (ξ=0
∞
W ( ξ ) ( e jθ )−ξ =W(z)|ejθ
(xi=0
e−jθt W ( e jθ )
Note that the W(z) for z=ejθ it is in fact to be treasured most rare.
the radius of convergence of W(z) coincides with the maximum of the moduli of the
eigenvalues. Assuming these are all less than 1, on the circumference of
the transfer function is well defined.
From the previous calculation, adopting the polar representation for W
W ( e jθ ) W ( e −jθ )
⇓ ⇓
M ( θ ) e jφ(θ) M ( − θ ) e jφ(−θ)
and observing that they are even and odd functions of theta
it is obtained
that is
yr (t) = M(θ)sin(θt + φ(θ))
t(k) (t(t-1)(t-2)...(t-k+1))
and let us consider the canonical input of order (k) for discrete-time systems
t(k)
u(t) =
k!
t
yr (t) = lim W(t−r)u(r)
t0 →−∞ (
r=t 0
t
r (k)
= lim W(t−r) t-r=θ
t0 →−∞ (
r=t 0
k!
∞
(t−θ)(k)
= W(θ)
(θ=0 k!
Since for the factorial polynomial of order (k) of a binomial the expansion holds
following
k
k
(t−θ)(k) = i (k−i)
+ i , -1 t (θ+i−1)
(i)
(i=0
5.7. Discrete-time systems: analysis in the complex variable domain 127
yes it has
k ∞
1 k
yr (t) = -1i t(k−i) W ( θ ) ( θ + i − 1 ) (i)
(i=0 k! + i , (θ=0
k i
-1 t(k−i)
= Mi
(i=0 i! (k−i)!
in which it is placed
dI W ( z )
Mi-1i
dzi |
z=1
in fact ∞ ∞
di W ( z ) di W(t) equals negative onei(t+i−1)(I)
W(t)
=
dz I dz (t=0 z t
i (t=0 z t+1
it is deduced that
Miequals negative oneI i!cI
which shows a close formal analogy with that obtained for the systems at
continuous time.
The steady-state response can also be calculated as the part
persistent forced response. In this case, it is useful to refer to the
domain of transforms.
A first aspect concerns the calculation of the transform of the input canon-
factorial icon
t[k] t(k−1)...(t−k+1)
=
k! k!
128 5. The stationary linear representations in the complex domain
results
t[k] z
Z" =
k! # (z−1)k+1
To verify what has been stated, it is necessary to premise a property of the
Zeta transformation. More precisely
d
Z tf(t) = -zF(z)
& ' dz
indeed
and multiplying by -z
d f(1) 2f(2)
-zF(z) + +. . .=Z t ·f(t)
dz = z z2 & '
We can now calculate the Z-transform.
d z z
Z[t] = -z =
dz z^{-1} (z−1)2
and then the transformation of the factorial polynomial
t[k] =t(t-1)...(t-k+1)
d 1 2z
Z[t [2]=Z[t·(t−1)] =−z =
dz(z−1)2 (z−1)3
t[2] t·(t−1) d 1 z
Z9 : =Z 9 : =−z · =
2! 2 dz(z−1) 22 (z−1)3
t(t−1)(t−2) z
Z9 :=
3! (z−1)4
and, in general:
t[k] z
Z9 :=
k! (z−1)k+1
5.7. Discrete-time systems: analysis in the complex variable domain 129
n
yf(z) c0 c1 ck rI
= + +. . .+ +
z (z−1)k+1(z−1)k (z−1) (i=1 z−pi
n
t(k) t(k−1)
yf (t) = c0 +c1 +. . .+c k + rI pit
k! (k−1)! (i=1
and infinite
t(k) t(k−1)
yr (t) = c0 +c1 +. . .+ck
k! (k−1)!
with
W(z)
c0=K= (z−1)k+1 - =W(1)
(z−1)k+1 -z=1
-
-
gain of the discrete system
d W(z)
c 1= ) (z−1)k+1 *-
dz(z−1) k+1 -z=1
-
-
.....
1di
c i=
W ( z )-
i!dzi -z=1
-
-
Even in the case of the discrete-time system, we can repeat the same with-
previous considerations; the coefficients are nothing but the coefficients
of the Taylor series development of W(z) around the point z = 1.
130 5. The stationary linear representations in the complex domain
We can repeat here the same considerations already made in the case of time.
continuo
r mi
φ(z) RI
=
z (i=1
( k=1
(z−λi )k
z
φ(t) = Rik Z −1 ) *
( (i k
(z−λI )k
in fact
λ λ2
Z[λ t f(t)] = f(0) + f(1) + f(z) + . . . 2
z z
from which
t(k−1) z/λ z zλ (k−1)
Zλ9t := = λ
=
(k−1)! (z/λ−1)k λk (z−λ)k (z−λ)k
that is
t(k−1) z
Zλ9t−k+1 :=
(k−1)! (z−λ)k
e
r mi
t(k−1)
φ(t) = RI λt−k+1
(i=1
( k=1
(k−1)! i
These are the laws of motion that appear; therefore no longer just laws of motion.
of the type (λt ), but also laws of motion with coefficients that are polynomials in
(t). It corresponds to the terms that characterize the power of a matrix A
when this is not regular (the diagonal form does not exist). The same applies
5.7. Discrete-time systems: analysis in the complex variable domain 131
z−1 W(s)
W(z) = Z 9L−1 ) *- :
z s -
-t=KT
-
To understand this expression, note that the system is- time-discrete.
linear criterion, it is sufficient to calculate W(z) to express the ratio between
a forced exit and the corresponding entry; moreover, the step response to
discrete time coincides with the sampling of the unit step response
continuous-time system behavior (that is because the response to a discrete step,
from a continuous step); therefore the sampled continuous impulse response:
W(s)
L−1 ) *-
s -
-kT
-
coincides with the system's impulse response in- discrete time. Making it
Z-transform
W(s)
Z 9L−1 ) *- :
s -
-kT
-
-
132 5. The stationary linear representations in the complex domain
and dividing by the Z-transform of the discrete-time step function, we obtain, for
as observed, the sought transfer function
W(s) t2
L−1 ) *
s 2
W(s) k2T 2
L−1 ) *-
s - 2
-t=kT
-
-
−1 W ( s ) T2
9 L )
From *- := z(k 2 )
s - 2
-t=kT
-
- z
z(k) =
(z−1)2
d z
z(k·k) = -z
dz(z−1)2
−1 W(s) T 2z(z + 1)
9 L )
From *- :=
s - 2 (z−1) 3
-t=kT
-
-
5.7. Discrete-time systems: analysis in the complex variable domain 133
zT
zDi =
∼e i i = 1, ..., m
1
i = m + 1, ... , n - 1 zDI discretized zero of sn−m
1
it is important to note that for (n−m)>2 the zeros of the discretized sn−m
W ( s ) = e −λs H(s)
λT = eT · mT 0<m<1
134 5. The stationary linear representations in the complex domain
W ( s ) −It is emT s a
=e
s s(s+a)
z−1 1 mT s emT s
W(z) = z 9 L−1 ) e *- :
z zl s − s+a -
-t=KT
-
-
z−1
= l+1 zδ) −1 (kT+mT)−e−a(kT+mT) *
z
z−1 z e−amT z
W(z) = ) *
z l+1z−1z−e −aT −
Thus, one obtains a W(z) which is a ratio of polynomials, in reality one obtains
that the degree in the denominator is equal to (λ + 1). Therefore, there is a function of
transfer that is a ratio of polynomials and therefore a representation with
the state that is of finite dimension, which in reality would not occur in this
case for the continuous-time system.
The forced model coincidence in the complex domain means that the
the same solutions identified for continuous-time systems can be
apply AI systems to discrete time.
5.8. Appendix 4
5.8. Appendix 4 135
=h(s−λ
5i i ) (s2-2α 2 2
5 k s+αk+ωk)
k
2zk s2
p(s) = h(-λi ) (ω2nk ) (1 + τI s) ) 1+ s+ 2 *
5i 5k 5i 5k ωnk ωnk
αk2+ω2k=ω2nk
ωnk ζk=−αk
Starting from the form of the generic factorization with real coefficients of a
polynomial(s), it is obtained
2zk′ s2
h′ (1 +τi′ s) ) 1+ ′
ωnk s+ ′2
ωnk
*
4i 4k
W(s) =
2zk s2
h (1 +τi s) ) 1+ ωnk s+ ω2nk *
4i 4k
2zk′ ′
s2
K (1 + τi′ s) ) 1+ ′
ωnk s+ '2
ωnk
*
4i 4k
W(s) =
2zk s2
sr(1 +τi s) ) 1+ ωnk s+ ω2nk *
4i 4k
Starting from such a factorized form with real coefficients, the so-called form
Bode's factored form, it is relatively easy to understand what the behavior is
in frequency of the well-represented system, as noted from the response ar-
monica, W(jω). It is indeed necessary to remember that its modulus and its phase
they represent, in the case of all eigenvalues having negative real parts, the
modification in amplitude and phase that a pure periodic signal undergoes while passing through
5.8. Appendix 4 137
2z s2
K s 1 + τs 1+ s+ 2
ωn ωn
there are the four different types of terms that characterize the W(s); the term
s2
costanteK, monomios, binomio 1 +τs, trinomio 1 +2ζ s+ . The ωn product
ω2n
of this type, in the numerator and denominator, gives rise to a
genericW(s).
The Bode diagrams
Post
|W(jω)|=M(ω)
phaseW(jω) = φ(ω)
M(ω)
7
φ(ω)
The justification for such a choice is provided in the following steps where we
note that the diagrams of a product transfer function, this that
responds to the cascade of two systems, obtained by summing the diagrams
the indBe modules of the phases
W ( s ) = W 1 (s)W2 (s)
1 2
W=Mejφ =M1 ejφ ·M2 ejφ
1 2
W=M1·M2 ej(φ +φ )
So, starting from the graphic representations of modules and phases, the rep
The overall graphic presentation is obtained by summing for.
Based on the previous consideration, the opportunity of choosing the methods
the representation is validated by the fact that, as already noted, the W(s)
the product of four terms, constant, monomial, binomial and trinomial.
KdB= 20 log10|K|
K= 100
K=−1000
5.8. Appendix 4 139
100
K = -1000
50
K = 100
0
-50
-100
10-2 10-1 100 101 102
pulsation
200
100
K = 100
0
-100
K = - 1000
-200
10-2 10-1 100 101 102
pulsation
Figure 5.1
that is the module is a line with a slope of 20dB per decade, with a semi-
π
pre-increasing. The phase, however, is constantly equal to. A
2
monomial factor at
the denominator is represented by opposing trends. The figure shows
the trends
140 5. The stationary linear representations in the complex domain
100
AT THE DENOMINATOR
50
-50
A NUMERATOR
-100
10-2 10-1 100 101 102
pulsation
A DENOMINATOR
-100
Figure 5.2
MdB= 20 log101 + ω2 τ 2
.
1 √
ω= ⇒MdB= 20 log102 = +3dB
|τ|
1 √
ω≪ =
⇒MdB20 ∼
log 10 1=0
|τ|
1
ω≫ = omega +1020 log tau
∼
⇒MdB20 log 10
|τ|
1
Perωmuch less than |τ| the asymptotic diagram coincides with the axis of
abscissa; a line with a slope of 20dB per decade represents the value of
module perω≫|τ | Thsinilein1etsrecst atω= 1
|τ| the x-axis. The
The approximate representation made up of a broken line on two sides is called
asymptotic diagram.
5.8. Appendice 4 141
1 1
The phase perω≪ |τ| it is practically null, for ω>>> |τ| e p r a t ic a m e n t e
π
the fairy conτ>0
2 is equal to− 2conτ<0. π
The modulus and phase of binomial terms are represented in the figure.
numerator and denominator for positive values of τ.
100
At the numerator
50
At the denominator
-50
-100
10-2 10-1 100 101 102
pulsation
200
AT THE NUMERATOR
100
-200
10-2 10-1 100 101102
pulsation
Figure 5.3
1
With reference to a binomial term in the denominator, in ω = |τ| the module
presents a variation of (1/2)dB−3d√ relative to the zero value of
asymptotic diagram
As observed in the figure, the trend corresponding to values of τ
negatives do not change regarding the module, while it becomes opposite
regarding the phase. With constant term tracking
and of the term binomial we are able to characterize the behavior in
frequency of a system that, for example, has a transfer function type
K
W(jω) =
1 + jωτ
Based on what has been stated, the trend of the module and the phase is characterized by
from the sum of the trends modulo - phase of the two factors that in this case
142 5. The stationary linear representations in the complex domain
1
soneKe1+jωτ. The frequency ω= |τ| and that pulsation corresponds
of which the value of the gain of the module is −3dB compared to the initial value,
this pulsation is normally referred to as the bandwidth, B3 .
We have yet to address the representation of a trinomium, which
we will support the denominator, corresponding therefore to the transfer function
mind of a dynamic system characterized by a pair of eigenvalues com-
conjugated pairs. This will allow assigning a physical meaning to parameters
specify what the natural frequency and the damping are
1
(jω)2
1 +2ζωjω+
n ω2n
ω2 4z2
MdB-20 log10 ? + 1− , + 2 ω2
ωn
2 ωn
At ω=ωn The value of the module depends on the second:
MdB-20 log102|z|
and for the unit module, that which corresponds to the cancellation of the part
imaginary and the reduction of the trinomial square of a term
binomial, the value of the modulus results in −20 log 102 = −6
yes, therefore a line that has a slope of -40dB per decade, which crosses
the axis of literally at the pointn . It could be verified that the trend
effectively coincident with the asymptotic one at a decade of
distance from ωn .
60
40
SMORZAMENTO = 0.01
20
0
SMORZAMENTO = 1
-20
-40
-60
10-1 100 101
pulsation
Figure 5.4
144 5. The stationary linear representations in the complex domain
0
SMORZAMENTO = 0.01
-50
-100
SMORZAMENTO = 1
-150
-200
10-2 10-1 100 101 102
Figure 5.5
2 live
ωn
φ = arctan
1− ωω2n
2
100
50 A NUMERATOR
AT THE DENOMINATOR
0
-50
-100
10-2 10-1 100 101 102
pulsation
200
A NUMERATOR
100
AT THE DENOMINATOR
0
-100
-200
10-2 10-1 100 101 102
pulsation
Figure 5.6
Obviously, the trends of the phase for negative values of ζ are opposite.
to those indicated.
It is interesting to examine the frequency behavior of one of the systems.
physicists considered in the previous lessons.
Bode diagrams of the mechanical damper
1/K
W(s) = Ms2
k + bsk + 1
1 M
=
ω2n k
k b
ω n= / ζ= √
M 2kM
Seb >2kMlo√ damping would result in 1 which corresponds to eigenvalues
reality.
We evaluate such parameters as the mass and the elastic constant vary.
molla and the coefficient of dynamic friction.
146 5. The stationary linear representations in the complex domain
104 1
b = 104⇒ζ= √ = = 0.05
2kM 20
102
b= 102⇒ζ= √ 0.0005
2kM
Let's examine what happens in the two cases when the system is stressed with a
own periodic pulsation signal = ωn= 10 rad/sec
1
1
- (j10)2 10
6
|W(j10)|= -= = 10−1−1b
- +j+ b
1 - 10b
- 100 k -
- -
- -
That is, if we stimulate with a pulse of 10 radians per second with a
a force of 100 Newtons results, at steady state, in an oscillation of 10−3 m, 1 mm
(millimeter), in the first case, 10 centimeters, in the second.
(jω)2
100+ 1 +2ζ
ωn(jω) + jω2n ,
W(jω) =
jω(1 + jωτ1 (1 + jωτ2 )
50
-50
-100
10-2 10-1 100 101 102
pulse
200
100
-100
-200
10-2 10-1 100 101 102
pulsation
Figure 5.7
148 5. The stationary linear representations in the complex domain
100
50
-50
-100
10-2 10-1 100 101 102
pulse
200
100
-100
-200
10-2 10-1 100 101 102
pulsation
Figure 5.8
The considerations made regarding the plotting of Bode diagrams are ...
important not for the actual tracking of the diagrams themselves, what can
to be entrusted to a very simple calculating tool, but to understand
for inspection on the factors of the transfer function what will be the behavior
tamento in frequenza e, quindi nel tempo, del sistema stesso. Innanzitutto
it is necessary to remember that the corresponding elements depend on various factors
standard: to a monomial term, which corresponds to an integral action-
In the perspective, a fixed trend is associated that translates into concrete terms in a
constant decrease in gain at different frequencies and a phase shift
constant; to a binomial term corresponds a standard trend that
translate into a maintenance of profit and a subsequent constant decrease
decrease with a corresponding increase in phase shift; the time constant does not
it modifies the device's behavior, but only its 'placement' in fre-
similar considerations apply for a trinomial term that also
can exert a resonance or a filtering action (think of a term
trinomio a numeratore).
5.8. Appendix 4 149
1
W(s) =
s2
) 1 + 2ζωns+ ω2n *
ζ= 0.1 ω n= 2
2
Y(t)
1
0
U(t)
-1
-2
-3
-4
-5
20 21 22 23 24 25 26 27 28 29 30
Figure 5.9
150 5. The stationary linear representations in the complex domain
INGRESSO = SIN( 2t )
5
4
Y(t)
3
1
U(t)
0
-1
-2
-3
-4
-5
20 21 22 23 24 25 26 27 28 29 30
Figure 5.10
INGRESSO = SIN( 5t )
5
1 U(t)
Y(t)
0
-1
-2
-3
-4
-5
20 21 22 23 24 25 26 27 28 29 30
Figure 5.11
5.8. Appendix 4 151
We are around the natural pulsation of the system as shown in the diagrams.
in Bode reported below
BODE DIAGRAM
20
10
-10
-20
10-1 100 101
OMEGA
-50
-100
-150
-200
10-1 100 101
OMEGA
Figure 5.12
152 5. The stationary linear representations in the complex domain
BODE DIAGRAM
6
0
10-1 100 101
OMEGA
-50
-100
-150
-200
10-1 100 101
OMEGA
Figure 5.13
In the module diagram, there are two parameters that are often
assumptions to characterize the frequency behavior: the passband B3 ,
that pulsation from which there is an attenuation in the response to
permanent regime higher than 0.707 compared to the value in ω= 0; the modulus at
resonanceMr that is, the maximum value of the modulus divided by the value of the modulus
in zero.
Pronounced resonances highlight the fact that there are frequencies
which are exalted, while a wide passband characterizes the ra-
small time constants with which the system responds.
parameters and the previous simple considerations find a correspondence in-
mediated in the behavior of the response to the unit step. In fact, a wide
passband, response speed, translates into a rapid tendency to the val-
regime hours (low response time), while the enhancement of signals in a
certain frequency range translates into the presence of an enhancement of the
response in the presence of the initial discontinuous stimulus and for this reason rich
in harmonic content.
5.8. Appendix 4 153
ωn = 5 z = 0.5
the result is shown in the following figure. The previous considerations
They intend to convey the qualitative trend in relation to the different values.
the parametersn .
FORCED RESPONSE WITH INPUT = SQUARE WAVE
1.5
0.5
EXIT
0
-0.5
ENTRY
-1
-1.5
0 1 2 3 4 5 6 7 8 9 10
TIME IN SECONDS
Figure 5.14
154 5. The stationary linear representations in the complex domain
2
EXIT
1.5
0.5
-0.5
-1 ENTRY
-1.5
-2
-2.5
0 1 2 3 4 5 6 7 8 9 10
TIME IN SECONDS
Figure 5.15
10
τ 1= 1
(1 + jωτ1 )
20
10
0
10-1 100 101
Frequency (rad/sec)
-50
-100
-150
-200
10-1 100 101
Frequency (rad/sec)
Figure 5.16
156 5. Stationary linear representations in the complex domain
-1
-2
-3
-4
-5
-6
-7
-2 0 2 4 6 8 10
REAL PART
Figure 5.17
K
τ1 τ2>0
(1 + τ1 s)(1
+τ2 s)
5.8. Appendix 4 157
20
10
-10
-20
-30
10-2 10-1 100 101
Frequency (rad/sec)
-50
-100
-150
-200
10-2 10-1 100 101
Frequency (rad/sec)
Figure 5.18
-1
-2
-3
-4
-5
-6
-7
-2 0 2 4 6 8 10
REAL PART
Figure 5.19
158 5. The stationary linear representations in the complex domain
10
τ1= 1 τ2= 2
(1 + jωτ1 (1 + jωτ2 )
10
τ 1= 1 τ 2= 2 τ 3= 2
(1 + jωτ1 (1 + jωτ2 (1 + jωτ3 )
20
-20
-40
-60
10-2 10-1 100 101
Frequency (rad/sec)
-100
-200
-300
10-2 10-1 100 101
Frequency (rad/sec)
Figure 5.20
5.8. Appendix 4 159
-1
-2
-3
-4
-5
-6
-7
-2 0 2 4 6 8 10
REAL PART
Figure 5.21
(s−0.5)
W(s) = 40
(s+ 1)(s+ 4)
160 5. The stationary linear representations in the complex domain
20
10
-10
10-2 10-1 100 101 102
Frequency (rad/sec)
200
100
-100
-200
10-2 10-1 100 101 102
Frequency (rad/sec)
Figure 5.22
6 OMEGA POSITIVES
-2
-4
-6 OMEGA NEGATIVES
-8
-6 -4 -2 0 2 4 6 8 10
REAL PART
Figure 5.23
5.8. Appendix 4 161
In conclusion, note that if there are singularities on the imaginary axis, the
The polar diagram is open, meaning there are branches that go to infinity.
(sI−A)aT
(sI−A)−1 =
|sI−A|
En−1I
En−2=AEn−1+an−1 I
En−3=AEn−2+an−2 I
..
.
E0=AE1+a1 I
0 =AE0+a0 I
162 5. The stationary linear representations in the complex domain
an−1=−s1
2an−2−(s2+an−1s1 )
3an−3=−(s 3+an−1s2+an−2s1 )
...
Furthermore, since the eigenvalues of the power of a matrix are the eigen-
the eigenvalues themselves raised to the power, and the sum of the eigenvalues is equal to the trace
of the matrix, it follows that
si=tr(Ai )
It follows that the i−th element of Newton's relations can be expressed in the
following mode
1
ai=− tr(AEI )
n−i
⇓
EiI
⇓
Ci=AEi
⇓
1
ai= tr(AEi )
n−1
⇓
Ei−1=Ci+ai I
⇓
PRINT(i, E i , ai )
i= 0 ⇒AE+a I 0Stop
0 =
0 ∼
⇓
i = i - 1 ⇒ ITERA
subsection Appendix 9
164 5. The stationary linear representations in the complex domain
ρf>|z| orfndeidfeebalivroacm
xehplt ofoniuntcfasi tI
(ρfradius of convergence). The class of functions for which such an operation of
The transformation is defined and extremely broad, and it is indeed intuitive to understand.
It is said that for sufficiently large values of the modulus, the series is convergent.
Being the samples of the function multiplied by z. t
Very often, to indicate that the function which is being transformed is defined
pert≥0 if you use the convention to represent the function conf(t)δ−1(t)
5.8. Appendix 4 165
being−1the unit step, the even function equals one for t≥0. A prime
The aspect we consider is the calculation for elementary functions. Among these
let's consider the temporal expressions of natural modes that, as we know,
characterize the explicit description of the system. The law of motion of a
aperiodic or alternating mode of type λtand if there is
λ λ2 1 z
Z[λ t δ−1 (t)] = Z[λt ] = 1 + + 2 +. . .= λ =
z z 1− z
z−λ
defined by
λ
- z -<1⇔|z|>|λ|
- -
- -
- -
The verification is simple if one observes that the function in question is the sum of
a geometric series with reason and λzfor a geometric series with ratio
post
Sn= 1 + a + a2+a3+...+an
yes
1−an+1
S n=
1−a
and therefore
1-an+1 1
S=S∞ = limSn= lim = |a|<1
n approaches infinity 1−a
n approaches infinity 1−a
scent= ejθt-−jθt
2j
166 5. The stationary linear representations in the complex domain
jθt −jθt 1
Z[sinθt] = Z " e −e jθt −jθt
#= 2j Z[e ]−Z[e ] =
2j
1 z z 2 −jθ 2 jθ
= ) − z-e−jθ * = 1z -ze -z +ze =
2j z -ejθ 2j z −z(e +e ) + 1
2 jθ −jθ
zsenθ
Z[senθt] = |z|>1
z 2-2zcosθ +1
where the property of linearity was used.
ρf=max{ρf,ρf} 1 2
Right translation
1
Z[f(t−1)] = F(z)
z
In fact, from the definition:
∞
f(0)f(1)
Z[f(t−1)] = f ( t − 1 ) z −t = + 2 +. . .
(t=0 z z
it can be calculated
f(0)f(1)
zZ[f(t−1)] =z ) *
+ 2 +. . . =F(z)
z z
5.8. Appendix 4 167
and therefore
1
Z[f(t−1)] = F(z)
z
that is, shifting to the right by one step in time is equivalent to multiplying by in 1zthe
domain of the complex variable.
Left translation
f(2)f(3)
Z[f(t+ 1)] = f(1) + + 2 +. . .
z z
it can be calculated
1 f(1)f(2)f(3)
Z[f(t+ 1)] =−f(0) + f(0) + 2 + +. . .
z z z z3
and therefore
Note the analogy with the theorem of differentiation in the Laplace domain.
Lf(t)˙
+ =, sF(s)−f(0)
Convolution
t
Z" f(τ)g(t−τ) =
# F(z)·G(z)
(
t=0
t
Z) f(τ)g(t−τ) =
,
(
t=0
168 5. The stationary linear representations in the complex domain
1 1
=f(0)·g(0)+f(0)·g(1)+f(1)·g(0) + f(0)g(2)+f(1)g(1)+f(2)g(0) + . .,.
z+ , z2 +
g(1)g(2) f(1) g(1)g(2)
=f(0)g(0)
+ + + 2 +. . . +, g(0) + + 2 +. . . ,
z z z + z z
f(2) g(1)
+ + g(0) + +. . . ,
z 2 z
f(1) f(2)
=f(0)G(z) + G(z) + G(z). . .
z z2
f(1)f(2)
=f(0)
+ + + 2 G(z) =, F(z) · G(z)
z z
t−1
H(t−τ)u(τ)
(
τ=0
and parallel convolution in the field of time translates into the product of
transform yourself.
f(0) = limF(z)
|z|→∞
In fact
f(1)f(2)
F(z) = f(0) + + 2 +. . .
z z
e
limF(z) = f(0)
|z|→∞
f(∞) = lim(z−1)F(z)
z→1
5.8. Appendix 4 169
In conclusion:
lim(z−1)F(z)−f(0) =f(∞)−f(0)
z→1
To conclude, let's calculate the Z-transform of [Link] calculation follows the steps
turned in the case of scaling
A A2 A3
Z[A t ] =I+ + 2 + 3 +. . .
z z z
and place
A An
Sn=I+ + + . . .+
z zn ,
results
A A A An An+1
+ I− z,
Sn=I− +
z ,+
I+ + . . .+
z zn ,
I−
z n+1
that is
-1
A An+1
SnI− + + I−
z, z n+1 ,
and, in conclusion
-1 -1
A zI−A
Z[A t ] =limn⇒∞ Sn=I− + = + =z(zI-A)−1
z, z ,
which converges for
A
@z <1⇒|z|>∥A∥
@
@ @
@ @
@ @
6. The geometric properties of
state space
This At(xa-xb ) = 0
So, xa`e in d is t in g u ib ile d a xbit implies that (xa -x b ) it goes out identically
Nothing is in free evolution, that is to say, it is unobservable.
It is easy to show that the inverse also holds; in fact, if there exists an observation x¯
able, then ∀xa , xbthat
xa-xbx bar
6.1. The properties of the internal structure: Observability 173
xaindistinguishableb
Therefore, it is enough to study the existence of statix¯ that give free evolution nothing.
in output, unobservable states, to obtain indistinguishable states.
The questions to which we want to respond are:
Under what conditions is an assigned state x¯ observable?
How to characterize the set of unobservable states?
Sotto quali condizioni il sistema`e tutto osservabile (assenza di stati inosserv-
able) ?
If we denote with I the set of unobservable states, a first result
the following is important.
The set of unobservable states
dk At
This xt=0CAk eAtt=0x=CAk x k = 0, 1, ...
dtk . .
. .
That is ifAtx = 0 then CAk x = 0 for each flame, and in particular
C
⎛ CA ⎞
.. x=0
⎜ ⎟.
⎜ CA n−1 ⎟
⎝ ⎠
174 6. The geometric properties of state space
Or vice versa,
CA k x = 0 k= 0, . . . , n−1
from the Cayley-Hamilton theorem it follows that:
tk
This Atx= CA k x= 0
/
k≥0
k!
C
⎛ CA ⎞
ker ..
⎜ ⎟
.
⎜ CA n−1 ⎟
⎝ ⎠
C
⎛ CA ⎞
ρ= .. m<n
⎜ ⎟.
⎜ CA n−1 ⎟
⎝ ⎠
if this is true it can be demonstrated that:
A 0 B1
TAT −1 = 0 11 1 TB= 0 B2 1 C T −1=10 )
A21A22
176 6. The geometric properties of the state space
with
m
z1
z=Tx 0 z2 1
m−n
ż1=A11 z1+B1 u
⎧
⎪⎨ ż2=A21 z1+A22 z2+B2 u
y=C1 z1
⎪⎩
and it is evident the composition in two under system is1 ,S2with state variables
z1this2 each of the subsystems is influenced by the subsystem S1no
the influenza today2while influences2the exit is characterized only by the
subsystemS1 .
So in these coordinates the internal structure is that of the figure
6.1. The properties of the internal structure: Observability 177
z 10
S
z 1
S 1
u
z y
2
S 2
z 20
Figure 6.1
where it appears that only the subsystemS1influence the exit. Moreover, since the
ternaA11, B1 C1S available as to t u t e1it seems to be in a state of use.
this breakdown clearly highlights the internal structure of the system, which is
the interconnection of the two subsystems, one fully observable and the other not
The observable has a dimension equal to the rank of the observability matrix.
m n−m
..
T −1 = . base in I 7
6 ∗
178 6. The geometric properties of state space
1
A= 0 -2 C=(01)
0−1 1
C 0 1
0 CA 1 = 0 0−1 1
0 1 1
0 0−1 1
x= 0 implies I=gen 80 1 9
0
and choosing
..
T −1 = , 0.1 -
.
1.0 .
TAT matrices are obtained−1eCT−1of the desired form.
Another aspect that is wanted to be highlighted is that the bond
forced input-output is characterized by the only observable subsystem. This
It is evident from the figure, but it is directly obtained through calculation. We have
in fact
¯ eAt¯ B¯
W ( t ) = C eAt B = C
k
t k A11 0 B1
(C10 ) 1 0 B2 1 =
/ k! 0 A21A22
t k A11
k 0 B1
= 10 ) 0 k 10 B 1 =
/ k! ∗ A22 2
=C1 eA tB111
All and only the observable modes characterize the observable subsystem.
In fact, if we remember that the observable modes are those whose laws are...
they belong to AceAt, that is λoi
¯ At¯ (C1 eA
Ce= 11 t
0) =ψ(t)
T
*
R(T) = x:x= ! eA(T−τ ) Bu(τ)dτ, for someone0 , u [t0 ,T) +
t0
T
*
R(T) = x:x= ! eA(T−τ ) Bu(τ)dτ, u(−∞,T ) +
−∞
∗ ⊥
Having stated what is asserted in the lemma, being R(L) and Ker(L) are...
tospazi, is equivalent to
∗
R(L)⊥ =Ker(L)
what is true in that nonetheless presobinR(L)⊥ and for eachainA
The result now recalled can be applied to the case under study by observing
that the function
∗
! eA(T−τ ) Bu(τ)dτ
−∞
It is also observed that the adjoint operator of a matrix is the matrix ag-
matrix (which in the real case coincides with the transpose) and it is observed that for the
previous theorem theKerCeAtcoincides with the kernel of the observation matrix
Ability, it is not difficult to verify that the following result holds.
6.2. The properties of the internal structure: Reachability 181
Theorem 2.
T ⊥
∗
R(T) = R , ! - B∗ eA(T−τ )
eA(T−τ ) Bu(τ)dτ=Ker x =
−∞ 6 7
⊥
B∗
⎛ B ∗ A∗ ⎞
Ker .. =R(B AB· · · An−1 B)
⎜ . n−1 ⎟
⎜ B ∗ A∗ ⎟⎠
⎝
1−2 -1
00 =A 1 =B
31 0
-1
..
R=Im ⎛ -1 . 1 ⎞
..
⎝ -1 . 3 ⎠
In this case, the system is fully accessible and therefore one can go anywhere.
any state.
The verification of the reachability of a state has become very simple, because
In reality, given a state: is x reachable? It was necessary to see if x ∈ R (if not.
182 6. The geometric properties of state space
R=Im -1-1
0 0 01
in this case the set of reachable states is a space of one dimension
sione, this means that for this system all and only the states of the type:
α
0 01 =xR
R=Im .. ..
n−1
6 B.. . . .A B7
∥
6.2. The properties of the internal structure: Accessibility 183
¯
R(n)xR∋t≤n
For a continuous-time system, S, entirely reachable if we denote with x Rhe
generic state, this can be achieved starting from, statox(0) = 0 in
an arbitrary time, T. In fact, it can be verified that with the entry
′
uR (t) = B′ e−AtC −1 e−ATxR
in which
T
′
C= ! e−AtBB ′ e−Atdt
0
T
′
eAT ( ! e−Aτ BB ′ e−Aτ dτ)C−1 e−ATxR
0
eATe−ATxR=xR
therefore the state reached at time Ts under the action of the indicated input is
proprioxR .
Being arbitrary, it is said that the property of reachability is different -
exercises for continuous time systems. Obviously, the less time it takes to
go from 0 axRThe entrance to be applied will be wider.
Another observation concerns the fact that the states reachable from
zeros can also be forced to zero (controllability). A state, x, is
controllable aTse there exists an input that brings it to zero at time T. For
it can be shown that linear and stationary representations are controllable
all and only the reachable states.
If indeed a state, xTand reachable at T, being the set of states
reachable a subspace of the state space invariant with respect to A,
184 6. The geometric properties of the state space
even the state−e ATxTit is dragging; when you are in a bad mood,Ttale
what
T
eATxT+ ! eA(T−τ ) ThisT(τ)d(τ) = 0
0
θ
l
u
Figure 6.2
¨
mgsinθ = mlθ + ucosθ
θ
0 θ̇1 = 0
x1=θ
*
x2=θ ˙
=
xẋ 1 2
⎧ g 1
⎨ cosx1 u
dx2ot= sin x1−
I ml
⎩
calculate the linear approximation and the set of controllable states.
0 1 0
A= , g - . =J(xe ) B= , 1 cosx - ..
cosx10. . − ml 1
l . xe =0 . xe =0
. .
.
0 1 0
˙ , g
ξ= - ξ+ , 1 - u
0 − ml
l
0 1
− ml ⎞
R= = ⎛
6B AB 7 1 0
⎝ − ml ⎠
This matrix has a rank of two, so the set of reachable states is everything.
space, therefore also the zero state with zero angular velocity is controllable.
What implies the zero controllability of the physical system for small displacements
thoughts around balance (for more insight, refer to what has been said)
with stability analysis through linearization).
186 6. The geometric properties of the state space
u
Figure 6.3
..
ρ [Link].. ..
6 . . n−1A B 7 =m < n
A A B1
TAT −1 = 0 11 12 TB= C T −1= (C1C2)
0 A22 1 0 0 1
6.2. The properties of the internal structure: Reachability 187
with
A11(m×m) B1(m×p) C1→(q×m)
Furthermore, the triplet11 , B1 , C1a ggingittuert
This block structure highlights the fact that the system in
realt`a`e consists of two subsystems, one fully reachable and the other not reachable.
reachable.
z 10
S
u z 1
S 1
+
y
+
z 2
S 2
z 20
Figure 6.4
The structure of the matrices in the new coordinates corresponds to the structure
of the system highlighted in the figure. In fact, to a partition
m
z1
z=Tx= 0 z2 1
n−m
ρ(B· · ·An−1 B) =m
n−1
11 B1
¯ ¯ ¯ A ¯n−1 B̄) = , B1A11 B1· · · A
ρ(B|AB|... - =m
0 0 ··· 0
with this it is demonstrated that the pair (A11, B1 )and it is all reachable. So if from
given overall system (A, B, C), the triplet (A) is extracted11, B1 , C1 the new
representation thus obtained,
A,B,C(n)
⎧
⎪⎪ T
⎪⎪
⎪⎪ ¯ ¯ ¯
A,B,C(n)
⎪⎨
⎪⎪ ⇓
⎪⎪ A11 B1 C1 (m < n)
⎪⎪
⎪⎩
it is degrading.
As for the calculation of the transformation T
To this end, it is sufficient to choose the prime new coordinates a basis in the set
the reachable ones, and the remaining ones - any, but obviously such as to generate
independent vectors.
base ∀
T −1 = 0 .. 1
R .
m-column n-m
6.2. The properties of the internal structure: Accessibility 189
1 -1
B= 0 1 A= 0
1 0 1 1
..
⎛ 1.−1 ⎞ =gen 1
R=Im .. : 1;
⎝ 1 - 1 ⎠
..
1 0
T −1 = , 1.0 . - T=
0 1
1 .. 1 −1 1
1
TAT −1 = 0 −1−3 TB=
0 -1 1 0 01
dx1ot-x1-3x2+u
⎧
⎪⎨ ẋ−x
2 2
y=C1 x1+C2 x2
⎪⎩
Another observation concerns the fact that the forced link entry - us-
The cited is represented by the only reachable subsystem. This is also what
suggests the intuition in observing the figure. The calculation of W(t) confirms
such intuition. In fact
¯
¯ , ¯B ,¯ C ) ⇒
W ( t ) = C eAt B = ( A ¯ At
C e¯ B
ma
∞ k
t k A11 A12 B1
(C1 , C2 ) =
/
k=0
k! 0 0 A22 1 0 0 1
∞ k
t k A11
k
∗ B1
(C1 , C2 ) 0 k 1 0 0 1 =
/ k! 0 A22
k=0
establishes a precise link between the reachability of the states and their excitability
natural methods.
As already mentioned, they are excitable with input impulses those natural modes whose
exponential laws appear in:
H(t) = e AT B
If one takes into account the fact that the exponential laws that characterize the
two subsystems1eS2 what inTAT1 they are the same as those of the matrix
A
A A
A→λI TAT −1 = 0 11 12
0 A22 1 →λI
It is understood that the eigenvalues of A are partly in A.11 , in number dim, e
part inA22 , in number of (n−m); furthermore from the calculation
¯ eA 11tB1
H(t) = e AtB ¯= 0 1
0
cio`e inH(t) ci stanno soltanto le leggi esponenziali che stanno ineA tB1 11
those that correspond to the eigenvalues of A11 , that is to say those that correspond
to the subsystems1 It could be shown that since the subsystem
S1everything is reachable all its ways are excitable. This allows to affir-
but the reachable subsystem is characterized by all and only the modes
excitable.
A final observation is needed to clarify what is meant when one states
that the properties of reachability and observability are dual to each other.
This stems from the fact that if one considers the two systems
A
=ẋ+
B
xu
*
y=Cx
̇′
=zAx+C′u
*
y=B′ u
the duality of each other, the complete reachability of one is equivalent
the complete observability of the other and vice versa. Therefore, the existence of a
6.3. The Kalman decomposition 191
Property data on one of them implies the existence of dual property on the
dual system.
A A B1
TAT −1 = 0 11 12 TB=
0 A22 1 0 0 1
C
⎛ CA ⎞
I=ker ..
⎜ .⎟
⎜ CA n−1 ⎟
⎝ ⎠
if the system is not fully observable, we know that there exists a transformation
in coordinates, T, such that
A 0
TAT −1 = 0 11 1 C T −1=10 )
A21A22
192 6. The geometric properties of state space
It will be shown below that it is possible to select the coordinates in such a way that
The properties of the states should be highlighted.
3
X = R
I
Figure 6.5
z S
u 1
S 1
z 2
S 2
+
z 3
S 3 y
+
z 4
S 4
z
Figura 6.6
The subsystemS1and influence today2 , daS3and also toS4and infinite from the entrance
u. The subsystemS2and influenza today4 , and from the entry. The subsystemS3
and influence today4 , edS4The autonomous modal evolves into the method of achievements.
The output is conditioned by S2eS4 .
In conclusion, this is the structure that appears from the outside in the coordinates.
S1It is a system that can be reached, but it is not servable;2 the suns are shining
theme characterized by reachable and observable states; The subsystem S3no
reachable and unobservable; finally, S4it is a non-reachable system
but observable.
It appears from this decomposition that well manifests the internal structure of the
196 6. The geometric properties of the state space
state, that the forced connection input output is characterized by only the sub-dose
themes2 .
This can be verified directly with the calculation. In fact, in the coordinates
z
W ( t ) = C e AtB=C2 eA tB222
and this for the block structure of matrices. Therefore, the impulse response
of the overall system is the response of the subsystem S2what is a system
everything is reachable and observable.
It is easily understood, and it would not be difficult to prove, that the trio
(A22 , B2 , C2 ), the one that characterizes the subsystem S2 , has the largest dimensions
small among all the triples of matrices that have fixed impulse response and
equal to W(t). Taking advantage of the connections that exist between property of rag-
reachability and observability of states, and excitability and observability of modes,
despite the considerations made regarding the implementation, the subsystem
S2extract a minimal achievement from a possible realization.
A u = iR1+iC 1 C1 V˙C=−i
1 R+u 1
6.3. The Kalman decomposition 197
1 VC1=VR+V1C⇒VR=V 3
C+VC1 1 3
VR1
iR 1 = = VC1−VC 3
R1 R1
VC1 VC3
C1 V˙C=−
1 + +u
R1 R1
˙ 1 VC1 VC3 u
→VC=− + +
C1 R1 C1 R1 C1
+
+ i +
i v
i
- +v - -
- y
+
i v
- + -
Figure 6.8
B iR=i
2 C+u 2
˙ 2
C2 VC=iR −u 2
VR2
VC3=VR+V2C 2 ⇒ = VC3-VC 2
R2 R2
VC2 VC3
C2 V˙C=− +
2
R2 R2 -u
˙ 2 VC2 VC3 u
→VC=− +
C2 R2 C2 R2 − C2
Infinite perx3you can write the current balance at node C:
C iR=i
1 C+iR 3 2
i C3 =−i+iR2 R1
198 6. The geometric properties of state space
˙
-iR+i2R=C3 V1 C 3
VR1
iR 1 = = VC1-VC 3
R1 R1
iR 2 = VC3-VC 2
R2
VC1 VC3 u
V =− + +
⎧ C1 C1 R1 C1 R1 C1
⎪⎪ VC2 VC3 u
⎪⎪ V =−
⎪⎪ C2 +
⎪⎨ C2 R2 C2 R2 − C2
VC1 VC2 VC3 1 1
⎪⎪ VC3 = + , + -
⎪⎪ C3 R1 C3 R2 − C3R1 R2
⎪⎪
⎪⎩
y=VC-V1 C 2
1 1
− R1 C1 0 R1 C1
⎛ 0 1 1 ⎞
− R2 C2 R2 C2
A=
⎜ ⎟
⎜ 1 1 1
, 1
+ 1
- ⎟⎟
⎜ R1 C3 R2 C3 − C3 R1 R2
⎝ ⎠
1
C1
C= ( 1 −1 0 ) B= ⎛ − C12 ⎞
⎝ 0 ⎠
6.4. The properties of the state space 199
⎪⎪ χ4=0
⎪⎪
⎩
200 6. The geometric properties of state space
6.4.a. Unobservability
6.4.b. Accessibility
¯ 0 , x0 ,u¯¯ )[t=x
ϕ(t−t 0 ,t) R t¯ < ∞
so:
¯
φ(t,0,u¯ [t 0t̄,) ) =xR
6.4. The properties of state space 201
so it can go, varying(0), in all the states that are in the image
diB; still
by varying (0) I (1) in all possible ways, we obtain all the states that
they are in the image of Be of AB, essentially:
n×p n×p
..
R(2) = Im 6 B . AB 7
That is, the matrix that is obtained when in the characteristic polynomial to the variable
If a matrix A is replaced, it is equal to the null matrix.
1 1 1
A= 0 , B = 0 1
1 11 1
that characterize a system certainly not entirely reachable, but control-
unstable. It is indeed evident that the input(0) = x1 (0) +x2 (0) brings the state to
zero in a single moment in time. Even in this case, the two properties are
equivalents are in the solar.
7. Linear input-output models
and representations with the state
under what conditions does this functional connection coincide with the forced response
of a finite-dimensional, stationary linear system, causal?
Since we know the constrained bond of a linear, stationary system to
finite dimension, characterized by:
t
yf(t) = ! A(t−τ ) B + Dδ(t - τ)u(τ)dτ
" This #
0
the formulated problem admits a solution if the nucleus, K(t), coincides with the
impulse response matrix of a system.
One thus arrives at the following reformulations of the problem in the domain
of time and the complex variable, respectively.
Problem of realization
AssignedK(t) (q×p) find (n,A,B,C,D) such that
The constructive proof, that is, it will be shown how starting from a matrix
of proper rational functions, it is possible to construct, by inspection of the coefficients
of the matrices that describe it, a realization of it.
Assume, without loss of generality, that K(s) is a matrix of functions.
strictly proper rational functions are those rational functions when the degree
the numerator is lower in degree than the denominator. It is said without
loss of identity because one can always shift from a proper form to
a strictly proper, in fact, a rational eigenmatrix can always
scrivere come
K(s) = K 0+K(s)
conK0constant matrix eK(s) strictly proper rational matrix. Given
D=K0 the problem of realization reduces to the calculation of din, A, B, C for
a strictly proper rational matrix of the type
D K(s)
1 1
s+1 0 s+1
K(s) = " s+1 # = $ %+ $ 1
%
s 1 s
0 I 0 ··· 0 0
⎛ .. .. .. .. ⎞
. . I . . ⎛ ⎞..
.. .
.. ..
AR= ⎜⎜ . . . 0 ⎟
⎟ B=R ⎜ ⎟...
⎜ ⎟ ⎜ ⎟
⎜ 0 0 I ⎟ ⎜ 0⎟
⎜ ······ ⎟ ⎜ ⎟
⎜ ⎜ I⎟
⎝ -a0 I ··· −an−1 I ⎟⎠ ⎝ ⎠
The identity matrices, I, are matrices that have dimensions (p×p). and, of
consequence, the matricesAR ,BR eCRthey have size (np×np), (np×p)
e (q×np), respectively.
Example:
s 0 1
1 " s+ 1 # " 1# +" #1 s
s+1
K(s) = " 1 # = s2+s =
s s2+s
0 1 0 0 1
AR= " # BR= " 1# CR= " 1
0-1 1#
in fact, it is easy to verify that
Another example:
1 0 1 0
" # + " 1 2#
s
K(s) = -1 1
1 + 3s + s2
208 7. Linear input-output models and state representations
in which
p= 2 n= 2
0 0 1 0
0 I ⎛ 0 0 0 1⎞
AR= " # =
-a0 I -a1 I ⎜ -1 0 -3 0⎟
⎜ ⎟
⎝ 0−1 0−3 ⎠
0 0
⎛ 0 0⎞ 1 0 1 0
B=R CR = " #
1 0 -1 1 1 2
⎜ ⎟
⎝ 0 1⎠
The demonstration of the correctness of the procedure followed to calculate
The implementation consists of verifying that the matrix of transfer functions
The measure associated with such representation coincides with the assigned K(s).
In this regard, it is necessary to calculate
T
(sI−AR )a
CR (sI−AR )−1 BR=CR BR
(sI−AR )
-1
s -1 0 ··· 0
⎛ .. ⎞ 0
.. ⎛ ⎞..
0 . . .
⎜ .. ⎟
= (B0 . . . B n−1 ) ⎜ .. ⎟ ⎜ ⎟...
⎜ . . 0 ⎟ ⎜ ⎟
⎜ ⎟ ⎜ 0⎟
⎜ ⎟ ⎜ ⎟
⎜ 0 -1 ⎟ ⎜ 1⎟
⎜ ⎟ ⎝ ⎠
⎜ a0an−2s+an−1 ⎟
⎝ ⎠
We refer to case p=q=1 for simplicity of notation. So
observe, first of all that the calculation of the determinant det(I−A) Rd'a
What is always true for a matrix that has the structure of AR the
The elements of the last row are precisely the coefficients of the characteristic polynomial.
7.1. The problem of realization 209
∗ ∗ 1 0
⎛ ∗ ∗ s ⎞ ⎛ ⎞..
.. .
(B0. . . B n−1 )
⎜ ∗ ∗ . ⎟ ⎜ ⎟... =
a0+a1 s+. . .+sn ⎜ .. ⎟ ⎜ ⎟
⎜ . ⎟ ⎜ 0⎟
⎜ ∗ ∗ ⎟ ⎜ ⎟
⎜ ⎟ ⎜ 1⎟
⎝ ∗ ∗ sn−1 ⎠ ⎝ ⎠
0
⎛ ⎞..
.
B0+B1 s+. . .+Bn−1 sn−1
$∗ ∗ % ⎜ ⎟... =
a0 +a1 s+. . .+sn ⎜ ⎟
⎜ 0⎟
⎜ ⎟
⎜ 1⎟
⎝ ⎠
and so
B0+B1 s+. . .+Bn−1 sn−1
K(s) =
a0+a1 s+. . .+sn
0 1 0
AR= " # BR= " 1# CR= (b0. . . b 1)
-a0−a1
ẋ=x
1 2
⎧
⎪⎨ dx2ot=−a 0 x1−a1 x2+u
y=b0 x1+b1 x2
⎪⎩
210 7. Linear input-output models and state representations
b1
. .
u x2 x2 =x1 x1 + y
+ - ! ! b0
+
-
a1
a2
Figure 7.1
Note that the achievement identified according to this technique has diminished.
sione (np×np); therefore every time the number of inputs is one, in-
depending on the number of outputs, the size of this achievement
ne.
The realization (AR , BR , CR this in canonical form reachable in
how much the state space associated with it is all reachable. In fact, to
start from the expression
0 1 0 0
⎛ .. .. ⎞ ⎛ ⎞...
AR= . . BR =
⎜ 0 1 ⎟ ⎜ 0⎟
⎜ an−1 ⎟⎠ ⎜ 1⎟
⎝ -a0 ··· ⎝ ⎠
CR= (b0. . . b n−1 )
the reachability matrix is easily calculated
0 0 0 1
⎛ .. .. .. ⎞
. . .
(B|. . . A n−1 B) = ⎜ ... 0 1 ⎟
⎜ ⎟
⎜ 0 1* ⎟
⎜ ⎟
⎜ 1* * ⎟
⎝ ⎠
that has a determinant different from zero and therefore full rank and even adn. This
prove that the realization (AR , BR , CR ) it is all reachable.
7.1. The problem of realization 211
Note that if the acquired results are linked to what is shown in the
study of reachability and observability properties, we have a
calculation procedure that allows us, given a generic realization of
a given nucleus, to calculate a new one that satisfies the property of being entire
reachable. If in fact (A, B, C) is a not fully reachable realization,
0· · · · · · -a
0 0I
⎛ .. ⎞ B0
I . ⎛ .. ⎞
Ao .. Bo .
= ⎜⎜ 0 . . . .
⎟
⎟ = ..
(nq×nq) ⎜ . ⎟ (nq×p) ⎜ ⎟
⎜ .. .. .. ⎟ ⎜ . ⎟
⎜ . . ⎟ ⎜ ⎟
⎜ 0· · · 0 ⎝ Bn−1 ⎠
⎝ I -an−1 I ⎟⎠
Co
(0 · · · · · · 0I)
(q×nq)
aT
s 0 ··· 0 +a0
⎛ .. .. .. ⎞ B0
s . . .
.. ⎛ .. ⎞
(0···01) ⎜-1 0 .. .. ⎟ .
. . 0 . ⎟ ..
a0+· · ·+sn ⎜ .. .. .. .. .. ⎟ ⎜ ⎟
⎜ . . . . . ⎟ ⎜ . ⎟
⎜ ⎟ ⎜ ⎟
⎜ ⎟ ⎝ Bn−1 ⎠
⎜ ⎟
⎝ 0· · · 0−1s+a n−1 ⎠
by virtue of the structure of the matrix Co the only thing we want to calculate
and the last term, as the preceding ones multiplied by zeros.
B0
∗ ∗ ∗ ⎛ .. ⎞
(0···01) ⎛ ⎞ .
a0+. . .+sn ∗ ∗ ∗ n−1 ⎜ .. ⎟
⎝ 1 · · · s ⎠ ⎜ . ⎟
⎜ Bn−1 ⎟
⎝ ⎠
0 0 -3 0 1-1
⎛ 0 0 0−3 ⎞ ⎛ 0 1⎞
Ao= Bo=
⎜ 1 0 -2 0⎟ ⎜ 2 0 ⎟
⎜ ⎟ ⎜ 1 0 ⎟⎠
⎝ 0 1 0−2 ⎠ ⎝
0 0 1 0
Co = "
0 0 0 1#
(q×p) q= 1
( 1 −1 ) + ( 2 0 )s
K(s) =
3 + 2s + s2
0−3 1−1
Ao= " # Bo= " Co( 0 1 )
1−2 2 0#
7.1. The problem of realization 213
0−a0 b
Ao = " # Bo = " 0 # Co= ( 0 1 )
1−a1 b1
dx1ot-a0 x2+b0 u
⎧
⎪⎨ dx2ot=x1−a1 x2+b1 u
y=x2
⎪⎩
u
b0 b1
. .
+ x1 x1 + x2 x2 y
! !
- -
a0 a1
Figure 7.2
0 0 1
Co ⎛ 0 0 1 ∗⎞
⎛ Co Ao ⎞ 0
.. = ⎜ ⎟
⎜ . ⎟ ⎜ ⎟
⎜ ⎟ ⎜ ⎟
⎝ Co Aon−1 ⎠ ⎜ ⎟
⎜ ⎟
⎝ 1 ∗ ⎠
The observability matrix is upper triangular and is therefore non-singular.
the realization is therefore entirely observable.
214 7. Linear input-output models and state representations
0 I 0
.. .. 0
AR ⎛ ⎞
= . . BR= ⎛ ⎞..
(np×np) ⎜ 0 I ⎟ ⎝.⎠
⎜ ⎟ I
⎝ -a0 I · · · · · · ⎠ I
−an−1
0· · · · · · -a
0 0I
⎛ .. .. ⎞
I . . B0
Ao .. .. .. ⎞
= ⎜⎜ 0
..
. . .
⎟
⎟ Bo = ⎛ .
(nq×nq) ⎜ .. .. .. ⎟
⎜ . . 0 . ⎟ ⎝ Bn−1 ⎠
⎜ ⎟
⎜ 0· · · 0 I -an−1 I ⎟⎠
⎝
having these realizations with dimensions equal to when one goes to realize a
nucleus that characterizes an input-output relationship with a single input p = 1
or with a single output.q= 1
Further elementary examples. The supplement:
1
K(s) = ⇒A=
s 0 B= 1 C= 1
=ẋyux
K K/τ 1 K
K(s) = = 1 ⇒A=− B= C= 1
1 + τs τ +s
τ τ
1 K
=ẋ−x+=
yux
τ τ
the double integrator:
1
K(s) =
s2
0 1 0
A= " B= " # C= ( 1 0 )
0 0# 1
dx1ot=x22ẋ=u y =x1
1 0 1
A= " # B= " # C= ( 1 0 )
-1-2 -1
the transfer function is:
1
W(s) =
s−1
and the construction of its realization according to one of the two techniques already
you see from
Ā= 1 B̄= 1 C̄= 1
A first observation concerns the proposed techniques which, as already pointed out,
highlighted provide realizations of dimensions in particular conditions. In
in this case such achievements are also minimal.
7.2. The minimum achievements 217
q= 1⇒(Ao , Bo , Co minimum
p = 1 ⇒ (AR , BR , CR minima
s2+ 2s 0
1
0 " 2s + 2s2+ 1 #
s+1
K(s) = " 2 1 # = =
s(s+2) s+2 s3+ 3s2+ 2s
0 0 2 0 1 0 2
" 2 + s+ " s
1# " 2 0# 0 1#
2s + 3s +s
2 3
p=q= 2
0 >N= 6
n= 3
1 0
s+1
⎛ ⎞ K (s)
K(s) = = " 1
2 1 K2 (s) #
⎜ s(s + 2) s+ 2 ⎟⎠
⎝
There are therefore two connections
(10)
K1 (s) =
s+ 1
realizable through matrices
1 ..
0 .
⎛ s+1 ⎞ ⎛ . ⎞
K(s) =
2 1
= K1 (s) .K2. (s)
⎜ s(s + 2) ⎜ ⎟
⎝ s+ 2 ⎟⎠ ⎜ .. ⎟
⎝ . ⎠
To each of these columns, I can associate a realization of size
3 perK1
s(s + 2)
" 2(s + 1) #
K1 (s) =
s(s + 1)(s + 2) →n1= 3
7.2. The minimum achievements 219
and 1 perK2
0
" 1#
K2 (s) =
s + 2 ⇒n2= 1
the overall realization will have dimension N=n1+n24 ed`e of the type
A1 0 B1 0
A= " o B= " o C= (Co1Co2)
0 A2o # 0 Bo2 #
It has been taken into account that the system results from the connection of two
subsystems that have the same outputs and each has a single input.
With reference to the example of multiple inputs and outputs under examination, the technique of
The implementation by rows is more convenient than that by columns.
In conclusion, with the information known to us, we are able to build
of achievements of a size smaller than what would turn out to be the
size of the realization built directly on the K(s).
greater than 1*
1 0 0 1 0 0 0
0 " 1 0# " 0 0# " #
s+ 1
K(s) = ⎛ 2 1
⎞ = + + -1 1
s s+1 s+ 2
⎝ s(s + 2) s+2 ⎠
r
Ki
K(s) =
1i=1 s−pi
Ki=Ci Bi
(q×ρi )(ρi×p)
A realization that turns out to be minimal is that of dimension N =
r
i=1ρiwith
2 Diagonal to elements coinciding with the poles of the expansion in fractions
partial with multiplicities each equal to the rank of the related residue
BeC, the aggregate of matrices BieCi
B1
⎛ ··· ⎞
. .. ..
B= ⎜ .. ⎟ C=
⎜ ⎟ 3 C1.· · · Cr .4
⎜ ··· ⎟
⎜ B ⎟
⎝ r ⎠
In the case of the previous example:
0 0 0
K 1= " = " # (10)
1 0# 1
1 0 1
K 2= " = " # (10)
0 0# 0
0 0 0
K 3= " # = " 1# (−1 1)
-1 1
and therefore
0 0 0
A= 0⎛−1 0⎞
⎝ 0 0-2 ⎠
1 0
0 1 0
B= ⎛ 1 0⎞ C= "
1 0 1#
⎝ -1 1 ⎠
This technique allows us, in the case of simple poles of the expansion in
partial fractions of K(s), to identify a minimal realization for another
7.3. Introduction to the study of Interconnected Systems 221
One last aspect that wants to be highlighted is that the achievements min-
they have dimensions equal to the rank of the following matrix, H, called the matrix of
Hankel. The result, which would not be difficult to prove, is important in the
applications as it allows to detect the order of the minimum realizations
eventually starting from the examination of entrance - exit behavior.
K(t)
Si=q×p
(q×p)
S 0S 1S 2· · ·
⎛ ⎞
S1S2S3· · ·⎟
H= S S S · · ·⎟ ..
⎜2 3 4 ⎟
S 3 S4 .
⎜ . .. .. ⎟
⎜ . ⎟
⎝ . . . ⎠
infinite matrix, in blocks, where
I
dK(t)
S i= 5
dtI 5t=0
5
5
The Hankel matrix from a realization (A, B, C) of a system
dynamic, can be obtained by replacing at K(t) the expression
of the impulsive response, obtaining for the coefficients:
Si→CAi B
The interconnected systems, that is, made up of several parts, of several subsystems among
They connected in various ways are the subject of this lesson.
The first problem that arises is the calculation of the model that describes how
specifically the interconnected system starting from the knowledge of the models
222 7. Linear input-output models and state representations
of the individual parts that make it up. The following introduces the
elementary connections: in series, in parallel and feedback.
Consider two systems, the first S1who has entered1the exit1 ; it
seconds2conu2hey2 A first connection that can be considered is that
the exit of the first becomes the entrance of the second; thus there is an interconnection
specified by the following conditions:
u2=y1 u=u1y=y2
ẋ=A
1 1 x1+B1 u12
ẋ=A2 x2+B2 y1
0 0
y1=C1 x1 y2=C2 x2
Therefore, it is possible to find the overall model by introducing into the descriptions of
each of the two systems the condition of interconnection, in this case:
ẋ1 A1 0 x1 B1
" ẋ2 # = " B2 C1A2 # " x2 # + " 0 # u
x1
y = ( 0 C2) " x2 #
It can therefore be noted that the aggregate of the two systems is still a system, and it
state space can be considered through the aggregate of the components
of the state space of the first system and the second. This is a connection
important the cascaded connection or series connection is reported in
figure
u1 y1 u2 y2
S1 S2
Figure 7.3
7.3. Introduction to the Study of Interconnected Systems 223
u1=u2=u y=y1+y2
ẋ=A
1 1 x1+B1 ux˙2=A2 x2+B2 u
0 0
y1=C1 x1 y2=C2 x2
y=y1+y2
ẋ1 A1 0 x1 B1
" ẋ2 # = " 0 +
A2 # " x2 # " B2 #
u
x1
y = (C1 C2) " x2 #
u + u1 y1 y
S1
-
y2
S2 u2
Figure 7.5
It is about connecting the two systems.1 edS2in a different way. One takes the
systems1with entered1he/she/it exited1and the systems2conu2hey2in sense
reverse; start from a defined external entry regardless and proceed
to subtract the exit from it2of the second system:
u1=u−y 2
Then consider the overall output of the system as the input of the system.
S2the exit1 , that is:
u2=y1 y=y1
An interconnection is built in which the first system is acted upon by the ...
the difference between an external input and its output brought back and subtracted,
after passing through the system2 at the entrance. There is a cycle that
it is called counterreaction, because the negative sign is introduced, meaning we subtract
draw the exit1upon entry. This feedback or cycle system
closed, or closed-loop is very important especially in the description of
natural phenomena, as well as in the technological and application fields.
The overall system, starting from the descriptions of the two systems.1 , S2 ,
here are the equations:
ẋ=A
1 1 x1+Bu−B1 C2 x2
0
y=C1 x1
dx2ot=A2 x2+B2 C1 x1
0
y2=C2 x2
7.3. Introduction to the study of Interconnected Systems 225
ẋ1 A1 -B1 C2 x1 B1
" 2ẋ # = " B2 C1 +
A2 # " x2 # " 0 #
u
x1
y= (C10 ) " x2 #
Three emblematic cases have been examined, a system with two parts in
cascata o in serie, un sistema con due parti in parallelo, e uno con due parti in
a counter-reaction scheme, it has been verified that assuming it as a state
The aggregate of the states of the component parts can describe a system.
of the same nature, with the finite, linear, stationary, regular dimension.
Ultimately, it has been found how the representation of an interconnected system
more parts can be identified starting from the models of the parts that constitute it
of course the interconnections can also be more complex,
So perhaps it is appropriate to introduce simpler graphic notations of
which of the block diagrams used so far.
Representations of interconnected systems through flow graphs
In the block diagrams considered so far, the variables are fast-
presented by arches that connect blocks; the block represents the system that
connect the two input - output variables; the connections occur for col-
ligament, derivation, and sum of arches. There are four symbols, four
notations through which the block diagram is constructed: arcs, blocks, sum-
mothers and letters.
Flow diagrams represent an alternative tool for describing
the interconnected systems. In them, the variable is represented by a point, a
nodou, the system through a segment that links two variables u−y, there is
a certain dualism this time the trait represents the system, ultimately the fact
that the same variable can be input or output of multiple systems, represents
simply the branching from a generic variable. The sum instead
is assumed following the convergence on a node, establishing that the vari-
able to be found at the point of convergence is the sum of the effects of the two
systems represented by the two segments. The notation is therefore very simple.
and allows for quick representation of interconnections. It has nothing of
conceptually different from the symbolism of the block diagram, but it is very
sleeker.
226 7. Linear input-output models and representations with state
The symbols have well-established names, the points are called nodes and represent
variable feelings, the traits are called branches, and the set of nodes and branches
it constitutes a flow graph, that is, the graphical representation of the system. In
In some cases there are open paths, parallel paths, however there are -
three characteristic configurations, namely rings, closed loops. Open paths
the rings that in a flow graph allow for the development of a theory for
describe the overall system in a very simple way.
The transfer function of the overall system
Cascade system case: be W1the transfer function that connects
entry - exit, eW2entrance - exit of the second system.
y1=W1 u1 y2=W2 u2
y=W2 W1 u=Wu
W=W1 W2
y= (W1+W2 )u=Wu
W=W1+W2
y=W1 (u−W2 y)
W1
y= u=Wu
1 +W1 W2
W1
W=
1 +W1 W2
This last expression is very interesting because1and the function of
transfer of the direct path connecting entrance - exit; 1 +W1 W2the
The denominator can be interpreted as 1 minus the product of the functions of
transfer of the ring (−W1 W2 Therefore, two concepts can be introduced.
important:
P=W1
L = -W1 W2
D = 1 - L = 1 + W1 W2
W1 P P
W= = =
1 +W1 W2 1−L D
The function has a direct interpretation. Suppose we divide the node,
from the case to counter-reaction, in two parts, a node where there are starting branches
in this case one, and the other with branches arriving, that is to cut the knot where
there is the loop, moreover it is assumed to insert a signal equal to 1 in this loop
interrupted. What is the signal at the end of the interrupted cycle?
1·(-W1 W2 ) =L
Hy=u,
i∆i Pi
S= 2 ∆
in which
PItransfer of′ I walk from u to y
∆⇒1−(−1)k+1 Pjk
1k1 j
S4
u 1 S1 S2 S3 1 y
S5
S6
Figure 7.6
P1=S1 S2 S3
0
P2=S3 S4
∆= 1−S2 S5−S2 S3 S6
in this case one must stop, because the rings touch, they have a branch in
common, but just one shared node is enough to no longer be considered
the product of the pairs, which therefore becomes zero, just like all those that
they follow.
∆ideprived of all those products of transfer functions that
involve some transfer function that is already in PI
∆1 = 1
0
∆21 - S2 S5
7.3. Introduction to the study of Interconnected Systems 231
P1 ∆1+P2 ∆2
S= = S1 S2 S3+S3 S4-2 S3 S4 S5
∆ 1−S2 S5−S2 S3 S6
Another example that is discussed is shown in the figure
S6 S8
u one S1 S2 S3 S4 S5 1 y
S7
S9
Figure 7.7
ẋi=Ai xi+Bi ui
S i: 0
yi=Ci xI
the calculation boils down to imposing the topological constraints established by the connections between
the various subsystems. As already seen, in the case of the cascading connection, one
ha
232 7. Linear input-output models and state representations
ẋ=A
1 1 x1+B1 u
⎧
⎪⎨ ẋ=A
2 2 x2+B2 C1 x1
y=C2 x2
⎪⎩
in which the state is the aggregate of the states of the subsystems. The heuristic procedure
which consists of enforcing compliance with constraints is simple in elementary cases,
but when connections with multiple subsystems are considered, with multiple
reverse reaction connections, it is necessary to refer to a procedure
systematic based on the analysis of the flow diagram of the interconnected system
connection. This procedure is based on the fact that in the absence of cycles, the calculation of
representation does not present difficulties. It is therefore a matter of releasing some
connections in order to obtain an acyclic graph, calculate its model and
then restore the open connections.
Consider the flow graph and make, possibly in number
at least, separations of nodes so that there are no more in the resulting graph
cycles. The separation must be done in such a way as to leave all branches
incoming on one side and outgoing branches on the other. For a graph of this type,
without loops, the representation with the state is immediately calculated, to
starting from that of the subsystems.
In fact, the state vector is the aggregate of the state of the individual subsystems.
the representation can be simply written by inspection on the graph
of flow. This will be highlighted by referring to the example in the figure
S4
u 1 S1 S2 S3 1 y
S5
S6
Figure 7.8
p= well, the part of the node separated where all the branches arrive
7.3. Introduction to the study of Interconnected Systems 233
A=
+
B
xẋu1 1 1 1
⎧
⎪⎪ ẋ=A 2 2 x2+B2 s
⎪⎪
⎪⎪ dx3ot=A3 x3+B3 C2 x2+B3 C4 x4
⎪⎪
⎪⎪
⎨Axuẋ 4
=
+
B 4 4 4
⎪⎪ ẋ=A
⎪⎪ 5 5 x5+B5 C2 x2+B5 C4 x4
⎪⎪
⎪⎪ ẋ=A x +B C x
⎪⎪ 6 6 6 6 3 3
⎩
y=C3 x3 p=C1 x1+C5 x5+C6 x6
⎪⎪ ẋ=A
⎪⎪ 5 5 x5+B5 C2 x2+B5 C4 x4
⎪⎪
⎪⎪ ẋ=A x +B C x
⎪⎪ 6 6 6 6 3 3
⎩
y= ( 0 0C30 0 0 )x
p= (C10 0 0C5C6)x
To restore the original graph, it is necessary to impose that it coincides with p.
By doing this, the representation with the state of the system is obtained.
detergent.
In conclusion, the overall representation is obtained: separating
some nodes of the graph, calculating the associated representation, imposing the
coincidence constraint of separate nodes
So two methods: one for calculating the transfer function, the other
for the calculation of the representation with the state, they put, among other things, in evi-
it is assumed that the state of the overall system is the aggregate of the states of the sub-
systems.
234 7. Linear input-output models and state representations
1
W(s) =
s−1
+
=
uxẋ
0
y=x
A system of dimension one has been obtained while the size of the system
comprehensive, interconnection of two subsystems of size one and two.
This fact can be explained in terms of loss of reachability or
of system observability, because the transfer function of the system
overall and influenced only by the fully reachable subsystem and everything os-
servable.
In reality, the case being considered corresponds to a loss of reason.
accessibility of the system. This can be verified by realizing the two subsystems,
finding the representations with the state, the overall representation and
Finally verifying that the system is not all reachable.
An opposite cancellation will correspond to a loss of observability.
p,(7.1)
)ut(,)xt(=
f)tẋ(
! "
y(t) = ηx(t), u(t), p (7.2).
! "
0 = f(xe ,0)
A
=
+
B)xutẋ(
The equilibrium states are those that satisfy the following equation:
0 = Axe
xandand it is stable
∀ϵ∃δ(ϵ) :
$ x2i %
#i=1 |xi | #i=1 max|x
i i|
A
=)xtẋ(
z(t) = αx(t)
ż()=
tαẋ()=
tαAx()=
tAz()t
therefore the new variable satisfies the same dynamic equation that it satisfied
lax, then having the same dynamic properties, I will be satisfied, for example,
the same condition of local asymptotic stability namely:
δa
∥x(t0 )−xe ∥< α ⇒ limα∥x(t)−x e= 0
t approaches infinity
Another aspect concerns the stability of a motion: let us suppose, again for
graphic representation requirements, the two-dimensional state space1 , x2
one starts from an initial state x(t)0 ), and imagine fixing an entrance; one will have
a motion that satisfies the equation:
)ut(,)xt(=
f)tẋ(
! "
if the initial state is disturbed, a different trajectory is born and a motion for-
disturbed,xp (t), which satisfies the same equation:
ẋ(t)
p = f xp (t), u(t)
! "
it indicates the evolution of the deviation between the two trends
ξ(t) = x p (t)−x(t)
it is obtained
Naturally, in the general case now being discussed, the system is not stationary.
because it depends on the particular input u(t) and the motor x(t). But in the case of the
stationary linear systems do not happen. In fact:
A
=
+
B)xutẋ(
7.3. Introduction to the Study of Interconnected Systems 223
u1=u2=u y=y1+y2
ẋ=A
1 1 x1+B1 ux˙2=A2 x2+B2 u
0 0
y1=C1 x1 y2=C2 x2
y=y1+y2
ẋ1 A1 0 x1 B1
" ẋ2 # = " 0 +
A2 # " x2 # " B2 #
u
x1
y = (C1 C2) " x2 #
||x(t)|| ≤ k||x(t0 )∥
ϵ
||x(t0 )∥< k
ϵ
||x(t)|| ≤ k k
So if the matrix φ is limited to evolution free limited; furthermore,
se∥φ(t−t 0 if it were not limited, there would be at least one element of the matrix
not limited to the changing of:
With the same arguments, it is easy to verify that a necessary condition is...
aware of asymptotic stability and that
lim∥φ(t−t 0 )∥= 0
t approaches infinity
A
=)xtẋ(
lim|φ(t−t 0 = 0
t→∞
244 8. Elements of stability theory
The imaginary axis is the boundary that divides the eigenvalues of the matrix.
There are two categories, those that give rise to converging natural ways and those
which give rise to divergent natural ways. In conclusion, the condition of
stability can be expressed by saying that the eigenvalues of A must not
positive, as long as they are simple; for asymptotic stability, they must all be with
Negative real part. Therefore, the imaginary axis is the stability boundary.
So given the dynamic matrix that characterizes the system, the stability
asymptotic or not, it is linked to the sign of the real part of the eigenvalues of A. This
and the condition that must be verified through some methods, as can be seen
they do not require the calculation of eigenvalues, but only the sign of their part
real.
The stability criteria, which we will discuss in the upcoming lessons, are used to verify
such allocation.
where it is indicated with the prime 1 or >1 the unitary or greater multiplicity
of an eigenvalue.
8.2. The Routh criterion 245
The Routh criterion allows us to determine if all the zeros have a real part.
negative and potentially how many are the zeros in the positive real part. The criterion
it consists of constructing a table and then verifying the coefficients
degli elementi della prima colonna della tabella. Il criterio in esame estende
the Cartesian rule on the correspondence of solutions with positive real parts and
negative to variations and permanences of the signs of the coefficients of a polynomial
Of second degree. As it will be remembered, a variation corresponds to a zero.
a positive real part, to a permanence one with a negative real part.
ax2+bx+c= 0
n an an−2 an−4
n−1 an−1 an−3 an−5
bn−2 bn−3 ...
cn−3 cn−4 ...
anan−2
' an−1an−3 '
bn−2 '
= ' '
'
' -an−1 '
anan−4
' an−1an−5 '
bn−3 = '' '
'
' -an−1 '
The following line, of order -3, is still constructed with the same procedure.
but with reference to the two previous lines; and so on
Example:
λ5+ 3λ4+ 2λ3-2λ2+ 2λ + 4
5 1 2 2
4 3 -2 4
3 8 2 0
3 3
2 11 4
− 4
1 50
11
0 16
In the construction of the Routh table, the result does not change if it is modified
you can multiply the rows by the same positive number; this allows for
simplify the calculations in the construction of the table.
Once the construction of the table is completed, the first column is considered.
table and they count the number of sign changes that occur in the transition
from one coefficient to another. The number of variations is equal to the number of zeros
the positive real part of the given polynomial.
8.2. The Routh criterion 247
In the case under study, the fifth degree polynomial has two zeros with real parts.
positive, the other three have negative real parts.
It is important to note that one may encounter obstacles in construction.
of the table; this occurs when an element of the first column becomes zero.
This can happen for only the first element of the row or for all the elements.
of the line. When one of these situations arises, stability is not achieved.
asymptotic.
The impediment can be overcome in three different ways as explained in the se-
guito.
Replacing the symbol ϵ with zero to represent a positive number and 'small-
color
1 2 3
2 4
ϵ 3
-4ε-6
ϵ
3
Once the table is constructed, the sign changes are counted. In the case
There are two sign changes in the exam, so it is concluded that there are two.
zero positive real part.
Multiplying the given polynomial by a binomial with a negative zero;
d1 (λ) =d(λ)(λ+ 1)
1 3 4
2 6
0 4
from which
1 5 10
3 9 4
6 26
-1 1
32
1
A different situation that can occur when building the table is nullification.
say an entire line. This can only happen in correspondence with a line of
odd order, in fact the two previous rows must be proportional and
8.2. The Routh Criterion 249
so they must have the same number of elements (note that in passing from
an odd line to an even one underlying the non-changing elements.
In this circumstance, it can be concluded that the polynomial under consideration is the
product of two polynomials; the first will have zeros with a real part characterized
given the sign variations of the elements in the first column of the table
so far built (the zeros of)1the positive real part are as many as the variables
sign actions that appeared in the first column of the constructed table
up to that moment); the second polynomial is of even degree and equal to the index
of the line preceding the line that was canceled, it only has even degree powers
and its coefficients are in order, from the one of maximum degree to that of
zero degree, the coefficients of the row that precedes the one that has been annulled. If
for example consider
5 1 3 1
4 1 3 1
3 0 0
The given polynomial is therefore equal to the product of two polynomials.
the first of first degree with a negative zero, the second of degree four
characterized by the presence of powers of λ only of even order
For such a polynomial, we can calculate the zeros, but in the case of the construction of the
if the table had stopped at a higher level, there could be some
250 8. Elements of stability theory
difficulty in seeing the sign of zeroes did2 (λ). Then the construction of the
Table can be resumed in another way. The polynomial is derived.2
d 4
(λ + 3λ2+ 1) = 4λ3+ 6λ
dλ
and the coefficients of the null row (that of order three in the case under study) are
sostituiscano i coefficienti di questo nuovo polinomio. Si pu`o or a c on t in u ar e
the construction of the table.
5 1 3 1
4 1 3 1
3 2 3
2 3 2
5
1 3
0 2
Upon completing the table, the analysis of the coefficients is carried out. In the case under study, it does not
there are sign variations, which guarantees stability in the case in question
simple, as it could be verified by calculating the zeros. It should be noted
the zeros of the polynomial2 (λ), with powers of only even degree, have a
double symmetry, with respect to each axis of the complex plane. This axis-
care that if there are no zeros on the positive real axis, all are found on the axis
imaginary.
Example:
λ2+ 9λ+ 20 = 0
8.3. The general method of Lyapunov* 251
It has a zero on the negative real part as a result of the elementary application of
Routh's criterion or, which is the same in the case of a second-degree polynomial
degree, of the Cartesian criterion. Are the zeros with real part less than -3?
A tale question it is still possible to respond by applying the Routh criterion to
polynomial obtained by substituting aλ, (λ−3)
The latter has zeros with negative real parts, which implies zeros with real parts.
less than -3 in the original polynomial.
and applying the Routh criterion to the following polynomial with real coefficients
the number of sign changes divided by two indicates the zeros that are located at
outside of our sector.
252 8. Elements of stability theory
)xt(=
f)tẋ( f ( x e) = 0
xe1
.. ⎞
xe= ⎛ .
⎝ xin ⎠
n
V(x) = 2
#i=1 (xI-xhi)
In the plane, n = 2,
(x1-xand)2+1 (x2-xe )2 2
V(x) ≥ 0
d ∂V dx
V(x(t)) =
dt ∂x dt·
ẋ1 n
∂V∂V ∂V . ∂V
$ ... % ⎛ .. ⎞ V̇(x) = fi (x)<0
∂x 1∂x 2 ∂x n · #i=1 ∂x I
⎝ dxnot ⎠
The Lyapunov criterion says that essentially, if you are able to find a
function V(x) defined positive in Ixsuch that
e its derivative along the motion is
definite negativity in the same neighborhood then one can conclude about stability
local asymptotic of the equilibrium state.
If you have talked about asymptotic stability, in reality the criterion also allows for
establish under what conditions the equilibrium state is locally simple
˙
stable mind, in this case it is sufficient that (V) is negative semidefinite.
ẋ-x 2
1 1 x2
&
ẋ=−x 3
2 2
↓
254 8. Elements of stability theory
x
xe= , -
0
0
, 0- =xestate of equilibrium
V(x) = x12+x22
-2x21x22-2x42=−2x22(x12+x22)≤0
A positive semidefinite function is obtained, therefore the equilibrium state xe
baitssiehretytilaernI .ebaltsm
uym
ispilbrilqeuiofetatsasi ti
simple global of this state, because it is a generalization of the criterion
following: if the function V(x) is positive defined throughout the space, also the
function V(x) is radially unbounded, that is, when x tends to infinity
the function tends to infinity, and the function (V) is˙ defined as negative or semi-negative
negative in all space, then there is also simple, non-local stability
my global.
Another physical example we will address concerns the dynamics of the configuration.
a rigid body in the absence of gravity. This model is assumed to describe the
variation of a satellite's orientation with respect to an inertial reference.
We denote with ω1the angular velocity around the axis x, ω2that around
to the assembly, eω3quella attorno all’asseze immaginiamo che l’orientamento del
satellite can change as a result of pairs, τ1 ,τ2 ,τ3 , which are generated
around the principal axes of inertia, which we will assume to coincide with the axes
of the reference system. The equation of dynamics can be written
J1 ω̇=
1 (J2-J3 )ω2 ω3+τ1
J2 ω̇=(J
2 3−J1 )ω3 ω1+τ2
J3 ω̇(J
3 1-J2 )ω1 ω2+τ3
oveJ1 J2 J3These are the moments of inertia around the axes x, y, z. This mod-
Hello generalizes orientation with respect to an axis: one can imagine having a
8.3. The general Lyapunov method* 255
metallic square in the plane, free from friction and that can be operated with a
pair around their center of mass, the equation of dynamics is
τ=ω̇
τ = -kω
yes it has
−
=ω̇
kω
& ω(t)→0
ω(t) = e −kt ω(0)
So it is enough to apply a torque proportional to the angular velocity,
A dynamic friction would be sufficient.
The same result occurs in the general case. Indeed, if one sets
ω→0. The verification of what is stated is easy if one substitutes into the equations of the
dynamics of the rigid body
τi=−k i ωi
1 1 1
V(x) = J1 ω21+J2 ω22+J3 ω23
2 2 2
Taking the derivative, it is easily verified that it is a defined function.
negative throughout the space with reference to the state ω = 0. In conclusion
V(x) is positive definite in the whole space with respect to the state ω = 0 and is
radially unlimited, and the (V) ˙that is calculated is a negative defined function
in all space, and this is enough to guarantee global asymptotic stability
of the zero state. Therefore, the body tends to stop.
256 8. Elements of stability theory
A
=)xtẋ( xe= 0
V(x) = x′ P x > 0
In reality, when considering quadratic forms, if one does not assume the matrix
symmetric and it can be easily verified that the same result is obtained that
it would obtain if the symmetric matrix obtained fromPaccording to the
so-called symmetrization procedure:
P+P T
P̃=
2
in fact
P+P T
x′ P x = x ′ x
2
As an example, it is immediately verified that:
1
1 1 1 −2 x
x′ , - x=x
′
,
-2 2 −
1
2 2 -
This explains why one limits oneself, in quadratic forms, to assuming matrices.
symmetric. The positivity is then verified with Sylvester's criterion.
With the previous premises on quadratic forms, we can state the cri-
Lyapunov's theorem for stationary linear systems that provides a condition
necessary and sufficient for the verification of asymptotic stability.
The condition requires that the matrix equation
A′ Q + QA = -P
the symmetric and positive definite unknown has a unique solution for each
A symmetric and positive definite matrix.
So fixed to referee P, that in the applications of the criterion is assumed
the identity I, the equation in the unknown Q must be solved; such a solution must
to be unique, symmetric, and positive definite. If so, it can be affirmed
that the system is stable, otherwise it can be asserted that the system is not stable.
The consideration is therefore stronger than in the non-linear context in which
the impossibility of concluding from a positive definite function does not
allows declaring the absence of ownership.
258 8. Elements of stability theory
V(x) = x′ Qx > 0
′
V̇(x) = x˙Qx + x Qx˙′
V̇(x) = x A′ Qx′ + x QAx′
V̇(x) = x(A′ Q +′ QA)x
-x′ P x < 0
z=eAT x
The criterion provides a necessary and sufficient condition based on the resolution.
n(n+1)
of a system of 2 equations. There are as many equations as there are
unknown in a symmetric matrix of dimension n. For example, with
0 1
A= , -
-5-6
the matrix equation
−1 =−5q12-5q12
0 = -5q22+q11-6q12
-1 = 2q12-12q22
These are three equations, in general.n(n+1), pern=
2 2.
post
ξ(t) = x(t) - x e
provides
∂f 1 ∂f 1
⎛ ∂x 1 ··· ∂x n ⎞
˙ = Aξ(t)
ξ(t) A=J(x) =e ..
⎜ . ⎟
⎜ ∂f n ∂f n ⎟
⎜ ··· ⎟
⎝ ∂x 1 ∂x n ⎠
As we will see, the asymptotic stability of the linear approximation implies the
local asymptotic stability of the statee . It follows that it is necessary to verify
preliminarily the state of equilibriumehe is not late.
It is possible to demonstrate that the state of equilibriumefor the nonlinear system
he is this or that, so he
|J(xe )|̸= 0
This condition is not necessary as shown by the following simple example:
-3ẋ
x=
love, despitee = 0, the equilibrium state is unique (in reality the solution
it is the algebraic multiplicity of the eigenvalue.
To show that the asymptotic stability of the linear system implies stability
local asymptotics of the stateeof the given nonlinear system, let V(ξ) = ξ′ Qξ,
with solution of A′ Q + QA = -P, a Lyapunov function for the system
linear approximating. The proof of the local asymptotic stability of the system
nonlinear is based on the use of the same function
positively defined function with respect to xe , which has a definitively negative derivative in
an appropriate contextx. e
This means that the asymptotic stability of the linear system implies that one
local asymptotic stability of the equilibrium state for the nonlinear system. Furthermore, it
it can show that the instability of the linear system due to the presence of at least one
8.4. Stability through linearization 261
Eigenvalue with a positive real part implies instability of the equilibrium state.
in the non-linear system, while nothing can be said in the case where the system
linear and in stability limit conditions or unstable due to mere presence
of eigenvalues with zero real part and geometric multiplicity greater than one.
This result concludes, as far as we are concerned, sufficient conditions.
stability coefficients. For example, if the linear approximation had a matrix
0
A= , -1
1−3 -
-1 0 0 0
J ( x e ) =A= , orJ(xe ) =A=
1 0- , 1 0-
nothing could be said about the stability of the equilibrium state regarding
to nonlinear dynamics.
Here are some simple application examples.
Dynamics of the pendulum
post
θ=x1 θ̇=x2
ẋ=x
1 2
262 8. Elements of stability theory
k g
ẋ=−x
2 2+ sin x1
m l
The calculation of equilibrium states of a
0
xe= ,
hπ -
It is noted that the equilibrium states 2hπ correspond to the pendulum in the
vertical position downward, while the equilibrium states (2h+ 1)π
correspond to the pendulum in the upper vertical position in a condition of
´ that the two types of equilibrium states have
initial velocity zero. It is clear
different properties: the vertical one, when disturbed, tends to be abandoned
Nato, the other conversely tends to be maintained if disturbed. The method
the linear approximation provides this result to support intuition.
The linear approximation corresponding to each of these states of
equilibrium can be calculated from the Jacobian matrix
0 1
J ( x ) = $ -g %'
cosx1 ' -k
l ' xe m
'
'
In the equilibrium states corresponding to the upward vertical position, one
ha
0 1
A= , g
-k -
l m
in the others, however:
0 1
Ā= $ -g %
-k
l m
With the dynamic matrix of this form, the study of stability is immediate.
as the elements of the last row represent the coefficients of the
characteristic polynomial changed signs. Therefore, in the first case
0 1 k g
A= $ g % →|λI−A|=λ 2+λ− characteristic polynomial
-k m l
l m
8.4. Stability through linearization 263
0 1 k g
Ā= $ -g ¯
% →|λI−A|=λ+λ+ 2
-k m l
l m
with two eigenvalues with negative real parts, i.e., asymptotic stability of the approximation
linear and thus local asymptotic stability of the equilibrium position
responding to the pendulum downward.
Logistic equation
The following nonlinear differential equation
2-)=
kx(t)
axt()tẋ(
x(a−kx) = 0
xe10
⎧ a
⎨ xe2 =
k
⎩
The stability analysis method through linearization can be used.
to study the stability of these two equilibrium states. The Jacobian matrix
in this case it is the scaling:
J ( x e ( ) = (a−2kx)x e
264 8. Elements of stability theory
e
per xe2if you want(xe2 ) =−a
The result is now evident: being a positive constant, the first state of
a
The equilibrium is unstable, while the state of equilibrium
k is stably asymptotic.
The previous considerations find further confirmation in the example that
next.
Predator-prey model
As already seen, the model suitable for representing the interaction between two species.
ypethet
ẋ(t)
1 = ax1 (t)−bx1 (t)x2 (t)
&
dx2ot(t) =−cx2 (t) +dx1 (t)x2 (t)
dovex1ex2they represent the density of the two populations and it has been assumed that
the prey population,x1 , grow, in the hypothesis of unlimited resources according to a
exponential growth with exponent a, while the predator population,
x2 , since the prey is its only nourishment, it decreases in the absence of interaction.
with the prey itself according to an exponential trend characterized by
coefficient - c. The coefficients a, b, c are all positive.
x1 (a - bx2 ) = 0
&
x2 (-c+dx1 ) = 0
The application of the linearization method does not allow one to conclude,
how the eigenvalues of the approximation can be easily verified
lines have zero real part. Yes, therefore in a condition where one cannot
8.4. Stability through linearization 265
infer nothing about the stability of the equilibrium state of the nonlinear system.
& x1 (a−bx2-kx1 ) = 0
x2 (-c+dx1 ) = 0
⎧ (a−bx2-kx1 ) = 0 c
⎨ (−c+dx1 ) = 0→x1=
d
⎩
c
, a−bx2-k =0
d-
a kc
x2= −
b db
c
d
, a kc - =xe
b − bd
therefore, a new state of equilibrium appears which is positive when ad - kc > 0
or, what is the same, ka > cdThsi conoidnti si expersed by saynig htat hte
the carrying capacity of the prey species must be sufficiently high to
can react to the predator's prey capacity. This state of balance is
positive and has a physical meaning, in the context being examined, on
the reader is invited to reflect.
The application of the method with respect to this state of equilibrium allows for verifying
about its local asymptotic stability. The steps are left for exercise.
mathematically explains well a phenomenon that occurs in the presence of
Competition between two species, prey - predator, occurs in reality.
Following a disturbance from a given equilibrium condition
an oscillatory phenomenon arises that converges to the state of equilibrium.
266 8. Elements of stability theory
At the completion of the study conducted, after having developed criteria and
internal stability methods, let's consider the input behavior
exit.
Under what conditions on the limited input system does it correspond
a limited exit? This property is also called external stability
A
=
+
B)xutẋ(
&
y(t) = Cx(t) x(0) = x 0
8.5. The Entry-Exit Stability 267
For all M,
0 there exists an N.
x ,M:|u(t)|< M⇒0 ∥y(t)∥< Nx ,M
two cases are distinguished: the first of stability in the zero state, related to
forced response when that ex0 = 0; a second of external stability in
any state, if this property is true ∀x0 this simply stability
external.
The conditions
t
2. 0 ∥W(τ)∥dτ< k2For all t
0
The sufficiency of the indicated conditions is easily demonstrated starting from the expression
of the answer
t
y(t) = ψ(t)x 0+ 0 W(t−τ)u(τ)dτ
0
||a+b|| ≤ ||a||+||b||
268 8. Elements of stability theory
it is obtained:
t
Still,
The norm of the dot product of a and b is less than or equal to the product of the norms of a and b.
implies
t
||y(t)|| ≤ ||ψ(t)|| ||x||0 ∥+ ∥W(t−τ)u(τ)∥dτ
0
0
t
||y(t)|| ≤ ||ψ(t)|| ||x||0 ∥+ ∥W(t−τ)∥∥u(τ)∥dτ
0
0
remembering that
||u(τ)|| < M
it concludes
t
||y(t)|| < ||ψ(t)|| ||x0 ∥+ ∥W(t−τ)∥dτ·M
0
0
one obtains the limitation of the response for each initial state. If the state
initial is null the first part is not there and the same applies to the limitation of the y, only
for hypothesis 2.. Therefore condition 2 is sufficient for external stability
in a zero state.
The necessity is demonstrated with similar arguments, fundamentally the
The demonstration consists of verifying that, if conditions 1 and 2 are not met.
they are satisfied, unlimited output evolutions can be obtained. As far as
Note that if the matrix ψ(t) does not have a bounded norm, some among
its elements are not limited to what corresponds to the existence of evolutions
unlimited free. Therefore condition 1.`is also necessary. Similarly
If it is assumed that hypothesis 2 is not satisfied, it can be demonstrated that there exists
a limited entrance, and it is that entrance that manages to take those terms
which are possibly not limited, corresponding to an evolution
8.5. The Stability of Entrance - Exit 269
forced in unlimited exit. Essentially, the two conditions are necessary and
sufficient. The formulation is given in terms of matrix norms. It is possible
understand how such conditions reflect on the characteristics of the system
particularly on the eigenvalues, this allows for establishing a connection between the different
types of stability studied.
ΛEo ⊂C−
=
uẋ
&
y=x
In fact: ∞
W ( t ) = 1 →0 dt=∞
0
This responds to the fact that if a constant input is allowed into the integrator,
the integral tends to infinity as the transfer function has a zero
270 8. Elements of stability theory
in zero.
Λ1o⊂Ce− o ⊂C−
Λ>1
e
ΛEo ⊂C−
oveCe− indicate the closed negative complex half-plane and the subscript 1 e>1 the
one is greater than one respectively. The second condition has already been
Commentary corresponds to 2.. The first corresponds to 1. How easy it is.
check remembering that in free evolution out comes the auto-
values associated with observable exponential laws.
Λ={λ1 , . . . ,λ n }
This ATΛo={λo1, . . . ,λ or }
eATBΛE={λE1 , . . . ,λ Es }
This ATBΛEo =ΛE∩Λo
To understand how the first of the given conditions corresponds to the [Link]
note that the boundedness of the norm of ψ(t) corresponds to the boundedness
of its elements. These elements are a linear combination of exponential functions.
essentially, possibly with polynomial coefficients in the presence of eigenvalues
a geometric multiplicity greater than one. They are limited only if all the
the eigenvalues that appear to you, those that are observable, are with strictly real part
negative, if the geometric multiplicity is greater than one, with a negative real part
or nothing for others.
So condition 1 is equivalent to the following
Λ1o⊂Ce− o ⊂C−
Λ>1
8.5. Entry-Exit Stability 271
L is a subset of C−
ΛE0⊂C−
Λ10⊂Ce− Λ>1
0 ⊂C− Λ0 ⊂C−
E
ΛE0≡Λ
That is, all modes are excitable and observable. In fact, in such a case, the eigenvalues
which appear in W(t), which by hypothesis have a negative real part, are all
those of the dynamic matrix which is, therefore, asymptotically stable.
External stability0= 0 implies external stability if
ΛE≡Λ
272 8. Elements of stability theory
ΛE≡Λ
ΛE0=Λ∩Λ0
ΛE0≡Λ
that is, if all the modes are excitable and observable. If we limit ourselves to considering
the implication between external stability and non-asymptotic internal stability, then
the condition of observability of the modes is sufficient. In fact, if all the modes are
observables may eventually exist in the spectrum of eigenvalues at
real part zero, and therefore of constant or periodic motion laws. And this
it implies simple internal stability. In conclusion, external stability implies
the asymptotic internal stability is
ΛE0≡Λ
Λ0≡Λ
-1 2 0 1
A= ⎛ 0 0 0⎞ B= ⎛ 0⎞ C= ( 1 1 0 )
⎝ -1 1 1 ⎠ ⎝ 0⎠
8.5. The Stability of Ingress - Egress 273
λ1-1
λ2= 0
λ3equals 1
0 2 0
(A−λ1 I)u 1 ⎛ 0 1 0 ⎞ u1= 0
⎝ -1 1 2 ⎠
0 2 0 2 2
⎛ 0 1 0⎞⎛⎞ 0 0 u1= ⎛ 0⎞
⎝ -1 1 2 ⎠⎝⎠ 1 ⎝ 1⎠
0 2 0
v′1 ⎛ 0 1 0 ⎞= 0 v′1(+21-1 0 )
⎝ -1 1 2 ⎠
perλ2
−a+ 2b= 0
-a + b + c = 0
-1 2 0
⎛ 0 0 0 ⎞ u 2= 0 a = 2b
⎝ −1 1 1 ⎠ -2b + b + c = 0
c=b
2
-1 2 0 u2= ⎛ 1⎞
v′2 ⎛ 0 0 0⎞
⎝ 1⎠
⎝ -1 1 1 ⎠ v′2( 0 1 0 )
perλ3
-2 2 0 0
⎛ 0−1 0 u⎞3= 0 u3= ⎛ 0⎞
⎝ -1 1 0⎠ ⎝ 1⎠
-2 2 0
v′3 ⎛ 0−1 0 = 0⎞ v3′ = (−20 1
1 )
⎝ -1 1 0⎠
274 8. Elements of Stability Theory
So:
Λ={−1; 0; 1}
ΛE{-1; 1}
Λ0 ={−1; 0}
ΛE0-1
A
=ẋ+
B
x+
Bs
us
⎧
⎪⎨ y=Cx
p=Cp x
⎪⎩
therefore the overall representation with the state is:
A
=ẋ+
B
x+
BC
us p x
&
y = Cx
8.6. The stability of interconnected systems 275
A
=
(ẋ+BC
s p x + Bu
&
y = Cx
Starting from the dynamic matrices of the subsystemsi that is starting from the
(AI , Bi , Ci ), the matrix A is first identified and then, starting from this,
the dynamic matrix of the interconnected system is (A+B s Cp ).
Considering that an interconnected system is characterized
A directed acyclic flow graph has as eigenvalues the aggregate of the eigenvalues.
of the subsystems that make it up, which is equivalent to saying that the property of
equivalence of stability to the stability property of subsystems, is understood
that the change in stability is connected to the transition from A to A (A+B s Cp ).
A
=ẋ+
Bu
xs
&
y = -Cp x
u=v-y
A
=ẋ+
B
xu
S: &
y = Cx
postou=−ysi obtains
A
=
(ẋ−BCx)+BV
S f: &
y = Cx
and it is put
The overall system will be stable when the second polynomial has no zeros.
the positive real part.
Observation 1.
dC H (s)
=|1 +F(s)|=|I+C(sI−A) −1 B|.
dAP(s)
|I−A+BC|=|sI−A||1 + (sI−A)|−1 BC
=|sI−A||1 +C(sI−A) −1 B|
dC H (s) =|sI−A+BC|
dAP(s) = |sI - A|
dCH(s) −1
= 1 + C(sI−A) B
dAP(s) ! "
it is understood that
sdiolshetnihncagseephsaofonsitaivrhaow
sngiofsoecprhetnieram
sesnetaousveiprhetnad
we must distinguish two situations: the case of a zero with a real part
positive and the case of a zero with a negative real part.
N=∆ϕ (1 + F)
278 8. Elements of stability theory
Note that the criterion allows for evaluating the stability of the overall system.
starting from the knowledge of the stability properties of the chain system
direct and from the graphic examination of the frequency behavior.
In the particular case where in the direct chain system there isp= 0,
if it has the criterion in its so-called reduced form
N= 0.
Some simple examples and the extension of the criterion to the case where
All eigenvalues have zero real part.
Consider the system described by the following transfer function
with gain, k, and time constants, τi , positive.
k
F(s) =
(1 + τ1 s)(1 +τ2 s)(1 +τ3 s)
The examination of the polar representation of F(jω) shown in the following figure
it allows to understand that the unit feedback system is stable if the point
8.6. The stability of interconnected systems 279
(-1, j0) +∞ 0
-∞
Figure 8.1
Seen in the complex plane, there will be a pole at zero and one at−1if one
τ builds
the image according to the polar diagram, as it approaches zero the curve
to infinity. But how do you calculate the number of turns around the point
Is the diagram open at (−1; j0)?
280 8. Elements of stability theory
-
0
(-1, j0) -∞
+∞
+
0
Figure 8.2
So the diagram needs to be closed and the result is different depending on the direction.
the diagram is closed. To overcome this ambiguity, one must proceed
as follows:
the polar diagram is no longer considered, but the Nyquist diagram, that is to say the
diagram showing how to construct the image according to the function F no longer
of the imaginary axis, but of the so-called Nyquist path: a path un-
it coincides with the imaginary axis at almost all points, but that around
At the singular points, the poles of the function are located on the imaginary axis, it leaves
of hooks that is to say it avoids them, and by convention it is established that this path
you leave the poles that have zero real part to the left, that is, you "consider" them
a negative real part. With this convention, even the eigenvalues have a part
Nothing will be considered in the negative complex half-plane. Note that
that, furthermore, with such an agreement,p= 0 does not correspond to a system in
a directly stable chain and not at the limit of stability.
induce conformal transformations, that is, those that preserve angles and directions
distance. So if in the domain to go from point A to point B, you
a hook traverses and a hook corresponds to a phase shift of π,
the same thing will happen in the image as well, so a rotation must be done
Moreover, the travel lines are maintained, which in this case leave
on the left the singular point and thus also the corresponding point, that is the
the image point must be left on the left. The image point is the point
improper that goes to infinity, this means that the closure of the diagram of
Nyquist must be done leaving the improper point on the left which is the image.
from the singular point.
Let's go back to the example; in this case, after the closure, it can be seen that the
Point (−1;j0) remains external, regardless of the gain value.
N=pp= 0
N=pp
the number of rotations the function F(jω) performs around the point
(−1;j0), when ω varies from −∞ to +∞, expresses the necessary and sufficient condition
asymptotic stability coefficient of the unity feedback system that has
in a direct chain.
We also saw: how the criterion allows for evaluation
stability as the gain of the direct chain system varies; in what way
It is possible to study stability in cases of non-unitary feedback.
282 8. Elements of stability theory
Example:
SystemS:
k(s + 1)
F(s) =
s(s−1)(s+10)
and, obviously
pp = 1
(-1, j0) -∞
+∞
-
0
Figure 8.3
8.6. The stability of interconnected systems 283
recalling the considerations made in the last lesson regarding the study
As the profit varies, it is understood that regarding stability
the unit feedback system will be unstable for gain values
lower than a critical value, k,∗ , stable for higher values.
How to evaluate∗ ?
For this purpose, the Routh criterion can be applied to the characteristic polynomial.
closed-loop system:
k(s + 1)
F(s) =
s(s−1)(s+10)
F(s)
W(s) =
1 + F(s)
k(s + 1)
W(s) =
s3+ 9s2+ (k−10)s+k
If a preliminary check is noted >10 it is necessary. The critical value may be
10 or any other potential greater value corresponding to the sign change
of some coefficient of the first column of the table. It is:
1 k−10
9 k
8k−90
9
k
90
k>0 k>
8
what it implies
90
k∗ =
8
The unit feedback system for this value is at the limit of stability.
s+ 1
F(s) = unstable system in direct chain p p= 0
s2
284 8. Elements of stability theory
s+1
F(s) =
s 2+ 1
Applying the Nyquist criterion, it will be checked that for every positive value,
both unit feedback systems are asymptotically stable.
The following considerations are useful for deciding on stability in the case where
there are eigenvalues with zero real partAPe/o indCH. Sia jω0 a point
of the imaginary axis and indicate withAPedmCHthe multiplicity with which
that point and possible zero of the aforementioned polynomials. It is also assumed to be known
mAP. If 1 + F(jω0 ) =c̸= 0,mCH=mAPit will be simply established if
mAP= 1, instability itselfAP1. If F(jω0 ) =−1,m CHmAPand if it has
instability weAP≥0. If, infinite, F(jω) = ∞, mAP> mCHand it has been established
if the closure involves a rotation dimAPπ.
So far it has been seen how to decide on stability starting from the dia-
Nyquist diagram, but this diagram has been deduced from Bode's.
So some information about the stability of the feedback system is already
contained in the Bode diagram of the direct chain system. In fact, it is pos-
you can evaluate the stability conditions on the Bode diagrams. This will
we will refer to a direct chain system free of poles except for
positive real.
With reference to the case where the feedback system is stable, and therefore
that the point (−1;j0) is to the left of the crossing; it is possible to introduce
two parameters that account for the stability of the system.
One of them is the module of the F when the phase is -π; if it is less than one,
Negative dB, the system will be stable. The smaller it is than one, the greater.
it goes to stability. It is customary to give this parameter the value of the gain when
the phase e−π, the name margin of gain as it tells us how much we
can modify the system's gain to reach a situation that is
stability limit.
The other parameter that is perhaps more common to refer to is the so-called
phase margin corresponds to the phase variation necessary to lead to
stability limit.
mϕ =π−phaseF |F|=1
8.7. The stability of discrete-time systems 285
It is understood that the definitions remain unchanged, the study conducted so far
It allows us to understand that the condition of asymptotic stability of systems
discrete-time linear and that the eigenvalues have a modulus less than one.
The verification can be conducted using different criteria; among these, those here
below exposed.
Stability conditions for discrete-time systems.
xand=Axe
(A−I)xe = 0
286 8. Elements of stability theory
||φ(t)|| < k
e
lim∥φ(t)∥= 0
t approaches infinity
φ(t) = At
in this case the limitation and the convergence to zero of the norm impose
that the eigenvalues belong to the unit circle.
So for the eigenvalues that are outside, there are divergent motions.
for those inside, there are converging motions, on the boundary there are modes
which can have constant amplitude if the multiplicity is one, if instead it is
greater than one, divergent motions occur again. The stability condition
Then it will be that the eigenvalues are all inside the circle, at most.
Simple eigenvalues are allowed at the border. For asymptotic stability
they must be strictly inside. In conclusion, even in this case
the stability condition corresponds to the allocation of the eigenvalues in a
fixed region of the complex plane. simple stability
1
& |λi|≤1
|λi1>|<1
asymptotic stability
|λI |<1
where it is indicated with the apostrophe 1 or >1 the unitary multiplicity or greater than one
of the self-value.
A first criterion is represented by Routh's criterion itself.
With reference to the polynomial
variable substitution
s+ 1
z←
s−1
it corresponds to performing a transformation from the unit circle to the half-plane
negative of the complex plane. On the polynomial it is possible to apply the
Routh's criterion and obtaining information on stability; in fact, the zeros in
Zeros with real part greater than zero correspond to zcon modulo>1.
A second criterion is the so-called Jury criterion, which consists, just like
Routh's, in verifying that among the coefficients of a table constructed at
This requires the verification of precise conditions.
a0 a1 ... ... an
⎛ an an−1. . . . . . a0 ⎞
b0 b1 . . . . . . bn−1
⎜ ⎟
⎜ bn−1 ... ... b0 ⎟
⎜ ⎟
⎜ c0 . . . . . . c n−2 ⎟
⎜ ⎟
⎝ cn−2 ... ... c0 ⎠
..
.
t 0t 1t 2
a0 a1
b0← ' anan−1 '
' '
' '
' '
aa
bn−1 ← ' 0 n '
' n a0
a '
' '
in general ' '
a0ak+1
bk← '' anan−k−1 '
'
' '
' '
The statement of Jury's criterion is as follows: necessary conditions and
sufficient for asymptotic stability are
288 8. Elements of stability theory
d(1) > 0
-1n d(-1) >0
|an |>|a 0 |
|b0 |>|bn−1 |
..
.
|t0 |>|t2 |
Example:
z 2+z+ 0.5
The first three conditions are preliminarily verified that do not depend on
coefficients of the table; thus, the table is constructed characterized by the
single line
1 1 0.5
e d e fi n it a n e g a t iv a in
e Ix.
The criterion is therefore stated in the same way and this time one must find
a function whose first derivative is defined (negative semidefinite). An-
once again the criterion of life under conditions that are only sufficient, in the
I feel that if in correspondence with V(x) defined as positive, the ∆V(x) is not the case.
Negative definite (semidefinite) cannot decide the stability of the state.
of balance.
Example:
x2 (t)
x (t + 1) =
⎧ 1 1 + x2 (t)2
⎪⎪
⎨ x1 (t)
x2 (t+ 1) =
⎪⎪ 1 +x2 (t)2
⎩
the state of equilibriume= 0 and posted:
V(x) = x12+x22
1
∆V(x) = $ % ≤0
-1V(x)
(1 + x22)2
This last quantity is negative semi-definite and therefore the equilibrium state.
qeuiebaltsyloblaglasi tiytilaerni ,etatsebaltsasi ti
the stated conditions of unlimited validity of the conditions apply throughout
the space.
Also with reference to discrete-time systems in the linear case, the cri-
terio is specified by creating necessary and sufficient conditions that are easy
verify.
290 8. Elements of stability theory
discrete time
A′ QA−Q=−P
V(x) =x′ Qx
For the necessity, one proceeds in a formally similar manner to what has been seen.
in the case of continuous time showing that
∞
Q= A′k PA k
#0
The comparison between the Routh criterion and the Lyapunov criterion is interesting in
due to the complexity of calculations in terms of the number of operations. Obvious-
Unfortunately, everything depends on the structure of the dynamic matrix; in the presence of
A high number of null elements in the matrix can be advantageous.
apply the Routh criterion as it should not be forgotten that for ap-
To apply the Routh criterion, it is necessary to first calculate the polynomial.
characteristic.
You have seen in this paragraph how the Lyapunov criterion is applied.
larizzi in the linear context, continuous time and discrete time, giving life to
necessary and sufficient conditions. This does not happen in the non-linear context where
the possibility of achieving a positive verification of stability is linked to
calculation of a function that satisfies the known conditions. This aspect makes
the application of the method is difficult in cases where the verification has not occurred
positive outcome using the defined positive elementary quadratic functions.
What to do when you cannot directly study stability
of a nonlinear system applying Lyapunov's criterion? Are there others?
methods?
In reality, another method, moreover with a complementary role to that
carried out from the Lyapunov criterion, it is based on the study of stability
of the nonlinear system based on the linear approximation. This will be
the subject of the next lesson.
292 8. Elements of stability theory
xe = √ a
⇓
0 < a < 1 if and only if |1 - 2a|√< 1
From this analysis it is derived that for 0 < a < 1 the solution xeto be √
certainly calculated with the given algorithm.
8.7. The stability of discrete-time systems 293
The obtained result has a local validity and the maximum remains undefined.
Deviation allowed from the initial condition with respect to the true value (that one
of equilibrium). More precise information about any deviation,
stability region, can be obtained using the Lyapunov criterion with
V(x) = |a - x|2 |
The calculation of V(x(t+ 1)) and ∆(V(x)) allows for easy verification.
if we assume 0 < a < 1 that
0≤x≤1
defines the region of stability.
Two examples to be carried out
Leslie model
price dynamics
The second example generalizes the example of price dynamics in...
Equilibrium conditions in the presence of a single good. Let us consider a market
characterized by the presence of two assets
d(t+ 1) = Ap(t+ 1) + d0
&
o(t+ 1) = Ep(t) + o0
doveA,E,d0eo0they are matrices and vectors of dimension two whose coefficients
they allow for simple interpretations in terms of production and demand.
equilibrium conditions it is possible to establish a dynamics equation of
price according to a standard report.