一、源程序(SZ-XY-SZ)
"1.SZ => XY"："2. XY => SZ"：NU"X0"V"Y0"O"S0"G"F0"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：X"XS"
= U+Rec(C,F)+ZcosE◢
Y"YS"=V+J+ZsinE◢
Goto 1←┙
Lbl 2：W=G+90Q：E=U+Rec(R,W)：F=V+J←┙
Lbl 3：{XY}：XY：D=Pol(X-E,Y-F):T=J：J<0=>T=J+360△B=T-W-180M：P=Abs B：S"S"=O+
DKcosP+ MπRP÷180◢
Z"Z"=DKsinB+MQ(R-D)◢
Goto 3←┙
二、使用说明
1.由里程、中边桩计算坐标
启动程序后按以下提示输入数据：
N ？1 （N=1，由里程中边桩计算坐标；N≠1，由坐标反算里程及中边桩)
X0 ？142.353 (直线或圆曲线起点的X坐标）
Y0 ？368.166 (直线或圆曲线起点的Y坐标）
S0 ？155.252 (直线或圆曲线起点的里程）
F0 ？60°54′03″(直线或圆曲线起点的切线方位角）
Q ？-1 (当线元为直线时Q=0；当线元为圆曲线且左偏时Q=-1；当线元为圆曲线且
右偏时Q=1)
R ？ 100 (圆曲线的半径)
S ？ 194.016 (中线里程)
Z ？ -5.251 (中线左侧5.251m的边桩；在中线上Z=0，在左侧<0，在右侧Z>0)
显示计算结果
XS=172.7467767
YS=395.2316807
2.由坐标反算里程、中边桩
启动程序后按以下提示输入数据：
N ？2 （N=1，由里程中边桩计算坐标；N≠1，由坐标反算里程及中边桩)
X0 ？142.353 (直线或圆曲线起点的X坐标）
Y0 ？368.166 (直线或圆曲线起点的Y坐标）
S0 ？155.252 (直线或圆曲线起点的里程）
F0 ？60°54′03″(直线或圆曲线起点的切线方位角）
Q ？-1 (当线元为直线时Q=0；当线元为圆曲线且左偏时Q=-1；当线元为圆曲线且
右偏时Q=1)
R ？ 100 (圆曲线的半径)
X ？ 172.7467767 (待求里程及中边桩的点的X坐标)
Y ？ 395.2316807 (待求里程及中边桩的点的Y坐标)
显示计算结果
S=194.016
Z=-5.251
三、注意事项
1.本程序主要是针对市政道路工程中的直线和圆曲线而编制的，未考虑缓和曲线等
曲线的计算；
2.程序中“R=1E-45”的意思为10的-45次方，输入程序时切莫输错。
