发明名称 一种复合分层抗干扰控制器
摘要 一种复合分层抗干扰控制器,涉及一类多源干扰控制系统的复合分层抗干扰控制器设计;该发明针对一类含有外部模型描述干扰、随机干扰和能量有界干扰的多源干扰控制系统。首先,构造降阶干扰观测器估计并抵消多源干扰控制系统中的外部模型描述干扰;其次,构造具有H<sub>∞</sub>和保成本性能指标的控制器,其中H<sub>∞</sub>性能指标抑制系统的能量有界干扰,保成本性能指标抑制系统的随机干扰,并能优化方差的上界;基于降阶干扰观测器、鲁棒H<sub>∞</sub>和保成本状态反馈控制器,构造复合分层抗干扰控制器;最后,基于凸优化算法求解复合分层抗干扰控制器增益阵列。本发明具有抗干扰性强、便于设计等优点,可用于航空航天器导航与控制系统等。
申请公布号 CN101571704A 申请公布日期 2009.11.04
申请号 CN200910086897.3 申请日期 2009.06.18
申请人 北京航空航天大学 发明人 郭雷;文新宇;曹松银;刘华;全伟;杨照华
分类号 G05B13/00(2006.01)I 主分类号 G05B13/00(2006.01)I
代理机构 北京科迪生专利代理有限责任公司 代理人 成金玉;卢 纪
主权项 1、一种复合分层抗干扰控制器,其特征在于包括以下步骤:首先,构造降阶干扰观测器估计并抵消多源干扰控制系统中的外部模型描述干扰;其次,构造具有H<sub>∞</sub>和保成本性能指标的控制器,抑制多源干扰控制系统中随机干扰和能量有界干扰;基于降阶干扰观测器、鲁棒H<sub>∞</sub>和保成本状态反馈控制器,构造复合分层抗干扰控制器;最后,基于凸优化算法求解复合分层抗干扰控制器增益阵列。具体步骤如下:<maths num="0001"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><mover><mi>x</mi><mo>.</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>Ax</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mi>Ff</mi><mrow><mo>(</mo><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>+</mo><msub><mi>B</mi><mn>1</mn></msub><mo>[</mo><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>d</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>]</mo><mo>+</mo><msub><mi>B</mi><mn>2</mn></msub><msub><mi>d</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>d</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>y</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>Cx</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>D</mi><mn>1</mn></msub><msub><mi>d</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>D</mi><mn>2</mn></msub><msub><mi>d</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>d</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced></math>]]></maths>其中,多源干扰控制系统状态变量x(t)∈R<sup>n</sup>,R<sup>n</sup>表示n维实向量空间,系统阵A∈R<sup>n×n</sup>,R<sup>n×n</sup>表示n×n维实矩阵空间,f(x(t))为多源干扰控制系统的非线性项,其增益阵为F,<maths num="0002"><![CDATA[<math><mrow><msub><mi>d</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&Element;</mo><msup><mi>R</mi><msub><mi>q</mi><mn>1</mn></msub></msup></mrow></math>]]></maths>为外部模型描述干扰,其增益阵为<maths num="0003"><![CDATA[<math><mrow><msub><mi>B</mi><mn>1</mn></msub><mo>&Element;</mo><msup><mi>R</mi><mrow><mi>n</mi><mo>&times;</mo><msub><mi>q</mi><mn>1</mn></msub></mrow></msup></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msub><mi>d</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&Element;</mo><msup><mi>R</mi><msub><mi>q</mi><mn>2</mn></msub></msup></mrow></math>]]></maths>为能量有界干扰,其增益阵为<maths num="0005"><![CDATA[<math><mrow><msub><mi>B</mi><mn>2</mn></msub><mo>&Element;</mo><msup><mi>R</mi><mrow><mi>n</mi><mo>&times;</mo><msub><mi>q</mi><mn>2</mn></msub></mrow></msup><mo>,</mo></mrow></math>]]></maths>d<sub>3</sub>(t)∈R<sup>n</sup>为随机干扰,y(t)∈R<sup>m</sup>为系统输出变量,C∈R<sup>m×n</sup>为系统输出阵,输出系统多源干扰增益阵分别为D<sub>1</sub>,D<sub>2</sub>;(1)构造降阶干扰观测器为:<maths num="0006"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><mover><mi>r</mi><mo>.</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><mi>W</mi><mo>+</mo><msub><mi>KB</mi><mn>1</mn></msub><mi>V</mi><mo>)</mo></mrow><mi>r</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mo>[</mo><mi>KA</mi><mo>-</mo><mrow><mo>(</mo><mi>W</mi><mo>+</mo><msub><mi>KB</mi><mn>1</mn></msub><mi>V</mi><mo>)</mo></mrow><mi>K</mi><mo>]</mo><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>KB</mi><mn>1</mn></msub><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mi>KFf</mi><mrow><mo>(</mo><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mover><mi>w</mi><mo>^</mo></mover><mo>=</mo><mi>r</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mi>Kx</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><msub><mover><mi>d</mi><mo>^</mo></mover><mn>1</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>V</mi><mover><mi>w</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced></math>]]></maths>其中,r(t)为参考状态变量,<img file="A2009100868970002C7.GIF" wi="97" he="60" />为外部干扰模型∑<sub>2</sub>状态变量w(t)∈R<sup>p</sup>的估计值,R<sup>p</sup>表示p维实向量空间,x(t)∈R<sup>n</sup>为多源干扰控制系统∑<sub>1</sub>的可测状态变量,A为系统阵,u(t)∈R<sup>q</sup>为控制输入,控制输入通道增益阵为<maths num="0007"><![CDATA[<math><mrow><msub><mi>B</mi><mn>1</mn></msub><mo>&Element;</mo><msup><mi>R</mi><mrow><mi>n</mi><mo>&times;</mo><msub><mi>q</mi><mn>1</mn></msub></mrow></msup><mo>,</mo></mrow></math>]]></maths><img file="A2009100868970002C9.GIF" wi="87" he="41" />表示n×q<sub>1</sub>维实矩阵空间,f(x(t))为多源干扰控制系统∑<sub>1</sub>的非线性项,F为其增益阵,K∈R<sup>p×m</sup>为待定观测器增益阵,W∈R<sup>p×p</sup>、<maths num="0008"><![CDATA[<math><mrow><mi>V</mi><mo>&Element;</mo><msup><mi>R</mi><mrow><msub><mi>q</mi><mn>1</mn></msub><mo>&times;</mo><mi>p</mi></mrow></msup></mrow></math>]]></maths>分别表示外部干扰模型∑<sub>2</sub>的系统矩阵、输出矩阵,<img file="A2009100868970002C11.GIF" wi="106" he="76" />为干扰观测器输出变量;m、n、p、q<sub>1</sub>均为自然数;(2)构造鲁棒H<sub>∞</sub>和保成本状态反馈控制器为:u<sub>f</sub>(t)=Lx(t)其中,u<sub>f</sub>(t)为状态反馈控制器,L为待定状态反馈控制器增益阵;(3)基于降阶干扰观测器、鲁棒H<sub>∞</sub>和保成本状态反馈控制器,构造复合分层抗干扰控制器为:<maths num="0009"><![CDATA[<math><mrow><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><msub><mover><mi>d</mi><mo>^</mo></mover><mn>1</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>u</mi><mi>f</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></math>]]></maths>其中,<img file="A2009100868970003C2.GIF" wi="92" he="67" />为外部模型描述干扰d<sub>1</sub>(t)的估计值,u<sub>f</sub>(t)为状态反馈控制器;(4)基于凸优化算法得降阶干扰观测器增益阵为<maths num="0010"><![CDATA[<math><mrow><mi>K</mi><mo>=</mo><msubsup><mi>P</mi><mn>2</mn><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msub><mi>R</mi><mn>2</mn></msub><mo>,</mo></mrow></math>]]></maths>状态反馈控制器增益阵为<maths num="0011"><![CDATA[<math><mrow><mi>L</mi><mo>=</mo><msub><mi>R</mi><mn>1</mn></msub><msubsup><mi>P</mi><mn>1</mn><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo>,</mo></mrow></math>]]></maths>其中P<sub>1</sub>、R<sub>1</sub>、P<sub>2</sub>、R<sub>2</sub>由以下凸优化问题求得:<maths num="0012"><![CDATA[<math><mrow><mi>min</mi><mfenced open='(' close=')' separators=' '><mtable><mtr><mtd></mtd></mtr></mtable><mfenced open='(' close=')'><mtable><mtr><mtd><msup><mi>x</mi><mi>T</mi></msup><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></mtd><mtd><msubsup><mi>e</mi><mi>w</mi><mi>T</mi></msubsup><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mfenced open='(' close=')'><mtable><mtr><mtd><msub><mi>P</mi><mn>1</mn></msub></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msub><mi>P</mi><mn>2</mn></msub></mtd></mtr></mtable></mfenced><msup><mfenced open='(' close=')'><mtable><mtr><mtd><msup><mi>x</mi><mi>T</mi></msup><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></mtd><mtd><msubsup><mi>e</mi><mi>w</mi><mi>T</mi></msubsup><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mi>T</mi></msup></mfenced></mrow></math>]]></maths><maths num="0013"><![CDATA[<math><mrow><mi>&Phi;</mi><mo>=</mo><mtable><mtr><mtd><mrow><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>&Phi;</mi><mn>11</mn></msub></mtd><mtd><msub><mi>B</mi><mn>1</mn></msub><mi>V</mi></mtd><mtd><msub><mi>B</mi><mn>2</mn></msub></mtd><mtd><msub><mi>P</mi><mn>1</mn></msub><msubsup><mi>C</mi><mn>1</mn><mi>T</mi></msubsup></mtd><mtd><mi>F</mi></mtd><mtd><msub><mi>P</mi><mn>1</mn></msub><msubsup><mi>C</mi><mn>3</mn><mi>T</mi></msubsup></mtd><mtd><mi>&lambda;</mi><msub><mi>P</mi><mn>1</mn></msub><msup><mi>U</mi><mi>T</mi></msup></mtd></mtr><mtr><mtd><mo>*</mo></mtd><mtd><msub><mi>&Phi;</mi><mn>22</mn></msub></mtd><mtd><msub><mi>R</mi><mn>2</mn></msub><msub><mi>B</mi><mn>2</mn></msub></mtd><mtd><msubsup><mi>C</mi><mn>2</mn><mi>T</mi></msubsup></mtd><mtd><mn>0</mn></mtd><mtd><msubsup><mi>C</mi><mn>4</mn><mi>T</mi></msubsup></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd><mtd><mo>-</mo><msup><mi>&gamma;</mi><mn>2</mn></msup><mi>I</mi></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd><mtd><mo>-</mo><mi>I</mi></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd><mtd><mo>-</mo><msup><mi>&lambda;</mi><mn>2</mn></msup><mi>I</mi></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd><mtd><mo>-</mo><mi>I</mi></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd><mtd><mo>*</mo></mtd><mtd><mo>-</mo><mi>I</mi></mtd></mtr></mtable></mfenced><mo>&lt;</mo><mn>0</mn></mrow></mtd></mtr></mtable></mrow></math>]]></maths>其中:x(0)、e<sub>w</sub>(0)为给定初始值,[C<sub>1</sub> C<sub>2</sub>]∈R<sup>(n+p)×(n+p)</sup>、[C<sub>3</sub> C<sub>4</sub>]∈R<sup>(n+p)×(n+p)</sup>为估计误差闭环系统∑<sub>3</sub>可调输出阵,R<sup>(n+p)×(n+p)</sup>表示(n+p)×(n+p)维实矩阵空间,λ为非线性权重参数,γ为干扰抑制度;Φ<sub>11</sub>=(AP<sub>1</sub>+B<sub>1</sub>R<sub>1</sub>)+(AP<sub>1</sub>+B<sub>1</sub>R<sub>1</sub>)<sup>T</sup>,U为∑<sub>1</sub>中非线性项f(x(t))的Lipschitz参数阵,Φ<sub>22</sub>=(P<sub>2</sub>W+R<sub>2</sub>B<sub>1</sub>V)+(P<sub>2</sub>W+R<sub>2</sub>B<sub>1</sub>V)<sup>T</sup>,<maths num="0014"><![CDATA[<math><mrow><msub><mi>B</mi><mn>2</mn></msub><mo>&Element;</mo><msup><mi>R</mi><mrow><mi>n</mi><mo>&times;</mo><msub><mi>q</mi><mn>2</mn></msub></mrow></msup></mrow></math>]]></maths>为能量有界干扰增益阵;符号*表示对称矩阵中相应部分的对称块,q<sub>2</sub>为自然数。
地址 100191北京市海淀区学院路37号
您可能感兴趣的专利