自由设站程序5800485097509860cg20.doc

基于边角联合后方交会原理的自由设站程序 后方交会for 5800
（大歪哥2D后方交会改编本网站原4850版之程序）
程序编制原理：基于边角联合后方交会法，按最小二乘间接条件严密平差原理
参考文献：
1、《测绘工程》2000年9月（第3期）——李全信《边角后方交会的精度分析及布设方案选择》
2、《城市道桥与防洪》2005年第1期——宋润嘉等《几种典型自由设站法测量的解算及精度分析》

源程序代码：
程序名：HJ
Lbl 1
Clr
"XA_"？X
"YA"？Y
"XB"？A
"YB"？B
"HD1"？C
"HD2"？D
"_ng Wu Cha"?G
"HD-1 Zhong Wu Cha"?E
"HD-2 Zhong Wu Cha"?F
"+-1"?Q
F->K:R->F
F->R:K->F
Pol(X-A,Y-B
360Frac ((J+360)/360->J
cos-1 ((C^<2>+I^<2>-D^<2>)÷(2CI))->O
O->V:C->K:D->C:K->D:cos^-1 ((C^<2>+I^<2>-D^<2>)÷(2CI))->O
O->W:C->K:D->C:K->D
3600(180-W-V-R)->Z
Csin W->H
206265->P
H^<2>G^<2>Z/(H^<2>G^<2>+P^<2>((cos W)^<2>F^<2>))->K
KP/G^<2>/H->N
NE^<2>cos V->L
NF^<2>cos W->M
R+K/3600->R
Int (Exp3(C+L)+0.5)/Exp3->C
Int (Exp3(D+M)+0.5)/Exp3->D
(C^<2>+I^<2>-D^<2>)/2/I^<2>->L
Sqrt(C^<2>/I^<2>-L^<2>)->M
A-X->V
B-Y->W
X+VL-QWM->U
Y+WL+QVM->S
Sqrt((K/Z)((E^<2>+F^<2>)/(sin R)^<2>+P^<2>E^<2>F^<2>/H^<2>/G^<2>))->M
ClrText
"_＜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的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、80.669
5、YP=97212.292
6、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 :J<0＝>J+360→J:cos-1（（I的平方+B的平方-A的平方）÷（2IB））→K:J-K→W:
W<0=>W+360→W:“FWJ（OP2）=”:W◆DMS ◢
“X＝”:X-Bcos(W)→M◢
“Y＝”:Y-Bsin(W)→N◢
Goto  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

一、程序编制原理：基于边角联合后方交会法，按最小二乘间接条件严密平差原理
二、源程序：
1、主程序：JIAOHUI-BJ1
Lbl 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):J≤0=>J=J+360△Prog“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=H^2G^2Z÷(H^2G^2+
P^2((cosW)^2F^2)):N=KP÷G^2÷H:L=NE^2cosV:M=NF^2cosW:R=R+K÷3600:"a0=":R→DMS▲
“CP=”：C=Int(1E3(C+L)+0.5)÷1E3:Pause 0:"DP=":D=Int(1E3(D+M)+0.5)÷1E3▲
L=(C^2+I^2-D^2)÷2÷I^2:M=√(C^2÷I^2-L^2):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=√((K÷Z)((E^2+F^2)÷(sinR)^2+P^2E^2F^2÷H^2÷G^2)):Int(1E7M+0.5)÷1E4▲Goto 1
2、子程序：COS
O=cos-1‎((C^2+I^2-D^2)÷(2CI))
3、子程序：DMS2
F=Int F+Frac 100F÷36+(Int 100F÷100-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)是三角形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,144°
43′40″，现求设站P点的坐标及点位精度。
运行主程序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=144°43′32.63″
2、CP=2221.986
3、DP=1053.104
4、XP=60980.669
5、YP=97212.292
6、MP=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是如何输入？
二、Int(1E3(D+M)+0.5)÷1E3中的1E3是不是10的三次方还是这个E是代表E"MS1"中的E？
三、就是4800与4850计算器中的输出部分更改问题：
（1）、R→DMS▲
（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=√((K÷Z)((E^2+F^2)÷(sinR)^2+P^2E^2F^2÷H^2÷G^2)):Int(1E7M+0.5)÷1E4▲
以上四句在4800中如何更改输入？
我很想用你的程序，可我用的是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.a0
2、1E3是10的3次方；
3、改成C“CP=”。。。。M“MP=”
请这位朋友多多看说明书，这些问题都比较菜，多努力哦！

回3楼的朋友：1、其中的“r“是输入有误，应为”R“。
2、测角与测距中误差请看仪器说明书或测量书
注：这个数学建模是比较新的，在现有的测量教科书上是找不到的，请查阅我上面所提的论文。
谢谢各位的支持！
6 楼：'我很认可
09-10-13 20:42
评：+10 分

楼主,在4800里,下面几句是不是这样改?
（1）、R→DMS▲改为R▲
（2）、“CP=”：C=Int(1E3(C+L)+0.5)÷1E3:Pause 0:"DP=":D=Int(1E3(D+M)+0.5)÷1E3▲
改为:C"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=√((K÷Z)((E^2+F^2)÷(sinR)^2+P^2E^2F^2÷H^2÷G^2)):Int(1E7M+0.5)÷1E4▲
改为:M"MP"=√((K÷Z)((E^2+F^2)÷(sinR)^2+P^2E^2F^2÷H^2÷G^2)):Int(1E7M+0.5)÷1E4▲
请楼主务必帮我看一下,实在不好意思,麻烦你了.谢谢!!!!!!!!!我会努力的.
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:22

9楼的朋友，请看清这是基于两边夹一角的联合后方交会法，按最小二乘间接条件严密平差原理编制的程序，现场自由设站时仅需两个已知点即可，而大多全站仪则需三个以上已知点，这在现场条件有限的情况下，此程序更显得可取之处。Good Luck!