2024年4月8日发(作者：咸丝祺)
-------------
附闭合导线计算
１、 源程序
F1 A1
L1 Defm 4N-2
L2 N:A:B:Pol(C-A,D-B):W<0=>W=W+360?T=W
L3 K=0=>M=T+180:E=C:F=D:GOTO 0: ≠>E:F:Pol(G-E,H-F):W<0=>W=W+3
60?M=W
L4 Lbl 0:L=0:U=0:I=0:R=2:Z[1]=T
L5 Lbl 1:{J}:Z[R]+360: ?R=N+1=>GOTO 2: ≠>R=R+1:GOTO 1
L6 Lbl 2:P”JB”=(Z[N+1]-M) ?Q”JL”=40√N?R=2
L7 Lb1 3:{S}:Z[N+R]=S:L=L+S?
L8 Z[2N-1+R]=Rec(S,(Z[R]-P(R-1)/N)):U=U+V
L9 Z[3N-2+R]=W:I=I+W:N=R=>GOT 4: ≠>R=R+1:GOTO 3
L10 Lbl 4:P=U+C-E?Q=I+D-F?
L11 G”1:M”=L/Pol(P,Q) ?R=2
L12 Lbl 5:X”XI”=C+Z[2N-1+R]-PZ[N+R]/L?Y”YI”=D+Z[3N-2+R]-QZ[N+R]/L?
L13 R=N=>GOTO 6: ≠>R=R+1:C=X:D=Y:GOTO 5
L14 Lbl 6:”END”
２、 说明
（1）、本程序可计算附和导线和闭合导线的坐标，计算的坐标系经过角度闭合
差及坐标增量闭合差分配后的结果，能显示角度闭合差、增量闭合差及导线全长
的相对精度；
（2）、输入的观测角为导线的左角。
3、程序代号注释
N?导线观测角的折角数；
A、B?导线起始点所后视的已知点的坐标x,y；
C、D?导线起始点（即设站点）的坐标x,y；
E、F?导线终点（已知点）的坐标x,y；
G、H?在导线终点设站观测前视已知点的坐标x,y；
T?起始站后视至起始点的方位角；
M?终点站至前视已知点的方位角；
J?观测的左角值；
JB?角度闭合差；
JL?允许的角度闭合差，程序中是以40√n计算的，如和要求的不一致，可改一
下L6语句中的有关部分。
S?所测导线的边长；
L?边长的累计数；
U?△x的累计数；
I?△y的累计数；
P?x坐标的闭合差；
Q?y坐标的闭合差；
K?转换符，当K=0时为计算闭合导线，当K≠0（任意数）时为计算附和导线。
-------------
-------------
面积计算（多边形法）
1、源程序
F1 A2
L1 N:P=A:Q=B:S=0:I=2
L2 Lbl 0:{C,D}:F=(A+C)(B-D):S=S+F
L3 A=C:B=D:I=I+1
L4 I≦N=>GOTO 0⊿
L5 F=(C+P)(D-Q):S=S+F:S”W”=S/2◢
３、 说明：
（1）、本程序适用于所测断面为多边形闭合图形的面积计算。
（2）、折点坐标按顺时针方向输入，得出的面积为正，否则为负，绝对值是一样
的。
４、 程序代号注释
A、B—计算面积起始点纵横坐标；
C、D—各转折点的纵横坐标；
S—代表计算过程中的有关面积；
S“W”—为图形最后需要的计算面积。
N—多边形的折点个数。
体积计算
1、源程序
F1 A3
L1 J=0:H=0:WG
L2 Lbl 0:{NAB}:NAB:P=A:Q=B:S=0:I=1
L3 Lbl 1:{CD}:S=S+(A+C)(B-D)/2:A=C:B=D:I=I+1
L4 I
L5 J≠1=>GOTO 2: ≠>L=G-H:V=(R+S+√（R*S))*L/3◢⊿W=W+V◢
L6 Lbl 2:R=S:H=G:J=1:{G}:G:GOTO 0
2、说明
（1）、程序可自动计算每一断面面积，当进行到第二个断面时就会显示出1~2断
面间的体积，而后再进行第三断面面积计算，并累计出1~3断面之间的体
积。。。。。。，直到最后得出需算断面间的总体积。
（2）、坐标输入时，应按顺时针方向逐个输入折点坐标，这样得出的面积为正
值，一个桩号折点输入完后，程序自动进入下一桩号的输入状态。
3、程序代号注释
G—断面桩号；
A、B—断面起算折点的坐标；
C、D—断面上其他折点坐标；
S—断面面积；
L—断面间距；
-------------
2024年4月8日发(作者：咸丝祺)
-------------
附闭合导线计算
１、 源程序
F1 A1
L1 Defm 4N-2
L2 N:A:B:Pol(C-A,D-B):W<0=>W=W+360?T=W
L3 K=0=>M=T+180:E=C:F=D:GOTO 0: ≠>E:F:Pol(G-E,H-F):W<0=>W=W+3
60?M=W
L4 Lbl 0:L=0:U=0:I=0:R=2:Z[1]=T
L5 Lbl 1:{J}:Z[R]+360: ?R=N+1=>GOTO 2: ≠>R=R+1:GOTO 1
L6 Lbl 2:P”JB”=(Z[N+1]-M) ?Q”JL”=40√N?R=2
L7 Lb1 3:{S}:Z[N+R]=S:L=L+S?
L8 Z[2N-1+R]=Rec(S,(Z[R]-P(R-1)/N)):U=U+V
L9 Z[3N-2+R]=W:I=I+W:N=R=>GOT 4: ≠>R=R+1:GOTO 3
L10 Lbl 4:P=U+C-E?Q=I+D-F?
L11 G”1:M”=L/Pol(P,Q) ?R=2
L12 Lbl 5:X”XI”=C+Z[2N-1+R]-PZ[N+R]/L?Y”YI”=D+Z[3N-2+R]-QZ[N+R]/L?
L13 R=N=>GOTO 6: ≠>R=R+1:C=X:D=Y:GOTO 5
L14 Lbl 6:”END”
２、 说明
（1）、本程序可计算附和导线和闭合导线的坐标，计算的坐标系经过角度闭合
差及坐标增量闭合差分配后的结果，能显示角度闭合差、增量闭合差及导线全长
的相对精度；
（2）、输入的观测角为导线的左角。
3、程序代号注释
N?导线观测角的折角数；
A、B?导线起始点所后视的已知点的坐标x,y；
C、D?导线起始点（即设站点）的坐标x,y；
E、F?导线终点（已知点）的坐标x,y；
G、H?在导线终点设站观测前视已知点的坐标x,y；
T?起始站后视至起始点的方位角；
M?终点站至前视已知点的方位角；
J?观测的左角值；
JB?角度闭合差；
JL?允许的角度闭合差，程序中是以40√n计算的，如和要求的不一致，可改一
下L6语句中的有关部分。
S?所测导线的边长；
L?边长的累计数；
U?△x的累计数；
I?△y的累计数；
P?x坐标的闭合差；
Q?y坐标的闭合差；
K?转换符，当K=0时为计算闭合导线，当K≠0（任意数）时为计算附和导线。
-------------
-------------
面积计算（多边形法）
1、源程序
F1 A2
L1 N:P=A:Q=B:S=0:I=2
L2 Lbl 0:{C,D}:F=(A+C)(B-D):S=S+F
L3 A=C:B=D:I=I+1
L4 I≦N=>GOTO 0⊿
L5 F=(C+P)(D-Q):S=S+F:S”W”=S/2◢
３、 说明：
（1）、本程序适用于所测断面为多边形闭合图形的面积计算。
（2）、折点坐标按顺时针方向输入，得出的面积为正，否则为负，绝对值是一样
的。
４、 程序代号注释
A、B—计算面积起始点纵横坐标；
C、D—各转折点的纵横坐标；
S—代表计算过程中的有关面积；
S“W”—为图形最后需要的计算面积。
N—多边形的折点个数。
体积计算
1、源程序
F1 A3
L1 J=0:H=0:WG
L2 Lbl 0:{NAB}:NAB:P=A:Q=B:S=0:I=1
L3 Lbl 1:{CD}:S=S+(A+C)(B-D)/2:A=C:B=D:I=I+1
L4 I
L5 J≠1=>GOTO 2: ≠>L=G-H:V=(R+S+√（R*S))*L/3◢⊿W=W+V◢
L6 Lbl 2:R=S:H=G:J=1:{G}:G:GOTO 0
2、说明
（1）、程序可自动计算每一断面面积，当进行到第二个断面时就会显示出1~2断
面间的体积，而后再进行第三断面面积计算，并累计出1~3断面之间的体
积。。。。。。，直到最后得出需算断面间的总体积。
（2）、坐标输入时，应按顺时针方向逐个输入折点坐标，这样得出的面积为正
值，一个桩号折点输入完后，程序自动进入下一桩号的输入状态。
3、程序代号注释
G—断面桩号；
A、B—断面起算折点的坐标；
C、D—断面上其他折点坐标；
S—断面面积；
L—断面间距；
-------------