4800、4850计算程序坐标转换

高速MC公路测量CASIO4800&4850万能坐标计算程序(完整版)

高速公路测量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、编制说明本程序是运用复化辛普生公式根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距，对该曲线段范围内任意里程中边桩坐标进行计算，以及对卡西欧扩充变量的灵活应用，实现了真正意义上的的全线贯通及曲线要素输入程序化（在不修改程序内容的情况下可通过运行程序输入任意多段曲线元要素）。

4850 4500 4800正反算程序

1．中边桩计算程序<1>主程序名—XYLbl 1:{KM}：K：M：Prog”3”:X=C+Mcos(W-90:Y=D+Msin(W-90:”X=”:X:Pause 0:”Y=”:Y▲Goto1<2>子程序—3PrOg”QX”：Prog”20”<3>曲线要素—QX数据实例：某高速公路中线设计如下图所示K<2000=>B:E:F:O:P:Q:R: K≤4570=>B=25:E=1000:F=1000:O=100P=66:Q=4000:R=800: K≤5040=>B=-17:E=1086.4378:F=1555.1451:O=90:P=91:Q=4570:R=800: K≤5670=>B=35:E=1168.2556:F=2016.0909:O=120:P=74:Q=5040:R=1000: B=-25:E=1180.9927:F=2636.2334:O=150:P=109:Q=5670:R=1500注：无缓和曲线时，O=0.00001当输入里程小于2000时，为新曲线要素输入，可修改值域范围<4> 计算过程—20Z=O：Prog”8”：A=B/AbsB:L=ABπR/180+O:T=(D+RcosN)tan(AB/2)+C-RsinN：K≤Q-T=>C=(K-Q)cosP+E：D=(K-Q)sinP+F:W=P： K≤Q-T+O=>Z=K-Q+T：Prog”9”：K≤Q-T+L-O=>Z=O: Prog”9”:V=90(K-Q+T-O)/ πR:S=2R sinV:N=W+AV:C= C+S cosN:D=D+S sinN:W=W+2AV: Z=Q-T+L-K: Prog”8”:V=P+B-180-AV:W=P+B-AN :C=S cosV+E+T cos(P+B:D=S sinV+F+T sin(P+B: Z=K-Q+2T-L:W=P+B:C=E+ZcosW:D=F+Z sinW△△△△N=Z[3]<5> 缓和曲线方程式—8V=RO:C=Z-Z^5/40V2+Z^9/3456V^4:D=Z^3/6V-Z^7/336V^3+Z^11/42240V^5:N= 90Z^2/πV:V=tan-1(D/C:S=(D^2+C^2)^0.5<6> 子程序—9Prog”8”:C=Scos(P+AV)+E-TcosP:D=Ssin(P+AV)+F-TsinP:W=P+AN2.桩号、边桩距返算程序—KM方式一：试算法（速度慢，程序简单）Fixm:K=Q-T+L-O:Prog”20”Lbl 1:{XY}Lbl 2:POL(X-C,Y-D:K=Isin(J-W+90)+K:Abs sin(J-W+90)<0.0001＝>”K=”:K:Pause0:”M=”:M=Icos(J-W+90▲≠>Prog”20”:Goto2△方式二：解算法（程序复杂，计算速度快）Defm5: Lbl 1:Fixm:{XY}:Pol(X-E,Y-F:Z[1]=I:Z[2]=J:cos(J-P<0＝>K=Q+Icos(J-P: ≠>K=Q+Icos(J-P-B)+2T+L△K≤Q-T＝>M=Icos(J-P+90:Goto5≠>Goto2△△Lbl 2:Z=O:Porg”8”:U=(D+RcosN)+cos(B/2):Z=P+B/2+90A:C=UcosZ+E:D=UsinZ+F:Pol(X-C,Y-D:Asin(J-P-AN+90 A<0＝>K=K-AZ[1]cos(Z[2]-P+90)tan(90 (K-Q+T)^2/πRO):Goto3:≠>Asin(J-P-B+AN+90A>0 ＝>K=K-AZ[1]cos(Z[2]-P-B-90 )tan(90(K-Q+T-L)^2/πRO):Goto4:≠>K=Q-T+O+Asin-1(sin(J-P-AN+90 A))×πR/180 :M=A(I-R:Goto5△△Lbl 3:Z=K-Q+T:Prog”9”Pol(X-C,Y-D:K=K+Isin(J-W+90 :Abs sin(J-W+90<.001＝>M=Icos(J-W+90 :Goto5: ≠>Goto3△Lbl 4:Z=Q-T+L-K:Prog”8”:V=P+B-180-AV:W=P+B-AN:C=cosV+E+Tcos(P+B:Pol(X-C,Y-D:K=K+Isin(J-W+90 :Abs sin(J-W+90 <.001＝>M=Icos(J-W+90:Goto5: ≠>Goto4△Lbl 5:”K=”:K:Pause0:”M=”:M▲注：程序“20”为中桩计算程序Q—直缓点里程；T—切线长，L—曲线长；O—缓和曲线长；W—切线方位角Pol—极坐标转换，为涵数特别注意数字0和字母O的区别一、程序功能本程序由一个主程序(TYQXjs)和两个子程——正算子程序(SUB1)、反算子程序(SUB2)序构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。

测量程序4800和4850

1.主程序(ZCX1)"1.SZ => XY"："2.XY => SZ"：N：U"X0"：V"Y0"：O"S0"：G"F0"：H"LS"：P"R0"：R"RN"：Q：C=1÷P：D=(P-R)÷(2HPR)：E=180÷π：N=1=>Goto 1：≠>Goto 2Δ换行Lbl 1：{SZ}：SZ：W=Abs(S-O)：Prog "SUB1"：X"XS"=X◢Y"YS"=Y◢F"FS"=F-90◢Prog“ZDS”2.Lbl 2：{XY}：XY：I=X：J=Y：Prog "SUB2"：S"S"=O+W◢Z"Z"=Z◢Prog“ZDS”2. 正算子程序(SUB1)A=0.1739274226：B=0.3260725774：K=0.0694318442：L=0.3300094782：F=1-L：M=1-K：X=U+W(Acos(G+QEKW (C+KWD))+Bcos(G+QELW(C+LWD))+Bcos(G+QEFW(C+FWD)) +Acos(G+QEMW(C+MWD)))：Y=V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+QEFW(C+FWD))+ Asin(G+QEMW(C+MWD)))：F=G+QEW(C+WD)+90：X=X+ZcosF：Y=Y+ZsinF换行3. 反算子程序(SUB2)T=G-90：W=Abs((Y-V)cosT-(X-U)sinT)：Z=0：Lbl 0：Prog "SUB1"：L=T+QEW(C+WD)：Z=(J-Y)cosL-(I-X)sinL：AbsZ<1E-6=>Goto1：≠>W=W+Z：Goto 0Δ换行Lbl 1：Z=0：Prog "SUB1"：Z=(J-Y)÷sinF换行Lbl 0：S≤终点桩号=>O=起点桩号：U=起点X：V=起点Y：G=起点在曲线上的方位角：H=曲线或直线长：P=起点半径（直线段和缓和曲线为无穷大）：R=终点半径（直线段和缓和曲线为无穷大）：Q=路线转向（左“—”右“+”直线为“0”）：Prog “ZCX1”：Goto 0⊿↙Lbl 1：S≤终点桩号=>O=起点桩号：U=起点X：V=起点Y：G=起点在曲线上的方位角：H=曲线或直线长：P=起点半径（直线段和缓和曲线为无穷大）：R=终点半径（直线段和缓和曲线为无穷大）：Q=路线转向（左“—”右“+”直线为“0”）：Prog “ZCX1”：Goto 0⊿↙。

卡西欧4800和4850计算程序

卡西欧4800\4850程序一：坐标反算：(ZBFS)①程序步骤：C”X1”:D”Y1”:E”X2”:F”Y2”:Fixm:Pol(E-C,F-D:I”S1-2=”▲J≤Ｏ＝＞J=J+360△J”A1-2=”②操作过程：ZBFS→EXE→输入X1值→EXE→输入Y1值→EXE→输入X2值→EXE→输入Y2值→EXE→EXE→EXE注：①：S1-2:计算得出的距离。

②：A1-2:计算得出的角度。

（按shift°′″转换为60进制的角度）③：此程序可循环计算。

③案例：已知：D1坐标（7811.23，606.136），D2坐标（7805.915，431.910）求解： a D1-D2(A1-2)=268°15′09.56 ″S D1-D2(S1-2)=174.3071二：坐标正算：(ZBZS)①程序步骤：C”X1”:D”Y1”:S”S1-2”:R”A1-2”:Fixm:X”X2”=C+Rec(S,R) ▲Y”Y2”=D+J②操作过程：ZBZS→EXE→输入X1值→EXE→输入Y1值→EXE→输入S1-2距离值→EXE→输入A1-2角度（例268°15′09.56 ″）值→EXE→EXE注：①：X2计算得出坐标，Y2计算得出坐标。

②：此程序可循环计算。

③案例：已知：D1坐标（7811.23，606.136），两点距离：S D1-D2(S1-2)=174.3071方位角：a D1-D2(A1-2)=268°15′09.56 ″求解：D2坐标（7805.915，431.910三：线路曲线计算程序:（XLQXJS）(万能公式)①程序步骤：1.Lbl 0:{E G}：A“XA”：B“YA”：C“CA”：D“1/RA”：E“1/RB”：F“DKA”：G“DKA”2.Lbl 1:{H O R}:H“DKI”：O“DL”：R“DR”：H＞G=﹥Goto 3⊿3.P=(E-D)/Abs(G-F):Q=Abs(H-F):I=P×Q:T=D+I4.J=C+(I+2D)×Q×90/π ▲5.M=C+(I/4+2D)Q×45/(2π):N=C+(3I/4+2D)Q×135/(2π):6.K=C+(I/2+2D)Q×45/π7.X=A+Q(cosC+4(cosM+cosN)+2cosK+cosJ)/12 ▲8.Y=B+Q(sinC+4(sinM+sinN)+2sinK+sinJ)/12 ▲9.U“XL”=X+Ocos(J-90) ▲10.V“YL”=Y+Osin(J-90) ▲11.W“XR”=X+Rcos(J+90) ▲12.Z“YR”=Y+ Rsin(J+90) ▲13. A=X：B=Y：C=J：D=T：F=H：Goto 1Lbl 3: A=X：B=Y：D=E：F=G：C=J：Goto 0程序说明：A- 曲线元起点A的坐标；B- 曲线元起点B的坐标；C- 曲线元起点A的切线坐标方位角；F- 曲线元起点A的里程；G- 曲线元起点B的里程；H- 曲线上待求点i的里程；D- 曲线元起点A的曲率；E- 曲线元终点B的曲率；XL-左边线点位X坐标；YL-左边线点位Y坐标；XR-右边线点位X坐标；YR-右边线点位Y坐标；X- 中线点位纵坐标；Y- 中线点位横坐标；DL-左边线距中线平距；DR-右边线距中线平距；该程序需要输入的数据为：（1）.曲线元起点A的坐标及切线坐标方位角，计算器上用“XA”，“YA”，“CA”显示；（2）.曲线元起点A和B的曲率，计算器上用I÷RA，I÷RB显示（曲线左偏时取“-”）；（3）.曲线元起点A和终点B的里程，计算器上用“DKA”，“DKB”显示；（4）.输入待求点里程和该点距左右的水平距离，计算器上用“DKI”，“DL”，“DR”显示；每算完一个待求点的中线及边线坐标，程序又让输入下一点的“DKI”，“DL”，“DR”当输入的“DKI”大于“DKB”时程序中显示“1÷RB”和“DKB”，此时输入下一个曲线元起点的曲率和里程，然后重复步骤（4），即可计算下一个曲线中线及边线点位坐标。

4850计算器使用说明

4850计算器程序使用说明一、XY-DA主程序—坐标反算（已知两点坐标，计算距离及方位角）XCZ？测站点X坐标；YCZ？测站点Y坐标；XHS？后视点X坐标；YHS？后视点Y坐标；DHS=两点间距离；AHS=两点间方位角。

二、GLCL主程序—计算中桩及边桩坐标XCZ？测站点X坐标；YCZ？测站点Y坐标；XHS？后视点X坐标；YHS？后视点Y坐标；DHS=两点间距离；AHS=两点间方位角。

JD？交点号；P？待算点桩号；X=，Y=待算点中桩X、Y坐标；D=待算中桩点到测站点距离；A=待算中桩点到测站点方位角；AO=测站中桩点到待算点与测站点到后视点的夹角；D(L—，R+)？待算里程距离中桩距离（即边距，左为-1，右为+1，输入0为重新进入里程输入）X=，Y=待算点边桩X、Y坐标；D=待算边桩点到测站点距离；A=待算边桩点到测站点方位角；AO=测站边桩点到待算点与测站点到后视点的夹角。

三、BZFY主程序—已知点坐标，计算该点里程桩号和距离该里程中线的距离（即边距）XCZ？测站点X坐标；YCZ？测站点Y坐标；XHS？后视点X坐标；YHS？后视点Y坐标；DHS=两点间距离；AHS=两点间方位角。

JD？交点号；ZHS？后视点高程；VHS？后视点棱镜高；△HHS？后视高差；X、Y测量点X、Y坐标；P=测量点里程桩号；D(L-，R+)测量点边距；△H？前视棱镜高；V？前视高差。

四、ZBJS主程序—已知点坐标及距离方位角计算另一点坐标A？B?已知点X，Y坐标；S？待算点到已知点距离；F？已知点到待算点方位角；X=，Y=待算点X，Y坐标。

五、ZX主程序—直线段计算中、边桩坐标Z？直线段起点桩号；A？，B? 直线段起点X、Y坐标；F？线路方位角；U？待算点里程桩号；ZZX=、ZZY=待算里程中桩X、Y坐标；BJ？待算里程边距；I(L-1，R+1)？-1为左边距，+1为右边距；BZX=、BZY=待算里程边桩X、Y坐标。

六、YQX主程序—计算圆曲线中、边桩坐标Z？圆曲线起点桩号；A？，B? 圆曲线起点X、Y坐标；R？圆曲线半径；F？圆曲线起点方位角；K(L-1，R+1)？曲线左偏为-1，右偏为+1；U？待算点里程桩号；ZZX=、ZZY=待算里程中桩X、Y坐标；I(L-1，R+1)？-1为左边距，+1为右边距；BZX=、BZY=待算里程边桩X、Y坐标。

4850和4800曲线坐标正算程序

★★★缓和曲线、圆曲线、直线段坐标计算程序★★★一：程序适用于计算器CASVO fx-4800P和CASVO fx-4850。

本程序可计算与线路中心成任意夹角的缓和曲线、圆曲线、直线段中、边桩坐标及切线方位角。

二：程序（在CASVO fx-4800P中输入，所有“：“后面都没有空格。

）1 A：R：C“LS”：D“JD”2 P=C∧2/24/R-C∧4/2688/R∧33 Q=C/2-C∧3/240/R∧24 B=90C/兀/R5 T=(R+P)tan(AbsA/2)+Q◢若在“4850”中输入则改为“T=”：T=(R+P)tan(AbsA/2)+Q◢6 W=(R+P)/cos(A/2)-R◢若在“4850”中输入则改为“W=”：W=(R+P)/cos(A/2)-R◢7 L=((AbsA)-2B)兀R/180+2C◢若在“4850”中输入则改为“L=”：L=((AbsA)-2B)兀R/180+2C◢8 G=D-T◢若在“4850”中输入则改为“G=”：G=D-T◢9 H=G+C◢若在“4850”中输入则改为“H=”：H=G+C◢10 I=G+L/2◢若在“4850”中输入则改为“I=”：I=G+L/2◢11 K=G+L-C◢若在“4850”中输入则改为“K=”：K=G+L-C◢12 M=G+L◢若在“4850”中输入则改为“M”：M=G+L◢13 N：E：F：J14 A＜0=>S=-1：≠=>S=1⊿15 U=F+A/2+90S16 V=W+R17 B=N+VcosU18 O=E+VsinU19 Lbl 120 {Z}21 Z≤G=>L=T+G-Z22 V=F+18023 U=F24 Goto 2⊿25 Z≤H=>L=Z-G26 V=L-L∧5/(90R∧2C∧2)27 L=30L∧2S/(兀RC)28 P=F+18029 Q=F+L30 U“FW”=F+3L◢若在“4850”中输入则改为“FW=”：U=F+3L◢31 Goto 4⊿32 Z≤K=>L=F+A/2+90S+180+180(Z-I)S/R/兀33 U“FW”=L+90S◢若在“4850”中输入则改为“FW=”：U= L+90S◢34 Goto 5⊿35 Z≤M=>L=M-Z36 V=L-L∧5/(90R∧2C∧2)37 L=30SL∧2/(兀RC)38 P=F+A39 Q=F+A+180-L40 U=F-3L+A41 Goto4⊿42 Z>M=>L=Z-M+T43 U=F+A44 V=U45 Goto 246 Lbl 247 X=N+LcosV◢48 Y=E+LsinV◢49 Goto 6⊿50 Lbl 351 {W}：若在“4850”中输入则改为{W}：{O}52 P“XL”=X+Wcos(U-90) ◢若在“4850”中输入时不改变，计算结果不显示XL=53 Q“YL”=Y+Wsin(U-90) ◢若在“4850”中输入则改为“YL=”：Q=Y+Wsin(U-O) ◢54 P“XR”=X+Wcos(U+90) ◢若在“4850”中输入则改为“XR=”：P=X+Wcos(U+180-O) ◢55 Q“YR”=Y+Wsin (U+90)◢若在“4850”中输入则改为“YR=”：Q=Y+Wsin (U+180-O)◢56 Goto 157 Lbl 458 X=N+TcosP+VcosQ◢若在“4850”中输入则改为”X=”：X=N+TcosP+VcosQ◢59 Y=E+TsinP+VsinQ◢若在“4850”中输入则改为”Y=”：Y=E+TsinP+VsinQ◢60 Goto 661 Lbl 562 X=B+R cosL◢63 Y=O+RsinL◢64 Goto 665 Lbl 666 J=1=> Goto 3⊿67 Goto 1~~完~~直线正算程序（放样程序）。

线路任意点坐标计算、及任意点对应桩号,左右偏距计算程序.tmp

线路任意点坐标计算、及任意点对应桩号，左右偏距计算程序卡西欧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；在缓和曲线上等于缓和段长度；在直线上为零。

卡西欧4800和4850放样坐标计算程序

已知任意点的里程、半宽，计算该点平面坐标的程序前言1、本计算程序适用范围仅限于某ZH点至下一相邻ZH点间的线路范围。

2、使用本程序前需进行导线点坐标计算并求出各个曲线的曲线要素如第一切线方位角、转向角、曲线半径、缓和曲线长度、切线长度、曲线长度等。

3、使用过程中要注意长短链对里程传递的影响关系，否则将导致计算结果错误。

4、往计算器中输入程序符号时将⊿改为。

5、本程序仅供参考,实际使用时请先进行验证。

6、本程序适用的计算器型号为卡西欧4800/4850。

“wangDongping” (程序名称)Defm 13↙XYVWRLZTC：Fixm↙Z[4]=0.5×L－L3÷(240×R2)↙Z[2]=Z[4]+(R+L2÷24÷R)Abs tan(T÷2)：G=Z＋Ans：Fix3：G＂JDDK=＂◢Z[9]=R Abs Tπ÷180＋L：Z[9]＂QXL=＂◢Z[13]=Z＋Ans：Z[13]＂HZDK=＂◢D=C＋T：D＜0⇒D=D＋3600⊿D＞3600⇒D=D－3600⊿D◢Z[11]＂HZX＂=V＋Z[2]×(CosC＋CosD)◢Z[12]＂HZY＂=W＋Z[2]×(SinC＋SinD)◢Norm：Prog＂WEIXONG2＂(程序“WaangDongping”运行到此时，自动转到程序“WangDangping2”进行下一步运行)“WangDongping2”（程序名称）Lb1 0：{K,S,Q,E,F,N,U,H}：P=Q＋Abs(E－F)×N：K≤Z⇒Goto1⊿Z＜K⇒K≤Z＋L⇒Goto2⊿Z＋L＜K⇒K≤Z＋Z[9]÷2⇒Goto3⊿Z＋Z[9]÷2＜K⇒K＜Z[13]－L⇒Goto4⊿Z[13]－L≤K⇒K＜Z[13]⇒Goto5⊿ K≥Z[13]⇒Goto6⊿Goto0↙Lb1 1：Z[1]=V＋(Z－K)×cos(C＋1800)＋Pcos(C＋S)＋Ucos(C＋H)：Z[3]=W＋(Z －K)×sin(C＋1800)＋Psin(C＋S)+Usin(C+H)：Goto7⊿Goto0↙Lb1 2：Z[5]=K－Z－(K－Z)5÷(40R2L2)：Z[6]=（K－Z）3÷（6RL）：A=T÷AbsT：Z[7]=(K－Z)2×180÷(6πRL)：Z[1]=V＋√((Z[5])2＋(Z[6])2)×cos(C＋A×Z[7])＋Pcos(C＋3×A×Z[7]＋S)＋Ucos(C＋3×A×Z[7]＋H)：Z[3]=W＋√((Z[5])2＋(Z[6])2)×sin(C＋A×Z[7])＋Psin(C＋3×A×Z[7]＋S)＋Usin(C＋3×A×Z[7]＋H)：Goto7⊿Goto0↙Lb1 3：B=L×180÷(2πR)：Z[8]=B＋(K－Z－L)×180÷(πR)：Z[5]=Z[4]+R×sinZ[8]：Z[6]=R+L2÷(24R)－RcosZ[8]：Z[7]=tan-1(Z[6]÷Z[5])：Z[10]=(2×(K－Z)－L)×180÷(2πR)：A=T÷AbsT：Z[1]=V+ √((Z[5])2+(Z[6])2)×cos(C+A×Z[7])+Pcos(C+A×Z[10]+S)＋Ucos(C＋A×Z[10]＋H)：Z[3]=W+√((Z[5])2＋(Z[6])2)×sin(C+A×Z[7])+Psin(C+A×Z[10]+S)＋Usin(D＋A×Z[10]＋H)：Goto7⊿Goto0↙Lbl 4：B=L×180÷(2πR)：Z[8]=B+(Z[13]－K－L)×180(πR)：Z[5]=Z[4]+R×sinZ[8]：Z[6]=R+L2÷(24R)－RcosZ[8]：Z[7]=tan-1(Z[6]÷Z[5])：Z[10]=(2×(Z[13]－K)－L)×180÷(2πR)：A=－T÷AbsT：Z[1]=Z[11]+ √((Z[5])2+(Z[6])2)×cos(D+180+A×Z[7])+Pcos(D+A×Z[10]+S)＋Ucos(D＋A×Z[10]＋H)：Z[3]=Z[12]+√((Z[5])2+(Z[6])2)×sin(D+180+A×Z[7])+Psin(D+A ×Z[10]+S)＋Usin(D＋A×Z[10]＋H)：Goto7⊿Goto0↙Lbl 5：Z[5]=Z[13]－K－(Z[13]－K) 5÷(40R2L2)：Z[6]=(Z[13]－K) 3÷(6RL)：A=－T÷AbsT：Z[7]=(Z[13]－K)2×180÷(6πRL)：Z[1]=Z[11]+ √((Z[5])2+(Z[6])2)×cos(D+180+A×Z[7])+Pcos(D+3×A×Z[7]+S)+Ucos（D+3×A ×Z[7]+H）：Z[3]=Z[12]+√((Z[5])2+(Z[6])2)×sin(D+180+A×Z[7])+Psin(D+3×A×Z[7]+S)＋Usin(D+3×A×Z[7]+H)：Goto7⊿Goto0↙Lb1 6：Z[1]=Z[11]+(K－Z[13])×CosD+Pcos(D+S)＋Ucos(D+H)：Z[3]=Z[12]+(K －Z[13])×sinD+Psin(D+S)＋Usin(D+H)：Goto7⊿Goto0↙Lb1 7：M=Z[1]－X：O=Z[3]－Y：Pol(M,O)：J＜0 J=J+3600⊿Fix3：J◢I◢P◢E－F◢M+X◢O+Y◢Norm：Goto0↙坐标程序符号说明1、X，Y：测站坐标2、V，W：ZH点X,Y坐标3、R：曲线半径4、L：缓和曲线长度5、Z：ZH点里程6、T：曲线转向角（左负右正）7、C：第一切线方位角8、K：放样点里程9、Q：放样半宽（平行于S角非切线边在里程点左右侧的宽度，若已知实际半宽P，则Q=P；若放路基边桩则输入变坡点处的半宽P）10、S：放样点与其对应的中线点连线和该中线点切线的夹角（左负右正，即放样点在线路左侧S为负，同理在右侧为正）11、F：放样点地面标高12、N：路基边坡坡度13、E：变坡点设计标高14、P=Q+Abs（E－F）×N（Abs为绝对值），若P为已知，则Q的值直接输入P的值即可，E、F、N均为零。

4800P~48500P主线坐标计算程序

4800P~4850P全线坐标计算程序子程序“ZB”D≤C=>Goto 1△←┘D＞C=>D≤C+H=>Goto 2△←┘D＞C+H=>D≤C+H+I => Goto 3△←┘D＞C+H+I=>D＜C+H+I+J=>Goto 4△←┘D≥C+H+I+J=>Goto 5△←┘Lbi 1←┘M＞180=>M=M-180∶≠>M=M+180△←┘X=A+(C-D)CosM∶Y=B+（C-D）SinM←┘M＞180=>M=M-180∶≠>M=M+180△←┘Z=M←┘Goto 6△←┘Lbi 2←┘V=D-C：W=RH←┘L=V-V^5/（40W²）+V^9/（3456W^4）-V^13/（599040W^6）←┘N=V^3/（6W）-V^7/（336W^3）+V^11/（42240W^5）-V^15/（9676800W^7）←┘O=tan-1 (N/L):P=√(L ²+N ²)←┘K=1=>Q=M+O△K=0=>Q=M-O:Q＜0=>Q=Q+360△△←┘X=A+P CosQ:Y=B+P SinQ←┘T=90V ²/(∏w)←┘K=1=>Z=M+T:Z＞360=>Z=Z-360△△K=0=>Z=M-T:Z＜0=>Z=Z+360△△←┘Goto 6△←┘Lbi 3←┘L=H-H ^3/(40R ²)+H ^5/(3456R ^4)-H ^7/(599040R^6)←┘N=H²/(6R)-H^4/(336R^3)+H^6/(42240R^5)-H^8/(9676800R^7)←┘O= tan-1(N/L):P=√(L²+N²)←┘K=1=>Q=M+O△K=0=>Q=M-O:Q＜0=>Q=Q+360△△←┘X=A+P Cos Q:Y=B+P SinQ←┘T=90H/(∏R)←┘K=1=>Z=M+T:Z＞360=>Z=Z-360△△K=0=>Z=M-T:Z＜0=>Z=Z+360△△←┘V=D-C-H←┘U=180V/(∏R):W=2R Sin (U/2)←┘K=1=>Z=Z+U/2△K=0=>Z=Z-U/2:Z＜0=>Z=Z+360△△←┘X=X+W CosZ:Y=Y+W SinZ←┘K=1=>Z=Z+U/2:Z＞360=>Z=Z-360△△K=0=>Z=Z-U/2:Z＜0=>Z=Z+360△△←┘Goto 6△←┘Lbi 4←┘V=C+ H+I+J-D:W=RJ←┘L=V-V^5/(40W²)+V^9/(3456W^4)-V^13/(599040W^6)←┘N=V^3/(6W)-V^7/(336W^3)+V^11/(42240W^5)-V^15/(9676800W^7)←┘O=tan-1(N/L):P=√(L²+N²)←┘G＞180=>Q=G-180: ≠>G＜180=>Q=G+180△←┘K=1=>U=Q-O:U＜0=>U=U+360△△K=0=>U=Q+O:U＞360=>U=U-360△△←┘X=E+P CosU:Y=F+P SinU←┘T=90V²/(∏W)←┘K=1=>Z=G-T:Z＜0=>Z=Z+360△△K=0=>Z=G+T:Z＞360=>Z=Z-360△△←┘Goto 6△←┘Lbi 5←┘V=D-C-H-I-J←┘X=E+V CosG:Y=F+V SinG←┘Z=G←┘Goto 6△←┘Lbi 6←┘S＜0=>Z＞90=>Z=Z-90: ≠>S＜0=>Z＜90=>Z=Z+270△△△△←┘S＞0=>Z=Z+90:Z＞360=>Z=Z-360△△←┘X=X+Abs S CosZ◢ Y=Y+Abs SinZ◢ Z=Z◢子程序“H2”V=C-D:W=RJ←┘L=V-V^5/(40W²)+V^9/(3456W^4)-V^13/(599040W^6)←┘N=V^3/(6W)-V^7/(336W^3)+V^11/(42240W^5)-V^15/(9676800W^7)←┘O= tan-1（N/L）：P=√（L²+N²）←┘G＞180=>Q=G-180：≠>Q=G+180△←┘K=1=>M=Q-O：M＜0=>M=M+360△△K=0=>M=Q+O：M＞360=>M=M-360△△←┘X=E+P CosM:Y=F+P SinM←┘T=90V²/(∏W)←┘K=1=>Z=G-T:Z＜0=>Z=Z+360△△K=0=>Z=G+T:Z＞360=>Z=Z-360△△←┘S＜0=>Z＞90=>Z=Z-90: ≠>S＜0=>Z＜90=>Z=Z+270△△△△←┘S＞0=>Z=Z+90△←┘X=X+Abs S Cos Z◢Y=Y+Abs S SinZ ◢Z=Z◢数据厍：名“自起“Lbi 0←┘{DS}:D“ZK”←┘D≤第二缓和曲起点桩号=>A=第一缓和曲起点桩号X坐标:B=第一缓和曲起点桩号Y坐标:C=第一缓和曲线起点:M=第一缓和曲起点方位角:K=路线左转时为0 右转时为1:R=半径:H=缓和曲线长没有缓和曲线时为1E-20:I=圆曲线长:Porg“ZB”:Goto 0:≠>D＞第二缓和曲起点桩号=>D≤第二缓和曲线止点桩号=>E=第二缓和曲线止点桩号X 坐标: F=第二缓和曲线止点桩号Y坐标:C=第二缓和曲止点桩号:G=第二缓和曲止点桩号方位角:K=路线左转时为0 右转时为1:R=半径: J=缓和曲线长: Porg“H2”:Goto 0: ≠>D＞第二缓和曲线止点桩号=>D≤下一交点第二缓和曲起点桩号=>A=第一缓和曲起点桩号X坐标: B=同上………………………………………………………………………………………………………………………………………终点是直线时≠>D＞上一交点第二缓和曲线止点桩号=>D≤路线终点桩号=>A=路线终点桩号X坐标 :B=路线终点桩号X坐标:C= 终点桩号:M=终点桩号方位角:K=路线左转时为0 右转时为1 :Porg“ZB”:Goto 0△←┘。

4800,4850编写程序详解

14
综合上述关于直线，圆曲线，缓和曲线缓和段程序的介绍，道路测量就是这三种线型的组合，计算方法就是这些。希望做的同事在回去工作中多看关于这方面的书籍，了解其原理，能够自己写出更多更好实用性更强的程序。
2020/8/27
15
道路测量直线，缓和曲线。圆曲线综合程序
路线程序（LX） D“JD”：A：﹛R﹜：C“LS” P=C2÷24÷R-C4÷2688÷R3+C6÷506880÷R5-
V=F+ZA：圆曲线上任意点切线方位角
2020/8/27
9
第4步： X=M+Hcos(F+ZI) ▲ 所求任意点X坐标 Y=N+Hsin(F+ZI) ▲ 所求任意点Y坐标第5步： J”J=” 距中桩距离。左-，右+。
X=X+Jcos(V+90) ▲
Y=Y+Jsin(V+90) ▲ 边桩X,Y坐标结果。 Goto 0 同上
经多年的实践经验，自己的总结，编程出了一个通过偏角法进行计算坐标的程序，适用于各种圆曲线，包括匝道的卵形曲线中的圆曲线，希望大家在以后的工作中能够用上。
2020/8/27
7
第1步： M “X0=”:N”Y0=”:R”R=“: F”F=”: 0”0=”:Z”Z=” 同直线赋予值一样，只介绍 R—半径，Z—转角，左-右+
5
第4步
X=X+J*cos(F+90) ▲ Y=Y+J*sin(F+90) ▲
所求的XY坐标就是左右边桩的坐标，90度为正交。
第5步 Goto 0 Goto为无条件转移命令， Goto 0即为无条件转移到
标识符为0的程序运行段落中来

4800及4850计算器程序

4800及4850计算器常用程序ZBZS（坐标正算）E“X1=”:F“Y1=”:G“T1=LbI 0：Fixm：{AS}：G>360=>G=G-360 “T2=”:G DMSE=E+ScosG：“X2=”：F=F+SsinG: “Y2=”X1：起算点的X坐标Y1：起算点的Y坐标T1：起算边指向计算前进方向的坐标方位角A：计算前进方向的左角ZBFS（坐标反算）A“X1=”：B“Y1=”：C“X2=”：D“Y2=”：POL（C-A，D-B）：Fixm：J>0=>W=J:≠“W（1-2）=”“S（1-2）=”：XLYH（线路优化程序）A“XA=”：B“YA=”：C“TA=”：E“1÷RA，Z-，Y+”：P“KA=LbI1: ｛FQ｝:F“1÷RB，Z-，Y+”：Q“KB=LbI2:｛VDZ｝:V“K=”:Z“JIA-JIAO=”:D“WAIYI=(Z-,Y+)”: O=V:O<P=>Goto 3: ≠LbI3: “K-SMALL”LbI4: O>Q=>Goto 5: ≠LbI5: “JISUAN:A-TXY”LbI6: T=C+90+3.1415926×(2E+(F-E)(V-P)÷(Q-P))(V-P):Rad:W=A+∫(cos(C0+0.5×(2E+(F-E)(X-P)÷(Q-P))(X-P), P,V,4)+D×cos(T0+ Z0)Y=B+∫(sin(C0+0.5×(2E+(F-E)(X-P)÷(Q-P))(X-P)),P,V,4)+D×sin(T0+ Z0):LbI 7: O≤Q=>Goto 9:≠LbI 9: W“X=Y“Y=”T“T=”XA：起算点的X坐标YA：起算点的Y坐标TA：起算点的切线方位角1÷RA，Z-，Y+：起算点的曲率（1/R），左偏曲线为正，右偏曲线为负KA：起算点的里程1÷RB，Z-，Y+：终算点的曲率（1/R），左偏曲线为正，右偏曲线为负KB：终算点的里程JIA-JIAO：边桩计算时所要输的与线路切向的夹角，计算点与线路正交时输90度，斜交时则输与线路切向右侧所夹的斜交角度WAIYI=(Z-,Y+)：边桩计算时所要输的计算点到中线点的直线距离（即边桩的外移距离），计算点位于线路左侧时为负，位于线路右侧时为正程序说明：本程序计算线路坐标时要将线路分段（一般分直线段、第一缓和曲线段、圆曲线段、第二缓和曲线段、直线段以此循环）进行计算，例如：缓和曲线加圆曲线计算时，必须将曲线分三段进行计算，先输（起点要素）ZH点的坐标XA、YA，ZH到JD的方位角（即ZH点的切线方位角）TA，ZH 点的曲率1÷RA、里程KA，再输（终点要素）HY点的曲率1÷RB、里程KB，然后输入ZH到HY间任意点的里程便可计算；当所输计算点里程大于HY点里程时，程序显示“JISUAN:A-TXY”表示所要计算点超过第一缓和曲线计算范围，此时程序自动默认HY点的曲线要素为圆曲线计算段的起点要素，用户只需再输入（终点要素）YH点的曲率1÷RB、里程KB，就可计算圆曲线上任意点；当计算点里程大于YH 点里程时，程序显示“JISUAN:A-TXY”表示所要计算点超过圆曲线计算范围，程序自动默认YH点的曲线要素为第二缓和曲线计算段的起点要素，用户只需再输入（终点要素）HZ点的曲率1÷RB、里程KB，就可计算第二缓和曲线上任意点。

4800计算器程序

4800-4850直线、曲线、螺旋线、匝道中边桩坐标计算程序直线、曲线、螺旋线、匝道中边桩坐标计算程序主程序名称（4800-4850）A“XA”：B“YB”：C“CA”：D“RA”：E“RB”：F“DKA”：G“DKB”换行D≠0=>D=1÷D：换行E≠0=>E=1÷E：换行Lb1 1:｛H，O，R｝换行H“DKI”：O“DL”：R“DR” 换行Lb1 0 换行P=（E－D）÷Abs（G－F）换行Abs（H－F）＞Abs（G－F）=>Q= Abs（G－F）：≠=>Q= Abs（H－F）：换行I=P×Q 换行J=C＋（I＋2D）×Q×90÷π 换行M=C＋（I÷8＋2D）×Q×45÷（4π）换行N=C＋（3I÷8＋2D）×Q×135÷（4π）换行U=C＋（5I÷8＋2D）×Q×225÷（4π）换行V=C＋（7I÷8＋2D）×Q×315÷（4π）换行K=C＋（I÷4＋2D）×Q×45÷（2π）换行W=C＋（I÷2＋2D）×Q×45÷π 换行Z=C＋（3I÷4＋2D）×Q×135÷（2π）换行X=A＋Q÷24×（cosC＋4×（cosM＋cosN＋cosU＋cosV）＋2×（cosK＋cosW＋cosZ）＋cosJ）换行Y=B＋Q÷24×（sinC＋4×（sinM＋sinN＋sinU＋sinV）＋2×（sinK＋sinW＋sinZ）＋sinJ）换行U=X＋O×cos（J－90）换行V=Y＋O×sin（J－90）换行W=X＋R×cos（J＋90）换行Z=Y＋R×sin（J＋90）换行Abs（H－F）＞Abs（G－F）=> Goto 3：换行H“DKI=” ◢ fx4850输入改为“DKI=”：H ◢X“X=” ◢ fx4850输入改为“X=”：X ◢Y“Y=” ◢ fx4850输入改为“Y=”：Y ◢U“XL=” ◢fx4850输入改为“XL=”：U◢V“YL=” ◢fx4850输入改为“YL=”：V◢W“XR=” ◢fx4850输入改为“XR=”：W◢Z“YR=” ◢fx4850输入改为“YR=”：Z◢Goto 1 换行Lb1 3 换行A=X：B=Y：D=E：F=G：C=J 换行｛E，G｝换行E“RB”：G“DKB”换行E≠0=>E=1÷E：换行Goto 0计算器显示及输入方法：XA？线型起算点的X坐标。

CASIO(4850-4800)交点法计算平曲线要数各线元起点坐标及切线方位角(更新)1

平曲线要数、线元各起点坐标及切线方位角计算程序CASIO（4850-4800）原创版单位：中铁二十局集团第二工程有限公司编程：辜建君Defm12:LbI0:{ABCFΟZREK}:A”JD”:B”JDX”:C”JDY”:F”FJ”:Ο”ZJ:Z-,Y+”:Z:R:E”LS1”:K”LS2”:E<1=>E=1E-9⊿K<1=>K=1E-9⊿Z[2]=E÷2-EEE÷240RR:Z[6]=K÷2-KKK÷240RR: Z[7]= E-E^3÷40R^2+E^5÷3456R^4-E^7÷599040R^6+E^9÷17542600R^8:Z[8]= E^2÷6R-E^4÷336R^3+E^6÷42240R^5-E^8÷9676800R^7+E^10÷3530096640R^9: Z[9]= K-K^3÷40R^2+K^5÷3456R^4-K^7÷599040R^6+K^9÷17542600R^8:Z[10]= K^2÷6R-K^4÷336R^3+K^6÷42240R^5-K^8÷9676800R^7+K^10÷3530096640R^9: G=sin-1((Z[7]-Z[2])÷R): U= sin-1((Z[9]-Z[6])÷R): Z[1]=RcosG+Z[8]-R: Z[5]=RcosU+Z[10]-R ：V=(EE-KK)÷24R÷SinAbsΟ: ”T1”:Z[3]=(R+Z[1])tan(AbsΟ÷2)+Z[2]-V◢”T2”:Z[4]=(R+Z[5])tan(AbsΟ÷2)+Z[6]+V◢I=tan-1((R+Z[1])÷(Z[3]-Z[2])):“L”:L=AbsΟлR÷180+(E+K)÷2◢”E0”:U=(R+Z[1])÷sinI-R◢X=A-Z[3]:Y=X+E:E<1=>GOTO 1: ≠> GOTO 2 ⊿LbI 1:”ZY”:X◢“ZYX=”:N=B+Z[3]Cos(F+180) ◢“ZYY=”:Q=C+Z[3]sin(F+180) ◢“FWJ=”:F->DMS◢GOTO 3LbI 2:”ZH”:X◢“ZHX=”:N=B+Z[3]cos(F+180) ◢“ZHY=”:Q=C+Z[3]sin(F+180) ◢“FWJ=”:F->DMS◢“HY”:Y◢”HYX=”:N=B+(Z[3]-Z[7])cos(F+180)+Z[8]cos(F+90Z) ◢“HYY=”:Q=C+(Z[3]-Z[7])sin(F+180)+Z[8]sin(F+90Z) ◢J=F+ZG：J>360=>J=F+ZG-360:⊿J<0=>J=F+ZG+360: ⊿“FWJ=”:J->DMS◢LbI 3:”QZ”:S=X+(L-K-E)÷2+E◢D=Rsin(90 * (L-K-E)÷ЛR+sin-1((Z[7]-Z[2])÷R))+Z[2]:H=R(1-cos(*90 (L-K-E)÷ЛR+sin-1((Z[7]-Z[2])÷R)))+Z[1]“QZX=”:N= B+(Z[3]-D)cos(F+180)+Hcos(F+90Z) ◢“QZY=”:Q= C+(Z[3]-D)sin(F+180)+Hsin(F+90Z) ◢J=F+ZG+90Z(L-K-E)÷ЛR:J>360=>J= F+ZG+90Z(L-K-E)÷ЛR -360: ⊿J<0=>J=F+ZG+90Z(L-K-E)÷ЛR+360:⊿“FWJ=”:J >DMS◢LbI 4:M=X+L-K:P=X+L:K<1=>GOTO 5: ≠> GOTO 6⊿LbI 5:”YZ”:P◢“YZX=”:N=B+Z[4]com(F+Ο) ◢“YZY=”:Q=C+Z[4]sin(F+Ο) ◢W=F+Ο:W>360=>W=F+Ο-360: ⊿W<0=>W=F+Ο+360: ⊿”FWJ=”:W->DMS◢G0TO 0LbI 6:”YH”:M◢”YHX=”:N=B+(Z[4]-Z[9])cos(F+Ο)+Z[10]cos(F+Ο+90Z) ◢“YHY=”:Q=C+(Z[4]-Z[9])sin(F+Ο)+Z[10]sin(F+Ο+90Z) ◢J=F+ZG+180Z(L-K-E)÷ЛR:J>360=>J= F+ZG+180Z(L-K-E)÷ЛR -360: ⊿J<0=>J=F+ZG+180Z(L-K-E)÷ЛR+360:⊿“FWJ=”:J >DMS◢“HZ”:P◢“HZX=”:N=B+Z[4]com(F+Ο) ◢“HZY=”:Q=C+Z[4]sin(F+Ο) ◢W=F+Ο:W>360=>W=F+Ο-360: ⊿W<0=>W=F+Ο+360: ⊿”FWJ=”:W->DMS◢GOTO 0说明：一、便于区分数字零于字母O，在程序中用Ο代替字母O二、本程序适合CASIO4850、4800、4500计算器使用。

4800--4850隧道断面程序坐标正反算

4800--4850隧道断面程序坐标正反算1.主程序(TYQXJS)"1.SZ => XY"："2.XY => SZ"：N：Prog"SJK":C=1÷P：D=(P-R)÷(2HPR)：E=180÷π：N=1=>Goto 1：≠>Goto 2Δ←┘Lbl 1：{SZ}：SZ：W=Abs(S-O)：Prog "SUB1"：X"XS"=X◢Y"YS"=Y◢F"FVJ"=F-90:F◢Goto 1←┘Lbl 2：{XYH}：XY：I=X：J=Y：Prog "SUB2"：S"S"=O+W◢Z"Z"=Z◢A=407.514里程+0.0053坡度(S-136694.623里程)-(S-136694.623里程)^2/3200+高差:←┘Z[1]=A+√((半径6.49^2-(Z-中线5.901)^2)◢H=√((Z-5.901)^2+(H-A)-6.49◢Goto 22. 正算子程序(SUB1)A=0.1739274226：B=0.3260725774：K=0.0694318442：L=0.3300094782：F=1-L：M=1-K：X=U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Bcos(G+QEFW (C+FWD))+Acos(G+QEMW(C+MWD)))：Y=V+W(Asin(G+QEKW(C+KWD))+Bsin(G+ QELW(C+LWD))+Bsin(G+QEFW(C+FWD))+Asin(G+QEMW(C+MWD)))：F=G+QEW(C+ WD)+90：X=X+ZcosF：Y=Y+ZsinF3. 反算子程序(SUB2)T=G-90：W=Abs((Y-V)cosT-(X-U)sinT)：Z=0：Lbl 0：Prog "SUB1"：L=T+QEW(C+ WD)：Z=(J-Y)cosL-(I-X)sinL：AbsZ<1E-6=>Goto1：≠>W=W+Z：Goto 0Δ←┘Lbl 1：Z=0：Prog "SUB1"：Z=(J-Y)÷sinF本程序"SUB1"和"SUB2"同上,新加"SJK"子程序为曲线要素数据库,也修正了主程序里的几个问题,H为断面高程计算抄欠挖.本程序"SJK"数据库为:S≤136750.876(HY)＝> O=136620.876(ZH):U=90246.898(136620.876的X坐标)：V=327.519(136620.876的Y坐标):G=198°12′37.8″：H=130：P=1960：R=1960：Q=-1:⊿←┘S≤****(YH)＝> O=****(HY):U=*******(HY的X坐标)：V=******(HY的Y坐标):G=***°**′**″：H=***：P=***：R=***：Q=-1:⊿←┘S≤****(HZ)＝> O=****(YH):U=*******(YH的X坐标)：V=******(YH的Y坐标):G=***°**′**″：H=***：P=***：R=***：Q=-1:⊿←┘依次往下推,本人昨天晚上才研究出来的,希望给以参考,并想给与高手修正,因本人的能力有限..这个程序只适合反算,正算的时候会出现很多问题,不适合正算,想给与高手指点,希望把正算的也给休整好.谢谢,本人QQ:8696847,希望给予联系,共同克服里面的程序的困难.2、输入与显示说明输入部分：1. SZ => XY2. XY = > SZN ? 选择计算方式，输入1表示进行由里程、边距计算坐标；输入2表示由坐标反算里程和边距。

卡西欧4800 4850计算器程序

卡西欧4800\4850计算器计算程序及说明一：坐标反算：(ZBFS)1 程序步骤：C”X1”:D”Y1”:E”X2”:F”Y2”:Fixm:Pol(E-C,F-D:I”S1-2=”▲J≤Ｏ＝＞J=J+360△J”A1-2=”2 操作过程：ZBFS→EXE→输入X1值→EXE→输入Y1值→EXE→输入X2值→EXE→输入Y2值→EXE→EXE→EXE注：①：S1-2:计算得出的距离。

②：A1-2:计算得出的角度。

（按shift°′″转换为60进制的角度）③：此程序可循环计算。

3 案例：已知：D1坐标（7811.23，606.136），D2坐标（7805.915，431.910）求解： a D1-D2(A1-2)=268°15′09.56 ″S D1-D2(S1-2)=174.3071二：坐标正算：(ZBZS)1 程序步骤：C”X1”:D”Y1”:S”S1-2”:R”A1-2”:Fixm:X”X2”=C+Rec(S,R) ▲Y”Y2”=D+J2 操作过程：ZBZS→EXE→输入X1值→EXE→输入Y1值→EXE→输入S1-2距离值→EXE→输入A1-2角度（例268°15′09.56 ″）值→EXE→EXE注：①：X2计算得出坐标，Y2计算得出坐标。

②：此程序可循环计算。

3 案例：已知：D1坐标（7811.23，606.136），两点距离：S D1-D2(S1-2)=174.3071方位角：a D1-D2(A1-2)=268°15′09.56 ″ 求解：D2坐标（7805.915，431.910三：线路曲线计算程序:（XLQXJS）(万能公式)1 程序步骤：1.Lbl 0:{E G}：A“XA”：B“YA”：C“CA”：D“1/RA”：E“1/RB”：F“DKA”： G“DKA”2.Lbl 1:{H O R}:H“DKI”：O“DL”：R“DR”：H＞G=﹥Goto 3⊿3.P=(E-D)/Abs(G-F):Q=Abs(H-F):I=P×Q:T=D+I4.J=C+(I+2D)×Q×90/π ▲5.M=C+(I/4+2D)Q×45/(2π):N=C+(3I/4+2D)Q×135/(2π):6.K=C+(I/2+2D)Q×45/π7.X=A+Q(cosC+4(cosM+cosN)+2cosK+cosJ)/12 ▲8.Y=B+Q(sinC+4(sinM+sinN)+2sinK+sinJ)/12 ▲9.U“XL”=X+Ocos(J-90) ▲10.V“YL”=Y+Osin(J-90) ▲11.W“XR”=X+Rcos(J+90) ▲12.Z“YR”=Y+ Rsin(J+90) ▲13. A=X：B=Y：C=J：D=T：F=H：Goto 1Lbl 3: A=X：B=Y：D=E：F=G：C=J：Goto 0程序说明：A- 曲线元起点A的坐标；B- 曲线元起点B的坐标；C- 曲线元起点A的切线坐标方位角；F- 曲线元起点A的里程；G- 曲线元起点B的里程；H- 曲线上待求点i的里程；D- 曲线元起点A的曲率；E- 曲线元终点B的曲率；XL-左边线点位X坐标；YL-左边线点位Y坐标；XR-右边线点位X坐标；YR-右边线点位Y坐标；X- 中线点位纵坐标；Y- 中线点位横坐标；DL-左边线距中线平距；DR-右边线距中线平距；该程序需要输入的数据为：（1）.曲线元起点A的坐标及切线坐标方位角，计算器上用“XA”，“YA”，“CA”显示；（2）.曲线元起点A和B的曲率，计算器上用I÷RA，I÷RB 显示（曲线左偏时取“-”）；（3）.曲线元起点A和终点B的里程，计算器上用“DKA”，“DKB”显示；（4）.输入待求点里程和该点距左右的水平距离，计算器上用“DKI”，“DL”，“DR”显示；每算完一个待求点的中线及边线坐标，程序又让输入下一点的“DKI”，“DL”，“DR”当输入的“DKI”大于“DKB”时程序中显示“1÷RB”和“DKB”，此时输入下一个曲线元起点的曲率和里程，然后重复步骤（4），即可计算下一个曲线中线及边线点位坐标。

4850计算器坐标正反计算程序

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

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

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

2．计算原理：利用Gauss-Legendre 5点通用公式计算线路中边桩坐标并计算放样数据。

二、源程序1、程序(XYZ)Defm 2:F=1：Z[1]=90Z“KM“:A=第一段线终点桩号：B=第二段线终点桩号：C=第三段线终点桩号：E=第四段线终点桩号：F=第五段线终点桩号：……（有多少段就加多少段;超过十五段，要另加子程序转过来）Z<A=>Prog”C1”:≠>Z<B=>Prog”C2”: ≠>Z<C=>Prog”C3”: ≠>Z<E=>Prog”C4”: ≠> Z<F=>Prog”C5”ΔΔΔΔΔ……（有几段线，则打几个三角）C=1÷P：S=(P-R)÷(2HPR)：E=180÷π：F=1=>Goto 1：≠>Goto 2Δ←Lbl 1：D：W=Z-O：Prog "A"：X"XZ"=X◢Y"YZ=":Y=Y◢Lbl 2：{XY}：XY：I=X：J=Y：Prog "B"：Z"KM"=O+W◢D=D◢2. 正算子程序(A)A=0.1184634425:B=0.2393143352:N=0.2844444444:K=0.0469100770:L=0.2307653449:M =0.5:X=U+W(Acos(G+QEKW(C+KWS))+Bcos(G+QELW(C+LWS))+Ncos(G+QEMW(C+MWS))+Bcos( G+QE(1-L)W(C+(1-L)WS))+Acos(G+QE(1-K)W(C+(1-K)WS))): Y=V+W(Asin(G+QEKW(C+KWS ))+Bsin(G+QELW(C+LWS))+Nsin(G+QEMW(C+MWS))+Bsin(G+QE(1-L)W(C+(1-L)WS))+Asin( G+QE(1-K)W(C+(1-K)WS)))：Z[2]=G+QEW(C+WS)+Z[1]：X=X+Dcos Z[2]：Y=Y+Dsin Z[2]3. 反算子程序(B)T=G-90：W=Abs((Y-V)cosT-(X-U)sinT)：D=0：Lbl 0：Prog "A"：L=T+QEW(C+WS)：Z=(J -Y)cosL-(I-X)sinL：AbsZ<1E-6=>Goto1：≠>W=W+Z：Goto 0Δ←┘Lbl 1：Z=0：Prog "A"：Z=(J-Y)÷sinZ[2]4.曲线元要素数据库：C1~C5…….O=xxxx.xxx:U=xxxx.xxx:V=xxxx.xxx:G=xxxx.xxx:H=xxxx.xxx:P=xxxx.xxx:R=xxxx.xxx :Q=x (注：一段线一个子程序，分别提供这几个要素数据。