卡西欧5800交点法坐标计算程序.doc

交点法坐标计算JDFZBJS
第一种：单一交点法计算任意里程位的坐标。
程序如下:
待求段的交点X坐标→V:
待求段的交点Y坐标→W:
前一段的交点X坐标→M:
前一段的交点Y坐标→N:
直缓点ZH的里程→Q:
缓直点HZ的里程→Z:
待求段圆曲线半径→R:
待求段缓和曲线长度→U:
待求段的切线长度→T:
待求段转向方向（左转取-1，右转取+1）→P:
待求段转向偏角→A:
10→DimZ:
Cls:
“KM=”?L:待求点里程
“L-1,R+1”?K:待求点位于前进方向左侧为-1,右侧为+1
“PJ=”?S:待求点到线路设计中心的平距
Lbl0:
Pol(V-M,W-N):
J→Z[1]:
Rec(T,Z[1]+180):
V+I→Z[2]:
W+J→Z[3]:
Lbl1:
IfL≤Q:
ThenRec(Q-L,Z[1]+180):
ElseGoto2:
IfEnd:
Z[2]+I→Z[4]:
Z[3]+J→Z[5]:
IfK=1:
ThenRec(S,Z[1]+90):
ElseRec(S,Z[1]+270):
IfEnd:
Goto6:
Lbl2:
Abs(Q-L)→Z[6]:
IfU≥Z[6]:
ThenPol(Z[6]-Z[6]ˆ(5)÷(40R2U2)+Z[6]^(9)÷(3456(RU)^(4)),Z[6]^(3)÷(6RU)-Z[6]^(7)÷(336(RU)^(3))+Z[6]^(11)÷(42240(RU)^(5))):
ElseGoto3:
IfEnd:
J→Z[7]:
Rec(I,Z[1]+Z[7]P):
Z[2]+I→Z[4]:
Z[3]+J→Z[5]:
Z[1]+Z[7]P+60PZ[6]2÷(RUπ)+90→Z[10]:
IfK=1:
ThenRec(S,Z[10]):
ElseRec(S,Z[10]+180):
IfEnd:
Goto6:
Lbl3:
Abs(Q-L)→Z[6]:
IfL≤Z-U:
ThenPol(Rsin(90(2Z[6]-U)÷(Rπ))+U-Rsin(90U÷(Rπ))-U^(3)÷(40R2),R-Rcos(90(2Z[6]-U)÷(Rπ))+U2÷(24R)):
ElseGoto4:
IfEnd:
J→Z[7]:
Rec(I,Z[1]+Z[7]P):
Z[2]+I→Z[4]:
Z[3]+J→Z[5]:
Z[1]+90P(2Z[6]-U)÷(Rπ)+90→Z[10]:
IfK=1:
ThenRec(S,Z[10]):
ElseRec(S,Z[10]+180):
IfEnd:
Goto6:
Lbl4:
IfL≥Z-UAndL≤Z:
ThenZ[1]+AP→Z[1]:
ElseGoto5:
IfEnd:
Rec(T,Z[1]):
V+I→Z[2]:
W+J→Z[3]:
Z→Q:
-P→P:
-K→K:
Z[1]+180→Z[1]:
Goto2:
Lbl5:
IfL≥Z:
ThenZ[1]+AP→Z[1]:
IfEnd:
Rec(T,Z[1]):
V+I→Z[2]:
W+J→Z[3]:
Rec(L-Z,Z[1]):
Z[2]+I→Z[4]:
Z[3]+J→Z[5]:
IfK=1:
ThenRec(S,Z[1]+90):
ElseRec(S,Z[1]+270):
IfEnd:
Goto6:
Lbi6:
Z[4]+I→Z[8]:
Z[5]+J→Z[9]:
Cls:
Fix3:
“Xp=”:计算的X坐标
“Yp=”:计算的Y坐标
Locate7,1,Z[8]:X坐标显示在第一行
Locate7,2,Z[9]◢Y坐标显示在第二行
Stop

注解:
在本程序中考虑了ZH点前和HZ点后的直线情况。程序最后用Stop结束一个里程位的计算,如果要计算的点多,可以在程序最前面(即程序开始处)加个语句行LblA:,在程序最后面将Stop改写成GotoA。但这样在停止计算时转出来较麻烦。
第二种：单一交点法计算任意里程位的坐标。
主程序JD
24→DimZ:
Cls:
交点里程桩号→K:
交点X坐标→X:
交点Y坐标→Y:
第一缓和曲线长度→B:
第二缓和曲线长度→C:
圆曲线半径→R:
ZH点方位角→M:
偏角(带符号)→O:
M+O→N:HZ点方位角
Prog“JDA”
Cls:
“T1=”:第一切线长度
“T2=”:第二切线长度
“L=”:曲线总长度
“LY=”:圆曲线长度
Locate9,1,S:
Locate9,2,T:
Locate9,3,L:
Locate9,4,Q◢
Cls:
“E=”:外距
“K(ZH)=”:直缓点里程
Locate11,1,E:
Locate8,2,Z[1]◢
Cls:
“K(HY)=”:缓圆点里程
“K(QZ)=”:曲中点里程
“K(YH)=”:圆缓点里程