Research Article  Open Access
DataDriven Photovoltaic System Modeling Based on Nonlinear System Identification
Abstract
Solar photovoltaic (PV) energy sources are rapidly gaining potential growth and popularity compared to conventional fossil fuel sources. As the merging of PV systems with existing power sources increases, reliable and accurate PV system identification is essential, to address the highly nonlinear change in PV system dynamic and operational characteristics. This paper deals with the identification of a PV system characteristic with a switchmode power converter. Measured inputoutput data are collected from a real PV panel to be used for the identification. The data are divided into estimation and validation sets. The identification methodology is discussed. A HammersteinWiener model is identified and selected due to its suitability to best capture the PV system dynamics, and results and discussion are provided to demonstrate the accuracy of the selected model structure.
1. Introduction
The modern power system is increasingly taking advantage of renewable energy sources entering the marketplace. Traditional central power stations with their pollution related problems will likely be replaced with cleaner and smaller power plants closer to the loads. The energy generated by the sun is one of the most promising, nonpolluting, free sources of energy [1]. Among their benefits, solarpowered systems are easily expanded. Despite their still relatively high cost, photovoltaic (PV) systems installed worldwide show a nearly exponential increase [2]. A PV cell directly converts sunlight into electricity, and the basic elementary device of PV systems is the PV cell [3]. PV systems have proven that they can generate power to very small electronic devices up to utilityscale PV power plants. The basic building block for PV systems is a PV panel consisting of a number of prewired cells in series [4]. Panels are then connected in series to increase voltage and in parallel to increase current; the product is power. A PV array is formed by series and parallel combinations of panels [5]. Figure 1 represents a generic PV array structure. The performance of a PV system is normally evaluated under the standard test condition (STC), where an average solar spectrum at 1.5 Air Mass (AM) is used, the irradiance is normalized to 1000 W/m^{2}, and the cell temperature is defined as 25°C [6]. However, under real operating conditions (i.e., varying irradiance as well as significant temperature changes), most commercial panels do not necessarily behave as in the specifications given by the manufacturers [7, 8]. In addition, PV panels perform differently according to the location, time of day, and season of the year.
For PV system, the relationship between environmental conditions and electrical output parameters (current and voltage) is highly nonlinear. For this reason, the process of modeling the dynamics of the PV system and identifying the model structure that captures reallife behavior is extremely essential for the purpose of controlling the output power. Another important purpose is to predict future performance of the system for maintenance and troubleshooting.
Modeling and simulation of PV systems have been the subject of many research studies [3, 9–12]. However, the focus of the studies was on the PV panel/array stage without including the power electronics of the complete system—that is, identifying only the nonlinear IV characteristic of the PV system. In this work, the PV system identification incorporated the entire setup considering a PV system, power converter, maximum power point tracker (MPPT), and a load.
Identification of linear/nonlinear systems became a hot topic in the 1960s, probably because proper models were needed to design good controllers [13]. At this stage, the system identification field was rather immature. The first attempt to put the field into order was the published pioneering work of Eykhoff [14]. Afterward, the field became completely structured when the books of Ljung [15] and Soderstrom and Stoica [16] laid out a complete theoretical framework and practical methodology to the identification process.
The core of the system identification process is to construct a mathematical model from observed inputoutput data. It is widely applied in different engineering disciplines to help understand the studied process, predict the system responses, and create better design with new specifications.
This work is organized as follows. Section 2 presents a description of the studied PV system. Section 3 details the system identification methodology. Section 4 discusses the results associated with implementing the PV system identification from applied inputoutput data, dividing the data into estimation and validation sets. Finally, Section 5 concludes the work.
2. PV System Description
The PV system illustrated in Figure 2 is considered for this study. The PV panel performance is affected by several environmental conditions; however, irradiance has the strongest impact on the panel output power. The panel is a multipurpose module consisting of 36 polycrystalline silicon cells connected in series. Under STC, the open circuit voltage () is 21.6 V, the short circuit current () is 5.16 A, the maximum power voltage () is 17.3 V, and the maximum power current () is 4.63 A. Figure 3 presents the IV and PV characteristics under different irradiances. These data are obtained experimentally [17], and they are not the data that will be used for the system identification procedure; nevertheless, prior information about the system is beneficial. What we are interested in is the current measurement as input and the power measurement as output with respect to time, as indicated in Figure 2.
The IV and PV characteristic curves show the standard PV panel/array behavior; as the irradiance increases, the panel current is directly proportional to the solar intensity per unit area. Cutting the irradiance in half leads to a drop in current by half. Decreasing irradiance also reduces voltage, but it does so following a logarithmic relationship that results in a relatively modest change of voltage.
In general, the operating point of the PV panel is not the maximum power point (MPP) along the IV curve, at which the panel operates with maximum efficiency and produces maximum power. Thus, a switchmode power converter (DC/DC buck converter, in our case) with a MPPT is utilized to maintain the PV panel operating point at the MPP. The MPPT controller retains this operating point by controlling the PV panel’s voltage or current independently of those of the load.
Depending on the application, the DC/DC converter stage could be a buck converter, boost converter, or buckboost converter. Several MPPT control algorithms have been proposed in different publications. Perturbation and observation (P&O) or the hill climbing algorithm is the most commonly employed method in commercial PV MPPTs [18–22]. Among other control mechanisms for MPPT are fuzzy logic [23–27], neural network [28–32], and incremental conductance (IncCond) control [33–37]. The IncCond algorithm is adapted in this study because it is more convenient and lends itself well to DSP and microcontroller control.
3. Identification Methodology
In this section, the principle of system identification is described. Dynamic models depend heavily on the amount of a priori knowledge about the dynamic process that is to be incorporated. Generally, modeling any dynamical system can be categorized into two approaches: firstprinciple modeling and datadriven modeling.
Firstprinciple modeling uses an understanding of the system’s physics to derive a mathematical representation, whereas datadriven modeling involves using empirical data to construct a model for the system. The two approaches classify the modeling process in terms of known parameters and structure into black box, grey box, and white box modeling as suggested by Figure 4.
According to [38], white box modeling is when a model is perfectly known; it is possible to construct the model entirely from prior knowledge and physical insight. Grey box modeling is when some physical insight is available, but several parameters remain to be determined from observed data. In black box modeling, no physical insight is available or used, but the chosen model structure belongs to candidate models that are known to have good flexibility and have been successful in the past.
The black box modeling is used for this work. It is usually a trialanderror process where one estimates the parameters of various structures and compares the results. The PV system constructed and implemented in Figure 2 is used for data collection purposes. The measured inputoutput data will then be utilized to select model structure and identify the system.
3.1. The PV System Identification Process
Figure 5 illustrates the measured data collected to be used for the system identification. The measured inputoutput data were taken from the PV system setup shown in Figure 2. When the irradiance is changing with time, the output maximum power will vary according to the maximum power values in Figure 3. The STC (1000 W/m^{2}) power is 82.5 W for the selected panel. The irradiance fluctuates, thereby causing the input current outputted from the panel to change. This results in continuous change of the output power of the system.
The flowchart in Figure 6 summarizes the identification process using the measured data. The inputoutput data are a sampled time domain signal. The next step is to select a suitable structure for the model. Selecting the model structure is the most difficult step in the identification process because there are a multitude of possibilities. Next, one can use the Matlab System Identification Toolbox and compare models to choose the model with the best performance. The final step is to evaluate the resulting model with a validation data set. It is worth mentioning that system identification is not a straightforward technique; rather, it is an iterative procedure involving several decisions to be made during the process.
3.2. HammersteinWiener Model
A nonlinear blockoriented model is frequently applied for an adequate description of the nonlinear behavior of a system over a range of operating conditions. The identified system is generally subdivided into a linear dynamic block and a nonlinear static block. The HammersteinWiener model is a nonlinear model that is used in many domains for its simplicity. The model is popular because it has convenient block representation, has a transparent relationship to linear systems, and is easier to implement than heavyduty nonlinear models such as neural networks and Volterra models. It can be used as a black box model structure because it offers flexible parameterization for nonlinear models. The model describes a dynamical system using one or two static nonlinear blocks in series with a linear block. Only the linear block contains dynamic elements. The linear block is a discretetime transfer function, and the nonlinear blocks are implemented using nonlinearity estimators such as saturation, wavelet, and dead zone. Figure 7 depicts the structure of the nonlinear HammersteinWiener model. The input signal passes through the first nonlinear block, a linear block, and a second nonlinear block to produce the output signal [39].
The HammersteinWiener structure can be described by the following general equations:where and are the inputs and outputs for the system, respectively. and are nonlinear functions that correspond to the input and output nonlinearities, respectively. For multiple inputs and multiple outputs, and are defined independently for each input and output channel. and are internal variables that define the input and output of the linear block, respectively. has the same dimension as . has the same dimension as . and in the linear dynamic block are linear functions. For outputs and inputs, the linear block is a transfer function matrix containing entries in the following form:where and . If only the input nonlinearity is present, the model is called a Hammerstein model. If only the output nonlinearity is present, the model is called a Wiener model.
The available nonlinearity estimators to be used in the identification process, by estimating the parameters of the input and output blocks, are dead zone, piecewise linear, saturation, sigmoid network, and wavelet network [40].
3.2.1. Dead Zone Function
The dead zone function generates zero output within a specified region that is called the dead zone. The lower and upper limits of the dead zone are specified as the start and the end of the dead zone parameters, respectively. The dead zone defines a nonlinear function , where is a function of . There are three intervals, which can be identified as follows: when has a value between and ; this zone is called the “zero interval” zone. See Figure 8.
3.2.2. Piecewise Linear Function
The piecewise linear function is defined as a nonlinear function , where is a piecewise linear (affine) function of and there are breakpoints where and . is linearly interpolated between the breakpoints. and are scalars.
3.2.3. Saturation Function
The saturation function can be defined as a nonlinear function , where is a function of . There are three intervals as shown in Figure 9 and they can be identified as follows:
3.2.4. Sigmoid Network Function
The sigmoid network nonlinear estimator uses neural networks comprising an input layer, an output layer, and a hidden layer employing sigmoid activation functions as represented by Figure 10.
It combines the radial basis neural network function using a sigmoid as the activation function. The estimator is based on the following expression: where is the input and is the output. is the regressor. is a nonlinear subspace and is a linear subspace. is a linear coefficient. is an output offset. is a dilation coefficient, is a translation coefficient, and is an output coefficient. is the sigmoid function, given by the following equation:
3.2.5. Wavelet Network Function
The wavelet estimator is a nonlinear function combining wavelet theory and neural networks. Wavelet networks are feedforward neural networks using a wavelet as an activation function, based on the following expression: where is the input and is the output, is a nonlinear subspace and is a linear subspace. is a linear coefficient. is an output offset. is a scaling coefficient and is a wavelet coefficient. is a scaling dilation coefficient and is a wavelet dilation coefficient. is a scaling translation coefficient and is a wavelet translation coefficient. The scaling function and the wavelet function are both radial functions and can be written as follows:In the system identification process, the wavelet coefficient , the dilation coefficient , and the translation coefficient are optimized during model learning steps to obtain the best performance model.
4. Results and Discussion
As mentioned earlier in Section 3, the measured inputoutput data collected in Figure 5 are the system identification data to be used for the estimation. Table 1 shows the values of irradiance, current, and power with time for the PV system.

Once the identification process results in finding the best model in a structure, another data set (validation data) is used to validate the model. The identification procedure follows the flowchart steps in Figure 6.
4.1. Model Estimation
Table 1 represents the variation in power for the data set collected for the estimation. The time domain data in Figure 5 are imported to the system identification toolbox in Matlab. Such a toolbox integrates techniques for nonlinear and linear models so that the complex problem of estimating and analyzing nonlinear models appears simple and systematic.
In our case, a model is to be identified for the PV system from the current data as input to the system, and the recorded power is to be identified as the system output. A nonlinear autoregressive exogenous (NARX) model was selected as a reference point. There are several models that can be selected from the toolbox. Figure 11(a) demonstrates the available classifications of system identification; however, discussing the detail of each one is beyond the scope of this paper.
(a)
(b)
(c)
(d)
Figure 11(b) shows the result of a nonlinear ARX model type with wavelet network nonlinearity. The figure indicates the model output matches the measured output with 77.48% accuracy. However, this match does not best describe the dynamics of the PV system. The input and output data each contain 100001 samples with sampling interval of 0.0001 seconds.
After the iterative process—by selecting different nonlinear models with their associated inputoutput nonlinearity, searching methods, and manipulating numerous estimation configurations—HammersteinWiener was found to be the best model that captures the main system dynamics, with a 94.51% fit. See Figure 11(c).
4.2. Model CrossValidation
Similarly, the model performance is evaluated using crossvalidation. In crossvalidation, the HammersteinWiener model is confronted with a new data set that is different from the data used to estimate the model. Figure 11(d) provides the new measured data set collected from the PV system, this time with different irradiances, namely, 800 W/m^{2}, 600 W/m^{2}, 400 W/m^{2}, and 800 W/m^{2}. This situation emulates an effect of a cloud blocking the irradiance from the PV system. In the same figure, it can be seen that the output of the HammersteinWiener model matches the validation data set well, with 93.98% accuracy. This shows that the estimation process is robust enough to handle different input and that the model estimation is successful. Table 2 lists the results of the best fit with different input and output nonlinearity showing the 93.98% accuracy. The table as well displays an assessment of the estimation in terms of final prediction error (FPE) and loss function.
 
FE: failed estimation. 
5. Conclusions
The results of this work showed that the nonlinear HammersteinWiener model was able to provide an accurate description of the PV system dynamics. The model was selected after an iterative process involving trial and error. The model must produce an accurate fit to the data, and the HammersteinWiener model showed an accuracy of 93.98% applying estimation and validation data sets. Developing such a black box model from measured inputoutput data can contribute to the design and implementation of nonlinear control strategies. In addition, researchers and engineers will be able to predict future performance of PV systems for maintenance and troubleshooting.
Competing Interests
The authors declare that they have no competing interests.
References
 W. T. da Costa, J. F. Fardin, L. D. V. B. M. Neto, and D. S. L. Simonetti, “Identification of photovoltaic model parameters by differential evolution,” in Proceedings of the IEEE International Conference on Industrial Technology (ICIT '10), pp. 931–936, IEEE, March 2010. View at: Publisher Site  Google Scholar
 D. Sera, R. Teodorescu, and P. Rodriguez, “PV panel model based on datasheet values,” in Proceedings of the IEEE International Symposium on Industrial Electronics (ISIE '07), pp. 2392–2396, Vigo, Spain, June 2007. View at: Publisher Site  Google Scholar
 M. G. Villalva, J. R. Gazoli, and E. R. Filho, “Comprehensive approach to modeling and simulation of photovoltaic arrays,” IEEE Transactions on Power Electronics, vol. 24, no. 5, pp. 1198–1208, 2009. View at: Publisher Site  Google Scholar
 G. M. Masters, Renewable and Efficient Electric Power Systems, John Wiley & Sons, Hoboken, NJ, USA, 2004.
 R. Ramaprabha and B. L. Mathur, “MATLAB based modeling to study the influence of shading on series connected SPVA,” in Proceedings of the 2nd International Conference on Emerging Trends in Engineering and Technology (ICETET '09), pp. 30–34, IEEE, Nagpur, India, December 2009. View at: Publisher Site  Google Scholar
 W. Xiao, W. G. Dunford, and A. Capel, “A novel modeling method for photovoltaic cells,” in Proceedings of the IEEE 35th Annual Power Electronics Specialists Conference (PESC '04), vol. 3, pp. 1950–1956, Aachen, Germany, June 2004. View at: Publisher Site  Google Scholar
 E. Saloux, M. Sorinand, and A. Teyssedou, “Explicit model of photovoltaic panels to determine voltages and currents at the maximum power point,” CETC Number 2010156, CanmetENERGY, Ottawa, Canada, 2010. View at: Google Scholar
 A. H. Alqahtani, M. S. Abuhamdeh, and Y. M. Alsmadi, “A simplified and comprehensive approach to characterize photovoltaic system performance,” in Proceedings of the IEEE Energytech Conference, pp. 1–6, Cleveland, Ohio, USA, May 2012. View at: Publisher Site  Google Scholar
 J. J. Soon and K.S. Low, “Photovoltaic model identification using particle swarm optimization with inverse barrier constraint,” IEEE Transactions on Power Electronics, vol. 27, no. 9, pp. 3975–3983, 2012. View at: Publisher Site  Google Scholar
 Y. A. Mahmoud, W. Xiao, and H. H. Zeineldin, “A parameterization approach for enhancing PV model accuracy,” IEEE Transactions on Industrial Electronics, vol. 60, no. 12, pp. 5708–5716, 2013. View at: Publisher Site  Google Scholar
 S. A. Rahman, R. K. Varma, and T. Vanderheide, “Generalised model of a photovoltaic panel,” IET Renewable Power Generation, vol. 8, no. 3, pp. 217–229, 2014. View at: Publisher Site  Google Scholar
 L. Cristaldi, M. Faifer, M. Rossi, and S. Toscani, “An improved modelbased maximum power point tracker for photovoltaic panels,” IEEE Transactions on Instrumentation and Measurement, vol. 63, no. 1, pp. 63–71, 2014. View at: Publisher Site  Google Scholar
 J. Schoukens, “System identification: what does it offer to electrical engineers,” in Proceedings of the 60th Anniversary of the Foundation of the Faculty of Electrical Engineering and Informatics of the Technical, University Budapest, December 2009. View at: Google Scholar
 P. Eykhoff, System IdentificationParameter and State Estimation, Wiley, London, UK, 1974.
 L. Ljung, System Identification—Theory for the User, PrenticeHall, Englewood Cliffs, NJ, USA, 1987.
 T. Soderstrom and P. Stoica, System Identification, PrenticeHall, Hemel Hempstead, UK, 1989.
 C. Osorio, ModelBased Design for Solar Power Systems, The MathWorks, Natick, Mass, USA, 2009.
 K. Ishaque and Z. Salam, “A review of maximum power point tracking techniques of PV system for uniform insolation and partial shading condition,” Renewable and Sustainable Energy Reviews, vol. 19, pp. 475–488, 2013. View at: Publisher Site  Google Scholar
 A. R. Reisi, M. H. Moradi, and S. Jamasb, “Classification and comparison of maximum power point tracking techniques for photovoltaic system: a review,” Renewable and Sustainable Energy Reviews, vol. 19, pp. 433–443, 2013. View at: Publisher Site  Google Scholar
 M. A. G. de Brito, L. Galotto, L. P. Sampaio, G. De Azevedo Melo, and C. A. Canesin, “Evaluation of the main MPPT techniques for photovoltaic applications,” IEEE Transactions on Industrial Electronics, vol. 60, no. 3, pp. 1156–1167, 2013. View at: Publisher Site  Google Scholar
 E. B. Youssef, P. Stephane, E. Bruno, and A. Corinne, “New P&O MPPT algorithm for FPGA implementation,” in Proceedings of the 36th Annual Conference of the IEEE Industrial Electronics Society (IECON '10), pp. 2868–2873, IEEE, Glendale, Calif, USA, November 2010. View at: Publisher Site  Google Scholar
 E. Setiawan, I. Hodaka, Y. Yamamoto, and I. Jikuya, “Improvement of maximum power point tracking based on nonlinear control of boost converter,” in Proceedings of the 1st IEEE Conference on Control, Systems & Industrial Informatics (ICCSII '12), pp. 180–184, Bandung, Indonesia, September 2012. View at: Publisher Site  Google Scholar
 C.S. Chiu and Y.L. Ouyang, “Robust maximum power tracking control of uncertain photovoltaic systems: a unified TS fuzzy modelbased approach,” IEEE Transactions on Control Systems Technology, vol. 19, no. 6, pp. 1516–1526, 2011. View at: Publisher Site  Google Scholar
 N. Patcharaprakiti and S. Premrudeepreechacharn, “Maximum power point tracking using adaptive fuzzy logic control for gridconnected photovoltaic system,” in Proceedings of the IEEE Power Engineering Society Winter Meeting, pp. 372–377, IEEE, January 2002. View at: Google Scholar
 S. A. Khan and M. I. Hossain, “Design and implementation of microcontroller based fuzzy logic control for maximum power point tracking of a photovoltaic system,” in Proceedings of the International Conference on Electrical and Computer Engineering (ICECE '10), pp. 322–325, IEEE, Dhaka, Bangladesh, December 2010. View at: Publisher Site  Google Scholar
 J. Li and H. Wang, “Maximum power point tracking of photovoltaic generation based on the fuzzy control method,” in Proceedings of the 1st International Conference on Sustainable Power Generation and Supply (SUPERGEN '09), pp. 1–6, Nanjing, China, April 2009. View at: Publisher Site  Google Scholar
 S.J. Kang, J.S. Ko, J.S. Choi et al., “A novel MPPT Control of photovoltaic system using FLC algorithm,” in Proceedings of the 11th International Conference on Control, Automation and Systems (ICCAS '11), pp. 434–439, IEEE, Gyeonggido, Republic of Korea, October 2011. View at: Google Scholar
 W.M. Lin, C.M. Hong, and C.H. Chen, “Neuralnetworkbased MPPT control of a standalone hybrid power generation system,” IEEE Transactions on Power Electronics, vol. 26, no. 12, pp. 3571–3581, 2011. View at: Publisher Site  Google Scholar
 M. Syafaruddin, E. Karatepe, and T. Hiyama, “Artificial neural networkpolar coordinated fuzzy controller based maximum power point tracking control under partially shaded conditions,” IET Renewable Power Generation, vol. 3, no. 2, pp. 239–253, 2009. View at: Publisher Site  Google Scholar
 T. Esram, Modeling and control of an alternatingcurrent photovoltaic module [Ph.D. thesis], Department of Electrical Engineering, University of Illinois at UrbanaChampaign, Urbana, Ill, USA, 2010.
 Y.H. Liu, C.L. Liu, J.W. Huang, and J.H. Chen, “Neuralnetworkbased maximum power point tracking methods for photovoltaic systems operating under fast changing environments,” Solar Energy, vol. 89, pp. 42–53, 2013. View at: Publisher Site  Google Scholar
 A. B. G. Bahgat, N. H. Helwa, G. E. Ahmad, and E. T. El Shenawy, “Maximum power point traking controller for PV systems using neural networks,” Renewable Energy, vol. 30, no. 8, pp. 1257–1268, 2005. View at: Publisher Site  Google Scholar
 B.J. Jeon, K.H. Kim, K. Lee, J.H. Im, G.B. Cho, and Y.O. Choi, “Battery controller design of standalone photovoltaic system using IncCond method,” in Proceedings of the 15th International Conference on Electrical Machines and Systems (ICEMS '12), pp. 1–5, October 2012. View at: Google Scholar
 C. Osorio, Webinar: ModelBased Design for Solar Power Systems, The MathWorks, Natick, Mass, USA, 2009.
 H. Kumar and R. K. Tripathi, “Simulation of variable incremental conductance method with direct control method using boost converter,” in Proceedings of the Students Conference on Engineering and Systems (SCES '12), March 2012. View at: Publisher Site  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
 A. Safari and S. Mekhilef, “Implementation of incremental conductance method with direct control,” in Proceedings of the IEEE Region 10 Conference (TENCON '11), pp. 944–948, Bali, Indonesia, November 2011. View at: Publisher Site  Google Scholar
 J. Sjöberg, Q. Zhang, L. Ljung et al., “Nonlinear blackbox modeling in system identification: a unified overview,” Automatica, vol. 31, no. 12, pp. 1691–1724, 1995. View at: Publisher Site  Google Scholar  MathSciNet
 L. Ljung, System Identification Toolbox 7 User's Guide, The MathWorks, Natick, Mass, USA, 2008.
 N. Patcharaprakiti, K. Kirtikara, K. Tunlasakun et al., “Modeling of photovoltaic grid connected inverters based on nonlinear system identification for power quality analysis,” in Electrical Generation and Distribution Systems and Power Quality Disturbances, InTech, Rijeka, Croatia, 2011. View at: Google Scholar
Copyright
Copyright © 2016 Ayedh Alqahtani et al. 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.