PID Control System Analysis & Design
PID Control System Analysis & Design
Control System
System
Analysis and Design
PROBLEMS, REMEDIES, AND FUTURE DIRECTIONS
IMAGESTATE
ith its three-term functionality offering treatment of both transient and steady-state responses,
proportional-integral-derivative (PID) control provides a generic and efficient solution to realworld control problems [1][4]. The wide application of PID control has stimulated and sustained research and development to get the best out of PID [5], and the search is on to find
the next key technology or methodology for PID tuning [6].
This article presents remedies for problems involving the integral and derivative terms. PID design objectives, methods, and future directions are discussed. Subsequently, a computerized, simulation-based approach
is presented, together with illustrative design results for first-order, higher order, and nonlinear plants. Finally,
we discuss differences between academic research and industrial practice, so as to motivate new research
directions in PID control.
FEBRUARY 2006
1066-033X/06/$20.002006IEEE
Authorized licensed use limited to: St. Xavier's Catholic College of Engineering. Downloaded on June 26, 2009 at 00:52 from IEEE Xplore. Restrictions apply.
GPID (s) =
1
U(s)
= KP 1 +
+ TD s ,
E(s)
TI s
where GPD (s) and GPI (s) are the factored PD and PI parts of
the PID controller, respectively, and
(1)
where U(s) is the control signal acting on the error signal E(s),
KP is the proportional gain, TI is the integral time constant, TD
is the derivative time constant, and s is the argument of the
Laplace transform. The control signal can also be expressed in
three terms as
1 4TD /TI
> 0.
2
1
U(s) = KP E(s) + KI E(s) + KD sE(s)
s
= UP (s) + UI (s) + UD (s),
(2)
1
1
= 1 + 2 2 > 1, for all ,
jTI
TI
(5)
1+
1
jTI
= tan1
1
TI
< 0, for all .
(6)
Hence, both gain margin (GM) and phase margin (PM) are
reduced, and the closed-loop system becomes more oscillatory
and potentially unstable.
Automatic Reset
If TI 4TD so that the series form (3) exists, antiwindup can be
achieved implicitly through automatic reset. The factored PI
part of (3) is thus implemented as shown in Figure 1 [8], [9].
Explicit Antiwindup
In nearly all commercial PID software packages and hardware
modules, however, antiwindup is implemented explicitly
through internal negative feedback, reducing UI (s) to [8][10]
Rise Time
Overshoot
Settling Time
Steady-State Error
Stability
Increasing K P
Decrease
Increase
Small Increase
Decrease
Degrade
Increasing K I
Small Decrease
Increase
Increase
Large Decrease
Degrade
Increasing K D
Small Decrease
Decrease
Decrease
Minor Change
Improve
FEBRUARY 2006
Authorized licensed use limited to: St. Xavier's Catholic College of Engineering. Downloaded on June 26, 2009 at 00:52 from IEEE Xplore. Restrictions apply.
UPD(s)
U(s)
KP
Actuator Model
U(s)
1 + jTD = 1 + 2 T2 > 1, for all ,
D
+
1
1 + TIs
to U(s)
is thus 1 + 1/(TI s), the same as the last factor of (3).
1
U(s) U(s)
UI (s) =
KP E(s)
,
TI s
(9)
G(s) =
(11)
(12)
1 + TD2 2
TD
,
KK
min
1,
P
1 + T2 2
T
(13)
G( j)GPD ( j) = KKP
For all , the gain satisfies
(10)
(7)
K Ls
e ,
1 + Ts
KKP
1/2
is
(14)
D tan1 T L.
G( j)GPD ( j) = tan1 T
1
1
TD
1 + jTD = tan1
[0, /2] for all ,
1
(8)
(15)
Hence, if TD > T/KKP and KP > 1/K, then by the Nyquist criterion, the closed-loop system is unstable. This analysis also confirms that some PID mapping formulas, such as the
Ziegler-Nichol (Z-N) formula obtained from the step-response
method, in which KP = (1.2(T/L)) (1/K) and TD is proportional
to L, are valid for only a limited range of values of the T/L ratio.
As an example, consider plant (10) with K = 10, T = 1 s, and
L = 0.1 s [7]. Control by means of a PI controller with KP =
0.644 > 1/K and TI = 1.03 s yields reasonable stability margins
and time-domain performance, as seen in Figures 2 and 3 (Set 1,
red curves). However, when a differentiator is added, gradually
FEBRUARY 2006
Authorized licensed use limited to: St. Xavier's Catholic College of Engineering. Downloaded on June 26, 2009 at 00:52 from IEEE Xplore. Restrictions apply.
TD s
1+
TD
s
(16)
Velocity Feedback
Because a lowpass filter does not completely remove, but
rather averages, impulse derivative signals caused by sudden changes of the setpoint or disturbance, modifications of
the unity negative feedback PID structure are of interest [8].
To block the effect of sudden changes of the setpoint, we
consider a variant of the standard feedback. This variant uses
Nichols Chart
40
30
Gain (dB)
Set 1:
Gain Margin: 7.75757
Phase Margin: 53.37157
Set 1
PIDeasy
Set 2
Set 3
PIDeasy:
Gain Margin: 9.16494
Phase Margin: 64.72558
Set 3: TD = 0.2
20
Set 2:
Gain Margin: 3.41323
Phase Margin: 86.39067
Set 2: TD = 0.1
10
Set 3:
Gain Margin: 2.26496
Phase Margin:
0
10
20
PIDeasy: TD = 0.0303
3.5
Set 1: TD = 0
3.0
2.5
2.0
1.5
Phase (deg.)
1.0
0.5
0.0
102
FIGURE 2 Destabilizing effect of the derivative term, measured in the frequency domain by GM and PM. Adding a derivative term increases
both the GM and PM, although raising the derivative gain further tends to reverse the GM and destabilize the closed-loop system. For
example, if the derivative gain is increased to 20% of the proportional gain (TD = 0.2 s), the overall open-loop gain becomes greater than
2.2 dB for all . At = 30 rad/s, the phase decreases to while the gain remains above 2.2 dB. Hence, by the Nyquist criterion, the
closed-loop system is unstable. It is interesting to note that MATLAB does not compute the frequency response as shown here, since MATLAB handles the transport delay factor ejL in state space through a Pad approximation.
FEBRUARY 2006
Authorized licensed use limited to: St. Xavier's Catholic College of Engineering. Downloaded on June 26, 2009 at 00:52 from IEEE Xplore. Restrictions apply.
the process variable instead of the error signal for the derivative action [14], as in
u(t) = KP e(t) + KI
e( ) d KD
d
y(t),
dt
(17)
where y(t) is the process variable, e(t) = r(t) y(t) is the error
signal, and r(t) is the setpoint or reference signal. The last term
of (17) forms velocity feedback and, hence, an extra loop that
is not directly affected by a sudden change in the setpoint.
However, sudden changes in disturbance or noise at the plant
output can cause the differentiator to produce a theoretically
infinite control signal.
Setpoint Filter
To further reduce sensitivity to setpoint changes and avoid
overshoot, a setpoint filter may be adopted. To calculate the
proportional action, the setpoint signal is weighted by a factor
b < 1, as in [8] and [14]
t
d
u(t) = KP b r(t) y(t) + KI 0 e( ) d KD y(t).
dt
Prefilter
For setpoint tracking applications, an alternative to using a Type
B or C structure is to cascade the setpoint with a prefilter that
has critically damped dynamics. When a step change in the setpoint occurs, continuous output of the prefilter helps achieve
soft start and bumpless control [8], [19]. However, a prefilter
does not solve the problem caused by sudden changes in the
disturbance since it is not embedded in the feedback loop.
(18)
Step Response
Set 1:
Kp: 0.6439
Ti:
1.0278
Td: 0
ITAE: 228.93
1.4
Set 1: TD = 0
1.2
PIDeasy:
Set 3
Kp: 0.6439
Ti:
1.02781
Td: 0.03025
ITAE: 104.86
Output
1.0
0.8
Set 2:
0.6
Kp: 0.6439
Ti:
1.0278
Td: 0.1
ITAE: 247.69
Set 2: TD = 0.1
0.4
0.2
Set 3:
0.0
0.0
0.2
0.4
PIDeasy: TD = 0.0303
0.6
0.8
Time (Sec.)
1.0
1.2
1.4
Kp: 0.6439
Ti:
1.0278
Td: 0.2
ITAE: 52,547.83
Set 3: TD = 0.2
FIGURE 3 Destabilizing effect of the derivative term, confirmed in the time domain by the closed-loop step response. Although increasing the
derivative gain initially decreases the oscillation, this trend soon reverses and the oscillation grows into instability.
FEBRUARY 2006
Authorized licensed use limited to: St. Xavier's Catholic College of Engineering. Downloaded on June 26, 2009 at 00:52 from IEEE Xplore. Restrictions apply.
Heuristic Methods
Heuristic methods evolve from empirical tuning (such as the ZN tuning rule), often with a tradeoff among design objectives.
Heuristic search now involves expert systems, fuzzy logic,
neural networks, and evolutionary computation [19], [22].
Mode
Action
PB
PI
PB
Value
1.084
L
T
0.977
L
1.164K
T
0.68
T
L
40.44 T
0.947
L
0.7369K
T
0.738
T
L
51.02 T
0.995
T
L
157.5 T
0.947
L
0.5438K
T
0.995
T
L
157.5 T
2.04K
TI
PB
PID
TI
TD
PD
PB
Analytical Methods
Because of the simplicity of PID control, parameters can be
derived analytically using algebraic relations between a plant
TD
FEBRUARY 2006
Authorized licensed use limited to: St. Xavier's Catholic College of Engineering. Downloaded on June 26, 2009 at 00:52 from IEEE Xplore. Restrictions apply.
Phase Margin ()
15
10
5
103
102
101
100
101
102
103
102
101
100
101
102
103
80
70
60
50
103
L/T
FIGURE 5 Gain and phase margins resulting from PIDeasy designs for firstorder delayed plants with various L/T ratios. While requirements of fast
transient response, no overshoot, and zero steady-state error are accommodated by time-domain criteria, multiobjective design goals provide
frequency-domain margins in the range of 911 dB and 6566 .
PIDeasy [7] is a software package that uses automatic simulations to search globally for controllers that meet all five
design objectives in both the time and frequency domains.
The search is initially performed offline in a batch mode
[19] using artificial evolution techniques that evolve both
TABLE 3 Multioptimal PID settings for normalized typical high-order plants. Since PIDeasys search priorities are time-domain
tracking and regulation, the corresponding gain and phase margins are given to assess frequency-domain properties.
Kp
Plants
1
G1 (s) =
(s + 1)
G2 (s) =
G3 (s) =
1
(s + 1)(1 + s)(1 + 2 s)(1 + 3 s)
1 s
(s + 1)3
1
es
G4 (s) =
(1 + s)2
PID Coefficients
Ti (s)
Td (s)
Resultant Margins
GM (dB)
PM ( )
=1
=2
=3
=4
=8
92.1
1.95
1.12
0.83
0.50
1.0
1.61
2.13
2.61
4.31
0.0022
0.14
0.28
0.43
1.01
26.8
13.9
9.05
102
62.4
60.7
61
58.9
= 0.1
= 0.2
= 0.5
5.53
2.87
1.19
1.03
1.08
1.36
0.04
0.07
0.17
52.8
38.6
19.1
68.7
66.3
62.6
= 0.1
= 0.2
= 0.5
= 1.0
= 2.0
= 5.0
1.03
0.96
0.79
0.63
0.48
0.36
2.15
2.18
2.23
2.30
2.39
2.58
0.31
0.33
0.39
0.47
0.57
0.72
19.4
16.6
13
7.52
7.45
2.69
61.2
61.6
62.4
50.9
58.6
40.4
= 0.1
= 0.2
= 0.5
= 2.0
= 5.0
= 10
0.23
0.30
0.49
1.04
1.42
1.65
0.43
0.59
1.07
3.49
8.32
16.35
0.12
0.17
0.26
0.49
0.92
1.59
10.4
10.4
10.5
15
24.2
32.8
66
65.8
65.6
62.4
62.1
62.1
FEBRUARY 2006
Authorized licensed use limited to: St. Xavier's Catholic College of Engineering. Downloaded on June 26, 2009 at 00:52 from IEEE Xplore. Restrictions apply.
controller parameters and their associated structures. For practical simplicity and reliability, the standard PID structure is maintained as much as possible, while allowing augmentation with
either lowpass or median filtering for the differentiator and with
explicit antiwindup for the integrator. The resulting designs are
then embedded in the PIDeasy package. Further specific tuning
can be continued by local, fast numerical optimization if the
plant differs from its model or data used in the initial design.
0.9
Output (mol/l )
0.8 0.74
0.7
0.6
dmax
0.4
0.3
0.31
1
0.1
1
,
(s + 1)
0.49
0.5
0.2
= 1, 2, 3, 4, 8,
(19)
1
,
(s + 1)(1 + s)(1 + 2 s)(1 + 3 s)
= 0.1, 0.2, 0.5,
(20)
1 s
,
= 0.1, 0.2, 0.5, 1, 2, 5,
(21)
G3 (s) =
(s + 1)3
1
es,
= 0.1, 0.2, 0.5, 2, 5, 10. (22)
G4 (s) =
(1 + s)2
G2 (s) =
3
Input (l/h)
where
y(t) = concentration in the outlet stream (mol/),
u(t) = flow rate of the feed stream (/h),
K = rate of reaction (/mol-h),
V = reactor volume ( ),
d = concentration in the inlet stream (mol/ ).
The setpoint, equilibrium, or steady-state operating trajectory
of the plant is governed by
Ky2 +
1
d
u y u = 0.
V
V
(24)
S3
0.8
(23)
Weighting
dy(t)
1
= Ky2 (t) +
d y(t) u(t),
dt
V
S2
S1
0.6
0.4
0.2
0
0.49
0.31
0.1
0.74
0.8 0.9
FEBRUARY 2006
Authorized licensed use limited to: St. Xavier's Catholic College of Engineering. Downloaded on June 26, 2009 at 00:52 from IEEE Xplore. Restrictions apply.
0.8
u(t) = [S1
0.6 0.53
0.4
0.2
0
0.5
1.5
2.5 3
Time (h)
3.5
4.5
0.5
1.5
2.5 3
Time (h)
3.5
4.5
6
4
2
0
S2
9.82 1.22 0.0376
1
S3 ] 15.6 0.784 0.0241
p1 e(t),
p
28.6 0.481 0.0137
(25)
0.6
0.5
0.3
0.9
4
5
Time (h)
4
5
Time (h)
10
8
6
4
2
0
FIGURE 9 Performance of the pseudolinear PID network applied to the nonlinear chemical process (23) at multiple operating levels that are not originally
used in the design process. The network tracks these setpoint changes accurately without oscillation. It can be seen that the control effort increases disproportionally to the setpoint change along the nonlinear trajectory, compensating
for the decreasing gain of the plant when the operating level is raised.
FEBRUARY 2006
Authorized licensed use limited to: St. Xavier's Catholic College of Engineering. Downloaded on June 26, 2009 at 00:52 from IEEE Xplore. Restrictions apply.
in the feedback loop, a median filter also outperforms a prefilter in dealing with disturbances.
Over the past half century, researchers have sought the
next key technology for PID tuning and modular realization.
Many design methods can be computerized and, with simulation packages widely used, the trend of computerizing simulation-based designs is gaining momentum. Computerizing
enables simulations to be carried out automatically, which
facilitates the search for the best possible PID settings for the
application at hand. A simulation-based approach requires no
artificial minimization of the control amplitude and helps
improve sluggish transient response without windup.
In tackling PID problems, it is desirable to use standard
PID structures for a reasonable range of plant types and operations. Modularization around standard PID structures
should also help improve the cost effectiveness of PID control
and maintenance. This way, robustly optimal design methods
such as PIDeasy can be developed. By including system identification techniques, the entire PID design and tuning process
can be automated, and modular code blocks can be made
available for timely application and real-time adaptation.
ACKNOWLEDGMENTS
This article is based on [25]. Kiam Heong Ang and Gregory
Chong are grateful to the University of Glasgow for a postgraduate research scholarship and to Universities UK for an
Overseas Research Students Award. The authors thank Prof.
Hiroshi Kashiwagi of Kumamoto University and Mitsubishi
Chemical Corp., Japan, for the nonlinear reaction process
model and data.
AUTHOR INFORMATION
Yun Li (Y.Li@elec.gla.ac.uk) is a senior lecturer at the University of Glasgow, United Kingdom, where he has taught and
conducted research in evolutionary computation and control
engineering since 1991. He worked in the U.K. National Engineering Laboratory and Industrial Systems and Control Ltd,
Glasgow, in 1989 and 1990. In 1998, he established the IEEE
CACSD Evolutionary Computation Working Group and the
European Network of Excellence in Evolutionary Computing
(EvoNet) Workgroup on Systems, Control, and Drives. He
was a visiting professor at Kumamoto University, Japan. He is
currently a visiting professor at the University of Electronic
Science and Technology of China. His research interests are in
parallel processing, design automation, and discovery of engineering systems using evolutionary learning and intelligent
search techniques. He has advised 12 Ph.D. students and has
140 publications. He can be contacted at the Department of
Electronics and Electrical Engineering, University of Glasgow,
Glasgow G12 8LT, U.K.
Kiam Heong Ang received a First-Class Honors B.Eng. and
a Ph.D. degree in electronics and electrical engineering from
the University of Glasgow, United Kingdom, in 1996 and
2005, respectively. From 19972000, he was a software engineer with Advanced Process Control Group, Yokogawa Engi-
neering Asia Pte. Ltd., Singapore. Since 2005, within the same
company, he has been working on process industry standardization and new technology development. His research
interests include evolutionary, multiobjective learning, computational intelligence, control systems, and engineering
design optimization.
Gregory Chong received a First-Class Honors B.Eng. degree
in electronics and electrical engineering from the University of
Glasgow, United Kingdom, in 1999. He is completing his
Ph.D. at the same university in evolutionary, multiobjective
modeling, and control for nonlinear systems.
REFERENCES
[1] J.G. Ziegler and N.B. Nichols, Optimum settings for automatic controllers, Trans. ASME, vol. 64, no. 8, pp. 759768, 1942.
[2] W.S. Levine, Ed., The Control Handbook. Piscataway, NJ: CRC Press/IEEE
Press, 1996.
[3] L. Wang, T.J.D. Barnes, and W.R. Cluett, New frequency-domain design
method for PID controllers, Proc. Inst. Elec. Eng., pt. D, vol. 142, no. 4, pp.
265271, 1995.
[4] J. Quevedo and T. Escobet, Eds., Digital control: Past, present and future
of PID control, in Proc. IFAC Workshop, Terrassa, Spain, Apr. 5, 2000.
[5] I.E.E. Digest, Getting the best out of PID in machine control, in Digest
IEE PG16 Colloquium (96/287), London, UK, Oct. 24, 1996.
[6] P. Marsh, Turn on, tune inWhere can the PID controller go next, New
Electron., vol. 31, no. 4, pp. 3132, 1998.
[7] Y. Li, W. Feng, K.C. Tan, X.K. Zhu, X. Guan, and K.H. Ang, PIDeasy and
automated generation of optimal PID controllers, in Proc. 3rd Asia-Pacific
Conf. Control and Measurement, Dunhuang, P.R. China, 1998, pp. 2933.
[8] K.J. strm and T. Hgglund, PID Controllers: Theory, Design, and Tuning.
Research Triangle Park, NC: Instrum. Soc. Amer., 1995.
[9] F.G. Shinskey, Feedback Controllers for the Process Industries. New York:
McGraw-Hill, 1994.
[10] C. Bohn and D.P. Atherton, An analysis package comparing PID antiwindup strategies, IEEE Contr. Syst. Mag., vol. 15, no. 2, pp. 3440, Apr.
1995.
[11] K.J. strm and T. Hgglund, The future of PID control, Contr. Eng.
Pract., vol. 9, no. 11, pp. 11631175, 2001.
[12] Techmation Inc., Techmation [Online], May 2004. Available: http://protuner.com
[13] J.P. Gerry and F.G. Shinskey, PID controller specification, white paper
[Online]. May 2004. Available: http://www.expertune.com/PIDspec.htm
[14] BESTune, PID controller tuning [Online]. May 2004. Available:
http://bestune.50megs.com
[15] Honeywell International Inc. [Online]. May 2004. Available: http://
www.Acs.Honeywell.Com/Ichome/Rooms/DisplayPages/LayoutInitial
[16] Y. Li, K.H. Ang, and G. Chong, Patents, software, and hardware for
PID control, IEEE Contr. Syst. Mag., vol. 26, no. 1, pp. 4254, 2006.
[17] J.P. Gerry, A comparison of PID control algorithms, Contr. Eng., vol.
34, no. 3, pp. 102105, Mar. 1987.
[18] A. Kaya and T.J. Scheib, Tuning of PID controls of different structures,
Contr. Eng., vol. 35, no. 7, pp. 6265, July 1988.
[19] W. Feng and Y. Li, Performance indices in evolutionary CACSD
automation with application to batch PID generation, in Proc. 10th IEEE Int.
Symp. Computer Aided Control System, Hawaii, Aug. 1999, pp. 486491.
[20] R. Gorez, A survey of PID auto-tuning methods, Journal A, vol. 38, no. 1,
pp. 310, 1997.
[21] A. ODwyer, Handbook of PI and PID Controller Tuning Rules. London:
Imperial College Press, 2003.
[22] Y. Li, K.H. Ang, G. Chong, W. Feng, K.C. Tan, and H. Kashiwagi,
CAutoCSDEvolutionary search and optimisation enabled computerautomated control system design, Int. J. Automat. Comput., vol. 1, no. 1,
pp. 7688, 2004.
[23] Specification Data File of Commander 355, ABB, SS/C355, Issue 3, 2001.
[24] K.J. strm and T. Hgglund, Benchmark systems for PID control, in
Proc. IFAC Workshop, Terrassa, Spain, 2000, pp. 165166.
[25] K.H. Ang, G. Chong, and Y. Li, PID control system analysis, design, and
technology, IEEE Trans. Contr. Syst. Tech., vol. 13, no. 4, pp. 559576, 2005.
FEBRUARY 2006
Authorized licensed use limited to: St. Xavier's Catholic College of Engineering. Downloaded on June 26, 2009 at 00:52 from IEEE Xplore. Restrictions apply.