发明名称 具有随机丢包的非线性双时标系统网络化控制方法
摘要 一种具有随机丢包的非线性双时标系统网络化控制方法,属于复杂系统网络化控制技术领域。该方法融合模糊逻辑、奇异摄动技术以及随机系统理论,在统一模型框架下,综合考虑被控系统自身的非线性双时标并存特性和网络通信引发的随机丢包问题,建立网络化非线性双时标系统的离散模糊奇异摄动模型,提出随机模糊状态反馈鲁棒控制方法。优点在于,解决现有网络控制方法无法消除复杂双时标系统快变量引起的系统失稳或稳态误差问题,大幅提高NTTSSs的网络化控制性能。
申请公布号 CN104570727A 申请公布日期 2015.04.29
申请号 CN201410650629.0 申请日期 2014.11.14
申请人 冶金自动化研究设计院 发明人 陈金香
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 北京华谊知识产权代理有限公司 11207 代理人 刘月娥
主权项 一种具有随机丢包的非线性双时标系统网络化控制方法,其特征在于:工艺步骤如下:步骤1、根据被控NTTSSs的动力学方程,建立NTTSSs的连续时间模糊奇异摄动模型;将NTTSSs的变化缓慢或能够直接测量的状态变量看作为慢变量,小参数相关或变化较快的状态变量看作为快变量,采用扇区非线性方法,建立NTTSSs的连续时间模糊奇异摄动模型:规则i:如果ξ<sub>1</sub>(t)是φ<sub>i1</sub>,...,ξ<sub>g</sub>(t)是φ<sub>ig</sub>,那么<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>E</mi><mi>&epsiv;</mi></msub><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>A</mi><mi>ci</mi></msub><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>B</mi><mi>ci</mi></msub><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000608375790000011.GIF" wi="549" he="71" /></maths>y(t)=C<sub>ci</sub>x(t)            (1)其中,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>E</mi><mi>&epsiv;</mi></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>I</mi><mrow><mi>n</mi><mo>&times;</mo><mi>n</mi></mrow></msub></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>&epsiv;</mi><msub><mi>I</mi><mrow><mi>m</mi><mo>&times;</mo><mi>m</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>,</mo><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>x</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>x</mi><mi>f</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000608375790000012.GIF" wi="865" he="135" /></maths>x(t)∈R<sup>(n+m)</sup>为状态变量,x<sub>s</sub>(t)∈R<sup>n</sup>为慢状态变量,x<sub>f</sub>(t)∈R<sup>m</sup>为快状态变量,u(t)∈R<sup>q</sup>为控制输入,y(t)∈R<sup>l</sup>为系统输出,φ<sub>i1</sub>,...,φ<sub>ig</sub>(i=1,2,...,r)均为模糊集合,ξ<sub>1</sub>(t),...,ξ<sub>g</sub>(t)为可测量的系统变量,A<sub>ci</sub>,B<sub>ci</sub>,C<sub>ci</sub>为合适维数矩阵,ε是奇异摄动参数;步骤2、建立被控NTTSSs的离散时间模糊奇异摄动模型控制系统中的传感器和执行器均采用时间驱动方式,且二者采用相同的采样时间T<sub>s</sub>,在零阶保持器的作用下,将以上连续时间模型(1),离散化为如下离散时间模糊奇异摄动模型:规则i:如果ξ<sub>1</sub>(k)是φ<sub>i1</sub>,...,ξ<sub>g</sub>(k)是φ<sub>ig</sub>,那么x(k+1)=E<sub>ε</sub>A<sub>i</sub>x(k)+E<sub>ε</sub>B<sub>i</sub>u(k)y(k)=C<sub>i</sub>x(k)               (2)其中,x(k)∈R<sup>(n+m)</sup>为状态变量,<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>x</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>x</mi><mi>f</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000608375790000013.GIF" wi="388" he="132" /></maths>x<sub>s</sub>(k)∈R<sup>n</sup>为慢状态变量,x<sub>f</sub>(k)∈R<sup>m</sup>为快状态变量,u(k)∈R<sup>q</sup>为控制输入,y(k)∈R<sup>l</sup>为系统输出,C<sub>i</sub>=C<sub>ci</sub>,A<sub>i</sub>,B<sub>i</sub>的值可采用如下matlab命令获得:sys=ss(A<sub>ci</sub>,B<sub>ci</sub>,C<sub>ci</sub>,0)             (3)c2d(sys,T<sub>s</sub>,′zoh′)               (4)给定[x(k);u(k)],应用标准模糊推理方法,得到全局离散时间模糊奇异摄动模型:x(k+1)=E<sub>ε</sub>A(μ)x(k)+E<sub>ε</sub>B(μ)u(k)y(k)=C(μ)x(k)                (5)其中,隶属度函数<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>&mu;</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>&xi;</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>w</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>&xi;</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>r</mi></munderover><msub><mi>w</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>&xi;</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow></mrow></mfrac><mo>,</mo><msub><mi>w</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>&xi;</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Pi;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>g</mi></munderover><msub><mi>&phi;</mi><mi>ij</mi></msub><mrow><mo>(</mo><msub><mi>&xi;</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>,</mo><mo>,</mo></mrow>]]></math><img file="FDA0000608375790000014.GIF" wi="1004" he="198" /></maths>φ<sub>ij</sub>(ξ<sub>j</sub>(k))为ξ<sub>j</sub>(k)在φ<sub>ij</sub>中的隶属度,设w<sub>i</sub>(ξ(k))≥0,for i=1,2,…,r,r为规则数,ξ(k)表示包含ξ<sub>1</sub>(k),...,ξ<sub>g</sub>(k)的向量,μ<sub>i</sub>(ξ(k))≥0,<img file="FDA0000608375790000021.GIF" wi="341" he="116" />为了便于记录我们令μ<sub>i</sub>=μ<sub>i</sub>(ξ(k)),<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>A</mi><mrow><mo>(</mo><mi>&mu;</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>r</mi></munderover><msub><mi>&mu;</mi><mi>i</mi></msub><msub><mi>A</mi><mi>i</mi></msub><mo>,</mo><mi>B</mi><mrow><mo>(</mo><mi>&mu;</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>r</mi></munderover><msub><mi>&mu;</mi><mi>i</mi></msub><msub><mi>B</mi><mi>i</mi></msub><mo>,</mo><mi>C</mi><mrow><mo>(</mo><mi>&mu;</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>r</mi></munderover><msub><mi>&mu;</mi><mi>i</mi></msub><msub><mi>C</mi><mi>i</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000608375790000022.GIF" wi="1244" he="157" /></maths>步骤3、随机模糊状态反馈鲁棒控制器设计在实际网络化控制系统中,因网络通信的不稳定,常常在传感器‑控制器通道和控制器‑执行器通道中丢失数据,导致传感器的输出不等于控制器输入与控制器输出不等于执行器输入的情况。为了描述上述数据丢包问题,本发明采用x<sub>c</sub>(k)∈R<sup>(n+m)</sup>描述控制器输入,<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>x</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>x</mi><mi>cs</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>x</mi><mi>cf</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000608375790000023.GIF" wi="418" he="133" /></maths>x<sub>cs</sub>(k)∈R<sup>n</sup>为慢状态变量,x<sub>cf</sub>(k)∈R<sup>m</sup>为快状态变量;u<sub>c</sub>(k)∈R<sup>q</sup>表示控制器输出,并将数据丢包描述为服从Bernoulli的取值为0和1随机切换序列,即x<sub>c</sub>(k)=α<sub>h</sub>x(k)          (7)u(k)=β<sub>h</sub>u<sub>c</sub>(k)         (8)其中,随机变量α<sub>h</sub>∈R和β<sub>h</sub>∈R分别为描述传感器‑控制器通道和控制器‑执行器通道数据丢包现象的相互独立的Bernoulli随机过程;α<sub>h</sub>=0表示传感器‑控制器通道通信失败数据丢包,α<sub>h</sub>=1表示传感器‑控制器通道通信成功;β<sub>h</sub>=0表示控制器‑执行器通道通信失败数据丢包,β<sub>h</sub>=1表示控制器‑执行器通道通信成功;为描述整体闭环系统数据丢包情况,定义另一Bernoulli过程δ<sub>h</sub>=α<sub>h</sub>β<sub>h</sub>,δ<sub>h</sub>=0表示闭环系统通信失败数据丢包,δ<sub>h</sub>=1表示闭环系统通信成功;随机变量δ<sub>h</sub>的分布律为P{δ<sub>h</sub>=h}=p<sup>h</sup>(1‑p)<sup>1‑h</sup>,h=0,1(0<p<1)     (9)其中,p表示闭环系统通信成功的概率,h表示δ<sub>h</sub>的取值;随机变量δ<sub>h</sub>的均值为:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mfenced open='' close=''><mtable><mtr><mtd><mi>E</mi><mrow><mo>(</mo><msub><mi>&delta;</mi><mi>h</mi></msub><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>h</mi><mo>=</mo><mn>0</mn></mrow><mn>1</mn></munderover><mi>hP</mi><mo>{</mo><msub><mi>&delta;</mi><mi>h</mi></msub><mo>=</mo><mi>h</mi><mo>}</mo></mtd></mtr><mtr><mtd><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>h</mi><mo>=</mo><mn>0</mn></mrow><mn>1</mn></munderover><msup><mi>hp</mi><mi>h</mi></msup><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>p</mi><mo>)</mo></mrow><mrow><mn>1</mn><mo>-</mo><mi>h</mi></mrow></msup></mtd></mtr><mtr><mtd><mo>=</mo><mi>p</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000608375790000024.GIF" wi="1067" he="389" /></maths>基于离散时间模糊奇异摄动模型(5),对被控对象NTTSSs设计随机模糊状态反馈控制器:控制器规则i:如果ξ<sub>1</sub>(k)是φ<sub>i1</sub>,...,ξ<sub>g</sub>(k)是φ<sub>ig</sub>,那么u<sub>c</sub>(k)=F<sub>i</sub>x<sub>c</sub>(k)         (11)其中,F<sub>i</sub>为控制器增益;采用标准的模糊推理方法‑即单点模糊化,乘积推理,和加权平均清晰化,得全局模糊控制器为u<sub>c</sub>(k)=F(μ)x<sub>c</sub>(k)         (12)其中,<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mi>F</mi><mrow><mo>(</mo><mi>&mu;</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>r</mi></munderover><msub><mi>&mu;</mi><mi>i</mi></msub><msub><mi>F</mi><mi>i</mi></msub><mo>;</mo></mrow>]]></math><img file="FDA0000608375790000031.GIF" wi="349" he="150" /></maths>将式(7)、(12)代入式(8)得被控系统控制输入u(k)=δ<sub>h</sub>F(μ)x(k)       (13)步骤4、推导网络化NTTSSs的闭环系统模型将式(13)代入式(5)得网络化NTTSSs的闭环系统模型:x(k+1)=E<sub>ε</sub>[A(μ)+δ<sub>h</sub>B(μ)F(μ)]x(k)y(k)=C(μ)x(k)          (14)步骤5、求解控制器增益采用谱范数方法和线性矩阵不等式方法,推导出鲁棒随机模糊状态反馈控制器存在的充分条件,下面是求解控制器增益的线性矩阵不等式:<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><mfenced open='[' close=']'><mtable><mtr><mtd><mo>-</mo><msup><mi>&gamma;</mi><mn>2</mn></msup><mi>Q</mi></mtd><mtd><mo>*</mo></mtd></mtr><mtr><mtd><msub><mi>A</mi><mi>i</mi></msub><mi>Q</mi><mo>+</mo><mi>p</mi><msub><mi>B</mi><mi>i</mi></msub><msub><mi>W</mi><mi>i</mi></msub></mtd><mtd><mo>-</mo><msubsup><mi>E</mi><mi>&epsiv;</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mi>Q</mi><msubsup><mi>E</mi><mi>&epsiv;</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup></mtd></mtr></mtable></mfenced><mo>&lt;</mo><mn>0</mn><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mi>r</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>15</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000608375790000032.GIF" wi="1462" he="128" /></maths><maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><mfenced open='[' close=']'><mtable><mtr><mtd><mo>-</mo><mn>2</mn><msup><mi>&gamma;</mi><mn>2</mn></msup><mi>Q</mi></mtd><mtd><mo>*</mo></mtd></mtr><mtr><mtd><mrow><mo>(</mo><msub><mi>A</mi><mi>i</mi></msub><mo>+</mo><msub><mi>A</mi><mi>j</mi></msub><mo>)</mo></mrow><mi>Q</mi><mo>+</mo><mi>p</mi><msub><mi>B</mi><mi>i</mi></msub><msub><mi>W</mi><mi>j</mi></msub><mo>+</mo><mi>p</mi><msub><mi>B</mi><mi>j</mi></msub><msub><mi>W</mi><mi>i</mi></msub></mtd><mtd><mo>-</mo><mn>2</mn><msubsup><mi>E</mi><mi>&epsiv;</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mi>Q</mi><msubsup><mi>E</mi><mi>&epsiv;</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup></mtd></mtr></mtable></mfenced><mo>&lt;</mo><mn>0,1</mn><mo>&le;</mo><mi>i</mi><mo>&lt;</mo><mi>j</mi><mo>&le;</mo><mi>r</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>16</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000608375790000033.GIF" wi="1545" he="134" /></maths>Q>0              (17)其中,Q∈R<sup>(n+m)×(n+m)</sup>为对称正定矩阵,W<sub>i</sub>∈R<sup>q×(n+m)</sup>为矩阵,q为控制输入u(k)的维数,0<γ≤1,p为网络化系统收到有效数据包的概率;控制器增益:F<sub>i</sub>=W<sub>i</sub>*Q<sup>‑1</sup>,i=1,2,…,r          (18)步骤6、将所得控制器Matlab代码传化为C语言代码,植入控制器;传感器和执行器均置于被控对象端,并通过采用相同的物理时钟,实现两者之间的时钟同步;执行器由缓冲区和零阶保持器组成,传感器采用时间驱动方式,控制器采用事件驱动方式,当采样数据包到达控制器节点时,控制器立刻采用该采样数据包进行控制信号计算,并将控制信号与采样数据包携带的时间戳封装成控制数据包,通过网络传送给执行器;当控制数据包到达执行器节点后,执行器将其携带的时间戳与缓冲区内控制信号的时间戳进行比较,并判断新到达的控制数据包是否“新”;“是”则将新到达的控制信号及其时间戳保存在缓冲区中,“否”则丢弃该控制数据包,零阶保持器采用时间驱动方式,即零阶保持器按照固定的采样周期,从缓冲区读取控制信号,并生成控制输入作用于被控对象,从而实现网络化NTTSSs的高精度控制。
地址 100071 北京市丰台区西四环南路72号