AMaximum Power Point Tracker with Automatic Step Size Tuning Scheme for Photovoltaic Systems

The purpose of this paper is to study on a novel maximum power point tracking (MPPT) method for photovoltaic (PV) systems. First, the simulation environment for PV systems is constructed by using PSIM software package. A 516 W PV system established with Kyocera KC40T photovoltaic modules is used as an example to finish the simulation of the proposed MPPT method. When using incremental conductance (INC) MPPT method, it usually should consider the tradeoff between the dynamic response and the steady-state oscillation, whereas the proposed modified incremental conductance method based on extension theory can automatically adjust the step size to track the maximum power point (MPP) of PV array and effectively improve the dynamic response and steady-state performance of the PV systems, simultaneously. Some simulation and experimental results are made to verify that the proposed extension maximum power point tracking method can provide a good dynamic response and steady-state performance for a photovoltaic power generation system.


Introduction
The characteristic curves representing the output currentvoltage (I-V ) and output power-voltage (P-V ) relations of PV modules will change in a nonlinear fashion, due to differing irradiation and ambient temperature conditions.Thus, maximum power tracking control must be imposed on the PV modules to make sure that power output stays at maximum any time of the day.For this reason, many experts and scholars have put forward all manner of proposals for a variety of MPPT mechanisms [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15][16][17].
One of the existing MPPT methods is the Perturb and Observe (P & O) method [1,2], which has a simple structure and requires very little parameter measurements.The direction of the next perturbation is determined by periodically increasing or decreasing the size of the output voltage of the PV module, and comparing it with the previous output voltage and power prior to any change.This is the most popular MPPT method used to track maximum power.However, even after P & O tracking approaches MPP, its perturbation does not stop there but oscillates around the MPP, resulting in the loss of energy.Though this phenomenon may be compensated for by setting a smaller perturbation with a smaller oscillation magnitude, if significant irradiation change takes place, it will cause the transient response time of the newly tracked MPP to increase.Yet, even if a larger perturbation can accelerate the response speed of tracking, it will result in a large increase in the amount of steady-state oscillation.Therefore, perturbation quantity needs to be appropriately selected as a performance trade-off between steady-state and transient responses.The maximum power tracking logic formula for the Incremental Conductance Method (INC) [3,4] is dI/dV = −I/V , where V is the output voltage, I is the output current, dI/dV is the dynamic conductance, and I/V is the static conductance.Traditional INC uses a fixed value for step adjustment of perturbation quantity during the tracking process.Similar to the P & O method, when the perturbation quantity is set to a smaller amount, the oscillation magnitude near the vicinity of the MPP will be smaller, making it easier to track the actual MPP.However, the tracking time will become relatively long.If perturbation quantity increases, International Journal of Photoenergy the tracking speed will accelerate.But that would increase the steady-state oscillation, resulting in a loss of energy and a reduction in the PV module's power output.Therefore, perturbation quantity depends on the equilibrium value between tracking accuracy and speed.Due to the complexity of the tracking process of INC method, the cost of the system will increase.Additionally, quantity measurement errors of the sensing element can further prevent the MPP from being tracked.However, to mitigate the defects of INC method, some experts and scholars have proposed a modified version of INC method [5] that uses a variable step size method of tracking and adds a fixed voltage tracking method to the tracking process.Thus, this method can be based on the characteristic curve of the PV module array to automatically adjust the size of each step.But the process involved in this method is highly complex, thus increasing the cost of the system.Factors such as switching the tracking system from a fixed-voltage tracking method to a modified INC method, or wrongfully setting the starting point of tracking, can result in system instability.Though the building of an intelligent fuzzy method for MPPT [6][7][8][9] can provide better performance for maximum power tracking of PV modules, these algorithms must be based on the output characteristics of PV modules to create rules for maximum power tracking control.For this reason, adaptive fuzzy logic [10] and parameter optimization techniques such as genetic algorithm [11] and particle swam optimization [12,13] have been applied to overcome the problem in MPPT algorithm.Due to the extremely complex nature of the computation process, the practicality of these control rules is quite limited.
A number of studies on MPPT have concentrated on the application of artificial neural network (ANN) [14,15].In most of these ANN-based methods, large number of field data considering atmospheric conditions are required to train the ANN.Moreover, the main problem of ANNbased methods is that it is system dependent and cannot be implemented for PV arrays with different characteristics.The combined use of fuzzy logic and ANN to track maximum power point in PV systems can be found in [16,17].In this method, ANN is trained offline using experimental data to desfine a reference voltage, which is the voltage at the maximum point according to the PV array characteristic.The reference voltage is then compared to the instantaneous array to generate a signal error.The signal error and change of the error are considered as the fuzzy logic controller (FLC) inputs.The FLC generates a duty cycle value for the pulse width modulation (PWM) generator.The PWM trigger signal is then applied to the switching of the boost converter connected to a PV array.A drawback of this method is that it needs much data for offline training.
In this paper, a novel intelligent technique based on extension theory is proposed and used together with INC based MPPT controller in a PV system.The proposed MPPT controller can adaptively tune the tracking step size of tradition INC MPPT method to obtain superior dynamic response and steady-state performance, simultaneously.The less constructed data utilized, no learning procedures needed, and easy implementation are the good features of the proposed MPPT method.

Extension Theory
In 1983 Professor Cai, a scholar from China, first proposed the concept of extension theory, which was designed to study things in terms of their extensity.This theory explores qualitative and quantitative solutions for contradictions among things [18].In addition, matter-element theory and extension mathematics are two core areas of extension theory.Matter-element theory mainly describes the extensity and transformation characteristics of matter-elements, while extension mathematics mainly focuses on the core calculus of the extension set and the correlation function [18].
Extension theory expresses information about things through the matter-element model, using the matterelement transformation to represent the changing correlation between the quality and quantity of the characteristics of matter.Then, the results from this correlation function are used to better understand the effects that these qualities and quantities have on the matter, so as to clearly express the level of impact from the characteristics of the matter.

The Concept of Extension Matter
Element.Extension theory handles problems through the use of the matter-element model.If such a model is represented by a mathematical function, it can be expressed asfollows: where R represents the description of the basic elements of matter, or matter element.N, C, and V represent the three elements that constitute the matter element.These three elements are N for the name of the matter, C for the characteristic value of matter, and V for the characteristic magnitude of matter.
In extension matter-element theory, when a characteristic of matter element is one of many, it is expressed as the x th characteristic and the x th corresponding characteristic value, respectively.
The vector form of this characteristic value is expressed as C = [c 1 , c 2 , . . ., c x ], and the characteristic magnitude is expressed as V = [v 1 , v 2 , . . ., v x ].Therefore, the extension matter element function of (1) can be rewritten as If the characteristic magnitude is an interval, then it is called a classical domain and is contained within a neighborhood domain.Assuming that an interval is set to between F 0 = a, b and F = d, e , with F 0 ∈ F, and that point f is any point on the interval F, then the corresponding matter element of F 0 = a, b can be expressed as C i is the characteristic value of F 0 , and V i is the characteristic magnitude of C i in the classical domain.The corresponding matter element of F, called R F , can be expressed as (4); where C j is the characteristic value of F, and V j is the characteristic magnitude of C j in the neighborhood domain:

Distance and Rank
Value.Classical mathematics measures the distance relationship between two points.Extension theory measures the distance relationship between any point on the actual domain and the interval, which is expressed as a function in (5): In addition to considering the correlation between points and intervals, there is also the need to consider the correlation between points and two intervals, or between intervals.Therefore, if we let F 0 = v a , v b and F = v d , v e be two separate intervals in the actual domain, with interval F 0 residing within interval F, then the rank value between point f and interval F 0 and interval F can be expressed as follows:

Correlation Function.
The correlation function is a function formed by dividing the rank value by distance, such as in the following: where when f = (v a + v b )/2, the value generated by the correlation function is at a maximum.Such a correlation function is called an elementary correlation function, as the diagram in Figure 1 shows.In addition, when λ g (f ) < −1, it means that point f is outside the interval F; when λ g (f ) > 0, it means that point f is inside the interval F 0 ; when −1 < λ g ( f ) < 0, it means that point f falls within the extension domain.

Extension Theory-Based MPPT Method
This study proposed an extension theory-based MPPT method designed to analyze the I-V characteristic curve of a photovoltaic module array using variation in slope errors as the characteristic of the proposed extension MPPT method and established 12 categories in the construction of the matter-element model.Calculations of degree of correlation related to the 12 categories were used to adjust the duty cycle of the MPPT boost converter circuit by category, to achieve maximum power tracking.

Characteristic Selection for the Proposed MPPT Method.
Figure 2 shows the schematic diagram of the I-V characteristic curve of a PV module array at various time points corresponding to its output current and voltage.Tradition incremental conductance MPPT method [3,4] is based on differentiation of PV power to its voltage and on condition of zero slope of P-V curve in maximum power point (MPP).Differencing PV power, especially, arises (8): From ( 8) and taking into account the basic condition of zero slope of P-V curve (9) is deduced: -I pv /V pv represents the opposite of static conductance of PV system and dI pv /dV pv the dynamic conductance.According to (9) these two quantities must be equal in MPP.Moreover, in the right of MPP is dI pv /dV pv < −I pv /V pv ; thus a reduction in PV's voltage is essential to achieve MPP.Similarly, in the left of MPP is dI pv /dV pv > −I pv /V pv ; thus an increase in PV's voltage is essential to achieve MPP.In this paper, to let the MPPT method possess adaptive capability, the step size of the INC MPPT method of the PV arrays is adaptively tuned by extension error tuning scheme.The tuning scheme is driven by a sum of static conductance I pv (k)/V pv (k) and dynamic conductance dI pv (k)/dV pv (k) which are defined as e(k) at kth sampling interval, respectively.The major purpose of this MPPT controller is to let the resulting dI pv /dV pv + I pv /V pv closely follow the reference dI pv /dV pv + I pv /V pv = 0 as shown in Figure 3 for driving the operation point toward the maximum power point.Thus, the general model e and ė trajectory can be predicted and plotted in Figure 4.In each category, the slope error e, change in the slope error ė, and change in the duty cycle ΔD were as shown in Table 1.In addition, the proximity between each category and the MPP was converted into a wave form to express the relationship between the 12 categories and the MPP; the results are displayed in Figure 5.As learned from Figures 4 and 5, Category numbers 3, 6, 9, and 12 showed larger amplitude swings because they were furthest away from the MPP.Category numbers 2, 5, 8, and 11 were closer to the MPP; category numbers 1, 4, 7, and 10 were closest to the MPP, so their swings in amplitude were smaller.

Establishment of a Matter-Element Model for the Characteristics of the Proposed MPPT Method.
Based on matter-element extension theory, the results from the matter-element model for a 12-category classical domain are shown in Table 2.The neighborhood domain constructed from the maximum and minimum values of each individual characteristic of the classical domain is displayed as follows: The following is the tracking control procedure using the proposed extension MPPT method.
Step 1. From slope error e and the amount of change in slope error ė in each category, construct a matter-element model: Step 2. Entering the to-be-classified slope error e and the amount of change in slope error ė, the matter-element model is Step 3. Based on the slope error e and the amount of change in slope error Δe, use (7) to calculate the correlation function between them and each individual category.
Step 4. Select weights W 1 and W 2 for each characteristic to represent the significance level of each characteristic.This study used a rule of thumb to select W 1 = 0.85 and W 2 = 0.15.Step 5. Calculate the degree of correlation with each category: W j K g j , g = 1, 2, . . ., 12.
Step 6. Select the maximum value from the calculated degree of correlation to identify the classified slope error e and the amount of change in slope error Δe into their respective categories.Setting the duty cycle variation ΔD according to the category to which it belongs, recalculate the new value of the duty cycle: where D old is the calculated duty cycle for the previous cycle, P is the polarity of error change ė, and λg is the maximum correlation degree of the category.

Simulation and Experimental Results
Figure 6 shows the system architecture of the proposed extension MPPT.Its maximum power tracking circuit architecture is a boost converter circuit; the design value of its circuit components is shown in Table 3 [19].With respect to simulation, the PSIM simulation software [20] in this study used Kyocera KC40T [21] as a reference module to construct a 516 W PV power generation system that applies variation in irradiation to simulate the proposed extension MPPT method and the traditional MPPT method and carried out a performance comparison.In terms of measurements, the PSoC chip made by Cypress Company was used to achieve the proposed extension MPPT strategy and performed the actual testing of the physical MPPT circuit under variable irradiation conditions.

Simulation Results
. This study carried out simulations under two different linear variations in irradiation conditions: one with an irradiation that dropped from 1,000 W/m 2 to 700 W/m 2 , then increased back to 1,000 W/m 2 ; and another with a irradiation that increased from 450 W/m 2 to 750 W/m 2 and then dropped back to 450 W/m 2 .These linear variations in irradiation are shown in Figures 7 and 8.Then, the proposed extension MPPT, traditional INC, and variable  step size INC methods were simulated and compared in terms of performance.The step variation (ΔD) of the traditional INC method was fixed at a value of 0.02, while the step variation of the variable step size INC method was set to 0.02 during the MPPT process.But when the power tracking was close to the MPP, step size variation changed from 0.02 to 0.01 to reduce the amount of steady-state oscillation.The change in the amount of step size change was based on the calculated e(k) Δ = dI pv (k)/dV pv (k) + I pv (k)/V pv (k) value; whether or not it was lower than the set value of 0.005 S formed the basis for changing the step variable quantity.The step changeability of the extension INC method was adjusted automatically based on the position of the working point as shown in Table 1.Figures 9,10,and 11 show the output waveforms of a PV module array using the traditional INC, variable step size INC, and proposed extension MPPT methods.As can be observed from Figures 9 to 11, the steady-state oscillation amplitude of the proposed extension MPPT method was the smallest of the three tracking methods.In addition, when irradiation changed from the linear interval to the constant value interval or changed from constant interval to linear interval, the proposed extension MPPT method has better  dynamic response than the traditional INC and variable step size INC methods.It can be seen that in addition to the good steady-state response of the proposed extension MPPT method, its transient response also has a good speed.Though the variable step size INC method can improve the amount of change in the steady-state response, its transient response slowed down.
In addition, to display excellence in maximum power tracking performance of the proposed extension MPPT method under various conditions of irradiation, this study followed a linear change in irradiation that increased from 450 W/m 2 to 750 W/m 2 and then dropped back down to 450 W/m 2 (irradiation changes as shown in Figure 8) to simulate the proposed tracking method.As can be seen from Figures 12, 13, and 14, regardless of whether it is under steady-state or transient response, the proposed extension MPPT method possesses similarly impressive performance, proving that the proposed tracking method has high robustness.

Measurement Results.
The PV module array used during the measurement was a solar simulator [22] made by    This shows that the proposed extension MPPT method has better performance in transient response in terms of tracking speed than the tradition INC and variable step size INC methods.In addition, the amount of steady-state oscillation is at a minimum after the proposed extension MPPT method has tracked the MPP.Thus, it also has a better steady-state performance than the traditional INC and the variable step size INC methods.That is, the power loss of the proposed extension MPPT method at the MPP is the least of the three MPPT methods.The output waveforms of the three tracking methods still oscillate upon achieving the MPP, due to maximum power tracking control still tracking back and forth in the vicinity of MPP.
In addition, to verify that the proposed extension MPPT method conforms with the simulated results, such that it still has good tracking performance for maximum power under various conditions of irradiation, the developed maximum power tracking circuit was used to carry out experimental testing under levels of irradiation that rose linearly from 450 W/m 2 to 750 W/m 2 and then dropped back down to 450 W/m 2 .As can be seen in Figures 18, 19 Figures 12 to 14.During the maximum power tracking process, it had better performance, regardless of transient or steady-state response.Table 4 shows the steady-state output power under various irradiation conditions using the extension MPPT method, the tradition INC method, and the variable step size INC method.To more clearly show that the steady-state performance of the proposed extension MPPT method is superior to that of the traditional INC and variable step size INC methods, this study postulated a single day irradiation of 1,000 W/m 2 , 800 W/m 2 , and 450 W/m 2 at sunshine times of 2 hrs, 3 hrs, and 4 hrs, respectively, for a comparison of the output power values from the three MPPT methods; the results of which are shown in Table 5.Table 5 shows that the output power value of the extension MPPT method was the highest of the three MPPT methods, which means that the power loss of the proposed extension MPPT method at the MPP was the least of the three MPPT methods.

Conclusion
This paper proposed an intelligent extension theory-based maximum power point tracking method for a 516 W PV

Figure 1 :
Figure 1: Schematic diagram of elementary correlation function.

Figure 2 :
Figure 2: Schematic diagram of P-V characteristic curve of output voltage and corresponding output power of a PV module array.

Figure 3 :Figure 4 :
Figure 3: Change diagram of I-V characteristic curve of slope error e and the change in slope error ė.

Figure 5 :
Figure 5: Dynamic slope analysis diagram of I-V characteristic curve.

Table 1 :−
Slope error e, change in slope error Δe, and change in duty cycle ΔD of the 12 categories.Category number Slope error category e Slope error change category ė Duty cycle step size ΔD 1 0 < e ≤ 0.00488 −125 < ė ≤ 0

Figure 6 :
Figure 6: Structural diagram of the proposed extension MPPT system.

Figure 7 :
Figure 7: Simulation waveform of irradiation that dropped from 1,000 W/m 2 to 700 W /m 2 and then increased back up to 1,000 W/m 2 .

Figure 8 :Figure 9 :
Figure 8: Simulation waveform of irradiation that increased from 450 W/m 2 to 750 W/m 2 and then dropped back to 450 W/m 2 .

Figure 10 :Figure 11 :
Figure 10: Simulated output power, dynamic conductance, static conductance, and duty cycle waveforms of a PV module array using the variable step size INC method under irradiation that dropped from 1,000 W/m 2 to 700 W/m 2 and then increased back up to 1,000 W/m 2 .

Figure 12 :
Figure 12: Simulated output power, dynamic conductance, static conductance, and sduty cycle waveforms of a PV module array using the traditional INC method under irradiation that increased from 450 W/m 2 to 750 W/m 2 and then dropped back to 450 W/m 2 .

Figure 13 :
Figure 13: Simulated output power, dynamic conductance, static conductance, and duty cycle waveforms of a PV module array using the variable step size INC method under irradiation that increased from 450 W/m 2 to 750 W/m 2 and then dropped back to 450 W/m 2 .

Figure 14 :Figure 15 :
Figure 14: Simulated output power, dynamic conductance, static conductance, and duty cycle waveforms of a PV module array using the proposed extension method under irradiation that increased from 450 W/m 2 to 750 W/m 2 and then dropped back to 450 W/m 2 .

Figure 16 :Figure 17 :
Figure 16: Measured output power, dynamic conductance, static conductance, and duty cycle waveforms of a PV module array using the variable step size INC method under irradiation that dropped from 1,000 W/m 2 to 700 W/m 2 and then increased back up to 1,000 W/m 2 .

Figure 18 :Figure 19 :
Figure 18: Measured output power, dynamic conductance, static conductance, and duty cycle waveforms of a PV module array using the traditional INC method under irradiation that increased from 450 W/m 2 to 750 W/m 2 and then dropped back to 450 W/m 2 .

Figure 20 :
Figure 20: Measured output power, dynamic conductance, static conductance, and duty cycle waveforms of a PV module array using the proposed extension method under irradiation that increased from 450 W/m 2 to 750 W/m 2 and then dropped back to 450 W/m 2 .

Table 2 :
The classical domains of categories.

Table 3 :
Component parameters of the implemented MPPT circuit., then rises again to 1,000 W/m 2 , as simulated by the solar simulator.As can be seen in Figures 15 to 17, at the 15 s mark, the irradiation began to change.The proposed extension MPPT method can achieve tracking of the new maximum power point (MPP) after about 1.1 s.As can be seen in Figures15 and 16, the tradition INC method and the variable step size INC method required 1.4 s and 1.6 s, respectively, to achieve tracking of the new MPP.

Table 4 :
Average value of steady-state power output under three types of irradiation using various MPPT methods.This study carried out simulation and measurement of maximum power tracking performance under variable irradiation conditions to verify the effectiveness of the proposed method.The simulated and measured results show that the proposed extension MPPT method possesses quick dynamic response to rapidly changing irradiation conditions.Its steady-state response at the MPP also has better performance than the traditional INC and variable step size INC MPPT methods.Meanwhile, the power loss of the proposed extension MPPT method at the MPP is far less

Table 5 :
Output energy during different hours of operation using various MPPT methods under three types of single day irradiation.