5800交点法正反算坐标(自己修改)
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,Ｙ：“FWJ”: Locate 6,4,Z[2]◢
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+Ssi n(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.44 1,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:
5800计算器程序下载
阅览室馆友我的图书馆 kaixin100 | | 分享
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+Zcos（F）→X：Y+Zsin（F）→Y 反算子程序SUB2 G-90→T (Y-V)cosT-(X-U)sin(T) →W Abs(W)→W:0→Z Lbl6：Prog "SUB1" T+QEW(C+WD) →L：(J-Y)cos（L）-(I-X)sin（L）→Z IF Abs（Z）<1E-6：Then0→Z:Prog "SUB1":(J-Y)÷sin（F）→Z:Else W+Z→W：Goto6:IfEnd 数据库子程序SUB0 Goto 1(线元可输入多条，分离式可在前多加一位，匝道一样。例：左幅为K129+500,右幅输线元参数里程为1129+500，其他不变，前面 1为任意数字，计算机便于区分) Lbl 1：IF S<线元终点里程：Then@@@→O(线元起点里程) ：@@@ →U(线元起点X坐标)：@@@→V(线元起点Y坐标)：@@@→G(线元起点计算方位角)：@@@→P(线元起点半径)：@@@→R(线元止点半径)：@@@→H(线元长度)：@@@→Q(线元左、右偏标志，左偏-1，右偏1，直线为0)：Return：IfEnd IF S<线元终点里程：Then@@@→O(线元起点里程) ：@@@ →U(线元起点X坐标)：@@@→V(线元起点Y坐标)：@@@→G线元(起点计算方位角)：@@@→P(线元起点半径)：@@@→R(线元止点半径)：@@@→H(线元长度)：@@@→Q(线元左、右偏标志，左偏-1，右偏1，直线为0)：Return：IfEnd 。。。。。。。。。。。。。。。。。。。。。。。。。。。 一程序功能 本程序由一个主程序(ZBJS)和3个子程——正算子程序(SUB1)、反 算子程序( SUB2) 、数据库子程序（SUB0）构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线 元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲 率半径、止点曲 率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标 进行正反算。另 外也可以将本程序中核心算法部分的两个子程序移植到其它相关的 程序中，用于对曲 线任意里程中边桩坐标进行正反算。本程序也可以在CASIO fx-4500P计算器及CASIO fx-4850P计算器上运行。 二、使用说明
FX5800P计算器坐标正反算程序Word文档
(以下程序是专业人士编写，本店铺不对程序负责，仅供您参考使用。) 卡西欧fx5800p计算器坐标正反算程序
一、程序功能 本程序由 6 个主程序、 5 个次子程序及 5 个参数子程序组成。主要用于公路测量中坐标正反算，设计任意点高程及横坡计算 , 桥涵放样，路基开挖口及填方坡脚线放样。程序坐标计算适应于任何线型 . 二、源程序 1. 主程序 1 ：一般放样反算程序(① 正算坐标、放样点至置仪点方位角及距离；② 反算桩号及距中距离 ) 程序名 :1ZD-XY Lb1 0:Norm 2 F=1 ： ( 正反算判别， F=1 正算， F=2 反算 , 也可以改 F 前加？，改 F 为变量 ) Z[1]=90 （与路线右边夹角） Prog ＂ THB ＂： F=1=>Goto 1:F=2=>Goto 2 Lb1 1: Ｆ ix 3: ＂Ｘ = ＂： Locate 6,4, Ｘ◢ ＂Ｙ＝＂： Locate 6,4, Ｙ◢ Ｐ rog ＂3JS”:Goto 0: Lb1 2:Fix 3: ＂ＫＭ＝＂： Locate 6,4, Ｚ◢ ＂Ｄ＝＂： Locate 6,4, Ｄ◢ Ｇ oto 0 ２．主程序２：高程序横坡程序 ( 设计任意点高程及横坡 ) 程序名： 2GC LbI 0:Norm 2 “KM”?Z:?D: Prog”H”:Fix 3:” H=”:Locate 6,4,H◢ “ I=”: Locate 6,4,I◢ Goto 0 3. 主程序 3 ：极坐放样计算程序 ( 计算放样点至置仪点方位角及距离 ) 程序名： 3JS X ： Y ： 1268 ．123→K( 置仪点 X 坐标 ) 2243 ．545→L （置仪点 Y 坐标，都是手工输入 , 也可以建导线点数据库子程序 , 个人认为太麻烦） Y-L→E ： X-K→F ： Pol(F,E):IF J<0:Then J+360→J:Int(J)+0.01Int(60Frac(J))+0.006Frac(60Frac(J)) →J:( 不习惯小数点后四位为角度显示的，也可以用命令J◢DMS◢ 来直接显示) Fix 4:” FWJ=”: Locate 6,4,J◢( 不习惯小数点后四位为角度显示的，也可以用命令 J◢DMS◢ 来直接显示 ) Fix 3:” S=”:Locate 6,4,I◢ 4 ．主程序 4 ：涵洞放样程序（由涵中心桩号计算出各涵角坐标、在主程序 3 中输入置仪点坐标后计算放样点至置仪点方位角及距离 ) 程序名： 4JH-XY LbI 0:Norm 2 90→Z[1]( 涵洞中心桩与右边夹角，手工输入，也可以修改成前面加？后变为变量 )
CASIO5800ZXZFS直线正反算坐标计算程序
[原创] CASIO fx-5800P ZXZFS 直线正反算坐标计算程序 原理：采用正余弦定理求待求点垂距进行确定里程和偏距。 公式：tan-1((yb-ya)/(xb-xa)) 、√(xB-xA)2+(yB-yA)2、sinα×D、Q÷tan(α) 说明：该程序适用于计算器 CASIO fx-5800P，针对直线正反算设计，正算可算直线中边桩坐标及其斜交、反算可推算直线距线路中心的里程偏距。 1、DK(QD)? 输入起点里程 2、X(QD)? 输入起点坐标X 3、Y(QD)? 输入起点坐标Y 4、FWJ? 输入计算方位角 5、J? 输入1为正算、输入2为反算 6、JSDK? 输入计算里程 主程序名：ZXZFS Deg：Fix 4↙ ＂DK(QD)＂?A：＂X(QD)＂？B：＂Y(QD)＂？C：＂FWJ＂？F↙ ＂1=>XY.2=>DK＂：＂J＂？W↙ If W=1：Then Goto 1：Else If W=2：Then Goto 2：IfEnd↙ Lb1 1↙ ＂JSDK＂？E：＂PJ＂?G：＂PA＂？H↙ ＂X=＂：B+(E-A) ×cos(F)+G×cos(F+H)→X◢ ＂Y=＂：C+(E-A) ×sin(F)+G×sin(F+H)→Y◢ ＂ZJD(X)＂？S：＂ZJD(Y)＂？T↙ ＂PJ=＂：Pol(X-S，Y-T)：I→M◢ ＂FWJ=＂：If J<0：Then J+360→N：Else J→N：N▲DMS◢ Goto 1↙ Lb1 2↙ ＂X=＂？K：＂Y=＂？L↙ Pol(K-B，L-C)：I→D：If J<0：Then J+360→P：Else J→P：IfEnd↙ If F M = (1.0/Re-1.0/Rs)/Ls; x=∫｛cos(Ta + L/Rs + 0.5*M *L*L)，0，L｝; y=∫｛sin(Ta +L/Rs + 0.5*M *L*L)，0，L｝; a(i)= Ta +L/Rs + 0.5*M *L*L Rs:缓和曲线起点半径 Re:缓和曲线止点半径 Rs，Re （NE坐标系下，右偏为正，左偏为负） Ta:缓和曲线起点的真北方位角 Ls:不完整缓和曲线长度。 此公式为缓和曲线在坐标系下任意位置的通用积分公式，能完全适应缓和曲线左偏、右偏、Rs >Re 、Rs CASIO fx-5800P线元法坐标正反算程序V2.0 说明：本程序适用于卡西欧计算器 CASIO fx-5800P,可对全线贯通坐标正反算、竖曲线高程计算。该程序可计算任意线型，包含（直线、圆曲线、缓和曲线、卵形曲线）等,还可以能通过坐标反推该点里程和距中线距离,适用测量员专用。 主程序名：ABCYT 第1步Deg：Fix 3：10→DimZ 第2步Lbl 3：＂1.DK=>XY＂：＂2.XY=>DK＂：＂Q＂?W：＂K0+＂?S：Prog＂ABCYTSJ＂：If P=0：Then 10^(45)→P：IfEnd：If R=0：Then 10^(45)→R：IfEnd 第3步1÷P→C：(P-R)÷(2HPR)→D：180÷π→E：If W=1：Then Goto 1：Else Goto2：IfEnd 第4步Lbl 1：＂E＂?Z 90=N：Abs(S-O)→W：Prog ＂ABCYTZ＂ 第5步Cls：＂F=＂：Locate 3,1,F°：＂X=＂：Locate 3,2,X：＂Y=＂：Locate 3,3,Y◢ 第6步Prog＂ABCYTSQX＂：Cls：＂H=＂：Locate 3,1,H◢ 第7步1→W：90→N：Goto 3 第8步Lbl 2：?X：?Y：X→I：Y→J：Prog＂ABCYTF＂：O+W→S 第9步Cls：＂K=＂：Locate 3,1,S：＂E=＂：Locate 3,2,Z◢ 第10步2→W：Goto 3 正算子程序名：ABCYTZ 第1步0.1739274226→A：0.3260725774→B：0.0694318442→K：0.3300094782→L 第2步1-L→F：1-K→M 第3步U+W×(A×cos(G+Q×E×K×W×(C+K×W×D))+B×cos(G+Q×E×L×W×(C+L×W× D))+B×cos(G+Q×E×F×W×(C+F×W×D))+A×cos(G+Q×E×M×W×(C+M×W×D)))→X 第4步V+W×(A×sin(G+Q×E×K×W×(C+K×W×D))+B×sin(G+Q×E×L×W×(C+L×W× D))+B×sin(G+Q×E×F×W×(C+F×W×D))+A×sin(G+Q×E×M×W×(C+M×W×D)))→Y 第5步G+Q×E×W×(C+W×D)→F：F+ N→Z[1] 第6步X+Z×cos（Z[1]）→X：Y+Z×sin（Z[1]）→Y 反算子程序名：ABCYTF 第1步Lbl 2：（S-O）→W：0→Z：Prog ＂ABCYTZ＂：F-90→Z[9]：(J-Y)×cos(Z[9])-(I-X)×sin(Z[9]) 道路中边桩坐标放样正反算CASIO fx-5800P程序（全线贯通） 编辑 | 删除 | 权限设置 | 更多▼ 设置置顶推荐日志转到私密记事本 转载自王中伟转载于2009年08月12日 17:34 阅读(1) 评论(0) 分类：技术交流权限: 公开 一、前言 本程序是《CASIO fx-5800P计算与道路坐标放样计算》中道路坐标放样计算程序的升级改进版本。原道路坐标放样计算程序只基于道路的单个基本型曲线，有效计算范围仅包括平曲线部分和前后的两条直线段，使用时需要输入平曲线设计参数，无坐标反算桩号功能。 改进后的程序名称为：道路中边桩坐标放样正反算程序（全线贯通），增加了可实现全线贯通的数据库功能和坐标反算桩号功能，主要是： 1．使用道路平面数据库子程序，可将一段或若干段道路的交点法格式平面参数（可容易从直线、曲线及转角表中获得）以数据库子程序形式输入计算器，程序在计算时省却了输入原始数据的麻烦； 2．坐标正算方面，输入桩号即可进行道路的中、边桩坐标计算，若输入了测站坐标，还可同时计算全站仪极坐标放样数据（拨角和平距）； 3．坐标反算方面，输入平面坐标，即可计算对应的桩号和距中距离（含左右信息）； 4．对于存在断链的道路，可分段分别编写数据库子程序，然后在主程序中添加一个路段选择的功能即可实现（可参照立交匝道程序中匝道的选择）。 程序的特点： 1．可进行中桩坐标的正、反算，程序代码简洁，便于阅读和改写； 2．主程序通过调用数据库子程序，省却了使用时输入平面参数的繁琐； 3．使用数据库子程序，换项目只需改写数据库子程序，程序通用性强。 二、道路示例项目基本资料 基本资料同《CASIO fx-5800P计算与道路坐标放样计算》第6章HY高速公路第2合同段（合同段起止桩号：K4+800~K9+600）。这里摘取直线、曲线及转角表资料如下（若图片不清晰，请参见参见教材P161附录1）： . 3、交点法、线元法坐标计算 坐标计算是根据图纸中“直线及曲线转角一览表”提供的数据计算道路中桩坐标，然后和图纸提供的“逐桩坐标表”比对，如果一样则说明输入平曲线参数输入正确，可以计算边桩坐标和其他结构物坐标了；如果中桩坐标不一样，一般是平曲线参数输入有误，需要重新检查输入，另一种结果是图纸有错，这种情况少见，但不代表没有。“直线及曲线转角一览表”和“逐桩坐标表”见附件1、附件2。 线元法是以路线的起点坐标、方位角、起终点桩号等节点元素来计算出要求的坐标；交点法是以路线的交点要素和路线的主要要素来求得坐标。 ①交点法 交点：路线的转折点，路线改变方向是相邻两直线的延长线相交的点。用JD表示， 有些图 纸上用 IP表示。 看下图： 交 点是针对曲线的（包含圆曲线和缓和曲线），一段曲线就有一个交点。交点参数有：坐标（X,Y）、交点桩号、转角值、圆曲线半径R、缓和曲线长度。 教学提供软件（轻松测量、双心软件、测量工具）交点法曲线要素输入说明： 1、QD起点坐标： 起点坐标必须在直线段上，或填写前一交点的坐标。 2、JD交点曲线要素： （1）交点桩号 （2）交点坐标（X,Y） （3）曲线半径R 始点的话，起始里程有时候需要校正，当然，并不是每个图纸给出的起点里程都需要校正，大多数图纸的起点里程已经被设计院校正过，我们输入平曲线的时候需要验证一下。如果我们按照图纸给出的起点里程输入，发现后面的交点里程都和图纸相差一个相同的值，这就表明我们输入的起点里程需要校正。 起始点里程正常输入，第二、三个交点输入完成后，检查第二个交点的切线长和交点 里程是否和图纸一样，如果切线长正确，交点里程不正确，说明起点里程需要校正，将第二个交点的里程与正确里程的差值，应用到起点里程中，从而使第二个交点里程和后面交点的里程与图纸吻合。 注意：交点法计算坐标适用的平曲线为对称或不对称缓和曲线、圆曲线。对于非普通的三单元曲线，交点法不适用。非普通的三单曲线例如下页的JD18及JD19处的平曲线， 的输入是否正确，有的图纸给的方位角数据较少，需要每隔几个线元才能检验方位角。 卡西欧5800公路坐标正反算程序 目录 一、坐标正算基本公式 (02) 二、坐标反算原理 (04) 三、高程数据库录入变换 (05) 四、计算器程序 (07) 01、ZBZS(坐标正算) (07) 02、ZBFS（坐标反算） (08) 03、GCJF(高程积分) (09) 04、PJFY（坡脚放样） (10) 05、JFCX（积分程序） (11) 06、ZBFY（坐标放样） (11) 07、DT（递推） (12) 08、HP（横坡） (13) 09、LK（路宽） (14) 10、SJK1（平面数据库） (14) 11、SJK2（纵面数据库） (14) 12、SJK3（左路宽度数据库） (15) 13、SJK4（右路宽度数据库） (15) 14、SJK5（横坡数据库） (16) 15、SJK6（下边坡数据库） (16) 16、SJK7（左上边坡数据库） (17) 17、SJK8（右上边坡数据库） (18) 五、后记 (19) CASIO 5800计算器公路工程测量程序 一、正算所涉及的计算公式 X R i d X αβB d Y d l d βI 图表 1 在图1中，A 点为回旋曲线起点，B 点为回旋曲线止点，I 点为所求坐标 点。设： A 点的X 坐标为X A ，Y 坐标为Y A ，A 点的切线方位角为α，A 点的曲率为 ρA ,A 点的里程为L A ,B 点的曲率为ρB ,B 点的里程为L B ,I 点的曲率为ρI ,I 点的 里程为L I 。I 点的切线角为β。 由于回旋线上各点曲率半径R i 和该点至曲线起点的距离L 成反比。故此任 意点的曲率为; C L R i i ==1ρ (c 为常数) （1） 由式(1)可知，回旋曲线任意点的曲率按线性变化，由此回旋曲线上里程为L i 点的曲率为; A B A i A B A i L L L L --?-+=)(ρρρρ （2） 当曲线右偏时ρB 、ρA 取正值，反之取负值。设： 本程序由一个主程序JD和三个子程序（JDA、JDB、JDC）构成，运行时只需运行主程序即可！ 本程序适用于单交点对称型、不对称型、无缓和曲线单圆曲线型一个交点范围内（含交点前后有直线段时）的曲线要素核对和坐标计算，手工输入要素，对设计图纸的“直线、曲线转角表”中交点数据进行复核验证，并为线元法程序提供起点坐标起点切线方位角等数据！当然本程序也可单独逐交点输入进行放样计算用！鉴于5800计算器的空间和以上所述本程序的主要目的，故此程序不修改为数据库版本！需要的自行修改结合XY框架自己修改为数据库反算程序等！ 主程序名：JD 24→Dimz↙ Cls :＂XC＂?U :＂YC＂?V :＂K（JD）＂?K :＂X（JD）＂?X : ＂Y（JD）＂?Y :＂LS1＂?B :＂LS2＂?C : ?R : ＂（ZH）FWJ°＂?M : ＂α（Z-,Y+）°＂?O : M+O→N : Prog ＂JDA＂↙ Cls :＂T1=＂:＂T2=＂:＂L=＂:＂LY=＂: Locate 4,1,S : Locate 4,2,T : Locate 4,3,L : Locate 4,4,Q◢ Cls :＂E=＂:＂K（ZH）=＂: Locate 7,1,E : Locate 7,2,Z[1]◢ Cls : ＂K（HY）=＂:＂K（QZ）=＂:＂K（YH）=＂:＂K（HZ）=＂: Locate 7,1, Z[2] : Locate 7,2, Z[3] : Locate 7,3, Z[4] : Locate 7,4, Z[5]◢ LbI 0 : ＂K×+×××＂?P : ＂Z＂?D : If D≠0 :Then ＂RJ＂?H : IfEnd : Prog ＂JDB＂↙ If D＜0 :Then Cls : ＂X(L)=＂:＂Y(L)=＂: Locate 6,1,F : Locate 6,2,G◢Pol(F-U,G-V : Cls : ＂S(L)=＂: Locate 6,1,I : ＂F(L)=＂: 360Frac((J+360)÷360▼DMS◢ Goto 0 : IfEnd↙ If D=0 :Then Cls : ＂X(Z)=＂:＂Y(Z)=＂: Locate 6,1,F : Locate 6,2,G : ＂QXFWJ(Z)=＂: Z▼DMS◢ Pol(F-U,G-V : Cls : ＂S(Z)=＂: Locate 6,1,I : ＂F(Z)=＂: 360Frac((J+360)÷360▼DMS◢ Goto 0 : IfEnd↙ If D＞0 :Then Cls : ＂X(R)=＂:＂Y(R)=＂: Locate 6,1,F : Locate 6,2,G◢Pol(F-U,G-V : Cls : ＂S(R)=＂: Locate 6,1,I : ＂F(R)=＂: 360Frac((J+360)÷360▼DMS◢ Goto 0 : IfEnd↙ 子程序1名： JDA If O＜0 :Then -1→W : Else 1→W : IfEnd : WO→A ↙ B2 ÷24÷R-B^(4)÷2688÷R ^(3) →Z[6] ↙ C2 ÷24÷R-C^(4)÷2688÷R ^(3) →Z[7] ↙ B÷2-B^(3)÷240÷R2 →Z[8] ↙ C÷2-C^(3)÷240÷R2 →Z[9] ↙ Z[8]+((R+Z[7]-(R+Z[6])cos(A))÷sin(A))→S↙ Z[9]+((R+Z[6]-(R+Z[7])cos(A))÷sin(A))→T↙ RAπ÷180+(B+C) ÷2→L↙ RAπ÷180-(B+C) ÷2→Q↙ 颜色文字Casio5800交点法与线元法（积木法）匝道坐标正反算放样程序（XUFENG 2011.2.14） 本人一直以来想找一个交点法与线元法相结合的坐标正反算程序，在网上找了很久很久，没能找到一个较为满意的，有幸在测量空间看到大歪哥的《Casio5800交点法程序》与《线元法（积木法）匝道坐标正反算放样程序》，根据歪哥意见“需要的自行修改结合XY框架自己修改为数据库反算程序等”，本人不才，采用最笨的办法将两个程序综合了一下，使之能既能进行交点法正反算，又能进行线元法正反算。在此特别感谢大歪哥！ 将程序发上来，愿与大家一同交流学习欢迎大家吐口水，只要能进步就行！ 程序由一个主程序ZBZFS和8个子程序（JS、XY-A、XY-B、JDYS、1、2、3、4）构成，运行时只需运行主程序即可！ 本程序适用于单交点对称型、不对称型、无缓和曲线单圆曲线型一个交点范围内（含交点前后有直线段时）的曲线要素核对和坐标正反算，手工输入要素，对设计图纸的“直线、曲线转角表”中交点数据进行复核验证，并能对单一线元进行坐标正反算。 1主程序名：ZBZFS（功能：进入计算主程序） 65→Dimz↙ Deg:Fix 3↙ ＂1.JD ZFS 2. ZHADAO ZFS＂? I: I→Z[61]: ＂1.ZHONG SHU JS 2. JS＂? I ↙ If I=1: Then Goto1: Else Goto2:IfEnd↙ LbI 1 :If Z[61]=1: Then Prog＂JDYS＂:Else Cls:＂K0＂?A:＂KN＂?L :＂X0＂?U :＂Y0＂?V :＂F0＂?W :＂R0＂?P :＂RN＂?Q:＂ZX:-1,+1,0＂?G:IfEnd↙ LbI 2 :Prog＂JS＂ 2子程序名：JS（功能：选择正算或反算模式） Cls:＂XC＂?H:＂YC＂?Z↙ 5800全线任意坐标计算程序 1. 正算主程序(ZHCX) （不运行） 8→DimZ 1÷P→Z[4 ]：(P-R)÷(2HPR)→D: 180÷π→E “Z=”?Z:”YJJ=”?A:Abs(S-O)→W 0.26→Z[1 ]: 0.74→B: 0.02→K: 0.82→Z[3 ]: 1-Z[3 ]→F:1-K→Z[2 ] U+W(Z[1 ]cos(G+QEKW(Z[4 ]+KWD))+Bcos(G+Z[3 ]QEW(Z[4 ]+ Z[3 ]WD))+Bcos(G+QEFW (Z[4 ]+FWD))+ Z[1 ]cos(G+ Z[2 ]QEW(Z[4 ]+ Z[2 ]WD)))→X: V+W(Z[1 ] sin (G+QEKW(Z[4 ]+KWD))+B sin(G+ Z[3 ]QEW(Z[4 ]+ Z[3 ]WD))+B sin(G+QEFW (Z[4 ]+FWD))+ Z[1 ] sin(G+ Z[2 ]QEW(Z[4 ]+ Z[2 ]WD)))→Y： G+QEW(Z[4 ]+WD)→F：X+Zcos(F+A)→X：Y+Zsin(F+A)→Y:If F≧360:Then F-360→F:IfEnd ”X=”:X→X◢ ”Y=”:Y→Y◢ If F﹤0:Then F+360→F:IfEnd ”QX FWJ=”:F▼DMS◢ “C=1=>XX: C=2=>XZ”: ”C=”?C: ”QHJU=”?L: If C=1：Then Goto 1:Else Goto 2: IfEnd 可以计算斜交斜做或斜交正做的桥涵坐标 Lbi 1 X+L cos(F)→X：Y+Lsin(F)→Y: Goto 3 Lbi 2 X+L cos(F+A-90)→X：Y+Lsin(F+A-90)→Y: Goto 3 Lbi 3 “QH-X=”: X →X◢ “QH-Y=”: Y →Y◢ Prog “FY” 2 . 参数子程序（直接运行） M(主线) 一条线路一个名称 “S=”?S If S≦线元终点：Then 线元起点X值→U: 线元起点Y值→V:线元起点切线方位角→G:线元起点桩号→O:线元长度→H：线元起点半径→P：线元终点半径→R：（左偏-1，或右偏 1）→Q:Goto 1:IfEnd … … If S≦线元终点：Then 线元起点X值→U: 线元起点Y值→V:线元起点切线方位角→G:线元起点桩号→O:线元长度→H：线元起点半径→P：线元终点半径→R：（左偏-1，或右偏 1）→Q:Goto 1:IfEnd Lbi 1 Prog “ZBJS” 3. 放样程序（FY）（不运行） “X0=”?M:“Y0=”？N Pol((X-M, Y-N) Casiofx-5800P交点法坐标计算程序 XC U 测站X坐标 YC V 测站Y坐标 K(JD) K 交点桩号 X(JD) X 交点X坐标 Y(JD) Y 交点Y坐标 LS1 B 第一缓和曲线长度 LS2 C 第二缓和曲线长度 R 圆曲线半径 FWJ M 起始边切线方位角 α(Z-,Y+) O 本交点处线路转角(左转为负，右转为正，度分秒输入) K×+××× 待求桩号 Z 待求桩号距中距离(左负值，右正值，中为0) RJ 斜交右角(线路切线前进方向与边桩右侧夹角) →? 主程序名：JD 24→Dimz ”XC”?U:”YC”?V:”K(JD)”?K:”X(JD)”?X:“Y(JD)”?Y:”LS1”?B:”LS2”?C:?R “FWJ”?M:“α(Z-,Y+)°”?O:M+O→N Prog“JDA” Cls:”E=“:”K(ZH)=“:Locate7,1,E:Locate7,2,Z[1]◢ Cls:“K(HY)=“:”K(QZ)=“:”K(YH)=“:”K(HZ)=“:Locate7,1,Z[2]:Locate7,2,Z[3]:Locate7,3,Z[4]:Locate7,4,Z[5]◢LbI0:“K×+×××”?P:“Z”?D:If D≠0:Then “RJ”?H:IfEnd:Prog“JDB” If D<0:Then Cls:“X(L)=“:”Y(L)=“:Locate6,1,F:Locate6,2,G◢ Pol(F-U,G-V:Cls:“S(L)=“:Locate6,1,I:“F(L)=“: 360Frac((J+360)÷360?DMS◢ Goto0:IfEnd If D=0:Then Cls:“X(Z)=“:”Y(Z)=“:Locate6,1,F:Locate6,2,G:“QXFWJ(Z)=“:Z?DMS◢ Pol(F-U,G-V:Cls:“S(Z)=“:Locate6,1,I:“F(Z)=“: 360Frac((J+360)÷360?DMS◢ Goto0:IfEnd If D>0:Then Cls:“X(R)=“:”Y(R)=“:Locate6,1,F:Locate6,2,G◢ Pol(F-U,G-V:Cls:“S(R)=“:Locate6,1,I:“F(R)=“: 360Frac((J+360)÷360?DMS◢ Goto0:IfEnd 子程序1名：JDA If O<0:Then -1→W:Else 1→W:IfEnd:WO→A B2÷24÷R-B^(4)÷2688÷R^(3)→Z[6] C2÷24÷R-C^(4)÷2688÷R^(3)→Z[7] B÷2-B^(3)÷240÷R2→Z[8] C÷2-C^(3)÷240÷R2→Z[9] Z[8]+((R+Z[7]-(R+Z[6])cos(A))÷sin(A))→S Z[9]+((R+Z[6]-(R+Z[7])cos(A))÷sin(A))→T RAπ÷180+(B+C)÷2→L RAπ÷180-(B+C)÷2→Q (R+(Z[6]+Z[7])÷2)÷cos(A÷2)-R→E K-S→Z[1] 坐标正反算程序 程序变量说明: D:选线。E、F:置仪点坐X、Y。K:待算点桩号。B:待算点距离。Q:待算点与线路交角（以顺时针方向）L:线元起点至待算点长度。0:待算点方位角。S:线元长度。U:线元起点半径。V:线元终点半径。 X、Y:正算结果。Z、计算序列。C、I、J T:内部运算。A、G、H、M、N、P、R、W未用变量。 ZB- Deg:Norm 2 Z[1]T Z:”两个空格ZS=1 Or FS=2Z If Z M 1 And Z M2:Then Cls:Locate 5,2,”NO 10r 2”：Locate 6,3 !Stop!”： Stop:Else Z Z[1]:lfEnd Z[10]T D:四个空格XL T 1 T 2T3九个空格XL N=N”D If D M 1 And D M 2 And D M 3:Then Cls:Locate 5,2,'XL r 1 T2T3”:Locate 5,3 !Stop!”： Stop: Else D T Z[9]: IfEnd If D=Z[10]:The n Cls:Do:Locate 3,1,"[EXE] GO!GO!':Locate 8,2,9-I ntg(Z/7) Locate 4,3,"10S T DAT":Locate 12,3,D:Isz Z:Z=70=>Goto 1:Lpwhile Getkey M 47:Goto 2:IfE nd Lbl 1 D=1=>Prog” ”:D=2=>Prog” ”:D=3=>Prog” ”:D=X=>Prog‘ D T Z[7]:For 1 T Z To Z[7]-1 Cls: Locate 9,2,Z[7]-Z Z[10Z+11]-Z[10Z+1]T L Prog" TQSJ:Next Lbl 2 Z[9] T Z[10] Cls:Z[1]=2=>Prog”FS' Z[2]T E: (X)=”E:E T Z[2]:Z[3] T F：”(Y)仝F:F T Z[3]:Prog ”ZS' FX5800计算器测量程序集版 一、程序功能 主要功能：采用交点法方式计算多条线路坐标正反算，可算任意复杂线型及立交匝道，包括C型，S型、卵型、回头曲线等；极坐标放样，全线路基边坡开挖口及坡脚放样计算、路基任意点设计高程、横坡、设计半幅宽度.隧道欠超挖放样计算。 新版本优化：1、优化程序语句、2、以复数形式输入变量及做数据库，取消原矩阵数据库;3、修改隧道超欠挖程序为通用形，不受圆心个数限制、4、新增测量资料表计算 二、源程序（绿色为程序名；蓝色为输入计算器内容）紫色为新版改动处（可以根据自己标段情况用相关主程序及子程序，再在0程序中汇总）0.汇总程序（1、坐标计算放样程序(1XY、A、AB、HX、JS、DX、QX、F、XY、X1)；2、坐标反算程序（2ZD、A、B、AB、HX、QX、F、ZD、X1)；3、高程计算查阅程序（3GC、H、I、QX、S1、I1）；4、路基半幅标准宽度查阅程序（4GD、C、QX、G1）；5、路基边坡及开挖口放样程序（5BP、 A、B、AB、HX、H、I、C、JS、DX、QX、F、ZD、X1、S1、I1、G1、W1）； 6、路基标准距离放样（6FM、A、AB、HX、H、I、JS、DX、QX、F、XY、X1、S1、I1)； 7、桥梁锥坡计算放样程序（7ZP、A、AB、HX、C、JS、DX、QX、F、XY、X1、G1)； 8、极坐标计算程序（8JS、JS、DS）； 9、隧道超欠挖计算程序（A、B、AB、HX、H、I、QX、S、SD、F、ZD、X1、S1、I1、SD1）运行后按1~9数子约半秒，则选择1至9的程序，返回时，在桩号输入-1,返回选择计算类型。输入-2,返回选择线路。 程序名：0（数子0） ClrMat:ClrVar:12→DimZ：Norm 2:Do:＂(XY=1，ZD=2 ，GC=3，GD=4，BP=5，FM=6，ZP=7，JS=8，SD=9)===>QING AN 1-9＂：Getkey→Z[3]:While Z[3]=35:Prog＂1XY＂:WhileEnd:While Z[3]=36:Prog＂2ZD＂：WhileEnd: While Z[3]=37:Prog＂3GC＂：WhileEnd: While Z[3]=21:Prog＂4GD＂：WhileEnd: While Z[3]=22:Prog＂5BP＂：WhileEnd: While Z[3]=23:Prog＂6FM＂：WhileEnd: While Z[3]=31:Prog＂7ZP＂：WhileEnd: While Z[3]=32:Prog＂8JS＂： 一．使用说明： 1.规定： （1）以道路中线的前进方向区分线路的左右偏转方向，线路左偏Q=-1，线路右偏Q=1，直线Q=0。 （2）当所求点位于线路中线时，Z=0；当位于左侧时，Z取负值；当位于右侧时，Z取正值。 （3）当线元为直线时，其起终点的曲率半径为无穷大，取10的45次方代替。 （4）当线元为圆曲线时，起终点曲率半径为圆曲线半径。 （5）当线元为缓和曲线时，起终点曲率半径与直线相接取10的45次方代替，与圆曲线相接取圆曲线半径。 2.输入与显示说明： （1）输入部分： N？：选择计算方式，输入1表示由里程边距计算坐标，输入2表示由 坐标计算里程边距。 X0？线元起点X坐标。 Y0？线元起点Y坐标。 S0？线元起点里程。 F0？线元起点切线方位角。 LS？线元长度。 R0？线元起点曲率半径。 RN?线元终点曲率半径。 Q？线元左右偏标志（左偏=-1，右偏=1，直线=0） S?正算时所求点里程。 Z？正算时所求点距中线边距。位于线路左侧输入负值，在右侧输入正 值。 AR?斜交角度。（为线路右角，在0-180度之间） X?反算时所求点X坐标。 Y？反算时所求点Y坐标。 （2）显示部分： XS=计算结果，所求点X坐标。 YS=计算结果，所求点Y坐标。 FS=计算结果，所求点切线方位角。 S=计算结果，反算时所求点里程 Z=计算结果，反算时所求点边距。 二．程序清单1（用于任意曲线元的坐标计算） 1.主程序(TYQXJS) ″1.SZ=>XY″:″2.XY=>SZ″:?N:″X0″?U:″Y0″?V:″S0″?O: ″F0″?G:″LS″?H:″ R0″?P:″RN″?R:?Q:1÷P→C:(P-R) ÷(2HPR) →D:180÷π→E:If N=1: THEN Goto 1:Else Goto 2:IfEnd:Lbl 1:?S:”D(Z-，Y+)=”?Z:90→T:”O=”?T:Abs(S-O) →W:Prog ″SUB1″:″XS″:X◢″YS″:Y◢″FS″:F◢ Goto 1:Lbl 2:?X:?Y:I→X:J→Y:Prog ″SUB2″:″S″:O+W→S:◢″Z″:Z◢Goto 2 2.正算子程序(SUB1) 0.1739274226→A:0.3260725774→B:0.0694318442→K:0.3300094782→L:1-L→F:1-K→ 目录 一、坐标正算基本公式 (02) 二、坐标反算原理 (04) 三、高程数据库录入变换 (05) 四、计算器程序 (07) 01、ZBZS(坐标正算) (07) 02、ZBFS（坐标反算） (08) 03、GCJF(高程积分) (09) 04、PJFY（坡脚放样） (10) 05、JFCX（积分程序） (11) 06、ZBFY（坐标放样） (11) 07、DT（递推） (12) 08、HP（横坡） (13) 09、LK（路宽） (14) 10、SJK1（平面数据库） (14) 11、SJK2（纵面数据库） (14) 12、SJK3（左路宽度数据库） (15) 13、SJK4（右路宽度数据库） (15) 14、SJK5（横坡数据库） (16) 15、SJK6（下边坡数据库） (16) 16、SJK7（左上边坡数据库） (17) 17、SJK8（右上边坡数据库） (18) 五、后记 (19) CASIO 5800计算器公路工程测量程序 一、正算所涉及的计算公式 X R i d X αβB d Y d l d βI 图表 1 在图1中，A 点为回旋曲线起点，B 点为回旋曲线止点，I 点为所求坐标点。设： A 点的X 坐标为X A ，Y 坐标为Y A ，A 点的切线方位角为α，A 点的曲率为ρ A ,A 点的里程为L A , B 点的曲率为ρB ,B 点的里程为L B ,I 点的曲率为ρI ,I 点的里 程为L I 。I 点的切线角为β。 由于回旋线上各点曲率半径R i 和该点至曲线起点的距离L 成反比。故此任 意点的曲率为; C L R i i ==1ρ (c 为常数) （1） 由式(1)可知，回旋曲线任意点的曲率按线性变化，由此回旋曲线上里程为L i 点的曲率为; A B A i A B A i L L L L --?-+=)(ρρρρ （2） 当曲线右偏时ρB 、ρA 取正值，反之取负值。设： 5800计算器全线坐标计算放样程序（修改第三版） 5800计算器全线坐标计算放样程序（修改版） “XLZBJSCX” ◢ LB1 0 ↙ CLS : FIX 4 : 30→DIM Z ↙ “XHS="?G(后视点X):"YHS="?L(后视点Y):"XZJ="?M(置镜点X):"YZJ="?N(置镜点Y)0l(G-M,L-N):"DH=":I(后视距)◢J<0=>J+360→J:"FH=":J→DMS◢(后视 方位角) LB1 1 ↙ “K=”?K ◢（计算里程） IF K<本段曲线终点里程 AND K≥上段曲线终点里程：THEN 本段终点里程→Z[1] : 上段曲线终点里程→Z[2] ：1→0 （注：左偏曲线输入-1→0,右偏曲线输入1→0）: 偏角→A：半径→R : 第一缓和曲线→Z[6] : 第二缓和曲线 →Z[7] : 交点X→B :交点Y→C : 小里程向交点方位角→E : 交点向大里程方位角→F : G0T0 2 : IFEND↙ …………（曲线段分段输入） 补充直线段输入如下 IF K<本段直线终点里程 AND K≥本段直线起点里程:THEN 1→0:本段直线终点里程→Z[2]:终点X→Z[16]:终点Y→Z[11]:方位角→E:G0T0 4:IFEND LB1 2 ↙（曲线要素计算） Z[6]/2- Z[6]^3/(240*R^2)+ Z[6]^5/(34560*R^4) →Z[8] ↙（M1） Z[7]/2- Z[6]^3/(240*R^2)+ Z[7]^5/(34560*R^4) →Z[9] ↙（M2） Z[6]^2/(24*R)- Z[6]^4/(2688*R^3) →Z[10] ↙（P1） Z[7]^2/(24*R)- Z[7]^4/(2688*R^3) →Z[11] ↙（P2） π*A*R/180+0.5*( Z[3]+ Z[2])→W ↙（曲线总长） 90* Z[6]/(R*π) →Z[14] ↙（第一缓和曲线总偏角） 90* Z[7]/(R*π) →Z[15] ↙（第二缓和曲线总偏角,可以省略） Z[8]＋（R+Z[10])TAN(A/2)-(Z[11]-Z[11] )/SIN A→Z[11]↙ (切线T1) Z[9]＋（R+Z[12])TAN(A/2)+(Z[10]-Z[12] )/SIN A→Z[12]↙ (切线T2) B+ Z[12]*C0S (E+180)→ Z[13] ↙（ZH点X） C+ Z[12]*SIN(E+180)→ Z[15] ↙（ZH点Y） Z[1]-S→Z[3] ↙ (ZH点里程) Z[3]+ Z[6]→Z[4] ↙ (HY点里程) Z[1]- Z[7]→Z[5] ↙ (YH点里程) G0T0 3 ↙ LB1 3 ↙(判断里程点与曲线关系) IF K≤Z[3] AND K> Z[2] : THEN G0T0 4 : IFEND ↙ IF K≤Z[4] AND K> Z[3] : THEN G0T0 5 : IFEND ↙ IF K≤Z[5] AND K> Z[4] : THEN G0T0 6 : IFEND ↙ 5800P编程计算器坐标正反算程序 一，主程序文件名（TYQXJS） LbI 4:“1.SZ=> XY’’:’’2.XY=>SZ’’︰?N:?S:Prog“SUB0” 1÷P→C：(P－R) ÷(2HPR) →D：180÷π→E：N=1=> Goto1:Goto2 LbI 1: ？Z:Abs(s－O)→W︰Prog“SUB1”︰“XS=”:X?“YS=”:Y? F－90→F:”FS=”:F?DMS? Goto 4 LbI 2:?X:?Y:X→I:Y→J:Prog”SUB2”:O＋W→S﹕“S=”S?“Z=”：Z? Goto 4 二:子程序（文件名SUB1） 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＋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＋QEMW(C＋QEMW(C＋MWD)))→Y:G＋QEW(C＋WD)＋90→F:X＋Zcos(F)→X:Y＋Zsin(F)→Y 三:反算子程序(文件名SUB2) G－90→T:Abs((Y－V)Cos(T)－(X－U)Sin(T)→W:0→Z:LbI 0:Prog“SUB1”:T＋QEW(C＋WD)→L:(J－Y)cos(L)－(I－X)sin(L)→Z:IfAbs(Z)＜0.000001:ThenGoto1:Else W＋Z→W:Goto0 LbI 1:0→Z:Prog“SUB1”:(J－Y)÷sin(F)→Z: 四:数据库(文件名) LbI 1:If S＜止点里程:Then起点方位角→G:起点里程→O:起点X坐标→X:起点Y坐标→Y:起点半径→P:止点半径→R:曲线长→H:左转(-1)右转(1)直线(0)→Q:Return:IfEnd 数据库说明: 每增加一段曲线就增加以上一句上面的语句，当半径为无穷大时以10的45次方代替，左转输入负1，右转输入正1 五:程序输入及显示说明 1:输入部分 1，SZ=>XY 2，XY=>SZ N? 选择计算方式输入1表示由里程边距计算坐标(正算)输入2表示由坐标计算里程及边距(反算) S? 正算时输入所求点的里程，反算时为所求点的近似里程(注意反算时的近似里程不能跨曲线，如果跨曲线需要进行两次反算，并以第一的里程结果作为近似里程，否则，反卡西欧5800计算器坐标正反算程序
5800曲线道路全线坐标正反算(已验证)
FX5800道路路线测量程序
交点法线元法坐标计算
卡西欧5800公路坐标正反算程序教程文件
交点法坐标计算
5800坐标程序
5800简单全线坐标计算程序
Casio fx-5800P 交点法坐标计算程序
5800公路坐标正反算程序
Fx5800计算器公路测量程序设计
卡西欧5800测量坐标计算程序
卡西欧5800公路坐标正反算程序
5800计算器全线坐标计算放样程序(修改第三版)
5800P编程计算器坐标正反算程序