发明名称 一种适用于新能源发电机的不确定大时滞的补偿方法
摘要 本发明公开了一种适用于新能源发电机的不确定大时滞的补偿方法。通过Fibonacci序列和时间尺度的参数自适应调节方法建立线性ADRC控制器,输出并网逆变器的控制量;利用互相关方法和滑动窗数字信号处理方法,由并网逆变器的控制量和输出电流计算得到不确定大时滞;采用Smith预估器计算得到新能源发电机并网逆变器控制量,实现不确定大时滞的补偿和新能源发电机的自适应鲁棒控制。本发明可准确估计新能源发电机中的不确定大时滞,降低新能源发电机输出电流的超调量,减小调节时间,并提高新能源发电机的运行稳定性和鲁棒性。
申请公布号 CN104578143B 申请公布日期 2016.06.29
申请号 CN201410775666.4 申请日期 2014.12.15
申请人 浙江大学 发明人 张国月;齐冬莲;吴越;张建良
分类号 H02J3/38(2006.01)I 主分类号 H02J3/38(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 林超
主权项 一种适用于新能源发电机的不确定大时滞的补偿方法,其特征在于包括以下步骤:1)通过Fibonacci数列和时间尺度的参数自适应调节方法建立线性ADRC控制器,输出并网逆变器的控制量;2)利用互相关方法和滑动窗数字信号处理方法,由并网逆变器的控制量和输出电流计算得到不确定大时滞;3)采用Smith预估器计算得到新能源发电机的并网逆变器的控制量,实现不确定大时滞的补偿和新能源发电机的自适应鲁棒控制;所述的步骤1)中建立的线性ADRC控制器采用以下公式(1):<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>f</mi><mi>s</mi><mo>=</mo><mo>-</mo><msup><mi>r</mi><mn>3</mn></msup><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>-</mo><mi>v</mi><mo>)</mo></mrow><mo>-</mo><mn>3</mn><msup><mi>r</mi><mn>2</mn></msup><msub><mi>x</mi><mn>2</mn></msub><mo>-</mo><mn>3</mn><msub><mi>rx</mi><mn>3</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>x</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msub><mi>x</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mi>h</mi><mo>&CenterDot;</mo><msub><mi>x</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>x</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msub><mi>x</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mi>h</mi><mo>&CenterDot;</mo><msub><mi>x</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>x</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msub><mi>x</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mi>h</mi><mo>&CenterDot;</mo><mi>f</mi><mi>s</mi></mrow></mtd></mtr><mtr><mtd><mrow><mi>e</mi><mo>=</mo><msub><mi>z</mi><mn>1</mn></msub><mo>-</mo><mi>y</mi></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>z</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msub><mi>z</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mi>h</mi><mo>&CenterDot;</mo><mo>&lsqb;</mo><msub><mi>z</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>&beta;</mi><mn>01</mn></msub><mi>e</mi><mo>&rsqb;</mo></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>z</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msub><mi>z</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mi>h</mi><mo>&CenterDot;</mo><mo>&lsqb;</mo><msub><mi>z</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>&beta;</mi><mn>02</mn></msub><mi>e</mi><mo>&rsqb;</mo></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>z</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msub><mi>z</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mi>h</mi><mo>&CenterDot;</mo><mo>&lsqb;</mo><msub><mi>z</mi><mn>4</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>&beta;</mi><mn>03</mn></msub><mi>e</mi><mo>+</mo><mi>u</mi><mo>&rsqb;</mo></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>z</mi><mn>4</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msub><mi>z</mi><mn>4</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mi>h</mi><mo>&CenterDot;</mo><mrow><mo>(</mo><mo>-</mo><msub><mi>&beta;</mi><mn>04</mn></msub><mi>e</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>e</mi><mn>1</mn></msub><mo>=</mo><msub><mi>x</mi><mn>1</mn></msub><mo>-</mo><msub><mi>z</mi><mn>1</mn></msub><mo>,</mo><msub><mi>e</mi><mn>2</mn></msub><mo>=</mo><msub><mi>x</mi><mn>2</mn></msub><mo>-</mo><msub><mi>z</mi><mn>2</mn></msub><mo>,</mo><msub><mi>e</mi><mn>3</mn></msub><mo>=</mo><msub><mi>x</mi><mn>3</mn></msub><mo>-</mo><msub><mi>z</mi><mn>3</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><mi>U</mi><mo>=</mo><msub><mi>&beta;</mi><mn>1</mn></msub><mo>&CenterDot;</mo><msub><mi>e</mi><mn>1</mn></msub><mo>+</mo><msub><mi>&beta;</mi><mn>2</mn></msub><mo>&CenterDot;</mo><msub><mi>e</mi><mn>2</mn></msub><mo>+</mo><msub><mi>&beta;</mi><mn>3</mn></msub><mo>&CenterDot;</mo><msub><mi>e</mi><mn>3</mn></msub><mo>-</mo><msub><mi>z</mi><mn>4</mn></msub><mo>-</mo><msub><mi>f</mi><mrow><mi>k</mi><mi>d</mi></mrow></msub></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000971607650000011.GIF" wi="1125" he="933" /></maths>其中:v为新能源发电机的并网逆变器的输出电流的参考值,x<sub>1</sub>、x<sub>2</sub>、x<sub>3</sub>分别为v的跟踪值、一阶微分信号和二阶微分信号,h为采样步长,y为新能源发电机的并网逆变器的输出电流,r为快速因子,k表示第k步运算,z<sub>1</sub>、z<sub>2</sub>、z<sub>3</sub>分别为y的跟踪值、一阶微分信号、二阶微分信号,z<sub>4</sub>为新能源发电机受到的内部和外部的扰动总和w的估计值,β<sub>01</sub>、β<sub>02</sub>、β<sub>03</sub>、β<sub>04</sub>分别为第一、第二、第三、第四待调参数;e为新能源发电机的并网逆变器的输出电流跟踪值z<sub>1</sub>和新能源发电机的并网逆变器的输出电流y的差值,e<sub>1</sub>、e<sub>2</sub>、e<sub>3</sub>分别为x<sub>1</sub>与z<sub>1</sub>、x<sub>2</sub>与z<sub>2</sub>、x<sub>3</sub>与z<sub>3</sub>的差值,β<sub>1</sub>、β<sub>2</sub>、β<sub>3</sub>分别为第五、第六、第七待调参数;U为虚拟控制量;f<sub>kd</sub>为新能源发电机的并网逆变器已建模动态,u为新能源发电机的并网逆变器的控制量,fs表示三阶线性跟踪微分器。
地址 310058 浙江省杭州市西湖区余杭塘路866号