Next Article in Journal
GSRec: A Graph-Sequence Recommendation System Based on Reverse-Order Graph and User Embedding
Next Article in Special Issue
Existence and Hyers–Ulam Stability of Stochastic Delay Systems Governed by the Rosenblatt Process
Previous Article in Journal
Empirical-Likelihood-Based Inference for Partially Linear Models
Previous Article in Special Issue
Combination of Functional and Disturbance Observer for Positive Systems with Disturbances
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Distributed Interval Observers with Switching Topology Design for Cyber-Physical Systems

The School of Mechanical and Electrical Engineering, Soochow University, Suzhou 215325, China
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(1), 163; https://doi.org/10.3390/math12010163
Submission received: 30 November 2023 / Revised: 25 December 2023 / Accepted: 29 December 2023 / Published: 4 January 2024
(This article belongs to the Special Issue Dynamical System and Stochastic Analysis)

Abstract

:
In this paper, the distributed interval estimation problem for networked Cyber-Physical systems suffering from both disturbances and noise is investigated. In the distributed interval observers, there are some connected interval observers built for the corresponding subsystems. Then, due to the communication burden in Cyber-Physical systems, we consider the case where the communication among distributed interval observers is switching topology. A novel approach that combines L methodology with reachable set analysis is proposed to design distributed interval observers. Finally, the performance of the proposed distributed interval observers with switching topology is verified through a simulation example.

1. Introduction

Cyber-Physical Systems (CPSs) are the combinations of physical procedures, high-efficiency computation, communication, and effective control defined by [1]. Architecturally, from [2], a typical CPS can be divided into three layers, which are composed of the sensing layer, the network layer, and the control layer. The development of distributed sensing and networking technologies such as [3,4] has enabled omnipresent sensing and computing capabilities. This has led to the implementation of CPSs in large-scale networks. CPSs are widely used in industrial informatics [5] manufacturing [6], healthcare [7], electrical grids [8], and so on. State estimation and observer design are crucial research areas in CPSs. Ref. [9] used a sliding mode observer and integrated the event-triggered mechanism to estimate the state of CPSs from sensor measurements. Ref. [10] introduced a security estimator combined with a Kalman filter to improve the practical performance of state estimation for CPSs. Ref. [11] accomplished state estimation and resilient control of CPSs using finite time observer techniques and switching schemes. It should be noted that the state estimation for CPSs with bounded disturbance and noise has not been investigated sufficiently.
On the other hand, disturbances and noise always exist in real systems, and the interval observer serves as a powerful estimator of upper and lower bounds for uncertain systems with disturbances and noise. In [12], the concept, as well as the framework of interval observer, were presented. Using the monotone system theory, Refs. [13,14] proposed the approach of coordinate transformation that serves as an efficient strategy to reduce the strict conditions for interval observer design. In recent years, the set-membership estimation method was applied effectively in interval observer design. A two-step interval observer design methodology that combines reachability analysis with robust observer was first presented in [15]. In [15], the H method and reachable set analysis are combined to design interval observers that eliminate the effects of perturbations and noise on the system. At the same time, the L (or L 2 ) method is also a powerful, robust property, and it is widely used in control and observation fields, such as [16,17,18]. It has been recently shown by [19,20,21] that the reachability analysis estimation method can not only enhance the accuracy of estimation but also increase the design freedom. Concurrently, in the context of distributed systems, several recent studies have been conducted on distributed interval observers, such as [22,23,24]. Ref. [22] designed distributed interval observers based on the monotone system approach for multiagent systems. Ref. [23] considered a distributed interval observers design problem for a class of linear time-invariant systems with uncertainty. At the same time, Ref. [24] improved distributed interval observers by using a set-membership estimation approach. However, the topology of interval observer in the aforementioned work [22,23,24] is fixed, and it is usually switching with respect to time in practice.
In light of the above discussion, for the problem of state estimation for uncertain CPSs, we apply distributed interval observers to such systems. Since the state estimation of CPSs with bounded disturbances and noise has not been investigated sufficiently, it is meaningful to study the state estimation problem for uncertain CPSs. Each interval observer has two types of observer gain: one is obtained by using the traditional observer design method, and the other one is determined by employing neighborhood information. Considering the practical communication problem of the network layer in CPSs, we suppose that the communication among distributed interval observers is described by switching topology. There are three main challenges: one is to design L observers with optimal performance for networked CPSs, another is to construct a reachable set analysis framework for CPSs, which then gives upper and lower bounds on the state of the system, and the last is to solve the switching topology problem among the observers. The contributions of this paper are summarized in aspects below.
(1)
A distributed interval observer methodology for CPSs is proposed. Compared with the monotone system method, the estimation accuracy is greatly improved by using the two-step method. The L technique is used to deal with the effects of uncertainty in observer design.
(2)
The switching topology with average dwell time (ADT) among distributed interval observers is taken into account and is more closely aligned with the actual system. It can also reduce the communication burden of CPSs.
This paper is structured as follows, with the rest of the paper starting with the graph theory, system model, and some basics presented in Section 2. In Section 3, the optimal robust observer is designed using the L technique. In order to complete the interval estimation, a reachability analysis methodology is used to design the distributed interval observer. In Section 4, the paper simulates a networked CPS with four Unmanned Aerial Vehicles (UAV) models to illustrate the effectiveness of the distributed interval observer. Finally, Section 5 concludes the paper.
Notation: For a matrix of real symmetry E R N × N , E 0 demonstrates that E is positive definite, while E 0 demonstrates that E is negative and H e ( E ) = E + E T . The maximal (minimum) eigenvalue of the matrix Q is denoted by λ m a x ( Q ) ( λ m i n ( Q ) ) . The norm L 2 of the vector v is represented by v 2 . In other words, v 2 = v T v . Similarly, the norm L of the vector v is denoted by v . The symbol * means the term that can arise due to symmetry in the symmetric matrices.

2. Preliminaries

2.1. Graph Theory

For a digraph G which has N vertices, A = [ a i j ] R N × N is called the adjacency matrix. The weight associated with the edge ( i , j ) S that connects node i to node j is called a i j , and A is provided by a i j = 0 . The length path from vertex i to vertex j is made up of t + 1 distinct vertices with consecutive vertices adjacent to each other. If there is a path connecting any two vertices of the graph G , the graph G is considered connected. L = D A is defined as Laplacian matrix of a graph G . D is referred to as a degree matrix of G . The Laplacian matrix L of a connected graph has a single zero eigenvalue, and 1 N is the associated eigenvector. Furthermore, 0 = λ 1 ( G ) λ N ( G ) if G is connected, where λ i ( G ) ( i = 1 , 2 , , N ) is the eigenvalue of a Laplacian matrix L .
Lemma 1
([25]). For a strongly connected graph G , denote t i ( i = 1 , 2 , , N ) as the left eigenvector with a 0 eigenvalue, and T = d i a g { t 1 , t 2 , , t N } , then, we can obtain T L + L T T 0 .
Lemma 2
([26]). Suppose that graph G is a strongly connected and balanced graph, and the algebraic connectivity of G is defined by a ( L ) = m i n t T x = 0 , x 0 x T ( T L + L T T ) x 2 x T T x where T = d i a g { t 1 , t 2 , , t N } . Then, we can obtain a ( L ) = λ m i n ( H e ( L ) 2 ) .

2.2. System Model

