卫星的Β1频率双差载波相位整周模糊度已经固定，即为已知值，对于公式(23) 的解算，在双差组合的情况下忽略的影响，则有：
[0195] 对公式(24)的参数X和丨、进行解算，计算得到位置参数X和 BDS卫星的Β3频率双差载波相位整周模糊度参数:Δ\τ/ν3;、： AWf、AVAf，然后取整确定 B3频率双差载波相位整周模糊度参数Δντν37、ΔΥ/Vf、AVTVf的整数数值；
[0196] 利用BDS卫星三频双差载波相位整周模糊度和三频双差载波相位观测值数据计算 BDS测站坐标值。图5(a)~（c)分别为计算的不同方向上的BDS测站坐标计算值与已知精确 坐标的差值时间序列，其中（a)是N方向（北方向）BDS测站坐标计算值与已知精确坐标差值 的时间序列；（b)是E方向（东方向）BDS测站坐标计算值与已知精确坐标差值的时间序列； (c)是U方向（高程方向）BDS测站坐标计算值与已知精确坐标差值的时间序列。N、E、U三个方 向的定位结果差值都处于厘米级，其精度分别为N方向0.0154米，E方向0.0143米，U方向 0.0191米，可以表明K)S卫星三频双差载波相位整周模糊度已经准确固定，能够得到高精度 的测站坐标。
【主权项】
1. 一种K)S卫星载波相位整周模糊度快速解算方法，其特征在于，包括以下步骤： 步骤1、利用IGSO卫星和MEO卫星的B1频率双差伪距观测值和B3频率双差伪距观测值， 计算得到IGSO卫星和MEO卫星的B1频率双差载波相位与B3频率双差载波相位的宽巷整周模 糊度； 步骤2、利用IGSO卫星和MEO卫星的B1频率双差载波相位与B3频率双差载波相位的宽巷 整周模糊度，解算GE0卫星的B1频率双差载波相位与B3频率双差载波相位的宽巷整周模糊 度； 步骤3、利用IGSO卫星、ME0卫星以及GE0卫星的B1频率双差载波相位与B3频率双差载波 相位的宽巷整周模糊度，计算确定IGSO卫星、ME0卫星以及GE0卫星的B1频率双差载波相位 整周模糊度、B2频率双差载波相位整周模糊度、B3频率双差载波相位整周模糊度，进而计算 得到BDS测站的坐标。2. 根据权利要求1所述的BDS卫星载波相位整周模糊度快速解算方法，其特征在于，所 述步骤1按如下步骤进行： 步骤1-1、BDS卫星向BDS测站接收机播发三频非差载波相位观测值数据和三频非差伪 距观测值数据； 步骤1-2、将BDS接收机接收三频非差载波相位观测值和三频非差伪距观测值进行双差 组合，得到BDS卫星的三频双差载波相位观测值和三频双差伪距观测值； 步骤1-3、利用IGSO卫星和ME0卫星的B1频率双差伪距观测值和B3频率双差伪距观测 值，计算IGSO卫星和ME0卫星的双差电尚层延迟误差残差、双差对流层延迟误差残差、双差 卫星轨道误差残差。3. 根据权利要求2所述的BDS卫星载波相位整周模糊度快速解算方法，其特征在于，所 述步骤1-3按如下步骤进行： 步骤1-3-1、利用IGSO卫星和ME0卫星的B1频率双差伪距观测值和B3频率双差伪距观测 值，计算IGSO卫星和ME0卫星的B1频率双差电离层延迟误差残差、B3频率双差电离层延迟误 差残差； 步骤1-3-2、利用IGSO卫星和ME0卫星的B1频率双差伪距观测值和B3频率双差伪距观测 值，计算IGSO卫星和ME0卫星的双差对流层延迟误差残差、双差卫星轨道误差残差。4. 根据权利要求1所述的BDS卫星载波相位整周模糊度快速解算方法，其特征在于，所 述步骤2按如下步骤进行： 步骤2-1、利用IGSO卫星和ME0卫星的B1频率双差载波相位观测值和B3频率双差载波相 位观测值、IGSO卫星和ME0卫星的双差电离层延迟误差残差、双差对流层延迟误差残差以及 双差卫星轨道误差残差，计算IGSO卫星和ME0卫星的B1频率双差载波相位与B3频率双差载 波相位的宽巷整周模糊度； 步骤2-2、根据IGSO卫星和ME0卫星的B1频率双差载波相位与B3频率双差载波相位的宽 巷整周模糊度，计算GE0卫星的B1频率双差载波相位与B3频率双差载波相位的宽巷整周模 糊度。5. 根据权利要求4所述的BDS卫星载波相位整周模糊度快速解算方法，其特征在于，所 述步骤2-1按如下步骤进行： 步骤2-1-1、根据IGSO卫星和ME0卫星的双差载波相位观测方程，得到IGSO卫星和ME0卫 星的B1频率载波相位整周模糊度与B3频率载波相位整周模糊度； 步骤2_1_2、利用步骤3计算得到的IGSO卫星和MEO卫星的双差电尚层延迟误差残差和 双差对流层延迟误差残差、双差卫星轨道误差残差，修正IGSO卫星和MEO卫星的B1频率双差 载波相位整周模糊度与B3频率双差载波相位整周模糊度； 步骤2-1-3、将修正后的IGSO卫星和MEO卫星的B1频率双差载波相位整周模糊度与B3频 率双差载波相位整周模糊度作差取整，确定B1频率双差载波相位与B3频率双差载波相位的 宽巷整周模糊度，或将计算出的多个历元的宽巷整周模糊度的数值取平均值再取整确定B1 频率双差载波相位与B3频率双差载波相位的宽巷整周模糊度。6. 根据权利要求4所述的BDS卫星载波相位整周模糊度快速解算方法，其特征在于，所 述步骤2-2按如下步骤进行： 步骤2-2-1、将IGSO卫星和MEO卫星的B1频率双差载波相位与B3频率双差载波相位的宽 巷整周模糊度带入包含位置参数的定位观测方程组； 步骤2-2-2、通过IGSO卫星和ME0卫星以及GE0卫星的包含位置参数的定位观测方程组， 计算确定GE0卫星的B1频率双差载波相位与B3频率双差载波相位的宽巷整周模糊度。7. 根据权利要求4所述的BDS卫星载波相位整周模糊度快速解算方法，其特征在于，所 述步骤3按如下步骤进行： 步骤3-1、将IGSO卫星、ME0卫星以及GE0卫星的B1频率双差载波相位与B3频率双差载波 相位的宽巷整周模糊度，带入包含位置参数以及IGSO卫星、ME0卫星和GE0卫星的B1频率双 差载波相位整周模糊度参数的定位观测方程组，计算确定IGSO卫星、ME0卫星以及GE0卫星 的B1频率双差载波相位整周模糊度； 步骤3-2、将IGSO卫星和ME0卫星以及GE0卫星的B1频率双差载波相位整周模糊度，带入 包含位置参数以及IGSO卫星、ME0卫星和GE0卫星的B2频率双差载波相位整周模糊度参数的 观测方程组，计算确定IGSO卫星、ME0卫星以及GE0卫星的B2频率双差载波相位整周模糊度； 步骤3-3、将IGSO卫星、ME0卫星以及GE0卫星的B1频率双差载波相位整周模糊度，带入 包含位置参数以及IGSO卫星、ME0卫星和GE0卫星的B3频率双差载波相位整周模糊度参数的 观测方程组，计算确定IGSO卫星、ME0卫星以及GE0卫星的B3频率双差载波相位整周模糊度。
【专利摘要】一种BDS卫星载波相位整周模糊度快速解算方法，包括：计算IGSO卫星和MEO卫星的B1频率双差载波相位与B3频率双差载波相位的宽巷整周模糊度；利用IGSO卫星和MEO卫星的B1频率双差载波相位与B3频率双差载波相位的宽巷整周模糊度，解算GEO卫星的B1频率双差载波相位与B3频率双差载波相位的宽巷整周模糊度；利用各卫星的B1频率双差载波相位与B3频率双差载波相位的宽巷整周模糊度，计算各卫星的B1频率双差载波相位整周模糊度、B2频率双差载波相位整周模糊度、B3频率双差载波相位整周模糊度，进而计算得到BDS测站的坐标。本发明使用伪距观测值计算MEO卫星、IGSO卫星的载波相位宽巷整周模糊度，克服了GEO卫星伪距观测噪声较大的影响，解决了BDS高精度快速定位的关键问题。
【IPC分类】G01S19/44
【公开号】CN105510946
【申请号】CN201510999033
【发明人】祝会忠, 高猛, 徐宗秋, 徐爱功, 葛茂荣, 徐彦田
【申请人】辽宁工程技术大学
【公开日】2016年4月20日
【申请日】2015年12月28日