FX5800计算器交点法正反算程序.doc

FX5800计算器交点法正反算程序

输入图纸中的整个标段交点要素，输入桩号及距中桩距离正算本标段任意点坐标，及极坐标放样的方位角及距离;由输入大概桩号及实测坐标，反算准确桩号及距中桩距离。(绿色为程序名，蓝色为输入计算器里的内容)
1、 主程序
程序名：1XY-KM
LbI 0：7→DimZ：Norm 2:1→A：＂A:X,Y=1，KM,D=2＂?A:
A=1=>Goto 1:A=2=>Goto 2:LbI 1:LbI A:Prog＂QX＂:90→B:＂PJ＂?B:B→Z[1]：LbI B:＂KM＂?Z:Z= -1=>Goto 0:Z=-2=>Goto A: Prog＂X1＂:?D: Prog "Z":I+D×COS(Z[2]+Z[1]→X: J+D×Sin(Z[2]+Z[1]→Y: Ｆix 3:＂   Ｘ=＂：Locate 6,4,Ｘ：＂   Ｙ＝＂：Locate 6,4,Ｙ：◢
Goto B:
LbI 2：Prog＂QX＂:＂KM＂?Z：LbI C: Prog＂X1＂: ＂XO＂?X:X=-1=>Goto 0:X=-2=>Goto 2: ＂Y0＂?Y:Prog ＂ZX＂: Fix 3:＂　ＫＭ＝＂：Locate 6,4,Ｚ：＂　　Ｄ＝＂：Locate 6,4,Ｄ◢
Goto C
2. 交点法正算子程序(Z)
程序名：Z
H2÷R÷24-H∧(4)÷2688÷R∧(3)→A:H÷2-H∧(3)÷240÷R2→B:（（H2-N2）÷24÷R÷Sin(Abs(P))-((H∧(4)-N∧(4))÷2688÷R∧(3))÷Sin(Abs(P→E:(R+A)tan(Abs(P)÷2)+B-E→T:P÷Abs(P→W:0→M:H→C: (R+N2÷R÷24-N∧(4)÷2688÷R∧(3))tan(Abs(P)÷2)+N÷2-N∧(3)÷240÷R2+E→Z[6]: O-T+ΠR×Abs(P)÷180+H÷2-N÷2→Z[7]:If Z≤O-T:Then Z-O→S:G→Z[2]:Goto 1: IfEnd: If Z≤O-T+H:Then Z-O+T→S:Prog “HX”:G+WK→Z[2]:A-T→A:G→E:Goto 2:IfEnd:If Z≤Z[7]: Then 180(Z-O+T-0.5H)÷R÷Π→S : A+R(1-Cos(S))→B:H÷2-H∧(3) ÷240÷R2+Rsin(S)→A:R→M:G+WS→Z[2]: A-T→A:G→E :Goto 2: IfEnd:
If Z≤Z[7]+N:Then Z[7]+N-Z→S:Z[6]→T :N→H:Prog “HX”:G+P→E:E-WK→Z[2]:T-A→A:C→H:Goto 2:IfEnd:
Z-Z[7]-N+Z[6]→S:G+P→Z[2]: LbI 1:U+Scos(Z[2])→I:V+Ssin(Z[2]→J: Return:LbI 2:U+Acos(E)-WBsin(E→I:V+ASin(E)+WBcos(E→J
3. 交点法缓和段子程序(HX)
程序名：HX
S-S∧(5)÷40÷R2÷H2+S∧(9)÷3456÷R∧(4) ÷H∧(4→A:S∧(3) ÷6÷R÷H-S∧(7)÷336÷R∧(3)÷H∧(3)+S∧(11)÷42240÷R∧(5) ÷H∧(5→B:90S2÷Π÷R÷H→K:RH÷S→M
4. 交点法反算子程序(ZX)
程序名：ZX
0→D:Do:Z+D→Z:Prog“Z”:Pol(X-I,Y-J+.000000001:J-Z[2]→J:Isin(J→S:Icos(J→D: If M≠0:Then Pol(M-WS,D:JMΠ÷180→D:IfEnd: LpWhile Abs(D)>0.001:Z+D→Z:S→D
5．线路选择子程序(线路选择输0时。则输曲线参数（QX）
程序名：QX
1→A:”QX:？=0，X1=1,X2=2”?A:A→Z[3]:
If A=0:Then ”JDKM”?O:”XJ”?U:”YJ”?V:”FW”?G:”LS1”?H:“LS2”？N：“JDZJ”？P:?R:IfEnd
6．坐标计算参数矩阵调出子程序（F）
程序名：F
Mat F[1,1]→O: Mat F[1,2]→U: Mat F[1,3]→V: Mat F[1,4]→G: Mat F[1,5]→H: Mat F[1,6]→N: Mat F[1,7]→P: Mat F[1,8]→R
7．数据子程序(根据自己标段设计数据修改红色部分)
程序名：X1（坐标计算要素程序,暂例子只给二条线路参数）

Z[3]=0=＞Return:Z[3]=1=＞Goto 1:Z[3]=2=＞Goto 2:
LbI 1:If  Z≤122269.856:Then [[121656.952, 2603848.019, 508351.147, 98°53′22.6″,660,300, 41°47′44.1″, 1415.097]]→Mat F: Prog”F”: Return: IfEnd:
If Z≤125115.9:Then [[123984.982, 2601974.464, 509885.441, 140°41′06.8″,700,400, －39°31′43.7″,2800]]→Mat F: Prog”F”: Return:IfEnd:
LbI 2:If  Z≤396.166:Then [[213.698, 2600569.493, 514256.323, 238°49′17.1″,107.584,120.474, 64°39′40″, 250]]→Mat F: Prog”F”: Return: IfEnd:
说明：If Z《交点HZ点桩号(或者直线段桩号)：Then[[交点桩号，交点X坐标，交点Y坐标，前交点至止交点方位角，交点前段缓和曲线长度，交点后段缓和曲线长度，交点转角，圆半径]] →Mat F: Prog”F”: Return: IfEnd:
注：转角分正负，左偏为负，右偏为正.加粗O为子母..Z[2]值为所计算点切线方位角。
输入部分：
1、＂A:X,Y=1，KM,D=2＂?  选择正反算，输1为正算，输2为反算；
3、”QX:？=0，X1=1，X2=2”? 选择手工输入参数及调入参数，输0则手工输入交点线元参数，输1则调入子程序X1中第一条线路，输入计算第二条线路。
4、PJ?  路线方向中线向右夹角，默认为90度；
5、JDKM？XJ？YJY？FW？LS1？LS2？JDZJ？R？分别输入：交点桩号，交点X坐标，交点Y坐标，前交点至止交点方位角，交点前段缓和曲线长度，交点后段缓和曲线长度，交点转角，圆半径
6、KM？ 输入计算段任意点桩号（主程序中输-1，返回选择正反算，输-2，返回选择输入选择线路或输入参数.
7、D？ 距中桩距离，左负，右正，中桩输0。
输出部分：
1、X=***    正算时，计算得出的所求点的X坐标
2、 Y=***    正算时，计算得出的所求点的Y坐标
5、ＫＭ=***    反算时，计算得出的所求点的里程
6、Ｄ=***    反算时，计算得出的所求点的边距
CASIO fx5800p全线高程计算程序
GAOCHEN  主程序

Lbl 1
“KM=，<0,Stop”:?K：K<0=>Stop:“PY=”？L：Prog”GK”
C-D→E:Abs(RE/2)→T：R(Abs(E)/E)→R
If　K≤B-T:Then 0→H:Else:If　K≥B+T
Then 0→H:D→C:Else K-B+T→H:Ifend:Ifend
A-(B-K)C-H2/(2R)→G：Cls
“KM=”：Locate 4，1，K：Locate  10，1，“PY=”：Locate 13，1，L：Fix 3
“H=”：Locate 4，2，G
Prog “PODU”：（E-B）/（D-A）（K-A）+B→I：(F-C)/（D-A）（K-A）+C→J
“HL=”：G+IL→X：Locate 4，3，X：Locate 11，3，“I=”：Locate 13，3，I*100
“HR=”：G+JL→Y：Locate 4，4，Y：Locate 11，4，“I=”：Locate 13，4，J*100◢  显示中边桩高程
Cls：Norm 2：“BM+HS≤0,Goto 1”？Z：Z≤0=> Goto 1：Cls    (输入视线高)
“KM=”：Locate 4，1，K：Locate  10，1，“PY=”：Locate 13，1，L：Fix 3
“QSM=”: Locate 6，2，Z-G                         (显示中桩读数)
“QSL=”: Locate 6，3，Z-X                          (显示左桩读数)
“QSR=”: Locate 6，4，Z-Y◢                         (显示右桩读数)
Norm 2：Cls：Goto1
（后面可加已知视线高计算读数部分，不想计算读数则视线高输入0或负数如不想显示麻烦，可将Locate语句去掉）
以下两个子程序不需运行，只是两个独立的数据库赋值程序，字母重复不影响计算结果
GK  数据库子程序

If K≤第二曲线起点桩号：Then 第一曲线交点高程→A：第一曲线交点桩号→B：第一曲线前坡→C：第一曲线后坡→D：第一曲线半径→R：Return：Ifend
…………….（有几个变坡点编几个If语句）

PODU  计算坡度子程序

If K≤第一变（非变）坡段终点：Then 第一曲线起点桩号→A：第一曲线起点左坡→B：第一曲线起点右坡→C：第一曲终点桩号→D：第一曲终点左坡→E：第一曲终点右坡→F： Return：Ifend
………………（每一个超高变化线元一个If语句）

结果显示：
高程显示
KM=0000.000  PY=0.000
H=  00.000     
HL=00.000    I=-1.5
HR=00.000    I=-1.5
读数显示
KM=0000.000  PY=0.000
QSM=  00.000
QSL=   00.000
QSR=  00.000