工程测量5800程序

5800通用线路测量程序及编程原理关于坐标计算内核1、程序采用数据保存在扩充存储的方式，更换或补充数据不用修改程序。

2、线路采用线路号来区分，程序内部和使用时都用自然里程参加计算。

3、平曲线不用区分匝道和主线，所有平曲线都转化成交点资料，此交点可以缓和曲线不对称、不完整，主线通常由交点起算、匝道均由曲线起点起算。

每条平曲线交点格式为9个数据。

4、直线、圆曲线、缓和曲线分别用数学模型不同的计算子程序，保证了计算速度。

直线、圆曲线子程序为精确公式，缓和曲线计算子程序为原创的复化段数可调的复化高斯2节点公式，计算精度可自由控制。

（现暂改为固定4次复化）************************************************************************************** *************数据库编制规则基本规则：圆曲线半径右转为正，坡度上坡为正，角度单位为度，边桩转角右转为正,线路号为非负整数,数据库中的直圆里程应该在起点和缓圆各输入一个互差0.0001的里程。

圆直点同样处理。

一、目录部分：格式为：起点里程、终点里程、平曲线资料起始位置、竖曲线资料起始位置、坡度资料起始位置起点里程、终点里程、平曲线资料起始位置、竖曲线资料起始位置、坡度资料起始位置……各条线路依次连续排列。

二、平曲线部分：格式为：圆半径起点里程、缓圆里程、圆缓里程、终点里程、切线长、交点X、交点Y、切线方位角、圆半径起点里程、缓圆里程、圆缓里程、终点里程、切线长、交点X、交点Y、切线方位角、圆半径……终点里程、圆半径各交点依次连续排列。

主线按图纸设计数据整理，匝道的线元数据组合成交点组合时，每条曲线由直线、缓和曲线、圆曲线、缓和曲线、直线组成，对不完整缓曲要对数据进行简单处理，不完整缓曲可以划归任一侧，不完整缓曲所在的起（终）点里程改为负，交点坐标输入为起点坐标，切线长输入0即可。

5800全线直曲计算程序通用经实践在5800上可行，反算运行时间就是慢点有4秒左右才出结果。

5800测量程序加入了数据库，综合了网上的一些而成。

非本人原创，希望能与广大朋友交流。

5800坐标正反程序(加入了数据库) 1. 主程序(TYQXJS)Lbl 4："1.SZ => XY"："2.XY => SZ"：？N：？S：Prog“SUB0”↙1÷P→C: (P-R)÷(2HPR) →D:180÷∏→E:N=1 => Goto1：Goto2：↙Lbl 1:? Z：Abs(S-O) →W：Prog "SUB1"："XS="：X◢"YS="：Y◢F-90→F：“FS=”：F▲DMS◢Goto4↙Lbl 2：？X：？Y：X→I：Y→J：Prog“SUB2”：O+W→S：“S=”：S◢“Z=”：Z◢Goto4↙2. 正算子程序(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+LWD))+ Bcos(G+QEFW(C+FWD))+Acos(G+QEMW(C+MWD)))→X：V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+B sin(G+QEFW(C+FWD))+Asin(G+QEMW(C+MWD))) →Y：G+QEW(C+WD)+90→F：X+Zcos（F）→X：Y+Zsin（F）→Y2. 反算子程序(SUB2)G-90→T：Abs((Y-V)cos(T)-(X-U)sin(T)) →W：0→Z：Lbl 0：Prog "SUB1"：T+QEW(C+WD) →L：(J-Y)cos（L）-(I-X)sin（L）→Z：IF Abs （Z）<0.000001：Then0→Z:Prog "SUB1" :(J-Y)÷sin（F）→Z Else W+Z→W:Goto4:IfEndSUB0 数据库子程序Goto1↙同时保存多个曲线时的指针Lbl 1：IF S<***（线元终点里程）:Then***→G（线元起点方位角）:***→O（线元起点里程）:***→U（线元起点X）:***→V（线元起点Y）:***→P（线元起点曲率半径）:***→R（线元终点曲率半径）: ***→H（线元起点至终点长度）：0或1、-1→Q：Return:IfEnd↙Lbl 1：IF S<***（线元终点里程）:Then***→G（线元起点方位角）:***→O（线元起点里程）:***→U（线元起点X）:***→V（线元起点Y）:***→P（线元起点曲率半径）:***→R（线元终点曲率半径）: ***→H（线元起点至终点长度）：0或1、-1→Q：Return:IfEnd ……………..为了便于解读，每增加一个线元增加一行语句，每增加一条曲线增加一个Lbl，每增加一个工程增加一个文件。

FX5800计算器公路测量常用程序集2.2 版一、程序功能主要功能：采用线元法与交点法相结合计算多条线路坐标正反算，可算任意复杂线型及立交匝道，包括C型，S型、卵型、回头曲线等；极坐标放样，全线路基边坡开挖口及坡脚放样计算、路基任意点设计高程、横坡、设计半幅宽度.隧道欠超挖放样计算。

新版本优化：1、调整一次显示结果；2、交点法中考虑了不对称缓和曲线；3、修改原版本部分地方笔误.(红色为修改处)二、源程序1.总主程序（1、坐标计算放样程序2、坐标反算程序；3、高程计算查阅程序；4、路基半幅标准宽度查阅程序；5、路基边坡及开挖口放样程序；6、路基标准距离放样；7、桥梁锥坡计算放样程序；8、极坐标计算程序；9、隧道超欠挖计算程序）运行后输入1~9数子则选择1至9的程序，返回时，在桩号输入-1,返回选择选择计算类型。

输入-2,返回选择线路。

坐标计算中输入-3,则显示本段曲线要素。

程序名：0ZCXLbl 0: 24→DimZ：Norm 2:1→A：＂A:XY=1，ZD=2 ，GC=3，GD=4， BP=5，FM=6，ZP=7，JS=8，SD=9＂?A:A=1=>Goto 1:A=2=>Goto 2:A=3=>Goto 3:A=4=>Goto 4:A=5=>Goto 5:A=6=>Goto 6:A=7=>Goto 7:A=8=>Goto 8:A=9=>Goto 9Lbl 1:Prog＂DX＂:LbI A:Prog＂QX＂:90→B:＂PJ1＂?B:B→C:＂PJ2＂?C:B→Z[1]:C→Z[8]:Lbl B:1→F:＂KM＂?Z:Z= -1=>Goto 0:Z=-2=>Goto A:Z=-3=>GotoX:Prog＂KM＂:?D:Porg＂THB＂:0→L:＂L0＂?L:Z[2]+Z[1]-Z[8] →E:X+L cos(E)→X:Y+Lsin(E) →Y:Prog＂XY＂:Prog＂JS＂:Goto BLbl 2:2→F:90→Z[1]:Prog＂QX＂: Lbl C:＂DKLC＂?Z:Z=-1=>Goto 0:Z=-2=>Goto 2:Z=-3=>Goto X:Prog＂KM＂: ＂X0＂?X: ＂Y0＂?Y:Prog＂THB＂:Porg＂ZD＂:Goto CLbI 3:Prog＂QX＂:→B:＂H-B＂?B:B→Z[9]:＂SXG＂?L:Lbl D:＂KM＂?Z:Z=-1=>Goto 0:Z=-2=>Goto 3:?D:＂ZS＂?P: Porg＂H＂:Fix 3:＂SJG=＂:Lcoate 6,4,H-Z[9]:＂SCG=＂:Lcoate 6,4,F:＂I=＂:Locate 6,4,I×100: If U<0:Then ＂W=＂: Locate 6,4, U: Else＂T=＂: Locate 6,4, U : IfEnd◢Goto DLbl 4:Prog＂QX＂:LbI E:＂KM＂?Z:Z=-1=>Goto 0:Z=-2=>Goto 4:?D:Prog＂GD＂:Fix 3:＂ SJGD=＂:Locate7,4,L: Goto ELbl 5:Prog＂QX＂:0.5→B:＂TH-GD＂?B:B→Z[23]:Lbl F:2→F:90→Z[1]:＂KM＂?Z:Z=-1 =>Goto 0:Z=-2 =>Goto 5:Prog＂KM＂:＂X0＂?X:＂Y0＂?Y:0→M:＂M0＂?M: M→Z[4]:Prog＂FBX＂:Fix 2: P→D ：Abs（D）-S→O: ＂ L0=＂:Locate 6,4,O:Prog＂ZD＂: Z[4]→T:＂ TW=＂:Lcoate 6,4,T◢ Goto FLbl 6:Prog＂DX＂:LbI G:Prog＂QX＂:LbI H:1→F:90→Z[1]:＂KM＂?Z:Z=-1=>Goto 0:Z=-2=>Goto G:Prog＂KM＂:?D:Prog＂THB＂:Prog＂XY＂:Prog＂JS＂:Prog＂H＂:0→M: ＂M0＂?M:Fix 2:H-M→T:＂ TW=＂:Locate 6,4,T◢ Goto HLbl 7:Prog＂DX＂:LbI I:Prog＂QX＂: ＂Z0＂?Z:0→D: ＂LD:Z-,Y+＂? D:Abs(D) →R:＂LR＂?R: D→Z[6]:R→Z[7]:Z→Z[8]:LbI J:0→Y: ＂L0:DZH+,XZH-＂?Y:Y=-1=>Goto 0:Y=-2=>Goto I:Prog＂ZP＂:Goto JLbl 8:Prog＂DX＂:Lbl K:?X:X=-1=>Goto 0:X=-2=>Goto 8:?Y:Prog＂JS＂:Goto KLbl 9:Porg＂DX＂:Lbl L:Prog”QX”:Lbl M:2→F:90→Z[1]: ＂KM＂?Z:Z=-1=>Goto 0:Z=-2=>Goto L:Z=-3=>Goto X: Prog＂SD1＂: 0.5→A:＂CQHD＂?A:“X0＂?X:＂Y0＂?Y:0→M:＂M0＂?M:M→Z[8]:A→Z[15]:Prog＂SD＂: Goto M: Lbl X:Mat F◢ Goto 02、次程序:路基开挖边线及填方坡脚线放样程序（输入填方放宽值、大概桩号及测量坐标、地面标高计算出偏移距离、桩号、距中距离、填挖高度）程序名：FBXProg“THB”:D→P:Prog“GD”:If D<0:Then 0.75-L→D:Goto H:Else L-0.75→D:Goto H:IfEnd:Lbl H:Prog“H”:H-0.03-Z[4] →Z[4]: Prog “W-1”:If Z[4]<0:Then –Z[4]→G:Goto W:Else Z[4] →G:Goto T:IfEnd:Lbl W: Z[10]+Z[11] +Z[12] +Z[13]→A: Z[10]+Z[11]+Z[12]→B:Z[10]+Z[11]→C: If G>A:Then Goto 1:Else If G>B:Then Goto 2:Else If G>C:Then Goto 3:Else If G>Z[10]:Then Goto 4:Else Goto 5:IfEnd:Lbl 1: L+Z[14]+Z[15]+Z[16]+Z[17]+Z[18]+(G-A)×Z[9]+Z[13]×Z[8]+Z[12]×Z[7]+Z[11]×Z[6]+Z[10]×Z[5]→S:Goto Z:Lbl 2: L+Z[14]+Z[15]+Z[16]+Z[17]+(G-B)×Z[8]+Z[12]×Z[7]+Z[11]×Z[6]+Z[10]×Z[5]→S:Goto Z:Lbl 3:L+Z[14]+Z[15]+Z[16]+(G-C)×Z[7]+Z[11]×Z[6]+Z[10]×Z[5]→S:Goto Z: Lbl 4:L+Z[14]+Z[15]+(G-Z[10])×Z[6]+Z[10]×Z[5]→S:Goto Z:Lbl 5:L+Z[14]+G×Z[5] →S:Goto Z:Lbl T:L+Z[23]→T:If G>Z[21]:Then T+Z[22] +Z[21]×Z[19]+(G-Z[21])×Z[20]→S:Goto Z:Else T+G×Z[19]→S:Goto Z: IfEnd:Lbl Z3.次程序3.隧道超欠挖值计算放样程序（输入隧道线路，大概桩号、输入衬砌厚度、测量三维坐标，计算准确桩号及位置、计算欠超挖值）程序名：SDProg “KM”:Prog“THB”: D→Z[13]:0→D:Prog“H”:Cos(Z[10])×Z[4] →E: H+Z[5] →Z[5]: H+Z[7] →Z[7]: E+Z[5] →Z[10]:If Z[8]>Z[10]:Then Goto R:Else Goto S:IfEnd:Lbl R:√((Z[13]-Z[9]) 2+(Z[8]-Z[5]) 2 )- Z[4]-Z[15] →L:Goto L:Lbl S:Z[4]-Z[6] →S:√(S 2-(Z[7]-Z[5]) 2 )×S÷Abs(S)→T:Abs(Z[13]-Z[9])-T→T:√(T 2+(Z[8]-Z[7]) 2 )- Z[6]-Z[15] →L:Goto L:Lbl L: Fix 2: “L0=”: Locate 6,4,L:Prog“ZD”:“H0=”: Locate 6,4,Z[8]-H◢4：极坐放样计算程序(计算放样点至置仪点方位角及距离)程序名：JS(一式)X：Y：Z[11]→K:Z[12]→L:Pol(X-K, Y-L):IF J<0:Then J+360→J:IFEnd:Int(J)+0.01Int(60Frac(J))+0.006Frac(60Frac(J)) →J: Fix 4: “FWJ=”: Locate 6,4,J:Fix 3:“JL=”:Locate 6,4,I◢JS(二式)X：Y：Z[11]→K:Z[12]→L:Pol(X-K, Y-L):IF J<0:Then J+360→J:IFEndFix 3:“JL=”:Locate 6,4,I: J▲DMS◢5．路基宽度子程序：路基标准半幅宽度计算程序(对于设计有加宽渐变的有用，如路基宽度无变化，则把此程序直接输入半幅宽度值至L)程序名： GD1→S: Prog “G-1”:Z-C→E: E ×(B-A)/S+A→L:6．坐标计算次程序程序名：THBIf Q>1:Then Goto J: Else 1÷P→C:(P-R)÷(2HPR)→S:180÷π→E:IfF=1:Then Abs(Z-O) →W:Prog "A":Ｇoto 2:Else X→I:Y→J:Prog ＂B＂:O＋Ｗ→Z:D→D:Ｇoto 2: IfEnd:LbI J: If F=1:Then Prog "Z":Ｇoto 1:Else Prog ＂ZX＂:Ｇoto 2: IfEnd: LbI 1:I+D×COS(Z[2]+Z[1]) →X: J+D×Sin(Z[2]+Z[1]) →Y: LbI 27. 线元法正算子程序(Ａ)程序名：A0.1184634425→A: 0.2393143352→B：0.2844444444→N：0.0→K：0.2307653449→L：0.5→M：U+W(Acos(G+QEKW(C+KWS))+Bcos(G+QELW(C+LWS))+Ncos(G+QEMW(C+MWS))+Bcos(G+Q E(1-L)W(C+(1-L)WS))+Acos(G+QE(1-K)W(C+(1-K)WS)))→X:V+W(Asin(G+QEKW(C+KWS))+Bsin(G+QELW(C+LWS))+Nsin(G+QEMW(C+MWS))+Bsin(G+Q E(1-L)W(C+(1-L)WS))+Asin(G+QE(1-K)W(C+(1-K)WS))) →Y:G+QEW(C+WS)→Z[2]:X+Ｄcos（Z[2]+Z[1]）→Ｘ: Y+Ｄsin（Z[2]+Z[1]）→Ｙ8. 线元法反算子程序(Ｂ)程序名： B G-90→T：Abs((Y-V)cos(T)-(X-U)sin(T))→W：0→D：Lbl 0：Prog "A"：T+QEW(C+WS)→L：(J-Y)cos(L)-(I-X)sin(L)→D：IF Abs(D)<0.01:Then Goto 1：Else W+D→W:Goto 0: IfEnd: Lbl 1:0→D：Prog "A"：(J-Y)÷sin(Z[2]+90) →D9. 交点法正算子程序(Z)程序名：ZH2÷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:If Z≤O -T:Then Z-O→S:G→Z[2]:Goto 2: IfEnd:If Z≤O -T+H:Then Z-O +T→S:Prog “HX”:G+WK→Z[2]:Goto 4:IfEnd:If Z≤O -T+πR×Abs(P)÷180+H÷2-N÷2: 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]:Goto 4: IfEnd:O -T+πR×Abs(P)÷180+H÷2+N÷2-Z→S:(R+N2÷R÷24-N∧(4)÷2688÷R∧(3))tan(Abs(P)÷2)+N÷2-N∧(3)÷240÷R2+E→T :N→H:Prog “HX”:G+P→S:S-WK→Z[2]:U+(T-A)Cos(S)-WBSin(S)→I:V+(T-A)Sin(S)+WBcos(S)→J:Goto 3:Lbl 4:U+(A-T)cos(G)-WBsin(G)→I:V+(A-T)Sin(G)+WBcos(G)→J: Goto 3:Lbl 2:U+Scos(Z[2])→I:V+Ssin(Z[2])→J: LbI 3:C→H10. 交点法缓和段转化子程序(HX)程序名：HXS-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→M11. 交点法反算子程序(ZX)程序名：ZXZ:0→D:Lbl 0:Prog“Z”:Pol(X-I,Y-J):J-Z[2]→J:Isin(J)→S:Icos(J) →I:If Abs(I)<0.1:Then Z+I→Z:S→D:Goto 2:Else Goto 1: LbI 1:If M=0:Then Z+I→Z:Goto 0:Eles Pol(M-WS,I):(JMπ)/180→I:Z+I→Z:Goto 0:IfEnd:LbI 212．高程计算子程序（H）程序名：HProg “S-1”:N-I→A:A÷Abs(A→W:R×Abs(A÷2→T:If Z<C-T:Then G+I(Z-C→H:Goto 0:IfEnd:If Z<C+T:Then G-TI+(I+(Z-C+T)W÷2÷R)×(Z-C+T→H:Goto 0: IfEnd:N→I: G+I(Z-C→H:Lbl 0:If D=0:Then H→H: I→I:Goto I:Else Prog “I”:H+V→H:Abs(L-13.高程超高计算程序（I）程序名：IIf Z[3]=1:Then Prog“I-1”:Goto 1:Else If Z[3]=2:Then Prog“I-2” :Goto1: Else If Z[3]=3:Then Prog“I-3” :Goto 1: Else If Z[3]=4:Then Prog“I-4” :Goto 1: Else If Z[3]=5:Then Prog“I-5” :Goto 1: Else IfZ[3]=6:Then Prog“I-6” :Goto 1: Else Prog“I-7”:Goto 1: IfEnd:Lbl 1: If W=1:Then Goto Z:Else Goto X: IfEnd :Lbl Z:If S=0:Then Abs(D)×M→V:Goto 2:Else Abs(D)×((N-M)×(Z-C)÷S+M)→V:Goto 2:IfEnd:Lbl X:If S=0:Then Abs(D)×M→V:Goto 2:Else Abs(D)×(((3((Z-C)÷S)2-2((Z-C)÷S)∧(3))×(N-M))+M)→V:Goto 2:IfEnd:Lbl 2:Abs(D)→E:V÷E→I:I(E-K)→V14．导线点子程序（DX）程序名：DXZ[11]→K:Z[12]→L:“XZ”?K:“YZ”?L:K→Z[11]:L→Z[12]15．线路选择子程序(线路选择输0时,则输曲线参数)（QX）程序名：QX1→A:“LX:XX-1,XA-2,XB-3,XC-4,XD-5,XE-6,XF-7”?A:A→Z[3]:If A=0:Then Goto 1:Else Goto 2: IfEnd:Lbl 1:“JDO”?O:“XO”?U:“YO”?V:“FW”?G:“LS”?H:?P:?R:“J=2,3;X=(0,1,-1)”?Q:If Q=2:Then H→N:Goto 2:Else If Q=3:Then“LS”?N:Lbl 216．线路中线元段判别子程序（KM）程序名：KMIf Z[3]=0:Then Goto 0:Else If Z[3]=1:Then Prog “XX-1”:Goto 0: Else If Z[3]=2:Then Prog “XA-2”:Goto 0: Else If Z[3]=3:Then Prog “XB-3”:Goto 0: Else If Z[3]=4:Then Prog “XC-4”:Goto 0: Else IfZ[3]=5:Then Prog “XD-5”:Goto 0: Else If Z[3]=6:Then Prog “XE-6”:Goto 0: Else Prog “XF-7”:Goto 0: IfEnd: LbI 017．坐标计算参数转化子程序（F）程序名：FMat 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]→P:Mat F[1,7]→R: Mat F[1,8]→Q:Q=2=>H→N:Q=3=>Mat F[1,9]→N18.桥梁锥坡放样子程序（ZP）程序名：ZP√(Z[6]2（1-Y2/Z[7]2）) →X:Z[8] →Z:Prog“GD”:If D>0:Then L+X→D:Else –L-X→D:IfEnd:Z[8]+Y→Z:D:1→F:90→Z[1]:Prog“KM”:Prog“THB”:Prog“XY”:Prog“JS”19．显示子程序（XY）和（ZD）①程序名：XYFix 3:＂XS=＂:Locate 6,4,Ｘ：＂YS=＂:Locate 6,4,Ｙ：②程序名：ZDFix 3:＂KM=＂:Locate 6,4,Z：＂ D=＂:Locate 6,4,Z：20．数据子程序(附后示例)①-1程序名：X X-1（线路1坐标计算要素程序）If Z≤26615.555:Then [[25900,11587.421,1847.983, 101。

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◢GotoA3．程序说明3.1 输入已知变量R—圆曲线半径，显示R？L—缓和曲线长，显示L？A—曲线转向角，显示A？E—直缓点纵坐标，显示ZH：X=？N—直缓点横坐标，显示ZH：Y=？F—第一切线方位角，即ZH至JD的方位角，显示FWJ=？K—直缓点里程，显示ZH：LC=？3.2 计算待求量T—切线长度，显示T= …S—曲线全长，显示S= …Z[3]—外矢距，即JD到QZ的距离，显示E0= …3.3 输入待求变量K—输入待求（测设）点的里程，显示LC=？V—横向偏距，即测设点左、右侧外移距的偏移量，若为中线点输入0；右侧输入“+”值，左侧输入“-”值。

卡西欧5800测量程序一、坐标计算源程序1.主程序(QXjs-000)Lbl 0："1.SZ => XY"："2.XY => SZ"：?N："S="?S：Prog " QXjs –SUB0"："D="?D: Prog " B2"：IfN=1：ThenGoto1：ElseGoto 2：←┘Lbl 1："S ="?S："Z ="?Z：Abs(S-O) →W：Prog "SUB1"："XS="：X◢ "YS="：Y◢"FWJ="：F-90DMS ◢Prog " B1":"H="：Z[8]◢Goto 0←┘Lbl 2："X="?X："Y="?Y： X→I：Y→J：Prog " SUB2"："S"：O+W◢"Z"：Z◢Abs (O+W)→S :Prog " B1":"H="：Z[8]◢Abs (Z-2.4)→Z: "CG="?A: "NX="?B: Prog " RRRR":"L="：L◢ Goto0←┘2.平曲线数据库子程序(QXjs-SUB0)：IfS≤***(本段线元终点里程): Then***→G(线元起点方位角)：***→O(线元起点里程)：***→U(线元起点X)：***→V(线元起点Y)：***→P(线元起点曲率半径)：***→R(线元终点曲率半径)：***→H(线元起点至终点长度)：***→Q(线元左为-1右+1)：Return:IfEnd←┘IfS≤***: Then***→G：***→O：***→U：***→V：***→P：***→R：***→H：***→H ：Return:IfEnd←┘IfS≤……带长短链的线元，输入时从长短链处分开输线元要素，3. 正算子程序(SUB1)1÷P→C：(P-R)÷(2HPR)→D：180÷π→E：0.1739274226→A：0.3260725774→B：0.0694318442→K：0.3300094785→L：1-L→F：1-K→M：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（想计算任意偏角时，可以将A值直接写成数值，再最后一步里F 后面加-90+A，就可以了，不过反算A值要输90，左偏-角右+角)4. 反算子程序(SUB2)G-90→T： Abs((Y-V)cos(T)-(X-U)sin(T)) →W： 0→Z←┘Lbl 0：Prog "SUB1"： T+QEW(C+WD) →L： (J-Y)cos(L)-(I-X)sin(L)→Z：IfAbsZ<10∧(-6) ：ThenGoto1：Else W+Z→W：IfEnd：Goto 0←┘Lbl 1：0→Z：Prog "SUB1"： (J-Y)÷sin(F )→Z←┘5高程数据库（B2）IfD=**: Then ***→Z[1](竖曲线变坡点里程) :***→Z[2](交点高程)：***→Z[3](前段坡率)：***→Z[4](后段坡率)：***→Z[5](切线长): ***→Z[6](竖曲线半径)：***→Z[7](圆弧在上面为-1，圆弧在下面为+1): ***→Z[9] (本段线元起点里程): ***→Z[10](本段线元终点里程)：6 高程计算程序(B1)IfS＜Z[9]:Then Cls:Stop:IfEnd←┘IfS≤Z[1]-Z[5]: Then Z[2]-(Z[1]-S)Z[3]→Z[8] :Return:IfEnd←┘IfS≤Z[1]:Then Z[2]-(Z[1]-S)Z[3]–Z[7](Z[5]-Z[1]+S)^2÷(2Z[6])→Z[8]: Return:IfEnd←┘IfS≤Z[1]+Z[5]:Then Z[2]+(S-Z[1])Z[4]- Z[7](Z[5]+Z[1]-S)^2÷(2Z[6])→Z[8]: Return:IfEnd←┘IfS＞Z[1]+Z[5]: Then Z[2]+(S-Z[1])Z[4]→Z[8] :IfEnd←┘IfS＞Z[10]:Then Cls:Stop:IfEnd←┘7半径计算（RRRR）围岩类型IfB=3: Then IfA≥Z[8]+1.3：Then6.98-√￣（（Z-0.00）^2+(A-Z[8]-2.27)^2）→L: Return:IfEnd←┘IfA≥Z[8]-0.78: Then√￣(6.98^2-0.97^2)-0.28(Z[8]+1.3-A)÷2.08-Z→L: Return:IfEnd←┘IfZ≥4.827: Then2.91-√￣((Z-3.711)^2+(Z[8]+1.0715-A)^2)→L: Return:IfEnd←┘IfZ＜ 4.827: Then15.27-√￣((Z[8]+12.87-A)^2+Z^2) →L: Return:IfEnd: Return:IfEnd←┘IfB=4: Then IfA≥Z[8]+1.3：Then7.21-√￣（（Z-0.00）^2+(A-Z[8]-2.27)^2）→L: Return:IfEnd←┘IfA≥Z[8]-0.78: Then√￣(7.21^2-0.97^2)-0.28(Z[8]+1.3-A)÷2.08-Z→L: Return:IfEnd←┘IfZ≥4.875: Then3.06-√￣((Z-3.711)^2+(Z[8]+1.0715-A)^2)→L: Return:IfEnd←┘IfZ＜ 4.875: Then15.42-√￣((Z[8]+12.87-A)^2+Z^2) →L:Return:IfEnd: Return:IfEnd←┘IfB=5: Then IfA≥Z[8]+1.3：Then7.39-√￣（（Z-0.00）^2+(A-Z[8]-2.27)^2）→L: Return:IfEnd←┘IfA≥Z[8]-0.78: Then√￣(7.39^2-0.97^2)-0.28(Z[8]+1.3-A)÷2.08-Z→L: Return:IfEnd←┘IfZ≥5.042: Then3.49-√￣((Z-3.711)^2+(Z[8]+1.0715-A)^2)→L: Return:IfEnd←┘IfZ＜ 5.042: Then15.85-√￣((Z[8]+12.87-A)^2+Z^2) →L: Return:IfEnd: Return:IfEnd←┘IfB=2: Then IfA≥Z[8]+0.3: Then6.46-√￣（（Z-0.00）^2+(A-Z[8]-2.27)^2）→L: Return:IfEnd←┘IfA≤Z[8]+0.3: Then2.3-√￣((Z-3.958)^2+(Z[8]+0.9916-A)^2)→L: Return:IfEnd: Return:IfEnd←┘IfB=1: Then IfA≥Z[8]+5.902:Then7.66-√￣（（Z+0.00）^2+(A-Z[8]-2.07)^2）→L: Return:IfEnd←┘IfA≥Z[8]+1.3: Then7.26-√￣((Z-0.35)^2+( A-Z[8]-2.27)^2）→L: Return:IfEnd←┘IfA≥Z[8]-0.78:Then√￣(7.26^2-0.97^2)-0.285(Z[8]+1.3-A)÷2.08-Z+0.35→L: Return:IfEnd←┘二说明（1）以线路前进方向区分左右，当曲线在左时P、R取负值，当曲线半径在右时P、R取正值，当曲线半径无穷大（即直线）时，P、R以10的45次代替。

卡西欧5800程序（完整版）Lbl 31If A>P-T+L:Then P-T+S-B-I^Z: Else P-T+B+I-Z:IfEndlIf K=l:Then If A<P-T+S-L:Then D-E -> B: Else E-D-B:IfEnd:IfEndlIf K=-l:Then If A<P-T+S-L:Then E・D - B: Else D ・E- B:IfEnd:IfEndFSHHQX1Z[5]->X:Z[6]->Y:Z[l]+90K->JlProg " HHQX u IX—Z[3]:Y—Z[4]lRYJ+2Fr+90K—JlProg " XZBH " IX—G:Y—Z1Z[3]->X:Z[4]->Y IProg、' DOTLINE u JDZBProg ” QDZB IZ[5]-Z[3]:Z[6]-Z[4] IIf O=-l:Then P-T+S—A: Prog n QDZB n :Z[3]^Z[5]:Z[4]^ Z[6]: Else Z[l]+90K+K(L^(6R))r^J:L-L A34-(9eR2)^C: Prog ”XZBH ” "End I(Z[6]-Y+Xtan(Z[2])-Z[5]tan(Z[l]))^(tan(Z[2])-tan(Z[l]))->G ItanZ[2](G-X)+Y->ZlProg n FWJ HZFXZ If M=6:Then Prog " ZSPD " : Else Prog n FSTJ n :IfEndFSLINE Prog ” QDZB IX—G:Y—Z:J—NlProg M FWJ H IN—J:F-J—FlIcos(F)-GlIf G<D:Then If A>P-T+S-L:Then P・T+S+G—Z:Isin(F)-B:Else P-T-G^Z:-Isin(F)^B:IfEnd :1-A:2-O: Else 0-A:IfEndXZBHX+Ccos(J)^XlY+Csin(J)-YHENG PO If A<P-T+L-Z[10] :Then Z[8]->V:Gotol:IfEndlE:Z[8] + （A-P+T-L+Z[10]）E->V: Goto l:IfEndlIf A<P-T+S-L:Then Z[9]^V:Goto l:IfEndlIf A<P-T+S-L+Z[10]:Then （Z[8]-Z[9]）vZ[10]->E:Z[9] + （A-P+T-S+L）E->V: Goto l:IfEndlZ[8]->V ILbl 1GAO CHENG……（交点桩号）-0：……（交点高程）-D:……（小桩号坡度）-N:……（切线长度）-G:……（圆曲线半径）-R I If A<O-G:Then 0^K: Goto l:IfEndlIf A<O+G:Then 1-K: Goto l:IfEndl（K=0为直线、K=1为凸圆、K二・l为凹圆）……（同上）Lbl l:Prog " HENG PO " IK （ A6G ） 2+（2R）-P ID+N（A-O）-P+WV^Fl FWJJS （测回法测坐标）A 1 XA u:B ” YA n :C H XB ':D " YB 1: K H CHANG DU ” :G11JIAO DU 11I=0:J=0:Pol（C-A,D-B）:J<0 J=J + 360zlE=J+GlJ " FWJ 1X"X " =A+KcosE，Y ' Y n二B+KsinE，”END M说明：H JIAO DU ”顺时钟为正、逆时针为负,“ XA “、' YA ”为测站点坐标.”XB“. ”YB”为后视点坐标，“X”. “Y”为前视即待测点坐标，" CHANG DU ”为测站点到待测点距离。

主程序JD—NELBI 1："JDk="？C："JDμ="？D："Li"？L："R="？R："JDX="？U："JDY="？V："JDA="？A："（+1R，-1L）"？M："Z（+R，—L）"？B："K="？P：Prog："JD—NE1"："X="？X"▲Y="？Y▲"F="？F▲Goto 1↙子程序JD—NE1（（R+L²÷（24R））tan（D÷2）+L÷2—L³÷（240R²）→T↙C—T→I↙I+L→J↙J+πDR÷180—L→H↙H+L→K↙If P＜I：Then U+（C-P）cos（A+180°）+Bcos（A+90°）→X↙V+（C-P）sin（A+180°）+Bsin（A+90°）→Y↙A→F：Return：IfEnd↙If P＞I And P＜J ：Then A+M×30°×（P—I）²÷（πRL）→O↙（P—I）—（P—I）²÷（90R²）→G↙U+Tcos（A+180°）+Gcos（O）→X↙V+Tsin（A+180°）+Gsin（O）→Y↙90°×（P—I）²÷（πRL）→W↙X+Bcos（A+MW+90°）→X↙Y+Bsin（A+MW+90°）→Y↙A+MW→F：Return：IfEnd↙If P＞J And P＜H ：Then A+M×90°×（L+P—J）÷（πR）→O↙2Rsin（（P—J）×90°÷（πR））→G↙U+Tcos（A+180°）+（L—L³÷（90R²））cos（A+M（30°×L）÷（πR））+Gcos（O）→X↙V+Tsin（A+180°）+（L—L³÷（90R²））sin（A+M（30°×L）÷（πR））+Gcsin（O）→X↙90°×（P—J）÷（πR）→W↙X+Bcos（O+MW+90°）→X↙Y+Bsin（O+MW+90°）→Y↙O+MW→F：Return：IfEnd↙If P＞HAnd P＜K ：Then A+MD+180°—M×30°×（P —K）²÷（πRL）→O↙K—P—（K—P）³÷（90R²）→G↙U+Tcos（A+MD）+Gcos（O）→X↙V+Tsin（A+MD）+Gsin（O）→Y↙（P—K）²×90°÷（πRL）→W↙X+Bcos（A+MD—MW+90°）→X↙Y+Bsin（A+MD—MW+90°）→Y↙A+MD—MW→F：Return：IfEnd↙If P＞K ：Then U+（T+P—K）cos（A+MD）+Bcos （A+MD+90°）→X↙V+（T+P—K）sin（A+MD）+Bsin（A+MD+90°）→Y↙A+MD→F：Return：IfEnd↙变量说明：JDk——交点桩号CJDμ——交点转角DLi——缓和曲线长度LR——圆曲线半径RJDX——交点X坐标UJDY——交点Y坐标VJDA——交点方位角A（+1R，—1L）——偏向MZ（+R，—L）——（边桩到中桩距离）BK——待求桩号PX——待求点X坐标Y——待求点Y坐标F——待求点的方位角本程序已经测试通过，第二缓和曲线起点附近有时会有1公分误差，是公式简化过程中产生的。

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

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

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

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

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

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

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

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

一．使用说明：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 22.正算子程序(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+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) →F:X+Zcos(F+T) →X:Y+Zsin(F+T) →Y3.反算子程序(SUB2)G-90→T:Abs((Y-V)cos(T)-(X-U)sin(T)) →W:0→Z:Lbl 0:Prog ″SUB1″:T+QEW(C+WD) →L:(J-Y)cos(L)-(I-X)sin(L)→Z:IF Abs(Z)<1E-6 :Then Goto 1:Else W+Z→W:Goto 0 : IfEnd:Lbl 1:0→Z:Prog ″SUB1″:(J-Y)÷sin(F)→Z 三．用于提前将曲线元要素输入子程序，计算时只输入里程桩号与边距。

公路测量坐标计算主程序（正算）一、SZ-XY（主程序）Lbl 0:4→Dim Z：ClsLbl 1:Cls:?SZ340000≤S And S＜341480.882Goto A341480.882≤S And S＜341900.82Goto B341900.82≤S And S＜343465.805Goto C343465.805≤S And S＜343885.805Goto D343885.805≤S And S＜344240.152Goto E344240.152≤S And S＜344590.152Goto F344590.152≤S And S＜344913.770Goto G344913.770≤S And S＜345263.770Goto H345263.770≤S And S＜349500Goto I349500≤S And S＜349690.625Goto J349690.625≤S And S＜350080.625Goto K350080.625≤S And S＜351523.359Goto L351523.359≤S And S＜351913.359Goto M351913.359≤S And S＜352048Goto NLbl A:Prog“P1”Prog“M1”Goto 0Lbl B:Prog“P2”Prog“M1”Goto 0Lbl C:Prog“P3”Prog“M1”Goto 0Lbl D:Prog“P4”Prog“M1”Goto 0Lbl E:Prog“P5”Prog“M1”Goto 0Lbl F:Prog“P6”Prog“M1”Goto 0Lbl G:Prog“P7”Prog“M1”Goto 0Lbl H:Prog“P8”Prog“M1”Goto 0Lbl I:Prog“P9”Prog“M1”Goto 0Lbl J:Prog“P10”Prog“M1”Goto 0Lbl K:Prog“P11”Prog“M1”Goto 0Lbl L:Prog“P12”Prog“M1”Goto 0Lbl M:Prog“P13”Prog“M1”Goto 0Lbl N:Prog“P14”Prog“M1”Goto 0M1（子程序）1/P→C:(P-R)/(2HPR) →D:180÷π→E Abs(S-O) →WProg “SUB1”:Cls“XS=”:X⊿“YS=”:Y⊿“FS=”:F-90►DMS⊿公路测量坐标计算主程序（反算）二、XY-SZ（主程序）Lbl 0:4→Dim Z：ClsLbl 2:Cls:?XY：X→I：Y→J“BH1-14”？→NN=1Prog“P1”N=2Prog“P2”N=3Prog“P3”N=4Prog“P4”N=5Prog“P5”N=6Prog“P6”N=7Prog“P7”N=8Prog“P8”N=9Prog“P9”N=10Prog“P10”N=11Prog“P11”N=12Prog“P12”N=13Prog“P13”N=14Prog“P14”Prog “SUB2”：O+W→S：Cls“S=”：S⊿“Z=”: Z⊿Goto 0公路线路主点要素数据子程序P1（原始数据，已修改）548503.126→U: 506857.907→V: 340000→O: 123°1’4.19’→G:1480.882→H:1*1045→P1*1045→R:0→QP2（原始数据）U=547696.193:V=508099.628:O=341480.882:G=123°1’4.19’:H=420:P=1*1045R=5500:Q=-1P3（原始数据）U=547471.850:V=508454.659:O=341900.882:G=120°49’48.62’:H=1564.923:P=5500：R=5500:Q=-1P4（原始数据）U=546870.501:V=509893.716:O=343465.805:G=104°31’39.79’:H=420:P=5500：R=1*1045:Q=-1 P5（原始数据）U=546775.533:V=510302.810:O=343885.805:G=102°20’24.18’:H=354.347:P=1*1045R=1*1045:Q=0P6（原始数据）U=546699.804:V=510648.971:O=344240.152:G=102°20’24.18’:H=350:P=1*1045R=7000:Q=-1P7（原始数据）U=546627.858:V=510991.487:O=344590.152:G=100°54’27.56’:H=323.618:P=7000R=7000:Q=-1P8（原始数据）U=546573.987:V=511310.561:O=344913.770:G=98°15’31.70’:H=350:P=7000R=1*1045:Q=-1P9（原始数据）U=546529.492:V=511657.712:O=345263.770:G=96°49’35.11’:H=4236.23:P=1*1045R=1*1045:Q=0P10（原始数据）U=546066.740:V=515963.181:O=349500:G=96°58’58.04’:H=190.625:P=1*1045R=1*1045:Q=0P11（原始数据）U=546043.566:V=475046.967:O=349690.625:G=96°58’58.04’:H=390:P=1*1045R=6000:Q=-1P12（原始数据）U=546000.351:V=475434.547:O=340000:G=95°7’14.43’:H=1442.734:P=6000R=6000:Q=-1P13（原始数据）U=546044.753:V=476873.123:O=351523.359:G=81°20’36.89’:H=390:P=6000R=1*1045:Q=-1P14（原始数据）U=546111.787:V=477257.300:O=351913.359:G=79°28’53.31’:H=524.641:P=1*1045R=1*1045:Q=0公路线路竖曲线计算程序GMGC(轨面高程)Lbl 1：“S”→？S340000≤S＜34022975.622+（S-340000）*0.0058→H340229≤S＜34051977.791-（S-340374）^2/(2*25000)-25000*0.0058*0.0058/2→H 340519≤S＜346018.7576.95-(S-340519)*0.0058→H346018.75≤S＜346101.2544.812+((S-346060)-25000*(0.0058+0.0025)/2)^2/(2*25000)- 25000*0.0058*0.0025/2→H346101.25≤S＜346628.7544.709-(S-346101.25)*0.0025→H346628.75≤S＜346691.2543.312+((S-346660)-25000*0.0025/2)^2/(2*25000)→H 346691.25≤S＜348882.543.312→H348882.5≤S＜348977.543.312-((S-348930)-25000*0.0038/2)^2/(2*25000)→H 348977.5≤S＜350878.7543.132-(S-348977.5)*0.0038→H350878.75≤S＜351121.2535.446+((S-351000)+25000*(0.0059-0.0038)/2)^2/(2*25000)+ 25000*0.0038*0.0059/2→H351121.25≤S＜35200036.161+(S-351121.25)*0.0059→H352000≤S＜35300041.346+(S-352000)*0.0050→HH⊿:Goto 1三、FX-5800P万能程序(TYQXJS)LB1 0:4→Dim Z:Cls"1.SZ→XY""2.XY→SZ"?→N:Cls"X0"?U:"Y0"?V:"S0"?O:"F0"?G:"LS"?H:"R0"？P："RN"?R:?Q1÷P→C:(P-R)÷(2HPR)→D:180÷π→EN=2Goto 2Lb1 1:Cls:?SZ:Abs(S-O)→WProg"SUB1":Cls"XS=":X◢"YS="Y◢"FS":F-90►DMS◢Goto 1Lb1 2:Cls:?XY:X→I:Y→JProg"SUB2":O+W→S:Cls"S=":S◢"Z=":Z◢Goto 2四、方位角、距离计算程序（P-P）Deg:Fix 3:"XO="?A:"YO="?BLb1 0:"XN=(<0→END)"?CWhile C>0"YN="?DPol(C-A,D-B):ClsIf J<0:Then J+360→F:Else J→F:IfEnd"DISTO→N(M)=":I◢"BEARO→N(DMS)=":F►DMS◢Goto 0:WhileEnd"END"五、子程序①SUB1(正算子程序)1÷P→ C： (P-R)÷(2HPR) → D： 180÷π→ E： 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+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： Abs((Y-V)cos(T)-(X-U)sin（T）) → W： 0→ Z：Lbl 0：Prog "SUB1-ZS"：T+QEW(C+WD)→L： (J-Y)cos(L)-(I-X)sin(L) → Z：If Abs(Z)<1E-6(注：此处1E-6可输入0.000001)：then Goto1：Else W+Z → W：Goto 0：IfEndLbl 1： 0→ Z：Prog "SUB1-ZS"： (J-Y)÷sin(F)→ Z③SUB3（高程子程序）O-S→L:A-B→W:Abs(R*W÷2) →T:O-T→M:O+T→PIf S≤M:ThenG-L*B→H:Goto5:Else If S≤O:Then Goto3:Else If S≤P:ThenGoto4:IfEnd: IfEnd: IfEndLbi3If W＞O：Then G+(M-S)2÷2÷R-L*B→H:Goto5: Else If W＜O: ThenG-(M-S)2÷2÷R-L*B→H:Goto5: IfEnd: IfEndLbi4If W＞O：Then G+( S-P)2÷2÷R-L*A→H:Goto5: Else If W＜O: ThenG-( S-P)2÷2÷R-L*A→H:Goto5: IfEnd: IfEndLbi5H→H:Return。