Next Article in Journal
On the Integration of Stiff ODEs Using Block Backward Differentiation Formulas of Order Six
Next Article in Special Issue
On a Semigroup Problem II
Previous Article in Journal
Securing Fingerprint Template Using Blockchain and Distributed Storage System
Previous Article in Special Issue
Critically-Finite Dynamics on the Icosahedron
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Observability and Symmetries of Linear Control Systems

by
Víctor Ayala
1,*,
Heriberto Román-Flores
1,
María Torreblanca Todco
2 and
Erika Zapana
2
1
Instituto de Alta Investigación, Universidad de Tarapacá, Casilla 7D, Arica 1000000, Chile
2
Departamento Académico de Matemáticas, Universidad, Nacional de San Agustín de Arequipa, Calle Santa Catalina, Nro. 117, Arequipa 04000, Peru
*
Author to whom correspondence should be addressed.
Symmetry 2020, 12(6), 953; https://doi.org/10.3390/sym12060953
Submission received: 7 April 2020 / Revised: 28 April 2020 / Accepted: 6 May 2020 / Published: 4 June 2020
(This article belongs to the Special Issue Symmetry and Dynamical Systems)

Abstract

:
The goal of this article is to compare the observability properties of the class of linear control systems in two different manifolds: on the Euclidean space R n and, in a more general setup, on a connected Lie group G. For that, we establish well-known results. The symmetries involved in this theory allow characterizing the observability property on Euclidean spaces and the local observability property on Lie groups.

1. Introduction

For general facts about control theory, we suggest the references, [1,2,3,4,5], and for general mathematical issues, the references [6,7,8].
In the context of this article, a general control system Σ can be modeled by the data Σ = ( M , D , h , N ) . Here, M is the manifold of the space states, D is a family of differential equations on M , or if you wish, a family of vector fields on the manifold, h : M N is a differentiable observation map, and N is a manifold that contains all the known information that you can see of the system through h. Normally, dim ( N ) dim ( M ) . One of the main problems here is to find topological, algebraic, and differentiable properties to impose on Σ , to recover the dynamics on M, through the partial information given by h .
The indistinguishable equivalence relation associated with the data permit decomposing the manifold into subsets where two points in the same set can not be distinguishable for D and h in N . More precisely, given two distinguishable states x , y M , there exist positive times t 1 , , t l > 0 with t = t 1 + + t l and a concatenation φ t = φ t 1 φ t l of flows associated with vector fields X 1 , , X l D such that:
h ( φ t ( x ) ) h ( φ t ( y ) ) .
The system is said to be observable if any indistinguishable class is a singleton set. In other words, any two points of M are distinguishable by Σ .
The goal of this article is just to compare the observability properties of the class of linear control systems in two different manifolds: on the Euclidean space R n and, in a more general setup, on a connected Lie group G. The symmetries involved in this theory allow:
  • characterizing the observability property on R n and also
  • characterizing the local observability property on G. After that, it is possible to get the global observability property on G through the fixed points of the drift vector field of the system.
In order to compare both systems, we establish well-known observability results in both cases. We cite the corresponding references.
Linear control systems on R n appear in many applications [9,10,11,12,13,14,15,16,17,18] for a related control system. Markus in [12] extended this notion from R n to a matrix group. Finally, in [19], the authors introduced the class of linear control system on any Lie group. In the references [18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35], the authors addressed relevant problems for this class of systems.
Furthermore, in [26], the author showed that each affine control system that generates a finite-dimensional Lie algebra on a connected manifold is diffeomorphic to a linear control system on a Lie group, or on a homogeneous space. This result shows the relevance of this class of systems on applications.
The article is organized as follows. In Section 2, we recall the notion of linear control systems Σ R n on Euclidean spaces. We establish the main observability results for Σ R n and the special symmetry contained in the data of the system, which completely determine this property. We finish the section with a couple of concrete examples. One of them is motivated by an example appearing in the book by Pontryagin [4]. Section 3 introduces the linear control system on connected Lie groups G. We show the face of the solution of this system and the notion of local observability and the (global) observability, and its main results appear in [27,28]. Before giving examples, we establish an algorithm to compute the relevant algebraic tools to decide the observability for the linear control system with observation Σ G O . We finish the section with a couple of examples. Finally, in the Appendix A, we give a short review on Lie theory to make the article more self-contained, [29].

2. Linear Control Systems on Euclidean Spaces

To introduce the definition of the linear control system on Euclidean spaces in this first section, we introduce elementary dynamic system notions about two classes of dynamics: linear and invariant.
Let us denote by R the set of real numbers. The Euclidean space R n = x = ( x 1 , , x n ) : x i R is generated by the basis vectors:
e i = ( 0 , , 1 ( i ) , , 0 ) = ( x i ) 0 , i = 1 , , n .
As an analytical differentiable manifold, the tangent space T x R n of R n at the state x R n is given by the vector space generated by:
( x i ) x : i = 1 , , n ,
where the vector ( x i ) x is obtained by the translation of e i by x .
A linear control system on R n essentially controls the behavior of a linear differential equation (linear vector field) through a number of constant vector fields, i.e., invariant (by translation) vector fields on R n .
In the sequel, we consider gl ( n , R ) the vector space of real matrices of order n . Any element A gl ( n , R ) determines a linear differential equation x · ( t ) = A x ( t ) with x ( 0 ) = x 0 as the initial condition. The corresponding solution reads as x ( t ) = e t A x 0 , where e A = Σ j = 0 A j j ! and A 0 = I d is the identity map on R n . Therefore, A induces the linear vector field:
x A x · ( t ) = A x ( t ) T x ( t ) R n
Then, the flow of the linear vector field X A defined by A reads as:
X t A : R n R n , X t A ( x ) = e t A x , t R .
On the other hand, any vector b = ( b 1 , , b i , , b n ) R n defines an invariant vector field X b on R n by:
X b ( x ) = Σ i = 0 n b i ( x i ) x T x R n ,
the solution of which is determined by its flows X t b ( x ) = x + t b , parallel lines to the line generated by b .
A linear control system on the Euclidean space R n is defined by the family of ordinary differential equations:
Σ R n : x · ( t ) = A x ( t ) + B u , u U ,
parametrized by the controls u U . Here, x R n , A gl n , R is a square real matrix of order n, B is an n × m real matrix, and U = L l o c 1 R , Ω is the family of admissible control functions. The set Ω R m is closed, and 0 i n t ( Ω ) . In other words, U is the space of locally integrable measure functions u : R Ω .
Consider a state x 0 R n and a control u U . The associated solution of Σ R n with these two parameters reads as:
φ t u x 0 = e t A x 0 + 0 t e τ A B u τ d τ .
In fact, φ t u x 0 is the only absolutely continuous function satisfying the Cauchy problem: [30]
φ t u x 0 · = A φ t u x 0 + B u , φ 0 u ( x 0 ) = x 0 ,
The solution set φ t u x 0 : t R describes a curve in R n in which elements are reached by x 0 through the specific dynamics u of Σ R n in positive and negative time. More generally, the positive orbit of Σ R n starting on any initial condition x 0 is given by:
S x 0 = φ t u x 0 : 0 t , u U .
Thus, S x 0 is the set of the states of R n that can be reached from x 0 in positive time by using all the control functions u U .

2.1. Observability

A linear control system with observation Σ R n O is a linear control system on R n given by the data contained in Formula ( 1 ) and the observation map h : R n R s . Here, s n and h is a linear transformation [30,31].
Our interest is in analyzing the observability property of Σ R n O and a special symmetry property of the pair ( A , h ) determining this notion. All the results in this section are very well known [30]. However, we need to include them in order to compare with the Lie group case.
In the sequel, we introduce the Σ R n O -indistinguishable relationship in R n . This equivalence relation allows us to decompose the space of states into equivalence classes, which contain those elements that cannot be distinguished from each other in R s by the observation map h : R n R s through every positive dynamic of Σ R n O . Precisely, we denote by C the matrix of h in the canonical basis. The linear map h does not distinguish the state x 1 from x 0 through any positive dynamic of Σ R n O , meaning that for any control u U and for each time t 0 :
C φ t u x 1 = C φ t u x 0 .
This condition is equivalent to:
C ( e t A x 0 + 0 t e τ A B u τ d τ ) = C ( e t A x 1 + 0 t e τ A B u τ d τ )
C e t A x 0 = C e t A x 1 , for each t 0
in particular, when t = 0 , C x 0 = C x 1 . i Recall that the kernel and the image of h , given respectively by:
K e r ( h ) = x R n : h ( x ) = 0 , Im ( h ) = y R s : x R n with y = h ( x )
are subspaces of R n .
Definition 1.
Let Σ R n O be a linear control system with observation:
1. 
Two states x 0 , x 1 R n are indistinguishable for Σ R n O , denoted as x 0 I x 1 , if:
x 1 x 0 K e r C e t A , t 0
2. 
Σ R n O is said to be observable from x R n if the equivalence class I ( x ) of x is x .
3. 
Σ R n O is said to be observable if it is observable from x for each x R n .
The main results of this section show that the observability property of linear control systems with observation totally depends on a symmetric matrix built with A and C . The next result shows some very well-known properties of I 0 , which determine completely the observability property of the system; see [30,31]. Since the proof is simple, we include it because it will inspire obtaining symmetric properties on the Lie group case.
Proposition 1.
Let Σ R n O be a linear control system with observation. Then,
1. 
I is an equivalence relation;
2. 
I 0 = t 0 K e r C e t A ;
3. 
I x = x + I 0 .
Proof. 
To prove the first item, we use the fact that for any t 0 , the set K e r C e t A R n is a vector subspace. It follows that x I x and also x I y imply y I x . Thus, I is reflexive and symmetric. Now, assume x 0 I x 1 and x 1 I x 2 . The equality:
x 2 x 0 = ( x 1 x 0 ) + ( x 2 x 1 ) K e r C ( e t A ) , t 0
shows that I is also transitive.
To prove 2, consider x I 0 . By definition, the analytical curve γ ( t ) = C ( e t A ) x = 0 for any t 0 . Taking the first derivative of γ , we obtain:
( d d t ) t C ( e t A ) x = C A ( e t A ) x = 0 .
When t = 0 , we get x K e r ( C A ) . By taking the derivative of C A ( e t A ) x = 0 at t = 0 , we obtain C A 2 x = 0 . Continuing this process, we obtain:
I ( 0 ) j = 0 n 1 K e r C A j .
On the other hand, if x j = 0 n 1 K e r C A j , it follows immediately that:
x K e r C e t A , t 0 .
In fact, from the Cayley–Hamilton theorem [30], A l S p a n A 1 , A 2 , , A n 1 for l n . Since, K e r C A j x = 0 , for j = 1 , , n 1 , the conclusion follows from the definition of the exponential map e t A = Σ j = 0 t j j ! A j and A 0 = I d .
Finally, let x R n be arbitrary by fixing y I ( x ) . By definition:
y x K e r C ( e t A ) , t 0 .
Thus,
y x I 0 y x + I 0 I ( x ) x + I ( 0 ) .
Let y I ( 0 ) . Since, y = ( x + y ) x I ( 0 ) ; then ( x + y ) + I x ; therefore, x + y I ( x ) , concluding the proof.
From the previous proposition, any indistinguishable equivalence class I ( x ) depends just on x and I 0 . Therefore, we get the main results of this section:
Theorem 1.
Let Σ R n O be a linear control system with observation:
Σ R n O i s o b s e r v a b l e I ( 0 ) = 0 .
 □

2.2. The Symmetric Matrix Determining Observability on R n

