Combination of Data-Driven Active Disturbance Rejection and Takagi-Sugeno Fuzzy Control with Experimental Validation on Tower Crane Systems

In this paper a second-order data-driven Active Disturbance Rejection Control (ADRC) is merged with a proportional-derivative Takagi-Sugeno Fuzzy (PDTSF) logic controller, resulting in two new control structures referred to as second-order data-driven Active Disturbance Rejection Control combined with Proportional-Derivative Takagi-Sugeno Fuzzy Control (ADRC–PDTSFC). The data-driven ADRC–PDTSFC structure was compared with a data-driven ADRC structure and the control system structures were validated by real-time experiments on a nonlinear Multi Input-Multi Output tower crane system (TCS) laboratory equipment, where the cart position and the arm angular position of TCS were controlled using two Single Input-Single Output control system structures running in parallel. The parameters of the data-driven algorithms were tuned in a model-based way using a metaheuristic algorithm in order to improve the efficiency of energy consumption.


Introduction
In this paper the authors propose to merge two techniques, the second-order data-driven Active Disturbance Rejection Control (ADRC) and the proportional-derivative Takagi-Sugeno Fuzzy (PDTSF) logic control in two ways resulting in two control structures referred to as second-order data-driven Active Disturbance Rejection Control combined with Proportional-Derivative Takagi-Sugeno Fuzzy Control (ADRC-PDTSFC).The purpose of these new combinations is to improve the overall performance of control system (CS) structures.
Two ADRC-PDTSFC structures are proposed in the next sections.For the first ADRC-PDTSFC structure (ADRC-PDTSFC1) the proportional-derivative (PD) component in data-driven ADRC is fuzzified using a PDTSF logic controller.For the second ADRC-PDTSFC structure (ADRC-PDTSFC2) a PDTSF logic controller is added to the ADRC algorithm.To determine if the proposed ADRC-PDTSFC structures perform better than the ADRC structure, several comparison case studies are considered, where the performance of the CS with ADRC algorithm is compared with the performance of the CS with ADRC-PDTSFC1 and ADRC-PDTSFC2.To have a fair comparison between these CSs, the parameters of ADRC, ADRC-PDTSFC1 and ADRC-PDTSFC2 are tuned in a model-based way using the Grey Wolf Optimizer (GWO) algorithm, which solves an optimization problem characterized by an objective function expressed by the weighted sum of squared control errors and squared control signals.This formulation of the objective functions ensures a reduced energy consumption and thus energy efficiency.The experimental validation is conducted on a nonlinear Multi Input-Multi Output (MIMO) three-degrees-of-freedom tower crane system (TCS).Three values of the penalty factor in the weighted objective function are considered, resulting in three sets of experimental scenarios.
ADRC [1][2][3][4] is one of the most popular data-driven techniques along with Model-Free Adaptive Control [5][6][7][8][9], Model-Free Control [10][11][12][13][14] or Virtual Reference Feedback Tuning [15][16][17][18].The main advantage that made data-driven techniques so popular is that they use only the input/output data from the process.This advantage of data-driven control (in the tuning) can be very useful when the accurate mathematical model of the process is impossible or very hard to be determined or when the plant has strong nonlinearities.
One of the most successful and commonly used control techniques is fuzzy logic.Some TCS applications with fuzzy logic control are H ∞ based on adaptive fuzzy control technique [29,30], anti-swing controller designed using a time-delayed feedback of the load swing angle and an anti-swing fuzzy logic controller [31], gain-scheduling anti-swing controller that employs fuzzy clustering techniques [32], fuzzy logic for selecting the best crane [33], anti-swing combined with fuzzy control [34], fuzzy logic with sensorless payload deflection feedback [35] and Mamdani fuzzy logic controller [36].
The literature reports several combinations of ADRC and fuzzy logic, namely experimental validation on a brushless DC motor in [37], simulation results for a ship course in [38], simulation results for a permanent magnet synchronous motor in [39], permanent magnet synchronous motor sensorless control system for electric vehicle validated by simulation results in [40], simulation results for a flying boat in [41], simulated and experimental validation for permanent magnet synchronous motor servo system in [42], simulated and experimental validation of three motor synchronous system in [43], simulated validation of attitude maneuver control of a spacecraft in [44], simulated validation of insulin infusion in type 1 diabetic patients in [45] or electrical actuation system validated through simulation results in [46].
The TCS nonlinear laboratory equipment is very practical because it is actually a small-scale representation of an industrial tower crane used in civil engineering.However, for the moment the state of the art is very poor when discussing about controllers designed for TCS.Some TCS applications with different types of controllers include simulation and experimental results using a tensor product model transformation-based nonlinear feedback controller, with an additional neural network-based friction compensator [47], simulation results using recurrent neural networks evolved by evolutionary algorithms [48], simulation and experimental results using gain scheduling feedback control [49], simulation and experimental results using nonlinear model predictive control [50], experimental results with input-shaping control to suppress double-pendulum payload oscillations [51], experimental results with adaptive control [52], experimental results for a nonlinear anti-swing feedback controller with an integral term [53], experimental results with hierarchical artificial neural network-based adaptive fuzzy logic [54].
The new contributions of this paper are:  Optimal tuning of ADRC, ADRC-PDTSFC1 and ADRC-PDTSFC2 parameters in terms of using more efficient the energy consumption.
The design of the proposed controllers in the following sections is accomplished by following the steps:

