Virtual Reference Feedback Tuning of Model-Free Control Algorithms for Servo Systems

: This paper proposes the combination of two data-driven techniques, namely virtual reference feedback tuning (VRFT) and model-Free Control (MFC) in terms of the VRFT of MFC algorithms dedicated to servo systems. VRFT ensures the automatic optimal computation of the parameters of three MFC algorithms represented by intelligent proportional (iP), intelligent proportional-integral (iPI), and intelligent proportional-integral-derivative (iPID) controllers. The combination of MFC and VRFT leads to a novel mixed MFC-VRFT approach. The approach is validated by experimental results related to the angular speed control of modular servo system laboratory equipment. The performance of the control systems with the MFC algorithms (iP, iPI, and iPID controllers) tuned by the mixed MFC-VRFT approach is compared with that of control systems with MFC algorithms tuned by a metaheuristics gravitational search algorithm (GSA) optimizer, and of control systems with I, PI and PID controllers optimally tuned by VRFT and GSA in the same optimization problem.


Introduction
The feature of data-driven model-free control (MFC) [1] is that it uses only the input/output (I/O) data in order to control unknown processes.MFC is also known in the literature as model-free tuning, as it approximates the process model with a so called ultra-local model, and a rapidly adaptable estimator is used in the process approximation.The control system (CS) stability is guaranteed in MFC if certain constraints are fulfilled.The control algorithms specific to MFC, i.e., the MFC algorithms, are usually referred to as intelligent controllers and the most popular ones are built around proportional (P), proportional-integral (PI), or proportional-integral-derivative (PID) controllers, and are called intelligent proportional (iP), intelligent PI (iPI), and intelligent PID (iPID) controllers.The MFC algorithms have been successfully applied to various processes, such as thermal processes [2], glycemia of type-1 diabetes [3], oscillation damping [4], wheeled autonomous vehicles [5], microalgae cultivation [6], experimental greenhouses [7], surface mounted permanent magnet synchronous motor drive systems [8], and twin rotor aerodynamic systems (TRASs) [9,10].Recent applications of MFC are reported in [11][12][13].Other model-free tuning techniques similar to MFC are active disturbance rejection control [14][15][16] and model-free adaptive control (MFAC) [17][18][19].
Virtual reference feedback tuning (VRFT) [20] is another data-driven technique that uses the I/O-collected data of the initial open-loop process, and automatically computes the parameters of linear or nonlinear controllers.The main disadvantage of VRFT is that it does not directly guarantee the Machines 2017, 5, 25 2 of 15 closed-loop CS stability.The control algorithms specific to VRFT, i.e., the VRFT algorithms, have been applied recently to TRASs [9,10], two-tank systems [21], quadrotors [22], unmanned gyroplanes [23], non-minimum phase systems [24], robotics [25], and power systems [26].
VRFT has been combined with other data-driven techniques, such as MFC [11], iterative feedback tuning [27], MFAC [28], and reinforcement Q-learning [29].According to [11] both MFC and VRFT complement each other through their features and they can exhibit improved CS performance by merging them in the so-called mixed MFC-VRFT approach [11].
One major shortcoming of MFC algorithms is that in all previously reported applications the parameters have been chosen empirically and not automatically.The problem is challenging if no process model is available, and it is mitigated here by using VRFT to compute MFC algorithms' parameters.
Building upon the previous discrete-time implementations of MFC reported in [1][2][3][4][5][6][7][8][9][10][11][12][13] and the mixed MFC-VRFT approach applied to TRASs in [10], this paper proposes a thorough experimental investigation of the new approach, spanning several MFC structures and comparisons with other model-based and model-free tuned controller structures.The main advantages of the proposed approach with respect to the state-of-the-art are: -It is a widely applicable one-shot data-driven tuning technique that ensures the experiment-based optimal tuning of MFC algorithms.-As shown in the applications of MFC algorithms [1][2][3][4][5][6][7][8][9][10][11][12][13], which have mainly achieved the discrete-time implementation of iP, iPI, and iPID controllers, the stability around a reference trajectory is obtained by appropriate gain tuning.-It can be easily generalized to other model-free control techniques.
This paper applies the MFC-VRFT approach to a nonlinear modular servo system (MSS) laboratory equipment and compares the performance of CSs with MFC algorithms (iP, iPI, and iPID controllers) tuned by the mixed MFC-VRFT approach with that of control CSs with MFC algorithms tuned by a metaheuristics gravitational search algorithm (GSA) optimizer, of CSs with I, PI, and PID controllers optimally tuned by VRFT, and CSs with I, PI, and PID controllers tuned by GSA.
The main contributions of this paper with the current literature in the field [1][2][3][4][5][6][7][8][9][10][11][12][13] are: -The mixed MFC-VRFT approach is applied here to an MSS and is experimentally validated by real-time results.-Three different MFC structures are tuned in this paper using the new approach in order to investigate how the controller structure complexity influences the CS performance.
The VRFT-based iP, iPI, and iPID controllers are compared against three other model-based and model-free tuned controllers to support the further use of the new model-free tuning approach.
The fair controller comparison is ensured by using the same optimization problem to optimally tune the parameters of both MFC algorithms and linear controllers.Three experimental case studies have been considered, (1), (2), and (3): (1) The performance of a CS with an iP controller optimally tuned by VRFT is compared with that of a CS with a model-based optimally tuned iP controller.Then the performance of the CSs with iP controllers is compared with that of a CS with an I controller, which is optimally tuned both in a model-based setting and using VRFT (i.e., in a data-driven model-free setting).(2) The same comparisons and controller tunings are carried out as in the experimental case study (1) but two iPI controllers and two PI controllers are involved.(3) The experimental case study (1) is applied, but two iPID controllers and two PID controllers are involved.
The paper is organized as follows: the MFC, VRFT, and MFC-VRFT approaches are briefly presented in Section 2, the application of the MFC-VRFT approach and the experimental results on the MSS laboratory equipment are given in Section 3, and the conclusions are highlighted in Section 4.

Overview of MFC
The discrete-time first-order ultra-local model of the unknown single input-single output process is [1]: where The control laws specific to iP, iPI, and iPID controllers are [1]: where , and y * (t d ) ∈ is the reference input to the feedback CS regarded as a desired trajectory to be tracked.By abuse of notation, K 1 , K 2 , K 3 ∈ are called the P, integral (I), and derivative (D) gains of the iP, iPI, and iPID controllers as they are easily computable in terms of their classical counterparts that appear in the recurrent equations in the error term ε(t d ), such as those in the right-hand side of Equation (2).ε(t d ) indicates the tracking error, ε(t d ) ∈ {ε iP (t d ), ε iPI (t d ), ε iPID (t d )} and, by case, the subscripts iP, iPI, and iPID will be added to ε(t d ) and y(t d ): F(t d ) is computed from the I/O data and Equation (1) using: where the subscripts can be iP, iPI, and iPID can be added to F(t d ), y(t d ), and u(t d ) depending on the MFC algorithm.A noise-corrupted y(t d ) can compromise the precision of the estimate F(t d ).Adding a lowpass filter on the measured output will remedy this problem and will also extend the process I/O dynamics, whose bandwidth will decrease.The extended process dynamics are not an issue since unmodeled dynamics are accounted for within the estimated uncertainty F(t d ).Subsequently, the closed-loop bandwidth will decrease and the CS performance requirements from y * (t d ) to y(t d ) have to be re-matched with the dynamics of the extended process.The difference between F(t d ) and F(t d ) is expressed as: where δ(t d ) ∈ is considered a negligible random disturbance in the controller design and tuning [9,10], Substituting the iP, iPI, and iPID control laws given in Equation (2) in the model of Equation ( 1) leads to the dynamics equations of the closed-loop CSs with iP, iPI, and iPID controllers: The CS structure with MFC algorithms, namely iP, iPI, and iPID controllers, is illustrated in Figure 1.Substituting the iP, iPI, and iPID control laws given in Equation (2) in the model of Equation (1) leads to the dynamics equations of the closed-loop CSs with iP, iPI, and iPID controllers: The CS structure with MFC algorithms, namely iP, iPI, and iPID controllers, is illustrated in Figure 1.The substitution of ) ( d t δ defined in Equation ( 5) into Equation ( 6) and in Equation ( 3) results in the dynamics equation of the tracking error for the CSs with iP, iPI, and iPID controllers: As shown in [9] and [10], MFC guarantees the CS stability if the roots of the characteristic polynomials resulted from Equation (7) are inside the unit circle and the controlled output will track the reference input if the tracking error is converging to zero.

Overview on Nonlinear VRFT
The goal of VRFT is to automatically tune the parameters of a linear or nonlinear controller such that the closed-loop CS tracks a reference model [30,31].With this respect, only the I/O data )} ( ), ( { collected after an open-loop experiment is used as follows: first, a rich frequency spectrum signal must be applied as input on the open-loop process; second, the designer must determine a linear or nonlinear reference model that should specify the behavior of the closed-loop CS.
The performance index used in VRFT is the model reference objective function [31]: where N is the number of samples, is the output of the nonlinear process modeled by the recurrent equation: )  The substitution of δ(t d ) defined in Equation ( 5) into Equation ( 6) and in Equation (3) results in the dynamics equation of the tracking error for the CSs with iP, iPI, and iPID controllers: As shown in [9,10], MFC guarantees the CS stability if the roots of the characteristic polynomials resulted from Equation (7) are inside the unit circle and the controlled output will track the reference input if the tracking error is converging to zero.

Overview on Nonlinear VRFT
The goal of VRFT is to automatically tune the parameters of a linear or nonlinear controller such that the closed-loop CS tracks a reference model [30,31].With this respect, only the I/O data {u(t d ), y(t d )} collected after an open-loop experiment is used as follows: first, a rich frequency spectrum signal must be applied as input on the open-loop process; second, the designer must determine a linear or nonlinear reference model that should specify the behavior of the closed-loop CS.
The performance index used in VRFT is the model reference objective function [31]: )), C θ is the nonlinear controller parameterized by vector θ, and the detailed form of the nonlinear controller used in VRFT is: with n uc and n ec are the known orders of the nonlinear controller parameterized by the vector θ.
The control (or feedback) error of the closed-loop CS is: where r(t d ) is the reference input applied to the closed-loop CS.The output of the nonlinear reference model is: where the reference model m is assumed to be invertible and the orders of output and reference input are n ym and n rm .The unknown linear/nonlinear process is assumed to be I/O stable.The virtual reference input is r(t d ) = m −1 (y(t d )), and it is computed assuming that the output trajectories of the closed-loop CS and the reference model are similar.The virtual control error e(t d ) is next calculated as e(t d ) = r(t d ) − y(t d ).
According to [31], if e(t d ) is filtered through a richly-parameterized nonlinear controller with parameter vector θ, the resulting signal should be similar to the open-loop control signal u(t d ).
The controller parameter θ is determined in terms of minimizing the objective function: Unlike in the classical linear [20] and nonlinear [30] VRFT where J MR (θ) and J VRFT (θ) can be made approximately equal in the vicinity of their minimizing arguments with the use of a prefilter, [31] suggests that, without using the prefilter, J MR (θ) and J VRFT (θ) can be minimized by the same argument in the case of a richly parameterized controller.
The design approach of the VRFT structure consists of the following steps: Step 1.1.Set the initial rich spectrum frequency signal u(t d ) that must be applied as input on the open-loop process and then collect the I/O data {u(t d ), y(t d )} after an open-loop experiment.
Step 1.2.Set a linear or nonlinear reference model m that should specify the behavior of the closed-loop CS and leads to the virtual reference vector r(t d ) = m −1 (y(t d )), such that the output trajectories of the closed-loop CS and the reference model are similar.
Step 1.3.Determine the parameter vector θ by minimizing the objective function in (13).
The CS structure controlled with VRFT is presented in Figure 2.
Machines 2017, 5, 25 5 of 15 with uc n and ec n are the known orders of the nonlinear controller parameterized by the vector θ .
The control (or feedback) error of the closed-loop CS is: where is the reference input applied to the closed-loop CS.The output of the nonlinear reference model is: where the reference model m is assumed to be invertible and the orders of output and reference input are ym n and rm n .
The unknown linear/nonlinear process is assumed to be I/O stable.The virtual reference input is )), ( ( ) ( and it is computed assuming that the output trajectories of the closed-loop CS and the reference model are similar.The virtual control error is filtered through a richly-parameterized nonlinear controller with parameter vector θ, the resulting signal should be similar to the open-loop control signal . The controller parameter θ is determined in terms of minimizing the objective function: Unlike in the classical linear [20] and nonlinear [30] VRFT where can be made approximately equal in the vicinity of their minimizing arguments with the use of a prefilter, [31] suggests that, without using the prefilter, can be minimized by the same argument in the case of a richly parameterized controller.The design approach of the VRFT structure consists of the following steps: Step 1.1.Set the initial rich spectrum frequency signal that must be applied as input on the open-loop process and then collect the I/O data )} ( ), ( { after an open-loop experiment.
Step 1.2.Set a linear or nonlinear reference model m that should specify the behavior of the closed-loop CS and leads to the virtual reference vector , such that the output trajectories of the closed-loop CS and the reference model are similar.
Step 1.3.Determine the parameter vector θ by minimizing the objective function in (13).
The CS structure controlled with VRFT is presented in Figure 2.

Mixed MFC-VRFT Approach
This sub-section presents how the VRFT approach is used to automatically compute the parameters of the MFC iP, iPI, and iPID controllers.Unlike [10], where VRFT solely tunes an iP controller, here VRFT also tunes the iPI and iPID controllers, to investigate the increased controller structure complexity over the CS performance.
The equation specific to the MFC law using Equation ( 4) is:

Mixed MFC-VRFT Approach
This sub-section presents how the VRFT approach is used to automatically compute the parameters of the MFC iP, iPI, and iPID controllers.Unlike [10], where VRFT solely tunes an iP controller, here VRFT also tunes the iPI and iPID controllers, to investigate the increased controller structure complexity over the CS performance.The equation specific to the MFC law using Equation ( 4) is: where } is the notation of the block of the MFC law that processes past samples of the feedback error (in similar fashion to P, PI, or PID controllers), and its detailed form is: where . Substituting the estimate of F(t d ) in u(t d ), the detailed control law becomes: Adding and subtracting the output value y(t d ) and the one step behind sample of the desired trajectory y * (t d − 1) in Equation ( 16) and then using Equations ( 3) and ( 15), the expressions of the control laws specific to the iP, iPI, and iPID controllers are: The control laws given in Equation ( 17) are also expressed in the nonlinear forms: where g is linear affine function of its arguments, being parameterized by T regarded as tunable vector parameters (with T indicating matrix transposition) and T are regarded as the non-tunable but known vector parameters and can be treated as either time-varying parameters or exogenous disturbances.The linear affine form of g is a particular case of a nonlinear controller structure with integral character.It fits the form of Equation (10), motivating the VRFT.
The VRFT of MFC algorithms is carried out by the mixed MFC-VRFT approach, and the merge of MFC and VRFT is briefly described as follows.The reference input y * (t d ) specific to MFC and the virtual reference input r(t d ) specific to VRFT are assumed to be equivalent.The virtual error e(t d ) specific to VRFT and the tracking error ε(t d ) specific to MFC are assumed to be equivalent using the fact that e(t d ) is a part of the nonlinear controller equation given in Equation (10), where, depending by the scenario, the parameter vector with the general notation θ ∈ {θ iP , θ iPI , T .The detailed forms of the objective function (13) using these expressions of θ are: Machines 2017, 5, 25

of 15
The CS structure with mixed MFC-VRFT algorithms is presented in Figure 3.The MFC algorithms tuned by VRFT in terms of minimizing the first, second and third objective functions in Equation (19) are referred to as VRFT-iP, VRFT-iPI, and VRFT-iPID controllers, respectively.
In summary, the designer must select a linear or nonlinear reference model that specifies the behavior of the closed-loop CS.The I/O data are next collected from the process in a dedicated open-loop experiment.An iP, iPI, or iPID controller structure is selected and the parameter α is chosen by the user, then the other tunable parameters are automatically determined by the GSA-based minimization of the objective function in Equation (13).The obtained parameters are used in the closed-loop CSs to track the reference model.
The design approach of the MFC-VRFT structure consists of the following steps: Step 2. Step 2.3.Set a linear or nonlinear reference model m that should specify the behavior of the closed-loop CS and leads to the virtual reference vector r(t d ) = m −1 (y(t d )), such that the output trajectories of the closed-loop CS and the reference model are similar.
Step 2.4.Determine the parameter vector θ by minimizing the objective function in Equation (13).
Machines 2017, 5, 25 7 of 15 The CS structure with mixed MFC-VRFT algorithms is presented in Figure 3.The MFC algorithms tuned by VRFT in terms of minimizing the first, second and third objective functions in Equation ( 19) are referred to as VRFT-iP, VRFT-iPI, and VRFT-iPID controllers, respectively.
In summary, the designer must select a linear or nonlinear reference model that specifies the behavior of the closed-loop CS.The I/O data are next collected from the process in a dedicated open-loop experiment.An iP, iPI, or iPID controller structure is selected and the parameter α is chosen by the user, then the other tunable parameters are automatically determined by the GSA-based minimization of the objective function in Equation (13).The obtained parameters are used in the closed-loop CSs to track the reference model.
The design approach of the MFC-VRFT structure consists of the following steps: Step after an open-loop experiment.
Step 2.3.Set a linear or nonlinear reference model m that should specify the behavior of the closed-loop CS and leads to the virtual reference vector , such that the output trajectories of the closed-loop CS and the reference model are similar.
Step 2.4.Determine the parameter vector θ by minimizing the objective function in Equation ( 13).

The Modular Servo System Equipment
The VRFT-iP, VRFT-iPI, and VRFT-iPID controllers tuned by the mixed MFC-VRFT approach presented in Section 2 are validated on single input-single output MSS laboratory equipment in order to control its velocity.The MSS is illustrated in Figure 4.The state-space model that describes the nonlinear MSS laboratory equipment and the MSS process parameters are given in [32] and [33].The equipment is able to collect discrete-time I/O data from the process in real-time.Therefore, the mixed MFC-VRFT approach will be applied to control this nonlinear process.

The Modular Servo System Equipment
The VRFT-iP, VRFT-iPI, and VRFT-iPID controllers tuned by the mixed MFC-VRFT approach presented in Section 2 are validated on single input-single output MSS laboratory equipment in order to control its velocity.The MSS is illustrated in Figure 4.
The CS structure with mixed MFC-VRFT algorithms is presented in Figure 3.The MFC algorithms tuned by VRFT in terms of minimizing the first, second and third objective functions in Equation ( 19) are referred to as VRFT-iP, VRFT-iPI, and VRFT-iPID controllers, respectively.
In summary, the designer must select a linear or nonlinear reference model that specifies the behavior of the closed-loop CS.The I/O data are next collected from the process in a dedicated open-loop experiment.An iP, iPI, or iPID controller structure is selected and the parameter α is chosen by the user, then the other tunable parameters are automatically determined by the GSA-based minimization of the objective function in Equation (13).The obtained parameters are used in the closed-loop CSs to track the reference model.
The design approach of the MFC-VRFT structure consists of the following steps: Step , such that the output trajectories of the closed-loop CS and the reference model are similar.
Step 2.4.Determine the parameter vector θ by minimizing the objective function in Equation ( 13).

The Modular Servo System Equipment
The VRFT-iP, VRFT-iPI, and VRFT-iPID controllers tuned by the mixed MFC-VRFT approach presented in Section 2 are validated on single input-single output MSS laboratory equipment in order to control its velocity.The MSS is illustrated in Figure 4.The state-space model that describes the nonlinear MSS laboratory equipment and the MSS process parameters are given in [32] and [33].The equipment is able to collect discrete-time I/O data from the process in real-time.Therefore, the mixed MFC-VRFT approach will be applied to control this nonlinear process.The state-space model that describes the nonlinear MSS laboratory equipment and the MSS process parameters are given in [32,33].The equipment is able to collect discrete-time I/O data from the process in real-time.Therefore, the mixed MFC-VRFT approach will be applied to control this nonlinear process.

Experimental Validation
Three experimental case studies, (1), ( 2) and (3), have been considered.Each case study consists of four experimental control scenarios where: (1) the performance of the model-free VRFT-iP and VRFT-I controllers is compared with that of model-based tuned iP and I controllers; (2) the performance of the model-free VRFT-iPI and VRFT-PI controllers is compared with that of model-based tuned iPI and PI controllers; and (3) the performance of the model-free VRFT-iPID and VRFT-PID controllers is compared with that of model-based tuned iPID and PID controllers.
In all experimental case studies, the sampling time has been set to T s = 0.1s, the number of samples to N = 900, and the reference input y * (t d ) to a combination of sine, step and ramp signals: (20) the parameter α to α = 80 has been set in step 2.1, and the discrete transfer function of the linear reference model m(z −1 ) in all experimental case studies has been set to according to steps 1.2 and 2.3 after applying steps 1.1 and 2.2, respectively: The performance of CSs with controllers tuned by the mixed MFC-VRFT approach is conditioned by both the choice of the reference model [10,28], and the choice of the initial rich spectrum frequency signal that is applied as input to the open-loop process [31].
The optimization problem defined for the sake of fair comparisons of the performance of VRFT-iP, VRFT-iPI, and VRFT-iPID controllers, VRFT-I, VRFT-PI, and VRFT-PID controllers, and I, PI, and PID controllers, is: where the vector variable (the parameter vector) θ 1 depends on the selected experimental case study as follows: T in experimental case study (2), and The tunable parameters of VRFT-iP, VRFT-iPI, and VRFT-iPID controllers are automatically determined by minimizing the objective function in Equation ( 13) using GSA, with the GSA parameters and implementations given in [34,35], and the tunable parameters computed using only the I/O data collected from the nonlinear MSS of the initial open-loop experiment.On the other hand, the parameters of the iP, iPI and iPID controllers and those of the I, PI, and PID controllers are determined by solving the optimization problem in (22) in a model-based setting using the nonlinear MSS process.The problem in Equation ( 22) is solved via GSA, with the objective function evaluated by the simulation of the CS behavior using the reference input illustrated as follows.The parameters of the VRFT-I, VRFT-PI, and VRFT-PID controllers are computed by minimizing the objective function in Equation ( 13) through a linear least-squares regression.Therefore, the parameters of VRFT-iP, VRFT-iPI, VRFT-iPID, VRFT-I, VRFT-PI, and VRFT-PID controllers are tuned in a model-free VRFT-based setting applying step 1.3.
Experimental case study (1) consists of four experimental scenarios, where the MSS laboratory equipment is controlled using the optimally-tuned parameters of two model-free VRFT-iP and VRFT-I = K 1 = 0.0119, and the model-based I controller parameter is θ opt 1 = K 1 = 0.0066.The obtained results have been averaged ten times in order to eliminate the effects of random disturbances.The average and the variance of J E in the experimental case study (1) are presented in Table 1. Figure 5 illustrates the control signal and the controlled output (the angular velocity) of the CSs with VRFT-iP, iP, VRFT-I and I controllers.Experimental case study (2) also consists of four experimental scenarios, where the MSS laboratory equipment is controlled using the optimally tuned parameters of VRFT-iPI, iPI, VRFT-PI and PI controllers.The VRFT-iPI controller parameters are θ opt The obtained results have been averaged ten times in order to eliminate the effects of random disturbances.The average and the variance of E J in the experimental case study (1) are presented in Table 1. Figure 5 illustrates the control signal and the controlled output (the angular velocity) of the CSs with VRFT-iP, iP, VRFT-I and I controllers.
, the parameters of the model-free VRFT-PI controller are As in the previous case study, the results in the experimental case study (2) have been averaged ten times.Table 2 gives the average and the variance of E J in experimental case study (2).The control signal and the controlled output of the CSs with VRFT-iPI, iPI, VRFT-PI and PI controllers are presented in Figure 6.As in the previous case study, the results in the experimental case study (2) have been averaged ten times.Table 2 gives the average and the variance of J E in experimental case study (2).The control signal and the controlled output of the CSs with VRFT-iPI, iPI, VRFT-PI and PI controllers are presented in Figure 6.Experimental case study (3) also sums four experimental scenarios, where the MSS laboratory equipment is controlled using the optimally tuned parameters of VRFT-iPID, iPID, and VRFT-PID controllers, and a model-based PID controller.The model-free VRFT-iPID controller parameters are The obtained results given in Table 3 and Figure 7 for experimental case study (3) have also been averaged ten times.The average and the variance of E J are presented in Table 3.The control signal and the controlled output of the CSs with VRFT-iPID, iPID, VRFT-PID and PID controllers are shown in Figure 7.The resulted characteristic polynomials of the VRFT-iP, iP, VRFT-iPI, iPI, VRFT-iPID, and iPID controllers together with the roots of the characteristic polynomials resulted from (7) that guarantee the MFC CS stability are presented in Table 4.  3 and Figure 7 for experimental case study (3) have also been averaged ten times.The average and the variance of J E are presented in Table 3.The control signal and the controlled output of the CSs with VRFT-iPID, iPID, VRFT-PID and PID controllers are shown in Figure 7.The resulted characteristic polynomials of the VRFT-iP, iP, VRFT-iPI, iPI, VRFT-iPID, and iPID controllers together with the roots of the characteristic polynomials resulted from (7) that guarantee the MFC CS stability are presented in Table 4.

Discussion of Experimental Results
The results obtained in experimental case study (1) show that that the CSs with VRFT-iP and iP controllers exhibit similar performance with the iP controller producing an insignificant advantage.Inserting the controller abbreviations as subscripts in the objective function, Table 1 . The responses given in Figure 5 illustrate the same, namely, the controlled outputs the CSs with VRFT-iP and iP controllers almost overlap with the reference trajectory signal.Unlike the model-free controllers, the CSs with I controllers tuned via model-free VRFT or GSA exhibit low performance, resulting in some small oscillations at the output, whereas this controller has a single parameter compared with the iP controller.
The results obtained in experimental case study (2) reveal that VRFT-iPI controller show lower performing than the iPI one, with small differences, namely . This is a somehow normal result, whereas the parameters of the iPI controller have been obtained using the reference trajectory and the mathematical model of the MSS, while the parameters of the VRFT-iPI controller were computed using only the I/O data from the initial open-loop experiment.In this experimental case study, the CS with model-based optimized linear PI controller has better performance than the CS with nonlinear model-free controller, reflected in the objective function ratio, Unlike the PI controller, the performance of the CS with VRFT-PI controller is somewhere between those of the CSs with VRFT-iPI and iPI controllers as shown in Table 2 and Figure 6. Figure 6 illustrates that all four controlled outputs overlap the reference trajectory.
The results given in Table 3 and Figure 7 for experimental case study (3) outline that the CS with iPID controller tracks better the reference trajectory than the CS with VRFT-iPID controller, i.e.,

Discussion of Experimental Results
The results obtained in experimental case study (1) show that that the CSs with VRFT-iP and iP controllers exhibit similar performance with the iP controller producing an insignificant advantage.Inserting the controller abbreviations as subscripts in the objective function, Table 1 indicates that J E VRFT−iP ≈ 1.04J E iP .The responses given in Figure 5 illustrate the same, namely, the controlled outputs the CSs with VRFT-iP and iP controllers almost overlap with the reference trajectory signal.Unlike the model-free controllers, the CSs with I controllers tuned via model-free VRFT or GSA exhibit low performance, resulting in some small oscillations at the output, J E I ≈ 14.55J E VRFT−I , whereas this controller has a single parameter compared with the iP controller.
The results obtained in experimental case study (2) reveal that VRFT-iPI controller show lower performing than the iPI one, with small differences, namely J E VRFT−iPi ≈ 1.58J E iPI .This is a somehow normal result, whereas the parameters of the iPI controller have been obtained using the reference trajectory and the mathematical model of the MSS, while the parameters of the VRFT-iPI controller were computed using only the I/O data from the initial open-loop experiment.In this experimental case study, the CS with model-based optimized linear PI controller has better performance than the CS with nonlinear model-free controller, reflected in the objective function ratio, J E iPI ≈ 2.02J E PI .Unlike the PI controller, the performance of the CS with VRFT-PI controller is somewhere between those of the CSs with VRFT-iPI and iPI controllers as shown in Table 2 and Figure 6. Figure 6 illustrates that all four controlled outputs overlap the reference trajectory.
The results given in Table 3 and Figure 7 for experimental case study (3) outline that the CS with iPID controller tracks better the reference trajectory than the CS with VRFT-iPID controller, i.e., J E VRFT−iPID ≈ 1.54J E iPID .This is an expected result as the controlled output of the CS with iPID controller overlaps the reference trajectory.The comparison of the best performance of the model-free CSs, namely, the CS with iPID controller, with that of the CS with model-based PID controller, points out that the controlled output of the CS with model-based PID controller exhibits better compared to that of the CS with iPID controller, reflected in J E iPID ≈ 1.93J E PID .The performance of the CS with VRFT-PID controller is slightly better than that of the CS with VRFT-iPID controller, J E VRFT−iPID ≈ 1.13J E VRFT−PID , but not than that of the CSs with iPID and PID controllers obtained in a model-based setting.Figure 7 illustrates that the differences are almost imperceptible, whereas the controlled outputs almost overlap.
The analysis of the performance of CSs with MFC algorithms tuned by VRFT shows that the best results are obtained for VRFT-iP, followed by VRFT-iPID and VRFT-iPI, and the comparison of the objective function values is J E VRFT−iPI ≈ 1.30J E VRFT−iP ≈ 1.07J E VRFT−iPID .
The comparison of the performance of CSs with MFC algorithms tuned by GSA highlights that the differences between the objective function values for iP, iPI, and iPID controllers are insignificant, as illustrated by J E iP ≈ 1.22J E iPID ≈ 1.16J E iPI .Concluding, the comparison of MFC algorithms tuned by model-free VRFT and model-based GSA proves that no serious improvement is achieved if the intelligent MFC adaptive structure complexity expressed in number of parameters is increased.
The comparison of the performances of CSs with linear I, PI, and PID controllers tuned by VRFT shows that the VRFT-PI and the VRFT-PID controllers perform eight times better than the VRFT-I controller and the objective function comparison is J E VRFT−I ≈ 8.07J E VRFT−PI ≈ 8.23J E VRFT−PID .
The comparison of CSs with linear I, PI and PID controllers tuned in a model-based setting reveals that the I controller has lower performance than the PI and PID ones.The CSs with PI and PID controllers perform ten times better than the CS with I controller, namely, J E I ≈ 10.47J E PI ≈ 10.52J E PID .
The experimental results presented in case studies (1), (2), and (3) indicate that, on one hand, better reference model tracking performance is achieved if the model-based tuning of iP, iPI, and iPID controllers is carried out versus the model-free MFC-VRFT approach, although the differences between model-free VRFT tuning and model-based optimal tuning (using GSA) of MFC algorithms are not high.This is considered normal because using the process model knowledge is expected to give superior performance.
On the other hand, the model-based optimal tuning results show that the classical linear PI and PID controllers outperform the iPI and iPID controllers, respectively, except the I versus the iP controller.In the case of I and iP controllers it is in fact normal that the latter is better as it actually has an integral component in the first expression in Equation ( 17), although it is only called iP.As such, even if it has only one tunable parameter, the iP controller more accurately resembles a pure integral controller.Thus, the complexity of MFC algorithms is not really justified in the iPI and iPID cases.According to the obtained results, in the model-free case, better performance is achieved using the VRFT-iP controller versus the VRFT-iPI, VRFT-iPID, iPID, VRFT-I, VRFT-PI, and VRFT-PID controllers, although the VRFT-iP controller has only one tunable parameter.As specified in the previous section, the MFC algorithms have an additional parameter, α, which, in this case, was not subjected to tuning, but it has been set by the designer, indicating a further shortcoming of MFC algorithms.
Since no major differences exist between model-based and model-free tuning in this application, it is more advantageous to tune the controller parameters using model-free VRFT.Our new mixed MFC-VRFT approach saves time as it automatically tunes the parameters of the MFC algorithms.
Arguably, the VRFT approach translates the selection of MFC parameters into the selection of reference model parameters, which may be considered just a problem reformulation.However, it is always preferable to prescribe performance in terms of the easily-interpretable indices of reference model responses, such as response time and overshoot.Nevertheless, this required performance has to be matched with the process dynamics.

Figure 2 .
Figure 2. The CS structure controlled with VRFT.

Figure 2 .
Figure 2. The CS structure controlled with VRFT.

1 = 1 =
two model-based iP and I controllers.The parameter of the VRFT-iP controller is θ opt θ iP1 = K 1 = −0.9018, the parameter of the iP controller is θ opt θ iP1 = K 1 = −0.6501, the parameter of the model-free VRFT-I controller is θ opt 1 and the parameters of the model-based PI controller are θopt 1 = [K 1 K 2 ] T = [0.0653− 0.0599] T .

Figure 5 .
Figure 5. Experimental results of case study (1): control signal u of CS with VRFT-iP (black), CS with iP (red), CS with VRFT-I (green), and CS with I (blue) controllers in (a); reference trajectory d y (solid black dashed), controlled output y of CS with VRFT-iP (black), CS with iP (red), CS with VRFT-I (green), and CS with I (blue) controllers in (b).

Figure 5 .
Figure 5. Experimental results of case study (1): control signal u of CS with VRFT-iP (black), CS with iP (red), CS with VRFT-I (green), and CS with I (blue) controllers in (a); reference trajectory y d (solid black dashed), controlled output y of CS with VRFT-iP (black), CS with iP (red), CS with VRFT-I (green), and CS with I (blue) controllers in (b).

Figure 6 .
Figure 6.Experimental results of case study (2): control signal u of CS with VRFT-iPI (black), CS with iPI (red), CS with VRFT-PI (green), and CS with PI (blue) controllers in (a); reference trajectory d y (solid black dashed), controlled output y of CS with VRFT-iPI (black), CS with iPI (red), CS with VRFT-PI (green), and CS with PI (blue) controllers in (b).

.
This is an expected result as the controlled output of the CS with iPID controller overlaps the reference trajectory.The comparison of the best performance of the
y(t d ) ∈ is the controlled output in the case of iP, iPI, and iPID controllers, respectively, with the notation y(t d ) ∈ {y iP (t d ), y iPI (t d ), y iPID (t d )}; u(t d ) ∈ is the control signal, u(t d ) ∈ {u iP (t d ), u iPI (t d ), u iPID (t d )}; F(t d ) ∈ are the results through the I/O-estimated values and account for the unmodeled dynamics and unknown disturbances, with the notation F(t d ) ∈ {F iP (t d ), F iPI (t d ), F iPID (t d )}; t d ∈ Z, t d ≥ 0 is the discrete time argument; and α ∈ is a constant chosen by the CS designer such that ∆y(t d + 1) = y(t d + 1) − y(t d ) and α u(t d ) have the same order of magnitude.
1. Set the design parameter α > 0 such that ∆y(t d + 1) = y(t d + 1) − y(t d ) and αu(t d ) have the same order of magnitude.Step 2.2.Set the initial rich spectrum frequency signal u(t d ) that must be applied as input on the open-loop process and then collect the I/O data {u(t d ), y(t d )} after an open-loop experiment.