测量坐标计算程序K

隧道测量---快速坐标正反算程序(4800-4850计算器）2008-11-23 20:47:00点击：212ZB （坐标正算）1. Lb1 12. ｛L｝3.SMNARCL4. K=L-S5. R=0=>X“X”=M+KcosA▲6. Y“Y”=N+KsinA▲7. Z“FW”=A▲8.Goto 4△R≠0=>Abs K≤C=>Goto 2△R≠0=>K＞C=>Goto 3△R≠0=>K<0=>Abs K＞C=>Goto 1△9. Lb1 210. O=90K²÷（πRC）11. G=K³÷（6RC）-Kˆ7÷（336（RC）³＋Kˆ11÷（42240（RC）ˆ5-Kˆ15÷9676800（RC）ˆ7）＋Kˆ19÷（3530096640（RC）ˆ9）12. F=K－Kˆ5÷（（40（RC）²）－Kˆ9÷（3456（RC）ˆ4）＋Kˆ13÷（599040（RC）ˆ6）－Kˆ17÷（40320×4532（RC）ˆ8）13. J=0=>X“X”=M＋GsinA＋cosA▲14. Y“Y”=N＋FsinA- cosA▲15. Z“FW”=A-O▲Goto 4△16．J≠0=>X“X”=M＋F cosA－GsinA▲17．Y“Y”=N＋G cosA＋F sinA▲18．Z“FW”=A＋O▲Goto 4△19．LbI 320．B=90（2K-C）÷（лR）21.U=C÷2－C³÷（240R²）＋R sinB＋Cˆ5÷（34560Rˆ4）-Cˆ7÷（599040 Rˆ6）22．V=C²÷（24R）＋R（1－cosB）－Cˆ4÷（2688R³）＋Cˆ6÷（42240 Rˆ5）23．J=0=>X“X”=M＋UcosA＋VsinA▲Y“Y”=N＋U sinA－V cosA▲Z“FW”=A-B▲Goto 4△24．J≠0=>X“X”=M＋U cosA－VsinA▲Y“Y”=N＋U sinA＋V cosA▲Z“FW”=A＋B▲Goto 4△25．Lb1 426．｛D｝：｛T｝27．D=0 => Goto 1△28．D≠0=>X“LX”=X＋Dcos（T＋Z）▲Y“LY”=Y＋Dsin（T＋Z）▲Goto 1△SUB（坐反算）“1.SZ=>UV”: “2.UV =>SZ”: N（=2): X“X0”: Y“Y0”: O“S0”=L: G“F0”(=Z): Q：Q=0=>R=1E-45：≠> R△J=0：M=AbsQ：K=1-M：N=1=>Goto 1：≠>Goto 2△←┙Lbl 1：{SZ}：SZ：L=S-O：B=90QL÷π÷R：F=G+B：E=F+B+90：C=LK+2MRsinQB：U"US"= X+Rec(C,F)+ZcosE◢V "VS" =Y+J+ZsinE◢Goto 1←┙Lbl 2：W=G+90Q：E=X+Rec(R,W)：F=Y+J：W=W+180M：W>360=>W=W-360△←┙Lbl 3：{UV}：UV：D=Pol(U-E,V-F):T=J：J<0=>T=J+360△B=T-W：P=KB+MAbs B：S"S"=O+ DKcosP+ MπRP÷180◢H"H"=DKsinB+MQ(R-D)◢Goto 3←┙正算由下面输入：S： ZH或HZ点里程M： ZH或HZ点X坐标N： ZH或HZ点Y坐标A：切线方位角R：曲线半径C：缓和曲线长度L：求点里程X：求点X坐标Y：求点Y坐标FW：求点切线方位角D：左、右边桩距离T：边桩与线路夹角LX：边桩X坐标LY：边桩Y坐标J：右偏输“1”，左偏输“0”反算由下面输入：X0: 线路中心X坐标Y0: 线路中心Y坐标L: 相对应里程Z: 相对应方位角Q: ＝０U: 求点X坐标V: 求点Y坐标S= 求点里程H= 求点边距说明：１.可用于缓和曲线坐标正反算，反算时X0,Y0,L,Z利用正算储存的数据，速度快．２.反算时先输入近似里程正算．３.隧道断面炮位放样２０多个点位在操作熟练时只需15分钟左右，提高工程进度．。

