[发明专利]单频GNSS相位稳定性监测的高频逐历元相位差方法

摘要

本发明公开了一种单频GNSS相位稳定性监测的高频逐历元相位差方法，在测站点上利用单频GNSS接收机进行高频观测，计算相邻历元的相位差，利用最小二乘方法，计算连续历元下相邻历元间隔内的测站三维(NEU)坐标改正数时间序列，通过时间序列统计分析，分别取得NEU坐标的平均值和标准偏差，计算异常比率。通过提出的方法解算连续观测异常比率，以此来评估单频GNSS接收机相位观测的稳定性。本发明利用单频GNSS数据，进行单点实时高频GNSS相位数据质量监测；无需高精度钟差、各种模型改正即能进行GNSS相位稳定性的评估。

主权项

一种单频GNSS相位稳定性监测的高频逐历元相位差方法，其特征在于，包括以下步骤：第一步，用单频GNSS接收机进行高频观测，计算相邻两个历元之间的相位差；利用单频GNSS接收机在某一测站P上进行高频观测，在观测卫星多于3颗的条件下，高频采集相位数据；对某一GNSS卫星j在历元t1上对应的相位为在GNSS信号不失锁的情况下，在下一个历元t2GNSS卫星j对应的相位为式(1)和(2)中，ρj(t1)为历元t1时测站P至卫星j的几何距离，即，ρj(t1)=[xj(t1)-xP]2+[yj(t1)-yP]2+[zj(t1)-zP]2]]>ρj(t2)为历元t2时测站P至卫星j的几何距离，即ρj(t2)=[xj(t2)-xP]2+[yj(t2)-yP]2+[zj(t2)-zP]2]]>[xj(t1),yj(t1),zj(t1)]为历元t1时GNSS卫星j的三维空间坐标；[xj(t2),yj(t2),zj(t2)]为历元t2时GNSS卫星j的三维空间坐标；(xP,yP,zP)为测站P在WGS84下的坐标；Sclock(j)和Rclock分别为GNSS卫星j和接收机的钟差；Shardware(j)和Rhardware分别为GNSS卫星j和接收机硬件延迟误差；Dtroposphere(j)和Dionosphere(j)分别为对流层和电离层延迟误差；N(j)为整周模糊度；Emultipath(j)为多路径效应；Ephase center(j)为天线相位中心偏差及其变化影响；ε1是历元t1时GNSS卫星j的残余误差；ε2是历元t2时GNSS卫星j的残余误差；对高频GNSS信号来说，在相邻的两个历元上，钟差、大气延迟、硬件延迟、多路径效应、整周模糊度和天线相位中心偏差及其变化影响基本上保持不变，式(1)减去式(2)，可得式中，是历元t1和历元t2时GNSS卫星j的相位差，v(j)＝ε1‑ε2；通过式(3)得出，对高频GNSS信号来说，上述钟差、大气延迟、硬件延迟、多路径效应、整周模糊度和天线相位中心偏差及其变化影响通过逐历元相位差移去或有效降低；相邻历元相位差的精度为其中，σ是相位精度，测站P的最初坐标是[x0,y0,z0]，坐标改正是[Δx,Δy,Δz]；将式(3)在[x0,y0,z0]处进行泰勒展开，忽略高阶项，获得线性形式为其中和是测站P到GNSS卫星j在历元t1和t2时的初始距离，Δx,Δy,Δz为相邻测站坐标改正数的三个时间序列；第二步，解算Δx,Δy,Δz设在历元t1和t2，测站P同时都观测到的GNSS卫星数为m，那么式(5)用矩阵表达为V＝AX+L  (6)式中A为系数矩阵，A=(x2(1)-x0)ρ02(1)-(x1(1)-x0)ρ01(1)(y2(1)-y0)ρ02(1)-(y1(1)-y0)ρ01(1)(z2(1)-z0)ρ02(1)-(z1(1)-z0)ρ01(1).........(x2(j)-x0)ρ02(j)-(x1(j)-x0)ρ01(j)(y2(j)-y0)ρ02(j)-(y1(j)-y0)ρ01(j)(z2(j)-z0)ρ02(j)-(z1(j)-z0)ρ01(j).........(x2(m)-x0)ρ02(m)-(x1(m)-x0)ρ01(m)(y2(m)-y0)ρ02(m)-(y1(m)-y0)ρ01(m)(z2(m)-z0)ρ02(m)-(z1(m)-z0)ρ01(m),]]>m＞3，利用最小二乘法就可以对式(6)进行解算，获得X的最优估值；式(6)的最小二乘解为X＝(ATPA)‑1ATPL  (7)方差‐协方差矩阵为DX=σ02(ATPA)-1---(8)]]>其中，单位权方差是σ02=VTPVm-3---(9)]]>式中，P是权矩阵，初始权是P0＝I，I为单位矩阵，通过式(9)可以得到新的单位权σ0，p(j)=σ0|v(j)|+c---(10)]]>其中，c是一个大于零的常数，通过公式(7)即可计算X；第三步，利用Δx,Δy,Δz序列评估相位稳定性相邻测站坐标改正数的三个时间序列Δx,Δy和Δz通过式(7)构建，从这三个时间序列中，得到平均值和对应的标准偏差，即MEANΔx,MEANΔy,MEANΔz,STDΔx,STDΔy和STDΔZ，并定义上界值和下界值分别为UB＝MEAN+kSTD  (11)LB＝MEAN‑kSTD  (12)其中，k是一个常数，如果Δx,Δy和Δz的一个解大于对应的上界值或小于对应的下界值，则被认为是非正常值；在一个连续时间段内，通过连续高频观测获得n个异常解；连续观测异常比率AP被定义为AP=nq-1×100%---(13)]]>其中，q是连续观测中的观测历元数；当这个连续时间段非常长，为10天或超过10天时，计算连续观测异常比率AP，其结果作为连续观测异常比率的临界值，把临界值作为判断相位稳定性的一个临界点，用来判断平时GNSS测量时某一段时间内的相位稳定性情况；即在一个时间段内，如果连续观测异常比率AP大于这个临界值，则表明相位是不稳定的，反之则相位是稳定的；在一个GNSS测站上，某一段时间内，根据计算出的连续观测异常比率，结合连续观测异常比率的临界值就可以评估单频相位数据在该某一段时间内的稳定性情况。