卡西欧5800测量坐标计算程序

一．使用说明：

1.规定：

（1）以道路中线的前进方向区分线路的左右偏转方向，线路左偏Q=-1，线路右偏Q=1，直线Q=0。

（2）当所求点位于线路中线时，Z=0；当位于左侧时，Z取负值；当位于右侧时，Z取正值。

（3）当线元为直线时，其起终点的曲率半径为无穷大，取10的45次方代替。

（4）当线元为圆曲线时，起终点曲率半径为圆曲线半径。

（5）当线元为缓和曲线时，起终点曲率半径与直线相接取10的45次方代替，与圆曲线相接取圆曲线半径。

2.输入与显示说明：

（1）输入部分：

N？：选择计算方式，输入1表示由里程边距计算坐标，输入2表示由坐

标计算里程边距。

X0？线元起点X坐标。

Y0？线元起点Y坐标。

S0？线元起点里程。

F0？线元起点切线方位角。

LS？线元长度。

R0？线元起点曲率半径。

RN?线元终点曲率半径。

Q？线元左右偏标志（左偏=-1，右偏=1，直线=0）

S?正算时所求点里程。

Z？正算时所求点距中线边距。位于线路左侧输入负值，在右侧输入正

值。

AR?斜交角度。（为线路右角，在0-180度之间）

X?反算时所求点X坐标。

Y？反算时所求点Y坐标。

（2）显示部分：

XS=计算结果，所求点X坐标。

YS=计算结果，所求点Y坐标。

FS=计算结果，所求点切线方位角。

S=计算结果，反算时所求点里程

Z=计算结果，反算时所求点边距。

二．程序清单1（用于任意曲线元的坐标计算）

1.主程序(TYQXJS)

″1.SZ=>XY″:″2.XY=>SZ″:?N:″X0″?U:″Y0″?V:″S0″?O: ″F0″?G:″LS″?H:″

R0″?P:″RN″?R:?Q:1÷P→C:(P-R) ÷(2HPR) →D:180÷π→E:If N=1: THEN Goto 1:Else Goto 2:IfEnd:Lbl 1:?S:”D(Z-，Y+)=”?Z:90→T:”O=”?T:Abs(S-O) →

W:Prog ″SUB1″:″XS″:X◢″YS″:Y◢″FS″:F◢

Goto 1:Lbl 2:?X:?Y:I→X:J→Y:Prog ″SUB2″:″S″:O+W→S:◢″Z″:Z◢Goto 2

2.正算子程序(SUB1)

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+QEF

W(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) →F:X+Zcos(F+T) →X:Y+Zsin(F+T) →Y

3.反算子程序(SUB2)

G-90→T:Abs((Y-V)cos(T)-(X-U)sin(T)) →W:0→Z:Lbl 0:Prog ″

SUB1″:T+QEW(C+WD) →L:(J-Y)cos(L)-(I-X)sin(L)→Z:IF Abs(Z)<1E-6 :Then Goto 1:Else W+Z→W:Goto 0 : IfEnd:Lbl 1:0→Z:Prog ″SUB1″:(J-Y)÷sin(F)→Z 三．用于提前将曲线元要素输入子程序，计算时只输入里程桩号与边距。

1.主程序(TYQXJS)(本程序可将线路要素提前输入计算器)

Lbl 1:?S:

IF S>=起点桩号AND S<终点桩号:THEN PROG “YS1”: IFEND:

IF S>=起点桩号AND S<终点桩号:THEN PROG “YS2”: IFEND:

IF S>=起点桩号AND S<终点桩号:THEN PROG “YS3”: IFEND:

……

1÷P→C:(P-R) ÷(2HPR) →D:180÷π→E:”D(Z-，Y+)=”?Z:90→T:”O=”?T:Abs(S-O) →W:Prog ″SUB1″:″XS″:X◢″YS″:Y◢″FS″:F◢

Goto 1

2.曲线要素子程序（第一段曲线）(YS1)

第一段曲线起点X坐标数值→U：第一段曲线起点Y坐标数值→V：第一段曲线起点里程数值→O：第一段曲线起点切线方位角数值→G：第一段曲线长度数值→H：第一段曲线起点半径数值→P：第一段曲线终点半径数值→R：第一段曲线偏转标志数值→Q

3. 曲线要素子程序（第二段曲线）(YS2)

第二段曲线起点X坐标数值→U：第二段曲线起点Y坐标数值→V：第二段曲线起点里程数值→O：第二段曲线起点切线方位角数值→G：第二段曲线长度数值→H：第二段曲线起点半径数值→P：第二段曲线终点半径数值→R：第二段曲线偏转标志数值→Q

4. 曲线要素子程序（第三段曲线）(YS3)

第三段曲线起点X坐标数值→U：第三段曲线起点Y坐标数值→V：第三段曲线起点里程数值→O：第三段曲线起点切线方位角数值→G：第三段曲线长度数值→H：第三段曲线起点半径数值→P：第三段曲线终点半径数值→R：第三段曲线偏转标志数值→Q

……

Lbl 1:?S:

IF S>=起点桩号AND S<终点桩号:THEN PROG “ZS1”: IFEND:

IF S>=起点桩号AND S<终点桩号:THEN PROG “ZS2”: IFEND:

IF S>=起点桩号AND S<终点桩号:THEN PROG “ZS3”: IFEND:

……