1、基于边角联合后方交会原理的自由设站程序 后方交会 for 5800（大歪哥 2D后方交会改编本网站原 4850版之程序）程序编制原理：基于边角联合后方交会法，按最小二乘间接条件严密平差原理参考文献：1、测绘工程2000 年 9月（第 3期）李全信边角后方交会的精度分析及布设方案选择2、城市道桥与防洪2005 年第 1期宋润嘉等几种典型自由设站法测量的解算及精度分析源程序代码：程序名：HJLbl 1Clr“XA_“？X“YA“？Y“XB“？A“YB“？B“HD1“？C“HD2“？D“_K:R-FF-R:K-FPol(X-A,Y-B360Frac (J+360)/360-Jcos-1 (C+I-D
2、)(2CI)-OO-V:C-K:D-C:K-D:cos-1 (C+I-D)(2CI)-OO-W:C-K:D-C:K-D3600(180-W-V-R)-ZCsin W-H206265-PHGZ/(HG+P(cos W)F)-KKP/G/H-NNEcos V-LNFcos W-MR+K/3600-RInt (Exp3(C+L)+0.5)/Exp3-CInt (Exp3(D+M)+0.5)/Exp3-D(C+I-D)/2/I-LSqrt(C/I-L)-MA-X-VB-Y-WX+VL-QWM-UY+WL+QVM-SSqrt(K/Z)(E+F)/(sin R)+PEF/H/G)-MClrText“_P=
3、“: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、J+360J:co
4、s-1（I 的平方+B 的平方-A 的平方）（2IB）K:J-KW:WW+360W:“FWJ（OP2）=”:WDMS “X”:X-Bcos(W)M“Y”:Y-Bsin(W)NGoto 1程序说明：X1 表示目标 1纵坐标，Y1 表示目标 1的横坐标；X2 表示目标 2的纵坐标，Y2 表示目标 2横坐标；D1 表示置仪点 O 到目标 1的距离D2 表示置仪点 O 到目标 2的距离显示 FWJ=？（表示 O 到目标 2的方位角）显示置仪点 O 的纵横坐标：X=？；Y=？基于边角联合后方交会原理的自由设站程序 for 4850P基于边角联合后方交会原理的自由设站程序 for 4850P一、程序编制原
5、理：基于边角联合后方交会法，按最小二乘间接条件严密平差原理二、源程序：1、主程序：JIAOHUI-BJ1Lbl 1:XYABCDRGEFQ:X“XA“:Y“YA“:A“XB“:B“YB“:C“D1“:D“D2“:R“a0“:G“MR“:E“MS1“:F“MS2“:Q“K(1,-1)“K=F:F=R:Prog“DMS2“:R=F:F=K:J=0:I=Pol(X-A,Y-B):J0=J=J+360Prog“COS”:V=O:K=C:C=D:D=K:Prog“COS“:W=O:K=C:C=D:D=K:Z=3600(180-W-V-R):H=CsinW:P=206265:K=H2G2Z(H2G2+P2
6、(cosW)2F2):N=KPG2H:L=NE2cosV:M=NF2cosW:R=R+K3600:“a0=“:RDMS“CP=”：C=Int(1E3(C+L)+0.5)1E3:Pause 0:“DP=“:D=Int(1E3(D+M)+0.5)1E3L=(C2+I2-D2)2I2:M=(C2I2-L2):V=A-X:W=B-Y:U=X+VL-QWM:S=Y+WL+QVM:“XP=“:Int(1E3U+0.5)1E3:Pause 0:“YP=“:Int(1E3S+0.5)1E3“MP=“:M=(KZ)(E2+F2)(sinR)2+P2E2F2H2G2):Int(1E7M+0.5)1E4Goto 1
7、2、子程序：COSO=cos-1(C2+I2-D2)(2CI)3、子程序：DMS2F=Int F+Frac 100F36+(Int 100F100-Int F)0.6三、说明：1、XA 是已知点 A的 x坐标；2、YA 是已知点 A的 y坐标；3、XB 是已知点 B的 x坐标；4、YB 是已知点 B的 y坐标；5、D1 是自由设站点 P至已知点 A的测距；6、D2 是自由设站点 P至已知点 B的测距；7、a0 是 PA与 PB的夹角，即APB；8、MR 是测角中误差（秒），如：MR=5；9、MS1 是 PA测距中误差（mm);10、MS2 是 PB测距中误差（mm);11、K(1,-1)是三角
8、形 PAB的编号顺序，逆时针编号输入-1，逆时针编号输入 1;12、a0=是经过平差后的夹角APB；13、CP=是经过平差后的 PA测距；14、DP=是经过平差后的 PB测距；15、XP=是经过平差后的自由设站点 P的 x坐标；16、YP=是经过平差后的自由设站点 P的 y坐标；17、MP=是经过平差后的点位中误差(mm);四、算例：已知点 A（61749.436,95127.533)，B（60112.617,97808.541);在 P点设置 GTS-310D型全站仪测得 PA、PB 及夹角 a0值分别为 2221.99m,1053.106m,1444340，现求设站 P点的坐标及点位精度。
9、运行主程序 JIAOHUI-BJ1：输入数据：1、XA？61749.436 回车2、YA？95127.533 回车3、XB？60112.617 回车4、YB？97808.541 回车5、D1？2221.99 回车6、D2？1053.106 回车7、a0？144.4340 回车8、MR？5 回车9、MS1？7.4*1E-3 回车10、MS2？5.1*1E-3 回车11、K（-1，1）？-1 回车(注：三角形 PAB按逆时针编号输入，则输入-1）显示：1、a0=1444332.632、CP=2221.9863、DP=1053.1044、XP=60980.6695、YP=97212.2926、MP=
10、14.577五、参考文献：1、测绘工程2000 年 9月（第 3期）李全信边角后方交会的精度分析及布设方案选择2、城市道桥与防洪2005 年第 1期宋润嘉等几种典型自由设站法测量的解算及精度分析谢谢发贴：+16 分 【测量空间】0查看楼主其它文章 返回 回复内容 只查看楼主与本人回复 倒序沙发：hanniu_1973 楼主 09-10-12 15:58基于这样精良的数学建模和详细算例的自由设站程序，居然没人顶？2 楼：我很认可 09-10-12 18:41评：+10 分楼主，你好，这方面的程序实在是太少，所以很高兴能有你提供共享，我是新手，有几个问题想请教你：一、R“a0“中的 a0是如何输入
11、？二、Int(1E3(D+M)+0.5)1E3 中的 1E3是不是 10的三次方还是这个 E是代表 E“MS1“中的 E？三、就是 4800与 4850计算器中的输出部分更改问题：（1）、RDMS（2）、“CP=”：C=Int(1E3(C+L)+0.5)1E3:Pause 0:“DP=“:D=Int(1E3(D+M)+0.5)1E3（3）、XP=“:Int(1E3U+0.5)1E3:Pause 0:“YP=“:Int(1E3S+0.5)1E3（4）、“MP=“:M=(KZ)(E2+F2)(sinR)2+P2E2F2H2G2):Int(1E7M+0.5)1E4以上四句在 4800中如何更改输入？
12、我很想用你的程序，可我用的是 4800计算器，苦于是新手，所以只能在这里求救于您了，望能得到楼主的解答，最后再次感谢你的好程序！3 楼：我也来顶 09-10-12 18:44评：+10 分楼主，我也是新手，对了，你的程序中的 K=F:F=r中的 r在计算器中是如何输入的，小弟不太懂，希望得到你的指导，谢谢！还有测角与测距的中误差都得输入？这些中误差是如何来确定的？能给个说明吗？4 楼：测量 09-10-12 19:21看测量学吧，学习的过程中了解公式的推导过程很重要5 楼：hanniu_1973 楼主 09-10-13 15:0回 2楼的朋友：1、a0 在按 functon-4.const-3
13、.a02、1E3 是 10的 3次方；3、改成 C“CP=”。M“MP=”请这位朋友多多看说明书，这些问题都比较菜，多努力哦！回 3楼的朋友：1、其中的“r“是输入有误，应为”R“。2、测角与测距中误差请看仪器说明书或测量书注：这个数学建模是比较新的，在现有的测量教科书上是找不到的，请查阅我上面所提的论文。谢谢各位的支持！6 楼：我很认可 09-10-13 20:42评：+10 分楼主,在 4800里,下面几句是不是这样改?（1）、RDMS改为 R（2）、“CP=”：C=Int(1E3(C+L)+0.5)1E3:Pause 0:“DP=“:D=Int(1E3(D+M)+0.5)1E3改为:C“
14、CP“=Int(1E3(C+L)+0.5)1E3:D“DP“=Int(1E3(D+M)+0.5)1E3（3）“XP=“:Int(1E3U+0.5)1E3:Pause 0:“YP=“:Int(1E3S+0.5)1E3改为:“XP“=Int(1E3U+0.5)1E3:“YP“=Int(1E3S+0.5)1E3（4）、“MP=“:M=(KZ)(E2+F2)(sinR)2+P2E2F2H2G2):Int(1E7M+0.5)1E4改为:M“MP“=(KZ)(E2+F2)(sinR)2+P2E2F2H2G2):Int(1E7M+0.5)1E4请楼主务必帮我看一下,实在不好意思,麻烦你了.谢谢!我会努力的.
15、7 楼：小测量 09-10-14 12:4评：+10 分我只有 4800计算器,希望楼主,能具体帮助改一下,我实在没辙了.谢谢!8 楼：紫云 09-10-15 17:53 评：+8 分我输了 程序但结果和老大的差几个毫米 不知道为什么9 楼：water3000 09-10-15 21:58那个全站仪都有后方交会程序吧，还用计算器后方交会？10 楼：hanniu_1973 楼主 09-10-16 16:229楼的朋友，请看清这是基于两边夹一角的联合后方交会法，按最小二乘间接条件严密平差原理编制的程序，现场自由设站时仅需两个已知点即可，而大多全站仪则需三个以上已知点，这在现场条件有限的情况下，此程序更显得可取之处。Good Luck!
