casio 5800p 計算程式
對程式進行部分更新,匝道計算比以前的公式精度高,程式公式簡潔,缺點是計算速度稍慢。
2010/5/24
匝道座標計算程式
a03zdzb 主程式名
lbl 0:」a=1,b=2,c=3,d=4,e=5,h=6,j=7」?a:」k=」?k
if a=1 then prog 「a」:goto 1:ifend a匝道資料庫
if a=2 then prog 「b」:goto 1:ifend b匝道資料庫
if a=3 then prog 「c」:goto 1:ifend c匝道資料庫
if a=4 then prog 「d」:goto 1:ifend d匝道資料庫
if a=5 then prog 「e」:goto 1:ifend e匝道資料庫
if a=6 then prog 「h」:goto 1:ifend h匝道資料庫
if a=7 then prog 「j」:goto 1:ifend j匝道資料庫
lbl 1:prog 「zdjs」:goto 0 主程式完。zdjs為計算過程子程式
zdjs子程式
2→dimz:k-f→i:(d-c)/(z-f)*i→l:v+(2c)*i*90/π→p:if p<0:then
360+p→p:else if p>360:then p-360→p:ifend:ifend:「hr=」:pdms
fix 3:「x=」:n+i*(cos(v)+4∑(cos(v+7.
5(2h+1)*i*((2h+1)*l/12+2c)/ π),h,0,5)+2∑(cos(v+15hi(hl/6+2c)/ π),h,1,5)+cos(p))/36→x
「y=」:e+i*(sin(v)+4∑(sin(v+7.5(2h+1)*i*((2h+1)*l/12+2c)/ π),h,0,5)+2∑(sin(v+15hi(hl/6+2c)/ π),h,1,5)+sin(p))/36→y
norm 1:「hj=」?g:「dl=」?m:fix 3:「xl=」:x+mcos(p-g)→z[1]
「yl=」:y+msin(p-g)→z[2]
norm 1:「dr=」?m:fix 3: 「xr=」:x+mcos(180+p-g)→x
「yr=」:y+msin(180+p-g)→y
norm 1: 「10x」?s:k+10s→k 程式結束
匝道資料庫輸入範例:
a 直線、圓曲線、緩和曲線(包括卵形曲線)一種線形為一組要素
if k≤167.43:then 0:→f:167.43→z:4218460.504→n:
415641.464→e:39°48′47.5″→v:-1/62.75→c:0→d:ifend
if k>167.43: then 167.43→f:334.684→z:4218611.934→n:
415610.526→e:323°22′28.4″→v:0→c:1/300→d:ifend
if k>334.684: then其他匝道資料輸入同此
程式說明:
a= 計算a匝道輸入1,b輸入2,c輸入3.........
k=待求點里程
hr=待求點的方位角
x=待求點中線x座標
y=待求點中線y座標
hj=中線與左邊夾角
dl=中線到左邊的距離
dr=中線到右邊的距離
10x=樁號自動遞增或遞減,如果是計算20m遞增的里程,可以輸入2,計算迴圈後,k自動遞增,輸入-2自動遞減。
資料庫字母含義:
f=線元起點里程
z=線元終點里程
n=起點座標x
e=起點座標y
v=起算點方位角
c=起點曲線半徑的倒數,直線為0,曲線左轉為負,右轉為正
d=終點曲線半徑的倒數,直線為0,曲線左轉為負,右轉為正
主線座標程式
a01zxzb主程式名
24->dimz :norm 1:lbl 9: 「li cheng=」?k:prog 「b01zb」:
t^2/r→z[1]:u^2/r→z[2]:z[1]^2/24/r→z[3]:z[1]/2-z[1]^(3)/240/r^2
→z[4]:z[2]^2/24/r-z[5]:z[2]/2-z[2]^(3)/240/r^2→z[6]:
(r+z[3])tan(abs(w)/2)+z[4]→z[7]:
(r+z[5])tan(abs(w)/2)+z[6]→z[8]:180z[1]/2/r/π→z[9]:
180z[2]/2/r/π→z[10]:r*(abs(w)-z[9]-z[10])*π/180→z[11]:
f-z[7]→a:a+z[1]→b:b+z[11]→ c:c+z[2]→d
if w<0:then -1→j:else 1→j:ifend
if k:else v+w→z[12]:z[8]→z[13]:z[5]→z[22]:z[6]→z[23]:1→z:ifend
n+zz[13]cos(z[12])→z[14]:e+zz[13]sin(z[12])→z[15]
kgoto 1
kgoto 2
kgoto 3
kgoto 4
k≥d=>goto 5
lbl 1:0→z[16]:z[12]+jz[16]→z[17]:(a-k)cos(z[17])→z[20]:
(a-k)sin(z[17])→z[21]:goto 6
lbl 2:30(k-a)^2/r/z[1]/π→z[16]:z[12]+3jz[16]→z[17]:
(k-a)-(k-a)^(5)/90/r^2/z[1]^2→z[18]:
z[18]cos(z[12]+jz[16])→z[20]:
z[18]sin(z[12]+jz[16])→z[21]:goto 6
lbl 3:90z[1]/r/π+180(k-b)/r/π→z[16]:z[12]+jz[16]→z[17]:
rsin(z[16])+z[23]→z[18]:r(1-cos(z[16]))+z[22]→z[19]:
z[18]cos(z[12])-jz[19]sin(z[12])→z[20]:
z[18]sin(z[12])+jz[19]cos(z[12])→z[21]:goto 6
lbl 4:30(k-d)^2/r/z[2]/π→z[16]:z[12]-3jz[16]→z[17]:
(d-k)-(d-k)^(5)/90/r^2/z[2]^2→z[18]:-z[18]cos(z[12]-jz[16])→
z[20]:-z[18]sin(z[12]-jz[16])->z[21]:goto 6
lbl 5:0→z[16]:z[12]-jz[16]→z[17]:
(k-d)cos(z[12])→z[20]:(k-d)sin(z[12])→z[21]:goto 6
lbl 6:norm 1: 「hr=」:z[17] dms
fix 3: 「x=」:z[14]+z[20]→x
「y=」:z[15]+z[21]→y
norm 1: 「hj=」?g: 「dl=」?m:fix 3: 「xl=」:x+mcos(z[17]-g)→i
「yl=」:y+msin(z[17]-g)→l
norm 1: 「dr=」?m : fix 3 : 「xr=」 : x+mcos(z[17]-g+180)→x
「yr=」 : y+msin(z[17]-g+180)→y
norm 1: 「10x」?s:k+10s->k:goto 9 程式結束
資料庫輸入範例:
b01zb 資料庫
if k≤72905.727:then 71847.992→f: 4219571.89→n:421323.274
→e: 0→t:0→u: 5000→r: 278°02′47.63″→v: -12°11′22.16″→w:ifend
if k>72905.727:then 73528.28→f: 4219450.207→n:419643.357
→e: 0→t:0→u: 4000→r: 265°51′25.48″→v:17°41′34.34″→w:ifend
if k
程式說明:
k=待求點里程
hr=待求點的方位角
x=待求點中線x座標
y=待求點中線y座標
hj=中線與左邊夾角
dl=中線到左邊的距離
dr=中線到右邊的距離
10x=樁號自動遞增或遞減
資料庫字母含義:
f=交點里程
n=起點座標x
e=起點座標y
t=前緩和曲線引數a
u=後緩和曲線引數a
r=半徑
v=前切線方位角
w=轉角,左轉輸入負值,右轉正值
Casio5800交點法與線元法
casio5800交點法與線元法 積木法 匝道座標正反算放 樣程式顏色文字casio5800交點法與線元法 積木法 匝道座標正反算放樣程式 xufeng 2011.2.14 本人一直以來想找乙個交點法與線元法相結合的座標正反算程式,在網上找了很久很久,沒能找到乙個較為滿意的,有幸在測量空間看到大歪哥...
5800應用程式
一 座標計算源程式 1.主程式 tyqxjs xy sz n u x0 v y0 o s0 g f0 h ls p r0 r rn q c 1 p d p r 2hpr e 180 n 1 goto 1 go to 2 lbl 1 sz w abs s o prog sub1 x xs x y ys...
CASIO 4800常見程式序集
序集1 由乙個已知點座標和乙個已知方位角計算支導線座標 程式名 zbjs 程式清單 x y t lb1 0 fixm t t a 180 t 0 t t 360 t 360 t t 360 t x x scos t y y ssin t goto0 程式執行說明 執行 zbjs 每執行一步按exe鍵...