CN114675311B - 高性能pvt解算方法及导航方法 - Google Patents

高性能pvt解算方法及导航方法 Download PDF

Info

Publication number
CN114675311B
CN114675311B CN202210597377.4A CN202210597377A CN114675311B CN 114675311 B CN114675311 B CN 114675311B CN 202210597377 A CN202210597377 A CN 202210597377A CN 114675311 B CN114675311 B CN 114675311B
Authority
CN
China
Prior art keywords
satellite
state information
time
calculation
current
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
Application number
CN202210597377.4A
Other languages
English (en)
Other versions
CN114675311A (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.)
Changsha Jinwei Information Technology Co ltd
Original Assignee
Changsha Jinwei Information 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 Changsha Jinwei Information Technology Co ltd filed Critical Changsha Jinwei Information Technology Co ltd
Priority to CN202210597377.4A priority Critical patent/CN114675311B/zh
Publication of CN114675311A publication Critical patent/CN114675311A/zh
Application granted granted Critical
Publication of CN114675311B publication Critical patent/CN114675311B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/35—Constructional details or hardware or software details of the signal processing chain
    • G01S19/37—Hardware or software details of the signal processing chain
    • 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/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining 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
    • 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/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining 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/42—Determining position
    • 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/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining 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/52—Determining velocity

Abstract

本发明公开了一种高性能PVT解算方法，包括获取PVT解算所需要的输入数据；采用第一线程，以设定的调度频度、卫星号周期性方式和时间间隔抽样方式进行卫星状态信息计算；同时采用第二线程，以二次曲线拟合的方式进行卫星状态信息外推和最小二乘解算；根据最小二乘解算结果得到最终的定位结果并完成PVT解算。本发明还公开了一种包括所述高性能PVT解算方法的导航方法。本发明的两个线程独立运算，两者之间进行卫星状态信息的交互，降低了系统负担，减少了解算耗时，提高了解算效率，而且定位精度较高、可靠性较好。

Description

