A Novel Hybrid MPPT Technique to Maximize Power Harvesting from PV System under Partial and Complex Partial Shading

: Photovoltaic (PV) system has been extensively used over the last few years because it is a noise-free, clean, and environmentally friendly source of energy. Maximum Power Point (MPP) from the PV energy systems is a challenging task under modules mismatching and partial shading. Up till now, various MPP tracking algorithms have been used for solar PV energy systems. Classical algorithms are simple, fast, and useful in quick tracing the MPP, but restricted to uniform weather conditions. Moreover, these algorithms do not search the Global Maxima (GM) and get stuck on Local Maxima (LM). However, bio-inspired algorithms help ﬁnd the GM but their main drawback is that they take more time to track the GM. This paper addresses the issue by using the combination of conventional Incremental Conductance (InC) with variable step size and bio-inspired Dragonﬂy Optimization (DFO) algorithms leading to a hybrid (InC-DFO) technique under multiple weather conditions, for instance, Uniform Irradiance (UI), Partial Shading (PS), and Complex Partial Shading (CPS). To check the robustness of the proposed algorithm, a comparative analysis is done with six already implemented techniques. The results indicate that the proposed technique is simple, efﬁcient with a quicker power tracking capability. Furthermore, it reduces undesired oscillation around the MPP especially, under PS and CPS conditions. The proposed algorithm has the highest efﬁciencies of 99.93%, 99.88%, 99.92%, and 99.98% for UI, PS1, PS2, and CPS accordingly among all techniques. It has also reduced the settling time of 0.75 s even in the case of the CPS condition. The performance of the suggested method is also veriﬁed using real-time data from the Beijing database.


