附和导线近似平差计算程序.doc

附和导线近似平差计算程序
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⊿
N＜P=>Goto0⊿
F=M-B：F＞360=>F=F-360⊿F▲
AbsF＜W=>F=-F÷P：≠=>GotoE⊿
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
N＜P=>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＜P=>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：≠=>GotoE⊿
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＜P=>Goto3⊿
LbIE
说明：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为导线未知点坐标平差值