Next, we show that the observability property of a linear system with observation depends strongly on the symmetric Gramian observability matrix, [30]
G 0 T = 0 T e τ A * C * C e τ A d τ ,
where P * denotes the transpose of the matrix P.
Theorem 2.
Let Σ R n O be a linear control system with observation. Therefore, for any T > 0 :
I 0 = K e r G 0 T .
Proof. 
The proof is independent of T. If x K e r G 0 , hence x * G 0 x = 0 . Thus,
x * ( 0 T e τ A * C * C e τ A d τ ) x = 0 T C e τ A x 2 d τ = 0 .
Since the application C e τ A x : 0 , T R n is analytical, the continuous function C e τ A x 2 must be null on the interval 0 , T . In fact, if there exists τ with C e τ A x 2 > 0 , by continuity, there exits ϵ > 0 and a neighborhood ( τ ε , τ + ε ) where C e t A x 2 > 0 . Hence,
0 T C e τ A x 2 d τ > 0
and we have a contradiction. Therefore,
γ τ = C e τ A x = 0 , 0 τ T .
However, γ is an analytical application, which vanishes on the open set ( 0 , T ) . Therefore, γ must be null on the whole real line. In particular, C e t A x = 0 , for each t R . As before, we obtain x I 0 .  □
Reciprocally, if x is indistinguishable from the origin,
G 0 x = 0 T e τ A * C * C e τ A d τ = 0 T e τ A * C * 0 d τ = 0 ,
and the proof is complete.
Remark 1.
Since for any T > 0 , the matrix G 0 T is symmetric, it is possible to consider its orthogonal space:
K e r G 0 T = I m G 0 T .
In particular, I m G 0 T is the space of the observable states of the system. On the other hand, any x R n has a unique decomposition:
x = x 0 + x I , x 0 x I
where x 0 is the observable component of x and x I is the non-observable part.
To analyze observability through the Gramian matrix, you need to compute an integral. However, it is possible to avoid this calculus through the following algebraic characterization.
Theorem 3.
Let Σ R n O be a linear control system with observation. Denote by O M s n × n R the matrix:
O = ( C C A C A 2 C A n 1 ) T r a n s p o s e d .
Then, I 0 = K e r ( O ) .
Proof. 
The proof follows directly from:
x I 0 C A j x = 0 , j = 0 , 1 , . . , n 1 C e t A x = 0 , t 0 .
As a consequence of the previous analysis, we obtain the following directly: □
Corollary 1.
Let Σ R n O be a linear control system with observation and T > 0 . The following are equivalent:
1. 
Σ R n O is observable;
2. 
G 0 T is invertible for any T > 0 ;
3. 
R a n k ( O ) = n .

2.3. Examples

The first example shows that it is possible to build trajectories through the observation function when the system is observable. This means, in this case, that it is possible to determine the solutions of the system on R n just with s < n parameters.
Example 1.
From an observation function, recover the optimal trajectories of the following problem:
S t o p p i n g a t r a i n a t t h e r a i l w a y s t a t i o n i n t h e s h o r t e s t t i m e .
This very famous optimal problem appears in the book [4]. We extend this model with a couple of observation maps.
Let us consider a railway modeled by the real line and without friction. We denote by x ( t ) the distance between the train and the origin 0 R considered as the station. Then, we get on R 2 the linear control system:
Σ R 2 : x ˙ ( t ) y ˙ ( t ) = 0 1 0 0 x ( t ) y ( t ) + 0 1 u ( t ) , u U .
In fact, from Newton’s law, F = M u , where the function u is the acceleration. Of course, we can consider M = 1 , then we get:
x · ( t ) = y ( t ) is the velocity and y · ( t ) = u ( t ) the acceleration .
The train is controlled by the admissible class of control U = L l o c 1 R , Ω = 1 , 1 , and the set
U = u : R Ω , u is locally integrable .
Here, 1 u ( t ) 1 with 1 and 1 representing the maximum and minimum possible accelerations and breaking, respectively. It is interesting to observe the geometric meaning of the problem: to solve it, we need first to associate with any initial state z 0 = ( x 0 , y 0 ) R 2 a curve solution of the system Σ R 2 transporting z 0 to the origin ( 0 , 0 ) .
R a n k ( O ) = 0 1 1 0 = 2 . On the other hand, the Lyapunov spectrum S p e c L y ( A ) of A, built with the real parts of the eigenvalues of A is equal to 0 . According to Example 3.2.16 in the book of Colonius and Kliemann [1], the system is controllable, which means that given any two points of the plane, each of them can be reached from the other through a solution of the system in positive time. In particular, the optimal problem is well posed. Hence, from a trivial computation, we know that for any state z, there exists u U and a +time t > 0 such that φ 0 u ( z 0 ) = z 0 and φ t u z 0 = ( 0 , 0 ) . However, we do not know what is the curve that does the job. Worst, among all the curves transporting z 0 to the origin, how doe we compute the trajectory with the shortest time? Fortunately, by the Pontryagin maximum principle [4], the optimal control u * exists, and its optimal value belongs to the boundary Ω = 1 , 1 . Furthermore, according to this principle, the optimal trajectories are concatenations with at most one change (i.e., the dimension of the manifold n minus one), of two main curves: the solutions of u = 1 and u = 1 . Each of the differential equations induced by these special controls decomposes the plane into (parallel) parabolas with respect to the y-axis, solving the problem completely. It is worth saying that this principle [4] won a relevant prize in Russia. Just observe that you reduce the possible values of the optimal trajectory from the interval 1 , 1 to the discrete set 1 , 1 ; from non-numerable infinity values to two possibilities.
Now, we apply the observability results to show that for the train problem, it is more relevant to know the first variable rather than the second one, i.e., knowing the distances gives you more information than the velocities.
  • Consider Σ R 2 and h = π 1 : R 2 R , h ( x , y ) = x the projection on the first variable. A simple computation shows that R a n k ( O ) = 2 . In fact,
    C = ( 1 0 ) , A = 0 1 0 0 , C A = ( 0 1 ) then O = 1 0 0 1 .
    According to Theorem 3, the system is observable.
  • Consider now the same system Σ R 2 and h = π 2 : R 2 R with h ( x , y ) = y the projection on the second variable. In this case:
    C = ( 0 1 ) , C A = ( 0 0 ) then O = 0 0 1 0 .
    and R a n k ( O ) = 1 . Therefore, the system is not observable. Furthermore, I ( 0 , 0 ) coincides with the x-axis. For any z R 2 , I ( z ) = z + R .
Geometrically, this situation can be explained as follows. In the case of π 1 , any two states with different first coordinates can be distinguished directly from their own observable function. Therefore, take any two states with the first coordinate equal to zero. By applying simultaneously the solutions coming from u = 1 to these two initial conditions, the projection of both curves will coincide. Now, by considering u = 1 , the projections of the curves will give different trajectories on the x-axis. Therefore, the system is observable since I ( 0 ) = 0 .
In the case of π 2 , consider two different states on the x-axis. The projections to the y-axis of the mentioned parabolas with u = 1 and u = 1 , will coincide. Therefore, I ( 0 ) = x -axis, and the system is not observable. In particular, through the knowledge of the velocities, it is not possible to recover the solution in the manifold. More precisely, contrary to the first case, here we cannot determine the initial condition and therefore all the trajectory!
Example 2.
The next example corresponds to a linearization of a non-linear control system modeling the attitude of a satellite in orbit around the Earth (see [2]) with different observation maps.
Let Σ R n be a linear control system with:
A = 0 1 0 0 3 β 2 0 0 2 β 0 0 0 1 0 2 β 0 0 , B = 0 0 1 0 0 0 0 1 , β 0 .
  • Consider the linear transformation h with matrix C = 1 0 0 0 0 0 1 0 . A simple computation shows that O is invertible. Hence, the system is observable.
  • If C = 1 0 0 0 : R 4 R , R a n k ( O ) = 3 , and the system is not observable. Therefore, by knowing just the first variable, it is not possible to rebuild the system.
  • If C = 0 0 1 0 : R 4 R , R a n k ( O ) = 4 , and the system is observable. Hence, the third variable contains enough information to recover any Σ -trajectories from a given initial condition.

3. Linear Control Systems on Lie Groups

The Appendix at the end of the article contains a short review of Lie theory; mainly, a basic notion about Lie groups, Lie algebras, and some examples. The reference [29] is a very good starting point.
We begin this section by extending the definition of a linear control system from the Euclidean space G = R n to a connected arbitrary Lie group G . From a dynamic point of view, a linear control system on the Abelian group R n is determined by:
the drift A gl ( n , R ) , and the B - column vectors b j R n , j = 1 , , m .
We start by the drift vector field A. For that, we need to extend the notion of the linear differential equation x · ( t ) = A x ( t ) from vector spaces to groups.
The first idea comes from the flow e t A : R n R n , t R , induced by matrix A , which is a one-parameter group of R n -automorphism. In fact, for any real number t, the linear map e t A is the inverse of e t A . As a first approach, we impose on X the same property. This means that the flow φ t of the drift X :
φ t : t R A u t ( G )
must be a subgroup of A u t ( G ) , the Lie group of the automorphism of G . We call X an infinitesimal automorphism. Therefore, the drift is determined by the formula:
X ( g ) = ( d d t ) t = o φ t ( g ) .
The Lie bracket between a linear and an invariant dynamic on R n reads as:
A x , b = A b R n .
Thus, the linear vector field determined by A leaves invariant the Lie algebra T 0 R n = R n of R n . Therefore, we impose that X leaves invariant the Lie algebra g of G under the Lie brackets. In [19], the authors introduced the following notion:
Definition 2.
A vector field X on G is said to be linear if:
X , Y g , f o r a n y Y g , a n d X ( e ) = 0 .
In [25], the author proved the following:
Theorem 4.
Let X be a linear vector field on G. The following are equivalent:
1. 
X is linear;
2. 
X is an infinitesimal automorphism.
On the other hand, any vector b R n determines by translation an invariant vector field X b on R n , as follows:
X b ( x ) = Σ i = 0 n b i ( x i ) x T x R n
Since,
B u = Σ j = 1 m u j b j
the linear control system Σ R n on R n controls the matrix A through m invariant vector fields X b 1 , , X b m .
Let G be a connected Lie group. In [19], the authors introduced the class of linear control systems Σ G on G through the controlled family of differential equations:
Σ G : g ˙ ( t ) = X ( g ( t ) ) + j = 1 m u j ( t ) X j ( g ( t ) ) , u U ,
parametrized as u U . Here, the drift X is linear, X 1 , , X m are invariant vector fields, and u = ( u 1 , , u m ) belongs to the class of admissible control functions U = L l o c 1 R , Ω .
If G = R n , we recover the definition of the linear control system in Euclidean spaces. In other words, Formula ( 2 ) is a perfect generalization of Formula ( 1 ) .
Remark 2.
Any linear vector field X determines a derivation D = D X of g , defined as:
D : g g s u c h t h a t D X = [ X , X ] ( e ) , f o r a l l X g .
Recall that [ A x , b ] = A b as in the Euclidean case. As a matter of fact, for each t R , ( d φ t ) e = e t D . In particular, from the commutative diagram, we obtain:
g ( d φ t ) e g exp exp G φ t G
φ t ( exp X ) = exp ( d φ t ) e X = exp ( e t D X ) , f o r a l l t R , X g .
Furthermore, assume G is also a simply connected and nilpotent Lie group. It is well known that in this case, the exponential map is a diffeomorphism. Consider a derivation D , and for g G , take X g such that log ( g ) = X . Therefore,
φ t ( g ) = exp ( e t D log ( g ) ) , for all t R , g G .
Finally, consider the inner derivations D g , induced by an element Y of g by the formula D ( · ) = · , Y . Since it flows:
φ t ( g ) = exp ( t Y ) g exp ( t Y ) , for all t R
is given by conjugation, and it is easy to determine the value of X = X D on G through the formula:
X = ( d d t ) t = 0 φ t .

3.1. The Solution of a Linear Control System on Lie Groups

Let Σ G be a linear control system on the connected Lie group G . For a control u U , the solution of:
g ( t ) · = X ( g ( t ) ) + j = 1 m u j ( t ) X j ( g ( t ) )
with arbitrary initial condition x G is given by:
γ t ( x ) = φ t ( β ( t , u ) x ) ,
where β ( t ) is a differentiable curve in g . Precisely, this curve satisfies the differential equation: [19]
β ( t , u ) · = ( φ t ) * ( j = 1 m u j ( t ) X X ( β ( t ) , u ) j ) .

3.2. Observability Properties

