大地主题解算.docVIP

大地主题解算 实验目的： 提高运用计算机语言编程开发的能力； 加深对大地主题解算计算公式及辅助参数的理解并掌握计算步骤； 通过编程语言实现大地主题解算。 工具： Windows XP Mode 环境下的Microsoft Visual C++ 6.0 注意事项： 计算所需变量多，容易混淆； 正反算函数的编写； 函数调用； 弧度与角度之间的转化。 实验要求： 提交报告，实验总结，编写代码； 独立编程，调试运行； 上交成果：编写思想，编写过程，问题分析，源代码，计算结果； 编程过程实现： 对白塞尔法大地主题解算有一定的了解，并参考教材P148-P150； 由于参数较多，而在C语言环境下很多符号无法定义，需要符合要求的定义符号替代书本上那些无法直接在C语言环境下定义的符号来达到实现实验的目的； 程序中采用弧度与度分秒之间转换的函数定义与调用，减轻一定的实验麻烦； 在C语言环境下，数学函数fabs代替abs起绝对值作用，atan代替arctan起反函数作用； 程序中尤其注意弧度与角度之间转换，在C语言环境下电脑默认为弧度。 六、源程序代码： #includestdio.h #includemath.h double hudu(double,double,double); /*度分秒转换为弧度*/ double du(double); /*弧度转换为度*/ double fen(double); /*弧度转换为分*/ double miao(double); /*弧度转换为秒*/ #define PI 3.1415926 void main (void) { int k; printf(请选择执行正算或者反算，若执行正算，请输入1；若执行反算，请输入2。\n); scanf(%d,k); /*正算*/ if(k==1) { double bz,lz,az,S,bz2,lz2,az2,B1,L1,A1,B2,L2,A2,bx,by,lx,ly,ax,ay; int bx2,by2,lx2,ly2,ax2,ay2; double e2,W1,sinu1,cosu1,sinA0,coto1,sin2o1,cos2o1,sin2o,cos2o,A,B,C,r,t,o0,o,g,sinu2,q; /*以度分秒顺序输入数据*/ printf(请输入大地线起点纬度度分秒\n); scanf(%lf%lf%lf,bx,by,bz); printf(请输入大地线起点经度度分秒\n); scanf(%lf%lf%lf,lx,ly,lz); printf(请输入大地方位角度分秒\n); scanf(%lf%lf%lf,ax,ay,az); printf(请输入大地线长度\n); scanf(%lf,S); /*调用函数*/ B1=hudu(bx,by,bz); L1=hudu(lx,ly,lz); A1=hudu(ax,ay,az); /*白塞尔大地主题解算*/ e2=0.006693421622966; W1=sqrt(1-e2*sin(B1)*sin(B1)); sinu1=sin(B1)*(sqrt(1-e2))/W1; cosu1=cos(B1)/W1; sinA0=cosu1*sin(A1); coto1=cosu1*cos(A1)/sinu1; sin2o1=2*coto1/(coto1*coto1+1); cos2o1=(coto1*coto1-1)/(coto1*coto1+1); A=6356863.020+(10718.949-13.474*(1-sinA0*sinA0))*(1-sinA0*sinA0); B=(5354.469-8.978*(1-sinA0*sinA0))*(1-sinA0*sinA0); C=(2.238*(1-sinA0*sinA0))*(1-sinA0*sinA0)+0.006; r=691.46768-(0.58143-0.00144*(1-sinA0*sinA0))*(1-sinA0*sinA0); t=(0.2907-0.0010*(1-sinA0*sinA0))*(1-sinA0*sinA0); o0=(S-(B+C*cos2o1)*sin2o1)/A; sin2o=sin2o1*cos(2*o0)+cos2o1*sin(2*o0); cos2o=cos2o1*cos(2*o0)