水准网平差程序的创新设计.pptVIP

系数矩阵B构建算法描述： 1、定义一个存储系数B的数组b()，其行数为观测测段数，列数为所选必要参数个数；定义一个存储常数阵L的数组L()，其行数为测段个数，列数为1； 2、取出第j个测段（测段数据存储在测段信息数组中的），首先判断该测段起点是否为已知点（也即其在已知点数组a()中的下标位于起算高程点个数内，如：起算高程点个数为5，起点在已知数组中的下标为i，当i小于等于5时，说明该起点为已知点，若i大于5，则说明该起点为未知点，同样终点是什么类型的点也通过该方法判断），记下其在a()中的下标i1（原始已知点数为i0）；判断该测段终点是什么类型的点记下其在a()中的下标i2； a)、若该测段的起点为未知点，且终点也为未知点，则矩阵B的第j行、第（i1-i0）列的元素为-1，第j行、第（i2 –i0）列的元素为1；那么第j行其它元素都为 0； b)、若该测段的起点为已知点，且终点为未知点，则矩阵B的第j行、第（i2-i0）列的元素为1，该行其它元素都为0； c)、若该测段的起点为未知点，且终点为已知点，则矩阵B的第j行、第（i1-i0）列的元素为-1，该行其它元素都为0； 3、根据每测段起点和终点点的类型，就可以计算出这个误差方程的常数项，把计算出的数值放到L()中去就可以了； 4、遍历完所有的测段，即可构建出矩阵B，即其所有的元素都存储在b()中了。 七、法方程的组建 有了B和L，基于最小二乘原理，对所选参数求偏导后的式子，推出法方程，即： 直接计算法方程系数阵和常数阵（累加法） 组成法方程是间接平差的关键步骤，由间接平差模型可知 ，在误差方程系数矩阵B，误差方程常数项以及观测值的权P确定时，法方程系数矩阵和自由项向量分别为BTPB和BTPL，组成法方程就是矩阵的乘积。但在水准网平差程序设计时，为了节省存储单元，可不必为误差方程矩阵B开辟存储数组，权阵也只保存对角线上的元素，组成法方程的计算就不是矩阵乘积运算了，而是根据观测值的起始点号、终点号、高差值和观测值的权直接计算出法方程的系数阵BTPB和自由项向量BTPL。 应将pk、-pk累加到法方程系数矩阵的什么位置？若第K个观测高差的起点点号和终点点号分别为i、j，则需将pk累加到法方程系数矩阵的第i行第i列和第j行第j列，将-pk累加到法方程系数矩阵的第j行和第i列。同理，应将pklk和-pklk分别累加到法方程常数项的第j行和第i行中。 组成法方程的过程就是以观测值总数为循环界，逐个观测值计算累加项的过程，在第k 个循环中，所做的工作如下： （1）获得高差的起点点号i和终点号j； （2）获得起点和终点的高程值Hi和Hj； （3）计算误差方程的自由项lk=Hj-Hi-hk； （4）将pk、pklk累加到法方程系数阵和法方程自由项的相应单元中。 可以将法方程系数矩阵和法方程常数项矩阵的计算定义成一个函数Ca_ATPA()。 void CLevelingAdjust::ca_ATPA() { int t=m_Pnumber; for(int i=0; it*(t+1)/2; i++) ATPA[i]=0.0; for(i=0; it; i++) ATPL[i]=0.0; for(int k=0; km_Lnumber; k++) { int i=StartP[k]; int j=EndP[k]; double Pk=P[k]; double Lk=L[k]-(Height[j]-Height[i]); ATPL[i]-=Pk*Lk; ATPL[j]+=Pk*Lk; ATPA[ij(i,i)]+=Pk; ATPA[ij(j,j)]+=Pk; ATPA[ij(i,j)]-=Pk; }} 八、水准网平差定权 水准网平差中还有涉及到权阵P的确定，有了权阵P ，和我们上面所构建出来的B和l ，就可以很容易的进行参数改正值和观测值改正值的解算了。 我们一般怎样来定权? 设有观测值Li(i=1,2,…n)它们的方差 为 若选定一常数 , 定义 ，并称P为观测值Li的权。对于一组观测值，其权之比等于相应方差的倒数之比。这表明，方差越小，其权越大；或者说精度越高，其权越大。因此权可以作为比较观测值之间的精度高低的一种指标。 结论: （1） 选择了一个 ，即对应一组权； （2） 一组观测值的权，其大小是随 的不同而异，但不论选何值，权之间的比例关系始终不变； （3）为了使全能够起到比较精度的作用，在同一问题中只能选取一个 值，不能同时选取几个 ，否则