Controllability

From HandWiki
Short description: Dynamic system property

Controllability is an important property of a control system and plays a crucial role in many control problems, such as stabilization of unstable systems by feedback, or optimal control.

Controllability and observability are dual aspects of the same problem.

Roughly, the concept of controllability denotes the ability to move a system around in its entire configuration space using only certain admissible manipulations. The exact definition varies slightly within the framework or the type of models applied.

The following are examples of variations of controllability notions which have been introduced in the systems and control literature:

  • State controllability
  • Output controllability
  • Controllability in the behavioural framework

State controllability

The state of a deterministic system, which is the set of values of all the system's state variables (those variables characterized by dynamic equations), completely describes the system at any given time. In particular, no information on the past of a system is needed to help in predicting the future, if the states at the present time are known and all current and future values of the control variables (those whose values can be chosen) are known.

Complete state controllability (or simply controllability if no other context is given) describes the ability of an external input (the vector of control variables) to move the internal state of a system from any initial state to any final state in a finite time interval.[1]:737

That is, we can informally define controllability as follows: If for any initial state [math]\displaystyle{ \mathbf{x_0} }[/math] and any final state [math]\displaystyle{ \mathbf{x_f} }[/math] there exists an input sequence to transfer the system state from [math]\displaystyle{ \mathbf{x_0} }[/math] to [math]\displaystyle{ \mathbf{x_f} }[/math] in a finite time interval, then the system modeled by the state-space representation is controllable. For the simplest example of a continuous, LTI system, the row dimension of the state space expression [math]\displaystyle{ \dot{\mathbf{x}} = \mathbf{A}\mathbf{x}(t) + \mathbf{B}\mathbf{u}(t) }[/math] determines the interval; each row contributes a vector in the state space of the system. If there are not enough such vectors to span the state space of [math]\displaystyle{ \mathbf{x} }[/math], then the system cannot achieve controllability. It may be necessary to modify [math]\displaystyle{ \mathbf{A} }[/math] and [math]\displaystyle{ \mathbf{B} }[/math] to better approximate the underlying differential relationships it estimates to achieve controllability.

Controllability does not mean that a reached state can be maintained, merely that any state can be reached.

Controllability does not mean that arbitrary paths can be made through state space, only that there exists a path within the prescribed finite time interval.

Continuous linear systems

Consider the continuous linear system [note 1]

[math]\displaystyle{ \dot{\mathbf{x}}(t) = A(t) \mathbf{x}(t) + B(t) \mathbf{u}(t) }[/math]
[math]\displaystyle{ \mathbf{y}(t) = C(t) \mathbf{x}(t) + D(t) \mathbf{u}(t). }[/math]

There exists a control [math]\displaystyle{ u }[/math] from state [math]\displaystyle{ x_0 }[/math] at time [math]\displaystyle{ t_0 }[/math] to state [math]\displaystyle{ x_1 }[/math] at time [math]\displaystyle{ t_1 \gt t_0 }[/math] if and only if [math]\displaystyle{ x_1 - \phi(t_0,t_1)x_0 }[/math] is in the column space of

[math]\displaystyle{ W(t_0,t_1) = \int_{t_0}^{t_1} \phi(t_0,t)B(t)B(t)^{T}\phi(t_0,t)^{T} dt }[/math]

where [math]\displaystyle{ \phi }[/math] is the state-transition matrix, and [math]\displaystyle{ W(t_0,t_1) }[/math] is the Controllability Gramian.

In fact, if [math]\displaystyle{ \eta_0 }[/math] is a solution to [math]\displaystyle{ W(t_0,t_1)\eta = x_1 - \phi(t_0,t_1)x_0 }[/math] then a control given by [math]\displaystyle{ u(t) = -B(t)^{T}\phi(t_0,t)^{T}\eta_0 }[/math] would make the desired transfer.

Note that the matrix [math]\displaystyle{ W }[/math] defined as above has the following properties:

  • [math]\displaystyle{ W(t_0,t_1) }[/math] is symmetric
  • [math]\displaystyle{ W(t_0,t_1) }[/math] is positive semidefinite for [math]\displaystyle{ t_1 \geq t_0 }[/math]
  • [math]\displaystyle{ W(t_0,t_1) }[/math] satisfies the linear matrix differential equation
[math]\displaystyle{ \frac{d}{dt}W(t,t_1) = A(t)W(t,t_1)+W(t,t_1)A(t)^{T}-B(t)B(t)^{T}, \; W(t_1,t_1) = 0 }[/math]
  • [math]\displaystyle{ W(t_0,t_1) }[/math] satisfies the equation
