卡西欧fx-5800P直线圆曲线坐标正反算程序
CASIO5800坐标正反算计算程序
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 0K——计算点的里程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÷π→F4. 数据库（命名为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中，程序中所有公式和部分函数结果均存储在ListAns列表数组中，要想多次调用最好随公式取出结果，并赋给变量。
卡西欧fx-5800p计算器测量程序
卡西欧fx-5800p计算器测量程序（作者：水江华）1.坐标反算程序"ZBFS QH1-4"◢Deg：Fix 3←┘"X+YiS(m）="?A←┘Lb1 0:"X+YiE(m)，0=>END="?B←┘Abs（B）=0=>Goto E←┘Arg（B-A)→JJ<0=>J+360→J"DIST(m）="：Abs（B-A)◢"α（DMS)=":J▶DMS◢Goto 0←┘Lb1 E:"ZBFS QH1-4=>END"2.坐标正算程序"ZBZS QH1-3"◢Deg：Fix 3←┘"XS(m)= "？A:"YS(m)= "？B←┘Lb1 1:"DIST(m)= "？L←┘If L≤0:Then Goto 2:IfEnd←┘"α（Deg)="？R←┘A+Lcos（R）→C←┘B+Lsin（R）→D←┘Fix 3："XE(m)= ":C◢"YE(m)= ":D◢Goto 1←┘Lb1 2:"ZBZS QH1-3=>END"3.单一闭、附合导线近似平差程序"BFHDXPC QH1-8"◢Deg:Norm 1:FreqOn←┘n→N←┘If List Y[N]=0:Then N-1→DimZ:Goto 0←┘Else N→DimZ:Goto 1: IfEnd←┘Lb1 0: "CLOSE Or CONECT TRA VE"◢"ANGLE NUM=":N◢"SIDE NUM=":N-1◢"UNKNOWN P NUM=":N-2◢"CLOSE(0)，CONNECT(≠0)="?→Z←┘Fix 3："XA(m)，≤0=>αA→B(Deg)="？A←┘If A>0：Then "YA(m)= "？B：Else "αA→B(Deg)="？R：IfEnd←┘"XB(m)= "？C:"YB(m)= "？D←┘C+Di→U：U→V←┘If A>0：Then Pol（C-A，D-B）：Cls←┘J<0=>J+360→J：J→R←┘"DIST A→B(m)="：I◢"αA→B(DMS)="：R▶DMS◢IfEnd←┘If Z≠0：Then "XC(m)= "？E:"YC(m)= "？F←┘E+Fi→V←┘"XD(m)，≤0=>αC→D(Deg)="？G←┘If G>0：Then "YD(m)= "？H←┘Pol（G-E，H-F）：Cls←┘J<0=>J+360→J：J→S←┘"DIST C→D(m)="：I◢"αC→D(DMS)="：S▶DMS◢Else "αC→D(Deg)="？S：IfEnd←┘Else If R>180：Then R-180→S：Else R+180→S：IfEnd←┘IfEnd←┘∑y→M←┘For 1→I To N←┘If I=1：Then R+List X[I]→L：Else L+List X[I]→L：IfEnd←┘If L>180：Then L-180→L：Else L+180→L：IfEnd←┘Nent←┘3600(L-S)→T←┘"α CLOSE ERROR(S)="：T◢40→W←┘Abs(T)>W=>"α CLOSE ERROR OVRE!"←┘-T÷N÷3600→P：0→O←┘For 1→I To N←┘If I=1：Then R+List X[I]+P→L←┘Else L+List X[I] +P→L：IfEnd←┘If L>180：Then L-180→L：Else L+180→L：IfEnd←┘If I<N：Then List Y[I]<L→Z[I] ：Z[I]+O→O：IfEnd←┘Nent←┘3600(L-S)→T←┘"CHECK α CLOSE ERROR(S)="：T◢U+O-V→O←┘Int(M÷Abs(O))→K←┘"DELTA X(m)= "：ReP(O)◢"DELTA Y(m)= "：ImP(O)◢"RELAT CLOSE ERROR=1÷"：K◢K<4000=>"RELAT CLOSE ERROR OVER! "←┘For 1→I To N-1←┘Z[I]-OList Y[I]÷M→X←┘If I=1：Then U+X→Y：Else Y+X→Y：IfEnd←┘Norm 1："POINT n="：I◢Fix 3："X(m)= "：ReP(Y)◢"Y(m)= "：ImP(Y)◢Nent←┘Y-V→Q←┘"CHECK X (m)= "：ReP(Q)◢"CHECK Y (m)= "：ImP(Q)◢Goto E←┘Lb1 1："NO DIRECT TRA VE"◢"ANGLE NUM="：N-1◢"SIDE NUM="：N◢"UNKNOWN P NUM="：N-1◢"XA(m)= "？A:"YA(m)= "？B←┘"XB(m)= "？C:"YB(m)= "？D←┘A+Bi→U:C+Di→V←┘0→H:U+List Y[1]<H→Z[1]←┘For 1→I To N-1←┘H+List X[1]→H←┘If H>180：Then H-180→H：Else H+180→L：IfEnd←┘以下未完程序由读者完成。
卡西欧fx5800p计算器 坐标正反算程序
(以下程序是专业人士编写，本店铺不对程序负责，仅供您参考使用。
)卡西欧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.1184634425→A ：0.2393143352→B ：0.2844444444→N 0.046910077→K ：0.2307653449→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←┘L bl 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 。
CASIO5800线元法程序
卡西欧5800P 线元法正反算程序1-JIN YU(主程序)Lbl 4: “1→ZS,2→FS?N(选择计算模式，1为正算，2为反算，3为超欠挖)N=1=>Goto 1:N=2=>Goto 2Lbl 1: “K=”?S：“P=”?Z：Prog “SJ-PM”： Abs(S-O) → W：Prog “SUB1-ZS”：“X=”：Locate4,4,X:“Y=”：Locate4,4,Y:F-90→F: “W=”: Locate4,4,F°◢Goto 4(正算-输入待求点里程K=、输入待求点偏距P=、显示待求点里程X=、显示待求点里程Y=、显示待求点里程设计高程H=)Lbl 2:“X=”?X：“Y=”？Y:Prog“SJ-PM”： X→ I： Y→J：Prog "SUB2-FS"：O+W→S: “K="：Locate4,4, S：“P="：Locate4,4, Z：S→ K:Prog“SJ-GC”：“H=”：Locate4,4, H:F-90→F: “W=”: Locate4,4,F°◢Goto 4(反算-输入实测点X=、输入实测点Y=、显示实测点里程K=、显示实测点偏距P=、显示实测点的设计高程H=)SJ-PM（子程序名-平面线形数据库）If S ≥59227.681（线元起点里程）:Then 2599818.013→U（线元起点X坐标）：496887.918→V（线元起点Y坐标）：59227.681→O（线元起点里程）：208028’39”→G（线元起点方位角）：100 →H（线元长度）：1×1045→P（线元起点曲率半径）：1×1045→R（线元终点曲率半径）：0→Q（线元左右偏标志：左-1右1直0）：IfEndIf S ≥59327.681（线元起点里程）:Then 2599730.112→U（线元起点X坐标）：496840.237→V（线元起点Y坐标）：59327.681→O（线元起点里程）：208028’39”→G（线元起点方位角）：90 →H（线元长度）：1×1045→P（线元起点曲率半径）：750→R（线元终点曲率半径）：-1 →Q（线元左右偏标志：左-1右1直0）：IfEnd(注：有多个平面线形，依照上面的依次变更, 每多一个，就增加一个。
卡西欧fx-5800 坐标正算程序
卡西欧fx-5800 坐标正算程序“SO”算方位Lbl 0“D”：Pol(X-A，Y-B) →D◢If J≥0:Then J◢Else “J”:J+360◢IfEnd“ZZ”正算直线坐标A+Rec(Abs(P-K),Q) →XB+J→Y“X”：X+Rec(D,Q+N×90) →X◢“Y”:Y+J→Y◢P?所求桩号 A，B：起点x,y K:起点桩号 Q：方位角D：边桩距离 N？左偏“-1”右偏：“+1”X= Y= 所求点坐标“YY”正算圆曲线坐标Q+C×（P-K）÷S×180÷∏→Z“X”：A+C×S×（Sin(Z)-Sin(Q)）+Rec(D,Z+N×90) →X◢“Y”:B-C×S×（Cos(Z)-Cos(Q)）+J→Y◢P:所求点桩号 A，B：起点x,y K：起点桩号Q：方位角 S：曲线半径 D？边桩距离N？左偏“-1”，右偏“+1”X= y= 所算点坐标 C：线路左转“-1”右转“+1”“H1”正算缓和曲线坐标Abs(P-K) →LT2÷L→HL2×180÷∏÷2÷T2→VQ+C×V→ZL-L∧（3）÷40÷H2+L∧（5）÷3456÷H∧（4）-L∧（7）÷H∧（6）÷599040→EL2÷6÷H- L∧（4）÷336÷H∧（3）+ L∧（6）÷42240÷H∧（5）- L∧（8）÷H∧（7）÷（6×1612800）→FIf C<0：Then-F→F：IfEndA+E×Cos(Q)-F×Sin(Q) →X“X”：X+Rec(D,Z+M×N×90) →X◢B+E×Sin(Q)+F×Cos(Q) →Y“Y”：Y+J◢P？所求点桩号 A，B：起点x,y K:起算点里程桩号Q：前切线方位角 T（是圆曲线要素的A1）=√（根号）缓曲总长×半径C=切线线路方向左转“-1”，右转“+1“M=第一缓曲“+1”，第二缓曲”-1“N：算边桩左偏“-1”，右偏“+1”D：边桩距离坐标正算数据库例子“肇兴隧道右线数据库”Lbl 0“P”?P:“N”？N：“D”？DIf P≥16836.5 And P≤17009.624:Then 66881.803→A：14635.782→B：16836.5→K:280032030.50→Q:902.044→T:1→C: 1→M:Prog“H1”：IfEnd (第一缓曲线)If P≥17009.624 And P≤17529.681:Then 66914.52→A：14465.781→B：17009.624→K:281035049.370→Q:4700→S：Prog“YY”:IfEnd （圆曲线）If P≥17539.681 And 17702.805≥P:Then 67102.343→A：13799.117→B：17702.805→K:108059031.50→Q:902.044→T:-1→C: -1→M:Prog“H1”：IfEnd (第二缓曲线)If P≥17702.805 And P≤19550:Then 67102.343→A：13799.17→B：17702.805→K:288059031.50→Q: Prog“ZZ”：IfEnd （直线）。
CASIO fx-5800计算器测量坐标程序
主程序坐标正反算程序名称：001XYLB1 A↙Cls : Fix 4 : 30→Dim Z :"v=1,v=0"?Z ↙If Z=1:Then Goto B:IfEnd↙(Z=1进入里程点坐标正算）If Z=0:Then Prog"FSLCBZ ": Goto A:IfEnd↙(Z=0进入反算里程边桩）Lb1 B ↙“K=”?K ：（计算里程）Prog"SJK1 ":Prog”ZBFY” :Goto A↙子程序数据库 SJK1IF 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 : Prog”XLZBJSCX”:Return: IFEND↙IF…………Prog”XLZBJSCX”:Return: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] :E →T:Return:IFEND子程序坐标正算，名称：XLZBJSCXLB1 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 : ELSE E-T →T : T<0=>360+T→T : IFEND ↙GOTO 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 : ELSE 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*S IN 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 ↙LLb1 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↙子程序反算里程边桩，名称：FSLCBZLb1 1↙"POINT=X"?X:" POINT=Y"?Y↙（输入任意测点的XY坐标）:“K=”?K （试算里程，任意输入标段内里程点即可，也可以输入大致的估算里程加快速度）↙Lb1 2↙Prog"SJK1 ":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"SJK1":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]:"W(-Z,Y+)=":I◢(偏距）（由于该程序不能准确判断边桩左右方向，暂作修改，取消左右边判定）"K=":K+S→K◢（里程）Return↙子程序坐标放样：ZBFYLB1 0 ↙"P(N)="?M:"P(E)="?N ↙“a=”?P◢(输入边桩与线路夹角，左-右+)“W=”?D◢（输入边桩距）Z[18]+D*COS(T+P) →Z[20] ↙Z[19]+D*SIN(T+P) →Z[21] ↙“XP=”: Z[20] ◢（边桩X）“YP=”: Z[21] ◢（边桩Y）Pol(Z[20]-M,Z[21]-N):"HD=":I◢（边桩放样距）J<0=>J+360→J:"HAR=":J→DMS◢（边桩放样方位角）Return↙直线短链要将该直线分成二条直线来对待编制数据库文件即可，曲线长链要单独建判定语句子程序并然后编制二个数据库（按来向里程和去向里程分别推出曲线起终点里程，以此编制曲线段的数据库）。
线路测量坐标正、反算原理及卡西欧5800计算器程序说明_secret..
线路测量坐标正、反算计算原理 及卡西欧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基本曲线要素及坐标计算各类测量书均对曲线要素及坐标计算有详细推到过程，此处仅列出计算公式，便于程序编制说明。
卡西欧5800计算器坐标正反算程序
卡西欧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↙。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
[精] fx 5800 直线圆曲线坐标正反算程序
(2010-07-01 21:50:11)
标签：
杂谈
直线计算程序
0→I:0→J:”X 0”?D:”Y0”?E:”X1”?B:”Y1”?C:Pol(B-D,C-E):J→A:If A<0: Then A+360→A:Else A→A:IfEnd:Lbl 0:?O:?S:If O≠0:Then Goto 1:IfEnd: D+Scos(A) →X:”X=”:X◢
E+Ssin(A) →Y:”Y=”:Y◢
Goto 0:Lbl 1:D+Scos(A)+Ocos(A+90) →X:”X=”:X◢
E+Ssin(A)+Osin(A+90) →Y:”Y=”:Y◢
Goto 0
输入程序时注意区别字母O与数字0
程序运行时符号说明
X0? Y0?分别输入直线起点的XY坐标值
X1? Y1?分别输入直线终点的XY坐标值
O? 输入边桩与中桩的距离(左边桩为负值,右边桩为正值),如计算中桩坐标输入0
S? 输入所求点到直线起点的距离
圆曲线计算程序
0→I:0→J: ”X0”?C:”Y0”?D:”X1”?E:”Y1”?F:?R:”L:-1 R:1”?N:”ZY”?W:Lbl 0:”LN”?T:
T-W→O:Pol(E-C,F-D):I→S:J→A: If A<0:Then A+360→A:Else A→A:IfEnd: sin-1(S÷(2R)) →K:2∏RK÷180→L:180O÷(2∏R) →G:(2R)sin(G)
→H:C+Hcos(A-KN+GN) →X:”X=”:X◢
D+Hsin(A-KN+GN) →Y:”Y=”:Y◢
R-0.5√(4R2-H2) →Q:”Q=”:Q◢
“S=”?V:If V=0:Then Goto 0:IfEnd:”L:-90 R:90”?U:X+Vcos(A-KN+2GN+U) →X:”X=”:X◢
Y+Vsin(A-KN+2GN+U) →Y:”Y=”:Y◢
Goto 0
输入程序时注意区别字母O与数字0
程序运行时符号说明
X0? Y0?分别输入直线起点的XY坐标值
X1? Y1?分别输入直线终点的XY坐标值
R? 输入圆曲线半径
L:-1 R:1?圆曲线向左转弯时输入-1,向右转弯时输入1
ZY? 输入起点桩号
LN? 输入所求点桩号(按下exe键后显示中桩坐标)
S? 输入所要求的边桩与中桩的距离
L:-90 R:90要求左边边桩输入-90,右边边桩输入90(按下exe键后显示边桩坐标)。