Next Article in Journal
Automatic Superpixel-Based Clustering for Color Image Segmentation Using q-Generalized Pareto Distribution under Linear Normalization and Hunger Games Search
Next Article in Special Issue
Homogenization of Smoluchowski Equations in Thin Heterogeneous Porous Domains
Previous Article in Journal
Novel Design Framework for Dual-Band Frequency Selective Surfaces Using Multi-Variant Differential Evolution
Previous Article in Special Issue
Characterization of Traveling Waves Solutions to an Heterogeneous Diffusion Coupled System with Weak Advection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Robust Stabilization and Observer-Based Stabilization for a Class of Singularly Perturbed Bilinear Systems

1
Department of Computer Science and Information Engineering, Southern Taiwan University of Science and Technology, 1, Nan-Tai St., Yongkang District, Tainan 71005, Taiwan
2
Department of Electrical Engineering, Southern Taiwan University of Science and Technology, 1, Nan-Tai St., Yongkang District, Tainan 71005, Taiwan
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(19), 2380; https://doi.org/10.3390/math9192380
Submission received: 16 August 2021 / Revised: 19 September 2021 / Accepted: 22 September 2021 / Published: 25 September 2021
(This article belongs to the Special Issue Asymptotic Analysis and Homogenization of PDEs)

Abstract

:
An infinite-bound stabilization of a system modeled as singularly perturbed bilinear systems is examined. First, we present a Lyapunov equation approach for the stabilization of singularly perturbed bilinear systems for all ε∈(0, ∞). The method is based on the Lyapunov stability theorem. The state feedback constant gain can be determined from the admissible region of the convex polygon. Secondly, we extend this technique to study the observer and observer-based controller of singularly perturbed bilinear systems for all ε∈(0, ∞). Concerning this problem, there are two different methods to design the observer and observer-based controller: one is that the estimator gain can be calculated with known bounded input, the other is that the input gain can be calculated with known observer gain. The main advantage of this approach is that we can preserve the characteristic of the composite controller, i.e., the whole dimensional process can be separated into two subsystems. Moreover, the presented stabilization design ensures the stability for all ε∈(0, ∞). A numeral example is given to compare the new ε-bound with that of previous literature.

1. Introduction

As is well known, a large class of physical systems in biology, engineering, ecology, socio-economics, nuclear, thermal, and chemical processes can be properly described by bilinear models [1,2,3,4]. The above-mentioned systems are usually described by higher-order dimensional types. Due to the high-order dimensions of these factories, a large amount of computer memory and considerable operating time are required to process such a large-scale system. The singular perturbation method [5,6,7,8,9,10,11,12,13,14,15,16] is very useful for solving the problems of high dimensionality and stiffness when dealing with large-scale systems. Although the singularly perturbed systems have so many applications and combinations with other methodologies, apparently many basic issues remain to further address. Especially, stability analysis and controller design are fundamental issues for singularly perturbed systems. It is well known that stabilizing the slow and fast subsystems will not guarantee the stability of the resulting closed-loop system for all singular parameters. Therefore, it is of great importance to find the singular parameter bound for the stability of reduced and closed-loop systems. Many papers have tried various ways to calculate the bound, but only sufficiently small bound for ε are available and some of them are rather conservative. The controller design of this type of model is costly to implement, and sometimes even impossible to implement.
However, thus far, there are few literatures on the stabilization of singularly perturbed bilinear systems [17,18,19,20]. Literature [17] uses Lyapunov’s theory to study the stabilization of singularly perturbed strictly bilinear systems. Asamoah and Jamshidi [18] extended the results of [19] to singularly perturbed bilinear systems, but some of the results in [17] were proven wrong by [19]. These papers [17,18,19] only guarantee the stability for ε∈(0, 1). The above results [17,18,19] seem to be limited to continuous time. For discrete system, Ref. [20] presented two kinds of robust controllers for stabilizing singularly perturbed discrete bilinear systems. That is, all stable solutions can be implemented under the assumption that all states are available for measurement. This assumption is usually not true in practice, that is, the state of the system may not be able to achieve state feedback. Therefore, it is usually necessary to design an observer and a controller, which leads to the asymptotic stability of a given equilibrium point [21].
In this paper, an infinite ε-bound stabilization of a system modeled as singularly perturbed bilinear systems is examined. First, we present a Lyapunov equation for the stabilization of singularly perturbed bilinear systems for all ε∈(0, ∞). The method is based on the Lyapunov stability theorem. The state feedback constant gain can be determined from the admissible region of the convex polygon. Second, we extend this technique to study the observer and observer-based controller of singularly perturbed bilinear systems for all ε∈(0, ∞). Where two different methods are proposed to design the observer and observer-based control gains, one is the estimator’s gain can be calculated with known bounded input, and the other is the state feedback gain can be calculated with known observer’s gain. The main advantage of the developed schemes is that the whole dimensional design can be directly separated into two lower dimensions without the decomposition and composition of the slow and fast subsystems. Moreover, the presented design insures the stability for all ε∈(0, ∞). We show by example that the new ε-bound can yield the maximum robust stability bound than that of the result [18].
The following notations will be used throughout the paper: the identity matrix with dimension n is denoted by In, λ(A) stands for the eigenvalues of matrix A, A denotes the norm of matrix A, i.e., A = M a x [ λ ( A T A ) ] 1 / 2 .

