数据结构课程设计——导线网平差程序的设计与实现.docxVIP

设计二：导线网平差程序的设计与实现 一、设计目的 立足于《数据结构与测绘软件开发》这一课程的课堂教学及其实验课程设计，为着实提高学生基于计算机辅助的方式切实解决工程实际问题的动手能力，通过本实习，一方面，使学生深入了解课堂所学知识，另一方面，通过实践掌握测绘行业软件设计与开发的基本方法，深刻掌握矩阵运算、曲线/曲面拟合的数值解法，掌握不同类型的典型测绘软件设计方法，使得学生初步具备编写测绘软件常用算法的能力以及开发中小规模测绘专业软件的能力。 二、设计内容 有导线网如图，观测了14条边长和16个转折角，已知测角精度，测边精度为。已知A、B、C、D、E、F点的坐标（无误差），如下表： 表1 已知点数据 点号X(m)Y(m)A5256.9534520.068B5163.7524281.277C3659.3713621.210D4119.8793891.607E4581.1505345.292F4851.5545316.953表2 角度观测值 编号角度观测值（° ′ ″）编号角度观测值（° ′ ″）1163 45 049169 10 30264 58 371098 22 043250 18 111194 53 504103 57 3412111 14 23583 08 051379 20 186258 54 1814268 06 047249 13 1715180 41 188207 32 3416103 23 08表3 边长观测值 编号边长观测值（m）编号边长观测值（m）1238.6198241.5602170.7599224.9963217.86910261.8264318.17311279.8405245.63512346.4436215.51413312.1097273.82914197.637试编程实现该导线网的平差，给出计算结果，并对其精度进行评价。 三、关键问题描述 3.1 未知点近似坐标计算 平面控制网进行平差计算时需要计算未知点的近似坐标 1.坐标计算公式 1、2点的坐标已知，并观测了1-2、1-3的夹角，根据这些数据可以求出3号点坐标 根据1、2两点的坐标，可以反算出1、2方向的方位角T12,3号点的坐标为 式子中S13为观测边长，为观测角度 2.计算流程 从读入的数据循环计算未知点的坐标，已计算出的坐标当做已知坐标的点处理参加下次计算，以此类推，逐步计算出未知点的坐标 3.实现算法 CMatrix CPlaneNetAdjust::XYJS() { CMatrix _XYJS(Pnumber,2); double T12; for(int i=0;iTnumber;i++) { int k1=cezhan[i]; double x1=xy[k1].X; double y1=xy[k1].Y; if(x1=0) { continue; } int k2=huoshi[i]; int k3=qianshi[i]; if(xy[k2].X0xy[k2].Y0) { T12=GetT12(k1,k2); } double s12=Gets12(k1,k2); double s13=Gets12(k1,k3); double T13=T12+guancejiao[i].Guancezhi; double dx=s13*cos(T13); double dy=s13*sin(T13); xy[k3].X=xy[k1].X+dx; xy[k3].Y=xy[k1].Y+dy; } for(int i=0;iPnumber;i++) { double temp1=xy[i].X; double temp2=xy[i].Y; _XYJS.setValue(i,0,temp1); _XYJS.setValue(i,1,temp2); } return _XYJS; } 3.2 误差方程列立 1.理论分析 平面控制网的误差方程都是非线性方程，必须引入参数近似值将误差方程线性化， 取的充分近似值 ，是微小量，在按台劳公式展开时可以略去二次和二次以上的项，而只取至一次项，于是可对非线性平差值观测方程式线性化，于是有如下的式子 对于观测角的改正数有 对于边长观测值的改正数有 2.实现算法如下： CMatrix C