Research Article  Open Access
BICO MPPT: A Faster Maximum Power Point Tracker and Its Application for Photovoltaic Panels
Abstract
This paper develops a maximum power point tracking (MPPT) algorithm to optimize photovoltaic (PV) array performance and make it more compatible to rapidly varying weather conditions. In particular, a novel extremum seeking controller (ESC), which uses a Bode ideal cutoff (BICO) filter in its structure, is designed and tested on a simulated PV array. The new algorithm is compared against the commonly used ESC MPPT algorithm with firstorder filters. The BICO extremum seeking controller achieves transient rise to the MPP faster than the common extremum seeking MPPT, which is the faster and more robust method among all other methods. This claim has been discussed and proved mathematically in this paper, in addition to simulation illustrations. This faster extremum seeking algorithm enables PV systems to detect rapid variations in the environmental conditions like irradiation and temperature changes.
1. Introduction
The PV cells exhibit a nonlinear  characteristics as shown in Figure 1 and their power output mainly depends on the nature of the connected load. Since connecting the load directly to the PV system results in a poor overall efficiency, to minimize the life cycle cost of any PV system, increasing the efficiency by changing the operating point of the system using an intermediate maximum power point tracker (MPPT) is desirable.
MPPT controls the output current and voltage and consequently output power of the PV panel adaptively to maintain maximum efficiency and better performance in the presence of environmental variations. Typically MPPT algorithms are implemented on a solar array using a switching power converter; for instance, in a gridtied inverter, the solar array charges a capacitor and then the current is switched out of the capacitor at an optimal varying duty cycle in order to extract maximum power from the PV array.
A number of solar power converter architectures with MPPT are discussed in the literature [1, 2]. As discussed in these works, convergence speed is one of the most important features among all different MPPT algorithms. Brunton et al. have pointed out in their paper that, “as irradiance decreases rapidly, the curve shrinks and the MPV and MPI decrease. If the MPPT algorithm does not track fast enough, the control current or voltage will fall off the curve” [3]. Consequently, any improvement in the rise time of MPPT improves the reliability of the system and increases the power extraction and efficiency of the whole system.
1.1. Maximum Power Point Tracking Algorithms
There are many different maximum power point tracking techniques for photovoltaic systems which are well established in the literature [4–11]. These techniques vary in many aspects as simplicity, convergence speed, digital or analogical implementation, sensors required, cost, range of effectiveness, and other aspects. In analog world, short current (SC), open voltage (OV), and CV are good options for MPPT; otherwise with digital circuits that require the use of microcontroller, Perturb and Observe (P&O), Incremental Conductance (IC), and temperature methods are easy to implement [12]. Table 1 and Figure 2 present the comparison among different MPPT methods considering the costs of sensors, microcontroller, and additional power components and also their efficiencies. In this table, A means absence, L means low, M means medium, and H means high.

Currently, the most popular and workhorse MPPT algorithm is P&O, because of its balance between performance and simplicity. However, it suffers from the lack of speed and adaptability which is necessary for tracking the fast transients under varying environmental conditions [3]. Recently, another adaptive algorithm, called extremum seeking control, has been developed [13] to overcome these weaknesses.
1.2. Extremum Seeking Control
A promising new robust MPPT algorithm is the method of extremum seeking (ES) control, which carries all P&O’s benefits like simplicity and performance, and in addition, improves its weaknesses [3]. Figures 3 and 4 present the comparison between ESC and P&O, where the inverter controls the current and voltage according to the setpoints which are provided by these two MPPT algorithms.
(a)
(b)
(c)
(a)
(b)
(c)
The power, current, and voltage are plotted versus time for the ESC and P&O algorithms as well as the true maximum power of the system.
P&O and ESC methods oscillate closely around the real maximum power voltage, as seen in the power versus time plot. Obviously, the ESC method rises to the MPP orders of magnitude more rapidly than the P&O.
ESC MPPT has some advantages from hardware implementation point of view. Brunton et al. have mentioned in their paper that “the ripplebased ES algorithm has good MPPT performance over a range of inverter capacitor sizes. Typically, the choice of capacitor is expensive because it must be well characterized and large enough to maintain a small ripple. However, because the ES control signal exploits the natural inverter ripple, a smaller capacitor allows the tracking of rapid irradiance changes.”
ESC for peak power point tracking method has been successfully applied to biochemical reactors [14, 15], ABS control in automotive brakes [16], variable cam timing engine operation [17], electromechanical valves [18], axial compressors [19], mobile robots [20], mobile sensor networks [21, 22], optical fibre amplifiers [23], and so on. A good survey of the literature on this topic prior to 1980 can be found in [24] and a more recent overview can be found in [25]. Astrom and Wittenmark rated extremum seeking as one of the most promising adaptive control methods [26].
Since extremum seeking control has better features and performance compare to P&O which is the best known MPPT algorithm, in this paper, the improvement of better than the best MPPT algorithm has been investigated.
2. Basic Regular Extremum Seeking Algorithm
As shown in Figure 5, ESC algorithm employs a slow periodic perturbation, , which is added to the estimated signal, . If the perturbation is slow enough, the plant appears as a static map and its dynamics do not interfere with the peak seeking scheme [13]. If is on either side of , which is the optimal point, the perturbation signal , will create a periodic response of which is either in phase or out of phase with . The highpass filter eliminates the “DC component" of . Thus and highpass filter will be approximately two sinusoids which are in phase if or out of phase if .
The integrator approximates the gradient update law which tunes to [13]. System of Figure 5 can be summarized in mathematical equations as follows: where is the convolution operator and is the inverse Laplace transform operator. The transfer functions for H_{HPF} and H_{LPF} in the regular SISO ESC scheme are and respectively, where [25]. This model will be used for stability analysis in this paper.
In the following sections, after introducing BICO filter [27, 28], the advantages of using this filter in the ESC algorithm from the stability and robustness point of view will be discussed.
3. BICO Extremum Seeking Control MPPT
In this section we are going to introduce a filter which was strongly favored by Bode [29], called Bode’s ideal cutoff characteristic (BICO) filter [27, 28]. The general transfer function of lowpass BICO filter is,
Its corresponding time response in the special case (rad/s) was derived by Oberhettinger and Baddi [30], but the time response of (2) is presented as, where is the th order Bessel function.
By replacing by , highpass BICO filter is obtained. Figure 6 compares the frequency response of highpass BICO filter and regular highpass filter () with the same cutoff frequency of rad/s.
As seen in this figure, BICO filter has a sharp edge in its cutoff frequency. This great feature causes almost no attenuation for frequencies higher than and a large attenuation in the lower frequencies. Therefore the behavior of this filter is close to an “ideal” filter. This sharp edge presents in the lowpass BICO filter. By combining highpass and lowpass BICO filters, the bandpass BICO filter with sharp edges in both sides can be obtained as well.
4. Stability Analysis of Extremum Seeking Control
4.1. Mathematical Modeling of ESC Scheme
The stability analysis of ESC algorithm has been investigated in [13, 31–34]. In these literatures traditional extremum seeking control with regular firstorder filters has been considered. In this paper, we have considered the similar stability analysis approach as [34] to compare the stability and robustness of BICOESC and regular SISO ESC. The difference between this work and [34] is that in this paper, the stability analysis has been done by considering both lowpass and highpass filters in the ESC structure but in [34] authors have investigated the stability of simplified ESC with only lowpass filter in its structure.
According to the previous discussion, the nonlinear map in the ESC scheme is considered to be concave and is assumed to have only one extremum point. Since in the PV system applications, MPPT is employed to extract maximum amount of power from PV panels, therefore, extremum point in this case is maximum point (), where at this point , and also . In ESC algorithm, the output of the nonlinear map is, where and are the amplitude and angular frequency of perturbation signal. Since the perturbation signal assumes to be small, therefore the Taylor expansion of (4) is where H.O.T. stands for higher order terms and is the approximation of . The highpass filtered signal will be,
Highpass filter acts as a derivative opereator in series with a lowpass filter (). By multiplying the modulation signal to the resulted signal from the highpass filter and passing the modulated signal through the lowpass filter and since , the output signal of integrator will be
Under the assumptions that the amplitude of the sinusoidal perturbation is small and the harmonics of highpass filter are attenuated by lowpass filter, output of lowpass filter is proportional to the gradient of the nonlinear map with respect to its input and time. Therefore in the neighborhood of the extremum point, the amplitude of the estimated signal is small, since the gradient is small. It can be seen that this amplitude depends on and .
4.2. Stability Analysis of Averaged ESC Scheme
The averaging method is typically used to analyze the periodic steady state solutions of weakly nonlinear systems. Since the amplitude of perturbation in the ESC scheme is small, this system can be evaluated by its averaged model. The averaged form of signal is, where . Therefore the averaged model of ESC scheme is
On the other hand, in the neighborhood of the extremum point, can be approximated as,
If the difference between the extremum point and averaged point is defined by and since , thus,
By defining , then, . Therefore, . Substituting this relationship in (9) gives,
Without loss of generality and by assuming that is the output signal from a highpass filter, H_{HPF}, which its cutoff frequency is higher than cutoff frequency of the lowpass filter, H_{LPF}, then (12) can be rewritten as
This system can be considered as a feedback system as shown in Figure 7, which can be considered as a perturbation. The loop gain of this system depends on the demodulation gain, , the integral gain, , and the curvature of nonlinear map, . This result completely matches with the results obtained from other analysis method in the simplified case [13, 31–33].
5. Analysis of ESC MPPT on PV Panels
5.1. PV Cell Model
The PV cell  curve is modeled using the light emitting diode equations [3]:
Values and definitions for further simulations are shown in Table 2.

