Research Article  Open Access
Angang Chen, Zhengyun Ren, Zhiping Fan, Xue Feng, "Combined Integrating Control Based on Dynamic Optimization Estimation", Mathematical Problems in Engineering, vol. 2020, Article ID 5917606, 14 pages, 2020. https://doi.org/10.1155/2020/5917606
Combined Integrating Control Based on Dynamic Optimization Estimation
Abstract
In this paper, a new type of controller based on the combined integrating systems is proposed, which is called a combined integrating controller. The controller’s current output is determined only by the current input and the average value of the controller’s previous output. The innovation is to introduce a timedelay positive feedback in the controller. The memory characteristics of the timedelay term keep accumulating error information and repeatedly learning so that the system can stably track and restrain any input signal without error. Simultaneously, the combined integrating controller is applied to conventional process systems based on dynamic optimization estimation in the case study to show absolute superiority over the nonpredictive control method (such as the classical PID control method).
1. Introduction
The ProportionalIntegralDerivative (PID) control algorithms and their variants are widely used in the process industry because they can meet most control requirements for various processes. Most PID tuning rules are based on conventional process models in the industry, such as the firstorder process with timedelay (FOPDT), secondorder process with timedelay (SOPDT), integrating processes, and highorder processes. The widespread use of FOPDT models, especially in PID controllers, is due to their simplicity and ability to capture the essential dynamics of multiple industrial processes [1–3]. On the other hand, many industrial processes can be successfully approximated by a secondorder timedelay transfer function. Especially in the chemical process industry, most control loops are PIDtype secondorder modeling plants or temperature control, as shown in [4, 5]. Compared to stable process models, the number of PID adjustment rules reported for the integrating process is less in the literature. The integrating process is challenging to control because load disturbances easily disrupt the balance between the input and output. In [6], controllers of the modified Smith predictor are designed for pure integrating, integrating plus first order and double integrating processes with large timedelay. However, if there is a larger timedelay in the process, the control performance obtained using the PID controller will be limited [7]. It may cause the dynamic performance of the control system to decline or the stability of the entire system to deteriorate and even cause system oscillation or divergence. Numerous methods of PID controller synthesis and tuning have been developed over the past few decades, such as ZieglerNicholas based tuning methods [8], internal model control (IMC) design methods [9, 10], frequencydomain method [11–13], and decomposition method based on multiscale control (MSC) scheme [14, 15].
With the indepth research on the process industry model, many complex processes may not be able to adequately describe with conventional models, such as fractional integral and differential express the dynamics of the actual process more flexible and accurate [16]. The distillation column bottom process, which Morari had referred to in [17], has a transfer function , which is a nonminimum phase timedelay process. Morari designed the controller of this kind of process based on the IMC and verified the control performance.
At the same time, through the analysis of the mechanism of the productive process, we get a new type of system transfer function that contains the term and call it a combined integrating system. This kind of system has been widely used in the fields of steel, chemical industry, and tobacco through our engineering practice. The most combined integrating systems are regarded as a loworder timedelay (LOTD) using the conventional control strategies, such as the PID control algorithm, Smith predictor algorithm, or IMC algorithm. To some extent, these control algorithms satisfy the needs of industrial production, which is why these processes are not taken seriously.
Nevertheless, with the rapid development of modern industry, product quality, energy consumption, and environmental protection put higher requirements on precision control. The traditional control has gradually been not suitable for combined integrating systems. Therefore, it needs to design new control strategies and methods to satisfy higher requirements. In this paper, a new type of controller based on the combined integrating systems is proposed, which is called a combined integrating controller. The purpose is to provide a deadtime compensating controller that can handle multitimedelay systems and manually adjust control parameters in the same way as PID controllers.
In the classical control period, the timedelay term is often regarded as a fixed part of the system, to be stabilized rather than as a control parameter. In [18], the timedelay term is introduced into the controller, and the timedelay coefficient is taken as the control parameter. In this paper, the timedelay and the combined integrating term are introduced into the controller and as control parameters. The effective combination of timedelay term and combined integrating term with “slow computing” characteristics is used to construct a controller with a simple structure and reduced computational complexity. In the following simulation, we will give the theoretical feasibility. The advantages of introducing the timedelay term and combined integrating term as control parameters in the control system are as follows: (1) the controller is simple, and the parameters depend on the timedelay and combined integrating term of the model. (2) The controller’s current output is only determined by the current input and the average value of the controller’s previous output. (3) The asynchronous operation of the control system is realized, which is easily converted into a digital controller, similar to a hybrid system.
This paper is organized as follows. In Section 2, dynamic modeling and estimation are introduced. In Section 3, based on the internal model principle and predictive control, we propose a combined integrating controller and give a detailed derivation process. In Section 4, the combined integrating controller mechanism analysis is introduced. In Section 5, the combined integrating controller robustness analysis is presented. In Section 6, we will compare the different control algorithms for the typical industrial processes based on dynamic optimization estimation. Finally, we give conclusions and further study in Section 7.
2. Dynamic Modeling and Estimation
The critical aspect of any controller is getting a sufficiently accurate model. The model form can include adjustable parameters that cannot be directly measured but can be adjusted to match steadystate and dynamic data. Models can be linear or nonlinear, empirical, or based on primary forms produced by materials and energy balances, reaction kinetics mechanisms, or other predefined model structures. The model structure can also include constraints such as fixed gain ratio, constraints on components, or other boundaries reflecting physical reality. Equation (1) is a statement of the general model form, which can include differential, algebraic, continuous, and integer variables:
The target value of equation (1) is determined by the initial state , a set of parameters , a trajectory of disturbance values , and a sequence of control moves . The equation includes differential , algebraic , and inequality constraints that are included to model physical phenomena. In addition to the model form, the objective function is vital to ensure desirable results. A common objective form is a leastsquares form: . A suitable initial guess for the state trajectories is often required for solver convergence.
The direct application of DAE to nonlinear optimization faces many challenges. In this paper, it is assumed that the dynamic model can be described by a set of ordinary differentialalgebraic equations (DAEs). The discretization of continuoustime representation allows largescale nonlinear programming (NLP) solvers to find solutions at specified intervals over time. There are many names and related techniques to find the mathematical relationship between the derivative and nonderivative. These techniques include orthogonal collocation on finite elements, direct transcription, Gauss pseudospectral method, Gaussian quadrature, Lobatto quadrature, Radau collocation, Legendre polynomials, Chebyshev polynomials, Jacobi polynomials, Laguerre polynomials, and so on [19–23]. Here, we use the orthogonal collocation method on the finite element with Lobatto quadrature to solve the DAEs (the differential term in equation (2)):
Euler’s method is based on twonode per horizon step linear discretization by
The orthogonal collocation method on the finite element with Lobatto quadrature uses nonlinear polynomial fitting. In engineering applications, each horizontal step can be linear within 3–6 nodes. Here, we give four nodes orthogonally by Lobatto quadrature: , , , :
In this derivation, a thirdorder polynomial approximates the solution at the four points on the horizon. The increase in the number of collocation points will increase the corresponding polynomial order. For the initial value problem, when the initial time is arbitrarily defined as zero, the coefficient is equal to . In order to determine the coefficients , , and , equation (4) is substituted into equation (5), and the results are obtained. Note that the coefficient in equation (4) is offset by to the right of equation (5). Finally, the rearrangement and solution of give the solutions of the three nodes in the future which are given as follows:
A simultaneous solution demonstrates the application of orthogonal collocation. The firstorder system is solved at 3–5 points from to using orthogonal collocation on finite elements with Lobatto quadrature. In this case, , , and the initial condition is specified at . For this problem, the time points for and are selected by Lobatto quadrature. The value of is specified at due to the initial condition. As shown in Figure 1, as the number of collocation points increases, the accuracy increases and the relative error decreases.
(a)
(b)
3. Combined Integrating Systems and Combined Integrating Controller
Definition 1. The system transfer function consists of a combined integral term, integrating plus timedelay process, and a general process system transfer function. It can be expressed as the following general equation:We define these openloop stable systems as combined integrating systems. The combined integrating term is an essential part of equation (6); is a general process system transfer function without integrating. In recent years, it is found that four typical controlled processes are widely used in the industry. Their transfer functions areIn the industrial process, the firstorder Padé approximation is often used for the timedelay. So we use this approximate method for the combined integrating systems, and the combined integrating systems are transformed into the conventional industrial systems (LODT). However, its characteristics also have changed. The combined integrating systems are nonminimum phase systems, and the approximation results are transformed into minimum phase systems. The corresponding controllers are designed by the traditional control method; though it can get the control performance satisfying the production requirements, it is difficult to improve the control performance further.
Due to the specific characteristics of the combined integrating systems, using the traditional PID algorithm, IMC, prediction PI (PPI) algorithm, and Smith predictive control algorithm cannot reach good control effect, rapidity, and stability. It needs to design a new control algorithm. The controller design of the combined integrating systems (7a) and (7b) is analyzed below.
4. Combined Integrating System (7a)
By choosing the appropriate desired closedloop response (see equation (8)) and using either a Padé approximation or a powerseries approximation for the timedelay, PI/PID controllers can be derived for process models commonly used in industrial applications [24, 25]:
According to the internal model controller's design principle, a twodegreeoffreedom controller is obtained, which can achieve the setpoint tracking performance, and the transfer function of the controller is
The relationship between the input and the output of the controller is
Equation (10) interprets that the current output of the new controller is determined only by the current input and previous controller output (see Figure 2).
The controller itself is a timedelay positive feedback controller, and its structure and realization are convenient and straightforward. The controller has good antidisturbance ability, response speed, and robustness. However, after the industrial application of the controller, it is found that the output oscillation of the controller is relatively large, and the wear of the actuator is relatively large, which severely limits its popularization in the industry.
To overcome the shortcomings of the controller output oscillation, it needs to smooth the positive feedback control part. Considering the combined integrating term average filter nature, the timedelay in Figure 3 is replaced and the is reduced to , so we get the following improvement of the controller structure.
In Figure 3, is the Laplace transform of the average value of output in the prior period ; the relationship between input and output of the controller in Figure 3 is
The transfer function of the controller is
The closedloop transfer function of the whole system is
The closedloop transfer function is consistent with the openloop transfer function. Assuming discrete time is , , , and are integers. Considering the convenience of implementation, transform equation (11) into an incremental form:
As can be seen from equation (14), the increment of controller output relates only to the current and previous inputs as well as the controller output in the prior period .
5. Combined Integrating System (7b)
The system (7b) has a pure timedelay more than the first system, following the front controller design ideas. The current controlled variable is only related to the controller output in the prior period . The controller’s positive feedback term is replaced by , and the following controller structure is obtained, as shown in Figure 4.
The relationship between the input and the output of the controller is
The input and output of equation (15) are discretized and transformed into an incremental form (assuming discretetime is , , , , , and are integers):
According to the above derivation, the definition of the combined integrating controller is a class controller with the combined integrating term , the controller output by the current input and the average value of previous controller output.
Table 1 lists the combined integrating controllers of various types of combined integrating systems.
6. Combined Integrating Controller Structure Analysis
From equation (15) and Table 1, it can be seen that combined integrating controller structure contains a combined integrating term and a timedelay term. The timedelay always is regarded as an unfavorable factor. More and more researchers show that the impact on timedelay on system performance is not always negative. By judiciously introducing timedelay into the system, the system performance and stability can be improved [26], which translates the unfavorable timedelay into a favorable one.
The current output of the combined integrating controller consists of the current input and the average controller output over the prior period. Due to the existence of , the combined integrating controller has the predictive function and eliminates the system steadystate residual error function. According to the definition of equation (15), we give the controller’s specific parameters, that is, , and simulate the step response of the controller (at t = 10 s plus step input). The controller response curve can be obtained, as shown in Figure 5.
From Figure 5, the output of the combined integrating controller consists of three parts: proportional action, predictive action (antidelay), and integral action. Compared with the output of the traditional PID controller, there is an essential difference. Firstly, the prediction function is different. If the controller’s input remains the same, the predictive effect will keep the controller unchanged in the pure timedelay. Secondly, the integral action is different; the strength of the integral is different at different periods. The most crucial point is that there is a different mechanism for achieving a steadystate residual error: the PID is achieved by direct integration of the inputs, and the combined integrating controller is obtained by averaging the output over a period. When the input remains constant, the average output value of the controller increases continuously due to the positive feedback structure of the combined integrating controller, to eliminate the steadystate residual.
7. Combined Integrating Controller Robustness Analysis
Robustness is the ability of the closedloop system to be insensitive to component changes. It is one of the most useful properties of feedback. The presence of large timedelay and parameter uncertainties in the process model poses a severe challenge to the design of robust feedback controls for chemical processes [27]. It is well known that a welldesigned control system should meet the following requirements besides nominal stability: disturbance attenuation, setpoint tracking, robust stability, and robust performance [28–30].
The tradeoff between stability and steadystate tracking accuracy is an important aspect of the performance of the combined integrating system. The stability analysis of the control system with a combined integrating controller is essentially stability analysis of a class of systems with multitimedelay. In recent years, many researchers have published relevant results and concluded that the stability space of multitimedelay processes is discontinuous [29]. It is methodologically feasible to apply these research results to the stability research of the subject [31]. Such systems are timedelay systems with , which bring great difficulties to the stability analysis of control systems. Zheng [32, 33] uses the stability analysis method based on Kronecker’s theory and Routh array to analyze the multitimedelay system with coefficients and delay uncertainties.
The research of robust stability does not need to obtain the mathematical expression of stability. However, it needs to get the relationship between the stability and the controller parameters and provides the theoretical basis for the design and parameter tuning of the combined integrating controller. Therefore, to control more applications in the actual industry, we use the most practical industrial field debugging method to analyze the relative stability.
From Figure 6, the closedloop system characteristic equation is as follows:where , , and are all positive.
Theorem 1. Assume and is a monotone decreasing function at .
However,
Then, if ,. Assume , and then only can satisfy . Therefore, the relative stability is obtained by setting the range of controller parameter ; then, we get the more accurate range by the relevant approximation:
According to Routh array stability,
The combined integrating controller is used to stabilize the feedback system, which can achieve relative stability and a certain degree of parameter tuning. However, the robust stability of the combined integrating controller will be illustrated by the concrete simulation examples.
Assume that the parameters of the combined integrating system (7B) are as follows: , , and , which is a kind of large timedelay process, using equation (15) to control. The step response curve and antijamming characteristics are shown in Figure 7 with a step increase at t = 10 s and a disturbance amplitude increase 1 at t = 45 s).
When the model of controlled process is mismatched, the parameters are , , and . While the controller remains the same, the system can still be well controlled as shown in Figure 8.
For the model structure mismatch of controlled process, it changes to the firstorder largetimedelay process (with a step increase at t = 10 s and a disturbance amplitude increase 1 at t = 100 s). Similarly, the controller remains unchanged, the response speed of the system slows down to a certain extent, and the smoothness is lacking, but the system still has excellent control performance: there are no overshoot and no steadystate error and antijamming performance, and the controller output is relatively smooth, as shown in Figure 9.
8. Case Study
Above, we just verify that the proposed controller has excellent control characteristics for the combined integrating systems. Next, we compare it with PID and Predictive PI control algorithms for common industrial processes to verify the effectiveness of the control algorithm. Various performance indicators were selected for testing and comparison (e.g., response time, settling time, overshoot, and integral error (ITAE)) [28–30, 34].
In the modern sense, cigarette products’ production process roughly goes through seven major production processes: tobacco leaf precuring, tobacco leaf recuring, tobacco leaf fermentation, cigarette formula, cigarette shreds, cigarette rolls, and cigarette packaging. The tobacco leaf precuring cannot be directly supplied to the cigarette factory but must be recured, that is, the second moisture adjustment of the tobacco leaf based on the tobacco leaf precuring, making it a real raw material for cigarette production. Before the raw tobacco is recured, the tobacco slices and the stems are separated by the leafbeating device. Tobacco leaf and stems were rebaked separately, packaged, and stored separately. Tobacco leaf precuring has uneven moisture content during grading and binding, and the climate and storage conditions are different in different regions.
The role of tobacco leaf recuring: first, adjust the moisture to prevent mildew; second, remove impurities and purify aroma; third, insecticide and sterilization and beneficial storage; fourth, maintain color, and facilitate production. After reheating and drying, the physical and chemical characteristics of the recured tobacco leaf are further optimized. The quality of the tobacco leaf is improved, and the hygroscopicity is weakened, which is conducive to the natural alcoholization of the tobacco leaf during storage. To meet the requirements of threshing and reduce the shattering during threshing, it is necessary to increase the temperature and humidity of the tobacco leaf before threshing to improve the toughness and processing resistance of the tobacco leaf.
WF3 series hotair leaf moisturizing machine needs to carry out secondary tobacco leaf moisture processing on the tobacco leaf after the moisture regaining machine so that the tobacco leaf’s moisture increases from about 12% to 16% after the vacuum moisture regains, and the temperature reaches about . The main body of the hot air leaf moisturizing machine is a slanting slewing body. Driven by the power device during operation, the cylinder continuously rotates, and the tobacco leaves first enter the drum from the upstream vibrating conveyor. On the inner surface of the drum, there are spirally distributed stirring rods. The rotation of the drum, the action of the stirring rod, and the tobacco leaf make a spiral movement with the cylinder. The drum rotates round and round, and the tobacco leaves are gradually sent out of the drum outlet.
During the flow of tobacco leaf in the cylinder, due to the relatively sealed state of the cylinder, the air jets from the three nozzles of steam, watersteam mixer, and hot air at the inlet form a hightemperature and highhumidity mixed airflow that quickly fills the relatively tight cylinder space, thus forming a high temperature and high humidity environment in the cylinder. The tobacco leaves are in full contact with the hightemperature and highhumidity steam and heated by the ambient humid and hot gas and absorb the moisture in the humid airflow. The heat and mass transfer between the tobacco leaf and the humid and hot gas causes the tobacco leaf to heat up and absorb moisture quickly. Under the joint suction action of the moisture discharge fan and the circulation fan at the outlet end of the drum, the hot and humid airflow in the drum flows from the inlet end to the outlet end of the drum. Figure 10 shows valves control the inlet flow of mixed water vapor and hot air. and are moisture sensor and transmitter from imported and exported tobacco leaf, is the flow sensor and transmitter of hot air, and is the flow sensor and transmitter of mixed water vapor.
Most of the tuning techniques are proposed for FOPDT. This paper presents the technique for obtaining the FOPDT model using dynamic optimization estimation. Most process engineers know how to adjust the different parameters of the controller to obtain desired closedloop behaviors. They also know how the three parameters , , and influence the control. The empirical model building aims to find a relationship between the response (output) variable and the input .
The controller requires an estimate of the static process gain, the dominating time constant, and the dead time. As seen above, these parameters are easily obtained from a step response experiment. Graphical method to obtain the static process gain, the apparent dead time, and the dominating time constant from a step response is given in Figure 11. However, this graphical method can only get an imprecise result. For the highprecision control system, it is impossible to use this identification method. Therefore, we choose to use the above dynamic optimization estimation method, at the same time, based on the results of the graphical method, as the initial value of dynamic optimization estimation, to increase the convergence and solution speed of optimization.
(a)
(b)
(c)
(d)
We get the system process model and input pressure disturbance model. The optimization parameters are as shown in Table 2.

