Next Article in Journal
Prescribed Performance Control with Sliding-Mode Dynamic Surface for a Glue Pump Motor Based on Extended State Observers
Previous Article in Journal
Hysteresis Modeling of a PAM System Using ANFIS
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Model-Based Observer Design Considering Unequal Measurement Delays

Department of Mechanical Engineering, University of Alberta, Edmonton, AB T6G 2V4, Canada
*
Author to whom correspondence should be addressed.
Actuators 2021, 10(11), 281; https://doi.org/10.3390/act10110281
Submission received: 8 September 2021 / Revised: 19 October 2021 / Accepted: 20 October 2021 / Published: 22 October 2021
(This article belongs to the Section Control Systems)

Abstract

:
State observers are essential components of a modern control system. It is often designed based on a mathematical model of the process, thus requiring detailed process knowledge. However, in the existing state estimation methods, equal delays are commonly assumed for all communication lines, which is unrealistic and poses problems such as instability and a degraded performance of observers when unequal time delays exist. In this paper, a design of observers considering the measurement delays is presented. To deal with this problem, a chain-based observer has been proposed in which each chain deals with one output delay, performs prediction for the unavailable output value, and passes it to the next chain. Convergence of each chain observer as well as overall state estimation were proven. To illustrate the performance of the proposed scheme, simulation studies were performed on a benchmark continuous stirred tank heater (CSTH) process.

1. Introduction

State observers play a key role in process monitoring, control, and automation by providing estimates of internal state variables using hardware output measurements. Given a process model, observers can be designed by a large number of available methods such as high gain [1], Kalman filter [2], sliding mode [3], and so on. Various types of observers have also been used in the feedback control design in the presence of delays [4,5]. One of the profound challenges in state estimation is that sensor measurements are often time delayed from the actual process outputs. The control problem of dynamic systems with delays has received growing attention in recent years due to novel application areas such as control over networks [6]. This might also occur due to the slow dynamics of physical sensors as well as communication delays from measurements transmitted between different processing units. If these delays are not compensated, they can lead to the performance degradation of the observers and can have a strong impact on their stability and robustness properties [7,8]. When designing an observer or controller for a system affected by such delays, it is often expedient to lump the measurement delay together with any system and actuation delay, as they are indistinguishable in terms of the ideal closed-loop dynamics of the system [9]. However, in some cases, it can be advantageous to model the measurement delays separately from the other delays.
For linear systems, when the system model is available, there exists literature to deal with equal output delays in state estimation [10,11,12]. Assuming equal delays for all output channels is unrealistic due to different communication line lengths or different sensors and technology types. Hence, the involvement of the delay could be in the communication line, be multiple, and different for each channel, which, however, has rarely been studied in the literature. Supposing equal delays will simplify the observer designing step, but can cause similar problems to ignoring them [13]. Others have assumed the delays on outputs [9,14] or states [15], which again, is simpler than considering them in the communication line (out of the process model). Furthermore, an extension of the available model-based methods to data-driven methods is not straightforward (if even it is possible). This paper intends to extend the Luenberger observer design proposed in [16] to account for equal and unequal measurement delays. Then, the subspace-based approach [17] can similarly be extended to design a data-driven observer in the presence of equal and unequal communication delays.
In this paper, we assumed that the output variables were measured but with delays. Furthermore, the delays might be different among the different output variables. Thus, the receiver node needs to be reconstructed to be compatible with the available data. To deal with this problem, a chain observer has been proposed in which each chain collects data, performs prediction, and passes it to the next chain. The chain observer is a common method to deal with delays in nonlinear systems. The idea of achieving the convergence of the state estimate by using a cascade of two observers (elementary chain observer) for a nonlinear system was first proposed in [18], while the idea of using more than two observers in the chain to deal with large measurement delays was proposed in [19]. In [18], it was shown that a chain of two observers is sufficient for asymptotic state reconstruction as long as the measurement delay is below a given threshold, which depends on the Lipschitz constants of the system. When the time delay exceeds such a threshold, more links must be added to the chain [19]. Thus, each observer in the chain is in charge of predicting the system state for a suitable fraction of the total delay. The chain observer method has also been used for dealing with large measurement delays [7].
In the proposed method, the chain observer was used to deal with the multiple un-equal measurement delays. It was proven that the proposed chain observer had the desired stability properties with delay-free measurements, and the design followed a two-step strategy. First, it obtained the delayed state trajectories using delayed measurements with the delay-free observer. These delayed state estimates were then used as predictors that can compensate for the presence of delay. This strategy was shown to converge to true states. In this paper, the convergence of each chain as well as the overall observer was proven. This paper shows the development of the chain observer, which holds stability as well as convergence, along with a new prediction method that is used in each chain.
The main contributions of this paper are as follows:
  • Designing an observer by considering equal measurement delays.
  • Designing a chain observer to deal with unequal measurement delays.
  • Proving the convergence of each chain as well as the overall observer.
The rest of this paper is organized as follows. In Section 2, the preliminaries are presented. The problem formulation is discussed in Section 3. Section 4 presents the theory of the proposed Luenberger-based observer for the case of equal and unequal measurement delays. In Section 5, a simulation through the continuous stirred tank heater (CSTH) process is used to illustrate the performance of the proposed scheme. The paper ends with concluding remarks in the final section.

2. Preliminaries of Model-Based (Luenberger) Observer Design

In this section, we present the preliminaries of designing a model-based observer approach following the work of [17]. Consider a discrete-time LTI (linear time invariant) system, which is described by
x ( k + 1 ) = A x ( k ) + B u ( k ) y ( k ) = C x ( k ) + D u ( k )
where u(k) ∈ Rl, y(k) ∈ Rm, and x(k) ∈ Rn represent process input, output, and state variable vectors, respectively. Suppose that the pair (C, A) with A ∈ Rn × n, C ∈ Rm × n is observable. The s ( = n ) -order observer with the following state and output equations
z(k + 1) = Gz(k) + Hu(k) + Ly(k) ∈ Rs
yˆ(k) = Wz(k) + V y(k) + Qu(k) ∈ Rm
is asymptotically stable, provided that G is stable (all of its eigenvalues lie within the circle one) and the observer parameters H , L , Q , W , V satisfy the following Luenberger equations:
T A − GT = LC, H = T B − LD
C = WT + VC, Q = D − VD
G ∈ Rs × s, T ∈ Rs × n, W ∈ Rm × s
where z ( k ) is the observer states related to the process states through the transformation matrix T ∈ Rs × n so that z ( k ) = T x ^ ( k ) , where x ^ ( k ) is an estimation of true state x ( k ) .

3. Problem Formulation

Sensors and communication networks are an integral part of a data transmission system. One of the main issues regarding measurements and network transmission is the presence of delays in each output measurement (say delay τ i for output y i ). Figure 1 illustrates the communication network and corresponding delays. Due to these delays, output signals y ( k ) required for state estimation (2) are not available instantaneously. Therefore, the objective of this work was to design an observer that utilizes delayed measurements (i.e., y(kτ) = [y1(kτ1) … ym(kτm)]) to estimate z ( k ) . As a special case, these delays may be equal in all outputs (i.e., τ i = d , for i = 1 , , m .) First, we developed an efficient observer scheme in the presence of constant equal measurement delays. The proposed observer, which deals with equal measurement delays, will have the form of:
z ( k ) = G z ( k d ) + i = 0 d 1 H i u ( k 1 i ) + L y ( k d ) T x ^ ( k ) = z ( k )
where it utilizes available data until the ( k d ) t h instant to estimate x ( k ) . The next section will also propose a novel method to deal with unequal measurement delays using a chain–observer structure.

4. The Proposed Observer Design

4.1. Observer Design in the Case of Equal Measurement Delays

In this section, it was assumed that the delays were equal in all outputs (i.e., τ i = d , for i = 1 , , m ). The following theorem designs an asymptotically stable observer (5) with the consideration of equal measurement delays when the system model is available. It is an extension of the Luenberger condition (4) for the case of equal measurement delays.
Theorem 1.
Suppose system model (1) is known with equal measurement delay d in all measured outputs. Furthermore, it is assumed that the pair (C,A) is observable. The observer with the following state and output equations
z ( k ) = G z ( k d ) + i = 0 d 1 H i u ( k 1 i ) + L y ( k d ) R s T x ^ ( k ) = z ( k )
y ^ ( k ) = W z ( k d ) + V y ( k d ) + i = 0 d Q i u ( k i )  
is asymptotically stable, namely l i m k ( T x ( k ) z ( k ) ) 0 , and has unbiased estimation for y ( k ) , that is l i m k ( y ^ ( k ) y ( k ) ) = 0 , provided the observer parameters defined in Equations (6) and (7) satisfy the following delayed-Luenberger equations:
State convergence condition:
T A d G T = L C ,   H d 1 = T A d 1 B L D , H i = T A i B ,   i = 0 , , d 2
Output convergence condition:
CAd = WT + VC, Q0 = D − VD,
Qj = CAj−1B, j = 1,..., d
G ∈ Rs × s, T ∈ Rs × n
where G should be chosen so that it is stable (i.e., every eigenvalue of G lies within the circle one). In addition, the closed loop dynamics of estimation error for state and output is reduced to
e ( k ) = T x ( k ) z ( k ) e ( k ) = G e ( k d ) y ( k ) y ^ ( k ) = W e ( k d ) .
 Proof. 
