# Differentiable Families of Planar Bimodal Linear Control Systems

Authors: Josep Ferrer; M. Dolors Magret; Marta Peña
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

We consider bimodal linear control systems consisting of two subsystems acting on each side of a given hyperplane, assuming continuity along it. For a differentiable family of planar bimodal linear control systems, we obtain its stratification
diagram and, if controllability holds for each value of the parameters, we construct a differentiable family of feedbacks which
stabilizes both subsystems for each value of the parameters.


---

## Body

## 1. Introduction

Piecewise linear control systems (in particular, the bimodal ones: see, for example, [1–3]) have attracted the interest of the researchers in recent years, as a special class of switched systems (see, e.g., [4–6]), by their wide range of applications, as well as by the possible theoretical approaches, even in the planar case (see, e.g., [7]).Bimodal linear control systems (BLCS) consist of two subsystems acting on each side of a given hyperplane, assuming continuity along the separating hyperplane. These systems present a complex dynamical behaviour, even for low dimensions, as has been shown in several works. For example, in [8], it is proved that a planar bimodal linear system is stable if each subsystem is stable, but this does not hold for a bimodal linear system with three-state variables. On the other hand, since typically the number of state variables of systems describing elementary circuits is two or three (see [9]), we devote a special attention to the planar case. Here we tackle two problems concerning parameterized families of planar BLCS.Firstly, obtaining of its stratification diagram with regard to the natural equivalence relation is defined by change of basis in the state space (which preserve the hyperplanes parallel to the separating one). Previously, it is necessary to list the possible equivalence classes and to obtain a complete set of classifying invariant parameters (Theorem6).By the way, Arnold's theory allows us to restrict this study to the so-called “miniversal” deformation families. Indeed, the equivalence classes are just the orbits of a certain group action, so that they are differentiable manifolds and Arnold's machinery is applicable. Moreover we remark that by joining the orbits according to the discrete classifying invariants one obtains differentiable “strata” (each one formed by the union of classes differing only on continuous classifying invariants). We list the dimension of each orbit and the corresponding strata (Proposition7). As an application of the previous results, we present the unobservable bifurcation diagram of a miniversal deformation (Example 10).Secondly we consider parameterized families of controllable BLCS. It is known (see [3]) that for each value of the parameter there is a feedback which stabilizes the corresponding system. Thus, we lead to the quite general question of whether pointwise solvability implies the existence of a nicely parameterized solution [10]. This parameterized family of pointwise stabilizers may not be differentiable (not even continuous). Our results allow constructing a differentiable family of feedbacks which stabilizes the corresponding system for each value of the parameter (Theorem 15).We point out that when dealing with parameterized families of BLCS, the nongeneric case of unobservable ones appears in a natural way. See, for example, the circuit modeling the Fitzghugh-Nagumo equations in [9], where the unobservable case appears if R
3
=
R
4.Finally, notice that, as in previous works concerning single control systems, we will use geometrical techniques: reducing bases, stratifications, and miniversal deformations In this sense, we expect that the geometrical approach in [11] could be translated to BLCS in the future.Throughout the paper,ℝ will denote the set of real numbers and M
n
×
m
(
ℝ
) the set of matrices having n rows and m columns and entries in ℝ (in the case where n
=
m, we will simply write M
n
(
ℝ
)).


## 2. Planar Bimodal Linear Control Systems

