CN110133636A - 一种基于区域相关度的抗差航迹关联方法 - Google Patents

一种基于区域相关度的抗差航迹关联方法 Download PDF

Info

Publication number
CN110133636A
CN110133636A CN201910426229.4A CN201910426229A CN110133636A CN 110133636 A CN110133636 A CN 110133636A CN 201910426229 A CN201910426229 A CN 201910426229A CN 110133636 A CN110133636 A CN 110133636A
Authority
CN
China
Prior art keywords
track
gray area
correlation
gray
region
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201910426229.4A
Other languages
English (en)
Inventor
衣晓
周威
关欣
赵志勇
周正
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Naval Aeronautical University
Original Assignee
Naval Aeronautical University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Naval Aeronautical University filed Critical Naval Aeronautical University
Priority to CN201910426229.4A priority Critical patent/CN110133636A/zh
Publication of CN110133636A publication Critical patent/CN110133636A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/50Systems of measurement based on relative movement of target
    • G01S13/58Velocity or trajectory determination systems; Sense-of-movement determination systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/66Radar-tracking systems; Analogous systems
    • G01S13/72Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar
    • G01S13/723Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar by using numerical data
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/415Identification of targets based on measurements of movement associated with the target

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明针对系统误差下的航迹关联问题,公开了一种基于区域相关度的抗系统误差多传感器航迹关联方法。该方法分析了时变系统误差对航迹的影响,将其用灰区域来描述;定义了区域间的绝对距离以及区域测度的概念,通过分析区域间的相对位置关系,进而求得区域间的相关度,以此来实现航迹关联。该方法抗差性能好,步骤简单,耗时少,具有很强的适用性。

Description

一种基于区域相关度的抗差航迹关联方法
技术领域
本发明涉及一种基于区域相关度的抗差航迹关联方法。
背景技术
在分布式多目标多传感器跟踪系统中,航迹关联是一项关键技术,旨在判断来自不同传感器的航迹哪些来自同一目标,然后合并同源航迹,提高航迹的准确性。在现实中,雷达探测过程受到系统误差的影响,得到的目标航迹量测值与目标真实位置存在偏差,导致航迹关联的正确率大大降低。
对此,目前主要方法是先对航迹进行误差配准,之后再进行航迹关联。准确的航迹关联是以精准的误差配准为前提,而高精度的误差配准又需以准确的航迹关联为基础,它们两个互为前提。因此研究在带有系统误差的情况下如何进行航迹关联是很有必要的。目前已有的抗系统误差航迹关联方法,大多通过多次迭代计算,方法耗时较为严重;同时方法的适用条件一般都比较苛刻,适用性较差。因此迫切需要一种更好的方法解决大系统误差条件下航迹关联问题。
发明内容
为了克服现有技术的缺陷,本发明公开了一种基于区域相关度的抗差航迹关联方法。本发明分析了时变系统误差对航迹的影响,将其用灰区域来描述。定义了区域间的绝对距离以及区域测度的概念,通过分析区域间的相对位置关系,进而求得区域间的相关度,以此来实现航迹关联。
为了实现本发明的目的,本发明提供了一种基于区域相关度的抗差航迹关联方法。所述航迹关联方法,包括以下步骤:
步骤一,根据目标真实位置存在的扇形灰区域,计算不同雷达之间的灰区域中心、灰区域之间的绝对距离d1和灰区域测度;所述灰区域中心为扇形灰区域角平分线的中点,所述灰区域之间的绝对距离d1为不同灰区域中心之间的欧式距离,所述灰区域测度为所述灰区域内点到灰区域中心的欧式距离最大值,所述扇形灰区域由雷达量测值和雷达系统误差范围推算得到;
步骤二,根据所述灰区域中心、灰区域之间绝对距离d1、灰区域测度,计算灰区域之间的区域相关度其中,σa、σb分别为不同灰区域的灰区域测度;
步骤三,根据航迹之间的区域相关度,计算航迹间的区域灰关联度ζlj,进而实现航迹关联;航迹间区域灰关联度的计算方法为:其中,ρ为灰关联分辨系数,κij为雷达1的航迹i和雷达2的航迹j间的区域相关度。
与现有技术相比,本发明具有如下有益效果:
①抗系统误差性能好。本发明中的关联方法将系统误差对航迹的影响区域化描述,并且根据区域间的相对位置关系求得航迹间的相关度,对系统误差抗性好。
②适用性强。本发明中,方法在任何条件下均能保持很高的性能,同时方法步骤简单,耗时短,因此在实际应用中适用性更强。
附图说明
图1是系统误差灰区域示意图。
图2是区域间相对位置关系示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步详细描述。
为了解决现有技术中存在的技术问题,发明人分析了时变系统误差对航迹的影响,将其用灰区域来描述。定义了区域间的绝对距离以及区域测度的概念,通过分析区域间的相对位置关系,进而求得区域间的相关度,以此来实现航迹关联。
假设异地配置的两部2D雷达A(0,0)、B(xs,0)同时对目标区域内Ns批目标进行跟踪探测。同一处理周期内,雷达探测数目均为Ns,k时刻雷达A的探测航迹集合Γa(k)为:
式中,表示雷达A在k时刻对目标i探测得到的航迹量测值。
k时刻雷达B的探测航迹集合Γb(k)为:
式中,表示雷达B在k时刻对目标j探测得到的航迹量测值。
步骤一,推算目标真实位置存在区域:
雷达对目标进行观测时,由于受到时变系统误差以及随机误差的影响,雷达量测值与目标真实值相比存在偏移和旋转。以雷达A为例,雷达量测值即:
式中,为目标量测值,r、θ为目标真实值,Δra、Δθa为雷达测距、测角系统误差,其范围为(0,Δram)、(0,Δθam),为随机误差。
随机误差引起的目标位置变化较小,在此处忽略随机误差对航迹的影响。那么根据雷达量测值和时变系统误差的范围,可以推算得到目标真实位置存在的区域,即:
目标真实位置区域为一扇形区域,称其为系统误差灰区域,如图1所示。目标真实位置必定在该系统误差灰区域内。
步骤二,根据所述灰区域值,计算不同雷达之间的灰区域中心、灰区域之间的绝对距离和灰区域测度:
取系统误差灰区域中心,即:
定义两系统误差灰区域中心距离为区域间绝对距离,即:
式中,xoa、yoa为雷达A的系统误差灰区域中心,xob、yob为雷达B的系统误差灰区域中心。
当区域大小一定时,区域间绝对距离越大,两区域就越远;而当区域大小不固定时,绝对距离也就不能完全表示两区域位置关系上的远近。因此引入区域测度的概念,以此来描述区域的大小。
系统误差灰区域是一扇形区域,取区域内到区域中心点距离的最大值为测度。区域测度越大,区域也就越大。要求扇形区域内到区域中心的最远点,即求一个以区域中心位圆心的最小覆盖圆的半径,而对于覆盖圆来说,其弧度必定远远大于扇形边界的弧度。以雷达A为例,灰区域的测度即:
式中,A、B、C、D是灰区域四个端点,O是灰区域中心。
步骤三,计算区域相关度矩阵:
将雷达A探测得到的Ns条航迹看作已知航迹,将雷达B探测的Ns条航迹看作待识别航迹,取k时刻雷达A探测的Ns条航迹值与雷达B探测的第j条航迹值组成Ns+1行航迹关联决策矩阵,即:
式中,表示雷达A在k时刻探测的第i条航迹值所形成的系统误差灰区域的中心坐标;表示雷达A在k时刻探测的第i条航迹值所形成的系统误差灰区域的测度。
定义区域相关度为:
区域相关度能够表示两区域的位置关系。当两区域中心重合时,此时两区域相关度最高,为1;随着两目标区域间绝对距离的增大,区域相关度不断减小;同时,当区域间绝对距离一定时,随着区域测度的增大,区域相关度也就越大。
将Ns条已知航迹与待识别航迹j进行比较,运用上述方法,计算航迹之间的区域相关度,形成区域相关度矩阵Τ,即:
式中,κij为航迹i的系统误差灰区域和航迹j的系统误差灰区域之间的区域相关度。
步骤四,计算航迹间的区域灰关联度:
计算已知航迹l与待识别航迹j间的区域灰关联度ζlj,即:
式中,ρ为灰关联分辨系数,这里取0.5。
ζlj直接体现了已知航迹l与待识别航迹j直接的相似度,ζlj越大,两航迹来自同一目标的概率也就越大。航迹关联采取最大区域灰关联度原则,即若则认为已知航迹l与待识别航迹j关联。
实施例
假设两部异地配置的2D雷达的坐标为(0,0)、(100km,0)。目标初始位置在50km×50km的矩形区域内随机产生,目标做匀速直线运动,初始速度方向和大小在0~2πrad和200~400m/s随机产生。雷达采样周期为1s,目标批次为20。雷达A的测距、测角系统误差为0.5km、0.5rad,雷达B的测距、测角系统误差为-0.5km、-0.5rad。采用本专利方法计算相似度,对来自两部雷达的航迹进行关联。
采用本发明提出的关联方法,上述设计要求可按如下技术措施实施。
首先读取来自两雷达的数据:
其中,分别表示雷达A和雷达B的第i个航迹文件,雷达探测目标航迹数为20。
对雷达探测数据进行数据预处理,如图2所示,得到每个航迹值的系统误差灰区域中心坐标以及区域测度,即:
系统误差灰区域中心:
雷达A灰区域四个端点A、B、C、D的坐标:
Xa=Ra·Φa1 Ya=Ra·Φa2 (17)
式中,
雷达B灰区域四个端点E、F、G、H的坐标:
Xa=Ra·Φa1 Ya=Ra·Φa2 (18)
式中,
雷达A、B的灰区域测度为:
在第k个融合周期内,取雷达B的第20条航迹与雷达A的20条航迹的区域中心和区域测度组成关联决策矩阵,即:
区域相关度的定义:
式中,d1为两区域中心间绝对距离,σa、σb为两区域测度。
将第21行与前20行分别比较,计算它们之间的区域相关度,得到区域相关度矩阵T,即:
T=[0.136 0.157 0.048 0.081 0.067 0.063 0.039 0.148 0.212 0.151 0.1260.083 0.063 0.334 0.052 0.469] (22)
根据区域相关度矩阵计算航迹间的灰关联度ζlj,即:
根据上述方法计算得到关联度矩阵,即:
Z=[0.760 0.855 0.413 0.557 0.482 0.728 0.457 0.544 0.487 0.518 0.3560.687 0.759 0.692 0.710 0.524 0.476 0.818 0.427 1] (24)
航迹关联采取最大区域灰关联度原则,即若因此则认为来自雷达航迹B的第20条航迹与来自雷达航迹A的第20条航迹关联。

Claims (4)

1.一种基于区域相关度的抗差航迹关联方法,其特征在于,包括以下步骤:
步骤一,根据目标真实位置存在的扇形灰区域,计算不同雷达之间的灰区域中心、灰区域之间的绝对距离d1和灰区域测度;所述灰区域中心为扇形灰区域角平分线的中点,所述灰区域之间的绝对距离d1为不同灰区域中心之间的欧式距离,所述灰区域测度为所述灰区域内点到灰区域中心的欧式距离最大值;
步骤二,计算灰区域之间的区域相关度其中,σa、σb分别为不同灰区域的灰区域测度;
步骤三,根据航迹之间的区域相关度,计算航迹间的区域灰关联度ζlj,进而实现航迹关联。
2.根据权利要求1所述的基于区域相关度的抗差航迹关联方法,其特征在于,所述扇形灰区域由雷达量测值和雷达系统误差范围推算得到。
3.根据权利要求1或2所述的基于区域相关度的抗差航迹关联方法,其特征在于,步骤三中航迹间区域灰关联度的计算方法为:其中,ρ为灰关联分辨系数,κij为来自雷达1的航迹i和来自雷达2的航迹j间的区域相关度。
4.根据权利要求3所述的基于区域相关度的抗差航迹关联方法,其特征在于,ρ取0.5。
CN201910426229.4A 2019-05-21 2019-05-21 一种基于区域相关度的抗差航迹关联方法 Pending CN110133636A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910426229.4A CN110133636A (zh) 2019-05-21 2019-05-21 一种基于区域相关度的抗差航迹关联方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910426229.4A CN110133636A (zh) 2019-05-21 2019-05-21 一种基于区域相关度的抗差航迹关联方法

Publications (1)

Publication Number Publication Date
CN110133636A true CN110133636A (zh) 2019-08-16

Family

ID=67572372

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910426229.4A Pending CN110133636A (zh) 2019-05-21 2019-05-21 一种基于区域相关度的抗差航迹关联方法

Country Status (1)

Country Link
CN (1) CN110133636A (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5430445A (en) * 1992-12-31 1995-07-04 Raytheon Company Synthetic aperture radar guidance system and method of operating same
CN103983269A (zh) * 2014-01-26 2014-08-13 衣晓 一种基于区实混合序列的异步航迹关联算法
CN106443664A (zh) * 2016-09-13 2017-02-22 中国人民解放军海军航空工程学院 系统误差下基于拓扑信息的雷达与esm航迹关联方法
CN106646450A (zh) * 2016-12-29 2017-05-10 中国人民解放军海军航空工程学院 基于距离分步聚类的雷达航迹抗差关联方法
CN107728140A (zh) * 2017-11-22 2018-02-23 中国电子科技集团公司第二十八研究所 一种警戒雷达多目标多通道并行跟踪处理方法
CN108168564A (zh) * 2017-12-04 2018-06-15 上海无线电设备研究所 一种基于lhd灰色关联度的航迹关联方法
US20180289982A1 (en) * 2017-04-05 2018-10-11 Accuray Incorporated Sequential monoscopic tracking
CN109298413A (zh) * 2018-09-01 2019-02-01 哈尔滨工程大学 一种针对性解决复杂电磁环境下的多目标跟踪数据关联问题的方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5430445A (en) * 1992-12-31 1995-07-04 Raytheon Company Synthetic aperture radar guidance system and method of operating same
CN103983269A (zh) * 2014-01-26 2014-08-13 衣晓 一种基于区实混合序列的异步航迹关联算法
CN106443664A (zh) * 2016-09-13 2017-02-22 中国人民解放军海军航空工程学院 系统误差下基于拓扑信息的雷达与esm航迹关联方法
CN106646450A (zh) * 2016-12-29 2017-05-10 中国人民解放军海军航空工程学院 基于距离分步聚类的雷达航迹抗差关联方法
US20180289982A1 (en) * 2017-04-05 2018-10-11 Accuray Incorporated Sequential monoscopic tracking
CN107728140A (zh) * 2017-11-22 2018-02-23 中国电子科技集团公司第二十八研究所 一种警戒雷达多目标多通道并行跟踪处理方法
CN108168564A (zh) * 2017-12-04 2018-06-15 上海无线电设备研究所 一种基于lhd灰色关联度的航迹关联方法
CN109298413A (zh) * 2018-09-01 2019-02-01 哈尔滨工程大学 一种针对性解决复杂电磁环境下的多目标跟踪数据关联问题的方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Y. HUANG, Y. ZHOU AND H. ZHANG: "《Heterogeneous Sensors Track-to-Track Correlation Algorithm Based on Gray Correlative Degree》", 《2008 INTERNATIONAL SYMPOSIUM ON COMPUTER SCIENCE AND COMPUTATIONAL TECHNOLOGY》 *
衣晓; 张怀巍; 曹昕莹; 何友: "《基于区间灰数的分布式多目标航迹关联算法》", 《航空学报》 *
郑庆利: "《三种视角的灰色关联度性质、建模与应用研究》", 《中国优秀硕士学位论文全文数据库 (基础科学辑)》 *

Similar Documents

Publication Publication Date Title
CN106646450B (zh) 基于距离分步聚类的雷达航迹抗差关联方法
CN103729859B (zh) 一种基于模糊聚类的概率最近邻域多目标跟踪方法
CN103792532B (zh) 一种单脉冲高精度测角系统及其测角方法
CN106249232A (zh) 基于目标运动态势信息数据关联策略的目标跟踪方法
CN110132284B (zh) 一种基于深度信息的全局定位方法
CN103047982B (zh) 基于角度信息的目标自适应跟踪方法
CN109738902B (zh) 一种基于同步信标模式的水下高速目标高精度自主声学导航方法
CN110738275B (zh) 基于ut-phd的多传感器序贯融合跟踪方法
CN103592620B (zh) 长基线时差定位系统中高重频信号解定位模糊方法
CN109521420B (zh) 基于多特征匹配的多目标跟踪方法
CN109782290B (zh) 一种防止跟踪偏离的自动水声目标方位跟踪方法
CN110196409B (zh) 一种基于区域集合相对距离的抗差异步航迹关联方法
CN110133636A (zh) 一种基于区域相关度的抗差航迹关联方法
CN109490868B (zh) 一种基于分布式垂直线列阵的海上目标运动分析方法
CN115166785B (zh) 一种基于三接收机钟差单差的导航欺骗式干扰检测方法
CN103605980B (zh) 基于局部特征匹配的isar目标部件检测方法
CN116449298A (zh) 一种基于改进拓扑的航迹关联方法
CN102353958A (zh) 超短基线垂直运动目标测量方法
Manyika et al. A tracking sonar sensor for vehicle guidance
CN109084752B (zh) 一种基于全连通约束的地磁导航定位方法
CN111104985B (zh) 一种异步航迹关联的加权滑窗方法
Lou et al. Indoor localization and map building for autonomous mobile robot
CN111856412B (zh) 基于船摇补偿数据融合的抗欺骗式假目标干扰方法
CN111398896B (zh) 基于基准线成对交叉定位点距离的虚假点去除预处理方法
Yi et al. Research on underwater multi-clutter data association based on fuzzy clustering

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20190816

WD01 Invention patent application deemed withdrawn after publication