Performance Portrait Method: An Intelligent PID Controller Design Based on a Database of Relevant Systems Behaviors

The article deals with a computer-supported design of optimal and robust proportional-integral-derivative controllers with two degrees of freedom (2DoF PID) for a double integrator plus dead-time (DIPDT) process model. The particular design steps are discussed in terms of intelligent use of all available information extracted from a database of control tracking and disturbance rejection step responses, assessed by means of speed and shape-related performance measures of the process input and output signals, and denoted as a performance portrait (PP). In the first step, the performance portrait method (PPM) is used as a verifier, for whether the pilot analytical design of the parallel 2DoF PID controller did not omit practically interesting settings and shows that the optimality analysis can easily be extended to the series 2DoF PID controller. This is important as an explicit observer of equivalent input disturbances based on steady-state input values of ultra-local DIPDT models, while the parallel PID controller, allowing faster transient responses, needs an additional low-pass filter when reconstructed equivalent disturbances are required. Next, the design efficiency and conciseness in analyzing the effects of different loop parameters on changing the optimal processes are illustrated by an iterative use of PPM, enabled by the visualization of the dependence between the closed-loop performance and the shapes of the control signals. The main contributions of the paper are the introduction of PPM as an intelligent method for controller tuning that mimics an expert with sufficient experience to select the most appropriate solution based on a database of known solutions. In doing so, the analysis in this paper reveals new, previously undiscovered dimensions of PID control design.


Introduction
The ever-increasing possibilities of computer technology allow an ever-increasing part of the processes related to the automatic control of systems to be carried out using programmable devices. Their design usually attempts to emulate important characteristics of human intelligence, leading to a discipline called intelligent control.
Intelligent control methods are successfully promoted, especially in areas where conventional control design is associated with the need to solve complex problems, such as, e.g., arising in nonlinear systems with incomplete and inadequate representation and under incomplete specifications of how to do this in an uncertain environment toward a not definitely specified control goal [4,5]. Their use also tends to reflect the personality of the designer, who prefers intuitive and heuristic design methods and returns to the analytical model-based approach only to verify and validate his intuitive results.
A large amount of work in the field of AI has already been devoted to optimizing the settings of PI and PID controllers. However, tied to the intrinsic limitations of this type of control, it has not always led to significant progress in the achieved performance. The reasons could be partially caused by the incomplete understanding of the functionality of these controllers. Just recently, they have been discussed as structures for system stabilization with an explicit or implicit disturbance observer (DOB) for the reconstruction of the input disturbance from steady-state input related to ultra-local integral models of controlled systems [27,28]. Although this type of DOB represented by a first-order lowpass filter is interesting due to its exceptional simplicity, it inherently places significant limitations on the dynamics of the reconstruction processes and on the acceptable settings of the controllers. In borderline, but not entirely rare, cases, these intrinsic limitations of PID control can finally lead to the abuse of the concept of intelligent control, when the adjective "intelligent" is rather a sign of the pursuit of a cheap and effective attribute. Ultimately, such approaches can lead to pejorative attitudes to misleading overly simplified "intelligent" designs.
The presented work brings a more consistent emphasis on two aspects of PID controller design-system modeling using ultra-local integral models and possible DOB-based PID interpretation-which belong to the basic indebtedness of traditional methods of automatic control and their innovative modifications used within intelligent control. Whereas the use of ultra-local models is today frequently preferred in such areas as model-free control (MFC) [29,30], or active disturbance rejection control (ADRC) [31,32], in terms of use in the field of PID control (despite frequent use, as in [33][34][35]), this is a neglected issue. Furthermore, the DOB-based view on PID control represents a relatively new and not yet sufficiently and systematically explored topic [27,28].
In addition, the paper develops another new dimension of an intelligent approach to the design of optimal and robust controllers based on the performance portrait method (PPM). After several initial steps [36,37], the PPM was fully developed, especially in the works [38,39]. It could be simply stated that this contribution represents an extension of the work dealing with optimal and robust control design of systems approximated by integrator plus dead-time (IPDT) models using 2DoF PI controllers [38] to the case of controlling DIPDT models by 2DoF PID controllers. However, such a formulation of the focus of the work is only very indicative because the control of DIPDT models brings a number of new moments that have influenced the scenario of work at many points.
With respect to this, the paper is structured as follows. A brief overview of the plant modeling from the point of view of positional control of autonomous vehicles is discussed in Section 2, together with the basic problems of an analytical 2DoF PD and PID controller design. Section 3 discusses the basic limitations of an analytical design of optimal 2DoF PID control and introduces the ideal shapes of transient responses, together with the related performance measures enabling us to evaluate deviations of achieved transients from their ideal wave-forms. Section 4 deals with the introduction of the performance portrait (PP) and its use in the performance portrait method (PPM), allowing optimal nominal parallel and series PID controller tuning. Contrary to [28], it offers an extension of PPM with the possibility of weighting the emphasis on the resulting behavior in terms of setpoint tracking and disturbance rejection while also considering the optimization of the reference setpoint signal weighting by optimal settings of pre-filter parameters.
Discussion of the main results of the paper in Section 5, together with their summary in the Conclusions, show that due to strong limitations of the up-to-now available mathematical apparatus, just the intelligent control based on the PPM is able to analyze and exploit all (even up to now hidden) dimensions of PID control.

PD and PID Controller Schemes for Analysis and Implementation
In order to be able to use illustrative physical interpretations of the considered steps when working with PPM and also to stress the differences with respect to the traditional approaches to PID control typical, e.g., in process control, let us assume the application of the method for controlling the position of vehicles with one degree of freedom of movement. Due to the existing analogies of the description of dynamic systems, the method can, of course, also be fully used for the control of any other dynamic systems, for the reliable stabilization of which the proportional-derivative (PD) controller must be used. As the main specific feature of the chosen domain of interpretation, we will highlight the fact that in terms of reference setpoint tracking, the need for smooth monotonic setpoint step responses without over-regulation will be emphasized.

Double-Integrator and Double-Integrator-Plus-Dead-Time Models
Due to the forces changing nonlinearly during a vehicle motion, whereby several system parameters are not fully known and can vary in time, control of a moving body can be treated by several possible approaches to nonlinear uncertain systems. All the possible approaches can be unified by application of Newton's second law when the vehicle acceleration is proportional to the resulting force applied (F = ma). When we leave the designations typical of physics in Figure 1 and instead denote the position of the vehicle by the variable x, its velocity asẋ and acceleration asẍ, and after introducing the gain K s = 1/m instead of the mass m, the system equation can be written in the form: Thereby, u r ≡ F corresponds to the system's manipulated variable (control signal), F r (x,ẋ) represents the resulting force acting against the control signal (consisting of friction and gravitational forces), and the input disturbance d i summarizes the influence of external or internal forces not taken into account by the considered model.
When considering repeated movements in a known environment and with a known load, it probably makes sense to apply nonlinear methods to the control design as exact linearization [40], or to consider linearization around the selected fixed operating point, replacing nonlinear dependencies F r (x,ẋ) with "local" linear approximations obtained, e.g., by the Taylor series expansion. However, when moving in an unknown environment, it is usually impossible to work on the control design with linear or nonlinear approximations F r (x,ẋ) depending on x andẋ. There is nothing left but to merge the unknown value F r (x 0 ,ẋ 0 ) in the current state x = (x 0 ,ẋ 0 ) t with the external disturbance d i into an equivalent disturbance: as typically used in ADRC.
In case of output disturbances d o , the output equation can be written as: As for the output disturbance, it can be, for example, an inaccuracy in the location of the vehicle or the drift of ships by sea currents. However, with regard to the inertia of the movement, the output disturbances of a jumping nature are certainly not possible, and therefore we will limit ourselves to considering the jump changes of the required position and the input disturbance (load) and consider y = x. In addition, when neglecting the possible control saturation by considering u(t) = u r (t) and taking into account that all other delays arising in the transmission and processing of information, in the actuators, filters, and sensors can be denoted as a model dead-time T m , the double integrator plus dead time (DIPDT) model [28,41] can be written as: In the case of a non-measurable output derivative, its value has to be reconstructed by a derivation combined with a binomial low-pass filter Q n (s): Thereby, the sampling period T s used for implementation must be chosen to fulfill T s << T f . The total loop dead-time T d must include, in addition to an estimate of the loop delay T m , an equivalent filter delay estimate T e : where T e (as proposed in [42]) approximates the filter dynamics in terms of a dead time that can simply be added to the total dead time T d , where the coefficient N can be given by values in the range from N = 0.5 (equivalence based on "the half rule") to N = 1 (equivalence based on "the average residence time").

Stabilizing PD Controller Tuning by the Triple Real Dominant Pole Method
Today, even people without an education in automatic control, but with the experience of a car driver, know that to stabilize the movement of the vehicle, its position, speed, load, and the driver's reaction time must be taken into account. Thus, to control the DIPDT system (4) with a constant input disturbance d i , a proportional-derivative (PD) controller expressed with the help of the Laplace transform as C s (s) = K p (1 + T D s) and supplemented by a disturbance feedforward counteracting the equivalent disturbances with the equation: has to be used, where y denotes the actual and w the reference setpoint output values. Here, d e represents a measurement, or estimate of the equivalent disturbance d e = d i + F r (x 0 ). Analysis of an optimal controller tuning for a DIPDT system by the triple real dominant pole method can be found in [41,43]. In [44], the PD controller has been proposed by modification of PI control for an IPDT system based on the relative time-delay margin. Furthermore, ref. [45] starts with the modification of PI control for an IPDT system, whereby the PD controller design is based on pole placement by using the Lambert W function.
DIPDT model stabilization by an ideal PD controller: having at its input the (tracking) control error: yields the closed loop transfer function: Its characteristic quasi-polynomial: has a triple real pole s o , when the loop parameters fulfill the requirements: i.e., when: Here, T o denotes the time constant associated with the real pole s o . As pointed out by one of the anonymous reviewers, according to [46], for the given closed-loop model, the maximum pole multiplicity is three. At the same time, to determine two unknown controller parameters and the position of the corresponding multiple pole, one has to solve a system (12) with minimally three solutions.
Although the analytical design of controllers of systems with delays and the transcendent nature of the mathematical description is associated with limitations, it will help to clarify several other possibilities of numerical solutions of the problem. The fastest nonoscillatory responses corresponding to a triple-real-dominant-pole (TRDP) s o = −0.5858/T d , which yields K p = 0.079/(K s T 2 d ); K d = 0.461/(K s T d ), can finally be shown to be guaranteed by dimensionless controller parameters: Due to the zero of F wy in (10), the corresponding setpoint step responses typically have overshoot. For its elimination, a pre-filter must be added to the PD controller (8): Its denominator cancels the numerator zero −1/T D . The numerator coefficient b can be determined to cancel one of the triple poles of the closed-loop s o = −1/T o .

From the PD Controller to the Series and Parallel PID Controller Design
In addition to Minorsky's work [1], which influenced the development of the simplest systems control theory and the development of relevant terminology by designing the three-term control concept, there was also a practical branch of research in which pneumaticbased controllers prevailed for a long time [47]. Although there has been significantly less information published in this area (in terms of patent protection), over time, elements of an intelligent approach that sought to emulate and replace the work of experts can also be traced. Thus, proportional-integration (PI) controllers were first created, in which the need to set the additional controller output signal (offset) required to achieve zero permanent control error under the influence of constant disturbances was replaced by the introduction of positive feedback from the controller output [27,28]. This created the controllers, originally called "automatic reset". Their essence can be explained by the observation that in integral systems, the steady state-output of the controller is equal to the negative value of the input disturbance d i . Because a counteracting signal must be added to the controller output to compensate for the constant d i acting on the plant input, a double inversion of the sign gives positive feedback. To measure the steady-state values of the controller output u r (see Figure 2), it was sufficient to take a given signal using a low-pass filter with a time constant T i substantially longer than the time constant of transients in the circuit with the stabilizing controller (denoted in previous section as T o ). Relatively short transients caused by input step changes affect the output of such a filter only negligibly, and the controller steady-state output values dominate in it. At this moment, it is necessary to stress that only when neglecting the constraints put on the control action in Figure 2 (only in the proportional control band), can the corresponding series PID controller be described by the transfer function: where the controller parameters considered with the indices "s" are: K c the controller gain, T i the integral, and T D the derivative time constant.
In such a situation, it is then also possible to get the parallel PID controller: with parameters denoted with the index "p". To get in both structures the same controller transfer function, it must hold: Whereas the re-calculation from a series to a parallel controller is always possible according to: the replacement of a parallel PID controller by a series one according to: is only possible for T ip ≥ 4T Dp .

Optimal Parallel PID Tuning by the Quadruple Real Dominant Pole Method
In [48], the PID controller design has been based on SIMC design [49,50]. In [33], it has been modified, based on optimization using the exact gradient method. In [44], the controller optimization was based on the dead-time margin.
Next, the analytical controller tuning will be based on the quadruple real dominant pole (QRDP) method [28,41,42]. If, for simplicity, we omit the index "p" from the parameter designation, application of the parallel PID controller (17) leads to the closed-loop transfer functions: Similar to above, from the characteristic quasi-polynomial: with four unknown parameters yields the requirement of a quadruple real dominant pole (QRDP) s o , which requires the conditions: to be met, the dimensionless parameters of the parallel PID controller: Overshooting of the setpoint step responses due to zeros of F wy in (21) can now be eliminated by a pre-filter: with the trivial numerator tuning The setpoint responses can be accelerated by cancelling a single dominant pole s o , and even more by cancelling two dominant poles s 2 o , when:

Basic Limitations of Analytical Tuning of PD and PID Controller
If we want to evaluate the basic pros and cons of the above analytical setup of PID controllers for the DIPDT system, we must first state the performance measures appropriate to the problem and its applications in (vehicle) control.

Evaluation of the Speed of Transient Response
For the tracking error (9) defined in the time domain as e(t) = w(t) − y(t), the speed of transients can be evaluated in terms of the integral of the absolute error (I AE): For the unit setpoint step responses and d i = d o = 0 it will be denoted as I AE w . For w(t) = 0, d o = 0 and unit input disturbance (load) step as I AE i . Some works (for example, [33,44]) used in derivation of Pareto-optimal (PO) PID controller with maximal sensitivity constraints in the cost function: where I AE o corresponds to unit step responses of the output disturbance d o . However, they do not give examples of applications with the possibility of step changes of d o . Of course, when applied to vessel control (as in [44]), output disturbances caused by drift due to sea currents can be expected. However, in no case will these be jump disturbances, and it is, therefore, questionable whether such optimization will lead to concise results. Similarly, when driving autonomous vehicles, output disturbances may occur due to inaccurate GPS location. However, such disturbances will be rather random in nature, and it will be more appropriate to classify them as measurement noise. In the following, therefore, the cost function: similar to [38,51] will be used to evaluate the optimality of control. Unit setpoint step responses lead under the TRDP PD controller (14) to optimal I AE values: When considering the parallel PID controller, depending on the prefilter (25) numerator order, the QRDP PID tuning (24) yields I AE values corresponding to unit setpoint and disturbance steps: Significant differences in the optimal values I AE w and I AE i of the controller's analytical design show why standardization using minimum achievable values I AE w,min and I AE i,min is necessary when optimizing general designs with an emphasis on setpoint tracking or disturbance rejection. They also document that, with the introduction of I-action, it is possible to compensate for the disturbances acting, but at the cost of prolonging the transients (in the case given by I AE w2 /I AE w by 5.513/4.1213 = 1.34 times, 1.92 times for I AE w1 , and by 2.5 times for I AE w0 ) with the current paradoxical increase in proportional gain (K c /K p = 0.1248/0.079 = 1.58 times).
When someone thinks about ways to improve the performance of the loops using a PID controller, it is also important to remember that: Detailed research based on the analysis of several other optimal settings of PI and PID controllers [27,28] has shown that this fundamental limitation of the disturbance reconstruction rate and thus the transient speed resulting from the nature of the disturbance observer used (reconstruction from steady-state controller output values) cannot be removed even using artificial intelligence methods.
Remark 2 (The inability to achieve QRDP optimal responses by a series PID controller). It follows from the transformation relations (20) for the conversion of the series controller parameters from the parallel PID parameters that due to: it is impossible to achieve QRDP optimal responses by a series PID controller with real parameters. Instead, optimal series PID controllers always have (see, for example, [33]): However, such a selection requires that their value be appropriately determined and the value of K cs added.

