# Full and Reduced-Order Unknown Input Observer Design for Linear Time-Delay Systems with Multiple Delays

Authors: Seifeddine Ben Warrad; Olfa Boubaker; Mihai Lungu; Saleh Mobayen
Journal: Mathematical Problems in Engineering
 (2018)

---

## Abstract

In this work, we consider the design problem of a full-order and also reduced-order unknown input observers for a particular class of time-delay systems. Asymptotic stability and existence conditions for the designed observers are established. The quadruple-tank process is used as a benchmark to prove the efficiency of the proposed algorithms.


---

## Body

## 1. Introduction

Time delays emerge in many industrial applications such as hydraulic systems, automotive applications, and telecommunication networks [1–7]. Commonly, delays describe propagation phenomena, energy transfer, or data transmission. Furthermore, it is recognized that delays are the key sources of instability and poor dynamics [8]. Since the sixteen’s century, the research theme of time-delay systems has concerned lots of interest and has been considered as one of the essential research areas in control theory for which many significant research works have been devoted. The topical survey [9] and the recent books [10–12] as well as their related references prove the richness and maturity of this field. Time delays can appear in either the state variables, control inputs, or measurement outputs and the negligence of time delays in the controller design can lead to undesirable dynamics such as oscillations, bifurcations, chaos dynamics, and instability [8]. In this trend, stabilization of linear systems with delays has obtained a huge interest (see, for example, recent works [13–15] and related references). As reported in [16], for stabilization of continuous linear systems subject to time delays, most research works only focus in systems with input delays and solve the problem in the absence of state delays. To the best of the authors’ knowledge, few research works have considered the stabilization problem with multiple delays in states and inputs where the problem is generally solved by assuming fully available state variables [1] or by resorting to output-feedback controllers [17]. The design problem of observer-based stabilization for systems with multiple delays in states and inputs is rarely solved [16, 18]. Until now, this issue of theoretical and practical importance has been remained a challenging problem for systems subject to unknown inputs, uncertainties, and input saturation.In control theory, modeling and controlling benchmark systems is another fundamental issue. Such systems are considered as challenging control problems, yet, in spite of their uncomplicated arrangement. In this outline, several systems are commonly used in control theory for testing novel control algorithms. We can cite in this framework the three-tank process [19], mass-spring-damper system [20], bouncing ball [21], TORA system [22], hard-disk drive system [23], magnetic levitation system [24], cart-inverted pendulum [25], Furuta pendulum [26], reaction wheel pendulum [27], beam-and-ball [28], two-link flexible manipulator [29], and so on. In the same direction, the well-known quadruple-tank benchmark [30–36] has attracted many attention since it can show elementary notions in estimation and control theory, particularly performance limitation due to the nonminimum-phase zeros and their output directions for multivariable systems. This system is composed of four coupled tanks, two pumps, and two valves [35]. Many research papers have been devoted to the problem of observer and controller design of this benchmark system without consideration of time delays [37–43] and to the best of our information, no researches have been investigated for unknown input observer design of time-delay quadruple-tank systems.Motivated by the authors’ early works for linear systems without delays [44], in this paper, we expand the problem of full-order unknown input observer (FOUIO) andreduced-order unknown input observer(ROUIO) design for linear systems with multiple time delays in states and inputs affected by unknown inputs. In this framework, it is important to note that a FOUIO is designed in [45–47] for linear systems with a single delay in the states and two delays in the inputs whereas another FOUIO is proposed in [48] for linear systems with a single delay in the states and a single delay in the inputs. For the last works, only the state information is used to design observers in open loops and to the best of our knowledge, there are no related works for which both the state vector and the input vectors are used to design the observer. Furthermore, in our best knowledge, there is not in the vast literature ROIU design approach that attempts to solve such a problem.The design procedures of both reduced-order and full-order observers are considered here. The key contributions of this work can be recapitulated as follows:(i)
The designed observers’ approaches do not require approximation of the infinite dimensional model of the delayed system in a lumped one as described in our previous works [45, 46].(ii)
Different from [46] where only the state vector is used to construct observers in open loop, in this paper, we construct the novel observers by using both state and input vectors. The structure of the proposed observers allows us to introduce more observer gains that can accomplish better control performance.(iii)
Different from the observer designed in [47] where the delayed system introduces a single delay in the states and a single delay in the inputs, in this paper, multiple time delays in the inputs are considered.The rest of this paper is organized as follows: the problem statement is considered in Section2. The full-order observer is designed in Section 3 whereas the reduced-order observer is developed in Section 4. Section 5 is devoted to the application of the proposed approaches to the quadruple-tank benchmark. Finally, conclusions are provided in Section 6.


