基于ORBCOMM 卫星机会信号的定位技术
2021-01-08秦红磊谭滋中丛丽赵超
秦红磊,谭滋中,丛丽,赵超
(北京航空航天大学 电子信息工程学院,北京100083)
尽管全球导航卫星系统(Global Navigation Satellite System,GNSS)随着位置服务需求的日益增长得到了快速的发展,但其存在信号弱、容易受到干扰等缺点,越来越多的备用导航定位技术被提出,以摆脱对GNSS系统的依赖。利用广泛存在的无线电机会信号导航(Opportunistic Navigation,Op Nav)[1-2]及协作机会信号导航(Collaborative Opportunistic Navigation,COpNav)近年来得到了广泛的关注[3]。典型系统如美国提出的全源导航(All Source Positioning and Navigation,ASPN)技术,该技术提出利用一切可以利用的信号资源实现融合导航定位;类似的还有英国国防部研发的Navsop(Navigation via Signals of Opportunity)系统。
机会信号功率强、占用频谱宽、广泛存在,因此具有抗干扰能力强且投入低的优点[4],而将低轨卫星无线电信号作为机会信号具有更广的覆盖范围。目前,成熟的典型低轨卫星系统有铱星(IRIDIUM)、轨道通信卫星(ORBCOMM)和全球星(Globalstar)等[5-7]。Globalstar卫星系统和ORBCOMM卫星系统本身均具备位置信息服务,铱星公司2016年宣布在第二代卫星 IRIDIUM NEXT系统中增加了定位服务,然而以上低轨卫星定位服务需授权才可使用。
国外计划未来建造大量低轨卫星系统,2015年OneWeb公司宣布在未来构建一个含有648颗卫星的低轨卫星系统[8];同年SpaceX和Google宣布合作构建一个含有更多卫星的低轨卫星系统,共含有4 000颗卫星[9];三星公司也计划在未来建立含有4 600颗低轨卫星的卫星系统[10];Boeing公司也启动了构建低轨卫星系统计划[11]。国内也计划发射大量的低轨卫星,如鸿雁全球卫星星座、虹云工程计划。现有低轨卫星系统和未来计划发射的大量低轨卫星为基于低轨卫星机会信号的定位技术提供了基础。
本文研究将低轨卫星作为机会信号源,在不需要授权下即可利用低轨卫星实现定位。在前期研究中,利用铱星(IRIDIUM NEXT)机会信号实现了定位[12]。本文主要研究将另一种低轨卫星(ORBCOMM)信号作为机会信号实现定位。由于ORBCOMM 卫星机会信号仅部分信息已知,重点分析信号体制并给出定位观测信息获取方法。提出利用瞬时多普勒定位方法实现定位,给出了详细的定位数学模型,对多普勒定位几何精度因子的物理意义进行了解释。通过搭建卫星信号采集平台,实现了实际卫星信号的定位。首先,对ORBCOMM 卫星机会信号通信体制进行了深入分析,为定位信息观测量的获取提供了基础。其次,介绍了基于ORBCOMM 卫星机会信号的定位原理,给出了定位观测信息多普勒频率的估计方法,以及瞬时多普勒定位原理及数学模型,并阐述了多普勒定位几何精度因子的物理意义。最后,利用实际接收的ORBCOMM第二代卫星(OG2)信号实现了定位,验证了本文方法的有效性和正确性。
1 ORBCOMM 卫星轨道及信号体制
ORBCOMM卫星星座分布决定了接收机可见卫星数量,而定位观测信息提取方法与信号体制有关,本节主要分析ORBCOMM 卫星星座结构和信号体制特性。
ORBCOMM卫星通信系统是能够实现数据全球通信的卫星系统[13]。ORBCOMM 空间星座最初计划如图1所示。主星座4个轨道平面(A,B,C,D),每个平面均匀配置8颗星,轨道高度为825 km,倾角为45°;辅助星座2个轨道平面(G,F),每个平面配置2颗星,倾角分别为108°和70°,轨道高度为775 km;赤道面(E)曾设6颗高度为975 km的卫星[14]。ORBCOMM 公司于2008年宣布部署第二代卫星(OG2)星座。目前,共12颗OG2卫星在轨运行,均匀分布在4个主轨道平面。OG2卫星轨道高度为620 km,轨道倾角为47°,轨道周期为97 min。本文采用OG2卫星实现接收机定位。
ORBCOMM卫星下行链路占用频带为137~138 MHz,共包括13个信道,频带分配如图2所示。12个带宽为25.0 kHz的信道用于与用户终端通信(采用FDMA多址连接),另外1个带宽为50.0 kHz的信道用于与关口站通信。所有ORBCOMM卫星共用12个用户链路下行信道,频率复用次数为4次[15]。本文主要采用卫星至用户下行链路信号实现定位。
图1 ORBCOMM系统星座计划Fig.1 ORBCOMM system constellation plan
图2 ORBCOMM卫星下行链路频带分配Fig.2 ORBCOMM down link frequency band assignment
用户下行链路采用SDPSK(Symmetrical Differential Phase Shift Keying)调制方式,数据率为4 800 bit/s。调制采用0数据状态和1数据状态分别引起相位负90°和正90°变化。
采集的实际信号显示每颗OG2卫星最多同时在2个信道内发射信号,所占信道频带在几天内保持不变。ORBCOMM 卫星在用户下行链路内持续发射SDPSK调制的信号,为提取定位观测信息并最终实现接收机定位提供了基础。
2 ORBCOMM 定位方法
接收机定位需要观测信息、卫星轨道和定位解算,本文利用ORBCOMM OG2卫星的多普勒频率作为定位观测信息,采用公开的TLE数据及轨道预测模型[16]计算OG2卫星轨道,定位方法采用瞬时多普勒定位技术。本节主要给出了OG2卫星多普勒频率的估计方法,并介绍了瞬时多普勒定位原理及数学模型。
2.1 ORBCOMM 卫星机会信号多普勒频率估计
OG2卫星机会信号到达地面功率一般高于噪声,然而相比于直接通过快速傅里叶变换(Fast Fourier Transform,FFT)处理并根据频谱对称性或峰值谱线对应的频率获得多普勒频率粗测量的其他信号,SDPSK调制的OG2信号无法通过其频谱直接获得多普勒测量值。另外,采用最大似然估计(Maximum Likelihood Estimation,MLE)方法[17]进行多普勒精确测量时,根据多普勒粗测量值产生的本地载波与OG2信号进行相关时,相关峰会受到数据位的影响。综上,本文先对OG2卫星机会信号平方处理获得2倍频信号,再利用FFT和MLE方法估计多普勒频率。
SDPSK调制信号表达式为
可以看出,平方后的信号码元相位变化为±π,是一种基带信号为0/1交替的特殊BPSK调制信号。基带信号不断重复0/1交替过程中,信号频谱分别在中心频率的左右两侧对称出现峰值谱线,其间距为码速率大小。令平方后的信号频谱中的2个峰值谱线对应频率分别为f1和f2且f1<f2,则原始信号中心频率f表达式为
式中:I(f)为最优估计函数;f为载波频率估计;x[n]为平方信号。
2.2 ORBCOMM 卫星多普勒定位方法
接收机获得ORBCOMM 卫星的多普勒频率测量后,可采用瞬时多普勒定位方法实现定位解算。本节主要介绍利用ORBCOMM 卫星机会信号实现接收机定位解算的原理及数学模型。
接收机在某时刻接收到ORBCOMM 卫星机会信号并测量出多普勒频率,其位置一定位于以ORBCOMM卫星为顶点的圆锥面上,且该圆锥面上所有点对应的多普勒频率均相等。当接收机同时得到多个多普勒频率时,通过多个等多普勒圆锥面的交点即可获得接收机位置[18],此即为ORBCOMM卫星机会信号多普勒定位原理的几何解释。一般情况下,同时可见ORBCOMM 卫星个数至多2颗,对于位置保持不变的静止接收机,可采用多历元定位方法[19],即采用多个不同时刻的ORBCOMM卫星多普勒频率实现接收机定位解算。
下面给出ORBCOMM 卫星多普勒定位的数学模型。静止接收机在某一时刻测量得到的卫星多普勒频率表达式为
式中:˙ρ为ORBCOMM卫星的多普勒频率,m/s;ru为接收机位置向量;r和v分别为ORBCOMM 卫星位置矢量和速度矢量;fu为接收机钟偏;ε˙ρk为多普勒频率误差。
式中:rk和vk分别为第k个多普勒频率测量对应的ORBCOMM 卫星的位置矢量和速度矢量;E3×3为单位矩阵。
根据式(8)~式(10)得到ORBCOMM卫星多普勒定位线性导航状态更新方程为
式中:Δu=[ΔruΔfu]为接收机状态更新向量;G为多普勒定位雅可比矩阵;ε为测量和线性化误差向量。
利用牛顿迭代运算,令ui和Δui分别表示第i次接收机状态量估计和修正量,则第i+1次迭代解算的接收机状态量估计ui+1的表达式为
多普勒频率测量等效为速度时,DOPD的单位为s,物理意义代表接收机真实位置以ORBCOMM 卫星多普勒测量误差的等效速度到达定位解算的位置所用的时间。
接收机持续获得ORBCOMM 卫星多普勒频率测量并利用TLE数据及轨道预测模型获得卫星轨道后,即可利用上述多普勒定位方法实现基于ORBCOMM卫星机会信号的接收机定位。
3 实验验证
3.1 ORBCOMM 实际信号分析及多普勒频率估计
本节主要对实际采集的ORBCOMM 卫星(OG2)机会信号进行时频域分析,并估计信号多普勒频率。采集平台主要包括天线、射频前端和数据处理单元,采集系统框图如图3所示。
图4给出了OG2卫星机会信号的实时频谱。此时仅有一颗OG2卫星划过接收机头顶,该颗卫星同时在信道11和信道9发送信号,分别对应图中A信号和B信号。2个信道的标准频率分别为137.287 5 MHz和137.737 5 MHz。在137.5 MHz频带附近的C信号为硬件采集平台射频前端引起的噪声信号。
图5给出了OG2卫星机会信号的三维瀑布图。可以看出,在9 min内可持续接收到OG2卫星信道9和信道11的信号,信号功率明显高于噪声。
图3 ORBCOMM卫星机会信号采集系统框图Fig.3 Block diagram of ORBCOMM signals of opportunity collection system
图4 ORBCOMM卫星用户链路信道9和信道11信号功率谱Fig.4 Signal power spectrum of ORBCOMM user links channel 9 and channel 11
图5 ORBCOMM卫星用户链路信道9和信道11信号瀑布图Fig.5 Signal waterfall plot of ORBCOMM user links channel 9 and channel 11
上述给出了OG2卫星机会信号实时分析结果,接下来将采集器中心频率设置为138 MHz,持续采集数据5 min。OG2卫星信道9和信道11信号频率经下变频变为27.7375 MHz和27.2875 MHz。对采集的信号进行滤波处理得到相应信道信号数据,并估计信号多普勒频率。
取100 ms的OG2卫星信道11信号进行平方处理并进行FFT变换后,可以获得2个峰值谱线,如图6所示。2个峰值谱线对应的频率f1和f2分别为54 577234.8 Hz和54582034.8 Hz,2个谱线的距离为4 800 Hz,其值等于OG2卫星机会信号的数据率。根据式(4)可以获得原始信号的中心频率粗略估计结果为27 289 817.4 Hz,对应的多普勒频率粗测量值为2 317.4 Hz。采用式(5)得到的信号中心频率精测量值为27 289 816.1 Hz,则卫星信号数据的多普勒频率精测量值为2 316.1 Hz。对OG2卫星信道9信号采用上述方法同样可以获得信号多普勒频率。
实验结果证明,OG2卫星持续在用户下行链路发送信号,并且可采用本文方法估计信号多普勒频率。
图6 ORBCOMM卫星信道11信号平方处理后的频谱Fig.6 Spectrum of ORBCOMM channel 11 signal after square processing
3.2 ORBCOMM 多普勒频率定位
本节对OG2卫星机会信号进行长时间采集,测量卫星多普勒频率,并结合卫星TLE轨道信息实现接收机静态定位。
持续接收OG2卫星机会信号48 min,期间接收机保持静止。将数据分为连续的多个360 ms数据块,采用匹配滤波器对每个数据块进行检测,当存在卫星机会信号时,估计信号中心频率和多普勒频率。图7给出此段时间内可见卫星的信号中心频率估计结果。
可以看出,48 min内可见卫星共3颗:Sat1、Sat2和Sat3,每颗卫星同时在2个信道内发射信号,共6个信号中心频率曲线,同一卫星不同信道信号的中心频率曲线形状大体相同。此段时间内任意时刻最多可见卫星1颗,Sat1和Sat3卫星星下点轨迹距离接收机较近且属于同一轨道,在信道9和信道11发送信号,Sat2卫星在信道6和信道10发送信号。
Sat1和Sat3卫星的信道11信号及Sat2卫星的信道6的信号多普勒曲线如图8所示。可以看到,ORBCOMM 卫星的多普勒范围在±3 kHz之间。Sat1和Sat3卫星的信道9及Sat2卫星的信道10的信号多普勒曲线与图8中对应卫星的多普勒曲线的形状大致相同。
本文利用卫星多个历元的多普勒频率实现定位,分为2颗卫星和3颗卫星定位2种情况。2颗卫星定位时,观测信息采用Sat1卫星信道11信号和Sat2卫星信道6信号的多普勒频率。3颗卫星定位时,增加了Sat3卫星信道11信号的多普勒频率观测。
图7 ORBCOMM卫星机会信号中心频率曲线Fig.7 Center frequency curves of ORBCOMM signals of opportunity
图8 ORBCOMM卫星用户链路信道6和信道11信号多普勒曲线Fig.8 Doppler curves of ORBCOMM userlinks channel 6 and channel 11 signal
首先,将每颗卫星所有不同时刻的多普勒频率测量值均匀分组;其次,在每组中随机选取一个测量值,不同卫星的共25个观测历元的多普勒频率组合成定位观测信息,并利用高程辅助[20]瞬时多普勒定位进行定位解算;最后,利用25个不同观测历元多普勒频率进行组合800次,分别进行定位解算并统计定位误差均值和RMS值。下面分别给出2颗卫星和3颗卫星情况下定位误差的统计结果和分析。
表1给出了高程辅助情况下的定位误差统计结果,此时高程辅助信息固定为高度真值。2颗卫星定位时,定位误差主要分布在北向,误差均值和RMS值均大于100 m,在东向上误差具有较大的波动。当增加卫星个数时,东向误差均值和RMS值变化较小,但北向误差明显有所改善,由于增加了卫星,改善了卫星几何分布,使得东向和北向误差RMS值都在80 m左右。
表2给出了无高程辅助和有高程辅助情况下的定位误差对比结果。可以看出,高程辅助可以有效改善定位精度,2颗卫星情况下,三维定位(3D定位)误差RMS值可以达到224 m,水平定位(2D定位)误差RMS值可以达到174 m,3颗卫星情况下,三维定位(3D定位)误差RMS值可以达到200 m,水平定位(2D定位)误差RMS值可以达到120 m。
表1 高程辅助多普勒定位误差统计结果Table 1 Error statistic r esults of Doppler positioning with height aid
表2 有无高程辅助情况下的定位误差结果对比Table 2 Comparison of error results of positioning with and without height aid
将高程辅助信息误差固定在40 m范围内,重复上述定位实验。实验结果显示,此时高程辅助信息误差对2颗卫星和3颗卫星定位的水平定位精度的影响均不超过15 m。综上所述,利用实际ORBCOMM 卫星机会信号可以实现接收机定位,高程辅助误差在40 m范围内时,接收机水平定位精度可以达到140 m。
4 结束语
本文提出利用ORBCOMM 卫星机会信号实现接收机定位,并利用实际信号进行实验验证。
对SDPSK调制的ORBCOMM 卫星机会信号进行平方处理,并利用FFT和MLE方法有效测量多普勒频率。利用瞬时多普勒定位方法,并结合轨道预测模型计算的卫星轨道信息可以实现接收机定位。实验结果表明,在高程辅助情况下水平定位精度可以达到140 m,验证了本文方法的正确性和有效性。为GNSS受到强干扰导致不能正常工作的环境下,提供一种有效的定位手段。
针对卫星机会信号定位技术,下一步主要研究联合采用ORBCOMM 和IRIDIUM 两种卫星机会信号实现接收机定位。