GPS整周模糊度解算方法

GPS整周模糊度解算方法探讨一、为什么要解算GPS整周模糊度？整周模糊度的确定是载波相位测量中的关键问题，这是因为：（1）精确的、不足1周的相位观测值rF和修复周跳后的正确的整周计数Int只有与正确的整周模糊度配合使用才有意义。模糊度参数一旦出错，就将导致大量的卫地距出现系统性的粗差，从而严重损害定位的精度和可靠性。正确确定整周模糊度N是获得高精度定位结果的必要条件。（2）在一般精度的GPS定位中，定位所需的时间实际上就是正确确定整周模糊度所需要的时间。快速确定整周模糊度对提高GPS定位的作业效率具有极其重要的作用；对开拓GPS定位技术的应用领域，将其推广应用到低等级控制测量和一般的工程测量等领域也具有极其重要的作用。二、GPS整周模糊度解算方法1、LAMBDA法1993年荷兰Delft大学的Teunissen教授提出了最小二乘模糊度降相关平差法，简称LAMBDA法。该方法可缩小搜索范围，加快搜索过程，是目前快速静态定位中最成功的一种模糊度搜索方法。LAMBDA法的基本原理：（1）整数变换在LAMBDA法中，并不直接对整数模糊度参数N进行搜索，而是先对初始解中的实数模糊度参数nNNNN,......,,21及其协因数阵NQ进行整数变换：NZzTZQZQNTz式中Z为整数变换矩阵。整数变换具有以下特点：当N为整数时，变换后的参数z也为整数；反之，当z为整数时，经逆变换后所得的zZNT1也为整数。整数变换并不是唯一的。我们希望整数变换后所得到的新参数nzzzz,......,,21之间的相关性能显著减小，其协因数阵zQ中的非对角线元素5.0，模糊度参数的方差也能大幅度减小。注意，整数变换指的是具有上述特性的一种数学变换方法，但并非只能对整数进行变换。在LAMBDA法的正变换中，是在对实数模糊度进行变换。（2）搜索算法欲寻求经整数变换后的新参数z的整数最小二乘解，实际上就是要寻找能满足下式的整数组合nzzzz,......,,21：min1zzQzzzT由于上式无法直接求解，故一般都采用搜索算法从备选组中将满足上式的整数组合z挑选出来。由于变换后的新参数的方差及参数间的互相关性均较前大大减小，故搜索工作将更为简便、迅速。求得最优的整数组合z后再进行逆变换：zZNT1变换后的参数N满足下列公式：min1NNQNNNT逆变换后求得的参数N就是最初要寻找的最佳的整数模糊的向量。2、快速模糊度解算法1990年E.Frei和G.Beutler提出了快速模糊度解算法（FARA—FastAmbiguityResolutionApproach）。（1）基本概念进行快速定位时虽然观测时间较短，但只要能正确确定整周模糊度，仍能获得相当好的结果，因此快速定位的关键在于快速确定整周模糊度。我们知道，用短时间的观测资料所建立的方程的状态很差，方程几乎是线性相关的，在这种情况下所求得的实数模糊度参数的中误差必然很大。整数模糊度向量N的备选组中只有一组整数模糊度组合是完全正确的，如果我们能将这组模糊度组合挑选出来取用，那么快速定位就能取得很好的结果。（2）搜索原理将备选组中的整数模糊度组合一一代入法方程中进行计算，那么能使观测残差的平方和为最小的这组整数模糊度组合就是最终的正确解。只有当所有的整周模糊度皆取正确值时，观测值得残差才会与载波相位测量的正确精度相对应，其他组的代入由于卫地距出现粗差，从而使观测值残差的平方和迅速增大。在未知参数必须为整数的情况下求最小二乘解的方法称为整数最小二乘法。最小二乘解的一种形式：min1NNQNNNT式中：N为初始解中求得的一组实数模糊度解；1NQ为这组实数模糊度的协因数阵；N为整数模糊度组合。能满足式子的这组整数模糊度就是我们所寻求的最优的整数模糊度组合。（3）FARA法FARA法的实质就是在上式进行计算前，先对备选组进行数理统计检验，把大量的显然不合理备选组先剔除掉，以减少计算工作量。统计检验的标准是：任意两个整数模糊度参数iN和jN之差ijN是否位于这两个模糊度差值的置信区间内。FARA法充分利用了初始解协因数阵中的非对角线元素所提供的模糊度间的相互关信息，对参数作进一步的数理统计检验。通过统计检验，可以把大量的不合理的整数组合迅速予以剔除。然后求出相应的miiV12（m为观测值总数）及单位权中误差miinmV12/（n为未知数的个数）。从原则上讲，能使取最小值的那组整数模糊度组合就是我们所寻求的最优的整数模糊度组合。（4）确认最优解需进行的三项统计检验①整数解与初始解所求得得基线向量的一致性检验。②整数解和初始解的单位权中误差的一致性检验。③整数解中最小单位权中误差与次最小单位权中误差次间的显著性检验。3、GPS变形监测中整周模糊度解算的新方法利用变形监测网中监测点坐标已知的特点，提出了一种新的解算整周模糊度的方法——DC（directcalculation）算法。该方法不需要组成和解算法方程，更不需要搜索和确认，而是直接计算整周模糊度。在GPS变形监测中，可采用单历元计算整周模糊度，单历元解是根据GPS单历元观测值解算基线向量，从而获得变形信息。（1）DC算法的原理如图所示，设j号卫星为参考卫星，则可以得到单差观测方程为：RRttjjjVVcN21212121（1）RRttkkkVVcN21212121（2）由式（1）、式（2）可得双差观测方程为：jkjkjkNN212121212121（3）由式（3）可以解出整周模糊度为：jkjkjkNN212121212121/可见，当已知卫星的位置和监测点的位置时，就可以直接计算出整周模糊度，上式就是解算整周模糊度的DC算法。（2）监测点的变形量对整周模糊度解算的影响由图所示，卫星到监测点间的距离为：222pspspsZZYYXX式中，sssZYX,,为卫星s的坐标，pppZYX,,为监测点p的坐标。载波相位的实际观测值tji与卫星和地球的距离的关系为：Ntji于是有：tNji由上式对求微分得：pppdzzdyydxxddN000应用协方差传播律得：22022202220222zyxpNzyx取zyx，得：222022222xxpNzyx若要求N0.5周，及2/1Lp，因为1L=0.1903m，所以有：mzyx09515.0当监测点的位移mdxzyx1648.032222，它对整周模糊度的影响小于等于半周。