专利名称：用于处理gps变形监测数据的方法
技术领域：
本发明涉及一种用于处理数据的方法，特别涉及一种用于处理GPS变形监测数据的方法。
背景技术：
随着GPS技术的发展，GPS已越来越多地用于变形监测。但在变形监测环境中，通常多路径或衍射信号较强，导致周跳发生频率较高，而常用的最小二乘搜索法(LS)、快速模糊度搜索法(FARA)、降低模糊度之间相关性的最小二乘非相关平差法(LAMBDA)、Cholesky分解搜索法等都需要利用一段时间内的多个历元的观测数据信息，且这多个历元的观测数据不能出现周跳，否则将很难确定整周模糊度。因此在变形监测应用中，GPS数据的处理一般采用单历元解算方法以避开周跳探测，即根据监测点位置变化不大的特点，利用变形量约束进行单历元模糊度求解。同样由于监测环境的限制，可视GPS卫星通常较少，导致卫星分布图形结构强度较弱，加上多路径效应的影响，使得所获取的结果精度不高，在变形序列中会出现跳变的情况，从而难以根据监测结果进行预警。针对以上现有技术中的用于处理GPS变形监测数据的方法的不足，研究精度和可靠性都得到提升的新的GPS数据处理方法，对变形监测技术的发展和其在工程中的实际应用具有积极的意义。

发明内容
针对现有技术中的用于处理GPS变形监测数据的方法无法回避跳变，导致难以根据监测结果精确预警的不足，本发明提出了一种新颖的用于处理GPS变形监测数据的方法。本发明提出了一种用于处理GPS变形监测数据的方法，包括:第一步:获取基准站和监测站当前历元观测数据及广播星历数据，利用相位观测值组成双差观测方程；第二步:设置两历元间的最大变形量；第三步:采用Cholesky分解构建模糊度搜索空间；第四步:搜索模糊度，并且使方差比大于3 ;第五步:获取监测点的单历元整周模糊度固定解坐标X、Y、Z和其协方差矩阵；第六步:利用平均间隙法对当前历元和前一历元间的变形值进行变形检验；第七步:经第六步的检验后，如果两历元间未发生位移，则将两个历元观测方程叠加的抗差最小二乘静态解作为当前历元的结果；若第六步的检验结果为两历元间有位移发生，则根据抗差序贯平差获取当前历元的结果。在一个实施例中，在第三步中:
权利要求
1.一种用于处理GPS变形监测数据的方法，包括: 第一步:获取基准站和监测站当前历元观测数据及广播星历数据，利用相位观测值组成双差观测方程； 第二步:设置两历元间的最大变形量； 第三步:采用Cholesky分解构建模糊度搜索空间； 第四步:搜索模糊度，并且使方差比大于3 ; 第五步:获取监测点的单历元整周模糊度固定解坐标X、Y、Z和其协方差矩阵； 第六步:利用平均间隙法对当前历元和前一历元间的变形值进行变形检验； 第七步:经第六步的检验后，如果两历元间未发生位移，则将两个历元观测方程叠加的抗差最小二乘静态解作为当前历元的结果；若第六步的检验结果为两历元间有位移发生，则根据抗差序贯平差获取当前历元的结果。
2.根据权利要求1所述的方法，其特征在于，在第三步中:
3.根据权利要求1或2所述的方法，其特征在于，在第六步中还包括如下子步骤: 第一子步骤:通过两历元观测值改正数获取经验方差:
4.根据权利要求1到3中任一项所述的方法，其特征在于，在第四步中，利用最小二乘法来搜索模糊度。
5.根据权利要求1到4中任一项所述的方法，其特征在于，在第五步中，利用最小二乘法来获取监测点的单历元整周模糊度固定解坐标X、Y、Z和其协方差矩阵。
全文摘要
本发明提出了一种用于处理GPS变形监测数据的方法，包括第一步获取基准站和监测站当前历元观测数据及广播星历数据，利用相位观测值组成双差观测方程；第二步设置两历元间的最大变形量；第三步采用Cholesky分解构建模糊度搜索空间；第四步搜索模糊度，并且使方差比大于3；第五步获取监测点的单历元整周模糊度固定解坐标X、Y、Z和其协方差矩阵；第六步利用平均间隙法对当前历元和前一历元间的变形值进行变形检验；第七步经第六步的检验后，将两个历元观测方程叠加的抗差最小二乘静态解作为当前历元的结果或者根据抗差序贯平差获取当前历元的结果。本方法很准确。
文档编号G01B15/06GK103148813SQ201310037829
公开日2013年6月12日 申请日期2013年1月31日 优先权日2013年1月31日
发明者朱自强, 戴吾蛟, 杨文韬, 曾凡河, 周勇, 谢依胜 申请人:湖南致力地质资源环境工程科技有限公司