2. Problem Formulation and Design Procedure

Consider the linear time-invariant singularly perturbed bilinear system
[ x ˙ 1 ( t ) ε x ˙ 2 ( t ) ] = [ A 11 A 12 A 21 A 22 ] [ x 1 ( t ) x 2 ( t ) ] + [ N 11 N 12 N 21 N 22 ] [ x 1 ( t ) x 2 ( t ) ] u ( t ) + [ B 1 B 2 ] u ( t )
y ( t ) = C 1 x 1 ( t ) + C 2 x 2 ( t )
where x 1 R n 1 , x 2 R n 2 , n(= n1 + n2) is the order of the whole system, uR is the control, y R m is the output and the matrices A i j , N i j , B i and C i are constant matrices with appropriate dimensions for i, j = 1, 2, and ε is the singular perturbation parameter.

2.1. Nonlinear State Feedback

First, the main objective of this subsection is to design a nonlinear state feedback controller for the given system (1) to be stable with ε∈(0, ∞). Here, we will apply a control of the form
u ( t ) = k 1 D 1 x 1 1 + x 1 T D 1 T D 1 x 1 k 2 D 2 x 2 1 + x 2 T D 2 T D 2 x 2
to stabilize the system (1), where D i R 1 × n i can be arbitrarily designed, and K i are constants. To investigate the stability of system (1), we choose the Lyapunov function candidate as
V ( x 1 ,   x 2 ) = x 1 T P 1 x 1 + ε x 2 T P 2 x 2
where P i R n i × n i are unique real symmetric positive-definite matrices satisfying the Lyapunov equations,
A i i T P i + P i A i i = 2 I n i
for i = 1, 2, and A i i are assumed to be stable.
Remark 1.
The system’s input can be deterministic signals with unknown upper bounds, i.e., u ( t ) Ω and Ω { u ( t ) | | u ( t ) | | k 1 | + | k 2 | } .
Substitution of (2) into (1) yields an alternative form
[ x ˙ 1 x ˙ 2 ] = [ A ¯ 11 A ¯ 12 A ¯ 21 A ¯ 22 ] [ x 1 ( t ) x 2 ( t ) ]
where
A ¯ i i = A i i k i N i i sin θ i k j N i i sin θ j k i B i D i cos θ i
A ¯ i j = A i j k i N i j sin θ i k j N i j sin θ j k j B i D j cos θ j
sin θ i = x i T D i T 1 + x i T D i T D i x i , · cos θ i = 1 1 + x i T D i T D i x i   and   90 ° θ i 90 °   for   i = 1 ,   2 .
Before deriving the main result, one preliminary lemma is reviewed firstly.
Lemma 1 ([22]).
For any matrices X and Y with appropriate dimensions, we have
X T Y + Y T X γ X T X + Y T Y
for any constant γ >0.
Theorem 1.
Consider the singularly perturbed bilinear system (1) and assume that Aii is Hurwitz and P i A i j A i j T P i < 1 , if k 1 and k 2 are selected inside the region of the convex polygon as shown in Figure 1, the nonlinear state feedback controller (2) stabilizes the system (1) in the entire state space for all ε∈(0, ∞), where the four intersections of convex polygon with k 1 –axis and with k 2 –axis are described follows,
k ¯ 1 = ( 1 P 1 A 12 A 12 T P 1 ) / ( 2 N 11 T P 1 + 2 D 1 T B 1 T P 1 + P 1 N 12 N 12 T P 1 + 2 )
k ^ 1 = ( 1 P 2 A 21 A 21 T P 2 ) / ( 2 N 22 T P 2 + P 2 N 21 N 21 T P 2 + P 2 B 2 D 1 D 1 T B 2 T P 2 + 1 )
k ¯ 2 = ( 1 P 1 A 12 A 12 T P 1 ) / ( 2 N 11 T P 1 + P 1 N 12 N 12 T P 1 + P 1 B 1 D 2 D 2 T B 1 T P 1 + 1 )
k ^ 2 = ( 1 P 2 A 21 A 21 T P 2 ) / ( 2 N 22 T P 2 + 2 D 2 T B 2 T P 2 + P 2 N 21 N 21 T P 2 + 2 )
Proof. 
From (3), (5) and Lemma 1, the Lyapunov function’s derivative becomes
V ˙ = x ˙ 1 T P 1 x 1 + x 1 T P 1 x ˙ 1 + ε x ˙ 2 T P 2 x 2 + ε x 2 T P 2 x ˙ 2 = i = 1 2 [ ( A ¯ i i x i + A ¯ i j x j ) T P i x i + x i T P i ( A ¯ i i x i + A ¯ i j x j ) ]   ( for   j = 1 , 2 ,   and   j i ) = i = 1 2 [ 2 x i T ( A i i T P i k i N i i T P i sin θ i k j N i i T P i sin θ j k i D i T B i T P i cos θ i ) x i + 2 x j T ( A i j T P i k i N i j T P i sin θ i k j N i j T sin θ j k j D j T B i T P i cos θ j ) x i ] i = 1 2 x i T [ 1 + P i A i j A i j T P i + | k i | ( 2 N i i T P i + 2 C i T B i T P i + P i N i j N i j T P i + 2 ) + | k j | ( 2 N i i T P i + P i B i D j D j T B i T P i + P i N i j N i j T P i + 1 ) ] x i
It is easily verified that the following inequality holds such that V ˙ < 0
1 + P i A i j A i j T P i + | k i | ( 2 N i i T P i + 2 D i T B i T P i + P i N i j N i j T P i + 2 ) + | k j | ( 2 N i i T P i + P i B i D j D j T B i T P i + P i N i j N i j T P i + 1 ) < 0
From the above inequality (8), the eight pinnacles on the coordinate axes of k 1 k 2 plane form a convex polygon as shown in Figure 1. Hence, the controller (2), whose gains are selected inside the admissible region of Figure 1, can stabilize the system (1) for all ε∈(0, ∞). □
Remark 2.
Singular perturbation technique is usually a method to deal with two-time-scale systems, where the analysis of the whole system can be split, for “small” values of ε (or the limit ε 0 ), into the analysis of two separate systems presented in the slow and fast time scales. One tries to design a controller for the whole system working on reduced order control problems for the two limit systems. In this paper, we do not utilize the time-scale properties of the system and this is another reason why asymptotic stability is guaranteed also for “large” ε, though the two-time-scale structure of the system may be destroyed. In fact, no information is given to determine whether or not the two-time-scale separation is preserved in the closed-loop case as ε approaches the maximum upper bound ε*. Note that we seek the maximum upper bound ε* for the closed-loop system to retain stable, the value of ε* with which we end up might turn out to be not so small (sometimes even infinite) [23]. As such, the calculation of ε* serves only one purpose, namely it guarantees that the system (1) is well behaved in the generalized stability sense for all ε∈(0, ε*).
Remark 3.
The hypotheses employed by Theorem 1, i.e., asymptotic stability of matrices A 11   a n d   A 22 , and inequalities P 1 A 12 A 12 T P 1 < 1 ,   P 2 A 21 A 21 T P 2 < 1 , guarantee that the open-loop system is asymptotically stable for any value of ε∈(0,∞). One can notice that the polygon of admissible gains given in Theorem 1 includes the origin of the k 1 k 2 plane. One query may raise. What is the reason to apply a closed-loop control to stabilize a stable system (1)? There is at least a main reason to show the necessity of the feedback control: It is well known that an open-loop system is sensitive for some parameter variations or external input disturbances. In other words, if u(t) = 0, the bilinear system becomes an open-loop linear system, which may be unstable easily under some external disturbances. But by the given nonlinear feedback control (2), the robust stability of the system will increase [2].
Remark 4.
There have been some papers as well [24,25,26] discussing the stabilization problem of the bilinear system with the open-loop stable system.

