发明名称 基于带遗忘因子的递推最小二乘光伏电池参数辨识方法
摘要 本发明公开了一种基于带遗忘因子的递推最小二乘光伏电池模型参数辨识方法,包括以下步骤:建立光伏电池的递推最小二乘模型形式,确定待辨识参数;初始化参数估值、遗忘因子及协方差矩阵;获取光伏电池的实时输出电压电流,更新参数估值和目标函数值;目标函数满足预设阈值时输出此时的参数估值,即为光伏电池参数最优取值。本发明将光伏电池的四参数模型转化为递推最小二乘模型形式,屏蔽了光伏电池内部结构及系统误差对模型的影响,通过迭代得到最优参数取值。它实现简单,能减少计算量,减少数据在计算机中占用的内存,提高了辨识速度。遗忘因子能够强调新数据的作用,逐渐遗忘旧数据的作用,使模型具有较高的精度,数值稳定性好。
申请公布号 CN106100582A 申请公布日期 2016.11.09
申请号 CN201610531729.0 申请日期 2016.07.07
申请人 国网青海省电力公司;国网青海省电力公司电力科学研究院;华北电力大学(保定) 发明人 杨军;徐岩;靳伟佳
分类号 H02S50/10(2014.01)I 主分类号 H02S50/10(2014.01)I
代理机构 北京科亿知识产权代理事务所(普通合伙) 11350 代理人 汤东凤
主权项 一种基于带遗忘因子的递推最小二乘光伏电池参数辨识方法,所述方法用于辨识光伏并网发电系统中m串n并的光伏阵列中的光伏电池参数;其特征在于:包括以下步骤:步骤1:建立光伏电池的递推最小二乘模型,确定待辨识参数θ:所述光伏电池的递推最小二乘模型为:<maths num="0001"><math><![CDATA[<mrow><mi>&Delta;</mi><mi>I</mi><mo>=</mo><mfenced open = "(" close = ")"><mtable><mtr><mtd><msub><mi>U</mi><mi>L</mi></msub></mtd><mtd><msubsup><mi>U</mi><mi>L</mi><mn>2</mn></msubsup></mtd><mtd><msubsup><mi>U</mi><mi>L</mi><mn>3</mn></msubsup></mtd></mtr></mtable></mfenced><mfenced open = "(" close = ")"><mtable><mtr><mtd><msub><mi>a</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>a</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><msub><mi>a</mi><mn>3</mn></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001043832090000011.GIF" wi="1270" he="263" /></maths>且ΔI=I<sub>sc</sub>‑I<sub>L</sub>                       (2)式中,θ<sup>T</sup>=(a<sub>1</sub>,a<sub>2</sub>,a<sub>3</sub>)为待辨识的参数,I<sub>L</sub>、U<sub>L</sub>分别为光伏电池的实时输出电压和实时输出电流,I<sub>sc</sub>为光生电流;步骤2:初始化:由以下具体分步骤组成:步骤2‑1:定义输入观测量的长度n和遗忘因子λ,0<λ≤1;步骤2‑2:设置待辨识的参数估值的初始时刻取值<img file="FDA0001043832090000012.GIF" wi="234" he="87" />设置协方差矩阵P(t)初始时刻的取值P(0)=aI<sub>n</sub>,a实数,I<sub>n</sub>为n×n的单位矩阵;步骤3:获取光伏电池的实时输出电压U<sub>L</sub>和实时输出电流I<sub>L</sub>:U<sub>L</sub>=U<sub>L_array</sub>/m                      (3)I<sub>L</sub>=I<sub>L_array</sub>/n                       (4)式中,U<sub>L_array</sub>和I<sub>L_array</sub>F分别为所述光伏阵列的实时输出电压实时输出电流;步骤4:产生观测矩阵<img file="FDA0001043832090000014.GIF" wi="109" he="80" />及其转置矩阵<img file="FDA0001043832090000015.GIF" wi="170" he="78" /><img file="FDA0001043832090000016.GIF" wi="1238" he="103" />其中,y(t)是光伏电池输出电压当前时刻的观测值,y(t‑1)为光伏电池输出电压前一时刻的观测值,y(t‑n)为光伏电池输出电压前n时刻的观测值;步骤5:计算当前时刻观测矩阵<img file="FDA0001043832090000017.GIF" wi="106" he="79" />的协方差矩阵P(t):<img file="FDA0001043832090000013.GIF" wi="1550" he="214" />步骤6:计算上一时刻的增益矩阵F(t‑1),<img file="FDA0001043832090000026.GIF" wi="1564" he="104" />步骤7:更新光伏电池的参数估值<img file="FDA0001043832090000021.GIF" wi="139" he="87" /><img file="FDA0001043832090000022.GIF" wi="1542" he="109" />步骤8:计算当前时刻的预测残差:<img file="FDA0001043832090000023.GIF" wi="1262" he="103" />并计算当前时刻的目标函数:<img file="FDA0001043832090000024.GIF" wi="1382" he="150" />步骤9:判断当前时刻的目标函数J<sub>t</sub>(θ)是否小于预设阈值,如果否,转向步骤3,如果是,转向步骤10;步骤10:输出参数估值<img file="FDA0001043832090000025.GIF" wi="135" he="94" />将其作为光伏电池参数θ的辨识结果。
地址 810008 青海省西宁市城西区胜利路89号