BDS接收机天线相位中心标定.docx

BDS接收机天线相位中心标定.docx

《BDS接收机天线相位中心标定.docx》由会员分享，可在线阅读，更多相关《BDS接收机天线相位中心标定.docx（26页珍藏版）》请在冰豆网上搜索。

BDS接收机天线相位中心标定

摘要：

在高精度GNSS定位中，接收机天线相位中心偏差（PCO）和天线相位中心变化（PCV）的影响不可忽略。

目前，IGS发布的绝对天线相位模型文件中包含了GPS/GLONASS系统的标定值，但是没有发布北斗系统（BDS）的标定值。

本文借助机械臂可以控制天线自由旋转，在数小时内实现全方位GNSS观测的特性，采用历元间差分的方法对接收机天线包括GPSL1/L2和BDSB1I/B2I/B3I等多个频点的PCO和PCV分别进行标定和拟合。

标定结果表明，比较最小二乘估计的GPSPCO与IGS发布值，其STD和RMS在L1/L2上均小于1mm；BDSPCO估计值的STD在B1I/B2I/B3I上分别为0.5、0.3、0.3mm。

利用球谐函数拟合的GPSPCV格网值与IGS发布值相比，其偏差在天顶距小于75°时均小于1.5mm。

BDSPCV拟合值范围均在-5~8mm，且随天顶距变化曲线呈现波谷状。

BDSPCV在低高度角处拟合值波动较大，随方位角变化曲线峰值-峰值最大达到了5.6mm。

关键词：

GNSS    天线相位中心    天线相位中心偏差    天线相位中心变化    北斗系统    机械臂标定    

BDSreceiverantennaphasecentercalibration

Abstract:

Inhigh-preciseGNSSapplications,theeffectsofreceiverantennaphasecenteroffset（PCO）andantennaphasecentervariation（PCV）arenotnegligible.Atpresent,theantennaphasecorrectiongridfilereleasedbyIGSonlycontainsthecalibrationvaluesofGPS/GLONASS,buttherearenothecalibrationvaluesoftheBeiDousystem（BDS）released.Inthispaper,weusetime-differencetocalibrateandfitGPSL1/L2&BDSB1I/B2I/B3IPCOsandPCVsofunknownantenna,dependingonthecharacteristicsthatrobotcanbeusedtocontrolthefreerotationoftheantennaandrealizeomnidirectionalGNSSobservationinfewhours.ComparedwiththeIGSpublishedvalues,thecalibrationresultsshowthattheestimatedGPSPCOsusingtheleastsquaresonL1/L2arelessthan1mmforSTDandRMS.TheSTDoftheBDSPCOsestimationvaluesare0.5,0.3,0.3mmonB1I/B2I/B3I,respectively.TheGPSPCVgridvaluesfittedbythesphericalharmonicfunctionarelessthan1.5mmbelow75degreeszenithcomparedtotheIGSpublishedvalue.TheBDSPCVsfittingrangesfrom-5to8mmwithzenith,andthecurveistroughwiththezenithvariation.Asazimuthchanging,theBDSPCVshavealargefluctuationinfittingvaluesatlowelevation,andthepeaktopeakhasamaximumof5.6mm.

Keywords:

 GNSS    antennaphasecenter    PCO    PCV    BDS    robotcalibration    

随着GNSS接收机硬件技术和软件算法的不断提高，使用载波相位观测值获得厘米级甚至亚厘米级的定位精度成为现实[1]，并在形变监测、地壳运动监测及精密工程测量中得到了广泛应用。

GNSS载波相位观测值是以天线相位中心（antennaphasecenter,APC）为测量点获得的。

而在定位数据处理观测模型中的坐标是以天线参考点（antennareferencepoint,ARP）为基准[2]。

由于天线本身的特性,天线相位中心在水平方向、垂直方向上与ARP都存在偏差，即天线相位偏差，在垂直方向上最大可达分米量级，直接影响定位结果中的高程分量。

如果不顾及此项改正，基线高程方向误差最大可达分米量级[3-5]。

目前，IGS发布最新的IGS14.ATX相位中心模型文件中给出了大部分测地型天线GPS/GLONASS等系统的改正值，但尚未发布BDS改正值。

虽然不同的载波频率，其相位中心不同，但是BDS的B1/B2频率和GPSL1/L2频率较为接近，近些年针对BDS的天线相位中心改正通常采用GPS的公布值来近似代替。

随着我国北斗三代在2018年底完成基本系统建设[6-7]，BDS用户越来越广泛，有必要对BDS进行天线相位中心标定。

本文借助机械臂可以控制天线自由旋转、并在短时间内实现全方位GNSS观测的特性，采用历元间差分的方法对TRM59800天线包括GPSL1/L2和BDSB1I/B2I/B3I等多个频点的相位中心进行了标定。

1天线相位中心偏差（PCO）和天线相位中心变化（PCV）概述

天线相位中心（APC）分为两个部分：

天线相位中心偏差（PCO）和天线相位中心变化（PCV）[8]。

如图1所示，PCO定义为由ARP指向接收信号的平均相位中心在天线固定坐标系（antenna-fixedcoordinatesystem，AFCS）中的三维坐标矢量，用NEU（北东天）表示。

PCV定义为瞬时相位中心到PCO的距离标量，一般随入射高度角θ和方位角α变化[3]。

几乎所有高精度GPS天线的PCV受方位角影响很小，它主要受高度角变化影响。

图1 天线相位中心偏差和天线相位中心变化Fig.1 Antennaphasecenteroffsetandantennaphasecentervariation

图选项 

天线相位中心改正对观测值的影响Δ由PCO在信号方向e上的投影E与PCV之和组成[9]

（1）

式中，负号是因为根据IGS天线模型文件交换格式中的定义[10]，PCV对距离量的改正为正,其与PCO在卫星视线上的投影符号相反所致。

PCO测量方法有旋转天线法、交换天线法、完全流动观测法和相对参考天线的双差相位观测法等[11]。

这些方法需要观测时间较长，流程较为烦琐，而且只能求解PCO，无法解算随高度角和方位角的变化而变化的PCV。

PCV的测量方法相比于PCO要复杂很多，因为PCV不仅与信号频率有关，而且还随信号的强弱以及方向而变化。

因此精确标定PCV必须多次重复测量。

对于GNSS天线，目前有效的标定方法分别是微波暗室测量法、相对现场标定和机器人标定方法[12],其对比见表1。

微波暗室测量法设备昂贵，测量过程复杂，不宜推广；相对现场标定方法假定参考天线的PCV为零，PCO为已知值，因此会产生系统误差；传统机器人标定方法是在室外开阔无遮挡的环境中标定，在短基线的一端安置天线相位中心己知的天线，另一端在自动机器人上安装待测天线，第1天机器人通过按照一定的步长进行旋转测量采集数据。

根据GPS卫星星座在连续两个恒星日内相同时刻的空间分布相同的特性，第2天在同一恒星时将被测天线初始方向旋转某角度（如果接收机天线在两个恒星日内保持天线方向不变，则多路径和PCV均被消除）进行同样的测量。

然后利用这些观测值,解算出用于拟合PCV模型的球谐函数系数。

在IGS相位中心改正模型文件中，超过一半的接收机天线是用机器人校准得到的[9-10]。

但是，由于BDS系统GEO/IGSO/MEO卫星轨道周期不同，该方法对北斗系统并不适用，每次标定并不能充分利用每颗卫星的数据，且标定时间长。

针对上述问题，本文采用历元间差分的方法[13]，借助机械臂可以快速旋转方向并能在短时间内获取天线半球空间的GPS/BDS观测数据的特点，采集均匀覆盖半球的数千个方向观测值进行测量。

由于不同方向的观测值蕴含了天线旋转所引起的PCO和PCV变化信息，因此采用历元间差分模型既可以削弱多路径又能实现对PCO和PCV的同步标定。

表1 PCV标定方法比较Tab.1 ComparisonofPCVcalibrationmethods

类型

微波暗室测量

相对现场校准

机器人绝对标定

实际GNSS信号

否

是

是

快速旋转天线

是

否

是

多路径

不存在

存在

历元差分消除

天线覆盖范围

覆盖半球

极有限区域

覆盖半球

表选项 

2机械臂相位中心绝对标定估计模型

在设计观测模型时，本文采用已知基线法[14]，在已知坐标的超短基线上进行同步观测，主站放置型号已知的参考天线，待测站上由机械臂连接待标定天线。

2.1数学模型

记基准天线B对i星的单频载波相位观测方程

（2）

式中, ϕBi为卫星的载波相位观测值；λ为波长；（ρBi）0为天线ARP到卫星的几何距离；δti为卫星钟差；δtB为接收机钟差；IonBi、TropBi分别为电离层延迟和对流层延迟；NBi为B测站对应i卫星的整周模糊度；WB为地基相位缠绕误差[15]；e为卫星信号矢量在AFCS下的单位矢量；εB为多路径和观测噪声。

同理待标定天线C的观测方程为

 （3）

如图2，与基准站不同的是，待测站方程中还包括了机械臂的固定旋转轴点F到待测天线ARP点的转动臂长在卫星矢量上的投影l=len·eCi改正。

图2 双差模型Fig.2 Doubledifferencemodel

图选项 

为了消除卫星和接收机钟差，在星间和测站间组成双差模型。

在超短基线下，对流层和电离层误差项也将被消除，并将基准站天线的已知PCO和PCV改正（或当基准站天线未知时其PCO和PCV的影响可以通过历元间（<10s）差分消除[9]）。

则双差方程为

 （4）

式中，∇ΔϕC－Bj－i=（ϕCj－ϕBj）－（ϕCi－ϕBi）；∇ΔNC－Bj－i为双差模糊度；ΔPCVCj－i=PCVCj－PCVCi；eCj-i=eCj-eCi；f=B1, B2, B3。

将已知坐标的几何距离、机械臂长、相位缠绕改正后，鉴于在天线旋转过程中卫星升降频繁，因此本模型中的模糊度采用单历元固定的策略[16-19]，避免了天线快速转动引起的周跳对参数估计的影响。

在式（4）中，PCV量级在毫米至厘米量级，测地型天线PCO根据天线尺寸不同，其U分量从数厘米至分米级不等。

利用短时间的静态基线解算与已知基线比较可以得到初始PCO分量，故可以用最小二乘解算得到精度较高的浮点解，然后基于LAMBDA算法单历元固定双差模糊度。

在分别固定3个频点的双差模糊度后，则单个频率的相位双差残差为

 （5）

由于天线在快速转动和摆动，在相邻的两个历元间（<10s），多路径噪声可以视为相同，因此通过历元间求差可以有效消除多路径。

最终得到历元差分观测方程为

 （6）

对于待标定天线C，天线在相邻观测历元进行了旋转，因此∇Δv（t2, t1）的变化包含了PCO和PCV的变化信息。

第1阶段求解PCO。

根据前文所述PCO可以视为是APC的平均值，因此在解算PCO的过程中首先忽略PCV项，记∇Δv（t2）－∇Δv（t1）为L（t2, t1）,得到历元间差分观测方程为

 （7）

根据大量卫星矢量双差和历元差数据，将式（7）中的矢量展开为分量形式，利用最小二乘求解得到PCO的NEU3个分量。

记∇ΔeCji（t2, t1）的N、E、U分量分别为∇ΔeNji、∇ΔeEji、∇ΔeUji，则误差方程的矩阵形式为

 （8）

 （9）

式中，A、L矩阵的下标代表参考星；P阵为权阵，其形式见2.2节随机模型。

最后根据最小二乘原理求解

得到PCO3个分量。

第2阶段估计PCV。

待得到PCO估计值后，将其值作为已知值消去，则历元差残差可写为

 （10）

为了获得均匀覆盖半球的PCV模型，可以采用多项式拟合和球谐函数拟合。

当只考虑卫星高度角的影响时，采用四阶多项式拟合效果最佳[3]。

当同时考虑高度角和方位角时，可采用文献[20]提出的相位中心变化球谐函数（式11），利用最小二乘反算系数Anm和Bnm

 （11）

式（10）是4项球谐函数的代数和，因此消除了0阶项，在构造设计矩阵时，需要从n=1开始。

此外，为了满足PCV在天顶处值为零的定义，需要附加约束条件：

。

相关研究表明，针对与高度角和方位角均相关的PCV模型，将式（11）展开成8阶5次（n=8, m=5）即可达到拟合精度要求[9]。

在求解完系数后，采用IGS发布的标准天线ATX文件格式[10]，方位角从0~360°以5°为步长，高度角从0~90°以5°为步长，分别计算上述格网点的PCV改正值。

2.2随机模型

对于原始相位观测量，通常采用卫星高度角定权[21-22]，这是基于低高度角卫星受大气影响较大、天线极化较弱、信噪比较低故而观测量的先验方差设置较大的先验假设。

在本文天线旋转时，低高度卫星虽然受大气影响较大，但可能会接近天线的天顶方向，此时极化强度可能会大于高度角较高的卫星，因此本文采用高度角和信噪比混合模型[23]，且将截止高度角提高到25°来削弱大气的影响。

设载波相位观测量为Φ，E为单位阵, σ0为单位权中误差。

对于BDS系统的3种类型卫星GEO/IGSO/GEOMEO，σ0分别取为0.003、0.003、0.009m[24]，则m颗卫星非差观测值的方差协方差阵为Dϕ=σ02Em。

设TSD为单差算子，TDD为双差算子，在超短基线中，可认为同一卫星对A、B两站等权，即DϕA=DϕB,故根据误差传播律可得到一个历元两个测站同时观测m颗卫星的单差观测量先验方差协方差矩阵DΔϕ和双差观测量先验方差协方差矩阵D∇Δϕ分别为

 （12）

 （13）

若不考虑不同观测历元观测量的相关性，测站A、B同时观测m颗卫星n个历元的双差观测量的方差协方差阵为

 （14）

设三差算子为TTD，则m颗卫星n个历元的三差观测量的方差协方差阵为[25]

则2.1节数学模型中的权阵P=Dδ∇Δϕ-1。

3标定策略

在2.1节的数学模型中，需要对所有项进行时空统一转换。

对于接收机时钟和机械臂控制器时间需要进行同步；对于在不同坐标框架下的矢量，需要统一归算到天线坐标系下。

（1）时间同步。

机械臂控制器由PC控制，其旋转到每一个位置均会输出此位置的起始时间、结束时间、天线高度角和方位角等信息。

但是此输出时间是PC系统时间，精度较差，而接收机观测数据的时间是GPST时。

为了实现接收机时间和PC输出时间同步，本文采用串口通信的方式，通过接收机串口读取时间戳来校正PC系统时间,校正间隔为30s，校正精度<0.1s。

（2）测站坐标。

为了得到基准站天线ARP和机械臂旋转轴点的毫米级先验坐标，利用已知天线先进行24h同步观测。

对参考站和测试站进行PPP和基线解算，最终得到毫米级的先验坐标矢量。

（3）当地水平坐标系与天线固定坐标系转换。

在上述模型中，卫星矢量是以天线固定坐标系（AFCS）为参考。

在WGS-84坐标系下得到的站星矢量经转化后可以归算到当地水平坐标系（locallevelsystem，LLS）下。

由LLS再转换到AFCS需要3步旋转。

如图3所示，第1个旋转角是机械臂安装固定后其Y轴与LLS北方向的夹角，记为偏航角y1。

如图4所示，第2个旋转角是天线的高度角，记为俯仰角p。

第3个旋转角是机械臂控制的天线旋转方位角，在图4中是天线北刻度与N″的夹角y2。

经过3步旋转即可得到卫星矢量在AFCS系下的矢量e。

上述旋转涉及的旋转矩阵如下

图3 机械臂旋转Fig.3 Robotrotation

图选项 

图4 坐标旋转Fig.4 Coordinaterotation

图选项 

（4）机械臂运动方式。

鉴于PCV由球谐函数表达，其自变量为方位角和高度角，因此天线应该在这两个维度分别旋转。

第1步，外循环为高度角，机械臂绕其固定轴点倾斜旋转；内循环为方位角，天线绕其自转轴旋转；此过程为分离出方位角变化对PCV的影响。

第2步，外循环为方位角，内循环为高度角。

此过程为分离出高度角变化对PCV的影响。

天线在旋转过程中，某些时刻卫星不可避免地会接近0°角入射。

为避免信号在经过物体边缘处发生衍射效应，此时应将该颗卫星剔除。

（5）历元差分消除多路径效果。

为了验证历元间差分对多路径消除效果，本文采集了采样率为1Hz的3h静态数据。

固定先验坐标后G29的双差残差如图5（a）。

从图5（b）的时间序列可以看出，历元间差分有效削弱了图5（a）中2:

00的波动趋势，呈现白噪声分布。

图5 多路径时间序列Fig.5 Multipathtimeseries

图选项 

4实例与结果4.1标定过程

试验选取了两台TrimbleR9接收机和两个Trimble59800扼流圈天线。

首先进行了24h的水平静态测量，由于天线类型相同，因此在未获得天线PCO改正信息的情况下利用Bernese解算同样可以获得毫米级的基线坐标矢量，解算得到的基线长度约为3.597m,验后中误差ENU分别为0.2、0.2、0.4mm。

然后在2018年DOY132—136进行了标定试验（见图6）。

每次标定过程采样率设置为1Hz，持续时间约为3h。

天线固定坐标系截止高度角根据北斗/全球卫星导航系统（GNSS）测量型天线性能要求及测试方法[26]中的相位中心偏差测试方法设定为5°。

图6 标定现场Fig.6 Calibrationsite

图选项 

分析图7可以发现，相较于静态观测卫星轨迹只能覆盖有限的区域，使用精密旋转机械臂控制待测天线其卫星轨迹分布均匀覆盖天线半球，能在数小时内得到几千个不同方向上的观测值分布。

图7 静态观测情况下和利用机械臂观测的数据分布对比Fig.7 Comparisonofdatadistributionbetweenstaticobservationandobservationwithroboticarm

图选项 

4.2标定结果

4.2.1PCO估计结果

为了验证标定流程和算法，本文首先对Trimble59800的L1/L2两个频点进行标定，并与IGS发布的结果进行比较。

连续5d的观测数据计算得到的Trimble59800的GPSPCO估计值见表2。

其中横向两组NEU是基准站分别用已知和未知型号天线的标定结果。

综合表2数据可以看出，待测天线PCO标定结果几乎不依赖于主站天线PCO改正信息。

GPSPCO两个频率标定值STD小于0.5mm，其与IGS发布值的RMS小于1mm。

表2 GPSPCOL1/L2多天标定值统计Tab.2 GPSPCOL1/L2multi-daycalibrationstatistics 

mm

GPSL1频率PCO

项目

基准站天线型号已知

基准站天线型号未知

North

East

Up

North

East

Up

DOY132

0.63

1.18

90.39

0.56

1.09

90.33

DOY133

0.43

0.90

89.90

0.48

0.95

89.84

DOY134

1.20

0.62

89.63

1.11

0.64

89.67

DOY135

0.95

0.98

90.34

0.90

1.05

90.17

DOY136

0.76

0.25

89.68

0.72

0.26

89.63

IGS标定值

1.04

0.70

89.51

1.04

0.70

89.51

平均值

0.79

0.79

89.99

0.75

0.80

89.93

RMS

0.36

0.33

0.58

0.37

0.33

0.50

STD

0.26

0.32

0.32

0.23

0.31

0.28

GPSL2频率PCO

项目

基准站天线型号已知

基准站天线型号未知

North

East

Up

North

East

Up

DOY132

0.14

-0.21

116.56

0.10

-0.28

116.49

DOY133

0.02

0.18

116.26

-0.01

0.13

116.27

DOY134

-0.25

0.59

116.89

-0.22

0.64

116.88