Robust Tuning of The Simplest Dead Time Compensators: Mikul A S Huba
Robust Tuning of The Simplest Dead Time Compensators: Mikul A S Huba
Robust Tuning of The Simplest Dead Time Compensators: Mikul A S Huba
Compensators
Mikuláš Huba
Slovak University of Technology in Bratislava, Slovakia
Email: mikulas.huba@stuba.sk
Abstract—This paper rediscovers and modifies the oldest and with a gain estimation error in the range 10-20% of the nominal
the simplest dead-time compensator (DTC) introduced originally value is shown. So, in total, the DTC performance is used to be
by [1] and shows new possibility for its robust tuning by the demonstrated in 3-6 different working points. This is, however,
performance portrait method. This method comprising computer still far from being enough for fully characterizing such a
based analysis, compressed characterization and visualization complex closed loop performance, as the DTC have, from a
of the closed loop properties achieved by evaluating set of
performance measures over chosen grid of loop parameters with
potential user view. In this paper the loop performance will
an optimization based on sweeping all possible stored data may be characterized by information extracted from hundreds, or
be successfully used for an optimal nominal or robust tuning thausends of experiments each one carried out for a different
[2]–[4] based either on the known plant model with parameters working point. The interesting shape and time related closed
taking values upon uncertainty intervals, or directly on measured loop properties will be mapped, stored, visualized and used
experimental data [5]. In order to enable an easy comparing with for the robust control design based on a compressed form
other possible solutions new robustness characteristics will be of information representation called here as the closed-loop
derived showing how a model uncertainty influences decrease of performance portrait (PP). The PP method (PPM) is then used
the loop performance. for comparing and optimizing properties of different DTCs.
will be chosen Such a task has been firstly treated for n = 1 A. Performance portrait method (PPM)
and Tp = Tn in [14]. Experience gained in dealing with DO
The second type of above mentioned limitations will be
based controllers for the first order plants gives motivation to
analyzed by the performance portrait method (PPM). To be
come back to this controller tuning with a question, how its
able to use the generated PP for any plant (5), the setpoint
closed loop robustness and noise attenuation depend on the
step responses will be mapped by using normed coordinates
order of the filter
1 K Tn Td
Qn (s) = (10) κ= ; τf = ; τd = ; p = Tm s (15)
(1 + Tn s)n Km Tm Tm
that may be used both in the disturbance reconstruction and which yields
as a prefilter in the setpoint filtration that prevents possible κe−τd p
Fs1 (p) = ((τn p+1)n −e−p )+κe−τd p
overshotings. (16)
Km κe−τd p ((τn p+1)n −e−p )
Fi1 (p) = κ((τn p+1)n −e−p )+κe−τd p
Definition 1 (Predictive I0 controller - P rI0 ): Under the
P rI0 controller we will understand a static feedforward control
with a gain 1/K0 extended by the input or output disturbance The robustness analysis by the PPM will be based on theo-
reconstruction by means of a ”predictive disturbance observer” rems similar to the Theorem 2 proven in [2] stating that when
(PDO) and its compensation at the static feedforward output considering PP including items IAE(ȳs ), IAE(ȳd ), T V 0 (ȳs ),
compensation (Figure 1) with the PDO filter time constant Tn T V 1 (ȳd ), T V 0 (ūs ) and T V 0 (ūd ) generated for the considered
and by the prefilter with the time constant Tp (in the simplest plant (5) with K = 1, Td = 1 over chosen grid of points
case with Tp = Tn ). Tm , n, Tn from the setpoint responses (r = 1, di = 0) with
output and input ȳs (τ ), ūs (τ ) and from the input disturbance
For Tp = Tn , the setpoint-to-output closed loop transfer responses (r = 0, di = 1) with ȳd (τ ), ūd (τ ) stored and
function of the PrI0 controller corresponding to the output y1 expressed over grid of normalized variables
and a plant-model parameter mismatch may be derived as
κ = K/Km ∈ [κP P min , κP P max ]
Ke−Td s τd = Td /Tm ∈ [τdP P min , τdP P max ] (17)
Fs1 (s) = (11)
Km ((Tn s + 1)n − e−Tm s ) + Ke−Td s τn = Tn /Tm ∈ [τnP P min , τnP P max ]
3 6
n=1
n=2 5.5
2.5 n=3
n=4
5
n=5
2
4.5
∆IAEΣ
IAEΣ
1.5 4
3.5
1
n=1
3
n=2
0.5 n=3
2.5 n=4
n=5
0 2
0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1
∆ =(T −T )/T ∆ =(T −T )/T
T dmax dmin dmin T dmax dmin dmin
Figure 2. Relative dead time robustness characteristics for = 0.001 and Figure 3. Absolute dead time robustness characteristics for = 0.001 and
different filter orders different filter orders
then, the loop performance measures corresponding to any trace out a line segment denoted as an Uncertainty Line
plant parameters K, Td and controller parameters Km , Tm , Tn Segment (ULS) that represents a special case of an uncertainty
belonging to the range (17) with the corresponding responses curve segment (UCS) [2]. It has been shown there that the aim
ys (t), us (t), yd (t) and ud (t) may be calculated according to of an optimal robust controller tuning satisfying given shape
IAE(ys ) = Td IAE(ȳs ) ; IAE(yd ) = KTd IAE(ȳd ) constraints may be interpreted as to optimally locate an UCS
(21) to lie in an area with admissible shape parameter values
T V0 (ys ) = T V 0 (ȳs ) ; T V1 (yd ) = K T V 1 (ȳd )
1 and, according to the Theorem 3 in [2], by locating an UCS
T V0 (us ) = K T V 0 (ūs ) ; T V1 (ud ) = T V 1 (ūd ) into areas with zero deviations from ideal shapes one gets
(18)
full performance invariance against dead time perturbations.
The aim is to keep for all possible values of an uncertain
B. Robustness characteristics for a dead time uncertainty parameter (19) and a chosen controller tuning the tolerable
In [14] it has been shown, how this controller may be tuned shape related deviations (3),(4) (9) yielding a minimal value
for n = 1 by the PPM by considering uncertainty of both of the cost function (2).
model parameters (6) in a 3D coordinate space. However, in
To simplify the evaluation of the performance changes, it
order to be able to evaluate impact of the filter order n on
will be characterized by the sum of the IAE values for the
the overall closed loop robustness, it is more appropriate to
setpoint and disturbance steps of the output y1 corresponding
derive new robustness characteristics expressing influence of a
to the dead time in the feedforward path
particular uncertain parameter (approximated by a fixed value
Km , or Tm ) on the loop performance. Thereby, we will require
to keep the shape related constraints (3) and (4) and to evaluate IAEΣ = IAEs + IAEd (22)
the corresponding increase of the IAE (2) values of the setpoint
and disturbance step responses. For such situations it is enough For an increasing ∆T (20), one gets by the PPM the deadtime
to generate 2D performance portraits corresponding either to robustness characteristics showing for particular filter orders n
κ = 1, or τd = const ≈ 1. the corresponding increase of the IAE (22) on ∆T (Figure 2).
It is to note that the minimal relative performance decrease
Thus, to characterize robustness aspects of the closed loop corresponds to n = 2. However, at the same time, the
control, the corresponding properties will now be checked by absolutely lowest performance decrease due to the dead time
the PPM [2], [3] for a single uncertain parameter. Firslty, the uncertainty corresponds to n = 1 (Figure 3). Thus, when
closed loop dead-time wishing to discuss about impact of the filter order on the
loop robustness, one has strictly define its meaning. Since
Td ∈ hTdmin , Tdmax i (19) the discussion of robustness have been mostly oriented on a
will be considered. Thereby, the dead time uncertainty may be robust stability [15], this aspect has been neglected up to now.
characterized by a relative uncertainty value In interpreting these characteristics by the PP (Figure 4) one
has to note that by increasing the dead time uncertainty (19)
∆T = (Tdmax − Tdmin )/Tdmin (20) the length of an uncertainty line segment (corresponding to
Then, in a normed 2D parameter subspace τd , τn , all possible all possible working points with the first coordinate Td /Tm
working points increases and thus, to be located in regions with the shape
related deviations below the chosen value (9) it must be shifted
κ = K/Km = 1 ; τd = Td /Tm ; τn = Tn /Tm (21) to higher Tn
TV0(us) TV1(ud) TV0(us) TV1(ud)
1 1
0.001
τn=Tn/Tm
τn=Tn/Tm
0.001 4 4
0.1 0.1
0. 1
0.00.1
1
0.5 1 0.5 1 0 1
0.00
0.00
0.1 0.101 2 0 2
10 10
1 0.1
1 1 1
0.6 0.8 1 0.6 0.8 1 0.5 1 1.5 2 0.5 1 1.5 2
TV0(ys) TV1(yd) TV0(ys) TV1(yd)
1 1
τn=Tn/Tm
τn=Tn/Tm
0.001 0.001 4 4
0.1 0.1
1 1
0.00.1 0.00.1
1
1
0.5 1 0.5 1
0.00
0.00
0.1 0.1 2 0 2 0
10 10
1 1 1 1
0.6 0.8 1 0.6 0.8 1 0.5 1 1.5 2 0.5 1 1.5 2
IAE(ys) IAE(yd) IAE(ys) IAE(yd)
1 2 1 2 3 3
τn=Tn/Tm
τn=Tn/Tm
4 4
4
4
3 3
8
0.5 0.5 6 2 6 2
4 1.5 4 1.5 2 4 3
2 4 3
65
18
30 65
18
30
23 23 2 1.5 234 2 1.5 234
0.6 0.8 1 0.6 0.8 1 0.5 1 1.5 2 0.5 1 1.5 2
τ =T /T τ =T /T κ =K/K κ =K/K
d d m d d m d m d m
Figure 4. The closed loop performance portrait for n = 1 and an uncertainty Figure 5. The closed loop performance portrait for n = 1 and an uncertainty
line segment (ULS) corresponding to an uncertain dead time (19), κ = 1 line segment (ULS) corresponding to an uncertain gain (23), τd = 1
τn≥τn,min=2
C. Robustness characteristics for a gain uncertainty 1.4
n=1
Let us now consider a plant gain uncertainty n=2
1.2
n=3
K ∈ hKmin , Kmax i (23) n=4
1 n=5
that may also be characterized by a relative value
∆K = (Ksmax − Ksmin )/Ksmin (24) 0.8
∆IAEΣ