In this section, we state the main ideas and results appearing in [28]. First, we introduce the notion of a linear control system with an observation map on a Lie group. Through the shape of the solution of Σ G given in the previous section, we are willing to show that as in the Euclidean space, it is possible to characterize the observability property without the use of the control vectors. Actually, the algebraic condition just depends on the drift vector field and the observation map. Therefore, we have here the first symmetric situation between both classes of systems.
In this section, we will not show the proof of the main results. In fact, the mathematics involved are hard and out of the scope of this article.
For any homomorphism h : G H between Lie groups, the kernel, and the image of h:
K e r ( h ) = x G : h ( x ) = e H and Im ( h ) = y H : x G with h ( x ) = y
are subgroups of G and H, respectively. Here, e H denotes the identity of H .
Definition 3.
Let G be a connected Lie group. A linear control system with observation Σ G O on G is given by the data contained in Formula ( 2 ) and the observation map h : G h o m H , where H is a Lie group with dim ( G ) dim ( H ) .
We start with the indistinguishable and observable notions of Σ G O .
Definition 4.
Let Σ G O be a linear control system with observation on G:
1. 
Two elements g 1 and g 2 are indistinguishable, denoted by g 1 I g 2 , if:
φ t ( g 1 g 2 1 ) K e r ( h ) , t 0
2. 
Σ G O is said to be observable from g G if the equivalence class I ( g ) is g .
3. 
Σ G O is said to be observable if it is observable from g G , for every g G .
4. 
Σ G O is said to be locally observable from g G if U ( g ) , a neighborhood of g , with:
I ( g ) U ( g ) = g .
5. 
Σ G O is said to be locally observable if it is locally observable from g , g G .
Remark 3.
In the Euclidean case, the equivalence class of the origin I ( 0 ) is a subspace of R n . This means I ( 0 ) is trivial or at least contains a line. In the non-trivial case, the local notion above does not have meaning. However, in the group case, I ( e ) could be discrete and even finite, here e = e G . Thus, the locally observable notion has perfect meaning. This is a fundamental difference between both Lie groups R n and G .
The formula φ t ( g 1 g 2 1 ) K e r ( h ) , t 0 , associated with the indistinguishable notion comes from the following analysis. Recall that φ t is an infinitesimal automorphism and h is a homomorphism. Therefore, given two state g 1 and g 2 and any control u U , we have:
g 1 I g 2 h ( φ t ( β ( t , u ) g 1 ) = h ( φ t ( β ( t , u ) g 2 )
h ( φ t ( g 1 ) ) = h ( φ t ( g 2 ) ) h ( φ t ( g 1 g 2 1 ) = e H .
In particular, by taking g 2 = e , we obtain:
I = I ( e ) = g G : φ t ( g ) K e r ( h ) , t 0 .
Furthermore, for any g G , I ( g ) = I g . In fact,
g 1 I g 2 g 2 I g 1 .
In the sequel, we denote Σ G O also by ( G , X , h , H ) .
Proposition 2.
Consider the system Σ G O = ( G , X , h , H ) . Then [28]:
1. 
I is a topologically closed and normal subgroup of G;
2. 
I = g G : φ t ( g ) K e r ( h ) , t R .
Coming back to the Euclidean case, certainly, any subspace of R n is topologically closed. Since R n is an Abelian group, any subgroup is normal, actually Abelian as well. On the other hand, Proposition 1 implies that we have a complete symmetry between the indistinguishable class of the identity element e of G and the neutral element 0 R n ,
I = x R n : φ t ( x ) K e r ( h ) , t R .

3.3. Local Observability

As we mentioned, a continuous Lie group could contain discrete subgroups. For instance, for any natural number k, the set:
R ( k ) = z C : z k = 1
is a finite subgroup of the circle S 1 . For example, R ( 4 ) = π 2 , π , 3 π 2 , 1 .
In the sequel, we show an algebraic version of Proposition 17. Let I and K be the Lie algebras (the tangent space at the identity element e) of I and the kernel K = K e r ( h ) of h, respectively. Define a d i ( X ) : g g inductively:
a d 0 ( X ) = I d , a d 1 ( X ) ( Y ) = X , Y , , a d i ( X ) = a d i 1 ( X ) ( X , Y ) .
Proposition 3.
Consider the system Σ G O = ( G , X , h , H ) . Then [28],
1. 
a d i ( X ) ( I ) I , for any natural number i 0 ;
2. 
I = i = 0 n 1 a d i ( X ) ( K ) .
The first condition means that the Lie bracket of any order between the drift X with elements in I will remain in I . Precisely, let Y I , then:
X , Y I , X , X , Y I , X , X , X , Y I , and so on .
The idea of the proof of 2 is as follows. Since I K , we need to prove:
i = 0 n 1 a d i ( X ) ( K ) K ,
which follows from 1 .
Now, we are in a condition to state the main results about the observable system from the local point of view.
Theorem 5.
Consider the system Σ G O = ( G , X , h , H ) . Then [28],
Σ G O i s l o c a l o b s e r v a b l e Σ G O i s l o c a l o b s e r v a b l e f r o m e I = 0 .
This result says that the observability property of Σ R n O on R n and the local observability property of Σ G O on G can be characterized in a symmetric way, i. e., with analogous mathematical tools. In fact,
I = 0 is equivalent to I = e .
In fact, in the Euclidean case e = 0 . Therefore, we recover Theorem 3 ,
Σ R n O is observable I = 0 .

3.4. Observability

The local observability property of Σ G O just depends on the indistinguishable equivalence class of e . To obtain an observability result, we need to consider something more: the set of fixed points of the drift vector field X , given by:
F i x ( X ) = g G : φ t ( g ) = e , for any t R .
Theorem 6.
Consider the system Σ G O = ( G , X , h , H ) . Then [28],
Σ G O i s o b s e r v a b l e I = 0 a n d F i x ( X ) K e r ( h ) = e .
Proof. 
Idea of the proof: Since I = e , we get I = 0 . If g F i x ( X ) K e r ( h ) ,
φ t ( g ) = g K e r ( h ) , t R .
By definition, g I , then g = e .
Reciprocally, consider the map:
Φ : R × G G defined by Φ ( t , g ) = φ t ( g ) .
We observe that Φ : R × I I is well defined. Pick g G , and define the application Φ g : R × g G . The set:
Φ g ( R × g ) = Im ( Φ g ) = φ t ( g ) : t R
is connected since the domain R × g is connected and Φ g is continuous. By taking t = 0 , we get g Im ( Φ g ) I .
I = 0 implies that I is a discrete subgroup of G . Therefore,
g = Im ( Φ g ) so g F i x ( X ) .
However, g I K e r ( h ) . It turns out that g = e , ending the proof. □

3.5. An Observability Algorithm

Let us denote by g the dual space of g . In order to compute the Lie algebra I , it is suitable to use a general algorithm provided by Isidori [2]. Starting from the Lie algebra K , we build a finite sequence of invariant g -subspaces, which converges to I ; see [28] for details. Let Σ G O = ( G , X , h , H ) be a linear control system with observation. Consider the following steps:
  • Choose a basis B = { Y 1 , , Y p } of K ;
  • Compute the subspace B = { w 1 , , w n p } of g ;
  • Compute the a d ( X ) ( B ) -basis to I , i.e.,
    a d ( X ) ( B ) = { a d i ( X ) ( w j ) 0 i , 1 j n p } ,
    a d 0 ( X ) = I d , a d ( X ) ( w ) = w a d ( X ) ,
    a d i ( X ) ( w ) = a d ( a d i 1 ( X ) ( w ) ) , i 2 .
The previous algorithm allows computing the relevant algebraic tool I .
Theorem 7.
Let us consider the system Σ G O = ( G , X , h , H ) . Therefore, [28]
I = ( S p a n a d ( X ) ( B ) )
Proof. 
It follows from Isidori‘s theorem in [2]. In fact, in this particular case, we obtain:
S p a n ( a d ( X ) ( B ) ) = I .
Recall that I can be computed by integration of the Lie algebra I .  □

3.6. Examples

In the sequel, we use the previous algorithm to compute some examples.
Example 3.
The Heisenberg–Lie group G of dimension three reads as:
G = g = 1 a c 0 1 b 0 0 1 : a , b , c R .
Let us denote by e i j * the coefficient of a matrix such that the only non-null element is e i j = 1 . The Lie algebra g of G is:
g = S p a n X 1 = e 12 * , X 2 = e 23 * , X 3 = e 13 * ,
where the only non-null bracket is X 1 , X 2 = X 3 . We consider the infinitesimal automorphism X induced by the matrix X 1 by conjugation:
X ( x ) = d d t t = 0 exp ( t X 1 ) g exp ( t X 1 ) = b X 3 .
With these data, we develop two examples:
  • Let Σ G O = ( G , X , h = π 2 , R ) be a linear control system on G, where:
    π 2 1 a c 0 1 b 0 0 1 = b .
    We have,
    K e r ( h ) = 1 a c 0 1 0 0 0 1 : a , c R and K = S p a n X 1 , X 3 .
    Thus, the β basis is X 1 , X 3 . It follows that the orthogonal basis β is generated by w = ( X 2 ) . Simple calculus shows that:
    X , X 1 = 0 , X , X 2 = X 3 and X , X 3 = 0 .
    Since a d ( X ) ( w ) = w a d ( X ) , we obtain:
    a d ( X ) ( w ) ( X 1 ) = 0 , a d ( X ) ( w ) ( X 2 ) = ( X 2 ) ( X 3 ) = 0 and a d ( X ) ( w ) ( X 3 ) = 0 .
    Thus, a d ( X ) ( w ) = 0 and a d i ( X ) ( w ) for any i . It follows that:
    0 = a d i ( X ) ( β ) K I = K .
    Consequently, Σ G O is not observable since it is not locally observable.
  • Let Σ G O = ( G , X , h , G / exp ( t X 1 ) be a linear control system on G, where:
    h : G G / exp ( t X 1 ) is the canonical projection .
    We have,
    K e r ( h ) = 1 0 0 0 1 b 0 0 1 : b R and K = S p a n X 2 .
    Therefore, β = X 2 and β = w 1 = ( X 1 ) , w 2 = ( X 3 ) . However, a d ( X ) ( w ) = ( X 2 ) . Then, I I = g I and I = 0 , which implies by Theorem 19 that the system Σ G O is locally observable.
On the other hand, the set of fixed points on the drift:
F i x ( X ) = 1 a c 0 1 0 0 0 1 : a , c R , in fact
φ t ( g ) = exp ( t X 1 ) g exp ( t X 1 ) = g b = 0 .
It turns out that the only element in F i x ( X ) and K e r ( h ) is the identity. From Theorem 6, the system Σ G O = ( G , X , h , G / exp ( t X 1 ) is observable.

4. Conclusions

In this paper, we first showed that the symmetric property of the Gramian matrix determines totally the observability property of any linear control system on a Euclidean space.
Second, the natural extension of a linear control system from R n to a Lie group G contains many symmetries. In fact, the drift A of Σ R n with flow e t A A u t ( R n ) , t R , can be generalized to the drift X of Σ G through the infinitesimal automorphism notion. Precisely, its flow φ t A u t ( G ) , for any t R . Furthermore, the column vectors of the matrix B of Σ R n determine m constant vector fields on R n , which can be generalized as m invariant vector fields on the Lie group G.
On the other hand, we proved that there are many symmetries in the given data, which allow characterizing the observability properties of Σ R n O and Σ G O . For instance,
  • The equivalent class I ( 0 ) of 0 R n and I ( e ) of e G have symmetric constructions;
  • The algebraic tool I allows characterizing the observability property of Σ R n O . The natural extension of I = I ( 0 ) from a subspace to a subgroup I = I ( e ) characterizes the locally observable property of Σ G O .
It is worth saying that contrary to the Abelian group R n , a Lie group G could contain discrete subgroups. To approach this relevant difference, we considered the sets of fixed points of the drift vector fields X . With that, we characterized the observability of Σ G O .

Author Contributions

Supervision, V.A.; writing, original draft, M.T.T. and H.R.-F.; writing, review and editing, E.Z. All authors have read and agreed to the published version of the manuscript.

Funding

Se agradece a la Universidad Nacional de San Agustín de Arequipa, UNSA, por el financiamiento del proyecto de investigación según Contrato N° IAI-014-2018-UNSA.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. A Short Review of Lie Theory

To define the class of linear control systems in a more general setting, we introduce some basic notion of Lie theory; mainly, matrix Lie groups, Lie algebras, and two special classes of dynamics: invariant and linear vector fields. For general facts of this matter, we mention the following references [6,7,29,32].
A group G is said to be a topological group if the applications μ : G × G G , defined by μ ( g , h ) = g h and ι : G G , determined by ι ( g ) = g 1 are continuous.
According to the solved Hilbert’s fifth problem, G has a structure of differential manifolds, and μ and ι are analytical applications. We call this structure a Lie group. In particular, it is possible to define the differential equation on G (vector fields) and to compute its solution through its flows.
Example A1.
The following sets are Lie groups with the natural structure:
1. 
The Euclidean space ( R n , + ) .
2. 
The multiplicative set G L ( n , R ) = A : det ( A ) 0 of invertible real matrices of order n. We denote by G L + ( n , R ) the connected component, which contains the identity map I d .
3. 
The n-dimensional torus T n = S 1 × . × S 1 , n-times, where S 1 = z C : z = 1 is the circle group.
4. 
The special linear group:
S L ( n , R ) = det 1 ( 1 )
and the special orthogonal group:
S O ( n , R ) = { A O ( d ) A A t = I d a n d det ( A ) = 1 }
5. 
A matrix representation of the Heisenberg–Lie group reads as:
G = g = 1 a c 0 1 b 0 0 1 : a , b , c R .
Define the right and the left translations maps on G by:
R g : G G , R g ( x ) = x g and L g : G G , L g ( x ) = g x
respectively. Their inverses are given by R g 1 and L g 1 .
Definition A1.
An invariant vector field X on G is determined by the formula X g = ( R g ) * ( X e ) , g G .
Here, ( R g ) * denotes de derivative of R g at the identity element e of G . Just observe that any invariant vector field is determined by its value at e .
The notion of Lie brackets is a measure of the non-commutative behavior between any two dynamics. In the case of a matrix group, the Lie bracket corresponds to the commutator, A , B = A B B A .
Here is an example coming from real life. Imagine an airplane flying always at the same altitude and moving 1000 km in every direction: first to the north, second to the west, then to the south, and finally, to the east. Except for some special states, the initial and ending points of the trip are always different! The plane moves through the solution of two vector fields on the sphere S 2 = x R 3 : x = 1 induced by rotational matrices A , B S O ( d , R ) . In this case, A , B 0 .
Denote by g the vector space of right invariant vector fields on G. As a matter of fact,
X , Y g X , Y g .
Therefore, g is isomorphic to the tangent space T e G , and it is called the Lie algebra g of G. The elements of g satisfy the skew-symmetric property:
[ X , Y ] = [ Y , X ] , for any X , Y g
and the Jacobi identity:
a d ( X ) ( [ Y , Z ] ) + a d ( Z ) ( [ X , Y ] ) + a d ( Y ) ( [ Z , X ] ) = 0 , if X , Y , Z g
Moreover, any subspace V g is a subalgebra if:
X , Y V X , Y V
V is an ideal if:
X V and Y g X , Y V .
Next, we show the algebra g of some groups through T e G :
  • T 0 R n = R n ;
  • T I d G L + ( n , R ) = gl ( n , R ) ;
  • T e S 2 = R 2 ;
  • T I d S O ( n , R ) = { A gl ( n , R ) A + A t = 0 } , the skew-symmetric matrices;
  • sl ( n , R ) = { A gl ( n , R ) t r ( A ) = 0 } , the trace zero matrices;
  • The Lie algebra of the Heisenberg group has the basis X 1 , X 2 , X 3 such that X 1 , X 2 = X 3 .
The exponential map exp : g G defined on any invariant vector field X g reads as:
exp ( X ) = X t ( e ) t = 1 = X 1 ( e ) .
In the classical case of the matrix group, we recover:
exp : gl ( d , R ) G L + ( d , R ) , exp A = ( e t A ) ( t = 1 ) = e A .
Since d ( exp ) 0 = e , V G a neighborhood of e where exp is a diffeomorphism (local). This fact comes from the inverse map theorem [32].
The set Aut ( G ) of G-automorphisms is a Lie group with Lie algebra aut ( G ) when G is connected [32,33,34,35]. Furthermore, if G is also simply connected, aut ( G ) is isomorphic to the Lie algebra g of g -derivations, i.e., the set of linear maps D : g g , which satisfy the Leibniz rule D [ X , Y ] = [ D X , Y ] + [ X , D Y ] .
A remarkable relationship between a homomorphism and its derivative is given by the commutative diagram:
g d ( φ ) e The derivative at e h exp g exp h G φ The homomorphism H
As an example, we recover the very well-known formula involving the determinant and the trace of a matrix:
e t r A = det ( exp A ) , A gl ( d , R ) .
where G = G L + ( n , R ) , g = gl ( n , R ) , H = R , and h = R .
Definition A2.
A Lie algebra g is said to be:
1. 
Abelian if every bracket is null;
2. 
Nilpotent if the sequence a d j ( g ) = a d j 1 ( g ) , a d 1 ( g ) , j = 1 , 2 , stabilizes at zero, i.e.,
k 1 : a d k ( g ) = a d k 1 ( g ) , a d 1 ( g ) = 0
3. 
Solvable if the sequence a d ( j ) ( g ) = a d ( j 1 ) ( g ) , a d ( j 1 ) ( g ) j = 1 , 2 , stabilizes at zero, i.e.,
k 1 : a d ( k ) ( g ) = a d ( k 1 ) ( g ) , a d ( k 1 ) ( g ) = 0
4. 
Simple if g is not Abelian and does not contain non-proper ideals;
5. 
Semisimple if the solvable radical r ( g ) (the largest solvable Lie subalgebra of g ) is null.
We classify a Lie group with the same name that is inherited from the classification of its Lie algebra.
Example A2.
1. 
R n is Abelian;
2. 
T m = S 1 × . × S 1 is Abelian and compact;
3. 
R n × T n , n , m N is the general Abelian group shape;
4. 
The three-dimensional Heisenberg group is nilpotent;
5. 
The affine group A y 0 1 : A G L ( n , R ) , y R n is solvable;
6. 
S O ( n , R ) is compact:
S O ( n , R ) = { A G L ( n , R ) A A t = e }
7. 
S O ( 4 , R ) is compact and semisimple;
8. 
S L ( n , R ) is non-bounded and semisimple.

References

  1. Colonius, F.; Kliemann, W. The Dynamics of Control; Birkhäuser: Boston, MA, USA, 2000. [Google Scholar]
  2. Isidori, A. Nonlinear Control Systems; Springer: Cham, Switzerland, 1998. [Google Scholar] [CrossRef]
  3. Jurdjevic, V. Geometric Control Theory; Cambridge University Press: Cambridge, UK, 1997; ISBN1 0521058244. ISBN2 9780521058247. [Google Scholar]
  4. Pontryagin, L.S.; Boltyanskii, V.G.; Gamkrelidze, R.V.; Mishchenko, E.F. The Mathematical Theory of Optimal Processes; English translation; Interscience: Hoboken, NJ, USA, 1962. [Google Scholar] [CrossRef]
  5. Sachkov, Y. Control Theory on Lie Groups; Lecture Notes SISSA: Trieste, Italy, 2006. [Google Scholar] [CrossRef]
  6. Helgason, S. Differential Geometry, Lie Groups and Symmetric Spaces; Academic Press: New York, NY, USA, 1978. [Google Scholar]
  7. San Martin, L.A.B. Algebras de Lie; Editorial UNICAMP: Campinas, Brazil, 1999. [Google Scholar]
  8. Varadarajan, V. Lie Groups, Lie Algebras and their Representations; Prentice Hall, Inc.: Upper Saddle River, NJ, USA, 1974. [Google Scholar] [CrossRef]
  9. Axelby, G.S. Round-Table Discussion on the Relevance of Control Theory. Automatica 1973, 9, 279–281. [Google Scholar] [CrossRef]
  10. Lee, D.H.; Milroy, I.P.; Tyler, K. Application of Pontryagin’s Maximum Principle to the Semi-automatic Control of Rail Vehicles. In Proceedings of the Second Conference on Control Engineering 1982: Merging of Technology and Theory to Solve Industrial Automation Problems, Newcastle, Australia, 25–27 August 1982; pp. 233–236. [Google Scholar]
  11. Leitmann, G. Optimization Techniques with Application to Aerospace Systems; Academic Press Inc.: London, UK, 1962. [Google Scholar]
  12. Markus, L. Controllability of Multi-Trajectories on Lie Groups. Lect. N.a. Math. 1980. [Google Scholar] [CrossRef]
  13. Pesch, H.J.; Plail, M. The Cold War and the Maximum Principle of Optimal Control. Documenta Mathematica; Extra Volume ISMP. 2012, pp. 331–343. Available online: https://www.math.uni-bielefeld.de/documenta/vol-ismp/48_pesch-hans-josef-cold-war.pdf (accessed on 1 January 2012).
  14. Shell, K. Applications of Pontryagin’s Maximum Principle to Economics. In Mathematical Systems Theory and Economics I and II, Volume 11/12; The Series Lecture Notes in Operations Research and Mathematical Economics; Springer: Berlin/Heidelberg, Germany, 1969; pp. 241–292. [Google Scholar]
  15. Naidu, D.S. Optimal Control Systems; CRC Press: Washington D.C., USA, 2002; ISBN 0-8493-0892-5. [Google Scholar]
  16. Ledzewicz, U.; Shattler, H. Optimal Control for a Bilinear Model with Recruiting Agent in Cancer Chemotherapy. In Proceedings of the MTNS, Kyoto, Japan, 24–28 July 2006. [Google Scholar]
  17. Ledzewick, U.; Shattler, H. Optimal controls for a two compartment model for cancer chemotherapy with quadratic objective. In Proceedings of the MTNS, Kyoto, Japan, 24–28 July 2006. [Google Scholar]
  18. Sussmann, H.; Willems, C. 300 Years of Optimal Control: From the Brachystochrone to the Maximum Principle. IEEE Control. Hist. Perspect. 1997, 17, 32–44. [Google Scholar] [CrossRef] [Green Version]
  19. Ayala, V.; Tirao, J. Linear Control Systems on Lie Groups and Controllability. AMS Proc. Symp. Pure Math. 1999, 64, 47–64. [Google Scholar]
  20. Ayala, V.; Da Silva, A. On the characterization of the controllability property for linear control systems on nonnilpotent, solvable three-dimensional Lie groups. J. Differ. Equ. 2019, 266, 8233–8257. [Google Scholar] [CrossRef] [Green Version]
  21. Ayala, V.; Da Silva, A. The control set of a linear control system on the two dimensional Lie group. J. Differ. Equ. 2020, 268, 6683–6701. [Google Scholar] [CrossRef]
  22. Ayala, V.; Da Silva, A.; Jouan, P.; Zsigmond, G. Control sets of linear systems on semi-simple Lie groups. J. Differ. Equ. 2020, 269, 449–466. [Google Scholar] [CrossRef] [Green Version]
  23. Ayala, V.; Da Silva, A. A Semigroup associated to a linear control system on a Lie group. Syst. Control. Lett. 2016, 98, 33–36. [Google Scholar] [CrossRef] [Green Version]
  24. Ayala, V.; Jouan, P. Almost Riemannian Geometry on Lie Groups. Siam J. Control Optim. 2016, 54, 2919–2947. [Google Scholar] [CrossRef] [Green Version]
  25. Jouan, P. Equivalence of Control Systems with Linear Systems on Lie Groups and Homogeneous Spaces. ESAIM Control Optim. Calc. Var. 2010, 16, 956–973. [Google Scholar] [CrossRef]
  26. Jouan, P. Controllability of Linear Systems on Lie groups. J. Dyn. Control Syst. 2011, 17, 591–616. [Google Scholar] [CrossRef]
  27. Ayala, V.; Kizil, E.; Hacibekiroglu, A. Observability of General Linear Pairs. Comput. Math. Appl. 2001, 39, 35–43. [Google Scholar] [CrossRef] [Green Version]
  28. Ayala, V.; Hacibekiroglu, A.K. Observable Linear Pairs. Comput. Appl. Math. 1997, 16, 205–214. [Google Scholar]
  29. Curtis, L.M. Matrix Groups; Springer: New York, NY, USA, 1984. [Google Scholar] [CrossRef]
  30. Wonham, M.W. Linear Multivariable Control: A Geometric Approach; Applications of Mathematics; Springer: New York, NY, USA, 1979; Volume 10, p. 326. ISBN 978-1-4612-1082-5. [Google Scholar]
  31. Kalman, R. Lecture Notes on Controllability and Observability; Springer International Publishing: Cham, Switzerland, 1968; ISBN 978-3-642-11063-4. [Google Scholar]
  32. Warner, F.W. Foundations of Differential Manifolds and Lie Groups; Scott Foreman: Glenview, IL, USA, 1971. [Google Scholar]
  33. Ayala, V.; San Martín, L. Controllability Properties of a Class of Control Systems on Lie Groups; Lectures Notes in Control and Information Science; Springer: London, UK, 2001. [Google Scholar] [CrossRef]
  34. Ayala, V.; Da Silva, A. Controllability of Linear Control Systems on Lie Groups with Semisimple Finite Center. SIAM J. Control. Optim. 2017, 55, 1332–1343. [Google Scholar] [CrossRef] [Green Version]
  35. Kalman, R.; Ho, Y.; Narendra, K. Controllability of Linear Dynamical Systems. Contrib. Diff. Equ. 1962, 1, 189–213. [Google Scholar]

Share and Cite

MDPI and ACS Style

Ayala, V.; Román-Flores, H.; Torreblanca Todco, M.; Zapana, E. Observability and Symmetries of Linear Control Systems. Symmetry 2020, 12, 953. https://doi.org/10.3390/sym12060953

AMA Style

Ayala V, Román-Flores H, Torreblanca Todco M, Zapana E. Observability and Symmetries of Linear Control Systems. Symmetry. 2020; 12(6):953. https://doi.org/10.3390/sym12060953

Chicago/Turabian Style

Ayala, Víctor, Heriberto Román-Flores, María Torreblanca Todco, and Erika Zapana. 2020. "Observability and Symmetries of Linear Control Systems" Symmetry 12, no. 6: 953. https://doi.org/10.3390/sym12060953

APA Style

Ayala, V., Román-Flores, H., Torreblanca Todco, M., & Zapana, E. (2020). Observability and Symmetries of Linear Control Systems. Symmetry, 12(6), 953. https://doi.org/10.3390/sym12060953

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop