In this paper, aspects of analytical design of PID controllers are studied, by combining pole placement technique with symmetrical optimum criterion. The proposed method is based on low-order plant model with pure integrator, and it can be used for both fast and slow processes. Starting from the desired closed-loop transfer function, which contains a second-order oscillating system and a lead-lag compensator, it is shown that the zero value depends on the real-pole value of closed-loop transfer function. In addition, there is only one pole value, which satisfies the assumptions of symmetrical optimum criterion imposed to open-loop transfer function. In these conditions, by combining the pole placement technique with symmetrical optimum criterion, the analytical expressions of the controller parameters can be simplified. For simulations, PID autopilot design for heading control problem of a conventional ship is considered.
1. Introduction
PID controllers are widespread in industry, being the most popular ones for decades. Hence, even small improvements in PID design could have major impact worldwide [1]. There are many theoretical and practical papers on PID tuning algorithms, in time and frequency domain. They are using analytical [1–5], graphical [6–8], or empirical methods, including artificial intelligence [9, 10].
Analytical methods rely on low-order plant models characterized by small number of parameters. Two reduced plant models are mostly used: pure integrator plus first-order model, like those for thermal and electromechanical processes, and the first-order plus dead-time model, like in chemical processes, respectively [11].
The proposed analytical method in the paper is based on low-order plant model with pure integrator, and it can be used for both fast and slow processes. An application example is discussed for autopilot design of conventional ships in course-keeping or course-changing control problems.
Pole placement technique (PPT) [12–14] and symmetrical optimum criterion (SOC) [15–18] were intensively studied as separated methods on both types of reduced plant models mentioned above.
If PPT is used to synthesize the PID controller, the first step is to specify some performance conditions of the closed-loop system, which lead to the analytical expression of the closed-loop transfer function (CLTF). The PID parameters depend on the CLTF and plant model parameters. In general, the closed-loop system is approximated by a second order reference model. If the plant model has pure integrator element and PID type is desired for the controller, then the CLTF must contain an additional lead-lag compensator. In this case, the reference model has four parameters, two from the second order reference model and the other two from the compensator. All these parameters are present in the analytical expressions of PID controller. In addition, if the well known symmetrical optimum criterion [6] is used, the parameters of compensator have particular values, as stated in the paper, which simplify the PID parameters.
By combining the pole placement technique with symmetrical optimum criterion, the analytical expressions of the PID parameters depend only on the parameters of the plant and second order reference models.
In this paper, PID autopilot design for heading control problem of a conventional ship is considered in simulations. The ship model is obtained from the equations describing the horizontal motion of the ship, expressing conservation of hydrodynamic forces and moments. For ship dynamics modeling, which characterize a slow process, the first type of reduced model is used, with pure integrator plus first-order model.
The autopilot must assure performance conditions in both course-keeping and course-changing control problems. It must be PID type to assure null stationary error for both step and ramp variations on reference and disturbances inputs [19, 20]. Thus, in course-keeping control, the PID autopilot can compensate the stationary error (e.g., drift effect), generated by constant or slowly time-variable disturbances (e.g., wind, sea currents). Also, in course-changing control problems, PID autopilot can act as a tracking system to a desired trajectory, consisting of line segments [21].
As a result, the reference model contains a second-order oscillating system and a lead-lag compensator, with real negative values for pole-zero pair, and the open-loop transfer function (OLTF) contains a double-integrator element.
The goal of this paper is to find the simplified analytical expressions of PID controller parameters, which satisfy two simultaneous conditions: the desired behavior of closed-loop system and symmetrical characteristics of the open-loop transfer function. Applying PPT, the controller parameters depend on the parameters of the plant and closed-loop transfer functions. The resulting open-loop transfer function contains a double-integrator element and a pole-zero pair, with real negative values. Hence, the symmetrical optimum criterion can be used. Imposing symmetrical characteristics of the open-loop transfer function, which can be obtained in certain conditions stated in the paper, the analytical expressions of the controller parameters can be simplified.
The paper is organized as follows. Section 2 provides mathematical models for the plant, PID controller, and reference closed-loop system used in simulations, for both fast and slow processes. In Section 3, the PID controller synthesis using PPT is presented. The analytical expressions of controller parameters are simplified in Section 4, by imposing symmetrical characteristics of OLTF. Section 5 describes the simulation results, using a PID autopilot for heading control of a ship. Conclusions are presented in Section 6.
2. Mathematical Models
Mathematical models for plant, PID controller, and reference closed-loop system used in simulations are discussed in this section. The classical factorized forms of transfer functions are considered, to express poles and zeros, based on which pole placement technique and symmetrical optimum criterion can be used.
In the paper, a low-order plant model with pure integrator is considered, for both fast and slow processes. The controller is of PID type, being imposed by performance conditions of closed-loop system, which are detailed at the end of this section.
The classical structure of the control loop, for a linear SISO process without disturbances, is illustrated in Figure 1.
Classical structure of the control loop.
The plant has a low-order model, which contains a pure integrator. It is characterized by a dominant time constant (TP) and a gain coefficient (kP). The analytical expressions of plant model and PID controller depend on the process type, fast or slow [22].
(a) If the process is fast, then the small time constants can not be neglected. An equivalent small time constant (TΣ) is added to the model, corresponding to the sum of parasitic time constants:
(1)HP(s)=kPs·(sTP+1)·(sTΣ+1),
where TΣ<TP.
In this case, the PID controller is of the following form:
(2)HC(s)=kCsTC·(sTC+1)·(sTC′+1),
where TΣ<TC′<TC.
The open-loop transfer function, H(s)=HC(s)·HP(s), is as follows:
(3)H(s)=kCsTC·(sTC+1)·(sTC′+1)·kPs·(sTP+1)·(sTΣ+1).
Using pole cancellation, the nonzero dominant pole of the plant model is cancelled by choosing TC′=TP. Thus, only two controller parameters must be determined: kC and TC.
It can be observed that if the process does not have any non-zero dominant pole (model without time constant TP), then the controller is of PI type (without time constant TC′), and the same parameters must be determined (kC and TC).
(b) If the process is slow, the equivalent small time constant can be neglected (TΣ≪TP) and the plant model is as follows:
(4)HP(s)=kPs·(sTP+1).
In this case, the PID controller contains a supplementary degree of freedom and it is of the following form:
(5)HC(s)=kCsTC·(sTC+1)·sTC′+1sT1+1,
where T1<TC′<TC.
The open-loop transfer function for slow process is as follows:
(6)H(s)=kCsTC·(sTC+1)·sTC′+1sT1+1·kPs·(sTP+1).
Again, the non-zero dominant pole of the plant model is cancelled, using pole cancellation, resulting TC′=TP. In this case, three controller parameters must be determined: kC, TC, and T1. This is a more general case, due to the time constant T1, which is not imposed by the process and it can be chosen. Therefore, in this paper, the plant model given in (4) is used for simulations, but discussions include both plant models given in (1) and (4).
The open-loop transfer functions, given in (3) and (6), have similar expressions, after pole cancellation:
(7)H(s)=HC(s)·HP(s)=kC·kP·(sTC+1)s2TC·(sT+1),
where the time constant T has different meanings, depending on the type of process considered. In the first case for fast processes, it represents the equivalent small time constant imposed by the process (T=TΣ), while in the second case it is a controller parameter (T=T1).
In general, the behavior of the closed-loop system is approximated by a second order reference model:
(8)H0(s)=ω02s2+2ζω0s+ω02,
where ω0>0 is the natural frequency and ζ>0 is the damping coefficient. The performance conditions of the closed-loop system can be specified, by imposing the model parameters (ζ and ω0).
Starting from the plant model with pure integrator and considering the reference model from (8), it results a real-PD controller. The corresponding OLTF is as follows:
(9)H(s)=HC(s)·HP(s)=H0(s)1-H0(s)=ω02s(s+2ζω0).
It can be observed that the open-loop transfer function, corresponding to CLTF presented in (8), has one pure integrator, which assures null stationary error only for step variations of reference input, denoted r in Figure 1.
If the controller must assure null stationary error for both step and ramp variations on reference and disturbances inputs, then the controller must be PID type. In this case, the open-loop transfer function contains a double-integrator element, one being included from the plant model, and the other one from the PID controller. The resulted OLTF has the general form expressed in (7), which cannot be generated from the reference model given in (8).
To obtain the expression of OLTF given in (7), the reference model must be completed with a lead-lag compensator, which contains a pole-zero pair, with real and negative values [22]:
(10)H0(s)=k0·(s+z)(s2+2ζω0s+ω02)·(s+p),
where z>0, p>0, and k0=ω02·p/z, so that |H0(0)|=1.
Particular values for pole-zero pair of lead-lag compensator simplify the analytical expressions of the controller parameters, as stated in the next sections. In addition, symmetrical characteristics of OLTF are obtained.
3. Aspects of PID Controller Design Using Pole Placement Technique
In this section, existence conditions of parameters of CLTF are discussed, and PID parameters are determined based on PPT [19]. Consider the control structure shown in Figure 1, with the OLTF given in (7) and the desired CLTF in (10).
Proposition 1.
For every value of the pole (s=-p) in the desired closed-loop transfer function given in (10), there is only one zero value (s=-z) for which the open-loop transfer function has a double-pole in origin, and in addition, the zero frequency is smaller than the pole frequency: z<p.
Proof.
The proposition demonstration includes the next lemma results.
Lemma 2.
The necessary and sufficient condition, for the existence of a double-pole in origin for the open-loop transfer function of a control system illustrated in Figure 1, starting from a desired CLTF of the form given in (10), is as follows:
(11)2ζpz+ω0z-ω0p=0.
Proof.
The transfer function of the open-loop system can be computed starting from the desired CLTF:
(12)H(s)=HC(s)·HP(s)=H0(s)1-H0(s)=ω02·p/z·(s+z)s3+s2(p+2ζω0)+s[2ζω0p+ω02(1-p/z)].
From (12) it can be observed that a double-pole in origin is obtained for every frequency, if and only if the coefficient of the third term of denominator is null:
(13)2ζω0p+ω02(1-pz)=0,
which is equivalent with equation in (11) for ω0≠0.
As a result, the unique zero value results, whose expression depends on the selected pole value:
(14)s=-z=-ω0p2ζp+ω0,∀p,ω0,ζ>0.
For every frequency (p) of the pole, the corresponding frequency (z) of the zero given in (14) is smaller:
(15)z=ω0p2ζp+ω0<p,∀p,ω0,ζ>0.
In this case, the real values of pole-zero pair and conjugate complex poles, of the desired CLTF given in (10), are illustrated in the complex plane in Figure 2.
Poles and zero of desired CLTF given in (10).
Using (13) and (15) in (12), the expression of the open-loop transfer function is obtained:
(16)H(s)=ω0·(2ζp+ω0)·(s+(ω0p/(2ζp+ω0)))s2·[s+(p+2ζω0)].
The zero and pole frequencies of the open-loop transfer function are denoted by ωz and ωp, respectively:
(17)ωz=ω0p2ζp+ω0=z,ωp=2ζω0+p.
The open-loop transfer function can be rewritten as follows:
(18)H(s)=ω0(2ζp+ω0)·(s+ωz)s2·[s+ωp].
Putting into evidence the time constants, the open-loop transfer function can be rewritten, like in (7):
(19)H(s)=ω02p·(s((2ζp+ω0)/ω0p)+1)s2(2ζω0+p)(s(1/(2ζω0+p))+1).
Equaling OLTF without pole cancellation, from (3) or (7), with (19), an equality of two polynomials of 3rd order in s variable is obtained:
(20)kP·kC·(2ζω0+p)·(sTC+1)·(sTC′+1)·(s12ζω0+p+1)=TC·ω02p·(sT+1)·(sTP+1)·(s2ζp+ω0ω0p+1),
where T=TΣ or T=T1, depending on the plant model.
Because the equality must be true for every frequency, it results a four equation system:(21a)kP·kC·(2ζω0+p)=TC·ω02p,(21b)TC·TC′·12ζω0+p=T·TP·2ζp+ω0ω0p,(21c)TC·TC′+(TC+TC′)·12ζω0+p=T·TP+(T+TP)·2ζp+ω0ω0p,(21d)TC+TC′+12ζω0+p=T+TP+2ζp+ω0ω0p.
The solutions of the equation system are the PID controller parameters [19]:(22a)kC=ω0kP·2ζp+ω02ζω0+p,(22b)TC=2ζp+ω0ω0p=1ωz,(22c)TC′=TP,(22d)T=12ζω0+p=1ωp.
It can be observed that the solution (22c) represents the pole cancellation condition, and it does not depend on the p pole value.
If a fast process is considered, then the time constant T given in (22d) represents the equivalent small time constant of the plant model (T=TΣ). If the process is slow, then the same time constant T is a controller parameter (T=T1).
The time constants must satisfy the inequalities:
(23)T<TC′<TC.
This can be transposed into frequency domain:
(24)ωz<1TP<ωp.
Using (17) in (24), a system of 2 inequalities is obtained:(25a)2ζω0+p>1TP,(25b)2ζp+ω0>pω0TP.
The reference model must be chosen so that the parameters of closed-loop transfer function (ω0, ζ, and p) satisfy this system of inequalities.
4. Design Aspects by Imposing Symmetrical Characteristics of OLTF
After applying PPT, the parameters of PID controller depend on the CLTF parameters (ω0, ζ, and p). In general, ω0 and ζ characterize the desired behavior of the closed-loop system, and they have fixed values, while the pole value can be chosen. Specific pole values can be imposed by using supplementary conditions. In this paper, the conditions for choosing the pole value refer to the symmetrical optimum criterion, which simplify the expressions of PID parameters.
The goal is to find that pole value of the CLTF, which satisfies the assumptions of symmetrical optimum criterion around natural frequency ω0, for the transfer function of open-loop system given in (16). Using this value, the expressions of PID parameters given in (22a), (22b), (22c), and (22d) are simplified.
Proposition 3.
There is only one admissible value for the pole (s=-p) of CLTF given in (10), so that the corresponding OLTF given in (16) has symmetrical characteristics around ω0, and this value is as follows:
(26)p=ω0,∀p,ω0,ζ>0.
Proof.
For the specified open-loop transfer function, only the symmetry of magnitude-frequency characteristic around natural frequency ω0 must be imposed, as it implies also the symmetry of phase-frequency characteristic.
The general form of the symmetrical optimum criterion imposes two conditions for magnitude-frequency characteristic of OLTF:
(a) cross-over frequency ω0 must be equally placed between zero and pole frequencies on the 10-base logarithmic scale:
(27)ω0ωz=ωpω0;
(b) for cross-over frequency ω0, the magnitude-frequency characteristic of OLTF must have 0 dB:
(28)|H(jω0)|=1.
Using (17) in (27), the first condition in p variable is as follows:
(29)ω0p2ζp+ω0·(2ζω0+p)=ω02⟺p2=ω02.
From (29), the solution results are as follows:
(30)p=ω0,∀p,ω0,ζ>0.
Therefore, the condition given in (27) is satisfied if and only if p=ω0.
For the second condition in (28), the magnitude of open-loop transfer function in frequency ω0 is computed from (18):
(31)|H(jω0)|=(2ζp+ω0)·ωz2+ω02ω0·ωp2+ω02.
The frequencies ωz and ωP are replaced with their expressions from (17), resulting the following:
(32)|H(jω0)|=p2+ω02+4ζpω0+4ζ2p2p2+ω02+4ζpω0+4ζ2ω02.
Using (32) in (28), the same solution in (30) is obtained: p2=ω02⇔p=ω0,∀p,ω0,ζ>0.
Using (30) in (15), it results the following:
(33)p=ω0⟹z=ω02ζ+1.
The expression of CLTF becomes
(34)H0(s)=ω02·(2ζ+1)·(s+(ω0/(2ζ+1)))(s2+2ζω0s+ω02)·(s+ω0).
From (17), the zero and pole frequencies of OLTF are
(35)ωz=ω02ζ+1,ωp=ω0·(2ζ+1).
The open-loop transfer function can be rewritten as follows:
(36)H(s)=ω02(2ζ+1)·(s+(ω0/2ζ+1))s2·[s+ω0(2ζ+1)].
The real values of pole-zero pair and complex conjugate poles, of the CLTF given in (34), are illustrated in Figure 3.
Poles and zero of CLTF for p=ω0.
The position of the zero, s=-z, depends on ζ parameter. Three situations are possible:
if ζ∈(0,1/2), then -ω0<-z<-ζω0 and the zero is placed between the two points: s=-ω0 and s=-ζω0, respectively;
if ζ=1/2, then -z=-ζω0. This is the particular case of the Kessler’s symmetrical optimum method [6];
if ζ>1/2, then 0>-z>-ζω0 and the zero is placed to the right of the point s=-ζω0. This is the case illustrated in Figure 3.
Knowing the pole value of CLTF, p=ω0, the simplified parameters of PID controller result from (22a), (22b), (22c), and (22d): (37a)kC=ω0kP,(37b)TC=2ζ+1ω0,(37c)TC′=TP,(37d)T=1(2ζ+1)·ω0.
In this case, the system of inequalities from (25a) and (25b) is simpler, including only 2 variables from reference model (ω0 and ζ):(38a)(2ζ+1)·ω0>1TP,(38b)2ζ+1ω0>TP.
Proposition 4.
In the case of symmetrical characteristics of the OLTF given in (36) around the natural frequency ω0, the phase margin and the distance between the frequency points on the 10-base logarithmic scale depend only on the parameter ζ.
Proof.
The distance between the frequency points on the logarithmic scale can be easily obtained, using (35) in (27):
(39)ω0ωz=ωpω0=2ζ+1.
Using the OLTF given in (36), the phase margin results as follows:
(40)φm=π+arg(H(jω0))=arctg(2ζ+1)-arctg(12ζ+1).
For particular case of the Kessler’s symmetrical optimum method (ζ=0.5), the distance between frequency points is equal with an octave and the phase margin is φm=36.87[deg].
5. Simulation Results
For simulations, the heading control problem of a ship is considered, using a PID autopilot.
The ship yaw angle (ψ) is the model output, and rudder command (δ) generated by autopilot is the control input. The ship model is linear with parametric uncertainties, depending on ship loading conditions and forward speed of the ship, while the wave characteristics change frequently.
It is a first order Nomoto model of the form given in (4), being identified for a ship speed of 22 knots [21]:
(41)HP(s)=ψ(s)δ(s)=kPs·(sTP+1),
where ψ(s) and δ(s) represent the Laplace transforms of yaw angle and rudder angle, respectively.
The ship model parameters are as follows:
(42)kP=-0.0834[s-1],TP=5.98[s].
The autopilot model is given in (5) and the desired CLTF is given in (10). The parameters ω0 and ζ are chosen from performance conditions [20]:
(43)ζ=0.9,ω0=0.1[rad/s].
Using the values from (42) and (43), it can be easily verified that the conditions in (38a) and (38b) are satisfied.
From the desired closed-loop transfer function with pole-zero pair specified in (30) and (33), and imposing symmetrical characteristics of the OLTF, the expressions in (34) and (36) are obtained. From (37a), (37b), (37c), and (37d), the autopilot parameters are obtained:
(44)kC=-1.2,TC=28[s],TC′=5.98[s],T1=3.57[s].
The step response of the closed-loop transfer function is illustrated in Figure 4. The step reference input is represented with dotted line.
Step response of the closed-loop system.
Considering the second order reference model given in (8), the resulting autopilot is of real-PD type. The response is over-dumped, being represented with dashed line. For PID autopilot based on reference model with lead-lag compensator given in (10), the response is overshot and it is drawn with continuous line.
Similarly, the ramp response of the CLTF is illustrated in Figure 5. The ramp reference input is represented with dotted line. Considering the second order reference model given in (8) with real-PD type, the response is shown with dashed line. An important stationary error can be observed. For reference model with lead-lag compensator given in (10) with PID autopilot, the response is represented with continuous line. In this case, a null stationary error is obtained.
Ramp response of the closed-loop system.
PID autopilot assures null stationary error for both step and ramp variations on reference inputs.
For reference model with compensator and PID autopilot, the symmetrical characteristics of the open-loop transfer function are illustrated in Figure 6. It can be observed that the two conditions for magnitude-frequency characteristic of OLTF generate the symmetry of phase-frequency characteristic. In this case, the phase margin is φm=50.69[deg].
Symmetrical characteristics of OLTF.
6. Conclusions
By combining the pole placement technique with symmetrical optimum criterion, the analytical expressions of PID parameters are simplified. In addition, it was demonstrated that there is only one possible pair for the pole-zero values of closed-loop transfer function, so that the PID controller satisfied two simultaneous conditions: the desired behavior of closed-loop system and symmetrical characteristics of the open-loop transfer function. For simulations, PID autopilot design for heading control problems of a conventional ship is considered. PID autopilot assures null stationary error for both step and ramp variations on reference inputs.
SuvaG. J.DattaA.BhattacharyyaS. P.New results on the synthesis of PID controllers20024722412522-s2.0-003647503210.1109/9.983352YuzJ. I.SalgadoM. E.From classical to state-feedback-based controllers200323458672-s2.0-004252754510.1109/MCS.2003.1213604SujitjornS.WiboonjaroenW.State-PID feedback for pole placement of LTI systems20112011202-s2.0-8155519526510.1155/2011/929430929430HerjólfssonG.Soffía HauksdóttirA.Direct computation of optimal PID controllersProceedings of the 42nd IEEE Conference on Decision and ControlDecember 2003112011252-s2.0-154229915610.1109/CDC.2003.1272718AstromK. J.HagglundT.1995Research Triangle Park, NC, USAInstrument Society of AmericaKesslerC.Das symmetrische optimum19586395400 and 432–436AstromK. J.HagglundT.Automatic tuning of simple regulators with specifications on phase and amplitude margins198420564565110.1016/0005-1098(84)90014-1TangY.OrtegaR.Adaptive tuning to frequency response specifications1993296155715632-s2.0-002770182010.1016/0005-1098(93)90021-KAb MalekM. N.Mohamed AliM. S.Evolutionary tuning method for PID controller parameters of a cruise control system using metamodeling2009200982-s2.0-7764925620110.1155/2009/234529234529WangL.FangX.DuanS.LiaoX.PID controller based on memristive CMAC network20132013651023810.1155/2013/510238DattaA.HoM. T.BhattacharyyaS. P.2000New York, NY, USASpringerWangQ.-G.ZhangZ.AstromK. J.ChekL. S.Guaranteed dominant pole placement with PID controllers20091923493522-s2.0-5874909485210.1016/j.jprocont.2008.04.012ValasekM.OlgacN.Efficient pole placement technique for linear time-variant SISO systems199514254514582-s2.0-002937090510.1049/ip-cta:19951959de la SenM.ON pole-placement controllers for linear time-delay systems with commensurate point delays2005200511231402-s2.0-1764440072610.1155/MPE.2005.123ZieglerJ. G.NicholsN. B.Optimum settings for automatic controllers194264759768PapadopoulosK. G.MermikliK.MargarisN. I.Optimal tuning of PID controllers for integrating processes via the symmetrical optimum criterionProceedings of the 19th Mediterranean Conference on Control and Automation (MED '11)June 2011grc128912942-s2.0-8005235349710.1109/MED.2011.5983018LoronL.Tuning of PID controllers by the non-symmetrical optimum method19973311031072-s2.0-0030817471PreitlS.PrecupR.-E.Extension of tuning relations after symmetrical optimum method for PI and PID controllers19993510173117362-s2.0-003259876110.1016/S0005-1098(99)00091-6NicolauV.2004Galati, RomaniaUniversity of Galati(Romanian)FossenT. I.1994New York, NY, USAJohn Wiley & SonsNomotoK.Response analysis of manoeuvrability and its application to ship design196611CeangaE.NichitaC.ProtinL.CutululisN. A.2001Bucharest, RomaniaTehnica