第三章 水准网平差程序设计.pptVIP

第三章 水准网平差程序设计 本章难点： 1、近似高程计算 2、最短路线的计算 3、误差方程及法方程的构建 4、直接计算出法方程的系数矩阵BTPB和常数矩阵BTPL 该课程中所采用的平差模型为间接平差，即所选的独立参数的个数等于必要观测数，这样可以将每个观测值表示成这t个参数的函数，组成观测方程。 间接平差的函数模型： 采用的平差准则为最小二乘法: 从而平差值结果为： 4.解算法方程，求取未知参数； 5.精度评定。 平差程序设计不仅仅是“写程序”，测量程序设计还包括程序功能设计、平差模型选择、算法选择、数据接口设计等内容。 一般应按数据处理和计算功能的划分，将网平差程序分为若干独立函数，每个函数完成特定的计算或操作，当需进行某种平差时，再调用这些函数即可。 要整洁、简练、实用、方便有关数据的输入或输出，必要的话要考虑和用户的交互，尽量不要画蛇添足。 水准网平差输入界面一 水准网平差输入界面二 水准网平差输入界面三 四、高差闭合差的计算与检查 为了检查水准高差的质量，水准网平差前一般进行附合水准路线闭合差计算与多边形闭合差计算，并进行闭合差的检核。 1.附合水准路线高差闭合差的计算 a）最短路线搜索算法 “邻接点”：若一个点借助 另一个点与目标点发生了 联系，称另一个点是这个点 的“邻接点”。 如果一个点与目标点有直接联系，我们也把目标点称为这个点的邻接点。 显然邻接点不是唯一的，但是我们规定每个点只有一个邻接点；当网中每个点的邻点都唯一确定时，顺着邻接点又可找到邻点的邻接点，这样一层一层地找下去，最终找到目标点，每个点到目标点的路线是唯一确定的。 根据邻接点的概念可知，任何一点沿邻点到目标点的路线长度是其邻接点到目标点的路线长度再加上该点与其邻接点之间的路线长度。当一个点有多个点可以作为邻接点时，只选择路线最短的一点作为邻接点，这样就可保证沿邻接点到目标点的路线是最短路线。 设p1选择p2作为邻接点，P2到目标点的路线长度为S2, p1到p2的路线长度为S12，则p1到目标点的路线长度S1为 S1=S2+S12 若还有一点P3也是P1的邻接点，则经由P3到目标点的路线长度为 S’1=S3+S13 为了寻求最短路线，显然当S1S’1时，应选p2作为p1的邻接点，反之应选择p3作为p1的邻接点。因此寻求最短路线的过程转化为给每个点寻找邻接点的过程，当每个点都有了邻接点，且这个点连接的路线又是最短路线时，搜素工作即告结束。 综上所述，设目标点点号为k1，搜索各点到k1的最短路线的计算过程如下： （1）定义两个数组。一个数组Neighbor用来存储每个点的邻接点点号；另一数组S用来存储每个点到目标点的路线长度。 （2）将S(k1)赋值为0，数组的其它元素赋值为1.0e30（每点到目标点的初始路线长度为无穷大），Neighbor(k1)赋值k1，Neighbor数组的其它元素赋值为-1（表示没有邻接点）。 （3）循环访问每一个观测值，取出每个观测值的起点点号（假定为p1）和终点点号（假定为p2 ） 根据p1、p2得到p1到目标点的路线长度S1和p2到目标点的长度S2，再结合p1、p2的路线长度S12，判断p1是否可以作为p2的邻接点，以及p2 是否可以作为p1的邻接点。 那么p1可作为p2的邻接点的条件为： S2S1+S12 …(条件1) p2可作为p1的邻接点的条件为： S1S2+S12 …(条件2) 若条件1成立时,Neighbor(p2)=p1,S(p2)=S1+S12； 若条件2成立时,Neighbor(p1)=p2,S(p2)=S1+S12。 （4）当水准网的规模较大时，在循环访问每一个观测值后某些点到目标点的路线还不是最短路线，甚至有些点还没有找到邻接点，需要转（3）继续为寻找路线最短的邻接点，知道所有的观测边不再满足条件1和条件2，这时表明每个点都已经找到了到目标点最短的邻接点，搜素工作结束。 最短路线搜索函数FindShortPath（p as integer, exclude as integer, int neighbor()as integer,diff() as double, S() as double） 从一个已知点出发用观测高差依次推算其它各点的高程，最后闭合到另外一个已知点上，闭合点上推算高程和已知高程值之差称之为附合水准路线闭合差，简称路线闭合差。设k1、k2为两个已知点，已知高程值分别为H1、H2，两点之间的推算路线由高差h1,h2……hn，各测段路线长度分别为S1,S2 ……Sn，路线闭合差为 式中，正负号取决于高差起始点到终点的方向与推算路线是否一