Next Article in Journal
A Multi-Communication-Based Demand Response Implementation Structure and Control Strategy
Previous Article in Journal
Simplified Regional Prediction Model of Long-Term Trend for Critical Frequency of Ionospheric F2 Region over East Asia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Discrete-Time First-Order Plus Dead-Time Model-Reference Trade-off PID Control Design

1
Department of Mechanical Engineering, Graduate School of Engineering, University of Hyogo, Hyogo 671-2280, Japan
2
Department of Telecommunications and Systems Engineering, Universitat Autònoma de Barcelona, 08193 Barcelona, Spain
*
Author to whom correspondence should be addressed.
Appl. Sci. 2019, 9(16), 3220; https://doi.org/10.3390/app9163220
Submission received: 19 July 2019 / Revised: 3 August 2019 / Accepted: 5 August 2019 / Published: 7 August 2019
(This article belongs to the Section Electrical, Electronics and Communications Engineering)

Abstract

:
The present study proposes a novel proportional-integral-derivative (PID) control design method in discrete time. In the proposed method, a PID controller is designed for first-order plus dead-time (FOPDT) systems so that the prescribed robust stability is accomplished. Furthermore, based on the control performance, the relationship between the servo performance and the regulator performance is a trade-off relationship, and hence, these items are not simultaneously optimized. Therefore, the proposed method provides an optimal design method of the PID parameters for optimizing the reference tracking and disturbance rejection performances, respectively. Even though such a trade-off design method is being actively researched for continuous time, few studies have examined such a method for discrete time. In conventional discrete time methods, the robust stability is not directly prescribed or available systems are restricted to systems for which the dead-time in the continuous time model is an integer multiple of the sampling interval. On the other hand, in the proposed method, even when a discrete time zero is included in the controlled plant, the optimal PID parameters are obtained. In the present study, as well as the other plant parameters, a zero in the FOPDT system is newly normalized, and then, a universal design method is obtained for the FOPDT system with the zero. Finally, the effectiveness of the proposed method is demonstrated through numerical examples.

1. Introduction

Proportional-integral-derivative (PID) [1,2,3,4,5,6,7,8,9] control has few tuning parameters: proportional gain, integral time, and derivative time, and its structure is simple. Hence, PID control has been widely used in industry, and numerous tuning methods have been proposed.
In the model-based approach, optimal tracking and robust stability are achieved. However, the control system must be redesigned whenever the controlled system is changed. As a simple model-based optimal design method, Ziegler and Nichols proposed the step response method (ZN method) [10]. Using the ZN method, the PID parameters are decided based on the step response trajectory so that the tracking performance is optimized. However, the robust stability is not taken into account [11,12]. Therefore, robust control designs have been proposed [13,14]. Using robust design methods, robust stability is obtained, but the tracking performance is not optimized. In [15], H 2 optimal design was proposed using internal model control (IMC) design, and hence, stability was assured. Although robust PID control systems are also designed in the discrete time domain [16,17], the stability margin cannot be prescribed.
The relationship between the tracking performance and the robust stability is a trade-off relationship [18]. Therefore, trade-off PID design methods have been proposed [19,20]. In the trade-off design methods, the PID parameters are decided such that the tracking performance is optimized subject to the prescribed robust stability, where servo or regulator mode is selected for performance optimization. A trade-off optimization approach has been designed in the continuous time domain that consists of a first-order plus dead-time (FOPDT) system [21,22], a second-order plus dead-time (SOPDT) system [23,24], and a two degrees-of-freedom (2DOF) system [25,26].
For a discrete time system, trade-off approaches have been proposed [27,28]. Discrete-time methods are useful for controlling discrete time systems, in which the controller is implemented with a digital computer. However, the robust stability was not assigned in [27], and the controlled plant was restricted to a non-zero system [28]. Therefore, in the present study, a new discrete time domain approach is proposed, in which the constraint is relaxed and the prescribed robust stability is accomplished. Specifically, the FOPDT model with a zero is normalized, and a universal design method is obtained. As a result, the trade-off design strategy is available for a large class of discrete time systems. Finally, the effectiveness of the proposed method is demonstrated through numerical examples.

2. Control System