In traditional control, there are two control methods for disturbance suppression: feedforward & feedback and cascade control. Cascade control is one of the most popular complex control structures that can be found in the process industries, implemented to improve the disturbance rejection properties of the controlled system [9, 10, 35]. As shown in Figure 12, the system’s goal is to control the outlet moisture of the tobacco leaf. The manipulated variable is the flow rate of the mixed water vapor, and the disturbance variable is the inlet pressure of the mixed water vapor. is the primary controller controlling the outlet moisture of tobacco leaf, and is the subcontroller controlling the mixed water vapor inlet pressure. is the inlet pressure sensor and transmitter of mixed water vapor. Here, we use cascade control with the combined integrating controller (CIC), PPI, IMC PID method, ZieglerNichols method, and Lambda method, respectively, in Figure 12. We compare the performance of the five algorithms based on three situations.
8.1. Case 1
For the nominal process model and disturbance model that we have established, we control the process object through the combined integrating control (CIC), PPI, IMC PID, ZieglerNichols method, and Lambda method. From the setpoint tracking response in Figure 13, we can see that the CIC, PPI, and IMC PID control algorithms reach the setpoint value within 2000 seconds, while the ZieglerNichols and Lambda control algorithms fail to track the setpoint value, a meaningless control algorithm. From Table 3, the overall performance indexes of the three effective control algorithms are counted.

