<?xml version="1.0" encoding="UTF-8"?>
<article xmlns:xlink="http://www.w3.org/1999/xlink">
<front><article-meta>
<title-group><article-title>Multi-Objective Optimization of a Fuzzy Tracking Controller Using Differential Evolution for Nonlinear Three-Dimensional Tracking</article-title></title-group>
<contrib-group>
<contrib contrib-type="author"><name><surname>Najafi</surname><given-names>Seyed Esmaeil</given-names></name><xref ref-type="aff" rid="aff1"/></contrib>
<contrib contrib-type="author"><name><surname>Poorahangryan</surname><given-names>Fereshteh</given-names></name><xref ref-type="aff" rid="aff2"/></contrib>
<contrib contrib-type="author"><name><surname>Jamali</surname><given-names>Ali</given-names></name><xref ref-type="aff" rid="aff3"/></contrib>
<aff id="aff1">Department of Industrial Engineering, S.R.C., Islamic Azad University, Tehran, Iran.</aff>
<aff id="aff2">Department of Electrical Engineering, Ayandegan University, Tonekabon, Iran.</aff>
<aff id="aff3">Department of Engineering, RMIT University, Melbourne, Australia.</aff>
</contrib-group>
<abstract><p>The issue of three-dimensional tracking in a system involving nonlinear dynamics, couplings, uncertainties, and varying working conditions has become a challenge in current control engineering practice. It is of vital importance that in such a case, the precision of tracking is maximized while the control effort and energy consumed are minimized. In this paper, the optimized Fuzzy Logic Controller (FLC) has been designed for the three-dimensional tracking of a missile (satellite) through application of the Differential Evolution (DE) algorithm. The fuzzy membership function parameters are considered as design parameters, and the minimization of both the tracking error and energy is performed simultaneously as conflicting objectives. The parameters of the membership functions have been optimized through the use of the DE algorithm, and the performance of the controller has been enhanced. The results obtained through simulation prove that the fuzzy control based on DE yields an ideal compromise between tracking performance and energy utilization as compared to Proportional, Integral, and Derivative (PID) and Genetic Algorithm (GA)-based controllers in terms of overshoot and settling time. The proposed method serves as an efficient technique for the optimal design of fuzzy controllers for tracking problems in three-dimensional space.</p></abstract>
<kwd-group kwd-group-type="author"><kwd>Fuzzy logic control</kwd><kwd>Fuzzy controller</kwd><kwd>Three-dimensional tracking</kwd><kwd>Missile tracking</kwd><kwd>Differential</kwd></kwd-group>
</article-meta></front>
<body>
<sec id="sec-1">
<p>evolution, Multi-objective optimization, Genetic algorithm.</p>
</sec>
<sec id="sec-2"><label>1</label><title>Introduction</title>
<p>The control of trajectories is a central task in contemporary control engineering that is applicable in mobile robots, automatic cars, Uncrewed Air Vehicles (UAVs), aerospace systems, marine vessels, and intelligent transport systems. The key purpose of trajectory tracking control is for the actual state of the dynamical system to accurately track the desired trajectory with minimal tracking error while assuring stability and satisfactory performance. However, in real-life systems, the controller performance is not solely judged by the ability to perform accurate tracking. Control efforts, energy usage, transient time, disturbance rejection, and actuator constraints are other important factors to take into account. They become even more challenging when designing a control system for a nonlinear dynamical system under uncertainty, disturbances, and changing operating conditions. Thus, the creation of intelligent controllers for accurate tracking with low control efforts is a current research direction.</p>
<p>Traditional model-based control approaches have found extensive application in trajectory tracking due to their sound theoretical background. However, the performance of these approaches is largely dependent on the accuracy of the mathematical model and proper tuning of the controller parameters. For complex and highly nonlinear systems, developing an accurate model is not always feasible, and errors in the modeling process and unmodeled dynamics could greatly impact the performance of the controller. In addition, a controller that is developed for a particular operating point might fail if changes occur in the dynamics of the system or operating conditions.</p>
<p>Fuzzy Logic Control (FLC) is one of the most mature intelligent control techniques designed for nonlinear and uncertain processes. The theoretical basis of fuzzy set theory was developed by Zadeh <xref ref-type="bibr" rid="ref-1">[1]</xref>, whereas the applicability of linguistic IF-THEN rules in control was shown by Mamdani and Assilian <xref ref-type="bibr" rid="ref-2">[2]</xref>. Fuzzy control systems can model relationships between control variables in terms of linguistic rules and take advantage of human expertise in designing control. Such properties make FLC a promising approach to nonlinear systems, where it is hard to build an accurate mathematical model. However, the effectiveness of fuzzy control greatly depends on the choice of parameters, such as membership functions, scaling factors, and rule-base parameters. Classical manual tuning and trial-and-error methods become impractical with the increase in the number of these parameters.</p>
<p>This is where evolutionary optimization techniques come into play, as they can effectively search through nonlinear parameter spaces without the need to know gradients. Among these techniques, Differential Evolution (DE), proposed by Storn and Price <xref ref-type="bibr" rid="ref-3">[3]</xref>, has been widely used due to its simplicity and ability to efficiently search through continuous problem domains having a small number of control parameters. The usage and development of DE in various optimization tasks have been extensively studied <xref ref-type="bibr" rid="ref-4">[4]</xref>, <xref ref-type="bibr" rid="ref-5">[5]</xref>, <xref ref-type="bibr" rid="ref-6">[6]</xref>. Thus, DE seems to be applicable to optimizing continuous fuzzy membership function parameters.</p>
<p>Several works were dedicated to the application of fuzzy control systems in combination with evolutionary optimization algorithms. Specifically, Pishkenari et al. <xref ref-type="bibr" rid="ref-7">[7]</xref> used DE and Genetic Algorithm (GA) to optimize a fuzzy controller for trajectory tracking of the CEDRA mobile robot. They showed the superiority of these optimized membership functions in terms of performance compared to manually designed ones. Later, Pishkenari et al. <xref ref-type="bibr" rid="ref-8">[8]</xref> utilized DE for optimal design of a fuzzy trajectory-tracking controller and compared its performance with that of GA.</p>
<p>The ability of DE to perform optimization in fuzzy controller design has also been exhibited in nonlinear control systems. Chen et al. <xref ref-type="bibr" rid="ref-9">[9]</xref> introduced an adaptive mutation-based DE technique for TSK-type fuzzy controller optimization where membership function parameters and consequent weight parameters are optimized together. This technique was tested on several nonlinear control problems, and it showed superior search ability and lower control errors. Ochoa et al. <xref ref-type="bibr" rid="ref-10">[10]</xref> introduced a parameter adaptation DE approach to optimize the Type-1 and Interval Type-2 fuzzy controllers. Another work by Ochoa et al. <xref ref-type="bibr" rid="ref-11">[11]</xref> utilized the Interval Type-3 fuzzy system to dynamically adapt the crossover parameter of DE in order to optimize the fuzzy controller of a direct current motor.</p>
<p>DE has also been used in solving multi-objective fuzzy control problems. For instance, Marinaki et al. <xref ref-type="bibr" rid="ref-12">[12]</xref> proposed Multi-Objective Differential Evolution (MODE) to optimize a fuzzy control system for suppressing vibrations in smart structures. From their experiments, DE was shown to be effective in dealing with control design problems that have many and possibly conflicting objectives. The problem of multiple objectives is significant in trajectory tracking since minimizing tracking errors might demand more control activities, hence more energy expenditure. In another study, Subraj et al. <xref ref-type="bibr" rid="ref-14">[13]</xref> proposed a fuzzy-controlled diversity-based DE for noisy multi-objective optimization problems, showing the effectiveness of fuzzy approaches in enhancing DE performance.</p>
<p>Other evolutionary algorithms have been used for optimization-based fuzzy control in trajectory-tracking systems. Tolossa et al. <xref ref-type="bibr" rid="ref-15">[14]</xref> carried out research on FLC for trajectory tracking of a two-wheeled differential-drive mobile robot through PSO for optimizing parameters of the membership function and GA for automated rule base reduction. The results obtained through their simulation and experimentation proved that through optimization, more accurate tracking can be achieved while at the same time reducing the control efforts. Several other recent studies have been able to show how DE-based fuzzy control can be applied in complicated engineering systems. For instance, Ouda and Alavifard <xref ref-type="bibr" rid="ref-16">[15]</xref> developed an adaptive fuzzy-DE method for optimizing damping and inertia coefficients of virtual synchronous generators. Wang et al. <xref ref-type="bibr" rid="ref-17">[16]</xref> designed a robust fuzzy control method for trajectory tracking of autonomous vehicles using multi-objective DE. Chen et al. <xref ref-type="bibr" rid="ref-18">[17]</xref> presented a modified constrained DE algorithm to optimize an interval Type-2 fuzzy Proportional, Integral, And Derivative (PID) controller for hydraulic actuators of humanoid robots.</p>
<p>Even with all these advances, there still exist some research areas that need to be looked into. Quite a number of studies have concentrated their attention on mobile robots and planar trajectory tracking; however, complex three-dimensional tracking systems have not been given much attention. Besides, quite a lot of work has concentrated on minimizing tracking error without giving much thought to control effort or energy expenditure, which should also be taken into account. Even though DE has proven to be effective at solving optimization problems and tuning fuzzy controllers, its use in the field of fuzzy controller optimization for missile three-dimensional tracking has not yet received adequate attention.</p>
<p>Therefore, the current research aims to design an optimal fuzzy controller for three-dimensional missile tracking through the DE algorithm. Fuzzy membership function parameters are considered as design variables, and they are optimized automatically using the DE algorithm. The objective function evaluates both tracking error and control effort so that the tradeoff is made between the two aspects. Performance of the developed controller is analyzed based on tracking error, control effort, and dynamic response, and compared with other optimization and control methods, such as GA.</p>
<p>The rest of this paper outlines the systematic design and analysis of the proposed fuzzy control scheme. First, the mathematical formulation and state-space model of the nonlinear tracking system are derived. The performance of a traditional PID controller is analyzed as a baseline, and afterwards a Mamdani fuzzy controller with optimized control rules using a GA is designed. Finally, a MODE method is used to find the optimum output torque membership functions of the fuzzy control system by taking into account both the tracking error and the control effort simultaneously. These Pareto optimal results are studied in terms of their dynamic response and are compared to the conventional PID and GA-based fuzzy controllers.</p>
</sec>
<sec id="sec-3"><label>2</label><title>Dynamic Modeling of the Tracking System</title>
<sec id="sec-4"><label>2.1</label><title>System Structure and Reference Frame</title>
<p>The system that is taken into account in this research is a satellite with two degrees of rotational freedom, the motion of which is determined by two rotations along two different axes. The configuration structure of the satellite and the reference frame used are shown in <xref ref-type="fig" rid="fig1">Figs. 1</xref> and 2. Two generalized coordinates of the system are referred to as (θ) and (α), while two control torques are denoted as (M1) and (M2). The center of mass of the satellite is situated at the origin of the reference frame.</p>
<fig id="fig1"><label>Figure 1</label><caption><p>Schematic of the satellite and its rotation directions.</p></caption><graphic content-type="pdf-region" xlink:href="pdf-region:3:66,631.4,537.1,780.3"/><p content-type="note">This figure is shown in the PDF version of the article.</p></fig>
<fig id="fig2"><label>Figure 2</label><caption><p>Schematic of the satellite and the body-fixed coordinate system.</p></caption><graphic content-type="pdf-region" xlink:href="pdf-region:3:66,461.2,537.1,616.9"/><p content-type="note">This figure is shown in the PDF version of the article.</p></fig>
</sec>
<sec id="sec-5"><label>2.2</label><title>Derivation of the Equations of Motion</title>
<p>In order to obtain the equations of motion, the Lagrange equations are used. They have the following general form:</p>
<p>d ∂L ∂L</p>
<p>( ) − = δQ, (1)</p>
<p>dt ∂q̇ ∂q where L = T - U, T is the kinetic energy, U is the potential energy, Q is the differential work, and q is the generalized coordinate. Due to the immobility of the center of mass, the variations of the potential energy are equal to zero (U = 0), and hence L = T.</p>
</sec>
<sec id="sec-6"><label>2.3</label><title>Calculation of Kinetic Energy</title>
<p>For calculating the kinetic energy, an equivalent rod structure is assumed as in <xref ref-type="fig" rid="fig2">Fig. 2</xref>. Considering an observer on the first rod, there is only the angular velocity θ̇, and the rotational movement about the e₃ axis does not affect it. Since the inertia matrix of the first rod is expressed as I₁ = diag(I₁ₓ, I₁ᵧ, I₁₂), its kinetic energy is as follows:</p>
<p>T1 = Iyy1θ̇ . (2)</p>
<p>For the second rod, the observer observes the angular velocity α̇ and also θ̇. Since the inertia matrix of the second rod is expressed as I₂=diag(I₂ₓ, I₂ᵧ, I₂₂), its kinetic energy can be calculated from the following equation:</p>
<p>1 2 2 2 2 T2 = [Ixx2(θ̇ cos α) + Iyy2(θ̇ sin α) + Izz2α̇ + (Iyx2 + Ixy2)θ̇ cos α sinα</p>
<p>2 (3)</p>
<p>+((Izx2 + Ixz2)θ̇ cos α + (Izy2 + Iyz2)θ̇ sinα)α̇ ].</p>
<p>Hence, the total kinetic energy of the system will be calculated by adding the kinetic energies of the two rods:</p>
<p>T = T1 + T2. (4)</p>
<p>The change in the potential energy of the system is equal to zero since the center of mass of the system is fixed at the origin of the coordinate system, and therefore there is no change in altitude. The Lagrangian of the system is equal to the kinetic energy of the system. Upon substitution of (T) and (U) in the Lagrange equation, the governing dynamic equations of the system can be found.</p>
</sec>
<sec id="sec-7"><label>2.4</label><title>Derivation of Equations of Motion by Using the Lagrange Equations</title>
<p>By substitution of the generalized coordinates q = θ and q = α in Eq. (1), the equations of motion for each degree of freedom will be found.</p>
<sec id="sec-8"><label>2.4.1</label><title>First case: Substitution of q = θ</title>
<p>Upon substitution of q = θ in Eq. (1):</p>
<p>d (∂L) = 1 d [2Ixx θ̇ cos2 α + 2Iyy θ̇ sin2 α + 2(Iyx + Ixy )θ̇ cos α sinα,</p>
<p>dt ∂q̇ 2 dt 2 2 2 2 (4)</p>
<p>((Izx + Ixz ) cos α + (Izy + Iyz ) sin α)α̇ + 2Iyy θ̇].</p>
<p>2 2 2 2 1</p>
<p>= 0. (5)</p>
<p>δQ = M1θ, (6)</p>
<p>where</p>
<p>Ixx (θ̈ cos2 α − 2θ̇α̇ cosα sin α) + Iyy (θ̈ sin2 α + 2θ̇α̇ cos α sin α), +(Iyx + Ixy )(θ̈ cos α sinα + θ̇(α̇ cos2 α − α̇ sin2 α))</p>
<p>1 1 (7)</p>
<p>+( (Izx2 + Ixz2)α̇ cosα + (Izy2 + Iyz2)α̇ sinα)α̇</p>
<p>+(1 (Izx + Ixz ) cos α + 1 (Izy + Iyz ) sin α)α̈ + Iyy θ̈ = M1.</p>
<p>2 2 2 2 2 2 1</p>
</sec>
<sec id="sec-9"><label>2.4.2</label><title>Second case: Substituting q=α</title>
<p>By substituting q = α into Eq. (1), we have:</p>
<p>d ∂L 1 d</p>
<p>( ) = [2Izz2α̇ + (Izx2 + Ixz2)θ̇ cosα + (Izy2 + Iyz2)θ̇ sinα]. (8)</p>
<p>dt ∂q̇ 2 dt ∂L 1 2 2 2 2 2 = [−2Ixx θ̇ cosα sinα + 2Iyy θ̇ cos α sinα + (Iyx + Ixy )θ̇ (cos α − sin α)</p>
<p>∂q 2 2 2 2 2 (9)</p>
<p>+ ((−(Izx + Ixz )θ̇ sin α) + (Izy + Iyz )θ̇ cos α) α̇ ].</p>
<p>2 2 2 2</p>
<p>δQ = M2, (10)</p>
<p>where</p>
<p>Izz α̈ + (Izx + Ixz )(θ̈ cos α − θ̇α̇ sinα), 2 2 2 2 1 1 2</p>
<p>+ (Izy2 + Iyz2)(θ̈ sinα + θ̇α̇ sin α) + ((−Ixx2θ̇ cos α sinα) (11)</p>
<p>2 1 2 2 2 +Iyy θ̇ cosα sinα) + (Iyx + Ixy )θ̇ (cos α − sin α) 2 2 2 2 + 1 ((−(Izx + Ixz )θ̇ sinα) + 1 (Izy + Iyz )θ̇ cos α)α̇ = M2.</p>
<p>2 2 2 2 2 2</p>
</sec>
</sec>
<sec id="sec-10"><label>2.5</label><title>Conversion of Second Order Equations to State Space Equations</title>
<p>To carry out numerical solution and optimize the system, the second-order Eqs. (4) and (6) need to be expressed as first-order equations through:</p>
<p>θ = y1 , θ̇ = y2 , α = y3,α̇ = y4. (12)</p>
<p>The state-space equations are obtained as follows:</p>
<p>y1 = θ ⇒ ẏ1 = θ̇. (13)</p>
<p>y2 = θ̇ ⇒ ẏ 2 = θ̈. (14)</p>
<p>y3 = α ⇒ ẏ3 = α̇ . (15)</p>
<p>y4 = α̇ ⇒ ẏ4 = α̈ . (16)</p>
<p>The set of Eqs. (6)-(9) is the last state space representation of the system, and the same equations will be used for the remaining analysis in this study.</p>
<p>(7854.81 cos2 α + 6461.67 sin2 α + 0.915632)θ̈ − 2786.28θ̇α̇ cosα sin α = M1</p>
<p>(17)</p>
<p>M1 + 2786.28θ̇α̇ cos α sinα ⇒ θ̈ = .</p>
<p>7854.81 cos2 α + 6461.67 sin2 α + 0.915632 1 2 2 6379.30α̈ + ((−7854.81θ̇ cos α sin α) + 6461.67θ̇ cos α sin α) = M2</p>
<p>(18)</p>
<p>M2+696.57θ̇ 2 cos α sin α ⇒ α̈ = .</p>
<p>6379.30 Now, by combining Eqs. (13)-(16), the following equations are obtained:</p>
<p>ẏ1 = y2. (19)</p>
<p>M1 + 2786.28y2y4 cos y3 sin y3</p>
<p>ẏ2 = 2 2 . (20)</p>
<p>7854.81 cos y3 + 6461.67 sin y3 + 0.915632</p>
<p>ẏ3 = y4. (21) ẏ4 = M2+696.57y22 cos y3 siny3. (22)</p>
<p>6379.30</p>
</sec>
</sec>
<sec id="sec-11"><label>3</label><title>Control of the System Using PID Controller</title>
<p>At first, the PID controller is applied to the set of equations for the system under consideration in this thesis in order to see the response of the system when an ordinary and suitable controller controls it.</p>
<p>A PID controller is made up of a proportional gain, an integrator, and a differentiator having constant coefficients. With respect to the coefficients of the PID in <xref ref-type="table" rid="tab1">Table 1</xref>, the response of the system for the unit step input will be shown in two modes.</p>
<table-wrap id="tab1"><label>Table 1</label><caption><p>PID coefficients.</p></caption><graphic content-type="pdf-region" xlink:href="pdf-region:6:66,706.3,537.1,757.6"/><p content-type="note">This table is shown in the PDF version of the article.</p></table-wrap>
<p>In <xref ref-type="fig" rid="fig3">Fig. 3</xref>, the horizontal axis shows the time while the vertical axis shows the response (error) for the first mode. The first mode includes the values of θ and θ̇.</p>
<fig id="fig3"><label>Figure 3</label><caption><p>Response of the system to unit step input for the first mode.</p></caption><graphic content-type="pdf-region" xlink:href="pdf-region:6:66,500,537.1,660.7"/><p content-type="note">This figure is shown in the PDF version of the article.</p></fig>
<p>Again, in <xref ref-type="fig" rid="fig4">Fig. 4</xref>, the horizontal axis shows the time while the vertical axis shows the response for the second mode. The second mode includes the values of α and α̇.</p>
<fig id="fig4"><label>Figure 4</label><caption><p>System response to unit step input for the second mode.</p></caption><graphic content-type="pdf-region" xlink:href="pdf-region:6:66,288.3,537.1,450.2"/><p content-type="note">This figure is shown in the PDF version of the article.</p></fig>
<p>As seen above, the response of the first variable has about 50% overshoot and a settling time of 10 seconds. In the case of the second variable, the response has 60% overshoot and a settling time of 15 seconds. This kind of response is not at all acceptable for a military system where the requirement is to track a target very quickly. Thus, this particular controller, despite all its merits, cannot be used for this system. With the analysis of the response of the PID controller being completed now, the next control strategy to discuss is the fuzzy controller. In a fuzzy controller, some rules as well as membership functions are required, which are normally provided by people who know about the system very well.</p>
<p>However, there is another way of doing it. The system rules can be obtained through optimization, where an optimal chromosome will be generated instead of having to depend on an expert. In fact, the optimization technique, after some time and through many generations, will result in the generation of a chromosome that is able to perform at the same level as the expert designer of the fuzzy controller. The only problem with this approach is that the membership functions of the system have to be initialized. This means that the membership functions do not have any other way other than depending on the initial knowledge of a person familiar with the system.</p>
</sec>
<sec id="sec-12"><label>4</label><title>Optimization using Genetic Algorithm</title>
<p>Before optimizing the torque, it is essential to distinguish the membership functions related to that torque. Torque M₁ has the two variables y1 = θ and y2 = θ̇. In the case of optimizing torque M₁, only the membership functions of those two parameters will be used, as shown in <xref ref-type="fig" rid="fig5">Fig. 5</xref>.</p>
<fig id="fig5"><label>Figure 5</label><caption><p>Relation between torque M₁ and the input membership functions.</p></caption><graphic content-type="pdf-region" xlink:href="pdf-region:7:66,550.2,537.1,698.4"/><p content-type="note">This figure is shown in the PDF version of the article.</p></fig>
<p>The same thing applies to torque M₂, which uses the variable y₃ = α y4 = α̇. This can be seen from <xref ref-type="fig" rid="fig6">Fig. 6</xref>.</p>
<fig id="fig6"><label>Figure 6</label><caption><p>Relation between torque M₂ and the input membership functions.</p></caption><graphic content-type="pdf-region" xlink:href="pdf-region:7:66,359.1,537.1,514.5"/><p content-type="note">This figure is shown in the PDF version of the article.</p></fig>
<p>In order to choose the right response, one should use an algorithm that solves one problem only, and in which the cost function is designed to minimize the error due to the step response of each of the states α and θ. In the cost function, the error of each state variable is calculated compared to its target value by using numerical integration using the trapezoidal rule. In addition to the complexity of this problem and the lengthy processing time needed to solve it, a small initial population size was selected. However, a mutation operator is included in the algorithm to optimize the current response.</p>
<sec id="sec-13"><label>4.1</label><title>Fuzzy Control System Design with Optimized Rules</title>
<p>After analyzing the performance of the PID control system, this part will consider the design of a fuzzy control system. The fuzzy system used is of Mamdani type, having five triangular membership functions for each input and two output torques, namely M₁ and M₂. The complete fuzzy control system has four inputs such that y1 = θ, y = θ̇, y = α and y₄ = α̇.</p>
<sec id="sec-14"><label>4.1.1</label><title>Design of membership functions</title>
<p>These membership functions have been defined with the help of previous knowledge of the system, and the objective is to make sure that the entire input range is covered (condition of completeness). Due to the same pattern of membership functions for each of the input variables, <xref ref-type="table" rid="tab2">Table 2</xref> shows the membership functions along with GA parameters.</p>
<table-wrap id="tab2"><label>Table 2</label><caption><p>Membership functions (for variables y1,y2,y3,y4) and GA parameters.</p></caption><graphic content-type="pdf-region" xlink:href="pdf-region:8:66,587.3,537.1,757.6"/><p content-type="note">This table is shown in the PDF version of the article.</p></table-wrap>
<p>Description of <xref ref-type="table" rid="tab2">Table 2</xref>: The first column depicts the linguistic labels used. The mathematical formulation for each of the membership functions in terms of the variable x is given in the second column. The third and fourth columns show the GA tuning parameter names and their respective values. These functions are identical for all four input variables y₁ to y₄; the only difference is the variable name used. The output membership functions of the torque M₁ and M₂ are taken to be symmetric from –10 to +10 kN·m.</p>
</sec>
</sec>
<sec id="sec-15"><label>4.2</label><title>Comparison of Controller Performance</title>
<p>In order to perform a quantitative comparison of the optimized fuzzy controller, a comparison with the performance of the classical PID controller is conducted. The coefficients of the PID controller are shown in <xref ref-type="table" rid="tab3">Table 3</xref> below:</p>
<table-wrap id="tab3"><label>Table 3</label><caption><p>PID controller coefficients.</p></caption><graphic content-type="pdf-region" xlink:href="pdf-region:8:66,364.8,537.1,416"/><p content-type="note">This table is shown in the PDF version of the article.</p></table-wrap>
<p>Membership functions for the four input variables (y1) to (y4) are depicted in <xref ref-type="fig" rid="fig7">Fig. 7</xref>. The normalized input space is considered over the range ([-1,1]) and consists of five fuzzy sets, namely, BN (Big-Negative), LN (Low-Negative), ZE (Zero), LP (Low-Positive), and BP (Positive). Triangular membership functions with proper overlap have been used to ensure a smooth transition from one fuzzy set to another. The same membership function structure has been considered for all four inputs, with the only difference being in their</p>
<fig id="fig7"><label>Figure 7</label><caption><p>Membership functions for the four input variables (y_1), (y_2), (y_3), and (y_4).</p></caption><graphic content-type="pdf-region" xlink:href="pdf-region:8:65.7,69.4,537.1,283.1"/><p content-type="note">This figure is shown in the PDF version of the article.</p></fig>
</sec>
<sec id="sec-16"><label>4.3</label><title>Optimization of Fuzzy Rules by Genetic Algorithm</title>
<p>Now, utilizing the values of the membership functions defined in <xref ref-type="table" rid="tab4">Table 4</xref>, the fuzzy rules are optimized by means of a GA. The parameters of the GA are given in <xref ref-type="table" rid="tab4">Table 4</xref> below.</p>
<table-wrap id="tab4"><label>Table 4</label><caption><p>GA optimization parameters.</p></caption><graphic content-type="pdf-region" xlink:href="pdf-region:9:66,624.6,537.1,698.3"/><p content-type="note">This table is shown in the PDF version of the article.</p></table-wrap>
<p>In the course of optimization, the objective function is designed as the minimization of the Integrated Squared Error (ISE) simultaneously for both states. Being a stochastic method, GA optimization was performed several times with different initial populations of solutions to choose the one providing the minimal objective function value. Dynamic characteristics of the suggested GA-optimized fuzzy controller were calculated and compared with those of the standard PID controller. Step responses of the GA-optimized fuzzy control system for the first and the second operating mode are shown in <xref ref-type="fig" rid="fig8">Fig. 8</xref>. As can be seen, the proposed controller demonstrates much better dynamic characteristics with much faster transient response, less overshoot, and shorter settling time in both operating modes.</p>
<fig id="fig8"><label>Figure 8</label><caption><p>Step response of the system for the first and second mode (θ and α).</p></caption><graphic content-type="pdf-region" xlink:href="pdf-region:9:66,330.5,537.1,496.5"/><p content-type="note">This figure is shown in the PDF version of the article.</p></fig>
</sec>
<sec id="sec-17"><label>4.4</label><title>Comparison of PID and Fuzzy Controllers</title>
<p><xref ref-type="fig" rid="fig9">Fig. 9</xref> compares the transient responses of the PID and GA-based fuzzy controllers for the second operating mode (α). It is seen that the classical PID controller has a relatively high overshoot, equal to about 60%, as well as noticeable oscillations before attaining the steady-state value. The required settling time is 15 s. However, the new fuzzy controller designed using the GA greatly improves the dynamic response of the closed-loop control system. Specifically, the maximum overshoot is minimized to less than 5%, whereas the settling time is less than 5s. Also, the oscillations, which occur with the classical PID controller, are almost entirely suppressed. In other words, the proposed fuzzy controller provides faster and smoother convergence to the reference value. It should be noted that the overshoot reduction provided by the new controller is 55% points, and the settling time reduction is greater than 60% compared to the PID controller.</p>
<p>The same improvement can be seen from the first operating mode (θ) as shown in <xref ref-type="fig" rid="fig10">Fig. 10</xref> below. The optimized fuzzy controller has a much smaller overshoot and converges to the steady state much faster than the PID controller while still having good tracking performance. This indicates that the proposed control scheme is able to achieve better transient behavior in both operating modes.</p>
<p>From the comparative results shown in <xref ref-type="fig" rid="fig9">Figs. 9</xref> and 10 above, it is clear that the proposed GA-based fuzzy controller can effectively improve the dynamic behavior of the tracking system. The optimized fuzzy controller has improved settling time, decreased overshoot, and improved damping performance, without needing the exact mathematical model of the plant. Thus, it can be concluded that the proposed GA-based fuzzy control scheme is an effective solution for the satellite attitude tracking problem.</p>
<fig id="fig9"><label>Figure 9</label><caption><p>Fuzzy versus PID controller for first mode (θ).</p></caption><graphic content-type="pdf-region" xlink:href="pdf-region:10:66,556.5,537.1,698.3"/><p content-type="note">This figure is shown in the PDF version of the article.</p></fig>
<fig id="fig10"><label>Figure 10</label><caption><p>Fuzzy versus PID controller for second mode (α).</p></caption><graphic content-type="pdf-region" xlink:href="pdf-region:10:66,409.3,537.1,542"/><p content-type="note">This figure is shown in the PDF version of the article.</p></fig>
</sec>
</sec>
<sec id="sec-18"><label>5</label><title>Multi-objective Optimization of Output Torque Membership Functions Using DE</title>
<p>In this work, the input membership functions of the fuzzy controller are supposed to have been optimally designed, which is why only the output membership functions corresponding to the control torques are optimized employing the DE technique.</p>
<p>In this optimization problem, both of the following objective functions are minimized: the tracking error and the control effort (torque). After the completion of the evolutionary process, the solutions that are not dominated in the last generation are considered Pareto optimal solutions. The parameters of MODE algorithm used in this work are listed in <xref ref-type="table" rid="tab5">Table 5</xref>.</p>
<table-wrap id="tab5"><label>Table 5</label><caption><p>Parameters of the DE algorithm.</p></caption><graphic content-type="pdf-region" xlink:href="pdf-region:10:66,108.9,537.1,203"/><p content-type="note">This table is shown in the PDF version of the article.</p></table-wrap>
<p>Using the values given in <xref ref-type="table" rid="tab5">Table 5</xref>, the optimization has been done, and the result obtained has been illustrated in the form of the Pareto front in <xref ref-type="fig" rid="fig11">Fig. 11</xref>.</p>
<fig id="fig11"><label>Figure 11</label><caption><p>Pareto front obtained using the DE algorithm.</p></caption><graphic content-type="pdf-region" xlink:href="pdf-region:11:66,667,537.1,780.3"/><p content-type="note">This figure is shown in the PDF version of the article.</p></fig>
<p>The step responses of the two operating modes of the system have been analyzed according to the obtained Pareto optimal solutions. These are shown in <xref ref-type="fig" rid="fig12">Figs. 12</xref> and 13, for the first and the second operating modes,</p>
<fig id="fig12"><label>Figure 12</label><caption><p>Step response of the first operating mode using the Pareto-optimal solution for the output torque (M1).</p></caption><graphic content-type="pdf-region" xlink:href="pdf-region:11:65.7,492.5,537.1,618"/><p content-type="note">This figure is shown in the PDF version of the article.</p></fig>
<fig id="fig13"><label>Figure 13</label><caption><p>Step response of the second operating mode using the Pareto-optimal solution for the output torque (M2).</p></caption><graphic content-type="pdf-region" xlink:href="pdf-region:11:66,353.1,537.1,465.2"/><p content-type="note">This figure is shown in the PDF version of the article.</p></fig>
<p>One Pareto-optimal solution has been chosen among all such solutions to compare with the previous controllers, such as the traditional PID controller and GA-optimized fuzzy controller. The optimized membership functions for the selected Pareto-optimal solution have been shown in <xref ref-type="fig" rid="fig14">Figs. 14</xref> and 15 for the output torques (M1) and (M2), respectively.</p>
<fig id="fig14"><label>Figure 14</label><caption><p>Optimized torque membership function of the output torque (M_1).</p></caption><graphic content-type="pdf-region" xlink:href="pdf-region:11:66,76.8,537.1,262.6"/><p content-type="note">This figure is shown in the PDF version of the article.</p></fig>
<fig id="fig15"><label>Figure 15</label><caption><p>Optimized torque membership function of the output torque (M2).</p></caption><graphic content-type="pdf-region" xlink:href="pdf-region:12:66,589.1,537.1,779.1"/><p content-type="note">This figure is shown in the PDF version of the article.</p></fig>
<p>The step responses of the closed-loop system are obtained based on the optimized torque membership functions and compared with the responses of the PID controller and GA-optimized fuzzy controller. The comparative results of the responses for the two operating modes are illustrated in <xref ref-type="fig" rid="fig16">Figs. 16</xref> and 17,</p>
<fig id="fig16"><label>Figure 16</label><caption><p>Comparative responses of the three controllers for the first operating mode.</p></caption><graphic content-type="pdf-region" xlink:href="pdf-region:12:65.7,322.5,537.1,525.9"/><p content-type="note">This figure is shown in the PDF version of the article.</p></fig>
<fig id="fig17"><label>Figure 17</label><caption><p>Comparative responses of the three controllers for the second operating mode.</p></caption><graphic content-type="pdf-region" xlink:href="pdf-region:12:66,122.1,537.1,308"/><p content-type="note">This figure is shown in the PDF version of the article.</p></fig>
<p>It can be seen from the figures above that the selected initial torque membership functions ensure adequate control performance. Nevertheless, the optimization with the DE method results in further improvement of the transient response through reduction of the overshoot and settling time. Hence, the proposed DE-based optimization improves the control system performance better than the GA-based one does.</p>
</sec>
<sec id="sec-19"><label>6</label><title>Conclusion</title>
<p>In this work, a FLC approach is proposed for the attitude tracking control of a nonlinear coupled system with two rotational degrees of freedom. The primary aim here was to enhance the tracking performance of the system while minimizing unwanted transients and effort exerted by the control action. Initially, a Mamdani-type FLC system was designed, and its performance was evaluated as compared to the classical PID controller. From the simulation results, it could be seen that the classical PID controller showed considerable overshoots and slow response time for both operating modes. The performance of this controller was thus unsatisfactory for the intended application of attitude tracking. For the first and second operating modes, the overshoot was 50% and 60% with response times of 10 and 15 seconds, respectively.</p>
<p>To address these constraints, a GA was then used to tune the fuzzy control rules. The GA tuned fuzzy controller showed considerable enhancement in the transient response in comparison to the PID controller. In the first operating mode, the overshoot value decreased from around 50% to below 20%, while the settling time dropped from about 10 s to below 7 s. In the second operating mode, the overshoot value decreased from around 60% to below 5%, while the settling time dropped from about 15 s to below 5 s. These results highlight the efficiency of the evolutionary tuning of fuzzy controllers.</p>
<p>Finally, a MODE algorithm was used to optimize the membership functions of the outputs corresponding to the control torques (M_1) and (M_2). Contrary to the previous single-objective optimization, the DE-based optimization considered both the tracking error and the control effort objectives at the same time. Therefore, the obtained non-dominated solutions form the Pareto front, representing various tradeoffs between the tracking error and control torque.</p>
<p>Furthermore, the comparative study of the chosen Pareto-optimal point with respect to the conventional PID and GA optimized fuzzy control methods showed that additional benefits could be attained through the output torque membership functions&apos; tuning process. In addition, the method based on DE improved such control performance characteristics as overshoot and settling time compared to the initially designed membership functions as well as GA based fuzzy control system. Thus, one may conclude that the DE technique represents an efficient tool for fuzzy controller parameters&apos; fine-tuning process while resolving contradictory control objectives at the same time.</p>
<p>On the whole, the suggested multi-objective fuzzy control system based on DE provides an effective approach for solving nonlinear 3D tracking problems in cases when the accuracy and control effort are the control objectives. The obtained simulation results prove that combining fuzzy control and the evolutionary algorithm leads to a significant increase in dynamic performance in comparison to conventional PID controllers and GA-based fuzzy optimization. Therefore, the suggested approach may be considered as an efficient framework for intelligent controller design for aerospace and other nonlinear tracking systems.</p>
</sec>
<sec id="sec-20"><title>Conflict of Interest</title>
<p>The authors declare no conflict of interest.</p>
</sec>
<sec id="sec-21"><title>Data Availability</title>
<p>All data are included in the text.</p>
</sec>
<sec id="sec-22"><title>Funding</title>
<p>This research received no specific grant from funding agencies in the public, commercial, or not-for-profit sectors.</p>
</sec>
</body>
<back><ref-list><title>References</title>
<ref id="ref-1"><label>1</label><mixed-citation>Zadeh, L. A. (1965). Fuzzy sets. Information and control, 8(3), 338–353. <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1016/S0019-9958(65)90241-X">https://doi.org/10.1016/S0019-9958(65)90241-X</ext-link></mixed-citation></ref>
<ref id="ref-2"><label>2</label><mixed-citation>Mamdani, E. H., &amp; Assilian, S. (1993). An experiment in linguistic synthesis with a fuzzy logic controller. In Readings in fuzzy sets for intelligent systems (pp. 283–289). Morgan Kaufmann. <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1016/B978-1-4832-1450-4.50032-8">https://doi.org/10.1016/B978-1-4832-1450-4.50032-8</ext-link></mixed-citation></ref>
<ref id="ref-3"><label>3</label><mixed-citation>Storn, R., &amp; Price, K. (1997). Differential evolution – a simple and efficient heuristic for global optimization over continuous spaces. Journal of global optimization, 11(4), 341–359. <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1023/A:1008202821328">https://doi.org/10.1023/A:1008202821328</ext-link></mixed-citation></ref>
<ref id="ref-4"><label>4</label><mixed-citation>Vesterstrom, J., &amp; Thomsen, R. (2004). A comparative study of differential evolution, particle swarm optimization, and evolutionary algorithms on numerical benchmark problems. Proceedings of the 2004 congress on evolutionary computation (IEEE cat. no.04th8753) (Vol. 2, pp. 1980-1987 Vol.2). IEEE. <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1109/CEC.2004.1331139">https://doi.org/10.1109/CEC.2004.1331139</ext-link></mixed-citation></ref>
<ref id="ref-5"><label>5</label><mixed-citation>Price, K., Storn, R., &amp; Lampinen, J. (2005). The differential evolution algorithm. In differential evolution: A practical approach to global optimization (pp. 37–134). Berlin, Heidelberg: Springer Berlin Heidelberg. <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1007/3-540-31306-0_2">https://doi.org/10.1007/3-540-31306-0_2</ext-link></mixed-citation></ref>
<ref id="ref-6"><label>6</label><mixed-citation>Das, S., &amp; Suganthan, P. N. (2011). Differential evolution: A survey of the state-of-the-art. IEEE transactions on evolutionary computation, 15(1), 4–31. <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1109/TEVC.2010.2059031">https://doi.org/10.1109/TEVC.2010.2059031</ext-link></mixed-citation></ref>
<ref id="ref-7"><label>7</label><mixed-citation>Pishkenari, H. N., Mahboobi, S. H., &amp; Meghdari, A. (2004). On the optimum design of fuzzy logic controller for trajectory tracking using evolutionary algorithms. IEEE conference on cybernetics and intelligent systems, 2004. (Vol. 1, pp. 660–665 vol.1). IEEE. <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1109/ICCIS.2004.1460494">https://doi.org/10.1109/ICCIS.2004.1460494</ext-link></mixed-citation></ref>
<ref id="ref-8"><label>8</label><mixed-citation>Pishkenari, H. N., Mahboobi, S. H., &amp; Alasty, A. (2011). Optimum synthesis of fuzzy logic controller for trajectory tracking by differential evolution. Scientia Iranica, 18(2), 261–267. <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1016/j.scient.2011.03.021">https://doi.org/10.1016/j.scient.2011.03.021</ext-link></mixed-citation></ref>
<ref id="ref-9"><label>9</label><mixed-citation>Chen, C. H. (2013). Design of TSK-type fuzzy controllers using differential evolution with adaptive mutation strategy for nonlinear system control. Applied mathematics and computation, 219(15), 8277–8294. <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1016/j.amc.2013.02.026">https://doi.org/10.1016/j.amc.2013.02.026</ext-link></mixed-citation></ref>
<ref id="ref-10"><label>10</label><mixed-citation>Ochoa, P., Castillo, O., &amp; Soria, J. (2020). Optimization of fuzzy controller design using a Differential Evolution algorithm with dynamic parameter adaptation based on Type-1 and Interval Type-2 fuzzy systems. Soft computing, 24(1), 193–214. <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1007/s00500-019-04156-3">https://doi.org/10.1007/s00500-019-04156-3</ext-link></mixed-citation></ref>
<ref id="ref-11"><label>11</label><mixed-citation>Ochoa, P., Castillo, O., Melin, P., &amp; Castro, J. R. (2023). Interval type-3 fuzzy differential evolution for parameterization of fuzzy controllers. International journal of fuzzy systems, 25(4), 1360–1376. <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1007/s40815-022-01451-4">https://doi.org/10.1007/s40815-022-01451-4</ext-link></mixed-citation></ref>
<ref id="ref-12"><label>12</label><mixed-citation>Marinaki, M., Marinakis, Y., &amp; Stavroulakis, G. E. (2015). Fuzzy control optimized by a multi-objective differential evolution algorithm for vibration suppression of smart structures. Computers &amp; structures, 147, 126–</mixed-citation></ref>
<ref id="ref-13"><label>137</label><mixed-citation><ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1016/j.compstruc.2014.09.018">https://doi.org/10.1016/j.compstruc.2014.09.018</ext-link></mixed-citation></ref>
<ref id="ref-14"><label>13</label><mixed-citation>Subburaj, B., Maheswari, J. U., Ibrahim, S. P. S., &amp; Kavitha, M. S. (2024). Population diversity control based differential evolution algorithm using fuzzy system for noisy multi-objective optimization problems. Scientific reports, 14(1), 17863. <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1038/s41598-024-68436-1">https://doi.org/10.1038/s41598-024-68436-1</ext-link></mixed-citation></ref>
<ref id="ref-15"><label>14</label><mixed-citation>Tolossa, T. D., Gunasekaran, M., Halder, K., Verma, H. K., Parswal, S. S., Jorwal, N., … &amp; Hote, Y. V. (2024). Trajectory tracking control of a mobile robot using fuzzy logic controller with optimal parameters. Robotica, 42(8), 2801–2824. <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1017/S0263574724001140">https://doi.org/10.1017/S0263574724001140</ext-link></mixed-citation></ref>
<ref id="ref-16"><label>15</label><mixed-citation>Awda, Y., &amp; Alowaifeer, M. (2024). Adaptive optimization of virtual synchronous generator based on fuzzy logic control and differential evolution. Ain Shams engineering journal, 15(4), 102606. <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1016/j.asej.2023.102606">https://doi.org/10.1016/j.asej.2023.102606</ext-link></mixed-citation></ref>
<ref id="ref-17"><label>16</label><mixed-citation>Wang, H., Wang, L., Zhang, Y., &amp; Niu, X. (2026). Optimized fuzzy feedback control for autonomous vehicle trajectory tracking via adaptive multi-objective differential evolution. Swarm and evolutionary computation, 104, 102365. <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1016/j.swevo.2026.102365">https://doi.org/10.1016/j.swevo.2026.102365</ext-link></mixed-citation></ref>
<ref id="ref-18"><label>17</label><mixed-citation>Chen, X., Dong, H., Shen, C., Li, H., &amp; Li, D. (2026). Optimizing interval type-2 fuzzy logic PID controller with an improved constraint differential evolution algorithm. Scientific reports, 16(1), 10426. <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1038/s41598-026-41203-0">https://doi.org/10.1038/s41598-026-41203-0</ext-link></mixed-citation></ref>
</ref-list></back>
</article>