CN110988922B - 一种基于ads-b的三维椭球交叉gnss干扰源定位方法 - Google Patents
一种基于ads-b的三维椭球交叉gnss干扰源定位方法 Download PDFInfo
- Publication number
- CN110988922B CN110988922B CN201910996596.8A CN201910996596A CN110988922B CN 110988922 B CN110988922 B CN 110988922B CN 201910996596 A CN201910996596 A CN 201910996596A CN 110988922 B CN110988922 B CN 110988922B
- Authority
- CN
- China
- Prior art keywords
- ads
- equation
- point
- interference source
- gnss
- 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
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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/21—Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于ADS‑B的三维椭球交叉GNSS干扰源定位方法,包括如下步骤:步骤1:获取ADS‑B接收机的所有飞机的ADS‑B数据信息;步骤2:搜索ADS‑B数据信息中有位置信息丢失的航线;步骤3:分析每条位置信息丢失航线的位置丢失点和位置重获取点;步骤4,利用线性化降维及最小二乘法实现三维椭球GNSS干扰源定位。本发明所公开基于ADS‑B的三维椭球交叉GNSS干扰源定位方法,利用机场逐步推广的ADS‑B(广播式自动相关监视)系统中的GNSS信息,不需要增加额外设备以及传输链路,成本低、易实现,为大范围的GNSS干扰源定位提供可能。
Description
技术领域
本发明属于GNSS电磁干扰定位技术领域,特别涉及该领域中的一种基于ADS-B的三维椭球交叉GNSS干扰源定位方法。
背景技术
全球导航卫星系统(GNSS)的应用十分广泛,小至个人的使用,大到国家关键基础设施的应用。但是GNSS信号到达地面的功率很低,同时民用信号的格式是公开的,因此极易受到各种无意和人为故意干扰的影响,很小的干扰也会对GNSS定位和授时的精度造成影响。
近些年来GNSS干扰引起了国际的广泛关注,及时发现GNSS干扰源并排除掉干扰源,才能够保障GNSS在各个方面应用的性能。人们对GNSS干扰源定位做了大量的研究,目前常用的GNSS干扰源定位方法包括到达时差(TDOA)定位、测向交汇定位、到达频率差定位技术等。这些技术都需要布设专用设备才能够实现,但是很小的GNSS干扰源就能够影响到GNSS,要对很小的GNSS干扰源进行定位,需要比较密集的布设专用设备,布设费用十分昂贵。
发明内容
本发明所要解决的技术问题就是提供一种成本低、易实现的基于ADS-B的三维椭球交叉GNSS干扰源定位方法。
本发明采用如下技术方案:
一种基于ADS-B的三维椭球交叉GNSS干扰源定位方法,其改进之处在于,包括如下步骤:
步骤1:获取ADS-B接收机的所有飞机的ADS-B数据信息;
步骤2:搜索ADS-B数据信息中有位置信息丢失的航线;
步骤3:分析每条位置信息丢失航线的位置丢失点和位置重获取点;
步骤4,利用线性化降维及最小二乘法实现三维椭球GNSS干扰源定位:
对于第k条有位置丢失点和重获取点的航线,建立以下方程:
其中,xJ、yJ、zJ分别为干扰源的三维坐标,xLk、yLk、zLk分别为第k条有位置丢失点和重捕获点航线的位置丢失点的三维坐标,xRk、yRk、zRk分别为第k条有位置丢失点和重捕获点航线的位置重获取点的三维坐标;
步骤51:对方程(1)进行去根号,可得:
步骤52:对方程(2)进行线性化降维,记:
通过寻找GNSS干扰源的初始位置(xJ,0,yJ,0,zJ,0)开始线性化,更新的干扰源的坐标为:
其中,ΔxJ1、ΔyJ1、ΔzJ1为增量,f(xJ,1,yJ,1,zJ,1)的泰勒展开式为:
经过这种变化后,就实现了方程(2)的线性化方程:
步骤53:对于每一条有位置丢失点和重捕获点航线的位置丢失点的航线,都可以建立方程(6),当多于三条及三条以上这样的航线时,可建立求解干扰源位置的方程组:
对于该方程组,采用最小二乘进行求解,最小二乘是给出一个无解的方程Ax=b,A有m行和n列,m>n;其中观测值b1,...,bm的个数大于求解参数的个数,称之为最优解应使误差向量的长度最小,则等价于使得||e||2=(b-Ax)T(b-Ax)最小,最小化该二次方程得到:
采用最小二乘可求得方程组(7)的解ΔxJ1、ΔyJ1、ΔzJ1;
步骤54:求得ΔxJ1、ΔyJ1、ΔzJ1后,可由方程(4)求得干扰源的新位置,将该新位置作为新的初始位置(xJ,0,yJ,0,zJ,0),从方程(5)到方程(7)循环计算,直到增量ΔxJ1、ΔyJ1、ΔzJ1收敛到设定阈值为止。
本发明的有益效果是:
本发明所公开基于ADS-B的三维椭球交叉GNSS干扰源定位方法,利用机场逐步推广的ADS-B(广播式自动相关监视)系统中的GNSS信息,不需要增加额外设备以及传输链路,成本低、易实现,为大范围的GNSS干扰源定位提供可能。采用ADS-B接收机接收的飞机数据,基于有位置丢失和位置重获取的航线中的位置丢失点和位置重获取点,结合传播模型,能够给出较为准确的干扰源位置,为干扰源的精确查找及排除提供位置参考。
附图说明
图1是本发明实施例1所公开方法的流程示意图;
图2是ADS-B信息中提取的受GNSS干扰影响的航线示意图;
图3是某一条受GNSS干扰影响的航线所确定干扰源的椭球示意图;
图4是基于ADS-B的三维椭球交叉GNSS干扰源定位结果示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图和实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
实施例1,如图1所示,本实施例公开了一种基于ADS-B的三维椭球交叉GNSS干扰源定位方法,包括如下步骤:
步骤1:获取ADS-B接收机的所有飞机的ADS-B数据信息;
步骤2:搜索ADS-B数据信息中有位置信息丢失的航线;
步骤3:分析每条位置信息丢失航线的位置丢失点和位置重获取点;
步骤4,利用线性化降维及最小二乘法实现三维椭球GNSS干扰源定位:
对于第k条有位置丢失点和重获取点的航线,建立以下方程:
其中,xJ、yJ、zJ分别为干扰源的三维坐标,xLk、yLk、zLk分别为第k条有位置丢失点和重捕获点航线的位置丢失点的三维坐标,xRk、yRk、zRk分别为第k条有位置丢失点和重捕获点航线的位置重获取点的三维坐标;
步骤51:对方程(1)进行去根号,可得:
步骤52:对方程(2)进行线性化降维,记:
通过寻找GNSS干扰源的初始位置(xJ,0,yJ,0,zJ,0)开始线性化,更新的干扰源的坐标为:
其中,ΔxJ1、ΔyJ1、ΔzJ1为增量,f(xJ,1,yJ,1,zJ,1)的泰勒展开式为:
经过这种变化后,就实现了方程(2)的线性化方程:
步骤53:对于每一条有位置丢失点和重捕获点航线的位置丢失点的航线,都可以建立方程(6),当多于三条及三条以上这样的航线时,可建立求解干扰源位置的方程组:
对于该方程组,采用最小二乘进行求解,最小二乘是给出一个无解的方程Ax=b,A有m行和n列,m>n;其中观测值b1,...,bm的个数大于求解参数的个数,称之为最优解应使误差向量的长度最小,则等价于使得||e||2=(b-Ax)T(b-Ax)最小,最小化该二次方程得到:
采用最小二乘可求得方程组(7)的解ΔxJ1、ΔyJ1、ΔzJ1;
步骤54:求得ΔxJ1、ΔyJ1、ΔzJ1后,可由方程(4)求得干扰源的新位置,将该新位置作为新的初始位置(xJ,0,yJ,0,zJ,0),从方程(5)到方程(7)循环计算,直到增量ΔxJ1、ΔyJ1、ΔzJ1收敛到设定阈值为止。
图2是ADS-B信息中提取的受GNSS干扰影响的航线示意图;图3是某一条受GNSS干扰影响的航线所确定干扰源的椭球示意图。
如图4所示,为了验证该方法在GNSS干扰源定位方面的性能,本实施例选取了10条航线进行仿真验证,其中干扰源的坐标为(20km,20km,0km),干扰源的功率为30dBm,随机生成10条航线,加入1dBm的随机误差,生成各航线的位置丢失点和位置重获取点,经过计算影响的范围为50km,采用本实施例提供的方法对干扰源进行定位,经过100次仿真测试,定位干扰源的平均位置为(19.9999km,19.9999km,-0.1007km),定位误差为100.6890m。
综上所述,本实施例提供的基于ADS-B的三维椭球交叉GNSS干扰源定位方法,可以实现干扰源的较准确定位,为干扰源的进一步排查和确认提供参考,对大范围GNSS干扰源的检测定位具有重要价值。
Claims (1)
1.一种基于ADS-B的三维椭球交叉GNSS干扰源定位方法,其特征在于,包括如下步骤:
步骤1:获取ADS-B接收机的所有飞机的ADS-B数据信息;
步骤2:搜索ADS-B数据信息中有位置信息丢失的航线;
步骤3:分析每条位置信息丢失航线的位置丢失点和位置重获取点;
步骤4,利用线性化降维及最小二乘法实现三维椭球GNSS干扰源定位:
对于第k条有位置丢失点和重获取点的航线,建立以下方程:
其中,xJ、yJ、zJ为干扰源的三维坐标,xLk、yLk、zLk为第k条有位置丢失点和重获取点航线的位置丢失点的三维坐标,xRk、yRk、zRk为第k条有位置丢失点和重获取点航线的位置重获取点的三维坐标;
步骤41:对方程(1)进行去根号,可得:
步骤42:对方程(2)进行线性化降维,记:
通过寻找GNSS干扰源的初始位置(xJ,0,yJ,0,zJ,0)开始线性化,更新的干扰源的坐标为:
其中,ΔxJ1、ΔyJ1、ΔzJ1为增量,fk(xJ,1,yJ,1,zJ,1)的泰勒展开式为:
经过这种变化后,就实现了方程(2)的线性化方程:
步骤43:对于有位置丢失点和重获取点的航线,都建立方程(6),当多于三条及三条以上这样的航线时,可建立求解干扰源位置的方程组:
对于该方程组,采用最小二乘进行求解,最小二乘是给出一个无解的方程Ax=b,A有m行和n列,m>n;其中观测值b1,...,bm的个数大于求解参数的个数,称之为最优解应使误差向量的长度最小,则等价于使得||e||2=(b-Ax)T(b-Ax)最小,最小化该二次方程得到:
采用最小二乘可求得方程组(7)的解ΔxJ1、ΔyJ1、ΔzJ1;
步骤44:求得ΔxJ1、ΔyJ1、ΔzJ1后,由方程(4)求得干扰源的新位置,将该新位置作为新的初始位置(xJ,0,yJ,0,zJ,0),从方程(5)到方程(7)循环计算,直到增量ΔxJ1、ΔyJ1、ΔzJ1收敛到设定阈值为止。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910996596.8A CN110988922B (zh) | 2019-10-19 | 2019-10-19 | 一种基于ads-b的三维椭球交叉gnss干扰源定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910996596.8A CN110988922B (zh) | 2019-10-19 | 2019-10-19 | 一种基于ads-b的三维椭球交叉gnss干扰源定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110988922A CN110988922A (zh) | 2020-04-10 |
CN110988922B true CN110988922B (zh) | 2021-09-17 |
Family
ID=70082176
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910996596.8A Active CN110988922B (zh) | 2019-10-19 | 2019-10-19 | 一种基于ads-b的三维椭球交叉gnss干扰源定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110988922B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112558113A (zh) * | 2020-12-04 | 2021-03-26 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种利用ads-b的基于网格概率遍历的gnss干扰源定位方法 |
CN113759392B (zh) * | 2021-08-23 | 2023-08-18 | 中国民航大学 | 基于飞行大数据的鲁棒gnss干扰源定位方法 |
CN113568010B (zh) * | 2021-08-31 | 2024-03-26 | 中国电子科技集团公司第五十四研究所 | 一种基于星载ads-b的干扰定位方法 |
CN115480271A (zh) * | 2022-09-13 | 2022-12-16 | 北京邮电大学 | Gps干扰源定位方法、装置、电子设备和可读存储介质 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103024894A (zh) * | 2011-09-21 | 2013-04-03 | 华为技术有限公司 | 多基站干扰源定位方法及设备 |
EP2846174A1 (en) * | 2013-09-05 | 2015-03-11 | MBDA UK Limited | Geo-location of Jamming Signals |
CN104570022A (zh) * | 2014-12-26 | 2015-04-29 | 北京理工雷科电子信息技术有限公司 | 一种基于dsp的多系统兼容gnss接收机及其接收系统 |
CN205607406U (zh) * | 2016-05-17 | 2016-09-28 | 怀化市东成电子有限公司 | 一种具有陀螺仪和加速度计的车载定位终端 |
CN106959453A (zh) * | 2017-05-18 | 2017-07-18 | 北京理工大学 | 一种用于辅助任务无人机获取卫星信号的辅助无人机 |
CN109471061A (zh) * | 2018-09-18 | 2019-03-15 | 宁波大学 | 一种鲁棒处理模型参数误差的接收信号强度差定位方法 |
CN109541538A (zh) * | 2018-12-12 | 2019-03-29 | 华东理工大学 | 基于机会参考源实现分布式ads-b无源定位防欺骗功能的方法 |
CN109581425A (zh) * | 2018-12-29 | 2019-04-05 | 南京天际易达通信技术有限公司 | 一种基于多接收机的卫星导航欺骗信号检测方法 |
CN109633695A (zh) * | 2019-01-10 | 2019-04-16 | 浙江大学 | 一种无人机对卫导干扰机的主动定位方法 |
-
2019
- 2019-10-19 CN CN201910996596.8A patent/CN110988922B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103024894A (zh) * | 2011-09-21 | 2013-04-03 | 华为技术有限公司 | 多基站干扰源定位方法及设备 |
EP2846174A1 (en) * | 2013-09-05 | 2015-03-11 | MBDA UK Limited | Geo-location of Jamming Signals |
CN104570022A (zh) * | 2014-12-26 | 2015-04-29 | 北京理工雷科电子信息技术有限公司 | 一种基于dsp的多系统兼容gnss接收机及其接收系统 |
CN205607406U (zh) * | 2016-05-17 | 2016-09-28 | 怀化市东成电子有限公司 | 一种具有陀螺仪和加速度计的车载定位终端 |
CN106959453A (zh) * | 2017-05-18 | 2017-07-18 | 北京理工大学 | 一种用于辅助任务无人机获取卫星信号的辅助无人机 |
CN109471061A (zh) * | 2018-09-18 | 2019-03-15 | 宁波大学 | 一种鲁棒处理模型参数误差的接收信号强度差定位方法 |
CN109541538A (zh) * | 2018-12-12 | 2019-03-29 | 华东理工大学 | 基于机会参考源实现分布式ads-b无源定位防欺骗功能的方法 |
CN109581425A (zh) * | 2018-12-29 | 2019-04-05 | 南京天际易达通信技术有限公司 | 一种基于多接收机的卫星导航欺骗信号检测方法 |
CN109633695A (zh) * | 2019-01-10 | 2019-04-16 | 浙江大学 | 一种无人机对卫导干扰机的主动定位方法 |
Non-Patent Citations (5)
Title |
---|
Detection and Localization of GNSS Radio Interference using ADS-B Data;Petr Jonáš*,et al;《International Conference on Military Technologies (ICMT) 2019》;20190531;p1-5 * |
Interference Localization for Satellite Navigation Systems;Andrew G. Dempster, et al;《Proceedings of the IEEE》;20160630;p1318-1326 * |
一种网格化的GNSS干扰源定位方法研究;刘睿等;《全球定位系统》;20151231;第16-20页 * |
干扰信号强度补偿的传感器网络干扰源定位;王胜等;《传感器与微系统》;20161231;第60-65页 * |
空管设备电磁干扰源定位算法研究;夏涛;《交通世界》;20181231;第14-15页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110988922A (zh) | 2020-04-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110988922B (zh) | 一种基于ads-b的三维椭球交叉gnss干扰源定位方法 | |
CN109298388B (zh) | 基于方位信息的超视距目标地理坐标直接估计方法 | |
CN106886039B (zh) | 基于城市三维地图的地面数字电视广播与导航卫星混合粒子滤波定位方法 | |
CN109061774B (zh) | 一种雷暴核关联性处理方法 | |
CN103270801B (zh) | 用于定位用户设备的方位的方法和系统 | |
CN104507050A (zh) | 一种WiFi室内定位中概率型指纹匹配方法 | |
CN110673087B (zh) | 基于ads-b广播信号飞机航路监测方法、装置及计算机存储介质 | |
CN110988929B (zh) | 电离层影响下的gbas系统性能评估方法及装置 | |
CN107607907B (zh) | 一种基于天基随机无线电信号的定位导航系统 | |
CN113740882B (zh) | 利用qar/ads-b数据基于到达功率差的稳健gnss干扰源定位方法 | |
Tobie et al. | Hybrid navigation filters performances between GPS, Galileo and 5G TOA measurements in multipath environment | |
CN117031519A (zh) | 一种多目标三维协同定位方法及系统 | |
KR102520764B1 (ko) | Uwb 측위 시스템에서의 시각 동기화 방법 및 위치 추정 방법 | |
Xu et al. | Assimilation of high frequency radar data into a shelf sea circulation model | |
CN112558113A (zh) | 一种利用ads-b的基于网格概率遍历的gnss干扰源定位方法 | |
Vicenzo et al. | Experimental Investigation of GNSS Direct Position Estimation in Densely Urban Area | |
CN113759392B (zh) | 基于飞行大数据的鲁棒gnss干扰源定位方法 | |
Jan et al. | Evaluation of positioning algorithms for wide area multilateration based alternative positioning navigation and timing (APNT) using 1090 MHz ADS-B signals | |
CN105607038A (zh) | 基于gpu实现的运动卫星干扰源定位方法 | |
Thevenon et al. | Pseudo-range measurements using OFDM channel estimation | |
Chu et al. | Multi-receiver direct position estimation tested on a full-scale fixed-wing aircraft | |
Lee et al. | Generation of Klobuchar coefficients for ionospheric error simulation | |
Zeller et al. | Relative train localization for cooperative maneuvers using GNSS pseudoranges and geometric track information | |
Guanli et al. | New algorithms on the solution to drifting problem of GPS positioning | |
Seo et al. | Evaluation of deep signal fading effects due to ionospheric scintillation on GPS aviation receivers |
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 |