道路中(高程)边桩坐标计算程序A解析
高等级公路中桩边桩坐标计算方法
线路工程测量
14.7 线路逐桩坐标计算
2、坐标反算
根据直线起点和终点的坐标，计算直线的边长和坐 标方位角，称为坐标反算。
AB
arctan
YAB X AB
DAB (XAB )2 (YAB )2
线路工程测量
14.7 线路逐桩坐标计算
三、中桩坐标计算
1、直线上点的坐标计算
xp xJDi1 DK p DK JDi1 cosi1,i y p yJDi1 DK p DK JDi1 sin i1,i
(1)第一缓和曲线及圆曲线上点的坐标计算 当P点位于第一缓和曲线(ZH-HY)上，按切线支距法 公式：
xP
l
l5 40R 2l02
JDi
yP
l3 6 Rl0
JDi1
HY ZH
YH HZ
JDi1
线路工程测量
14.7 线路逐桩坐标计算
由坐标转换公式，P点在线路坐标下坐标：
X P X ZH xP cosi1,i KyP sin i1,i YP YZH xP sin i1,i KyP cosi1,i
曲线右偏时K=1；曲线左偏时K=-1；
JDi
JDi1
HY ZH
YH HZ
JDi1
线路工程测量
14.7 线路逐桩坐标计算
(2)圆曲线上点的坐标计算 当P点位于HY-YH圆曲线上，则：
xP m R sin
yP R P R cos
0
DK P
DK HY
R
1800
JDi
HY ZH
JD
YH HZ
线路工程测量
14.7 线路逐桩坐标计算
§14.7 线路中桩、边桩坐标计算
一、引言
公路逐桩坐标计算程序
2 )sinαAB +(
P 点为顺时针方向时，其方位角为 αAB+900 P 点为逆时针方向时，其方位角为 αAB−900
第一段缓和曲线时：以直缓点（ZH）为起点计算，αAB 为 ZH 点的坐标方位角，L 为 P 点距 ZH 点的距离。 第二段缓和曲线时：以缓直点（HZ）为起点计算，αAB 为 HZ 点坐标方位角的反 方向即 HZ 点方位角加 180 度，L 为 P 点距 HZ 点的距离，加减 90 度刚好与第一 段缓和曲线相反。图如下：
缓和曲线转角公式：β
=
L2
2RLs
（2） 边桩坐标计算公式
左侧
XA=XP+T1cos(αAB± β-900) YA=YP+T1sin(αAB± β-900)
右侧
XB=XP+T2cos(αAB± β+900) YB=YP+T2sin(αAB± β+900)
*第二段缓和曲线计算边桩坐标时， 注意加减 90 度时， 与第一段缓和曲线相反。
2、 缓和曲线段 （1） 中桩坐标计算公式
Xp=X1+(L − Yp=Y1+(L −
L5 40R 2 Ls L5 40R 2 Ls
2 )cosαAB+(
L3 6RLs L3 6RLs
− −
L7 336R 3 ������s 3 L7 336R 3 ������s 3
)sin(αAB±900) )sin(αAB±900)
右侧
XB=XP+T2cos(αAB± β+900) YB=YP+T2sin(αAB± β+900)
第二章
公路导线测量计算
1
公路导线测量为附合导线测量，按路线前进方向测量右角。β
路任意点(中,边桩)坐标高程程序
公路任意点（中，边桩）坐标高程程序500字节（4800，4850）作者：周新全主程序（坐标和高程） ZCX程序为公路坐标高程计算的最最简单适用的程序。
而且可以计算非对称型，不完整的曲线，复曲线，卵形曲线，匝道等多种功能。
利用高等数学（积分和复数）自己编写。
计算精度非常高。
需要加我QQ：179122447。
或者EMAIL：laoye5403@.01 Prog "DATE 1":L=Z-D:V=1:R=90/πBC:D=L(F+.02)/C（ 1.主程序中i为复数的i）,02 L＞L-S=＞L=S-L:V=-1:C=U:G=G+HA:T=OΔ计算第二缓和曲线03 L＜0=＞I=L:J=0:P=0:W=0:O=.02:Y=.02:Goto1Δ计算直线04 L＜C＝＞I=∫(cos(RX^2,0,L:J=∫(sin(RX^2,0,L:P=RL^2:O=H(D-.02:Y=-O:W=L:L＜2*.02/(F+.02＝＞O=.02-.5D(1-H:Y=.02-.5D(1+H:Goto1Δ△计算缓和曲线05 L＜S-C＝＞P=RC(2L-C):I=BsinP+C/2-C^3/240B^2:J=B(1-cosP) + C^2/24B:W=C:O=HF:Y=-OΔ计算圆曲线06 Lbl 1:I=VI:J=HJ:R=(I+J i)(cosG+i sinG) 坐标转换07 I=E-VT(cosG+i sinG):I=R+I:P=G+HVP 方位角和坐标计算08 V=5+WQ(1+H)/2C:U=5+WQ(1-H)/2C 计算路基宽(假定半幅路基宽为5米）09 Prog "DATE 2":S=Z-A+T:S＞T=＞S=2T-S:D=-E△10 B=B-DT+DS11 S≥0=＞B=B+FS^2/2R△计算中桩高程。
12 "ZHONG X Y P H":I◢ P->DMS◢B◢"YOU X Y H":M=I+V(cos(P+90)+i sin(P+90◢ O=B-O V◢"ZUO X Y H":N=I+U(cos(P-90)+i sin(P-90◢ Y=B-YU◢数据库程序（DATE 1)(平曲线数据）Z<缓直点=>A=xxx:B=xxx:C=xxx:D=xxx:E=xxx:G=xxx:F=xxx:Q=xxx:U=xxx:Goto 2...Lbl 1:U=CLbl 2: H=A/AbsA：A=AbsAT=0.5C-C^3/240B2+(B+U^2/24B-(B+C^2/24BcosA)/sinAO=0.5U-U^3/240B2+(B+C^2/24B-(B+U^2/24BcosA)/sinAS=πAB/180+0.5(C+U)数据库程序（DATE 2)（竖曲线数据）Z<（变坡点+切线长）桩号=>A=xxx:B=xxx:R=xxx:E=xxx:Goto 1...Lbl 1E-D＜0=＞F=-1：Goto2△E-D＞0=＞F=1△Lbl 2:T=Abs(R(E-D)/2说明：使用范围，加宽为渐变加宽，超高加宽在全缓和曲线范围内，超高绕中线旋转。
道路中桩边桩计算程序(附件)
附件：道路中桩、边桩计算公式(fx-5800P程序)主程序：“QXZBJS”“QZH”?B: （线路起点桩号，前一个曲线的HZ或YZ,或是）“LS1”?C:“LS2”?I: （第一缓和曲线与第二缓和曲线，可以不等）”JDZH”?N:”JDX”?G:”JDY”?H: (交点桩号与坐标)”XZH”?M: （下一个曲线的ZH或ZY桩号）“T1”?S:”T2”?J: （第一、二切线长）“QXC”?Z: （曲线全长，含缓和曲线长）“JDPJ”?A: （本交点的转向角）“QDFWJ”?O: （起始点到交点的方位角）“R”?R: （本曲线的圆曲线半径）“W”?W: （曲线偏转信息，左为－１，右为＋１）Lbl 0:“K”?K: （所求断面的桩号）If K<B: （比较所求桩号是否小于起点桩号）Then Goto6: （条件为真，转到Lbl 6）IfEnd:If K>N-S:（判定所求桩号是否大于ZH或ZY点，即所求点是否在曲线段）Then Goto 1:（条件为真，转到Lbl 1）IfEnd:(条件为假时运行下例程序，即所求桩号在第一段直线上)G+(N-K)COS(O+180)→X ▲（以上条件都为假是运行该段，所求桩号在直线段，H+(N-K)Sin(O+180)→Y ▲求其中桩坐标）O→T: （起始方位角赋值与T）Prog “ZI-1”Goto 0Lbl 1:If K>N-S+C: （判定所求桩号是否大于HY点）Then Goto 2: (条件为真，转到Lbl 2，即为所求桩在圆曲线或第二缓和曲线上) IfEnd:(条件为假时运行下例程序，即所求桩号在第一缓和曲线上)((K-N+S)2/(6RC))*(180/π) →Q ▲（所求桩中心所占缓和曲线长度的角度）O+3WQ→T:（ZH到所求点的方位角）√(((K-N+S)-(K-N+S)5/(40R2C2))2+(K-N+S)3/(6RC))2)→D:（所求点与ZH点旋长）G+Scos(180+O)+Dcos(O+WQ)→X ▲（根据ZH点坐标和旋长计算中心点坐标）H+Ssin(180+O)+Dsin(O+WQ)→Y ▲Prog “ZI-1”Goto 0Lbl 2:If K>N-S+Z-I:（判定所求桩号是否大于YH点）Then Goto 3: (条件为真，转到Lbl 3，即为所求桩在第二缓和曲线上)IfEnd:(条件为假时运行下例程序，即所求桩号在圆曲线段上)(90(K-N+S-C))/(πR) →Q:（所求占的圆心角的一半，旋的外角）2Rsin(Q) →D:（所求点与HY点之间的旋长）O+W((C/(2R))*(180/π))+2WQ →T:（HY到所求点的方位角）G+Scos(180+O)+ √((C-C5/(40R2C2))2+(C3/(6RC))2)cos(O+W((C2/(6RC))*(180/π)))+Dcos(O+WQ+W((C/(2R))*(180/π)))→X ▲（根据HY点坐标和旋长计算中心点坐标）H+Ssin(180+O)+ √((C-C5/(40R2C2))2+(C3/(6RC))2)sin(O+W((C2/(6RC))*(180/π)))+Dsin(O+WQ+W((C/(2R))*(180/π)))→Y ▲Prog “ZI-1”Goto 0:Lbl 3:If K>N-S+Z:（判定所求桩号是否大于HZ或YZ点）Then Goto 4: (条件为真，转到Lbl 4，即为所求桩在第二段直线上)TfEnd:(条件为假时运行下例程序，即所求桩号在第二缓和曲线段上)((N-S+Z-K)2/(6RI))*(180/π) →Q（所求桩中心所占缓和曲线长度的角度）O+WA-180-3WQ→T（HZ到所求点的方位角）√((N-S+Z-K-( N-S+Z-K)5/(40R2I2))2+( N-S+Z-K)3/(6RI))2) →D（所求点与HZ 点旋长）G+Jcos(O+WA)+Dcos(OWA-180-WQ) →X ▲（根据HZ点坐标和旋长计算中心点坐标）H+Jsin(O+WA)+Dsin(OWA-180-WQ) →Y ▲Prog “ZI-2”:Goto 0:Lbl 4:If K>M:（判定所求桩号是否大于本次计算的桩号范围）Then Goto 6: (条件为真，转到Lbl 6，即为所求桩超出范围)IfEnd:(条件为假时运行下例程序，即所求桩号在第二段直线上)G+Jcos(O+WA)+(K-(N-S+Z))cos(O+WA) →X ▲（中心坐标）H+Jsin(O+WA)+(K-(N-S+Z))sin(O+WA) →Y ▲O+WA→T:Prog “ZI-1”Goto 0:Lbl 6:“END”：子程序“ZI-1”(求边桩坐标)“L”?L:（输入中心至左右边桩的距离L，左为负，右为正）X+Lcos(T+90) →E ▲Y+Lsin(T+90) →F ▲“ZI-2”(求边桩坐标)“L”?L:（输入中心至左右边桩的距离L，左为负，右为正）X+Lcos(T-90) →E ▲Y+Lsin(T-90) →F ▲。
中边桩高程及加宽计算程序
中边桩高程及加宽计算程序（绕中轴旋转高次抛物线超高加宽）2008年4月23日徐锋GCJKJS（主程序）A“P1”B“P2”RD“JD ZHANG”I“JD-H”S“LH”G“LJK”M“SJJK”P“LGHP”Q“CGHP”U“ZH”X“HZ”Z“(Z－，Y+)”：Z[1]=Z：LBI0：{Z}：Z：W=B-A：T=AbsWR÷2：E=TW÷4：C=D-T：F=D+T：W≥0=>K=（Z-C）2÷（2R）：≠> K= -（Z-C）2÷2R：Goto1 Lbi1：Z≤D-T=>Goto5：≠>Z＞D-T=>Z≤D+T=>Goto3：≠>Z＞D+T=>Goto4：Lbi5：H=I-A （D-Z）GotoD：Lbi3：H=I-A（D-Z）+KGotoD：Lbi4：H=I-B（D-Z）GotoDLbiDZ≤U=>L＝0：Goto2⊿Z≤U+S=>L＝Z－U：Goto2⊿Z≤X－S=>L＝S：Goto2⊿Z≤X=>L ＝X－Z：Goto2⊿Z>X=>L＝0：Goto2⊿Lbi2：J“JK”=（4×（L÷S）3-3×（L÷S）4）×MY“WC”=（L÷S×（P+Q））×G÷2-P×G÷2O=2P÷（P+Q）×S：L≤O=>V“LC”=-P×（G÷2+J）GotoC：≠=>V“LC”=-（（L÷S）×（P+Q）-P）×（G÷2+J）GotoCLbiCZ[1]>0=> Prog“A”：Goto 0⊿Prog “B”：Goto0⊿A（子程序）“G=”：H=H◢L<0.001=>Y=-Y: “ZG=”：H=H+Y◢“YG=”：H=H◢Goto 1 ≠=>Y=Y:“ZG=”：H=H+Y◢“YG=”：H=H-Y+V◢Goto 1Lbi1“ZK=”：Z[2]=G/2◢“YK=”：Z[3]=G/2+J◢B（子程序）“G=”：H=H◢L<0.001=>V=V: “ZG=”：H=H+V◢“YG=”：H=H◢Goto 1≠=>Y=Y:“ZG=”：H=H+V◢“YG=”：H=H+Y－V◢Goto 1Lbi1“ZK=”：Z[3]=G/2+J◢“YK=”：Z[4]=G/2◢说明：一、程序中字母及符号意义：P1（A）：第一条纵坡坡度（i %，如0.02、-0.02）P2（B）：第二条纵坡坡度（i %，如0.02、-0.02）R：竖曲线半径JD ZHANG（D）？变坡点桩号I（JD-H）？变坡点高程S“LH”？缓和曲线长（当没有缓和曲线时，请输为0.0001）G“LJK”？标准路基全幅宽度M“SJJK”？弯道设计最大加宽值P“LGHP”？标准路拱横坡坡值Q“CGHP”？弯道设计最大超高值U“ZH”？直缓桩桩号X“HZ”？缓直桩桩号Z“(Z－，Y+)”？左偏输入负值，右偏输为正值Z？待求点桩号“G=”中桩高程“ZG=”左边桩高程“YG=”右边桩高程“ZK=”左幅路基宽度“YK=”右幅路基宽度二、输入、计算要点1．一次只能输入一条平竖曲线参数2．纵坡坡度须代入正负号进行运算特别说明：1、本程序单下划线部份是中桩高程计算部份2、本程序双下划线部份是超高加宽计算部份3、高程计算部份引用了杨小杰的程序，在此特别感谢优点：将高程、超高、加宽一次性计算出来，减少计算的工作量中边桩高程及加宽计算程序（绕中轴旋转高次抛物线超高加宽）2008年4月23日徐锋GCJKJS（主程序）LbiP：A“P1”B“P2”RD“JD ZHANG”I“JD-H” S“LH” G“LJK” M“SJJK” P“LGHP”Q“CGHP” U“ZH” X“HZ” Z“(Z－，Y+)”：Z[1]=Z：Goto0回车LBI0：{Z}：Z：W=B-A：T=AbsWR÷2：E=TW÷4：C=D-T：F=D+T：W≥0=>K=（Z-C）2÷（2R）：≠> K= -（Z-C）2÷2R：Goto1回车Lbi1：Z≤D-T=>Goto5：≠>Z＞D-T=>Z≤D+T=>Goto3：≠>Z＞D+T=>Goto4：Lbi5：H=I-A（D-Z）回车GotoD：Lbi3：H=I-A（D-Z）+K回车GotoD：Lbi4：H=I-B（D-Z）回车GotoD回车LbiD回车Z≤U=>L＝0：Goto2⊿Z≤U+S=>L＝Z－U：Goto2⊿Z≤X－S=>L＝S：Goto2⊿Z≤X=>L＝X－Z：Goto2⊿Z>X=>L＝0：Goto2⊿回车Lbi2：J“JK”=（4×（L÷S）3-3×（L÷S）4）×M回车Y“WC”=（L÷S×（P+Q））×G÷2-P×G÷2回车O=2P÷（P+Q）×S：L≤O=>V“LC”=-P×（G÷2+J）回车GotoC：≠=>V“LC”=-（（L÷S）×（P+Q）-P）×（G÷2+J）回车GotoC回车LbiC回车Z[1]>0=> Prog“A”：Goto 0⊿Prog “B”：Goto0⊿A（子程序）“G=” ：H=H◢L<0.001=>Y=Y: “ZG=” ：H=H+Y◢“YG=” ：H=H◢Goto 1 ≠=>Y=Y:“ZG=” ：H=H+Y◢“YG=” ：H=H-Y+V◢Goto 1Lbi1“ZK=” ：Z[2]=G/2◢“YK=” ：Z[3]=G/2+J◢B（子程序）“G=” ：H=H◢L<0.001=>V=V: “ZG=” ：H=H+V◢“YG=” ：H=H◢Goto 1≠=>Y=Y:“ZG=” ：H=H+V◢“YG=” ：H=H+Y－V◢Goto 1Lbi1“ZK=” ：Z[3]=G/2+J◢“YK=” ：Z[4]=G/2◢说明：一、程序中字母及符号意义：P1（A）：第一条纵坡坡度（i %，如0.02、-0.02）P2（B）：第二条纵坡坡度（i %，如0.02、-0.02）R：竖曲线半径JD ZHANG（D）？变坡点桩号I（JD-H）？变坡点高程S“LH”？缓和曲线长（当没有缓和曲线时，请输为0.0001）G“LJK”？标准路基全幅宽度M“SJJK”？弯道设计最大加宽值P“LGHP”？标准路拱横坡坡值Q“CGHP” ？弯道设计最大超高值U“ZH”？直缓桩桩号X“HZ”？缓直桩桩号Z“(Z－，Y+)”？左偏输入负值，右偏输为正值Z？待求点桩号“G=”中桩高程“ZG=”左边桩高程“YG=”右边桩高程“ZK=”左幅路基宽度“YK=”右幅路基宽度二、输入、计算要点1．一次只能输入一条平竖曲线参数2．纵坡坡度须代入正负号进行运算特别说明：1、本程序单下划线部份是中桩高程计算部份2、本程序双下划线部份是超高加宽计算部份3、高程计算部份引用了杨小杰的程序，在此特别感谢优点：将高程、超高、加宽一次性计算出来，减少计算的工作量本文可共享补充点击修改共享补充New!查看楼主其它文章<<返回→回复内容只查看楼主与本人回复沙发：xufeng 楼主08-4-23 16:37回欢迎大家使用，并提出宝贵的意见，以便进一步提高，相信对你的测量计算工作有所帮助2楼：李洪08-4-24 18:19回能算匝道不?3楼：'geergirl 08-4-25 19:45回评：-4分代码错误很多，根本无法正常运行！4楼：XUFENG 楼主08-4-26 11:53回3楼:geergirl3楼的朋友切莫乱加评论，你用了吗，没用过就别乱糟蹋别人的心血，我的程序给你这样的人看实在很不幸5楼：XUFENG 楼主08-4-26 11:55回2楼:李洪只算主线，匝道一般用内插计算6楼：XUFENG 楼主08-4-26 12:8回程序的好坏，要用了才知道，如果真有问题请指出不对之处，大家共同进步，不然就不要妖言祸众，让别人不敢用7楼：'123 08-4-26 12:25回谢谢，用后评。
4800公路中边桩高程计算程序
4800公路中边桩高程计算程序４８００公路中边桩高程计算程序fx-4800高等级公路中边桩标高程计算程序（适用超高方式为有分隔带绕左内侧旋转，无分隔绕线路中心旋转）随着高等级公路建设的不断发展，对路基施工高程控制要求不断提高，经常在现场要对路基中边桩高程进行测量，由于施工线路较长，高程数据较多，而且只提供设计桩号的高程，对非设计断面的横坡及高程计算相当复杂，因此利用fx-4800计算器的编程功能及扩展变量的储存功能，编制一个程序，可以预先输入整条线路的曲线要素，使用时，只要输入任意点桩号就可以计算中边桩高程。
程序由中桩高程计算主程序和横坡计算子程序组成，程序如下：1、中桩高程计算程序文件名（ZZGC）程序说明Lb1 0Fix 3（输入变坡点数） N=?（输入扩展变量位置数） V=?M=V（输入待求桩号） {K}Lb1 1K>Z[2+M]-Z[3N+3+M]=> （确定待求桩号在哪个竖曲线内）M=M+1：Goto 1（读取该竖曲线半径） R=Z[2N+2+M]（读取该竖曲线切线长） T=Z[3N+2+M]（读取该变坡点前坡率）I=Z[4N+2+M]÷100（读取该变坡点后坡率）L=Z[4N+3+M]÷100（计算待求桩号与该竖曲线起点距离）S=K-Z[1+M]+Z[3N+2+M]M=V=> （如果桩号在起点与第1竖曲线起点之间，转Lb1 5 ）Goto 5⊿（计算坡率差） A=I-L（判断竖曲线凹凸）J=AbsA÷AS （如果距离在前半曲线内，计算改正值）Y=S^2÷2R ：Goto 2 ⊿S<2T=> （如果距离在后半曲线内，计算改正值）Y=(2T-S)^2÷2R ：Goto 3 ⊿Goto4Lb1 2H=Z[N+2+M]-I(T-S)- YJD"H" =D ◢Goto 6Lb1 3H= Z[N+2+M]+L(S-T)- YJD"H" =D ◢Goto 6Lb1 4（在直线内，计算高程） H=Z[N+2+M]+L(K-Z[1+M])D"H" =D ◢Goto 6Lb1 5I=Z[V+4N+3] ÷100H=Z[V+N+3]-I(Z[V+2]-K）D"H" =D ◢Lb1 6(输入边桩距离，右正，左负) {B}Prog "HP js" (转横坡计算子程序，应与横坡计算文件名相同)Fix 3B<0=> D=D-QAbsB÷100 ：Goto 7 ⊿B>0=> D=D-PAbsB÷100 ：Goto 8 ⊿Lb1 7D"ZBH" =D ◢（显示左边桩高程）Goto 0Lb1 8D"YBH" =D ◢（显示右边桩高程）Norm：1Goto 0中桩高程计算程序扩展变量设置表起点第1变坡点第2变坡点第3变坡点第n变坡点终点里程Z[V+1]Z[V+2]Z[V+3]Z[V+4]Z[V+n+1] Z[V+n+2]高程Z [V+n+3]Z[V+n+4]Z[V+n+5]Z[V+2n+2]R 竖曲线半径Z[V+2n+3] Z[V+2n+4] Z[V+2n+5] Z[V+3n+2] T（竖曲线切线长）Z[V+3n+3] Z[V+3n+4] Z[V+3n+5] Z[V+4n+2]I （坡率）Z[V+4n+3] Z[V+4n+4] Z[V+4n+5] Z[V+5n+2] Z[V+5n+3]2、横坡计算程序文件名（HP js）程序说明Lb1 0（输入平曲线交点数） N=?（输入标准横坡值） I=?（输入LC值） T=?（输入扩展变量位置） V=?M=V(输入桩号,如作为子程序使用则无此行) {K}Lb1 1K>Z[2+M] => M=M+1 ：Goto 1 ⊿（读取平曲线全长） L=Z[N+2+M]（读取LS1） E=Z[2N+2+M]（读取LS2） F=Z[3N+2+M]（读取标准超高值） H=Z[4N+2+M]（读取线路偏向，右=1，左=-1） A=Z[5N+2+M]S=K-Z[M+1]M=V=> P=I ：Q=I ：Goto 7 ⊿E=0=> （无LS则不设超高） P=I ：Q=I ：Goto 7 ⊿S≤T=> （如果该点在第1 Lc内，转Lb2） S=S ：Goto 2 ⊿S≤E=> （如果该点在 Lc-Ls1终点，转Lb3） S=S-T ：Goto 3 ⊿S≤L-F=> （如果该点在圆曲线内，转Lb4）Goto 4 ⊿S≤L-T=> （如果该点在Ls2起点-Lc，转Lb 5）S=S-(L-F) ：Goto 5 ⊿S≤L=> （如果该点在第2Lc内，转Lb6） S=T-(L-S) ：Goto 6 ⊿（在直线） P=I ：Q=I ： Goto 7Lb1 2A=1=> （右偏） P=I ：Q=I-2IS÷T ：Goto 7 ⊿A=-1=> （左偏） P= I-2IS÷T ：Q=I ：Goto 7 ⊿Lb1 3A=1=> P=I+S(H-I) ÷(E-T) ：Q=-P ：Goto 7 ⊿A=-1=> P=-I- S(H-I) ÷(E-T) ：Q=-P ：Goto 7 ⊿Lb1 4A=1=> P=H ：Q=-H ：Goto 7 ⊿A=-1=> P=-H ：Q=H ：Goto 7 ⊿Lb1 5A=1=> P=H-S(H-I) ÷(F-T)：Q=-P ：Goto 7 ⊿A=-1=> P=-H+S(H-I) ÷(F-T)：Q=-P ：Goto 7 ⊿Lb1 6A=1=> P=I：Q=-I+2IS÷T ：Goto 7 ⊿A=-1=> P=-I+2IS÷T ：Q=I ：Goto 7 ⊿Lb1 7Fix 3P"YI"=P ◢（显示右横坡，如不需显示则无“◢”）Q"ZI"=Q ◢（显示左横坡，如不需显示则无“◢”）Norm ：1(如作为子程序使用则无此行) Goto 0横坡计算程序扩展变量设置表起点第1交点第2交点第3交点第Z交点终点线路起点、曲线起点、线路终点桩号Z[V+1] Z[V+2] Z[V+3] Z[V+4] Z[V+Z+1] Z[V+Z+2]L (平曲线全长) Z[V+Z+3] Z[V+Z+4] Z[V+Z+5] Z[V+2Z+2]LS1 Z[V+2Z+3] Z[V+2Z+4] Z[V+2Z+5] Z[V+3Z+2]LS2 Z[V+3Z+3] Z[V+3Z+4] Z[V+3Z+5] Z[V+4Z+2]i (标准超高值) Z[V+4Z+3] Z[V+4Z+4] Z[V+4Z+5] Z[V+5Z+2]偏转方向(右偏=1，左偏=-1) Z[V+5Z+3] Z[V+5Z+4] Z[V+5Z+5] Z[V+6Z+2]3、使用说明把上述程序输入到fx-4800内，然后扩展足够的变量个数，即至少为5N+3+6Z+2个，N为变坡点个数，Z 为平曲线交点数。
公路逐桩坐标计算程序
公路逐桩坐标计算程序(可以计算对称、不对称缓和曲线)Lb1 0Z=？V=？W=V+2：Fixm｛K｝Lb1 1K>Z[W+5Z+4]=>W=W+1:Goto 1⊿(判断桩号在哪个交点范围，就是该交点曲线起点至下一交点曲线起点) S=K-Z[W+5Z+3] （计算该桩号与曲线起点的距离）R=Z[W+2Z+2]：L=Z[W+3Z+2]:E=Z[W+4Z+2] （读取该交点曲线要素R、Ls1 、Ls2）Pol（Z[W]-Z[W-1],Z[W+Z+2]-Z[W+Z+1]）（计算该交点与下一交点直线方位角）J<0=>J=J+360⊿A=JPol（Z[W-1]-Z[W-2],Z[W+Z+1]-Z[W+Z]）（计算该交点与上一交点直线方位角）J<0=>J=J+360⊿C=A-J:A=J （计算偏角）W=V+2=>Goto2⊿（如果桩号在起点与第一交点曲线起点之间，则转Lb1 2 ）I=Abs（tan（c÷2））M=L÷2-L^3÷240R^2:N=E÷2-E^3÷240R^2P=L^2÷6R-L^4÷336R^3-R（1-cos（90L÷πR））Q=E^2÷6R-E^4÷336R^3-R（1-cos（90E÷πR））D=（P-Q）I÷2 : F=（P+Q+2R）I÷2M=F+M-D:Q=F+N+DN=πRAbsC÷180+（L+E）÷2X=Z[W-1]-McosAY=Z[W+Z+1]-MsinAM=Z[W-1]+Qcos（A+C）V=Z[W+Z+1]+Qsin（A+C）Q=AbsC÷CS≤L=>P=0:Goto3⊿（如果桩号在第一缓和曲线内，则转Lb1 3）S≤N-E=>S=S-L:Goto4⊿（如果桩号在圆曲线内，则转Lb1 4）S≤N=>S=N-SQ=-Q:A=A+C-180:X=M:Y=V:L=E：P=180:Goto3 ⊿（如果桩号在第二缓和曲线内，则转Lb1 3）P=A+C:S=S-N:D=M+ScosP:F=V+SsinPGoto6 （如果桩号在直线内，则转Lb1 6）Lb1 2P=A+CD=Z[W-1]+ScosPF=Z[W+Z+1]+SsinP:Goto6Lb1 3I=S-S^5÷40R^2÷L^2+S^9÷3456R^4÷L^4J=Q（S^3÷6RL-S^7÷336R^3÷L^3）P=P+A+90Q S^2÷πRL:Goto5Lb1 4M=90（2S+L）÷πRI=RsinM+L÷2-L^3÷240R^2J=Q（L^2÷24R+R（1-cosM））P=A+QMLb1 5D=X+IcosA-jsinA:F=Y+JcosA+IsinALb1 6D″X=″◢（结果显示X坐标）F″Y=″◢（结果显示Y坐标）P″AT=″◢（结果显示该桩号方位角）{BO}：B″S″O″⊿″ （输入边桩距离，交角）P=P+OL″XB″=D+BcosP◢（结果显示边桩X坐标）M″YB″=F+Bs inP◢（结果显示边桩Y坐标）以上是坐标计算程序，括号内是程序计算的大致原理及说明，中间部分为直线、圆曲线、缓和曲线计算的各种公式，大家也知道，书上也有。
道路中边桩坐标计算
内移距和切线增长距则可取:
(２－6)
2、局部坐标计算
（1）、缓和曲线段。缓和曲线段上各待定点坐标按缓和曲线参数方程计算,即
（2-7）
（2）、圆曲线线段。圆曲线段上各待定点坐标,可按图2.2写出
（２-８)
图2.2圆曲线局部坐标
注：式中 为圆曲线上的点到圆曲线起点的弧长(里程差）
2.3 中桩坐标计算
平面曲线按其半径的不同分为圆曲线和缓和曲线。圆曲线上任意一点的曲率半径处处相等。缓和曲线是在直线与圆曲线，圆曲线与圆曲线之前设置的曲率半径连续渐变的一段过渡曲线;缓和曲线上任意一点曲率半径处处在变化。当缓和曲线作为直线与圆曲线之间的介曲线时，其半径变化范围自无穷大至圆曲线半径R，若用以连接半径为Ｒ1和Ｒ2的圆曲线时,缓和曲线的半径便自Ｒ1向R２过渡。
２、切线角公式,如图2.1所示,可知切线角公式为：
（2-2）
3、回旋线参数方程式为:
(2-3）
注:当圆曲线半径较大时,一般略去高次项,x只取前一、二项，y取前一项即可。缓和曲线终点HY（或YH)的坐标即为:
（2-4)
2.2.2缓和曲线局部坐标计算
1、如图2．1当半径较小时应取更多的项,实际计算取前五项即可,其中Ａ为回旋线参数，以下为回旋线参数方程取前五项的计算公式:
(2-1５）
故AP方向的坐标方位角和Ｐ点坐标为：
(２-16）
3、采用坐标转换计算中桩坐标，利用坐标转换公式将xｏy坐标系下的坐标的切线角,曲线右偏为正,左偏为负，再将 坐标系下的坐标转换到测量坐标系下则可得：
(2-18）
其中 为A点的切线在测量坐标系下的坐标方位角。
按曲线的连接方式不同,可分为：
a、单圆曲线，亦称为单曲线，即具有单一半径的曲线