The continuous time FOPDT transfer function is described as follows:
P ( s ) = K T s + 1 e L s
where K is the gain, T is the time constant, and L is the dead time. In the present study, the discrete time control model is designed with sampling interval T s , and hence, the continuous time model is expressed as follows:
P d ( z 1 ) = b 0 + b 1 z 1 1 a 1 z 1 z ( d + 1 )
where a 1 , b 0 , and b 1 are the coefficient parameters, d is the dead-time in the discrete time, and z 1 denotes the backward shift operator. The discrete time plant parameters correspond to the continuous time plant parameters as follows:
a 1 = e T s T
b 0 = K 1 a 1 e L 0 T
b 1 = K a 1 e L 0 T a 1
d = L 1
where a non-negative integer L 1 satisfies the following equation:
L = L 1 T s + L 0 ( L 0 < T s )
Therefore, the input/output relationship in discrete time is given as follows:
y ( k ) = P d ( z 1 ) u ( k )
where y ( k ) is the system output (plant output) and u ( k ) is the control input.
The present study proposes a new trade-off design method for the following PID control law:
u ( k ) = C e ( z 1 ) e ( k ) C y ( z 1 ) y ( k ) C e ( z 1 ) = K p 1 + T s T i ( 1 z 1 ) C y ( z 1 ) = K p T d ( 1 z 1 ) T s e ( k ) = r ( k ) y ( k )
where r ( k ) is the reference input and K P , T i , and T d are the proportional gain, the integral time, and the derivative time, respectively. A block diagram of the control system is shown in Figure 1, in which the control input is disturbed by disturbance d ( k ) . Equation (8) is the discrete version of the next continuous time control law:
U ( s ) = K p 1 + 1 T i s E ( s ) T d s Y ( s ) E ( s ) = R ( s ) Y ( s )
where L [ · ] means the Laplace transform, U ( s ) = L [ u ( t ) ] , Y ( s ) = L [ y ( t ) ] , and R ( s ) = L [ r ( t ) ] .
The closed-loop systems from r ( k ) and d ( k ) to u ( k ) and y ( k ) , respectively, are obtained as follows:
u ( k ) = C e ( z 1 ) 1 + C d ( z 1 ) P d ( z 1 ) r ( k ) C d ( z 1 ) P d ( z 1 ) 1 + C d ( z 1 ) P d ( z 1 ) d ( k )
y ( k ) = C e ( z 1 ) P d ( z 1 ) 1 + C d ( z 1 ) P d ( z 1 ) r ( k ) + P d ( z 1 ) 1 + C d ( z 1 ) P d ( z 1 ) d ( k ) C d ( z 1 ) = C e ( z 1 ) + C y ( z 1 )
From the above equations, the relationship between the reference tracking and the disturbance rejection is a trade-off relationship.

3. Design Objective

In the proposed method, the PID parameters of the discrete time control law are designed based on a constrained optimization problem. The constrained optimization problem consists of the constraint condition as the prescribed stability margin and the objective function defined by the index with respect to the servo and regulator performances, respectively. The constraint condition and the objective function are defined in Section 3.1 and Section 3.2, respectively.

3.1. Constraint Condition

The constraint condition is used to obtain the stability margin. In the proposed method, the stability margin is designed using the following sensitivity function:
S f ( z 1 ) = 1 1 + C d ( z 1 ) P d ( z 1 )
Using the sensitivity function, the constraint condition is defined as follows:
| M s M s d | = 0
M s = max ω | S f ( e j ω ) |
where M s denotes the maximum value of the sensitivity function S f ( e j ω ) and M s d is the designed value of M s .
The design range of M s d is recommended to be from 1.4 2.0 [4]. The relationship between M s d and the stability margin is in inverse proportion.

3.2. Objective Function

The objective function is defined based on the sum of absolute errors (SAE) as follows:
J d = k = 1 | r ( k ) y ( k ) |
In a one-degree-of-freedom (1DOF) system, the relationship between the servo and regulation performances is a trade-off relationship. Therefore, the objective functions in the reference tracking optimization and the disturbance rejection optimization are given by Equations (16) and (17), respectively:
J d s = k = 1 r s ( k ) C e ( z 1 ) P d ( z 1 ) 1 + C d ( z 1 ) P d ( z 1 ) r ( k )
J d r = k = 1 r d ( k ) P d ( z 1 ) 1 + C d ( z 1 ) P d ( z 1 ) d ( k )
where r s ( k ) = 1 , r ( k ) = 1 , r d ( k ) = 0 , and d ( k ) = 1 .
In the discussed control system, the stability and the tracking performance are simultaneously decided because they are in a trade-off relationship.

4. Proportional-Integral-Derivative Parameter Design

In order to obtain a universal design method for arbitrary FOPDT models, first, in Section 4.1, the plant model and the PID control law are normalized. Second, in Section 4.2, optimal PID parameters are obtained subject to the pre-established stability margin. Third, in Section 4.3, the parameter tuning rule is derived based on the optimal parameters calculated in Section 4.2. Finally, the algorithm for the proposed design method is summarized in Section 4.4.

4.1. Normalization

The plant model and the control law are normalized to obtain a universal design method. The normalization parameters for the plant model, τ 0 and τ a , are given as follows:
τ 0 = d log a 1 + log b 0 a 1 + b 1 a 1 ( b 0 + b 1 )
τ a = log a 1
Moreover, the normalization parameters for the control law, κ p , τ i , and τ d , are given as follows:
κ p = b 0 + b 1 1 a 1 K p
τ i = T i log a 1 T s
τ d = T d log a 1 T s
The derivation of the normalization parameters is shown in Appendix A.
The expression of the plant model using the normalization parameters has merits for optimization analysis and controller design. Since the optimization problem is stated in terms of just two parameters τ 0 and τ a , the analysis and elaboration of coefficient parameters are simplified. Furthermore, the controller parameters are tuned independent of the process gain and time constant.

4.2. Optimization

The controller parameters are optimized for the servo and regulation operations, respectively, subject to the established stability margin.
As the optimization tool, the fmincon function in MathWorks MATLAB software was used, where the prescribed robust stability is set: M s d { 1.4 , 1.6 , 1.8 , 2.0 } , and the optimal controller parameters are calculated based on the normalization parameters for the controlled plant: τ 0 { 0.3 , , 1.7 } and τ a { 0.01 , , 0.1 } . As an example, the optimized controller parameters for M s d = 1.4 in the servo are shown in Figure 2.

4.3. Controller Parameter Decision

The tuning rule for the controller parameters is proposed in Section 4.3.1, and the prescribed robust stability using the proposed parameter tuning rule is evaluated in Section 4.3.2.

4.3.1. Tuning Rule

On the pair of plant parameters τ 0 and τ a with the established robust stability M s d , the controller parameters are expressed in terms of the normalization parameters for the plant model as:
κ p = α 0 + α 1 τ 0 α 2
τ i = β 0 + β 1 τ 0 + β 2 τ 0 2 + β 3 τ 0 3
τ d = γ 0 + γ 1 τ 0 + γ 2 τ 0 2
where α i ( i = 0 , 1 , 2 ) , β j ( j = 0 , 1 , 2 , 3 ) , and γ k ( k = 0 , 1 , 2 ) are defined as follows:
α 0 = a 00 + a 01 τ a
α 1 = a 10 + a 11 τ a
α 2 = a 20 + a 21 τ a
β 0 = b 00 + b 01 τ a
β 1 = b 10 + b 11 τ a
β 2 = b 20 + b 21 τ a
β 3 = b 30 + b 31 τ a
γ 0 = c 00 + c 01 τ a
γ 1 = c 10 + c 11 τ a
γ 2 = c 20 + c 21 τ a
where the coefficient parameters in Equation (26)–Equation (35) in the servo optimization are shown in Table 1, and those in the regulation optimization are also shown in Table 2.
The controller parameters are optimized in the limited range. Although the range can be expanded, the tuning rule is redesigned when the range is changed.

4.3.2. Evaluation of the Prescribed Robust Stability

Non-obtained controller parameters are interpolated, where using the tuning rule proposed in Section 4.3.1, the controller parameters are calculated on the pair of the normalization parameters for the controlled plant: τ 0 { 0.3 , 0.31 , , 1.7 } and τ a { 0.01 , 0.011 , , 1.0 } , where M s d is varied in { 1.4 , ⋯, 2.0 }. The calculated controller parameters for M s d = 1.4 , merged into the preliminary calculated parameters shown in Figure 2, are shown in Figure 3. The calculated parameters using the proposed tuning rule are sufficiently close to the optimized parameters.
Next, M s values calculated using the calculated controller parameters in the servo and regulation optimization are shown in Figure 4 and Figure 5, respectively. The minimum and maximum calculated M s values are summarized in Table 3. Since the errors between M s d and M s are within ± 5 % , the proposed tuning rule has sufficient precision.

4.4. Algorithm

The proposed design procedure is summarized as the following algorithm:
  • The plant parameters and the controller parameters are normalized.
  • For the normalized systems, the constrained optimization problem is solved such that the reference and disturbance responses are optimized, respectively, and the optimal controller parameters are obtained.
  • Based on the obtained optimized parameters, the tuning rule for the controller parameters is designed.
  • Using the tuning rule, the practical PID parameters are decided.
When the tuning rule is obtained once, the PID parameters are calculated without solving the constrained optimization problem.

5. Numerical Examples

First, in Section 5.1, the trade-off design between the servo and regulator design is confirmed. The accomplished robust stability is then shown in Section 5.2. Finally, the proposed method is compared with two conventional discrete time trade-off design methods [27,28] in Section 5.3 and Section 5.4, respectively.

5.1. Trade-off Tracking Performance Comparison

As a controlled plant, the following transfer function is used:
P d ( z 1 ) = 0.0231 + 0.0114 z 1 1 0.9753 z 14
The transfer function is the discrete time representation of the following continuous time system with T s = 0.03 s:
P ( s ) = 1.4 1.2 s + 1 e 0.4 s
Using the proposed method, the PID parameters are decided based on Equation (36) with M s d { 1.4 , 1.6 , 1.8 , 2.0 } . The simulation results obtained using the PID parameters are shown in Figure 6, where the reference input is 1.0 and the control input is disturbed by a unit step signal after 15 s. Figure 6 shows that the servo design is superior to the regulator design with respect to reference tracking performance. On the other hand, the regulator design is superior to the servo design with respect to disturbance rejection.
The control performance is evaluated using J s and J r , where J s denotes the SAE value from the start until 15 s, and J r also denotes the SAE value from 15 s until the end. Here, J s and J r , as well as the PID parameters and M s in the servo and regulator, respectively, are summarized in Table 4. In the servo design, J s is smaller than J r . In the regulator design, J r is smaller than J s .
From both Figure 6 and Table 4, the larger M s d , the better the control performance. Therefore, the trade-off design is achieved using the proposed method.

5.2. Robust Stability

