CN106802153B - 基于单频导航原始观测量地面处理的高精度测定轨方法 - Google Patents
基于单频导航原始观测量地面处理的高精度测定轨方法 Download PDFInfo
- Publication number
- CN106802153B CN106802153B CN201710055239.2A CN201710055239A CN106802153B CN 106802153 B CN106802153 B CN 106802153B CN 201710055239 A CN201710055239 A CN 201710055239A CN 106802153 B CN106802153 B CN 106802153B
- Authority
- CN
- China
- Prior art keywords
- satellite
- frequency
- observed quantity
- real time
- perturbation
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/24—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for cosmonautical navigation
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Signal Processing (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Astronomy & Astrophysics (AREA)
- Automation & Control Theory (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明提供了一种基于单频导航原始观测量地面处理的高精度测定轨方法,其包括步骤一、获取原始观测量并下传、步骤二、准实时轨道动力学平滑处理、步骤三、事后单频精密定轨处理,星载导航接收机将实时获取的导航卫星原始观测量下传地面,地面采用准实时轨道动力学平滑处理和事后单频精密定轨处理方法,将卫星测定轨精度从10米量级提高到0.2米,成功解决了卫星自主测定轨精度低从而限制了卫星应用的问题,从地面处理上取得突破,取得了降低卫星复杂度、有效提高卫星测定轨精度、提高了卫星应用能力等有益效果,对于我国后续配置单频BDS、GPS接收机的卫星创造了一个新的提高测定轨精度的方法,该发明在本领域内应用将十分广泛。
Description
技术领域
本发明涉及一种测定轨方法,具体地,涉及一种基于单频导航原始观测量地面处理的高精度测定轨方法。
背景技术
卫星上搭载高动态单频GPS(美国全球定位***)、BDS(中国北斗卫星导航***)接收机,能够实时获取GPS和BDS导航信号在第一频率上的测距码、载波相位和多普勒等观测值,同时根据接收到的导航电文信息利用原始观测量自主确定卫星位置和速度,但是,由于单频观测值中电离层延迟难以精确模型化,导航电文中发布的导航卫星轨道和钟差精度有限,星载接收机计算能力有限,实时定轨算法中难以考虑精确的观测模型和动力学模型,因而星上定轨获取的卫星位置和速度精度较差,难以满足卫星应用的要求,因此,必须在提高卫星测定轨精度上提供一种新的方法解决上述问题,目前没有发现同本发明类似技术的说明或报道,也尚未收集到国内外类似的资料。
发明内容
针对现有技术中的缺陷,本发明的目的是提供一种基于单频导航原始观测量地面处理的高精度测定轨方法,星载导航接收机将实时获取的导航卫星原始观测量下传地面,地面采用准实时轨道动力学平滑处理和事后单频精密定轨处理方法,将卫星测定轨精度从10米量级提高到0.2米,成功解决了卫星自主测定轨精度低从而限制了卫星应用的问题,从地面处理上取得突破,取得了降低卫星复杂度、有效提高卫星测定轨精度、提高了卫星应用能力等有益效果,对于我国后续配置单频BDS、GPS接收机的卫星创造了一个新的提高测定轨精度的方法,该发明在本领域内应用将十分广泛。
根据本发明的一个方面,提供一种基于单频导航原始观测量地面处理的高精度测定轨方法,其特征在于,其包括以下步骤:
步骤一,获取原始观测量并下传;
步骤二,准实时轨道动力学平滑处理;
步骤三,事后单频精密定轨处理。
优选地,所述步骤一获取原始观测量并下传中,卫星装载高动态单频BDS、GPS接收机,以1秒的高采样率实时获取BDS和GPS导航信号在第一频率上的测距码、载波相位和多普勒观测值,同时根据接收到的导航电文信息利用原始观测量自主确定卫星位置和速度,导航接收机按照约定的格式将原始观测量、实时定位数据组帧,存储到固态存储器里,通过高速数据传输通道下传地面,同时下传的还有卫星姿态数据。
优选地,所述步骤二准实时轨道动力学平滑处理中,准实时轨道动力学平滑处理采用轨道拟合方法,具体是将星上确定的卫星位置作为观测量来拟合卫星运动方程,求解卫星在参考时刻的位置速度和指定的动力学模型参数,在卫星动力学运动方程中考虑作用在卫星上的力有地球球形二体引力、地球非球形引力摄动、太阳系天体引力摄动、相对论效应摄动、地球固体潮汐摄动、地球海洋潮汐摄动、固体极移潮汐摄动、海洋极移潮汐摄动、大气阻力摄动、太阳光压摄动,此外,为了吸收未精确考虑、难以模型化、未知的其它力,在卫星轨道切向和法向上附加正余弦周期性的经验加速度。
优选地,所述步骤三事后单频精密定轨处理中,事后单频精密定轨处理采用星载单频BDS、GPS确定低轨卫星精密轨道技术,将星上下传的导航原始观测量中的伪距、载波相位作为观测值建立观测方程,通过轨道积分获取卫星动力学轨道作为运动方程,同时求解观测方程和运动方程达到观测值和卫星运动状态的最佳匹配,在卫星动力学运动方程中考虑作用在卫星上的力与准实时轨道动力学平滑中考虑的一致,观测方程中考虑的因素为导航卫星钟误差、接收机钟误差、电离层延迟、相对论效应、天线相位缠绕、天线相位中心偏差和变化,同一时刻同一导航卫星同一频率的伪距和载波相位观测值中的电离层延迟一阶项大小相同符号相反,因而采用伪距与载波相位求和再取一半的组合观测消除电离层一阶项,忽略微小的电离层二阶项影响,接收机钟误差在每一个观测时刻附加一个参数估计,组合观测值的距离偏差在每一个导航卫星跟踪弧段内估计一个参数,大气阻力规格化系数每3小时估计一个,周期性经验加速度分别在轨道切向和法向上每3小时估计一个余弦振幅和一个正弦振幅,同时估计接收机天线在Z方向上的相位中心偏差,参数估计中根据待估参数的特性通过带有参数递推、消去、恢复的最小二乘方法加快数据处理过程。
与现有技术相比,本发明具有如下的有益效果:星载导航接收机将实时获取的导航卫星原始观测量下传地面,地面采用准实时轨道动力学平滑处理和事后单频精密定轨处理方法,将卫星测定轨精度从10米量级提高到0.2米,成功解决了卫星自主测定轨精度低从而限制了卫星应用的问题,从地面处理上取得突破,取得了降低卫星复杂度、有效提高卫星测定轨精度、提高了卫星应用能力等有益效果,对于我国后续配置单频BDS、GPS接收机的卫星创造了一个新的提高测定轨精度的方法,该发明在本领域内应用将十分广泛。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1为基于单频导航原始观测量地面处理的高精度测定轨方法过程示意图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进。这些都属于本发明的保护范围。
如图1所示,本发明基于单频导航原始观测量地面处理的高精度测定轨方法包括如下的步骤:
步骤一,获取原始观测量并下传:卫星装载高动态单频BDS、GPS接收机,以1秒的高采样率实时获取BDS和GPS导航信号在第一频率上的测距码、载波相位和多普勒等观测值,同时根据接收到的导航电文信息利用原始观测量自主确定卫星位置和速度,实现的测定轨精度是10米量级,导航接收机按照约定的格式将原始观测量、实时定位数据组帧,存储到固态存储器里,通过高速数据传输通道下传地面,同时下传的还有卫星姿态数据等。
步骤二,准实时轨道动力学平滑处理:准实时轨道动力学平滑处理采用轨道拟合方法,具体是将星上确定的卫星位置作为观测量来拟合卫星运动方程,求解卫星在参考时刻的位置速度和指定的动力学模型参数,在卫星动力学运动方程中考虑作用在卫星上的力有地球球形二体引力、地球非球形引力摄动、太阳系天体引力摄动、相对论效应摄动、地球固体潮汐摄动、地球海洋潮汐摄动、固体极移潮汐摄动、海洋极移潮汐摄动、大气阻力摄动、太阳光压摄动,此外,为了吸收未精确考虑、难以模型化、未知的其它力,在卫星轨道切向和法向上附加正余弦周期性的经验加速度,通过上述方法,将卫星测定轨精度从10米量级提高到3米量级。
步骤三,事后单频精密定轨处理:事后单频精密定轨处理采用星载单频BDS、GPS确定低轨卫星精密轨道技术,将星上下传的导航原始观测量中的伪距、载波相位作为观测值建立观测方程,通过轨道积分获取卫星动力学轨道作为运动方程,同时求解观测方程和运动方程达到观测值和卫星运动状态的最佳匹配,在卫星动力学运动方程中考虑作用在卫星上的力与准实时轨道动力学平滑中考虑的一致,观测方程中考虑的因素为导航卫星钟误差、接收机钟误差、电离层延迟、相对论效应、天线相位缠绕、天线相位中心偏差和变化等,同一时刻同一导航卫星同一频率的伪距和载波相位观测值中的电离层延迟一阶项大小相同符号相反,因而采用伪距与载波相位求和再取一半的组合观测消除电离层一阶项,忽略微小的电离层二阶项影响,接收机钟误差在每一个观测时刻附加一个参数估计,组合观测值的距离偏差在每一个导航卫星跟踪弧段内估计一个参数,大气阻力规格化系数每3小时估计一个,周期性经验加速度分别在轨道切向和法向上每3小时估计一个余弦振幅和一个正弦振幅,同时估计接收机天线在Z方向上的相位中心偏差,参数估计中根据待估参数的特性通过带有参数递推、消去、恢复的最小二乘方法加快数据处理过程,通过上述方法,将卫星测定轨精度提高到0.2米。
从上述方案的描述中可以看出:本发明基于单频导航原始观测量地面处理的高精度测定轨方法,星载导航接收机将实时获取的导航卫星原始观测量下传地面,地面采用准实时轨道动力学平滑处理和事后单频精密定轨处理方法,将卫星测定轨精度从10米量级提高到0.2米,成功解决了卫星自主测定轨精度低从而限制了卫星应用的问题,从地面处理上取得突破,取得了降低卫星复杂度、有效提高卫星测定轨精度、提高了卫星应用能力等有益效果,对于我国后续配置单频BDS、GPS接收机的卫星创造了一个新的提高测定轨精度的方法,该发明在本领域内应用将十分广泛。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变形或修改,这并不影响本发明的实质内容。
Claims (1)
1.一种基于单频导航原始观测量地面处理的高精度测定轨方法,其特征在于,其包括以下步骤:
步骤一,获取原始观测量并下传;
步骤二,准实时轨道动力学平滑处理;
步骤三,事后单频精密定轨处理;
所述步骤一获取原始观测量并下传中,卫星装载高动态单频BDS、GPS接收机,以1秒的高采样率实时获取BDS和GPS导航信号在第一频率上的测距码、载波相位和多普勒观测值,同时根据接收到的导航电文信息利用原始观测量自主确定卫星位置和速度,单频BDS、GPS接收机按照约定的格式将原始观测量、实时定位数据组帧,存储到固态存储器里,通过高速数据传输通道下传地面,同时下传的还有卫星姿态数据;
所述步骤二准实时轨道动力学平滑处理中,准实时轨道动力学平滑处理采用轨道拟合方法,具体是将星上确定的卫星位置作为观测量来拟合卫星运动方程,求解卫星在参考时刻的位置速度和指定的动力学模型参数,在卫星动力学运动方程中考虑作用在卫星上的力有地球球形二体引力、地球非球形引力摄动、太阳系天体引力摄动、相对论效应摄动、地球固体潮汐摄动、地球海洋潮汐摄动、固体极移潮汐摄动、海洋极移潮汐摄动、大气阻力摄动、太阳光压摄动,此外,为了吸收未精确考虑、难以模型化、未知的其它力,在卫星轨道切向和法向上附加正余弦周期性的经验加速度;
所述步骤三事后单频精密定轨处理中,事后单频精密定轨处理采用单频BDS、GPS接收机确定低轨卫星精密轨道技术,将星上下传的导航原始观测量中的伪距、载波相位作为观测值建立观测方程,通过轨道积分获取卫星动力学轨道作为运动方程,同时求解观测方程和运动方程达到观测值和卫星运动状态的最佳匹配,在卫星动力学运动方程中考虑作用在卫星上的力与准实时轨道动力学平滑中考虑的一致,观测方程中考虑的因素为导航卫星钟误差、接收机钟误差、电离层延迟、相对论效应、天线相位缠绕、天线相位中心偏差和变化,同一时刻同一导航卫星同一频率的伪距和载波相位观测值中的电离层延迟一阶项大小相同符号相反,因而采用伪距与载波相位求和再取一半的组合观测消除电离层一阶项,忽略微小的电离层二阶项影响,接收机钟误差在每一个观测时刻附加一个参数估计,组合观测值的距离偏差在每一个导航卫星跟踪弧段内估计一个参数,大气阻力规格化系数每3小时估计一个,周期性经验加速度分别在轨道切向和法向上每3小时估计一个余弦振幅和一个正弦振幅,同时估计接收机天线在Z方向上的相位中心偏差,参数估计中根据待估参数的特性通过带有参数递推、消去、恢复的最小二乘方法加快数据处理过程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710055239.2A CN106802153B (zh) | 2017-01-24 | 2017-01-24 | 基于单频导航原始观测量地面处理的高精度测定轨方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710055239.2A CN106802153B (zh) | 2017-01-24 | 2017-01-24 | 基于单频导航原始观测量地面处理的高精度测定轨方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106802153A CN106802153A (zh) | 2017-06-06 |
CN106802153B true CN106802153B (zh) | 2019-09-17 |
Family
ID=58988239
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710055239.2A Active CN106802153B (zh) | 2017-01-24 | 2017-01-24 | 基于单频导航原始观测量地面处理的高精度测定轨方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106802153B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109782313A (zh) * | 2019-01-17 | 2019-05-21 | 上海卫星工程研究所 | 地面精密单点定位数据处理方法 |
CN112444828A (zh) * | 2019-08-30 | 2021-03-05 | 河北地质大学 | 一种仅基于伪距观测值确定gnss卫星轨道的方法 |
CN110988932B (zh) * | 2019-12-11 | 2022-02-01 | 武汉大学 | 一种提高星载gps接收机实时钟差解算精度的方法 |
CN112305575B (zh) * | 2020-09-25 | 2024-04-02 | 北京空间飞行器总体设计部 | 一种高轨sar卫星精密定轨*** |
CN117724128B (zh) * | 2024-02-07 | 2024-04-30 | 中南大学 | 一种低轨卫星轨道预报方法、***、终端及介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102305630A (zh) * | 2011-05-17 | 2012-01-04 | 哈尔滨工业大学 | 基于扩展卡尔曼滤波的sar卫星自主定轨方法 |
CN103363994A (zh) * | 2012-04-09 | 2013-10-23 | 陈刘成 | 仅基于无线电载波相位观测的卫星精密定轨技术 |
CN103592653A (zh) * | 2013-10-17 | 2014-02-19 | 中国科学院光电研究院 | 用于局部地区单频卫星导航用户的电离层延迟修正方法 |
CN103675861A (zh) * | 2013-11-18 | 2014-03-26 | 航天恒星科技有限公司 | 一种基于星载gnss多天线的卫星自主定轨方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6463366B2 (en) * | 2000-03-10 | 2002-10-08 | Schafer Corp | Attitude determination and alignment using electro-optical sensors and global navigation satellites |
-
2017
- 2017-01-24 CN CN201710055239.2A patent/CN106802153B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102305630A (zh) * | 2011-05-17 | 2012-01-04 | 哈尔滨工业大学 | 基于扩展卡尔曼滤波的sar卫星自主定轨方法 |
CN103363994A (zh) * | 2012-04-09 | 2013-10-23 | 陈刘成 | 仅基于无线电载波相位观测的卫星精密定轨技术 |
CN103592653A (zh) * | 2013-10-17 | 2014-02-19 | 中国科学院光电研究院 | 用于局部地区单频卫星导航用户的电离层延迟修正方法 |
CN103675861A (zh) * | 2013-11-18 | 2014-03-26 | 航天恒星科技有限公司 | 一种基于星载gnss多天线的卫星自主定轨方法 |
Non-Patent Citations (2)
Title |
---|
星载GPS伪距观测量重建及定轨精度分析;彭冬菊等;《遥感学报》;20101231;第14卷(第6期);正文第1-5节 |
用单频GPS数据实现低轨卫星动力学法定轨研究;张德成等;《导航定位学报》;20140331;第2卷(第1期);正文第1-3节 |
Also Published As
Publication number | Publication date |
---|---|
CN106802153A (zh) | 2017-06-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106802153B (zh) | 基于单频导航原始观测量地面处理的高精度测定轨方法 | |
CN108919634B (zh) | 一种北斗三频非差非组合观测值时间传递***及方法 | |
WO2022007211A1 (zh) | 一种基于gnss的实时高精度波浪测量方法及装置 | |
CN102033236B (zh) | 一种卫星导航位置速度联合估计方法 | |
CN104316947B (zh) | Gnss/ins超紧组合导航装置及相对导航*** | |
CN110376621A (zh) | 一种基于北斗三号B2b信号的卫星定位方法及装置 | |
CN104597471A (zh) | 面向时钟同步多天线gnss接收机的定向测姿方法 | |
CN204101732U (zh) | 一种gnss新体制信号捕获装置 | |
CN110986879A (zh) | 一种电力线塔倾斜实时监测方法及*** | |
CN104280746A (zh) | 一种惯性辅助gps的深组合半实物仿真*** | |
CN113253314B (zh) | 一种低轨卫星间时间同步方法及*** | |
CN104536016A (zh) | 一种gnss新体制信号捕获装置及方法 | |
CN106842263A (zh) | 基于导航信号空时处理与矢量跟踪相结合的卫星导航方法 | |
CN107544075A (zh) | 基于精密单点定位与调和分析估计海潮负荷位移参数方法 | |
CN103529482B (zh) | 一种高精度确定载体动态加速度的方法 | |
CN103454665A (zh) | 一种双差gps/sins组合导航姿态测量方法 | |
CN106767922B (zh) | 一种海洋卫星星载高度计在轨绝对定标***及定标方法 | |
CN102253396A (zh) | 一种高动态gps载波环跟踪方法 | |
CN107193023A (zh) | 一种具有闭式解的高精度北斗卫星***单点定位方法 | |
CN109613582B (zh) | 一种车载实时单频米级伪距定位方法 | |
CN106768179A (zh) | 基于连续运行gnss站信噪比数据的潮位的测量方法 | |
CN106597492A (zh) | 卫星导航接收机及其抗远近效应的方法和室内定位方法 | |
CN103543454A (zh) | 一种嵌入在移动通讯网中的卫星定轨*** | |
CN107132423B (zh) | 一种探测电离层电子密度总数的方法及装置 | |
CN210242756U (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 |