平差程序设计剖析.docxVIP

中国地质大学（武汉）信息工程学院-《数据结构》课程实习报告 2015年秋 PAGE \* MERGEFORMAT- 19 - 115141-01 姓名实习三二叉树及其应用 PAGE \* MERGEFORMAT- 20 - + 实习题目 【题目一】角度、弧度、方位角 【问题描述】 （1）从键盘输入一个角度，判断其合法性； （2）从键盘输入一个角度值将其转化为对应的弧度值输出； （3）从键盘输入一个弧度值将其转化为对应的角度值输出； （4）从键盘输入两个点的坐标的值（x1,y1）（x2,y2），求方位角并输出。 【算法分析与设计】 问题分析与设计思想 题目一包含了四个部分，因此可以定义四个函数实现上述功能。 函数（一）主要应用if-else的分支结构来实现，首先是取出输入值的度位判断其是否位于[0，360)，接下来取其分位判断是否位于[0,60)，最后判断秒位是否位于[0,60)，这样就可以实现输入角度合法性判断的功能了。 函数（二）的设计就比较简单了，主要是要先取出输入的角度的度、分、秒对应的值a1,a2,a3，然后利用利用数学公式a1*PI/180+a2*PI/(180*60)+a3*PI/(180*60*60)式中（PI=3.1415926）；输出求得的结果即可。 函数（三）设计思路和函数二的整体思路基本一致，只不过首先要先利用数学公式a*180/PI求出输入弧度对应的以度为单位的值，然后就是取出改值的度、分、秒对应的值a1,a2,a3，最后用公式a1+a2/100+a3/10000；求得我们平差中常用的角度值。 函数（四）整体设计思路是运用反三角函数atan，但是要考虑一些细节问题，因为求得的(y2-y1)/(x2-x1)值可能为负数为了后面的计算方便要先调用一个求绝对值的函数fabs()，然后调用atan()函数求得的结果是弧度值，因此要调用之前的弧度转环为度的函数，接下来运用if-else if-else的分支语句分不同的象限求值，最后输出求得的结果。 出现的问题及解决办法 第一次设计函数四的功能时，在调用函数atan（）是由于疏忽是用下面这种方式调用的atan(y2-y1)/(x2-x1)，细心观测会发现少了一个括号，因此导致输出的结果不正确，然后就加断点分步运行调试最终找出问题所在。 3效率分析与改进设想 因为四个函数的代码较短，函数功能比较简单因此效率很高，暂时还没想出改进设想，期待日后更加深入的学习。 【源代码】 见附件。// HzJ.cpp : 定义控制台应用程序的入口点。 // #include stdafx.h #includeiostream #includecmath using namespace std; const double PI=3.1415926; void PD_angle(double a){ //传入一个角度 double c; int b,d; b=(int)a; if(b0||b360) //判断度位是否合法 cout输入角度不合法endl; else { c=(a-b)*100; d=(int)c; if(d=60) //判断分位是否合法 cout输入角度不合法endl; else{ c=(c-d)*100; if(c=60) //判断秒位是否合法 cout输入角度不合法endl; else cout输入角度合法endl; } } } double JzH(int a){ double b,d; int c1,c2,c3; c1=(int)a; //取出来度位 b=(a-c1)*100; c2=(int)b; //取出来分位 c3=(int)(b-c2)*100; //取出来秒位 d=c1*PI/180+c2*PI/(180*60)+c3*PI/(180*60*60); //计算弧度 cout输入的角度转为弧度是dendl; return d; } double HzJ(double a){ double c1,c2,c3,d; int b1,b2,b3; c1=a*180/PI; //将弧度转化为以度为单位 b1=(int)c1; //取出度位 c2=(c1-b1)*60; b2=(int)c2; //取出分位