Optimal Design of Fuzzy Systems Using Di ﬀ erential Evolution and Harmony Search Algorithms with Dynamic Parameter Adaptation

: This paper presents a study of two popular metaheuristics, namely di ﬀ erential evolution (DE) and harmony search (HS), including a proposal for the dynamic modiﬁcation of parameters of each algorithm. The methods are applied to two cases, ﬁnding the optimal design of a fuzzy logic system (FLS) applied to the optimal design of a fuzzy controller and to the optimization of mathematical functions. A fuzzy logic controller (FLC) of the Takagi–Sugeno type is used to ﬁnd the optimal design in the membership functions (MFs) for the stabilization problem of an autonomous mobile robot following a trajectory. A comparative study of the results for two modiﬁed metaheuristic algorithms is presented through analysis of results and statistical tests. Results show that, statistically speaking, optimal fuzzy harmony search (OFHS) is better in comparison to optimal fuzzy di ﬀ erential evaluation (OFDE) for the two presented study cases.


Introduction
Recently, the harmony search algorithm (HS) has proven to be an interesting method in solving complex problems in intelligent computing, and several authors are focused on implementing this algorithm in different fields on computer science. For example, advances for this algorithm are presented by Kim et al. in [1], HS-based management of distributed energy and storage systems in microgrids is described by Ceylan et al. in [2], an HS is used for neural networks to improve fraud detection in banking system by Daliri in [3], a discrete HS for flexible job shop scheduling problem with multiple objectives is explained by Gao et al. in [4], a HS for optimal design of PID controller is described by Kayabekir et al. in [5], an HS method to solve the vehicle routing problem by Liu et al. in [6], an improved HS is described by Ouyang et al. in [7], an HS with dynamic adaptation of Parameters for problem control presented by Peraza et al. in [8], an HS for feature selection for facial emotion recognition using cosine similarity by Saha et al. in [9], a novel HS and its application to data clustering presented by Talaei et al. in [10], a case study to test a fuzzy HS is described by Valdez et al. in [11], and an improved differential-based HS with linear dynamic domain is presented by Zhu et al. in [12].
Another important meta-heuristic that has proven to be an excellent algorithm based on its implementation results in different areas of intelligent computing is the differential evolution (DE) algorithm, and some important works using DE include an Adaptive DE with novel mutation strategies in multiple sub-populations presented by Cui et al. in [13], a novel DE for solving constrained engineering optimization problems presented by Mohamed in [14], a hybrid real-code population-based incremental learning and DE for many-objective optimization of an automotive floor-frame described by Pholdee et al. in [15], time series forecasting for building energy consumption using weighted support vector regression with DE optimization technique presented by Zhang et al. in [16], an improved adaptive DE for continuous optimization outlined by Yi et al. in [17], and an adaptive DE with sorting crossover rate for continuous optimization problems presented by Zhou et al. in [18].
To validate the efficiency in the performance of the metaheuristics algorithms, two cases are used. The first one is in the area of stabilization of non-linear plants, and the second one is the area of mathematical functions. For the first case, some interesting related works can be mentioned-a Particle Swarm Optimization-Modified Frequency Bat (PSO-MFB) is implemented for stable path planning of an autonomous mobile robot (AMR) as presented by Ajeil et al. [19], the use of multiple sensors for AMR navigation is presented by Hoang et al. [20], a mobile robot path-planning using oppositional-based improved firefly algorithm under cluttered environment is presented by Panda et al. [21], fuzzy sets in dynamic adaptation of parameters of a bee colony optimization for controlling the trajectory of an AMR are presented by Amador-Angulo et al. [22], a robot path planning using modified artificial bee colony algorithm is presented by Nayyar et al. [23], and a design and implementation of a fuzzy path optimization system for omnidirectional AMR in real-time is presented by Cuevas et al. [24]. An interesting work where DE is used in the field of non-linear plants is for robot path planning as presented by Jain et al. [25]. In this work, the optimal design of the MFs is search for a fuzzy logic system (FLS) [26,27], specifically a Takagi-Sugeno fuzzy controller is used in the experiments, this inference mechanisms in the FLS is an excellent tool a comparative of fuzzy sets, for example; a Takagi-Sugeno fuzzy logic controller (FLC) for a Liu-Chen four-scroll chaotic system is presented by Vaidyanathan et al. in [28], a Whale optimization algorithm-based Sugeno FLC for fault ride-through improvement of grid-connected variable speed wind generators is presented by Qais et al. in [29], a Fault tolerant trajectory tracking control design for interval-type-2 Takagi-Sugeno fuzzy logic system is outlined by Maalej et al. in [30], a Sugeno-Mamdani fuzzy system based soft computing approach toward sensor node localization with optimization is presented by Kumar et al. in [31], a hybrid technique of Mamdani and Sugeno-based fuzzy interference system approach is explained by Devi et al. in [32], an optimization of fuzzy logic (Takagi-Sugeno) blade pitch angle controller in wind turbines by genetic algorithm is presented by Civelek in [33], a control and balancing of two-wheeled mobile robots using Sugeno fuzzy logic in the domain of AI techniques is presented by Chouhan et al. in [34], and a stable Takagi-Sugeno fuzzy control designed by optimization is presented by Vrkalovic et al. in [35].
The second important field in intelligent computing to evaluate metaheuristic algorithms is in mathematical functions, and some interesting related works include Hussien et al. presenting a New binary whale optimization algorithm for discrete optimization problems [36]; Ochoa et al. presenting a DE with a fuzzy logic approach for dynamic parameter adjustment using benchmark functions [37]; Peraza et al. presenting an HS with dynamic adaptation of parameters for the optimization of a benchmark set of functions [38]; Rao presents the Rao algorithms-three metaphor-less simple algorithms for solving optimization problems [39]; Sulaiman et al. presenting a barnacles mating optimizer-a new bio-inspired algorithm for solving engineering optimization problems [40]; Yue et al. presenting a hybrid grasshopper optimization algorithm with invasive weed for global optimization [41]; and Perez et al. presenting a bat algorithm comparison with genetic algorithm using benchmark functions [42]. Finally, an interesting work is presented by Mulo et al. as a combination of modified HS and DE optimization techniques in economic load dispatch [43].
Therefore, it is important to analyze the study cases that have been implemented in the works previously presented [37,38]. In this paper, two metaheuristic algorithms are implemented, which are harmony search (HS) and differential evolutional (DE); the main purpose is to find the optimal parameters for each algorithm that allows the user to minimize the fitness functions in each of the study cases. Thus, an improvement on the original algorithm is proposed with the dynamic adaptation of parameters through the fuzzy logic system, specially to find the optimal design in the values for HS for harmony memory accepting (HMR) and pitch adjustment rate (PArate) to represent the optimized fuzzy harmony search (OFHS) and the dynamic parameters for mutation rate (F) and crossover rate (CR) for DE to represent the optimized fuzzy differential evolution (OFDE). The two proposed methods are applied for optimizing fuzzy controllers and mathematical functions, a statistical test is presented with the obtained results, and the main objective is to measure the efficiency in the performance of metaheuristic algorithms for the two study cases.
The organization of the sections in this paper is described in the following. Section 2 describes two metaheuristics implemented in this paper; Differential Evolution and Harmony Search Algorithms, Section 3 describes the optimal design of fuzzy system approach for the two study cases. Section 4 describes two study cases to analyze-the first is the stabilization of trajectory for an autonomous mobile robot (AMR), and the second is the optimization in mathematical functions. Section 5 describes a comparative of the results for the metaheuristic algorithms, and a statistical test is also presented. Section 6 presents some relevant conclusions and future works for this paper.