高性能PVT解算方法及导航方法
技术领域
本发明属于导航技术领域，具体涉及一种高性能PVT解算方法及导航方法。
背景技术
随着经济技术的发展和人们生活水平的提高，导航技术已经广泛应用于人们的生产和生活当中，给人们的生产和生活带来了无尽的便利。因此，保证导航过程的精度和实时性，就成为了研究人员的研究重点。
GNSS导航接收机解算的主要目的是为了求解出接收机的位置（Position）、速度（Velocity）和时间（Time），这个过程称作PVT解算。为了满足多系统、多频点的应用场景，GNSS导航接收机需要具备多系统、多频点同时解算的特性。参与PVT解算的有效卫星数目和信号分量越多，就越能提高接收机定位的可靠性和定位精度；但是，相应也会增大PVT解算的运算量和耗时，增加系统负担；这将使得GNSS导航接收机的应用很难得到本质的提升和突破。
目前，在传统的接收机的解算方法中，PVT解算过程通常包括以下几个步骤：（1）获取观测量数据和导航电文数据；（2）计算卫星状态信息（包括卫星位置、速度和加速度等）；（3）最小二乘解算；（4）定位结果处理及输出。该过程的典型方法流程图如图1所示。
在上述传统解算的过程中，四个步骤顺序执行，每一个步骤的计算结果均要作为下一个步骤的输入。整个PVT求解运算过程的耗时主要集中在计算卫星状态信息和最小二乘解算两个部分；而其中计算卫星状态信息耗时更多，占整个PVT解算耗时的60~70%。
在多系统、多卫星、多频点的应用场景下，为了提高定位精度，同一颗卫星的不同频点往往会使用当前频点的导航电文数据参与解算，这就需要基于当前频点的导航电文数据再次计算该卫星的更高精度的卫星状态信息；这使得运算量成倍增加。另外，高频度解算时，高频的每一拍PVT解算都需要完成所有有效卫星的状态信息计算，这使得高频度解算的接收机压力更大；此外，现有技术在进行PVT解算的过程中，其计算流程繁复，效率也较低。因此，传统的PVT解算过程，解算耗时巨大，限制了高频度、低成本等应用需求。
发明内容
本发明的目的之一在于提供一种精度较高、可靠性较好且解算耗时更低的高性能PVT解算方法。
本发明的目的之二在于提供一种包括了所述高性能PVT解算方法的导航方法。
本发明提供的这种高性能PVT解算方法，包括如下步骤：
S1. 获取PVT解算所需要的输入数据；
S2. 根据步骤S1获取的输入数据，采用第一线程，以设定的调度频度、卫星号周期性方式和时间间隔抽样方式，将所有卫星状态信息的计算工作分摊到不同的时隙，并逐个卫星进行卫星状态信息的计算；
S3. 在步骤S2进行的同时，采用第二线程，结合二次运动模型理论，以二次曲线拟合的方式进行卫星状态信息外推，并进行最小二乘解算；
S4. 根据步骤S3得到的最小二乘解算结果，得到最终的定位结果，完成PVT解算。
采用第一线程进行计算的优先级低于采用第二线程进行计算的优先级。
步骤S2计算得到卫星状态信息的计算结果后，将计算结果更新到全局卫星状态信息缓存中，供步骤S3进行卫星状态信息外推使用；在对全局卫星状态信息缓存进行任何操作时，均采用互斥量对全局卫星状态信息缓存进行防护，从而保证卫星状态信息的完整性、操作的可靠性和安全性。
步骤S1所述的获取PVT解算所需要的输入数据，具体包括obs观测量数据和nav导航电文数据；其中，nav导航电文数据包括eph星历数据，satclk卫星钟差数据，iono电离层数据和tgd群延迟数据；然后将获取的输入数据按照设定的数据结构和格式进行组装。
步骤S2所述的设定的调度频度，具体为100Hz。
步骤S2所述的卫星号周期性方式，具体包括如下步骤：
总卫星数N为
Figure 92953DEST_PATH_IMAGE002
，其中N 1为GPS卫星的数量，N 2为GLO卫星的数量，N 3为GAL卫星的数量，N 4为北斗卫星的数量；在设定的调度频度下，每一拍仅计算一颗卫星的状态信息；当计算完最后一颗卫星的状态信息后，再从第一颗卫星开始周期性的进行计算。
步骤S2所述的时间间隔抽样方式，具体包括如下步骤：
已经计算过第i颗卫星的状态信息后，当再次轮到对第i颗卫星进行计算时，只有当前时刻与上次进行状态信息计算的时刻的时间间隔大于设定值时，才对第i颗卫星进行状态信息的计算；否则，直接跳过第i颗卫星的状态信息计算过程。
步骤S2所述的进行卫星状态信息的计算，具体包括如下步骤：
A. 获取卫星的星历数据信息；
B. 对全局卫星状态信息缓存中的当前卫星的状态信息有效标志进行判断：
若当前卫星的状态信息有效标志无效，则进行后续步骤；
若当前卫星的状态信息有效标志有效，则按照时间间隔抽样方式进行判断；
C. 对当前卫星的星历信息进行检测：
若当前卫星的星历信息有效，则进行后续步骤；
若当前卫星的星历信息无效，则设置当前卫星的状态信息有效标志为无效，将结果保存到全局卫星状态信息缓存中，并等待下一拍的调度计算；
D. 采用当前卫星的星历信息，计算当前时刻的卫星位置；
E. 采用当前卫星的星历信息，计算当前时刻的卫星速度；
F. 采用当前卫星的星历信息，计算当前时刻的卫星加速度；
G. 计算完成后，设置当前卫星的状态信息有效标志为有效，并将计算结果进行保存，保存到全局卫星状态信息缓存中。
步骤A所述的卫星星历的数据信息，具体包括星历参考时间t oe 、卫星轨道长半轴a s 的平方根值
Figure 603569DEST_PATH_IMAGE004
、轨道偏心率e s 、星历参考时间t oe 时刻的轨道倾角i 0、周内时为0时的轨道升交点赤经
Figure 117727DEST_PATH_IMAGE006
、轨道近地角距
Figure 981777DEST_PATH_IMAGE008
、星历参考时间t oe 时刻的平近点角M 0、平均运动角速度校正值
Figure 655204DEST_PATH_IMAGE010
、轨道倾角对时间的变化率
Figure 946508DEST_PATH_IMAGE012
、轨道升交点赤经对时间的变化率
Figure 807017DEST_PATH_IMAGE014
、升交点角距余弦调和校正振幅C uc 、升交点角距正弦调和校正振幅C us 、轨道半径余弦调和校正振幅C rc 、轨道半径正弦调和校正振幅C rs 、轨道倾角余弦调和校正振幅C ic 和轨道倾角正弦调和校正振幅C is 。
所述的步骤D，具体包括如下步骤：
D1. 计算归一化时间t k 为
Figure 740338DEST_PATH_IMAGE016
，式中t为当前时刻；计算平均角速度n为
Figure 143637DEST_PATH_IMAGE018
，式中
Figure 733406DEST_PATH_IMAGE020
为地心引力常数；计算平近点角M k 为
Figure 222156DEST_PATH_IMAGE022
，式中n为计算得到的平均角速度；
D2. 以
Figure 428010DEST_PATH_IMAGE024
为初始值，采用算式
Figure 810450DEST_PATH_IMAGE026
进行迭代，直至精度满足设定的要求，得到最终的偏近点角E；E j 为第j次计算得到的偏近点角；
D3. 计算真近点角v k 为
Figure 974715DEST_PATH_IMAGE028
；计算升交点角距
Figure 950761DEST_PATH_IMAGE030
为
Figure 960305DEST_PATH_IMAGE032
；
D4. 采用如下算式计算摄动校正项：
Figure 462831DEST_PATH_IMAGE034
Figure 266839DEST_PATH_IMAGE036
Figure 730181DEST_PATH_IMAGE038
式中
Figure 668050DEST_PATH_IMAGE040
为升交点角距修正量；
Figure 900448DEST_PATH_IMAGE042
为卫星地心向径修正量；
Figure 406516DEST_PATH_IMAGE044
为卫星轨道倾角修正量；
D5. 计算升交点角距u k 为
Figure 950630DEST_PATH_IMAGE046
；计算卫星矢量长度r k 为
Figure 567556DEST_PATH_IMAGE048
；计算轨道倾角i k 为
Figure 654461DEST_PATH_IMAGE050
；
D6. 采用如下算式计算卫星轨道平面的位置
Figure 331430DEST_PATH_IMAGE052
和升交点赤经
Figure 628419DEST_PATH_IMAGE054
：
Figure 517877DEST_PATH_IMAGE056
Figure 724868DEST_PATH_IMAGE058
Figure 697372DEST_PATH_IMAGE060
式中
Figure 357023DEST_PATH_IMAGE062
为地球自转角速度；
D7. 采用如下算式计算地心地固直角坐标系中的卫星坐标
Figure 315752DEST_PATH_IMAGE064
为
Figure 111670DEST_PATH_IMAGE066
，
Figure 255075DEST_PATH_IMAGE068
和
Figure 136443DEST_PATH_IMAGE070
。
所述的步骤E，具体为将地心地固直角坐标系中的卫星坐标
Figure 898863DEST_PATH_IMAGE064
对时间求导，从而得到当前时刻的卫星速度
Figure 676851DEST_PATH_IMAGE072
为：
Figure 600944DEST_PATH_IMAGE074
Figure 94243DEST_PATH_IMAGE076
Figure 129195DEST_PATH_IMAGE078
式中
Figure 24338DEST_PATH_IMAGE080
为升交点赤经
Figure 119333DEST_PATH_IMAGE082
对时间的变化率；
Figure 240873DEST_PATH_IMAGE084
为轨道倾角i k 对时间的变化率；
Figure 345095DEST_PATH_IMAGE086
为卫星轨道平面的位置
Figure 829166DEST_PATH_IMAGE088
对时间的变化率。
所述的步骤F，具体为将步骤E得到的当前时刻的卫星速度对时间求导，从而得到当前时刻的卫星加速度
Figure 360642DEST_PATH_IMAGE090
为：
Figure 703898DEST_PATH_IMAGE092
Figure 736445DEST_PATH_IMAGE094
Figure 481548DEST_PATH_IMAGE096
式中
Figure 183924DEST_PATH_IMAGE098
为第一系数，且
Figure 14477DEST_PATH_IMAGE100
；
Figure 585136DEST_PATH_IMAGE102
为第二系数，且
Figure 919165DEST_PATH_IMAGE104
；
Figure 792443DEST_PATH_IMAGE106
为轨道倾角i k 对时间的二次求导；
Figure 234926DEST_PATH_IMAGE108
为卫星轨道平面的位置
Figure 219062DEST_PATH_IMAGE110
对时间的二次求导。
所述的步骤S3，具体包括如下步骤：
a. 从全局卫星状态信息缓存中获取当前卫星的状态信息有效标志：
若当前卫星的状态信息有效标志为无效，则本步骤结束；并获取下一颗卫星的状态信息有效标志；
若当前卫星的状态信息有效标志为有效，则计算并得到当前卫星的外推时间间隔
Figure 673177DEST_PATH_IMAGE112
为
Figure 576411DEST_PATH_IMAGE114
，式中t pos 为当前PVT解算时刻，t s 为当前卫星的最近一次状态信息计算的时刻；
b. 设置时间阈值t th1：
若当前卫星的外推时间间隔
Figure 381556DEST_PATH_IMAGE112
大于t th1时，判定当前卫星异常，不再外推当前卫星的状态信息，且不再使用当前卫星参与后续的解算过程；
若当前卫星的外推时间间隔
Figure 169384DEST_PATH_IMAGE112
小于或等于t th1时，采用如下二次曲线模型外推当前卫星在t pos 时刻的位置和速度：
Figure 605569DEST_PATH_IMAGE116
Figure 820650DEST_PATH_IMAGE118
Figure 706566DEST_PATH_IMAGE120
Figure 298084DEST_PATH_IMAGE122
Figure 726792DEST_PATH_IMAGE124
Figure 971828DEST_PATH_IMAGE126
式中上标s表示为当前卫星s；
Figure 833DEST_PATH_IMAGE128
为当前卫星s在t s 时刻的位置；
Figure 130463DEST_PATH_IMAGE130
为当前卫星s在t s 时刻的速度；
Figure 272731DEST_PATH_IMAGE132
为当前卫星s在t s 时刻的加速度；
Figure 95194DEST_PATH_IMAGE134
为当前卫星s外推到t pos 时刻的位置；
Figure 48979DEST_PATH_IMAGE136
为当前卫星s外推到t pos 时刻的速度；
c. 重复步骤a~b，直至遍历完成所有卫星；
d. 将所有正常卫星的状态信息外推后，采用外推得到的状态信息进行后续的最小二乘解算。
本发明还提供了一种包括了所述高性能PVT解算方法的导航方法，包括如下步骤：
（1）实时获取导航所需的数据信息；
（2）采用上述的高性能PVT解算方法进行实时PVT解算，得到实时的定位结果；
（3）根据步骤（2）得到的实时定位结果，完成最终的导航。
本发明提供的这种高性能PVT解算方法及导航方法，将运算量较大的卫星状态信息计算放在优先级较低的第一线程中，采用设定的调度频度、按卫星号周期性的方式、按时间间隔抽样方式执行，将原本需要在一拍内完成的所有卫星状态信息的计算工作，分摊到不同的时隙，并逐个卫星开始计算；把运算量较少的卫星状态信息外推和最小二乘解算放在优先级较高的第二线程中执行，结合二次运动模型理论，采用二次曲线拟合方式进行卫星状态信息外推，并采用系统配置的解算频度进行计算；对传统的PVT解算流程进行优化，将卫星状态信息计算过程进行剥离和优化；两个线程独立运算，两者之间进行卫星状态信息的交互，降低了系统负担，减少了解算耗时，提高了解算效率，而且定位精度较高、可靠性较好。
附图说明
图1为现有的PVT解算方法的方法流程示意图。
图2为本发明的解算方法的方法流程示意图。
图3为本发明的解算方法的时序示意图。
图4为本发明的导航方法的方法流程示意图。
具体实施方式
如图2所示为本发明的解算方法的方法流程示意图，其对应的时序图如图3所示：本发明提供的这种高性能PVT解算方法，包括如下步骤：
S1. 获取PVT解算所需要的输入数据；具体包括obs观测量数据和nav导航电文数据；其中，nav导航电文数据包括eph星历数据，satclk卫星钟差数据，iono电离层数据和tgd群延迟数据等；然后将获取的输入数据按照设定的数据结构和格式进行组装；
S2. 根据步骤S1获取的输入数据，采用第一线程，以设定的调度频度、卫星号周期性方式和时间间隔抽样方式，将所有卫星状态信息的计算工作分摊到不同的时隙，并逐个卫星进行卫星状态信息的计算；
具体实施时，设定的调度频度优选为100Hz，即10ms为一拍，每一拍仅计算一颗卫星的状态信息；
卫星号周期性方式具体为：
总卫星数N为
Figure 310196DEST_PATH_IMAGE002
，其中N 1为GPS卫星的数量，N 2为GLO卫星的数量，N 3为GAL卫星的数量，N 4为北斗卫星的数量；在设定的调度频度下，每一拍仅计算一颗卫星的状态信息；当计算完最后一颗卫星的状态信息后，再从第一颗卫星开始周期性的进行计算；具体实施时，N 1为32，N 2为28，N 3为30，N 4为37，N为127；
时间间隔抽样方式具体为：
已经计算过第i颗卫星的状态信息后，当再次轮到对第i颗卫星进行计算时，只有当前时刻与上次进行状态信息计算的时刻的时间间隔大于设定值（优选为5s）时，才对第i颗卫星进行状态信息的计算；否则，直接跳过第i颗卫星的状态信息计算过程；
根据以上的设定的调度频度、卫星号周期性方式和时间间隔抽样方式，将原本需要在一拍内完成的所有卫星状态信息的计算工作，分摊到不同的时隙，并逐个卫星开始计算；
进行卫星状态信息的计算，则具体包括如下步骤：
A. 获取卫星的星历数据信息；具体包括星历参考时间t oe 、卫星轨道长半轴a s 的平方根值
Figure 309901DEST_PATH_IMAGE004
、轨道偏心率e s 、星历参考时间t oe 时刻的轨道倾角i 0、周内时为0时的轨道升交点赤经
Figure 303265DEST_PATH_IMAGE137
、轨道近地角距
Figure 526436DEST_PATH_IMAGE008
、星历参考时间t oe 时刻的平近点角M 0、平均运动角速度校正值
Figure 122502DEST_PATH_IMAGE138
、轨道倾角对时间的变化率
Figure 380308DEST_PATH_IMAGE012
、轨道升交点赤经对时间的变化率
Figure 278994DEST_PATH_IMAGE014
、升交点角距余弦调和校正振幅C uc 、升交点角距正弦调和校正振幅C us 、轨道半径余弦调和校正振幅C rc 、轨道半径正弦调和校正振幅C rs 、轨道倾角余弦调和校正振幅C ic 和轨道倾角正弦调和校正振幅C is ；
B. 对全局卫星状态信息缓存中的当前卫星的状态信息有效标志进行判断：
若当前卫星的状态信息有效标志无效，则进行后续步骤；
若当前卫星的状态信息有效标志有效，则按照时间间隔抽样方式进行判断；
C. 对当前卫星的星历信息进行检测：
若当前卫星的星历信息有效（包括星历已经收齐、星历健康、星历龄期有效等），则进行后续步骤；
若当前卫星的星历信息无效，则设置当前卫星的状态信息有效标志为无效，将结果保存到全局卫星状态信息缓存中，并等待下一拍的调度计算；
D. 采用当前卫星的星历信息，计算当前时刻的卫星位置；具体包括如下步骤：
D1. 计算归一化时间t k 为
Figure 114095DEST_PATH_IMAGE016
，式中t为当前时刻；计算平均角速度n为
Figure 389218DEST_PATH_IMAGE139
，式中
Figure 767110DEST_PATH_IMAGE020
为地心引力常数；计算平近点角M k 为
Figure 836697DEST_PATH_IMAGE140
，式中n为计算得到的平均角速度；
D2. 以
Figure 424673DEST_PATH_IMAGE024
为初始值，采用算式
Figure 237909DEST_PATH_IMAGE141
进行迭代，直至精度满足设定的要求，得到最终的偏近点角E；E j 为第j次计算得到的偏近点角；
D3. 计算真近点角v k 为
Figure 470307DEST_PATH_IMAGE142
；计算升交点角距
Figure 710795DEST_PATH_IMAGE030
为
Figure 520488DEST_PATH_IMAGE143
；
D4. 采用如下算式计算摄动校正项：
Figure 137414DEST_PATH_IMAGE034
Figure 224319DEST_PATH_IMAGE036
Figure 25922DEST_PATH_IMAGE038
式中
Figure 198277DEST_PATH_IMAGE144
为升交点角距修正量；
Figure 353315DEST_PATH_IMAGE145
为卫星地心向径修正量；
Figure 560306DEST_PATH_IMAGE044
为卫星轨道倾角修正量；
D5. 计算升交点角距u k 为
Figure 736072DEST_PATH_IMAGE046
；计算卫星矢量长度r k 为
Figure 395723DEST_PATH_IMAGE146
；计算轨道倾角i k 为
Figure 479086DEST_PATH_IMAGE050
；
D6. 采用如下算式计算卫星轨道平面的位置
Figure 540583DEST_PATH_IMAGE052
和升交点赤经
Figure 559354DEST_PATH_IMAGE054
：
Figure 706302DEST_PATH_IMAGE147
Figure 330706DEST_PATH_IMAGE058
Figure 981130DEST_PATH_IMAGE148
式中
Figure 170803DEST_PATH_IMAGE062
为地球自转角速度；
D7. 采用如下算式计算地心地固直角坐标系中的卫星坐标
Figure 929680DEST_PATH_IMAGE064
为
Figure 230212DEST_PATH_IMAGE149
，
Figure 722DEST_PATH_IMAGE068
和
Figure 361296DEST_PATH_IMAGE150
；
E. 采用当前卫星的星历信息，计算当前时刻的卫星速度；具体为将地心地固直角坐标系中的卫星坐标
Figure 607469DEST_PATH_IMAGE151
对时间求导，从而得到当前时刻的卫星速度
Figure 446112DEST_PATH_IMAGE152
为：
Figure 336708DEST_PATH_IMAGE153
Figure 727238DEST_PATH_IMAGE076
Figure 70495DEST_PATH_IMAGE154
式中
Figure 978408DEST_PATH_IMAGE080
为升交点赤经
Figure 723510DEST_PATH_IMAGE082
对时间的变化率；
Figure 550521DEST_PATH_IMAGE084
为轨道倾角i k 对时间的变化率；
Figure 381073DEST_PATH_IMAGE086
为卫星轨道平面的位置
Figure 561519DEST_PATH_IMAGE088
对时间的变化率；
F. 采用当前卫星的星历信息，计算当前时刻的卫星加速度；具体为将步骤E得到的当前时刻的卫星速度对时间求导，从而得到当前时刻的卫星加速度
Figure 426707DEST_PATH_IMAGE155
为：
Figure 424619DEST_PATH_IMAGE156
Figure 476888DEST_PATH_IMAGE157
Figure 461025DEST_PATH_IMAGE158
式中
Figure 305353DEST_PATH_IMAGE098
为第一系数，且
Figure 349532DEST_PATH_IMAGE159
；
Figure 154677DEST_PATH_IMAGE102
为第二系数，且
Figure 942505DEST_PATH_IMAGE104
；
Figure 641339DEST_PATH_IMAGE106
为轨道倾角i k 对时间的二次求导；
Figure 856420DEST_PATH_IMAGE160
为卫星轨道平面的位置
Figure 883282DEST_PATH_IMAGE110
对时间的二次求导；
G. 计算完成后，设置当前卫星的状态信息有效标志为有效，并将计算结果进行保存，保存到全局卫星状态信息缓存中；
S3. 在步骤S2进行的同时，采用第二线程，结合二次运动模型理论，以二次曲线拟合的方式进行卫星状态信息外推，并进行最小二乘解算；
本步骤所采用的“以二次曲线拟合的方式进行卫星状态信息外推”的技术方案，是一种全新的计算方式；该方案基于卫星轨道理论模型和大量测试数据的统计结果，结合二次运动模型理论，实现了PVT解算的快速进行和优化；
具体实施时，步骤S3具体包括如下步骤：
a. 从全局卫星状态信息缓存中获取当前卫星的状态信息有效标志：
若当前卫星的状态信息有效标志为无效，则本步骤结束；并获取下一颗卫星的状态信息有效标志；
若当前卫星的状态信息有效标志为有效，则计算并得到当前卫星的外推时间间隔
Figure 325066DEST_PATH_IMAGE112
为
Figure 753773DEST_PATH_IMAGE161
，式中t pos 为当前PVT解算时刻，t s 为当前卫星的最近一次状态信息计算的时刻；
b. 设置时间阈值t th1（优选为10s）：
若当前卫星的外推时间间隔
Figure 139755DEST_PATH_IMAGE112
大于t th1时，判定当前卫星异常（包括失锁、星历无效等情况），不再外推当前卫星的状态信息，且不再使用当前卫星参与后续的解算过程；
若当前卫星的外推时间间隔
Figure 919492DEST_PATH_IMAGE112
小于或等于t th1时，采用如下二次曲线模型外推当前卫星在t pos 时刻的位置和速度：
Figure 173756DEST_PATH_IMAGE162
Figure 456970DEST_PATH_IMAGE163
Figure 279432DEST_PATH_IMAGE164
Figure 139941DEST_PATH_IMAGE165
Figure 807683DEST_PATH_IMAGE166
Figure 210982DEST_PATH_IMAGE126
式中上标s表示为当前卫星s；
Figure 204346DEST_PATH_IMAGE128
为当前卫星s在t s 时刻的位置；
Figure 286571DEST_PATH_IMAGE130
为当前卫星s在t s 时刻的速度；
Figure 23583DEST_PATH_IMAGE132
为当前卫星s在t s 时刻的加速度；
Figure 281389DEST_PATH_IMAGE134
为当前卫星s外推到t pos 时刻的位置；
Figure 304709DEST_PATH_IMAGE136
为当前卫星s外推到t pos 时刻的速度；
c. 重复步骤a~b，直至遍历完成所有卫星；
d. 将所有正常卫星的状态信息外推后，采用外推得到的状态信息进行后续的最小二乘解算；
本步骤在进行状态信息外推时，当外推时间间隔
Figure 280755DEST_PATH_IMAGE168
在几分钟之内时，采用二次曲线拟合计算所得的卫星位置和速度的误差分别小于10cm和1mm/s，这对定位结果产生的误差极其微小，可以忽略不计；而它的计算比直接用卫星星历参数来计算卫星位置和速度的方法要快约20倍，极大的提升了PVT解算的效率；
S4. 根据步骤S3得到的最小二乘解算结果，得到最终的定位结果，完成PVT解算。
在具体实施时，采用第一线程进行计算的优先级低于采用第二线程进行计算的优先级；同时，步骤S2计算得到卫星状态信息的计算结果后，将计算结果更新到全局卫星状态信息缓存中，供步骤S3进行卫星状态信息外推使用；在对全局卫星状态信息缓存进行任何操作时，均采用互斥量对全局卫星状态信息缓存进行防护，从而保证卫星状态信息的完整性、操作的可靠性和安全性。
以下结合一个实施例，对本发明的解算方法进行进一步说明：
导航接收机使用ARM A7处理器600MHz主频进行四系统PVT定位解算，四系统分别为北斗B1，GPS的L1，GLO的L1和GAL的E1，第一线程调度频度100Hz（要求每10ms计算一颗卫星的状态信息），四系统卫星最大数N为127，PVT定位解算频度为10Hz（要求1秒内解算10次，即每100ms解算一次），运行步骤如下：
（1） T1+100*K1 ms时刻，接收机跟踪了15颗北斗卫星、10颗GPS卫星、7颗GLO卫星、8颗GAL卫星，这些卫星的观测量数据组装到obs数据结构中；同时，所有卫星的导航电文数据全部接收完成，已组装到nav数据结构中。
（2） T1+100*K1 ms时刻，第一线程和第二线程都满足执行时机，各自启动执行。系统优先执行高优先级的第二线程，秒内第一拍定位解算使用之前已经计算并存储好的卫星状态信息集合进行外推，第一拍解算完成耗时3.53ms。
（3） T1+100*K1+3.53 ms时刻，开始执行第一线程，运行了0.16ms完成卫星s的状态信息计算。
（4） T1+100*K1+10 ms时刻，再次执行第一线程，运行了0.01ms跳过卫星s+1的计算，因为该卫星号没有被跟踪。
（5） T1+100*K1+20 ms时刻，再次执行第一线程，运行了0.17ms完成卫星s+2的状态信息计算。
……
（6） T1+100*K1+90 ms时刻，再次执行第一线程，运行了0.16ms完成卫星s+9的状态信息计算。
（7） T1+100*K2 ms时刻，秒内第二拍定位解算启动执行，再次组装obs数据，并使用已经更新了卫星s、卫星s+2、…、卫星s+9的状态信息的集合进行外推，第二拍解算完成耗时3.54ms。
（8）上述步骤（7）中，K2为K1的下一拍解算，把K2当作K1，继续重复步骤（1）。
（9） 步骤（2）中，在接收机开机的时候卫星状态信息集合是空的，第二线程定位解算得不到定位结果，导航接收机的第一个定位结果输出会在步骤（2）中K1迭代到某一个KN时给出。
（10） 步骤（6）和步骤（7）会存在同时操作全局卫星状态信息缓存的可能性，接收机引入互斥量防护技术消除了多线程同时操作的灾难性后果。
以下对传统的PVT解算方法和本发明所述的高性能PVT解算方法进行性能对比：
导航接收机使用ARM A7处理器600MHz主频进行四系统PVT定位解算，定位解算频度分别配置为10Hz、20Hz、50Hz和100Hz，对比传统方法和本发明所述方法的实验结果数据如下表1所示：
表1实验数据对比示意表
Figure 555879DEST_PATH_IMAGE170
从上表数据可以看出：使用传统解算方法最大只能支持50Hz解算频度，对于100Hz解算频度已经无法支持了。而使用本发明所述的高性能PVT解算方法，接收机系统的CPU使用率和解算总耗时都大大降低。在20Hz解算的情况下，本发明方法系统负担下降为传统方法的60%左右，PVT解算效率提升了62%左右，且能够支持到更高的100Hz频度解算，改善效果非常明显。
如图4所示为本发明的导航方法的方法流程示意图：本发明提供的这种包括了所述高性能PVT解算方法的导航方法，包括如下步骤：
（1）实时获取导航所需的数据信息；
（2）采用上述的高性能PVT解算方法进行实时PVT解算，得到实时的定位结果；
（3）根据步骤（2）得到的实时定位结果，完成最终的导航。

Claims (10)

1.一种高性能PVT解算方法，其特征在于包括如下步骤：
S1. 获取PVT解算所需要的输入数据；
S2. 根据步骤S1获取的输入数据，采用第一线程，以设定的调度频度、卫星号周期性方式和时间间隔抽样方式，将所有卫星状态信息的计算工作分摊到不同的时隙，并逐个卫星进行卫星状态信息的计算；
S3. 在步骤S2进行的同时，采用第二线程，结合二次运动模型理论，以二次曲线拟合的方式进行卫星状态信息外推，并进行最小二乘解算；
S4. 根据步骤S3得到的最小二乘解算结果，得到最终的定位结果，完成PVT解算。
2.根据权利要求1所述的高性能PVT解算方法，其特征在于采用第一线程进行计算的优先级低于采用第二线程进行计算的优先级；步骤S2计算得到卫星状态信息的计算结果后，将计算结果更新到全局卫星状态信息缓存中，供步骤S3进行卫星状态信息外推使用；在对全局卫星状态信息缓存进行任何操作时，均采用互斥量对全局卫星状态信息缓存进行防护，从而保证卫星状态信息的完整性、操作的可靠性和安全性。
3.根据权利要求2所述的高性能PVT解算方法，其特征在于步骤S1所述的获取PVT解算所需要的输入数据，具体包括obs观测量数据和nav导航电文数据；其中，nav导航电文数据包括eph星历数据，satclk卫星钟差数据，iono电离层数据和tgd群延迟数据；然后将获取的输入数据按照设定的数据结构和格式进行组装。
4.根据权利要求3所述的高性能PVT解算方法，其特征在于步骤S2所述的卫星号周期性方式，具体包括如下步骤：
总卫星数N为
Figure 608575DEST_PATH_IMAGE002
，其中N 1为GPS卫星的数量，N 2为GLO卫星的数量，N 3为GAL卫星的数量，N 4为北斗卫星的数量；在设定的调度频度下，每一拍仅计算一颗卫星的状态信息；当计算完最后一颗卫星的状态信息后，再从第一颗卫星开始周期性的进行计算。
5.根据权利要求4所述的高性能PVT解算方法，其特征在于步骤S2所述的时间间隔抽样方式，具体包括如下步骤：
已经计算过第i颗卫星的状态信息后，当再次轮到对第i颗卫星进行计算时，只有当前时刻与上次进行状态信息计算的时刻的时间间隔大于设定值时，才对第i颗卫星进行状态信息的计算；否则，直接跳过第i颗卫星的状态信息计算过程。
6.根据权利要求5所述的高性能PVT解算方法，其特征在于步骤S2所述的进行卫星状态信息的计算，具体包括如下步骤：
A. 获取卫星的星历数据信息；
B. 对全局卫星状态信息缓存中的当前卫星的状态信息有效标志进行判断：
若当前卫星的状态信息有效标志无效，则进行后续步骤；
若当前卫星的状态信息有效标志有效，则按照时间间隔抽样方式进行判断；
C. 对当前卫星的星历信息进行检测：
若当前卫星的星历信息有效，则进行后续步骤；
若当前卫星的星历信息无效，则设置当前卫星的状态信息有效标志为无效，将结果保存到全局卫星状态信息缓存中，并等待下一拍的调度计算；
D. 采用当前卫星的星历信息，计算当前时刻的卫星位置；
E. 采用当前卫星的星历信息，计算当前时刻的卫星速度；
F. 采用当前卫星的星历信息，计算当前时刻的卫星加速度；
G. 计算完成后，设置当前卫星的状态信息有效标志为有效，并将计算结果进行保存，保存到全局卫星状态信息缓存中。
7.根据权利要求6所述的高性能PVT解算方法，其特征在于步骤A所述的卫星星历的数据信息，具体包括星历参考时间t oe 、卫星轨道长半轴a s 的平方根值
Figure 943741DEST_PATH_IMAGE004
、轨道偏心率e s 、星历参考时间t oe 时刻的轨道倾角i 0、周内时为0时的轨道升交点赤经
Figure 266138DEST_PATH_IMAGE006
、轨道近地角距
Figure 79374DEST_PATH_IMAGE008
、星历参考时间t oe 时刻的平近点角M 0、平均运动角速度校正值
Figure 311772DEST_PATH_IMAGE010
、轨道倾角对时间的变化率
Figure 676894DEST_PATH_IMAGE012
、轨道升交点赤经对时间的变化率
Figure 96374DEST_PATH_IMAGE014
、升交点角距余弦调和校正振幅C uc 、升交点角距正弦调和校正振幅C us 、轨道半径余弦调和校正振幅C rc 、轨道半径正弦调和校正振幅C rs 、轨道倾角余弦调和校正振幅C ic 和轨道倾角正弦调和校正振幅C is 。
8.根据权利要求7所述的高性能PVT解算方法，其特征在于所述的步骤D，具体包括如下步骤：
D1. 计算归一化时间t k 为
Figure 837934DEST_PATH_IMAGE016
，式中t为当前时刻；计算平均角速度n为
Figure 659259DEST_PATH_IMAGE018
，式中
Figure 70649DEST_PATH_IMAGE020
为地心引力常数；计算平近点角M k 为
Figure 370568DEST_PATH_IMAGE022
，式中n为计算得到的平均角速度；
D2. 以
Figure 791185DEST_PATH_IMAGE024
为初始值，采用算式
Figure 732596DEST_PATH_IMAGE026
进行迭代，直至精度满足设定的要求，得到最终的偏近点角E；E j 为第j次计算得到的偏近点角；
D3. 计算真近点角v k 为
Figure 580467DEST_PATH_IMAGE028
；计算升交点角距
Figure 364752DEST_PATH_IMAGE030
为
Figure 323481DEST_PATH_IMAGE032
；
D4. 采用如下算式计算摄动校正项：
Figure 384978DEST_PATH_IMAGE034
Figure 403749DEST_PATH_IMAGE036
Figure 144172DEST_PATH_IMAGE038
式中
Figure 906592DEST_PATH_IMAGE040
为升交点角距修正量；
Figure 681650DEST_PATH_IMAGE042
为卫星地心向径修正量；
Figure 871323DEST_PATH_IMAGE044
为卫星轨道倾角修正量；
D5. 计算升交点角距u k 为
Figure 505566DEST_PATH_IMAGE046
；计算卫星矢量长度r k 为
Figure 806098DEST_PATH_IMAGE048
；计算轨道倾角i k 为
Figure 701241DEST_PATH_IMAGE050
；
D6. 采用如下算式计算卫星轨道平面的位置
Figure 796236DEST_PATH_IMAGE052
和升交点赤经
Figure 652197DEST_PATH_IMAGE054
：
Figure 881053DEST_PATH_IMAGE056
Figure 771648DEST_PATH_IMAGE058
Figure 303124DEST_PATH_IMAGE060
式中
Figure 646380DEST_PATH_IMAGE062
为地球自转角速度；
D7. 采用如下算式计算地心地固直角坐标系中的卫星坐标
Figure 413348DEST_PATH_IMAGE064
为
Figure 892871DEST_PATH_IMAGE066
，
Figure 595248DEST_PATH_IMAGE068
和
Figure 550434DEST_PATH_IMAGE070
；
所述的步骤E，具体为将地心地固直角坐标系中的卫星坐标
Figure 996459DEST_PATH_IMAGE064
对时间求导，从而得到当前时刻的卫星速度
Figure 861647DEST_PATH_IMAGE072
为：
Figure 734925DEST_PATH_IMAGE074
Figure 660898DEST_PATH_IMAGE076
Figure 910614DEST_PATH_IMAGE078
式中
Figure 630308DEST_PATH_IMAGE080
为升交点赤经
Figure 799121DEST_PATH_IMAGE082
对时间的变化率；
Figure 338687DEST_PATH_IMAGE084
为轨道倾角i k 对时间的变化率；
Figure 126514DEST_PATH_IMAGE086
为卫星轨道平面的位置
Figure 700715DEST_PATH_IMAGE088
对时间的变化率；
所述的步骤F，具体为将步骤E得到的当前时刻的卫星速度对时间求导，从而得到当前时刻的卫星加速度
Figure 40429DEST_PATH_IMAGE090
为：
Figure 332871DEST_PATH_IMAGE092
Figure 658810DEST_PATH_IMAGE094
Figure 353096DEST_PATH_IMAGE096
式中
Figure 598133DEST_PATH_IMAGE098
为第一系数，且
Figure 112291DEST_PATH_IMAGE100
；
Figure 241921DEST_PATH_IMAGE102
为第二系数，且
Figure 915348DEST_PATH_IMAGE104
；
Figure 472231DEST_PATH_IMAGE106
为轨道倾角i k 对时间的二次求导；
Figure 473685DEST_PATH_IMAGE108
为卫星轨道平面的位置
Figure 407006DEST_PATH_IMAGE110
对时间的二次求导。
9.根据权利要求8所述的高性能PVT解算方法，其特征在于所述的步骤S3，具体包括如下步骤：
a. 从全局卫星状态信息缓存中获取当前卫星的状态信息有效标志：
若当前卫星的状态信息有效标志为无效，则本步骤结束；并获取下一颗卫星的状态信息有效标志；
若当前卫星的状态信息有效标志为有效，则计算并得到当前卫星的外推时间间隔
Figure 934939DEST_PATH_IMAGE112
为
Figure 662724DEST_PATH_IMAGE114
，式中t pos 为当前PVT解算时刻，t s 为当前卫星的最近一次状态信息计算的时刻；
b. 设置时间阈值t th1：
若当前卫星的外推时间间隔
Figure 151474DEST_PATH_IMAGE112
大于t th1时，判定当前卫星异常，不再外推当前卫星的状态信息，且不再使用当前卫星参与后续的解算过程；
若当前卫星的外推时间间隔
Figure 888486DEST_PATH_IMAGE112
小于或等于t th1时，采用如下二次曲线模型外推当前卫星在t pos 时刻的位置和速度：
Figure 270925DEST_PATH_IMAGE116
Figure 169611DEST_PATH_IMAGE118
Figure 145658DEST_PATH_IMAGE120
Figure 420781DEST_PATH_IMAGE122
Figure 657727DEST_PATH_IMAGE124
Figure 992894DEST_PATH_IMAGE126
式中上标s表示为当前卫星s；
Figure 925078DEST_PATH_IMAGE128
为当前卫星s在t s 时刻的位置；
Figure 131456DEST_PATH_IMAGE130
为当前卫星s在t s 时刻的速度；
Figure 98275DEST_PATH_IMAGE132
为当前卫星s在t s 时刻的加速度；
Figure 604342DEST_PATH_IMAGE134
为当前卫星s外推到t pos 时刻的位置；
Figure 289402DEST_PATH_IMAGE136
为当前卫星s外推到t pos 时刻的速度；
c. 重复步骤a~b，直至遍历完成所有卫星；
d. 将所有正常卫星的状态信息外推后，采用外推得到的状态信息进行后续的最小二乘解算。
10.一种包括了权利要求1~9之一所述的高性能PVT解算方法的导航方法，其特征在于包括如下步骤：
（1）实时获取导航所需的数据信息；
（2）采用权利要求1~9之一所述的高性能PVT解算方法进行实时PVT解算，得到实时的定位结果；
（3）根据步骤（2）得到的实时定位结果，完成最终的导航。
CN202210597377.4A 2022-05-30 2022-05-30 高性能pvt解算方法及导航方法 Active CN114675311B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210597377.4A CN114675311B (zh) 2022-05-30 2022-05-30 高性能pvt解算方法及导航方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210597377.4A CN114675311B (zh) 2022-05-30 2022-05-30 高性能pvt解算方法及导航方法

Publications (2)

Publication Number Publication Date
CN114675311A CN114675311A (zh) 2022-06-28
CN114675311B true CN114675311B (zh) 2022-08-12

Family

ID=82081085

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210597377.4A Active CN114675311B (zh) 2022-05-30 2022-05-30 高性能pvt解算方法及导航方法

Country Status (1)

Country Link
CN (1) CN114675311B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116011561B (zh) * 2023-03-28 2023-06-30 武汉大学 基于神经网络的信息外推方法、装置、设备及存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109001771A (zh) * 2018-06-04 2018-12-14 北京未来导航科技有限公司 导航卫星和低轨卫星实时钟差确定及预报方法和系统
CN112235033A (zh) * 2020-09-24 2021-01-15 中国科学院国家授时中心 一种基于rdss通信链路的gnss共视系统
CN113359510A (zh) * 2021-06-04 2021-09-07 北京理工大学 北斗卫星导航系统信号模拟器数据实时仿真系统及方法
CN113534221A (zh) * 2021-09-16 2021-10-22 长沙海格北斗信息技术有限公司 高精度高频度接收机解算方法、解算系统及导航接收机

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9057780B2 (en) * 2013-04-18 2015-06-16 California Institute Of Technology Real-time and post-processed orbit determination and positioning
WO2018009088A1 (en) * 2016-07-04 2018-01-11 Llc "Topcon Positioning Systems" Gnss positioning system and method using multiple processing threads

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109001771A (zh) * 2018-06-04 2018-12-14 北京未来导航科技有限公司 导航卫星和低轨卫星实时钟差确定及预报方法和系统
CN112235033A (zh) * 2020-09-24 2021-01-15 中国科学院国家授时中心 一种基于rdss通信链路的gnss共视系统
CN113359510A (zh) * 2021-06-04 2021-09-07 北京理工大学 北斗卫星导航系统信号模拟器数据实时仿真系统及方法
CN113534221A (zh) * 2021-09-16 2021-10-22 长沙海格北斗信息技术有限公司 高精度高频度接收机解算方法、解算系统及导航接收机

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
A Parallel Approach for Multi-GNSS Ultra-Rapid Orbit Determination;Chunhua Jiang 等;《Remote Sensing》;20210901;第13卷(第17期);1-18 *
A parallel GPS satellite positioning algorithm based on broadcast ephemeris;Si Zhang 等;《2015 8th International Congress on Image and Signal Processing (CISP)》;20160218;7-11 *
卫星轨道递推的GPU集成式并行加速方法;孔繁泽 等;《哈尔滨工业大学学报》;20210630;第53卷(第6期);13-20 *

Also Published As

Publication number Publication date
CN114675311A (zh) 2022-06-28

Similar Documents

Publication Publication Date Title
WO2020233158A1 (zh) 一种基于智能手机的高精度单点定位方法与装置
CN114675311B (zh) 高性能pvt解算方法及导航方法
Stamp et al. Comparison of oral 25-hydroxycholecalciferol, vitamin D, and ultraviolet light as determinants of circulating 25-hydroxyvitamin D
CN101950025B (zh) 用于局域增强系统的数据质量监测方法
CN1135265A (zh) 一种确定的完整性监视的外推导航装置
JP5359187B2 (ja) 長期予測軌道データの信頼性判定方法、長期予測軌道データの提供方法及び情報提供装置
CN104215974A (zh) 一种卫星导航系统的完好性监测可用性确定方法
CN110553653B (zh) 基于多源数据驱动的航天器轨道确定方法
CN110007317A (zh) 一种选星优化的高级接收机自主完好性监测方法
CN110275799B (zh) 账务系统不停机日切点快照余额的方法
Mullin et al. Validation of the REVEAL Prognostic Equation and Risk Score Calculator in Incident Systemic Sclerosis–Associated Pulmonary Arterial Hypertension
JP5359188B2 (ja) 長期予測軌道データの信頼性判定方法、長期予測軌道データの提供方法及び情報提供装置
CN110501733A (zh) 一种自适应格网vrs生成与服务方法
Banville et al. Ionospheric monitoring using" integer-levelled" observations
CN108873023A (zh) 一种提高定位精度的观测量的处理方法
CN107238848B (zh) 卫星导航系统可用性的评估方法及装置
CN101866021A (zh) 并行化Abel变换大气参数数据处理方法
CN103033827B (zh) 一种卫星位置解算方法
CN114002713A (zh) 一种卫星轨道参数递推处理预报系统
CN111781641B (zh) 一种被动震源定位方法、系统、终端及可读存储介质
CN110011719A (zh) 一种利用gps定轨数据的星上轨道获取控制方法
CN113504557A (zh) 面向实时应用的gps频间钟差新预报方法
CN116299561A (zh) 一种降低电离层不规则体对精密单点定位影响的方法
CN113341686B (zh) 一种单星多历元校时方法、装置、计算机设备和存储介质
CN111999750B (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