5800任意里程坐标计算 支持反算里程 绝对能用

曲线任意里程中边桩坐标正反算(CASIO fx-5800P计算器)程序曲线任意里程中边桩坐标正反算(CASIO fx-5800P计算器)程序一、程序功能本程序由一个主程序(TYQXJS)和两个子程——正算子程序(SUB1)、反算子程序(SUB2)序构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。

另外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲线任意里程中边桩坐标进行正反算。

本程序也可以在CASIO fx-4500P计算器及CASIO fx-5800 P计算器上运行。

特别申明：(1). 适用于弧长小于2倍半径的各种线元坐标正反算，精度优于1mm；(2). 在引用该核心计算部分时，请注明来源。

（3）该程序可计算标准、非标准曲线，但不可以计算与线路成斜交角的坐标二、源程序主程序(QXJS)"1.SZ => XY"："2.XY => SZ"?N:"X0"?U: "Y0"?V: "S0"?O: "FWJ0"?G: "LS"?H:"R0"?P?"RS"?R: "PX "?Q：1÷P→C：(P-R)÷(2HPR)→D：180÷π→E：If N=1:Then Goto 1：Else Goto 2:IfEndLbl 1: "KS"?Z:"ES"?S：Abs(Z-O)→W：Prog "QXZS":"XS=" :X→X◢"YS=": Y→Y◢"FWJS=":(F-90)DMS◢Prog "GC"Goto1Lbl 2："XS"?X:"YS"?Y：X→I：Y→J：Prog "QXFS":O+W→Z"KS=":Z→Z◢"ES=": S→S◢"FWJS=":(F-90)DMS◢Prog "GC"Goto 21.全线计算主程序(XLJS)"1.SZ => XY"："2.XY => SZ"："1 or 2"?N: Lbl 0："KS"?Z:"ES"?S：Prog "XXX":1÷P→C:(P-R)÷(2HP R)→D:180÷π→E：If N=1:Then Goto 1: Else Goto 2:IfEndLbl 1:Abs(Z-O)→W：Prog "QXZS":"XS=" :X→X◢"YS=": Y→Y◢"FWJS=":(F-90)DMS◢Prog "GC"Goto 0Lbl 2："XS"?X:"YS"?Y：X→I：Y→J：Prog "QXFS":O+W→Z"KS=":Z→Z◢"ES=": S→S◢"FWJS=":(F-90)DMS◢Prog "GC"Goto 22. 正算子程序(QXZS)0.1739274226→A:0.3260725774→B:0.0694318442→K:0.3300094782→L:(1-L)→F:(1-K)→M:(U+W(Ac os(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Bcos(G+QEFW(C+FWD))+Acos(G+QEMW(C+MW D))))→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+Scos(F→X:Y+Ssin(F→YReturn3. 反算子程序(QXFS )G-90→T:Abs((Y-V)cos(T)-(X-U)sin(T))→W:0→S：Lbl 0：Prog " QXZS ":T+QEW(C+WD)→L:(J-Y)co s(L)-(I-X)sin(L)→S:If Abs(S)<10^(-6):Then Goto1:Else S+W→W:Goto 0:IfEndLbl 1:0→S：Prog " QXZS ":(J-Y)÷sin(F)→SReturn数据库：XXXLbl 1:If Z≤下一起点里程: Then本线路起点里程→O：本线路X→U：本线路Y→V:本线路方位角→G：本线路长度→H：本线路起点半径→P：本线路止点半径→R:±1→Q (左偏Q=-1，右偏Q=1，直线段Q=0）：Else If Z≤下一起点里程: Then本线路起点里程→O：本线路X→U：本线路Y→V:本线路方位角→G：本线路长度→H：本线路起点半径→P：本线路止点半径→R:±1→Q :IfEnd: IfEndReturn四、数据库1、线路平面数据库数据名：“XXX”可根据自己须要更改Lbl 1:If Z≤80(计算段终点里程) : Then 0→O (计算段起点里程): 62450.391→U(起点里程纵坐标): 56853.282→V (起点里程横坐标): 220°31′33″→G(计算段起点里程切线方位角): 80→H(缓和线): 10^45→P曲率半: 782.716→R(曲率半): -1→Q (线路左右偏): Else If Z≤xxx: Then0→O (计算段起点里程): 62450.391→U(起点里程纵坐标): 56853.282→V (起点里程横坐标): 220°31′33″→G(计算段起点里程切线方位角): 80→H(缓和线): 10^45→P曲率半: 782.716→R(曲率半): -1→Q: IfEnd: IfEndReturnO=80:U=62388.712:V=56802.348:G=²²²²²²²²²²²²注：数据库内数据为示意数据，编两曲线中直线段时“H=”为该段长；半径“P=”、“R=”详以上使用说明；坐标、方位角为计算段起点的。

一．使用说明：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 三．用于提前将曲线元要素输入子程序，计算时只输入里程桩号与边距。

5800公路线路计算万能程序本程序以高斯四节点法为计算内核，多条线路正反算，中桩，边桩高程计算，超高计算，超高缓和计算，加宽计算，加宽缓和计算,边坡开口线计算，挡土墙坡脚线计算，考虑了中间绿化带的影响，适用与国家高速公路至乡村四级公路计算和放样，路基路面工程可以直接得出中边桩的设计三维坐标，去掉加宽和超高影响的计算困难，边坡和坡脚线计算可直接在边坡上提取坐标带入程序，经过所有设计因素的综合，得出更改边距并显示修改偏差后的坐标.反算速度明显提高,只需要3秒.程序无错，可放心输入，另本程序可以增加隧道超欠挖计算子程序，非常方便，因为每个隧道的断面数据不一样，所以在此没有明确写出. 计算器主程序：ZHU-CHENG-XULbiO: “1, ZS=FS, 2ZS, 3FS, 4XY=>SG. 5。

” ？U: U=l = > Prog" ZS-FS” :进入公路三维程序l>2=>Porg” ZS” :进入坐标正算程序U=3 = > Porg" FS” :进入坐标反算程序U二4 = >Porg” XY =>SG” :进入大地坐标转施工坐标程序????GotoO:说明：计算器总的主程序，进入选择各种分支计算程序。

1为公路三维计算，2为普通正算，3为普通反算，4为大地坐标转施工坐标。

°。

选择错误重新选择。

此程序可以不用输入，只为给大家一个思路，可以把计算器所有程序集中到一个主程序内管理.公路三维部分正反算选择程序：ZS-FS Deg: //设置角度模式20-DimZ: //扩展变量“1LZ=>XY, 2XY=>LZ, 3BIANP0-FY” ？ U：//正反算选择，正算选1,反算选2, 坡口坡脚选3If U=l: ThenProg” ZS~XH M : If End: //进入正算循环主程序IfU=2:Then Prog” FS~ XH” : If End: //进入反算循环主程序IfU=3:ThenProg” BIANPO-XH” : If End: 〃进入边坡开挖主程序正算循环主体程序ZS-XH“1PT-2SJ”？W：〃普通计算和设计边距计算选择“XL-XZ “？ U：〃选择线路1〜NIf W=l: Then 〃运行普通计算模式边踞自由输入LbiO:” L二”?L:” Z二” ?Z: 〃输入桩号和边距Porg" ZS-XH-1 “： GotoO:IfEnd:运行正算循环子程序If W=2:Then //运行设计计算模式边距按设计输入Lbil: ” L二” ? L：〃输入桩号Porg v Oo 5LUKUAN” : //录入公路加宽后,左右侧设计边距Porg” ZS-XH-1 “：Gotol:IfEnd: //运行正算循环子程序正算循环子程序ZS-XH-1Prog” PM-XL-XZ” : 〃进入平面线路选择程序对线路线元进行选择Prog" CHA0XIAN- CL” ： //里程超限处理Prog” GL-ZS” ：〃进入坐标正算程序Prog” GL-BZ-ZS” : 〃进入公路边桩正算程序Prog” GC-XL-XZ” ：〃进入高程线路选择程序对线路线元进行选择Prog” GC” : // 进入高程部分提取中桩高程Prog\进入超高数据库提取超高数据If W=l:ThenProg" B-H” : 〃普通计算时进入边桩高程部分提取边桩高程ElseZ[3] -Z: Porg “B-H” : //按设计计算时候代入减去半边绿化带后的路基边距, 提取边桩高程Ifend:Prog” ZS-XS” : //进入正算显示程序，提取正算三维结果反算循环主体程序FS-XH0-Z: “XL-XZ” ? U：〃选择线路 1 〜N LbiO：“L” ？ L：“X二” ? T：“Y二” ？ S：〃输入起算假定桩号，实测X Y Prog” GL-FS” : /进入公路反算程序Prog” FS-XS” : 〃进入反算显示程序GotoO:坐标正算程序GL-ZS.1739274226-1:. 3260725774-J:.0694318442-*K:. 3300094782-E : //四节点法计算系数录入1-E-F： 1-K-M：Q+KrO(C+KD) -K: //数据库要用X、Y、Q作为起点坐标方位角Q+ErO(C+ED) -E://r 是弧度，输入方法FUNCTION 选 5 Q+FrO(C+FD) -*F: Q-MrO(C+MD) -*M:X+0(ICos (K) +JCos (E) + JCos (F) +ICos (M)) -»X: //计算出X 坐标Y+0(ISin (K) +JSin(E) +JSin (F) +ISin(M)) -Y: //计算出Y 坐标Q+Or (C+D) -Q：〃计算出算点方位角,0是字母,r是弧度：边桩坐标计算程序GL-BZ-ZS:X+Zcos (Q+90) -*G: 〃计算出边桩X坐标Y+ZSin(Q+90) -M: 〃计算出边桩Y坐标公路反算程序GL-FS20-V: //V值设置为20,为第一次反算控制Do：Porg” PM-XL-XZ” ：〃进入平面线路选择程序Porg” CHAOXIAN-CL” : 〃里程超限处理程序If V>10: // V越小反算速度可能越慢，但是半径过小有可能算不出来，10为建议值Then Porg" GL~JY~ZS M : //进入公路简易正算二节点程序Else Porg" GL_ZS n : // 进入公路正算4节点程序IfEnd:T-XT:S-Y-J: //试算点与测点X, Y距离差Pol(I, J):J-Q-J: //方位角转换为接近试算点到测点的值：Rec(I, J): //计算出X 差和Y差(X里程，J边)I/(1-CJ-2DJ) -V: //计算出里程大致差距(此公式为加速反算)L+V-L：〃修改后的里程Abs(V) -V: //修改值的绝对值LpWhile V>0. 001: 〃小于一亳米通过，大于一亳米循环计算中桩高程计算程序GCIf I-A<O:Thenl-*H:IfEnd: 〃开口向上部分赋值IfI-A=0：ThenGotol: If End: // 无竖曲线部分赋值IfI-A>O:Then-l-*H:IfEnd: //开口向下时赋值(L-C) I+K+ (L-C) *2/(2E)H-1: //计算得出竖曲线部分中桩高程Goto2:Lbil:K+(L-C)I-1: 〃计算得出无竖曲线部分中桩高程Lbi2:边桩高程计算程序B-H(L-C) / (D-C) ->P: //计算点占本线元长比例3P*(2* (S))-2P*(37S)) -P: //计算得线性和高次抛物线系数F+P (T-F) -E：// 初步计算计算点横坡If V=l:Then If N= T:Then //直线和直缓,缓直，曲线左偏时If Z<0： ThenIfE>Z[4] :ThenZ[4]-*E: If End : //计算左侧时,横坡大于一般横坡取一般横坡I+Abs ⑵ E-0: If End : 〃左侧高程IfZ>0:Then I+ZF-0: IfEnd: If end: //右侧高程得数If N=1:Then IfZ<0:Then I+Abs (Z)E-O: IfEnd：//曲线右偏时，左侧高程得得IfZ>O:Then If E>Z[4]:Then Z[4]-E:IfEnd //右侧高程横坡大于一般横坡,取一般横坡I+ZE-*O:IfEnd： If End: If End: //右侧高程得数IfV二2:ThenIf Z>O:Then I+ZE-0: IfEnd://圆曲线全超高和S曲线超高渐变部分右侧高程得数If Z<0: Thenl-Abs (Z) E^O: If End: If End: //圆曲线全超高和S曲线超高渐变部分左侧高程得数备注:如果有更多的坡级，大家可以使用增加变量.完全满足大家边坡放样需求设计标准路面宽度子程序CHANGSHU12-ZE11：〃录入12米为路面半宽度，根据设计修改1. 5->Z[21：〃录入1. 5米为绿化带半宽，根据设计修改，无绿化带输入0 -0. 02 -Z[4]: 〃录入-0。

FX5800计算器交点法正反算程序1、主程序程序名：1XY-KMLbI 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,Ｘ：＂Ｙ＝＂：Loca te 6,4,Ｙ：◢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,Ｚ：＂Ｄ＝＂：Loca te 6,4,Ｄ◢Goto C2. 交点法正算子程序(Z)程序名：ZH2÷R÷24-H∧(4)÷2688÷R∧(3)→A:H÷2-H∧(3)÷240÷R2→B:（（H2-N2）÷24÷R÷S in(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:T hen 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 :G oto 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+Ssin(Z[2]→J: Return:LbI 2:U+Acos(E)-WBsin(E→I:V+ASin(E)+WBcos(E→J3. 交点法缓和段子程序(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 →M4. 交点法反算子程序(ZX)程序名：ZX0→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时。

卡西欧5800P计算器曲线任意⾥程中边桩坐标正反算卡西欧5800P计算器曲线任意⾥程中边桩坐标正反算⼀、程序功能本程序由⼀个主程序(ZHUCHENXU)和⼏个⼦程序——正算⼦程序(SUB1)、反算⼦程序(SUB2)、数据库(SUB3)、算⽅位⾓程序(FA)、使⽤⽅位⾓算坐标(NE)、边坡放样程序(BIANPO)、隧道圆⼼放样程序（SDYX）构成，可以根据直线、圆曲线、缓和曲线（完整或⾮完整型）的线元要素（起点坐标、起点⾥程、起点切线⽅位⾓、线元长度、起点曲率半径、⽌点曲率半径）及⾥程边距或坐标，对该曲线段范围内任意⾥程中边桩坐标进⾏正反算。

⼆、源程序增加变量16→DimZ1.主程序(ZHUCHENXU)"1.SZ → XY，2.XY → SZ，3.SDYX，4. BIANPO "？N：Lbl 1：”ZHUANGHAO=”?S：Prog "SUB3"：N≥2=>Goto 2：Abs(S-O)→W：”BIANZHU=”?→Z:Prog "SUB1"："XS=”:X→X ◢"YS=”:Y→Y◢"FS=”:(F-90) →F:F►DMS◢Prog "FA"：Goto 1:Lbl 2：”CX=”?X:X→I：”CY=”?Y: Y→J：”DMG=”?→Z[7]：Prog "SUB2"："S=":(O+W) →S◢"Z="Z→Z◢If N=3:ThenProg”SDYX”:IfEnd ：If N=4:Then Prog”BIANPO”:IfEnd:Goto 22. 正算⼦程序(SUB1)0.1739274226 →Z[1]：0.3260725774→B：0.0694318442→K：0.3300094782→L：(1-L) →F：(1-K)→M：U+W(Z[1]cos(G+QEK W(C+KWD))+Bcos(G+QELW(C+LWD))+Bcos(G+QEFW (C+FWD))+Z[1]cos(G+QEMW(C+MWD)))→X：V+W(Z[1]sin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+QEFW(C+FWD))+Z[1]sin(G+QEMW(C+MWD))) →Y：G+QEW(C+WD)+90→F：X+ZcosF→X：Y+ZsinF→Y3. 反算⼦程序(SUB2)G-90 →Z[2]：Abs((Y-V)cosZ[2]-(X-U)sinZ[2]) →W：0→Z：Lbl 0：Prog "SUB1"：Z[2]+QEW(C+WD) →L：(J-Y)cosL-(I-X)sinL→Z：If Abs(Z)<1E-6:Then Goto1：Else W+Z→W：Goto 0: IfEndLbl 1：0→Z：Prog "SUB1"：(J-Y)÷sinF→Z4. 数据库(SUB3)If S<226100:Then Cls:Stop: IfEndIf S<226255.833:Then 226100→O：99037.9736→U：4282.3590→V：201.3266304→G：1045→P：1045→R：253.543→H：0→Q: Goto 1:IfEnd…………If S>236005.651:Then Cls:Stop: IfEndLbl 1: 1÷P→C：(P-R)÷(2HPR) →D：180÷π→E：5. 算⽅位⾓程序(FA)LbI0：“XC=”? A : “YC=”? T :X：Y0→J：POI（X-A，Y-T）→I：“BC=”：I◢If J＜0 Then 360+J→J：ElseJ→J :IfEnd“FW=”：J►DMS◢ Goto06. 使⽤⽅位⾓算坐标(NE)LbI0：“X=”? X: “Y=”? Y: “L=”? L:“F0”?F：L×Cos（F）+X→N: “N=”N◢L×Sin（F）+Y→E: “E=”E◢Goto07. 边坡放样程序(BIANPO)“QPDJL=”? → Z[3]: “QPDGD=”? → Z[4]: “PG=” ? → Z[5]：2→Z[6]：“P1=” ? → Z[8]：“P2=” ? → Z[9]：“P3=” ? → Z[10] (没有的话就删)：Z[7]-Z[4] → Z[11]：If Z[11]〈0 Then AbsZ[11] →Z[11]: “T=”：Z[11] ◢Else“W=”：Z[11] ◢IfEndZ[11]÷Z[5] →Z[12]：If Z[12]〈1 Then Z[8]×Z[11] →Z[13]：Else Z[8]×Z[5]+Z[9]×Z[5]×（Z[12]-1）→Z[13]：IfEnd:AbsZ-（Z[13]+Z[6]×Int Z[12]+Z[3]）→Z[14]:“BP+C-J=”: Z[14] ◢说明：Z[3]= QPDJL ………………起坡点离中桩距离Z[4]=“QPDGD” ……………起坡点设计标⾼Z[5]=“PG” ………………坡⾼Z[6]=“PTK” ……………… 平台宽度Z[7]=“DMG” ………………实测地⾯⾼度Z[8]=“P1” ………………⼀级边坡坡度Z[9]=“P2” ………………⼆级边坡坡度Z[10]=“P3” ……………… 三级边坡坡度(没有的话就删)Z[14]=“BP+C-J=” ………………边坡超⽋挖宽度（+表超，-表⽋）8. 隧道圆⼼放样程序（SDYX）IF S<229000:THEN 243.997+0.0205×(S-226520)-(S-227318)^2/2/16000+1.503 →Z[10]:“NGC=”：Z[10] ◢ELSE 286.385-0.023×(S-230500)- (S-230593.753)^2/2/165000+1.503→Z[10]:“YLJ=”：Z[10] ◢IfEndZ[10]- Z[7] →Z[15]:“+SG,-JD”：Z[15] ◢√（Z[15]^2+(ABS(Z)-6.8)^2）→Z[16]: “R=”:Z[16] ◢三、使⽤说明1、规定(1) 以道路中线的前进⽅向（即⾥程增⼤的⽅向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。

5800坐标正反算程序“1.If SZThen XY”：“2. If XYThen SZ”：“N=”？→N：起点X→U：起点Y →V：起点里程→O：起点方位角→G：曲线长度→H：起点曲率半径→P：终点曲率半径→R：（左偏－1；直线0；右偏1）→Q：1÷P→C：（P－R）÷(2HPR)→D：180÷π→E：If 1→N：Then Goto 1：Else Goto 2：IfEnd←Lbl 1：“S＝”？→S：“Z＝”？→Z：Abs(S－O)→W：Prog “SUB1”：X→X▲Y→Y▲Lbl 2：“X＝”？→X：“Y＝”？→Y：X→I：Y→J：Prog “SUB2”：O+W→S▲Z→Z▲Goto 2正算子程序SUB10.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))+A cos(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+ZcosF→X：Y+ZsinF→Y←反算子程序SUB2G－90→T：Abs((Y－V)cosT－(X－U)sinT)→W：0→Z：Lbl 0：Prog “SUB1”：T+QEW(C+WD) →L：（J－Y）cosL－(I－X)sinL→Z：If AbsZ＜0.000001：Then Goto1：Else W+Z→W：Goto0：IfEnd←Lbl 1：0→Z：Prog“SUB1”：（J－Y）÷sinF→Z5800坐标反算程序Lbl 0：“X=”?→X：“Y＝”？→Y：焦点X→N：焦点Y→E：转向角→A：半径→R：缓和曲线长→C：切线长→T：直缓点里程→G：缓直点里程→M：方位角→F：If G＜M：Then 1→P：F→U：Else －1→P：F+A→U：IfEnd←If A＞0：Then 1→Z：Else －1→Z：IfEnd←Pol(X－N，Y－E) ←G+PT+Rec(I,J－U)→S：J→D←If G＞M：Then Goto3：IfEnd←If S≤G ：Then S→K▲ZD→Q▲Goto 0←Else Goto4：IfEnd←Lbl 3←If S≥G：Then S→K▲ZD→Q▲Goto 0←Lbl 4←C2÷24÷R－C4÷2688÷R3→W：C÷2－C3÷240÷R2→V：90C÷R÷π→B：←Pol(P(S－G) －V，ZD－R－W) ←If J≥B－90：Then Goto 5：IfEnd←P(S－G)+P(S－G)5÷40÷R2÷C2→Z←Z3÷(6RC) →H←If A＞0：Then D－H→W：Else －D－H→ W：IfEnd←If RC≤2WC：Then RC÷W→Z：Else(RC－√（R2C2－2RCWZ）) ÷W→Z：IfEnd←G+PZ→K▲K－PZ5÷(40R2C2) →V：Z3÷(6RC) →H←If A＜0：Then －1→Z：Else 1→Z：IfEnd←√((S－V)2+(ZD－H)2) →W←If ZD≥H：Then W→Q▲Else －W→Q▲IfEnd←Goto 0←Lbl 5←(90+J)π÷180－C÷2R→J←If RJ≥P(M－G) －2C：Then Goto 6：IfEnd←G+P(C+JR) →K▲R－I→Q▲Goto 0←Lbl 6←Pol(P(S－G) －T，ZD) ←M－P(T－Rec(I，J－AbsA)) →O←If G＞M：Then Goto 7：IfEnd←If O＞M：Then O→K▲J→O▲Goto 0←Else Goto 8：IfEnd←Lbl 7←If O≤M：Then O→K▲J→Q▲IfEnd←Goto 0←Lbl 8←P(M－O)+P(M－O)5÷(40R2C2) →Z←Z3÷(6RC)) →H←J－H→W←If RC≤2WZ：Then RC÷W→Z：Else(RC－√(R2C2－2RCWZ)) ÷W→Z：IfEnd ←M－PZ→K▲K+PZ5÷(40R2C2) →V：Z3÷(6RC) →H：√（（O－V）2+(J－H)2）→W←If J≥H：Then W→Q▲Else －W→QIfEnd←Goto 0←。

卡西欧5800计算器坐标正反算程序主程序坐标正反算程序名称：ZBZFSLb1 A↙30→Dim Z :"0=ZS,1=FS"?Z ↙If Z=0:Then Goto B:IfEnd↙ (Z=0进入里程点坐标正算）If Z=1:Then Prog"FSLCZ "↙ (Z=1进入反算里程边桩）Lb1 B ↙Prog”ZBZS” ↙子程序反算里程桩名称：FSLCZLb1 1↙"XK="?X:"YK="?Y↙（输入任意测点的X、Y坐标）“K=”?K ↙（试算里程，任意输入标段内里程点即可，也可以输入大致的估算里程加快速度）Lb1 2↙Prog"SJK ":T-90 →W:W<0=>360+W→W:Abs((Y-Z[19])*Cos(W)-(X-Z[18])*Sin(W))→S↙If S<0.0001:Then Goto 4:Else Goto 3:Ifend↙Lb1 3↙K+S→K: Prog"SJK":T-90 →W:W<0=>360+W→W:Abs((Y-Z[19])*Cos W-(X-Z[18])*Sin (W)→Q↙If Q<0.0001 :Then Q→S: Goto 4:Else if Q<S: Then K+Q→K:Goto 2 :Else if Q>S :Then K-Q→K:Goto 2:Ifend:Ifend:Ifend↙Lb1 4↙Pol (X-Z[18],Y-Z[19]:"DP(-Z+Y)=":I◢(偏距）"K=":K+S→K◢（里程）Goto 1↙子程序坐标正算名称：ZBZSLb1 0 ↙“XHS="?G(后视点X):"YHS="?L(后视点Y):"XZJ="?M(置镜点X):"YZJ="?N(置镜点Y):Pol(G-M,L-N):"DH=":I(后视距)◢J<0=>J+360→J:"FH=":JDMS◢(后视方位角)Lb1 1↙“K=”?K ：（输入所需计算里程）Prog"SJK "↙XI : Z[18]↙YI: Z[19] ↙Pol(Z[18]-M,Z[19]-N): J<0=>J+360→J↙“PJ=”?P↙ (输入桩与线路夹角)“PD=”?D↙ （输入桩距中线的距离）Z[18]+D*Cos(T+P) →Z[20] ↙Z[19]+D*Sin(T+P) →Z[21] ↙“X=”: Z[20] ◢ （放样坐标X）“Y=”: Z[21] ◢ （放样坐标Y）Pol(Z[20]-M,Z[21]-N):"D=":I◢（放样距）J<0=>J+360→J:"F=":JDMS◢（放样方位角）Goto 1↙子程序数据库名称：SJKif K<本段曲线终点里程And K≥上段曲线终点里程：Then 本段曲线终点里程→Z[1] : 上段曲线终点里程（第一段曲线输起点的里程）→Z[2] ：1→O （注：左偏曲线输入-1→O,右偏曲线输入1→O）: 半径→R :曲线偏角→A:第一缓和曲线→Z[6] : 第二缓和曲线→Z[7] : 交点X→B :交点Y→C : 小里程向交点方位角→E : 交点向大里程方位角→F : Pr og”JSPB”:Return: Ifend↙if…………Prog”JSPB”:Retur n:Ifend（曲线段分段输入）↙补充直线段输入如下 (只需输线路的最后一段直线数据）if K<本段直线终点里程And K≥本段直线起点里程:Then 本段直线终点里程→Z[3]:终点X→Z[16]:终点Y→Z[17]:方位角→E: Z[16]+ (K- Z [3])*Cos(E)→Z[18]:Z[17]+ (K- Z[3])*Sin(E)→Z[19] : Return:Ifend ↙子程序计算判别名称：JSPBLb1 2 ↙（曲线要素计算）Z[6]/2- Z[6]^3/(240*R^2)+ Z[6]^5/(34560*R^4) →Z[8] ↙ （M1）Z[7]/2- Z[7]^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[6]+ Z[7])→Z[25] ↙ （曲线总长）90* Z[6]/(R*π) →Z[14] ↙ （第一缓和曲线总偏角）90* Z[7]/(R*π) →Z[15] ↙（第二缓和曲线总偏角,可以省略）Z[8]＋（R+Z[10])Tan(A/2)-(Z[10]-Z[11] )/Sin( A)→Z[12]↙ (切线T1)Z[9]＋（R+Z[11])Tan(A/2)+(Z[10]-Z[11] )/Sin (A）→Z[13]↙ (切线T2)B+ Z[12]*Cos (E+180)→ Z[16] ↙ （ZH点X）C+ Z[12]*Sin(E+180)→ Z[17] ↙ （ZH点Y）Z[1]- Z[25]→Z[3] ↙ (ZH点里程)Z[3]+ Z[6]→Z[4] ↙ (HY点里程)Z[1]- Z[7]→Z[5] ↙ (YH点里程)Goto 3 ↙LB1 3 ↙(判断里程点与曲线关系)if K≤Z[3] And K> Z[2] : Then Goto 4 : Ifend ↙if K≤Z[4] And K> Z[3] : Then Goto 5 : Ifend ↙if K≤Z[5] And K> Z[4] : Then Goto 6 : Ifend ↙if K≤Z[1] And K> Z[5] : Then Goto 7 : Ifend ↙Lb1 4 ↙（里程小于直缓点直线独立坐标）K- Z[3] →Z[23] : 0→Z[24] : E→T : Goto 8↙Lb1 5 ↙（第一缓和曲线独立坐标）K- Z[3] →H ↙H-H^5/(40*R^2* Z[6]^2)+H^9/(3456*R^4* Z[6]^4) →Z[23] ↙H^3/(6*R* Z[6])-H^7/(336*R^3* Z[6]^3) →Z[24] ↙90*H^2/( R*π* Z[6]) →T ↙if O>0 ：Then T +E→T : Eles E-T →T : T<0=>360+T→T : Ifend ↙Got o 8 ↙Lb1 6 ↙（圆曲线独立坐标）K- Z[4] →H ↙H*180/( R*π)+ Z[14]→T ↙R*Sin(T)+ Z[8]→Z[23] ↙R*(1-Cos(T))+ Z[10]→Z[24] ↙if O>0 ：Then T +E→T : Eles E-T →T : T<0=>360+T→T : Ifend ↙Goto 8 ↙Lb1 7 ↙（第二缓和曲线独立坐标）Z[1] -K →H ↙H-H^5/(40*R^2* Z[7]^2)+H^9/(3456*R^4* Z[7]^4) →U↙H^3/(6*R* Z[7])-H^7/(336*R^3* Z[7]^3) →V ↙90*H^2/( R*π* Z[7]) →T ↙Z[13]Cos(A)+ Z[12]-U*Cos(A)-V*Sin(A)→Z[23] ↙Z[13]*Sin(A)-U*Sin(A)+V*Cos(A)→Z[24] ↙if O>0 ：Then F-T→T : T<0=>360+T→T : Else F+T →T : Ifend ↙Goto 8 ↙Lb1 8if O<0 : Then - Z[24]→Z[24] : Ifend ↙Z[16]+Z[23]*Cos(E) -Z[24]*Sin（E）→Z[18] ↙Z[17]+Z[23]*Sin(E)+Z[24]*Cos(E)→Z[19] ↙Return↙此程序大约占计算器2900字节！此程序本人已亲自验算，无误！！子程序反算里程桩名称：FSLCZLb1 1↙"XK="?X:"YK="?Y↙（输入任意测点的X、Y坐标）“K=”?K ↙（试算里程，任意输入标段内里程点即可，也可以输入大致的估算里程加快速度）Lb1 2↙Prog"SJK ":T-90 →W:W<0=>360+W→W:Abs((Y-Z[19])*Cos(W)-(X-Z[18])* Sin(W))→S↙If S<0.0001:Then Goto 4:Else Goto 3:Ifend↙Lb1 3↙K+S→K: Prog"SJK":T-90 →W:W<0=>360+W→W:Abs((Y-Z[19])*Cos W-(X-Z [18])*Si n (W)→Q↙If Q<0.0001 :Then Q→S: Goto 4:Else if Q<S: Then K+Q →K: Goto 2 :Else if Q>S :Then K-Q→K:Goto 2:Ifend:Ifen d:Ifend↙Lb1 4↙Pol (X-Z[18],Y-Z[19]:"DP(-Z+Y)=":I◢(偏距）"K=":K+S→K◢（里程）Goto 1↙Lb1 4↙Pol (X-Z[18],Y-Z[19]:------------差个撒样"DP(-Z+Y)=":I◢(偏距）"K=":K+S→K◢（里程）Goto 1↙。

CASIO 9860/5800 坐标正反算计算程序1. 坐标正算主程序（命名为ZBZS）第1 行：Lbl 0: ” K= K: ” BIAN=” Z: ” a =” B第2 行：Prog “ A”第3 行：” X仝：N+Zcos（F+B）丄第4 行：” 丫仝:E+Zsin（F+B）丄第5行：” F=” :F丄第6 行：Goto 0K――计算点的里程BIAN --- 计算点到中桩的距离（左负右正）a ---- 取前右夹角为正2. 坐标反算桩号和偏距主程序（命名为ZBFS）第1 行：” X1=” C: ” 丫1=” D: ” K1=” K第2 行：Lbl 0:Prog“A”第3 行：Pol(C-N,D-E)第4 行：List Ans[1]T I第5 行：List Ans[2]T J第6 行：Icos(F-J) T S:K+S T K第7 行：Abs(S)>=>Goto 0第8行：” K仁”：K丄第9 行：” BIAN=” ：lsin（J-F）Z JX1——取样点的X坐标Y1——取样点的Y坐标K1――输入时为计算起始点（在线路内即可），输出时为反算点的桩号Z――偏距（左负右正）3. 计算坐标子程序（命名为XYF）为了简洁，本程序由数据库直接调用，上述中的正反算主程序不直接调用此程序第1 行：K-A T S：（Q-P）- L T I第2N+/ (cos(F+X(2P+XI) X 90-n ),0,S) T N行：第3E+/ (sin(F+X(2P+XI) X 90-n ),0,S) T E行：第4F+S(2P+SI) X 90-nT F行：4. 数据库(命名为A)第1 行：K< =>Stop第2 行A: T N: T E:92 ° 26' 40 F:0 f P:1/240 f Q: f L:KC A+L =>GoTo 1 (第一缓和曲线，圆半径为240)第3 T A: T N:T E: 100 °50'"T F: 1/240 T P:1/240 T Q:T L: K C A+L =>Goto 1 (第圆曲线，半径为240)第4 T A: T N:T E: 118 °15'"T F: 1/240 T P: 0 T Q:f L: K C A+L =>Goto 1 （第二缓和曲线，圆半径为240）第5 T A: T N:T E:124 °50'"T F:0 T P:-1/180 T Q:T L:K C A+L=>Goto 1:Stop （第一缓和曲线，圆半径为180）第6 行：Lbl 1:Prog “ XYF”A――曲线段起点的里程N 曲线段起点的x坐标E 曲线段起点的y坐标F ---- 曲线段起点的坐标方位角P――曲线段起点的曲率（左负右正）Q曲线段终点的曲率（左负右正）L――曲线段长度（尽量使用长度，为计算断链方便）说明：（1 ）在9860中，程序中所有公式和部分函数结果均存储在List Ans列表数组中,要想多次调用最好随公式取出结果，并赋给变量。

反算原理在图中，A点为已知坐标而待求对应中桩桩号及边距的点。

B点为假定的A点对应中桩桩号点。

显然，B点并不对应于A点。

做出B点的切线，过A点做辅助线垂直于B点的切线，相交于C点。

设：B点的切线方位角为α，B点的桩号为K B,B点的坐标分别为X B、Y B，A点的桩号为K A,A点的坐标分别为X A、Y A，“B-A”的方位角为β，“B-A”的距离为N，“B-C”的距离为L，“C-A”的距离为Z。

根据前面的坐标正算的公式可以得到α，X B、Y B值。

根据计算器内置的Pol（X A-X B，Y A-Y B）公式（直角坐标转换为极坐标）能得到β，N值。

（1）（2）当L=0时，B点是对应于A点的，K B=K A，Z即为A点的距中桩的距离。

当L≠0时，则采用K B=K B+L,对B点进行新的假定，进而再次对L进行解算，直至L=0，或则L值在容许误差范围之内。

坐标正反算通用程序(终极篇)1. 坐标正算主程序(命名为KP-XY)第1行：Lbl 0“DK=”?K:“PJ=”? P:“α=”?W第2行：Prog “GP”(GP为线路名称)第3行：“X=”:X+Pcos(F+W)→X◢第4行：“Y=”:Y+Psin(F+W)→Y◢第5行：“FWJ=”:F►DMS◢第6行： Goto 0K——计算点的里程PJ——计算点到中桩的距离（左负右正）α——取线路前进方向右夹角为正X、Y--为计算坐标、F为方位角GP--为线路平曲线名称2. 坐标反算桩号和偏距主程序(命名为XY-KP)第1行：“DK=”?K:“X:=”? U:“Y:=”?V第2行：Lbl 0Prog “GP”第3行：Pol(U-X,V-Y):Icos(F-J)→Z[1]:K+Z[1]→K 第4行：Abs(Z[1])>0.0001=>Goto 0第5行：“DK=”:K◢第6行：“PJ=”:Isin(J-F)→P◢X:= —取样点的X坐标Y:= —取样点的Y坐标DK= 输入时为计算起始点(在线路内即可)，输出时为反算点的桩号P—偏距（左负右正）3. 计算坐标子程序(命名为XYF)为了简洁，本程序由数据库直接调用，上述中的正反算主程序不直接调用此程序第1行：K-G→S：(D-C)÷L→I第2行：X+∫(cos(F+X(2C+XI)×90÷π),0,S)→X第3行：Y+∫(sin(F+X(2C+XI)×90÷π),0,S)→Y第4行：F+S(2C+SI)×90÷π→F第5行：F<0=>F+360→F: F>360=>F-360→F4. 数据库（命名为A）第1行：0→G:3925525.975→X:502796.176→Y:3°47′30.2″→F:1/996→C:1/180→K≤G+L =>GoTo1（第一缓和曲线）第2行：245.607→G: 3925566.736→X:502801.006→Y:11°30′39.4″→F: 1/180→C:1/180→D:91.027→L:K≤G+L =>Goto1（圆曲线）第3行：318.522→G: 3925647.682→X:502804.484→X:40°29′08.7″→F: 1/180→C: 1/295.54→D: 70.37→L:K≤G+L =>Goto1（第二缓和曲线）第4行：Lbl 1:Prog “XYF”G——曲线段起点的里程X——曲线段起点的x坐标Y——曲线段起点的y坐标F——曲线段起点的坐标方位角C——曲线段起点的曲率(半径倒数，直线为0，左负右正)D——曲线段终点的曲率(半径倒数，直线为0，左负右正)L——曲线段长度（尽量使用长度，为计算断链方便）说明：（1）适用于任意线形：直线（0→C、0→D）、圆曲线（圆半径倒数→C、圆半径倒数→D）、缓和曲线（0或圆半径倒数→C、圆半径倒数或0→D）、卵形曲线（接起点圆的半径倒数→C、接终点圆的半径倒数→D），曲线左转多加一负号。

CASIO 5800计算器程序一、程序功能：本程序由一个主程序（AK）及两个子程序：坐标正算程序（SUB1）、坐标反算程序（SUB2）平曲线数据库（PQX）构成，可以根据曲线段的组成部分：直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、止点坐标、起点切线方位角、线元长度、起点曲率半径、止点曲率半径、里程、边距或坐标），对该曲线段范围内任意里程中边桩坐标进行正反算。

二、源程序：1、主程序（AK）A排在前面，方便运行程序，不用往下翻页了。

Deg：“0→ZX,1→A,2→B”?Z ：“1→JDF,2→XYF”?E：“X(CZ)”?A：“Y”?B 说明：设置角度模式，选择线路0为主线，1为A匝道。

可以自己再增加。

选择计算方法，交点法还是线元法计算坐标，1为交点法，2为线元法。

输入测站坐标，如不需要计算极坐标，即放样的时候需要输入坐标到仪器里面去，则X(CZ)输0，Y不用管可以为任何值，建议把测站坐标输进去，省了放一个点输一次坐标数据到仪器里面去。

只有一条线路可以不需要黄色部分。

“1.Ki→XY,2.XY→Ki”?N：Lbl 0：“Ki”?K：Prog“PQX”IfE=1：Then Prog“JD”：IfEndN=1＝＞Goto 1：N=2＝＞Goto 2Lbl 1?LK-O→W：Prog“SUB1”说明：1为正算，即知道桩号和偏距（注意只能计算边桩，不能计算斜交的点，构造物可以通过后面的ZHIXIAN程序计算），2为反算，即知道坐标，求桩号和偏距。

Ki为桩号，正算直接输，反算需要输入附近桩号，第一次计算最好输好桩号之后，把反算出来的桩号再输进去，重新反算下，如果和第一次反算的桩号一样，那就说明反算的桩号正确。

调用平曲线数据库，输入偏距L（偏距左输负数，右输正数），如果选择的是交点法计算，那需要输入是数据库是交点数据，调用正算子程序SUB1。

If A≠0：Then Pol（X-A,Y-B）：“JULI=”：I◢J<0=>J+360→J：“FWJ=”：J ▶DMS◢Else “Xi=”：X◢“Yi=”：Y◢“FWJ=”：F ▶DMS◢If End ：Goto0说明：正算，如果前面X(CZ)输0，则直接显示X、Y坐标和方位角，否则显示仪器到放样点的距离和方位角。

(以下程序是专业人士编写，本店铺不对程序负责，仅供您参考使用。

)卡西欧fx5800p计算器坐标正反算程序一、程序功能本程序由 6 个主程序、 5 个次子程序及 5 个参数子程序组成。

主要用于公路测量中坐标正反算，设计任意点高程及横坡计算 , 桥涵放样，路基开挖口及填方坡脚线放样。

程序坐标计算适应于任何线型 .二、源程序1. 主程序 1 ：一般放样反算程序(① 正算坐标、放样点至置仪点方位角及距离；② 反算桩号及距中距离 )程序名 :1ZD-XYLb1 0:Norm 2F=1 ： ( 正反算判别， F=1 正算， F=2 反算 , 也可以改 F 前加？，改 F 为变量 )Z[1]=90 （与路线右边夹角）Prog ＂ THB ＂： F=1=>Goto 1:F=2=>Goto 2Lb1 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２．主程序２：高程序横坡程序 ( 设计任意点高程及横坡 )程序名： 2GCLbI 0:Norm 2“KM”?Z:?D:Prog”H”:Fix 3:” H=”:Locate 6,4,H◢“ I=”: Locate 6,4,I◢Goto 03. 主程序 3 ：极坐放样计算程序 ( 计算放样点至置仪点方位角及距离 )程序名： 3JSX ： Y ：1268 ．123→K( 置仪点 X 坐标 )2243 ．545→L （置仪点 Y 坐标，都是手工输入 , 也可以建导线点数据库子程序 , 个人认为太麻烦）Y-L→E ： X-K→F ： Pol(F,E):IF J<0:ThenJ+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-XYLbI 0:Norm 290→Z[1]( 涵洞中心桩与右边夹角，手工输入，也可以修改成前面加？后变为变量 )1→F:Prog”THB”:?L:Z[2]-Z[1] →E:X+Lcos(E) →X:Y+Lsin(E) →Y:Fix 3: ＂Ｘ = ＂： Locate 6,4, Ｘ◢＂Ｙ＝＂： Locate 6,4, Ｙ◢Ｐ rog ＂3JS”:Goto 0:5 ．主程序 5: 路基开挖边线及填方坡脚线放样程序（输入大概桩号及测量坐标、地面标高计算出偏移距离、桩号、距中距离、填挖高度）程序名： 5FBXLbI 0:Norm 2: 18→DimZ:2→F:90→Z[1]:Prog“THB”:Z:D:”M0”?M:M→Z[4]:D→Z[3]:Prog”6GD”:L→Z[6]:If D<0:Then 0.75-L→D:Goto H:Else L-0.75→D:Goto H:IfEndLbI H:Prog”H”:H-0.03-Z[4] →Z[5]:Z[6] →L:If Z[5]<0:Then –z[5]→G:Goto W:Else Z[5] →G:Goto T:LbI W:Prog “W0”:Z[10]+Z[11] →A: If G>A:Then Goto 1:Else IfG>Z[10]:Then Goto 2:Else Goto 3:IfEnd:LbI 1:L+Z[12]+Z[13]+Z[14]+(G-A)×Z[9]+Z[11]×Z[8]+Z[10]×Z[7]:Goto Z: LbI 2:L+Z[12]+Z[13]+(G-Z[10])×Z[8]+Z[10]×Z[7]:G oto Z:LbI 3:L+Z[12]+G×Z[7]:Goto z:LbI T:L+0.5→N:If G>Z[17]:Then(N+Z[18]+(G-Z[17])×Z[16]+Z[17]×Z[15])→S:Goto Z:Else(N+G×Z[15])→S:Goto z:LbI Z:Z[3]→D:Fix 2:Abs(D)-S→T:” L0=”:L Locate 6,4,T◢＂ＫＭ＝＂： Locate 6,4, Ｚ◢＂Ｄ＝＂： Locate 6,4, Ｄ◢“ TW=”: Locate 6,4,Z[5]◢Goto 06 ．主程序 6 ：路基标准半幅宽度计算程序 ( 对于设计有加宽渐变的有用，如路基宽度无变化，则把此程序直接输入半幅宽度值至 L)程序名： 6GDProg “G0”Z-C→E:(B-A)×E/S+A→L:L:7 ．坐标计算次程序（ＴＨＢ）程序名： THB18→DimZ ：＂ＫＭ＂ ?Z ： Prog ＂Ｘ０＂1÷P→ Ｃ： (P-R)÷(2HPR) → Ｓ：180÷π→ Ｅ：Ｆ =1=>Goto 1 ：Ｆ=2=>Goto 2←┘Lbl 1 ： ? Ｄ： Abs( Ｚ－O) →W ： Prog " Ａ " ： X ：Ｙ：Ｇ oto 3 LbI 2: Ｘ：Ｙ：Ｘ→I ：Ｙ→J ： Prog ＂Ｂ＂：Ｏ＋Ｗ→Z ：Ｄ→ Ｄ：Ｇ oto 3LbI 3:IF Ｆ＝１ Then Ｘ：Ｙ： Else Ｚ：Ｄ8. 正算子程序 ( Ａ )程序名： A0.25→A ：0.52→B ：0.44→N 0.0→K ：0.49→L ：0.5→M ：U+W(Acos(G+QEKW(C+KW Ｓ ))+Bcos(G+QELW(C+LW Ｓ ))+Ncos(G+QEMW(C+MWＳ ))+Bcos(G+QE(1-L)W(C+(1-L)WS))+Acos(G+QE(1-K)W(C+(1-K)WS))) →X ：V+W(Asin(G+QEKW(C+KW Ｓ ))+Bsin(G+QELW(C+LWＳ ))+Nsin(G+QEMW(C+MWS))+Bsin(G+QE(1-L)W(C+(1-L)WS))+Asin(G+QE(1-K)W (C+(1-K)W Ｓ))) →Y ：G+QEW(C+W Ｓ )+ Ｚ［１］→ Ｚ［２］： X+ Ｄ cos （Ｚ［２］）→ Ｘ： Y+ Ｄ sin （Ｚ［２］）→ Ｙ9. 反算子程序 ( Ｂ )程序名： BG-90→ Ｔ： Abs((Y-V)cos(T)-(X-U)sin(T)) → Ｗ：0→ Ｄ： Lbl 0 ： Prog " Ａ " ： T+QEW(C+W Ｓ) → Ｌ： (J-Y)cos(L)-(I-X)sin(L)→ Ｄ： IF Abs( Ｄ )<0.01:Then Goto1 ： Else W+ Ｄ→W:Goto 0←┘Lbl 1:0→D ： Prog " Ａ " ： (J-Y)÷sin( Ｚ［２］) →D:10 ．高程计算子程序（ H ）程序名： HProg “S0”:R:T:C:G:I:C-T→F:Z-F→L:C+T→E:G-TI→Q:If T=O:ThenQ+LI→H:Goto 0:Else If Z<F:Then Q+LI→H:Goto 0:Else If Z≤E:ThenQ+LI+L2÷2÷R→H:Goto 0:LbI 0:H:If D=0:Then Goto I:Else Prog “I”:H+V→H:Goto I:LbI I:H:I: 11. 高程超高计算程序（ I ）程序名： IProg”I0”:W=1=> Goto 0:W=2=>Goto 1:LbI 0:If L=0:Then Abs(D)×M→V:Goto 2:ElseAbs(D)×((N-M)×(Z-C)÷L+M)→V:Goto 2:IfEnd:LbI 1:If L=0:Then Abs(D)×M→V:Goto 2:ElseAbs(D)×(((3((Z-C)÷L)2-2((Z-C)÷L)∧(3))×(N-M))+M)→V:Goto 2:IfEnd: LbI 2:Abs(D)→E:V÷E→I:I(E-K)→V:12 ．数据子程序 ( 附后示例 )① 程序名： X0 （坐标计算要素程序）If Z≥25900 And Z≤26615.555:Then25900→O:11587.421→U:1847.983→V:101 。

卡西欧5800任意线元坐标正反算程序（包括直线，圆曲线，缓和曲线正算及反算）“A=?,NEW=0，OLD≠0”?A: A≠0=>Goto S:“O=”?O:“U=”?U:“V=”?V:“G=”?G:“H=”?H:“P=”?P:“R=”?R:“Q=”?Q: Lbl S ：Fix 3:7→DimZ:1÷P→C:(P-R)÷(2HPR)→D:180÷∏→E:“1=>ZS,2=>FS”?N:If N=1：Then“Z-H=”?S: Z[2] →Z:“-L,+R=”?Z: Z→Z[2]：S-O→W:If W<0 or W>H:Then“=>ERROR!”◢ Goto S:IfEnd:Else “X=”?X:“Y=”?Y:X →I:Y→J:G-90→T:Abs((Y-V)cos(T)-(X-U)sin(T)→W:0→Z:IfEnd:Lbl A: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)))→Z[5]:V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+QEFW(C+FWD))+Asin(G+QEMW(C +MWD)))→Z[6]:G+QEW(C+WD)+90→F:Z[5]+Zcos(F)→Z[5]:Z[6]+Zsin(F)→Z[6]:T+QEW(C+WD)→L:(J-Z[6])cos(L)-(I-Z[5])sin(L)→Z:If N=1:Then“FWJ=”:F-90◣DMS ◢“X=”:Z[5]→X◢“Y=”:Z[6]→Y◢Goto S:IfEnd:If Abs(Z)＜0.00000001:Then(J-Z[6])÷sin(F)→Z:O+W→S:Else W+Z→W:Goto A:IfEnd:If S<O or S>O+H:Then“=>ERROR!”◢Goto S:IfEnd:“Z-H=”:S◢“-L,+R=”:Z◢Goto S说明：A：等于0为输入新的线元参数数据，反之进入坐标正反算。

卡西欧5800公路坐标正反算程序(总19页)-CAL-FENGHAI.-(YICAI)-Company One1-CAL-本页仅作为文档封面，使用请直接删除目录一、坐标正算基本公式 (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αβBd 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 成反比。

故此任意点的曲率为;CL R i i ==1ρ (c 为常数) （1） 由式(1)可知，回旋曲线任意点的曲率按线性变化，由此回旋曲线上里程为L i 点的曲率为;AB A i A B A i L L L L --⨯-+=)(ρρρρ （2） 当曲线右偏时ρB 、ρA 取正值，反之取负值。

曲线任意里程中边桩坐标正反算(CASIO fx-5800P计算器)程序一、程序功能及原理1.功能说明：本程序由一个主程序(TYQXJS)和五个子程——正算子程序(SUB-ZS)、反算子程序(SUB-FS)等构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。

本修改版程序既可实现正算全线贯通，亦可实现反算全线贯通。

本程序在CASIO fx-5800P计算器运行。

2．计算原理：利用Gauss-Legendre 5点通用公式正算线路中边桩坐标、线外测点至曲线元起点和终点的垂距的符号是否相异（即Dca×Dcb<0=>该测点在其线元内）进行判断并利用该线元要素反算中桩里程、支距，最后计算出放样数据。

二、源程序1.主程序(TYQXJS)Deg:fix 3119→DimZ“INPUT(0) Or DATA(Else)”?ILbl 0:“1.SZ=>XY,2.XY=>SZ，3.TF=>CK,4.SD=>FY，5.TW=>FY”?NIf N=1 Or N=5:Then Goto 1Else If N=2 Or N=3 Or N=4:Then Goto 2Else Goto 3IfEnd:IfEndLbl 1:“K(m)=”?SIf S<0:Then Goto 0:IfEnd“JL(m)=”?ZIf Z≠0:Then “ANGLE→R(Deg)=”?M:IfEndIf I=0:Then Prog “DAT1”:Else Prog “DAT2”:IfEndS-O→W:If W<0:Then Goto 0:Else If W>H:Then Goto 0:IfEnd:IfEndProg “SUB-ZS”:Prog “SUB-GC”If Z<0:Then“XL(m)=”:X◢“YL(m)=”:Y◢If N=5:Then Prog “SUB-TW”:IfEndElse If Z>0:Then “XR(m)=”:X◢“YR(m)=”:Y◢If N=5:Then Prog “SUB-TW”:IfEndElse “X(m)=”:X◢“Y(m)=”:Y◢“Hs(m)=”:L◢“FWJ=”: F►DMS◢IfEnd:IfEndGoto 1Lbl 2:“X(m)=”?X:If X<0:Then Goto 0:IfEnd“Y(m)=”?YIf N=3 Or N=4:Then “H(m)=”?→Z[2]:IfEndX→Z[4]:Y→Z[5]:90→MIf I=0:Then Prog “DAT1”G+QEH(C+HD)→F:1→J:U→Z[8J]:V→Z[8J+1]G-M→A:F-M→B:Prog “SUB-ZX1”:Goto BG+M→A:F+M→B:Prog “SUB-ZX1”:Goto B:IfEnd0→J:Lbl A: Isz JZ[8J+3]-M→A:Z[8(J+1)+3]-M→B:Prog “SUB-ZX1”If Z[6]Z[7]<0:Then J→J:Prog“DAT1” :Goto B:IfEndZ[8J+3]+M→A:Z[8(J+1)+3]+M→B:Prog “SUB-ZX1”If Z[6]Z[7]<0:Then J→J:Prog “DAT1” :Goto B:Else Goto A:IfEndLbl B:Prog“SUB-FS”O+W→S:Prog “SUB-GC”“K(m)=”:S◢“Hs(m)=”:L◢“JL(m)=”:Z◢If N=3:Then Prog “WIDE”:IfEndIf N=4:Then Prog “SDFY”:IfEndGoto 2Lbl 3:“TYQXJS→END”2. 正算子程序(SUB-ZS)0.1184634425→A: 0.2393143352→B:0.2844444444→Z[1]:0.0469100770→K:0.2307653449→L:0.5→Z[3]U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Z[1]cos(G+QEZ[3]W(C+Z[3]WD))+Bcos(G+QE(1 -L)W(C+(1-L)WD))+Acos(G+QE(1-K)W(C+(1-K)WD)))→XV+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Z[1]sin(G+QEZ[3]W(C+Z[3]WD))+Bsin(G+QE(1-L)W(C+(1-L)WD))+Asin(G+QE(1-K)W(C+(1-K)WD)))→YG+QEW(C+WD)+M→FX+Zcos(F)→X:Y+Zsin(F)→YF-M→F:If F<0:Then F+360→F:IfEndReturn3. 反算子程序(SUB-FS)G-M→T:Abs((Y-V)cos(T)-(X-U)sin(T))→W:0→ZLbl C:Prog “SUB-ZS”T+QEW(C+WD)→L:(Z[5]-Y)cos(L)-(Z[4]-X)sin(L)→ZIf Abs（Z）<10-6:Then Goto D:Else W+Z→W：Goto C:IfEndLbl D:0→Z:Prog “SUB-ZS”:(Z[5]-Y)÷cos(F)→Z:Return4．高程计算子程序（SUB-GC）Prog “DAT3”0．5RAbs(0.01(E-D))→TIf E>D:Then 1→C:Else -1→C:IfEndIf S<A:Then -0.01D→W:Else 0.01E→W:IfEndAbs(A-S)→LIf L>T:Then B+WL→L:Else B+WL+C(T-L)2÷2÷R→L:IfEndReturn5. 垂距计算子程序(SUB-ZX1)(Z[5]-Z[8J+1])cos(A)-(Z[4]-Z[8J])sin(A)→Z[6](Z[5]-Z[8(J+1)+1])cos(B)-(Z[4]-Z[8(J+1)])sin(B)→Z[7]Return6.曲线元要素数据库：DAT1(两条线路，可扩充为多条线路)If I=0:Then Goto E: Else If I=2:Then Goto F:IfEnd:IfEndJ=1=>***→Z[8]:***→Z[9]:***→Z[10]:***→Z[11]:***→Z[12]:***→Z[13]:***→Z[14]:*** →Z[15]：IfEndJ=2=>***→Z[16]:***→Z[17]:***→Z[18]:***→Z[19]:***→Z[20]:***→Z[21]:***→Z[22]:***→Z[23]:IfEndJ=3=>***→Z[24]***→Z[25]:***→Z[26]:***→Z[27]:***→Z[28]:***→Z[29]:***→Z[30]:***→[31]:IfEnd………………………………………………..Goto HLbl F: J=1=>***→Z[8]:***→Z[9]:***→Z[10]:***→Z[11]:***→Z[12]:***→Z[13]:***→Z[14]:*** →Z[15]：IfEndJ=2=>***→Z[16]:***→Z[17]:***→Z[18]:***→Z[19]:***→Z[20]:***→Z[21]:***→Z[22]:***→Z[23]:IfEndJ=3=>***→Z[24]***→Z[25]:***→Z[26]:***→Z[27]:***→Z[28]:***→Z[29]:***→Z[30]:***→[31]:IfEnd(注：如有多个曲线元要素继续添加入数据库DAT1中)Goto GLbl E:“X0=”?U:“Y0=”?V:“S0=”?O:“F0=”?G:“LS=”?H:“R0=”?P:“RN=”?R:“Q=”?Q Lbl G:1÷P→C:(P-R)÷(2HPR)→D:180÷π→E:Return7．隧道放样子程序：SDFY5.5→A:8.5→B:L+1.5→CIf I=1:Then -1→K:Else 1→K:IfEnd 给左右隧道符号赋值1.25K→F:4.25K→L:7.25K→K 给各圆心距路线的偏距赋值Z[2]-C→D 计算测点至圆心的高差D，If D>A:Then Goto 1:Else If D<0:Then Goto 2:IfEnd:IfEndZ-L→E:Abs(E)→E:√(D2+E2)→W:W-A→T: 计算测点至圆心的平距E，实测半径W“R(m)=”:W◢“DR(m)=”:T◢If E≤A:Then D-√(A2-E2)→C:E-√(A2-D2)→L 计算测点水平偏移量L、垂直偏移量C “Dh(m)=”:C◢显示测点垂直偏移量C“DL(m)=”:L◢显示测点水平偏移量LElse E-√(A2-D2)→L“DL(m)=”:L◢IfEnd:ReturnLbl 1Z-L→E:Abs(E)→E:√(D2+E2)→W:W-A→T: 计算测点至圆心的平距E，实测半径W“R(m)=”:W◢“DR(m)=”:T◢If E≤A:Then D-√(A2-E2)→C“Dh(m)=”:C◢IfEnd:ReturnLbl 2:Abs(D)→D:Abs(Z)→Z:Abs(K)→K:Abs(F)→FIf Z>K:Then Z-F→E:Else Z+K→E:IfEnd√(D2+E2)→W:W-B→T“R(m)=”:W◢“Rc(m)=”:T◢If E≤B:Then D-√(B2-E2)→C:E-√(B2-D2)→L 计算测点水平偏移量L、垂直偏移量C “Dh(m)=”:C◢显示测点垂直偏移量C“DL(m)=”:L◢显示测点水平偏移量LElse E-√(B2-D2)→L“DL(m)=”:L◢显示测点水平偏移量LIfEnd:Return8．曲线元要素判断数据库：DAT2(两条线路，可扩充为多条线路)If I=2:Then Goto A:IfEndIf S<第一线元终点里程:Then 1→J:Prog “DAT1”:ReturnElse If S<第二线元终点里程:Then 2→J:Prog “DAT1”:ReturnElse If S<第三线元终点里程:Then 3→J:Prog “DAT1”:ReturnElse If S<第四线元终点里程:Then 4→J:Prog “DAT1”:ReturnIfEnd:IfEnd:IfEnd:IfEndIf S≥第四线元起点里程And S<第五线元终点里程：Then 5→J：Prog “DAT1”：Return .................................Else If S<第n线元里程：Then n→J：Prog “DAT1”：ReturnIfEnd:IfEnd:IfEnd:IfEndLbl AIf S<第一线元终点里程:Then 1→J:Prog “DAT1”:ReturnElse If S<第二线元终点里程:Then 2→J:Prog “DAT1”:ReturnElse If S<第三线元终点里程:Then 3→J:Prog “DAT1”:ReturnElse If S<第四线元终点里程:Then 4→J:Prog “DAT1”:ReturnIfEnd:IfEnd:IfEnd:IfEndIf S≥第四线元起点里程And S<第五线元终点里程：Then 5→J：Prog “DAT1”：Return Else If S<第n线元里程：Then n→J：Prog “DAT1”：ReturnIfEnd:IfEnd:IfEnd:IfEnd .................................9．高程数据库子程序：DAT3（两条线路，可扩充为多条线路）If I=0:Then Goto A:Else If I=2:Then Goto B:IfEnd:IfEndIf S≤46500:Then 46150→A:361.26→B:40000→R:-0.3→D:-1.2→E:ReturnElse If S≤47200:Then 46800→A:353.46→B:20000→R:-1.2→D:1→E:ReturnElse If S≤48400:Then 47500→A:360.46→B:20000→R:1→D:-1.6→E:ReturnElse If S≤49600:Then 48700→A:341.26→B:10000→R:-1.6→D:2.45→E:Return IfEnd:IfEnd:IfEnd:IfEndIf S>49600 And S≤50300:Then 49840→A:369.19→B:100000→R:2.45→D:2.8→E:Return Else If I=2:Then49800→A:368.21→B:150000→R:2.45→D:2.7→E:ReturnIfEnd:IfEndLbl A:“C PEG(m)=”?→A“C H(m)=”?→B“R(m)=”?→R“I1=”?→D“I2=”?→EReturnLbl B49800→A:368.21→B:150000→K:2.45→D:2.7→E:Return10.填方路基宽度检查子程序（WIDE）“WIDTH CHECK”◢“H→MID(m)=”?→A 起坡点与中桩高差常数“DIST→MID(m)=”?→B 起坡点与中桩距离常数“i=”?→C 路基填筑坡度CAbs(Z[2]-L+A)→DAbs(Z)-B-D→E“DL(m)=”:E◢显示实有宽度与应有宽度之差值Return11.坡口坡脚线坐标计算子程序:SUB-TW“TIANWA LINE”◢“DIST→MID(m)=”?A 起坡点与中桩距离常数“H→MID(m)=”?B 起坡点与中桩高差常数“T(1) Or W(Else)=”?C 填方边坡为1，其它为挖方边坡“i=”?E 填方或挖方边坡坡度Abs(Z)-A→D: DE→KIf C=1:Then L-B-K→W:Else L+B+K→W:IfEnd“JS H(m)=”:W◢计算出应有高程“SC H(m)=”？→D 输入实测高程D-W→T“GC(m)=”:T◢显示实测高程与计算高程的差异Return三、使用说明1、规定(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。

CASIO fx-5800P线元法坐标正反算程序说明：本程序适用于卡西欧计算器 CASIO fx-5800P,可对全线贯通坐标正反算、竖曲线高程计算。

该程序可计算任意线型，包含（直线、圆曲线、缓和曲线、卵形曲线）等,还可以能通过坐标反推该点里程和距中线距离,适用测量员专用。

主程序名：ABCYT第1步Deg：Fix 3：10→DimZ第2步Lbl 3：＂1.DK=>XY＂：＂2.XY=>DK＂：＂Q＂?W：＂DK＂?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：＂W＂?Z：＂α＂?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：＂S=＂：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])→Z[10]第2步If Abs(Z[10])>0.001：Then S+Z[10]→S：Goto 2：Else Goto 1：IfEnd第3步Lbl 1：(Y-J)÷sin(Z[9])→Z第4步Pol(X-I，Y-J)：If Z<0：Then -1×I→Z：Else 1×I→Z：IfEnd数据库名：ABCYTSJ第1步Goto1第2步Lbl 1：If S<7586.707 Or S>13346.96：Then Cls：Locate 2,2,＂PQX＂：Locate 4,3,＂CHAOXIAN＂：Locate 10,4,＂→Stop＂◢第3步Stop：IfEnd第4步Lbl 1：IF S<7946.707：Then 98°56′56″→G：7586.707→o：3378605.445→U：453648.704→V：0→P：4500→R：360→H：1→Q：Return：IfEnd第5步Lbl 1：IF S<11766.03：Then 101°14′26″→G：7946.707→o：3378544.714→U：454003.518→V：4500→P：4500→R：3819.323→H：1→Q：Return：IfEnd第6步Lbl 1：IF S<12126.03：Then 149°52′11″→G：11766.03→o：3376389.890→U：457018.324→V：4500→P：0→R：360→H：1→Q：Return：IfEnd第7步Lbl 1：IF S<13346.96：Then 152°09′41.68″→G：12126.03→o：3376073.846→U：457190.654→V：0→P：0→R：1220.93→H：0→Q：Return：IfEnd第n步……………………………………………………数据输入说明：第1步Goto1第2步Lbl 1：If S<本条线路起点里程S>本条线路止点里程Cls：Locate 2,2,＂PQX＂：Locate 4,3,＂CHAOXIAN＂：Locate 10,4,＂→Stop＂◢第3步Stop：IfEnd第4步Lbl 1：If S<本线元止点里程：Then线元起点切线方位角→G：线元起点桩号→O：线元起点坐标X→U：线元起点坐标Y→V：线元起点半径（直线为0、曲线为半径）→P：线元止点半径（直线为0、曲线为半径）→R：线元长度→H：线元转向（左转为-1、右转为1、直线为0）→Q：Return：IfEnd第n步Lbl 1：每增加一行则为增加一个线元要素。

卡西欧5800程序完整版＂ZS=1,FS=2＂？N: 1÷P →C:(P-R)÷(2HPH)→D:180÷π→E: If N=1: Then Goto 1: Else Goto 2: IfE nd←┚LbI 1:1 →I:0 →Z：Abs(K-Ｏ) →W: Prog＂ZS＂←┚F-90→F←┚＂JJ＂?K:＂PJ＂?L: ＂ZX=＂：X+Lcos(F+K) →Z[3]◢＂ZY=＂: Y+Lsin(F+K) →Z[4]◢＂QXFWJ=＂:F DMS◢Abs(Z) →Z:Goto 8←┚LbI 2:＂X＂?I:＂Y＂?J: Prog＂FS＂:＂FSZH=＂:Ｏ+W→K◢＂FSJL=＂:Z◢Lbl 8正算子程序（ZS）（不运行）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+ZcosF→X：Y+ZsinF→Y3反算子程序（FS）（不运行）G-90→T：Abs（（Y-V））cos（T）-（X-U）sin（T）→W：0→Z←┚Lbl 0←┚prog＂ZS＂：T+QEW（C+WD）→L：（J-Y）cosL-（I-X）sinL→Z：If Abs(Z)＜0.001: Then Goto1：E lse W+Z→W：IfEnd :Goto 0←┚Lbl 1：0→Z：Prog＂ZS＂：（J-Y）÷sinF→ZZX-WN（运行主程序，该程序适用于匝道，对全段完整非对称曲线计算太过麻烦）＂CZX＂?S(仪站X坐标):＂CZY＂?W（仪站Y坐标）LbI 0: ＂K＂?K: If K≤126.9（第1条单一曲线的终点下面同上）: Then Goto A：Else If K≤176.9 And K＞126.9: Then Goto B:Else If K≤259.562 And K＞176.9: Then Goto C:Else If K≤309.562 And K ＞259.562: Then Goto D：IfEnd:IfEnd:IfEnd:IfEnd←┚LbI A:92094.336（起点X坐标）→U:8287.962（起点Y坐标）→V:0起点桩号→O:145°11＇54＂（方位角）→G:126.9（曲线长度）→H:10^(45) （起点半径）→P:10^(45) （止点半径）→R:+1（路线左右偏，左负右正，直线输0）→Q:Prog ＂WN＂:Goto WLbI B:91990.134→U:8360.388→V:126.9→O:145°11＇54＂→G:50→H:10^(45) →P:150 →R:+1→Q:Pr og ＂WN＂:Goto WLbI C:91947.609→U:8386.570→V:176.9→O:154°44＇51＂→G:82.662→H:150 →P:150 →R:+1→Q:Pr og ＂WN＂:Goto WLbI D:91867.101→U:8399.992→V:259.562→O:186°19＇19＂→G:50→H:150 →P: 10^(45) →R:+1→Q: Prog ＂WN＂:Goto WLbI W:＂PJ1＂?T←┚（该项是只左右边桩再左右偏距，输0为不再左右偏，可以直接出方位角、距离，适用于桥梁桥台、盖梁、涵洞八字墙等放样，以路线前进方向，前-，后+）If T=0：Then Z[3]→Z: Z[4]→M: Goto P: Else Goto V: IfEnd←┚LbI V:＂JJ1＂? O←┚（左右边桩再偏夹角，用180°减去路线夹角）＂X1=＂: Z[3]+Tcos（F+N+ O）→U◢＂Y1=＂:Z[4]+Tsin（F+N+ O）→V◢U→Z: V→M: GotoP←┚LbI P:Pol(Z-S,M-W): If J≤0：Then J+360°→J: IfEnd←┚＂JL=＂:I◢（仪站至测点的距离）＂FWJ=＂:J DMS◢（仪站至测点的方位角）Goto0备注：1．规定以道路中线前进方向（即里程增大的方向）区分左右；当线元往左偏时Q=-1；当元线往右偏时，Q=1：当元线为直线时，Q=0。