The CIC control algorithm proposed in this paper is superior to PPI and IMC in terms of response time dynamic performance index and steadystate performance index (rise time, settling time, and steadystate error). The ITAE performance index from the system error also confirms that the CIC control performance is improved twice as much as the PPI and IMC.
8.2. Case 2
As another disturbance of the system, the hot air flow affects the moisture at the tobacco leaf outlet. When the hot air flow fluctuates positively, it will enhance the tobacco leaf humidification process and increase the tobacco leaf outlet moisture. We assume that, at , a disturbance model is added. The input of the first scheme is a signal with a positive step amplitude of 1, accompanied by Gaussian random noise (the mean of 0, the variance of 0.5). The input of the second scheme is a signal with a negative step amplitude of 1, accompanied by Gaussian random noise (the mean of 0, the variance of 0.5).
From the setpoint tracking response in Figure 14, we can see that the CIC, PPI, and IMC PID control algorithms reached the setpoint value within the first 2000 seconds, while the ZieglerNichols and Lambda control algorithms failed to track the setpoint value. At , a positive disturbance and a negative disturbance are added. The five control algorithms actively suppress the disturbance and retrack the setpoint value. After adding disturbance to the ZieglerNichols and Lambda control algorithms, the control performance has been improved, but there are still steadystate errors. Table 4 summarizes the system error performance index ITAE of three effective control algorithms in positive disturbance and negative disturbance.
(a)
(b)

