Skip to main content

ORIGINAL RESEARCH article

Front. Aerosp. Eng., 10 May 2023
Sec. Navigation, Guidance, and Control
Volume 2 - 2023 | https://doi.org/10.3389/fpace.2023.1158718

Interval observers design for systems with ostensible Metzler system matrices

www.frontiersin.orgDušan Krokavec* www.frontiersin.orgAnna Filasová
  • Department of Cybernetics and Artificial Intelligence, Faculty of Electrical Engineering and Informatics, Technical University of Košice, Košice, Slovakia

This paper attempts to resolve the problem concerning the interval observers design for linear systems with ostensible Metzler system matrices. Because system dynamics matrices are partially different from strictly Metzler structures, a solution is achieved by constructing a composed system matrix representation, which combines pre-compensated interval matrix structures fixed with a prescribed region of D-stability and the reconstructed strictly Metzler matrix structure, related to the original interval system matrix parameter definition. A novel design procedure is presented, which results in a strictly positive observer gain matrix and guarantees that the lower estimates of the positive state variables are non-negative when considering the given system structure and the non-negative system state initial values. The design is computationally simple since it is reduced to the feasibility of the set of linear matrix inequalities.

1 Introduction

Interval observers have appeared as an alternative technique for robust state estimation (Moisan et al., 2009). Whilst, when using the technique based on classical observers, only the initial condition is assumed to be unknown (Luenberger, 1971), interval observers structures are constructed assuming that the upper and lower bounds of the initial conditions are known (Raïssi and Efimov, 2018; Khan et al., 2020). The main limitation to the interval observers theory is that the trajectories of the system that start from an internally bounded initial condition will enclose the stable system trajectory only if the system is positive that its system matrix is Metzler and Hurwitz and that other matrix parameters are non-negative (Farina and Rinaldi, 2000). Thus, the positivity of interval estimation error dynamics is one of the most restrictive assumptions for interval observers design. When restricted to the Metzler structure of system matrices, as well as to non-negative input and output matrices, such systems are referred to as Metzler systems (Nikaido, 1968; Smith, 1995; Liu et al., 2011), with a stringent approach that reflects the diagonal stabilization principle. Although a certain class of systems can be transformed through a change of coordinates into positive cooperative systems (Mazenc and Bernard, 2011; Mazenc and Bernard, 2014), no general technique exists for such construction.

When maintaining platforms for positive systems with nonnegative states (Nikaido, 1968; Smith, 1995; Moisan et al., 2009), the theory of Metzler matrices (Berman et al., 1989) implies some additional parametric constraints to reflect the system positiveness (Shorten et al., 2009) and to construct the system representation (Son and Hinrichsen, 1996; Gao et al., 2005; Liu et al., 2017; Ito and Dinh, 2020). Since the linear time-invariant system theory cannot be directly used for linear positive systems, various combinations of linear programming and linear matrix inequalities (LMI) are generally used to represent Metzler systems (Ait Rami and Tadeo, 2006; Shu et al., 2008; Anderson and Murray, 2018; Guo et al., 2020). The benefits of a potential unification are presented (Krokavec and Filasová, 2018) when reflecting diagonal stabilization and associated Metzler system matrix parametric representations by a specific set of LMIs.

The system matrix parametric constraints give rise to substantially complex design methods when applied to positive systems with interval-defined model parameters (Ganesan, 2007). To demarcate the object of study in this field, Metzler matrix transforms are reflected for interval observers analysis (Mason, 2012; Chambon et al., 2015). Interval observers design for linear time-varying (LTV) systems, as well as for a class of non-linear time-varying systems with output specifications, exploits static coordinate transformation (Efimov et al., 2013; Raïssi and Efimov, 2018) when translating a stable LPV system to another stable and cooperative LPV system. The LMI-based conditions applicable in interval observers design for positive Metzler systems have also been studied (Krokavec and Filasová, 2020b). The utilization of interval observers for interconnected schemes is often applied in relation to distributed interval estimation and distributed feedback control (Wang et al., 2020; Wang X L et al., 2022; Zhang et al., 2022); this also reflects that their application for continuous linear large-scale systems is limited due to the system’s complexity (Wang T et al., 2022). These problems are still open in distributed applications since it is difficult to ensure that the system state will be enclosed by the cooperative estimated upper and lower bounds of the observed system state (Huong, 2022; Li et al., 2022), as well as in interval estimation strategy for anti-disturbance control of drones (Yong et al., 2020).

This paper contributes to the properties of the interval state estimation for linear systems with ostensible Metzler system matrices. It outlines a new LMI-based approach to determine interval observers with positive observer gains using a combined representation of the ostensible Metzler system matrix (Krokavec and Filasová, 2022). Design conditions are formulated using LMIs, respecting the diagonal stabilization principle, Metzler system matrix parametric constraints, and given interval matrix bounds. Because linear systems with ostensible Metzler matrices are not positive, even if their matrix parameters are non-negative, the main limitation to the solution is that the interval estimation only works for system state variables whose trajectories for a non-negative initial state are non-negative.

The outline of this paper is as follows. Following an introduction in Section 1, the basic preliminaries are discussed in Section 2. Section 3 presents the LMI structures necessary for observer stability and the positive gain, and the design method of interval observers for a given class of positive systems is presented in Section 4. To illustrate the design task, its efficiency is demonstrated by numerical solutions in Section 5; in Section 6, conclusions are briefly presented.

Throughout the paper, X ≺ 0 conveys briefly that a real square matrix X is a symmetric and negative definite, notations xT and XT identify the transpose of a vector or a matrix, In indicates the nth order unit matrix, function ρ() reflects the eigenvalue spectrum of a real square matrix, diag[ ] enters a block diagonal matrix, the symbol * is used as ellipsis in a symmetric matrix, L is the set of real function z(t) with the property ‖z(t)‖ < , the relations x1x2 and X1X2 operate on corresponding elements element-wise, R(R+) qualifies the set of (non-negative) real numbers, Rn×r(R+n×r) refers to the set of n × r (non-negative) real matrices, and M+n×n denotes the set of strictly Metzler square matrices.