Meta-Heuristic Background
Metaheuristics are used for problem solving. This section describes the basic concepts of two metaheuristics that are used later in the paper, which have had great impact on the solution of different problems and are presented in more detail in the following subsections.

Differential Evolution Algorithm
In 1994, Storn and Price proposed the differential evolution (DE) algorithm [44], and the basic structure of the algorithm is based mainly on the following operations:

Population Structure
The DE contains Np D-dimensional vectors of real-valued parameters, where P x : The current population, x i,g : is composed of those vectors, and indices start with 0 to simplify working with arrays and modular arithmetic. The index, g = 0, 1, . . . , g max , indicates the generation to which a vector belongs. Once initialized, DE mutates randomly chosen vectors to produce an intermediary population P v,g of Np mutant vectors, v i,g . Each vector in the current population is recombined with a mutant vector to produce a trial population, P u , the NP, mutant vector u i,g . P x,g = x i,g , i = 0, 1, . . . , Np − 1, g = 0, 1, . . . , g max (1) Appl. Sci. 2020, 10, 6146 4 of 21

Initialization
The upper and lower limits of each parameter must be established, and these 2D values can be collected by two initialized vectors-D-dimensional, b L y b U , where the subscripts L and U indicate the lower and upper limits, respectively.
2.1.3. Mutation Three vectors are selected at random which multiplied by a scale factor F create the mutated vector, the scale factor, F ∈ (0, 1) is a positive real number that controls the rate at which the population evolves.

