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

CASIO fx4850线路正反算测量程序（全站仪伴侣）CL-测量Norm :Lbl 1:Cls:{U}:U"1 ZS 2 FS"=1 =>Prog "ZS":≠>U=2 =>Prog "FS":≠>Goto 1△△D-读Lbl 1:K[B]=Z[G+B]:Dsz B:Goto 1 FS-反算Lbl 1:Cls:{KVW}:V"X"W"Y":Lbl 2:Prog "P":Pol(V-X,W-Y):J=J-A:I=Rec(I,J):K=K+I:AbsI≥E-3=>Goto 2:≠>Prog "WZ":Goto 1△HQ-缓曲N=I-I^5÷40S2+I^9÷3456S^4-I^13÷599040S^6:M=I^3÷6S-I^7÷336S^3+I^11÷42240S^5:Z=90IAbsI÷πS:J=90-(90-tan-1(M÷N))AbsI÷I:I=√(M2+N2)HY-缓圆I=O:Prog "HQ":Prog "XY":I=K-L-O HZ-缓直I=P: Prog "Y": Prog "XY":I=Q:S=RQ:Prog “HQ”:J=Z-J:Prog "XY":I=K-L-O-P-QML-目录Cls:Fix 3:C=Z[100]:D=Z[101] :Lbl 1:B=6:G=C-6D:Prog "D":K≥L =>K≤M =>C=N:D=O:E=P:F=Q: Goto 2△△Dsz D:Goto 1:Lbl 2P-平Prog"ML": Prog"PM"PM-平面Lbl 1:B=9:G=C-8D-1:Prog "D":K≤T =>Prog "PQ":Goto 2△Dsz D:Goto 1:Lbl 2 PQ-平曲Prog "ZH":K≤L =>Prog "ZX":≠>K≤L+O =>Prog "HQ":≠>Prog "HY":K≤L+O+P =>Prog "Y":≠> Prog "HZ":K≤L+O+P+Q => Prog "HQ" :≠>Prog "ZX"△△△△Prog "XY"PY-偏移Lbl 1:I=0:J=90:{IJ}:I"L"≠0 =>Z=J"<":Prog "XY":Prog "ZB":Goto 1△PZ-平纵Prog "ML":Prog "PM":Prog "ZM"SQ-竖曲Z=N+(K-M)L:K≤M-T =>Z:≠>Z=Z- (K-M+T)2÷2O△WZ-位置"KJ":K:Pause1:J▲XY-XYJ=A+J: A=A+Z:X=X+Rec(I,J):Y=Y+JY-圆J=90I÷πR: I=2RSinJ: Z=2JZ-纵Prog"ML": Prog"ZM"ZB-坐标A=360Frac((A+360)÷360):"ZB":X:Pause 1:Y:Pause 1ZH-直缓X=M:Y=N:A=S:I=K-L:S=ORZM-纵面Lbl 1:B=5:G=E-4F-1:Prog "D":T=Abs(O(P-L)÷2):K≤M+T =>Prog "SQ":Goto2△Dsz F:Goto 1:Lbl 2ZS-正算{U}:Lbl 1:{K}:U"1 B"=1=> Prog "P": Prog "PY" :Goto 1△Prog "PZ":Prog "ZB": Z▲Goto 1ZX-直线J=0:Z=0使用说明运行CL程序首先程序提示正反算1 ZS2 FS？输1正算、输2反算，其它重选。

