卡西欧9750公路线元法计算程序
时间:2010-07-21 09:51:32来源:本站作者:sun2001124我要投稿我要收藏投稿指南
.新版程序把线元法和交点法已经集成在一个模块中了，用户只需修改JD程序和ZA程序中的数据部分即可,其余不需作任何的改动。 2.因为每条路高程计算不尽相同，且比较复杂，现在可利用PC机EXCEL计算好打印成表格带到工地上使用，所以本版程序未对线路高程序进行专门的编程计算，而是利用统计计算模式中来输入桩号（第一列X）及左、右高程（第
二、三列Y，Freq），这种输入数据的方式最为直观，易发现错误，也易修改，输入完毕后运行S程序对数据按桩号进行排序，在程序中通过调用GG程序来进行内插计算，SG=-1得左标高，SG=1得右标高（若SG输入0,则可进行一般的线性内插计算）。 3.在JD程序和XY程序中，先将一个计算单元的数据置入矩阵F中（1行8列或1行9列），这样程序可读性极好。 4．相比原CASIO4850程序操作习惯，作了一点小小的改动，测站坐标存在Z[10]，N中，X坐标原存在M中容易被误操作修改，而设计标高存在M中，这样易于修改，因为CASIO5800没有IN，OUT功能，很不方便。
4.程序利用Z[2]变量值来判断是采用交点法还是线元法模型计算，Z[2]＝0为线元法，否则为交点法。
一、PQX程序：计算中边桩坐标及近似的桩号反算，在运行模式直接调用。① Z[10]→S:”XO”?S:S→Z[10]:”YO”?N:Prog “AU”② Lbl 2:?L:Prog “Z”:Prog “E”:1n→O:90→S
③ Lbl 4:”JJ”?S:”YC”?O:SO=0 =>Goto 2‘原来lbl 后没有标号４的。④ O=-1 =>Goto 6
⑤ “X,Y”:R+OCos(Z+S)→X▲U+OSin(Z+S)→Y▲Prog “D”:Goto 4⑥ Lbl 6:Z[7]→X:Z[8]→Y:”XF”?X:”YF”?Y:XY=0 =>Goto 4⑦
X→Z[7]:Y→Z[8]:Pol(X-R,Y-U+1p):Z+S-J→J:”YC,DL,L”:ICos(J)→O▲ISin(J)→I▲L+I▲Goto 6
二、P程序：在程序中提供一个自由运算的模式。
①Lbl 1:”TMP”?I:If I≠0:Then “RST”:I▲Goto 1:IfEnd
二、LYC程序：进行桩号反算及边坡放样，在运行模式直接调用。①Prog “AU”
②Lbl 1:Z[7]→X: Z[8]→Y:
Z[6]→S: ”XF”?X :X→Z[7]:”YF”?Y:Y→Z[8]: ”ZF”?S: S→Z[6]③Lbl 2:Prog “Z”:Y=U =>Y+1p→Y
④Pol(X-R,Y-U):J-Z→J:Isin(J)→O:Icos(J)→I
⑤If Abs(I)≤0.1:Then Prog “E”:”L,YC”:L+I→L▲O▲Goto 3:IfEnd⑥If Z[9]≠0:Then Pol(Z[9]-SO,I):πJZ[9]÷180→I:IfEnd⑦”DL”:I▲L+I→L:Goto 2
⑧Lbl 3: Z[6]→S:If S=0:Then Goto 1:IfEnd
⑧ M→Z
⑨ Lbl 4:”SG”?Z:Z→M:If Abs(Z)=1: Then Prog “GG”:Y→Z:If X=1:Then
X→Z:IfEnd:Z→M:IfEnd
⑾If Z=0:Then Z[11]→C:Z[12]→D: Z[13]→E:
Z[14]→F:”LA”?C:C→Z[11]: ”HA”?D:D→Z[12]: ”LB”?E:E→Z[13]: ”HB”?F:F→Z[14]:D+(E-C)-1(F-D)(L-C)→Z:Z→M:IfEnd
⑿Lbl 5:”SG”:Z:”TW”:Z-S→C▲:Z[3]→P:
Z[4]→Q: :Z[5]→T:”YC0”?P:P→Z[3]: :”HC0”?Q:Q→Z[4] :”M”?T:T→Z
[5]:O÷Abs(O)→I:”BL”:I(P+TC(C+Q)÷Abs(C))→U▲”DB”:IU-IO→J▲Ans→I:Prog “P”:Goto 1
四、I程序：通过手工方式输要素
①Z[2]=0 =>Goto 2
②Z[1]→R:Z[2]→U:”JL”?K:”JX”?G:”JY”?H:”JA”?V:”JJ”?W:”R”?A:”LS”?B:”L1”?R:”L2”?U:R→Z[1]:U→Z[2]:ReTurn
③Lbl 2:”L0”?K:”X0”?G:”Y0”?H:”A0”?V”:”L1”?F:”L-0 Z=-1 Y=1”?W:W=0 =>ReTurn
④”R0”?A:”R1”?B
五、D程序：根据测站点的坐标及放样点的坐标计算水平距离及方位角①Y=N =>Y+1n→Y
②Pol(X-Z[10],Y-N):JJ+360→J
③”A,D”:J→DMS▲I▲
六、E程序：桩号出界判断
①Z[2]=0 =>Goto 2
②If L>Z[2]: Then “>L2”:L▲IfEnd
③If L
④Goto 3
⑤Lbl 2:If L
⑥If L>F:Then “>L1”:L-F▲IfEnd
⑦Lbl 3
七、GG程序：用统计数据中桩号，左右高程数据进行内插计算。
①FreqOff : n→C:FreqOn:If LmaxX: Then “L>maxX”:C-1→I:Goto 8:IfEnd
③For 1→I To C-1
④If L≤List X: Then Break:IfEnd:Next
⑤Lbl 8 : (L-List X)÷(List X- List X)→J : List Y+J(List Y-List Y)→X : List Freq +J(List Freq -List Freq )→Y
八、S程序:对统计数据列表中数据按桩号进行排序。
①FreqOff : n→C:FreqOn“WAIT„n=”:Locate 10,1,C
②If C
③For 1→D To C: D→I:
④For D+1→E To C:
⑤If List X[E]
⑥If I≠D: Then List X[D]→F:List X→List X[D]:F→List X: List
Y[D]→F:List Y→List Y[D]:F→List Y: List Freq[D]→F:List Freq→List Freq [D]:F→List Freq :IfEnd:Next
⑦Locate 1,2,”minX:”Locate 6,2,minX
⑧Locate 1,3,”minX:”Locate 6,3,maxX
九、X程序:计算缓和曲线坐标增量xx,yy
①AB→U:U÷O→Z[9]:0.5O2÷U→R:O-0.1OR2→I:6-1OR(2-7-1R2)→J:Rr→R
十、JS程序:内插计算的接口程序，在运行模式直接调用。
①Prog “S”
②Lbl 1:”L=”?L:L=0 =>Return
③Prog “GG”:”X,Y”:X▲Y▲Goto 1
十一、GZW程序：构造物坐标放样计算，在运行模式直接调用。
①Z[10]→S:”XO”?S:S→Z[10]:”YO”?N：0→T
②Lbl 1;”X0”?R:”Y0”?U:”A0”?Z:IF Z=-1 :Then
“A0,DD”:J→Z▲I▲IfEnd
③If Z=0:Then R→X:U→Y:Prog “D”:Goto 1:IfEnd
\'在④后增加一句，解决了程序再次运行时，交角的初值被置为测站点的Ｘ坐标的BUG.
④O-Z→S:”J0”?S:Z+S→O:”0-XX,YY→X,Y”:”1-X,Y→XX,YY”?T:If T≠0: Then Goto 3:IfEnd
⑤Lbl 2:”XX”?P:”YY”?Q:If PQ=0:Then Goto 1:IfEnd
⑥”X,Y”:R+Pcos(Z)+Qcos(O)→X▲U+Psin(Z)+Qsin(O)→Y▲Prog “D”:Goto 2
⑦Lbl 3:”XF”?X:”YF”?Y:If XY=0: Then Goto 1:IfEnd
⑧”XX,YY”:((X-R)sin(O)-(Y-U)cos(O))÷sin(S)→P▲((Y-U)cos(Z)-(X-R)sin(Z))÷sin(S)→Q▲
⑨Goto 3
十二、Y程序：线元法计算匝道坐标
①L-K→C:If L>F:Then F-K→C:IfEnd
②If W=0 or C
④If B=0:Then (A-1C)r→C:Asin(C)→I:A(1-cos(C))→J:A→Z[9]:Goto 7:IfEnd
⑤A-B→S:B(F-K)÷Abs(S)→R:AR→U:S÷Abs(S)→S:0.5R2÷U→D:R+SC→T:⑥0.5T2÷U→E:U÷T→Z[9]:T(1-0.1E2+E22÷216-E^(6)÷9360)-R(1-0.1D2+D22÷216-D^(6)÷9360)→I:
⑦3-1(TE(1-E2÷14+E22÷440-E^(6) ÷25200)-RD(1-D2÷14+D22÷440-D^(6)÷25200))→J:Pol(1n+Abs(I),Abs(J):Rec(I,Abs(J-Dr):Abs(E-D)r→C⑧Lbl 7:G+Icos(V)-WJsin(V)→R:H+Isin(V)+WJcos(V) →U:V+WC→Z⑨Lbl 8:If L>F:Then R+Rec(L-F,Z) →R:U+J→U:0→Z[9]:IfEnd
十三、Z程序：交点法计算中桩坐标
①If Z[2]=0:Then Prog “Y”:Return:IfEnd
②B2÷A÷24→P:0.1B(5-P÷A) →Q:Q+(A+P)tan(0.5Abs(W))
→T:K-T→C:C+B→D:D+πAAbs(W) ÷180→F:F-B→E
③W÷Abs(W)→S:0→Z[9]:If L≤C:Then L-K→O:V→Z:Goto 2:IfEnd④If L≥F:Then L-F+T→O:V+W→Z:Goto 2:IfEnd
⑤If L≥E:Then
F-L→O:Prog ”X”:V+W→O:O-SR→Z:G+(T-I)cos(O)-SJsin(O)→R:H+(T-I)sin
(O)+Sjcos(O) →U:Goto 3:IfEnd
⑥If L>D:Then ((L-D+0.5B)÷A)r→O:Q+Asin(O)→I:P+A(1-cos(O))
→J:A→Z[9]:V+SO→Z:Goto 4:IfEnd
⑦L-C→O:Prog “X”:V+SR→Z
⑧Lbl 4:G+(I-T)cos(V)-Sjsin(V)→R:H+(I-T)sin(V)+Sjcos(V)→U:Goto 3⑨Lbl 2:G+Ocos(Z)→R:H+Osin(Z)→U
⑩Lbl 3
十四、AU程序: 注：L0-输测段桩号近似值，”AU”子程序根据桩号近似值自动调用曲线要素（包括线元法和交点法要素），Z[2]=0则调用线元法数据，否则为交点法数据
（L0=0-默认为原曲线要素 L0=-1-手工输入曲线要素）
①Cls:” L0= 0-Conti.”:” L0=-1-Input”:If Z[2]=0:Then Locate 13,2,”(XY)”:Else Locate 13,2,”(JD)”:IfEnd
②Z[2]→S:”0－XY E.－JD”?S:S→Z[2]:”L0”?L:If L=0:Then ReTurn:IfEnd③If L=-1:Then Prog “I”:ReTurn:IfEnd
④If Z[2]=0:Then Prog “ZA”:Else Prog “JD”:IfEnd
⑤Mat F[1,1]→K: Mat F[1,2]→G: Mat F[1,3]→H: Mat F[1,4]→V: Mat F[1,5]→W: Mat F[1,6]→A: Mat F[1,7]→B: If Z[2]=0:Then Mat
F[1,8]→F:Else Mat F[1,8]→Z[11]: Mat F[1,9]→Z[12]:IfEnd
十五、JD程序：交点法矩阵每一行数据依次为交点桩号，交点X坐标,交点Y坐标,起始边方位角，偏角，圆曲线半径,缓和曲线长度,计算单元起点，计算单元终点
①If L＜24900: Then “L＜minX”:L▲Stop:IfEnd
②If L＞30801.27: Then “L＞maxX”:L▲Stop:IfEnd
③IF L≤26227.359:Then [[25676.236,
142225.377,504270.949,6°24°59.6°,-70°50°10.6°,770,300.007,24900,26227.359]]→Mat F:Return:IfEnd
④IF L≤27318.276:Then [[27026.458,
142873.016,502918.028,295°34°49.°,-15°31°18.5°,1500,180.003,26227.359,27318.276]]→Mat F: Return:IfEnd
⑤IF L≤28522.182:Then [[28148.818,
143069.516,501810.234,280°03°30.5°,41°45°02.1°,780,206.532,2731
8.276,28522.182]]→Mat F: Return:IfEnd
⑥IF L≤29246.112:Then [[28895.742,
143666,501341,321°48°32.6°,-22°33°08.8°,1250,215.74,28538.355,29246.112]]→Mat F: Return:IfEnd
⑦IF L≤30801.27:Then [[30431.525,
144420,499995,299°15°23.8°,49°29°52.7°,720,160.598,29246.112,30801.27]]→Mat F: Return:IfEnd
十六、ZA程序
注2：线元法矩阵每一行数据依次为起点桩号，起点X坐标，起点Y坐标，起点方位角，左右转标志，起点曲率半径，终点曲率半径，终点桩号
①If L＜0: Then “L＜minX”:L▲Stop:IfEnd
②If L＞1479.85: Then “L＞maxX”:L▲Stop:IfEnd
③If L＞1000: Then Goto 2:IfEnd
④If L＞553.49: Then “L＞maxX”:L▲Stop:IfEnd
⑤IF L≤82.699:Then
[[0,142872.513,502902.111,290°03°00°,-1,1509.5,1509.5,82.699]]→Mat F:Return:IfEnd
⑥IF L≤111.519:Then
[[82.699,142898.724,502823.687,286°54°40°,-1,1509.5,1T,111.519]]→Mat F:Return:IfEnd
⑦IF L≤155.022:Then
[[111.519,142906.932,502796.061,286°21°51°,-1,1T,1000,155.022]]→Mat F:Return:IfEnd
⑧IF L≤269.54:Then
[[155.022,142918.885,502754.233,285°07°04°,-1,1000,1000,269.54]]→Mat F:Return:IfEnd
⑨IF L≤359.475:Then
[[269.54,142942.363,502642.211,278°33°23°,-1,1000,1T,359.475]]→Mat F:Return:IfEnd
⑩IF L≤434.421:Then
[[359.475,142953.072,502552.924,275°58°48°,1,1T,1200,434.421]]→Mat F:Return:IfEnd
⑾IF L≤463.597:Then
[[434.421,142961.655,502478.474,277°46°09°,1,1200,1200,463.597]]→Mat F:Return:IfEnd
⑿IF L≤501.14:Then
[[463.597,142965.95,502449.617,279°09°44°,1,1200,1T,501.14]]→Mat F:Return:IfEnd
⒀IF L≤553.488:Then
[[501.14,142972.314,502412.618,280°03°31°,0,1T,1T,553.488]]→Mat F:Return:IfEnd
⒁Lbl 2:IF L≤1169.959:Then
[[1000,142957.323,502388.342,102°55°07°,0,1T,1T,1169.959]]→Mat F:Return:IfEnd
⒂IF L≤1215.026:Then
[[1169.959,142919.325,502554,102°55°07°,-1,1T,1500,1215.026]]→Mat F:Return:IfEnd
⒃IF L≤1261.836:Then
[[1215.026,142909.471,502597.975,102°03°28°,-1,1500,1500,1261.836]]→Mat F:Return:IfEnd
⒄IF L≤1283.436:Then
[[1261.836,142900.408,502643.898,100°16°11°,-1,1500,1T,1283.436]]→Mat F:Return:IfEnd
⒅IF L≤1307.44:Then
[[1283.436,142896.66,502665.17,99°51°26°,-1,1T,700,1307.44]]→Mat F:Return:IfEnd
⒆IF L≤1346.219:Then
[[1307.44,142892.686,502688.842,98°52°30°,-1,700,700,1346.219]]→Mat F:Return:IfEnd
⒇IF L≤1363.61:Then
[[1346.219,142887.767,502727.303,95°42°03°,-1,700,1T,1363.61]]→Mat F:Return:IfEnd
(21)IF L≤1387.958:Then
[[1363.61,142886.183,502744.622,94°59°20°,1,1T,500,1387.958]]→Mat F:Return:IfEnd
(22)IF L≤1451.374:Then
[[1387.958,142883.869,502768.858,96°23°02°,1,500,500,1451.374]]→Mat F:Return:IfEnd
(23)IF L≤1477.462:Then
[[1451.374,142872.845,502831.266,103°39°03°,1,500,1490.5,1477.462]]→Mat F:Return:IfEnd
(24)IF L≤1479.848:Then
[[1477.462,142866.175,502856.486,105°38°50°,1,1490.5,1490.5,1479.848]]→Mat F:R