•
In the first step an extended Luenberger state observer was developed to estimate the not directly observable effects of the neglected dynamics as a disturbance.By the use of this estimation an ADRC controller was developed within the frames of the LTI systems with proportional and derivative feedback terms in discrete time approximation.

•
In the next step the ADRC controller was further improved by fuzzifying the PD terms in the first version ("ADRC-PDTSFC1 structure") in a Two Inputs-Single Output Fuzzy Controller with a relatively simple rule-base using 2 × 3 fuzzy membership functions.

•
The second version ("ADRC-PDTSFC2 structure") is a modification of the same initial ADRC controller, too, by adding further PD feedback terms to the original ones before fuzzifying them.This latter structure has more independent parameters than the previous one.

•
In the case of both versions the so-obtained "control models" (containing the parameters of the membership functions) were inserted into an optimal controller in which the sum of the squares of the tracking errors plus the sum of the squares of the control signals weighted with a "penalty factor" were minimized for three different penalty factors: (a) a big value that made the precise trajectory tracking impossible; (b) a rather reasonable value that allowed acceptable tracking by consuming up more energy by the controller; and (c) a small value that only slightly prohibited the energy consumption of the control process.

•
In the last step experimental investigation of the operation of these controllers was carried out.
The paper is organized as follows: a short overview of second-order data-driven ADRC is presented in the next section.Section 3 offers the design approach of ADRC-PDTSFC structures.Section 4 describes the nonlinear TCS laboratory equipment used for experimental validation.The real time experiments that validate the ADRC, ADRC-PDTSFC1 and ADRC-PDTSFC2 structures are given in Section 5.The conclusions with the advantages and the disadvantages of the proposed structure are highlighted in Section 6.

Data-Driven Second-Order ADRC Structure
The design of the second order ADRC algorithm starts with the second-order Single Input-Single Output (SISO) process model [1][2][3]: where τ ∈ R is the natural period of oscillation of the process, ω ∈ R is the damping factor of the process, κ ∈ R is the process gain, u(t) ∈ R is the continuous-time control signal, y(t) ∈ R is the continuous-time controlled output and in η(t) ∈ R is the unknown disturbance that also includes the unmodeled dynamics of the process.
Using the notation a = κ τ 2 , a ∈ R and splitting it into a = a 0 + ∆a, with a 0 ∈ R as a known estimate and ∆a ∈ R as the unknown modelling error, the second-order SISO process in (1) is rewritten as: ..
Considering a generalized disturbance term f (t) ∈ R that may contain terms that convey inconvenient coupling with y(t), .y(t), .y(t) that are "hidden" by η(t), f (t) depends on the unknown disturbances and unmodeled dynamics of the process and the unknown modelling error, with the notation the process model resulted after substituting the generalized disturbance term in (3) into the second-order SISO process in (2) results as a disturbed double integrator: ..

y(t) and f (t).
The gain matrix L = [l 1 l 2 l 3 ] T ∈ R 3×1 (column matrix, where T indicates matrix transposition) will be used in order to provide the estimate of the state variables in (5).The dynamic model of ELSO is [1][2][3]: where A ∈ R 3×3 , B ∈ R 3×1 , C ∈ R 1×3 are ELSO matrices.Substituting x1 (t) with ŷ(t), (6) becomes: The elements of gain L can be determined by arbitrarily choosing the eigenvalues of the matrix (A − LC) in the left half-pane.The discrete-time model of ELSO accounting for the zero-order hold element is The discrete-time second-order control law specific to ADRC is: where ∆y * (k) = y * (k) − y * (k − 1), y * (k) ∈ R is the desired reference trajectory, ê(k) ∈ R is the approximate tracking error, ∆ê(k) = ê(k) − ê(k − 1), and K 1 and K 2 are the parameters of the discrete-time second-order ADRC control law.Equation ( 9) indicates that we deal with a modified PD controller, so K 1 and K 2 can also be seen as the parameters of a PD component in the control law [1][2][3].
As shown in [55], the stability is guaranteed for a second-order system if there exists a Lyapunov function that satisfies an inequality, where the partial derivatives of the Lyapunov function is smaller than a user-chosen parameter.Then if under a disturbance another inequality is satisfied, the disturbed system will still be asymptotically stable.
The control system structure with second-order discrete-time data-driven ADRC algorithm is presented in Figure 1.where , and 1 K and 2 K are the parameters of the discrete-time second-order ADRC control law.Equation ( 9) indicates that we deal with a modified PD controller, so 1 K and 2 K can also be seen as the parameters of a PD component in the control law [1][2][3].
As shown in [55], the stability is guaranteed for a second-order system if there exists a Lyapunov function that satisfies an inequality, where the partial derivatives of the Lyapunov function is smaller than a user-chosen parameter.Then if under a disturbance another inequality is satisfied, the disturbed system will still be asymptotically stable.
The control system structure with second-order discrete-time data-driven ADRC algorithm is presented in Figure 1.

Second-Order Data-Driven ADRC-PDTSFC1 Structure
The discrete-time second order control law specific to data-driven ADRC-PDTSFC1 is a combination of the discrete-time second order data-driven ADRC control law in ( 9) and a PD Takagi-Sugeno Fuzzy logic controller.The PD component in ADRC-PDTSFC1 structure is fuzzified using a Takagi-Sugeno Fuzzy logic controller leading to the following control law: where the component is replaced by the PDTSFC1 component, and it is designed using the recurrent equation of the PD component in (9) [56]: This PDTSFC1 component is built around a Two Inputs-Single Output Fuzzy Controller (TISO-FC) as shown in Figure 2, where e B ˆ and e B Δ are the parameters of the membership functions, and the SUM and PROD operators are used in the inference engine of the fuzzy component [56].The rule base consists of nine rules presented in Table 1, with the rule consequents:

Second-Order Data-Driven ADRC-PDTSFC1 Structure
The discrete-time second order control law specific to data-driven ADRC-PDTSFC1 is a combination of the discrete-time second order data-driven ADRC control law in ( 9) and a PD Takagi-Sugeno Fuzzy logic controller.The PD component in ADRC-PDTSFC1 structure is fuzzified using a Takagi-Sugeno Fuzzy logic controller leading to the following control law: where the component χ(k) is replaced by the PDTSFC1 component, and it is designed using the recurrent equation of the PD component in ( 9) [56]: Energies 2019, 12, 1548 This PDTSFC1 component is built around a Two Inputs-Single Output Fuzzy Controller (TISO-FC) as shown in Figure 2, where B ê and B ∆ê are the parameters of the membership functions, and the SUM and PROD operators are used in the inference engine of the fuzzy component [56].The rule base consists of nine rules presented in Table 1, with the rule consequents: where the component is replaced by the PDTSFC1 component, and it is designed using the recurrent equation of the PD component in ( 9) [56]: This PDTSFC1 component is built around a Two Inputs-Single Output Fuzzy Controller (TISO-FC) as shown in Figure 2, where e B ˆ and e B Δ are the parameters of the membership functions, and the SUM and PROD operators are used in the inference engine of the fuzzy component [56].The rule base consists of nine rules presented in Table 1, with the rule consequents:  Table 1.The rule base of the ADRC-PDTSFC controller.
Table 1, ( 11) and ( 12) lead to the conclusion that the PDTSFC1 component is practically a nonlinear combination of three discrete-time PD controllers placed in the rule consequents, which change according to the fired rules.The main purpose of the parameter γ i is used to adjust the overshoot of the CS with the ADRC-PDTSFC1 technique for TCS control [56].
The control system structure with second-order discrete-time ADRC-PDTSFC1 algorithm is presented in Figure 3.
Table 1.The rule base of the ADRC-PDTSFC controller.11) and ( 12) lead to the conclusion that the PDTSFC1 component is practically a nonlinear combination of three discrete-time PD controllers placed in the rule consequents, which change according to the fired rules.The main purpose of the parameter i γ is used to adjust the overshoot of the CS with the ADRC-PDTSFC1 technique for TCS control [56].
The control system structure with second-order discrete-time ADRC-PDTSFC1 algorithm is presented in Figure 3.