一．文件名:GSZFS（高斯-勒让德）Deg:I=0:I“JS.1,FS.2,SJK.3”：I=1=>GotoA⊿I=2=>GotoC⊿I=3=>GotoB⊿↙LbI A:{KZ}:QXYFLPROKZ“―0+”:Prog“A”:G=G-90:G<0=>G=G+360:G>36 0=>G=G-360⊿⊿“XYF=”:U:Pause0:V: Pause0:G→DMS▲(Q+J)=K=>Goto0: ≠> GotoA⊿↙LbI 0:Q=K:X=U:Y=V:F=G:P=R: “-----AC-----”: ▲ Prog“GS” ↙LbI B: {KZ}:K:Z“―0+”: Prog“SJK”: Prog“A”:C“C.X”:T“C.Y”:I=0:J=0:PoI(C-U.T-V: “FY=”:J=J+180→DMS I▲GotoB↙LbI C:{MN}:M“(X)”:N“(Y)”:K=0: C=0：I=0:J=0:PoI(M-X,N-Y:K=Q+Icos(F-J) ↙LbI 2:K=K-C: Prog“A”:PoI(U-M,V-N:T=G-90-J:C=IcosT:Int(1000C)=0=>Goto3:≠> Goto2⊿↙LbI 3:sinT<0=>A= -1:≠>A=1⊿“(K)=”:K+C: Pause0:“(―0+)=”:IA▲ GotoC↙二．文件名:A（内核）W=K-Q:D=OrW2(P-R)/2LPR:E=OrWP-1:A=.1184634425:B=.2393143352 :S=.046910077:H=.2307653449:Z[1]=F+SE+S2D:Z[2]=F+HE+H2D:Z[3]= F+.5E+.25D:Z[4]=F+(1-H)2D+E(1-H):Z[5]=F+(1-S)2D+E(1-S):G=F+E+D +90+(斜交) ↙U=X+ZcosG+W(AcosZ[1]+BcosZ[2]+.2844444444cosZ[3]+BcosZ[4]+Aco sZ[5]):V=Y+ZsinG+W(AsinZ[1]+BsinZ[2]+.2844444444sinZ[3]+BsinZ[4]+AsinZ[5])三.文件名:SJK（数据库）K≤下个起点桩号=>Q=起点桩号:X=起点（X）:Y=起点（Y）:F=前切线方位角:L=线型长度:P=起点（R）:R=终点（R）:O=转向:≠>K≤。

CASIO 4850 道路测量正反算程序变量意义（凡是坐标均按复数a+bi输入，所有的值均按数学坐标系考虑）以下所说的线元为直线、圆曲线、缓和曲线A 每个线元的起点坐标（由设计图提供，或者计算出来）B 直线线元的终点（方向）坐标（由设计图提供，或者计算出来）C 正反算时计算点的坐标结果（中间结果，不做事前计算）D 反算时计算点的坐标（现场实测值,手动输入）E 反算时计算桩号的向量（前一点）（中间结果，不做事前计算）,显示结果时为正东坐标。

F 反算时计算桩号的向量（后一点）（中间结果，不做事前计算）G 反算时计算桩号的步长（程序设定为20米,可自定义设定，一般10--50）H 反算时计算桩号前一点的坐标结果（中间结果，不做事前计算）I (D-C)与E的夹角（中间结果，不做事前计算）J (D-H)与F的夹角（中间结果，不做事前计算）K 缓和曲线的曲线总长（由设计图提供，或者计算出来）L 正算时为偏移计算点的长度，反算时为0（正算时手动输入，左正右负）M 正反算判断，等于0时进入正算，等于1时进入反算N 反算时起点桩号，显示结果时为正北坐标。

O 圆曲线的圆心坐标（由设计图提供，或者计算出来）P 曲线左偏时为-1，右偏时为+1（根据设计图事先确定）Q 缓和曲线起点的切线角（由设计图提供，或者计算出来）R 圆曲线的半径，缓和曲线的终点半径（由设计图提供，或者计算出来）S 计算点桩号（正算时手动输入）T 每个线元的起点桩号（由设计图提供，或者计算出来）U 反算时终点桩号V 缓和曲线在标准位置时的坐标计算结果（中间结果，不做事前计算）W 与道路中线的偏角（手动输入，左偏转范围0度--180度，没有右偏）XYZ 缓和曲线在标准位置时的切线角计算结果（中间结果，不做事前计算）主程序清单(每个语句后均有回车键)说明：→条件判断△条件判断结束▲显示结果≠→条件判断的否则Lb1 0M=0:F=0:H=0:G=20:N=?:U=? 设定中间参数(N和U是起点和终点的桩号)Lb1 1"ZHENG SUAN=0""FAN SUAN=1":{M} 根据变量M的值确定正反算M=0→Goto 2△M=0,进入正算程序M=1→Goto 3△M=1,进入反算程序Goto 0 M=其他，回到起点Lb1 2Prog "计算程序":Goto 0 进入正算程序，Lb1 3 进入反算程序S=N 起点桩号是整个平曲线的起点，此时输入程序中，下同（可自定义起点）“JI SUAN DIAN”:{D} 输入计算点坐标（a+bi）S=S-GLb1 5S=S+GS>U→“GUO ZHONG DIAN”▲GOTO 0△终点桩号是整个平曲线的终点，下同.Pause 4,暂停2秒。

CASIO4850计算器公路测量计算程序程序设计:魏加训2009.3.10Casio 4850计算器数据库型全线任意位置高程计算程序一、主程序：HLbl 0: X “BM-H”：U “HS”：“SXG=”: U+X◢输入水准点高程及后视，显示视线高程Lbl 1：{ P K }: “QS”K “K××+×××”: Prog“SJK”↙（换路修改此处为对应数据库名即可）C= C÷100: D= D÷100:↙D>C=> W=1:≠>W = -1⊿K<A=> I= -C:≠> I=D⊿↙Abs (A-K):L> R Abs(D-C)÷2 =>H=Z+IL: "H(Z-SJ)=": H◢Goto2↙≠>H= Z+IL+W(R Abs(D-C)÷2- Abs (A-K))2÷2÷R: "H(Z-SJ)=": H◢Lbl 2: {ONJP}:O “JGC(h)” : N“Z”:J “HP”:N=Abs N : B=H-O+N J÷100: "H(C) =" : U+X-P :Pause 0 : “H(S)=”: B◢“⊿h=”: B-U-X+P◢输入前视后显示实测高、设计高及高差（+填，-挖）Goto 1二、竖曲线数据库格式：K<下一竖曲线起点里程=>C=本竖曲线变坡点之前坡度(%前分子有正负): D=本竖曲线变坡点之后坡度(%前分子有正负): R=本竖曲线半径: A=本竖曲线变坡点桩号: Z=本竖曲线变坡点高程: Goto 1⊿↙三、数据库样例: SJKK<下一竖曲线起点里程=> C=本竖曲线变坡点之前坡度(%前分子有正负):D=本竖曲线变坡点之后坡度(%前分子有正负):R=本竖曲线半径: A=本竖曲线变坡点桩号: Z=本竖曲线变坡点高程: Goto 1⊿↙K<下一竖曲线起点里程=> C=本竖曲线变坡点之前坡度(%前分子有正负):D=本竖曲线变坡点之后坡度(%前分子有正负):R=本竖曲线半径: A=本竖曲线变坡点桩号: Z=本竖曲线变坡点高程: Goto 1⊿↙…………………………………………………………………………………………………………………..(注：如一条路有多个竖曲线要素则按顺序排列依次添加入数据库S.Z中)K≤下一竖曲线起点里程=> C=本竖曲线变坡点之前坡度(%前分子有正负): D=本竖曲线变坡点之后坡度(%前分子有正负): R=本竖曲线半径: A=本竖曲线变坡点桩号: Z=本竖曲线变坡点高程: Goto 1⊿↙Lbl 1↙四、输入与显示说明BM-H ？输入水准点的高程HS？输入水准点上的后视读数SXG= ×××显示视线高程“QS”？输入所求点前视读数K××+×××? 输入所求的桩号H(Z-SJ)= ×××显示所求桩号的路面中桩设计高程JGC(h) ? 输入结构层厚度或从设计路面高程点处垂直向下的厚度Z ？输入所求点距中线的边距(注:不分正负号均可，有中央分隔带的输入时需自行扣除中央分隔带半宽)HP？输入所求桩号处的横坡(%前分子有正负).H(C) = 显示计算得出的所求点的实测高程H(S) = 显示计算得出的所求点的设计高程⊿h= 显示设所求点的计与实测的高差（+填，-挖）注：（①当Z和JGC(h)为0时，则H=×××为路面中桩设计高程等于前面显示的H(S) ；②当Z为0，JGC(h)不为0时，则H=×××为路面垂直向下至某层的中桩设计高程；③当Z、JGC(h)、HP三者同时不为0时，则H=×××为路面垂直向下至某层的边桩对应Z位置的设计高程；④当Z不为0、JGC(h)为0、HP不为0时，则H=×××为路面边桩对应Z位置的设计高程；）附：数据库中用到的竖曲线基本计算公式：竖曲线起点桩号=本竖曲线变坡点桩号+本竖曲线之切线长T竖曲线切线长T=RAbs（i2-i1）÷2。

程序一、(YX)道坐标正反算引导程序，主程序LBi 1:N“1--2”←N=1=> Prog " YXSJ（线路平面数据库引导程序)": Prog"BZZB边桩坐标"： Prog “HYXSJ （ZX高程数据库”: Prog"H"（高程计算程序）←“X0=”：X◢“Y0=”：Y◢“H0=”：H◢（显示计算点的设计高程）≠> Lbl 1 Goto 1:≠＝＞K=3000:LBi 4：Prog " YX（线路平面数据库引导程序)":Prog " ZZZB ":U “X1”: V “Y1”:Poi(U-X,V-Y)(A点至B点坐标反算）： A=tan-1((V-Y)/(U-X)):(U-X)≤0＝＞A=A+1800：J≤0=>J=J+360∠（B点至A的方位角）≠＝＞A=A+3600△R=A-T：R=J-T：（计算角CAB） I=ㄏ（（V-Y）2+（U-X）2）：F=I(i为AB距离）*CosR：（计算AC距离）K=K+F：Abs（绝对值）F≥0.001=>Goto 1:(设置反算结果的精度，达不到就返回再运算) ≠>Z=I*SinR△“K=”：K◢“Z=”：Z◢（说明：输入1时进入线路坐标正算,输入2时进入坐标反算,当进入坐标反算时会先提示输入里程，这是因为这个坐标反算程序需要先输入一个离计算点相近的里程作为引导参考计算，逐渐接近，可以设置精度，F≥0.001，即计算结果≤0.001m停止计算，并显示结果。

）程序二、(BZZB)边桩坐标计算Prog "ZZZB"：X=X+ Cos(T+G)*Z(G为计算边桩与中线法切线方向的夹角，正为右偏，负为左偏，Z为中桩到边桩的距离，G，Z为0时计算中桩坐标)：Y=Y+Sin（T+G）*Z：“R0=”：T▲程序三、(ZZZB)核心计算程序，线路中桩坐标通用计算公式S=绝对值Abs（K-A）:(计算点到线元起点距离恒为正值)P= (D-C)/绝对值（Abs(B-A)) ：(计算弧长恒为正值)T=R+E( (C+P*S)+ C)*S/2*180/∏:(计算点之切线方位角)T＞360=>T=T-360△T＜0=>T=T+360△A=4 (Cos(R+E((C+P* S/4)+ C)* S/23*180/∏)+ Cos(R+ E ( (C+P* (S/4+ S/2))+ C)* (S/4+ S/2)/2*180/∏))：B=4 (Sin(R+E( (C+P* S/4)+ C)* S/23*180/∏)+ Sin (R+ E ((C+P* (S/4+ S/2))+ C)* (S/4+ S/2)/2*180/∏))←X=X+S/2/6*｛A+2(Cos(R+E((C+P* S/2)+ C)* S/4*180/∏)) +CosT + CosR｝←Y=Y+S/2/6*｛B+2(Sin (R+E((C+P* S/2)+ C)* S/4*180/∏)) +Sin T + Sin R｝：为了尽量简化，此方程式N取值为2等分，如果用在卵形曲线上小于50的半径精度不够。

这是本人自编的Casio-4850测量程序,实用!请测友提出保贵意见! 愿和测友共享!符号说明-------------------A ？交点桩号B ？转角(左-右+)C ？缓和曲线长D = Li (距ZH点长度)E ？X [起算点]F ？Y [起算点]G ？X [交点]H ？Y [交点]I = POL（）：r极距J = POL（）：θ极角K ？DK [桩号］L = 曲线长（ZH～HZ）M ？X[后视点]N ？Y[后视点]O= 方位角[ZH～JD]P = 方位角[JD～HZ]Q = Li(距HZ点长度）R ？半经S = HD(水平距离)T = T (切线长)U ？X[测站]V ？Y[测站]W = 测站至后视方位角X = X坐标Y = Y坐标Z = 测站至测点方位角Z[1]= β0Z[2]= mZ[3]= pZ[4]= [ZH]点X坐标Z[5]= [ZH]点Y坐标Z[6]= [HZ]点X坐标Z[7]= [HZ]点Y坐标Z[8]= 圆曲线所夹圆心角Z[9]= 待测点X坐标Z[10]= 待测点Y坐标=============================================================== ================程序-------------------------------------------------------------------------------DEG: Defm10：A“DK[JD]”：B“JD°°°Z[-]Y[+]”→DMS： R：C“L0”：E“X[QSD]”：F“Y[QSD]”：G“X[JD]”：H“Y[JD]”：U“X[CE ZHAN]”：V“Y[CE ZHAN]”：M“X[HOU SHI]”：N“Y[HOU SHI]”： Z[1]=90C/(πR)：Z[2]=C/2-C^3/(240R^2): Z[3]=C^2/24R-C^4/(2688R^3):T=（R+Z[3]）tg（AbsB/2）+Z[2]：L=πR（AbsB-2Z[1]）/180+2C：I=0：J=0： Pol（M-U，N-V）：J<0=>J= J+360：△：“FWJ[CZ～HS]=”：J→DMS▲LbI A｛K｝：K“DK”：I=0：J=0：Pol（G-E，H-F）：J＜0=>J=J+360：△：P=J+B： P＜0=>P=P+360：△：Z[4]=G+Tcos（J+180）：Z[5]=H+Tsin（J+180）：Z[6]=G+TcosP：Z[7]=H+Tsi nP：D=K-A+T：D＜0=>X=﹣D: Y=0：O=J+180：Goto 1：△D＜C =>X=D- D^5/(40R^2C^2) + D^9/(3456R^4C^4)-D^13/(599040R^6C^6): Y=D^3/(6RC) -D^7/(336R^3C^3) +D^11/(42240R^5C^5):B＜0=>Y=-Y:△: O=J: Goto1:△:D≤L-C=>Z[8]=180(D-0.5C)/(πR): X=RsinZ[8]+Z[2]: Y=R-RcosZ[8]+Z[3]: B＜0=>Y=-Y:△: O=J: Goto1: △:Q=K-A+T-L: Q＜0 => Q=-Q：X=Q- Q^5/(40R^2C^2) + Q^9/(3456R^4C^4)-Q^13/(599040R^6C^6):Y=Q^3/(6RC) -Q^7/(336R^3C^3) +Q^11/(42240R^5C^5):B＞0=> Y=–Y：△：P=P+180：Goto2：△：Q≥0=>X=Q：Y=0：P=P： Goto2：△LbI 1 :Z[9]=Z[4]+XcosO-YsinO: Z[10]=Z[5]+XsinO+YcosO: goto 3LbI 2 :Z[9]=Z[6]+XcosP-YsinP: Z[10]=Z[7]+XsinP+YcosP: Goto 3LbI 3 :I=0: J=0: Pol(Z[9]-U,Z[10]-V): J<0=>J=J+360:△“X=”:Z[9]▲“Y=”:Z[10]▲“FWJ[CZ～DK]°°°=”: J→DMS“HD=”:I▲Goto A=============================================================== =================注:①本程序用于CASIO fx-4850p型计算器。

.CASIO fx-4850p公路工程计算器程序集编号主程序名（其中文含义）程序计算功能一FANG YANG （放样）计算路线中、边桩坐标，测站点与路线中、边桩坐标间的距离和方位角等二CHAO GAO (超高) 计算设计高程、左右两侧路基横坡度、宽度、左中右处设计标高之差值及设计高程三YYLDFWJ (任意两点方位角) 计算任意两点间的距离和方位角四SQXBG（竖曲线标高）计算竖曲线设计标高五TIAN WA G-C(填挖高差) 计算左、中、右处填挖高差六 BZZB （边桩坐标）主要进行桥梁及涵洞各点放样计算七 ZHUI PO （锥坡）进行锥坡放样计算编号子程序名（其中文含义）程序计算功能1 XY （坐标计算子程序）利用交点参数数据库进行计算2 P-Q-X J-S （平曲线计算子程序）利用交点坐标进行计算曲线要素3 YAO-SU (平曲线要素计算子程序) 进行平曲线要素计算4 ZB (平曲线的坐标计算子程序) 进行平曲线的坐标计算5 X1 （曲线要素数据库子程序）曲线要素数据库（按桩号先后顺序依次输入）6 ZHONG ZHUANG（中桩坐标计算之子程序）计算任意点中桩坐标及其切线方位角7 ZH HY（直缓点缓圆点处中桩坐标计算之子程序）计算ZH至HY间中桩坐标及其切线方位角8 HY YH（缓圆点圆缓点坐标计算之子程序）计算HY至YH间中桩坐标及其切线方位角9 YH HZ（圆缓点直缓点坐标计算之子程序）计算YH至HZ间中桩坐标及其切线方位角10 S0（缓和曲线长度等于零坐标计算之子程序）计算S=0的圆曲线中桩坐标及其切线方位角11 ZXZB（直线段坐标计算之子程序）计算直线段中桩坐标及其切线方位角12 CHAO GAO-1（超高计算之子程序）计算左右两侧路基横坡度、宽度、左中右处设计标高之差值13 JK--B[2]（加宽计算之子程序）计算内侧路基宽度14 JU-LI FWJ(任意两点间距离、方位角计算之子程序) 计算任意两点间的距离和方位角15 SQXBG-1（竖曲线标高计算之子程序）计算竖曲线设计标高16 SQX CAN SHU（竖曲线参数）输入各个竖曲线相关参数（按桩号先后顺序依次输入）17 CAN SHU（圆曲线参数）输入各个平曲线相关参数（按桩号先后顺序依次输入）18 S244-QIAO（S244 线桥梁坐标参数） S244 线桥梁坐标参数数据库（按桩号先后顺序依次输入）19 S244-BM（S244 线桥梁附近的导线点参数） S244 线桥梁附近的导线点参数数据库（按桩号先后顺序依次输入）常用字母的一般含义：（输入时）1. 判断参数M，当M≥0时，调用子程序相应参数；当M＜0时，直接输入相应参数。

CASIO4850程序（一）坐标反算：（ZBFS）（1）A“XA”:B“YA”:C“XB”:D“YB”S=POL（（C-A）,（D-B））▲E=Frac（J÷360+1）×360“FWJ=”:E▲注：（1）计算XA？→第一点X坐标，YA? →第一点Y坐标，XB? →第二点X坐标，YB? →第一点Y坐标（2）结果r→两点间距离，J→象限角，FWJ→坐标方位角（二）直线坐标：（ZXZB）（1）Lbl 0：｛K,D｝:S=K-〖起点里程如:1188000〗：A=〖坐标方位角（十进制）如：163.3395411〗：X=〖起点X坐标〗+S×cos A:Y〖起点Y坐标〗+S×sin A:D“BZ”:D=1＝＞｛L,G｝:L“ZL”:A[1]=A-90:N=X+L×cosA[1]:E=Y+L×sinA[1]:G=“YL”:A[2]=A+90:X[2]=X+G×cosA[2]:Y[2]=Y+G×sinA[2]:“ZC”:“N=”N=N▲“E=”:E=E▲“YC”:“X[2]=”:X[2]=X[2] ▲“Y[2]=”:Y[2]=Y[2] ▲Goto 0:≠＝＞“X=”:X=X▲“Y=”:Y=Y▲Goto 0注：（1）计算K? →里程BZ? →任意数中桩坐标，1边桩坐标：ZL? →左边桩，RL?→右边桩（2）结果N=左X坐标，E＝左Y坐标，X[2]=右X坐标，Y[2]=右Y坐标注明：程序修改方位角163°20′22.35″=（十进制）163.3395411（三）曲线坐标｛缓和曲线，圆曲线｝（1）主程序（STAKEOUT）Prog “DA.SKIP”Lbl 1:｛H,Y｝H“（CE SHE）K”Y“（CE SHE）E”S=H-Q:Z=FYProg “P”Pol（V-K,W-N）T＜0＝＞T=T-360⊿U=T:K=IU“（CE SHE）X”=A+KcosT“（CE SHE）X=”:U=U▲K“（CE SHE）Y”=Ksin T+B“（CE SHE）Y=”:K=K▲Goto 1（2）子程序（3） 1.（DA.SKIP）DegC“DATA IN”C=0＝＞Fixm:≠＝＞Prog“DA.IN”△2.（DA.IN）R“（BAN JING）R”L“（HUAN HE）L”F“（QU XIAN）FX”Q“（ZH）K”X“（HZ）K”S=（X-Q）÷2:Z=0Prog “XY.BAS”Z[1]=V:Z[2]=W:Z[3]=KD“（CE ZHAN）K”O“（CE ZHAN）E”S=D-Q:Z=FOProg “P”M=V:N=WE“（HOU SHI）K”G“（HOU SHI）FW”S=E-Q:Z=OProg“P”Pol（V-M,W-N）Fixm:P=J3.（P）S≤X-Q-L＝＞Prog “XY.BAS”:≠＝＞S=X-Q-S:Prog“XY.BAS”:Pol（V-Z[1],W-Z[2]）:Fixm:V=Z[1]-Icos （2FZ[3]-J）:W=Z[2]-ISin （2FZ[3]-J）:｛E,U,G｝:K=2Z[3]-K△4.（XY.BAS）A=A:B=B:S≤0＝＞K=0:V=S:W=FZ:≠＝＞S≤L＝＞K=90S^2÷πRL：V=S-S^5÷40R^2 L^2-Zsin K:W=FS^3÷6RL-FS^7÷336R^3L^3＋FZcosK△△S＞L＝＞K=90（2S-L）÷ΠR:V=L-L^3÷40R^2-Rsin（90L÷πR）+（R-Z）sinK:W=FL^2÷6R-FL^4÷336R^3+FRcos（90L÷πR）+F（Z-R）cosK△注：计算：①DATA IN? →模式输0以外的任意数修改曲线要素要素，输入0继续上次的计算。

主程序：ZBJSM:fixm:LbI0:M≥0=>{KS}:KS: ≠>{NE}:NEΔLbI1:F=0:G=0:L=0:Prog”R”:G=0=>G=LΔPol(Q-V,T-W):Z=J:Pol(V-X,W-Y):Q=Z-J:Z=J:Q<180=>Q=Q+360ΔQ>180=>Q=Q-360ΔT=R^-1=(R+TL^2÷24)tg(AbsQ÷2)+L÷2-T^2L^3÷240:U=K-U:U≤L=>X=Up rog”HU”:Goto2ΔU=U-L:J=兀RabsQ÷180-L:U≤J=>X=L:Prog”HU”:U=90U÷兀R:X=X+Rec(2RsinU,U+L):Y=Y+J:L=L+2U:Goto2ΔT=-T:Z=Z+Q:X=L-U+J:Prog”HU”:X=-X:L =-L:LbI2:Q<0=>Y=-Y:L=-L:J=G:G=F:F=IΔX=V+Rec(Pol(X-T,Y),J+Z):Y=W+J:L=L+Z:M<0=>K=K+Rec(Pol(N-X,E-Y),J-L):AbsI≥0.1=>Goto1Δ”K=”:K▲”S=”:S=J▲Goto3?≠>X=X+Rec(S,L+90):”X=”:X▲”Y=”:Y=Y+J▲LbI3:Prog”H”:H=(K-U)(W-V)÷(Z-U)+V:Z=Z-O:U=U+J:K>Z=>H=H+(K-Z)^2÷2RΔK<U=>H=H+(U-K)^2÷2IΔ”H=”:H= H-0.000▲Goto0子程序：HUX≤0=>F=2:G=2:Y=0:L=0:Goto2ΔF=0=>F=2:G=2:Goto1ΔX>G=>G=F:F=-F:Goto1ΔY=F:F=X(-Y-2)÷G+Z:G=4G÷(Y+2):X≥G=>G=-F:≠>G=ZΔLbI1:I=R^-1L^-1:Y=IX^3÷6-I^3X^7÷336:L=90IX^2÷兀:X=X-I^2X^5÷40:LbI2?子程序：DQ=V:T=W:V=X:W=Y:Z=U?子程序：CZ=U:W=V:R=I:O=J子程序：R(曲线要素数据库)P10 X=……:Y=……:Prog”D”:(交点坐标，后向前输)P9 X=……:Y=……:Prog”D”:P8 X=……:Y=……:U=……(P9交点的第一缓和曲线起点里程)：K≥U=>R=……(平曲线半径)：L=……(缓和曲线缓和段长度)：Goto1ΔProg”D”:PN (输入方法同P8)P1 X=……:Y=……:U=……：K≥U=>R=……：L=……:LbI1?子程序：H(竖曲线要素)P10 U=……(竖曲线变坡点里程，从后向前输)：V=……(设计高程)：I=……(竖曲线半径，凸为“-”、凹为“+”)J=……(切线长):Prog”C”:P9 U=……:V=……:I=……：J=……：K≥U=> Goto1ΔProg”C”:PN (输入方法同P9)P1 U=……:V=……:I=……：J=……：LbI1?注： K=所求点里程 S=边距，左“-”、右“+” M＝正算为“+”反算为“-”此程序可把全线路要素输入数据库，能执行直线、圆曲线、LS1=LS1的缓和曲线的坐标，高程放线正反算。