Next Article in Journal
Data-Driven Adaptive Control for Laser-Based Additive Manufacturing with Automatic Controller Tuning
Next Article in Special Issue
Trajectory/Path-Following Controller Based on Nonlinear Jerk-Level Error Dynamics
Previous Article in Journal
Cereal Stem Stress: In Situ Biomechanical Characterization of Stem Elasticity
Previous Article in Special Issue
Particle Swarm Optimization of a Passivity-Based Controller for Dynamic Positioning of Ships
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Observer Design for Nonlinear Invertible System from the View of Both Local and Global Levels

1
Guizhou Provincial Key Laboratory of Internet Plus Collaborative Intelligent Manufacturing, Electrical Engineering School, Guizhou University, Guiyang 550025, China
2
University de Toulouse, UPS, LAAS, F-31400 Toulouse, France
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(22), 7966; https://doi.org/10.3390/app10227966
Submission received: 22 September 2020 / Revised: 28 October 2020 / Accepted: 4 November 2020 / Published: 10 November 2020
(This article belongs to the Special Issue Nonlinear Control: Theory and Applications)

Abstract

:
This paper emphasizes the importance of the influences of local dynamics on the global dynamics of a control system. By considering an actuator as an individual, nonlinear subsystem connected with a nonlinear process subsystem in cascade, a structure of interconnected nonlinear systems is proposed which allows for global and local supervision properties of the interconnected systems. To achieve this purpose, a kind of interconnected observer design method is investigated, and the convergence is studied. One major difficulty is that a state observation can only rely on the global system output at the terminal boundary. This is because the connection point between the two subsystems is considered unable to be measured, due to physical or economic reasons. Therefore, the aim of the interconnected observer is to estimate the state vector of each subsystem and the unmeasurable connection point. Specifically, the output used in the observer of the actuator subsystem is replaced by the estimation of the process subsystem observer, while the estimation of this interconnection is treated like an additional state in the observer design of the process subsystem. Expression for this new state is achieved by calculating the derivatives of the output equation of the actuator subsystem. Numerical simulations confirm the effectiveness and robustness of the proposed observer, which highlight the significance of the work compared with state-of-the-art methods.

1. Introduction

A modern system often consists of a series of interconnected dynamical units (sensors, actuators and system components) and, therefore, very complicated dynamics are exhibited [1]. Technological advances mean that these units themselves are dynamic systems and exhibit complicated dynamics. Therefore, a modern control system can be viewed as composed of dynamic subsystems connected in a series. In all situations, the global plant can be analyzed at different levels—down to the component level—to estimate the reliability of the whole plant.
In recent years, the topic of observer design for separate nonlinear systems has been widely discussed in the literature, like high-gain observers [2,3,4], sliding mode observers [5,6], adaptive observers [7,8,9], and unknown input observers (UIOs) [10,11,12]. These methodologies are typically centralized monitoring systems where intelligence is either at the system level or at the field device level of the processing plant. For the former, it aims at monitoring plant dynamics from the viewpoint of a global system, like in [13,14,15]. In these methods, the dynamics of subcomponents (i.e., actuators) are often neglected. They are, generally speaking, treated as constants in the input or output coefficient matrix (function) of the process system model. For the latter, it focuses on the field device level, aiming at analyzing the internal dynamics of a specific subcomponent while the influences of local internal dynamics on the global dynamics are neglected, like in [3,16].
However, centralized observers may not be suitable for modern control systems. On one hand, a modern control system is in fact an interconnected system, while the centralized observer just enables an individual component to monitor internal dynamics locally. However, the dynamics of the field devices can cause significant disturbances to the global process and influence the quality of the final product [17,18]. On the other hand, due to uneconomical measurement costs or physical environment factors such as high temperature, it is impossible to measure the state or partial state, such as in [19,20]. To overcome these difficulties, an effective way is to decompose the system into several interconnected subsystems so that the observer can be decentralized in each subsystem. In this way, it allows the analysis of less complex subcomponents to study the characteristics of interconnected systems.
The last few decades have also witnessed significant improvements in dynamic networks which consist of a number of interconnected units, like in [1,21,22,23,24]. A typical approach is to design distributed observers for each subsystem using the internal information of each subsystem, and then all the observers are aggregated to form the total estimator [25,26,27]. For instance, in [26], an observer was designed for the whole system from the separate synthesis of observers for each subsystem, assuming that for each of these separate designs, the states from the other subsystem were available. In [27], an observer for each subsystem was proposed, using the state estimation of the previous subsystem. In addition, a quasi-input-to-state stability and input-to-state dynamical stability (ISS/ISDS) reduced-order observer for the whole system was designed, considering the interconnections of quasi-ISS/ISDS reduced-order observers for each subsystem. A major challenge for these methods is the availability of the measurement of the interconnections between subsystems. Therefore, it is interesting to consider the problem of whether we can prove that, under some conditions, the effect in lower subsystems can be distinguished from higher subsystems, thus avoiding full measurements of the local subsystem. For example, in this work, the interconnection is the output of the actuator, and it is not economical or realistic to measure its output. Contributions dealing with the state observation problem for interconnected systems subjected to unknown interconnections have received less extensive treatment in the literature. In [28], a promising method to solve the state observation problem of nonlinear systems modeled by ODE-PDE series was proposed. A similar problem was also studied in [29], where the interconnected system was composed of a nonlinear system and a linear system.
In this paper, the problem of state estimation for interconnected nonlinear dynamic systems is studied. An interconnected system consists of two nonlinear dynamic subsystems, and the interconnection point is unknown. Thus, one major difficulty is that state observation can only rely on the output of the terminal subsystem, making existing observers useless. Therefore, the problem considered here is that the output of the nonlinear system cannot be measured directly, while part of the state measurement of the second nonlinear system can be obtained. Two issues are highlighted here. Firstly, it is assumed that the measurement value used by the observer of the former subsystem is unmeasurable, and the solution is to replace it with the estimated value of the observer of the latter subsystem. Secondly, in the latter subsystem, the estimated interconnection provided for the former subsystem is regarded as an additional state to form a new, extended subsystem. The expression of the new state is obtained by calculating the derivative of the output equation of the former subsystem.
The contribution of this paper mainly lies in its emphasis of the importance of the influences of local internal dynamics (actuator) on the global dynamics of a control system. A method is proposed to distinguish the influence of low subsystems in higher subsystems, even if the full measurement of local subsystems cannot be realized. Thus, the goal of the design methodology is to enable or simplify observer design for systems that are otherwise difficult to handle by allowing the designer to focus on a smaller, nonlinear subsystem. That is to say, we mainly focus on observing, for example, how the change of an internal parameter at the local level affects the global output at the global level. As a result, both local and global dynamics supervision is achieved, as well as analysis of the influences of local internal dynamics on the global dynamics.
The rest of the paper is organized as follows. Problem formulation is introduced in Section 2, where the type of dynamic units of the interconnected system is explained, and the main objective is introduced. Section 3 contains all the results for the observer design, with respect to the interconnected systems. Some numerical simulation examples are given to illustrate the effectiveness of the proposed methods in Section 4. Finally, a conclusion is made in Section 5.

2. Motivations and Problem Formulations