For a networked CPS, consider a network with N subsystems. The following is the ith subsystem with disturbances and noise:
x i ( k + 1 ) = A x i ( k ) + B u i ( k ) + p i ( k ) , y i ( k ) = C x i ( k ) + q i ( k ) ,
where x i ( k ) R n is the state, u i ( k ) R n is the control input, y i ( k ) R m is the output, p i ( k ) R n is the disturbance, q i ( k ) R m is the noise. A R n × n , B R n × n and C R m × n are constant matrices.
In the following, it is supposed that the communication topology of the subsystem is strongly connected. The global dynamic system of (1) can be given:
x ( k + 1 ) = A ¯ x ( k ) + B ¯ u ( k ) + p ( k ) , y ( k ) = C ¯ x ( k ) + q ( k ) ,
where x = [ x 1 T , , x N T ] T , u = [ u 1 T , , u N T ] T , y = [ y 1 T , , y N T ] T , q = [ q 1 T , , q N T ] T , p = [ p 1 T , , p N T ] T , A ¯ = d i a g { A , , A N } , B ¯ = d i a g { B , , B N } and C ¯ = d i a g { C , , C N } .
Since the information can be received by a single subsystem from its neighborhood, we consider the case of switching topology of the observer system and then present ϕ ( k ) , a stepwise constant function that takes values from the finite collection S = 1 , 2 , , N . The observer dynamics of ith subsystem are:
x ^ i ( k + 1 ) = A x ^ i ( k ) + B u i ( k ) + L i ( y i ( k ) C x ^ i ( k ) ) + κ ϕ ( k ) M i j = 1 N a i j ϕ ( k ) ( x ^ j ( k ) x ^ i ( k ) ) ,
where κ ϕ ( k ) is the coupled gain that needs to be designed and M i and L i are observer gains of the ith subsystem, and a i j ϕ ( k ) represents the connectivity weight from subsystem i to subsystem j at moment k.
By subtracting (1) from (3), we obtain the error dynamics of a single subsystem:
e i ( k ) = x i ( k ) x ^ i ( k ) , e i ( k + 1 ) = ( A L i C ) e i ( k ) κ ϕ ( k ) M i j = 1 N a i j ϕ ( k ) ( x ^ j ( k ) x ^ i ( k ) ) + D i d i ( k ) ,
with D i = I L i and d i ( k ) = p i ( k ) q i ( k ) .
Then we can obtain the dynamic system of the global observer
x ^ ( k + 1 ) = A ¯ x ^ ( k ) + B ¯ u ( k ) + L ¯ ϕ ( k ) ( y ( k ) C ¯ x ^ ( k ) ) + κ ϕ ( k ) M ¯ ϕ ( k ) ( L ϕ ( k ) I N ) ( x ( k ) x ^ ( k ) ) ,
with L ¯ = d i a g { L 1 , , L N N } and M ¯ = d i a g { M 1 , , M N N } .
The global error system e ( k ) = [ e 1 T ( k ) , , e N T ( k ) ] T is as follows
e ( k + 1 ) = ( A ¯ L ¯ ϕ ( k ) C ¯ κ ϕ ( k ) M ¯ ϕ ( k ) ( L ϕ ( k ) I N ) ) e ( k ) + D ¯ d ( k ) ,
and the compact form of (6) can be writen as
e ( k + 1 ) = Γ ϕ ( k ) e ( k ) + D ¯ d ( k ) ,
where D ¯ = d i a g { D 1 , , D N N } and Γ ϕ ( k ) = A ¯ L ¯ ϕ ( k ) C ¯ κ ϕ ( k ) M ¯ ϕ ( k ) ( L ϕ ( k ) I N ) .
Definition 1
([27]). If the following condition holds,
e ( k ) 2 ϰ d 2 + V ( 0 ) ω k ,
where ϰ > 0 and 0 < ω < 1 , V ( 0 ) = e T ( 0 ) P ϕ ( k ) e ( 0 ) and P ϕ ( k ) 0 . Then the observer (5) is a L robust observer for system (2).
Definition 2
([28]). Let N ψ ( k 1 , k 2 ) be the switching times of ψ ( k ) across the range [ k 1 , k 2 ) . If
N ψ ( k 1 , k 2 ) N 0 + k 2 k 1 τ ˜ ,
for given τ ˜ > 0 and N 0 0 , τ ˜ is the average dwell time (ADT) of the switching signal ψ ( k ) . In this paper, we let N 0 = 0 .
Definition 3
([29]). The definition of an α-dimensional zonotope is as follow
Ω = ν H Y α = ν + H z , z Y α ,
where ν R ι represents a given vector, H R ι × α represents a given matrix, Y α is a unitary box made up of α unitary intervals and Y α is a unitary interval. In the sequel, the zonotope Ω is represented as ν , H for the sake of simplicity.
Lemma 3
([30]). The following equation is satisfied for a zonotope defined in (10):
ν 1 , H 1 ν 2 , H 2 = ν 1 + ν 2 , [ H 1 , H 2 ] ,
W ν , H = W ν , W H ,
ν , H ν , H ¯ ,
where H 1 and H 2 represent the shape matrices of each of the zonotopes, and ν 1 and ν 2 are their centers. H ¯ R N × N means a diagonal matrix with H ¯ M , M = i = 1 α | H M , i | , M = 1 , . . . , ρ .
Remark 1.
H ¯ can be expressed in the form shown below:
H ¯ = i = 1 α | H 1 , i | 0 0 i = 1 α | H ρ , i | .
If α > ρ , then ν , H ν , H ¯ is applied to reduce the order of high order zonotopes.
Remark 2.
For zonotopes Ω M R N , M = 1 , . . . , ρ , the Minkowski sum of them is
M = 1 ρ Ω M = Ω 1 Ω 2 Ω ρ .
Definition 4
([19]). For an α-order zonotope, there is an interval hull Ω that could contain Ω in its entirety:
Ω B o x ( Ω ) = [ a , b ] ,
with a = [ a 1 , , a α ] T , b = [ b 1 , , b α ] T and B o x ( · ) stands for the interval hull. For any zonotope, the interval hull is the smallest interval vector.
Lemma 4
([29]). If Ω = ν , H , the components of its interval hull are
a i = ν i j = 0 α H i j , i = 1 , , ι , b i = ν i + j = 0 α H i j , i = 1 , , ι .
Lemma 5
([29]). Given zonotopes Ω i , i = 1 , 2 , , m
B o x ( m = 1 ρ Ω m ) = m = 1 ρ ( B o x ( Ω m ) ) .
Lemma 6
([31]). For the given symmetric matrix A B B T C , the following inequalities are equivalent:
( 1 ) A B B T C 0 ,
( 2 ) C 0 ; A B A 1 B T 0 ,
( 3 ) A 0 ; C B T C 1 B 0 .
Assumption 1.
The initial state of the ith subsystem satisfies the following condition
x ̲ i ( 0 ) x i ( 0 ) x ¯ i ( 0 ) .
Assumption 2.
The disturbances and output noise in system (6) are bounded, which are:
d ( k ) 2 d ,
where d is a constant.
Assumption 3.
The initial state, disturbance and noise, and initial error are represented by x ( 0 ) , d ( 0 ) and e ( 0 ) , which can be wrapped as in the equations below:
x ( 0 ) ν 0 , H 0 = X ( 0 ) , d ( 0 ) 0 , E d = D ( 0 ) , e ( 0 ) 0 , H 0 = E ( 0 ) ,
where H 0 is given matrix and E d = d i a g { d i } .

