作者：檸檬樹???? QQ：249910569

1.坐標正算主程序(命名為ZBZS)

第1行：Lbl 0:{W,T,B}:W”K=”:T”LEN=”:B

第2行：Prog “A”

第3行：X=D+Tcos(F+B)◢

第4行：Y=E+Tsin(F+B)◢

第5行：F=F◢

第6行：Goto 0

K——計算點的里程

LEN——計算點到中樁的距離（左負右正）

B——取前右夾角為正

2.坐標反算樁號和偏距主程序(命名為ZBFS)

第1行：{U,V,K}:U”X1”:V”Y1”: W”K1”

第2行：I=0:J=0

第3行：Lbl 0:Prog “A”

第4行：Pol(U-D,V-E):S=Icos(F-J):W=W+S

第5行：Abs(S)>0.0001=>Goto0△

第6行：T=Isin(J-F)

第7行：K=W◢

第8行：T”LEN”=T◢

X1——取樣點的X坐標

Y1——取樣點的Y坐標

K1——輸入時為計算起始點(在線路內即可)，輸出時為反算點的樁號

T——偏距（左負右正）

3.計算坐標子程序(命名為XYF)

為了簡潔，本程序由數據庫直接調用，上述中的正反算主程序不直接調用此程序

第1行：S=W-A：I=(Q-P)÷L

第2行：D=D+∫(cos(F+X(2P+XI)×90÷π),0,S,4)

第3行：E=E+∫(sin(F+X(2P+XI)×90÷π),0,S,4)

第4行：F=F+S(2P+SI)×90÷π

4.數據庫（命名為A）

第1行：W≤175.191=> GoTo2△

第2行：A=175.191:D=428513.730:E=557954.037:F=92°26′40″:P=0:Q=1/240:L=70.417:

W≤A+L=>GoTo 1△

第3行：A=245.607: D=428507.298:E=558024.092:F=100°50′59.4″:P=1/240:Q=1/240:

L=72.915: W≤A+L=>Goto 1△

第4行：A=318.522: D=428482.988:E=558092.538:F=118°15′25.2″:P=1/240:Q=0:L=

55.104: W≤A+L=>Goto 1△

第5行：A=373.627:D=428453.283:E=558138.912:F=124°50′4.5″0:P=0:Q=-1/180:L=67.222:

W≤A+L=>Goto 1△Goto 2

第6行：Lbl 1:Prog “XYF”: Goto 3

第7行：Lbl 2:D=0:E=0:F=0

第8行：Lbl 3

A——曲線段起點的里程

D——曲線段起點的x坐標

E——曲線段起點的y坐標

F——曲線段起點的坐標方位角

P——曲線段起點的曲率(左負右正)

Q——曲線段終點的曲率(左負右正)

L——曲線段長度（盡量使用長度，為計算斷鏈方便）

說明：

（1）正算主程序可以計算一般邊樁的坐標，如要計算類似涵洞端墻的坐標需增加兩個變量，具體方法參考本程序集中的第1篇辛普生公式的坐標計算通用程序

（2）程序規定，左偏曲線曲率（半徑倒數）輸入負值，右偏曲線曲率輸入正值，直線上點曲率輸入0。

（3）本程序精度極高，不受曲線半徑大小影響，即使極小半徑的螺旋曲線等誤差僅為萬分之一（0.1mm），可以忽略。

（4）若是從大里程向小里程的反方向計算，則曲率取正方向時的負值，方位角減去(或加上)180度。

（5）有多個匝道的項目，可隨時更改正反算主程序中的紅色字體部分來調用其它線路的數據

（6）反算樁號偏差為1mm

（7）可以計算任意線型的任意點坐標

5.坐標計算通用程序(命名為ZB)

本程序與上述中的幾個程序無關，是在計算坐標時采用輸入曲線元參數模式下的程序，為臨時使用方便。

第1行：A”K0“:D”X0“:E”Y0“:F”F0“:P”1÷R1“:Q”1÷R2“:L“L0“

第2行：Lbl 0:{K,T,B}:K”K”:T”LEN”:B

第3行：S=Abs(K-A)：I=(Q-P)÷L

第4行：D=D+∫(cos(F+X(2P+XI)×90÷π),0,S,4)

第5行：E=E+∫(sin(F+X(2P+XI)×90÷π),0,S,4)

第6行：F=F+S(2P+SI)×90÷π

第7行：X=D+Tcos(Z+B)◢

第8行：Y=E+Tsin(Z+B)◢

第9行：F=F◢（可省掉）

第10行：Goto 0