本发明属于光伏阵列最大功率跟踪研究技术领域，具体涉及一种基于bezier多项式的光伏阵列最大功率跟踪方法。
背景技术
光电转换效率不高是当前制约光伏发电技术进一步发展的主要瓶颈之一，为提高光伏阵列的发电效率，需以最大功率跟踪技术对不同工况下光伏阵列的输出进行控制，以保证工况变化时光伏电站能运行在最大功率点附近，提高光伏电站的发电能力，具有良好的现实意义。
国内外相关学者对光伏阵列的最大功率跟踪技术进行了较为深入的研究，并给出了一些比较典型的算法，按照控制策略不同，可分为按照某一规律预先给定电压或电流的数值，以使光伏阵列工作在最大功率点附近的静态模型跟踪策略；通过实时检测光伏阵列运行时电压和电流等参量变化，进而确定最大功率跟踪方向的动态模型跟踪策略；利用自动寻优算法对光伏阵列的最大功率进行优化搜索的基于智能算法的跟踪策略等。在上述典型的最大功率跟踪研究策略中，多以matlab软件为平台，通过编写s函数和simulink模型或以迭代计算等方式对光伏阵列进行建模，该过程不仅实现较为繁琐，也存在着计算时间长的弊端，提高了最大功率跟踪技术在工程中推广和使用的难度，阻碍了光伏发电技术的进一步发展。
技术实现要素：
针对现有技术中的上述不足，本发明提供的基于bezier多项式的光伏阵列最大功率跟踪方法解决了现有技术中实现光伏阵列最大功率跟踪实现过程繁琐，难以在工程中广泛推广等问题。
为了达到上述发明目的，本发明采用的技术方案为：基于bezier多项式的光伏阵列最大功率跟踪方法，包括以下步骤：
s1、测试光伏阵列的实际辐照度和温度值；
s2、计算实际测试辐照度和温度下，光伏模块的开路电压vocb，短路电流iscb和最大功率点(vmb，imb)；
其中，vmb为最大功率点电压，imb为最大功率点电流；
s3、基于开路电压vocb，短路电流iscb和最大功率点(vmb，imb)，利用bezier函数的光伏模块建模方法，得到给定条件下光伏模块的输出特性曲线；
s4、对输出特性曲线进行多项式拟合，并给出对应输出特性曲线的多项式拟合结果；
s5、利用电流测试装置检测光伏阵列的输出电流，同时利用电压检测装置检测光伏模块的输出电压；
s6、根据检测到的输出电压和电流，基于增量电导法，对多项式拟合结果得出的最大功率点进行逐一微调；
s7、将检测到的输出电压值作为自变量，使光伏模块工作输出特性在微调后的最大功率点电压处，使光伏模块工作在最大输出功率，实现光伏阵列的最大功率跟踪。
进一步地，所述步骤s2中：
最大功率点vmb的计算公式为：
vmb＝vmln(e+βδir)(1-γδt)
最大功率点电流imb的计算公式为：
短路电流iscb的计算公式为：
开路电压vocb的计算公式为：
vocb＝vocln(e+βδir)(1-γδt)
其中，δt为实际工况下的电池温度与标准测试条件下的电池温度差；
δir为实际工况下的辐照度与标准测试条件下辐照度的差；
α、β和γ均为拟合常数，且α＝0.0025/℃，β＝0.0005/(w/m2)，γ＝0.00288/℃；
irb为标准测试条件下的辐照度，且irb＝1000w/m2，
t标准测试条件下的电池温度，且t＝25℃；
vm、im是标准测试工况下光伏阵列最大功率点处的电压和电流；
isc、voc是标准测试工况下光伏阵列短路电流和开路电压；
ir是实际测试的辐照度。
进一步地，所述步骤s3具体为：
s31、建立平面直角坐标系，在建立的平面直角坐标系中确定短路电流点、最大功率点以及开路电流点的位置，并根据其位置，确定bezier函数的两个控制点在平面直角坐标系中的位置；
s32、根据所述短路电流点、最大功率点、开路电流点的位置，根据bezier函数特性，确定最大功率点左侧和右侧的bezier函数曲线上点的坐标；
s33、根据最大功率点左侧和右侧的点的坐标，绘制出光伏模块的输出特性曲线。
进一步地，所述步骤s4中的多项式拟合结果为一个以电压为自变量，以电流为函数的8次多项式。
进一步地，
所述步骤s5中k时刻和k-1时刻检测到的电压分别为v(k)和v(k-1)，电流分别为i(k)和i(k-1)；
所述步骤s6具体为：
s61、计算δvk＝v(k)-v(k-1)和δik＝i(k)-i(k-1)；
其中，δvk为相邻时刻检测到电压差值，δik为相邻时刻检测到的电流差值；
s62、判断|δvk|是否小于δ1，若是则进入步骤s621，否则进入步骤s622；
其中，δ1预先设定的相邻时刻检测到电压差值的阈值；
s621、判断|δik|是否小于δ2，若是则进入s63，否则进入步骤s6211；
其中，δ2预先设定的相邻时刻检测到电流差值的阈值；
s6211、判断δik是否大于0，若是则按照一定步长增加pwm占空比，并进入步骤s63；否则按照一定步长减少pwm占空比，并进入步骤s63；
s622、计算g＝i(k)/v(k)和δgk＝δik/δvk；
其中，g为外部负载的电导，δgk为相邻时刻光伏阵列的微变电导；
s623、判断∣g+δgk∣是否小于δ3，若是进入步骤s63；否则进入步骤s6231；
s6231、判断g+δgk是否大于0；若是则以一定步长增加电压，并进入步骤s63；否则以一定步长减少电压，并进入步骤s63；
s63、完成本次最大功率点的搜索。
进一步地，所述输出特性曲线进行多项式拟合前，需要对利用bezier函数的建立光伏模块模型的两组计算结果进行数据整合，并去除其中的重复点数据；
所述重复点数据为最大功率点的电压和电流数据。
本发明的有益效果为：本发明提供的基于bezier多项式的光伏阵列最大功率跟踪方法能够进行有效的最大功率跟踪，并且具有实现简单且跟踪精度好的特点。
附图说明
图1为本发明提供的实施例中基于bezier多项式的光伏阵列最大功率方法实现流程图。
图2为本发明提供的实施例中利用bezier函数进行光伏模块建模方法实现流程图。
图3本发明提供的实施例中利用增量电导法进行最大功率点微调方法实现流程图。
图4为本发明提供的实施例中bn-20m输出特性曲线图。
图5为本发明提供的实施例光伏阵列最大功率跟踪模型框图。
图6为本发明提供的实施例中光伏阵列输出功率示曲线图。
图7为本发明提供的实施例中光伏阵列输出电流曲线图。
图8为本发明提供的实施例中光伏阵列输出电压曲线图。
具体实施方式
下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。
如图1所示，基于bezier多项式的光伏阵列最大功率跟踪方法，包括以下步骤：
s1、测试光伏阵列的实际辐照度和温度值；
s2、计算实际测试辐照度和温度下，光伏模块的开路电压vocb，短路电流iscb和最大功率点(vmb，imb)；
其中，vmb为最大功率点电压，imb为最大功率点电流；
由光生伏打效应可得，入射光粒子的能量对光伏半导体材料中的电子和空穴等载流子的运动具有重要影响，而温度同样直接关系到载流子的扩散和漂移的过程，进而影响到光伏阵列的实际出力特性，因此需考虑付照度和温度对光伏阵列实际输出特性的影响效果。因此，在给定了温度和辐照度变化的非标准测试条件工况条件下时，上述步骤s2中，
最大功率点vmb的计算公式为：
vmb＝vmln(e+βδir)(1-γδt)
最大功率点电流imb的计算公式为：
短路电流iscb的计算公式为：
开路电压vocb的计算公式为：
vocb＝vocln(e+βδir)(1-γδt)
其中，δt为实际工况下的电池温度与标准测试条件下的电池温度差；
δir为实际工况下的辐照度与标准测试条件下辐照度的差；
α、β和γ均为拟合常数，且α＝0.0025/℃，β＝0.0005/(w/m2)，γ＝0.00288/℃；
irb为标准测试条件下的辐照度，且irb＝1000w/m2，
t标准测试条件下的电池温度，且t＝25℃；
vm、im是标准测试工况下光伏阵列最大功率点处的电压和电流；
isc、voc是标准测试工况下光伏阵列短路电流和开路电压；
ir是实际测试的辐照度。
利用上述计算公式可对辐照度和温度变化时光伏阵列最大功率、开路电压和短路电流等参数进行计算，进而可用bezier函数方法对不同辐照度和温度下的光伏阵列输出最大功率点左右两侧的输出曲线进行分别建模。
s3、基于开路电压vocb，短路电流iscb和最大功率点(vmb，imb)，利用bezier函数的光伏模块建模方法，得到给定条件下光伏模块的输出特性曲线；
如图2所示，上述步骤s3具体为：
s31、建立平面直角坐标系，在建立的平面直角坐标系中确定短路电流点、最大功率点以及开路电流点的位置，并根据其位置，确定bezier函数的两个控制点在平面直角坐标系中的位置；
s32、根据所述短路电流点、最大功率点、开路电流点的位置，根据bezier函数特性，确定最大功率点左侧和右侧的bezier函数曲线上点的坐标；
s33、根据最大功率点左侧和右侧的点的坐标，绘制出光伏模块的输出特性曲线。
s4、对输出特性曲线进行多项式拟合，并给出对应输出特性曲线的多项式拟合结果；
步骤s3中，基于bezier函数建模方法，可得到光伏阵列的输出特性曲线数据，为实现简单的最大功率跟踪，利用matlab软件可以很容易地对bezier函数的计算结果进行多项式拟合。但在多项式拟合之前，需预先对两组bezier函数的结果数据进行整合，并去除其中的重复点数据，其重复数据是指最大功率点的电压和电流数据，产生重复点数据的根本原因在于以bezier函数计算光伏阵列的输出特性时，分别将短路电流点，最大功率点和开路电压点分别作为二阶bezier函数的起点和终点，导致最大功率点数据重复出现两次，因而，需要在多项式拟合前对重复点数据进行处理。并且，大量统计规律表明，利用8次多项式对光伏阵列的拟合效果最好。因此，该拟合多项式结果是一个以电压为自变量，以电流为函数的8次多项式。
s5、利用电流测试装置检测光伏阵列的输出电流，同时利用电压检测装置检测光伏模块的输出电压；
s6、根据检测到的输出电压和电流，基于增量电导法，对多项式拟合结果得出的最大功率点进行逐一微调；
根据光伏阵列的确定温度和辐照度时的电压/电流和电流/功率曲线，可知光伏阵列的电压功率曲线是一个单峰的曲线，在输出功率最大点，功率对电压的导数为零，要寻找最大功率点，只要在功率对电压的导数大于零的区域增加电压，在功率对电压的导数小于零的区域减小电压，在导数等于零或非常接近于零的时候，电压保持不变即可；当电压不变电流增加时，增加工作电压，在电压不变时电流减小时，减小工作电压。
上述步骤s5中k时刻和k-1时刻检测到的电压分别为v(k)和v(k-1)，电流分别为i(k)和i(k-1)；
如图3所示，上述步骤s6具体为：
s61、计算δvk＝v(k)-v(k-1)和δik＝i(k)-i(k-1)；
其中，δvk为相邻时刻检测到电压差值，δik为相邻时刻检测到的电流差值；
s62、判断|δvk|是否小于δ1，若是则进入步骤s621，否则进入步骤s622；
其中，δ1预先设定的相邻时刻检测到电压差值的阈值；
s621、判断|δik|是否小于δ2，若是则进入s63，否则进入步骤s6211；
其中，δ2预先设定的相邻时刻检测到电流差值的阈值；
s6211、判断δik是否大于0，若是则按照一定步长增加pwm占空比，并进入步骤s63；否则按照一定步长减少pwm占空比，并进入步骤s63；
s622、计算g＝i(k)/v(k)和δgk＝δik/δvk；
其中，g为外部负载的电导，δgk为相邻时刻光伏阵列的微变电导；
s623、判断∣g+δgk∣是否小于δ3，若是进入步骤s63；否则进入步骤s6231；
s6231、判断g+δgk是否大于0；若是则按照一定步长增加电压，并进入步骤s63；否则按照一定步长减少电压，并进入步骤s63；
s63、完成最大功率点的微调。
上述增量电导法通过设定一些很小的变化阈值δ1、δ2和δ3，时光伏阵列最后稳定在最大点附近的某个点，而不是来回跳动，当光伏阵列模型系统从一个稳态过度到另一个稳态时，增量电导法能根据电流的变化能够做出正确的判断，实现对计算最大功率点的微调。
s7、将检测到的输出电压值作为自变量，使光伏模块工作输出特性在微调后的最大功率点电压处，使光伏模块工作在理论的最大输出功率，实现光伏阵列的最大功率跟踪。
在发明的一个实施例中，利用步骤s31～s33所述流程，对光伏阵列pm两侧的输出特性分别进行拟合计算，以bn20-m光伏阵列为例，选取bezier曲线控制点距最大功率点(左、右)两侧的距离分别为lp1p2长度的0.113倍和0.149倍。结果如图4所示，图4中实线为实测结果，虚线为利用bezier函数计算的结果。
表1bn-20m计算结果及误差
表1中，电流为实测结果，δi为实测数据与bezier函数拟合电流之间的绝对误差，可知，bn-20m的最大绝对误差为126.3ma，最大相对误差为13.44％，且相对误差和绝对误差均出现在最大功率点的右侧，导致其误差较大的原因主要由以下几个方面，首先，最大功率点右侧输出特性曲线比较陡峭，表现出准电压源的特性，在较小的电压变化范围内，电流变化的尺度非常大，导致这部分拟合精度低于最大功率点左侧准电流源特性部分；其次，实测数据精度较低也是导致其误差较大的原因之一。进一步计算其平均绝对误差为37.9ma，平均相对误差为3.5％，拟合精度仍然能够较好地满足工程需求，证明了利用bezier函数进行光伏输出曲线建模方法的有效性和准确性。如表2所示，
表2不同工况参数表
对不同辐照度和温度下(表2中的前两列)的光伏阵列的电流i和u进行多项式拟合，其结果如下：
i1＝1.069e-08*u^8-8.804e-07*u^7+2.833e-05*u^6-0.000464*u^5+0.00417*u^4-0.02037*u^3+0.0496*u^2-0.0514*u+2.753
i2＝9.133e-09*u^8-7.315e-07*u^7+2.289e-05*u^6-0.000365*u^5+0.00318*u^4-0.0151*u^3+0.0358*u^2-0.0362*u+1.882
i3＝6.808e-09*u^8-5.258e-07*u^7+1.587e-05*u^6-0.000244*u^5+0.00205*u^4-0.0094*u^3+0.0215*u^2-0.0209*u+1.049
表2中(后3列)为不同工况下在给定的光伏阵列辐照度和温度下根据公式(1)～(4)得到的最大功率点的计算值，根据表1中的计算结果可对光伏阵列最终最大功率跟踪结果进行对比和验证。
在本发明的一个实施例中，提供了利用增量电导法建立的基于bezier多项式的最大功率跟踪模型如图5所示，基于bezier的光伏阵列最大功率跟踪方法的实现是在升压斩波电路(boost)电路中进行的，因为boost电路可简单地实现增压输出的功能，对输出电压较低的光伏阵列比较适用。分析图4中boost电路的工作原理，可知在不考虑功率损失的情况下，该模型系统的输入功率和输出功率相等，即：
v1i1＝v2i2
式中，v1，i1分别为boost输入(光伏输出)的电压和电流；v2，i2分别为boost输出电压和电流；
并且，输入电压、电流和输出电压、电流之间满足：
式中，d为boost电路的占空比。
光伏阵列和负载侧的电阻r1与r2满足
综上可知，采用的boost电路可以通过调节占空比d，达到调节电池侧等效电阻r1的目的，进而调节boost电路的工作点，以使光伏阵列运行在最大功率点附近。
上述最大功率跟踪模型利用受控电流源等效实际光伏阵列的输出，通过检测受控电流源两端的电压作为bezier函数多项式y1(u)，y2(u)，y3(u)的输入信号，y1(u)，y2(u)，y3(u)分别为按照表2中给定系数的光伏阵列输出功率多项式的情景a、情景b和情景c一致。
上述最大功率跟踪模型中，step模块用来确定不同辐照度和温度下光伏阵列输出曲线拟合多项式的接入时间，step以三个阶跃函数确定光伏输出曲线的接入时间，具体为在0s时，光伏阵列工作在情景a的工况下，0.5s时工作在情景b的工况下，而在1s时工作在情景c的工况下。对该模型进行仿真条件为ode23tb算法，仿真时，用增量电导对最大功率点进行微调时，最小步长为1e-7；最大步长为1e-6；相对误差为1e-5，c1＝c2＝300uf，l1＝80uh，r＝200ω。相应的仿真结果分别如图6～8所示。
将图6～图8中的仿真结果与表1中的bezier理论计算结果相比可知，在时间分别为0s，0.5s，1s时，光伏阵列的输出功率均可靠地运行在实际最大功率点188.7w，118.3w和61.1w处，所对应的输出电流和电压的一致性也较好。表明所构建的bezier函数多项式的最大功率跟踪模型能够进行有效的最大功率跟踪，并且具有实现简单且跟踪的精度较好的特点。
本发明的有益效果为：本发明提供的基于bezier多项式的光伏阵列最大功率跟踪方法能够进行有效的最大功率跟踪，并且具有实现简单且跟踪精度好的特点。