北斗三号卫星多频多通道差分码偏差估计与分析.docx
- 文档编号:8940588
- 上传时间:2023-02-02
- 格式:DOCX
- 页数:18
- 大小:951.71KB
北斗三号卫星多频多通道差分码偏差估计与分析.docx
《北斗三号卫星多频多通道差分码偏差估计与分析.docx》由会员分享,可在线阅读,更多相关《北斗三号卫星多频多通道差分码偏差估计与分析.docx(18页珍藏版)》请在冰豆网上搜索。
北斗三号卫星多频多通道差分码偏差估计与分析
摘要:
差分码偏差(DCB)是电离层建模与导航定位授时的主要误差源,北斗多频多通道信号衍生出一系列新的DCB。
本文首先分析了北斗三号卫星的码观测值组合及可估的DCB类型,建立了北斗三号卫星多频码偏差估计的数学模型,利用IGS实测数据首次估计得到了22种不同类型的北斗DCB。
在此基础上,全面比较分析了各类DCB的内符合精度、外符合精度及月稳定度。
结果表明,北斗三号卫星各类DCB的闭合差基本都在0.2ns以内,具有较好的内符合精度;估计结果与中科院(CAS)、德国宇航中心(DLR)提供的DCB产品具有一致性,与CAS的6种DCB偏差基本在0.1ns以内,与DLR的4种DCB偏差基本在0.2ns以内;由于误差传递的影响,通过线性转换得到DCB值的精度和可靠性不及DCB直接估计量;北斗三号卫星各类DCB的月平均标准差为0.083ns,具有较好的中长期稳定性;相较于北斗二号卫星,北斗三号卫星的DCB稳定性相对更优。
关键词:
北斗三号全球导航卫星系统 多频多通道 差分码偏差 一致性 稳定性
Estimationandanalysisofthemulti-frequencyandmulti-channelDCBforBDS-3
Abstract:
Differentialcodebias(DCB)isoneofthemajorerrorsinionosphericmodeling,satellitenavigation,positioning,andtiming.AnewseriesofDCBsisderivedfromBDSmulti-frequencyandmulti-channelsignals.Firstly,thispaperanalyzesthecodeobservationcombinationandestimableDCBtypeforBDS-3,establishesthemathematicalmodelofmulti-frequencyandmulti-channelDCBestimation,andestimatesmorethan20typesofDCBbyusingIGSdata.Onthisbasis,theprecision,accuracy,andmonthlystabilityofvariousDCBsarecomparedandanalyzedcomprehensively.TheresultsindicatethattheclosureerrorsofBDS-3DCBsarebasicallywithin0.2ns,whichshowsgoodprecision.TheestimatedresultshaveagoodagreementwiththeDCBproductsprovidedbyCASandDLR.SixtypesofDCBdifferenceswithCASarebasicallywithin0.1ns.FourtypesofDCBdifferenceswithDLRarebasicallywithin0.2ns.Duetotheinfluenceoferrorpropagation,theaccuracyandreliabilityofDCBobtainedbylineartransformationarenotasgoodasDCBestimateddirectly.ThemonthlymeanSTDofBDS-3DCBsis0.083ns,showinggoodmedium-andlong-termstability.ComparedwiththatofBDS-2,theDCBstabilityofBDS-3isrelativelybetter.
Keywords:
BDS-3 multi-frequencyandmulti-channel differentialcodebias consistency stability
GNSS导航信号差分码偏差(differentialcodebias,DCB)是指由于硬件延迟导致同一时刻不同频率或同一频率上不同测码信号之间的时延差异,包括卫星端差分码偏差和接收机端差分码偏差[1]。
差分码偏差是影响电离层总电子含量(totalelectroncontent,TEC)监测和建模的主要误差源[2-4],忽略卫星和接收机DCB会导致TEC计算误差达到数十个纳秒[4-5]。
差分码偏差也直接影响利用伪距进行导航定位与授时的精度,其误差可达数米[6-8]。
因此,在电离层延迟估计与建模、导航定位与授时等应用中必须准确分离卫星和接收机DCB。
自1998年以来,国际GNSS服务组织(InternationalGNSSService,IGS)启动了电离层工作组,将GPS卫星DCB作为其电离层模型的副产品。
从2003年起,IGS分析中心开始提供GLONASS卫星的DCB[9]。
文献[10-12]基于GPS观测数据对GPS卫星和接收机的DCB参数进行了一系列的研究,解算结果与IGS发布的结果差异为亚纳秒级。
近年来,随着我国北斗卫星导航系统的发展,文献[3]提出采用“两步法”,利用区域电离层模型通过较少的跟踪站即可精确估计北斗二号卫星的DCB。
在此基础上,文献[13]在MGEX框架(multi-GNSSexperiment,MGEX)下实现了多模GNSS的差分码偏差估计,取得了与欧洲定轨中心(CentreforOrbitDeterminationinEurope,CODE)、德国宇航中心(DeutschesZentrumfürLuft-undRaumfahrt,DLR)等机构相当的精度。
2020年6月23日,随着最后一颗北斗三号组网卫星成功送入预定轨道,北斗三号全球导航卫星系统星座部署全面完成。
相较于北斗二号卫星,北斗三号卫星提供的频率资源和调制方式更加丰富。
它不仅继承了北斗二号卫星的频率,还增加了两个互操作频点B1C和B2a,且同一频率具有多个支路,其频点与伪距通道信息见表1[14-19]。
北斗三号更多的可用卫星和更加丰富的频率资源为导航定位带来了新的发展机遇。
与此同时,多频率多通道信号也衍生出一系列新的偏差亟需解决,而差分码偏差就是其中之一。
早在北斗三号系统建设之初,已有学者对北斗三号卫星的DCB估计方法及其特性进行了研究[2, 15, 20-21]。
但受限于当时可用的北斗三号卫星数量较少,且地面接收设备观测值类型的限制,绝大多数研究只聚焦于B1Ⅰ、B2Ⅰ和B3Ⅰ3个频点的少数几类DCB。
截至目前,全球可持续稳定提供DCB产品的机构中国科学院(ChineseAcademyofScience,CAS)和DLR分别只提供8种、6种北斗卫星DCB,具体的DCB类型见表2。
尽管通过线性组合方式可重构计算部分其他类型的DCB,但仍无法涵盖所有类型的DCB。
当前,随着北斗三号系统的全面组网和地面接收机的更新和升级,使得北斗三号多频多通道的DCB研究成为可能。
因此,本文针对北斗三号多频多通道的差分码偏差估计研究,对于提升北斗导航定位与授时服务性能具有重要的现实意义。
表1 BDS卫星信号频率与伪距类型[14-19]Tab.1 BDSsatellitessignalsfrequencyandcodeobserva-tions[14-19]
信号
频率/MHz
伪距
BDS-2
BDS-3
B1Ⅰ
1561.098
C2ⅠC2QC2X
B3Ⅰ
1268.52
C6ⅠC6QC6X
B2Ⅰ
1207.140
C7Ⅰ
-
C7Q
-
C7X
-
B1C
1575.42
-
C1D
-
C1P
-
C1X
B2a
1176.45
-
C5D
-
C5P
-
C5X
B2b
1207.140
-
C7D
-
C7P
-
C7Z
B2(B2a+
B2b)
-
C8D
1191.795
-
C8P
-
C8X
表选项
表2 CAS和DLR北斗卫星DCB类型Tab.2 TheBDSsatelliteDCBtypesofCASandDLR
机构
CAS
DLR
DCB类型
C1P-C5PC1P-C6ⅠC1X-C5X
C1X-C6ⅠC1X-C7ZC1X-C8X
C2Ⅰ-C7ⅠC2Ⅰ-C6Ⅰ
C1X-C2ⅠC2Ⅰ-C5XC2Ⅰ-C8X
C2Ⅰ-C6ⅠC2Ⅰ-C7ⅠC2Ⅰ-C7Z
表选项
1差分码偏差估计方法1.1DCB估计数学模型
忽略伪距观测噪声和多路径误差,任意两种类型的码观测值的观测方程可以表示为[15]
(1)
式中,P表示码观测值;上标s和下标r分别表示卫星端和接收机端;Li和Lj分别表示频率i和频率j;x和y分别表示码观测值为x类型和y类型;ρrs表示卫星与测站间的几何距离;IsLi和IsLj分别表示Li和Lj频率上的电离层延迟;Trs表示对流层延迟;br, Li, x和bsLi, x分别表示卫星端和接收机端PLi, xs码观测值的硬件延迟;br, Lj, y和bLj, ys分别表示卫星端和接收机端PLj, ys码观测值的硬件延迟。
将式
(1)中两种码观测值的观测方程作差,形成无几何距离(geometry-free,GF)组合观测量
(2)
由于电离层延迟与频率有关,当i=j时,即两种码观测值在同一频率上,电离层延迟相同,IsLi-IsLj=0。
当i≠j时,即两种码观测值在不同频率上,电离层延迟无法消除。
此时DCB估计方法主要有两种[13]。
一种方法是在电离层建模时,同步估计频间DCB参数。
另一种方法是采用经验或已知的电离层模型修正电离层延迟,然后再估计频间DCB参数。
CAS和DLR分别采用第一种和第二种估计方法。
对于全球电离层TEC模型,第一种方法计算量大且只能估计参与计算的信号间DCB;对于区域电离层TEC模型,卫星DCB解算受区域模型的精度影响较大。
由于北斗观测站数量和分布的限制,仅使用北斗观测数据难以确定精确的电离层模型[22]。
因此本文采用第二种方法,利用已有的高精度电离层格网产品(globalionospheremaps,GIM)改正电离层延迟,进而估计频间DCB。
通常认为一天内的DCB为常数,将式
(2)中无几何组合观测量在一天内取平均,削弱观测噪声和多路径误差,得到综合的DCB观测值
(3)
式中,DCBr, Li, xLj, y和DCBsLi, xLj, y分别为待估的接收机端和卫星端DCB参数;
;VTECrs为电离层穿刺点处垂直方向总电子含量;M(Ζrs)为投影函数;Ζ为接收机位置的卫星天顶距。
由于卫星端和接收机端DCB参数线性相关,所形成的法方程秩亏。
通常将该类型所有卫星端DCB之和约束为0,即添加“零基准”约束,其约束方程可以表示为
(4)
式中,Ns为该类型DCB所包含的卫星总数。
约束方程的不同会导致DCB估计结果的不同,而不同机构采用的约束方程可能不一致,同一机构不同时间受观测卫星数量变化影响其基准也会产生变化。
因此,不同机构或不同时间之间的DCB要进行基准统一才能进行比较[23-26]。
1.2DCB估值基准统一
假设A基准有m颗卫星,B基准有A的m颗卫星中的n颗卫星(n < m)。
A基准下,对m颗卫星施加“零基准”,其约束方程可以表示为
(5)
式中,
为A基准下的卫星端DCB参数。
若要将A基准下的卫星端DCB参数转换到B基准下,则对这n颗卫星施加“零基准”,其约束方程可以表示为
(6)
式中,Sm为B基准下约束方程系数向量,
为B基准下的卫星端DCB。
两个基准间待估参数转换公式可以表示为
(7)
式中,I为m阶单位矩阵。
如果n > m,则采用相同的方式,把B基准下的DCB参数调整到A基准下。
如果要分析一段时间跨度内DCB的稳定性,由于会存在某些天某些卫星值缺失的情况,选取所有天都有效的卫星施加零基准约束,再将每天的基准统一到此零基准,只有这样,才能对一段时间内的DCB进行稳定性分析。
2试验数据与处理策略2.1试验数据获取
当前,大部分支持北斗的测站均能接收到C2Ⅰ、C6Ⅰ、C7Ⅰ这3种码观测值,三者之间的DCB容易估计。
为了确保其余各类DCB的可估性,首先需要搜集尽可能多的北斗多频多通道码观测值。
其次,为了客观评价本文DCB估值与参考值之间的一致性和稳定性,还需要选取与DLR、CAS尽可能相同的测站数及分布。
表3为利用2020年4月1日至2020年4月30日(对应DOY92-DOY121)一个月内所有IGS站网数据,统计得到的可估DCB类型及可用测站数。
将这些测站取并集,综合考虑测站的空间分布与观测质量,得到本文试验所用的58个测站,如图1所示,图中不同颜色代表不同类型的接收机,这与DLR和CAS所选的测站基本一致。
表3 北斗三号卫星可估的DCB类型及可用测站数(不含C2Ⅰ-C7Ⅰ、C2Ⅰ-C6Ⅰ、C6Ⅰ-C7Ⅰ)Tab.3 BDS-3satelliteestimableDCBtypesandavailablestationnumbers(excludeC2-C7Ⅰ,C2Ⅰ-C6Ⅰ,C6Ⅰ-C7Ⅰ)
DCB类型
测站数
DCB类型
测站数
C1P-C2Ⅰ
25
C2Ⅰ-C5X
34
C1P-C5P
25
C5X-C6Ⅰ
34
C1P-C6Ⅰ
25
C5X-C7Z
27
C1X-C2Ⅰ
28
C5X-C8X
23
C1X-C5X
28
C2Ⅰ-C8X
23
C1X-C6Ⅰ
28
C6Ⅰ-C8X
23
C1X-C7Z
21
C7Z-C8X
22
C1X-C8X
17
C2Ⅰ-C7Z
28
C2Ⅰ-C5P
25
C6Ⅰ-C7Z
27
C5P-C6Ⅰ
25
表选项
图1 测站分布Fig.1 Distributionofstations
图选项
本文使用CODE发布的GIM产品改正电离层延迟。
CODE采用基于球谐函数的全球电离层TEC建模方法计算GIM,其内符合精度约为2.4TECU[27-28]。
2.2数据处理策略
图2给出了DCB估计的数据处理流程。
其中,质量控制部分主要包括设置卫星截止高度角为20°、剔除波动大(标准差大于1)和观测数量少(时长小于50min)的GF组合时间序列;最小二乘中单位权中误差σ0=0.1,依据标准差给综合DCB观测值定权。
图2 DCB估计Fig.2 DCBestimation
图选项
由表3可知,可估的北斗三号DCB类型共有19种,各类码组合的可用测站数为17~34个。
为了与北斗二号进行比较,同时还估计了C2Ⅰ-C6Ⅰ、C2Ⅰ-C7Ⅰ、C6Ⅰ-C7Ⅰ这3种DCB类型,总共估计得到了22种类型的DCB。
其中C2Ⅰ-C7Ⅰ、C6Ⅰ-C7Ⅰ类型的DCB只存在于北斗二号卫星,C2Ⅰ-C6Ⅰ类型的DCB同时存在于北斗二号和北斗三号卫星,其余类型的DCB则只存在于北斗三号卫星。
为了评估北斗三号多频多通道DCB估计的精度和可靠性,从内符合精度、外符合精度、估值稳定性等方面进行分析。
需要注意的是,在比较分析之前需要采用1.2节中的方法消除由不同机构的基准差异、同一机构不同时间卫星数量变化引起的基准不一致。
3试验结果与分析3.1内符合精度分析
闭合差是反映内符合精度的重要指标,给定任意3个码观测值,可构造3个DCB,但是,理论上只有两个独立的DCB,即3个DCB之间的闭合差为零[29-30]。
由于各类DCB是基于不同观测值独立解算的,且平差计算时受观测噪声和模型误差影响,导致其闭合差通常并不等于零。
闭合差的大小反映了DCB估值与模型的吻合程度,闭合差越小,内符合精度越高。
以北斗三号卫星3个新频点之间形成的3类DCB(C1X-C5X、C1X-C8X、C5X-C8X)为例,图3给出了北斗三号各卫星的每日的DCB闭合差(DCBC1X-C5X+DCBC5X-C8X+DCBC8X-C1X)时间序列。
由图3可知,各卫星每日的DCB闭合差在零值附近波动,无显著系统性偏差,且绝大多数闭合差分布在0.3ns以内。
少数卫星如C28、C32、C35卫星分别在DOY113、DOY120、DOY97的闭合差达到1ns,这可能与当天可用的观测数量和数据质量有关。
统计北斗三号各卫星DCB月平均闭合差及其标准差,如图4所示。
从图4中可以看出,北斗三号卫星DCB估值的内符合精度较好,平均闭合差都在0.2ns以内,且大部分优于0.1ns,其中C40卫星平均闭合差最小为0.008ns,C36卫星的平均闭合差最大为0.2ns。
各卫星闭合差的标准差都在0.4ns以内,且大部分优于0.2ns,这表明C1X-C5X、C1X-C8X、C5X-C8X3个DCB估值之间具有较好的一致性,内部符合精度较高。
其他类型的DCB闭合差具有类似的统计特性,限于篇幅,不单独罗列。
图3 北斗三号各卫星DCB日闭合差时间序列(C1X-C5X、C1X-C8X、C5X-C8X)Fig.3 TimeseriesofBDS-3satelliteDCBdailyclosureerrors(C1X-C5X、C1X-C8X、C5X-C8X)
图选项
图4 北斗三号各卫星DCB月平均闭合差及标准差(C1X-C5X、C1X-C8X、C5X-C8X)Fig.4 MonthlymeanclosureerrorsandSTDsofBDS-3satelliteDCBs(C1X-C5X、C1X-C8X、C5X-C8X)
图选项
3.2外符合精度分析
为了验证北斗三号卫星DCB估计结果的外符合精度,将本文的估计结果(下文简称为WHU)与CAS和DLR机构的DCB产品进行比较。
除C2Ⅰ-C6Ⅰ之外,CAS和DLR的DCB产品中分别有6种、4种北斗三号卫星DCB类型,将其依次与WHU求差,得到北斗三号各卫星DCB估值的外符合精度,如图5、图6所示。
从图5中可以看出WHU和CAS的平均偏差基本都在0.1ns以内,这说明WHU与CAS的DCB产品之间一致性较好。
其中C25卫星C1X-C8X类型DCB平均偏差较大,达到0.3ns,这与C1X-C8X类型DCB的数据质量有关,见表3,用于估计C1X-C8X类DCB的测站数量仅为17,明显少于其他DCB类型。
从图6中可以看出,除C45卫星之外,WHU与DLR之间DCB的平均偏差基本在0.2ns以内,说明二者之间也具有较好的一致性。
至于C45卫星,DLR和WHU之间C1X-C2Ⅰ、C2Ⅰ-C5X两种类型的DCB差值分别为-0.88、1.15ns。
这与该卫星的数据、质量和这两类DCB的稳定性有关,分析DLR和WHU的DCB估值方差-协方差矩阵发现,所有卫星C1X-C2Ⅰ的月平均标准差分别为0.054、0.051ns,而C45卫星的标准差分别为0.308、0.155ns;类似的,C45卫星C2Ⅰ-C5X的月平均标准差分别为0.135、0.213ns,显著大于其他卫星。
图5 WHU与CAS估计的6种北斗三号卫星DCB之间的平均偏差Fig.5 The6typesofBDS-3satelliteDCBmeandifferencesofWHUandCAS
图选项
图6 WHU与DLR估计的4种北斗三号卫星DCB之间的平均偏差Fig.6 The4typesofBDS-3satelliteDCBmeandifferencesbetweenWHUandDLR
图选项
尽管CAS和DLR只提供少数几种DCB产品类型,但是根据DCB之间的线性关系,可利用已有的DCB类型经过转换得到一些新的DCB类型,比如DCBC1X-C8X=DCBC1X-C5X+DCBC5X-C8X。
因此,还可将WHU的DCB估值与CAS和DLR的部分DCB转换值进行间接比较,如图7、图8所示。
从图7和图8中可以看出,相比于图4和图5中DCB估值的直接比较法,由于受误差传播的影响,间接法获得的各机构之间的DCB差值明显更大,一致性较差,其差值分布在1ns之内,这也说明了直接估计多频多通道DCB的必要性。
图7 WHU与CAS的13种北斗三号卫星DCB(转换值)之间的平均偏差Fig.7 The13typesofBDS-3satelliteDCB(transformvalues)meandifferencesbetweenWHUandCAS
图选项
图8 WHU与DLR的7种北斗三号卫星DCB(转换值)之间的平均偏差Fig.8 The7typesofBDS-3satelliteDCB(transformvalues)meandifferencesbetweenWHUandDLR
图选项
3.3估值稳定性分析
为分析北斗三号卫星DCB估值的稳定性,首先选取3个新增频点上的码观测值C1X、C5X和C8X构成的DCB类型C1X-C5X、C1X-C8X为例,将2020年4月内各卫星对应的DCB每日估值的时间序列绘于图9。
可以看出,大多数卫星的DCB估值分布在±20ns,少数卫星的DCB数值较大,其中C33卫星的DCB估值达到-78.2ns左右,在导航定位与授时中若不考虑,将带来严重的系统偏差。
在一个月的时间内,各卫星的DCB估值变化较小,表明它们具有较好的天稳定度。
图9 BDS-3各卫星C1X-C5X和C1X-C8XDCB时间序列Fig.9 TimeseriesofBDS-3satellitesC1X-C5XandC1X-C8XDCB
图选项
进一步分析各类DCB的稳定性,表4给出WHU、CAS和DLR直接估计或通过线性转换得到的DCB值的月平均标准差(STD)。
其中,黑体表示线性转换值。
由表4可知,WHU的DCB估值稳定性普遍优于CAS和DLR,反映本文估计的DCB稳定性良好。
DLR的平均STD虽优于CAS,但其产品提供的DCB类型较少,且涵盖的码观测值类型也较少,导致某些DCB类型无法通过线性转换得到。
利用CAS的DCB产品,虽然可通过线性转换恢复所有DCB类型,但因受误差传递和累积的影响,转换得到的DCB稳定性明显较差,进一步说明了直接估计多频多通道DCB的必要性。
表4 北斗卫星DCB的月平均STDTab.4 MonthlymeanSTDofBDSsatelliteDCBs ns
DCB类型
WHU
CAS
DLR
C1P-C2Ⅰ
0.026
0.069
-
C1P-C5P
0.
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 北斗 三号 卫星 多频多 通道 差分码 偏差 估计 分析