The Design of Frequency Filters of Iterative Feedback Tuning Using Particle Swarm Optimization

Iterative feedback tuning (IFT) is a data-based tuning approach that minimizes a quadratic performance index using some closedloop experimental data. A control weighting coefficient, known as lambda, and two frequency filters are the most important parameters which can significantly improve the performance of the method. One of the major problems in IFT is tuning these parameters. This paper presents a new approach to tune frequency filters using particle swarm optimization (PSO). At the end, the performance of the proposed method is evaluated by two case study simulations.


Introduction
Iterative feedback tuning (IFT) is a data-based method for the tuning of controllers with restricted complexity which was proposed by Hjalmarsson et al. [1] in 1994.In this method, the problem of model bias could be avoided by replacing the information carried by the model with information achieved directly from the system itself.This leads to an iterative method where the controller parameters are successively updated using information from closed-loop experiments with the most recent controller in the loop.IFT has proved to be very effective in practice and is now greatly used in process control [2].Since 1994 many experiences have been achieved by IFT algorithm and lots of improvements have been observed [3][4][5][6][7][8][9][10][11][12][13].
In order to fully take advantage of IFT, a new scheme needs to be introduced to tune the efficient parameters.In earlier studies, the importance of IFT parameters was confirmed and a few tuning key points were presented [8], but the tuning of these parameters has not been treated in much detail.However, some approaches have been proposed to tune similar parameters.Shridhar and Cooper [14] derived an analytical expression in which the suppression coefficients are calculated as a function of the plant model parameters.
Kai et al. [15] proposed a min-max algorithm to design tuning parameters.Al-Ghazzawi et al. [16] presented an approach to tuning MPC on-line based on sensitivity equations derived from a step response model with linear constraints.The IFT's parameters are mainly a scalar, named lambda, and two frequency filters.This paper provides a method to tune the frequency filters which can emphasize or suppress specific frequency bands of the output and control signals.These filets affect the IFT's performance greatly.Therefore, convenient design of these filters seems critical.In order to tune these filters, the particle swarm optimization (PSO) is proposed, which tunes frequency filters using the inputoutput data achieved from IFT algorithm.This method is proposed for the first order filters with two parameters.
The paper is organized as follows: a brief explanation of IFT algorithm is presented in Section 2. In Section 3, the particle swarm optimization method is briefly described and implemented to tune the frequency filters.Consequently, we shed some light on the performance of the proposed method for tuning the frequency filters, by applying the improved IFT to some case study systems in Section 4. Finally, the paper is concluded in Section 5 with a summary of key points and results.
2 Advances in Electrical Engineering

Iterative Feedback Tuning Approach
Consider the SISO closed-loop system as presented in Figure 1.From the figure, it can be followed that where  is an unknown single input-single output (SISO) operator which is linear time invariant (LTI), () is the process output, () is the corresponding process input, and V() is an immeasurable disturbance which is assumed to be stochastic.  is an external deterministic reference signal which is independent of {V()}. = {  (),   ()} are linear time invariant transfer functions parameterized by parameter vector .
Assuming   as a desired output response to the reference signal , for closed-loop system, the error between the achieved and desired response is ( The control design objective can be formulated as where [⋅] denotes expectation with respect to the disturbance V. Scalar  is the weighting coefficient expressing the relative importance of the penalty on the control effort, and the symbols   and   are some frequency filters which can emphasize or suppress specific frequency bands of the output and control signals.These filters can be used for some purposes such as (i) emphasizing or suppressing specific frequency bands of the outputs and control signals, for instance, to prevent unwanted oscillations in these signals: (ii) using as notch filters in the frequency bands where the measurement noise dominates; (iii) meeting specific frequency domain performance specifications, such as constraints on the sensitivities.
For minimizing , IFT estimates the gradient / only using some special closed-loop experiments.The detailed procedure of estimating the gradient can be seen in [8].
After estimating the gradient /, the parameters of the controller are updated by the following iterative algorithm: Here,   is an appropriate positive definite matrix which is typically a Gauss-Newton approximation of the Hessian of  while   is a real positive scalar which defines the step size of the algorithm.

Particle Swarm Optimization
Particle swarm optimization (PSO) is a population-based optimization approach first proposed in 1995 [17,18].This method is urged by the observation of social interaction and animal behaviours such as fish schooling and bird flocking [19][20][21].Similar to most optimization techniques, PSO requires a cost function evaluation function relevant to the particle's position.As   and   affect the output signal and the input signal, the proposed cost function consists of the integrated absolute error (IAE) defined as follows [22]: And the second part of the cost function is total variation (TV) of the manipulated input.This criterion is used to evaluate the required control effort which is a convenient measure of "smoothness" of control input and should be as small as possible [22]: Thus, the cost function can be defined as All solutions in PSO can be represented as particles in a swarm.Each particle has a position and velocity vector and each position coordinate represents a parameter value.The current position of th particle of the swarm is denoted by and  gbest  are the personal best (Pbest) position and global best (Gbest) position of the th particle.Each particle is initialized with a random position and velocity.The velocity of each particle is accelerated toward the global best and its own personal best based on the following equation: Here  1 and  2 are two random numbers generated from a uniform distribution in the range [0, 1];  1 and  2 are the acceleration constants and  is the inertia weight factor.Suitable selection of these parameters provides a balance between global and local explorations and helps the particles  converge to Gbest.In this paper, the constricted version of PSO [23] is used which sets the parameters as follows: where The positions are updated based on their movement over a discrete time interval (Δ) as follows: where Δ usually is set to 1. Then the cost function at each position is reevaluated to update Gbest and Pbest.In this paper, the frequency filters are considered as the first order filter as follows: Therefore, the proposed PSO is employed to tune four parameters  1 ,  1 ,  2 , and  2 .In each iteration of PSO, the IFT algorithm obtains the optimal control parameter with the current filters parameters.Using the updated controller parameters, cost function ( 7) is changed and it leads to different parameters of frequency filters achieved in PSO algorithm.

Simulation Study
Example 1.A gas boiler system is considered where the gas valve position  2 is used to control the temperature  of the water.The following transfer function is used to represent the relationship between the variations in  2 and  when the flow of water is constant [24]: This system is controlled by a PI controller tuned by IFT algorithm.The frequency filters tuned by the proposed PSO method are as follows: The simulation results in Figure 2 and Table 1 indicate that the IFT algorithm tuned by the proposed method gives an acceptable response which can substitute the trial-and-error tuning of the filters.As mentioned before, the frequency filters play a significant role in IFT performance.By many simulations it was observed that IFT is very sensitive to filters parameters.By a random selection of the filters parameters, it was observed that poorly tuned filters can lead to unstable responses.To show the filters effect on the responses, the filters parameters are obtained randomly around the optimal values by a 30% tolerance.From Figure 2 and Table 1, it can be seen that small variations of filters parameters can lead to large changes in output response.Example 2. Consider the control system of a heated tank.The dynamics relating the control action and the temperature can be modelled with an FOPDT model as follows [24]: The output response and control signal of the system are depicted in Figure 3 in which the IFT tuned by the proposed method has been compared to IFT with the randomly tuned filters.Here, more freedom was assigned to select the random filters.Therefore, the effect of these filters has been shown better.The simulation results can be seen in Table 2.

Conclusion
The earlier attempts to tune the frequency filters were based on trial-and-error.In this study, these filters were designed using the PSO method.The proposed method was applied to a case study system.It was observed that the performance of IFT modified by the proposed methods was greatly improved.Therefore, IFT has become more capable with a better performance.The proposed method for designing frequency filters in this study is not merely limited to IFT algorithm; therefore, future studies are recommended to generalize the proposed approach to other methods with similar parameters.Also, a further study with focus on high order filters is suggested.

Figure 2 :
Figure 2: Simulation results of Example 1 for evaluating the effect of proposed filters.

Figure 3 :
Figure 3: Simulation results of Example 2 for evaluating the effect of proposed filters.

Table 1 :
Controller setting for Example 1 (evaluating the proposed frequency filters effect).

Table 2 :
Controller setting for Example 2 (evaluating the proposed frequency filters effect).