Next Article in Journal
Comparison of Different Numerical Interface Capturing Methods for the Simulation of Faraday Waves
Previous Article in Journal
Degradation of Landfill Leachate Using UV-TiO2 Photocatalysis Combination with Aged Waste Reactors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimization of a Fuzzy Automatic Voltage Controller Using Real-Time Recurrent Learning

1
Facultad de Ingeniería, Universidad Distrital Francisco José de Caldas, Bogotá 11021-110231588, Colombia
2
Departamento de Ingeniería Eléctrica, Electrónica de Computadores y Sistemas, Campus de Viesques, Universidad de Oviedo, 33204 Gijón/Xixón, Spain
*
Author to whom correspondence should be addressed.
Submission received: 26 April 2021 / Revised: 17 May 2021 / Accepted: 21 May 2021 / Published: 27 May 2021
(This article belongs to the Section Automation Control Systems)

Abstract

:
The automatic voltage regulator is an important component in energy generation systems; therefore, the tuning of this system is a fundamental aspect for the suitable energy conversion. This article shows the optimization of a fuzzy automatic voltage controller for a generation system using real-time recurrent learning, which is a technique conventionally used for the training of recurrent neural networks. The controller used consists of a compact fuzzy system based on Boolean relations, designed having equivalences with PI, PD, PID, and second order controllers. For algorithm implementation, the training equations are deduced considering the structure of the second order compact fuzzy controller. The results show that a closed-loop fuzzy control strategy was successfully implemented using real-time recurrent learning. In order to implement the controllers optimization, different weighting values for error and control action are used. The results show the behavior of the configurations used and its performance considering the steady state error, overshoot, and settling time.

1. Introduction

The power grid is a non-linear complex system consisting of various interconnected systems or control areas. In general, the distribution network is employed to control various nodes in the supply network with power supply stabilizers, power systems, and Automatic Voltage Regulators (AVRs) [1,2].
It is also noteworthy that the quick expansion of energetic resources distributed along the low voltage network induces violations in the voltage limits, especially in rural zones with a low electrical charge. Therefore, the use of devices for line voltage regulation arises as an approach to solve such issue and at low costs for the grid [3].
Following [4], excitation systems for synchronous generators contribute to the effective control of voltage and the stability of the energy system. Even if such an excitation system is adjusted to allow the proper operation of the electric energy scheme in a wide range of operating conditions, it may be necessary to tune it again when aiming at the improvement of the system’s stability amid disregarded operation conditions.
In a generator, the excitation system helps keep the energy and the control of power flux employing an automatic voltage regulator. The job of an AVR consists of maintaining the magnitude in the energy of a synchronous generator at a certain level; hence, the stability in an AVR may affect the security in an energy system [5]. Furthermore, an AVR can be found in a Direct Current (DC) Microgrid system to regulate the Alternating Current (AC) power produced by the generators before rectification and then supplied to DC grid [6].

1.1. Automatic Voltage Regulator Design

Regarding related works to a voltage controller design, AVR voltage responses and frequencies for different gains in a specified limit for output voltage are studied in [7]. This work proposes a Proportional Integral (PI) controller achieving a response with minimal overshoot and small settling time.
Reference [2] displays a study of predictive Proportional Integral Derivative (PID) automatic voltage regulator used in a multi-machine feeding system. This work carries out the comparison with a traditional PID controller observing that the proposed method may provide better attenuation of the system’s oscillators in the power grid with small and large disturbances.
Meanwhile, an automatic voltage regulator AVR based on series voltage compensation with an AC chopper is proposed in [8], this consists of an AC chopper with Pulse Width Modulation (PWM), and a transformer for the series voltage compensation. The AC chopper provides AC-AC direct energy conversion with no energy storage elements; thus, the size and cost of the AVR are reduced. The proposed AVR can compensate the voltage drops and the increment of voltage in the power input.

1.2. Automatic Voltage Regulator Optimization

Given the relevance of the AVR systems in the generation of energy, the literature allows observing different works performing the optimization of controllers employed with the AVRs. In [5], is displayed an approach with evolutionary computation to establish the optimal values of the gains of a PID employing evolutionary algorithms. This work shows the relevance of achieving the controller’s adjustment.
Regarding the employment of other bio-inspired algorithms, in [9] the optimization of a Fractional-Order Proportional-Integral-Derivative (FOPID) controller is performed, since such controller includes two additional parameters compared to a conventional PID and its adjustment process is more complex. Therefore, it is used the Salp Swarm Optimization Algorithm (SSA) for the selection of the optimized parameters of the controller FOPID to achieve the optimal dynamic response and stability; besides, an analysis of stability is shown utilizing zero pole stability criteria and bode diagrams.
A method for the optimal adjustment of a FOPID to an automatic voltage regulator is presented in [10]. The method is based on the Yellow Saddle Goatfish Algorithm (YSGA). The performance of the obtained FOPID controller is verified when compared to various FOPID controllers adjusted by other metaheuristic algorithms.
In addition, in [11] is presented a technique to determine optimal value gains of a PID controller for an AVR using Cuckoo Search (CS) evolutionary algorithm. The dynamic performance of the proposed controller is evaluated considering the transitory response characteristics like the rise time, settling time, overshoot, and steady state error.
Finally, in [12] is displayed a hybrid metaheuristic method for the optimal adjustment of four different types of PID controllers for an automatic voltage regulator system. The method is based on the optimization algorithm of Simulated Annealing - Manta Ray Foraging Optimization (SA-MRFO). The performance of the ideal PID, real PID, fractional order PID, and PID obtained from second-order derived controllers is verified through the comparison with the tuned controllers by different algorithms presented in the literature.