2 Basic preliminaries

To explain the technique used, the main question can be illustrated using the linear Metzler systems given as follows:

q̇t=Aqt+But+Ddt,(1)
yt=Cqt,(2)

where q(t)R+n, u(t)Rr, and y(t)R+m are non-negative vectors of the system, input and output, disturbance d(t)L, d(t)R+rd is norm-bounded and non-negative, and BR+n×r, DR+n×rd, and CR+m×n are non-negative matrices.

Definition 1. Berman et al. (1989). A square matrix AM+n×n is a strictly Metzler matrix if all its diagonal elements are negative and all its off-diagonal elements are positive.

This study also uses the term “a purely Metzler matrix AMn×n” if all diagonal elements of A are negative and all its off-diagonal elements are non-negative, and “an ostensible Metzler matrix AMn×n” if all diagonal elements of A are negative and at least one off-diagonal element is negative while the number of non-negative off-diagonal elements is prevalent.

2.1 Positive continuous-time linear systems

The following assumptions make it possible to cover constraints on the parameters of system (1) and (2) when studying the system positivity and the conditions of its diagonal stabilizability.

Lemma 1. Tanaka and Langbort (2011). Disturbance-free system (1), (2) is internally positive if and only if A is (strictly, purely) Metzler and B, C are entry-wise non-negative.

In consequence, any solution of an autonomous and disturbance-free linear system with a Metzler matrix AMn×n is element-wise non-negative for all t ≥ 0, provided that q (0) ≥0. The output solution y(t) for such a defined solution is non-negative if CR+m×n. Such dynamical systems are called non-negative only if initial conditions in R+n are considered.

Remark 1. A strictly Metzler matrix A={aij}I,j=1n makes of n2 constraints

aii<0i=1,,n,aij,ij>0i,j=1,,n.(3)

These parametric constraints imply the strict application of the diagonal stabilization principle (Shorten et al., 2009; Mason, 2012) in analysis. If a strictly Metzler AM+n×n is represented (in relation to the observer design) in the following rhombic form, where the diagonal items reflected by the column index define multiple circular shifts of elements of the columns of A as (Krokavec and Filasová, 2020b)

AΘ=a11a21a22a31a32a33an1an2an3anna12a13a1na23a2nan1,n,(4)

the diagonal stabilization principle can be appropriately respected using the derived diagonal matrix structures related to AΘ as

Al,l=diaga11a22ann0,(5)
Al+h,l=diaga1+h,1an,nha1,nh+1ah,n0,(6)

for h = 0, 1, , n−1.

Remark 2. Defining the matrix LRn×n in the circulant permutation form (Horn and Johnson, 1995)

L=0T1In10,L1=LT(7)

and considering a diagonal matrix Z=diagz11z22znn, where ZRn×n, then

LTZL=diagz22znnz11.(8)

The aforementioned results can be combined and reflected by the following lemma.

Lemma 2. Krokavec and Filasová (2020a) If a positive matrix JR+n×m forces the strictly Metzler matrix Ae=AJCR+n×n, where AR+n×n is strictly Metzler and CR+m×n is non-negative, then Ae is parameterized as

Ae=h=0n1Al+h,lk=0rJjhCkLhT,(9)
Al,lk=0rJkCk0,(10)
Al+h,lk=0rJkhCkLhT0,(11)

where h = 1, , n − 1 and the diagonal matrices Jk,Jkh,CkR+n×n are composed in the following ways:

C=c1TcmT,Ck=diagckT,(12)
J=j1jm,Jk=diagjk,Jkh=LhTJkLh.(13)

The proof of the aforementioned lemma is based on the fact that only diagonal matric representations are applicable for the diagonal stabilization of positive systems.

2.2 Ostensible Metzler matrices

Given a system with the dynamical model (1) and (2) and considering that ARn×n is an ostensible Metzler matrix, then inclusion of the negative off-diagonal elements of A into the design task is built on the following basic facts from the theory of matrices.

Definition 2. Shores (2007) Matrix XRn×n is similar to the matrix ΛRn×n if there exists an invertible similarity transform matrix SRn×n such that

S1XS=Λ.(14)

If X and Λ are similar, then they have the same eigenvalues, their algebraic multiplicities are the same, their characteristic polynomials are the same, and their determinants and traces are the same.

Remark 3. Let {vkCn}k=1n be the set of eigenvectors for a matrix XRn×n and {λkC}k=1n is the associated set of eigenvalues of X such that eigenvalues are all distinct, then (14) implies

V=v1v2vn,Λ=diagλ1λ2λn(15)

and {vkCn}k=1n are linearly independent.

Theorem 1. Shores (2007) If for X,YRn×n it can be set Y = cX + dIn with scalars c,dR, c ≠ 0, and InRn×n, then the eigenvalues of Y are

ξk=cλk+d,(16)

where λk runs over ρ(X) with k = 1, , n, and the eigenvectors of X and Y are identical.Supposing that A is ostensible Metzler, then the proposed idea means decoupling the system matrix A so that A = Ap + Am, where Ap is strictly Metzler and Am is entry-wise negative and Hurwitz.Lemma 3. Krokavec and Filasová (2022) A strictly Metzler ApM+n×n and an entry-wise negative and Hurwitz AmRn×n to the composed form of the ostensible Metzler matrix A=Ap+AmMn×n exist if there exist positive scalars η,δR+ such that with

λo=maxkλk+|λk+=realλk>0,λkρA°m,(17)
p=λo+δ,Ad+pIn0,Ad=diaga11ann,A=aiji,j=1n,(18)

it yields

