自由设站程序5800485097509860cg20.doc

基于边角联合后方交会原理的自由设站程序 后方交会for 5800（大歪哥2D后方交会改编本网站原4850版之程序）程序编制原理：基于边角联合后方交会法，按最小二乘间接条件严密平差原理参考文献：1、《测绘工程》2000年9月（第3期）——李全信《边角后方交会的精度分析及布设方案选择》2、《城市道桥与防洪》2005年第1期——宋润嘉等《几种典型自由设站法测量的解算及精度分析》源程序代码：程序名：HJLbl 1ClrXA_？XYA？YXB？AYB？BHD1？CHD2？D_P?RCeJiao Zhong Wu Cha?GHD-1 Zhong Wu Cha?EHD-2 Zhong Wu Cha?F+-1?QF-K:R-FF-R:K-FPol(X-A,Y-B360Frac ((J+360)/360-Jcos-1 ((C^2+I^2-D^2)÷(2CI))-OO-V:C-K:D-C:K-D:cos^-1 ((C^2+I^2-D^2)÷(2CI))-OO-W:C-K:D-C:K-D3600(180-W-V-R)-ZCsin W-H206265-PH^2G^2Z/(H^2G^2+P^2((cos W)^2F^2))-KKP/G^2/H-NNE^2cos V-LNF^2cos W-MR+K/3600-RInt (Exp3(C+L)+0.5)/Exp3-CInt (Exp3(D+M)+0.5)/Exp3-D(C^2+I^2-D^2)/2/I^2-LSqrt(C^2/I^2-L^2)-MA-X-VB-Y-WX+VL-QWM-UY+WL+QVM-SSqrt((K/Z)((E^2+F^2)/(sin R)^2+P^2E^2F^2/H^2/G^2))-MClrText_＜P=:Locate 5,1,R:R黑三角DMS黑三角PA=:PB=:XP=:YP=:MP=:Locate 5,1,C:Locate 5,2,D:Locate 5,3,Int (Exp3U+0.5)/Exp3:Locate 5,4,Int (Exp3S+0.5)/Exp3:Locate 5,5,M黑三角Goto 1三、说明：1、XA是已知点A的x坐标；2、YA是已知点A的y坐标；3、XB是已知点B的x坐标；4、YB是已知点B的y坐标；5、HD1是自由设站点P至已知点A的测距；6、HD2是自由设站点P至已知点B的测距；7、P是PA与PB的夹角，即∠APB；8、CeJiao Zhong Wu Cha是测角中误差（秒），如：MR=5；9、HD-1 Zhong Wu Cha是PA测距中误差（mm);10、HD-2 Zhong Wu Cha是PB测距中误差（mm);11、+-是三角形PAB的编号顺序，顺时针编号输入1，逆时针编号输入-1;12、＜P=是经过平差后的夹角∠APB；13、PA=是经过平差后的PA测距；14、PB=是经过平差后的PB测距；15、XP=是经过平差后的自由设站点P的x坐标；16、YP=是经过平差后的自由设站点P的y坐标；17、MP=是经过平差后的点位中误差(m);算例：已知点A（61749.436,95127.533)，B（60112.617,97808.541);在P点设置GTS-310D型全站仪测得PA、PB及夹角a0值分别为2221.99m,1053.106m,144°43′40″，现求设站P点的坐标及点位精度。运行主程序JIAOHUI-BJ1：输入数据：1、XA？61749.436????回车2、YA？95127.533????回车3、XB？60112.617????回车4、YB？97808.541????回车5、HD1？2221.99????回车6、HD2？1053.106????回车7、P？144.4340????回车8、CeJiao Zhong Wu Cha？5???????? 回车9、HD-1 Zhong Wu Cha？7.4E-3????回车10、HD-2 Zhong Wu Cha？5.1*E-3????回车11、+ -1？-1????回车(注：三角形PAB按逆时针编号输入，则输入-1）显示：1、P=144.726144°43′31.86″2、PA=2221.9863、PB=1053.1044、XP=60980.6695、YP=97212.2926、MP=0.0015 测边交会Lbl1:“X1”？C:“Y1”?O：“X2”？X：“Y2”?Y：“D1”?A：“D2”?B：0→I：1→J:Pol（X-C，Y-?O?:J0＝J