# Optimal Fuzzy Proportional-Integral-Derivative Control for a Class of Fourth-Order Nonlinear Systems using Imperialist Competitive Algorithms

Authors: S. Hadipour Lakmesari; Z. Safipour; M. J. Mahmoodabadi; Yousef Ibrahim; Saleh Mobayen
Journal: Complexity
 (2022)

---

## Abstract

The proportional integral derivative (PID) controller has gained wide acceptance and use as the most useful control approach in the industry. However, the PID controller lacks robustness to uncertainties and stability under disturbances. To address this problem, this paper proposes an optimal fuzzy-PID technique for a two-degree-of-freedom cart-pole system. Fuzzy rules can be combined with controllers such as PID to tune their coefficients and allow the controller to deliver substantially improved performance. To achieve this, the fuzzy logic method is applied in conjunction with the PID approach to provide essential control inputs and improve the control algorithm efficiency. The achieved control gains are then optimized via the imperialist competitive algorithm. Consequently, the objective function for the cart-pole system is regarded as the summation of the displacement error of the cart, the angular error of the pole, and the control force. This control concept has been tested via simulation and experimental validations. Obtained results are presented to confirm the accuracy and efficiency of the suggested method.


---

## Body

## 1. Introduction

### 1.1. Background and Motivation

The proportional-integral-derivative (PID) control method is one of the most technical approaches in the control field, which has experienced remarkable development with regard to its use in mechanical systems [1–5]. In fact, this control approach is renowned for its simple structure. However, tuning its parameters is a challenge that requires carrying out some strategies such as classical tuning rules, adaptation laws, optimal techniques, and fuzzy systems [6–11]. Through these methods, fuzzy logic, which was initially introduced by Zadeh [12], is a very effective way to adjust and determine the PID control gains [13].


### 1.2. Literature Review

The improvement of fuzzy-PID methods for diverse mechanical systems has attracted attention in the last couple of years. Wang et al. designed a novel fuzzy-PID on a model of a coke furnace and compared its efficiency with conventional PID control [14]. Lapa and Cpalka proposed a nature-inspired technique to automate the determination of the parameters and structure of a flexible fuzzy-PID technique [15]. Kumar employed a new application of the genetic algorithm to optimize the parameters of type-2 fuzzy-PID controllers used for five-degrees of freedom robot manipulators [16]. Sa-ngiamvibool applied an optimal fuzzy-PID controller for hydrothermal systems [17]. Barik et al. [18] used a fuzzy proportional, derivative integral controller to investigate the frequency fluctuations due to wind and solar-thermal based renewable power [19]. Acharya and Das [20] introduced a fuzzy-PID controller to enhance the pressure tracking profile of pressure-controlled ventilators [16].In this respect, the heuristic parameters of such controllers should be defined by a proper method. One of the effective methods to determine these parameters is to take advantage of evolutionary algorithms such as ICA [21]. For instance, Ahmadi et al. used a novel scenario using the composition of the ICA and fuzzy-PID to solve the speed matter in an electric motor [22]. Cornejo et al. developed a tuning method for PID schemes based on evolutionary algorithms, particularly the ICA, the genetic algorithm (GA), and the active-set algorithm [23]. Kusama et al. used various artificial intelligence methods such as the ICA and particle swarm optimization to tune the PID controller gains [24]. Eventually, in [25], the ICA was utilized for a fractional-order fuzzy PI method to optimize the control gains.


### 1.3. Contributions

This paper develops some previous studies [26–30] by proposing a new optimal fuzzy-PID method that is introduced and executed on a given test-bed system. The innovation of this research is the novel control strategy presented, specifically the design of its fuzzy systems. In the proposed strategy, two fuzzy inference engines are applied to enhance the performance of the control approach. Next, the imperialist competitive algorithm is employed to optimize the proposed controller gains. Obtained results confirm the better time response and high stabilization accuracy for the presented approach in comparison to the other methods.


### 1.4. Organization of Paper

The study is structured and presented in seven sections. Section2 defines the dynamic model of the 2-DOF cart-pole system. Section 3 explains the fuzzy-PID approach. Section 4 presents the ICA algorithm. Section 5 demonstrates how to apply the proposed optimal fuzzy-PID control to the cart-pole system. Section 6 demonstrates the simulation and experimental outcomes for the cart-pole system. Section 7 presents discussion and conclusions.


## 1.1. Background and Motivation