3. Main Results

In this part, we provide sufficient conditions for the observer with L property for CPS. Then, a reachability analysis methodology is used to design the distributed interval observer.
Theorem 1.
Let η and θ be two given constants with 0 < η < 1 and θ > 1 . If there exists a constant χ and matrices P n 0 R N × N , P m 0 R N × N such that
m i n χ 2 , s u b j e c t t o : η P n 0 χ 2 I N P n A ¯ W n C ¯ κ n Q n ( L n I N ) P n D ¯ P n 0 , P n θ P m , 1 τ ˜ + l n η l n θ < 0 , κ n > 1 a ( L n ) ,
where W n = P n L ¯ n , Q n = P n M ¯ n , n m , n , m S and τ ˜ satisfying ADT. Then (5) is a robust L observer for system (3).
Proof. 
Please see the Appendix A.    □
Remark 3.
In Theorem 1, the parameter κ n depends on a ( L n ) . The algebraic connectivity a ( L n ) of a graph tends to increase with graph stability. As the stability of the graph increases, there will be a greater range of options for κ n .
Remark 4.
In practice, the ADT τ ˜ and the disturbance attenuation level χ 2 stand for the performance of the observer. Owing to the fact that the ADT 1 τ ˜ + l n η l n θ < 0 depends on η and θ, it is necessary to select suitable values for η and θ to minimize the ADT τ ˜ .
After completing the design of the optimal L observer, we need to construct the interval observer by designing the interval hull that can completely wrap the system disturbances, noise, and errors. The real-time error from (5) can be wrapped by the zonotope as
e ( k ) ν 0 , H ( k ) = E ( k ) .
Then, we add the reachability analysis methodology to the distributed L observer designed in Theorem 1 and then give the following interval observer:
x ¯ = x ^ + e ¯ x ̲ = x ^ + e ̲
where e ¯ and e ̲ are the upper and lower bounds of e ( k ) .
Theorem 2.
An interval estimate of the system state is provided by (22), if Assumption 3 holds, then the error e ( k ) has the following upper and lower bounds:
e ̲ ( k ) , e ¯ ( k ) = B o x ( ( n = 0 k 1 Γ n ) E ( 0 ) ) n = 0 k 1 B o x ( ( m = 1 k 1 Γ m ) Γ n 1 D D ) ,
where B o x ( E ( 0 ) ) = e ̲ ( 0 ) , e ¯ ( 0 ) .
Proof. 
Based on Assumption 3, we can obtain E ( 1 ) ,
E ( 1 ) = Γ ϕ ( 0 ) E ( 0 ) D D .
When k = 1 , then we can obtain
E ( 2 ) = Γ ϕ ( 1 ) E ( 1 ) D D = Γ ϕ ( 1 ) Γ ϕ ( 0 ) E ( 0 ) Γ ϕ ( 1 ) D D D D .
Iterating the above process yields
E ( k + 1 ) = ( n = 0 k 1 Γ n ) E ( 0 ) n = 0 k 1 ( m = 1 k 1 Γ m ) Γ n 1 D D .
Then the interval hull below describes the set E ( k )
E ( k ) = e ̲ ( k ) , e ¯ ( k ) = B o x ( ( n = 0 k 1 Γ n ) E ( 0 ) ) n = 0 k 1 B o x ( ( m = 1 k 1 Γ m ) Γ n 1 D D ) .
   □
