fx-5800P 坐标正反算程序

坐标正反算程序

Lb1 0:Fix 3:“X0”?X:“Y0”?Y:“ZS2,FS1”?N?If N=1:Then Goto 1:Else Goto 2:IfEnd?

Lb1 1:“X1”?A:“Y1”?B:Pol(A-X,B-Y)?

If J<0:Then J+360→J:IfEnd?

“JIAODU=”J?DMS◢

“R=”:I◢

Goto 0?

Lb1 2:“JIAODU=”?C:“R=”?D:X+Dxcos(C)→E:Y+Dxsin(C)→F?

“X=”:E◢

“Y=”:F◢

Goto 0?

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◣

工程测量坐标正反算通用程序(终极篇)

第五篇坐标正反算通用程序(终极篇) 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?DMS◢ 第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):Icos(F-J)→S:K+S→K 第4行：Abs(S)>0.0001=>Goto 0 第5行：”K1=”:K◢ 第6行：”BIAN=”:Isin(J-F)→Z◢ X1——取样点的X坐标 Y1——取样点的Y坐标 K1——输入时为计算起始点(在线路内即可)，输出时为反算点的桩号 Z——偏距（左负右正） 注：在9860或9960中需将第3行替换为Pol(C-N,D-E): List Ans[1]→I ：List Ans[2]→J：Icos(J-F)→S:K+S→K，正反算主程序所有输入赋值多加一赋值符号（→），其他所有除数据库外的程序均保持不变 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+S I)×90÷π→F 第5行：F<0=>F+360→F: F>360=>F-360→F 4. 数据库（命名为A） 第1行：K≤175.191=>Stop（超出后显示Done） 第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（第一缓和曲线） 第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（圆曲线）

线元法线路坐标正反算程序

经苦心钻研，奋战多日，终于编写出了代码短，速度快，精度高，功能全的线路坐标正反算程序，欢迎试用并提出宝贵意见。 功能简介及特点： 1、选用高斯-勒让德公式作计算内核，保证精度，模块化设计，便于扩充功能。 2、线元数据可自动从数据库调用，也可手工输入。 3、可管理多条线路，如里程不在线路或线元范围，将警告里程偏大、偏小。 4、边桩计算设计为导线式递推方式，可用于由一个中桩推出结构物所有角点坐标。 5、反算实现了智能化操作，只需输入线路号（或手工输线元资料）、坐标，不需近似里程，即可自动从起点向后开始试算出里程、位置，如对算出里程、位置表示怀疑，还可以让计算器从终点起再向前试算下一个可能的位置（匝道、回头曲线同一坐标可能会有一个以上结果）。第三次及以后试算才要求输入近似里程。 6、程序代码规范简洁，便于阅读、理解。 完整程序清单: ZFS %正反算主程序 B=.1739274226:C=.5-B: Lbl 1:U"0 ZS 1 FS"=0=>Prog "ZS": ≠>U=1=>Prog"FS":≠>Goto 1

ZS %正算子程序 {K}:Prog"ZZ":I=0:{I}:I"L"≠0=>"Prog"WY":≠>Prog"ZB" FS %反算子程序 {KVW}:V"XC"W"YC":Lbl 2:Prog "ZZ":I=V-S:J=W-T:Pol(I,J: J=J-F:K=K+Rec(I,J:AbsI<1m=>Prog"WZ":≠>Goto 2Δ M=0:{M}:M"0 NEXT"=0=>U=U+1:Goto 2:≠>U=1 ZZ %高斯法中桩子程序(4节点) Prog"XL":M=K-L:O=(P-R)÷2PQR: D=.0694318442:E=.3300094782:F=1:G=1-E:H=1-D: I=5:Lbl 1:C[I]=A+MrC[I](1÷P+OMC[I]:Dsz I:Goto 1: S=X+M(BcosD+CcosE+CcosG+BcosH: T=Y+M(BsinD+CsinE+CsinG+BsinH WY %外移点计算子程序 Lbl 1:J=90:{J}:J=F+J"<":F=J:S=S+Rec(I,J:T=T+J: Prog"ZB":I=0:{I}:I"L"≠0=>Goto 1 WZ %位置显示子程序 "KJ":K:Pause 1:J◢ ZB %坐标显示子程序 "XY":S:Pause 1:T◢ YC %异常处理子程序 U=1=>K=L:U=2Δ U=3=>K=M:U=4Δ

