发明名称 基于串并联估计模型的柔性机械臂系统模糊控制方法
摘要 一种基于串并联估计模型的柔性机械臂系统模糊控制方法,包括:建立柔性机械臂伺服系统的动态模型,初始化系统状态、采样时间以及控制参数;采用模糊神经网络控制的自学与自适应能力,补偿机械臂系统中的未知非线性不确定项;进行通过动态面技术,在每一步设计中引入虚拟控制量,并依次通过一阶低通滤波器,避免传统反演控制法所带来的复杂度爆炸问题;同时,基于串并联估计模型,定义变量预测值并设计相关变化律,提高控制系统的鲁棒性。本发明提供一种能够有效改善柔性机械臂伺服系统控制性能的基于串并联估计模型的模糊自适应控制方法,实现系统的精确快速跟踪。
申请公布号 CN106113040A 申请公布日期 2016.11.16
申请号 CN201610569377.8 申请日期 2016.07.19
申请人 浙江工业大学 发明人 陈强;高灵捷;余梦梦
分类号 B25J9/16(2006.01)I 主分类号 B25J9/16(2006.01)I
代理机构 杭州斯可睿专利事务所有限公司 33241 代理人 王利强
主权项 一种基于串并联估计模型的柔性机械臂系统模糊控制方法,其特征在于:包括以下步骤:步骤1:建立柔性机械臂伺服系统动态模型,初始化系统状态、采样时间以及控制参数,过程如下:1.1柔性机械臂伺服系统动态模型的运动方程表达式为<maths num="0001"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>I</mi><mover><mi>q</mi><mo>&CenterDot;&CenterDot;</mo></mover><mo>+</mo><mi>K</mi><mrow><mo>(</mo><mi>q</mi><mo>-</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>+</mo><mi>M</mi><mi>g</mi><mi>L</mi><mi>s</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><mi>q</mi><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><mi>J</mi><mover><mi>&theta;</mi><mo>&CenterDot;&CenterDot;</mo></mover><mo>-</mo><mi>K</mi><mrow><mo>(</mo><mi>q</mi><mo>-</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>=</mo><mi>u</mi></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001052969080000011.GIF" wi="844" he="135" /></maths>其中,q与θ分别为机械臂连杆和电机的转动角度;<img file="FDA0001052969080000012.GIF" wi="30" he="53" />与<img file="FDA0001052969080000013.GIF" wi="35" he="55" />分别为机械臂连杆和电机的转动角加速度;g为重力加速度;I为连杆的惯量;J为电机的惯量;K为弹簧的刚度系数;M与L分别为连杆的质量与长度;u为控制信号;1.2定义:x<sub>1</sub>=q,<img file="FDA0001052969080000014.GIF" wi="166" he="54" />x<sub>3</sub>=θ,<img file="FDA0001052969080000015.GIF" wi="170" he="63" />式(1)改写为<maths num="0002"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mn>1</mn></msub><mo>=</mo><msub><mi>x</mi><mn>2</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mn>2</mn></msub><mo>=</mo><msub><mi>x</mi><mn>3</mn></msub><mo>+</mo><msub><mi>f</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mn>3</mn></msub><mo>=</mo><msub><mi>x</mi><mn>4</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mn>4</mn></msub><mo>=</mo><mfrac><mn>1</mn><mi>J</mi></mfrac><mi>u</mi><mo>+</mo><msub><mi>f</mi><mn>4</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mi>y</mi><mo>=</mo><msub><mi>x</mi><mn>1</mn></msub></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001052969080000016.GIF" wi="534" he="342" /></maths>其中,y为系统位置输出轨迹;<img file="FDA0001052969080000017.GIF" wi="29" he="62" />与<img file="FDA0001052969080000018.GIF" wi="31" he="56" />分别为机械臂连杆和电机的转动角速度;<maths num="0003"><math><![CDATA[<mrow><msub><mi>f</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><mfrac><mrow><mi>M</mi><mi>g</mi><mi>L</mi></mrow><mi>I</mi></mfrac><mi>sin</mi><mi> </mi><msub><mi>x</mi><mn>1</mn></msub><mo>-</mo><mfrac><mi>K</mi><mi>I</mi></mfrac><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>-</mo><msub><mi>x</mi><mn>3</mn></msub><mo>+</mo><mfrac><mi>I</mi><mi>K</mi></mfrac><msub><mi>x</mi><mn>3</mn></msub><mo>)</mo></mrow><mo>,</mo><msub><mi>f</mi><mn>4</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mi>K</mi><mi>J</mi></mfrac><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>-</mo><msub><mi>x</mi><mn>3</mn></msub><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0001052969080000019.GIF" wi="1358" he="95" /></maths>步骤2:针对式(1),计算控制系统位置跟踪误差,利用模糊系统逼近复杂非线性项,设计系统状态预测误差以及预测变量的变化律,设计虚拟控制量,并通过一阶低通滤波器输出,最后更新模糊系统权值,过程如下:2.1定义系统的位置跟踪误差s<sub>1</sub>=x<sub>1</sub>‑y<sub>r</sub>  (3)其中,y<sub>r</sub>为二阶可导期望位置轨迹;2.2设计虚拟控制量<maths num="0004"><math><![CDATA[<mrow><msubsup><mi>x</mi><mn>2</mn><mi>d</mi></msubsup><mo>=</mo><mo>-</mo><msub><mi>k</mi><mn>1</mn></msub><msub><mi>s</mi><mn>1</mn></msub><mo>+</mo><msub><mover><mi>y</mi><mo>&CenterDot;</mo></mover><mi>r</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00010529690800000110.GIF" wi="486" he="67" /></maths>其中,k<sub>1</sub>为常数且满足k<sub>1</sub>>0,<img file="FDA00010529690800000111.GIF" wi="50" he="56" />为期望速度轨迹;2.3定义一个新的变量<img file="FDA00010529690800000112.GIF" wi="75" he="55" />让虚拟控制量<img file="FDA00010529690800000113.GIF" wi="60" he="62" />通过时间常数为τ<sub>2</sub>的一阶低通滤波器,得<maths num="0005"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>&tau;</mi><mn>2</mn></msub><msubsup><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mn>2</mn><mi>c</mi></msubsup><mo>+</mo><msubsup><mi>x</mi><mn>2</mn><mi>c</mi></msubsup><mo>=</mo><msubsup><mi>x</mi><mn>2</mn><mi>d</mi></msubsup></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>x</mi><mn>2</mn><mi>c</mi></msubsup><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow><mo>=</mo><msubsup><mi>x</mi><mn>2</mn><mi>d</mi></msubsup><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001052969080000021.GIF" wi="479" he="135" /></maths>2.4定义滤波误差<img file="FDA0001052969080000022.GIF" wi="307" he="59" />为消除滤波误差对控制效果的影响,定义第一补偿信号z<sub>1</sub>,其变化律表达式为<maths num="0006"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mover><mi>z</mi><mo>&CenterDot;</mo></mover><mn>1</mn></msub><mo>=</mo><mo>-</mo><msub><mi>k</mi><mn>1</mn></msub><msub><mi>z</mi><mn>1</mn></msub><mo>+</mo><msub><mi>z</mi><mn>2</mn></msub><mo>+</mo><msub><mi>&chi;</mi><mn>2</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>z</mi><mn>1</mn></msub><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001052969080000023.GIF" wi="582" he="119" /></maths>其中,z<sub>2</sub>为第二补偿信号;2.5定义跟踪误差补偿信号v<sub>1</sub>=s<sub>1</sub>‑z<sub>1</sub>    (7)2.6定义误差变量<maths num="0007"><math><![CDATA[<mrow><msub><mi>s</mi><mn>2</mn></msub><mo>=</mo><msub><mi>x</mi><mn>2</mn></msub><mo>-</mo><msubsup><mi>x</mi><mn>2</mn><mi>c</mi></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001052969080000024.GIF" wi="398" he="63" /></maths>2.7为了逼近复杂的非线性不确定项f<sub>2</sub>(X),定义以下模糊系统<img file="FDA0001052969080000025.GIF" wi="550" he="64" />其中,<img file="FDA0001052969080000026.GIF" wi="53" he="55" />为理想权重;<img file="FDA0001052969080000027.GIF" wi="422" he="63" />ε<sub>2</sub>为模糊逼近误差,ε<sub>N2</sub>为逼近误差上界,满足|ε<sub>2</sub>|≤ε<sub>N2</sub>;<img file="FDA0001052969080000028.GIF" wi="1238" he="63" />的表达式为<img file="FDA0001052969080000029.GIF" wi="629" he="103" />其中,μ<sub>l</sub>(x<sub>j</sub>)为隶属度函数,其表达式为<img file="FDA00010529690800000210.GIF" wi="914" he="99" /><img file="FDA00010529690800000211.GIF" wi="54" he="43" />为常数,exp(·)为指数函数;2.8设计虚拟控制量<img file="FDA00010529690800000212.GIF" wi="814" he="65" />其中,k<sub>2</sub>为常数且满足k<sub>2</sub>>0,<img file="FDA00010529690800000213.GIF" wi="52" he="63" />为<img file="FDA00010529690800000214.GIF" wi="51" he="62" />的估计值;2.9定义一个新的变量<img file="FDA00010529690800000215.GIF" wi="75" he="59" />让虚拟控制量<img file="FDA00010529690800000216.GIF" wi="59" he="63" />通过时间常数为τ<sub>3</sub>的一阶低通滤波器,可得<maths num="0008"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>&tau;</mi><mn>3</mn></msub><msubsup><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mn>3</mn><mi>c</mi></msubsup><mo>+</mo><msubsup><mi>x</mi><mn>3</mn><mi>c</mi></msubsup><mo>=</mo><msubsup><mi>x</mi><mn>3</mn><mi>d</mi></msubsup></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>x</mi><mn>3</mn><mi>c</mi></msubsup><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow><mo>=</mo><msubsup><mi>x</mi><mn>3</mn><mi>d</mi></msubsup><mo>(</mo><mn>0</mn><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00010529690800000217.GIF" wi="478" he="135" /></maths>2.10定义滤波误差<img file="FDA00010529690800000218.GIF" wi="306" he="67" />为消除滤波误差对控制效果的影响,定义第二补偿信号z<sub>2</sub>,设计其变化律表达式<maths num="0009"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mover><mi>z</mi><mo>&CenterDot;</mo></mover><mn>2</mn></msub><mo>=</mo><mo>-</mo><msub><mi>k</mi><mn>2</mn></msub><msub><mi>z</mi><mn>2</mn></msub><mo>+</mo><msub><mi>z</mi><mn>3</mn></msub><mo>-</mo><msub><mi>z</mi><mn>1</mn></msub><mo>+</mo><msub><mi>&chi;</mi><mn>3</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>z</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001052969080000031.GIF" wi="710" he="118" /></maths>其中,z<sub>3</sub>为第三补偿信号;2.11定义预测误差<maths num="0010"><math><![CDATA[<mrow><msub><mi>z</mi><mrow><mn>2</mn><mi>N</mi><mi>N</mi></mrow></msub><mo>=</mo><msub><mi>x</mi><mn>2</mn></msub><mo>-</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mn>2</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>14</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001052969080000032.GIF" wi="451" he="63" /></maths>其中,<img file="FDA0001052969080000033.GIF" wi="54" he="55" />为系统状态x<sub>2</sub>的预测值;2.12设计串并联估计模型<img file="FDA0001052969080000034.GIF" wi="720" he="134" />其中,β<sub>2</sub>为常数且满足β<sub>2</sub>>0;2.13定义跟踪误差补偿信号v<sub>2</sub>=s<sub>2</sub>‑z<sub>2</sub>  (16)2.14设计模糊系统权重估计值<img file="FDA0001052969080000035.GIF" wi="48" he="64" />的调节规律为<img file="FDA0001052969080000036.GIF" wi="901" he="79" />其中,δ<sub>2</sub>与r<sub>z2</sub>为常数,且δ<sub>2</sub>>0,r<sub>z2</sub>>0,γ<sub>2</sub>为对称正定矩阵;2.15定义误差变量<maths num="0011"><math><![CDATA[<mrow><msub><mi>s</mi><mn>3</mn></msub><mo>=</mo><msub><mi>x</mi><mn>3</mn></msub><mo>-</mo><msubsup><mi>x</mi><mn>3</mn><mi>c</mi></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>18</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001052969080000037.GIF" wi="398" he="61" /></maths>2.16设计虚拟控制量<maths num="0012"><math><![CDATA[<mrow><msubsup><mi>x</mi><mn>4</mn><mi>d</mi></msubsup><mo>=</mo><msubsup><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mn>3</mn><mi>c</mi></msubsup><mo>-</mo><msub><mi>k</mi><mn>3</mn></msub><msub><mi>s</mi><mn>3</mn></msub><mo>-</mo><msub><mi>s</mi><mn>2</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>19</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001052969080000038.GIF" wi="564" he="64" /></maths>其中,k<sub>3</sub>为常数,且k<sub>3</sub>>0;2.17定义一个新的变量<img file="FDA0001052969080000039.GIF" wi="74" he="55" />让虚拟控制量<img file="FDA00010529690800000310.GIF" wi="59" he="62" />通过时间常数为τ<sub>4</sub>的一阶低通滤波器,可得<maths num="0013"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>&tau;</mi><mn>4</mn></msub><msubsup><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mn>4</mn><mi>c</mi></msubsup><mo>+</mo><msubsup><mi>x</mi><mn>4</mn><mi>c</mi></msubsup><mo>=</mo><msubsup><mi>x</mi><mn>4</mn><mi>d</mi></msubsup></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>x</mi><mn>4</mn><mi>c</mi></msubsup><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow><mo>=</mo><msubsup><mi>x</mi><mn>4</mn><mi>d</mi></msubsup><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>20</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00010529690800000311.GIF" wi="480" he="135" /></maths>2.18定义滤波误差<img file="FDA00010529690800000312.GIF" wi="307" he="63" />为消除滤波误差对控制效果的影响,定义第三补偿信号z<sub>3</sub>,设计其变化律表达式<maths num="0014"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mover><mi>z</mi><mo>&CenterDot;</mo></mover><mn>3</mn></msub><mo>=</mo><mo>-</mo><msub><mi>k</mi><mn>3</mn></msub><msub><mi>z</mi><mn>3</mn></msub><mo>+</mo><msub><mi>z</mi><mn>4</mn></msub><mo>-</mo><msub><mi>z</mi><mn>2</mn></msub><mo>+</mo><msub><mi>&chi;</mi><mn>4</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>z</mi><mn>3</mn></msub><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>21</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00010529690800000313.GIF" wi="710" he="120" /></maths>2.19定义跟踪误差补偿信号v<sub>3</sub>=s<sub>3</sub>‑z<sub>3</sub>  (22)步骤3:设计控制器输入,过程如下:3.1定义误差变量<maths num="0015"><math><![CDATA[<mrow><msub><mi>s</mi><mn>4</mn></msub><mo>=</mo><msub><mi>x</mi><mn>4</mn></msub><mo>-</mo><msubsup><mi>x</mi><mn>4</mn><mi>c</mi></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>23</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001052969080000041.GIF" wi="398" he="63" /></maths>3.2为了逼近复杂的非线性不确定项f<sub>4</sub>(X),定义以下模糊系统<maths num="0016"><math><![CDATA[<mrow><msub><mi>f</mi><mn>4</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mi>&theta;</mi><mn>4</mn><mo>*</mo></msubsup><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>&epsiv;</mi><mn>4</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>24</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001052969080000042.GIF" wi="542" he="63" /></maths>其中,<img file="FDA0001052969080000043.GIF" wi="55" he="55" />为理想权重;ε<sub>4</sub>为模糊逼近误差,ε<sub>N4</sub>为逼近误差上界,满足|ε<sub>4</sub>|≤ε<sub>N4</sub>;3.3设计控制器输入为u<img file="FDA0001052969080000044.GIF" wi="838" he="69" />其中,k<sub>4</sub>为常数且满足k<sub>4</sub>>0,<img file="FDA0001052969080000045.GIF" wi="52" he="64" />为<img file="FDA0001052969080000046.GIF" wi="51" he="55" />的估计值;3.4定义补偿信号z<sub>4</sub>,设计其变化律表达式<maths num="0017"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mover><mi>z</mi><mo>&CenterDot;</mo></mover><mn>4</mn></msub><mo>=</mo><mo>-</mo><msub><mi>k</mi><mn>4</mn></msub><msub><mi>z</mi><mn>4</mn></msub><mo>-</mo><msub><mi>z</mi><mn>3</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>z</mi><mn>4</mn></msub><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>26</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001052969080000047.GIF" wi="502" he="119" /></maths>3.5定义预测误差<maths num="0018"><math><![CDATA[<mrow><msub><mi>z</mi><mrow><mn>4</mn><mi>N</mi><mi>N</mi></mrow></msub><mo>=</mo><msub><mi>x</mi><mn>4</mn></msub><mo>-</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mn>4</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>27</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001052969080000048.GIF" wi="454" he="63" /></maths>其中,,<img file="FDA0001052969080000049.GIF" wi="57" he="55" />为系统状态x<sub>4</sub>的预测值;3.6设计串并联估计模型<img file="FDA00010529690800000410.GIF" wi="701" he="127" />其中,β<sub>4</sub>为常数且满足β<sub>4</sub>>0;3.7定义跟踪误差补偿信号v<sub>4</sub>=s<sub>4</sub>‑z<sub>4</sub>  (29)3.8设计模糊系统权重估计值<img file="FDA00010529690800000411.GIF" wi="52" he="63" />的调节规律为<img file="FDA00010529690800000412.GIF" wi="894" he="79" />其中,δ<sub>4</sub>,r<sub>z4</sub>为常数,且δ<sub>4</sub>>0,r<sub>z4</sub>>0,γ<sub>4</sub>为对称正定矩阵;步骤4:设计李雅普诺夫函数<maths num="0019"><math><![CDATA[<mrow><mi>V</mi><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>4</mn></msubsup><msubsup><mi>v</mi><mi>i</mi><mn>2</mn></msubsup><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mrow><mo>(</mo><msubsup><mover><mi>&theta;</mi><mo>~</mo></mover><mn>2</mn><mi>T</mi></msubsup><msubsup><mi>&gamma;</mi><mn>2</mn><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msub><mover><mi>&theta;</mi><mo>~</mo></mover><mn>2</mn></msub><mo>+</mo><msubsup><mover><mi>&theta;</mi><mo>~</mo></mover><mn>4</mn><mi>T</mi></msubsup><msubsup><mi>&gamma;</mi><mn>4</mn><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msub><mover><mi>&theta;</mi><mo>~</mo></mover><mn>4</mn></msub><mo>+</mo><msub><mi>r</mi><mrow><mi>z</mi><mn>2</mn></mrow></msub><msubsup><mi>z</mi><mrow><mn>2</mn><mi>N</mi><mi>N</mi></mrow><mn>2</mn></msubsup><mo>+</mo><msub><mi>r</mi><mrow><mi>z</mi><mn>4</mn></mrow></msub><msubsup><mi>z</mi><mrow><mn>4</mn><mi>N</mi><mi>N</mi></mrow><mn>2</mn></msubsup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>31</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00010529690800000413.GIF" wi="1454" he="92" /></maths>对式(31)进行求导得:<maths num="0020"><math><![CDATA[<mrow><mover><mi>V</mi><mo>&CenterDot;</mo></mover><mo>=</mo><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>4</mn></msubsup><msub><mi>v</mi><mi>i</mi></msub><msub><mover><mi>v</mi><mo>&CenterDot;</mo></mover><mi>i</mi></msub><mo>+</mo><mrow><mo>(</mo><mo>-</mo><msubsup><mover><mi>&theta;</mi><mo>~</mo></mover><mn>2</mn><mi>T</mi></msubsup><msubsup><mi>&gamma;</mi><mn>2</mn><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msub><mover><mover><mi>&theta;</mi><mo>^</mo></mover><mo>&CenterDot;</mo></mover><mn>2</mn></msub><mo>-</mo><msubsup><mover><mi>&theta;</mi><mo>~</mo></mover><mn>4</mn><mi>T</mi></msubsup><msubsup><mi>&gamma;</mi><mn>4</mn><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msub><mover><mover><mi>&theta;</mi><mo>^</mo></mover><mo>&CenterDot;</mo></mover><mn>4</mn></msub><mo>+</mo><msub><mi>r</mi><mrow><mi>z</mi><mn>2</mn></mrow></msub><msub><mover><mi>z</mi><mo>&CenterDot;</mo></mover><mrow><mn>2</mn><mi>N</mi><mi>N</mi></mrow></msub><msub><mi>z</mi><mrow><mn>2</mn><mi>N</mi><mi>N</mi></mrow></msub><mo>+</mo><msub><mi>r</mi><mrow><mi>z</mi><mn>4</mn></mrow></msub><msub><mover><mi>z</mi><mo>&CenterDot;</mo></mover><mrow><mn>4</mn><mi>N</mi><mi>N</mi></mrow></msub><msub><mi>z</mi><mrow><mn>4</mn><mi>N</mi><mi>N</mi></mrow></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>32</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00010529690800000414.GIF" wi="1621" he="95" /></maths>如果<img file="FDA00010529690800000415.GIF" wi="146" he="56" />则判定系统是稳定的。
地址 310014 浙江省杭州市下城区潮王路18号浙江工业大学科技处