The problem of state observation is investigated for an interconnected nonlinear system, modeled by two cascaded nonlinear dynamical subsystems: the process and the actuator subsystems. As shown in Figure 1, an interconnected nonlinear system structure is proposed by considering both the actuator and the process as individual dynamic subsystems connected in cascade. The aim is to accurately estimate the state vector of both subsystems, as well as the interconnection. As a result, both local and global dynamics supervision are realized, and analysis of the influence of local internal dynamics on global dynamics is achieved.
A dynamic process subsystem is proposed in an input affine form:
Σ p : { x ˙ = f ( x ) + g ( x ) u a , x ( t 0 ) = x 0 y = x  
where x n is the state of the process subsystem, y p is the output of the global system, which is also the output of the process subsystem, u a m is the input of the process subsystem, which is also the output of the actuator subsystem, u a is assumed to be inaccessible, and f and g are smooth vector fields on n .
The actuator subsystem is described as
Σ a : { x ˙ a = f a ( x a , u ) x a ( t 0 ) = x a 0 u a = h a ( x a , u )
where x a n is the state, u   l is the input and constant parameter to be monitored, u a R m is the output of the actuator subsystem, which is also the input of the process subsystem, f a is the smooth vector field on   n , and h a is the smooth vector field on m .
Considering the interconnected system described by Equations (1) and (2), it is required to monitor the performance of the interconnected system from the perspective of a single subsystem and the whole system; that is, it is required to describe the cause and effect relationships between the subsystem variables and global system output y, thus providing advanced predictive maintenance techniques in an operating plant. The left invertibility of the interconnected system is then required for ensuring that the impact of local variables on the global level is distinguishable. The property of distinguishability of the two inputs or parameters refers to their capacity to generate different output signals for a given input signal.
One way to achieve this purpose is to propose observers for each of the subsystems and the whole network. However, the main difficulty is that the connection point between the two subsystems cannot be measured. This is because the interconnection point is the output of the actuator subsystem, and online measurement is difficult to achieve due to physical or uneconomical reasons. In addition, the measured value could be unreliable, due to its rough operation environment.
Therefore, the state observation in this work can only rely on the global system output (i.e., the process state at the terminal boundary). As shown in Figure 1, the particular aim of our design is to accurately estimate the state vectors x and x a of each subsystem online, as well as the unmeasured interconnection vector u a .

3. Observer Design

The structure of the proposed interconnected observer is shown in Figure 2. It is a two-stage interconnected observer system, consisting of an actuator and a process state estimator. The actuator state estimator deals with estimating the states of the actuator subsystem, where the major challenge is that the output is inaccessible. Aiding the actuator, the process state estimator is a coordinator that extends the interconnection as an additional state of the process subsystem. This process state estimator generates an input sequence which is applied to the actuator subsystem. Then, the overall observer estimates the states and interconnections of the interconnected system by using the estimates of the two estimators.
The main idea of the interconnected observer design is as follows. In the first aspect, the unknown interconnection is extended as new states of the process subsystem, where the expression can be achieved via derivatives of the output expression of the actuator subsystem, thus forming a new process subsystem. Then, for the state estimation of each subsystem of the interconnected system, the state estimation of one subsystem is realized by the state estimation of the other subsystem, and the global estimator is formed by the set of both observers. Specifically, it is assumed that an existing observer is already available for the actuator subsystem a , where the measured output is u a , while the observer is implemented using an estimate of u a , denoted by u ˜ a . In order to obtain this estimate, the state space of the process subsystem p is extended to include u a as an additional state. By calculating the derivatives of the value u a of the actuator subsystem, the expression of the time derivatives of u a is obtained; it is a function of u , with derivatives of u and x a . In summation, for the studied interconnected nonlinear systems, an interconnected observer design method is proposed by combining both actuator and process subsystem state estimators.

3.1. Observer Design for the Interconnected System

3.1.1. Interconnected System Extension

For the interconnected system described by Equations (1) and (2), in order to facilitate analysis, the unknown interconnection u a is extended as a new state x u :
{ x u : u a x ˙ u : = u ˙ a
To get a function for x ˙ u , inspired by the work proposed in [29], let us derivate the output u a in Equation (2) to get the following function:
x ˙ u : = u ˙ a = ε ( u , u ˙ , x a ) = h a u ( u , x a ) u ˙ + h a x a ( u , x a ) f a ( u , x a )
where the function ε ( u , u ˙ , x a ) is with respect to the time derivative of the output u a in Equation (2).
We can define Assumption 1 as follows: For any function u U ( t , x ^ a ) A . C ( + , ) , there exists a real constant satisfied by γ 1 :
ε ( u , u ˙ , x ^ a ) ε ( u , u ˙ , x a ) γ 1 x a x ^ a
Assumption 1.
It refers to the global Lipchitz-type condition of function ε , although this condition seems restrictive, becomes much lower since u and u a are bounded, which is usually the case in physical situations.
Moreover, this boundedness can be found by introducing saturation in the argument of h a . According to [29], if u and h a belong to a compact set U , Y , then the global Lipchitz-type condition of function ε can be replaced by local smoothness by using saturations.
Thus, a new interconnected system is constituted of Σ : = Σ p + Σ a + x u :
Σ : { x ˙ = f ( x ) + g ( x ) x u x ˙ u = ε ( u , u ˙ , x a ) x ˙ a = f a ( x a , u ) y = x x ( t 0 ) = x 0 ;   x a ( t 0 ) = x a 0 ;   x u = h a ( x a , u )
where the input of the system is u , the output is y , and x u is an unmeasured state.
Let ξ 1 = [ ξ 11 ξ 12 ] = [ x x u ] , ξ 2 = x a . Thus, the above system becomes
Σ : { ξ ˙ 1 = f ¯ 1 ( ξ 1 , ξ 2 , u ) ξ ˙ 2 = f ¯ 2 ( ξ 2 , u ) ( y 1 y 2 ) = ( h 1 ( ξ 11 ) h 2 ( ξ 2 ) )
where f ¯ 1 ( ξ 1 , ξ 2 , u ) : = ( f ( ξ 11 ) + g ( ξ 11 ) ξ 12 ε ( u , u ˙ , ξ 2 ) ) , f ¯ 2 ( ξ 2 , u ¯ ) : = f a ( ξ 2 , u ) , and ( h 1 ( ξ 11 ) h 2 ( ξ 2 ) ) : = ( ξ 11 h a ( x a , u ) ) .
The above system can be divided into two subsystems:
Σ O i : { ξ ˙ i = f ¯ i ( ξ i , ξ i ¯ , u ) y i = h i ( ξ i )
where i = { 1 , 2 } and i ¯ denotes the complementary index of i (i.e., { i , i ¯ } = { 1 , 2 } ).

3.1.2. State Estimator Design for the New Process Subsystem

The following system can be viewed as a transformed form of the process subsystem described in Equation (1):
Σ P : { ξ ˙ 1 = f ¯ 1 ( ξ 1 , ξ 2 , u ) y 1 = ξ 11
where f ¯ 1 ( ξ 1 , ξ 2 , u ) : = ( f ( ξ 11 ) + g ( ξ 11 ) ξ 12 ε ( u , u ˙ , ξ 2 ) ) .
The new process subsystem in Equation (9) can be expressed as follows:
Σ P : { ξ ˙ 1 = G ( ξ 11 ) ξ 1 + F ( ξ 11 ) + ε ¯ ( u , u ˙ , ξ 2 ) y 1 = C ξ 1
where G ( ξ 11 ) = ( 0 g ( ξ 11 ) 0 0 ) ,   F ( ξ 11 ) = ( f ( ξ 11 ) 0 ) , C ( I n 0 ) , ε ¯ ( u , u ˙ , ξ 2 ) = [ 0 ε ( u , u ˙ , ξ 2 ) ] T , and I n is the n × n identity matrix.
As demonstrated in [30], supposing that the following assumptions related to the boundedness of the states, signals, and functions are satisfied, an extended, high-gain observer for the system in Equation (10) can be formed.
Assumption 2.
It states that there exist finite real numbers ϱ ,   τ with 0 < ϱ τ , and that ϱ 2 I n F T ( ξ 11 ) F ( ξ 11 ) τ 2 I n .
Assumption 3.
Is that F ( ξ 11 ) is a global Lipchitz, with respect to ξ 11 , locally and uniformly with respect to u.
Assumption 4.
It states that g ( ξ 11 ) is a global Lipchitz with respect to ξ 11 .
Then, an extended, high-gain observer for the system in Equation (10) can be given as
Σ PO : { ξ ^ ˙ 1 = G ( ξ ^ 11 ) ξ ^ 1 + F ( ξ ^ 11 ) + ε ¯ ( u , u ˙ , ξ ^ 2 ) Λ 1 ( ξ ^ 1 ) S θ 1 C T ( C ξ ^ 1 y 1 ) y 1 = C ξ ^ 1
where H ( ξ ^ 1 ) : = Λ 1 ( ξ ^ 1 ) S θ 1 C T is the gain function and Λ ( ξ ^ 1 ) : = [ I 0 0 g ( ξ ^ 11 ) ] , S θ is the unique symmetric positive definite matrix, satisfying the following algebraic Lyapunov equation:
θ S θ + A T S θ + S θ A C T C = 0
where A = [ 0 I 0 0 ] , θ > 0 is a parameter defined by Equation (12), and the solution is
S θ = [ 1 θ I 1 θ 2 I 1 θ 2 I 2 θ 3 I ]  
Then, the gain of the estimator can be given by
H ( ξ ^ 1 ) = Λ 1 ( ξ ^ 1 ) S θ 1 C T = [ 2 θ I θ 2 g 1 ( ξ ^ 11 ) ]
The state estimation error is expressed as
e p ( t ) = ξ ^ 1 ( t ) ξ 1 ( t )
Then, by subtracting the corresponding Equations (10) and (11), the following error dynamics can be obtained:
e ˙ p = G ( ξ ^ 11 ) ξ ^ 1 + F ( ξ ^ 11 ) + ε ¯ ( u , u ˙ , ξ ^ 2 ) + H ( ξ ^ 1 ) ( C ξ ^ 1 y ) G ( ξ 11 ) ξ 1 F ( ξ 11 ) ε ¯ ( u , u ˙ , ξ 2 )
Assumption 5.
Is that for any u U ( t , e p ) A . C ( + , ) , there exists a continuously differentiable function V p , and the positive constants α 1 , β 1 , γ 11 , γ 12 satisfy the following:
{ ( a ) γ 11 e p 2 V p ( t , e p ) γ 12 e p 2 ( b ) V p t ( t , e p ) + V p e p ( t ,   e p ) e ˙ p α 1 e p 2 ( c ) V p e p ( t , e p ) β 1 e p .
Theorem 1.
It states that if Assumption 1 and Assumption 5 are satisfied by properly choosing a relatively high-gain tuner parameter θ such that the following conditions are met: (1) if ξ 2 ξ ^ 2 converges to 0, then make θ > η 1 α 1 , and (2) if ξ 2 ξ ^ 2 is bounded by e ˜ p , then choose a value for θ such that ( θ * α 1 η 1 η 2 θ * ) 0 .
Then, the system in Equation (11) becomes a converging observer for the system described in Equation (10), which is a transformed form of the process subsystem described in Equation (1).
The proof is given in Appendix A.

3.1.3. State Estimator Design for the New Actuator Subsystem

Equation (18) can be viewed as a transformed form of the actuator subsystem described in Equation (2):
Σ a : { ξ ˙ 2 = f ¯ 2 ( ξ 2 , u ) y 2 = h 2 ( ξ 2 , u )
where f ¯ 2 ( ξ 2 , u ) : = f a ( ξ 2 , u ) .
A convergence observer can be designed for the system in Equation (18) as follows:
Σ aO : { ξ ^ ˙ 2 = f ¯ 2 ( ξ ^ 2 , u ) + κ 2 ( 2 , ξ ^ 2 ) ( h 2 ( ξ ^ 2 ) y 2 ) ˙ 2 = Θ 2 ( ξ ^ 2 , u , 2 )
where κ 2 and Θ 2 are smooth gain functions, with respect to their arguments, the state variable ( 2 , ξ ^ 2 ) belongs to ( n × Θ 2 ) , and Θ 2 is a subset of n , which is positively invariant by the second equation of (19).
The state estimation error is defined as
e a : = ξ ^ 2 ( t ) ξ 2 ( t )
Then, by subtracting the corresponding Equations (18) and (19), we get the following error dynamics:
e ˙ a ( t , e a ) = f ¯ 2 ( ξ 2 , u ) f ¯ 2 ( ξ ^ 2 , u ) κ 2 ( 2 , ξ ^ 2 ) ( h 2 ( ξ ^ 2 ) y 2 )
where Κ 2 ( u , ξ ^ 2 , y 2 ) : = κ 2 ( 2 , ξ ^ 2 ) ( h 2 ( ξ ^ 2 ) y 2 ) .
In order to formulate a solution to the convergence of the above observer, we need to follow Assumption 6, with respect to the error Lyapunov function introduced in [26]. This error Lyapunov function shows the equivalence of the existence of an error Lyapunov function and the existence of a converging observer.
Assumption 6.
Is that for any u U ( t , e a ) A . C ( + , ) , there exists a continuously differentiable function V a and positive constants α 2 ,   β 2 ,   γ 21 ,   γ 22 to satisfy
{ ( a ) γ 21 e a 2 V a ( t , e a ) γ 22 e a 2 ( b ) V a t ( t , e a ) + V a e a ( t , e a ) e ˙ a α 2 e a 2 ( c ) V a e a ( t , e a ) β 2 e a .
The observer defined by Equation (19) is a converging observer if Assumption 6 is satisfied. However, the observer in Equation (19) can only be realized when the output y 2 is measurable, which is not the case. Due to this, y 2 in our design represents the output of the actuator subsystem. It is assumed to be unmeasured and, therefore, y 2 must be replaced by an estimate y ˜ 2 using the available measurements.
Fortunately, that estimation of y ˜ 2 is available in the process subsystem observer in Equation (11). By substituting u a for y ˜ 2 , we can now implement the observer in Equation (19) for an actuator subsystem as
Σ aO : { ξ ^ ˙ 2 = f ¯ 2 ( ξ ^ 2 , u ) + κ 2 ( 2 , ξ ^ 2 ) ( h 2 ( ξ ^ 2 ) y ˜ 2 ) ˙ a = Θ 2 ( ξ ^ 2 , u , 2 )
where Κ 2 ( u , ξ ^ 2 , y ˜ 2 ) : = κ 2 ( 2 , ξ ^ 2 ) ( h 2 ( ξ ^ 2 ) y ˜ 2 ) .
The estimation error is produced again by subtracting the corresponding equation in Equations (18) and (23), and the new error dynamics are achieved as follows:
e ˜ ˙ a ( t , e ˜ a ) = f ¯ 2 ( ξ 2 , u ) f ¯ 2 ( ξ ^ 2 , u ) Κ 2 ( u , ξ ^ 2 , y ˜ 2 ) = f ¯ 2 ( ξ 2 , u ) f ¯ 2 ( ξ ^ 2 , u ) Κ 2 ( u , ξ ^ 2 , y 2 ) + Κ 2 ( u , ξ ^ 2 , y 2 ) Κ 2 ( u , ξ ^ 2 , y ˜ 2 ) = e ˙ a ( t , e a ) + Κ 2 ( u , ξ ^ 2 , y 2 ) Κ 2 ( u , ξ ^ 2 , y ˜ 2 )
In order to ensure the stability of the error dynamics in Equation (24), an assumption is required with respect to the sensitivity of Κ 2 ( u , ξ ^ 2 , y 2 ) , with changes of y 2 .
Assumption 7.
It provides a sufficient condition for achieving this purpose, stating that for any u U , ( t , ξ ^ 2 , y ˜ 2 ) A . C ( + , ) , there exists a real constant γ 2 to satisfy
Κ 2 ( u , ξ ^ 2 , y ˜ 2 ) Κ 2 ( u , ξ ^ 2 , y 2 ) γ 2 y 2 y ˜ 2 .
Similar to Assumption 1, Assumption 7 implies a global Lipchitz-type condition on function Κ 2 such that, in a physical problem, u , y 2 are bounded. Therefore, it can also be replaced by a local smoothness condition.
In addition to asking that the state estimation error e a converge to 0 in the absence of disturbances, we want it to still converge to 0 if a disturbance is present, but converge to 0 and remain bounded if the disturbance is bounded. Therefore, Assumption 7 implies that the definition of e ˜ ˙ a ( t , e ˜ a ) in (24) is not affected.
In particular, since output y ˜ 2 , used in the observer in Equation (23), is in fact a virtual measurement which is estimated by the output of the process subsystem, an estimation error becomes unavoidable. This estimation error can be viewed as a bounded disturbance to the real output of the actuator subsystem y 2 . Therefore, the basic problem addressed in this work is the design of nonlinear observers that possesses robustness to the disturbance affecting the real output.
Theorem 2.
It says that if Assumptions 6 and 7 are satisfied, then the observer described in Equation (23) is a converging observer for the actuator subsystem described in Equation (18).
The proof is given in Appendix B.

3.1.4. Interconnected Observer

The interconnected observer for the studied interconnected system by the system in Equations (11) and (23) is constituted as follows:
Σ O : { ξ ^ ˙ 1 = G ( ξ ^ 11 ) ξ ^ 1 + F ( ξ ^ 11 ) + ε ¯ ( u , u ˙ , ξ ^ 2 ) Λ 1 ( ξ ^ 1 ) S θ 1 C T ( C ξ ^ 1 y 1 ) ξ ^ ˙ 2 = f ¯ 2 ( ξ ^ 2 , u ) + κ 2 ( 2 , ξ ^ 2 ) ( h 2 ( ξ ^ 2 ) ξ ^ 12 )
where the virtual measurement y ˜ 2 in Equation (23) is replaced by the estimation ξ ^ 12 . The observer estimation errors satisfy the following equation:
e p ( t ) = ξ ^ 1 ( t ) ξ 1 ( t ) , e ˜ a ( t ) = ξ ^ 2 ( t ) ξ 2 ( t )

3.2. Interconnected Observer Analysis

The observer in Equation (26) has been designed so that the dynamics of the corresponding error system in Equation (27) are governed as follows:
{ e ˙ p ( t , e p ) = G ( ξ ^ 11 ) ξ ^ 1 + F ( ξ ^ 11 ) + ε ¯ ( u , u ˙ , ξ ^ 2 ) + H ( ξ ^ 1 ) ( C ξ ^ 1 y ) G ( ξ 11 ) ξ 1 F ( ξ 11 ) ε ¯ ( u , u ˙ , ξ 2 ) e ˜ ˙ a ( t , e ˜ a ) = e ˙ a ( t , e a ) + Κ 2 ( u , ξ ^ 2 , h 2 ( ξ ^ 2 ) ) Κ 2 ( u , ξ ^ 2 , ξ ^ 12 )
To analyze the system in Equation (28), our purpose is to study the stability of the error dynamics.
Theorem 3.
If Assumptions 1–7 are satisfied, then a relatively high value of θ can be chosen such that 1 1 2 α 2 ( θ η 1 ) ( η 3 + η 2 θ ) , and the error dynamics governed in Equation (28) are convergent.
Proof. 
The objective is to analyze the stability of the error dynamics. To achieve this purpose, by using V a and V p , defined in the proof of Theorems 1 and 2, the following Lyapunov function candidate is constructed:
V ( t , e ˜ p , e ˜ a ) = V p ( t , e ˜ p ) + V a ( t , e ˜ a )
Then, the time derivation of V ( t , e ˜ p , e ˜ a ) yields
V ˙ ( t , e ˜ p , e ˜ a ) = V p t ( t , e ˜ p ) + V p e p ( t , e ˜ p ) e ˜ ˙ p ( t , e ˜ p ) term 1 + V a t ( t , e ˜ a ) + V a e a ( t , e a ) e ˙ a ( t , e a ) + V a e a ( t , e a ) ( Κ ( u , ξ ^ 2 , y 2 ) Κ ( u , ξ ^ 2 , ξ ^ 12 ) ) term 2
Let us analyze the different terms on the right side of Equation (30), starting with term 1 and using results in the proof of Theorem 1:
V ˙ p ( t , e ˜ p ) = V p t ( t , e ˜ p ) + V p e p ( t , e ˜ p ) e ˜ ˙ p ( t , e ˜ p ) ( θ + η 1 ) V p + η 2 θ V p ξ 2 ξ ^ 2
In turn, by using results in the proof of Theorem 2, then term 2 on the right side of Equation (30) develops as follows:
V ˙ a ( e ˜ a ) = V a t ( t , e ˜ a ) + V a e a ( t , e a ) e ˙ a ( t , e a ) + V a e a ( t , e a ) ( Κ 2 ( u , ξ ^ 2 , y 2 ) Κ 2 ( u , ξ ^ 2 , ξ ^ 12 ) ) α 2 V a + η 3 V a ξ 12 ξ ^ 12
Then, the overall inequality yields
V ˙ ( t , e ˜ p , e ˜ a ) α 2 V a + η 3 V p V a + ( θ + η 1 ) V p + η 2 θ V p V a
Now, set V a * = α 2 V a , V p * = ( θ η 1 ) V p , and V * = V a * + V p * .
Let us assume that
ς = min { α 2 , ( θ η 1 ) }
In this case,
V * 2 ς V
It should also be noted that
V p * + V a * 2 V p * V a * = 2 α 2 ( θ η 1 ) V p V a
Thus,
V p V a 1 2 α 2 ( θ η 1 ) V *
It is easy to get that inequality in Equation (33) to yield the following:
V ˙ ( e ˜ p , e ˜ a ) V * + 1 2 α 2 ( θ η 1 ) ( η 3 + η 2 θ ) V * ( 1 1 2 α 2 ( θ η 1 ) ( η 3 + η 2 θ ) ) V * 2 ς ( 1 1 2 α 2 ( θ η 1 ) ( η 3 + η 2 θ ) ) V
Now, it suffices to choose a value of θ such that ( 1 1 2 α 2 ( θ η 1 ) ( η 3 + η 2 θ ) ) 0 .
This ends the proof. □

4. Simulations

Numerical simulations were performed to validate that the interconnected observer given by Equation (26) can be implemented for monitoring the performance of an interconnected system. A case study was developed on an intensified heat exchanger (IHEX). The pilot consisted of three process plates sandwiched between five utility plates. Two pneumatic control valves were used to control the utility and process fluid. More relative information could be found in [17]. Moreover, the outlet fluid flow rates of the control valves were assumed to be unmeasured to ensure a realistic simulation. Therefore, during the course of the simulation work, the proposed observers were designed for estimating unmeasured inlet fluid flows and monitoring the performance of the IHEX.

4.1. System Modelling

4.1.1. Actuator Subsystem Modelling

The pneumatic control valve is used to act as an actuator in this system. By applying Bernoulli’s continuous flow law of incompressible fluids, we have
F = C v f ( X ) Δ P sg
where F is the flow rate ( m 3 s 1 ), Δ P is the fluid pressure drop across the valve (Pa), sg is the specific gravity of the fluid and equals 1 for pure water, X is the valve opening percentage, C v is the valve coefficient, and f(X) is the flow characteristic which is defined as the relationship between the valve capacity and fluid traveling through the valve. In [3], a pneumatic control valve had a dynamic model as follows:
m d 2 X dt + μ dX dt + kX = p c A a
dp c dt = ( p c + p a ) ( C V P A a V 0 + A a X dX dt )
where A a is the diaphragm area on which the pneumatic pressure acts, p c is the pneumatic pressure, m is the mass of the control valve stem, μ is the friction of the valve stem, k is the spring compliance, X is the stem displacement or percentage opening of the valve, C V P is the command pressure, and p a is the air pressure. The following definitions also apply:
x a T = [ x a 1 x a 2 x a 3 x a 4 x a 5 x a 6 ] = [ X 1 dX 1 dt p c 1 X 2 dX 2 dt p c 2 ]
u T = [ u 1 u 2 ] = [ C V P 1 C V P 2 ] , u a T = [ F 1 F 2 ] = [ C v Δ P 1 sg X 1 C v Δ P 2 sg X 2 ]
C = [ c 1 c 2 c 3 c 4 c 5 c 6 ] = [ C v Δ P 1 sg 0 0 C v Δ P 2 sg 0 0 ]
where p c 1 , p c 2 , X 1 , X 2 , C V P 1 , C V P 2 , Δ P 1 , and Δ P 2 correspond to p c , X , C V P , and Δ P in Equations (36)–(38), respectively, and subscripts one and two represent two different control valves. The actuator subsystem is then described as
{ x ˙ a 1 = x a 2 x ˙ a 2 = k 1 m x a 1 μ 1 m x a 2 + A a m x a 3 x ˙ a 3 = x a 3 u 1 A a V 0 + A a x a 1 x a 3 x a 2 + p a u 1 A a p a V 0 + A a x a 1 x a 2 x ˙ a 4 = x a 5 x ˙ a 5 = k 2 m x a 4 μ 2 m x a 5 + A a m x a 6 x ˙ a 6 = x a 6 u 2 A a V 0 + A a x a 4 x a 5 x a 6 + p a u 2 A a p a V 0 + A a x a 4 x a 5 y = C x a

4.1.2. Process Subsystem Modelling

The IHEX can be modeled based on the mass and energy balances, which describe the evolution of characteristic values such as temperature, mass, composition, and pressure. Considering the heat exchanger system taken from [17], the dynamic equation governing the heat balance of the fluids is given by
T ˙ p = UA ρ p c pp V p ( T u T p ) + 1 V p ( T pi T p ) F p
T ˙ u = UA ρ u c pu V u ( T p T u ) + 1 V u ( T ui T u ) F u
where ρ p , ρ u are the densities of the fluids (in kg · m 3 ), V p , V u are the volumes of the fluids (in m 3 ), c pp , c pu are the specific heats of the fluids (in J · kg 1 · K 1 ), U is the overall heat transfer coefficient (in J · m 2 · K 1 · s 1 ), A is the reaction area (in m 2 ), F p , F u are the mass flow rates of the fluids (in kg · s 1 ), and T pi , T ui are the inlet temperatures of the fluids.
If the state vector is defined as x T = [ x 1 , x 2 ] T = [ T p , T u ] T , the control input as u a T = [ u a 1 , u a 2 ] T = [ F p , F u ] T , and the output vector as y T = [ y 1 , y 2 ] T = [ T p , T u ] T , then the above two equations can be rewritten as
{ x ˙ = f ( x ) + i = 1 2 g i ( x ) u a y = h ( x , u a )
where f ( x ) = ( f 1 ( x ) f 2 ( x ) ) = ( h p A ρ p C p p V p ( T p T u ) h u A ρ u C p u V u ( T u T p ) ) , g ( x ) = ( g 1 , g 2 ) = ( ( T pi T p ) V p 0 0 ( T ui T u ) V u ) , and y 1 = x 1 , y 2 = x 2 .
By using Equation (34), a function for the derivatives for u a is obtained:
u ˙ a = ε ( u , u ˙ , x a ) = h a u ( u , x a ) u ˙ + h a x a ( u , x a ) f a ( u , x a ) = ( C v Δ P 1 sg 0 C v Δ P 2 sg 0 ) x a + ( A a m C v Δ P 1 sg A a m C v Δ P 2 sg ) u
If the state vector is defined as x 1 T = [ x 11 , x 12 ] T = [ T p , T u ] T , the unmeasured state as x 2 T = [ x 21 , x 22 ] T = [ u a 1 , u a 2 ] T = [ F p , F u ] T , and the output vector as y T = [ y 1 , y 2 ] T = [ T p , T u ] T , then Equations (44) and (45) can be rewritten as
{ x ˙ 1 = G 1 ( x 1 ) x 2 + g 1 ( x 1 , u ) x ˙ 2 = ε ( u , u ˙ , x a ) y = x 1
where G 1 ( x 1 ) = ( ( T pi x 11 ) V p 0 0 ( T ui x 12 ) V u ) and f 1 ( x ) = ( h p A ρ p C p p V p ( x 11 x 12 ) h u A ρ u C p u V u ( x 12 x 11 ) ) .

4.2. Observer Design

4.2.1. Observer 1 for the Actuator Subsystem

In this model, outputs were considered as unmeasured and were substituted by its estimation proposed in Observer 2, then an extended high-gain observer of the form in Equation (26) for the system in Equation (41) is given by
{ x ^ ˙ a 1 = x ^ a 2 k 1 ( C v Δ P 1 sg x ^ a 1 x ^ 21 ) x ^ ˙ a 2 = k 1 m x ^ a 1 μ 1 m x ^ a 2 + A a m x ^ a 3 k 2 ( C v Δ P 1 sg x ^ a 1 x ^ 21 ) x ^ ˙ a 3 = x ^ a 3 u 1 A a V 0 + A a x ^ a 1 x ^ a 3 x ^ a 2 + p a u 1 A a p a V 0 + A a x ^ a 1 x ^ a 2 k 3 ( C v Δ P 1 sg x ^ a 1 x ^ 21 ) x ^ ˙ a 4 = x ^ a 5 k 4 ( C v Δ P 1 sg x ^ a 4 x ^ 22 ) x ^ ˙ a 5 = k 2 m x ^ a 4 μ 2 m x ^ a 5 + A a m x ^ a 6 k 5 ( C v Δ P 1 sg x ^ a 4 x ^ 22 ) x ^ ˙ a 6 = x ^ a 6 u 2 A a V 0 + A a x ^ a 4 x ^ a 5 x ^ a 6 + p a u 2 A a p a V 0 + A a x ^ a 4 x ^ a 5 k 6 ( C v Δ P 1 sg x ^ a 4 x ^ 22 )

4.2.2. Observer 2 for the Process Subsystem

It should be noted that the original system in Equation (45) has been augmented with the differential equation u ˙ a = ε ( u , u ˙ , x a ) ; that is to say, the unknown inputs are treated like an unmeasured state. Then, it is possible to design an observer of the form in Equation (26) for the system in Equation (46) as follows:
{ x ^ ˙ 1 = ( ( T pi x ^ 11 ) V p 0 0 ( T ui x ^ 12 ) V u ) x ^ 2 + ( h p A ρ p C p p V p ( x ^ 11 x ^ 12 ) h u A ρ u C p u V u ( x ^ 12 x ^ 11 ) ) ( 2 θ 2 θ ) ( y ^ y ) x ^ ˙ 2 = ( C v Δ P 1 sg 0 C v Δ P 2 sg 0 ) x ^ a + ( A a m C v Δ P 1 sg A a m C v Δ P 2 sg ) u ( θ 2 h p A ρ p C p p V p ( x ^ 11 x ^ 12 ) θ 2 h u A ρ u C p u V u ( x ^ 12 x ^ 11 ) ) ( y ^ y ) y ^ = x ^ 1

4.3. Numerical Simulations

In order to test the performance of the proposed observers, two numerical simulations were carried out. Considering the actuator and process models given by Equations (44) and (46), Observer 1 in Equation (47) and Observer 2 in Equation (48) were designed for estimating unmeasured inlet flows F p , F u and monitoring the performance of the final products T p , T u . Available measurements involve the inlet–outlet temperatures of the hex reactor and the pneumatic pressure p c 1 , p c 2 of the actuators.
Two cases were considered. In Case 1, constant inlet flows F p , F u were considered in both fluids. By contrast, in Case 2, F p , F u were considered to be time varying. Observers 1 and 2 were simulated with respect to the actuator and process subsystems, using the values corresponding to an IHEX system from [17]. The parameters in the process subsystem were as follows: hA = 214.8   W · K 1 , V p = 2.685 × 10 5   m 3 ,   V u = 1.141 × 10 4   m 3 , and the inlet temperatures T pi and T ui were 76 °C and 15 °C, respectively. The parameters in the actuator subsystem were as follows: m = 2 kg, A a = 0.029 m2, μ = 1500 Ns/m, k = 6089 Ns/m, Pc for the utility fluid was 1 MPa (1.2 Mpa for the process fluid), and the pressure drop Δ P in the utility fluid was 0.6 MPa (60 KPa in the process fluid).
In order to illustrate the robustness, external disturbances and measurement noise were considered. Suppose the output measurement is corrupted by a colored noise. In addition to noise, no error is assumed in measuring T p   and   T u . Moreover, to compare the effectiveness of the proposed method with other existing ones, we employed a typical unknown input observer (UIO) in [31] to show the differences.

4.3.1. Case 1: Both Fluid Flow Rates F p , F u Are Constant

The objective of this simulation was to prove the convergence of the observers in a common situation where both fluid flow rates remained constant over a long time. The computed inlet flow rate of the utility fluid F u was 4.22 × 10 5   m 3 s 1 , and the inlet flow rate of the process fluid F p was a constant 4.17 × 10 6   m 3 s 1 . The computed value meant the expected true values of the actuators. The initial conditions of the process model were T p 0 = 80 °C and T u 0 = 20 °C, respectively, and for the observers they were T ^ p 0 = T ^ u 0 = 30 . The discrepancies between the initial conditions of the process and those of the observers were reasonable and realistic, considering that the temperature was a process variable that could be easily measured. In order to evaluate the observer performance against uncertainties of the knowledge of the fluid flow rate, the initial value of the estimates were F ^ p 0 = F ^ u 0 = 0 in both observers. This assumption represented a relatively rough situation in the practical engineering world. However, simulation results showed encouraging results. The tuning parameters were k 1 = k 3 = 100 , k 2 = k 4 = 0.15 (for Observer 1) and θ = 80 (for Observer 2).
The results are reported in Figure 3 and Figure 4. In Figure 3, the dashed curves correspond to the estimates using Observer 1, and the solid lines are the measured temperatures. It can be seen that, whether noise-free or noise-corrupted, the convergence of the estimated T ^ p and T ^ u values proved to be fast (in several seconds). It is not surprising because, actually, T p and T u were the measured outputs of the overall system.
The main contribution of the proposed method is the capacity that, besides states, it can estimate the unknown connection of an interconnected system by the final measured outputs, where the unknown connection represents the inlet fluid flow rate F p , F u in the IHEX system. Figure 5 and Figure 6 give the encouraging results. Figure 5 shows the computation and estimation of the process fluid flow rate F p , using a UIO and Observer 2. Figure 6 is the results for the utility fluid flow rate F u . As expected, F ^ p and F ^ u follow different trajectories before they converge toward the true value (computed value) in a relatively short transient period, whether or not noise exists. However, the convergence speed of the proposed method is obviously faster than that of the UIO. In addition, the proposed observer was more robust with noise. For F p , as shown in Figure 5, if no noise was present, after less than 5 s, the three curves overlapped. Compared with the curve of the UIO on the dashed, dotted line, the curve of the proposed Observer 2 in dashed, dotted line only needed less than 1 s to track the simulated curve of the solid line. Moreover, from Figure 5, we can see that the curve of the UIO in the dashed line was significantly more affected by noise than that of Observer 2 in the dashed, dotted line. The similar results are shown in Figure 6 with respect to F u . According to Figure 6, it took about 1.5 s for the three curves to be overlapped, while for the curve of Observer 2 on the dashed line, less than 0.5 s was taken. The impact of noise is relatively obvious on the curve of the UIO, while for the curve of Observer 2, the influence was less significant.
In both fluids, these differences were caused by their varied initial values. However, it is proven that if adequate values of the tuning parameters k and θ are selected, no matter the degree of deviation of the initial value of F ^ p 0 and F ^ u 0 from the simulated values in the system model, convergences are guaranteed. Larger values of these tuning parameters ensure a smaller convergence time, while smaller values have the opposite effect. However, large tuning values should be avoided, since the observer may become too sensitive to measurement noise in real-time applications. According to the above analysis, we can readily conclude that the proposed interconnected observer works effectively and robustly for the designed purposes with proper tuning parameters.

4.3.2. Case 2: The Fluid Flow Rates Vary Due to Parameter Changes

The present computations were executed to get an accurate screening of the variation of the observer estimate by corroborating if they were in agreement with the simulated fluid flow rates which underwent either an abrupt change or a gradual variation due to, for example, aging or erosion.
The parameter effects were taken into account in the following way. Initial values of F u = 4.22 × 10 5   m 3 s 1 and F p = 4.17 × 10 6   m 3 s 1 are considered, followed by an abrupt change of F u at 80 s. The reason for this change is the variation of parameter Δ P , with the value changing from 0.6 MPa to 0.4 Mpa. Several factors can be attributed to this kind of variation, such as valve clogging or an unexpected pressure drop across the control valves. After that, at t = 150 s, F p began to deteriorate due to an increase of the spring compliance k 2 in the process fluid actuator. One main reason contributing to this change is erosion. Because of erosion, the gland packing of the valve may loosen, which leads to stem vibration. In the simulation, a value of 1000 nm 1 was added to the spring compliance k 2 . This simulation was carried out using the same constants used in the previous simulation and the same values of T pi and T ui , as well as T p 0 and T u 0 , were used. The initial conditions of both observers, as well as the observer parameters ( k 1 , k 2 , k 3 , k 4 , and θ ) were the same as the previous ones. These variations are illustrated in Figure 7, Figure 8, Figure 9 and Figure 10.
As shown in Figure 7 and Figure 8, whether in a noise-free or noise-corrupted situation, the estimated outlet fluid temperatures T ^ p and T ^ u on the dashed line can better track the curve of the measurements T p and T u on the solid line after a short transient time. At 80 s, both curves decreased unexpectedly and finally stabilized at a new level. A drop of 0.2 °C was observed. Shortly after, at t = 150 s, another drop happened, and a new stable level was expected with a 0.9 °C reduction. These decreases imply the influences of parameter changes in the fluid actuators, and no further variations illustrate the occurrence of additional changes. A similar result was obtained in the estimated T ^ u of the utility fluid in Figure 8. It is shown that, due to changes of Δ P 1 and k 2 , the measured T u dropped 0.2 °C and 0.5 °C at 80 s and 150 s, respectively. The estimated T ^ u on the dashed line tracks T u after the observer converges.
The simulation curves indicate that the proposed observer was effective at tracking system performance. However, it is shown that, because of noise, the drops caused by local parameter changes at 80 s could not be visually observed on the global output directly. In order to analyze the influence of these changes, we have to detect it through the local dynamics supervision. Therefore, it is very meaningful to monitor both the local and global dynamics of a control system.
It can be seen from Figure 9 that, in a noise-free situation, two estimate curves on the dashed line and the dashed, dotted line converge to the simulated value F p of the solid line quickly after a transient response. The dashed line is the estimated process fluid flow rate F ^ p produced by the UIO, and the dashed, dotted line was generated by Observer 2. Later, at 150 s, there was an unexpected drop by the simulated F p on the solid line. Fortunately, both estimate curves responded quickly to this variation, and the curve of the UIO took 1.5 s to track F p again, while the curve of Observer 2 converged more quickly than that of the UIO. The decrease implies parameter changes in the process fluid actuator, which satisfy the assumption that k 2 changes at t = 150 s. When it comes to the noise-corrupted situation, the curve of the UIO was significantly affected by the noise, although a drop at 150 s can still be observed. Luckily, the curve of Observer 2 was relatively robust to the noise, as it converged to the simulated value of the solid line with a minor difference compared with the noise-free situation.
Figure 10 demonstrates the results for the utility fluid. The same results were obtained in a noise-free situation at a time of 80 s. As expected, the simulated F u of the solid line jumped due to the change of Δ P 1 . The estimated value of F ^ u of the dashed line by the UIO and the dashed, dotted line by Observer 2 track well with the value of F u on the solid line. Just as with the process fluid situation, the convergence speed of Observer 2 was faster than that of the UIO. If measurements were corrupted by noise, the curve of Observer 2 was not obviously influenced. Conversely, the curve of the UIO was greatly impacted.
To sum up, obviously, the proposed interconnected observer is effective, even when the unknown connection is varying in time simultaneously. Therefore, the proposed observer proves its ability to monitor the performance of interconnected systems and to estimate unknown interconnections.

5. Conclusions

The goal of the design methodology presented in this paper was to enable or simplify observer design for systems that are otherwise difficult to handle by allowing the designer to focus on a smaller, nonlinear subsystem. That is to say, we mainly focused on observing, for example, how the change of an internal parameter at the local level affects the global output at the global level.
An interconnected observer is designed to estimate both the state and unmeasured interconnection at the local and global levels. As a result, both local and global dynamics can be observed, as well as the influence of local dynamics on global dynamics. In particular, the interconnection is not supposed to be accessible to measurement. In order to achieve this purpose, firstly, an existing observer is extended to estimate the states of the actuator subsystem. Particularly, the information of the actuator subsystem output is substituted by their estimates, achieved by the observer of the process subsystem. Secondly, a kind of extended, high-gain observer is produced to estimate the states of the process subsystem, which is subjected to a precise unknown input. The unknown input is considered a new state of the process subsystem, and it is expressed as a function of the inputs, derivatives of the inputs, and the states of the actuator subsystem. Thus, an interconnected observer is proposed by using the estimates of the states and the unmeasured interconnection, and the convergence is investigated. Finally, satisfactory simulation results are obtained to confirm the effectiveness and robustness of the proposed method.
In this paper, it is clear that the physical motivation for the decomposition of a control system into the actuator and the process parts is physically motivated. In this respect, for control analysis purposes, the condition for decomposition of an independent control system may be the target of future investigations, like the inverted pendulum on a cart. Another open question worth addressing is the demonstration of stability and sensitivity of the estimation error, like the use of ISS to investigate the stability of the estimation error in [32].

Author Contributions

M.Z. and Z.L. conceived and designed the study. M.Z. and Q.W. performed the observer design, carried out simulations, and wrote the original draft. B.D. and X.C. reviewed and edited the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, Grant No. 62003106 and No. 51867006, Talent Project of GZU (2018) 02. Key Lab construction project of Guizhou Province (2016) 5103.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
ISSinput-to-state stability
ISDSinput-to-state dynamical stability
UIOunknown input observer
IHEXintensified heat exchanger

Appendix A

Proof of Theorem 1.
Before preceding the convergence proof of the observer, one introduces the following notations:
S θ = 1 θ Δ θ S 1 Δ θ
where
S 1 = S θ | θ = 1   and   Δ θ = [ I n 0 0 1 θ I n ]
To proof Theorem 1, the estimation error is introduced as e p ( t ) = ξ ^ 1 ( t ) ξ 1 ( t ) .
Then, subtracting corresponding Equations (10) and (11), one gets the following error dynamics:
e ˙ p = G ( ξ ^ 11 ) ξ ^ 1 + F ( ξ ^ 11 ) + ε ¯ ( u , u ˙ , ξ ^ 2 ) + Λ 1 ( ξ ^ 1 ) S θ 1 C T ( C ξ ^ 1 y ) G ( ξ 11 ) ξ 1 F ( ξ 11 ) ε ¯ ( u , u ˙ , ξ 2 )
= ( G ( ξ ^ 11 ) + Λ 1 ( ξ ^ 1 ) S θ 1 C T C ) e p + ( G ( ξ ^ 11 ) G ( ξ 11 ) ) ξ 1 + F ( ξ ^ 11 ) F ( ξ 11 ) + e ε ¯ ( t , e ˜ a )
where e ε ¯ ( t , e ˜ a ) = ε ¯ ( u , u ˙ , ξ ^ 2 ) ε ¯ ( u , u ˙ , ξ 2 ) .
By setting e ˜ p = Δ θ Λ ( ξ ^ 11 ) e p , one can then get
e ˜ ˙ p ( t ) = θ ( A S 1 1 C T C ) e ˜ p + Λ ˙ ( ξ ^ 11 ) Λ 1 ( ξ ^ 11 ) e ˜ p + Δ θ Λ ( ξ ^ 11 ) ( F ( ξ ^ 11 ) F ( ξ 11 ) ) + Δ θ Λ ( ξ ^ 11 ) G ( ξ ^ 11 ) ξ 1 G ( ξ 11 ) ξ 1 Δ θ Λ ( ξ ^ 11 ) e ε ¯ ( t , e ˜ a )
To analyze the dynamics of the error system, the following positive Lyapunov function candidate is considered:
V p ( t , e ˜ p ) = e ˜ p T S 1 e ˜ p
Convergence of the observer is described by the time derivation of V p ( t , e ˜ p ) . Then, we obtain
V ˙ p ( t , e ˜ p ) = 2 e ˜ p T S 1 e ˜ ˙ p = θ ( 2 e ˜ p T S 1 A e ˜ p 2 e ˜ p T C T C e ˜ p ) + 2 e ˜ p T S 1 Λ ˙ ( ξ ^ 11 ) Λ 1 ( ξ ^ 11 ) e ˜ p + 2 e ˜ p T S 1 Δ θ Λ ( ξ ^ 11 ) ( F ( ξ ^ 11 ) F ( ξ 11 ) ) + 2 e ˜ p T S 1 Δ θ Λ ( ξ ^ 11 ) ( ( G ( ξ ^ 11 ) G ( ξ 11 ) ) ξ 1 ) 2 e ˜ p T S 1 Δ θ Λ ( ξ ^ 11 ) e ε ¯ ( t , e ˜ a ) = θ V p θ C e ˜ p 2 + 2 e ˜ p T S 1 Λ ˙ ( ξ ^ 11 ) Λ 1 ( ξ ^ 11 ) e ˜ p + 2 e ˜ p T S 1 Δ θ Λ ( ξ ^ 11 ) ( F ( ξ ^ 11 ) F ( ξ 11 ) ) + 2 e ˜ p T S 1 Δ θ Λ ( ξ ^ 11 ) ( ( G ( ξ ^ 11 ) G ( ξ 11 ) ) ξ 1 ) 2 e ˜ p T S 1 Δ θ Λ ( ξ ^ 11 ) e ε ¯ ( t , e ˜ a ) θ V p + 2 S 1 e ˜ p Λ ˙ ( ξ ^ 11 ) Λ 1 ( ξ ^ 11 ) e ˜ p + 2 S 1 e ˜ p Λ ( ξ ^ 11 ) Δ θ ( F ( ξ ^ 11 ) F ( ξ 11 ) ) + 2 S 1 e ˜ p Δ θ Λ ( ξ ^ 11 ) ( ( G ( ξ ^ 11 ) G ( ξ 11 ) ) ξ 1 ) + 2 θ S 1 e ˜ p G ( ξ ^ 11 ) e ε ¯ ( t , e ˜ a ) θ V p + 2 μ S 1 e ˜ p e ˜ p + 2 ρ S 1 e ˜ p Δ θ ( F ( ξ ^ 11 ) F ( ξ 11 ) ) + 2 S 1 e ˜ p Δ θ Λ ( ξ ^ 11 ) ( ( G ( ξ ^ 11 ) G ( ξ 11 ) ) ξ 1 ) + 2 τ γ 1 θ S 1 e ˜ p ξ 2 ξ ^ 2
where μ = sup t 0 Λ ˙ ( ξ ^ 11 ) Λ 1 ( ξ ^ 11 ) , ρ is the upper boundary of Λ ( ξ ^ 11 ) , e ε ¯ ( t , e ˜ a ) γ 1 ξ 2 ξ ^ 2 , as proposed in Assumption 1, and τ is a finite real number with 0 < ϱ τ , such that ϱ 2 I n F T ( ξ 11 ) F ( ξ 11 ) τ 2 I n , as given in Assumptions 2–4.
Then, we have
Δ θ ( F ( ξ ^ 11 ) F ( ξ 11 ) ) f ( ξ ^ 11 ) f ( ξ 11 ) σ e ˜ p
where σ denotes the Lipchitz constants of f ( ξ 11 ) .
Similarly, we have
Δ θ Λ ( ξ ^ 11 ) ( ( G ( ξ ^ 11 ) G ( ξ 11 ) ) ξ 1 ) = ( ( g ( ξ ^ 11 ) g ( ξ 11 ) ) ξ 12 ) ϵ e ˜ p
where the positive constant ϵ denotes the boundary of ξ 12 .
Thus, we get the following:
V ˙ p ( t , e ˜ p ) θ V p + η 1 V p + η 2 θ V p ξ 2 ξ ^ 2
( θ + η 1 ) V p + η 2 θ V p ξ 2 ξ ^ 2
where η 1 = ( 2 μ + 2 ρ σ + 2 ϵ ) ξ ( S 1 ) and ξ ( S 1 ) = λ max ( S 1 ) / λ min ( S 1 ) , η 2 = 2 τ γ 1 λ max ( S 1 ) . λ max ( S 1 ) (resp. λ max ( S 1 ) ) is the largest (resp. the smallest) eigenvalue of λ max ( S 1 ) .
Thus, the following hold true:
  • if ξ 2 ξ ^ 2 converges to 0, it results in V ˙ p ( t , e ˜ p ) ( θ + η 1 ) V p . Then, by taking θ > θ * = η 1 , the negative of the right side of the above inequality is obtained;
  • if ξ 2 ξ ^ 2 is bounded by e ˜ p , it results in V ˙ p ( t , e ˜ p ) ( θ + η 1 + η 2 θ ) ) V p . Then, by choosing θ > θ * such that ( θ * η 1 η 2 θ * ) 0 , the negative of the right side of the above inequality is obtained.
That ends the proof. □

Appendix B

Proof of Theorem 2.
In order to show that the system described in Equation (23) represents a converging observer for the system in Equation (18), we need to make its corresponding error dynamics in Equation (24) coincide with Assumption 6, which has been proven to be a condition for the existence of a converging observer.
e ˜ ˙ a ( t , e ˜ a ) = f ¯ 2 ( ξ 2 , u ) f ¯ 2 ( ξ ^ 2 , u ) Κ 2 ( u , ξ ^ 2 , y ˜ 2 ) = f ¯ 2 ( ξ 2 , u ) f ¯ 2 ( ξ ^ 2 , u ) Κ 2 ( u , ξ ^ 2 , y 2 ) + Κ 2 ( u , ξ ^ 2 , y 2 ) Κ 2 ( u , ξ ^ 2 , y ˜ 2 ) = e ˙ a ( t , e a ) + Κ 2 ( u , ξ ^ 2 , y 2 ) Κ 2 ( u , ξ ^ 2 , y ˜ 2 )
Therefore, by computing the time derivation of V a with respect to the trajectory e ˜ a in Equation (24), using Assumptions 6 and 7, it follows that
V ˙ a ( t , e ˜ a ) = V a t ( t , e a ) + V a e a ( t , e a ) e ˙ a ( t , e a ) + V a e a ( t , e a ) ( Κ 2 ( u , ξ ^ 2 , y 2 ) Κ 2 ( u , ξ ^ 2 , y ˜ 2 ) )
α 2 e a 2 + β 2 γ 2 e a y 2 y ˜ 2
α 2 e a 2 + η 3 e a y 2 y ˜ 2
where η 3 = β 2 γ 2 . Since the output y ˜ 2 used in the observer in Equation (23) is in fact a virtual measurement which is estimated by the output of the process subsystem, it is in fact ξ 12 ,   ξ ^ 12 .
V ˙ a ( t , e ˜ a ) α 2 e a 2 + η 3 e a ξ 12 ξ ^ 12
Thus, the right side of this inequality is negative if the following conditions are met:
  • α 2 β 2 γ 2 and ξ 12 ξ ^ 12 is bounded by e a , resulting in V ˙ a ( t , e ˜ a ) ( α 2 η 3 ) e a 2 ;
  • ξ 12 ξ ^ 12 converges to 0, resulting in V ˙ a ( t , e ˜ a ) α 2 e a 2 .
This ends the proof. □

References

  1. Yang, J.; Zhu, F.; Yu, K.; Bu, X. Observer-based state estimation and unknown input reconstruction for nonlinear complex dynamical systems. Commun. Nonlinear Sci. Numer. Simul. 2015, 20, 927–939. [Google Scholar] [CrossRef]
  2. Gauthier, J.; Hammouri, H.; Othman, S. A simple observer for nonlinear systems applications to bioreactors. IEEE Trans. Autom. Control. 1992, 37, 875–880. [Google Scholar] [CrossRef]
  3. Bartyś, M.; Patton, R.; Syfert, M.; Heras, S.D.L.; Quevedo, J. Introduction to the DAMADICS actuator FDI benchmark study. Control. Eng. Pr. 2006, 14, 577–596. [Google Scholar] [CrossRef]
  4. Moreno, J.A.; Alvarez, J. On the estimation problem of a class of continuous bioreactors with unknown input. J. Process. Control. 2015, 30, 34–49. [Google Scholar] [CrossRef]
  5. De Persis, C.; Isidori, A. A geometric approach to nonlinear fault detection and isolation. IEEE Trans. Autom. Control. 2001, 45, 853–865. [Google Scholar] [CrossRef] [Green Version]
  6. Zhu, F. State estimation and unknown input reconstruction via both reduced-order and high-order sliding mode observers. J. Process. Control. 2012, 22, 296–302. [Google Scholar] [CrossRef]
  7. Besançon, G. Remarks on nonlinear adaptive observer design. Syst. Control. Lett. 2000, 41, 271–280. [Google Scholar] [CrossRef]
  8. Vijay, P.; Tadé, M.O.; Ahmed, K.; Utikar, R.; Pareek, V.K. Simultaneous estimation of states and inputs in a planar solid oxide fuel cell using nonlinear adaptive observer design. J. Power Sources 2014, 248, 1218–1233. [Google Scholar] [CrossRef] [Green Version]
  9. Li, Z.; Dahhou, B. A new fault isolation and identification method for nonlinear dynamic systems: Application to a fermentation process. Appl. Math. Model. 2008, 32, 2806–2830. [Google Scholar] [CrossRef]
  10. Pan, S.; Xiao, D.; Xing, S.; Law, S.; Du, P.; Li, Y. A general extended Kalman filter for simultaneous estimation of system and unknown inputs. Eng. Struct. 2016, 109, 85–98. [Google Scholar] [CrossRef]
  11. Manaa, I.; Barhoumi, N.; M’Sahli, F. Unknown inputs observers design for a class of nonlinear switched systems. Int. J. Model. Identif. Control. 2015, 23, 45. [Google Scholar] [CrossRef]
  12. Besançon, G.; Munteanu, I. Control strategy for state and input observer design. Syst. Control. Lett. 2015, 85, 118–122. [Google Scholar] [CrossRef]
  13. Guo, S.; Zhu, F. Actuator fault detection and interval reconstruction based on interval observers. IFAC-PapersOnLine 2017, 50, 5061–5066. [Google Scholar] [CrossRef]
  14. Buciakowski, M.; Witczak, M.; Puig, V.; Rotondo, D.; Nejjari, F.; Korbicz, J. A bounded-error approach to simultaneous state and actuator fault estimation for a class of nonlinear systems. J. Process. Control. 2017, 52, 14–25. [Google Scholar] [CrossRef] [Green Version]
  15. Vijay, P.; Tadé, M.O. An adaptive non-linear observer for the estimation of temperature distribution in the planar solid oxide fuel cell. J. Process. Control. 2013, 23, 429–443. [Google Scholar] [CrossRef]
  16. Marzat, J.; Piet-Lahanier, H.; Damongeot, F.; Walter, É. Fault diagnosis for nonlinear aircraft based on control-induced redundancy. In Proceedings of the SysTol 2010: Conference on Control and Fault-Tolerant Systems, Nice, France, 6–8 October 2010; pp. 119–124. [Google Scholar]
  17. Théron, F.; Anxionnaz-Minvielle, Z.; Cabassud, M.; Gourdon, C.; Tochon, P. Characterization of the performances of an innovative heat-exchanger/reactor. Chem. Eng. Process. Process. Intensif. 2014, 82, 30–41. [Google Scholar] [CrossRef] [Green Version]
  18. Zhang, M.; Li, Z.-T.; Cabassud, M.; Dahhou, B. Root cause analysis of actuator fault based on invertibility of interconnected system. Int. J. Model. Identif. Control. 2017, 27, 256. [Google Scholar] [CrossRef] [Green Version]
  19. Zou, T.; Wu, S.; Zhang, R. Improved state space model predictive fault-tolerant control for injection molding batch processes with partial actuator faults using GA optimization. ISA Trans. 2018, 73, 147–153. [Google Scholar] [CrossRef]
  20. Kadlec, P.; Gabrys, B.; Strandt, S. Data-driven Soft Sensors in the process industry. Comput. Chem. Eng. 2009, 33, 795–814. [Google Scholar] [CrossRef] [Green Version]
  21. Boccaletti, S.; Latora, V.; Moreno, Y.; Chavez, M.; Hwang, D. Complex networks: Structure and dynamics. Phys. Rep. 2006, 424, 175–308. [Google Scholar] [CrossRef]
  22. Keliris, C.; Polycarpou, M.M.; Parisini, T. A robust nonlinear observer-based approach for distributed fault detection of input–output interconnected systems. Automatica 2015, 53, 408–415. [Google Scholar] [CrossRef] [Green Version]
  23. Gao, N.; Darouach, M.; Alma, M.; Voos, H. Decentralized dynamic-observer-based control for large scale nonlinear uncertain systems. In Proceedings of the 2015 American Control Conference (ACC), Chicago, IL, USA, 1–3 July 2015; pp. 4131–4136. [Google Scholar]
  24. Li, Y.; Sanfelice, R.G. Interconnected Observers for Robust Decentralized Estimation With Performance Guarantees and Optimized Connectivity Graph. IEEE Trans. Control. Netw. Syst. 2015, 3, 1–11. [Google Scholar] [CrossRef] [Green Version]
  25. Chakrabarty, A.; Sundaram, S.; Corless, M.J.; Buzzard, G.T.; Żak, S.H.; Rundell, A.E. Distributed unknown input observers for interconnected nonlinear systems. In Proceedings of the 2016 American Control Conference (ACC), Boston, MA, USA, 6–8 July 2016; pp. 101–106. [Google Scholar]
  26. Besançon, G.; Hammouri, H. On observer design for interconnected systems. J. Math. Syst. Estim. Control 1998, 8, 1–26. [Google Scholar]
  27. Dashkovskiy, S.; Naujok, L. Quasi-ISS/ISDS observers for interconnected systems and applications. Syst. Control. Lett. 2015, 77, 11–21. [Google Scholar] [CrossRef]
  28. Ahmed-Ali, T.; Giri, F.; Krstic, M.; Lamnabhi-Lagarrigue, F. Observer design for a class of nonlinear ODE–PDE cascade systems. Syst. Control. Lett. 2015, 83, 19–27. [Google Scholar] [CrossRef]
  29. Grip, H.F.; Saberi, A.; Johansen, T.A. Observers for interconnected nonlinear and linear systems. Automatica 2012, 48, 1339–1346. [Google Scholar] [CrossRef] [Green Version]
  30. Farza, M.; Busawon, K.; Hammouri, H. Simple nonlinear observers for on-line estimation of kinetic rates in bioreactors. Automatica 1998, 34, 301–318. [Google Scholar] [CrossRef]
  31. Zhang, M.; Li, Z.-T.; Cabassud, M.; Dahhou, B. Unknown input reconstruction: A comparison of system inversion and sliding mode observer based techniques. In Proceedings of the CCC2017 (2017 Chinese Control Conference), Daian, China, 26–27 July 2017; pp. 7172–7177. [Google Scholar]
  32. Alessandri, A.; Bagnerini, P.; Cianci, R. State observation for Lipchitz nonlinear dynamical systems basen on Lyapunov functions and functionals. Mathematics 2020, 8, 1424. [Google Scholar] [CrossRef]
Figure 1. Structure of an interconnected system.
Figure 1. Structure of an interconnected system.
Applsci 10 07966 g001
Figure 2. Structure of the proposed interconnected observer.
Figure 2. Structure of the proposed interconnected observer.
Applsci 10 07966 g002
Figure 3. Output of the process fluid temperature T p . The solid line is the measured value, and the dashed line is estimated by Observer 1. (a) represents Noise-free case while (b) represents Noise-corrupted case.
Figure 3. Output of the process fluid temperature T p . The solid line is the measured value, and the dashed line is estimated by Observer 1. (a) represents Noise-free case while (b) represents Noise-corrupted case.
Applsci 10 07966 g003
Figure 4. Output of the utility fluid temperature T u . The solid line is the measured value, and the dashed line is estimated by Observer 1. (a) represents Noise-free case while (b) represents Noise-corrupted case.
Figure 4. Output of the utility fluid temperature T u . The solid line is the measured value, and the dashed line is estimated by Observer 1. (a) represents Noise-free case while (b) represents Noise-corrupted case.
Applsci 10 07966 g004
Figure 5. The computation and estimation of the process fluid flow rate F p . The solid line is the computed value; the dashed line is estimated by the unknown input observer (UIO); and the dashed, dotted line is estimated by Observer 2. (a) represents Noise-free case while (b) represents Noise-corrupted case.
Figure 5. The computation and estimation of the process fluid flow rate F p . The solid line is the computed value; the dashed line is estimated by the unknown input observer (UIO); and the dashed, dotted line is estimated by Observer 2. (a) represents Noise-free case while (b) represents Noise-corrupted case.
Applsci 10 07966 g005
Figure 6. The computation and estimation of the utility fluid flow rate F u . The solid line is the computed value; the dashed line is estimated by the UIO; and the dashed, dotted line is estimated by Observer 2. (a) represents Noise-free case while (b) represents Noise-corrupted case.
Figure 6. The computation and estimation of the utility fluid flow rate F u . The solid line is the computed value; the dashed line is estimated by the UIO; and the dashed, dotted line is estimated by Observer 2. (a) represents Noise-free case while (b) represents Noise-corrupted case.
Applsci 10 07966 g006
Figure 7. Outlet temperature of the process fluid T p . The solid line is the measured value, and the dashed line is estimated by Observer 1. (a) represents Noise-free case while (b) represents Noise-corrupted case.
Figure 7. Outlet temperature of the process fluid T p . The solid line is the measured value, and the dashed line is estimated by Observer 1. (a) represents Noise-free case while (b) represents Noise-corrupted case.
Applsci 10 07966 g007
Figure 8. Outlet temperature of the utility fluid T u . The solid line is the measured value, and the dashed line is estimated by Observer 1. (a) represents Noise-free case while (b) represents Noise-corrupted case.
Figure 8. Outlet temperature of the utility fluid T u . The solid line is the measured value, and the dashed line is estimated by Observer 1. (a) represents Noise-free case while (b) represents Noise-corrupted case.
Applsci 10 07966 g008
Figure 9. The computation and estimation of the process fluid flow rate F p . The solid line is the computed value; the dashed line is estimated by the UIO; and the dashed, dotted line is estimated by Observer 2. (a) represents Noise-free case while (b) represents Noise-corrupted case.
Figure 9. The computation and estimation of the process fluid flow rate F p . The solid line is the computed value; the dashed line is estimated by the UIO; and the dashed, dotted line is estimated by Observer 2. (a) represents Noise-free case while (b) represents Noise-corrupted case.
Applsci 10 07966 g009
Figure 10. The computation and estimation of the utility fluid flow rate F u . The solid line is the computed value; the dashed line is estimated by the UIO; and the dashed, dotted line is estimated by Observer 2. (a) represents Noise-free case while (b) represents Noise-corrupted case.
Figure 10. The computation and estimation of the utility fluid flow rate F u . The solid line is the computed value; the dashed line is estimated by the UIO; and the dashed, dotted line is estimated by Observer 2. (a) represents Noise-free case while (b) represents Noise-corrupted case.
Applsci 10 07966 g010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, M.; Wu, Q.; Chen, X.; Dahhou, B.; Li, Z. Observer Design for Nonlinear Invertible System from the View of Both Local and Global Levels. Appl. Sci. 2020, 10, 7966. https://doi.org/10.3390/app10227966

AMA Style

Zhang M, Wu Q, Chen X, Dahhou B, Li Z. Observer Design for Nonlinear Invertible System from the View of Both Local and Global Levels. Applied Sciences. 2020; 10(22):7966. https://doi.org/10.3390/app10227966

Chicago/Turabian Style

Zhang, Mei, Qinmu Wu, Xiangping Chen, Boutaïeb Dahhou, and Zetao Li. 2020. "Observer Design for Nonlinear Invertible System from the View of Both Local and Global Levels" Applied Sciences 10, no. 22: 7966. https://doi.org/10.3390/app10227966

APA Style

Zhang, M., Wu, Q., Chen, X., Dahhou, B., & Li, Z. (2020). Observer Design for Nonlinear Invertible System from the View of Both Local and Global Levels. Applied Sciences, 10(22), 7966. https://doi.org/10.3390/app10227966

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