2.2. Observer and Observer-Based State Feedback

Furthermore, the objective of the subsection is to design a full state observer and observer-based feedback for the given system (1) to be stable with ε∈(0, ∞). Now, it is required to construct an observer dynamical equation for the system (1) as below:
x ^ ˙ 1 = A 11 x ^ 1 + A 12 x ^ 2 + N 11 x ^ 1 u + N 12 x ^ 2 u + B 1 u + L 1 ( y C 1 x ^ 1 C 2 x ^ 2 )
ε x ^ ˙ 2 = A 21 x ^ 1 + A 22 x ^ 2 + N 21 x ^ 1 u + N 22 x ^ 2 u + B 2 u + L 2 ( y C 1 x ^ 1 C 2 x ^ 2 )
where x ^ 1 R n 1 ,   x ^ 2 R n 2 are the observer states, and L 1 R n 1 × q and L 2 R n 2 × q are suitable observing gain matrices that will be designed so as to ensure the estimated state x ^ 1 x 1   and   x ^ 2 x 2 as t . It should also be added that the pair ( A i , C i ) is detectable.
Here, we will apply an observer-based feedback of the form
u ( t ) = k 1 D 1 x ^ 1 1 + x ^ 1 T D 1 T C 1 x ^ 1 k 2 D 2 x ^ 2 1 + x ^ 2 T D 2 T D 2 x ^ 2 = k 1 sin θ ^ 1 k 2 sin θ ^ 2 = k 1 cos θ ^ 1 D 1 x ^ 1 k 2 cos θ ^ 2 D 2 x ^ 2
to stabilize the system (1), where D i R 1 × n i can be arbitrarily designed, k i is constants, sin θ ^ i = x ^ i T D i T 1 + x ^ i T D i T D i x ^ i ,   cos θ ^ i = 1 1 + x ^ i T D i T D i x ^ i and 90 ° θ ^ i 90 ° for i = 1, 2.
Hence, we obtain the state and error ( e 1 = x 1 x ^ 1   and   e 2 = x 2 x ^ 2 ) equations which can be written as follows,
[ x ˙ 1 ε x ˙ 2 ] = A ^ 11 [ x 1 x 2 ] + A ^ 12 [ e 1 e 2 ]
[ e ˙ 1 ε e ˙ 2 ] = A ^ 21 [ x 1 x 2 ] + A ^ 22 [ e 1 e 2 ]
where
A ^ 11 = [ A ¯ 11 A ¯ 12 A ¯ 21 A ¯ 11 ] ,   A ^ 12 = [ k 1 B 1 D 1 cos θ ^ 1 k 2 B 1 D 2 cos θ ^ 2 k 1 B 2 D 1 cos θ ^ 1 k 2 B 2 D 2 cos θ ^ 2 ] ,   A ^ 21 = 0   and   A ^ 22 = [ A ˜ 11 A ˜ 12 A ˜ 21 A ˜ 22 ] ,
in which A ˜ i i = A ¯ i i k i N i i sin θ ^ i k j N i i sin θ ^ j , A ˜ i j = A ¯ i j k i N i j sin θ ^ i k j N i j sin θ ^ j , and A ¯ i i = A i i L i C i ,   A ¯ i j = A i j L i C j for i , j = 1 , 2 , i j .
For the stability of the system (11), the eigenvalues of A ^ 11   and   A ^ 22 are designed within the unit circle, i.e., | λ ( A ^ 11 ) | < 1 and | λ ( A ^ 22 ) | < 1 . The objective is to find the suitable estimator matrices L 1 and L 2 such that the estimated error approaches zero. This means that the estimated state x ^ 1 and x ^ 2 can approach the real state x 1 and x 2 for all ε∈(0, ∞), respectively.
In order to prove the stability of the error Equation (11), the Lyapunov function candidate is given by
V ( e 1 , e 2 ) = e 1 T P o 1 e 1 + ε e 2 T P o 2 e 2
where P o i R n i × n i are unique real symmetric positive-definite matrices satisfying the Lyapunov equations,
A ¯ i i T P o i + P o i A ¯ i i = 2 I n i
for i = 1, 2, and A ¯ i i are assumed to be stable.
Theorem 2.
For the singularly perturbed bilinear systems (1) with observer-based controller (10), in which k 1 and k 2 can be calculated from Theorem 1. If there exist L 1 and L 2 such that
P o i A ¯ i j A ¯ i j T P o i < 1 ( 2 N i i T P o i + P o i N i j N i j T P o i + 1 ) ( | k 1 | + | k 2 | )
Then the system (11) is asymptotically stable for all ε∈(0, ∞).
Proof. 
The proof can be divided into two parts:
(i)
The stability of the system’s state Equation (11).
In fact, this result is the same as Theorem 1, i.e., if x 1 R n 1 are selected inside the region of the convex polygon as shown in Figure 1, then A ^ 11 can be guaranteed to be stable.
(ii)
The stability of the error Equation (11).
From (11), (12) and Lemma 1, the Lyapunov function’s derivative becomes
V ˙ ( e 1 , e 2 ) = e ˙ 1 T P o i e 1 + e 1 T P o i e ˙ 1 + ε e ˙ 2 T P o i e 2 + ε e 2 T P o i e ˙ 2 = i = 1 2 [ ( A ˜ i i e i + A ˜ i j e j ) T P o i e i + e i T P o i ( A ˜ i i e i + A ˜ i j e j ) ] ( for   j = 1 , 2 ,   and   j i ) = i = 1 2 [ 2 e i T ( A ¯ i i T P o i k i N i i T P o i sin θ ^ i k j N i i T P o i sin θ ^ j ) e i + 2 e j T ( A ¯ i j T k i N i j T P o i sin θ ^ i k j N i j T P o i sin θ ^ j ) e i ] i = 1 2 e i T [ 1 + P o i A ¯ i j A ¯ i j T P o i + ( | k i | + | k j | ) ( 2 N i i T P o i + P i N i j N i j T P i + 1 ) ] e i
It is easily verified that the following inequality holds such that V ˙ ( e 1 , e 2 ) < 0
P o i A ¯ i j A ¯ i j T P o i < 1 ( | k i | + | k j | ) ( 2 N i i T P o i + P i N i j N i j T P i + 1 )
Hence, as k 1 and k 2 are to be calculated from Theorem 1, and matrices L 1 and L 2 can be chosen to satisfy the above inequality, then the state and error system (11) are both asymptotically stable for all ε∈(0, ∞). □
Remark 5.
It should be noted that in Theorem 2, the values | k 1 | + | k 2 | will affect the choices of L 1 and L 2 . In other words, the larger  | k 1 | + | k 2 | is, the more negative  A ¯ i i should be, and the smaller  P o i A ¯ i j A ¯ i j T P o i should be.
Remark 6.
There has been some literature [27,28,29,30] discussing only the observer problem of the bilinear systems with bounded input, while the observer-based controller design can be developed as well in this paper.
The controller and observer can be designed in sequence in Theorem 2, but, we do not guarantee there exists L 1 and L 2 such that inequality (14) holds. In the following theorem, we change the sequence of design. In other words, first, matrices L 1 and L 2 can be arbitrarily designed such that A ¯ 11 and A ¯ 22 are both stable; second, the range of values k 1   and   k 2 can be determined by means of matrices L 1 and L 2 in Theorem 3.
Theorem 3.
For the singularly perturbed bilinear systems (1) with observer-based feedback (10), where L i is designed such that  A ¯ i i is stable. If  k 1   a n d   k 2 are selected inside the region of the convex polygon as shown in Figure 2, then the observer-based feedback (10) stabilizes the system (1) in the entire state space for all ε∈(0,∞), where the six intersections of convex polygon with k1–axis and with k2–axis are described as Equation (7), and
k ¯ = m i n ( k ¯ 1 , k ¯ 2 )
in which
k ¯ 1 = ( 1 P o 1 A ¯ 12 A ¯ 12 T P o 1 ) / ( 2 N 11 T P o 1 + P o 1 N 12 N 12 T P o 1 + 1 ) k ¯ 2 = ( 1 P o 2 A ¯ 21 A ¯ 21 T P o 2 ) / ( 2 N 22 T P o 2 + P o 2 N 21 N 21 T P o 2 + 1 )
Proof. 
The proof can be also divided into two parts:
(i)
The stability of the system’s state Equation (11).
In fact, the inequalities (16) can be obtained from Theorem 1. The results guarantee that A ^ 11 is stable.
(ii)
The stability of the error Equation (11).
After some mathematical manipulations from inequality (14), it is clear that A ^ 22 is stable if inequality (16) holds. □
Remark 7.
For the structures of the observer (9) and observer-based feedback (10), there are two different procedures to find the parameter matrix L i and constant gain  k i . In Theorem 2, we can first design the feedback gain  k i  according to Theorem 1, then the  L i  can be calculated by means of inequality (14), but we cannot guarantee there exists   L i  such that inequality (14) holds. On the contrary, in Theorem 3, the parameter matrix  L i  can be designed such that  A ¯ i i  is stable, then the feedback gain  k i  can be determined; hence, it is clear that the controller and observer are both compatible in Theorem 3.

3. Example

Consider a fourth-order singularly perturbed bilinear system:
[ x ˙ 1 ( 1 ) ( t ) x ˙ 1 ( 2 ) ( t ) ε x ˙ 2 ( 1 ) ( t ) ε x ˙ 2 ( 2 ) ( t ) ] = [ 1 0 0.2 0.1 0.5 1.5 0.3 0 0 0.2 3 1 0.3 0.1 0 2 ] [ x 1 ( 1 ) ( t ) x 1 ( 2 ) ( t ) x 2 ( 1 ) ( t ) x 2 ( 2 ) ( t ) ] + [ 1 0 0.5 1 2 1 0 2 0.5 1 0 0 0 2 0 0 ] [ x 1 ( 1 ) ( t ) x 1 ( 2 ) ( t ) x 2 ( 1 ) ( t ) x 2 ( 2 ) ( t ) ] u ( t ) + [ 1 0.5 1 2 ] u ( t )
y ( t ) = [ 1 0 ] [ x 1 ( 1 ) ( t ) x 1 ( 2 ) ( t ) ] + [ 1 0 ] [ x 2 ( 1 ) ( t ) x 2 ( 2 ) ( t ) ]
where x 1 = [ x 1 ( 1 ) x 1 ( 2 ) ] T and x 2 = [ x 2 ( 1 ) x 2 ( 2 ) ] T . First, according to the Lyapunov Equation (4), we can obtain P 1 = [ 1.0667 0.1333 0.1333 0.6667 ] and P 2 = [ 0.3333 0.0667 0.0667 0.5333 ] . Let D i = [ 1 1 ] for i = 1, 2, and the four intersections, k ¯ 1 = 0.0927 ,   k ¯ 2 = 0.1129 ,   k ^ 1 = 0.2361   and   k ^ 2 = 0.1642 can be calculated from (7). Therefore, we conclude that system (17) can be stabilized by the following feedback control law for all ε∈(0, ∞):
u ( t ) = k 1 ( x 1 ( 1 ) + x 1 ( 2 ) ) 1 + ( x 1 ( 1 ) + x 1 ( 2 ) ) 2 k 2 ( x 2 ( 1 ) + x 2 ( 2 ) ) 1 + ( x 2 ( 1 ) + x 2 ( 2 ) ) 2
where k1 and k2 are selected inside the admissible region of the convex polygon as shown in Figure 3. Whereas the result of [14] promises the stability for ε∈(0, 1) only.
Furthermore, two different procedures for the design of observer and observer-based state feedback for the system (17) can be presented as follows.
(I)
First, by Theorem 1, let k1 = 0.05, k2 = 0.05, and choose L i = [ 1 0 ] T such that A ^ i i is stable, from the Lyapunov Equation (13), we can obtain P o 1 = [ 0.5238 0.0952 0.0952 0.6667 ] and P o 2 = [ 0.25 0.0417 0.0417 0.5208 ] . The inequality (14) can be held according to the above parameters, hence, an observer equation can be represented as below:
[ x ^ ˙ 1 ( 1 ) ( t ) x ^ ˙ 1 ( 2 ) ( t ) ε x ^ ˙ 2 ( 1 ) ( t ) ε x ^ ˙ 2 ( 2 ) ( t ) ] = [ 1 0 0.2 0.1 0.5 1.5 0.3 0 0 0.2 3 1 0.3 0.1 0 2 ] [ x ^ 1 ( 1 ) ( t ) x ^ 1 ( 2 ) ( t ) x ^ 2 ( 1 ) ( t ) x ^ 2 ( 2 ) ( t ) ] + [ 1 0 0.5 1 2 1 0 2 0.5 1 0 0 0 2 0 0 ] [ x ^ 1 ( 1 ) ( t ) x ^ 1 ( 2 ) ( t ) x ^ 2 ( 1 ) ( t ) x ^ 2 ( 2 ) ( t ) ] u ( t ) + [ 1 0.5 1 2 ] u ( t ) + [ 1 0 1 0 ] ( y ( t ) [ 1 0 1 0 ] [ x ^ 1 ( 1 ) ( t ) x ^ 1 ( 2 ) ( t ) x ^ 2 ( 1 ) ( t ) x ^ 2 ( 2 ) ( t ) ] )
Therefore, we conclude that system (17) can be stabilized by the following feedback control law for all ε∈(0, ∞);
u ( t ) = 0.05 ( x ^ 1 ( 1 ) + x ^ 1 ( 2 ) ) 1 + ( x ^ 1 ( 1 ) + x ^ 1 ( 2 ) ) 2 0.05 ( x ^ 2 ( 1 ) + x ^ 2 ( 2 ) ) 1 + ( x ^ 2 ( 1 ) + x ^ 2 ( 2 ) ) 2
(II)
First, we can choose L i = [ 1 1 ] T such that A ^ i i is stable, and matrices P o 1 = [ 0.7143 0.2857 0.2857 0.6667 ]   and   P o 2 = [ 0.2857 0.1429 0.1429 0.5714 ] can be also obtained. Hence, an observer of the system (17) for all ε∈(0, ∞) can be represented as below:
[ x ^ ˙ 1 ( 1 ) ( t ) x ^ ˙ 1 ( 2 ) ( t ) ε x ^ ˙ 2 ( 1 ) ( t ) ε x ^ ˙ 2 ( 2 ) ( t ) ] = [ 1 0 0.2 0.1 0.5 1.5 0.3 0 0 0.2 3 1 0.3 0.1 0 2 ] [ x ^ 1 ( 1 ) ( t ) x ^ 1 ( 2 ) ( t ) x ^ 2 ( 1 ) ( t ) x ^ 2 ( 2 ) ( t ) ] + [ 1 0 0.5 1 2 1 0 2 0.5 1 0 0 0 2 0 0 ] [ x ^ 1 ( 1 ) ( t ) x ^ 1 ( 2 ) ( t ) x ^ 2 ( 1 ) ( t ) x ^ 2 ( 2 ) ( t ) ] u ( t ) + [ 1 0.5 1 2 ] u ( t ) + [ 1 1 1 1 ] ( y ( t ) [ 1 0 1 0 ] [ x ^ 1 ( 1 ) ( t ) x ^ 1 ( 2 ) ( t ) x ^ 2 ( 1 ) ( t ) x ^ 2 ( 2 ) ( t ) ] )
Second, according to Theorem 3, the six intersections, k ¯ 1 = 0.0927 ,   k ¯ 2 = 0.1129 ,   k ^ 1 = 0.2361 ,   k ^ 2 = 0.1642 and k ¯ = 0.0927 can be calculated from (16). Therefore, we conclude that system (17) can be stabilized by the following feedback control law for all ε∈(0, ∞):
u ( t ) = k 1 ( x ^ 1 ( 1 ) + x ^ 1 ( 2 ) ) 1 + ( x ^ 1 ( 1 ) + x ^ 1 ( 2 ) ) 2 k 2 ( x ^ 2 ( 1 ) + x ^ 2 ( 2 ) ) 1 + ( x ^ 2 ( 1 ) + x ^ 2 ( 2 ) ) 2
where k 1   and   k 2 are selected inside the admissible region of the convex polygon, as shown in Figure 4.

4. Conclusions

According to the Lyapunov stability theorem, the designs of observer and observer-based stabilizers are presented for the singularly perturbed bilinear systems in this paper. The main advantage of this paper is that the whole n 1 + n 2 dimensional process can be separated into two parts ( n 1 and n 2 dimensional), which do not require the slow-fast decomposition. i.e., the main advantage of the proposed approach is that the whole dimensional design can be directly separated into two lower dimensions without the decomposition and composition of the slow and fast subsystems. Besides, this paper has three major contributions as well for the observer and controller design of the singularly perturbed bilinear systems. The first one is that the matrix N 22 = 0 is moved [18]. The second one is that the observer and closed-loop system are both stable for all ε∈(0, ∞), by those methods, the singular perturbation parameter is treated as a structured uncertainty. The last one is that we propose two different procedures for determining observer matrix L i and feedback gain k i in the structures of the observer (9) and the observer-based stabilizer (10), respectively.

Author Contributions

Methodology, J.-S.C.; validation, C.-T.C. and D.-H.C.; formal analysis, J.-S.C.; investigation, J.-S.C.; writing—original draft preparation, C.-T.C.; writing—review and editing, D.-H.C.; visualization, D.-H.C.; supervision, J.-S.C.; project administration, J.-S.C.; funding acquisition, J.-S.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research is supported by the Ministry of Science and Technology, Taiwan, Grant MOST 109-2221-E-218-020-MY2 and MOST 109-2622-E-218-006-CC2.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mohler, R.R. Nonlinear Systems, Vol. II, Applications to Bilinear Control; Prentice Hall: Englewood Cliffs, NJ, USA, 1991. [Google Scholar]
  2. Wang, W.-J.; Chiou, J.-S. A non-linear control design for the discrete-time multi-input bilinear systems. Mechatronics 1991, 1, 87–94. [Google Scholar] [CrossRef]
  3. Benesty, J.; Paleologu, C.; Dogariu, L.-M.; Ciochină, S. Identification of Linear and Bilinear Systems: A Unified Study. Electronics 2021, 10, 1790. [Google Scholar] [CrossRef]
  4. Nie, Z.; Gao, F.; Yan, C.-B. A Multi-Timescale Bilinear Model for Optimization and Control of HVAC Systems with Consistency. Energies 2021, 14, 400. [Google Scholar] [CrossRef]
  5. Kokotovic, P.V.; Khalil, H.K.; O’Reilly, J. Singular Perturbation Methods in Control: Analysis and Design; Academic Press: London, UK, 1986. [Google Scholar]
  6. Naidu, D.S. Singular Perturbation Methodology in Control Systems; IEE and Perter Peregrinus: Stevenage Herts, UK, 1988. [Google Scholar]
  7. Chao, C.-T.; Chen, D.-H.; Chiou, J.-S. Stabilization and the Design of Switching Laws of a Class of Switched Singularly Perturbed Systems via the Composite Control. Mathematics 2021, 9, 1664. [Google Scholar] [CrossRef]
  8. Bobodzhanov, A.; Kalimbetov, B.; Safonov, V. Generalization of the Regularization Method to Singularly Perturbed Integro-Differential Systems of Equations with Rapidly Oscillating Inhomogeneity. Axioms 2021, 10, 40. [Google Scholar] [CrossRef]
  9. Liu, W.; Wang, Y.; Tian, Y. Feedback Stabilization for Discrete-Time Singularly Perturbed Systems via Limited Information with Data Packet Dropout. IEEE Access 2021, 9, 40585–40594. [Google Scholar] [CrossRef]
  10. Cerpa, E.; Prieur, C. Singular Perturbation Analysis of a Coupled System Involving the Wave Equation. IEEE Trans. Autom. Control 2020, 65, 4846–4853. [Google Scholar] [CrossRef]
  11. Song, J.; Niu, Y. Dynamic Event-Triggered Sliding Mode Control: Dealing with Slow Sampling Singularly Perturbed Systems. IEEE Trans. Circuits Syst. II Express Briefs 2020, 67, 1079–1083. [Google Scholar] [CrossRef]
  12. Dragan, V. On the Linear Quadratic Optimal Control for Systems Described by Singularly Perturbed Itô Differential Equations with Two Fast Time Scales. Axioms 2019, 8, 30. [Google Scholar] [CrossRef] [Green Version]
  13. Abdullazade, N.N.; Chechkin, G.A. Perturbation of the Steklov problem on a small part of the boundary. J. Math. Sci. 2014, 196, 441–450. [Google Scholar] [CrossRef]
  14. Borisov, D.I.; Cardone, G.; Checkhin, G.A.; Koroleva, Y.O. On elliptic operators with Steklov condition perturbed by Dirichlet condition on a small part of boundary. Calc. Var. Partial. Differ. Equ. 2021, 60, 1–44. [Google Scholar] [CrossRef]
  15. Chechkina, A.G. On singular perturbation of a Steklov-type problem with asymptotically degenerate spectrum. Dokl. Math. 2011, 84, 695–698. [Google Scholar] [CrossRef]
  16. Chechkina, A.G. Homogenization of spectral problems with singular perturbation of the Steklov condition. Izvestiya: Mathematics 2017, 81, 199–236. [Google Scholar] [CrossRef]
  17. Tzafestat, S.G.; Anagonostou, K.E. Stabilization of singularly perturbed strictly bilinear systems. IEEE Trans. Automat. Contr. 1984, 29, 943–946. [Google Scholar] [CrossRef]
  18. Asamoah, F.; Jamshidi, M. Stabilization of a class of singularly perturbed bilinear systems. Int. J. Contr. 1987, 46, 1589–1594. [Google Scholar] [CrossRef]
  19. Li, T.-H.S.; Sun, Y.-Y. Comments on: Stabilization of a class of singularly perturbed bilinear systems. Int. J. Contr. 1988, 48, 1357–1358. [Google Scholar] [CrossRef]
  20. Chiou, J.-S.; Kung, F.-C.; Li, T.-H.S. Robust stabilization of a class of singularly perturbed discrete bilinear systems. IEEE Trans. Automat. Contr. 2000, 45, 1187–1191. [Google Scholar] [CrossRef]
  21. Pan, S.-J.; Tsai, J.S.-H.; Pan, S.-T. Application of Genetic Algorithm on Observer-Based D-Stability Control for Discrete Multiple Time-Delay Singularly Perturbation Systems. Int. J. Innov. Comput. Inf. Control 2011, 7, 3345–3358. [Google Scholar]
  22. Zhou, K.; Khargonekar, P.P. Robust stabilization of linear systems with norm-bounded time-varying uncertainty. Sys. Control Lett. 1988, 10, 17–20. [Google Scholar] [CrossRef]
  23. Saydy, L. New stability/performance results for singularly perturbed system. Automatica 1996, 32, 807–818. [Google Scholar] [CrossRef]
  24. Longchamp, R. Controller design for bilinear systems. IEEE Trans. Automat. Control 1980, 25, 547–548. [Google Scholar] [CrossRef]
  25. Genesio, R.; Tesi, A. The output stabilization of SISO bilinear systems. IEEE Trans. Automat. Control 1988, 33, 950–952. [Google Scholar] [CrossRef]
  26. Benallou, A.; Mellichamp, D.A.; Seborg, D.E. Optimal stabilizing controllers for bilinear systems. Int. J. Control 1988, 48, 1487–1501. [Google Scholar] [CrossRef]
  27. Hara, S.; Fumta, K. Minimal order state observers for bilinear systems. Int. J. Control 1976, 24, 705–718. [Google Scholar] [CrossRef]
  28. Funahashi, Y. Stable state estimator for bilinear system. Int. J. Control 1979, 29, 181–188. [Google Scholar] [CrossRef]
  29. Derese, I.; Noldus, E.J. Existence of bilinear state observers for bilinear system. IEEE Trans. Automat. Control 1981, 26, 590–592. [Google Scholar] [CrossRef]
  30. Derse, I.; Stevens, P.; Noldus, E.J. Observers for bilinear systems with bounded input. Int. J. Syst. Sci. 1979, 10, 649–668. [Google Scholar] [CrossRef]
