发明名称 不确定性时延双时标系统模糊时延状态反馈控制方法
摘要 一种不确定性时延双时标系统模糊时延状态反馈控制方法,复杂系统高精度控制技术领域。该方法基于不确定性标准离散模糊奇异摄动模型,设计模糊时延状态反馈控制器,实现复杂UTDNTTSSs的高精度控制。融合模糊逻辑、奇异摄动技术、时延及不确定性系统理论,建立被控UTDNTTSSs的不确定性标准离散模糊奇异摄动模型,结合模糊逻辑与时延系统理论,设计模糊时延状态反馈控制器,采用谱范数与LMIs方法,推导出控制器存在充分条件,通过CE150直升机姿态控制系统高精度控制仿真,说明了该方法的具体实施过程与有效性。优点在于,解决现有UTDNTTSSs建模与控制方法,无法消除时延、快模态及参数摄动引起的稳态误差问题,大幅提高UTDNTTSSs的控制性能。
申请公布号 CN104460322A 申请公布日期 2015.03.25
申请号 CN201410815382.3 申请日期 2014.12.23
申请人 冶金自动化研究设计院 发明人 陈金香
分类号 G05B13/04(2006.01)I;G06F19/00(2011.01)I 主分类号 G05B13/04(2006.01)I
代理机构 北京华谊知识产权代理有限公司 11207 代理人 刘月娥
主权项 一种时延非线性双时标系统的时延状态反馈控制方法,其特征在于:步骤一、根据被控UTDNTTSSs的动力学方程,建立其不确定性连续时延模糊奇异摄动模型规则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><mrow><mo>(</mo><msub><mi>A</mi><mi>ci</mi></msub><mo>+</mo><mi>&Delta;</mi><msub><mi>A</mi><mi>ci</mi></msub><mo>)</mo></mrow><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><msub><mi>A</mi><mi>cdi</mi></msub><mo>+</mo><mi>&Delta;</mi><msub><mi>A</mi><mi>cdi</mi></msub><mo>)</mo></mrow><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>&tau;</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="FDA0000641401920000011.GIF" wi="1245" he="67" /></maths>y(t)=Cx(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="FDA0000641401920000012.GIF" wi="867" he="135" /></maths>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>为系统输出,w(t)∈R<sup>q</sup>为外扰,φ<sub>i1</sub>,...,φ<sub>ig</sub>(i=1,2,...,r)均为模糊集合,ξ<sub>1</sub>(t),...,ξ<sub>g</sub>(t)为可测量的系统变量,A<sub>ci</sub>,A<sub>cdi</sub>,B<sub>ci</sub>,C为合适维数矩阵,ΔA<sub>ci</sub>andΔA<sub>cdi</sub>为合适维数的不确定性矩阵,ε是奇异摄动参数,τ(0<τ<τ<sub>m</sub>)为时延常数,τ<sub>m</sub>为上确界;步骤二、建立被控UTDNTTSSs的不确定性时延标准离散模糊奇异摄动模型控制系统中的传感器和执行器均采用时间驱动方式,且二者采用相同的采样时间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>+ΔA<sub>i</sub>)x(k)+E<sub>ε</sub>(A<sub>di</sub>+ΔA<sub>di</sub>)x(k‑τ)+E<sub>ε</sub>B<sub>i</sub>u(k)y(k)=Cx(k)                  (2)其中,ΔA<sub>i</sub>,ΔA<sub>di</sub>为适当维数的不确定性矩阵,A<sub>i</sub>,A<sub>di</sub>,B<sub>i</sub>为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>A</mi><mi>i</mi></msub><mo>=</mo><msubsup><mi>E</mi><mi>&epsiv;</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msup><mi>e</mi><mrow><msubsup><mi>E</mi><mi>&epsiv;</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msub><mi>A</mi><mi>ci</mi></msub><msub><mi>T</mi><mi>s</mi></msub></mrow></msup><mo>,</mo><msub><mi>A</mi><mi>di</mi></msub><mo>=</mo><msubsup><mi>E</mi><mi>&epsiv;</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msup><mi>e</mi><mrow><msubsup><mi>E</mi><mi>&epsiv;</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msub><mi>A</mi><mi>cdi</mi></msub><msub><mi>T</mi><mi>s</mi></msub></mrow></msup><mo>,</mo><msub><mi>B</mi><mi>i</mi></msub><mo>=</mo><msubsup><mi>E</mi><mi>&epsiv;</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msubsup><mo>&Integral;</mo><mn>0</mn><mi>h</mi></msubsup><msubsup><mi>E</mi><mi>&epsiv;</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msup><mi>e</mi><mrow><msubsup><mi>E</mi><mi>&epsiv;</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msub><mi>A</mi><mi>ci</mi></msub><mi>&tau;</mi></mrow></msup><mi>d&tau;</mi><msub><mi>B</mi><mi>ci</mi></msub></mrow>]]></math><img file="FDA0000641401920000013.GIF" wi="1255" he="109" /></maths>给定[x(k);u(k)],应用标准模糊推理方法,得到全局不确定性时延标准离散模糊奇异摄动模型:x(k+1)=E<sub>ε</sub>(A(μ)+ΔA(μ))x(k)+E<sub>ε</sub>(A<sub>d</sub>(μ)+ΔA<sub>d</sub>(μ))x(k‑τ)+E<sub>ε</sub>B(μ)u(k)y(k)=Cx(k)                   (3)其中,隶属度函数<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="FDA0000641401920000014.GIF" wi="1001" he="196" /></maths>φ<sub>ij</sub>(ξ<sub>j</sub>(k))为ξ<sub>j</sub>(k)在φ<sub>ij</sub>中的隶属度,ξ(k)表示包含ξ<sub>1</sub>(k),...,ξ<sub>g</sub>(k)的向量,设w<sub>i</sub>(ξ(k))≥0,for i=1,2,…,r,r为规则数,μ<sub>i</sub>(ξ(k))≥0,<img file="FDA0000641401920000015.GIF" wi="345" he="124" />为了便于记录,令μ<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><msub><mi>A</mi><mi>d</mi></msub><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>di</mi></msub></mrow>]]></math><img file="FDA0000641401920000021.GIF" wi="1061" he="151" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>&Delta;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><mi>&Delta;</mi><msub><mi>A</mi><mi>i</mi></msub><mo>,</mo><mi>&Delta;</mi><msub><mi>A</mi><mi>d</mi></msub><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><mi>&Delta;</mi><msub><mi>A</mi><mi>di</mi></msub><mo>,</mo></mrow>]]></math><img file="FDA0000641401920000022.GIF" wi="1062" he="149" /></maths>步骤三、时延模糊状态反馈鲁棒控制器设计基于模型(3),对被控对象UTDNTTSSs,设计时延模糊状态反馈鲁棒控制器;控制器规则i:如果ξ<sub>1</sub>(k)是φ<sub>i1</sub>,...,ξ<sub>g</sub>(k)是φ<sub>ig</sub>,那么u(k)=F<sub>i</sub>x(k)+F<sub>di</sub>x(k‑τ)          (4)其中,F<sub>i</sub>、F<sub>di</sub>均为控制器增益;采用标准的模糊推理方法‑即单点模糊化,乘积推理,和加权平均清晰化,得全局模糊控制器为u(k)=F(μ)x(k)+F<sub>d</sub>(μ)x(k‑τ)            (5)其中,<maths num="0007" id="cmaths0007"><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><msub><mi>F</mi><mi>d</mi></msub><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>di</mi></msub><mo>;</mo></mrow>]]></math><img file="FDA0000641401920000023.GIF" wi="752" he="154" /></maths>步骤四、推导出闭环系统模型将控制率(5)应用于被控系统模型(3),获得闭环系统模型:x(k+1)=E<sub>ε</sub>(A(μ)+B(μ)F(μ)+ΔA(μ))x(k)+E<sub>ε</sub>(A<sub>d</sub>(μ)+B(μ)F<sub>d</sub>(μ)+ΔA<sub>d</sub>(μ))x(k‑τ)                                                   (6)步骤五、采用谱范数方法和线性矩阵不等式方法,推导出时延模糊状态反馈鲁棒控制器存在的充分条件,推导过程不要求知道系统不确定性参数的上确界,下面是求解控制器增益的线性矩阵不等式组:Γ<sub>ii</sub><0 i=1,2,...,r.                          (7)Γ<sub>ij</sub>+Γ<sub>ji</sub><0 i<j<r(i=1,2,...,r;j=2,3,...,r)                 (8)P>0                          (9)Y>0                          (10)其中,γ(0<γ≤1)为给定常数,<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mi>P</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>P</mi><mn>11</mn></msub></mtd><mtd><mo>*</mo></mtd></mtr><mtr><mtd><msub><mi>P</mi><mn>21</mn></msub></mtd><mtd><msub><mi>P</mi><mn>22</mn></msub></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000641401920000024.GIF" wi="375" he="131" /></maths>(P<sub>11</sub>∈R<sup>(n+m)τ×(n+m)τ</sup>,P<sub>22</sub>∈R<sup>(n+m)×(n+m)</sup>为对称正定矩阵),<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><mi>Y</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>Y</mi><mn>11</mn></msub></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msub><mi>Y</mi><mn>22</mn></msub></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000641401920000025.GIF" wi="367" he="133" /></maths>(Y<sub>11</sub>,Y<sub>22</sub>∈R<sup>(n+m)×(n+m)</sup>为对称正定矩阵且<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><msub><mi>Y</mi><mn>11</mn></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>Y</mi><mn>111</mn></msub></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msub><mi>Y</mi><mn>112</mn></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000641401920000026.GIF" wi="442" he="135" /></maths>Y<sub>111</sub>∈R<sup>(n+m)×(n+m)(τ‑1)</sup>,Y<sub>112</sub>∈R<sup>(n+m)(τ‑1)×(n+m)</sup>),矩阵V<sub>i1</sub>∈R<sup>q×(n+m)(τ‑1)</sup>,W<sub>i</sub>∈R<sup>q×(n+m)</sup>,<img file="FDA0000641401920000027.GIF" wi="1158" he="242" /><img file="FDA0000641401920000031.GIF" wi="1166" he="244" /><maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><mi>S</mi><mo>=</mo><msub><mi>P</mi><mn>11</mn></msub><mo>-</mo><mi>&gamma;</mi><msub><mi>Y</mi><mn>11</mn></msub><mo>-</mo><mi>&gamma;</mi><msubsup><mi>Y</mi><mn>11</mn><mi>T</mi></msubsup><mo>,</mo></mrow>]]></math><img file="FDA0000641401920000038.GIF" wi="485" he="77" /></maths><maths num="0012" id="cmaths0012"><math><![CDATA[<mrow><mi>J</mi><mo>=</mo><msub><mi>P</mi><mn>22</mn></msub><mo>-</mo><mi>&gamma;</mi><msub><mi>Y</mi><mn>22</mn></msub><mo>-</mo><mi>&gamma;</mi><msubsup><mi>Y</mi><mn>22</mn><mi>T</mi></msubsup><mo>,</mo></mrow>]]></math><img file="FDA0000641401920000033.GIF" wi="485" he="80" /></maths>φ<sub>i</sub>=[A<sub>di</sub> 0<sub>(n+m)×(n+m)(τ‑1)</sub>],χ<sub>i</sub>=[B<sub>i</sub> 0<sub>(n+m)×(n+m)(τ‑1)</sub>],<maths num="0013" id="cmaths0013"><math><![CDATA[<mrow><msub><mi>V</mi><mi>i</mi></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>V</mi><mrow><mi>i</mi><mn>1</mn></mrow></msub></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msub><mn>0</mn><mrow><mrow><mo>(</mo><mi>n</mi><mo>+</mo><mi>m</mi><mo>)</mo></mrow><mrow><mo>(</mo><mi>&tau;</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>&times;</mo><mrow><mo>(</mo><mi>n</mi><mo>+</mo><mi>m</mi><mo>)</mo></mrow></mrow></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000641401920000034.GIF" wi="677" he="136" /></maths>Ω<sub>ij</sub>=A<sub>i</sub>Y<sub>22</sub>+B<sub>i</sub>W<sub>j</sub>,时延τ=1时Υ=0<sub>(n+m)×(n+m)</sub>,∏=I<sub>(n+m)×(n+m)</sub>,时延τ≥2时<img file="FDA0000641401920000035.GIF" wi="941" he="355" /><maths num="0014" id="cmaths0014"><math><![CDATA[<mrow><mi>&Pi;</mi><mo>=</mo><msubsup><mfenced open='[' close=']'><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mn>0</mn></mtd><mtd><mi>I</mi></mtd></mtr></mtable></mfenced><mrow><mrow><mo>(</mo><mi>n</mi><mo>+</mo><mi>m</mi><mo>)</mo></mrow><mi>&tau;</mi><mo>&times;</mo><mrow><mo>(</mo><mi>n</mi><mo>+</mo><mi>m</mi><mo>)</mo></mrow></mrow><mi>T</mi></msubsup><mo>,</mo></mrow>]]></math><img file="FDA0000641401920000036.GIF" wi="797" he="99" /></maths>控制器增益:<maths num="0015" id="cmaths0015"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><msub><mi>F</mi><mi>i</mi></msub><mo>=</mo><msub><mi>W</mi><mi>i</mi></msub><mo>*</mo><msubsup><mi>Y</mi><mn>22</mn><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo>,</mo><msub><mi>F</mi><mi>di</mi></msub><mo>=</mo><msub><mi>V</mi><mrow><mi>i</mi><mn>1</mn></mrow></msub><mo>*</mo><msubsup><mi>Y</mi><mn>111</mn><mrow><mo>-</mo><mn>1</mn></mrow></msubsup></mtd><mtd><mi>for</mi></mtd><mtd><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>r</mi><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000641401920000037.GIF" wi="1464" he="85" /></maths>步骤六、将所得控制器Matlab代码传化为C语言代码,植入控制器,控制器采用事件驱动方式,当采样数据到达控制器时,控制器立刻进行计算,并将控制信号传给执行器,执行器按照固定的采样周期读取控制信号,生成控制输入,作用于被控UTDNTTSSs,实现UTDNTTSSs的高精度控制。
地址 100071 北京市丰台区西四环南路72号