2022年测绘程序设计实验指导书 2.pdf

测绘程序设计实 验 指 导 书华北科技学院土木工程系测绘工程教研室2010 年 8 月名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 29 页 - - - - - - - - - 实验一：线性方程组解算一、实验目的与要求初步掌握线性方程组的算法；初步掌握线线性方程组在计算机上实现方法。二、实验安排实验共需用 2 学时。实验每个小组1 人，用 C/VB等语言在计算机上完成程序代码的编写，并调试通过，然后用算例印证程序的正确性。三、实验步骤及要点在课余时间，提前准备线性方程组解算程序代码；上机时写入代码，并调试运行；用算例进行印证。四、实验指导（一）Turbo C 工作环境一个 C语言程序的实施是从进入Turbo C 的集成环境开始的，而进入C 语言的环境，一般有两种途径：从DOS环境进入和从Windows环境进入。1、从 DOS 环境进入：在 DOS 命令行上键入：CCD TC （指定当前目录为TC子目录）CTC (进入 Turbo C 环境) 这时进入 Turbo C集成环境的主菜单窗口，屏幕显示如下图所示。2、从 Windows环境进入：在 Windows 95/98 环境中，如果本机中已安装了Turbo C，可以在桌面上建立一个快捷方式， 双击该快捷图标即可进入C语言环境。 或 者从开 始菜 单 中 找到“ 运行 ” ， 在 运 行对 话框 中键 入“C:TCTC” ， “确定”即可。刚进入 TC环境时，光带覆盖在“ File ”上，整个屏幕由四部分组成，依次为：主菜单、编辑窗口、信息窗口和功能提示行（或称快速参考行）。（1）主菜单名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 29 页 - - - - - - - - - 显示屏的顶部是主菜单条，它提供了8 个选择项：File 处理文件（装入、存盘、选择、建立、换名存盘、写盘） ，目录操作（列表、改变工作目录） ，退出 Turbo C，返回 DOS状态。Edit 建立、编辑源文件。Run 自动编辑、连接并运行程序。Compile 编辑、生成目标文件组合成工作文件。Project 将多个源文件和目标文件组合成工作文件。Option 提供集成环境下的多种选择和设置（如设置存储模式、选择编参数、诊断及连接任选项）以及定义宏；也可记录Include 、Output 及 Library文件目录，保存编译任选项和从配置文件加载任选项。Debug 检查、改变变量的值、查找函数，程序运行时查看调用栈。选择程序编译时是否在执行代码中插入调试信息。Break/Watch 增加、删除、编辑监视表达式， 及设置、清除、执行至断点。在主菜单中，Edit 选项仅仅是一条进入编辑器的命令。其他选项均为下拉式菜单，包含许多命令选项，使用方向键移动光带来选择某个选项时，按回车键，表示执行该命令，若屏幕上弹出一个下拉菜单，以提供进一步选择。（二）Turbo C 环境中运行 C语言源程序的步骤(1) 编辑源文件在主菜单下，直接按 Alt+F 键，或按 F10后将光带移到“File ”选项上，按回车键，在“ File ”下面出现一个下拉菜单，菜单中有以下选项：Load F3 表示加载或装入一个文件Pick Alt+F3 从指定的文件列表中选择文件装入编辑器。New 表示新文件，缺省文件名为NONAME.C。Save F2 将正在编辑的文件存盘。Write to 将正在编辑的文件写入指定的文件中，若文件名已存在，则重写。Directory 表示文件目录。Change Dir 改变驱动器及目录名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 29 页 - - - - - - - - - OS shell 进入 Turbo C 命令行模式，命令EXIT 可返回集成环境Quit Alt+x 退出 Turbo C，返回 DOS 状态。建立一个新文件，可用光标移动键将“File ”菜单中的光带移到“New ”处，按回车键，即可打开编辑窗口。此时，编辑窗口是空白的，光标位于编辑窗口的左上角，屏幕自动处于插入模式，可以输入源程序。 屏幕右上角显示缺省文件名为NONAME.C， 编辑完成之后，可用F2 或选择“ Save”或“ Write to”进行存盘操作，此时系统将提示用户将文件名修改成为所需要的文件名。(2) 源程序的编译、连接直接按 F9键，或将菜单“ Compile”中的光带移到“ Make EXE file”项上，按回车键，就可实现对源程序的编译、连接。若有错误，则在信息窗口显示出相应的信息或警告，按任意键返回编辑窗口，光标停在出错位置上，可立即进行编辑修改。修改后，再按 F9键进行编辑、连接。如此反复，直到没有错误为止，即可生成可执行文件。注意： C程序的连接是在编译后自动完成的。(3) 执行程序直接按 Ctrl+F9键，即可执行 .EXE 文件；或在主菜单中（按F10进入主菜单）将光带移到“ Run”选项，按回车键，弹出一个菜单，选择“ Run ”选项，回车。这时并不能直接看到输出结果。输出结果是显示在用户屏幕上，在 TC屏幕上看不到，直接按复合键Alt+F5 ，或选择“Run”菜单中的“ User Screen”选项，即可出现用户屏幕，查看输出结果。按任意键返回 TC集成环境。另外，选择“Run”菜单下的“Run ”项，或直接按 Ctrl+F9 键，可将 C程序的编译、连接、运行一次性完成，即第3 步中包含有第2 步的工作。如果程序需要输入数据，则在运行程序后，光标停留在用户屏幕上， 要求在用户输入数据，数据输入完成后程序继续运行，直至输出结果。 如果运行结果不正确或其他原因需要重新修改源程序，则需重新进入编辑状态。修改源程序，重复以上步骤，直到结果正确为止。(4) 退出 Turbo C 集成环境名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 29 页 - - - - - - - - - 退出 Turbo C环境， 返回操作系统状态。可在主菜单选择“File ”菜单的“ Quit ”选项，或者直接按Alt+x 。在执行退出 Turbo C环境时，系统将检查一下当前编辑窗口的程序是否已经存盘，若未存盘，系统将弹出一个提示窗口，提示是否将文件存盘，若按“Y”则将当前窗口内的文件存盘后退出；若按“N ”则不存盘退出。（三）线性方程组解算具体算法参见讲义。利用高斯消去法，根据算法，写出程序代码。五、思考题如何利用其他解法完成计算机数据处理。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 29 页 - - - - - - - - - 实验二：矩阵解算与角度运算一、实验目的与要求初步掌握矩阵的加减乘除及矩阵的转置、求逆等基本运算在计算机上的实现方法；掌握度分秒转换在计算机上实现的方法。二、实验安排本实验需要 2 学时。实验每个小组1 人，用 C语言在计算机上完成程序代码的编写，并调试通过，然后用算例印证程序的正确性。三、实验步骤及要点1在课下写出程序代码，重点是矩阵求逆和度分秒转换；2上课时在计算机上写入程序代码，调试运行，通过后用算例进行测试。四、实验指导1. 矩阵解法参见教材，下面给出测试数据。122113312A5647854531A2. 角度互化问题在计算机数据处理中经常遇到，计算机要求三角函数中的角度必须是弧度制， 这样引用计算机中的三角函数进行计算时，结果才是正确的，反过来，求出的角度值也是一弧度为单位的，所以还需把弧度转换为度分秒的形式。五、思考题在什么情况下需要使用度分秒转化；在什么情况下利用矩阵解算程序。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 29 页 - - - - - - - - - 实验三：控制测量常用程序设计一、实验目的和要求初步掌握控制测量数据处理过程中的数据结构类型；初步学会控制测量常用软件设计方法、过程、程序代码编写方法及程序调试过程。二、实验安排本实验需要 4 学时。实验每个小组1 人，用 C语言在计算机上完成程序代码的编写，并调试通过，然后用实际算例印证程序的正确性。三、实验步骤及要点课余时间完成预先布置的控制测量数据处理问题，完成程序代码的纸上编写和设计工作； （上级题目可以根据自己的情况，任意选做一个即可）在计算机上写入程序代码，进行调试；用预先准备好的算例进行印证， 要求每人至少准备2 组典型算例，并打印测试结果，附在实验报告中。四、实验指导控制网数据结构： 在计算机中，带有结构性质的数据元素的集合为数据结构。 控制网数据结构主要研究控制网的数据存储格式和存储方法。 主要包括图的数据结构和控制网数据结构，常用的存储模式主要有混合结构、邻接表结构。在网图数据中，常用邻接表方式，平面控制网概算程序设计指导：1）流程图：名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 7 页，共 29 页 - - - - - - - - - 2）测算实例：如图所示平面控制网，已知点4 个，未知点3个，15 个方向观测值。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 8 页，共 29 页 - - - - - - - - - 程序代码自己独立完成。3高斯投影坐标换带程序设计指导1）坐标正算公式如下：2）坐标反算公式名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 9 页，共 29 页 - - - - - - - - - 五、思考题1思考如何编写高程控制网严密平差程序。2如何编写高程测量概算程序。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 10 页，共 29 页 - - - - - - - - - 实验四：工程测量程序设计一、实验目的和要求初步掌握工程测量中常用软件设计方法和技巧。二、实验安排本实验需要 4 学时。实验每个小组1 人，用 C语言在计算机上完成程序代码的编写，并调试通过，然后用实际算例印证程序的正确性。三、实验步骤及要点课余时间完成预先布置的控制测量数据处理问题，完成程序代码的纸上编写和设计工作；在计算机上写入程序代码，进行调试；用预先准备好的算例进行印证，要求每人至少准备12组典型算例，并打印测试结果，附在实验报告中。四、实验指导施工放样地面标定要素计算程序设计公式请参考工程测量学教材。 独立研究和设计其存储的数据格式，写出程序代码，上机调试通过运行。变形观测程序设计公式计算法请参考工程测量学教材。独立研究和设计其存储的数据格式，写出程序代码，上机调试通过运行。五、思考题编写变形观测程序时， 请思考数据处理与图形绘制为一体的程序设计实现方法。实验五： GPS程序设计一、实验目的和要求初步掌握 GPS 数据处理软件中的数据结构格式；初步掌握 GPS 平差软件的编写方法和思路；初步了解 GPS 数据通讯软件设计方法。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 11 页，共 29 页 - - - - - - - - - 二、实验安排本实验需要 4 学时。实验每个小组1 人，用 C语言在计算机上完成程序代码的编写，并调试通过，然后用实际算例印证程序的正确性。三、实验步骤及要点课余时间完成预先布置的控制测量数据处理问题，完成程序代码的纸上编写和设计工作；在计算机上写入程序代码，进行调试；用预先准备好的算例进行印证，要求每人至少准备12组典型算例，并打印测试结果，附在实验报告中。四、实验指导各种 GPS 数据格式参见教材；卫星位置计算程序设计指导：计算步骤和过程如教材所示，下面给出程序代码，同学们据此进行上机调试。卫星位置计算程序代码：名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 12 页，共 29 页 - - - - - - - - - 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 13 页，共 29 页 - - - - - - - - - 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 14 页，共 29 页 - - - - - - - - - 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 15 页，共 29 页 - - - - - - - - - 3数据格式转换程序设计（程序略）4GPS 接收机与计算机通讯程序设计（选作）五、思考题名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 16 页，共 29 页 - - - - - - - - - 如何编写实用的GPS 后处理软件。如何从 GPS 星历数据文件中提取数据。实验六： Excel 在测绘工程中的应用一、实验目的和要求初步了解电子表格在测绘中的应用方法。二、实验安排本实验需要 2 学时。实验每个小组1 人，用Excel在计算机上绘制电子表格，利用此表格实现测量外业记录和数据处理，然后用实际算例印证其可行性。三、实验步骤及要求1. 用 VBA完成前方交会自定义函数；用Excel 绘制电子表格，写入数据，用函数完成计算，输出计算结果，打印成果表。2. 准备一个算例，如水准测量内业计算，点号距离（km）观测高差(m) 改正数(mm) 改正后高差(m) 高程(m) 03 34.612 1.5 +8.237 +6 +8.243 N1 42.855 2.5 1.425 +9 1.416 N2 41.439 2.0 3.691 +8 3.683 04 37.756 6.0 +3.121 +23 +3.144 fh=h (H04 H03 )=+3.121 3.144= 0.023m fh容= 40 =40 98mm=0.098m 或准备一个经纬仪导线计算算例， 闭合或附合导线均可， 数据自备。四、实验指导L06.名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 17 页，共 29 页 - - - - - - - - - Excel 常用函数1、ABS函数函数名称： ABS 主要功能：求出相应数字的绝对值。使用格式： ABS(number) 参数说明： number代表需要求绝对值的数值或引用的单元格。应用举例：如果在B2单元格中输入公式： =ABS(A2)，则在 A2单元格中无论输入正数（如100）还是负数（如 -100） ，B2 中均显示出正数（如 100） 。特别提醒：如果number 参数不是数值，而是一些字符（如A等） ，则 B2中返回错误值“ #VALUE ！”。2、AND 函数函数名称： AND 主要功能：返回逻辑值：如果所有参数值均为逻辑“真（TRUE ）”，则返回逻辑“真（TRUE ）”，反之返回逻辑“假（FALSE ）”。使用格式： AND(logical1,logical2, .) 参数说明： Logical1,Logical2,Logical3：表示待测试的条件值或表达式，最多这30 个。应用举例：在 C5单元格输入公式： =AND(A5=60,B5=60) ，确认。如果 C5中返回 TRUE ，说明 A5和 B5中的数值均大于等于60，如果返回 FALSE ，说明 A5和 B5中的数值至少有一个小于60。特别提醒： 如果指定的逻辑条件参数中包含非逻辑值时，则函数返回错误值“ #VALUE! ”或“#NAME ”。3、AVERAGE 函数函数名称： AVERAGE 主要功能：求出所有参数的算术平均值。使用格式： AVERAGE(number1,number2, ) 参数说明： number1,number2,：需要求平均值的数值或引用单元格（区域），参数不超过 30 个。应用举例：在 B8单元格中输入公式：=AVERAGE(B7:D7,F7:H7,7,8)，确认后，即可求出B7至 D7区域、F7 至 H7区域中的数值和7、8 的平均值。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 18 页，共 29 页 - - - - - - - - - 特别提醒： 如果引用区域中包含“ 0”值单元格，则计算在内；如果引用区域中包含空白或字符单元格，则不计算在内。4、COLUMN 函数函数名称： COLUMN 主要功能：显示所引用单元格的列标号值。使用格式： COLUMN(reference) 参数说明： reference为引用的单元格。应用举例：在 C11单元格中输入公式： =COLUMN(B11)，确认后显示为 2（即 B列） 。特别提醒：如果在B11 单元格中输入公式： =COLUMN() ，也显示 出2 ； 与 之 相 对 应 的 还 有 一 个 返 回 行 标 号 值 的 函 数ROW(reference)。5、CONCATENATE函数函数名称： CONCATENATE 主要功能： 将多个字符文本或单元格中的数据连接在一起，显示在一个单元格中。使用格式： CONCATENATE(Text1，Text) 参数说明： Text1、Text2为需要连接的字符文本或引用的单元格。应用举例：在 C14单元格中输入公式：=CONCATENATE(A14,B14,.com) ，确认后，即可将 A14单元格中字符、 、B14单元格中的字符和 .com 连接成一个整体，显示在C14单元格中。特别提醒：如果参数不是引用的单元格，且为文本格式的，请给 参 数 加 上英 文 状 态 下的 双 引 号， 如 果 将上 述 公 式改 为 ：=A14&B14&.com ，也能达到相同的目的。6、COUNTIF 函数函数名称： COUNTIF 主要功能：统计某个单元格区域中符合指定条件的单元格数目。使用格式： COUNTIF(Range,Criteria) 参数说明： Range代表要统计的单元格区域； Criteria表示指定的条件表达式。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 19 页，共 29 页 - - - - - - - - - 应 用 举 例 ： 在C17单 元 格 中 输 入 公 式 ：=COUNTIF(B1:B13,=80)，确认后，即可统计出 B1至 B13单元格区域中，数值大于等于80 的单元格数目。特别提醒：允许引用的单元格区域中有空白单元格出现。7、DATE 函数函数名称： DATE 主要功能：给出指定数值的日期。使用格式： DATE(year,month,day) 参数说明： year 为指定的年份数值（小于9999） ；month为指定的月份数值（可以大于12） ；day 为指定的天数。五、思考题更为复杂的计算问题如何用Excel 实现。实验七： Mat lab在测绘工程中的应用一、实验目的和要求初步掌握用 Mat lab 进行数据处理的方法和过程；初步掌握 Mat lab 软件设计方法。二、实验安排本实验需要 24学时。实验每个小组1 人，用 Mat lab 在计算机上实现各种测量数据处理问题，学会开发简单的测绘程序。三、实验步骤及要点准备算例；用 Matlab 在计算机上实现数据处理；开发简单的测绘程序，调试运行，用实例印证。4. 用 Math lab 完成实验一中的方程组结算工作。四、实验指导MATLAB 的主要线性代数运算如表 A1 所示为常用的矩阵和线性代数运算函数，用户可以用 help matfun获得更多内容。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 20 页，共 29 页 - - - - - - - - - 表 A1 常用线性代数函数B=A矩阵转置C=A+B矩阵相加C=A*B矩阵相乘C=Ak矩阵幂C=A.*B矩阵点乘，即两维数相同的矩阵各对应元素相乘expm(A)指数矩阵，也就是eAinv(A)矩阵的逆矩阵det(A)矩阵的行列式的值rank(A)计算矩阵的秩eig(A)矩阵的特征值X,D=eig(A)矩阵的特征向量X 和以特征值为元素的对角阵Dp=poly(A)矩阵的特征多项式r=roots(p)特征多项式方程的根名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 21 页，共 29 页 - - - - - - - - - conv(p1,p2)两多项式相乘 常用的控制系统处理函数1） TF2SS 传递函数转换到状态空间表达式A,B,C,D = TF2SS(NUM,DEN) 将系统：转换成：其中: NUM=bm,bm-1, ,b1,b0 ，DEN=1,an-1,an-2, ,a1,a0 2） ZP2SS 零极点型传递函数转换到状态空间表达式A,B,C,D = ZP2SS(Z,P,K) 除了以外，其它与 TF2SS 相同。3）SS2TF 状态空间表达式转换到传递函数NUM,DEN = SS2TF(A,B,C,D,iu) 即求第 iu 个输入信号对输出y(t) 的传递函数，即：名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 22 页，共 29 页 - - - - - - - - - SS2TF 的调用返回值为 W(S)的分子多项式的系数矩阵NUM 和分母多项式的系数向量 DEN 。4）SS2ZP . 状态空间表达式到零极点形式的传递函数的转换Z,P,K = SS2ZP(A,B,C,D,iu) 5）TF2ZP 一般传递函数转换到零极点型传递函数Z,P,K = TF2ZP （NUM,DEN) 6）ZP2TF 零极点型传递函数转换到一般传递函数NUM,DEN = ZP2TF(Z,P,K) 7）SS2SS 状态空间表达式的线性变换A1,B1,C1,D1=ss2ss(A,B,C,D,T) 其中 T为变换矩阵。注意变换方程为：X1=TX ，而不是常见的 X=TX1 。所以要与用户习惯的变换方程一致，我们必须用T 的逆代入上式，即：A1,B1,C1,D1=ss2ss(A,B,C,D,inv(T) 8）CANON 求状态空间表达式的对角标准型As,Bs,Cs,Ds,Ts=canon(A,B,C,D,mod) 其中 Ts为变换矩阵，注意变换方程为：Xs=TsX 。9）CTRB 计算系统的可控判别矩阵M M=ctrb(A,B) 10）OBSV 计算系统的可观判别矩阵N N=obsv(A,C) 11）IMPULSE 求系统的单位脉冲响应y,x=impulse(A,B,C,D,in,t) y,x=impulse(num,den,t) 12）STEP 求系统的单位阶跃响应y,x=step(A,B,C,D,in,t) y,x=step(num,den,t) 13）LSIM 求系统对任意输入函数u(t) 的响应y,x=lsim(A,B,C,D,u,t) y,x=lsim(num,den,u,t) 14）C2D 连续系统状态方程转换为离散状态方程,T 为采样周期G,H=c2d(A,B,T) 相关的函数还有 D2C,D2D 。15 LYAP 求解如下形式的李雅普诺夫方程：名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 23 页，共 29 页 - - - - - - - - - 所以求解用户习惯的李雅普诺夫方程，我们必须用 A的转置A代入，即：P=lyap(A,Q) 16）BODE 绘制系统的伯德图mag,phase=bode(num,den,w) 17）NIQUIST 绘制系统的奈奎斯特图re,im=nyquist(num,den,w) 18）PLOT 绘图函数plot(X,Y,str ) 可以用不同颜色、不同符号绘制曲线，其中str 可以是下列三组选项的任意组合。y -yellow . -point - -solidm -magentao-circle: -dottedc -cyanx-x-mark-. -dashdot r -red + -plus- -dashed g -green* -starb -blues -squarew -white d -diamondk -black v -triangle (down) -triangle (up) -triangle (right)p -pentagramh -hexagram五、思考题 如何用 Matlab 开发测绘软件。 如何实现 Matlab 与 C语言通讯。实验八：测绘应用程序设计一、实验目的与要求1对以上理论内容的验证与应用。2通过学习掌握测绘软件开发过程与方法，初步具备测绘软件开发基本技能。二、实验安排1实验时数 12 学时。2每实验小组可以由34人组成，或独立完成。若由几个人完成程序设计，应进行合理的分工。三、实验步骤和要点1熟悉程序设计任务书的基本内容，调查了解软件需求状况，进行需求分析；2进行总体设计。根据所调查收集的资料和任务书的要求，对系统的硬件资源进行初步设计，提出硬件配置计划；进行软件总体设计，设计出软件程序功能的模块；3根据总体设计的结果，进行详细设计，进行数据存储格式设计、算法等，写出逻辑代码；4编写程序代码，调试运行；5程序试运行。最后同学们可根据自己的选题，写出软件开发名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 25 页，共 29 页 - - - - - - - - - 设计书一份，打印程序代码和运行结果。四、实验指导1同学们可从下列软件设计题目中人选一题进行开发设计，也可以结合实际自己拟题。2题目如下：序号题目要求1 导线严密平差程序设计适合于三、四等导线或一级导线测量平差计算，按严密平差算法编制程序，界面美观，最好用可视化程序设计。2 三四等水准测量平差计算采用严密平差方法，包括水准测量记录及计算、水准测量概算、平差计算等内容。并用程序验证图3 测边网平差计算程序设计适合于工程三、四等测边网的平差计算，要求界面美观，数据处理快、准，精度满足规范要求。4 摄影测量空间后交- 前交方法要求可以利用一定的控制点和像点，实现空间前交和后交法，计算出像片的外方位元素、待定点的地面坐标等（已知数据见下）5 水准网间接平差可以选择闭合或者复合书准网，按照条件平差/ 间接平差的原理进行程序设计。6 水准网条件平差7 道路曲线放样计算可以进行道路圆曲线、缓和曲线的的必要放样元素的计算。输出中庄桩号、转角、坐标等。7 沉降观测数据处理程序设计要求能够进行数据录入、数据计算、沉降分析、绘制下沉曲线图、形成沉降观测分析报告名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 26 页，共 29 页 - - - - - - - - - 等功能，界面美观，可视化，程序实用。3每人上交一份程序设计书并附运行测试结果和程序代码。报告按教师指定格式撰写。摄影测量空间后交 -前交的计算数据实例：四个角控制点地面坐标及对应的像点坐标：点号控制点像点（左片）X Y Z x y 1 36589.41 25273、 32 2195.17 -86.15 -68.99 2 37631.08 31324、 51 728、69 -53.40 82.21 3 39100.97 24934.98 2386.50 -14.78 -76.63 4 40426.54 30319.81 757.31 10.46 64.43 航摄比例尺：1：50000 航高： 7583 像片(左)的内方位元素： 0 ，0 ，153.24 待定点的像点坐标点号左片右片x y x y 1 4363 960 1660 796 2 6920 944 4273 876 3 6676 2265 3985 2201 4 7685 4860 4944 4857 5 5979 4449 3214 4387 6 3936 4040 1129 3903 7 3544 5239 696 5104 8 6233 6069 3426 6033 9 5135 7816 2239 7760 10 4656 6267 1810 6180 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 27 页，共 29 页 - - - - - - - - - 五、思考题1复杂的或任意图形的控制网如何用计算机实现平差计算，研究其数据结构存储格式。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 28 页，共 29 页 - - - - - - - - - - 1 -附件一：封面华北科技学院实 验 报 告课程名称：测绘程序设计专业：测绘工程班级学号：姓名：任课教师：名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 29 页，共 29 页 - - - - - - - - -