The proportional-integral-derivative (PID) control method is one of the most technical approaches in the control field, which has experienced remarkable development with regard to its use in mechanical systems [1–5]. In fact, this control approach is renowned for its simple structure. However, tuning its parameters is a challenge that requires carrying out some strategies such as classical tuning rules, adaptation laws, optimal techniques, and fuzzy systems [6–11]. Through these methods, fuzzy logic, which was initially introduced by Zadeh [12], is a very effective way to adjust and determine the PID control gains [13].


## 1.2. Literature Review

The improvement of fuzzy-PID methods for diverse mechanical systems has attracted attention in the last couple of years. Wang et al. designed a novel fuzzy-PID on a model of a coke furnace and compared its efficiency with conventional PID control [14]. Lapa and Cpalka proposed a nature-inspired technique to automate the determination of the parameters and structure of a flexible fuzzy-PID technique [15]. Kumar employed a new application of the genetic algorithm to optimize the parameters of type-2 fuzzy-PID controllers used for five-degrees of freedom robot manipulators [16]. Sa-ngiamvibool applied an optimal fuzzy-PID controller for hydrothermal systems [17]. Barik et al. [18] used a fuzzy proportional, derivative integral controller to investigate the frequency fluctuations due to wind and solar-thermal based renewable power [19]. Acharya and Das [20] introduced a fuzzy-PID controller to enhance the pressure tracking profile of pressure-controlled ventilators [16].In this respect, the heuristic parameters of such controllers should be defined by a proper method. One of the effective methods to determine these parameters is to take advantage of evolutionary algorithms such as ICA [21]. For instance, Ahmadi et al. used a novel scenario using the composition of the ICA and fuzzy-PID to solve the speed matter in an electric motor [22]. Cornejo et al. developed a tuning method for PID schemes based on evolutionary algorithms, particularly the ICA, the genetic algorithm (GA), and the active-set algorithm [23]. Kusama et al. used various artificial intelligence methods such as the ICA and particle swarm optimization to tune the PID controller gains [24]. Eventually, in [25], the ICA was utilized for a fractional-order fuzzy PI method to optimize the control gains.


## 1.3. Contributions

This paper develops some previous studies [26–30] by proposing a new optimal fuzzy-PID method that is introduced and executed on a given test-bed system. The innovation of this research is the novel control strategy presented, specifically the design of its fuzzy systems. In the proposed strategy, two fuzzy inference engines are applied to enhance the performance of the control approach. Next, the imperialist competitive algorithm is employed to optimize the proposed controller gains. Obtained results confirm the better time response and high stabilization accuracy for the presented approach in comparison to the other methods.


## 1.4. Organization of Paper

The study is structured and presented in seven sections. Section2 defines the dynamic model of the 2-DOF cart-pole system. Section 3 explains the fuzzy-PID approach. Section 4 presents the ICA algorithm. Section 5 demonstrates how to apply the proposed optimal fuzzy-PID control to the cart-pole system. Section 6 demonstrates the simulation and experimental outcomes for the cart-pole system. Section 7 presents discussion and conclusions.


## 2. Cart-Pole System Formulation

The 2-DOF cart-pole system is regarded as a nonlinear under-actuated system. Indeed, it is a reputable test-bed to analyze control strategies [31, 32]. Generally, such systems include practical applications such as hypersonic vehicles, overhead cranes, rockets, and robots [33, 34]. The tracker’s design and control require thorough consideration. Therefore, stabilization of under-actuated systems has recently attracted many researchers’ attention. Accordingly, various control methods such as feedback linearization, sliding mode control, H∞ control, neural networks, and other robust control have been designed to apply proper controllers for these systems [26–30]. In fact, the under-actuated systems include many benefits, such as decreasing the construction cost, having a minimal number of actuators, and lightening the system.Consider a fourth-order nonlinear cart-pole system shown as follows [25]:(1)x˙1=x2x˙2=n1+h1u+rx˙3=x4x˙4=n2+h2u+rn1=M+mgsinx1−mLsinx1cosx1x22L4/3M+m−mcos2x1n2=mgsinx1cosx1−4/3mLx22sinx14/3M+m−mcos2x1h1=cosx1L4/3M+m−mcos2x1h2=4/34/3M+m−mcos2x1,where, x=x1,x2,x3,x4T is the state vector, nx,t and hx,t are nonlinear functions displaying system dynamics, and rt defines the uncertainty. x1,x2,x3,andx4 respectively represent the position and the velocity of the cart, the angle and the angular velocity of the pole. In addition, g,m,M, and L respectively show the gravity acceleration, the mass of the pole, the mass of the cart, and the half-length of the pole. The general schematic of the cart-pole system is demonstrated in Figure 1.Figure 1
Configuration of the considered system.In the next sections, we define general concepts of an optimal fuzzy-PID approach and then derive a controller for the 2-DOF cart-pole system based on fuzzy logic and an imperialist competitive algorithm.


## 3. Fuzzy-PID Controller

PID control is the most prevalent method among all feedback controls because of its simplicity and wide applications. As the method with a simple control structure that is straightforward to adjust, it was utilized by plant operators. The control inputs of the control scheme can be obtained through the following procedure [28]:(2)x¨=n3+h3u,θ¨=n4+h4u,where, u is the control input. Also, n3, n4,h3, and h4 are nonlinear functions. Consider the state-space formulation of the nonlinear system with the following: (3)x˙1=x2x˙2=x3x˙3=n3+h3ux˙4=x5x˙5=x6x˙6=n4+h4u,in which(4)x=x1,x2,x3,x4,x5,x6=∫x,x,x˙∫θ,θ,θ˙,Is the state vector while its desired value is as follows:(5)xd=∫xd,xd,xd˙∫θd,θd,θd˙,Accordingly(6)ex=x−xd,and(7)eθ=θ−θd,As the control inputs, the common form of the PID controller is written according to the following equations:(8)uPIDx=KPxext+Kix∫0τexτdτ+Kdxddtext,(9)uPIDθ=KPθeθt+Kiθ∫0τeθτdτ+Kdθddteθt,where, KPx and KPθ are the proportional gains, Kix and Kiθ denote the integral gains, and Kdx and Kdθ represent the derivative gains. The selection and tuning of the gain values are crucial issues in the design of PID controllers. Fuzzy rules have the capability to be combined with controllers such as PID to tune their coefficients and allow the controller to provide more precise answers. Thus, to provide adaptive gains for the PID controller, the fuzzy logic technique as a regulator method is used. The following equation defines fuzzy-PID control:(10)uFuzzy−PIDx=Kixfw1+Kpxfw2+Kdxfw3+Kiθfw4+Kpθfw5+Kdθfw6,where, wii = 1, 2, …, 6 represent the fuzzy variables containing inputs ∫x,x,x˙, ∫θ,θ,θ˙ achieved by first fuzzy system (FFS). Kixf,Kpxf,Kdxf,Kiθf,Kpθf, and Kdθf are obtained with respect to the following procedure:(11)Kixf=Kixc+ε1Δγ1Kpxf=Kpxc+ε2Δγ2Kdxf=Kdxc+ε3Δγ3Kiθf=Kiθc+ε4Δγ4Kpθf=Kpθc+ε5Δγ5Kdθf=Kdθc+ε6Δγ6,where, i = 1, 2, …, 6 denote input items ∫x,x,x˙ , ∫θ,θ,θ˙ achieved by second fuzzy system (SFS). εi,i=1,2,…,6aswellasKixc,Kpxc,Kdxc,Kiθc,Kpθc,andKdθc are weight constants that can be acquired by trial and error. Meanwhile, one effective technique to determine such parameters is evolutionary algorithms. Accordingly, in this paper, the imperialist competition algorithm (ICA) is used to determine the parameters of the proposed controller.


## 4. Imperialist Competition Algorithm

Considering the complicated world where we live, it can be divided into local places controlled by local powers and superpowers, which can have global scale influence. This pattern is more obvious when we look at that part of history when countries were colonized by more powerful imperialists. Generally, the more powerful country, the more colonies it could have. Therefore, there was a strict competition between powerful countries to conquest weaker countries and extended their territory. This model leads to an evolutionary algorithm termed the imperialist competition algorithm (ICA). Although most optimization algorithms are based upon modeling of natural processes, the ICA is inspired by imperialist competition and begins with initial populations, which are called countries. In this regard, several of the best countries would be chosen as the imperialist states, while the rest would be their colonies. If the optimization problem isN-dimensional, a country is a 1×N array, illustrated by the following equation:(12)Country=p1,p2,p3,…,pNvar,where, pi are variables that are optimized. Since our optimization problem is designing constant parameters of the fuzzy-PID method, the initial countries are those variables such as Kixc,Kpxc,Kdxc,Kiθc,Kpθc, and Kdθc:(13)Country=country11country12country21..countryyl.The function of the ICA is defined as looking for a country with the greatest socio-political specifications such as culture, religion, and economic. From an optimization standpoint, it gives rise to acquire an optimal answer of a problem, which has the minimum cost value. In this regard, each country’s cost can be determined as follows:(14)fcountry=fp1,p2,p3,…,pNvar.In this paper, there are 50 countries, which will result in 50 initial cost functions. Moreover, a country which has less cost is regarded as a powerful one; thus,Nimp of the strongest countries would be selected from empires while the remains are the colonies of these imperialists. With regard to initial cost functions, Nimp is determined as 10 in this paper.This competition is simulated by selecting some of the feeblest empires to start a rivalry out of all empires to govern the weak colonies. Starting the rivalry, imperialists tend to possess colonies. The most powerful empires are more likely to possess these colonies. This movement is simulated by shifting all colonies to imperialist ones. For random variablex, the following relation can be written as:(15)x∼U0,β×d.In whichd represents the distance of colony and imperialist while β denotes a number which is greater than 1. More precisely, β≫1 leads to keep out colonies from the imperialist state, and an approaching value to one for β decreases the seeking capability of the ICA. It is worth noting that β is regarded as 2 in this study that results to approach countries toward the global minimum. Meanwhile, approaching the imperialist, a colony may obtain a location with smaller cost in comparison to that. In such a situation, the colony and imperialist alter their locations. Subsequently, the ICA would proceed by the imperialist in the fresh location so that the colonies would be incorporated by the imperialist in the fresh location. This movement is depicted in Figure 2.Figure 2
Changing position of colonies to their related imperialist.Although the overall vigor of an empire is chiefly influenced by the vigor of an imperialist country, the vigor of colonies of a given empire affects the entire vigor of that empire too. The process can be formulated as follows:(16)T.C.n=Costimperialistn+ξmeanCostcoloniesofempiren.In whichT.C.n represents the overall cost of the n-th empire, as well as ξ denotes a positive constant less than 1. Through imperialistic competition, gradually the weaker empires’ power is reduced and the power of stronger ones is enhanced. Therefore, weak empires would miss out their colonies so that eventually, the most powerful empire would remain while others would collapse. On the other hand, one idea to end the process can be taken advantage of a number of maximum iterations that is regarded as 200 in this study. The suggested ICA is briefly explained in Figure 3 as a flowchart.Figure 3
Flowchart of the ICA.


## 5. Optimal Fuzzy-PID Control on the Cart-Pole System

This section describes the design of an optimal fuzzy-PID method based on two fuzzy inference engines for stabilization of the cart-pole system. The controller includes six inputs so that each of them contains one first fuzzy system (FFS) and second fuzzy system (SFS). These FFSs and SFSs are exclusively designed for the stabilization of the cart position and the pole angular. Regarding Equation (1), u as a control input is considered as u=11×uFuzzy−PID, and other parameters are regarded as L=0.5m,M=1kg,m=0.1kg,andg=9.8m/s2. The block diagram relevant to the fuzzy-PID controller is illustrated in Figure 4 indicates that each of the state variables x1,x2,x3,x4,x5,x6 related to the cart-pole system is compared with its desired value. These state variables (∫x,x,x˙, ∫θ,θ,θ˙) are normalized with control correction factors (CCFs) and then are fed into the fuzzy system. The values of correction factors related to ∫x,x,x˙, ∫θ,θ,θ˙ in FFS blocks are shown in Table 1 as they are considered a constant value for SFSs as π/3. Moreover, the absolute value of x5 is selected as the only input variable for all SFS blocks. In more details, in the first step, output values of FFS-i and SFS-i are multiplied, and then, by addition of all output values, the fuzzy-PID control output would be obtained. Eventually, the CCF block converts the output value of the suggested controller to the cart-pole control force u=ufuzzy−PID×CCF. In what follows, two fuzzy inference engines utilized in this paper is explained.Figure 4
Block diagram of the designed control scheme.Table 1
Values of correction factors corresponding to∫x,x,x˙, ∫θ,θ,θ˙ in FFS.
VariablesCorrection factorsNormalized factors∫x1X1x2.4X2x˙1X3∫θ1X4θ30°X5θ˙100°/sX6


### 5.1. First Fuzzy System (FFS)

Generally, every input of the fuzzy system includes one FFS and SFS. To calculate each FFS-i, fuzzy rules are regulated according to Table 2. The membership functions of FFS -i input including Sij are defined as triangular and trapezoidal, as Figure 5. Further, αi and βi shown in Table 2 are fuzzy outputs, which should be optimized. FFS -i outputs (fi) for the cart and pole are the same and calculated as follows:(17)fi=∑j=13Sijf^ij∑j=13Sij,i=1,2,3.Table 2
Fuzzy rules related to FFSs of the fuzzy system inputs for the cart-pole system.
IfXii=1,2,…,6ThenfiSi1αi<0Si20Si3βi>0Figure 5
Membership functions related to FFSs of fuzzy system inputs for the cart-pole system.


### 5.2. Second Fuzzy System (SFS)

As mentioned above, FFSs affect directly the output while for the determination of output valuesΔwi, each input section indirectly relates to one SFS. The inference the dynamic variables of SFSs, fuzzy rules are given in Table 3. After normalization, the absolute value of the pole angle would be determined as the single input value (x5). Figure 6 shows membership functions related to Dij which are defined triangular. Coefficients γi,Δi,δi,φi shown in Tables 3 and 4 are fuzzy system outputs that should be optimized. In addition, SFS-i outputs (Δwi) for the pole and cart are respectively calculated as follows:(18)Δwi=∑j=13Dijw^ij∑j=13Dij,i=1,2,3,(19)Δwi=∑j=13Dijw^ij∑j=13Dij,i=4,5,6.Table 3
Fuzzy rules of SFSs related to three dynamical variables of the pole.
IfX5ThenΔwii=1,2,3Di10Di20≤γi≤1Di30.5≤Δi≤1Figure 6
Membership functions of SFS related to system inputs.Table 4
Fuzzy rules of SFSs related to three dynamical variables of the cart.
IfX5ThenΔwii=4,5,6Di10≤δi≤0.5Di20≤φi≤1Di30Hence, after calculation of it is possible to obtain the SFS variable, and with respect tofi0 for the FFS variable, ufuzzy−PID as an output control force would be obtained as follows:(20)uFuzzy−PID=Kˇixf1+Kˇpxf2+Kˇdxf3+Kˇiθf4+Kˇpθf5+Kˇdθf6=∑i=1NwiDfi0,Finally,ufuzzy−PID would perform as a control effort on the cart-pole system after applying the control correction factor on it, as follows:(21)u=ufuzzy−PID×CCF.In the next step, the optimum design variables would be determined by means of the ICA algorithm mentioned in Section4. For this purpose, objective functions should be determined as the first step by the following equations:(22)Of1=∫0txt−xdtdt=∫0textdtOf2=∫0tθt−θdtdt=∫0teθtdtOf3=∫0tFtdt,where, xdt = θdt=0 denotes the desired setpoints while Of1,Of2, and Of3 respectively represent the error of the cart and angle position, as well as the system force control. To optimize all three objective functions, the final objective function is expressed as the summation of them as follows:(23)Of=Of1+Of2+Of3.


## 5.1. First Fuzzy System (FFS)

Generally, every input of the fuzzy system includes one FFS and SFS. To calculate each FFS-i, fuzzy rules are regulated according to Table 2. The membership functions of FFS -i input including Sij are defined as triangular and trapezoidal, as Figure 5. Further, αi and βi shown in Table 2 are fuzzy outputs, which should be optimized. FFS -i outputs (fi) for the cart and pole are the same and calculated as follows:(17)fi=∑j=13Sijf^ij∑j=13Sij,i=1,2,3.Table 2
Fuzzy rules related to FFSs of the fuzzy system inputs for the cart-pole system.
IfXii=1,2,…,6ThenfiSi1αi<0Si20Si3βi>0Figure 5
Membership functions related to FFSs of fuzzy system inputs for the cart-pole system.


## 5.2. Second Fuzzy System (SFS)

As mentioned above, FFSs affect directly the output while for the determination of output valuesΔwi, each input section indirectly relates to one SFS. The inference the dynamic variables of SFSs, fuzzy rules are given in Table 3. After normalization, the absolute value of the pole angle would be determined as the single input value (x5). Figure 6 shows membership functions related to Dij which are defined triangular. Coefficients γi,Δi,δi,φi shown in Tables 3 and 4 are fuzzy system outputs that should be optimized. In addition, SFS-i outputs (Δwi) for the pole and cart are respectively calculated as follows:(18)Δwi=∑j=13Dijw^ij∑j=13Dij,i=1,2,3,(19)Δwi=∑j=13Dijw^ij∑j=13Dij,i=4,5,6.Table 3
Fuzzy rules of SFSs related to three dynamical variables of the pole.
IfX5ThenΔwii=1,2,3Di10Di20≤γi≤1Di30.5≤Δi≤1Figure 6
Membership functions of SFS related to system inputs.Table 4
Fuzzy rules of SFSs related to three dynamical variables of the cart.
IfX5ThenΔwii=4,5,6Di10≤δi≤0.5Di20≤φi≤1Di30Hence, after calculation of it is possible to obtain the SFS variable, and with respect tofi0 for the FFS variable, ufuzzy−PID as an output control force would be obtained as follows:(20)uFuzzy−PID=Kˇixf1+Kˇpxf2+Kˇdxf3+Kˇiθf4+Kˇpθf5+Kˇdθf6=∑i=1NwiDfi0,Finally,ufuzzy−PID would perform as a control effort on the cart-pole system after applying the control correction factor on it, as follows:(21)u=ufuzzy−PID×CCF.In the next step, the optimum design variables would be determined by means of the ICA algorithm mentioned in Section4. For this purpose, objective functions should be determined as the first step by the following equations:(22)Of1=∫0txt−xdtdt=∫0textdtOf2=∫0tθt−θdtdt=∫0teθtdtOf3=∫0tFtdt,where, xdt = θdt=0 denotes the desired setpoints while Of1,Of2, and Of3 respectively represent the error of the cart and angle position, as well as the system force control. To optimize all three objective functions, the final objective function is expressed as the summation of them as follows:(23)Of=Of1+Of2+Of3.