Second-Order Data-Driven ADRC-PDTSFC2 Structure
The discrete-time second-order control law specific to data-driven ADRC-PDTSFC2 is also a combination of the discrete-time second order ADRC control law in (9) and a PDTSF logic controller.A PD component is added to the control law as follows: ), ( ) ( ) ( and it is next fuzzified using a Takagi-Sugeno fuzzy logic controller and added to the ADRC law in (9) leading to the following control law: where the component is replaced by the PDTSFC2 component, resulting in the combined ADRC-PDTSFC2 control law.As in the previous case, the PDTSFC2 component is built around a TISO-FC respecting the representation given in Figure 2. Defining again the parameters e B ˆ and e B Δ of the input membership functions, using the same inference engine of the fuzzy component that makes use of the SUM and PROD operators [56], the rule base is the same, it consists of nine rules presented in Table 1 with i γ in ( 12) used to adjust the overshoot of the CS with the ADRC-PDTSFC2 technique for TCS control [56].
The control system structure with second-order discrete-time ADRC-PDTSFC2 algorithm is

Second-Order Data-Driven ADRC-PDTSFC2 Structure
The discrete-time second-order control law specific to data-driven ADRC-PDTSFC2 is also a combination of the discrete-time second order ADRC control law in (9) and a PDTSF logic controller.A PD component is added to the control law as follows: and it is next fuzzified using a Takagi-Sugeno fuzzy logic controller and added to the ADRC law in (9) leading to the following control law: where the component χ(k) is replaced by the PDTSFC2 component, resulting in the combined ADRC-PDTSFC2 control law.
As in the previous case, the PDTSFC2 component is built around a TISO-FC respecting the representation given in Figure 2. Defining again the parameters B ê and B ∆ê of the input membership functions, using the same inference engine of the fuzzy component that makes use of the SUM and PROD operators [56], the rule base is the same, it consists of nine rules presented in Table 1 with γ i in ( 12) used to adjust the overshoot of the CS with the ADRC-PDTSFC2 technique for TCS control [56].
The control system structure with second-order discrete-time ADRC-PDTSFC2 algorithm is presented in Figure 4.

The Nonlinear TCS Laboratory Equipment
All ADRC, ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers are validated experimentally on TCS, where the cart position and the arm angular position are controlled using two SISO loops running in parallel.All experiments are conducted starting with zero initial conditions.The nonlinear state-space model that describes the behavior of the TCS laboratory equipment is [ x y where ω x is the cart position, θ is the arm angular position, 8067 .9 = g is the gravitational , α and β are the angles that describe the payload position in x-z plane, L is the lift's line length , parameter L is not the same with gain matrix , and the outputs are the cart position

The Nonlinear TCS Laboratory Equipment
All ADRC, ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers are validated experimentally on TCS, where the cart position and the arm angular position are controlled using two SISO loops running in parallel.All experiments are conducted starting with zero initial conditions.The nonlinear state-space model that describes the behavior of the TCS laboratory equipment is [57]: . x 1 = x 5 , .
where x ω is the cart position, θ is the arm angular position, g = 9.8067 is the gravitational acceleration, x 1 = β, x 2 = α, α and β are the angles that describe the payload position in x-z plane, L is the lift's line length L = 0.8, parameter L is not the same with gain matrix L ∈ R 3×1 defined in Section 2, u c and u θ are the control signals for cart position control and arm angular position control, respectively, that is, the Pulse Width Modulation (PWM) duty cycle of the cart position and arm angular position direct current (DC) motors, the control limit values of u c and u θ are limited by the PWM duty cycle between [−1, 1], As Figure 5 shows, the communication between the personal computer and the TCS laboratory equipment is ensured through a power interface unit, a real-time data acquisition and a control universal bus called RT-DAC/USB2 board, where the PWM control logic is stored in a XILINX chip [57].The software pack provided by [57] called TowerCrane Toolbox also ensures the communication between the personal computer and the TCS laboratory equipment.
The block diagram of the MIMO control system structure that consists of two SISO control systems for separate cart position control and arm position control is given in Figure 6.As Figure 5 shows, the communication between the personal computer and the TCS laboratory equipment is ensured through a power interface unit, a real-time data acquisition and a control universal bus called RT-DAC/USB2 board, where the PWM control logic is stored in a XILINX chip [57].The software pack provided by [57] called TowerCrane Toolbox also ensures the communication between the personal computer and the TCS laboratory equipment.

Experimental Results
The block diagram of the MIMO control system structure that consists of two SISO control systems for separate cart position control and arm position control is given in Figure 6.As Figure 5 shows, the communication between the personal computer and the TCS laboratory equipment is ensured through a power interface unit, a real-time data acquisition and a control universal bus called RT-DAC/USB2 board, where the PWM control logic is stored in a XILINX chip [57].The software pack provided by [57] called TowerCrane Toolbox also ensures the communication between the personal computer and the TCS laboratory equipment.
The block diagram of the MIMO control system structure that consists of two SISO control systems for separate cart position control and arm position control is given in Figure 6.

Experimental Results
This section is dedicated to the presentation of three case studies and also experimental scenarios to validate the new controllers on TCS laboratory equipment in order to control the cart position and the arm angular position.In all experimental scenarios the parameters of the CS with ADRC controller, the CS with ADRC-PDTSFC1 controller and the CS with ADRC-PDTSFC2 controller are optimally tuned in a model-based way using a metaheuristic GWO algorithm [58][59][60] that solves the optimization problem: where Λ = [Λ c Λ θ ] T is a vector variable (column matrix) that contains the tunable parameters of the ADRC, ADRC-PDTSFC1 or ADRC-PDTSFC2 control laws, with the subscripts c and θ related to cart position and arm angular position, respectively, Λ is inserted as an additional argument to the variables in ( 16) in order to point out that all variables depend on the controller parameters gathered in Λ, Λ * = [Λ * c Λ * θ ] T is a vector variable (column matrix) that contains the optimal parameters of ADRC, ADRC-PDTSFC1 or ADRC-PDTSFC2 control laws, Λ contains the tunable parameters of ADRC control law expressed as: the tunable parameters of ADRC-PDTSFC1 control law are grouped in the matrices: the tunable parameters of ADRC-PDTSFC2 control law are grouped in the matrices: J e,u is the objective function, the sampling time is T s = 0.01 s, the duration of an experiment is set to 79 s, N = 7900 is the number of samples and it is the result of the division between the duration of the experiment and the sampling time, the control errors are the same with the ones illustrated in Figure 6, and λ is a penalty factor applied to the control law in order to reduce the energy consumption.The membership function parameters B êc , B ∆êc , B êθ , B ∆êθ related to ADRC-PDTSFC1 and ADRC-PDTSFC2 control laws are set by the designer such that all fuzzy rules are fired during the dynamic regime in all experiments conducted on TCS; this will contribute to the improvement of CS tracking performance.The energy consumption alleviation aspect is studied through three case studies, where the penalty parameter λ is set to 0.1, 0.01 and 0.001 in the first, second and third case study, respectively.The experiments in the three case studies are conducted using the same reference trajectory [y * c y * θ ] T :  ( In order to determine the parameters of ELSO for both cart position control and arm angular position control by ADRC, ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers, two up-chirp signals, one for each control signal, are applied to TCS in an open-loop experiment in order to determine whether ELSO outputs track the controlled outputs y c (k) and y θ (k).The parameters of ELSO, namely gains l 1 , l 2 , l 3 , are obtained by arbitrarily choosing a set of eigenvalues of (A − LC) in the half plane (to ensure the ELSO stability), after validating the ELSO modeled in (8) on the TCS equipment in the dynamic regime characterized by chirp control signals and setting the initial state vector of the ELSO to [0.45 0 0].We chose the eigenvalues {−0.11, −0.24, −0.69} for cart position and {−0.2, −0.17, −0.75} for arm angular position, and the following ELSO gains were obtained: A GWO algorithm with several details, including parameterization and vector operations given [58][59][60], is applied to solve the optimization problem in ( 16) using the nonlinear mathematical model of TCS in (15) and the dynamic regime characterized by the reference trajectory in (20), zero initial conditions and no additional disturbances, which is kept in all experiments in order to evaluate the objective function values.A specific feature of GWO algorithms is: where Λ A GWO ( ), Λ B GWO ( ), Λ Γ GWO ( ) are first three best vector solutions, also called alpha, beta, and delta, obtained after an iteration of GWO algorithm, is the iteration index, and max is the maximum number of iterations.The updated vector solution is: where i is the index of the current agent, namely the grey wolf, and the total number of agents is N GWO .The GWO algorithm is mapped onto the optimization problem defined in ( 16) in terms of: The parameters of the GWO algorithm that solves the optimization problem defined in ( 16) are N GWO = 20 and max = 100.They were chosen such that to reach an acceptable tradeoff to converge and computational resources.The search domains for all parameters of ADRC, ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers was set to the interval [−100, 100].
The first case study is characterized by setting the penalty factor to λ = 0.1.The GWO-based solving of the optimization problem in ( 16) leads to the following values of controller parameters: and the membership function parameters set to B êc = 0.1400, B ∆êc = 3.6616, B êθ = 0.1401, B ∆êθ = 0.0196 for the CS with ADRC-PDTSFC2 controller.
The experimental results for the three control systems are presented in Figure 7 in terms of the following plots: in Figure 7a) the control signal u c for ADRC is with red line, u c for ADRC-PDTSFC1 is with blue line and u c for ADRC-PDTSFC2 is with green line; in Figure 7b) the reference trajectory y * c is with thick black dashed line, the controlled output y c for ADRC is with red line, y c for ADRC-PDTSFC1 is with blue line and y c for ADRC-PDTSFC2 is with green line; in Figure 7c) the control signal u θ for ADRC is with red line, u θ for ADRC-PDTSFC1 is with blue line and u θ for ADRC-PDTSFC2 is with green line; in Figure 7d) the reference trajectory y * θ is with thick black dashed line, the controlled output y θ for ADRC is with red line, y θ for ADRC-PDTSFC1 is with blue line and y θ for ADRC-PDTSFC2 is with green line.
In order to have a fair comparison between the obtained results of the CSs with ADRC, ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers, and to mitigate the effects of random disturbances that occur unexpectedly during the experiments, the graphical results presented in Figure 7 were averaged ten times.with thick black dashed line, the controlled output θ y for ADRC is with red line, θ y for ADRC-PDTSFC1 is with blue line and θ y for ADRC-PDTSFC2 is with green line.
In order to have a fair comparison between the obtained results of the CSs with ADRC, ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers, and to mitigate the effects of random disturbances that occur unexpectedly during the experiments, the graphical results presented in Figure 7 were averaged ten times.2. .In this regard, considering the variance of the objective functions, the ranking of the objective functions in Table 2 is The obtained results representing the average and the variance of the objective functions J e,u of CSs with ADRC, ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers are presented in Table 2. Analysis of the results given in Figure 7 and Table 2 leads to the conclusion that the CSs with ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers have performance in the same range and both outperform the CS with ADRC controller.The following relationships between the average values of the objective functions of the three CSs are expressed: J e,u ADRC 1.127J e,u ADRC−PDTSFC1 and J e,u ADRC 1.127J e,u ADRC−PDTSFC2 .In this regard, considering the variance of the objective functions, the ranking of the objective functions in Table 2 is J e,u ADRC−PDTSFC2 < J e,u ADRC−PDTSFC1 < J e,u ADRC .
The second case is characterized by modifying the penalty factor to λ = 0.01.After executing the GWO in order to solve the optimization problem in (16), the following values of controller parameters were obtained: and the membership function parameters set to B êc = 0.1400, B ∆êc = 3.1161, B êθ = 0.1400, B ∆êθ = 0.0123 for the CS with ADRC-PDTSFC2 controller.
As in the previous case study, the experimental results for the control systems, namely ADRC, ADRC-PDTSFC1 and ADRC-PDTSFC2, are presented in Figure 8 in terms of the following plots: in Figure 8a) the control signal u c for ADRC is with red line, u c for ADRC-PDTSFC1 is with blue line and u c for ADRC-PDTSFC2 is with green line; in Figure 8b) the reference trajectory y * c is with thick black dashed line, the controlled output y c for ADRC is with red line, y c for ADRC-PDTSFC1 is with blue line and y c for ADRC-PDTSFC2 is with green line; in Figure 8c) the control signal u θ for ADRC is with red line, u θ for ADRC-PDTSFC1 is with blue line and u θ for ADRC-PDTSFC2 is with green line; in Figure 8d) the reference trajectory y * θ is with thick black dashed line, the controlled output y θ for ADRC is with red line, y θ for ADRC-PDTSFC1 is with blue line and y θ for ADRC-PDTSFC2 is with green line.for the CS with ADRC-PDTSFC2 controller.
As in the previous case study, the experimental results for the control systems, namely ADRC, ADRC-PDTSFC1 and ADRC-PDTSFC2, are presented in Figure 8 in terms of the following plots: in Figure 8   In order to have a fair comparison between the results obtained of the CSs with ADRC, ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers, and to mitigate the effects of random disturbances that occur unexpectedly during the experiments, the graphical results presented in Figure 8 were averaged ten times to eliminate this type of disturbance.
The results obtained that represent the average and the variance of the objective functions u e J , of CS with ADRC, ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers are presented in Table 3.In order to have a fair comparison between the results obtained of the CSs with ADRC, ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers, and to mitigate the effects of random disturbances that occur unexpectedly during the experiments, the graphical results presented in Figure 8 were averaged ten times to eliminate this type of disturbance.
The results obtained that represent the average and the variance of the objective functions J e,u of CS with ADRC, ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers are presented in Table 3.As in the first case study, the analysis of the results given in Figure 8 and Table 3 leads to the conclusion that the CSs with ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers have similar performances and they overperform the CS with ADRC controller.The following relationships between the average values of the objective functions of the three CSs are expressed: J e,u ADRC 1.385J e,u ADRC−PDTSFC1 and J e,u ADRC 1.437 e,u ADRC−PDTSFC2 .Considering the variance of the objective functions, the ranking of the objective functions in Table 3 is J e,u ADRC−PDTSFC2 < J e,u ADRC−PDTSFC1 < J e,u ADRC .
The third case study is characterized by setting the penalty factor this time to λ = 0.001.The optimization problem in (16) and the membership function parameters set to B êc = 0.1400, B ∆êc = 3.2774, B êθ = 0.1400, B ∆êθ = 0.0240 for the CS with ADRC-PDTSFC2 controller.As in the case studies above, the experimental results for the CSs with ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers are presented in Figure 9 in terms of the following plots: in Figure 9a) the control signal u c for ADRC is with red line, u c for ADRC-PDTSFC1 is with blue line and u c for ADRC-PDTSFC2 is with green line; in Figure 9b) the reference trajectory y * c is with thick black dashed line, the controlled output y c for ADRC is with red line, y c for ADRC-PDTSFC1 is with blue line and y c for ADRC-PDTSFC2 is with green line; in Figure 9c) the control signal u θ for ADRC is with red line, u θ for ADRC-PDTSFC1 is with blue line and u θ for ADRC-PDTSFC2 is with green line; in Figure 9d) the reference trajectory y * θ is with thick black dashed line, the controlled output y θ for ADRC is with red line, y θ for ADRC-PDTSFC1 is with blue line and y θ for ADRC-PDTSFC2 is with green line.
To have a fair comparison between the obtained results of the CSs with ADRC, ADRC-PDTSFC1 and ADRC-PDTSFC2 controller, and mitigate the effects of random disturbances that occur unexpectedly during the experiments, the graphical results presented in Figure 9 were averaged ten times.Using the same presentation style in the first two case studies, the results obtained in the third case study representing the average and the variance of the objective functions u e J , of the CSs with ADRC, ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers are presented in Table 4.As in the previous case studies, the analysis of the results given in Figure 9 and Table 4 allows for the conclusion that the CSs with ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers have similar performances, and both perform better than the CS with ADRC controller.The following relationships between the average values of the CSs with ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers are expressed: Considering the above, the variance and the ranking of the objective functions in Table 4 is

