Mathematical Problems in Engineering

Volume 2013, Article ID 104702, 6 pages

http://dx.doi.org/10.1155/2013/104702

## Model Predictive Control Based on Kalman Filter for Constrained Hammerstein-Wiener Systems

^{1}Institute of Advanced Control Technology, Dalian University of Technology, Dalian 116024, China^{2}Dalian Institute of Science and Technology, Dalian 116052, China

Received 6 August 2013; Accepted 22 September 2013

Academic Editor: William Guo

Copyright © 2013 Man Hong and Shao Cheng. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

To precisely track the reactor temperature in the entire working condition, the constrained Hammerstein-Wiener model describing nonlinear chemical processes such as in the continuous stirred tank reactor (CSTR) is proposed. A predictive control algorithm based on the Kalman filter for constrained Hammerstein-Wiener systems is designed. An output feedback control law regarding the linear subsystem is derived by state observation. The size of reaction heat produced and its influence on the output are evaluated by the Kalman filter. The observation and evaluation results are calculated by the multistep predictive approach. Actual control variables are computed while considering the constraints of the optimal control problem in a finite horizon through the receding horizon. The simulation example of the CSTR tester shows the effectiveness and feasibility of the proposed algorithm.

#### 1. Introduction

Temperature control of the batch reactor during chemical processes [1] is the key to producing qualified chemical products. The control requirement is to change the temperature in the reactor with the given temperature process curve set up in advance for the entire batch reaction process. The accuracy of following the temperature change is the scale used to judge the performance of the designed controller. Two factors affect temperature change. First is the change of the interlayer cooling capacity, where the size reflects the input change of the control system. Second is the heat production issue of the chemical reaction in the reactor with a variation in size that cannot be measured generally because of uncertain, strong, and nonlinear characteristics, raising higher requirements for the design of the controller.

The Hammerstein-Wiener model can describe the processes in the continuous stirred tank reactor (CSTR), pH neutralization reaction process, and other nonlinear chemical processes. Many research results were achieved in the nonlinear model predictive control (NMPC) [2–4]. Considering that the state of the actual system is usually unpredictable, the literature [5] presents a new two-step NMPC algorithm that constrains the output feedback of the Hammerstein system. The literature [6] designs the predictive controller of the batch reaction based on the Kalman filter; the DMC is used as the framework and is combined with the size of the interference evaluated by the Kalman filter. Its impact on the system output for the problems regarding uncertainty and strong nonlinear characteristics is reflected by the system during the heat production of the batch reaction itself. However, no reasonable method was found to resolve the complex relationship among the states of the practical nonlinear system, the constraints of the intermediate mass, and the impact of the interference signal on the output of the system.

Considering the Hammerstein-Wiener system of state, output, and intermediate mass constraints, this paper proposes the predictive control algorithm of the batch reaction based on the Kalman filter for the constrained Hammerstein-Wiener system by combining the literature [5, 6]. The paper designs the output feedback optimal control law of the linear subsystem with the state observer, and it evaluates the size of interference and its impact on the system output with the linear Kalman filter. By predicting the solving error of the nonlinear algebraic equation and by roll optimizing the error performance indicator of the finite horizon, the relationship between the outputs calculates the actual control variables. The simulation study of the CSTR proves that the control method in this paper overcomes the impact of the heat production in the reaction in terms of the accuracy of temperature tracking and improves the control effect to a certain extent.

#### 2. Description of Batch Reaction Process and the Hammerstein-Wiener System

##### 2.1. Description of Batch Reaction Process [1]

The batch reactor process is generally operated under a closed reactor condition so that the reactants are added to the reactor once to react in the given condition. The intensifying reaction will cause the heat to accumulate and the temperature to rise and further intensify. Obviously, this process, as well as the interlayer temperature, is unsteady. As the control process is complex, the batch reactor control system is generally required to design the cascade structure, in which the main circuit of the outer ring takes the interlayer temperature as the control input and the reaction temperature as the control output, using and , respectively. The secondary circuit takes the cooling water flow rate as the control input and the interlayer temperature as the control output, using and , respectively. The value of the secondary circuit is set up by the main circuit and the secondary circuit, is generally controlled with a simple PID. The external structure diagram of the reactor is shown in Figure 1.

##### 2.2. Description of the Hammerstein-Wiener System

The Hammerstein-Wiener model consists of two static, nonlinear modules and one dynamic linear module. Its model structure is shown in Figure 2.

Consider the following discrete-time multivariable Hammerstein-Wiener absolute system: where , and are the state variable, intermediate variable, and the input variable of the system, respectively; is the output variable in the linear part; is the constant coefficient matrix with the corresponding dimension; are the continuously differentiable nonlinear function vectors, indicating the static nonlinear link between the system input and the intermediate variable and between the linear part and the intermediate variable and meeting .

