测绘程序设计来自测绘魏志刚.pptxVIP

§4-3 平面控制网平差程序设计 控制网的平差方法基本上分为两类：条件平差法和间接平差法。 在手工计算时，主要考虑的是选择哪种法方程使得阶数较少，又便于解算的方法，因此平差方法的选择依控制网的具体情况而定。 然而对于电算来说，主要考虑的是程序编制的难易度和通用性。于是用条件平差方法编制通用程序就必须顾及所有可能的条件式，而要做到这一点，即使只要求程序适合于测角网，程序编起来也是相当困难。 但是对于间接平差来说，情况就完全不同了。因为无论任何一个平面控制网，通常只包含两类观测值，即：方向观测值和边长观测值，;;ajk、bjk称j、k方向的方向系数;对于某i测站其误差方程式为; 3）方向系数aij、bij的计算;6)法方程式的组成; 2．边长观测值误差方程式 1）边长观测值误差方程式的一般形式; 二、方位角条件与边长条件的列立 1．方位角条件; 三、法方程式的组成与解算 1．法方程式的组成 2．法方程式的解算 高斯消去法; §4-4 高程控制网平差程序设计 这里主要介绍以测段的高差为平差元素、待定点的高程为未知数的间接平差法的解算过程。 一、高程控制网间接平差的步骤及总框图 1．高程控制网间接平差计算步骤 (1)计算待定点的近似高程X＝(x1，x2 ．．．)T (2)列出误差方程 (3)组成法方程 (4)解算法方程，求得dX＝(dx1，dx2．．．)T (5)求得平差后的高程X＝X十dX (6)如果dX中的绝对值最大的一个小于给定的值，则平差结束，否则由新的近似值再一次平差， 即重复第(2)、(3)、(4)、(5)步。 (7)精度评定。; 2．高程控制网间接平差流程图 高程控制网间接平差是以测段高差为平差元素。一般为三角高程网和水准高程网的平差计算。要能够评定网中任意点高程和任意两点间高差的精度。主要包括以下三个主模块。 (1)数据输入模块 (2)平差计算模块 (3)精度评定及成果输出模块; (2)简变约定(WORD为无符号整型) WORD IKP：已知点点数 WORD IUP：未知点点数 WORD IQ：总点数。IQ??IKP+IUP WORD IGL：已知概略高程点数 WORD IOH：高差观测个数 WORD IZ：未知数总个数。IZ＝2*IUP WORD UD：评定精度的个数 Double EM：验前单位权中误差(mm); (3)数组(为了叙述方便，忽略掉各一维数组的零单元，因此说明数组时下标加1) Double DH[10H+1]：高差观测值数组 Double PDH[10H+1)：定权元素数组 double H[IQ]：高程数组 WORD ZDHA[I0H+1]，ZDHB[I0H+1]：高差的两个端点点号数组。高差观测值是有方向的，ZDHA[I0H+1]为起点，ZDHB [IOH+1]为终点。 DWORDZJD[2*IQ+1]评定精度元素数组，以XXXXXXXX方式存放，前四位为一个点号，后四位为另一个点号。评定高程精度时只有一个点号。 Char Name[IQ+1][18]：以点号为序存放各点的点名。不得大于18个字符。 Double RM[IZ*(IZ+1)／2+1]，RW[IQ+1]：法方程系数和法方程常数项数组。;试按间接平差法求B、C、D点高程的平差值。 ; 例4—2 如图所示的水准高程控制网，已知点2个， 未知点3个，总点数5，已知概略高程点数0，高差观测值个数7，评定精度的个数3。网形编号如图4—2所示。;各测段高差及距离; 二、待定点高程计算 1．原理 待定点的高程计算是为了得到未知数的高程近似值。在端点号数组ZDHA[IOH+1]和ZDHB[I0H+1]中，如果该两点有一点是已知或已算出的，而另一点未知，则可由其中已知高程的点与未知点间的观测高差运算得到未知点的高程。为了区别哪些点已知或已算出，哪些点待计算，需要设立标志数组NameBz[IQ]，当NameBz[I]＝1时，说明第1点高程已知或已算出，NameBz[I]=0时，说明第1点待计算。为了防止由于编号不当或观测值排列顺序的限制而一次不能完全推算出所有的待定点，采用多次推算的方法来解决。 2．程序;;