Design of Distributed Interval Observers for Multiple Euler–Lagrange Systems

Abstract

This paper investigates the problem of distributed interval estimation for multiple Euler–Lagrange systems. An interconnection topology is supposed to be strongly connected. To design distributed interval observers, the coordinate transformation method is employed. The construction of the distributed interval observer is given by the monotone system theory, and the stability is analyzed by the Lyapunov stability theory. Unlike the current works, each sub-interval observer has its own gain; in addition to this, additional observer gains are used to reduce the conservatism of design. The gains of all sub-interval observers are determined by both the monotone system theory and the Lyapunov stability theory. Finally, a simulation example verifies the feasibility of the presented distributed interval observers.

1. Introduction

Euler–Lagrange systems (ELSs) have been studied by a wide range of scholars as a common class of systems for describing many real mechanical models. In fact, some of the variables in a specific mechanical system are not directly available. Therefore, a number of scholars have addressed the observer design problem for ELSs. A speed observer design method for general ELSs was presented in [1]. The authors in [2] designed a nonlinear disturbance observer for ELSs and realized the estimation of unknown disturbances. Ref. [3] used a sliding-mode disturbance observer for ELSs to design a semi-global exponential controller. With the deepening of research and the development of multi-agent systems, scholars investigated the problem of the state estimation of multiple ELSs (MELSs) in recent years. In [4], the authors designed distributed observers (DOs) for a class of MELSs without external disturbances. Ref. [5] studied the output-based tracking control problem for MELSs, in which the authors designed a novel nonlinear DO. In [6], a DO-based fault-detection method for MELSs was introduced, and a DO-based controller was also designed. In addition to this, the investigation on MELSs with uncertainties had received considerable attention. Ref. [7] addressed the problem of leader-following consensus of uncertain MELSs under switching network topology and proposed a DO design method to recover the states. Moreover, the distributed adaptive observers were designed for MELSs with uncertainties in [8]. There are also other interesting works on MELSs [9,10,11,12] herein.

As an effective estimation method, interval observers (IOs) can recover the upper and lower boundaries of uncertain systems. The concept of IOs was introduced in [13]. The IOs are further divided into centralized IOs (CIOs) and distributed IOs (DIOs). There are many studies on CIOs [14,15,16,17,18,19,20,21]. For continuous systems, an IO design method was introduced for linear time-invariant systems with disturbances in [14]. The authors in [15] gave an effective design method of IO for time-varying systems. For discrete-time systems, [16] formulated a design framework of IO, and [17] investigated the IO design problem for time-varying discrete-time systems. In [18], the IO design method for discrete-time switched systems was studied. Over the past decade, the research of functional IO has also attracted the attention of some researchers [19,20,21]. On the other hand, the investigation on DIOs has also been reported by some works recently such as [22,23,24]. The authors of [22] presented an method to design DIOs, which was applied to fault detection for multi-agent systems. The authors of [23] considered the design method of DIOs for a class of linear time-invariant systems with uncertainties. Finally, The authors of [24] extended the DIOs design method to a class of fractional-order systems. The problem of estimation of MELSs is always complicated due to the uncertainty contained in the system itself. Most of the existing works add constraints to achieve state estimation; however, this problem can be avoided if an interval estimation approach is used and state estimation can also be achieved. Therefore, it is necessary to design DIOs for MELSs. To the authors’ best knowledge, there are few reports on DIOs for MELSs.

Based on the above discussion, this paper studies the DIOs design problem of MELSs. First of all, the coordinate transformation method must necessarily be introduced. Then, the framework of DIOs is constructed based on the monotone system theory. Additionally, the stability of DIO is guaranteed by Lyapunov stability theory. The rest of the paper is structured as follows. Section 2 gives the preliminaries, which include the graph theory and system model. Section 3 provides the DIO design method based on monotone system theory. In Section 4, a numerical example is used to show the validity of the presented design method of DIO. Finally, Section 5 is the conclusion of the paper.

Notaion: For a matrix E, $E^{+}$ denotes $max\{0,E\}$, $E^{-}$ = $E^{+}-E$. For a real symmetric matrix $O\in R^{n\times n},O\succ 0(O\prec 0)$ indicates that O is positive (negative) definite, $He\left(O\right)=O^T+O$. ⊗ represents the Kronecker product, and the diag{·} denotes a block diagonal matrix.

2. Preliminaries

2.1. Graph Theory

For a digraph G with N vertices, the adjacency matrix $\mathcal{A}\in R^{n\times n}=\left[a_{ij}\right]$ is given by $a_{ij}=0$, if $(i,j)\in$S, where S is the set of edges, and $a_{ij}=1$, otherwise. The length path from a vertex i to a vertex j is a sequence of $r+1$ distinct vertices starting at i and ending at j, and consecutive vertices are adjacent to each other. A graph G is said to be connected if there is a path between any two vertices of the graph G. The Laplacian matrix $\mathcal{L}$ of graph $\mathcal{G}$ is defined as $\mathcal{L}:=\mathcal{D}-\mathcal{A}$. $\mathcal{D}$ is called the degree matrix of $\mathcal{G}$. For a connected graph, the Laplacian matrix has a single zero eigenvalue and the corresponding eigenvector is $1_N$. Additionally, if G is connected then $0={\lambda }_1\left(G\right)\le {\lambda }_2\left(G\right)\le \cdots \le {\lambda }_N\left(G\right)$, where ${\lambda }_i\left(G\right)$ is the eigenvalue of $\mathcal{L}$.

Lemma 1 

