大地测量实验报告.docx

袁大地测量实验报告蚈莆节罿膈肇芄姓名:孙志勇莁学号: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); 肇肄/*计算球面长度*/薄 a0=(S-(B+C*cos2a1)*sin2a1)/A;蚀 m=sin2a1*cos(