发明名称 一种基于量测存储的航迹起始方法
摘要 本文公开了一种基于历史量测信息的航迹起始方法。该算法针对窄发多收边扫描边跟踪雷达系统和被动跟踪系统在实际场景中的目标跟踪。现有的跟踪算法大都假设目标能够被连续观测到且不同目标的量测信息同步,本发明特点:考虑了边扫描边跟踪雷达系统中目标被观测到的时间间隔不定的情况,充分利用历史量测信息,具有更快的航迹起始速度和成功跟踪概率;本发明同时还可以解决集中式跟踪系统和被动跟踪系统中多传感器异步通信时航迹起始的问题。优点:充分利用了历史信息,实现简单,成功跟踪概率更高;有效解决了属于不同目标的量测信息异步时的航迹起始问题。
申请公布号 CN106199584A 申请公布日期 2016.12.07
申请号 CN201610519551.8 申请日期 2016.07.05
申请人 电子科技大学 发明人 易伟;陈方园;姜萌;方梓成;王佰录;李溯琪;崔国龙;孔令讲
分类号 G01S13/72(2006.01)I;G01S13/58(2006.01)I 主分类号 G01S13/72(2006.01)I
代理机构 电子科技大学专利中心 51203 代理人 张杨
主权项 一种基于量测存储的航迹起始方法,该方法包括如下步骤:步骤1、初始化相关参数,初始化参数主要包括:数据关联波门大小<img file="FDA0001040744910000011.GIF" wi="98" he="78" />速度波门大小(v<sub>min</sub>,v<sub>max</sub>),时间帧k=1,未关联量测矩阵<img file="FDA0001040744910000012.GIF" wi="145" he="61" />步骤2、利用恒虚警检测器处理雷达回波信息,得到量测信息集合<img file="FDA0001040744910000013.GIF" wi="418" he="78" />其中,n为检测到的量测个数;步骤3、如果k=1,将所有量测z<sub>1</sub>∈Z<sub>1</sub>保存到矩阵<img file="FDA0001040744910000014.GIF" wi="59" he="61" />中,并令k=k+1并跳转至步骤2;步骤4、若当前时刻无航迹,将所有量测z<sub>k</sub>∈Z<sub>k</sub>保存到矩阵<img file="FDA0001040744910000015.GIF" wi="91" he="63" />中,并跳转至步骤7,若存在航迹,进行步骤5;步骤5、第k帧量测与航迹进行数据关联:5.1.预测目标的状态、估计误差协方差矩阵:<maths num="0001"><math><![CDATA[<mrow><msubsup><mi>x</mi><mi>k</mi><mo>-</mo></msubsup><mo>=</mo><msub><mi>F</mi><mi>k</mi></msub><msub><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow>]]></math><img file="FDA0001040744910000016.GIF" wi="214" he="62" /></maths><maths num="0002"><math><![CDATA[<mrow><msubsup><mi>P</mi><mi>k</mi><mo>-</mo></msubsup><mo>=</mo><msub><mi>F</mi><mi>k</mi></msub><msub><mi>P</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><msup><msub><mi>F</mi><mi>k</mi></msub><mo>&prime;</mo></msup><mo>+</mo><msub><mi>Q</mi><mi>k</mi></msub></mrow>]]></math><img file="FDA0001040744910000017.GIF" wi="390" he="55" /></maths>其中,<img file="FDA0001040744910000018.GIF" wi="318" he="63" />为k‑1时刻目标的状态估计值,<img file="FDA0001040744910000019.GIF" wi="52" he="54" />为k时刻状态一步预测值,(x,y)和<img file="FDA00010407449100000110.GIF" wi="123" he="63" />分别为目标的坐标和速度,P<sub>k‑1</sub>为k‑1时刻的状态估计误差协方差矩阵,<img file="FDA00010407449100000111.GIF" wi="58" he="55" />为k时刻一步状态预测误差协方差矩阵,第Q<sub>k</sub>为当前时刻过程噪声,F<sub>k</sub>为状态转移矩阵;5.2.第k帧量测与已有航迹进行数据关联,如果满足下式,则认为关联成功:v'(k)S<sup>‑1</sup>(k)v(k)&lt;γ<sup>2</sup>其中,<img file="FDA00010407449100000112.GIF" wi="403" he="71" />表示量测与预测位置的差,S(k)表示v(k)的协方差矩阵,H表示观测矩阵;步骤6、利用第k帧关联到的量测z<sub>k</sub>更新目标状态:<maths num="0003"><math><![CDATA[<mrow><msub><mover><mi>x</mi><mo>^</mo></mover><mi>k</mi></msub><mo>=</mo><msubsup><mi>x</mi><mi>k</mi><mo>-</mo></msubsup><mo>+</mo><msub><mi>G</mi><mi>k</mi></msub><mrow><mo>(</mo><msub><mi>z</mi><mi>k</mi></msub><mo>-</mo><msub><mi>H</mi><mi>k</mi></msub><msubsup><mi>x</mi><mi>k</mi><mo>-</mo></msubsup><mo>)</mo></mrow></mrow>]]></math><img file="FDA00010407449100000113.GIF" wi="509" he="79" /></maths><maths num="0004"><math><![CDATA[<mrow><msub><mi>P</mi><mi>k</mi></msub><mo>=</mo><mrow><mo>(</mo><mi>I</mi><mo>-</mo><msub><mi>G</mi><mi>k</mi></msub><msub><mi>H</mi><mi>k</mi></msub><mo>)</mo></mrow><msubsup><mi>P</mi><mi>k</mi><mo>-</mo></msubsup></mrow>]]></math><img file="FDA00010407449100000114.GIF" wi="381" he="62" /></maths>其中,G<sub>k</sub>为卡尔曼增益,H<sub>k</sub>为量测矩阵,I为单位矩阵,<img file="FDA00010407449100000115.GIF" wi="49" he="61" />为k时刻目标的状态估计值,P<sub>k</sub>为k时刻状态估计误差协方差矩阵;步骤7、未关联到航迹的量测<img file="FDA0001040744910000021.GIF" wi="45" he="63" />与历史量测<img file="FDA0001040744910000022.GIF" wi="203" he="62" />起始新航迹:7.1.若满足下式,则起始新航迹:<maths num="0005"><math><![CDATA[<mrow><msub><mi>v</mi><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow></msub><mo>&lt;</mo><mo>|</mo><mo>|</mo><mfrac><mrow><msubsup><mi>z</mi><mi>k</mi><mi>u</mi></msubsup><mo>-</mo><msubsup><mi>z</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>u</mi></msubsup></mrow><mrow><msub><mi>t</mi><mi>k</mi></msub><mo>-</mo><msub><mi>t</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow></mfrac><mo>|</mo><mo>|</mo><mo>&lt;</mo><msub><mi>v</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub></mrow>]]></math><img file="FDA0001040744910000023.GIF" wi="470" he="149" /></maths>其中,t<sub>k</sub>为当前帧量测的观测时间,t<sub>k‑1</sub>为历史量测的观测时间,不同历史量测的观测时间可能不同,v<sub>min</sub>,v<sub>max</sub>为速度门限,该速度门限根据事情情况设定;7.2.更新量测矩阵<img file="FDA0001040744910000024.GIF" wi="107" he="63" />用来起始航迹的量测<img file="FDA0001040744910000025.GIF" wi="72" he="63" />从<img file="FDA0001040744910000026.GIF" wi="78" he="62" />中删除,同时,未起始航迹的量测<img file="FDA0001040744910000027.GIF" wi="168" he="63" />保存到集合<img file="FDA0001040744910000028.GIF" wi="80" he="63" />中:<maths num="0006"><math><![CDATA[<mrow><msubsup><mi>Z</mi><mi>k</mi><mi>u</mi></msubsup><mo>=</mo><mrow><mo>(</mo><msubsup><mi>Z</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>u</mi></msubsup><mo>-</mo><msubsup><mi>z</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>u</mi></msubsup><mo>)</mo></mrow><mo>&cup;</mo><msup><msubsup><mi>z</mi><mi>k</mi><mi>u</mi></msubsup><mo>&prime;</mo></msup></mrow>]]></math><img file="FDA0001040744910000029.GIF" wi="411" he="77" /></maths>步骤8、令k=k+1,如果k&lt;=K,其中K为总的观测帧数,返回步骤2。
地址 611731 四川省成都市高新区(西区)西源大道2006号