《(vc++ mfc)高斯平均引数大地主题正反算》由会员分享，可在线阅读，更多相关《(vc++ mfc)高斯平均引数大地主题正反算（18页珍藏版）》请在金锄头文库上搜索。
1、1地球科学与环境工程学院实 验 报 告 书课 程 名 ： 学 号： 姓 名： 指导老师： 日 期： 1目 录一、目的与要求 .1二、实验内容 .1三、计算公式整理 .1四、 程序代码.4五、计算结果 .15六、实验体会 .161一、目的与要求参考椭球面是大地测量计算的基准面。大地坐标是椭球面上的基本坐标系，根据大地测量的观测成果（如距离与方向） ，从大地原点出发，逐点计算在椭球面上的大地坐标；或根据两点的大地坐标，计算它们之间的大地线长度和大地方位角，这类计算称为大地问题解算（或称为大地主题解算） 。大地问题解算的用途是多方面的，随着现代空间技术和航空航天、航海等领域的发展，大地问题解算（尤其
2、是大地反算）有着更为重要的作用，因此需要熟练掌握其计算。二、实验内容在大地测量学基础教材中，介绍了高斯平均引数法与白塞尔方法的计算过程、步骤。鉴于此，需要熟练掌握高斯平均引数法与白塞尔方法解大地主题问题的基本方法与原理。采用所熟悉的计算机语言编程计算。计算时采用克拉索夫椭球参数，至少完成其中一种方法正反算，按照数据序号选取不同的已知数据，在计算结果中注明所选取的数据序号，选取其它数据作为无效数据处理。三、计算公式整理3.1、高斯平均引数正算计算公式（SAngleTrans:AngleTrans(void)AngleTrans:AngleTrans(void)/度分秒转换为弧度double An
3、gleTrans:trans1(double DFM) D=floor(DFM);F=floor(DFM-D)*100);M=(DFM-D-F/100)*10000;Ten=D+F/60+M/3600;Rad=Ten/180*Pi;return Rad;5/弧度转换为度分秒double AngleTrans:trans2(double Rad) Ten=Rad/Pi*180;D=floor(Ten);F=(Ten-D)*60;M=(F-floor(F)*60;F=floor(F);DFM=D+F/100+M/10000;return DFM;/十进制度转化为弧度double AngleTran
4、s:trans3(double D) Rad=D/180*Pi;return Rad;4.3、正反算类的头文件：#pragma onceclass ZhengFanSuanpublic:ZhengFanSuan(void);ZhengFanSuan(void);double zB1,zL1,zA12,zS,fB1,fL1,fB2,fL2;double ZhengSuanB(double zB1,double zL1,double zA12,double zS), ZhengSuanL(double zB1,double zL1,double zA12,double zS),ZhengSuanA
5、(double zB1,double zL1,double zA12,double zS);double FanSuanA12(double fB1,double fL1,double fB2,double fL2),FanSuanS(double fB1,double fL1,double fB2,double fL2),FanSuanA21(double fB1,double fL1,double fB2,double fL2);4.3、正反算类的源文件：#include StdAfx.h#include ZhengFanSuan.h#includeAngleTrans.h#include
6、ZhengFanSuan:ZhengFanSuan(void)6ZhengFanSuan:ZhengFanSuan(void)AngleTrans _AngleTrans;const double e1= 0.0066934216622966, e2=0.006738525414683, a=6378245.0000, b= 6356863.01877,temp=pow(10.0, -10);/精度要求double Calc_M(double z) /计算Mm double x=a*(1-pow(e1,2)/sqrt(pow(1-pow(e1,2)*pow(sin(z),2),3);retur
7、n x;double Calc_N(double z) /计算Nm double x=a/sqrt(1-pow(e1,2)*pow(sin(z),2);return x;double Calc_t(double z) /计算tm double x=tan(z);return x;double Calc_yita(double z) /计算yitam double x=pow(e2,2)*pow(cos(z),2);return x;/正算纬度double ZhengFanSuan:ZhengSuanB(double zB1,double zL1,double zA12,double zS) d
8、ouble M10000,N10000,t10000,B10000,Bm10000,L10000,Lm10000,A10000,Am10000,yita10000;double _zB1=_AngleTrans.trans1(zB1),_zL1=_AngleTrans.trans1(zL1),_zA12=_AngleTrans.trans1(zA12),7final;t1=tan(_zB1);yita1=pow(e2,2)*pow(cos(_zB1),2);N1=a/sqrt(1-pow(e1,2)*pow(sin(_zB1),2);M1=a*(1-pow(e1,2)/sqrt(pow(1-p
9、ow(e1,2)*pow(sin(_zB1),2),3);B0=_AngleTrans.trans3(206265/M1*zS*cos(_zA12)/3600),L0=_AngleTrans.trans3(206265*zS*sin(_zA12)/(N1*cos(_zB1)/3600),A0=L0*sin(_zB1);Bm1=_zB1+1/2*B0;Lm1=_zL1+1/2*L0;Am1=_zA12+1/2*A0;int i=1;B1=_AngleTrans.trans3(206265/M1*zS*cos(Am1)*(1+pow(zS,2)/(24*pow(N1,2)*(pow(sin(Am1
10、),2)*(2+3*pow(t1,2)+2*pow(yita1,2)+3*pow(cos(Am1),2)*pow(yita1,2)*(pow(t1,2)-1-pow(yita1,2)-4*pow(yita1,2)*pow(t1,2)/3600);Li=_AngleTrans.trans3(206265/(N1*cos(Bm1)*zS*sin(Ami)*(1+pow(zS,2)/(24*pow(Ni,2)*(pow(ti,2)*pow(sin(Ami),2)-pow(cos(Ami),2)*(1+pow(yitai,2)-9*pow(yitai,2)*pow(ti,2)/3600);Ai=_An
11、gleTrans.trans3(206265/Ni*zS*sin(Ami)*ti*(1+pow(zS,2)/(24*pow(Ni,2)*(pow(cos(Ami),2)*(2+7*pow(yitai,2)+9*pow(yitai,2)*pow(ti,2)+5*pow(yitai,4)+pow(sin(Ami),2)*(2+pow(ti,2)+2*pow(yitai,2)/3600);while(Bi-Bi-1=temp)Bi+1=_AngleTrans.trans3(206265/Mi*zS*cos(Ami)*(1+pow(zS,2)/(24*pow(Ni,2)*(pow(sin(Ami),2
12、)*(2+3*pow(ti,2)+2*pow(yitai,2)+3*pow(cos(Ami),2)*pow(yitai,2)*(pow(ti,2)-1-pow(yitai,2)-4*pow(yitai,2)*pow(ti,2)/3600);Li+1=_AngleTrans.trans3(206265/(Ni*cos(Bmi)*zS*sin(Ami)*(1+pow(zS,2)/(24*pow(Ni,2)*(pow(ti,2)*pow(sin(Ami),2)-pow(cos(Ami),2)*(1+pow(yitai,2)-9*pow(yitai,2)*pow(ti,2)/3600);Ai+1=_A
13、ngleTrans.trans3(206265/Ni*zS*sin(Ami)*ti*(1+pow(zS,2)/(24*pow(Ni,2)*(pow(cos(Ami),2)*(2+7*pow(yitai,2)+9*pow(yitai,2)*pow(ti,2)+5*pow(yitai,4)+pow(sin(Ami),2)*(2+pow(ti,2)+2*8pow(yitai,2)/3600);Bmi+1=_zB1+1/2*Bi;Lmi+1=_zL1+1/2*Li;Ami+1=_zA12+1/2*Ai;Mi+1=Calc_M(Bi);Ni+1=Calc_N(Ni);yitai+1=Calc_yita(y
