[发明专利]一种GNSS接收机高采样钟差解算方法

权利要求书

1.一种GNSS接收机高采样钟差计算方法，所述GNSS接收机为地面静态GNSS接收机，所述方法包括：

使用低采样GNSS接收机观测数据进行精密单点定位解算，获得低采样时刻的第一接收机钟差，以及低采样时刻的接收机位置参数、对流层天顶延迟参数和模糊度参数；

获取上述三类参数在高采样时刻对应的高采样参数值，利用三类高采样参数值计算高采样时刻的第二接收机钟差；

利用第一接收机钟差和第二接收机钟差，计算最优化的高采样接收机钟差；

所述利用第一接收机钟差和第二接收机钟差，计算最优化的高采样接收机钟差；具体包括：

两个相邻低采样时刻的第一接收机钟差为dT^s(t)、dT^s(t+k)，并且这两个时刻第二接收机钟差为dT^m(t)、dT^m(t+k)，则两个时刻内高采样接收机钟差序列为dT^m(t)，dT^m(t+l)…dT^m(t+p*l)…dT^m(t+k)；

则在[t,t+k]时间区间内，高采样接收机钟差为：

其中，dT^f(t)为t时刻的高采样接收机钟差，k表示低采样的时间间隔，l表示高采样率的时间间隔，则在低采样时间区间[t,t+k]内，高采样时间序列个数对于第p个高采样时刻：t+p*l,p＝1,2…P-2。

2.根据权利要求1所述的GNSS接收机高采样钟差计算方法，其特征在于，低采样的时间间隔为30秒～15分钟，高采样的时间间隔为50赫兹～5秒。

3.根据权利要求2所述的GNSS接收机高采样钟差计算方法，其特征在于，所述获取上述三类参数在高采样时刻对应的高采样参数值，具体包括：

对于相邻的两个低采样时刻t和(t+k)，t时刻的第一接收机钟差为dT^s(t)；接收机位置参数为[x_t,y_t,z_t]^T；对流层天顶湿延迟为dZPD_t,wet，模糊度参数为j为观测到的GNSS卫星的序号，n为观测到的GNSS卫星的总数；(t+k)时刻的第一接收机钟差为dT^s(t+k)；接收机位置参数为[x_t+k,y_t+k,z_t+k]^T；对流层天顶湿延迟为dZPD_t+k,wet，模糊度参数为

对于第p个高采样时刻：t+p*l,p＝1,2…P-2，其接收机位置参数为(t+k)时刻的接收机位置参数[x_t+k,y_t+k,z_t+k]^T；

利用t时刻和(t+k)时刻的对流层天顶延迟参数值dZPD_t,wet和dZPD_t+k,wet，进行线性内插得到第p个高采样时刻的对流层天顶延迟参数值；

搜索第p个高采样时刻观测到的GNSS卫星在输出的模糊度文件中所在的起始和结束时间，如果搜索到，则取该范围内模糊度值作为该时刻的模糊度值，否则，模糊度值标识为NAN。

4.根据权利要求3所述的GNSS接收机高采样钟差计算方法，其特征在于，所述利用三类高采样参数值计算高采样时刻的第二接收机钟差；具体包括：

对于第p个高采样时刻：t+p*l,p＝1,2…P-2，利用该时刻的接收机位置参数、对流层天顶延迟参数和模糊度参数，根据观测方程计算第j个GNSS卫星对应的接收机钟差由此得到n个接收机钟差估值为

则第p个高采样时刻的第二接收机钟差为：

其中，Elev_j为第j颗观测到的GNSS卫星相对于接收机的高度角。