View source for Controllability

'''Controllability''' is an important property of a [[Control system|control system]] and plays a crucial role in many regulation problems, such as the stabilization of [[BIBO stability|unstable systems]] using feedback, tracking problems, obtaining [[Optimal control|optimal control]] strategies, or, simply prescribing an input that has a desired effect on the state.

Controllability and [[Observability|observability]] are [[Duality (mathematics)|dual]] notions. Controllability pertains to regulating the state by a choice of a suitable input, while observability pertains to being able to know the state by observing the output (assuming that the input is also being observed).

Broadly speaking, the concept of controllability relates to the ability to steer a system around in its configuration space using only certain admissible manipulations. The exact definition varies depending on the framework or the type of models dealt with.

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

* State controllability: the ability to steer the system between states
* Strong controllability: the ability to steer between states over any specified time window
* Collective controllability: the ability to simultaneously steer a collection of dynamical systems
* Trajectory controllability: the ability to steer along a predefined trajectory rather than just to a desired final state
* Output controllability: the ability to steer to specified values of the output
* Controllability in the behavioural framework: a compatibility condition between past and future input and output trajectories

== State controllability ==

The state of a [[Deterministic system|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.

''{{dfn|Complete state controllability}}'' (or simply ''{{dfn|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.<ref name="Ogata97">{{cite book |last= Ogata |first= Katsuhiko |date= 1997 |title= Modern Control Engineering |edition= 3rd |publisher= Prentice-Hall |location= Upper Saddle River, New Jersey |isbn= 978-0-13-227307-7}}</ref>{{rp|737}}

That is, we can informally define controllability as follows:
If for any initial state {{math|'''x'''{{sub|0}}}} and any final state {{math|'''x'''{{sub|f}}}} there exists an input sequence to transfer the system state from {{math|'''x'''{{sub|0}}}} to {{math|'''x'''{{sub|f}}}} in a finite time interval, then the system modeled by the [[State-space representation|state-space representation]] is controllable. For the simplest example of a continuous, linear time-invariant (LTI) system, the row dimension of the state space expression {{math|1='''ẋ''' = ''A'''''x'''(''t'') + ''B'''''u'''(''t'')}} 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|'''x'''}}, then the system cannot achieve controllability. It may be necessary to modify {{mvar|A}} and {{mvar|B}} 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 a finite time interval. When the time interval can also be specified, the [[Dynamical system|dynamical system]] is often referred to as being strongly controllable.

== Continuous linear systems ==

Consider the continuous linear system <ref group="note">A [[Linear time-invariant system|linear time-invariant system]] behaves the same but with the coefficients being constant in time.</ref>
<math display=block>\begin{align}
\dot{\mathbf{x}}(t) &= A(t) \mathbf{x}(t) + B(t) \mathbf{u}(t) \\
\mathbf{y}(t) &= C(t) \mathbf{x}(t) + D(t) \mathbf{u}(t).
\end{align}</math>

There exists a control {{mvar|u}} from state {{math|'''x'''{{sub|0}}}} at time {{math|''t''{{sub|0}}}} to state {{math|'''x'''{{sub|1}}}} at time {{math|''t''{{sub|1}} > ''t''{{sub|0}}}} if and only if {{math|'''x'''{{sub|1}} − ''ϕ''(''t''{{sub|0}}, ''t''{{sub|1}})'''x'''{{sub|0}}}} is in the [[Column space|column space]] of
<math display=block>W(t_0,t_1) = \int_{t_0}^{t_1} \phi(t_0,t)B(t)B(t)^\mathsf{T}\phi(t_0,t)^\mathsf{T} \mathrm{d}t</math>
where {{math|''ϕ''(''t{{sub|0}}'', {{itco|''t''}})}} is the [[State-transition matrix|state-transition matrix]], and {{math|''W''(''t''{{sub|0}}, ''t''{{sub|1}})}} is the [[Controllability Gramian]].

In fact, if {{math|''η''{{sub|0}}}} is a solution to {{math|1=''W''(''t''{{sub|0}}, ''t''{{sub|1}})''η'' = ''x''{{sub|1}} − ''ϕ''(''t''{{sub|0}}, ''t''{{sub|1}})''x''{{sub|0}}}} then a control given by {{math|1=''u''(''t'') = −''B''(''t''){{sup|T}}''ϕ''(''t''{{sub|0}}, ''t''){{sup|T}}''η''{{sub|0}}}} would make the desired transfer.

Note that the matrix {{math|''W''(''t''{{sub|0}}, ''t''{{sub|1}})}} defined as above has the following properties:
{{blist
| {{mvar|W}} is [[Symmetric matrix|symmetric]]
| {{mvar|W}} is positive semidefinite for {{math|''t''{{sub|1}} ≥ ''t''{{sub|0}}}} 
| {{mvar|W}} satisfies the linear [[Matrix differential equation|matrix differential equation]]
 <math display=block>\frac{\mathrm{d}}{\mathrm{d}t}W(t,t_1) = A(t)W(t,t_1)+W(t,t_1)A(t)^\mathsf{T}-B(t)B(t)^\mathsf{T}, \; W(t_1,t_1) = 0</math>
| {{mvar|W}} satisfies the equation<ref>{{cite book|first=Roger W.|last=Brockett|title=Finite Dimensional Linear Systems|publisher=John Wiley & Sons|year=1970|isbn=978-0-471-10585-5}}</ref>
 <math display=block>W(t_0,t_1) = W(t_0,t) + \phi(t_0,t)W(t,t_1)\phi(t_0,t)^\mathsf{T}</math>
}}

== Rank condition for controllability==
The [[Controllability Gramian]] involves integration of the [[State-transition matrix|state-transition matrix]] of a system. A simpler condition for controllability is a [[Rank (linear algebra)|rank]] condition analogous to the Kalman rank condition for time-invariant systems.

Consider a continuous-time linear system {{math|Σ}} smoothly varying in an interval {{math|{{bracket|''t''{{sub|0}}, ''t''}}}}:
<math display=block>\begin{align}
\dot{\mathbf{x}}(t) &= A(t) \mathbf{x}(t) + B(t) \mathbf{u}(t) \\
\mathbf{y}(t) &= C(t) \mathbf{x}(t) + D(t) \mathbf{u}(t).
\end{align}</math>

The state-transition matrix {{mvar|ϕ}} is also smooth. Introduce the {{times|{{mvar|n}}|{{mvar|m}}}} matrix-valued function {{math|1={{itco|''M''}}{{sub|0}}(''t'') = ''ϕ''(''t''{{sub|0}}, ''t'')''B''(''t'')}} and define
<math display=block>M_k(t) = \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 matrix {{mvar|M}} ({{mvar|M{{sub|i}}}}, for {{math|1=''i'' = 0, 1, ... , ''k''}}):
<math display=block>M^{(k)}(t) := \left[M_0(t), \ldots, M_k(t)\right] .</math>

If there exists a {{math|1=''t̅'' ∈ {{bracket|''t''{{sub|0}}, ''t''}}}} and a nonnegative integer {{mvar|k}} such that {{math|1=rank {{itco|''M''}}{{sup|(''k'')}}({{itco|''t̅''}}) = ''n''}}, then {{math|Σ}} is controllable.<ref name=":0">Eduardo D. Sontag, [https://books.google.com/books?id=f9XiBwAAQBAJ Mathematical Control Theory: Deterministic Finite Dimensional Systems].</ref>

If {{math|Σ}} is also analytically varying in an interval {{math|{{bracket|''t''{{sub|0}}, ''t''}}}}, then {{math|Σ}} is controllable on every nontrivial subinterval of {{math|{{bracket|''t''{{sub|0}}, ''t''}}}} [[If and only if|if and only if]] there exists a {{math|1=''t̅'' ∈ {{bracket|''t''{{sub|0}}, ''t''}}}} and a nonnegative integer {{mvar|k}} such that {{math|1=rank {{itco|''M''}}{{sup|(''k'')}}(''t{{sub|i}}'') = ''n''}}.<ref name=":0" />

The above methods can still be complex to check, since it involves the computation of the state-transition matrix {{mvar|ϕ}}. Another equivalent condition is defined as follow. Let {{math|1=''B''{{sub|0}}(''t'') = ''B''(''t'')}}, and for each {{math|''i'' ≥ 0}}, define
<math display=block>B_{i+1}(t) = A(t)B_i(t) - \frac{\mathrm{d}}{\mathrm{d} t}B_i(t). </math>
In this case, each {{mvar|B{{sub|i}}}} is obtained directly from the data {{math|(''A''(''t''), ''B''(''t''))}}. The system is controllable if there exists a {{math|''t̅'' ∈ {{bracket|''t''{{sub|0}}, ''t''}}}} and a nonnegative integer {{mvar|k}} such that {{math|1=rank({{bracket|''B''{{sub|0}}({{itco|''t̅''}}), ''B''{{sub|1}}({{itco|''t̅''}}), ... , ''B{{sub|k}}''({{itco|''t̅''}})}}) = ''n''}}.<ref name=":0" />

=== Example ===

Consider a system varying analytically in {{math|(−∞, ∞)}} and matrices
<math display=block>\begin{align}
A(t) &= \begin{bmatrix}
t & 1 & 0\\ 
0 & t^{3} & 0\\ 
0 & 0 & t^{2} 
\end{bmatrix}, \\
B(t) &= \begin{bmatrix}
0 \\ 
1 \\ 
1
\end{bmatrix}.
\end{align}</math>
Then 
<math display=block> [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 {{nowrap|<math>\mathbb{R}</math>.}}

=== Continuous linear time-invariant (LTI) systems ===

Consider the continuous linear [[Time-invariant system|time-invariant system]]
<math display=block>\begin{align}
\dot{\mathbf{x}}(t) &= A \mathbf{x}(t) + B \mathbf{u}(t) \\
\mathbf{y}(t) &= C \mathbf{x}(t) + D \mathbf{u}(t)
\end{align}</math>
where 
* {{math|'''x'''}} is the {{times|{{mvar|n}}|1}} state vector;
* {{math|'''y'''}} is the {{times|{{mvar|m}}|1}} output vector;
* {{math|'''u'''}} is the {{times|{{mvar|r}}|1}} input (or control) vector;
* {{mvar|A}} is the {{times|{{mvar|n}}|{{mvar|n}}}} state matrix;
* {{mvar|B}} is the {{times|{{mvar|n}}|{{mvar|r}}}} input matrix;
* {{mvar|C}} is the {{times|{{mvar|m}}|{{mvar|n}}}} output matrix; and
* {{mvar|D}} is the {{times|{{mvar|m}}|{{mvar|r}}}} feedthrough (or feedforward) matrix.

The {{times|{{mvar|n}}|{{mvar|nr}}}} controllability matrix is given by
<math display=block>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 (linear algebra)|rank]] (i.e. {{math|1=rank(''R'') = ''n''}}.

== Discrete linear time-invariant (LTI) systems ==

For a discrete-time linear state-space system (i.e. time variable {{nowrap|<math>k\in\mathbb{Z}</math>)}} the state equation is
<math display=block>\textbf{x}(k+1) = A\textbf{x}(k) + B\textbf{u}(k)</math>
where {{mvar|A}} is an {{times|{{mvar|n}}|{{mvar|n}}}} matrix and {{mvar|B}} is a {{times|{{mvar|n}}|{{mvar|r}}}} matrix (i.e. {{math|'''u'''}} is {{mvar|r}} inputs collected in a {{times|{{mvar|r}}|1}} vector). The test for controllability is that the {{times|{{mvar|n}}|{{mvar|nr}}}} matrix
<math display=block>\mathcal{C} = \begin{bmatrix}B & AB & A^{2}B & \cdots & A^{n-1}B\end{bmatrix}</math>
has full row [[Rank (linear algebra)|rank]] (i.e., {{math|1=rank({{itco|{{mathcal|C}}}}) = ''n''}}). That is, if the system is controllable, {{mathcal|C}} will have {{mvar|n}} columns that are linearly independent; if {{mvar|n}} columns of {{mathcal|C}} are linearly independent, each of the {{mvar|n}} states is reachable by giving the system proper inputs through the variable {{math|''u''(''k'')}}.

===Derivation===

Given the state {{math|'''x'''(0)}} at an initial time, arbitrarily denoted as {{math|1=''k'' = 0}}, the state equation gives {{math|1='''x'''(1) = {{itco|''A''}}'''x'''(0) + {{itco|''B''}}'''u'''(0)}},  then {{math|1='''x'''(2) = {{itco|''A''}}'''x'''(1) + {{itco|''B''}}'''u'''(1) =}} {{math|''A''{{sup|2}}'''x'''(0) + {{itco|''AB''}}'''u'''(0) + {{itco|''B''}}'''u'''(1)}}, and so on with repeated back-substitutions of the [[State variable|state variable]], eventually yielding
<math display="block">\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 display=block>\textbf{x}(n)-A^n\textbf{x}(0)= [B \, \, AB \, \, \cdots \, \, A^{n-1}B] [\textbf{u}^\mathsf{T}(n-1) \, \, \textbf{u}^\mathsf{T}(n-2) \, \, \cdots \, \, \textbf{u}^\mathsf{T}(0)]^\mathsf{T}.</math>

Imposing any desired value of the state vector {{math|'''x'''(''n'')}} 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|1=''n'' = 2}} and {{math|1=''r'' = 1}} (i.e. only one control input). Thus, {{mvar|B}} and {{mvar|AB}} are {{times|2|1}} vectors. If {{math|{{bracket|''B'' {{itco|''AB''}}}}}} has rank 2 (full rank), and so {{mvar|B}} and {{mvar|AB}} are linearly independent and span the entire plane. If the rank is 1, then {{mvar|B}} and {{mvar|AB}} are [[Line (geometry)|collinear]] and do not span the plane.

Assume that the initial state is zero.

At time {{math|1=''k'' = 0}}:
<math display=block>x(1) = A\textbf{x}(0) + B\textbf{u}(0) = B\textbf{u}(0)</math>
At time {{math|1=''k'' = 1}}:
<math display=block>x(2) = A\textbf{x}(1) + B\textbf{u}(1) = AB\textbf{u}(0) + B\textbf{u}(1)</math>

At time {{math|1=''k'' = 0}} all of the reachable states are on the line formed by the vector {{mvar|B}}.

At time {{math|1=''k'' = 1}} all of the reachable states are linear combinations of {{mvar|AB}} and {{mvar|B}}.
If the system is controllable then these two vectors can span the entire plane and can be done so for time {{math|1=''k'' = 2}}.

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 {{mvar|n}}, but the case of {{math|1=''n'' = 2}} is easier to visualize.

===Analogy for example of {{mvar|n}} = 2===
Consider an [[Engineering:Car analogy|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 {{mvar|C}} 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 {{mvar|C}} is 2.

If you change this example to {{math|1=''n'' = 3}} then the analogy would be flying in space to reach any position in 3D space (ignoring the orientation of the [[Engineering:Aircraft|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>\dot{\mathbf{x}} = \mathbf{f(x)} + \sum_{i=1}^m \mathbf{g}_i(\mathbf{x})u_i</math>

are locally accessible about <math>x_0</math> if the accessibility distribution <math>R</math> spans <math>n</math> space, when <math>n</math> equals the dimension of <math>x</math> and R is given by:<ref>Isidori, Alberto (1989). ''Nonlinear Control Systems'', p. 92–3. Springer-Verlag, London. {{ISBN|3-540-19916-0}}.</ref>

:<math>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>[\mathrm{ad}^k_{\mathbf{f}}\mathbf{\mathbf{g}}]</math> is the repeated [[Lie bracket of vector fields|Lie bracket]] operation defined by

: <math>[\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.

== Controllability via state feedback ==
When control authority on a [[Linear dynamical system|linear dynamical system]] is exerted through a choice of a time-varying feedback gain matrix <math>K(t)</math>, the system

: <math>\dot{\mathbf{x}} = (A - BK(t))\mathbf{x}</math>

is nonlinear, in that products of control parameters and states are present. The accessibility distribution <math>R</math> is, as before,

: <math>R= \begin{bmatrix} B & AB & \cdots & A^{n-1}B \end{bmatrix}. </math>

It is clear that for the system to be controllable, it is necessary that <math>R</math> has full column rank. It turns out that this condition is also sufficient. However, the (optimal) control strategy explained earlier needs to be modified so that the trajectory when applying an optimal input to steer the system between the specified states, does not pass through the origin, else the regulating input cannot be written in feedback form <math>u=-K(t)\mathbf{x}</math>. The controllability, as well as the strong controllability of this bilinear system
was proven in 
<ref name="abdelgalil_georgiou">{{cite journal |author1=Mahmoud Abdelgalil |author2=Tryphon T. Georgiou | title = Collective Steering in Finite Time: Controllability on <math>GL^+(n,\mathbb{R})</math> | journal = IEEE Transactions on Automatic Control | doi = 10.1109/TAC.2025.3574186 | year = 2025 }}</ref>.

== Collective controllability -- Control of the state transition via feedback ==

Collective controllability represents the ability to steer <math>n</math> linear dynamical systems that obey identical dynamics

: <math>\dot{\mathbf{x}}^{(i)}(t) = A \mathbf{x}^{(i)}(t) + B \mathbf{u}^{(i)}(t)</math>

where <math>n</math> equals the dimension of <math>\mathbf{x}</math>, 
between specified starting and ending configurations by way of a common state feedback gain matrix <math>K(t)</math>, and thereby, each instantiating a control input

: <math>\mathbf{u}^{(i)}(t)=K(t){\mathbf{x}}^{(i)}(t) </math>

for <math>i\in\{1,\ldots,n\}</math>, respectively.

The accessibility distribution <math>R</math> having full column rank is trivially a necessary condition. It is also sufficient, and in fact, the collective is strongly controllable, in that it can be steered from an initial
configuration

: <math>\Phi(0)= \begin{bmatrix} \mathbf{x}^{(1)}(0)\ldots \mathbf{x}^{(n)}(0) \end{bmatrix} </math>

to any specified terminal configuration

: <math>\Phi(T)= \begin{bmatrix} \mathbf{x}^{(1)}(T)\ldots \mathbf{x}^{(n)}(T) \end{bmatrix}, </math>
provided <math>\det(\Phi(0)\Phi(T))>0 </math>, over any specified time interval <math>[0,T] </math> through a choice of a common time-varying feedback gain matrix <math>K(t)</math> provided <math>R</math> has full column rank <ref name="abdelgalil_georgiou"/>.

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

This can be easily shown by controllable-uncontrollable decomposition.

== Output controllability ==

''{{dfn|Output controllability}}'' is the related notion for the output of the system (denoted {{math|'''y'''}} in the previous equations); the output controllability describes the ability of an external input to move the output from any [[Initial condition|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 {{math|1=''D'' = 0}} and matrix {{mvar|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 {{math|1=''D'' = 0}} and a {{mvar|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|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 {{mvar|A}}, {{mvar|B}}, {{mvar|C}}, and {{mvar|D}}, the {{times|{{mvar|m}}|{{math|({{itco|''n''}}+1)''r''}}}} ''{{dfn|output controllability matrix}}''
<math display=block>\begin{bmatrix} CB & CAB & CA^2B & \cdots & CA^{n-1}B & D\end{bmatrix}</math>
has full row rank (i.e. rank {{mvar|m}}) if and only if the system is output controllable.<ref name="Ogata97" />{{rp|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 [[Engineering:Actuator|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<ref>{{cite journal |author1=Claire J. Tomlin |author2=Ian Mitchell |author3=Alexandre M. Bayen |author4=Meeko Oishi|author4-link=Meeko Oishi | title = Computational Techniques for the Verification of Hybrid Systems | journal = Proceedings of the IEEE |volume=91 |issue=7 |pages=986–1001 | doi = 10.1109/jproc.2003.814621 | year = 2003 |bibcode=2003IEEEP..91..986T | url = http://www.cs.ubc.ca/~mitchell/Papers/publishedIEEEproc03.pdf | access-date = 2012-03-04
|citeseerx=10.1.1.70.4296 }}</ref> and [[Viability theory|viability theory]].<ref>{{cite book|author=Jean-Pierre Aubin |title= Viability Theory |year=1991 |publisher=Birkhauser |isbn=978-0-8176-3571-8 }}</ref>

== 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&ndash;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.<ref name="Polderman98">{{cite book|author1=Jan Polderman |author2=Jan Willems |title=Introduction to Mathematical Systems Theory: A Behavioral Approach|edition=1st|year=1998|publisher=Springer Verlag|location=New York|isbn=978-0-387-98266-3}}</ref>{{rp|151}}

== Stabilizability ==
A slightly weaker notion than controllability is that of '''stabilizability'''. A system is said to be [[Hautus lemma#Hautus Lemma for stabilizability|stabilizable]] when all uncontrollable state variables can be made to have [[Stability theory|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.<ref name="Anderson+Moore/book:1989">{{cite book|author1=Brian D.O. Anderson|author2=John B. Moore|title=Optimal Control: Linear Quadratic Methods|year=1990|publisher=Prentice Hall|location=Englewood Cliffs, NJ|isbn=978-0-13-638560-8}}</ref>

== Reachable set ==
Let {{math|''T'' ∈ '''''Т'''''}} and {{math|''x'' ∈ ''X''}} (where {{mvar|X}} is the set of all possible states and {{mvar|'''Т'''}} is an interval of time). The reachable set from {{mvar|x}} in time {{mvar|T}} is defined as:<ref name=":0" />
<math display=block>R^T{(x)} = \left\{ z \in X : x \overset{T}{\rightarrow} z \right\}, </math>
where {{math|''x'' {{overset|''T''|→}} ''z''}} denotes that there exists a state transition from {{mvar|x}} to {{mvar|z}} in time {{mvar|T}}.

For autonomous systems the reachable set is given by:
<math display=block>\mathrm{Im}(R)=\mathrm{Im}(B)+\mathrm{Im}(AB)+....+\mathrm{Im}(A^{n-1}B),</math>
where {{mvar|R}} is the controllability matrix.

{{math theorem |name= Claim |math_statement= In terms of the reachable set, the system is controllable if and only if {{nowrap|<math>\mathrm{Im}(R)=\mathbb{R}^n</math>.}}}}
{{math proof |title= Proof |proof=
We have the following equalities:
<math display=block>\begin{align}
R &= [B\ AB ....A^{n-1}B] \\
\mathrm{Im}(R) &= \mathrm{Im}([B\ AB ....A^{n-1}B]) \\
\mathrm{dim(Im}(R)) &= \mathrm{rank}(R)
\end{align}</math>
Considering that the system is controllable, the columns of {{mvar|R}} should be linearly independent. So:
<math display=block>\begin{align}
\mathrm{dim(Im}(R)) &= n \\
\mathrm{rank}(R) &= n \\
\mathrm{Im}(R) &= \R^{n}\qquad \blacksquare
\end{align}</math>
}}

A related set to the reachable set is the controllable set, defined by:
<math display=block>C^T{(x)} = \left\{ z \in X : z \overset{T}{\rightarrow} x \right\}.</math>
The relation between reachability and controllability is presented by Sontag:<ref name=":0" />
{{olist |type= lower-alpha
| An {{mvar|n}}-dimensional discrete linear system is controllable if and only if:
<math display=block>R(0)=R^k{(0)=X}</math>
(where {{mvar|X}} is the set of all possible values or states of {{mvar|x}} and {{mvar|k}} is the time step).
| A continuous-time linear system is controllable if and only if:
<math display=block>\begin{align}
R(0) &= R^e{(0)=X} \\
C(0) &= C^e{(0)=X}
\end{align}</math>
for all {{math|''e'' > 0}}.
}}

===Example===
Let the system be an n dimensional discrete-time-invariant system from the formula:
<math display=block>\phi(n,0,0,w)=\sum\limits_{i=1}^n A^{i-1}Bw(n-1),</math>
where {{mvar|ϕ}} (final time, initial time, state variable, restrictions) is defined as the transition matrix of a state variable {{mvar|'''x'''}} from an initial time {{math|0}} to a final time {{mvar|n}} with some restrictions {{mvar|w}}).

It follows that the future state is in {{math|''R''{{sup|''k''}}(0)}} if and only if it is in {{nowrap|<math>\mathrm{Im}(R)</math>,}} the image of the linear map {{mvar|R}}, defined as:
<math display=block>R(A,B)\triangleq [B\ AB ....A^{n-1}B],</math>
which maps
<math display=block>u^{n}\mapsto X.</math>

When {{math|1=''u'' = ''K''{{i sup|''m''}}}} and {{math|1=''X'' = ''K''{{i sup|''n''}}}} we identify {{math|''R''(''A'', ''B'')}} with a {{times|{{mvar|n}}|{{mvar|nm}}}} matrix whose columns are {{math|''B'', ''AB'', ... , ''A''{{i sup|''n''−1}}''B''}} in that order. If the system is controllable the rank of {{math|{{bracket|''B'' ''AB'' ... ''A''{{i sup|''n''−1}}''B''}}}} is {{mvar|n}}. If this is true, the image of the linear map {{mvar|R}} is all of {{mvar|X}}. Based on that, we have:
<math display=block>R(0)=R^k{(0)=X}</math> with <math>X\in\R^{n}.</math>

== See also ==
* [[Observability]]
* [[Hautus lemma]]

==Notes==
{{Reflist|group=note}}

==References==
{{reflist}}

==External links==
* [http://www.mathworks.com/help/toolbox/control/ref/ctrb.html MATLAB function for checking controllability of a system] {{Webarchive|url=https://web.archive.org/web/20120210033123/http://www.mathworks.com/help/toolbox/control/ref/ctrb.html |date=2012-02-10 }}
* [http://reference.wolfram.com/mathematica/ref/ControllableModelQ.html Mathematica function for checking controllability of a system]


[[Category:Classical control theory]]

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

{{Sourceattribution|Controllability}}