最小半径处坐标X=494380.677 Y=477971.271 最小半径R=40
最小半径处方位角：W=358°28′45.9″ 缓和曲线旋转常数P=65 与之相连的圆曲线长度I=39.216 最小半径处桩号J= W=358°28′45.9″
此段数据计算范围可以从K0+78.504到匝道结束。
当桩号L位于K0+000~ K0+78.504时，曲线参数按第一段编写，当L>78.504时，按第二段缓和曲线编写，当然第二段缓和曲线计算范围也可以从圆曲线起点K0+039.289到K0+78.504中间任意桩号开始。下面是D匝道的数据程序文件（文件名“DZD”）：
3→DimZ：X0”?A:”Y0”?B:“Ln”?L: Abs(L)-Int(Abs(1000L))/1000→Z[3]:Lbl 0:If Z[3]≠0.0001: Then ?L: Else “Xp”?X:”Yp”?Y:X→Z[1]：Y→Z[2]：IfEnd:Lbl 1:If L≥0 :Then 494347.926→X：477989. 873→Y：40→R：302°18°26°→W：55→P:39.216→I: 39.289→J：IfEnd: If L≥78.504：Then 494380.677→X：477971.271→Y：40→R：358°28°45.9°→W：-65→P:39.216→I: 78.504→J：IfEnd:Prog”FDCH”: Prog”XLJS”: If Z[3]≠0.0001:Then Goto 0:Else (Z[1]－X)cos(Ｏ)＋(Z[2]－Y)sin(Ｏ)→N:L＋N→L:－(Z[1]－X)sin(O)+(Z[2]－Y)cos（O）→K：If Abs（N）≥0.001：Then Goto 1：Else “L=“:L◢ “K=”:K◢ IfEnd:Goto 0: IfEnd
在上面程式中，每个缓和曲线的线元参数放在一个If ? Then ?? IfEnd循环语句中，If 后面紧跟计算的起点范围，终点范围从下段的起点范围结束。
变量定义：
1、X、Y 输入缓和曲线半径最小处的X、Y坐标
2、R 缓和曲线最小半径，当线路左转时，半径数据前面加“-”号，以判断线路转向。 3、W 输入缓和曲线半径最小处的方位角
4、P 缓和曲线旋转常数，第一缓和曲线长度为正数，第二缓和曲线常数为负数。当不设缓和曲线时，P=0
5、I 与缓和曲线最小半径端连接的圆曲线长度。 6、J缓和曲线小半径段的桩号。
再看看E匝道数据文件编辑
E匝道线元设计数据表
直线 EK0+000.00 EK0+071.296 EK0+143.296 EK0+203.366 EK0+248.491 EK0+313.333 右偏 右偏 右偏 左偏 左偏 线元名称 起点桩号 偏向 半径或缓和曲线旋转常数A 120 200 95 75 86.75 线元起点坐标、方位角 线元长度(m) 71.296 72 60.070 45.125 64.841 17.800 X(m) 494360.619 494373.205 494381.631 494372.527 494355.878 494337.373 Y(m) 477706.293 477776.470 477847.870 477907.018 477948.932 478010.658 方位角 79°49′57.0″ 79°49′57.0″ 90°08′44.6″ 107°21′16.3″ 113°49′05.6″ 92°24′19.0″
缓和曲线 圆曲线 缓和曲线 缓和曲线 圆曲线 在E匝道线元表中，我们可以看出，E匝道式由一段直线段→缓和曲线段→圆曲线段→缓和曲线段→缓和曲线段→圆曲线段连接起来的。值得注意的是，在位于K0+248.491处，是两段反向的缓和曲线连接点，在编写数据程序文件时，分三段编写，第一段从K0+000到K0+203.366，第二段从K0+203.366到K0+248.919，第三段从K0+248.919到线路终点K0+331.13。具体文件如下：（文件名“EZD”）
3→DimZ：X0”?A:”Y0”?B:“Ln”?L:Abs(L)-Int(Abs(1000L))/1000→Z[3]:Lbl 0:If Z[3]≠0.0001: Then ?L: Else “Xp”?X:”Yp”?Y:X→Z[1]：Y→Z[2]：IfEnd:Lbl 1:If L≥0 :Then 494381.6316→X：477847.870→Y：200→R：90°08°44.6°→W：120→P:60.07→I: 143.296→J：IfEnd: If L≥203.366：Then 494372.527→X：477948.932→Y：200→R：107°21°16.3°→W：-95→P:60.07→I: 203.366→J：IfEnd: If L≥248.491：Then 494337.373→X：478010.658→Y：86.75→R：92°24°19°→W：75→P:17.8→I: 313.333→J：IfEnd: Prog”FDCH”: Prog“XLJS”：If Z[3]≠0.0001:Then Goto 0:Else (Z[1]－X)cos(Ｏ)＋(Z[2]?－Y)sin(Ｏ)→N:L＋N→L:Int(1000L)÷1000+0.0001→L:－(Z[1]－X)sin(O)+(Z[2]－Y)cos（O）→K：If Abs（N）≥0.001：Then Goto 1：Else “L=“:L◢
“K=”:K◢ IfEnd:Goto 0: IfEnd
为了让使用者熟悉掌握数据程序文件的编写线元法的数据程序文件，现在将C、F、G等匝道的线元数据列出来，供使用者练习。 C匝道数据
半径或线元起点坐标、方位角 线元名缓和曲线元长 称 起点桩号 偏向 线旋转度(m) X(m) Y(m) 方位角 常数A 圆曲线 CK0+000.00 右偏 32 56.876 495005.991 478056.677 56°19′42.3″ 直线 CK0+056.876 21.277 494991. 262 478104. 123 158°09′54.4″ 缓和曲线 CK0+078.153 右偏 65 30.179 494971. 511 478112.036 158°09′54.4″ 圆曲线 CK0+108.332 右偏 140 30.431 494943.127 478122.242 164°20′25.7″ 缓和曲线 CK0+138.763 右偏 65 30.179 494913.167 478127.219 176°47′40.6″ 直线 CK0+168.941 右偏 85 422.023 494883.008 478126.740 182.°58′11.8″
缓和曲线 CK0+590.964 右偏 85 85 494461.552 478104.874 182.°58′11.8″ 圆曲线 CK0+675.964 右偏 85 72.739 494379.485 478086.681 211°37′04.2″ F匝道数据
半径或线元起点坐标、方位角 线元名起点桩号 偏向 缓和曲线元长称 线旋转度(m) X(m) Y(m) 方位角 常数A 直线 FK0+000.00 111.525 494591.311 478420.033 252°32′00.2″ 缓和曲线 FK0+111.505 右偏 100 90.909 494557.835 478313.650 252°32′00.2″ 圆曲线 FK0+202.434 右偏 110 101.114 494542.812 478224.689 276°12′33.6″ 缓和曲线 FK0+303.548 右偏 120 130.909 494595.310 478142.420 328°52′35.7″ 数据
G匝道
三、路计算程序
道路计算程序是核心程序，但是道路计算程序是固定的，使用者不需要做任何修改，也不直接运行道路计算程序分为两部分，第一部分是常规的对成型线型组合计算的程序（文件名：XLJS），利用交点法编辑的数据文件可以直接调用；第二部分是适用于非对称的线型组合计算文件（文件名：FDCH），在数据文件调用完“FDCH”文件后，还要调用“XLJS”文件。 线路计算程序（XLJS）
Deg:Fix 4:G÷Abs(G)→V:F2÷(24R)-F^(4)÷(2688R^(3))+F^(6)÷(506880R^(5))-F^(8)÷(154828800R^(7))→P:F÷2-F^(3)÷（240R2）+F^（5）÷(34560R^（4）)-F^（7）÷(8386560R^（6）)+F^（9）÷(3158507520R^（8）)→Q:Q+(R+P)Vtan(G÷2)→T:πRGV÷180+F→I:IF L≤H-T:Then L-H→J:Goto 1:IfEnd:If L Lbl 3:180(J-0.5F) ÷πR→O:Q+Rsin (O)-T→M:(P+R(1-cos (O)))V→N:C+VO→O：Goto B Lbl B:D+Mcos(C)-Nsin(C)→X:E+Msin(C)+Ncos(C)→Y :Goto C:Lbl C:If Z[3]≠0.0001:Then ?K: \=>180-U →U:O+U×Abs(K)÷(K+10^(-10)) →Z:X+Abs(K)×cos(Z)→X◢ Y+Abs(K)×sin(Z)→Y◢ tan-1((Y-B)÷(X-A+10^(-10))→W 线元名称 缓和曲线 圆曲线 缓和曲线 起点桩号 偏向 半径或缓和曲线旋转常数A 175 230 190 线元起点坐标、方位角 线元长度(m) 133.152 149.124 149.72 X(m) 494718.802 494587.600 494472.335 Y(m) 478112.963 478093.368 478002.890 方位角 182°58′11.8″ 199°33′17.5″ 236°42′21.1″ GK0+000.000 GK0+133.152 GK0+282.286 右偏 右偏 右偏 道 X-A<0=>W+180→W:W<0=>W+360→W:√((X-A)2+(Y-B)2)→S \S◢ \W>DMS◢ IfEnd 非对称线型组合计算程序（FDCH） Abs(P+10^(-10)）÷(P+10^(-10))→S：Abs (R)÷R→V:Abs(R)→R:P2÷R→F:180(I+F)V÷(πR)→G:90FSV÷(πR)→O:If S=1:Then W-O→C:Else If S=-1:Then W-O-G→C :IfEnd:IfEnd:If C<0:Then C=C+360→C:IfEndIf C>360:Then C-360→C:IfEnd:F2÷(24R)-F^(4)÷(2688R^(3))+F^(6)÷(506880R^(5))-F^(8)÷(154828800R^(7))→P:F÷2-F^(3)÷(240R2)+F^5÷(34560R^(4))-F^(7)÷(8386560R^6)+F^9÷(3158507520R^(8))→Q:Q+(R+P)Vtan(G÷2)→T:If S=1:Then J-F+T→H:C→Z:If S=-1:Then J-I-F+T→H:C+G→Z:IfEnd:IfEnd:(F-F^(3)÷(40R2)+F^(5)÷(3456R^(4))-F^(7)÷(599040R^(6))+F^(9)÷(175472640R^(8))-T)S→M:(F2÷(6R)-F^(4)÷(336R^(3))+F^(6)÷(42240R^(5))-F^(8)÷(9676800R^(7))+F^(10)÷(3530097000R^(9)))V→N:X-Mcos(Z)+Nsin(Z)→D:Y-Msin(Z)-Ncos(Z)→E 四、程序执行 程序运行分为两种，第一种是根据根据使用者输入的桩号和距中距离计算出所求点的坐标和测站点到该点的方位角和水平距离，可以直接利用计算器提供的距离和方位角进行极坐标放样。另一种是根据使用者输入的实测坐标计算出该点的桩号和距中距离，可以根据桩号和距中距离进行高边坡、隧道断面点的设计高程、以及进行路面高程控制计算方面很方便。 （一）、根据桩号和距中距离计算所求点的坐标、及方位角和水平距离。 要进行程序的线路计算，直接运行数程序文件就可以了，不需要理会计算程序。如果计算某桩号L距离中线K处的坐标X，Y和该点到测站点的距离和方位角，运行提示如下： X0? 提示输入测站点X坐标 Y0？ 提示输入测站点Y坐标 （测站点的位数和数据程序文件里坐标取位要一致）