International Journal of Computational Intelligence Systems

Volume 14, Issue 1, 2021, Pages 1042 - 1052

Slime Mould Algorithm-Based Tuning of Cost-Effective Fuzzy Controllers for Servo Systems

Authors
Radu-Emil Precup1, *, ORCID, Radu-Codrut David1, ORCID, Raul-Cristian Roman1, ORCID, Emil M. Petriu2, ORCID, Alexandra-Iulia Szedlak-Stinean1, ORCID
1Department of Automation and Applied Informatics, Politehnica University of Timisoara, Bd. V. Parvan 2, 300223 Timisoara, Romania
2School of Electrical Engineering and Computer Science, University of Ottawa, 800 King Edward, Ottawa, ON K1N 6N5, Canada
*Corresponding author. Email: radu.precup@aut.upt.ro
Corresponding Author
Radu-Emil Precup
Received 17 August 2020, Accepted 5 March 2021, Available Online 16 March 2021.
DOI
10.2991/ijcis.d.210309.001How to use a DOI?
Keywords
Low-cost fuzzy control; Optimal tuning; Position control; Servo systems; Slime Mould Algorithm
Abstract

This paper suggests five new contributions with respect to the state-of-the-art. First, the optimal tuning of cost-effective fuzzy controllers represented by Takagi–Sugeno–Kang proportional-integral fuzzy controllers (TSK PI-FCs) is carried out using a fresh metaheuristic algorithm, namely the Slime Mould Algorithm (SMA), and a fuzzy controller tuning approach is offered. Second, a relatively easily understandable formulation of SMA is offered. Third, a real-world application of SMA is given, focusing on the optimal tuning of TSK PI-FCs for nonlinear servo systems in terms of optimization problems that target the minimization of discrete-time cost functions defined as the sum of time multiplied by squared control error. Fourth, using the concept of improving the performance of metaheuristic algorithms with information feedback models, proposed by Wang and Tan, Improving metaheuristic algorithms with information feedback models, IEEE Trans. Cybern. 49 (2019), 542–555, Gu and Wang, Improving NSGA-III algorithms with information feedback models for large-scale many-objective optimization, Fut. Gen. Comput. Syst. 107 (2020), 49–69, and Zhang et al., Enhancing MOEA/D with information feedback models for large-scale many-objective optimization, Inf. Sci. 522 (2020), 1–16, new metaheuristic algorithms are introduced in terms of inserting the model F1 in SMA and other representative algorithms, namely Gravitational Search Algorithm (GSA), Charged System Search (CSS), Grey Wolf Optimizer (GWO) and Whale Optimization Algorithm (WOA). Fifth, the real-time validation of the cost-effective fuzzy controllers and their tuning approach is performed in the framework of angular position control of laboratory servo system. The comparison with other metaheuristic algorithms that solve the same optimization problem for optimal parameter tuning of cost-effective fuzzy controllers suggestively highlights the superiority of SMA. Experimental results are included.

Copyright
© 2021 The Authors. Published by Atlantis Press B.V.
Open Access
This is an open access article distributed under the CC BY-NC 4.0 license (http://creativecommons.org/licenses/by-nc/4.0/).

1. INTRODUCTION

Nonlinear models characterize the majority of physical systems and industrial processes, inspiring others to develop suitable approaches for the analysis of nonlinear systems. The rapidly increasing research on fuzzy control that copes with nonlinear systems has taken place in recent decades. This is important because fuzzy control is a robust and inexpensive mathematical approach to control highly complex nonlinear or nonanalytical systems in many industrial applications, and can be considered as a relatively simple initial nonlinear and convenient control approach. As shown by Precup et al. [1] and Precup and David [2], respectively, the systematic design and tuning of fuzzy control systems is supported by analyses that include stability, controllability, observability, sensitivity and robustness. The optimal tuning of fuzzy controllers is combined with these analyses. These analyses are viewed in the context of one of the hot research topic in the recent years, namely the development of the improvement of fuzzy systems performance, owing to their extensively growing applications. To improve the performance of fuzzy control systems, type-1 fuzzy control systems were also extended and generalized to type-2 and type-3 ones especially focusing on interval ones in this regard.

Metaheuristic algorithms are successfully applied because they ensure higher performance and require lower computing capacity and time versus deterministic algorithms in several optimization problems. Such challenging optimization problems are those specific to the optimal (parameter) tuning of fuzzy (logic) controllers, where both the process and the controller are nonlinear and deterministic algorithms are not successful. The following metaheuristic algorithms have been applied most recently to the optimal tuning of fuzzy controllers in representative examples: adaptive weight Genetic Algorithm (GA) for gear shifting control [3], GA-based multiobjective optimization for electric vehicle powertrain control [4], GA for hybrid power systems control [5], engines control [6], energy management in hybrid vehicles [7], servo system control [2], wellhead back pressure control systems [8], micro-unmanned helicopter control [9], Particle Swarm Optimization (PSO) algorithm with compensating coefficient of inertia weight factor for filter time constant adaptation in hybrid energy storage systems control [10], set-based PSO algorithm with adaptive weights for optimal path planning of unmanned aerial vehicles [11], PSO algorithm for zinc production [12] and inverted pendulum control [13], hybrid PSO-Artificial Bee Colony algorithm for frequency regulation in microgrids [14], Imperialist Competitive Algorithm for human immunodeficiency control [15], Grey Wolf Optimizer (GWO) algorithms for sun-tracker systems [16] and servo system control [2], PSO, Cuckoo Search and Differential Evolution (DE) for gantry crane systems position control [17], Whale Optimization Algorithm (WOA) for vibration control of steel structures [18], Grasshopper Optimization Algorithm for load frequency control [19], DE for electro-hydraulic servo system control [20], Gravitational Search Algorithm (GSA) and Charged System Search (CSS) for servo system control [2].

Slime Mould Algorithm (SMA) is a fresh metaheuristic algorithm proposed in [21]. It mimics the oscillation mode of slime mould (i.e., the Physarum polycephalum) in nature when producing positive and negative propagation feedback in the path toward the food. Compared to other metaheuristic algorithms, it is proved in [21] that SMA exhibits improved exploratory and exploitation features.

Other popular metaheuristic algorithms with fresh results are hybrid PSO-GA [22] and hybrid GSA-GA [23] for constrained solutions, water cycle [24] and bat [25] algorithms for combinatorial optimization, GWO with hierarchical fuzzy operator [26], hierarchical GA multiobjective optimization of neural networks [27], Cross-Entropy algorithm for manufacturing processes [28], Water Circle algorithm for traveling salesman problem [29], Conflict Monitoring algorithm [30], Jaya optimization algorithm for biped robots [31] and load forecast [32] and Chemical Reaction algorithm for community detection [33].

Some of the most recently proposed metaheuristic algorithms concern both the improvement of other algorithms and the development of new ones. Suggestive approaches are the improvement of Nondominated Sorting Genetic Algorithm-III (NSGA-III) algorithm with adaptive mutation operator and the application of this improved NSGA-III to Big Data tuning [34], the analysis of the behavior of crossover operators in NSGA-III algorithm and the operator improvement such that to work in large-scale optimization [35], the improvement of DE algorithm by a selection mechanism and its application to fuzzy job-shop scheduling problems [36], and population extremal optimization algorithms applied to continuous optimization and nonlinear controller optimal tuning problems [3739].

An important and general improvement of metaheuristic algorithms in terms of modifying their structure by including information feedback models is proposed and formulated generally by Wang and Tan [40]. Six types of information feedback models are defined in [40], where individuals from previous iterations are selected in either a fixed or random manner and embedded in the update process of the algorithms. A detailed exemplification on PSO is offered in [40], but experimental applications to other representative algorithms are included as well. This concept is included successfully by Gu and Wang in NSGA-III [41] and Zhang et al. in multiobjective evolutionary algorithms based on decomposition in [42], and has a big potential to be incorporated in other algorithms.

The optimal tuning of fuzzy controller can also be seen, discussed and investigated in the more general and generous context of fuzzy modeling. The estimation capability of fuzzy systems depends on their rule parameters and structure of fuzzy sets. Up to now, many learning methods have been developed for optimizing rule parameters, i.e., which belongs to topic of optimal tuning of parameters of fuzzy model structures. Moreover, evolving fuzzy systems are also a popular topic as they also learn and optimally build the structure of fuzzy systems along with their parameters. Also, research is also focused on the tuning rules, which can be extracted such that to adapt the structure of the fuzzy sets, but these approaches usually also require that the output derivative of fuzzy systems to be computed with respect to the parameters which are actually learned and specifically tuned.

This paper is built upon authors’ recent papers on the optimal tuning of Takagi–Sugeno–Kang proportional-integral fuzzy controllers (TSK PI-FCs) [1,2,4345] applied to servo systems control such that to obtain a reduced parametric sensitivity (with respect to process gain and time constants), and suggests a novel SMA-based tuning approach. The approach is focused on the position control of nonlinear servo systems viewed as controlled processes. SMA is involved in solving a minimization-type optimization problem, which involves a cost function that is equal to the sum of time multiplied by squared control error, representing the discrete-time version of the Integral of Time Multiplied by Squared Error (ITSE) used in continuous-time control.

The paper suggests the following new contributions:

  • a fuzzy controller tuning approach,

  • a relatively easily understandable formulation of SMA,

  • a real-world application of SMA,

  • using the concept of improving the performance of metaheuristic algorithms with information feedback models [4042], new metaheuristic algorithms are proposed by inserting the model F1 in SMA and other representative algorithms, namely GSA, CSS, GWO and WOA,

  • the real-time validation of cost-effective fuzzy controllers and their tuning approach.

These contributions are significant and also advantageous in the context of the state-of-the-art briefly discussed in this section as the superiority of SMA versus other optimization algorithms, namely PSO, GSA, CSS, GWO and WOA, is proved by means a comparison included in the paper. The comparison is supported by experimental results on a real-world application, namely angular position control of nonlinear servo system; real data from authors’ lab is included.

The next sections in the sequel are structured as follows: the optimization problem is defined in Section 2, and the associated models of the process and fuzzy controller are specified. SMA and, based on this algorithm, the new tuning approach that produces optimal TSK PI-FCs are presented in Section 3. Starting with the introduction of the information feedback model F1 in PSO in [40], details on the new metaheuristic algorithms obtained by inserting the information feedback model F1 in GSA, CSS, GWO and WOA are given in Section 4. The tuning approach is validated in Section 5 in terms of offering experimental results and a comparison as well. The conclusions are drawn in Section 6.

2. PROCESS MODELS, CONTROLLER MODELS AND OPTIMIZATION PROBLEM

The control system structure is illustrated in Figure 1. Figure 1 represents an extended version of the block diagram and description of sub-systems given in [46], highlighting that a continuous-time process is controlled (t0 indicates the continuous time) by a digital control (tdZ,td0 is the discrete-time index).

Figure 1

Setpoint filter type two-degree-of-freedom (2-DOF) fuzzy control system structure.

The sub-systems and variables involved in Figure 1 are [46]: F – the setpoint (or the reference input) filter, FC – the fuzzy controller, which is represented by the TSK PI-FC, PWM-DAC – the Pulse Width Modulated (PWM) Digital to Analog Converter that produces the analog voltage applied to the process P (in fact, the actuator), ADC – the Analog to Digital Converter needed for digital process control, r – the reference input, r1 – the reference input filtered through F, y – the controlled output, u – the control signal, and e=ry – the control error. The sub-system F is included in Figure 1 because this control system is a setpoint filter type two-degree-of-freedom (2-DOF) structure, which usually operates with Proportional-Integral (PI) or Proportional-Integral-Derivative (PID) linear controllers, but nonlinear TSK PI-FCs will be used in this paper. Disturbance inputs are not included in Figure 1 because the integral component of the controller copes with load-type disturbance inputs. The presentation in this section is based on authors’ past papers as, e.g., [2,4346], therefore figures and equations are reused in this regard.

The servo system that plays the role of P in Figure 1 is characterized by the continuous-time state space model

m(t)=1,if u(t)ub,u(t)+uc/ucuc,if uc<u(t)<uc,0,if uc|u(t)|ua,u(t)ua/ubua,if ua<u(t)<ub,1,if u(t)ub,x˙1(t)x˙2(t)=0101/TΣx1(t)x2(t)+0kP/TΣm(t),y(t)=10x1(t)x2(t)T,(1)
where the dynamics of both measurement instrumentation and actuator are included in P, the parameters and variables represent [46]: kP>0 – the servo system (or process) gain, TΣ>0 – the small time constant of the process, u(t) – the control signal as a PWM signal within the interval (due to the PWM duty cycle) 1u(t)1, x1(t)=α(t) (rad) – the angular position, x2(t)=ω(t) (rad/s) – the angular speed, m(t) – the output of the saturation and dead zone static nonlinearity, T – matrix transposition, and the parameters of the nonlinearity fulfill the conditions 0<ua<ub and 0<uc<ub.

As shown in [2,4446], the nonlinearity in (1) is neglected in order to enable the cost-effective linear and nonlinear (including fuzzy) controller design. The transfer function of this simplified model of P, which is next used in the controller design and tuning, is

P(s)=kEP/s1+TΣs,(2)
where kEP is the equivalent process gain:
kEP=kP/(ubuu),if ub<u(t)<uu,kP/(ubua),if ua<u(t)<ub.(3)

The servo system process parameters are: kEP=170,TΣ=0.76 s,ua=0.15,ub=1 and uu=0.15.

PI controllers are recommended in [47,48] for the process modeled in (2). The transfer function of the PI controller (to be replaced with FC in Figure 1) is

C(s)=kc1+sTi/s=kC1+1/sTi,kC=kcTi,(4)
where kC>0 or kc>0 is the controller gain and Ti>0 is the controller integral time constant.

The Extended Symmetrical Optimum (ESO) method suggested in [47,48] is recommended to be applied in order to carried out the PI controller tuning as it guarantees a trade-off to a set of empirical control system quality indices: percent overshoot, settling time, rise time and phase margin. That trade-off is ensured by a single design parameter with the notation β and the recommended domain 1<β20. The PI tuning conditions specific to the ESO method are [47,48]

kc=1/ββkEPTΣ2,Ti=βTΣ,kC=1/βkEPTΣ,(5)
and the transfer function of the setpoint filter F, which additionally helps the control system performance enhancement by canceling a zero in the closed-loop transfer function (in the linear case) with respect to the setpoint, is
F(s)=1/1+βTΣs.(6)

The quality of the control system is increased by inserting the cost-effective nonlinear controller represented by TSK PI-FC in terms of the structure given in Figure 1. Incidentally, the design and tuning start with employing the knowledge and experience gained in terms of experimenting PI controllers in various applications. The structure and input membership functions of TSK PI-FC are presented in Figure 2, where [46]: q−1 – the backward shift operator, TISO-FC – the Two Inputs-Single Output Fuzzy Controller sub-system (it is actually a nonlinear sub-system without dynamics as it is introduced outside), e(td) – control error, u(td) – control signal, Δe(td)e(td) variation, and Δu(td)u(td) variation.

Figure 2

Takagi-–Sugeno-–Kang proportional-integral fuzzy controllers (TSK PI-FC) configuration.

The expressions of the two increments illustrated in Figure 2 are obtained by Tustin’s discretization method that leads to the recurrent equation of the incremental discrete-time PI controller

Δu(td)=KPΔe(td)+μe(td),(7)
where KP and µ are
KP=kcTiTs/2,μ=2T s/2TiTs,(8)
where Ts>0 is the sampling period, which has to be set such that to have a quasi-continuous digital control implementation.

As considered in [47], TISO-FC in Figure 2 uses for defuzzification a weighted average based approach, while the inference engine employs SUM and PROD operators. The rule base of this sub-system consists of only two rules for the sake of cost-effective fuzzy control design and implementation [2,4346], and it is expressed in Table 1, where the parameter η,0.25<η<0.75, is inserted for overshoot diminution.

Δe(td)e(td)

N ZE P
P KPΔe(td)+μe(td) KPΔe(td)+μe(td) ηKPΔe(td)+μe(td)
ZE KPΔe(td)+μe(td) KPΔe(td)+μe(td) KPΔe(td)+μe(td)
N ηKPΔe(td)+μe(td) KPΔe(td)+μe(td) KPΔe(td)+μe(td)
Table 1
Rule base of TISO-FC.

The parameter η makes the difference in the rule consequents, which consist of two different linear discrete-time PI controllers. This structure of the fuzzy controller makes it behave as a bumpless interpolator between the two linear PI controllers in the rule consequents.

The modal equivalence principle [49] is applied to express one tuning equation for TSK PI-FC

BΔe=μBe,(9)
so the parameters of TSK PI-FC are grouped in the parameter vector ρ:
ρ=βBeηT(10)
of only three elements (the first one from the linear part and the last two ones from the TISO-FC part), illustrating that cost-effective fuzzy control is targeted.

As specified in Section 1, the optimization problem is defined as

ρ*=argmin ρDρJ(ρ),J(ρ)=td=0tde2(td,ρ),(11)
where J is the cost function, which is expressed as the sum of discrete times multiplied by squared control errors, ρ* is the optimal parameter vector
ρ*=β*Be*η*T,(12)
and Dρ is the feasible domain of ρ. Several recommendations are given in [2,4346] to set to set the domain Dρ, but the most important ones concern technical and economical issues including stability constraints [5055]. Although the stability analysis appears to be the source of many other synthesis issues, it is always important so less conservative stability conditions need to be established. Therefore, one of the first objectives of the stability investigation is to evaluate the maximum permissible parameter variations in order to sustain the delay model’s stability. In addition, that should be analyzed and discussed in relation to the accuracy of the controlled process model in (1) and also the accuracy of the cost function model in (11). Different approaches are used in case of model-free data-driven control. Following, SMA is inserted in solving (11).

3. SMA AND FUZZY CONTROLLER TUNING APPROACH

The operating mechanism of SMA described in accordance with the standard formulation introduced in [21] starts with the random initialization of the population of agents, i.e., agent positions in the slime mould, such that to belong to the search domain Ds. A total number of N agents (i.e., elements of the slime mould) is used, and each agent is assigned to a position vector Xi(k)

Xi(k)=xi1(k)xif(k)xiq(k)TDs,i=1N,(13)
where [46] xif(k) is the position of ith agent in fth dimension, f=1q, k is the index of the current iteration, k=1kmax, and kmax is the maximum number of iterations. The expression of the search domain is
Ds=l1,u1××l f,u f××lq,uqq,(14)
where lf are the lower bounds, uf are the upper bounds, f=1q, and
xif(k)l f,u f,i=1N,f=1q.(15)

The search process specific to SMA continues with approaching, wrapping and grabbling the food. This is modeled by several equations. Using the notation Si(k) for the fitness (i.e., cost function in the context of (11)) of ith agent with the position vector Xi(k), the population is ranked in the ascending order of the fitness function values (i.e., cost function values in the context of (11)) for optimization problems that target the minimization of the cost functions. The following sets of position vectors are obtained:

{X1(k),X2(k),,XH(k),XH+1(k),,XN(k)}=SetFHSetSH,SetFH=X1(k),X2(k),,XH(k),SetSH=XH+1(k),,XN(k),(16)
where
S1(k)S2(k)SH(k)SH+1(k)SN(k),(17)
H=N/2 is the integer part of N/2, SetFH is the ranked set of agents in the first half of the population, and SetSH is the ranked set of agents in the second half of the population.

The vector Wi(k) of weights of slime mould consists of the elements (weights) wif(k), f=1q

Wi(k)=wi1(k)wif(k)wiq(k)T,i=1N,(18)
which are computed in terms of
wif(k)=1+rifloggi(k),if Xi(k)SetFH,1+rifloggi(k),otherwise,gi(k)=Sb(k)Si(k)/Sb(k)Sw(k)+ε+1,i=1N,f=1q,(19)
where log is the notation for the decimal logarithm, rif are random numbers, 0rif1, Sb(k) is the best fitness obtained at the current iteration, Sw(k) is the worst fitness obtained at the current iteration, with the values (using (17))
Sb(k)=mini=1NSi(k)=S1(k),Sw(k)=maxi=1NSi(k)=SN(k),(20)
and ε,ε=const>0, is relatively small to avoid zero denominator in (19).

The parameters involved in the position update equation are computed in terms of

a=arctanh1k/kmax,b=1k/kmax,(21)
which allow the computation of the uniformly distributed random numbers vaf and vbf
vafa,a,vbfb,b,f=1q.(22)

Using the notation Xb(k) for the best agent, i.e., the solution obtained in all iterations so far, with the elements xbf(k), f=1q:

Xb(k)=xb1(k)xbf(k)xbq(k)T,(23)
and introducing two agents XA(k) and XB(k) randomly selected in the population:
XA(k)=xA1(k)xA f(k)xAq(k)T,XB(k)=xB1(k)xB f(k)xBq(k)T,(24)
with randomly selected A and B, A,B=1N, the position update equation is
xi f(k+1)=rfu fl f+l f,if r f<z,xb f(k)+vaf(wi f(k)    xA f(k)xB f(k)),if r fz and ri f<pi,vb fxi f(k),otherwise,i=1N,(25)
where rf are random numbers, 0rf1, z is a parameter with recommended values in [21] in the interval 0z0.1, and the parameters pi are computed in terms of [21]
pi=tanh|Si(k)Smin(k)|,i=1N,(26)
and Smin(k) is the best fitness obtained in all iterations so far:
Smin(k)=mini=1N,kb=0...kSi(kb).(27)

SMA consists of steps SMA1 to SMA7. These steps are given as follows after the revision of the algorithm presented in [21]:

  • SMA1. The initial random population, which consists of the position vectors Xi(k)Ds of N agents, i.e., i=1N, is generated such that to fulfill (13)(15). The maximum number of iterations is set to kmax. The iteration index is initialized to k=0. The best fitness Sbest is initialized to Sbest= because the optimization problem defined in (11) is a minimization-type one.

  • SMA2. The performance of each element of the population of agents (i.e., the slime mould) is evaluated; this evaluation involves simulations and/or experiments conducted on the fuzzy control systems if the optimal tuning of fuzzy controllers is carried out. The evaluation leads to the fitness function values of all agents.

  • SMA3. The best fitness Smin(k) obtained in all iterations is computed using (27). If Smin(k)<Sbest, then the best fitness is Sbest will be updated to Smin(k) and the corresponding agent will represent the updated best agent (vector) Xb(k).

  • SMA4. The vector Wi(k) of weights of slime mould is calculated using (18)(20).

  • SMA5. Agents’ positions are updated using (25), which is supported by (21)(24), (26) and (27).

  • SMA6. The iteration index k is incremented with one. The algorithm goes on with step 2 until kmax is reached.

  • SMA7. The solution to the optimization problem is the last best agent obtained so far, i.e., the vector Xb(kmax), and the minimum value of the cost function is Sbest.

SMA is used in the optimization problem defined in (11) in terms of the following relationships:

Xi(k)=ρ,i=1N,Si(k)=J,i=1N,Ds=Dρ,Xb(kmax)=ρ*,q=3.(28)

Using the aspects presented in this section and the previous one, the SMA-based tuning approach dedicated to TSK PI-FCs consists of proceeding the steps 1–4. These steps are briefly described as follows:

Step 1. As pointed out in [46], the sampling period Ts is set in accordance with the requirements of quasi-continuous digital control.

Step 2. The feasible domain Dρ is set in order to include all constraints imposed to the elements of ρ.

Step 3. SMA is integrated using steps SMA1 to SMA7 and (28) in solving the optimization problem defined in (11) leading to the optimal parameter vector ρ* in (12) and three of the optimal parameters of TSK PI-FC, namely β*, Be* and η*.

Step 4. The optimal parameter BΔe of TSK PI-FC, with the notation BΔe*, results from (9) calculated for the optimal parameters in (12) substituted in (8). The setpoint filter with the transfer function defined in (6) is used in the fuzzy control system structure using β* instead of β.

4. INSERTING INFORMATION FEEDBACK MODEL F1 IN METAHEURISTIC ALGORITHMS

Using the notation in (13) for the position vector Xi(k) of ith agent at iteration k of the metaheuristic algorithm without information feedback model (i.e., the initial algorithm), the notation Si(k) for its fitness, which is also equal to the cost function value J in terms of (28), and inserting the notation Zi(k) for the position vector of ith agent at iteration k of the metaheuristic algorithm with information feedback model F1

Zi(k)=zi1(k)zif(k)ziq(k)TDs,i=1N,(29)
where zif(k) is the position of ith agent in fth dimension, f=1q, the information feedback model F1 is characterized by the recurrent equation [4042]
Zi(k+1)=αF1Xi(k)+βF1Zi(k),(30)
where the expressions of the weighting factors αF1 and βF1, which satisfy αF1+βF1=1, are [4042]
αF1=Si(k)/(Si(k)+Si(k+1)),βF1=Si(k+1)/(Si(k)+Si(k+1)).(31)

Equations (30) and (31) are successfully integrated in PSO algorithms [40], in NSGA-III algorithms [41], and in multiobjective evolutionary algorithms based on decomposition [42]. The agent (or particle) velocity and position update equations employed in PSO algorithms are [2]

Vi(k+1)=w(k)Vi(k)+c1r1(Pg,bestXi(k))+c2r2(Pi,bestXi(k)),Xi(k+1)=Xi(k)+Vi(k+1),i=1N,(32)
where Vi(k) is ith agent particle (agent) velocity vector at iteration k,
Vi(k)=vi1(k)vif(k)viq(k)T,i=1N,(33)

vif(k) is the velocity of ith agent in fth dimension, f=1q,r1,r2 are random variables with uniform distribution between 0 and 1, c1,c2>0 are weighting factors of the stochastic acceleration terms that pull each particle toward their end position, w(k)>0 is the inertia weight, Pi,best is the best position (vector) achieved by ith agent (particle), and Pg,best is the best position (vector) explored by the entire swarm so far. The PSO algorithm with information feedback model F1 will be referred to as PSO-F1 algorithm.

The position and velocity updates in GSA are governed by the recurrent equations [2]

vi f(k+1)=ρivi f(k)+ai f(k),xi f(k+1)=xi f(k)+vi f(k+1),i=1N,f=1q,(34)
where ρi is a uniform random variable, 0ρi1, and aif(k) is the acceleration of ith agent at iteration k. GSA with information feedback model F1 will be referred to as GSA-F1 algorithm.

The new position and velocity vectors of ith agent (or charged particle) in CSS algorithms are calculated in terms of [2]

Xi(k+1)=ri1ka(k)(Fi/mi)(Δk)2+ri2kv(k)Vi(k)Δk+Xi(k),Vi(k+1)=Xi(k+1)Xi(k)/Δk,i=1N,(35)
where ka(k) is the acceleration parameter, kv(k) is the velocity parameter, ri1,ri2, are uniformly distributed random numbers within [0, 1], mi is the mass of ith agent, and Δk is the time step set to 1. The CSS algorithm with information feedback model F1 will be referred to as CSS-F1 algorithm.

The updated agent (i.e., grey wolf) position vector in GWO algorithms is obtained in terms of [2]

Xi(k+1)=(Xα(k+1)+Xβ(k+1)+ Xδ(k+1))/3,i=1N,(36)
where Xα(k), Xβ(k) and Xδ(k) are the top three agents positions vectors computed in terms of [2]
J(Xα(k))=mini=1N{J(Xi(k))|Xi(k)Ds},J(Xβ(k))=mini=1N{J(Xi(k))|Xi(k)Ds\{Xα(k)}},J(Xδ(k))=mini=1N{J(Xi(k))|Xi(k)Ds\{Xα(k),Xβ(k)}}.(37)

The GWO algorithm with information feedback model F1 will be referred to as GWO-F1 algorithm.

Using the notation Xbest(k) for the best agents’ position (i.e., the position of the prey tracked by the whales, which are the agents) vector in the population at iteration k, the position update equation in WOA is [45]

Xi(k+1)=Xbest(k)fi(k),if C1,Xrand(t)gi(k),if C2,hi(k)eb lcos(2πl)+Xbest(k),otherwise,C1:p<0.5 and |aif(k)|1f=1q,C2:p<0.5 and |aif(k)|>1f=1q,i=1N,(38)
where p is random number within 0p1, the parameter b = const defines the shape of the logarithmic spiral of whales’ path, set to b=1 for simplifying the shape, l is a uniformly distributed random number, 1l1, rand is an arbitrary agent index, and fi(k),gi(k) are distance vectors. WOA with information feedback model F1 will be referred to as WOA-F1 algorithm.

5. VALIDATION AND COMPARISON

SMA with information feedback model F1 (SMA-F1) and the tuning approach presented in Section 3 are validated in this section by the optimal tuning of TSK PI-FCs for the angular position control of the experimental setup described in [56], which belongs to the Intelligent Control Systems Laboratory of the Politehnica University of Timisoara, Romania. A photo of the experimental setup is given in Figure 3. The parameter values of the servo system models given in (1) and (2) are ua=0.15, ub=1, uu=0.15, kP=kEP=140 and TΣ=0.92 s. The application of the three steps of the tuning approach is described as follows.

Figure 3

Experimental setup.

The sampling period was set to Ts=0.01 s in step 1 of the tuning approach. The vector variables ρ of the cost function belong to the feasible (and also search) domain Dρ set in step 2 as

Dρ=β|3β17×Be|20Be38×η|0.25η0.75.(39)

The dynamic regimes considered in the optimization problem are characterized by the r=40 rad step type modification of the setpoint and zero initial conditions. The upper limit of the sum in (11) was set to 2000 instead of . Therefore, the time horizon used in the evaluation of the cost function by experiments conducted on the fuzzy control system is 20 s.

The parameters of SMA were set in step 3 using the recommendations given in [21], namely z=0.03 and ε=0.001. In addition, authors’ experience in other metaheuristic algorithms applied to optimal tuning of fuzzy controllers and fuzzy models as well is reflected in setting the parameters N=20 and kmax=20 to achieve an acceptable trade-off to convergence and computational resources. Fast convergence is aimed in this paper versus our recent papers [2,4346], where kmax=100 was imposed.

The simulations and experiments were run in Matlab R2007b environment on a personal computer with Windows 10 operating system and the following hardware configuration: CPU Intel Core i5-7500, 3.4 GHz quad-core, 16 GB of DDR3 RAM memory at 1600 MHz, 256 SSD. These details are convenient for other researchers to redo the experiments, thus making the work in this paper relatively easily acceptable.

The optimal controller parameters and the corresponding minimum cost function Jmin in relation with (11) are presented in Table 2. Table 2 includes suggestively the comparison with other metaheuristic algorithms, namely PSO, GSA, CSS, GWO, WOA, PSO with information feedback model F1 (PSO-F1), GSA with information feedback model F1 (GSA-F1), CSS with information feedback model F1 (CSS-F1), GWO with information feedback model F1 (GWO-F1) and WOA with information feedback model F1 (WOA-F1). The same numbers of agents N=20 and maximum number of iterations kmax=20 as in case of SMA and SMA-F1 were used for all algorithms in the comparison. In addition, also to ensure a fair comparison of the algorithms, in order to reduce the effects of random parameters in the algorithms (always present in metaheuristic optimization algorithms), the results are presented as averaged values after the first Nt=10 runs (trials) of the algorithms; the comparisons included in the past papers were based on the best five runs of algorithms. The values of the other nonrandom parameters of the algorithms used in the comparison are presented in [2,4446]. The readers are invited to contact the authors if they are interested in the implementations of all these algorithms.

Algorithm

Be*

η*

Be*

β*

kC*

Ti*(s)

Jmin

SMA 38 0.75 0.00857633 4.821519 0.0035384 4.435799 2797210
PSO 36.5966 0.5935440 0.0819451 4.913505 0.00352454 4.520424 3523353
GSA 33.0287 0.7414017 0.0786853 4.770794 0.00360137 4.398443 3693850
CSS 23.1039 0.6662344 0.0715622 3.776415 0.00407048 3.4743 4935620
GWO 37.9457 0.6994276 0.0857001 4.910214 0.00352702 4.517395 3315304
WOA 35.0330 0.6513980 0.0748931 5.09787 0.0034414 4.682604 3297086
SMA-F1 38 0.75 0.0857638 4.821494 0.00353585 4.435774 2797210
PSO-F1 27.5769 0.6315154 0.0601471 5.161258 0.00341579 4.748359 3905199
GSA-F1 32.6942 0.7219313 0.0685116 5.194281 0.00340791 4.778737 3444169
CSS-F1 20.8791 0.7095039 0.0464305 4.962937 0.00351122 4.565903 4213287
GWO-F1 34.9617 0.7367405 0.0770037 4.030243 0.00350845 4.550747 3355471
WOA-F1 34.2792 0.6915126 0.1020510 3.799553 0.00407587 3.495585 5236392

CSS, Charged System Search; CSS-F1, Charged System Search with information feedback model F1; GSA, Gravitational Search Algorithm; GSA-F1, Gravitational Search Algorithm with information feedback model F1; GWO, Grey Wolf Optimizer; GWO-F1, Grey Wolf Optimizer with information feedback model F1; PSO, Particle Swarm Optimization; PSO-F1, Particle Swarm Optimization with information feedback model F1; SMA, Slime Mould Algorithm; SMA-F1, Slime Mould Algorithm with information feedback model F1; WOA, Whale Optimization Algorithm; WOA-F1, Whale Optimization Algorithm with information feedback model F1.

Table 2
Average values of controller parameters and minimum cost function J.

The performance of SMA and other algorithms considered in Table 2 was assessed using two quality indices. The first index reflects the convergence speed cs by measuring the number of evaluations required for the cost function J(ρ) until obtaining the optimal solution ρ*. The second index reflects the accuracy rate ar defined as the percent value of standard deviation of J(ρ) obtained by running a certain metaheuristic algorithm divided to the average value of the solution, AvgJmin:

ar=StDev%Jmin=100StDevJminAvgJmin,AvgJmin=1Nth=1NtJmin(h),StDevJmin=1Nt1h=1NtJmin(h)AvgJmin2,(40)
where Jmin is the value of J(ρ) measured after running the algorithm, and the superscript h indicates the observed trial, h=1Nt. Since Table 2 gives contains values, the values in the last column correspond actually to Avg(Jmin) instead of Jmin.

The search performance of all metaheuristic algorithms analyzed in this section is synthesized in Table 3 in terms of the values of the quality indices cs and ar.

The results presented in Table 2 indicate that the best performance of the control system as far as the average value of the cost function is concerned is obtained by SMA and SMA-F1 followed by WOA, and the worse performance is obtained by CSS and WOA-F1. The comparison of accuracy and resource usage on the basis of the results given in Table 3 reveals that the best accuracy rate is obtained by SMA and SMA-F1 followed by WOA, and the worse is obtained by GSA and GWO-F1. The best convergence speed is obtained by SMA-F1 followed by SMA-F1 and GWO, and the worse by PSO-F1 and PSO. Concluding, the best overall performance is exhibited by SMA and SMA-F1, which is followed by WOA, justifying the application of SMA in this paper.

Algorithm cs ar
SMA 19.7 0
PSO 311.4 18.67
GSA 167.2 18.86
CSS 97.9 12.07
GWO 20 18.52
WOA 18.7 6.24
SMA-F1 14.3 0
PSO-F1 269.9 11.77
GSA-F1 204.4 12.82
CSS-F1 162.9 11.33
GWO-F1 180.4 20.49
WOA-F1 11.9 14.9

CSS, Charged System Search; CSS-F1, Charged System Search with information feedback model F1; GSA, Gravitational Search Algorithm; GSA-F1, Gravitational Search Algorithm with information feedback model F1; GWO, Grey Wolf Optimizer; GWO-F1, Grey Wolf Optimizer with information feedback model F1; PSO, Particle Swarm Optimization; PSO-F1, Particle Swarm Optimization with information feedback model F1; SMA, Slime Mould Algorithm; SMA-F1, Slime Mould Algorithm with information feedback model F1; WOA, Whale Optimization Algorithm; WOA-F1, Whale Optimization Algorithm with information feedback model F1.

Table 3
Values of quality indices cs and ar.

The results presented in Table 2 show that adding the information feedback model F1 is beneficial for SMA, GSA and CSS as far as the reduction of cost function values is concerned. The results presented in Table 3 highlight that adding the information feedback model F1 offers the reduction of the convergence speed of SMA, PSO and WOA. But the main advantage of adding the information feedback model, as outlined in Table 3, is the reduction of the accuracy rate of all algorithms except GSO and WOA.

These conclusions are drawn for the process, fuzzy controller, optimization problem and dynamic regime considered in this paper. Other processes, controllers, optimization problems and dynamic regimes are expected to lead to different conclusions; challenging processes in this regard are those considered in [5764].

Two typical fuzzy control system responses are presented in Figure 4 in terms of controlled output and control signal versus time considering the average controller after the first iteration of SMA (with the parameters Be=33.2784, η=0.5268 and BΔe=0.0676), and the average controller after the application of the 20 iterations of SMA (with the parameters given in the first row in Table 2). Figure 4 clearly shows the improvement of empirical quality indices after the application of SMA. Both plots in Figure 4 also highlight the effects of the process nonlinearity in the initial and final part of fuzzy control system responses.

Figure 4

Real-time experimental results expressed as fuzzy control system responses y and u with Takagi–Sugeno–Kang proportional-integral fuzzy controllers (TSK PI-FC) after one iteration of Slime Mould Algorithm (SMA) (dotted line) and after 20 iterations of SMA (continuous line).

It is certainly obvious from Figure 4 that the output trajectories blend to reference input value as time approaches the end of the time horizon. As a consequence, with the suggested controller, the control system addressed can be stabilized. However, as pointed out in [44], the stability cannot be proved.

6. CONCLUSIONS

This paper proposed an approach to the SMA-based tuning of cost-effective fuzzy controllers for servo systems. This approach is advantageous because of two reasons, (1) and (2), which represent the pros of the approach proposed in this paper in the framework of metaheuristic optimization algorithms: (1) it is computational efficient as proved by the values of two quality indices related to accuracy and resource usage, (2) SMA has a small number of parameters. The validation and comparison on real-world servo system position control shows the superior performance of SMA versus other similar metaheuristic algorithms reflected by smaller value of the cost function, defined, as specified in Section 1, as the sum of time multiplied by squared control error, which is the discrete-time version of ITSE used in continuous-time control.

The limitations of SMA, which are also the cons of the approach proposed in this paper, are also twofold, (i) and (ii): its current implementation is offline, (ii) it requires a relatively large number of evaluations of the cost function (or the fitness function). But these two cons are also met at other metaheuristic algorithms.

Future research will be focused on the mitigation of these limitations by modifying SMA in order to work online in adaptive fuzzy control system structure, and solving other optimization problems including those specific to optimal tuning of controller parameters in optimization problems expressed as reference tracking control problems.

Moreover, except the algorithms used in this paper, some of the most representative computational intelligence algorithms can be used to solve the problems, like Monarch Butterfly Optimization (MBO) [65], Earthworm Optimization Algorithm (EWA) [66], Elephant Herding Optimization (EHO) [67] and Moth Search (MS) algorithm [68]. Another group of representative algorithms is discussed in [69].

Nevertheless, a viable way to apply the results given in this paper is represented by networked control systems. With the increasing progress in advanced data communication technology, sampled-based event-triggered control protocols have been introduced in such systems and especially in complex networked control systems to effectively alleviate the communication burden by reducing computation loads and data transmission rates. Therefore this motivates the attempt to further enhance the performance of such control systems in terms of using fuzzy control, and the optimal tuning of tuning parameters is a convenient way to cope with such systems.

CONFLICT OF INTEREST

The authors declare that they do not have any conflict of interest.

AUTHOR’S CONTRIBUTION

The contribution of the authors to this paper is equal.

ACKNOWLEDGMENTS

This work was supported by grants of the Romanian Ministry of Education and Research, CNCS - UEFISCDI, project numbers PN-III-P4-ID-PCE-2020-0269, PN-III-P1-1.1-PD-2019-0637, PNIII-P1-1.1-PD-2016-0331, PN-III-P1-1.1-TE-2019-1117, PN-III-P2-2.1-PTE-2019-0694, within PNCDI III, and by the NSERC of Canada.

REFERENCES

29.J. Pimentel, C.J. Ardila Hernandez, E. Niño, D.J. Molinares, and J. Ruiz-Rangel, Water cycle algorithm: implementation and analysis of solutions to the bi-objective travelling salesman problem, Int. J. Artif. Intell., Vol. 17, 2019, pp. 105-136.
30.M. Moattari and M.H. Moradi, Conflict monitoring optimization heuristic inspired by brain fear and conflict systems, Int. J. Artif. Intell., Vol. 18, 2020, pp. 45-62.
47.S. Preitl and R.-E. Precup, On the algorithmic design of a class of control systems based on providing the symmetry of open-loop Bode plots, Sci. Bull. UPT Trans. Autom. Control Comput. Sci., Vol. 41, 1996, pp. 47-55.
56.Modular Servo System, User’s Manual, Inteco Ltd., Krakow, Poland, 2007.
Journal
International Journal of Computational Intelligence Systems
Volume-Issue
14 - 1
Pages
1042 - 1052
Publication Date
2021/03/16
ISSN (Online)
1875-6883
ISSN (Print)
1875-6891
DOI
10.2991/ijcis.d.210309.001How to use a DOI?
Copyright
© 2021 The Authors. Published by Atlantis Press B.V.
Open Access
This is an open access article distributed under the CC BY-NC 4.0 license (http://creativecommons.org/licenses/by-nc/4.0/).

Cite this article

TY  - JOUR
AU  - Radu-Emil Precup
AU  - Radu-Codrut David
AU  - Raul-Cristian Roman
AU  - Emil M. Petriu
AU  - Alexandra-Iulia Szedlak-Stinean
PY  - 2021
DA  - 2021/03/16
TI  - Slime Mould Algorithm-Based Tuning of Cost-Effective Fuzzy Controllers for Servo Systems
JO  - International Journal of Computational Intelligence Systems
SP  - 1042
EP  - 1052
VL  - 14
IS  - 1
SN  - 1875-6883
UR  - https://doi.org/10.2991/ijcis.d.210309.001
DO  - 10.2991/ijcis.d.210309.001
ID  - Precup2021
ER  -