戴凯阳 邵 搏 熊 帅 张 键 雷哲哲
1 中国电子科技集团公司第二十研究所,西安市白沙路1号,710068
北斗系统在不断发展与完善的同时也需要与国际接轨,得到国际社会的认可。目前国际民航组织正在开展双频多星座(dual frequency and multi constellation,DFMC)星基增强系统(space based augmentation system,SBAS)标准与建议措施(standard and recommended practices,SARPs)的制定工作,BDS作为被增强对象、北斗星基增强系统(Beidou navigation satellite system,BDSBAS)作为服务提供商均被写入DFMC SBAS SARPs中[1-3]。在DFMC SBAS SARPs修订过程中需要对其中的相关内容进行验证,其中北斗系统新频点多路径误差特性验证工作是一个意义重大的议题。美国、欧洲已经分别完成了GPS和Galileo系统的多路径误差特性验证工作[4-6],北斗系统急需开展相关研究并形成工作文件。因此,开展北斗系统多路径误差特性的验证工作不仅可以弥补国内在该方面研究的不足,而且可以助力北斗系统写入国际民航组织标准,为推动北斗系统走向国际迈出坚实一步。
国际民航组织针对导航卫星系统的多路径误差特性制定了相关规定,BDS多路径误差特性需要被JWGs/6提出的DFMC多路径误差模型[7]所包络。另外,导航卫星接收机数据所使用的接收机天线群延时需要满足航空无线电技术委员会发布的DO-373标准[8]。因此,本文为了验证BDS多路径误差特性符合DFMC多路径误差模型的要求,首先利用符合DO-373天线群延时特性的天线接收BDS在B1c和B2a上的观测数据;然后提出一种北斗三号新频点的多路径误差特性验证方法,对观测数据进行处理并绘制多路径误差随高度角的变化曲线;最后分析其与国际民航多路径误差模型的符合性。
本次实验测试的接收天线为NovAtel GNSS-750扼流圈天线,微波暗室实验设备主要包括矢量网络分析仪(Agilent E5071B)、标准天线(CR125B)和直流耦合器(Mini-Circuits ZFBT-4R2G-FT+)。实验设备连接示意图如图1所示,标准天线与被测天线处于同一水平方向上。
图1 微波实验设备连接示意图Fig.1 Test equipment connection diagram
为了得到NovAtel GNSS-750扼流圈天线在北斗三号新频点(B1c和B2a)上的群延时参数,进行天线群延时实验。测试步骤如下:
1)搭建测试环境,将NovAtel GNSS-750扼流圈天线通过直流耦合器与矢量网络分析仪端口2连接,将标准天线与矢量网络分析仪端口1连接,并使2个天线正面法向对齐。
2)设置矢量网络分析仪测试频点,在B1c频点1 575.42 MHz(或B2a频点1 176.45 MHz)校准矢量网络分析仪。
3)启动直流电源,设置输出电压值为5 V。
4)将天线的高度角调整为0°,控制转台水平方向转动1周,同时记录矢量网络分析仪在指定方位角Azi(0°、10°、20°、30°…、330°、340°、350°)的天线群延时值τ(Azi,Ele),其中Azi为方位角,Ele为高度角。
5)控制转台转动以使被测天线在高度角方向转动5°,重复步骤4),逐步增大高度角数值,直至完成垂直方向90°的测试,同时记录指定高度角Ele(5°、10°、15°、20°、…、80°、85°、90°)的数据。
(1)
图2给出NovAtel GNSS-750扼流圈天线在B1c和B2a频点的天线群延时方位图,可以看出,各频点的天线群延时主要跟高度角相关。
图2 B1c和B2a频点天线群延时方位图Fig.2 Antenna group delay on B1c and B2a
图3给出在B1c和B2a频点的最大天线群延时随高度角的变化曲线,图中虚线为DO-373天线群延时要求包络线,实线为实测各高度角对应的最大天线群延时曲线。可以看出,NovAtel GNSS-750扼流圈天线各高度角最大天线群延时被要求曲线包络,证明其在B1c和B2a频点上的天线群延时特性符合DO-373的要求。
图3 B1c和B2a频点的最大天线群延时随高度角变化曲线Fig.3 Maximum antenna group delay varies with elevation angles on B1c and B2a
利用NovAtel GNSS-750扼流圈天线采集原始观测数据,在西安某单位楼顶架设接收机天线,接收机型号为Septentrio PolaRx5。采集2020-09-01~08的观测数据,包括BDS B1c/B2a频点以及GPS L1/L5频点。
基于采集的BDS B1c/B2a和GPS L1/L5频点上的伪距观测量和载波相位观测量,给出本文验证方法的具体实施步骤(图4)。
图4 多路径误差特性验证流程Fig.4 Verification flow chart of mutipath error characteristics
北斗三号监测站采集到北斗三号卫星的观测数据和导航电文信息,其新频点B1c、B2a基本观测方程如下[9]:
ρB1c=r+c(dtu-dts)+T+
IB1c+MPρ,B1c+ερ,B1c
(2)
φB1c=r+c(dtu-dts)+T-IB1c+
NB1cλB1c+MPφ,B1c+εφ,B1c
(3)
ρB2a=r+c(dtu-dts)+T+
IB2a+MPρ,B2a+ερ,B2a
(4)
φB2a=r+c(dtu-dts)+T-IB2a+
NB2aλB2a+MPφ,B2a+εφ,B2a
(5)
式中,ρB1c和φB1c分别为B1c频点伪距和载波相位观测值,c为光速,dtu和dts分别为接收机钟差和卫星钟差,T为对流层延迟,IB1c为B1c频点电离层延迟,MPρ,B1c和MPφ,B1c分别为B1c频点伪距和载波相位上的多路径误差,ερ,B1c和εφ,B1c分别为B1c频点伪距和载波相位上与接收机相关的噪声误差,NB1c为B1c频点载波相位整周模糊度,λB1c为B1c频点波长,ρB2a和φB2a分别为B2a频点伪距和载波相位观测值,IB2a为B2a频点电离层延迟,MPρ,B2a和MPφ,B2a分别为B2a频点伪距和载波相位上的多路径误差,ερ,B2a和εφ,B2a分别为B2a频点伪距和载波相位上与接收机相关的噪声误差,NB2a为B2a频点载波相位整周模糊度,λB2a为B2a频点波长。
1)在φB1c和φB2a上进行周跳探测和标记,并剔除信噪比较低的数据。周跳探测主要使用无几何距离组合(geometry-free,GF)周跳探测方法,检测量b为:
b=φB1cλB1c-φB2aλB2a
(6)
对检测量前后历元作差,此处差值阈值设置为0.05 m,当差值大于0.05 m时,则认为当前历元发生周跳,对观测数据进行标记。
剔除信噪比低的历元观测值,B1c频点信噪比门限为29 dBHz,B2a频点上信噪比门限为27 dBHz。
2)利用B1c和B2a频点上的伪距观测值和载波相位观测值计算伪距与载波的差值(code minus carrier,CMC)。
针对单频B1c或B2a频点,利用双频载波观测值消除单频CMC数据中的电离层延迟:
CMCB1c=ρB1c-φB1c=2IB1c+MPρ,B1c+
ερ,B1c-MPφ,B1c-εφ,B1c-NB1cλB1c
(7)
CMCB2a=ρB2a-φB2a=2IB2a+MPρ,B2a+
ερ,B2a-MPφ,B2a-εφ,B2a-NB2aλB2a
(8)
(9)
(10)
CMCI_free,B1c=CMCB1c-2IB1c=
MPρ,B1c+ερ,B1c-NB1cλB1c
(11)
CMCI_free,B2a=CMCB2a-2IB2a=
MPρ,B2a+ερ,B2a-NB2aλB2a
(12)
式中,CMCB1c和CMCB2a分别为B1c和B2a频点上的CMC,CMCI_free,B1c和CMCI_free,B2a分别为B1c和B2a频点消除单频电离层影响后的CMC值。因为伪距上的多路径误差与噪声远大于载波上的,因此式(7)和式(8)中-MPφ,B1c-εφ,B1c≈0和-MPφ,B2a-εφ,B2a≈0,式(11)和式(12)中-NB1cλB1c和-NB2aλB2a为常数。
双频组合计算方法为:
(13)
(14)
CMCI_free=ρI_free-φI_free
(15)
式中,ρI_free和φI_free分别为无电离层组合的伪距和载波相位观测值,fB1c和fB2a分别为B1c和B2a频点的频率,CMCI_free为无电离层组合的CMC。
3)利用CMC数据的平均值来消除整周模糊度的影响,即消除步骤2)中的常数项。遇到步骤1)中探测出现的周跳标记,则重启该过程:
(16)
(17)
按照上述方法对2020-09-01~09-08的BDS和GPS双频观测数据进行处理,得到B1c/L1、B2a/L5以及双频的多路径误差σMP&AGDV随高度角变化曲线,如图5~7所示。
图5 BDS B1c和GPS L1频点σMP&AGDV随高度角变化曲线Fig.5 The curves of σMP&AGDV varies with elevation angles on BDS B1c and GPS L1
图6 BDS B2a和GPS L5频点σMP&AGDV随高度角变化曲线Fig.6 The curves of σMP&AGDV varies with elevation angles on BDS B2a and GPS L5
图7 BDS和GPS双频组合频点σMP&AGDV随高度角变化曲线Fig.7 The curves of σMP&AGDV varies with elevation angles on BDS and GPS dual-requency combination
从图5~7可以看出:
1)基于2 d和8 d数据得到的BDS/GPS多路径误差曲线的变化趋势几乎保持一致,即多路径误差曲线的变化趋势基本上不受天数影响。
2)BDS与GPS分别在B1c/L1、B2a/L5和双频组合情况下的多路径误差曲线的变化趋势基本一致,即BDS与GPS多路径误差特性基本一致。
3)BDS在B1c、B2a以及双频组合情况下的多路径误差特性符合DFMC多路径误差模型的要求。
本文首先基于微波暗室实验分析了NovAtel GNSS-750天线的天线群延时特性,验证该天线的群延时特性符合DO-373的要求;然后提出一套可用于分析验证北斗三号新频点多路径误差特性的方法,并给出明确的处理流程和实施步骤;最后基于此天线,进行BDS和GPS原始双频观测数据的采集和处理。结果表明,BDS的多路径误差特性与GPS基本一致,并且能够满足国际民航提出的DFMC多路径误差模型要求。本文能够为北斗系统多路径误差特性分析提供理论依据和实施思路,并且为北斗系统写入国际民航组织标准增砖添瓦。
致谢:感谢欧洲DLR课题组提供参考资料。