Next Article in Journal
Progressive Development and Challenges Faced by Solar Rotary Desiccant-Based Air-Conditioning Systems: A Review
Previous Article in Journal
Evaluation of the Thermal and Morphological Properties of γ-Irradiated Chitosan-Glycerol-Based Polymeric Films
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonlinear Offset-Free Model Predictive Control based on Dynamic PLS Framework

1
School of Information and Control Engineering, Liaoning Petrochemical University, Fushun 113001, China
2
Institute of Intelligence Science and Engineering, Shenzhen Polytechnic, Shenzhen 518055, China
*
Author to whom correspondence should be addressed.
Processes 2021, 9(10), 1784; https://doi.org/10.3390/pr9101784
Submission received: 19 August 2021 / Revised: 28 September 2021 / Accepted: 5 October 2021 / Published: 7 October 2021
(This article belongs to the Section Automation Control Systems)

Abstract

:
A nonlinear offset-free model predictive control based on a dynamic partial least square (PLS) framework is proposed in this paper. A multi-output multi-input system is projected into latent variable space by a PLS outer model. For each latent variable model, the T–S fuzzy model is used to describe the nonlinear characteristics of the system; while the state-space model is used in T–S fuzzy model consequent parameters to describe the dynamic characteristics. A disturbance model is introduced in the state-space model. For model state variables, a state observer is used to compensate for the mismatch of the model. The case study results for the pH neutralization process show that the MPC controller based on this method can guarantee the tracking performance of the nonlinear system without static error.

1. Introduction

Nonlinearity is a common characteristic in industrial processes. Many control algorithms are proposed for nonlinear systems [1]. A nonlinear model predictive control (MPC) is one of its significant branches [2]. A nonlinear MPC can be used to solve the control problem of nonlinear process. However, a mismatch between model and actual processes exists [3,4], which is intensified by unmeasurable disturbance. This leads to a decline in MPC control performance. Conventional nonlinear MPC cannot guarantee a zero-static error performance [5].
Many researchers have proposed a variety of solutions to the offset-free control problem for nonlinear systems. Jap et al. [6] presented a stochastic-tube MPC strategy that systematically handles plant–model mismatch, guarantees stability in the presence of changing operating conditions, and ensures offset-free tracking for all reachable operating conditions. Wang et al. [7] proposed an offset-free robust model predictive control for the current control of a three-phase induction motor where an incremental model is adopted to improve the parameter robustness and to realize offset-free prediction, and an improved current observer is designed to realize time-delay compensation in control variable optimization. Awha et al. [8] proposed an offset-free control scheme based on a multiple model predictive controller (MMPC) and an adaptive integral action controller for nonlinear processes; the multiple model is used to capture the essence of the nonlinear process, a min–max approach is utilized to counter the effect of parametric uncertainty between the linear models and the nonlinear process, and an adaptive integral action controller is run in parallel to the proposed controller.
The neural network and fuzzy models are two nonlinear models that are widely used in nonlinear modeling and control [9]. Benrabah et al. [10] proposed a feed-forward multilayer neural network to approximate nonlinear mapping, and a teaching–learning-based optimization strategy is used to resolve the optimization problem in MPC. Wu et al. [11] introduced a hybrid modeling method to integrate the first-principles model and the recurrent neural network. With the obtained model, a nonlinear MPC is designed and applied to a chemical process. L. Kong and J. Yuan [12] proposed a fuzzy model predictive controller for a boiler-turbine system, which guarantees asymptotic stability and satisfies the input constraints. Escano et al. [13] proposed a complexity reduction Takagi–Sugeno fuzzy model that is obtained by finding a function basis via a functional principal component analysis and used for MPC. It is applied in a low-cost microcontroller system. Han et al. [14] used an adaptive fuzzy neural network identifier to match the nonlinear behaviors of wastewater treatment processes with process data and proposed a multi-objective control strategy that had low computational cost.
Partial least square (PLS) proved to be rewarding in the data-driven modeling area and has been applied to many areas such as quality prediction, process monitoring, and chemometrics [15]. The principle of PLS is to maximize the covariance between the input and output data sets in a reduced space. Its main advantages include: (1) original data dimensionality reduction; (2) multivariate auto-pairing; (3) modeling of non-square system. In order to handle the dynamic modeling problem, a variety of methods combing the dynamic model with PLS have been proposed in the past decades. Yining Dong and S. Joe Qin [16] proposed a dynamic inner PLS model where the explicit dynamic inner model was given, and the inner and outer models were made consistent at the same time. A method that combined the autoregressive exogenous structure and PLS was proposed by Kaspar and Ray [17]. Qinghua Chi et al. [18] extended this method and discussed the relevant identification method in the inner PLS. Junghui Chen et al. [19] proposed another dynamic PLS framework with the ARX model. Besides these modeling methods, many attempts have been made to put forward new control strategies that compromise the merits of PLS. Junghui Chen et al. [19] proposed a novel decoupling PID strategy with PLS. Lü and Liang [20] proposed a multi-loop constrained MPC scheme. Jianhua Zhang [21] proposed a robust control based on PLS. The studies on nonlinear MPC are relatively fewer. There are two main difficulties. One is that PLS is a linear regression method. It is unable to deal with nonlinear problems. The other is that conventional PLS could not describe the dynamic characteristics of the system. In order to address the first shortcoming, some researchers proposed some nonlinear PLS methods by introducing a nonlinear modeling method [22,23,24,25]. Some researchers extended these methods to nonlinear control problems and proposed nonlinear PLS control methods. Zhao et al. [26] used neural networks with an ARX structure as internal models to describe the nonlinear characteristics of the system and designed an internal model controller. Chi et al. [27] proposed a multi-model predictive control strategy, which divided the process into several blocks according to the working area and used the area models to describe the system and design the controller.
This paper attempts to solve the nonlinear dynamic PLS modeling problem. The main work of this paper is as follows: a brief review of the T–S fuzzy model and the fuzzy model used in this paper is presented in Section 2. In Section 3, a T–S fuzzy model is used in the inner model of a dynamic PLS framework to describe the nonlinear character of the system. Based on the model in Section 3, a nonlinear MPC in a dynamic PLS framework is given in Section 4, and an offset-free control method based on it is discussed. In Section 5, a case study based on the PH process is described to demonstrate the merit of the proposed method. Finally, conclusions are drawn in Section 6.

2. T–S-Fuzzy Model-Overview

The fuzzy model is a modeling method based on the fuzzy combination rule. The fuzzy rule is expressed in the form of “IF...THEN...”, where “IF...” is called the fuzzy rule premise and “THEN...” is called the fuzzy rule consequent. Fuzzy models can be divided into three categories according to the different consequents: (1) Mamdani type (fuzzy language model); (2) Fuzzy relational type; (3) Takagi–Sugeno (T–S) type.
The T–S model expresses the global nonlinear system with a local linearized model under each fuzzy rule. Such a model form is advantageous for the design of the corresponding controller of each fuzzy rule and for the use of the principle of parallel distributed compensation (PDC) to obtain the overall controller output. In this section, the T–S model is used as the inner model of PLS to describe the dynamic nonlinear characteristic of the system and to design the corresponding controller. The inner model of PLS is a series of SISO sub-models; however, only the SISO T–S model is introduced here.
In the T–S model [28], the form of the fuzzy rule describing the system dynamics is shown in Equation (1):
Li: IF z(k) is Ωi, THEN yi(k) = Ai y(k-1) + Bi x(k), i = 1,…,K.
where Li is the i-th fuzzy rule; z(k) is the premise input vector composed of input and output data. Ωi is the fuzzy set. Here, the ARX model is used as the consequent parameter of the fuzzy rule. Ai and Bi are polynomials, and K is the number of fuzzy rule. The fuzzy set Ωi is defined as:
μ i ( z ( k ) ) : R [ 0 , 1 ] .
where μ i ( z ( k ) ) is the membership function of the fuzzy set Ωi, which characterizes the degree that z ( k ) belongs to the fuzzy rule Li.
With Equations (1) and (2), the total output of the T–S model is obtained as follows:
y ( k ) = i = 1 K μ i ( z ( k ) ) y i ( k ) i = 1 K μ i ( z ( k ) ) .
The fuzzy model can be derived by using the first principle model of the nonlinear system. But this method is not universal. A more common method of establishing a fuzzy model is by identifying the input and output data. The following is an introduction to the identification of the T–S model premise (the membership function) and the T–S model consequent (the linear equation parameters).
The membership function is obtained by clustering the algorithm. Commonly used clustering algorithms include K-means clustering, fuzzy, C-means (FCM) clustering, and Gustafson–Kessel (GK) clustering.
The K-means theory was first seen in the E. W. Forgy article [29], which is the most widely used clustering algorithm because of its simplicity and efficiency. In the K-means algorithm, the distance between the sample and the cluster center is determined by Euclidean distance. The goal is to minimize the distance in the same cluster and to maximize the distance between different clusters by minimizing the objective function, as shown in Equation (4):
J = i = 1 K k = 1 l Z i k z ( k ) v i 2 .
where l denotes the number of samples, v i represents the cluster center, Z i k represents the membership function and satisfies Z i k = 1 , z ( k ) cluster i 0 ,   else .
The disadvantage of the K-means clustering is that it is not global optimum but local optimum, and is only suitable for clustering numerical data. A direct fuzzy K-means clustering method of classifying the data is by a flexible partitioning method where the membership function is not in the form of 0 or 1, but a value between 0 and 1. Bezdek [30] introduced the fuzzy index m into the objective function Equation (4) and rewrote it as follows:
J ( U , V ) = i = 1 K k = 1 l ( μ i ( k ) ) m z ( k ) v i 2 .
where U = [ μ 1 ( k ) , , μ K ( k ) ] is the fuzzy partition of z, μ i ( k ) is the membership degree of z ( k ) corresponding to cluster center vi, which satisfies μ i ( k ) [ 0 , 1 ] and i = 1 K μ i ( k ) = 1 , V = [ v 1 , , v K ] is the cluster center vector, and m is the fuzzy index.
A fuzzy K-means clustering method based on the objective function Equation (5) is called fuzzy C-means clustering. The solution of Equation (5) based on the Lagrangian multiplier method is:
μ i ( k ) = 1 j = 1 K d i k d j k 2 m 1 , where   d i k 2 = z ( k ) v i 2 .
The equation of the cluster center is
v i = k = 1 l μ i m ( k ) z ( k ) k = 1 l μ i m ( k )
The GK clustering method [31] introduces an additional symmetric positive definite parameter matrix in each cluster, which is similar to the multivariate normal distribution covariance matrix, so that it can be used to deal with elliptic clustering problems. This means that the solution to this positive parameter matrix must be included in the optimization step of Equation (5) to solve the cluster center. The objective function of GK clustering is defined as follows:
J ( U , V , A ) = i = 1 K k = 1 l ( μ i ( k ) ) m d i k , A i 2
The Euclidean distance in K-means clustering and fuzzy C-means clustering is replaced by the Mahalanobis distance in GK clustering, which is defined as:
d i k , A i 2 = ( z ( k ) v i ) T A i ( z ( k ) v i )
where A i = ρ i det ( F i ) 1 / p F i 1 , ρ i > 0 i ,
F i = k = 1 l μ i m ( k ) ( z ( k ) v i ) T ( z ( k ) v i ) k = 1 l μ i m ( k )
F = [ F 1 , , F K ] is the fuzzy covariance matrix of GK clustering.
In this paper, GK clustering is used to obtain the membership function. The steps of the solution are as follows:
Step 1: Parameter initialization. Select the number of clusters K, the fuzzy index m, and the convergence error ε > 0 . The fuzzy partition matrix U is then randomly initialized;
Step 2: Calculate the cluster center V according to Equation (7);
Step 3: Calculate the covariance matrix F according to Equation (10);
Step 4: Calculate the Mahalanobis distance d i k , A i 2 according to Equation (9);
Step 5: Update the fuzzy partition according to Equation (11);
μ i ( k ) = 1 j = 1 K d i k d j k 2 m 1 , d i k > 0 μ i ( k ) = 1 , d i k = 0
Step 6: Repeat Step 2 up to Step 5 until the sum of the squared error of two iterations of the fuzzy partition matrices is less than the given convergence error.
The identification of the consequent parameters can be performed after U, V, and F are obtained. The ARX model shown in Equation (1) is the consequent, and the least squares method can be used for identification.
For the subsequent controller design, the ARX structure consequent is rewritten as a state space based form. Rewrite Equation (1) as:
L i :   IF   z ( k )   is   Ω i ,   THEN x i ( k + 1 ) = A i x i ( k ) + B i u ( k ) y i ( k ) = C i x i ( k ) ,   i = 1 , , K
where Ai, Bi, and Ci are state equation parameters and are satisfied A i = a i , 1 1 0 0 a i , 2 0 1 0 a i , n a 1 0 0 1 a i , n a 0 0 0 , B i = b i , 1 b i , n b , C i = 1 0 0 T . Ai, Bi, and Ci can be converted from the ARX model of Equation (1).
According to the fuzzy inference method [32], the fuzzy inference model and the single-point fuzzification theory generated by the central average defuzzification can be expressed as the global fuzzy model (shown in Equation (13)).
x ( k + 1 ) = A x ( k ) + B u ( k ) y ( k ) = C x ( k )
where A = i = 1 K μ i ( k ) A i , B = i = 1 K μ i ( k ) B i , C = i = 1 K μ i ( k ) C i .
Equation (13) replacing Equation (12) is also known as the dynamic fuzzy global state space model.

