《基于边角联合后方交会原理的自由设站程序 for 4850p》由会员分享，可在线阅读，更多相关《基于边角联合后方交会原理的自由设站程序 for 4850p（4页珍藏版）》请在金锄头文库上搜索。
1、基于边角联合后方交会原理的自由设站程序 for 4850P一、程序编制原理：基于边角联合后方交会法，按最小二乘间接条件严密平差原理二、源程序：1、主程序：JIAOHUI-BJ1Lbl 1:XYABCDRGEFQ:XXA:YYA:AXB:BYB:CD1:DD2:Ra0:GMR:EMS1:FMS2:QK(1,-1)K=F:F=R:ProgDMS2: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:ProgCOS:W=O:K=C:C=D:D=K:Z=3600(180-W-V-R):H=CsinW:P=206265:K=H2G
2、2Z(H2G2+P2(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)1E3MP=:M=(KZ)(E2+F2)(sinR)2+P2E2F2H2G2):Int(1E7M+0.5)1E4Goto 1
3、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 测距中误差（
4、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.1
5、06m,1444340，现求设站 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=1444332.632、CP=2221.9863、DP=1053.1
6、044、XP=60980.6695、YP=97212.2926、MP=14.577五、参考文献：1、测绘工程2000 年 9 月（第 3 期）李全信边角后方交会的精度分析及布设方案选择2、城市道桥与防洪2005 年第 1 期宋润嘉等几种典型自由设站法测量的解算及精度分析沙发：hanniu_1973 楼主 09-10-12 15:58 回基于这样精良的数学建模和详细算例的自由设站程序，居然没人顶？2 楼：我很认可 09-10-12 18:41 回评：+10 分楼主，你好，这方面的程序实在是太少，所以很高兴能有你提供共享，我是新手，有几个问题想请教你：一、Ra0中的 a0 是如何输入？二、Int(
7、1E3(D+M)+0.5)1E3 中的 1E3 是不是 10 的三次方还是这个 E 是代表 EMS1中的 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 中如何更改输入？我很想用你的程序，
8、可我用的是 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-
9、3.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)1E
10、3改为:CCP=Int(1E3(C+L)+0.5)1E3:DDP=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改为:MMP=(KZ)(E2+F2)(sinR)2+P2E2F2H2G2):Int(1E7M+0.5)1E4请楼主务必帮我看一下,实在不好意思,麻烦你了.谢谢!我会努力的.7 楼：小测量 09-10-14 12:4 回评：+10 分我只有 4800 计算器,希望楼主,能具体帮助改一下,我实在没辙了.谢谢!
