平差程序设计

+

实习题目

【题目一】角度、弧度、方位角

【问题描述】

（1）从键盘输入一个角度，判断其合法性；

（2）从键盘输入一个角度值将其转化为对应的弧度值输出；

（3）从键盘输入一个弧度值将其转化为对应的角度值输出；

（4）从键盘输入两个点的坐标的值（x1,y1）（x2,y2），求方位角并输出。

【算法分析与设计】

1、问题分析与设计思想

题目一包含了四个部分，因此可以定义四个函数实现上述功能。

函数（一）主要应用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的分支语句分不同的象限求值，最后输出求得的结果。

2、出现的问题及解决办法

第一次设计函数四的功能时，在调用函数atan（）是由于疏忽是用下面这种方式调用的

atan(y2-y1)/(x2-x1)，细心观测会发现少了一个括号，因此导致输出的结果不正确，然后就加断点分步运行调试最终找出问题所在。

3效率分析与改进设想

因为四个函数的代码较短，函数功能比较简单因此效率很高，暂时还没想出改进设想，期待日后更加深入的学习。

【源代码】

见附件。// HzJ.cpp : 定义控制台应用程序的入口点。//

#include"stdafx.h"

#include

#include

usingnamespace std;

constdouble PI=3.1415926;

void PD_angle(double a){ //传入一个角度

double c;

int b,d;

b=(int)a;

if(b<0||b>360) //判断度位是否合法 cout<<"输入角度不合法<

else

{

c=(a-b)*100;

d=(int)c;

if(d>=60) //判断分位是否合法 cout<<"输入角度不合法<

else{

c=(c-d)*100;

if(c>=60) //判断秒位是否合法 cout<<"输入角度不合法<

else

cout<<"输入角度合法<

}

}

}

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<<"输入的角度转为弧度是"<

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; //取出分位

c3=(c2-b2)*60;

b3=(int)c3; //取出秒位

d=b1+b2/100+b3/10000;

cout<<"输入的弧度转为角度为a："<

return d;

}

double FwJ(double x1,double y1,double x2,double y2){ double a,b,d,x,y;

if(x1==x2 &&y1*y2>0){

d=90;

cout<<"两点间的方位角为a："<

}

elseif(x1=x2 && y1*y2<0){

d=270;

cout<<"两点间的方位角为a："<

}

else{

x=x2-x1;

y=y2-y1;

a=atan(y/x);

a=(a>0)? a:-a;

d=HzJ(a);

//x=x2-x1;

//y=y2-y1;

if(x>0&&y>0)

d=d;

elseif(x<0&&y>0)

d=180-d;

elseif(x<0&&y<0)

d=180+d;

else

d=360-d;