CN109814651B - 基于粒子群的光伏电池多峰值最大功率跟踪方法及系统 - Google Patents

Abstract

本发明公开了一种基于粒子群的光伏电池多峰值最大功率跟踪方法及系统，本发明通过获得局部阴影条件下光伏电池输出的电压以及电流数据值，将粒子初始化在可能的电压值处，通过迭代的方式不断更新粒子最大功率点，寻找个体历史最优位置和全局最优位置，从而确定复杂局部阴影条件下的多峰值最大功率点，达到相对应的结束条件之后搜寻最大功率点结束。与传统的最大功率跟踪方法相比，本发明公开的最大功率跟踪方法可以避免搜寻陷入局部最优，保证跟踪的精度和速度，减小了跟踪过程中的震荡损耗，更符合实际。

Description

基于粒子群的光伏电池多峰值最大功率跟踪方法及系统

技术领域

本发明涉及一种基于粒子群的光伏电池多峰值最大功率跟踪方法及系统，属于光伏发电工程领域。

背景技术

随着矿石能源的逐步消耗以及环境污染的逐步加剧，可再生能源越来越多的受到人们的青睐。作为一种环保、可再生的能源，太阳能光伏发电技术已经被广泛应用，光伏(PV)装置被越来越多地投入到各种应用中。但是光电转换效率问题一直是阻碍光伏发电发展的一个技术瓶颈。为了更大限度地发展和利用太阳能资源，对光伏发电系统的最大功率点进行跟踪是提高光电转换效率十分有效的途径，然而最大功率跟踪的主要挑战是要应对它随温度和日照变化的非线性输出特性，并且如果整个阵列没有得到均匀光照，例如在分布式光伏发电的城市中，遇到树木的遮挡、云层的变化以及灰尘的不均匀积累，该特性变得更为复杂，出现多峰特性。多峰的存在降低了现有最大功率点跟踪(MPPT)方案的有效性，因为它们不能辨别局部最高点和整体最高点。

因此对部分遮蔽情况下光伏系统最大功率点跟踪的研究至关重要。现有传统技术多为单峰值最大功率跟踪方法，而传统方法跟踪速度慢，且在稳态时容易发生功率震荡，增加功率损耗，更为主要的是，极易发生误判陷入局部最优。

基于以上分析，本发明提出了一种基于粒子群算法的最大功率跟踪方法，可以有效避免陷入局部最优获得局部阴影条件下光伏电池的组的最大功率点。

发明内容

本发明要解决的技术问题在于，针对现有技术多为单峰值最大功率跟踪方法，而传统方法跟踪速度慢，且在稳态时容易发生功率震荡，增加功率损耗，极易发生误判陷入局部最优的技术缺陷，提供了一种基于粒子群的光伏电池多峰值最大功率跟踪方法。

本发明解决其技术问题，所采用的基于粒子群的光伏电池多峰值最大功率跟踪方法包括：

S1、根据光伏电池阵列的开路电压U_oc，在[0,U_oc]上选取m个点，作为粒子i的初始电压V_i(k)，i∈[1,2,3,…,m]，k表示迭代的次数，初始电压时k＝0,m为大于1的正整数；初始化粒子种群的参数，包括了粒子权重ω的取值，自我学习因子c₁与社会学习因子c₂的设定范围；

S2、计算各个粒子的适应度值，适应度值等于光伏电池阵列的输出功率；根据各个粒子的适应度值寻找种群的个体最优位置与全局最优位置；

S3、更新各个粒子的位置与速度，得到新的下一代粒子种群的适应度值；

S4、分别比较下一代粒子种群的个体最优位置与全局最优位置与本次更新之前的个体最优位置与全局最优位置，以分别更新个体最优位置与全局最优位置，更新的准则：计算当前位置的适应度值后与本次更新之前的最优适应度值比较，将较大的适应度值对应的位置作为该粒子的新的个体最优位置，将所有的新的个体最优位置粒子中最大的适应度值对应的位置作为全局最优位置；

S5、返回步骤S3，直至达到收敛条件或者达到最大迭代次数，将最后的全局最优位置作为光伏系统的最大功率点所对应的工作电压。

进一步地，在本发明的基于粒子群的光伏电池多峰值最大功率跟踪方法中，所述的在[0,U_oc]上选取m个点的选取的要求是线性均匀选取或者是将粒子均匀分布在理论值处，所述的理论值数值等于0.8*U_oc。

进一步地，在本发明的基于粒子群的光伏电池多峰值最大功率跟踪方法中，所述的学习因子c₁及c₂的设定范围是(0,2)，学习因子c₁及c₂在(0,2)区间内随机取值。

进一步地，在本发明的基于粒子群的光伏电池多峰值最大功率跟踪方法中，所述的粒子权重ω取值是线性递减权重，计算公式为：

其中c为介于0至1之间的调节系数，k为当前迭代的次数，maxgen为最大迭代次数，ω_min为最小惯性权重，ω_max为最大惯性权重，f_i ^k为第k代中的第i个粒子的适应度值，

和

分别为第k代中的第i个粒子的最小和最大适应度值，

为第k代中的第i个粒子的权重。

进一步地，在本发明的基于粒子群的光伏电池多峰值最大功率跟踪方法中，所述更新各个粒子的位置与速度度具体是指：

在n维的搜索空间中，由m个粒子组成的种群记为X＝(x₁,...,x_i，...,x_m)，第i个粒子所在的位置为x_i＝(x_i1，x_i2，...x_in)^T，速度为v_i＝(v_i1，v_i2，...v_in)^T，第i个粒子在空间搜索过程中搜索到的自身最优位置P_i＝(p_i1，p_i2，...p_id)^T，全局最优位置P_g＝(p_g1，p_g2，...p_gd)^T，每个粒子代表一个解，粒子的位置与速度以迭代的方式进行更新，公式为：

是粒子i在第k次迭代的速度，是

第d维的速度分量；

是粒子i在第k次的位置，

是

第d维的位置分量；

是粒子i在第k次迭代的个体最优位置，

是

第d维的个体最优位置分量；

是粒子群在第k次迭代的全局最优位置，

是

第d维的最优位置分量，

为第k代中的第i个粒子的权重，r₁与r₂是在[0,1]区间内均匀分布的随机数。

进一步地，在本发明的基于粒子群的光伏电池多峰值最大功率跟踪方法中，所述达到收敛条件是指：

判断所有粒子所对应的适应度值最大与最小的差值，如果差值小于预设差值，表示达到收敛条件，此时粒子群算法执行完毕，否则未达到收敛条件。

进一步地，在本发明的基于粒子群的光伏电池多峰值最大功率跟踪方法中，适应度值的计算公式为：

是粒子i在第k次迭代的速度，代表光伏电池输出的电压值；

代表光伏电池当前输出电压测得的电流；

是粒子的目标函数值代表当前光伏电池输出电压对应的输出功率，即适应度值。

本发明为解决其技术问题，所提供的一种基于粒子群的光伏电池多峰值最大功率跟踪系统包括：

初始化模块，用于根据光伏电池阵列的开路电压U_oc，在[0,U_oc]上选取m个点，作为粒子i的初始电压V_i(k)，i∈[1,2,3,…,m]，k表示迭代的次数，初始电压时k＝0,m为大于1的正整数；初始化粒子种群的参数，包括了粒子权重ω的取值，自我学习因子c₁与社会学习因子c₂的设定范围；

最优位置计算模块，用于计算各个粒子的适应度值，适应度值等于光伏电池阵列的输出功率；根据各个粒子的适应度值寻找种群的个体最优位置与全局最优位置；

粒子种群更新模块，用于更新各个粒子的位置与速度，得到新的下一代粒子种群的适应度值；

最优位置更新模块，用于分别比较下一代粒子种群的个体最优位置与全局最优位置与本次更新之前的个体最优位置与全局最优位置，以分别更新个体最优位置与全局最优位置，更新的准则：计算当前位置的适应度值后与本次更新之前的最优适应度值比较，将较大的适应度值对应的位置作为该粒子的新的个体最优位置，将所有的新的个体最优位置粒子中最大的适应度值对应的位置作为全局最优位置；

最优位置更新模块，用于返回粒子种群更新模块，直至达到收敛条件或者达到最大迭代次数，将最后的全局最优位置作为光伏系统的最大功率点所对应的工作电压。

进一步地，在本发明的基于粒子群的光伏电池多峰值最大功率跟踪系统中，所述的粒子权重ω取值是线性递减权重，计算公式为：

其中c为介于0至1之间的调节系数，k为当前迭代的次数，maxgen为最大迭代次数，ω_min为最小惯性权重，ω_max为最大惯性权重，f_i ^k为第k代中的第i个粒子的适应度值，

和

分别为第k代中的第i个粒子的最小和最大适应度值，

为第k代中的第i个粒子的权重。

进一步地，在本发明的基于粒子群的光伏电池多峰值最大功率跟踪系统中，所述更新各个粒子的位置与速度度具体是指：

在n维的搜索空间中，由m个粒子组成的种群记为X＝(x₁,...,x_i，...,x_m)，第i个粒子所在的位置为x_i＝(x_i1，x_i2，...x_in)^T，速度为v_i＝(v_i1，v_i2，...v_in)^T，第i个粒子在空间搜索过程中搜索到的自身最优位置P_i＝(p_i1，p_i2，...p_id)^T，全局最优位置P_g＝(p_g1，p_g2，...p_gd)^T，每个粒子代表一个解，粒子的位置与速度以迭代的方式进行更新，公式为：

是粒子i在第k次迭代的速度，是

第d维的速度分量；

是粒子i在第k次的位置，

是

第d维的位置分量；

是粒子i在第k次迭代的个体最优位置，

是

第d维的个体最优位置分量；

是粒子群在第k次迭代的全局最优位置，

是

第d维的最优位置分量，

为第k代中的第i个粒子的权重，r₁与r₂是在[0,1]区间内均匀分布的随机数。

实施本发明的基于粒子群的光伏电池多峰值最大功率跟踪方法及系统，具有以下有益效果：

1.本发明利用基于微生物行为机制的粒子群算法的全局搜索特点应用到光伏电池的最大功率跟踪方面，避免了陷入局部最优，提升了搜寻速度，减小了震荡，降低了损耗。

2.在粒子群算法的参数设置中，将传统的权重保持不变，优化为一个深化到粒子群的每一代和代中的每一个粒子的自适应调整的惯性权重参数提高了精度，从而更好的实现最大功率跟踪。

附图说明

下面将结合附图及实施例对本发明作进一步说明，附图中：

图1为本发明实施例的光伏最大功率跟踪的电路结构示意框图；

图2为本发明实施例的一种光伏电池多峰值粒子群最大功率跟踪控制方法的流程图；

图3为本发明实施例的一种光伏电池多峰值粒子群最大功率跟踪控制方法的电路结构图；

图4为本发明实施例的一种光伏电池多峰值粒子群最大功率跟踪控制方法的仿真模型；

图5为采用本发明方法控制后的光伏阵列输出功率曲线。

具体实施方式

为了对本发明的技术特征、目的和效果有更加清楚的理解，现对照附图详细说明本发明的具体实施方式。

参考图1，其为图1为本发明实施例的光伏最大功率跟踪的电路结构示意框图，该电路工作结构包括：

光伏电池阵列模块、DC/DC变换器模块、负载或逆变模块，其中控制单元包括数据采集单元、粒子群算法以及控制器单元；

数据采集单元采集光伏电池阵列得到粒子群算法所需要的输出电流与输出电压值，利用粒子群算法进行最大功率的计算与确定；

控制器单元将粒子群得到的占空比输出给DC/DC变换器，稳定工作在粒子群算法得到的工作电压点，供给负载电能。

参考图2，本发明的基于粒子群的光伏电池多峰值最大功率跟踪方法包括：

S1、根据光伏电池阵列的开路电压U_oc，在[0,U_oc]上选取m个点，作为粒子i的初始电压V_i(k)，i∈[1,2,3,…,m]，k表示迭代的次数，初始电压时k＝0,m为大于1的正整数；初始化粒子种群的参数，包括了粒子权重ω的取值，自我学习因子c₁与社会学习因子c₂的设定范围；在[0,U_oc]上选取m个点的选取的要求是线性均匀选取或者是将粒子均匀分布在理论值处，所述的理论值数值等于0.8*U_oc。学习因子c₁及c₂的设定范围是(0,2)，学习因子c₁及c₂在(0,2)区间内随机取值。

S2、计算各个粒子的适应度值，适应度值等于光伏电池阵列的输出功率；根据各个粒子的适应度值寻找种群的个体最优位置与全局最优位置。

适应度值的计算公式为：

是粒子i在第k次迭代的速度，代表光伏电池输出的电压值；

代表光伏电池当前输出电压测得的电流；

是粒子的目标函数值代表当前光伏电池输出电压对应的输出功率，即适应度值。

S3、更新各个粒子的位置与速度，得到新的下一代粒子种群的适应度值。

为了加快算法的收敛速度，在初始阶段，设置较大的ω值，可以使算法不易陷入局部最优；粒子权重ω取值是线性递减权重，计算公式为：

其中c为介于0至1之间的调节系数，k为当前迭代的次数，maxgen为最大迭代次数，ω_min为最小惯性权重，ω_max为最大惯性权重，fi_i ^k为第k代中的第i个粒子的适应度值，

和

分别为第k代中的第i个粒子的最小和最大适应度值，

为第k代中的第i个粒子的权重。

所述更新各个粒子的位置与速度度具体是指：

在n维的搜索空间中，由m个粒子组成的种群记为X＝(x₁,...,x_i，...,x_m)，第i个粒子所在的位置为x_i＝(x_i1，x_i2，...x_in)^T，速度为v_i＝(v_i1，v_i2，...v_in)^T，第i个粒子在空间搜索过程中搜索到的自身最优位置P_i＝(p_i1，p_i2，...p_id)^T，全局最优位置P_g＝(p_g1，p_g2，...p_gd)^T，每个粒子代表一个解，粒子的位置与速度以迭代的方式进行更新，公式为：

是粒子i在第k次迭代的速度，是

第d维的速度分量；

是粒子i在第k次的位置，

是

第d维的位置分量；

是粒子i在第k次迭代的个体最优位置，

是

第d维的个体最优位置分量；

是粒子群在第k次迭代的全局最优位置，

是

第d维的最优位置分量，

为第k代中的第i个粒子的权重，r₁与r₂是在[0,1]区间内均匀分布的随机数。

S4、分别比较下一代粒子种群的个体最优位置与全局最优位置与本次更新之前的个体最优位置与全局最优位置，以分别更新个体最优位置与全局最优位置，更新的准则：计算当前位置的适应度值后与本次更新之前的最优适应度值比较，将较大的适应度值对应的位置作为该粒子的新的个体最优位置，将所有的新的个体最优位置粒子中最大的适应度值对应的位置作为全局最优位置。即：计算每个粒子当前位置的适应度值，与其经历的最优适应度值比较，哪个好就把哪个粒子位置更新为p_besti(第i个粒子的最优)；在所有的粒子当中找适应度值最好的，更新g_best(全局最优)。

粒子的飞行速度由动量部分、认知部分和社会部分组成；粒子位置由上一次迭代的位置加上运动的速度；目标函数f是衡量粒子所在位置优劣的适应度函数。粒子的飞行速度v即是ΔU，v与最大功率点距离成正比，r1r2使得电压步长具有随机性，更大的概率趋近最大功率点；粒子的位置x对应直流侧Vdc；适应度函数即公式(2)判断光伏电池组件是否运行在最大功率点。本发明当中n为1即一维搜索空间，种群粒子数量m设置为10，最大迭代步数maxgen为30，输出的功率作为适应度值，加入了自适应权重，并设置了终止条件，即当最大粒子的功率值与最小粒子的功率值差值小于0.04W的时候，如果没有迭代到最大步数，则MPPT搜索结束。

S5、返回步骤S3，直至达到收敛条件或者达到最大迭代次数，将最后的全局最优位置作为光伏系统的最大功率点所对应的工作电压。达到收敛条件是指：判断所有粒子所对应的适应度值最大与最小的差值，如果差值小于预设差值，表示达到收敛条件，此时粒子群算法执行完毕，否则未达到收敛条件。

图3最大功率跟踪控制方法的电路中最大功率控制器(MPPT控制器)将测量得到的电流电压数值，通过算法以及控制器之后得到可输出到开关管的占空比，控制开关管的通断使系统工作在期望的电压点。当MPPT控制器输出高电平，开关管T导通时，二极管D连接到电池板正极承受反压而截止，电容C向负载供电，电池板的电压全部加到电感L两端，电感电流线性增长，储存的磁场能量也渐渐增加，直到开关管被截止；当MPPT控制器输出低电平，开关管T截止时，由于L的自身特性，会改变L的电压极性，这样电池板与电感串联经二极管D流向负载，电感电流线性减小，直到开关管T被导通。

根据图2的流程图以及图3的电路图，搭建了图4的仿真模型图，其中S-Function部分就是对粒子群算法的代码实现。

图4光伏阵列输出功率曲线图中可以得到最终的输出功率为122W，并且最终光伏电池输出功率稳定在了这个最大功率点处，输出稳定震荡较小，可以基本得出本发明方法可以有效解决多峰值最大功率跟踪的问题，克服传统方法陷入局部优的缺点，实现最大功率的跟踪问题。

本发明为解决其技术问题，所提供的一种基于粒子群的光伏电池多峰值最大功率跟踪系统，该系统与上述的方法对应，具体可参照上述方法。该系统包括：

初始化模块，用于根据光伏电池阵列的开路电压U_oc，在[0,U_oc]上选取m个点，作为粒子i的初始电压V_i(k)，i∈[1,2,3,…,m]，k表示迭代的次数，初始电压时k＝0,m为大于1的正整数；初始化粒子种群的参数，包括了粒子权重ω的取值，自我学习因子c₁与社会学习因子c₂的设定范围；

最优位置计算模块，用于计算各个粒子的适应度值，适应度值等于光伏电池阵列的输出功率；根据各个粒子的适应度值寻找种群的个体最优位置与全局最优位置；

粒子种群更新模块，用于更新各个粒子的位置与速度，得到新的下一代粒子种群的适应度值；

最优位置更新模块，用于分别比较下一代粒子种群的个体最优位置与全局最优位置与本次更新之前的个体最优位置与全局最优位置，以分别更新个体最优位置与全局最优位置，更新的准则：计算当前位置的适应度值后与本次更新之前的最优适应度值比较，将较大的适应度值对应的位置作为该粒子的新的个体最优位置，将所有的新的个体最优位置粒子中最大的适应度值对应的位置作为全局最优位置；

最优位置更新模块，用于返回粒子种群更新模块，直至达到收敛条件或者达到最大迭代次数，将最后的全局最优位置作为光伏系统的最大功率点所对应的工作电压。

进一步地，在本发明的基于粒子群的光伏电池多峰值最大功率跟踪系统中，所述的粒子权重ω取值是线性递减权重，计算公式为：

其中c为介于0至1之间的调节系数，k为当前迭代的次数，maxgen为最大迭代次数，ω_min为最小惯性权重，ω_max为最大惯性权重，f_i ^k为第k代中的第i个粒子的适应度值，

和

分别为第k代中的第i个粒子的最小和最大适应度值，

为第k代中的第i个粒子的权重。

进一步地，在本发明的基于粒子群的光伏电池多峰值最大功率跟踪系统中，所述更新各个粒子的位置与速度度具体是指：

在n维的搜索空间中，由m个粒子组成的种群记为X＝(x₁,...,x_i，...,x_m)，第i个粒子所在的位置为x_i＝(x_i1，x_i2，...x_in)^T，速度为v_i＝(v_i1，v_i2，...v_in)^T，第i个粒子在空间搜索过程中搜索到的自身最优位置P_i＝(p_i1，p_i2，...p_id)^T，全局最优位置P_g＝(p_g1，p_g2，...p_gd)^T，每个粒子代表一个解，粒子的位置与速度以迭代的方式进行更新，公式为：

是粒子i在第k次迭代的速度，是第d维的速度分量；

是粒子i在第k次的位置，

是

第d维的位置分量；

是粒子i在第k次迭代的个体最优位置，

是

第d维的个体最优位置分量；

是粒子群在第k次迭代的全局最优位置，

是

第d维的最优位置分量，

为第k代中的第i个粒子的权重，r₁与r₂是在[0,1]区间内均匀分布的随机数。

上面结合附图对本发明的实施例进行了描述，但是本发明并不局限于上述的具体实施方式，上述的具体实施方式仅仅是示意性的，而不是限制性的，本领域的普通技术人员在本发明的启示下，在不脱离本发明宗旨和权利要求所保护的范围情况下，还可做出很多形式，这些均属于本发明的保护之内。

Claims (8)

1.一种基于粒子群的光伏电池多峰值最大功率跟踪方法，其特征在于，包括：

S1、根据光伏电池阵列的开路电压U_oc，在[0,U_oc]上选取m个点，作为粒子i的初始电压V_i(k)，i∈[1,2,3,…,m]，k表示迭代的次数，初始电压时k＝0,m为大于1的正整数；初始化粒子种群的参数，包括了粒子权重ω的取值，自我学习因子c₁与社会学习因子c₂的设定范围；

S2、计算各个粒子的适应度值，适应度值等于光伏电池阵列的输出功率；根据各个粒子的适应度值寻找种群的个体最优位置与全局最优位置；

S3、更新各个粒子的位置与速度，得到新的下一代粒子种群的适应度值；

S4、分别比较下一代粒子种群的个体最优位置与全局最优位置与本次更新之前的个体最优位置与全局最优位置，以分别更新个体最优位置与全局最优位置，更新的准则：计算当前位置的适应度值后与本次更新之前的最优适应度值比较，将较大的适应度值对应的位置作为该粒子的新的个体最优位置，将所有的新的个体最优位置粒子中最大的适应度值对应的位置作为全局最优位置；

S5、返回步骤S3，直至达到收敛条件或者达到最大迭代次数，将最后的全局最优位置作为光伏系统的最大功率点所对应的工作电压；

所述的粒子权重ω取值是线性递减权重，计算公式为：

其中c为介于0至1之间的调节系数，k为当前迭代的次数，maxgen为最大迭代次数，ω_min为最小惯性权重，ω_max为最大惯性权重，

为第k代中的第i个粒子的适应度值，

和

分别为第k代中的第i个粒子的最小和最大适应度值，

为第k代中的第i个粒子的权重。

2.根据权利要求1所述的基于粒子群的光伏电池多峰值最大功率跟踪方法，其特征在于，所述的在[0,U_oc]上选取m个点的选取的要求是线性均匀选取或者是将粒子均匀分布在理论值处，所述的理论值数值等于0.8*U_oc。

3.根据权利要求1所述的基于粒子群的光伏电池多峰值最大功率跟踪方法，其特征在于，所述的学习因子c₁及c₂的设定范围是(0,2)，学习因子c₁及c₂在(0,2)区间内随机取值。

4.根据权利要求1所述的基于粒子群的光伏电池多峰值最大功率跟踪方法，其特征在于，所述更新各个粒子的位置与速度具体是指：

在n维的搜索空间中，由m个粒子组成的种群记为X＝(x₁,...,x_i，...,x_m)，第i个粒子所在的位置为x_i＝(x_i1，x_i2，...x_in)^T，速度为v_i＝(v_i1，v_i2，...v_in)^T，第i个粒子在空间搜索过程中搜索到的自身最优位置P_i＝(p_i1，p_i2，...p_id)^T，全局最优位置P_g＝(p_g1，p_g2，...p_gd)^T，每个粒子代表一个解，粒子的位置与速度以迭代的方式进行更新，公式为：

是粒子i在第k次迭代的速度，是

第d维的速度分量；

是粒子i在第k次的位置，

是

第d维的位置分量；

是粒子i在第k次迭代的个体最优位置，

是

第d维的个体最优位置分量；

是粒子群在第k次迭代的全局最优位置，

是

第d维的最优位置分量，

为第k代中的第i个粒子的权重，r₁与r₂是在[0,1]区间内均匀分布的随机数。

5.根据权利要求1所述的基于粒子群的光伏电池多峰值最大功率跟踪方法，其特征在于，所述达到收敛条件是指：

判断所有粒子所对应的适应度值最大与最小的差值，如果差值小于预设差值，表示达到收敛条件，此时粒子群算法执行完毕，否则未达到收敛条件。

6.根据权利要求1所述的基于粒子群的光伏电池多峰值最大功率跟踪方法，其特征在于，适应度值的计算公式为：

是粒子i在第k次迭代的速度，代表光伏电池输出的电压值；

代表光伏电池当前输出电压测得的电流；

是粒子的目标函数值代表当前光伏电池输出电压对应的输出功率，即适应度值。

7.一种基于粒子群的光伏电池多峰值最大功率跟踪系统，其特征在于，包括：

初始化模块，用于根据光伏电池阵列的开路电压U_oc，在[0,U_oc]上选取m个点，作为粒子i的初始电压V_i(k)，i∈[1,2,3,…,m]，k表示迭代的次数，初始电压时k＝0,m为大于1的正整数；初始化粒子种群的参数，包括了粒子权重ω的取值，自我学习因子c₁与社会学习因子c₂的设定范围；

最优位置计算模块，用于计算各个粒子的适应度值，适应度值等于光伏电池阵列的输出功率；根据各个粒子的适应度值寻找种群的个体最优位置与全局最优位置；

粒子种群更新模块，用于更新各个粒子的位置与速度，得到新的下一代粒子种群的适应度值；

最优位置更新模块，用于分别比较下一代粒子种群的个体最优位置与全局最优位置与本次更新之前的个体最优位置与全局最优位置，以分别更新个体最优位置与全局最优位置，更新的准则：计算当前位置的适应度值后与本次更新之前的最优适应度值比较，将较大的适应度值对应的位置作为该粒子的新的个体最优位置，将所有的新的个体最优位置粒子中最大的适应度值对应的位置作为全局最优位置；

最优位置更新模块，用于返回粒子种群更新模块，直至达到收敛条件或者达到最大迭代次数，将最后的全局最优位置作为光伏系统的最大功率点所对应的工作电压；

所述的粒子权重ω取值是线性递减权重，计算公式为：

其中c为介于0至1之间的调节系数，k为当前迭代的次数，maxgen为最大迭代次数，ω_min为最小惯性权重，ω_max为最大惯性权重，

为第k代中的第i个粒子的适应度值，

和

分别为第k代中的第i个粒子的最小和最大适应度值，

为第k代中的第i个粒子的权重。

8.根据权利要求7所述的基于粒子群的光伏电池多峰值最大功率跟踪系统，其特征在于，所述更新各个粒子的位置与速度具体是指：

在n维的搜索空间中，由m个粒子组成的种群记为X＝(x₁,...,x_i，...,x_m)，第i个粒子所在的位置为x_i＝(x_i1，x_i2，...x_in)^T，速度为v_i＝(v_i1，v_i2，...v_in)^T，第i个粒子在空间搜索过程中搜索到的自身最优位置P_i＝(p_i1，p_i2，...p_id)^T，全局最优位置P_g＝(p_g1，p_g2，...p_gd)^T，每个粒子代表一个解，粒子的位置与速度以迭代的方式进行更新，公式为：

是粒子i在第k次迭代的速度，是

第d维的速度分量；

是粒子i在第k次的位置，

是

第d维的位置分量；

是粒子i在第k次迭代的个体最优位置，

是

第d维的个体最优位置分量；

是粒子群在第k次迭代的全局最优位置，

是

第d维的最优位置分量，

为第k代中的第i个粒子的权重，r₁与r₂是在[0,1]区间内均匀分布的随机数。

Images