匝道上回旋线任意点坐标方位角计算
源程序：
“ZDHQ ”
C“A2”: H=1=＞ M=0 ◣ H=2 =＞M=C/R ◣ H=1=＞ G=N“ZH（G）HY”: ≠＞ G=N“ZH（G）HY”+180D/Rл ◣ G <360 G“ZH（G）YH”◢ ≠＞ G“ZH（G）YH”=G-360◢ ◣ Lbl 1：{ZLAB}：Z=0 Goto2 ◣ P=1=＞ L=M+L-K ◣ P=-1= ＞L=M-L+K ◣ A：B： T=G-90PM2/Cл： E=T：E=E+90PL2/Cл：E <360 E◢ ≠＞ E=E-360 ◢ ◣ O=（L∧3-M∧3）/6C-（L∧7-M∧7）/336C∧3+（L∧11-M∧11）/42240C∧5-（L∧15-M∧15）/9676800C∧7：F=L-M-（L∧5-M∧5）/40C2+（L∧9-M∧9）/3456C∧4：U“XN”=W+PF cosT-OsinT ◢ Z“YN”=V+PFsinT+OcosT ◢ X=U： Y=Z： Prog“F”： Prog“BZ”：Goto 1：Lbl 2：Z=0=＞W=O：V=F：P=＞1= K=L-M+K ◣ P=-1=＞K=M-L+K ◣ F=G：“END”
子程序
“XY”
U“XP”=X+ Rec(S，Q ◢ Z“YP”=Y+J
“ F ”
Pol（X-A，Y-B：J＜0=＞ J=J+360 ◢ ≠＞ J ◢ ◣ I ◢
“BZ ”
X=U: Y=Z： O=U: F=Z: {S}: T=1 ◢ Q=E-90 : Prog“XY”◢ X=U: Y=Z:
Prog“F”: J◢:I”S”◢{S}: T=T+1 ◢ X=O: Y=F: Q=E+90：Prog“XY”◢ X=U: Y=Z: Prog“F”:J◢:I”S”◢
说明：
A2 = L0R1R2 /（R1-R2） 当R1为∞时则R1=1000000R2
A2 — 曲率半径变化率的平方
H — H= 1计算第一缓和曲线 H=2 计算第二缓和曲线
ZH（G）HY— H = 1时过ZH H = 2时过HY起算方位角
P — 起点RK半径＞终点RL半径时输入1，反之输入- 1
L — 回旋线待求点里程
K — 回旋线起点里程
E — 所求点切线方位角
A、B — 测站点坐标
W、V — 回旋线起点坐标
S — 中线至边线距离
R — 已算过的前圆曲线半径，（左负右正）
D — 已算过的前圆曲线长度
Z = 0 结束本曲线计算，显示END
求：
ZH（G）HY — 过本曲线ZH或
过YH点切线方位角
E — 所求点切线方位角
XN、YN、J、I — 所求桩号点中桩及左、右边桩坐标 及 测站至所求点方位角、边长
T=1 为所求桩号点左边桩坐标
T=2 为所求桩号点右边桩坐标
算例：
计算第一缓和曲线：
已知： A2= 902 H= 1 ZH = 134.6108889
P= -1 L= 50 K= 47.86
W= 933.866 V= 312.967 S=20
A= 900 B= 300
成果：
ZH（G）YH =134.6108889 过本曲线ZH、YH切线方位角
E = 134.5946919 （所求点切线方位角）
中桩：XP= 932.363 =24°07 ′13.34 ″测站至所求点方位角
YP=314.491 S = 35.459 测站至所求点边长
左桩：XP= 946.605 =31°28 ′32.17 ″测站至所求点方位角
YP= 328.532 S = 54.645 测站至所求点边长
右桩：XP=918.121 = 1°25 ′07.93 ″ 测站至所求点方位角
YP =300.449 S = 18.127 测站至所求点边长
计算第二缓和曲线：
A2= 1202 H= 2 R= -200
HY（G）= 75.15295764 D= 58.78
P= 1 L= 280 K= 272.985
W= 905.365 V= 514.788 S= 20
A= 900 B= 500
成果：
ZH（G）YH = 58.31372804 过本曲线YH切线方位角
E = 56.40197921 （所求点切线方位角）
中桩：XP= 909.150 = 66°08 ′45.38 ″ 测站至所求点方位角
YP= 520.694 S = 22.626 测站至所求点边长
左桩：XP= 925.809 = 20°27 ′17.09 ″ 测站至所求点方位角
YP= 509.626 S = 27.546 测站至所求点边长
右桩：XP= 892.492 = 103°18 ′03.15 ″ 测站至所求点方位角
YP=531.761 S = 32.636 测站至所求点边长
- 本文固定链接: https://www.chtfs.com/1723/
- 转载请注明: dawang 于 测绘途夫 发表
