Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Recent Developments of the Lauricella String Scattering Amplitudes and Their Exact SL(K + 3,C) Symmetry
Next Article in Special Issue
Hardware in the Loop Topology for an Omnidirectional Mobile Robot Using Matlab in a Robot Operating System Environment
Previous Article in Journal
Monotone Iterative Technique for the Periodic Solutions of High-Order Delayed Differential Equations in Abstract Spaces
Previous Article in Special Issue
A Heuristic Approach to Shelf Space Allocation Decision Support Including Facings, Capping, and Nesting
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Least Squares Support Vector Machine-Based Multivariate Generalized Predictive Control for Parabolic Distributed Parameter Systems with Control Constraints

1
Department of Automation, Harbin University of Science and Technology, Harbin 150086, China
2
Key Laboratory of Advanced Manufacturing and Intelligent Technology, Ministry of Education, Harbin University of Science and Technology, Harbin 150086, China
3
School of Electrical Engineering, Computing and Mathematical Sciences, Curtin University, Perth, WA 6845, Australia
4
School of Mathematical Sciences, Sunway University, Selangor 47500, Malaysia
*
Author to whom correspondence should be addressed.
Symmetry 2021, 13(3), 453; https://doi.org/10.3390/sym13030453
Submission received: 13 February 2021 / Revised: 8 March 2021 / Accepted: 8 March 2021 / Published: 10 March 2021
(This article belongs to the Special Issue Symmetry in Optimization and Control with Real World Applications)

Abstract

:
This manuscript addresses a new multivariate generalized predictive control strategy using the least squares support vector machine for parabolic distributed parameter systems. First, a set of proper orthogonal decomposition-based spatial basis functions constructed from a carefully selected set of data is used in a Galerkin projection for the building of an approximate low-dimensional lumped parameter systems. Then, the temporal autoregressive exogenous model obtained by the least squares support vector machine is applied in the design of a multivariate generalized predictive control strategy. Finally, the effectiveness of the proposed multivariate generalized predictive control strategy is verified through a numerical simulation study on a typical diffusion-reaction process in radical symmetry.

1. Introduction

Modeling and control issues in the process industry involving distributed parameter systems (DPSs), which are infinite-dimensional systems, are challenging. The identification of unknown parameters of DPSs is much more complex than that of the lumped parameter systems (LPSs). In general, an infinite-dimensional problem is required to be approximated by finite-dimensional LPSs. Then, the model reduction methods are used to reduce the dimension of the approximate LPSs. For further details, see [1].
A system governed by parabolic partial differential equations (PDEs) is a typical DPS. It can be approximated by a set of low-dimensional ordinary differential equations through the use of Proper Orthogonal Decomposition (POD) technique [2], which is also called Karhunen-Loève expansion (K-L expansion) [3] or Principal Component Analysis (PCA) [4]. The POD technique is a popular technique for carrying out the task of model order reduction of parabolic systems. The spatial basis functions constructed from those slow modes of the system can capture the dominant dynamics while those fast modes can be ignored [5]. For the identification of the unknown parameters of the approximate LPSs constructed from the low-dimensional time-series, conventional approaches, such as projection transformation, can be used. The Singular Value Decomposition and K-L expansion is applied to construct a low-dimensional LPS with a simple structure. However, the accuracy of the approximation is low for a nonlinear DPS [6,7]. The K-L expansion combined with the neural network method can be applied to construct a sufficiently accurate approximation for a nonlinear DPS [8,9,10,11,12,13,14]. The support vector machine (SVM) is a general identification approach, and the least-squares support vector machine (LS-SVM) is its modified version. The LS-SVM, which simplifies the calculation by replacing inequality constraints with equality constraints, inherits the characteristics of the global optimization property and high generalization capabilities of SVM [15]. The model reduction approach based on the K-L expansion combined with LS-SVM is applied to regression problems that do not have a large amount of sample data. It is potentially applicable to a wide range of DPSs [16,17]. Several improved PCA methods are available in the literature, such as the incremental POD method [18], the adaptive POD method [19], the adaptive K-L/Fuzzy method [20], and the nonlinear-PCA method [21], for dealing with more complex nonlinear DPSs. For the nonlinear parabolic DPSs with time-dependent spatial domains, a temporally-local model order-reduction technique [22] and sparse POD-Galerkin methodology [23] are proposed and have been successfully applied to handle the moving boundary problem for a hydraulic fracturing process.
Generalized predictive control (GPC) strategy is introduced in [24,25]. Due to its simplicity for implementation, and excellent performance and robustness, it has been successfully implemented in a diversity of industrial applications from single-variable to multi-variable, and linear to nonlinear [26].
In this paper, we present a new multivariate GPC strategy used in conjunction with the POD technique to reduce the model order of weakly nonlinear DPSs of parabolic type. Then, the LS-SVM method is used to identify the unknown parameters of the low-dimensional time-series. Detailed methodology used to construct the generalized predictive temperature control in the tubular chemical reactor [27] will be discussed, and the control constraints and stability analysis will also be considered in the paper. The GPC based on SVM method has been successfully applied to LPSs with nonlinear characteristics [28,29,30]. The predictor obtained using the SVM method combined with the POD technique is found to have higher accuracy. It is applicable to a wider range of nonlinear DPSs. On the other hand, the GPC strategy used in conjunction with POD and recursive least square method has been applied to DPSs [31,32,33]. Thus, it is advantageous to combine GPC with SVM for the modeling and control of DPSs.
The remainder of this paper is structured as follows. First, the system equations are described in Section 2. Then, the POD and LS-SVM for regression are briefly introduced, and the POD and LS-SVM-based predictive model is derived in Section 3. The POD and LS-SVM-based multivariate GPC strategy is presented in Section 4. In Section 5, the numerical simulation of a diffusion-reaction process in radical symmetry is carried out to examine the effectiveness of the proposed algorithm in Section 3. The summary and future work of our paper are included in Section 6.

2. System Description

We consider a class of weakly nonlinear PDEs of parabolic type as given below:
y ( x , t ) t = A 2 y ( x , t ) 2 x + B u ( t ) + f ( x , t , y ( x , t ) , y ( x , t ) x ) ,
where t is the time variable, x Ω is the spatial variable, and Ω is the spatial domain, y ( x , t ) is the controlled output variable, u ( t ) is the control input variable, A and B are the linear operators in Hilbert space H , and f ( · ) is an unknown weakly nonlinear function. System (1) is applicable to model a variety of physical and biochemical processes, such as catalytic reaction rods, steel casting, and the snap oven [5].
For system (1), it contains an unknown weakly nonlinear function f ( · ) . To obtain accurate information of such a system, a sufficient number of sensors need to be placed along the spatial location. Due to actual physical conditions, only a small number of actuators is allowed to be mounted to observe the state. The input-output data are obtained from the actual production process under random signal excitation. The control algorithm is developed in three stages: In Section 3, the output { y ( x m , t n ) } m = 1 , n = 1 M , N excited by the input { u ( k , t n ) } k = 1 , n = 1 K , N are used by the POD technique to reduce the dimensionality of the approximate model, and M, N, and K are, respectively, the number of the sensors, the sample duration, and the actuators. A low-dimensional time-series obtained after model order reduction is identified by LS-SVM. In Section 4, the POD used in conjunction with the LS-SVM-based multivariate GPC strategy is utilized to design the required controller through solving a receding horizon optimization problem at every sample period.

3. Model Reduction by POD and LS-SVM

The POD and LS-SVM-based model will be established in this section. Since a high dimensional output is not ideal for the design of the controller, the POD technique will be applied first to reduce the dimension of the model. Let the input u ( · , t ) K × N , and let y ( x m , t n ) M × N . The POD basis of rank can be calculated by Algorithm 1 proposed in [2].
Algorithm 1 Proper Orthogonal Decomposition (POD) basis of rank
Require: Snapshots { y ( · , t ) } j = 1 M M , P O D r a n k M , where is the rank of the POD basis.
  • Step 1. Set Y = [ y 1 , , y M ] T M × N .
  • Step 2. Determine R = Y Y T M × M .
  • Step 3. Compute eigenvalue decomposition [ Ψ , Λ ] = e i g ( R ) .
  • Step 4. Set ψ i = Ψ · , i M , a n d λ i = Λ i i , λ 1 λ 2 λ λ M , f o r i = 1 , , M .
  • Step 5. Compute E ( ) = i = 1 λ i / i = 1 M λ M .
  • return POD basis Ψ = { ψ i } i = 1 and eigenvalues { λ i } i = 1 .
Choose the minimum such that E ( ) 99 % . Then,
y ( x , t ) = Ψ ( x ) Y ( t ) M × N ,
Y ( t ) = Ψ T ( x ) y ( x , t ) × N ,
where Ψ ( i ) ( x ) and Y ( i ) ( t ) , i = 1 , , , represent, respectively, the i-th basis function and the corresponding low-dimensional time-series, Y ( t ) = [ Y ( 1 ) ( t ) , , Y ( ) ( t ) ] T . Due to the orthogonality of Ψ ( i ) ( x ) , the Y ( i ) ( t ) are independent. Through the use of the LS-SVM approach proposed in [15], the i-th estimated value of the time-series is:
Y ^ ( i ) ( t n ) = w ( i ) T ϕ ( i ) ( X ( i ) ( t n ) ) + b ( i ) ,
where X ( i ) ( t n ) = [ Y ( i ) ( t n 1 ) , , Y ( i ) ( t n n y ) , u 1 ( t n 1 ) , , u 1 ( t n n u ) , , u M ( t n 1 ) , , u M ( t n n u ) ] T , i = 1 , , , w ( i ) is the scale factor, b ( i ) is the bias, ϕ ( i ) ( · ) is a function which maps the vector X ( i ) from function space into feature space. The optimization problem can be expressed as:
m i n w ( i ) , b ( i ) , ξ ( i ) J ( w ( i ) , ξ ( i ) ) = 1 2 w ( i ) T w ( i ) + 1 2 γ ( i ) j = 1 N ξ ( i ) 2 , s . t . Y ( i ) ( t n ) = w ( i ) T ϕ ( X ( i ) ( t n ) ) + b ( i ) + ξ ( i ) , j = 1 , 2 , , N ,
where γ ( i ) is the regularization factor, and ξ ( i ) is the relaxation factor. According to Lagrange multiplier method, the Lagrange function is:
L ( w ( i ) , b ( i ) , ξ ( i ) , α ( i ) ) = 1 2 w ( i ) T w ( i ) + 1 2 γ ( i ) j = 1 N ξ ( i ) 2 j = 1 N α ( i ) ( j ) ( w ( i ) T ϕ ( X ( i ) ) + b ( i ) + ξ ( i ) Y i ) ,
where α ( i ) is the Lagrange operator. Let L / w ( i ) = L / b ( i ) = L / e ( i ) = L / α ( i ) = 0 , the coefficients α ( i ) and b ( i ) are obtained through solving the following linear algebraic equations:
0 1 1 1 X ( i ) ( t 1 ) T · X ( i ) ( t 1 ) + 1 γ ( i ) X ( i ) ( t 1 ) T · X ( i ) ( t N ) 1 X ( i ) ( t N ) T · X ( i ) ( t 1 ) X ( i ) ( t N ) T · X ( i ) ( t N ) + 1 γ ( i ) b ( i ) α 1 ( i ) α N ( i ) = 0 Y ( i ) ( t 1 ) Y ( i ) ( t N ) ,
The time-series Y ( i ) can be obtained as:
Y ^ ( i ) ( t n ) = τ = 1 N α τ ( i ) · K ( i ) X ( i ) ( t n ) , X ( i ) ( τ ) + b ( i ) = j = 1 n y τ = 1 N α τ ( i ) X ( i ) ( τ ) [ j ] z j · Y ( i ) ( t n ) + k = 1 K j = 1 n u τ = 1 N α τ ( i ) X ( i ) ( τ ) [ n y + ( k 1 ) · n u + j ] z j · u k ( t n ) + b ( i ) ,
where K ( i ) ( · , · ) is the kernel function. The linear kernel function K ( i ) X ( i ) ( t n ) , X ( i ) ( τ ) = X ( i ) ( t n ) T · X ( i ) ( τ ) is chosen, and X ( i ) ( t n ) [ j ] is the j-th term of X ( i ) ( t n ) , z j is the j-step backward shift operator.
The matrix form of the time-series autoregressive exogenous (ARX) model obtained by the LS-SVM approach is:
A ( z 1 ) Y ( t ) = B ( z 1 ) U ( t 1 ) + b ,
where A ( z 1 ) = I A , A = d i a g ( A ( 1 ) , , A ( ) ) , A ( i ) = j = 1 n y τ = 1 N α τ ( i ) X ( i ) ( τ ) [ j ] z j , B ( z 1 ) = [ B ( 1 ) , , B ( ) ] , B ( i ) = [ B 1 ( i ) , , B K ( i ) ] , B k ( i ) = k = 1 K j = 1 n u τ = 1 N α τ ( i ) X ( i ) ( τ ) [ n y + ( k 1 ) · n u + j ] z j , i = 1 , , , and b = [ b ( 1 ) , , b ( ) ] T , k = 1 , , K .

4. POD and LS-SVM-Based Multivariate GPC

GPC strategy proposed in [23,24] is a class of predictive control algorithm developed for adaptive control due to its excellent performance and simplicity for implementation. It has been applied in a diversity of actual industrial plants during the past three decades. A multivariate GPC strategy will be proposed for parabolic DPSs with unknown coefficients in this section.

4.1. Multivariate GPC Strategy

Consider the following ARX model containing a zero-mean white noise:
A ( z 1 ) Y ( t ) = B ( z 1 ) U ( t ) + C ( z 1 ) e ( t ) Δ ,
where the output Y ( t ) × 1 , the input U ( t ) K × 1 , and the white noise with zero mean e ( t ) × 1 . Δ = 1 z 1 , C ( z 1 ) is the identity matrix. Minimizing the quadratic cost function on a finite horizon to obtain the optimal control increment as follows:
J ( N p , N u ) = j = 1 N p Y ^ ( t + j | t ) w ( t + j ) 2 + j = 1 N u Δ u ( t + j 1 ) λ 2 ,
where Y ^ ( t + j | t ) is the output predictive j-step ahead value based on previous state up to time t. λ is the weighted coefficient of Δ u . N p and N u are, respectively, the prediction horizons and control horizons. The reference trajectory w ( t + j ) is defined as:
w ( t + j ) = α w w ( t + j 1 ) + ( 1 α w ) Y s p ( t + j ) ,
Y s p ( t + j ) = Ψ ( x ) y s p ( x , t + j ) ,
where y s p is the set point value. Consider the Diophantine equation:
I × = E j ( z 1 ) A ˜ + z j F j ( z 1 ) ,
where A ˜ = A Δ . The order of the polynomial matrix E j ( z 1 ) is j 1 , and that of the polynomial matrix F j ( z 1 ) is n y . Since b is constant, Δ b is zero.
Let E j ( z 1 ) B ( z 1 ) = G j ( z 1 ) + z j G j p ( z 1 ) and f j = G j p ( z 1 ) Δ u ( t 1 ) + F j ( z 1 ) Y ( t ) , the j-step prediction is:
Y ^ ( t + j | t ) = G j ( z 1 ) Δ u ( t + j 1 ) + f j ,
By solving E j ( z 1 ) recursively, the N p -step prediction value Y ^ can be written as:
Y ^ ( t + 1 | t ) Y ^ ( t + 2 | t ) Y ^ ( t + N u | t ) Y ^ ( t + N p | t ) = G 0 0 0 G 1 G 0 0 G N u 0 G 0 G N p 1 G N p 2 G N p N u Δ u ( t ) Δ u ( t + 1 ) Δ u ( t + N u 1 ) + f 1 f 2 f N u f N p ,
where G j ( z 1 ) = i = 0 j 1 G i z i . The matrix form of the prediction value Y ^ is:
Y = G Δ u + f .
Rewrite the quadratic criterion (11) as:
J = ( G Δ u + f w ) T ( G Δ u + f w ) + Δ u T λ Δ u ,
If there are no constraints, the optimum Δ u without constraint can be obtained as:
Δ u = H G T ( w f ) ,
where H = ( G T G + λ I ) 1 . It is difficult to avoid control constraints in practical systems. The constrained incremental control law directly truncated at the boundary is not optimal. Considering the control constraints Θ Δ u θ , the quadratic criterion using a Lagrange multiplier is given as follows [31]:
J ( N p , N u ) = j = 1 N p Y ^ ( t + j | t ) w ( t + j ) 2 + j = 1 N u Δ u ( t + j 1 ) λ 2 + Θ Δ u θ ε 2 ,
where ε is the weighting marix of control constraint which is positive definite. The optimum Δ u ˜ with constraint can be obtained as:
Δ u ˜ = Δ u H Θ T ( Θ H Θ T + ε 1 ) 1 ( Θ Δ u θ ) ,
The quadratic criterion is solved repeatedly at the subsequent sampling times, only the first term of Δ u ˜ is applied to the actual system, so the control action is:
u ( k , t j ) = u ( k , t j 1 ) + Δ u ˜ ( k , t j ) ,
where k, k = 1 , , K , is defined in Section 2. The constrained control action is handled by the truncated algorithm:
u ( k , t j ) = u m i n , u ( k , t j ) u m i n u ( k , t j ) , u m i n u ( k , t j ) u m a x u m a x , u ( k , t j ) u m a x ,
where u m i n u ( k , t j ) u m a x .
The proposed POD and LS-SVM-based multivariate GPC algorithm can be summarized in Algorithm 2.
Algorithm 2 The POD and least-squares support vector machine (LS-SVM)-based multivariate generalized predictive control (GPC).
Require: A set of output y ( x m , t n ) is derived by appropriate excitation signals.
  • Step 1. POD basis { ψ i } i = 1 and eigenvalues { λ i } i = 1 are obtained by Algorithm 1.
  • Step 2. The predictor Y ^ and the ARX model are built by LS-SVM with a linear kernel.
  • Step 3. Compute Δ u by Equation (16). If the constrained condition is not satisfied, replaced Δ u with Δ u ˜ using Equation (18).
  • Step 4. Use the control law u ( k , t j ) to the actual system at instant t j .
  • Step 5. Let j = : j + 1 , Repeat Step 3 and Step 4 at instant t j + 1 .

4.2. The Stability Analysis

We assumed that the following assumptions are satisfied throughout:
Assumption 1.
The reference trajectory w ( t + j ) is bounded, and the prediction error of model (8) is bounded.
Assumption 2.
System (1) is open-loop stable, so that the influence of the control action exerted at the past time tends to zero with the increase of time.
Assumption 3.
The fast modes of the dynamic of system (1) are ignored, so that it suffices to consider the stability conditions of the nominal system (9).
Assumption 4.
The future disturbances are equal to that at time t. The mismatch and any uncertainty of model (15) can be estimated as:
D ^ ( t + 1 ) = [ d ^ ( t + 1 ) T d ^ ( t + N p ) T ] T ,
d ^ ( t + 1 ) = ϵ ( Y ( t + 1 ) Y ^ ( t + 1 ) ) ,
where Y ( t + 1 ) = Ψ T ( x ) y ( x , t + 1 ) , and 0 ϵ 1 .
The real value of the nominal system (10) at time t + 1 are
Y ( t + 1 ) = Y ^ ( t + 1 ) + d ^ ( t + 1 ) = G 0 Δ u ( t ) + f 1 + ϵ ( Y ( t ) Y ^ ( t ) ) .
Assume that only Y ( t + 1 ) and u ( t ) are included in the closed-loop system. Then, Y ˜ ( t + 1 ) = [ Y ( t + 1 ) Δ u ( t ) ] . Let D = H G T , Δ u ( t ) = D 1 F Y ( t ) + D 1 where D 1 is the first columns of D . Substitute Δ u ( t ) into Equation (26). Then,
Y ˜ ( t + 1 ) = A + B D 1 F 0 D 1 F I Y ˜ ( t ) + D 1 D 1 ω + e ( t ) 0 .
Lemma 1.
(Unconstrained case) [6] If the eigenvalue λ i of the matrix A satisfy | λ i | < 1 , then the closed-looped system is asymptotically stable, where
A ˜ = A + B D 1 F 0 D 1 F I .
Proof. 
By Assumptions 1 and 2, the output Y ˜ ( t + 1 ) = A ˜ Y ˜ ( t ) + f 1 ( ω ) + f 2 ( e ( t ) ) is bounded. Thus, the closed-looped system is asymptotically stable if that max | λ i | < 1 and f 1 ( ω ) + f 2 ( e ( t ) ) δ , where δ is a sufficient small number. By Assumption 3, system (1) is stable if the nominal system (10) is stable. □
Let D ˜ = ( I H Θ T ( Θ H Θ 1 + ϵ 1 ) 1 Θ ) H G T , Δ u ( t ) ˜ = D 1 ˜ F Y ( t ) + D 1 ˜ where D 1 ˜ is the first columns of D ˜ .
Lemma 2.
(Constrained case) [6] If the eigenvalue λ ˜ i of the matrix A ˜ satisfy | λ ˜ i | < 1 , then the closed-looped system is asymptotically stable, where
A ˜ = A + B D 1 ˜ F 0 D 1 ˜ F I .
Proof. 
The proof is similar to that given for Lemma 1. □

5. Case Study

The proposed multivariate GPC strategy is applied in the simulation study of an exothermic catalytic reaction in a long thin rod in radical symmetry, which has the form A B [34] of transport reaction as shown in Figure 1.
The goal is to produce as much product B as possible while consuming the least amount of energy by adjusting the axial cooling medium temperature. Assume that the density and the heat capacity of the rod are constant, and that the reaction rate of the rod is uniform. Due to radial symmetry, the dynamical distribution of the axial temperature on the rod is governed by the following PDE of parabolic type:
T t = 2 T 2 x + β T e T 1 + T e γ ¯ + β u ( T c T ) ,
subject to the first kind of boundary conditions, and initial condition:
T ( 0 , t ) = 0 , T ( L , t ) = 0 , T ( x , 0 ) = T 0 ( x ) ,
where t, x, and L are, respectively, the dimensionless time, axial position, and the length of reactor, and the controlled variables T ( x , t ) and the control variables T c ( x , t ) are, respectively, the axial temperature and the cooling medium temperature of the reactor. In the simulation, let the length of the rod L = π , the initial temperature T 0 ( x ) = 0.5 , the heat of reaction β T = 50 , the heat transfer coefficient β u = 2 , and the activation energy γ ¯ = 4 . The 21 temperature sensors and 4 cooling medium temperatures T c ( x , t ) are equally located along the long rod. Each sample period Δ t = 0.01 , and the total period is t s = 2 . All the actual simulation values are calculated using Matlab PDE toolbox. The random excitation input signal that varies randomly between 4 to 1 are shown in Figure 2 and the corresponding output is shown in Figure 3.
From Algorithm 1, the minimum number of the POD basis for ensuring E ( ) 99 % is 4. The first four basis functions ϕ ( x ) are shown in Figure 4.
The first 160 groups of data are used as training set, while the latter 40 groups of data are used as testing set. Then, the low-dimensional time-series Y ^ T ( i ) ( t n ) are identified by LS-SVM:
Y ^ T ( i ) ( t n ) = j = 1 n y τ = 1 N α τ ( i ) X T ( i ) ( τ ) [ j ] z j · Y T ( i ) ( t n ) + k = 1 K j = 1 n u τ = 1 N α τ ( i ) X T ( i ) ( τ ) [ n y + ( k 1 ) · n u + j ] z j · u k ( t n ) + b ( i ) , ( 31 )
Different structures of the low-dimensional time-series model are considered for our numerical experiments, and the orders of n y = 4 and n u = 3 are finally selected for which the desired accuracy of the corresponding reduced dimensional model is achieved. The parameters N = 160 , and K = 4 . The 10-fold cross-validation is used to obtained the regularization factor γ . The LS-SVM-based model can approximate the low-dimensional time-series satisfactorily which is shown in Figure 5. The error of the predictive model is between [ 0.02 , 0.015 ] which is shown in Figure 6, and the mean of the root mean square error is 0.0001 . The average training times of the model are [16.4, 5.5, 20.5, 2.8] ms, respectively.
For expressions (11) and (12), choose λ = 1 , α w = 0.3 , N p = 3 , N u = 2 as the parameters of the controller. The unconstrained control increment and control action are calculated by expressions (19) and (22), and the first column of the control variable T c is used in the actual system. The process is repeated at every sample period. The control variables T c and the output response curves for the unconstrained case are shown in Figure 7 and Figure 8, respectively.
In Figure 7 and Figure 8, the controller can quickly drive the system dynamics back to the steady-state in 30 sample periods, and the steady-state error is below 1 × 10 5 . From the simulation results, we can see that the proposed controller has excellent tracking performance and can be effectively applied to parabolic DPSs involving unknown coefficients. The fluctuation range of control variables between 8 to 1.5 is large.
We now consider the case of control constraints. More specifically, 0.5 Δ u 0.5 , and 4 u 1 . The constrained control increment and control action are updated by expressions (21) and (23). The control variables T c and the output response curves under the constrained case are shown in Figure 9 and Figure 10, respectively.
In Figure 9 and Figure 10, we can see that the controller is working well in the case of control constraints. Due to the presence of the control constraints, the settling time is now extended to 50 sample periods. The time taken for control action calculation in each sampling period is given in Figure 11.
In Figure 11, we can see that the time taken for unconstrained case is in milliseconds, and it is a bit longer for the constrained case.

6. Conclusions

A new control strategy, named POD and LS-SVM-based multivariate GPC, was developed for weakly nonlinear parabolic DPSs with unknown coefficients. The POD was used for model order reduction, resulting in low-dimensional lumped parameter systems. The desired control was designed based on this low-dimensional model. The LS-SVM with linear kernel function is found effective in the identification issues for linear and weakly nonlinear plants. An interesting question to ask is: Are there kernel functions for general nonlinear systems? The POD used in conjunction with the LS-SVM-based multivariate GPC strategy is easy to implement and is applicable to weakly nonlinear parabolic systems with unknown coefficients. As an illustration, a diffusion-reaction process was considered and solved using the proposed methods. The controller obtained has achieved the desired good performance. For future research, the online LS-SVM and radial basis kernel function will be applied to improve the modeling accuracy and to a more complex problem with nonlinearity. Furthermore, the method will also be extended to handle moving boundary problems, where the single POD basis function is replaced by a temporally-local basis function.

Author Contributions

Conceptualization, L.A. and K.L.T.; formal analysis, L.A.; funding acquisition, L.D.; investigation, L.A.; methodology, L.A. and K.L.T.; project administration, L.D. and K.L.T.; Software, L.A.; Supervision, K.L.T.; Validation, L.A.; visualization, Y.X.; writing—original draft, L.A. and Y.X.; writing—review and editing, L.A., L.D. and K.L.T.. All authors have read and agreed to the published version of the manuscript.

Funding

Part of this project is funded by National Natural Science Foundation of China No. 61673141, 61806060, Chinese Scholarship Council (CSC) under grant 201608230319, Natural Science Foundation of Heilongjiang Province(LH2019F024), Harbin Program of Science and Technology (2017RAQXJ006).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The author is also thankful to National Natural Science Foundation of China No. 61673141, 61806060, Chinese Scholarship Council (CSC) under grant 201608230319, Natural Science Foundation of Heilongjiang Province(LH2019F024), Harbin Program of Science and Technology (2017RAQXJ006).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
ARXAutoregressive exogenous
DPSDistributed parameter systems
GPCGeneralized predictive control
K-LKarhunen-Loève
LPSLumped parameter systems
LSLeast squares
PCAPrincipal component analysis
PDEPartial differential equation
PODProper orthogonal decomposition
SVMSupport vector machine

References

  1. Li, H.X.; Qi, C.K. Spatio-Temporal Modeling of Nonlinear Distributed Parameter Systems—A Time/Space Separation Based Approach; Springer: New York, NY, USA, 2007. [Google Scholar]
  2. Volkwein, S. Lecture: Proper Orthogonal Decomposition: Theory and Reduced-Order Modelling; Universtity of Konstanz: Konstanz, Germany, 2013. [Google Scholar]
  3. Mallat, S. A Wavelet Tour of Signal Processing-The Sparse Way, 3rd ed.; Elsevier: Singapore, 2009. [Google Scholar]
  4. Jolliffe, I. Principal Component Analysis, 2nd ed.; Springer: New York, NY, USA, 2002. [Google Scholar]
  5. Smyshlyaev, A.; Krstic, M. Adaptive Control of Parabolic PDEs; Princeton University Press: London, UK, 2010. [Google Scholar]
  6. Zheng, D.; Hoo, K.A. System identification and model-based control for distributed parameter systems. Comput. Chem. Eng. 2004, 28, 1361–1375. [Google Scholar] [CrossRef]
  7. Jing, X.T.; Li, H.X.; Lei, J.; Yang, N. KL-SVD based One-dimensional Modal Control combining with Diameter Control Subsystem as State Estimator in Czochralski Process. In Proceedings of the 9th International Conference on Intelligent Human-Machine Systems and Cybernetics, Hangzhou, China, 26–27 August 2017; pp. 189–192. [Google Scholar]
  8. Aggelogiannaki, E.; Sarimveis, H. Nonlinear model predictive control for distributed parameter systems using data driven artificial neural network models. Comput. Chem. Eng. 2008, 32, 1225–1237. [Google Scholar] [CrossRef]
  9. Li, H.X.; Deng, H.; Zhong, J. Model-based integration of control and supervision for one kind of curing process. IEEE Trans. Electron. Pack. Manuf. 2004, 27, 177–186. [Google Scholar] [CrossRef]
  10. Deng, H.; Li, H.X.; Chen, J.R. Delay-robust stabilization of a hyperbolic PDE-ODE system. IEEE Trans. Contr. Syst. Technol. 2005, 13, 686–700. [Google Scholar]
  11. Romijn, R.; Özkan, L.; Weiland, S.; Ludlage, J.; Marquardt, W. A grey-box modeling approach for the reduction of nonlinear systems. J. Process. Contr. 2008, 18, 906–914. [Google Scholar] [CrossRef]
  12. Bartecki, K. PCA-based approximation of a class of distributed parameter systems classical vs neural network approach. Bull. Pol. Acad. Sci. 2012, 60, 651–660. [Google Scholar] [CrossRef] [Green Version]
  13. Wang, M.L.; Shi, H.B. Spatiotemporal prediction for nonlinear parabolic distributed parameter system using an artificial neural network trained by group search optimization. Neurocomputing 2013, 113, 234–240. [Google Scholar] [CrossRef]
  14. Wang, M.L.; Shi, H.B. An adaptive neural network prediction for nonlinear parabolic distributed parameter system based on block-wise moving window technique. Neurocomputing 2014, 133, 67–73. [Google Scholar] [CrossRef]
  15. Suykens, J.A.K.; Vandewalle, J. Least squares support vector machine classifiers. Neural Process. Lett. 1999, 6, 293–300. [Google Scholar] [CrossRef]
  16. Qi, C.K.; Li, H.X. A LSSVM Modeling Approach for Nolinear Distributed Parameter Processes. In Proceedings of the 7th World Congress on Intelligent Control and Automation, Chongqing, China, 25–27 June 2008; pp. 569–574. [Google Scholar]
  17. Qi, C.K.; Li, H.X.; Zhang, X.X.; Zhao, X.C.; Li, S.Y.; Gao, F. Time/Space-Separation-Based SVM Modeling for Nonlinear Distributed Parameter Processes. Ind. Eng. Chem. Res. 2011, 50, 332–341. [Google Scholar] [CrossRef]
  18. Xu, C.; Schuster, E. Model Order Reduction for High Dimensional Linear Systems based on Rank-1 Incremental Proper Orthogonal Decomposition. In Proceedings of the 2011 American Control Conference, San Francisco, CA, USA, 29 June–1 July 2011; pp. 2975–2981. [Google Scholar]
  19. Pourkargar, D.B.; Armaou, A. APOD Based Control of Linear Distributed Parameter Systems Under Sensor Controller Communication Bandwidth Limitations. AIChE J. 2013, 61, 434–447. [Google Scholar] [CrossRef]
  20. Lu, X.J.; Zou, W.; Huang, M.H. An adaptive modeling method for time-varying distributed parameter processes with curing process applications. Nonlinear Dynam. 2015, 82, 865–876. [Google Scholar] [CrossRef]
  21. Wang, M.L.; Qi, C.K.; Yan, H.C.; Shi, H.B. Hybrid neural network predictor for distributed parameter system based on nonlinear dimension reduction. Neurocomputing 2016, 171, 1591–1597. [Google Scholar] [CrossRef]
  22. Narasingam, A.; Siddhamshetty, P.; Kwon, J.S. Temporal clustering for order reduction of nonlinear parabolic PDE systems with time-dependent spatial domains: Application to a hydraulic fracturing process. AIChE J. 2017, 63, 3818–3831. [Google Scholar] [CrossRef]
  23. Sidhu, H.S.; Narasingam, A.; Siddhamshetty, P.; Kwon, J.S. Mohtadi and P. S. Tuffs, Model order reduction of nonlinear parabolic PDE systems with moving boundaries using sparse proper orthogonal decomposition: Application to hydraulic fracturing. Comput. Chem. Eng. 2018, 112, 92–100. [Google Scholar] [CrossRef] [Green Version]
  24. Clarke, D.W.; Mohtadi, C.; Tuffs, P.S. Generalized predictive control–Part I:the basic algorithm. Automatica 1987, 23, 137–148. [Google Scholar] [CrossRef]
  25. Clarke, D.W.; Mohtadi, C.; Tuffs, P.S. Generalized predictive control–Part II:extensions and interpetations. Automatica 1987, 23, 149–160. [Google Scholar] [CrossRef]
  26. Camacho, E.F.; Bordons, C. Model Predictive Control, 2nd ed.; Springer: New York, NY, USA, 2007. [Google Scholar]
  27. Ai, L.; Zhang, D.S.; Teo, K.L.; Deng, L.W. Generalized Predictive Temperature Control in Tubular Chemical Reactors by means of POD and LSSVM. In Proceedings of the 39th Chinese Control Conference, Shenyang, China, 27–29 July 2020; pp. 2407–2412. [Google Scholar]
  28. Iplikci, S. Support vector machines-based generalized predictive control. Int. J. Robust. Nonlin. 2006, 16, 843–862. [Google Scholar] [CrossRef]
  29. Iplikci, S. Online trained support vector machines-based generalized predictive control of non-linear systems. Int. J. Adapt. Control. 2006, 20, 599–621. [Google Scholar] [CrossRef]
  30. Guo, Z.K.; Guan, X.P. Nonlinear generalized predictive control based on online least squares support vector machines. Nonlinear Dynam. 2015, 79, 1163–1168. [Google Scholar] [CrossRef]
  31. Li, N.; Hua, C.; Wang, H.F.; Li, S.Y.; Ge, S.S. Time-Space Decomposition Based Generalized Predictive Control of a Transport-Reaction Process. Ind. Eng. Chem. Res. 2011, 50, 11628–11635. [Google Scholar] [CrossRef]
  32. Hua, C.; Li, N.; Li, S.Y. Time-space ARX modeling and predictive control for distributed parameter system. Control Theory Appl. 2011, 28, 1711–1716. [Google Scholar]
  33. Liu, P.F.; Tao, J.L.; Wang, X.G. Decoupling ARX Time-space model based GPC for distributed parameter system. In Proceedings of the 34th Chinese Control Conference, Hangzhou, China, 28–30 July 2015; pp. 612–617. [Google Scholar]
  34. Christofides, P.D. Nonlinear and Robust Control of PDE Systems: Methods and Applications to Transport-Reaction Processes; Birkhäuser: Boston, MA, USA, 2001. [Google Scholar]
Figure 1. The schematic of the tubular reactor.
Figure 1. The schematic of the tubular reactor.
Symmetry 13 00453 g001
Figure 2. The random excitation input signal.
Figure 2. The random excitation input signal.
Symmetry 13 00453 g002
Figure 3. The corresponding output.
Figure 3. The corresponding output.
Symmetry 13 00453 g003
Figure 4. The first four Proper Orthogonal Decomposition (POD) basis functions.
Figure 4. The first four Proper Orthogonal Decomposition (POD) basis functions.
Symmetry 13 00453 g004
Figure 5. The performance of the least-squares support vector machine (LS-SVM)-based model.
Figure 5. The performance of the least-squares support vector machine (LS-SVM)-based model.
Symmetry 13 00453 g005
Figure 6. The prediction error of the POD used in conjunction with LS-SVM based model.
Figure 6. The prediction error of the POD used in conjunction with LS-SVM based model.
Symmetry 13 00453 g006
Figure 7. The control variables T c for the unconstrained case.
Figure 7. The control variables T c for the unconstrained case.
Symmetry 13 00453 g007
Figure 8. The output response curves for the unconstrained case.
Figure 8. The output response curves for the unconstrained case.
Symmetry 13 00453 g008
Figure 9. The control variables T c under the constrained case.
Figure 9. The control variables T c under the constrained case.
Symmetry 13 00453 g009
Figure 10. The output response curves under the constrained case.
Figure 10. The output response curves under the constrained case.
Symmetry 13 00453 g010
Figure 11. The time consuming of control law calculation.
Figure 11. The time consuming of control law calculation.
Symmetry 13 00453 g011
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ai, L.; Xu, Y.; Deng, L.; Teo, K.L. Least Squares Support Vector Machine-Based Multivariate Generalized Predictive Control for Parabolic Distributed Parameter Systems with Control Constraints. Symmetry 2021, 13, 453. https://doi.org/10.3390/sym13030453

AMA Style

Ai L, Xu Y, Deng L, Teo KL. Least Squares Support Vector Machine-Based Multivariate Generalized Predictive Control for Parabolic Distributed Parameter Systems with Control Constraints. Symmetry. 2021; 13(3):453. https://doi.org/10.3390/sym13030453

Chicago/Turabian Style

Ai, Ling, Yang Xu, Liwei Deng, and Kok Lay Teo. 2021. "Least Squares Support Vector Machine-Based Multivariate Generalized Predictive Control for Parabolic Distributed Parameter Systems with Control Constraints" Symmetry 13, no. 3: 453. https://doi.org/10.3390/sym13030453

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