## 2. Problem Statement

Consider a class of time-delay systems described by(1)x˙t=A0xt+A1xt-τ0+∑i=0kBiut-τi+1+Ddt,yt=Cxt,where xt∈Rn is the state vector, ut∈Rm  is the input vector, dt∈Rq  is the unknown input vector, and yt∈Rp  is the output vector; A0∈Rn×n,A1∈Rn×n, Bi∈Rn×m,C∈Rp×n, and D∈Rn×q are constant matrices with appropriate dimensions, while τ0andτi+1 for i=0, …, k are known and constant time delays.The objective of this paper is to design two different approaches of unknown input observers (UIOs) for system (1):( 1 ) An FOUIO is described by(2)ζ˙t=Nζt+Lyt+Jyt-τ0+∑i=0kGiut-τi+1,x^t=ζt-Eyt,where ζt∈Rn and x^t∈Rn are the observer state vector and the estimated state vector, respectively. N,L,J,E, and Gi (i=0…k) are matrices of appropriate dimensions to be calculated such that x^t converges asymptotically to xt under the following assumptions: (A1) rankD=q; (A2) rankC=p; (A3) the pair A0,C is observable (detectable).( 2 ) AnROUIOis also designed in this paper for the estimation of the state vector xt by obtaining the estimation of the vector: ζ-t=Fxt,F∈Rs×n, where ζ-t∈Rs is the vector enclosing a part of the system’s states or a combination of them. The dynamics of the ROUIO is(3)Zζ-^˙t=N-ζ-^t+L-yt+J-yt-τ0+∑i=0kG-iut-τi+1,x^t=Hζ-^t+E-yt,where ζ-^t∈Rs and x^t∈Rn are the observer state vector and estimated state vector, respectively. N-∈Rs×s,L-∈Rs×p,J-∈Rs×p,G-i∈Rs×m for (i=0...k),H∈Rn×s,E-∈Rn×p,F∈Rs×n and Z∈Rs×s are matrices of appropriate dimensions to be calculated such that ζ-^t and x^t converge asymptotically to ζ-t and xt, respectively, under the assumptions as(A4)rankD=q; (A5) rankC=p; (A6) the pair A0,C is observable; (A7) n>s>p≥q; (A8) n=p+q.Remark 1.
For the design procedure of ROUIO, a particular form for the matrix C can be chosen asC=C10p×q, where C1∈Rp×p is a full rank (nonsingular) matrix. This is not restrictive as long as the matrix C is full row rank (conditions A2≡A5), and there will always be an orthogonal transformation which leads to the equation yt=C10xt.
Either for theFOUIO (2) or theROUIO (3), system (1) is controlled using the state-feedback controller as(4)u=Kx^twhere K∈Rm×n is designed using one of the appropriate approaches for system (1), for example, the approach proposed in [1].


## 3. Design of the Full-Order Observer