[math]\displaystyle{ W(t_0,t_1) = W(t_0,t) + \phi(t_0,t)W(t,t_1)\phi(t_0,t)^{T} }[/math][2]

Rank condition for controllability

The Controllability Gramian involves integration of the state-transition matrix of a system. A simpler condition for controllability is a rank condition analogous to the Kalman rank condition for time-invariant systems.

Consider a continuous-time linear system [math]\displaystyle{ \Sigma }[/math] smoothly varying in an interval [math]\displaystyle{ [t_0,t] }[/math] of [math]\displaystyle{ \mathbb{R} }[/math]:

[math]\displaystyle{ \dot{\mathbf{x}}(t) = A(t) \mathbf{x}(t) + B(t) \mathbf{u}(t) }[/math]
[math]\displaystyle{ \mathbf{y}(t) = C(t) \mathbf{x}(t) + D(t) \mathbf{u}(t). }[/math]

The state-transition matrix [math]\displaystyle{ \phi }[/math] is also smooth. Introduce the n x m matrix-valued function [math]\displaystyle{ M_0(t) = \phi(t_0,t)B(t) }[/math] and define

[math]\displaystyle{ M_k(t) }[/math] = [math]\displaystyle{ \frac{\mathrm{d^k} M_0}{\mathrm{d} t^k}(t), k\geqslant 1 }[/math].

Consider the matrix of matrix-valued functions obtained by listing all the columns of the [math]\displaystyle{ M_i }[/math], [math]\displaystyle{ i = 0,1, \ldots, k }[/math]:

[math]\displaystyle{ M^{(k)}(t) := \left[M_0(t), \ldots, M_k(t)\right] }[/math].

If there exists a [math]\displaystyle{ \bar{t} \in [t_0,t] }[/math] and a nonnegative integer k such that [math]\displaystyle{ \operatorname{rank}M^{(k)}(\bar{t})=n }[/math], then [math]\displaystyle{ \Sigma }[/math] is controllable.[3]

If [math]\displaystyle{ \Sigma }[/math] is also analytically varying in an interval [math]\displaystyle{ [t_0,t] }[/math], then [math]\displaystyle{ \Sigma }[/math] is controllable on every nontrivial subinterval of [math]\displaystyle{ [t_0,t] }[/math] if and only if there exists a [math]\displaystyle{ \bar{t} \in [t_0,t] }[/math] and a nonnegative integer k such that [math]\displaystyle{ \operatorname{rank}M^{(k)}(t_i)=n }[/math].[3]

The above methods can still be complex to check, since it involves the computation of the state-transition matrix [math]\displaystyle{ \phi }[/math]. Another equivalent condition is defined as follow. Let [math]\displaystyle{ B_0(t) = B(t) }[/math], and for each [math]\displaystyle{ i \geq 0 }[/math], define

[math]\displaystyle{ B_{i+1}(t) }[/math]= [math]\displaystyle{ A(t)B_i(t) - \frac{\mathrm{d}}{\mathrm{d} t}B_i(t). }[/math]

In this case, each [math]\displaystyle{ B_i }[/math] is obtained directly from the data [math]\displaystyle{ (A(t),B(t)). }[/math] The system is controllable if there exists a [math]\displaystyle{ \bar{t} \in [t_0,t] }[/math] and a nonnegative integer [math]\displaystyle{ k }[/math] such that [math]\displaystyle{ \textrm{rank}( \left[ B_0(\bar{t}), B_1(\bar{t}), \ldots, B_k(\bar{t}) \right]) = n }[/math].[3]

Example

Consider a system varying analytically in [math]\displaystyle{ (-\infty,\infty) }[/math] and matrices

[math]\displaystyle{ A(t) = \begin{bmatrix} t & 1 & 0\\ 0 & t^{3} & 0\\ 0 & 0 & t^{2} \end{bmatrix} }[/math], [math]\displaystyle{ B(t) = \begin{bmatrix} 0 \\ 1 \\ 1 \end{bmatrix}. }[/math] Then [math]\displaystyle{ [B_0(0),B_1(0),B_2(0),B_3(0)] = \begin{bmatrix} 0 & 1 & 0 &-1\\ 1 & 0 & 0&0\\ 1 & 0 & 0&2 \end{bmatrix} }[/math] and since this matrix has rank 3, the system is controllable on every nontrivial interval of [math]\displaystyle{ \mathbb{R} }[/math].