Crossover
A uniform crossover, also known as discrete (dual) recombination is used to search for differential mutations, where Cr ∈ (0, 1) is a positive real number:

Selection
In this operation, the test vector is analyzed, u i,g , in case it has a value equal to or less than its target vector, x i,g . The target vector is replaced in the next generation, otherwise the target retains its place in the population until the condition is met: The aforementioned operations are repeated until the stopping criterion is satisfied.

Harmony Search Algorithm
Harmony search (HS) was originally proposed by Zong Woo Geem in 2001 [45]. This algorithm is inspired by jazz music composition. It is made up of the following five operations: Step 1: Initialize the problem and parameters, LB and UB indicate the lower and upper limits, respectively: Step 2: Initialize the harmony memory (HM): Appl. Sci. 2020, 10, 6146 5 of 21 Step 3: Improvise a new harmony: improvisation is the most important part of the algorithm because the memory considerations (HMR), pitch adjustment (PAR), and randomization operations are generated, as indicated in the following equation: Step 4: Update the harmony memory: the objective function is used to evaluate the HM update with a new solution vector x_ (new). It is compared if the new vector solution is better than the worst solution of the historical vector, and then the worst historical solution is excluded and replaced by a new one.
Step 5: Check the stopping criteria: the aforementioned operations are repeated until the stop criterion is satisfied.

Optimal Design of Fuzzy Systems Approach
The proposed method is based on the concepts and equations of the original DE and HS algorithms. Both original algorithms have the problem of parameters being fixed throughout the iterations. In previous work, this problem was eliminated using fuzzy logic to perform the dynamic parameter adjustment. This method consisted of a Mamdani fuzzy system with two inputs (iterations and diversity) and two outputs for each method, as shown by Castillo et al. in [46]. Therefore, the main contribution in this paper is that the proposed method uses the original DE and HS algorithms to find the best architecture of the previous fuzzy differential evolution (FDE) and fuzzy harmony search (FHS) fuzzy systems shown by Castillo et al. in [46], which are responsible for adjusting parameters along the iterations. The proposed methods are optimized fuzzy differential evolution (OFDE) and optimized fuzzy harmony search algorithm (OFHS). Once the best architectures of the proposed OFDE and OFHS algorithms have been found, two types of problems are considered for optimization: benchmark mathematical functions and a non-linear control problem. Figure 1 shows the general diagram of the proposed method in this paper.

