New Fusion Algorithm-Reinforced Pilot Control for an Agricultural Tricopter UAV

: Currently, fuzzy proportional integral derivative (PID) controller schemes, which include simpliﬁed fuzzy reasoning decision methodologies and PID parameters, are broadly and e ﬃ caciously practiced in various ﬁelds from industrial applications, military service, to rescue operations, civilian information and also horticultural observation and agricultural surveillance. A fusion particle swarm optimization (PSO)–evolutionary programming (EP) algorithm, which is an improved version of the stochastic optimization strategy PSO, was presented for designing and optimizing controller gains in this study. The mathematical calculations of this study include the reproduction of EP with PSO. By minimizing the integral of the absolute error (IAE) criterion that is used for estimating the system response as a ﬁtness function, the obtained integrated design of the fusion PSO–EP algorithm generated and updated the new elite parameters for proposed controller schemes. This progression was used for the complicated non-linear systems of the attitude-control pilot models of a tricopter unmanned aerial vehicle (UAV) to demonstrate an improvement on the performance in terms of rapid response, precision, reliability, and


Introduction
Evolutionary programming (EP) is one of the four primary evolutionary algorithm paradigms and was first proposed by Lawrence J. Fogel in the United States in 1966 [1]. The development of EP, which is parallel to global search techniques, was continued by his son D. B. Fogel. In contrast to the genetic algorithm (GA), the EP algorithm does not include a crossover operator but only contains a mutation operator and a specific selection strategy; however, its parameters are allowed to evolve. In EP progress, the prior knowledge of dynamics system or environmental information is not necessary. Those optimization methodologies are alike in terms of their natural evolution theory and differ mostly in terms of their parameter representation approach.
Particle swarm optimization (PSO) was proposed by Kennedy and Eberhart [2] to solve optimization problems, and belongs to the class of swarm intelligence methods. This algorithm is based on stochastic and population-based searches and tries to imitate the behavior of a flock of birds and shoal of fish. The original PSO method only has a few adjusted parameters and has higher computation speed, higher accuracy, and smaller memory size compared with the methods such as neural networks, genetic computation, and machine learning. However, a method for determining the proper values of PSO parameters at the optimal level is yet to be found. This challenging problem has received considerable attention in [3][4][5]. Hence, at present, in order to solve the local optima trap and improve the optimization process, many different methodologies of hybrid PSO algorithms have been enhanced to overcome this downside. For example, Juang proposed a hybrid PSO algorithm known as the hybrid version of a genetic algorithm (GA) with PSO [6], which incorporates the elite strategy for evolutionary operations of a GA. In this strategy, the upper half of optimal-performance individuals in a population is regarded as elites. Bergh and Engelbrecht [7] introduced a cooperation of PSO. Huang and Mohan [8] proposed a hybrid condition to obtain PSO robustness. To estimate the parameters of non-linear systems, Li et al. [9] presented the combination of PSO and simulated annealing for the jumping property. Next, Araújo and Uturbey [10] successfully applied a hybrid PSO algorithm to the generation and demand of the dispatch problem. Later, Ali and Tawhid solved many engineering optimization problems by using a hybrid PSO-differential evolution method in [11]. Recently, hybrid PSO algorithms have been exploited to solve unmanned aerial vehicle (UAV) path planning in [12,13]. The combination of PSO and support vector regression (SVR) for tourist arrivals forecasting is utilized by [14]. Hybrid PSO and the direct search method (DSM) for optimization of large-scale economic dispatch with valve-point effects is exploited in [15]. The hybrid swarm intelligence approach for blog success prediction is developed in [16]. All the aforementioned techniques can improve the parameters evolution speed. The use of fusion and hybrid algorithms has been increasing in recent decades because these methods combine the advantages of different algorithms and also eliminate the drawbacks of single optimization algorithms. Therefore, this study employs the fusion PSO-EP algorithm to reinforce the pilot control of an agricultural tricopter UAV.
Due to the effortless control formation, ease of design, and low cost of traditional proportional integral derivative (PID) controllers, they have been currently employed in real-world systems. Moreover, a novel hybrid control approach involving a fuzzy-PID controller was successfully used for an autonomous mini helicopter [17]. The model of fuzzy logic rule is useful when the physical progression of obtaining information is ambiguous, and limited data is available. The use of a fuzzy logic control (FLC) is proven as one of the most efficient and systematic approaches for solving non-linear practical systems by emulating human logic rather than using an accurate mathematical model.
The membership functions of fuzzy controller design are independent of extracted rules. Thus, a high system performance cannot be achieved, especially when there are many input variables in a complex system. Nevertheless, in most of circumstances, the fuzzification and de-fuzzification and also the membership functions could be adjusted to improve the performance of systems [18][19][20][21].
UAVs, which are able to autonomously operate, are considerably complicated non-linear systems with uncertainties. The operating conditions of UAVs are also unpredictable, for example, changes in direction and speed of wind, results in variations in the flight attitude, direction, speed and precision. Therefore, the control approaches of an aerial robot, especially in attitude controller design issues [22][23][24][25][26][27][28][29][30][31][32][33] find many difficulties to provide high performance, absolute precision and robustness. Thus, a tricopter, which can be successfully applied in horticultural observation and agricultural surveillance [34][35][36][37][38] is selected as a best model to present the efficiency of proposed controller design.
In general, this study proposed the combination of a fuzzy controller that design by the triangular membership functions and a conventional PID controller to compose the fuzzy-PID controller scheme. Then, the proposed controllers utilize the advantage of the fusion evolutionary PSO algorithm by using a fitness function associated with a performance index of the system-the integral of absolute error (IAE) [39,40]-to optimize the controller parameters. The designed controllers demonstrate the benefits of the improving performance of Tricopter pilot control models that could be fruitfully operated on agriculture fields [41][42][43][44].