## 6. Results and Discussion

In what follows, the system performance is considered in various conditions to challenge the proposed strategy. Hence, the results are obtained for three different initial values. For all cases, the specifications of the cart-pole system and the ICA parameters are written according to Tables5 and 6.Table 5
Specifications of the cart-pole system.
ParametersValueM0.1 kgL0.5 mM1 kgG9.8 m/s2dmax0Table 6
ICA parameters.
ParametersDescriptionValueMaxiterMaximum iteration200NcountryInitial population50Ββ coefficient2NimpThe number of imperialists10


### 6.1. Case Study 1

In the first case study, the initial condition isx0=∫x0,x0,x˙0,∫θ0,θ0,θ˙0=0,2,0,0,0,0 and compared with three existing studies [28–30] including Yi and Yabazaki, 2000; Mahmoodabadi et al. 2016; and Mahmoodabadi and Jahanshahi, 2016. In Figure 7, the settling time of the proposed method for the cart position is shorter than other methods. According to Figure 8, the settling time for the pole angle is also shorter than all three methods while it has a smaller undershoot in comparison with the two other approaches. Figure 9 displays the comparison of the control efforts between the proposed methods and other techniques.Figure 7
Cart position based on the suggested fuzzy-PID controller and the techniques presented in references [28–30].Figure 8
Pole angle based on the suggested fuzzy-PID controller and the techniques presented in references [28–30].Figure 9
Control efforts based on the suggested fuzzy-PID controller and the techniques presented in references [28–30].


### 6.2. Case Study 2

In the second case study, the initial conditioni is compared with Yi and Yabazaki’s approach. According to Figure 10, the settling time of the proposed approach for the cart position is 1.5s shorter than Yi and Yabazaki’s method. Similarly, Figure 11 indicates that the settling time of the suggested method for the pole angle also has a better value while it has greater overshoot. In addition, the highest absolute value of the control input of the presented method is 30N smaller than Yi and Yabazaki’s scheme, as is seen in Figure 12.Figure 10
Cart position based on the suggested fuzzy-PID controller and Yi and Yabazaki’s technique with initial condition0,0,0,0,π/6,0.Figure 11
Pole angle based on the suggested fuzzy-PID controller and Yi and Yabazaki’s technique with0,0,0,0,π/6,0 initial condition.Figure 12
Control efforts based on the suggested fuzzy-PID controller and Yi and Yabazaki’s technique with0,0,0,0,π/6,0 initial condition.


### 6.3. Case Study 3

In the third case study, the initial condition is and compared with Yi and Yabazaki’s method. According to Figures13 and 14, the proposed method has a better time response in comparison to those techniques. Figure 15 indicates the comparison of the control inputs between the two methods. Further, the optimum design variables of the proposed control approach for three case studies are given in Table 7.Figure 13
Cart position based on the suggested fuzzy-PID controller and Yi and Yabazaki’s technique with0,2,0,0,π/6,0 initial condition.Figure 14
Pole angle based on the suggested fuzzy-PID controller and Yi and Yabazaki’s technique with0,2,0,0,π/6,0 initial condition.Figure 15
Control efforts based on the suggested fuzzy-PID controller and Yi and Yabazaki’s technique with0,2,0,0,π/6,0 initial condition.Table 7
Design variables determined by means of the ICA for three case studies.
Design variablesValues1st case study2nd case study3rd case studyα1−1.06−1.08−0.99β10.91.10.9α20.99−1.07−1.06β20.911.020.92α3−0.9−0.9−1.01β30.991.11.07α4−1−0.95−1.1β40.971.091.1α5−0.96−1.1−0.99β51.031.11.1α6−0.9−1.06−1.1β61.010.90.9δ10.970.90.9ϕ10.460.550.55Kixc0.30.030.03ε1−0.03−0.03−0.03δ21.081.041.04of6.5813.2818.73CCF1112.112.1ϕ20.470.450.45Kpxc0.580.640.64ε20.260.320.26δ30.981.10.9ϕ30.550.550.55Kdxc0.120.120.12ε30.140.150.16γ10.500.550.55Δ111.071.1Kiθc2.93.173.17ε44.164.174.17γ20.530.450.45Δ21.030.90.9Kpθc2.572.872.87ε5−1.35−1.35−1.35γ30.450.450.45Δ30.90.90.9Kdθc2.111.731.73ε64.704.504.50


### 6.4. Case Study 4

In this section, the authors intend to examine the suggested approach through an experimental setup. The practical tests are executed by using the Simulink® toolboxes and MATLAB® Real-Time®. Figure16 displays the practical cart-pole system, which is created in University of Zanjan, in the hanging and inverted modes. The angular position of the pendulum and position of the cart are computed via E40S Autonics Company encoders. The PCI-1751 data acquisition card (PCI-1751 DAQ) is used for the collection of input-output data and formation of interrelation amid computer and practical cart-pole system. The PCI-1751 DAQ includes the digital-to-analog and analog-to-digital converters. The cart position and the angular position of the pole controlled by the fuzzy-PID technique are respectively shown in Figures 17 and 18. In more detail, the considered system is in the swing-up mode in the first 4.5 seconds. Then, the control signal switches to the balancing controller. As a result, the pendulum was in the upright situation after 4.5 seconds. As observed in Figure 17, the proposed controller is succeeded in bringing the pole from radian to 0 in about 8 seconds. The pendulum remains almost equal to zero afterward. Moreover, as it can be observed from Figure 18, the cart position changes rapidly in the swing-up period. Moreover, Figure 19 depicts the control effort obtained by the proposed control method during the experimental test. These experimental outcomes are compatible with the simulation ones that approve the efficiency of the offered method.Figure 16
Configuration of the practical cart-pole system ((a) the pendulum in hanging mode, (b) the pendulum in inverted mode).Figure 17
Pole angle based on the suggested fuzzy-PID controller.Figure 18
Cart position based on the suggested fuzzy-PID controller.Figure 19
Control effort based on the suggested fuzzy-PID controller.


## 6.1. Case Study 1

In the first case study, the initial condition isx0=∫x0,x0,x˙0,∫θ0,θ0,θ˙0=0,2,0,0,0,0 and compared with three existing studies [28–30] including Yi and Yabazaki, 2000; Mahmoodabadi et al. 2016; and Mahmoodabadi and Jahanshahi, 2016. In Figure 7, the settling time of the proposed method for the cart position is shorter than other methods. According to Figure 8, the settling time for the pole angle is also shorter than all three methods while it has a smaller undershoot in comparison with the two other approaches. Figure 9 displays the comparison of the control efforts between the proposed methods and other techniques.Figure 7
Cart position based on the suggested fuzzy-PID controller and the techniques presented in references [28–30].Figure 8
Pole angle based on the suggested fuzzy-PID controller and the techniques presented in references [28–30].Figure 9
Control efforts based on the suggested fuzzy-PID controller and the techniques presented in references [28–30].


## 6.2. Case Study 2

In the second case study, the initial conditioni is compared with Yi and Yabazaki’s approach. According to Figure 10, the settling time of the proposed approach for the cart position is 1.5s shorter than Yi and Yabazaki’s method. Similarly, Figure 11 indicates that the settling time of the suggested method for the pole angle also has a better value while it has greater overshoot. In addition, the highest absolute value of the control input of the presented method is 30N smaller than Yi and Yabazaki’s scheme, as is seen in Figure 12.Figure 10
Cart position based on the suggested fuzzy-PID controller and Yi and Yabazaki’s technique with initial condition0,0,0,0,π/6,0.Figure 11
Pole angle based on the suggested fuzzy-PID controller and Yi and Yabazaki’s technique with0,0,0,0,π/6,0 initial condition.Figure 12
Control efforts based on the suggested fuzzy-PID controller and Yi and Yabazaki’s technique with0,0,0,0,π/6,0 initial condition.


## 6.3. Case Study 3

In the third case study, the initial condition is and compared with Yi and Yabazaki’s method. According to Figures13 and 14, the proposed method has a better time response in comparison to those techniques. Figure 15 indicates the comparison of the control inputs between the two methods. Further, the optimum design variables of the proposed control approach for three case studies are given in Table 7.Figure 13
Cart position based on the suggested fuzzy-PID controller and Yi and Yabazaki’s technique with0,2,0,0,π/6,0 initial condition.Figure 14
Pole angle based on the suggested fuzzy-PID controller and Yi and Yabazaki’s technique with0,2,0,0,π/6,0 initial condition.Figure 15
Control efforts based on the suggested fuzzy-PID controller and Yi and Yabazaki’s technique with0,2,0,0,π/6,0 initial condition.Table 7
Design variables determined by means of the ICA for three case studies.
Design variablesValues1st case study2nd case study3rd case studyα1−1.06−1.08−0.99β10.91.10.9α20.99−1.07−1.06β20.911.020.92α3−0.9−0.9−1.01β30.991.11.07α4−1−0.95−1.1β40.971.091.1α5−0.96−1.1−0.99β51.031.11.1α6−0.9−1.06−1.1β61.010.90.9δ10.970.90.9ϕ10.460.550.55Kixc0.30.030.03ε1−0.03−0.03−0.03δ21.081.041.04of6.5813.2818.73CCF1112.112.1ϕ20.470.450.45Kpxc0.580.640.64ε20.260.320.26δ30.981.10.9ϕ30.550.550.55Kdxc0.120.120.12ε30.140.150.16γ10.500.550.55Δ111.071.1Kiθc2.93.173.17ε44.164.174.17γ20.530.450.45Δ21.030.90.9Kpθc2.572.872.87ε5−1.35−1.35−1.35γ30.450.450.45Δ30.90.90.9Kdθc2.111.731.73ε64.704.504.50


## 6.4. Case Study 4

In this section, the authors intend to examine the suggested approach through an experimental setup. The practical tests are executed by using the Simulink® toolboxes and MATLAB® Real-Time®. Figure16 displays the practical cart-pole system, which is created in University of Zanjan, in the hanging and inverted modes. The angular position of the pendulum and position of the cart are computed via E40S Autonics Company encoders. The PCI-1751 data acquisition card (PCI-1751 DAQ) is used for the collection of input-output data and formation of interrelation amid computer and practical cart-pole system. The PCI-1751 DAQ includes the digital-to-analog and analog-to-digital converters. The cart position and the angular position of the pole controlled by the fuzzy-PID technique are respectively shown in Figures 17 and 18. In more detail, the considered system is in the swing-up mode in the first 4.5 seconds. Then, the control signal switches to the balancing controller. As a result, the pendulum was in the upright situation after 4.5 seconds. As observed in Figure 17, the proposed controller is succeeded in bringing the pole from radian to 0 in about 8 seconds. The pendulum remains almost equal to zero afterward. Moreover, as it can be observed from Figure 18, the cart position changes rapidly in the swing-up period. Moreover, Figure 19 depicts the control effort obtained by the proposed control method during the experimental test. These experimental outcomes are compatible with the simulation ones that approve the efficiency of the offered method.Figure 16
Configuration of the practical cart-pole system ((a) the pendulum in hanging mode, (b) the pendulum in inverted mode).Figure 17
Pole angle based on the suggested fuzzy-PID controller.Figure 18
Cart position based on the suggested fuzzy-PID controller.Figure 19
Control effort based on the suggested fuzzy-PID controller.


## 7. Conclusions

In this research paper, a fuzzy proportional-integral-derivative controller was designed to handle the cart-pole system as a test-bed. More precisely, the objective was to tune and update the proportional-integral-derivative coefficients using fuzzy rules. In addition, the constant parameters of the fuzzy proportional-integral-derivative controller were optimized via the imperialist competitive algorithm. The objective function was regarded as the summation of the angle error of the pole, the error of the cart position, and the system control force. Furthermore, the proposed control approach was further examined by comparison with existing studies with three different initial conditions. The reported results confirmed that the introduced strategy is effectively able to control the given under-actuated system. In case study 1, the proposed method proposed the minimum settling time of the cart position and pole angle compared to References [26, 27]. As shown in case study 2, the suggested control strategy could reduce the settling time of the cart position around 50%. With respect to case study 3, the settling time of the cart position decreased around 43%. The results exhibited the appropriate performance of the proposed controller by showing an optimal smooth control input. Although the presented scheme is appropriate, there are some developments that can be still made. To design the fuzzy system, other membership functions, fuzzifiers, and defuzzifiers such as Gaussian membership function, pseudo-exponential membership function, triangular fuzzifier, Gaussian fuzzifier, and center of gravity defuzzifier can be also used. The suggested strategy can be applied for stabilizing other under-actuated systems such as ball and beam, seesaw, and ball and wheel systems. The introduced scenario can be combined with other linear and nonlinear control strategies to stabilize different higher-order dynamical systems. Other evolutionary algorithms such as bacterial foraging algorithm, whale optimization algorithm, particle swarm optimization, and ant colony among others can be utilized to select the controller parameters.



---

*Source: 3554897-2022-10-17.xml*