1.3. Real-Time Recurrent Learning

The Real-Time Recurrent Learning (RTRL) algorithm is usually employed for identification and control of dynamic systems using recurrent neural networks [13,14]. This algorithm corresponds to the application of the chain rule when using the descending gradient method for the adjustment of parameters in a recurrent neural network. Even though this is mostly used for training of neural networks, it can be also applied in other optimization applications. This algorithm calculates the derivatives of states and outputs concerning all weights as the network processes the sequence (during the forward step) by which is not necessary the neural network unfolding [13]. The derivative of the states regarding the weights at the moment n is calculated starting from the states, and the derivatives at the moment n 1 , and the input at the moment n. In this way, a set of recurrence equations is employed for the neural network training.
In this paper, is used the RTRL algorithm for the optimization of a fuzzy controller in close-loop, however, to present the RTRL theory is considered the classical application for a Recurrent Network (RN). An example of a fully connected RN is shown in Figure 1, where z 1 represents a memory element (delay).
Considering the presented in [15], for a step time n, the parameters of a fully connected recurrent network are:
  • x k [ n ] : signal applied to input node k.
  • y k [ n ] : output of processing node k.
  • d k [ n ] : target output of processing node k.
  • s k [ n ] : activation of processing node k.
  • w i j : weight values connecting node j to node i.
  • I: set of indexes associated with the input nodes.
  • U: set of indexes associated with the processing nodes.
The neural connections for the output of a processing node k in the step time n + 1 is y k [ n + 1 ] = f k ( s k [ n ] ) , which is calculated using Equation (1).
y k [ n + 1 ] = f k p I w k p x p [ n ] + q U w k q y q [ n ] k U
In Equation (1) f k is the activation function. On the other hand, taking T [ n ] as the set of indexes belonging to U representing the existence of a target value for a processing node k in a step time n, then, the total instantaneous squared error at the step time n is given by Equation (2).
E [ n ] = 1 2 k T [ n ] d k [ n ] y k [ n ] 2
In Equation (2), the error signal for each node is e k [ n ] = d k [ n ] y k [ n ] . Using the gradient descent method where η is the learning rate, the weights can be updated using Equation (3).
w i j [ n ] = w i j [ n 1 ] η d E [ n ] d w i j
In this way is observed that:
d E [ n ] d w i j = k T [ n ] e k [ n ] d y k [ n ] d w i j = k T [ n ] e k [ n ] p i j k [ n ] i U , j U I
Calculating the derivative of y k [ n + 1 ] (given in Equation (1)) with respect to w i j is obtained:
p i j k [ n + 1 ] = f k ( s k [ n ] ) z j [ n ] δ k i + q U w k q p i j q [ n ] i U , k U , j U I
In Equation (5) δ k i is the Kronecker delta and z j [ n ] is defined as presented in Equation (6).
z j [ n ] = x j [ n ] i f j I , y j [ n ] i f j U .

1.4. Proposal Approach and Document Organization

This article shows the optimization of a fuzzy automatic voltage controller for a generation system employing the algorithm real-time recurrent learning. The controller’s optimization is performed in discrete time domain. The architecture of this controller consists of a fuzzy compact structure based on Boolean relations. The main purpose of the document is to deduce the equations of the algorithm for the controller’s training, as well as its application for the AVR system using the weighting values for error and action control obtaining different behaviors. In this way, for the application of generation systems the real-time recurrent learning is a suitable option for training neuro-fuzzy controllers, which is the central aspect of this article.
The main contributions of this work are described as follows.
  • Considering a close-loop control system, it is presented the deduction of the training equations employed to implement the training algorithm using the AVR discrete time model. Furthermore, the algorithm steps used for controller optimization are presented.
  • Fuzzy controller schemes are proposed making analogies with linear controllers and the result of the training process is shown.
The document is organized as follows. Section 2 details the model of an automatic voltage regulator. Section 3 displays the architecture of the fuzzy controller (compact fuzzy based on Bolean relations). Section 4 presents the deduction of training equations and the algorithm used for parameter adaptation. Next, Section 5 presents the implementation of compact fuzzy controllers: Second Order (SO), PI (Proportional Integral), Proportional Derivative (PD), and Proportional Integral Derivative (PID). Section 6 describes the experimental results; finally, Section 7 presents the conclusions.

2. Automatic Voltage Regulator

In energy systems, the automatic voltage regulator is a crucial part for achieving energy exchange and regulating the energy for the consumers. The AVR in the power system allows supplying quality reliable energy aiming at reducing the steady-state error to zero in an interconnected power system [7]. In the design of an AVR, the main requisites are quick response, low overshoot and steady-state error to zero at deviation from reference voltage [2]; all of this is to keep the reliability in the operations of the power system. Generally, a real power station has more than one generator connected to a busbar and each one has an AVR [16,17].
In Figure 2, is observed the schematic diagram of an AVR system. A simple AVR includes four main components: an amplifier, an exciter, a generator, and a sensor [18,19]. In this order, the controller allows having a suitable behavior of AVR.
As mentioned, the purpose of designing an AVR is to keep the output voltage at a certain level in a generation system. An electric power grid includes more than one generator connected to a similar busbar and each one has an AVR. Since the purpose of this system is to control the energy of the power grid to which the generator is connected by a power transformer, the level of energy is constantly measured as a feedback signal using an energy sensor that may consist of a transformer. After being rectified and filtered, this signal is compared to the reference voltage in the comparator to obtain a voltage error signal. The error signal passes through the controller, then it is amplified to feed the exciter for the adjustment of the voltage/current of the field winding in the generator’s in a way that any voltage deviation may be compensated [20].
According to [5], to obtain the mathematical model and transfer function of the AVR, the four main components are linearized taking into account their time constants. Considering the presented in [2,5,21], an AVR can have a representation like the one in Figure 3, where:
  • k A : Amplifier gain (10–40), nominal 10.
  • k E : Excitation gain (1–10), nominal 1.
  • k G : Generator gain (0.7–1.0), depends on the load.
  • k S : Sensor gain with nominal value of 1.
  • τ A : Amplifier time constant (0.02 s–0.1 s), nominal 0.1 .
  • τ E : Exciter time constant (0.4 s–1.0 s), nominal 0.4 .
  • τ G : Generator time constant (1.0 s–2.0 s), depends on the load.
  • τ S : Sensor time constant (0.001 s–0.06 s), nominal 0.01 .

