浙江大学学报(工学版), 2020, 54(3): 589-596 doi: 10.3785/j.issn.1008-973X.2020.03.020

航空航天技术

基于北斗B3频点的低轨卫星实时定轨性能评估

余鑫,, 金小军,, 莫仕明, 张伟, 徐兆斌, 金仲和

Performance assessment of BeiDou B3-based real-time orbit determination for LEO satellites

YU Xin,, JIN Xiao-jun,, MO Shi-ming, ZHANG Wei, XU Zhao-bin, JIN Zhong-he

通讯作者: 金小军,男,副教授. orcid.org/0000-0002-8804-9890. E-mail: axemaster@zju.edu.cn

收稿日期: 2019-01-28  

Received: 2019-01-28  

作者简介 About authors

余鑫(1995—),男,硕士生,从事星载GPS/BDS定轨技术研究.orcid.org/0000-0002-6991-6549.E-mail:21724023@zju.edu.cn , E-mail:21724023@zju.edu.cn

摘要

基于自主研制的搭载于浙江大学皮星三号(ZDPS-3)任务的星载GPS/北斗双模接收机,在北斗B3频点上开展低轨卫星实时定轨仿真试验,对基于北斗(特别是B3频点)的实时定轨性能进行评估. 结果表明:B3频点的测距伪码码速率高,其抗噪声性能优于B1频点,有利于提升北斗导航系统下的低轨卫星定轨精度. 在建立观测模型和定轨算法模型并检验接收机实测数据质量的基础上,利用导航信号模拟器建立半物理仿真实验平台. 实验结果表明,在仅使用北斗二号14颗导航卫星的条件下,利用B3频点能获得明显优于B1频点的实时定轨精度,并且与基于处于完全运行状态的GPS的定轨精度相当.

关键词: 星载GNSS接收机 ; 北斗B3频点 ; 低轨卫星实时定轨 ; 半物理仿真 ; 皮纳卫星

Abstract

Real-time orbit determination simulations were conducted using the space-borne GPS/BDS (BeiDou navigation satellite system) receiver developed for the ZDPS-3 mission. The orbit determination performance based on BDS, especially the B3 frequency, was evaluated. As results, with the higher chip rate, the anti-noise performance of pseudo-noise code of B3 was obviously better than that of B1, which can be utilized to achieve a better orbit determination accuracy. Based on the establishment of the observation model and the orbit determination algorithm model, a hardware-in-loop simulation experimental platform was established after the quality of the receiver’s field observation data was evaluated. The simulation results show that, under the same conditions, the real-time orbit determination accuracy of the B3 frequency is better than that of B1, and comparable to that of GPS in full operation.

Keywords: space-borne GNSS receiver ; BeiDou B3 frequency ; real-time orbit determination ; hardware-in-loop simulation ; pico-satellite

PDF (1687KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

余鑫, 金小军, 莫仕明, 张伟, 徐兆斌, 金仲和. 基于北斗B3频点的低轨卫星实时定轨性能评估. 浙江大学学报(工学版)[J], 2020, 54(3): 589-596 doi:10.3785/j.issn.1008-973X.2020.03.020

YU Xin, JIN Xiao-jun, MO Shi-ming, ZHANG Wei, XU Zhao-bin, JIN Zhong-he. Performance assessment of BeiDou B3-based real-time orbit determination for LEO satellites. Journal of Zhejiang University(Engineering Science)[J], 2020, 54(3): 589-596 doi:10.3785/j.issn.1008-973X.2020.03.020

长期以来,基于星载导航接收机的低轨卫星定轨主要依赖GPS系统[1-2]. 随着北斗卫星导航系统(BeiDou navigation satellite system,BDS)的稳步建设和不断完善,单独使用BDS进行定轨已经成为可能. 截至2018年12月,北斗三号完成19颗卫星发射组网,预计到2020年底完成30颗卫星的发射与调试,全面建成北斗三号导航系统. 与GPS仅采用MEO轨道的导航卫星相比,北斗系统采用MEO、GEO、IGSO三种轨道卫星组成的混合星座,几何分布更优,可视时间更长. 采用独立自主的北斗导航系统对低轨卫星进行定轨,对于保障我国航天任务具有重要意义.

北斗导航系统提供多个频点的导航信号. 2018年2月,北斗官方公开了B3信号的空间信号接口控制文件,为全球用户提供公开服务. B3频点测距码的码速率达到10.23 Mcps,为GPS C/A码的10倍,BDS B1频点的5倍. 理论上更高的码速率具有更优的抗噪声性能,从而能获得更高的伪距观测精度,有利于提升低轨卫星的定轨精度. 基于BDS进行低轨卫星定轨目前已有研究,Li等[3-8]对搭载GPS/BDS双模接收机的风云三号卫星定轨进行研究,在接收机最多只能同时跟踪6颗北斗导航卫星的条件下,可以实现事后高精度定轨,在实时定轨方面,佘丽丽[9]使用风云三号C星数据,在亚太区域内使用BDS可以达到小于15 m的定轨精度;熊超等[10]采用模拟器仿真数据,仅使用BDS进行实时定轨仿真测试,得到的三轴误差在5 m以内. 已有研究尚未将B3频点信号用于低轨卫星定轨,定轨精度偏低,不足以满足高精度实时定轨的需求.

本文基于自主研制的搭载于浙江大学皮星三号(ZDPS-3)任务的星载GPS/BDS接收机,对基于北斗(特别是B3频点)的实时定轨性能进行评估,以期充分发挥B3频点伪码观测精度高的优势,提升低轨卫星实时定轨性能.

1. 定轨理论模型

1.1. 基本观测模型

伪距与载波相位是GNSS导航定位中最主要的观测量. 综合考虑各项误差后,可以得到伪距观测量ρ与载波相位观测量Ф的观测方程:

$\rho = \rho _{\rm{r}}^{\rm{s}} + c\left[ {\delta \left( {{t_{\rm{r}}}} \right) - \delta \left( {{t^{\rm{s}}}} \right)} \right] + {I_\rho } + {\rho _{{\rm{pco}}}} + {\varepsilon _\rho },$

$\varPhi = \rho _{\rm{r}}^{\rm{s}} + c\left[ {\delta \left( {{t_{\rm{r}}}} \right) - \delta \left( {{t^{\rm{s}}}} \right)} \right] + {I_\varPhi } + B + {\varPhi _{{\rm{pco}}}} + {\varepsilon _\varPhi }.$

式中:ρrs为接收机与导航星之间的几何距离,c为光速,δtr)为接收机钟差,δts)为导航卫星钟差,IρIФ分别为伪距、载波相位的电离层延迟,Фpcoρpco分别为天线相位中心偏差值(phase center offset,PCO)的误差,B为载波相位模糊度,ερεФ分别为伪距与载波相位的测量噪声. 对于低轨卫星,其运行轨道位于对流层之上,因此可以忽略低轨卫星受到的对流层延迟影响.

1.2. 伪距观测误差

接收机的码跟踪环常采用非相干DLL,非相干前减后功率是最常用的鉴相器之一. 在不考虑多径和其他干扰的情况下,码环的测量误差主要包括由热噪声所致的码相位抖动和动态应力误差两部分. 若σtDLL代表由热噪声所致的码相位测量误差均方差,则对于非相干前减后功率鉴相器,以伪码码片数为单位的σtDLL值可通过下式[11-12]进行计算:

$ \sigma _{\rm{tDLL}}^2 = \left\{ {\begin{array}{*{20}{c}} {\dfrac{{{B}_{\text{L}}}}{2C/{{N}_{0}}}D\left[ 1+\dfrac{2}{\left( 2-D \right){{T}_{\text{coh}}}C/{{N}_{0}}} \right], D \geqslant \dfrac{{\text{π}}}{{{B_{\rm{\rm{fe}}}}{T_{\rm{C}}}}};}\\ \dfrac{{{B}_{\text{L}}}}{2C/{{N}_{0}}}\left[ \dfrac{1}{{{B}_{\rm fe}}{{T}_{\rm{C}}}}+\dfrac{{{B}_{\rm fe}}{{T}_{\rm{C}}}}{\text{π} -1}{{\left( D-\dfrac{1}{{{B}_{\rm fe}}{{T}_{\rm{C}}}} \right)}^{2}} \right]\times \\ \left[ 1+\dfrac{2}{\left( 2-D \right){{T}_{\rm{coh}}}C/{{N}_{0}}} \right],\dfrac{1}{{{B}_{\rm fe}}{{T}_{\rm{C}}}} < D <\dfrac{\text{π} }{{{B}_{\rm fe}}{{T}_{\rm{C}}}};\\ {\dfrac{{{B_{\rm{L}}}}}{{2C/{N_0}}}\dfrac{1}{{{B_{{\rm{\rm{fe}}}}}{T_{\rm{C}}}}}\left( {1 + \dfrac{1}{{{T_{\rm{coh}}}C/{N_0}}}} \right),D \leqslant \dfrac{1}{{{B_{\rm{fe}}}{T_{\rm{C}}}}}.} \end{array}} \right. $

式中:Bfe为射频前端带宽,TC为伪码码宽,D为前后相关器间距,BL为环路噪声带宽,C/N0为载噪比,Tcoh为相干积分时间.

在相同参数条件下,B3码由热噪声引起的测量误差均方差理论上比B1小5倍.例如在本文所用的接收机中,B3频点的射频前端带宽Bfe为20 MHz,B1频点为4 MHz,前后相关器间距D为1个码片,在环路噪声带宽BL、载噪比C/N0与相干积分时间Tcoh相同的条件下,B1、B3频点由热噪声导致的以码片相位为度量的跟踪误差相同. 而B3的码片周期是B1的1/5,故换算为距离单位后B3频点由热噪声导致的跟踪误差小于B1. 需要指出的是,除了热噪声,决定DLL跟踪精度的误差源还应包括动态应力,两者共同影响最终的环路跟踪精度.

1.3. 消电离层组合观测模型

在低轨卫星定轨中,电离层延迟是影响定轨精度的主要误差来源之一. 对于双频观测数据,星载接收机常采用双频伪距组合的方式,可以有效消除电离层一阶项的影响;对于单频接收机数据,常采用Klobuchar等电离层模型[13],但是由于低轨卫星运行高度往往处于电离层峰值之上,且电离层模型精度受到太阳活动强弱的影响,直接将这些模型应用到低轨卫星定轨上,往往会带入不同程度的模型误差[14-15]. 因此,单频点定轨较为常用的方法是采用组合观测数据,以消除电离层效应带来的误差,常用的有GRAPHIC消电离层组合(group and phase ionospheric correction,GRAPHIC).

1.3.1. 双频伪距组合

伪距在通过电离层时的延迟量与其载波频率的平方成反比,利用这一特性,可构建组合观测值以消除电离层的影响:

${\rho _{{\rm{IF}}}} = \frac{{f_1^2}}{{f_1^2 - f_2^2}}{\rho _1} - \frac{{f_2^2}}{{f_1^2 - f_2^2}}{\rho _2}.$

式中:ρIF为伪距组合观测量,ρ1ρ2分别为载波频率f1f2上的伪距观测量,组合后的伪距观测噪声标准差为

${\sigma _{{\rm{IF}}}} = \sqrt {\frac{{{\gamma ^2}}}{{{{\left( {\gamma - 1} \right)}^2}}}\sigma _1^2 + \frac{1}{{{{\left( {\gamma - 1} \right)}^2}}}\sigma _2^2} .$

式中:γ=(f1/ f22σ1σ2分别为各自频点上的伪距观测噪声标准差. 可见,伪距组合的方式虽然能消除电离层延迟,却不可避免地放大了观测噪声,组合后的定轨精度甚至有可能低于单频的定轨精度.

1.3.2. GRAPHIC组合

根据同一导航信号频率上伪距和载波相位观测值的电离层延迟一阶项大小相等、符号相反的特性,可通过GRAPHIC组合的方式消除电离层影响[16],其组合方式为伪距观测量ρ与载波相位观测量Ф的数学平均,观测方程为

$\begin{split} {\rho _{\rm{G}}} =\left( {\varPhi + \rho } \right)/2 = & \rho _{\rm{r}}^{\rm{s}} + c\left[ {\delta \left( {{t_{\rm{r}}}} \right) - \delta \left( {{t^{\rm{s}}}} \right)} \right] + \\ & {B_{\rm{G}}} + \rho _{{\rm{G}}, {\rm{pco}}} + {\varepsilon _{\rm{G}}}, \end{split} $

式中:εG为GRAPHIC组合后的测量噪声, ${\rho _{{\rm{G}},{\rm{pco}}}}$BG分别为组合后的天线相位中心偏差值和模糊度.

GRAPHIC组合测量噪声的标准差如下式:

$\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!{\sigma _{\rm{G}}} = 0.5{\left( {\sigma _\rho ^2 + \sigma _\varPhi ^2} \right)^{1/2}}.$

由于载波相位观测值噪声远小于伪距噪声,式(7)可近似为σG≈0.5σρ,即观测噪声约为伪距噪声的一半. 然而GRAPHIC组合方式会引入载波相位的模糊度,本文的处理方法是将模糊度BG作为扩展卡尔曼滤波(extended kalman filter,EKF)[17-20]的一个状态量来估计,在EKF滤波算法中更新,其初始值设置为BG=(Фρ)/2. 由于天线相位PCO的误差值较小,对BG赋初值时可忽略不计.

1.4. 实时定轨算法

实时定轨常以序贯或者递推方法进行数据处理和参数估计,递推方法常采用EKF滤波算法. 基于星载BDS的低轨卫星定轨模型和观测模型,与GPS没有实质差异,区别在于观测模型中的时间系统、坐标系统和观测量精度. EKF的状态量包括卫星位置rsf、速度vsf和接收机钟差δtr). 对于单频定轨,EKF状态量还包括GRAPHIC组合的模糊度BG. EKF算法流程包括时间更新和测量更新,其时间更新方程为

${{Y}}_i^ - = f\left( {{t_i};{{Y}}_{i - 1}^ + } \right),$

${{P}}_i^ - = {{{\varPhi }}_i}{{P}}_{i - 1}^ - {{\varPhi }}_i^{\rm{T}} + {{{Q}}_i}.$

式中:下标ii−1为相应的观测历元;ti为观测历元i的时间; ${{Y}}_{i - 1}^ + $为EKF状态量的后验估计,即i−1历元的最终计算结果; ${{Y}}_{i }^ - $为EKF状态量的先验估计;f为EKF状态量的递推函数,与时间ti和上一历元的EKF状态量 ${{Y}}_{i - 1}^ + $有关; ${{P}}_{i - 1}^ - $为EKF状态量协方差的先验估计;Фi为EKF状态转移矩阵,也是f函数的雅各比矩阵;ФiTФi的转置;Qi:为EKF状态量的过程噪声. 式(8)为EKF状态量的递推公式,当前历元接收机钟差值与上一历元的值相等,模糊度的递推使用模糊度更新算法实现. EKF的测量更新方程为

${{{K}}_i} = {{P}}_i^ - {{H}}_i^{\rm{T}}{\left( {{{{W}}_i} + {{{H}}_i}{{P}}_i^ - {{H}}_i^{\rm{T}}} \right)^{ - 1}},$

${{Y}}_i^ + = {{Y}}_i^ - + {{{K}}_i}\left[ {{{{z}}_i} - h\left( {{{Y}}_i^ - } \right)} \right],$

${{P}}_i^ + = \left( {{{I}} - {{{K}}_i}{{{H}}_i}} \right){{P}}_i^ - .$

式中:i为相应的观测历元;Ki为卡尔曼增益矩阵;zi为观测量;h为观测量关于EKF状态量的函数;Hih的雅克比矩阵;Wi为观测噪声的协方差矩阵;Yi+为EKF状态量的后验估计,即i历元位置、速度、接收机钟差和模糊度的值;I为单位矩阵;Pi+为EKF状态量的后验协方差矩阵. 对于单频数据,zi为GRAPHIC组合,即 ${{{z}}_i} = \left[ {\rho _{\rm{G}}^1,\rho _{\rm{G}}^2, \cdots ,\rho _{\rm{G}}^{{n}}} \right]$,对于双频数据,zi为伪距消电离层组合 ${{{z}}_i} = $ $ \left[ {\rho _{{\rm{IF}}}^1,\rho _{{\rm{IF}}}^2, \cdots ,\rho _{{\rm{IF}}}^{{n}}} \right]$.

2. 原始观测数据实测及精度评估

浙江大学皮星三号(ZDPS-3)是浙江大学研制的新型高功能密度应用型皮纳卫星,主要用于对地面测控系统进行标定和校准[21-22]. 星载GNSS导航接收机为该卫星的重要载荷,由浙江大学自主研制,除了支持GPS L1和L2频点,还支持北斗B1和B3双频模式,其样机如图1所示,支持的伪码信号如表1所示. 其中,f 为载波频率,fc 为码速率.

图 1

图 1   搭载于ZDPS-3的星载GNSS接收机样机实物图及原理框图

Fig.1   Photo and block diagram of GNSS receiver prototype for ZDPS-3


表 1   ZDPS-3接收机支持的频点

Tab.1  Supported frequencies of the receiver for ZDPS-3

信号 f/MHz fc/ Mcps 信号 f/MHz fc/ Mcps
B1I 1561.098 2.046 L1 C/A 1575.42 1.023
B3I 1268.52 10.23 L2C 1227.60 1.023

新窗口打开| 下载CSV


零基线法是工程上常用的评估接收机原始观测数据性能的方法. 实验方法是通过特制的L波段信号分路器,将接收机天线收到的导航信号分成完全相等的两路,馈送给2台接收机. 理论上两路信号的对流层、电离层等主要测量误差来源相同,通过将观测到的伪距、载波相位数据在接收机之间和导航卫星之间进行两级差分,得到双差观测值,消除主要误差来源,因此双差观测值的期望值为0. 测量噪声为原始观测噪声的2倍,进而对接收机的原始观测数据质量作出评估.

为保证信号质量,实验地点选在楼顶空旷地带,时间为2018年7月1日22点至次日5点,接收机采用上文所述ZDPS-3 GPS/BDS双模接收机. 统计载噪比高于42 dBHz的卫星数据,以仰角最高的卫星为参考卫星,计算其他卫星相对于参考卫星的双差观测值. 1 h内BDS及GPS的外场实测伪距双差数据如图23所示,ρd为伪距的双差观测值. 统计观测时段所有卫星数据,B1伪距双差噪声标准差为0.43 m,B3伪距双差噪声标准差为0.17 m,而L1、L2分别为0.60与0.62 m. 如前所述,码速率的提高能够降低热噪声的影响,提升码环跟踪精度. 本文重点是利用B3伪码速率高的优势验证提升定轨精度的可能性,尚未对B3频点对应的环路参数进行优化设计. 根据实测结果,B3实际的伪距观测精度达到B1频点的2.53倍,且显著优于L1、L2频点.

图 2

图 2   BDS伪距双差观测值

Fig.2   Double difference pseudorange observation for BDS


图 3

图 3   GPS伪距双差观测值

Fig.3   Double difference pseudorange observation for GPS


统计载波相位观测数据,B1和B3载波相位双差值的标准差分别为3.00 mm和3.60 mm,B3频点的载波相位精度略差于B1,这是由于B3载波波长大于B1频点,从而造成绝对精度的差异;如果换算成载波周数,则两者精度较为接近. 同理,L1频点外场实测载波相位双差值标准差为2.80 mm,略优于L2频点的3.90 mm.

3. 半物理仿真验证

由于目前暂时没有搭载北斗B3频点接收机的在轨卫星数据,本文通过半物理仿真实验进行验证. BDS测试使用率先支持B3频点的卫导NSS8000模拟器,半物理仿真平台如图4所示. 如图5所示为仿真实验的流程图. 本文的BDS测试数据采用2018年2月6日的仿真观测数据,低轨卫星的轨道根数设置值如表2所示,其中,a 为轨道半长轴,e 为偏心率,i 为轨道倾角,Ω 为升交点经度,ω 为近日点幅角,M 为平近点角. 由于模拟器尚未支持北斗三号系统,基于星座卫星数为14颗的北斗二号系统进行仿真测试. 由于北斗二号主要覆盖亚太地区,当低轨卫星运行到覆盖区域外时,仅能观测到少数卫星而不足以定位,因此可稳定定轨时长占轨道周期的1/3. 为此,本文主要评估可定轨时段的性能. 为对照评估本文所用接收机上的GPS定轨性能,使用英国思博伦GSS9000模拟器在相同条件下对GPS定轨进行仿真,GPS导航星座的卫星数为32颗,仿真时长为低轨卫星的一个轨道周期(5 400 s).

表 2   低轨卫星轨道根数

Tab.2  Orbit parameters of low-orbit satellite

参数 设置值 参数 设置值
$a$ 6 907 395.5 m $\varOmega $ 115.5611°
$e$ 0.001163 $\omega $ 128.1038°
$i$ 97.4304° $M$ 76.7532°

新窗口打开| 下载CSV


图 4

图 4   半物理仿真实验平台

Fig.4   Hardware-in-loop simulation platform


图 5

图 5   星载GPS/BDS定轨仿真实验流程图

Fig.5   Flow chart for simulation of GPS/BDS orbit determination


3.1. 双频伪距组合

消电离层组合观测模型采用双频伪距组合,以模拟器设置的精密卫星轨道为参考,统计定轨结果的位置误差并表示在地心地固(earth centered earth fixed,ECEF)坐标系下. 如图6(a)所示为BDS定轨位置误差,如图6(b)所示为GPS定轨位置误差. 如表3所示为定轨位置误差均方差(root mean square,RMS)统计结果.

表 3   位置误差均方差统计表

Tab.3  Root mean square statistics of orbit determination

双频伪距组合 RMS/m
X Y Z 三轴
BDS 1.88 10.04 3.88 10.92
GPS 3.49 7.45 3.79 9.06

新窗口打开| 下载CSV


图 6

图 6   双频伪距组合实时定轨结果

Fig.6   Real-time orbit determination results of dual frequency pseudorange combination


由以上结果可见,GPS的总体定轨精度略优于BDS,这是由于GPS星座为完全运行状态,可用星数更多,仿真时段的平均可用星数达到9颗. 而BDS星座为部分运行状态(总共仅14颗星),平均可用星数仅为6颗. 统计仿真时段的几何精度因子(geometric dilution precision, GDOP)如图7所示,GPS的平均GDOP为1.78;而BDS的GDOP值波动较大,在亚太区域内的平均GDOP值为3.60,亚太区域外由于星数太少,几何分布急剧变差甚至导致不能定位. 如前所述,双频伪距组合会放大伪距噪声,定轨精度甚至有可能劣于单频伪距定轨精度. 表3为定轨位置误差均方差(root mean square, RMS)统计结果。其中,RMSX、RMSY、RMSZ、和RMSS分别为ECEF坐标系下X轴、Y轴、Z轴以及三轴方向的位置误差均方差. 如表3所示,无论是GPS还是BDS,采用双频伪距方式定轨的三轴位置精度接近10 m量级,不足以满足高精度实时定轨的需求.

图 7

图 7   GDOP统计

Fig.7   Statistics of GDOP


3.2. GRAPHIC组合

GRAPHIC组合方式引入了载波相位观测量,能够大幅降低观测噪声. 在相同测试条件下,采用GRAPHIC消电离层组合观测模型,统计仿真测试定轨结果的位置误差并表示在ECEF坐标系下. 如图8所示为BDS定轨位置误差,如图9所示为GPS定轨位置误差. 如表45所示分别为BDS、GPS为定轨位置误差均方差的统计结果.

表 4   BDS位置误差均方差统计表

Tab.4  Root mean square statistics of orbit determination for BDS

双频伪距组合 RMS/m
X Y Z 三轴
B1 0.72 3.77 1.45 4.10
B3 0.83 2.47 0.99 2.79

新窗口打开| 下载CSV


表 5   GPS位置误差均方差统计表

Tab.5  Root mean square statistics of orbit determination for GPS

双频伪距组合 RMS/m
X Y Z 三轴
L1 1.21 1.33 1.72 2.49
L2 1.34 1.47 1.73 2.64

新窗口打开| 下载CSV


图 8

图 8   BDS GRAPHIC组合实时定轨结果

Fig.8   Real-time orbit determination results of GRAPHIC for BDS


图 9

图 9   GPS GRAPHIC组合实时定轨结果

Fig.9   Real-time orbit determination results of GRAPHIC for GPS


由以上结果可见,B1频点的定轨位置精度达到4.10 m,较双频伪距组合有显著提升;B3频点由于具有码速率优势,精度提升至2.79 m. 如前所述,GPS由于目前在轨星数更多,可以达到更优的GDOP值,故L1、L2仿真结果稍优. 待北斗三代建设完成,根据前文观测数据实测评估对比结果,预期可以获得优于GPS的定轨性能. 总体而言,采用GRAPHIC组合可以获得明显优于双频伪码组合的定轨精度,而B3频点由于测距码抗噪声性能优势显著,可以在GRAPHIC组合下获得相对较好的定轨精度.

伪距测量噪声是GRAPHIC组合观测精度的主要误差源. 根据外场实测及在轨场景仿真实验结果,B3伪距噪声低的优势可以明显地体现在定轨结果中,在实际低轨卫星应用中,采用B3频点GRAPHIC组合进行低轨卫星实时定轨具有较好的工程应用价值. 需要指出的是,本文的主要贡献为基于B3频点进行低轨卫星实时定轨试验,并评估B3频点实时定轨性能及其相对于常用单频导航频点(GPS和BDS B1)的优势,因此并未将提升观测量的绝对精度和定轨算法本身作为研究重点. 通过接收机环路的优化设计可提高观测性能,同时采用简化动力学定轨代替运动学定轨等措施,可进一步提升实时定轨精度.

4. 结 语

本文在对星载北斗接收机实测数据进行评估的基础上,利用导航信号模拟器进行了半物理仿真实验,重点对基于星载北斗接收机B3频点的低轨卫星实时定轨性能进行了评估. 实验结果表明:采用B1/B3双频伪距组合方式的定轨精度较低,位置精度为10 m左右,而采用GRAPHIC消电离层观测组合能实现定轨精度的大幅提升. 特别是在仅使用北斗二号14颗导航卫星的条件下,基于B3频点GRAPHIC组合的实时定轨精度达到2.79 m,优于B1频点的4.10 m,并与处于完全运行状态的GPS系统的定轨精度相当,可满足高精度实时定轨的需求. 在北斗三代全球导航系统建设完成之后,低轨卫星的定轨性能将得到进一步提升,可有力保障我国的航天任务.

参考文献

MONTENBRUCK O, RAMOS-BOSCH P

Precision real-time navigation of LEO satellites using global positioning system measurements

[J]. GPS Solutions, 2008, 12 (3): 187- 198

DOI:10.1007/s10291-007-0080-x      [本文引用: 1]

BOCK H, HUGENTOBLER U, SPRINGER T A, et al

Efficient precise orbit determination of LEO satellites using GPS

[J]. Advances in Space Research, 2002, 30 (2): 295- 300

DOI:10.1016/S0273-1177(02)00298-3      [本文引用: 1]

LI M, LI W W, SHI C, et al

Precise orbit determination of the Fengyun-3C satellite using onboard GPS and BDS observations

[J]. Journal of Geodesy, 2017, 91 (11): 1313- 1327

DOI:10.1007/s00190-017-1027-9      [本文引用: 1]

SHI C, FAN L, LI M, et al

An enhanced algorithm to estimate BDS satellite's differential code biases

[J]. Journal of Geodesy, 2016, 90 (2): 161- 177

DOI:10.1007/s00190-015-0863-8     

XIONG C, LU C F, ZHU J, et al

Orbit determination using real tracking data from FY3C-GNOS

[J]. Advances in Space Research, 2017, 60 (3): 543- 556

DOI:10.1016/j.asr.2017.04.013     

YI B, GU D F, CHANG X, et al

Integrating BDS and GPS for precise relative orbit determination of LEO formation flying

[J]. Chinese Journal of Aeronautics, 2018, 31 (10): 2013- 2022

DOI:10.1016/j.cja.2018.06.001     

ZENG T, SUI L F, JIA X L, et al

Results and analysis of BDS precise orbit determination with the enhancement of Fengyun-3C

[J]. Acta Geodaetica Et Cartographica Sinica, 2017, 46 (7): 824- 833

CAI Y R, BAI W H, WANG X Y, et al

In-orbit performance of GNOS on-board FY3-C and the enhancements for FY3-D satellite

[J]. Advances in Space Research, 2017, 60 (12): 2812- 2821

DOI:10.1016/j.asr.2017.05.001      [本文引用: 1]

佘丽丽. 基于星载GNSS的LEO实时精密定轨算法研究[D]. 北京: 中国科学院大学, 2018: 56.

[本文引用: 1]

SHE Li-li. Research on LEO real-time precise orbit determination based on onboard GNSS [D]. Beijing: Chinese Academy of Sciences, 2018: 56.

[本文引用: 1]

熊超, 贾小林, 朱俊, 等

星载BDS/GPS低轨卫星定轨精度分析

[J]. 测绘科学, 2017, (6): 121- 128

[本文引用: 1]

XIONG Chao, JIA Xiao-lin, ZHU Jun, et al

Precision analysis of orbit determination with space-borne BDS and GPS measurements for LEO satellite

[J]. Science of Surveying and Mapping, 2017, (6): 121- 128

[本文引用: 1]

谢钢. GPS原理与接收机设计[M]. 北京: 电子工业出版社, 2009.

[本文引用: 1]

BETZ J W, KOLODZIEJSKI K R

Generalized theory of code tracking with an early-late discriminator part I: lower bound and coherent processing

[J]. IEEE Transactions on Aerospace and Electronic Systems, 2009, 45 (4): 1538- 1556

DOI:10.1109/TAES.2009.5310316      [本文引用: 1]

WANG N, YUAN Y B, LI Z, et al

Improvement of Klobuchar model for GNSS single-frequency ionospheric delay corrections

[J]. Advances in Space Research, 2016, 57 (7): 1555- 1569

[本文引用: 1]

王军. 星载GPS单频实时定轨电离层改正方法研究[D]. 武汉: 武汉大学, 2017: 21.

[本文引用: 1]

WANG Jun. Novel methods for improving ionospheric delay correction in orbit determination with single-frequence space-borne GPS observation [D]. Wuhan : Wuhan University, 2017: 21.

[本文引用: 1]

MONTENBRUCK O

Kinematic GPS positioning of LEO satellites using ionosphere-free single frequency measurements

[J]. Aerospace Science and Technology, 2003, 7 (5): 396- 405

[本文引用: 1]

阮仁桂, 吴显兵, 冯来平

单频精密单点定位观测模型和电离层处理方法比较

[J]. 武汉大学学报: 信息科学版, 2013, 38 (9): 1023- 1028

[本文引用: 1]

RUAN Ren-gui, WU Xian-bing, FENG Lai-ping

Comparison of observation models and ionospheric elimination approaches for single frequency precise point positioning

[J]. Geomatics and Information Science of Wuhan University, 2013, 38 (9): 1023- 1028

[本文引用: 1]

王瞧, 金小军, 阳芳, 等

浙江大学皮星二号实时精密自主定轨算法与半实物仿真

[J]. 传感技术学报, 2016, 29 (8): 1193- 1199

DOI:10.3969/j.issn.1004-1699.2016.08.012      [本文引用: 1]

WANG Qiao, JIN Xiao-jun, YANG Fang, et al

Precision real-time orbit determination algorithm and hardware-in-the-loop simulation for ZDPS-2 satellite

[J]. Chinese Journal of Sensors and Actuators, 2016, 29 (8): 1193- 1199

DOI:10.3969/j.issn.1004-1699.2016.08.012      [本文引用: 1]

DAI X L, LOU Y D, DAI Z Q, et al

Real-time precise orbit determination for BDS satellites using the square root information filter

[J]. GPS Solutions, 2019, 23 (2): 1- 14

GUO R, CHEN J P, ZHU L F, et al

Kinematic orbit determination method optimization and test analysis for BDS satellites with short-arc tracking data

[J]. Acta Geodaetica Et Cartographica Sinica, 2017, 46 (4): 411- 420

CHEN P, SHU L Z, DING R, et al

Kinematic single-frequency relative positioning for LEO formation flying mission

[J]. GPS Solutions, 2015, 19 (4): 525- 535

DOI:10.1007/s10291-014-0410-8      [本文引用: 1]

周阳, 王春晖, 金小军, 等

新型皮卫星星载接收机设计

[J]. 传感技术学报, 2011, 24 (5): 694- 699

DOI:10.3969/j.issn.1004-1699.2011.05.014      [本文引用: 1]

ZHOU Yang, WANG Chun-hui, JIN Xiao-jun

A new receiver structure of TT&C transponder for pico-satellite

[J]. Chinese Journal of Sensors and Actuators, 2011, 24 (5): 694- 699

DOI:10.3969/j.issn.1004-1699.2011.05.014      [本文引用: 1]

DEND S C, MENG T, WANG H, et al

Flexible attitude control design and on-orbit performance of the ZDPS-2 satellite

[J]. Acta Astronautica, 2017, 130 (1): 147- 161

[本文引用: 1]

/