Figure 1. The admissible region k1 and k2 for the controller (2).
Figure 1. The admissible region k1 and k2 for the controller (2).
Mathematics 09 02380 g001
Figure 2. The admissible region k1 and k2 for the observer-based controller (10).
Figure 2. The admissible region k1 and k2 for the observer-based controller (10).
Mathematics 09 02380 g002
Figure 3. The admissible region k1 and k2 for the controller of example 1.
Figure 3. The admissible region k1 and k2 for the controller of example 1.
Mathematics 09 02380 g003
Figure 4. The admissible region k1 and k2 for the observer-based controller of example 1.
Figure 4. The admissible region k1 and k2 for the observer-based controller of example 1.
Mathematics 09 02380 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, D.-H.; Chao, C.-T.; Chiou, J.-S. Robust Stabilization and Observer-Based Stabilization for a Class of Singularly Perturbed Bilinear Systems. Mathematics 2021, 9, 2380. https://doi.org/10.3390/math9192380

AMA Style

Chen D-H, Chao C-T, Chiou J-S. Robust Stabilization and Observer-Based Stabilization for a Class of Singularly Perturbed Bilinear Systems. Mathematics. 2021; 9(19):2380. https://doi.org/10.3390/math9192380

Chicago/Turabian Style

Chen, Ding-Horng, Chun-Tang Chao, and Juing-Shian Chiou. 2021. "Robust Stabilization and Observer-Based Stabilization for a Class of Singularly Perturbed Bilinear Systems" Mathematics 9, no. 19: 2380. https://doi.org/10.3390/math9192380

APA Style

Chen, D. -H., Chao, C. -T., & Chiou, J. -S. (2021). Robust Stabilization and Observer-Based Stabilization for a Class of Singularly Perturbed Bilinear Systems. Mathematics, 9(19), 2380. https://doi.org/10.3390/math9192380

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