GB/T 39396.1-2020.Quality requirements for international GNSS monitoring and assessment system (iGMAS)-Part 1: Observation data.
1范围
GB/T 39396的本部分规定了全球连续监测评估系统(iGMAS)跟踪站GNSS观测数据质量指标要求、评定项目及观测数据要求、评定方法。
GB/T 39396.1适用于iGMAS测量型GNSS接收机研制与指标检定、跟踪站选址、观测数据质量评定与分析等工作。
2规范性引用文件
下列文件对于本文件的应用是必不可少的。凡是注日期的引用文件,仅注日期的版本适用于本文件。凡是不注日期的引用文件,其最新版本(包括所有的修改单)适用于本文件。
GB/T 39267-2020北斗 卫星导航术语
3术语、定义和缩略语
3.1 术语和定义
GB/T 39267-2020界定的以及下列术语和定义适用于本文件。
3.1.1
观测数据 observation data
iGMAS跟踪站GNSS接收机记录的测量数据。
注:包含伪距、载波相位、多普勒频移和信号强度等数据。
3.1.2
观测数据完整率 observation data integrity rate
接收机观测到卫星的实际历元数据量与理论历元数据量的比值。
3.1.3
粗差 outlier
离群的误差。
注:在相同观测条件下作一系列的观测,其绝对值超过限差的测量偏差。
3.1.4
周跳比 observations per slip
在某时间段内,接收机观测数据的实际历元数据量与发生周跳历元数据量的比值。
注:反映了周跳发生的平均观测历元数。
3.1.5
接收机钟跳 receiver clock jump
接收机厂商为保持接收机内部时钟与GNSS时同步精度在一定范围之内,控制接收机钟差不超过某一阈值而插人的时钟跳跃。

ICS49.020
中华人民共和国国家标准
GB/T39396.1—2020
全球连续监测评估系统（iGMAS）质量要求第1部分：观测数据
Quality requirements for international GNSS monitoring and assessmentsystem (iGMAS)—Part 1:Observation data2020-11-19发布
国家市场监督管理总局
国家标准化管理委员会
2021-06-01实施
规范性引用文件
术语、定义和缩略语·
术语和定义
缩略语
观测数据质量指标要求
数据完整率
周跳比
多路径误差
伪距噪声
载波相位噪声
载噪比
评定项目及观测数据要求
评定方法
数据完整率
周跳比
多路径误差
伪距噪声
载波相位噪声
载噪比
附录A（规范性附录）
附录B（资料性附录）
GNSS卫星信号载波频率
双频组合频率选择
GB/T39396.1—2020
GB/T39396《全球连续监测评估系统(iGMAS)质量要求》分为两个部分：第1部分：观测数据；
第2部分：产品
本部分为GB/T39396的第1部分
本部分按照GB/T1.1—2009给出的规则起草。本部分由中央军委装备发展部提出。本部分由全国北斗卫星导航标准化技术委员会（SAC/TC544)归口。GB/T39396.1—2020
本部分起草单位：战略支援部队信息工程大学、中国卫星导航工程中心、长安大学、卫星导航系统与装备技术国家重点实验室、中国航天标准化研究所。本部分主要起草人：李建文、焦文海、李涌涛、刘莹、黄观文、苏牡丹、蔚保国、潘林、伍蔡伦、庞鹏、盛传贞、陈晨、王凯、谢松、干维嘉、张楷时、顾晨钟。m
1范围
GB/T39396.1—2020
全球连续监测评估系统（iGMAS）质量要求第1部分：观测数据
GB/T39396的本部分规定了全球连续监测评估系统（iGMAS)跟踪站GNSS观测数据质量指标要求、评定项目及观测数据要求、评定方法本部分适用于iGMAS测量型GNSS接收机研制与指标检定、跟踪站选址、观测数据质量评定与分析等工作。
2规范性引用文件
下列文件对干本文件的应用是必不可少的。凡是注日期的引用文件，仅注日期的版本适用干本文件。凡是不注日期的引用文件，其最新版本（包括所有的修改单)适用于本文件。GB/T39267—2020北斗卫星导航术语术语、定义和缩略语
术语和定义
GB/T39267一2020界定的以及下列术语和定义适用于本文件。3.1.1
观测数据
observationdata
iGMAS跟踪站GNSS接收机记录的测量数据。注：包含伪距、载波相位、多普勒频移和信号强度等数据3.1.2
观测数据完整率observationdataintegrityrate接收机观测到卫星的实际历元数据量与理论历元数据量的比值。3.1.3
粗差outlier
离群的误差。
注：在相同观测条件下作一系列的观测，其绝对值超过限差的测量偏差。3.1.4
周跳比
observationsperslip
在某时间段内，接收机观测数据的实际历元数据量与发牛周跳历元数据量的比值。注：反映了周跳发生的平均观测历元数。3.1.5
receiverclockjump
接收机钟跳
接收机厂商为保持接收机内部时钟与GNSS时同步精度在一定范用之内，控制接收机钟差不超过某一阈值而插入的时钟跳跃。
GB/T39396.1—2020
多路径误差multipatherror
由非直达导航信号引人的测距误差。3.1.7
pseudorangenoise
伪距噪声
伪距测量随机误差。
注：泛指导航信号产生端、接收端产生或引人的各种随机误差对伪距测量的综合影响值、3.1.8
载波相位噪声carrierphasenoise载波相位测量随机误差。
注：泛指导航信号产生端、接收端产生或引人的各种随机误差对载波相位测量的综合影响值。3.1.9
载噪比
carrier-to-noiseratio;CNR
载波信号功率与噪声功率谱密度之比。注：单位为dBHz。
缩略语
下列缩略语适用于本文件。
BDS：北斗卫星导航系统（BeiDouNavigationSatelliteSystem）Galileo：伽利略卫星导航系统（GalileoNavigationSatelliteSystem）GLONASS：全球导航卫星系统（GlobalNavigationSatelliteSystem）GNSS：全球卫星导航系统（GlobalNavigationSatelliteSystem）GPS：全球定位系统（GlobalPositioningSystem）iGMAS：全球连续监测评估系统（InternationalGNSSMonitoringandAssessmentSystem）4观测数据质量指标要求
1数据完整率
任一GNSS系统任一可观测频率（信号）的观测数据完整率应不小干90%（采样率30s）。任一GNSS系统的系统观测数据完整率应不小于90%（采样率30s）。周跳比
任一GNSS系统的周跳比应不小于500（采样率30s)。4.3
多路径误差
任一GNSS系统任一频率（信号）所有卫星的多路径误差的平均值应不大于0.5m（采样率30s）。伪距噪声
采样率为1s的观测数据，单个信号所有卫星的伪距噪声的平均值应为（载波名称见附录A)：a）0.511MCPS码速率信号（L2C、G1C、G2C)≤18cm；1.023MCPS码速率信号（B1C、L1C/A、L2C/A、L1C、E1)≤12cm；b)
2.046MCPS码速率信号（B11)≤10cm；c）
5.11MCPS、5.115MCPS码速率信号(G1P、G2P、E6)≤8cm；e）10.23MCPS码速率信号（B2a、L1P、L2P、L5C、E5a、E5b、B3I)≤6cm。4.5载波相位噪声
GB/T39396.1—2020
任一GNSS系统任一频率（信号）所有卫星的载波相位噪声的平均值应≤0.01周（采样率1s）。4.6载噪比
任一GNSS系统任一频率（信号）所有卫星的载噪比的平均值应不小于35dBHz（采样率30s）。评定项目及观测数据要求
对于任一GNSS系统，评定项目、观测数据要求及适用的评估对象应符合表1规定。涉及的GNSS卫星信号载波频率见附录A，涉及的双频组合频率选择参见附录B。表1评定项目、观测数据要求及适用的评估对象观测数据要求
数据完整率
周跳比
多路径误差
伪距噪声
载波相位噪声
载噪比
6评定方法
数据完整率
天文件
天文件
天文件
天文件
天文件
天文件
高度角
观测频率
(信号)个数
滤波与平滑
非滤波非平滑
非滤波非平滑
非滤波非平滑
非滤波非平滑
非滤波非平滑
非滤波非平滑
适用的评估对象
按公式（1)计算任一系统s任一频率（信号）的观测数据完整率，按公式（2)计算任一系统的观测数据完整率。
式中：
B.）×100%
ZB:)×100%
GNSS系统s的频率（信号）f观测数据完整率，%；n
在观测时间段内，观测的卫星总数；观测卫星序号，j=1,2，，n；
在观测时间段内，GNSS系统s卫星j在频率（信号）f的实际观测历元总数；在观测时间段内，GNSS系统s卫星i在频率（信号）f的理论历元总数：系统s观测数据完整率，%；
.（1)
·（2）
在观测时间段内，GNSS系统s卫星j所有观测频率（信号）均有观测数据的实际观测历3
GB/T39396.1—2020
元总数；
6.2周跳比
在观测时间段内，GNSS系统s卫星的理论历元总数。按公式（3)计算观测时段内任一GNSS系统s的周跳比。C,=
式中：
Os,total
O s,elips
观测时段内，GNSS系统s的周跳比；-GNSS系统s，在观测时段内观测的实际历元总数：-GNSS系统s，在观测时段内发生的周跳历元总数（周跳数）。·（3)
周跳数osclips采用先MW（Melbourne-Wubbena）组合探测，再GF（Geometry-Free）组合探测的方法：
MW组合探测周跳方法
按公式（4），得到历元t时刻的MW组合观测量。Pmw(t.)=Φi(t,）—2(t.)
式中：
p2(t:)
P2(t,)
fi-f[e(t)
历元t：时刻的MW组合观测量，单位为周；02(t.)
·（4)
历元t时刻的MW组合观测量第一个频率（信号）载波相位噪声观测量，单位为周；
历元t时刻的MW组合观测量第一个频率（信号）载波相位噪声观测量，单位为周；
第一个频率（信号）的载波频率（频率选择方法参见附录B）)，单位为兆赫（MHz)；第二个频率（信号）的载波频率（频率选择方法参见附录B），单位为兆赫（MHz)；历元t；时刻的MW组合观测量第一个频率（信号）伪距观测量，单位为米（m）；第一个频率（信号)的载波的波长，单位为兆赫（MHz)；历元t：时刻的MW组合观测量第二个频率（信号）伪距观测量，单位为米（m）；第二个频率（信号)的载波的波长，单位为兆赫（MHz)。按公式（5)构造MW组合周跳检验量DMw（t,），当DMw（t,)≥1时，在排除粗差的前提下，该历元存在周跳。
Dmw(t,)=mw(t,)—mw(ti-1)
式中：
Dmw(t,)
$pmw(t)
Pmw(ti-,)
历元t，时刻MW组合周跳检验量，单位为周；历元t，时刻MW组合观测量，单位为周：历元t;-时刻MW组合观测量，单位为周。GF组合探测周跳方法
按公式（6），得到历元t，时刻的GF观测量。PGF(t）=(t)
式中：
PGF(t,)
历元t；时刻的GF组合观测量，单位为周；（5）
（6）
历元t，时刻的GF组合观测量第一个频率（信号）载波相位观测量，单位为周；入2
P2(t,)
第二个频率（信号）的载波的波长，单位为兆赫（MHz）；第一个频率（信号）的载波的波长，单位为兆赫（MHz)；GB/T39396.1—2020
-历元t，时刻的GF组合观测量第一个频率（信号）载波相位观测量，单位为周。按公式（7)构造GF组合周跳检验量Dcr（t,），当Dcr（t)≥1时，在排除粗差和接收机钟跳的前提下，该历元存在周跳。
式中：
DGF(t,)
PGF(t:)
PGF(ti-i)
6.3多路径误差
DGF(t)=PGF(t,）—GF(t-)
历元t：时刻GF组合周跳检验量，单位为周；历元t，时刻GF组合观测量，单位为周；历元t,-,时刻GF组合观测量，单位为周。（7）
本部分提供的多路径误差检验方法仅适用于GNSS双频、多频观测数据。按公式（8)计算观测时段内任一GNSS系统、任一频率（信号）、任一卫星的多路径误差值SZAC
式中：
S[M,(t)-M]
观测时段内，观测频率(信号)f的观测值多路径误差值，单位为米（m）；观测时段内，观测历元总数；
历元序号，i=1，2，.，N；
观测频率（信号）f，在历元t，时刻的多路径计算值（含整周模糊度影响），单位为米(m);
观测频率（信号)f，在观测时段内多路径计算平均值（含整周模糊度影响），单位为米(m)。
历元t时刻的多路径计算值M，即M,（t），按公式（9)计算得到。Mn=pi
式中：
fl频率（信号）任一观测历元多路径计算值（含整周模糊度影响），单位为米（m）；0
第一个频率（信号)对应历元的伪距观测量，单位为米（m）；第一个频率（信号）)的载波频率（频率选择方法参见附录B)，单位为兆赫（MHz)；第二个频率（信号）的载波频率（频率选择方法参见附录B)，单位为兆赫（MHz)；第一个频率（信号）对应历元的载波相位观测量，单位为米（m）；第二个频率（信号）对应历元的载波相位观测量，单位为米（m）；f2频率（信号)任一观测历元多路径计算值（含整周模糊度影响），单位为米（m)；第二个频率（信号)对应历元的伪距观测量，单位为米（m）。M,按公式(10)计算得到。
式中：
观测时段内，观测历元总数；
（9）
...（10)
GB/T39396.1—2020
历元序号，i=1,2，...，N。
在多路径误差计算时，应首先在观测数据中剔除粗差、修复接收机钟跳和周跳等的影响。6.4伪距噪声
本部分伪距噪声检验方法仅适用干采样间隔为1s的伪距观测数据。按公式（11)计算观测时段内，任一GNSS系统任一频率（信号）所有观测卫星的伪距噪声（平均值）。XApks
式中：
观测时段内，任一GNSS系统任一频率（信号）伪距噪声（平均值），单位为米（m）；观测时段内，任一GNSS系统任一频率（信号）观测的卫星总数；观测卫星序号，j=1,2，，n；
（11）
观测时段内，任一GNSS系统任一频率（信号）第颗卫星的伪距噪声，单位为米（m）。A0kMs按公式(12)计算得到。
式中：
Ao(t)—免费标准下载网bzxz
C[（t)2
观测时段内，任一GNSS系统任一频率（信号）第i颗卫星的观测历元总数；历元序号，i=1,2，，N；
.（12)
任一GNSS系统任一频率（信号）第i颗卫星，在历元t，时刻的伪距噪声估值，单位为米(m)。
△o（t,)按公式(13)计算。
Ap(t)=p(t)-pi(t,)
式中：
（13）
任一GNSS系统任一频率（信号）第i颗卫星，在历元t，时刻的伪距观测值，单位为米(m);
任一GNSS系统任一频率（信号）第颗卫星，在历元t，时刻的伪距二次多项式拟合值，oi(t,)
单位为米（m）。
按公式（14)，进行伪距二次多项式拟合；拟合计算得到二次多项式系数后，可得到各历元的伪距二次多项式拟合值（t，）。在观测时段内，从起始历元开始，每120个历元为一个拟合窗口，拟合窗口不重叠。靠近观测时段结尾时，当剩余历元数不小于3个历元观测数据时，可按剩余历元数拟合；当剩余历元数小干3个历元观测数据时，可舍弃不用。当观测数据有中断现象时，可分段分别处理。f(r)=aa2+br+c
式中：
伪距二次多项式拟合函数；
拟合窗口内的历元编号，二次多项式自变量，=0.1，，118，119；二次多项式二次项系数；
二次多项式一次项系数；
二次多项式常数项。
6.5载波相位噪声
.（14）
本部分载波相位噪声检验方法仅适用于采样间隔为1s的载波相位观测数据。按公式（15)计算观测时段内，任一GNSS系统任一频率（信号）所有观测卫星的载波相位噪声（平均值）。
小提示：此标准内容仅展示完整标准里的部分截取内容，若需要完整标准请到上方自行免费下载完整标准文档。