CN103886196B - 一种抗大***误差的航迹关联方法 - Google Patents
一种抗大***误差的航迹关联方法 Download PDFInfo
- Publication number
- CN103886196B CN103886196B CN201410095990.1A CN201410095990A CN103886196B CN 103886196 B CN103886196 B CN 103886196B CN 201410095990 A CN201410095990 A CN 201410095990A CN 103886196 B CN103886196 B CN 103886196B
- Authority
- CN
- China
- Prior art keywords
- radar
- flight path
- delta
- association
- systematic error
- 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.)
- Expired - Fee Related
Links
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种抗大***误差的雷达组网航迹关联方法。该方法首先以雷达某一待关联航迹为中心,在其可能的***误差范围内进行航迹预关联假设,得到多个预关联关系,然后根据各个预关联关系,利用航迹的偏差信息估计出相应的***误差,并根据***误差对该雷达所有航迹进行修正,然后进行全局匹配相似度计算,最后根据最优匹配关系来确定目标的关联关系。本发明的优点在于航迹关联前不需要对雷达***误差进行校正,在大***误差下能够较好解决虚增批与错误关联,这是传统概率统计关联算法不能做到的。该方法能够实现大***误差下航迹的正确关联,可应用于雷达组网的信息融合处理***。
Description
技术领域
本发明属于雷达组网航迹关联技术领域,特别是一种抗大***误差的航迹关联方法。
背景技术
雷达组网的航迹关联,一般是建立在***误差配准的基础上,然后采用基于统计理论的方法,有加权法、独立序贯法、经典分配法、独立双门限法、修正最近邻法和基于模糊数学的方法,通常认为关联门限阈值的大小直接依赖于网内各雷达的状态估计误差协方差,并依据经验将关联门限放大,使其能够容忍配准后仍存在的***误差残差,但如果雷达的***误差未配准,或者配准后残差过大,采用基于统计理论的方法关联,将会出现大量增批。
发明内容
发明目的:本发明所要解决的技术问题是针对现有技术的不足,提供一种抗大***误差的航迹关联方法。
为了解决上述技术问题,本发明公开了一种抗大***误差的航迹关联方法,包括如下步骤:
步骤(1)以待关联的雷达航迹为中心,根据可能的最大***误差,计算其关联区域范围;
步骤(2)计算处于该区域范围内的所有已存在的***航迹,并假设该雷达航迹与处于该区域***航迹一一关联,形成预关联假设集;
步骤(3)根据预关联假设,计算该雷达航迹与预关联***航迹的航迹偏差;
步骤(4)根据航迹偏差,进一步计算该雷达的***误差;
步骤(5)根据计算的***误差,对该雷达所有航迹进行修正;
步骤(6)将修正后航迹与***航迹进行全局匹配相似度计算;
步骤(7)根据全局匹配相似度,确定关联关系。
本发明中,所述大***误差,可以是大的测距***误差和大的测向***误差。
本发明步骤(1)中的关联区域范围R(x,y)是根据最大***误差计算得到,满足:
其中(x,y)为关联区域范围的二维平面直角坐标,为当前待关联雷达航迹Ti的最新位置坐标,(x0,y0)为该雷达站基点位置坐标,Δθmax为最大的测向***误差,Δρmax为最大的测距***误差。为当前待关联雷达航迹Ti的最新位置,(x0,y0)为该雷达站基点,Δθmax,Δρmax分别为最大的测向***误差和测距***误差。
本发明步骤(2)中,假定该雷达航迹与处于该区域***航迹一一关联,形成预关联假设集表述如下:
Hj:Ti→TSj,j=1,2…n
其中Hj为第j个假设,Ti表示第i条雷达航迹,处于关联区域的***航迹总数为n条,TSj为第j条***航迹,符号→表示其左右两航迹关联。
本发明步骤(3)中的航迹偏差可采用多周期点迹的均值统计偏差,描述如下:
其中(Δxj,Δyj)分别表示第j个假设下的航迹偏差,m表示多周期航迹点数,表示航迹Ti的第k个坐标点,表示***航迹j的第k个坐标。
本发明步骤(4)中的***误差也将形成多个假设值,计算如下:
其中为原航迹Ti位置点(xT,yT)对应的方位与距离。Δθj和Δρj分别为假设Hj下的雷达***方位误差和距离误差。
本发明步骤(5)中的航迹修正要对该雷达所有航迹的当前位置点进行修正,描述如下:
其中l为该雷达航迹的总数。
本发明步骤(6)中的全局匹配相似度表示该雷达所有航迹与***航迹匹配的相似度总和,相似度的计算可由航迹坐标、速度因素构成:
ζbest=min(ζj),j=1,2…m
其中为航迹Ti的x-y平面速度分量,是***航迹TSj的x-y平面速度分量,wpos是相似度中位置项部分的权重,wv是相似度中速度项部分的权重,本发明中设置位置项权重wpos0.8,速度项权重0.2,ζi,j是雷达航迹Ti在假设Hj下的最优匹配相似度,ζj是该雷达所有航迹在假设Hj的全局匹配相似度,ζbest是所有假设中最优匹配相似度,即最优假设。
本发明步骤(7)确定关联关系可采用双门限(C1,C2,满足C1<C2)判定方法,若ζbest<C1,best假设成立,判定航迹Ti与Tbest为关联,C1≤ζbest≤C2为模糊,C2≤ζbest为不关联。
本发明针对大***误差情况采用了假设关联,全局匹配印证的关联方式,能够有效实现大***误差下的航迹正确关联。
有益效果:本发明的显著优点在于1)航迹关联前不需要进行***误差配准;2)在***误差较大,或者配准后***误差残差较大的情况下,仍然能够保持正确关联,具有较好的工程应用前景。
附图说明
下面结合附图和具体实施方式对本发明做更进一步的具体说明,本发明的上述和/或其他方面的优点将会变得更加清楚。
图1是本发明一种抗大***误差的航迹关联方法的信息流程图。
具体实施方式
结合图1,本发明一种雷达组网***误差配准方法,步骤如下:
步骤(1)以待关联的雷达航迹为中心,根据可能的最大***误差,计算其关联区域R(x,y),满足:
其中(x,y)为二维平面直角坐标,为当前待关联雷达航迹Ti的最新位置,(x0,y0)为该雷达站基点,Δθmax,Δρmax分别为最大的测向***误差和测距***误差。
步骤(2)计算处于该区域范围内的所有已存在的***航迹,并假设该雷达航迹与处于该区域***航迹一一关联,形成预关联假设集:
Hj:Ti→TSj,j=1,2…n
其中Hj为第j个假设,处于关联区域的***航迹总数为n条,TSj为第j条***航迹,符号→表示其左右两航迹关联。
步骤(3)根据假定的关联关系,计算该雷达航迹与***航迹的多周期航迹偏差,采用多周期点迹的均值统计偏差,描述如下:
其中(Δxj,Δyj)分别表示第j个假设下的航迹偏差,m表示多周期航迹点数,表示航迹Ti的第k个坐标点,表示***航迹j的第k个坐标。
步骤(4)根据航迹偏差,进一步计算各假设下该雷达的***误差:
其中为原航迹Ti位置点(xT,yT)对应的方位与距离。
步骤(5)根据计算的***误差,对各假设下该雷达所有航迹进行修正:
其中l为该雷达航迹的总数。
步骤(6)将修正后航迹与***航迹进行全局匹配相似度计算,全局匹配相似度表示该雷达所有航迹与***航迹匹配的相似度总和,相似度的计算可由航迹坐标、速度因素构成,如下:
ζbest=min(ζj),j=1,2…m
其中为航迹Ti的x-y平面速度分量,是***航迹TSj的x-y平面速度分量,ζi,j是航迹Ti在假设Hj下的最优匹配相似度,ζj是该雷达所有航迹在假设Hj的全局匹配相似度,ζbest是所有假设中最优匹配相似度,也就是最优假设。
步骤(7)根据全局匹配相似度,可采用双门限(C1,C2,满足C1<C2)判定方法确定其关联关系,若ζbest<C1,best假设成立,判定航迹Ti与Tbest为关联,C1≤ζbest≤C2为模糊,等待下一周期再判别,C2≤ζbest为不关联,为新航迹。
实施例1
下面结合附图1,说明本发明一个实施例:
本实施例对雷达组网中某部存在大***误差雷达的航迹进行关联处理,该雷达***方位***误差为Δθ=3°,距离***误差Δρ=0.3km,站基点坐标为(120.345km,33.249km),最大可能方位***误差Δθmax=5°,Δρmax=1km。该雷达共探测到4个目标(T1,T2,T3,T4),多雷达融合***已通过组网的其它雷达形成了本区域的***航迹,与之对应的4个目标分别为(TS1,TS2,TS3,TS4),某周期时空对准后的各航迹最新位置点及其速度如表1所示。由于***误差的存在,T1,T2,T3,T4与对应的TS1,TS2,TS3,TS4航迹之间有较大的偏差,而此时T1与TS3反而更为接近,传统的概率统计关联方法容易造成增批或错误关联。
表1航迹数据
下面以T1为当前需判关联航迹进行说明。
步骤(1)以T1为中心,根据Δθmax=5°,Δρmax=1km,计算其关联区域R(x,y)范围,满足:
即:
步骤(2)计算处于该区域范围内的所有已存在的***航迹,T1与TS1,TS2,TS3,TS4的距离分别为:10.166、170.8612、4.2447和54.5290,那么处于该区域范围内的***航迹为TS1和TS3,假设T1分别与TS1和TS3一一关联,形成预关联假设集:
H1:T1→TS1
H2:T1→TS3
步骤(3)根据假定的关联关系,计算该雷达航迹与***航迹的航迹偏差:
步骤(4)根据航迹偏差,进一步计算各假设下该雷达的***误差:
步骤(5)根据计算的***误差,对各假设下该雷达所有航迹进行修正,假设1下的航迹修正为:
即:
假设2下的航迹修正为:
即:
步骤(6)将修正后航迹与***航迹进行全局匹配相似度计算,取wpos=0.8,wv=0.2,在假设1、假设2下分别计算修正后航迹(T1,T2,T3,T4)的最优匹配相似度。
假设1情况下,T1与TS1假定匹配,剩余航迹T2,T3,T4分别与TS2,TS3,TS4最匹配,相似度为:
那么:
假设2情况下,T1与TS3假定匹配,剩余航迹T2,T3,T4分别与TS2,TS1,TS4最匹配,相似度为:
那么:
最优假设的相似度为:
ζbest=min(ζ1,ζ2)=ζ1=0.91
步骤(7)根据全局匹配相似度,采用双门限判定方法确定其关联关系,令C1=5(表示5km偏差范围内),C2=10(表示10km偏差范围内),当前最优假设的相似度满足ζbest<C1,即假设H1成立,判定T1与TS1关联。
本发明提供了一种抗大***误差的航迹关联方法,具体实现该技术方案的方法和途径很多,以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。本实施例中未明确的各组成部分均可用现有技术加以实现。
Claims (8)
1.一种抗大***误差的航迹关联方法,其特征在于,包括如下步骤:
步骤(1)以待关联的雷达航迹为中心,根据最大***误差,计算雷达的关联区域范围;
步骤(2)计算处于所述关联区域范围内的所有已存在的***航迹,并假设该雷达航迹与处于该区域***航迹一一关联,形成预关联假设集;
步骤(3)根据预关联假设,计算雷达航迹与预关联***航迹的航迹偏差;
步骤(4)根据航迹偏差,进一步计算该雷达***误差;
步骤(5)根据计算的***误差,对该雷达所有航迹进行修正;
步骤(6)将修正后航迹与***航迹进行全局匹配相似度计算;
步骤(7)根据全局匹配相似度,确定关联关系;
步骤(1)中的关联区域范围R(x,y)描述如下:
其中(x,y)为关联区域范围的二维平面直角坐标,为当前待关联雷达航迹Ti的最新位置坐标,(x0,y0)为雷达站基点位置坐标,Δθmax为最大的测向***误差,Δρmax为最大的测距***误差。
2.根据权利要求1所述的抗大***误差的航迹关联方法,其特征在于,所述大***误差为测距***误差和测向***误差。
3.根据权利要求1所述的抗大***误差的航迹关联方法,其特征在于,步骤(2)中,假设该雷达航迹与处于该区域***航迹一一关联,形成预关联假设集表示如下:
Hj:Ti→TSj,j=1,2…n,
其中Hj为第j个假设,Ti表示第i条雷达航迹,处于关联区域的***航迹总数为n条,TSj为第j条***航迹,符号→表示其两个航迹关联。
4.根据权利要求3所述的抗大***误差的航迹关联方法,其特征在于,步骤(3)中的航迹偏差采用多周期点迹的均值统计偏差,公式如下:
其中(Δxj,Δyj)分别表示第j个假设下的航迹偏差坐标,m表示多周期航迹点数,表示雷达航迹Ti的第k个坐标点,表示***航迹j的第k个坐标。
5.根据权利要求4所述的抗大***误差的航迹关联方法,其特征在于,步骤(4)中的***误差包括1个以上假设值,计算如下:
其中为原雷达航迹Ti位置点(xT,yT)对应的方位与距离,Δθj和Δρj分别为假设Hj下的雷达***方位误差和距离误差。
6.根据权利要求5所述的抗大***误差的航迹关联方法,其特征在于,步骤(5)中的航迹修正要对该雷达所有航迹的当前位置点进行修正,公式如下:
其中l为该雷达航迹的总数。
7.根据权利要求6所述的抗大***误差的航迹关联方法,其特征在于,步骤(6)中的全局匹配相似度表示该雷达所有航迹与***航迹匹配的相似度总和,相似度的计算公式为:
ζbest=min(ζj),j=1,2…n,
其中为雷达航迹Ti的x-y平面速度分量,是***航迹TSj的x-y平面速度分量,wpos是相似度中位置项部分的权重,wυ是相似度中速度项部分的权重,ζi,j是雷达航迹Ti在假设Hj下的最优匹配相似度,ζj是该雷达所有航迹在假设Hj的全局匹配相似度,ζbest是所有假设中最优匹配相似度,即最优假设。
8.根据权利要求7所述的抗大***误差的航迹关联方法,其特征在于,步骤(7)确定关联关系采用双门限C1和C2判定方法,且C1<C2,若ζbest<C1,best假设成立,判定航迹Ti与Tbest为关联,若C1≤ζbest≤C2判定为模糊,若C2≤ζbest判定为不关联。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410095990.1A CN103886196B (zh) | 2014-03-14 | 2014-03-14 | 一种抗大***误差的航迹关联方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410095990.1A CN103886196B (zh) | 2014-03-14 | 2014-03-14 | 一种抗大***误差的航迹关联方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103886196A CN103886196A (zh) | 2014-06-25 |
CN103886196B true CN103886196B (zh) | 2016-08-24 |
Family
ID=50955086
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410095990.1A Expired - Fee Related CN103886196B (zh) | 2014-03-14 | 2014-03-14 | 一种抗大***误差的航迹关联方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103886196B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104268865B (zh) * | 2014-09-18 | 2017-06-16 | 中国人民解放军海军航空工程学院 | 基于模板匹配的集中式多传感器编队目标跟踪方法 |
CN104318094B (zh) * | 2014-10-15 | 2017-04-05 | 中国电子科技集团公司第二十八研究所 | 一种编队飞行的雷达组网航迹关联方法 |
CN105510896B (zh) * | 2015-12-03 | 2019-07-02 | 成都能通科技有限公司 | 一种集中式多雷达数据处理的加权最近邻域数据关联方法 |
CN106934324A (zh) * | 2015-12-30 | 2017-07-07 | 南京理工大学 | 基于简化多假设算法的雷达数据关联方法 |
CN107066806B (zh) * | 2017-02-15 | 2019-07-26 | 中国人民解放军海军航空工程学院 | 航迹关联方法及装置 |
CN107063259B (zh) * | 2017-03-08 | 2020-06-09 | 四川九洲电器集团有限责任公司 | 一种航迹关联方法及电子设备 |
CN107271974B (zh) * | 2017-06-08 | 2020-10-20 | 中国人民解放军海军航空大学 | 一种基于稳定角点的时空误差求取方法 |
CN108490432B (zh) * | 2018-05-18 | 2020-03-27 | 中船重工鹏力(南京)大气海洋信息***有限公司 | 一种导航雷达***中全局目标跟踪的方法 |
CN109856649A (zh) * | 2019-02-28 | 2019-06-07 | 百度在线网络技术(北京)有限公司 | 导航定位的误差分析方法、装置及存储介质 |
CN111157969B (zh) * | 2020-01-20 | 2022-01-25 | 四川九洲电器集团有限责任公司 | 一种多传感器***误差在线估计和数据融合方法 |
CN113124871B (zh) * | 2021-03-17 | 2023-05-02 | 北京理工大学 | 一种基于数据质量评估的自适应航迹关联方法 |
CN116520307B (zh) * | 2023-02-20 | 2023-10-27 | 扬州宇安电子科技有限公司 | 一种基于目标追踪的雷达协同检测***及方法 |
CN116380148B (zh) * | 2023-04-06 | 2023-11-10 | 中国人民解放军93209部队 | 多传感器目标跟踪***的两级时空误差标校方法及装置 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102194332A (zh) * | 2011-03-24 | 2011-09-21 | 中国船舶重工集团公司第七○九研究所 | 一种自适应航迹数据相关方法 |
-
2014
- 2014-03-14 CN CN201410095990.1A patent/CN103886196B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102194332A (zh) * | 2011-03-24 | 2011-09-21 | 中国船舶重工集团公司第七○九研究所 | 一种自适应航迹数据相关方法 |
Non-Patent Citations (3)
Title |
---|
《一种抗差型航迹关联算法研究》;沈蕾 等;《舰船电子工程》;20081031;第28卷(第10期);135-138 * |
《基于多假设的航迹关联方法研究》;文云峰 等;《舰船科学技术》;20110228;第33卷(第2期);116-120 * |
《雷达情报数据融合***的误差校正和航迹关联技术研究》;陈列;《中国优秀硕士学位论文全文数据库-信息科技辑》;20080215;第2008年卷(第2期);摘要,第3.3、3.4.1、3.4.5、3.5、4.4.2节 * |
Also Published As
Publication number | Publication date |
---|---|
CN103886196A (zh) | 2014-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103886196B (zh) | 一种抗大***误差的航迹关联方法 | |
CN108536851B (zh) | 一种基于移动轨迹相似度比较的用户身份识别方法 | |
CN106793087B (zh) | 一种基于aoa和pdoa的阵列天线室内定位方法 | |
CN107909600A (zh) | 一种基于视觉的无人机实时运动目标分类与检测方法 | |
CN106055091A (zh) | 一种基于深度信息和校正方式的手部姿态估计方法 | |
CN108182413A (zh) | 一种矿井移动目标检测与跟踪识别方法 | |
CN105828289A (zh) | 一种基于信道状态信息的无源室内定位方法 | |
CN104602342A (zh) | 一种基于iBeacon设备的高效室内定位方法 | |
CN104050368B (zh) | ***误差下基于误差补偿的群航迹精细关联算法 | |
CN105954712A (zh) | 联合无线电信号复包络和载波相位信息的多目标直接定位方法 | |
CN104268866A (zh) | 基于运动信息与背景信息相结合的视频序列配准方法 | |
CN110488222B (zh) | 一种nlos条件下svm与重心坐标相结合的uwb定位方法 | |
CN103135117A (zh) | 一种分散式多机器人协同定位算法 | |
CN104331078B (zh) | 一种基于位置映射算法的多机器人协作定位方法 | |
CN106842121A (zh) | 视距与非视距混合环境中基于到达时间差的稳健定位方法 | |
CN105913452A (zh) | 空间碎片实时检测与跟踪方法 | |
CN112929826A (zh) | 一种基于uwb三角定位与指纹信息的室内混合定位方法 | |
CN104573713A (zh) | 一种基于图像纹理特征的互感器红外图像识别方法 | |
CN103792515B (zh) | 一种异平台2维雷达与红外传感器量测数据合成方法 | |
CN106154221A (zh) | 一种基于wlan的半监督定位方法 | |
CN105096341A (zh) | 基于三焦张量和关键帧策略的移动机器人位姿估计方法 | |
CN108737952A (zh) | 基于rssi测距改进的多边形加权质心定位方法 | |
CN110412504A (zh) | 基于角度与时差信息的无源多站多目标关联与定位方法 | |
CN103491627A (zh) | 一种集成多种算法的近距离实时精确定位方法 | |
CN106295395A (zh) | 基于图划分的不确定轨迹隐私保护方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160824 |