Group Teaching Optimization Algorithm Based MPPT Control of PV Systems under Partial Shading and Complex Partial Shading

The most cost-effective electrical energy is produced by photovoltaic (PV) systems, and with the smallest carbon footprint, making it a sustainable renewable energy. They provide an excellent alternative to the existing fossil fuel-based energy systems, while providing 4% of global electricity demand. PV system efficiency is significantly reduced by the intrinsic non-linear model, maximum power point (MPP), and partial shading (PS) effects. These two problems cause major power loss. To devise the maximum power point tracking (MPPT) control of the PV system, a novel group teaching optimization algorithm (GTOA) based controller is presented, which effectively deals with the PS and complex partial shading (CPS) conditions. Four case studies were employed that included fast-changing irradiance, PS, and CPS to test the robustness of the proposed MPPT technique. The performance of the GTOA was compared with the latest bio-inspired techniques, i.e., dragon fly optimization (DFO), cuckoo search (CS), particle swarm optimization (PSO), particle swarm optimization gravitational search (PSOGS), and conventional perturb and observe (P&O). The GTOA tracked global MPP with the highest 99.9% efficiency, while maintaining the magnitude of the oscillation <0.5 W at global maxima (GM). Moreover, 13–35% faster tracking times, and 54% settling times were achieved, compared to existing techniques. Statistical analysis was carried out to validate the robustness and effectiveness of the GTOA. Comprehensive analytical and statistical analysis solidified the superior performance of the proposed GTOA based MPPT technique.


Introduction
Solar energy is utilized in commercial and standalone projects for electricity generation using concentrated solar and photovoltaic principles. Solar energy provides numerous merits, such as being In the literature many MPPT techniques have been proposed, and depending upon the accuracy, techniques vary significantly by the type of control strategies. Voltage, current, and duty cycle are utilized, along with a mathematical model of PV cells, to devise fundamental techniques, and are utilized to predict optimized operating points on the P-V curve. In open circuit voltage (OCV), is used similarly, and is utilized to model MPPT control [10]. In this model mathematical approximations are utilized to predict the operating point, and no real-time monitoring of output power is involved, so the accuracy cannot be reliable. Moreover, the operating point is not reliable because the I-V and P-V curves show the dependence upon operating conditions, and hence efficiency is compromised [11].
The simplest and most famous classical MPPT technique is the P&O algorithm. In the perturb and observe (P&O) method, a predefined small change is introduced to perturb the system, consequently, the observed values are calculated and a power comparison yields the direction of change on the P-V curve. It is a simple and effective technique, and has been implemented in standalone and hybrid modes. However, like other techniques, it does also offer certain critical drawbacks, such as the oscillations caused by the differential increase in duty cycle. To minimize these oscillations, the step size has to be small, which consequently causes a dip in tracking time. In the literature many MPPT techniques have been proposed, and depending upon the accuracy, techniques vary significantly by the type of control strategies. Voltage, current, and duty cycle are utilized, along with a mathematical model of PV cells, to devise fundamental techniques, and are utilized to predict optimized operating points on the P-V curve. In open circuit voltage (OCV), V OC is used similarly, and I SC is utilized to model MPPT control [10]. In this model mathematical approximations are utilized to predict the operating point, and no real-time monitoring of output power is involved, so the accuracy cannot be reliable. Moreover, the operating point is not reliable because the I-V and P-V curves show the dependence upon operating conditions, and hence efficiency is compromised [11].
The simplest and most famous classical MPPT technique is the P&O algorithm. In the perturb and observe (P&O) method, a predefined small change is introduced to perturb the system, consequently, the observed values are calculated and a power comparison yields the direction of change on the P-V curve. It is a simple and effective technique, and has been implemented in standalone and hybrid modes. However, like other techniques, it does also offer certain critical drawbacks, such as the oscillations caused by the differential increase in duty cycle. To minimize these oscillations, the step size has to be small, which consequently causes a dip in tracking time. Hence, there is a tradeoff between tracking time and efficiency in this technique. Incremental conductance (IC) is yet another technique which utilizes the change in instantaneous conductance with respect to current operating point to decide the direction of the operating point. In both P&O and IC, if the power (P) is improved, the direction of increment is retained and vice versa [12]. These techniques are easy to implement and are low cost. Under rapidly changing ambient conditions their draw backs become prominent and they fail to track actual MPP, due to their lack in properties, in terms of adaptively and robustness. Another technique known as the hill climbing (HC) mechanism utilizes direct control action. It suffers from similar drawbacks. Work done to improve the step size of the aforementioned techniques, improved the performance and efficiency but increased the complexity, and required accurate sensors. A comprehensive mathematical model integration with operating PV conditions is cumbersome, moreover, a partial shading scenario, and intrinsic oscillations could not be minimized significantly [13].
In an instance when PV system does not receive equal levels of irradiance, PS happens. The mismatch of produced current and voltage activates the PS effect. Isolated under shaded modules of arrays, bypass diodes in the antiparallel configuration are introduced. Initiation of bypass diodes produces multiple local maxima in I-V and P-V characteristic curves. It is observed that the number of peaks on a P-V curve is equal to the number of partially shaded modules connected in a series string of arrays. Local maxima (LM) and global maxima (GM) are two categories of maxima. There is only one GM, and the point of operation is called global maximum power point (GMPP). Classical techniques cannot differentiate among multiple maxima, hence considerable power is lost [14][15][16].
To share information regarding the best position, the PSO based MPPT technique uses the social interaction of swarm particles to share information. This information sharing is carried out using a weighted vector, which influences the position update of swarm particles, and helps the swarm particles to converge to the global solution. To break the LM trap, random numbers are embedded in velocity and position update equations. Random numbers produce unwanted oscillations and affect the settling time. The particle swarm optimization gravitational search (PSOGS) incorporates the effective exploration of PSO, with the accurate exploitation of GS, and improvises the working of PSO significantly. Although improvement in steady-state oscillations is observed, settling time is not significantly minimized. GS integration with PSO increases the number of parameters to be fine-tuned, the complexity of the technique, and the computational burden. ACO uses pheromones to share the information, which evaporate with time. Individually every ant is a member of a set representing a probable solution, so a cautious selection mechanism is required. Scout ants are also utilized to explore the search space [29]. In CSA random values are assigned using Levy flight, which causes unwanted fluctuations in output [20]. Lower population size makes these shortcomings meaningfully higher. A higher population number is utilized, which takes a toll on computation power and the time required to track GM. The main qualities required from bio-inspired techniques are quick convergence, least random oscillations, GM detection, robustness, and high power convergence. Although CSA, GHO, and PSOGS show reasonable tracking ability, the logical flow of these techniques involves Electronics 2020, 9,1962 4 of 24 randomness, and is more complex than the ordinary PSO algorithm, which makes their real-world implementation challenging [30].
Hybrid MPPT control takes advantage of robust conversion of classical techniques, and the SI techniques can be based upon the physical laws of nature, such as the swarming behavior of communal species. To overcome the drawbacks of the existing techniques by improving the performance over the iterative time, a new bio-inspired meta-heuristic algorithm group teaching optimization algorithm (GTOA) [31] was utilized to implement direct control of PV systems for the MPPT problem. Its ability to track GM was phenomenal. A comprehensive statistical and analytical study validated the superior performance of the proposed GTOA based MPPT for various weather conditions. A wide-ranging study was made to authenticate the effectiveness and robustness of the GTOA in this study. The major contribution of the proposed GTOA based MPPT techniques are 1.
The introduction of a novel implementation of GTOA based MPPT technique for the partial shading problem of the PV system. The exploration and exploitation of search space are controlled by parameter "b" and "c" for maximizing the GM identification and quicker tracking of GM.

2.
GTOA is effective in the optimization of uni-modal, multi-modal, and complex mathematical problems. The proposed GTOA has increased the overall tracking efficiency.

