当前位置：文档之家 › 线路坐标正算及反算程序

线路坐标正算及反算程序

  • 格式：doc
  • 大小：49.00 KB
  • 文档页数：3

一.程序清单:
1.主程序(TYQXJS)
"K0="？O："X0="？U："Y0="？V："FWJ="？G："LS="？H："+1，0，-1="？Q：If 0= Q：Then 1e45→P：1e45→R：Goto 0：Else "R1="？P："R2="？R：Goto 0：←┘
Lbl 0：(P-R)÷(2HPR) →D："Z=1,F=2"? N：If 1=N：Then Goto 1：Else Goto 2：←┘
Lbl 1："SK="？S："JZ="？Z：Abs(S-O) →W：Prog "SUB1"：F-90→F："X="：X◢"Y="：Y◢"FWJ="：F◢Goto 1◢
Lbl 2："CX="？X："CY="？X→I：Y→J：Prog "SUB2"："O+W→S：Cls：Locate 1，1，"SK"←┘Locate 6，1，S←┘Locate 1，3，"JZ" ←┘Locate 6，3，Z←┘Goto 2◢
2. 正算子程序(SUB1)
0.1739274226→A：0.3260725774→B：0.0694318442→K：0.3300094782→L：1-L→F：1-K→M：
U+W(Acos(G+57.29577951QKW(1/P+KWD))+Bcos(G+57.29577951QLW(1/P+LWD))+Bco s(G+57.29577951QFW(1/P+FWD))+Acos(G+57.29577951QMW(1/P+MWD))) →X：
V+W(Asin(G+57.29577951QKW(1/P+KWD))+Bsin(G+57.29577951QLW(1/P+LWD))+Bsin (G+57.29577951QFW(1/P+FWD))+Asin(G+57.29577951QMW(1/P+MWD))) →Y：
G+57.29577951QW(1/P+WD)+90→F：X+ZcosF→X：Y+ZsinF→Y：Return
3. 反算子程序(SUB2)
G-90→T：Abs((Y-V)cosT-(X-U)sinT)→W：0→Z：Lbl 0：Prog "SUB1"：T+Q57.295779 51W(1/P+WD) →L：(J-Y)cosL-(I-X)sinL→Z：If AbsZ<1E-6：Then Goto1：Else W+Z→W：Goto 0←┘
Lbl 1：0→Z：Prog "SUB1"：(J-Y)÷sinF→Z：Return
二.增设数据库程序（ＳＪＫ）
Lbo0：｛S｝：If S∠下一线元起点里程Then O=起点里程：U=起点X：V=起点Y：G=起算方位角：H=元长度：P=起点半径：R=终点半径：Q=0或1、-1：Prog“TYQX
JS”：If End←┘
三、使用说明
1、规定
(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，
Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。

(2) 当所求点位于中线时，Z=0；当位于中线左铡时，Z取负值；当位于中线中线右侧时，Z取正值。

(3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。

(4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆
弧的半径。

(5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45 次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。

止点与直线相接时，曲率半
径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。

(6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。

止点与直线相接时，曲率半径等
于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。

2、输入与显示说明
输入部分：
1. SZ => XY
2. XY = > SZ
N ? 选择计算方式，输入1表示进行由里程、边距计算坐标；输入2表示由坐标反算里程和边距。

X0 ？线元起点的X坐标
Y0 ？线元起点的Y坐标
S0 ？线元起点里程
F0 ？线元起点切线方位角
LS ？线元长度
R0 ？线元起点曲率半径
RN ？线元止点曲率半径
Q ？线元左右偏标志(左偏Q=-1，右偏Q=1，直线段Q=0）
S ？正算时所求点的里程
Z ？正算时所求点距中线的边距(左侧取负，值右侧取正值，在中线上取零)
X ？反算时所求点的X坐标
Y ？反算时所求点的Y坐标
显示部分：
XS=×××正算时，计算得出的所求点的X坐标
YS=×××正算时，计算得出的所求点的Y坐标
FS=×××正算时，所求点对应的中线点的切线方位角
S=×××反算时，计算得出的所求点的里程
Z=×××反算时，计算得出的所求点的边距
说明：正算和反算子程序单独输入，不与主程序输在一起。

合集下载

坐标正反算程序

坐标正反算程序

1.坐标反算(ZBFS)主程序XY：N=X: T=Y:Pros"DAT1"W"X0"：S"Y0"：O"K0"：G"F0"：H"KN"：P"R0"：R"RN"：Q”Q(-Z +Y)” :D=(P-R)÷(2 Abs (H-O)PR):Z=Abs((T-S)cos(G-90)－(N-W)sin(G-90))：L=0：M”M(YJJ)”=90: （注：此处若不给M赋值，则可计算斜交点）Lbl 0：Prog " SUB1 "：L=(T-Y)cos(G-90+QZ(1÷P+ZD)×180÷π)-(N-X)sin(G-90+QZ(1÷P +ZD) ×180÷π)：AbsL<1E-6=>Goto1：≠>Z=Z+L：Goto 0Δ←┘Lbl 1：L=0：Prog " SUB1 "：L=(T-Y)÷sinF:”K=”:K=O+Z◢”L=”:L=L2.正算主程序 ( ZBZS)W"X0"：S"Y0"：O"K0"：G"F0"：H"KN"：P"R0"：R"RN"：Q”Q(-Z +Y)” :D=(P-R)÷(2Abs(H-O)PR)：L”L(-Z +Y)” :M”ANG”=90:（注：此处若不给M赋值，则可计算斜交点）Z=Abs(K-O)：Prog"SUB1"：”FWJ=”:F=F-M:”X=”:X=X◢”Y=”:Y=Y3.正算子程序(SUB1)Defm4:A=0.1184634425:B=0.2393143352:Z[4]=0.2844444444:C=0.046910 0770:E=0.2307653449:Z[1]=0.5:X=W+Z(Acos(G+QCZ(1÷P+CZD)×180÷π)+Bcos(G+QEZ(1÷P+EZD)×180÷π)+Z[4]cos(G+QZ[1]Z(1÷P+Z[1]ZD)×180÷π)+Bcos(G+Q(1-E)Z(1÷P+(1-E)Z D)×180÷π)+Acos(G+Q (1-C)Z(1÷P+(1-C)ZD) ×180÷π)):Y=S+Z(Asin(G+QCZ(1÷P+CZD)×180÷π)+Bsin(G+QEZ(1÷P+EZD)×180÷π) +Z[4]sin(G+QZ[1]Z(1÷P+Z[1]ZD)×180÷π)+Bsin(G+Q(1-E)Z(1÷P+(1-E)ZD)×180÷π)+Asin(G+Q (1-C)Z(1÷P+(1-C)ZD) ×180÷π))：F=G+QZ(1÷P+ZD) ×180÷π+M：X=X+LcosF：Y=Y+LsinF4.曲线元要素数据库：DAT-01K≥O=>K＜H=> W=**:S=**:O=**:G=**:H=**:P=**:R=**:Q=**⊿⊿←┘K≥O=>K＜H=> W=**:S=**:O=**:G=**:H=**:P=**:R=**:Q=**⊿⊿←┘K≥O=>K＜H=> W=**:S=**:O=**:G=**:H=**:P=**:R=**:Q=**⊿⊿←┘K≥O=>K＜H=> W=**:S=**:O=**:G=**:H=**:P=**:R=**:Q=**⊿⊿←2、输入与显示说明（1）输入部分：X0 ？线元起点的X坐标Y0 ？线元起点的Y坐标K0 ？线元起点里程F0 ？线元起点切线方位角KN ？线元终点里程R0 ？线元起点曲率半径RN ？线元止点曲率半径Q ？线元左右偏标志(左偏Q=-1，右偏Q=1，直线段Q=0)K ？正算时所求点的里程L ？正算时所求点距中线的边距(左侧取负值，右侧取正值，在中线上取零) ANG？正算边桩时左右边桩连线与线路中线的右交角X ？反算时所求点的X坐标Y ？反算时所求点的Y坐标M ? 斜交右角线元要素数据库中K≥O=>K＜H=>中的O和H分别为该段线元起点里程和终点里程A、B、Z[4] 是Gauss-Legendre求积公式中的插值系数C 、E、Z[1] 是Gauss-Legendre求积公式中的求积节点（2）显示部分：X=×××正算时，计算得出的所求点的X坐标Y=×××正算时，计算得出的所求点的Y坐标K=×××反算时，计算得出的所求点的里程L=×××反算时，计算得出的所求点的边距。

5800公路坐标正反算程序(线元法数据输入)

5800公路坐标正反算程序(线元法数据输入)

坐标正反算程序程序变量说明：D：选线。

E、F:置仪点坐X、Y。

K：待算点桩号。

B：待算点距离。

Q：待算点与线路交角（以顺时针方向）。

L：线元起点至待算点长度。

O：待算点方位角。

S：线元长度。

U：线元起点半径。

V：线元终点半径。

程序运行内容说明：开机运行程序后，1、进入主程序，选择正算或反算，1为正算，2为反算。

2、选择线路，1为第线。

2为第线。

3为第线等等。

3、如果上一次计算过的话提示是不更新数线路参数，十秒钟不操作更新，按EXE不更新。

如果上一次不是计算现在要计算的线路的话直接更新数据库（数据库更新会倒计时提示）。

3、正算。

提示入入置仪点坐标。

提示输入待算点桩号。

提示输入离中桩距离。

提示输入与线路方向夹角。

结果：坐标X、Y。

距置仪点距离和方位角。

再次提示输入待算点桩号。

反算。

提示输入实测点坐标X、Y。

结果：待算点桩号程距中桩距离。

再次提示输入实测点坐标X、Y。

程序特点：1、输入参数少，只要输入线路起点的桩号、半径、方位角、XY坐标以及各线元的桩号和半径(左+右-直线×1045)即可。

2、调用变量少，本程还有A、G、H、M、N、P、R、W字母变量未使用，可以用于其它程序而不会和本程序产生冲突。

3、数据库独立，更换工程时数据更换方便。

4、界面友好，尽量做到人人都容易上手，操作一学就会。

5、程序尽量运用5800计算器各种指令，使程序速度提高。

并能够起到帮助大家学习5800各指令工作原理，以便更好掌握编程。

5800P公路路线坐标正反算及高程计算程序

5800P公路路线坐标正反算及高程计算程序

公路路线座标正反算及高程计算程序Fx-5800计算机程序1.FYJS 主程序Fix 3: “1.SZ=>NE:2.NE=>SZ:3.S=>H”:?Q:Q=1 => Goto1:Q=2 => Goto2: Q=3 => Goto3←┘Lbl 1: ＂K＂?S: ＂-B,0,+B＂?Z: “PJ”?M:Prog“5.DAT-SJK”:Prog“2.SUB1”: N+Zcos (F+M) →N:E+Zsin(F+M) →E←┘Prog“4.FY”：Goto 1←┘Lbl 2: “X”?B: “Y”?C:＂K＂?S: Pro g“3.SUB2”: “K=”:S◢“-B,0,+B=”:Z◢Goto 2←┘Lbl 3: ＂K＂?S ：Prog“5.1GCJSCX”: Goto 3＂K＂?S: 所计算点桩号＂-B,0,+B＂?Z: 计算点左右边距左- 右+ 中0“PJ”?M 所计算点斜距的右夹角“X”?B: 反算实测X坐标“Y”?C:反算实测Y坐标＂K＂?S: 反算点附近大概桩号2.SUB1 正算子程序0.5（1÷R-1÷P）÷L→D:S-O(字母 )→X←┘U+∫(cos(A+(X÷P+DX^(2))×180÷π,0,X)→N←┘V+∫(sin(A+(X÷P+DX^(2))×180÷π,0,X)→E←┘A+(X÷P+DX^(2))×180÷π→F←┘3.SUB2 反算子程序Lbl 1:0→Z：1→Q：Prog“5.DAT-SJK”: Prog“2.SUB1” :Pol(N-B+10^(-46), E-C+10^(-46)):Isin(F-90-J) →W:S+W→S:Abs(W)>0.0001 => Goto1←┘Lbl 2: 0→Z：Prog“2.SUB1”:(C-E) ÷sin(F+90) →Z4.FY 放样子程序"HS:X"?N:"HS:Y"?E :"CE:X"?G:"CE:Y"?H:Pol(N-G,E-H): J<0=>J+360→J: " FWJ=": J►DMS◢"JL=": I◢"HS:X"?N: 后视点（计算点）X坐标"HS:Y"?E : 后视点（计算点）Y坐标"CE:X"?G: 测站点X坐标"CE:Y"?H: 测站点Y坐标5.1高程计算子程序（5.1GCJSCX）Prog“ 6.SQX-DAT”:Z[6] Abs(Z[2]- Z[1]) ÷2→Z[3]:if Z[2]- Z[1]>θ :Then 1→Z[4]:Else -1→Z[4]:ifendif S<K:Then - Z[1]→T:Else Z[2]→T:ifendAbs (K-S)→Z[5]if Z[5]- Z[3]>θ:Then Z[7]+T Z[5]→Z[8]:“HS=”：Z[8]◢Else Z[7]+T Z[5]+ Z[4] (Z[3]- Z[5])2÷2÷Z[6]→Z[8]: “HS=”：Z[8]◢ifend5.DAT-SJK 数据库子程序Goto1←┘同时保存多个曲线时的指针Lbl 1←┘IF S<476533.296（线元终点里程）:Then31°24'29.07″→A（线元起点方位角）:476088→O（线元起点里程）:4689534.997→U（线元起点X）:479874.003→V（线元起点Y）:1045→P（线元起点曲率半径）:1045→R（线元终点曲率半径）:445.296→L（线元起点至终点长度）: R eturn:IfEnd←┘IF S<476803.296:Then31°45'48.37″→A:476533.296→O:4689913.601→U:480108.413→V:1045→P:12000→R: 270→L: Return:IfEnd←┘IF S<477106.771:Then32°24'28.85→A:476803.296→O:4690142.627→U:480251.404→V:12000→P:12000→R: 303.475→L: Return:IfEnd←┘IF S<477376.771:Then33°51'25.2″→A:477106.771→O:4690396.754→U:480417.272→V:12000→P:1045→R: 270→L: Return:IfEnd←┘IF S<488923.776:Then34°30'05.68→A:477376.771→O:4690619.834→U:480569.372→V:1045→P:1045→R: 11547.005→L: Return:IfEnd←┘IF S<489353.776:Then34°30'05.68″→A:488923.776→O:4700135.844→U:487109.930→V:1045→P:-10000→R:430→L: Return:IfEnd←┘IF S<490817.499:Then33°16'10.99″→A:489353.776→O:4700491.941→U:487350.944→V:-1000 0→P:-10000→R: 1463.723→L: Return:IfEnd←┘IF S<491247.499:Then24°52'59.52″→A:490817.499→O:4701770.051→U:488061.643→V:-1000 0→P:-1045→R: 430→L: Return:IfEnd←┘IF S<509336.037:Then23°39'04.83″→A:491247.499→O:4702162.678→U:488236.961→V:-1045→P:-1045→R: 18088.538→L: Return:IfEnd←┘程序名：(6.SQX-DAT)(高程竖曲线要素子程序)if S<332.599:Then 0.02793→Z[1]:-0.03309→Z[2]:1800→Z[6]:235→K:32.7→Z[7]: goto 1:ifendif S <734.085:Then -0.03309→Z[1]:-0.00485→Z[2]:6000→Z[6]:460→K:25.254→Z[7]: goto 1:ifend-- Lbl 1Z[1]—前坡: Z[2]—后坡: Z[6]—竖曲线半径 K—变坡点桩号: Z[7]—变坡点设计高三、使用说明1、规定(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当曲线半径在左时，P、R取负值，当曲线半径在右时，P、R取正值，当曲线半径为无穷大（即直线）时，P、R以10的45次代替。

公路路线坐标正反算通用程序

公路路线坐标正反算通用程序

公路路线坐标正反算通用程序（万能程序）此程序为缓和曲线在坐标系下任意位置的通用正反算程序，能完全适应缓和曲线左偏、右偏、Rs>Re、Rs<Re等各种情况，不必先凑成完整缓和曲线，降低算法的复杂程度。

虽然此公式是由缓和曲线推导出来，也适于直线和圆曲线，可降低计算机编程的复杂程度主程序：QXJS-0Lbl4:”1.KS=>XY”:”2. XY=> KS”?Q:”K”?K: Prog”QXJS-SUB0”Lbl0:Q=1=>Goto1: Q=2=> Goto2Lbl1:”S”?S:Prog”QXJS-SUB1”:IfF﹤0:thenF+360→F:Ifend”X=”:Locate5,4,N:”Y=”:Locate5,4,E:”FWJ=”:F▲Lbl2:”X”?B:”Y”?C:B →N:C →E: Prog”QXJS-SUB2”:”K=”:K ▲”S=”:S▲Goto4正算子程序：QXJS-SUB10.5（1÷R-1÷P）÷L→D:K-O→XU+∫(cos（A+(X÷P+DX²)×180÷∏,0,X）→NV+∫(sin(A+(X÷P+DX²)×180÷∏,0,X）→EA+ (X÷P+DX²)×180÷∏→FN+Scos(F+90) →N:E+Ssin（F+90) →E反算子程序：QXJS-SUB2Lbl1:0→S:1→Q: Prog”QXJS-SUB0”: Prog”QXJS-SUB1”Pol(N-B+10∧(-46),E-C+10∧(-46)):Isin(F-90-J) →W:K+W→KAbs（W）>0.0001=>Goto1Lbl2：0→S: Prog”QXJS-SUB1”:(C-E)÷sin(F+90) →S数据库子程序：QXJS-SUB0Lbl1:If K<(线元终点里程): then ***→A（线元起点方位角）：***→O(线元起点里程)：***→U（线元起点X）：***→V（线元起点Y）: ***→P（线元起点曲率半径）: ***→R（线元终点曲率半径）: ***→L（线元起点至终点长度）:Return：Ifend…………………………………每增加一个线元增加一行语句，每增加一条曲线增加一个Lbl数据库子程序：QXJS-A,C,E,G,H,J,T,。