基本观测方程在GPS长基线测量的一种应用.docVIP

基本观测方程在GPS长基线测量的一种应用 Hiroshi Isshiki, Pegasusnet公司, 静冈,日本 E-mail:isshiki@dab.hi-ho.ne.jp 作者简介 Dr.H.Isshiki是Pegasusnet公司日本静冈的首席科学家(Isshiki et al.(2000a,b),(2001,a,b))。然而，如果整周模糊度刚好需要，那就必须要完整的解算出来。 为了解决整周不定性问题，对于短基线距离测量来说，人们已经提出很多解决方法，并且这些方法得到了成功的应用。本文提出了一种方法，这种方法利用了四个基本观测方程（即用一台双频接收机得到的两个测码伪距观测方程和两个测相伪距观测方程）。根据这种方法，当接收机钟差、卫星钟差和卫星轨道误差都充分小时，这两个初始相位整周模糊度、电离层延迟、距离真值的总和和电离层延迟在理论上完全可以得到解决。如果卫星和接收机钟差可以在双差模型中得到消除，那么在双差模型中也可以解决两个初始相位整周模糊度、电离层延迟、距离真值的总和和电离层延迟问题。 行驶在海洋中距离陆地几千公里远的移动物体的精密定位是很重要的，但是定位中的问题并没有完全得到解决。单点定位（不需要参考站）可能是一种解决这些问题的方法。上述方法运用到这个方向可能非常有用。基于这个目的，我们必须提出一种用来消除卫星钟差、卫星轨道误差和对流层延迟的实用方法。 基本观测方程 定义 和分别为历元，测码伪距，测相伪距，真实距离，载波L1的电离层延迟和对流层延迟，上标i和下标分别表示卫星和接收机，下标对应载波L1,L2。基本观测方程如下所示： 用 * 标记的量表示观测值，没有用*标记的量表示真值。 和 是噪声部分。一般来说，前者比后者大十倍。 和 的表达式如下： 其中 噪声在测码伪距中，是由于卫星产生的 噪声在测码伪距中，是由于接收机产生的 噪声（不同于电离层和对流层延迟）在测码伪距中，是由于GPS信号在卫星和接收机之间的传播路径上产生的 杂噪声在测码伪距中 噪声在测相伪距中，是由于卫星产生的 噪声在测相伪距中，是由于接收机产生的 噪声（不同于电离层和对流层延迟在测相伪距中，是由于GPS信在卫星和接收机之间的传播路径上产生的 杂噪声在测相伪距中 对于卫星和接收机每一种组合，都可以组成五个方程包含五个未知量 和。所以，这个方程组在双频观测情况下是不可解。然而，如果把和合并为，则四个方程中包含四个未知量和。所以这个方程中肯定是可以解的，但是，和不能单独解出，而只能解出的值。假设三频观测，未知量就会增加到五个， 分别为： 和一个新的未知量，其中是在中的初始相位整周模糊度。此时，方程增加到六个，但是仍然没法分别解出和。下面，我们将继续讨论双频观测。 定义 和 为新的观测量，通过忽略方程组（1）的噪声影响而获得的。 公式如下: 方程组（3）解算结果如下： 把方程（4c）和（4d）分别求和并求差，初始相位整周模糊度和用如下公式表达： 把和代入，是的65倍。因此，的精度要比的精度高的多。 观测量的线性组合做如下总结。定义和为窄巷观测值 窄巷观测值的观测方程可化为如下形式： 下面，定义和为宽巷观测值的虚拟频率，波长和初始相位整周模糊度： 宽巷观测值的观测方程可化为如下形式： 最后，消去电离层误差的观测方程可化为如下形式： 在这个方程组合中包含了电离层延迟的影响。 观测值拟合 假设方程组（1）的误差相互独立,由此可以知道：由方程组（1）的平均值导出如下方程： 其中表示在某一观测历元t的观测值的平均值。在下面的算例中，从开始观测到观测结束的时间跨度平均值。 根据方程组（11），可推导出 ，的表达式如下： 方程组（12）可以由方程（4c）和方程（4d）取平均值获得。也就是说如果时间跨度充分长的情况下，噪声误差可以得到消除， ，的真值也可以得到。 如果把方程（1c）和方程（1d）代入方程（12a）和方程（12b），则可以推导出，的表达式如下： 把方程（13a）和方程（13b）右边最后部分忽略不计，则剩下部分，表达式为： 根据方程组（14），用和可以计算出和。这就意味着前者与后者的精度相同了。因此，方程组（14）能显著提高和的精度。 但是的求解需要的估值。 遗憾的是，是未知的。如果表示表达式部分的最大值，则可以近似等于： 其中为卫星和接收机之间的高度角。 当假设干分量部分为2.1m和湿分量部分为0.4m，的近似值为2.5m。如果使用水汽辐射计会得到一个更为精确的值。 算例分析 宽巷观测值的初始相位整周模糊度基线长度在和之间双差模型的收敛性如图1.红色线条所示。蓝色线条表示基于方程(4c,d)和（12a,b）计算出的结果。基线长度有数百米。基于方程（12a，b）计算出的结果有很好的收敛性。当基线长度很长时，卫星钟差、卫星轨道误差、卫星信