坐标正反算计算公式

坐标正反算公式

一、GPS数据处理相关术语 1、三维无约束平差 三维无约束平差是以基线解算所得到的三维静态基线向量为观测值，待定参数主要为GPS 网中点的坐标；同时，利用基线解算时随基线向量一同输出的基线向量的方差阵，形成平差的随机模型，最终形成平差完整的数学模型。随后对所形成的数学模型进行求解，根据平差结果来确定观测值中是否存在粗差，数学模型是否有需要改进的部分，若存在问题，则采用相应的方法进行处理并重新进行求解；若未发现问题，则输出最终结果，并进行后续的数据处理。 2、三维约束平差 三维约束平差是以基线解算所得到的三维静态基线向量为观测值，在平差过程中引入会使GPS 网的尺度、方向和位置发生变化的外部起算数据，从而实现GPS 网成果由基线解算时GPS 卫星星历所采用的参照系（WGS84 ）到特定参照系的转换，得到在特定参照系下的经过用户约束条件约束的点三维空间坐标。 二、南方GPS数据处理软件的平差方式

三维约束平差是指在基线解算后，WGS84坐标系下的三维平差，在三维平差中是不需要当地平面直角坐标系下的已知点坐标，当需要用到WGS84经纬度或空间直角坐标的用户可加载已知点的WGS84空间坐标（如果只有经纬度时，可采用COORD4.1软件进行转换，本站免费提供）进行三维约束平差，即可得到与已知点相匹配的WGS84坐标。 一般情况下，在“已知点坐标录入”窗口中，我们都没有输入WGS8坐标，而只输入当地坐标系下的已知坐标，此时GPS处理软件会自动识取一个坐标点的WGS84坐标进行约束平差。如下图：

如果在某些控制测量中,需要得到精确的WGS84经纬度或空间坐标时,让系统自动识取显然是不行的,此时我们只要为参与平差的已知点的WGS84空间坐标输入后再进行三维平差即可 在这里,我们加入了两个已知点的WGS84空间坐标,三维平差后,列表中会显示两个"固定"字样的点,说明,在进行三维平差中,我们把这两个点做为起算点,进行平差别的未知点。

5800曲线道路全线坐标正反算(已验证)

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])

坐标正反算定义及公式

坐标正反算定义及公式 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-9860计算器坐标正反算通用程序

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）可以计算任意线型的任意点坐标

5800公路坐标正反算程序

坐标正反算程序 程序变量说明: 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'

5800交点法正反算程序

5800 交点法正反算程序 自编5800 支距及偏角法正反算程序 该程序使用支距偏角法，正算速度一般，反算速度较辛普森公式反算速度稍快 计算范围从直缓到缓直，如若缓直后为直线则无需输入变量直接继续计算 此程序在正算循环里程超限上占用字节较多，适用公路路基临时输入变量放样 运行只需从主程序进入，请各位高手多多指正 主程序“ZS FS” 20→Dim Z “1→NEW,2→OLD”？Z ↙ IF Z=1:Then Goto 1: Else Goto 2 if end↙选择输入选一手动输入曲线要数选二提取数据库 输入交点X,Y LbI 1:”JDX”?T:“JDY”?B 输入切线1，2的长 “T1”?C:”T2”?D 输入前后缓和曲线的长 “L1”?E: “L2”?F “ZHL”?G 输入直缓里程 “FWJ”Q 输入直缓点切线方位角 “-PJ+”?S 输入偏角左负右正 ?R: 输入圆曲线半径 Goto 4↙ ↙ LbI2:prog “SJK”:goto 4 ↙输入观测点里程输入N值1为正算其他为反算 LbI4: “KCD”?K: ?N ： If N=1: T hen Goto 5: Else prog “ZSXH” Goto 6: If end↙