Conclusions
Using data-driven control as a starting point, this paper proposed a fusion of control methods which are meant to improve energy consumption for a TCS, resulting in three data-driven controllers, namely ADRC, ADRC-PDTSFC1 and ADRC-PDTSFC2.The new ADRC-PDTSFC structures improve the overall CS performance by compensating for shortcomings of ADRC.The data-driven controllers were validated through real-time experiments to control the cart position and the arm angular position of the MIMO nonlinear TCS laboratory equipment.The data-driven controllers were implemented as two independent SISO control system structures that are running in parallel.Using the same presentation style in the first two case studies, the results obtained in the third case study representing the average and the variance of the objective functions J e,u of the CSs with ADRC, ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers are presented in Table 4.As in the previous case studies, the analysis of the results given in Figure 9 and Table 4 allows for the conclusion that the CSs with ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers have similar performances, and both perform better than the CS with ADRC controller.The following relationships between the average values of the CSs with ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers are expressed: J e,u ADRC 1.375J e,u ADRC−PDTSFC1 and J e,u ADRC 1.336J e,u ADRC−PDTSFC2 .Considering the above, the variance and the ranking of the objective functions in Table 4 is J e,u ADRC−PDTSFC1 < J e,u ADRC−PDTSFC2 < J e,u ADRC .

Conclusions
Using data-driven control as a starting point, this paper proposed a fusion of control methods which are meant to improve energy consumption for a TCS, resulting in three data-driven controllers, namely ADRC, ADRC-PDTSFC1 and ADRC-PDTSFC2.The new ADRC-PDTSFC structures improve the overall CS performance by compensating for shortcomings of ADRC.The data-driven controllers were validated through real-time experiments to control the cart position and the arm angular position of the MIMO nonlinear TCS laboratory equipment.The data-driven controllers were implemented as two independent SISO control system structures that are running in parallel.
This paper also suggested the improvement of the energy consumption for TCS, therefore penalty factor λ was inserted in the proposed objective function in (16) in order to reduce the energy consumption.Three case studies with three values of the penalty factor λ were defined discussed in order to determine how the energy consumption influences the output tracking performance of TCS.Considering the experimental results given in Section 5 and illustrated in Figures 6-8 and Tables 2-4, the conclusion is that if the practitioner wants to minimize the energy consumption, a large value of the penalty factor should be used, but this can seriously affect the tracking performance defined in cart and arm angular positions.If the penalty factor in the objective function is set to a small value, the energy consumption will increase, but also the CS tracking performance will be improved.Therefore, the practitioner must achieve a compromise to energy consumption and output tracking performance of TCS.
The main advantage of the proposed ADRC-PDTSFC structures in comparison with the ADRC structure is that both ADRC-PDTSFC1 and ADRC-PDTSFC2 improve the overall tracking performances of the TCS for any penalty parameter λ used in the objective function.In the first case study the ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers are 1.375 times better than the ADRC controller as far as the average value of the minimum objective function is concerned.In the second case study the ADRC-PDTSFC1 is 1.385 times better than the ADRC controller and the ADRC-PDTSFC2 controller is 1.437 times better than the ADRC controller.And in the third case study the ADRC-PDTSFC1 is 1.375 times better than the ADRC controller and the ADRC-PDTSFC2 controller is 1.336 times better than the ADRC controller.In all three case studies the CSs with ADRC-PDTSFC exhibit improved performance over the CS with ADRC controller.
The main shortcoming of the proposed ADRC-PDTSFC structures is that the number of parameters is higher compared to the initial ADRC structure.The ADRC controller applied to a SISO process has 5 parameters.The ADRC-PDTSFC1 controller applied to a SISO process has 10 parameters, that is, 5 parameters of the ADRC controller plus 3 parameters of the rule base of the fuzzy component plus 2 parameters of the input membership functions.The ADRC-PDTSFC2 controller applied to a SISO process has 12 parameters, namely 5 parameters of the ADRC controller plus 5 parameters of the rule base of the fuzzy component including the PD component plus 2 parameters of the input membership functions.