5.2. BICO ESC vsersus Regular SISO ESC
In order to compare the qualitative behavior of BICO MPPT with the regular SISO ESC MPPT, the RootLocus (RL) analysis has been chosen. Clearly, other stability analysis methods are applicable at this point.
Since RootLocus (RL) method is serving for linear systems, it is important to point out that the curvature constant, , is an uncertain parameter which depends on different factors like the PV panel manufacturer and weather conditions. According to [33], .
As seen in Figure 7, the characteristic polynomial of averaged ESC system is , and to analyze the behavior of this system, RootLocus (RL) analysis method can be employed. Since there is no command in MATLAB to plot RootLocus (RL) for BICO transfer function, this equation has been solved for different and the roots have been plotted. To compare the behavior of regular SISO ESC with BICO ESC, the roots for this system has been plotted with the same method instead of using “rlocus” command.
Figure 8 shows the comparison between the root locus of averaged ESC using BICO and regular firstorder filters. In this root locus plot, the constant gain is assumed to be and cutoff frequencies in both filters are assumed to be rad/s. Clearly, by using firstorder filters in the averaged ESC scheme, for some values of , system has complex poles which cause an oscillatory behavior in the response of the system. On the other side, by using BICO filters, roots of the characteristic polynomial (system poles) are always real numbers for any value of .
Besides the location of the poles, as can be seen in this figure, for the same value of , BICO filter has farther poles with respect to the origin compared to the firstorder filter case. Therefore, the bandwidth of BICO filter is higher than the firstorder filter with similar cutoff frequency. Higher bandwidth means faster response for BICO ESC. This feature is justifiable by looking at the frequency response of both filters. Since the bandwidth of BICO is close to the bandwidth of an ideal filter, BICO ESC becomes the fastest achievable MPPT algorithm.
6. Simulation Illustrations
Figures 9 and 10 show the maximum power point of a PV panel with the defined parameters in Table 2. As can be seen in these figures, environmental conditions and especially variations in the sun irradiation will change the nonlinear behavior of PV panels. Shadows, cloudy or dusty weather, and temperature variations cause moving of optimal operating point in PV systems. When temperature increases, the maximum output power of PV panels decreases and vice versa. In Figure 9 the variations of optimal operating point by changing the environmental temperature from 20°C to 30°C are illustrated. As Figure 10 presents, variation of peak output power happens in the wider range when sun irradiation changes. Generally, in PV panels, when sun irradiation increases, the output power increases, but when temperature increases, the maximum extractable power reduces.
To compare the proposed MPPT method which is called BICO MPPT with the ESC MPPT, the working conditions of both algorithms have been defined to be same. For all simulations, the ambient temperature is 25°C and the irradiation is assumed to be . Also the cutoff frequency of highpass filters is rad/s, and for lowpass filter, this frequency is rad/s. Under these conditions, from Figure 9, the maximum amount of power which can be extracted from the simulated PV panel is 48 Watts and this maximum happens around 17 Volts. To implement the BICO MPPT, the discrete approximation of this filter in MathWorks, Inc.’s website has been used.
Figure 11 shows the outputs of extremum seeking algorithms resulted from two different MPPT methods. Figure 12 represents the maximum voltage tracking in BICO MPPT algorithm and ESC with firstorder filters. As expected and proved before, BICO MPPT converges to the peak power point two times faster than the regular ESC MPPT algorithm. Faster convergence speed is due to the higher bandwidth of the averaged BICO system.
Figures 13 and 14 show the performance of the proposed MPPT approach compare to the ESC in the presence of white noise. As can be seen in these figures, in the presences of a white noise (noise power ), which is considered as the variations in the nonlinear map behavior, , BICO MPPT preforms better than ESC MPPT from attenuation and tracking point of view.
Another parameter which is considered in these simulations is the robustness of these two methods to the variations of the system gain. Figures 15(a) and 15(b) illustrate the performance of BICO MPPT and regular ESC MPPT to the gain variations, respectively. From these figures, it can be concluded that, BICO MPPT tolerates higher gain variations, but by increasing the gain of the integrator in ESC MPPT, it starts oscillating.
(a)
(b)
7. Conclusions
The ubiquity of using MPPT in the renewable energy systems has stimulated the persistent development of various MPPT algorithms. According to the comparison in this paper, one of the most popular methods of maximum power tracking method is the P&O method. However, this method fails in those situations when there is a need to track the MPPs rapidly.
ESC is a new robust MPPT algorithm which carries all P&O’s benefits and improves its weaknesses. In this paper, a novel ES algorithm has been presented which is called BICO MPPT. This algorithm employs BICO filter in its structure to improve the ESC rise time response even further. BICO filter is a forgotten part of Bode’s research that for the first time in this paper was used for PV applications.
BICO filter is the closest filter to the ideal filter and this feature of BICO improves the performance of ESC algorithm. Also, the advantages of BICO has been discussed in this paper. As discussed in this paper, using BICO in the ESC structure increases the bandwidth of the system, which improves the system response.
In addition, applying BICO filter to the ESC algorithm moves all the roots of characteristic polynomial of averaged BICO ESC system to the real axis. Consequently, this system has no pole with imaginary part, and therefore theoretically the system has no oscillatory response by increasing its integrator gain. These poles, however, can have imaginary part by increasing the gain of the system in the regular ESC. This feature increases the robustness of the BICO MPPT compared to regular ESC.
As can be seen in the results, BICO MPPT not only can follow the maximum power point faster than regular ESC, but also it shows more robustness in the presence of disturbance in the system or gain variations.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
References
 B. Bose, P. Szczesny, and R. Steigerwald, “Microcomputer control of a residential photovoltaic power conditioning system,” IEEE Transactions on Industry Applications, vol. 21, no. 5, pp. 1182–11191, 1985. View at: Google Scholar
 R. A. Mastromauro, M. Liserre, T. Kerekes, and A. Dell'Aquila, “A singlephase voltagecontrolled gridconnected photovoltaic system with power quality conditioner functionality,” IEEE Transactions on Industrial Electronics, vol. 56, no. 11, pp. 4436–4444, 2009. View at: Publisher Site  Google Scholar
 S. L. Brunton, C. W. Rowley, S. R. Kulkarni, and C. Clarkson, “Maximum power point tracking for photovoltaic optimization using ripplebased extremum seeking control,” IEEE Transactions on Power Electronics, vol. 25, no. 10, pp. 2531–2540, 2010. View at: Publisher Site  Google Scholar
 T. Esram and P. L. Chapman, “Comparison of photovoltaic array maximum power point tracking techniques,” IEEE Transactions on Energy Conversion, vol. 22, no. 2, pp. 439–449, 2007. View at: Publisher Site  Google Scholar
 J. J. Nedumgatt, K. B. Jayakrishnan, S. Umashankar, D. Vijayakumar, and D. P. Kothari, “Perturb and observe MPPT algorithm for solar PV systemsmodeling and simulation,” in Proceedings of the Annual IEEE India Conference: Engineering Sustainable Solutions (INDICON '11), pp. 1–6, December 2011. View at: Publisher Site  Google Scholar
 R. Alonso, P. Ibáñez, V. Martínez, E. Román, and A. Sanz, “An innovative perturb, observe and check algorithm for partially shaded PV systems,” in Proceedings of the 13th European Conference on Power Electronics and Applications (EPE '09), pp. 1–8, September 2009. View at: Google Scholar
 N. Femia, G. Petrone, G. Spagnuolo, and M. Vitelli, “Optimization of perturb and observe maximum power point tracking method,” IEEE Transactions on Power Electronics, vol. 20, no. 4, pp. 963–973, 2005. View at: Publisher Site  Google Scholar
 G. J. Kish, J. J. Lee, and P. W. Lehn, “Modelling and control of photovoltaic panels utilising the incremental conductance method for maximum power point tracking,” IET Renewable Power Generation, vol. 6, no. 4, pp. 259–266. View at: Google Scholar
 A. Safari and S. Mekhilef, “Simulation and hardware implementation of incremental conductance MPPT with direct control method using cuk converter,” IEEE Transactions on Industrial Electronics, vol. 58, no. 4, pp. 1154–1161, 2011. View at: Publisher Site  Google Scholar
 H. Malek, S. Dadras, Y. Q. Chen, R. Burt, and J. Cook, “Maximum power point tracking techniques for efficient photovoltaic microsatellite power supply system,” in Proceeding of the Small Satellite Conference, Logan, Utah, USA, 2012. View at: Google Scholar
 H. Malek, S. Dadras, and Y. Q. Chen, “A fractional order maximum power point tracker: Stability analysis and experiments,” in Proceeding of the IEEE 51st Annual Conference on Decision and Control (CDC '12), pp. 6861–6866, Maui, Hawaii, USA, 2012. View at: Google Scholar
 R. Faranda and S. Leva, “Energy comparison of MPPT techniques for PV systems,” WSEAS Transactions on Power Systems, vol. 3, no. 6, pp. 1–6, 2008. View at: Google Scholar
 M. Krstić and H.H. Wang, “Stability of extremum seeking feedback for general nonlinear dynamic systems,” Automatica, vol. 36, no. 4, pp. 595–601, 2000. View at: Publisher Site  Google Scholar
 M. Guay, D. Dochain, and M. Perrier, “Adaptive extremum seeking control of continuous stirred tank bioreactors with unknown growth kinetics,” Automatica, vol. 40, no. 5, pp. 881–888, 2004. View at: Publisher Site  Google Scholar
 G. Bastin, D. Nešić, Y. Tan, and I. Mareels, “On extremum seeking in bioprocesses with multivalued cost functions,” Biotechnology Progress, vol. 25, no. 3, pp. 683–689, 2009. View at: Publisher Site  Google Scholar
 S. Drakunov, U. Ozguner, P. Dix, and B. Ashrafi, “ABS control using optimum search via sliding modes,” IEEE Transactions on Control Systems Technology, vol. 3, no. 1, pp. 79–85, 1995. View at: Publisher Site  Google Scholar
 D. Popović, M. Janković, S. Magner, and A. Teel, “Extremum seeking methods for optimization of variable cam timing engine operation,” in Proceeding of the American Control Conference, pp. 3136–3141, Denver, Colo, USA, June 2003. View at: Google Scholar
 K. S. Peterson and A. G. Stefanopoulou, “Extremum seeking control for soft landing of an electromechanical valve actuator,” Automatica, vol. 40, no. 6, pp. 1063–1069, 2004. View at: Publisher Site  Google Scholar
 H.H. Wang, S. Yeung, and M. Krstić, “Experimental application of extremum seeking on an axialflow compressor,” IEEE Transactions on Control Systems Technology, vol. 8, no. 2, pp. 300–309, 2000. View at: Publisher Site  Google Scholar
 C. G. Mayhew, R. G. Sanfelice, and A. R. Teel, “Robust sourceseeking hybrid controllers for autonomous vehicles,” in Proceedings of the American Control Conference (ACC '07), pp. 1185–1190, New York, NY, USA, July 2007. View at: Publisher Site  Google Scholar
 E. Biyik and M. Arcak, “Gradient climbing in formation via extremum seeking and passivitybased coordination rules,” in Proceedings of the 46th IEEE Conference on Decision and Control (CDC '07), pp. 3133–3138, New Orleans, La, USA, December 2007. View at: Publisher Site  Google Scholar
 P. Ögren, E. Fiorelli, and N. E. Leonard, “Cooperative control of mobile sensor networks: adaptive gradient climbing in a distributed environment,” IEEE Transactions on Automatic Control, vol. 49, no. 8, pp. 1292–1302, 2004. View at: Publisher Site  Google Scholar
 P. M. Dower, P. M. Farrell, and D. Nešić, “Extremum seeking control of cascaded Raman optical amplifiers,” IEEE Transactions on Control Systems Technology, vol. 16, no. 3, pp. 396–407, 2008. View at: Publisher Site  Google Scholar
 J. Sternby, “Extremum control systems: an area for adaptive control,” in Proceedings of the International Symposium on Adaptive Systems, vol. 24, pp. 151–160, San Francisco, Calif, USA, 1980. View at: Google Scholar
 K. B. Ariyur and M. Krstic, RealTime Optimization by ExtremumSeeking Control, WileyInterscience, 2003.
 K. J. Astrom and B. Wittenmark, Adaptive Control, AddisonWesley, 2nd edition, 1995.
 Y. Luo, T. Zhang, B. Lee, C. Kang, and Y. Q. Chen, “Disturbance observer design with bode's ideal cutoff filter in harddiscdrive servo system,” Mechatronics, vol. 23, no. 7, pp. 856–862, 2013. View at: Google Scholar
 L. E. Olivier, I. K. Craig, and Y. Q. Chen, “Fractional order and BICO disturbance observers for a runofmine ore milling circuit,” Journal of Process Control, vol. 22, no. 1, pp. 3–10, 2012. View at: Publisher Site  Google Scholar
 T. T. Hartley and C. F. Lorenzo, “Optimal fractionalorder damping,” in Proceedings of the ASME Design Engineering Technical Conferences and Computers and Information in Engineering Conference, pp. 737–743, Chicago, Ill, USA, September 2003. View at: Google Scholar
 Oberhettinger and Baddi, Tables of Laplace Transforms, Springer, Berlin, Germany, 1973.
 Y. Tan, D. Nešić, and I. Mareels, “On stability properties of a simple extremum seeking scheme,” in Proceedings of the 45th IEEE Conference on Decision and Control (CDC '06), pp. 2807–2812, December 2006. View at: Google Scholar
 Y. Tan, D. Nešić, and I. Mareels, “On nonlocal stability properties of extremum seeking control,” Automatica, vol. 42, no. 6, pp. 889–903, 2006. View at: Publisher Site  Google Scholar
 R. Leyva, C. Olalla, H. Zazo et al., “Mppt based on sinusoidal extremumseeking control in pv generation,” International Journal of Photoenergy, vol. 98, no. 4, pp. 529–542, 2011. View at: Google Scholar
 R. Leyva, P. Artillan, C. Cabal, B. Estibals, and C. Alonso, “Dynamic performance of maximum power point tracking circuits using sinusoidal extremum seeking control for photovoltaic generation,” International Journal of Electronics, vol. 98, no. 4, pp. 529–542, 2011. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Hadi Malek and YangQuan Chen. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.