The proof of Theorem 2 is now completed.
Based on Theorems 1 and 2, the design of a distributed interval observer with switching topology can be implemented by Algorithm 1.
Remark 5.
According to Theorem 1, the robust L observer is designed to reduce the effect of outside disturbance and output noise. It is obvious that a bounded interval hull exists based on the result of Theorem 1, which can include errors and disturbances. Then, using the given interval hull E ( 0 ) as a starting point, Theorem 2 gives the interval hull E ( k ) . We propose a reachable set analysis technique by combining Theorem 1 with Theorem 2.
Algorithm 1: Algorithm for designing the distributed interval observer with switching topology.
(1)
Model CPSs with given disturbances and noise.
(2)
Design distributed observers for subsystems.
(3)
Select the appropriate κ n according to the switching topology.
(4)
Solve the LMI problem in Theorem 1 using the information of the bounds of disturbances and noise.
(5)
Calculate the gains of the observers by M ¯ n = P n 1 Q n , L ¯ n = P n 1 W n .
(6)
Determine the ADT by given η and θ.
(7)
Obtain the zonotopes of the initial value according to (18)
(8)
Transform the zonotope at the moment k = n into an interval hull starting from n = 0 .
(9)
Iterate the interval hull in step (6).
(10)
Obtain the interval hull at the moment k = n + 1 .
(11)
The interval observer is obtained through (21).
Remark 6.
If the modeling uncertainty is taken into consideration, the new model of this paper is as follows:
x i ( k + 1 ) = ( A + Δ A ( k ) ) x i ( k ) + ( B + Δ B ( k ) ) u i ( k ) + p i ( k ) , y i ( k ) = ( C + Δ C ( k ) ) x i ( k ) + q i ( k ) ,
where Δ A ( k ) , Δ B ( k ) and Δ C ( k ) represent the uncertainty of the system, respectively. Then, we may use the norm-bounded condition on Δ A ( k ) , Δ B ( k ) and Δ C ( k ) to design the L interval observer. However, it is not an easy task to construct the interval hull of the corresponding error system since the time-varying terms are contained. In the near future, we will deal with systems with model uncertainty with the method proposed in this paper.
Remark 7.
In Theorems 1 and 2, limited by the current knowledge of authors, this paper only gives sufficient conditions for observer design. In the future, we will study the necessary conditions for the design of interval observers, and in conjunction with this paper, we will give the necessary and sufficient conditions for the design of interval observers.

4. Simulation

Among CPSs, UAVs have recently achieved widespread application. In this section, we simulate through a networked CPSs with four UAVs.
Referring to [11], the dynamical system of each UAV is as follows
α ˙ i ( t ) β ˙ i ( t ) = A α i ( t ) β i ( t ) + B u i ( k ) + p i ( k ) , y i ( k ) = C α i ( t ) β i ( t ) + q i ( k ) ,
where the pitch rate and angle of attack of each UAV are indicated by β i ( t ) and α i ( t ) . The schematic of each UAV is shown in Figure 1, and the detailed derivation of the dynamics model is omitted here.
Table 1 displays the output noise and the external disturbance borrowed from [11], and the partial matrix values are as follows:
A = 0.8825 0.0987 0.8458 0.9122 , B = 0.0194 0.0036 1.9290 0.3803 , C = 1 0.2 .
For illustrative purposes, Figure 2 shows a switching communication topology G 1 , G 2 and G 3 with four UAVs. Figure 3 displays the change in the switching signal ϕ ( k ) . Then, from Figure 2, we can give the corresponding matrix L n :
L 1 = 1 1 0 0 1 3 1 1 0 1 2 1 0 1 1 2 , L 2 = 1 1 0 0 0 2 1 1 0 0 1 1 1 1 0 2 , L 3 = 1 0 0 1 0 1 0 1 0 0 1 1 1 1 1 3 .
From Lemmas 1 and 2, we can obtain a ( L 1 ) = a ( L 2 ) = a ( L 3 ) = 1 , then we have κ > 1 . In this simulation, κ 1 = 1.8 , κ 2 = 2.3 , κ 3 = 3.2 are chosen. We can determine L n and M n by solving (19), the observer gains for 1th subsystem are listed below
L 1 = 47.8095 46.2386 , L 2 = 9.5528 47.7273 , L 3 = 17.8244 11.9289 .
M 1 = 41.7480 36.5506 1.8874 6.9585 , M 2 = 10.3302 8.4984 79.5140 71.1991 , M 3 = 11.7589 7.5478 6.1219 9.0834 .
The peak estimation error of Angle of Attack and Pitch Rate of 1th subsystem are 0.4184 and 0.7668 , respectively. The disturbance attenuation level χ 2 and ADT τ ˜ are as follows
χ 2 = 0.4984 , τ ˜ 0.1763 .
Below are the results of the numerical simulation. The states of the original systems and the upper and lower observers are depicted in Figure 4 and Figure 5. x i j are the original states of the subsystems. u i j and v i j reflect the bounds produced via the interval hull technique used in this paper and the monotone system method in [32], where i denotes the ith subsystem and j denotes the jth state. Figure 6 shows the error system for a single subsystem, e u i j , and e v i j reflect the observation error through the method used in this paper and the monotone system method in [32].
It is evident that the states of CPS are completely surrounded by those of the upper and lower observers. From Figure 4, Figure 5 and Figure 6, it can be seen that the distributed interval observer designed in this paper has higher estimation accuracy compared to the traditional monotone system approach in [32]. We design the optimal robust observer using the L technique, which reduces the design requirements of the observer, unlike the H technique applied in [24]. Different from [23] and most of the work, we consider for the first time the case of switching topology among distributed interval observers. Thus, the proposed distributed interval observers design method for CPSs is effective and feasible.
Remark 8.
For networked CPSs with perturbations and noise, we propose a class of distributed interval observer design methods with switching topology that combine the design of L observers with interval hulls. This class of methods eliminates the need to consider the error system to be Schur and, therefore, eliminates the need to use coordinate transformation methods, significantly reducing the conservatism of the estimation. It can be seen from the observer form (3), sufficient conditions (19), and the proof of Theorem 1. In addition, the estimation accuracy is better than that in [32], as it can also be seen from (22).

