文档介绍：Forpersonaluseonlyinstudyandresearch;mercialuse薆大地测量实验报告薁羁薆蚆羂荿虿姓名:孙志勇螆学号:20104220莃专业班级:测绘工程10级(2)膀莈袆螃薈膆袆袀芀羅羆芁螈实验名称:白塞尔大地主题解算(正算和反算)羈实验目的:。,继而在球面上进行大地主题解算,最后再将球面上的计算结果换算到椭球面上的基本方法和步骤。。蒀实验环境:MicrosoftVisualC++螇注意事项:。。肃实验步骤:正算:,B,C及d,。蕿反算:、球面长度及经差。,B,。莂程序源代码:芈正算:#include<>蒅#include<>肂##(double,double,double);蒅voidmain(void)蒂{薁 doubleA1,B1,L1,S,A2,B2,L2;袅 doublex1,x2,x3,y1,y2,y3,z1,z2,z3;薅 doubleW1,sinu1,sinu2,cosu1,sinA0;袃 doublecota1,cos2a1,sin2a1,cosA0A0;罿 doubleA,B,C,d,e,a0,a1,m;袈 doublen,a,Q,R;蚅 printf("请输入数据B1=");羀 scanf("%lf%lf%lf",&x1,&x2,&x3);蚁 B1=F(x1,x2,x3);蚇 printf("请输入数据L1=");螄 scanf("%lf%lf%lf",&y1,&y2,&y3);莁 L1=F(y1,y2,y3);腿 printf("请输入A1=");蒆 scanf("%lf%lf%lf",&z1,&z2,&z3);袄A1=F(z1,z2,z3);螂 printf("请输入S=");袁 scanf("%lf",&S);葿printf("B1=%f\n",B1);羄 printf("L1=%f\n",L1);膃 printf("A1=%f\n",A1);艿 printf("S=%f\n",S);芈肄/*计算起点的规划纬度*/薄 W1=sqrt(1-ee*sin(B1)*sin(B1));肁 sinu1=sin(B1)*sqrt(1-ee)/W1;羇 cosu1=cos(B1)/W1;肄螁 蒈 printf("W1=%f\n",W1);螆 printf("sinu1=%f\n",sinu1);膄 printf("cosu1=%f\n",cosu1);膁膀/*计算辅助函数值*/螈 sinA0=cosu1*sin(A1);芄 cota1=cosu1*cos(A1)/sinu1;薂 sin2a1=2*cota1/(cota1*cota1+1);蚈 cos2a1=(cota1*cota1-1)/(cota1*cota1+1);薇 莄 printf("sinA0=%f\n",sinA0);羃 printf("cota1=%f\n",cota1);莀 printf("sin2a1=%f\n",sin2a1);莆 printf("cos2a1=%f\n",cos2a1);蒃莄/*计算系数ABC及de*/袈 cosA0A0=1-sinA0*sinA0;荿 A=+(-(*cosA0A0))*cosA0A0;薃 B=(-*cosA0A0)*cosA0A0;蒁 C=(*(cosA0A0))*cosA0A0+;薀 d=-(-*cosA0A0)*cosA0A0;膈 e=(-cosA0A0*)*cosA0A0;薃袂 printf("cosA0A0=%f\n",cosA0A0);节 printf("A=%f\n",A);袇 printf("B=%f\n",B);蚃 printf("C=%f\n",C);芃 printf("d=%f\n",d);螀 printf("e=%f\n",e); 蚆螃/*计算球面长度*