Let us consider a bimodal linear control system (BLCS) given by(1)
x
˙
(
t
)
=
A
1
x
(
t
)
+
B
1
u
(
t
)
,
y
(
t
)
=
C
x
(
t
)
,
if
y
(
t
)
≤
0
,
x
˙
(
t
)
=
A
2
x
(
t
)
+
B
2
u
(
t
)
,
y
(
t
)
=
C
x
(
t
)
,
if
y
(
t
)
≥
0
,
where A
1
,
A
2
∈
M
n
(
ℝ
); B
1
,
B
2
∈
M
n
×
1
(
ℝ
); C
∈
M
1
×
n
(
ℝ
). One assumes that the dynamics is continuous along the separating hyperplane H
=
{
x
∈
ℝ
n
:
C
x
=
0
}; that is to say, both subsystems coincide for y
(
t
)
=
0.By means of a linear change in the state variablex
(
t
), one can consider C
=
(
1
0
⋯
0
)
∈
M
1
×
n
(
ℝ
). Hence H
=
{
x
∈
ℝ
n
:
x
1
=
0
} and continuity along H is equivalent to
(2)
B
2
=
B
1
,
A
2
e
i
=
A
1
e
i
,
2
≤
i
≤
n
.
We will write from now on B
=
B
1
=
B
2.Definition 1.
In the above conditions, one says that the triple of matrices(
A
1
,
A
2
,
B
) defines a bimodal linear control system (BLCS). Throughout the paper, 𝒳 will denote the set of these triples:
(3)
𝒳
=
{
(
A
1
,
A
2
,
B
)
∈
M
n
(
ℝ
)
×
M
n
(
ℝ
)
×
M
n
×
1
(
ℝ
)
|
A
2
e
i
=
A
1
e
i
,
2
≤
i
≤
n
}
which is obviously a (
n
2
+
2
n
)-differentiable manifold.
The system is calledobservable if
(4)
rank
⁡
(
C
C
A
i
⋯
C
A
i
n
-
1
)
=
n
,
i
=
1,2
.A natural goal is simplifying the matricesA
1
,
A
2, and B by means of changes in the variables x
(
t
) which preserve the qualitative behavior of the system. So, one considers linear changes in the state variables space preserving the hyperplanes x
1
(
t
)
=
k.Definition 2.
One callsadmissible basis changes those given by the matrices
(5)
𝒮
:
=
{
S
∈
G
l
n
(
ℝ
)
|
S
=
(
1
0
U
T
)
,
T
∈
G
l
n
-
1
(
ℝ
)
,
U
∈
M
n
×
1
(
ℝ
)
(
1
0
U
T
)
}
.
Then,(
A
1
,
A
2
,
B
)
,
(
A
1
′
,
A
2
′
,
B
′
)
∈
𝒳 are said to be equivalent if there exists a matrix S
∈
𝒮 (representing an admissible basis change) such that (
A
1
′
,
A
2
′
,
B
′
)
=
(
S
-
1
A
1
S
,
S
-
1
A
2
S
,
S
-
1
B
).Notice that the matrixC is not involved in this definition since C
S
=
C for any S
∈
𝒮.When considering canonical forms, it is necessary that the coefficients appearing in them as well as the conditions used to distinguish the different types do not depend on the admissible basis which one considers; that is to say, they are preserved under admissible basis changesS
∈
𝒮. It is wellknown that tr
⁡
⁡
A
1, tr
⁡
⁡
A
2, det
⁡
A
1, and det
⁡
A
2 are invariant under any basis change S
∈
G
l
n
(
ℝ
). We focus on the additional invariants when only admissible basis changes S
∈
𝒮 are considered.Definition 3.
A real number (resp., a property) associated with a triple(
A
1
,
A
2
,
B
) is called 𝒮
-invariant if it is preserved by admissible basis changes; that is to say, it has the same value (resp., it is also true) for any other triple (
A
1
′
,
A
2
′
,
B
′
)
𝒮-equivalent to the given one.For example, it is obvious that they areS-invariant: the top coefficient b
1 in B, the matrix C, and the condition of (
A
1
,
A
2
,
B
) being observable.We introduce another𝒮-invariant that will be used under additional hypotheses.Definition 4.
Given a triple(6)
A
1
=
(
a
1
a
3
a
2
a
4
)
,
A
2
=
(
γ
1
a
3
γ
2
a
4
)
,
B
=
(
b
1
b
2
)
one writes(7)
Δ
0
=
det
⁡
(
a
3
b
1
a
4
b
2
)
=
a
3
b
2
-
a
4
b
1
,
Δ
12
=
a
2
(
a
4
-
γ
1
)
-
γ
2
(
a
4
-
a
1
)
,
Δ
1
=
b
1
a
2
+
(
a
4
-
a
1
)
b
2
,
Δ
2
=
b
1
γ
2
+
(
a
4
-
γ
1
)
b
2
.Lemma 5.
The above triple is unobservable if and only ifa
3
=
0. In this case one has (1)
det
⁡
(
(
a
1
-
γ
1
a
2
-
γ
2
)
|
A
i
(
a
1
-
γ
1
a
2
-
γ
2
)
)
=
(
a
1
-
γ
1
)
Δ
12, i
=
1,2,
det
⁡
(
B
|
A
i
B
)
=
b
1
Δ
i, i
=
1,2;
(2)
the action ofS
∈
𝒮 transforms Δ
1, Δ
2, and Δ
12, respectively, into
(8)
1
det
⁡
S
Δ
1
,
1
det
⁡
S
Δ
2
,
1
det
⁡
S
Δ
12
.
In particular, it is𝒮-invariant the sign (positive, negative, or zero):
(9)
sign
⁡
(
Δ
1
Δ
2
)
.Proof.
Clearly,(
*
)
(
C
C
A
1
)
=
(
1
0
a
1
a
3
)
,
(
C
C
A
2
)
=
(
1
0
γ
1
a
3
)
do not have maximal rank when a
3
=
0. Then(1)
it is a straightforward computation;
(2)
ifa
3
=
0, then a
1 and γ
1 are eigenvalues of A
1 and A
2.
The action ofS transforms the matrices in (
*
) into their left product by S
-
1.Theorem 6.
With the above notation:(1)
Table1 summarizes some 𝒮-invariant numbers and properties, as well as the hypotheses for each one;
(2)
Table2 lists the possible canonical forms and the classification criteria.Table 1
𝒮-Invariant numbers.
Hypotheses
Numbers
Properties
(1)
Δ
0
,
b
1
a
3
=
0
(2)
b
1
=
0
b
2
=
0
(3)
a
3
=
0
a
1
,
γ
1
,
a
4
Δ
12
=
0
,
Δ
1
=
0
,
Δ
2
=
0
(3′)
a
3
=
0
,
Δ
12
≠
0
Δ
1
/
Δ
12
,
Δ
2
/
Δ
12
(4)
a
3
=
0
,
a
1
=
a
4
a
2
=
0
(4′)
a
3
=
0
,
γ
1
=
a
4
γ
2
=
0
(5)
b
1
=
0
,
a
3
=
0
,
a
1
=
a
4
b
2
/
a
2
(5′)
b
1
=
0
,
a
3
=
0
,
γ
1
=
a
4
b
2
/
γ
2Table 2
Canonical forms.
Label
Classification criteria
Canonical forms
CF1
a
3
≠
0
(
tr
⁡
A
1
1
-
det
⁡
A
1
0
)
,
(
tr
⁡
A
2
1
-
det
⁡
A
2
0
)
,
(
b
1
Δ
0
)
CF2
a
3
=
0
,
a
1
≠
a
4
,
γ
1
≠
a
4
,
Δ
12
≠
0
(
a
1
0
0
a
4
)
,
(
γ
1
0
1
a
4
)
,
(
b
1
-
Δ
1
/
Δ
12
)
CF3
a
3
=
0
,
a
1
≠
a
4
,
γ
1
≠
a
4
,
Δ
12
=
0
,
Δ
1
≠
0
(
a
1
0
0
a
4
)
,
(
γ
1
0
0
a
4
)
,
(
b
1
1
)
CF4
a
3
=
0
,
a
1
≠
a
4
,
γ
1
≠
a
4
,
Δ
12
=
0
,
Δ
1
=
0
(
a
1
0
0
a
4
)
,
(
γ
1
0
0
a
4
)
,
(
b
1
0
)
CF5
a
3
=
0
,
a
1
=
a
4
,
γ
1
≠
a
4
,
a
2
≠
0
(
a
4
0
1
a
4
)
,
(
γ
1
0
0
a
4
)
,
(
b
1
Δ
2
/
Δ
12
)
CF5′
a
3
=
0
,
a
1
≠
a
4
,
γ
1
=
a
4
,
γ
2
≠
0
(
a
1
0
0
a
4
)
,
(
a
4
0
1
a
4
)
,
(
b
1
-
Δ
1
/
Δ
12
)
CF6
a
3
=
0
,
a
1
=
a
4
,
γ
1
≠
a
4
,
a
2
=
0
,
Δ
2
≠
0
(
a
4
0
0
a
4
)
,
(
γ
1
0
0
a
4
)
,
(
b
1
1
)
CF6′
a
3
=
0
,
a
1
≠
a
4
,
γ
1
=
a
4
,
γ
2
=
0
,
Δ
1
≠
0
(
a
1
0
0
a
4
)
,
(
a
4
0
0
a
4
)
,
(
b
1
1
)
CF7
a
3
=
0
,
a
1
=
a
4
,
γ
1
≠
a
4
,
a
2
=
0
,
Δ
2
=
0
(
a
4
0
0
a
4
)
,
(
γ
1
0
0
a
4
)
,
(
b
1
0
)
CF7′
a
3
=
0
,
a
1
≠
a
4
,
γ
1
=
a
4
,
γ
2
=
0
,
Δ
1
=
0
(
a
1
0
0
a
4
)
,
(
a
4
0
0
a
4
)
,
(
b
1
0
)
CF8
a
3
=
0
,
a
1
=
a
4
=
γ
1
,
a
2
≠
0
,
γ
2
≠
0
,
b
1
≠
0
(
a
4
0
1
a
4
)
,
(
a
4
0
γ
2
/
a
2
a
4
)
,
(
b
1
0
)
CF9
a
3
=
0
,
a
1
=
a
4
=
γ
1
,
a
2
≠
0
,
γ
2
≠
0
,
b
1
=
0
(
a
4
0
1
a
4
)
,
(
a
4
0
γ
2
/
a
2
a
4
)
,
(
0
b
2
/
a
2
)
CF10
a
3
=
0
,
a
1
=
a
4
=
γ
1
,
a
2
≠
0
,
γ
2
=
0
,
b
1
≠
0
(
a
4
0
1
a
4
)
,
(
a
4
0
0
a
4
)
,
(
b
1
0
)
CF10′
a
3
=
0
,
a
1
=
a
4
=
γ
1
,
a
2
=
0
,
γ
2
≠
0
,
b
1
≠
0
(
a
4
0
0
a
4
)
,
(
a
4
0
1
a
4
)
,
(
b
1
0
)
CF11
a
3
=
0
,
a
1
=
a
4
=
γ
1
,
a
2
≠
0
,
γ
2
=
0
,
b
1
=
0
(
a
4
0
1
a
4
)
,
(
a
4
0
0
a
4
)
,
(
0
b
2
/
a
2
)
CF11′
a
3
=
0
,
a
1
=
a
4
=
γ
1
,
a
2
=
0
,
γ
2
≠
0
,
b
1
=
0
(
a
4
0
0
a
4
)
,
(
a
4
0
1
a
4
)
,
(
0
b
2
/
γ
2
)
CF12
a
3
=
0
,
a
1
=
a
4
=
γ
1
,
a
2
=
0
,
γ
2
=
0
,
b
1
≠
0
(
a
4
0
0
a
4
)
,
(
a
4
0
0
a
4
)
,
(
b
1
0
)
CF13
a
3
=
0
,
a
1
=
a
4
=
γ
1
,
a
2
=
0
,
γ
2
=
0
,
b
1
=
0
,
b
2
≠
0
(
a
4
0
0
a
4
)
,
(
a
4
0
0
a
4
)
,
(
0
1
)
CF14
a
3
=
0
,
a
1
=
a
4
=
γ
1
,
a
2
=
0
,
γ
2
=
0
,
b
1
=
0
,
b
2
=
0
(
a
4
0
0
a
4
)
,
(
a
4
0
0
a
4
)
,
(
0
0
)Proof.
(i) Concerning Table1
(1)
the𝒮-action on A
1 and B can be formulated as
(10)
S
-
1
(
A
1
,
B
)
(
S
∣
0
0
∣
1
)
=
(
1
0
u
t
)
-
1
(
a
1
a
3
b
1
a
2
a
4
b
2
)
(
1
0
0
u
t
0
0
0
1
)
=
(
*
*
|
(
1
0
u
t
)
-
1
(
a
3
b
1
a
4
b
2
)
(
t
0
0
1
)
)
.
Therefore,Δ
0 is 𝒮-invariant:
(11)
det
⁡
(
(
1
0
u
t
)
-
1
(
a
3
b
1
a
4
b
2
)
(
t
0
0
1
)
)
=
det
⁡
(
a
3
b
1
a
4
b
2
)
.
We have seen thata
3
≠
0 if and only if
(12)
rank
⁡
(
C
C
A
i
)
=
2
,
i
=
1,2
,
which is𝒮-invariant.
(2)
Ifb
1
=
0, then
(13)
S
-
1
(
0
b
2
)
=
1
t
(
t
0
-
u
1
)
(
0
b
2
)
=
1
t
(
0
b
2
)
.
(3)
Ifa
3
=
0, then a
1
,
a
4, and γ
1 are the eigenvalues of A
1 and A
2. Then
(14)
Ker
⁡
(
A
1
-
a
1
I
)
=
Ker
⁡
(
A
2
-
γ
1
I
)
if and only if(15)
rank
⁡
(
a
2
a
4
-
a
1
γ
2
a
4
-
γ
1
)
=
1
or, equivalently,(16)
0
=
det
⁡
(
a
2
a
4
-
a
1
γ
2
a
4
-
γ
1
)
=
Δ
12
.
In a similar way,Δ
1
=
0 if and only if
(17)
(
b
1
b
2
)
∈
Ker
⁡
(
A
1
-
a
1
I
)
andΔ
2
=
0 if and only if(18)
(
b
1
b
2
)
∈
Ker
⁡
(
A
2
-
γ
1
I
)
.
(3′)
This case follows from (3) and the above lemma.
(4)
Clearly, ifa
3
=
0 and a
1
=
a
4, then a
2
=
0 if and only if A
1 diagonalizes.
(4′)
Analogously than (4) forγ
2
=
0.
(5)
Returning to the formulation in (1):(19)
(
1
0
u
t
)
-
1
(
a
1
0
0
a
2
a
1
b
2
)
(
1
0
0
u
t
0
0
0
1
)
=
(
a
1
0
0
a
2
t
a
1
b
2
t
)
.
(5′)
Analogously than (5) forγ
1
=
a
4.
(ii) This case follows from [12], bearing in mind the 𝒮-invariants in (i).


## 3. Stratification Diagrams

In the previous section we have partitioned the set of BLCS into equivalence classes, characterized by the reduced forms in Theorem6. In order to study, for example, the changes when a BLCS is perturbed, some natural questions arise about the geometric structure of this equivalence partition. We will see in a moment that each equivalence class is a manifold, as well as the “strata” obtained by joining the classes that differ only in the continuous classification parameters. Their dimensions are listed in Proposition 7.Concerning perturbations, small changes in the coefficients of the matrices defining the system may give rise to nonequivalent systems. Then in order to explain the behavior of the system under small perturbations, it is necessary to know the nearby equivalence classes. The “stratification diagram” (“bifurcation diagram” in Arnold's terminology) of a parameterized family of systems is the partition of the parameter space according to the equivalence class. Arnold's theory [13] shows that they are induced by the ones of the so-called versal deformations. In particular, the miniversal deformations are the simplest parameterized families which provide all the information about which equivalence classes are near a given one, that is to say, which canonical forms appear when the given one is perturbed.The starting point is that the above equivalence classes are actually the orbits with regard to the action of the Lie group𝒮 on the differentiable manifold 𝒳:
(20)
α
:
𝒮
×
𝒳
→
𝒳
defined by
(21)
α
(
S
,
𝒳
)
=
(
S
-
1
A
1
S
,
S
-
1
A
2
S
,
S
-
1
B
)
.
Given any triple of matrices (
A
1
,
A
2
,
B
)
∈
𝒳, we will denote by 𝒪
(
A
1
,
A
2
,
B
) its orbit (or equivalence class).As an application of the closed orbit lemma (see [14]), we deduce that equivalence classes are differentiable manifolds. Namely, any equivalence class is a locally closed differentiable submanifold of 𝒳 and its boundary is a union of equivalence classes or orbits of strictly lower dimension. In particular, equivalence classes or orbits of minimal dimension are closed.Moreover, orbits having the same discrete classification parameters (but differing in the continuous ones) can be joined in a finite number of “strata” which in our case are as well differentiable manifolds (see [12]). For the commodity of the reader we adapt the results there.Proposition 7.
Table3 lists the dimensions of the orbits (i.e., the equivalence classes) and the corresponding strata (i.e., the union of the orbits of the same type when the parameters appearing in the canonical form vary).Table 3
Dimension of orbits and strata for each case in Table2.
Canonical form
Dimension of the orbit
Dimension of the stratum
CF1
2
8
CF2
2
7
CF3
2
6
CF4
1
5
CF5, CF5′
2
6
CF6, CF6′
2
5
CF7, CF7′
1
4
CF8
2
5
CF9
1
4
CF10, CF10′
2
4
CF11, CF11′
1
3
CF12
1
3
CF13
1
2
CF14
0
1To illustrate the above considerations and as an application of the results in Section2, we will present (see Figure 1) the unobservable bifurcation diagram of the miniversal deformation of a system of type CF10′. The main definitions and results about deformations and versality can be found in [13, 15]. Here we rewrite them down, adapted to our particular case.Figure 1
Stratification diagram.Definition 8.
Adeformation of (
A
1
,
A
2
,
B
)
∈
𝒳 is a differentiable map φ
:
U
→
𝒳, with U an open neighbourhood of the origin ℝ
d, such that φ
(
0
)
=
(
A
1
,
A
2
,
B
).
A deformationφ
:
U
→
𝒳 of (
A
1
,
A
2
,
B
) is called versal at 0 if, for any other deformation of (
A
1
,
A
2
,
B
), ψ
:
V
→
𝒳, there exists a neighbourhood V
′
⊆
V with 0
∈
V
′, a differentiable map γ
:
V
′
→
U with γ
(
0
)
=
0, and a deformation of the identity I
∈
𝒮, θ
:
V
′
→
𝒮, such that ψ
(
μ
)
=
α
(
θ
(
μ
)
,
φ
(
γ
(
μ
)
)
) for all μ
∈
V
′.
A versal deformation with minimal number of parametersd is called miniversal deformation.A miniversal deformation can be obtained from the normal space to the orbit with regard to some scalar product.Proposition 9 (see [16]).
We consider the following scalar product in𝒳:
(22)
〈
(
A
1
,
A
2
,
B
)
,
(
A
1
′
,
A
2
′
,
B
′
)
〉
=
tr
⁡
(
A
1
t
A
1
′
)
+
tr
⁡
(
A
2
t
A
2
′
)
+
tr
⁡
(
B
t
B
′
)
.
(i)
The normal space to the orbit of
(
A
1
,
A
2
,
B
)
at
(
A
1
,
A
2
,
B
)
N
(
A
1
,
A
2
,
B
)
𝒪
(
A
1
,
A
2
,
B
)
∩
𝒳
is the vector subspace consisting of triples
(
X
1
,
X
2
,
Y
)
∈
𝒳
such that
(23)
A
1
X
1
t
-
X
1
t
A
1
+
A
2
X
2
t
-
X
2
t
A
2
-
B
Y
t
∈
𝔸
,
where
𝔸
is the set
(24)
𝔸
=
{
M
=
(
m
i
j
)
|
m
i
j
=
0,2
≤
i
≤
n
,
1
≤
j
≤
n
}
.
(ii)
Then the mapping
(25)
ℝ
d
→
𝒳
(
η
1
,
…
,
η
d
)
→
(
A
1
,
A
2
,
B
)
+
η
1
V
1
+
⋯
+
η
d
V
d
,
where
{
V
1
,
…
,
V
d
}
is any basis of the vector space
N
(
A
1
,
A
2
,
B
)
𝒪
(
A
1
,
A
2
,
B
)
, is a miniversal deformation of
(
A
1
,
A
2
,
B
)
.Normal spaces of two equivalent triples can be obtained one from the other. Thus, it is always possible to restrict ourselves to the case where the triple is in its canonical form.Here, as an application of the previous results, we present the unobservable bifurcation diagram of a miniversal deformation: Figure1 shows the geometrical configuration of the unobservable strata near a given system of type CF10′.Example 10.
Consider a bimodal linear dynamical system of type CF10′ whose canonical form is(26)
A
1
=
(
a
4
0
0
a
4
)
,
A
2
=
(
a
4
0
1
a
4
)
,
B
=
(
b
1
0
)
.
Then, N
(
A
1
,
A
2
,
B
)
𝒪
(
A
1
,
A
2
,
B
)
∩
𝒳 is the vector subspace consisting of triples (
X
1
,
X
2
,
Y
)
∈
𝒳
(27)
X
1
=
(
x
1
x
3
x
2
x
4
)
,
X
2
=
(
x
5
x
3
x
6
x
4
)
,
Y
=
(
y
1
y
2
)
such that
(28)
x
6
=
0
a
4
x
5
+
b
1
y
2
=
0
.
Moreover, parameterx
3 must be zero to avoid observable perturbations and parameters x
4
,
y
1 give orbits in the initial stratum.
Then the unobservable perturbations in the normal space to the stratum of(
A
1
,
A
2
,
B
) are parameterized by
(29)
φ
(
x
1
,
x
2
,
x
5
)
=
(
(
a
4
+
x
1
0
x
2
a
4
)
,
(
a
4
+
x
5
0
1
a
4
)
,
(
b
1
-
a
4
b
1
x
5
)
)
.
We denote byE
i the set of all triples of matrices having canonical form of type (CFi), i
=
1
,
…
,
14.
Clearly, if onlyx
1 (resp., x
2) is nonzero, it lies in E
5
′ (resp., E
8). But for only x
5, the strata E
6 and E
7 are possible in principle, depending on the value of Δ
2. In our case
(30)
Δ
2
=
b
1
γ
2
+
(
a
4
-
γ
1
)
b
2
=
b
1
+
(
-
x
5
)
(
-
a
4
b
1
x
5
)
=
1
b
1
(
b
1
2
+
a
4
x
5
2
)
.
Hence, it belongs toE
7 for x
5
2
=
-
b
1
2
/
a
4, and to E
6 otherwise.
In a similar way, ifx
1
,
x
5
≠
0 only E
2, E
3, and E
4 are possible. We have Δ
0
=
-
x
2
x
5
+
x
1. Hence, x
2
=
0 implies Δ
0
≠
0, which corresponds to E
2. If x
2
≠
0, it gives again E
2 except on the hyperbolic paraboloid x
1
=
x
2
x
5. When it happens,
(31)
Δ
1
=
b
1
x
2
+
x
1
a
4
b
1
x
5
=
x
2
b
1
(
b
1
2
+
a
4
x
5
2
)
.
Hence, it lies inE
4 for x
5
2
=
-
b
1
2
/
a
4, and in E
3 otherwise.
Finally, it is straightforward that one obtainsE
5 for x
1
=
0, x
2
,
x
5
≠
0, and E
5
′ for x
5
=
0, x
1
,
x
2
≠
0. In summary (see Figure 1),(i)
ifx
2
,
x
5
=
0, x
1
≠
0, then φ
(
x
1
,
x
2
,
x
5
)
∈
E
5
′;
(ii)
ifx
1
,
x
5
=
0, x
2
≠
0, then φ
(
x
1
,
x
2
,
x
5
)
∈
E
8;
(iii)
ifx
1
,
x
2
=
0, x
5
2
=
-
b
1
2
/
a
4, then φ
(
x
1
,
x
2
,
x
5
)
∈
E
7;
(iv)
ifx
1
,
x
2
=
0, x
5
≠
0, x
5
2
≠
-
b
1
2
/
a
4, then φ
(
x
1
,
x
2
,
x
5
)
∈
E
6;
(v)
ifx
5
=
0, x
1
,
x
2
≠
0, then φ
(
x
1
,
x
2
,
x
5
)
∈
E
5
′;
(vi)
ifx
2
=
0, x
1
,
x
5
≠
0, then φ
(
x
1
,
x
2
,
x
5
)
∈
E
2;
(vii)
ifx
1
=
0, x
2
,
x
5
≠
0, then φ
(
x
1
,
x
2
,
x
5
)
∈
E
5;
(viii)
ifx
1
,
x
2
,
x
5
≠
0, x
1
=
x
2
x
5, x
5
2
=
-
b
1
2
/
a
4, then φ
(
x
1
,
x
2
,
x
5
)
∈
E
4;
(ix)
ifx
1
,
x
2
,
x
5
≠
0, x
1
=
x
2
x
5, x
5
2
≠
-
b
1
2
/
a
4, then φ
(
x
1
,
x
2
,
x
5
)
∈
E
3;
(x)
ifx
1
,
x
2
,
x
5
≠
0, x
1
≠
x
2
x
5, then φ
(
x
1
,
x
2
,
x
5
)
∈
E
2.


## 4. Controllability and Families of Stabilizers

We have seen that in a differentiable family of BLCS different equivalence classes can appear for different values of the parameters. Let us see that, however, some global treatments are possible. Indeed, we will prove that a differentiable family of stabilizers exists if each BLCS in the given family is controllable.The notion of controllability of a single system is extended to bimodal ones in a natural way.Definition 11.
A BLCS is(completely) controllable if for any pair of states (
x
0
,
x
f
) there exists a locally integrable input u such that the solution x
x
0
,
u passes through x
f; that is, x
x
0
,
u
(
T
)
=
x
f for some T
>
0.A well-known remarkable fact is that a single linear systemx
˙
=
A
x
+
B
u is controllable if and only if its “controllability matrix” (
B
A
B
⋯
A
n
-
1
B
) has maximal rank. For planar BLCS we recall the characterization of controllability of planar BLCS obtained in [1] for observable systems and generalized in [17] to unobservable ones.Proposition 12.
Let one consider a planar BLCS defined by(
A
1
,
A
2
,
B
). One writes C
1, C
2 the controllability matrices of both subsystems
(32)
C
1
=
(
B
A
1
B
)
,
C
2
=
(
B
A
2
B
)
.
Then, it is controllable if and only if
(33)
det
⁡
C
1
det
⁡
C
2
>
0
.Remark 13.
(1) Notice that, in particular, both subsystems must be controllable, but it is not a sufficient condition.
(2) Whereas for single systems the subset of controllable ones is open and dense, the above proposition shows that it is not for BLCS systems: controllability is an open, but not generic, property.If the control function is a so-called “feedback” of the typeu
(
t
)
=
f
(
x
(
t
)
), one obtains a new dynamical system (“in closed loop”). For single linear systems x
˙
=
A
x
+
B
u, a feedback u
=
F
x gives x
˙
=
(
A
+
B
F
)
x. A remarkable fact is that it is stable for some suitable F, provided that the initial control system is controllable.As a natural generalization, in [3] any controllable BLCS is proved to be feedback stabilizable. Hence, if a differentiable parameterized family (
A
1
(
s
)
,
A
2
(
s
)
,
B
(
s
)
) is pointwise controllable (observable or not), then it is also pointwise stabilizable; that is to say, for any s
∈
ℝ there is a common feedback F
(
s
) such that both closed-loop systems A
1
(
s
)
+
B
(
s
)
F
(
s
), A
2
(
s
)
+
B
(
s
)
F
(
s
) are stable. However, the family F
(
s
) may not be differentiable (not even continuous). Here we prove that differentiable families of stabilizer feedbacks exist for n
=
2.As we have pointed out in the Introduction, the unobservable case appears generically in parameterized families of bimodal systems. A typical case is considered in the following example. As an application of the above proposition, we characterize when this family is pointwise controllable.Example 14.
Let us consider the parameterized family of planar BLCS(34)
A
1
(
s
)
=
(
a
1
s
a
2
a
4
)
,
A
2
(
s
)
=
(
γ
1
s
γ
2
a
4
)
,
B
=
(
b
1
b
2
)
,
where s
∈
ℝ. Obviously, the systems defined by these matrices are observable except for s
=
0. Let us see that the family is pointwise controllable (i.e., for any s
∈
ℝ the corresponding system is controllable) if and only if b
1
≠
0 and  (i)
a
2
γ
2
>
0, if b
2
=
0,
(ii)
det
⁡
(
a
1
b
1
a
2
b
2
)
=
det
⁡
(
γ
1
b
1
γ
2
b
2
), otherwise.From Proposition12, for any s
∈
ℝ (including the case where s
=
0) the corresponding system is controllable if and only if
(35)
(
b
1
2
a
2
+
b
1
b
2
a
4
-
b
1
b
2
a
1
-
b
2
2
s
)
×
(
b
1
2
γ
2
+
b
1
b
2
a
4
-
b
1
b
2
γ
1
-
b
2
2
s
)
>
0
.
In particular b
1
≠
0 (it suffices to take s
=
0).Ifb
2
=
0, the above inequality is
(36)
(
b
1
2
a
2
)
(
b
1
2
γ
2
)
>
0
that is to say,
(37)
a
2
γ
2
>
0
.
Assume now b
2
≠
0. In general, two polynomials of degree 1 have the same sign at any point if and only if they have the same root and the slopes have the same sign. In our case both slopes are -
b
2
2, so that the above inequality holds if and only if
(38)
b
1
2
a
2
+
b
1
b
2
a
4
-
b
1
b
2
a
1
b
2
2
=
b
1
2
γ
2
+
b
1
b
2
a
4
-
b
1
b
2
γ
1
b
2
2
which is equivalent (recall b
1
≠
0) to
(39)
b
1
a
2
-
b
2
a
1
=
b
1
γ
2
-
b
2
γ
1
.Finally, we prove the existence of differentiable families of stabilizers for differentiable families of planar controllable bimodal systems.Theorem 15.
Let(40)
(
A
1
(
s
)
,
A
2
(
s
)
,
B
(
s
)
)
,
s
∈
ℝ
be a differentiable family of planar BLCS. If it is pointwise controllable, then there is a differentiable family of feedbacks F
(
s
), s
∈
ℝ, such that
(41)
A
1
(
s
)
+
B
(
s
)
F
(
s
)
,
A
2
(
s
)
+
B
(
s
)
F
(
s
)
are stable for any s
∈
ℝ.
More explicitly, if(42)
A
1
(
s
)
=
(
a
1
a
3
a
2
a
4
)
,
A
2
(
s
)
=
(
γ
1
a
3
γ
2
a
4
)
,
B
(
s
)
=
(
b
1
b
2
)
,
C
1
=
(
B
(
s
)
A
1
(
s
)
B
(
s
)
)
,
C
2
=
(
B
(
s
)
A
2
(
s
)
B
(
s
)
)
,
where all the coefficients are assumed to be differentiably depending on s
∈
ℝ, one can take
(43)
F
=
(
f
1
f
2
)
=
(
x
y
)
1
det
⁡
C
1
(
det
⁡
(
b
1
a
1
b
2
a
2
)
det
⁡
(
b
1
a
3
b
2
a
4
)
-
b
2
-
b
1
)
with
(44)
x
<
-
tr
⁡
A
1
,
-
tr
⁡
A
2
,
y
<
det
⁡
A
1
,
y
<
-
det
⁡
(
b
1
γ
1
-
a
1
b
2
γ
2
-
a
2
)
det
⁡
(
b
1
a
3
b
2
a
4
)
1
det
⁡
C
2
x
+
det
⁡
A
2
det
⁡
C
1
det
⁡
C
2
.Proof.
By hypothesis, we assume(45)
det
⁡
C
1
det
⁡
C
2
>
0
for any s
∈
ℝ.
We look forF
(
s
)
=
(
f
1
f
2
) where again we assume the coefficients depending on s
∈
ℝ, such that the eigenvalues of
(46)
(
a
1
+
b
1
f
1
a
3
+
b
1
f
2
a
2
+
b
2
f
1
a
4
+
b
2
f
2
)
,
(
γ
1
+
b
1
f
1
a
3
+
b
1
f
2
γ
2
+
b
2
f
1
a
4
+
b
2
f
2
)
have negative real part for any s
∈
ℝ or, equivalently, the matrices have negative trace and positive determinant; that is to say,
(47)
b
1
f
1
+
b
2
f
2
<
-
a
1
-
a
4
,
b
1
f
1
+
b
2
f
2
<
-
γ
1
-
a
4
,
f
1
(
a
3
b
2
-
a
4
b
1
)
+
f
2
(
a
2
b
1
-
a
1
b
2
)
<
a
1
a
4
-
a
2
a
3
,
f
1
(
a
3
b
2
-
a
4
b
1
)
+
f
2
(
γ
2
b
1
-
γ
1
b
2
)
<
γ
1
a
4
-
γ
2
a
3
.
We change the variables (
f
1
,
f
2
) by (
x
,
y
) defined by
(48)
x
=
b
1
f
1
+
b
2
f
2
,
y
=
(
b
2
a
3
-
b
1
a
4
)
f
1
+
(
b
1
a
2
-
b
2
a
1
)
f
2
,
which is a change of variables, because (by hypothesis)
(49)
det
⁡
(
b
1
b
2
a
3
-
b
1
a
4
b
2
b
1
a
2
-
b
2
a
1
)
=
det
⁡
C
1
≠
0
.
Then
(50)
(
f
1
f
2
)
=
(
x
y
)
1
det
⁡
C
1
(
det
⁡
(
b
1
a
1
b
2
a
2
)
det
⁡
(
b
1
a
3
b
2
a
4
)
-
b
2
-
b
1
)
.
With this change of variables, the desired inequalities become
(51)
x
<
-
a
1
-
a
4
=
-
tr
⁡
A
1
,
x
<
-
γ
1
-
a
4
=
-
tr
⁡
A
2
,
y
<
a
1
a
4
-
a
2
a
3
=
det
⁡
A
1
,
(
a
3
b
2
-
a
4
b
1
)
(
b
1
a
2
-
b
2
a
1
)
x
-
b
2
y
det
⁡
C
1
-
(
γ
2
b
1
-
γ
1
b
2
)
(
b
2
a
3
-
b
1
a
4
)
x
-
b
1
y
det
⁡
C
1
<
γ
1
a
4
-
γ
2
a
3
.
It is straightforward that the last inequality can be rewritten:
(52)
det
⁡
(
b
1
γ
1
-
a
1
b
2
γ
2
-
a
2
)
det
⁡
(
b
1
a
3
b
2
a
4
)
1
det
⁡
C
1
x
+
y
det
⁡
C
2
det
⁡
C
1
<
det
⁡
A
2
.Example 16.
For the family in Example14, when b
1
≠
0, b
2
=
0, a
2
γ
2
>
0, differentiable families of feedbacks are given by
(53)
(
f
1
f
2
)
=
(
x
y
)
1
b
1
a
2
(
a
2
-
a
4
0
-
1
)
,
x
<
-
a
1
-
a
4
,
-
γ
1
-
a
4
,
y
<
a
1
a
4
-
a
2
s
,
y
<
(
a
2
γ
2
-
1
)
a
4
x
+
a
2
γ
2
γ
1
a
4
-
a
2
s
.
For example,
(54)
x
=
min
⁡
{
-
a
1
-
a
4
,
-
γ
1
-
a
4
}
-
1
≡
α
,
y
=
β
-
a
2
2
4
s
2
,
β
=
min
⁡
{
a
1
a
4
,
(
a
2
γ
2
-
1
)
a
4
α
+
a
2
γ
2
γ
1
a
4
}
-
1
.
That is,
(55)
F
(
s
)
=
(
α
b
1
-
a
4
α
b
1
a
2
-
β
b
1
a
2
+
a
2
4
b
1
s
2
)
.


## 5. Conclusion

In this work we consider planar bimodal linear control systems (BLCS) consisting of two subsystems acting on each side of a given hyperplane, assuming continuity along it. The set of BLCS is partitioned into equivalence classes by reducing each triple of matrices by means of a suitable change of basis. For a differentiable family of such systems (for example, perturbations of a given one) we study its stratification diagram, that is to say, the different equivalence types appearing for different values of the parameters. On the other hand, in spite of these different classes (even nonobservable ones), if pointwise controllability holds, we construct a differentiable family of feedbacks which stabilizes both subsystems for each value of the parameters.Some extensions of this work could be the application of the same techniques to tackle the case of piecewise linear control systems composed of a different partition of the state space, for example, the ones composed of three regions, with the same subsystem acting on the outer ones (see, for example, [9]). Another possible work would be the extension to bimodal linear systems with three-state variables, starting from our results in [17].



---

*Source: 292813-2014-01-22.xml*