A New Control Performance Evaluation Based on LQG Benchmark for the Heating Furnace Temperature Control System

: Temperature control systems are a series of processes with large time ‐ delay and non ‐ linear characteristics. Research shows that using fractional ‐ order modeling and corresponding control strategies can better control these processes. At the same time, the existing studies for control performance assessment are almost committed to the integer order control systems, and the methods used in few literatures on performance assessment of fractional order systems are also one ‐ sided. This paper applies the linear quadratic Gaussian (LQG) evaluation benchmark to the performance evaluation of fractional ‐ order control systems for the first time, starting with the LQG evaluation benchmark considering the input and output performance. The LQG benchmark can be obtained by the analytical algorithm, which simplifies the complexity of LQG solution. Finally, taking the application of the fractional predictive function control (FO ‐ PFC) controller in the experiment of industrial heating furnace temperature control as an example, the effectiveness of the LQG benchmark is verified.


Introduction
In the industry production process, the high-performance controller is the essential element of advanced control systems.Under the complex industrial environment, the performance of a loop controller may decline after a long time of operation [1,2].Therefore, it is great practical significance to make use of the control performance assessment technology to evaluate and monitor loop controllers online in real time, in order to ensure the high-performance operation of the loop controllers and the safe and efficient operation of the industrial control system [3,4].
The research about control performance assessment technologies has a history of more than 40 years.In 1978, Devries et al. [5] first proposed the concept of performance evaluation and analyzed the effectiveness of quantitative control of paper machine by using multiple time series technology, which greatly stimulates the research in the field of subsequent evaluation of controller performance [6][7][8][9][10][11][12][13][14][15][16][17].When a two-dimensional Roesser model was known, Wei et al. [18] proposed a method for evaluating the control effect of iterative learning controller (ILC), which extended the LQG benchmark from onedimensional to two-dimensional firstly.When parameters of the batch process were unknown, Wang et al. [19] developed a model-free control performance assessment method based on date-driven control in 2015.Different from the finite-horizon LQG benchmark in reference [18,19], Wang et al. [20] adopted an infinite-horizon quadratic performance cost function to design LQG benchmark and extended the results to a multiple input multiple output (MIMO) system.For disposing the inconsistency of the minimum variance criterion when evaluating non-Gaussian perturbation systems, Zhang et al. [21] proposed an improved Renyi entropy benchmark for the performance evaluation of common cascade control systems in 2019.By the performance index combined entropy with the output mean, the running state of the control system was evaluated.
Due to the widespread use of proportional integral derivative (PID) controllers and model predictive control, the control performance assessment of them is a hot research spot.For achieving the minimum variance to evaluate control loop performance, Ko et al. [22] and Sendjaja et al. [23] adopted a minimum variance control (MVC) benchmark for evaluation which can get excellent consequences by calculating the proportion integration differentiation regulator.Horton et al. [24] showed the optimum performance of the PI regulator, which was as the widely-used controller.Fang et al. [25] adopted LQG benchmark to assess running state of the IMC-PID controller under the circumstance of model matching and mismatch in experiment of industrial furnace temperature control.By considering the tradeoff between LQG cost and average sampling rate, an interpretation of the optimal sampled-data controller and event-based sampling policy were obtained in terms of PID control [26].This method helped us understand what implementation characteristic was important for obtaining excellent event-based PID performance.However, in industrial production process, integer order model is insufficient to precisely describe information of a system.Additional fractional integration and differentiation can make the process model more flexible and accurate, and reflect the dynamic characteristics of the actual process [27][28][29].In order to expand the traditional concept of integral order controller, a fractional order controller is emerged.According to previous studies, the control effect of fractional order controller is more excellent than integral order controller [30][31][32][33].Consequently, the methods and techniques of performance evaluation on fractional order controller are also worthy of studying.Škarda et al. [34,35] proposed a new performance assessment indicator by using the ideal shape based on the sensitivity function.By calculating the fractional order process, the optimal performance control for all PID control models was archived.The evaluation method of minimum variances based on the classical theory that the normalization process dead time great impact on maximum performance, so only a fixed structure for the evaluation process controller.However, there are some limitations.For solving the tradeoff between servo controller performance and the adjustment controller performance, a combined performance evaluation index is developed by Meneses et al. [36].The method carried out the theoretical research that may be implemented on the basic of the performance of the process delay time normalized and constrained by the robustness but also for evaluating the fractional order PID controller.However, considering that there are drawbacks of insufficient application scope on the existing performance assessment methods for fractional order controllers, these indicators structured performance appraisal are one-sided and the angle output is from the process point of view.
Although LQG benchmark is a more practical industrial process evaluation criterion, which considers both of the variance of the system input and the variance of the system output, its application is limited to the order of the process model.Therefore, it is significance to extend the LQG benchmark to assess the running status of fractional order system with fractional order controllers.This paper uses an analytical algorithm to solve the LQG benchmark.Taking the heating furnace temperature controller as an example, which is called FO-PFC controller, the validity of the LQG benchmark is verified.In the meantime, the effectiveness of FO-PFC controller is studied whether the model matches or not.
The main contribution of this paper is applying the LQG benchmark to the performance assessment study of a fractional order controller for the first time.Section 2 covers solution algorithm of LQG benchmark.The block diagram of heating furnace control process is presented in Section 3. In Section 4, LQG benchmark is used to assess the FO-PFC controller of the heating furnace temperature control system.

