对称复曲线（直线、缓曲、圆曲、缓曲、直线）中边桩坐标计算

该程式主要适用于对称复曲线（直线、缓曲、圆曲、缓曲、直线）中边桩坐标计算，单次半算范围为前一交点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