Evaluating the Excessive Controller Effort
In application of PID control, a high attention has to be devoted to measurement noise impact, which in combination with not appropriate controller tuning can lead to an unacceptable excessive controller effort. In [33,44,48], it has been evaluated in terms of total variation: Although TV(u) plays an important role in circuit assessment, the identification of this measure with excessive controller effort is incorrect because the total amount of absolute increments also includes useful control signal changes used to achieve the desired output change [50]. In other words, reaching the "ideal" value of TV(u) = 0 would mean that the output of the DIPDT system if it was initially in a steady state, would not change at all during the transient.
Although the value of TV (35) is related to excessive control efforts, it does not in itself have an immediate physical or mathematical interpretation. However, appropriate TV modifications can be useful in evaluating the deviations of the time responses from optimal shapes based on the concept of monotonicity. Monotonic output variable changes mean that all its increments have contributed to the output change from an initial value y 0 to the final value y ∞ without causing some counteractive movement. Thus, the output signal samples y i viewed with the sampling period T s yield a measure [38]: that can be interpreted as a deviation from monotonicity. TV 0 (y) > 0 signalizes a nonmonotonic transient with some increments not contributing to the requested change, e.g., for a smooth response with a single overshoot amplitude ∆, TV 0 (y) = 2∆. Therefore, TV 0 (y) finds application mainly in the evaluation of setpoint step responses at the system output. The situation will be more complicated when evaluating the output waveforms after a disturbance step. Due to the delayed response of the controller, which occurs as a result of to the transport delay of the system, the deviation of the output from the setpoint will increase in the first phase of the input disturbance step coming after steady state. When the response of the controller finally manifests itself, ideally the output of the system monotonically returns to the desired value: the output thus consists of two monotonic sections. Such a course is hereinafter referred to as one pulse (1P) shape. The deviation of the real course from its ideal 1P shape is calculated by summing the deviations from the monotonicity on these two individual sections according to: where y m / ∈ (y 0 , y ∞ ) lying outside the interval formed by the initial and final output values y 0 and y ∞ represents the point separating the monotonic sections of the 1P response.
Feldbaum's principle of time-optimal control of a vehicle described by the secondorder differential equation expressed the need for two intervals with constrained control signal alternating individual limits more than half a century ago [52]. It is true that in the literature focused on smoother PID controllers, Feldbaum's theorem has been practically forgotten (with exception of few authors, as in [53,54]). However, the requirement to account for active control pulses applies not only to relaying rectangular control interventions, but also to their smoother modifications, in which case some of the limit values may not even be reached. Especially in the field of vehicle motion control, it is certainly not possible to optimize their movement by minimizing the value of TV, but it must be adjusted for the values of useful interventions needed to accelerate and brake the vehicle. Thus, for a double integrator [37,38], the optimal shape of the control signal u(t), t ∈ [0 − , ∞) associated with an ideal monotonic setpoint changes and 1P disturbance responses of the output (given by the inverse plant dynamics) will be specified by two extreme points u m1 , u m2 occurring at times t 1 , t 2 ∈ (0, ∞), lying between the initial and final values u 0 and u ∞ and satisfying the condition: They form three monotonic control intervals of the two-pulse (2P) optimal control signal shape and the excessive control effort (summing the deviations from monotonicity over these three intervals) has to be evaluated according to: It means that for more complex controller output shapes (as, for example, in Figure 3), the total variations (35) has to be decreased by a "useful" contribution term combining the acceleration-deceleration amplitudes 2(u m1 − u m2 ) with the total control signal change u ∞ − u 0 considered with the polarity depending on the first amplitude value with respect to the final steady-state control signal.

