Comparison between Genetic Algorithms of Proportional– Integral–Derivative and Linear Quadratic Regulator Controllers, and Fuzzy Logic Controllers for Cruise Control System

: One of the most significant and widely used features currently in autonomous vehicles is the cruise control system that not only deals with constant vehicle velocities but also aims to optimize the safety and comfortability of drivers and passengers. The accuracy and precision of system responses are responsible for cruise control system efficiency via control techniques and algorithms. This study presents the dynamic cruise control system model, then investigates a genetic algorithm of the proportional–integral–derivative (PID) controller with the linear quadratic regulator (LQR) based on four fitness functions, the mean squared error (MSE), the integral squared error (ISE), the integral time squared error (ITSE) and the integral time absolute error (ITAE). Then, the response of the two controllers, PID and LQR, with the genetic algorithm was compared to the response performance of the fuzzy and fuzzy integral (Fuzzy-I) controllers. The MATLAB 2024a program simulation was employed to represent the system time response of each proposed controller. The output simulation of these controllers shows that the type of system stability response was related to the type of controller implemented. The results show that the Fuzzy-I controller outperforms the other proposed controllers according to the least Jmin function, which represents the minimum summation of the overshoot, settling time, and steady-state error of the cruise control system. This study demonstrates the effectiveness of driving accuracy, safety, and comfortability during acceleration and deceleration due to the smoothness and stability of the Fuzzy-I controller with a settling time of 5.232 s and when converging the steady-state error to zero.


Introduction
In recent years, autonomous vehicles have made great strides in progress and development in the field of transportation at the individual and societal levels, especially in terms of comfort, safety, energy consumption, and maneuvering accuracy [1][2][3][4][5].One of the most significant contributing factors that played a prominent role in the stages of development of autonomous vehicles is cruise control systems, which not only originally contributed to the stability of the vehicle speed but were also incorporated into many features in the gradual evaluation of control in autonomous vehicles, such as modifying speed automatically to maintain a safe distance between vehicles to avoid collision [6][7][8][9].To achieve the optimal and most accurate performance of vehicles, the cruise control system uses various control techniques and methods such as the proportional-integral-derivative (PID) control, fuzzy logic control (FLC), the linear quadratic regulator (LQR), and model predictive control (MPC) [10][11][12][13][14][15].
Many studies related to the cruise control system have been published and utilized various control methods.Pradhan et al. proposed a PID controller to improve the stability of the cruise control system based on the Ant Lion algorithm to optimize the PID parameters [16].They showed better performance compared to the other proposed control approaches.Similarly, Wu et al. proposed a PID controller with the chaotic ant swarm algorithm to enhance the time response for an automobile cruise system [17].They demonstrated increased driving comfortability and safety as well as decreased fuel consumption.Singhal et al., in [18], applied the GNA (global neighborhood algorithm) to specify PID controller parameters to mitigate the system response error that also led to increased driving comfortability.Other researchers also adopted a PID controller by applying different scenarios to increase control parameter optimization such as particle swarm optimization, the genetic algorithm, Fractional Order PID (FOPID), and the model predictive controller (MPC) [19][20][21][22].
From another point of view, several publication studies focused on using a fuzzy logic controller for an autonomous vehicle control system.Abdurohman et al. utilized a fuzzy control system to select the optimal control parameters of a PID controller [23].Their proposed method showed the safety, comfortability, and stability of the cruise control system by adjusting the accelerator pedal under different drive conditions.On the other hand, Mao et al. proposed a fuzzy model predictive control algorithm (Fuzzy-MPC) to achieve better performance than the traditional algorithm under different road conditions by considering the safety, comfortability, and stability of the cruise control system [24].Likewise, Shojaeefard et al. presented fuzzy logic control to avoid the collision between vehicles by modifying the vehicle speed and acceleration dynamically to meet the minimum requirement distance based on changing weather conditions [25].From a medical standpoint, Yang et al. introduced fuzzy control to analyze driver comfortability according to the driver's heat rate (HRV) [26].In order to decrease driver discomfort physiologically related to the acceleration of braking, they extracted electrocardiography (ECG) sensor data, and stabilized the cruise control system by applying fuzzy control.
In addition to the controller methods mentioned above, some researchers have used the linear quadratic regulator (LQR) for a more complex cruise control system.Varma et al. designed a LQR controller for autonomous vehicle path tracking in consideration to roll and yaw stability control, where most car accidents occur [27].Li et al. applied a variable weight of LQR to minimize the fluctuation of longitudinal acceleration that keeps cars in line and improves driving safety comparing to constant weight [28].From another point of view, Zhang et al. adopted the Fuzzy-LQR method to design a cruise control system to modify the accelerating speed during lane change to avoid vehicle collision and ensure driving comfort and safety [29].In this paper, three different controllers for the cruise control system were presented with genetic algorithm (GA) optimization to obtain the appropriate control parameters.A performance analysis comparison of the PID, FLC, and LQR controllers was conducted and discussed in terms of using various fitness functions, which make the system more accurate and stable using the MATLAB simulation program during acceleration and deceleration (mean squared error (MSE), integral squared error (ISE), integral time squared error (ITSE), and integral time absolute error (ITAE)).
This paper is divided in to five sections as follows: the first section includes the background and a literature review in the introduction.In the second section, the linearization of the cruise control system model and mathematical expression are derived.In section three, the three proposed controllers are designed using a genetic algorithm and fitness functions.Then, the simulation and results section is presented using the MATLAB program to show the best performance among controllers on applying various fitness functions.Finally, a performance system analysis of the proposed controllers will be summarized in the conclusion.

Model of the Cruise Control System
The main concern of using the cruise control system is the stability of vehicle speed according to driver input instruction.Rather than using the accelerator or brake pedal to maintain a constant vehicle speed, the cruise control system utilizes the feedback signals from the speed sensor to converge or match the speed reference (v ref ) set by the driver due to modifications to the gate angle of engine throttle actuator (u), and then specifies the appropriate driving force of the engine (F d ).The block diagram of the cruise control system is represented by Figure 1 [30].

Model of the Cruise Control System
The main concern of using the cruise control system is the stability of vehicle speed according to driver input instruction.Rather than using the accelerator or brake pedal to maintain a constant vehicle speed, the cruise control system utilizes the feedback signals from the speed sensor to converge or match the speed reference (vref) set by the driver due to modifications to the gate angle of engine throttle actuator (u), and then specifies the appropriate driving force of the engine (Fd).The block diagram of the cruise control system is represented by Figure 1 [30].By taking into consideration the aerodynamic drag force (Fa), the climbing resistance force (Fg) and, according to Newton's second law, the nonlinear longitudinal dynamic of the system can be represented by the following equation [31]: where M represents the overall passengers and vehicle mass, while     is the inertia force.All forces in Equation (1) can be found in the following expression [32]: =  sin  (3) where Ca is the coefficient of aerodynamic drag, v and vw are the actual vehicle and wind speed, respectively, g is the gravity acceleration, and θ is the incline angle of road.Equation (4) has the maximum and minimum limitations of brake and acceleration forces, respectively, where C1 is the value added to the constant of the actuator force saturation,  − is the delay factor, τ is the reaction time of driver, and Ts is a time constant.Figure 2 shows the above three force equations in a block diagram of vehicle dynamics [31].
To simplify the control equations of the system shown in Figure 2, a value of zero will be set to all the initial condition due to disturbances that eliminated the wind and grading effects.Due to the simplification assumption, Equation (1) should be rearranged into the following equation: By taking into consideration the aerodynamic drag force (F a ), the climbing resistance force (F g ) and, according to Newton's second law, the nonlinear longitudinal dynamic of the system can be represented by the following equation [31]: where M represents the overall passengers and vehicle mass, while M d dt v is the inertia force.All forces in Equation (1) can be found in the following expression [32]: where C a is the coefficient of aerodynamic drag, v and v w are the actual vehicle and wind speed, respectively, g is the gravity acceleration, and θ is the incline angle of road.Equation (4) has the maximum and minimum limitations of brake and acceleration forces, respectively, where C 1 is the value added to the constant of the actuator force saturation, e −τs is the delay factor, τ is the reaction time of driver, and T s is a time constant.Figure 2 shows the above three force equations in a block diagram of vehicle dynamics [31].By applying inverse Laplace for Equation (4), Equation ( 1) is rewritten into the following equation: To simplify the control equations of the system shown in Figure 2, a value of zero will be set to all the initial condition due to disturbances that eliminated the wind and grading effects.Due to the simplification assumption, Equation (1) should be rearranged into the following equation: By applying inverse Laplace for Equation (4), Equation ( 1) is rewritten into the following equation: To overcome the nonlinearity of Equation ( 5), Osman et al. differentiated both sides of the equation [31], whereas Izci et al. used Taylor's series expansion [30].Thus, the transfer function of vehicle speed output related to the pedal actuator input of Equation ( 6) can be written as follows [30,31]: By substituting the parameters and specification values in Table 1 [30,31], the linearized transfer function of the derived mathematical expression of the system can be found in Equation ( 8), while Equations ( 9) and ( 10) represent the state space of the system used by the proposed controller, as outlined in the next section.

Controller Design
In this section, the analysis and design of the presented controllers will be presented and discussed.This section is divided into two parts to ensure a comprehensive understanding of the methodologies employed.The first part will focus on the design and implementation of the proportional-integral-derivative (PID) and linear quadratic regulator (LQR) controllers.These controllers will be optimized using the genetic algorithm (GA) optimization method.This part will delve into the theoretical foundations of the PID and LQR controllers, followed by a detailed explanation of how the genetic algorithm is applied to enhance their performance and efficiency.The second part will focus on the design and development of the fuzzy logic controller and the fuzzy integral controller.This section will cover the principles of fuzzy logic, explaining its application in control systems and how it can be integrated with traditional control methods to improve system responsiveness and robustness.The fuzzy integral controller will be introduced as an advanced technique that combines fuzzy logic with integral control to achieve superior control outcomes.

PID and LQR Controllers with the GA Optimization Method
The genetic algorithm (GA) method is considered one of the most prominent and efficacious optimization techniques utilized across numerous domains, including control system design.It has been extensively employed to ascertain the optimal values for the parameters of controllers, serving as a critical step in the controller design process.This method belongs to the category of evolutionary algorithms, which are inspired by the principles of natural selection and genetics.By emulating the process of natural selection, genetic algorithms iteratively evolve solutions to optimization problems, ultimately converging on the most optimal solutions.
In the context of controller design, the genetic algorithm plays an integral role in augmenting the performance and efficiency of control systems.It commences with the population of potential solutions, each represented as a set of chromosomes or individuals.These individuals are evaluated based on a predefined fitness function, which assesses the extent to which each solution satisfies the desired control objectives.The algorithm then applies genetic operators, such as selection, crossover, and mutation, to generate a new population of potential solutions.Through selection, the most proficient individuals are chosen to transmit their attributes to subsequent generations, analogous to the concept of survival of the fittest.Crossover enables the amalgamation of characteristics from disparate individuals, facilitating the creation of potentially superior offspring.Mutation introduces stochastic variations in the offspring's chromosomes, ensuring genetic diversity and aiding the algorithm in evading local optima.By iteratively applying these operators over multiple generations, the genetic algorithm incrementally refines the solutions, steering them toward optimal parameter values for the controllers.The outcome is a set of optimized controller parameters that enhance the overall performance, stability, and robustness of the control system.
To use GA, firstly, a population size that represents the potential solutions of the optimal problem should be selected; this population will be represented as chromosome structures, which are typically a binary string.This chromosome structure includes all the controller's elements that must be optimized.In this paper, it will be used for the PID and LQR controllers.For the PID controller, it is required to optimize the values of the K p , K i , and K d , while for the LQR controller, it is required to optimize the Q and R matrices.
After forming the chromosome, the next step is the use of the fitness function in order to obtain the performance of each one.The next step is the crossover; this is between a pair of individuals, so the information between the two chromosomes will be exchanged to create totally new chromosomes.After crossover, mutation will occur randomly to only some of the individuals in order to increase the diversity in the population.
Figure 3 shows the flowchart of the GA.
In the domain of control system design utilizing genetic algorithms, the selection of an appropriate fitness function is of the utmost importance.The fitness function serves as a quantitative metric that evaluates the performance of a given set of controller parameters concerning the specified objectives.The choice of fitness function influences the optimization process's outcome, as different functions can produce different results.Thus, identifying the most suitable fitness function for the model in use is crucial for achieving optimal performance where each fitness function has its strengths and weaknesses, and their applicability is contingent upon the specific control objectives and evaluation criteria.In the domain of control system design utilizing genetic algorithms, the selection of an appropriate fitness function is of the utmost importance.The fitness function serves as a quantitative metric that evaluates the performance of a given set of controller parameters concerning the specified objectives.The choice of fitness function influences the optimization process's outcome, as different functions can produce different results.Thus, identifying the most suitable fitness function for the model in use is crucial for achieving optimal performance where each fitness function has its strengths and weaknesses, and their applicability is contingent upon the specific control objectives and evaluation criteria.
Commonly utilized fitness functions include mean squared error (MSE), integral squared error (ISE), integral time squared error (ITSE), and integral time absolute error (ITAE).Each of these functions is designed to address specific control evaluation scenarios by emphasizing performance aspects such as the transient response, steady-state response, or both.
The mean squared error (MSE) function is frequently employed due to its straightforward nature and efficacy in minimizing the average of the squared deviations between the desired and actual system outputs.This function is particularly beneficial when the primary objective is to reduce the overall deviation in the system's response from the desired trajectory.Integral squared error (ISE), on the other hand, integrates the square of the error over time, thereby placing emphasis on mitigating substantial deviations throughout the entire response period.This characteristic renders ISE suitable for systems minimizing both transient and steady-state errors.Integral time squared error (ITSE) further refines this approach by weighting errors that occur later in the response, making it advantageous for applications where the transient response is critically important and later-stage precision is essential.At the same time, the integral time absolute error (ITAE) emphasizes the minimization of time-weighted absolute errors, prioritizing rapid error correction and minimal overshoot.This makes ITAE particularly well-suited for scenarios where a balance between transient and steady-state performance is desirable.
The equations of these fitness functions are represented in Equations (11)(12)(13)(14).Commonly utilized fitness functions include mean squared error (MSE), integral squared error (ISE), integral time squared error (ITSE), and integral time absolute error (ITAE).Each of these functions is designed to address specific control evaluation scenarios by emphasizing performance aspects such as the transient response, steady-state response, or both.
The mean squared error (MSE) function is frequently employed due to its straightforward nature and efficacy in minimizing the average of the squared deviations between the desired and actual system outputs.This function is particularly beneficial when the primary objective is to reduce the overall deviation in the system's response from the desired trajectory.Integral squared error (ISE), on the other hand, integrates the square of the error over time, thereby placing emphasis on mitigating substantial deviations throughout the entire response period.This characteristic renders ISE suitable for systems minimizing both transient and steady-state errors.Integral time squared error (ITSE) further refines this approach by weighting errors that occur later in the response, making it advantageous for applications where the transient response is critically important and later-stage precision is essential.At the same time, the integral time absolute error (ITAE) emphasizes the minimization of time-weighted absolute errors, prioritizing rapid error correction and minimal overshoot.This makes ITAE particularly well-suited for scenarios where a balance between transient and steady-state performance is desirable.
The equations of these fitness functions are represented in Equations ( 11)- (14).
To use the GA with the PID controller, it is crucial to understand the parts of the PID controller and their effect on the response.It consists of three parts (proportional, integral, and derivative), and the sum of their effect forms the controller output.When increasing the proportional gain K p , the response of the system will be faster, but it will increase the overshoot of the response while increasing the derivative gain K d and reduce the overshoot and any sudden changes in the response.Finally, the use of the integral part will eliminate the steady-state error, even though adding it will have a negative effect on the transient response.
Equation (15) shows the representation of the control single of the PID controller, while Figure 4 shows the block diagram of the GA with the PID controller.To use the GA with the PID controller, it is crucial to understand the parts of controller and their effect on the response.It consists of three parts (proportional, and derivative), and the sum of their effect forms the controller output.When in the proportional gain Kp, the response of the system will be faster, but it will incr overshoot of the response while increasing the derivative gain Kd and reduce the o and any sudden changes in the response.Finally, the use of the integral part will e the steady-state error, even though adding it will have a negative effect on the t response.
Equation (15) shows the representation of the control single of the PID co while Figure 4 shows the block diagram of the GA with the PID controller.For the LQR controller, the optimization of the weighting matrices (Q and R to minimize the quadratic cost function where they play a crucial role in detecting of each state and control input in the controller, as the Q matrix is the state-cos matrix, while the R matrix is the control weight matrix.

𝑢(𝑡) = 𝐾 𝑝 . 𝑒(𝑡) + 𝐾 𝑖 ∫ 𝑒(𝑡)𝑑𝑡 + 𝐾 𝑑 𝑑𝑒(𝑡) 𝑑𝑡
Figure 5 shows the block diagram of the LQR controller which is a state f control technique that uses the quadratic cost function to determine the optima input to minimize the performance index of the system.The role of the feedbac to measure the state variables of the system to calculate the optimal control inpu For the LQR controller, the optimization of the weighting matrices (Q and R) is used minimize the quadratic cost function where they play a crucial role in detecting the role of each state and control input in the controller, as the Q matrix is the state-cost weight matrix, while the R matrix is the control weight matrix.
Figure 5 shows the block diagram of the LQR controller which is a state feedback control technique that uses the quadratic cost function to determine the optimal control input to minimize the performance index of the system.The role of the feedback loop is to measure the state variables of the system to calculate the optimal control inputs.As the LQR controller used the state space representation shown the quadratic cost function to be minimized is given by Equation ( 14).

𝑥 ̇(𝑡) = 𝐴 𝑥(𝑡) + 𝐵 𝑢(𝑡)
Where x(t) is the state vector, u(t) is the control input vector, A is the state B is the input matrix.As the LQR controller used the state space representation shown in Equation ( 16), the quadratic cost function to be minimized is given by Equation ( 14).
x (t) = A x(t) + B u(t) (16) where x(t) is the state vector, u(t) is the control input vector, A is the state matrix, and B is the input matrix.
where Q is the state cost matrix and R is the control cost matrix.By using the optimal values of the Q and R matrices, the controller gain K can be computed using Equations (18)(19)(20).
where P is the solution to the continuous-time algebraic Riccati equation.
Finally, the block diagram represents the LQR controller with the GA is shown in Figure 6.
B is the input matrix.
where Q is the state cost matrix and R is the control cost matrix.By using values of the Q and R matrices, the controller gain K can be computed using (18)(19)(20).
+  −  −1    +  = 0 Where P is the solution to the continuous-time algebraic Riccati equation.
Finally, the block diagram represents the LQR controller with the GA i Figure 6.

Fuzzy Logic Controller
The fuzzy logic controller (FLC) is a knowledge-based controller that use logic principle and utilizes linguistic variables and fuzzy rules to represent complex, nonlinear systems.For this reason, it has been widely applied, esp systems that have uncertainty and nonlinearity.One drawback of the FLC is sponse still has steady-state errors, even though it diminishes the overshoot.T

Fuzzy Logic Controller
The fuzzy logic controller (FLC) is a knowledge-based controller that uses the fuzzy logic principle and utilizes linguistic variables and fuzzy rules to represent and control complex, nonlinear systems.For this reason, it has been widely applied, especially with systems that have uncertainty and nonlinearity.One drawback of the FLC is that its response still has steady-state errors, even though it diminishes the overshoot.To overcome this problem, an integral part (I) will be added to design a Fuzzy-I controller.The FLC has two inputs, proportional error PE and the change in error CE, and one output, which is the output of the controller U. Table 2 list the linguistic rules that will be used for the fuzzy and Fuzzy-I controllers where both inputs and the output have the following seven membership functions:

•
Negative big (NB);  The method that is used to determine the fuzzy characteristics of the output is called the fuzzy inference, which takes into account every rule in the rule base for the fuzzy inputs.The next step is finding the crisp control output using the defuzzification process.Figure 7 shows the membership functions of the inputs and the output.It is important to use trapezoidal functions for both ends of the inputs and the output to cover all extreme values while using gaussian functions for a smooth change in membership.The method that is used to determine the fuzzy characteristics of the the fuzzy inference, which takes into account every rule in the rule ba inputs.The next step is finding the crisp control output using the defuzz Figure 7 shows the membership functions of the inputs and the output.I use trapezoidal functions for both ends of the inputs and the output to c values while using gaussian functions for a smooth change in membersh

Simulation and Results
In this section, firstly, the PID and LQR controllers will be run ind four fitness functions to choose the controller with the best fitness MATLAB Simulink will be used to compare the results for all controlle the PID controller gain values using the GA optimization method for the u tions shown in Table 3, the first comparison is made between the PID different four fitness functions, as shown in Figure 8.

Simulation and Results
In this section, firstly, the PID and LQR controllers will be run individually for the four fitness functions to choose the controller with the best fitness function.Then, MATLAB Simulink will be used to compare the results for all controllers.After finding the PID controller gain values using the GA optimization method for the used fitness functions shown in Table 3, the first comparison is made between the PID controllers with different four fitness functions, as shown in Figure 8.At the same time, the same test is conducted to the LQR; the values of the K matrix are shown in Table 4, and the step response is shown in Figure 9.At the same time, the same test is conducted to the LQR; the values of the K matrix are shown in Table 4, and the step response is shown in Figure 9.As the response specifications for both controllers have very similar responses, the Jmin equation was used to find the controller with the best fitness function by comparing the transient and steady-state response specifications, which are the overshoot percent As the response specifications for both controllers have very similar responses, the Jmin equation was used to find the controller with the best fitness function by comparing the transient and steady-state response specifications, which are the overshoot percent (%OS), rise time (Tr), and settling time with a threshold of 2% (Ts), as shown in Table 5.As the PID-ITSE and LQR-ITAE have the minimum Jmin value, they will be used in the MATLAB Simulink to be compared with the fuzzy and Fuzzy-I controllers.The block diagram for all controllers with the model of the cruise control system is shown in Figure 10.  Figure 12 shows the time response in the case of the acceleration of speed from 60 km/h (16.667 m/s) to 90 km/h (25 m/s) and the deceleration of speed back to 60 km/h.
To compare the time response specifications of the designed controllers, their values are shown in Table 6.It is important to use the Jmin function to compare the responses.Step response of PID, LQR, fuzzy, and Fuzzy-I controllers with speed acceleration and deceleration.
To compare the time response specifications of the designed controllers, their values are shown in Table 6.It is important to use the Jmin function to compare the responses.As shown in Table 6, the Fuzzy-I controller has the smallest Jmin value, and it be concluded that it is the best controller for the cruise control system.Moreover, the fuzzy controller has more advantages over the LQR controller, like its ability to handle nonlinearity, robustness, ease of implementation, flexibility, adaptability, and no requirements for precise modeling.

Conclusions
Based on our extensive research, we concluded that the Fuzzy-I controller has the best time response compared to the PID-GA, LQR, fuzzy, and Fuzzy-I controllers for cruise control systems.Before reaching this conclusion, firstly, the performance of the PID-GA and LQR-GA were tested separately using four different fitness functions, which were the MSE, ISE, ITSE, and ITAE.Then, the best controller with the best fitness functions was chosen for a later comparison with the fuzzy and Fuzzy-I controllers, where the best con-  As shown in Table 6, the Fuzzy-I controller has the smallest Jmin value, and it be concluded that it is the best controller for the cruise control system.Moreover, the fuzzy controller has more advantages over the LQR controller, like its ability to handle nonlinearity, robustness, ease of implementation, flexibility, adaptability, and no requirements for precise modeling.

Conclusions
Based on our extensive research, we concluded that the Fuzzy-I controller has the best time response compared to the PID-GA, LQR, fuzzy, and Fuzzy-I controllers for cruise control systems.Before reaching this conclusion, firstly, the performance of the PID-GA and LQR-GA were tested separately using four different fitness functions, which were the MSE, ISE, ITSE, and ITAE.Then, the best controller with the best fitness functions was chosen for a later comparison with the fuzzy and Fuzzy-I controllers, where the best controller for the PID controller was with the ITSE fitness function, while for the LQR controller, it was ITAE.Secondly, the responses of the PID-ITSE and LQR-ITAE were compared with the fuzzy and Fuzzy-I controllers.As it was difficult to determine which controller had the best response, the Jmin function was used to find it; the Fuzzy-I has the smallest Jmin, so it is the best controller for the cruise control system.

Figure 1 .
Figure 1.Block diagram of cruise control system.

Figure 1 .
Figure 1.Block diagram of cruise control system.

World 15 Figure 2 .
Figure 2. Block diagram of vehicle dynamics.

Figure 2 .
Figure 2. Block diagram of vehicle dynamics.

Figure 5 .
Figure 5. Block diagram of the LQR controller.

Figure 5 .
Figure 5. Block diagram of the LQR controller.

Figure 7 .
Figure 7. Input and output membership functions.

Figure 7 .
Figure 7. Input and output membership functions.

Figure 8 .
Figure 8.Time response of the PID controller for all fitness functions.

Figure 9 .
Figure 9.Time response of the LQR controller for all fitness functions.

Figure 9 .
Figure 9.Time response of the LQR controller for all fitness functions.

15 Figure 10 .
Figure 10.Block diagram of the cruise control system with controllers.

Figure 10 .
Figure 10.Block diagram of the cruise control system with controllers.The time response of the PID-ITSE, LQR-ITAE, fuzzy, and Fuzzy-I controllers are shown in Figure 11.Figure12shows the time response in the case of the acceleration of speed from 60 km/h (16.667 m/s) to 90 km/h (25 m/s) and the deceleration of speed back to 60 km/h.To compare the time response specifications of the designed controllers, their values are shown in Table6.It is important to use the Jmin function to compare the responses.

Figure 10 .
Figure 10.Block diagram of the cruise control system with controllers.The time response of the PID-ITSE, LQR-ITAE, fuzzy, and Fuzzy-I controllers are shown in Figure 11.

Figure 12
Figure12shows the time response in the case of the acceleration of speed from 60 km/h (16.667 m/s) to 90 km/h (25 m/s) and the deceleration of speed back to 60 km/h.

Figure 12 .
Figure 12.Step response of PID, LQR, fuzzy, and Fuzzy-I controllers with speed acceleration and deceleration.

Author Contributions:
The authors have the following contribution.A.M.: controller design, simulation and results, conclusion, writing-original draft preparation, software.K.Y.A.A.-b.: writingoriginal draft preparation, Abstract, Literature review, Introduction, Methodology, References.R.S.: supervision, investigation, reviewing and editing.All authors have read and agreed to the published version of the manuscript.

Figure 8 .
Time response of the PID controller for all fitness functions.

Table 5 .
Time response specifications of PID and LQR controllers for all fitness functions.

Table 6 .
Controller time response specifications.

Table 6 .
Controller time response specifications.