Ap=Ad+A++ηΣ+pIn=A°p+pIn,Am=AηΣpIn=A°mpIn,(19)

where

Σ=011110111110(20)

and

A=0a12a1,n1a1na210a2,n1a2nan1,1an1,20an1,nan1an2an1,n10,A+=0a12+a1,n1+a1n+a21+0a2,n1+a2n+an1,1+an1,2+0an1,n+an1+an2+an1,n1+0,(21)
aij+=aijifaij>0,0ifaij<0,aij=aijifaij<0,0ifaij>0,(22)

are defined for i, j ∈ ⟨1, n⟩, ij, whilst ρ(A°m) is the set of eigenvalues of the matrix A°m.Remark 4. Structure (21) implies, since the sum of the eigenvalues of a matrix equals its trace,

trA°m=k=1nλkm=0(23)

and so the set {λkmC}k=1n consists of stable and unstable eigenvalues. For repeated eigenvalues, one must add them according to their multiplicity, but this does not change the existence of a real eigenvalue with a maximal positive value to be compensated by the construction used.Remark 5. An upper-bound parameter ηR+ must be chosen such that, in the final result, it must be set λo > η. Since (19) is constructed as a strictly Metzler matrix, all the elements on its main diagonal must be negative. This implies the boundary condition in defining a stable D-stability region with δR+ such that

maxiaii+p:aii+pAd+pIn<0,p=λo+δ.(24)

2.3 Intervally defined ostensible Metzler matrices

In this case, is assumed that q(0) and the ostensible Metzler system parameter A are unknown but bounded by constant bounding vectors and constant bounding matrices of appropriate dimensions in such a way that (these inequalities being understood element-wise (Jaulin et al., 2001))

0q̲0q0q̄0,A̲AĀ.(25)

Since the main goal is the design of an interval observer of the state, it is considered that BR+n×r, DR+n×rd, and CR+m×n are known non-negative matrices, the system input u(t) is norm-bounded, and the following assumption is adopted.

Assumption 1. The function bounds d̲,d̄R+rd and d̲,d̄L are given such that

0d̲dtd̄.(26)

This assumption states that the disturbance is known up to some interval error ed=d̄d̲.

Corollary 1. Strictly Metzler A̲p,ĀpM+n×n and an entry-wise negative and Hurwitz A̲m,ĀmRn×n to the composed forms of the ostensible Metzler matrices A̲=A̲p+A̲m,Ā=Āp+ĀmMn×n exist if there exist positive scalars η̲,η̄,δ̲,δ̄R+ such that, with

λ̲o=maxkλ̲k+|λ̲k+=realλ̲k>0,λ̲kρA̲°m,(27)
λ̄o=maxkλ̄k+|λ̄k+=realλ̄k>0,λ̄kρĀ°m,(28)
p̲=λ̲o+δ̲,A̲d+p̲In0,A̲d=diaga̲11a̲nn,A̲=a̲iji,j=1n,(29)
p̄=λ̄o+δ̄,Ād+p̄In0,Ād=diagā11ānn,Ā=āiji,j=1n,(30)

it yields

A̲p=A̲d+A̲++η̲Σ+p̲In=A̲°p+p̲In,A̲m=A̲η̲Σp̲In=A̲°mp̲In,(31)
Āp=Ād+Ā++η̄Σ+p̄In=Ā°p+p̄In,Ām=Āη̄Σp̄In=Ā°mp̄In,(32)

where Σ is from (20) and A̲+, Ā+, A̲, and Ā are constructed by the rules defined in (21) and (22).

3 General interval observers structure

Under these introduced assumptions, the interval observers equations for systems with intervally given ostensible Metzler matrices can be defined as follows:

q̲̇et=A̲q̲et+But+Jyty̲et=A̲eq̄et+But+Jyt,(33)
q̇̄et=Āq̄et+But+Jytȳet=Āeq̄et+But+Jyt,(34)

where q̲e(t)Rn and q̄e(t)Rn are, respectively, the lower and upper interval estimates for the state q(t) and

Āe=ĀJC,A̲=A̲JC,(35)
y̲t=Cq̲t,ȳt=Cq̄t,y̲et=Cq̲et,ȳet=Cq̄et.(36)

Using the observation errors

e̲t=qtq̲et,ēt=qtq̄et,(37)

it follows from (1), (33), and (34) that

ė̲t=A̲ee̲t+Ddt,ė̄t=Āee̲t+Ddt.(38)

To construct a Hurwitz stable A̲e,ĀeRn×n, guaranteeing also strictly Metzler and Hurwitz matrices A̲pe,ĀpeM+n×n when implementing for ostensible Metzler A̲,ĀMn×n, A̲=A̲p+A̲mMn×n, and Ā=Āp+ĀmMn×n, then (38) can be rewritten as

ė̲t=A̲pee̲t+A̲me̲t+Ddt,ė̄t=Āpee̲t+Āme̲t+Ddt,(39)

where

A̲pe=A̲pJC,Āpe=ĀpJC,A̲e=A̲pe+A̲m,Āe=Āpe+Ām.(40)

To apply the parametrization principle in designing this class of observer, the following corollary is objective.

Corollary 2. State observation error dynamics (40) entail the parameterizations of the strictly Metzler matrices A̲p and Āp as follows:

A̲pl,l=diaga̲p11a̲pnn,Āpl,l=diagāp11āpnn,(41)
A̲pl+h,l=diaga̲p,1+h,1a̲p,n,nha̲p,1,nh+1a̲ihn,(42)
Āpl+h,l=diagāp,1+h,1āp,n,nhāp,1,nh+1āphn,(43)
A̲pe=h=0n1A̲pl+h,lj=0rJjhCjLhT,(44)
Āpe=h=0n1Āpl+h,lj=0rJjhCjLhT,(45)

while the parameterizations (12) and (13) stay unchanged.

Provided that (25) is satisfied, then for all tR+, the estimates q̲e(t) and q̄e(t) are bounded with the limit properties, illustrated by the following remark.

Remark 6. Performing an inner adjustment for (38) as

q̇tq̲̇et=A̲eqtq̲et+Ddt,(46)
q̲̇et=q̇tA̲eq̲t+A̲eq̲etDdt,(47)

respectively, and substituting (1) in (47) yields

q̲̇et=AA̲JCqt+A̲eq̲et+But=AA̲q̲t+JCq̲t+A̲eq̲et+But=AA̲q̲t+JCq̲t+A̲ep+A̲mq̲et+But,(48)

and, if A̲e=A̲ep+A̲m is Hurwitz, CR+m×n is nonnegative, JR+n×m is positive, and A̲A, then the lower system state estimate produced by the interval observer constructed on the system model with the ostensible Metzler matrix converges to a non-negative trajectory if JCq̲(t)>0. Consequently, provided that q̲(0)q(0)q̄(0), then for all tR+, the estimates q̲(t) and q̄(t) given by (33) and (34) produce the interval bounds only to those system state variables qi(t), i = 1, , n which are non-negative.

4 Interval observers design

The design goals are Hurwitz stable matrices A̲eRn×n and ĀeRn×n and strictly Metzler and Hurwitz matrices A̲peM+n×n and ĀpeM+n×n when implementing for ostensible Metzler A̲Mn×n and ĀMn×n. A solution method, resulting in positive matrix gain JR+n×m, is given in Theorem 2.

Theorem 2. The matrices A̲ep,Āep,R+n×n are strictly Metzler and Hurwitz and the matrices Āe and A̲eRn×n are Hurwitz if, for the given ostensible Metzler matrices A̲,ĀRn×n and non-negative CR+m×n, there exist positive definite diagonal matrices P,RkR+n×n and positive scalars μ̲,μ̄R+ that for h = 1, , n − 1, lT=11 and the parameters from Corollary 2 satisfy the LMIs

P0,Rk0,(49)
Ω̲**DTPμ̲Ird*C0μ̲Im0,Ω̄**DTPμ̄Ird*C0μ̄Im0,(50)
PĀpl,lk=1mRkCk0,PA̲pl,lk=1mRkCk0,(51)
PLhA̲pl+h,lLhTk=1mRkLhCkLhT0,PLhĀpl+h,lLhTk=1mRkLhCkLhT0,(52)
Ω̲=PA̲p+A̲pTP+PA̲m+A̲mTPk=1rRkllTCk+CkllTRk,(53)
Ω̄=PĀp+ĀpTP+PĀm+ĀmTPk=1rRkllTCk+CkllTRk.(54)

Confirming the feasible task, the interval observer gain is given as

Jk=P1Rk,jk=Jkl,J=j1jm.(55)

Proof. To respect the diagonal stabilization principle, v(e̲(t)) is served as a Lyapunov function for (37) using a symmetric positive definite matrix PR+n×m and a positive scalar µ̲R+ such that

ve̲t=e̲TtPe̲t+μ̲10te̲yTτe̲yτμ̲2dTτdτdτ>0,(56)

whose time-derivative for the observer error trajectory must satisfy

v̇e̲t=ė̲TtPe̲t+e̲TtPė̲t+μ̲1e̲yTte̲ytμ̲dTtdt<0.(57)

Applying in inequality (57) the observer error dynamics (37) gives the following:

v̇e̲t=e̲TtA̲eTP+PA̲ee̲t+e̲TtPDdt+dTtDTPe̲t+μ̲1e̲TtCTCe̲tμ̲dTtdt<0.(58)

Thus, constructing a common notation e̲d(t) that is readily representable for the used variables as

e̲dTt=e̲TtdTt,(59)

then there is reasonable grounds to conclude that

v̇e̲dt=e̲dTtΩ̲°e̲dt<0,(60)

where, for the covered systematization,

Ω̲°=A̲eTP+PA̲e+μ̲1CTCPDDTPμ̲Ird0.(61)

Therefore, the new form of LMI after applying the property of the Schur complement is

PA̲e+A̲eiTP**DTPμ̲Ird*C0μ̲Im,0(62)

and, using (13) and 40, it can be set as

PA̲pe=PA̲pJC=PA̲pk=1mPjkckT=PA̲pk=1mPJkllTCk,(63)

where the column vector l is used to uncover the diagonal matrix structures. Thus, (62) implies (50) and (53) when substituting

PJk=Rk,A̲e=A̲pe+A̲m.(64)

Separating h = 0 from (44) diagonal part and multiplying its left side by P yields

PA̲pl,lj=1mPJjCj0(65)

and using notation (64) then (65) implies (51). Analogously, it can be obtained when taking from (44) a component for h ≠ 0 and multiplying its left side by PLh (since LhLhT = In) that

PLhA̲pl,l+hLhTk=1mPLhLhTJkLhCkLhT0(66)

and, using notation (64), then (66) implies (52).Analogously, all this can be carried out for the upper bound parameters. This concludes the proof.

5 Illustrative examples

In this section, two examples are presented to demonstrate the effectiveness of the interval observers design.

Example 1. To illustrate the proposed design principles, the stable interval ostensible strictly Metzler systems (1) and (2) are constructed on the matrices

A̲=0.1721.941.450.1421.960.380.1000.172.91,Ā=0.1582.061.550.1421.640.320.2000.172.55,D=0.121.090.21,CT=100100.

To apply Theorem 2 conditions, the derived design parameters are selected as

A̲d=diag0.1721.962.91,Ād=diag0.1581.642.55,

and the related matrix structures are constructed from the system matrix bounds as follows:

A̲=0000.14200.38000,A̲+=01.941.450000.100.170,Σ=011101110,
Ā=0000.14200.32000,Ā+=02.061.550000.200.170.

Thus, using Σ and η = 0.005 yields for A°m that

A̲°m=A̲ηΣ=00.0050.0050.14700.3250.0050.0050,ρA̲°m=0.05410.00500.0491,λ̲0=0.0491,
Ā°m=ĀηΣ=00.0050.0050.14700.3250.0050.0050,ρĀ°m=0.05110.00500.0461,λ̄0=0.0461.

Setting δ̲=δ̄=0.003 means p̲=λ̲0+δ̲=0.0521 and p̄=λ̄0+δ̄=0.0491, and so A̲m=A̲°mp̲In, Ām=Ā°mp̄In take the values

A̲m=0.05210.00500.00500.14700.05210.38500.00500.00500.0521,ρA̲m=0.10620.04710.0030,
Ām=0.07610.00500.00500.14700.07610.32500.00500.00500.0761,ρĀm=0.12720.07110.0300.

Furthermore, A̲p=A̲d+A̲++ηΣ+p̲In and Āp=Ād+Ā++ηΣ+p̄In are computed as

A̲p=0.11991.94501.45500.00501.90790.00500.10500.17502.8579,ρA̲p=0.05961.91332.9128,
Āp=0.06832.06501.55500.00501.55030.00500.20500.17502.4603,ρĀp=0.06521.55722.5869,

which are strictly Metzler, and their rhombic representations imply the diagonal matrices for the observer synthesis

A̲pl,l=diag0.11991.90792.8579,Āpl,l=diag0.06831.55032.4603,
A̲pl+1,l=diag0.00500.17501.4550,Āpl+1,l=diag0.00500.17501.5550,
A̲pl+2,l=diag0.10501.94500.0050,Āpl+2,l=diag0.20502.06500.0050,

whilst straightforward calculations give

C1=diag100,C2=diag010,lT=1111,L=0T1I30.

Using LMIs defined by Theorem 2, the feasible matrix variables result in the non-negative gain matrix when applying the SeDuMi package (Peaucelle et al., 2002)

P=diag2.42343.12892.71250,R1=diag3.08450.00480.09880,R2=diag2.36271.15040.14880,J=1.27280.97500.00150.36770.03640.0548,μ̲=4.0252,μ̄=4.1147.

This infuses the strictly Metzler and Hurwitz matrices A̲pe=A̲pJC and Āpe=ĀpJC as

A̲pe=1.39270.97001.45500.00352.27560.00500.06860.12022.8579,ρA̲pe=1.32352.27942.9233,
Āpe=1.34111.09001.55500.00351.91800.00500.16860.12022.4603,ρĀpe=1.13661.92362.6591,

where A̲peĀpe. In addition, it can be seen that, due to the structure of matrix C, the elements on the third columns of matrices A̲pe and Āpe have not changed compared to A̲p and Āp.

Applying the same gain matrix to the ostensible Metzler matrices yields

A̲e=A̲JC=1.44480.96501.45000.14352.32770.38000.06360.11522.9100,
Āe=ĀJC=1.43081.08501.55000.14352.00770.32000.16360.11522.5500.

It should be noted that the positions of the negative off-diagonal elements in A̲ and A̲e, as well as in Ā and Āe, have been preserved. In addition, in the considered case, A̲eĀe.

By simulating the response of the autonomous system with considered interval ostensible Metzler parameters to better illustrate the ostensible Metzler phenomena, the dynamics of the system were

A=0.1652.001.500.1421.800.350.1500.172.73,A̲AĀ,

the initial system state was set as q(0) = [ 0.5 7.5 0 ]T, q̲e(0)=q̄e(0)=0, and σd2=0.04. The simulation is executed in the MATLAB framework using Simulink.

Figure 1 depicts the time responses of the first system state variable and its upper and lower estimations; Figure 2 shows the time responses for the third system state variable. Although the given system is not positive, it can be seen from Figures 1 and 2 that the behaviors of these state variables are correctly intervally estimated by the proposed interval observer if the components A̲mq(t) and Āmq(t), indicated in (48), are compensated by suitably choosing the observer initial states, satisfying conditions q̲e(0)q(0)q̄e(0). Since the state variable q2(t) is undefined in sign, its interval estimation is also undefined in sign. This case is trivial and is not presented.

FIGURE 1
www.frontiersin.org

FIGURE 1. Estimation of the first state variable in the autonomous mode.

FIGURE 2
www.frontiersin.org

FIGURE 2. Estimation of the third state variable in the autonomous mode.

Moreover, considering the effect of the fixed uncompensated part with prescribed D-stability region related to A̲m, Ām, the proposed approach leads to a structure that has the properties of a stable system. By using the tuning parameter δ, the D-stability region can be analytically continued.

Example 2. To demonstrate the application validity of the suggested interval observer, the second example is presented on the linearized dynamic model of a U.S. Navy F-404 engine which powers the F/A-18 aircraft (Kwon et al., 1999). The corresponding dynamic model is a stable interval ostensible purely Metzler system (1), (2), written as

A̲=1.460002.42800.83572.40.37880.310702.1300,Ā=1.460002.42800.33571.40.37880.310702.1300,B=0.41825.20300.39010.12450.51860.0236.

Since the interval matrices of the system are purely Metzler, due to the structure of their second column, it is advantageous if the measurement system corresponds the following matrix C

C=100001,D=111.

Applying analogously as the aforementioned conditions of Theorem 2, the resulting matrix representations are

A̲m=0.07880.00500.00500.84070.07880.38380.00500.00500.0788,ρA̲m=0.15960.07380.0030,
Ām=0.08770.00500.00500.34070.08770.38380.00500.00500.0877,ρĀm=0.15040.08270.0300,
A̲p=1.38120.00502.43300.00502.32120.00500.31570.00502.1512,ρA̲p=0.80892.32132.7234,
Āp=1.37230.00502.43300.00501.31230.00500.31570.00502.1423,ρĀp=0.79991.31242.7145,

