1÷P→C:(P-R)÷(2HPR) →D: 180÷π→E←┚
0.1739274226→A: 0.3260725774→B: 0.0694318442→K: 0.3300094782→L: 1-L→F: 1-K→M:U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Bcos(G+QEFW(C+FWD))+Acos(G+QEMW(C+MWD)))→X: V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+QEFW(C+FWD))+Asin(G+QEMW(C+MWD)))→Y: G+QEW(C+WD)+90→F: X+ZcosF→X: Y+ZsinF→Y
3反算子程序(FS)(不运行)
G-90→T: Abs((Y-V))cos(T)-(X-U)sin(T)→W: 0→Z←┚
Lbl 0←┚
prog"ZS":T+QEW(C+WD)→L:(J-Y)cosL-(I-X)sinL→Z:If Abs(Z)<0.001: Then Goto1:Else W+Z→W:IfEnd :Goto 0←┚
Lbl 1: 0→Z:Prog"ZS":(J-Y)÷sinF→Z
ZX-WN(运行主程序,该程序适用于匝道,对全段完整非对称曲线计算太过麻烦)
"CZX"?S(仪站X坐标):"CZY"?W(仪站Y坐标)
LbI 0: "K"?K: If K≤126.9(第1条单一曲线的终点下面同上): Then Goto A:Else If K≤176.9 And K>126.9: Then Goto B:Else If K≤259.562 And K>176.9: Then Goto C:Else If K≤309.562 And K>259.562: Then Goto D: IfEnd:IfEnd:IfEnd:IfEnd←┚
LbI A:92094.336(起点X坐标)→U:8287.962(起点Y坐标)→V:0起点桩号→O:145°11'54"(方位角)→G:126.9(曲线长度)→H:10^(45) (起点半径)→P:10^(45) (止点半径)→R:+1(路线左右偏,左负右正,直线输0)→Q:Prog "WN":Goto W
LbI B:91990.134→U:8360.388→V:126.9→O:145°11'54"→G:50→H:10^(45) →P:150 →R:+1→Q:Prog "WN":Goto W
LbI C:91947.609→U:8386.570→V:176.9→O:154°44'51"→G:82.662→H:150 →P:150 →R:+1→Q:Prog "WN":Goto W
LbI D:91867.101→U:8399.992→V:259.562→O:186°19'19"→G:50→H:150 →P: 10^(45) →R:+1→Q:Prog "WN":Goto W
LbI W:"PJ1"?T←┚(该项是只左右边桩再左右偏距,输0为不再左右偏,可以直接出方位角、距离,适用于桥梁桥台、盖梁、涵洞八字墙等放样,以路线前进方向,前-,后+)
If T=0:Then Z[3]→Z: Z[4]→M: Goto P: Else Goto V: IfEnd←┚
LbI V:"JJ1"? O←┚(左右边桩再偏夹角,用180°减去路线夹角)
"X1=": Z[3]+Tcos(F+N+ O)→U◢
"Y1=" :Z[4]+Tsin(F+N+ O)→V◢
U→Z: V→M: GotoP←┚
LbI P:Pol(Z-S,M-W): If J≤0:Then J+360°→J: IfEnd←┚
"JL=":I◢(仪站至测点的距离)
"FWJ=":J DMS◢(仪站至测点的方位角)
Goto0
备注:
1. 规定以道路中线前进方向(即里程增大的方向)区分左右;当线元往左偏时Q=-1;当元线往右偏时,Q=1:当元线为直线时,Q=0。
2. 当线元为直线时,其起点、止点的曲率半径为无穷大,以10的45次方代替。
3. 当线元为圆曲线时,无论其起点、止点与什么线元相接,其曲率半径均等于
说明:本程序适用于计算器 CASIO fx-5800P,可计算:第一夹直线段、第一缓和曲线段、圆曲线段、第二缓和曲线段、第二夹直线段的线路中、边桩坐标及切线方位角,反算里程和偏距。数据库名:YTGoto 1Lbl 1:If Z<10000:Then 3378226.731→U:456053.721→V:3373351.385→N:458628.379→E:9557.84→D:2040→R:260→S:280→T:-14°18ˊ10″→A:Return:IfEndLbl 1:If Z<………………………………………………………………………………………………………………Lbl 1:If Z<本交点HZ里程或下个交点ZH里程:Then 前交点坐标X→U:前交点坐标Y→V:本交点坐标X→N:本交点坐标Y→E:交点里程→D:曲线半径→R:第一缓和曲线长→S:第二缓和曲线长→T:转角→A:Return:IfEndLbl 1每增加一行则为增加一个交点要素,……表示下一个交点要素。注:1、◢ 为显示指令。2、 ↙为回车(EXE)标识。3、()括号内为注解不用输入。计算参数:1=>XY2=>DK
J? 2(输入1表示调用数据库,输入2表示手动输入数据)
1、QJD(X)=? 3378226.731(前交点坐标X)
2、QJD(Y)=? 456053.721(前交点坐标Y)
3、JD(X)=? 3373351.385(本交点坐标X)
4、JD(Y)=? 458628.379(本交点坐标Y)
5、JD(DK)=? 9557.84(本交点里程)
6、R=? 2040(曲线半径)
7、LS1?= 260(第一缓和曲线长)
8、LS2?= 280(第二缓和曲线长)
9、α?= -14°18ˊ10″(转角)
10、JSDK?= 9500(计算里程)
11、W?= 10(偏距:左为负,有为正数)
12、α?= 90(偏角)
………………………………
13、FWJ=146°34ˊ4.38″(切线方位角)
14、X= 3373402.498(计算结果坐标X)
15、Y= 458602.6783(计算结果坐标Y)
16、ZH=09+170.838,X=3373693.598,Y=458447.6569,FWJ=152°9′41.68″
17、HY=09+430.838,X=3373466.361,Y=458573.9060,FWJ=148°30′37.36″
18、QZ=09+550.461,X=3373366.244,Y=458639.3431,FWJ=145°9′2.25″
19、HY=09+670.084,X=3373270.134,Y=458710.5365,FWJ=141°47′27.15″
20、HZ=09+950.084,X=3373058.317,Y=458893.5689,FWJ=137°51′31.68″资料地址: https://skydrive.live.com/?cid=23F98B07A5CC3A18&id=23F98B07A5CC3A18%213143
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)