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.

1. Introduction

Time delays emerge in many industrial applications such as hydraulic systems, automotive applications, and telecommunication networks [17]. 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 [1012] 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 [1315] 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 [3036] 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 [3743] 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) and reduced-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 [4547] 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 Section 2. 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 bywhere is the state vector,   is the input vector,   is the unknown input vector, and   is the output vector; , , and are constant matrices with appropriate dimensions, while 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):

An FOUIO is described bywhere and are the observer state vector and the estimated state vector, respectively. , and (i=0…k) are matrices of appropriate dimensions to be calculated such that converges asymptotically to under the following assumptions: (A1) ; (A2) ; (A3) the pair is observable (detectable).

An ROUIO is also designed in this paper for the estimation of the state vector by obtaining the estimation of the vector: , where is the vector enclosing a part of the system’s states or a combination of them. The dynamics of the ROUIO iswhere and are the observer state vector and estimated state vector, respectively. for and are matrices of appropriate dimensions to be calculated such that and converge asymptotically to and , respectively, under the assumptions as

(A4) ; (A5) ; (A6) the pair is observable; (A7) ; (A8) .

Remark 1. For the design procedure of ROUIO, a particular form for the matrix C can be chosen as , where is a full rank (nonsingular) matrix. This is not restrictive as long as the matrix C is full row rank (conditions A2A5), and there will always be an orthogonal transformation which leads to the equation .
Either for the FOUIO (2) or the ROUIO (3), system (1) is controlled using the state-feedback controller aswhere 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., as with an arbitrary convergence rate for any x(0), d(t), and u(t), if there exists a positive-definite matrix verifying the matrix inequalityand if the matrix fulfills the following conditions:

Proof. Define the observer reconstruction error vector as the difference between the estimated state described by (2) and the state vector related to system (1) asUsing (1) and (2), the dynamics of the estimation error becomesorUsing , the error dynamics (9) is written in the following form:If conditions (6) are satisfied, the error dynamics of the observer is expressed as . 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 notationthe first term of (6) becomes

Now, using (11) and (12), the dynamics (2) can be written as

If the pair is not observable, the calculation of matrix is made such that the observer is asymptotically stable if and only if 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   ;   Ł,

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 is observable. In [2], it is proved that the above condition is equivalent with the pair observable or at least detectable (first existence condition of the full-order observer). The second constraint is the condition 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

To calculate N, L, J, P, G0, and G1, we partition matrices A0, A1, C, N, P, L, and J as follows:where

The observer’s design is concentrated into the following theorem.

Theorem 4. If the rank conditions in Theorem 3 are satisfied, the UIO (2) can be designed for system (1) as follows:wherewhile and are arbitrarily chosen; should be a nonsingular matrix.

Proof. Since , i.e., the matrix (CD) has full column rank, we can use the pseudoinverse of , given by [17]: By means of this pseudoinverse matrix, we can solve the fourth equation of (6); one successively getsHence, one gets Replacing P into the third equation of (6), we obtain . Now, substituting (14) into the first equation of (6), we obtainBy using (17), system (19) is divided into the following systems:Working in the hypothesis: the matrix is full column rank. By solving systems (20) and (21), we obtainNow, solving the second equation of (6), we getThe unknown matrices of the FOUIO (2), i.e., N, L, J, P, and , 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: there exists a symmetric and positive-definite matrix which verifies the LMI (5), i.e., the matrix N is Hurwitz; the pair 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 as 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 is required in our design approach. Therefore, the unique existence condition of the observer is The theorem’s proof is now complete.

The algorithm associated to the delay-dependent observer is summarized below.

Algorithm 5. Design of FOUIO

Step 1. Given system (1), verify the related assumptions A1-A3.

Step 2. If the existence conditions given in Theorem 3 hold, one calculates the generalized inverse matrix 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 the matrix N is Hurwitz, i.e., there exists a symmetric and positive-definite matrix verifying the LMI (5); the pair 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 verifying the LMI:with , and if there exists a vector such that the next conditions are satisfied:

Proof. Define the observer error vector defined as the difference between the estimated vector described by (3) and the state vector related to system (1) asUsing (1), (3), and (4), the dynamics of the estimation error becomesIf the error dynamics (27) can be written aswhere
Let the state estimation error be defined as . From (1) and (3), one obtains . If , 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 matrices , and , let us partition the matrices A0, A1, C, , and D as follows:where

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: there exists the matrix verifying the LMI (24); the matrices , and are given bywhere the matrices and are arbitrarily chosen, and are submatrices of the matrix , and .

Proof. Using (29) and (30), the fourth equation of (25) becomesTo compute the matrix from this equation, one can choosewith an unknown matrix to be computed. With the notation , one gets , where is the generalized pseudoinverse of , given by According to the assumption A5 and to the fact that the matrix is full rank, we can conclude that , and the generalized pseudoinverse of can be always obtained.
Equation (33) can be then written asand the matrix is chosen arbitrarily.
Now, using the expression of the matrix Z, we obtain the matrices in the third equation of (25). We arbitrarily choose the matrices and ; as a consequence, the matrix is arbitrarily chosen. Using this as well as the form of matrix C given in (29) and (30), i.e., , the first equation (25) can be written in a form of matrix system as follows:by solving system (35), we obtainbearing in mind the fact that , it results that is a submatrix of the nonsingular matrix ; from this, one deduces that , i.e., is a full column rank matrix.
Now, solving the second equation of (25), we getUsing (29) and (30), the sixth equation (25) can be transformed into a matrix system asfrom the first and second equations of (38), we obtain and . On the other hand, the matrix is calculated from the third equation of (35) as follows:where is an unknown matrix to be computed. Because is a nonsingular matrix and , one gets , i.e., is a full column rank matrix. Taking this into account, one successively gets and . From the fourth equation of system (38), we obtain ; by concatenation of the matrices and , 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 matrix of from (29) is chosen; the matrices , and are partitioned with (29) and (30). The matrices and and thus the matrix are arbitrarily chosen; the matrices Z1, , are obtained via (31).

Step 3. The matrices and thus the matrix are arbitrarily chosen; by solving system (35) with respect to the matrices , or, directly, by means of (36), these four matrices are obtained; the matrices and are then calculated with (31). Matrices and are computed by solving system (38).

Step 4. Using the above calculated matrix , one checks if there exists a symmetric and positive-definite matrix verifying the LMI (24); if so, the matrix 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 Figure 1. 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 (I =1…4) is measured by means of a pressure sensor. The process inputs are the input voltages of the pumps ( and ) and the output variables are the tank levels and

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]with the state vector, the input vector, and the output vector, respectively, and with

Using the following numerical valuesthe matrices from (1) are given by

5.3. Software Implementation of the Full-Order Observer

For simulation results, the following numerical values are chosen: , , , , , , and , , The control gain matrix is computed using the approach proposed in [1]

For the FOUIO software implementation, one obtains

5.4. Software Implementation of the Reduced-Order Observer

During the first step of the Algorithm 8, 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 , and We obtain

Figure 2 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.

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.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.