Hopfield Neural Network Optimized Fuzzy Logic Controller for MaximumPower Point Tracking in a Photovoltaic System

This paper presents a Hopfield neural network (HNN) optimized fuzzy logic controller (FLC) for maximum power point tracking in photovoltaic (PV) systems. In the proposed method, HNN is utilized to automatically tune the FLC membership functions instead of adopting the trial-and-error approach. As in any fuzzy system, initial tuning parameters are extracted from expert knowledge using an improved model of a PV module under varying solar radiation, temperature, and load conditions. The linguistic variables for FLC are derived from, traditional perturbation and observation method. Simulation results showed that the proposed optimized FLC provides fast and accurate tracking of the PV maximum power point under varying operating conditions compared to that of the manually tuned FLC using trial and error.


Introduction
Because of the demand for electric energy and environmental issues such as pollution and the effects of global warming, renewable energy sources are considered as an option for generating clean energy technologies [1].Among them, the photovoltaic (PV) energy from solar radiation has received great attention, as it appears to be one of the most promising renewable energy sources in the world [2].PV systems have been developed to supply clean energy for fulfilling the energy demand required by the modern society.However, the widespread use of PV systems poses several challenges such as increasing the efficiency of energy conversion, ensuring the reliability of power electronic converters, and meeting the requirements for grid connection [3].One step to overcome the problem of low energy conversion efficiency of PV modules and to get the maximum possible power is to extract maximum power from the PV system at every instant of time.To achieve this, it is mandatory to operate the PV systems close to its maximum power point (MPP).Maximum power point tracking (MPPT) system is an electronic system that plays a vital role in operating the PV modules in a manner that it produces it is maximum power according to the situation [4].Many MPPT control strategies have been elaborated in the literature, starting with simple techniques such as voltage and current feedbackbased MPPT to more improved power-based MPPT such as the perturbation and observation (P&O) technique and the incremental conductance technique [5].
Recently, fuzzy logic has been applied for tracking the MPP of PV systems in [6][7][8][9] because it has the advantages of being robust, design simplicity, and minimal requirement for accurate mathematical model.It is found that fuzzy logicbased P&O and hill climbing MPPT methods perform batter due to optimized perturbation.However, the fuzzy methods depend on careful selection of parameter, definition of membership function, and the fuzzy rules table.Developing fuzzy method also involves expert knowledge and experimentation in selecting parameters and membership functions.For this reason, adaptive fuzzy logic control [10] and parameter optimization techniques such as genetic algorithm [11] and particle swam optimization [12,13] have been introduced to overcome the problem in MPPT algorithms.
A number of studies on MPPT have concentrated on the application of artificial neural network (ANN) [14,15].In most of these ANN-based methods, large number of field data considering atmospheric conditions are required to train the ANN.Moreover, the main problem of ANNbased methods are that it is system dependent and cannot be International Journal of Photoenergy implemented for PV arrays with different characteristics.In a related work, a voltage-based MPPT using ANN has been developed in which an optimal instantaneous voltage factor was determined from a trained ANN [16].The inputs of the ANN consist of temperature module and solar irradiation.
The combined use of fuzzy logic and ANN to track maximum power point in PV systems can be found in [17,18].In this method, ANN is trained offline using experimental data to define a reference voltage, which is the voltage at the maximum power point according to the PV array characteristic.The reference voltage is then compared to the instantaneous array voltage to generate a signal error.The signal error and change of the error are considered as the FLC inputs.The FLC generates a duty cycle value for the pulse width modulation (PWM) generator.The PWM is then  applied to the switching of the boost converter connected to a PV array.A drawback of this method is that it needs much data for offline training.
In this paper, a new variant of intelligent technique is proposed and used together with fuzzy logic-based MPPT controller in a PV system.Here, the fuzzy logic MPPT controller is integrated with the Hopfield neural network (HNN) to optimize the membership function of the fuzzy system.The HNN has been applied to various fields since Hopfield proposed the model [19,20].In optimization problems, the HNN has a well-demonstrated capability of 0 1 g(λu) Figure 5: Curve of a sigmoid function.finding solutions to complex tasks.HNN has been applied to solve optimization problems based on convergence of the energy function and coefficients weighting [21][22][23].

