先发现周跳，并给出   让其浮动迭代           N凑整

           测站坐标作为近似值   计       算

   以上过程许多型号的接收机都可利用随机计算软件完成，因其计算的数学模型大同小异，不再逻列。它的360个历元全部进行预处理后，它的输出成果基线向量如附表五所列。

这里，我将1~360个历元的数据进行分历元计算，即如表一、二、三、四所列成果按1~90历元、90~180历元、180~270历元、270~360历元分历元按相对静态预处理模式计算。这些分历元计算的成果我们可以分析出具有以下的特点：

   （A）因为分历元的导航文件及其测站文件等都是按照整测段的数据计算，避免了按小段观测（比如，只观测30分钟，即1~90历元）时，整周模糊度不能确定成准确计算，并且单点定位成果甚至于由于数据量的影响而不能计算的情况。

   （B）如果需要第90历元（即开机后90历元×15秒/60秒=22.5分钟）时刻的监测点的最或然位置，则选取（90+45历元）~（90-45历元）历元的间隔区段进行计算。（因为GPS天线也随着变形监测点的位移而移动）。

   如果需要知道第180历元（即开机后180历元×15秒/60秒=45分钟）时刻的变形点最或然位置，则选取（180-45）~（180+45）历元间隔区段进行计算。

   所计算的第90历元和第180历元的数据成果列于附表6和附表7。

   从以上我们可以归纳出有以下的计算经验公式：

   A、计算，从开机时刻起算，第X分钟的变形观测量，选取的历元为：

      （X×4-45）~（X×4+45）

                    其中：（X×4-45）>0

B、计算，从开机时刻起算，第Y历元时的变形观测量，选取的历元为：

      （Y-45）~（Y+45）    其中（Y-45）>0

    （2）观测数据成果的分析

    我们可以推导出如下的变形量计算公式：

   其中：j变形点号，i为某时刻（或历元），o为基准点，

        A为某一认定的值或固定历元时的值。

       X、j、i为i时刻j点的向量值（变形点相对于基准点）。

   这样，我们针对算例，2号点为基准点，我们可把1号点及3号点所计算的相对2号点的瞬时（或然）值，绘成如图一、图二所示的变形点相对于参考点的位置变化图。当然，为了能更为准确地显示它的变化规律，我们可以把数据分成若干个历元进行处理，则它的运动轨迹将更能较好地反映出来，在此，为了明了起见，只选择了如附表所列数据进行了展点（计算展点时，须对同步闭合差分配）。

   由图一和图二，我们可以清楚地看出该变形点是基本上呈什么样的运动形式位移的。

5 结论

   综合以上所述，可得出GPS技术在工程变形监测中的应用，有以下的步骤：（1）根据监测的目的，在图上选点，然后到野外踏勘，以保证所选点位满足布网的要求和野外观测所具备的条件，最后得到要施测的概略点位；（2）按照每个点发出三条独立基线且边长分布较为均匀的原则并根据接收机台数的多少和布网原则，设计网的观测图形，并选定可能要追加施测的路线；（3）给定所需的可区分度指标（或精度指标）进行计算，直到达到给定的要求为止，最后得到增加独立观测基线后的最终施测方案；（4）观测时段和周期的设计；（5）采用连续观测分历元数据处理，或间断观测，整段数据处理方法进行数据计算；（6）利用观测值的统计，成因分析，一元或多元线性回归，逐步回归，图表等形式进行数据分析。

   由以上方法可知，利用这一监测方式，可以解决同一变形体，由于受阳光、风力、外界因素干扰（如大桥上行驶的车辆对大桥的影响）等，在某一时刻时，它们之间最或然的相对关系（相对位移及扰度等）及它们相对于参考基准点的绝对位移。分析它们各项因子对变形体影响的显著性，为工程及设计，运营管理单位提供最可靠的最直观数据及分析资料。

本文在撰写过程中，得到了上海市测绘院朱妙珍高工等的悉心指导，在此，深表感谢。

     表1：                 对 应 图 中 位 置1                  1~90历元

     表2：                  对 应 图 中 位 置2               90~180历元

     表3：                 对 应图 中 位 置3              180~270历元

     表4：                 对 应图 中 位 置4               270~360历元

     表5：                                                      1~360历元

      表6：               对 应图 中 A、C 位 置             45~135历元

      表7：                 对 应图 中 B、D 位 置          135~225历元

                                            Y(mm)

   注：此图中所采用数据均以WGS—84坐标系中的x、y分量所计算的，按二维方式的示意图，为了更确切地反映实际，应按地方所采用的投影平面进行坐标转化和高斯投影，故为能说明问题明了起见，在此未做此项，特此说明，图二同理。