Figure 1 .
Figure 1.The control system structure with second-order discrete-time Active Disturbance Rejection Control (ADRC) algorithm.

Figure 1 .
Figure 1.The control system structure with second-order discrete-time Active Disturbance Rejection Control (ADRC) algorithm.

Figure 2 .
Figure 2. The block diagram of the proportional-derivative Takagi-Sugeno Fuzzy control (PDTSFC)1 component and the input membership functions.

Figure 2 .
Figure 2. The block diagram of the proportional-derivative Takagi-Sugeno Fuzzy control (PDTSFC)1 component and the input membership functions.

Figure 3 .
Figure 3.The control system structure with second-order discrete-time ADRC algorithm.

Figure 3 .
Figure 3.The control system structure with second-order discrete-time ADRC algorithm.

Figure 4 .
Figure 4.The control system structure with second-order discrete-time ADRC algorithm.
c u and θ u are the control signals for cart position control and arm angular position control, respectively, that is, the Pulse Width Modulation (PWM) duty cycle of the cart position and arm angular position direct current (DC) motors, the control limit values of c u and θ u are limited by the PWM duty cycle between The diagram which shows the communication between the personal computer (for digital control) and the TCS laboratory equipment and the graphical definition of main variables in(15) are presented in Figure5.

Figure 4 .
Figure 4.The control system structure with second-order discrete-time ADRC algorithm.
Figure 6 illustrates the two SISO controllers, that is, cart position controller and arm angular position controller, with ADRC, ADRC-PDTSFC1 or ADRC-PDTSFC2 algorithms, the parameter vector c Λ of cart position controller, the parameter vector θ Λ of arm angular position controller, the reference input (m) * c y for cart position control (the set-point for c y ), the reference input (rad) * θ y for angular position control (the set-point for θ y ), the estimated control error control are used in the next section in the definition of the objective function as part of the optimization problem for optimal tuning of controller parameters.

Figure 6 .
Figure 6.The Multi Input-Multi Output (MIMO) control system block diagram.
Figure 6 illustrates the two SISO controllers, that is, cart position controller and arm angular position controller, with ADRC, ADRC-PDTSFC1 or ADRC-PDTSFC2 algorithms, the parameter vector c Λ of cart position controller, the parameter vector θ Λ of arm angular position controller, the reference input (m) * c y for cart position control (the set-point for c y ), the reference input control are used in the next section in the definition of the objective function as part of the optimization problem for optimal tuning of controller parameters.

Figure 6 .
Figure 6.The Multi Input-Multi Output (MIMO) control system block diagram.

Figure 6 .
Figure 6.The Multi Input-Multi Output (MIMO) control system block diagram.

Figure 7 .
Figure 7. Experimental results in the first case study for control systems (CSs) with ADRC, ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers.The obtained results representing the average and the variance of the objective functions u e J ,

Table 2 .⋅
The measured objective functions for the CS with ADRC, ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers.Analysis of the results given in Figure7and Table2leads to the conclusion that the CSs with ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers have performance in the same range and both outperform the CS with ADRC controller.The following relationships between the average values of the objective functions of the three CSs are expressed:

Figure 7 .
Figure 7. Experimental results in the first case study for control systems (CSs) with ADRC, ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers.

Energies 2019 ,
12, x FOR PEER REVIEW 12 of 19 for the CS with ADRC controller, a) the control signal c u for ADRC is with red line, c u for ADRC-PDTSFC1 is with blue line and c u for ADRC-PDTSFC2 is with green line; in Figure 8 b) the reference trajectory * c y is with thick black dashed line, the controlled output c y for ADRC is with red line, c y for ADRC-PDTSFC1 is with blue line and c y for ADRC-PDTSFC2 is with green line; in Figure 8 c) the control signal θ u for ADRC is with red line, θ u for ADRC-PDTSFC1 is with blue line and θ u for ADRC-PDTSFC2 is with green line; in Figure 8 d) the reference trajectory * θ y is with thick black dashed line, the controlled output θ y for ADRC is with red line, θ y for ADRC-PDTSFC1 is with blue line and θ y for ADRC-PDTSFC2 is with green line.

Figure 8 .
Figure 8. Experimental results in the second case study for CSs with ADRC, ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers.

Figure 8 .
Figure 8. Experimental results in the second case study for CSs with ADRC, ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers.

Figure 9 .
Figure 9. Experimental results in the third case study for CSs with ADRC, ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers.

Table 4 .
The measured objective functions for the CS with ADRC, ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers.

Figure 9 .
Figure 9. Experimental results in the third case study for CSs with ADRC, ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers.

Table 2 .
The measured objective functions for the CS with ADRC, ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers.

Table 3 .
The measured objective functions for the CS with ADRC, ADRC-PDTSFC1 and ADRC-PDTSFC2 controllers.