IntegralType EventTriggered Receding Horizon Control of Nonlinear Systems with Additive Disturbance
Abstract
This paper studies the eventtriggered receding horizon control (RHC) of continuoustime nonlinear systems. An integraltype eventtriggered mechanism is proposed to save the communicational resource, and a less conservative robustness constraint is introduced to the RHC scheme for compensating the additive disturbance. Based on these formulations, the designed eventtriggered algorithm is shown to have better performance on avoiding unnecessary communication. Furthermore, the feasibility of the integraltype eventtriggered RHC scheme and the stability of the closedloop system are rigorously investigated. Several sufficient conditions to guarantee these properties are established, which indicate that a tradeoff exists for designing parameters such as the prediction horizon, the disturbance bound, the triggering level, and the contraction rate for the robustness constraint. The effectiveness of the proposed algorithm is illustrated by a numerical example.
I Introduction
Recent research interests and efforts have been directed towards reducing the computation or/and communication load in control systems by using the eventtriggered scheme. This control paradigm employs the socalled eventtriggered mechanism (ETM) to avoid unnecessary communication. Compared with the conventional periodic control, eventtriggered control treats sampling instants as a design parameter while the periodic scheme samples the states at fixed time instants. The eventtriggered paradigm is mainly to design a scheduling mechanism with the capability of determining when the states or sensor outputs should be sampled. The benefit of using this method is that it can provide a smaller communication rate and thus reduce the communication load, especially in cases where the computation and communication resources are scarce. Therefore the eventtriggered control scheme has been widely studied under such resourcelimited scenarios [1, 2, 3, 4, 5, 6, 7].
Several pioneering works have been devoted to build the fundamental blocks of the eventtriggered control [1, 2, 3, 4]. The works [1, 2, 3] used a constant thresholdbased eventtriggered scheme, which enable the next sampling when the norm of state or estimation error exceeds a certain constant bound. Compared with the conventional periodic control, this method is shown to have great advantages in terms of reducing communication rate [1]. In another early work [4], the proposed eventtriggered scheme adopted a socalled relative threshold policy for a class of nonlinear systems, where it allows transmissions when the norm of measurement error exceeds a weighted norm of the state. In addition, a lower bound for interexecution time is proved to exist for avoiding the Zeno behavior. The study of statebased eventtriggered control can be found in [3, 8, 5], and outputbased eventtriggered control has been reported in [9]. In [7], the authors proposed a periodic eventtriggered scheme for reducing the communication rate, where the eventtriggering condition is checked periodically with a fixed time interval. To further reduce the communication rate, the authors in [10] proposed a novel integralbased eventtriggered scheme by incorporating the integral of estimated errors to the eventtriggering condition.
In control research area, Receding Horizon Control (RHC) has been one of the most successful control methodologies. The basic idea of RHC framework is to solve optimization problems online at each sampling instant, and apply the corresponding control action to the plant. In particular, introducing eventtriggered scheme to RHC is of great importance in term of alleviating communication load, thus receives many research studies [11, 12, 13, 14, 15, 16]. By using the eventtriggered scheme, the eventtriggered RHC can avoid unnecessary communication, and thus reduce the frequency of solving heavy computational optimization problems. Therefore, employing eventtriggered scheme on RHC is preferable in term of saving computational cost. Some research has been concentrated on eventtriggered RHC of linear systems [11, 12]. The authors in [12] studied the eventtriggered RHC of linear systems with additive disturbance by using a tubebased approach. There have been also research interests focused on nonlinear systems [13, 14, 15]. In [13], an eventtriggered scheme was proposed for nonlinear systems with additive disturbance by continuously measuring the error between the actual and the predicted trajectories. In order to acquire the benefit of avoiding Zeno behavior, the authors in [14] proposed an eventtriggered mechanism design which can guarantee that the interexecution time is lower bounded. Eventtriggered RHC can also be found in decentralized and networked control systems [17, 15]. In [15], an eventtriggered RHC framework was proposed for stabilizing the distributed nonholonomic systems. However, such RHCbased eventtriggered schemes only use the instantaneous information of the actual and predicted state. The integral of the error between actual and predicted trajectories can be used in the eventtriggered RHC framework for further reducing the communication rate compared with the existing results.
In this paper, we investigate the robust integraltype eventtriggered RHC problem for the continuoustime nonlinear systems with additive disturbance, aiming at alleviating the computational load while ensuring the feasibility of proposed RHC problem and the stability of the closedloop system. The main contributions of this work are threefold:

An integraltype eventtriggered RHC algorithm has been designed for the continuoustime nonlinear system in the presence of additive disturbance. This ETM is proposed by using the integral of the error between the actual and predicted states. By using the ETM, the optimization problem will be only solved when the accumulated error reaches the designed triggering level. The triggering level is also designed to avoid the Zeno behavior;

A novel robustness constraint is proposed to compensate the additive disturbance for the closedloop system. The nominal state in the optimization problem is required to satisfy a timevarying constraint, which is decreasing proportionally to time. Moreover, this constraint will shrink into an ellipsoidal terminal region after a prediction horizon. By using this unique configuration, the optimization problem admits a less conservative initial feasible region.

The feasibility of the integral eventtriggered RHC and the stability of the closedloop systems are thoroughly studied. Sufficient conditions for ensuring the feasibility and stability are provided, respectively. It is also shown that the feasibility and stability conditions are subject to the prediction horizon, the bound of disturbance, the triggering level, and the contraction rate for the robustness constraint. Moreover, there exists a design tradeoff for the parameters when the control performance and the computational load are considered.
The paper is organized as follows: Section II describes the problem formulation of the proposed eventtriggered RHC algorithm. Section III states the feasibility and stability results of the proposed algorithm. Section IV gives a simulation example to verify the effectiveness and computational efficiency of the proposed algorithm.
Notations: The real space is denoted by and the symbol represents the set of all positive integers. For a given matrix , we use and to denote its transpose and inverse. For a symmetric matrix , we write or if is positive definite (PD) or positive semidefinite (PSD). The largest and smallest eigenvalues of are denoted by and , respectively. Given a column vector , represents its Euclidean norm, and is the weighted norm. We also use the notation . Given a continuously differentiable vectorvalued function on , we use to represent its Jacobian matrix.
Ii Problem Formulation
Iia System Dynamics and Optimization Problem
We consider a continuoustime nonlinear system with additive disturbance as follows
(1) 
where is the state variable, is the control input, and is the additive disturbance. The system satisfies and has a Lipschitz constant . The control input , where is a compact set containing the origin. Moreover, the disturbance is also in a compact set containing the origin, which is bounded by . By linearizing the nonlinear system (1) at the equilibrium , we can obtain the linearized statespace model:
(2) 
where and .
In the following, we introduce a conventional assumption for the linearized model (2), which is necessary for analyzing the closedloop performance of the nonlinear system (1).
Assumption 1.
There exists a feedback control law such that the closedloop system matrix is Hurwitz.
In addition, we also make use of a conventional result of the control invariant property of the nonlinear system (1).
Lemma 1.
[18] If is twice continuously differentiable, , is piecewise rightcontinuous and suppose that Assumption 1 holds, then given a stabilizable , and two symmetric positivedefinite matrices and , there exists a constant such that: (1) The Lyapunov equation admits a unique solution , where and satisfies ; (2) and is control invariant by the feedback control law for the nonlinear system (1).
The nominal system of (1) can be defined as
(3) 
where and are the predicted states and control sequence, respectively. Note that the nominal system dynamics (3) will be used for constructing the equality constraint of the following optimization problem. In order to avoid ambiguity, we take explicit notations and for the predicted state and control trajectory at the eventtriggered instant . It should be noticed that the eventtriggered instants are generated by using the integraltype ETM, which will be elaborated after introducing the optimization control problem.
Then the nonlinear optimization problem can be formulated as
(4)  
s.t.  (5)  
(6)  
(7) 
The cost function is defined as follows
(8) 
where , , is the identity matrix, is defined by using the method from Lemma 1, is the prediction horizon, is the designed parameter for defining the terminal set, is the scaling ratio, and is the contraction rate for the robustness constraint (7). It should be noted that the terminal constraint is .
By using the aforementioned system dynamics and optimization problem, we construct the following closedloop system as
(9) 
where is the eventtriggered instant. In this control framework, we use the optimal control sequence as the control input generated by solving the nonlinear optimization problem (4). In order to reduce the communication between the controller and the actuator, we take an integraltype ETM for determining when the next optimization should be conducted, i.e. solving the optimization problem and transmitting the optimal control sequence to the actuator.
Remark 1.
In the optimization problem (4), the robustness constraint (7) is used for compensating the additive disturbance when applying the optimal control sequence to the closedloop system (9). The shrinking constraintbased method for robust RHC has been firstly used in [19]. In their configuration, the robustness constraint shrinks very fast to the terminal constraint as time evolves. In our approach, a less conservative constraint is proposed, where the shrinking rate is a constant. Intuitively, our proposed robustness constraint will provide a larger feasible set for solving the optimization problem.
IiB Integraltype Eventtriggered Mechanism
For the nonlinear system (1), we assume that the eventtriggered instants set is with , where is determined by the proposed eventtriggered scheme (10). We also make the assumption that there is no time delay and inaccuracy when transmitting the sensor measurements to the digital controller. By using this configuration, an integraltype ETM is introduced for scheduling and implementing the sampling tasks, i.e. determining the eventtriggered instants . As shown in Fig. 1, the ETM produces an ON/OFF signal to the sensor by measuring the error between the actual state and the predicted optimal state. Based on this configuration, we can acquire the benefit of reducing both the computational and communicational load. The integraltype eventtriggering condition is designed as
(10) 
where is a minimum time instant satisfying the triggering condition (10). This ETM can take account of accumulated error between the measured state and optimal predicted state generated by the RHC algorithm (4) over the current period , which is different from the eventtriggered setting in [14]. The following theorem shows some important properties of the proposed integraltype ETM.
Theorem 1.
Proof.
This proof can be done by two steps.
Step 1: The upper bound of interexecution intervals is . From the design of the integraltype ETM, it can be directly deduced that all the intervals is less than or equal to the prediction horizon .
Step 2: The lower bound for interexecution intervals can be designed as . To prove this result, we firstly consider the upper bound for at . We assume here that the sensor measurements of the states are accurate, thus it follows that . By using the triangle inequality, we have Then by applying the integral form of GronwallBellman inequality, it can be obtained that Substituting the previous inequality to (10), we can deduce that Since is strictly larger than zero for , we can choose as (11) by setting . Therefore, it can be guaranteed that the lower bound of the triggered time interval is . The proof is completed. ∎
Remark 2.
Note that the designed ETM is based on the integral of the error between actual states and predicted states. The main difference of the integraltype ETM is that the accumulated error between two consecutive eventtriggered instants is taken into account. It is worthwhile to point out that this eventtriggered scheme is only valid for the system in the presence of disturbance.
Iii Main Results
Iiia Feasibility Analysis
Following a conventional setup for RHC framework, we construct a classical feasible control sequence for the optimization problem (4). This same control policy has been widely exploited by [20, 21], and it can be given as
(12) 
Then the suboptimal feasible control and state trajectory evolves as
(13) 
Before presenting the result of this section, we would like to introduce a lemma which will be used in the following analysis.
Lemma 2.
Given a continuously differentiable vectorvalued function on , then the following inequality holds
(14) 
Proof.
For every , we have two results: Subtracting the aforementioned two equations, it yields By employing the triangle inequality, we can deduce from the above equality that Since every defined on the closed interval is equal or less than the right side of the above inequality, thus the result (14) holds. ∎
To simplify the analysis of integraltype eventtriggered configuration, we make use of a term for , which can be found continuously differentiable. Note that is defined as the candidate state trajectory generated by the nominal system (13), and is the solution of the optimization problem (4). In addition, we let the suboptimal feasible state by sampling the state at the triggered time instant . For convenience, we denote , where is the square root of the . Then we can propose the following result.
Corollary 1.
Given and defined on , the following inequality holds
(15) 
Proof.
By Lemma 2, we can obtain that
(16) 
Note that , then it can be deduced from the above inequality that
(17) 
Next we show that the upper bound of is related to the triggering level . By using the result from Theorem 1, the upper bound of is for . Thus, it can be obtained that for . Note that we also have for , where the triggering level is designed as the integral of from to . By following simple calculation, we can obtain that
(18) 
and consequently it follows that
(19) 
where . Since the function gets its maximum at , the above inequality becomes
(20) 
According to GronwallBellman inequality, one can obtain (15) by substituting (20) to (17). The proof is thus completed. ∎
Now we can analyze the iterative feasibility of the RHC problem (4), implying that if the RHC problem admits a solution at current time instant then a feasible solution exists for the next time instant. To prove this result, we use a conventional feasible control sequence candidate at time instant defined in (12), where for and for . In the following theorem, we will show that the designed control sequence candidate can steer the feasible state into if some conditions can be satisfied. In addition, it is also necessary to show that the candidate state will remain in the designed state constraint (7).
Assumption 2.
The optimization problem (4) admits a feasible solution for the initial time .
Theorem 2.
Proof.
First, we show that the designed control sequence for drives into , i.e. . Let us construct an error norm for . By using Corollary 1, we can obtain that
(25) 
Then it follows that . By using the Triangle inequality, we have
(26) 
which implies that .
Note from Theorem 1 that the designed triggering level . In order to steer the candidate state trajectory into , one can simply deduce that the following inequality must holds
(27) 
From (27), it can be also obtained that the maximum bound for disturbance satisfies .
Second, we consider the candidate trajectory for . By using Lemma 1, we can verify that is an invariant set for the closedloop system . Consequently, we can deduce that . By the virtue of comparison principle for , it follows that
(28) 
which indicates that . By using Theorem 1, we can have . To obtain , it is equivalent to show that . With some calculation, one can obtain to guarantee the previous inequality holds. Similar argument can be found in [14].
Third, we show that will satisfy the state constraint (7). For , one can get
(29) 
which can be easily derived from (15). Then we need to prove
(30) 
By some calculation, it can be obtain that . For , it can be derived from (28) that
(31) 
In order to prove , it is equivalent to show
(32) 
For brevity, we denote , and it turns out that . By evaluating the derivative of , it can be verified that is nonpositive for , which indicates . Finally, the designing parameter should be set as for guaranteeing the satisfaction of the proposed robustness constraint. The proof is completed. ∎
Remark 3.
Note from Theorem 2 that the feasibility can be affected by the prediction horizon , the Lipschitz constant , and the disturbance bound . In order to achieve the recursive feasibility, the prediction horizon should be lower bounded, and the design parameter in (7) should be lower bounded as well. Specifically, a lager leads to a larger initial feasible region. It should be also noted that the maximum allowable disturbance can be decided as (27), which shows that the allowable disturbance and the prediction horizon are correlated to each other. Thus the tradeoff between these two parameters must be taken into account when designing the algorithm.
IiiB Stability Analysis
In this part, we investigate the closedloop stability of the proposed integraltype eventtriggered RHC. In the following theorem, we mainly analyze the nonincreasing properties of the cost function in (8). Due to the existence of the disturbance, it is worthwhile to point out that the closedloop stability can be achieved to converge to an invariant set. Since we use the RHC configuration, the analysis for stability can be divided into two steps: One is to ensure that the optimal trajectory will enter the terminal set in finite time; the other is to prove that the closedloop system is stable after the state enters the terminal set .
Theorem 3.
Proof.
This theorem will be proved by two steps.
Step 1: For all initial state , we aim to show the state trajectory enters in finite time. In this situation, we construct an error term of two Lyapunov functions as follows . Expanding this term yields
(34) 
Substituting (12) to the above equation, we can obtain that
(35) 
Note from Lemma 1 that . Taking integral from to of the above inequality yields
(36) 
Applying this fact to , it can be shown that
(37) 
To analyze the above inequality, we firstly consider the term