It is assumed that the system controllability, state observability, and the origin as the equilibrium point of the system, the parameters of the system states, intermediate variables, and the input constraints are defined as follows: where , and represent the compact convex set, including the origin as the interior point.

#### 3. The Constrained Predictive Controller Design Based on the Kalman Filter

##### 3.1. The Constrained NMPC Algorithm [5]

The Hammerstein-Wiener model can be divided into the Hammerstein model [7] and the Wiener model [8]. Different approaches are adopted for these two models to design the controller of the system and to deal with the impact of system output and interference signal on the system. This section mainly describes the designs of the output feedback optimal control law of the Hammerstein linear subsystem based on the state observer.

The considered output signal is assumed to have no direct control relationship with the input signal in formula (1). Therefore, in its corresponding linear subsystem , , ; if the observer status is , the optimal control law of the linear subsystem can be described as , where , is the positive definite weighting matrix of the intermediate mass, is a symmetric positive definite solution to the Riccati equation , and is the semipositive definite state weighting matrix.

The actual control variable can be calculated by solving the nonlinear equations , inevitably introducing the calculation error, and thereby affecting the closed-loop stability and optimality. Next, the actual control is determined through the prediction and the scroll optimization of the calculation error.

Considering the system in formula (1), the calculation error objective function is defined as follows: where the integer predicts the step length and the positive definite matrix is the input weighting matrix.

Thus, the NMPC issue of the constrained Hammerstein system is defined as follows:

The sequential quadratic programming is applied to calculate online the optimization problem in formula (4) and obtain the optimal control solution:

According to the principle of rolling optimization, the NMPC law based on the observation state is defined as follows:

In formula (6), the observation state can be obtained through the following all-dimensional observer: where is the gain matrix of the observer.

Control law (6) observes the system status with the state observer (7) at each sampling time. It then updates the initial conditions of the optimal control problem (4) on this basis and scrolls the optimization to calculate the predictive control variables of the current.

The observation error of the defined state is ; thus, the NMPC closed-loop Hammerstein system is where is the calculation error of the control variables; the controller gain is given by .

The NMPC closed-loop structure of the whole system is shown in Figure 3.

##### 3.2. The Predictive Controller Design Based on the Kalman Filter

This section mainly introduces the interference model [6] based on the Kalman filter and evaluates the size of the interference signal as well as its impact on the system output based on the Wiener model. Meanwhile, with the above constrained NMPC output, the relationship between their outputs solves the error through the predictive nonlinear algebraic equation, and roll optimizes the performance indicators in the finite horizon to calculate the actual control variables.

The input interference model is the key to achieving unbiased control and to improving the tracking performance. To improve the tracking performance and to realize the unbiased control, the system is characterized by introducing the interference model. The model is built in the input interference along with the influence of heat produced by the materials in the chemical reactor on the output. Davison and Smith [8] regarded this method as the standard in eliminating the steady-state error in the design of the quadratic regulator:

After the establishment of the input interference model, the input interference is evaluated. If the interference is the state vector, the following augmented-state space model is obtained: where and are the process noises when the mean value is zero.

According to the measured value of the current output, the system state and the interference of the current can be determined using the following formula: where and represent the corresponding Kalman gains of the state variable and the interference , respectively.

Then, the predictive value of the augmented state at the next moment is

For the sake of discussion, the following are set: , , , and . The Kalman gains of each cycle can then be calculated by , where is the covariance of the estimation error of the current; are the process incentive noise covariance matrix and the observation noise covariance matrix, respectively, and their values directly affect the performance of the Kalman filter. Greg Welch and Gary Bishop (2006) discussed it in detail in their writings. Interference has no direct relationship with input , but the observability of interference ensures passing the Kalman filter and then evaluates . When designing the controller, the evaluation value of is used to remove its influence on the system output. After evaluating the size of using the Kalman filter in each sampling period, feedforward compensation is made on this part, and its influence on the system output is included in the prediction model. The real-time control variable is where ; ; and deducts the influence of interference on the output from the expected value to form the new one.

The predictive control algorithm flowchart of the constrained Hammerstein-Wiener model with the feedforward compensation is shown in Figure 4.

#### 4. Simulation Analysis

The experimental platform is the CSTR tester that can simulate in real time the entire process of the batch reaction. The main parameters selected include the volume of the reactants: 200 L, the spin speed of the stirrer: 90 , the power of the stirrer: 1.5 Kw, the feeding of the reactant : 50 kg, and the specific heat: . The valve operation of the heating steam comprises manual and automatic modes. The initial stage of the reaction mainly involves the evoked response within the reactor. In other words, the temperature of the reactor is improved to 60°C with steam. The steam valve is manually operated in this process. The temperature in the reactor continues to rise because of the heat production in the reaction. The steam valve is switched to automatic mode to implement the control strategies. The process curve of the reaction temperature is set as follows: rising at 0.2°C/s, remaining constant at 120°C, and the overshoot of the temperature being less than 1.5°C. The result is compared with that of the traditional PID control, as shown in Figure 5.

Compared with the PID control, the PID control method cannot make the reactor temperature follow the process curve well because of the inherent nonlinearity and time lag characteristics of the object. However, the method in this paper achieves a good control effect.

In practical production, different products are produced through different processes, requiring corresponding control tasks to be changed, that is, selecting different process curves while keeping the algorithm unchanged. In the same simulation platform, the process curve is set to rise at 0.3°C/s and be constant at 120°C, with other parameters unchanged. The control result is shown in Figure 6(a), whereas the traditional PID algorithm control result is shown in Figure 6(b).

For batch production, the different natures of the reactants result in different exothermic characteristics in the entire reaction process, requiring the adopted control algorithm to have good robustness. In other words, the original control result remains even when the material property changes. Now, the reactant is used with the specific heat of . The reaction temperature process curve is set as follows: rising at 0.2°C/s and remaining constant at 120°C. The result is compared with that of the traditional PID control, as shown in Figure 7.

The traditional PID control algorithm is difficult to use in controlling the reaction temperature in a satisfactory range because of the specific heat change of the reactant. However, the method proposed in this paper shows good robustness, indicating practical significance to the complex and changeable strong, nonlinear, chemical process system.

#### 5. Conclusions

With the temperature tracking issue in the batch reaction process as the study object, as well as with the condition that the value of the heat production cannot be measured and its influence law on the output is unknown in the whole process of the reactant reaction, this paper proposes the predictive control algorithm of the batch reaction based on the Kalman filter for the constrained Hammerstein-Wiener system. Combined with state observer, the paper designs the output feedback optimal control law of the linear subsystem and characterizes the unmodeled dynamic system with the output interference model, which compensates for the impact of the heat production from the reaction on the system output and improves the accuracy of temperature tracking. The simulation experiment is conducted to verify the control results of this algorithm on different process curves and determine its robustness. The proposed algorithm, when compared with the traditional PID algorithm, fully demonstrates effectiveness and feasibility.

#### Acknowledgments

This work is partially supported by the National Natural Science Foundation of China (no. 61074020). The authors also wish to thank the associate editor and the anonymous reviewers for their valuable and constructive criticisms and suggestions.

#### References

- A. Jutan and A. Uppal, “Combined feedforward-feedback servo control scheme for an exothermic batch reactor,”
*Industrial & Engineering Chemistry Process Design and Development*, vol. 23, no. 3, pp. 597–602, 1984. View at Google Scholar · View at Scopus - H. Song, W. Gui, and C. Yang, “Identification of Hammerstein-Wiener model with least squares support vector machine,” in
*Proceedings of the 26th Chinese Control Conference (CCC '07)*, pp. 260–264, July 2007. View at Publisher · View at Google Scholar · View at Scopus - M. Hong and S. Cheng, “Nonlinear model predictive control based on LS-SVM Hammerstein-wiener model,”
*Journal of Computational Information Systems*, vol. 8, no. 4, pp. 1373–1381, 2012. View at Google Scholar · View at Scopus - J. C. Gómez, A. Jutan, and E. Baeyens, “Wiener model identification and predictive control of a pH neutralisation process,”
*IEE Proceedings*, vol. 151, no. 3, pp. 329–338, 2004. View at Publisher · View at Google Scholar · View at Scopus - H. E. De-feng, Y. U. Li, and Z. O. U. Tao, “Output feedback nonlinear predictive control for constrained Hammerstein systems,”
*Control Engineering of China*, vol. 16, no. 4, pp. 416–418, 2009. View at Google Scholar - X. F. Zhu and D. E. Seborg, “Nonlinear predictive control based on Hammerstein models,”
*Control Theory & Applications*, vol. 11, no. 5, pp. 564–575, 1994. View at Google Scholar · View at MathSciNet - A. L. Cervantes, O. E. Agamennoni, and J. L. Figueroa, “A nonlinear model predictive control system based on Wiener piecewise linear models,”
*Journal of Process Control*, vol. 13, no. 7, pp. 655–666, 2003. View at Publisher · View at Google Scholar · View at Scopus - E. J. Davison and H. W. Smith, “Pole assignment in linear time-invariant multivariable systems with constant disturbances,”
*Automatica*, vol. 7, pp. 489–498, 1971. View at Google Scholar · View at MathSciNet