LQG Benchmark
The LQG benchmark is based on MVC benchmark, which extends and takes the output and input variances of the system into consideration.Consequently, a more nuanced state of the control system is excavated by LQG benchmark.This provides the optimal performance limits of the controlled system in respect of input and output variance under any linear controller.The tradeoff curve is obtained by solving an LQG problem where the cost function of LQG is defined as: In Equation ( 1), Q is the output weighting matrix and R is the input weighting matrix.The LQG problem can be solved by changing the input weighting  , and a class of LQG optimal discrete points can be obtained.The horizontal axis represents the best input variance and the vertical axis represents the best output variance, then the LQG compromise tradeoff curve [25,37] is obtained by regressing these discrete points as described in Figure 1, which indicates the relationship between the optimal input variance and the optimal output variance.Through the relationship, the minimum performance value and input-output variance of the linear controller can be obtained.However, because the disturbance is always existed in industry, the top of the curve in Figure 1 shows what may happen in the actual operation.Var Y are the optimal input and output variances.
Then the corresponding performance index can be defined as:

LQG Solution via State Space Model
Consider a state space model: where A B C 、 、 are coefficient matrix.k  is Gaussian white noise and H is the Kalman filter gain matrix.k X and k Y are the state vector and output vector at time instant k .k U is the input vector at time instant k .The process equation and Kalman estimator are integrated as follows: where ˆk X is the estimated state at time instant k .

=( ) ( )
which can be obtained through the MATLAB optimization toolbox.S can be gained by calculating the Riccati equation described as is weighting matrix state in the quadratic index.N is the operation weight of the state matrix, which is set to zero.Then we can solve for ( ) Var Y as follows: Due to the limited memory of main control cards in distributed control systems (DCS), the amount of computation must be allocated reasonably.Therefore, the rapidness and simplicity of the LQG solution must be considered to save the memory for the use of controllers and security interlock programs in the process units.By combining a subspace method, a more accurate modelfree LQG solution can be achieved (If anyone wants to know more about subspace identification, please refer to Appendix B. The way of getting data repository can be seen in Appendix A).But the model-free LQG solution process is too complicated, which will take up more memory of main control cards.From the perspective of the application, the proposed LQG benchmark based on a known model subspace is simpler and the calculation is smaller, which can be directly embedded into an industrial system.

FO-PFC Controller Principle
As one of the advanced control methods, PFC can get more excellent control performance in some high-precision processes.Traditional PFC, which is based on the integer order model and the corresponding control method, has certain limitations in controlling fractional order systems.Studies have shown that fractional order modality can be added to the objective function of conventional PFC for optimization.Better process tuning performance can be achieved than conventional PFC because the process model is described more accurately.In this paper, the process model of industrial heating furnace temperature is built through the fractional order model, and the control method of the system is FO-PFC algorithm.The block diagram of control process is displayed in Figure 2. According to literature [38], the optimal control law of FO-PFC is shown as follows: where s T is sample time,  is the cost function of integral coefficients.
, , Y U E are corresponding matrix constructed from input, output and error respectively, , , B C D are corresponding coefficients of matrix (see more details for reference [38]).P is the predictive horizon.

