隧道测量-快速坐标正反算程序
(4800-4850计算器）
ZB （坐标正算）
1. Lb1 1 2. ｛L｝
3.SMNARCL 4. K=L-S
5. R=0=>X“X”=M+KcosA▲ 6. Y“Y”=N+KsinA▲ 7. Z“FW”=A▲
8.Goto 4△R≠0=>Abs K≤C=>Goto 2△R≠0=>K＞C=>Goto 3△
R≠0=>K<0=>Abs K＞C=>Goto 1△ 9. Lb1 2
10. O=90K²÷（πRC）
11. G=K³÷（6RC）-Kˆ7÷（336（RC）³＋Kˆ11÷（42240（RC）ˆ5-Kˆ15÷9676800（RC）ˆ7）＋Kˆ19÷（35300960（RC）ˆ9）
12. F=K－Kˆ5÷（（40（RC）²）－Kˆ9÷（3456（RC）ˆ4）＋Kˆ13÷（599040（RC）ˆ6）－Kˆ17÷（40320×4532（RC）ˆ8）
13. J=0=>X“X”=M＋GsinA＋cosA▲
14. Y“Y”=N＋FsinA- cosA▲ 15. Z“FW”=A-O ▲ Goto 4△
16．J≠0=>X“X”=M＋F cosA－GsinA▲ 17．Y“Y”=N＋G cosA＋F sinA▲ 18．Z“FW”=A＋O▲ Goto 4△ 19．LbI 3
20．B=90（2K-C）÷（лR）
21.U=C÷2－C³÷（240R²）＋R sinB＋Cˆ5÷（34560Rˆ4）-Cˆ7÷（599040 Rˆ6）
22．V=C²÷（24R）＋R（1－cosB）－Cˆ4÷（2688R³）＋Cˆ6÷（42240 Rˆ5）
23．J=0=>X“X”=M＋UcosA＋VsinA▲ Y“Y”=N＋U sinA－V cosA▲ Z“FW”=A-B▲ Goto 4△
24．J≠0=>X“X”=M＋U cosA－VsinA▲ Y“Y”=N＋U sinA＋V cosA▲ Z“FW”=A＋B▲ Goto 4△ 25．Lb1 4
26．｛D｝：｛T｝ 27．D=0 => Goto 1△
28．D≠0=>X“LX”=X＋Dcos（T＋Z）▲ Y“LY”=Y＋Dsin（T＋Z）▲ Goto 1△
SUB（坐反算）
“1.SZ=>UV”: “2.UV =>SZ”: N（=2): X“X0”: Y“Y0”: O“S0”=L: G“F0”(=Z):
Q：Q=0=>R=1E-45：≠> R△J=0：M=AbsQ：K=1-M：N=1=>Goto 1：≠>Goto 2△←┙
Lbl 1：{SZ}：SZ：L=S-O：B=90QL÷π÷R：F=G+B：E=F+B+90：C=LK+2MRsinQB：U\"US\" = X+Rec(C,F)+ZcosE◢ V \"VS\" =Y+J+ZsinE◢ Goto 1←┙
Lbl 2：W=G+90Q：E=X+Rec(R,W)：F=Y+J：W=W+180M：W>360=>W=W-360△←┙
Lbl 3：{UV}：UV：D=Pol(U-E,V-F):T=J：J<0=>T=J+360△B=T-W：P=KB+MAbs B：S\"S\"=O+ DKcosP+ MπRP÷180◢ H\"H\"=DKsinB+MQ(R-D)◢ Goto 3←┙
正算由下面输入： S： ZH或HZ点里程 M： ZH或HZ点X坐标 N： ZH或HZ点Y坐标 A： 切线方位角 R： 曲线半径 C： 缓和曲线长度 L： 求点里程 X： 求点X坐标 Y： 求点Y坐标 FW： 求点切线方位角 D： 左、右边桩距离 T： 边桩与线路夹角 LX： 边桩X坐标 LY： 边桩Y坐标
J： 右偏输“1”，左偏输“0” 反算由下面输入： X0: 线路中心X坐标 Y0: 线路中心Y坐标 L: 相对应里程 Z: 相对应方位角 Q: ＝０
U: 求点X坐标 V: 求点Y坐标 S= 求点里程 H= 求点边距
说明：１.可用于缓和曲线坐标正反算，反算时X0,Y0,L,Z利用正算储存的数据，速度快． ２.反算时先输入近似里程正算．
３.隧道断面炮位放样２０多个点位在操作熟练时只 需15分钟左右，提高工程进度．