The state equation of Equation (1) can be rewritten as:
x ( k ) = A x ( k 1 ) + B u ( k 1 )
Substituting x ( k ) from (11) in (1) yields
( k + 1 ) = A 2 x ( k 1 ) + B u ( k ) + A B u ( k 1 )
Similarly, x ( k 1 ) can be replaced by A x ( k 2 ) + B u ( k 2 ) , which yields
x ( k + 1 ) = A 3 x ( k 2 ) + B u ( k ) + A B u ( k 1 ) + A 2 B u ( k 2 )
Continuing the above procedure, after d step reaches
( k + 1 ) = A d + 1 x ( k d ) + i = 0 d A i B u ( k 1 i )
or
x ( k ) = A d x ( k d ) + i = 0 d 1 A i B u ( k 1 i )
The state estimation error is defined as
( k ) = T x ( k ) z ( k )
Substituting Equations (15) and (6) in Equation (16) gives
e ( k ) = T A d x ( k d ) + i = 0 d 1 T A i B u ( k 1 i ) G z ( k d ) i = 0 d 1 H i u ( k 1 i ) L y ( k d )
Now by substitution of T A d G T = L C , H d 1 = T A d 1 B L D , H i = T A i B , and y ( k d ) = C x ( k d ) + D u ( k d ) , the above equation is reduced to
( k ) = G e ( k d )
The estimation error (18) will converge to zero, provided G is a stable matrix. Performing similar steps for y ^ ( k ) y ( k ) gives
y ( k ) y ^ ( k ) = W e ( k d )
Converging e ( k ) to zero proves relation (10). Thus, the proof of Theorem 1 is completed. □
Remark 1.
The matrix G can be defined as follows:
G = [ G 0 g ] ,           g = [ g 1 g s 1 g s ] ,     G 0 = [ 0 0 0 1 0 0 0 1 0 0 0 1 ]  
then, the designed observer poles are roots of polynomial function xs−gsxs−1−···−g1 = 0. Therefore, to have the desired performance, one can select proper poles and determine the corresponding polynomial function coefficients (gi for i = 1,..., s).

4.2. Observer Design in the Case of Unequal Delays

This section considers the extension of the proposed observers (6) and (7) to deal with unequal measurement delay. The proposed observer consists of m chain–observers, each one being in the form of (6) and (7). Without loss of generality, let τ 1 τ 2 τ m and define d 1 = τ 1 ,   d i = τ i τ i 1 ,   i = 2 , , m . Before explaining the proposed approach of designing the observer considering unequal delay, a theorem that provides stability conditions for a system with multiple delays is presented. This theorem is an extension to the single delay case, which was presented in [20]. The result of this theorem enables us to design an observer with the convergence property.
Theorem 2.
Suppose the system model with the assumption A 0 2 0 , is given by
x ( k + 1 ) = A 0 x ( k ) + j = 1 N A j x ( k d j )  
where N is the value of time delays. Then, it is asymptotically stables if there exists a real symmetric matrix P   >   0 so that:
( 1 + ε m ) A 0 T P A 0 + N ( 1 + ε m ε m ) j = 1 N A j T P A j P   <   0 ε m = ( N j = 1 N A j 2 2 ) 1 2 A 0 2 1
where Ψ 2 = λ max ( Ψ T Ψ )
 Proof. 
See Appendix A. □
Remark 2.
When A 0 2 = 0 , system model (21) can be rewritten as
x ( k + 1 ) = j = 1 N A j x ( k d j )       = A 1 x ( k d 1 ) + j = 2 N A j x ( k d j )
Now, define a new variable as x ( k ) = x ( k d 1 ) . Then, the augmented model becomes
[ x ( k + 1 ) x ´ ( k + 1 ) ] = [ 0 0 0 A 1 ] [ x ( k ) x ´ ( k ) ] + j = 2 N [ A j 0 0 0 ] [ x ( k d j ) x ´ ( k d j ) ]
or
x ¯ ( k + 1 ) = A ¯ 0 x ¯ ( k ) + j = 2 N A ¯ j x ¯ ( k d j )
where x ¯ ( k ) = [ x ( k ) x ´ ( k ) ] ,     A ¯ 0 = [ 0 0 0 A 1 ] , and  A ¯ j = [ A j 0 0 0 ] . Then, Theorem2 can be applied for stability analysis of (25).
Suppose that due to different sensors or measurement techniques or communication topologies, each output measurement is received with a delay that may be different from delays in other output channels. Therefore, in order to estimate the state x ( k ) at sample time k , the output y i is only available up to k τ i for i = 1 , , m . The observer gains can be parameterized by the set of the available measurements. In order to handle the unequal measurement delays, we propose constructing a new observer scheme. The proposed observer scheme is shown in Figure 2, in which y i ( 1 : s ) is the output i from sample 1 to s . The proposed observer gains can be calculated in a composition of m chain–observers, each of which predicts the unavailable output due to measurement delay as well as updates the values of the estimated states and then passes them to the next chain.
The following theorem provides the design of the asymptotically stable chain–observer with the consideration of unequal measurement delays provided that the system model is available.
Theorem 3.
Suppose that A, B, C, D, in Equation (1) are known, andτ = [ τ1 τm] represents unequal measurement delays, and define d 1 = τ 1 ,   d i = τ i τ i 1 ,   i = 2 , , m . The design of m -chain–observer is achieved by solving the chain-Luenberger equations for i = 1 , , m chain–observers. The observer i is designed using the estimated output y ^ j + 1 ( k τ i ) for j = i , , m 1 provided by the previous chain–observers from 1 to i 1 as
z ( k τ i 1 ) = G i z ( k τ i ) + L i y ¯ ( k τ i )                                 + i = 0 d i 1 H i u ( k 1 i τ i 1 )  
where
y ^ ( k τ i 1 ) = [ y 1 ( k τ i )     y i ( k τ i )   y ^ i + 1 ( k τ i )   y ^ m ( k τ i )   y ^ m ( k τ i ) ] T i x ^ ( k τ i 1 ) = z ( k τ i 1 )  
and
y ^ ( k τ i 1 ) = W i z ( k τ i ) + V i y ¯ ( k τ i )                                 + i = 0 d i Q i u ( k i τ i 1 )  
Each chain–observer satisfies the following convergence conditions:
T i A d i G j T i = L i C ,     H d i 1 = T i A d i 1 B L i D  
where G i for i = 1 , , m should be chosen so that it is stable and should satisfy the stability condition Theorem 2 as
( 1 + ε i ) G i T P G i + ( m i ) ( 1 + ε i ε i ) j = i + 1 m Φ j T P Φ j P < 0 ε i = ( ( m i ) j = i + 1 m Φ j 2 2 ) 1 2 G i 2 1
where P > 0 , Φ j = L i p = i + 1 m V p W j I j . It follows then, that the constructed observer is asymptotically convergent, and the observer output equation delivers unbiased estimation for y ( k ) , that is,
l i m k ( y ^ ( k ) y ( k ) ) 0
 Proof. 
See Appendix B. □
To design an observer with unequal output delay, Theorem 3 can be utilized to construct each chain observer provided that g i (i.e., G i ) in each chain must be selected in such a way that the stability condition (29) is satisfied. Therefore, the chain observer can be designed using the following algorithm (Algorithm 1):
Algorithm 1. Chain observer design algorithm
Step 1: set i = m
Step 2: Select a proper G i that satisfies stability conditions (29).
Step 3: Design the chain-observer i using Theorem 1 with d = d i .
Step 4: Produce the complementary output y ¯ ( k τ i 1 ) using available output and input, then estimate output data by Equation (27).
Step 5: Go to step 2 and repeat till i = 1 .

5. Simulation

In this section, the effectiveness of the proposed approach was demonstrated on a simulation setup of a continuous stirred tank heater (CSTH). The system is written as a linear state-space model to generate input–output data with the following parameters [16]:
x = [ V T H T T h j ] ,         y = [ h T T T T h j ]
A = [ 0 0 0 626.4371 5.9406 × 10 3 36.55 0 0 1.2019 × 10 3 ]  
B = [ 1 0 0 0 0 3041.5241 ]
C = [ 31.831 0 0 0 3.8578 × 10 5 0 0 0 1 ]
where V ˙ i n V ˙ o u t u 1 ,   P h u 2 are two input variables, and the water level h T , the temperature of the water in the tank T T as well as in the heating jacket T h j are output variables. The physical meanings of the process variables and parameters used above are listed in Table 1.
The output variables had the delay vector of τ = [2 5 7]. The order of the observer was s = 3. Set g1 = [10−2, −3 × 10−8, 3 × 10−4] so that all three poles were located in 10−4 at each chain. Then, we designed a chain–observer using Algorithm 1 when m = 3.
The following is the design procedure:
Chain 1-
Step 1.1 : Set d = 7 5 = 2 and s = 3 ;
Step 1.2 : Set
g = [ 10 12 3 × 10 8 3 × 10 4 ]
G = [ 0 0 10 12 1 0 3 × 10 8 0 1 3 × 10 4 ]  
which results in three poles at 10 4 .
Step 1.3: Calculate L , T , H i as
L = [ 3 × 10 7 2.3 × 10 6 8.3 × 10 12 1.9 × 10 4 6.7 × 10 2 5.7 × 10 6 0.02 0.0172 1.2 × 10 6 ]
H 1 = [ 1.6 × 10 3 2.1 × 10 4 4.1 × 10 4 4.5 × 10 5 1.3 × 10 4 2.4 × 10 5 ]  
T = [ 6.2 × 10 3 2.6 × 10 6 5.7 × 10 6 0.652 6.6 × 10 7 1.2 × 10 6 0.732 2.1 × 10 7 6.7 × 10 7 ]  
H 0 = [ 6.21.6 × 10 3 0.175 0.652 0.038 0.7323 0.02 ]  
Step 1.4 : Calculate V , W , Q i as
W = [ 0.32 1 0.89 162.21 507.3 452.98 246.22 769.4 687.04 ]
V = [ 2 × 10 12 1.1 × 10 7 2.3 × 10 15 4.5 × 10 6 3.5 × 10 5 0 1.4 × 10 9 8.2 × 10 5 1.4 × 10 6 ]
V = [ 2 × 10 12 1.1 × 10 7 2.3 × 10 15 4.5 × 10 6 3.5 × 10 5 0 1.4 × 10 9 8.2 × 10 5 1.4 × 10 6 ]
Q 0 = [ 0 0 0 0 0 0 ]  
Q 1 = [ 31.83 0 0 0 0 30411.5 ]  
Q 2 = [ 0 0 0.0241 0 0 36.55 ]  
Chain 2-
Step 2.1 : Set d = 5 2 = 3 and s = 3 ;
Step 2.2 : Set
g = [ 10 12 3 × 10 8 3 × 10 4 ]
G = [ 0 0 10 12 1 0 3 × 10 8 0 1 3 × 10 4 ]  
which results in three poles at 10 4 .
Step 2.3 : Calculate L , T , H i
L = [ 2.5 × 10 8 2 × 10 4 1.4 × 10 11 9.2 × 10 4 0.998 8.27 × 10 11 1.6 × 10 4 0.0236 1 × 9 ]
H 2 = [ 1.4 × 10 4 3.1 × 10 7 3 × 10 6 4 × 10 7 6.4 × 10 7 3.3 × 10 7 ]                                                 T = [ 2.9 × 10 2 3.8 × 10 5 8.2 × 10 11 5.3 × 10 3 9.1 × 10 7 1 × 10 9 4.5 × 10 5 8.7 × 10 8 1.2 × 10 9 ] H 1 = [ 0.0241 3.2 × 10 8 5.7 × 10 4 3.8 × 10 7 5.4 × 10 5 4.6 × 10 7 ] H 0 = [ 0.0295 2.4 × 10 6 5.3 × 10 3 3 × 10 5 4.6 × 10 5 3.6 × 10 5 ]
Step 2.4 : Calculate V , W , Q i as
W = [ 0.1 0.53 6.35 0.007 0.036 0.428 1.48 7.51 87.86 ]
V = [ 1 × 10 11 3 × 10 6 2.3 × 10 15 2.6 × 10 8 2.1 × 10 7 5 × 10 16 1.5 × 10 10 4.3 × 10 5 1.7 × 10 9 ] Q 0 = [ 0 0 0 0 0 0 ]                                     Q 1 = [ 31.83 0 0 0 0 30411.5 ]                                           Q 2 = [ 0 0 0.0241 0 0 36.55 ]                                           Q 3 = [ 0 0 1.4 × 10 4 0 0 0.044 ]  
Chain 3-Same as Chain 1.
The designed chain–observer (Method 1) was compared with the classic observer (Method 3), which ignores the measurement delays, and the observer (Method 2) [21], which supposes (one) equal measurement delay (delay = 2 in this test). The method proposed in [21] deals with both state and output delays; in this test, state delay is supposed to be equal to zero. Figure 3, Figure 4 and Figure 5 show the estimated values for the system states x1, x2, and x3 by the above-mentioned observers, respectively. These figures show that the proposed observer estimation error converges to zero in 10 s. Observer #2 (considered equal delays) had a bias in estimation and was slow in convergence, and observer #3 (ignored delays) diverged. The input signal to the system is depicted in Figure 6.
Remark 3.
Based on the results presented above, we can conclude that
  • By changing the observer poles, the convergence speed can be tuned.
  • Each chain is an observer in the form of (6).
  • Increasing the communication delays may increase the convergence time and transient estimation error magnitude, but according to theorem 3, it will not make the observer unstable.

6. Conclusions

Delays are an inherent part of any system. Ignoring it could cause instability or degradation in the performances of the control or estimation blocks. Measurement and communication lines produce some delays, which may be unequal due to different line distances or sensor types and technologies. Supposing equal measurement/communication delays causes a similar problem to ignoring them. In this paper, we proposed a new approach to constructing an observer, consisting of a bank of chain–observers to deal with equal and unequal measurement delays. Analytical results for the construction of the observers are provided, and convergence of the chain–observers is shown. Simulations demonstrate the efficacy of the proposed design scheme through a CSTH system in the presence of measurement delays. Having a chain-based observer design facilitates the extension of the proposed approach to a data-driven method. Therefore, the advantage of the proposed approach is that it can be used for designing data-driven (subspace-based) observers [17] in the presence of equal and unequal communication delays. Designing such a data-driven observer is in our future research scope.

Author Contributions

Conceptualization, Y.A. and L.Z.; Methodology, Y.A. and L.Z.; Software, Y.A.; Validation, Y.A. and L.Z.; Formal analysis, Y.A.; Investigation, Y.A.; Resources, Y.A. and L.Z.; Data curation, Y.A.; Writing—original draft preparation, Y.A.; Writing—review and editing, Y.A. and L.Z.; Visualization, Y.A.; Supervision, L.Z.; Project administration, L.Z.; Funding acquisition, L.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research is supported by Alberta Innovates with award number 202102493.

Data Availability Statement

Data sharing is not applicable (no new data were created or analyzed in this study, therefore, data sharing was not applicable to this article).

Acknowledgments

The authors gratefully acknowledge the financial support from Alberta Innovates.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Proof of Theorem 2

