Nguyen and Nguyen Environ Syst Res (2015) 4:24

DOI 10.1186/s40068-015-0047-9

RESEARCH Open Access

Mathematical modeling of photovoltaic

cell/module/arrays with tags in Matlab/Simulink
Xuan Hieu Nguyen1* and Minh Phuong Nguyen2

Abstract 
Background:  Photovoltaic (PV) array which is composed of modules is considered as the fundamental power con-
version unit of a PV generator system. The PV array has nonlinear characteristics and it is quite expensive and takes
much time to get the operating curves of PV array under varying operating conditions. In order to overcome these
obstacles, common and simple models of solar panel have been developed and integrated to many engineering
software including Matlab/Simulink. However, these models are not adequate for application involving hybrid energy
system since they need a flexible tuning of some parameters in the system and not easily understandable for readers
to use by themselves. Therefore, this paper presents a step-by-step procedure for the simulation of PV cells/modules/
arrays with Tag tools in Matlab/Simulink. A DS-100M solar panel is used as reference model. The operation characteris-
tics of PV array are also investigated at a wide range of operating conditions and physical parameters.
Result:  The output characteristics curves of the model match the characteristics of DS-100M solar panel. The output
power, current and voltage decreases when the solar irradiation reduces from 1000 to 100 W/m2. When the tempera-
ture decreases, the output power and voltage increases marginally whereas the output current almost keeps con-
stant. Shunt resistance has significant effect on the operating curves of solar PV array as low power output is recorded
if the value of shunt resistance varies from 1000 ohms to 0.1 ohms.
Conclusion:  The proposed procedure provides an accurate, reliable and easy-to-tune model of photovoltaic array.
Furthermore, it also robust advantageous in investigating the solar PV array operation from different physical parame-
ters (series, shunt resistance, ideality factor, etc.) and working condition ( varying temperature, irradiation and espe-
cially partial shadow effect) aspects.
Keywords:  Photovoltaic array, Matlab/Simulink, P–V and I–V curve, Simulation, Tag

Background and (Gonzalez-Longatt 2005). Here, the effect of tem-

Mathematical modeling of PV module is being con- perature, solar irradiation, and diode quality factor and
tinuously updated to enable researchers to have a better series resistance is evaluated. A difficulty of this method
understanding of its working. The models differ depend- is to require readers programming skills so it is not easy
ing on the types of software researchers used such as to follow. Another method which is the combination
C-programming, Excel, Matlab, Simulink or the tool- between Matlab m-file and C-language programming is
boxes they developed. even more difficult to clarify (Gow and Manning 1999).
A function in Matlab environment has been developed Among other authors, a proposed model is based on
to calculate the current output from data of voltage, solar solar cell and array’s mathematical equations and built
irradiation and temperature in the study of (Walker 2001) with common blocks in Simulink environment in (Salmi
et  al. 2012), (Panwar and Saini 2012), (Savita Nema and
Agnihotri 2010), and (Sudeepika and Khan 2014). In
*Correspondence: [email protected]; [email protected] these studies, the effect of environmental conditions
1
Department of Electric Power System, Faculty of Engineering, Vietnam (solar insolation and temperature), and physical param-
National University of Agriculture, Trau Quy town, Gia Lam district,
Hanoi 10000, Vietnam eters (diode’s quality factor, series resistance Rs, shunt
Full list of author information is available at the end of the article

© 2015 Nguyen and Nguyen. This article is distributed under the terms of the Creative Commons Attribution 4.0 International
License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any
medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons
license, and indicate if changes were made.
Nguyen and Nguyen Environ Syst Res (2015) 4:24 Page 2 of 13

resistance Rsh, and saturation current, etc.) is investi- Methods

gated. One disadvantage of these papers is lack of pre- Mathematical equivalent circuit for photovoltaic array
senting simulation procedure so it causes difficulties for The equivalent circuit of a PV cell is shown in Fig. 1. The
readers to follow and simulate by themselves later. This current source Iph represents the cell photocurrent. Rsh
disadvantage is filled in by (Jena et al. 2014), (Pandiarajan and Rs are the intrinsic shunt and series resistances of
and Muthu 2011). A step-by-step procedure for simulat- the cell, respectively. Usually the value of Rsh is very large
ing PV module with subsystem blocks with user-friendly and that of Rs is very small, hence they may be neglected
icons and dialog in the same approach with Tarak Salmi to simplify the analysis (Pandiarajan and Muthu
and Savita Nema is developed  by Jena, Pandiarajan and 2011).  Practically, PV cells are grouped in larger units
Muthu et al. However, the biggest gap of the studies men- called PV modules and these modules are connected in
tioned above is shortage of considering the effect of par- series or parallel to create PV arrays which are used to
tially shading condition on solar PV panel’s operation. generate electricity in PV generation systems. The equiv-
In other researches, authors used empirical data and alent circuit for PV array is shown in Fig. 2.
Lookup Table or Curve Fitting Tool (CFtool) to build The voltage–current characteristic equation of a
P–V and I–V characteristics of solar module (Banu and solar cell is provided as (Tu and Su 2008; Salmi et  al.
Istrate 2012). The disadvantage of this method is that it is 2012):  Module photo-current Iph:
quite challenging or even unable to collect sufficient data
Iph = [Isc + Ki (T − 298)] × Ir/1000 (1)
if no experimental system be available so  that modeling
curves cannot be built and modeled. Here, Iph: photo-current (A); Isc: short circuit current (A)  ;
From the work of (Ibbini et  al. 2014) and (Ven- Ki: short-circuit current of cell at 25 °C and 1000 W/m2;
kateswarlu and Raju 2013), a solar cell block which has T: operating temperature (K); Ir: solar irradiation (W/m2).
already been built in Simscape/Simulink environment is Module reverse saturation current Irs:
employed. With this block, the input parameters such as
short circuit current, open circuit voltage, etc. is provided
Irs = Isc /[exp(qVOC /NS knT ) − 1] (2)
by manufacturers. The negative point of this approach is Here, q: electron  charge,  =  1.6  ×  10 C; Voc: open −19

that some parameters including saturation current, tem- circuit voltage  (V); Ns: number of cells connected in
perature, and so on cannot be evaluated. series; n: the ideality factor of the diode; k: Boltzmann’s
Solar model developed with Tag tools in Simulink envi- constant, = 1.3805 × 10−23 J/K.
ronment is recorded in the research of (Varshney and Tariq
2014), (Mohammed 2011), etc. In these papers, only two
aspects (solar irradiation and temperature) are investigated
without providing step-by-step simulation procedure.
In overall, although having advantages and disadvan-
tages, different methods have similar gaps as follows:

••  The proposed models are not totally sufficient to

study all parameters which can significantly affect
to I–V and P–V characteristics of solar PV array,
including physical parameters such as saturation Fig. 1  PV cell equivalent circuit (Salmi et al. 2012)
current, ideality factor, series and shunt resistance,
etc. and environmental working conditions (solar
insolation, temperature and especially shading
effect).
••  Lack of presenting step-by-step simulation procedure
and this causes difficulties for readers and research-
ers to follow and do simulation by themselves.

Therefore, the study proposes a robust model built with

Tag tools in Simulink environment. The proposed model
shows strength in investigating all parameters’ influence
on solar PV array’s operation. In addition, a unique step-
by-step modeling procedure shown allows readers to fol-
low and simulate by themselves to do research. Fig. 2  Equivalent circuit of solar array (Tu and Su 2008)
Nguyen and Nguyen Environ Syst Res (2015) 4:24 Page 3 of 13

The module saturation current I0 varies with the cell Reference model
temperature, which is given by: The 100 W solar power module is taken as the reference
module for simulation and the detailed parameters of
module is given in Table 1.
 3   
T q × Eg0 1 1
I0 = Irs exp − (3)
Tr nk T Tr

Here, Tr: nominal temperature = 298.15 K; Eg0: band gap

energy of the semiconductor, = 1.1 eV; The current out-
put of PV module is:
Table 1 Electrical characteristics data of  DS-100  M PV
I = NP × Iph − NP × I0 module
   
V /NS + I × Rs /NP Name DS-100M
× exp − 1 − Ish (4)
n × Vt
Rated power (Vmp) 100 W
With
Voltage at maximum power (Vmp) 18 V
k ×T Current at maximum power (Imp) 5.55 A
Vt = (5) Open circuit voltage (VOC) 21.6 V
q
Short circuit current (ISC) 6.11 A
and
Total number of cells in series (NS) 36
V × NP /NS + I × RS Total number of cells in parallel (NP) 1
Ish = (6) Maximum system voltage 1000 V
Rsh
Range of operation temperature −40 °C to 80 °C
Here: Np: number of PV modules connected in paral-
lel; Rs: series resistance  (Ω); Rsh: shunt resistance  (Ω); The electrical specifications are under test conditions of irradiance of 1 kW/m2 ,
spectrum of 1.5 air masses and cell temperature of 25 °C
Vt: diode thermal voltage (V).

Fig. 3  Input parameters for simulation model

Nguyen and Nguyen Environ Syst Res (2015) 4:24 Page 4 of 13

Step by step procedure for modeling of photovoltaic arrays current at 25  °C and 1000  W/m2  =  6.11 A; Voc is PV
with tags module open circuit voltage at 25  °C and 1000  W/
A mathematical model of PV array including fundamen- m2  =  0.6  V; Eg0 is the band gap energy for sili-
tal components of diode, current source, series resistor con  =  1.1  eV. Rs is series resistor, normally the value
and parallel resistor is modeled with Tags in Simulink of this one is very small,  =  0.0001  Ω; Rsh is shunt
environment (http://mathwork.com). The simulation  of resistor, the value of this is so large, = 1000 Ω (Fig. 3).
solar module is based on equations given in the section
Step 2
above and done in the following steps.
Module photon-current is given in Eq.  (1)  and modeled
Step 1
as Fig. 4 (Ir0 = 1000 W/m2).
Provide input parameters for modeling:
Iph = [Isc + Ki (T − 298)] × Ir/1000 (7)
Tr is reference temperature  =  298.15  K; n
is ideality factor  =  1.2; k is Boltzmann con- Step 3
stant  =  1.3805  ×  10−23 J/K; q is electron Module reverse saturation current is given in Eq. (2) and
charge  =  1.6  ×  10−19; Isc is PV module short circuit modeled as Fig. 5.

Fig. 4  Modeled circuit for Eq. (1)

Fig. 5  Modeled circuit for Eq. (2) and (5)

Nguyen and Nguyen Environ Syst Res (2015) 4:24 Page 5 of 13

Fig. 6  Modeled circuit for Eq. (3)

Fig. 7  Modeled circuit for Eq. (6)

Nguyen and Nguyen Environ Syst Res (2015) 4:24 Page 6 of 13

Fig. 8  a Modeled circuit for Eq. (4). b Solar subsystem simulation model

Nguyen and Nguyen Environ Syst Res (2015) 4:24 Page 7 of 13

Fig. 9  Simulation model of solar PV array

Step 4 and a solar system of two DS-100M panels connected

in series, each panel has the key specifications listed in
Module saturation current I0 is given in Eq. (3) and mod-
Table 1.
eled as Fig. 6.
Step 5 Result and discussion
Simulation scenario
Modeled circuit for Eq. (6) (Fig. 7).
With the developed model, the PV array characteris-
Step 6 tics are estimated as follows.
Modeled circuit for Eq. (4) (Fig. 8).
(i) I–V and P–V characteristics under varying irradia-
Step 7 tion with constant temperature are given in Fig. 11a
and b. Here, the solar irradiation changes with val-
The solar module simulation procedure is shown from Fig. 3
ues of 100, 500 and 1000  W/m2 while temperature
to Fig. 8b. The solar PV array includes six modules and each
keeps constant at 25 °C.
module has six solar cells connected in series. Therefore, the
Summary when the irradiation increases, the cur-
proposed model of solar PV array is given in Fig. 9.
rent and voltage output increase. This results in rise
in power output in this operating condition.
Experimental test
(ii) I–V and P–V characteristics under varying tem-
In order to validate the Matlab/Simulink model, the PV
perature and constant irradiation are obtained in
test system of Fig.  10 is installed. It consists of a rheo-
Fig.  12a and b. Here, the temperature varies from
stat, a solar irradiation meter, two digital multi-meters
Nguyen and Nguyen Environ Syst Res (2015) 4:24 Page 8 of 13

slightly and this results in slight net reduction in

power output. However, a significant decrease in
current, voltage and power output is recorded when
the value of Rsh is 0.1 Ω.
(iv) 
I–V and P–V characteristics under varying Ns
and Np are obtained in Fig.  14a and b. In practice,
PV cells are connected in series into PV module
and these PV modules then are connected in series
or parallel to form PV array for generating more
electricity from sunlight. The reference model is
36-series-connected-cell array so two cases are
studied: two modules are connected in series and
two modules are connected in parallel.
Summary
Fig. 10  Setup of the experimental solar system •  With two modules connected in series (Ns  =  72,
Np = 1), the value of current output is similar to that
of it in case of one module (Ns = 36, Np = 1) but the
25 to 50 and 75 °C respectively whereas the irradia- voltage output doubles so the power output doubles.
tion level keeps constant at 1000 W/m2. •  In term of two modules connected in paral-
Summary when the operating temperature increases, lel (Ns = 36, Np = 2), the value of voltage output is
the current output raises marginally but the voltage similar to that of it in case of one module (Ns = 36,
output decreases drastically. This leads to net reduc- Np = 1) but the current output doubles so the power
tion in power output with rise in temperature. output doubles. Similar value of power output is
(iii) I–V and P–V characteristics under varying shunt/paral- experienced in both cases of two modules despite dif-
lel resistance Rsh, constant temperature and irradiation ferent ways in module connection (parallel or series).
are shown in Fig.  13a and b. In this case, Rsh changes
with three values of 0.1, 1 and 1000 Ω, respectively. The proposed model has advantages not only in studying
 Summary when Rsh varies between 1000 and 1 Ω, effect of physical parameters such as series resistance Rs,
the current output and voltage output decreases shunt resistance Rsh, etc. but also in investigating impact

Fig. 11  a Output I–V curve with varying irradiation. b Output P–V curve with varying irradiation
Nguyen and Nguyen Environ Syst Res (2015) 4:24 Page 9 of 13

Fig. 12  a Output I–V curve with varying temperature. b Output P–V curve with varying temperature

of environmental condition like varying temperature, Case Description

solar irradiation and especially shading effect. In this
3 Two shaded modules (receive irradi-
study, the evaluation of shading effect on solar PV array’s ation of 500 W/m2), others receive
operation is carried out through following cases. The full irradiation of 1000 W/m2
simulation results are given in Fig. 15a and b. 4 Two shaded modules (receive
irradiation of 500 and 250 W/
m2), others receive irradiation of
Case Description 1000 W/m2
1 No shaded PV module (full irradia-
tion on solar PV array): 1000 W/m2
2 One shaded module (receives irradi-
Summary
ation of 500 W/m2), others receive
full irradiation of 1000 W/m2
•  The power output of PV array reduces  noticeably
when it works under partial shading condition.
Nguyen and Nguyen Environ Syst Res (2015) 4:24 Page 10 of 13

Fig. 13  a Output I–V curve with varying Rsh. b Output P–V curve with varying Rsh

Fig. 14  a I–V curve with varying Ns and Np. b P–V curve with varying Ns and Np
Nguyen and Nguyen Environ Syst Res (2015) 4:24 Page 11 of 13

Fig. 15  a I–V curves under partial shading condition. b P–V curves under partial shading condition

•  The I–V curve experiences multiple steps whereas Conclusion

the P–V curve gives many local peaks along with A step-by-step procedure for simulating a PV array with
the maximum power point (the global peak). In Tag tools, with user-friendly icons and dialogs in Matlab/
addition, more shaded modules are higher number Simulink block libraries is shown. This modeling proce-
of power output peaks is shown. dure serves as an aid to help people to closer understand
of I–V and P–V operating curves of PV module. In addi-
Experimental results and validation tion, it can be considered as a robust tool to predict the
The Matlab/Simulink model is evaluated for the experi- behavior of any solar PV cells, modules and arrays under
mental test system (two DS-100M panels are connected varying environmental conditions (temperature, irradia-
in series). The results are shown in Fig. 16. On the other tion and partially shading condition) and physical param-
hand, the empirical results with a solar irradiation of eters (series resistance, shunt resistance, ideality factor
520 W/m2 and operating temperature of 40 °C are given and so on). This research is the first step to study a hybrid
in Fig. 17. The I–V and P–V simulation and experimental system where a PV power generation connecting to other
results show a good agreement in terms of short circuit renewable energy production sources like wind or bio-
current, open circuit voltage and maximum power output. mass energy systems.
Nguyen and Nguyen Environ Syst Res (2015) 4:24 Page 12 of 13

Fig. 16  Solar Matlab/SIMULINK results

Fig. 17  Solar system experimental results

Authors’ contributions Authors’ information

XHN initiated, proposed model developed in Matlab/Simulink and analyzed. Xuan Hieu Nguyen is a lecturer in Department of Electric Power System, Fac-
He also prepared a draft manuscript for publication. MPN assisted in designing, ulty of Engineering, Vietnam National University of Agriculture, Hanoi, Vietnam.
data collection, analysis and reviewed the manuscript and edited many times He had bachelor degree in electric power system, Hanoi University of Science
and added his inputs. Finally, XHN decided finally the content of the research and Technology, Vietnam in 2008. He received master degree in electrical
revised final manuscript. All authors read and approved the final manuscript. engineering in University of Wollongong, Australia in 2012. Minh Phuong
Nguyen and Nguyen Environ Syst Res (2015) 4:24 Page 13 of 13

Nguyen is graduate student in Faculty of Engineering, Vietnam National Ibbini MS et al (2014) Simscape solar cells model analysis and design. In:
University of Agriculture. Zaharim A, Sopian K, Bulucea A, Niola V, Skala V (eds) 8th International
conference on renewable energy sources (RES’ 14), 2nd International con-
Author details ference on environmental informatics (ENINF’14), Kuala Lumpur, Malaysia,
1
 Department of Electric Power System, Faculty of Engineering, Viet- 23–25 April 2014, WSEAS Press
nam National University of Agriculture, Trau Quy town, Gia Lam district, Jena C, Das A, Paniigrahi CK, Basu M (2014) Modelling and simulation of
Hanoi 10000, Vietnam. 2 Faculty of Engineering, Vietnam National University photovoltaic module with buck-boost converter. Int J Adv Eng Nano
of Agriculture, Trau Quy town, Gia Lam district, Hanoi 10000, Vietnam. Technol 1(3):4
Mohammed SS (2011) Modeling and simulation of photovoltaic module using
Acknowledgements Matlab/Simulink. Int J Chem Environ Eng 2(5):6
The authors are grateful to the support by this work through the project Pandiarajan N, Muthu R (2011) Mathematical modeling of photovoltaic mod-
“Study, design and manufacture a solar PV system using SPV technology ule with Simulink. International Conference on Electrical Energy Systems
served for chicken farms in Faculty of Animal Science, Vietnam National Uni- (ICEES 2011), p 6
versity of Agriculture”, Vietnam (2014–2017). Panwar S, Saini RP (2012) Development and simulation photovoltaic model
using Matlab/Simulink and its parameter extraction. International confer-
Competing interests ence on computing and control engineering (ICCCE 2012)
The authors declare that they have no competing interests. Salmi T, Bouzguenda M, Gastli A, Masmoudi A (2012) Matlab/simulink based
modelling of solar photovoltaic cell. Int J Renew Energy Res 2(2):6
Received: 23 July 2015 Accepted: 15 October 2015 Savita Nema RKN, Agnihotri Gayatri (2010) Matlab/Simulink based study of
photovoltaic cells/modules/arrays and their experimental verification. Int
J Energy Environ 1(3):14
Sudeepika P, Khan GMG (2014) Analysis of mathematical model of PV cell
module in Matlab/Simulink environment. Int J Adv Res Electr Electr
Instrum Eng 3(3):7
References Tu H-LT, Su Y-J (2008) Development of generalized photovoltaic model using
Banu I-V, Istrate M (2012) Modeling and simulation of photovoltaic arrays. MATLAB/SIMULINK. Proc World Congr Eng Comput Sci 2008:6
World energy system conference, p 6 Varshney A, Tariq A (2014) Simulink model of solar array for photovoltaic
Gonzalez-Longatt FM (2005) Model of photovoltaic module in Matab. II power generation system. Int J Electr Electr Eng 7(2):8
CIBELEC 2005, vol 2005, p 5 Venkateswarlu G, Raju PS (2013) Simscape model of photovoltaic cell. Int J Adv
Gow JA, Manning CD (1999) Development of a photovoltaic array model for Res Electr Electr Instrum Eng 2(5):7
use in power-electronics simulation studies. IEEE Proc Electr Power Appl Walker G (2001) Evaluating MPPT converter topologies using a Matlab PV
146(2):8 model. J Electr Electr Eng Aust 21(1):7