[发明专利]零/短基线GLONASS相位通道间偏差校准方法及系统

权利要求书

1.零/短基线GLONASS相位通道间偏差校准方法，其特征是，包括：

S100读取单历元的GLONASS观测数据，计算站星几何距离，并组建双差观测量；

S200对频间偏差IFB进行粗校准，本步骤进一步包括：

S210判断当前历元对应的双差相位观测量中是否存在特定通道数之差的双差相位观测量，若存在，执行步骤S220；否则，返回步骤S100，读取下一历元的GLONASS观测数据；

S220选取特定通道数之差的双差相位观测量，直接取整固定双差模糊度，获得双差整周模糊度；其中，当GLONASS观测数据是单频数据时，特定通道数之差为±1；当GLONASS观测数据是双频数据时，特定通道数之差为±1、±2和±3；

当GLONASS观测数据是单频数据时，双差整周模糊度为L1双差整周模糊度

其中，[·]表示四舍五入取整符号；表示在测站m、n和卫星p、q间的L1双差整周模糊度；表示在测站m、n和卫星p、q间以周为单位的L1双差相位观测量；和分别表示卫星q和p在测站m、n间的单差站星几何距离；和分别表示卫星q和p的L1波长；k^pq表示卫星p、q间的通道数之差；a₁表示频分多址中第1个频段通道数的系数；dt_mn表示测站m、n间接收机的相对钟差；

当GLONASS观测数据是双频数据时，双差整周模糊度为宽巷双差整周模糊度

其中，[·]表示四舍五入取整符号；表示测站m、n和卫星p、q间的宽巷双差整周模糊度；表示测站m、n和卫星p、q间的宽巷双差相位观测量；和分别表示卫星p和q在测站m、n间的单差站星几何距离；和分别表示卫星q和p的宽巷波长；k^pq表示卫星p、q间的通道数之差；dt_mn表示测站m、n间接收机的相对钟差；a₁和a₂分别表示频分多址中第1个和第2个频段通道数的系数；

S230将固定的双差整周模糊度代回双差观测方程，反算相位ICB值；

子步骤S230中，相位ICB值的反算公式如下：

其中，ICB^pq表示卫星p、q间的相位通道间偏差；表示在测站m、n和卫星p、q间的双差站星几何距离项；j表示不同的解算模式，可取1或WL；

对于单频数据，j取1，则和分别表示卫星p和q的L1波长；和分别表示L1波段下卫星q和p在测站m、n间的站间单差相位观测值；表示在测站m、n和卫星p、q间的L1双差整周模糊度；表示由伪距计算的参考卫星p在测站m、n间的L1站间单差模糊度；

对于双频数据，j取WL，则和分别表示卫星q和p的宽巷波长；和分别表示卫星q和p在测站m、n间的宽巷站间单差相位观测值；表示测站m、n和卫星p、q间的宽巷双差整周模糊度；表示由伪距计算的参考卫星p在测站m、n间的宽巷站间单差模糊度；

S240对步骤S230反算的相位ICB值进行线性拟合，获得IFB粗校准值；线性拟合公式为IFB₀＝ICB^pq/k^pq，其中，IFB₀为IFB粗校准值，ICB^pq为步骤S230反算的相位ICB值；

S300重复S100～S200，直至在所有历元下均完成IFB粗校准，取IFB粗校准值的平均值，即IFB校准初值IFB₀；

S400频间偏差IFB进行精校准，本步骤进一步包括：

S410读取一单历元的GLONASS观测数据，并选取该历元所有的双差相位观测量；

S420利用IFB校准初值改正子步骤S410所选取的双差相位观测量中的相位ICB，并取整固定双差整周模糊度；

子步骤S420中，取整固定双差整周模糊度的公式如下：

其中，f_i⁰表示频分多址中第i个频率的常数项，IFB₀表示粗校准IFB值，c表示光速；i表示信号波段编号，对单频数据，i＝1，上式表示L1双差整周模糊度的固定公式，相关的变量为L1频率所对应变量；对双频数据，i＝1,2，上式分别表示L1和L2双差整周模糊度的固定公式，相关的变量为L1和L2频率对应变量；

S430将子步骤S420固定的双差整周模糊度代回双差观测方程，反算相位ICB值；

子步骤S430中，相位ICB值的反算公式如下：

S440对子步骤S430反算的相位ICB值进行线性拟合，获得IFB精校准值；线性拟合公式为IFB₁＝ICB^pq/k^pq，其中，IFB₁为IFB精校准值：ICB^pq为步骤S430反算的相位ICB值；

S500重复步骤S400，直至对所有历元完成IFB精校准。