Before proving Theorem 2, the following lemma is presented.
Lemma A1.
For any two matrices F and G of dimension n × m and for any square matrix P = P T > 0 of dimension n , the following matrix inequality holds,
( F + G ) T P ( F + G ) ( 1 + ε ) F T P F + ( 1 + ε 1 ) G T P G  
where ε is a positive constant.
Now, to prove Theorem 2, consider the following Lyapunov function:
V ( x ( k ) ) = x ( k ) T P x ( k ) + j = 1 N l = 1 d j x T ( k l ) S j x ( k l ) P = P T > 0 S j = S j T 0
Taking forward the difference from the above equation gives:
V ( x ( k ) )   V ( x ( k + 1 ) ) V ( x ( k ) )                                                         = [ A 0 x ( k ) + j = 1 N A j x ( k d j ) ] T P [ A 0 x ( k ) + j = 1 N A j x ( k d j ) ]                                                         x ( k ) T P x ( k ) + j = 1 N l = 1 d j x T ( k + 1 l ) S j x ( k + 1 l )                                                         j = 1 N l = 1 d j x T ( k l ) S j x ( k l ) P = P T > 0
S j = S j T 0
Applying Lemma A1 on inequality (A3) yields
Δ V ( x k ) ( 1 + ε ) x T ( k ) A 0 T P A 0 x ( k ) +   + ( 1 + ε 1 ) j = 1 N x T ( k d j ) A j T P j = 1 N A j x ( k d j )   x T ( k ) P x ( k ) + x T ( k ) j = 1 N S j x ( k )   j = 1 N x T ( k d j ) S j x ( k d j )
Now, based on Cauchy-Schwarz inequality [22], namely for any real vector v i , the inequality ( i = 1 m v i ) T ( i = 1 m v i ) m i = 1 m v i T v i holds, and we have
Δ V ( x k ) x T ( k ) [ ( 1 + ε ) A 0 T P A 0 + j = 1 N S j P ] x ( k )   + ( 1 + ε 1 ) N j = 1 N x T ( k d j ) A j T P A j x ( k d j )   j = 1 N x T ( k d j ) S j x ( k d j )
Taking two factors x T ( k d j ) and x ( k d j ) out of the left- and right-hand side of two last terms in Equation (A5) yields
Δ V ( x k ) x T ( k ) [ ( 1 + ε ) A 0 T P A 0 + j = 1 N S j P ] x ( k )   + j = 1 N x T ( k d j ) [ N ( 1 + ε 1 ) A j T P A j S j ] x ( k d j )
Now, select S j = N ( 1 + ε 1 ) A j T P A j . We have
Δ V ( x k ) x T ( k ) [ ( 1 + ε ) A 0 T P A 0   + N ( 1 + ε 1 ) j = 1 N A j T P A j P ] x ( k )       = φ ( x , ε )
Since matrices A 0 T P A 0 and A j T P A j are symmetric and positive semi-definite, based on [23], we have
φ ( x , ε ) x T ( k ) [ ( 1 + ε ) λ max ( A 0 T P A 0 )   + N ( 1 + ε 1 ) j = 1 N λ max ( A j T P A j ) P ] x ( k )           = g ( ε ) λ max ( P ) x ( k ) 2 2
where
g ( ε ) = ( 1 + ε ) σ max 2 ( A 0 ) + N ( 1 + ε 1 ) j = 1 N σ max 2 ( A j )
The minimum of the function g can be found by taking its derivative with respect to the parameter ε , which yields
d g ( ε ) d ε = 0 σ max 2 ( A 0 ) N ε 2 j = 1 N σ max 2 ( A j ) = 0
Therefore, the optimum value for ε is
ε m = ( N j = 1 N σ max 2 ( A j ) ) 1 2 σ max 1 ( A 0 )
We can conclude
Δ V ( x k ) φ ( ε m , x k )   = x T ( k ) [ ( 1 + ε m ) λ max ( A 0 T P A 0 )   + N ( 1 + ε m 1 ) j = 1 N λ max ( A j T P A j ) P ] x ( k )
Now, if the condition ( 1 + ε m ) A 0 T P A 0 + N ( 1 + ε m ε m ) j = 1 N A j T P A j P < 0 is satisfied, then considering Equation (A8), we have
Δ V ( x k ) β x ( k ) 2 2
The proof is completed.

Appendix B. Proof of Theorem 3