Maximum Power Point Tracking and Fuzzy Logic
PV is not a constant DC energy source but has variation of output power, which depends strongly on the current drawn by the load.Besides, PV characteristic also changes with temperature and irradiation variation [4].Thus, the output voltage (V ) of PV varies with the current (I). Figure 1 shows the characteristic of a 200 W Sanyo PV module [24].The module can be used as a single system or it can be connected to other similar modules to increase the voltage and current.
In multimodule systems, PV modules are wired in series and parallel to form a PV array.From Figure 1, it can be seen that the PV module voltage varies from 0 V until the open circuit voltage (V oc ) of the module.Similarly, the current varies from 0 A until the short circuit current (I sc ) of the module.The V oc and I sc of a PV module also depend on temperature and solar irradiation.For any PV system, the output power is increased by tracking the maximum power point (MPP) of the system.To achieve this, an MPPT controller is required to track the optimum power of the PV system and it is usually connected to a boost converter located between the PV panel and load as shown in Figure 2 [6][7][8].Many different control techniques such as the P&O and fuzzy logic are used in the MPPT controller.
The main components of a fuzzy logic controller (FLC) are fuzzification, fuzzy rules inference, and defuzzification as shown in Figure 3.The input variables to the FLC are the change in power of PV (ΔP k pv ) array and direction of change in duty cycle (ΔD k ) of boost converter whereas the output of the FLC is the change of the duty cycle that must be applied to control boost converter (ΔD k+1 ref ).The formulation for ΔP k pv and ΔD k can be expressed as follows: where, P k pv : PV array output power at kth iteration, V k pv : PV array output voltage at kth iteration, I k pv : PV array output current at kth iteration, D k : measured duty cycle to control switch of DC-DC converter at kth iteration, D k+1 ref : duty cycle that must be applied to control switch of DC-DC converter at (k + 1) th iteration, Δ: a small change.
The universe of discourse for the first input variable (ΔP k pv ) is assigned in terms of its linguistic variable by using seven fuzzy subsets which are denoted by negative large (NL), negative medium (NM), negative small (NS), zero (Z), positive small (PS), positive medium (PM), and positive large (PL).The universe of discourse for the second input variable (ΔD k ) defines the changes in direction of duty cycle which is classified into three fuzzy sets, namely, negative (N), zero (Z), and positive (P).The output variable (ΔD k ref ) is assigned in terms of its linguistic variable by using nine fuzzy subsets which are denoted by negative double large (NLL), negative large (NL), negative medium (NM), negative small (NS), zero (Z), positive small (PS), positive medium (PM), positive large (PL), and positive double large (PLL).Then, the fuzzy rules are generated as shown in Table 1 with ΔP k pv and ΔD k as inputs while ΔD k ref as the output.This table is also known as fuzzy associative matrix (FAM).The fuzzy inference of the FLC is based on Mamdani's method, which is associated with the max-min composition.The defuzzification technique is based on the centroid method, which is used to compute the crisp output, ΔD k ref .

Design of Optimized FLC with HNN
The proposed design is to develop optimal membership function of the FLC especially for MPPT in PV systems application.The HNN representation and the integration of HNN and FLC is described in the following subsections.

Subsystem
Step 5 Step 3 Step 2 Step

Hopfield Neural Network Representation.
The HNN network is useful for associative memory and optimization in a symmetrical structure.The basic structure of the HNN is shown in Figure 4 [25].
The HNN uses a two-state threshold "neuron" that follows a stochastic algorithm where each neuron, or processing element, N i has two states with values either 0 or 1.The inputs of each neuron come from two sources; external inputs, I i , and inputs from other neurons, N j .The total input to neuron N i is given by  where in − N i : total input to neuron i, w ji : synaptic interconnection strength from neuron N i to neuron N j , I i : external input to neuron N i , o − N j : output of neuron N j .Each neuron samples its input at random times.It changes the value of its output or leaves it fixed according to a threshold rule with thresholds θ i : Then the energy function of the HNN is defined as The change in E due to the changing state of neuron N i by Δo − N i is given by where Δo − N i is the change in the output of neuron N i .
The continuous and deterministic model of the HNN is based on continuous variables and responses but retains all of the significant behaviors of the original model described above.The output variable o − N i for neuron N i has values in the range of 0 ≤ o N i ≤ 1 and the input-output function is a continuous and monotonically increasing function of the input in − N i to neuron N i .The dynamics of the neurons is defined as [18] du i dt where τ: a constant, u i : input of HNN.A typical output of neuron N i is a sigmoid function as shown in Figure 5. Mathematically it is given by where λ is the gain that determines the shape of the sigmoid function.
The energy function of the continuous HNN is similarly defined as [19]  and it is change in energy is given by dE/dt is always less than zero because g is a monotonically increasing function.Therefore, the network solution moves in the same direction as the decrease in energy.The solution seeks out a minimum of E and comes to a stop at stability point.

3.2.
Integrating HNN and FLC.Despite using expert knowledge in the formulation of the inference rules and the membership functions of FLC, there are still some defects such as center of fuzzification and range of the fuzzification.To improve these defects, the proposed FLC uses HNN to find the optimal membership functions which is achieved by considering the following steps.
(1) Defining neuron for the HNN.In the design of the proposed optimal FLC, two inputs, ΔP k pv and direction ΔD k , and one output, ΔD k ref , are used as described before.For simplicity, the design is only based on membership functions of ΔP k pv and ΔD k ref .
ΔP k pv is described with seven membership functions, as illustrated in Figure 6, and ΔD k ref is described with nine membership functions, as illustrated in Figure 7.In Figure 6, the centers of ΔP k pv membership function are x 1 , x 2 , x 3 , x 4 , x 5 , x 6 , and x 7 while in Figure 7   When comparing the values of neurons in Figures 6 and  7, the following constraints should be satisfied.
(2) Defining objective function for the optimization problem.The goal of MPPT is to achieve ΔP k pv = 0 and ΔD k ref = 0. Therefore, the quadratic criterion to be minimized is where, E: energy function to be minimized, A, B: constants.From (11), the first part of E, which is defuzzification of the universe of ΔP k pv using a centroid function as: where μ(ΔP k pv ) i is a membership value of x i .E 1 depends only on neurons N i (i = 1, 2, 3, 4, 5, 6, 7).Knowing that the left side of ΔP k pv in the membership function of the first input is μ(ΔP k pv ) i / = 0 for x 1 , x 2 , x 3 , then, ΔP k pv can be rewritten as Then, the first half of energy function E 1 can be rewritten as E 1a : Considering that the right side and center of ΔP k pv in the membership function of the first input of FLC is μ(ΔP k pv ) i = 0 for x 1 , x 2 , x 3 , hence ΔP k pv can be rewritten as: Similarly, the second half of energy function E 1 is rewritten as E 1b : Since the first part of E 1 is the summation of left side (E 1a ) and right side (E 1b ) of the ΔP k pv membership function, E 1 can be expressed as: where o N = 0 for i = 4, j = 4.The second part, E 2 = (1/2)B(ΔD k ref ) 2 is related to the output of FLC and depends only on neurons N i (i = 8,9,10,11,12,13,14,15,16).The ΔD k ref can be defined by defuzzification by using the centroid method and is written as: Similar to the equations shown in obtaining E 1 , E 2 can be expressed as where N = 0 for i = 12, j = 12 and μ(ΔD k ref ) is the membership value of z n .
Finally, the total energy function E is expressed as where N = 0 for i = 4, 12; j = 4, 12.
By comparing (8) with (20), the weight matrix of neurons N i to N j in the HNN is derived and given as: where, , for i and j = 1, 2, 3, 4, 5, 6, 7, , for i, j = 8, 9, 10, 11, 12, 13, 14, 15, 16,   w i j = 0, for other. ( (3) Design for Physical Implementation.The physical implementation of MPPT for PV systems using Hopfield optimized FLC (HFLC) is described in terms of a block diagram as shown in Figure 8.As shown in Figure 8, the system consists of PV array, DC-DC converter, load, and control system.The control system consist of voltage and current measurement system, controlled pulse with modulation (PWM) generator with HFLC, and a switching driver circuit to gain the PWM.At the initial step (k = 0), the control system generates square wave signals with a small duty cycle (D), of value 10%.In the next step (k = 1), the value of D is increased by ΔD k

Simulation Results
The performance of the proposed HFLC under different operating conditions is validated using the MAT-LAB/Simulink software.In the PV model shown in Figure 9, there are two groups of PV arrays connected in parallel.Every group consists of 3 PV modules connected in series.
The PV module parameters are obtained from the Sanyo HIP-200BA3 PV technical datasheet [24].In the simulations, first the characteristics of the PV module are validated and then the performance of the HFLC under various conditions is evaluated to investigate the effectiveness of the HFLC method.

Validation of PV Module Simulation. Figures 10 and 11
show the results of the I-V characteristics of the simulated PV module as a function of irradiation and temperature, respectively.It can be observed from the above figures that the I-V curves of the simulated PV module are quite similar to the I-V curves of the Sanyo HIP-200BA3 PV module provided by the Sanyo manufacturer in Figure 1.Therefore, it is quite reasonable to use the PV module model to verify the performance of the proposed HFLC-based MPPT controller under simulation environment.

Performance of MPPT by Using HFLC. Figures 12 and 13
show the performance of the HFLC in finding the maximum power point (MPP) of the PV system shown in Figure 9 under varying irradiations and temperatures, respectively.From the figures, the MPP obtained from HFLC is compared with the theoretical MPP.The results of the MPP clearly show that both MPPs are very close to each other.
To further demonstrate the performance of the HFLC MPPT controller, simulations were performed under the following test conditions.(i) Constant temperature at 25 • C and changing the solar radiation slowly and drastically.
(ii) Constant solar radiation at 1000 W/m 2 and changing the temperature slowly and drastically.
(iii) Constant temperature at 25 • C and considering partial shading and change in solar radiation.
The MPPT controller was also tested using the conventional FLC and the P&O MPPT methods.

Effect of Changing the Solar Radiation.
To analyze the effect of solar radiation, simulations were carried out under various solar irradiations but at constant temperature of 25 • C. Figure 14 shows the change in solar irradiation from 1000 W/m 2 to 900 W/m 2 .Figure 15 shows the PV output power when subjected to the changing solar irradiations.From Figure 15, it can be seen that the HFLC, and FLC-based MPPT gives greater PV output powers than the P&O-based MPPT.
Figure 16 shows the sudden change in solar irradiation from 1000 W/m 2 to 950 W/m 2 while Figure 17 shows the response of the MPPT controller in terms of PV output From Figure 19, it can be noted that for slow temperature changes, the MPPT controlled by HFLC and FLC gives higher PV output power than the P&O method especially at the transient state.While in the case of drastic change in temperature, the MPPT controlled by HFLC achieve the highest PV output at the transient state as shown in Figure 21.

Effect of Partial Shaded Solar Irradiation.
Simulation is also performed to illustrate the effectiveness of the MPPT of PV systems under some partial shading case.In this case, it is assumed that a half of the PV array receives constant solar irradiation of 1000 W/m 2 and the other half with shading solar irradiation which changes from 1000 W/m 2 to 950 W/m 2 .This condition is depicted in Figure 22.The power harvested from the PV array for this case is shown in Figure 23.From Figure 23, it can be seen that performance of MPPT controlled by HFLC is the best among the other compared methods.
Another case describes a situation of solar irradiation changing from 400 W/m 2 to 1000 W/m 2 , to a half of the PV array while the other half receiving constant 400 W/m 2 under shading as shown in Figure 24.The PV output power for this case is shown in Figure 25.From the figure, it can be noted that the FLC and P&O failed to track the MPP correctly.However, the proposed HFLC MPPT method successfully finds the MPP around 590 W as shown in Figure 25.
Figures 26, 27, 28, and 29 show the characteristic curves of voltage versus power (V -P) of the modeled PV array under various partial shading conditions described in Table 2.
From Figures 26 and 27, it can be seen that the HFLC method is accurate in finding the MPP (590 W) while the conventional FLC and P&O methods failed to do so.Generally, the conventional FLC and P&O methods just find local maximum power point.However, for case as depicted in Figures 28 and 29, all of the MPPT methods correctly tracked the MPP.

Conclusion
A new Hopfield optimized FLC for MPPT of PV system is proposed in which improvement is made by applying HNN to find the optimal width of each fuzzification input and output of the FLC.A complete PV system with HFLC MPPT controller was modeled and implemented in Matlab/Simulink to simulate various irradiation and temperature conditions so as to verify the performance of the proposed MPPT method.Simulation results show that the proposed HFLC MPPT method is robust and accurate compared to the other conventional MPPT methods.The HFLC MPPT method successfully tracks the global maximum power point of a PV module even under partial shading.

Figure 9 :
Figure 9: MATLAB simulation model of a PV MPPT controller.

Figure 10 :
Figure 10: Effects of solar radiation at constant temperature on the PV module.

Figure 11 :
Figure 11: Effects of temperature at constant solar radiation on the PV module.

Figure 12 :Figure 13 :
Figure 12: Performance of the HFLC-based MPPT under various irradiation at constant temperature T = 25 • C.
ref that is defined by HFLC as discussed in the previous section.The value of D is always updated by the increment ΔD k ref .The value of ΔD k ref tends to change in either positive or negative direction as |ΔD k ref | decline towards zero.

Table 1 :
Fuzzy rules for the proposed FLC.
, the centers of ΔD k ref membership function are z 1 , z 2 , z 3 , z 4 , z 5 , z 6 , z 7 , z 8 , and z 9 .Based on the number of centers of ΔP k pv and ΔD k ref , the proposed HNN consists of 16 neurons with variables given

Table 2 :
Various partial shaded solar irradiation.Effect of Change in Temperature.This simulation is carried out to illustrate the performance of the MPPT methods under constant solar irradiation of 1000 W/m 2 and changes in temperature.Figures 18 and 20 depict the slow and sudden changes in temperature, respectively.Figures 19 and 21 show the corresponding PV output powers during slow and sudden changes in temperature, respectively.