大地测量主题解算白赛尔大地主题解算
大地测量学白塞尔大地主题解算实习报告(模板)1
大地测量上机实习报告
一、实习时间：第9周星期三
二、实习地点：系办三楼机房
三、实习内容：白塞尔大地主题解算
四、实习目的：（1）尝试编写白塞尔主题解算的程序
（2）通过实例检验解算的结果
五、实习步骤：
1、先了解书本上关于白塞尔大地主题解算的步骤；
2、双击打开白塞尔大地主题解算的程序，如下：
3、在“选择椭球”中选择“克氏椭球”，“选择正反算”中选择“正算”；
4、然后在输入框中输入相应的数据，如下：
5、接下来单击菜单栏中的“数据正确”，记得正算结果，如下所示：
6、重复第3步，选择反算，并输入相应数据：
7、点击“数据正确”，即得反算结果：
六、实习总结：通过白塞尔大地主题结算方法的上机操作，懂得了白塞尔法结算大地主题是将椭球面上的大地元素按照白塞尔投影条件投影到辅助球面上，继而在球面上进行大地主题解算，再将球面上的计算结果换算到椭球面上的原理，并熟悉白塞尔大地主题解算的程序。
白塞尔大地主题算法
1．球面上大地主题正解方法
此时已知量：φ 1，α1及σ；要求量：φ 2，α2及λ. 首先按(i)式计算 ，继而用下式计算φ 2球面上坐标关系式
如图13—4所示，在椭球 面极三角形PP1P2中．用B． L．S及A分别表示大地线上某 点的大地坐标，大地线长及其 大地方位角。在球面极三角形 P'P1'P2'中，与之相应，用φ , λ，σ及α分别表示球面大圆弧 上相应点的坐标，弧长及方位 角。
M-子午圈曲率半径 N-卯酉圈曲率半径
M=N=1
为简化计算，白塞尔提出如下三个投影条件: 1 . 椭球面大地线投影到球面上为大圆弧； 2．大地线和大圆弧上相应点的方位角相等; 3．球面上任意一点的纬度等于椭球面上相应 点的归化纬度.
u-归化纬度，
r-平行圈半径，a-长半轴
（三）白塞尔微分方程的积分
(一)在球面上进行大地主题解算
如图13—3示，在球面上 有两点P1和P2，其中P1点的大 地纬度φ 1，大地经度λ1,P2点 大地纬度φ 2，大地经度λ2；P1 和P2点间的大圆弧长为σ， P1P2的方位角为α1，其反方位 角为α2，球面上大地主题正算 是已知φ 1，α1，σ，要求φ 2, α2及经差λ；反算问题是已知 φ 1,φ 2及经差λ，要求σ，α1及 α2。
（四）白塞尔法大地主题正算步骤
（五）白塞尔法大地主题反算步骤
白塞尔大地主题解算方法
白塞尔法解算大地主题的基本思想是将椭球面 上的大地元素按照白塞尔投影条件投影到辅助球 面上，继而在球面上进行大地主题解算，最后再 待球面上的计算结果换算到椭球面上。由此可见 ，这种方法的关键问题是找出椭球面上的大地元 素与球面上相应元素之间的关系式。同时也要解 决在球面上进行大地主题解算的方法。
白塞尔大地主题解算的基本思想
白塞尔大地主题解算的基本思想
首先，白塞尔大地主题解算的基本思想之一是建立地壳变形的力学模型。
地壳变形是地球表面的一项重要现象，是由于地质作用和地球内部力
学过程的结果。
地壳变形的力学模型是研究地壳变形的重要方法。
常见的
地壳变形的力学模型有弹性模型、弹塑性模型和拟弹性模型等。
这些模型
可以描述地壳的变形过程，通过对地壳的变形过程进行建模，可以更好地
理解地壳变形的机制和动力学过程。
其次，白塞尔大地主题解算的思想之一是研究现今地球表面的动力学
过程。
地球表面的动力学过程包括板块运动、地震活动、火山喷发等。
这
些过程在地球的长期演化中起着重要的作用。
通过对现今地球表面动力学
过程的研究，可以揭示地球内部的结构和动力学机制，进而更好地理解地
壳变形的成因和发展变化的规律。
第三，白塞尔大地主题解算的基本思想之一是研究地壳变形的成因。
地壳变形的成因是地壳运动的基本原因，也是地球科学研究的一个重要问题。
地壳变形的成因包括构造运动、地壳应力状态的改变、地震活动等。
通过研究地壳变形的成因，可以更好地了解地壳变形的机制和规律，进而
为地震预测和地壳运动的控制提供科学依据。
总之，白塞尔大地主题解算的基本思想是通过建立地壳变形的力学模型，研究现今地球表面的动力学过程，探究地壳变形的成因，揭示地壳变
形与地球动力学过程之间的相互关系，以推动地壳运动的理论和应用研究。
这一思想对于研究地壳运动的机制和规律，了解地壳变形的成因和动力学
过程，提高地震预测和地壳运动控制的水平具有重要意义。
白塞尔大地主题解算课程设计
白塞尔大地主题解算课程设计大地测量学课程设计设计题目:白塞尔大地主题解算学院: 矿业学院专业: 测绘工程班级: 测绘122学号: 1208010255学生姓名: 郑英智指导教师: 张俊2014年12月26日目录1. 白塞尔法大地主题正算步骤 12. 白塞尔法大地主题反算步骤 33. 同一平行圈弧长、子午线弧长与大地线比较大小 44. 程序代码 55. 演算示例 96. 参考文献 127. 心得体会 13教师评语 14白塞尔大地主题解算摘要:白塞尔法解算大地主题的基本思想是将椭球面上的大地元素按照白塞尔投影条件(,椭球面大地线投影到球面上为大圆弧;,大地线和大圆弧上相应点的方位角相等;,球面上任意一点的纬度等于椭球面上相应点的归化纬度。
)投影到辅助球面上，继而在球面上进行大地主题解算，最后再将球面上的计算结果换算到椭球面上。
1.白塞尔法大地主题正算步骤2已知、、、、()、S，计算、、()。
BLABLAaAAe1112221221(1)将椭球面元素投影到球面上2utanu,1,etanB求:,由B1111A,计算辅助量和 ,01sinA,cosusinAtan,,tanusecA , 011111, ,计算球面长度，将S化为,,,S，,sin,cos(2,，,)，,sin2,cos(4,，2,) 11式中系数分别为:''''246k7k15k,,,,(1,，,，？) ,bAb464256''246,Bkk37k'',,,,(,，，？) A48512''46,Ckk'',,,,(,，？) A12812822'2k,ecosA 024635kkk ,1，,,，？A464256246kk15kB,,，，？41651246kk3C,,？128512,,,,S上式右端含有，因此需要迭代计算。
第一次迭代取近似值，第二次计算取 0,,,S，,sin,cos(2,，,)，,sin2,cos(4,，2,) 01010- 1 -,以后计算用代换,代入上式迭代计算，直到所要求的精度为止。
大地主题解算深度干货超精
大地主题解算(正算)代码与白塞尔大地主题解算大地主题解算(正算)代码：根据经纬度和方向角以及距离计算另外一点坐标新建模块->拷贝下面的大地主题(正算)代码，调用方法示例：起点经度：116.235(度)终点纬度：37.435(度)方向角：50（度）长度：500（米）终点经纬度（"经度,纬度"）=Computation(37.435,116.235,50,500)Const Pi = 3.14Private a, b, c, alpha, e, e2, W, V As Double'a 长轴半径'b 短轴'c 极曲率半径'alpha 扁率'e 第一偏心率'e2 第二偏心率'W 第一基本纬度函数'V 第二基本纬度函数Private B1, L1, B2, L2 As Double'B1 点1的纬度'L1 点1的经度'B2 点1的纬度'L2 点2的经度Private S As Double '''''大地线长度Private A1, A2 As Double'A1 点1到点2的方位角'A2 点2到点1的方位角Function Computation(STARTLAT, STARTLONG, ANGLE1, DISTANCE As Double) As StringB1 = STARTLATL1 = STARTLONGA1 = ANGLE1S = DISTANCEa = 6378245b = 6356752.3142c = a ^ 2 / balpha = (a - b) / ae = Sqr(a ^ 2 - b ^ 2) / ae2 = Sqr(a ^ 2 - b ^ 2) / bB1 = rad(B1)L1 = rad(L1)A1 = rad(A1)W = Sqr(1 - e ^ 2 * (Sin(B1) ^ 2))V = W * (a / b)Dim W1 As DoubleE1 = e ''''第一偏心率'// 计算起点的归化纬度W1 = W ''Sqr(1 - e1 * e1 * Sin(B1 ) * Sin(B1 )) sinu1 = Sin(B1) * Sqr(1 - E1 * E1) / W1cosu1 = Cos(B1) / W1'// 计算辅助函数值sinA0 = cosu1 * Sin(A1)cotq1 = cosu1 * Cos(A1)sin2q1 = 2 * cotq1 / (cotq1 ^ 2 + 1)cos2q1 = (cotq1 ^ 2 - 1) / (cotq1 ^ 2 + 1)'// 计算系数AA,BB,CC及AAlpha, BBeta的值。
白塞尔大地主题正反算
地球椭球 CGCS2000 6378137 6356752.314 6399593.626 1/298.257222101 0.00669438 0.006739497 第一偏心率e2= 第一偏心率e‘2= sinu1= cotσ1= A0= B= β= cos2(σ1+σ0)= sinu2= 大地方位角A21= 第一偏心率e2= u1= u2= cos(2σ‘1+Δσ‘)= ω2 = ω5 = s2= s5= 大地方位角A21= u1= u2= cos(2σ‘1+Δσ‘)= ω2 = ω5 = s2= s5=
2
克拉索夫斯基椭球 BJ54 6378245 6356863.019 6399698.902 1/298.3 0.006693422 0.006738525 BJ54 35°00′00.000000″ 35°00′00.000000″ 0.820055451 1.9591E-09 0.002206928 0.241484809 0.002515577 0.002514192 34°59′59.544529″ BJ54 35°00′00.000000″ 35°00′00.000000″ 0.002508527 4.64788E-09 -4.15588E-09 0.003065877 9.19076E-10 16000.00241 34°59′59.544529″ 35°10′30.957713″ 0.002508527 4.64788E-09 -4.15588E-09 0.003065877 9.19076E-10
IUGG--1979椭球 WGS84 6378137 6356752.314 6399593.626 1/298.257223563 0.00669438 0.006739497 6356863.019 89°59′59.999859″ 16000 0.820055451 -1 6360368.854 0.003351407 -0.005031132 6.90795E-06 35°10′30.957713″ 6378245 34°59′59.544529″ 35°10′30.957713″ 0.820055451 0.003358976 9.24261E-07 0.999451116 -0.000548959 90°00′00.000141″ 35°00′00.000000″ 35°00′00.000000″ 0.820055451 0.003358976 9.24261E-07 0.999451116 -0.000548959
大地测量学复习资料#(精选.)
1.垂线偏差：地面一点上的重力向量g和相应椭球面上的法线向量n之间的夹角定义为该点的垂线偏差。
2.参考椭球：具有确定参数（长半径a和扁率α），经过局部定位和定向，同某一地区大地水准面最佳拟合的地球椭球，叫参考椭球。
3.大地线：椭球面上两点间的最短程曲线叫做大地线。
4.力高：水准面在纬度45度处的正常高。
5.大地主题解算：已知某些大地元素推求另一些大地元素的计算工作叫大地主题解算。
6.大地主题正算：已知P１点的大地坐标(L１，Ｂ１），Ｐ１至P２的大地线长S及其大地方位角，计算P２点的大地坐标(L２，Ｂ２）和大地线S在P２点的反方位角Ａ２１，这类问题叫做大地主题正算。
7.大地基准：是指能够最佳拟合地球形状的地球椭球的参数及椭球定位和定向8.高斯投影：横轴椭圆柱等角投影（假象有一个椭圆柱横套在地球椭球体外，并与某一条子午线相切，椭球柱的中心轴通过椭球体中心，然后用一定投影方法，将中央子午线两侧各一定范围内的地区投影到椭圆柱上，再将此柱面展开成投影面）。
9.大地测量学：是指在一定的时间与空间参考系中，测量和描绘地球形状及其重力场并监测其变化，为人类活动提供关于地球的空间信息的一门科学。
10.理论闭合差：由水准面不平行而引起的水准环线闭合差，称为理论闭合差。
11.地心坐标系：地心坐标系是在大地体内建立的O-XYZ坐标系。
原点O设在大地体的质量中心，用相互垂直的X，Y，Z三个轴来表示，X轴与首子午面与赤道面的交线重合，向东为正。
Z轴与地球旋转轴重合，向北为正。
Y 轴与XOZ平面垂直构成右手系。
12.高斯投影正、反算公式进行换带计算的步骤。
这种方法的实质是把椭球面上的大地坐标作为过度坐标。
首先把某投影带内有关点的平面坐标（x，y）1利用高斯投影反算公式换算成椭球面上的大地坐标（B，l），进而得到L=L0+l，然后再由大地坐标(B,l)，利用投影正算公式换算成相邻带的平面坐标（x，y）2在计算时，要根据第2带的中央子午线来计算经差l，亦即此时l=L-L0。
白塞尔大地主题解算(正算和反算)
大地测量实验报告实验名称：白塞尔大地主题解算（正算和反算）实验目的：1.通过编写白塞尔大地主题电算程序进一步掌握白塞尔法解算大地主题的基本思想。
2.熟练掌握将椭球面上的大地元素按照白塞尔投影条件投影到辅助球面上，继而在球面上进行大地主题解算，最后再将球面上的计算结果换算到椭球面上的基本方法和步骤。
3.学会掌握计算机编程的基本能力。
实验环境：Microsoft Visual C++注意事项：1.在编写程序的过程当中要注意代码的前后统一和重复。
2.注意数值类型的转换和度分秒的换算。
实验步骤：正算：1.计算起点的规划纬度2.计算辅助函数值3.计算系数A,B,C及d,e.4计算球面长度5.计算经度差改正数6.计算终点大地坐标及大地方位角。
反算:1.进行计算辅助函数值2.用逐次趋近法同时计算起点大地方位角、球面长度及经差。
3.计算系数A,B,C及大地线长度S.4.计算反方位角及确定符号。
程序源代码：正算：#include<stdio.h>#include<math.h>#define ee 0.006694384999588#define I 3.141592653double F(double,double,double);void main(void){double A1,B1,L1,S,A2,B2,L2; double x1,x2,x3,y1,y2,y3,z1,z2,z3; double W1,sinu1,sinu2,cosu1,sinA0;doublecota1,cos2a1,sin2a1,cosA0A0;double A,B,C,d,e,a0,a1,m;double n,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=6356755.288+(10710.341-(13.534*cosA0A0))*cosA0A0;B=(5355.171-9.023*cosA0A0)*cosA0A0;C=(2.256*(cosA0A0))*cosA0A0+0.006;d=691.46768-(0.58143-0.00144*cosA0A0)*cosA0A0;e=(0.2907-cosA0A0*0.0010)*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(2*a0)+cos2a1*sin(2*a0);n=(cos2a1)*(cos(2*a0))-(sin2a1)*(sin(2*a0));a=a0+((B+5*C*n))*m/A;printf("a0=%f\n",a0);printf("m=%f\n",m);printf("n=%f\n",n);printf("a=%f\n",a);/*计算经度差改正数*/Q=(d*a+(e*(m-sin2a1))/3600/180*I)*sinA0;printf("Q=%f\n",Q);/*计算终点大地坐标及大地方位角*/sinu2=sinu1*cos(a)+cosu1*cos(A1)*sin(a);B2=180*atan(sinu2/((sqrt(1-ee))*(sq rt(1-sinu2*sinu2))))/I;R=180*atan(sin(A1)*sin(a)/(cosu1*co s(a)-sinu1*sin(a)*cos(A1)))/I;printf("sinu2=%f\n",sinu2);printf("B2=%f\n",B2);printf("R=%f\n",R*180/I);/*确定R的值*/if(sin(A1)>0 && tan(R)>0)R=abs(R);else if(sin(A1)>0 && tan(R)<0)R=I-abs(R);else if(sin(A1)<0 && tan(R)<0)R=-abs(R);elseR=abs(R)-I;/*确定L2A2的值*/L2=(L1*180/I+R-(Q/206265*180/I));A2=atan(cosu1*sin(A1)/(cosu1*cos(a) *cos(A1)-sinu1*sin(a)));if(sin(A1)<0&&tan(A2)>0)A2=(fabs(A2))*180/I;else if(sin(A1)<0&&tan(A2)<0)A2=(I-fabs(A2))*180/I;else if(sin(A1)>0&&tan(A2)>0)A2=(I+fabs(A2))*180/I;elseA2=(2*I-fabs(A2))*180/I;printf("A2=%3f\n B2=%3f\nL2=%3f\n",A2,B2,L2);}double F(double a2,double b2,doublec2){double d2;d2=(double)(a2+1.0*b2/60+1.0*c2/3600);d2=(d2/180)*I;return (d2);}注：A1，B1，L1，S分别为大地线起点的大地方位角、纬度、经度、大地线长;B2,L2,A2为大地线终点纬度、经度及方位角。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
(4 192) ( 4 194)
三阶导数
d 3B V5 2 2 2 2 2 2 2 2 2 2 2 cos A [sin A ( 1 3 t 9 t ) 3 cos A ( 1 t 5 η t )] 3 3 dS c d 3 L 2V 2 2 t sec B sin A cos A 3 dS c
11
Fundation of Geodesy
(1)建立级数展开式:
S MP , 2 2 S MP 1 2
dB S 1 d 2 B S 2 1 d 3 B S 3 B2 BM ( ) M ( 2 ) ( 3 ) dS 2 2 dS 4 6 dS 8
dB S 1 d 2 B S 2 1 d 3 B S 3 B1 BM ( ) M ( 2 ) M ( 3 )M dS 2 2 dS 4 6 dS 8
(
dB dB ) ( ) dB ( ) M f ( Bm , Am ) ( dS )( BM Bm ) ( dS )( AM Am ) + dS Bm Am
S 2 d 2B BM Bm ( 2 )M 8 dS S 2 d 2B ( 2 )m 8 dS
二阶导数：
d 2 B dB dB dB dA V4 2 2 2 ( ) ( ) t ( 3 cos A sin A) 2 2 dS B dS dS A dS dS c ( 4 190)
6
Fundation of Geodesy
d 2 L dL dB dL dA 2V 2 ( ) ( ) 2 t sec B sin A cos A 2 dS B dS dS A dS dS c d 2 A dA dB dA dA V 2 2 2 ( ) ( ) sin A cos A ( 1 2 t ) 2 2 dS B dS dS A dS dS c
cos A dS M sin A dS N cos B tan B sin AdS N
Fundation of Geodesy
2.以白塞尔大地投影为基础
1) 按椭球面上的已知值计算球面相应值，即实现椭球面 2) 3) 按球面上得到的数值计算椭球面上的相应数值，即实 现从圆球向椭球的过渡。 典型解法：白塞尔大地主题解算 特点：解算精度与距离长短无关，它既适用于短距离 解算，也适用于长距离解算。可适应20 000km 或更长的 距离，这对于国际联测，精密导航，远程导弹发射等都 具有重要意义。
3
Fundation of Geodesy
4L2 L( S ), A21 A( S ) B(0) B1 , L(0) L1 , A(0) A12
d nB S n dB d 2B S 2 d 3B S 3 B2 B1 B ( n )1 ( )1 S ( 2 )1 ( 3 )1 dS n! dS dS 2! dS 3!
4 2 2 5m ) sin 2 Am ( 2 t m 2m )]} 5次
"
18
Fundation of Geodesy
B2 B1 B,
L2 L1 L,
A21 A12 A 180
B2 B1 1 1 Bm ( B2 B1 ) B1 B1 B 2 2 2
d 3 L 2V 3 2 2 2 2 2 sec B [sin A cos A ( 1 + 3 t ) t sin A] 3 3 dS c
7 Fundation of Geodesy
u S cos A1
v S sin A1
8
Fundation of Geodesy
9
Fundation of Geodesy
1 Am A12 A 2
注意： 从公式可知，欲求Δ Ｌ，Δ Ｂ及Δ Ａ，必先有Ｂｍ及Ａ ｍ。但由于Ｂ2和Ａ21未知，故精确值尚不知，为此须 用逐次趋近的迭代方法进行公式的计算。 除此之外，此方法适合与200公里以下的大地问题解算， 其计算经纬计算精度可达到0.0001”, 方位角计算精度 可达到0.001”。
主题解算分为: 短距离(<400km) 中距离(<1000km) 长距离(1000km以上)
1
Fundation of Geodesy
1.以大地线在大地坐标系中的微分方程为基础，直 接在地球椭球面上进行积分运算。
• 主要特点：解算精度与距离有关，距离越长， 收敛越慢，因此只适用于较短的距离 典型解法：高斯平均引数法
( 4 190)
d 2 B dB dB dB dA V4 2 2 2 ( ) ( ) t ( 3 cos A sin A) dS 2 B dS dS A dS dS c2 d 2 L dL dB dL dA 2V 2 ( ) ( ) 2 t sec B sin A cos A 2 dS B dS dS A dS dS c d 2 A dA dB dA dA V 2 ( ) ( ) 2 sin A cos A(1 2t 2 2 ) 2 dS B dS dS A dS dS c
19 Fundation of Geodesy
4.7.4 高斯平均引数反算公式 高斯平均引数反算公式可以依正算公式导出:
S sin Am 2 2 L " 2 S sin Am N m cos Bm [ S t sin Am m 2 " 24 N m
2 2 2 S 2 cos 2 Am (1 m 9m tm )]
为了计算 B, L, A 的级数展开式，关键问题是推求 各阶导数。
4 Fundation of Geodesy
dB dS dL dS dA dS
cos A V 3 cos A M c sin A V sec B sin B N cos B c tan B V sin A tan B sin A N c
5 Fundation of Geodesy
(4 192)
( 4 194)
一阶导数：
dB cos A V 3 dS M c cos A sin A V dL sec B sin B N cos B c dS V dA tan B dS N sin A c tan B sin A
13
Am AM
Fundation of Geodesy
dB ( ) M f ( BM , AM ) F ( Bm BM Bm , Am AM Am ) + dS
( dB f f ) M f ( Bm , Am ) ( )( BM Bm ) ( )( AM Am ) + dS Bm Am
cos A dB dS M sin A dL N cos B dS tan B dA sin A dS N
2
P 2 B2 B1 P 1 P 2 L2 L1 P 1 P 2 A A 1 2 P 1
dA d 3 A A21 A12 A ( ) M S ( 3 ) M S 3 dS 24 dS
(2) BM , AM Bm , Am
Bm 1 1 ( B2 B1 ), Am ( A21 A12 180 ) 2 2
Bm BM ,
14
S2 d2A AM Am ( ) 2 M 8 dS S2 d2A ( ) 2 m 8 dS
Fundation of Geodesy
(3)由大地线微分方程依次求偏导数:
3 2 cos Am Vm Vm dB ( )m cos Am cos Am dS Mm c Nm
dnL Sn dL d 2L S2 d 3L S3 L2 L1 L ( n )1 ( )1 S ( 2 )1 ( 3 )1 dS n! dS dS 2! dS 3!
d n A S n dA d 2 A S2 d 3A S3 A2 A1 180 A ( n )1 ( )1 S ( 2 )1 ( 3 )1 dS n ! dS dS 2 ! dS 3!
dB d 3 B B2 B1 B ( ) M S ( 3 )M S 3 dS 24 dS
12 Fundation of Geodesy
(4 200)
(4 201)
(4 202)
同理可得:
dL d 3 L L2 L1 L ( ) M S ( 3 )M S 3 dS 24 dS
3 Vm dB ( ) ( cos Am ) 3 2 c ( dS )m t mm cos Am B B Nm
3 Vm dB ( ) ( cos Am ) 2 V c ( dS )m m sin Am A A Nm
2 Vm S 3 d 3B 2 2 2 2 2 ( 3 )M cos A [sin A ( 1 3 t 9 t m m m m m m ) + 3 24 dS 24 N m 2 2 2 2 2 3m cos 2 Am (1 t m +m 5t m m ) S 3 + 5次
S cos Am
B " N m S cos Am 2 2 2 2 [ S sin A (2 3 t 2 ) m m m 2 2 " Vm 24 N m
2 2 2 2 2 2 3m S cos 2 Am (t m 1 m 4m t m )]
4.7 大地测量主题解算概述
4.7.1 大地主题解算的一般说明
正算 : 已知 P 1( B 1, L 1 ), S12 , A 12 求 P 2 ( B2 , L2 ), A 21 反算 : 已知 P ，求 S12 , A12 , A21 1( B 1, L 1 ), P 2 ( B2 , L2 )