1.一种顾及到GLONASS频间偏差的时差监测方法，其特征在于，包括以下步骤：S1、实时输入接收机原始观测数据，包括原始伪距和广播星历；S2、对第一步输入的原始观测数据进行预处理，剔除原始观测值中的异常值，删除重复的电文；S3、选择BDS系统和GLONASS系统组成时差监测方程组成时差监测方程，确定待估参数；S3中，对于任一观测站的接收机观测到某一卫星j，其在第i个频点的伪距观测量表示为：式中， 为BDS系统在观测站r的接收机观测到的某一卫星m，其在第i个频点的原始伪距观测量，以米为单位； 为GLONASS系统在观测站r的接收机观测到的某一卫星n，其在第i个频点的原始伪距观测量，以米为单位；下标r为观测站编号；上标m,n为卫星编号；C,R分别代表BDS系统和GLONASS系统； 和 分别为为信号发射时刻的观测站r到卫星m和n间的几何距离；c为真空中光速； 为以北斗时为基准的接收机钟差； 为BDS系统与GLONASS系统的系统时差；IFB i,0 为GLONASS频间偏差；dt C,j 、dt R,j 分别表示BDS系统、GLONASS系统的卫星钟差； 和 分别为BDS系统中的卫星m和GLONASS系统中的卫星n在L i 频点上到观测站r路径上的电离层延迟， 和 分别为BDS系统中的卫星m和GLONASS系统中的卫星n到观测站r路径上的对流层延迟， 为BDS伪距观测值观测噪声和其他未模型化误差， 为GLONASS伪距观测值的观测噪声和其他未模型化误差，未模型化误差包括对流层残余误差，电离层二阶项误差；将方程(1)进行线性化，可得BDS系统与GLONASS系统的时差监测方程为式中， 为无电离层组合后的从观测站r到BDS系统中的卫星m的伪距观测值的残差， 为无电离组合后的从观测站r到GLONASS系统中的卫星n的伪距观测值的残差； 为BDS系统中的卫星m到观测站r的单位矢量， 为GLONASS系统中的卫星n到观测站r的单位矢量，[Δx,Δy,Δz]为坐标改正数， 为BDS系统无电离层组合后的伪距观测值的残差及未模型化误差， 为GLONASS系统无电离层组合后的伪距观测值的残差及未模型化误差，未模型化误差包括对流层残余误差，电离层二阶项误差；此时方程(2)中的待估参数为待估参数包含有位置参数即坐标改正数[Δx,Δy,Δz]、接收机钟差参数即以北斗时为基准的接收机钟差 系统时差参数即BDS系统与GLONASS系统的系统时差 以及GLONASS系统的13个频间偏差参数，即IFB 1,0 ,…,IFB 13,0 ；S4、利用卡尔曼滤波对待估参数进行参数估计；利用卡尔曼滤波对待估参数进行参数估计，位置参数、BDS系统与GLONASS系统的系统时差 和GLONASS系统的频间偏差参数在一天的时间间隔内是保持不变的，时差解算的状态方程和过程噪声为式中，I 3×3 为3阶的单位矩阵， 和 分别为位置参数的过程噪声方差、接收机钟差参数的过程噪声方差以及GLONASS频间偏差参数的过程噪声方差；过程噪声的大小设置为伪距观测噪声的十分之一；GLONASS频间偏差参数的初值设置为0.23，GLONASS系统频间偏差参数在滤波的过程中会逐渐收敛到一个稳定的值。
2.根据权利要求1所述的顾及到GLONASS频间偏差的时差监测方法，其特征在于，S2中，预处理包括包括判断伪距观测值是否在正常范围内以及依据卫星健康标示判断伪距观测值所对应的卫星状态是否正常，在进行上述判断分析后将不合格的观测数据作为异常值进行剔除。