From Equations (1) and (27), we have
y ( k ) = C x ( k ) + D u ( k )
and
y ^ ( k ) = W i z ( k d ) + V i y ¯ ( k d )                                 + i = 0 d i Q i u ( k i )
Substituting x ( k ) from Equation (A14) in (A15), gives
y ( k ) = C A d x ( k d ) + i = 0 d 1 C A i B u ( k 1 i ) + D u ( k )
Now, subtracting (A16) from (A15), gives
y ( k ) y ^ ( k ) = C A d x ( k d ) + i = 0 d 1 C A i B u ( k 1 i ) + D u ( k ) W i z ( k d ) + V i y ¯ ( k d ) i = 0 d i Q i u ( k i )
or
y ( k τ i 1 ) y ^ ( k τ i 1 )                                                         = C A d x ( k d τ i 1 ) + i = 0 d 1 C A i B u ( k 1 i τ i 1 )                                                         + D u ( k τ i 1 ) W i z ( k d τ i 1 ) V i y ¯ ( k d τ i 1 )                                                         i = 0 d i Q i u ( k i τ i 1 )
Now, recalling that τ i = d + τ i 1 , the above equation can be rewritten as
y ( k τ i 1 ) y ^ ( k τ i 1 )                                                         = C A d x ( k τ i ) + i = 0 d 1 C A i B u ( k 1 i τ i 1 )                                                         + D u ( k τ i 1 ) W i z ( k τ i ) V i y ¯ ( k τ i )                                                         i = 0 d i Q i u ( k i τ i 1 )
Now, according to Equation (10), for chain i , we can write
y ^ ( k τ i 1 ) = y ( k τ i 1 ) W i e ( k d i )   f o r   i = 1 , , m
Then, we have
y ¯ ( k τ i 1 ) = y ( k τ i 1 ) W i I i e ( k d i ) .  
From (A19), for each chain i, we get
y ( k τ i 1 ) y ^ ( k τ i 1 ) = C A d i x ( k τ i )   + j = 0 d i 1 C A j B u ( k j τ i 1 ) + D u ( k τ i 1 )   W i 1 z ( k τ i ) V i 1 y ¯ ( k τ i 1 )   j = 0 d i 1 Q j u ( k j τ i 1 ) .  
Substituting Equation (A21) in the above equation yields
y ( k τ i 2 ) y ^ ( k τ i 2 ) = C A d i 1 x ( k τ i 1 )   + j = d i 1 d i 2 1 C A j B u ( k j τ i 1 ) + D u ( k τ i 2 )   W i 1 z ( k τ i 1 ) V i 1 y ( k τ i 1 )   V i 1 W i I i e ( k τ i ) j = d i 1 d i 2 1 Q i u ( k j τ i 2 )
Now, selecting C A d i 1 = W i 1 T i 1 + V C , Q d i 1 = D V i 1 D , Q j = C A j B , j = d i + 1 , , , d i 1 gives
y ( k τ i 2 ) y ^ ( k τ i 2 ) = V i 1 W i e ( k τ i )                   + W i 1 e ( k τ i 1 )
Then,
y ¯ ( k τ i 1 ) = y ( k τ i 1 ) V i 1 W i I i e ( k τ i )               W i 1 I i 1 e ( k τ i 1 )
Comparing Equations (A20) and (A25), and following similar steps from (A20) to (A23) on (A25), we will reach:
y ( k τ i 3 ) y ^ ( k τ i 3 ) = V i 2 V i 1 W i e ( k τ i ) + V i 1 W i 1 e ( k τ i 1 ) + + W i 2 e ( k τ i 2 )                  
Repeating this sequence until chain m , gives
y ( k τ i ) y ^ ( k τ i ) = j = i + 1 m p = m j V p W j e ( k τ j )
and
y ¯ ( k τ i 1 ) = y ( k τ i 1 ) j = i + 1 m p = m j V p W j I j e ( k τ j )
According to Equation (A27), the convergence of states brings
l i m k y ( k τ i ) y ^ ( k τ i ) = 0
Now, the estimation error can be written as
e ( k + 1 ) = T i x ( k + 1 ) z ( k + 1 )
Substituting Equations (A22) and (26) in Equation (A30) gives
e ( k + 1 ) = T m A d m + 1 x ( k τ m )   + j = 0 d i T i A j B u ( k j τ m ) G m z ( k τ m )   j = 0 d i H j u ( k j τ m ) L i y ¯ ( k τ m )
Now, substituting T i A d i G i T i = L i C , H d i 1 = T i A d i 1 B L i D , H j = T i A j B , j = 0 , , d i 2 and considering Equation (26) to (29), the above equation is reduced to
e ( k + 1 ) = G m e ( k τ m ) L i j = i + 1 m p = m j V p W j I j e ( k τ j )
From Equation (18) in Theorem 1 and Remark 2, we know that
e ( k ) = G e ( k d )
Considering the stability of G, we have lim k e ( k ) = 0 . Then, considering all matrices in Equation (A32) are constant, we can conclude
lim k y ( k τ i ) y ^ ( k τ i ) = lim k e ( k τ i )   0
Therefore, the estimation error will converge to zero, provided that the error dynamics (A32) is stable (i.e., condition (29) is satisfied).

References

  1. Boizota, N.; Busvelleb, E.; Gauthier, J.-P. An adaptive high-gain observer for nonlinear systems. Automatica 2010, 46, 1483–1488. [Google Scholar] [CrossRef]
  2. Youngjoo, K.; Bang, H. Introduction to Kalman filter and its applications. In Introduction and Implementations of the Kalman Filter; Govaers, F., Ed.; IntechOpen: Rijeka, Croatia, 2019. [Google Scholar]
  3. Spurgeon, S.K. Sliding mode observers: A survey. Int. J. Syst. Sci. 2008, 39, 751–764. [Google Scholar] [CrossRef]
  4. Regaieg, M.A.; Kchaou, M.; Bosche, J.; El-Hajjaji, A.; Chaabane, M. Robust dissipative observer-based control design for discrete-time switched systems with time-varying delay. IET Control Theory Appl. 2019, 13, 3026–3039. [Google Scholar] [CrossRef]
  5. Kchaou, M.; Al Ahmadi, S. Robust control for nonlinear uncertain switched descriptor systems with time delay and nonlinear input: A sliding mode approach. Complexity 2017, 2017, 1027909. [Google Scholar] [CrossRef] [Green Version]
  6. Karafyllis, I.; Malisoff, M.; Mazenc, F.; Pepe, P. Recent Results on Nonlinear Delay Control Systems; Springer International Publishing: Cham, Switzerland, 2016. [Google Scholar]
  7. Nagmani, G.; Karthik, C.; Rajan, G.S. Observer-Based Exponential Stabilization for Time-Delay Systems via Augmented Weighted Integral Inequality. J. Frankl. Inst. 2019, 356, 9023–9042. [Google Scholar] [CrossRef]
  8. Gnaneswaran, N.; Joo, Y.H.; Kim, H.S. A linear matrix inequality-based extended dissipativity criteria for linear systems with additive time-varying delays. IFAC J. Syst. Control 2019, 10, 100070. [Google Scholar] [CrossRef]
  9. Won, J.K. An Observer for State Estimation in the Presence of Multiple Measurement Delays. Master’s Thesis, California State University, Long Beach, CA, USA, 2019. [Google Scholar]
  10. Olbrot, A. Observability and observers for a class of linear systems with delays. IEEE Trans. Autom. Control 1981, 26, 513–517. [Google Scholar] [CrossRef]
  11. Zheng, G.; Bejarano, F.J. Observer design for linear singular time-delay systems. Automatica 2017, 80, 1–9. [Google Scholar] [CrossRef]
  12. Sanz, R.; Garcia, P.; Krstic, M. Observation and stabilization of LTV systems with time-varying measurement delay. Automatica 2019, 103, 573–579. [Google Scholar] [CrossRef]
  13. Cacace, F.; Germani, A.; Manes, C. A chain observer for nonlinear systems with multiple time-varying measurement delays. SIAM J. Control Optim. 2014, 52, 1862–1885. [Google Scholar] [CrossRef]
  14. Subbarao, K.; Muralidhar, P.C. A state observer for lti systems with delayed outputs: Time-varying delay. In Proceedings of the 2008 American Control Conference, Seattle, WA, USA, 11–13 June 2008. [Google Scholar]
  15. Cao, Y.Y.; Sun, Y.X.; Cheng, C. Delay-dependent robust stabilization of uncertain systems with multiple state delays. IEEE Trans. Autom. Control 1998, 43, 1608–1612. [Google Scholar]
  16. Ding, S.X. Data-Driven Design of Fault Diagnosis and Fault-Tolerant Control Systems; Springer: London, UK, 2014. [Google Scholar]
  17. Ding, S.X.; Yin, S.; Wang, Y.; Wang, Y.; Yang, Y.; Ni, B. Data-driven design of observers and its applications. In Proceedings of the 18th IFAC World Congress, Milano, Italy, 28 August–2 September 2011. [Google Scholar]
  18. Germani, A.; Manes, C.; Pepe, P. State observation of nonlinear systems with delayed output measurements. In Proceedings of the 2nd IFAC Workshop on Time Delay Systems, Ancona, Italy, 11–13 September 2000; pp. 58–63. [Google Scholar]
  19. Kahelras, M.; Ahmed-Ali, T.; Giri, F.; Lamnabhi-Lagarrigue, F. Sampled-Data Chain-Observer Design for a Class of Delayed Nonlinear Systems. Int. J. Control 2018, 91, 1076–1090. [Google Scholar] [CrossRef]
  20. Stojanovic, S.B.; Debeljkovic, D.L. On the asymptotic stability of linear discrete time delay systems. Eng. Syst. Des. Anal. 2004, 41731, 805–811. [Google Scholar]
  21. Dinh Cong, H. A fresh approach to the design of observers for time-delay systems. Trans. Inst. Meas. Control 2018, 40, 477–503. [Google Scholar]
  22. Callebaut, D. Generalization of the Cauchy-Schwarz inequality. J. Math. Anal. Appl. 1965, 12, 491–494. [Google Scholar] [CrossRef]
  23. Meyer, C.D. Matrix Analysis and Applied Linear Algebra; Siam: Philadelphia, PA, USA, 2000; Volume 71. [Google Scholar]
Figure 1. Schematic showing the transmission lines with measurement delays.
Figure 1. Schematic showing the transmission lines with measurement delays.
Actuators 10 00281 g001
Figure 2. The proposed m chain–observer (yi (1: s) is the output i from sample 1 to sample s).
Figure 2. The proposed m chain–observer (yi (1: s) is the output i from sample 1 to sample s).
Actuators 10 00281 g002
Figure 3. The state 1 estimation results of the three introduced observers.
Figure 3. The state 1 estimation results of the three introduced observers.
Actuators 10 00281 g003
Figure 4. The state 2 estimation results of the three introduced observers.
Figure 4. The state 2 estimation results of the three introduced observers.
Actuators 10 00281 g004
Figure 5. The state 3 estimation results of the three introduced observers.
Figure 5. The state 3 estimation results of the three introduced observers.
Actuators 10 00281 g005
Figure 6. Shows inputs signals of the CSTH system.
Figure 6. Shows inputs signals of the CSTH system.
Actuators 10 00281 g006
Table 1. The process variables and parameters.
Table 1. The process variables and parameters.
SymbolDescriptionUnit
VTWater volume in the tankL
HTEnthalpy in the tankJ
ThjTemperature in the heating jacket°C
V ˙ i n ,   V ˙ o u t Water flows in and out of the tank1/s
PhElectrical heater powerW
hTWater level in the tankm
TTWater temperature in the tank°C
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Alipouri, Y.; Zhong, L. Model-Based Observer Design Considering Unequal Measurement Delays. Actuators 2021, 10, 281. https://doi.org/10.3390/act10110281

AMA Style

Alipouri Y, Zhong L. Model-Based Observer Design Considering Unequal Measurement Delays. Actuators. 2021; 10(11):281. https://doi.org/10.3390/act10110281

Chicago/Turabian Style

Alipouri, Yousef, and Lexuan Zhong. 2021. "Model-Based Observer Design Considering Unequal Measurement Delays" Actuators 10, no. 11: 281. https://doi.org/10.3390/act10110281

APA Style

Alipouri, Y., & Zhong, L. (2021). Model-Based Observer Design Considering Unequal Measurement Delays. Actuators, 10(11), 281. https://doi.org/10.3390/act10110281

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop