主权项 |
一种光伏阵列最大功率点滑模跟踪控制方法,包括有多块太阳能电池板的光伏阵列经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>‾</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>‾</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>Δi</mi><mo>′</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Δi</mi><mo>′</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>‾</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>‾</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>Δu</mi><mo>′</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Δu</mi><mo>′</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>‾</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>‾</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的导数。 |