International Conference on Renewable Energies and Power Quality (ICREPQ’17)

Malaga (Spain), 4th to 6th April, 2017

Renewable Energy and Power Quality Journal
(RE&PQJ) ISSN 2172-038 X, No.15 April 2017

The Newton Raphson Method in the Extraction of Parameters of PV Modules

L. R. D. Reis, J. R. Camacho, D. F. Novacki.

School of Electrical Engineering

Universidade Federal de Uberlândia
(UFU)
Campus Santa Mônica- Av. João Naves de Ávila, 2121. ZipCode:38400.902 – Uberlândia – MG - Brazil
Phone: +55 34 3239 4734, email: [email protected], [email protected], d
[email protected].

Abstract. This research is justified by the expansion of electricity; And (ii) heliothermic, which consists of the
electrical systems and the large-scale adoption of small generation of electricity through the heating of a fluid
electric power generation systems, in this particular case and will function in a similar way to a conventional
the semiconductor photovoltaic(PV) systems connected to thermoelectric plant [3]. Among the two types of solar
the grid. In order to analyse these systems it is imperative power plants, photovoltaic systems have proved to be
that the parameters of these generation systems are of more viable and therefore more competitive.
reasonable reliability. This research aims to study The efficiency of photovoltaic systems is influenced
mathematical methods to obtain parameters of by three questions: the efficiency of the photovoltaic
photovoltaic modules for static and dynamic studies in module, which depends on the way it is made, the
medium voltage electrical energy systems. Therefore, it is materials used and the constructive aspects of the
necessary to use equations that characterize this converter used, and lastly the efficiency of the algorithm,
generation system, depending on the parameters of which controls and traces the point of maximum power
interest. These model parameters are: generated photo for the system, the maximum power point tracking
current, reverse saturation current, ideality factor, series (MPPT) algorithm.
resistance and parallel resistance. These parameters are Improving the efficiency of the module and the
of fundamental importance for the simulation of this converter is a more labor-intensive task, which requires
generation system and the reconstruction of V.I and P.V research and development by the manufacturers,
curves of photovoltaic systems in power systems. Other however, improving the efficiency of the MPPT
genetic algorithm methods have already been used in this algorithm is easier and cheaper [4].
research in the past to obtain these parameters when To work on the development and improvement of the
totally unknown [1], [2]. The Newton Raphson method is control algorithm it is necessary to first extract the
a mathematical method widely used in power systems parameters of the module or the photovoltaic system that
and can be used when one has a rough idea of the one wishes to work. Therefore, this article aims to obtain
solution. This article clearly shows the possibility of a method that is able to extract such parameters through
using this fast method. traditional fast numerical methods applied to engineering
problems, able to characterize the behavior of the module
Keywords. Parameters extraction, Newton Raphson in working conditions.
Method, Photovoltaic Module, Modeling. The extraction of parameters of a photovoltaic
module has already been studied by another researcher at
1. Introduction the Center of Research and Extension in Alternative
Energies (NUPEA). Although, the research developed in
Among all renewable energy sources, the most abundant is this case used a genetic algorithm based method, the
solar energy. This energy reaches the earth through the Trust-Region Optimization Method [1], [2], since the
solar rays, which are electromagnetic waves that carry parameters were completely unknown.
light and heat necessary for terrestrial life.
The conversion of solar energy into electrical energy 2. PV Module Modelling
occurs through the solar irradiation that reaches the surface
of the earth. The energy transformation can be harnessed Some technical characteristics of the photovoltaic module
through two different types of power plants: (i) are provided by the manufacturer, they are: The open
photovoltaic, which consists of the arrangement of circuit voltage VOC, obtained from the measurement of
photovoltaic modules installed in a relatively large area, voltage in the module when it has no load; the current of
usually made of semiconductors (silicon, germanim, the module when the device is shorted I SC (short circuit
etc…), capable of converting solar irradiation directly into
https://doi.org/10.24084/repqj15.416 1 RE&PQJ, Vol.1, No.15, April 2017
current) and the maximum power (PMP) of a photovoltaic
module, reached when are obtained the maximum power In the case of PV modules with various series connected
current (IMP) and the maximum power voltage (V MP)[5], cells, equation 1 become:
illustrated in Figure 1.
q(V+I Rs) V+I
I = IL − ID (e Nsm K T − 1) − Rs (2)
Rp

Being Ns the number of series connected cells.

The reverse saturation current of the diode (I S) is
produced by the component when the cell is not under
sunlight, i.e. it is in the dark. It depends on the local
temperature and is generated thermally.
The series and parallel resistors presented are
responsible for internal dissipations in the cell. The series
resistance (RS) originates in the resistance of the
semiconductor material itself, in the metal contacts and in
metal-semiconductor junctions, characterizes the
material's ohmic losses, the series resistance contributes
to reduce the Isc of the cell, but practically does not affect
Fig, 1: Tipical V-I curve of a PV cell[6]. its Voc.
The parallel resistance (Rsh) is caused by impurities
A. Equivalent Circuit of a Cell Photovoltaics: Single and defects in the structure, especially near the edges,
Diode Model which produce an internal path for a leakage current, thus
reducing the current effectively produced by the device,
A photovoltaic cell is simply a P-N junction so Rsh consists of reducing VOC without influence in ISC
semiconductor. When exposed to sunlight, a current is [6].
generated and varies proportionally to the irradiance at that
time [7]. Figure 2 below shows the equivalent circuit of a 3. The Newton Raphson Method
photovoltaic cell with one diode.
The Newton Raphson Method is an iterative method that
consists in estimate of a given function f(x) with an initial
guess.
The method is obtained through the Taylor series
expansion in (x − xO) given below.
1 2
f (x)  f (x )  f '(x )(x  x )  f ''(x )(x  x )  ...  0
0 0 0 0 0
2

Suppose the initial guess is very close to the real root

of the equation. Then (x-x0) is small enough and only the
first terms are important for estimating the value of the
root, given x0. By truncating the series in the second term
Fig. 2: Equivalent circuit f a PV cell with one diode [8]. we obtain the general formula of Newton Raphson's
method [9].
Equation describing the circuit is given by:
f (x 0 )
q (V+IRs) V + IRs x1  x0  (3)
I = Iph − ID (e mKT − 1) (1) f ' (x0
− Rsh
)
where:
Therefore, given xn, the point x(n+1) will be obtained
I - cell current (A), by intersecting the tangent line at f(x) in xn with the x
Iph - photo generated current (A), axis. You can illustrate this mathematically, according to
ID - diode reverse saturation current (A), equation (3) and graphically, according to figure (3).
q - electron elementary charge (1,602 ×
10–19coulombs), Rshor Rp o voltag
K – the Boltzmann constant (1,381 × 10–23 – cell u e (V).
J/K), shunt t
T - cell temperature (K), resistanc p
m – diode ideality fator, e (Ω), u
Rs – cell series resistance (Ω), V – cell t
https://doi.org/10.24084/repqj15.416 2 RE&PQJ, Vol.1, No.15, April 2017
 xn1  x n
 f (x n )
f '(x ) (
n
4
)

https://doi.org/10.24084/repqj15.416 3 RE&PQJ, Vol.1, No.15, April 2017

 relationship of voltage and power, highlighting the point
of maximum power.

Figure 3: Illustration of the Newton Raphson Method. Fig. 4: Curve 1 – Relationship of voltage and current.

The convergence of the Newton Raphson Method is

guaranteed for a certain interval [a,b] containing the root
of f(x), provided that f(x) and f'(x) are continuous in this
interval and that f(α) = 0, where α is the root of f(x)[10].

4. Problem Formulation
The data were obtained through measurements performed
on a photovoltaic module. The technical characteristics of
the module used are shown in Table I.

Table I. - Feature of the PV Module

Fig. 5: Curve 1 – Relationship of voltage and power.
Value
Power 40 W
VNp 16.6 V
INp 2.45 A
Voc 20.5 V
IS€ 2.80 A
Ns 36
Material Type monocrystalline silicon

The voltage and current curves of the photovoltaic

module were obtained through the mini KLA meter. The
meter operates with a temperature range of -20 to 100° C
and an accuracy of ∓0.4%. In addition, the meter provides
the internal temperature of the module and the solar Fig. 6: Curve 2 – Relationship of voltage and current
irradiance to which it is subjected.
Two tests were carried out with the same photovoltaic
module at different times. The data of the two curves
obtained are shown in Table II.

Table II. - Used curves.

Irradiance Temperature
Curve 1 225 W/m² 25°C
Curve 2 596.8W/m² 35°C

With the knowledge of data to obtain the curves, a

program was developed in the Octave computational tool
to generate the two curves, according to the figures below.
Figures four and six show the relationship of voltage and Fig. 7: Curve 2 – Relationship of voltage and power.
current obtained, highlighting the five points used in the
developed program. Figures 5 and 7 illustrate the In order to apply the Newton Raphson numerical
method to the proposed problem it is necessary to use
several points of the sampled curve. In the case of a diode
model, it is necessary to extract five parameters, and it is essential to formulate five equations so that the solution of
  (V1  I1 x4 )  (V1  I1 x4 )
  (V1 1 4I x ) 6

the problem
1 becomes possible. Therefore, it10x (V  I x211)
is necessary
36 x3
to4e being close 10
to V I x. 
6 OC
36 x3
e second
The and Ifourth  I11chosen
Vwere x4 
curve fivepoints 1 2 1
36 x3
extract from the10 real
6 (e 1) that will be used in
36 x 32 arbitrarily, considering 36 x(I32, V2) = L/2 xand (I4, V4x) 2= 3.L/2,
solving the problem. In this work, such points were used:  (V2  I so 2 x4 )that L represents the total
5 5 
 (V2 number
I 2 x4 ) of points of the real
The firstcorresponds to the 2 2 4 initial point of the curve, being
 
10  6 x (V 4  I x )  e 36curve.
(V I x ) 
x3 The five
10I xpoints
6
e sampled are Ihighlighted
36 x 3
V 22I 4inx VxI
I1 near ISC1; The
10 third
6 (e corresponds
36 x3
 1)to the maximum222power 22
curves of figures (4) and (6). Through  2 equation 2, it is
point PMP, with I3 = IMP and V3 = VMP; The36 x 32
fifth then possible to generate 36 x 3the five equations x5 x 52 in the
used 
  (V 3 3I4x )
system, as shown below. 36 x3
 (V  3I 3x4)
J 1   VIx 
 (V3  I 3 x 4 )
36 x3 10  6 x (V  I233x ) e 36 x3
10  6 I x e I
1 10  6 (e  1) 4
3 2 3 33 4 

 36 x 32 36 x 3 x5 x 2 5
 (V4  I 4 x 4 )
  (V44 4I x )
 (V4  I 4 x4 ) 
10  6 x (V  I244x ) e 36 x3
10  6 I x e 36 x3 I VI x 
 10  6 (e
1
36 x3
 1) 4
4 2 4 44 4 
36 x 32 36 x x x
 52
  (V5  I 5 x4 )
3 5
 (V5  I 5 x4 )
  (V55 4I x )
10  6 x (V 255 I 4x )  e 36 x3 VIx
10  6 I x e 36 x3 I
1
 10  6 (e 36 x3
 1)   
5 2 5 55 4
 36 x 32 36 x 3 x5 x 52 

corresponds to the last sampled point of the curve, V5

Fig. 8: Inverse of the Jacobian Matrix

Being the inverse of the Jacobian matrix given by

 (V1  I1x4 ) (V  I x ) Figure 8, where,
f (x)  x 106 x (e 36  1)  1 1 4  I  0
x3
1 1 2 1
x5 x I 
 1 p
 (V2  I 2 x4 ) (V  I x )    h 
1)  2 24
 I  02 x I
f 2 (x)  x1 10 x 2(e 6 36
x3 x5 q  2  d 
 and x    m  .

T  x3   R 
 1)  (V 3 I x3 4) I  0
 (V3  I3 x4 )
f (x)  x 106 x (e 36
x3 4 s
3 1 2
1   
x5
 (V  I x )   sh 
 R
 x
 5

444 (V  I x )
f 4 (x)  x1 106 x 2(e 36
x3
1)  4 44
 I 4 0 5. Results Analysis
x5
 (V5  I 5 x4
6 ) 36 x3 (V5  I 5 x4 ) After the tests made un the program developed using
f 5(x)  x 10
1 x (e2 1)   I 5 0 Octave, the following results were obtained:
x5
Table III. - Initial Guess Used
Given the equations, it is possible to assemble the INITIAL GUESS (IG)
system to solve the problem. 1º x0 = [1.5; 0.5; 1.0; 0.5; 50]
1
2º x0 = [2.0; 0.1; 1.0; 1.0; 100]
 f1 f1 f1 f1 f  3º x0 = [2.0; 1.0; 1.0; 1.0; 100]
x x2 x3 x4 1x5
4º x0 = [2.0; 1.0; 1.0; 1.0; 50]
 1 
 x1,k 1  x1, k   f2 f2 f2 f2 f2   f1 ( xk )  5º x0 = [0.6; 2.5; 1.5; 2.0; 100]
 
x  x   x1 x2 x3 x4 x5   f (x )
 f f f f f   
2k
During the simulations executed several guesses were
 x    
2, k 1 2,k
x  3 3 3 3 3
 . f ( x )
3k
3,k 1 x x x  used, some of them converged to solve the problem,
x  3, k
 x   x
1
x 2 3 4 5 f (x ) others did not. From those that converged, shown in
 4, k 1   
4,k f4 f4 f4
f4
 x f4x   4 k     f (x ) Table III, it was observed that the generated values were
 5, k  5,k  x 2 x3 x4 identical. In addition, it was noted that the initial guess
 5   5 k
1  x 
x  f1
 x 5 f5 f5 f5 f5  when closer to the root generated fewer iterations, as
x x x x
 1 2 3 4 5  shown in Tables IV and V.
 Table IV. - Obtained Parameters of the Curve 1
CURVE 1
IG Iph(A) Is.10-6 (µA) m Rs (Ώ) Rsh(Ώ) Iters
1º 0.6687 2.006 1.4288 1.1686 120.58 9
2º 0.6687 2.006 1.4288 1.1686 120.58 7
3º 0.6687 2.006 1.4288 1.1686 120.58 9
4º 0.6687 2.006 1.4288 1.1686 120.58 9
5º 0.6687 2.006 1.4288 1.1686 120.58 4

Table V. - Obtained parameters of curve 2

CURVE 2
IG Iph(A) Is. 10-6 (µA) m Rs (Ώ) Rsh(Ώ) Iters
1º 1.8462 0.0441 1.0245 0.1979 58.92 5
2º 1.8462 0.0441 1.0245 0.1979 58.92 7 Fig. 11: Curve P x V, T=25°C e G=225.0 W/m².
3º 1.8462 0.0441 1.0245 0.1979 58.92 8
4º 1.8462 0.0441 1.0245 0.1979 58.92 8

With the results obtained in tables V and VI, it was

possible to simulate the corresponding curves. Figures 9
and 10 show a relationship of current and voltage,
comparing a real curve with a simulated curve. It is
possible to observe that the simulated curves go through
the five points used in the problem formulation, thus
obtaining a good approximation of the real curve.
Figures 11 and 12 show the relation of power and
voltage. It is possible to observe that the simulated curves
were able to reach the points of maximum power and that
their shape is close to the real curve.
Fig. 12: Curve P x V, T=35°C e G=596.8 W/m².

When using the Newton Raphson method, it is

necessary to have a knowledge of the approximate value
that one wishes to obtain, since the initial guess should be
close to the solution of the problem, otherwise the
method is not able to find the solution. This knowledge is
necessary because, the Jacobian matrix calculated in each
iteration must be nonsingular, that is, with the
determinant other than zero.
The study of photovoltaic systems has evolved a lot in
the last years, this has generated in the researchers of the
subject an improvement of the knowledge, thus allowing
to estimate values close to the real parameters of
photovoltaic cells.
However, tests were performed with values distant
Fig. 9: Curve I x V, T=25°C e G=225.0 W/m². from the solution, as an example, when the initial guess n
curve 1 was x0 = [1.5; 1e-3; 1.5; 1.5; 100], and in curve 2
when x0 = [1.5; 1.0; 1.5; 1.5; 100], the result was not
obtained because the method did not show convergence.
Comparing the results obtained in this article with
those from [1], the parameters Iph (photogenerated
current), m (diode ideality factor) and Rsh (parallel
resistance) approximated to those found by that author.
The parameters Id (diode reverse current) and Rs (series
resistance) diverged considerably, as shown in Table VI.

Table VI. - Obtained Parameters fromTrust-Region Methods [7]

Fig. 10: Curve I x V, T=35°C e G=596.8 W/m².

 Curve Iph(A) Is.10-6 m Rs (Ώ) Rsh(Ώ)
(A)
1º 0.6736 2.4258 1.3390 2.0364 129.8172
2º 1.8498 0.1632 1.0411 0.5282 69.1844
6. Conclusions photovoltaic modules through the use of single and double
exponential models”, in International Conference on Renewable
The parameters from VxI curve using the method proposed Energies and Power Quality (ICREPQ), 2011.
in this work presented satisfactory results when compared
to the actual curve obtained from the photovoltaic module. [2] P. Rodrigues, “The extraction of parameter of photovoltaic
panels from the solution of a non-linear equation system
The Newton Raphson method can be used for various through thrust region techniques”. MSc Dissertation.
types of problems with nonlinear equations, being able to Universidade Federal de Uberlândia, Uberlândia, Minas Gerais,
cover several areas of science. However, it is fundamental Brazil, 2012, 117p. (In portuguese)
to have the knowledge of a range of values close to the
solution, otherwise convergence is not obtained, as was [3] FGV ENERGIA. “Complementary Renewable Energies”.
shown in this work. Energy Notebooks from Fundação Getúlio Vargas. December,
The results found in this work approximated those 2015, Year 2, N° 4.
obtained in [1], even using different methods, the
[4] Shivananda Pukhrem,“A Photovoltaic Panel Model in
divergences in the calculated parameters were small,
Matlab/Simulink”, Wroclaw University of Technology (2013).
allowing a comparison between the two studies.
The method is very reliable to get parameters of [5] R. Zilles, W. N. Macêdo, M. A. B. Galhardo and S. H. F.
modules when the user has some knowledge of the PV Oliveira, “Photovoltaic systems connected to the grid”,
system under investigation, and the approximated values Workshop Texts, 2012, 208p. (In portuguese)
of parameters to be used as an initial guess. The Thrust
Region Method algorithm used in the past is very useful [6] J.T. Pinho, and M.A. Galdino, “Photovoltaic Systems
when the user has no idea on the initial guess of Engineering Manual”, CEPEL-CRESESB, Rio de Janeiro
parameters. (2014), 530p. (In portuguese)
Finally, this work presents a simple and fast
[7] Jacob Jamesm Nedumgatt, Jayakrishnan K. B., Umashankar
mathematical method that can be applied in practical S., Vijayakumar D., “Perturb and Observe MPPT Algorithm for
electrical engineering problems of photovoltaic systems. Solar PV Systems-Modeling and Simulation”, in Proc. School
of Electrical Engineering, VIT University.
Acknowledgements
[8] F. Kininger, “Photovoltaic Systems Technology”, Kassel
The authors are grateful to CNPQ (Brazilian National University (2003), Germany, 155p.
Council for Scientific and Technological Development) for [9] M. D. Smith, “Newton-Raphson Technique”, 1998,
the financial support of this research and to UFU Electronic document available in:
(Universidade Federal de Uberlândia) for the support and http://web.mit.edu/10.001/Web/Course_Notes/NLAE/node6.ht
opportunity to conduct the research in this institution. ml. Accessed in October, 02, 2016.

References [10] E. F. Cota, A. F. Bastos, S. R. Silva e H. A. Pereira,

“Application of the Newton-Raphson Method in the Control of
[1] P. Rodrigues, J. R. Camacho, and F. B. Matos “The the External Resistance of Wind Generators During Voltage
application of Trust Region Method to estimate the parameters of Sags”, in Proc. IV SBSE (2012), Goiânia, Goiás, Brasil. (In
portuguese)