发明名称 一种采用自抗扰控制技术的网络化无刷直流电机时延补偿和控制方法
摘要 一种采用自抗扰控制技术的网络化无刷直流电机时延补偿和控制方法,包括以下步骤:1)建立含有时变网络诱导时延的无刷直流电机控制系统模型。将网络化无刷直流电机控制系统描述为一个具有一步输入时滞的离散时间线性时变系统,进而将时变时延引起的系统不确定动态部分描述为系统的加性噪声;2)扩张状态观测器设计。本发明采用扩张状态观测器可将时变时延引起的不确定动态作为总和扰动的一部分进行估计;3)对网络化无刷直流电机控制系统中的时变时延项的补偿过程。本发明可以有效地将时变网络诱导时延引起的不确定动态用扩张状态观测器实时估计并补偿,此方法对时延引起的不确定性和系统内外扰动以及模型不确定性均具有很强的抑制能力。
申请公布号 CN104142627A 申请公布日期 2014.11.12
申请号 CN201410314689.5 申请日期 2014.07.03
申请人 浙江工业大学 发明人 张文安;刘凯;俞立
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 杭州斯可睿专利事务所有限公司 33241 代理人 王利强
主权项 一种采用自抗扰控制技术的网络化无刷直流电机时延补偿和控制方法,其特征在于:所述方法包括以下步骤:步骤1)建立含有时变网络诱导时延的无刷直流电机控制系统模型,将网络化无刷直流电机控制系统描述为一个具有一步输入时滞的离散时间线性时变系统,进而将时变时延引起的系统不确定动态部分描述为系统的加性噪声,包括以下过程:1.1)建立无刷直流电机控制系统的线性化传递函数模型无刷直流电机控制系统的二阶系统模型的传递函数为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfrac><mrow><mi>N</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mrow><mrow><mi>U</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mrow></mfrac><mo>=</mo><mfrac><mi>b</mi><mrow><msup><mi>s</mi><mn>2</mn></msup><mo>+</mo><mi>as</mi><mo>+</mo><mi>c</mi></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000532111380000016.GIF" wi="1136" he="132" /></maths>其中,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>a</mi><mo>=</mo><mfrac><mi>r</mi><mrow><mi>L</mi><mo>-</mo><mi>M</mi></mrow></mfrac><mo>,</mo><mi>b</mi><mo>=</mo><mfrac><mrow><mi>&pi;</mi><msub><mi>K</mi><mi>e</mi></msub></mrow><mrow><mn>30</mn><mi>J</mi><mrow><mo>(</mo><mi>L</mi><mo>-</mo><mi>M</mi><mo>)</mo></mrow></mrow></mfrac><mo>,</mo><mi>c</mi><mo>=</mo><mfrac><mrow><mn>2</mn><msubsup><mi>K</mi><mi>e</mi><mn>2</mn></msubsup></mrow><mrow><mi>J</mi><mrow><mo>(</mo><mi>L</mi><mo>-</mo><mi>M</mi><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000532111380000011.GIF" wi="959" he="150" /></maths>r是定子相绕组的电阻,L是绕组的自感,M是两相绕组间的互感,K<sub>e</sub>是电动势系数,J为电机转动惯量,N(s)是电机转速的拉普拉斯变换,U(s)为导通两相的支路电压的拉普拉斯变换;将式(1)所示的传递函数模型转化为如下的状态空间模型:<img file="FDA0000532111380000012.GIF" wi="1175" he="216" />其中,x<sub>1</sub>为直流电机的转速,x<sub>2</sub>为直流电机加速度,u为控制量,即导通两相的支路电压;1.2)获得时变网络诱导时延影响下的电机控制系统模型数据包在网络中传输时,存在从传感器到控制器和控制器到执行器之间的时延,用<img file="FDA0000532111380000013.GIF" wi="62" he="77" />和<img file="FDA0000532111380000014.GIF" wi="70" he="78" />分别表示测量信号从传感器传输到控制器所经历的时延和控制量从控制器传输到执行器的时延,那么控制回路总的网络诱导时延是<img file="FDA0000532111380000015.GIF" wi="285" he="77" />由于时延小于一个采样周期,直流电机在一个周期内的控制输入电压u(t)由两部分构成,一部分是由上一周期计算得到的控制输入电压u(k‑1),另一部分是当前周期计算得到的控制输入电压u(k),且具有以下形式:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mi>u</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo></mtd><mtd><msub><mi>t</mi><mi>k</mi></msub><mo>&le;</mo><mn>1</mn><mo>&lt;</mo><msub><mi>t</mi><mi>k</mi></msub><mo>+</mo><msub><mi>&tau;</mi><mi>k</mi></msub></mtd></mtr><mtr><mtd><mi>u</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo></mtd><mtd><msub><mi>t</mi><mi>k</mi></msub><mo>+</mo><msub><mi>&tau;</mi><mi>k</mi></msub><mo>&le;</mo><mi>t</mi><mo>&lt;</mo><msub><mi>t</mi><mi>k</mi></msub><mo>+</mo><mi>T</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000532111380000021.GIF" wi="1321" he="150" /></maths>其中,T是采样周期,t<sub>k</sub>表示第k个采样时刻;因此,根据式(2)和(3),离散化后的含有时变网络诱导时延的无刷直流电机控制系统模型为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><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>T</mi><msub><mi>x</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>x</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>aT</mi></mrow></msup><msub><mi>x</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mfrac><mn>1</mn><mi>a</mi></mfrac><mrow><mo>(</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>a</mi><mrow><mo>(</mo><mi>T</mi><mo>-</mo><msub><mi>&tau;</mi><mi>k</mi></msub><mo>)</mo></mrow></mrow></msup><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mrow><mo>(</mo><mi>bu</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>cx</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>-</mo><mfrac><mn>1</mn><mi>a</mi></mfrac><mrow><mo>(</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>aT</mi></mrow></msup><mo>-</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>a</mi><mrow><mo>(</mo><mi>T</mi><mo>-</mo><msub><mi>&tau;</mi><mi>k</mi></msub><mo>)</mo></mrow></mrow></msup><mo>)</mo></mrow><mrow><mo>(</mo><mi>bu</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><msub><mi>cx</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000532111380000022.GIF" wi="1512" he="354" /></maths>将e<sup>‑aT</sup>用1‑aT近似后,将式(4)化为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><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>T</mi><msub><mi>x</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><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>T</mi><mrow><mo>(</mo><mo>-</mo><msub><mi>ax</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mfrac><mn>1</mn><mi>aT</mi></mfrac></mrow><mrow><mo>(</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>a</mi><mrow><mo>(</mo><mi>T</mi><mo>-</mo><msub><mi>&tau;</mi><mi>k</mi></msub><mo>)</mo></mrow></mrow></msup><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mrow><mo>(</mo><mi>bu</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>cx</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>-</mo><mfrac><mn>1</mn><mi>a</mi></mfrac><mrow><mo>(</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>aT</mi></mrow></msup><mo>-</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>a</mi><mrow><mo>(</mo><mi>T</mi><mo>-</mo><msub><mi>&tau;</mi><mi>k</mi></msub><mo>)</mo></mrow></mrow></msup><mo>)</mo></mrow><mrow><mo>(</mo><mi>bu</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><msub><mi>cx</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000532111380000023.GIF" wi="1617" he="351" /></maths>将式(5)中由时变时延τ<sub>k</sub>引起的时变动态用一个新的状态变量x<sub>3</sub>(k)表示,即<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>x</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><msub><mi>ax</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mfrac><mi>b</mi><mi>aT</mi></mfrac><msup><mi>e</mi><mrow><mo>-</mo><mi>a</mi><mrow><mo>(</mo><mi>T</mi><mo>-</mo><msub><mi>&tau;</mi><mi>k</mi></msub><mo>)</mo></mrow></mrow></msup><mi>u</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mfrac><mn>1</mn><mi>aT</mi></mfrac><mrow><mo>(</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>a</mi><mrow><mo>(</mo><mi>T</mi><mo>-</mo><msub><mi>&tau;</mi><mi>k</mi></msub><mo>)</mo></mrow></mrow></msup><mo>-</mo><mn>1</mn><mo>)</mo></mrow><msub><mi>cx</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mfrac><mn>1</mn><mi>aT</mi></mfrac><mrow><mo>(</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>aT</mi></mrow></msup><mo>-</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>a</mi><mrow><mo>(</mo><mi>T</mi><mo>-</mo><msub><mi>&tau;</mi><mi>k</mi></msub><mo>)</mo></mrow></mrow></msup><mo>)</mo></mrow><mrow><mo>(</mo><mi>bu</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><msub><mi>cx</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000532111380000024.GIF" wi="1930" he="135" /></maths>,并令<img file="FDA0000532111380000025.GIF" wi="505" he="127" />由此将由式(5)表示的网络化无刷直流电机控制系统模型扩张成如下的三阶系统模型:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><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><msub><mi>Tx</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><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>T</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mfrac><mi>b</mi><mi>aT</mi></mfrac><mi>u</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow></mtd></mtr><mtr><mtd><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>Tw</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000532111380000026.GIF" wi="1350" he="295" /></maths>其中,x<sub>1</sub>(k+1)、x<sub>2</sub>(k+1)、x<sub>3</sub>(k+1)分别为电机转速x<sub>1</sub>(k)、电机加速度x<sub>2</sub>(k)、新扩张状态量x<sub>3</sub>(k)的下一采样时刻的值;步骤2)设计扩张状态观测器,用于估计无刷直流电机控制系统中由网络诱导时延引起的不确定性;步骤3)设计带有扩张状态观测器的网络化无刷直流电机自抗扰控制器,实现时变网络诱导时延的实时补偿和电机转速实时控制。
地址 310014 浙江省杭州市下城区朝晖六区潮王路18号