发明名称 一种基于非线性扩张状态观测器的永磁同步电机混沌控制方法
摘要 基于非线性扩张状态观测器的永磁同步电机混沌控制方法,包括:建立永磁同步电机系统的混沌模型,初始化系统状态以及相关控制参数;通过坐标变换,将永磁同步电机混沌模型转变为更适宜非线性扩张状态观测器设计的Brunovsky标准形式;设计非线性扩张状态观测器,用于估计不可测系统状态和参数扰动;根据非线性扩张状态观测器估计的系统状态和参数扰动,设计自适应滑模变结构控制器,改善滑模控制中的抖振问题,并保证系统的混沌状态快速稳定收敛至零点。
申请公布号 CN104201941A 申请公布日期 2014.12.10
申请号 CN201410306896.6 申请日期 2014.06.30
申请人 浙江工业大学 发明人 陈强;翟双坡;王晶
分类号 H02P6/00(2006.01)I;H02P23/12(2006.01)I;H02P25/02(2006.01)I 主分类号 H02P6/00(2006.01)I
代理机构 杭州天正专利事务所有限公司 33201 代理人 王兵;黄美娟
主权项 基于非线性扩张状态观测器的永磁同步电机混沌控制方法,包括以下步骤:步骤1,建立如式(1)所示的永磁同步电机系统的混沌模型,初始化系统状态以及相关控制参数;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><mfrac><mrow><mi>d</mi><msub><mover><mi>i</mi><mo>~</mo></mover><mi>d</mi></msub></mrow><mi>dt</mi></mfrac><mo>=</mo><mo>-</mo><msub><mover><mi>i</mi><mo>~</mo></mover><mi>d</mi></msub><mo>+</mo><mover><mi>&omega;</mi><mo>~</mo></mover><msub><mover><mi>i</mi><mo>~</mo></mover><mi>q</mi></msub><mo>+</mo><msub><mover><mi>u</mi><mo>~</mo></mover><mi>d</mi></msub></mtd></mtr><mtr><mtd><mfrac><mrow><mi>d</mi><msub><mover><mi>i</mi><mo>~</mo></mover><mi>q</mi></msub></mrow><mi>dt</mi></mfrac><mo>=</mo><mo>-</mo><msub><mover><mi>i</mi><mo>~</mo></mover><mi>q</mi></msub><mo>-</mo><mover><mi>&omega;</mi><mo>~</mo></mover><msub><mover><mi>i</mi><mo>~</mo></mover><mi>d</mi></msub><mo>+</mo><mi>&gamma;</mi><mover><mi>&omega;</mi><mo>~</mo></mover><mo>+</mo><msub><mover><mi>u</mi><mo>~</mo></mover><mi>q</mi></msub></mtd></mtr><mtr><mtd><mfrac><mrow><mi>d</mi><mover><mi>&omega;</mi><mo>~</mo></mover></mrow><mi>dt</mi></mfrac><mo>=</mo><mi>&sigma;</mi><mrow><mo>(</mo><msub><mover><mi>i</mi><mo>~</mo></mover><mi>q</mi></msub><mo>-</mo><mover><mi>&omega;</mi><mo>~</mo></mover><mo>)</mo></mrow><mo>-</mo><msub><mover><mi>T</mi><mo>~</mo></mover><mi>L</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000530314150000011.GIF" wi="1121" he="444" /></maths>其中,<img file="FDA0000530314150000012.GIF" wi="239" he="78" />为状态变量,分别表示直轴和交轴定子电流以及转子角频率;<img file="FDA0000530314150000013.GIF" wi="58" he="70" />和<img file="FDA0000530314150000014.GIF" wi="53" he="74" />表示直轴和交轴的定子电压;<img file="FDA0000530314150000015.GIF" wi="58" he="81" />为外部扭矩;σ和γ为常值参数;步骤2,通过坐标变换,将永磁同步电机混沌模型转变为更适宜非线性扩张状态观测器设计的Brunovsky标准形式,具体是:2.1,令<img file="FDA0000530314150000016.GIF" wi="464" he="84" />则式(1)可以等效为<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mn>1</mn></msub><mo>=</mo><mi>&sigma;</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>2</mn></msub><mo>-</mo><msub><mi>x</mi><mn>1</mn></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mn>2</mn></msub><mo>=</mo><mo>-</mo><msub><mi>x</mi><mn>2</mn></msub><mo>-</mo><msub><mi>x</mi><mn>1</mn></msub><msub><mi>x</mi><mn>3</mn></msub><mo>+</mo><mi>&gamma;</mi><msub><mi>x</mi><mn>1</mn></msub><mo>+</mo><mi>u</mi></mtd></mtr><mtr><mtd><msub><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mn>3</mn></msub><mo>=</mo><mo>-</mo><msub><mi>x</mi><mn>3</mn></msub><mo>+</mo><msub><mi>x</mi><mn>1</mn></msub><msub><mi>x</mi><mn>2</mn></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000530314150000017.GIF" wi="1196" he="235" /></maths>其中,x<sub>1</sub>,x<sub>2</sub>,x<sub>3</sub>为系统状态且x<sub>2</sub>,x<sub>3</sub>不可测,σ和γ为未知参数,u为控制信号,<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mover><mi>u</mi><mo>~</mo></mover><mi>d</mi></msub><mo>=</mo><msub><mover><mi>u</mi><mo>~</mo></mover><mi>q</mi></msub><mo>=</mo><msub><mi>T</mi><mi>L</mi></msub><mo>=</mo><mn>0</mn><mo>;</mo></mrow>]]></math><img file="FDA0000530314150000018.GIF" wi="353" he="70" /></maths>2.2,为便于控制器设计,将式(1)分解为如下两个子系统:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mn>1</mn></msub><mo>=</mo><mi>&sigma;</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>2</mn></msub><mo>-</mo><msub><mi>x</mi><mn>1</mn></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mn>2</mn></msub><mo>=</mo><mo>-</mo><msub><mi>x</mi><mn>2</mn></msub><mo>-</mo><msub><mi>x</mi><mn>1</mn></msub><msub><mi>x</mi><mn>3</mn></msub><mo>+</mo><mi>&gamma;</mi><msub><mi>x</mi><mn>1</mn></msub><mo>+</mo><mi>u</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000530314150000019.GIF" wi="1195" he="158" /></maths>和<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mn>3</mn></msub><mo>=</mo><mo>-</mo><msub><mi>x</mi><mn>3</mn></msub><mo>+</mo><msub><mi>x</mi><mn>1</mn></msub><msub><mi>x</mi><mn>2</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00005303141500000110.GIF" wi="1121" he="81" /></maths>其中,式(4)可以认为是式(2)的内动态方程,即:当x<sub>1</sub>,x<sub>2</sub>收敛至零点时,有<img file="FDA00005303141500000111.GIF" wi="176" he="68" />成立,从而x<sub>3</sub>也可以渐近收敛至零点;因此,本发明的控制目的为:设计控制器u,使得式(3)中的两个状态x<sub>1</sub>和x<sub>2</sub>收敛至零点;2.3,设<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>y</mi><mn>1</mn></msub><mo>=</mo><msub><mi>x</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>y</mi><mn>2</mn></msub><mo>=</mo><mi>&sigma;</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>2</mn></msub><mo>-</mo><msub><mi>x</mi><mn>1</mn></msub><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="FDA0000530314150000021.GIF" wi="1108" he="157" /></maths>则式(3)可转变为如下所示的Brunovsky标准形式:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mover><mi>y</mi><mo>&CenterDot;</mo></mover><mn>1</mn></msub><mo>=</mo><msub><mi>y</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><msub><mover><mi>y</mi><mo>&CenterDot;</mo></mover><mn>2</mn></msub><mo>=</mo><mi>a</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>+</mo><mi>bu</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000530314150000022.GIF" wi="1097" he="157" /></maths>其中,a(x)=σ[‑x<sub>2</sub>‑x<sub>1</sub>x<sub>3</sub>+γx<sub>1</sub>‑σ(x<sub>2</sub>‑x<sub>1</sub>)],b=σ;2.4,令a<sub>0</sub>=a(x)+Δbu,Δb=b‑b<sub>0</sub>,其中b<sub>0</sub>为b的估计值,可根据经验给定;基于扩张状态观测器的设计思想,通过定义扩展状态y<sub>3</sub>=a<sub>0</sub>,则式(6)可以改写为以下等效形式:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mover><mi>y</mi><mo>&CenterDot;</mo></mover><mn>1</mn></msub><mo>=</mo><msub><mi>y</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><msub><mover><mi>y</mi><mo>&CenterDot;</mo></mover><mn>2</mn></msub><mo>=</mo><msub><mi>y</mi><mn>3</mn></msub><mo>+</mo><msub><mi>b</mi><mn>0</mn></msub><mi>u</mi></mtd></mtr><mtr><mtd><msub><mover><mi>y</mi><mo>&CenterDot;</mo></mover><mn>3</mn></msub><mo>=</mo><mi>h</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000530314150000023.GIF" wi="1047" he="235" /></maths>其中,<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><mi>h</mi><mo>=</mo><msub><mover><mi>a</mi><mo>&CenterDot;</mo></mover><mn>0</mn></msub><mo>;</mo></mrow>]]></math><img file="FDA0000530314150000024.GIF" wi="140" he="71" /></maths>步骤3,设计非线性扩张状态观测器,用于估计不可测系统状态和参数扰动;令z<sub>i</sub>,i=1,2,3,分别为式(7)中状态变量y<sub>i</sub>的观测值,定义观测误差为e<sub>oi</sub>=z<sub>i</sub>‑y<sub>i</sub>,则设计非线性扩张状态观测器表达式为:<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mover><mi>z</mi><mo>&CenterDot;</mo></mover><mn>1</mn></msub><mo>=</mo><msub><mi>z</mi><mn>2</mn></msub><mo>-</mo><msub><mi>&beta;</mi><mn>1</mn></msub><msub><mi>e</mi><mrow><mi>o</mi><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mover><mi>z</mi><mo>&CenterDot;</mo></mover><mn>2</mn></msub><mo>=</mo><msub><mi>z</mi><mn>3</mn></msub><mo>-</mo><msub><mi>&beta;</mi><mn>2</mn></msub><mi>fal</mi><mrow><mo>(</mo><msub><mi>e</mi><mrow><mi>o</mi><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>&alpha;</mi><mn>1</mn></msub><mo>,</mo><mi>&delta;</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>b</mi><mn>0</mn></msub><mi>u</mi></mtd></mtr><mtr><mtd><msub><mover><mi>z</mi><mo>&CenterDot;</mo></mover><mn>3</mn></msub><mo>=</mo><mo>-</mo><msub><mi>&beta;</mi><mn>3</mn></msub><mi>fal</mi><mrow><mo>(</mo><msub><mi>e</mi><mrow><mi>o</mi><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>&alpha;</mi><mn>2</mn></msub><mo>,</mo><mi>&delta;</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000530314150000025.GIF" wi="1172" he="235" /></maths>其中,β<sub>1</sub>,β<sub>2</sub>,β<sub>3</sub>&gt;0为观测器增益;fal(·)为原点附近具有线性段的连续幂次函数,表达式为:<maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><mi>fal</mi><mrow><mo>(</mo><msub><mi>e</mi><mrow><mi>o</mi><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>&alpha;</mi><mi>i</mi></msub><mo>,</mo><mi>&delta;</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mfrac><msub><mi>e</mi><mrow><mi>o</mi><mn>1</mn></mrow></msub><msup><mi>&delta;</mi><mrow><mn>1</mn><mo>-</mo><msub><mi>&alpha;</mi><mi>i</mi></msub></mrow></msup></mfrac></mtd><mtd><mo>|</mo><msub><mi>e</mi><mrow><mi>o</mi><mn>1</mn></mrow></msub><mo>|</mo><mo>&le;</mo><mi>&delta;</mi></mtd></mtr><mtr><mtd><msup><mrow><mo>|</mo><msub><mi>e</mi><mrow><mi>o</mi><mn>1</mn></mrow></msub><mo>|</mo></mrow><msub><mi>&alpha;</mi><mi>i</mi></msub></msup><mi>sign</mi><mrow><mo>(</mo><msub><mi>e</mi><mrow><mi>o</mi><mn>1</mn></mrow></msub><mo>)</mo></mrow></mtd><mtd><mo>|</mo><msub><mi>e</mi><mrow><mi>o</mi><mn>1</mn></mrow></msub><mo>|</mo><mo>></mo><mi>&delta;</mi></mtd></mtr></mtable></mfenced><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2,3</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000530314150000026.GIF" wi="1438" he="218" /></maths>其中,δ&gt;0表示线性段的区间长度,0&lt;α<sub>i</sub>&lt;1;通过选择合适的参数β<sub>i</sub>,fal(·)函数可以保证z<sub>i</sub>→y<sub>i</sub>,i=1,2,3.即:观测误差可以收敛到|y<sub>i</sub>‑z<sub>i</sub>|≤d<sub>i</sub>,其中d<sub>i</sub>&gt;0为很小的正数.步骤4,根据非线性扩张状态观测器估计的系统状态和参数扰动,设计自适应滑模变结构控制器;。4.1,为将系统状态x<sub>1</sub>和x<sub>2</sub>稳定到原点,设计基于滑模变结构方法设计自适应控制器u,其中滑模面设计如式(10)所示:s=y<sub>2</sub>+λ<sub>1</sub>y<sub>1</sub>.        (10)s的一阶导数为:<maths num="0012" id="cmaths0012"><math><![CDATA[<mrow><mfenced open='' close=''><mtable><mtr><mtd><mover><mi>s</mi><mo>&CenterDot;</mo></mover><mo>=</mo><msub><mover><mi>y</mi><mo>&CenterDot;</mo></mover><mn>2</mn></msub><mo>+</mo><msub><mi>&lambda;</mi><mn>1</mn></msub><msub><mover><mi>y</mi><mo>&CenterDot;</mo></mover><mn>1</mn></msub></mtd></mtr><mtr><mtd><mo>=</mo><msub><mi>y</mi><mn>3</mn></msub><msub><mrow><mo>+</mo><mi>b</mi></mrow><mn>0</mn></msub><mi>u</mi><mo>+</mo><msub><mi>&lambda;</mi><mn>1</mn></msub><msub><mi>y</mi><mn>2</mn></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000530314150000031.GIF" wi="1058" he="147" /></maths>其中,λ<sub>1</sub>&gt;0为控制参数;4.2,由式(11),基于扩张状态观测器(8)的普通滑模控制器(SMC+ESO)设计为<maths num="0013" id="cmaths0013"><math><![CDATA[<mrow><msup><mi>u</mi><mo>*</mo></msup><mo>=</mo><mfrac><mn>1</mn><msub><mi>b</mi><mn>0</mn></msub></mfrac><mrow><mo>(</mo><mo>-</mo><msub><mi>z</mi><mn>3</mn></msub><mo>-</mo><msub><mi>&lambda;</mi><mn>1</mn></msub><msub><mi>z</mi><mn>2</mn></msub><mo>-</mo><msup><mi>k</mi><mo>*</mo></msup><mi>sign</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000530314150000032.GIF" wi="1221" he="141" /></maths>其中,k<sup>*</sup>&gt;0满足k<sup>*</sup>≥d<sub>3</sub>+λ<sub>1</sub>d<sub>2</sub>;4.3,由于式(12)中观测误差上界d<sub>2</sub>和d<sub>3</sub>难以准确获得,因此,k<sup>*</sup>往往无法精确得到;为解决上述问题,设计基于扩张状态观测器的自适应滑模控制器(ASMC+ESO),其具体表达形式为:<maths num="0014" id="cmaths0014"><math><![CDATA[<mrow><mi>u</mi><mo>=</mo><mfrac><mn>1</mn><msub><mi>b</mi><mn>0</mn></msub></mfrac><mrow><mo>(</mo><mo>-</mo><msub><mi>z</mi><mn>3</mn></msub><mo>-</mo><msub><mi>&lambda;</mi><mn>1</mn></msub><msub><mi>z</mi><mn>2</mn></msub><mo>-</mo><mi>ksign</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000530314150000033.GIF" wi="1221" he="140" /></maths>其中,k=k(t)为控制器参数,其参数自适应律如下所示:<maths num="0015" id="cmaths0015"><math><![CDATA[<mrow><mover><mi>k</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>k</mi><mi>m</mi></msub><mo>|</mo><mi>s</mi><mo>|</mo><mi>sign</mi><mrow><mo>(</mo><mo>|</mo><mi>s</mi><mo>|</mo><mo>-</mo><mo>&Element;</mo><mo>)</mo></mrow><mo>,</mo></mtd><mtd><mi>k</mi><mo>></mo><mi>&mu;</mi></mtd></mtr><mtr><mtd><mi>&mu;</mi><mo>,</mo></mtd><mtd><mi>k</mi><mo>&le;</mo><mi>&mu;</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>14</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000530314150000034.GIF" wi="1233" he="156" /></maths>其中,k<sub>m</sub>&gt;0,μ&gt;0为很小的正常数,用于保证k&gt;0;4.4,设计李雅普诺夫函数<img file="FDA0000530314150000041.GIF" wi="373" he="134" />则可以证明式(6)中的所有信号均是一致有界的;同时,系统跟踪误差e可以在有限时间内收敛至平衡点e=0。
地址 310014 浙江省杭州市下城区潮王路18号