Case Study
In this part, the process of the heating furnace temperature is selected as an example and a fractional order model is selected as the temperature process model.The fractional order model is developed through corresponding measurement data and the controller of the temperature is FO-PFC controller.Finally, the performance of control system is tested by LQG benchmark through MATLAB simulation.

System Description
The temperature controlled object is sxf-4-10 industrial electric heating furnace in this paper.The principle diagram is given in Figure 3.The voltage during normal operation is 220 V and the power during normal operation is 4 KW.The control flow of the furnace heating system is as follows: heating furnace temperature is sampled by thermocouple and then transferred into the computer through signal amplifier and signal acquisition card.The FO-PFC controller runs on the computer and the output of the FO-PFC algorithm exports to the signal output card to control on and off state of the solid state relay.When the solid state relay is on state, the state of the heating furnace is heating.Otherwise, the heating loop is in the off-state.

Process Model Establishment
An input of 20% duty-ratio is given to the heating furnace and the heating furnace begins to heat.During the heating period, real-time data of furnace temperature is collected every two seconds until the output trend of the heating furnace temperature is stable with no obvious change.The actual temperature responses of this process are shown in Figure 4.For general industrial process control systems, the first order plus dead time (FOPDT) model is generally chosen to describe the relationship between inputs and outputs, which is beneficial for the convenience of subsequent controller design.Here, in order to describe the actual process model more accurately, the fractional order model is chosen, which the new luus-jaakola (NLJ) method will be used to identify.
The obtained model is: Since the FO-FPC controller to be evaluated is designed as a nominal model, the fractional order model is replaced by the high-order integer order model, which is obtained by the Outsaloup approximation method [39].According to the Outsaloup approximation method, some extra parameters must be chosen in the appropriate value.The order of the derivative operator is 4, and the higher and lower frequency approximation interval are where: The transfer function model of Equation ( 8) is further discretized.The sampling interval is 20 s and the discrete model of Equation ( 8) can be obtained as shown below: where: The practical temperature data, the responses curve of the fractional order model and highorder integer order model are focused on Figure 4.The goodness of fit of the two models is superior even if the allowable error still exists in the Oustaloup approximation method.In this paper, the perturbation z function is assumed as follows: To sum up, the process model of the furnace is consists of two parts, which can be expressed as:

Computation of LQG Benchmark
In order to test the stability of the temperature control system, a Gaussian white noise is added to the system, which has average of zero and the variance is 0.1.
The process transfer function is converted into the state space type like Equation ( 3) for convenient calculation.Through the above method shown from Equation (3) to Equation ( 5), the LQG benchmark can be computed.At the same time, since the interval  must be sufficient to delegate the overall trend of the LQG curve, the weighted value  changes within [0.000000001, 160000].Furthermore, the corresponding variance of the optimal input-output is calculated.After that, the LQG performance curve is obtained which is shown in Figure 5.

