主权项 |
一种系统误差下基于拓扑信息的雷达与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>[</mo><msub><mi>x</mi><mi>A</mi></msub><mi>cot</mi><mrow><mo>(</mo><msubsup><mi>θ</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>]</mo><mo>-</mo><mo>[</mo><msub><mi>x</mi><mi>B</mi></msub><mi>cot</mi><mrow><mo>(</mo><msubsup><mi>θ</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>]</mo></mrow><mrow><mi>cot</mi><mrow><mo>(</mo><msubsup><mi>θ</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>θ</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>[</mo><msub><mi>y</mi><mi>A</mi></msub><mi>tan</mi><mrow><mo>(</mo><msubsup><mi>θ</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>]</mo><mo>-</mo><mo>[</mo><msub><mi>y</mi><mi>B</mi></msub><mi>tan</mi><mrow><mo>(</mo><msubsup><mi>θ</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>]</mo></mrow><mrow><mi>tan</mi><mrow><mo>(</mo><msubsup><mi>θ</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>θ</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>ξ</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>ρ</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>ρ</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>γ</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>Σ</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><munderover><mo>Σ</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></munderover><munderover><mo>Σ</mo><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mn>2</mn></munderover><msubsup><mi>ξ</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航迹关联情况。 |