Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Nonlinear adaptive NeuroFuzzy feedback linearization based MPPT control schemes for photovoltaic system in microgrid

  • Muhammad Awais,

    Roles Conceptualization, Data curation, Investigation, Methodology, Software, Writing – original draft

    Affiliation Department of Electrical and Computer Engineering, COMSATS University Islamabad, Abbottabad, Pakistan

  • Laiq Khan ,

    Roles Supervision

    laiq@cuiatd.edu.pk

    Affiliation Department of Electrical and Computer Engineering, COMSATS University Islamabad, Abbottabad, Pakistan

  • Saghir Ahmad,

    Roles Formal analysis, Writing – review & editing

    Affiliation Department of Electrical and Computer Engineering, COMSATS University Islamabad, Abbottabad, Pakistan

  • Sidra Mumtaz,

    Roles Validation, Visualization

    Affiliation Department of Electrical and Computer Engineering, COMSATS University Islamabad, Abbottabad, Pakistan

  • Rabiah Badar

    Roles Data curation, Resources

    Affiliation Department of Electrical and Computer Engineering, COMSATS University Islamabad, Abbottabad, Pakistan

Abstract

Renewable energy resources connected to a single utility grid system require highly nonlinear control algorithms to maintain efficient operation concerning power output and stability under varying operating conditions. This research work presents a comparative analysis of different adaptive Feedback Linearization (FBL) embedded Full Recurrent Adaptive NeuroFuzzy (FRANF) control schemes for maximum power point tracking (MPPT) of PV subsystem tied to a smart microgrid hybrid power system (SMG-HPS). The proposed schemes are differentiated based on structure and mathematical functions used in FRANF embedded in the FBL model. The comparative analysis is carried out based on efficiency and performance indexes obtained using the power error between the reference and the tracked power for three cases; a) step change in solar irradiation and temperature, b) partial shading condition (PSC), and c) daily field data. The proposed schemes offer enhanced convergence compared to existing techniques in terms of complexity and stability. The overall performance of all the proposed schemes is evaluated by a spider chart of multivariate comparable parameters. Adaptive PID is used for the comparison of results produced by proposed control schemes. The performance of Mexican hat wavelet-based FRANF embedded FBL is superior to the other proposed schemes as well as to aPID based MPPT scheme. However, all proposed schemes produce better results as compared to conventional MPPT control in all cases. Matlab/Simulink is used to carry out the simulations.

Introduction

The energy demand of the globe is mainly fulfilled by fossil fuel. Increasing energy demand and limitation of fossil fuel supplies boost the cost of electricity. The environmental hazards due to greenhouse emission and scarcity of fossil fuel supplies diverted the focus towards renewable energy resources. Renewable power is clean, sustainable, green, economical, and durable. The productivity of renewable energy resources depends on meteorological conditions. Thus, a single stand-alone source is unable to supply continuous reliable energy. Therefore several renewable and non-renewable resources are integrated to form a single HPS. Solar energy is a huge reservoir of green energy blessed to this planet [1, 2]. To utilize it most reliably and efficiently needs its conversion into electrical energy [3]. Many researchers struggle for the extraction of maximum power from photovoltaic (PV) arrays as its demand increases in both stand-alone and grid-connected modes [4]. The grid-connected mode is trending these days due to drawbacks like high cost, heavy batteries, and regular maintenance related to traditional stand-alone PV systems [5, 6].

Different topologies, circuits, and control algorithms are required for the extraction of maximum power from PV. Maximum power tracking from a PV system is of interest for a long time and is usually preferred using a boost converter due to its advantages over other techniques [7, 8]. PV and IV characteristics of the PV system vary according to the atmospheric conditions, thus maximum power tracking becomes difficult [5, 8, 9]. Numerous MPPT techniques are practiced and proposed for the extraction of maximum power from the PV system that is then supplied to the grid through inverters [10].

Literature studies show that various conventional and intelligent methods are used for the MPPT of the PV system. The conventional techniques include perturb and observe (P&O) [1114], fractional open-circuit voltage (OCV) [15], Hill Climbing algorithm (HCA) [16], incremental conductance (IC) [17], and ripple correlation (RC) [18] method. P&O, IC, and HCA have the drawback of oscillations around maximum power point (MPP), thus loses power and also show degraded energy conversion efficiency. Under fast varying conditions, P&O and IC have a slow response. RC algorithm needs prior knowledge of the phase relationship between the DC power ripple and DC voltage, to determine the MPP. Also, the implementation of the RC method is quite complex. The accuracy and efficiency of both OCV and SCC are low.

Advanced MPPT control techniques include artificial neural network (ANN) [19], fuzzy logic control (FLC) [20], feedback linearization [21], adaptive NeuroFuzzy [22], wavelet-based NeuroFuzzy [8, 23], and adaptive feedback linearization based NeuroFuzzy [9]. The major drawback of FLC is the exponential growth in the number of membership functions and the fuzzy if-then rules with an increase in variables. The ANN model requires periodic training to ensure convergence to the accurate MPP. Feedback linearization is a recent technique used for the MPPT problem [5, 9]. This technique decouples a nonlinear system into linear subsystems, thus linear control laws can be implemented on these decoupled subsystems. Feedback linearization has its applications in the control of robotics [24, 25], HVDC link [26], motor control [27], etc. However, for un-modeled system dynamics, the classical feedback linearization control has limited robustness and stability [25]. Adaptive feedback linearization-based NeuroFuzzy algorithms are used to nullify the drawbacks of classical feedback linearization. The NeuroFuzzy algorithm combines the reasoning method of humans in fuzzy systems with the learning abilities of neural networks. This hybrid adaptive scheme can deal with system nonlinearities, uncertainties, and fluctuations. Therefore, in an adaptive feedback linearization-based NeuroFuzzy control scheme, the adaptive feedback linearization control is applied to the nonlinear system identified by the hybrid NeuroFuzzy inference systems. In [8, 9] and [23] NeruoFuzzy schemes are presented for the identification and control of grid-connected PV subsystems. A feedforward NeuroFuzzy scheme is presented having static mapping thus time-domain dynamical responses require a larger number of neurons [28]. It limits the number of inputs and can explore the characteristics of local structure [2830]. On the other hand, recurrent NeuroFuzzy control schemes consist of feedback and feedforward connections between layers [28]. The neurons in this scheme form complicated dynamics and use their natural temporal operation to deal with time-varying input/output [28, 29].

Mismatching of PV cells along with the weather and environmental changes like shadows of buildings, dust, moving clouds, and trees give rise to partial shading conditions (PSC), which result in power loss and produce numerous local MPPs (LMPPs) [3133]. Numerous conventional and non-conventional techniques were developed to tackle PSC in the near past but they have drawbacks of lager power fluctuation, lower power output, and complexity of control design in some cases [3237].

This paper presents a comparison of the performance of four different adaptive feedback linearization (FBL) techniques incorporated with full recurrent adaptive NeuroFuzzy (FRANF) based controllers for a PV system in a grid-integrated SMG-HPS for three different cases.

Seven-layered full recurrent adaptive NeuroFuzzy structure embedded with four different mathematical functions and wavelets is used to estimate the nonlinear functions of FBL control. FRANF structures are based on Standard Additive Model, Fourier Series, Mexican Hat Wavelet, and Chebyshev Wavelet for the estimation purpose. An online learning algorithm based on the gradient-decent method is applied to update all the parameters of the FRANF structure for each proposed control scheme adaptively. Three different cases are used on the same system to test the performance of proposed controllers. The extraction of maximum power under varying conditions from the PV system in these scenarios is challenging. The comparison of the proposed control schemes is based on the power quality of the SMG-HPS as well as extracted power and performance indexes obtained during the simulation process. The performance of the proposed control schemes is evaluated against the adaptive PID (aPID) control scheme.

The article comprised of four major sections. Section 1 describes the testbed for this research work. Section 2 gives detail of the proposed control schemes for the PV subsystem. Section 3 and 4 provide details of results and conclusion respectively.

1 SMG-HPS and PV subsystem

The sources in the testbed for this research work consists of wind-turbine (WT), electrolyzer, micro-turbine (MT), solid oxide fuel cell (SOFC), PV, Micro-hydro (MH), Bio-mass (BM), super-capacitor (SC), batteries, and utility grid (UG). The loads are prosumer, residential load, charging station (CS), and plug-in hybrid-electric vehicles (PHEVs).

Renewable resources like WT, electrolyzer, PV, SOFC, SC, and batteries are connected to DC link and deliver power to the connected converters. These converters are connected to the AC bus that has UG, MT, MH, CS, BM, along with all types of loads attached to it as shown in Fig 1. Specifications of all the sources are described in section 3 Table 1. All the source converters except PV are being controlled by aPID.

1.1 PV cell model

A simple pn-junction that converts solar irradiation into electric energy is known as PV cell [32]. It is comprised of current source, parallel diode and a series resistor, which is further connected to boost converter [8, 32]. PV cells are combined in clusters which are then connected to one-another is series and parallel fashion to obtain the desired power level [9, 32, 33]. The relationship between output voltage and current of PV cell is given as: [6, 8, 19, 32, 33] (1) where all the symbols are defined in Nomenclature. Photo-current, ip can be determined by solar irradiation as: (2) Beside this, saturation current, is of PV cell varies with temperature and can be related as: (3) The Eqs (1), (2) and (3) are used to design PV system and also show that the output of a PV array depends on solar irradiance and temperature of the environment [32].

2 Adaptive feedback linearization embedded NeuroFuzzy MPPT for PV system

2.1 Control law design

Feedback linearization is a tool that transforms nonlinear system dynamics into linear ones algebraically, either fully or partly, hence linear control techniques can be applied. nonlinearities can be eliminated from any nonlinear system represented in a companion form as [38]: (4) where y ∈ ℜ is plant output, n is relative degree of system, f(x) and g(x) are unknown nonlinear functions, uMPPT ∈ ℜ is control input, and is the state-space vector. The control problem is to find uMPPT that assures y(t) follows the desired trajectory yd(t). If a new input v represents the plant’s input then: (5) Using the control law of Eq (5) in Eq (4), nonlinear terms will be canceled and input-output integral form is obtained as: (6) The nonlinear functions in Eq (5) are estimated using NeuroFuzzy systems. represents the estimated f(x) and represents g(x). The adaptive feedback control law can be rewritten as: (7) The control law of Eq (7) is based on control input v = −Kv℘ − YD and identified nonlinear functions. where Kv is constant and ℘ the tracking error is given as [38]: (8) where ∧ = [λ1 λ2 ⋯  λn-1] is constant weight vector and e = xxd represents error matrix.

Online adoption of ∧ ensures occupation of poles of sn−1 + λn−1sn−2 + ⋯ + λ1 in left half of complex plane. To bring the tracking error to zero, the following control law is entertained by identification of and through FRANF, where is obtained by using saturation block and identifier is set to give 0.1 as output if during identification. (9) To generate appropriate control law, FBL control coefficients ∧ are updated using n LMS algorithm, as: where , ζ = e(n−1) and is estimation of ∧. Fig 2 shows structure of Adaptive Feedback Linearization embedded with FRANF. The control law requires estimated and of Eq (7) which are identified using NeuroFuzzy identification.

Fig 3 shows the internal control system and the closed-loop adaptive feedback control strategy.

thumbnail
Fig 3. Adaptive FBL embedded FRANF control for MPPT of PV system.

https://doi.org/10.1371/journal.pone.0234992.g003

2.2 Full recurrent adaptive NeuroFuzzy identification techniques based on different NeuroFuzzy architectures

A variety of FRANF identifiers are used in order to identify the nonlinear and functions for PV system in SMG-HPS. The seven-layered FRANF system uses NeuroFuzzy concept for estimation.

Fuzzy logic uses IF-THEN rules for approximation of unknown functions using standard fuzzy model. The unknown functions and can be identified by the standard fuzzy model using a set of rules as;

Rm: IF x1 is and xn is THEN y is .

Let fuzzy logic controller has q inputs, ρ1, ρ2, …, ρq. The output of NeuroFuzzy system is given as: (10) where is the membership function, ρj and βl are adjustable parameters. It is the point in R at which achieves its maximum value. m is the number of fuzzy rules used to construct the identifier, is the jth fuzzy set corresponding to the lth fuzzy rule, and βl is centroid of the lth fuzzy set corresponding to identifier output, and . Eq (10) can be written for and using fuzzy-basis function vector ξ(x), as (11) and (12) where (13) and (14) and ξ(x) is given as (15) A number of mathematical relations and functions are available for designing a fast and robust NeuroFuzzy identifier. The following variants are used to design antecedent and the consequent part of the fuzzy logic system for this research work.

2.2.1 Antecedent part.

The transformation of continuous input variables into linguistic variables is fuzzification. A membership function is always required for the transformation. The importance of the membership function is based on its shape that translates complete information of the plant (uncertainties and nonlinearities) in fuzzy inference system. The membership function chosen for this research work is the Gaussian membership function due to the following properties:

  • local and nonlinear nature
  • smooth output
  • relation between the radial basis functions Neural Networks (NNs) and fuzzy system

Gradient-based techniques are highly suitable for use due to the continuous differentiable nature of Gaussian membership function. It is expressed as: (16) where mij and σij are the mean and variance of the ith input and jth membership function.

2.2.2 Variants of consequent part.

The consequent part generates weights based on different mathematical functions like Fourier series function, wavelet networks and polynomial NNs. The operation of consequent part takes place in parallel to antecedent part and produces final output of identifier at defuzzification layer. The variants of consequent part used for this research work are Standard Additive Model (SAM), Fourier series function and wavelet networks (Mexican hat wavelet and Chebyshev wavelet).

2.2.2.1 Standard Additive Model (SAM). SAM is an important case of additive fuzzy systems that can estimate any uninterrupted function uniformly over a closed space. In SAM the fuzzy rules are given as:

If xi = C Then yi = D such that C and D are one to one mapping of input and output spaces. Fig 4 shows the architecture of SAM. Following theorem defines SAM.

Theorem 1.1 Consider the fuzzy system is SAM: F(x) = Centeroid(D) = Centeroid Then F(x) represents convex sum of n then-part set centroids; (17) where, (18) are the convex coefficients or discrete probability weights, ck is the centroid of the then-part set and Vk is the finite positive volume.

2.2.2.2 Fourier Series Neural Networks (FsNNs). Fast convergence and accurate modeling capabilities of Fourier series NNs are known with gradient descent algorithm. Mutually orthogonal sine and cosine are the basis functions of FsNNs that guarantee better estimation and convergence.

In FsNN, the jth input signal initiates neurons with equal number of sine and cosine as activation function. Adaptation of output weights occur during learning. Overall output of the FsNN is given as: (19) where, ⊗ denotes Kronecker product, n shows the total number of inputs and . Fig 5 shows structure of FsNN. Sine and cosine are the basis functions of FsNNs that span over the infinite interval of time. They are not localized in time and have infinite energy.

2.2.2.3 Fuzzy Wavelet Neural Networks (NNs). For a better estimation of nonlinear functions, wavelet NNs were proposed as a substitute to feedforward NNs. Due to enormous neurons, NNs may get struck on the local minima that result in slower convergence of the network. To get rid of it, wavelet functions (WFs) can be used in the structure. Wavelets are waves having a limited duration and zero mean value. The localization characteristics of wavelets and the fast learning abilities of NNs results in better outcomes for complex nonlinear system modeling. The schematic diagram of wavelet NN is given in Fig 6. Two wavelet activation functions, Mexican hat and Chebyshev wavelet are used in this research work.

  • Mexican hat wavelet (MHW) is a negative normalized, non-orthogonal second derivative of Gaussian function. MHW function is expressed as; (20) where Ψi(xi) is the family of wavelets obtained by single ψ(xi) function with parameters dilation dij and tij respectively.
  • Chebyshev wavelet (CW) have four arguments; k ϵ , r = 1, 2, …, 2k−1, and ; s is the degree of the first kind Chebyshev polynomial and t ϵ [0, 1] is the normalized time. They can be expressed on the interval [0, 1] as: (21) where (22) s = 0, 1, …N − 1, where N is a fixed positive integer. Orthonormality of the system is given by the coefficients in Eq (22). Here nth degree Chebyshev polynomials orthogonal to weight function on the interval [−1, 1] is represented by that satisfy the following recursive formula: (23) (24) The weight functions of CW, generates orthogonal wavelets on dilation and translation as .

2.2.3 Proposed FRANF identifier.

The FRANF has seven layers as shown in Fig 7. The antecedent part consists of the first three layers, whereas the rest four layers are consequent part layers. The n number of input signals in the first layer is equivalent to the m number of nodes and these nodes are used for input distribution. Let and represents input and output of ith node in kth layer. The operation function of nodes and the signal propagation in each layer is given as under:

  • Layer 1 (Input Layer): It takes input variables and these inputs are transmitted to the next layer by the nodes. Feedback connections are the part of this layer that add temporal relationship in the network. (25) (26) where, i = 1, 2, 3, ……n; k is number of iterations.
  • Layer 2 (Membership Layer): Each node represents one linguistic term and computes membership degree and fuzzy set for all input signals entering into the system. Linguistics terms in each node are computed using Gaussian membership function. (27) (28) Where ij subscript shows the jth term of the ith input , where and the superscript (2) represents second layer. Also, i = 1, 2, 3, ……n; k is number of iterations, is the past value of membership function and is the recurrent weight.
  • Layer 3 (Rule Layer): The rule layer is just the product of membership functions. The number of rules, i.e. R1, R2, …Rn in this layer are equal to number of nodes and each node corresponds to one fuzzy rule. Min operator is used to compute the output signal’s value in each rule. Each is the input for the next (consequent) layer. (29) (30)
  • Layer 4 (Consequent Layer): It determines the difference in the proposed control techniques. The general description of this layer is; (31) where Hi is the output of the hidden layer (previous values of mathematical function), and Fi is the feedback weight and the superscript (4) denotes the layer 4. Each mathematical function is multiplied by the weight, wi of the NNs in fourth layer. (32) where represents one of the used mathematical function discussed in section 2.2.2.
  • Layer 5 (Defuzzification Layer 1): The product of the outputs of antecedent and consequent parts for each input is taken in this layer and then added to each other. (33)
  • Layer 6 (Defuzzification Layer 2): In this layer the summation of rules (output of antecedent) part is calculated. (34)
  • Layer 7 (Output Layer): The final output of the FRANF identifier is estimated in the seventh layer as given below: (35)
