Next Article in Journal
Minimization Problems for Functionals Depending on Generalized Proportional Fractional Derivatives
Next Article in Special Issue
Finite Difference–Collocation Method for the Generalized Fractional Diffusion Equation
Previous Article in Journal
Multivalent Functions and Differential Operator Extended by the Quantum Calculus
Previous Article in Special Issue
Finite Time Stability of Fractional Order Systems of Neutral Type
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fractional-Order Interval Observer for Multiagent Nonlinear Systems

The School of Mechanical and Electrical Engineering, Soochow University, Suzhou 215131, China
*
Author to whom correspondence should be addressed.
Submission received: 31 May 2022 / Revised: 20 June 2022 / Accepted: 23 June 2022 / Published: 25 June 2022
(This article belongs to the Special Issue Fractional-Order System: Control Theory and Applications)

Abstract

:
A framework of distributed interval observers is introduced for fractional-order multiagent systems in the presence of nonlinearity. First, a frame was designed to construct the upper and lower bounds of the system state. By using monotone system theory, the positivity of the error dynamics could be ensured, which implies that the bounds could trap the original state. Second, a sufficient condition was applied to guarantee the boundedness of distributed interval observers. Then, an extension of Lyapunov function in the fractional calculus field was the basis of the sufficient condition. An algorithm associated with the procedure of the observer design is also provided. Lastly, a numerical simulation is used to demonstrate the effectiveness of the distributed interval observer.

1. Introduction

Fractional calculus is the definition of the differential and integral of arbitrary order systems. During the last 20 years, fractional calculus grew into a hot topic and attracted increasing interest [1,2]. It plays an important role when people are faced with natural dynamics problems. It works well when it is applied in describing the memory and hereditary properties of manifold materials. In addition, the applications of fractional calculus are rapidly expanding, such as non-Fickian dynamics [3], fractional boundary value problems [4], and variable-order thermostat models [5].
Due to the inaccessibility of the system state in many physical backgrounds and the availability of output information, research on observation problems is very significant. With the indepth study of observation problems, interval observers are widely accepted as an efficient tool in reconstructing the system state with nonlinearity or bounded uncertainty [6,7]. The concept of interval observers was first proposed by Gouzé, and it was successfully applied to uncertain biological systems [8]. After that, research on interval observers can be divided into two parts. The first is based on monotone system theory, and the second was developed from set-membership estimation. The monotone-system-theory-based method requires researchers to design suitable observer gains in order to ensure that the error dynamics is positive and bounded. In this case, the bounds from the interval observer converge to the original system. Inspired by this idea, there are many new techniques for various systems emerging. For linear time-invariant systems, the time-varying coordinate transformation technique was used in interval observer design [9]. Efimov et al. also extended the interval estimation technique to nonlinear time-varying systems [10]. An interval observer for switched systems was also introduced in [11]. On the other hand, a set-membership estimation-based method combined robust observer design with reachability analysis [12]. Interval observers in this framework obviously improved estimation accuracy. Additionally, scholars presented a new integrated version of interval observers that achieved an expected tradeoff between robust estimation conservatism and computational complexity [13]. In general, research on interval observers for integer-order systems is fruitful, which is instructive for us to design interval observers of fractional-order systems.
A tremendous expansion for multiagent systems (MASs) has been witnessed, which mainly regards output regulation [14], observer design [15] and consensus control [16]. It started from simple two-order systems and attracted more attention on various complex MASs, such as switched [17], nonlinear [18], and time-delaying [19] systems. The observers designed for MASs are called distributed observers if the communication topology is taken into account. A distributed observer was designed for leader-following MASs, and an actual vehicle model was used to verify the function of the distributed observer [20]. In [21], the existence condition of the distributed observer was established to recover the state of nonlinear MASs. There are some works for fractional-order MASs in [22] and [23]. The consensus control problem for fractional-order MASs under a fixed topology was studied by stability theory and the Lyapunov method [22]. For unknown nonlinear dynamics, an adaptive control protocol was introduced that was used in a trajectory tracking problem [23]. There were also some interesting results about the observation problem in [24,25,26,27]. Traditional observers and distributed controllers were both used in fractional-order MASs [24]. Compared with [24], the time-varying formation control was investigated, and the consensus protocol was designed to match a switching topology [25]. For the fractional-order heterogeneous nonlinear MASs, the consensus control problem was converted into the stabilisation problem using the distributed control technique [26]. In [27], the object was the fractional-order system with unknown orders, and a robust observer-based controller was formulated to solve the consensus problem.
For these works [24,25,26,27], the introduced observers were all Luenberger-like. When the original system suffers from uncertain disturbance or nonlinearity, Luenberger-like observers are unable to recover the system state. It is challenging to recover the system state when the general fractional-order system has a complex communication with its neighborhood. Therefore, a distributed interval observer is a better choice for fractional-order MASs. In many physical problems such as trajectory tracking or output regulation, distributed interval observers provide a more accurate state information, which is a great step forward for the consensus control protocol design. However, there is not any paper focusing on distributed interval observers for fractional-order MASs. In the field of fractional calculus, there are some interesting works about fractional-order systems, including state reconstruction and algorithm programming [28,29,30]. Motivated by the interval observer mentioned above, we developed a framework of distributed interval observers for the general fractional-order MASs. The specific contributions are outlined as follows:
1.
By investigating the fractional differential problem, a novel fractional-order Lyapunov function is proposed for the boundedness problem in the observer deign, which is an approach to prove that a matrix is Hurwitz.
2.
Different from linear MASs, nonlinear MASs are considered, and the solution for Lipschitz functions is combined with the distributed interval observer design.
3.
For fractional-order MASs, the communication topology is applied to the observer design. Each observer of the corresponding agent could accept the information from its adjacent observers. A novel distributed interval observer was first designed for the fractional-order MASs.
The rest of the paper is structured as follows. Section 2 mainly presents some preparation work, including the fractional differential, graph theory and fractional-order systems. Section 3 proposes the distributed interval observer design method for fractional-order MASs. An illustrative example is given in Section 4 to verify the validity of the designed observer. Lastly, conclusions and future work are drawn in Section 5.
Notations: For two vectors x R n and y R n , x < ( ) y are understood with x i < ( ) y i , i { 1 , , n } . For  A 0 and A 0 , symbol ( ) means that matrix A is positive(negative) definite. For a matrix B, there exist three properties: B + = m a x ( B , 0 ) , B = B + B and | B | = B + + B . Matrix B T denotes the transpose of B, and H e ( B ) is defined as H e ( B ) = B + B T . λ m i n ( A ) represents the minimal nonzero eigenvalue of A. ⊗ is the Kronecker product used in MASs. 1 N means a N-order column vector, and all elements are 1.