高速公路测量CASIO4800&4850万能坐标计算程序(完整版）程序特点：真正的全线贯通坐标正反计算、任意斜角计算！！！程序中加入测站点，真正的实现了“坐标法”与“极坐标法”两种放样方法的同时显示的功能，使得放样操作方法选择时更加灵活！！！在曲线元要素输入时仅需要输入第一段全部曲线元要素，后面曲线元要素除起点半径、终点半径、曲线长、转向需输入外其他要素均从前一曲线按辛普森8等分计算得出，解决了主线坐标计算无法获得第二段及其以后曲线元起点参数的问题；辛普森公式任意等分，满足所有精度要求；全线曲线元数据一次性程序化输入，参数存储采用扩充变量数据库，无需修改程序内容；多功能采用单程序编程，避免频繁调用子程序，提高运算速度。

一、程序：ZBJSW“1.ZS 2.FS 3.SZ”:W=1=>Z[2]=0:V=0：Goto 1 ΔW=2=> Goto 4ΔW=3=> O “KOU LING”:O≠123456=>O=0: “OUT”◢Goto CΔO=0: V=0:Z[1]=0:Goto 0←┘Lbi 0←┘”N0.”：Z[1]+1 ◢Z[1]=0=>{ABCREFGUKO}:A“X0”:B“Y0”:C“F0”:R“R0”:E“RN”:F“D0”:G “LS”:U“G”:K“X（00）”: O“Y（00）”: Z[Z[1]×8+3]=A:Z[Z[1]×8+4]=B:Z[Z[1]×8+5]=C:Z[Z[1]×8+6]= R-1:Z[Z[1]×8+7]= E-1:Z[Z[1]×8+8]=F: Z[Z[1]×8+9]=F+G: Z[Z[1]×8+10]=U:“NEXT”◢Isz Z[1]: Goto 0ΔZ[1]=1=>D=Z[9]:Z=0:Z[2]=0:GOTO 2ΔD=Z[(Z[1]-1)×8+9]:Z=0:Z[2]=Z[1]-1:GOTO 2←┘Lbi A←┘Z[Z[1]×8+3]=X:Z[Z[1]×8+4]=Y:Z[Z[1]×8+5]=J: Z[Z[1]×8+8]=D: {REGU}:R“R0”：E “RN”: G“LS”:U“G”: Z[Z[1]×8+6]=R-1：Z[Z[1]×8+7]=E-1: Z[Z[1]×8+9]=D+G: Z[Z[1]×8+10]=U:“NEXT”◢Isz Z[1]: Goto 0←┘Lbi 1←┘{DZT }:D：Z：T“RJ”：Z[2]=0:Goto 2←┘Lbi 2←┘V≠1=>Z[2]>Z[1] =>GoToCΔΔD≤Z[Z[2]×8+9]=> A=Z[Z[2]×8+3]:B=Z[Z[2]×8+4]: C =Z[Z[2]×8+5]:R=Z[Z[2]×8+6]: E=Z[Z[2]×8+7]: F=Z[Z[2]×8+8]: G=Z[Z[2]×8+9]: U=Z[Z[2]×8+10]: Goto3ΔIsz Z[2]:Goto 2←┘Lbi 3←┘W=3 =>N=8：≠P＝U（E-R）÷Abs（G-F）：Q＝Abs（D-F）÷N：S=90Q÷π：J＝C+(NPQ+2UR)NS：L=1←┘X＝A+Q÷6×(Cos C+Cos J +4∑（Cos (C+((L+0.5)PQ+2UR)×(L+0.5)S),L,0,(N-1)）+2∑（Cos (C+((LPQ+2UR)LS,L,1,(N-1)))+ZCos（J+ T）←┘Y＝B+Q÷6×(Sin C+Sin J +4∑（Sin (C+((L+0.5)PQ+2UR)×(L+0.5)S),L,0,(N-1)）+2∑（Sin (C+((LPQ+2UR)LS,L,1,(N-1)))+Z Sin（J+T）:V=1=>Goto6ΔV=2=>Goto9ΔV=3=> GOTO CΔW=3=>GOTO AΔZ=0=>“X（Z）=”:X:Pause 0: “Y（Z）=”:Y◢Pol（（X－K），（Y－O））←┘“S（Z）=”:I ◢J<0=> J=J+360Δ“F（Z）=”: J→DMS◢Goto 1ΔZ<0=>“X（L）=”:X:Pause 0: “Y（L）=”:Y◢Pol（（X－K），（Y－O））←┘fx4850①“S（L）=”:I ◢J<0=> J=J+360Δ“F（L）=”: J→DMS◢Goto 1ΔZ>0=>“X（R）=”:X:Pause 0: “Y（R）=”:Y ◢Pol（（X－K），（Y－O））←┘“S（R）=”:J ◢J<0=> J=J+360Δ“F（R）=”: J→DMS◢Goto 1 ←┘Z=0=> X “X（Z）=”◢Y “Y（Z）=”◢Pol（（X－K），（Y－O））←┘I“S（Z）=”◢J<0=> J=J+360ΔJ“F（Z）=”◢Goto 1ΔZ<0=> X “X（L）=”◢Y “Y（L）”◢Pol（（X－K），（Y－O））←┘fx4800②I“S（L）=”◢J<0=> J=J+360ΔJ“F（L）=”◢Goto 1ΔZ>0=> X “X（R）=”◢Y “Y（R）=”◢Pol（（X－K），（Y－O））←┘I“S（R）=”◢J<0=> J=J+360ΔJ“F（R）=”◢Goto 1 ←┘Lbi 4←┘{MH} ：M“X”：H“Y”：Z[2]=0：GOTO 5←┘Lbi 5←┘V=1：D= Z[Z[2]×8+9]：Z=0：T=90：GOTO 2←┘Lbi 6←┘K=（（H -B）Cos（C-90）-（M-A）Sin（C-90））×（（H -Y）Cos（J-90）-（M-X）Sin（J-90））：K≤0=> Goto 7ΔIsz Z[2]：Goto5←┘Lbi 7←┘D=F+Abs((H -B)Cos(C-90)-(M-A)Sin（C-90）)：D>G=> Isz Z[2]: Goto5ΔGoto 8←┘Lbi 8←┘V=2 ：GOTO 3←┘Lbi 9 ←┘K=(H -Y)Cos(J-90)-(M-X)Sin(J-90)：Abs K<（1÷E）^3=>Goto BΔD=D+K ：GOTO 8←┘Lbi B←┘V=3 ：Z=0：Goto 3←┘Lbi C←┘Z=(H-Y) ÷Sin(J+90)：“D”:D:Pause 0: “Z”: Z◢4850输出(Z=(H-Y) ÷Sin(J+90)：D“D”◢Z “Z”◢4800输出)GOTO 4←┘Lbi C←┘二、说明a、编制说明本程序是运用复化辛普生公式根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距，对该曲线段范围内任意里程中边桩坐标进行计算，以及对卡西欧扩充变量的灵活应用，实现了真正意义上的的全线贯通及曲线要素输入程序化（在不修改程序内容的情况下可通过运行程序输入任意多段曲线元要素）。

4800公路测量放样计算程序CASIO fx-4800P计算器程序DA DI——N E (公路弯道坐标计算)K“JD”：R：S“LS”：F“PJ——L－R＋”：A“ZH——N”：B“ZH——E”：C“JD——N”：D“JD——E”M“ZH”=K－（R＋S2÷（24R））tan（Abs F÷2）－S÷2＋S^3÷（240R2）▲N“HY”=M＋S▲O“QZ”=M＋（πRAbs F÷180＋S）÷2▲Z[1]“YH”=M＋πRAbs F÷180▲Q“HZ”= Z[1]＋S▲T=K－M▲L=Q－M▲E=（R＋S2÷（24R））÷cos（F÷2）－R▲Lb1 0：{G}：G“ZHUANG HAO”：G＜O＝=＞Z=G－M≠=＞Z=Q－G：△{U}：U“I——B OUT＋IN－”Z＜S＋0.005＝=＞H=tan-1（24Z2S2R2－Z^6）÷（48Z^3R^3－6Z^4SR））：V=Z－Z^5÷（40S2R2）＋Usin H：W=Z^3÷（6SR）－Z ^7÷（336Z^3R^3）－Ucos H：≠=＞H=180（Z－S÷2）÷（πR）：V=Rsin H＋S÷2－S^3÷（240R2）＋Usin H：W=R－Rcos H＋S2÷（24R）－Ucos H：△F＜0＝=＞G＜O＝=＞X=V：Y=W：≠=＞X=T＋Tcos F－Vcos F－Wsin Abs F：Y=Tsin Abs F －Vsin Abs F＋Wcos F：△≠=＞G＜O＝=＞X=V：Y=－W：≠=＞X=T＋Tcos F－Vcos F－Wsin Abs F：Y=Vsin Abs F－Tsin Abs F－Wcos F：△△C=A＝=＞D≥B＝=＞P=90：≠=＞P=270：△≠=＞C＜A＝=＞P= tan-1（（D－B）÷（C－A））＋180：≠=＞P= tan-1（（D－B）÷（C－A））：△△I“COM——N”=A＋√（X2＋Y2）cos（P－tan-1（Y÷X））▲J“COM——E”=B＋√（X2＋Y2）sin（P－tan-1（Y÷X））▲G=G＋20Goto 0程序运行输入：JD？——交点桩号R？——圆曲线半径LS？——缓和曲线长度PJ——L－R＋？——偏角右偏为正ZH——N？——ZH点N坐标ZH——E？——ZH点E坐标JD——N？——JD点N坐标JD——E？——JD点E坐标输出：ZH= ——ZH点桩号HY= ——HY点桩号QZ= ——QZ点桩号YH= ——YH点桩号HZ= ——HZ点桩号T= ——切线长L= ——曲线长E= ——外矢距再输入：ZHUANG HAO？——输入任一点桩号I——B OUT＋IN－？——距中线距离（外侧为正）再输出：COM——N= ——计算出N坐标COM——E= ——计算出E坐标以此循环——再输入——再输出程序中各字母所代表的含义：A——ZH点N坐标B——ZH点E坐标C——JD点N坐标D——JD点E坐标E——外矢距F——偏角（右偏为正）G——弯道内任一点桩号H——弯道内任一点切线角I——计算出的N坐标J——计算出的E坐标K——交点桩号L——曲线长M——ZH点桩号N——HY点桩号O——QZ点桩号P——ZH--JD方位角Q——HZ点桩号R——圆曲线半径S——缓和曲线长度T——切线长U——距中线距离（外侧为正）V——支距W——支距X——弯道内坐标Y——弯道内坐标Z——任一点曲线长度Z[1]——YH点桩号。

786635957WOSHIchunzi1314CASIOfx-4500PA计算器程序和操作一、水准测量高程计算程序1、程序名：H-（1）第一行：L1 lbl 0 第二行：L2 Z：A：｛BK｝第三行：L3 H=Z＋A－B◢第四行：L4 V=K-H ◢第五行：L5 Goto 02、操作步骤①开机：按AC键：②呼出文件名：按FILE键，显示H-（1）文件名为止：③按EXE键，显示Z？，输入水准点BM后视点高程：如129.919；④按EXE键，显示A？，输入后视水准尺读数：如1.633；⑤按EXE键，显示B？，输入前视水准尺读数：如4.186；⑥按EXE键，显示H=127.366，此结果是K3＋300中桩地面髙程；⑦按EXE键，显示K？，输入K3＋300中桩设计高程127.032；⑧按EXE键，显示V=K-H，-0.334此结果是K3＋300中桩挖出数据；⑨按EXE键，显示B？，连续算下去……（－挖；＋填）。

二、坐标反算程序①1、程序名；1ZBFS第一行：L1 AB：Fixm：｛CD｝第二行：L2 Pol(C-A，D-B）◢第三行: L3 W<0 =〉W=W+360第四行：L4 lntW＋0.01lnt（）＋0.006 Frac（60 Frac W）◢2、操作步骤①开机：按AC键：②呼出文件名：按FILE键，显示1ZBFS文件名为止：③按EXE键，显示A？，输入第一点x坐标：如363.567（也可当测站）；④按EXE键，显示B？，输入第一点y坐标：如814.454（也可当测站）；⑤按EXE键，显示C？，输入第二点x坐标：如406.260 （也可作后视点）；⑥按EXE键，显示D？，输入第二点y坐标：如1029.145（也可作后视点）；⑦按EXE键，显示Pol(C-A，D-B）：218.89476为两点间距离；⑧按EXE键，显示lntW＋0.01lnt（60 Frac W）＋0.006 Frac（60 Frac W）78.45108两点方位角；(即是78°45′11″)三、坐标正算程序1、程序名；2ZBZS第一行：L1 AB：Lbl5：｛ST｝第二行：L2 C=A＋Rce（S，T）◢第三行：L3 D=B＋W◢第四行：L4 Goto52、操作步骤①开机：按AC键：②呼出文件名：按FILE键，显示2ZBZS文件名为止：③按EXE键，显示A？，输入第一点x坐标：如363.567；④按EXE键，显示B？，输入第一点y坐标：如814.454；⑤按EXE键，显示S？，输入两点间距离：如218.89476；⑥按EXE键，显示T？，输入两点方位角：如78°45′10.8″；⑦按EXE键，显示C=A＋Rce（S，T），为第二点x坐标：如406.260；⑧按EXE键，显示D=B＋W，为第二点y坐标：如1029.14499；四、坐标反算程序②1、程序名；ZF第一行：L1 Lbl 0第二行：L2 A：B：｛CD｝：C≤0 G=〉oto2第三行：L3 X=C－A第四行：L4 Y=D－B第五行：L5 S=Pol（X，Y）◢第六行：L6 T=W第七行：L7 T＜0=〉第八行：L8 T=W ◢第九行：L9 Goto 0第十行：L10 Lbl 1第十一行：L11 T″T″=360＋T◢第十二行：L12 Goto 0第十三行：L3 lbl 2第十四行：L14 ｛AB｝第十五行：L15 Goto 02、操作步骤①开机：按AC键：②呼出文件名：按FILE键，显示ZF文件名为止：③按EXE键，显示A？，输入第一点x坐标：如363.567 ；④按EXE键，显示B？，输入第一点y坐标：如814.454；⑤按EXE键，显示C？，输入第二点x坐标：如354.618；⑥按EXE键，显示D？，输入第二点y坐标：如553.341；⑦按EXE键，显示Pol(C-A，D-B）S=Pol（X，Y）：261.2663为两点间距离；⑧按EXE键，显示T=268.0371接着按SHIFT °′″键显T=268°2′13.55″为两点的方位角；⑨按EXE键，显示C？，输入下一点x值；⑩重复以上操作方法（略）五、线路直线段点坐标计算程序(1) ㈠本可计算前切线YZ至ZY点(或后切1、程序名；Z-X,Y(1) 线ZY至YZ)间直线上任点中、边桩坐标。

测量程序1.ZBFS（坐标反算）LbI0：/M“CZ:X”N“CZ:Y”/{XY}/U=X-M:V=Y-N/Prog“FW”：I“DIST”▲J“FWJ”▲Goto 0/2.DXJS（导线计算）LbI0：/{BD}/F“FWJ”M“CZX”N“CZY”A“JC”H“ZS”P“XC”Q“YC”L“CD”C/J=F+B-180-A/H:Prog“JD”:J“FWJ”▲U=DcosJ▲V=DsinJ▲X=M+U-P*D/L▲Y=N+V-Q*D/L▲C=0＝＞Goto1△/Goto 0/Lb1 1：/{BD}/J=J+B-180-A/H:Prog“JD”:J“FWJ”▲U=DcosJ ▲V=DsinJ ▲X=X+U-P*D/L▲Y=Y+V-Q*D/L▲Goto 1/3.QXJS（曲线坐标计算）C“C=-1＝＞L,C=1＝＞R”/Prog“QXYS”/Prog“QXZD”/LbI0：/{KE}/Prog“QXZB”/Goto 0/4.LICHPJ（里程计算）C“ZJL:C =-1,R:C=1”/Prog“QXYS”/o“JD:DK”M“JD:X”N“JD:Y”/Prog“QXZD”/Z[15]=Z[10]+Z[3]CosF+(R+Z[1])Cos(F+90C)/Z[16]=Z[11]+Z[3]SinF+(R+Z[1])Sin(F+90C) /W=L: Prog“HHQX”/J=F+JC: Prog“JD”:Z[17]=Z[10]+I CosJ:Z[18]=Z[11]+ISinJ/U=Z[17]-Z[15]:V=Z[18]-Z[16]/Prog“FW”:Z[19]=J:J=J+(A-Z[5]-Z[6])C:Prog “JD”/Z[20]=J: “STAR”▲LbI0：/{XY}:X“POINT:X”:Y“POINT:Y”/U=X-Z[15]:V=Y-Z[16]: Prog“FW”:W=J/ J=C(J-Z[19]):Prog“JD”:Z[21]=J/J=C(W-Z[20]):Prog“JD”:Z[22]=J/Z[21]>180=> Goto1△/Z[22]≤180=> Goto2:≠=>“HY-YH”▲Z[23]=O-Z[7]+L+πRZ[21]÷180▲Z[24]=-C(I-R)▲Goto0/LbI 1:/U=X-Z[10]:V=Y-Z[11]: Prog“FW”:J=J-F: Prog“JD”/Z[25]=ICosJ/Z[25] ≤0=>“≤ZH”▲Z[23]=O-Z[7]+Z[25]▲Z[24]=ISinJ▲△/Z[25]>0=>Z=0: Goto3△/Goto0/LbI 2:/U=X-Z[13]:V=Y-Z[14]:Prog“FW”:W=J:J=Z[12]+180: Prog“JD”/J=W-J: Prog“JD”:Z[25]=ICosJ:Z[25]≥0=>“≥HZ”▲Z[23]=O-Z[7]+Z[9]+Z[25]▲Z[24]=ISinJ▲△/Z[25]<0=>Z=0: Goto4△/Goto0/LbI 3:/Z=Z+Z[25]:W=Z: Prog“HHQX”:Q=J/J=F+QC:Prog”JD”:Z[23]=Z[10]+ICosJ:Z[ 24]=Z[11]+ISinJ:J=J+2QC:Prog“JD”:W=J:U=X-Z[23]:V=Y-Z[24]:Prog“FW”/J=J-W:Prog”JD”:Z[25]=ICosJ:AbsZ[25]>0 .001=> Goto3△/“ZH-HY”▲Z[23]=O-Z[7]+Z+Z[25]▲Z[24]=ISinJ▲Goto0/LbI 4:/Z=Z-Z[25]:W=Z:Prog“HHQX”/Q=J:J=Z[12]-QC: Prog“JD”/Z[23]=Z[13]+I CosJ/Z[24]=Z[14]+I SinJ/J=J-2QC: Prog“JD”:J=J+180: Prog “JD”:W=J/U=X-Z[23]:V=Y-Z[24]: Prog“FW”:J=J-W: Prog“JD”/Z[25]= I CosJ/Abs Z[25]>0.001=> Goto4△/“YH-HZ”▲Z[23]=O-Z[7]+Z[9]-Z+Z[25]▲Z[24]=I SinJ▲Goto0/5.ZBXTZHLbI0:/ASTKLMN/{XY}:V=T-L:U=S-K/Prog“FW”:G=J/X:U=X-K:V=Y-L: Prog“FW”/J=J-G: Prog“JD”/J=A+J-180: Prog“JD”/O=M+ICosJ▲P=N+ISinJ▲Goto0/6.ZXDKFYU“CZX”V“CZY”/Z[26]=U:Z[27]=V/C“C=-1=>L,C=1=>R”/Prog“QXYS”/Prog“QXZD”/LbI0:/{KE}/Prog“QXZB”/U=X-Z[26]:V=Y-Z[27]/Prog“FW”/I“DIST”▲J“FWJ”▲Goto0/9.QXYS（曲线要素）RL“L01”S“L02”A“PJ”/Z[1]=L2÷24÷R-L4÷2688÷R3/Z[2]=S2÷24÷R-L4÷2688÷R3/Z[3]=L÷2-L3÷240÷R2/Z[4]=S÷2-S3÷240÷R2/Z[5]=90L÷π÷R/Z[6]=90S÷π÷R/Z[7]=(R+ Z[2]-( R+Z[1]) cosA)÷sinA+ Z[3]/Z[8]=(R+ Z[1]-( R+ Z[2]) cosA)÷sinA+ Z[4]▲Z[9]=L+S+πR(A-Z[5])-Z[6])÷180▲10.JD（角度处理）J<0=>J=J+360△/J≥360=>J=J-360△/11.QXZD（曲线主点坐标）M“JD:X”N“JD:Y”F“ZH-JD:FWJ”O “JD:DK”/Z[10]=M-Z[7]CosF:Z[10] “ZH:X”▲Z[11]=N-Z[7]SinF:Z[11] “ZH:Y”▲J=F+CA: Prog“JD”:J=J-180: Prog “JD”:Z[12]=J/Z[13]=M-Z[8]CosJ:Z[13] “HZ:X”▲Z[14]=N-Z[8]SinJ:Z[14] “HZ:Y”▲12.HHQX（缓和曲线）U=W-W5÷(40R2L2)/V= W3÷(6RL)- W7÷(336R3L3)/Prog“FW”/13.FW（方位角反算）I=0:J=0/Pol(V,U):J=90-J/Prog“JD”/14.QXZB（曲线坐标）P=K-O+Z[7]/P≤0=> Goto1: ≠=>P≤L=> Goto2: ≠=>P≤Z[9]-S=> Goto3: ≠=>P≤Z[9] => Goto4:≠=> Goto5△△△△/LbI1:/“≤ZH”▲E:J=F+90:Prog“JD”:G=J:/X=M-(Z[7]-P)CosF+ECosG▲Y=N-(Z[7]-P)SinF+ESinG▲Goto0/LbI2:/“ZH-HY”▲E/W=P: Prog“HHQX”:Q=J:J=F+3QC-90C: Prog“JD”:G=J/J=F+QC: Prog“JD”:T=J/X=Z[10]+I CosJ-C×ECosG▲Y=Z[11]+I SinJ-C×ESinG▲G▲Goto0/LbI3:/“HY-YH”▲E/J=F+90C: Prog“JD”:T=J/J=T+180+Z[5]C: Prog“JD”/J=J+180(P-L)C÷R÷π: Prog“JD”:G=J/ X=Z[10]+Z[3]CosF+(R+Z[1])CosT+(R-CE) CosG▲Y=Z[11]+Z[3]SinF+(R+Z[1])SinT+(R-CE) SinG▲G▲Goto0/LbI4:/“YH-HZ”▲E/W=Z[9]-P:B=L:L=S:Prog“HHQX”:L=B:Q=J/J=Z[12]-3QC+90C: Prog“JD”:G=J/J=Z[12]-QC: Prog“JD”:T=J/X=Z[13]+I CosT-C×ECosG▲Y=Z[14]+I SinT-C×ESinG▲G▲Goto0/LbI5:/“≥HZ”▲E:J=Z[12]-90: Prog“JD”:G=J/X=M-(P-Z[9]+Z[8])CosZ[12]-CECos(Z[12] +90C)▲Y=N-(P-Z[9]+Z[8])SinZ[12]-CESin(Z[12]+9 0C)▲LbI0:/说明：“/”为回车键。

线路任意点坐标计算、及任意点对应桩号，左右偏距计算程序卡西欧4800、4850系列计算器测量计算程序一、字母含义；K: 表示拟计算位置的线路桩号;H: 表示计算位置距路线中心线的偏距，左偏为正，右偏为负。

T: 各段线路上作为起算点处的切线方位角。

S: 拟计算点到起算点的曲线长。

L: 在圆曲线上表示曲线长，在缓和曲线上表示缓和段长度，在直线上为零。

R: 表示曲线半径，左偏为正，右偏为负E、F: 起算点的坐标值。

M”X1”N”Y1”: 已知点坐标，求其对应位置桩号及左右偏距。

二、程序1、坐标计算(COORD){K,H}:KH“L+，R-”:Prog “DATA”:”(X,Y)=”:X=X+HSinW :Pause 1:Y=Y-HcosW:2、坐标反算线路桩号(FS ZH)Fix 4:M”X1”N”Y1”:LbI 1:Prog “DATA”:PoI (M-X,N-Y):Fixm:J<0=>J=J+360:≠>J=J⊿Abs(Sin(W-J ))=1=>”K=”:K: Pause 1 :“L+,R- =”: H=ISin(W-J):≠>K=K+Icos(W-J): Goto 13、子程序Prog “DATA”K<（第１段与第２段线路分界处的路线桩号）＝> T=(第１段起算点处的切线方位角值):Ｓ＝Ｋ－（第１段起算点处的路线桩号）:Ｌ＝（在圆曲线上等于Ｓ；在缓和曲线上等于缓和段长度；在直线上为零。

）:R=（曲线半径）:E=（第1段起算点的X坐标值）:F=（第1段起算点的Y 坐标值）≠ >K<(第2段与第3段线路分界处的路线桩号):=>T=(第2段起算点处的切线方向角值):S=K-（第2段起算点处的路线桩号）:L=（在圆曲线上等于S；在缓和曲线上等于缓和段长度；在直线上为零。

）:R=（曲线半径）:E=（第2段起算点的X坐标值）F=（第2段起算点的Y坐标值）≠ >K<(第3段与第4段线路分界处的路线桩号):=>T=(第3段起算点处的切线方向角值):S=K-（第3段起算点处的路线桩号）:L=（在圆曲线上等于S；在缓和曲线上等于缓和段长度；在直线上为零。

计算细那么1、坐标计算：X 1=X+Dcosα,Y1=Y+Dsin α。

式中Y 、 X 为坐标， D 为两点之间的距离，Α 为方位角。

2、方位角计算：1〕、方位角 =tan=两坐标增量的比值，然后用计算器按出他们的反三角函数〔±号判断象限〕。

2〕、方位角： arctan〔 y2- y1)/(x2-x 1)。

加减 180〔大于 180 就减去 180〔还大于 360 就在减去 360〕、小于 180 就加 180 如果 x 轴坐标增量为负数，那么结果加 180°。

如果为正数，那么看 y 轴的坐标增量，如果 Y 轴上的结果为正，那么算出来的结果就是两点间的方位角，如果为负值，加360°。

S=√（y2- y1)+(x2-x 1),1)、当 y2- y1>0，x2-x 1>0 时；α =arctan〔 y2- y1)/(x2-x 1)。

2)、当 y2- y1<0，x2-x 1>0 时；α =360° +arctan〔y2- y1)/(x2-x 1)。

3)、当 x2-x 1<0 时；α =180° +arctan〔y2- y1)/(x2-x 1)。

再用两点之间的距离公式可算距离(根号下两个坐标距离差的平方相加〕。

拨角： arctan〔y2- y1)/(x2-x 1)1、例如：两条巷道要互相平行掘进的话，求它们的拨角：方法〔前视边方位角减后视边方位〕在此后视边方位要加减 180°，假设拨角结果为负值为左偏“逆时针〞〔 +360°就可化为右偏，正值为右偏“顺时针〞。

2、在图上标识方位的方法：就是导线边与Y 轴的夹角。

3、高程计算：目标高程 =测点高程 +?h〔高差〕 +仪器高—占标高。

4、直角坐标与极坐标的换算：〔直角坐标用坐标增量表示；极坐标用方位角和边长表示〕1〕、坐标正算〔极坐标化为直角坐标〕一个点的坐标及该点至未知点的距离和方位角，计算未知点坐标方位角，知A(Xa,Ya) 、Sab、αab，求 B(Xa,Ya)解： ?Xab=Sab×COSαab 那么有 Xb=Xa+?Xab ?Yab=Sab × SIN αab Yb=Ya+?Yab2)、坐标反算，两点的坐标，求两点的距离〔称反算边长〕和方位角(称反算方位角〕的方法A(Xa,Ya) 、 B(Xb,Yb), 求α ab、 Sab。

5800计算器全线坐标计算放样程序（修改版）“XLZBJSCX” ◢ （第一个程序）Lb1 0 ↙Cls : Fix 4 : 30→Dimz ↙“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=":J►DMS ◢ ( 后视方位角 ) Lbl 1 ↙ ( If （如果的意思） And （和字的意思） Then （然后的意思） )“K=”?K ◢ （计算里程） ( 下面有色的是数据库 )If （如果） K< 51760.052 本曲线缓直点桩号 And （和）K≥ 51048.785 上一个曲线缓直点桩号： Then （然后）本曲线缓直点桩号51760.052 →Z[1] : 上一个曲线缓直点桩号51048.785 →Z[2] ： 1 -1 →O（注：左偏曲线输入 - 1→O, 右偏曲线输入1→O） : 偏角12 ’ 23 ’ 19.5 ’ →A ：半径3289.486 →R : 第一缓和曲线长度0 →Z[6] : 第二缓和曲线长度0 →Z[7] : 交点 X 坐标→B : 交点 Y 坐标→C : 小里程向交点方位角→E : 交点向大里程方位角→F : Goto 2 : IfEnd ↙ ( 重兰字是输入的数字 ) ………… （曲线段分段输入）补充直线段输入如下If （如果） K< 本段直线终点里程 And （和）K≥ 本段直线起点里程 :Then （然后）1→O: 本段直线终点里程→Z[3]: 终点坐标X→Z[16]: 终点坐标Y→Z[17]: 方位角→E:Goto 4:IfEnd ↙Lb1 2 ↙ （曲线要素计算）Z[6] ÷ 2- Z[6]^ 3 ÷ (240R^2)+ Z[6]^ 5 ÷ (34560*R^4) →Z[8] ↙ （ M1 不输）Z[7 ] ÷ 2- Z[7]^ 3 ÷ (240R^2)+ Z[7]^ 5 ÷ (34560 R^4) →Z[9] ↙ （ M2 不输）Z[6]^ 2 ÷ (24R)- Z[6]^ 4 ÷ (2688R^3) →Z[10] ↙ （ P1 不输）Z[7]^ 2 ÷ (24R)- Z[7]^ 4 ÷ (2688R^3) →Z[11] ↙ （ P2 不输）(πAR ) ÷ 180+0. 5 × ( Z[6]+ Z[7])→S ↙ （曲线总长）90 × Z[6 ] ÷ ( R × π) →Z[14] ↙ （第一缓和曲线总偏角）( ×÷以后自己改 )9 0 × 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]-S→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] A nd K> Z[3] : T hen Goto 5 : IfEnd ↙If K≤Z[5] A nd K> Z[4] : T hen Goto 6 : IfEnd ↙If K≤Z[1] A nd K> Z[5] : T hen Goto 7 : IfEnd ↙LB1 4 ↙ （里程小于直缓点直线独立坐标）K- Z[3] →X : 0→Y : E→T : Prog“TYZBCX” ：Goto 1 ↙Lb1 5 ↙ （第一缓和曲线独立坐标）K- Z[3] →H ↙H-H^5/(40*R^2* Z[6]^2)+H^9/(3456*R^4* Z[6]^4) →X ↙H^3/(6*R* Z[6])-H^7/(336*R^3* Z[6]^3) →Y ↙90*H^2/( R*π* Z[6]) →T ↙IF O >0 ：Then T +E→T : Else E-T →T : T<0=>360+T→T : IfEnd ↙PRrog“TYZBCX” ：Goto 1 ↙Lb1 6 ↙ （圆曲线独立坐标）K- Z[4] →H ↙H*180/( R*π)+ Z[14]→T ↙R*SIN( T)+ Z[8]→X ↙R*(1-COS (T))+ Z[10]→Y ↙IF O >0 ：Then T +E→T : Else E-T →T : T<0=>360+T→T : IfEnd ↙Prog“TYZBCX” ：Goto 1 ↙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)→X ↙Z[13]*SIN( A)-U*SIN( A)+V*COS (A)→Y ↙IF O >0 ： Then F-T→T : T<0=>360+T→T : Else F+T →T : IfEnd ↙Prog“TYZBCX” ：Goto 1 ↙子程序：“TYZBCX” ↙ （统一坐标计算）（第 2 个程序）IF O<0 : Then -Y→Y : IfEnd ↙“QXJ=” :T ◢ （计算里程点切线方位角，可以不显示）Z[16]+X*COS (E)-Y*SIN( E)→Z[18] ↙Z[17]+X*SIN (E ） +Y*COS (E ）→Z[19] ↙“XI=” : Z[18] ◢ （ XI 中线 X ）“ YI =” : Z[19] ◢ （ YI 中线 Y ）Pol(Z[18]-M,Z[19]-N):"DI=":I ◢ （中桩放样距）J<0=>J+360→J:"FI=": J ►DMS ◢ （中桩放样方位角）“ PJ =”?P ◢ ( 输入边桩与线路夹角 PJ ，左偏– 90 右 +90 )“ PD =”?D ◢ （输入边桩距 PD ）Z[18]+D*COS(T+P) →Z[20] ↙Z[19]+D*SIN(T+P) →Z[21] ↙“XP=”: Z[20] ◢ （ XP 边桩 X ）“YP=”: Z[21] ◢ （ YP 边桩 Y ）Pol(Z[20]-M,Z[21]-N):"DP=":I ◢ （边桩放样距）J<0=>J+360→J:"FP=":J ►DMS ◢ （边桩放样方位角）Return ↙ （以上都要输进计算器）注解不输K 里程 XI 中线 X YI 中线 Y PD 输入边桩距 XP 边桩 X 坐标YP 边桩 Y 坐标 I 边桩放样距 PJ 输入边桩与线路夹角，左偏– 90 右 +90 ) I 中桩放样距 T 计算里程点切线方位角，可以不显示卡西欧FX5800全线贯通万能正、反算程序FX5800计算器的积分程序（正反算、全线贯通、新线路）终极版ZHUCHENGXU 主程序"1.ZS,2.FS" ?→Q输入1正算，输入2反算“NEW=0,OLD≠0”?ZIf Z=0:Then “X0=”?A:“Y0=”?B:“C0=”?C:“1/R0=”?D:“1/RI=”?E:“SP=”?F:“EP=”?G:Ifend:Q=2=>Goto 2Lbl 1 :“KM=,<0Stop”?H:H<0=>Stop:“PJ=”?O:“PY=”?LLbl Z:Z=1=> Prog“01”：Z=2=> Prog“02”选择数据库文件，可增加H- F→X:0.5(E-D)÷(G-F)→NC+(XD+NX2)*180÷π→P：P<0=>P+360→P：P>360=>P-360→PA+∫(cos(C+(XD+NX2)*180÷π),0,X)+Lcos(P+O)→UB+∫(sin(C+(XD+NX2)*180÷π),0,X)+Lsin(P+O)→VQ=2=>Goto 4：Cls：Fix 3"Xn=":Locate 4，1，U："Yn=": Locate 5，2，V：“FWJ=”：P▶DMS◢Norm 2：Cls：Goto 1Lbl 2:“XD=,<0,STOP”？R：R<0=>Stop:“YD=”？S“KMDG=”？H ：90→O:0→L：Goto Z （H线路范围内的任意桩号）Lbl 4:Pol（R-U，S-V）：J<0 => J+360→JWhile abs（Icos（J-P））≤0.001:P-J>180=> J+360→J： P-J<-180=> P+360→P:IF P-J>0:then -I→L:else I→L ifendGoto 3: Whileend：H+Icos（J-P）→H:Goto ZLbl 3:Cls：Fix 3“KM=”: Locate4，1，H：“PY=”: Locate4，2，L◢Norm 2：Cls：Goto 201（数据库子程序）If H<=第一曲线终点桩号：then 第一曲线起点X→A：第一曲线起点Y→B：第一曲线起点方位角→C：起点曲率→D：终点曲率→E：起点桩号→F：终点桩号→G：return：ifend……………程序说明：1、该程序可以计算任意线形（直线、圆曲线、缓和曲线、不完整曲线）任意桩号的坐标（正算，输入1），也可根据坐标计算该点到线路的距离及垂足桩号（反算，输入2）；2、（NEW=0，OLD≠0）？如果要计算的点为数据库线路中的点，则输入数据库编号（以整数1、2、3…代替输入）；如果在数据库中没有要计算线路的数据，则输入曲线要素X0：曲线起点X坐标；Y0：曲线起点Y坐标；C0:曲线起点方位角；R0-1、 RI-1:曲线起点、终点曲率，直线为0，曲线左偏输入负值，右偏输入正值；SP、 EP：曲线起点桩号，终点桩号；KM：待求点桩号；PJ：正斜交的设定；PY：偏中距离，线路上的点输入0，右偏输入+值，左偏输入-值；3、正算显示坐标及切线方位角；反算输入线路的任意桩号（此桩号越接近真实值计算速度越快）、待求点坐标，显示待求点桩号及偏中距离；4、正算子程序为积分公式编写而成；反算子程序为角度趋近的方法编写，计算速度有点慢。