附和导线近似平差计算程序附和导线近似平差计算程序P:A:B:I:J:K:L:Defm12 P=P+2:N=0:M=A W=60√P÷3600 LbI0 N=N+1 {C} Z[2N-1]=C:M=M+C M>180=>M=M-180:≠=>M=M+180⊿ NGoto0⊿ F=M-B:F>360=>F=F-360⊿F▲ AbsFF=-F÷P:≠=>Goto E⊿ N=0:M=A LbI1 N=N+1:M=M+Z[2N-1]+F M>180=>M=M-180▲ ≠=>M=M+180▲⊿ Z[2N-1]=M ...
Goto0⊿
F=M-B:F>360=>F=F-360⊿F▲
AbsF Goto1⊿
N=0:M=0:G=0:H=0:P=P-1
LbI2
N=N+1
{D}:Z[2N]=D:M=M+D
X=DcosZ[2N-1]:Y=DsinZ[2N-1]
G=G+X:H=H+Y
Z[2N-1]=X:Z[2N]=Y
N Goto2⊿
G=G+I-K▲H=H+J-L▲
T=√(G×G+H×H)▲
Q=M÷T▲
Q>2000=>G=-G÷M:H=-H÷M:≠=>Goto E⊿
N=0:X=I:Y=J
LbI3
N=N+1▲
D=√(Z[2N-1]×[2N-1]+Z[2N]×[2N])
V=Z[2N-1]+DG▲
U=Z[2N]+DH▲
X=X+V▲Y=Y+U▲
Z[2N-1]=X:Z[2N]=Y
N Goto3⊿
LbI E
说明:1程序中Defm12是按4个未知点设置的,计算时要进行修改。如未知点数为n,修改数为2(n+2)。
1:P输入未知点数
2:A输入起点方位角
3:输入终点方位角
4:I,J输入起点坐标
5:K,L输入终点坐标
6:C从导线起点依次观测角
7:F角度闭合差
8:M导线第一边依次方位角平差值
9:D导线第一边依次输边长
10:G,H为(⊿X,⊿Y)导线闭合差11:T为导线绝对误差
12:Q为导线相对闭合差分母,分子为1 13:N导线未知点数
14:V,U为导线第一边坐标增量平差值15:X,Y为导线未知点坐标平差值