2. Preliminaries

2.1. Fractional Calculus

There are several definitions used to regard the fractional differentiation operator, such as the Gr u ¨ nwald–Letnikov, Riemann–Liouville and Caputo definitions introduced in [31]. Among them, the Caputo definition is the most frequently used and has many applications in the engineering field. In the Caputo definition, the initial value of the differentiation is considered.
The Caputo fractional derivative for a function f ( t ) is
t 0 C D t α f ( t ) = 1 Γ ( n α ) t 0 t f ( n ) ( τ ) ( t τ ) α n + 1 d τ ,
where t 0 C D t α is the fractional integration differentiation operator, t 0 represents the initial time, and α [ n 1 , n ] is the order of the system. Γ ( · ) is a gamma function that is introduced in Definition 1. f ( t ) is a differentiable function with the nth derivative. In this paper, we mainly focus on the fractional-order system with n = 1 and t 0 = 0 , then we have 0 < α < 1 . When n = 1 , f ( t ) is only required to have the first derivative. For the sake of simplicity, the differentiation operator t 0 C D t α is replaced by D t α .
The Gr u ¨ nwald-Letnikov fractional derivative for a function f ( t ) is
t 0 G L D t α f ( t ) = lim N [ t t 0 N ] α j = 0 N 1 ( 1 ) j f ( t j [ t t 0 N ] ) .
The Riemann-Liouville fractional derivative for a function f ( t ) is
t 0 R L D t α f ( t ) = 1 Γ ( α ) t 0 t ( t τ ) α 1 f ( τ ) d τ , α < 0 f ( t ) , α = 0 D t n [ t 0 R L D t α n f ( t ) ] , α > 0
with n 1 α < n .
For a wide class of functions, the Gr u ¨ nwald–Letnikov and Riemann–Liouville definitions are equivalent [32]. However, it is difficult for the Gr u ¨ nwald-Letnikov definition to have a Laplace transform. The Laplace transform of the Riemann–Liouville definition is
L [ 0 R L D t α f ( t ) ] = s α F ( s ) , q 0 s α F ( s ) k = 0 n 1 s k 0 R L D t α k 1 f ( 0 ) , n 1 q < n
where F ( s ) is the Laplace transform of f ( t ) .
The Laplace transform of the Caputo definition is
L [ 0 C D t α ] = s α F ( s ) k = 0 n 1 s α 1 k f ( k ) ( 0 ) .
Comparing (4) and (5), it is obvious that the Riemann–Liouville fractional derivative is unsuitable for the Laplace transform technique because it requires the knowledge of noninteger-order derivatives of the function at t = 0 , while the Caputo fractional derivative only requires the knowledge of integer-order derivatives of the function. This is why Caputo fractional derivative was chosen.
Definition 1
([33]). Gamma function is an important element for the Caputo fractional differentiation operator, which is defined by
Γ ( a ) = 0 e t t a 1 d t ,
where value Γ ( a ) is in the convergence of the right-hand side of the complex plane.

2.2. Graph Theory

A communication topology G ( N , ξ , A ) contains three elements: N { 1 , , N } is the node set of MASs, ξ N × N represents an edge set for nodes in MASs, and A R N × N is an adjacent matrix of the graph. For A = [ a i j ] R N × N , a j i is the weight from node i to node j. In this paper, the graph has no self loops, i.e.,  a i i = 0 , i N . a j i 0 if and only if ( i , j ) ξ . For node i, its neighborhood is denoted as N i = { j | j N , a j i 0 , j i } . For G , it is defined as a strongly connected graph if and only if any node in it is mutually reachable. For adjacent matrix A , the corresponding Laplacian matrix is defined as L = [ l i j ] N × N , where l i i = Σ j = 1 N a i j and l i j = a i j , i j . If a topology is call as a balanced topology, the edges of it must be balanced, which means that edge ( i , j ) belongs to set ξ if edge ( j , i ) belongs to set ξ .
Lemma 1
([34]). for a Laplacian matrix L , zero is one of the eigenvalues and it has a fixed right eigenvector 1 N . The other nonzero eigenvalues are all positive. If there exists a directed spanning tree in G , zero is a simple eigenvalue of L .
Lemma 2
([35]). considering G as a strongly connected graph, suppose that r = [ r 1 , , r N ] is the left eigenvector connected with the eigenvalue zero. Then, we have R L + L T R 0 , where R = d i a g { r 1 , , r N } . For a balanced graph, we have r 1 = = r N .
Lemma 3
([35]). If G is a strongly connected graph; its generalized algebraic connectivity is define as a ( L ) = m i n r T x = 0 , x 0 x T ( R L + L T R ) x 2 x T R x . Due to the topology in this paper being defined as balanced, matrix R can be converted into R = r 1 I N . The generalized algebraic connectivity is equal to a ( L ) = λ m i n ( H e ( L ) 2 ) .

2.3. Fractional-Order Systems

Consider the following fractional-order system for the i-th agent.
D t α x i ( t ) = A x i ( t ) + f ( x i ( t ) ) , y i ( t ) = C x i ( t ) ,
where x i R n is the system state, y i R m is the output, and f ( x i ) R n is the Lipschitz function. A R n × n and C R m × n are matrices with suitable dimensions. Communication topology G for (7) is a strongly connected balanced graph.
Property 1.
Given a matrix M R n × n , M is a Metzler matrix if all its elements outside the main diagonal are non-negative. For example, the following matrix is Metzler:
M = 1 1 2 0.5 .
Property 2.
Given a matrix N R n × n , M is a Hurwitz matrix if its all real parts of the eigenvalues are negative. For example, the following matrix is Hurwitz:
N = 3 2 2 3 .
For nonlinearity f ( x i ) , we began the analysis with the following properties.
Property 3
([36]). the differentiable global Lipschitz function f ( x i ) can be divided into two increasing functions a ( x i ) and b ( x i ) ; their relationship is
f ( x i ) = a ( x i ) b ( x i ) .
Property 4
([36]). for a Lipschitz function f ( x i ) , there exists a differentiable global Lipschitz function f ˜ ( x i k 1 , x i k 2 ) , such that
  • f ˜ ( x i k 1 , x i k 2 ) = a ( x i k 1 ) b ( x i k 2 ) ,
  • f ˜ ( x i , x i ) = f ( x i ) ,
  • f ˜ a ( x i ) 0 and f ˜ b ( x i ) 0 .
