5800万能隧道程序

本文格式为Word版，下载可任意编辑——5800万能隧道程序

5800万能隧

道程序

1,主程序

Lbl0：\：\：？T：？S：Prog“SUB0〞↙

1÷P→C:(P-R)÷(2HPR)→D:180÷∏→E:T=1=>Goto1：Goto2：↙

Lbl1:“LICHENG=〞？K：“BANJNG=〞？P：“PIANJIN=〞？S：

：Abs(S-O)→W：Prog\↙Prog\：G+√ˉ(P2－S2)+1.66◢\：X◢\：Y◢\：Z◢

F-90→F：“FS=〞：F?DMS◢Goto0

Lbl2：Prog“SUB2〞↙Prog“YFS3〞Goto0↙

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

3.反算子程序(SUB2)

Lbi1∶“BANJING〞？P：“X〞？X：“Y〞？Y：“Z〞？Z：X→I：Y→J：G-90→T：

Abs((Y-V)cos(T)-(X-U)sin(T))→W：0→Z：Lbl4：Prog\：

T+QEW(C+WD)→L：(J-Y)cos（L）-(I-X)sin（L）→Z：IFAbs（Z）

Else“SHANG〞∶W◢Return∶Stop

11子程序(GC)高程

IfK≥93560AndK＜93805∶Then504.33+2.4%(K－93560)→G∶Return∶IfEnd↙

IfK≥93805∶Then510.220+2.4%(K－93805.399)－(K－93805.399)2÷2÷18504→G∶Return∶IfEnd↙

SUB0数据库子程序

Goto1↙同时保存多个曲线时的指针Lbl1：IFSXY2.XY=>SZ

N?选择计算方式，输入1表示进行由

里程、边距计算坐标；输入2表示由坐标反算里程和边距。S

？正算时所求点的里程，反算时为所求点的近似里程，反算时输入的近似里程不能夸线元，若夸线元的话，需进行其次次反算，并以第一次的结果做为近似里程，否则反算的里程及边距不对，Z？正算时所求点距中线的边距(左侧取负，值右侧取正值，在中线上取零)X？反算时所求点的X坐标Y？反算时所求点的Y坐标显示部分：

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

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

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

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

Z=×××反算时，计算得出的所求点的边距曲线任意里程中边桩坐标正反算(CASIO

fx-4800P计算器)程序一、程序功能

本程序由一个主程序(TYQXJS)和两个子程——正算子程序(SUB1)、反算子程序(

SUB2)序构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线

元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲

率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。另

外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲

线任意里程中边桩坐标进行正反算。本程序也可以在CASIOfx-4500P计算器及CASIOfx-4850P计算器上运行。

特别申明：(1).适用于弧长小于2倍半径的各种线元坐标正反算，精度优