where δ̲=δ̄=0.003, λ̲0=0.0758, λ̄0=0.0577, p̲=λ̲0+δ̲=0.0788, and p̄=λ̄0+δ̄=0.0607.

Analogously constructing the diagonal matrices for the interval observer synthesis from the rhombic representations of the interval system matrices and for the used matrix C, the feasible matrix variables resulting from the conditions defined by Theorem 2 are

P=diag2.26803.34532.40980,R1=diag1.44540.00620.24740,R2=diag2.88260.00601.15820,J=0.63731.27100.00180.00180.10270.4806,μ̲=4.2335,μ̄=4.4557.

This infuses the strictly Metzler and Hurwitz matrices A̲pe=A̲pJC and Āpe=ĀpJC as

A̲pe=2.01850.00501.16200.00322.32120.00320.21300.00502.6318,Āpe=2.00960.00501.16200.00321.31230.00320.21300.00502.6229,
ρA̲pe=1.74062.32132.9096,ρĀpe=1.31221.73192.9007,

where A̲peĀpe. In addition, it can be seen that, due to the structure of matrix C, the elements on the second columns of matrices A̲pe and Āpe have not changed compared to A̲p and Āp.

Using these ostensible purely Metzler matrices results in stable, purely Metzler structures

A̲e=A̲JC=2.097301.15700.83752.40.38060.208002.6106,Āe=ĀJC=2.097301.15700.33751.40.38060.208002.6106.

By simulating the response of the observer in the forced mode, it is set as A̲AĀ

A=1.460002.42800.58571.90.37880.310702.1300,ut=Wwt,W=2001,wt=0.3520.076,σd2=0.012,
q0=0.2503.7500.025T,q̲e0=0.1500T,q̄e0=0.400.1T,q̲e0q̄e0.

Figure 3 depicts the time responses of the first system state variable and its upper and lower estimations; Figure 4 shows the time responses for the third system state variable. Although the given system is not positive, it can be seen from Figure 3 and Figure 4 that the behaviors of these state variables are correctly intervally estimated by the proposed interval observer.

FIGURE 3
www.frontiersin.org

FIGURE 3. Estimation of the first state variable in the forced mode.

FIGURE 4
www.frontiersin.org

FIGURE 4. Estimation of the third state variable in the forced mode.

Note for both examples, since μ̲, μ̄ are scalar variables defined directly by a feasible solution of LMIs, they can be indicated as the values at the disturbance attenuation levels. The scalar variable p provides an additional degree of freedom in solving the problem of the dynamics of an interval observer, which should generally be faster than the dynamics of the system. Because the synthesis method is a two-step procedure, it is possible to sequentially define the locations of the stable regions first for the uncontrolled stable dynamics of A̲m and Ām by defining the D-region of stability using the parameter p > 0 (ν > 0 is just some small positive value by means of which Σ is regularized) and then, indirectly via LMIs, finding a solution that guarantees the required rate of estimation error convergence. Both tasks are parametrically dependent, while mutual interaction in the resulting dynamics of the observer is defined by the parameter p > 0, and its interactive setting is, as a rule, sufficient.

6 Concluding remarks

This paper presents new results concerning the interval state estimation of intervally defined ostensible Metzler systems. It proposes how this problem can be formulated using a positive parametric representation and how a constructive procedure based on LMIs can be used respecting the diagonally stabilization principle. It is therefore proven that the gain matrix of the interval observer can be constructed for strict positivity when the stability of the interval observer is defined for a strictly Metzler approximation of the ostensible Metzler system matrix in combination with its stable complement, having a prescribed region of D-stability. The intention was to define the synthesis conditions based only on the quadratic Lyapunov function and to suppress the influence of disturbance in the state estimation by setting the upper bounds of the H norm of its transfer function matrix. The proposed synthesis conditions are not singular, ensuring fast enough convergence of estimation errors, and do not require prior knowledge of the disturbance boundary. With a constant output matrix and the fact that only the upper and lower bounds of the system dynamics matrix are required, such interval observers have relatively high robustness to changes in system parameters. No comparable results in the field of interval estimators for systems with Metzler dynamics seem to have been published so far.

The use of the class of application models was strictly limited by the occurrence of the description of dynamics in the form of Metzler matrices, a class which also includes models of turbo engines applied in the field of networked aircraft fault tolerant control and diagnosis (Jin and Chen, 2014; Li et al., 2020). The goal of the idea was to derive a method for application in the context of interval observer-based methodology for aircraft engine diagnosis and fault-tolerant control (Lamouchi et al., 2022). It is still left as an open question.

This approach requires further theoretical investigation, especially if the considered continuous-time systems have ostensible Metzler system matrices that have a dominant number of negative and zero elements outside the main diagonal. Further research is thus envisaged on both theoretical and applied aspects in anti-disturbance tracking control for unmanned aerial vehicles and drones considering ostensible Metzler and Hurwitz model parameter setting (Yong, 2022; Song et al., 2023).

Data availability statement

The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.

Author contributions

AF elaborated the principles of the observer parameter synthesis and implemented their numerical validation. DK addressed the design and constraint principle assembling into a set of LMIs in design for ostensible Metzler continuous-time linear MIMO systems.

Funding

The research covering the work field presented in this paper was funded by VEGA, the Grant Agency of the Ministry of Education and Academy of Science of the Slovak Republic, under Grant No. 1/0483/21. This support is very gratefully acknowledged.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Ait Rami, M., and Tadeo, F. (2006). “Linear programming approach to impose positiveness in closed-loop and estimated states,” in Proceedings of the 17th international symposium on mathematical theory of networks and systems MTNS 2006, 2470–2477.

Google Scholar

Anderson, J., and Murray, R. (2018). “Structured feedback optimization for Metzler dynamics,” in Proceedings of the 57th IEEE conference on decision and control CDC 2018, 4417–4424.

CrossRef Full Text | Google Scholar

