4800平曲线坐标反算桩号程序（对称与非对称缓和曲线）

4800-对称与非对称缓和曲线坐标计算程序

主程序：XY=NH

\？M ：\？N ： 0→T

LbI 1：Prog\： B∧2/24/A→P：B/2-B∧3/240/A∧2→Q：K- Z[1]→C ：C+B→D ：D+(AbsW-90(B+U)/π/A) πA /180→E：E+U→F If W<0 : Then S= -1：Else 1→S:IfEnd

LbI 6：If L≤C: Then K-L→O： G-Ocos V→R： H-Osin V→U： V→Z：Goto 3:IfEhd

IF L≥F: Then L-F→O: V+W→Z： G+Ocos Z→R： H+Osin Z→U：Goto 3: IfEND

IF L

IF L>E:Then F-L→O ：U→B:Prog\：V+W→O ：O-SR→Z： G+（Z[2]-I）cos O-SJsin O→R ：H+（Z[2]-I）sin O+JScos O→U：Goto 3:Ifend 180（L-D+B/2）/ π/A→O： Asin O+Q→I： P+A（1-cos O）→J： V+SO→Z：G-（Z[1]-I）cos V-SJsin V→R ：H-（Z[1]-I）sin V+SJcos V→U LbI 3：？T：If T=1:Then Goto 5 ：Else 0→O : Ifend \◢ \◢ \◢

LbI 4：\\？O：If O=0 : Then Goto 1：IFEnd If O=-1:Then Goto 5 : Ifend

R-Osin Z→X ： U+Ocos Z→Y ：0→J :0→I：Pol（X-M，Y-N）：J<0=> J+360→J:Ifeng \\：X◢ \◢ \◢ \◢ Goto 4:Ifend

LbI 5： \\？X：\ ？Y If X=-1 : Then 0→T：Goto 4 : Ifend

Pol（X-R，Y-U） ： Icos（J-Z）→O ：If AbsO≤0.001: Then L+ O→L◢ Isin（J-Z）→Z[3] \\◢

Goto 5 : Else T=1：L=L+O : Goto 1 : Ifend

子程序：ZH--HY

AB→U： O-O∧5/40/U∧2→I： O∧3/6/U-O∧7/336/U∧3→J： 90O∧2/πU→R

数据库子程序：JD-X

？L:L≥45250.23 And L＜46500.29: Then 45635.294→K : 3303516.326→G: 34635482.418→H : 154°23′25.3″→V : -23°48′12.4″→W : ？？？？→ z[1] : ？？？→z[2] : 1100→A : 240→B ：250→U : Ifend ？L＞.......... ：Ifend：Ifend

？

使用过程“显示”说明：

1、 M\、N\为输入测站点坐标。

2、 K、G、H、V、W、Z[1]、Z[2]、A、B、U分别为：输入交点桩号、交点X坐标、交点Y坐标、起始方位角、偏角（左-右+）、第一切线长度、第二切线长度、半径、第一缓和曲线长度、第二缓和曲线长度。

3、 L、X0、Y0、I分别为：输入桩号，输出中桩X坐标、中桩Y坐标、即时方位角。

4、 S、X1、Y1、G、L分别为输入边桩偏距（左-右+）、输出边桩X坐标、边桩Y坐标、边桩放样方位角与距离。

（注:当S= -1 时，程序为实测点坐标反推该点里程桩号计算；当XF= -1 时，程序为返回中、边桩坐标计算；当S= 0 时，程序返回里程推算坐标计算。）

5、 XF、YF、分别为输入实测点X坐标、Y坐标。 6、 L+O、SF为输出实测点桩号、偏中距离。