Optimal Design of Fuzzy Systems Approach
The proposed method is based on the concepts and equations of the original DE and HS algorithms. Both original algorithms have the problem of parameters being fixed throughout the iterations. In previous work, this problem was eliminated using fuzzy logic to perform the dynamic parameter adjustment. This method consisted of a Mamdani fuzzy system with two inputs (iterations and diversity) and two outputs for each method, as shown by Castillo et al. in [46]. Therefore, the main contribution in this paper is that the proposed method uses the original DE and HS algorithms to find the best architecture of the previous fuzzy differential evolution (FDE) and fuzzy harmony search (FHS) fuzzy systems shown by Castillo et al. in [46], which are responsible for adjusting parameters along the iterations. The proposed methods are optimized fuzzy differential evolution (OFDE) and optimized fuzzy harmony search algorithm (OFHS). Once the best architectures of the proposed OFDE and OFHS algorithms have been found, two types of problems are considered for optimization: benchmark mathematical functions and a non-linear control problem. Figure 1 shows the general diagram of the proposed method in this paper. The proposed fuzzy systems have triangular membership functions at the inputs and outputs in both algorithms. The triangular membership function depends on three parameters which are a, b, and c as given by Equation (14), where a and c represents the lower parts of the triangle and b locates the peak. This triangular membership function is represented graphically in Figure 2. The proposed fuzzy systems have triangular membership functions at the inputs and outputs in both algorithms. The triangular membership function depends on three parameters which are a, b, and c as given by Equation (14), where a and c represents the lower parts of the triangle and b locates the peak. This triangular membership function is represented graphically in Figure 2.
Appl. Sci. 2020, 10, x 6 of 23 Tables 1 and 2 summarize the representation of the knowledge of the inputs and outputs of the fuzzy system with triangular membership functions.      Low.
Figures 3 and 4 illustrate the parameters that form each input and output triangular membership function.
The fuzzy system was optimized with DE and HS; in the DE, it is necessary to use a chromosome, and in HS, a harmony to optimize the membership functions (MFS) is necessary, and in this case, it is called a vector, as shown in Figure 5. The vector has 44 positions, of which six parameters are fixed and 38 are optimized.  The fuzzy system was optimized with DE and HS; in the DE, it is necessary to use a chromosome, and in HS, a harmony to optimize the membership functions (MFS) is necessary, and in this case, it is called a vector, as shown in Figure 5. The vector has 44 positions, of which six parameters are fixed and 38 are optimized. Figure 5. Vector for the optimization optimized fuzzy differential evolution (OFDE) and optimal fuzzy harmony search (OFHS).   The fuzzy system was optimized with DE and HS; in the DE, it is necessary to use a chromosome, and in HS, a harmony to optimize the membership functions (MFS) is necessary, and in this case, it is called a vector, as shown in Figure 5. The vector has 44 positions, of which six parameters are fixed and 38 are optimized. Figure 5. Vector for the optimization optimized fuzzy differential evolution (OFDE) and optimal fuzzy harmony search (OFHS).      The limits for the optimization of each triangular membership function are specified in Table 3. Table 3. Limits of the parameters of the membership functions.

MFs Parameters
Low MF: High MF: High MF:

Study Cases
Two types of problems are used in this work in order to validate the proposed methodology. The first one is benchmark mathematical functions, and the second one is a benchmark control problem of the Sugeno type. The explanation of each of the problems is shown in more detail in the following subsections.

Benchmark Mathematical Functions
A set of eight benchmark mathematical functions are selected for the first part of the experimentation. The idea of using mathematical functions is to explore the behavior of our proposal. In Table 4 we show the search domain, minimum and their respective equation for each benchmark function.
The set of used functions is diverse, as we have simple functions and at the same time complicated functions, including both levels of complexity helps in our experimentation to check that it works in simple functions, as well as in complicated functions.

Robot Mobile Controller
The second control problem is a unicycle mobile robot, and a Sugeno type controller was used, which has a higher level of complexity. The main goal of this controller is to make the robot follow a reference trajectory. The robot is composed of two drive wheels located on the same axis and a front free wheel that is used only for stability. The graphical representation of the robot mobile is illustrated in Figure 6, in which we can observe the two wheels that make up the robot, the two torques, and the angle with respect to the Y axis.
Appl. Sci. 2020, 10, x 9 of 23 The set of used functions is diverse, as we have simple functions and at the same time complicated functions, including both levels of complexity helps in our experimentation to check that it works in simple functions, as well as in complicated functions.

Robot Mobile Controller
The second control problem is a unicycle mobile robot, and a Sugeno type controller was used, which has a higher level of complexity. The main goal of this controller is to make the robot follow a reference trajectory. The robot is composed of two drive wheels located on the same axis and a front free wheel that is used only for stability. The graphical representation of the robot mobile is illustrated in Figure 6, in which we can observe the two wheels that make up the robot, the two torques, and the angle with respect to the Y axis. Equations (15) and (16) represent the operation of the mathematical model of the mobile robot, and for each equation a brief explanation is given of each of the variables that constitute the equations is presented.
where q = (x, y, θ) T is the vector of the configuration coordinates; v = (v, ) T is the vector of velocities; τ = �τ 1, τ 2 � is the vector of torques applied to the wheels of the robot where 1 and 2 denote the torques of the right and left wheel, respectively; ∈ 2 is the uniformly bounded disturbance vector; ( ) ∈ 2 2 is the positive-definite inertia matrix; Equations (15) and (16) represent the operation of the mathematical model of the mobile robot, and for each equation a brief explanation is given of each of the variables that constitute the equations is presented. M(q) . v + C q, where q = (x, y, θ) T is the vector of the configuration coordinates; v = (v,) T is the vector of velocities; τ = (τ 1, τ 2 ) is the vector of torques applied to the wheels of the robot where τ 1 and τ 2 denote the torques of the right and left wheel, respectively; P ∈ R 2 is the uniformly bounded disturbance vector; M(q) ∈ R 2X2 is the positive-definite inertia matrix; q ϑ is the vector of centripetal and Coriolis forces; D ∈ R 2X2 is a diagonal positive-definite damping matrix. The kinematic system is determined by Equation (16): where (x,y) is the position in the X − Y (world) reference frame; θ is the angle between the heading direction and the x-axis; v and ω are the linear and angular velocities, respectively. Equation (17) represents the non-holonomic constraint, which this system has, which corresponds to a non-slip wheel condition preventing the robot from moving sideways.
One of the membership functions used in the controller is the trapezoidal one, which has four scalar parameters-a, b, c, and d as given by the Equation (18), where a and d represents the lower corners of the trapezoid and the b and c locate the shoulders. To represent the graphical form of the trapezoidal membership function, we illustrate it in Figure 7.
Equation (17) represents the non-holonomic constraint, which this system has, which corresponds to a non-slip wheel condition preventing the robot from moving sideways.
̇cos −̇sin = 0 (17) One of the membership functions used in the controller is the trapezoidal one, which has four scalar parameters-a, b, c, and d as given by the Equation (18)  The structure of the fuzzy system of the Sugeno-type controller is composed of two inputs: the first one is the linear velocity error, and the second is the angular velocity error, which are granulated into three membership functions (negative, zero, positive). Both are composed of trapezoidal functions on the extremes and a triangular one in the center, as can be appreciated in Figure 8. The structure of the fuzzy system of the Sugeno-type controller is composed of two inputs: the first one is the linear velocity error, and the second is the angular velocity error, which are granulated into three membership functions (negative, zero, positive). Both are composed of trapezoidal functions on the extremes and a triangular one in the center, as can be appreciated in Figure 8.
The structure of the fuzzy system of the Sugeno-type controller is composed of two inputs: the first one is the linear velocity error, and the second is the angular velocity error, which are granulated into three membership functions (negative, zero, positive). Both are composed of trapezoidal functions on the extremes and a triangular one in the center, as can be appreciated in Figure 8.   Tables 5 and 6 show the mathematical knowledge of each membership function for each of the inputs. Table 5. Mathematical representation of the input 1 membership function.

Angular Velocity Error (ew)
The controller used is of the Takagi-Sugeno type and is formed by nine fuzzy rules that govern the trajectory of the robot based on the controller. Table 7 contains more detail of these fuzzy rules. The Sugeno coefficients of this controller are presented in Table 8. Table 7. Rules for the fuzzy logic controller. The surface of this controller is shown in detail in Figure 9.
The proposed OFDE and OFHS methods are used to optimize the parameter values of the membership functions of the mobile robot controller in order to minimize the root mean square error (RMSE) described in Equation (19). Figure 10 represents the vector with the points of the membership functions of both controller inputs. It has 22 total points, where eight are fixed and 14 are optimized. These points that we call optimized are the points (parameter values) that the original algorithm provides in order to create a new structure of the fuzzy system. The limits for the optimization of each membership function are specified in Table 9. The proposed OFDE and OFHS methods are used to optimize the parameter values of the membership functions of the mobile robot controller in order to minimize the root mean square error (RMSE) described in Equation (19). Figure 10 represents the vector with the points of the membership functions of both controller inputs. It has 22 total points, where eight are fixed and 14 are optimized. These points that we call optimized are the points (parameter values) that the original algorithm provides in order to create a new structure of the fuzzy system.
The proposed OFDE and OFHS methods are used to optimize the parameter values of the membership functions of the mobile robot controller in order to minimize the root mean square error (RMSE) described in Equation (19). Figure 10 represents the vector with the points of the membership functions of both controller inputs. It has 22 total points, where eight are fixed and 14 are optimized. These points that we call optimized are the points (parameter values) that the original algorithm provides in order to create a new structure of the fuzzy system. The limits for the optimization of each membership function are specified in Table 9. The limits for the optimization of each membership function are specified in Table 9. Table 9. Limits of the parameters of the membership functions.

MFs Parameters
Negative MF:

Comparison of Results
Experiments were performed with the proposed OFDE and OFHS algorithms applied to benchmark mathematical functions and benchmark control cases to validate the proposed methodology. This section shows the experimentation and the obtained results.

Results for the Benchmark Mathematical Functions
Thirty experiments were performed with the proposed OFDE and OFHS methods applied to each of the eight benchmark mathematical functions shown in Section 4.1. The simulations were obtained as a result of software implemented in Matlab 2013 and 2015, from Mathworks, Los Angeles, CA, USA. The parameters used for the experiment for the OFDE method are shown in Table 10 and for the OFHS method are presented in Table 11.  Table 10 shows the values of the population, dimension, generation, run, mutation rate (F), and crossover rate (CR) parameters used in the experiment for the proposed OFDE. Table 11 shows the values of the harmony, dimension, iteration, run, harmony memory accepting (HMR), and pitch adjustment rate (PArate) parameters used in the experiment for the proposed OFHS.
The results obtained by using the proposed OFDE method are shown in Table 12, and the results obtained by using the proposed OFHS method are shown in Table 13. These tables indicate the execution time of each mathematical function, where we can note lower times for OFHS.  Tables 12 and 13 show the minimum and maximum obtained errors, the averages of the 30 experiments, and their standard deviations for all the mathematical functions.
A statistical z-test was performed to validate the efficacy of the proposed methods. The claim is that the OFDE method has achieved lower values compared to the OFHS method. The parameters of the statistical test are as follows: Alpha 0.05, confidence level 95%, sample size is 30 and the critical value is −1,645. where µ 1 = mean o f the OFDE algorithmµ 2 = mean o f the OFHS algorithm (20) The hypotheses would be as follows: H o : µ 1 ≥ µ 2 and H a : µ 1 < µ 2 (Claim). We can say that there is enough statistical evidence for all values below the critical value Z or −1.645. Table 14 shows the results of the mean, standard deviation, and the Z-values obtained; in three functions, significant evidence is obtained, and in five, no significant evidence is obtained. We can indicate that the OFDE method is not better than OFHS in these mathematical benchmark functions.

Experiments and Results for the Benchmark Control Problem
Thirty experiments were performed with the proposed OFDE and OFHS methods applied to mobile robot controller. The parameters used for the experiment with the OFDE method are shown in Table 10 and for the OFHS method are shown in Table 11.
The results obtained by optimizing the parameters of the mobile robot controller, with the OFDE and OFHS methods are presented in Tables 15 and 16, respectively. The noise applied to this controller is 0.5 (Gaussian random number). The execution times of each algorithm are also shown. We can summarize that, in both optimized algorithms, with our proposal, the improvement on average is observed, and in the same way, we can observe that the best result obtained is by using our new proposal. The comparison between the optimized method with noise and without noise also shows that our proposal improves on problems with uncertainty, and this is observed for both algorithms. Figures 11 and 12 represent the best simulations obtained with each of the methods; for Figure 11, the OFDE method with the best paths without noise and noise are illustrated, and Figure 12 illustrates the OFHS method with the best paths without noise and with noise. To corroborate the results from both algorithms, Table 17 is presented, which includes some results using different algorithms and different ways of optimizing applied to this control problem. Table 17 contains the minimum, maximum, average, standard deviation, and the number of experiments for other metaheuristics that were applied to the same problem.  To corroborate the results from both algorithms, Table 17 is presented, which includes some results using different algorithms and different ways of optimizing applied to this control problem. Table 17 contains the minimum, maximum, average, standard deviation, and the number of experiments for other metaheuristics that were applied to the same problem.  To corroborate the results from both algorithms, Table 17 is presented, which includes some results using different algorithms and different ways of optimizing applied to this control problem. Table 17 contains the minimum, maximum, average, standard deviation, and the number of experiments for other metaheuristics that were applied to the same problem.
In order to corroborate all the experimentation carried out with our proposal, three methods are chosen from Table 17 to perform statistical tests, and these methods are the only ones that meet the characteristics to perform a Z test.
The following explains how the statistical tests are applied: Table 18 shows the OFDE results with noise which represents our µ 1 and we perform a statistical test for each of the three methods (GAs optimization FLC Type-2, GA optimization Type-1, and GA optimization FLC Type-2), which is represented with µ 2 in each of the statistical tests.  The same dynamics is applied in Table 19, where three statistical tests are performed, but in a comparison with the OFHS results with noise. The claim is that the OFDE or OFHS method has lower (better) values compared to the GA optimization FLC Type-2 or GA optimization Type-1 or GA optimization FLC Type-2 method. The parameters of the statistical test are as follows: Alpha 0.05, confidence level 95%, sample size is 30 and the critical value is −1645. where µ 1 = OFDE or OFHS algorithm µ 2 = DA optimization FLC Type-2 or GA optimization Type-1 or GA optimization FLC Type-2 The hypotheses would be as follows: H o : µ 1 ≥ µ 2 and H a : µ 1 < µ 2 (Claim). In this case, we can say that there is enough evidence for all values below the critical value Z or −1.645.
Analyzing the results of the statistical tests for both methods used (OFDE and OFHS), we can note that, in all cases, our proposal has significant statistical evidence of being the better method when compared to the three methods of literature.

Conclusions
The main conclusions in this paper consist of highlighting the minimization of errors when the original algorithms are modified using the techniques of fuzzy logic. In both study cases in this paper, the proposed method demonstrates a stabilization and improvement with respect to the original algorithms. In this case, the proposed modification statistically shows in the mathematical functions a minimization of the error and, for the second problem of the mobile robot controller, stabilization in the trajectory even when noise is added in the model.
A comparison was presented between the two proposed methods in the case studies. In the first case of benchmark mathematical functions, OFHS achieves better results than OFDE. In the second case of control, both methods achieve good results by optimizing the robot controller with and without noise; in both cases, it is possible to minimize the error by applying noise to the controller. By comparing the results obtained with other existing methods in the literature, the effectiveness of the proposed methodology can be validated. The results obtained show the efficiency of the fuzzy logic system to find the optimal design in the membership functions (MFs) that allows the user to identify the values of the HMR and PArate parameters for the HS algorithm and the values of the F and CR parameters for the DE algorithm, and these values allow the user to obtain the modified algorithms that demonstrate a better performance in error minimization for the two study cases.
As to future work, we envision the implementation of this proposed method using interval type-2 fuzzy logic systems, as well as the experimentation with other non-linear models and being able to visualize the efficiency that this proposed method can demonstrate on the fuzzy control field.