发明名称 系统误差下基于拓扑信息的雷达与ESM航迹关联方法
摘要 本发明公开了一种系统误差下基于拓扑信息的雷达与ESM航迹关联算法。对于异地配置的雷达与ESM,方位线相交会产生若干交叉定位点,其中包含错误的鬼点和正确关联的点。在得到所有交叉定位点后,本发明首先基于距离信息进行局部关联,排除大量不可能关联的情况,减少了计算量,节约了时间。然后对剩下的交叉定位点划分组合,基于拓扑信息分别计算不同组合与雷达量测的关联程度,关联程度最大的组合所对应的交叉定位点依次与相应的雷达量测关联,从而得到所有雷达与ESM目标的关联关系。相比于传统的航迹关联算法,本发明充分利用目标间的拓扑信息,能够很好的抑制系统误差对航迹关联造成的影响,提高正确关联率,减少误关联率和漏关联率。
申请公布号 CN106443664A 申请公布日期 2017.02.22
申请号 CN201610820652.9 申请日期 2016.09.13
申请人 中国人民解放军海军航空工程学院 发明人 关欣;彭彬彬;衣晓;孙贵东
分类号 G01S13/72(2006.01)I 主分类号 G01S13/72(2006.01)I
代理机构 代理人
主权项 一种系统误差下基于拓扑信息的雷达与ESM航迹关联方法,其特征在于,包含如下步骤:步骤一、获取交叉定位点:第i个雷达目标与第j个ESM目标的角度交叉定位点x<sub>ij</sub>和y<sub>ij</sub>坐标<maths num="0001"><math><![CDATA[<mrow><mfenced open = "" close = "}"><mtable><mtr><mtd><mrow><msub><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mo>&lsqb;</mo><msub><mi>x</mi><mi>A</mi></msub><mi>cot</mi><mrow><mo>(</mo><msubsup><mi>&theta;</mi><mi>i</mi><mi>A</mi></msubsup><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow><mo>-</mo><msub><mi>y</mi><mi>A</mi></msub><mo>&rsqb;</mo><mo>-</mo><mo>&lsqb;</mo><msub><mi>x</mi><mi>B</mi></msub><mi>cot</mi><mrow><mo>(</mo><msubsup><mi>&theta;</mi><mi>j</mi><mi>B</mi></msubsup><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow><mo>-</mo><msub><mi>y</mi><mi>B</mi></msub><mo>&rsqb;</mo></mrow><mrow><mi>cot</mi><mrow><mo>(</mo><msubsup><mi>&theta;</mi><mi>i</mi><mi>A</mi></msubsup><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow><mo>-</mo><mi>cot</mi><mrow><mo>(</mo><msubsup><mi>&theta;</mi><mi>j</mi><mi>B</mi></msubsup><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow></mrow></mfrac></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>y</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mo>&lsqb;</mo><msub><mi>y</mi><mi>A</mi></msub><mi>tan</mi><mrow><mo>(</mo><msubsup><mi>&theta;</mi><mi>i</mi><mi>A</mi></msubsup><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow><mo>-</mo><msub><mi>x</mi><mi>A</mi></msub><mo>&rsqb;</mo><mo>-</mo><mo>&lsqb;</mo><msub><mi>y</mi><mi>B</mi></msub><mi>tan</mi><mrow><mo>(</mo><msubsup><mi>&theta;</mi><mi>j</mi><mi>B</mi></msubsup><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow><mo>-</mo><msub><mi>x</mi><mi>B</mi></msub><mo>&rsqb;</mo></mrow><mrow><mi>tan</mi><mrow><mo>(</mo><msubsup><mi>&theta;</mi><mi>i</mi><mi>A</mi></msubsup><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow><mo>-</mo><mi>tan</mi><mrow><mo>(</mo><msubsup><mi>&theta;</mi><mi>j</mi><mi>B</mi></msubsup><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow></mrow></mfrac></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001114193610000011.GIF" wi="1549" he="366" /></maths>其中,r<sub>j</sub>(k)和<img file="FDA0001114193610000012.GIF" wi="126" he="71" />分别是k时刻第i个雷达目标的距离与方位测量值,(x<sub>A</sub>,y<sub>A</sub>)为雷达的坐标,(x<sub>B</sub>,y<sub>B</sub>)为ESM位的坐标,<img file="FDA0001114193610000013.GIF" wi="129" he="76" />为k时刻ESM第j个目标的角度测量值;步骤二、基于距离信息的局部关联:计算第i个雷达目标与第j个ESM目标的角度交叉定位点距雷达的距离d<sub>ij</sub>(k)<maths num="0002"><math><![CDATA[<mrow><msub><mi>d</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msqrt><mrow><msup><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>-</mo><msub><mi>x</mi><mi>A</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>y</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>-</mo><msub><mi>y</mi><mi>A</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001114193610000014.GIF" wi="1454" he="103" /></maths>令D<sub>ij</sub>(k)=d<sub>ij</sub>(k)‑r<sub>i</sub>(k),由于雷达角度与距离和ESM的角度之间的测量误差相互独立,所以D<sub>ij</sub>(k)的方差为<img file="FDA0001114193610000015.GIF" wi="396" he="83" />基于距离信息的局部粗关联处理规则为:当<img file="FDA0001114193610000016.GIF" wi="453" he="78" />时,即认为第i个雷达目标与第j个ESM目标来自同一个目标,即两航迹试验关联;当不满足时,认为第i个雷达目标与第j个ESM目标不可能来自同一个目标;步骤三、对剩余交叉定位点划分组合:按照每个雷达目标只与一个ESM目标关联的原则,对A进行拆分时遵循:从航迹关联矩阵A的每一行,选出唯一的一个1元素,作为航迹关联事件矩阵B<sub>l</sub>在该行和该列唯一的非零元素;其中,A为经过基于距离信息的局部关联之后,得到的雷达与ESM可能航迹关联矩阵;B<sub>l</sub>为航迹关联事件矩阵,是由按照上述原则拆分得到的,<img file="FDA0001114193610000017.GIF" wi="1405" he="302" />其中,l=1,2,…,N,<img file="FDA0001114193610000018.GIF" wi="47" he="70" />为航迹关联事件矩阵B<sub>l</sub>第i行第j列元素,为布尔型变量;步骤四、基于拓扑信息的全局灰关联:第l个航迹关联事件矩阵第i个目标与第i个雷达点迹的空间拓扑相关系数序列为<maths num="0003"><math><![CDATA[<mrow><msubsup><mi>&xi;</mi><mrow><mi>l</mi><mi>i</mi></mrow><mi>m</mi></msubsup><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mi>l</mi></munder><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mi>j</mi></munder><mo>|</mo><msubsup><mi>V</mi><mi>i</mi><mi>m</mi></msubsup><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><msubsup><mi>U</mi><mrow><mi>l</mi><mi>i</mi></mrow><mi>m</mi></msubsup><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>|</mo><mo>+</mo><mi>&rho;</mi><munder><mi>max</mi><mi>l</mi></munder><munder><mi>max</mi><mi>j</mi></munder><mo>|</mo><msubsup><mi>V</mi><mi>i</mi><mi>m</mi></msubsup><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><msubsup><mi>U</mi><mrow><mi>l</mi><mi>i</mi></mrow><mi>m</mi></msubsup><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>|</mo></mrow><mrow><mo>|</mo><msubsup><mi>V</mi><mi>i</mi><mi>m</mi></msubsup><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><msubsup><mi>U</mi><mrow><mi>l</mi><mi>i</mi></mrow><mi>m</mi></msubsup><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>|</mo><mo>+</mo><mi>&rho;</mi><munder><mi>max</mi><mi>l</mi></munder><munder><mi>max</mi><mi>j</mi></munder><mo>|</mo><msubsup><mi>V</mi><mi>i</mi><mi>m</mi></msubsup><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><msubsup><mi>U</mi><mrow><mi>l</mi><mi>i</mi></mrow><mi>m</mi></msubsup><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>|</mo></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001114193610000021.GIF" wi="1581" he="185" /></maths>其中,U<sub>li</sub>为第l个航迹关联事件矩阵第i个目标的空间拓扑描述为序列,V<sub>i</sub>为k时刻第i个雷达点迹的拓扑信息为其与所有邻居的位置差向量序列,i=1,2,…n,l=1,2,…,N,j=1,2,…,n‑1,m=1,2;由于第l个航迹关联事件矩阵有n个交叉点,每个交叉点包含n‑1个邻居,且每一个交叉定位点坐标序列都包含横纵轴两个相关系数序列分量,定义第l个航迹关联事件矩阵与雷达目标的空间拓扑的相似度的灰关联度为:<maths num="0004"><math><![CDATA[<mrow><msub><mi>&gamma;</mi><mi>l</mi></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></munderover><munderover><mo>&Sigma;</mo><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mn>2</mn></munderover><msubsup><mi>&xi;</mi><mrow><mi>l</mi><mi>i</mi></mrow><mi>m</mi></msubsup><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001114193610000022.GIF" wi="1390" he="135" /></maths>计算得到γ<sub>s</sub>=max{γ<sub>l</sub>|l=1,2,…,N}即第s个航迹关联事件矩阵对应的交叉定位点迹与雷达目标点迹空间拓扑相似程度最高,然后依据第s个航迹关联事件矩阵的非零元素判定多目标情况下的各雷达与ESM航迹关联情况。
地址 264001 山东省烟台市芝罘区二马路188号科研部