Continuous linear time-invariant (LTI) systems

Consider the continuous linear time-invariant system

[math]\displaystyle{ \dot{\mathbf{x}}(t) = A \mathbf{x}(t) + B \mathbf{u}(t) }[/math]
[math]\displaystyle{ \mathbf{y}(t) = C \mathbf{x}(t) + D \mathbf{u}(t) }[/math]

where

[math]\displaystyle{ \mathbf{x} }[/math] is the [math]\displaystyle{ n \times 1 }[/math] "state vector",
[math]\displaystyle{ \mathbf{y} }[/math] is the [math]\displaystyle{ m \times 1 }[/math] "output vector",
[math]\displaystyle{ \mathbf{u} }[/math] is the [math]\displaystyle{ r \times 1 }[/math] "input (or control) vector",
[math]\displaystyle{ A }[/math] is the [math]\displaystyle{ n \times n }[/math] "state matrix",
[math]\displaystyle{ B }[/math] is the [math]\displaystyle{ n \times r }[/math] "input matrix",
[math]\displaystyle{ C }[/math] is the [math]\displaystyle{ m \times n }[/math] "output matrix",
[math]\displaystyle{ D }[/math] is the [math]\displaystyle{ m \times r }[/math] "feedthrough (or feedforward) matrix".

The [math]\displaystyle{ n \times nr }[/math] controllability matrix is given by

[math]\displaystyle{ R = \begin{bmatrix}B & AB & A^{2}B & ...& A^{n-1}B\end{bmatrix} }[/math]

The system is controllable if the controllability matrix has full row rank (i.e. [math]\displaystyle{ \operatorname{rank}(R)=n }[/math]).

Discrete linear time-invariant (LTI) systems

For a discrete-time linear state-space system (i.e. time variable [math]\displaystyle{ k\in\mathbb{Z} }[/math]) the state equation is

[math]\displaystyle{ \textbf{x}(k+1) = A\textbf{x}(k) + B\textbf{u}(k) }[/math]

where [math]\displaystyle{ A }[/math] is an [math]\displaystyle{ n \times n }[/math] matrix and [math]\displaystyle{ B }[/math] is a [math]\displaystyle{ n \times r }[/math] matrix (i.e. [math]\displaystyle{ \mathbf{u} }[/math] is [math]\displaystyle{ r }[/math] inputs collected in a [math]\displaystyle{ r \times 1 }[/math] vector). The test for controllability is that the [math]\displaystyle{ n \times nr }[/math] matrix

[math]\displaystyle{ \mathcal{C} = \begin{bmatrix}B & AB & A^{2}B & \cdots & A^{n-1}B\end{bmatrix} }[/math]

has full row rank (i.e., [math]\displaystyle{ \operatorname{rank}(\mathcal C) = n }[/math]). That is, if the system is controllable, [math]\displaystyle{ \mathcal C }[/math] will have [math]\displaystyle{ n }[/math] columns that are linearly independent; if [math]\displaystyle{ n }[/math] columns of [math]\displaystyle{ \mathcal C }[/math] are linearly independent, each of the [math]\displaystyle{ n }[/math] states is reachable by giving the system proper inputs through the variable [math]\displaystyle{ u(k) }[/math].

Derivation

Given the state [math]\displaystyle{ \textbf{x}(0) }[/math] at an initial time, arbitrarily denoted as k=0, the state equation gives [math]\displaystyle{ \textbf{x}(1) = A\textbf{x}(0) + B\textbf{u}(0), }[/math] then [math]\displaystyle{ \textbf{x}(2) = A\textbf{x}(1) + B\textbf{u}(1)= A^2\textbf{x}(0)+AB\textbf{u}(0)+B\textbf{u}(1), }[/math] and so on with repeated back-substitutions of the state variable, eventually yielding

[math]\displaystyle{ \textbf{x}(n)=B\textbf{u}(n-1) + AB\textbf{u}(n-2) + \cdots + A^{n-1}B\textbf{u}(0) + A^n\textbf{x}(0) }[/math]

or equivalently

[math]\displaystyle{ \textbf{x}(n)-A^n\textbf{x}(0)= [B \, \, AB \, \, \cdots \, \, A^{n-1}B] [\textbf{u}^T(n-1) \, \, \textbf{u}^T(n-2) \, \, \cdots \, \, \textbf{u}^T(0)]^T. }[/math]

Imposing any desired value of the state vector [math]\displaystyle{ \textbf{x}(n) }[/math] on the left side, this can always be solved for the stacked vector of control vectors if and only if the matrix of matrices at the beginning of the right side has full row rank.

Example

For example, consider the case when [math]\displaystyle{ n=2 }[/math] and [math]\displaystyle{ r=1 }[/math] (i.e. only one control input). Thus, [math]\displaystyle{ B }[/math] and [math]\displaystyle{ A B }[/math] are [math]\displaystyle{ 2 \times 1 }[/math] vectors. If [math]\displaystyle{ \begin{bmatrix}B & AB\end{bmatrix} }[/math] has rank 2 (full rank), and so [math]\displaystyle{ B }[/math] and [math]\displaystyle{ AB }[/math] are linearly independent and span the entire plane. If the rank is 1, then [math]\displaystyle{ B }[/math] and [math]\displaystyle{ AB }[/math] are collinear and do not span the plane.

Assume that the initial state is zero.

At time [math]\displaystyle{ k=0 }[/math]: [math]\displaystyle{ x(1) = A\textbf{x}(0) + B\textbf{u}(0) = B\textbf{u}(0) }[/math]

At time [math]\displaystyle{ k=1 }[/math]: [math]\displaystyle{ x(2) = A\textbf{x}(1) + B\textbf{u}(1) = AB\textbf{u}(0) + B\textbf{u}(1) }[/math]

At time [math]\displaystyle{ k=0 }[/math] all of the reachable states are on the line formed by the vector [math]\displaystyle{ B }[/math]. At time [math]\displaystyle{ k=1 }[/math] all of the reachable states are linear combinations of [math]\displaystyle{ AB }[/math] and [math]\displaystyle{ B }[/math]. If the system is controllable then these two vectors can span the entire plane and can be done so for time [math]\displaystyle{ k=2 }[/math]. The assumption made that the initial state is zero is merely for convenience. Clearly if all states can be reached from the origin then any state can be reached from another state (merely a shift in coordinates).

This example holds for all positive [math]\displaystyle{ n }[/math], but the case of [math]\displaystyle{ n=2 }[/math] is easier to visualize.

Analogy for example of n = 2

Consider an analogy to the previous example system. You are sitting in your car on an infinite, flat plane and facing north. The goal is to reach any point in the plane by driving a distance in a straight line, come to a full stop, turn, and driving another distance, again, in a straight line. If your car has no steering then you can only drive straight, which means you can only drive on a line (in this case the north-south line since you started facing north). The lack of steering case would be analogous to when the rank of [math]\displaystyle{ C }[/math] is 1 (the two distances you drove are on the same line).

Now, if your car did have steering then you could easily drive to any point in the plane and this would be the analogous case to when the rank of [math]\displaystyle{ C }[/math] is 2.

If you change this example to [math]\displaystyle{ n=3 }[/math] then the analogy would be flying in space to reach any position in 3D space (ignoring the orientation of the aircraft). You are allowed to:

  • fly in a straight line
  • turn left or right by any amount (Yaw)
  • direct the plane upwards or downwards by any amount (Pitch)

Although the 3-dimensional case is harder to visualize, the concept of controllability is still analogous.

Nonlinear systems

Nonlinear systems in the control-affine form

[math]\displaystyle{ \dot{\mathbf{x}} = \mathbf{f(x)} + \sum_{i=1}^m \mathbf{g}_i(\mathbf{x})u_i }[/math]

are locally accessible about [math]\displaystyle{ x_0 }[/math] if the accessibility distribution [math]\displaystyle{ R }[/math] spans [math]\displaystyle{ n }[/math] space, when [math]\displaystyle{ n }[/math] equals the rank of [math]\displaystyle{ x }[/math] and R is given by:[4]

[math]\displaystyle{ R = \begin{bmatrix} \mathbf{g}_1 & \cdots & \mathbf{g}_m & [\mathrm{ad}^k_{\mathbf{g}_i}\mathbf{\mathbf{g}_j}] & \cdots & [\mathrm{ad}^k_{\mathbf{f}}\mathbf{\mathbf{g}_i}] \end{bmatrix}. }[/math]

Here, [math]\displaystyle{ [\mathrm{ad}^k_{\mathbf{f}}\mathbf{\mathbf{g}}] }[/math] is the repeated Lie bracket operation defined by

[math]\displaystyle{ [\mathrm{ad}^k_{\mathbf{f}}\mathbf{\mathbf{g}}] = \begin{bmatrix} \mathbf{f} & \cdots & j & \cdots & \mathbf{[\mathbf{f}, \mathbf{g}]} \end{bmatrix}. }[/math]

The controllability matrix for linear systems in the previous section can in fact be derived from this equation.

Null Controllability

If a discrete control system is null-controllable, it means that there exists a controllable [math]\displaystyle{ u(k) }[/math] so that [math]\displaystyle{ x(k_0) = 0 }[/math] for some initial state [math]\displaystyle{ x(0) = x_0 }[/math]. In other words, it is equivalent to the condition that there exists a matrix [math]\displaystyle{ F }[/math] such that [math]\displaystyle{ A+BF }[/math] is nilpotent.

This can be easily shown by controllable-uncontrollable decomposition.

Output controllability

Output controllability is the related notion for the output of the system (denoted y in the previous equations); the output controllability describes the ability of an external input to move the output from any initial condition to any final condition in a finite time interval. It is not necessary that there is any relationship between state controllability and output controllability. In particular:

  • A controllable system is not necessarily output controllable. For example, if matrix D = 0 and matrix C does not have full row rank, then some positions of the output are masked by the limiting structure of the output matrix, and therefore unachievable. Moreover, even though the system can be moved to any state in finite time, there may be some outputs that are inaccessible by all states. A trivial numerical example uses D=0 and a C matrix with at least one row of zeros; thus, the system is not able to produce a non-zero output along that dimension.
  • An output controllable system is not necessarily state controllable. For example, if the dimension of the state space is greater than the dimension of the output, then there will be a set of possible state configurations for each individual output. That is, the system can have significant zero dynamics, which are trajectories of the system that are not observable from the output. Consequently, being able to drive an output to a particular position in finite time says nothing about the state configuration of the system.

For a linear continuous-time system, like the example above, described by matrices [math]\displaystyle{ A }[/math], [math]\displaystyle{ B }[/math], [math]\displaystyle{ C }[/math], and [math]\displaystyle{ D }[/math], the [math]\displaystyle{ m \times (n+1)r }[/math] output controllability matrix

[math]\displaystyle{ \begin{bmatrix} CB & CAB & CA^2B & \cdots & CA^{n-1}B & D\end{bmatrix} }[/math]

has full row rank (i.e. rank [math]\displaystyle{ m }[/math]) if and only if the system is output controllable.[1]:742

Controllability under input constraints

In systems with limited control authority, it is often no longer possible to move any initial state to any final state inside the controllable subspace. This phenomenon is caused by constraints on the input that could be inherent to the system (e.g. due to saturating actuator) or imposed on the system for other reasons (e.g. due to safety-related concerns). The controllability of systems with input and state constraints is studied in the context of reachability[5] and viability theory.[6]

Controllability in the behavioral framework

In the so-called behavioral system theoretic approach due to Willems (see people in systems and control), models considered do not directly define an input–output structure. In this framework systems are described by admissible trajectories of a collection of variables, some of which might be interpreted as inputs or outputs.

A system is then defined to be controllable in this setting, if any past part of a behavior (trajectory of the external variables) can be concatenated with any future trajectory of the behavior in such a way that the concatenation is contained in the behavior, i.e. is part of the admissible system behavior.[7]:151

Stabilizability

A slightly weaker notion than controllability is that of stabilizability. A system is said to be stabilizable when all uncontrollable state variables can be made to have stable dynamics. Thus, even though some of the state variables cannot be controlled (as determined by the controllability test above) all the state variables will still remain bounded during the system's behavior.[8]

Reachable set

Let T ∈ Т and x ∈ X (where X is the set of all possible states and Т is an interval of time). The reachable set from x in time T is defined as:[3]

[math]\displaystyle{ R^T{(x)} = \left\{ z \in X : x \overset{T}{\rightarrow} z \right\} }[/math], where xTz denotes that there exists a state transition from x to z in time T.

For autonomous systems the reachable set is given by :

[math]\displaystyle{ \mathrm{Im}(R)=\mathrm{Im}(B)+\mathrm{Im}(AB)+....+\mathrm{Im}(A^{n-1}B) }[/math],

where R is the controllability matrix.

In terms of the reachable set, the system is controllable if and only if [math]\displaystyle{ \mathrm{Im}(R)=\mathbb{R}^n }[/math].

Proof We have the following equalities:

[math]\displaystyle{ R=[B\ AB ....A^{n-1}B] }[/math]
[math]\displaystyle{ \mathrm{Im}(R)=\mathrm{Im}([B\ AB ....A^{n-1}B]) }[/math]
[math]\displaystyle{ \mathrm{dim(Im}(R))=\mathrm{rank}(R) }[/math]

Considering that the system is controllable, the columns of R should be linearly independent. So:

[math]\displaystyle{ \mathrm{dim(Im}(R))=n }[/math]
[math]\displaystyle{ \mathrm{rank}(R)=n }[/math]
[math]\displaystyle{ \mathrm{Im}(R)=\R^{n}\quad \blacksquare }[/math]

A related set to the reachable set is the controllable set, defined by:

[math]\displaystyle{ C^T{(x)} = \left\{ z \in X : z \overset{T}{\rightarrow} x \right\} }[/math].

The relation between reachability and controllability is presented by Sontag:[3]

(a) An n-dimensional discrete linear system is controllable if and only if:

[math]\displaystyle{ R(0)=R^k{(0)=X} }[/math] (Where X is the set of all possible values or states of x and k is the time step).

(b) A continuous-time linear system is controllable if and only if:

[math]\displaystyle{ R(0)=R^e{(0)=X} }[/math] for all e>0.

if and only if [math]\displaystyle{ C(0)=C^e{(0)=X} }[/math] for all e>0.

Example Let the system be an n dimensional discrete-time-invariant system from the formula:

Φ(n,0,0,w)=[math]\displaystyle{ \sum\limits_{i=1}^n A^{i-1}Bw(n-1) }[/math] (Where Φ(final time, initial time, state variable, restrictions) is defined is the transition matrix of a state variable x from a initial time 0 to a final time n with some restrictions w).

It follows that the future state is in [math]\displaystyle{ R^k{(0)} }[/math] ⇔ it is in the image of the linear map:

Im(R)=R(A,B)≜ Im([math]\displaystyle{ [B\ AB ....A^{n-1}B] }[/math]),

which maps,

[math]\displaystyle{ u^{n} }[/math]→X

When [math]\displaystyle{ u=K^{m} }[/math] and [math]\displaystyle{ X=K^{n} }[/math] we identify R(A,B) with a n by nm matrix whose columns are the columns of [math]\displaystyle{ B,\ AB, ....,A^{n-1}B }[/math] in that order. If the system is controllable the rank of [math]\displaystyle{ [B\ AB ....A^{n-1}B] }[/math] is n. If this is truth, the image of the linear map R is all of X. Based on that, we have:

[math]\displaystyle{ R(0)=R^k{(0)=X} }[/math] with XЄ[math]\displaystyle{ \R^{n} }[/math].

See also

Notes

  1. A linear time-invariant system behaves the same but with the coefficients being constant in time.

References

  1. 1.0 1.1 Katsuhiko Ogata (1997). Modern Control Engineering (3rd ed.). Upper Saddle River, NJ: Prentice-Hall. ISBN 978-0-13-227307-7. 
  2. Brockett, Roger W. (1970). Finite Dimensional Linear Systems. John Wiley & Sons. ISBN 978-0-471-10585-5. 
  3. 3.0 3.1 3.2 3.3 3.4 Eduardo D. Sontag, Mathematical Control Theory: Deterministic Finite Dimensional Systems.
  4. Isidori, Alberto (1989). Nonlinear Control Systems, p. 92–3. Springer-Verlag, London. ISBN:3-540-19916-0.
  5. Claire J. Tomlin; Ian Mitchell; Alexandre M. Bayen; Meeko Oishi (2003). "Computational Techniques for the Verification of Hybrid Systems". Proceedings of the IEEE 91 (7): 986–1001. doi:10.1109/jproc.2003.814621. http://www.cs.ubc.ca/~mitchell/Papers/publishedIEEEproc03.pdf. Retrieved 2012-03-04. 
  6. Jean-Pierre Aubin (1991). Viability Theory. Birkhauser. ISBN 978-0-8176-3571-8. 
  7. Jan Polderman; Jan Willems (1998). Introduction to Mathematical Systems Theory: A Behavioral Approach (1st ed.). New York: Springer Verlag. ISBN 978-0-387-98266-3. 
  8. Brian D.O. Anderson; John B. Moore (1990). Optimal Control: Linear Quadratic Methods. Englewood Cliffs, NJ: Prentice Hall. ISBN 978-0-13-638560-8. 

External links

fr:Représentation d'état#Systèmes linéaires