立交匝道坐标放样正反算casiofxp程序带数据库功能secret

可以算任意斜交涵洞轴线的坐标，增加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Δ←┘(第二线元数据要素)S∠下一线元起点里程＝＞O=本线元起点里程：U=本线元起点X：V=本线元起点Y：G=本线元起算方位角：H=本线元长度：P=起点曲率半径：R=终点曲率半径：Q=0或1、-1：Goto0Δ←┘(第三线元数据要素)。

曲线任意里程中边桩坐标正反算及高程、隧道断面超欠挖程序(CASIO fx-5800P计算器)程序(加数据库版)1.正算主程序[KP-XY]"DK="?K:"PJ="?P:Prog"SJK""XS=":X+P cos(Z[4]+90)→X◢"YS=":Y+P sin(Z[4]+90)→Y◢"QX-FWJ=":Z[4]►DMS◢Prog"SQX":O→I:O→J:"XC="?A:"YC="?B:Pol(X-A,Y-B)J<0＝＞360+J→J"FWJ=":J►DMS◢"P=":I◢2.反算主程序[XY-KP]"DK="?K:"X:"?U:"Y:"?VLbl 0Prog"SJK":O→I:O→J:Pol(U-X,V-Y):J<0＝＞J+360→JIcos(J-Z[4])→Z[5]If Abs(Z[5])≥0.001:Then K+Z[5]→K:Goto 0Else Isin(J-Z[4])→P:Ifend"DK=":K◢"PJ=":P◢Prog"SQX"Return3.隧道超欠挖主程序[CQW]Prog"XY-KP"："Q="?QIf Q≥H+6.471:Then "X=":Abs(P）-√(4.26²-(Q-H-4.34)²)→X◢"Y=":Abs(Q-H-4.34)-√(4.26²-P²)→Y◢Eles If Q≥H+2.17:Then "X=":Abs（P）-(√(7.41²-(H+2.77-Q)²)-2.73）→X◢Eles If Q≥H-0.415：Then "X=":Abs（P）-(4.65-(H+2.17-Q)×0.22÷2.58)→X◢"Y=":H-0.42-Q→Y◢Eles If Q≥H-3.2：Then "Y(YG)=":(H+5.2-Q)-√(7.1²-P²)→Y◢IfEnd:IfEnd:IfEnd:IfEnd4.综合洞超欠挖主程序（ZHD)Prog"XY-KP":K-302480→D:"Q="?QIf Q≥H＋4.8:Then "X=":Abs(D)-√(2.5²-(Q-H-4.8)²)→X◢"Y=":Abs(Q-H-4.8)-√(2.5²-D²)→Y◢Else If Q≥H:Then "X=":2.5-Abs(D)→X◢"Y=":H-Q→Y◢说明：302480综合洞里程，2.5半径，4.8总高度5.竖曲线主程序[SQX]0.003(6.984/1000)→I:0.006(-6.4/1000)→J:20000(15000)→R:301420(42350)→L:R Abs(I-J)÷2→Z[2]:RAbs(I-J)÷(I-J)→RIf K≤L-Z[2]:Then 0→Z:Else If K≥L+Z[2]:Then 0→Z:J→I:Else K-L+Z[2]→Z IfEnd:IfEnd"H=":200.621(164.958 )-I(L-K)-Z²÷2÷R→H◢说明：0.003前坡坡度，0.006后坡坡度，20000竖曲线半径，301420竖曲线交点里程，200.621竖曲线交点高程（有竖曲线使用考虑后的高程）0.003→I:0.006→J:0→R:301420→L:RAbs(I-J)÷2→Z[2]:RAbs(I-J)÷(I-J)→RIf K≤L-Z[2]:Then 0→Z:Else If K≥L+Z[2]:Then 0→Z:J→I:Else K-L+T→Z: IfEnd:IfEnd"H=":200.621-I(L-K)→H◢未有竖曲线6.数据库程序[SJK]4500（4500）→R:340（120）→G:157°25′52.55″（220.350722°）→Z:300412.156（43227.29）→F:-23°37′53.48″（6.23477）→O1→Z[1]: 0<0＝＞-1→Z[1]K-F→SIf S≤0:Then S→I:Z+180→J:Z→Z[4]:2953202.491（3678.283）+Icos(J)→X:497703.0998（1696.445）+Isin(J)→Y:Else If S≤G:Then 30Z[1]S²÷(πRG)→W:√((S-S^5÷(40R²G²))²+(S^3÷(6RG))²)→I:Z+W→J:Z+3W→Z[4]:2953202.491（3678.283）+Icos(J)→X:497703.0998（1696.445）+Isin(J)→Y:Else If S≤2196.015（622.3879）-G:Then S-G→S:90Z[1]S÷(πR)→W:2Rsin(Abs(W))→I:Z+90Z[1]G÷π÷R+W→J:J+W→Z[4]:2952890.216（3587.499）+Icos(J)→X:497837.5235（1617.971）+Isin(J)→Y:Else If S≤2196.015:Then 2196.015-S→S:30Z[1]S²÷(πRG)→W:√((S-S^5÷(40R²G²))²+(S^3÷(6RG))²)→I:Z+0+180-W→J:Z+0-3W→Z[4]:2951406.6187（3229.258）+Icos(J)→X:498932.0473（1266.092）+Isin(J)→Y:Else If S＞2196.015:Then S-2196.015→I:Z+O→Z[4]:2951406.6187（3229.258）+Icos(Z[4])→X:498932.0473（1266.092）+Isin(Z[4])→Y:IfEnd:IfEnd:IfEnd:IfEnd:IfEndReturn说明：4500→R曲线半径，340→G缓和曲线长度，157°25′52.55″→Z线路方位角，-23°37′53.48″→O转向角（沿线路方向左负右正），300412.156→F直缓点里程（ZH），2953202.491，497703.0998沿线路方向直缓点坐标（ZH），2196.015曲线长度，2952890.216，497837.5235沿线路方向缓圆点坐标（HY），2951406.6187，498932.0473沿线路方向缓直点坐标（HZ）。

线路测量坐标正、反算计算原理 及卡西欧fx-5800P 程序说明一、计算原理在各测量书中对于坐标正算的相关计算式均有说明，故在此不做详解，仅对正算过程中需要用到的原理及公式做一汇总。

对于坐标反算，虽然都采用无限趋近原理进行计算，但计算方式各有差别，本文仅对其中一种自认为相对简单易懂并便于操作的原理进行详解。

1.1 坐标转换[1]如图1，设X P 、Y P 为P 点在国家控制网坐标系中的坐标；x P 、y P 为P 点在工程独立控制网坐标系中的坐标。

X O 、Y O 为工程独立坐标系原点o 在国家坐标系中的坐标，Δα为两坐标系纵坐标轴的交角。

如果一条边在国家坐标系中的坐标方位角为A ，而在工程独立坐标系中的坐标方位角为α，则：(1-1)当由工程独立坐标系坐标换算至国家坐标系坐标时，换算公式为：(1-2)当由国家坐标系坐标换算至工程独立坐标系坐标时，也可使用式(1-2)，此时应将X 、Y 与x 、y 互换，且。

1.2 坐标方位角关系计算 1.2.1 正、反坐标方位角[2]一条直线的坐标方位角与直线的前进方向有关，沿直线前进方向的坐标方位角称为正坐标方位角，与其相反方向的坐标方位角称为反坐标方位角。

如图2，由于轴子午线之间是互相平行的，因此同一直线的正、反坐标方位角相差180°，即：(2-1)当时，取“+”号；当时，取“-”号。

1.2.2 坐标方位角的推算[3]1.2.2.1 转折角为右角如图3(a)，为已知边坐标方位角，为推算边的坐标方位角，为该两边所夹的右角，则：(2-2)1.2.2.2 转折角为左角如图3(b)，为已知边坐标方位角，为推算边的坐标方位角，(2-3)无论用右角还是左角推算，如遇出现负数的情形，应加上。

当推算出的坐标方位角超过，则应减去。

如果已知导线起始边的坐标方位角和导线各右角或左角时，可直接用以下公式计算第n边的坐标方位角：按右角：(2-4)按左角：(2-5)1.3基本曲线要素及坐标计算各类测量书均对曲线要素及坐标计算有详细推到过程，此处仅列出计算公式，便于程序编制说明。

C A S I O f x-5800P计算与道路坐标放样计算
一、前言
本程序是《CASIO fx-5800P计算与道路坐标放样计算》中道路坐标放样计算程序的升级改进版本。

原道路坐标放样计算程序只基于道路的单个基本型曲线，有效计算范围仅包括平曲线部分和前后的两条直线段，使用时需要输入平曲线设计参数，无坐标反算桩号功能。

改进后的程序名称为：道路中边桩坐标放样正反算程序（全线贯通），增加了可实现全线贯通的数据库功能和坐标反算桩号功能，主要是：
1．使用道路平面数据库子程序，可将一段或若干段道路的交点法格式平面参数（可容易从直线、曲线及转角表中获得）以数据库子程序形式输入计算器，程序在计算时省却了输入原始数据的麻烦；
2．坐标正算方面，输入桩号即可进行道路的中、边桩坐标计算，若输入了测站坐标，还可同时计算全站仪极坐标放样数据（拨角和平距）；
3．坐标反算方面，输入平面坐标，即可计算对应的桩号和距中距离（含左右信息）；
4．对于存在断链的道路，可分段分别编写数据库子程序，然后在主程序中添加一个路段选择的功能即可实现
程序的特点：
1．可进行中桩坐标的正、反算，程序代码简洁，便于阅读和改写；
2．主程序通过调用数据库子程序，省却了使用时输入平面参数的繁琐；
3．使用数据库子程序，换项目只需改写数据库子程序，程序通用性强。

程序代码
计算器操作正算说明记住要输入正算1
C匝道计算结果。