Next Article in Journal
Extending the Convergence Domain of Methods of Linear Interpolation for the Solution of Nonlinear Equations
Next Article in Special Issue
From Symmetry to Asymmetry: The Use of Additional Pulses to Improve Protection against Ultrashort Pulses Based on Modal Filtration
Previous Article in Journal
Solution and Interpretation of Neutrosophic Homogeneous Difference Equation
Previous Article in Special Issue
Artificial Intelligence: An Energy Efficiency Tool for Enhanced High performance computing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Overshoot Elimination for Control Systems with Parametric Uncertainty via a PID Controller

1
Department of Automation and Robotics, National Research Tomsk Polytechnic University, Tomsk 634050, Russia
2
University Administration Services, Tomsk State University of Control Systems and Radioelectronics, Tomsk 634050, Russia
*
Author to whom correspondence should be addressed.
Symmetry 2020, 12(7), 1092; https://doi.org/10.3390/sym12071092
Submission received: 3 June 2020 / Revised: 22 June 2020 / Accepted: 28 June 2020 / Published: 1 July 2020
(This article belongs to the Special Issue Information Technologies and Electronics)

Abstract

:
One of the key performance requirements for different control systems is non-overshooting step response, so that the controllable value should not overcome the reference value within a transient process. The problem of providing a non-overshooting step response was examined in this paper. Despite much scientific research being dedicated to the overshoot elimination problem, there are little to no results regarding parametric uncertainty for the discussed problem. Consideration of parametric uncertainty, particularly in the form of interval-given parameters, is essential, since in many physical processes, electronic devices and control systems parameter values can be obtained with acceptable error, and they can vary under different conditions. The main result of our research is the development of a proportional-integral-derivative (PID)-controller tuning approach for systems with interval-given parameters that provides a non-overshooting step response for such classes of control systems. This paper investigates analytical conditions and constraints for linear time invariant (LTI) systems in order to have no overshoot, enhances them with respect to parametric uncertainty, and formulates rules for tuning choices of parameters.

1. Introduction

Regarding industrial process control, it is essential for control systems to meet the technological requirements. One of the key performance indices is overshoot. Particularly, it is essential for positioning control systems, machinery, and several thermal and chemistry processes to show no overshoot. Considering linear time invariant (LTI) systems, there has been much research worldwide that has been dedicated to the problem of overshoot. Papers [1,2] present approaches for achieving acceptable overshoot percentages on the poles domination theory basis. Several research papers, such as [3,4,5,6], have been focused on full overshoot elimination using Vandermonde-like matrices for error function representation. Other ways were proposed in [7], where the problem was considered with respect to impulse response behavior, and in [8], where minimizing of integrated absolute error (IAE) was performed. In the context of industrial process control, it is worth noting different approaches based on canonical proportional-integral-derivative (PID)-controller design methods. In particular, several results in [9,10,11,12,13] were obtained using the Ziegler–Nichols tuning method. More theoretical techniques take place as well, such as [14,15] that state and prove theorems for necessary and sufficient pole-zero configuration of the transfer function, providing a non-overshooting step response. As opposed to several analytical approaches, numeric optimization for the purpose of PID-controller tuning was regarded as well in [16,17]. Another option for controller design and tuning is the initial choice of a desired stability degree for each pole of a closed-loop control system.
Considering real-world practical applications, it is worth mentioning that systems parameters tend to be uncertain, due to measurement errors and varying conditions of functioning. Even slight parameters variations can lead to significant changes in system dynamics. Parametric uncertainty makes controller design significantly more difficult, since the controller has to meet the requirements for varying states. In addition, a robust controller tuning procedure is challenging, since many existing well-known methods become inapplicable under parametric uncertainty. Despite the high level of complexity for a robust controller design, there have been several results presented, such as [18,19] that are based on minimization of oscillatory degree for systems with interval-given parameters. Regarding systems with interval-given parameters, one of the key techniques for analysis and control design is the Kharitonov theorem. For example, the authors of [20,21,22] address vertex polynomial properties for a robust PID-control design, but in many cases, vertex polynomials do not indicate all behavior aspects of the system with interval-given parameters. Another example of a robust PID-controller design was presented in [23,24] in the context of DC converters on the mixed approach that involves Kharitonov polynomials and phase portraits.

2. Materials and Methods

The presented research is conducted on the basis of the canonical closed-loop control scheme, depicted in Figure 1.
In Figure 1, C ( s ) is the controller, G ( s ) is the plant, and u , ε , and y are the reference signal, error signal, and output signal respectively.
The controller, C ( s ) in this case, is a classic PID controller, which is one of the most distributed control solutions. The PID-controller transfer function is
C ( s ) = D s 2 + K s + I s ,
where D , K and I are the coefficients of the PID controller. The plant G ( s ) is represented with the second-order transfer function
G ( s ) = K P s 2 + 2 α s + α 2 + ω 2 ,
with poles s = α ± j ω and gain K P . In cases of plants with initially real poles, parameter ω is supposed to be complex with a null real part.
At the first stage of the research, a stationary LTI system was examined. The task was to figure out the rule for PID-controller coefficients, providing a non-overshoot step response. In previous research [14], one of the necessary conditions for this was that the poles of the transfer function are strictly real.

2.1. Providing Real Closed-Loop Pole Configuration

Regarding control systems structure, plant and controller transfer functions, the closed loop transfer function is
W ( s ) = K P ( D s 2 + K s + I ) s 3 + ( 2 α + D K P ) s 2 + ( α 2 + ω 2 + K K P ) s + I K P
Conducting changing D K P = K , K K P = K , and I K P = I , a characteristic equation of (2) can be represented as:
s 3 + ( 2 α + D ) s 2 + ( α 2 + ω 2 + K ) s + I
According to the theorem in [14], it is necessary for transfer function poles to be exclusively real. Since the roots type is defined with a discriminant sign, let us write a discriminant for (3) in the following form [25]
Δ C E ( α , ω , K , I , D ) = ( 2 α + D ) 2 ( α 2 + ω 2 + K ) 2 27 I 2 4 I ( 2 α + D ) 3 4 ( α 2 + ω 2 + K ) 3 + 18 I ( 2 α + D ) ( α 2 + ω 2 + K )
Substituting arbitrary real positive values of α , ω , and K in (4), one can plot the surface of Δ C E ( D ,   I ) . Since only positive values of Δ C E are of interest, in the same coordinates, let us plot plane Δ C E = 0 to visualize constraints of the region of interest.
Figure 2 illustrates that region where Δ C E > 0 , denoted as Ω , is bounded with two curves (white dashed lines). Expressing I from (4) yields two functions describing bounding curves for Ω :
{ I 1 ( α , ω , K , D ) = 18 α K + 9 D K + 18 ω 2 α + 9 ω 2 D 12 α D 2 15 α 2 D + 2 α 3 2 D + 2 ( α 2 3 ω 2 + 4 α D + D 2 3 K ) 3 27 ; I 2 ( α , ω , K , D ) = 18 α K + 9 D K + 18 ω 2 α + 9 ω 2 D 12 α D 2 15 α 2 D + 2 α 3 2 D 2 ( α 2 3 ω 2 + 4 α D + D 2 3 K ) 3 27 .
For further investigation, let us plot region Ω in the D I plane as it is shown in Figure 3.
Lines of constraints for region Ω have an intersection point ξ ( D ξ , I ξ ) . Point ξ ( D ξ , I ξ ) defines peak values for parameters I and D .
For generalization purposes, let us get analytical expressions for I ξ and D ξ using (5)
D ξ = 2 α + 3 ( α 2 + ω 2 + K )
In the same way, substitution of (6) in (5) gives the expression [25]
I ξ = 3 ( α 2 + ω 2 + K ) 3 9 .
Let us investigate region Ω in more detail. It is known that for the third-order polynomial, its roots can be analytically calculated with the Cardano formula. Regarding the form of (2), the dominating pole can be calculated as follows
S D O M = u + v 3 + u v 3 a ,
where:
u = 9 ( 2 α + D ) ( α 2 + ω 2 + K ) 27 I 2 ( 2 α + D ) 3 54 ,
v = 27 ( 27 I 9 ( 2 α + D ) ( α 2 + ω 2 + K ) + 2 ( 2 α + D ) 3 ) + 4 ( 3 ( α 2 + ω 2 + K ) ( 2 α + D ) 2 ) 3 54 ,
a = ( 2 α + D ) 3 .
Investigating derivatives S D O M I and S D O M D in the field of real values states that (8) is a monotone, non-increasing function and thus, the stability degree η ( α , ω , K , I , D ) = | Re ( S D O M ( α , ω , K , D , D ) ) | will increase according to I . Regarding the I range of values, one can assume that the maximum stability degree can be reached in ξ point. Substituting (5) and (7) in (8) gives the following expression
η ξ M A X ( α , ω , K ) = 3 ( α 2 + ω 2 + K ) 3 .
Function (9) defines the maximum stability degree that can be reached for a given plant, with the chosen PID-controller proportional coefficient with strictly real poles for a closed-loop transfer function. In addition, it should be noticed that in ξ point, there is a triple pole T 1 = T 2 = T 3 = T .

2.2. Non-Overshoot Step Response Condition

Regarding special point ξ , let us check non-overshooting conditions that were formulated in [15]. According to [15], it is necessary and sufficient that if at least one of the following conditions holds, then system step response has no overshoot.
{ T 2 ( K I T ) T D I 0 , T ( D I T 2 ) 0 ;
{ T 2 ( K I T ) T D I 0 T 3 ( 2 T 2 ( K I 5 2 T ) ) + ( T ( D I 2 T 2 ) ) 2 0
For a more convenient robust controller synthesis procedure, let us find such a K value that guarantees a non-overshoot step response for any p ( D * , I * ) Ω . Investigating (10) and (11), it is clear that expression
φ ( T i , D , K , I ) = T i 2 ( K I T i ) T i D I ,
where T i , i , i 1 , 3 ¯ , is a part of each of the conditions (10) and (11). For generalization purposes, let us solve (12) with respect to T i . The solution is an expression that can be written as follows
T i = K ± K 2 4 D I 2 I .
Regarding the highest coefficient sign and (13), one can infer that
{ φ ( T i , D , K , I ) 0 ,   T i [ K K 2 4 D I 2 I , K + K 2 4 D I 2 I ] φ ( T i , D , K , I ) < 0 ,   T i ( ,   K K 2 4 D I 2 I ) , (   K + K 2 4 D I 2 I , + ) .
Since K , D , and I are strictly positive, regarding (14), one can assume that the condition φ ( T i , D , K , I ) < 0 will always hold, and it is a monotonic decreasing function of T i . Next, substitution of (6, 7, 9) into (11) yields that
K 2 α 2 + 6 α ω + 2 ω 2
Regarding the case of real poles of plant transfer function and, thus, complex value of ω with a null real part, the aforementioned substitution and expressing K gives:
K 16 7 α 2 + 30 7 α ( 5 7 α + 1 7 4 α 2 7 ω 2 )
In other words, a proper choice of K provides a non-overshooting step response for a closed-loop system within every point p ( D * , I * ) Ω .

2.3. Plant with Interval-Given Parameters

Regarding practical applications, the exact parameters values are unknown and, basically, can be represented as a confidence limit. Moreover, parameters tend to vary due to temperature, humidity, pressure, and mechanical deterioration. In addition, in the field of outdoor mobile robotics, conditions of functioning are highly heterogeneous and demand special approaches regarding varying parameters [26,27,28]. Thus, transfer function parameters basically have to be represented as an interval value. Since further mathematical representation of the plant contains interval-given values, the controller synthesis should be conducted accordingly.
Let us consider second-order transfer function with parametric uncertainty as a plant model
W ( s ) = K P a s 2 + b s + c ,
where a [ a _ ,   a ¯ ] , b [ b _ ,   b ¯ ] , c [ c _ ,   c ¯ ] , and K P [ K P _ ,   K P ¯ ] are given intervals.
According to [14,15], a system with interval-given parameters is basically a set of stationary LTI systems, and its poles’ location can be represented as a region with borders defined by ranges of interval parameters. For the purpose of analysis and control design for systems with interval-given parameters, it is sufficient only to consider external borders of poles localization region (see Figure 4). A typical well-known representation of the poles localization region is multiparametric interval root locus (MIRL). Since possible plant pole configurations are real poles or a complex–conjugate pair, MIRL is always symmetric with respect to the x-axis. The symmetry property of MIRL allows us to simplify further research, i.e., for the second-order transfer function, it is sufficient to investigate only one half of the MIRL. In other words, each symmetrical pair of points that belongs to MIRL is a poles pair for the same LTI system. In addition, since only half of MIRL is sufficient, computational load is reduced.
Since every point that belongs to MIRL is a single LTI system that has its own region Ω with ξ point, the region of poles localization forms a corresponding region of ξ points on D I plane. With accordance to (6) and (7), I ξ and D ξ values are defined by K = K K P , and due to the interval nature of K P , for every K P [ K P _ , K P ¯ ] mapping, M ( K P * ) ,   K P * [ K P _ , K P ¯ ] can be obtained (see Figure 5). According to (5), regions of positive discriminant values Ω can be obtained for ξ M ( K P * ) ,   K P * [ K P _ , K P ¯ ] . The main aim of the research is to obtain PID-controller coefficient values such that a closed-loop system has a non-overshooting step response for any plant parameter value variations within given ranges. Since every point with coordinates ( D , I ) Ω provides strictly real closed-loop poles for (2) for corresponding point ξ , then some pair of controller coefficient values ( D * , I * ) gives { s ,   s < 0   | ξ i K p M ( K P * ) , K P * [ K P _ , K P ¯ ] } in case ( D * , I * ) Ω * , such that Ω * = Ω i K p for every ξ i K p i .
Similarly to each Ω region, for resultant region Ω * , constraints are defined with (5) for particular values of ξ . Thus, in order to choose ( D * , I * )   properly, it is sufficient to define points ξ 1 and ξ 2 , for which (5) forms constraints for the desired set Ω * , which contains I and D values that provide a non-overshooting step response under parameter variation.

2.4. Constraints Clarification for Ω *

The next problem in the research is to generalize constraints for Ω * . It can be suggested that Ω * is constrained with (5) for arguments providing the I ξ value to be minimum according to (6) maximum D I 0 —the value of D that turns I 2 ( α , ω , k , D ) into zero. The value of D I 0 is defined as follows:
D I 0 ( α , ω , K ) = 2 α + 2 ( α 2 + ω 2 + K ) .
It should be noted that desired points ξ 1 and ξ 2 can either belong to vertices or edges mapping of MIRL. For further convenience, (7) and (18) can be rewritten as functions of plant model parameters. Since (1) and (17) are equivalent representations of the plant transfer function, (7) and (18) can be represented as follows:
I ξ ( a , b , c , K ) = 1 a 3 3 ( 4 a 4 + 4 K a 2 + 4 c a b 2 ) 3 72 ,
D I 0 ( a , b , c , K ) = 2 c a + K b a .
Regarding the assumption that desired points ξ 1 and ξ 2 belong to edges of MIRL mapping M ( K P i ) , it should be noted that the only possible way of edge location for ξ 1 and ξ 2 is a nonmonotonic behavior of (19) or (20). It is known that each edge of MIRL and, consequently, its mapping M ( K P ) are formed by the varying of single interval parameters with others fixed in their limit values, and it turns (19) and (20) into one-variable functions forming each edge.
In order to define whether (19) and (20) are monotonic or not along each interval parameter, one can find partial derivatives for (19) and (20) with respect to a , b , c , and K for checking their monotonicity property. Partial derivation for (19) yields:
I ξ ( a , b , c , K ) a = 3 ( 4 a 4 2 c a + b 2 ) ( 4 a 4 + 4 K a 2 + 4 c a b 2 ) 2 24 a 4 ( 4 a 4 + 4 K a 2 + 4 c a b 2 ) 3
I ξ ( a , b , c , K ) b = 3 b ( 4 a 4 + 4 K a 2 + 4 c a b 2 ) 2 24 a 3 ( 4 a 4 + 4 K a 2 + 4 c a b 2 ) 3
I ξ ( a , b , c , K ) c = 3 ( 4 a 4 + 4 K a 2 + 4 c a b 2 ) 2 12 a 2 ( 4 a 4 + 4 K a 2 + 4 c a b 2 ) 3
I ξ ( a , b , c , K ) c = 3 ( 4 a 4 + 4 K a 2 + 4 c a b 2 ) 2 12 a ( 4 a 4 + 4 K a 2 + 4 c a b 2 ) 3
Investigating (21)–(24), it can be concluded that if 4 a 4 + 4 K a 2 + 4 c a b 2 and 4 a 4 2 c a + b 2 possess a non-zero value at any a [ a _ ,   a ¯ ] , b [ b _ ,   b ¯ ] , c [ c _ ,   c ¯ ] , and K [ K _ , K ¯ ] , then (19) has a monotonic behavior. Alternatively, for solutions for one of (21), (22), (23), or (24) with respect to changing variables, one can get the exact investigated parameter value that gives ξ 1 .
Similarly, partial derivation for (20) gives:
D I 0 ( a , b , c , K ) a = b c c a + K 2 a 2 ,
D I 0 ( a , b , c , K ) b = 1 a ,
D I 0 ( a , b , c , K ) c = 1 a c a + K 2 ,
D I 0 ( a , b , c , K ) K = 1 c a + K .
Functions (26)–(28) are always non-zero, so that it is sufficient to check the behavior of (25). In addition, one has to check (25) only in the case of parameter a variation.
Regarding the aforementioned results, if (21)–(28) are monotonic functions within interval parameters bounds, in special points ξ 1 ( α L ,   ω L ) and ξ 2 ( α U ,   ω U ) that define regions in which intersection gives Ω * , only two members of interval family form the set containing I and D values that satisfies non-overshoot conditions for all the interval family. The two special members of the family are defined as follows:
W ξ 1 ( s ) = K P _ a ¯ s 2 + b _ s + c ¯ ,
W ξ 2 ( s ) = K P ¯ a _ s 2 + b _ s + c ¯ .
In case some of (21)–(28) are non-monotonic, one has to find the exact solution with respect to varying parameters and substitute it to (29) or (30), instead of limiting the value of the parameter.

2.5. PID-Controller Coefficient Choice

The final problem within the present research is to find the exact values for the PID-controller coefficients D and I . One has to note that D = K P D and I = K P I are interval values, since K P is a given interval. Thus, it is essential to find D and I values such that:
{ D Ω * ,   I Ω *   |   K P [ K P _ ; K P ¯ ] } .
In order to satisfy (31), one can check two points denoted as P U ( D K P ¯ , I K P ¯ ) and P L ( D K P _ , I K P _ ) , since any point ( D K P * , I K P * ) ,   K P [ K P _ ; K P ¯ ] belongs to the linear interval P L U . Let us investigate the generalized condition for P L U to lay within the Ω * region. In other words, the following problem is to choose D and I in such a way that P L U Ω * is as depicted in Figure 6.
According to Figure 6, one can observe that with a known range of gain K P and known expression for constraints (5), D and I values providing P L U Ω * can be found from inequality of the form:
I U ( α U , ω U , K U , D U ) K P ¯ I L ( α L , ω L , K L , D L ) K P _
Regarding the limit case, i.e., P L I 2 and P U I 1 , the aforementioned inequality can be represented as the following equality:
I U ( α U , ω U , K U , D U ) K P ¯ = I L ( α L , ω L , K L , D L ) K P _ .

3. Example

In order to verify the obtained results, let us consider arbitrary second-order transfer functions with interval-given parameters that are described with the expression:
W ( s ) = [ 16 , 18 ] [ 1.2 , 1.4 ] s 2 + [ 8.3 ,   9.7 ] s + [ 28 ,   36 ]
Using MATLAB, let us plot the MIRL of (33). The MIRL for (33) is shown in Figure 7.
Regarding the results obtained in [19], one can obtain the external border of the MIRL presented in Figure 7 by vertices numbered as 3 ,   4 ,   2 ,   6 ,   5 ,   7 , and corresponding edges that connect external vertices. The candidate transfer functions according to (29) and (30) for the plant are:
W ξ 1 ( s ) = 16 1.4 s 2 + 8.3 s + 36
W ξ 1 ( s ) = 18 1.2 s 2 + 8.3 s + 36
The investigation according to (21)–(28) yields that candidate functions are satisfied in the aforementioned conditions and, thus, define Ω * . W ξ 1 and W ξ 2 parameters are: α L = 3.4643 , ω L = 2.8282 , α U = 3.4526 , and ω U = 4.2461 . With respect to W ξ 1 and W ξ 2 as proportional coefficients for PID controllers, the lowest values should be chosen, and according to (14), one can obtain K = 6.5808 , and thus, K = K K P = [ 75.2091 ,   98.712 ] .
For the investigated plant, regarding the external border presented in Figure 7, mapping M ( K P ) for different K P values within its range was calculated and plotted (see Figure 8).
Substituting α L = 3.4643 , ω L = 2.8282 , α U = 3.4526 , and ω U = 4.2461 into (5), one can obtain analytical expression for Ω * constraints and that gives:
I U ( D ) = 32.2913 D 2 ( D 2 + 13.8104 D 338.3036 ) 3 27 1.54 D 2 0.074 D 3 + 271.756
I L ( D ) = 21.0685 D + 2 ( D 2 + 13.8527 D 237.6221 ) 3 27 1.53 D 2 0.074 D 3 + 195.251
Next, solving (32) with respect to variable D yields P U ( 19.9695 , 90.8278 ) and P L ( 17.1167 ,   77.8524 ) , calculation of the PID-controller coefficients gives the following results: D = 1.331 , K = 6.5808 , and I = 6.055 . For the obtained values, the condition P L U Ω * holds and provides a non-overshooting closed-loop step response for (2). The obtained closed-loop transfer function has the following form
W C L ( s ) = [ 11.43 ;   15 ]   ( 1.331 s 2 + 6.581 s + 6.055 ) s 3 + [ 21.1379 ;   28.049 ]   s 2 + [ 95.2234 ;   128.712 ]   s + [ 69.1985 ;   90.8226 ]
For confirmation results, let us obtain step responses for W C L ( s ) in different realizations, i.e., with different parameters value combinations. In Figure 9, fifty different members of the interval family are presented.
In Figure 9, there are a set of step responses for W C L ( s ) that are obtained with regular steps of each parameter combination, and the overshoot value is zero for each of the responses.

4. Conclusions

This paper presents the approach for PID-controller tuning providing a non-overshoot step response for a second-order transfer function with interval-given parameters. Designed for a stationary LTI-system approach based on necessary conditions of closed-loop poles, a real configuration with further application is necessary, and sufficient conditions within the region Ω of PID-controller parameters guarantees strictly real closed-loop poles. Next, the rule for proportional coefficient calculation was introduced while holding necessary and sufficient conditions in every point that belong to region Ω . In addition, another key performance index in many control applications is the stability degree [29,30]. In this paper, (9) defines the maximum stability with respect to plant parameters and chosen proportional PID-controller coefficients and thus allows us to tune the PID controller according to desired stability degree.
Considering parametric uncertainty and its interval representation, the approach was enhanced and applied to this class of systems. Since in systems with interval-given parameters, poles move and form regions of localization, the approach has to regard a certain set of cases. Using (5) and (6), one can get mapping of poles’ localization regions into a PID-controller parameter plane. The next step was getting the region of PID-controller coefficients that allows maintaining the aforementioned conditions for each member of interval family. Finally, regarding all constraints and conditions, the equation whose solution gives the exact PID-controller coefficient was formulated.
One of the key disadvantages is difficulty of enhancing of the presented results for higher-order systems due to lack of analytical methods for a higher-order equation solution. Using the Ferrari method for a fourth-order equation, plant model objects can probably be third-order. On the other hand, a wide class of typical industrial plants can be described with second- or third-order transfer function [31] that makes the presented approach applicable for many problems.
One of the main advantages of the approach is a high degree of generality, since all basic functions and expression are presented in general form, and one can investigate how exactly the parameters of the controller or plant affect the design procedure and, thus, the final result. The generalization of design techniques allows us to eliminate so-called “loss of insight”, in which [32] corresponded to a numerical calculation. As a direction of further research, generalizing results for transfer functions of higher order with zeros should be examined. Another topic of further research can be PID-controller design with consideration of settling time, energy, and more complex kinds of uncertainties as well. Moreover, regarding future directions of research, one thought is to enhance it to linear time-varying (LTV) systems, since despite being in the class of LTV systems, its parameters are not just regular intervals, but functions of time. Due to fundamental physical constraints, ranges of definition for such kind of functions are constrained as well. The presented research forms the minimal necessary basis for providing a required performance index within the given interval, and the possible next step is to advance control quality with respect to law of parameters variation and other conditions [33].

Author Contributions

Conceptualization, A.T. and S.E.; methodology, A.T.; validation, S.E.; investigation, A.T. and S.E.; software, A.T.; visualization, A.T.; writing—original draft preparation, A.T.; writing—review and editing, S.E.; formal analysis, S.E. and S.Z.; data curation, S.E.; supervision, S.Z.; project administration, S.Z.; funding acquisition, S.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research and APC was funded by Russian Ministry of Science and Higher Education grant number FEWM-2020-0036.

Acknowledgments

Authors would like to thank National Research Tomsk Polytechnic University for administrative and technical support. Additionally, authors would like to thank Shir-Kuan Lin and Chang-Jia Fang for providing additional comments on their work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Moore, K.L.; Bhattacharyya, S.P. A technique for choosing zero locations for minimal overshoot. IEEE Trans. Autom. Control 1990, 35, 577–580. [Google Scholar] [CrossRef]
  2. Efimov, S.V.; Kurgankin, V.V.; Zamyatin, S.V. Designing transfer function with the required direct performance measures based on the Laplace transform. Optoelectron. Instrum. Data Process. 2014, 50, 348–353. [Google Scholar] [CrossRef]
  3. Hauksdottir, A.S. A sufficient and necessary condition for extrema-free step responses of single-zero continuous-time systems with real poles. In Proceedings of the 38 Conference on Decision & Control, Phoenix, AZ, USA, 7–10 December 1999; pp. 499–500. [Google Scholar]
  4. Hauksdottir, A.S. Analytic Expression of Transfer Function Reponses and Choice Numerator Coefficients (Zeros). IEEE Trans. Autom. Control 1996, 41, 1482–1488. [Google Scholar] [CrossRef]
  5. Hauksdottir, A.S.; Hjaltadottir, H. Closed-form expressions of transfer function responses. In Proceedings of the American Control Conference, Denver, CO, USA, 4–6 June 2003; pp. 3234–3239. [Google Scholar]
  6. Herjolfsson, G.; Hauksdottir, A.S. Direct computation of optimal PID controllers. In Proceedings of the 42nd IEEE Conference on Decision and Control, Maui, HI, USA, 9–12 December 2003; pp. 1120–1225. [Google Scholar]
  7. Jayasuriya, S.; Dharne, A.G. Necessary and sufficient conditions for non-overshooting step responses for LTI systems. In Proceedings of the 2002 American Control Conference, Anchorage, AK, USA, 8–10 May 2002; pp. 500–510. [Google Scholar]
  8. Nguyen, N.H.; Nguyen, P.D. Overshoot and settling time assignment with PID for first-order and second-order systems. IET Control Theory Appl. 2018, 12, 2407–2416. [Google Scholar] [CrossRef]
  9. Prakash, R.; Anburaja, S.A.; Rishivanth, S.; Govinda Kumar, E. Non-overshoot time response of third order system using cascade PID-Lead compensator controller. In Proceedings of the 2017 Second International Conference on Electrical, Computer and Communication Technologies (ICECCT), Coimbatore, Tamil Nadu, India, 22–24 February 2017; pp. 1–6. [Google Scholar]
  10. Visioli, A. Practical PID Control, 1st ed.; Springer: London, UK, 2006; p. 314. [Google Scholar]
  11. Anbarasan, M.; Suji Prasad, S.J.; Meenakumari, R.; Balakrishnan, P.A. Modified PID controller for avoiding overshoot in temperature of barrel heating system. In Proceedings of the 2013 International Conference on Emerging Trends in VLSI, Embedded System, Nano Electronics and Telecommunication System (ICEVENT), Tiruvannamalai, Tamil Nadu, India, 7–9 January 2013; pp. 1–6. [Google Scholar]
  12. Souza, M.R.S.B.; Murofushi, R.H.; Tavares, J.J.P.Z.; Ribeiro, J.F. PID Tuning for the Pitch Angle of a Two-Wheeled Vehicle. In Proceedings of the 2015 12th Latin American Robotics Symposium and 2015 3rd Brazilian Symposium on Robotics (LARS-SBR), Uberlândia, Brazil, 28 October–1 November 2015; pp. 371–375. [Google Scholar]
  13. Wu, H.; Su, W.; Liu, Z. PID controllers: Design and tuning methods. In Proceedings of the 2014 9th IEEE Conference on Industrial Electronics and Applications (ICIEA), Hangzhou, China, 9–11 June 2014; pp. 808–813. [Google Scholar]
  14. Kobayashi, H. Output overshoot and pole-zero configuration. In Proceedings of the 12th IFAC World Congress Automat. Contr., IFAC, Sydney, Australia, 18–23 July 1993; pp. 529–532. [Google Scholar]
  15. Lin, S.K.; Fang, C.J. Nonovershooting and monotone nondecreasing step responses of a third-order SISO linear system. IEEE Trans. Autom. Control 1997, 42, 1299–1303. [Google Scholar]
  16. Zhmud, V.; Yadrishnikov, O.; Poloshchuk, A.; Zavorin, A.A. Modern key technologies in automatics: Structures and numerical optimization of regulators. In Proceedings of the 2012 7th International Forum on Strategic Technology (IFOST), Tomsk, Russia, 18–21 September 2012. [Google Scholar]
  17. Prokopiev, A.P.; Ivanchura, V.I.; Emelyanov, R.T. Synthesis PID Controller for Objects Second Order with Regard to the Location Poles. J. Sib. Fed. Univ. Eng. Technol. 2016, 1, 50–60. [Google Scholar] [CrossRef]
  18. Nesenchuk, A.A. A method for synthesis of robust interval polynomials using the extended root locus. In Proceedings of the 2017 American Control Conference (ACC), Seattle, WA, USA, 24–26 May 2017; pp. 1715–1720. [Google Scholar]
  19. Sukhodoev, M.S.; Gayvoronskiy, S.A.; Zamyatin, S.V. Parametric synthesis of linear regulator in interval system with guaranteed root quality indices. Bull. Tomsk Polytech. Univ. 2007, 311, 9–12. [Google Scholar]
  20. Ouyang, H.; Yue, J.; Su, Y. Design and Application of PID Controllers Based on Interval Computing Theory. In Proceedings of the 2011 International Conference on Electrical and Control Engineering, Yichang, China, 16–18 September 2011; pp. 1505–1510. [Google Scholar]
  21. Chapellat, H.; Bhattacharyya, S.P. A generalization of Kharitonov’s theorem; Robust stability of interval plants. IEEE Trans. Autom. Control 1989, 34, 306–311. [Google Scholar] [CrossRef]
  22. Zhmud, V.; Zavorin, A. The design of the control system for object with delay and interval-given parameters. In Proceedings of the 2015 International Siberian Conference on Control and Communications (SIBCON), Omsk, Russia, 21–23 May 2015. [Google Scholar]
  23. Marcillo, K.E.L.; Guingla, D.A.P.; Barra, W.; De Medeiros, R.L.P.; Rocha, E.M.; Benavides, D.A.V.; Nogueira, F.G. Interval robust controller to minimize oscillations effects caused by constant power load in a DC multi-converter buck-buck system. IEEE Access 2019, 7, 26324–26342. [Google Scholar] [CrossRef]
  24. Marcillo, K.E.L.; Rocha, E.M.; Junior, W.B.; Benavides, D.A.V.; de Medeiros, R.L.P.; Alves, M.S.; Leão, M.P.V. Robust control for DC–DC buck converter under parametric uncertainties. In Proceedings of the 22nd Congresso Brasileiro de Automatica (CBA), Paraiba, Brazil, 9–12 September 2018; pp. 1–8. [Google Scholar]
  25. Tsavnin, A.V.; Efimov, S.V.; Zamyatin, S.V. Providing real closed-loop transfer functions poles for plant with interval-given parameters for overshoot elimination. In Proceedings of the 11th International Congress on Ultra Modern Telecommunications and Control Systems and Workshops (ICUMT), Dublin, Ireland, 28–30 October 2019; pp. 1–7. [Google Scholar]
  26. Andrakhanov, A.; Belyaev, A. Navigation learning system for mobile robot in heterogeneous environment: Inductive modeling approach. In Proceedings of the 12th International Scientific and Technical Conference on Computer Sciences and Information Technologies (CSIT 2017), Lviv, Ukraine, 5–8 September 2017; pp. 543–548. [Google Scholar]
  27. Andrakhanov, A.; Belyaev, A. GMDH-Based Learning System for Mobile Robot Navigation in Heterogeneous Environment. In Advances in Intelligent Systems and Computing; Shakhovska, N., Stepashko, V., Eds.; Springer: Cham, Switzerland, 2018; Volume 689, pp. 1–20. [Google Scholar]
  28. Fadeev, A.S.; Zarnitsyn, A.Y.; Tsavnin, A.V.; Belyaev, A.S. Cyber-physical system prototype development for control of mobile robots group for general mission accomplishment. In Proceedings of the Second International Conference on Material Science, Smart Structures and Applications: ICMSS-2019, Erode, India, 21–22 November 2019. [Google Scholar]
  29. Pushkarev, M.I.; Gaivoronsky, S.A.; Efimov, S.V.; Zamyatin, S.V. Parametric synthesis of maximum stability degree and specified accuracy linear automatic control system PI-controller. In Proceedings of the 18th International Conference on Soft Computing, MENDEL 2012, Brno, Czech Republic, 27–29 June 2012; pp. 344–349. [Google Scholar]
  30. Pushkarev, M.I.; Efimov, S.V.; Gayvoronsky, S.A.; Chenkova, O.A. A single-loop DC motor control system design with a desired aperiodic degree of stability. In Proceedings of the International Conference on Mechanical Engineering, Automation and Control Systems 2015 MEACS 2015, Tomsk Polytechnic University, Tomsk, Russian, 1–4 December 2015. [Google Scholar]
  31. Astrom, K.J.; Hagllund, T. PID Controllers Theory, Design and Tuning, 2nd ed.; Instrument Society of America: Pittsburgh, PA, USA, 1994. [Google Scholar]
  32. Barmish, B.R. New Tool for Robustness of Linear Systems; Macmillan Publishing Company: New York, NY, USA, 1994. [Google Scholar]
  33. Xiang, W.; Xiao, J. Stabilization of switched continuous-time systems with all modes unstable via dwell time switching. Automatica 2014, 3, 940–945. [Google Scholar] [CrossRef]
Figure 1. Typical closed-loop structure.
Figure 1. Typical closed-loop structure.
Symmetry 12 01092 g001
Figure 2. Surface of discriminant values with zero plane.
Figure 2. Surface of discriminant values with zero plane.
Symmetry 12 01092 g002
Figure 3. Region Ω in a plane of controller coefficients D I .
Figure 3. Region Ω in a plane of controller coefficients D I .
Symmetry 12 01092 g003
Figure 4. Multiparametric interval root locus for second-order plant with parametric uncertainties.
Figure 4. Multiparametric interval root locus for second-order plant with parametric uncertainties.
Symmetry 12 01092 g004
Figure 5. Set of mappings M ( K P i ) of multiparametric interval root locus (MIRL) with respect to the K P range.
Figure 5. Set of mappings M ( K P i ) of multiparametric interval root locus (MIRL) with respect to the K P range.
Symmetry 12 01092 g005
Figure 6. P L and P U localization within Ω * and corresponding D and I ranges.
Figure 6. P L and P U localization within Ω * and corresponding D and I ranges.
Symmetry 12 01092 g006
Figure 7. Interval root locus for (33).
Figure 7. Interval root locus for (33).
Symmetry 12 01092 g007
Figure 8. Set of experimental mappings M ( K P i ) of the MIRL external border for (33) and obtained ξ 1 and ξ 2 .
Figure 8. Set of experimental mappings M ( K P i ) of the MIRL external border for (33) and obtained ξ 1 and ξ 2 .
Symmetry 12 01092 g008
Figure 9. Experimental set of step responses for obtained closed-loop transfer function.
Figure 9. Experimental set of step responses for obtained closed-loop transfer function.
Symmetry 12 01092 g009

Share and Cite

MDPI and ACS Style

Tsavnin, A.; Efimov, S.; Zamyatin, S. Overshoot Elimination for Control Systems with Parametric Uncertainty via a PID Controller. Symmetry 2020, 12, 1092. https://doi.org/10.3390/sym12071092

AMA Style

Tsavnin A, Efimov S, Zamyatin S. Overshoot Elimination for Control Systems with Parametric Uncertainty via a PID Controller. Symmetry. 2020; 12(7):1092. https://doi.org/10.3390/sym12071092

Chicago/Turabian Style

Tsavnin, Alexey, Semen Efimov, and Sergey Zamyatin. 2020. "Overshoot Elimination for Control Systems with Parametric Uncertainty via a PID Controller" Symmetry 12, no. 7: 1092. https://doi.org/10.3390/sym12071092

APA Style

Tsavnin, A., Efimov, S., & Zamyatin, S. (2020). Overshoot Elimination for Control Systems with Parametric Uncertainty via a PID Controller. Symmetry, 12(7), 1092. https://doi.org/10.3390/sym12071092

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