The CIC control algorithm proposed in this paper is superior to PPI and IMC in terms of timedomain dynamic performance index and steadystate performance index (perturbation amount, recovery time, and steadystate error) that suppress disturbance response. The system error performance index ITAE in Table 4 confirms that the CIC control performance is 1.5 times higher than PPI and IMC. The ZieglerNichols and Lambda control performance have been improved, indicating that the disturbance does not always deteriorate the control performance. The above example proves that the control performance can be improved.
8.3. Case 3
There is model adaptation when the system gets a wrong model due to data acquisition error or sensor equipment damage. Here, we assume that there is no mismatch in the system’s disturbance model and that the process model of the system changes from FOPDT model to the combined integrating system model, which is a large timedelay system; that is, , , .
When the system model does not match, the ZieglerNichols and Lambda control algorithms cannot stabilize the system; from the setpoint value tracking response of Figure 15, we can see that the CIC, PPI, and IMC PID control algorithms can all stabilize the system; the timedomain dynamic performance of the system (rise time, settling time, and overshoot) has declined. Table 5 summarizes the system error performance index ITAE of three effective control algorithms under model mismatch.

The CIC control algorithm proposed in this paper is superior to PPI and IMC in suppressing the timedomain dynamic performance index and steadystate performance index (rise time, settling time, and steadystate error) of the model mismatch response. ITAE from the system error performance index confirms that CIC control performance is improved two times as PPI and IMC.
The statistical data in the case analysis verifies that the CIC control algorithm proposed in this paper is superior to the traditional control algorithm in terms of timedomain dynamic performance, steadystate performance, and system error performance under setpoint tracking, disturbance suppression, and model mismatch. On the other hand, it is worth noting that modeling errors and measurement errors and disturbances do not always deteriorate the control performance. The above examples prove that the control performance can be improved.
9. Conclusion
This paper proposes a combined integrating controller claiming its applications and characteristics for the industry’s linear timeinvariant system. Extensive simulations were carried out to compare the performance of the proposed controller with popular controllers such as PID and Predictive PI control algorithms for conventional industrial processes. From the simulations, we see that the proposed controller has some superior characteristics and indexes. The proposed controller provides a more flexible framework due to the controller’s output only being determined by the system input and the average value of the controller’s output over a period. Further, the combined integrating term and the timedelay term are introduced into the controller structure, which breaks the disadvantage of the timedelay term. Overall, the proposed controller provides an attractive control alternative to the conventional controller.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare no conflicts of interest.
Authors’ Contributions
Angang C and Zhengyun R contributed to conceptualization, methodology, original draft preparation, and financial support. Zhiping F and Xue F contributed to validation. Angang C contributed to final edition check.
Acknowledgments
This research work was supported by the Fundamental Research Funds of the Central Universities (CUSFDHD2018100) and the China Scholarship Council (201806630067).
References
 Astrom, J. Karl, and H. Tore, PID Controllers: Theory, Design, and Tuning, ISAThe Instrumentation, Systems and Automation Society, Durham, NC, USA, 1995.
 G. Fedele, “A new method to estimate a firstorder plus time delay model from step response,” Journal of the Franklin Institute, vol. 346, no. 1, pp. 1–9, 2009. View at: Publisher Site  Google Scholar
 B. M. Sarif, D. V. A. Kumar, and M. V. G. Rao, “Comparison study of PID controller tuning using classical analytical methods,” International Journal of Applied Engineering Research, vol. 13, pp. 5618–5625, 2018. View at: Google Scholar
 M. Parada, R. A. Borges, D. Sbarbaro, and P. L. D. Peres, “Robust PID design for secondorder processes with timedelay and structured uncertainties,” IFAC Proceedings Volumes, vol. 44, no. 1, pp. 4614–4619, 2011. View at: Publisher Site  Google Scholar
 Q. H. Seer and J. Nandong, “Stabilization and PID tuning algorithms for secondorder unstable processes with timedelays,” ISA Transactions, vol. 67, pp. 233–245, 2017. View at: Publisher Site  Google Scholar
 M. Ajmeri and A. Ali, “Simple tuning rules for integrating processes with large time delay,” Asian Journal of Control, vol. 17, no. 5, pp. 2033–2040, 2015. View at: Publisher Site  Google Scholar
 P. Larsson and T. Hägglund, “Comparison between robust PID and predictive PI controllers with constrained control signal noise sensitivity,” IFAC Proceedings Volumes, vol. 45, no. 3, pp. 175–180, 2012. View at: Publisher Site  Google Scholar
 K. J. Åström and T. Hägglund, “Revisiting the ZieglerNichols step response method for PID control,” Journal of Process Control, vol. 14, no. 6, pp. 635–650, 2004. View at: Publisher Site  Google Scholar
 M. R. Cesca and J. L. Marchetti, “IMC design of cascade control,” Computer Aided Chemical Engineering, Elsevier, Amsterdam, The Netherlands, 2005. View at: Publisher Site  Google Scholar
 C. V. N. Rao and R. P. Sree, “IMC based controller design for integrating systems with time delay,” Indian Chemical Engineer, vol. 52, no. 3, pp. 194–218, 2012. View at: Publisher Site  Google Scholar
 Y. Gohda, S. Masuda, and Y. Matsui, “A direct PID gain tuning based on FRIT method using optimal filter in the frequency domain,” IFAC Proceedings Volumes, vol. 46, no. 11, pp. 349–354, 2009. View at: Publisher Site  Google Scholar
 Š. Bucz and A. Kozáková, “The sine wave tuning method: robust PID controller design in the frequency domain,” Journal of Electrical Systems and Information Technology, vol. 2, no. 3, pp. 358–367, 2015. View at: Publisher Site  Google Scholar
 E. Grassi and K. Tsakalis, “PID controller tuning by frequency loopshaping: application to diffusion furnace temperature control,” IEEE Transactions on Control Systems Technology, vol. 8, no. 5, pp. 842–847, 2000. View at: Publisher Site  Google Scholar
 G. Fedele, “Independent design of multiloop PI/PID controllers for interacting multivariable processes,” Journal of Process Control, vol. 20, no. 8, pp. 922–933, 2010. View at: Google Scholar
 J. Nandong and Z. Zang, “Multiloop design of multiscale controllers for multivariable processes,” Journal of Process Control, vol. 24, no. 5, pp. 600–612, 2014. View at: Publisher Site  Google Scholar
 R. Zhang, Q. Zou, Z. Cao, and F. Gao, “Design of fractional order modeling based extended nonminimal state space MPC for temperature in an industrial electric heating furnace,” Journal of Process Control, vol. 56, pp. 13–22, 2017. View at: Publisher Site  Google Scholar
 G. S. Virk, M. Morari, and E. Zafiriou, “Robust process control,” Chemical Engineering Science, vol. 45, no. 4, p. 1149, 1990. View at: Publisher Site  Google Scholar
 F. Xie and H. Lu, “Delaydependent robust Hinfinity stability for uncertain discretetime Markovian jump singular system,” in Proceedings of the Chinese Control and Decision Conference, Guilin China, June 2009. View at: Publisher Site  Google Scholar
 L. T. Biegler, G. S. Virk, M. Morari, and E. Zafiriou, “Solution of dynamic optimization problems by successive quadratic programming and orthogonal collocation,” Computers & Chemical Engineering, vol. 8, no. 34, pp. 243–247, 1984. View at: Publisher Site  Google Scholar
 M. Cizniar, D. Salhi, M. Fikar, and E. Abderrazak, “A MATLAB package for orthogonal collocations on finite elements in dynamic optimisation,” in Proceedings of the 15th International Conference Process Control, Strbske Pleso, Slovakia, June 2015. View at: Google Scholar
 L. T. Biegler, “An overview of simultaneous strategies for dynamic optimization,” Chemical Engineering and Processing: Process Intensification, vol. 46, no. 11, pp. 1043–1053, 2007. View at: Publisher Site  Google Scholar
 S. Kameswaran and L. T. Biegler, “Convergence rates for direct transcription of optimal control problems using collocation at Radau points,” Computational Optimization and Applications, vol. 41, no. 1, pp. 81–126, 2008. View at: Publisher Site  Google Scholar
 J. D. Hedengren, R. A. Shishavan, K. M. Powell, and T. F. Edgar, “Nonlinear modeling, estimation and predictive control in APMonitor,” Computers & Chemical Engineering, vol. 70, pp. 133–148, 2014. View at: Publisher Site  Google Scholar
 D. Chen and D. E. Seborg, “PI/PID controller design based on direct synthesis and disturbance rejection,” Industrial & Engineering Chemistry Research, vol. 41, no. 19, pp. 4807–4822, 2002. View at: Publisher Site  Google Scholar
 N.S. Pai, S.C. Chang, and C.T. Huang, “Tuning PI/PID controllers for integrating processes with deadtime and inverse response by simple calculations,” Journal of Process Control, vol. 20, no. 6, pp. 726–733, 2010. View at: Publisher Site  Google Scholar
 A. A. Khan, J. R. Moyne, and D. M. Tilbury, “Favorable effect of time delays on tracking performance of typeI control systems,” IET Control Theory & Applications, vol. 2, no. 3, pp. 210–218, 2008. View at: Publisher Site  Google Scholar
 G. Szita and C. K. Sanathanan, “Robust design for disturbance rejection in time delay systems,” Journal of the Franklin Institute, vol. 334, no. 4, pp. 611–629, 1997. View at: Publisher Site  Google Scholar
 W. Tan, H. J. Marquez, and T. Chen, “Performance assessment of PID controllers,” Control and Intelligent Systems, vol. 32, no. 3, 2004. View at: Publisher Site  Google Scholar
 X.N. Zhang and G.H. Yang, “Performance analysis for multidelay systems in finite frequency domains,” International Journal of Robust and Nonlinear Control, vol. 22, no. 8, pp. 933–944, 2012. View at: Publisher Site  Google Scholar
 C. Dalen and D. Di Ruscio, “Performance optimal PI controller tuning based on integrating plus time delay models,” Algorithms, vol. 11, no. 6, pp. 86–109, 2018. View at: Publisher Site  Google Scholar
 T. N. L. Vu and M. Lee, “Multiloop PI controller design based on the direct synthesis for interacting multitime delay processes,” ISA Transactions, vol. 49, no. 1, pp. 79–86, 2010. View at: Publisher Site  Google Scholar
 D. Zheng, Z. Ren, and J.a. Fang, “Stability analysis of multiple timedelayed system,” ISA Transactions, vol. 47, no. 4, pp. 439–447, 2008. View at: Publisher Site  Google Scholar
 D. Zheng, Z. Ren, and J. Fang, “Stability analysis of time delayed system with coefficient uncertainty and time delay uncertainty,” European Journal of Control, vol. 16, no. 1, pp. 5–13, 2010. View at: Publisher Site  Google Scholar
 Q. Zhang, F.F. Lee, YG Wang, W. Zhang, and Q. Chen, “Performance assessment of cascade control system with nonGaussian disturbance based on minimum entropy,” Symmetry, vol. 11, no. 3, pp. 379–394, 2019. View at: Publisher Site  Google Scholar
 O. Arrieta, R. Vilanova, and P. Balaguer, “Procedure for cascade control systems design: choice of suitable PID tunings,” International Journal of Computers Communications & Control, vol. 3, no. 3, pp. 235–248, 2008. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Angang Chen et al. 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.