CN111308504A - 一种基于相位观测值实时探测北斗卫星轨道机动的方法 - Google Patents

一种基于相位观测值实时探测北斗卫星轨道机动的方法 Download PDF

Info

Publication number
CN111308504A
CN111308504A CN201911173999.9A CN201911173999A CN111308504A CN 111308504 A CN111308504 A CN 111308504A CN 201911173999 A CN201911173999 A CN 201911173999A CN 111308504 A CN111308504 A CN 111308504A
Authority
CN
China
Prior art keywords
observation
beidou satellite
value
satellite
maneuver
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
CN201911173999.9A
Other languages
English (en)
Other versions
CN111308504B (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.)
National Time Service Center of CAS
Original Assignee
National Time Service Center of CAS
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 National Time Service Center of CAS filed Critical National Time Service Center of CAS
Priority to CN201911173999.9A priority Critical patent/CN111308504B/zh
Publication of CN111308504A publication Critical patent/CN111308504A/zh
Application granted granted Critical
Publication of CN111308504B publication Critical patent/CN111308504B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/03Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
    • G01S19/08Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing integrity information, e.g. health of satellites or quality of ephemeris data

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

一种基于相位观测值实时探测北斗卫星轨道机动的方法
技术领域
本发明涉及星轨道机动实时探测技术领域,具体涉及一种基于相位观测值实时探测北斗卫星轨道机动的方法。
背景技术
北斗卫星导航星座由地球静止轨道、倾斜地球轨道和中圆轨道三种星座组成。在日月等星体的引力作用下,卫星的轨道会发生变化。因此,为了保证卫星在设计的轨道上运行,需要通过轨道机动方式进行轨道调整。尤其是地球静止轨道卫星,轨道机动更为频繁。当前,轨道机动的探测方法存在精度低,精准性差,不易实时操作或工作方式复杂等诸多不足。
如何提高轨道机动探测的准确性,同时方便实时操作,简单方面,能高效探测出轨道机动,减小卫星机动对精密定轨的影响,并延长用户端卫星的可用观测时间,对于提升精密导航定位授时服务具有重要价值。
发明内容
针对上述存在的技术问题,本发明提供了一种基于观测值历元差分的原理,采用相位观测值历元差分测速模型求解测站的观测值矩阵、三维速度和观测值残差,并确定观测值矩阵和观测值残差STD的经验阈值,通过观测值残差的STD来判断卫星是否发生机动,通过观测值矩阵来识别机动卫星的PRN号码的实时探测北斗卫星轨道机动的方法。
本发明的技术方案为:一种基于相位观测值实时探测北斗卫星轨道机动的方法,包括以下步骤:
步骤一、数据获取;
收集测站上北斗卫星的双频观测值以及北斗数据处理需要的辅助数据,辅助数据包括北斗卫星广播星历、天线相位中心和地球自转参数;
步骤二、数据预处理;
根据步骤一获取的北斗卫星广播星历,对北斗卫星的双频观测值进行数据质量检查和粗差剔除,删除无北斗卫星星历和观测值不完整的数据,得到干净的数据;
步骤三、误差修正;
对步骤二得到的干净的数据进行相对论、潮汐、天线相位中心、对流层和地球自转误差的修正;
步骤四、建立方程;
首先对步骤二得到的北斗卫星干净的双频相位观测值分别进行无电离层组合,形成无电离层组合观测值;同时根据北斗卫星广播星历和测站初始位置,建立北斗卫星历元间差分的观测方程,如式(1):
Figure BDA0002289491320000021
根据北斗卫星位置和测站概略位置计算卫星的高度角,并根据北斗卫星高度角和观测噪声确定其对应的随机模型如式(2):
Figure BDA0002289491320000022
其中,ΦIF为无电离层组合的相位观测值,角标s代表卫星,角标r代表测站,e为卫星和接收机天线间的单位矢量,Δ代表历元间差分,t代表观测时间,c代表光速,δTr为接收机钟差,
Figure BDA0002289491320000023
为其他非模型化误差,包括星历残差、大气残差和多路径残差;θ为卫星高度角;
观测方程线性化后,其对应的误差方程和权阵如式(3);
Figure BDA0002289491320000024
其中,V为观测值残差矢量,A为未知参数的系数矩阵,
Figure BDA0002289491320000025
为未知参数矢量,L为观测值扣除误差后的常数项矢量;
步骤五、参数解算;
根据步骤四的误差方程(3)和随机模型(2)进行最小二乘参数估计,获取三维速度值
Figure BDA0002289491320000031
以及对应的观测值残差(V)和观测值常数值矩阵(L);
步骤六、计算经验阈值;
根据步骤五获取的观测值残差V,计算残差V的STD序列,并将无轨道机动条件下一天的STD序列平均值的3倍作为判断北斗卫星机动是否发生的经验阈值;
步骤七、判断轨道机动;
若当前历元计算出的观测值残差STD超出步骤六得到的经验阈值,并且持续5分钟一直超出阈值,则判定为北斗卫星轨道机动发生;若当前历元计算出的观测值残差STD未超出步骤六得到的经验阈值,则判定为北斗卫星轨道机动未发生;
步骤八、判断机动卫星号;
若探测出北斗卫星轨道机动发生,则根据步骤五得到的观测值常数矩阵L,判定最大的L值对应的北斗卫星为即轨道机动卫星。
进一步地,步骤一中北斗卫星的双频观测值包括双频伪距观测值和双频载波相位观测值,能够有效降低外界环境因素对北斗卫星的双频相位观测值的影响,提高北斗卫星的双频相位观测值的数据质量。
进一步地,步骤二的具体操作为:通过双频伪距观测值定位计算基准站的坐标,与测站坐标进行对比,若对比结果超过预设阈值,则通过残差分析算法剔除删除无北斗卫星星历和观测值不完整的数据,得到干净的数据,同时剔除该粗差北斗卫星。
进一步地,步骤3中相对论和潮汐修正使用IERS Conventions 2010中指定的模型修正,天线相位中心修正采用igs14.atx模型修正,对流层修正采用Saastamoinen模型修正,地球自转误差修正使用IERS EOP C04模型修正。
进一步地,步骤三还包括对天线相位缠绕、固体潮、多径误差和观测噪声的修正,天线相位缠绕使用Wu模型进行改正,固体潮采用IERS2010模型进行修正,多径误差,观测噪声在数据处理中忽略不计。
进一步地,步骤四中形成无电离层组合观测值后,对无电离层组合观测值应用滤波器过滤,以获得表示对无电离层组合的载波-相位整数模糊度估计的状态值和关联统计信息。
与现有技术相比,本发明的有益效果:
第一,直接相位观测值历元间差分处理,得到高精度的速度和对应参数信息,可以实现精准探测。本发明将北斗相位观测值直接进行历元间差分处理,既消除了模糊度参数,并可直接获取速度信息和对应的参数求解信息。因采用高精度的相位观测值,因此精准性高。
第二,有效减弱共性误差,提高参数求解的性能。历元间差分原理可以消除或减弱一些诸如多路径、坐标系统、大气,星历等共性误差,提高了参数求解的误差修正水平和参数估计性能。
第三,方法简单可靠,方便实时实施。相对于其他轨道机动探测方法一方面需要辅助信息,另一方面需要多站联合解算,本发明在用户端即可以简单实施,而且只需简单的历元差分解算就可以完成轨道机动探测,方便实时应用。
附图说明
图1是本发明的流程框图;
具体实施方式
实施例1:一种基于相位观测值实时探测北斗卫星轨道机动的方法,包括以下步骤:
步骤一、数据获取;
收集测站上北斗卫星的双频观测值以及北斗数据处理需要的辅助数据,辅助数据包括北斗卫星广播星历、天线相位中心和地球自转参数;
步骤二、数据预处理;
根据步骤一获取的北斗卫星广播星历,对北斗卫星的双频观测值进行数据质量检查和粗差剔除,删除无北斗卫星星历和观测值不完整的数据,得到干净的数据;
步骤三、误差修正;
对步骤二得到的干净的数据进行相对论、潮汐、天线相位中心、对流层和地球自转误差的修正;
步骤四、建立方程;
首先对步骤二得到的北斗卫星干净的双频相位观测值分别进行无电离层组合,形成无电离层组合观测值;同时根据北斗卫星广播星历和测站初始位置,建立北斗卫星历元间差分的观测方程,如式(1):
Figure BDA0002289491320000051
根据北斗卫星位置和测站概略位置计算卫星的高度角,并根据北斗卫星高度角和观测噪声确定其对应的随机模型如式(2):
Figure BDA0002289491320000052
其中,ΦIF为无电离层组合的相位观测值,角标s代表卫星,角标r代表测站,e为卫星和接收机天线间的单位矢量,Δ代表历元间差分,t代表观测时间,c代表光速,δTr为接收机钟差,
Figure BDA0002289491320000053
为其他非模型化误差,包括星历残差、大气残差和多路径残差;θ为卫星高度角;
观测方程线性化后,其对应的误差方程和权阵如式(3);
Figure BDA0002289491320000054
其中,V为观测值残差矢量,A为未知参数的系数矩阵,
Figure BDA0002289491320000061
为未知参数矢量,L为观测值扣除误差后的常数项矢量;
步骤五、参数解算;
根据步骤四的误差方程(3)和随机模型(2)进行最小二乘参数估计,获取三维速度值
Figure BDA0002289491320000062
以及对应的观测值残差(V)和观测值常数值矩阵(L);
步骤六、计算经验阈值;
根据步骤五获取的观测值残差V,计算残差V的STD序列,并将无轨道机动条件下一天的STD序列平均值的3倍作为判断北斗卫星机动是否发生的经验阈值;
步骤七、判断轨道机动;
若当前历元计算出的观测值残差STD超出步骤六得到的经验阈值,并且持续5分钟一直超出阈值,则判定为北斗卫星轨道机动发生;若当前历元计算出的观测值残差STD未超出步骤六得到的经验阈值,则判定为北斗卫星轨道机动未发生;
步骤八、判断机动卫星号;
若探测出北斗卫星轨道机动发生,则根据步骤五得到的观测值常数矩阵L,判定最大的L值对应的北斗卫星为即轨道机动卫星。
实施例2:一种基于相位观测值实时探测北斗卫星轨道机动的方法,包括以下步骤:
步骤一、数据获取;
收集测站上北斗卫星的双频观测值以及北斗数据处理需要的辅助数据,辅助数据包括北斗卫星广播星历、天线相位中心和地球自转参数;北斗卫星的双频观测值包括双频伪距观测值和双频载波相位观测值,能够有效降低外界环境因素对北斗卫星的双频相位观测值的影响,提高北斗卫星的双频相位观测值的数据质量;
步骤二、数据预处理;
根据步骤一获取的北斗卫星广播星历,通过双频伪距观测值定位计算基准站的坐标,与测站坐标进行对比,若对比结果超过预设阈值,则通过残差分析算法剔除删除无北斗卫星星历和观测值不完整的数据,得到干净的数据,同时剔除该粗差北斗卫星;
步骤三、误差修正;
对步骤二得到的干净的数据进行相对论、潮汐、天线相位中心、对流层和地球自转误差的修正;相对论和潮汐修正使用IERS Conventions 2010中指定的模型修正,天线相位中心修正采用igs14.atx模型修正,对流层修正采用Saastamoinen模型修正,地球自转误差修正使用IERS EOP C04模型修正;还包括对天线相位缠绕、固体潮、多径误差和观测噪声的修正,天线相位缠绕使用Wu模型进行改正,固体潮采用IERS2010模型进行修正,多径误差,观测噪声在数据处理中忽略不计;
步骤四、建立方程;
首先对步骤二得到的北斗卫星干净的双频相位观测值分别进行无电离层组合,形成无电离层组合观测值;同时根据北斗卫星广播星历和测站初始位置,建立北斗卫星历元间差分的观测方程,如式(1):
Figure BDA0002289491320000071
根据北斗卫星位置和测站概略位置计算卫星的高度角,并根据北斗卫星高度角和观测噪声确定其对应的随机模型如式(2):
Figure BDA0002289491320000072
其中,ΦIF为无电离层组合的相位观测值,角标s代表卫星,角标r代表测站,e为卫星和接收机天线间的单位矢量,Δ代表历元间差分,t代表观测时间,c代表光速,δTr为接收机钟差,
Figure BDA0002289491320000081
为其他非模型化误差,包括星历残差、大气残差和多路径残差;θ为卫星高度角;
观测方程线性化后,其对应的误差方程和权阵如式(3);
Figure BDA0002289491320000082
其中,V为观测值残差矢量,A为未知参数的系数矩阵,
Figure BDA0002289491320000083
为未知参数矢量,L为观测值扣除误差后的常数项矢量;形成无电离层组合观测值后,对无电离层组合观测值应用滤波器过滤,以获得表示对无电离层组合的载波-相位整数模糊度估计的状态值和关联统计信息;
步骤五、参数解算;
根据步骤四的误差方程(3)和随机模型(2)进行最小二乘参数估计,获取三维速度值
Figure BDA0002289491320000084
以及对应的观测值残差(V)和观测值常数值矩阵(L);
步骤六、计算经验阈值;
根据步骤五获取的观测值残差V,计算残差V的STD序列,并将无轨道机动条件下一天的STD序列平均值的3倍作为判断北斗卫星机动是否发生的经验阈值;
步骤七、判断轨道机动;
若当前历元计算出的观测值残差STD超出步骤六得到的经验阈值,并且持续5分钟一直超出阈值,则判定为北斗卫星轨道机动发生;若当前历元计算出的观测值残差STD未超出步骤六得到的经验阈值,则判定为北斗卫星轨道机动未发生;
步骤八、判断机动卫星号;
若探测出北斗卫星轨道机动发生,则根据步骤五得到的观测值常数矩阵L,判定最大的L值对应的北斗卫星为即轨道机动卫星。

Claims (7)

1.一种基于相位观测值实时探测北斗卫星轨道机动的方法,其特征在于,包括以下步骤:
步骤一、数据获取;
收集测站上北斗卫星的双频观测值以及北斗数据处理需要的辅助数据,所述辅助数据包括北斗卫星广播星历、天线相位中心和地球自转参数;
步骤二、数据预处理;
根据步骤一获取的北斗卫星广播星历,对所述北斗卫星的双频观测值进行数据质量检查和粗差剔除,删除无北斗卫星星历和观测值不完整的数据,得到干净的数据;
步骤三、误差修正;
对步骤二得到的干净的数据进行相对论、潮汐、天线相位中心、对流层和地球自转误差的修正;
步骤四、建立方程;
首先对步骤二预处理后得到的北斗卫星干净的双频相位观测值分别进行无电离层组合,形成无电离层组合观测值;同时根据所述北斗卫星广播星历和测站初始位置,建立北斗卫星历元间差分的观测方程,如式(1):
Figure FDA0002289491310000011
根据北斗卫星位置和测站概略位置计算卫星的高度角,并根据北斗卫星高度角和观测噪声确定其对应的随机模型如式(2):
Figure FDA0002289491310000012
其中,ΦIF为无电离层组合的相位观测值,角标s代表卫星,角标r代表测站,e为卫星和接收机天线间的单位矢量,Δ代表历元间差分,t代表观测时间,c代表光速,δTr为接收机钟差,
Figure FDA0002289491310000021
为其他非模型化误差,包括星历残差、大气残差和多路径残差;θ为卫星高度角;
观测方程线性化后,其对应的误差方程和权阵如式(3);
Figure FDA0002289491310000022
其中,V为观测值残差矢量,A为未知参数的系数矩阵,
Figure FDA0002289491310000023
为未知参数矢量,L为观测值扣除误差后的常数项矢量;
步骤五、参数解算;
根据步骤四所述的误差方程(3)和随机模型(2)进行最小二乘参数估计,获取三维速度值
Figure FDA0002289491310000024
以及对应的观测值残差(V)和观测值常数值矩阵(L);
步骤六、计算经验阈值;
根据步骤五获取的观测值残差V,计算残差V的STD序列,并将无轨道机动条件下一天的STD序列平均值的3倍作为判断北斗卫星机动是否发生的经验阈值;
步骤七、判断轨道机动;
若当前历元计算出的观测值残差STD超出步骤六得到的经验阈值,并且持续5分钟一直超出阈值,则判定为北斗卫星轨道机动发生;若当前历元计算出的观测值残差STD未超出步骤六得到的经验阈值,则判定为北斗卫星轨道机动未发生;
步骤八、判断机动卫星号;
若探测出北斗卫星轨道机动发生,则根据步骤五得到的观测值常数矩阵L,判定最大的L值对应的北斗卫星为即轨道机动卫星。
2.根据权利要求1所述的一种基于相位观测值实时探测北斗卫星轨道机动的方法,其特征在于,所述步骤一中北斗卫星的双频观测值包括双频伪距观测值和双频载波相位观测值。
3.根据权利要求2所述的一种基于相位观测值实时探测北斗卫星轨道机动的方法,其特征在于,所述步骤二的具体操作为:通过双频伪距观测值定位计算基准站的坐标,与测站坐标进行对比,若对比结果超过预设阈值,则通过残差分析算法剔除删除无北斗卫星星历和观测值不完整的数据,得到干净的数据,同时剔除该粗差北斗卫星。
4.根据权利要求1所述的一种基于相位观测值实时探测北斗卫星轨道机动的方法,其特征在于,所述步骤3中相对论和潮汐修正使用IERS Conventions2010中指定的模型修正,天线相位中心修正采用igs14.atx模型修正,对流层修正采用Saastamoinen模型修正,地球自转误差修正使用IERS EOP C04模型修正。
5.根据权利要求1所述的一种基于相位观测值实时探测北斗卫星轨道机动的方法,其特征在于,所述步骤三还包括对天线相位缠绕、固体潮、多径误差和观测噪声的修正,所述天线相位缠绕使用Wu模型进行改正,所述固体潮采用IERS2010模型进行修正,所述多径误差,观测噪声在数据处理中忽略不计。
6.根据权利要求1所述的一种基于相位观测值实时探测北斗卫星轨道机动的方法,其特征在于,所述步骤四中形成无电离层组合观测值后,对所述无电离层组合观测值应用滤波器过滤。
7.根据权利要求1所述的一种基于相位观测值实时探测北斗卫星轨道机动的方法,其特征在于,所述步骤一中北斗卫星的双频观测值包括双频载波相位观测值和双频伪距观测值。
CN201911173999.9A 2019-11-26 2019-11-26 一种基于相位观测值实时探测北斗卫星轨道机动的方法 Active CN111308504B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911173999.9A CN111308504B (zh) 2019-11-26 2019-11-26 一种基于相位观测值实时探测北斗卫星轨道机动的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911173999.9A CN111308504B (zh) 2019-11-26 2019-11-26 一种基于相位观测值实时探测北斗卫星轨道机动的方法

Publications (2)

Publication Number Publication Date
CN111308504A true CN111308504A (zh) 2020-06-19
CN111308504B CN111308504B (zh) 2023-10-31

Family

ID=71154275

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911173999.9A Active CN111308504B (zh) 2019-11-26 2019-11-26 一种基于相位观测值实时探测北斗卫星轨道机动的方法

Country Status (1)

Country Link
CN (1) CN111308504B (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112198532A (zh) * 2020-09-30 2021-01-08 哈尔滨工程大学 一种顾及参数灵敏度的北斗广播星历机动型故障监测方法
CN113406678A (zh) * 2021-05-12 2021-09-17 中铁第四勘察设计院集团有限公司 滤波方法、装置、设备和存储介质
CN113568011A (zh) * 2021-07-16 2021-10-29 航天科工海鹰集团有限公司 一种ppk多径消除技术
CN114254568A (zh) * 2022-02-28 2022-03-29 浙江国遥地理信息技术有限公司 基于人工智能决策树的gps遥感洪水预警方法
CN114355418A (zh) * 2021-12-07 2022-04-15 中国科学院国家授时中心 基于北斗地基增强系统的事后数据质量评估方法和系统
CN114690210A (zh) * 2022-04-18 2022-07-01 山东大学 一种基于多普勒观测值的北斗卫星机动探测方法
CN114895328A (zh) * 2022-05-12 2022-08-12 中国科学院国家授时中心 基于多普勒观测值的北斗卫星轨道机动识别方法及系统
CN115096315A (zh) * 2022-06-07 2022-09-23 哈尔滨工业大学 针对稀疏数据的航天器目标机动检测方法
CN115856982A (zh) * 2023-02-22 2023-03-28 广州导远电子科技有限公司 一种相对位置获取方法、装置、存储介质及电子设备
CN115902981A (zh) * 2022-11-14 2023-04-04 中南大学 列车定位优化方法、系统及轨道交通车辆

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106501815A (zh) * 2016-09-26 2017-03-15 中国人民解放军63920部队 一种仅天基测角跟踪的空间目标轨道机动融合检测方法
CN108535746A (zh) * 2018-02-27 2018-09-14 中国科学院测量与地球物理研究所 一种探测gnss卫星轨道机动的方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106501815A (zh) * 2016-09-26 2017-03-15 中国人民解放军63920部队 一种仅天基测角跟踪的空间目标轨道机动融合检测方法
CN108535746A (zh) * 2018-02-27 2018-09-14 中国科学院测量与地球物理研究所 一种探测gnss卫星轨道机动的方法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
R.TU ET AL.: "Real time monitoring ground motion using GPS with real time corrections", 《SURVEY REVIEW》 *
R.TU ET AL.: "Real time monitoring ground motion using GPS with real time corrections", 《SURVEY REVIEW》, 31 December 2014 (2014-12-31), pages 2 *
戴小蕾: "基于平方根信息滤波的GNSS导航卫星实时精密定轨理论与方法", 《中国博士学位论文全文数据库 基础科学辑》 *
戴小蕾: "基于平方根信息滤波的GNSS导航卫星实时精密定轨理论与方法", 《中国博士学位论文全文数据库 基础科学辑》, 15 June 2017 (2017-06-15), pages 2 *
李浩军 等: "基于历元间差分技术的精密单点定位研究", 《宇航学报》 *
李浩军 等: "基于历元间差分技术的精密单点定位研究", 《宇航学报》, 31 March 2010 (2010-03-31), pages 748 - 749 *
魏振军: "《概率论与数理统计》", 30 April 2009, pages: 170 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112198532B (zh) * 2020-09-30 2023-07-21 哈尔滨工程大学 一种顾及参数灵敏度的北斗广播星历机动型故障监测方法
CN112198532A (zh) * 2020-09-30 2021-01-08 哈尔滨工程大学 一种顾及参数灵敏度的北斗广播星历机动型故障监测方法
CN113406678A (zh) * 2021-05-12 2021-09-17 中铁第四勘察设计院集团有限公司 滤波方法、装置、设备和存储介质
CN113406678B (zh) * 2021-05-12 2022-06-10 中铁第四勘察设计院集团有限公司 滤波方法、装置、设备和存储介质
CN113568011A (zh) * 2021-07-16 2021-10-29 航天科工海鹰集团有限公司 一种ppk多径消除技术
CN114355418A (zh) * 2021-12-07 2022-04-15 中国科学院国家授时中心 基于北斗地基增强系统的事后数据质量评估方法和系统
CN114254568A (zh) * 2022-02-28 2022-03-29 浙江国遥地理信息技术有限公司 基于人工智能决策树的gps遥感洪水预警方法
CN114690210A (zh) * 2022-04-18 2022-07-01 山东大学 一种基于多普勒观测值的北斗卫星机动探测方法
CN114895328A (zh) * 2022-05-12 2022-08-12 中国科学院国家授时中心 基于多普勒观测值的北斗卫星轨道机动识别方法及系统
CN114895328B (zh) * 2022-05-12 2024-05-14 中国科学院国家授时中心 基于多普勒观测值的北斗卫星轨道机动识别方法及系统
CN115096315A (zh) * 2022-06-07 2022-09-23 哈尔滨工业大学 针对稀疏数据的航天器目标机动检测方法
CN115902981A (zh) * 2022-11-14 2023-04-04 中南大学 列车定位优化方法、系统及轨道交通车辆
CN115902981B (zh) * 2022-11-14 2024-01-30 中南大学 列车定位优化方法、系统及轨道交通车辆
CN115856982A (zh) * 2023-02-22 2023-03-28 广州导远电子科技有限公司 一种相对位置获取方法、装置、存储介质及电子设备

Also Published As

Publication number Publication date
CN111308504B (zh) 2023-10-31

Similar Documents

Publication Publication Date Title
CN111308504B (zh) 一种基于相位观测值实时探测北斗卫星轨道机动的方法
CN107102346B (zh) 一种基于北斗系统的多天线测姿方法
CN104714244B (zh) 一种基于抗差自适应Kalman滤波的多系统动态PPP解算方法
CN109581452B (zh) 一种gnss参考站载波相位整周模糊度解算方法
CN110045407B (zh) 一种分布式伪卫星/gnss优化定位方法
CN108120994B (zh) 一种基于星载gnss的geo卫星实时定轨方法
CN110208836B (zh) 基于卡尔曼滤波的gnss高适应性周跳探测与修复方法
CN110007326B (zh) 一种用于星基增强系统的双频测距误差参数生成方法
CN111983641B (zh) 一种用于实时生成北斗星基增强系统完好性参数的方法
CN112305574A (zh) 北斗gnss卫星实时定位定向数据预处理系统及方法
Ke et al. Integrating GPS and LEO to accelerate convergence time of precise point positioning
CN104865587A (zh) 基于geo星座粗定位的北斗接收机快速定位方法
CN110988917B (zh) 一种卫星轨道机动状态实时监测方法
CN112285745A (zh) 基于北斗三号卫星导航系统的三频模糊度固定方法及系统
CN111290004A (zh) 伪距差分定位方法、装置、电子设备及存储介质
CN113253314A (zh) 一种低轨卫星间时间同步方法及系统
CN111505693A (zh) 一种载波相位历元差分和多普勒观测融合测速方法
CN114935770B (zh) 一种多历元加快精密单点定位收敛速度的方法及装置
CN111505677A (zh) 一种基于地面参考站观测的geo卫星轨道机动修复方法
CN103675858A (zh) 北斗系统b1与gps系统l1载波相位混频差分方法
CN111562594A (zh) 基于站间单差和模糊度固定的北斗精密时间传递方法
CN110058274B (zh) 一种卫星导航系统间的时差监测方法及系统
CN104991265B (zh) 一种北斗卫星导航系统用户统一性定位方法
CN114895328B (zh) 基于多普勒观测值的北斗卫星轨道机动识别方法及系统
CN111650612A (zh) 一种基于双差模式的实时精密时间传递方法

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