白塞尔大地主题解算.docxVIP

白塞尔大地主题解算 方向: 学号: 姓名: 一.基本思路： 基本思想：将椭球而上的大地元素按照白塞尔投影条件投影到辅助球而上，继而在球而 上进行大地主题解算，最后在将球而上的汁算结果换算到椭球而上。其关键问题是找出椭球 而上的大地元素与球而上相应元素之间的关系式，同时解决在球而上进行大地主题解算的方 法。 正算流程： 计算起点的归化纬度 计算辅助函数值，解球而三角形 按公式计算相关系数A,B,C以及(】，B 计算球面长度 巧=丄{3 -sin2cr1(B + Ccos2a1)} A sin 2(er】+ = sin 2? sin 2q)+ cos 2cr} cos 2q) cos 2( b] + 6) = cos 2b] cos 一 sin 2r: sin 2t{) b = 6 + 丄[B + 5C cos 2( + 巧)]sin 2( + cr0) A 5 ?计算纬度差改正数 A-L = 3 = sinAJab + 0(shi2(bi +tr(J—sin2bj] 6.汁算终点大地坐标及大地方位角 sin u2 = sin u} cos b + cos 绚 cos sin a \j\-e2 sin = sin B. W2 cos 心=——cosB? ~ vv2 ? tanw2 = Jl一，tan优 ,2 n 1 siir 心 ? tair = = arctan 1-e^ 1 一 cos* 冲 r sino-sinA , ?I = arctan[ J cosz/j costr-sinz/, sinacosA, A. =arctan cos u{ sin A】 cos u} cos er cos Ax - sin ux sin a 1 J]./ sin I/? Jl-cos w2 反算流程： 1.辅助计算 2 ?用逐次趋近法同时讣算起点大地方位角、球而长度及经差，第一次趋近时，取8 = 0. tanAj = sin 2 cos 心 cos “i sin u2 - sin u} cos u2 cos A p = sin 2 cos m. q = b\-b、cos A A} = arctan — q ps\nA} +qcosA tan b = ￡ ? COSb cosflr = t/l +a2 cos A 讣算下式，重复上述汁算过程2. sin cr = psinA{ +qcosA} (sincr) (T = arctan {COSffJ sin Ao = cos Wj sin A} tan b】=tan ? sec k-L = 8 -sin4)[eb + 0(sin2b2 -sin2b【)] 3?计算大地线长度S 刁=s+b A = L+S S = 4b + sin2b](〃 + Ccos2b])-shi26(B + Ccos26) 4 ?讣算反方位角 cosw. sin 2 A. = arctan ! S cosA-b, 二.已知数据 序号 Bl ( LI ? A12 S12(m) 1 130. 8000 三.源代码: #include #include #define e 〃克拉索夫斯基椭球体第一偏心率 void main() { int kBlO’Bll丄 10丄X,A10,AlbB20,B2bL20丄2bA20,A21; double B12丄 1 乙A12，B22,L22,A22; double Bi丄丄2，A2丄，pi; ( double A,B,C,afa，beta; double al,a2,bl,b2,p,c\,x,y; double Wl/W2,si nul，si nu2,cosubcosu2,si nA0,cotsigmabsin2sigmahcos2sigma：Usigma0，sin2,cos2,sigma ，sins，coss，deltaO,delta」amda; pi=4*atan(l); printfC1白塞尔大地主题正算请输入l\n白塞尔大地主题反算请输入2\n)； scanf(吆 cT,k); if(k==l) { printf(H请输入大地线起点纬度B经度L,大地方位角A及大地线长度S:\nH)； scanf(,,%d%d%lf%d%d%lf%d%d%lf%lf,,,B10/BlbB12,L20,LlhL：L2,A10,Al：LAl 2,S); Bl=(B10+(float)B11^0+B1^3600)*pi/180; Ll=(L10+(float)L11^0+L1^3600)*pi/180; Al=(A10+(float)All/60+A1^3600)*pi/180; 〃讣算起点规划纬度〃计算起点规划纬度 〃讣算起点规划纬度 〃计算起点规划纬度 〃计算起点规划纬度 〃计算辅