发明名称 一种机载雷达地面动静目标稳定跟踪方法
摘要 本发明涉及一种机载雷达地面动静目标稳定跟踪方法,根据动、静目标特点,采用同时发射非线性调频和步进频信号分别对动、静目标进行检测。非线性调频信号主要用于检测动目标,步进频信号主要用于检测静目标。在雷达的一个扫描帧可以同时得到地面动目标和静目标的一次点迹。与传统的方法相比,本发明能够同时关注地面动目标和静目标,反应地面目标的整体态势,当目标由动到静或由静到动时能持续跟踪,保证地面动、静目标的连续稳定跟踪。
申请公布号 CN104502907A 申请公布日期 2015.04.08
申请号 CN201410775618.5 申请日期 2014.12.15
申请人 西安电子工程研究所 发明人 段芳芳;张明;周凯;王伟;唐尧
分类号 G01S13/66(2006.01)I 主分类号 G01S13/66(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 一种机载雷达地面动静目标稳定跟踪方法,其特征在于步骤如下:步骤1:雷达发射非线性调频信号用于检测动目标,发射步进频信号检测静目标,在雷达的一个扫描帧同时得到地面动目标和静目标的一次点迹;步骤2:分别进行动目标和静目标的一次点迹进行去旁瓣和加权凝聚处理;步骤3:采用基于两级Hough变换的航迹起始方法将动和静目标的凝聚点分别建立各自的备选航迹;将前一帧备选航迹进行运动补偿,补偿到当前帧载机所在位置;步骤4:预测当前帧航迹的位置:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mover><mi>X</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>/</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mi>F</mi><mover><mi>X</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>/</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000633320040000011.GIF" wi="763" he="85" /></maths>其中:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>F</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mi>T</mi></mtd><mtd><msup><mi>T</mi><mn>2</mn></msup><mo>/</mo><mn>2</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mi>T</mi></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000633320040000012.GIF" wi="465" he="197" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>X</mi><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><mi>x</mi></mtd><mtd><mover><mi>x</mi><mo>&CenterDot;</mo></mover></mtd><mtd><mover><mi>x</mi><mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mrow></mover></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>,</mo></mrow>]]></math><img file="FDA0000633320040000013.GIF" wi="407" he="81" /></maths>x、<img file="FDA0000633320040000014.GIF" wi="45" he="62" />和<img file="FDA0000633320040000015.GIF" wi="38" he="69" />分别表示当前帧航迹的位置,速度和加速度分量;步骤5:重复步骤1~步骤2,得到下一个点迹的凝聚点;步骤6:以航迹的预测位置为中心,将落入波门内的量测数据和航迹相关联,得到相关组合;如果量测数据没有落在任何一个航迹的波门内,以该数据建立一个新的航迹;如果在航迹的波门内没有观测点,分两种情况:情况1:当动目标航迹与所有动目标凝聚点都不相关,则进行两假设跟踪,得到相关组合;假设1:目标未被检出,则将航迹按照之前的运动规律进行外推;假设2:目标静止,则进行动目标航迹与静目标凝聚点之间的相关,如果有相关点,则用其进行航迹更新,如果没有相关点,则保持航迹原来的位置不动,速度为零;情况2:当静目标航迹与所有静目标凝聚点都不相关,则进行两假设跟踪;假设1:目标未被检出,则保持航迹原来的位置不动,速度为零,目标仍然静止;假设2:目标转为运动,则进行静目标航迹与动目标凝聚点之间的相关,如果有相关点,则用其进行航迹更新,并估计航迹的速度,如果没有相关点,则保持航迹原来的位置不动,等待下一帧的相关,得到相关组合;将所有动目标航迹和静目标航迹的相关组合按照下式参与运算,找到一个全局代价最小的最优分配,用其进行航迹的最终更新;a<sup>*</sup>(k)=arg min<sub>a(k)</sub>C(k|a(k))其中:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>C</mi><mrow><mo>(</mo><mi>k</mi><mo>|</mo><mi>a</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>=</mo><msubsup><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><msub><mn>0</mn><mn>1</mn></msub><mo>,</mo><msub><mn>0</mn><mn>2</mn></msub><mo>,</mo><mn>1</mn></mrow><mrow><mi>M</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></msubsup><msubsup><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></msubsup><mi>a</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mi>c</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000633320040000021.GIF" wi="1183" he="108" /></maths>其中:c(k,m,n)是分配a(k,m,n)的代价值:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>c</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><mi>ln</mi><mrow><mo>(</mo><mfrac><mrow><mi>&Phi;</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow></mrow><mrow><mi>&Phi;</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>m</mi><mo>,</mo><mn>0</mn><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000633320040000022.GIF" wi="636" he="125" /></maths>当情况1时:Ф(k,m,n)计算如下:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>&Phi;</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>P</mi><mi>D</mi></msub><mi>P</mi><mo>{</mo><mo>|</mo><msub><mi>v</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo><mo>&GreaterEqual;</mo><msub><mi>v</mi><mn>0</mn></msub><mo>}</mo><mi>&Lambda;</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow></mtd><mtd><mi>m</mi><mo>></mo><mn>0</mn><mo>,</mo><mi>n</mi><mo>></mo><mn>0</mn></mtd></mtr><mtr><mtd><mi>&lambda;</mi></mtd><mtd><mi>m</mi><mo>></mo><mn>0</mn><mo>,</mo><mi>n</mi><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>P</mi><mi>D</mi></msub><mo>)</mo></mrow><mi>P</mi><mo>{</mo><mo>|</mo><msub><mi>v</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo><mo>&GreaterEqual;</mo><msub><mi>v</mi><mn>0</mn></msub><mo>}</mo><mi>&lambda;</mi></mtd><mtd><mi>m</mi><mo>=</mo><msub><mn>0</mn><mn>1</mn></msub><mo>,</mo><mi>n</mi><mo>></mo><mn>0</mn></mtd></mtr><mtr><mtd><mi>P</mi><mo>{</mo><mo>|</mo><msub><mi>v</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo><mo>&lt;</mo><msub><mi>v</mi><mn>0</mn></msub><mo>}</mo><mi>&lambda;</mi></mtd><mtd><mi>m</mi><mo>=</mo><msub><mn>0</mn><mn>2</mn></msub><mo>,</mo><mi>n</mi><mo>></mo><mn>0</mn></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000633320040000023.GIF" wi="1412" he="285" /></maths>其中:m&gt;0,n&gt;0为测量与航迹相关;m&gt;0,n=0为测量与已存在的航迹不相关,为虚警;m=0<sub>1</sub>,n&gt;0为假设测量1与航迹相关;m=0<sub>2</sub>,n&gt;0为假设测量2与航迹相关;Λ(k,m,n)是交互多模型估计的组合概率,v<sub>0</sub>是最小可检测速度,v<sub>r</sub>(k)是第k个扫描帧的径向速度(如果第k个扫描帧没有测量值则v<sub>r</sub>(k)是由第k‑1个扫描帧预测得来)。假设速度服从高斯分布,则可以计算出P{|v<sub>r</sub>(k)|≥v<sub>0</sub>}和P{|v<sub>r</sub>(k)|&lt;v<sub>0</sub>};虚警的空间密度λ为:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mi>&lambda;</mi><mo>=</mo><mfrac><msub><mi>P</mi><mi>FA</mi></msub><mrow><msup><mn>12</mn><mrow><mn>3</mn><mo>/</mo><mn>2</mn></mrow></msup><msub><mi>&sigma;</mi><mi>r</mi></msub><msub><mi>r&sigma;</mi><mi>&theta;</mi></msub><msub><mi>&sigma;</mi><mover><mi>r</mi><mo>&CenterDot;</mo></mover></msub></mrow></mfrac></mrow>]]></math><img file="FDA0000633320040000024.GIF" wi="380" he="122" /></maths>其中P<sub>FA</sub>是虚警概率,距离、方位、距离速度三维的测量精度分别为<img file="FDA0000633320040000025.GIF" wi="374" he="83" />和<img file="FDA0000633320040000026.GIF" wi="178" he="85" />当情况2时,Ф(k,m,n)计算如下:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mi>&Phi;</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>P</mi><mi>D</mi></msub><mi>P</mi><mo>{</mo><mo>|</mo><msub><mi>v</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo><mo>=</mo><mn>0</mn><mo>}</mo><mi>&Lambda;</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow></mtd><mtd><mi>m</mi><mo>></mo><mn>0</mn><mo>,</mo><mi>n</mi><mo>></mo><mn>0</mn></mtd></mtr><mtr><mtd><mi>&lambda;</mi></mtd><mtd><mi>m</mi><mo>></mo><mn>0</mn><mo>,</mo><mi>n</mi><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>P</mi><mi>D</mi></msub><mo>)</mo></mrow><mi>P</mi><mo>{</mo><mo>|</mo><msub><mi>v</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo><mo>=</mo><mn>0</mn><mo>}</mo><mi>&lambda;</mi></mtd><mtd><mi>m</mi><mo>=</mo><msub><mn>0</mn><mn>1</mn></msub><mo>,</mo><mi>n</mi><mo>></mo><mn>0</mn></mtd></mtr><mtr><mtd><mi>P</mi><mo>{</mo><mo>|</mo><msub><mi>v</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo><mo>!</mo><mo>=</mo><mn>0</mn><mo>}</mo><mi>&lambda;</mi></mtd><mtd><mi>m</mi><mo>=</mo><msub><mn>0</mn><mn>2</mn></msub><mo>,</mo><mi>n</mi><mo>></mo><mn>0</mn></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000633320040000031.GIF" wi="1402" he="315" /></maths>情况2需要计算P{|v<sub>r</sub>(k)|=0}和P{|v<sub>r</sub>(k)|!=0};所述波门包括距离波门、方位波门、俯仰波门和多普勒波门;步骤7、滤波:按照α‑β‑γ的滤波方法,对航迹进行滤波<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><mover><mi>X</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>/</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mover><mi>X</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>/</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><mi>K</mi><mo>[</mo><mi>y</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mi>H</mi><mover><mi>X</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>/</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>]</mo></mrow>]]></math><img file="FDA0000633320040000032.GIF" wi="1119" he="85" /></maths><maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><mover><mi>X</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>/</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mi>F</mi><mover><mi>X</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>/</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000633320040000033.GIF" wi="735" he="80" /></maths>其中增益矩阵<maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><mi>K</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mi>&alpha;</mi></mtd></mtr><mtr><mtd><mi>&beta;</mi><mo>/</mo><mi>T</mi></mtd></mtr><mtr><mtd><mi>&gamma;</mi><mo>/</mo><msup><mi>T</mi><mn>2</mn></msup></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000633320040000034.GIF" wi="299" he="196" /></maths>观测矩阵H=[1 0 0];步骤8、航迹管理:当一条航迹多次没有观测点造成多次外推时,将航迹删除;重复步骤1~步骤8,进行动静目标稳定跟踪。
地址 710100 陕西省西安市长安区凤栖东路