The effectiveness of the accomplished robust stability is shown. The PID controller is designed based on the nominal model Equation (36) in Section 5.1 in the servo and regulation, respectively. The scenario is that the plant is Equation (36) from the start to 30 s and is changed to Equation (38) after 30 s as the model perturbation.
P d ( z 1 ) = 0.0530 1 0.9788 z 21
The transfer function is the discrete time version of the following continuous time model with the sampling interval T s = 0.03 s:
P ( s ) = 2.5 1.4 s + 1 e 0.6 s
The simulation results for the servo and regulator designs are shown in Figure 7, where the reference input is given by a unit step function, and the control input is disturbed by a unit step function after 15 s. The simulation results show that the effect of the model perturbation is suppressed by small M s d . However, note that the control performance is superior before the model perturbation when the value of M s d is large.

5.3. Comparison with the IMC-Based Method

Consider the following discrete time system:
P d ( z 1 ) = 0.0100 1 0.99 z 1 z 25
The system is the discrete time representation of the following continuous time system with a sampling interval of 0.01 s:
P ( s ) = 1 s + 1 e 0.25 s
In the conventional discrete time trade-off design method [27], the IMC controller is designed as:
C I M C ( z 1 ) = ( 1 a 1 ) ( 1 λ c ) b 0 ( 1 λ c z 1 )
where λ c is the trade-off design parameter. The obtained IMC controller is approximated by the following discrete time PID control law:
C P I D ( z 1 ) = K p 1 + 1 T i ( 1 z 1 ) + ( 1 α T d ) T d ( 1 z 1 ) 1 α T d z 1
where λ c is designed in the range of [0.8,0.99]. However, the designed control systems are unstable when λ c 0.96 , and thus, the conventional control law is designed such that λ c is set to 0.97 , 0.98 , and 0.99 , respectively, where α = 0.1 . The designed PID parameters for the conventional method and the proposed method are shown in Table 5.
The simulations were conducted using the conventional and proposed methods, where the reference input was set to 1.0 , and the control input was disturbed by a unit step function after 30 s. The conventional method is compared with the proposed servo and regulator optimization methods in Figure 8. The obtained M s values and the evaluated values J s and J r are also summarized in Table 5, where J s denotes the SAE value from the start until 30 s, and J r also denotes the SAE value from 30 s until the end. Table 5 shows that the conventional method provides a trade-off design by selecting λ c even though no value is assigned to M s . Moreover, the tracking performances obtained using the proposed method are superior to those obtained using the conventional method.

5.4. Comparison with the Conventional Discrete-Time Method

The conventional discrete time design method [28] is compared with the proposed method. Here, two scenarios are conducted, in which a non-zero plant and a zero-included plant, respectively, are controlled.
In the first simulation, we consider the following non-zero discrete time plant:
P d 1 ( z 1 ) = 0.05127 z 1 1 0.9486 z 1 z 10
where Equation (44) is the discrete time representation of Equation (45) with a sampling interval of T s = 0.05 s.
P 1 ( s ) = 1 0.95 s + 1 e 0.5 s
Equation (44) has no zero since the continuous time dead-time is an integer multiple of the sampling interval.
In the second simulation, the controlled discrete time system is given as follows:
P d 2 ( z 1 ) = 0.0201 + 0.02473 z 1 1 0.9552 z 1 z 7
where Equation (46) is the discrete time version of the continuous time system given by Equation (47) with a sampling interval of T s = 0.061 s.
P 2 ( s ) = 1 1.33 s + 1 e 0.4 s
Equation (46) has a zero since the dead-time in the continuous time model is not an integer multiple of the sampling interval. Since the discrete time system has a zero, the conventional method is not directly used. Therefore, Equation (46) is hereby approximated by the next discrete time system, and the conventional method is used:
P d 2 ( z 1 ) = 0.04483 z 1 1 0.9552 z 1 z 7
In the simulations, the reference input is set to 1.0 . Furthermore, the control input is disturbed by a unit step function signal after 10 s. The obtained PID parameters are shown in Table 6, Table 7, Table 8 and Table 9. Using the obtained parameters, the discrete time models Equations (44) and (46) are controlled, respectively, and the output results are plotted in Figure 9 and Figure 10. Furthermore, the obtained M s value and index values J s and J r are also shown in Table 6, Table 7, Table 8 and Table 9, where J s denotes the SAE value while the control is not disturbed, and J r denotes the SAE value while the control input is disturbed.
In the case of the non-zero system Equation (44), J s , J r , and M s obtained using the conventional method were close to those obtained using the proposed method. On the other hand, when the controlled plant was a zero-included system, Equation (46) was out of range of the conventional method. Therefore, the tracking performances using the conventional method were inferior to those using the proposed method, and hence, the SAE values using the conventional method were larger than those using the proposed method. Furthermore, the prescribed robust stabilities M s d were achieved using the proposed method, while on the other hand, the M s values obtained using the conventional method were insufficient.
The simulation results showed that both the conventional and proposed methods were useful for the non-zero plants. However, when the zero-included plant was controlled, the design objective was achieved using the proposed method even though the conventional method was not available. Therefore, the proposed method is a more general method of the conventional discrete time design method.

