CN110906929A - 一种卫星异常数据识别方法及装置 - Google Patents

一种卫星异常数据识别方法及装置 Download PDF

Info

Publication number
CN110906929A
CN110906929A CN201911272942.4A CN201911272942A CN110906929A CN 110906929 A CN110906929 A CN 110906929A CN 201911272942 A CN201911272942 A CN 201911272942A CN 110906929 A CN110906929 A CN 110906929A
Authority
CN
China
Prior art keywords
data
satellite
abnormal
speed
columns
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
Application number
CN201911272942.4A
Other languages
English (en)
Other versions
CN110906929B (zh
Inventor
张霞
廉洁
马玉海
吴炜平
袁春贵
刘凯
杨毅强
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhongke Aerospace (Guangzhou) Equipment Industry Co.,Ltd.
Original Assignee
Beijing Zhongke Aerospace Exploration Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Zhongke Aerospace Exploration Technology Co Ltd filed Critical Beijing Zhongke Aerospace Exploration Technology Co Ltd
Priority to CN201911272942.4A priority Critical patent/CN110906929B/zh
Publication of CN110906929A publication Critical patent/CN110906929A/zh
Application granted granted Critical
Publication of CN110906929B publication Critical patent/CN110906929B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/10Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
    • G01C21/12Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
    • G01C21/16Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
    • G01C21/165Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation combined with non-inertial navigation instruments
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/20Instruments for performing navigational calculations
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/24Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for cosmonautical navigation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42Determining position
    • G01S19/45Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement
    • G01S19/47Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement the supplementary measurement being an inertial measurement, e.g. tightly coupled inertial

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Astronomy & Astrophysics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Navigation (AREA)

Abstract

本申请公开一种卫星异常数据识别方法及装置。所述方法包括获取惯导系统解算出的速度和位置,以及获取卫星定位系统输出的速度和位置;根据惯导系统解算出的速度和位置以及卫星定位系统输出速度和位置得到量测变量;根据量测变量预测新息数据,并由新息数据流计算均值和标准差,然后计算Grubbs检验统计量;根据设定的检出水平对检验统计量进行约束;若经检出水平约束后的检验统计量为异常值则剔除卫星定位系统数据,否则输入卡尔曼滤波器进行滤波估计。本申请可实现实时在线检测,具有简捷的递推关系,计算量小不会造成系统的额外负担;基于历史统计数据,具有一定弹性,能克服异常数据的不利影响,也使判别系统在异常数据的影响下快速恢复,避免死锁,提高卡尔曼滤波估计的可靠性。

Description

一种卫星异常数据识别方法及装置
技术领域
本申请涉及信息技术领域,尤其涉及一种卫星异常数据识别方法及装置。
背景技术
现有的工程上卫星异常数据判定方法基本采用测量范围判定,结合卫星的定位信息,如定位标志、DOP值等做层层状态判断。如《箭载GPS异常数据检测与剔除》,针对箭载GPS原始测量数据的不同特性,结合时间序列差分及门限检测等算法,给出了测量数据的状态位检测、时标错位检测、CodeTime检测和多普勒频率检测等异常数据检测与剔除算法;或者根据观测值的统计特性,通过计算标准差等来设置阈值对野值进行判断和剔除,通过某种信息处理和特征提取方法;或者采用自适应滤波算法和其它滤波算法相结合的算法,例如《结合自适应滤波和神经网络的GNSS/INS抗差组合导航算法》结合了自适应滤波和神经网络,提出了两种GNSS/INS抗差自适应组合导航解算方案。
对GNSS(全球导航卫星系统定位)接收机输出信息进行层层判别的方法是现行进行GNSS异常数据判别中较为常见的方法,但由于市场上的GNSS接收机型号众多,本身的状态信息输出难以统一和严格准确,且判断过程纷繁复杂,容易导致误判;计算标准差来设置阈值也需要设计者掌握一些先验信息,取得足够的经验才能确保所设计的阈值可靠安全;采用自适应滤波和其它滤波算法相结合的方法虽然适应性较强,但是算法繁杂,造成了系统的额外负担。
发明内容
本申请提供了一种卫星异常数据的识别方法,包括:
获取惯导系统解算出的速度和位置,以及获取卫星定位系统输出的速度和位置;
根据惯导系统解算出的速度和位置以及卫星定位系统输出速度和位置得到量测变量;
根据量测变量预测新息数据,并由预测出的新息数据流计算均值和标准差,然后根据均值和标准差计算Grubbs检验统计量;
设定检出水平,根据检出水平对Grubbs检验统计量进行约束;
判断经检出水平约束后的Grubbs检验统计量是否为异常值,如果是,则剔除该异常值,否则输入卡尔曼滤波器进行卡尔曼滤波估计。
如上所述的卫星异常数据的识别方法,其中将惯导系统解算出的速度和位置,与卫星定位系统输出的速度和位置作差值,将差值作为量测变量。
如上所述的卫星异常数据的识别方法,其中通过以下公式预测新息数据:
Figure BDA0002314712210000021
其中,k代表当前滤波周期,k-1代表上一滤波周期,ak代表预测的新息数据,Zk代表当前滤波周期的系统观测矢量,Hk代表当前滤波周期的系统观测矩阵,
Figure BDA0002314712210000022
代表上一周期对本周期的一步预测均值。
如上所述的卫星异常数据的识别方法,其中发射惯性坐标系下惯导系统/卫星定位系统的松组合导航系统的观测方程为
Figure BDA0002314712210000023
其中,Z为系统观测矢量;H为系统观测矩阵;v为系统观测噪声阵;Zv为发射惯性坐标系下速度观测方程,Zp为发射惯性坐标系下位置观测方程。
如上所述的卫星异常数据的识别方法,其中Zv为发射惯性坐标系下速度观测方程:
Figure BDA0002314712210000024
其中,
Figure BDA0002314712210000025
即Hv为三行十五列的矩阵,其中前三列为0,接下来三列为对角矩阵,后九列为0;vv=[δVgx δVgy δVgz]T为GNSS接收机速度测量误差。
如上所述的卫星异常数据的识别方法,其中Zp为发射惯性坐标系下位置观测方程:
Figure BDA0002314712210000026
其中,
Figure BDA0002314712210000027
即Hp为三行十五列的矩阵,其中前六列为0,接下来三列为对角矩阵,后六列为0;vp=[δXg δYg δZg]T为GNSS接收机速度测量误差。
如上所述的卫星异常数据的识别方法,其中通过对每个滤波周期的新息计算,可得到若干新息数据流为a1,a2…ak,取最新测量的n组新息数据ak-n+1,ak-n+2…ak,其均值(或中位值)为
Figure BDA0002314712210000031
其标准差为σ,计算Grubbs检验统计量
Figure BDA0002314712210000032
如上所述的卫星异常数据的识别方法,其中根据GNSS的滚转角速度和加加速度的适应范围和要求设定检出水平α,然后使用检出水平α对Grubbs检验统计量进行约束。
如上所述的卫星异常数据的识别方法,其中通过查询格拉布斯临界系统表,判断经检出水平约束后的Grubbs检验统计量是否超过格拉布斯临界值,如果是,则该Grubbs检验统计量为异常值,否则该Grubbs检验统计量为正常值。
本申请还提供一种卫星异常数据识别装置,所述装置执行上述卫星异常数据的识别方法。
本申请实现的有益效果如下:
(1)结合卡尔曼滤波观测值中所包含的信息,可实现实时在线检测,具有简捷的递推关系,计算量小不会造成系统的额外负担;
(2)具有良好的容错能力,该异常识别算法基于历史统计数据,具有一定弹性,不仅能使滤波器克服异常数据的不利影响,也使判别系统在异常数据的影响下快速恢复,避免死锁,提高卡尔曼滤波估计的可靠性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明中记载的一些实施例,对于本领域普通技术人员来讲,还可以根据这些附图获得其他的附图。
图1是本申请实施例提供的一种卫星异常数据的识别方法流程图;
图2为模拟弹道GNSS速度曲线;
图3为GNSS速度放大图;
图4为正确识别异常数据的示意图;
图5为本申请实施例提供的一种卫星异常数据的识别装置示意图。
具体实施方式
下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例一
本申请实施例一提供一种卫星异常数据的识别方法,如图1所示,包括:
步骤110、获取惯导系统解算出的速度和位置,以及获取卫星定位系统输出的速度和位置;
本申请提供的卫星异常数据识别方法适用于组合导航结算中,组合导航包括惯导系统(SINS,即惯性导航系统,是一种不依赖外部信息也不向外部辐射能量的自主式导航系统)和卫星定位系统(GNSS,即全球导航卫星定位系统,是一种能在地球表面或近地空间的任何地点为用户提供全天候的三维坐标和速度以及时间信息的空基无线电导航定位系统);
本申请中,惯导系统解算出的速度
Figure BDA0002314712210000041
位置为
Figure BDA0002314712210000042
卫星定位系统输出的速度
Figure BDA0002314712210000043
位置为
Figure BDA0002314712210000044
步骤120、根据惯导系统解算出的速度和位置以及卫星定位系统输出速度和位置得到量测变量;
具体地,将惯导系统解算出的速度和位置,与卫星定位系统输出的速度和位置作差值,将差值作为量测变量,然后对此次的量测变量进行步骤130-150的评估操作,验证本次量测变量是否为异常值,以此减少卫星定位的误判率;
其中,发射惯性坐标系下惯导系统/卫星定位系统的松组合导航系统的观测方程由速度观测方程、位置观测方程组成:
发射惯性坐标系下速度观测方程如下:
Figure BDA0002314712210000045
其中,
Figure BDA0002314712210000046
即Hv为三行十五列的矩阵,其中前三列为0,接下来三列为对角矩阵,后九列为0;vv=[δVgx δVgy δVgz]T为GNSS接收机速度测量误差;
发射惯性坐标系下位置观测方程如下:
Figure BDA0002314712210000051
其中,
Figure BDA0002314712210000052
即Hp为三行十五列的矩阵,其中前六列为0,接下来三列为对角矩阵,后六列为0;vp=[δXg δYg δZg]T为GNSS接收机速度测量误差;
由此可得发射惯性坐标系下惯导系统/卫星定位系统的松组合导航系统的观测方程如下:
Figure BDA0002314712210000053
其中,Z为系统观测矢量;H为系统观测矩阵;v为系统观测噪声阵。
步骤130、根据量测变量预测新息数据,新息数据为均方误差计算;
具体地,通过以下公式预测新息数据:
Figure BDA0002314712210000054
其中,k代表当前滤波周期,k-1代表上一滤波周期,ak代表预测的新息数据,Zk代表当前滤波周期的系统观测矢量,Hk代表当前滤波周期的系统观测矩阵,
Figure BDA0002314712210000055
代表上一周期对本周期的一步预测均值。
需要说明的是,步骤110和步骤120为惯导系统/卫星定位系统组合导航计算中的基本步骤,无需占用额外计算量。
步骤140、根据预测出的新息数据流计算均值和标准差,然后根据均值和标准差计算Grubbs检验统计量;
通过对每个滤波周期的新息计算,可得到若干新息数据流为a1,a2…ak,取最新测量的n组新息数据ak-n+1,ak-n+2…ak,其均值(或中位值)为
Figure BDA0002314712210000056
其标准差为σ,计算Grubbs(格拉布斯检验法,是一种统计检验方法)检验统计量
Figure BDA0002314712210000057
步骤150、设定检出水平,根据检出水平对Grubbs检验统计量进行约束;
具体地,根据GNSS的滚转角速度和加加速度的适应范围和要求设定检出水平α,优选地,由于越接近滚转角速度和加加速度边界则表示要求越严格,所以将检出水平α设定的小一些;然后使用检出水平α对Grubbs检验统计量进行约束。
步骤160、查询格拉布斯临界系统表,判断经检出水平约束后的Grubbs检验统计量是否超过格拉布斯临界值,如果是,则确定该测量值为异常值,剔除该异常值,否则将该测量值作为正常数据输入卡尔曼滤波器进行卡尔曼滤波估计;
其中,格拉布斯临界系统表如下表所示:
Figure BDA0002314712210000061
表1
具体地,通过查询格拉布斯临界系统表可以确定Grubbs检验统计量是否为异常值,如果检验出该统计量为异常值,则剔除该异常GNSS数据,防止异常数据进入卡尔曼滤波器量测更新过程,防止滤波器发散。
例如,图2为模拟弹道GNSS速度曲线,图3为GNSS速度放大图。从图上可知,速度Vy方向在500、1000、1500、2000数据点上存在异常值,正常GNSS速度误差在0.1m/s(即一倍的检出水平σ),而图2中的异常值偏离正常值约为1.2m/s,所以为了防止异常数据进入卡尔曼滤波器会出现滤波器发散的问题,通过本申请可成功实现异常值判别,由此剔除异常值在很大程度上减小误判率。图4为正确识别异常数据的示意图,从图中可以看出位于点数500、1000、1500、2000的GNSS速度异常值都被剔除出来,实现了异常值100%判别;对无异常值的数据判别出现了一次误判,误判率低于0.5‰,在卫星数据的异常判别中,对正常数据仅一次误判只会使滤波少一次量测更新过程,基本不影响整体的滤波效果,其影响可忽略。
此外,本申请还提供一种卫星异常数据识别装置,如图5所示,所述装置包括IMU(惯性测量)模块、导航解算模块、卫星定位模块、新息计算模块、检出水平计算模块、卡尔曼滤波模块和格拉布斯判别模块;所述装置执行上述卫星异常数据的识别方法。
以上所述实施例,仅为本申请的具体实施方式,用以说明本申请的技术方案,而非对其限制,本申请的保护范围并不局限于此,尽管参照前述实施例对本申请进行了详细的说明,本领域的普通技术人员应当理解:任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化,或者对其中部分技术特征进行等同替换;而这些修改、变化或者替换,并不使相应技术方案的本质脱离本申请实施例技术方案的精神和范围。都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应所述以权利要求的保护范围为准。

Claims (10)

1.一种卫星异常数据的识别方法,其特征在于:
获取惯导系统解算出的速度和位置,以及获取卫星定位系统输出的速度和位置;
根据惯导系统解算出的速度和位置以及卫星定位系统输出速度和位置得到量测变量;
根据量测变量预测新息数据,并由预测出的新息数据流计算均值和标准差,然后根据均值和标准差计算Grubbs检验统计量;
设定检出水平,根据检出水平对Grubbs检验统计量进行约束;
判断经检出水平约束后的Grubbs检验统计量是否为异常值,如果是,判定卫星定位系统数据为异常值,进行剔除,否则输入卡尔曼滤波器进行卡尔曼滤波估计。
2.如权利要求1所述的卫星异常数据的识别方法,其特征在于,将惯导系统解算出的速度和位置,与卫星定位系统输出的速度和位置作差值,将差值作为量测变量。
3.如权利要求1所述的卫星异常数据的识别方法,其特征在于,通过以下公式预测新息数据:
Figure FDA0002314712200000011
其中,k代表当前滤波周期,k-1代表上一滤波周期,ak代表预测的新息数据,Zk代表当前滤波周期的系统观测矢量,Hk代表当前滤波周期的系统观测矩阵,
Figure FDA0002314712200000012
代表上一周期对本周期的一步预测均值。
4.如权利要求3所述的卫星异常数据的识别方法,其特征在于,发射惯性坐标系下惯导系统/卫星定位系统的松组合导航系统的观测方程为
Figure FDA0002314712200000013
其中,Z为系统观测矢量;H为系统观测矩阵;v为系统观测噪声阵;Zv为发射惯性坐标系下速度观测方程,Zp为发射惯性坐标系下位置观测方程。
5.如权利要求4所述的卫星异常数据的识别方法,其特征在于,Zv为发射惯性坐标系下速度观测方程:
Figure FDA0002314712200000014
其中,Hv=[03×3 diag[1 1 1]03×9]3×15,即Hv为三行十五列的矩阵,其中前三列为0,接下来三列为对角矩阵,后九列为0;vv=[δVgx δVgy δVgz]T为GNSS接收机速度测量误差。
6.如权利要求4所述的卫星异常数据的识别方法,其特征在于,Zp为发射惯性坐标系下位置观测方程:
Figure FDA0002314712200000021
其中,Hp=[03×6 diag[1 1 1]03×6]3×15,即Hp为三行十五列的矩阵,其中前六列为0,接下来三列为对角矩阵,后六列为0;vp=[δXg δYg δZg]T为GNSS接收机速度测量误差。
7.如权利要求1所述的卫星异常数据的识别方法,其特征在于,通过对每个滤波周期的新息计算,可得到若干新息数据流为a1,a2…ak,取最新测量的n组新息数据ak-n+1,ak-n+2…ak,其均值(或中位值)为
Figure FDA0002314712200000022
其标准差为σ,计算Grubbs检验统计量
Figure FDA0002314712200000023
8.如权利要求1所述的卫星异常数据的识别方法,其特征在于,根据GNSS的滚转角速度和加加速度的适应范围和要求设定检出水平α,然后使用检出水平α对Grubbs检验统计量进行约束。
9.如权利要求1所述的卫星异常数据的识别方法,其特征在于,通过查询格拉布斯临界系统表,判断经检出水平约束后的Grubbs检验统计量是否超过格拉布斯临界值,如果是,则该Grubbs检验统计量为异常值,否则该Grubbs检验统计量为正常值。
10.一种卫星异常数据识别装置,其特征在于,所述装置执行权利要求1-9任意一项所述的卫星异常数据的识别方法。
CN201911272942.4A 2019-12-12 2019-12-12 一种卫星异常数据识别方法及装置 Active CN110906929B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911272942.4A CN110906929B (zh) 2019-12-12 2019-12-12 一种卫星异常数据识别方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911272942.4A CN110906929B (zh) 2019-12-12 2019-12-12 一种卫星异常数据识别方法及装置

Publications (2)

Publication Number Publication Date
CN110906929A true CN110906929A (zh) 2020-03-24
CN110906929B CN110906929B (zh) 2022-01-25

Family

ID=69825112

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911272942.4A Active CN110906929B (zh) 2019-12-12 2019-12-12 一种卫星异常数据识别方法及装置

Country Status (1)

Country Link
CN (1) CN110906929B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021213525A1 (zh) * 2020-04-24 2021-10-28 北京智行者科技有限公司 自适应的gps错误观测值识别方法、存储介质、产品及系统
CN114485725A (zh) * 2021-12-22 2022-05-13 深圳元戎启行科技有限公司 数据异常检测方法、自动驾驶平台以及计算机可读存储介质

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103487820A (zh) * 2013-09-30 2014-01-01 东南大学 一种车载捷联/卫星紧组合无缝导航方法
US20140058696A1 (en) * 2012-08-24 2014-02-27 Samsung Electro-Mechanics Co., Ltd. Apparatus and method for controlling gain automatically in inertia sensor
CN107121684A (zh) * 2017-05-02 2017-09-01 北京航空航天大学 一种基于残差卡方检验法的gps诱骗识别和阈值决策方法
CN108205432A (zh) * 2016-12-16 2018-06-26 中国航天科工飞航技术研究院 一种观测实验数据异常值的实时剔除方法
CN109084761A (zh) * 2018-08-07 2018-12-25 合肥正阳光电科技有限责任公司 一种基于惯性寻北装置的抗冲击寻北解算方法
CN109471143A (zh) * 2018-12-11 2019-03-15 北京交通大学 自适应容错的列车组合定位方法
CN109724598A (zh) * 2019-03-08 2019-05-07 哈尔滨工程大学 一种gnss/ins松组合时延误差的估计及补偿方法
CN109764892A (zh) * 2018-12-30 2019-05-17 华中光电技术研究所(中国船舶重工集团有限公司第七一七研究所) 近红外星等的校准方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140058696A1 (en) * 2012-08-24 2014-02-27 Samsung Electro-Mechanics Co., Ltd. Apparatus and method for controlling gain automatically in inertia sensor
CN103487820A (zh) * 2013-09-30 2014-01-01 东南大学 一种车载捷联/卫星紧组合无缝导航方法
CN108205432A (zh) * 2016-12-16 2018-06-26 中国航天科工飞航技术研究院 一种观测实验数据异常值的实时剔除方法
CN107121684A (zh) * 2017-05-02 2017-09-01 北京航空航天大学 一种基于残差卡方检验法的gps诱骗识别和阈值决策方法
CN109084761A (zh) * 2018-08-07 2018-12-25 合肥正阳光电科技有限责任公司 一种基于惯性寻北装置的抗冲击寻北解算方法
CN109471143A (zh) * 2018-12-11 2019-03-15 北京交通大学 自适应容错的列车组合定位方法
CN109764892A (zh) * 2018-12-30 2019-05-17 华中光电技术研究所(中国船舶重工集团有限公司第七一七研究所) 近红外星等的校准方法
CN109724598A (zh) * 2019-03-08 2019-05-07 哈尔滨工程大学 一种gnss/ins松组合时延误差的估计及补偿方法

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021213525A1 (zh) * 2020-04-24 2021-10-28 北京智行者科技有限公司 自适应的gps错误观测值识别方法、存储介质、产品及系统
US11598885B2 (en) 2020-04-24 2023-03-07 Beijing Idriverplus Information Technology Co., Ltd. Method for adaptive identification of erroneous GPS observed value
CN114485725A (zh) * 2021-12-22 2022-05-13 深圳元戎启行科技有限公司 数据异常检测方法、自动驾驶平台以及计算机可读存储介质
CN114485725B (zh) * 2021-12-22 2024-05-24 深圳元戎启行科技有限公司 数据异常检测方法、自动驾驶平台以及计算机可读存储介质

Also Published As

Publication number Publication date
CN110906929B (zh) 2022-01-25

Similar Documents

Publication Publication Date Title
CN110140065B (zh) Gnss接收机保护等级
CN103348261B (zh) 使用二重差分模糊度约束条件解析整周模糊度的导航系统和方法
US10041808B2 (en) Method of sensor data fusion
US8773303B2 (en) Position tracking device and method
US11442174B2 (en) Systems and methods for inertial measurement unit aided detection and exclusion against spoofing attacks
CN111060133B (zh) 一种用于城市复杂环境的组合导航完好性监测方法
CN110906929B (zh) 一种卫星异常数据识别方法及装置
CN102221365A (zh) 用于确定惯性导航系统故障的系统和方法
EP3012665B1 (en) Systems and methods for averaging satellite sigmas and readmitting excluded satellite measurements into differential corrections and integrity monitors
CN111708054B (zh) 基于粒子群优化算法的araim垂直保护级优化方法
EP2685214B1 (en) Multiple truth reference system and method
CN112130177A (zh) 一种基于稳定分布的地基增强系统完好性监测方法
Gao et al. Double-channel sequential probability ratio test for failure detection in multisensor integrated systems
CN114152958B (zh) 一种基于多数据源的机载卫星导航欺骗式干扰检测方法
CN113267792A (zh) 一种基于监测站网的电离层异常监测预警方法
JIANG et al. A robust fault detection algorithm for the GNSS/INS integrated navigation systems
CN105511481B (zh) 一种星载定轨优化方法
CN111856526B (zh) 非直达径卫星导航信号识别方法、系统及介质
CN115291253B (zh) 一种基于残差检测的车辆定位完好性监测方法及系统
CN112198533A (zh) 一种多假设下的地基增强系统完好性评估系统及方法
CN116756686A (zh) 一种飞行器的强抗扰高度状态估计方法及系统
CN111999750B (zh) 针对杆臂不准的实时单站周跳探测改进方法
CN107783166B (zh) Gps速度异常检测与修复的方法及其系统
CN112764059B (zh) 接收机自主完好性监测方法及装置
Grof et al. Voting based navigation system fault detection for aircraft position information during final approach

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
CP03 Change of name, title or address
CP03 Change of name, title or address

Address after: 511400 room 1702, No. 1, Cuiying street, Nansha District, Guangzhou, Guangdong Province

Patentee after: Guangzhou Zhongke Aerospace Exploration Technology Co.,Ltd.

Address before: 100176 1201, 12 / F, building 6, yard 18, Kechuang 13th Street, Beijing Economic and Technological Development Zone, Daxing District, Beijing

Patentee before: Beijing Zhongke Aerospace Exploration Technology Co.,Ltd.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20230411

Address after: No. 18 Zhixin Fourth Road, Nansha District, Guangzhou City, Guangdong Province, 511400

Patentee after: Zhongke Aerospace (Guangzhou) Equipment Industry Co.,Ltd.

Address before: 511400 room 1702, No. 1, Cuiying street, Nansha District, Guangzhou, Guangdong Province

Patentee before: Guangzhou Zhongke Aerospace Exploration Technology Co.,Ltd.