Remark 1.
In Property 3, a Lipschitz function is transformed into two increasing functions, which is just for us to introduce the function f ˜ ( · , · ) . Because a ( x i ) and b ( x i ) are both increasing functions, one can deduce that f ˜ a ( x i ) 0 and f ˜ b ( x i ) 0 easily. For  x ̲ x x ¯ , we have f ˜ ( x ̲ , x ¯ ) f ˜ ( x , x ) f ˜ ( x ¯ , x ̲ ) , which are the upper and lower bounds of the Lipschitz function f ( x ) in the structure of distributed interval observers.
Assuming that the bounds of the system state satisfy x ̲ i x i x ¯ i , on the basis of Properties 3 and 4, one can deduce that
f ˜ ( x ̲ i , x ¯ i ) f ˜ ( x i , x i ) = f ( x i ) f ˜ ( x ¯ i , x ̲ i ) .
By using the generalized Taylor formula, f ˜ ( x ¯ i , x ̲ i ) f ˜ ( x i , x i ) is written as
f ˜ ( x ¯ i , x ̲ i ) f ˜ ( x i , x i ) = 0 1 f ˜ δ 1 ( τ δ 1 + ( 1 τ ) δ 2 ) d τ ( δ 1 δ 2 ) ,
where δ 1 = x ¯ x ̲ and δ 2 = x x . It follows from Property 4 that
f ˜ ( x ¯ i , x ̲ i ) f ˜ ( x i , x i ) = [ 0 1 a x ( τ x ¯ + ( 1 τ ) x ) d τ 0 1 b x ( τ x + ( 1 τ ) x ̲ ) d τ ] ( δ 1 δ 2 ) = [ F 1 ( x ¯ , x ) F 2 ( x , x ̲ ) ] ( δ 1 δ 2 ) .
Similarly, we have
f ˜ ( x i , x i ) f ˜ ( x ¯ i , x ̲ i ) = [ F 3 ( x , x ̲ ) F 4 ( x ¯ , x ) ] ( δ 1 δ 2 ) ,
where matrices F i , i { 1 , , 4 } are non-negative and can be derived from (11). On the basis of the above discussion, the following property is presented.
Property 5
([36]). for f ( x i ) and f ˜ ( · , · ) defined in Property 4, if Jacobian matrix f ˜ δ 1 is bounded, there exist matrices F i , i { 1 , , 4 } such that
f ˜ ( x ¯ i , x ̲ i ) f ˜ ( x i , x i ) F 1 e ¯ i + F 2 e ̲ i , f ˜ ( x i , x i ) f ˜ ( x ̲ i , x ¯ i ) F 3 e ¯ i + F 4 e ̲ i ,
where e ¯ i = x ¯ i x i and e ̲ i = x i x ̲ i .
Example 1.
For nonlinear function f ( x ) = s i n x , corresponding functions f ˜ ( x ¯ , x ̲ ) and f ˜ ( x ̲ , x ¯ ) are defined as
f ˜ ( x ¯ , x ̲ ) = s i n ( x ¯ ) + x ¯ x ̲ , f ˜ ( x ̲ , x ¯ ) = s i n ( x ̲ ) + x ̲ x ¯ ,
where f ( x ) , f ˜ ( x ̲ , x ¯ ) and f ˜ ( x ¯ , x ̲ ) satisfy Properties 3 and 4.
The functions mentioned in Property 3 are defined as a ( x ) = s i n ( x ) + x and b ( x ) = x . a ( x ) and b ( x ) are both obviously increasing functions. Then, f ˜ ( x ¯ , x ̲ ) f ˜ ( x , x ) can be transformed into
f ˜ ( x ¯ , x ̲ ) f ˜ ( x , x ) = s i n ( x ¯ ) + x ¯ x ̲ s i n ( x ) = ( s i n ( x ¯ ) s i n ( x ) ) + ( x ¯ x ) + ( x x ̲ ) 2 ( x ¯ x ) + ( x x ̲ ) .
From Property 5, we have f ˜ ( x ¯ , x ̲ ) f ˜ ( x , x ) F 1 e ¯ + F 2 e ̲ . Combining (14) with Property 5, we chose F 1 = 2 I and F 2 = I . Similarly, we have F 3 = I and F 4 = 2 I .
From Example 1, Property 5 is feasible, and the result from Formulas (9)–(12) stands.
For System (7), nonlinear function f ( x i ( t ) ) was assumed to satisfy Properties 3–5; then, the interval observer for (7) was designed with
D t α x ¯ i ( t ) = A x ¯ i ( t ) + L ( y i ( t ) C x ¯ i ( t ) ) + γ M i = 1 N a i j ( x ¯ j ( t ) x ¯ i ( t ) ) + f ˜ ( x ¯ i , x ̲ i ) , D t α x ̲ i ( t ) = A x ̲ i ( t ) + L ( y i ( t ) C x ̲ i ( t ) ) + γ M i = 1 N a i j ( x ̲ j ( t ) x ̲ i ( t ) ) + f ˜ ( x ̲ i , x ¯ i ) ,
where M and L are interval observer gains.
The error dynamics of the interval observer is
D t α e ¯ i ( t ) = D t α x ¯ i ( t ) D t α x i ( t ) = ( A L C γ M j = 1 N L i j ) e ¯ i ( t ) + f ˜ ( x ¯ i , x ̲ i ) f ˜ ( x i , x i ) , D t α e ̲ i ( t ) = D t α x i ( t ) D t α x ¯ i ( t ) = ( A L C γ M j = 1 N L i j ) e ̲ i ( t ) + f ˜ ( x i , x i ) f ˜ ( x ̲ i , x ¯ i ) .
Then, e ( t ) , f ˜ ( x ¯ , x ̲ ) , f ˜ ( x ̲ , x ¯ ) and f ˜ ( x , x ) are defined as
e ( t ) = e 1 e N , f ˜ ( x ¯ , x ̲ ) = f ˜ ( x ¯ 1 , x ̲ 1 ) f ˜ ( x ¯ N , x ̲ N ) , f ˜ ( x ̲ , x ¯ ) = f ˜ ( x ̲ 1 , x ¯ 1 ) f ˜ ( x ̲ N , x ¯ N ) , f ˜ ( x , x ) = f ˜ ( x 1 , x 1 ) f ˜ ( x N , x N ) .
System (16) can be written in the following form:
D t α e ¯ ( t ) = ( I N ( A L C ) γ ( L M ) e ¯ ( t ) + f ˜ ( x ¯ , x ̲ ) f ˜ ( x , x ) , D t α e ̲ ( t ) = ( I N ( A L C ) γ ( L M ) ) e ̲ ( t ) + f ˜ ( x , x ) f ˜ ( x ̲ , x ¯ ) .

3. Main Results

Proving the boundedness of the distributed interval observer is equal to proving the convergence of the error dynamics. The Lyapunov function is a great choice for it. For error dynamics e ˙ ( t ) = A e e ( t ) , the following Lyapunov function is constructed:
V ( t ) = e T ( t ) P e ( t ) ,
where P is a positive symmetric matrix with suitable dimensions.
In the integer-order system, V ˙ ( t ) < 0 is the sufficient condition to prove the convergence of error dynamics. V ˙ ( t ) < 0 is equivalent to
A e T P + P A e 0 .
Equation (19) is the strict LMI that can be computed with MATLAB. Formula (19) can yield that A e is a Hurwitz matrix. By configuring matrix A e , the convergence of the system can be reached.
However, the above details are mainly about the integer-order system. There are no corresponding lemmas about the fractional-order system. Therefore, the fractional-order extension of a Lyapunov candidate function is introduced to demonstrate the convergence of the error dynamics by referring [37].
Lemma 4.
Consider error dynamics e ( t ) R n to be a continuous and derivable function. Then, for any time t 0 , the fractional derivative of the Lyapunov function is
D t α V ( t ) ( D t α e T ( t ) ) P e ( t ) + e T ( t ) P ( D t α e ( t ) ) ,
where V ( t ) = e T ( t ) P e ( t ) is the Lyapunov function connected with e ( t ) .
Proof. 
Denote that J = ( D t α e T ) P e + e T P ( D t α e ) D t α V .
According to the definition of fractional calculus, D t α V ( t ) is equivalent to
D t α V ( t ) = 1 Γ ( 1 α ) 0 t V ˙ ( τ ) ( t τ ) α d τ = 1 Γ ( 1 α ) 0 t e T ( τ ) P e ˙ ( τ ) + e ˙ T ( τ ) P e ( τ ) ( t τ ) α d τ .
Considering (21), J is rewritten as
J = 1 Γ ( 1 α ) 0 t e T ( t ) P e ˙ ( τ ) + e ˙ T ( τ ) P e ( t ) e T ( τ ) P e ˙ ( τ ) e ˙ T ( τ ) P e ( τ ) ( t τ ) α d τ = 1 Γ ( 1 α ) 0 t ( e T ( t ) e T ( τ ) ) P e ˙ ( τ ) + e ˙ T ( τ ) P ( e ( t ) e ( τ ) ) ( t τ ) α d τ = 1 Γ ( 1 α ) 0 t z ˙ T ( τ ) P z ( τ ) + z T ( τ ) P z ˙ ( τ ) ( t τ ) α d τ ,
where z ( τ ) = e ( t ) e ( τ ) .
Then, the proof for Lemma 1 is equivalent to proving that
1 Γ ( 1 α ) 0 t ( z ˙ T ( τ ) ) P z ( τ ) + z T ( τ ) P z ˙ ( τ ) ( t τ ) α d τ 0 .
For u = z T ( t ) P z ( t ) , one can deduce that
d u d t = z ˙ T P z ( t ) + z T ( t ) P z ˙ ( t ) .
Then, v ( τ ) = ( t τ ) α Γ ( 1 α ) is denoted, which yields
d v ( τ ) = α ( t τ ) α 1 Γ ( 1 α ) d τ .
On the basis of (24) and (25), (23) is equal to
z T ( τ ) P z ( τ ) Γ ( 1 α ) ( t τ ) α 0 t α Γ ( 1 α ) 0 t z T ( τ ) P z ( τ ) ( t τ ) α d τ 0 ,
with
z T ( τ ) P z ( τ ) Γ ( 1 α ) ( t τ ) α 0 t = z T ( τ ) P z ( τ ) Γ ( 1 α ) ( t τ ) α τ = t z T ( 0 ) P z ( 0 ) Γ ( 1 α ) ( t ) α .
When τ = t , the first term of (27), has nondeterminacy, and the analysis of its limitation is necessary:
lim τ t z T ( τ ) P z ( τ ) Γ ( 1 α ) ( t τ ) α = e ˙ T ( τ ) P e ( t ) + e ˙ T ( τ ) P e ( τ ) e T ( t ) P e ˙ ( τ ) + e T ( τ ) P e ˙ ( τ ) Γ ( 1 α ) α ( t τ ) α 1 = e ˙ T ( τ ) P e ( t ) + e ˙ T ( τ ) P e ( τ ) e T ( t ) P e ˙ ( τ ) + e T ( τ ) P e ˙ ( τ ) Γ ( 1 α ) α ( t τ ) 1 α .
Due to α ( 0 , 1 ) , it is obvious that z T ( τ ) P z ( τ ) Γ ( 1 α ) ( t τ ) α 0 when τ = t .
Expression (26) can be simplified to
z T ( 0 ) P z ( 0 ) Γ ( 1 α ) t α α Γ ( 1 α ) 0 t z T ( τ ) P z ( τ ) ( t τ ) α d τ 0 .
For P 0 and Γ ( · ) > 0 , it is certain to stand.    □
Theorem 1.
For System (15), there exist bounds x ¯ and x ̲ , such that x ̲ ( t ) x ( t ) x ¯ ( t ) if the following conditions are satisfied:
1.
Ω = I N ( A L C ) γ ( L M ) is Metzler;
2.
The initial condition of (7) satisfies x ̲ ( 0 ) x ( 0 ) x ¯ ( 0 ) ;
3.
Nonlinear function f ( x i ( t ) ) possesses the features introduced in Properties 3–5.
Proof. 
From Properties 3 and 4, the upper bound of the nonlinear function is equal to
f ˜ ( x ¯ i , x ̲ i ) = a ( x ¯ i ) b ( x ̲ i ) .
Then, it follows from (30) that
f ˜ ( x ¯ i , x ̲ i ) f ˜ ( x i , x i ) = ( a ( x ¯ i ) a ( x i ) ) ( b ( x ̲ i ) b ( x i ) ) .
Functions a ( · ) and b ( · ) are all increasing, which yields
f ˜ ( x ¯ i , x ̲ i ) f ˜ ( x i , x i ) 0 .
Similarly, one can deduce that
f ˜ ( x i , x i ) f ˜ ( x ̲ i , x ¯ i ) 0 .
From x ̲ ( 0 ) x ( 0 ) x ¯ ( 0 ) , the initial value of the error dynamics satisfies e ¯ ( 0 ) 0 and e ̲ ( 0 ) 0 . If matrix Ω is Metzler, and (32) and (33) are true, it follows that e ¯ ( t ) 0 and e ̲ ( t ) 0 , i.e.,  x ̲ ( t ) x ( t ) x ¯ ( t ) , for any t 0 .    □
Remark 2.
For Theorem 1, we constructed a frame for (7). The positivity of the error dynamics is guaranteed, which implies x ̲ ( t ) x ( t ) x ¯ ( t ) . However, an interval observer not only requires the error dynamics be positive, but also that the upper and lower errors are in a convergence of zero. Then, we give the following theorem.
Theorem 2.
On the basis of the result in Theorem 1, given constant τ > 0 and positive matrix P = P T , if there exists a solution such that
Π ^ = H e ( P A U C + P N 1 ) 2 τ I P N 2 + N 3 T P P N 3 + N 2 T P H e ( P A U C + P N 4 ) 2 τ I 0 ,
γ > τ a ( L ) ,
where γ is the coupling strength, L = P 1 U and M = P 1 are the observer gains, then System (15) is a distributed interval observer.
Proof. 
The Lyapunov function is chosen as follows:
V ( t ) = i = 1 N r i e ¯ i T ( t ) P e ¯ i ( t ) + i = 1 N r i e ̲ i T ( t ) P e ̲ i ( t ) .
From Lemma 1, the fractional derivative of V ( t ) is
D t α V ( t ) i = 1 N r i ( D t α e ¯ i T ) ( t ) P e ¯ i ( t ) + i = 1 N r i e ̲ i T ( t ) P ( D t α e ̲ i ( t ) ) + i = 1 N r i ( D t α e ̲ i T ) ( t ) P e ̲ i ( t ) + i = 1 N r i e ¯ i T ( t ) P ( D t α e ¯ i ( t ) ) .
Substituting the error dynamics into (35)
D t α V ( t ) = i = 1 N r i ( Ω i e ¯ i + Δ f ¯ ( x i ) ) T P e ¯ i + i = 1 N r i e ̲ i T P ( Ω i e ̲ i + Δ f ̲ ( x i ) ) + i = 1 N r i ( Ω i e ̲ i + Δ f ̲ ( x i ) ) T P e ¯ i + i = 1 N r i e ¯ i T P ( Ω i e ¯ i + Δ f ¯ ( x i ) ) i = 1 N r i ( Ω i e ¯ i + F 1 e ¯ i + F 2 e ̲ i ) T P e ¯ i + i = 1 N r i e ¯ i T P ( Ω i e ¯ i + F 1 e ¯ i + F 2 e ̲ i ) + i = 1 N r i e ̲ i T P ( Ω i e ¯ i + F 3 e ¯ i + N F e ̲ i ) + i = 1 N r i ( Ω i e ¯ i + F 3 e ¯ i + F 4 e ̲ i ) T P e ̲ i = ( ( R ( A L C ) γ R L M ) e ¯ + R F 1 e ¯ + R F 2 e ̲ ) T ( I P ) e ¯ + e ¯ T ( I P ) ( ( R ( A L C ) γ R L M ) e ¯ + R F 1 e ¯ + R F 2 e ̲ + ( ( R ( A L C ) γ R L M ) e ̲ + R N 3 e ¯ + R F 4 e ̲ ) T ( I P ) e ¯ + e ¯ T ( I P ) ( ( R ( A L C ) γ R L M ) e ¯ + R F 3 e ¯ + R N 4 e ̲ ) = e ¯ T ( R ( H e ( P A P L C + P F 1 ) ) γ ( R L + L T R ) P M ) e ¯ + e ̲ T ( R ( H e ( P A P L C + P F 4 ) ) γ ( R L + L T R ) P M ) e ̲ + e ¯ T ( P F 2 + F 3 T P ) e ̲ + e ̲ T ( P F 3 + F 2 T P ) e ¯ ,
where Ω i = A L C γ M j = 1 N L i j , Δ f ¯ ( x i ) = f ˜ ( x ¯ i , x ̲ i ) f ˜ ( x i , x i ) and Δ f ̲ ( x i ) = f ˜ ( x i , x i ) f ˜ ( x ̲ i , x ¯ i ) .
According to Lemma 3,  R L + L T R could be simplified as
2 a ( L ) e ¯ T R e ¯ e ¯ T ( R L + L T R ) e ¯ ,
2 a ( L ) e ̲ T R e ̲ e ̲ T ( R L + L T R ) e ̲ .
Taking M = P 1 , (37) and (38) into account, it follows from (36) that
D t α V ( t ) e ¯ T ( R ( H e ( P A P L C + P F 1 ) ) 2 γ a ( L ) R I n ) e ¯ + e ̲ T ( R ( H e ( P A P L C + P F 4 ) ) 2 γ a ( L ) R I n e ̲ + e ¯ T ( P F 2 + F 3 T P ) e ̲ + e ̲ T ( P F 3 + F 2 T P ) e ¯ .
Considering γ > τ a ( L ) , we have
D t α V ( t ) e ¯ T ( R ( H e ( P A P L C + P F 1 ) ) 2 τ R I n e ¯ + e ̲ T ( R ( H e ( P A P L C + P F 4 ) ) 2 τ R I n e ̲ + e ¯ T ( P F 2 + F 3 T P ) e ̲ + e ̲ T ( P F 3 + F 2 T P ) e ¯ = ϵ T ( t ) R Π ϵ ( t ) ,
where ϵ ( t ) = [ e ¯ T ( t ) , e ̲ T ( t ) ] T and
Π = H e ( P A P L C + P F 1 ) 2 τ I n P F 2 + F 3 T P P F 3 + F 2 T P H e ( P A P L C + P F 4 ) 2 τ I n .
To satisfy the LMI toolbox, U = P L is applied in Π , which leads to
Π ^ = H e ( P A U C + P F 1 ) 2 τ I n P F 2 + F 3 T P P F 3 + F 2 T P H e ( P A U C + P F 4 ) 2 τ I n .
Then, matrix Π ^ 0 is equal to D t α V ( t ) < 0 , which implies that lim t e ¯ ( t ) = 0 and lim t e ̲ ( t ) = 0 . The boundedness of the error dynamics could be guaranteed.    □
Remark 3.
Constant τ is simple without an actual effect. However, it is a parameter connected with γ. If we just compute γ, the LMI toolbox just gives one feasible solution. Nevertheless, if we compute τ, γ > τ a ( L ) would have more regions to select.
On the basis of Theorem 2, an algorithm was constructed to design distributed interval observers for fractional-order MASs.

4. Simulation

Considering a fractional-order MAS with nonlinearity, the state-space model is similar to (7), where
A = 0.5 2 0 1 , C = 0 1 , f ( x i ) = s i n ( x i ) .
For nonlinear function f ( x i ) , corresponding function f ˜ ( x ¯ i , x ̲ i ) is defined as
f ˜ ( x ¯ i , x ̲ i ) = s i n ( x ¯ i ) + x ¯ i x ̲ i .
For functions a ( x i ) = s i n ( x i ) + x i and b ( x i ) = x i , it is obvious that a ( x i ) and b ( x i ) are increasing functions. a ( x i ) = s i n ( x i ) + x i and b ( x i ) = x i are substituted into (13); then,
s i n ( x ¯ i ) + x ¯ i x ̲ i s i n ( x i ) = ( s i n ( x ¯ i ) s i n ( x i ) ) + ( x ¯ i x i ) + ( x i x ̲ i ) 2 ( x ¯ i x i ) + ( x i x ̲ i ) = F 1 e ¯ i + F 2 e ̲ i .
We chose F 1 = 2 I and F 2 = I . Similarly, we have F 3 = I and F 4 = 2 I .
Laplacian matrix L is
L = 3 1 1 1 1 2 1 0 1 1 2 0 1 0 0 1 .
By using Lemma 3, we obtained generalized algebraic connectivity a ( L ) = 1 . By solving the LMI in Theorem 2, the observer gains can be computed:
L = 1.1771 0.2282 , M = 0.1906 0.0757 0.0757 0.1470 , τ = 2.9151 .
Then we chose γ = 3 . The initial value of the original system state is defined as x ( 0 ) = [ 1 2 3 5 2 1 1 4 ] T . The initial value of the distributed interval observer is defined as x ¯ ( 0 ) = [ 6 7 8 10 7 4 4 9 ] T and x ̲ ( 0 ) = [ 4 3 2 0 4 6 6 1 ] T .
By performing Steps 5–8 in Algorithm 1, we can obtain Figure 1, Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6. Then, Figure 1 and Figure 2 show the original system state of the four agents. Figure 3 and Figure 4 display the bounds from a distributed interval observer. v i j means the upper bound of the jth state of the ith agent, while u i j means the lower bound of the jth state of the ith agent. From Figure 3 and Figure 4, the bounds of the distributed observer trap the state of the original system. Define that e i j + = x ¯ i j x i j and e i j + = x i j x ̲ i j . For Figure 5 and Figure 6, the error between the observer and the original system could be reduced to a bounded value, which implies that the distributed interval observer is feasible.
Algorithm 1 Distributed interval estimation for fractional-order MASs.
Step
1: Given the constant matrix L , compute the generalized algebraic connectivity a ( L ) ;
Step
2: For given matrices A R n × n , C R m × n , compute the observer gains L and M and the constant τ from Theorem 2;
Step
3: Select appropriate constant γ > τ a ( L ) ;
Step
4: Ensure that the matrix Ω is a Metzler matrix;
Step
5: For the total time t = 10 , choosing the step size h = 0.001 , the step N = T / h can be calculated.
Step
6: Based on Step 5, construct two loops. The first loop is for the original fractional-order system and the other loop is for the distributed interval observer.
Step
7: Construct an array with N + 1 volumes to store the output of the original system in step 6.
Step
8: Establish a distributed interval observer based on the output in Step 7 and the other loop proposed in Step 6.

5. Conclusions

In this paper, a distributed interval observer design methodology for linear fractional-order MASs with nonlinearity was proposed. A Lyapunov method that is useful for observer and controller design was first introduced for general fractional-order systems. For MASs, graph theory was applied to fractional-order systems, and the strict LMI and an effective algorithm were presented for observer design. Lastly, an example was given to demonstrate the effectiveness of the proposed method. In the future, by using the H technique, we aim to focus on research regarding the consensus control or formation control of fractional-order MASs.

Author Contributions

Conceptualization, H.Z. and J.H.; methodology, J.H.; software, H.Z.; validation, H.Z., J.H. and S.H.; formal analysis, H.Z.; investigation, S.H.; resources, H.Z.; data curation, H.Z.; writing—original draft preparation, H.Z.; writing—review and editing, H.Z.; visualization, H.Z.; supervision, J.H.; project administration, J.H.; funding acquisition, J.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research were funded by the Natural Science Foundation of Jiangsu province of China under grant BK2021-1309, and in part by the Open Fund for Jiangsu Key Laboratory of Advanced Manufacturing Technology under grant HGAMTL-2101 and by Undergraduate Training Program For Innovation and Entrepreneurship, Soochow University (202210285152Y).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to express our great appreciation to the editors and reviewers.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hartley, T.; Lorenzo, C. Fractional-Order system identification based on continuous order-distributions. Signal Process. 2003, 83, 2287–2300. [Google Scholar] [CrossRef]
  2. Li, M.; Li, D.; Wang, J.; Zhao, C. Active disturbance rejection control for fractional-order system. ISA Trans. 2013, 52, 365–374. [Google Scholar] [CrossRef] [PubMed]
  3. Lu, B.; Zhang, Y.; Reeves, D.M.; Sun, H.; Zheng, C. Application of tempered-stable time fractional-derivative model to upscale subdiffusion for pollutant transport in field-scale discrete fracture networks. Mathematics 2018, 6, 5. [Google Scholar] [CrossRef] [Green Version]
  4. Turab, A.; Mitrović, Z.; Savić, A. Existence of solutions for a class of nonlinear boundary value problems on the hexasilinane graph. Adv. Differ. Equ. 2021; to be published. [Google Scholar] [CrossRef]
  5. Rezapour, S.; Souid, M.S.; Bouazza, Z.; Hussain, A.; Etemad, S. On the fractional variable order thermostat model: Existence theory on cones via piece-wise constant functions. J. Funct. Space, 2022; to be published. [Google Scholar] [CrossRef]
  6. Zhang, K.; Jiang, B.; Yan, X.; Shen, J. Interval sliding mode observer based incipient sensor fault detection with application to a traction device in China railway high-speed. IEEE Trans. Veh Technol. 2019, 68, 2585–2597. [Google Scholar] [CrossRef]
  7. Huang, J.; Ma, X.; Che, H.; Han, Z. Further result on interval observer design for discrete-time switched systems and application to circuit systems. IEEE Trans. Circuits Syst. II-Express Briefs 2019, 67, 2542–2546. [Google Scholar] [CrossRef]
  8. Gouzé, J.; Rapaport, A.; Hadj-Sadok, M. Interval observers for uncertain biological systems. Ecol. Modell. 2000, 133, 45–56. [Google Scholar] [CrossRef]
  9. Mazenc, F.; Bernard, O. Interval observers for linear time-invariant systems with disturbances. Automatica 2011, 47, 140–147. [Google Scholar] [CrossRef] [Green Version]
  10. Raïssi, T.; Efimov, D.; Zolghadri, A. Interval state estimation for a class of nonlinear systems. IEEE Trans. Autom. Control 2011, 57, 260–265. [Google Scholar] [CrossRef]
  11. Dinh, T.N.; Marouani, G.; Raïssi, T.; Wang, Z.; Messaoud, H. Optimal interval observers for discrete-time linear switched systems. Int. J. Control 2020, 93, 2613–2621. [Google Scholar] [CrossRef]
  12. Huang, J.; Che, H.; Raïssi, T.; Wang, Z. Functional interval observer for discrete-time switched descriptor systems. IEEE Trans. Autom. Control 2021, 67, 2497–2504. [Google Scholar] [CrossRef]
  13. Xu, F.; Yang, S.; Wang, X. A novel set-theoretic interval observer for discrete linear time-invariant systems. IEEE Trans. Autom. Control 2021, 66, 773–780. [Google Scholar] [CrossRef]
  14. Cai, H.; Lewis, F.L.; Hu, G.; Huang, J. The adaptive distributed observer approach to the cooperative output regulation of linear multi-agent systems. Automatica 2017, 75, 299–305. [Google Scholar] [CrossRef]
  15. Han, W.; Trentelman, H.; Wang, Z.; Shen, Y. A simple approach to distributed observer design for linear systems. IEEE Trans. Autom. Control 2018, 64, 329–336. [Google Scholar] [CrossRef] [Green Version]
  16. Huang, J.; Yang, M.; Zhang, Y.; Zhang, M. Consensus control of multi-agent systems with P-one-sided Lipschitz. ISA Trans. 2022, 125, 42–49. [Google Scholar] [CrossRef]
  17. Li, C.; Yu, X.; Liu, Z.W.; Huang, T. Asynchronous impulsive containment control in switched multi-agent systems. Inf. Sci. 2016, 370, 667–679. [Google Scholar] [CrossRef]
  18. Hua, C.; You, X.; Guan, X. Leader-Following consensus for a class of high-order nonlinear multi-agent systems. Automatica 2016, 73, 138–144. [Google Scholar] [CrossRef]
  19. Chen, D.; Liu, G. A networked predictive controller for linear multi-agent systems with communication time delays. J. Frankl. Inst. 2020, 357, 9442–9466. [Google Scholar] [CrossRef]
  20. Hong, Y.; Chen, G.; Bushnell, L. Distributed observers design for leader-following control of multi-agent networks. Automatica 2008, 44, 846–850. [Google Scholar] [CrossRef]
  21. Liu, T.; Huang, J. A distributed observer for a class of nonlinear systems and its application to a leader-following consensus problem. IEEE Trans. Autom. Control 2018, 44, 1221–1227. [Google Scholar] [CrossRef]
  22. Yu, Z.; Jiang, H.; Hu, C. Leader-Following consensus of fractional-order multi-agent systems under fixed topology. Neurocomputing 2015, 149, 613–620. [Google Scholar] [CrossRef]
  23. Li, Z.; Gao, L.; Chen, W.; Xu, Y. Distributed adaptive cooperative tracking of uncertain nonlinear fractional-order multi-agent systems. IEEE-CAA J. Automatic. 2016, 47, 222–234. [Google Scholar] [CrossRef]
  24. Zhu, W.; Li, W.; Zhou, P.; Yang, C. Consensus of fractional-order multi-agent systems with linear models via observer-type protocol. Neurocomputing 2017, 230, 60–65. [Google Scholar] [CrossRef]
  25. Gong, Y.; Wen, G.; Peng, Z.; Huang, T.; Chen, Y. Observer-Based time-varying formation control of fractional-order multi-agent systems with general linear dynamics. IEEE Trans. Circuits Syst. II-Express Briefs 2019, 67, 82–86. [Google Scholar] [CrossRef]
  26. Wen, G.; Zhang, Y.; Peng, Z.; Yu, Y.; Rahmani, A. Observer-Based output consensus of leader-following fractional-order heterogeneous nonlinear multi-agent systems. Int. J. Control 2020, 93, 2516–2524. [Google Scholar] [CrossRef]
  27. Afaghi, A.; Ghaemi, S.; Ghiasi, A.; Badamchizadeh, M. Adaptive fuzzy observer-based cooperative control of unknown fractional-order multi-agent systems with uncertain dynamics. Soft Comput. 2020, 24, 3737–3752. [Google Scholar] [CrossRef]
  28. Danca, M.; Kuznetsov, N. Matlab code for Lyapunov exponents of fractional-order systems. Int. J. Bifurc. Chaos 2018, 28, 1850067. [Google Scholar] [CrossRef] [Green Version]
  29. Tian, X.; Yang, Z. Adaptive stabilization of a fractional-order system with unknown disturbance and nonlinear input via a backstepping control technique. Symmetry 2019, 12, 55. [Google Scholar] [CrossRef] [Green Version]
  30. Huong, D. Design of functional interval observers for non-linear fractional-order systems. Asian J. Control 2020, 22, 1127–1137. [Google Scholar] [CrossRef]
  31. Podlubny, I. Fractional derivatives and integrals. In Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications; Elsevier: San Diego, CA, USA, 1998; pp. 41–80. [Google Scholar]
  32. Tavazoei, M.; Haeri, M. A note on the stability of fractional order systems. Math. Comput. Simul. 2009, 79, 1566–1576. [Google Scholar] [CrossRef]
  33. Kamal, S.; Sharma, R.K.; Dinh, T.N.; Ms, H.; Bandyopadhyay, B. Sliding mode control of uncertain fractional-rder systems: A reaching phase free approach. Asian J. Control 2021, 23, 199–208. [Google Scholar] [CrossRef]
  34. Ren, W.; Beard, R. Consensus seeking in multiagent systems under dynamically changing interaction topologies. IEEE Trans. Autom. Control 2005, 50, 655–661. [Google Scholar] [CrossRef]
  35. Yu, W.; Chen, G.; Cao, M.; Kurths, J. Second-Order consensus for multiagent systems with directed topologies and nonlinear dynamics. IEEE Trans. Syst. Man Cybern. Syst. 2009, 40, 881–891. [Google Scholar]
  36. Moisan, M.; Bernard, O. Robust interval observers for global Lipschitz uncertain chaotic systems. Syst. Control Lett. 2010, 59, 687–694. [Google Scholar] [CrossRef]
  37. Aguila-Camacho, N.; Duarte-Mermoud, M.; Gallegos, J. Lyapunov functions for fractional order systems. Commun. Nonlinear Sci. Numer. Simul. 2014, 19, 2951–2957. [Google Scholar] [CrossRef]
Figure 1. The first state of the four-agents.
Figure 1. The first state of the four-agents.
Fractalfract 06 00355 g001
Figure 2. The second state of the four-agents.
Figure 2. The second state of the four-agents.
Fractalfract 06 00355 g002
Figure 3. The bounds of the first state for the four-agents.
Figure 3. The bounds of the first state for the four-agents.
Fractalfract 06 00355 g003
Figure 4. The bounds of the second state for the four-agents.
Figure 4. The bounds of the second state for the four-agents.
Fractalfract 06 00355 g004
Figure 5. The error of the first state for the four-agents.
Figure 5. The error of the first state for the four-agents.
Fractalfract 06 00355 g005
Figure 6. The error of the second state for the four-agents.
Figure 6. The error of the second state for the four-agents.
Fractalfract 06 00355 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, H.; Huang, J.; He, S. Fractional-Order Interval Observer for Multiagent Nonlinear Systems. Fractal Fract. 2022, 6, 355. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract6070355

AMA Style

Zhang H, Huang J, He S. Fractional-Order Interval Observer for Multiagent Nonlinear Systems. Fractal and Fractional. 2022; 6(7):355. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract6070355

Chicago/Turabian Style

Zhang, Haoran, Jun Huang, and Siyuan He. 2022. "Fractional-Order Interval Observer for Multiagent Nonlinear Systems" Fractal and Fractional 6, no. 7: 355. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract6070355

Article Metrics

Back to TopTop