LbI5: “X”?M: “Y”?Z Prog “FSXH”: Goto 7: If end↙选择反算后输入测点的X,Y 假设里程与正算输入一起，无需另输入 LbI6：“W”?W求边桩输入边距左-右+ X+WCos(P+90) →X: “X=”:X ◢显示所求X值 ◢显示所求Y值 Y+WSin(P+90) →Y: “Y=”:Y ◢显示所求点切线方位角 P<0 => P+360→P: “FW=”:P Goto 4 LbI7: “KCD=”:K ◢显示反算出的里程 “W=”:W◢显示反算出的边距 Goto 4↙ 子程序正算循环“ZSXH” ↙ R(Abs(S)-90E÷(∏R)-90F÷(∏R)) ∏÷180→Ｊ If K<=(说明：小于等于号)E+G+J: Then Goto 1: Else Goto 3: If end↙ LbI1: Q→H:T-C*Cos(H)→X:B-C*Sin(H)→Y: 根号(R∏) →Ａ↙ If K<=G+E:Then K-G→L:prog “B” :P→P:prog “C”:Goto 6: Else E→L: prog “B”:P→P: prog “C”:X→X:Y→Y:Goto 2: If end↙ ↙ LbI2: K-G-E→L: prog “D”:P→P: prog“C”: Goto 6 根号(RF)→Ａ↙ LbI3: Q+S→H:T+D*Cos(H)→X:B+D*Sin(H)→Y: If K<=G+E+J+F: Then G+E+J+F-K→L: prog “B”：

2021年坐标正反算计算公式

坐标正反算公式 欧阳光明（2021.03.07） 一、GPS数据处理相关术语 1、三维无约束平差 三维无约束平差是以基线解算所得到的三维静态基线向量为观测值，待定参数主要为 GPS 网中点的坐标；同时，利用基线解算时随基线向量一同输出的基线向量的方差阵，形成平差的随机模型，最终形成平差完整的数学模型。随后对所形成的数学模型进行求解，根据平差结果来确定观测值中是否存在粗差，数学模型是否有需要改进的部分，若存在问题，则采用相应的方法进行处理并重新进行求解；若未发现问题，则输出最终结果，并进行后续的数据处理。 2、三维约束平差 三维约束平差是以基线解算所得到的三维静态基线向量为观测值，在平差过程中引入会使 GPS 网的尺度、方向和位置发生变化的外部起算数据，从而实现 GPS 网成果由基线解算时 GPS 卫星星历所采用的参照系（ WGS84 ）到特定参照系的转换，得到在特定参照系下的经过用户约束条件约束的点三维空间坐标。二、南方GPS数据处理软件的平差方式

三维约束平差是指在基线解算后，WGS84坐标系下的三维平差，在三维平差中是不需要当地平面直角坐标系下的已知点坐标，当需要用到WGS84经纬度或空间直角坐标的用户可加载已知点的WGS84空间坐标（如果只有经纬度时，可采用 COORD4.1软件进行转换，本站免费提供）进行三维约束平差，即可得到与已知点相匹配的WGS84坐标。 一般情况下，在“已知点坐标录入”窗口中，我们都没有输入WGS8坐标，而只输入当地坐标系下的已知坐标，此时GPS处理软件会自动识取一个坐标点的WGS84坐标进行约束平差。如下图： 如果在某些控制测量中,需要得到精确的WGS84经纬度或空间坐标时,让系统自动识取显然是不行的,此时我们只要为参与平差的已知点的WGS84空间坐标输入后再进行三维平差即可 在这里,我们加入了两个已知点的WGS84空间坐标,三维平差后,列表中会显示两个"固定"字样的点,说明,在进行三维平差中,我们把这两个点做为起算点,进行平差别的未知点。

与控制点程序结合坐标正反算通用程序

GAUSSLE坐标正反算fx－4850程序 源程序 1.正算主程序GSZS I"X0"：S"Y0"：O"K0"：G"F0"：H"KN"：P"R0"：R"RN"：Q”Q(-Z +Y)” : D=(P-R)÷(2(H-O)PR)：KL”L(-Z +Y)” :M”ANG(YJJ)”: J=Abs(K-O)： Prog"SUB1"：”FWJ=”:F=F-M:”X=”:X=X:Pause0: ”Y=”:Y=Y◢ 2. 反算主程序GSFS XY：Z[2]=X：Z[3] =Y： I"X0"：S"Y0"：O"K0"：G"F0"：H"KN"：P"R0"：R"RN"：Q”Q(-Z +Y)” : D=(P-R)÷(2(H-O)PR): J=Abs((Y-S)cos(G-90)－(X-I)sin(G-90))： L=0：M”M(YJJ)”=90: Lbl 0：Prog "SUB1"： L=(Z[3]-Y)cos(G-90+QJ(1÷P+JD)×180÷π)-(Z[2]-X)sin(G-90+QJ(1÷P +JD) ×180÷π)：AbsL<1E-6=>Goto1：≠>J=J+L：Goto 0Δ←┘ Lbl 1：L=0：Prog "SUB1"：L=(Z[3]-Y)÷sinF:”K=”:K=O+J:Pause0:”L=”:L=L◢ 3. 正算子程序(SUB1) Defm 4: A=0.1184634425:B=0.2393143352:Z[4]=0.2844444444:C=0.0469100770:E=0.230 7653449:Z[1]=0.5: X=I+J(Acos(G+QCJ(1÷P+C JD)×180÷π)+Bcos(G+Q EJ(1÷P+E JD)×180÷π)+Z[4]cos(G +QZ[1]J(1÷P+Z[1]JD)×180÷π)+Bcos(G+Q(1-E)J(1÷P+(1-E)JD)×180÷π)+Acos(G+Q (1-C)J(1÷P+(1-C)JD) ×180÷π)): Y=S+J(Asin(G+QCJ(1÷P+C JD)×180÷π)+Bsin(G+Q EJ(1÷P+E JD)×180÷π)+Z[4]sin(G +QZ[1]J(1÷P+Z[1]JD)×180÷π)+Bsin(G+Q(1-E)J(1÷P+(1-E)JD)×180÷π)+Asin(G+Q (1-C)J(1÷P+(1-C)JD) ×180÷π))： F=G+QJ(1÷P+JD) ×180÷π+M：X=X+LcosF：Y=Y+LsinF 4. 曲线元要素数据库：DAT-M K≥O=>K＜H=> I=**:S=**:O=**:G=**:H=**:P=**:R=**:Q=**⊿⊿←┘ K≥O=>K＜H=> I=**:S=**:O=**:G=**:H=**:P=**:R=**:Q=**⊿⊿←┘ …………………………… K≥O=>K≤H=> I=**:S=**:O=**:G=**:H=**:P=**:R=**:Q=**⊿⊿←┘ (注：如有多个曲线元要素继续添加入数据库DAT-M中) 5、M线（坐标正算）组合程序MG-ZB Prog”DAT-M”:Prog”GSZS” 6、M线（坐标计算－放样）组合程序MG-FY Prog”MG-ZB”: Prog”FY” 7、M线（坐标反算）组合程序M-GSFB Prog”DAT-M”:Prog”GSFS”

坐标正反算

一、坐标正算与坐标反算 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］，

fx-5800p全线坐标正反算带高程计算程序(线元法)

曲线任意里程中边桩坐标正反算(CASIO fx-5800P计算器)程序 一、程序功能及原理 1.功能说明：本程序由一个主程序(TYQXJS)和五个子程——正算子程序(SUB-ZS)、反算子程序(SUB-FS)等构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。本修改版程序既可实现正算全线贯通，亦可实现反算全线贯通。本程序在CASIO fx-5800P计算器运行。 2．计算原理：利用Gauss-Legendre 5点通用公式正算线路中边桩坐标、线外测点至曲线元起点和终点的垂距的符号是否相异（即Dca×Dcb<0=>该测点在其线元内）进行判断并利用该线元要素反算中桩里程、支距，最后计算出放样数据。 二、源程序 1.主程序(TYQXJS)(A) Deg:fix 3 119→DimZ “INPUT(0) Or DATA(Else)”?I Lbl 0:“1.SZ=>XY,2.XY=>SZ，3.TF=>CK,4.SD=>FY，5.TW=>FY”?N If N=1 Or N=5:Then Goto 1 Else If N=2 Or N=3 Or N=4:Then Goto 2 Else Goto 3 IfEnd:IfEnd Lbl 1:“K(m)=”?S If S<0:Then Goto 0:IfEnd “JL(m)=”?Z If Z≠0:Then “ANGLE→R(Deg)=”?M:IfEnd If I=0:Then Prog “DAT1”:Else Prog “DAT2”:IfEnd S-O→W:If W<0:Then Goto 0:Else If W>H:Then Goto 0:IfEnd:IfEnd Prog “SUB-ZS”:Prog “SUB-GC” If Z<0:Then“XL(m)=”:X◢“YL(m)=”:Y◢ If N=5:Then Prog “SUB-TW”:IfEnd Else If Z>0:Then “XR(m)=”:X◢“YR(m)=”:Y◢ If N=5:Then Prog “SUB-TW”:IfEnd Else “X(m)=”:X◢“Y(m)=”:Y◢“Hs(m)=”:L◢“FWJ=”: F?DMS◢ IfEnd:IfEnd

卡西欧计算器坐标的正反算

可以算任意斜交涵洞轴线的坐标，增加T为斜交角度，规定T为涵轴右侧方向与“线路前进方向切线”之间的夹角，当涵轴与线路正交时，T=90,其他操作与原程序一样； 1. 正算子程序(SUB1) [color=Red]A=0.26：B=0.74：K=0.02：L=0.82：F=1-L： M=1-K：X=U+W(Acos(G+57.2958QKW(1/P+KWD))+Bcos(G+57.2958QLW(1/P+LW D))+Bcos(G+57.2958QFW (1/P+FWD))+Acos(G+57.2958QMW(1/P+MWD)))：Y=V+W(Asin(G+57.2958QKW(1/ P+KWD))+Bsin(G+ 57.2958QLW(1/P+LWD))+Bsin(G+57.2958QFW(1/P+FWD))+Asin(G+57.2958QMW (1/P+MWD)))：F=G+57.2958QW(1/P+ WD)+90：X=X+Zcos(F-90+T)：Y=Y+Zsin(F-90+T) 2. 反算子程序(SUB2) W=Abs((Y-V)cos(G-90)-(X-U)sin(G-90))：Z=0：Lbl 0：Prog "SUB1"：L=(G-90)+5 7.2958QW(1/P+ WD)：Z=(J-Y)cosL-(I-X)sinL：AbsZ<1E-6=>Goto1：≠>W=W+Z：Goto 0Δ←┘ Lbl 1：Z=0：Prog "SUB1"：Z=(J-Y)÷sinF 二.增设数据库程序（ＳＪＫ主程序） Lb1 4："1.SZ => XY"："2.XY => SZ"：｛NS｝：S∠下一线元起点里程＝＞O =本线元起点里程：U=本线元起点X：V=本线元起点Y：G=本线元起算方位角：H =本线元长度：P=起点曲率半径：R=终点曲率半径：Q=0或1、-1：Prog“TYQXJS”：Goto0Δ←┘(第一线元数据要素) S∠下一线元起点里程＝＞O=本线元起点里程：U=本线元起点X：V=本线元起点Y：G=本线元起算方位角：H=本线元长度：P=起点曲率半径：R=终点曲率半径：Q=0或1、-1：Goto0Δ←┘(第二线元数据要素)

线路坐标正算及反算程序

一.程序清单: 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