Open access peer-reviewed chapter

A New Methodology for Tuning PID‐Type Fuzzy Logic Controllers Scaling Factors Using Genetic Algorithm of a Discrete‐Time System

By Wafa Gritli, Hajer Gharsallaoui and Mohamed Benrejeb

Submitted: November 1st 2016Reviewed: May 5th 2017Published: August 30th 2017

DOI: 10.5772/intechopen.69604

Downloaded: 368

Abstract

In this chapter, a proportional‐integral derivative (PID)‐type fuzzy logic controller (FLC) is proposed for a discrete‐time system in order to track a desired trajectory generated using the flatness property. In order to improve the performance of the proposed controller, genetic algorithm (GA) based on minimizing the integral of the squared error (ISE) is used for tuning the input and output PID‐type FLC scaling factors online. The considered controller is applied to an electronic throttle valve (ETV). GA tuning shows a better and robust performance compared to Simulink design optimization (SDO) algorithm in terms of tracking a desired trajectory with disturbances rejection.

Keywords

  • PID‐type FLC
  • scaling factors
  • genetic algorithm
  • integral of the squared error
  • Simulink design optimization technique
  • flatness
  • electronic throttle valve

1. Introduction

Fuzzy logic control (FLC) has been widely used in many successful industrial applications. The first FLC algorithm was implemented by Mamdani in 1974. Unlike conventional control, which is based on mathematical model of a plant, an FLC usually embeds the intuition and experience of a human operator and may provide a nonlinear relationship induced by membership functions, rules and defuzzification. In that respect, FLC has been reported to be successfully used for a number of complex and nonlinear systems and are proved to be more robust and their performances are less sensitive to parametric variations than conventional controllers.

In the literature, various types such as proportional integral (PI), proportional derivative (PD) and proportional‐integral derivative (PID) of FLCs have been proposed. For example, PI‐type FLCs have been successfully implemented in many physical applications such the control of the temperature and pressure of a steam engine and control the steering and speed of an automobile. However, performance of PI‐type FLCs for higher order systems and nonlinear systems may be poor due to the large overshoot and the excessive oscillation. PD‐type FLCs are suitable for a limited class of systems and they are not recommendable in the presence of measurement noise and sudden load disturbances. Theoretically, PID‐type FLCs provide a good performance. However, there are difficulties associated with the generation of an efficient rule base and the tuning of parameters.

In the proposed PID‐type FLC, the design of parameters within two groups: structural parameters and tuning parameters. Basically, structural parameters include input/output (I/O) variables to fuzzy inference, fuzzy linguistic sets, membership functions, fuzzy rules, inference mechanism and defuzzification mechanism, which are usually determined during offline design. Tuning parameters include I/O scaling factors (SF) and parameters of membership functions (MF), which can be calculated during online adjustments of the controller in order to enhance the process performance [1].

The appropriate selection of input and output scaling factors is very important because they have significant effects on the dynamic of fuzzy controller. This leads researchers to explore the best method in searching optimum PID‐type FLC parameters. Various strategies or methods have been used up to now. In Ref. [2], Qiao and Mizumoto proposed a peak observer mechanism‐based method to adjust the PID‐type FLC parameters. This self‐tuning mechanism decreases the equivalent integral control component of the fuzzy controller gradually with the system response process time. Furthermore, Woo et al. [3] developed a method based on two empirical functions evolved with the system’s error information. In Ref. [1], Guzelkaya et al. proposed a technique that adjusts the scaling factors, corresponding to the derivative and integral components, using a fuzzy inference mechanism. However, the major disadvantages of all these PID‐type FLC tuning method are the difficult choice of their relative parameters and mechanisms. To overcome these difficulties, differential search algorithm (DSA) meta‐heuristic technique is proposed for systematically tuning the scaling factors of the PID‐type FLC in Ref. [4]. The fuzzy control design is formulated as a constrained optimization problem, which is efficiently solved based on an improved DSA. In this proposed technique, different optimization criteria such as integral square error (ISE) and maximum overshoot are considered in order to guarantee more robustness and performance control objectives.

In this chapter, a genetic algorithm (GA)‐based heuristic optimization technique has been implemented to obtain better performance compared to the Simulink design optimization (SDO) technique. GA is based upon minimizing the error between the output system and the desired trajectory starting from a flat output variable generated using the flatness property. Various performance indices can be used. In this study, the integral of squared error (ISE) index has been used in order to minimize the error between the output and the desired flat trajectory. The methods are applied in a discrete‐time framework to an electronic throttle valve as a case of study.

2. PID‐type fuzzy logic controller description

In this study, we will deal with fuzzy PID‐type controllers formed using one PD‐type FLC with an integrator at the output.

The PID‐type fuzzy logic controller structure is shown in Figure 1 [5], where Keand KdKe,KdR+are the input scaling factors are α and β α,βR+the output scaling factors.

Figure 1.

PID‐type FLC.

The inputs variables, well known as the error ek between the desired trajectory ykdand the measure yk, as well as the error variation Δek given by Eqs. (1) and (2) where Te is the sampling period.

E1
E2

The output variable Δukof such a controller is the variation of the control signal ukwhich can be defined as Eq. (3).

E3

The output of the PID‐type fuzzy is given by Eq. (4) [5]

E4

Thus, the equivalent control components of the PID‐type FLC are such that

Proportional gain: αKeP+βKdD

Integral gain: βKeP

Derivative gain: αKdD

where the terms Pand Dare given by Eqs. (5) and (6) [2].

E5
E6

The fuzzy controllers with a product‐sum inference method, centroid defuzzification method and triangular uniformly distributed membership functions for the inputs and a crisp output proposed in Refs. [2, 6] are used in our case of study.

Table 1 gives the linguistic levels, assigned to the variables ek, Δek and Δuk, as follows: NL: negative large; N: negative; ZR: zero; P: positive; PL: positive large.

ek\ΔekNZRP
NNLNZR
ZRNZRP
PZRPPL

Table 1.

Fuzzy rules‐base.

3. Scaling factors tuning using genetic algorithm

In this work, a new method is proposed for tuning the coefficients of PID‐type FLCs. This method adjusts the input scaling factor corresponding to the derivative coefficient and the output scaling factor corresponding to the integral coefficient of the PID‐type FLC using genetic algorithm, as shown in Figure 2. The integral of squared error (ISE) index has been used in order to minimize the error between the output and the desired flat trajectory.

Figure 2.

PID‐type FLC scaling factors tuning.

3.1. Genetic algorithm

Genetic algorithm was first proposed by Holland [7]. It is a heuristic optimization technique inspired by the mechanism of natural selection. It is used in order to solve highly complex problems. GA starts with an initial population containing a number of parameters, where each one is regarded as the genes of a chromosome and can be structured by a string of concatenated values. Each chromosome represents a solution of the problem and its performance is evaluated based on fitness function.

In the beginning, an initial chromosome population is randomly generated. The chromosomes are candidate solutions to the problem. Then, the fitness values of all chromosomes are evaluated by calculating the objective function. So, a group of the best chromosomes is selected based on the fitness of each individual. In this ‘surviving’ population, the genetic operators of crossover and mutation are applied in order to create the next population solution. The above steps are repeated until a specific termination criterion is found.

  • Reproduction: Create a part of the new population by simply copying without changing the selected individuals from the present population. Also, new population has the possibility of selection by already developed solutions [8].

  • Crossover: Create new individuals as offspring of two parents. It is a recombination operator that combines selected subparts called crossover points of two parent chromosomes. The individuals resulting in this way are the offspring [8].

  • Mutation: Create a new individual for the new population by randomly mutating a selected individual. The modifications can consist of changing one or more values in the representation or adding/deleting parts of the representation [8].

To compute the fitness of each chromosome, the objective functions are used. Many authors use integral of time multiplied by absolute error (ITAE), mean of the squared error (MSE), integral of absolute error (IAE) and integral of the squared error (ISE) as performance index [9, 10].

In this chapter, the method of tuning PID‐type FLC parameters using GA consists in finding the optimal I/O scaling factors, which minimize the defined objective function, chosen as the ISE in order to specify more performance in terms of tracking a desired trajectory.

If yd(t)is the desired trajectory and yis the output trajectory, then error e(t) is

E7

and the ISE can be defined by

E8

Fitness function is taken as inverse of error, i.e., performance index.

E9

3.2. Tuning procedure

The overall flowchart for optimization using GA is shown in Figure 3. Initially, a number of populations N have been generated for the scaling factors Ke,Kd,α and β. Each individual of these N sets in the current population is evaluated using the objective function ISE. Based on the values of the objective function, out of these N possible solutions, the good solutions are retained and the others are eliminated. A new population is formed by applying the genetic operators (reproduction, crossover and mutation) to these selected individuals. This process of production of a new generation and its evaluation is performed repetitively. The algorithm continues until the population converges to the stop criterion.

Figure 3.

Flowchart of the GA optimization algorithm.

4. Flatness and trajectory planning

The flat property has been introduced in Ref. [11] for continuous‐time nonlinear systems. It can be stated in a discrete‐time version which leads for the design of a control which ensures a tracking of a desired trajectory. One major property of differential flatness is that the state and the input variables can be directly expressed, without integrating any differential equation, in terms of the flat output and a finite number of its derivatives. The flatness approach will be used in this chapter in a discrete‐time framework.

The studied dynamic linear discrete system is described by Eq. (10).

E10

where q is the forward operator, ukand ykare the input and the output, respectively, and A(q)and B(q)are polynomials defined by

E11
E12

where the parameters aiand biare constants, i=0,1,…,n−1. The partial state of such a dynamic system can be considered as a discrete flat output zkwhich can be expressed as a function of input and output signals as following

E13
E14

Often, the real output signal ykto be controlled is not a flat output. Then, it is necessary to plan a desired trajectory for the flat output [11] and to consider thereafter the relation (14).

The open loop control law can be determined by the following relations [12].

E15
E16

where fand g are the vectorial functions. Then, it is sufficient to find a desired continuous flat trajectory tzd(t)that must to be differentiable at the (r+1)order.

The polynomial interpolation technique is used in order to plan the desired flat trajectory zd(t). Let consider the state vector zd(t)=(zd(t)z˙d(t)zd(r+1)(t))Tcontaining the desired continuous flat output and its successive derivatives. The expression of zd(t)can be given in Eq. (17); t0and tfare the two moments known in advance.

E17

where M1and M2are such as [12]

E18
E19

and the vectors c1and c2defined by

E20
E21

Then, the output desired trajectory ykdis defined. In the discrete‐time framework, the real output ykhas asymptotically to track this such as Eq. (22).

E22

In the following section of this chapter, the efficiency of the proposed methodology for tuning PID‐type FLC scaling factors has been validated on a discrete‐time system: an electronic throttle valve for a defined desired trajectory generated using the flatness property and compared to the Simulink design optimization (SDO) technique.

5. Case of study: electronic throttle valve (ETV)

Throttle valve is one of the most important devices in the engine management system. In conventional engine, the amount of airflow into the combustion system has been adjusted by the throttle valve, which is connected mechanically to an accelerator pedal [13]. The electronic throttle body (ETB) regulates air inflow into the car engine. Compared to the mechanical throttle, a well‐controlled ETB can reduce fuel consumption.

5.1. System modelling

The case of the ETV is described in Figure 4.

Figure 4.

Electronic throttle valve system.

The electrical part is modelled by Eq. (23)

E23

where L is the inductance R is the resistance u(t) and i(t) are the voltage and the armature current, respectively, kvis an electromotive force constant and ωm(t)is the motor rotational speed.

The mechanical part of the throttle is modelled by a gear reducer characterized by its reduction ratio γsuch as Eq. (24)

E24

where CLis the load torque and Cg is the gear torque. The mechanical part is modelled according to Eq. (25), such that [14, 15].

E25

and

E26

where θ(t)is the throttle plate angle, J is the overall moment of inertia, Ce=kei(t)is the electrical torque where Ke is a constant, Cf is the torque caused by mechanical friction, Cr is the spring resistive torque and Cais the torque generated by the airflow. The electronic throttle valve involves two complex nonlinearities due to the nonlinear spring torque Crand the friction torque Cf. They are given by their static characteristics [16]. The static characteristic of the nonlinear spring torque Cr is defined by

E27

For θminθθmax,kr, is the spring constant, Dis a constant, θ0is the default position and sgn(.) is the following signum function

E28

The friction torque function Cfof the angular velocity of the throttle plate can be expressed as

E29

where fv and fc are two constants. By substituting in Eq. (25) the expressions of Cg, Cf and Cr and by neglecting the torque generated by the airflow Ca, the two nonlinearities sgn(θθ0)and sgn(ω)and the two constants krγθ0and fv the transfer function of the linear model becomes (30) [15].

E30

with ks = (180 / π / γ2)kr and q as the Laplace operators.

5.2. Simulation results

The identified parameters of H(q)are given in Table 2 at 25°C temperature [15].

ParametersValues
R(Ω)2.8
L(H)0.0011
ke (N.m/A)0.0183
kv (v/rad/s)0.0183
J (kg.m2)4 × 10−6
γ16.95

Table 2.

Model’s parameters.

The corresponding discrete‐time transfer function is given by Eq. (31) for the sample time Te= 0.002s.

E31

The desired continuous time flat trajectory zd(t)can be computed according to the following polynomial form

E32

where cst1and cst2are constants, k0=3s,k1=6s,k2=10sandk315sare the instants of transitions, B(1)is the static gain between the flat output zkand the output signal ykfor each operating mode and poly1(k)andpoly2(k)are polynomials calculated using the technique of polynomial interpolation.

The desired trajectory is then given in Figure 5.

Figure 5.

Desired trajectory.

The obtained optimal I/O scaling factors (Ke,Kd,α,β)for Simulink design optimization technique and GA are summarized in Table 3.

SF\MethodSDOGA
Ke0.11440.0086
Kd1.59970.4612
α2.62390.1108
β0.00010.1934

Table 3.

PID‐type fuzzy scaling factors values.

The obtained results are given in Figures 69.

Figure 6.

System outputs using Simulink design optimization.

Figure 7.

System outputs using genetic algorithm.

Figure 8.

Tracking errors.

Figure 9.

Control signals.

Figures 8 and 9 show the responses with Simulink design optimization technique and GA tuning using ISE criterion. Based on a comparative analysis, better results were there obtained with the GA tuning method.

All results, for obtained scaling factors values, are acceptable and show the effectiveness of the proposed GA tuning method in terms of the tracking desired trajectory with disturbances rejection in comparison with the SDO technique.

6. Conclusion

In this chapter, an optimization technique was introduced to tune the parameters of PID‐type fuzzy logic controller (FLC). The idea is to use the genetic algorithm (GA)‐based heuristic optimization technique in order to solve highly complex problems. In order to specify more robustness and performance of the proposed GA‐tuned PID‐type FLC, optimization criteria such as integral square error (ISE) is considered.

The proposed controller is applied to an electronic throttle valve (ETV) in the discrete‐time framework in order to track a desired trajectory starting from a flat output generated using flatness property. The performance comparison with the Simulink design optimization (SDO) technique shows the efficiency of the proposed GA‐tuned approach in terms of tracking a desired trajectory with disturbances rejection.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Wafa Gritli, Hajer Gharsallaoui and Mohamed Benrejeb (August 30th 2017). A New Methodology for Tuning PID‐Type Fuzzy Logic Controllers Scaling Factors Using Genetic Algorithm of a Discrete‐Time System, Modern Fuzzy Control Systems and Its Applications, S. Ramakrishnan, IntechOpen, DOI: 10.5772/intechopen.69604. Available from:

Embed this chapter on your site Copy to clipboard

<iframe src="http://www.intechopen.com/embed/modern-fuzzy-control-systems-and-its-applications/a-new-methodology-for-tuning-pid-type-fuzzy-logic-controllers-scaling-factors-using-genetic-algorith" />

Embed this code snippet in the HTML of your website to show this chapter

chapter statistics

368total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Applications of the Fuzzy Logic to the Energy Conversion Systems on Board of UAVs

By Dinca Liviu and Corcau Jenica Ileana

Related Book

First chapter

A Real-Time Speech Enhancement Front-End for Multi-Talker Reverberated Scenarios

By Rudy Rotili, Emanuele Principi, Stefano Squartini and Francesco Piazza

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More about us