Example 1: Performance Evaluation of PD and PID Controllers
To illustrate the properties of derived TRDP PD and QRDP PID controllers, consider their evaluation using the proposed performance measures and comparison with the design of the PD controller according to [45].
For the plant model (4), implementation filter (5), and the delay equivalence (6) with parameters: it will be considered: 1.
The PD controller according to [45] with parameters K p = 0.12705/(KsT 2 d ) = 0.0813; K d = 0.53228/(KsT d ) = 0.4258, the prefilter (15), T D = K d /K p = 5.2369 with the numerator tuning b = 0 (GJ). Simulation with the step T s = 0.001 yields the setpoint step responses in Figure 4 with the performance measures in Table 1.  (14) with the prefilter (15), b = T 0 (TRD Pb) and b = 0 (TRDP0) and the PD controller according to Gerov and Jovanovic [45] with the prefilter (15), b = 0 (GJ). In terms of I AE w (Table 1), the fastest transients apparently correspond to the TRDPb controller. However, due to the zero relative degree of the prefilter, there is a discontinuity of the control action variable at the beginning. The undesirably high peak in u can be removed by using a saturation block, without significantly slowing down the transient that occurs in the case of TRDP0. Similar problems would arise when ingesting TRDP PID with a zero relative degree prefilter. Thanks to the reconstruction of disturbances in the TRDP PID from steady-state values of the controller output, this controller also yields the slowest transients. However, all responses with controllers set by the multiple real pole method are characteristic with a minimum number of monotonic sections at the input and output of the system, which corresponds to the zero values of TV 0 (y w ) and TV 2 (u w ).
When using the controller GJ according to [45], the output response is already overregulated even when using the pre-filter (15) with b = 0, which is signaled by the non-zero value TV 0 (y w ). The output overshooting ∆y = 0.021 < TV 0 (y w )/2 then also increases the overall controller activity signalized by the TV 2 (u w ) > 0 value.

