发明名称 基于目标拓扑信息的AIS与雷达航迹抗差关联算法
摘要 本发明公开了一种基于目标拓扑信息的AIS与雷达航迹抗差关联算法,该技术属于AIS与雷达信息融合技术领域。现有的AIS与雷达航迹关联算法一般都假设雷达不存在系统误差,然而这是与实际情况相违背的,因而在实际工程应用中现有算法的正确航迹关联率都比较低。为了有效解决实际应用中的AIS与雷达航迹抗差关联问题,本发明给出了一种雷达系统存在系统误差的情况与AIS航迹进行关联的算法。该算法首先根据目标的拓扑关系构建了模糊因素集,然后根据模糊因素的重要性分配权重,最后给出模糊航迹关联的抗差准则以及多义性处理,并最终实现航迹的正确关联。
申请公布号 CN102707276B 申请公布日期 2015.11.25
申请号 CN201210175557.X 申请日期 2012.05.23
申请人 中国人民解放军海军航空工程学院 发明人 何友;王国宏;朱洪伟;唐小明;熊伟
分类号 G01S13/58(2006.01)I;G01S13/60(2006.01)I;G01S7/40(2006.01)I 主分类号 G01S13/58(2006.01)I
代理机构 代理人
主权项 一种基于目标拓扑信息的AIS与雷达航迹抗差关联方法,包括模糊因素集的构建与抗差关联模型的构建,以及基于模糊信息处理理论的航迹关联方法,其特征在于:步骤(1):通过AIS与雷达上报的目标之间的拓扑关系构建航迹的模糊因素集,定义<img file="FSB0000138395730000011.GIF" wi="303" he="67" />分别表示AIS上报的第i条目标航迹以及雷达上报的第j条目标航迹在第k时刻对应的相对拓扑向量;AIS中1×M维的相对拓扑向量<img file="FSB0000138395730000012.GIF" wi="123" he="73" />设定如下<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mi>T</mi><mi>A</mi><mi>i</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mo>[</mo><munderover><mi>&Sigma;</mi><mrow><mi>s</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mn>1</mn></msub></munderover><msubsup><mi>&rho;</mi><mi>A</mi><msubsup><mi>ii</mi><mi>s</mi><mn>1</mn></msubsup></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><munderover><mi>&Sigma;</mi><mrow><mi>s</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>m</mi></msub></munderover><msubsup><mi>&rho;</mi><mi>A</mi><msubsup><mi>ii</mi><mi>s</mi><mi>m</mi></msubsup></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><munderover><mi>&Sigma;</mi><mrow><mi>s</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>M</mi></msub></munderover><msubsup><mi>&rho;</mi><mi>A</mi><msubsup><mi>ii</mi><mi>s</mi><mi>M</mi></msubsup></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSB0000138395730000013.GIF" wi="1567" he="138" /></maths>式中,<img file="FSB0000138395730000014.GIF" wi="960" he="105" />表示k时刻AIS的待关联目标i点迹与目标i<sub>s</sub>点迹间的欧氏距离,N<sub>m</sub>表示落入象限m中的目标数;定义模糊因素1为<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>u</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mn>1</mn><mo>-</mo><mi>Cor</mi><mrow><mo>(</mo><msubsup><mi>T</mi><mi>A</mi><mi>i</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo><msubsup><mi>T</mi><mi>R</mi><mi>j</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>;</mo><mi>i</mi><mo>&Element;</mo><msub><mi>U</mi><mi>A</mi></msub><mo>,</mo><mi>j</mi><mo>&Element;</mo><msub><mi>U</mi><mi>R</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSB0000138395730000015.GIF" wi="1452" he="86" /></maths>式中,<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>Cor</mi><mrow><mo>(</mo><msubsup><mi>T</mi><mi>A</mi><mi>i</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo><msubsup><mi>T</mi><mi>R</mi><mi>j</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msubsup><mi>T</mi><mi>A</mi><mi>i</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msubsup><mi>T</mi><mi>R</mi><mi>j</mi></msubsup><msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mi>T</mi></msup></mrow><msqrt><mo>|</mo><msubsup><mi>T</mi><mi>A</mi><mi>i</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo><mo>|</mo><msubsup><mi>T</mi><mi>R</mi><mi>j</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo></msqrt></mfrac><mo>;</mo></mrow>]]></math><img file="FSB0000138395730000016.GIF" wi="761" he="182" /></maths>对于各目标航速间的信息模糊因素,采用各目标航速间的欧氏距离,则模糊因素2为<maths num="0004" id="cmaths0004"><math><![CDATA[<mfenced open='' close='' separators=''><mtable><mtr><mtd><msub><mi>u</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mo>|</mo><msubsup><mi>v</mi><mi>A</mi><mi>i</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msubsup><mi>v</mi><mi>R</mi><mi>j</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo></mtd></mtr><mtr><mtd><mo>=</mo><mo>|</mo><msup><mrow><mo>(</mo><msup><mrow><mo>(</mo><msubsup><mover><mi>x</mi><mo>.</mo></mover><mi>A</mi><mi>i</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msubsup><mover><mi>y</mi><mo>.</mo></mover><mi>A</mi><mi>i</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow><mfrac><mn>1</mn><mn>2</mn></mfrac></msup><mo>-</mo><msup><mrow><mo>(</mo><msup><mrow><mo>(</mo><msubsup><mover><mi>x</mi><mover><mo>.</mo><mo>^</mo></mover></mover><mi>R</mi><mi>j</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msubsup><mover><mi>y</mi><mover><mo>.</mo><mo>^</mo></mover></mover><mi>R</mi><mi>j</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow><mfrac><mn>1</mn><mn>2</mn></mfrac></msup><mo>|</mo><mo>,</mo><mi>i</mi><mo>&Element;</mo><msub><mi>U</mi><mi>A</mi></msub><mo>,</mo><mi>j</mi><mo>&Element;</mo><msub><mi>U</mi><mi>R</mi></msub></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mfenced>]]></math><img file="FSB0000138395730000017.GIF" wi="1737" he="300" /></maths>对于各目标航向变化率间的信息模糊因素,采用各目标航向变化率间的欧氏距离,则模糊因素3为<maths num="0005" id="cmaths0005"><math><![CDATA[<mfenced open='' close='' separators=''><mtable><mtr><mtd><msub><mi>u</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mo>|</mo><mrow><mo>(</mo><msubsup><mi>&theta;</mi><mi>A</mi><mi>i</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msubsup><mi>&theta;</mi><mi>A</mi><mi>i</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mrow><mo>(</mo><msubsup><mi>&theta;</mi><mi>R</mi><mi>j</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msubsup><mi>&theta;</mi><mi>R</mi><mi>j</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>)</mo></mrow><mo>|</mo></mrow><mi>T</mi></mfrac></mtd></mtr><mtr><mtd><mfrac><mrow><mo>|</mo><mrow><mo>(</mo><msup><mi>tg</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>[</mo><mfrac><mrow><msubsup><mover><mi>y</mi><mo>.</mo></mover><mi>A</mi><mi>i</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mrow><msubsup><mover><mi>x</mi><mo>.</mo></mover><mi>A</mi><mi>i</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></mfrac><mo>]</mo><mo>-</mo><msup><mi>tg</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>[</mo><mfrac><mrow><msubsup><mover><mi>y</mi><mo>.</mo></mover><mi>A</mi><mi>i</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><mrow><msubsup><mover><mi>x</mi><mo>.</mo></mover><mi>A</mi><mi>i</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac><mo>]</mo><mo>)</mo></mrow><mo>-</mo><mrow><mo>(</mo><msup><mi>tg</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>[</mo><mfrac><mrow><msubsup><mover><mi>y</mi><mover><mo>.</mo><mo>^</mo></mover></mover><mi>R</mi><mi>j</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mrow><msubsup><mover><mi>x</mi><mover><mo>.</mo><mo>^</mo></mover></mover><mi>R</mi><mi>j</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></mfrac><mo>]</mo><mo>-</mo><msup><mi>tg</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>[</mo><mfrac><mrow><msubsup><mover><mi>y</mi><munder><mo>^</mo><mo>.</mo></munder></mover><mi>R</mi><mi>j</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><mrow><msubsup><mover><mi>x</mi><mover><mo>.</mo><mo>^</mo></mover></mover><mi>R</mi><mi>j</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac><mo>]</mo><mo>)</mo></mrow><mo>|</mo></mrow><mi>T</mi></mfrac><mo>,</mo><mi>i</mi><mo>&Element;</mo><msub><mi>U</mi><mi>A</mi></msub><mo>,</mo><mi>j</mi><mo>&Element;</mo><msub><mi>U</mi><mi>R</mi></msub></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mfenced>]]></math><img file="FSB0000138395730000018.GIF" wi="1908" he="374" /></maths>其中,<img file="FSB0000138395730000019.GIF" wi="512" he="97" />为k时刻雷达对目标j的状态估计;由模糊因素1、2、3组成的集合为模糊因素集;步骤(2):通过建立的模糊因素集构建雷达与AIS航迹之间的抗差关联模型,选取正态型模糊隶属度函数,基于第l个模糊因素判为两目标航迹相似的隶属度为<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>&mu;</mi><mi>l</mi></msub><mo>=</mo><mi>&mu;</mi><mrow><mo>(</mo><msub><mi>&mu;</mi><mi>l</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><msub><mi>&tau;</mi><mi>l</mi></msub><mrow><mo>(</mo><mfrac><msubsup><mi>u</mi><mi>l</mi><mn>2</mn></msubsup><msubsup><mi>&sigma;</mi><mi>l</mi><mn>2</mn></msubsup></mfrac><mo>)</mo></mrow><mo>)</mo></mrow><mo>,</mo><mi>l</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSB00001383957300000110.GIF" wi="1564" he="181" /></maths>式中,σ<sub>l</sub>是模糊因素集<img file="FSB00001383957300000111.GIF" wi="44" he="68" />中模糊因素l的展度,而τ<sub>l</sub>是调整度;采用加权平均法进行综合评价,通过下式计算两目标航迹间的综合相似度<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>f</mi><mi>ij</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>l</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>a</mi><mi>l</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msub><mi>&mu;</mi><mi>l</mi></msub><mo>;</mo><mi>i</mi><mo>&Element;</mo><msub><mi>U</mi><mi>A</mi></msub><mo>,</mo><mi>j</mi><mo>&Element;</mo><msub><mi>U</mi><mi>R</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSB0000138395730000021.GIF" wi="1488" he="137" /></maths>其中,a<sub>l</sub>为第l个模糊因素u<sub>l</sub>对应的权;由此,对k时刻AIS的n<sub>A</sub>条目标航迹和雷达的n<sub>R</sub>条目标航迹构建以下模糊关联矩阵<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><munder><mi>F</mi><mo>~</mo></munder><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mtext>=</mtext><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>f</mi><mn>11</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd><mtd><msub><mi>f</mi><mn>12</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>f</mi><mrow><mn>1</mn><msub><mi>n</mi><mi>R</mi></msub></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>f</mi><mn>21</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd><mtd><msub><mi>f</mi><mn>22</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>f</mi><mrow><mn>2</mn><msub><mi>n</mi><mi>R</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>f</mi><mrow><msub><mi>n</mi><mi>A</mi></msub><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd><mtd><msub><mi>f</mi><mrow><msub><mi>n</mi><mi>A</mi></msub><mn>2</mn></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>f</mi><mrow><msub><mi>n</mi><mi>A</mi></msub><msub><mi>n</mi><mi>R</mi></msub></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSB0000138395730000022.GIF" wi="1551" he="335" /></maths>获得模糊关联矩阵式(7)后,在关联矩阵<img file="FSB0000138395730000023.GIF" wi="124" he="67" />中找出最大元素,且最大元素f<sub>ij</sub>(k)>ε,判定航迹i和j为试验关联航迹,并划去该元素所对应的所有行列元素;由此得到新的关联矩阵<img file="FSB0000138395730000024.GIF" wi="157" he="74" />但原矩阵的行列号不变,对其重复进行上述过程,得到矩阵<img file="FSB0000138395730000025.GIF" wi="215" he="75" />直至<img file="FSB0000138395730000026.GIF" wi="123" he="63" />的所有元素均小于阈值ε为止,剩余的元素的行列号对应为不关联矩阵;根据得到的航迹试验关联矩阵采用双门限航迹关联算法进行航迹关联与多义性处理实现AIS与雷达的航迹抗差关联。
地址 264001 山东省烟台市二马路188号海军航空工程学院