发明名称 一种针对多源干扰系统的复合分层自适应滤波器
摘要 本发明涉及一种针对多源干扰系统的复合分层自适应滤波器,所述多源干扰,包括外系统输出变量、参数不确定性和范数有界噪声;首先,针对由外系统输出变量所描述的干扰,利用已知信息建立干扰模型,进而设计干扰观测器估计干扰并抵消;其次,针对参数不确定性,设计自适应更新律,对不确定参数进行估计,进而构造滤波器输入项进行补偿;再次,通过H<sub>∞</sub>性能优化,进一步抑制以上两步中产生的估计误差和模型中的范数有界噪声;最后,结合干扰观测器、自适应输入补偿项和H<sub>∞</sub>性能优化结果,构造复合分层滤波器;本发明具有精细抗干扰的特点,相对于传统的卡尔曼滤波器和H<sub>∞</sub>滤波器精度更高,适用面更广,可应用于车辆、船舶以及飞行器的导航系统。
申请公布号 CN105865441A 申请公布日期 2016.08.17
申请号 CN201610196187.6 申请日期 2016.03.31
申请人 北京航空航天大学 发明人 郭雷;李文硕;乔建忠;张培喜;王春
分类号 G01C21/00(2006.01)I 主分类号 G01C21/00(2006.01)I
代理机构 北京科迪生专利代理有限责任公司 11251 代理人 杨学明;顾炜
主权项 一种针对多源干扰系统的复合分层自适应滤波器,其特征在于包括以下步骤:首先,建立多源干扰系统的状态空间模型;其次,针对由外系统输出变量所描述的干扰,利用已知信息建立外系统模型,进而构造干扰观测器估计并抵消;再次,针对参数不确定性,设计自适应更新律,对不确定参数进行估计,进而构造滤波器输入项进行补偿;下一步,通过H<sub>∞</sub>性能优化,进一步抑制估计误差和范数有界噪声对参考性能的影响;最后,将干扰观测器、自适应输入补偿项和H<sub>∞</sub>滤波器进行复合,得到复合分层自适应滤波器;具体步骤如下:第一步,建立多源干扰系统的状态空间模型针对含有外部可建模干扰、参数不确定性干扰和范数有界噪声的多源干扰系统,建立系统状态空间模型如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>&Sigma;</mi><mn>1</mn></msub><mo>:</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mo>(</mo><mi>t</mi><mo>)</mo><mo>=</mo><mi>A</mi><mi>x</mi><mo>(</mo><mi>t</mi><mo>)</mo><mo>+</mo><mi>F</mi><mi>f</mi><mo>(</mo><mi>x</mi><mo>(</mo><mi>t</mi><mo>)</mo><mo>)</mo><mo>+</mo><msub><mi>B</mi><mn>1</mn></msub><msub><mi>d</mi><mn>1</mn></msub><mo>(</mo><mi>t</mi><mo>)</mo><mo>+</mo><msub><mi>B</mi><mn>2</mn></msub><msub><mi>d</mi><mn>2</mn></msub><mo>(</mo><mi>t</mi><mo>)</mo><mo>+</mo><mi>&omega;</mi><mo>(</mo><mi>t</mi><mo>)</mo></mtd></mtr><mtr><mtd><mi>y</mi><mo>(</mo><mi>t</mi><mo>)</mo><mo>=</mo><mi>C</mi><mi>x</mi><mo>(</mo><mi>t</mi><mo>)</mo><mo>+</mo><mi>&upsi;</mi><mo>(</mo><mi>t</mi><mo>)</mo></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000954861480000011.GIF" wi="1246" he="167" /></maths>其中,x(t)为多源干扰系统状态变量,y(t)为多源干扰系统输出变量,d<sub>1</sub>(t)为外部可建模干扰,d<sub>2</sub>(t)为参数不确定性干扰,ω(t)、υ(t)分别为系统方程和量测方程中的范数有界噪声,A,B<sub>1</sub>,B<sub>2</sub>和C为参数已知的矩阵,非线性项f(x(t))表示系统建模时的线性化误差项或神经网络逼近误差,满足Lipschitz条件,即存在已知Lipschitz参数阵U∈R<sup>n×n</sup>使得如下不等式成立:||f(x<sub>1</sub>(t))‑f(x<sub>2</sub>(t))||≤||F(x<sub>1</sub>(t)‑x<sub>2</sub>(t))||其中,x<sub>1</sub>(t)和x<sub>2</sub>(t)∈{x(t)|t为时间}为系统状态集合中的任意两个状态,F为非线性项的增益阵;外部模型描述干扰d<sub>1</sub>(t)由下列外部干扰模型∑<sub>2</sub>表示:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>&Sigma;</mi><mn>2</mn></msub><mo>:</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mover><mi>w</mi><mo>&CenterDot;</mo></mover><mo>(</mo><mi>t</mi><mo>)</mo><mo>=</mo><mi>W</mi><mi>w</mi><mo>(</mo><mi>t</mi><mo>)</mo><mo>+</mo><mi>&delta;</mi><mo>(</mo><mi>t</mi><mo>)</mo></mtd></mtr><mtr><mtd><msub><mi>d</mi><mn>1</mn></msub><mo>(</mo><mi>t</mi><mo>)</mo><mo>=</mo><mi>V</mi><mi>w</mi><mo>(</mo><mi>t</mi><mo>)</mo></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000954861480000012.GIF" wi="581" he="165" /></maths>其中,δ(t)为范数有界变量,w(t)为可建模干扰子系统的状态变量,W表示可建模干扰子系统的系统矩阵,V为可建模干扰子系统的输出矩阵;系统中的参数不确定性d<sub>2</sub>(t)满足d<sub>2</sub>(t)=ρ(y,t)θ,其中ρ(y,t)为已知的连续函数,θ为未知的待估计参数;第二步,设计干扰观测器补偿可建模干扰针对多源干扰系统∑<sub>1</sub>中的外部可建模干扰d<sub>1</sub>(t),设计干扰观测器对其进行估计,并求得估计值<img file="FDA0000954861480000021.GIF" wi="130" he="79" />干扰观测器结构如下:<maths num="0003" id="cmaths0003"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mover><mover><mi>w</mi><mo>^</mo></mover><mo>&CenterDot;</mo></mover><mo>=</mo><mi>W</mi><mover><mi>w</mi><mo>^</mo></mover><mo>(</mo><mi>t</mi><mo>)</mo><mo>+</mo><msub><mi>K</mi><mn>2</mn></msub><mo>(</mo><mi>y</mi><mo>(</mo><mi>t</mi><mo>)</mo><mo>-</mo><mover><mi>y</mi><mo>^</mo></mover><mo>(</mo><mi>t</mi><mo>)</mo><mo>)</mo></mtd></mtr><mtr><mtd><msub><mover><mi>d</mi><mo>^</mo></mover><mn>1</mn></msub><mo>(</mo><mi>t</mi><mo>)</mo><mo>=</mo><mi>V</mi><mover><mi>w</mi><mo>^</mo></mover><mo>(</mo><mi>t</mi><mo>)</mo></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000954861480000022.GIF" wi="670" he="191" /></maths>其中,<img file="FDA0000954861480000023.GIF" wi="84" he="55" />为干扰观测器状态变量,<img file="FDA0000954861480000024.GIF" wi="102" he="71" />为d<sub>1</sub>(t)的估计值,K<sub>2</sub>∈R<sup>p×m</sup>为待定观测器增益矩阵,R<sup>p×m</sup>表示p×m维实矩阵空间,通过对(W B<sub>1</sub>V)配置负实部极点得到,p和m均为自然数;y(t)为多源干扰系统输出变量,<img file="FDA0000954861480000025.GIF" wi="85" he="55" />为自适应鲁棒滤波器输出变量;第三步,构造自适应更新律补偿参数不确定性针对多源干扰系统∑<sub>1</sub>中的参数不确定性干扰d<sub>2</sub>(t),设计自适应输入项u(t)对参数不确定性干扰进行补偿;第四步,设计H<sub>∞</sub>滤波器抑制范数有界干扰经过干扰补偿与自适应补偿后,系统Σ<sub>1</sub>中存在的不确定干扰为:干扰估计误差e<sub>w</sub>(t)和参数不确定性补偿误差e<sub>θ</sub>(t),范数有界干扰ω(t)、υ(t)以及δ(t);以上干扰均通过优化滤波系统的H<sub>∞</sub>范数来进行抑制,构造如下H<sub>∞</sub>滤波器:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mover><mover><mi>x</mi><mo>^</mo></mover><mo>&CenterDot;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>A</mi><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mi>F</mi><mi>f</mi><mrow><mo>(</mo><mover><mi>x</mi><mo>^</mo></mover><mo>(</mo><mi>t</mi><mo>)</mo><mo>)</mo></mrow><mo>+</mo><msub><mi>K</mi><mn>1</mn></msub><mo>&lsqb;</mo><mi>y</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mi>C</mi><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0000954861480000026.GIF" wi="931" he="79" /></maths>其中,<img file="FDA0000954861480000027.GIF" wi="797" he="83" /><img file="FDA0000954861480000028.GIF" wi="93" he="71" />为滤波器的状态变量,<img file="FDA0000954861480000029.GIF" wi="174" he="70" />为滤波器的非线性项,用于补偿对象模型中的非线性项f(x(t)),K<sub>1</sub>为待定滤波器增益阵,通过极小化如下性能指标求得:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>J</mi><mo>=</mo><mi>&alpha;</mi><mo>|</mo><mo>|</mo><msub><mi>T</mi><mi>&infin;</mi></msub><mo>|</mo><msubsup><mo>|</mo><mi>&infin;</mi><mn>2</mn></msubsup></mrow>]]></math><img file="FDA00009548614800000210.GIF" wi="254" he="87" /></maths>式中:<img file="FDA00009548614800000213.GIF" wi="243" he="55" />其中e=[e<sub>w</sub>(t)<sup>T</sup> e<sub>θ</sub>(t)<sup>T</sup> ω(t)<sup>T</sup> υ(t)<sup>T</sup> δ(t)<sup>T</sup>]<sup>T</sup>,<img file="FDA00009548614800000211.GIF" wi="390" he="55" /><img file="FDA00009548614800000212.GIF" wi="622" he="71" />α为由设计者给定的正实数,表示范数有界干扰e的抑制水平;第五步,基于干扰观测器、自适应输入补偿项和H<sub>∞</sub>滤波器,构造具有抗多源干扰能力的复合分层自适应滤波器滤波器采用复合分层的结构设计,对系统中的外部可建模干扰d<sub>1</sub>(t)进行抵消,补偿参数不确定性干扰d<sub>2</sub>(t),同时对范数有界干扰ω(t)、υ(t)、δ(t),干扰估计误差e<sub>w</sub>(t),参数不确定性补偿误差e<sub>θ</sub>(t)进行抑制,复合分层自适应滤波器结构如下:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mover><mover><mi>x</mi><mo>^</mo></mover><mo>&CenterDot;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>A</mi><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mi>F</mi><mi>f</mi><mrow><mo>(</mo><mover><mi>x</mi><mo>^</mo></mover><mo>(</mo><mi>t</mi><mo>)</mo><mo>)</mo></mrow><mo>+</mo><msub><mi>K</mi><mn>1</mn></msub><mo>&lsqb;</mo><mi>y</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mi>C</mi><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>+</mo><msub><mi>B</mi><mn>1</mn></msub><msub><mover><mi>d</mi><mo>^</mo></mover><mn>1</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>.</mo></mrow>]]></math><img file="FDA0000954861480000031.GIF" wi="1329" he="87" /></maths>
地址 100191 北京市海淀区学院路37号