thumbnail
Fig 7. Architecture of Full Recurrent Adaptive NeuroFuzzy System.

https://doi.org/10.1371/journal.pone.0234992.g007

2.2.4 Optimization algorithm.

The training of FRANF is for the adjustment of input-output pairs or a given function by fine-tuning network parameters. Mean square error is used as cost function for the training purpose, given as [3944]; (36) Gradient descent method is used for fast cost function reduction and convergence [3944]. The general equation is given as follows: (37) where, gk is the gradient of cost function at kth iteration, γ > 0 is the learning rate and k is the iteration index. The update equations for antecedent part and consequent part are as follows:

2.2.4.1 Update equations for antecedent part parameters. The update equation for variants of Gaussian membership function is derived from following chain rule [3944]. (38) where, χ shows the variant like mean, variance and feedback weight of Gaussian membership function.

Update equation for mean, mi is; (39) Update equation for variance, σi is; (40) Update equation for recurrent weight, θi is; (41)

2.2.4.2 Update equations for consequent part. Chain rule for updating variants of consequent part is [3944]; (42) where ∃i represents one of the mathematical functions discussed in section 2.2.2 and shows the variants of mathematical function like centroid, dilation, translation and volume etc.

Update equation for centroid, ci of SAM is; (43) where, wi is the adaptive weight of consequent part obtained by Eq (51), wi(k) and wl(k) are the weights of SAM. Update equation for volume, Vi of SAM is; (44) Update equation for weight, wi(k) of FsNNs is; (45) where, wi is the weight of consequent part from Eq (51), wi(k) is the weight of FsNNs.

Chain rule for updating variants of MHW used in this research work is: (46) where here Ψi represents MHW, Zi defined below is an intermediate state variable and shows the variants of MHW like, translation, dilation, feedback weight etc. (47) Update equation for translation, ti is; (48) Update equation for dilation, di is; (49) Update equation for feedback weight, Fi of MHW is; (50) The update equations of CW are applied according to the equations given in section 2.2.2.3.

Weight of consequent part are being updated according to the following chain rule: (51) where wi represents the weight of the consequent layer. (52) The feedback weight of antecedent part of all proposed controllers is adaptive while feedback weight of consequent part of only Mexican hat wavelet is adaptive. The feedback weight of rest of the consequent parts of other controllers is closed-loop fixed gain.

2.3 Proposed controller Schemes

2.3.1 Adaptive Feedback Linearization (FBL) embedded with Full Recurrent Adaptive NeuroFuzzy (FRANF) Standard Additive Model (SAM) control.

This controller is based on the FBL scheme which is embedded with a FRANF-SAM identifier to identify unknown nonlinear functions for FBL. The antecedent part of FRANF is modeled using Gaussian membership function while the consequent part is modeled using SAM as discussed in section 2.2.2.1 and section 2.2.3. The overall explicit control law for this control scheme is given in Eq (53). (53)

2.3.2 Adaptive Feedback Linearization (FBL) embedded with Full Recurrent Adaptive NeuroFuzzy (FRANF) Fourier Series (FS) control.

This controller is based on the FBL scheme which is embedded with a FRANF-FS identifier to estimate unknown nonlinear functions for FBL. The antecedent part of FRANF is modeled using Gaussian membership function while the consequent part is modeled using FS as discussed in section 2.2.2.2 and section 2.2.3. The overall explicit control law for this control scheme is given in Eq (54). (54)

2.3.3 Adaptive Feedback Linearization (FBL) embedded with Full Recurrent Adaptive NeuroFuzzy(FRANF) Mexican Hat Wavelet (MHW) control.

This controller is based on the FBL scheme which is embedded with a FRANF-MHW identifier to identify unknown nonlinear functions for FBL. The antecedent part of FRANF is modeled using Gaussian membership function while the consequent part is modeled using MHW as discussed in section 2.2.2.3 and section 2.2.3. The overall explicit control law for this control scheme is given in Eq (55). (55)

2.3.4 Adaptive Feedback Linearization (FBL) embedded with Full Recurrent Adaptive NeuroFuzzy (FRANF) Chebyshev Wavelet (CW) control.

This controller is based on the FBL scheme which is embedded with a FRANF-CW identifier to identify unknown nonlinear functions for FBL. The antecedent part of FRANF is modeled using Gaussian membership function while the consequent part is modeled using CW as discussed in section 2.2.2.3 and section 2.2.3. The overall explicit control law for this control scheme is given in Eq (56). (56)

Results and discussion

The SMG system is developed in Matlab/Simulink R2015a for evaluation of the performance of the proposed controllers. The system is developed by using power generation sources like UG, PV, WT, electrolyzer, SOFC, MT, and the backup sources i.e. SC and batteries whose ratings are given in Table 1. The purpose of using multiple sources is to entertain the dynamic residential load and CS load.

Intelligent supervisory control is an essential part of the SMG system. It monitors the power generation and load variations during the simulation time. The supervisory control ensures the power consumption from renewable resources at a priority level during their availability period. It satisfies the load demand by shifting the load to other sources and UG during peak hours and in the absence of renewable power.

Intelligent control schemes like Adaptive FBL embedded FRANF-SAM, Adaptive FBL embedded FRANF-FS, Adaptive FBL embedded FRANF-MHW, Adaptive FBL embedded FRANF-CW are used to extract the maximum power from PV system connected in the microgrid. The performance of all the proposed control schemes is compared to aPID based MPPT control scheme.

3.1 Case studies

Three different cases are taken in this research work, e.g., (a) Step change in both solar irradiation and temperature; (b) Partial shading condition; and (c) Daily field data of solar irradiation and temperature in Islamabad, to evaluate the performance of proposed adaptive FBL embedded FRANF controllers under PSC compared to aPID control scheme.

3.1.1 Step change in both solar irradiation and temperature.

Step-changing solar irradiation and temperature profile are simulated for 24 seconds, where each second represents one hour. The solar irradiation is gradually increased in many steps to its maximum and then gradually decreased to zero. Figs 8 and 9 shows the step profile of solar irradiation and ambient temperature used for this case study.

Fig 10 shows the output power comparison of all intelligent proposed control schemes with aPID. Every step change is accurately determined by the proposed control schemes and the each proposed control scheme generates more output power compared to aPID. Furthermore, it is also clear that the power generated by three proposed controllers, i.e., adaptive FBL embedded FRANF-CW, adaptive FBL embedded FRANF-FS, and adaptive FBL embedded FRANF- SAM, seems to be very near to each other. However, the power generated by adaptive FBL embedded FRANF-MHW is maximum than other proposed schemes.

thumbnail
Fig 10. PV MPP tracked vs reference power for step change profile.

https://doi.org/10.1371/journal.pone.0234992.g010

Power error (Perror) can be described as: (57) where Pref is the reference power and PPV is tracked output power obtained from the PV system under the action of applied controllers. The maximum and average power error of all the control schemes is given in Table 2.

thumbnail
Table 2. Peak power error and average power error of all control schemes for step changing solar profile.

https://doi.org/10.1371/journal.pone.0234992.t002

Performance indexes including Integral Square Error (ISE), Integral Time Square Error (ITSE), Integral Absolute Error (IAE), and Integral Time Absolute Error (ITAE), calculated based on Perror in Eq (57) as shown in Figs 11, 12, 13 and 14. The comparison of performance indexes plots shows that the accumulative error in all schemes increases with time. The index of adaptive FBL embedded FRANF-MHW is least among all proposed controllers. Table 3 shows the values of various indexes of all the proposed controllers compared to aPID showing least to most performing from top to bottom.

The switching in the converter circuits arises harmonics in the load voltage and current that are not synchronized with the frequency of the system. It wastes power as heat and should be minimized. Voltage fluctuations and flickers are caused by higher frequency harmonics.

Fig 15 shows a comparison of the percentage change in total harmonic distortion (THD) for load current due to individual control scheme. The result shows that the percentage change in THD due to the Adaptive FBL embedded FRANF-MHW control scheme is the smallest of all and proves its better performance among other proposed controllers. The percentage change in frequency is also shown in Fig 16 for all proposed controllers. It can be observed that the percentage change in frequency due to the adaptive FBL embedded FRANF-MHW control scheme is almost flat and nearly zero.

thumbnail
Fig 15. %age THD change in load current for step changing solar profile.

https://doi.org/10.1371/journal.pone.0234992.g015

thumbnail
Fig 16. %age frequency change for step changing solar profile.

https://doi.org/10.1371/journal.pone.0234992.g016

3.1.2 Partial shading condition.

PSCs arise due to moving clouds, airplanes, dust, and shadows of the building. They cause a sudden drop in irradiation level as well as temperature ranging from a short interval of time to many hours. To evaluate the performance of proposed controllers for this type of phenomenon, PSCs are introduced in solar and temperature profiles by introducing multiplying factors at certain intervals of time as shown in Fig 17.

This suddenly changes the magnitude of irradiation level, as well as temperature, thus produce fluctuations in irradiation and temperature curve as shown in Figs 18 and 19.

thumbnail
Fig 18. Partial shading profile of solar irradiance level (W/m2).

https://doi.org/10.1371/journal.pone.0234992.g018

thumbnail
Fig 19. Partial Shading profile of ambient temperature (°C).

https://doi.org/10.1371/journal.pone.0234992.g019

Fig 20 shows the output power tracked by all the proposed control schemes compared with aPID. The sudden drops in irradiation and temperature level are successfully adopted by all the proposed control schemes, whereas, aPID produces the least accurate results. It is obvious to note that the output power produced using FRANF-CW, FRANF-SAM, and FRANF-FS is better than each other but the performance of FRANF-MHW is superior over all the other proposed control schemes. The maximum and average power errors of all the control schemes calculated according to Eq (57) are given in Table 4.

thumbnail
Fig 20. PV MPP tracked vs reference power for PSC profile.

https://doi.org/10.1371/journal.pone.0234992.g020

thumbnail
Table 4. Peak power error and average power error of all control schemes for partial shading condition.

https://doi.org/10.1371/journal.pone.0234992.t004

Performance indexes calculated based on Perror in Eq (57) are shown in the Figs 21, 22, 23, and 24. Comparative analysis based on the various performance indexes is shown in Table 5.

Fig 25 shows a comparison of the percentage change in total harmonic distortion (THD) for load current for all control schemes. It is obvious to note that the percentage change in THD due to the adaptive FBL embedded FRANF-MHW control scheme is the smallest compared to the rest of the proposed controllers. The percentage change in frequency is also shown in Fig 26 for all proposed controllers. It can be observed that the percentage change in frequency due to adaptive FBL embedded FRANF-MHW control scheme is the least among all the proposed schemes.

thumbnail
Fig 25. %age THD change in load current for partial shading solar profile.

https://doi.org/10.1371/journal.pone.0234992.g025

thumbnail
Fig 26. %age frequency change for partial shading solar profile.

https://doi.org/10.1371/journal.pone.0234992.g026

3.1.3 Daily field data of solar irradiation and temperature in Islamabad.

