CN109581359B - 一种esm无源信息与雷达有源信息关联的方法 - Google Patents
一种esm无源信息与雷达有源信息关联的方法 Download PDFInfo
- Publication number
- CN109581359B CN109581359B CN201811607966.6A CN201811607966A CN109581359B CN 109581359 B CN109581359 B CN 109581359B CN 201811607966 A CN201811607966 A CN 201811607966A CN 109581359 B CN109581359 B CN 109581359B
- Authority
- CN
- China
- Prior art keywords
- esm
- radar
- track
- measurement
- time
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 15
- 238000005259 measurement Methods 0.000 claims abstract description 50
- 230000004927 fusion Effects 0.000 claims abstract description 13
- 239000008186 active pharmaceutical agent Substances 0.000 claims description 4
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000003786 synthesis reaction Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 abstract description 4
- 238000009825 accumulation Methods 0.000 abstract description 2
- 238000005094 computer simulation Methods 0.000 abstract 1
- 230000005855 radiation Effects 0.000 description 5
- 230000002596 correlated effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Systems 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/86—Combinations of radar systems with non-radar systems, e.g. sonar, direction finder
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
本发明属于信息融合技术领域,涉及一种ESM无源信息与雷达有源信息关联的方法。本发明首先使用线性内插法,将ESM数据对准到雷达测量时刻;然后使用1NN最近邻分类建立关联判别函数,并利用类似纽曼‑皮尔逊准则,建立判别门限;再使用Dempster‑Shafer证据理论融合不同时刻的信息;最后使用硬判决确定关联矩阵。计算机仿真表明本发明经过一定时间积累具有较好的性能。
Description
技术领域
本发明属于信息融合技术领域,涉及一种ESM无源信息与雷达有源信息关联的方法。
背景技术
在多传感器信息融合***中,多平台之间异类传感器信息融合一直是一个亟待解决的问题。在雷达-电子支援措施(Electronic Support Measures,ESM)协同工作的多平台***中,ESM通过对目标辐射源信号的分析获得辐射源的特征信息(发射频率,脉冲宽度等),且抗干扰能力强。雷达的目标定位能力强,且精度较高,但抗干扰能力和识别能力较差。通过雷达与ESM数据融合,一方面可以综合不同传感器的数据实现更可靠的决策;另一方面可以利用雷达的高精度方位数据,对目标威胁综合判定、综合态势生成提供有力的数据支持。
雷达-ESM数据融合首先应将两类传感器的数据进行关联,即要确定ESM数据源于哪一个雷达目标。由于异类传感器之间可以利用的信息不同,例如雷达可以获得辐射源的方位、速度和航向,而ESM传感器只能提供方位信息。因而雷达和ESM之间可以利用的信息只有角度,且两者的数据率一般不同,使得航迹关联存在很多不确定性。对于运动的目标辐射源,雷达-ESM航迹关联结果应实时更新,如何将之前时刻与当前时刻的关联结果进行融合,并且控制计算量在可以接受的范围内,也是目前需要解决的问题。
发明内容
本发明提供了一种ESM无源信息与雷达有源信息关联的方法。该方法首先利用模糊综合函数,制定关联准则,将每一时刻雷达与ESM数据之间的马氏距离(MahalanobisDistance)转化为航迹关联信任度,再利用DS证据理论,将该时刻的航迹关联信任度与之前时刻的航迹关联信任度融合,最后根据当前融合结果,利用最大似然的方法做出硬判决。
本发明采用的技术方案是:一种多平台数据实时关联的方法,包括以下步骤:
第一步:雷达和ESM时间对准
当雷达和ESM数据率不同时,首先应将观测数据进行时间对准。由于雷达数据率一般比较高(1次/秒),对于第j条航迹,两次ESM数据之间有nj个雷达测量值。本发明使用线性内插法,将ESM数据对准到雷达测量时刻。
第二步:建立关联判别函数
利用贝叶斯公式,式(3)可以表示为:
如果各个假设的先验概率相等,则式(5)可以表示为:
对(6)取对数后,式(5)可简化为
如果令:
其中要使得式子(7)左边最大也就是使得(8)最小。所以在先验概率相等的情况下,正态分布的贝叶斯分类被简化为了1NN最近邻分类。由于随着nj的增大,Dj也有增大趋势,无法直接将Dj作为雷达和ESM航迹关联的判别函数。所以令:
第三步:选取关联门限
(1)TL的选取
假设在方位上与雷达航迹j最接近的雷达航迹为k,门限TL决定的是本应该与雷达航迹j相关的ESM航迹,被错误地判为与雷达航迹k相关。令雷达航迹k与雷达航迹j在方位上间隔λ,则:Hk:εki~N(μ,σ2)和Hk:Dk~χ2(nj,λ)其中,门限TL的选择使得ESM测量值被错判为与方位上最邻近雷达航迹k的概率为Pe,则
(2)TH的选取
门限TH决定的是正确关联的雷达航迹被错误判为不相关的概率,即漏关联概率,设漏关联概率为Pm,则:Hj:εji~N(0,1)和Hj:Dj~χ2(nj)。Dj服从中心卡方分布,根据自由度为nj和漏关联概率Pm得到中心卡方分布的上侧分位数即
第四步:时间维度的信息融合
使用Dempster-Shafer证据理论是为了融合不同时刻的信息,这一理论可以理解为概率论的推广,即将概率(置信度)分配给一系列事件,而不是相互排斥的单个事件。进行时间维度的信息融合,应按照以下步骤进行。
(1)建立识别框架
在测量时刻t,雷达航迹有n个,Xi为第i条雷达航迹在ESM测量空间的投影。对于ESM测量e,建立辨识框架其中Ti表示“e与目标航迹Xi相关”,表示“e与目标航迹Xi不相关”,T0表示“e与任何目标航迹都不相关”,θ表示“不知道e的关联情况”。置信度是t时刻证据对于识别框架中各个元素的赋值,分别记为mt(Ti),mt(θ)。
(2)计算统计距离量
利用第二步的方法计算e与所有雷达航迹Xi之间的统计距离。
(3)证据产生
根据εi和门限TL、TH,为当前识别框架赋予信任度。
最后对当前时间的识别框架进行归一化。
(4)证据合成
得到了当前时刻识别框架对所有雷达航迹的信任度之后,使用DS证据理论将前一时刻的ESM测量信息与当前时刻的ESM测量信息进行融合。融合按照以下规则进行:
在t时刻,如果不需要做出硬判决,则等待ESM上报下一时刻的测量数据,返回第一步开始执行。
在t时刻,如果需要做出硬判决,则应按照以下规则执行:
找到mt(Ti)=maxmt(Tj),j=1,2,...,n,如果mt(Ti)>0.5,则ESM数据与雷达航迹Ti相关。
若m(T0)>0.5,则ESM数据与任何雷达航迹不相关。
若m(θ)>0.5,则ESM测量与雷达航迹的相关情况未知。
本发明的有益效果是:。
附图说明
图1:本发明一种ESM无源信息与雷达有源信息关联的方法的一种具体实施方式流程图;
图2:本发明具体实施中,目标运动方位示意图;
图3:本发明具体实施中,性能曲线示意图,其中(a)是关联错误率,(b)是漏关联率;
具体实施方式
下面结合具体场景详细说明本发明的技术方案,仿真场景设置如下:
假设我方ESM和雷达同地布置,坐标为(0,0)。敌方有四架飞机向我方飞来,每架飞机上搭载两个辐射源。相邻敌机之间的方位间隔为λ=1.5σ2。雷达只对前三个目标进行跟踪,测角均方误差为σr 2=0.1°,数据率1秒/次。ESM对所有目标进行角度测量,测角均方误差为σe 2=3°,数据率为4秒/次。雷达和ESM分别对目标的角度测量值如图2所示。
第一步:雷达和ESM时间对准
当雷达和ESM数据率不同时,首先将观测数据进行时间对准。对于第j条航迹,两次ESM数据之间有nj=4个雷达测量值。本发明使用线性内插法,将ESM数据对准到雷达测量时刻。
第二步:建立关联判别函数
第三步:选取关联门限
(1)TL的选取
(2)TH的选取
第四步:时间维度的信息融合
(1)建立识别框架
在测量时刻ti,雷达航迹有n个,Xi为第i条雷达航迹在ESM测量空间的投影。对于ESM测量e,建立辨识框架其中Ti表示“e与目标航迹Xi相关”,表示“e与目标航迹Xi不相关”,T0表示“e与任何目标航迹都不相关”,θ表示“不知道e的关联情况”。置信度是t时刻证据对于识别框架中各个元素的赋值,分别记为mt(Ti),mt(θ)。
(2)计算统计距离量
对于ESM测量e,利用第二步的方法计算e与所有雷达航迹Xi之间的统计距离。
(3)证据产生
最后对ti时刻的识别框架进行归一化。
(4)证据合成
得到了当前时刻识别框架对所有雷达航迹的信任度之后,使用DS证据理论将之前时刻的信息与当前时刻的信息进行融合。融合按照以下规则进行:
在ti时刻,如果需要做出硬判决,则应按照以下规则执行:
找到mt(Ti)=maxmt(Tj),j=1,2,...,n,如果mt(Ti)>0.5,则ESM数据与雷达航迹Ti相关;若m(T0)>0.5,则ESM数据与任何雷达航迹不相关;若m(θ)>0.5,则ESM测量与雷达航迹的相关情况未知。
本发明方法性能用以下两个数据定义:
1.错误率:本应该与雷达航迹Ti关联的ESM航迹,被错误判定为与雷达航迹Tj相关。
2.漏关联率:本应该与雷达航迹Ti关联的ESM航迹,被判决为不与任何雷达相关。
做500次蒙特卡洛仿真实验,仿真结果如图3所示,由图3(a)可以看出,经过约10个ESM测量时刻,目标的关联错误率下降到5%;由图3(b)可以看出,经过约3个ESM测量时刻,目标的漏关联概率下降到5%。说明经过一定时刻的积累,本发明可以达到预期关联指标,完成ESM测量数据和雷达测量数据的关联。
表4和表5分别展示了其中两个时刻的关联信任度矩阵。
表4t1时刻的关联可信度矩阵
表5t10时刻的关联可信度矩阵
对比表4和表5可以发现,经过10个时刻的迭代,原来无法判定是否与雷达航迹相关的ESM航迹,能够正确地与雷达航迹关联。例如t1时刻,无法判定ESM测量1与雷达航迹的关联情况。经过几次数据融合,在t10时刻EMS测量1与雷达航迹1的关联概率为0.999,与雷达航迹2的关联概率为0.001,几乎可以判定ESM测量1与雷达航迹1相关。
Claims (1)
1.一种ESM无源信息与雷达有源信息关联的方法,其特征在于,包括以下步骤:
第一步:雷达和ESM时间对准;
使用线性内插法,将ESM数据对准到雷达测量时刻:
设在tk时刻,ESM上传观测数据θ(tk),ESM前一次上报数据θ(tk-1)的时刻与当前时刻[tk-1,tk]之间,雷达有n个测量值,定义ESM和雷达的时间对准公式为:
第二步:建立关联判别函数:
第三步:建立关联门限
(1)建立低门限TL
假设在方位上与雷达航迹j最接近的雷达航迹为k,门限TL决定的是本应该与雷达航迹j相关的ESM航迹,被错误地判为与雷达航迹k相关;令雷达航迹k与雷达航迹j在方位上间隔其中μ为雷达航迹j和雷达航迹k在方位上的差值;
(2)建立高门限TH
第四步:时间维度的信息融合:
(1)建立识别框架
在测量时刻t,雷达航迹有n个,Xi为第i条雷达航迹在ESM测量空间的投影;对于ESM测量e,建立辨识框架其中Ti表示“e与Xi相关”,表示“e与Xi不相关”,T0表示“e与任何目标航迹都不相关”,θ表示“不知道e的关联情况”;
(2)计算统计距离
对于ESM测量e,利用第二步的方法计算e与所有雷达航迹Xi之间的统计距离;
(3)证据产生
根据εi和门限TL、TH,为当前识别框架赋予信任度,最后对当前时间的识别框架进行归一化;具体方法为:
根据εi和门限TL、TH,为当前识别框架赋予信任度
最后对当前时间的识别框架进行归一化
(4)证据合成
得到了当前时刻识别框架对所有雷达航迹的信任度之后,使用DS证据理论将前一时刻的ESM测量信息与当前时刻的ESM测量信息进行融合:
第五步:判决输出
在t时刻,如果不需要做出硬判决,则等待ESM上报下一时刻的测量数据,返回第一步开始执行;
在t时刻,如果需要做出硬判决,则应按照以下规则执行:
找到mt(Ti)=max mt(Tj),j=1,2,...,n,如果mt(Ti)>0.5,则ESM数据与雷达航迹Ti相关;
若mt(T0)>0.5,则ESM数据与任何雷达航迹不相关;
若mt(θ)>0.5,则ESM测量与雷达航迹的相关情况未知。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811607966.6A CN109581359B (zh) | 2018-12-27 | 2018-12-27 | 一种esm无源信息与雷达有源信息关联的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811607966.6A CN109581359B (zh) | 2018-12-27 | 2018-12-27 | 一种esm无源信息与雷达有源信息关联的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109581359A CN109581359A (zh) | 2019-04-05 |
CN109581359B true CN109581359B (zh) | 2023-04-28 |
Family
ID=65932027
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811607966.6A Active CN109581359B (zh) | 2018-12-27 | 2018-12-27 | 一种esm无源信息与雷达有源信息关联的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109581359B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110018450A (zh) * | 2019-04-30 | 2019-07-16 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | Ais与雷达角度***偏差的关联校准方法 |
CN110109095B (zh) * | 2019-04-30 | 2022-10-28 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 目标特征辅助多源数据的关联方法 |
CN110675418B (zh) * | 2019-09-26 | 2023-04-18 | 深圳市唯特视科技有限公司 | 一种基于ds证据理论的目标轨迹优化方法 |
CN111257826B (zh) * | 2020-01-15 | 2021-12-28 | 西安交通大学 | 一种多源异构传感器复合跟踪方法 |
CN111723474B (zh) * | 2020-06-01 | 2022-02-18 | 四川九洲电器集团有限责任公司 | 一种基于最大似然准则的临近目标识别方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106443665A (zh) * | 2016-09-13 | 2017-02-22 | 中国人民解放军海军航空工程学院 | 基于辅助距离的雷达与esm航迹关联方法 |
CN106443664A (zh) * | 2016-09-13 | 2017-02-22 | 中国人民解放军海军航空工程学院 | ***误差下基于拓扑信息的雷达与esm航迹关联方法 |
-
2018
- 2018-12-27 CN CN201811607966.6A patent/CN109581359B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106443665A (zh) * | 2016-09-13 | 2017-02-22 | 中国人民解放军海军航空工程学院 | 基于辅助距离的雷达与esm航迹关联方法 |
CN106443664A (zh) * | 2016-09-13 | 2017-02-22 | 中国人民解放军海军航空工程学院 | ***误差下基于拓扑信息的雷达与esm航迹关联方法 |
Non-Patent Citations (4)
Title |
---|
"D-S证据理论在ESM-雷达航迹关联中的应用";卜卿;《指挥信息***与技术》;20110831;第2卷(第4期);第60-63页 * |
"Emitter Target Recognition Based on Multi-sensor Data Fusion of ESM and IR";WANG Jiegui;《ICSP2008 Proceedings》;20081231;第1508-1511页 * |
"Modified combination rule of D-S evidence theory";LAN Xu-hui et al.;《2010 IEEE》;20101231;第1-3页 * |
"雷达与ESM模糊相关的双门限判决方法";王国宏;《电子对抗技术》;19970131(第1期);第15-20页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109581359A (zh) | 2019-04-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109581359B (zh) | 一种esm无源信息与雷达有源信息关联的方法 | |
CN109508000A (zh) | 异构多传感器多目标跟踪方法 | |
CN114859339B (zh) | 一种基于毫米波雷达的多目标跟踪方法 | |
Houenou et al. | A track-to-track association method for automotive perception systems | |
CN110109095B (zh) | 目标特征辅助多源数据的关联方法 | |
CN106443664B (zh) | ***误差下基于拓扑信息的雷达与esm航迹关联方法 | |
Ma et al. | Target tracking system for multi-sensor data fusion | |
US20060092074A1 (en) | Process for sensor resources management | |
CN104297739B (zh) | 一种通航监视中对光电跟踪设备的引导方法 | |
Abdallah et al. | A machine learning approach for localization in cellular environments | |
CN110058205A (zh) | 一种基于迭代最近点算法的警戒雷达***误差校正方法 | |
CN113075648B (zh) | 一种无人集群目标定位信息的聚类与滤波方法 | |
CN111259332B (zh) | 一种杂波环境下的模糊数据关联方法及多目标跟踪方法 | |
CN108596117B (zh) | 一种基于二维激光测距仪阵列的场景监控方法 | |
CN110907903B (zh) | 一种基于多传感器的自适应跟踪处理方法 | |
KR20170103260A (ko) | 신호의 도래 시간차 정보 매칭을 이용한 위치 탐지 방법 | |
CN111487612A (zh) | 基于cpd的异地配置雷达/esm航迹抗差关联方法 | |
CN113687429B (zh) | 一种确定毫米波雷达监测区域边界的装置及方法 | |
CN113988228B (zh) | 一种基于rfid和视觉融合的室内监控方法及*** | |
CN112990125B (zh) | 一种辐射源雷达属于目标平台的判别方法 | |
Si et al. | Robust Poisson multi-Bernoulli filter with unknown clutter rate | |
RU2562616C1 (ru) | Способ получения радиотехнической информации и радиотехнический комплекс для его осуществления | |
CN117310642A (zh) | 一种基于多尺度聚类的多雷达密集群目标航迹关联方法 | |
CN111487613A (zh) | 基于分级聚类的雷达/esm航迹抗差关联方法 | |
CN116953692A (zh) | 一种主动雷达和被动雷达协同跟踪下的航迹关联方法 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20240319 Address after: No. 351 Tiyu Road, Xiaodian District, Taiyuan City, Shanxi Province 030000 Patentee after: NORTH AUTOMATIC CONTROL TECHNOLOGY INSTITUTE Country or region after: Zhong Guo Address before: 611731, No. 2006, West Avenue, hi tech West District, Sichuan, Chengdu Patentee before: University of Electronic Science and Technology of China Country or region before: Zhong Guo |