AVR Discrete Time Model

Considering that the fuzzy controller optimization takes place in discrete time, it is proceeding to obtain the respective AVR model in this domain. For this, the transfer function of the AVR corresponds to:
G ( s ) = k A 1 + τ A s · k E 1 + τ E s · k G 1 + τ E s · k S 1 + τ S s
G ( s ) = 10 0.1 s + 1 · 1 0.4 s + 1 · 0.8 1.5 s + 1 · 1 0.01 s + 1
G ( s ) = 8 0.0006 s 4 + 0.0679 s 3 + 0.81 s 2 + 2.01 s + 1
Transforming this transfer function into a discrete time representation considering a sampling time of T s = 0.05 s , and using the bilinear transformation method is obtained that:
G ( z ) = K 1.1 z 4 + 4.4 z 3 + 6.6 z 2 + 4.4 z + 1.1 z 4 2.021 z 3 + 0.9133 z 2 + 0.3293 z 0.2195
where K = 1 × 10 3 . In general, this transfer function can be written as:
G ( z ) = B 0 z 4 + B 1 z 3 + B 2 z 2 + B 3 z + B 4 z 4 + A 1 z 3 + A 2 z 2 + A 3 z + A 4
Using negative exponents:
Y ( z ) U ( z ) = B 0 + B 1 z 1 + B 2 z 2 + B 3 z 3 + B 4 z 4 1 + A 1 z 1 + A 2 z 2 + A 3 z 3 + A 4 z 4

3. Fuzzy Controller

According to [22], fuzzy logic has wide applicability in control systems, due to its flexibility for the implementation of control strategies, using knowledge of the system and linguistically proposing control strategies.
The architecture of the fuzzy controller is proposed considering its equivalence with a discrete time linear controller, in this way, the controller is optimized using Real-Time Recurrent Learning.
The controller architecture uses the structure of fuzzy logic systems based on Boolean relationships. In this regard, reference [23] presents the design methodology starting from a controller using Boolean sets which then become fuzzy. For the implementation of the fuzzy controller, two possible regions corresponding to negative and positive values in each universe of discourse are considered, therefore, for its implementation, membership functions of the sigmoidal type are used.
The fuzzy controller uses the compact structure presented in [24] that allows having an equivalence with a linear controller (in discrete time) using fuzzy sets to model the actions of direct inputs and the controller feedback. The general scheme of the control system can be seen in Figure 4.
To carry out the design of the fuzzy controller, a second order linear controller (in discrete time) is considered where the transfer function is:
C ( z ) = U ( z ) E ( z ) = b 0 + b 1 z 1 + b 2 z 2 1 + a 1 z 1 + a 2 z 2
The difference Equation (in discrete time) of this controller corresponds to:
u [ n ] = a 1 u [ n 1 ] a 2 u [ n 2 ] + b 0 e [ n ] + b 1 e [ n 1 ] + b 2 e [ n 2 ]
where the respective coefficients a i , b j are constant. For the fuzzy controller, these constants are replaced by non-linear relations given by fuzzy membership functions, such that:
u [ n ] = f 1 ( u [ n 1 ] ) + f 2 ( u [ n 2 ] ) + f 3 ( e [ n ] ) + f 4 ( e [ n 1 ] ) + f 5 ( e [ n 2 ] )
In general terms:
u [ n ] = f c ( u [ n 1 ] , u [ n 2 ] , e [ n ] , e [ n 1 ] , e [ n 2 ] )
For the fuzzy control system is used the membership functions shown in Figure 5, where fuzzy sets of sigmoidal type are employed to model the negative and positive values of the universe of discourse. In this way, the membership function μ i j ( x i ) is given by Equation (17).
μ i j ( x i ) = 1 + e σ i j ( x i γ i j ) 1
Considering the fuzzy sets of Figure 5, and the general structure of the controller given by the Equation (15) it is obtained the scheme of Figure 6, where the fuzzy controller used is shown.
Considering that x i { u [ n 1 ] , u [ n 2 ] , e [ n ] , e [ n 1 ] , e [ n 2 ] , } , the controller output can be calculated as:
u [ n ] = i = 1 5 j = 1 2 v i j μ i j ( x i ) = i = 1 5 f i
In this way, it is obtained Equation (19), being the set of controller parameters h i j { v i j , σ i j , γ i j } , where v i j is a scalar value and σ i j , γ i j the parameters of the membership function μ i j ( x i ) .
f i = v i 1 μ i 1 ( x i ) + v i 2 μ i 2 ( x i )

4. Fuzzy Controller Training Process

For a practical approach from the scheme of Figure 3, it is obtained the respective equivalent system in discrete time shown in Figure 4. In this model, the sensor is included as part of the plant, this being the variable to control.
In order to perform the adaptation of the fuzzy controller parameters h i j , it is used the learning rate η , if chosen with a very large value, it can cause that the algorithm does not converge, while choosing a very small value can cause the algorithm to require more time to converge [25]. Considering the presented in [25], the training method is based on the direction of maximum descent given by the gradient, therefore, the adaptation (optimization) of h i j is carried out in the following way:
h i j ( k + 1 ) = h i j ( k ) η d J d h i j
where J corresponds to the adjustment function defined by the error signal e [ n ] = r [ n ] y [ n ] then:
J = 1 2 P r [ n ] y [ n ] 2 + Q ( u [ n ] ) 2 = 1 2 P e 2 [ n ] + Q u 2 [ n ]
In Equation (21) the factor P weights the error signal and Q the control action, in a way that combinations of these values allow having different behaviors as desired.
The derivative of J depending on the adjustment parameters is:
d J d h i j [ n ] = d J d y [ n + 1 ] d y [ n + 1 ] d h i j [ n ]
On the other hand, from the transfer function (12), the corresponding recurrence equation for the plant is:
y [ n ] = A 1 y [ n 1 ] A 2 y [ n 2 ] A 3 y [ n 3 ] A 4 y [ n 4 ] + B 0 u [ n ] + B 1 u [ n 1 ] + B 2 u [ n 2 ] + B 3 u [ n 3 ] + B 4 u [ n 4 ]
Considering the controller architecture (Figure 6) and Equation (19) for f i the dynamics of the controller is given by:
u [ n ] = f 1 ( u [ n 1 ] ) + f 2 ( u [ n 2 ] ) + f 3 ( e [ n ] ) + f 4 ( e [ n 1 ] ) + f 5 ( e [ n 2 ] )
To implement the training Equation (20), the derivative of the plant output y [ n ] with respect to a controller parameter h i j is:
d y d h i j [ n ] = A 1 d y d h i j [ n 1 ] A 2 d y d h i j [ n 2 ] A 3 d y d h i j [ n 3 ] A 4 d y d h i j [ n 4 ] + B 0 d u d h i j [ n ] + B 1 d u d h i j [ n 1 ] + B 2 d u d h i j [ n 2 ] + B 3 d u d h i j [ n 3 ] + B 4 d u d h i j [ n 4 ]
In the same way, the derivative of the control action u [ n ] with respect to h i j is:
d u [ n ] d h i j = d f 1 ( u [ n 1 ] ) d h i j + d f 2 ( u [ n 2 ] ) d h i j + d f 3 ( e [ n ] ) d h i j + d f 4 ( e [ n 1 ] ) d h i j + d f 5 ( e [ n 2 ] ) d h i j
Meanwhile, the derivative of the error e [ n ] with respect to a controller parameter h i j is:
d e [ n ] d h i j = d r [ n ] d h i j d y [ n ] d h i j = d y [ n ] d h i j
In order to establish the respective derivatives using Equation (19) is obtained:
d f l ( x l ) d h i j = d d h i j ( v l 1 μ l 1 ( x l ) ) + d d h i j ( v l 2 μ l 2 ( x l ) )
where l = 1 , , 5 , i = 1 , , 5 and j = 1 , 2 , therefore, the following cases are presented:
Case 1: When l i and j can be equal to any value m = 1 , 2 , then:
d f l ( x l ) d h i j = d d x l ( v l 1 μ l 1 ( x l ) ) d x l d h i j + d d x l ( v l 2 μ l 2 ( x l ) ) d x l d h i j
d f l ( x l ) d h i j = d d x l ( v l 1 μ l 1 ( x l ) ) + d d x l ( v l 2 μ l 2 ( x l ) ) d x l d h i j
The respective derivatives are:
d d x l ( v l m μ l m ( x l ) ) = d d x l v l m 1 + e σ l m ( x l γ l m ) 1
d d x l ( v l m μ l m ( x l ) ) = v l m 1 + e σ l m ( x l γ l m ) 2 e σ l m ( x l γ l m ) σ l m
Case 2: When l = i and j = 1 , it is established that:
d f i ( x i ) d h i 1 = d d h i 1 ( v i 1 μ i 1 ( x i ) ) + d d x i ( v i 2 μ i 2 ( x i ) ) d x i d h i 1
Case 3: When l = i and j = 2 is obtained:
d f i ( x i ) d h i 2 = d d x i ( v i 1 μ i 1 ( x i ) ) d x i d h i 2 + d d h i 2 ( v i 2 μ i 2 ( x i ) )
To calculate the respective derivatives of cases 2 and 3, it is first obtained:
d d x i ( v i j μ i j ( x i ) ) = v i j 1 + e σ i j ( x i γ i j ) 2 e σ i j ( x i γ i j ) σ i j
For the other derivatives:
d d h i j ( v i j μ i j ( x i ) ) = d d h i j v i j 1 + e σ i j ( x i γ i j ) 1
In the case of the parameter h i j = v i j is determined:
d d v i j ( v i j μ i j ( x i ) ) = 1 + e σ i j ( x i γ i j ) 1 v i j 1 + e σ i j ( x i γ i j ) 2 e σ i j ( x i γ i j ) ( σ i j ) d x i d v i j
Taking the parameter h i j = σ i j it is established that:
d d σ i j ( v i j μ i j ( x i ) ) = v i j 1 + e σ i j ( x i γ i j ) 2 e σ i j ( x i γ i j ) γ i j x i σ i j d x i d σ i j
With the parameter h i j = γ i j the respective derivative is:
d d γ i j ( v i j μ i j ( x i ) ) = v i j 1 + e σ i j ( x i γ i j ) 2 e σ i j ( x i γ i j ) σ i j σ i j d x i d γ i j
For Equations (37)–(39) is considered that:
P v i j = 1 + e σ i j ( x i γ i j ) 1 P σ i j = v i j 1 + e σ i j ( x i γ i j ) 2 e σ i j ( x i γ i j ) γ i j x i P γ i j = v i j 1 + e σ i j ( x i γ i j ) 2 e σ i j ( x i γ i j ) σ i j
Furthermore, also:
R v i j = v i j 1 + e σ i j ( x i γ i j ) 2 e σ i j ( x i γ i j ) σ i j R σ i j = v i j 1 + e σ i j ( x i γ i j ) 2 e σ i j ( x i γ i j ) σ i j R γ i j = v i j 1 + e σ i j ( x i γ i j ) 2 e σ i j ( x i γ i j ) σ i j
Considering the above, Equations (37)–(39) can be written as:
d d v i j ( v i j μ i j ( x i ) ) = P v i j + R v i j d x i d v i j
d d σ i j ( v i j μ i j ( x i ) ) = P σ i j + R σ i j d x i d σ i j
d d γ i j ( v i j μ i j ( x i ) ) = P γ i j + R γ i j d x i d γ i j
Taking in a general way a controller parameter h i j it has:
d d h i j ( v i j μ i j ( x i ) ) = P h i j + R h i j d x i d h i j
In this way Equation (33) corresponding to the case when l = i and j = 1 can be written as:
d f i ( x i ) d h i 1 = P h i 1 + R h i 1 d x i d h i 1 + d d x i ( v i 2 μ i 2 ( x i ) ) d x i d h i 1
d f i ( x i ) d h i 1 = P h i 1 + R h i 1 + d d x i ( v i 2 μ i 2 ( x i ) ) d x i d h i 1
Meanwhile, for Equation (34) when l = i and j = 2 is obtained:
d f i ( x i ) d h i 2 = d d x i ( v l 1 μ i 1 ( x i ) ) d x i d h i 2 + P h i 2 + R h i 2 d x i d h i 2
d f i ( x i ) d h i 2 = P h i 2 + d d x i ( v i 1 μ i 1 ( x i ) ) + R h i 2 d x i d h i 2
Using the previous equations, the update of the controller parameters is carried out using the general Equation (20) having for each parameter the following equations:
v i j [ k + 1 ] = v i j [ k ] η P e [ n ] d e [ n ] d v i j + Q u [ n ] d u [ n ] d v i j
σ i j [ k + 1 ] = σ i j [ k ] η P e [ n ] d e [ n ] d σ i j + Q u [ n ] d u [ n ] d σ i j
γ i j [ k + 1 ] = γ i j [ k ] η P e [ n ] d e [ n ] d γ i j + Q u [ n ] d u [ n ] d γ i j