Theorem 2.
Consider system (1) with assumptions (A1-A3) presented above. The FOUIO (2) is convergent, i.e., et→0 as t→∞ with an arbitrary convergence rate for any x(0), d(t), and u(t), if there exists a positive-definite matrix R∈Rn×n verifying the matrix inequality(5)NTR+RN<0and if the matrix P=I+EC∈Rn×n fulfills the following conditions:(6)NP+LC-PA0=0,JC-PA1=0Gi-PBi=0∀i=0⋯k,PD=0.Proof.
Define the observer reconstruction error vector as the difference between the estimated statex^t described by (2) and the state vector related to system (1) as(7)et=x^t-xt=ζt-xt-Eyt=ζt-I+ECxt.
Using (1) and (2), the dynamics of the estimation error becomes(8)e˙t=ζ˙t-I+ECx˙t,or(9)e˙t=Net+NI+EC+LC-I+ECA0xt+JC-I+ECA1xt-τ0+∑i=0kGi-I+ECBiut-τi+1-I+ECDdt.
UsingP=I+EC, the error dynamics (9) is written in the following form:(10)e˙t=Net+NP+LC-PA0xt+JC-PA1xt-τ0+∑i=0kGi-PBiut-τi+1-PDdt.
If conditions (6) are satisfied, the error dynamics of the observer is expressed as e˙t=Net. FOUIO’s convergence is then achieved if the LMI (5) is verified.The determination of the matrices N, L, J, and P from the first, second, and fourth equations of (6) is a difficult task because we have to calculate four matrices by using only three equations. In order to use the well-known results obtained for the classical full-order observer without unknown inputs [2], by means of the notation(11)K0=L+NE,the first term of (6) becomes(12)N=PA0-K0C.Now, using (11) and (12), the dynamics (2) can be written as(13)ζ˙t=PA0-K0Cζt+Lyt+Jyt-τ0+∑i=0kGiut-τi+1,x^t=ζt-Eyt.If the pair(PA0,C) is not observable, the calculation of matrix K0 is made such that the observer is asymptotically stable if and only if (PA0,C) is detectable [2]. By extending the approach from [2], we can conclude that the necessary and sufficient conditions to design a stable observer are given by the following theorem.Theorem 3 (see [2]).
For system (1), the full-order observer (13) exists if and only if (1)  rankCD=q; (2)  ranksI-PA0C=n,∀s∈Ł, Res≥0.Proof.
The proof of the theorem is an extension of the approach from [2]. Because system (13) is in the form of a standard observer equation, then the matrix K0 can be calculated such that the observer (2) is asymptotically stable if and only if the pair PA0,C is observable. In [2], it is proved that the above condition (2) is equivalent with the pair (PA0,C)- observable or at least detectable (first existence condition of the full-order observer). The second constraint is the condition (1) from Theorem 3; bearing in mind the assumptions A1 and A2, as well as the dimensions of the matrices C and D, this condition is always fulfilled. Therefore, up to this point in the design of the unknown input observer, the only existence condition of FOUIO is related to the observability or at least detectability of the pair (PA0,C).To calculate N, L, J, P, G0, and G1, we partition matrices A0, A1, C, N, P, L, and J as follows:(14)A0=A011A012A021A022,A1=A111A112A121A122,CT=C1TC2T,N=N11N12N21N22,P=P11P12P21P22,L=L1L2,J=J1J2,where(15)A0∈Rn×n,A011∈Rn-p×n-p,A012∈Rn-p×p,A021∈Rp×n-p,A022∈Rp×p,A1∈Rn×n,A111∈Rn-p×n-p,A112∈Rn-p×p,A121∈Rp×n-p,A122∈Rp×p,N∈Rn×n,N11∈Rn-p×n-p,N12∈Rn-p×p,N21∈Rp×n-p,N22∈Rp×p,P∈Rn×n,P11∈Rn-p×n-p,P12∈Rn-p×p,P21∈Rp×n-p,P22∈Rp×p,L∈Rn×p,L1∈Rn-p×p,L2∈Rp×p,J∈Rn×p,J1∈Rn-p×p,J2∈Rp×p,C∈Rp×n,C1∈Rp×n-p,C2∈Rp×p.The observer’s design is concentrated into the following theorem.Theorem 4.
If the rank conditions in Theorem3 are satisfied, the UIO (2) can be designed for system (1) as follows:(16)E=DCD+C,P=In-DCD+C,Gi=PBifork=0⋯k,J=PA1C+,N=N11F12-F11-F12P22-1P21C1-C2P22-1P21+C2P22-1N21F22-F21-F22P22-1P21C1-C2P22-1P21+C2P22-1,L=F11-F12P22-1P21C1-C2P22-1P21+F21-F22P22-1P21C1-C2P22-1P21+,where(17)F11=P11A011+P12A021-N11P11,F12=P11A012+P12A022-N11P12,F21=P21A011+P22A021-N21P11,F22=P21A012+P22A022-N21P21,while N11 and N21 are arbitrarily chosen; P22 should be a nonsingular matrix.Proof.
SincerankCD=q, i.e., the matrix (CD) has full column rank, we can use the pseudoinverse of CD, given by [17]: CD+=CDTCD-1CDT. By means of this pseudoinverse matrix, we can solve the fourth equation of (6); one successively gets(18)I+ECD=0⇔D+ECD=0⇔E=-DCD+.
Hence, one getsP=In+EC=In-DCD+C. Replacing P into the third equation of (6), we obtain Gi=PBi=In-DCD+CBi. Now, substituting (14) into the first equation of (6), we obtain(19)N11P11+N12P21+L1C1-P11A011-P12A021=0n-p×n-p,N11P12+N12P22+L1C2-P11A012-P12A022=0n-p×p,N21P11+N22P21+L2C1-P21A011-P22A021=0p×n-p,N21P12+N22P22+L2C2-P21A012-P22A022=0p×p
By using (17), system (19) is divided into the following systems:(20)N12P21+L1C1=F11,N12P22+L1C2=F12,(21)N22P22+L2C2=F22,N22P21+L2C1=F21.
Working in the hypothesis: the matrixC1-C2P22-1P21∈Rp×n-p is full column rank. By solving systems (20) and (21), we obtain(22)L1=F11-F12P22-1P21C1-C2P22-1P21+,L2=F21-F22P22-1P21C1-C2P22-1P21+,N12=F12-F11-F12P22-1P21C1-C2P22-1P21+C2P22-1,N22=F22-F21-F22P22-1P21C1-C2P22-1P21+C2P22-1.
Now, solving the second equation of (6), we get(23)J=PA1C+.
The unknown matrices of the FOUIO (2), i.e., N, L, J, P, and Gi, have been successfully determined with respect to the arbitrarily chosen matrices N11 and N21. According to Theorem 2, the obtained matrix N should satisfy two constraints: (1) there exists a symmetric and positive-definite matrix R∈Rn×n which verifies the LMI (5), i.e., the matrix N is Hurwitz; (2) the pair PA0,C is observable. If these two conditions are satisfied, the matrix N has been obtained properly, the dynamics of the observer error has an homogeneous form, and et→0 as t→∞ with an arbitrary convergence rate for any x(0), d(t), and u(t). Otherwise, other matrices N11, N21 are chosen, the matrices from (17) are again calculated, and systems (20) and (21) are solved until the above presented constraints are fulfilled. In fact, no condition regarding the observability or detectability of the pairs PA0,C is required in our design approach. Therefore, the unique existence condition of the observer is rankCD=q. The theorem’s proof is now complete.The algorithm associated to the delay-dependent observer is summarized below.Algorithm 5.
Design of FOUIOStep 1.
Given system (1), verify the related assumptions A1-A3.Step 2.
If the existence conditions given in Theorem3 hold, one calculates the generalized inverse matrix CD+ related to the matrix CD. By using it, one computes the matrix E by means of the first equation of (16); then, one computes the matrices P, G0, and G1 with the second, third, and fourth equations of (16), respectively.Step 3.
One partitions the matrices A0, A1, C, N, P, L, and J according to (14).Step 4.
Matrices N11 and N21 are arbitrarily chosen; the matrices from (17) are determined.Step 5.
One solves systems (20) and (21); after that, the matrices N, L, and J are built.Step 6.
Using the obtained matrix N, if the required constraints are satisfied, then(1) the matrix N is Hurwitz, i.e., there exists a symmetric and positive-definite matrix R∈Rn×n verifying the LMI (5); (2) the pair PA0,C is observable or at least detectable; the matrix N has been obtained properly; otherwise, Steps 4–6 are repeated (in a “while” loop) until these conditions are fulfilled. Bearing in mind that there is no risk for infinite “while” loop, there can be concluded that the FOUIO’s design algorithm has no existence conditions, being characterized by lack of a priori restrictions on the class of systems to be considered.Step 7.
The observer described by (2) is completely designed and the time history of the system’s estimated states can be obtained.