([25]). Assume that $\mathcal{G}$ is a strongly connected graph. Denote $r_i(i=1,\cdots ,N)$ as the left eigenvector with eigenvalue 0, and $R=diag\{r_1,\cdots \cdots ,r_N\}$; then, $R\mathcal{L}+{\mathcal{L}}^{T}R\ge 0$ holds.

Lemma 2 

([26]). If graph $\mathcal{G}$ is strongly connected, the generalized algebraic connectivity of $\mathcal{G}$ could be represented as $a\left(\mathcal{L}\right)=\underset{r^{T}x=0,x\ne 0}{min}\frac{x^T(R\mathcal{L}+{\mathcal{L}}^{T}R)x}{2x^{T}Rx}$. If the topology is balanced, the matrix $R=r_1{I}_N$. Then, $a\left(\mathcal{L}\right)={\lambda }_{min}\left(\frac{He\left(\mathcal{L}\right)}{2}\right)$ holds.

2.2. System Model

In general, an ELS with second-order dynamics can be written as

$$M_i\left(q_i\right){\ddot{q}}_i+C_i(q_i,{\dot{q}}_i){\dot{q}}_i+G_i\left(q_i\right)=u_i,$$

(1)

where $q_i\in R^n$ is the vector of generalized coordinates, $M_i\left(q_i\right)\in R^{n\times n}$ represents the inertia matrix, $C(q_i,{\dot{q}}_i)\in R^{n\times n}$ denotes the vector of Coriolis and centripetal forces, $G\left(q\right)\in R^n$ is the gravity vector, and $u_i\in R^n$ stands for the vector of generalized control input forces. Next, we take the following transformations for system (1):