Training Algorithm

Algorithm 1 shows the steps used for controller training, in the first step is chosen the initial parameter settings for the controller. In the next step, the control system output is calculated using the plant model. Then it proceeds to perform the adjustment of the fuzzy controller parameters using the equations that involve the dynamics of the control system and the derivatives (Equations (25)–(27) and subsequent equations). It is important to note that the adjusted parameters are stored in an auxiliary variable since during this step the controller does not use these values. Then, it is returned to the step where the control system is evaluated with n = n + 1 until complete the simulation time N T . When completing the simulation time, the controller parameters are updated with the optimized values and it is returned to the step where is calculated the control system response for a new iteration k = k + 1 , until the objective function J T ( k ) given by Equation (53) is less than a defined ε , or until k is equal to a defined number K T .
It should be noted that in this algorithm the AVR discrete time model is used for controller training, therefore, it is not necessary to use additional AVR data. The controller and the plant model are integrated in a control close-loop, in this way, the control system input is the reference desired value r [ n ] .
J T ( k ) = n = 1 N T J [ k , n ]
Algorithm 1: Training controller.
Processes 09 00947 i001

5. Implementation of Compact Fuzzy Controllers

The training method is based on the direction of maximum descent given by the gradient, its convergence to the optimum point depends on the initial values of the parameters. Considering that the membership functions used allow representing the actions for positive and negative values, their initial configuration is taken from Figure 5. With this approach, it can be seen that fuzzy systems allow establishing an initial configuration that can be later optimized.
For the design of compact fuzzy controllers, linear control strategies are considered as a reference:
  • Second Order Controller (SO).
  • Proportional Integral Controller (PI).
  • Proportional Derivative Controller (PD).
  • Proportional Integral Derivative Controller (PID).
Considering the presented in [24], the compact fuzzy controllers are implemented considering Equation (24) associated with Figure 6, having the following expression:
u [ n ] = i = 1 5 f i ( x i )
where x i { u [ n 1 ] , u [ n 2 ] , e [ n ] , e [ n 1 ] , e [ n 2 ] } . In this way, to establish the equivalence with the discrete time equation of the linear controllers SO, PI, PD, PID, the functions f i with i = 1 , , 5 must be configured as presented below.

5.1. Second Order Controller Structure

The discrete linear controller considered consists of a second order system with transfer function:
U ( z ) E ( z ) = b 0 + b 1 z 1 + b 2 z 2 1 + a 1 z 1 + a 2 z 2
The respective difference equation is:
u [ n ] = a 1 u [ n 1 ] a 2 u [ n 2 ] + b 0 e [ n ] + b 1 e [ n 1 ] + b 2 e [ n 2 ]
Then for the fuzzy controller is considered the structure:
u [ n ] = f 1 ( u [ n 1 ] ) + f 2 ( u [ n 2 ] ) + f 3 ( e [ n ] ) + f 4 ( e [ n 1 ] ) + f 5 ( e [ n 2 ] )