## 4. Design of the Reduced-Order Observer

Theorem 6.
Consider the LTI multivariable system (1) with multiple delays under the five assumptions (A4-A8) presented above; the ROUIO (3) is asymptotically stable, if and only if there exist symmetric and positive-definite matrix R-∈Rs×s verifying the LMI:(24)N~TR-+R-N~<0,with N~=Z-1N-, and if there exists a vector ζ-t=Fxt such that the next conditions are satisfied:(25)N-F+L-C-ZFA0=0,J-C-ZFA1=0,G-i-ZFBi=0∀i=0⋯k,ZFD=0,E-C+HF=In;Proof.
Define the observer error vector defined as the difference between the estimated vectorζ-^t described by (3) and the state vector related to system (1) as(26)e-t=ζ-^t-ζ-t
Using (1), (3), and (4), the dynamics of the estimation error becomes(27)Ze-˙t=N-e-t+N-F-ZFA0+L-Cxt+J-C-ZFA1xt-τ0+∑i=0kG-i-ZFBiut-τi+1-ZFDdt.
IfN-F+L-C-ZFA0=0,J-C-ZFA1=0,G-i-ZFBi=0∀i=0⋯k,ZFD=0 the error dynamics (27) can be written as(28)e-˙t=N~e-t,where N~=Z-1N-.
Let the state estimation error be defined asext=x^t-xt. From (1) and (3), one obtains ext=HF-In+E-Cxt. If HF-In+E-C=0, ext=0 for any x(t) results.
Furthermore, if condition (24) is satisfied then the observer’s error dynamics (28) is asymptotically stable. Now, the proof of the Theorem 6 is complete.To compute the matricesF,Z,E-,N-,L-,J-,H,G-0, and G-1, let us partition the matrices A0, A1, C, N-,H,L-,J-,F,Z, and D as follows:(29)A0=A011A012A021A022,A1=A111A112A121A122,C=C10p×q=C~1C~2,D=D1D2,N-=N-11N-11N-21N-22=N-1N-2,H=H11H11H21H22,L-=L-1L-2,J-=J-1J-2,E-=E-1E-2,F=F11F12F21F22=0p×qIp×p0s-p×q0s-p×p,Z=Z11Z12Z21Z22=Z1Z2,Z1=Z11Z21,Z2=Z12Z22,where(30)A0∈Rn×n,A011∈Rq×q,A012∈Rq×p,A021∈Rp×q,A022∈Rp×p,A1∈Rn×n,A111∈Rq×q,A112∈Rq×p,A121∈Rp×q,A122∈Rp×p,N-∈Rs×s,N-11∈Rq×p,N-12∈Rq×s-p,N-21∈Rs-q×p,N-22∈Rs-q×s-p,H∈Rn×s,H11∈Rq×p,H12∈Rq×s-p,H21∈Rp×p,H22∈Rp×s-p,L-∈Rs×p,L-1∈Rq×p,L-2∈Rs-q×p,J-∈Rs×p,J-1∈Rq×p,J-2∈Rs-q×p,E-∈Rn×p,E-1∈Rq×p,E-2∈Rp×p,D∈Rn×q,D1∈Rq×q,D2∈Rp×q,Z∈Rs×s,Z11∈Rq×p,Z12∈Rq×s-p,Z21∈Rs-q×p,Z22∈Rs-q×s-p,Z1∈Rs×p,Z2∈Rs×s-p,N-1∈Rs×p,N-2∈Rs×s-p,C∈Rp×n,C1∈Rp×p,C~1∈Rp×q,C~2∈Rp×p.The ROUIO’s design is concentrated into the following theorem.Theorem 7.
Consider the LTI multivariable (1) under the assumptions A4-A8; the ROUIO (3) is asymptotically stable if and only if the following conditions are satisfied: (1) there exists the matrix R->0 verifying the LMI (24); (2) the matrices G-0,G-1,N-,L-,J-,E-, and H are given by(31)Z1=Is×p-D-2C1D2+C1,G-i=0s×qZ1Bi∀i=1⋯k,N-=Z1A014-A013C~+1C~2N-2,L-=Z1A021C~+1,J-=ZFA1C+,E-=C~+1IP-C~1C1C~1+C1,H=-C~1+C~20q×s-pIp-C~2+C~1C1C~1+C1C~20p×s-p,where the matrices Z2∈Rs×s-p and N-2∈Rs×s-p are arbitrarily chosen, C~1∈Rp×q and C~2∈Rp×p are submatrices of the matrix C, and D-2=Is×pD2.Proof.
Using (29) and (30), the fourth equation of (25) becomes(32)Z1D2=0s×q.
To compute the matrixZ1 from this equation, one can choose(33)Z1=Is×p+T1C1,with T1∈Rs×p an unknown matrix to be computed. With the notation D-2=Is×pD2, one gets T1=-D-2C1D2+, where C1D2+ is the generalized pseudoinverse of C1D2, given by C1D2+=C1D2TC1D2-1C1D2T. According to the assumption A5 and to the fact that the matrix C1 is full rank, we can conclude that rankC1D2=q, and the generalized pseudoinverse of C1D2 can be always obtained.
Equation (33) can be then written as(34)Z1=Is×p-D-2C1D2+C1,and the matrix Z2∈Rs×s-p is chosen arbitrarily.
Now, using the expression of the matrix Z, we obtain the matricesG-i in the third equation of (25). We arbitrarily choose the matrices N12 and N22; as a consequence, the matrix N-2 is arbitrarily chosen. Using this as well as the form of matrix C given in (29) and (30), i.e., C=C10p×q=C~1C~2, the first equation (25) can be written in a form of matrix system as follows:(35)-Z11A021+L-1C~1=0q×q,N-11-Z11A022+L-1C~2=0q×p,-Z21A021+L-2C~1=0s-q×q,N-21-Z21A022+L-2C~2=0s-q×p;by solving system (35), we obtain(36)N-11=Z11A022-A021C~1+C~2,L-1=Z11A021C~1+,N-21=Z21A022-A021C~1+C~2,L-2=Z21A021C~1+;bearing in mind the fact that q<p, it results that C~1 is a submatrix of the nonsingular matrix C1; from this, one deduces that rankC~1=q, i.e., C~1 is a full column rank matrix.
Now, solving the second equation of (25), we get(37)J-=ZFA1C+.
Using (29) and (30), the sixth equation (25) can be transformed into a matrix system as(38)E-1C~1=Iq,E-1C~2+H11=0q×p,E-2C~1=0p×q,E-2C~2+H21=Ip;from the first and second equations of (38), we obtain E-1=C~1+ and H11=-E-1C~2=-C~1+C~2. On the other hand, the matrix E-2 is calculated from the third equation of (35) as follows:(39)E-2=Ip+T2C1,where T2∈Rp×p is an unknown matrix to be computed. Because C1∈Rp×p is a nonsingular matrix and rankC~1=q, one gets rankC1C~1=q, i.e., C1C~1 is a full column rank matrix. Taking this into account, one successively gets T2=-C~1C1C~1+ and E-2=Ip-C~1C1C~1+C1. From the fourth equation of system (38), we obtain H12=Ip-C~2+C~1C1C~1+C1C~2; by concatenation of the matrices E-1,E-2,H11,H21,H12=0q×s-p and H22=0p×s-p, the last two equations (31) result. Now, the proof of the Theorem 7 is complete.Algorithm 8.
Design of the ROUIO.Step 1.
Given system (1), verify the five related assumptions A4-A8.Step 2.
The matrixF of from (29) is chosen; the matrices A0,A1,C,D,N-,H,L-,J-,E-, and Z are partitioned with (29) and (30). The matrices Z12 and Z22 and thus the matrix Z2 are arbitrarily chosen; the matrices Z1, G-i∀i=1⋯k, are obtained via (31).Step 3.
The matricesN-12,N-22 and thus the matrix N-2 are arbitrarily chosen; by solving system (35) with respect to the matrices N-11,N-21, L-1,L-2 or, directly, by means of (36), these four matrices are obtained; the matrices N-,J- and L- are then calculated with (31). Matrices E- and H are computed by solving system (38).Step 4.
Using the above calculated matrixN-, one checks if there exists a symmetric and positive-definite matrix R-∈Rn×n verifying the LMI (24); if so, the matrix N- has been properly obtained; otherwise, Steps 2–4 are repeated (in a “while” loop) until this condition is fulfilled.Step 5.
The observer described by (3) is completely designed and the time history of the system’s estimated states can be obtained.


