GPS变形监测-原与方法(王新洲).ppt

武汉大学灾害监测与防治研究中心 五、GPS天线阵列变形系统 4、 GPS天线阵列变形监测系统的核心技术 通信及远程控制 远程管理 本系统采用网络数据库管理技术，实现系统的远程操作、数据的远程处理。 五、GPS天线阵列变形系统 4、 GPS天线阵列变形监测系统的核心技术 相关理论与算法 DC算法 随着研究的深入，本人提出了一种全新的单历元解算法。该算法的核心是利用变形监测的特点和已有信息，提出了一种既准确又快速的整周模糊度解算方法——DC算法。该方法不需要组成和解算法方程，更不需要搜索和确认，而是直接计算整周模糊度。 DC算法的基本原理 如图，双差后有： 两式求差，得 ： 令： 五、GPS天线阵列变形系统 4、 GPS天线阵列变形监测系统的核心技术 相关理论与算法 DC算法 则可以解出双差整周模糊度： 当已知卫星的位置和监测点的位置时，由该式可直接计算出的整周模糊度的浮点解。 如此计算时，监测点的变形量对整周模糊度的影响有多大，在变形多大的时候才可以直接采用这种方法来解算呢？ 经推导，应用协方差传播律得： 若要求 。 五、GPS天线阵列变形系统 4、 GPS天线阵列变形监测系统的核心技术 相关理论与算法 DC算法 即当 =0.1648m时，位移对整周模糊度的影响小于等于半周。这种利用监测点的位移为约束条件直接解算单历元变形量的方法就称为DC算法。 位移小于0.1648m时整周模糊度的解算结果 采用TRIMBLE 5700接收机采集实测GPS数据,采样率为15秒，采集时间为2个小时，组成5个卫星对，以4号星为参考星，分别按DC算法和一般的单历元解算方法解算了卫星对4-7、4-8、4-11、4-20、4-28的L1波段50个历元的整周模糊度 。 计算结果比较 10.80 12.00 57.66 57.75 -9.32 -8.18 67.50 68.76 49.37 50.07 最小值 12.80 12.10 57.84 57.82 -6.04 -8.04 71.20 69.00 50.58 50.17 最大值 一般的单历元算法 DC算法 一般的单历元算法 DC算法 一般的单历元算法 DC算法 一般的单历元算法 DC算法 一般的单历元算法 DC算法 4-28 （准确周数12） 4-20 （准确周数58） 4-11 （准确周数-8） 4-8 （准确周数69） 4-7 （准确周数50） 卫星对 五、GPS天线阵列变形系统 4、 GPS天线阵列变形监测系统的核心技术 相关理论与算法 DC算法 。 五、GPS天线阵列变形系统 4、 GPS天线阵列变形监测系统的核心技术 相关理论与算法 DC算法 左图是卫星对4-28的双差模糊度实数解，正确的周数为12周，当位移从0cm变化到16cm时，双差整周模糊度的变化也越来越大，但与真实的双差整周模糊度之间的差值始终小于0.5周。 五、GPS天线阵列变形系统 4、 GPS天线阵列变形监测系统的核心技术 相关理论与算法 DC算法 当位移从16cm变化到60cm时，当位移为60cm时，浮点解与真实的双差整周模糊度的差值达到2周。 位移量大于0.1648m时整周模糊度的解算方法 由上图知，当监测点的位移大于0.1648m时，直接解算L1载波的整周模糊度是不准确的，而在实际工作中，存在位移量大于0.1648m的情况，如悬索桥的位移往往大于0.5m。在这种情况下，我们又导出了组合求解的方法。 当位移值大到0.7m时，如果还希望对整周模糊度的影响小于等于半周，并应用DC算法计算整周模糊度，则载波相位的波长应为： m 因为不同的载波相位组合观测值,可以认为是不同波长和精度的载波相位观测值。设和为一理想卫星发射的双频相位观测值，则它们组合的观测值： 适当选取n和m，可以保证组合载波的波长大于等于0.8083m。于是可应用DC算法计算组合载波的整周模糊度 ，使位移值对组合载波的整周模糊度的影响小于等于半周。 如此解出的组合载波的整周模糊度虽然比较准确，但不是我们所需要的最终结果。我们所需要的最终结果是 和 的整周模糊度 和 。 为了解决这个问题，我们再找一个波长大于0.8083m以上的载波相位组合观测值 ，并应用DC算法解得此组合的整周模糊度