Performance Portrait Method
Examining all possible meaningful solutions to a problem before deciding which one best meets the defined criteria is one of the basic aspects of "intelligent" problem solving. It is the basis of a method called brainstorming, but also of such legally established procedures as public procurement. Verification of assumed meaningful solutions is also an approach that is very popular in the design of optimal control by practitioners who prefer experimentation over a theoretically justified design. They return to the theoretical arguments only when they already know the optimal solution and use the theoretical justification to support its adoption.
The performance portrait method (PPM) [28,[36][37][38] assumes verification of loop properties at all possible relevant controller settings. For the DIPDT models, such meaningful intervals of optimal controller parameters were analytically investigated by selecting different groupings of real dominant poles in [28]. After choosing the appropriate quantization of particular dimensionless variables, which should be fine enough to capture all the important specifics of the problem under consideration, it enabled one to limit the control design to a grid of points and thus to spare the time and energy. The finer the quantization is choosen, the higher the total number of points at which one has to verify the properties of the system, either by simulation or real-time control. By concise approximation (modeling) of initial systems and work with dimensionless parameters (as in (24)), the applicability of information obtained by system analysis on a given grid can be extended to a significantly wider set of systems.
When analyzing closed loop properties at the controller parameters corresponding to the selected tuning point, several performance measures are evaluated, such as I AE w , I AE i , TV 0 (y w ), TV 1 (y i ), TV 2 (u w ), TV 2 (u i ), corresponding to the steps of the reference variable w and the load disturbance d i . This list is, of course, not exhaustive and may be extended or modified as necessary (for example, by the values of the sensitivity functions M s and M t preferred by many authors, maximum overshoot values, settling times, etc.). As we showed in the previous example, the amplitude of the maximum overshoot can also be estimated from the value of TV 0 (y w ). Since the values of M s and M t are important in relation to the limits of stability, but with higher demands on positioning accuracy and the shape of transients, they do not give as concise information as deviations from monotonicity and will not be used in this contribution.

Definition 1 (Performance Portrait (PP) and Performance Portrait Method (PPM)).
The information about the process is obtained, evaluated, and stored using computer technology. If necessary, it can be step-wise extended or reduced and reused for optimal and robust design of control circuits.
The performance portrait method (PPM) represents a digitization of controller settings inspired by the behavior of practitioners who first verify everything by experimenting. It is a reflection of four types of learning styles [55] in the field of controller design and by its nature falls within the trends of Industry 4.0.
Furthermore, if someone refers to the performance portrait method (PPM) as a "brute force" approach [34], it should be noted that traditional methods of optimizing controller parameters are much closer to that name. With PPM in general, there is no problem with convergence to the optimal solution, while the search for optimums can be repeated with minimal effort indefinitely, with the choice of new criteria of optimality (or using their various combinations).

Example 2: Generating PP of the Parallel PID Controller
To facilitate the possibility of repeating the next procedure to potential applicants, the generation of PP will be explained using simulation of the loop with a 2DoF parallel PID controller with an additional disturbance observer. Based on the experience of [28] and after adding another dimension needed to determine the appropriate prefilter coefficient b, it follows as a suitable choice: In order to be able to generate PP independently of the choice of the derivative action filter needed for the parallel PID controller implementation, a state-space scheme with available output derivative will be used ( Figure 5) with the prefilter: deleting the zero of F wy (s): The PP will be generated by simulating the unit setpoint and input disturbance step responses corresponding to the grid parameters (41). For each grid point, the performance measures I AE w , I AE i , TV 2 (u w ), TV 2 (u i ), TV 0 (y w ), TV 1 (y i ) will be stored in the corresponding performance measures matrices.
The simulations have been carried out in a total of 31 × 31 × 9 × 9 + 31 × 31 × 9 = 86,490 points. The duration of the PP calculation on a standard PC was 37 h. Because the calculation of performance measures at individual points is independent of each other, the whole problem can be advantageously parallelized. Of course, in the traditional evaluation of experiments by practitioners, it would not be realistic to work with such a high number of experiments and process characteristics using different performance measures. Thus, in terms of the detailed information considered in the controller design, PPM multiplies the capabilities of human experts.
In the next process of using PPM, we determine whether there is a point (or in the case of uncertain systems a compact group of points) in which a closed circuit would meet the specified limits on the values of dimensionless parameters. These must be chosen so that the requirements for the real system with specific values of the DIPDT model parameters are ultimately met.

Example 3: Optimal Nominal Tuning of the Parallel PID Controller
Choice of the "optimal" nominal solution will be based on weighted I AE values (30), allowing to take into account the weight on the results from the evaluation of unit setpoint and input disturbance step responses. As constraints in the optimization, the admissible shape-related constraints have been specified by the vector yw , yi , uw , ui as: Since the above analytical approach studied in Section Äl' yields ideal shapes of transient responses, it can be characterized by constraints: In practical applications, each of the permissible restrictions may be specified by a separate number. However, if we want to show general trends in this work regarding the impact of permissible restrictions, we will choose, for simplicity, a vector of "sufficiently" small figures: This choice is surely of an "ad hoc" nature, has a strong influence on the resulting transients, but represents a trade-off between practical usability and computational effort.
In a similar way as in [38], the following theorem can be derived for working with dimensionless quantities and quantities of a real circuit corresponding to the model (4) with T m = T d .

Theorem 1.
Let us consider PP, including items I AE(ȳ s ), I AE(ȳ i ), TV 0 (ȳ s ), TV 1 (ȳ i ), TV 2 (ū s ), and TV 2 (ū i ) generated for the IPDT plant with K s = 1, T d = 1 over the chosen grid of points K c , T i , T D , b by simulating setpoint responses (w = 1, d i = 0) with output and inputȳ s (τ), u s (τ) and input disturbance responses (w = 0, d i = 1) withȳ i (τ),ū i (τ). The properties of these responses are stored and expressed over a grid of dimensionless variables: Then, the PP items corresponding to any plant model parameters K s , T d , and prefilter tuning b belonging to the range (47) with the corresponding responses y s (t), u s (t), y i (t), and u i (t) can be calculated according to: This theorem shows that it is enough to carry out the optimal design just in the space of dimensionless parameters and to calculate the actual values of the performance measures just when wishing to show the relative character of some chosen admissible shape-related deviations (44)- (46), which have to be interpreted in the context of particular model parameters T d and K s .
For s w = 0.5 and s i = 1 − s w = 0.5, the search over the stored PP with performance measures satisfying the above requirements yields vectors of optimal controller tuning parameters: A comparison of transients with the QRDP analytical method in Figure 6 shows that PPM makes it possible to find faster transients. Thereby, the advantage of QRDP tuning (24) is that it works with lower controller gains (which is important in noisy loops) and, nominally, it gives ideal response shapes ( = 0) for any values of the model parameters K s and T d . The non-zero shape deviations specified for dimensionless quantities are subject to changes depending on the parameters of the model when moving to real values of the shape deviations, which means that they can ultimately be significantly different [28,38]. This transformation must be taken into account when specifying search constraints. The speed of transients can be increased by tolerating larger deviations from ideal courses (i.e., by increasing the values).
The location of the found optimal settings in the overall context of PP can be seen using the cross-sections in Figure 7. In the case of setpoint responses, it can be noticed that at higher values of b, the area corresponding to the transients with small deviations from the monotonicity disappears. Allowing larger deviations from ideal shapes of transients leads to acceleration (decreased I AE values) of both setpoint and disturbance responses.
Remark 3 (Impact of increased loop inertia on the ideal controller output shape). Allowable areas specified by TV 2 (u) ≤ 0.001 in Figure 7 are very small and difficult to display using contours. However, from the application of Feldbaum's theorem to a double integrator supplemented by the integral action of the controller and delays, which further increase the circuit order, the question arises whether it might not be interesting to look for solutions with higher allowable TV 2 (u) values or directly with higher number of allowable pulses of optimal system input, evaluated, for example, by TV 2 (u) [50]. Figure 6. Unit setpoint step responses (left) and unit input disturbance step responses (right) achieved with the QRDP PID controller (black), responses corresponding to POI-PID controllers [44] with a prescribed M s constraints and responses corresponding to optimal tuning calculated by the PPM with the parameter grid (41) under the performance specifications (46) and (44); T m = 1; K s = 1; T s = 0.001.

Parallel PID Controller Tuning Optimized under Consideration of Sensitivity Constraints
To get a picture of the differences in the approaches of traditional PID controller design, often based on I AE optimization supplemented by sensitivity constraints, the responses obtained using PPM and QRDP analytical tuning are compared with the optimal solutions from [44] given for three different M s values: κ = 0.0694; τ D = 5.7675; τ i = 13.3862; M s = 1.59 κ = 0.0974; τ D = 5.0837; τ i = 11.9603; M s = 1.80 κ = 0.1215; τ D = 4.6796; τ i = 11.2708; M s = 2.00 (50) To avoid overshooting of the setpoint step responses, the prefilter (42) with the numerator parameter b = 0 is used, whereas the output derivative signal is taken from the plant model. Figure 6 show these responses yet much slower than for the QRDP tuning. If the reason for slowing down the responses should be the argument that it is necessary to increase the robustness in case of uncertainty of the model, the answer is that in that case, it is possible to search for other controller tuning using PPM. However, it does not make sense every time, even if the parameters of the model are supposed to be known well enough.

Example 4: Optimal Nominal Tuning of the Series PID Controller
It is well known (see for example [33,44]) that higher dynamics can be achieved with parallel PID controllers than with series controllers. As mentioned in Remark 2, even a parallel QRDP PID controller cannot be converted to an equivalent series controller. However, the parallel controllers only contain the disturbance observer implicitly [27,28]. If we need to explicitly monitor the disturbance for some reason, the controller must be supplemented with an additional block of low-pass filter (as in Figure 5) with the time constant T i . From the point of view of disturbance reconstruction, but also with regard to the impact of control signal limitations, it is therefore more advantageous to work with a series PID controller, including an explicit disturbance observer.
Although the specification of the PP (41) calculation was designed with regard to the calculation of the parallel controller, the obtained PP can also be used to calculate the parameters of the series PID controller. At the same time, it can be expected that the fastest achievable transients will correspond to the borderline case in terms of parameter conversion (20), when: When we search in PP (41) using these constraints, we find that for s w = 0.5 and s i = 1 − s w = 0.5, it is not possible to find an optimal series PID controller that would provide transients with a deviation from ideal shapes ≤ 0.001.

Remark 4 (Influence of rounding errors).
We have already pointed out the small areas of permissible TV 2 (u) values in Remark 3. Due to κ and δ quantization in the τ D = δ/κ calculation and in further conversion steps, the situation will worsen further. The problems could be eliminated by generating another PP for the series PID design, which would work directly with κ, τ D , τ i , and b s corresponding to the prefilter: However, it is also possible to use the PP calculated for the parallel PID, but the values of the permissible values of TV 2 (u) need to be increased slightly. Figure 8, therefore, shows the optimal transients corresponding to the optimal series PID when selecting a modified vector of admissible shape deviations: = yw = yi = uw = ui ∈ {0.002, 0.01, 0.1} (53) Then, the search over the stored PP with performance measures satisfying the above requirements yields vectors of optimal controller tuning parameters (the indices "p" and "s" are used to distinguish the primarily found dimensionless values of the parallel controller and the recalculated values of the serial PID): = 0.002 : κ p = 0.13; δ p = 0.57; τ ip = 18; κ s = 0.065; τ Ds = 9; τ is = 9; b p = 0.7; = 0.01 : κ p = 0.15; δ p = 0.60; τ ip = 16; κ s = 0.075; τ Ds = 8; τ is = 8; b p = 0.7; = 0.1 : κ p = 0.20; δ p = 0.69; τ ip = 14; κ s = 0.100; τ Ds = 7; τ is = 7; b p = 0.7.
The responses were obtained using a series PID scheme with the output derivative available in Figure 9.
A detailed look at the location of the corresponding operating points in the PP (see Figure 10) shows that, in terms of the size of the dimensionless gains, influencing strongly the speed of transients, the limits on the permissible value of TV 2 (u) are decisive.
If, in order to accelerate the transients, we keep the allowable shape limits of the output yw and yi according to (53) and increase the allowable input (controller output) restrictions in all three cases to the value: Another way to accelerate transients would be to accept transients with a higher number of monotonic control intervals at the controller output. Figure 8. Unit setpoint step responses (left) and unit input disturbance step responses (right) achieved with the parallel QRDP PID controller (black) and responses corresponding to series PID controller (see Figure 9) with the optimal tuning calculated by the PPM over the parameter grid (41) under the performance specifications (53) and (44); T m = 1; K s = 1; T s = 0.001. The resulting waveforms of the individual circuit responses are presented in Figure 11. Even without numerical quantification of the obtained responses, there is an evident significant increase in speed, especially in terms of compensation of disturbance steps, but also in terms of the speed of reconstruction of the input disturbance. Since this can already be noticed in the series PID controller in Figure 9 compared to the reactions of the parallel PID controller Figure 5, it can be stated that while the parallel PID controller allows faster transients in terms of output, a series PID controller with an explicit disturbance observer is faster in terms of disturbance reconstruction.

Remark 5 (Limitations on T i selection ).
Understanding PID functionality, which is based on disturbance reconstruction from the controller output steady-state values, is crucial to understanding why it is not possible to arbitrarily reduce the integral time constant T i with respect to the transients of the loop stabilised by the PD controller. The existence of T i size limitations, mentioned already in Remark 1, is obvious from all the above analytical and numerical calculations of the optimal PID controller settings. Therefore, countless attempts to overcome this handicap by applying AI methods cannot be considered intelligent.

Example 5: Considering the Pareto-Front in Setpoint Tracking and Disturbance Rejection
It is well known that a PID controller designed with the best setpoint tracking in mind may not provide optimal disturbance rejection and vice versa [33,38,51]. Depending on the specific requirements of the considered application, it may therefore be interesting to look for the optimal setting of the controller at variable values of the weight parameters s w and s i = 1 − s w of the cost function (30). Figure 11. Unit setpoint step responses (left) and unit input disturbance step responses (right) achieved with the parallel QRDP PID controller (black) and responses corresponding to series PID controller (see Figure 9) with the optimal tuning calculated by the PPM over the parameter grid (41) under the output performance specifications (53) and (44) and the input admissible deviations (55); T m = 1; K s = 1; T s = 0.001. The admissible shape deviations will again be chosen as a mix of the allowable shape limits of the output yw = yi = 0.001 and the allowable input deviations (55) with uw = ui = 0.3. PP search with s w = 1 and s w = 0, whereby s i = 1 − s w , then provides controller parameters in Table 2. Table 2. Parameters of parallel (pPID) and series PID controllers (sPID) found for yw = yi = 0.001 and uw = ui = 0.3 in PP (41); κ p = K s K cp T 2 d , δ p = K s K dp T d , The setpoint and disturbance step responses corresponding to the specified parallel and series PID controllers are in Figure 12. As expected, with both controllers, setpoint step responses with s w = 1 are faster than with s w = 0. Here, too, the input disturbance is reconstructed more quickly in the case of series PIDs. However, it may come as a surprise that the optimal parallel PID with s w = 1 has a significantly worse disturbance response than both series PID controllers. Figure 12. Unit setpoint step responses (left) and unit input disturbance step responses (right) achieved with the parallel (p) and series (s) PID controllers (black) with the optimal tuning calculated by the PPM over the parameter grid (41) under the output performance specifications (55) and (44) and the cost function weights s w = 1 and s w = 0, s i = 1 − s w ; T m = 1; K s = 1; T s = 0.001.

Discussion
In terms of the need for the laborious initial generation of PP, required for its further use in the design of controllers, PPM represents an approach with high initial computational demands. However, the control design is, by nature, an interactive activity in which we frequently step-by-step strive to achieve the best possible performance by taking into account the often conflicting requirements arising from an infinitely diverse range of demands for real applications. Furthermore, with regard to very simple optimization options based on searching the stored properties of the waveforms corresponding to the individual settings, this is extremely advantageous with the use of PPM. Moreover, this method is also advantageous in terms of the support provided in improving the features of interactive designs (for example, it shows which constraints are dominant in the effort to improve dynamics). In terms of time savings compared to traditional optimization, whose processes must be repeated for each change in the cost function and given constraints, it is significantly more advantageous for interactive PPM design.
We have mentioned several times that the DOB contained in PID controllers is based on the properties of ultra-local integral models of the system. However, we will return to their importance once again. In the design, we try to approximate the practically unlimited range of various dynamical systems using suitably selected mathematical models, for which we implement the design of optional parameters of the controller. From the point of view of the universal applicability of the results achieved by means of individual models, ultra-local integral models [38,56] have proved to be important since the first known method [57] for setting controllers. Although over time it has become possible to achieve a significant improvement in the dynamic properties of the circuit with more complex structures (as, for example, [58]), due to their simplicity, PID controllers still maintain an exceptional place in practice. As mentioned in Section 2, in [33][34][35]44], and in many other ADRC [31,32] or MFC publications [29,30], the use of ultra-local models to approximate different dynamical processes is an important robust control tool in which the system's internal feedback is included in an equivalent disturbance compensated by an integral action of the controller. Their optimal setting for the control of DIPDT models using PPM thus naturally complements the work [38,39] devoted to the optimal and robust setting of PI controllers for IPDT models.
The digitization of the automatic control design, which is the mainstay of the Industry 4.0-5.0 campaigns, is of course also reflected in the design of the controllers themselves. From this point of view, it is far more advantageous to work with various computer-support systems of the PPM type than with traditional text overviews of optimal settings of [59] controllers. Although analytical methods as [33,49] provide a free parameter to take into account application requirements, the impact of that parameter may not be predictable enough to allow effective design. The optimization method with sensitivity constraints are, in general, not effectively and consistently applicable to unstable systems.
Regarding the development trends in the field of automatic control, we should note the increase in open solutions, in which the designer does not have to rely on a relatively narrow offer of industrially manufactured controllers with a relatively fixed structure, but implements its own solution using a far more universal base of programmable embedded controllers. This is also related to the development of ever new regulators and structures [58,60], which significantly push the boundaries of achievable circuit characteristics. It follows from this statement that the PPM will also need to be developed for the design of additional controllers with more advanced and efficient disturbance observers.

Conclusions
The paper discussed the optimal design of series and parallel PID controllers for a double integrator plus dead-time process using the performance portrait method. The proposed design imitates an expert who selects the most suitable solution based on appropriate criteria (performance measures) from a database of known solutions. The database of relevant responses, referred to as the "performance portrait", is evaluated in terms of the weighted sum of the IAE criterion in closed-loop tracking and disturbance rejection, and the shape of the input and output signals of the process (compared to their ideal shapes). This new intelligent control method has no convergence problems and allows interactive multi-criteria optimization of the control according to the required performance and robustness.  Acknowledgments: The first author would like to thank H. Hoyer and M. Gerke for creating excellent working conditions during the Mercator professorship at the Fern-Universität in Hagen, which led to the discovery of the performance portrait method. This work was also supported by E-Academia Slovaca, a non-profit organization, Sadmelijská 1, 831 06 Bratislava, Slovakia.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript: