发明名称 一种光伏阵列最大功率点滑模跟踪控制方法及系统
摘要 本发明为一种光伏阵列最大功率点滑模跟踪控制方法及系统。本方法根据光伏阵列的输出电压和电流得到开关控制信号的最大功率点跟踪控制方法具体为:步骤Ⅰ对电压、电流信号进行滑动平均滤波预处理,以减小采样测量干扰对最大功率点跟踪控制精度的影响,同时提高最大功率点的跟踪速度。步骤Ⅱ计算基于不完全偏微分的滑模函数S2的值,步骤Ⅲ对开关控制信号进行迟滞处理,根据S2的值与设定值δ关系,得到0或1的开关控制信号u。本系统单片机内存储基于不完全偏微分滑模面函数计算程序和迟滞处理程序,实现调整控制DC-DC变换电路的等效负载使之与光伏阵列内阻匹配,快速、精确跟踪光伏阵列的最大功率,结构简单,成本低,易于实现。
申请公布号 CN104460819B 申请公布日期 2016.04.13
申请号 CN201410814066.4 申请日期 2014.12.23
申请人 桂林电子科技大学 发明人 党选举;杨阳;姜辉;伍锡如;张向文;朱国魂;李爽
分类号 G05F1/67(2006.01)I 主分类号 G05F1/67(2006.01)I
代理机构 桂林市持衡专利商标事务所有限公司 45107 代理人 欧阳波
主权项 一种光伏阵列最大功率点滑模跟踪控制方法,包括有多块太阳能电池板的光伏阵列经DC‑DC降压变换电路连接负载,光伏阵列的输出端连接电压和电流采样电路,电压和电流采样电路接入单片机,单片机存储最大功率点跟踪控制方法的程序,单片机的信号输出接入DC‑DC降压变换电路;电压和电流采样电路监测光伏阵列的输出电压和电流,单片机根据实时电压和电流信号、按最大功率点跟踪控制方法的程序计算得到开关控制信号,单片机将所得开关控制信号送入DC‑DC降压变换电路,控制其中开关管的通断,调整脉宽调制波的占空比,控制DC‑DC降压变换电路的等效负载使之与光伏阵列内阻匹配,以跟踪控制光伏阵列的最大功率点;其特征在于根据光伏阵列的输出电压和电流得到开关控制信号的最大功率点滑模跟踪控制方法具体步骤如下:步骤Ⅰ、滑动平均滤波对电压、电流信号预处理对采集到的电压、电流信号u<sub>pv</sub>、i<sub>pv</sub>进行滑动平均滤波,即:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mover><mi>u</mi><mo>&OverBar;</mo></mover><mrow><mi>p</mi><mi>v</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><msub><mi>u</mi><mrow><mi>p</mi><mi>v</mi></mrow></msub><mo>(</mo><mrow><mi>k</mi><mo>-</mo><mi>N</mi><mo>+</mo><mn>1</mn></mrow><mo>)</mo><mo>+</mo><mo>...</mo><mo>.</mo><mo>+</mo><msub><mi>u</mi><mrow><mi>p</mi><mi>v</mi></mrow></msub><mo>(</mo><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo><mo>+</mo><msub><mi>u</mi><mrow><mi>p</mi><mi>v</mi></mrow></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow><mo>/</mo><mi>N</mi></mrow>]]></math><img file="FDA0000859135880000011.GIF" wi="1244" he="95" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mover><mi>i</mi><mo>&OverBar;</mo></mover><mrow><mi>p</mi><mi>v</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><msub><mi>i</mi><mrow><mi>p</mi><mi>v</mi></mrow></msub><mo>(</mo><mrow><mi>k</mi><mo>-</mo><mi>N</mi><mo>+</mo><mn>1</mn></mrow><mo>)</mo><mo>+</mo><mo>...</mo><mo>.</mo><mo>+</mo><msub><mi>i</mi><mrow><mi>p</mi><mi>v</mi></mrow></msub><mo>(</mo><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo><mo>+</mo><msub><mi>i</mi><mrow><mi>p</mi><mi>v</mi></mrow></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow><mo>/</mo><mi>N</mi></mrow>]]></math><img file="FDA0000859135880000012.GIF" wi="1182" he="95" /></maths>其中,k是离散采样时间、u<sub>pv</sub>(k)为第k次电压采样值、i<sub>pv</sub>(k)为第k次电流采样值电流、<img file="FDA0000859135880000013.GIF" wi="357" he="87" />为采样电压、电流第k次滑动平均滤波后的输出值,N为每次求平均的采样值个数;步骤Ⅱ、基于不完全偏微分的滑模面函数基于不完全偏微分的滑模面函数为<img file="FDA0000859135880000014.GIF" wi="566" he="95" />用步骤Ⅰ所得的<img file="FDA0000859135880000015.GIF" wi="358" he="87" />计算<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><msup><mi>&Delta;i</mi><mo>&prime;</mo></msup><mrow><mi>p</mi><mi>v</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msub><msup><mi>a&Delta;i</mi><mo>&prime;</mo></msup><mrow><mi>p</mi><mi>v</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>a</mi><mo>)</mo></mrow><mrow><mo>(</mo><msub><mover><mi>i</mi><mo>&OverBar;</mo></mover><mrow><mi>p</mi><mi>v</mi></mrow></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>-</mo><msub><mover><mi>i</mi><mo>&OverBar;</mo></mover><mrow><mi>p</mi><mi>v</mi></mrow></msub><mo>(</mo><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000859135880000016.GIF" wi="1110" he="94" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><msup><mi>&Delta;u</mi><mo>&prime;</mo></msup><mrow><mi>p</mi><mi>v</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msub><msup><mi>a&Delta;u</mi><mo>&prime;</mo></msup><mrow><mi>p</mi><mi>v</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>a</mi><mo>)</mo></mrow><mrow><mo>(</mo><msub><mover><mi>u</mi><mo>&OverBar;</mo></mover><mrow><mi>p</mi><mi>v</mi></mrow></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>-</mo><msub><mover><mi>u</mi><mo>&OverBar;</mo></mover><mrow><mi>p</mi><mi>v</mi></mrow></msub><mo>(</mo><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000859135880000017.GIF" wi="1166" he="95" /></maths>a=T/(T+Δt),T为惯性系数,Δt为采样周期,所得△i'<sub>pv</sub>和△u'<sub>pv</sub>值代入基于不完全偏微分的滑模面函数式,计算S2的值;步骤Ⅲ、迟滞处理开关控制信号对开关控制信号进行迟滞处理,即得到开关控制信号u如下式:<img file="FDA0000859135880000021.GIF" wi="886" he="240" />式中:δ为迟滞量,0.0004≤δ≤0.0006,<img file="FDA0000859135880000022.GIF" wi="75" he="86" />是S2的导数。
地址 541004 广西壮族自治区桂林市金鸡路1号