CN109991633A - 一种低轨卫星实时定轨方法 - Google Patents
一种低轨卫星实时定轨方法 Download PDFInfo
- Publication number
- CN109991633A CN109991633A CN201910164665.9A CN201910164665A CN109991633A CN 109991633 A CN109991633 A CN 109991633A CN 201910164665 A CN201910164665 A CN 201910164665A CN 109991633 A CN109991633 A CN 109991633A
- Authority
- CN
- China
- Prior art keywords
- epoch
- observed quantity
- error
- cycle slip
- carrier phase
- 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
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/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/27—Acquisition or tracking or demodulation of signals transmitted by the system creating, predicting or correcting ephemeris or almanac data within the receiver
-
- 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/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/30—Acquisition or tracking or demodulation of signals transmitted by the system code related
-
- 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
Abstract
本发明公开了一种低轨卫星实时定轨方法,包括以下步骤:从星载GNSS接收机获取待定位历元的双频测码伪距和载波相位观测量;对获取的双频测码伪距和载波相位观测量分别进行误差修正,包括电离层延迟、相对论效应、天线相位中心误差、地球自转效应;对修正误差后的载波相位观测量进行周跳和粗差的实时探测与修复;在使用经过处理的载波相位观测值对测码伪距观测量进行平滑;使用地面上注的IGS超快速星历计算导航星的位置和钟差信息;将上一历元的定位结果作为初值,使用扩展卡尔曼滤波器计算当前历元的低轨卫星定位结果;判断定位任务是否结束,若没结束,则返回第一步骤继续进行计算。本发明在适应星上处理器有限计算资源达到较高定位精度。
Description
技术领域
本发明涉及实时定位技术,更具体地说,尤其涉及一种低轨卫星实时定轨方法。
背景技术
当前,诸如高精度测绘、导航增强等领域都需要低轨卫星的高精度定位数据,随着需求的不断扩大,越来越多的应用对定位的实时性和高精度都提出了要求。目前国内外学者进行的研究大部分都是采用以下两种方法:方法一是以测码伪距或侧相位距作为观测量,使用广播星历计算导航星位置和钟差,再结合简化的动力学方程,使用扩展卡尔曼滤波器进行计算;方法二除了将观测量换成载波相位以外,其它步骤与方法一大体相同。但是它们都有各自的不足之处:方法一定位的精度不理想;方法二由于涉及到模糊度的解算,状态向量维度大大增加,在星载计算机上应用的难度很大。
本发明将有效地解决上述缺陷,面向低轨卫星的实际工程应用,充分考虑精度要求以及星载计算机计算能力的限制,解决兼顾计算精度与计算效率的星载算法关键问题。。
发明内容
本发明的目的在于克服现有技术的不足,提供一种兼顾计算精度与计算效率的低轨卫星实时定轨方法,在对星载GNSS接收机观测量进行数据预处理的同时,使用地面提前上注的IGS超快速预测星历代替精度不高的广播星历来计算导航卫星的位置和钟差,结合低轨卫星的动力学方程,使用扩展卡尔曼滤波器对目标卫星进行定位。
本发明是根据以下技术方案实现的:
一种低轨卫星实时定轨方法,其特征在于,包括以下步骤:
步骤S1:从星载GNSS接收机获取待定位历元的双频测码伪距和载波相位观测量;
步骤S2:对获取的双频测码伪距和载波相位观测量分别进行误差修正,包括:电离层延迟、相对论效应、天线相位中心误差、地球自转效应;
步骤S3:对修正误差以后的载波相位观测量进行周跳和粗差的实时探测与修复;
步骤S4:在使用经过处理的载波相位观测值对测码伪距观测量进行平滑;
步骤S5:使用地面上注的IGS超快速星历计算导航星的位置和钟差信息;
步骤S6:将上一历元的定位结果作为初值,使用扩展卡尔曼滤波器计算当前历元的低轨卫星定位结果;
步骤S7:判断定位任务是否结束,若没有结束,则返回步骤S1继续进行计算。
上述技术方案中,步骤S3具体包括:
步骤S301:使用M-W组合计算整周模糊度互差项Nδ
其中,N1,N2分别为f1,f2频率上的整周模糊度,L1,L2分别为f1,f2频率上的载波相位观测值,ρ1,ρ2分别为f1,f2频率上的测码伪距观测值;
步骤S302:计算前i个历元Nδ的平均值<Nδ>i及其均方根误差σi,为了保证方法的实时性,采用以下方式计算:
其中,Nδi为第i个历元的整周模糊度互差项;
步骤S303:使用上述公式计算历元i的整周模糊度互差项Nδi、历元i+1的整周模糊度互差项Nδ(i+1)以及前i-1个历元整周模糊度互差项的均值<Nδ>i-1以后,通过|Nδi-<Nδ>i-1|≥σ和|Nδ(i+1)-Nδi|<1探测周跳和粗差,如果这两个判断条件同时满足,认为在历元i-1和i之间已经发生了周跳;如果有|Nδi-<Nδ>i-1|≥σ和|Nδ(i+1)-Nδi|≥1同时成立,认为历元i时刻的相位突变是粗差造成的,公式中的σ为根据星载GNSS接收机的实际情况设定好的阈值;
步骤S304:如果判断出来有周跳存在,进行如下计算:
接着分别在周跳发生前后的两个历元联立Nδ和W两个方程,分别解出周跳发生前后历元在f1,f2两个频率上的整周模糊度N1,N2,进而解出周跳发生前后历元整周模糊度的变化量,也就是周跳,接着对周跳发生之后的连续观测历元中对观测量进行修正;
如果判断出有粗差存在,则将粗差对应的观测量进行剔除。
上述技术方案中,所述步骤S4进一步包括:如果出现意外情况星上无法获取超快速星历,则使用广播星历计算导航星位置和钟差。
上述技术方案中,步骤S6进一步包括:滤波器中解算的状态向量包括三维位置向量、三维速度向量以及接收机钟差,以平滑后的测码伪距作为观测量,首先使用上一历元的定位结果或事先设定的初值作为卡尔曼滤波器的初始状态向量;然后使用4阶龙格库塔法对低轨卫星的精密动力学方程进行数值积分得到位置、速度向量的一步预测值,同时把接收机钟差当作白噪声模型来计算一步预测值;最后使用经过载波相位平滑以后的测码伪距观测值对状态向量的一步预测值进行修正。
与现有技术相比,本发明具有如下的有益效果:
本发明可以在星上计算资源消耗较小的前提下,完成较高精度的低轨卫星定位工作,可以适用于低轨卫星星载GNSS接收机的实时定位。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1为本发明的低轨卫星实时定位方法的流程图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变化和改进。这些都属于本发明的保护范围。
如图1所示,一种低轨卫星实时定轨方法,它包括以下步骤:
步骤S1:从星载GNSS接收机获取待定位历元的双频测码伪距和载波相位观测量;
步骤S2:对获取的双频测码伪距和载波相位观测量分别进行误差修正,包括:电离层延迟、相对论效应、天线相位中心误差、地球自转效应;
步骤S3:对修正误差以后的载波相位观测量进行周跳和粗差的实时探测与修复;具体步骤为:
步骤S3具体包括:
步骤S301:使用M-W组合计算整周模糊度互差项Nδ
其中,N1,N2分别为f1,f2频率上的整周模糊度,L1,L2分别为f1,f2频率上的载波相位观测值,ρ1,ρ2分别为f1,f2频率上的测码伪距观测值;
步骤S302:计算前i个历元Nδ的平均值<Nδ>i及其均方根误差σi,为了保证方法的实时性,采用以下方式计算:
其中,Nδi为第i个历元的整周模糊度互差项;
步骤S303:使用上述公式计算历元i的整周模糊度互差项Nδi、历元i+1的整周模糊度互差项Nδ(i+1)以及前i-1个历元整周模糊度互差项的均值<Nδ>i-1以后,通过|Nδi-<Nδ>i-1|≥σ和|Nδ(i+1)-Nδi|<1探测周跳和粗差,如果这两个判断条件同时满足,认为在历元i-1和i之间已经发生了周跳;如果有|Nδi-<Nδ>i-1|≥σ和|Nδ(i+1)-Nδi|≥1同时成立,认为历元i时刻的相位突变是粗差造成的,公式中的σ为根据星载GNSS接收机的实际情况设定好的阈值;
步骤S304:如果判断出来有周跳存在,进行如下计算:
接着分别在周跳发生前后的两个历元联立Nδ和W两个方程,分别解出周跳发生前后历元在f1,f2两个频率上的整周模糊度N1,N2,进而解出周跳发生前后历元整周模糊度的变化量,也就是周跳,接着对周跳发生之后的连续观测历元中对观测量进行修正;
如果判断出有粗差存在,则将粗差对应的观测量进行剔除。
步骤S4:使用经过处理的载波相位观测值对测码伪距观测量进行平滑;
考虑到实时计算的需求,使用以下公式来进行处理:
初始条件为:
步骤S5:使用地面上注的IGS超快速星历计算导航星的位置和钟差信息;
特别需要说明的是,以GPS系统为例,IGS发布的精密星历有以下几种:事后精密星历以及钟差的精度最高,但是不能实时获得;快速星历及钟差同样不能实时获取;超快速星历分为观测和预报两种,其中观测的同样不能实时获得,预报的可以提供未来24h的星历和钟差信息,适用于本发明。
由于精密星历和钟差相邻两个点的时间间隔为15min,因此需要使用九阶拉格朗日插值法计算导航卫星信号发射时刻的位置和钟差。
步骤S6:将上一历元的定位结果作为初值,使用扩展卡尔曼滤波器计算当前历元的低轨卫星定位结果;其中,如果是初始历元则改为预先设定的初值。
状态向量为三维位置向量、三维速度向量以及接收机钟差。
观测向量为经过载波相位平滑以后的测码伪距观测量。
扩展卡尔曼滤波器的状态方程为低轨卫星的动力学方程,方程中考虑的几种摄动力模型如表1所示:
表1摄动力模型以及参数设置
在计算一步预测值时,使用4阶龙格库塔法对动力学方程进行数值积分。
为了减小计算量,关于位置和速度的状态转移矩阵为:
其中
对于钟差,则认为其在相邻历元间是不变的,因此对应到转移矩阵中的元素为1,并且与速度和位置无关。
滤波器中的观测方程为:
其中,xi,yi,zi分别对应导航卫星在信号发射时刻在ECEF坐标系下的坐标;x,y,z为接收机在ECEF坐标系下的坐标;cδt是接收机钟差(导航星钟差已经被修正)。
观测矩阵为:
其中
为接收机与GPS星之间的距离。
状态噪声协方差矩阵设置为:
其中Qy表示位置和速度的噪声协方差矩阵,在这里认为不同轴上的分量之间相互无关,但是同一轴上的位置和速度是相关的。同一轴上的位置和速度的噪声协方差矩阵形式为:
在对观测噪声R阵进行设置时,认为不同类型的观测值和不同GPS星之间的观测值是互不相关的,所以R阵表示为:
步骤S7:判断定位任务是否结束,若没有结束,则返回步骤S1继续进行计算。
如果在当前历元计算完成以后还要接着进行定位,则需要等待一个设定好的历元间隔才可以继续计算,这个历元间隔的设定必须保证星载计算机在这段时间里可以完成步骤S1~S7的所有工作,因此需要针对不同的硬件环境作出调整。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变化或修改,这并不影响本发明的实质内容。在不冲突的情况下,本申请的实施例和实施例中的特征可以任意相互组合。
Claims (4)
1.一种低轨卫星实时定轨方法,其特征在于,包括以下步骤:
步骤S1:从星载GNSS接收机获取待定位历元的双频测码伪距和载波相位观测量;
步骤S2:对获取的双频测码伪距和载波相位观测量分别进行误差修正,包括:电离层延迟、相对论效应、天线相位中心误差、地球自转效应;
步骤S3:对修正误差以后的载波相位观测量进行周跳和粗差的实时探测与修复;
步骤S4:在使用经过处理的载波相位观测值对测码伪距观测量进行平滑;
步骤S5:使用地面上注的IGS超快速星历计算GNSS星的位置和钟差信息;
步骤S6:将上一历元的定位结果作为初值,使用扩展卡尔曼滤波器计算当前历元的低轨卫星定位结果;
步骤S7:判断定位任务是否结束,若没有结束,则返回步骤S1继续进行计算。
2.根据权利要求1所述的一种低轨卫星实时定轨方法,其特征在于,步骤S3具体包括:
步骤S301:使用M-W组合计算整周模糊度互差项Nδ
其中,N1,N2分别为f1,f2频率上的整周模糊度,L1,L2分别为f1,f2频率上的载波相位观测值,ρ1,ρ2分别为f1,f2频率上的测码伪距观测值;
步骤S302:计算前i个历元Nδ的平均值<Nδ>i及其均方根误差σi,为了保证方法的实时性,采用以下方式计算:
其中,Nδi为第i个历元的整周模糊度互差项;
步骤S303:使用上述公式计算历元i的整周模糊度互差项Nδi、历元i+1的整周模糊度互差项Nδ(i+1)以及前i-1个历元整周模糊度互差项的均值<Nδ>i-1以后,通过|Nδi-<Nδ>i-1|≥σ和|Nδ(i+1)-Nδi|<1探测周跳和粗差,如果这两个判断条件同时满足,认为在历元i-1和i之间已经发生了周跳;如果有|Nδi-<Nδ>i-1|≥σ和|Nδ(i+1)-Nδi|≥1同时成立,认为历元i时刻的相位突变是粗差造成的,公式中的σ为根据星载GNSS接收机的实际情况设定好的阈值;
步骤S304:如果判断出来有周跳存在,进行如下计算:
接着分别在周跳发生前后的两个历元联立Nδ和W两个方程,分别解出周跳发生前后历元在f1,f2两个频率上的整周模糊度N1,N2,进而解出周跳发生前后历元整周模糊度的变化量,也就是周跳,接着对周跳发生之后的连续观测历元中对观测量进行修正;
如果判断出有粗差存在,则将粗差对应的观测量进行剔除。
3.根据权利要求1所述的一种低轨卫星实时定轨方法,其特征在于,所述步骤S4进一步包括:如果出现意外情况星上无法获取超快速星历,则使用广播星历计算导航星位置和钟差。
4.根据权利要求1所述的一种低轨卫星实时定轨方法,其特征在于,步骤S6进一步包括:滤波器中解算的状态向量包括三维位置向量、三维速度向量以及接收机钟差,以平滑后的测码伪距作为观测量,首先使用上一历元的定位结果或事先设定的初值作为卡尔曼滤波器的初始状态向量;然后使用4阶龙格库塔法对低轨卫星的精密动力学方程进行数值积分得到位置、速度向量的一步预测值,同时把接收机钟差当作白噪声模型来计算一步预测值;最后使用经过载波相位平滑以后的测码伪距观测值对状态向量的一步预测值进行修正。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910164665.9A CN109991633A (zh) | 2019-03-05 | 2019-03-05 | 一种低轨卫星实时定轨方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910164665.9A CN109991633A (zh) | 2019-03-05 | 2019-03-05 | 一种低轨卫星实时定轨方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109991633A true CN109991633A (zh) | 2019-07-09 |
Family
ID=67130139
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910164665.9A Pending CN109991633A (zh) | 2019-03-05 | 2019-03-05 | 一种低轨卫星实时定轨方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109991633A (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110515097A (zh) * | 2019-09-02 | 2019-11-29 | 江苏省测绘工程院 | 应用于基准站的gnss卫星观测粗差剔除方法和装置 |
CN110631588A (zh) * | 2019-09-23 | 2019-12-31 | 电子科技大学 | 一种基于rbf网络的无人机视觉导航定位方法 |
CN110673175A (zh) * | 2019-09-16 | 2020-01-10 | 西安空间无线电技术研究所 | 一种基于gnss广播星历的高轨卫星高精度自主定轨方法 |
CN110879407A (zh) * | 2019-12-12 | 2020-03-13 | 北京眸星科技有限公司 | 一种基于完好性风险模型的卫星导航观测量新息检测方法 |
CN110988932A (zh) * | 2019-12-11 | 2020-04-10 | 武汉大学 | 一种提高星载gps接收机实时钟差解算精度的方法 |
CN110988941A (zh) * | 2019-12-27 | 2020-04-10 | 北京遥测技术研究所 | 一种高精度实时绝对定轨方法 |
CN110988917A (zh) * | 2019-12-10 | 2020-04-10 | 中国科学院国家授时中心 | 一种卫星轨道机动状态实时监测方法 |
CN111045062A (zh) * | 2019-11-29 | 2020-04-21 | 航天恒星科技有限公司 | 一种基于电磁星的星基电离层反演方法 |
CN111123322A (zh) * | 2019-12-11 | 2020-05-08 | 上海卫星工程研究所 | 星载gnss接收机的观测值实时数据预处理方法、系统、介质及设备 |
CN111487660A (zh) * | 2020-04-24 | 2020-08-04 | 北京航空航天大学 | 一种高精度实时微纳卫星集群导航算法 |
CN111505679A (zh) * | 2020-04-20 | 2020-08-07 | 中国科学院国家空间科学中心 | 一种基于星载gnss的leo初轨确定方法 |
CN112666583A (zh) * | 2020-12-15 | 2021-04-16 | 上海卫星工程研究所 | 适应gnss接收机输出状态的单拍轨道递推方法及系统 |
CN112764067A (zh) * | 2020-12-23 | 2021-05-07 | 深圳创维数字技术有限公司 | Gps卫星星历数据的获取方法、装置、车载及可读存储介质 |
CN117289319A (zh) * | 2023-11-27 | 2023-12-26 | 北京李龚导航科技有限公司 | 基于阵列天线的载波相位定向方法、装置及电子设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008038919A1 (en) * | 2006-09-29 | 2008-04-03 | Electronics And Telecommunications Research Institute | Precise orbit determination system and method using gps data and galileo data |
CN101334458A (zh) * | 2008-06-03 | 2008-12-31 | 电子科技大学 | 一种卫星导航定位中载波相位的周跳修复方法 |
CN105203110A (zh) * | 2015-08-28 | 2015-12-30 | 中国科学院空间应用工程与技术中心 | 一种基于大气阻力模型补偿的低轨卫星轨道预报方法 |
CN107229061A (zh) * | 2017-07-18 | 2017-10-03 | 武汉大学 | 一种基于低轨卫星的星地差分实时精密定位方法 |
-
2019
- 2019-03-05 CN CN201910164665.9A patent/CN109991633A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008038919A1 (en) * | 2006-09-29 | 2008-04-03 | Electronics And Telecommunications Research Institute | Precise orbit determination system and method using gps data and galileo data |
CN101334458A (zh) * | 2008-06-03 | 2008-12-31 | 电子科技大学 | 一种卫星导航定位中载波相位的周跳修复方法 |
CN105203110A (zh) * | 2015-08-28 | 2015-12-30 | 中国科学院空间应用工程与技术中心 | 一种基于大气阻力模型补偿的低轨卫星轨道预报方法 |
CN107229061A (zh) * | 2017-07-18 | 2017-10-03 | 武汉大学 | 一种基于低轨卫星的星地差分实时精密定位方法 |
Non-Patent Citations (4)
Title |
---|
GEOFFREY BLEWITT: "An automatic editing algorithm for GPS data", 《GEOPHYSICAL RESEARCH LETTERS》 * |
吴显兵: "星载GPS低轨卫星几何法定轨及动力学平滑方法研究", 《中国优秀博硕士学位论文全文数据库(硕士)基础科学辑》 * |
涂先勤: "星载双频GPS精密定轨及精度评估方法研究", 《中国优秀硕士学位论文全文数据库基础科学辑》 * |
郑永艾等: "在轨卫星 GPS 导航实时高精度定轨控制仿真", 《计算机仿真》 * |
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110515097B (zh) * | 2019-09-02 | 2021-06-22 | 江苏省测绘工程院 | 应用于基准站的gnss卫星观测粗差剔除方法和装置 |
CN110515097A (zh) * | 2019-09-02 | 2019-11-29 | 江苏省测绘工程院 | 应用于基准站的gnss卫星观测粗差剔除方法和装置 |
CN110673175A (zh) * | 2019-09-16 | 2020-01-10 | 西安空间无线电技术研究所 | 一种基于gnss广播星历的高轨卫星高精度自主定轨方法 |
CN110631588A (zh) * | 2019-09-23 | 2019-12-31 | 电子科技大学 | 一种基于rbf网络的无人机视觉导航定位方法 |
CN111045062A (zh) * | 2019-11-29 | 2020-04-21 | 航天恒星科技有限公司 | 一种基于电磁星的星基电离层反演方法 |
CN110988917A (zh) * | 2019-12-10 | 2020-04-10 | 中国科学院国家授时中心 | 一种卫星轨道机动状态实时监测方法 |
CN110988917B (zh) * | 2019-12-10 | 2021-09-10 | 中国科学院国家授时中心 | 一种卫星轨道机动状态实时监测方法 |
CN110988932A (zh) * | 2019-12-11 | 2020-04-10 | 武汉大学 | 一种提高星载gps接收机实时钟差解算精度的方法 |
CN111123322A (zh) * | 2019-12-11 | 2020-05-08 | 上海卫星工程研究所 | 星载gnss接收机的观测值实时数据预处理方法、系统、介质及设备 |
CN111123322B (zh) * | 2019-12-11 | 2022-07-29 | 上海卫星工程研究所 | 星载gnss接收机的观测值实时数据预处理方法、系统、介质及设备 |
CN110879407A (zh) * | 2019-12-12 | 2020-03-13 | 北京眸星科技有限公司 | 一种基于完好性风险模型的卫星导航观测量新息检测方法 |
CN110879407B (zh) * | 2019-12-12 | 2023-08-18 | 北京眸星科技有限公司 | 一种基于完好性风险模型的卫星导航观测量新息检测方法 |
CN110988941A (zh) * | 2019-12-27 | 2020-04-10 | 北京遥测技术研究所 | 一种高精度实时绝对定轨方法 |
CN111505679A (zh) * | 2020-04-20 | 2020-08-07 | 中国科学院国家空间科学中心 | 一种基于星载gnss的leo初轨确定方法 |
CN111487660A (zh) * | 2020-04-24 | 2020-08-04 | 北京航空航天大学 | 一种高精度实时微纳卫星集群导航算法 |
CN112666583A (zh) * | 2020-12-15 | 2021-04-16 | 上海卫星工程研究所 | 适应gnss接收机输出状态的单拍轨道递推方法及系统 |
CN112764067A (zh) * | 2020-12-23 | 2021-05-07 | 深圳创维数字技术有限公司 | Gps卫星星历数据的获取方法、装置、车载及可读存储介质 |
CN112764067B (zh) * | 2020-12-23 | 2023-07-07 | 深圳创维数字技术有限公司 | Gps卫星星历数据的获取方法、装置、车载及可读存储介质 |
CN117289319A (zh) * | 2023-11-27 | 2023-12-26 | 北京李龚导航科技有限公司 | 基于阵列天线的载波相位定向方法、装置及电子设备 |
CN117289319B (zh) * | 2023-11-27 | 2024-04-05 | 北京李龚导航科技有限公司 | 基于阵列天线的载波相位定向方法、装置及电子设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109991633A (zh) | 一种低轨卫星实时定轨方法 | |
CN106168672B (zh) | 一种gnss多模单频rtk周跳探测方法及装置 | |
CN102736090B (zh) | 位置计算方法和位置计算装置 | |
CN104076348B (zh) | 一种雷达超视距基线无源协同定位方法 | |
CN109709591A (zh) | 一种面向智能终端的gnss高精度定位方法 | |
CN102230971A (zh) | Gps多天线测姿方法 | |
CN103199949B (zh) | 伪卫星定位系统中的零基线自闭环星地时间同步方法 | |
CN110412638A (zh) | 一种低成本三天线gnss rtk定位及测姿方法 | |
CN104459740A (zh) | 一种定位终端的高精度位置差分定位方法 | |
CN110058282A (zh) | 一种基于双频gnss智能手机的ppp高精度定位方法 | |
CN101750600A (zh) | 一种用于测绘与导航的实时高精度定位方法 | |
CN105891860A (zh) | 一种基于误差分离模式的gnss区域伪距差分增强定位方法 | |
CN101614802A (zh) | 一种导航卫星姿态测量方法 | |
CN105510942A (zh) | 一种基于卡尔曼滤波的gps单点定位系统 | |
CN108827322A (zh) | 一种多星协同测向定位观测系统优化设计与评估方法 | |
CN113253314A (zh) | 一种低轨卫星间时间同步方法及系统 | |
CN105093251A (zh) | Gnss接收机静态模式下的高精度相对定位方法 | |
CN104614739A (zh) | 基于抗干扰滤波器的北斗多频接收机信号联合跟踪方法 | |
Havyarimana et al. | Toward accurate intervehicle positioning based on GNSS pseudorange measurements under non-Gaussian generalized errors | |
CN116755126B (zh) | 一种基于三维模型映射匹配的北斗实时精准定位方法 | |
CN105204049A (zh) | 一种基于三维转一维的载波相位整周模糊度搜索的定位方法 | |
CN115343734B (zh) | 一种基于双线性内插半球模型的gnss变形监测方法 | |
CN109085626B (zh) | 一种定位方法及装置 | |
CN104765057A (zh) | 卫星定位方法、卫星定位装置与计算机可读取媒体 | |
Yue et al. | A Fusion Localization Algorithm with Adaptive Kalman Gain for Port Container Seamless 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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190709 |