卡西欧5800坐标正反算程序

卡西欧5800坐标正反算程序

1.主程序(TYQXJS)

Lb1 4:“1.SZ=>XY”：“2.XY=>SZ”:？N:？S:Prog“SUB0”↙

1÷P→C：（P-R）÷（2HPR）→D：180÷π→E:N=1=>Goto 1:Goto 2: ↙

Lb1 1:?Z:Abs(S-O) →W: Prog“SUB1”:“XS=”：X◢

“YS=”:Y◢

F-90→F:“FS=”F▲DMS◢ Goto 4↙

Lb1 2:?X:?Y:X→I:Y→J:Prog“SUB2”：O+W→S:“S=”：S◢

“Z=”:Z◢Goto 4↙

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+QEFW(C+FWD))+Acos(G+QEMW(C+MW D)))

→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↙

3. 反算子程序（SUB2）

G-90→T:Abs((Y-V)cos(T)-(X-U)sin(T)) →W:0→Z:Lb1 0:Prog“SUB1”:T+QEW(C+WD) →L:(J-Ycos(L)-(I-X)sin(L) →Z:If Abs(Z)〈 0.000001：Then 0→Z:Prog“SUB1”:(J-S)÷sin(F) →Z:Else W+Z→W:Goto 0:IfEnd↙

4. 数据库子程序（SUB0）

Goto1↙同时保存多个曲线时的指针

Lb1 1:If S〈 000.000（线元终点里程）：Then 000.000→G(线元起点方位角)：000.000→O （线元起点里程）：000.000→U(线元起点X坐标)：000.000→Y（线元起点Y坐标）：000.000→P(线元起点曲率半径)：000.000→R（线元终点曲率半径）：000.000→H（线元起点至终点长度）：0或1、-1→Q：Eeturn:IfEnd↙

Lb1 1:If S〈 000.000（线元终点里程）：Then 000.000→G(线元起点方位角)：000.000→O （线元起点里程）：000.000→U(线元起点X坐标)：000.000→Y（线元起点Y坐标）：000.000→P(线元起点曲率半径)：000.000→R（线元终点曲率半径）：000.000→H（线元起点至终点长度）：0或1、-1→Q：Eeturn:IfEnd↙

………………….

为了便于解读，每增加一个线元增加一行语句，每增加一条曲线增加一个Lb1，每增加一个工程增加一个文件。数字0和大写字母O分不清楚的，用查找命令即可。

三、使用说明

1、规定

（1）以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0.

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

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

（4）当线元为圆曲线时，无论其起点、止点于什么线元相接，其曲率半径均等于圆弧的半径。

（5）当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次方代替，与圆曲线相接时，曲率半径等于圆曲线的半径，止点于直线相接时，曲率半径为无穷大，以10的45次方代替，与圆曲线相接时，曲率半径等于圆曲线的半径。

（6）当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定值，与圆曲线相接时，曲率半径等于圆曲线半径。止点与直线相接时，曲率半径等于设计规定值，与圆曲线相接时，曲率半径等于圆曲线半径。

2、输入与显示说明

输入部分：

1.SZ=>XY

2.XY=>SZ

N?选择计算方式，输入1表示进行由里程、边距计算坐标；输入2表示由坐标反算里程和边距。

S?正算时所求点的里程，反算时为所求点的近似里程，反算时输入的近似里程不能跨线元，若跨线元的话，需进行第二次反算，并以第一次的结果作为近似里程，否则反算的里程及边距不对，Z？正算时所求点距中的边距（左侧为负，右侧为正，在中线上为零）

X？反算时所求点的X坐标

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

显示部分：

XS=000.000正算时，计算得出的所求点X坐标

YS=000.000正算时，计算得出的所求点Y坐标

FS=000.000正算时，所求点对应的中线点的切线方位角

S=000.000反算时，计算得出的所求点的里程

Z=000.000反算时，计算得出的所求点的边距