发明名称 |
一种基于小型风机的最大风能跟踪算法 |
摘要 |
本发明涉及风机控制算法领域。一种基于小型风机最大风能跟踪算法是改进的扰动观测算法与粒子群算法结合,改进的扰动算法不同于传统算法,根据电磁转矩与转速关系,对发电机电磁转矩提出了一种新的定义,即 <img file="dest_path_image002.GIF" wi="81" he="18" />。由这个公式可知,通过改变<img file="dest_path_image004.GIF" wi="8" he="11" />角可以找寻最大电磁转矩从而间接寻找最大功率。在对<img file="dest_path_image004a.GIF" wi="8" he="11" />角寻优的过程中利用了粒子群算法,提高了寻优的效率与精准度。这种方法特别适用于风速变化时能够迅速调整<img file="dest_path_image004aa.GIF" wi="8" he="11" />角,使输出功率恢复到最大值,从而实现最大风能追踪的快速性、精准性。 |
申请公布号 |
CN105069275A |
申请公布日期 |
2015.11.18 |
申请号 |
CN201510379656.3 |
申请日期 |
2015.06.29 |
申请人 |
广西大学 |
发明人 |
张涵;彭宇宁;李瀛;谢一飞;肖东裕;孙实超 |
分类号 |
G06F19/00(2011.01)I |
主分类号 |
G06F19/00(2011.01)I |
代理机构 |
成都九鼎天元知识产权代理有限公司 51214 |
代理人 |
陈科恒 |
主权项 |
一种基于小型风机的最大风能跟踪算法,其特征在于:按照以下步骤,步骤1,测定发电机初始运行状态下的发电机转速ω<sub>m</sub>,定义发电机电磁转矩T<sub>e</sub>满足公式:T<sub>e</sub>=tan(θ)ω<sub>m</sub><sup>2</sup>式1;步骤2,在式1参数θ稳态区域范围内任意选择微粒的初始位置<img file="FDA0000748534800000011.GIF" wi="246" he="96" />其中<img file="FDA0000748534800000012.GIF" wi="185" he="84" />初始速度V<sup>(0)</sup>=[V<sub>i</sub><sup>0</sup>],V<sub>θ</sub>范围为[V<sub>θmin</sub>,V<sub>θmax</sub>],初始速度V<sub>θ</sub>满足:V<sub>i</sub><sup>0</sup>=rand()·(V<sub>θmax</sub>‑V<sub>θmin</sub>)+V<sub>θmin</sub>式2;步骤3,采用所述发电机的电磁转矩T<sub>e</sub>的控制效果作为评价函数,适应值函数定义为f(x)=tanxω<sub>m</sub><sup>2</sup>式3;步骤4,定义所述微粒初始状态下局部最优值<img file="FDA0000748534800000013.GIF" wi="148" he="83" />和全局最优值<img file="FDA0000748534800000014.GIF" wi="175" he="85" />k=0;微粒位置<img file="FDA0000748534800000015.GIF" wi="66" he="77" />和速度<img file="FDA0000748534800000016.GIF" wi="68" he="76" />满足;速度方程组式3为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mi>V</mi><mi>i</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><mi>ω</mi><mo>·</mo><msubsup><mi>V</mi><mi>i</mi><mi>k</mi></msubsup><mo>+</mo><msub><mi>c</mi><mn>1</mn></msub><mrow><mo>(</mo><msubsup><mi>p</mi><mi>i</mi><mi>k</mi></msubsup><mo>-</mo><msubsup><mi>S</mi><mi>i</mi><mi>k</mi></msubsup><mo>)</mo></mrow><mo>+</mo><msub><mi>c</mi><mn>2</mn></msub><mrow><mo>(</mo><msubsup><mi>p</mi><mi>g</mi><mi>k</mi></msubsup><mo>-</mo><msubsup><mi>S</mi><mi>i</mi><mi>k</mi></msubsup><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000748534800000017.GIF" wi="848" he="90" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><mi>i</mi><mi>f</mi><mrow><mo>(</mo><msubsup><mi>V</mi><mi>i</mi><mi>k</mi></msubsup><mo>></mo><msub><mi>V</mi><mrow><mi>θ</mi><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>)</mo></mrow></mrow></mtd><mtd><mrow><msubsup><mi>V</mi><mi>i</mi><mi>k</mi></msubsup><mo>=</mo><msub><mi>V</mi><mrow><mi>θ</mi><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><mi>i</mi><mi>f</mi><mrow><mo>(</mo><msubsup><mi>V</mi><mi>i</mi><mi>k</mi></msubsup><mo><</mo><msub><mi>V</mi><mrow><mi>θ</mi><mi>m</mi><mi>i</mi><mi>n</mi></mrow></msub><mo>)</mo></mrow></mrow></mtd><mtd><mrow><msubsup><mi>V</mi><mi>i</mi><mi>k</mi></msubsup><mo>=</mo><msub><mi>V</mi><mrow><mi>θ</mi><mi>min</mi></mrow></msub></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000748534800000018.GIF" wi="606" he="180" /></maths>位置方程组式4为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>S</mi><mi>i</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><msubsup><mi>V</mi><mi>i</mi><mi>k</mi></msubsup><mo>+</mo><msubsup><mi>S</mi><mi>i</mi><mi>k</mi></msubsup></mrow>]]></math><img file="FDA0000748534800000019.GIF" wi="345" he="82" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><mi>i</mi><mi>f</mi><mrow><mo>(</mo><msubsup><mi>S</mi><mi>i</mi><mi>k</mi></msubsup><mo>></mo><msub><mi>S</mi><mrow><mi>θ</mi><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>)</mo></mrow></mrow></mtd><mtd><mrow><msubsup><mi>S</mi><mi>i</mi><mi>k</mi></msubsup><mo>=</mo><msub><mi>S</mi><mrow><mi>θ</mi><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><mi>i</mi><mi>f</mi><mrow><mo>(</mo><msubsup><mi>S</mi><mi>i</mi><mi>k</mi></msubsup><mo><</mo><msub><mi>S</mi><mrow><mi>θ</mi><mi>min</mi></mrow></msub><mo>)</mo></mrow></mrow></mtd><mtd><mrow><msubsup><mi>S</mi><mi>i</mi><mi>k</mi></msubsup><mo>=</mo><msub><mi>S</mi><mrow><mi>θ</mi><mi>min</mi></mrow></msub></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA00007485348000000110.GIF" wi="602" he="176" /></maths>其中,i∈[0,N‑1],N为微粒个数,i代表第i个微粒,K为迭代次数,C<sub>1</sub>、C<sub>2</sub>为加速常数,ω为惯性权重;步骤5,更新所述微粒的局部最优值和全局最优值;由步骤4确定出迭代后的微粒的位置<img file="FDA00007485348000000111.GIF" wi="93" he="77" />和速度<img file="FDA00007485348000000112.GIF" wi="127" he="77" />利用式5和式6更新微粒的局部最优值和全局最优值;<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msubsup><mi>p</mi><mi>i</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><mfenced open = '{' close = ''><mtable><mtr><mtd><msubsup><mi>p</mi><mi>i</mi><mi>k</mi></msubsup></mtd><mtd><mrow><mi>f</mi><mrow><mo>(</mo><msubsup><mi>s</mi><mi>i</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow><mo>≤</mo><mi>f</mi><mrow><mo>(</mo><msubsup><mi>p</mi><mi>i</mi><mi>k</mi></msubsup><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><msubsup><mi>s</mi><mi>i</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup></mtd><mtd><mrow><mi>f</mi><mrow><mo>(</mo><msubsup><mi>s</mi><mi>i</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow><mo>≥</mo><mi>f</mi><mrow><mo>(</mo><msubsup><mi>p</mi><mi>i</mi><mi>k</mi></msubsup><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000748534800000021.GIF" wi="581" he="243" /></maths>式5<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msubsup><mi>p</mi><mi>g</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><mfenced open = '{' close = ''><mtable><mtr><mtd><msubsup><mi>p</mi><mi>g</mi><mi>k</mi></msubsup></mtd><mtd><mrow><mi>f</mi><mrow><mo>(</mo><msubsup><mi>s</mi><mi>i</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow><mo>≤</mo><mi>f</mi><mrow><mo>(</mo><msubsup><mi>p</mi><mi>g</mi><mi>k</mi></msubsup><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><msubsup><mi>s</mi><mi>i</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup></mtd><mtd><mrow><mi>f</mi><mrow><mo>(</mo><msubsup><mi>s</mi><mi>i</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow><mo>≥</mo><mi>f</mi><mrow><mo>(</mo><msubsup><mi>p</mi><mi>g</mi><mi>k</mi></msubsup><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000748534800000022.GIF" wi="590" he="258" /></maths>式6步骤6:通过<img file="FDA0000748534800000023.GIF" wi="68" he="77" />和<img file="FDA0000748534800000024.GIF" wi="65" he="85" />由式1解得对应电磁转矩<img file="FDA0000748534800000025.GIF" wi="89" he="70" />通过<img file="FDA0000748534800000026.GIF" wi="94" he="83" />和<img file="FDA0000748534800000027.GIF" wi="96" he="92" />由式1解得对应电磁转矩<img file="FDA0000748534800000028.GIF" wi="114" he="75" />判断是否满足算法终止条件,定义<img file="FDA0000748534800000029.GIF" wi="309" he="91" />终止条件为:E<E<sub>limit</sub>;若满足,则算法停止,输出优化解θ=θ<sub>opt</sub>,否则将步骤5结果代入步骤4;步骤7,由所述参数θ=θ<sub>opt</sub>通过式1得到电磁转矩T<sub>e</sub>,通过矢量控制方法调节发电机输出最大的功率。 |
地址 |
530004 广西壮族自治区南宁市大学东路100号 |