>[0101] 步骤2-1-2、利用步骤3计算得到的IGSO卫星和MEO卫星的双差电离层延迟误差残 差和双差对流层延迟误差残差、双差卫星轨道误差残差，修正IGS0卫星和ΜΕ0卫星的B1频率 双差载波相位整周模糊度与B3频率双差载波相位整周模糊度；
[0102] 利用IGS0卫星和ΜΕ0卫星的双差电离层延迟误差残差AV/cwf和双差对流层延迟误 差残差A V0rbs、双差卫星轨道误差残差Δ VTrops，修正公式（12)、（13)，得到：
[0103]
[0105] 将公式(14)、（15)相减得到：
[0106]
[0107]步骤2-1-3、将修正后的IGS0卫星和ΜΕ0卫星的B1频率双差载波相位整周模糊度与 B3频率双差载波相位整周模糊度作差zWiVf - AWVf取整，确定B1频率双差载波相位与B3 频率双差载波相位的宽巷整周模糊度，或将计算出的多个历元的宽巷整周模糊度的数值取 平均值再取整确定B1频率双差载波相位与B3频率双差载波相位的宽巷整周模糊度。
[0108] 步骤2-2、根据IGS0卫星和ΜΕ0卫星的B1频率双差载波相位与B3频率双差载波相位 的宽巷整周模糊度，计算GE0卫星的B1频率双差载波相位与B3频率双差载波相位的宽巷整 周模糊度；
[0109] 步骤2-2-1、将IGS0卫星和ΜΕ0卫星的B1频率双差载波相位与B3频率双差载波相位 的宽巷整周模糊度带入包含位置参数的定位观测方程组；
[0110] 包含位置参数的定位观测方程组如公式(17)所示：
[0111]
[0112]
[0113]
[0114] 其中，λ13表示B1频率载波相位与B3频率载波相位的宽巷模糊度波长，Δ VN13 = Δ ▽Nr Δ VN3为BDS卫星的B1频率双差载波相位的宽巷整周模糊度与B3频率双差载波相位的 宽巷整周模糊度，X为BDS测站的位置参数，h为位置参数的系数，ξ为宽巷观测值电离层延迟 误差、对流层延迟误差和卫星轨道误差的双差残差综合影响，上标I表示IGSO卫星，上标Μ表 示ΜΕΟ卫星，上标G表示GEO卫星；
[0115] 步骤2-2-2、通过IGS0卫星和ΜΕ0卫星以及GE0卫星的包含位置参数的定位观测方 程组，计算确定GE0卫星的Β1频率双差载波相位与Β3频率双差载波相位的宽巷整周模糊度。
[0116] 通过公式（17)解算出Β1频率双差载波相位与Β3频率双差载波相位的宽巷整周模 糊度数值，并对其进行取整确定Β1频率双差载波相位与Β3频率双差载波相位的宽巷整周模 糊度；
[0117] IGS0卫星和ΜΕ0卫星的Β1、Β3频率的宽巷模糊度已经固定，即为已知值;对于公式 (17)的解算，在双差组合的情况下忽略ξ的影响，则有：
[0125] L1、LM、LG是分别利用载波相位观测值、初始几何距离以及固定的Β1频率、Β3频率载 波相位宽巷整周模糊度计算得到的已知项。
[0126] 对公式（18)的BDS测站的位置参数X和GE0卫星的B1频率双差载波相位与B3频率双 差载波相位的宽巷整周模糊度▲▽^进行解算，计算得到 BDS测站的位置参数X和GE0卫星 的B1频率双差载波相位与B3频率双差载波相位的宽巷整周模糊度参数AViVg，然后通过取 整确定GE0卫星的B1频率双差载波相位与B3频率双差载波相位宽巷整周模糊度AWV|的整 数数值。
[0127] 步骤3、利用IGS0卫星、ΜΕ0卫星以及GE0卫星的B1频率双差载波相位与B3频率双差 载波相位的宽巷整周模糊度，计算确定IGS0卫星、ΜΕ0卫星以及GE0卫星的B1频率双差载波 相位整周模糊度、B2频率双差载波相位整周模糊度、B3频率双差载波相位整周模糊度，进而 计算得到BDS测站的坐标。
[0128] 如图4所示，具体步骤如下：
[0129] 步骤3-1、将IGS0卫星、ΜΕ0卫星以及GE0卫星的B1频率双差载波相位与B3频率双差 载波相位的宽巷整周模糊度，带入包含位置参数以及IGS0卫星、ΜΕ0卫星和GE0卫星的B1频 率双差载波相位整周模糊度参数的定位观测方程组，计算确定IGS0卫星、ΜΕ0卫星以及GE0 卫星的B1频率双差载波相位整周模糊度；
[0130] 包含位置参数以及IGS0卫星、ΜΕ0卫星和GE0卫星的B1频率双差载波相位整周模糊 度参数的观测方程组如公式(19)所示：
[0131]
[0132]
[0133] (19)
[0134]
[0135]
[0136] 其中，人:为。.19203mm，表示B1频率载波相位模糊度的波长，〖:为价频率载波相位观 测值电离层延迟误差、对流层延迟误差和卫星轨道误差的双差残差综合影响，上标I表示 IGS0卫星，上标Μ表示ΜΕ0卫星，上标G表示GE0卫星；
[0137] BDS卫星的Β1频率双差载波相位与Β3频率双差载波相位的宽巷整周模糊度已经固 定，即为已知值，对于公式(19)的解算，在双差组合的情况下忽略ξΑ:的影响，则有：
[0151] 对公式(20)的位置参数X和Δν^Υ/、AWVf、AVTVf进行解算，计算得到位置参数 X和BDS卫星的B1频率双差载波相位整周模糊度参数AW/、麻沒;' Δν<，然后通过取 整确定BDS卫星的Β1频率双差载波相位整周模糊度AViV/、AViVf、AVAf的整数数值；
[0152] 步骤3-2、将IGS0卫星和ME0卫星以及GEO卫星的B1频率双差载波相位整周模糊度， 带入包含位置参数以及IGS0卫星、ΜΕ0卫星和GE0卫星的B2频率双差载波相位整周模糊度参 数的观测方程组，计算确定IGS0卫星、ΜΕ0卫星以及GE0卫星的B2频率双差载波相位整周模 糊度；
[0153] 包含位置参数以及IGS0卫星、ΜΕ0卫星和GE0卫星的B2频率双差载波相位整周模糊 度参数的观测方程组如公式(21)所示：
[0154]
[0155]
[0156] (21)
[0157]
[0158]
[0159] 其中，λΑ〇.24834πι，λ2表示B2频率载波相位模糊度的波长，ξ;*Β2频率载波相位 观测值电离层延迟误差、对流层延迟误差和卫星轨道误差的双差残差综合影响，其他符号 与公式（19)相同，上标I表示IGS0卫星，上标Μ表示ΜΕ0卫星，上标G表示GE0卫星；
[0160] BDS卫星的Β1频率双差载波相位整周模糊度已经固定，即为已知值，对于公式(21) 的解算，在双差组合的情况下忽略的影响，则有：
[0173] 对公式（22)的参数X和AW〗、AVTVf、AViVf进行解算，计算得到位置参数X和 BDS卫星的B2频率双差载波相位整周模糊度参数Δ▽况丨、AViVf、AViVf，然后取整确定B2 频率双差载波相位整周模糊度参数、AVAf、AWVf的整数数值；
[0174] 步骤3-3、将IGS0卫星、ΜΕ0卫星以及GEO卫星的B1频率双差载波相位整周模糊度， 带入包含位置参数以及IGS0卫星、ΜΕ0卫星和GE0卫星的B3频率双差载波相位整周模糊度参 数的观测方程组，计算确定IGS0卫星、ΜΕ0卫星以及GE0卫星的B3频率双差载波相位整周模 糊度。
[0175] 包含位置参数以及IGS0卫星、ΜΕ0卫星和GE0卫星的B3频率双差载波相位整周模糊 度参数的观测方程组如公式(23)所示：
[0181] 其中，λ3为0.23633m，表示B3频率载波相位模糊度的波长，ξ3为B3频率载波相位观 测值电离层延迟误差、对流层延迟误差和卫星轨道误差的双差残差综合影响，其他符号与 公式(21)相同，上标I表示IGS0卫星，上标Μ表示ΜΕ0卫星，上标G表示GE0卫星；
[0182] BDS