测绘程序设计课程实习报告2013.docVIP

测绘程序设计课程实习报告 一、实习目的 巩固之前所学的c++类定义、文件操作等各方面知识，掌握矩阵计算的基本算法，熟悉并熟练运用平差模型，掌握测量平差的基本原理和公式，并处理测量观测数据。 二、实习内容 该课程的实习主要为c语言及c++编程基础，矩阵操作算法，条件平差，间接平差，以及在水准网平差中的应用。 三、实习过程 实习1（文件的基本操作以及角度、弧度的转化） 准备 C++文件操作函数包含于头文件fstreamofstreamifstream，一般情况下使用fstream即可。在计算过程中需要文件输出格式控制，函数包含于iomainip。c++的流操作函数包含于iostream，相应的流函数在命名空间“std”。 角度转弧度的过程中，首先判断正负，对符号进行不同情况的处理。将度数的分，秒位分别除以60、3600，然后把得到的结果与度位绝对值相加，判断符号。 弧度转角度的过程中，同样判断正负，取小数位乘以60取整数得到相应角度的分，再取剩余的小数乘以60即得到相应角度的秒。 代码 #include iostream #include math.h #include fstream #include string using namespace std; const double pi = 3.14159; int main() { fstream file(a.txt); //定义文件操作对象 string aa; fileaa; coutaaendl; //输出文件中的内容 int a = 1; //定义重复输入使用的变量 while(a) //限制循环，重复输入 { float a1,a2,a3; coutendl; double deg; cout输入度endl; cina1; cout输入分endl; cina2; cout输入秒endl; cina3; cout角度为 a1° a2′ a3″ endl; 　　　　if(a1 = 0) //判断输入度数的正负 deg = a1 + a2/60 + a3/3600; else deg = a1 - a2/60 - a3/3600; double radian; 　　　　radian = deg/180*pi; cout转换得弧度为 radianendl; coutendl; cout输入弧度： endl; 　　　　cinradian; 　　　　bool tm = 0; deg = radian*180/pi; if(deg0) 　　　　 { tm = 1; deg = -deg; 　　　　 } a1 = int(deg); a2 = int( (deg-a1)*60 ); a3 = ( (deg-a1-a2/60)*3600 ); if(tm == 1) { a1 = -a1; tm = 0; } cout转换后角度为 a1° a2′ a3″endl; 　　　　 } return 0; } 运行结果 输入度 20 输入分 55 输入秒 32 角度为 20° 55′ 32″ 转换得弧度为 0.36522 输入弧度： 3.14 转化后角度为 179° 54′ 32.0387″ Press any key to continue 出现的问题 角度转换过程中需要考虑角度的正负号问题，角度转化可以添加是否超过59的判断。 实习2（概率计算） 准备 在误差检验、粗差探测等平差计算中经常需要计算某些特定分布函数的反函数值，为满足计算需要，使用CProbability类对相关计算进行封装。 代码及运行结果与书上相同。 实习3（矩阵基本运算） 准备 矩阵的运算有加法、减法、乘法、转置、求逆，封装在一个类中。为了方便之后的计算与使用，类中还添加了从文件中读入矩阵、输出矩阵，向控制台输入输出矩阵的函数。 代入运算的矩阵均用动态的一维数组表示。先定义double型指针，在确定矩阵大小后动态申请空间，row行col列的矩阵定义为：double *mat = new double[row*col]; 第i行j列的元素为mat[i*col+j]。 矩阵的加法、减法只需遍历两个矩阵的每个元素，相加之后放到结果矩阵中（如有必要先判断相加的两矩阵行数与列数是否相等）。 乘法计算中设相成的两个矩阵mat1、mat2的行列分别为row1、col,col、col2，则所求结果矩阵第i行j列的元素 用三个循环即可计