6. Conclusions

The present study proposed a new trade-off PID control design method for discrete time FOPDT systems including a zero. In the proposed method, the regulator- or servo-optimal PID controller was designed in discrete time. In the proposed method, since robust stability was a design parameter, it was adjustable depending on the model perturbation.
In the conventional discrete time design method [28], the designable class was restricted such that the dead-time in the continuous time system must be an integer multiple of the sampling interval, and hence, no zero appeared. On the other hand, in the proposed method, the constraint condition was relaxed.

Author Contributions

Conceptualization, R.K. and T.S.; methodology, R.K., T.S., and R.V.; software, R.K.; validation, T.S. and R.V.; formal analysis, R.K.; investigation, T.S. and R.V.; resources, R.K.; data curation, R.K. and T.S.; writing, original draft preparation, T.S.; writing, review and editing, R.V.; visualization, R.K., T.S., and R.V.; supervision, T.S. and R.V.; project administration, T.S., R.V., and Y.K.; funding acquisition, T.S., R.V., and Y.K.

Acknowledgments

The present study was supported by JSPS KAKENHI Grant Number 16K06425. This research is supported by the Catalan Government under Project 2017 SGR 1202 and also by the Spanish Government under Project DPI2016-77271-R co-funded with European Union ERDF funds.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Derivation of Normalization Parameters

In the same way as the conventional continuous time design, the normalization parameters are obtained.
Here, τ 0 is derived as follows:
τ 0 = L T = T s T d + L 0 T
In this equation, T s T is derived using Equation (3):
T s T = log a 1
and L 0 T is also derived by eliminating K from Equations (4) and (5):
L 0 T = log b 0 + a 1 ( b 0 + b 1 ) a 1
Equations (A2) and (A3) are substituted into Equation (A1), and Equation (18) is then obtained.
From Equation (A2), τ a is derived as:
τ a = T s T = log a 1
Since K is the steady-state gain, κ p is obtained as:
κ p = K K p = b 0 + b 1 1 a 1 K p
From Equation (A2), τ i and τ d are calculated as follows:
τ i = T i T = T i log a 1 T s
τ d = T d T = T d log a 1 T s

References

  1. O’Dwyer, A. Handbook of PI and PID Controller Tuning Rules; Imperial College Press: London, UK, 2003. [Google Scholar]
  2. Silva, G.J.; Datta, A.; Bhattacharyya, S. PID Controllers for Time-Delay Systems; Birkhäuser: Boston, MA, USA, 2005. [Google Scholar]
  3. Johnson, M.A.; Moradi, M.H. (Eds.) PID Control: New Identification and Design Methods; Springer-Verlag: London, UK, 2005. [Google Scholar]
  4. Åström, K.J.; Hägglund, T. Advanced PID Control; Instrumentation, Systems, and Automation Society: Pittsburgh, PA, USA, 2006. [Google Scholar]
  5. Visioli, A. Practical PID Control; Springer: London, UK, 2006. [Google Scholar]
  6. Yu, C.C. Autotuning of PID Controllers, 2nd ed.; Springer: London, UK, 2006. [Google Scholar]
  7. Vilanova, R.; Visioli, A. (Eds.) PID Control in the Third Millennium; Springer: London, UK, 2012. [Google Scholar]
  8. Boiko, I. Non-Parametric Tuning of PID Controllers; Springer: London, UK, 2013. [Google Scholar]
  9. Alfaro, V.; Vilanova, R. Model-Reference Robust Tuning of PID Controllers; Advances in Indstrial Control; Springer: Cham, Switzerland, 2016. [Google Scholar]
  10. Ziegler, J.; Nichols, N. Optimum settings for automatic controllers. Trans. ASME 1942, 64, 759–768. [Google Scholar] [CrossRef]
  11. Shinskey, F.G. How good are our controllers in absolute performance and robustness? Measur. Control 1990, 23, 114–121. [Google Scholar] [CrossRef]
  12. Shinskey, F.G. Process-Control Systems. Application, Design, and Tuning, 4th ed.; McGraw-Hill: New York, NY, USA, 1996. [Google Scholar]
  13. Chen, C.; Huang, H.; Hsieh, C. Tuning of PI/PID controllers based on specification on maximum closed-loop amplitude ratio. J. Chem. Eng. Jpn. 1999, 32, 783–788. [Google Scholar] [CrossRef]
  14. Matušů, R.; Prokop, R. Computation of robustly stabilizing PID controllers for interval systems. SpringerPlus 2016, 5, 702. [Google Scholar] [CrossRef] [PubMed]
  15. Jahanmiri, A.; Fallahi, H. New methods for process identification and design of feedback controller. Chem. Eng. Res. Des. 1997, 75, 519–522. [Google Scholar] [CrossRef]
  16. Ackermann, J.; Kaesbauer, D.; Bajcinca, N. Discrete-time robust PID and three-term control. IFAC Proc. Vol. 2002, 35, 127–132. [Google Scholar] [CrossRef]
  17. Hiroaki, M.; Takayoshi, U.; Yoshiyuki, T.; Toshio, T. Guaranteed cost PI control for uncertain discrete time systems with additive gain. In Proceedings of the 2009 European Control Conference (ECC), Budapest, Hungary, 23–26 August 2009; pp. 2319–2324. [Google Scholar]
  18. Garpinger, O.; Hägglund, T.; Åström, K.J. Criteria and trade-offs in PID design. IFAC Proc. Vol. 2012, 45, 47–52. [Google Scholar] [CrossRef]
  19. Alfaro, V.; Vilanova, R.; Méndez, V.; Lafuente, J. Performance/robustness tradeoff analysis of PI/PID servo and regulatory control systems. In Proceedings of the 2010 IEEE International Conference on Industrial Technology, Vina del Mar, Chile, 14–17 March 2010; pp. 111–116. [Google Scholar]
  20. Alcántara, S.; Pedret, C.; Vilanova, R. On the model matching approach to PID design: Analytical perspective for robust servo/regulator tradeoff tuning. J. Process Control 2010, 20, 596–608. [Google Scholar] [CrossRef]
  21. Arrieta, O.; Vilanova, R. Simple servo/regulation Proportional-Integral-Derivative (PID) tuning rules for arbitrary Ms-based robustness achievement. Ind. Eng. Chem. Res. 2012, 51, 2666–2674. [Google Scholar] [CrossRef]
  22. Kurokawa, R.; Inoue, N.; Sato, T.; Arrieta, O.; Vilanova, R.; Konishi, Y. Simple optimal PID tuning method based on assgined robust stability—Trade-off design based on servo/regulation performance. Int. J. Innov. Comput. Inf. Control 2017, 13, 1953–1963. [Google Scholar]
  23. Alfaro, V.M.; Méndez, V.; Vilanova, R. Robust-performance tuning of 2DoF PI/PID controllers for first- and second-order-plus-dead-time models. In Proceedings of the 2011 9th IEEE International Conference on Control and Automation (ICCA), Santiago, Chile, 19–21 December 2011; pp. 237–242. [Google Scholar]
  24. Sato, T.; Hayashi, I.; Horibe, Y.; Vilanova, R.; Konishi, Y. Optimal robust PID control for first- and second-order plus dead-time processes. Appl. Sci. 2019, 9, 1934. [Google Scholar] [CrossRef]
  25. Alfaro, V.; Vilanova, R. Robust tuning and performance analysis of 2DoF PI controllers for integrating controlled processes. Ind. Eng. Chem. Res. 2012, 51, 13182–13194. [Google Scholar] [CrossRef]
  26. Alfaro, V.; Vilanova, R. Model-reference robust tuning of 2DoF PI controllers for first- and second-order plus dead-time controlled processes. J. Process Control 2012, 22, 359–374. [Google Scholar] [CrossRef]
  27. Jiyao, C.; Yueling, C.; Tao, L. Discrete-time domain IMC-based PID control design for industrial processes with time delay. In Proceedings of the 2016 35th Chinese Control Conference (CCC), Chengdu, China, 27–29 July 2016; pp. 5946–5951. [Google Scholar]
  28. Tajika, H.; Sato, T.; Vilanova, R.; Konishi, Y. Optimal PID control in discrete time using sensitivity function. In Proceedings of the 2015 23rd Mediterranean Conference on Control and Automation (MED), Torremolinos, Spain, 16–19 June 2015; pp. 268–273. [Google Scholar]
Figure 1. Discrete-time PID control system.
Figure 1. Discrete-time PID control system.
Applsci 09 03220 g001
Figure 2. Relationships among optimized controller parameters and the normalized parameters τ 0 and τ a in servo design.
Figure 2. Relationships among optimized controller parameters and the normalized parameters τ 0 and τ a in servo design.
Applsci 09 03220 g002
Figure 3. Relationship between approximated controller parameters in servo design.
Figure 3. Relationship between approximated controller parameters in servo design.
Applsci 09 03220 g003
Figure 4. Relationships among the obtained M s , τ 0 , and τ a in servo design.
Figure 4. Relationships among the obtained M s , τ 0 , and τ a in servo design.
Applsci 09 03220 g004
Figure 5. Relationships among the obtained M s , τ 0 , and τ a in regulation design.
Figure 5. Relationships among the obtained M s , τ 0 , and τ a in regulation design.
Applsci 09 03220 g005
Figure 6. Output and input responses obtained using the proposed method in both servo and regulator design, where the system is disturbed by a unit step disturbance after 15 s.
Figure 6. Output and input responses obtained using the proposed method in both servo and regulator design, where the system is disturbed by a unit step disturbance after 15 s.
Applsci 09 03220 g006
Figure 7. Output and input responses obtained using the proposed method for each M s d with model perturbation, where the system is disturbed after 15 s; left: servo optimization; right: regulator optimization.
Figure 7. Output and input responses obtained using the proposed method for each M s d with model perturbation, where the system is disturbed after 15 s; left: servo optimization; right: regulator optimization.
Applsci 09 03220 g007
Figure 8. Results of the IMC-based method and the proposed servo and regulator optimization methods, where the system is disturbed by a unit step disturbance after 10 s; left: servo optimization; right: regulator optimization.
Figure 8. Results of the IMC-based method and the proposed servo and regulator optimization methods, where the system is disturbed by a unit step disturbance after 10 s; left: servo optimization; right: regulator optimization.
Applsci 09 03220 g008
Figure 9. Output and input trajectories for the proposed and conventional designs for Equation (44), where the system is disturbed by a unit step disturbance after 10 s; left: servo optimization; right: regulator optimization.
Figure 9. Output and input trajectories for the proposed and conventional designs for Equation (44), where the system is disturbed by a unit step disturbance after 10 s; left: servo optimization; right: regulator optimization.
Applsci 09 03220 g009
Figure 10. Output and input trajectories for the proposed and conventional designs for Equation (46), where the system is disturbed by a unit step disturbance after 10 s; left: servo optimization; right: regulator optimization.
Figure 10. Output and input trajectories for the proposed and conventional designs for Equation (46), where the system is disturbed by a unit step disturbance after 10 s; left: servo optimization; right: regulator optimization.
Applsci 09 03220 g010
Table 1. Coefficient parameters in Equations (26)–(35) in servo design.
Table 1. Coefficient parameters in Equations (26)–(35) in servo design.
M s d 1.4 1.6 1.8 2.0
a 00 0.21300.27780.32810.3098
a 01 −0.4643−0.6376−0.8185−0.7722
a 10 0.43610.58030.69320.8100
a 11 −0.3767−0.4236−0.3308−0.4577
a 20 −1.0067−1.0169−1.0150−0.9861
a 21 1.75091.79511.90031.8503
b 00 1.13681.14511.20971.3995
b 01 −1.6140−1.1310−0.7911−1.9403
b 10 −0.03940.31520.45160.1364
b 11 1.43930.0802−1.25932.0622
b 20 0.1724−0.0447−0.10940.1498
b 21 −0.92190.35211.6861−1.2358
b 30 −0.03260.02650.0354−0.0201
b 31 0.2070−0.1725−0.56770.2429
c 00 −0.01900.0000660.00470.0091
c 01 −0.1314−0.0898−0.0615−0.0129
c 10 0.31930.28190.33770.3596
c 11 0.33300.03810.03630.0514
c 20 0.0056−0.0100−0.0242−0.0090
c 21 −0.0527−0.01240.0078−0.0046
Table 2. Coefficient parameters in Equations (26)–(35) in regulator design.
Table 2. Coefficient parameters in Equations (26)–(35) in regulator design.
M s d 1.4 1.6 1.8 2.0
a 00 0.20850.27180.29990.3672
a 01 −0.6075−0.8871−0.6490−1.4148
a 10 0.44450.58970.72670.7914
a 11 −0.3597−0.3261−0.7568−0.1116
a 20 −1.0048−1.0010−0.9840−1.0107
a 21 2.42192.50222.17382.7688
b 00 0.21750.12080.16760.1793
b 01 1.01421.43500.51520.5668
b 10 1.30581.53591.44781.3845
b 11 −4.3025−4.9006−1.6551−1.4977
b 20 −0.7838−0.8310−0.6531−0.4397
b 21 3.78624.07340.99920.8169
b 30 0.22500.20670.15190.0589
b 31 −1.0977−1.1117−0.2245−0.1967
c 00 −0.00310.01390.01520.0314
c 01 0.08020.11030.07650.1761
c 10 0.44560.37830.36070.3006
c 11 0.33910.0800−0.0139−0.3791
c 20 −0.0467−0.0296−0.0374−0.0100
c 21 −0.1076−0.01070.01860.2333
Table 3. Maximum and minimum values of obtained M s .
Table 3. Maximum and minimum values of obtained M s .
M s d ServoRegulator
minmaxminmax
1.4 1.39231.40881.39041.4216
1.6 1.58361.61301.58191.6183
1.8 1.77251.82561.77381.8266
2.0 1.95182.03591.95272.0356
Table 4. Proportional-integral-derivative parameter, M s , and index values.
Table 4. Proportional-integral-derivative parameter, M s , and index values.
M s d K p T i T d M s J s J r
Servo Optimization1.4 1.0217 1.3331 0.1048 1.3998 0.9576 1.2253
1.6 1.3709 1.4633 0.1090 1.5964 0.7638 1.1531
1.8 1.6359 1.5879 0.1360 1.7937 0.7064 1.0688
2.0 1.8093 1.7116 0.1537 1.9936 0.6970 1.0274
Regulator Optimization1.4 1.0159 0.6876 0.1737 1.4052 1.3048 0.8667
1.6 1.3430 0.6641 0.1681 1.5944 1.0673 0.6466
1.8 1.6065 0.7020 0.1597 1.7913 0.9705 0.5302
2.0 1.8217 0.7174 0.1589 1.9922 0.9458 0.4565
Table 5. Proportional-integral-derivative parameters, M s , and sum of absolute errors (SAE) values using the proposed and internal model control (IMC)-based methods for Equation (40).
Table 5. Proportional-integral-derivative parameters, M s , and sum of absolute errors (SAE) values using the proposed and internal model control (IMC)-based methods for Equation (40).
K p T i T d M s J s J r
IMC-based method [27] λ c = 0.99 0.8248103.104.66961.16251.25001.2493
λ c = 0.98 1.3978104.837.53461.29370.75330.7500
λ c = 0.97 1.8184106.079.46581.56590.59190.5833
Proposed Servo M s d = 1.4 1.91201.12420.06061.40140.59360.5878
M s d = 1.6 2.57241.21070.06901.59190.47330.4705
M s d = 1.8 3.06411.30620.08711.78690.43670.4260
M s d = 2.0 3.38231.42770.09841.98510.43130.4216
Proposed Regulator M s d = 1.4 1.91930.50000.10701.39990.80790.3352
M s d = 1.6 2.53190.46050.10791.58580.76460.2432
M s d = 1.8 3.02340.49280.10381.77780.69510.1975
M s d = 2.0 3.44070.50120.10691.99400.66420.1689
Table 6. Proportional-integral-derivative parameters, M s , and SAE values in the servo design using Equation (44).
Table 6. Proportional-integral-derivative parameters, M s , and SAE values in the servo design using Equation (44).
M s d K p T i T d M s J s J r
Proposed method1.40.93731.04700.14451.40021.17371.1171
1.61.25291.18660.13461.60090.94890.9465
1.81.50081.29760.16491.80230.88330.8638
2.01.65951.37630.18672.00430.87030.8282
Conventional method1.40.97321.08880.13671.42241.14271.1186
1.61.29961.23980.13761.63110.95390.9531
1.81.53321.35850.16821.83630.89550.8846
2.01.68491.43740.18942.04420.88020.8513
Table 7. Proportional-integral-derivative parameters, M s , and SAE values in the regulator design using Equation (44).
Table 7. Proportional-integral-derivative parameters, M s , and SAE values in the regulator design using Equation (44).
M s d K p T i T d M s J s J r
Proposed method1.40.92390.66630.21901.40091.36800.8922
1.61.22540.68370.20201.60171.27680.6949
1.81.46600.72680.18861.80261.20420.5834
2.01.66140.75380.17962.00541.16530.5092
Conventional method1.40.96080.66520.21511.42271.35530.8618
1.61.26480.68650.20021.63221.26850.6745
1.81.48860.71540.18181.83141.21430.5733
2.01.70450.75800.17892.05591.15830.4945
Table 8. Proportional-integral-derivative parameters, M s , and SAE values in the servo design using Equations (46) and (48), respectively.
Table 8. Proportional-integral-derivative parameters, M s , and SAE values in the servo design using Equations (46) and (48), respectively.
M s d K p T i T d M s J s J r
Proposed method1.41.46641.43900.10091.40261.00880.9796
1.61.97251.57880.10661.60100.80200.7978
1.82.35771.71460.13541.80140.74120.7236
2.02.60431.84630.15502.00100.73200.7037
Conventional method [28]1.41.31441.50500.10751.35711.14491.1407
1.61.76191.68410.11941.50240.95530.9495
1.82.07451.82970.14731.65760.88130.8737
2.02.28041.94910.16621.80760.85370.8441
Table 9. Proportional-integral-derivative parameters, M s , and SAE values in the regulator design using Equation (46) and (48), respectively.
Table 9. Proportional-integral-derivative parameters, M s , and SAE values in the regulator design using Equation (46) and (48), respectively.
M s d K p T i T d M s J s J r
Proposed method1.41.43320.72740.17901.40261.29580.6457
1.61.89800.70080.17441.60301.21110.4786
1.82.27240.73520.16451.80281.12440.3915
2.02.57590.75270.16592.00761.07780.3375
Conventional method [28]1.41.29220.75350.19821.35711.33300.7262
1.61.69900.72950.18331.51241.22670.5433
1.82.01910.73280.17311.66201.16310.4495
2.02.29230.78100.16591.81111.08860.3925

Share and Cite

MDPI and ACS Style

Kurokawa, R.; Sato, T.; Vilanova, R.; Konishi, Y. Discrete-Time First-Order Plus Dead-Time Model-Reference Trade-off PID Control Design. Appl. Sci. 2019, 9, 3220. https://doi.org/10.3390/app9163220

AMA Style

Kurokawa R, Sato T, Vilanova R, Konishi Y. Discrete-Time First-Order Plus Dead-Time Model-Reference Trade-off PID Control Design. Applied Sciences. 2019; 9(16):3220. https://doi.org/10.3390/app9163220

Chicago/Turabian Style

Kurokawa, Ryo, Takao Sato, Ramon Vilanova, and Yasuo Konishi. 2019. "Discrete-Time First-Order Plus Dead-Time Model-Reference Trade-off PID Control Design" Applied Sciences 9, no. 16: 3220. https://doi.org/10.3390/app9163220

APA Style

Kurokawa, R., Sato, T., Vilanova, R., & Konishi, Y. (2019). Discrete-Time First-Order Plus Dead-Time Model-Reference Trade-off PID Control Design. Applied Sciences, 9(16), 3220. https://doi.org/10.3390/app9163220

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