5.2. PI Equivalent Structure

In this case, the controller has the form:
C ( z ) = k p + K i 1 z 1 = ( K p + K i ) K p z 1 1 z 1
considering the compact structure it can be expressed as:
U ( z ) E ( z ) = b 0 b 1 z 1 1 z 1
In this way, the difference equation of the PI controller corresponds to:
u [ n ] = u [ n 1 ] + b 0 e [ n ] b 1 e [ n 1 ]
Then for the PI fuzzy controller the structure is:
u [ n ] = f 1 ( u [ n 1 ] ) + f 3 ( e [ n ] ) + f 4 ( e [ n 1 ] )

5.3. PD Equivalent Structure

The PD controller can be represented as:
C ( z ) = K p + K d ( 1 z 1 ) = ( K p + K d ) K d z 1
then the expression for this is:
U ( z ) E ( z ) = b 0 b 1 z 1
The corresponding difference equation is:
u [ n ] = b 0 e [ n ] b 1 e [ n 1 ]
Then the structure for the PD fuzzy controller is:
u [ n ] = f 3 ( e [ n ] ) + f 4 ( e [ n 1 ] )

5.4. PID Equivalent Structure

In the case of the PID controller, it can be described as:
C ( z ) = K p + K i 1 z 1 + K d ( 1 z 1 )
performing the respective operations it is obtained:
C ( z ) = ( K p + K i + K d ) ( K p + 2 K d ) z 1 + K d z 2 1 z 1
in general, it can be written as:
U ( z ) E ( z ) = b 0 b 1 z 1 + b 2 z 2 1 z 1
The respective difference equation for this controller is:
u [ n ] = u [ n 1 ] + b 0 e [ n ] b 1 e [ n 1 ] + b 2 e [ n 2 ]
Then for the PID fuzzy controller the structure is:
u [ n ] = f 1 ( u [ n 1 ] ) + f 3 ( e [ n ] ) + f 4 ( e [ n 1 ] ) + f 5 ( e [ n 2 ] )

5.5. Configurations

Considering the general expression of the compact fuzzy controller given by Equation (54), the equivalences with PI, PD, PID, and second order linear controllers are obtained as shown in Table 1, this according to the respective discrete time Equations (57), (61), (65), and (70).

6. Experimental Results

In order to carry out the experimental tests, it is possible to have different combinations of P and Q, therefore, for the values selection one of these is left fixed, in this case P = 1 and Q is varied on a logarithmic scale in a way that in the results the effect of these combinations is observed. Furthermore, the learning rate is set to 2 × 10 5 . The values considered of P and Q are the following:
  • Case PQ1: P = 1 , Q = 0.01 .
  • Case PQ2: P = 1 , Q = 0.1 .
  • Case PQ3: P = 1 , Q = 1 .
  • Case PQ4: P = 1 , Q = 10 .
The performance index used corresponds to the one shown in Equation (71), where r is the desired output value, y the response of the control system, u the control signal, P and Q weighting values of the error and control action, finally, N T is the total number of data (total simulation steps).
J T = 1 2 n = 1 N T P ( r [ n ] y [ n ] ) 2 + Q ( u [ n ] ) 2
The process carried out to obtain the experimental results considers the different values of Q and the controllers, thus the steps are:
(1)
For each value of Q: 0.01 , 0.1 , 1, and 10, do the following:
(2)
For each controller: SO, PI, PD, and PID, do the following:
(3)
Perform the controller training using Algorithm 1.
(4)
Calculate the respective performance index given by Equation (71).
(5)
Store the results.
(6)
Return to step 2 until completing all controllers.
(7)
Return to step 1 until completing all values of Q.
In this way are obtained the results for the values of Q using different compact fuzzy controllers. It should be noted that since the structure and the initial configuration of the controllers given by each μ i j ( x i ) is established, only one run of the algorithm is required since the initial configuration is not random. In addition, the reference value r [ n ] = 1 since the controller seeks to regulate the normalized output voltage.
The optimization process is performed for the different cases of P and Q for each controller obtaining the results of the objective function shown in Table 2. In this table is seen that when increasing Q the value of the objective function also increases.
Regarding the performance of the control system, Table 3 shows the performance metrics corresponding to steady state error % e s s , overshoot % P O , and settling time T s . It is to be appreciated that the results with the suitable error in steady state are obtained for the combination of P and Q used in the configuration PQ3. Figure 7, Figure 8 and Figure 9 show the simulations results for the output y ( t ) , control signal u ( t ) , and the objective function J T ( k ) .
The results of the optimization process for y ( t ) are presented in Figure 7 showing the system’s output with each of the controllers after being optimized. Figure 8 describes the results for the control action u ( t ) observing that when the value of Q increases, the control action decrease. Finally, in Figure 9 is shown the evolution of the training process presenting the values of J T ( k ) for each controller considered and the respective setting employed for P and Q.
Considering case PQ3, Figure 10 displays the comparison results of controllers SO, PI, PD, and PID presenting the output signal y ( t ) and the control action u ( t ) . In these results, SO has the largest overshot and the PI controller a suitable behavior considering the overshot and setting time.
Finally, Figure 11 shows the response detail of the optimized fuzzy controllers, where the steady-state error and the overshoot are better observed for the SO, PI, PD, and PID controllers.

7. Conclusions

In this article, the process for the optimization of an automatic voltage controller using Real-Time Recurrent Learning was presented, making the deduction of the respective equations used for the implementation of this algorithm.
In order to carry out the experimental tests, different configurations for the P and Q values are considered; P weights the error signal and Q the control action, in a way that combinations of these values allow having different behaviors.
The results show the response of the different fuzzy controllers considered which were determined by the analogy of linear controllers, namely, base on preliminary knowledge of the controllers, therefore the controller topology considered was SO, PI, PD, PID.
The configurations display different dynamic behaviors according to values of P and Q. Checking the results, the control action decreases when the value of Q increases.
The record of the objective function J T ( k ) in each of the cases considered shows that the optimization process is successful using the deduced equations.
This scheme can be used to tune the controller in real time with an adaptive scheme when there are variations in the plant. Besides, in an additional work other fuzzy controller configurations could be used.

Author Contributions

Conceptualization, H.E., I.M. and H.L.; Methodology, H.E., I.M. and H.L.; Project administration, H.L.; Supervision, H.L.; Validation, I.M.; Writing—original draft, H.E.; Writing—review and editing, H.E., I.M. and H.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

In this work were not carried out test on individuals (humans).

Informed Consent Statement

The study presented in this work does not involve human beings.

Data Availability Statement

No external data were required for this work.

Acknowledgments

The authors express gratitude to the Universidad de Oviedo, and the Universidad Distrital Francisco José de Caldas.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mohamed, T.H.; Mohammed, A.; Abd-Eltawwab, A.; Hiyama, T. Wide-area Power System Oscillation Damping using Model Predictive Control Technique. IEEJ Trans. Power Energy 2011, 131, 536–541. [Google Scholar] [CrossRef]
  2. Liu, Q.; Mohamed, T.; Kerdphol, T.; Mitani, Y. PID-MPC Based Automatic Voltage Regultaor Design in Wide-Area Interconected Power System. Int. J. Emerg. Technol. Adv. Eng. 2014, 4, 412–417. [Google Scholar]
  3. Holt, M.; Rehtanz, C. Optimizing line-voltage-regulators with regard to power quality. Electr. Power Syst. Res. 2021, 190, 106654. [Google Scholar] [CrossRef]
  4. Calderón, J.G.; Castellanos, R.; Ramirez, M. Stability Enhancement of an Industrial Power System by AVR Gain Readjustment. IEEE Lat. Am. Trans. 2017, 15, 663–668. [Google Scholar] [CrossRef]
  5. Bhatt, V.; Bhongade, S. Design Of PID Controller In Automatic Voltage Regulator (AVR) System Using PSO Technique. Int. J. Eng. Res. Appl. (IJERA) 2013, 3, 1480–1485. [Google Scholar]
  6. D’Agostino, F.; Kaza, D.; Martelli, M.; Schiapparelli, G.P.; Silvestro, F.; Soldano, C. Development of a Multiphysics Real-Time Simulator for Model-Based Design of a DC Shipboard Microgrid. Energies 2020, 13, 3580. [Google Scholar] [CrossRef]
  7. Musa, B.U.; Kalli, B.M.; Shettima, K. Modeling and Simulation of Lfc and Aver with PID Controller. Int. J. Eng. Sci. Invent. 2013, 2, 54–57. [Google Scholar]
  8. Park, C.-Y.; Kwon, J.-M.; Kwon, B.-H. Automatic voltage regulator based on series voltage compensation with ac chopper. IET Power Electron. 2012, 5, 719–725. [Google Scholar] [CrossRef]
  9. Khan, I.A.; Alghamdi, A.S.; Jumani, T.A.; Alamgir, A.; Awan, A.B.; Khidrani, A. Salp Swarm Optimization Algorithm-Based Fractional Order PID Controller for Dynamic Response and Stability Enhancement of an Automatic Voltage Regulator System. Electronics 2019, 8, 1472. [Google Scholar] [CrossRef] [Green Version]
  10. Micev, M.; Calasan, M.; Oliva, D. Fractional Order PID Controller Design for an AVR System Using Chaotic Yellow Saddle Goatfish Algorithm. Mathematics 2020, 8, 1182. [Google Scholar] [CrossRef]
  11. Sikander, A.; Thakur, P. A new control design strategy for automatic voltage regulator in power system. ISA Trans. 2020, 100, 235–243. [Google Scholar] [CrossRef] [PubMed]
  12. Micev, M.; Ćalasan, M.; Ali, Z.M.; Hasanien, H.M.; Abdel Aleem, S.H.E. Optimal design of automatic voltage regulation controller using hybrid simulated annealing—Manta ray foraging optimization algorithm. Ain Shams Eng. J. 2021, 12, 641–657. [Google Scholar] [CrossRef]
  13. Du, K.L.; Swamy, M.N.S. Neural Networks and Statistical Learning; Springer: London, UK, 2013. [Google Scholar]
  14. Singh, M.; Singh, I.; Verma, A. Identification on Non Linear Series-Parallel Model Using Neural Network. MIT Int. J. Electr. Instrum. Eng. 2013, 3, 21–23. [Google Scholar]
  15. Mak, M.W.; Ku, K.W.; Lu, Y.L. On the improvement of the real time recurrent learning algorithm for recurrent neural networks. Neurocomputing 1999, 24, 13–36. [Google Scholar] [CrossRef]
  16. Khalil, A.; Peng, A.S. An Accurate Method for Delay Margin Computation for Power System Stability. Energies 2018, 11, 3466. [Google Scholar] [CrossRef] [Green Version]
  17. Machowski, J.; Lubosny, Z.; Bialek, J.W.; Bumby, J.R. Power System Dynamics: Stability and Control, 3rd ed.; John Wiley & Sons: West Sussex, UK, 2020. [Google Scholar]
  18. Shayeghi, H.; Dadashpour, J. Anarchic Society Optimization Based PID Control of an Automatic Voltage Regulator (AVR) System. Electr. Electron. Eng. 2012, 2, 199–207. [Google Scholar] [CrossRef] [Green Version]
  19. Modabbernia, M.; Alizadeh, B.; Sahab, A.; Moghaddam, M.M. Robust control of automatic voltage regulator (AVR) with real structured parametric uncertainties based on H and μ-analysis. ISA Trans. 2020, 100, 46–62. [Google Scholar] [CrossRef]
  20. Çelik, E.; Durgut, R. Performance enhancement of automatic voltage regulator by modified cost function and symbiotic organisms search algorithm. Eng. Sci. Technol. Int. J. 2018, 21, 1104–1111. [Google Scholar] [CrossRef]
  21. Gaing, Z.L. A particle swarm optimization approach for optimum design of PID controller in AVR system. IEEE Trans. Energy Convers. 2004, 19, 384–391. [Google Scholar] [CrossRef] [Green Version]
  22. Rajasekaran, S.; Pai, G.A.V. Neural Networks, Fuzzy Systems and Evolutionary Algorithms: Synthesis and Applications; PHI Learning Private Limited: Delhi, India, 2017. [Google Scholar]
  23. Espitia, H.; Soriano, J.; Machón, I.; López, H. Design Methodology for the Implementation of Fuzzy Inference Systems Based on Boolean Relations. Electronics 2019, 8, 1243. [Google Scholar] [CrossRef] [Green Version]
  24. Espitia, H.; Soriano, J.; Machón, I.; López, H. Compact Fuzzy Systems Based on Boolean Relations. Appl. Sci. 2021, 11, 1793. [Google Scholar] [CrossRef]
  25. Chow, T.W.S.; Cho, S.Y. Neural Networks and Computing: Learning Algorithms and Applications; Series in Electrical and Computer Engineering; Imperial College Press: London, UK, 2007; Volume 7. [Google Scholar]
