资源描述：
《大地测量学课程设计》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、大地测量学课程设计设计题目：白塞尔大地主题解算学院：矿业学院专业：测绘工程班级：********学号：120801024*目录1.基本原理及思想12.白塞尔法大地主题正算步骤23.白塞尔法大地主题反算步骤44.同一平行圈弧长、子午线弧长与大地线比较大小65.程序代码86.演算示例137.参考文献168.心得体会179.教师评语183白塞尔大地主题解算一：基本原理建立以椭球中心为中心，以任意长（或单位长）为半径的辅助球，按以下三个步骤计算。第一， 按一定条件将椭球面元素投影到辅助球面上。第二， 在球面
2、上解算大地问题。第三， 将求得的球面元素按投影关系换算到相应的椭球元素。关键：确定球面元素与椭球面元素的关系，即它们间的投影关系。二：白塞尔法解算大地主题的基本思想:以辅助球面为基础,将椭球面三角形转换为辅助球面的相应三角形,由三角形对应元素关系,将椭球面上的大地元素按照白塞尔投影条件投影到辅助球面上，然后在球面上进行大地主题解算，最后再将球面上的计算结果换算到椭球面上。三：在球面上进行大地主题解算球面上的大地主题正算：球面上的大地主题反算：18球面三角元素间的相互关系：四：正反算步骤1.白塞尔法大地
3、主题正算步骤已知、、、、()、S，计算、、()。（1）将椭球面元素投影到球面上由求：计算辅助量和,计算球面长度，将S化为18式中系数分别为：上式右端含有，因此需要迭代计算。第一次迭代取近似值，第二次计算取以后计算用代换代入上式迭代计算，直到所要求的精度为止。一般取。（2）解算球面三角形计算计算或计算（3）将球面元素换算到椭球面上由求18或将球面经差化为椭球面经差l，求l式中式中的最大值为，故在计算时通常可以略去不计。象限的判定符号++--符号+-+--l符号--++符号+-+-其中、为锐角。2.白塞尔法
4、大地主题反算步骤已知、、、、、，计算()、()、S。（1）将椭球面元素投影到球面上18由B求u，，l=,,采用逐次趋近法，由l计算在反算中，已知椭球面上经差l，球面经差上的对应经差未知，为了由l求，由下式可知还需计算、、，计算又还需量，故需要进行迭代计算。第一次趋近，取l；，或判断的象限p符号++--q符号+-+-判断象限+-18l+仿照上述计算步骤迭代计算，直到为止。（2）将球面元素换算到椭球面上或象限的判断与前面一致五：同一平行圈弧长、子午线弧长与大地线比较大小子午线弧长计算公式：式中：平行圈弧长公
5、式：（）不同纬度对应的一些弧长的数值B子午线弧长平行圈弧长1′1″1′1″0°110576m1842.94m30.716m111321m1855.36m30.923m15°1106561844.2630.7381075521792.5429.87630°1108631847.7130.795964881608.1326.8021845°1111431852.3930.873788481314.1421.90260°1114231857.0430.95155801930.0215.575°11162518
6、60.4231.00728902481.718.02890°1116961861.631.027000利用白塞尔大地主题反算求解大地线长S纬度为30°，经差为1°的平行圈弧长S=96488m，两点间大地线长为96487.595m经度为30°，纬度差为1°的子午线弧长X=110863m，两点间大地线长为110862.869m通过比较可知，同一平行圈或同一子午线两点间大地线长度与对应的平行圈弧长或子午线弧长相等。18六：程序代码#include#includedoubleh
7、udu(double,double,double);/*度分秒转换为弧度*/doubledu(double);/*弧度转换为度*/doublefen(double);/*弧度转换为分*/doublemiao(double);/*弧度转换为秒*/#definePI3.1415926voidmain(void){intk;printf("请选择大地主题算法，若执行正算，请输入1；若执行反算，请输入2。");scanf("%d",&k);/*大地主题正算*/if(k==1){doubleax,ay,az,
8、bx,by,bz,cx,cy,cz,S,dz,ez,fz,B1,B2,L1,L2,A1,A2;intdx,dy,ex,ey,fx,fy;doublee2,W1,sinu1,cosu1,sinA0,coto1,sin2o1,cos2o1,sin2o,cos2o,A,B,C,r,t,o0,o,g,sinu2,q;/*输入度分秒数据*/printf("请输入大地线起点纬度度分秒");scanf("%lf%lf%lf",&ax,&ay,&az)