Simulation Results
From the reference [38], the predictive horizon, coefficient of smoothness and the cost function of integral coefficients of the corresponding PFC controller are chosen as follows.However, some uncertainties and the immeasurable disturbances are inevitable in practical production, which lead to some deviations in establishing the model.To verify the stability of the controller, the cases of model match and mismatch need to be considered.Model 1 is the match model, and Model 2 and Model 3 are mismatch models.Three groups of the model are as follows: Because the running condition of the FO-PFC controller under the steady state is our main evaluation object, 600 °C and 605 °C are chosen as the temperature set-point.Figure 6 shows the performance evaluation consequences of actual temperature and duty cycle changes of the solidstate relay, which are obtained through multiple experiments.On the basis of experimental results, the corresponding parameter and calculation results can be obtained in Table 1 under  From the above results, it can be seen that the heating performance of the furnace can be greatly improved.By means of the evaluation and simulation experiments, the parameters of the FO-PFC controller, the integral and differential coefficients can be adjusted by constantly trying to use the wrong approach.Through a lot of experiments, an optimal set of controller parameters can be obtained, which can let control system get the best control performance.The group of optimal control parameters is:    Through the actual input variances and output variances comparison of improved parameter and current parameter, which are shown in Figure 8a,b, no matter in the condition of model matching or model mismatching, the improved parameter can obtain excellent control performance.Regardless of the model matching or the model mismatch, reducing the input and output of variances can be done by adjusting the FO-PFC parameters to a certain extent, which can be seen from Figure 7a to Figure 8b.After the expected parameter adjustment, the controller improves a lot in terms of control performance.Compared with the current parameter, the performance point of the improved parameter is closer to LQG performance tradeoff curve.In addition, whether the model matches or not, the gap of the input and output variances is obvious by contrasting improved parameter with current parameter, which are displayed in Figure 8a,b.When the temperature set-value is 600 °C, the input and output variances decrease by 32.09% and 7.56% respectively in Model 1; In Model 2, the input and output variances all decrease that the input variances decrease by 28.24% and the output variances decrease by 10.87%;In Model 3, the input and output variances all decrease that the input variances decrease by 19.69% and the output variances decrease by 9.07%.When the temperature set-value is 605 °C, the input and output variances all decrease that the input variances decrease by 16.04% and the output variances decrease by 10.34% in Model 1; In Model 2, the input and output variances all decrease that the input variances decrease by 18.94% and the output variance decrease 9.02%; In Model 3, the input and output variances all decrease so that the input variance decrease by 15.13% and the output variances decrease by 10.27%.Due to the smaller output variances, the stability of the controller in improved parameter is more superior to that of the current parameter.By redesigning the controller parameters after the evaluating the running status of the controller, the effect of the controller is better.The LQG benchmark provides a graphic way for manually adjusting controller parameters.

Conclusions
The proposed LQG benchmark is applied to the performance assessment of the FO-PFC controller in the heating furnace temperature control system.The analytic method is used to calculate the LQG benchmark, which simplifies the complexity of LQG solution.Through the solution of linear quadratic Gaussian problem, the optimum values of input variances and output variances under different weight coefficients can be obtained.The control effect and performance point under different parameters are obtained through simulation experiments.Finally, the LQG benchmark can be used to evaluate the control performance of the FO-PFC controller whether the model matches or not.When control performance still has great room for improvement, we can use LQG benchmark to seek for the more optimal control effect of the controller than before by adjusting the parameters.This benchmark enables us to assess the control effect of the controller in the control experiment of the heating furnace, and adjust its parameter in time if the expected value is not achieved.

Figure 1 .
Figure 1.Linear quadratic Gaussian (LQG) performance tradeoff curve.According to the LQG performance tradeoff curve, the current input variance is ( ) k Var U and

Figure 2 .
Figure 2. The block diagram of control process.
function model of the high-order integer order is obtained as follows:

Figure 4 .
Figure 4.The practical temperature data and the response curve of the two models.

Figure 6 .
Figure 6.Limit curve based on LQG performance and actual variance of the system.

.
the condition of model match and model mismatch.The performance of the FO-PFC controller has been evaluated against the performance indices.When the temperature set-point is set 600 °C: under the circumstance of match model shown as Model 1, the performance indicators are 34Under the condition of mismatch which are shown as Model 2 and Model 3, the performance indicators are 42temperature of the heating furnace can be set arbitrarily according to the demand.When the temperature is set to 605 °C: the corresponding performance indices of After adjusting and recalculating the relevant data and performance indicators, Figure7a-cfurther show the simulation results and the selected model parameters are three different parameters, which make the experiment more convincing.

Figure 7 .
Figure 7. (a) A comparison of the improved parameter and the current parameter (Model 1).(b) A comparison of the improved parameter and the current parameter (Model 2).(c) A comparison of the improved parameter and the current parameter (Model 3).

Figure 8 .
Figure 8.(a) The input variance comparison of the improved parameter and the current parameter.(b) The output variance comparison of improved parameter and current parameter.

Table 1 .
Statistical Results of Steady State Performance.

Model k Var(u ) k Var(y ) opt lqg Var(u ) u η (%) opt lqg Var(y ) y η (%)
When the temperature set-point is 600 °C, under the circumstance of model matching which is shown as Model 1, the performance indicators are improved to 66.98%

Table 2 .
Statistical results of steady state performance.