International Journal of Control Science and Engineering
p-ISSN: 2168-4952 e-ISSN: 2168-4960
2012; 2(5): 127-135
doi: 10.5923/j.control.20120205.05
Hicham Hihi
Electrical Engineering Department, LGECOS, Cadi Ayyad University, ENSA, Av Abdelkrim khattabi BP, 57540000, Marrakech, Morocco
Correspondence to: Hicham Hihi , Electrical Engineering Department, LGECOS, Cadi Ayyad University, ENSA, Av Abdelkrim khattabi BP, 57540000, Marrakech, Morocco.
| Email: | ![]() |
Copyright © 2012 Scientific & Academic Publishing. All Rights Reserved.
In this paper, a new methodology for analysis of structural observability of controlled switching linear systems modelled by bond graphs is proposed. Causal manipulations on the bond graph model enable to determine graphically the observable subspace. A novel definition of observability is proposed. Finally, two sufficient conditions of observability are derived. The proposed method, based on a bond graph theoretic approach, assumes only the knowledge of the systems structure. These conditions can be implemented by classical bond graph theory algorithms based on finding particular paths and cycles in a bond graph.
Keywords: Hybrid Systems, Switching Systems, Bond Graph, Structural Observability
Cite this paper: Hicham Hihi , "Structural Observability of Controlled Switching Linear Systems", International Journal of Control Science and Engineering, Vol. 2 No. 5, 2012, pp. 127-135. doi: 10.5923/j.control.20120205.05.
![]() | (1) |
is the state variable,
is the input variable, and
is the output variable.
is a piecewise constant switching function and
the hybrid state.If we consider this system in a particular mode i, the equation (1) can be written as :![]() | (2) |
,
,
,
and q the number of mode.Remark 1 System (2) can be considered as a linear time invariant system (LTI).Assumptions 11) We suppose that
,
and
matrices are constant on
,where
, and constant
is an arbitrarily small and independent of mode
. For instance, suppose that the dynamics in (1) are given by
over the finite time interval
. At time
the dynamic in interval
is given by
. 2) We assume that the state vector
does not jump discontinuously at
.Under these assumptions, the observability of (1) was defined:Definition 1[22] The CSLS (1) is observable, if there exist a time
and a switching path
, such that initial state
can be determined by the corresponding system input
and the system output
,
.![]() | (3) |
[9] of system (1) is given by equation (4),![]() | (4) |
is the observability matrix of the
mode.Theorem 1[9] The CSLS (1) is observable, if
.Remark 2 From this theorem, we can deduce that:1) The system (1) can be observable, if there is only one observable sub-system (mode). 2) However, it is possible that no sub-system is observable but that the system (1) is observable.
,
and![]() | (5) |
.This theorem is a geometric criterion, thus, it is easy to transform it into algebraic form.Definition 2The joint observability matrix of system (1) is defined as :![]() | (6) |
.Theorem 3[9]System (1) is observable, if and only if
.Proof.For
, we have :
For
, we have :
In a similar way one finds :
However
and
, then
Of another share, we have :


