发明名称 一种基于优化灰色预测补偿的永磁同步电机转速控制方法
摘要 本发明涉及一种基于优化灰色预测补偿的永磁同步电机转速控制方法,针对噪声和小扰动等不确定因素对永磁同步电机转速控制精度带来的影响,采用灰色系统理论建立灰色预测模型,对因干扰项的存在而产生的转速偏差进行超前预测,进而对实际控制过程予以实时补偿,提高永磁同步电机的转速跟踪性能和鲁棒性。而传统灰色方法在初值选取方面存在固有缺陷,首先提出了一种优化后的灰色预测模型,然后运用于电机的转速控制中。本发明通过推导灰色模型,得到了一种优化后的灰色预测模型,在矢量控制的基础上用此灰色模型预测永磁同步电机未来的转速偏差,进行超前预测补偿,提高了电机转速的控制精度和鲁棒性。
申请公布号 CN104242744B 申请公布日期 2016.08.24
申请号 CN201410479991.6 申请日期 2014.09.19
申请人 西北工业大学 发明人 骆光照;张蓉;张莎;涂文聪;韩复振;孟锦豪
分类号 H02P6/08(2016.01)I;H02P21/00(2016.01)I;H02P21/18(2016.01)I 主分类号 H02P6/08(2016.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 一种基于优化灰色预测补偿的永磁同步电机转速控制方法,其特征在于步骤如下:步骤1:建立永磁同步电机矢量控制模型<img file="FDA0000573118920000011.GIF" wi="657" he="128" />其中:i<sub>d</sub>,i<sub>q</sub>分别为d‑q坐标下的电流;ψ<sub>f</sub>为转子永磁体磁链;L<sub>d</sub>,L<sub>q</sub>为d‑q坐标系下的定子电感值;p为极对数;F为粘滞系数;J为转动惯量;ω为转子电角速度;T<sub>m</sub>为负载转矩;T<sub>e</sub>为电磁转矩;U为噪声及扰动条件下转速的不确定部分;所述电磁转矩方程:T<sub>e</sub>=p[ψ<sub>f</sub>i<sub>q</sub>‑(L<sub>d</sub>‑L<sub>q</sub>)i<sub>d</sub>i<sub>q</sub>];步骤2:将永磁同步电机运行过程中,理想条件和实际情况下的转速误差作为灰色预测序列的原始序列为Δω<sup>(0)</sup>,通过累加公式<img file="FDA0000573118920000012.GIF" wi="462" he="150" />得到累加序列Δω<sup>(1)</sup>={Δω<sup>(1)</sup>(1),Δω<sup>(1)</sup>(2),…,Δω<sup>(1)</sup>(n)};步骤3:建立微分方程:<img file="FDA0000573118920000013.GIF" wi="457" he="136" />通过最小二乘求解公式<maths num="0001"><math><![CDATA[<mrow><mfenced open='[' close=']'><mtable><mtr><mtd><mover><mi>a</mi><mo>^</mo></mover></mtd></mtr><mtr><mtd><mover><mi>b</mi><mo>^</mo></mover></mtd></mtr></mtable></mfenced><mo>=</mo><msup><mrow><mo>(</mo><msup><mi>B</mi><mi>T</mi></msup><mi>B</mi><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msup><mi>B</mi><mi>T</mi></msup><mi>Y</mi><mo>,</mo></mrow>]]></math><img file="FDA0000573118920000014.GIF" wi="424" he="225" /></maths>得到待估参数a,b,且a为发展系数,b为灰色作用量,其中,<maths num="0002"><math><![CDATA[<mrow><mi>B</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mo>-</mo><msup><mi>Z</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mtd><mtd><mn>1</mn></mtd></mtr><mtr><mtd><mo>-</mo><msup><mi>Z</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mtd><mtd><mn>1</mn></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>-</mo><msup><mi>Z</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000573118920000015.GIF" wi="426" he="336" /></maths>Z<sup>(1)</sup>(k+1)为预测模型的背景值,Y=[Δω<sup>(0)</sup>(2),Δω<sup>(0)</sup>(3),…,Δω<sup>(0)</sup>(n)]<sup>T</sup>,<maths num="0003"><math><![CDATA[<mrow><msup><mi>Z</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msup><mi>&Delta;&omega;</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><msup><mi>&Delta;&omega;</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><mn>2</mn></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000573118920000016.GIF" wi="691" he="139" /></maths>步骤4:根据估计值,得到灰色算法微分方程的解为:<img file="FDA0000573118920000017.GIF" wi="449" he="129" />所述C为常数项,为:<maths num="0004"><math><![CDATA[<mrow><mi>C</mi><mo>=</mo><mfrac><mrow><msup><mi>&Delta;&omega;</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mfrac><mi>b</mi><mi>a</mi></mfrac></mrow><msup><mi>e</mi><mrow><mo>-</mo><mi>at</mi></mrow></msup></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000573118920000018.GIF" wi="378" he="187" /></maths>步骤5:令t=k,并通过累减Δω<sup>(0)</sup>=Δω<sup>(1)</sup>(k)‑Δω<sup>(1)</sup>(k‑1),还原为原始序列,则:<maths num="0005"><math><![CDATA[<mrow><msup><mi>&Delta;&omega;</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><msup><mi>&Delta;&omega;</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mfrac><mi>b</mi><mi>a</mi></mfrac><mo>)</mo></mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>e</mi><mi>a</mi></msup><mo>)</mo></mrow><msup><mi>e</mi><mrow><mo>-</mo><mi>a</mi><mrow><mo>(</mo><mi>i</mi><mo>-</mo><mi>k</mi><mo>)</mo></mrow></mrow></msup><mo>,</mo></mrow>]]></math><img file="FDA00005731189200000110.GIF" wi="813" he="122" /></maths>此时取k=i,则:<maths num="0006"><math><![CDATA[<mrow><msup><mi>&Delta;&omega;</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msup><mo>=</mo><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><msup><mi>&Delta;&omega;</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>-</mo><mfrac><mi>b</mi><mi>a</mi></mfrac><mo>)</mo></mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>e</mi><mi>a</mi></msup><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000573118920000019.GIF" wi="648" he="124" /></maths>其中,i=1,2,…,n;步骤6:对原始序列的n个数的预测值,设预测序列为Δω',将预测的此序列进行累加得到:Δω'<sup>(1)</sup>={Δω'<sup>(1)</sup>(1),Δω'<sup>(1)</sup>(2),…,Δω'<sup>(1)</sup>(n)},取最近的Δω'<sup>(1)</sup>(n)值作为初值,同时a,b值也相应改变为a'、b',得到:<img file="FDA0000573118920000021.GIF" wi="853" he="134" />其中,j=n+1,n+2,…;步骤7:得到用一个分段函数表示的灰色预测模型为:<maths num="0007"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><msup><mi>&Delta;&omega;</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><msup><mi>&Delta;&omega;</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>-</mo><mfrac><mi>b</mi><mi>a</mi></mfrac><mo>)</mo></mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>e</mi><mi>a</mi></msup><mo>)</mo></mrow></mtd><mtd><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi></mtd></mtr><mtr><mtd><msup><mi>&Delta;&omega;</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><msup><mi>&Delta;&omega;</mi><mrow><mo>&prime;</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></msup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><mfrac><msup><mi>b</mi><mo>&prime;</mo></msup><msup><mi>a</mi><mo>&prime;</mo></msup></mfrac><mo>)</mo></mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>e</mi><msup><mi>a</mi><mo>&prime;</mo></msup></msup><mo>)</mo></mrow><msup><mi>e</mi><mrow><mo>-</mo><msup><mi>a</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><mi>j</mi><mo>-</mo><mi>n</mi><mo>)</mo></mrow></mrow></msup></mtd><mtd><mi>j</mi><mo>=</mo><mi>n</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>n</mi><mo>+</mo><mn>2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000573118920000022.GIF" wi="1400" he="281" /></maths>步骤8:以灰色预测模型,对未来转速误差进行超前预测,设此预测补偿量为Δω<sub>p</sub>,则得到的转速修正量为ω=ω<sub>c</sub>+Δω<sub>p</sub>;在矢量控制双闭环控制策略下,将此转速修正量作为控制输入量嵌入控制系统中,实现转速的超前预测补偿控制。
地址 710072 陕西省西安市友谊西路127号