基于边角联合后方交会原理的自由设站程序 for 4850P.doc

基于边角联合后方交会原理的自由设站程序 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=H2G2Z(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 12、子程序：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)是三角形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点的坐标及点位精度。运行主程序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=14.577五、参考文献：1、测绘工程2000年9月（第3期）李全信边角后方交会的精度分析及布设方案选择2、城市道桥与防洪2005年第1期宋润嘉等几种典型自由设站法测量的解算及精度分析1）、RDMS改为R（2）、“CP=”：C=Int(1E3(C+L)+0.5)1E3:Pause 0:DP=:D=Int(1E3(D+M)+0.5)1E3改为: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