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

一、实习目的 《测绘程序设计》是一门理论与实践并重的课程，课程设计是测量数据处理理论学习的一个重要实践环节，可以看做是在学习了专业基础理论课《误差理论与测量平差基础》课程后进行的一门实践课程， 3、在实现矩阵的运算之后，学会文件的读写操作，再加上一些程序实现间接平差/条件平差的运算。 4、在已实现的间接平差程序上实现水准网平差，将观测数据输入到文本，学会从观测数中读取相关矩阵。 三、实习过程 （一）角度与弧度的相互转化 C++程序设计中，关于角度的计算以弧度制为单位，而在测量以及具体工作中我们通常习惯以角度制为单位。这样，在数据处理中，经常需要在角度制与弧度制之间进行相互转化。这里，我们利用C++数学函数库math.h中的相关函数完成这两种功能。 这里，我们使用double类型数据表示角度制数和弧度制数。例如：123度44分58.445秒，用double类型表示为123.4458445，其中分、秒根据小数位确定。 在角度制与弧度制的转化中，涉及如下图所示的两个环节。 程序设计基本思想： 应用角度与弧度的比例关系，进行相互转化。 JD（角度）=HD（弧度）*180/π； HD（弧度）=JD（角度）*π/180 （二）文件读写操作 1.思路及方法: 头文件#includefstream，以test.text文件为例 (1)读取文件: ifstream fin；fin.open(“test.txt”)检查是否正常打开.is_open())；if(fin.is_open())；if(fin.fail())if(!fin.good())；if(!fin)；读取文件的下一行；fin.getline(buf,90);；getline(fin,line); （2）关闭文件: fin.close();将数据写入文件ofstream fout;foutdata;fout.close(); 使用之前要先检查文件是否正常打开，结束以后要记得关闭文件。 （三）动态创建数组和矩阵的运算 1.动态创建数组 ： （1）动态数组创建有两种方法： 方法一：用一位数组代替二维数组 方法二：二维数组 （2）主要思想是：使用new和delete两个运算符。new的功能是动态分配内存，语法格式为：new 类型名 T（初值列表）；该语句在程序运行过程中申请分配用于存放T类型数据的内存空间，并使用初值列表中的值进行初始化。如果内存申请成功，new运算返回一个指向新分配内存首地址的T类型的指针，如果失败，就返回空指针。delete的功能用来删除由new建立的对象，释放指针所指向的内存空间。 2.矩阵的运算（加、减、乘、转置及求逆） 在数据处理程序中，经常需要用到矩阵的基本计算，包括矩阵的加减、乘积、转置、求逆等。在测量数据处理中，还需要用到一些具体的、有规律和统一的矩阵计算。在此，为程序简单起见，我们用C++语言的数组表示矩阵，定义矩阵的最大维数为MAX（定义矩阵存储空间大小），函数中对同类计算使用C++重载的概念进行处理，从而使程序清晰、应用方便。 3.思路及方法： 平差计算过程中涉及到的都是矩阵的计算，因此我们要编程实现矩阵的各种计算。具体来说： 矩阵的加（减）法：首先要判断两个矩阵是否具有相同的行数和 数，若行数或列数不同，则输出“无法相加（减）！”。 （2）矩阵的乘法：对于乘法首先仍要判断两个矩阵的行数和列数，判断的准则是第一个矩阵的列数要等于第二个矩阵的行数，只有这样两个矩阵才能够相乘。 对n×m阶矩阵A和p×q阶矩阵B，A和B的乘法运算可定义为： 由上面的式子可以看出，乘积AB表示矩阵A的所有元素与B之间的乘积组合而成的较大的矩阵，BA则完全类似．AB和BA均为np×mq矩阵，但一般情况下ABBA．和普通矩阵的乘法不同，乘法并不要求两个被乘矩阵满足任何维数匹配方面的要求A^(-1)=(1/|A|)×A* ，其中A^(-1)表示A的逆矩阵，其中|A|为矩阵A的，A*为矩阵A的 逆矩阵的另外一种常用的求法： (A|E)经过得到(E|A^(-1)) 只有满足一下条件才能进行矩阵的求逆运算： a秩等于行数 行列式不为(或列向量)是组存在一个矩阵,与它的乘积是作为线性方程组的系数有唯一解满秩可以经过初等行变换化为单位矩阵伴随可以表示成初等矩阵的乘积它的转置可逆它去左(右)乘另一个矩阵,秩不变 间接平差法（参数平差法）是通过选定t个与观测值有一定关系的独立未知量作为参数，将每个观测值都分别表达成这t个参数的函数，建立函数模型，按最小二乘原理，用求自由极值的方法解出参数的最或然值，从而求得各观测