$$\left\{\begin{array}{c}\begin{array}{cc}\hfill & {\dot{x}}_1=x_2,\hfill \\ \hfill & {\dot{x}}_2=h\left(x_1\right)u+f(x_1,x_2),\hfill \end{array}\end{array}\right.$$

(2)

where $x_{i1}=q_i$, $x_{i2}={\dot{q}}_i$, $h\left(x_{i1}\right)=M{\left(x_{i1}\right)}^{-1}$, $f(x_{i1},x_{i2})=-M{\left(x_{i1}\right)}^{-1}(C(x_{i1},x_{i2})x_{i2}+G\left(x_{i1}\right))$. The joint displacements are measurable; under the above state reconstruction, we can obtain the following equation:

$$\left\{\begin{array}{c}\begin{array}{cc}\hfill & {\dot{x}}_i\left(t\right)=A{x}_i\left(t\right)+B{u}_i\left(t\right)+{\varphi }_i\left(x_i\left(t\right)\right),\hfill \\ \hfill & y_i\left(t\right)=C{x}_i\left(t\right),\hfill \end{array}\end{array}\right.$$

(3)

where $A=\left[\begin{array}{cc}{0}_{n\times n}& I_{n\times n}\\ 0_{n\times n}& 0_{n\times n}\end{array}\right]\in R^{2n\times 2n}$, $B=\left[\begin{array}{c}{0}_{n\times n}\\ h\left(x_{i1}\right)\end{array}\right]\in R^{2n\times n}$, $C=\left[\begin{array}{c}{I}_{n\times n}{0}_{n\times n}\end{array}\right]\in R^{n\times 2n}$, and ${\varphi }_i\left(x_i\left(t\right)\right)=\left[\begin{array}{c}{0}_{n\times 1}\\ f(x_{i1},x_{i2})\end{array}\right]\in R^{2n}$.

Lemma 3 

([27]). If $\varphi \left(x_i\right)$ is a Lipschitz function that is globally differentiable, then there are two increasing Lipschitz functions $f\left(x_i\right)$ and $g\left(x_i\right)$ such that

$$\varphi \left(x_i\right)=f\left(x_i\right)-g\left(x_i\right).$$

(4)

Lemma 4 

([27]). For the function $\varphi \left(x_i\right)$ in Lemma 3, a global Lipschitz function $\overline{\varphi }(x_{ia},x_{ib})$ exists such that:

⋆

${\displaystyle \overline{\varphi }(x_i,x_i)=\varphi \left(x_i\right)}$,

⋆

${\displaystyle \frac{\partial \overline{\varphi }}{x_{ia}}\ge 0}$, and ${\displaystyle \frac{\partial \overline{\varphi }}{x_{ib}}\le 0}$.

The above lemmas can help us obtain the boundaries of nonlinear function $\varphi (x_i,x_i)$:

$${\underline{x}}_i\le x\le {\overline{x}}_i\Rightarrow \overline{\varphi }({\underline{x}}_i,{\overline{x}}_i)\le \varphi (x_i,x_i)\le \overline{\varphi }({\overline{x}}_i,{\underline{x}}_i).$$

(5)

Lemma 5 

([24]). For $\varphi \left(x_i\right)$, $\overline{\varphi }({\overline{x}}_i,{\underline{x}}_i)$ and $\overline{\varphi }({\underline{x}}_i,{\overline{x}}_i)$ defined in Lemma 4, matrices $F_i\in \{1,2,3,4\}$ exist such that

$$\left\{\begin{array}{c}\begin{array}{cc}\hfill & \overline{\varphi }({\overline{x}}_i,{\underline{x}}_i)-\varphi \left(x_i\right))\le F_1{\overline{e}}_i+F_2{\underline{e}}_i,\hfill \\ \hfill & \varphi \left(x_i\right)-\overline{\varphi }({\underline{x}}_i,{\overline{x}}_i)\le F_3{\overline{e}}_i+F_4{\underline{e}}_i,\hfill \end{array}\end{array}\right.$$

(6)

where ${\underline{e}}_i=x_i-\underline{x_i}$ and ${\overline{e}}_i={\overline{x}}_i-x_i$.

Lemma 6 

([14]). Given a constant matrix $A\in R^{m\times n}$ and a vector $\zeta \in R^{n\times 1}$, if $\zeta \in [\underline{\zeta },\overline{\zeta }]$ holds, then

$$A^{+}\underline{\zeta }-A^{-}\overline{\zeta }\le A\zeta \le A^{+}\overline{\zeta }-A^{-}\underline{\zeta }.$$

(7)

In order to obtain the main results, we need the following assumptions.

Assumption 1. 

The nonlinear function ${\varphi }_i\left(x_i\left(t\right)\right)$ is considered as a global Lipschitz function.

Assumption 2. 

For every sub-system, the initial state of system (3) satisfies the following inequality:

$${\underline{x}}_i\left(0\right)\le x_i\left(0\right)\le {\overline{x}}_i\left(0\right).$$

(8)

Assumption 3. 

The topology graph is considered to be balanced and strongly connected.

Assumption 4. 

Matrices L, M, $\overline{Y}$, and $\underline{Y}$ exist such that Π is Metzler with

$$\begin{array}{c}\hfill \mathrm{\Pi }=\left[\begin{array}{cc}\mathrm{\Psi }+\overline{Y}& \overline{Y}\\ \underline{Y}& \mathrm{\Psi }+\underline{Y}\end{array}\right],\end{array}$$

where $\mathrm{\Psi }=\hat{A}-L\hat{C}-\gamma M(\mathcal{L}\otimes I_n)$.

Remark 1. 

The above assumptions have been commonly used and widely accepted in previous studies. Similar to [22], the topology graph is supposed to be balanced and strongly connected in this paper; based on this, we provide the following results.

3. Main Results

To design DIO for system (3), some preliminary preparations are needed. First, the equivalence transformation is performed for system matrix $A\in R^{n\times n}$. Using the linear transformation $z_i=H{x}_i$, the system (3) becomes

$$\left\{\begin{array}{c}\begin{array}{cc}\hfill & \dot{z_i}\left(t\right)=\overline{A}{z}_i\left(t\right)+\overline{B}{u}_i\left(t\right)+H{\varphi }_i\left(H^{-1}{z}_i\left(t\right)\right),\hfill \\ \hfill & y_i\left(t\right)=\overline{C}{z}_i\left(t\right),\hfill \end{array}\end{array}\right.$$

(9)

where $\overline{A}=HA{H}^{-1}$, $\overline{B}=HB$, and $\overline{C}=C{H}^{-1}$.

Remark 2. 

In order to estimate the bounds of the nonlinear function ${\varphi }_i\left(x\right)$, we define that

$$\left\{\begin{array}{c}\begin{array}{cc}\hfill & {\phi }_i({\overline{z}}_i,{\underline{z}}_i)={\varphi }_i(H_a^{+}{\overline{z}}_i-H_a^{-}{\underline{z}}_i,H_a^{+}{\underline{z}}_i-H_a^{-}{\overline{z}}_i),\hfill \\ \hfill & {\phi }_i({\underline{z}}_i,{\overline{z}}_i)={\varphi }_i(H_a^{+}{\underline{z}}_i-H_a^{-}{\overline{z}}_i,H_a^{+}{\overline{z}}_i-H_a^{-}{\underline{z}}_i),\hfill \end{array}\end{array}\right.$$

(10)

where ${\overline{z}}_i\left(t\right)$ and ${\underline{z}}_i\left(t\right)$ represent the estimated values of $z_i\left(t\right)$, and $H_a=H^{-1}$.

Under Lemma 6, the following inequality can be obtained:

$$\begin{array}{cc}\hfill {\mathrm{\Phi }}_i({\underline{z}}_i,{\overline{z}}_i)& =H^{+}{\phi }_i({\underline{z}}_i,{\overline{z}}_i)-H^{-}{\phi }_i({\overline{z}}_i,{\underline{z}}_i)\hfill \\ \hfill & \le H{\varphi }_i\left(H^{-1}{z}_i\right)\hfill \\ \hfill & \le H^{+}{\phi }_i({\overline{z}}_i,{\underline{z}}_i)-H^{-}{\phi }_i({\underline{z}}_i,{\overline{z}}_i)\hfill \\ \hfill & ={\mathrm{\Phi }}_i({\overline{z}}_i,{\underline{z}}_i).\hfill \end{array}$$

Remark 3. 

After the transformation of the coordinates, the conclusion of Lemma 5 takes the following form:

$$\left\{\begin{array}{c}\begin{array}{cc}\hfill & {\mathrm{\Phi }}_i({\overline{z}}_i,{\underline{z}}_i)-H{\varphi }_i\left(H^{-1}{z}_i\right)\le N_1{\overline{e}}_i+N_2{\underline{e}}_i,\hfill \\ \hfill & H{\varphi }_i\left(H^{-1}{z}_i\right)-{\mathrm{\Phi }}_i({\underline{z}}_i,{\overline{z}}_i)\le N_3{\overline{e}}_i+N_4{\underline{e}}_i,\hfill \end{array}\end{array}\right.$$

(11)

where $N_i=F_{i}H(i=1,2,3,4)$ are constant matrices.

To achieve $\underline{z}\le z\le \overline{z}$, the following DIOs for each sub-system are constructed:

$$\left\{\begin{array}{c}\begin{array}{cc}\hfill {\dot{\overline{z}}}_i\left(t\right)=& \overline{A}{\overline{z}}_i\left(t\right)+\overline{B}{u}_i\left(t\right)+{\mathrm{\Phi }}_i({\overline{z}}_i,{\underline{z}}_i)+L_i(y_i-\overline{C}{\overline{z}}_i\left(t\right))+{\overline{Y}}_i({\overline{z}}_i\left(t\right)-{\underline{z}}_i\left(t\right))\hfill \\ \hfill & +\gamma M_i\sum _{j=0}^N{a}_{ij}({\overline{z}}_j-{\overline{z}}_i),\hfill \\ \hfill {\dot{\underline{z}}}_i\left(t\right)=& \overline{A}{\underline{z}}_i\left(t\right)+\overline{B}{u}_i\left(t\right)+{\mathrm{\Phi }}_i({\underline{z}}_i,{\overline{z}}_i)+L_i(y_i-\overline{C}{\underline{z}}_i\left(t\right))-{\underline{Y}}_i({\overline{z}}_i\left(t\right)-{\underline{z}}_i\left(t\right))\hfill \\ \hfill & +\gamma M_i\sum _{j=0}^N{a}_{ij}({\underline{z}}_j-{\underline{z}}_i).\hfill \end{array}\end{array}\right.$$

(12)

Then, the error system for i-th sub-system is obtained

$$\left\{\begin{array}{c}\begin{array}{cc}\hfill {\dot{\overline{e}}}_i\left(t\right)=& {\dot{\overline{z}}}_i-{\dot{z}}_i\hfill \\ \hfill =& (\overline{A}-L_i\overline{C}-\gamma M_i\sum _{j=0}^N{\mathcal{L}}_{ij}){\overline{e}}_i\left(t\right)+{\mathrm{\Phi }}_i({\overline{z}}_i,{\underline{z}}_i)-H{\varphi }_i\left(H^{-1}{z}_i\left(t\right)\right)+{\overline{Y}}_i({\overline{e}}_i\left(t\right)-{\underline{e}}_i\left(t\right)),\hfill \end{array}\\ \begin{array}{cc}\hfill {\dot{\underline{e}}}_i\left(t\right)=& {\dot{\underline{z}}}_i-{\dot{z}}_i\hfill \\ \hfill =& (\overline{A}-L_i\overline{C}-\gamma M_i\sum _{j=0}^N{\mathcal{L}}_{ij}){\underline{e}}_i\left(t\right)+H{\varphi }_i\left(H^{-1}{z}_i\left(t\right)\right)-{\mathrm{\Phi }}_i({\overline{z}}_i,{\underline{z}}_i)+{\underline{Y}}_i({\overline{e}}_i\left(t\right)-{\underline{e}}_i\left(t\right)).\end{array}\end{array}\right.$$

(13)

Under Assumption 3, the dynamics of the global system are expressed as

$$\left\{\begin{array}{c}\begin{array}{cc}\hfill & \dot{z}\left(t\right)=\hat{A}z\left(t\right)+\hat{B}u\left(t\right)+\hat{\varphi }\left(z\left(t\right)\right),\hfill \\ \hfill & y\left(t\right)=\hat{C}z\left(t\right),\hfill \end{array}\end{array}\right.$$

(14)

where $z={[z_1^T,\cdots ,z_N^T]}^T$, $y={[y_1^T,\cdots ,y_N^T]}^T$, $\hat{A}=diag\{\underset{N}{\underbrace{\overline{A},\cdots ,\overline{A}}}\}$, $\hat{B}=diag\{\underset{N}{\underbrace{\overline{B},\cdots ,\overline{B}}}\}$, $u={[u_1^T,\cdots ,u_N^T]}^T$, $\hat{\varphi }={[{\left(H{\varphi }_1\right)}^T,\cdots ,{\left(H{\varphi }_N\right)}^T]}^T$, $\hat{C}=diag\{\underset{N}{\underbrace{\overline{C},\cdots ,\overline{C}}}\}$.

The dynamics of the global observer system are given by

$$\left\{\begin{array}{c}\begin{array}{cc}\hfill \dot{\overline{z}}\left(t\right)=& \hat{A}\overline{z}\left(t\right)+\hat{B}u\left(t\right)+\overline{\mathrm{\Phi }}+L(y-\hat{C}\overline{z}\left(t\right))+\overline{Y}(\overline{z}\left(t\right)-\underline{z}\left(t\right))\hfill \\ \hfill & -\gamma M(\mathcal{L}\otimes I_n)\overline{z}\left(t\right),\hfill \\ \hfill \dot{\underline{z}}\left(t\right)=& \hat{A}\underline{z}\left(t\right)+\hat{B}u\left(t\right)+\underline{\mathrm{\Phi }}+L(y-\hat{C}\underline{z}\left(t\right))-\underline{Y}(\overline{z}\left(t\right)-\underline{z}\left(t\right))\hfill \\ \hfill & -\gamma M(\mathcal{L}\otimes I_n)\underline{z}\left(t\right),\hfill \end{array}\end{array}\right.$$

(15)

where $\overline{\mathrm{\Phi }}=\left[\begin{array}{c}{\mathrm{\Phi }}_1({\overline{z}}_1,{\underline{z}}_1)\\ ⋮\\ {\mathrm{\Phi }}_N({\overline{z}}_N,{\underline{z}}_N)\end{array}\right]$, $\underline{\mathrm{\Phi }}=\left[\begin{array}{c}{\mathrm{\Phi }}_1({\underline{z}}_1,{\overline{z}}_1)\\ ⋮\\ {\mathrm{\Phi }}_N({\underline{z}}_N,{\overline{z}}_N)\end{array}\right]$, $L=\left[\begin{array}{ccc}{L}_1& & \\ & \ddots & \\ & & L_N\end{array}\right]$, $\overline{Y}=\left[\begin{array}{ccc}{\overline{Y}}_1& & \\ & \ddots & \\ & & {\overline{Y}}_N\end{array}\right]$, $\underline{Y}=\left[\begin{array}{ccc}{\underline{Y}}_1& & \\ & \ddots & \\ & & {\underline{Y}}_N\end{array}\right]$.

Theorem 1. 

If Assumptions 1–4 hold, the $\overline{z}\left(t\right)$ and $\underline{z}\left(t\right)$ given in (15) satisfy $\underline{z}\left(t\right)\le z\left(t\right)\le \overline{z}\left(t\right)$.

Proof. 

The dynamics of the global error system are given by

$$\left\{\begin{array}{c}\begin{array}{cc}\hfill \dot{\overline{e}}\left(t\right)=& \dot{\overline{z}}\left(t\right)-\dot{z}\left(t\right)\hfill \\ \hfill =& [\hat{A}-L\hat{C}-\gamma M(\mathcal{L}\otimes I_n)]\overline{e}\left(t\right)+\overline{\mathrm{\Phi }}-\hat{\varphi }\left(z\left(t\right)\right)+\overline{Y}(\overline{e}\left(t\right)-\underline{e}\left(t\right)),\hfill \end{array}\\ \begin{array}{cc}\hfill \dot{\underline{e}}\left(t\right)=& \dot{z}\left(t\right)-\dot{\underline{z}}\left(t\right)\hfill \\ \hfill =& [\hat{A}-L\hat{C}-\gamma M(\mathcal{L}\otimes I_n)]\underline{e}\left(t\right)+\hat{\varphi }\left(z\left(t\right)\right)-\underline{\mathrm{\Phi }}+\underline{Y}(\overline{e}\left(t\right)-\underline{e}\left(t\right)).\end{array}\end{array}\right.$$

(16)

Define $ϵ\left(t\right)=\left[\begin{array}{c}\overline{e}\left(t\right)\\ \underline{e}\left(t\right)\end{array}\right]$, then we can obtain the derivative of $ϵ\left(t\right)$

$$\dot{ϵ}\left(t\right)=\mathrm{\Pi }ϵ\left(t\right)+\tilde{\mathrm{\Phi }},$$

(17)

where

$$\begin{array}{c}\hfill \mathrm{\Pi }=\left[\begin{array}{cc}\mathrm{\Psi }+\overline{Y}& \overline{Y}\\ \underline{Y}& \mathrm{\Psi }+\underline{Y}\end{array}\right],\tilde{\mathrm{\Phi }}=\left[\begin{array}{c}\overline{\mathrm{\Phi }}-\hat{\varphi }\left(z\left(t\right)\right)\\ \hat{\varphi }\left(z\left(t\right)\right)-\underline{\mathrm{\Phi }}\end{array}\right].\end{array}$$

(18)

The proof of the validity of (12) and (15) is equivalent to the proof of the non-negativity of the error system. From Lemma 4 and Remark 2, we have $\overline{\mathrm{\Phi }}-\hat{\varphi }\left(z\left(t\right)\right)\ge 0$ and $\hat{\varphi }\left(z\left(t\right)\right)-\underline{\mathrm{\Phi }}\ge 0$. It is obvious that $\overline{e}\left(0\right)\ge 0$ and $\underline{e}\left(0\right)\ge 0$ from Assumption 2. By Assumption 4, we have $\mathrm{\Pi }$ is Metzler. Above all, considering the monotone system theory, $ϵ\left(t\right)\ge 0$ holds, which implies $\underline{z}\left(t\right)\le z\left(t\right)\le \overline{z}\left(t\right)$ for all $t\ge 0.$ Hence, the proof of Theorem 1 is completed.

After proving the boundedness of DIO, we proceed to design the observer gain $L_i,{\overline{Y}}_i,{\underline{Y}}_i$, and $M_i$ to guarantee the stability of the DIO (15).

Theorem 2. 

Given a positive definite matrix $P=P^T$ and a constant $\tau >0$, if a solution exists such that

$$\tilde{\mathrm{\Omega }}=\left[\begin{array}{cc}He(P\overline{A}-Q_i\overline{C}+P{N}_1+{\overline{W}}_i)-2\tau I_n& P{N}_2-{\overline{W}}_i+{N_3}^{T}P+{\underline{W}}_i^T\\ {N_2}^{T}P-{\overline{W}}_i^T+P{N}_3+{\underline{W}}_i& He(P\overline{A}-Q_i\overline{C}+P{N}_4-{\underline{W}}_i)-2\tau I_n\end{array}\right]\le 0,$$

$$\gamma >\frac{\tau }{a\left(\mathcal{L}\right)},$$

(19)

where $L_i=P^{-1}{Q}_i$, ${\overline{Y}}_i=P^{-1}{\overline{W}}_i$, ${\underline{Y}}_i=P^{-1}{\underline{W}}_i$, and ${M_i}^{-1}=P$ are the observer gains, γ is the coupling strength, then the system (15) is a DIO of (14).

Proof. 

A Lyapunov function candidate can be defined as $V\left(t\right)={\sum }_{i=0}^N{r}_i{\overline{e}}_i^{T}P{\overline{e}}_i+{\sum }_{i=0}^N{r}_i{\underline{e}}_i^{T}P{\underline{e}}_i$; then, the derivative of $V\left(t\right)$ can be written as □

$$\begin{array}{cc}\hfill \dot{V}\left(t\right)=& \sum _{i=0}^N{r}_i{(\mathrm{\Gamma }{\overline{e}}_i+{\overline{\chi }}_i+{\overline{Y}}_i({\overline{e}}_i-{\underline{e}}_i))}^{T}P{\overline{e}}_i+\sum _{i=0}^N{r}_i{(\mathrm{\Gamma }{\underline{e}}_i+{\underline{\chi }}_i+{\underline{Y}}_i({\overline{e}}_i-{\underline{e}}_i))}^{T}P{\underline{e}}_i\hfill \\ \hfill & +\sum _{i=0}^N{r}_i{\overline{e}}_i^{T}P(\mathrm{\Gamma }{\overline{e}}_i+{\overline{\chi }}_i+{\overline{Y}}_i({\overline{e}}_i-{\underline{e}}_i))+\sum _{i=0}^N{r}_i{\underline{e}}_i^{T}P(\mathrm{\Gamma }{\underline{e}}_i+{\underline{\chi }}_i+{\underline{Y}}_i({\overline{e}}_i-{\underline{e}}_i))\hfill \\ \hfill \le & \sum _{i=0}^N{r}_i{(\mathrm{\Gamma }{\overline{e}}_i+N_1{\overline{e}}_i+N_2{\underline{e}}_i+{\overline{Y}}_i({\overline{e}}_i-{\underline{e}}_i))}^{T}P{\overline{e}}_i\hfill \\ \hfill & +\sum _{i=0}^N{r}_i{(\mathrm{\Gamma }{\underline{e}}_i+N_3{\overline{e}}_i+N_4{\underline{e}}_i+{\underline{Y}}_i({\overline{e}}_i-{\underline{e}}_i))}^{T}P{\underline{e}}_i\hfill \\ \hfill & +\sum _{i=0}^N{r}_i{\overline{e}}_i^{T}P(\mathrm{\Gamma }{\overline{e}}_i+N_1{\overline{e}}_i+N_2{\underline{e}}_i+{\overline{Y}}_i({\overline{e}}_i-{\underline{e}}_i))\hfill \\ \hfill & +\sum _{i=0}^N{r}_i{\underline{e}}_i^{T}P(\mathrm{\Gamma }{\underline{e}}_i+N_3{\overline{e}}_i+N_4{\underline{e}}_i+{\underline{Y}}_i({\overline{e}}_i-{\underline{e}}_i))\hfill \\ \hfill =& ((R\otimes (\overline{A}-L_i\overline{C})-\gamma R\mathcal{L}\otimes M_i)\overline{e}+R\otimes (N_1+{\overline{Y}}_i)\overline{e}+R\otimes (N_2-{\overline{Y}}_i)\underline{e})(I\otimes P)\overline{e}\hfill \\ \hfill & +{((R\otimes (\overline{A}-L_i\overline{C})-\gamma R\mathcal{L}\otimes M_i)\underline{e}+R\otimes (N_3+{\underline{Y}}_i)\overline{e}+R\otimes (N_4-{\underline{Y}}_i)\underline{e})}^T(I\otimes P)\underline{e}\hfill \\ \hfill & +{\overline{e}}^T(I\otimes P)((R\otimes (\overline{A}-L_i\overline{C})-\gamma R\mathcal{L}\otimes M_i)\overline{e}+R\otimes (N_1+{\overline{Y}}_i)\overline{e}+R\otimes (N_2-{\overline{Y}}_i)\underline{e})\hfill \\ \hfill & +{\underline{e}}^T(I\otimes P)((R\otimes (\overline{A}-L_i\overline{C})-\gamma R\mathcal{L}\otimes M_i)\underline{e}+R\otimes (N_3+{\underline{Y}}_i)\overline{e}+R\otimes (N_4-{\underline{Y}}_i)\underline{e})\hfill \\ \hfill =& {\overline{e}}^T(R\otimes He(P\overline{A}-P{L}_i\overline{C}+P(N_1+{\overline{Y}}_i))-\gamma (R\mathcal{L}+{\mathcal{L}}^{T}R)\otimes PM)\overline{e}\hfill \\ \hfill & +{\overline{e}}^T(R\otimes (P(N_2-{\overline{Y}}_i)+{(N_3+{\underline{Y}}_i)}^{T}P))\underline{e}+{\underline{e}}^T(R\otimes ({(N_2-{\overline{Y}}_i)}^{T}P+P(N_3+{\underline{Y}}_i)))\overline{e}\hfill \\ \hfill & +{\underline{e}}^T(R\otimes He(P\overline{A}-P{L}_i\overline{C}+P(N_4-{\underline{Y}}_i))-\gamma (R\mathcal{L}+{\mathcal{L}}^{T}R)\otimes PM)\underline{e},\hfill \end{array}$$

(20)

where $\mathrm{\Gamma }=\overline{A}-L_i\overline{C}-\gamma M_i{\sum }_{j=0}^N{\mathcal{L}}_{ij}$, ${\overline{\chi }}_i={\mathrm{\Phi }}_i({\overline{z}}_i,{\underline{z}}_i)-H{\varphi }_i\left(H^{-1}{z}_i\left(t\right)\right)$, and ${\underline{\chi }}_i=H{\varphi }_i\left(H^{-1}{z}_i\left(t\right)\right)-{\mathrm{\Phi }}_i({\underline{z}}_i,{\overline{z}}_i)$.

It follows from Lemma 2 that

$$2a\left(\mathcal{L}\right){\overline{e}}^{T}R\overline{e}\le {\overline{e}}^T(R\mathcal{L}+{\mathcal{L}}^{T}R)\overline{e},$$

(21)

$$2a\left(\mathcal{L}\right){\underline{e}}^{T}R\underline{e}\le {\underline{e}}^T(R\mathcal{L}+{\mathcal{L}}^{T}R)\underline{e}.$$

(22)

Considering $M_i=P^{-1}$, (21) and (22), (20) can be written as

$$\begin{array}{cc}\hfill \dot{V}\left(t\right)\le & {\overline{e}}^T(R\otimes He(P\overline{A}-P{L}_i\overline{C}+P(N_1+{\overline{Y}}_i))-2\gamma a\left(\mathcal{L}\right)R\otimes I_n)\overline{e}\hfill \\ \hfill & +{\overline{e}}^T(R\otimes (P(N_2-{\overline{Y}}_i)+{(N_3+{\underline{Y}}_i)}^{T}P))\underline{e}+{\underline{e}}^T(R\otimes ({(N_2-{\overline{Y}}_i)}^{T}P+P(N_3+{\underline{Y}}_i)))\overline{e}\hfill \\ \hfill & +{\underline{e}}^T(R\otimes He(P\overline{A}-P{L}_i\overline{C}+P(N_4-{\underline{Y}}_i))-2\gamma a\left(\mathcal{L}\right)R\otimes I_n)\underline{e}.\hfill \end{array}$$

(23)

Noting that $\gamma >\frac{\tau }{a\left(\mathcal{L}\right)}$, we have

$$\begin{array}{cc}\hfill \dot{V}\left(t\right)\le & {\overline{e}}^T(R\otimes He(P\overline{A}-P{L}_i\overline{C}+P(N_1+{\overline{Y}}_i))-2\tau R\otimes I_n)\overline{e}\hfill \\ \hfill & +{\overline{e}}^T(R\otimes (P(N_2-{\overline{Y}}_i)+{(N_3+{\underline{Y}}_i)}^{T}P))\underline{e}+{\underline{e}}^T(R\otimes ({(N_2-{\overline{Y}}_i)}^{T}P+(P{N}_3+{\underline{Y}}_i)))\overline{e}\hfill \\ \hfill & +{\underline{e}}^T(R\otimes He(P\overline{A}-P{L}_i\overline{C}+P(N_4-{\underline{Y}}_i))-2\tau R\otimes I_n)\underline{e}\hfill \\ \hfill =& {ϵ}^T\left(t\right)R\otimes \mathrm{\Omega }ϵ\left(t\right),\hfill \end{array}$$

(24)

where $ϵ\left(t\right)={[{\overline{e}}^T,{\underline{e}}^T]}^T$, and

$$\Omega =\left[\begin{array}{cc}He(P\overline{A}-P{L}_i\overline{C}+P(N_1+{\overline{Y}}_i))-2\tau I_n& P(N_2-{\overline{Y}}_i)+{(N_3+{\underline{Y}}_i)}^{T}P\\ {(N_2-{\overline{Y}}_i)}^{T}P+P(N_3+{\underline{Y}}_i)& He(P\overline{A}-P{L}_i\overline{C}+P(N_4-{\underline{Y}}_i))-2\tau I_n\end{array}\right].$$

(25)

Then, $\dot{V}\left(t\right)<0$ and $\mathrm{\Omega }\prec 0$ are equal, which shows that $li{m}_{t\to \infty }\overline{e}\left(t\right)=0$ and $li{m}_{t\to \infty }\underline{e}\left(t\right)=0$. The stability of DIO can be guaranteed.

Remark 4. 

To satisfy the LMI toolbox, $Q_i=P{L}_i,{\overline{W}}_i=P{\overline{Y}}_i,{\underline{W}}_i=P{\underline{Y}}_i$ are applied to Ω, which results in

$$\tilde{\mathrm{\Omega }}=\left[\begin{array}{cc}He(P\overline{A}-Q_i\overline{C}+P{N}_1+{\overline{W}}_i)-2\tau I_n& P{N}_2-{\overline{W}}_i+{N_3}^{T}P+{\underline{W}}_i^T\\ {N_2}^{T}P-{\overline{W}}_i^T+P{N}_3+{\underline{W}}_i& He(P\overline{A}-Q_i\overline{C}+P{N}_4-{\underline{W}}_i)-2\tau I_n\end{array}\right].$$

4. Simulation

In this section, a team of 2-DOF manipulator models with an directed graph is applied to validate the feasibility of the presented DIOs for ELSs. The example runs on a CPU Inter Core i7-10750H CPU with 2.59 GHZ and 8 GB of RAM; the main software used is MATLAB 2018a.

The parameters of the model can be referred to [12], and the description of the parameters is shown in

Table 1. The descriptions of the parameters.

Parameter	Description $(\mathit{i}=1,2,3,4)$
$q_i$	The position of manipulator $i$
${\dot{q}}_i$	The velocity of manipulator $i$
${\ddot{q}}_i$	The acceleration of manipulator $i$
$M_i$	The inertia matrix of manipulator $i$
$C_i$	The vector of Coriolis and centrifugal force of manipulator $i$
$G_i$	The vector of Gravitational force of manipulator $i$
$u_i$	The control input of manipulator $i$

. First of all, the matrix $\mathcal{L}$ is

$$\mathcal{L}=\left[\begin{array}{cccc}3& -1& -1& -1\\ -1& 2& -1& 0\\ -1& -1& 2& 0\\ -1& 0& 0& 1\end{array}\right],$$

From Lemma 2, it follows that $\alpha \left(\mathcal{L}\right)=1$. Then, we chose the coordinate transformation matrix as

$$H=\left[\begin{array}{cccc}1& 0& 0& 0\\ 0& 1& 0& 0\\ -1& 0& 1& 0\\ 0& -1& 0& 1\end{array}\right].$$

By solving (19), the gain of the observer can be calculated.

$${\overline{W}}_i=\left[\begin{array}{cccc}1.7448& 0& 0.3703& 0\\ 0& 1.7448& 0& 0.3703\\ 0.3703& 0& 1.0920& 0\\ 0& 0.3703& 0& 1.0920\end{array}\right],{\underline{W}}_i=\left[\begin{array}{cccc}-1.7448& 0& -0.3703& 0\\ 0& -1.7448& 0& -0.3703\\ -0.3703& 0& -1.0920& 0\\ 0& -0.3703& 0& -1.0920\end{array}\right],$$

$$Q=\left[\begin{array}{cc}5.3976& 0\\ 0& 5.3976\\ 3.0250& 0\\ 0& 3.0250\end{array}\right],M_i=\left[\begin{array}{cccc}0.5235& 0& 0.0920& 0\\ 0& 0.5235& 0& 0.0920\\ 0.0920& 0& 0.8917& 0\\ 0& 0.0920& 0& 0.8917\end{array}\right],$$

and

$$P=\left[\begin{array}{cccc}1.9456& 0& -0.2008& 0\\ 0& 1.9456& 0& -0.2008\\ -0.2008& 0& 1.1422& 0\\ 0& -0.2008& 0& 1.1422\end{array}\right].$$

Noting Remark 4, we can obtain

$${\overline{Y}}_i=\left[\begin{array}{cccc}0.9474& 0& 0.2943& 0\\ 0& 0.9474& 0& 0.2943\\ 0.4908& 0& 1.0078& 0\\ 0& 0.4908& 0& 1.0078\end{array}\right],{\underline{Y}}_i=\left[\begin{array}{cccc}-0.9474& 0& -0.2943& 0\\ 0& -0.9474& 0& -0.2943\\ -0.4908& 0& -1.0078& 0\\ 0& -0.4908& 0& -1.0078\end{array}\right],$$

$$L_i=\left[\begin{array}{cc}3.1040& 0\\ 0& 3.1040\\ 3.1941& 0\\ 0& 3.1941\end{array}\right].$$

The initial value is defined as $x\left(0\right)={\left[0564415634616415\right]}^T$. The initial values of DIO are chosen as $\overline{x}\left(0\right)={\left[101097101013910107131010810\right]}^T$ and $\underline{x}\left(0\right)={[-4-4-4.5-3.5-4-4-5-1-4-4-4-6-4-4-1-4]}^T$, and $\gamma$ is given as $\gamma =3$.

The simulation results are displayed in the figures below. It is worth noting that i in $v_{ij}$ refers to the i-th manipulator and j represents the j-th state. Figure 1, Figure 2, Figure 3 and Figure 4 depict the original state trajectory of velocity and the state trajectory of DIO at the two joints of the four manipulators, respectively. It can be seen that the boundaries of the DIO recover the state of the original system. Figure 5, Figure 6, Figure 7 and Figure 8 show the errors in the DIO and the original system velocity for the four manipulators at the two joints, respectively. From the figures, we can see that the error between the original system and observer converge to a small bound, which means that the DIO is valid.

Remark 5. 

In the simulation, we first need to design a coordinate transformation matrix and then solve $\tilde{\mathrm{\Omega }}<0$ (25) by using the LMI toolbox, and the feasible solution is applied to the simulation model we have built. The complexity of the computation is reflected in solving the LMI, and the details are shown in

Table 2. Complexity of the calculation.

Parameter	Description	Time
${\overline{W}}_i$	Observer gain	0.436 s
${\underline{W}}_i$	Observer gain	0.436 s
Q	Observer gain	0.456 s
$M_i$	Weight matrix	0.318 s
P	Weight matrix	0.318 s
${\overline{Y}}_i$	Observer gain	0.624 s
${\underline{Y}}_i$	Observer gain	0.624 s
$L_i$	Observer gain	0.624 s

.

It is worth noting that the parameters in Table 2 all need to be recalculated if the topology or the coordinate transformation matrix changes.

Remark 6. 

Compared with [22], MELSs contain nonlinearities, which makes it more complicated to design the IO. That is, the design method proposed in the paper can also be applied to more general systems. Unlike [24], two observer gains ${\overline{Y}}_i$ and ${\underline{Y}}_i$ are considered to make the the sufficient condition less conservative while ensuring the existence of the interval observer. As shown in

Table 3. The solution under different $\tau$ [24].

Value of $\mathit{\tau }$	Theorem 2	Theorem 2 in [24]
−100	✓	-

, the conditions obtained using the method in this article are solvable. However, using the method in [24], there is no solution.

5. Conslusions

In this paper, a DIO design method for MELSs is studied. The design conditions are relaxed by using the coordinate transformation method. Two methods of designing observer gains are formulated under monotone system theory and Lyapunov stability theory. Meanwhile, the additional observer gains $\overline{Y}$ and $\underline{Y}$ are introduced in DIO, which also provide more feasibility for the existence of the observer. Finally, an example is given to verify the effectiveness of the presented DIO design method. In the future, we will try to relax the requirements of the topology graph and study the design of DIO-based controller algorithms.