Berman, A., Neumann, M., and Stern, R. (1989). Nonnegative matrices in dynamic systems. New York: John Wiley and Sons.

Google Scholar

Chambon, E., Apkarian, P., and Burlion, L. (2015). “Metzler matrix transform determination using a nonsmooth optimization technique with an application to interval observers,” in Proceedings of the SIAM conference on control and its applications CT 15, 205–211.

CrossRef Full Text | Google Scholar

Efimov, D., Raïssi, T., Chebotarev, S., and Zolghadri, A. (2013). Interval state observer for nonlinear time varying systems. Automatica 49 (1), 200–205. doi:10.1016/j.automatica.2012.07.004

CrossRef Full Text | Google Scholar

Farina, L., and Rinaldi, S. (2000). Positive linear systems. Theory and applications. New York: John Wiley and Sons.

Google Scholar

Ganesan, K. (2007). On some properties of interval matrices. Int. J. Math. Comput. Sci. 1 (1), 35–42.

Google Scholar

Gao, H., Lam, J., Wang, C., and Xu, S. (2005). Control for stability and positivity. Equivalent conditions and computation. IEEE Trans. Circuits Syst. II. Express Briefs 52 (9), 540–544. doi:10.1109/tcsii.2005.850525

CrossRef Full Text | Google Scholar

Guo, Z., Henry, D., Guo, J., Wang, Z., Cieslak, J., and Chang, J. (2020). Metzler matrix-based switching control scheme for linear systems with prescribed performance guarantees. IFAC-PapersOnLine 53 (2), 6428–6433. doi:10.1016/j.ifacol.2020.12.1784

CrossRef Full Text | Google Scholar

Horn, R. A., and Johnson, C. R. (1995). Matrix analysis. New York: Cambridge University Press.

Google Scholar

Huong, D. (2022). Secure interval estimations for time-varying delay interconnected systems using novel distributed functional observers. Int. J. Adapt. Control Signal Process. 36, 1373–1393. doi:10.1002/acs.3400

CrossRef Full Text | Google Scholar

Ito, H., and Dinh, T. D. (2020). Asymptotic and tracking guarantees in interval observer design for systems with unmeasured polytopic nonlinearities. IFAC-PapersOnLine 53 (2), 5010–5015. doi:10.1016/j.ifacol.2020.12.1099

CrossRef Full Text | Google Scholar

Jaulin, L., Kieffer, M., Didrit, O., and Walter, E. (2001). “Applied interval analysis with examples in parameter and state estimation,” in Robust control and robotics (London: Springer-Verlag).

Google Scholar

Jin, Y., and Chen, W. (2014). “Robust fault detection and estimation for turbofan engines subject to adaptive controllers via observer and ToMFIR techniques,” in Proceedings of the 9th IEEE conference on industrial electronics and applications, 672–677.

Google Scholar

Khan, A., Xie, W., Zhang, L., and Liu, L. W. (2020). Design and applications of interval observers for uncertain dynamical systems. IET Circuits, Devices Syst. 14, 721–740. doi:10.1049/iet-cds.2020.0004

CrossRef Full Text | Google Scholar

Krokavec, D., and Filasová, A. (2020a). “Control design for linear strictly Metzlerian descriptor systems,” in Proceedings of the 18th European control conference ECC 2020, 2092–2097.

CrossRef Full Text | Google Scholar

Krokavec, D., and Filasová, A. (2020b). “Interval observer design for uncertain linear continuous-time Metzlerian systems,” in Proceedings of the 28th mediterranean conference on control and automation MED 2020, 15–18.

CrossRef Full Text | Google Scholar

Krokavec, D., and Filasová, A. (2018). “LMI based principles in strictly Metzlerian systems control design,” in Mathematical problems in engineering 2018, 1–14. doi:10.1155/2018/9590253

CrossRef Full Text | Google Scholar

Krokavec, D., and Filasová, A. (2022). “State control of linear systems with potentially Metzler dynamics,” in CONTROLO 2022, lecture notes in electrical engineering. Editors L. B. Palma, R. Neves-Silva, and L. Gomes (Cham: Springer Nature), Vol. 930, 689–701.

CrossRef Full Text | Google Scholar

Kwon, W. H., Kim, P. S., and Park, P. G. (1999). A receding horizon Kalman FIR filter for linear continuous-time systems. IEEE Trans. Automatic Control 44 (11), 2115–2120. doi:10.1109/9.802927

CrossRef Full Text | Google Scholar

Lamouchi, R., Raissi, T., Amairi, M., and Aoun, M. (2022). Interval observer-based methodology for passive fault tolerant control of linear parameter-varying systems. Trans. Inst. Meas. Control 44 (5), 986–999. doi:10.1177/01423312211040370

CrossRef Full Text | Google Scholar

Li, D., Chang, J., and Chen, W. (2022). Event-triggered controller design for LTI systems. A distributed interval observer-based approach. ISA Trans. 131, 146–159. doi:10.1016/j.isatra.2022.04.049

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, T., Tang, X., Ge, J., and Fei, S. (2020). Event-based fault-tolerant control for networked control systems applied to aircraft engine system. Inf. Sci. 512, 1063–1077. doi:10.1016/j.ins.2019.10.039

CrossRef Full Text | Google Scholar

Liu, L. J., Zhao, X., Sun, X. M., and Wang, W. (2017). “New approaches to positive observer design of linear positive systems,” in Proceedings of the 2017 Chinese automation congress CAC 2017, 7195–7198.

CrossRef Full Text | Google Scholar

Liu, X., Xia, Y., and Xi, H. (2011). “Estimating adverse selection and moral hazard effects with hospital invoices data in a government-controlled healthcare system,” in Proceedings of the international conference on electronics, communications and control ICECC 21, 883–901. doi:10.1002/hec.1756

PubMed Abstract | CrossRef Full Text | Google Scholar

