CASIO5800坐标正反算计算程序
CASIO 9860/5800坐标正反算计算程序
1. 坐标正算主程序(命名为ZBZS)
第1行：Lbl 0:”K=”?K:”BIAN=”? Z:”α=”?B
第2行：Prog “A”
第3行：”X=”:N+Zcos(F+B)◢
第4行：”Y=”:E+Zsin(F+B)◢
第5行：”F=”:F◢
第6行：Goto 0
K——计算点的里程
BIAN——计算点到中桩的距离（左负右正）
α——取前右夹角为正
2. 坐标反算桩号和偏距主程序(命名为ZBFS)
第1行：”X1=”? C:”Y1=”?D: ”K1=”?K
第2行：Lbl 0:Prog “A”
第3行：Pol(C-N,D-E)
第4行：List Ans[1]→I
第5行：List Ans[2]→J
第6行：Icos(F-J)→S:K+S→K
第7行：Abs(S)>0.0001=>Goto 0
第8行：”K1=”:K◢
第9行：”BIAN=”:Isin(J-F)→Z◢
X1——取样点的X坐标
Y1——取样点的Y坐标
K1——输入时为计算起始点(在线路内即可)，输出时为反算点的桩号
Z——偏距（左负右正）
3. 计算坐标子程序(命名为XYF)
为了简洁，本程序由数据库直接调用，上述中的正反算主程序不直接调用此程序第1行：K-A→S：(Q-P)÷L→I
第2行：N+∫(cos(F+X(2P+XI)×90÷π),0,S)→N
第3行：E+∫(sin(F+X(2P+XI)×90÷π),0,S)→E
第4行：F+S(2P+SI)×90÷π→F
4. 数据库（命名为A）
第1行：K≤175.191=>Stop
第2行：175.191→A:428513.730→N:557954.037→E:92°26′40″→F:0→P:1/240→Q:70.417→L:
K≤A+L =>GoTo 1（第一缓和曲线，圆半径为240）
第3行：245.607→A: 428507.298→N:558024.092→E: 100°50′59.4″→F: 1/240→P:1/240→Q:
72.915→L: K≤A+L =>Goto 1（第圆曲线，半径为240）
第4行：318.522→A: 428482.988→N:558092.538→E: 118°15′25.2″→F: 1/240→P: 0→Q: 55.104
→L: K≤A+L =>Goto 1（第二缓和曲线，圆半径为240）
第5行：373.627→A:428453.283→N:558138.912→E:124°50′4.5″→F:0→P:-1/180→Q:67.222→L:
K≤A+L=>Goto 1:Stop（第一缓和曲线，圆半径为180）
第6行：Lbl 1:Prog “XYF”
A——曲线段起点的里程
N——曲线段起点的x坐标
E——曲线段起点的y坐标
F——曲线段起点的坐标方位角
P——曲线段起点的曲率(左负右正)
Q——曲线段终点的曲率(左负右正)
L——曲线段长度（尽量使用长度，为计算断链方便）
说明：
（1）在9860中，程序中所有公式和部分函数结果均存储在List Ans列表数组中，要想多次调用最好随公式取出结果，并赋给变量。
（2）正算主程序可以计算一般边桩的坐标，如要计算类似涵洞端墙的坐标需增加第二偏距和转角两个变量。
（3）程序规定，左偏曲线曲率（半径倒数）输入负值，右偏曲线曲率输入正值，直线上点曲率输入0，例如直线段，线元起点和终点均输入0，第一缓和曲线分别输入0和圆半径的倒数，圆曲线均输入半径倒数，第二缓和曲线分别输入圆半径倒数和0，卵形曲线分别输入对应圆半径的倒数
（4）若是从大里程向小里程的反方向计算，则曲率取正方向时的负值，方位角减去(或加上)180度。
（5）有多个匝道的项目，可随时更改正反算主程序中的红色字体部分来调用其它线路的数据
（6）反算桩号偏差为1mm
（7）可以计算任意线型的任意点坐标
5. 坐标计算通用程序(命名为ZB)
本程序与上述中的几个程序无关，是在计算坐标时采用输入曲线元参数模式下的程序，为临时使用方便。
第1行：”K0=“? A:”X0=“? N:”Y0=“?E:”F0=“?F:”1÷R1=“?P:”1÷R2=“?Q”:“L0=“?L
第2行：Lbl 0:”K=”?K:”BIAN=”? Z:”α=”?B
第3行：Abs(K-A)→S：(Q-P)÷L→I
第4行：N+∫(cos(F+X(2P+XI)×90÷π),0,S)→X
第5行：E+∫(sin(F+X(2P+XI)×90÷π),0,S)→Y
第6行：”F=”:F+S(2P+SI)×90÷π◢（可省掉）
第7行：”X=”:X+Zcos(Z+B)◢
第8行：”Y=”:Y+Zsin(Z+B)◢
第9行：Goto 0
说明：本程序测试于卡西欧9860，适用于卡西欧计算器5800以上版本，去掉反算程序中第4、5行后可用在5800计算器上。
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 一、 Lbl 3:"1→ZS,2→FS"?Q Q=1=>Goto 1:Q=2=>Goto 2 Lbl 1:"CE:X"?M:"CE:Y"?F:"JL"?L:"FWJ"?A:Rec(L,A):M+I→C:F+J→D Cls "X=":Locate 3,1,C:"Y=":Locate 3,2,D◢ Goto 3 Lbl 2::"CE:X"?G:"CE:Y"?H:"(HOU)FY:X"?N:"(HOU)FY:Y"?E Pol(N-G,E-H) If J<0:Then J+360→Y:Else J→Y:IfEnd Cls "FY JL=":Locate 10,1,I:"FY FWJ=":Y◆DMS◢ Goto 3 进入程序运行如下: 1→ZS,2→FS? 输入1为正算,2为反算. 以输入1为例: CE:X? 测站点X(5796.717) CE:Y? 测站点Y(5212.569) JL? 仪器测得的距离(321.889) FWJ? 仪器测得的方位角(193-41-07) 得到:X=5483.966 Y=5136.414 再按EXE,输2为例: CE:X? 测站点X(5796.717) CE:Y? 测站点Y(5212.569) (HOU)FY:X? 后视或放样的X(5483.966) (HOU)FY:Y? 后视或放样的Y(5136.414) 得到:FY JL=321.889 FY FWJ=193-41-6.79 二、 Deg : Fix 3 : “XZ→0:YZ→1”？A : If A = 1: Then Goto 1 : IfEnd ↙ If A = 0 : Then “BS→0:XY→1:AND→2:DK→3:L(I)→4 ”？O : IfEnd ↙ If O = 4: Then Goto 1 : IfEnd ↙ If O = 3: Then Prog “F.2 ”: If X= 0 : Then Goto 1 : IfEnd : IfEnd ↙ If O≠1: Then “X1 ”？X : “Y1”？Y : X→Z［11］: Y→Z［12］: “X2 ”？P : “Y2”？Q : Pol( P－X , Q－Y) : If J﹤0 : Then J + 360→J : IfEnd : Cls : “S12= ”: Locate 6 ,1, I : “B12= ”: J ?DMS◣ 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程序（全线贯通） 编辑 | 删除 | 权限设置 | 更多▼ 设置置顶推荐日志转到私密记事本 转载自王中伟转载于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）： . 坐标正反算公式 一、GPS数据处理相关术语 1、三维无约束平差 三维无约束平差是以基线解算所得到的三维静态基线向量为观测值，待定参数主要为GPS 网中点的坐标；同时，利用基线解算时随基线向量一同输出的基线向量的方差阵，形成平差的随机模型，最终形成平差完整的数学模型。随后对所形成的数学模型进行求解，根据平差结果来确定观测值中是否存在粗差，数学模型是否有需要改进的部分，若存在问题，则采用相应的方法进行处理并重新进行求解；若未发现问题，则输出最终结果，并进行后续的数据处理。 2、三维约束平差 三维约束平差是以基线解算所得到的三维静态基线向量为观测值，在平差过程中引入会使GPS 网的尺度、方向和位置发生变化的外部起算数据，从而实现GPS 网成果由基线解算时GPS 卫星星历所采用的参照系（WGS84 ）到特定参照系的转换，得到在特定参照系下的经过用户约束条件约束的点三维空间坐标。 二、南方GPS数据处理软件的平差方式 三维约束平差是指在基线解算后，WGS84坐标系下的三维平差，在三维平差中是不需要当地平面直角坐标系下的已知点坐标，当需要用到WGS84经纬度或空间直角坐标的用户可加载已知点的WGS84空间坐标（如果只有经纬度时，可采用COORD4.1软件进行转换，本站免费提供）进行三维约束平差，即可得到与已知点相匹配的WGS84坐标。 一般情况下，在“已知点坐标录入”窗口中，我们都没有输入WGS8坐标，而只输入当地坐标系下的已知坐标，此时GPS处理软件会自动识取一个坐标点的WGS84坐标进行约束平差。如下图： 如果在某些控制测量中,需要得到精确的WGS84经纬度或空间坐标时,让系统自动识取显然是不行的,此时我们只要为参与平差的已知点的WGS84空间坐标输入后再进行三维平差即可 在这里,我们加入了两个已知点的WGS84空间坐标,三维平差后,列表中会显示两个"固定"字样的点,说明,在进行三维平差中,我们把这两个点做为起算点,进行平差别的未知点。 一、前言本程序是《CASIO fx-5800P计算与道路坐标放样计算》中道路坐标放样计算程序的升级改进版本。原道路坐标放样计算程序只基于道路的单个基本型曲线，有效计算范围仅包括平曲线部分和前后的两条直线段，使用时需要输入平曲线设计参数，无坐标反算桩号功能。改进后的程序名称为：道路中边桩坐标放样正反算程序（全线贯通），增加了可实现全线贯通的数据库功能和坐标反算桩号功能，主要是： 1．使用道路平面数据库子程序，可将一段或若干段道路的交点法格式平面参数（可容易从直线、曲线及转角表中获得）以数据库子程序形式输入计算器，程序在计算时省却了输入原始数据的麻烦； 2．坐标正算方面，输入桩号即可进行道路的中、边桩坐标计算，若输入了测站坐标，还可同时计算全站仪极坐标放样数据（拨角和平距）； 3．坐标反算方面，输入平面坐标，即可计算对应的桩号和距中距离（含左右信息）； 4．对于存在断链的道路，可分段分别编写数据库子程序，然后在主程序中添加一个路段选择的功能即可实现（可参照立交匝道程序中匝道的选择）。程序的特点： 1．可进行中桩坐标的正、反算，程序代码简洁，便于阅读和改写； 2．主程序通过调用数据库子程序，省却了使用时输入平面参数的繁琐； 3．使用数据库子程序，换项目只需改写数据库子程序，程序通用性强。二、道路示例项目基本资料基本资料同《CASIO fx-5800P计算与道路坐标放样计算》第6章HY高速公路第2合同段（合同段起止桩号： K4+800~K9+600）。这里摘取直线、曲线及转角表资料如下. . . 三、程序代码 . . . . . 坐标正反算定义及公式 Corporation standardization office #QS8QHH-HHGX8Q8-GNHHJ8 第六章→第三节→导线测量内业计算 导线计算的目的是要计算出导线点的坐标，计算导线测量的精度是否满足要求。首先要查实起算点的坐标、起始边的方位角，校核外业观测资料，确保外业资料的计算正确、合格无误。 一、坐标正算与坐标反算 1、坐标正算 已知点的坐标、边的方位角、两点间的水平距离，计算待定点的坐标，称为坐标正算。如图6-6 所示，点的坐标可由下式计算： 式中、为两导线点坐标之差，称为坐标增量，即： 【例题6-1】已知点A坐标，=1000、=1000、方位角=35°17＇36.5"，两点水平距离=200.416，计算点的坐标？ 35o17＇36.5"=1163.580 35o17＇36.5"=1115.793 2、坐标反算 已知两点的坐标，计算两点的水平距离与坐标方位角，称为坐标反算。可知，由下式计算水平距离与坐标方位角。 （6-3） （6-4） 式中反正切函数的值域是-90°～+90°，而坐标方位角为0°～360°，因此坐标方位角的值，可根据、的正负号所在象限，将反正切角值换算为坐标方位角。 【例题6-2】=3712232.528、=523620.436、=3712227.860、=523611.598，计算坐标方位角计算坐标方位角 、水平距离。 =62°09＇29.4"+180°=242°09＇29.4" 注意：一直线有两个方向，存在两个方位角，式中：、的计算是过A点坐标纵轴至直线的坐标方位角，若所求坐标方位角为，则应是A点坐标减点坐标。 坐标正算与反算，可以利用普通科学电子计算器的极坐标和直角坐标相互转换功能计算，普通科学电子计算器的类型比较多，操作方法不相同，下面介绍一种方法。 【例题6-3】坐标反算，已知=2365.16、=1181.77、 =1771.03、=1719.24，试计算坐标方位角、水平距离。 键入1771.03-2365.16按等号键［=］等于纵坐标增量，按储存键［］， 键入1719.24-1181.77按等号键［=］等于横坐标增量，按［］键输入，按［］显示横坐标增量，按［］键输入，按第二功能键［2ndF］，再按［］键，屏显为距离，再按［］键，屏显为方位角。 【例题6-4】坐标正算，已知坐标方位角=294°42＇51"， =200.40，试计算纵坐标增量横坐标增量。 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) 5800-9860计算器坐标正反算通用程序 1. 坐标正算主程序(命名为ZBZS) 第1行：Lbl 0:”K=”?K:”BIAN=”? Z:”α=”?B 第2行：Prog “A” 第3行：”X=”:N+Zcos(F+B)◢ 第4行：”Y=”:E+Zsin(F+B)◢ 第5行：”F=”:F◢ 第6行：Goto 0 K——计算点的里程 BIAN——计算点到中桩的距离（左负右正） α——取前右夹角为正 2. 坐标反算桩号和偏距主程序(命名为ZBFS) 第1行：”X1=”? C:”Y1=”?D: ”K1=”?K 第2行：Lbl 0:Prog “A” 第3行：Pol(C-N,D-E) 第4行：List Ans[1]→I 第5行：List Ans[2]→J 第6行：Icos(F-J)→S:K+S→K 第7行：Abs(S)>0.0001=>Goto 0 第8行：”K1=”:K◢ 第9行：”BIAN=”:Isin(J-F)→Z◢ X1——取样点的X坐标 Y1——取样点的Y坐标 K1——输入时为计算起始点(在线路内即可)，输出时为反算点的桩号 Z——偏距（左负右正） 3. 计算坐标子程序(命名为XYF) 为了简洁，本程序由数据库直接调用，上述中的正反算主程序不直接调用此程序第1行：K-A→S：(Q-P)÷L→I 第2行：N+∫(cos(F+X(2P+XI)×90÷π),0,S)→N 第3行：E+∫(sin(F+X(2P+XI)×90÷π),0,S)→E 第4行：F+S(2P+SI)×90÷π→F 4. 数据库（命名为A） 第1行：K≤175.191=>Stop 第2行： 175.191→A:428513.730→N:557954.037→E:92°26′40″→F:0→P:1/240→Q:70.417→L: K≤A+L =>GoTo 1（第一缓和曲线，圆半径为240） 第3行：245.607→A: 428507.298→N:558024.092→E: 100°50′59.4″→F: 1/240→P:1/240→Q: 72.915→L: K≤A+L =>Goto 1（第圆曲线，半径为240） 第4行：318.522→A: 428482.988→N:558092.538→E: 118°15′25.2″→F: 1/240→P: 0→Q: 55.104 →L: K≤A+L =>Goto 1（第二缓和曲线，圆半径为240） 第5行：373.627→A:428453.283→N:558138.912→E:124°50′4.5″→F:0→P:- 1/180→Q:67.222→L: K≤A+L=>Goto 1:Stop（第一缓和曲线，圆半径为180） 第6行：Lbl 1:Prog “XYF” A——曲线段起点的里程 N——曲线段起点的x坐标 E——曲线段起点的y坐标 F——曲线段起点的坐标方位角 P——曲线段起点的曲率(左负右正) Q——曲线段终点的曲率(左负右正) L——曲线段长度（尽量使用长度，为计算断链方便） 说明： （1）在9860中，程序中所有公式和部分函数结果均存储在List Ans列表数组中，要想多次调用最好随公式取出结果，并赋给变量。 （2）正算主程序可以计算一般边桩的坐标，如要计算类似涵洞端墙的坐标需增加第二偏距和转角两个变量。 （3）程序规定，左偏曲线曲率（半径倒数）输入负值，右偏曲线曲率输入正值，直线上点曲率输入0，例如直线段，线元起点和终点均输入0，第一缓和曲线分别输入0和圆半径的倒数，圆曲线均输入半径倒数，第二缓和曲线分别输入圆半径倒数和0，卵形曲线分别输 入对应圆半径的倒数 （4）若是从大里程向小里程的反方向计算，则曲率取正方向时的负值，方位角减去(或加上)180度。 （5）有多个匝道的项目，可随时更改正反算主程序中的红色字体部分来调用其它线路的数据 （6）反算桩号偏差为1mm （7）可以计算任意线型的任意点坐标 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→ 坐标正反算公式 欧阳光明（2021.03.07） 一、GPS数据处理相关术语 1、三维无约束平差 三维无约束平差是以基线解算所得到的三维静态基线向量为观测值，待定参数主要为 GPS 网中点的坐标；同时，利用基线解算时随基线向量一同输出的基线向量的方差阵，形成平差的随机模型，最终形成平差完整的数学模型。随后对所形成的数学模型进行求解，根据平差结果来确定观测值中是否存在粗差，数学模型是否有需要改进的部分，若存在问题，则采用相应的方法进行处理并重新进行求解；若未发现问题，则输出最终结果，并进行后续的数据处理。 2、三维约束平差 三维约束平差是以基线解算所得到的三维静态基线向量为观测值，在平差过程中引入会使 GPS 网的尺度、方向和位置发生变化的外部起算数据，从而实现 GPS 网成果由基线解算时 GPS 卫星星历所采用的参照系（ WGS84 ）到特定参照系的转换，得到在特定参照系下的经过用户约束条件约束的点三维空间坐标。二、南方GPS数据处理软件的平差方式 三维约束平差是指在基线解算后，WGS84坐标系下的三维平差，在三维平差中是不需要当地平面直角坐标系下的已知点坐标，当需要用到WGS84经纬度或空间直角坐标的用户可加载已知点的WGS84空间坐标（如果只有经纬度时，可采用 COORD4.1软件进行转换，本站免费提供）进行三维约束平差，即可得到与已知点相匹配的WGS84坐标。 一般情况下，在“已知点坐标录入”窗口中，我们都没有输入WGS8坐标，而只输入当地坐标系下的已知坐标，此时GPS处理软件会自动识取一个坐标点的WGS84坐标进行约束平差。如下图： 如果在某些控制测量中,需要得到精确的WGS84经纬度或空间坐标时,让系统自动识取显然是不行的,此时我们只要为参与平差的已知点的WGS84空间坐标输入后再进行三维平差即可 在这里,我们加入了两个已知点的WGS84空间坐标,三维平差后,列表中会显示两个"固定"字样的点,说明,在进行三维平差中,我们把这两个点做为起算点,进行平差别的未知点。 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 ↙ 一、坐标正算基本公式 (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、SJKl （平面数据库） (14) IK SJK2 （纵面数据库） (14) 12、SJK3 （左路宽度数据库） (15) 13、SJK4 （右路宽度数据库） (15) 14、SJK5 （横坡数据库） (16) 15、SJK6 （下边坡数据库） (16) 16、SJK7 （左上边坡数据库） (17) 17、SJK8 （右上边坡数据库） (18) 五、后记 (19) CASIO 5800计算器公路工程测量程序 4 一、正算所涉及的计算公式 图表1 在图1中，A点为回旋曲线起点，B点为回旋曲线止点，I点为所求坐标 点。 设： A点的X坐标为XA , Y坐标为Y A, A点的切线方位角为α, A点的曲率为P A点的里程为L A, B点的曲率为p 8, B点的里程为L B, I点的曲率为p∣, I点的里A, 程为L I O I点的切线角为B o 由于回旋线上各点曲率半径Ri和该点至曲线起点的距离L成反比。故此任意点 的曲率为； Q =丄=土（C为常数）（1 ） R i C 由式（1）可知，回旋曲线任意点的曲率按线性变化，由此回旋曲线上里程为 L 点的曲率为； i 一、坐标正算与坐标反算 1、坐标正算 已知点的坐标、边的方位角、两点间的水平距离，计算待定点的坐标，称为坐标正算。如图6-6 所示，点的坐标可由下式计算： 式中、为两导线点坐标之差，称为坐标增量，即： 【例题6-1】已知点A坐标，=1000、=1000、方位角 =35°17＇36.5"，两点水平距离=200.416，计算点的坐标？ 35o17＇36.5"=1163.580 35o17＇36.5"=1115.793 2、坐标反算 已知两点的坐标，计算两点的水平距离与坐标方位角，称为坐标反算。如图6-6可知，由下式计算水平距离与坐标方位角。 （6-3） （6-4）式中反正切函数的值域是-90°～+90°，而坐标方位角为0°～360°，因此坐标方位角的值，可根据、的正负号所在象限，将反正切角值换算为坐标方位角。 【例题6-2】=3712232.528、=523620.436、 =3712227.860、=523611.598，计算坐标方位角计算坐标方位角 、水平距离。 =62°09＇29.4"+180°=242°09＇29.4" 注意：一直线有两个方向，存在两个方位角，式中：、 的计算是过A点坐标纵轴至直线的坐标方位角，若所求坐标方位角为，则应是A点坐标减点坐标。 坐标正算与反算，可以利用普通科学电子计算器的极坐标和直角坐标相互转换功能计算，普通科学电子计算器的类型比较多，操作方法不相同，下面介绍一种方法。 【例题6-3】坐标反算，已知=2365.16、=1181.77、 =1771.03、=1719.24，试计算坐标方位角、水平距离。 键入1771.03-2365.16按等号键［=］等于纵坐标增量，按储存键［］， 键入1719.24-1181.77按等号键［=］等于横坐标增量，按［］键输入，按［］显示横坐标增量，按［］键输入，按第二功能键［2ndF］， CASIO fx-4800P、fx-5800P型计算器用于线路施工 曲线中线点坐标的计算程序 中铁十局三建公司工程技术部 摘要：本文介绍了CASIO fx-4800P 、fx-5800P型计算器程序编制用于铁路、公路曲线线路内任意中线点的坐标计算程序及使用方法。本计算程序具有操作简便、计算快捷、应用广泛等特点、极大地减轻了测量工作者的内业工作量，对于测量工作者有较大的参考和指导作用。关键词：曲线线路施工测量计算程序 1．概述 过去，线路中线施工放样基本依靠经纬仪和钢尺了来进行角度及距离测量。对于曲线线路一般的测量方法是：经纬仪置于某一中线点上，采用偏角法拨角再用钢尺量距来定出中线点。随着电子技术进步和经济发展，测量仪器和测量方法的不断改进，目前，全站仪已广泛地应用于工程施工测量中，极大的提高了测量工作效率。但是，在进行铁路、公路工程的曲线线路施工测设时，需要在线路所在区域建立统一坐标系或独立坐标系，利用坐标变换的方法，将整个曲线的三个部分（第一缓和曲线、中间圆曲线、第二缓和曲线）统一到同一坐标系中。根据坐标系的建立，计算出整个曲线内任意点的坐标，再采用全站仪利用极坐标方法进行施工放样。前提是首先利用计算器计算出各中线点坐标，然后才能进行放样。而普通型计算器不仅计算速度慢，且要求计算者必须正确地记忆很多计算公式，计算繁琐而且容易出错，满足不了现场测设工作的要求。为了能够快速准确地为全站仪提供测设 数据，发挥全站仪快速测设的特点，提高测量工作效率，应采用可编程的计算器，编制计算程序。本文主要介绍应用CASIO fx-4800P型计算器的计算程序，供公司测量同行们参照使用。 2．计算程序 QXZBJS（文件名：曲线坐标计算fx-4800P） Defm2:R：L：A：N“ZH：X=”：E“ZH：Y=”：F：“FWJ=”：K“ZH：LC=”： P=L2/(24R)-L4/(2688R3):M=L/2-L3/(240R2):T“T”=（R+P）tng(A/2)+M ◢G=RAπ/180：“S”S=G+L◢ LbiA：｛C，V｝：C“CSDLC=”：V“HXPJ=”：D=C-K：D≤L＝＞I=D-D5/（40R2L2）：U=D3/(6RL)-D7/（336R3L3):J=√（I2+U2）：Goto1：≠＞D≤G＝＞O=90（2D-L）/( Rπ)：I=RsinO+M：U=R（1-cosO）+P：J=√（I2+U2）：Goto2：≠＞D=S-（C-K）： = D-D5/（40R2L2）：Z[2]=D3/6RL-D7/（336R3L3）：I=T+（T-Z[1]）cosA-Z[2]sinA：U=（T-Z[1]）sinA+Z[2]cosA：J=√（I2+U2）：Goto3：Lbi1:{Q}:Q“Z=1；Y=2”：Q=1＝＞Q=F-30D2/ （RLπ）：H=F-90D2/ （RLπ）：≠＞Q=F+30D2/ （RLπ）：H=F+90D2/（ RLπ）⊿ Goto4：Lbi2:{Q}:Q“Z=1；Y=2”：Q=1＝＞Q=F-tng-1(U/I)：H=F-O：≠＞Q= F+tng-1(U/I)：H=F+O⊿Goto4：Lbi3:{Q}:Q“Z=1；Y=2”：Q=1＝＞Q=F-tng-1(U/I):H=F-(A-90(S-(C-K))2/ （RLπ）): ≠＞Q= F+tng-1(U/I):H=F+(A-90(S-(C-K))2/ （RLπ）):⊿ Goto4： Lbi4:B=90+H:H＜0＝＞H“QXFWJ”=B+360◢≠＞H≥360＝＞H“QXFWJ”=H-360 ◢≠＞H“QXFWJ”=H◢⊿Goto5： Lbi5: X“CSD:X”=JcosQ+N+VcosB◢ Y“CSD:Y”=JsinQ+E+VsinB◢坐标正反算程序计算器
卡西欧5800计算器坐标正反算程序
FX5800道路路线测量程序
坐标正反算计算公式
公路测量卡西欧5800万能程序
坐标正反算定义及公式
5800简单全线坐标计算程序
5800-9860计算器坐标正反算通用程序
Fx5800计算器公路测量程序设计
卡西欧5800测量坐标计算程序
2021年坐标正反算计算公式
5800计算器全线坐标计算放样程序(修改第三版)
卡西欧5800公路坐标正反算程序
坐标正反算
CASIO fx-5800p测量程序