GPS整周模糊度解算方法.doc

《GPS整周模糊度解算方法.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、GPS整周模糊度解算方法探讨一、为什么要解算GPS整周模糊度？整周模糊度的确定是载波相位测量中的关键问题，这是因为：（1）精确的、不足1周的相位观测值和修复周跳后的正确的整周计数只有与正确的整周模糊度配合使用才有意义。模糊度参数一旦出错，就将导致大量的卫地距出现系统性的粗差，从而严重损害定位的精度和可靠性。正确确定整周模糊度是获得高精度定位结果的必要条件。（2）在一般精度的GPS定位中，定位所需的时间实际上就是正确确定整周模糊度所需要的时间。快速确定整周模糊度对提高GPS定位的作业效率具有极其重要的作用；对开拓GPS定位技术的应用领域，将其推广应用到低等级控制测量和一般的工程测量等领域也具有极

2、其重要的作用。二、GPS整周模糊度解算方法1、LAMBDA法1993年荷兰Delft大学的Teunissen教授提出了最小二乘模糊度降相关平差法，简称LAMBDA法。该方法可缩小搜索范围，加快搜索过程，是目前快速静态定位中最成功的一种模糊度搜索方法。LAMBDA法的基本原理：（1）整数变换在LAMBDA法中，并不直接对整数模糊度参数进行搜索，而是先对初始解中的实数模糊度参数及其协因数阵进行整数变换：式中为整数变换矩阵。整数变换具有以下特点：当为整数时，变换后的参数也为整数；反之，当为整数时，经逆变换后所得的也为整数。整数变换并不是唯一的。我们希望整数变换后所得到的新参数之间的相关性能显著减小，

3、其协因数阵中的非对角线元素，模糊度参数的方差也能大幅度减小。注意，整数变换指的是具有上述特性的一种数学变换方法，但并非只能对整数进行变换。在LAMBDA法的正变换中，是在对实数模糊度进行变换。（1）搜索算法欲寻求经整数变换后的新参数的整数最小二乘解，实际上就是要寻找能满足下式的整数组合：由于上式无法直接求解，故一般都采用搜索算法从备选组中将满足上式的整数组合挑选出来。由于变换后的新参数的方差及参数间的互相关性均较前大大减小，故搜索工作将更为简便、迅速。求得最优的整数组合后再进行逆变换：变换后的参数满足下列公式：逆变换后求得的参数就是最初要寻找的最佳的整数模糊的向量。2、快速模糊度解算法1990

4、年E.Frei和G.Beutler提出了快速模糊度解算法（FARA—FastAmbiguityResolutionApproach）。（1）基本概念进行快速定位时虽然观测时间较短，但只要能正确确定整周模糊度，仍能获得相当好的结果，因此快速定位的关键在于快速确定整周模糊度。我们知道，用短时间的观测资料所建立的方程的状态很差，方程几乎是线性相关的，在这种情况下所求得的实数模糊度参数的中误差必然很大。整数模糊度向量的备选组中只有一组整数模糊度组合是完全正确的，如果我们能将这组模糊度组合挑选出来取用，那么快速定位就能取得很好的结果。（2）搜索原理将备选组中的整数模糊度组合一一代入法方程中进行计算，那么

5、能使观测残差的平方和为最小的这组整数模糊度组合就是最终的正确解。只有当所有的整周模糊度皆取正确值时，观测值得残差才会与载波相位测量的正确精度相对应，其他组的代入由于卫地距出现粗差，从而使观测值残差的平方和迅速增大。在未知参数必须为整数的情况下求最小二乘解的方法称为整数最小二乘法。最小二乘解的一种形式：式中：为初始解中求得的一组实数模糊度解；为这组实数模糊度的协因数阵；为整数模糊度组合。能满足式子的这组整数模糊度就是我们所寻求的最优的整数模糊度组合。（3）FARA法FARA法的实质就是在上式进行计算前，先对备选组进行数理统计检验，把大量的显然不合理备选组先剔除掉，以减少计算工作量。统计检验的标准

6、是：任意两个整数模糊度参数和之差是否位于这两个模糊度差值的置信区间内。FARA法充分利用了初始解协因数阵中的非对角线元素所提供的模糊度间的相互关信息，对参数作进一步的数理统计检验。通过统计检验，可以把大量的不合理的整数组合迅速予以剔除。然后求出相应的（为观测值总数）及单位权中误差（为未知数的个数）。从原则上讲，能使取最小值的那组整数模糊度组合就是我们所寻求的最优的整数模糊度组合。（4）确认最优解需进行的三项统计检验①整数解与初始解所求得得基线向量的一致性检验。②整数解和初始解的单位权中误差的一致性检验。③整数解中最小单位权中误差与次最小单位权中误差间的显著性检验。3、GPS变形监测中整周模糊度

7、解算的新方法利用变形监测网中监测点坐标已知的特点，提出了一种新的解算整周模糊度的方法——DC（directcalculation）算法。该方法不需要组成和解算法方程，更不需要搜索和确认，而是直接计算整周模糊度。在GPS变形监测中，可采用单历元计算整周模糊度，单历元解是根据GPS单历元观测值解算基线向量，从而获得变形信息。（1）DC算法的原理如图所示，设号卫星为参考卫星，则可以得到单差观测方程为：（