With
,
,and
Then
.We exposed algebraic and geometric criteria of analysis of the properties of observability of CSLS. The next section is devoted to the graphic interpretation of these results by using the Bond graph approach.
on I elements and charges
on C elements), and the complementary state vector z(t)is composed of power variables (the efforts eon Celements and flows f on I elements); -
and
represent the variables going out of and into the R field; - the vector u(t)is composed of the sources; -
is composed of the zero valued variables imposed by the switches in this configuration; -
is composed of the complementary variables in the switches; - the vector y(t) is composed of the continuous outputs.![]() | Figure 1. Structure junction |
![]() | (7) |
.
is a positive matrix. Let assume that
is an invertible positive matrix. Then the second row leads to
.The third line of (7) gives :![]() | (8) |
![]() | (9) |
![]() | (10) |
.![]() | (11) |
,
,
,
,
and
.Therefore, for N switches and by considering
, we have
modes :![]() | (12) |
: the number of elements in integral causality in BGDi. i indicate the mode i.-
: the number of elements remaining in integral causality in BGDi, when a dualization of the maximum number of continuous output y and discrete output
is applied (in order to eliminate elements in integral causalities).To study structural observability of CSLS modelled by bond graph, graphical methods are proposed in the form of two sufficient conditions. In fact, formal representation of observability subspace is given for bond graph models. It is calculated through causal manipulations. The base of this subspace is used to propose a procedure to study the system observability.
.- BG-rank
. With
,
.Property 1 BG-rank
=rank
.Example 1 We consider the following acausal bond graph model. Shown in Figure 2: ![]() | Figure 2. Acausal bond graph model |
on
,
on
(
). The dimension of the system is
. We have one switch, then the number of possible configurations is
. The bond graph models in integral causality for these two configurations (modes) are given by Figure 3.![]() | Figure 3a. Bond graph model in integral causality (mode 1) : :![]() |
![]() | Figure 3b. Bond graph model in integral causality (mode 2) : : ![]() |
.The application of the derivative causality, for example on mode 1 (Figure 3.a), give the following BGD1 (Figure 4).![]() | Figure 4. BGD1 + dualization of detectors |
, this mode is observable by continuous outputs Df1 and Df2 and discrete output
, then the system is structurally observable. This result can be verified using formal calculation on the bond graph model in integral causality[3].To study the observability of system (12), it is necessary to apply this result to all modes; if one observable mode exists, the procedure is stopped. The case where no mode is observable, but when the system is observable, can be verified by formal calculation of combined matrix (4). This calculation can be formally effected :- by using the bond graph model in integral causality[3], or - by calculating the observability subspace from bond graph model in derivative causality.We chose to translate the latter in the form of a second sufficient condition. For that, formal representation of structural observability subspace, denoted as
, is given for BG model. It is calculated using causal manipulations. The base of this subspace is used to propose a procedure to study the structural observability of system.
elements remaining in integral causality and
elements in derivative causality.
algebraic equations can be written (Equation (13)) :![]() | (13) |
is either an effort variable
for
-element in integral causality or a flow variable
for
-element in integral causality,-
is either an effort variable
for
-element in derivative causality or a flow variable
for
- element in derivative causality,-
is the gain of the causal path between the 
or
-elements in integral causality and the 
or
-elements in derivative causality.Let us consider the
column vectors
whose components are the coefficients of the variables
and
in equation (13).Property 2The
column vectors
are orthogonal to the structural observability subspace vectors of the
mode. We write
and
.Procedure 1Calculation of
1) On the BGDi, dualize the maximum number of output detectors in order to eliminate the elements in integral causality,2) For each element in integral causality, write the algebraic relations with elements in derivative causality (equation 13),3) Write a column vector
for each algebraic relation with the causal path gains (Equation (13)).![]() | Figure 5a. Bond graph model in derivative causality (mode 1) |
![]() | Figure 5b. Bond graph model in derivative causality (mode 2) |
is in integral causality, so wecan write
. The coefficients of the algebraic relation are multiplied by the inductance or capacitance parameter, because of the form of the output matrix in the state equation. Thus we obtain
.In order to calculate a
basis,it is enough to find
independent row vectors
. These vectors are gathered in the matrix given by
.In the same manner, from the BGDi (and dualization of output detectors)
algebraic relations can be written (14).![]() | (14) |
is either a flow variable
for
-element in derivative causality or an effort variable
for
- element in derivative causality,-
is either a flow variable
for
-element in integral causality or a flow variable
for
- element in integral causality,-
is the gain of the causal path between the
element in derivative causality and the
element in integral causality. Suppose now the
row vectors
whose components are the coefficients of the variables
and
in equation (14).Procedure 2 Calculation of
1) On the BGDi, dualize the maximum number of output detectors in order to eliminate the elements in integral causality.2) For each element remaining in derivative causality, write the algebraic relation with elements in integral causality, (Equation (14)).3) Write a row vector
for each algebraic relation with the different gains of the causal paths, (Equation (14)).Property 3The
row vectors
compose a basis for the structural observability subspace of
mode.With
and
.Example 3 We implement procedure 2on the previous example. For mode 1, the two dynamic elements
and
are not causally connected with
, we can write
, the corresponding vectors are
and
. The algebraic equations corresponding to the elements
,
and
are given by : 

,
and
. Thus, we have
and
.Some calculation is carried out for mode 2. We obtain
, with
,
a, and
.The graphical calculation of structural observability subspaces and remark 2 lead to theorem 5.Theorem 5 If
, the CSLS system is structurally observable.Proof. We have shown for a given mode that the bond graph model in derivative causality is characterized by an algebraic equation of the form (14) from which, we construct a
base of structural observability subspace of the ith mode, denoted
.After commutation from ith mode to (i+1)th mode, implement a derivative causality on the bond graph model and dualization the maximum number of continuous and discrete outputs. We can write another algebraic relation (equation 15).![]() | (15) |
.However, the condition :
is sufficient for observability of the system, which implies that condition
is also sufficient.Example 4 Theorem 5 is now applied to the previous bond graph model, we have
, then, this system is structurally controllable.![]() | Figure 6. The acausal BG |
![]() | Figure 7. The BGIi, a) mode 1, b) mode 2, c) mode 3, d) mode 4 |
![]() | Figure 8. BGDi + dualization of detectors |
and
stays in integral causality. After dualization of the discrete output
associated to
, only one element
remaining in integral causality (figure 8.A), So this mode is not observable.In the same way, the other modes are not observable, therefore, step 1 is not verified.Step 2 : Verification of sufficient condition 2▪ Calculation of W1 (mode 1)For mode 1, the element
is in integral causality, we have one algebraic relation can be written
. The coefficients of the algebraic relation are multiplied by the inductance or capacitance parameter, because of the form of the output matrix in the state equation. Then
.The algebraic equations corresponding to the elements
and
are given by:
,
, then
,
. The dynamical element C is not causally connected with
, we can write
. The corresponding vector is
. Thus, we have
and
.Same calculation is carried out for the three other modes. We obtain
,
and
, with
,
,
,
,
,
and
.We have 

The system is structurally observable.Remark 3 If these conditions are not checked, it is necessary to use a necessary and sufficient condition. This result will be done in a future work.