Introduction
The demand for Electrical energy is increasing over time. Thus, the consumption of fossil fuels is also growing rapidly and their results greatly affect the environment, for instance, in terms of fossil fuel produces Carbon and Sulphur dioxide which causes the greenhouse effect. Presently, Renewable Energy Resources (RES) are playing a significant role in substituting the traditional energy resources of fossil fuel. Thus, many countries are switching to RES e.g., Iceland, Sweden, the U.S., Costa Rica, and Nicaragua. In the previous decade among all RES, solar energy has emerged as a promising source of energy and has gained more attention because it is a noise-free, pollution-free, increasingly cost-effective, and limitless source of energy. In the coming years based on recent advancements in solar Mechanism is incorporated. In this way, the unneeded motion of the particle is reduced causing fast convergence. An initial form of the SSJ proposed algorithm is evaluated in comparison with widely known MPPT algorithms i.e., modified InC, hybrid PSO, and modified CS. The effectiveness of these techniques are analyzed based on searching several required computation cycles, global maximum power point correctly, tracking time, and the scanned interval. Results demonstrate that the proposed technique significantly increases the accuracy and speed of finding the real global maximum power point. Moreover, it is irrespective of the complexity of PS conditions [10].
In [11], the authors describe the analysis of PS conditions for PV energy systems. In a PV energy system during partial shading conditions, if the PV modules increase in a string, then the 0.8 V oc model is used to raise the efficiency. Peaks in the case of PS are situated at the multiples of 0.8 V oc . For finding the correct position of peaks in the PV energy system for partial shading conditions, the authors [11] used a set of new equations model. Problems of partial shading in long-term PV systems are resolved using a model. However, this model slows the speed of the overall PV system.
Using the three-point-based method, the authors in [12] examine the parameters extraction of a single diode model. The pivot point is used to make points of critical data. The model forces to pass through these points after that Evolutionary Algorithm (EA) is applied to fine-tune the remaining section of the suitable curve. The EA technique shows substantial enhancement for fitting accuracy. From the MATLAB/Simulink results, it is concluded that the suggested method constantly converges to an optimum solution in less than 50 iterations but increases the system complexity [12].
The authors in [7] also proposed many techniques based on artificial intelligence like neural networks & fuzzy logic control. In the Artificial Neural Network (ANN), the input combination is considered as the surrounding condition of PV parameters, and output is assumed to be either reference voltage or duty cycle. As compared to the conventional neural network (NN) technique for MPPT, a multi-level neuro-fuzzy model has been proposed. The proposed model has been proven better efficient under a wide range of operating circumstances [13]. However, as the characteristic of PV array changes, the system needs to train periodically to obtain high accuracy [14].
The main benefits of Fuzzy Logic Control (FLC) is the acceptance of the noisy input signal and its ability to converge during a sudden change in surrounding conditions [15,16]. However, this technique determined the error based on a rule base table. The drawback of this method is that it highly depends on the availability of previous knowledge of the input signal [7]. In the literature, the Takagi-Sugeno model is used for further improvement in the stability and robustness of the system [17,18].
State-Space (SS) based MPPT technique uses the SS model to represent PV systems and dynamic feedback controllers to track MPP. Due to the ease of change in system parameters, this technique is very suitable even in fast-changing weather conditions and partially shaded conditions. However, no experimental verification is given [19,20]. The Slide Control Method (SCM) is achieving MPP that achieved the control of DC/DC converters. This technique normally employs the hysteric mode control of the DC/DC converter. The leading advantage of this method is its dynamic response and fast speed. Although, the unstable switching frequency is its major drawback [21].
Another MPP Flower Pollination Algorithm (FPA) has been introduced by the researchers in which five pollens or duty cycles are initially generated with an FPA between 0-1 and measure the power of each pollen. The pollen that gives the maximum power is chosen as the best power. To deciding which pollen will pass through global pollination or local pollination, a new number between 0-1 is randomly generated that has a 0.8 switching probability "P" for achieving the power. Researchers concluded that the best power gives the global best value after 25 iterations only [22]. The major disadvantages of FPA are high computation time, difficulty to tune parameters, complex structure, and procedural complications [23].
Modified FPA is described for the MPP of the PV system. The performance of the system is improved but the structure of FPA has become more complex [24]. Anderson et al. [25], examined the MPP using global optimization issues known as Random Search Method (RSM) as it tracks the MPP using a random solution. In defined limits, Anderson et al. produce the solution randomly and search for appropriate solutions that have a range between 0-1. To measure the power, these values send to the DC-DC converter against each input. Furthermore, these values are restructured continuously until the MPP is obtained after each iteration.
In [26] a Modified Firefly algorithm (FA) for the calculation of MPP is elaborated. To decrease the randomization of fireflies, randomization α and β parameters are introduced by varying irradiance and temperature. From the result, the authors of [26] concluded that efficiency and tracking speed is improved as compared to the conventional FA.
The bypass diodes are reverse biased and have little effect under UI levels. Whenever the PV array is partially shaded, the bypass diode is forward biased, and the current flows via the diode instead of through the PV module. Many MPP's occur in the P-V curve under partial shading due to the existence of bypass diodes. Traditional MPPT algorithms are looking for the peak of the P-V curve centered on a localized initial stage. As a result, they are unable to distinguish between LMPP and GMPP [27,28].
The authors in [29] improved the modified FA in terms of searching speed and decreases the number of computations of conventional FA. A Simplified FA is elaborated in which decreases the computational load and rises the MPPT in the PV energy systems. Steady-state oscillations also reduce additionally. Many other optimization techniques for MPP in PV systems such as Artificial Bee Colony (ABC) [30], Ant Colony Optimization (ACO) [31], Bacterial Fragmentation and aging (BFA) [32], Genetic Algorithm (GA) [33], Differential Evaluation (DE) [34], Grey Wolf Optimization (GWO) [35] and Simulated Annealing (SA) [36] algorithms have been discussed and evaluated based on response time, accuracy, efficiency, partial shading, and performance of the PV system under rapid response.
Up till now, under all circumstances, no technique gave the best results in terms of convergence and settling time, power tracking capability, and efficiency. There is still exists a research gap that how the PV system gives the response and addresses the PS condition according to the system requirements. These gaps still required to be eliminated by hybridizing the two techniques. So, We hybridize the conventional InC and softcomputing DFO techniques and use their advantages. If the uniform condition occurs then there is no benefit of using the soft computing techniques because they are slow with respect to change in any environmental condition. In the majority of cases, there is no existing PS condition and system response uniformly. However, if the PS does not occur in the PV system then InC is an ideal technique for use. In the world mostly P&O and InC techniques are used because they are fast and simple but cannot address the PS and CPS. On the other hand, a research gap exists in soft computing techniques where some tackle the PS condition and cannot respond to the CPS condition. Our contribution in this paper is as follows: The performance of the suggested method is also verified using real-time data from the Beijing database.
The rest of the paper is arranged in the following Sections. Section 2 described the partially shaded PV energy system. Sections 3 and 4 elaborate on the previous and proposed MPPT techniques for PV energy systems. Section 5 presented the simulation model. In Section 6 results and discussion are illustrated for all weather conditions. Finally, the conclusion of the paper is given in Section 7.

Description of the Power and Control Circuits
Ideally, whenever light hits the photovoltaic cell, it behaves similarly to a current source and is connected with a parallel diode. A PV module's characteristics may be acquired through an ideal model of a single diode. Practically, a single diode model consists of parallel and series resistance named as R p and R s respectively.The schematic diagram of a single diode model of PV is demonstrated in Figure 1.
where the values of I d and I sh are represented below in (2) and (3): where reverse saturation current is denoted by I o . The value of V T is equal to diode's ideality factor is presented by N s , 'K' stands for Boltzmann constant and 'T' is expressed as the temperature of the photovoltaic cell. The complete governing equation of a PV cell and its output current (I) is shown by DC-DC converter in PV energy systems is an integral part of the output power. There are several converters commonly used in a PV system to modify the value of output power like a buck, boost, and buck-boost converter. They decrease or increase the value of voltages according to the system requirements. Among all the other DC to DC converters, the boost converter has the advantages in terms of voltages regulation, power, and efficiency to achieve the desired results. PV energy system's output voltages should be controlled, varying the value of duty cycle of a DC-DC boost converter. For achieving this, the output voltage and connected load must be synchronized. Figure 2 presents the basic schematic diagram of the DC to DC boost converter. Initially, the reference voltage (V re f r ) at GM of the PV system is adjusted then the value of V re f r of the boost converter is modulated to modify PV voltage as a reference voltage. The value of the duty cycle is taken from 0 to 1. The maximum energy can be harvested from the PV energy system by implementing the control algorithm in a microcontroller having zero oscillation. The mathematical equations of output capacitance, duty cycle, input capacitance, output voltage, and inductance of the boost converter are represented below in (6)- (10).
where V in is the input voltage, f switching represents the switching frequency, I ripple denotes the ripple current and DC elaborates the duty cycle of a DC-DC boost converter. In PV systems conventional techniques may be stuck on LM and do not achieve GM for maximum energy harvesting. Due to this reason, bio-inspired techniques are used in this research for searching the MPP.

Description of UI, PS, and CPS under Study
Uniform irradiance case has the same irradiance values on the PV panels. In the partially shaded condition, it observes the various irradiance level. Due to the PS, a series of PV modules linked together may experience a low current and do not give high efficiency. Many PV modules are connected in series to increase voltages and it is called a string. Multiple strings are connected with a single inverter.
However, due to different irradiance levels, it may result in the non-uniform current that causes the mismatching and hot-spot effects. These effects can be reduced by connecting a bypass diode that gives an extra path to the flow of the high current of other PV modules. PV array of 4 × 1 of same and different irradiance levels is shown in Figure 3a,b respectively. PS faces the problem of MPPT that results in low efficiency of PV energy system. This problem is solved by conventional algorithms such as P&O and InC as their decisions are based on negative and positive gradients [5]. These techniques may not find any difference between a single GM and LM that operates on one point only. However, PS produces more than one peak on the curve, and in this way, by using these conventional algorithms the performance is compromised. Bio-inspired algorithms can locate GM successfully for PV energy systems. Whenever PV modules in large numbers are combined in series and give the output power, then CPS takes place [10]. Multiple close peaks are found on the P-V and I-V characteristics curves. The collection of these number of different peaks on power and voltage curves is known as a cluster. The CPS block diagram is demonstrated in Figure 4. There is a single GM in all clusters of CPS cases. This single maximum is also known as cluster head maxima. However, in a string, the number of partially shaded PV modules have an equal number of peaks obtained on a PV curve from the PV energy system. P-V characteristics curve of unique type is shown if the shading of the PV energy systems is broadly distributed. Moreover, this curve has multiple close peaks points. In CPS, finding GM is a big problem.

Incremental Conductance (InC) Technique
First, the InC method finds the slope of conduction then changes the duty cycle of the converter. This algorithm works for slope values. A duty cycle of converter rises only for positive slope values and declines for others. It will repeat the same process until found the MPP is on the curve [37] The MPP of the PV array is determined when the derivative of Power by voltage is equal to zero. Basic equations of InC that are used are given as follow: Equation (11) can be write down as: Therefore, the authors of [38] proposed a permitted error of approximately = 0.002, to eliminate the steady-state error. Due to the complexity, the convergence speed of the overall system decreases, and the cost of the system increase [39].

Dragon Fly Optimization (DFO) Technique
DFO is also a bio-inspired algorithm. Around the world, 3000 different species of dragonfly exist. Adult and nymph are the two milestones of a dragonfly's lifecycle. In life, DF spends their most of time in the nymph. During the metamorphism process, they become adults. Dragonfly examines small predators that hunt all the other species of insects. They have rare and unique swarming behavior. There are two main purposes of the DF swarm; the first is hunting and the second is migration. Migration and hunting are also known as dynamic and static swarms. In a static swarm, DF makes a small group over a small area and fly in back and forth motion to hunt the other species of flying prey e.g., mosquitoes and butterflies. Static swarms have two important characteristics named abrupt changes and local movements. On the other hand, in the dynamics swarm, a massive number of DF's make the group for migration from one place to another place. These are the two main behaviors for the inspiration of using the DF technique. These two behaviors are very similar to meta-heuristics two-phase such as exploitation and exploration. The primitive pattern of DFO is shown in Figure 5 in between the individuals. DF's motion in the swarm can be mathematically represented using the following five behaviors known as separation, alignment, cohesion, food, and enemy [40].

1.
Separation. In a static swarm, any DF does not collide with the other fly in the neighborhood. The separation S i of the i th individuals can be computed as S i = − ∑ N (k=1) (P − P k ). The present position of DF is denoted by P, P k represents the k th neighboring DF and N is the total number of individuals neighboring.

2.
Alignment. It indicates that the velocity of the other DF matches with the individuals in the same neighbourhood. Alignment represents with the symbol A i and it is . V k stands for velocity of the k th neighbouring individuals.

3.
Cohesion. It is the liability of individuals to move about the midpoint of the mass of neighborhoods. Cohesion represents as C i and its Equation is as C i = ∑ N (k=1) P k N − P. P denote the current position of the DF's individuals and P k describe the k th position of neighboring DF.

4.
Food. For the sake of survival, all the individuals move towards food. Attraction of food F i at position x is shown F i = P Food + P. The position of current food is represented by P Food & P stands for the current individual's position.

5.
Enemy. All the individual moves away from the enemy and their equation is given as E i = P Enemy + P. P Enemy shows the position of the current enemy and P presents the current individuals position. These five attributes affect the position of individuals DF's. Final upgraded position of the individuals DF's is calculated as Where the value of ∆P i is given in (16) below as follow: where s describes the weight of separation, S i shows the separation of i th DF, a represents the weight of alignment, A i stands for i th alignment of DF, c denotes the cohesion's weight, C i elaborates the i th cohesion of DF, f shows the weighting factor of food, F i describe the ith food of DF, e is the weighting factor of enemy and E th is the i th enemy of DF.

Proposed Hybrid (InC-DFO) Technique for PV Energy Systems
For the sake of making a new strong technique, the authors of this paper merge the powerful characteristics of InC with variable step size and DFO algorithm known as (InC-DFO) technique based on their features in a meta-heuristic context.
The proposed InC-DFO technique is a computational technique that has the advantages of both algorithms. These two techniques are sequentially executed to enhance the effectiveness of the optimization problem. The output of the InC technique with variable step size calculates the initialization of the DFO technique. However, the InC-DFO technique is applied for all weather conditions like UI, partial, and complex partial conditions for PV energy systems. InC-DFO algorithm is used in two stages, in the first stage, the InC technique with variable step size is used because it performs efficiently under UI values and their less convergence and settling time but this technique is trapped on LM and results in power loss.
However, practically it is not possible to achieve constant MPP so there is still some oscillation present around MPP. Thus, a satisfying compromise has to be made among the oscillations and the dynamics for MPPT with fixed step size in the conventional InC technique. The smaller step size can decrease the oscillations around the MPP with decelerating dynamics. On the other side, a larger step size gives the speedy tracking of MPP but it may results in steady-state unavoidable oscillations. For solving these challenges, the InC algorithm is used with variable step size is given as follows: where d(k − 1) and dk are the boost converter duty cycle at previous instant (k − 1) and at instant k. N represents the scaling factor to adjust the step size. Whereas dev × P , dev × V and dev × I are the output power, voltage and current of PV array respectively. It can be written as follow: The location of the DF's in global MPP searching of PV systems is the duty ratio at which dc to dc boost converter runs. It is a one-dimensional issue. The PV array's power output is the goal function to be optimized, and the DA is employed to maximize this power. The following procedures are followed to implement the DF algorithm in the second stage for global MPP tracking.

•
Step 1: The DF's are dispersed with a one search space within d min and d max that are the search space in contexts of duty cycle. The values of d min and d max in this article are selected to be and 2% & 98%, correspondingly. The population's value of the step P i is initialized. • Step 2: PV power of the DF's at each position (duty ratio) is calculated. The food supply and opponent positions are modified as a result of this. • Step 3: Modified the 's', 'a', 'c', 'f', and 'e' values. • Step 4: The Equations S i = − ∑ N (k=1) (P − P k ) and E i = P Enemy + P are used to compute 'S', 'A', 'C', 'F', and 'E' for each individual. • Step 5: The radius of the neighborhood has been modified. • Step 6: The DF's positions are changed with the help of equation P i = P i + ∆P i . If this equation produces a location outside of the d min and d max bounds, the flies are placed on the other side of the boundary. • Step 7: Exit the algorithm if the termination criteria are satisfied.

•
Step 8: If the perceived power changes, relaunch the tracking.
Any individual's neighborhood radius is expanded in this application, and the weights 'f', 'c', 'a', 's', and 'e' are lowered in accordance with the iteration calculation. This helps the swarm accomplish effective exploration in its early generations. The DF's are more impacted by each other in subsequent rounds as the algorithm approaches convergence. It would also ensure precise convergence and keep the DF's out of local maxima. The opponent and food sources are picked as their best and the worst solutions identified by the swarm up to that point. It enables the swarm to travel towards the best place it has identified thus far while avoiding less promising parts of the search field. The termination criteria for this task is set to be the circumstance in which the mean deviation of all DF's around the optimal location would be less than 5%. This number can be chosen by the user based on the necessary convergence precision. The pseudocode of DF algorithm [40] is shown in Algorithm 1. There are many factors in the suggested DF method used for hybridization that has to be fine-tuned. Various exploitative and exploitative behaviors may be created during optimization using separation 's', alignment 'a', cohesion 'c', food 'f', and enemy 'e' parameters. From these parameters, some of them have constant values and others have changing values according to the scenarios. DF neighbors are extremely essential, therefore each imitation DF is considered to have a neighborhood with a specific radius. In a dynamic swarm, DF tends to coordinate their flight while maintaining adequate separation and coherence. Therefore, in a static swarming, alignments are usually extremely low but cohesiveness is quite strong to assault prey.
As a result, while searching the search space, we give DF low cohesion and high alignment weights, whereas when exploiting the search space, we give DF's high cohesion and low alignment. The radius of neighborhoods is enhanced linearly with the number of iterations when transitioning from exploration to exploitation. The food supply and opponent are picked from among the best and worst solutions discovered thus far by the entire swarm. This results in convergence into promising portions of its optimal solution and divergence away from non-promising ones. Various exploratory or exploiting swarm behaviors can also be obtained with varying weight levels. For exploiting, the value of 'c' is kept large, but an is kept low because static swarming requires strong cohesion and low orientation. In a dynamic swarm, but on the other hand, 'c' is maintained modest while an is raised for better exploration qualities.
A heuristic method in MPP tracking aims to identify GM in its transient state with the least amount of time and reduce oscillations in the steady-state surrounding GM to decrease power loss. DF algorithm is used to improve MPP searching by considering these two features. In our approach, the duty cycle represents the search space, which accepts values between 0 and 1. The output duty cycle and the location of the DF are related. Within the search space, the solution is generated at random. Except for initiation, the method stages are repeated until the stopping condition is fulfilled or the maximum number of iterations is reached. The greatest value of PV power at GM is the fitness function.

Partial Shading Detection Method to Transfer from InC to DFO
To comprehend the PS detection mechanism of the InC-DFO technique, it is important to figure out the I-V curve relationship of a single diode PV cell. The complete P-V curves of the PV system are shown in Figure 6 that consist of UI, PS1, and PS2 conditions. For this aim, an equivalent circuit model is examined. We calculate the value of I pv from (5).
where 'K' stands for Boltzmann constant, 'T' is expressed as the temperature of the photovoltaic cell and G represents the irradiance. If the value of T is remaining constant and G is assumed to be varying then I pv rewritten as follow: I pv,STC is the photovoltaic current at standard testing condition i.e., 1000 W/m 2 at 25°C. From [41], taking into consideration that maximum current I mpp exists in the vicinity of 0.9I sc , Equation (22) The detection strategy of InC-DFO takes I mpp and I sc points on I-V curves of the Sun Power SPR-315E-WHT-D module and their specifications are given in Table 1. The irradiance values that are chosen from the I-V curve are demonstrated in Figure 7.  When values of G drops from 1000 W/m 2 to 700 W/m 2 , it is noticed that the ratio of I sc results as I sc I sc,STC = 4.202 6.14 = 0.684 ≈ 0.7 Respectively, the same outcomes could be drawn during the change in MPP current ratio when G falls from 1000 W/m 2 to 700 W/m 2 , i.e., It is observed from the (24), that the measured values of irradiance at I mpp and I sc are approximately equal. In [41], I mpp conceals near the 0.8 open-circuit voltage of a selected array V oc−array and I sc is in the region of 0.8 open-circuit voltage V oc of a module. Consequently, I mpp and I sc may be denoted as I 0.8Voc−array and I 0.8Voc , respectively. So, (24) could be written using these terms: MPPT records the corresponding currents as 0.8Voc − array and 0.8Voc when they take a sample at 0.8Voc − array and 0.8Voc respectively. Values of G are determined at 0.8Voc − array and 0.8Voc is specified as G-1 and G-2, respectively.
In UI condition, values of G-1 and G-2 are almost identical due to the slope current between I 0.8Voc , and I 0.8Voc−array is very small. This phenomenon is clarified in Figure 7 through the I-V characteristics curve. At 3 different values of G i.e., 1000 W/m 2 , 700 W/m 2 and 300 W/m 2 of PV array, corresponding currents are labeled at I 0.8Voc , and I 0.8Voc−array . So, the values of G-1 and G-2 are calculated using (27) at I 0.8Voc and I 0.8Voc−array sequentially as follows: It is concluded from (28) and (29) that G-1 and G-2 values are near to the real value of G like 1000 W/m 2 . Though, (27) is estimated on the basis such that I mpp is situated at 0.9 I sc . This consideration is not mandatory correct for all values of G because of creating a mismatch value. In above discussed computation, mismatch value is 23 W/m 2 among G-1 and G-2 at 1000 W/m 2 at 1000 W/m 2 . However, G mismatch value at 700 W/m 2 between G-1 and G-2 is reduced to 12 W/m 2 , whereas at 300 W/m 2 G mismatch value within G-1 and G-2 is decreased closely to zero.
Based on the mismatch values of G, an important conclusion could be constructed that in stage one, InC techniques with variable step size, act as UI conditions for a below specific threshold value of irradiance. This threshold value 40 is calculated by authors [42] on I-V curves using an extensive survey of ten poly and mono-crystalline different PV modules. Therefore, it could be deduced that if the value of mismatch is greater than 40 W/m 2 then the PS condition is prevailing and it acts as a UI condition for less than 40 W/m 2 .
To overcome the problem of power loss, the bio-inspired DFO technique is used in the second stage for searching the GM. InC technique with variable step size initiate in stage one and force the PV energy system to converge to 80 percent of I sc . once the InC algorithm achieved the I sc then it triggers the DFO algorithm which settles the PV array to exact MPP. Thus, InC-DFO not only removes the cons of both methods but also searches the MPP at PS and CPS conditions. The flowchart of the proposed InC-DFO algorithm is given in Figure 8.  Figure 9 illustrates the MPPT methods' implementation in a commonly used PV system design. The PV model presented by [43] is used to build the entire arrangement. Four modules are linked in series to form a PV array. P mpp = 200 W, V mpp = 26 V, V oc = 32.9 V, I mp = 7.6 A, I sc = 8.21 are the parameters for the PV module. The implementation of the developed Hybrid(InC-DFO) based MPPT is assessed using a PV array that is exposed to a variety of partial shade patterns. The PV array and the battery load are connected via a boost converter. The parameters of the DC to DC converter are C in = 10 uF, L = 200 uH, C out = 47 uF, and F sw = 20 kHz, as computed in [44]. As illustrated in Figure 9, the MPPT methods have been performed in embedded blocks. This block accepts I pv and V pv as inputs and outputs the duty cycle of a DC to DC boost converter.

Results and Discussion
Four cases are described in this section for all-weather conditions e.g., UI, PS, and CPS conditions. From four cases, one is taken for UI, two for PS, and one for CPS. After different results and trails, 4 distant cases in terms of different irradiance values are chosen for performance calculation of nature-inspired algorithms. A model of PV energy system has been developed on the MATLAB/Simulink software. In a PV energy system model, the rating of a single PV module is given in Table 1.
All previous algorithms like InC, CS, PSO, DFO, FA, and ACO, and the proposed MPPT Hybrid (InC-DFO) technique are tested on the developed PV model. Four cases are used for testing all MPPT techniques for PV energy systems in which UI condition is illustrated in case 1. PS condition is described in cases 2 and 3. CPS condition expresses in case 4. On the basis of settling and tracking time, current, transients of voltage, power, efficiency, and oscillation, results are analyzed. All four cases and their results are discussed as follows.

Case 1. Uniform Irradiance Condition
All the modules of the PV system received the same irradiance level in this case and it changes with respect to time. The irradiance pattern is described in Table 2. Figure 7 presents the PV energy systems UI curves of P-V and I-V. The pattern of uniform irradiance that is taken is elaborated in Table 2.
When comparing the power of different techniques, it is observed that the power at GM is 864 W and the highest power tracked at 700( KW m 2 ) irradiance level is 863.4 W in case of uniform irradiance condition for proposed hybrid (InC-DFO) technique. While comparing the power of all techniques, FA has the lowest power 852 W. The tracked power in the sequence is as follows: Hybrid > DFO > InC > CS > ACO > PSO > FA. CS algorithm has more oscillation and the highest setting time of 1.94 s. Among all the techniques, FA tracked the least power of 852 W. The proposed hybrid technique saves power loss by reducing the oscillation and enhancing the efficiency of 99.93% of PV energy systems. CS has the highest oscillations that cause power loss when compared to other techniques.

Case 2. Partial Shading Condition
In case 2 partial shading is discussed and their irradiance pattern is presented in Table 2. Partial shading occurs because of not receiving equal irradiance values by PV arrays. Bypass diodes are used for less-performing diodes due to PS. There are multiple peaks that appear on the I-V and the P-V characteristics curve that can result in complications for searching the peaks. Furthermore, multiple peaks are characterized as a number of LM and single GM. Conventional techniques cannot differentiate between LM and GM that may result in the reduction of power loss. The I-V and P-V characteristic curves of PS cases 2 and 3 are described in Figure 10. There are three LM and single GM observed on the I-V and P-V characteristics curve. The value of GM for PS1 is 512 W.

Case 3. Partial Shading Condition
Case three also represents another PS condition and their irradiance pattern is illustrated in Table 2. The value of GM is tracked in this case of PS. The Value of GM tracked is 781 W. In case 2 of partial shading, a 4 × 1 PV array is used. The power of the proposed hybrid technique is compared with other DFO, PSO, CS, FA ACO, and InC techniques. Among all the techniques PSO has the highest oscillation and ACO has no oscillation but they cannot find the GM and are stuck on LM. Power loss occurs due to oscillation in the PV system. It is noted that the proposed technique presents the lowest settling time that decreases the power losses and enhances the efficiency of the PV energy systems. CS has the highest settling time that may result in loss of power.

Case 4. Complex Partial Shading Condition
Condition of CPS presents in case 4 and their irradiance pattern is illustrated in Table 2. Figure 12 shows the P-V & I-V characteristics curve of CPS phenomena. The proposed method can remove the oscillation of the PSO and CS algorithm, thus increasing the efficiency of the PV energy systems. However, InC gives the lowest efficiency 87.52% because it is not able to find GM and is stuck on any value of LM.
In the case of CPS, the hybrid algorithm achieved the highest efficiency 99.98%. In the first stage, the InC technique is used because of less settling time but it is stuck in local maxima and cannot find the GM. So, in the second stage, we use the DFO that has the highest capability of tracking power. The proposed hybrid technique has zero oscillation and the highest convergence time that results in enhancing the efficiency of the PV system. Power, Current, and Voltage comparison of all techniques for Case 3 and Case 4 are represented in Figure 13. To validate the effectiveness of the proposed (InC-DFO) and its comparison with six other MPPT techniques, four different cases of PV array are taken for all weather conditions. Among four cases, one is for UI condition, two cases for partial shading condition (PSC), and one is for CPS condition. On the basis of these four cases, the performance valuation of the PV system is carried out and their numerical analysis is illustrated in Table 3.  Table 3, it is observed that it has the lowest convergence and settling time. The proposed algorithm tracks the GM efficiently and gives maximum efficiency of 99.93%, 99.88%, 99.92%, and 99.98% in all operating conditions due to their two-stage running phenomena.

Case 5. Field Atmospheric Data of Beijing
Several factors are used to demonstrate the PV system's feasibility. Solar PV systems are commonly used in open regions where they are exposed to adverse weather conditions. Case 5 is all about a field study on a Chinese city, namely Beijing. Field atmospheric data are used to assess the proposed InC-DFO MPP searching method in real-world situations. Seasonal temperatures and irradiation levels are available from the Chinese meteorological agency (CMA) [46]. Total direct sunlight temperature and irradiance of two seasons are considered throughout the course of a day. The PV system consists of 4 × 1 arrays, each rated at 315.072 watts, for an installed capacity of 1260.288 watts. The data provides information on weather factors such as temperature and irradiance of the PV system. Transients of The 24-h have a 10-min resolution. Maximum power obtained during a 24-h period is used to determine the efficiency of the MPP tracking InC-DFO algorithm. Table 4 shows the results of Case 5. Summer is probably the most productive season in Beijing, having peak production power of 897.2 W between 10 a.m. and 3 p.m. In the fall, the highest hour of average power is 811 W. Spring and winter are the least useful seasons in terms of power achievement. Beijing is located at 39.9042 North latitude and 116.4074 East latitude. According to the findings in all-weather situations, the InC-DFO algorithm gathers the most power. Figures 14 and 15 demonstrate a comprehensive examination of the InC-DFO algorithm collected power with, ACO, FA, DFO, InC, PSO, and CS in Beijing over the course of 24 h in the summer and spring. The results demonstrate that the proposed InC-DFO algorithm delivers better results and is more efficient.  for all weather conditions like UI, PS, and CPS. This proposed technique saved the power losses because of the fast searching ability of a GM and reduced the settling time that may enhance the performance of the PV energy systems in a transient state. Moreover, ripples and overshoot is eliminated because of zero oscillation in the transient state. It is thus concluded that although the conventional techniques are suitable for uniform irradiance, they do not find GM in the case of partial shading. Similarly, they require less time for tracking but they are stuck on LM. On the other side, bio-inspired methods are complex and take more time in finding the MPP. The proposed hybrid technique is the combination of InC and DFO, and it is used in two stages. In the first stage, the conventional InC technique with variable step size is adopted because of less convergence and settling time while in the second stage, a bio-inspired DFO technique is implemented for finding the GM quicker. The proposed technique outperforms the competition in terms of MPP monitoring and integration at the global MPP. The suggested MPPT technique efficiently achieves tracking time improvements of 10-29% and spontaneous oscillation reductions of more than 91%. The hybrid MPPT technique provides faster maximum power point (MPP) monitoring, lower computational pressure, and increased performance.