## 5. Numerical Simulation Results

### 5.1. The Quadruple-Tank Process

The quadruple-tank process is shown in Figure1. It is composed of four connected tanks and two pumps that split water into two tanks [30, 31, 49]. The inlet flow of each tank is measured by an electromagnetic flow-meter and regulated by a pneumatic valve whereas the level of each tank hi (I =1…4) is measured by means of a pressure sensor. The process inputs are the input voltages of the pumps (ϑ1 and ϑ2) and the output variables are the tank levels y1 and y2.Figure 1
The quadruple-tank process.


### 5.2. The Quadruple-Tank Process Model

For practical considerations, we assume that the operating regime of the process is well known that the quadruple-tank process is affected by delays and that the transport delays are perfectly symmetric. In such a case, the model of the process can be written as case study of model (1) and can be described by [30](40)x˙t=A0xt+A1xt-τ1+B0ut-τ2+B1ut-τ3+Ddt,yt=Cxt,with x=h1h2h3h4Tbeing the state vector, u=ϑ1ϑ2Tbeing the input vector, and y=y1y2Tbeing the output vector, respectively, and with(41)A0=-a1A1g2h100000-a2A2g2h200000-a3A3g2h300000-a4A4g2h40,B0=γ1k1A100γ2k2A20000,A1=00a3A1g2h300000a4A2g2h4000000000,B1=000001-γ2k2A31-γ1k1A40.Using the following numerical values(42)D1=D2=D3=D4=9.2cm,A1=A2=A3=A4=πD122,hmax=50cm,hi0=hmax2=25cm,g=981cm·s-2,d1=0.2cm,a1=a2=πd12,d3=0.9,a3=a4=πd322,k1=7.39,k2=6.92,the matrices from (1) are given by(43)A0=-0.00210000-0.00210000-0.04240000-0.0424,B0=0.1113×γ1000.1042×γ20000,A1=000.042400000.042400000000,B1=000001-γ2×0.10421-γ1×0.10420.


### 5.3. Software Implementation of the Full-Order Observer

For simulation results, the following numerical values are chosen:x0=-446-5, dt=0.3sin⁡2πft, DT=1020-1215-232830-36, γ1=0.333, γ2=0.307, α=0.4420, and τ1=5s, τ2=3s, τ3=4s. The control gain matrix is computed using the approach proposed in [1] (44)K=-0.1603-0.1765-0.0795-0.2073-0.1977-0.1579-0.2288-0.0772For the FOUIO software implementation, one obtains(45)E=-100-11.264-0.032-1.540.02,P=000000001.26-0.03210-1.54-0.02001,G0=00000.046-0.01-0.0570.06,G1=000000.0720.0740,N=-10000-10011-0.04240110-0.0424,L=0000-0.05-0.013-0.060.08,J=0000-0.0430.095-0.0250.017,


### 5.4. Software Implementation of the Reduced-Order Observer

During the first step of the Algorithm8, we find the matrices and vectors from (3). In this paper, the order is chosen as equal n=4, m=2, p=2, q=2, and s=3 and the constant time delays are chosen as τ1=5s,τ2=2s, and τ3=4s. We obtain(46)E-=10010000,H=000000100010,G-0=000.1000.07,G-1=00.0720.074000,N-=1013002.945000.46000,L-=0000-0.05-0.013-0.060.08,J-=0000-0.0430.095-0.0250.017,Z=10-140.2220-0.04390.52550.13770.01110.8992000.9452Figure2 shows the effectiveness of two designed observers for the quadruple-process and their performances for which the variables x1 and x2 (variation of water levels in tanks 1 and 2) are measured and only the variables x3 and x4 (variation of water levels in tanks 3 and 4) are estimated. The effectiveness of two observers is proved by the superpose of the three curves.Figure 2
FOUIO and ROUIO of the quadruple-tank process.


## 5.1. The Quadruple-Tank Process

The quadruple-tank process is shown in Figure1. It is composed of four connected tanks and two pumps that split water into two tanks [30, 31, 49]. The inlet flow of each tank is measured by an electromagnetic flow-meter and regulated by a pneumatic valve whereas the level of each tank hi (I =1…4) is measured by means of a pressure sensor. The process inputs are the input voltages of the pumps (ϑ1 and ϑ2) and the output variables are the tank levels y1 and y2.Figure 1
The quadruple-tank process.


## 5.2. The Quadruple-Tank Process Model

For practical considerations, we assume that the operating regime of the process is well known that the quadruple-tank process is affected by delays and that the transport delays are perfectly symmetric. In such a case, the model of the process can be written as case study of model (1) and can be described by [30](40)x˙t=A0xt+A1xt-τ1+B0ut-τ2+B1ut-τ3+Ddt,yt=Cxt,with x=h1h2h3h4Tbeing the state vector, u=ϑ1ϑ2Tbeing the input vector, and y=y1y2Tbeing the output vector, respectively, and with(41)A0=-a1A1g2h100000-a2A2g2h200000-a3A3g2h300000-a4A4g2h40,B0=γ1k1A100γ2k2A20000,A1=00a3A1g2h300000a4A2g2h4000000000,B1=000001-γ2k2A31-γ1k1A40.Using the following numerical values(42)D1=D2=D3=D4=9.2cm,A1=A2=A3=A4=πD122,hmax=50cm,hi0=hmax2=25cm,g=981cm·s-2,d1=0.2cm,a1=a2=πd12,d3=0.9,a3=a4=πd322,k1=7.39,k2=6.92,the matrices from (1) are given by(43)A0=-0.00210000-0.00210000-0.04240000-0.0424,B0=0.1113×γ1000.1042×γ20000,A1=000.042400000.042400000000,B1=000001-γ2×0.10421-γ1×0.10420.


## 5.3. Software Implementation of the Full-Order Observer

For simulation results, the following numerical values are chosen:x0=-446-5, dt=0.3sin⁡2πft, DT=1020-1215-232830-36, γ1=0.333, γ2=0.307, α=0.4420, and τ1=5s, τ2=3s, τ3=4s. The control gain matrix is computed using the approach proposed in [1] (44)K=-0.1603-0.1765-0.0795-0.2073-0.1977-0.1579-0.2288-0.0772For the FOUIO software implementation, one obtains(45)E=-100-11.264-0.032-1.540.02,P=000000001.26-0.03210-1.54-0.02001,G0=00000.046-0.01-0.0570.06,G1=000000.0720.0740,N=-10000-10011-0.04240110-0.0424,L=0000-0.05-0.013-0.060.08,J=0000-0.0430.095-0.0250.017,


## 5.4. Software Implementation of the Reduced-Order Observer

During the first step of the Algorithm8, we find the matrices and vectors from (3). In this paper, the order is chosen as equal n=4, m=2, p=2, q=2, and s=3 and the constant time delays are chosen as τ1=5s,τ2=2s, and τ3=4s. We obtain(46)E-=10010000,H=000000100010,G-0=000.1000.07,G-1=00.0720.074000,N-=1013002.945000.46000,L-=0000-0.05-0.013-0.060.08,J-=0000-0.0430.095-0.0250.017,Z=10-140.2220-0.04390.52550.13770.01110.8992000.9452Figure2 shows the effectiveness of two designed observers for the quadruple-process and their performances for which the variables x1 and x2 (variation of water levels in tanks 1 and 2) are measured and only the variables x3 and x4 (variation of water levels in tanks 3 and 4) are estimated. The effectiveness of two observers is proved by the superpose of the three curves.Figure 2
FOUIO and ROUIO of the quadruple-tank process.


## 6. Conclusion

This study considers the design of two approaches of UIOs for systems with multiple time-delays: a full-order delay-dependent UIO and a reduced-order delay-dependent UIO. The quadruple-tank process was used as a benchmark to prove the effectiveness of the new observer algorithms for the case study of multivariable nonminimum phase systems with multiple delays.



---

*Source: 1745734-2018-11-26.xml*