Fusion Particle Swarm Optimization-Evolutionary Programming (PSO-EP) Algorithm-Reinforced Precision Control Design
In this section, we proposed the optimization algorithm PSO-EP for tuning the five parameters of the fuzzy-PID controller that apply to the tricopter pilot control. This simultaneous tuning task is complicated since the final trade-off in multi-objective optimization [10] must find out as fast as possible, the fusion PSO and EP algorithm is employed to save design efforts, operation time, and also obtain better performance.

Fusion PSO-EP Algorithm
The fundamental PSO technique includes a number of particles moving in a multi-dimensional searching space. The PSO algorithm is initiated on a random population and searches the optimal point by continually updating this population. The PSO algorithm improvement are as follows: (1) the evolution operators-crossover and mutation-are involved and (2) the adjustment of various free parameters is not required. Each potential solution, represented by the set of particles, is assigned a random velocity to be "flown" in the problem area. The current data can be shared with related particles at the optimal fitness time. The position and velocity vectors of each particle are remembered and updated at the moment at which the vectors reached the highest fitness level through iterations. The particle-birds that fly through the searching swarm space are being attracted toward the optimal solution.
Assume that the searching space has a dimension D and m particles. The ith particle settle at position , and the "flying" particle velocity is presented by . Each particle adjusts its velocity vector based on the best tracking position P i (t) = [p i1 (t), p i2 (t), . . . , p iD (t)], and based on its neighbor P g (t) that is found as the best update particle of time t, where P g (t) = [p g1 (t), p g2 (t), . . . , p gD (t)]. Then, the essential PSO algorithm can be established in Equations (1) and (2): where ω is the inertia weighting factor, c 1 and c 2 as learning rates. The variables rand 1 and rand 2 are random distribution values [0-1].
The PSO algorithm has the advantage of a rapid convergence in 3D searching space. However, its weakness is it is easily trapped in a local minimum. Consequently, the combination of PSO and EP will cover this issue and furthermore create the new elite algorithm.
EP is one of the four primary evolutionary paradigms and does not include the crossover operator. However, in contrast to the GA, EP employs the mutation and the strategy of individual selection. The EP population members consider each of the species, thus a new offspring is rapidly created [1]. Hence, this study proposed the combination of PSO and EP to create the elite fusion PSO-EP algorithm.
In general, many problems in different scientific disciplines and real-time technical applications can be expressed as a constrained optimization problem and this is presented in Equation (3): where Ω is the constraint set, f (E) is utilized as the benchmark function and E j = [e 1 , e 2 , . . . , e pop ], j = {1, 2, . . . pop} represents the variable vectors.
Initially, an individual population is generated randomly. Each individual is represented by a set of PID parameter gains [k P , k I , k D ] and the fuzzy control gains [k e , k de ]. The membership function of fuzzy µ is verified by using: µ e(t) : U e(t) → −k e(t) , k e(t) , µ de(t) : U de(t) → −k de(t) , k de(t) .
The optimization problem is related to find the minimal of the performance index f (E) with five parameters: K p , K I , K D , K e , K de ∈ E. The progression is conducted as follows: Step 1. Velocity strategy updated. The fitness function value of each individual in the population is computed. The x i (t) value of each particle is compared with the resulting value of x i (t) from the contemporary evaluation. If the next value of x i (t) is superior than the former value, it will be set as the next location. Equations (1) and (2) are used to update the particle velocity.
Step 2. Position strategy updated. The Equation (4) is used to update particle's position and the mutation rate value mr is set. The positions p i (t) and p g (t) are updated using Equations (1) and (2).
Step 3. Mutation EP strategy. Each of x i (t) is mutated and then generated X j (t) = [x j1 (t), x j2 (t), . . . , x jD (t)] using Equation (5): where x jD represents the jth individual error and N(0, σ 2 ) denotes a normal distribution Step 4. Reproduction of EP. Apply roulette wheel selection technique to every particle j.
Step 5. Repeat until the optimized fitness function f (E) is achieved.

Optimal Control Design
To improve the dynamic response and also reduce/eliminate the steady state error of a system, the PID controllers are introduced. Moreover, in order to improve the quality performance of the controlled systems, the fuzzy controller is usually combined to the PID control. In fuzzy controller design, the tracking error e(t) and the differential tracking error de(t) are used as fuzzy inference system inputs, as displayed in Figure 1. The conventional triangular membership functions of e(t) and de(t) are used as inputs, and the output is e Fuzzy (t).
Mathematics 2020, 8, x FOR PEER REVIEW 4 of 16 contemporary evaluation. If the next value of xi(t) is superior than the former value, it will be set as the next location. Equations (1) and (2) are used to update the particle velocity.
Step 2. Position strategy updated. The Equation (4) is used to update particle's position and the mutation rate value mr is set. The positions pi(t) and pg(t) are updated using Equations (1) and (2).
Step 3. Mutation EP strategy. Each of xi(t) is mutated and then generated using Equation (5): where xjD represents the j th individual error and N(0, σ 2 ) denotes a normal distribution Step 4. Reproduction of EP. Apply roulette wheel selection technique to every particle j.
Step 5. Repeat until the optimized fitness function f(E) is achieved.

Optimal Control Design
To improve the dynamic response and also reduce/eliminate the steady state error of a system, the PID controllers are introduced. Moreover, in order to improve the quality performance of the controlled systems, the fuzzy controller is usually combined to the PID control. In fuzzy controller design, the tracking error e(t) and the differential tracking error de(t) are used as fuzzy inference system inputs, as displayed in Figure 1. The conventional triangular membership functions of e(t) and de(t) are used as inputs, and the output is eFuzzy (t). By using the set of linguistic fuzzy rules, its mechanism can produce the appropriate action for fuzzy control. The mechanism comprises one comparable control module and one input/output component: two inputs signal e(t) and de(t) and one output signal eFuzzy(t). To reduce the system error, ke(t) and kde(t) were employed. The fuzzification method: the center of the area is utilized to minimum values of the Fuzzy Mandani type, hence, 7 × 7 fuzzy rules are expressed in Table 1  By using the set of linguistic fuzzy rules, its mechanism can produce the appropriate action for fuzzy control. The mechanism comprises one comparable control module and one input/output component: two inputs signal e(t) and de(t) and one output signal e Fuzzy (t). To reduce the system error, k e(t) and k de(t) were employed. The fuzzification method: the center of the area is utilized to minimum values of the Fuzzy Mandani type, hence, 7 × 7 fuzzy rules are expressed in Table 1. Because this combination yields basic execution factors of the fuzzy control rules, seven partitions were formed: negative large (NL), negative medium (NM), negative small (NS), zero (ZE), positive small (PS), positive medium (PM), and positive large (PL).
The satisfactory performance requirement obtained by using fuzzy-PID parameter controllers has been presented in many empirical studies [16][17][18]29]. The formal structure of the fuzzy-PID controller can be expressed as in Equation (6): The gain coefficients used in the PSO-EP algorithm for discovering the P, I, and D values are evaluated as follows: k P , k I , k D in range [0, 50]. Based on the PID output range, the range of reference input values of e(t) and de(t) of fitness functions f(E) that are generated using Equation (7) as follows: In this study, the fitness function included five components: the rising time, the settling time, the peak of system response, the maximum overshoot, and the minimum fitness function criterion (α 5 , α 6 , α 7 ), as described in Equation (8): where α i represents weighting factors. The output performance was evaluated by using the fitness function. Then, the values of these parameters were recorded. Based on the high output performance of the controller design, the information of response including: rising time, settling time, peak time and overshoot (or undershoot) make the decision about the best performance. Fastest response means that the fastest rising time is also fastest settling time, related to α 1 and α 2 , the next is without overshoot or peak, relative to α 3 and α 4 and the final is the steady state error that is minimum, relative to α 5, α 6 and α 7 . We assigned appropriate the weights [α 1 , α 2 , α 3 , α 4 ] = [30,25,10,5] and the corresponding relative weights [α 5 , α 6 , α 7 ] = [7, 5, 3] were also computed. The error in the output performance is highlighted. Therefore, if the fitness function f(E) is rising then the hybrid PSO-EP algorithm will make the decision that this particle must be ignored in the group of particles. Thus, the particle groups that cover an enormous concealed error can be eliminated, and the system convergence speed can be increased. For evaluating the control design quality, the definite perform index is exploited. The fusion PSO-EP algorithm reinforced using fuzzy-PID controllers was used to assign optimal control gains for a tricopter to achieve a trade-off between flight pilot convenience and system stability, and this algorithm is explained in Figure 2. The system performance index uses the IAE criterion, that is [39,40] as follows: IAE = ∞ 0 e(t) .dt.

Flight Hovering Mode
The UAVs possessing the property of vertical take off and landing (VTOL-UAVs) will have the formation like a helicopter or drone/quadrotor, hexa-copter or octo-copter. Nevertheless, more rotors lead to a more risky operation and higher power consumption. Hence, this paper presents a novel tricopter UAV that may be applied successfully to agriculture activities at any environment. Agricultural Tricopter model, based on the Newton-Euler formulation was developed and obtained from previous studies [22][23][24]31,41,42], as shown in Figure 3. The UAV was described using two lefthand axes-one is the axes of the Earth coordinate system, and the other is a body frame. The positive x axis is aligned with the two front rotors (rotor 1 and 2) and the positive y axis is with the right (rotor 2). However, the positive z axis is directed downward. The roll φ, pitch θ, and yaw ψ triple pilot angles were with the positive x, y, and z axes, respectively. The tilt angle μ was determined by using

Flight Hovering Mode
The UAVs possessing the property of vertical take off and landing (VTOL-UAVs) will have the formation like a helicopter or drone/quadrotor, hexa-copter or octo-copter. Nevertheless, more rotors lead to a more risky operation and higher power consumption. Hence, this paper presents a novel tricopter UAV that may be applied successfully to agriculture activities at any environment. Agricultural Tricopter model, based on the Newton-Euler formulation was developed and obtained from previous studies [22][23][24]31,41,42], as shown in Figure 3. The UAV was described using two left-hand axes-one is the axes of the Earth coordinate system, and the other is a body frame. The positive x axis is aligned with the two front rotors (rotor 1 and 2) and the positive y axis is with the right (rotor 2). However, the positive z axis is directed downward. The roll ϕ, pitch θ, and yaw ψ triple pilot angles were with the positive x, y, and z axes, respectively. The tilt angle µ was determined by using the y-z coordinate axis. The dynamic of the agricultural tricopter has six degrees of freedom. One tilt angle and three rotor speeds are the four inputs to pilot control. The UAVs have commands similar to those of a traditional helicopter: collective, lateral, longitudinal, and yaw or pedal [22][23][24]31,41,42]. The commands are expressed as δcol (collection), δlat (roll control), δlon (pitch control), and δped (yaw control). The roll (φ) control was created because the two front rotors were maneuvered at different speeds. For example, when the speed of rotor 1 is low and the speed of rotor 2 is high, the UAV moves toward the left and vice versa. When the third tail rotor changes its velocity, pitch (θ) control is realized. Yaw (ψ) control is conducted by varying the tilt angle μ. The model notations are described in Table 2. The dynamics of the designed model of a UAV's rigid body were derived using Newton's laws and could be easily translated and rotated in a three-dimensional space. The motion of nonlinear equations is expressed in Equations (9) Rotation equations: The dynamic of the agricultural tricopter has six degrees of freedom. One tilt angle and three rotor speeds are the four inputs to pilot control. The UAVs have commands similar to those of a traditional helicopter: collective, lateral, longitudinal, and yaw or pedal [22][23][24]31,41,42]. The commands are expressed as δcol (collection), δlat (roll control), δlon (pitch control), and δped (yaw control). The roll (ϕ) control was created because the two front rotors were maneuvered at different speeds. For example, when the speed of rotor 1 is low and the speed of rotor 2 is high, the UAV moves toward the left and vice versa. When the third tail rotor changes its velocity, pitch (θ) control is realized. Yaw (ψ) control is conducted by varying the tilt angle µ. The model notations are described in Table 2. The dynamics of the designed model of a UAV's rigid body were derived using Newton's laws and could be easily translated and rotated in a three-dimensional space. The motion of non-linear equations is expressed in Equations (9)-(11): Force equations: Moment equations:

The Agricultural Tricopter Notation
The external moments (u, v, w) The translational velocities (p, q, r) The rotational velocities (angular velocities) (Φ, θ, Ψ) The rotational angles (Roll, Pitch, Yaw) (I xx , I YY , I ZZ ) The rotational inertias u 2 , u 3 , u 4 The channel inputs (Roll, Pitch, Yaw) The hovering mode is the flying mode where the object hangs in the air. This phenomenon occurs when a physical force against the attraction of the gravity is applied; thus the object is stable. The hovering state equilibrium point, named as trimpoint, must be found out because of the tricopter stability issue mentioned previously [41][42][43][44]. The Newton-Raphson method [45] that possesses the vital benefit is rapid quadratic convergence, is the superior selection to calculate the trimpoint of the tricopter non-linearization model. This model has four equations as demonstrated in Equation (12): where b is thrust factor, d is drag factor. Ω 1 , Ω 2 , Ω 3 are the rotor speeds of three rotorcrafts, respectively and µ is the tail tilt angle. Assume that the thrust T i = bΩ i 2 and torque Q i = dΩ i 2 is constant mean that b and d factor are fixed. U T is the total thrust of UAV in the z-axis. U X , U Y , U Z are roll, pitch and yaw moments, respectively. Consequently, the system is asymptotically stable. Therefore, the linear quadratic regulator (LQR) control, which is an appropriate state feedback controller, is exploited to warrant and to keep the system in stable state.

Linear Quadratic Regulator (LQR) Controller Design
The Tricopter system must be asymptotically stable to maintain the hovering mode at all the operation time. Hence, the continuous system exists a state feedback control with the input is u(t) = −k.x(t). The feedback gain matrix K is solved by using the Riccati equation A T P + PA + Q − PBR −1 B T P = 0. If k = −R −1 B T P, then the close loop system is asymptotically stable as shown in Figure 4.

The Agricultural Tricopter Notation FX, FY, FZ
The External Forces L, M, N The external moments (u, v, w) The translational velocities (p, q, r) The rotational velocities (angular velocities) The rotational angles (Roll, Pitch, Yaw) (Ixx, IYY, IZZ) The rotational inertias u2, u3, u4 The channel inputs (Roll, Pitch, Yaw) The hovering mode is the flying mode where the object hangs in the air. This phenomenon occurs when a physical force against the attraction of the gravity is applied; thus the object is stable. The hovering state equilibrium point, named as trimpoint, must be found out because of the tricopter stability issue mentioned previously [41][42][43][44]. The Newton-Raphson method [45] that possesses the vital benefit is rapid quadratic convergence, is the superior selection to calculate the trimpoint of the tricopter non-linearization model. This model has four equations as demonstrated in Equation (12): where b is thrust factor, d is drag factor. 1 2 3 , , Ω Ω Ω are the rotor speeds of three rotorcrafts, respectively and μ is the tail tilt angle. Assume that the thrust is constant mean that b and d factor are fixed. UT is the total thrust of UAV in the z-axis. UX, UY, UZ are roll, pitch and yaw moments, respectively. Consequently, the system is asymptotically stable. Therefore, the linear quadratic regulator (LQR) control, which is an appropriate state feedback controller, is exploited to warrant and to keep the system in stable state.

Linear Quadratic Regulator (LQR) Controller Design
The Tricopter system must be asymptotically stable to maintain the hovering mode at all the operation time. Hence, the continuous system exists a state feedback control with the input is u(t) = −k.x(t). The feedback gain matrix K is solved by using the Riccati equation

Numerical Simulation Results
The object simulation parameters that are listed in Table 3, were derived from previous studies [21,22,29]. However, the value of the tilt angle µ should be negative in the initial condition to maintain system stability. The simulation pilot controls are demonstrated in each channel RPY (roll, pitch, and yaw). The speed of each rotorcraft that was determined using the Newton-Raphson method is as follows: Ω1 = 8418 rpm, Ω2 = 8418 rpm, Ω3 = 8676 rpm. Moreover, the tilt angle µ is −0.0210 rad (−2.407°).

Numerical Simulation Results
The object simulation parameters that are listed in Table 3, were derived from previous studies [21,22,29]. However, the value of the tilt angle µ should be negative in the initial condition to maintain system stability. The simulation pilot controls are demonstrated in each channel RPY (roll, pitch, and yaw). The speed of each rotorcraft that was determined using the Newton-Raphson method is as follows: Ω 1 = 8418 rpm, Ω 2 = 8418 rpm, Ω 3 = 8676 rpm. Moreover, the tilt angle µ is −0.0210 rad (−2.407 • ).
According to the pilot angle control of a UAV in the roll ( Figure 5), Pitch (Figure 6), and Yaw ( Figure 7) axes of each channel, we exhibited a useful system model as follows: (a) The fitness function of the proposed controller design. The fitness function is smaller (or the inverse fitness is larger) mean that the optimal convergence is better.

Discussion
The comparison between the results obtained using five controller strategies integrated with the fusion EP-PSO algorithm is presented in Table 4. The five controller strategies that were implemented to pilot the tricopter attitude channels roll, pitch and yaw are: PSO-optimized PID gains, fusion PSO-EP optimized PID gains, fuzzy-PID controller, PSO-optimized fuzzy-PID gains, and the proposed methodology: fusion PSO-EP algorithm reinforced fuzzy-PID controller. Firstly, the PSO algorithm performs well when applied to find out the three optimal P, I and D control gains. The responses are just after 1 second. Secondly, the fusion EP-PSO is exploited to tune PID control, the responses of this method are faster than the original PSO. Next, the fuzzy-PID controller is designed by the Ziegler-Nichols method as mentioned above. This design display the advantage of the fuzzy-PID control method in faster response and stable. Later, The fuzzy-PID control with five parameters: e, de, P, I and D is successfully optimized by using the fundamental PSO algorithm. This means that the PSO algorithm works well in 5 dimensions. The system responses are faster than 5%. Finally, the proposed fusion EP-PSO algorithm is implemented well to show the fastest responses just after 0.2 s. However, the responses perform the control trade-off in jagged/sawtooth phenomenon to achieve the faster target but neither error nor overshoot. Moreover, the steady-state error is almost zero mean that the errors, between the input command and output response, of the proposed controllers design to tricopter models are tiny. Generally, all of the

Discussion
The comparison between the results obtained using five controller strategies integrated with the fusion EP-PSO algorithm is presented in Table 4. The five controller strategies that were implemented to pilot the tricopter attitude channels roll, pitch and yaw are: PSO-optimized PID gains, fusion PSO-EP optimized PID gains, fuzzy-PID controller, PSO-optimized fuzzy-PID gains, and the proposed methodology: fusion PSO-EP algorithm reinforced fuzzy-PID controller. Firstly, the PSO algorithm performs well when applied to find out the three optimal P, I and D control gains. The responses are just after 1 second. Secondly, the fusion EP-PSO is exploited to tune PID control, the responses of this method are faster than the original PSO. Next, the fuzzy-PID controller is designed by the Ziegler-Nichols method as mentioned above. This design display the advantage of the fuzzy-PID control method in faster response and stable. Later, The fuzzy-PID control with five parameters: e, de, P, I and D is successfully optimized by using the fundamental PSO algorithm. This means that the PSO algorithm works well in 5 dimensions. The system responses are faster than 5%. Finally, the proposed fusion EP-PSO algorithm is implemented well to show the fastest responses just after 0.2 s. However, the responses perform the control trade-off in jagged/sawtooth phenomenon to achieve the faster target but neither error nor overshoot. Moreover, the steady-state error is almost zero mean that the errors, between the input command and output response, of the proposed controllers design to tricopter models are tiny. Generally, all of the controller designs are achieved the high-quality performance without overshoot, high precision and high reliability.

Conclusions
Five controller designs were applied to a tricopter pilot system in this research work. The proposed fusion PSO-EP algorithm is reinforced to determining the parameters of the fuzzy-PID controllers, it then fruitfully implemented the tricopter attitude models in each pilot channel: roll, pitch, and yaw. The proposed controller design provides the optimal pilot performance in only after 0.2 s when compared with the other controller designs. The accomplished performance reveals that the proposed controller provides the fastest response, no overshoot, high reliability, and high stability. The efficient controller schemes for the non-linear models of an agricultural tricopter are successfully achieved. In further work, the disturbance attack to the system models would be also taken into consideration and the most imperative need is implementing the proposed algorithms to real-time systems.