Luenberger, D. G. (1971). An introduction to observers. IEEE Trans. Automatic Control 16 (6), 596–602. doi:10.1109/tac.1971.1099826

CrossRef Full Text | Google Scholar

Mason, O. (2012). Diagonal Riccati stability and positive time-delay systems. Syst. Control Lett. 61 (1), 6–10. doi:10.1016/j.sysconle.2011.09.022

CrossRef Full Text | Google Scholar

Mazenc, F., and Bernard, O. (2011). Interval observers for linear time-invariant systems with disturbances. Automatica 47 (1), 140–147. doi:10.1016/j.automatica.2010.10.019

CrossRef Full Text | Google Scholar

Mazenc, F., and Bernard, O. (2014). ISS interval observers for nonlinear systems transformed into triangular systems. Int. J. Robust Nonlinear Control 24, 1241–1261. doi:10.1002/rnc.2937

CrossRef Full Text | Google Scholar

Moisan, M., Bernard, O., and Gouzé, J. L. (2009). Near optimal interval observers bundle for uncertain bioreactors. Automatica 45 (1), 291–295. doi:10.1016/j.automatica.2008.07.006

CrossRef Full Text | Google Scholar

Nikaido, H. (1968). Convex structures and economic theory. New York: Academic Press.

Google Scholar

Peaucelle, D., Henrion, D., Labit, Y., and Taitz, K. (2002). User’s Guide for SeDuMi interface 1.04 (toulouse. LAAS-CNRS.

Google Scholar

Raïssi, T., and Efimov, D. (2018). Some recent results on the design and implementation of interval observers for uncertain systems. Automatisierungstechnik 66 (3), 213–224. doi:10.1515/auto-2017-0081

CrossRef Full Text | Google Scholar

Shores, T. S. (2007). Applied linear algebra and matrix analysis. New York: Springer.

Google Scholar

Shorten, R., Mason, O., and King, C. (2009). An alternative proof of the Barker, Berman, Plemmons (BBP) result on diagonal stability and extensions. Linear Algebra its Appl. 430 (1), 34–40. doi:10.1016/j.laa.2008.06.037

CrossRef Full Text | Google Scholar

Shu, Z., Lam, J., Gao, H., Du, B., and Wu, L. (2008). Positive observers and dynamic output-feedback controllers for interval positive linear systems. IEEE Trans. Circuits Syst. I. Regul. Pap. 55 (10), 3209–3222. doi:10.1109/tcsi.2008.924116

CrossRef Full Text | Google Scholar

Smith, H. L. (1995). Monotone dynamical systems. An introduction to the theory of competitive and cooperative systems. Providence: American Mathematical Society.

Google Scholar

Son, N. K., and Hinrichsen, D. (1996). Robust stability of positive continuous time systems. Numer. Funct. Analysis Optim. 17 (5-6), 649–659. doi:10.1080/01630569608816716

CrossRef Full Text | Google Scholar

Song, Y., Yong, K., and Wang, X. (2023). Disturbance interval observer-based robust constrained control for unmanned aerial vehicle path following. Drones 7, 90. doi:10.3390/drones7020090

CrossRef Full Text | Google Scholar

Tanaka, T., and Langbort, C. (2011). The bounded real lemma for internally positive systems and H-infinity structured static state feedback. IEEE Trans. Automatic Control 56 (9), 2218–2223. doi:10.1109/tac.2011.2157394

CrossRef Full Text | Google Scholar

Wang T, T., Li, Y., and Xiang, W. (2022). Design of interval observer for continuous linear large-scale systems with disturbance attenuation. J. Frankl. Inst. 359 (8), 3910–3929. doi:10.1016/j.jfranklin.2022.03.014

CrossRef Full Text | Google Scholar

Wang, X., Jiang, G. P., Yang, W., Su, H., and Wang, X. (2020). Neighborhood interval observer based coordination control for multi-agent systems with disturbances. IFAC-PapersOnLine 53 (2), 10994–10999. doi:10.1016/j.ifacol.2020.12.023

CrossRef Full Text | Google Scholar

Wang, X. L., Luo, H., Zhao, J., and Yin, Q. (2022). Consensus of time-varying interval uncertain multiagent systems via reduced-order neighborhood interval observer. Complexity 2022, 1–14. doi:10.1155/2022/5800097

CrossRef Full Text | Google Scholar

Yong, K., Chen, M., Shi, Y., and Wu, Q. (2020). Hybrid estimation strategy-based anti-disturbance control for nonlinear systems. IEEE Trans. Automatic Control 66 (10), 4910–4917. doi:10.1109/tac.2020.3043998

CrossRef Full Text | Google Scholar

Yong, K. (2022). Disturbance interval observer-based carrier landing control of unmanned aerial vehicles using prescribed performance. Sci. Sin. Inf. 52 (9), 1711–1726.

Google Scholar

Zhang, H., Huang, J., and He, S. (2022). Fractional-order interval observer for multiagent nonlinear systems. Fractal Fract. 6 (7), 355. doi:10.3390/fractalfract6070355

CrossRef Full Text | Google Scholar

Keywords: Metzler systems, parametric constraints, diagonal stabilization, linear matrix inequalities, applied interval analysis, interval observers

Citation: Krokavec D and Filasová A (2023) Interval observers design for systems with ostensible Metzler system matrices. Front. Aerosp. Eng. 2:1158718. doi: 10.3389/fpace.2023.1158718

Received: 04 February 2023; Accepted: 17 April 2023;
Published: 10 May 2023.

Edited by:

Robert Fonod, Swiss Federal Institute of Technology Lausanne, Switzerland

Reviewed by:

Shuyi Shao, Nanjing University of Aeronautics and Astronautics, China
Xiaodong Shao, Beihang University, China

Copyright © 2023 Krokavec and Filasová. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Dušan Krokavec, dusan.krokavec@tuke.sk

Download