CN103293950B - 一种基于lssvm的光伏最大功率点跟踪的控制方法 - Google Patents
一种基于lssvm的光伏最大功率点跟踪的控制方法 Download PDFInfo
- Publication number
- CN103293950B CN103293950B CN201310120738.7A CN201310120738A CN103293950B CN 103293950 B CN103293950 B CN 103293950B CN 201310120738 A CN201310120738 A CN 201310120738A CN 103293950 B CN103293950 B CN 103293950B
- Authority
- CN
- China
- Prior art keywords
- training
- model
- power point
- vector machine
- maximum power
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000012549 training Methods 0.000 claims abstract description 48
- 238000012360 testing method Methods 0.000 claims abstract description 23
- 230000005855 radiation Effects 0.000 claims abstract description 13
- 238000002790 cross-validation Methods 0.000 claims abstract description 11
- 238000004088 simulation Methods 0.000 claims description 18
- 238000012843 least square support vector machine Methods 0.000 claims description 10
- 238000012706 support-vector machine Methods 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 4
- 238000012795 verification Methods 0.000 claims description 4
- 241001497337 Euscorpius gamma Species 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 claims description 3
- 238000013178 mathematical model Methods 0.000 claims description 3
- 230000009514 concussion Effects 0.000 abstract 1
- 230000006870 function Effects 0.000 description 8
- 230000007547 defect Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000005457 optimization Methods 0.000 description 5
- 238000010248 power generation Methods 0.000 description 5
- 238000013528 artificial neural network Methods 0.000 description 3
- 230000010355 oscillation Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/50—Photovoltaic [PV] energy
- Y02E10/56—Power conversion systems, e.g. maximum power point trackers
Landscapes
- Photovoltaic Devices (AREA)
Abstract
本发明公开了基于LSSVM的光伏最大功率点跟踪的控制方法,其步骤有:(1)建立光伏电池模型;(2)确定模型中设定温度值和太阳辐射强度下MPPT对应的电压和电流;(3)根据温度、太阳辐射强度、电压及电流的序列,得到最小二乘支持向量机模型中的训练集和测试集样本;(4)利用交叉验证方法选择最小二乘支持向量机的训练集的训练参数;(5)利用训练集的样本数据,通过选择的训练参数对最小二乘支持向量机进行训练,得到训练模型;(6)根据训练模型对所述测试集进行测试并对结果进行评价。本发明通过挖掘已知数据的内在规律,能够有效实现最大功率点处的输出电压和电流的跟踪,并且有效避免系统在最大功率点处震荡运行的缺点。
Description
技术领域
本发明属于光伏发电领域,具体涉及一种基于LSSVM的光伏最大功率点跟踪的控制方法。
背景技术
随着社会的快速发展,能源危机和环境污染成为亟待解决的问题。可再生能源的利用成为近年来各个领域研究的热点。太阳能具有普遍的存在性、储量的无限性、使用的清洁性等优点,是人类社会很理想的未来能源。利用太阳能进行发电是可持续发展能源战略的重要决策。
为了实现光能的最大利用率,需要对光伏电池进行最大功率点跟踪。目前,常用的实现最大功率点跟踪的方法主要有以下几种:恒定电压法、扰动观察法、电导增量法、最优梯度法、神经网络法等。其中,恒定电压法在温度变化不大时具有良好的可靠性和稳定性,但实际中,一年四季或每天晨午温差较大,该方法会降低系统运行效率。扰动观察法结构简单,测量参数少。电导增量法在外界温度和光照等环境发生变化时,其动态特性和跟踪特性较好。但扰动观察法和电导增量法由于步长固定,都存在一个缺点,即步长过小,则导致光伏电池长时间停留在低功率工作状态;如果步长过大,则导致光伏发电系统输出振荡加剧。针对此缺点,有学者提出变步长方法,当离最大功率点较远时,采用较大步长,加快寻优速度;当离最大功率点较近时,采用较小步长,保证跟踪精度,利于系统稳定。最优梯度法是一种以梯度法为基础的多维无约束最优化问题的数值计算法,该方法中需要选择电压增量系数,相当于扰动观察法和电导增量法中的步长。神经网络具有大规模并行处理能力、容错性、自组织及自适应能力强等特点,但也有易陷入局部最小、神经网络结点个数难以确定等不可避免的缺点。
发明内容
针对现有技术的不足,本发明提出一种基于MATLAB/Simulink中LSSVM工具箱的光伏最大功率点跟踪的控制方法,能够有效避免系统在最大功率点处震荡运行。
本发明提供的一种基于LSSVM的光伏最大功率点跟踪的控制方法,其改进之处在于,所述方法包括如下步骤:
(1)建立光伏电池模型;
(2)确定模型中设定温度值和太阳辐射强度下MPPT对应的电压和电流;
(3)根据温度、太阳辐射强度、电压及电流序列,得到最小二乘支持向量机模型中的训练集和测试集样本;
(4)利用交叉验证方法选择最小二乘支持向量机的训练集的训练参数;
(5)利用训练集的样本数据,通过选择的所述训练参数对最小二乘支持向量机进行训练,并根据最优参数得到训练模型;
(6)根据所述训练模型对所述测试集进行测试并对结果进行评价。
其中,步骤(1)所述建立光伏电池模型是指根据光伏电池的工程数学模型,在MATLAB中建立其仿真模型。
其中,步骤(3)训练集和测试集样本的数据集为:
其中,
写成矩阵的形式即为:
式中,x1(n)为温度序列,x2(n)为太阳辐射强度序列;y1(n)为在最大功率点的输出电压序列,y2(n)为在最大功率点的输出电流序列。
其中,步骤(4)所述训练参数包括正则化参数γ和核函数参数σ2;其核函数的表达式如下:
其中,步骤(4)所述利用交叉验证方法选择最小二乘支持向量机的训练集的训练参数包括:
将训练集分割成K个子集,其中一个子集被保留作为验证模型的数据,其他K-1个子集用来训练,交叉验证重复K次,每个子集验证一次;
在训练前,给定正则化参数γ和核函数参数σ2的范围,即γ∈[γmin,γmax],σ2∈[σ2 min,σ2 max],使γ和σ2在这个范围内遍历取值,对每一个组合参数(γ,σ2)均进行10次交叉验证;
选取均方误差MSE(MeanSquareError)作为评价指标,将这10次测试结果的MSE取平均值作为该参数对应的指标,最后选择指标最高的一组参数组合作为最终的最优参数。
与现有技术比,本发明的有益效果为:
(1)本发明通过挖掘已知数据的内在规律,能够有效实现最大功率点处的输出电压和电流的跟踪。
(2)本发明能够有效避免系统在最大功率点处震荡运行的缺点。
(3)本发明能够避免算法中有步长时的缺点,即步长过小,则导致光伏电池长时间停留在低功率工作状态;如果步长过大,则导致光伏发电系统输出振荡加剧。
(4)本发明是一种基于结构风险最小化的学习算法,既有严格的理论基础,又能较好地解决过学习、欠学习、非线性、高维数和局部最优等问题。
(5)本发明通过选择支持向量机进行计算,是一种新的思路,为以后的理论知识提供基础。
附图说明
图1为本发明提供的光伏电池I-V特性曲线示意图。
图2为本发明提供的光伏电池P-V特性曲线示意图。
图3为本发明提供的利用LSSVM实现MPPT的流程图。
图4为本发明提供的测试集仿真结果示意图。
图5为本发明提供的测试集仿真结果误差示意图。
图6为本发明提供的测试集仿真结果相对误差示意图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步的详细说明。
本实施例提供的一种基于LSSVM的光伏最大功率点跟踪的控制方法,其流程图如图3所示,具体包括如下步骤:
步骤(1),在MATLAB仿真软件中建立光伏电池模型。根据光伏电池的工程数学模型,在MATLAB中建立其仿真模型,本实施例中,设置光伏电池的参数如下:Voc=42V,Isc=4.5A,Vm=34V,Im=4A,Tref=25°C,Sref=1kW/m2。该模型输出特性曲线的仿真结果如图1和图2所示。
步骤(2),根据步骤(1)中所述的光伏电池模型,得到某些温度和太阳辐射强度下MPPT对应的电压和电流。分别将温度在0~40°C之间、太阳辐射强度在0~4kW/m2之间间隔取值,经过仿真分析得到对应温度和太阳辐射强度条件下的最大功率点输出电压和电流。
步骤(3),根据步骤(2)中所述的电压和电流序列,得到最小二乘支持向量机模型中的训练集和测试集样本。
得到的数据集为
其中,
写成矩阵的形式即为:
式中,x1(n)为温度序列,x2(n)为太阳辐射强度序列;y1(n)为在最大功率点的输出电压序列,y2(n)为在最大功率点的输出电流序列。
具体样本数据如表1和表2所示。
表1训练集数据
表2测试集数据
根据已知某些典型温度和太阳辐射强度条件下的最大功率点以及工作在最大功率点处的输出电压和输出电流,可以找到一个映射函数,使得
Y(n)=F(X(n))(4)
最小二乘支持向量是在标准支持向量机的基础上发展而来的,它将标准支持向量机中约束条件中的不等式约束变为等式约束,从而使求解二次规划问题转化为求解一组线性方程问题回归模型的优化问题。描述如下:
设决策函数为y=(ω·x)+b,最优化问题变为:
引入Lagrange函数
根据KKT条件,最优化条件变为:
写成矩阵形式并消去ω和ξ,化简得到
求解此线性方程组就可以求得决策函数。对于非线性回归问题,引入核函数方法来求解即可。选用RBF核函数,如式(9)所示。故最小二乘支持向量机预测模型中有两个参数需要选择,即正则化参数γ和核函数参数σ2。
步骤(4),利用交叉验证方法对最小二乘支持向量机的训练参数进行选择。
K折交叉验证方法的思路如下:将训练集分割成K个子集,K一般大于等于2,其中一个子集被保留作为验证模型的数据,其他K-1个子集用来训练。交叉验证重复K次,每个子集验证一次。在训练前,给定正则化参数γ和核函数参数σ2的范围,即γ∈[γmin,γmax],σ2∈[σ2 min,σ2 max],使γ和σ2在这个范围内遍历取值,对每一个组合参数(γ,σ2)均进行10次交叉验证(K=10)。选取均方误差MSE(MeanSquareError)作为评价指标,将这10次测试结果的MSE取平均值作为该参数对应的指标,最后选择指标最高(即MSE最小)的一组参数组合作为最终的最优参数。
步骤(5),利用训练集样本数据,采用步骤(4)中所述方法选择到的参数对最小二乘支持向量机进行训练,得到训练模型。
步骤(6),利用得到的训练模型,对测试集进行测试并对结果进行评价。测试集仿真结果如图4所示,图中,方形符号代表在MPP处的电压、电流值,星形符号代表该方法控制下输出的电压、电流值;测试集仿真结果误差如图5所示,图中,方形符号代表仿真结果中电压的误差,星形符号代表仿真结果中电流的误差;测试集仿真结果相对误差如图6所示,图中,方形符号代表仿真结果中电压的相对误差,星形符号代表仿真结果中电流的相对误差。
根据图5、图6的仿真结果,MPP处的输出电压仿真结果的相对误差最大值为3.2%,绝大部分结果相对误差小于1%。MPP处的输出电流仿真结果的相对误差最大值为1%,绝大部分结果相对误差小于0.5%。可见,利用LSSVM实现MPPT的方法,能够较准确地实现最大功率点处的输出电压和电流的跟踪。本发明能够有效避免系统在最大功率点处震荡运行的缺点。
利用LSSVM方法仿真得到的电压或者电流值,可以应用到光伏发电系统后续的控制电路中。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细的说明,所属领域的普通技术人员应当理解:依然可以对本发明的具体实施方式进行修改或者等同替换,而未脱离本发明精神和范围的任何修改或者等同替换,其均应涵盖在本发明的权利要求范围当中。
Claims (4)
1.一种基于LSSVM的光伏最大功率点跟踪的控制方法,其特征在于,所述方法包括如下步骤:
(1)建立光伏电池模型;
(2)确定模型中设定温度值和太阳辐射强度下MPPT对应的电压和电流;
(3)根据温度、太阳辐射强度、电压及电流的序列,得到最小二乘支持向量机模型中的训练集和测试集样本;
(4)利用交叉验证方法选择最小二乘支持向量机的训练集的训练参数;
(5)利用训练集的样本数据,通过选择的所述训练参数对最小二乘支持向量机进行训练,得到训练模型;
(6)根据测试集对所述训练模型进行测试并对结果进行评价;
步骤(4)所述训练参数包括正则化参数γ和核函数参数σ2;其核函数的表达式如下:
步骤(4)所述利用交叉验证方法选择最小二乘支持向量机的训练集的训练参数包括:
将训练集分割成K个子集,其中一个子集被保留作为验证模型的数据,其他K-1个子集用来训练,交叉验证重复K次,每个子集验证一次;
在训练前,给定正则化参数γ和核函数参数σ2的范围,即γ∈[γmin,γmax],σ2∈[σ2 min,σ2 max],使γ和σ2在这个范围内遍历取值,对每一个组合参数(γ,σ2)均进行10次交叉验证;
选取均方误差MSE作为评价指标,将这10次测试结果的MSE取平均值作为该组合参数对应的指标,最后选择指标最高的一组参数组合作为最终的最优参数。
2.如权利要求1所述的控制方法,其特征在于,步骤(1)所述建立光伏电池模型是指根据光伏电池的工程数学模型,在MATLAB中建立其仿真模型。
3.如权利要求1所述的控制方法,其特征在于,步骤(3)训练集和测试集样本的数据集为:
其中,
写成矩阵的形式即为:
式中,x1(n)为温度序列,x2(n)为太阳辐射强度序列;y1(n)为在最大功率点的输出电压序列,y2(n)为在最大功率点的输出电流序列。
4.如权利要求1所述的控制方法,其特征在于,步骤(5)是根据最优参数建立训练模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310120738.7A CN103293950B (zh) | 2013-04-09 | 2013-04-09 | 一种基于lssvm的光伏最大功率点跟踪的控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310120738.7A CN103293950B (zh) | 2013-04-09 | 2013-04-09 | 一种基于lssvm的光伏最大功率点跟踪的控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103293950A CN103293950A (zh) | 2013-09-11 |
CN103293950B true CN103293950B (zh) | 2016-04-20 |
Family
ID=49095024
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310120738.7A Active CN103293950B (zh) | 2013-04-09 | 2013-04-09 | 一种基于lssvm的光伏最大功率点跟踪的控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103293950B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103823504B (zh) * | 2014-03-20 | 2016-01-20 | 常州工学院 | 一种基于最小二乘支持向量机的最大功率跟踪控制方法 |
CN106130006A (zh) * | 2016-06-07 | 2016-11-16 | 塔里木大学 | 一种基于matlab微电网光伏系统稳定性分析方法 |
CN106452355B (zh) * | 2016-10-17 | 2018-04-10 | 温州大学 | 一种基于模型辨识的光伏发电系统最大功率跟踪方法 |
CN108267951A (zh) * | 2016-12-30 | 2018-07-10 | 南京理工大学 | 一种基于核偏最小二乘的最大功率点跟踪控制方法 |
CN112446391A (zh) * | 2019-08-28 | 2021-03-05 | 辽宁省交通高等专科学校 | 基于最小二乘支持向量机的太阳高度角预测方法 |
CN110985289B (zh) * | 2019-12-04 | 2021-01-01 | 浙江大学 | 基于svr和smc的风电机组带预设性能的mppt方法 |
CN114510111B (zh) * | 2021-12-29 | 2023-09-12 | 北京华能新锐控制技术有限公司 | 部分遮阳光伏阵列的全局mppt控制方法与装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101877559A (zh) * | 2010-06-29 | 2010-11-03 | 苏州市职业大学 | 光伏发电系统最大功率点快速准确跟踪控制的复合控制方法 |
CN102817781A (zh) * | 2012-09-06 | 2012-12-12 | 华锐风电科技(集团)股份有限公司 | 风力发电机的风能捕获控制方法和装置及风力发电系统 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3318974B2 (ja) * | 1992-07-16 | 2002-08-26 | 日本電池株式会社 | 太陽電池の最大電力点追跡制御法 |
JPH08123561A (ja) * | 1994-10-20 | 1996-05-17 | Meidensha Corp | 太陽光発電システムの最大出力追従制御方法および装置 |
-
2013
- 2013-04-09 CN CN201310120738.7A patent/CN103293950B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101877559A (zh) * | 2010-06-29 | 2010-11-03 | 苏州市职业大学 | 光伏发电系统最大功率点快速准确跟踪控制的复合控制方法 |
CN102817781A (zh) * | 2012-09-06 | 2012-12-12 | 华锐风电科技(集团)股份有限公司 | 风力发电机的风能捕获控制方法和装置及风力发电系统 |
Non-Patent Citations (2)
Title |
---|
分布式发电系统储能优化配置;田军;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20110815(第8期);论文第24-31页 * |
基于MATLAB的光伏电池通用数学模型;王长江等;《电力科学与过程》;20090430;第25卷(第4期);第12-13页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103293950A (zh) | 2013-09-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103293950B (zh) | 一种基于lssvm的光伏最大功率点跟踪的控制方法 | |
Ge et al. | A hybrid model for short-term PV output forecasting based on PCA-GWO-GRNN | |
Ciulla et al. | Forecasting the cell temperature of PV modules with an adaptive system | |
CN103823504B (zh) | 一种基于最小二乘支持向量机的最大功率跟踪控制方法 | |
CN107562992B (zh) | 一种基于svm和粒子群算法的光伏列阵最大功率跟踪方法 | |
CN103258103A (zh) | 基于偏最小二乘回归的戴维南等值参数辨识方法 | |
CN106126887A (zh) | 一种基于参数辨识的多峰值光伏阵列最大功率点跟踪方法 | |
CN105512775A (zh) | 一种光伏发电系统的功率预测方法 | |
Bi et al. | An Enhanced 0.8$\text {V} _\text {OC} $-Model-Based Global Maximum Power Point Tracking Method for Photovoltaic Systems | |
CN103973221A (zh) | 一种基于实测数据的光伏阵列参数辨识方法 | |
Huang et al. | Parameter identification for photovoltaic models using an improved learning search algorithm | |
Liu et al. | A fault detection and diagnosis technique for solar system based on Elman neural network | |
Mellit et al. | A survey on the application of artificial intelligence techniques for photovoltaic systems | |
CN105913161A (zh) | 一种基于多目标优化的光伏系统最大功率点的获取方法 | |
CN103995559B (zh) | 一种基于环境参数模型的定电压mppt控制方法和系统 | |
Mao et al. | A convolutional neural network-based maximum power point voltage forecasting method for pavement PV array | |
Wang et al. | Digital twin based maximum power point estimation for photovoltaic systems | |
Hadji et al. | Experimental analysis of genetic algorithms based MPPT for PV systems | |
He | Solar Maximum Power Point Tracking Algorithm Based on Wavelet Neural Network | |
CN116258057A (zh) | 一种基于Transformer进行光伏发电预测分析的方法 | |
El Marghichi | A Solar PV Model Parameter Estimation Based on the Enhanced Self-Organization Maps | |
Liu et al. | FSR‐SSL: A fault sample rebalancing framework based on semi‐supervised learning for PV fault diagnosis | |
Bennagi et al. | Comprehensive study of the artificial intelligence applied in renewable energy | |
Yixin et al. | Improved SPSO-based parameter identification of solar PV cells IV model | |
CN108776843A (zh) | 基于波尔兹曼机的光伏超短期功率预测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |