发明名称 基于集中式MIMO雷达的多目标认知跟踪方法
摘要 本发明属于雷达系统多目标跟踪技术领域,公开了复杂环境下基于集中式MIMO雷达的多目标认知跟踪方法。该基于集中式MIMO雷达的多目标认知跟踪方法包括以下步骤:设定第q个目标在k时刻的状态<img file="DDA0000461152810000011.GIF" wi="61" he="78" />以及其概率密度函数;设定观测矩阵,得出观测矩阵的条件概率密度;计算k时刻的贝叶斯信息矩阵<img file="DDA0000461152810000012.GIF" wi="142" he="103" />;得出k时刻与k+1时刻BIM之间的递推关系式,在给定发射功率的假设下,计算k+1时刻第q目标跟踪误差的贝叶斯克拉美罗下界;建立功率分配模型,并进行求解,根据求解结果在k+1时刻向第q个目标发射具有对应功率的波束。
申请公布号 CN103760556B 申请公布日期 2016.03.23
申请号 CN201410032533.8 申请日期 2014.01.23
申请人 西安电子科技大学 发明人 刘宏伟;纠博;严俊坤;蒲文强;戴奉周;周生华
分类号 G01S13/66(2006.01)I 主分类号 G01S13/66(2006.01)I
代理机构 西安睿通知识产权代理事务所(特殊普通合伙) 61218 代理人 惠文轩
主权项 基于集中式MIMO雷达的多目标认知跟踪方法,其特征在于,包括以下步骤:S1:所述多目标为第1个目标至第Q个目标,Q为大于1的自然数;将第q个目标在k时刻的状态表示为<img file="FDA0000844799600000011.GIF" wi="92" he="76" />q取1至Q,k为大于0的自然数;将第q个目标在k时刻的状态的概率密度表示为<img file="FDA0000844799600000012.GIF" wi="174" he="102" />S2:设定观测矩阵Z<sub>q,k</sub>,<img file="FDA0000844799600000013.GIF" wi="382" he="118" />其中,<img file="FDA0000844799600000014.GIF" wi="93" he="86" />为观测向量,m<sub>q,k</sub>表示第q个目标在k时刻跟踪波门内检测到的点数;m<sub>q,k</sub>的概率密度为p(m<sub>q,k</sub>);根据Z<sub>q,k</sub>、p(m<sub>q,k</sub>)、以及<img file="FDA0000844799600000015.GIF" wi="189" he="102" />得出观测矩阵Z<sub>q,k</sub>的条件概率密度<img file="FDA0000844799600000016.GIF" wi="398" he="110" />S3:k时刻第q个目标状态的贝叶斯信息矩阵表示为<img file="FDA0000844799600000017.GIF" wi="174" he="100" />并用概率数据互联方法获取第q个目标在k时刻目标状态<img file="FDA0000844799600000018.GIF" wi="61" he="77" />的估计;S4:根据p(m<sub>q,k</sub>)、<img file="FDA0000844799600000019.GIF" wi="181" he="101" />以及<img file="FDA00008447996000000110.GIF" wi="398" he="110" />得出<img file="FDA00008447996000000111.GIF" wi="150" he="102" />与<img file="FDA00008447996000000112.GIF" wi="187" he="100" />之间的递推关系式;设置初始时刻每个波束的发射功率,并根据所述初始时刻第q个波束的发射功率得出初始时刻第q个目标的状态的贝叶斯信息矩阵<img file="FDA00008447996000000113.GIF" wi="181" he="102" />预测出第q个目标在k+1时刻的跟踪误差的贝叶斯克拉美罗下界<img file="FDA00008447996000000114.GIF" wi="401" he="125" /><img file="FDA00008447996000000115.GIF" wi="666" he="126" />P<sub>q,k+1</sub>表示k+1时刻集中式MIMO雷达第q个波束的发射功率;S5:建立如下功率分配模型:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><munder><mi>min</mi><msub><mi>P</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub></munder><mi>F</mi><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000844799600000021.GIF" wi="285" he="111" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><msub><mover><mi>P</mi><mo>&OverBar;</mo></mover><mrow><mi>q</mi><mi>m</mi><mi>i</mi><mi>n</mi></mrow></msub><mo>&le;</mo><msub><mi>P</mi><mrow><mi>q</mi><mo>,</mo><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>&le;</mo><msub><mover><mi>P</mi><mo>&OverBar;</mo></mover><mrow><mi>q</mi><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>,</mo><mi>q</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>Q</mi></mrow>]]></math><img file="FDA0000844799600000022.GIF" wi="902" he="93" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><munderover><mo>&Sigma;</mo><mrow><mi>q</mi><mo>=</mo><mn>1</mn></mrow><mi>Q</mi></munderover><msub><mi>P</mi><mrow><mi>q</mi><mo>,</mo><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mi>P</mi><mrow><mi>t</mi><mi>o</mi><mi>t</mi><mi>a</mi><mi>l</mi></mrow></msub></mrow>]]></math><img file="FDA0000844799600000023.GIF" wi="342" he="166" /></maths>其中,<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>F</mi><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>=</mo><munder><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mi>q</mi></munder><msqrt><mrow><mi>T</mi><mi>r</mi><mrow><mo>(</mo><mo>&lsqb;</mo><msub><mi>C</mi><mrow><mi>C</mi><mi>R</mi><mi>L</mi><mi>B</mi></mrow></msub><mo>(</mo><msub><mi>P</mi><mrow><mi>q</mi><mo>,</mo><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><msub><mo>|</mo><msubsup><mi>&xi;</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mi>q</mi></msubsup></msub><mo>&rsqb;</mo><mo>)</mo></mrow></msqrt><mo>,</mo></mrow>]]></math><img file="FDA0000844799600000024.GIF" wi="956" he="173" /></maths>F(P<sub>k+1</sub>)表示精度最差的目标的跟踪误差;P<sub>k+1</sub>表示k+1时刻各个波束的发射功率;<img file="FDA0000844799600000025.GIF" wi="117" he="92" />表示第q个波束的设定最小发射功率,<img file="FDA0000844799600000026.GIF" wi="124" he="92" />表示第q个波束的设定最大发射功率;P<sub>total</sub>为集中式MIMO雷达各波束设定的总发射功率;S6:通过求解所述功率分配模型,得出k+1时刻各个波束的发射功率P<sub>k+1</sub>;在k+1时刻向各目标发射具有对应功率的波束;S7:重复执行步骤S4至步骤S6,集中式MIMO雷达对每个目标进行持续的跟踪,直到目标超出雷达的跟踪范围;在步骤S1中,<img file="FDA0000844799600000027.GIF" wi="638" he="86" />其中,<img file="FDA0000844799600000028.GIF" wi="228" he="101" />表示第q个目标在k时刻的位置,<img file="FDA0000844799600000029.GIF" wi="227" he="102" />表示第q个目标在k时刻的速度;<img file="FDA00008447996000000210.GIF" wi="70" he="79" />为第q个目标在k时刻的雷达散射截面积;T表示观测时间间隔;目标运动模型表示为:<img file="FDA00008447996000000211.GIF" wi="389" he="85" />其中,<img file="FDA00008447996000000212.GIF" wi="70" he="86" />为与第q个目标对应的状态转移矩阵,<img file="FDA00008447996000000213.GIF" wi="100" he="78" />表示零均值的高斯白噪声,其协方差矩阵为<img file="FDA00008447996000000214.GIF" wi="165" he="86" />在步骤S2中,若跟踪波门内的观测值来源于目标时,<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msubsup><mi>z</mi><mrow><mi>q</mi><mo>,</mo><mi>k</mi></mrow><mi>j</mi></msubsup><mo>=</mo><msub><mi>h</mi><mrow><mi>q</mi><mo>,</mo><mi>k</mi></mrow></msub><mrow><mo>(</mo><msubsup><mi>&xi;</mi><mi>k</mi><mi>q</mi></msubsup><mo>)</mo></mrow><mo>+</mo><msub><mi>w</mi><mrow><mi>q</mi><mo>,</mo><mi>k</mi></mrow></msub></mrow>]]></math><img file="FDA00008447996000000215.GIF" wi="486" he="102" /></maths>其中,h<sub>q,k</sub>(·)为观测函数,w<sub>q,k</sub>为零均值的高斯白噪声,其协方差矩阵为Σ<sub>q,k</sub>;反之,当跟踪波门内的观测值来源于虚警时,<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msubsup><mi>z</mi><mrow><mi>q</mi><mo>,</mo><mi>k</mi></mrow><mi>j</mi></msubsup><mo>=</mo><msub><mi>&upsi;</mi><mrow><mi>q</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>;</mo></mrow>]]></math><img file="FDA00008447996000000216.GIF" wi="261" he="87" /></maths>虚警在跟踪波门内服从均匀分布,其个数服从泊松分布。
地址 710071 陕西省西安市太白南路2号