该程式主要适用于对称复曲线（直线、缓曲、圆曲、缓曲、直线）中边桩坐标计算，单次半算范围为前一交点HZ点到后一交点ZH点。
E“JDK” F“JDX” G“JDY” R H“LS” I“ALF”J“ZJ”
Ｔ＝（Ｒ＋Ｈ２÷２４Ｒ－Ｈ3÷２６８８Ｒ3）Ｔａｎ（ＡｂｓＪ÷２）＋Ｈ÷２－Ｈ3÷２４０Ｒ２
Ｌ＝Ｒ（πＡｂｓＪ÷１８０－Ｈ÷Ｒ）
Ｌｂｏ０：｛ＫＱ｝：Q=90:X＝F－TCosM：Y=G＋TSinM：Z＝Ｅ－Ｔ：Z＝K－O：K＞Z →ＧOTO１：Z=K-Z: M＝I： P＝M－Q：GOTO5
Ｌbo１：Ｊ＜０→Ｎ＝１：←Ｎ＝２△Ｋ＞Z＋Ｈ→ＧOto２：p＝（－１）Ｎ（Ｋ－Z）２÷Ｈ÷Ｒ÷６×１８０÷π：Ｍ＝Ｉ＋p：Ｚ＝Ｋ－Z－（Ｋ－Z）５÷９０÷Ｒ２Ｈ２：Ｐ＝Ｍ＋P－Ｑ：Ｇto５
Ｌbo２：z＝z＋Ｈ：Ｋ＞z＋Ｌ→Ｇoto３：o＝Ｉ＋（－１）ＮＨ÷２÷Ｒ×１８０÷π：p＝（－１）Ｎ（Ｋ－Z）÷Ｒ×１８０÷π：Ｍ＝Ｉ＋（－１）ＮＨ÷Ｒ÷６×１８０／π：Ｚ＝Ｈ－Ｈ３÷９０Ｒ２：Ｘ＝Ｘ＋ＸＣosＭ：Ｙ＝Ｙ＋ＳinＭ：Ｍ＝o＋p÷２：Ｚ＝２ＲＳin（（－１）ＮP÷２）：q＝o＋p－Ｑ：ＧotoＯ５
Ｌbo３：z＝zＨ＋Ｌ：Ｍ＝Ｉ＋Ｊ：Ｘ＝Ｆ＋ＴＣosＭ：Ｙ＝Ｇ＋ＴＳinＭ：ＫＬbo４：Ｊ<０→Ｎ＝２：←Ｎ＝２△p＝（－１）Ｎ（Ｋ－Z）２÷Ｈ÷Ｒ÷６×１８０÷π：Ｚ＝Ｋ－Z－（Ｋ－Z）５÷９０÷Ｒ２Ｈ２：Ｍ＝Ｍ＋p：q＝Ｍ＋２p－Ｑ：Ｇoto５
Ｌbo５：｛V｝：X＝Ｘ＋ＺＣosＭ＋V“WIDE”ＣosQ▲Y＝Y＋ＺSINＭ＋V“WIDE”SINQ▲
符号说明：
ＪＤK？：交点桩号；
ＪＤＸ？：交点桩Ｘ坐标；
ＪＤＹ？：交点Ｙ坐标；
Ｒ？：圆曲线半径；
ＬＳ？：缓和曲线长；
ＡＬＦ？：该交点与前一交点直线段的方位角（度）；
ＺＪ？：该曲线转角（左转为正、右转为负）；
Ｋ？：欲求点桩号；
WIDE？：中桩至边桩距离（左正右负，中桩为零）
X：欲求点X坐标；
Y欲求点Y坐标；
1、计算范围为前一交点HZ至后一交点ZH。
2、计算曲线为单圆曲线时，LS输入一无限小数值（如LS=1E45）。
3、程式中为使计算时输入元素最少， 边桩设为正交边桩，需计算任意交角边桩时，需把程式中第三行 Q=90删除即可在计算过程中提示要求输入边桩角度。
本人为半路出家的和尚，工民建专业毕业后从事公路施工测量。那时亦未能发现网站上可以下载相关计算程式。闭门造车，猛查资料三天方至成型。后一直没空未对其加以深究，现休息在家，将其贴于该站，一者请大家多多指教，二者解部分朋友当日我之困境。
FDS-1:我有一个输入切线长的程序,你可试用
Ｇ“Ｊ（Ｘ）”：Ｈ“Ｊ（Ｙ）”：Ｍ“Ｊ（Ｋ）”：Ｒ：Ｌ“ＬＳ１”：Ｗ“ＬＳ２”：Ｔ“Ｔ１”：Ｕ“Ｔ２”：Ｚ“Ｚ（Ｋ）”：Ｃ“Ｑ（ＦＷＪ）”：Ｄ“Ｚ（ＦＷＪ）”：Ｊ＝Ｃ＋１８０－Ｄ：Ｊ＜０＝＞Ｊ＝Ｊ＋３６０△Ｊ＞３６０＝＞Ｊ＝Ｊ－３６０△Ｊ＜１８０＝＞Ｐ＝１：≠＞Ｐ＝－１△Ｑ＝Ｍ－Ｔ：Ａ＝Ｇ＋Ｔcos(C+180)：Ｂ＝Ｈ＋Ｔsin(C+180)：Ｅ＝Ｇ＋ＵcosD：Ｆ＝Ｈ＋ＵsinD
Lbl 0：Ｆixm：{Ｋ｝：Ｋ≤Ｑ＝＞“ＦＷＪ＝”：Ｎ＝Ｃ→ＤＭＳ▲“Ｘ＝”：Ａ＋Rec(Q-K，N+180)▲“Ｙ＝”：Ｂ＋Ｊ▲Goto 0△Ｋ≥Ｚ＝＞“ＦＷＪ＝”：Ｎ＝Ｄ→ＤＭＳ▲“Ｘ＝”：Ｅ＋Rec(K-Z，N)▲“Ｙ＝”：Ｆ＋Ｊ▲Goto 0△Ｚ－Ｋ＜Ｗ＝＞Ｚ[3]＝Ｗ：Ｓ＝Ｚ－Ｋ：Prog“PQ1”：“ＦＷＪ＝”：Ｎ＝Ｄ－Ｚ[7]→ＤＭＳ▲“Ｘ＝”：Ｅ＋Ｉcos(Ｄ+180)＋Ｊsin(Ｄ+180)▲“Ｙ＝”：Ｆ－Ｊcos(Ｄ+180)＋Ｉsin(Ｄ+180)▲Goto 0△Ｚ[3]＝Ｌ：Ｋ－Ｑ＜Ｌ＝＞Ｓ＝Ｋ－Ｑ：Prog“PQ1”：“ＦＷＪ＝”：Ｎ＝Ｃ＋Ｚ[7]→ＤＭＳ▲“Ｘ＝”：Ａ＋ＩcosＣ－ＪsinＣ▲“Ｙ＝”：Ｂ＋ＪcosＣ＋ＩsinＣ▲Goto 0△Ｓ＝Ｋ－Ｑ－Ｌ：Ｚ[４]＝Ｃ＋90P(L+S)/(Rπ)
Z[5]=2Rsin(S÷R÷π×90):N=Z[4]+90PS÷R÷π:Z[6]=N:S=L:L=0＝＞I=0：J=0≠＞Prog“PQ1”：“FWJ=”：N=Z[6]→DMS▲X=A+IcosC－JsinC：Y=B+JcosC+IsinC：“X=”：X+Rec（Z[5]，Z[4]）▲“Y=”：Y+J▲Goto0△
PQ1：
I=S－S^5÷（40R^2Z[3]^2）+S^9/(3456R^4Z[3]^4)：J=P（S^3/(6RZ[3])－S^7/(336R^3Z[3]^3)）：Z[7]=90PS^2÷（Z[3]πR）
缓和段坐标计算（ZH---HY）
Z:Z←┘
Lbl1:{K}:K←┘
L=k-261301.35（ZH里程）:C=L-(L^5÷40÷7000^2÷200^2):B=(L^3÷6÷7000÷200)-(L^7÷336÷7000^3÷200^3)←┘
"X=":X=2571749.599(ZHX坐标)+CCos221'41'30.71"(ZH前方位角)-ZBSin221'41'30.71"◢
"Y=":Y=517518.910(ZHY坐标)+CSin221'41'30.71"(ZH前方位角)+ZBCon221'41'30.71"◢
{DR}:R"RA":D←┘
"N=":N=X+DCos(221'41'30.71"+R+ZL^2÷2÷200÷7000*180÷π)◢
"E=":E=Y+DSin(221'41'30.71"+R+ZL^2÷2÷200÷7000*180÷π)◢
Goto1
7000---为半径 不同的曲线修改下就OK 比如下一个半径为600，把7000该为600就OK
200---为缓和长 不同的缓和长修改下就OK 比如下一个半径为150，把200该为150就OK
（261301.35为ZH里程） (ZHX坐标) (ZHY坐标) (ZH前方位角) 都是根据你需要的要素而改动，我的着个上面只是个参数例子。不等长就需要修改
圆曲线坐标计算（HY---YH）
Z:Z←┘
Lbl1:{K}:K←┘
M=(200*180)÷2÷7000÷π:F=(((K-261501.35)*180)÷7000÷π)+M:B=(7000SinF)+(200÷2-200^3÷240÷7000^2):C=7000(1-CosF)+(200^2÷24÷7000)←┘
"X=":X=2571749.599(ZHX坐标)+BCos221'41'30.71"(ZH前方位角)-ZCSin221'41'30.71"◢
"Y=":Y=517518.910(ZHY坐标)+BSin221'41'30.71"(ZH前方位角)+ZCCon221'41'30.71"◢
{DR}:R"RA":D←┘
"N=":N=X+DCos(221'41'30.71"+ZM+R+Z(K-261501.35)÷7000*180÷π)◢
"E=":E=Y+DSin(221'41'30.71"+ZM+R+Z(K-261501.35)÷7000*180÷π)◢
Goto1
7000---为半径 不同的曲线修改下就OK 比如下一个半径为600，把7000该为600就OK
200---为缓和长 不同的缓和长修改下就OK 比如下一个半径为150，把200该为150就OK
（261501.35为HY里程） (ZHX坐标) (ZHY坐标) (ZH前方位角) 都是根据你需要的要素而改动，我的着个上面只是个参数例子。
后缓和段坐标计算（YH---HZ）
Z:Z←┘
Lbl1:{K}:K←┘
L=263047.451（HZ里程）-K:C=L-(L^5÷40÷7000^2÷200^2):B=(L^3÷6÷7000÷200)-(L^7÷336÷7000^3÷200^3)←┘
"X=":X=2570329.184(HZ坐标)-CCos（221'41'30.71"(ZH前方位角)+Z*12'39'18"(线路转角)）+Z(-B)Sin（221'41'30.71"(ZH前方位角)+Z*12'39'18"(线路转角)）◢
"Y=":Y=516510.803(HZ坐标)-CSin（221'41'30.71"(ZH前方位角)+Z*12'39'18"(线路转角)）-Z(-B)Con（221'41'30.71"(ZH前方位角)+Z*12'39'18"(线路转角)）◢
{DR}:R"RA":D←┘
"N=":N=X+DCos（221'41'30.71"(ZH前方位角)+Z*12'39'18"(线路转角)）+R-ZL^2÷2÷200÷7000*180÷π)◢
"E=":E=Y+DSin（221'41'30.71"(ZH前方位角)+Z*12'39'18"(线路转角)）+R-ZL^2÷2÷200÷7000*180÷π)◢
Goto1
7000---为半径 不同的曲线修改下就OK 比如下一个半径为600，把7000该为600就OK
200---为缓和长 不同的缓和长修改下就OK 比如下一个半径为150，把200该为150就OK
（263047.451为HZ里程） (HZ坐标) (HZ坐标) (ZH前方位角) (线路转角) 都是根据你需要的要素而改动，我的着个上面只是个参数例子。不等长就需要修改 缓和长 200 要有不等长的线路可以根据需要来修改
直线上就不用说的大家都会的哦
本程序说明
程序运行：
Z 为线路左偏或右偏 左偏时输入-1 右偏时输入1 切记
K 为所求的里程
X 为中桩坐标X
Y 为中桩坐标Y
RA 为线路算边桩左侧和右侧角度 任意角度 算左边距离时输入-90或-80度 算右边距离时输入+90或+80度
D 为距离
N 为边桩坐标X
E 为边桩坐标Y
下面是：例子 为线路左偏
缓和段坐标计算（ZH---HY）例子
里程：261501.35 中桩 X=2571599.622 Y=517386.5976 左侧-10度、距离为10 M N=2571591.037 E=517381.4661 右侧+10度、距离为10 M N=2571593.309 E=517378.84
圆曲线坐标计算（HY---YH）例子
里程：261600 中桩 X=2571524.572 Y=517322.5698 左侧-70度、距离为10 M N=2571515.909 E=517327.4613 右侧+70度、距离为10 M N=2571528.006 E=517313.178
后缓和段坐标计算（YH---HZ）例子
里程：262847.452 中桩 X=2570503.579 Y=516608.7076 左侧-90度、距离为10 M N=2570498.602 E=516617.3815 右侧+90度、距离为10 M N=2570508.558 E=516600.0358
作者：非非（feizai1126)
联系QQ：8696847
联系QQ：8696847
联系QQ：8696847