CN116840872A - 电离层扰动检测方法及装置 - Google Patents
电离层扰动检测方法及装置 Download PDFInfo
- Publication number
- CN116840872A CN116840872A CN202311098917.5A CN202311098917A CN116840872A CN 116840872 A CN116840872 A CN 116840872A CN 202311098917 A CN202311098917 A CN 202311098917A CN 116840872 A CN116840872 A CN 116840872A
- Authority
- CN
- China
- Prior art keywords
- ionospheric
- covariance matrix
- ionosphere
- observation
- threshold
- 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
Links
- 239000005433 ionosphere Substances 0.000 title claims abstract description 44
- 238000001514 detection method Methods 0.000 title claims abstract description 28
- 239000011159 matrix material Substances 0.000 claims abstract description 55
- 238000000034 method Methods 0.000 claims abstract description 33
- 238000009826 distribution Methods 0.000 claims description 18
- 230000001186 cumulative effect Effects 0.000 claims description 8
- 238000010276 construction Methods 0.000 claims description 4
- 238000007781 pre-processing Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000001934 delay Effects 0.000 claims 1
- 238000004590 computer program Methods 0.000 description 12
- 238000010586 diagram Methods 0.000 description 9
- 238000003860 storage Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 6
- 238000005259 measurement Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000003416 augmentation Effects 0.000 description 3
- 238000012937 correction Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000013179 statistical model Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
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/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
Landscapes
- Engineering & Computer Science (AREA)
- Signal Processing (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)
Abstract
本发明涉及一种电离层扰动检测方法及装置,所述方法包括根据GNSS卫星的观测数据获取电离层延迟;基于电离层延迟构建协方差矩阵;计算协方差矩阵的最大特征值,并基于预设误警率设定阈值,对比最大特征值和阈值,根据对比结果判断是否存在电离层扰动。本发明提供的方法对中小尺度电离层扰动敏感,当扰动发生时,实时提供报警信息,具有低漏检率;同时,当判断结果为不存在电离层扰动时,则不报警,具有低误警率的优点。
Description
技术领域
本发明属于电离层物理研究技术领域,具体涉及一种电离层扰动检测方法及装置。
背景技术
全球卫星导航系统(global navigation satellite system,GNSS)信号自卫星向地面传输时需穿越电离层,由于电离层折射效应导致信号传播速度变慢,因而引发定位中的电离层误差。星基增强系统通过向用户播发电离层误差的改正数和该改正数的误差区间,以减轻电离层误差对用户的定位精度和完好性的影响。但是星基增强系统所播发的改正数依赖于5°×5°经纬度的大尺度格网模型,无法涵盖几百公里以内的中小尺度电离层扰动引发的影响,导致当局域电离层扰动发生时,星基增强系统所播发的信息可能不再准确,此时应及时发出警报,宣告此时GNSS定位服务处于不可用或需谨慎使用的状态。
相关技术中,总电子含量变化指数法(rate of total electronic contentchange,ROTI)法被定义为5分钟内每颗卫星所反演的TEC时间序列的标准差。当某个时段的ROTI明显高于其他时段时,视为电离层扰动存在。在不同的地磁暴条件下,ROTI和电离层闪烁活动以及电离层时间变化之间存在着一致性。但是由于缺乏概率分布的先验统计模型,报警阈值难以设定。因此,ROTI法仅能定性分析电离层扰动现象,而不适用于星基增强系统中的检测。
一种快速傅里叶变换方法(fast fourier transform,FFT)被作为时间异常检测器。其将TEC时间序列数据表示成三角函数或者它们的积分的线性组合,进而分离出扰动存在的时间和周期。但仅可分析单独一个卫星的时间序列数据,在多个导航卫星信号路径上均存在电离层扰动时,基于傅里叶变换的方法失效。
一种最大最小特征值检测法(maximum and minimum eigenvalue detection,MME)是计算每颗卫星反演的总电子含量(total electronic content,TEC)时序的最大特征值与最小特征值,将最大特征值与最小特征值的比值作为度量因子,根据经验设置一个阈值。若该度量因子超出阈值,则视为此时存在电离层扰动。但其检测结果异常值常常过大,这是由于协方差矩阵的最小特征值过小,因而极易诱发超过阈值,这导致误警率高。
综上所述,现有的局域电离层扰动的检测方法存在漏检率高及误警率高的问题。
发明内容
有鉴于此,本发明的目的在于克服现有技术的不足,提供一种电离层扰动检测方法及装置,以解决现有技术中局域电离层扰动的检测方法存在漏检率高及误警率高的问题。
为实现以上目的,本发明采用如下技术方案:一种电离层扰动检测方法,包括:
根据GNSS卫星的观测数据获取电离层延迟;
基于所述电离层延迟构建协方差矩阵;
计算所述协方差矩阵的最大特征值,并基于预设误警率设定阈值,对比所述最大特征值和阈值,根据对比结果判断是否存在电离层扰动。
进一步的,所述观测数据包括双频载波相位观测值和伪距观测值,所述根据GNSS卫星的观测数据获取电离层延迟,包括:
对所述双频载波相位观测值和伪距观测值进行组合,得到组合观测量;
对所述组合观测量进行预处理后,计算得到每个GNSS卫星反演的TEC时间序列;
根据所述TEC时间序列获取电离层延迟。
进一步的,所述基于所述电离层延迟构建协方差矩阵,包括:
基于所述TEC时间序列确定观测周期,以及在所述观测周期内参与计算的卫星数目;
根据所述观测周期和卫星数据,得到观测矩阵;
计算所述观测矩阵的协方差矩阵。
进一步的,采用以下方式计算所述观测矩阵的协方差矩阵,
其中,x(n,t)为观测矩阵,x +(n,t)是x(n,t)的转置,Rx为协方差矩阵,t代表当前历元;n表示当前历元地面站观测到的卫星的总数目。
进一步的,所述对比所述最大特征值和阈值,根据对比结果判断是否存在电离层扰动,包括:
若所述最大特征值大于所述阈值,则存在电离层扰动;
若所述最大特征值小于或等于所述阈值,则不存在电离层扰动。
进一步的,所述基于预设误警率设定阈值,包括:
基于预设的累积概率分布模型根据预设误警率要求确定相应的阈值;其中,所述累积概率分布模型用于描述协方差矩阵最大特征值的概率密度分布。
进一步的,还包括:
当存在电离层扰动时,根据所述最大特征值确定最大阈值,并确定所述最大阈值对应的实际误警率。
进一步的,所述根据GNSS卫星的观测数据获取电离层延迟,包括:
获取GNSS卫星输出的TEC时间序列,根据所述TEC时间序列获取电离层延迟。
本申请实施例提供一种电离层扰动检测装置,包括:
获取模块,用于根据GNSS卫星的观测数据获取电离层延迟;
构建模块,用于基于所述电离层延迟构建协方差矩阵;
对比模块,用于计算所述协方差矩阵的最大特征值,并基于预设误警率设定阈值,对比所述最大特征值和阈值,根据对比结果判断是否存在电离层扰动。
进一步的,还包括:
反推模块,用于当存在电离层扰动时,根据所述最大特征值确定最大阈值,并确定所述最大阈值对应的实际误警率。
本发明采用以上技术方案,能够达到的有益效果包括:
本发明提供一种电离层扰动检测方法及装置,本申请通过GNSS卫星的观测数据获取电离层延迟,从而构建协方差矩阵,计算协方差矩阵的最大特征值,并基于预设误警率设定阈值,对比最大特征值和阈值,根据对比结果判断是否存在电离层扰动。通过本申请提供的方法能够对中小尺度电离层扰动敏感,当扰动发生时,可实时提供报警信息,具有低漏检率;同时具有低误警率,当电离层扰动不存在时,检测方法不应报警,提高GNSS的可用性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明电离层扰动检测方法的步骤示意图;
图2为本发明电离层扰动检测方法的流程示意图;
图3为本发明电离层扰动检测装置的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将对本发明的技术方案进行详细的描述。显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所得到的所有其它实施方式,都属于本发明所保护的范围。
下面结合附图介绍本申请实施例中提供的一个具体的电离层扰动检测方法及装置。
如图1所示,本申请实施例中提供的电离层扰动检测方法,包括:
S101,根据GNSS卫星的观测数据获取电离层延迟;
一些实施例中,本申请中的观测数据包括双频载波相位观测值和伪距观测值,所述根据GNSS卫星的观测数据获取电离层延迟,包括:
对所述双频载波相位观测值和伪距观测值进行组合,得到组合观测量;
对所述组合观测量进行预处理后,计算得到每个GNSS卫星反演的TEC时间序列;
根据所述TEC时间序列获取电离层延迟。
其中,对GNSS卫星的双频载波相位观测值和伪距观测值进行组合,得到组合观测量,对组合观测量进行预处理为,将观测量中包含的差分码偏差、模糊度和周跳等误差项去除,计算得到每颗卫星所反演的TEC时间序列数据。
如果待预测的GNSS卫星能够直接输出TEC时间序列,则可直接获取电离层延迟。
S102,基于所述电离层延迟构建协方差矩阵;
一些实施例中,所述基于所述电离层延迟构建协方差矩阵,包括:
基于所述TEC时间序列确定观测周期,以及在所述观测周期内参与计算的卫星数目;
根据所述观测周期和卫星数据,得到观测矩阵;
计算所述观测矩阵的协方差矩阵。
具体的,本申请中首先通过TEC时间序列确定观测周期τ和在这整个时段该GNSS站所能收到的卫星数目n。这样就可以将观测周期、卫星数目,列成一个τ行n列或者n行τ列的观测矩阵,观测矩阵如下:
其中,I v表示TEC或电离层延迟;t代表当前历元;τ是观测周期;s n代表在当前历元,地面站可以观测到的卫星的编号;n表示当前历元地面站观测到的卫星的总数目。
然后,采用以下方式计算观测矩阵x(n,t)的协方差矩阵Rx。
其中,x +(n,t)是x(n,t)的转置。
S103,计算所述协方差矩阵的最大特征值,并基于预设误警率设定阈值,对比所述最大特征值和阈值,根据对比结果判断是否存在电离层扰动。
一些实施例中,所述对比所述最大特征值和阈值,根据对比结果判断是否存在电离层扰动,包括:
若所述最大特征值大于所述阈值,则存在电离层扰动;
若所述最大特征值小于或等于所述阈值,则不存在电离层扰动。
作为一个具体的实施例,所述基于预设误警率设定阈值,包括:
基于预设的累积概率分布模型根据预设误警率要求确定相应的阈值;其中,所述累积概率分布模型用于描述协方差矩阵最大特征值的概率密度分布。
首先计算得到协方差矩阵的最大特征值λmax,本申请中根据Tracy-Widom分布和预设的误警率要求计算阈值。其中,Tracy-Widom分布是描述协方差矩阵最大特征值的概率密度分布的一种函数,其累积概率分布γ也就是阈值可见表1;误警率P fa可根据具体要求设定,例如,设置为10-3。当误警率P fa确定后,可根据该表查到累积概率分布值γ。
表1最大特征值的概率密度分布
最后,判断最大特征值λ max和累积概率分布值γ(阈值)的大小关系。若最大特征值λ max大于阈值γ,则视为此时存在电离层扰动。
一些实施例中,还包括:
当存在电离层扰动时,根据所述最大特征值确定最大阈值,并确定所述最大阈值对应的实际误警率。
具体的,当确定存在电离层扰动时,例如最大特征值是1,误警率P fa是10-1对应的阈值是0.45,但是比最大特征值大于0.98,阈值0.98对应的5*10-2,则确定实际误警率是5*10-2。
进而,可推导该方法的实际误警率P fa为:
电离层扰动检测方法的工作原理为:本申请采用的数据是地基全球卫星导航系统站所反演的总电子含量,将该时段所有卫星的TEC时间序列数据列入矩阵中,计算该矩阵的协方差矩阵的最大特征值,若最大特征值超出理论阈值,则视为此时发生电离层扰动。其中,阈值通过统计学理论严格推导得出。
具体的,由于协方差最大特征值可以表征该矩阵的最大特征量,即电离层的变化量。以最大特征值作为度量因子表征电离层扰动,避免了除以最小特征值的环节,最小特征值往往过小,它若作为分母往往导致度量因子过大,使得报警频繁,误警率较高,由此导致GNSS的可用性降低。本申请采用仅基于最大特征值的度量因子拥有已被多次验证的Tracy-Widom概率分布模型作为理论参考,使得阈值设立更加合理。
如图2所示,本申请实施例提供一种电离层扰动检测装置,包括:
获取模块201,用于根据GNSS卫星的观测数据获取电离层延迟;
构建模块202,用于基于所述电离层延迟构建协方差矩阵;
对比模块203,用于计算所述协方差矩阵的最大特征值,并基于预设误警率设定阈值,对比所述最大特征值和阈值,根据对比结果判断是否存在电离层扰动。
本申请提供的电离层扰动检测装置的工作原理为,获取模块201根据GNSS卫星的观测数据获取电离层延迟;构建模块202基于所述电离层延迟构建协方差矩阵;对比模块203计算所述协方差矩阵的最大特征值,并基于预设误警率设定阈值,对比所述最大特征值和阈值,根据对比结果判断是否存在电离层扰动。
一些实施例中,本申请提供的电离层扰动检测装置,还包括:
反推模块,用于当存在电离层扰动时,根据所述最大特征值确定最大阈值,并确定所述最大阈值对应的实际误警率。
本申请提供一种计算机设备,包括:存储器1和处理器2,还可以包括网络接口3,所述存储器存储有计算机程序,存储器可以包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM)。该计算机设备存储有操作系统4,存储器是计算机可读介质的示例。所述计算机程序被所述处理器执行时,使得所述处理器执行电离层扰动检测方法,图3中示出的结构,仅仅是与本申请方案相关的部分结构的框图,并不构成对本申请方案所应用于其上的计算机设备的限定,具体的计算机设备可以包括比图中所示更多或更少的部件,或者组合某些部件,或者具有不同的部件布置。
在一个实施例中,本申请提供的电离层扰动检测方法可以实现为一种计算机程序的形式,计算机程序可在如图3所示的计算机设备上运行。
一些实施例中,所述计算机程序被所述处理器执行时,使得所述处理器执行以下步骤:根据GNSS卫星的观测数据获取电离层延迟;基于所述电离层延迟构建协方差矩阵;计算所述协方差矩阵的最大特征值,并基于预设误警率设定阈值,对比所述最大特征值和阈值,根据对比结果判断是否存在电离层扰动。
本申请还提供一种计算机存储介质,计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光光盘(DVD)或其他光学存储、磁盒式磁带存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。
一些实施例中,本发明还提出了一种计算机可读存储介质,存储有计算机程序,所述计算机程序被处理器执行时,根据GNSS卫星的观测数据获取电离层延迟;基于所述电离层延迟构建协方差矩阵;计算所述协方差矩阵的最大特征值,并基于预设误警率设定阈值,对比所述最大特征值和阈值,根据对比结果判断是否存在电离层扰动。
综上所述,本发明提供一种电离层扰动检测方法及装置,方法包括根据GNSS卫星的观测数据获取电离层延迟;基于电离层延迟构建协方差矩阵;计算协方差矩阵的最大特征值,并基于预设误警率设定阈值,对比最大特征值和阈值,根据对比结果判断是否存在电离层扰动。本发明提供的方法对中小尺度电离层扰动敏感,当扰动发生时,实时提供报警信息,具有低漏检率;同时,当判断结果为不存在电离层扰动时,则不报警,具有低误警率的优点。
可以理解的是,上述提供的方法实施例与上述的装置实施例对应,相应的具体内容可以相互参考,在此不再赘述。
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器和光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令方法的制造品,该指令方法实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。
Claims (10)
1.一种电离层扰动检测方法,其特征在于,包括:
根据GNSS卫星的观测数据获取电离层延迟;
基于所述电离层延迟构建协方差矩阵;
计算所述协方差矩阵的最大特征值,并基于预设误警率设定阈值,对比所述最大特征值和阈值,根据对比结果判断是否存在电离层扰动。
2.根据权利要求1所述的方法,其特征在于,所述观测数据包括双频载波相位观测值和伪距观测值,所述根据GNSS卫星的观测数据获取电离层延迟,包括:
对所述双频载波相位观测值和伪距观测值进行组合,得到组合观测量;
对所述组合观测量进行预处理后,计算得到每个GNSS卫星反演的TEC时间序列;
根据所述TEC时间序列获取电离层延迟。
3.根据权利要求2所述的方法,其特征在于,所述基于所述电离层延迟构建协方差矩阵,包括:
基于所述TEC时间序列确定观测周期,以及在所述观测周期内参与计算的卫星数目;
根据所述观测周期和卫星数据,得到观测矩阵;
计算所述观测矩阵的协方差矩阵。
4.根据权利要求3所述的方法,其特征在于,采用以下方式计算所述观测矩阵的协方差矩阵,
其中,x(n, t)为观测矩阵,x +(n, t)是x(n, t)的转置,Rx为协方差矩阵,t代表当前历元;n表示当前历元地面站观测到的卫星的总数目。
5.根据权利要求1所述的方法,其特征在于,所述对比所述最大特征值和阈值,根据对比结果判断是否存在电离层扰动,包括:
若所述最大特征值大于所述阈值,则存在电离层扰动;
若所述最大特征值小于或等于所述阈值,则不存在电离层扰动。
6.根据权利要求1所述的方法,其特征在于,所述基于预设误警率设定阈值,包括:
基于预设的累积概率分布模型根据预设误警率要求确定相应的阈值;其中,所述累积概率分布模型用于描述协方差矩阵最大特征值的概率密度分布。
7.根据权利要求1所述的方法,其特征在于,还包括:
当存在电离层扰动时,根据所述最大特征值确定最大阈值,并确定所述最大阈值对应的实际误警率。
8.根据权利要求1所述的方法,其特征在于,所述根据GNSS卫星的观测数据获取电离层延迟,包括:
获取GNSS卫星输出的TEC时间序列,根据所述TEC时间序列获取电离层延迟。
9.一种电离层扰动检测装置,其特征在于,包括:
获取模块,用于根据GNSS卫星的观测数据获取电离层延迟;
构建模块,用于基于所述电离层延迟构建协方差矩阵;
对比模块,用于计算所述协方差矩阵的最大特征值,并基于预设误警率设定阈值,对比所述最大特征值和阈值,根据对比结果判断是否存在电离层扰动。
10.根据权利要求9所述的装置,其特征在于,还包括:
反推模块,用于当存在电离层扰动时,根据所述最大特征值确定最大阈值,并确定所述最大阈值对应的实际误警率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311098917.5A CN116840872A (zh) | 2023-08-30 | 2023-08-30 | 电离层扰动检测方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311098917.5A CN116840872A (zh) | 2023-08-30 | 2023-08-30 | 电离层扰动检测方法及装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116840872A true CN116840872A (zh) | 2023-10-03 |
Family
ID=88167427
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311098917.5A Pending CN116840872A (zh) | 2023-08-30 | 2023-08-30 | 电离层扰动检测方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116840872A (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007032947A1 (en) * | 2005-09-09 | 2007-03-22 | Trimble Navigation Limited | Ionosphere modeling apparatus and methods |
WO2013003636A2 (en) * | 2011-06-28 | 2013-01-03 | Nextnav, Llc | Coding in a wide area positioning system (waps) |
CN106093967A (zh) * | 2016-08-22 | 2016-11-09 | 中国科学院上海天文台 | 一种伪距相位综合的电离层延迟求解方法 |
CN108363084A (zh) * | 2018-01-18 | 2018-08-03 | 和芯星通科技(北京)有限公司 | 利用卫星定位的方法和装置、卫星导航接收机、存储介质 |
CN110687559A (zh) * | 2019-11-04 | 2020-01-14 | 中国电子科技集团公司第五十四研究所 | 适用于机载的gnss无缝高精度定位和完好性评估方法 |
CN111175796A (zh) * | 2020-01-20 | 2020-05-19 | 桂林电子科技大学 | 一种网络rtk中长基线模糊度快速解算方法 |
US20210173090A1 (en) * | 2019-12-06 | 2021-06-10 | Raytheon Company | Ionospheric gradient h1 monitor |
CN114690207A (zh) * | 2020-12-31 | 2022-07-01 | 复旦大学 | 差分电离层建模方法及系统 |
-
2023
- 2023-08-30 CN CN202311098917.5A patent/CN116840872A/zh active Pending
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007032947A1 (en) * | 2005-09-09 | 2007-03-22 | Trimble Navigation Limited | Ionosphere modeling apparatus and methods |
WO2013003636A2 (en) * | 2011-06-28 | 2013-01-03 | Nextnav, Llc | Coding in a wide area positioning system (waps) |
CN106093967A (zh) * | 2016-08-22 | 2016-11-09 | 中国科学院上海天文台 | 一种伪距相位综合的电离层延迟求解方法 |
CN108363084A (zh) * | 2018-01-18 | 2018-08-03 | 和芯星通科技(北京)有限公司 | 利用卫星定位的方法和装置、卫星导航接收机、存储介质 |
CN110687559A (zh) * | 2019-11-04 | 2020-01-14 | 中国电子科技集团公司第五十四研究所 | 适用于机载的gnss无缝高精度定位和完好性评估方法 |
US20210173090A1 (en) * | 2019-12-06 | 2021-06-10 | Raytheon Company | Ionospheric gradient h1 monitor |
CN111175796A (zh) * | 2020-01-20 | 2020-05-19 | 桂林电子科技大学 | 一种网络rtk中长基线模糊度快速解算方法 |
CN114690207A (zh) * | 2020-12-31 | 2022-07-01 | 复旦大学 | 差分电离层建模方法及系统 |
Non-Patent Citations (5)
Title |
---|
CHUNHUA JIANG等: "A parallel equivalence algorithm based on MPI for GNSS data processing", JOURNAL OF SPATIAL SCIENCE, pages 1 - 21 * |
JACEK PAZIEWSKI等: "Study on desirable ionospheric corrections accuracy for network-RTK positioning and its impact on time-to-fix and probability of successful single-epoch ambiguity resolution", SCIENCEDIRECT, pages 1098 - 1111 * |
TONG LIU等: "Maximum eigenvalue detection for local ionospheric anomaly - theory, cases, statistics, and potential application", MEASUREMENT, pages 1 - 9 * |
徐炜: "GNSS中长基线解算技术与程序实现", 中国优秀硕士论文全文数据库 基础科学辑, pages 008 - 165 * |
董毅: "GNSS/INS组合动对动精密相对定位方法研究", 中国博士学位论文全文数据库 信息科技辑, pages 136 - 232 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107064961B (zh) | 对卫星导航系统完好性监测性能进行测试的方法及装置 | |
Castaldo et al. | P‐RANSAC: An Integrity Monitoring Approach for GNSS Signal Degraded Scenario | |
US10416315B2 (en) | False alarm distribution in advanced receiver autonomous integrity monitoring | |
Zhao et al. | High-rate Doppler-aided cycle slip detection and repair method for low-cost single-frequency receivers | |
CN109085619B (zh) | 多模gnss系统的定位方法及装置、存储介质、接收机 | |
CN112068161B (zh) | 一种多路径误差削减方法及装置 | |
Pesonen | A framework for Bayesian receiver autonomous integrity monitoring in urban navigation | |
Li et al. | Integrity and continuity allocation for the RAIM with multiple constellations | |
CN111856513A (zh) | 卫星观测值获取方法、装置、计算机设备和存储介质 | |
CN115390095A (zh) | 一种获取电离层延迟的方法、装置及介质 | |
CN112731460A (zh) | 卫星导航系统校准数据的完好性风险的测试方法及其系统 | |
CN111142125B (zh) | 一种卫星完好性监测方法和系统 | |
CN112198533B (zh) | 一种多假设下的地基增强系统完好性评估系统及方法 | |
Seepersad et al. | Integrity monitoring in precise point positioning | |
Gao et al. | An improved detection method of GNSS faults with fractional information divergence | |
CN116840874A (zh) | 差分数据质量监测方法、装置、电子设备及存储介质 | |
CN115616622B (zh) | 一种故障检测方法、装置、设备及介质 | |
CN116840872A (zh) | 电离层扰动检测方法及装置 | |
Osechas et al. | Carrier‐Phase Acceleration RAIM for GNSS Satellite Clock Fault Detection | |
CN112733315B (zh) | 完好性风险的验证方法及系统 | |
Tu et al. | Evaluation of GNSS Receiver Autonomous Integrity Monitoring for multiple outliers with a smart random sample consensus strategy | |
CN112444825A (zh) | 一种基于北斗geo的电离层tec同化建模方法 | |
CN107783166B (zh) | Gps速度异常检测与修复的方法及其系统 | |
Ni et al. | An improved RAIM scheme for processing multiple outliers in GNSS | |
Gioia et al. | Inter-system Bias: Stability and impact on multi-constellation positioning |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20231003 |