CN112528211A - Method for fitting solar cell IV curve - Google Patents

Method for fitting solar cell IV curve Download PDF

Info

Publication number
CN112528211A
CN112528211A CN202011490703.9A CN202011490703A CN112528211A CN 112528211 A CN112528211 A CN 112528211A CN 202011490703 A CN202011490703 A CN 202011490703A CN 112528211 A CN112528211 A CN 112528211A
Authority
CN
China
Prior art keywords
fitting
data
curve
solar cell
function
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011490703.9A
Other languages
Chinese (zh)
Other versions
CN112528211B (en
Inventor
毛翌春
朱炬
年夫来
周康
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Cetc Instrumentation Anhui Co ltd
Original Assignee
Cetc Instrumentation Anhui Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Cetc Instrumentation Anhui Co ltd filed Critical Cetc Instrumentation Anhui Co ltd
Priority to CN202011490703.9A priority Critical patent/CN112528211B/en
Publication of CN112528211A publication Critical patent/CN112528211A/en
Application granted granted Critical
Publication of CN112528211B publication Critical patent/CN112528211B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • General Physics & Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Algebra (AREA)
  • Computational Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Testing Of Individual Semiconductor Devices (AREA)
  • Photovoltaic Devices (AREA)

Abstract

The invention discloses a solar cell IV curve fitting method, and belongs to the field of IV curve fitting. The method for obtaining the solar cell IV fitting curve realizes fitting of discrete I, V data originally collected by the solar cell through a plurality of sections of composite functions, obtains corresponding fitting functions through setting different data sections, and then combines the functions together to obtain an integrated smooth I-V data curve, thereby not only obtaining a test curve with high fitting precision, but also effectively avoiding the phenomenon of fitting curve oscillation.

Description

Method for fitting solar cell IV curve
Technical Field
The invention belongs to the field of IV curve fitting, and particularly relates to a solar cell IV curve fitting method.
Background
The model of the solar cell is a typical nonlinear transcendental function, and the characteristic parameters of the solar cell are difficult to solve through the function. The characteristic parameters of a solar cell are analyzed by engineering by collecting (I, V) data for a set of solar cells. However, the (I, V) data collected by the test circuit is a discrete set of test data, and noise data is inevitably included in the collected data due to the noise effect of the test circuit, resulting in deviation and fluctuation of the test data. In order to obtain a smooth I-V data curve, a curve fitting method is usually used to obtain a smooth curve, so that the smooth curve can be optimally fitted to the acquired data, and the fluctuation and noise of the acquired data can be avoided and reduced. The fitting algorithms adopted at present comprise an exponential function analytical method, a pseudo-Monte Carlo algorithm and the like, but because the algorithms have certain defects, or the fitting precision is low, or the algorithms are complex and the like, a polynomial fitting algorithm is mostly adopted in engineering application, but the polynomial fitting algorithm has the defects that the fitting precision is related to the polynomial order, and meanwhile, the phenomenon of fitting curve oscillation easily occurs.
Disclosure of Invention
In order to solve the problems, the invention provides a method for fitting the solar cell IV curve, which is reasonable in design, overcomes the defects of the prior art and has a good effect.
In order to achieve the purpose, the invention adopts the following technical scheme:
a fitting method of a solar cell IV curve comprises the steps of conducting segmentation analysis on a group of discrete solar cell I, V data collected by a test circuit, conducting polynomial fitting on different data segments to obtain corresponding polynomial fitting functions, integrating the polynomial fitting functions into a group of fitting functions, and finally obtaining a high-precision smooth I-V curve by utilizing the fitting function group.
Preferably, the method specifically comprises the following steps:
step 1: collecting discrete data of the solar cell by using a test circuit, wherein the number of collected points is n;
step 2: calculating I in collected datai*ViDividing the whole data acquisition segment into two segments from a first acquisition point to an mth acquisition point, and from an m +1 th acquisition point to an nth acquisition point, wherein i is 1-n;
and step 3: calculating a coefficient set f [ 1-15 ] of a 14-order polynomial fitting function according to the acquired data of the full data segment 1-n to obtain a fitting function 1, as shown in formula (1):
Figure BDA0002840600390000011
wherein, ViRepresenting the voltage value, Inew, in the raw acquired dataiRepresenting the current value after 14-order polynomial fitting, wherein the value range of i is 1-n;
and 4, step 4: calculating a coefficient set h [ 1-15 ] of a 14-order polynomial fitting function corresponding to the data segment according to the data of the first segment of acquired data segment to obtain a fitting function 2, as shown in formula (2):
Figure BDA0002840600390000021
wherein, VjRepresenting the voltage value, Inew, in the raw acquired datajRepresenting the current value after 14-order polynomial fitting, wherein the value range of j is 1-m;
and 5: the fitting function 1 and the fitting function 2 are combined into a function group, data fitting is carried out on the function group corresponding to different data sections respectively, and finally a group of complete (Inew, V) fitting data is integrated, wherein the function group is shown as a formula (3):
Figure BDA0002840600390000022
wherein, InewiFor the fitted current value, ViThe collected voltage value is m, the position of a data segmentation point is m, and the total number of collected data points is n;
step 6: and finally obtaining the solar cell I-V curve with high precision and smoothness according to the fitting function group.
The invention has the following beneficial technical effects:
the method realizes the fitting of discrete (I, V) data originally collected by the solar cell through a multi-section composite function, obtains corresponding fitting functions through setting different data sections, and then combines the functions together to obtain an integrated smooth I-V data curve; the invention can not only obtain the test curve with high fitting precision, but also effectively avoid the phenomenon of fitting curve oscillation.
Drawings
FIG. 1 is a schematic diagram of an equivalent circuit of a solar cell according to the present invention;
FIG. 2 is a diagram of originally collected discrete data according to example 1 of the present invention;
FIG. 3 is a schematic diagram of a 14 th order polynomial fit curve according to example 1 of the present invention;
FIG. 4 shows full data segment data according to example 2 of the present invention;
FIG. 5 is a schematic diagram of oscillation occurring by using a single 14 th order polynomial fitting curve according to embodiment 2 of the present invention;
FIG. 6 is a diagram illustrating a comparison between a multi-segment complex function fitted curve and a 14 th order polynomial fitted curve according to example 2 of the present invention.
Detailed Description
The method adopts a multi-section complex function method to fit the originally acquired (I, V) data to obtain a high-precision smooth I-V curve, and effectively solves the oscillation problem of the fitted curve.
The invention is described in further detail below with reference to the following figures and detailed description:
fig. 1 shows an equivalent circuit of a solar cell, and fig. 1 shows an equivalent circuit in the case of illumination.
Under the condition of illumination, the volt-ampere characteristic equation of the solar cell model is shown as the formula (4):
Figure BDA0002840600390000031
wherein, I and V are the output current and the output voltage of the solar cell respectively; i isLIs a photo-generated current; i isdIs a dark current; i is0Is a diode reverse saturation current; q is the electron charge constant; k is the boltzmann constant; t is the junction temperature at test time; rsIs a series resistor; rshIs a parallel resistor.
Example 1
The (I, V) data set of solar cells under light conditions collected by the test circuit is a discrete set of data, as indicated by the dots in fig. 3. In order to obtain a smooth I-V curve, polynomial fitting is usually used in engineering to achieve the curve fit. Equation (5) below is a 14-order polynomial fitting function for the discrete data points in fig. 3 (the 14-order polynomial fitting function is used because the 14-order fitting accuracy substantially meets the engineering requirements), and the fitting curve is as shown in fig. 3. The collected (I, V) data and the 14 th order polynomial fitted current data are listed in the table in fig. 2.
Figure BDA0002840600390000032
Wherein, ViThe value is the originally collected voltage value, Inewi is the fitted current value, and i is the total number of collected points from 1.
It can be seen from fig. 3 that the error between the 14 th-order fitting curve and the original collected data is small, and a smooth I-V fitting curve of the solar cell can be obtained.
Example 2
However, for some (I, V) data acquisitions, fitting the function with only a single polynomial tends to cause oscillations in the fitted curve, as shown in FIG. 5. In fig. 5, the dots are the original collected data, and the curve is a fitting curve obtained by applying a 14-order polynomial fitting function to the data of the whole data segment, and the fitting function is shown in formula (6). The collected (I, V) data and the 14 th order polynomial fitted current data are listed in the table in fig. 4.
Figure BDA0002840600390000033
Wherein, ViThe value is the originally collected voltage value, Inewi is the fitted current value, and i is the total number of collected points from 1.
It can be seen from fig. 5 that the curves fitted with a single polynomial show oscillations.
The curve obtained by piecewise fitting the raw collected data in fig. 4 by the multi-segment complex function fitting method of the present invention is shown in fig. 6. In fig. 6, the diamond points are the original collected data, the smooth curve is the fitted curve obtained by fitting the full data segment data with a 14-order polynomial fitting function, and the curve with points is the fitted curve obtained by fitting the full data segment data with a multi-segment complex function. The fitting function group of the multi-section complex function fitting method is shown as the formula (7):
Figure BDA0002840600390000041
where m is the position of the segmentation point, which has a value of 54; n is the total number of collected points, and the value is 399.
As can be seen from FIG. 6, the method of the present invention uses a multi-segment complex function to fit the original sampling data, so that the oscillation phenomenon of the fitting curve can be effectively solved, and the fitting accuracy is improved.
It is to be understood that the above description is not intended to limit the present invention, and the present invention is not limited to the above examples, and those skilled in the art may make modifications, alterations, additions or substitutions within the spirit and scope of the present invention.

Claims (2)

1. A fitting method of a solar cell IV curve is characterized in that a group of discrete solar cell I, V data collected by a test circuit is subjected to sectional analysis, polynomial fitting is carried out on different data sections to obtain corresponding polynomial fitting functions, the polynomial fitting functions are integrated into a group of fitting functions, and a high-precision smooth I-V curve is finally obtained by utilizing the group of fitting functions.
2. The method for fitting the solar cell IV curve according to claim 1, comprising the following steps:
step 1: collecting discrete data of the solar cell by using a test circuit, wherein the number of collected points is n;
step 2: calculating I in collected datai*ViThe point m with the maximum value divides the whole data acquisition segment into two segmentsThe first section is from the first acquisition point to the mth acquisition point, and the second section is from the m +1 th acquisition point to the nth acquisition point, wherein i is 1-n;
and step 3: calculating a coefficient set f [ 1-15 ] of a 14-order polynomial fitting function according to the acquired data of the full data segment 1-n to obtain a fitting function 1, as shown in formula (1):
Figure FDA0002840600380000011
wherein, ViRepresenting the voltage value, Inew, in the raw acquired dataiRepresenting the current value after 14-order polynomial fitting, wherein the value range of i is 1-n;
and 4, step 4: calculating a coefficient set h [ 1-15 ] of a 14-order polynomial fitting function corresponding to the data segment according to the data of the first segment of acquired data segment to obtain a fitting function 2, as shown in formula (2):
Inewj=h(1)+h(2)*Vj+h(3)*Vj 2+h(4)*Vj 3+h(5)*Vj 4+h(6)*Vj 5+h(7)*Vj 6+h(8)*Vj 7+h(9)*Vj 8+h(10)*Vj 9+h(11)*Vj 10+h(12)*Vj 11+h(13)*Vj 12+h(14)*Vj 13+h(15)*Vj 14 (2)
wherein, VjRepresenting the voltage value, Inew, in the raw acquired datajRepresenting the current value after 14-order polynomial fitting, wherein the value range of j is 1-m;
and 5: the fitting function 1 and the fitting function 2 are combined into a function group, data fitting is carried out on the function group corresponding to different data sections respectively, and finally a group of complete (Inew, V) fitting data is integrated, wherein the function group is shown as a formula (3):
Figure FDA0002840600380000012
wherein, InewiFor the fitted current value, ViThe collected voltage value is m, the position of a data segmentation point is m, and the total number of collected data points is n;
step 6: and finally obtaining the solar cell I-V curve with high precision and smoothness according to the fitting function group.
CN202011490703.9A 2020-12-17 2020-12-17 Method for fitting solar cell IV curve Active CN112528211B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011490703.9A CN112528211B (en) 2020-12-17 2020-12-17 Method for fitting solar cell IV curve

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011490703.9A CN112528211B (en) 2020-12-17 2020-12-17 Method for fitting solar cell IV curve

Publications (2)

Publication Number Publication Date
CN112528211A true CN112528211A (en) 2021-03-19
CN112528211B CN112528211B (en) 2022-12-20

Family

ID=75000809

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011490703.9A Active CN112528211B (en) 2020-12-17 2020-12-17 Method for fitting solar cell IV curve

Country Status (1)

Country Link
CN (1) CN112528211B (en)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3622770A (en) * 1968-08-21 1971-11-23 Hughes Aircraft Co Straight line segment function generator
CN103633937A (en) * 2013-12-10 2014-03-12 南昌航空大学 Solar cell parameter extracting method based on Lambert W function and polynomial fitting
CN104484703A (en) * 2014-12-30 2015-04-01 合肥工业大学 Sigmoid function fitting hardware circuit based on Remez approximating algorithm
CN104716903A (en) * 2015-03-17 2015-06-17 上海交通大学 Automatic segmentation polynomial fitting method applicable to photovoltaic array output characteristic curve
CN106199483A (en) * 2016-06-13 2016-12-07 电子科技大学 A kind of nonlinear error compensation method of digital multimeter
CN107612656A (en) * 2017-09-18 2018-01-19 北京航空航天大学 A kind of Gaussian approximation method for simplifying suitable for polarization code
WO2018112881A1 (en) * 2016-12-23 2018-06-28 深圳中兴力维技术有限公司 Rapid prediction method for battery charging performance and system thereof
CN111623898A (en) * 2020-06-24 2020-09-04 广东纽恩泰新能源科技发展有限公司 NTC thermistor temperature and resistance detection method and temperature sensor detection method

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3622770A (en) * 1968-08-21 1971-11-23 Hughes Aircraft Co Straight line segment function generator
CN103633937A (en) * 2013-12-10 2014-03-12 南昌航空大学 Solar cell parameter extracting method based on Lambert W function and polynomial fitting
CN104484703A (en) * 2014-12-30 2015-04-01 合肥工业大学 Sigmoid function fitting hardware circuit based on Remez approximating algorithm
CN107247992A (en) * 2014-12-30 2017-10-13 合肥工业大学 A kind of sigmoid Function Fitting hardware circuits based on row maze approximate algorithm
CN104716903A (en) * 2015-03-17 2015-06-17 上海交通大学 Automatic segmentation polynomial fitting method applicable to photovoltaic array output characteristic curve
CN106199483A (en) * 2016-06-13 2016-12-07 电子科技大学 A kind of nonlinear error compensation method of digital multimeter
WO2018112881A1 (en) * 2016-12-23 2018-06-28 深圳中兴力维技术有限公司 Rapid prediction method for battery charging performance and system thereof
CN107612656A (en) * 2017-09-18 2018-01-19 北京航空航天大学 A kind of Gaussian approximation method for simplifying suitable for polarization code
CN111623898A (en) * 2020-06-24 2020-09-04 广东纽恩泰新能源科技发展有限公司 NTC thermistor temperature and resistance detection method and temperature sensor detection method

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
CHEN Y: "Parameters extraction from commercial solar cells I-V characteristics and shunt analysis", 《APPLIED ENERGY》 *
CHEN Y: "Parameters extraction from commercial solar cells I-V characteristics and shunt analysis", 《APPLIED ENERGY》, vol. 88, no. 15, 31 December 2011 (2011-12-31) *
任明凡等: "太阳电池I-U特性曲线的拟合", 《安徽工业大学学报(自然科学版)》 *
任明凡等: "太阳电池I-U特性曲线的拟合", 《安徽工业大学学报(自然科学版)》, no. 02, 30 April 2007 (2007-04-30) *
肖文波: "基于Lambert W函数提取太阳电池五参数的方法及光强对参数影响的规律", 《科学通报》 *
肖文波: "基于Lambert W函数提取太阳电池五参数的方法及光强对参数影响的规律", 《科学通报》, vol. 59, no. 36, 31 December 2014 (2014-12-31) *

Also Published As

Publication number Publication date
CN112528211B (en) 2022-12-20

Similar Documents

Publication Publication Date Title
CN110956312B (en) Photovoltaic power distribution network voltage prediction method based on EMD-CNN deep neural network
CN110648249B (en) Annual power balance measuring and calculating method, device and equipment
CN113496311A (en) Photovoltaic power station generated power prediction method and system
CN111582555B (en) Photovoltaic power prediction method based on foundation cloud picture image characteristics
WO2023201552A1 (en) County-wide photovoltaic prediction method based on cluster division and data enhancement
JP6355491B2 (en) Pyrometer performance degradation state estimation device, psoriometer performance degradation state estimation system
CN112528211B (en) Method for fitting solar cell IV curve
CN115689055A (en) Short-term solar irradiance prediction method and device
CN107944466B (en) Rainfall deviation correction method based on segmentation idea
CN109193791B (en) Wind power convergence tendency state-based quantification method based on improved shape value
US20140122003A1 (en) Prediction method for sun-tracking type photovoltaic system
CN108876060B (en) Big data based prediction method for wind power output probability of sample collection
CN113344293B (en) Photovoltaic power prediction method based on NCA-fusion regression tree model
CN115600503A (en) Multi-junction compound battery model parameter optimization method based on simplex-simulated annealing method
CN106709587B (en) Direct radiation prediction method based on conventional weather forecast
CN111666691B (en) Statistical method for atmospheric optical turbulence parameters
CN114663781A (en) Carbon source inversion method for urban multi-source environmental data
CN111596384B (en) Inclined plane radiation prediction method based on weather type effective identification
CN109299431B (en) Photovoltaic module performance consistency evaluation method based on Hilbert marginal spectrum characteristics
CN112070319B (en) Ultra-short-term photovoltaic power generation power prediction method based on time sequence dynamic regression
Xia et al. Monthly calibration and optimization of Ångström-Prescott equation coefficients for comprehensive agricultural divisions in China
CN105975803A (en) Engineering model of high-precision silicon solar cell and calculating method
CN116260141B (en) Reconstruction method and system for power of photovoltaic power station and reconstruction terminal
Liu et al. A Self-Data-Driven Method for Lifetime Prediction of PV Arrays Considering the Uncertainty and Volatility
CN115900936A (en) Horizontal plane total radiation and scattered radiation prediction method based on weather classification

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information

Address after: 233010 no.1300, Yinghe Road, Yuhui District, Bengbu City, Anhui Province

Applicant after: CLP kesiyi Technology (Anhui) Co.,Ltd.

Address before: 233010 no.1300, Yinghe Road, Yuhui District, Bengbu City, Anhui Province

Applicant before: CETC INSTRUMENTATION (ANHUI) Co.,Ltd.

CB02 Change of applicant information
GR01 Patent grant
GR01 Patent grant