3.
The best solution is attained in the mechanism and used to converge the population. Meanwhile, the rest of the swarm particles iteratively search the space, due to which performance is not compromised. A comprehensive analysis was done using extensive case studies, and specifically CPS was elaborated.

PV Model
To model the PV electrical characteristics of a solar cell, an equivalent single diode model is presented in Figure 2. The resistances added by parallel and series connections are integrated in the model. performance over the iterative time, a new bio-inspired meta-heuristic algorithm group teaching optimization algorithm (GTOA) [31] was utilized to implement direct control of PV systems for the MPPT problem. Its ability to track GM was phenomenal. A comprehensive statistical and analytical study validated the superior performance of the proposed GTOA based MPPT for various weather conditions. A wide-ranging study was made to authenticate the effectiveness and robustness of the GTOA in this study. The major contribution of the proposed GTOA based MPPT techniques are 1. The introduction of a novel implementation of GTOA based MPPT technique for the partial shading problem of the PV system. The exploration and exploitation of search space are controlled by parameter "b" and "c" for maximizing the GM identification and quicker tracking of GM. 2. GTOA is effective in the optimization of uni-modal, multi-modal, and complex mathematical problems. The proposed GTOA has increased the overall tracking efficiency. 3. The best solution is attained in the mechanism and used to converge the population. Meanwhile, the rest of the swarm particles iteratively search the space, due to which performance is not compromised. A comprehensive analysis was done using extensive case studies, and specifically CPS was elaborated.

PV Model
To model the PV electrical characteristics of a solar cell, an equivalent single diode model is presented in Figure 2. The resistances added by parallel and series connections are integrated in the model. ℎ is the current produced which is directly proportional to the incidence of light. The current through diode is and the current ( ) flowing towards the output is shown in Equation (1).
The practical model of PV cell is modeled by adding intrinsic resistances, i.e., ℎ is the series resistance, and ℎ is the parallel resistance. The mathematical modeling is shown in Equation (3). I pvh is the current produced which is directly proportional to the incidence of light. The current through diode is I o and the current (I) flowing towards the output is shown in Equation (1).
The practical model of PV cell is modeled by adding intrinsic resistances, i.e., R sh is the series resistance, and R ph is the parallel resistance. The mathematical modeling is shown in Equation (3). A "Sunpreme SNPM-GX-315" solar panel was utilized in this study, and the electrical characteristics are shown in Table 1.

Partial Shading Condition
Under static irradiance, if all the PV panels are connected in series, the knee point of the I-V and P-V curve is the GMPP, as shown in Figure 3a,b, respectively. Electronics 2020, 9, x FOR PEER REVIEW 5 of 26 A "Sunpreme SNPM-GX-315" solar panel was utilized in this study, and the electrical characteristics are shown in Table 1.

Partial Shading Condition
Under static irradiance, if all the PV panels are connected in series, the knee point of the I-V and P-V curve is the GMPP, as shown in Figure 3a,b, respectively.
When light of different irradiance falls on the series of connected modules, partial shading occurs. In the PS condition, the PV panels produce a non-uniform current, which causes hot spots and mismatching effects. The bypass diode appears as a viable solution to reduce these effects, which produce peaks with different magnitudes in I-V and P-V characteristic curves, as shown in Figure  3c,d. Figure 4 shows the PV panels under static and non-static irradiance, with bypass and blocking diodes. The number of peaks in I-V and P-V characteristic curves depends upon the number of panels receiving non-uniform irradiance. Due to PS, there is one GM and multiple LMs.  When light of different irradiance falls on the series of connected modules, partial shading occurs. In the PS condition, the PV panels produce a non-uniform current, which causes hot spots and mismatching effects. The bypass diode appears as a viable solution to reduce these effects, which produce peaks with different magnitudes in I-V and P-V characteristic curves, as shown in Figure 3c,d. Figure 4 shows the PV panels under static and non-static irradiance, with bypass and blocking diodes. The number of peaks in I-V and P-V characteristic curves depends upon the number of panels receiving non-uniform irradiance. Due to PS, there is one GM and multiple LMs.

MPPT Using GTOA
A novel group teaching optimization algorithm was used for MPPT. The methodology of the proposed GTOA is explained in Section 3.

Framework of GTOA
A novel group teaching optimization algorithm was used for MPPT. By simulating a group teaching mechanism, the knowledge of the whole class (c) can be improved, which is the basic idea behind the proposed technique, i.e., GTOA. To implement the GTOA tor the optimization, a simple group of teaching model was presented based on the following rules.
a. The only difference between students is the capability to be acquiescent of knowledge. The greater challenge for the teacher to formulate the teaching plan depends upon the differences in ability to accept knowledge. b. The quality of a good teacher is to pay additional attention to the students who have a poor capability to accept knowledge. c. By self-learning, or by interacting with classmates, a student can improve his knowledge during the free time. d. To improve the knowledge of students, a good teacher allocation method is helpful. Four phases are proposed in this model, which are represented in Figure 5.

MPPT Using GTOA
A novel group teaching optimization algorithm was used for MPPT. The methodology of the proposed GTOA is explained in Section 3.

Framework of GTOA
A novel group teaching optimization algorithm was used for MPPT. By simulating a group teaching mechanism, the knowledge of the whole class (c) can be improved, which is the basic idea behind the proposed technique, i.e., GTOA. To implement the GTOA tor the optimization, a simple group of teaching model was presented based on the following rules. a.
The only difference between students is the capability to be acquiescent of knowledge. The greater challenge for the teacher to formulate the teaching plan depends upon the differences in ability to accept knowledge. b.
The quality of a good teacher is to pay additional attention to the students who have a poor capability to accept knowledge. c.
By self-learning, or by interacting with classmates, a student can improve his knowledge during the free time. d.
To improve the knowledge of students, a good teacher allocation method is helpful. Four phases are proposed in this model, which are represented in Figure 5.

Ability Grouping
To represent the knowledge of the whole class, a normal distribution function was used which is shown in Equation (5).
where x is value for which normal distribution function is required. μ is the mean and δ is the standard deviation. In group teaching optimization, all students are divided into two groups. An outstanding group is a group of students with a good ability to grasp knowledge, whereas the group with the poor capability to grasp knowledge is called the average group.

Teacher Phase
In this phase, students learn from the teacher, which is the second rule previously defined. In the GTOA the teacher makes separate plans for average and outstanding groups.

Teacher phase I:
The teacher focuses on improving the knowledge of the whole class due to the good ability of students to accept knowledge. Students belonging to an outstanding group can increase knowledge by Equation (6).
where the number of students is represented by , is the knowledge of each student, is the teacher's knowledge, is the mean groups' knowledge. Teacher factor is represented by , ℎ _ is the knowledge of student j learning from the teacher. Randomness is introduced by , , and in the range (0,1).

Teacher phase II:
Due to weak knowledge accepting ability, based on the second rule, the teacher provides more attention to the average group. An average group of students can gain knowledge by Equation (9).

Ability Grouping
To represent the knowledge of the whole class, a normal distribution function was used which is shown in Equation (5).
where x is value for which normal distribution function is required. µ is the mean and δ is the standard deviation. In group teaching optimization, all students are divided into two groups. An outstanding group is a group of students with a good ability to grasp knowledge, whereas the group with the poor capability to grasp knowledge is called the average group.

Teacher Phase
In this phase, students learn from the teacher, which is the second rule previously defined. In the GTOA the teacher makes separate plans for average and outstanding groups.

Teacher phase I:
The teacher focuses on improving the knowledge of the whole class due to the good ability of students to accept knowledge. Students belonging to an outstanding group can increase knowledge by Equation (6).
where the number of students is represented by N, X j is the knowledge of each student, T is the teacher's knowledge, M is the mean groups' knowledge. Teacher factor is represented by F, X teacher_ j is the knowledge of student j learning from the teacher. Randomness is introduced by a, b, and c in the range (0,1).

Teacher phase II:
Due to weak knowledge accepting ability, based on the second rule, the teacher provides more attention to the average group. An average group of students can gain knowledge by Equation (9).
Electronics 2020, 9, 1962 8 of 24 where d is randomness (0,1). Equation (10) addresses the problem where a student cannot get knowledge from the teacher phase.

Student Phase
In spare time students can gain knowledge by self-learning, or by interacting with fellow students, which can be represented mathematically as shown in Equation (11). The student phase links to the third rule by including student phases I and II.
where e and g are two randoms (0,1), X t+1 student,j is knowledge of the student j, and X t+1 teacher,j is knowledge of the student j learning from the teacher. Students cannot gain knowledge in this phase. He or she can be addressed by Equation (12).

Teacher Allocation Phase
In order to improve the student's knowledge, a good teacher allocation method is crucial, and which is defined by the fourth rule. Inspired by the hunting behavior of grey wolves, the top three best students are selected, as presented in Equation (13).
where X t f irst , X t second , and X t third are the top three students, respectively.

Working of GTOA as MPPT
GTOA works on the basis of group learning, which increases the efficiency and reduces the oscillations, due to the mean value of knowledge derived from the top three students. CPS conditions is explained in Figure 6. The tracking mechanism of GTOA's particles in search space over the P-V curve is shown in Figure 7. Figure 7a shows the power tracked, Figure 7b is the P-V curve under PS conditions, and Figure 7c represents the transients in the reference voltage generated by GTOA. The Boost convert provides the necessary control and load interface given in Appendix A.
Partial shading occurs due to the difference in irradiance levels on PV modules connected in series. Due to bypass diodes, the modules connected in series create complex multiple peaks in the form of groups, called a cluster. In every cluster, the MPP is called the cluster head maxima (CHM). CPS conditions are shown in Figure 6. The left half-plane contains cluster 1, whose maxima is 1117 W. From left to right, the values of peaks in cluster 1 are 792.9 W, 958.3 W, 1028 W, 1117 W, and 1080 W. In the right half of the plane, cluster 2 contains GMPP whose value is 1136 W. The values of the local peaks in cluster 2 from left to right are 1136 W, 1128 W, 1035 W, and 828.4 W. The CHM in cluster 1 and GMPP has a difference of 19 W which is less than 2% of GMPP. Therefore, the re-initialization of the code will not occur.
Electronics 2020, 9, x FOR PEER REVIEW 9 of 26 learning through grouping, proper re-initialization of code, and initialization of population give effective results in CPS.      Due to the complex partial shading, around 6% of power is lost due to the incorrect detection of GM. Due to incorrect initialization and mistuning of the parameter in bio-inspired techniques, the exploring populations gets stuck in the local maxima, and cause power loss. In GTOA, effective learning through grouping, proper re-initialization of code, and initialization of population give effective results in CPS.

Case 1: Fast Varying Irradiance
Case 1 represents the fast varying irradiance pattern, to check the robustness of MPPT techniques in the tracking of mechanism power and re-tracking with changes in the weather conditions. After every 2 s, a change is produced in irradiance levels to reinitialize the MPPT technique to re-track the MPP. At 2 s irradiance level changed from 1000 W/m 2 to 450 W/m 2 , and again at 4 s it changed to 750 W/m 2 , to test fast varying conditions. Irradiance levels are shown in Table 2. Figures 9 and 10 show the power tracked, and zoomed-in view of the power tracked. In the first interval at STC, the MP is 1264 W, and the power achieved by the GTOA, DFO, PSOGS, PSO, CS, and P&O, was 1263. 5

Case 1: Fast Varying Irradiance
Case 1 represents the fast varying irradiance pattern, to check the robustness of MPPT techniques in the tracking of mechanism power and re-tracking with changes in the weather conditions. After every 2 s, a change is produced in irradiance levels to reinitialize the MPPT technique to re-track the MPP. At 2 s irradiance level changed from 1000 W/m 2 to 450 W/m 2 , and again at 4 s it changed to 750 W/m 2 , to test fast varying conditions. Irradiance levels are shown in Table 2.        Large oscillations can be observed in P&O around the GM, which cause power loss. In the first interval under STC, the efficiency achieved by GTOA, DFO, PSOGS, PSO, and CS was 99.96%, 99.84%, 99.68%, 99.92%, 99.92%, and 98.02%, respectively. PSO achieved higher efficiency compared to PSOGS, at the cost of higher fluctuations. To measure the performance of MPPT techniques, in terms of tracking power, the average power was taken.
In case 1 the average power was 921.3 W, and the average power tracked by GTO, DFO, PSOGS, PSO, CS, and P&O was 916.8 W, 905.8 W, 905 W, 903 W, 913 W, and 888 W, respectively. Re-initialization of particles under fast changing irradiance can be observed in Figure 11, in the duty cycle comparison at 2 s and 4 s. A zoomed-in duty cycle comparison can also be seen in Figure 12.

Case 2: Partial Shading
This case dealt with the PS condition, as mentioned in Figure 6, which contains multiple LMs and one GM. PSC was used to check the robustness of the MPPT technique. The irradiance pattern was 650 W/m 2 , 520 W/m 2 , 250 W/m 2 , and 400 W/m 2 . For case 2, MPP was at 427 W. Irradiance pattern is shown in Table 2.
Randomly selected positions of swarm particles in initialization are a better solution to track GM under PS conditions, but take more time to converge at GM. PSO and CS have more random oscillations due to velocity vectors. Another solution to this problem is the use of a large number of

Case 2: Partial Shading
This case dealt with the PS condition, as mentioned in Figure 6, which contains multiple LMs and one GM. PSC was used to check the robustness of the MPPT technique. The irradiance pattern was Electronics 2020, 9,1962 13 of 24 650 W/m 2 , 520 W/m 2 , 250 W/m 2 , and 400 W/m 2 . For case 2, MPP was at 427 W. Irradiance pattern is shown in Table 2.
Randomly selected positions of swarm particles in initialization are a better solution to track GM under PS conditions, but take more time to converge at GM. PSO and CS have more random oscillations due to velocity vectors. Another solution to this problem is the use of a large number of populations, but this method drastically increases the computational time, which results in the use of more powerful and expensive hardware. PSO is only used at Pbest and Gbest to update the position, but in PSC only those techniques that successfully track the GM have avoidance capability of LM trap and cohesion. DFO also shows LM avoiding trap by using a cohesion property with less tracking efficiency. GTOA provides a solution to these problems. Tracked power is presented in Figure 15, and a zoomed-in comparison is shown in Figure 16. The duty cycle comparison of GTOA is presented in Figure 17, and a zoomed in comparison of duty cycle transients is presented in Figure 18. Voltage and current transient comparisons are shown in Figures 19 and 20, respectively. This case dealt with the PS condition, as mentioned in Figure 6, which contains multiple LMs and one GM. PSC was used to check the robustness of the MPPT technique. The irradiance pattern was 650 W/m 2 , 520 W/m 2 , 250 W/m 2 , and 400 W/m 2 . For case 2, MPP was at 427 W. Irradiance pattern is shown in Table 2.
Randomly selected positions of swarm particles in initialization are a better solution to track GM under PS conditions, but take more time to converge at GM. PSO and CS have more random oscillations due to velocity vectors. Another solution to this problem is the use of a large number of populations, but this method drastically increases the computational time, which results in the use of more powerful and expensive hardware. PSO is only used at Pbest and Gbest to update the position, but in PSC only those techniques that successfully track the GM have avoidance capability of LM trap and cohesion. DFO also shows LM avoiding trap by using a cohesion property with less tracking efficiency. GTOA provides a solution to these problems. Tracked power is presented in Figure 15, and a zoomed-in comparison is shown in Figure 16. The duty cycle comparison of GTOA is presented in Figure 17, and a zoomed in comparison of duty cycle transients is presented in Figure 18. Voltage and current transient comparisons are shown in Figures 19 and 20

Case 3: Partial Shading
Case 3 dealt with a PS condition whose irradiance pattern was 720 W/m 2 , 850 W/m 2 , 1000 W/m 2 , and 550 W/m 2 . A power comparison is shown in Figure 21, and a zoomed-in comparison is presented in Figure 22. The duty cycle comparison of all the techniques in GTOA is shown in Figure 23, and the zoomed-in comparison of duty cycles in Figure

Case 3: Partial Shading
Case 3 dealt with a PS condition whose irradiance pattern was 720 W/m 2 , 850 W/m 2 , 1000 W/m 2 , and 550 W/m 2 . A power comparison is shown in Figure 21, and a zoomed-in comparison is presented in Figure 22. The duty cycle comparison of all the techniques in GTOA is shown in Figure 23, and the zoomed-in comparison of duty cycles in Figure Table 3. The CPS P-V curve is explained in Figure 6, which has two clusters, and GM occurs at 1136 W. CSA and PSO trap the local maxima due to iterative decrease in velocity vectors, which reduces the exploration phase and settles at LM. PSOGS also gets trapped in LM due to the restricted movement of particles caused by the gravitational effect, which can be observed in the power comparison and the zoomed-in comparison figures, presented in Figures 27 and 28, respectively. A comparison of duty cycle transients, and a zoomed-in duty cycle comparison, are presented in Figures 29 and 30, respectively, and which show the exploration and exploitation behavior of GTOA in CPS conditions.                   Table 4 shows the comparative analysis of the MPPT techniques, which shows that GTOA mitigates the drawbacks of conventional and SI based MPPT techniques; which are oscillations around GM, trapping of LM, and high tracking time. The oscillations of GTOA around GM was less than 0.5 W, which is 97.9% less than P&O. Tracking and settling time also exhibited the performance  Table 4 shows the comparative analysis of the MPPT techniques, which shows that GTOA mitigates the drawbacks of conventional and SI based MPPT techniques; which are oscillations around GM, trapping of LM, and high tracking time. The oscillations of GTOA around GM was less than 0.5 W, which is 97.9% less than P&O. Tracking and settling time also exhibited the performance of the MPPT technique. GTOA tracked the GM in 140 ms, and settled within 230 ms. On average GTOA was 13% to 35% faster in tracking the GM. Table 4. Quantitative analysis of GTOA with dragon fly optimization (DFO), particle swarm optimization (PSO), perturb and observe (P&O), cuckoo search (CS), and particle swarm optimization gravitational search (PSOGS).

Statistical Analysis
In this section statistical analysis of GTOA is presented, in Figure 33. The sensitivity of techniques was inspected by mean absolute error (MAE) in Equation (15), relative error (RE) in Equation (14), and root means square error (RMSE) in Equation (16). STC values of the PV modules were tested against the values given in the datasheet to test the standard. Standard errors indicated that the performance of the proposed GTOA technique in the transient, as well as steady state, had a higher success rate, which indicated the effective detection of GM under all conditions.
Electronics 2020, 9, 1962 20 of 24 where P pvi represents the power of the panel at STC, P pv represents tracked power, and n represents the number of samples. Figure 33 presents RMSE, MAE, and RE. Average power is given in Appendix C.

Conclusions
A novel GTOA based MPPT control technique is presented in this paper. This technique was experimentally evaluated against well-known P&O, DFO, PSO, PSOGS, and CS technology. The experimental simulation was carried out in five cases. Rapidly changing irradiance and PS conditions were studied. Cases 2 and 3 for PS condition, and special case 4 for complex partial shading were included. The experimental results, analysis results, and statistical analysis results showed that the proposed method has strong robustness, and is superior to competitive technology in tracking time, settling time, power tracking efficiency, and oscillation reduction. Compared with DFO, PSOGS, PSO, CS, and P&O technology, the proposed GTOA technology could track GMPP within 140 ms, while the settling time of the GTOA against competing techniques was almost 14% faster. The power efficiency of the GTOA was greater than 99.9%. Compared with other SI techniques, the average increase was 5%, however as compared to P&O it increased by 46% under PS. Overshoot, oscillation reduction, and robustness are the promising characteristics of the proposed GTOA technique. The

Conclusions
A novel GTOA based MPPT control technique is presented in this paper. This technique was experimentally evaluated against well-known P&O, DFO, PSO, PSOGS, and CS technology. The experimental simulation was carried out in five cases. Rapidly changing irradiance and PS conditions were studied. Cases 2 and 3 for PS condition, and special case 4 for complex partial shading were included. The experimental results, analysis results, and statistical analysis results showed that the proposed method has strong robustness, and is superior to competitive technology in tracking time, settling time, power tracking efficiency, and oscillation reduction. Compared with DFO, PSOGS, PSO, CS, and P&O technology, the proposed GTOA technology could track GMPP within 140 ms, while the settling time of the GTOA against competing techniques was almost 14% faster. The power efficiency of the GTOA was greater than 99.9%. Compared with other SI techniques, the average increase was 5%, however as compared to P&O it increased by 46% under PS. Overshoot, oscillation reduction, and robustness are the promising characteristics of the proposed GTOA technique. The GTOA showed significant efficiency and redundancy under extreme MPPT conditions. Steady-state oscillation was reduced to less than 0.5 W, and the transition from exploration to exploitation phase was very effective. The proposed GTOA successfully overcomes the deficiencies of the existing SI and classical MPPT controllers. In light of these comprehensive results and discussion, it is concluded that the proposed GTOA based control outperforms the existing techniques by a significant margin.  Initialize the maximum number of iterations and initialize the population using Equation (A5) within bounds given Equation (A6). X i = {i = 1, 2, 3, . . . , N} (A5) where U is the upper boundary, L i is the lower boundary, r is a random number.
Step 2: Fitness evaluation Compute fitness of each particle in population and select the best solution G.
Step 3: Stopping criteria: Check if number of iteration is greater than the maximum iterations, stop algorithm and output the best solution G is governed by Equation (A7).
where P PV is the power of panel. Above condition will re-initialize the code and execute step 1.
Step 4: Teacher Allocation Phase: Select top three best particles and calculate the Teacher T using Equation (13).
Step 5: Ability Group Phase: Divide population in two groups depending upon fitness. Best half of the population in outstanding group X best and others in average group X avg .
Step 6: Teacher and Student phases: For the best group X best , teacher phase is calculated using Equations (6)- (8) and (10). After that student phase is conducted using Equations (11) and (12).
For the average group X avg , teacher phase is implemented using Equations (9) and (10). After that, create a new population in X best . After that student phase is conducted using Equations (11) and (12). Create a new population of X avg .
Step 7: Reconstruction of population X best and X avg constructs a new population.
Step 8: Re-evaluation of population Evaluate the population and select the optimal solution G and increase the current iteration by Equation (A8). iter = iter + 1 (A8) Go to step 3 and execute. Figure A1 shows the summary of tracked power by competing MPPT techniques as compared to GTOA. Highest power is achieved by the proposed MPPT controller in all case studies by effectively minimizing the tracking time and power loss caused by partial shading.

Appendix C
Electronics 2020, 9, x FOR PEER REVIEW 24 of 26 Step 6: Teacher and Student phases: For the best group , teacher phase is calculated using Equations (6)- (8) and (10). After that student phase is conducted using Equations (11) and (12).
For the average group , teacher phase is implemented using Equations (9) and (10). After that, create a new population in . After that student phase is conducted using Equations (11) and (12). Create a new population of .
Step 7: Reconstruction of population and constructs a new population. Step