高斯平均引数大地主题正反算

地球科学与环境工程学院实验报告书
课程名：
学号：
姓名：
指导老师：
日期：
目录
一、目的与要求 (1)
二、实验内容 (1)
三、计算公式整理 (1)
四、程序代码 (4)
五、计算结果 (15)
六、实验体会 (16)
一、目的与要求
参考椭球面是大地测量计算的基准面。

大地坐标是椭球面上的基本坐标系，根据大地测量的观测成果（如距离与方向），从大地原点出发，逐点计算在椭球面上的大地坐标；或根据两点的大地坐标，计算它们之间的大地线长度和大地方位角，这类计算称为大地问题解算（或称为大地主题解算）。

大地问题解算的用途是多方面的，随着现代空间技术和航空航天、航海等领域的发展，大地问题解算（尤其是大地反算）有着更为重要的作用，因此需要熟练掌握其计算。

二、实验内容
在《大地测量学基础》教材中，介绍了高斯平均引数法与白塞尔方法的计算过程、步骤。

鉴于此，需要熟练掌握高斯平均引数法与白塞尔方法解大地主题问题的基本方法与原理。

采用所熟悉的计算机语言编程计算。

计算时采用克拉索夫椭球参数，至少完成其中一种方法正反算，按照数据序号选取不同的已知数据，在计算结果中注明所选取的数据序号，选取其它数据作为无效数据处理。

三、计算公式整理
、高斯平均引数正算计算公式（S< 200 km）
、高斯平均引数正算计算公式（S< 200 km）
四、程序代码
、角度转换类的头文件：
#pragma once
const double Pi=;
class AngleTrans
{
public:
AngleTrans(void);
~AngleTrans(void);
double D,F,M,DFM,Rad,Ten;
double trans1(double DFM), f"),B2);
(_T("%.7f"),L2),
(_T("%.7f"),A21);
UpdateData(false);
}
、反算的计算按钮代码：
void C大地主题高斯引数正反算Dlg::OnBnClickedButton1() { UpdateData(true);
ZhengFanSuan _ZhengFanSuan;
double B1=_wtof(fB1),
L1=_wtof(fL1),
B2=_wtof(fB2),
L2=_wtof(fL2),
S=(B1,L1,B2,L2),
A12=(B1,L1,B2,L2),
A21=(B1,L1,B2,L2);
(_T("%.7f"),S);
(_T("%.7f"),A12);
(_T("%.7f"),A21);
UpdateData(false);
}
、清零按钮代码：
void C大地主题高斯引数正反算Dlg::OnBnClickedButton3() { UpdateData(true);
(_T("%.0f"),0);
(_T("%.0f"),0);
(_T("%.0f"),0);
(_T("%.0f"),0);
(_T("%.0f"),0);
(_T("%.0f"),0);
(_T("%.0f"),0);
UpdateData(false);
}
void C大地主题高斯引数正反算Dlg::OnBnClickedButton4()
{ UpdateData(true);
(_T("%.0f"),0);
(_T("%.0f"),0);
(_T("%.0f"),0);
(_T("%.0f"),0);
(_T("%.0f"),0);
(_T("%.0f"),0);
(_T("%.0f"),0);
UpdateData(false);
}
五、计算结果
数据组号：18
六、实习体会
实验一开始本来考虑用控制台做，后来发现用控制台做出来可视化效果不好，而且程序中要用很多cin和cout，写代码十分麻烦，所以最后选择了用MFC做。

其实技术上的难题并不多，主要在于公式上面：大量的乘法、开方和乘方，一不小心就会出错，而且在检查的过程中也不容易发现；在显示计算结果时出现了几次奇奇怪怪的数据，后来在检查代码时发现是由于某些原因（没经历循环，而目标值必须从循环过程中获得、if判断时漏掉情况等）导致要显示的变量没有被赋值，才会无法显示出需要的结果；最初没有加上清零按钮，在调试的过程中发现经常改数值会使得输入数据很麻烦，于是最后加上了清零按钮，以便使用。

总的来讲在实验过程中没有遇到什么大的问题，小的问题却是接二连三地来，好在最后都能够一一解决，也知道了一些C++中的一些报错是什么原因引起的，感觉收获不小。