5. Conclusions

In this paper, a distributed interval estimation method for uncertain CPSs is investigated. Due to the communication burden of networked CPSs, we consider the case of switching topology among distributed interval observers. To improve the accuracy of the estimation, a reachability analysis is introduced in conjunction with the L technique. Finally, the validity of the main results of this paper is verified by one example. In the future, we may focus our research on the distributed interval estimation of the attacked CPSs.

Author Contributions

Methodology, J.Z.; validation, C.L.; writing—original draft preparation, J.Z.; writing—review and editing, J.H.; supervision, J.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Proof. 
Define A = η P n 0 0 χ 2 I N , B = P n A ¯ W n C ¯ κ n Q n ( L n I N ) P n D ¯ , C = P n . Using Lemma 6 and the fact that C 0 , we can determine that
A B A 1 B T 0 .
Substituting Q n = P n M ¯ n , W n = P n L ¯ n into (5), we obtain
G = Γ n T P n Γ n η P n D ¯ T P n Γ n D ¯ T P n D ¯ χ 2 I N 0 .
where Γ n = A ¯ L ¯ n C ¯ κ n M ¯ n ( L n I N ) . It follows from (A2) that
e T ( k ) d T ( k ) Γ n T P n Γ n η P n D ¯ T P n Γ n D ¯ T P n D ¯ χ 2 I N e ( k ) d ( k ) 0 .
Thus, (A3) implies that
e T ( k ) ( Γ n T P n Γ n η P n ) e ( k ) + d T ( k ) ( D ¯ T P n D ¯ χ 2 I N ) d ( k ) + 2 d T ( k ) D ¯ T P n Γ n e ( k ) < 0 .
We then choose the following Lyapunov function
V n ( k ) = e T ( k ) P n e ( k ) .
Thus,
V n ( k ) = V n ( k + 1 ) V n ( k ) = e T ( k ) ( Γ n T P n Γ n η P n ) e ( k ) + d T ( k ) ( D ¯ T P n D ¯ χ 2 I N ) d ( k ) + 2 d T ( k ) D ¯ T P n Γ n e ( k ) .
Further simplification of (A6) yields
V n ( k ) < ( η 1 ) e T ( k ) P n e ( k ) + χ 2 d T ( k ) d ( k ) ,
which means
V n ( k + 1 ) < η V n ( k ) + χ 2 d T ( k ) d ( k ) .
Then consider the interval [ k σ , k ) . Iterating (A8) yields
V n ( k ) < η k k σ V n ( k σ ) + χ 2 ϱ = 0 k k σ 1 η ϱ d T ( k 1 ϱ ) d ( k 1 ϱ ) .
Suppose that ϕ ( k σ 1 ) = m , using P n θ P m , then
V n ( k ) < η k k σ V n ( k σ ) + χ 2 ϱ = 0 k k σ 1 η ϱ d T ( k 1 ϱ ) d ( k 1 ϱ ) < η k k σ + 1 θ V m ( k σ 1 ) + η k k σ θ χ 2 d T ( k σ 1 ) d ( k σ 1 ) + χ 2 ϱ = 0 k k σ 1 η l d T ( k 1 ϱ ) d ( k 1 ϱ ) .
Iterating (A10) yields the inequality below
V n ( k ) < η k θ N ϕ ( 0 , k ) V ϕ ( 0 ) ( 0 ) + χ 2 ϱ = 0 k k σ 1 η ϱ θ N ϕ ( k ϱ , k ) d T ( k 1 ϱ ) d ( k 1 ϱ ) .
In view of θ > 1 and N ϕ ( 0 , k ) N ϕ ( k ϱ , k ) , we have
χ 2 ϱ = 0 k k σ 1 η ϱ θ N ϕ ( k ϱ , k ) d T ( k 1 ϱ ) d ( k 1 ϱ ) < χ 2 ϱ = 0 k k σ 1 η ϱ θ N ϕ ( 0 , k ) d T ( k 1 ϱ ) d ( k 1 ϱ ) .
It follows from (A11) and (A12) that
V n ( k ) < θ N ϕ ( 0 , k ) [ η k V ϕ ( 0 ) ( 0 ) + χ 2 ϱ = 0 k k σ 1 η l d T ( k 1 ϱ ) d ( k 1 ϱ ) ] .
Due to V n ( k ) λ ̲ e ( k ) 2 2 , the error e ( k ) in (6) satisfies
e ( k ) 2 2 1 λ ̲ ( P n ) V n ( k ) , e ( k ) 2 2 θ N ϕ ( 0 , k ) λ ̲ ( P n ) ( η k V ϕ ( 0 ) ( 0 ) + χ 2 ϱ = 0 k 1 η ϱ d ( k ) 2 ) .
By Definition 1, the proof of Theorem 1 is completed. □

References

  1. Cao, X.; Cheng, P.; Chen, J.; Sun, Y. An online optimization approach for control and communication codesign in networked cyber-physical systems. IEEE Trans. Ind. Inform. 2012, 9, 439–450. [Google Scholar] [CrossRef]
  2. Gunes, V.; Peter, S.; Givargis, T.; Vahid, F. A survey on concepts, applications, and challenges in cyber-physical systems. KSII Trans. Internet Inf. Syst. 2014, 8, 4242–4268. [Google Scholar]
  3. Kim, K.D.; Kumar, P.R. Cyber–physical systems: A perspective at the centennial. Proc. IEEE 2012, 100, 1287–1308. [Google Scholar]
  4. Cai, S.; Lau, V.K. Zero MAC latency sensor networking for cyber-physical systems. IEEE Trans. Signal Process. 2018, 66, 3814–3823. [Google Scholar] [CrossRef]
  5. Sangaiah, A.K.; Siarry, P. Cognitive Brain-Inspired Cyber-Physical Systems in Industrial Informatics. Front. Neurorobot. 2022, 16, 926538. [Google Scholar] [CrossRef]
  6. Mitchell, R.; Chen, R. Modeling and analysis of attacks and counter defense mechanisms for cyber physical systems. IEEE Trans. Reliab. 2015, 65, 350–358. [Google Scholar] [CrossRef]
  7. Verma, R. Smart city healthcare cyber physical system: Characteristics, technologies and challenges. Wirel. Pers. Commun. 2022, 122, 1413–1433. [Google Scholar] [CrossRef]
  8. Hozdić, E.; Butala, P. Concept of socio-cyber-physical work systems for Industry 4.0. Tehnički Vjesn. 2020, 27, 399–410. [Google Scholar]
  9. Wu, C.; Hu, Z.; Liu, J.; Wu, L. Secure estimation for cyber-physical systems via sliding mode. IEEE Trans. Cybern. 2018, 48, 3420–3431. [Google Scholar] [CrossRef]
  10. Chang, Y.H.; Hu, Q.; Tomlin, C.J. Secure estimation based Kalman filter for cyber–physical systems against sensor attacks. Automatica 2018, 95, 399–412. [Google Scholar] [CrossRef]
  11. Zhang, C.L.; Yang, G.H.; Lu, A.Y. Resilient observer-based control for cyber-physical systems under denial-of-service attacks. Inf. Sci. 2021, 545, 102–117. [Google Scholar] [CrossRef]
  12. Gouzé, J.L.; Rapaport, A.; Hadj-Sadok, M.Z. Interval observers for uncertain biological systems. Ecol. Model. 2000, 133, 45–56. [Google Scholar] [CrossRef]
  13. 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]
  14. Mazenc, F.; Dinh, T.N.; Niculescu, S.I. Interval observers for discrete-time systems. Int. J. Robust Nonlinear Control 2014, 24, 2867–2890. [Google Scholar] [CrossRef]
  15. Tang, W.; Wang, Z.; Wang, Y.; Raïssi, T.; Shen, Y. Interval estimation methods for discrete-time linear time-invariant systems. IEEE Trans. Autom. Control 2019, 64, 4717–4724. [Google Scholar] [CrossRef]
  16. Zhou, X.; Wang, Z.; Wang, J. Automated Ground Vehicle Path-Following: A Robust Energy-to-Peak Control Approach. IEEE Trans. Intell. Transp. Syst. 2022, 23, 14294–14305. [Google Scholar] [CrossRef]
  17. Zhang, W.; Chen, Y.; Gao, H. Energy-to-peak control for seismic-excited buildings with actuator faults and parameter uncertainties. J. Sound Vib. 2011, 330, 581–602. [Google Scholar] [CrossRef]
  18. Feng, J.; Han, K. Robust full- and reduced-order energy-to-peak filtering for discrete-time uncertain linear systems. Signal Process. 2015, 108, 183–194. [Google Scholar] [CrossRef]
  19. 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]
  20. Guo, S.; Ren, W.; Ahn, C.K.; Wen, C.; Lam, H.K. Reachability analysis-based interval estimation for discrete-time Takagi–Sugeno fuzzy systems. IEEE Trans. Fuzzy Syst. 2021, 30, 1981–1992. [Google Scholar] [CrossRef]
  21. Ren, W.; Guo, S. State and Faults Interval Estimations for Discrete-time Linear Systems. Int. J. Control Autom. Syst. 2023, 21, 2303–2312. [Google Scholar] [CrossRef]
  22. Zhang, Z.H.; Yang, G.H. Distributed fault detection and isolation for multiagent systems: An interval observer approach. IEEE Trans. Syst. Man Cybern. Syst. 2018, 50, 2220–2230. [Google Scholar] [CrossRef]
  23. Li, D.; Chang, J.; Chen, W.; Raïssi, T. IPR-based distributed interval observers design for uncertain LTI systems. ISA Trans. 2022, 121, 147–155. [Google Scholar] [CrossRef] [PubMed]
  24. Huang, J.; Zhang, H.; Raïssi, T. Distributed interval estimation methods for multiagent systems. IEEE Syst. J. 2022, 17, 1843–1852. [Google Scholar] [CrossRef]
  25. 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. Part B (Cybern.) 2009, 40, 881–891. [Google Scholar]
  26. Li, Z.; Liu, X.; Fu, M.; Xie, L. Global H consensus of multi-agent systems with Lipschitz non-linear dynamics. IET Control Theory Appl. 2012, 6, 2041–2048. [Google Scholar] [CrossRef]
  27. Briat, C. Robust stability and stabilization of uncertain linear positive systems via integral linear constraints: L1-gain and L-gain characterization. Int. J. Robust Nonlinear Control 2013, 23, 1932–1954. [Google Scholar] [CrossRef]
  28. Zhao, X.; Zhang, L.; Shi, P.; Liu, M. Stability of switched positive linear systems with average dwell time switching. Automatica 2012, 48, 1132–1137. [Google Scholar] [CrossRef]
  29. Alamo, T.; Bravo, J.M.; Camacho, E.F. Guaranteed state estimation by zonotopes. Automatica 2005, 41, 1035–1043. [Google Scholar] [CrossRef]
  30. Combastel, C. Zonotopes and Kalman observers: Gain optimality under distinct uncertainty paradigms and robust convergence. Automatica 2015, 55, 265–273. [Google Scholar] [CrossRef]
  31. Boyd, S.; El Ghaoui, L.; Feron, E.; Balakrishnan, V. Linear Matrix Inequalities in System and Control Theory; SIAM: Philadelphia, PA, USA, 1994. [Google Scholar]
  32. Zhang, H.; Huang, J.; He, S. Fractional-order interval observer for multiagent nonlinear systems. Fractal Fract. 2022, 6, 355. [Google Scholar] [CrossRef]
