CASIO5800线元法程序

卡西欧5800P 线元法正反算程序

1-JIN YU(主程序)

Lbl 4: “1→ZS,2→FS?N(选择计算模式，1为正算，2为反算，3为超欠挖)

N=1=>Goto 1:N=2=>Goto 2

Lbl 1: “K=”?S：“P=”?Z：Prog “SJ-PM”： Abs(S-O) → W：Prog “SUB1-ZS”：“X=”：Locate4,4,X:“Y=”：Locate4,4,Y:F-90→F: “W=”: Locate4,4,F°◢

Goto 4

(正算-输入待求点里程K=、输入待求点偏距P=、显示待求点里程X=、显示待求点里程Y=、显示待求点里程设计高程H=)

Lbl 2:“X=”?X：“Y=”？Y:Prog“SJ-PM”： X→ I： Y→J：Prog "SUB2-FS"：O+W→S: “K="：Locate4,4, S：“P="：Locate4,4, Z：S→ K:Prog“SJ-GC”：“H=”：Locate4,4, H:F-90→F: “W=”: Locate4,4,F°◢

Goto 4

(反算-输入实测点X=、输入实测点Y=、显示实测点里程K=、显示实测点偏距P=、显示实测点的设计高程H=)

SJ-PM（子程序名-平面线形数据库）

If S ≥59227.681（线元起点里程）:Then 2599818.013→U（线元起点X坐标）：496887.918→V（线元起点Y坐标）：59227.681→O（线元起点里程）：208028’39”→G（线元起点方位角）：100 →H（线元长度）：1×1045→P（线元起点曲率半径）：1×1045→R（线元终点曲率半径）：0→Q（线元左右偏标志：左-1右1直0）：IfEnd

If S ≥59327.681（线元起点里程）:Then 2599730.112→U（线元起点X坐标）：496840.237→V（线元起点Y坐标）：59327.681→O（线元起点里程）：208028’39”→G（线元起点方位角）：90 →H（线元长度）：1×1045→P（线元起点曲率半径）：750→R（线元终点曲率半径）：-1 →Q（线元左右偏标志：左-1右1直0）：IfEnd

(注：有多个平面线形，依照上面的依次变更, 每多一个，就增加一个。每次只需要修改以上的数据或增加一个判断。)

SUB1-ZS(正算子程序)

1÷P→ C： (P-R)÷(2HPR) → D： 180÷π→ E： 0.1739274226→A： 0.3260725774→ B：0.0694318442→ K： 0.3300094782→ L： 1-L→ F：1-K→ M：

U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Bcos(G+QEFW(C+FWD))+Acos(G+QEMW(C+MWD) )) → X：

V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+QEFW(C+FWD))+Asin(G+QEMW(C+MWD) )) → Y： G+QEW(C+WD)+90→ F： X+Zcos(F)→ X： Y+Zsin(F)→ Y

SUB2-FS(反算子程序)

G-90→T： Abs((Y-V)cos(T)-(X-U)sin（T）) → W： 0→ Z：Lbl 0：Prog "SUB1-ZS"： T+QEW(C+WD)→L： (J-Y)cos(L)-(I-X)sin(L) → Z：If Abs(Z)<1E-6(注：此处1E-6可输入0.000001)：then Goto1：Else W+Z → W：Goto 0：IfEnd

Lbl 1： 0→ Z：Prog "SUB1-ZS"： (J-Y)÷sin(F)→ Z

说明：(正算1秒，反算和超欠挖5秒完成)

所有程序名。不同线路，只需改动SJ-PM，SJ-GC，SJ-CQW三个子程序里的内容，其它不变。

1为正算，2为反算，3为超欠挖.

K= 正算时，输入所求点里程。反算时得出里程结果

X=，Y= 正算时得出结果。反算时输入实测坐标

P= 正算时输入偏距。反算时得出偏距

H= 正算时得出高程。反算时得出高程

H1= 超欠挖时输入实测高程

CQ→Y?= 得出超欠挖值及第几圆心上