总的设计流程图:
输入已知点信息和观测信息 数据编辑 NO 闭合差检验合格 YES 计算近似坐标 组成边长观测误差方程组成方向观测方程 组成法方程并求逆 计算坐标改正值及验后单位权方差 判断参数改正数 精度评定 输出计算结果
代码书写与调试:
由于上面已经展示过几乎所有的代码(根据每一个模块),真正运行整一个程序的时候却是按照步骤一步步地一个模块地运行,但是问题一定会出现的,而且还不是很少,关键就是怎么样通过一步步地发现问题所在,并根据实际情况去试调程序,合理地按照平差理论去运行程序。一下列举了一些出现的主要问题。 方位角
由于所提供的数据时边角三角网的数据,这样就会涉及到方位角问题,而且还不止一次的出现方位角的运算(如待定点的近似坐标的求解过程中出现过),这就应该惊醒一个重要的细节:就是所求的方位角是否合理?怎么样去判断呢?就是通过if(elseif)语句去判断其与0和360度或者2*pi 的关系,在去改正方位角。原则上一旦出现方位角的求解出都是要用这个语句去判断,这样才能保证精密。和提高正确性。 单位问题
单位问题出现在两处,就是观测值的单位权的求解处,二是误差方程的系数求解处。
由于这两者是在后面的法方程求解时都要用上所以,在此之前要保证两者的单位相对统一,这个统一要分两处,边长观测量方程系数和方向观测方程系数,特别注意的是原始数据提供的观测精度(两种)的单位,和误差方程的系数的单位,而在实际中,我们组应用的是方向观测用秒为单位,而变成观测以米为单位,这样出来的结果最后的平差数据处理,如精度评定,误差椭圆计算等,也要注意单位问题。 比例参数
根据实际求出的误差椭圆参数(E F fi)的大小和所绘制的三角网的大小,这两者之间通过比例参数的调整才能更好地在显示误差椭圆。 弧度转角度
对于弧度转角度这个小问题,我不得不提,这是一个关键问题。为什么?因为我们组把其他问题都解决了后, 组的坐标改正数,误差椭圆参数等都很不正常的显示在电脑屏幕上,一次次地打击我们的信心。最后,我在仔细地对照书本,细心查询,观察异同发现,一个最不关紧要的差异,就是弧度转角度的小问题,我们原来采用的是MATLAB系统自带的弧度转角度的函数,结果我们发现这个问题直接就是导致待定点的近似坐标差异特别大,尤其在待定点3的x,y坐标上面,差了好几米(跟提供的平差结果数据相比),这样后面的就影响到误差方程系数的不正常,和平差计算得出的改正数,单位全中误差,点位误差和误差椭圆参数的异常。而这样我不得不自己再编两个弧度和角度相互转换的函数文件。于是重新采用转换函数程序之后,我们的数据就非常正常的出来的。
平差程序设计成果展示
进入界面
导入数据文件
数据文件的查询 查询结果:
近似坐标的计算 计算结果:
控制网点为分布
显示平差计算成果
百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读,免费范文网,提供经典小说综合文库MATLAB测量平差程序实习报告..(5)在线全文阅读。
相关推荐: