## PI controller design method with desired phase margin and settling time

14. Apríl, 2014, Autor článku: **Osuský Jakub**, Elektrotechnika

Ročník 7, číslo 4
Pridať príspevok

In this paper controller design method for SISO systems with performance specification in term of phase margin and desired settling time is presented. Mentioned method is based on modification of Neimark D-partition method which ensure desired phase margin of open loop and desired settling time by setpoint change. The developed frequency domain design technique is graphical, interactive and insightful. Theoretical results are demonstrated on examples.

**1. Introduction**

Frequency domain techniques for analysis and controller design dominate SISO control system theory. Bode, Nyquist, Nichols, and root locus are the usual tools for SISO system analysis.

Frequency methods are often used for controller design because it is easy to ensure performance of closed loop system through phase margin in open loop (Nagurka and Yaniv, 2003). To achieve the desired phase margin, controllers are usually designed using Bode characteristics (Fung, et al 1998),( Ho, et al 1995).

In this paper Neimarks method of D-partition is used (Neimark, 1992). This method is usually used for controller design which ensures closed loop stability with desired stability degree. In this paper is presented modification of this method aimed on phase margin instead of stability degree. Settling time is ensured using crossover frequency which has significant influence on this second performance specification.

**2. Preliminaries and problem formulation**

Consider classical feedback system depicted in Fig. 1.

*Figure 1: Classical feedback system*

The aim is to design PI controller for plant using Neimark D-partition method so, that not only stability will be ensured but performance in term of phase margin and desired settling time too.

**3. Theoretical results**

Consider closed loop feedback system with PI controller and SISO plant.

(1) |

Substituting (1) into characteristic function (2)

(2) |

and after small manipulations is possible to obtain:

(3) |

After other substitution and division into real and imaginary part:

(4) |

If is changing step by step in interval from real part of (4) is possible to calculate frequency dependent vector of complex numbers which plotted in complex plane create D-curve for parameter k. Similar is it with imaginary part of (4) from which is possible obtain k_{i}. D-curve plot for parameters k and k_{i} (PI controller) can be done in one step. This way is possible to calculate and paint marginal D-curves. If PI controller parameters will be chosen from this D-curves system with controller will have Nyquist plot intersecting point -1+0j in complex plane (Figure. 2).

*Figure. 2 Nyquist plot intersecting point -1+0j*

To obtain any phase margin it is necessary to move intersection point with unit circle form point -1+0j to point situated in lower left quadrant on the same circle. To obtain desired phase margin , (2) have to be modified as follows:

(5) |

where values of x and y depends on (Figure 3)

*Figure. 3 Nyquist plot with intersection ensuring phase margin .*

According to situation on Figure 3 and known radius of unit circle, values of x and y can be calculate as follows:

(6) |

For calculating D-curves of PI controller parameters k and k_{i}, ensuring desired phase margin (4) is modified to (7).

(7) |

Settling time of controlled system depends on plant gain and values of controller parameters. Higher values of k_{i} controller parameter provide shorter settling time for the price of higher controller output. D-curves of plants with possitive gain, starts in negative values of both PI controller parameters not far from origin. If plant is possible to stabilize with chosen phase margin using PI control structure, with increasing frequency D-curves continues to area where both PI controller parameters has possitive values (Fig. 4).

*Fig. 4 PI controller D-curves for different phase margin*

Chosing controller parameters from D-cruve also crossover frequency is chosen (frequency when open-loop system reaches unity circle). Controller parameters from D-curves maximum point M, yields shortest settling time. Consider parameters with higher frequency as maximum point frequency has no practical purpose. By increasing of crossover freqeuncy from (lowest frequency which yields possitive controller parameters) to , both controller parameters will increase and settling time T_{s} will decrease.

(8) |

According to (8) it is possible to define controller design algorithm, ensuring desired phase margin and settling time consisting of following steps:

- Chose desired phase margin and settling time T
_{s} - Calculate x and y according to (6)
- Calculate and frequency vectors from (7) and find as frequency where has maximum value
- Plot D-curve with k parameter on x-axis and k
_{i}on y-axis and design temporary PI controller for any crossover frequency - Measure settling time T
_{s_tmp}of system with designed temporary controller - Calculate frequency ensuring desired settling time T
_{s}as - If back to a) and chose higher T
_{s}or lower , else design final PI controller for crossover frequency

**4. Case study**

Consider stable 3-th order plant

a) For this plant PI controller ensuring phase margin and settling time T_{s}=6s will be designed.

b) According to (6) and desired phase margin:

c) D-curves for parameters k and k_{i} can be calculated in one step and plotted in one figure, because characteristic function is divided into real and imaginary part and each parameter is calculated from another part.

(8) |

d) In fig. 5 is plotted D-curve for parameters k and k_{i} calculated for . Parameter vector reaches the maximum value 1.5569 by frequency rad/s.

*Fig. 5 D-curve for phase margin *

Temporary controller has following parameters:

and from bode plot Fig. 6 it is possible to see that phase margin is really and crossover frequency rad/s.

*Fig. 6 Bode plot for open-loop system with temporary PI controller*

e) Settling time is the time required for the output to reach and remain within a 2% error band following input stimulus. This way for system with temporary controller was measured settling time T_{s_tmp}=7.3s (Fig. 7).

*Fig. 7 Step response for closed loop system with temporary PI controller*

f) Final crossover frequency is now calculated according to (9)

rad/s |

g) rad/s

From Fig. 8 is possible to see that phase margin is equal to desired phase margin and also crossover frequency is very close to calculated frequency .

*Fig. 8 Bode plot for open-loop system with final PI controller*

Settling time was measured from unity step response (Fig. 9) and the value is T_{s_fin}=6.15s what is very close to desired settling time T_{s}=6s.

*Fig. 9 Step response for closed loop system with final PI controller*

**5. Conclusions**

In this paper modification of Neimark D-partition method was presented. This controller design approach ensures not only stability but performance in term of phase margin and desired settling time too. The developed frequency domain design technique is graphical, interactive and insightful and it is useful for stable systems. Theoretical results have been verified on case study.

**Acknowledgments**

This work is supported by the Slovak Research and Development Agency (APVV) using financial support no. APVV-0772-12.

**References**

- Fung, H.W., Wang, Q.G., Lee T.H.: PI Tuning in terms of gain and phase margins, Automatica, 34, pp. 1145-1149, 1998.
- Ho, W.K., Hang, C.C.,Cao, L.: Tuning of PID controllers based on gain and phase margin specifications, Automatica, 31, pp. 497-502, 1995.
- Hudzovič, P.: Teória riadenia I., Alfa, 1980, Bratislava
- Nagurka M., Yaniv, O.: Robust PI Controller Design Satisfying Gain and Phase Margin Constraints. Proceedings of the American Control Conference, Denver, Colorado June 4-6, 2003
- Neimark, Y.I.: Robust stability and D-partition, Automation and Remote Control 53 (1992) (7), pp. 957–965

Coauthor of this paper is Stanislav Števo, Institute of Control and Industrial Informatics, Faculty of Electrical Engineering and Information Technology, Slovak University of Technology, Ilkovičova 3, 812 19 Bratislava, Slovak Republic