Figure 1. Longitudinal axis system of UAV.
Figure 1. Longitudinal axis system of UAV.
Mathematics 12 00163 g001
Figure 2. Three switching communication topology of UAVs.
Figure 2. Three switching communication topology of UAVs.
Mathematics 12 00163 g002
Figure 3. The change in the switching signal of distributed interval observers.
Figure 3. The change in the switching signal of distributed interval observers.
Mathematics 12 00163 g003
Figure 4. Angle of attack and interval estimates of UAVs.
Figure 4. Angle of attack and interval estimates of UAVs.
Mathematics 12 00163 g004
Figure 5. Pitch rate and interval estimates of UAVs.
Figure 5. Pitch rate and interval estimates of UAVs.
Mathematics 12 00163 g005
Figure 6. Observation error of angle of attack and pitch rate of UAVs.
Figure 6. Observation error of angle of attack and pitch rate of UAVs.
Mathematics 12 00163 g006
Table 1. The output noise and the external disturbance.
Table 1. The output noise and the external disturbance.
SubsystemOutput NoiseExternal Disturbance
1 0.1 cos ( k ) 0.2 + 0.1 cos ( 0.5 k )
2 0.1 sin ( k ) 0.1 cos ( 0.1 k )
3 0.01 cos ( k ) cos ( 0.2 π k ) + 0.3 sin ( 0.2 π k )
4 0.05 sin ( k ) cos ( 0.3 π k ) + 0.1 sin ( 0.3 π k )
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, J.; Huang, J.; Li, C. Distributed Interval Observers with Switching Topology Design for Cyber-Physical Systems. Mathematics 2024, 12, 163. https://doi.org/10.3390/math12010163

AMA Style

Zhang J, Huang J, Li C. Distributed Interval Observers with Switching Topology Design for Cyber-Physical Systems. Mathematics. 2024; 12(1):163. https://doi.org/10.3390/math12010163

Chicago/Turabian Style

Zhang, Junchao, Jun Huang, and Changjie Li. 2024. "Distributed Interval Observers with Switching Topology Design for Cyber-Physical Systems" Mathematics 12, no. 1: 163. https://doi.org/10.3390/math12010163

APA Style

Zhang, J., Huang, J., & Li, C. (2024). Distributed Interval Observers with Switching Topology Design for Cyber-Physical Systems. Mathematics, 12(1), 163. https://doi.org/10.3390/math12010163

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