卡西歐4800、4850系列計算器測量計算程序
壹、字母含義;
K: 表示擬計算位置的線路樁號 ;
H: 表示計算位置距路線中心線的偏距,左偏為正,右偏為負。
T: 各段線路上作為起算點處的切線方位角。
S: 擬計算點到起算點的曲線長。
L: 在圓曲線上表示曲線長,在緩和曲線上表示緩和段長度,在直線上為零。
R: 表示曲線半徑,左偏為正,右偏為負
E、F: 起算點的坐標值。
M”X1”N”Y1”: 已知點坐標,求其對應位置樁號及左右偏距。
二、程序
1、 坐標計算(COORD)
{K,H}:KH“L+,R-”:Prog “DATA”:”(X,Y)=”:
X=X+HSinW :Pause 1:Y=Y-HcosW:
2、坐標反算線路樁號(FS ZH)
Fix 4:M”X1”N”Y1”:LbI 1:Prog “DATA”:PoI (M-X,N-Y):Fixm:J<0=>J=J+360:≠>J=J⊿ Abs(Sin(W-J ))=1=>”K=”:K: Pause 1 :
“L+,R- =”: H=ISin(W-J):≠>K=K+Icos(W-J): Goto 1
3、子程序Prog “DATA”
K<(第1段與第2段線路分界處的路線樁號)=> T=(第1段起算點處的切線方位角值):S=K-(第1段起算點處的路線樁號):L=(在圓曲線上等於S;在緩和曲線上等於緩和段長度;在直線上為零。):R=(曲線半徑):E=(第1段起算點的X坐標值):F=(第1段起算點的Y 坐標值)
≠ >K<(第2段與第3段線路分界處的路線樁號):=>T=(第2段起算點處的切線方向角值):S=K-(第2段起算點處的路線樁號):L=(在圓曲線上等於S;在緩和曲線上等於緩和段長度;在直線上為零。):R=(曲線半徑):E=(第2段起算點的X坐標值)F=(第2段起算點的Y坐標值)
≠ >K<(第3段與第4段線路分界處的路線樁號):=>T=(第3段起算點處的切線方向角值):S=K-(第3段起算點處的路線樁號):L=(在圓曲線上等於S;在緩和曲線上等於緩和段長度;在直線上為零。):R=(曲線半徑):E=(第3段起算點的X坐標值)F=(第3段起算點的Y坐標值)
……
≠ >K<(第n-1段與第n段線路分界處的路線樁號):=>T=(第n-1段起算點處的切線方向角值):S=K-(第n-1段起算點處的路線樁號):L=(在圓曲線上等於S;在緩和曲線上等於緩和段長度;在直線上為零。):R=(曲線半徑):E=(第n-1段起算點的X坐標值)F=(第n-1段起算點的Y坐標值)
≠>T=(第n段起算點處的切線方向角值):S=K-(第n段起算點處的路線樁號):L=(在圓曲線上等於S;在緩和曲線上等於緩和段長度;在直線上為零。):R=(曲線半徑):E=(第n段起算點的X坐標值)F=(第n段起算點的Y坐標值)⊿⊿⊿⊿⊿⊿⊿⊿(⊿***n-1個)
L ≠ 0 = > C =90S^2/ЛRL : ≠ > C=0:⊿
L=S => V=2RsinC : Q=T-C : W=T-2C ≠> L=0 =>V=S : Q=T : W=T :≠> Pol(S-S^5/40R^2L^2+S^9/3456R^4L^4,S^3/6RL-S^7/336R^3L^3+S^11/42240R^5L^5):Fixm:V=I:Q=T-J:W=T-CΔΔ
X=E+VcosQ:
Y=F+VsinQ
三、計算示例(浙江龍麗高速公路某標段線路)
計算坐標:COORD
{K,H}:KH“L+,R+”:Prog “DATA”:”(X,Y)=”:X=X+HSinW Pause 1:
Y=Y-HcosW:
計算樁號、偏距:FS ZH
Fix 4:Deg:M”X1”N”Y1”:LbI 1:Prog “DATA”:PoI (M-X,N-Y):Fixm:J<0=>J=J+360:≠>J=J⊿
Abs(Sin(J- W ))=1=>”K=”:K: Pause 1 :
“L+,R- =”: H=ISin(J-W):≠> K=K+Icos(J-W):Goto 1
數據庫:DATA
K<74886.420=>T=151°44?18.1":S=K-73544.990:R=-2500:E=155009.949:F=439722.191: S<0=>L=S:≠>L=0Δ
≠>K<77616.980=>T=131°13?39" : S=K-76676.330:R=5000:E=152435.423: F=441466.078: S<0 =>L=S:≠>L=0Δ
≠>K<80984.800=>T=152°13?22.7" : S=K-79449.180:R=-5000:E=150385.156: F=443302.124:S<0=>L=S:≠>L=0Δ
≠>T=134°00?16.4":S=K-82574.660:R=5000:E=147760.219:F=444968.04: S<0=>L=S:≠>L=0ΔΔΔΔ
L ≠ 0 => C =90S^2/ЛRL : ≠> C=0: Δ
L=S => V=2RsinC : Q=T-C : W=T-2C ≠> L=0 =>V=S : Q=T : W=T :≠> Pol(S-S^5/40R^2L^2+S^9/3456R^4L^4,S^3/6RL-S^7/336R^3L^3+S^11/42240R^5L^5):Fixm:V=I:Q=T-J:W=T-CΔΔ
X=E+VcosQ: Y=F+VsinQ