CN1306280C - 双平台多辐射源测向测时交叉定位方法 - Google Patents
双平台多辐射源测向测时交叉定位方法 Download PDFInfo
- Publication number
- CN1306280C CN1306280C CNB031292720A CN03129272A CN1306280C CN 1306280 C CN1306280 C CN 1306280C CN B031292720 A CNB031292720 A CN B031292720A CN 03129272 A CN03129272 A CN 03129272A CN 1306280 C CN1306280 C CN 1306280C
- Authority
- CN
- China
- Prior art keywords
- platform
- radiation source
- data
- measurement data
- measurement
- 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
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明涉及一种双平台多辐射源测向测时交叉定位方法,主从平台独立采集数据,每组数据包括辐射源的方位角以及对应的时间标记,N个辐射源形成N*N个交叉点,两个平台形成N*N组测量数据和交叉点相对应。每个扫描周期结束从平台通过通信传输将数据传送到主平台进行融合,通过计算双平台间测量数据的相近度来进行最优匹配,去除虚假定位。本发明无需辐射源的先验信息,降低了平台间数据传输、处理和时间同步的复杂性,运算速度快,定位精度高,可广泛应用于航海、航空、宇航、侦察、测控、救援和地球物理学等各类军、民用系统。
Description
技术领域:
本发明涉及一种纯角度无源定位时双平台多辐射源测向测时交叉定位方法,是红外搜索与跟踪系统、无源定位系统、红外预警系统、大视场目标监视系统、空中交通管制系统等的一项核心技术,在各类军、民用系统中均可有广泛的应用。
背景技术:
对辐射源的准确定位,在诸多军用和民用系统中具有极其重要的意义,对军事系统来说,它有助于精确打击武器的使用,为最终摧毁消灭对方提供有力的保障;就民用系统来说,可以为辐射源提供可靠的服务,起到安全保障作用。在航海、航空、宇航、侦察、测控、救援和地球物理学中扮演着重要的角色。
有源定位是利用雷达、激光、声纳等有源设备来对辐射源进行定位,它具有全天候、高精度等优点。然而有源定位的使用,要发射大功率的信号,这样很容易暴露自己,被对方发现从而遭受破坏,使定位精度受到很大影响。通过对辐射源信号的截获、测量获得目标的位置和航迹,实现定位的方法,称为无源定位。它不发射自己的信号,具有作用距离远,隐蔽性好等优点。因而具有极强的生存能力和反隐身能力。和有源定位系统相比具有低耗费、低维修、对环境友好等优点。随着信号截获和处理技术的发展,无源定位技术的研究和应用越来越广泛。
在无源定位方法中,由于方向测量定位方法(测向交叉定位方法)具有方向角度变化慢,对平台间的时间统一要求低等特点,使得方向测量定位方法成为发展最早、研究最多和使用最广泛的无源定位技术。它是通过高精度的测向设备,在两个或两个以上的观测平台对辐射源测向,各个测向线(视线)的交叉点就是辐射源的地理位置。根据在各观测平台测得的辐射源的方向和各观测平台的位置,通过三角运算就可以确定交叉点的坐标。然而在多辐射源的情况下,此方法的主要缺点是存在虚假定位。对于空间存在N个辐射源,两个平台共有2N条方向线,存在N*N个交叉点,其中有N*(N-1)个交叉点是虚假定位。因此从实际应用的角度出发,要求去除这些虚假定位。
国内外研究人员针对虚假点的去除提出了一些方法。对于多平台(平台数≥3)无源定位,已经证明其平台间数据关联具有NP计算复杂性难题。目前研究人员主要应用图论、运筹学、组合数学等工具求解该问题,消除虚假定位的方法主要有最小距离法、最大似然法、谱相关法、整数规划法、拉格朗日松弛算法、神经网络等,但是平台过多导致平台间的数据传输、处理和时间的同步都变得十分复杂,限制了其应用。一般情况下,很难依据两个平台来去除虚假点,目前双平台无源定位主要是结合辐射源的先验信息或者冗余信息去除虚假点。信号处理的方法要求知道信号的一些先验信息,通过信号的谱相关理论、盲分离算法等进行定位,平台接收装置必须采用均匀阵列,算法复杂。
发明内容:
本发明的目的在于针对现有技术的上述不足,提出一种双平台多辐射源测向测时交叉定位方法,可以在两平台的情况下,无需先验信息即可根据平台间测量数据的相近度来进行最优匹配,去除虚假定位,满足国内实际系统的性能需求。
为实现这样的目的,本发明的技术方案中采用双平台进行定位,主从平台独立采集数据。每个扫描周期结束从平台通过通信传输将数据传送到主平台进行融合。在主平台融合中心进行数据处理,计算双平台间测量数据的相近度来进行最优匹配,去除虚假定位。对于N个辐射源,每个平台测量出N组数据,每组数据包括辐射源的方位角以及对应的时间标记。N个辐射源可以形成N*N个交叉点,两个平台形成N*N组测量数据和交叉点相对应。在一组数据中,根据两个方位角可以得到此交叉点和主平台的斜距,另外,通过测量时间差和主平台的测量方位角也可以得到此交叉点和主平台的斜距。根据两个斜距的误差程度,可以得到这组测量数据的相近度。根据相近度的数值可以完成虚假点的去除。
本发明的方法包括如下具体步骤::
1.测量数据采集:主平台从平台是两个独立的系统,每个平台都可以独立测量监视区域的辐射源。平台携带的红外搜索与跟踪系统(IRST)可以测量目标的方位角,测量精度可以达到3毫弧;平台采用基于全球定位系统GPS的统一时间基准。GPS的特点是可以为全球各地随时提供一个高精度的时钟,从而也就保证了各平台信号的同步。由GPS接收板、高稳定度晶体震荡器和时钟板组成的GPS钟可以提供高精度的时间标准,时间精度可以达到20纳秒。一个扫描周期结束,从平台通过通信传输将数据传送到主平台进行后续处理。
2.测量数据处理:融合主从平台的测量数据,对于N个辐射源,形成N*N组测量数据(即N*N个交叉点)。每组测量数据包括三个元素:(1)主平台的测量角度;(2)从平台的测量角度;(3)主从平台的测量时差。根据简单三角公式,由(1)(2)可以计算出辐射源到主平台的斜距。根据测角测时差定位公式,由(1)(3)可以计算出辐射源到主平台的斜距。如果不考虑测量误差的情况下,如果交叉点是辐射源,显然两个斜距应该相等。然而实际情况测量噪声肯定存在的,因此对于真实辐射源,两个斜距数值也不相等。根据斜距的误差,定义主从平台的测量数据的相近度P(P∈[0,1])。P的数值越大,表明平台间测量数据来自同一辐射源的可能性越大,即此组数据代表真实辐射源。重复上述过程,得到一个N*N的相近度矩阵
j1表示主平台的第j1个测量数据,j2表示从平台的第j2个测量数据。通过测量数据的处理,最终得到相近度矩阵A。
3.虚假点去除:因为平台的每个测量角度代表一个辐射源,不同的辐射源有不同的角度,因此从A中选出最大的N个不同行不同列元素,则此N个元素对应的测量数据就代表辐射源,其余都为虚假点。从而完成虚假点去除。
本发明可用于两个测量平台对多个辐射源进行定位,通过引入时刻量测量,提出相近度的概念,算法简单有效、实时性强,无需辐射源的先验信息,降低了平台间数据传输、处理和时间同步的复杂性,运算速度快,正确定位精度高,能够有效去除无源定位中的虚假点,可广泛应用于航海、航空、宇航、侦察、测控、救援和地球物理学等各类军、民用系统,具有广阔的市场前景和应用价值。
附图说明:
图1为测向交叉定位示意图。
如图1所示,E为辐射源,S1,S2为测量平台。已知测量平台的位置,测量的角度,以X轴为基准。
图2为多辐射源情况下测向交叉定位产生虚假点示意图。
如图2所示,监视区域存在3个辐射源,共产生9个交叉点,其中3个交叉点是真实辐射源,其余6个是虚假点。
图3为本发明实施例中辐射源水平编队态势图。
图4为本发明实施例中辐射源十字编队态势图。
具体实施方式:
为了更好地理解本发明的技术方案,以下结合附图对本发明的实施方式作进一步描述。为方便说明问题,在此作如下假设,这些假设并不影响算法实质,并且本发明很容易推广到三维实际情况中去。
a.观测平台和辐射源都在同一平面内,如XY平面内;
b.每个测量平台都能正常接收到所有辐射源辐射的信号,无杂波干扰;
c.信号没有重叠情况,即有N个目标,可以测到N组数据;
d.双平台接收数据具有相同的测量精度,并且已经时空配准。
1.测量数据采集
在这些假设前提下,每个测量站分别测得N个辐射源的方位角α以及对应的时间标记T。本系统采用基于全球定位系统GPS的统一时间基准。GPS的特点是可以为全球各地随时提供一个高精度的时钟,从而也就保证了各平台信号的同步。由GPS接收板、高稳定度晶体震荡器和时钟板组成的GPS钟可以提供高精度的时间标准。因此定义i平台对j辐射源的观测向量为:
(i=1,2;j=1,2....N)。平台1为主平台,2为从平台。平台2通过通信传输将数据传送到平台1进行融合,形成N*N的测量数据集。sj1,j2={(α1,j1,α2,j2,T1,j1,T2,j2)}(其中j=j1,…,jm≤N;为一个交叉点的排列,即第1个观测平台测到第j1个辐射源数据,第2个测量平台测到第j2个辐射源数据)。
2.测量数据处理
在图1中,根据简单的三角关系有:
设辐射源E到达S1、S2的时间差为ΔT,S1测量的角度为α1,则有:
联立(3)、(4)解得:
这里c为光速,
ES代表ES线段的距离。
定义时差Δtj1,j2=T1,j1-T2,j2,重新定义任一个交叉点的测量数据集sj1,j2={(α1,j1,α2,j2,ΔTj1,j2)},由(1)、(2)、(5)可以得出
如果测量组合sj1,j2来源于同一个辐射源,且没有测量噪声,则应满足如下表达式:
其中x1,y1为平台1的坐标。如果存在测量噪声,则(6)式左右不相等,定义dj1,j2为(6)式左边同(6)式右边距离的不等程度,
在dj1,j2基础上定义平台间测量数据的相近度Pj1,j2,
Pj1,j2=exp(-dj1,j2) (8)
Pj1,j2数值越大,表明平台1、2第j1组和第j2组数据来自同一辐射源概率越大,即此交叉点是真实辐射源概率越大。反之,Pj1,j2越接近0,表明平台1,2第j1组和第j2组数据来自同一辐射源概率越小,即此交叉点是虚假点概率越大。
对于N2个交叉点,分别计算其相近度P,组成一个N×N相近度矩阵
3.虚假点去除
从相近度矩阵中选出N个不同行不同列的元素,则对应着N个真实目标,其余为虚假点。从而实现虚假点去除。具体实现方法:首先从A中选出最大元素,则此元素为一真实辐射源,然后把此元素对应的行列去除,在剩下的元素中再次选择最大元素,然后再去除此元素对应的行列,依此类推,得到N个不同行不同列的元素。
以下给出两种典型辐射源编队计算实例。
设测向交叉定位系统由两平台组成,平台1在直角坐标系站址为(0km,0km),平台2为(5km,0km)。空间辐射源为5个,采用两种辐射源编队进行Monte Carlo仿真。红外测角精度为3毫弧,GPS时钟测时精度为20纳秒。定义正确定位概率为:
正确定位概率=正确定位次数/仿真次数 (9)
(1)水平编队
假设k时刻辐射源位置如下:辐射源1(100km,100km);辐射源2(100+d km,100km);辐射源3(100-d km,100km);辐射源4(100+2d km,100km);辐射源5(100-2d km,100km)。d为目标间距。图3为辐射源水平编队态势图。表1为一随机情况每个点的相近度计算。从表中可以看出,(1,1)、(2,2)、(3,3)、(4,4)、(5,5)组合最有可能是辐射源,其余都是虚假点。针对不同编队间距做1000次Monte Carlo仿真,正确定位概率如表2所示。
表1水平编队交叉点相近度(d=1km)
表2水平编队正确定位概率
间距(Km) | 0.5 | 1 | 1.5 | 2 | 2.5 | 3 |
正确定位概率 | 0.992 | 0.989 | 0.995 | 0.993 | 0.988 | 0.992 |
(2)十字编队
假设k时刻目标位置如下:辐射源1(100km,100km);辐射源2(100+d km,100km);辐射源3(100-d km,100km);辐射源4(100km,100-d km);辐射源5(100km,100+d km)。d为目标间距。图4为辐射源十字编队态势图。表3为一随机情况每个点的相近度计算。从表中可以看出,(1,1)、(2,2)、(3,3)、(4,4)、(5,5)最有可能是辐射源,其余都是虚假点。针对不同编队间距做1000次Monte Carlo仿真,正确定位概率如表4所示。
表3十字编队交叉点相近度(d=1km)
表4十字编队正确定位概率
间距(km) | 0.5 | 1 | 1.5 | 2 | 2.5 | 3 |
正确定位概率 | 0.343 | 0.546 | 0.749 | 0.836 | 0.885 | 0.891 |
通过Monte Carlo仿真可以看出,水平编队正确定位概率对辐射源间距离信息不敏感,在相同噪声情况下,正确定位概率基本上没有什么变化。十字编队正确定位概率对辐射源距离信息很敏感,随着距离的增大,正确定位概率显著增大。在此仿真条件下,对于水平编队飞行来说,角度标准差在3mrad内,时间标准差在20ns内,正确定位概率就可以达到99%,而对于十字编队来说,在相同测量精度条件下,间距d大于1.5km情况下正确定位率可以达到80%以上,都可以满足实际系统的需求。
Claims (1)
1、一种双平台多辐射源测向测时交叉定位方法,其特征在于包括如下具体步骤:
(1)测量数据采集:主平台从平台是两个独立的系统,平台携带的红外搜索与跟踪系统IRST用以测量目标的方位角,平台采用基于全球定位系统GPS的统一时间基准,一个扫描周期结束,从平台通过通信传输将数据传送到主平台进行后续处理;
(2)测量数据处理:融合主从平台的测量数据,对于N个辐射源,形成N*N组测量数据,即N*N个交叉点,根据每组测量数据分别计算出两个辐射源到主平台的斜距,根据斜距的误差,定义主从平台的测量数据的相近度P(P∈[0,1]),P的数值越大,表明平台间测量数据来自同一辐射源的可能性越大,即此组数据代表真实辐射源,重复上述过程,得到一个N*N的相近度矩阵
j1表示主平台的第j1个测量数据,j2表示从平台的第j2个测量数据;
(3)虚假点去除:从相近度矩阵A中选出最大的N个不同行不同列元素,则此N个元素对应的测量数据就代表真实辐射源,其余都为虚假目标,从而完成虚假点去除。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB031292720A CN1306280C (zh) | 2003-06-12 | 2003-06-12 | 双平台多辐射源测向测时交叉定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB031292720A CN1306280C (zh) | 2003-06-12 | 2003-06-12 | 双平台多辐射源测向测时交叉定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1477406A CN1477406A (zh) | 2004-02-25 |
CN1306280C true CN1306280C (zh) | 2007-03-21 |
Family
ID=34153465
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB031292720A Expired - Fee Related CN1306280C (zh) | 2003-06-12 | 2003-06-12 | 双平台多辐射源测向测时交叉定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1306280C (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101561493B (zh) * | 2009-05-26 | 2012-05-23 | 上海大学 | 基于双飞行机器人对有源目标的定位方法与系统 |
CN101672913B (zh) * | 2009-10-27 | 2012-09-19 | 湖南农业大学 | 激光三点动态定位方法及系统 |
CN101943750A (zh) * | 2010-08-18 | 2011-01-12 | 中华电信股份有限公司 | 3d红外线讯号源追踪装置 |
CN102175991B (zh) * | 2011-01-16 | 2012-10-31 | 哈尔滨工程大学 | 基于最大定位似然传感器配置的目标定位方法 |
CN102589549B (zh) * | 2011-12-22 | 2017-09-29 | 中国人民解放军海军航空工程学院 | 地球坐标系空间中三站测向交叉定位跟踪算法 |
CN103926569B (zh) * | 2014-01-09 | 2018-10-02 | 中国人民解放军海军航空大学 | 三坐标雷达网基于交叉定位点关联的集中式与分布式压制干扰鉴别方法 |
CN104239678B (zh) * | 2014-07-01 | 2018-03-02 | 中国电子科技集团公司第三十六研究所 | 一种实现干涉仪测向定位的方法和装置 |
CN106680772B (zh) * | 2016-12-08 | 2019-06-25 | 成都理想境界科技有限公司 | 一种信号识别系统、定位系统、定位基站及信号识别方法 |
CN108983192B (zh) * | 2018-07-27 | 2022-04-19 | 西安电子科技大学 | 基于gps辐射源的雷达运动目标参数估计方法 |
CN110673090B (zh) * | 2019-10-14 | 2022-08-05 | 电子科技大学 | 基于dbscan的无源多站多目标定位方法 |
CN111398896B (zh) * | 2020-04-02 | 2020-12-22 | 哈尔滨工程大学 | 基于基准线成对交叉定位点距离的虚假点去除预处理方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1030531A1 (en) * | 1999-02-17 | 2000-08-23 | Lucent Technologies Inc. | Method for combining multiple measurements to determine the position of a mobile transceiver |
CN1383282A (zh) * | 2001-04-20 | 2002-12-04 | Lg电子株式会社 | 判断移动通信终端位置的设备和方法 |
-
2003
- 2003-06-12 CN CNB031292720A patent/CN1306280C/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1030531A1 (en) * | 1999-02-17 | 2000-08-23 | Lucent Technologies Inc. | Method for combining multiple measurements to determine the position of a mobile transceiver |
CN1383282A (zh) * | 2001-04-20 | 2002-12-04 | Lg电子株式会社 | 判断移动通信终端位置的设备和方法 |
Also Published As
Publication number | Publication date |
---|---|
CN1477406A (zh) | 2004-02-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1306280C (zh) | 双平台多辐射源测向测时交叉定位方法 | |
Magruder et al. | Performance analysis of airborne photon-counting lidar data in preparation for the ICESat-2 mission | |
Dommety et al. | Potential networking applications of global positioning systems (GPS) | |
CN110873570B (zh) | 用于位置信息的定源、生成并且更新表示位置的地图的方法和装置 | |
CN104215239B (zh) | 基于视觉的无人机自主着陆导引装置实现的导引方法 | |
CN102411142B (zh) | 三维目标跟踪方法和系统、处理器、软件程序产品 | |
CN103196443B (zh) | 基于光流和附加信息的飞行体姿态测量方法与系统 | |
CN110988851B (zh) | 一种基于星位优化的异轨单星分时测频定位方法 | |
WO1995017685A1 (en) | Vehicle location and collision avoidance system | |
Chow et al. | Point-based versus plane-based self-calibration of static terrestrial laser scanners | |
CN102004244B (zh) | 多普勒直接测距法 | |
CN102749637A (zh) | 一种车载gps精确定位的实现方法 | |
CN106597416A (zh) | 一种地面GPS辅助的LiDAR数据高程差的误差修正方法 | |
CN102692621A (zh) | 一种ads-b与雷达的联合系统误差估计方法 | |
US20160041254A1 (en) | Method for the passive localization of radar transmitters | |
Tang et al. | The in-orbit calibration method based on terrain matching with pyramid-search for the spaceborne laser altimeter | |
CN103759725B (zh) | 一种基于六通道光电传感器的偏振方位角确定方法 | |
CN103065036A (zh) | 将传感器数据与先验数据结合的系统和方法 | |
CN113534184B (zh) | 一种激光感知的农业机器人空间定位方法 | |
CN109946648B (zh) | 一种车路协同下基于超宽带的车辆高精度定位方法 | |
CN109725339A (zh) | 一种紧耦合的自动驾驶感知方法及系统 | |
CN112462401B (zh) | 基于浮动车轨迹数据的城市峡谷快速探测方法及装置 | |
CN110456391A (zh) | 一种基于北斗geo卫星约束的单站多系统多星硬件延迟计算方法 | |
Sirikonda et al. | Integration of low-cost IMU with MEMS and NAVIC/IRNSS receiver for land vehicle navigation | |
JP7396472B2 (ja) | 位置計測装置、測位方法、及びプログラム |
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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20070321 Termination date: 20100612 |