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

一、实习目的 《测绘程序设计》是一门理论与实践并重的课程，课程设计是测量数据处理理论学习的一个重要实践环节，可以看做是在学习了专业基础理论课《误差理论与测量平差基础》课程后进行的一门实践课程，其目的是增强学生对测量平差基础理论的理解，牢固掌握测量平差的基本原理和公式，熟悉测量数据处理的基本原理和方法，灵活准确地应用于解决各类数据处理的实际问题，回顾所学习过的数据结构、C++高级语言程序设计、误差理论与平差基础、概率统计等课程知识，并能用所学的计算机知识，编制程序来计算实现矩阵（数组）、条件平差/间接平差以及水准网东平差的运算。 实习内容 1、 了解角度值与弧度制的关系，编写简单的程序，实现角度制与弧度制的相互转化； 2、结合C++高级语言程序设计和数据结构，回顾矩阵的运算（矩阵的加减法、矩阵的乘法、矩阵的转置、矩阵的求逆），矩阵的运算是进行平差的主要程序，无论是哪一种平差模型，都不可或缺。因此，编写程序实现矩阵的运算十分重要。 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．和普通矩阵的乘法不同，乘法并不要求两个被乘矩