3. Nonlinear Dynamic PLS Model

3.1. Dynamic PLS Model

PLS was first proposed by Herman Wold’s original nonlinear iterative partial least square (NIPALS) algorithm [33]. It requires two data matrices, input data matrix X = (xij)l×m and output data matrix Y = (yij) l×n, where l denotes the number of sample, while m and n denote the number of input variables and output variables, respectively. The raw data X and Y need to be scaled. In this paper, we suppose that X and Y are scaled by matrix W x and W y , respectively. In the NIPALS method, X and Y are projected into latent variable space as a sum of a series of vectors, as shown in Equation (14). This is called the outer model:
X = r = 1 R t r p r T + E * = T P T + E * Y = r = 1 R u r q r T + F * = U Q T + F *
where T = [ t 1 , t 2 , , t R ] and U = [ u 1 , u 2 , , u R ] are the score matrices of X and Y, respectively. P = [ p 1 , p 2 , , p R ] T and Q = [ q 1 , q 2 , , q R ] T are the loading matrices of X and Y, respectively. E* and F* are residual matrices of X and Y, respectively. R is the number of latent variables.
The relationship between X and Y is constructed as the relationship between tr and ur in latent variable space, which is called the inner model:
u r = b r t r
where b r = u r T t r t r T t r ; B = ( T T T ) 1 T T U ; U = T B .
Combing Equation (14) with Equation (15), the total regression model explained in Y can be written as follows:
Y = T B Q T + F * = r = 1 R b r t r q r T + F *

3.2. Nonlinear Dynamic Fuzzy PLS Model

The commonly used nonlinear PLS modeling methods usually have two forms. One is to modify the outer and inner models of PLS at the same time in order to fit the nonlinear characteristics. Such a method is adopted by Malthouse et al. [34], projection pursuit regression is used in the outer model to extract the principal component of input and output data, and the feedforward neural network is used in the inner model to regress the score matrix T and U. The other is to add nonlinearity in the inner model while keeping the outer model unchanged. There has been a lot of research on this method. Qin and McAvoy [35] used the neural network in the inner model. Bang et al. [22] proposed the fuzzy PLS method, which uses the fuzzy model as the inner model of PLS. Abdel-Rahman and Lim improved the fuzzy PLS model and proposed to use the type-two fuzzy model as the inner model of PLS, which is called QFPLS. These methods, no matter what form is used to improve the original PLS method, are all for static data modeling, just as the linear PLS method. In this paper, the T–S model introduced in Section 2 is used as the inner model of PLS to build a nonlinear dynamic PLS model. In the latent variable space, each sub-model is a SISO model. Therefore, the fuzzy model Equation (13) is used as the sub-model. The input and output of the model are the corresponding tr and ur scores.
The r-th inner fuzzy model in latent variable space is as follows:
L i :   IF   z ( k )   is   Ω i ,   THEN   x r ( k + 1 ) = A i x r ( k ) + B i t r ( k ) u r ( k ) = C i x r ( k ) ,   i = 1 , , K
where z ( k ) = [ u r ( k 1 ) , u r ( k n a ) , t r ( k 1 ) , t r ( k n b ) ] .
The global form of Equation (17) is the following:
x r ( k + 1 ) = A r x r ( k ) + B r t r ( k ) u r ( k ) = C r x r ( k ) or   x r ( k + 1 ) = i = 1 K μ i ( k ) A r , i x r ( k ) + B r , i t r ( k ) u r ( k ) = i = 1 K μ i ( k ) C r , i x r ( k )
where A r = i = 1 K μ i ( k ) A r , i , B r = i = 1 K μ i ( k ) B r , i , C r = i = 1 K μ i ( k ) C r , i .
The modeling process of dynamic fuzzy PLS can be summarized in the following six steps:
Step 1: The data matrix X and Y are normalized to X0 and Y0, and the normalized coefficient matrix Wx and Wy are obtained; r = 1;
Step 2: Let the first row of Y0 be ur; obtaining score vectors tr and ur by iterative solution using Equation (19) up to Equation (22) with the NIPALS method;
w r T = u r T X r 1 ( u r T u r )
t r = X r 1 w r X r 1 w r
q r = Y r 1 T t r Y r 1 T t r
u r = Y r 1 q r
Step 3: tr and ur are used as input and output data of the T–S fuzzy model. Identify the fuzzy model premise by using Equation (8) up to Equation (10) and the GK clustering steps, and obtain the fuzzy partition U, cluster center V, and covariance matrix F;
Step 4: Obtain the fuzzy model consequent of Equation (12) or Equation (13) by identifying consequent parameters with the least square. Then, obtain the latent variable output u r * with it;
Step 5: Calculate the loading vectors of X and Y with Equations (23) and (24);
p = X r 1 T t r ( t r T t r )
q r * = Y r 1 T u r * Y r 1 T u r *
Step 6: Let X r = X r 1 t r p r T and Y r = Y r 1 u r * q r * T . Xr and Yr are used as new data matrices X0 and Y0. Return to Step1 until the PLS model meets the accuracy requirements.

4. Offset-Free Model Predictive Control Based on a Nonlinear Dynamic PLS Framework

4.1. Dynamic PLS Control Framework

After the dynamic PLS is performed, the MIMO system is described as a series of SISO sub-models in latent variable space. Due to the orthogonality of the score vector pairs tr and vr, these sub-models are completely decoupled. Hence, many conventional SISO control strategies can be applied to it. This method was proposed by Kaspar and Ray [17], which is called the dynamic PLS control framework. This framework is shown in Figure 1, where Q+ is the appropriate inverse of Q. Setpoint Yset and system output Y are scaled by scaling matrix W y 1 , and projected into latent variable space by inverse loading matrix Q +. The resulting errors are used as the inputs of the controllers. The control law T is back projected into the original space by loading matrix P and anti-scaled by matrix W x , then they are used as the input of actual system.

4.2. Offset-Free Fuzzy Model Predictive Control Based on Dynamic PLS Framework

Based on the dynamic inner model Equation (18), the model predictive controller can be designed. At time k, the fuzzy model premise vector z(k) is known. The fuzzy partition μ i ( k ) at k can then be calculated. Following Equation (18), the r-th sub-model output u r ( k ) at time k can be obtained. It should be noted that μ i ( k ) is a time-varying variable. It needs to be calculated at every new moment. Hence, although Ai, Bi, and Ci in Equation (18) are time invariant models, Ar, Br, and Cr need to be calculated at every new moment. When new models are available, the control law could be solved.
Similar to linear systems, the feedback of control is the model state xr. When the state space model in the dynamic PLS framework matches the plant model accurately, xr can explain the system output Y well. Due to the decoupling scheme of PLS, the mismatch error F * ( k ) is unavoidably created. This leads to poor control performance or a steady-state error. In order to solve this problem, the output feedback is introduced into the controller through the observer. Define the error between the actual system output u Y , r ( k ) in latent variable space and the inner model output u r ( k k 1 ) with the following formula:
d r ( k ) = u Y , r ( k ) u r ( k )
Put the error into Equation (18) as:
x r ( k + 1 ) = A r x r ( k ) + B r t r ( k ) + E r d r ( k ) u r ( k ) = C r x r ( k ) + d r ( k )
where E r = i = 1 K μ i ( k ) E r , i , E r , i = a i , 1 a i , 2 a i , n a 1 a i , n a T .
Establish the fuzzy observer shown in Equation (27).
Li: IF z(k) is Ωi, THEN
x ^ r ( k + 1 ) = A i x ^ r ( k ) + B i t r ( k ) + E r , i d ^ r ( k ) + F r , i ( u r ( k ) u ^ r ( k ) ) u ^ r ( k ) = C i x ^ r ( k ) + d ^ r ( k ) ,   i = 1 , , K
where F r , i is the observer gain matrix, x ^ r k and y ^ r ( k ) are the estimation of state and output, respectively. Equation (27) can be augmented and rewritten in a global fuzzy model form as follows:
x ^ r ( k + 1 ) d ^ r ( k + 1 ) = A r E r 0 1 x ^ r ( k ) d ^ r ( k ) + B r 0 t r ( k ) + F r 0 ( u r ( k ) u ^ r ( k ) ) u ^ r ( k ) = C r x ^ r ( k ) + d ^ r ( k )
Bring the output equation in Equations (26) and (28) to the sate equation in Equation (28) as:
x ^ r ( k + 1 ) d ^ r ( k + 1 ) = i = 1 K j = 1 K μ i ( k ) μ j ( k ) A r , i E r , i 0 1 x ^ r ( k ) d ^ r ( k ) + B r , i 0 t r ( k ) + F r , i C r , j 0 ( x r ( k ) x ^ r ( k ) + d r ( k ) d ^ r ( k ) )
Equation (26) is subtracted from Equation (29). The error state of observer is as follows:
x ^ ˙ r ( k + 1 ) d ^ ˙ r ( k + 1 ) = i = 1 K j = 1 K μ i ( k ) μ j ( k ) A r , i E r , i 0 1 F r , i 0 [ C r , j 1 ] x ^ ˙ r ( k ) d ^ ˙ r ( k )
where x ^ ˙ r ( k ) = x r ( k ) x ^ r ( k ) , d ^ ˙ r ( k ) = d r ( k ) d ^ r ( k ) .
Thus, in the PLS framework, r-th MPC is based on the minimization of a quadratic cost function as:
J r = j = 1 N p , r δ r ( j ) u set , r ( k + j ) u ^ r ( k + j k ) 2 2 + j = 1 N u , r λ r ( j ) Δ t ^ r ( k + j 1 k ) 2 2
where Np,r and Nc,r denote prediction horizon and control horizon, respectively. u set , r is the setpoint in the latent space, which is transformed from the setpoint in original space.

4.3. Getting a Proper Observer Gain F r , i

Wang et al. [36] pointed out that the sub-model (Equation (27)) stability of each fuzzy rule in the fuzzy model stable cannot guarantee a global model (Equation (28)) stability. Therefore, the stability of the observer should be discussed before solving the observer gain F r , i . The stability of the T–S fuzzy model is mainly analyzed by the Lyapunov stability theory [37]. The global Lyapunov function is used for stability analysis in this paper.
When t r ( k ) = 0 , one can determine the Lyapunov function for the discrete fuzzy system of Equation (26) discussed in this paper as follows:
V ( x ) = x T P 1 x
where P is the positive definite symmetric matrix, x is the state variable. Similar conclusions to Theorem 1 can be obtained from Ref. [38].
Theorem 1.
If a symmetric positive definite matrix P exists, it satisfies the linear matrix inequality (LMI) Equation (33). The discrete fuzzy system described in Equation (26) is globally exponentially stable.
A i T P 1 A i P 1 < 0 ,   i = 1 , , K
Using the Schur complements [39], Equation (33) is equivalent to the following formula:
P A i P P T A i T P < 0 ,   i = 1 , , K  
For the fuzzy model of Equation (30), one can obtain the conclusion of Theorem 2.
Theorem 2.
If there are positive definite symmetric matrices P and a set of matrices Qi (i = 1,…,K) satisfying LMI such as Equations (35) and (36), then the observer error state equation of Equation (30) is globally exponentially stable. Moreover, observer gain is F ^ r , i = P 1 Q i , i = 1,…,K.
P P A ^ r , i Q i C ^ r , i A ^ r , i T P C ^ r , i T Q i T P < 0 ,   i = 1 , , K ;
P P A ^ r , i Q i C ^ r , j A ^ r , i T P C ^ r , j T Q i T P + P P A ^ r , j Q j C ^ r , i A ^ r , j T P C ^ r , i T Q j T P < 0
where  i < j K
Proof of Theorem 2.
Equation (30) can be rewritten as follows:
x ^ ˙ r ( k + 1 ) = i = 1 K μ i ( k ) μ i ( k ) A ^ r , i F ^ r , i C ^ r , i G i , j x ^ ˙ r ( k ) + 2 i = 1 K j > i μ i ( k ) μ j ( k ) A ^ r , i F ^ r , i C ^ r , j + A ^ r , j F ^ r , j C ^ r , i 2 x ^ ˙ r ( k )
where A ^ r , i = A r , i E r , i 0 1 , F ^ r , i = F r , i 0 , C ^ r , i = [ C r , i 1 ] , A ^ r , j = A r , j E r , j 0 1 , F ^ r , j = F r , j 0 , C ^ r , j = [ C r , j 1 ] .
According to Theorem 1, one can get the conclusion that the discrete fuzzy system described in Equation (37) is exponentially stable if there exists a symmetric positive definite matrix P satisfying the LMI Equations (38) and (39).
A ^ r , i F ^ r , i C ^ r , i T P 1 A ^ r , i F ^ r , i C ^ r , i P 1 < 0
A ^ r , i F ^ r , i C ^ r , j + A ^ r , j F ^ r , j C ^ r , i 2 T P 1 A ^ r , i F ^ r , i C ^ r , j + A ^ r , j F ^ r , j C ^ r , i 2 P 1 < 0
where i < j K
Equation (38) can be rewritten as follows:
A ^ r , i F ^ r , i C ^ r , i T P T P 3 P A ^ r , i F ^ r , i C ^ r , i P 1 < 0
P A ^ r , i P F ^ r , i C ^ r , i T P 3 P A ^ r , i P F ^ r , i C ^ r , i P 1 < 0
Use the Schur complements [39] to satisfy the following LMI establishment as:
P P A ^ r , i P F ^ r , i C ^ r , i A ^ r , i T P C ^ r , i T P F ^ r , i T P < 0
Let Q i = F ^ r , i C ^ r , i , then Equation (42) is rewritten as:
P P A ^ r , i Q i C ^ r , i A ^ r , i T P C ^ r , i T Q i T P < 0 ,   i = 1 , , K .
Empathy Equation (39) is rewritten as follows:
P P A ^ r , i Q i C ^ r , j A ^ r , i T P C ^ r , j T Q i T P + P P A ^ r , j Q j C ^ r , i A ^ r , j T P C ^ r , i T Q j T P < 0 , i < j K
By solving Equations (43) and (44), one can get the feasible P and Qi. Additionally, observer gain is F ^ r , i = P 1 Q i , i = 1,…,K. □
Based on Theorem 2, Equation (30) of global stability is obtained; it can then further ensure the stability of the proposed method.

5. Case Study

A typical nonlinear process, the PH neutralization titration process, is used as a case study to verify the effectiveness of the proposed method. The schematic illustration is shown in Figure 2. The process has three inputs: q1 (HNO3), q2 (NaHCO3), and q3 (NaOH). The three enter the reactor and are fully mixed and neutralized after mixing. In this case study, q1 and q2 are selected as the manipulation variables, and q3 is a constant value. The outputs are the PH value and the level of reactor h (cm). The model of the process consists of three nonlinear ordinary differential equations and a nonlinear equation, which is shown in Equation (45) [40]. The parameters of the model are shown in Table 1.
{ d h d t = 1 A ( q 1 + q 2 + q 3 C v h 0.5 ) d W a 4 d t = 1 A h [ ( W a 1 W a 4 ) q 1 + ( W a 2 W a 4 ) q 2 + ( W a 3 W a 4 ) q 3 ] d W b 4 d t = 1 A h [ ( W b 1 W b 4 ) q 1 + ( W b 2 W b 4 ) q 2 + ( W b 3 W b 4 ) q 3 ] W a 4 + 10 14 p H + W b 4 1 + 2 × 10 p H p K 2 1 + 10 p K 1 p H + 10 p H p K 2 10 p H = 0 W a 4 + 10 14 p H + W b 4 1 + 2 × 10 p H p K 2 1 + 10 p K 1 p H + 10 p H p K 2 10 p H = 0
In order to build a dynamic PLS model, a set of pseudo random signals is generated to excite the system and to generate identification data. This set of signals should include the range of manipulating variables as far as possible. The data set (1500 samples) for identification is shown in Figure 3. The sampling-time interval for this case study is 15 s.
In this case study, the number of latent variables in the proposed model R is 2. The fuzzy rule of each sub-latent variable model is 2. The model order of consequent is na = 1 and nb = 1, and the fuzzy exponent is 2. The other data set (1500 samples) generated is used for cross validation; the validation results are shown in Figure 4. Here, a new norm shown in Equation (46) is used to measure the accuracy of the model.
η = n = 1 sample   number Prediction   output n   Actual   output n Actual   output n
The η for the cross-validation results of the proposed method shown in Figure 4 are 68.99 and 124.70 for h and PH, respectively. From the results, one can conclude that the model of level h is more accurate. There is an error in the PH model, which is because the nonlinearity of PH is stronger.
A control simulation based on the proposed method with the model obtained (Equation (28)) and the conventional MPC with model Equation (18) is made; the results are shown in Figure 5. OMPC denotes the proposed method. The control horizon is 2 and the predictive horizon is 9 for both methods. One can conclude that although the conventional MPC method guarantees stability, there is static error. The proposed method can track the setpoints without static error. Shown in Figure 6 is the control results of the conventional method in latent variable space. The solid line represents the output of the dynamic PLS model in latent variable space. The dashed line represents the projection of the actual system output in latent variable space. Compared with the results of conventional MPC in original space (Figure 5), one can conclude that conventional MPC can make the sub-model in latent variable space track the projection of setpoints in this space. However, due to the model mismatch, the control law projected back into original space cannot make the system output track the setpoints without static error.
To verify the control effect when the parameters of the actual system change or the model parameters are uncertain, the method proposed in this paper is compared with the predictive function control based on the adaptive dynamic fuzzy PLS model introduced in reference [41] (denoted as PFC). In the following simulation, the model parameter Cv is modified to simulate the actual system change. A 500-step simulation is performed. At the 100th step, Cv is changed from 8.79 to 9.79. The results are shown in Figure 7.
In the method used in reference [41], the parameters of the PLS model are identified and modified adaptively. This method will increase the computational burden and cause a large fluctuation of the system because it involves updating the internal and external models of the dynamic PLS model at the sample time of the parameter change and some subsequent sampling moments. The method proposed in this paper can still make the output of the system track the setpoints through its feedback when the model parameters of the system change and the fluctuation is small. For the nonlinear PH variable, there are static errors in both methods. The static error of the proposed method is 0.027, less than that of the method in reference [41], which is 0.055. This is because the proposed model is not reidentified when the model parameters changed.

6. Conclusions

In order to solve the problem of offset-free control for nonlinear systems, a nonlinear offset-free MPC based on a dynamic fuzzy PLS model is proposed.
The T–S fuzzy model is used as the inner model of PLS to describe the dynamic nonlinear characteristics of the MIMO system. Moreover, the disturbance model and observer are used to compensate the model mismatch during the MPC controller design. The simulation results of the PH neutralization titration show that the proposed method can achieve the offset-free tracking to the setpoint. It also achieves satisfactory control performance when external disturbances and model parameters are changed in the system.

Author Contributions

Conceptualization, X.J. and S.L.; methodology, Q.Z.; software, H.Y.; validation, Q.Z., X.J. and S.L.; formal analysis, Q.Z. and X.J.; investigation, H.Y.; writing—original draft preparation, Q.Z.; writing—review and editing, X.J. and H.Y.; visualization, H.Y.; supervision, S.L.; project administration, X.J. and S.L.; funding acquisition, Q.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China, grant number 62073158.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wu, C.; Blaabjerg, F. Advanced control of power electronic systems—An overview of methods. In Control of Power Electronic Converters and Systems; Elsevier: Amsterdam, The Nehterlands, 2021; pp. 1–33. [Google Scholar]
  2. Pereira, D.F.; Lopes, F.; Watanabe, E.H. Nonlinear model predictive control for the energy management of fuel cell hybrid electric vehicles in real-time. IEEE Trans. Ind. Electron. 2021, 68, 3213–3223. [Google Scholar] [CrossRef]
  3. Xu, X.; Simkoff, J.M.; Baldea, M.; Chiang, L.H.; Castillo, I.; Bindlish, R.; Ashcraft, B. Data-driven plant–model mismatch estimation for dynamic matrix control systems. Int. J. Robust Nonlinear Control 2020, 30, 7103–7129. [Google Scholar] [CrossRef]
  4. Yuan, X.; Zhang, S.; Zhang, C. Improved model predictive current control for SPMSM drives with parameter mismatch. IEEE Trans. Ind. Electron. 2019, 67, 852–862. [Google Scholar] [CrossRef]
  5. Gerksic, S.; Strmcnik, S.; Boom, T. Feedback action in predictive control: An experimental case study. Control Eng. Pract. 2008, 16, 321–332. [Google Scholar] [CrossRef]
  6. Jap, A.; Tlmsa, B.; Am, A. Mixed stochastic-deterministic tube MPC for offset-free tracking in the presence of plant-model mismatch. J. Process Control 2019, 83, 102–120. [Google Scholar]
  7. Wang, Z.; Zheng, Z.; Li, Y.; Sun, J. An offset-free robust model predictive control with incremental model and improved current observer for induction motor. Int. Trans. Electr. Energy Syst. 2019, 29, 1–18. [Google Scholar] [CrossRef]
  8. Awha, B.; Ssc, B. An offset-free MPC formulation for nonlinear systems using adaptive integral controller. ISA Trans. 2019, 91, 66–77. [Google Scholar]
  9. Hamdia, K.M.; Zhuang, X.; Rabczuk, T. An efficient optimization approach for designing machine learning models based on genetic algorithm. Neural Comput. Appl. 2021, 33, 1923–1933. [Google Scholar] [CrossRef]
  10. Benrabah, M.; Kara, K.; Aitsahed, O.; Hadjili, M.L. Constrained nonlinear predictive control using neural networks and teaching–learning-based optimization. J. Control Autom. Electr. Syst. 2021, 32, 1228–1243. [Google Scholar] [CrossRef]
  11. Wu, Z.; Rincon, D.; Christofides, P.D. Process structure-based recurrent neural network modeling for model predictive control of nonlinear processes. J. Process Control 2020, 89, 74–84. [Google Scholar] [CrossRef]
  12. Kong, L.; Yuan, J. Generalized Discrete-time nonlinear disturbance observer based fuzzy model predictive control for boiler-turbine systems. ISA Trans. 2019, 90, 89–106. [Google Scholar] [CrossRef] [PubMed]
  13. Escano, J.M.; Bordons, C.; Witheephanich, K.; Gomez-Estern, F. Fuzzy model predictive control: Complexity reduction for implementation in industrial systems. Int. J. Fuzzy Syst. 2019, 21, 2008–2020. [Google Scholar] [CrossRef]
  14. Han, H.; Liu, Z.; Hou, Y.; Qiao, J. Data-driven multiobjective predictive control for wastewater treatment process. IEEE Trans. Ind. Inform. 2020, 16, 2767–2775. [Google Scholar] [CrossRef]
  15. Metz, M.; Lesnoff, M.; Abdelghafour, F.; Akbarinia, R.; Roger, J.M. A “big-data” algorithm for KNN-PLS. Chemometr. Intell. Lab. Syst. 2020, 203, 104076. [Google Scholar] [CrossRef]
  16. Dong, Y.; Qin, S.J. Dynamic-inner partial least squares for dynamic data modeling. IFAC-PapersOnLine 2015, 48, 117–122. [Google Scholar] [CrossRef]
  17. Kaspar, M.H.; Ray, W.H. Dynamic PLS modelling for process control. Chem. Eng. Sci. 1993, 48, 3447–3461. [Google Scholar] [CrossRef]
  18. Fei, Z.; Liu, K.; Chi, Q.; Liang, J. Latent-variable nonlinear model predictive control strategy for a pH neutralization process. Asian J. Control Affil. ACPA Asian Control Profr. Assoc. 2015, 17, 2427–2434. [Google Scholar]
  19. Chen, J.; Cheng, Y.C.; Yea, Y. Multiloop PID controller design using partial least squares decoupling structure. Korean J. Chem. Eng. 2005, 22, 173–183. [Google Scholar] [CrossRef]
  20. Lü, Y.; Liang, J. Multi-loop constrained iterative model predictive control using ARX-PLS decoupling structure. Chin. J. Chem. Eng. 2013, 21, 1129–1143. [Google Scholar] [CrossRef]
  21. Zhang, J.; Lin, M.; Chen, J.; Xu, J.; Li, K. PLS-based multi-loop robust H-2 control for improvement of operating efficiency of waste heat energy conversion systems with organic Rankine cycle. Energy 2017, 123, 460–472. [Google Scholar] [CrossRef] [Green Version]
  22. Bang, Y.H.; Chang, K.Y.; Lee, I.B. Nonlinear PLS modeling with fuzzy inference system. Chemometr. Intell. Lab. Syst. 2002, 64, 137–155. [Google Scholar] [CrossRef]
  23. Abdel-Rahman, A.I.; Lim, G.J. A nonlinear partial least squares algorithm using quadratic fuzzy inference system. J. Chemometr. 2010, 23, 530–537. [Google Scholar] [CrossRef]
  24. Li, X.; Zhao, Z.; Liu, F. Latent variable iterative learning model predictive control for multivariable control of batch processes. J. Process Control 2020, 94, 1–11. [Google Scholar] [CrossRef]
  25. Gao, T.; Luo, H.; Yin, S.; Kaynak, O. A recursive modified partial least square aided data-driven predictive control with application to continuous stirred tank heater. J. Process Control 2020, 89, 108–118. [Google Scholar] [CrossRef]
  26. Zhao, Z.; Bin, H.U.; Liang, J. Multi-loop adaptive internal model control based on a dynamic partial least squares model. J. Zhejiang Univ. Sci. A Appl. Phys. Eng. 2011, 22, 190–200. [Google Scholar] [CrossRef]
  27. Chi, Q.; Liang, J. A multiple model predictive control strategy in the PLS framework. J. Process Control 2015, 25, 129–141. [Google Scholar] [CrossRef]
  28. Takagi, T.; Sugeno, M. Fuzzy identification of systems and its applications to modeling and control. Read. Fuzzy Sets Intell. Syst. 1993, 15, 387–403. [Google Scholar]
  29. Forgy, E.W. Cluster analysis of multivariate data: eIciency versus interpretability of classifications. Biometrics 1965, 21, 768–780. [Google Scholar]
  30. Bezdek, J.C. Pattern Recognition with Fuzzy Objective Function Algorithms; Springer: Berlin/Heidelberg, Germany, 1981. [Google Scholar]
  31. Gustafson, D.E.; Kessel, W.C. Fuzzy Clustering with a Fuzzy Covariance Matrix. In Proceedings of the 1978 IEEE Conference on Decision and Control Including the 17th Symposium on Adaptive Processes, San Diego, CA, USA, 10–12 January 1978. [Google Scholar]
  32. Cao, S.G.; Rees, N.W. Identification of dynamic fuzzy models. Fuzzy Sets Syst. 1995, 74, 307–320. [Google Scholar] [CrossRef]
  33. Barker, M.; Rayens, W. Partial least squares for discrimination. J. Chemometr. 2003, 17, 166–173. [Google Scholar] [CrossRef]
  34. Malthouse, E.C.; Tamhane, A.C.; Mah, R. Nonlinear partial least squares. Comput. Chem. Eng. 1997, 21, 875–890. [Google Scholar] [CrossRef]
  35. Qin, S.J.; Mcavoy, T.J. Nonlinear PLS modeling using neural networks. Comput. Chem. Eng. 1992, 16, 379–391. [Google Scholar] [CrossRef]
  36. Wang, H.O.; Tanaka, K.; Griffin, M.F. An approach to fuzzy control of nonlinear systems: Stability and design issues. IEEE Trans. Fuzzy Syst. 2002, 4, 14–23. [Google Scholar] [CrossRef]
  37. Zhang, T.; Gang, F.; Zeng, X.J. Output tracking of constrained nonlinear processes with offset-free input-to-state stable fuzzy predictive control. Automatica 2009, 45, 900–909. [Google Scholar] [CrossRef]
  38. Tanaka, K.; Ikeda, T.; Wang, H.O. Fuzzy regulators and fuzzy observers: Relaxed stability conditions and LMI-based designs. IEEE Trans. Fuzzy Syst. 2002, 6, 250–265. [Google Scholar] [CrossRef] [Green Version]
  39. Su, H. Basic Theory of Robust Control; Science Press: Beijing, China, 2010. [Google Scholar]
  40. Hall, R.C.; Seborg, D.E. Modelling and Self-Tuning Control of a Multivariable pH Neutralization Process Part I: Modelling and Multiloop Control. In Proceedings of the American Control Conference, Pittsburgh, PA, USA, 21–23 June 1989. [Google Scholar]
  41. Lü, Y. Modeling and Predictive Controller Design Based on the Dynamic PLS Method; Zhejiang University: Hangzhou, China, 2013. [Google Scholar]
Figure 1. Dynamic PLS control framework proposed by Kaspar and Ray.
Figure 1. Dynamic PLS control framework proposed by Kaspar and Ray.
Processes 09 01784 g001
Figure 2. A schematic illustration of the process.
Figure 2. A schematic illustration of the process.
Processes 09 01784 g002
Figure 3. The data set for identification.
Figure 3. The data set for identification.
Processes 09 01784 g003
Figure 4. Cross-validation results of the proposed method.
Figure 4. Cross-validation results of the proposed method.
Processes 09 01784 g004
Figure 5. Comparison of the two control methods.
Figure 5. Comparison of the two control methods.
Processes 09 01784 g005
Figure 6. Results of the conventional method in latent variable space.
Figure 6. Results of the conventional method in latent variable space.
Processes 09 01784 g006
Figure 7. Results under the model disturbance.
Figure 7. Results under the model disturbance.
Processes 09 01784 g007
Table 1. The parameters of the process.
Table 1. The parameters of the process.
A = 207 cm2.
Cv = 8.75 mL·cm−1/2·s−1
pK1 = 6.35
pK2 = 10.25
Wa1 = 3 × 10−3 M
Wa2 = −3 × 10−2 M
Wa3 = −3.05 × 10−3 M
Wb3 = −5.05 × 10−5 M
q1 = 16.6 mL·s−1
q2 = 0.55 mL·s−1
q3 = 15.6 mL·s−1
[acid] = 0.003 M HNO3
[buffer] = 0. 03 M NaHCO3
[base] = 0.003 M NaOH
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhao, Q.; Jin, X.; Yu, H.; Lu, S. Nonlinear Offset-Free Model Predictive Control based on Dynamic PLS Framework. Processes 2021, 9, 1784. https://doi.org/10.3390/pr9101784

AMA Style

Zhao Q, Jin X, Yu H, Lu S. Nonlinear Offset-Free Model Predictive Control based on Dynamic PLS Framework. Processes. 2021; 9(10):1784. https://doi.org/10.3390/pr9101784

Chicago/Turabian Style

Zhao, Qiang, Xin Jin, Huapeng Yu, and Shan Lu. 2021. "Nonlinear Offset-Free Model Predictive Control based on Dynamic PLS Framework" Processes 9, no. 10: 1784. https://doi.org/10.3390/pr9101784

APA Style

Zhao, Q., Jin, X., Yu, H., & Lu, S. (2021). Nonlinear Offset-Free Model Predictive Control based on Dynamic PLS Framework. Processes, 9(10), 1784. https://doi.org/10.3390/pr9101784

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