文档介绍：《测绘程序设计》实****报告院系:姓名:班级:学号:2012年4月《测绘程序设计》实****报告一、总述:本学期我们《测绘程序设计》课程共有6周课程的实****程序设计的题目是基于上学期《误差理论与测量平差》课程的设计,主要有水准网平差、边角网平差、GPS网平差。由于时间的问题,本学期程序设计主要是水准网平差,有余力的同学继续向下做。二、实****时间地点及小组安排:1、实****时间:每周周二、周五下午四节课;2、实****地点;信息楼208机房;3、小组安排:以小组为单位,每组5-6个人,每个小组指定一位编程能力较好的同学作为组长,负责小组人员到课情况、作业完成程度、指导成员上机以及负责日常作业的检查。三、课程上机实****1、实****内容概述:第一次实****由老师系统的介绍了本学期上机实****的大体概况,讲解了测量平差的具体具体类型,给大家重新熟悉了间接平差的具体步骤。之后由我们编写了一个小程序:角度与弧度的相互转化;第二次实****大家已经完成上次的实****任务,之后我们开始间接平差的实****主要是矩阵间的操作;第三次课开始,每次上机操作之前先由老师简单介绍本次课程需要完成的任务,之后给大家看PPT讲解程序设计的具体操作。大约花费半个小时,剩下的时间我们自己编写代码。期间有不懂的问题可以随时问老师、学长或者与其他同学交流讨论;几次课之后,我们都完成了间接平差的实****由老师负责检查。检查完成之后,大家进行水准网平差的实****2、具体实****程序::程序设计基本思想:应用角度与弧度的比例关系,进行相互转化。JD(角度)=HD(弧度)*180/π;HD(弧度)=JD(角度)*π/180;::程序设计基本思想:采用指针数组,为矩阵的行列动态分配空间,然后根据间接平差的步骤,一步步求得最后的结果。:间接平差的步骤:①:根据题目要求,写出误差方程式v=Aδx+l;②:根据误差方程式系数阵A求N,N=AT*P*A(AT表示A的转置,P为权阵);③:根据误差方程式常数项l求C,C=AT*P*A;④:由高斯约化法或行列变化求N-1;⑤:求δx,δx=-N-1*C;⑥:将δx带回原误差方程,求误差改正数v,v=Aδx+l;:主要程序分析:①:矩阵相加:应用for循环控制条件,C[i][j]=A[i][j]+B[i][j];②:矩阵相乘:应用for循环控制条件,sum=sum+A[i][k]*B[k][j];③:矩阵求逆:采用行列变换法求出矩阵的逆矩阵。先初始化一个单位阵a[i][j],然后将N化为单位阵,相应的a[i][j]才用相同的变换,化为N-1;:程序注意事项:①:为指针数组动态分配空间,注意行列的大小。如A:m行n列,P:m行m列,l:m行1列,N:m行m列等等;②:一般的数据都采用double型,力求精确;③:应用srtprecision()函数设置精度,同时注意包含头文件#include<iomanip>。:程序设计收获:①:重新温****了间接平差的解题思路和方法,加深印象;②:掌握了一门技能,用程序实现间接平差的解题,简单、高级;③:掌握了C++函数的运用,熟悉了相应的C++语言技巧,如setpricision()函数的运用;④:与同学交流了思想,加深情谊,锻炼了团队合作意识与能力。::程序设计基本思想:已经编写了间接平差的程序,其中的矩阵操作可以直接运用。水准网平差主要是实现对误差方程式系数矩阵、常数项、权阵的输出,要采用相应的程序技巧将水准网的路径转化为矩阵中的数字,以此实现矩阵的输出。:水准网平差的步骤:①:采用程序设计技巧,控制相应的循环条件,将误差方程式的系数阵A、常数项l以及权阵P求出,组成误差方程式v=Aδx+l;②③④⑤⑥步同间接平差;⑦:求观测值改正之后的值L_L=L+v;⑧:求各未知点的高程值H_H=H+δx。:主要程序分析:①:系数矩阵的求法:◆先初始化系数矩阵A[m][n]为0矩阵;◆然后控制循环条件,依次输入各段水准路线的起点、终点;◆依次判断各水准路线的起点、终点在矩阵A中的位置,若是起点,则把原先的值0赋值为-1;若是终点,则把原先的值0赋值为1,其余的不变,求得A。②:常数项的求法:◆对应水准网点的编号,初始化标记数组a[P](P为总点数),已知点的赋值为0,未知点的赋值为1;◆在输入水准路线起始点的同时,判断,如果起点标记为0,终点标记为1或起点标记为1,终点标记为0,则求出未知点的近似高程H[],并求得常数项l,如果起始点分别标记为0和2或起始点分别标记为2和0,并把未知点点号标记为2;这样就把水准路线中带有已知点路线的未知高程点求得了近似高程H和常数项l;◆控制循环条件,判断标记数组,如果有1值,判断,如果起始点分别标记为