The wind speed (m/s), ambient temperature (°C), and solar irradiation (W/m2) are obtained from the Pakistan Meteorological Department (PMD) for a complete solar day at Islamabad station. In this simulation, each hour is modeled for one second of simulation time. The case study taken is the Defense Housing Authority (DHA), Islamabad, Pakistan. The irradiation varies with the appearance of the sun. The average irradiation level during day time is about 1000 W/m2, while the average temperature level is 20°C having a maximum peak of about 42.6°C during day time. Fig 27 shows the irradiance profile and Fig 28 shows the temperature profile used for this case study.

Fig 29 shows the complete power profile of all the sources and total load connected in the microgrid system during one day. Fig 30 shows the comparison of the power error of all the proposed controllers with aPID. The sudden spikes in the results are because of sudden changes in the irradiance level under PSCs that causes a control error. This initiates the controller which adjusts the duty cycle to minimize the power error.

Fig 30 also shows that the maximum power error generated by the least performing controller adaptive FBL embedded FRANF-SAM is not exceeding 31.3 kW at its highest peak, whereas the maximum power error generated by the best performing controller i.e. Adaptive FBL embedded FRANF-MHW is below 21.9 kW at its highest peak. However, aPID produced a maximum power error of 58kW at its highest peak. The maximum peak and average power error of all the control schemes starting from the least performance to the best is shown in Table 6.

thumbnail
Table 6. Peak power error and average power error of all proposed controllers.

https://doi.org/10.1371/journal.pone.0234992.t006

The output power comparison of all proposed controllers and reference power is shown in Fig 31. For every sharp change in reference power, there is also a sharp spike in tracked power at the same instant of time. This proves the ability of proposed controllers in dealing with sharp sudden changes in the reference signal. The only difference in the performance of the proposed controllers is the net output generated power at the same instant of time. It can be seen clearly that adaptive FBL embedded FRANF-SAM and adaptive FBL embedded FRANF-CW are producing low output power as compared to adaptive FBL embedded FRANF-FS and Adaptive FBL embedded FRANF-MHW. This is because of the different consequent structures used in FRANF. The PPVFBMHW acquires the PV system output power with steady-state error = 0.3 kW, undershoot = -9.19% and overshoot = 4.39%. The PPVFBFS extracts the PV system output power with steady-state error = 11.5 kW, undershoot = -20.95% and overshoot = 2.48%. The PPVFBSAM obtained the PV system output power with steady state-error = 18.1 kW, undershoot = -26.99% and overshoot = 1.16%. The PPVFBCW obtained the PV system output power with steady state-error = 18.4 kW, undershoot = -28.33% and overshoot = 1.49%.

The dynamic efficiency ηPV of various controllers on the basis of their tracked powers and power error are calculated as: (58) where; PPV = VPV × IPV, to = 0 h and tf = 24 h are initial and final intervals respectively. The dynamic efficiency, ηPV of the proposed controllers are shown in Fig 32. The peak ηPV of aPID is 75.56%, adaptive FBL embedded FRANF-SAM is 87.1%, adaptive FBL embedded FRANF-CW is 86.4%, adaptive FBL embedded FRANF-FS is 93.1% and adaptive FBL embedded FRANF-MHW is 95.9%. The efficiency plot of all the controllers shows a small variation. This is due to continuous change in the reference signal according to the solar irradiation at that instant of time and shows the high sensitivity of the proposed controllers. The peak and average efficiency of all controllers starting from the least to the most efficient are shown in Table 7.

thumbnail
Table 7. Peak and average efficiency of all proposed controllers.

https://doi.org/10.1371/journal.pone.0234992.t007

Performance indexes including Integral Square Error (ISE), Integral Time Square Error (ITSE), Integral Absolute Error (IAE), and Integral Time Absolute Error (ITAE), calculated based on Perror in Eq (57) as shown in Figs 33, 34, 35 and 36. The comparison of performance indexes plots shows that the accumulative error in all schemes increases with time. Again it is clear that the index of adaptive FBL embedded FRANF-MHW is least among all proposed controllers. Table 8 shows the values of various indexes of all the proposed controllers in contrast with aPID showing least to most performing from top to bottom.

thumbnail
Table 8. Performance indexes of all proposed controllers.

https://doi.org/10.1371/journal.pone.0234992.t008

Fig 37 shows a comparison of the percentage change in total harmonic distortion (THD) for load current due to individual control scheme. The percentage change in THD complies with the IEEE standard 1547 [45]. The result shows that the percentage change in THD due to the adaptive FBL embedded FRANF-MHW control scheme is the smallest of all and proves its better performance among other proposed controllers. The percentage change in frequency is also shown in Fig 38 for all proposed controllers. It can be observed that the percentage change in frequency due to the adaptive FBL embedded FRANF-MHW control scheme is almost flat and nearly zero. The percentage change in frequency for all the proposed controllers is within the acceptable range according to the IEEE standard 1547 [45].

Another factor of comparison is the percentage change in voltage (VRMS) produced in the AC-bus during the conversion process. Fig 39 shows the performance of the individual proposed scheme. It can be seen that the percentage change in VRMS due to adaptive FBL embedded FRANF-MHW control scheme has the least value among all proposed control scheme and thus ensures its superior performance. The overall performance of all the described controllers is compared based on power quality, average and peak power error in MPPT, efficiency, performance indexes, percentage THD change, percentage frequency change, and percentage VRMS in load current. For daily field data, a spider chart is plotted for all comparable performance parameters in Fig 40. Statistical values are scaled to a specific level for each individual parameter to obtain a clear view. The graphical view of IAE, ITAE, ISE, ITSE, maximum peak power error, and average power error in Fig 40 shows that adaptive FBL embedded FRANF-MHW control has the least values for these parameters as compared to other proposed controllers. Whereas, maximum peak efficiency and average efficiency values of Adaptive FBL embedded FRANF-MHW are the highest comparably. The spider chart shows the superior performance of the adaptive FBL embedded FRANF-MHW control scheme over other proposed controllers in a single glance.

Conclusions

In this article, four intelligent control schemes are presented for the MPPT problem of a grid-connected PV subsystem in an SMG-HPS for three different solar and temperature profiles. Results of all the control schemes are compared against one another and aPID for various parameters obtained through simulations. The performance of all the proposed control schemes is within the acceptable range and simply cannot be rejected at any ground. However, the overall analysis shows the performance of the adaptive FBL embedded FRANF-MHW is superior to all the other proposed control schemes. This is due to the use of a continuous signal wavelet that is the MHW in the consequent part of the identifier as well as the adaptive recurrent weights in antecedent and consequent parts of the same identifier. On the other hand, the CW is a discrete composite mathematical function and the recurrent weight of the consequent part of this and other proposed identifier schemes are fixed gains that are not being updated during the simulation.

Future studies include testing of adaptive FBL embedded FRANF hybrid wavelet control compared to those presented in this article.

References

  1. 1. Dolan DSL, Durago J, Taufik. Development of a photovoltaic panel emulator using Labview. In: 2011 37th IEEE Photovoltaic Specialists Conference; 2011. p. 001795–001800.
  2. 2. Khan M, Mathew L. Different Kinds of Maximum Power Point Tracking Control Method for Photovoltaic Systems: A Review. Archives of Computational Methods in Engineering. 2016;24.
  3. 3. Khan M, Shukla P, Mustafa R, Chatterji S, Mathew L. Different types of maximum power point tracking techniques for renewable energy systems: A survey. vol. 1715; 2016. p. 020015.
  4. 4. Subudhi B, Pradhan R. A Comparative Study on Maximum Power Point Tracking Techniques for Photovoltaic Power Systems. IEEE Transactions on Sustainable Energy. 2013;4(1):89–98.
  5. 5. Lalili D, Mellit A, Lourci N, Medjahed B, Berkouk EM. Input output feedback linearization control and variable step size MPPT algorithm of a grid-connected photovoltaic inverter. Renewable Energy. 2011;36(12):3282–3291. https://doi.org/10.1016/j.renene.2011.04.027.
  6. 6. Mellit A, Massi Pavan A. Performance prediction of 20kWp grid-connected photovoltaic plant at Trieste (Italy) using artificial neural network. Energy Conversion and Management. 2010;51:2431–2441.
  7. 7. Xiao W, Ozog N, Dunford WG. Topology Study of Photovoltaic Interface for Maximum Power Point Tracking. IEEE Transactions on Industrial Electronics. 2007;54(3):1696–1704.
  8. 8. Hassan SZ, Li H, Kamal T, Arifoaylu U, Mumtaz S, Khan L. Neuro-Fuzzy Wavelet Based Adaptive MPPT Algorithm for Photovoltaic Systems. Energies. 2017;10(3).
  9. 9. Mumtaz S, Ahmad S, Khan L, Ali S, Kamal T, Hassan SZ. Adaptive Feedback Linearization Based NeuroFuzzy Maximum Power Point Tracking for a Photovoltaic System. Energies. 2018;11(3).
  10. 10. IJMTST E. Adaptive Control Scheme for PV Based Induction Machine. 2016.
  11. 11. Kamran M, Mudassar M, Fazal MR, Asghar MU, Bilal M, Asghar R. Implementation of improved Perturb & Observe MPPT technique with confined search space for standalone photovoltaic system. Journal of King Saud University–Engineering Sciences. 2018; https://doi.org/10.1016/j.jksues.2018.04.006.
  12. 12. Abdourraziq MA, Ouassaid M, Maaroufi M, Abdourraziq S. Modified P O MPPT technique for photovoltaic systems. In: 2013 International Conference on Renewable Energy Research and Applications (ICRERA); 2013. p. 728–733.
  13. 13. Abdelsalam AK, Massoud AM, Ahmed S, Enjeti PN. High-Performance Adaptive Perturb and Observe MPPT Technique for Photovoltaic-Based Microgrids. IEEE Transactions on Power Electronics. 2011;26(4):1010–1021.
  14. 14. Ahmad J, Spertino F, Di Leo P, Ciocia A. A Variable Step Size Perturb and Observe Method Based MPPT for Partially Shaded Photovoltaic Arrays. In: PCIM Europe 2016; International Exhibition and Conference for Power Electronics, Intelligent Motion, Renewable Energy and Energy Management; 2016. p. 1–8.
  15. 15. Gokmen N, Karatepe E, Ugranli F, Silvestre S. Voltage band based global MPPT controller for photovoltaic systems. Solar Energy. 2013;98:322–334. https://doi.org/10.1016/j.solener.2013.09.025.
  16. 16. Bahari MI, Tarassodi P, Naeini YM, Khalilabad AK, Shirazi P. Modeling and simulation of hill climbing MPPT algorithm for photovoltaic application. In: 2016 International Symposium on Power Electronics, Electrical Drives, Automation and Motion (SPEEDAM); 2016. p. 1041–1044.
  17. 17. Putri RI, Wibowo S, Rifaaei M. Maximum Power Point Tracking for Photovoltaic Using Incremental Conductance Method. Energy Procedia. 2015;68:22–30. https://doi.org/10.1016/j.egypro.2015.03.228.
  18. 18. Trivedi A, Gupta A, Pachauri RK, Chauhan YK. Comparison of Perturb amp; Observe and Ripple correlation control MPPT algorithms for PV array. In: 2016 IEEE 1st International Conference on Power Electronics, Intelligent Control and Energy Systems (ICPEICES); 2016. p. 1–5.
  19. 19. Messalti S, Harrag AG, Loukriz AE. A new neural networks MPPT controller for PV systems. In: IREC2015 The Sixth International Renewable Energy Congress; 2015. p. 1–6.
  20. 20. Chekired F, Larbes C, Rekioua D, Haddad F. Implementation of a MPPT fuzzy controller for photovoltaic systems on FPGA circuit. Energy Procedia. 2011;6:541–549. https://doi.org/10.1016/j.egypro.2011.05.062.
  21. 21. Yang B, Yu T, Shu H, Qiu D, Zhang Y, Cao P, et al. Passivity-based linear feedback control of permanent magnetic synchronous generator-based wind energy conversion system: design and analysis. IET Renewable Power Generation. 2018;12(9):981–991.
  22. 22. Hassan S, Kamal T, Mumtaz S, Khan L. In: An Online Self Recurrent Direct Adaptive Neuro-Fuzzy Wavelet Based Control of Photovoltaic Systems; 2019. p. 233–250.
  23. 23. Mumtaz S, Khan L. Adaptive control paradigm for photovoltaic and solid oxide fuel cell in a grid-integrated hybrid renewable energy system. PLOS ONE. 2017;12(3):1–33.
  24. 24. Cambera JC, Feliu-Batlle V. Input-state feedback linearization control of a single-link flexible robot arm moving under gravity and joint friction. Robotics and Autonomous Systems. 2017;88:24–36. https://doi.org/10.1016/j.robot.2016.11.019.
  25. 25. Kim J, Croft EA. Full-State Tracking Control for Flexible Joint Robots With Singular Perturbation Techniques. IEEE Transactions on Control Systems Technology. 2019;27(1):63–73.
  26. 26. Chen Y, Damm G, Benchaib A, Lamnabhi-Lagarrigue F. Feedback linearization for the DC voltage control of a VSC-HVDC terminal. In: 2014 European Control Conference (ECC); 2014. p. 1999–2004.
  27. 27. Accetta A, Alonge F, Cirrincione M, Pucci M, Sferlazza A. Feedback Linearizing Control of Induction Motor Considering Magnetic Saturation Effects. IEEE Transactions on Industry Applications. 2016;52(6):4843–4854.
  28. 28. Chow TWS, Yong Fang. A recurrent neural-network-based real-time learning control strategy applying to nonlinear systems with unknown dynamics. IEEE Transactions on Industrial Electronics. 1998;45(1):151–161.
  29. 29. Shi Z, Shi M, Li C. The prediction of character based on recurrent neural network language model. In: 2017 IEEE/ACIS 16th International Conference on Computer and Information Science (ICIS); 2017. p. 613–616.
  30. 30. Williams RJ, Zipser D. A Learning Algorithm for Continually Running Fully Recurrent Neural Networks. Neural Computation. 1989;1(2):270–280.
  31. 31. Salem F, Awadallah MA. Detection and assessment of partial shading in photovoltaic arrays. Journal of Electrical Systems and Information Technology. 2016;3(1):23–32. https://doi.org/10.1016/j.jesit.2015.10.003.
  32. 32. Yang B, Yu T, Zhang X, Li H, Shu H, Sang Y, et al. Dynamic leader based collective intelligence for maximum power point tracking of PV systems affected by partial shading condition. Energy Conversion and Management. 2019;179:286–303. https://doi.org/10.1016/j.enconman.2018.10.074.
  33. 33. Yang B, Zhong L, Zhang X, Shu H, Yu T, Li H, et al. Novel bio-inspired memetic salp swarm algorithm and application to MPPT for PV systems considering partial shading condition. Journal of Cleaner Production. 2019;215:1203–1222. https://doi.org/10.1016/j.jclepro.2019.01.150.
  34. 34. Mohamed MA, Diab AAZ, Rezk H. Partial shading mitigation of PV systems via different meta-heuristic techniques. Renewable Energy. 2019;130:1159–1175. https://doi.org/10.1016/j.renene.2018.08.077.
  35. 35. Eltamaly AM, Farh HMH, Al-Saud MS. Grade point average assessment for metaheuristic GMPP techniques of partial shaded PV systems. IET Renewable Power Generation. 2019;13:1215–1231(16).
  36. 36. Mostafaee G, Ghandehari R. Power enhancement of photovoltaic arrays under partial shading conditions by a new dynamic reconfiguration method. Journal of Energy Management and Technology. 2020;4(1):46–51.
  37. 37. Gupta S, Saurabh K. Artificial mountain ape optimization algorithm for maximum power point tracking under partial shading condition. In: 2017 International Conference on Energy, Communication, Data Analytics and Soft Computing (ICECDS); 2017. p. 796–801.
  38. 38. Ahmad S, Khan L. Adaptive feedback linearization based HVDC damping control paradigm for power system stability enhancement. In: 2016 International Conference on Emerging Technologies (ICET); 2016. p. 1–6.
  39. 39. Werbos PJ. Generalization of backpropagation with application to a recurrent gas market model. Neural Networks. 1988;1(4):339–356. https://doi.org/10.1016/0893-6080(88)90007-X.
  40. 40. Tutschku K. Recurrent Multilayer Perceptions for Identification and Control: The Road to Applications; 1995.
  41. 41. W Sun YW. A recurrent fuzzy neural network based adaptive control and its application on robotic tracking control. Neural Information Processing–Letters and Reviews. 2004;5:19–26.
  42. 42. Song J, Shi H. Dynamic system modeling based on wavelet recurrent fuzzy neural network. In: 2011 Seventh International Conference on Natural Computation. vol. 2; 2011. p. 766–770.
  43. 43. Liu Y, Lin Y, Wu S, Chuang C, Lin C. Brain Dynamics in Predicting Driving Fatigue Using a Recurrent Self-Evolving Fuzzy Neural Network. IEEE Transactions on Neural Networks and Learning Systems. 2016;27(2):347–360. pmid:26595929
  44. 44. Pratama M, Lu J, Lughofer E, Zhang G, Er MJ. An Incremental Learning of Concept Drifts Using Evolving Type-2 Recurrent Fuzzy Neural Networks. IEEE Transactions on Fuzzy Systems. 2017;25(5):1175–1192.
  45. 45. IEEE Standard for Interconnection and Interoperability of Distributed Energy Resources with Associated Electric Power Systems Interfaces. IEEE Std 1547-2018 (Revision of IEEE Std 1547-2003). 2018; p. 1–138.