Figure 1. Example of a fully connected recurrent neural network using m inputs, and r processing nodes.
Figure 1. Example of a fully connected recurrent neural network using m inputs, and r processing nodes.
Processes 09 00947 g001
Figure 2. Schematic diagram of an AVR system.
Figure 2. Schematic diagram of an AVR system.
Processes 09 00947 g002
Figure 3. Block diagram of the AVR system.
Figure 3. Block diagram of the AVR system.
Processes 09 00947 g003
Figure 4. General diagram of the control system.
Figure 4. General diagram of the control system.
Processes 09 00947 g004
Figure 5. Membership functions used.
Figure 5. Membership functions used.
Processes 09 00947 g005
Figure 6. Fuzzy control system scheme.
Figure 6. Fuzzy control system scheme.
Processes 09 00947 g006
Figure 7. Optimization process results for the output y ( t ) .
Figure 7. Optimization process results for the output y ( t ) .
Processes 09 00947 g007
Figure 8. Optimization process results for the control action u ( t ) .
Figure 8. Optimization process results for the control action u ( t ) .
Processes 09 00947 g008
Figure 9. Optimization process results for the objective function J T ( k ) .
Figure 9. Optimization process results for the objective function J T ( k ) .
Processes 09 00947 g009
Figure 10. Results of y ( t ) and u ( t ) for case PQ3.
Figure 10. Results of y ( t ) and u ( t ) for case PQ3.
Processes 09 00947 g010
Figure 11. Detail results of y ( t ) for case PQ3.
Figure 11. Detail results of y ( t ) for case PQ3.
Processes 09 00947 g011
Table 1. Compact fuzzy controller f i configurations.
Table 1. Compact fuzzy controller f i configurations.
Controller f 1 f 2 f 3 f 4 f 5
SOYesYesYesYesYes
PIYesNoYesYesNo
PDNoNoYesYesNo
PIDYesNoYesYesYes
Table 2. Results of the optimization process.
Table 2. Results of the optimization process.
ConfigurationPQ1PQ2PQ3PQ4
SO7.15607.15727.221611.5680
PI7.85647.85767.919612.1120
PD7.71727.71927.788012.0156
PID7.08207.08327.148811.5252
Table 3. Performance metrics.
Table 3. Performance metrics.
ConfigurationPQ1PQ2PQ3PQ4
Metric % e ss % PO T s % e ss % PO T s % e ss % PO T s % e ss % PO T s
SO1.30162.922703.26241.14152.919703.26230.432092.89013.261013.82122.64183.2444
PI1.79710.715512.53431.63940.714332.53490.089320.70262.541113.15340.602142.599
PD1.55541.038202.39771.39941.036602.39830.136111.02152.404513.30990.88582.4640
PID1.20323.512503.29281.04293.509103.29280.532433.47473.293413.93993.1823.29510
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Espitia, H.; Machón, I.; López, H. Optimization of a Fuzzy Automatic Voltage Controller Using Real-Time Recurrent Learning. Processes 2021, 9, 947. https://0-doi-org.brum.beds.ac.uk/10.3390/pr9060947

AMA Style

Espitia H, Machón I, López H. Optimization of a Fuzzy Automatic Voltage Controller Using Real-Time Recurrent Learning. Processes. 2021; 9(6):947. https://0-doi-org.brum.beds.ac.uk/10.3390/pr9060947

Chicago/Turabian Style

Espitia, Helbert, Iván Machón, and Hilario López. 2021. "Optimization of a Fuzzy Automatic Voltage Controller Using Real-Time Recurrent Learning" Processes 9, no. 6: 947. https://0-doi-org.brum.beds.ac.uk/10.3390/pr9060947

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop