CN110907953A - 一种卫星故障识别方法、装置及软件接收机 - Google Patents
一种卫星故障识别方法、装置及软件接收机 Download PDFInfo
- Publication number
- CN110907953A CN110907953A CN201910996032.4A CN201910996032A CN110907953A CN 110907953 A CN110907953 A CN 110907953A CN 201910996032 A CN201910996032 A CN 201910996032A CN 110907953 A CN110907953 A CN 110907953A
- Authority
- CN
- China
- Prior art keywords
- satellite
- value
- matrix
- positioning
- false alarm
- 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.)
- Granted
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/20—Integrity monitoring, fault detection or fault isolation of space segment
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Security & Cryptography (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种卫星故障识别方法、装置及软件接收机,涉及卫星导航与定位技术领域;该方法包括:采集接收机位置数据和钟差,根据接收机位置数据、钟差构建观测矩阵并采用奇偶空间算法对观测矩阵进行分解,生成奇偶空间矩阵;对奇偶空间矩阵进行矩阵运算,得到参与定位的每个卫星的斜率并提取最大斜率值;基于最大斜率值和误警检测值计算定位卫星的相似径向误差保护值,当相似径向误差保护值不小于预设值时,根据奇偶空间矩阵对应的卫星偏差数据计算校验统计量;根据校验统计量、误警检测值和预置的阈值因子检测参与定位的卫星中是否存在故障;本发明具有时间复杂度低、空间复杂度低的特点,能够适应常态飞行、高动态飞行的卫星故障检测和识别。
Description
技术领域
本发明属于卫星导航与定位技术领域,更具体地,涉及一种卫星故障识别方法、装置及软件接收机。
背景技术
在软件接收机的用户位置、速度以及时间信息(Position Velocity Time,PVT)解算 中,需要实时地输出用户经度、纬度、高度和速度等测量信息。为了能够提高接收机的导航精度,必须实时地对卫星进行监测。然而复杂环境下多星座可见卫星的变化容易引 起系统故障,导致出现检测性能不稳定的问题。因此,需要设计一种接收机监测、识别 故障卫星的算法(RAIM,Receiver Autonomous Integrity Monitoring),当观测信息的误差 超过阈值后,将启动该算法使其能够成功剔除故障卫星,从而保证用户实时接收高精度、 高可靠性的定位信息。
启动RAIM算法至少需要5颗卫星参与定位,进行故障卫星剔除至少需要6颗卫星参与定位。
传统的RAIM算法主要包括最小二乘残差算法、奇偶空间校验的算法,二者在数学上是等价的,故二者具有相同的卫星伪距跳变识别度。但是二者的运算角度具有较大的 差别,最小二乘残差是在欧式距离上进行计算的,计算的空间复杂度高;而奇偶空间校 验是在奇偶空间上进行计算的,计算的时间复杂度高。由于弹载接收机内部DSP资源 有限,因此对算法的时间复杂度、空间复杂度都具有较高要求。
因此,研制一种时间复杂度低、空间复杂度低,并且能够适应常态飞行和高动态飞行,并且具有较高的卫星伪距跳变识度的算法至关重要。
发明内容
针对现有技术的至少一个缺陷或改进需求,本发明提供了一种卫星故障识别方法、 装置及软件接收机,结合最小二乘残差算法和奇偶空间校验算法的各自优点,具有时间 复杂度低、空间复杂度低的特点,并且具有相同的卫星伪距跳变识别度,使其能够适应常态飞行、高动态飞行的卫星故障检测和识别;其目的在于解决现有方法存在的时间复 杂度、空间复杂度高的问题。
为实现上述目的,按照本发明的第一个方面,提供了一种卫星故障识别方法,包括以下步骤:
S1:采集基于最小二乘算法进行卫星定位时收敛的接收机位置数据和接收机钟差, 根据所述接收机位置数据、接收机钟差构建观测矩阵并采用奇偶空间算法对所述观测矩 阵进行分解,生成奇偶空间矩阵;
S2:采用奇偶空间算法对所述标准正交矩阵进行分解,生成奇偶空间矩阵;
S3:对所述奇偶空间矩阵进行矩阵运算,得到参与定位的每个卫星的斜率并提取最 大斜率值;
S4:根据所述最大斜率值和误警检测值计算相似径向误差保护值,当所述相似径向 误差保护值不小于预设值时,根据所述奇偶空间矩阵对应的卫星偏差数据计算校验统计 量;
S5:根据所述校验统计量、误警检测值和预置的阈值因子检测参与定位的卫星中是 否存在故障。
优选的,上述卫星故障识别方法,其步骤S3中具体包括:
(1)根据所述奇偶空间矩阵构造第一矩阵和第二矩阵;
B=In-Qp TQp
A=H\B;
其中,Qp表示奇偶空间矩阵;B表示第一矩阵;In表示单位矩阵;A表示第二矩阵; H表示观测矩阵;n表示参与定位的卫星数量;
优选的,上述卫星故障识别方法,根据预置的正态分布方差、误报警率计算所述误警检测值:
其中,σγ表示误警检测值;表示正态分布方差;PFA表示误报警率;n表示参与 定位的卫星数量;icdf('ch2',1-PFA,n-4)表示求取以(n-4)为自由度,累计概率为 (1-PFA)的χ2分布的临界值。
优选的,上述卫星故障识别方法,其步骤S5中具体包括:
其中,FSSE表示校验统计量;σγ表示误警检测值;β表示阈值因子;n表示参与定 位的卫星数量。
优选的,上述卫星故障识别方法,还包括排除故障卫星的步骤,具体为:
(1)根据预置的误报警率计算检测门限值;
TD=Φ-1(1-PFA/2n)
其中,TD表示检测门限值;PFA表示误报警率;Φ-1表示正态分布的标准值;n表 示参与定位的卫星数量;
(2)若rimax=max(|PT·Qp(:,i)|/(σ0·norm(Qp(:,i))))>TD,表明第i颗卫星的门限值 大于检测门限值,则判定第i颗卫星存在故障并将其剔除;否则判定参与定位的卫星中 不存在故障卫星。
优选的,上述卫星故障识别方法,所述相似径向误差保护值为最大斜率值与误警检 测值的乘积。
优选的,上述卫星故障识别方法,所述预设值为中国地区飞机非精密进场值,当 n<=6,ARPm(n)=337;当n=7,ARPm(n)=343;当n>=7,ARPm(n)=347;
其中,ARPm(n)表示飞机非精密进场值;n表示参与定位的卫星数量。
优选的,上述卫星故障识别方法,所述校验统计量的计算方法为:
FSSE=PTP,P=QPε;
其中,FSSE表示校验统计量;P表示奇偶空间矢量,ε表示伪距噪声矢量。
按照本发明的第二个方面,还提供了一种卫星故障识别装置,包括:
创建单元,用于采集基于最小二乘算法进行卫星定位时收敛的接收机位置数据和接 收机钟差,根据所述接收机位置数据、接收机钟差构建观测矩阵并对所述观测矩阵进行正交三角分解,得到标准正交矩阵;采用奇偶空间算法对所述标准正交矩阵进行分解, 生成奇偶空间矩阵;
第一计算单元,用于对所述奇偶空间矩阵进行矩阵运算,得到参与定位的每个卫星 的斜率并提取最大斜率值;
第一计算单元,用于根据所述最大斜率值和误警检测值计算相似径向误差保护值, 当所述相似径向误差保护值不小于预设值时,根据所述奇偶空间矩阵对应的卫星偏差数 据计算校验统计量;
检测单元,用于根据所述校验统计量、误警检测值和预置的阈值因子检测参与定位 的卫星中是否存在故障。
按照本发明的第三个方面,还提供了一种软件接收机,其包括至少一个处理单元、以及至少一个存储单元,
其中,所述存储单元存储有计算机程序,当所述程序被所述处理单元执行时,使得所述处理单元执行上述任一项卫星故障识别方法的步骤。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,能够取得下列有益 效果:
(1)本发明提供的卫星故障识别方法、装置及软件接收机,通过结合最小二乘残差算法和奇偶空间校验算法的各自优点,使得新算法不仅与传统最小二乘算法,或者奇 偶空间校验算法相比,具有时间复杂度低、空间复杂度低的特点,并且具有相同的卫星 伪距跳变识别度,使其能够适应常态飞行、高动态飞行的卫星故障检测和识别。
(2)本发明提供的卫星故障识别方法、装置及软件接收机,在相同的误警率和漏警率的前提下,通过在卫星故障检测和识别时设置阈值因子,增大了卫星的最小伪距跳 变识别度。
附图说明
图1是本发明提供的卫星故障识别方法的流程图;
图2是本发明实施例提供的卫星故障识别方法的算法流程图;
图3是本发明实施例提供的软件接收机通过接收卫星信号模拟器发射的卫星信号, 通过模拟器上位机,控制GPS卫星均伪距不跳变的定位效果图;
图4是本发明实施例提供的软件接收机通过接收卫星信号模拟器发射的卫星信号, 通过模拟器上位机,控制BD卫星均伪距不跳变的定位效果图;
图5是本发明实施例提供的软件接收机通过接收卫星信号模拟器发射的卫星信号, 通过模拟器上位机,控制GPS+BD卫星均伪距不跳变的组合定位效果图;
图6是本发明实施例提供的软件接收机通过接收卫星信号模拟器发射的卫星信号, 通过模拟器上位机,控制GPS的27#卫星伪距跳变45m,使用卫星故障识别算法剔除故障卫星的效果图。
图7是本发明实施例提供的软件接收机通过接收卫星信号模拟器发射的卫星信号, 通过模拟器上位机,控制BD的103#卫星伪距跳变45m,使用卫星故障识别算法剔除故障卫星的效果图。
图8是本发明实施例提供的软件接收机通过接收卫星信号模拟器发射的卫星信号, 通过模拟器上位机,控制GPS+BD卫星组合定位时,GPS伪距跳变45m,使用卫星故 障识别算法剔除故障卫星的效果图。
图9是本发明实施例提供的软件接收机通过接收卫星信号模拟器发射的卫星信号, 通过模拟器上位机,控制GPS+BD卫星组合定位时,BD的103#卫星伪距跳变45m,使 用卫星故障识别算法剔除故障卫星的效果图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征 只要彼此之间未构成冲突就可以相互组合。
本发明所提供的一种卫星故障识别方法,联合最小二乘残差算法和奇偶空间校验算 法,减少算法的时间复杂度和空间复杂度;另外通过设置阈值因子,使识别算法具有较高的故障卫星伪距跳变识别度。
图1是本发明提供的卫星故障识别方法的流程图,参见图1所示,该方法包括以下步骤:
S1:采集基于最小二乘算法进行卫星定位时收敛的接收机位置数据和接收机钟差, 根据所述接收机位置数据、接收机钟差构建观测矩阵H并对所述观测矩阵H进行正交三角分解,得到标准正交矩阵Q;
假设某时间内有n颗卫星参与定位,若n<5,则不启动RAIM算法;若n>5,则 启动RAIM算法,同时构建采用最小二乘算法进行定位时收敛的观测矩阵H;对观测矩 阵H进行正交三角分解(QR分解),得到标准正交矩阵Q和上三角矩阵R。
S2:采用奇偶空间算法对所述标准正交矩阵Q进行分解,生成奇偶空间矩阵Qp;
本实施例中,提取标准正交矩阵Q的转置矩阵QT的5~n行数据,生成奇偶空间矩阵Qp。
S3:根据所述奇偶空间矩阵Qp计算参与定位的每个卫星的斜率并提取最大斜率值Slopemax;具体包括:
(1)根据奇偶空间矩阵Qp构造第一矩阵B,使得B=In-Qp TQp;然后构造第二矩 阵A,使得A=H\B;
(2)根据第一矩阵B和第二矩阵A计算参与定位的每个卫星的斜率并取其最大值Slopemax;
在软件接收机中,最难检测的参与卫星定位的故障卫星就是其能够产生最大斜率的 卫星。因此,本发明采用相似径向误差保护值ARP作为定位卫星在空间的几何分布作为衡量标准;首先计算参与定位的每颗卫星的斜率,并取其最大值Slopemax;
S4:根据最大斜率值Slopemax和误警检测值σγ计算相似径向误差保护值ARP,当相似径向误差保护值ARP不小于预设值时,根据奇偶空间矩阵Qp计算校验统计量FSSE;
其中,icdf('ch2',1-PFA,n-4)表示求取以(n-4)为自由度,累计概率为(1-PFA)的χ2分布的临界值。
相似径向误差保护值ARP为最大斜率值Slopemax与误警检测值σγ的乘积,即ARP=Slopemax·σγ;
然后进行RAIM算法可用性判断,若ARP=Slopemax·σr≥ARPm(n),则RAIM算法 可用;否则不可用。其中ARPm(n)是中国地区飞机非精密进场值,当n<=6,ARPm(n)=337。 当n=7,ARPm(n)=343。当n>=7,ARPm(n)=347;
当判定RAIM算法可用时,根据奇偶空间矩阵Qp计算校验统计量FSSE,具体的:
FSSE=PTP,P=QPε;
其中,P表示奇偶空间矢量,表示卫星的偏差数据;ε表示伪距噪声矢量。
S5:根据校验统计量FSSE、误警检测值σγ和预置的阈值因子检测参与定位的卫星中是否存在故障具体包括:
S6:排除故障卫星;
根据预置的误报警率PFA计算检测门限值TD;若第i颗卫星的门限值大于检测门限值TD,则判定第i颗卫星存在故障并将其剔除;否则判定参与定位的卫星中不存在故障 卫星;具体的,若rimax=max(|PT·Qp(:,i)|/(σ0·norm(Qp(:,i))))>TD,则表明第i颗卫星 存在故障并将其剔除;
其中,TD=Φ-1(1-PFA/2n),Φ-1表示正态分布的标准值;
考虑到软件接收机快速定位的需求,通常将可见卫星数n存储在数组中,接收机DSP程序在烧写到FLASH之前,就可计算出误警检测值数组σr'[n-4]、检测门限数组 TD'[n-4]中,有效避免重复计算,提高定位速度。
下面结合附图和具体的实例对本发明提供的卫星故障识别方法作更进一步的详细 说明。
实施例一
图2较为详细地描述了本发明提出的卫星故障识别方法的具体的算法流程,参见图 2所示,具体流程如下:
(1)软件接收机通过天线接收卫星信号后,首先启动接收机内FLASH中的存储程序,实现对卫星信号的位同步、帧同步后,解析有效卫星星历信息,进行卫星导航电文 译码,最后获取卫星编号、卫星伪距等信息。通过上述信息计算出卫星的俯仰角、方位 角,并删除俯仰角小于10度的卫星。
(2)待选取有效卫星后,采用最小二乘算法进行粗略定位后,选用收敛后的接收机位置x、y、z以及接收机钟差dt构建观测矩阵H;软件接收机上电一段时间,得到某一 时刻的观测矩阵H如下所示:
(3)从观测矩阵H可知,当前参与定位的卫星数n=7,启动卫星故障识别算法,同时对观测矩阵H进行QR分解,得到Q、R矩阵,获取奇偶空间矩阵Qp;
(4)构造矩阵B,使得B=In-Qp TQp;然后构造矩阵A,使得A=H\B;
TD=Φ-1(1-PFA/2n)
其中,n=7;
(7)卫星故障识别算法可用性判断:ARP=Slopemax·σr≥ARPm(7)=343,则卫星故障识别算法可用;
(8)根据步骤(3)中的奇偶空间矩阵Qp,计算校验统计量FSSE=PTP;其中,P代表 奇偶空间矢量,该奇偶空间矢量表示卫星的偏差数据;
(10)故障排除:rimax=max(|PT·Qp(:,i)|/(σ0·norm(Qp(:,i))))>TD,则说明第1颗卫 星存在故障。
(11)在剔除参与定位的故障卫星后,采用最小二乘算法重新利用剩余卫星进行定位, 输出精确的用户位置和用户钟差;
(12)采用VS2010编写MFC程序,将接收机的位置x、y、z,速度信息vx、vy、vz, UTC时间,参与定位卫星的伪距、伪距率、速度、跟踪状态、多普勒计数等有用信息输 出显示给用户。
图6~图9中的实线框表示在满足高动态场景下,卫星伪距发生跳变时,使用本发明 提供的卫星故障识别方法成功将其剔除,并且给出0X42F8的标志;状态掩码0X42F8 即表示识别出此颗卫星是故障卫星。
实施例二
本实施例提供了一种卫星故障识别装置,该装置包括:
创建单元,用于采集基于最小二乘算法进行卫星定位时收敛的接收机位置数据和接 收机钟差,根据上述接收机位置数据、接收机钟差构建观测矩阵H并对该观测矩阵H进行正交三角分解,得到标准正交矩阵Q;采用奇偶空间算法对上述标准正交矩阵Q进行 分解,生成奇偶空间矩阵Qp;
第一计算单元,用于对上述奇偶空间矩阵Qp进行矩阵运算,得到参与定位的每个卫星的斜率并提取最大斜率值Slopemax;
第一计算单元,用于根据上述最大斜率值Slopemax和误警检测值σγ计算相似径向误 差保护值ARP,当该相似径向误差保护值ARP不小于预设值时,根据奇偶空间矩阵Qp对应的卫星偏差数据计算校验统计量FSSE;
检测单元,用于根据上述校验统计量FSSE、误警检测值σγ和预置的阈值因子检测参与定位的卫星中是否存在故障。
关于卫星故障识别装置及其各功能单元的具体限定可以参见上文中对于卫星故障 识别方法的限定,在此不再赘述。上述卫星故障识别装置的各个功能单元可全部或部分通过软件、硬件及其组合来实现。上述各功能单元可以以硬件形式内嵌于或独立于计算 机装置的处理器中,也可以以软件形式存储于计算机装置的存储器中,以便于处理器调 用执行以上各个单元对应的操作。
实施例三
本实施例提供了一种软件接收机,其包括至少一个处理器、以及至少一个存储器,其中,存储器中存储有计算机程序,当计算机程序被处理器执行时,使得处理器执行上 述卫星故障识别方法的步骤。处理器和存储器的类型不作具体限制,例如:处理器可以 是微处理器、数字信息处理器、片上可编程逻辑系统等;存储器可以是易失性存储器、 非易失性存储器或者它们的组合等。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限 制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种卫星故障识别方法,其特征在于,包括以下步骤:
S1:采集基于最小二乘算法进行卫星定位时收敛的接收机位置数据和接收机钟差,根据所述接收机位置数据、接收机钟差构建观测矩阵并采用奇偶空间算法对所述观测矩阵进行分解,生成奇偶空间矩阵;
S2:对所述奇偶空间矩阵进行矩阵运算,得到参与定位的每个卫星的斜率并提取最大斜率值;
S3:基于所述最大斜率值和误警检测值计算定位卫星的相似径向误差保护值,当所述相似径向误差保护值不小于预设值时,根据所述奇偶空间矩阵对应的卫星偏差数据计算校验统计量;
S4:根据所述校验统计量、误警检测值和预置的阈值因子检测参与定位的卫星中是否存在故障。
5.如权利要求1所述的卫星故障识别方法,其特征在于,还包括排除故障卫星的步骤,具体为:
根据预置的误报警率计算检测门限值,若第i颗卫星的门限值大于检测门限值,则判定第i颗卫星存在故障并将其剔除;
TD=Φ-1(1-PFA/2n)
其中,TD表示检测门限值;PFA表示误报警率;Φ-1表示正态分布的标准值;n表示参与定位的卫星数量。
6.如权利要求1所述的卫星故障识别方法,其特征在于,所述相似径向误差保护值为最大斜率值与误警检测值的乘积。
7.如权利要求1所述的卫星故障识别方法,其特征在于,所述预设值为中国地区飞机非精密进场值,当n<=6,ARPm(n)=337;当n=7,ARPm(n)=343;当n>=7,ARPm(n)=347;
其中,ARPm(n)表示飞机非精密进场值;n表示参与定位的卫星数量。
8.如权利要求1所述的卫星故障识别方法,其特征在于,所述校验统计量的计算方法为:FSSE=PTP,P=QPε;
其中,FSSE表示校验统计量;P表示奇偶空间矢量,ε表示伪距噪声矢量。
9.一种卫星故障识别装置,其特征在于,包括:
创建单元,用于采集基于最小二乘算法进行卫星定位时收敛的接收机位置数据和接收机钟差,根据所述接收机位置数据、接收机钟差构建观测矩阵并采用奇偶空间算法对所述观测矩阵进行分解,生成奇偶空间矩阵;
第一计算单元,用于对所述奇偶空间矩阵进行矩阵运算,得到参与定位的每个卫星的斜率并提取最大斜率值;
第一计算单元,用于根据所述最大斜率值和误警检测值计算定位卫星的相似径向误差保护值,当所述相似径向误差保护值不小于预设值时,根据所述奇偶空间矩阵对应的卫星偏差数据计算校验统计量;
检测单元,用于根据所述校验统计量、误警检测值和预置的阈值因子检测参与定位的卫星中是否存在故障。
10.一种软件接收机,其特征在于,包括至少一个处理单元、以及至少一个存储单元,
其中,所述存储单元存储有计算机程序,当所述程序被所述处理单元执行时,使得所述处理单元执行权利要求1~8任一项所述方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910996032.4A CN110907953B (zh) | 2019-10-18 | 2019-10-18 | 一种卫星故障识别方法、装置及软件接收机 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910996032.4A CN110907953B (zh) | 2019-10-18 | 2019-10-18 | 一种卫星故障识别方法、装置及软件接收机 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110907953A true CN110907953A (zh) | 2020-03-24 |
CN110907953B CN110907953B (zh) | 2022-04-29 |
Family
ID=69814609
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910996032.4A Active CN110907953B (zh) | 2019-10-18 | 2019-10-18 | 一种卫星故障识别方法、装置及软件接收机 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110907953B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114114351A (zh) * | 2022-01-24 | 2022-03-01 | 浙江时空道宇科技有限公司 | 卫星导航测量故障的排除方法、装置、设备及存储介质 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070279279A1 (en) * | 2006-06-06 | 2007-12-06 | Honeywell International Inc. | System and method for enhancing the performance of satellite navigation receivers |
US20100033370A1 (en) * | 2003-12-02 | 2010-02-11 | Gmv Aerospace And Defence, S.A. | Gnss navigation solution integrity in non-controlled environments |
CN101806903A (zh) * | 2010-04-01 | 2010-08-18 | 北京航空航天大学 | 一种用于卫星导航系统的接收机自主完好性监测方法 |
EP2402785A1 (en) * | 2010-06-23 | 2012-01-04 | Astrium GmbH | An improved RAIM algorithm |
CA2751939A1 (en) * | 2010-09-10 | 2012-03-10 | Thales | Method for correcting position estimations by selecting pseudo-distance measurements |
CN102901971A (zh) * | 2012-09-27 | 2013-01-30 | 电子科技大学 | 基于奇偶矢量法的双星故障识别方法 |
CN103592657A (zh) * | 2013-09-30 | 2014-02-19 | 北京大学 | 一种基于钟差辅助的低可见星下单模raim实现方法 |
CN103592656A (zh) * | 2013-10-17 | 2014-02-19 | 航天恒星科技有限公司 | 一种适用于星载导航接收机的自主完好性监测方法 |
CN106707304A (zh) * | 2016-12-23 | 2017-05-24 | 天津七六四通信导航技术有限公司 | 一种卫星导航接收机故障星检测方法 |
-
2019
- 2019-10-18 CN CN201910996032.4A patent/CN110907953B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100033370A1 (en) * | 2003-12-02 | 2010-02-11 | Gmv Aerospace And Defence, S.A. | Gnss navigation solution integrity in non-controlled environments |
US20070279279A1 (en) * | 2006-06-06 | 2007-12-06 | Honeywell International Inc. | System and method for enhancing the performance of satellite navigation receivers |
CN101806903A (zh) * | 2010-04-01 | 2010-08-18 | 北京航空航天大学 | 一种用于卫星导航系统的接收机自主完好性监测方法 |
EP2402785A1 (en) * | 2010-06-23 | 2012-01-04 | Astrium GmbH | An improved RAIM algorithm |
CA2751939A1 (en) * | 2010-09-10 | 2012-03-10 | Thales | Method for correcting position estimations by selecting pseudo-distance measurements |
CN102901971A (zh) * | 2012-09-27 | 2013-01-30 | 电子科技大学 | 基于奇偶矢量法的双星故障识别方法 |
CN103592657A (zh) * | 2013-09-30 | 2014-02-19 | 北京大学 | 一种基于钟差辅助的低可见星下单模raim实现方法 |
CN103592656A (zh) * | 2013-10-17 | 2014-02-19 | 航天恒星科技有限公司 | 一种适用于星载导航接收机的自主完好性监测方法 |
CN106707304A (zh) * | 2016-12-23 | 2017-05-24 | 天津七六四通信导航技术有限公司 | 一种卫星导航接收机故障星检测方法 |
Non-Patent Citations (2)
Title |
---|
KEJING CAO等: ""Research On Improved RAIM Algorithm Based On Parity Vector Method"", 《2013 INTERNATIONAL CONFERENCE ON INFORMATION TECHNOLOGY AND APPLICATIONS》 * |
李靖宇: "组合导航系统完好性技术研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114114351A (zh) * | 2022-01-24 | 2022-03-01 | 浙江时空道宇科技有限公司 | 卫星导航测量故障的排除方法、装置、设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN110907953B (zh) | 2022-04-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11237276B2 (en) | System and method for gaussian process enhanced GNSS corrections generation | |
KR102479956B1 (ko) | 전역 항법 위성 시스템을 위한 다중 위성 측정 실패 검출 및 분리의 방법 | |
CN110007317B (zh) | 一种选星优化的高级接收机自主完好性监测方法 | |
US6166683A (en) | System and method for high-integrity detection and correction of cycle slip in a carrier phase-related system | |
CN111060133B (zh) | 一种用于城市复杂环境的组合导航完好性监测方法 | |
US9983314B2 (en) | System for excluding a failure of a satellite in a GNSS system | |
US20240142637A1 (en) | System and method for gaussian process enhanced gnss corrections generation | |
CN109085619B (zh) | 多模gnss系统的定位方法及装置、存储介质、接收机 | |
CN110927757B (zh) | 卫星观测数据的质控方法、装置及定位装置 | |
FR2963436A1 (fr) | Procede de determination d'un volume de protection dans le cas de deux pannes satellitaires simultanees | |
CN109307876B (zh) | 一种适用于gnss矢量跟踪的自主完好性监测方法 | |
CN111708054B (zh) | 基于粒子群优化算法的araim垂直保护级优化方法 | |
CN108088470A (zh) | 一种组合导航慢变斜坡故障完好性监测方法 | |
JIANG et al. | A robust fault detection algorithm for the GNSS/INS integrated navigation systems | |
CN105511481B (zh) | 一种星载定轨优化方法 | |
CN110907953B (zh) | 一种卫星故障识别方法、装置及软件接收机 | |
CN111913195B (zh) | 基于陆基无线电导航信息的gps接收机抗欺骗干扰处理方法 | |
CN116859415A (zh) | 一种快速稳健高精度的多故障卫星识别及定位方法 | |
CN112083463A (zh) | 模糊度是否固定正确的检测方法及装置、定位终端 | |
CN115561782A (zh) | 一种基于奇偶矢量投影的组合导航中卫星故障检测方法 | |
CN115062097A (zh) | 一种定位完好性监测方法、终端及计算机存储介质 | |
CN110146900B (zh) | 一种卫星导航低复杂度遍历式raim防欺骗方法及装置 | |
US20210356602A1 (en) | Single delta range differences using synthetic clock steering | |
CN115390099B (zh) | 基于奇偶矢量投影的最大值联合卡方故障排除方法 | |
CN113093248A (zh) | 一种范围一致性raim方法 |
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 |