下载文档原格式
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。
CASIO4850全线坐标正反算程序主程序：ZHANG LIANGLbI Ø ：Deg:T"WJM":C=2Ø+4Ø(T-1):T=Ø:N"1.SR,2.JS":N=2=> Goto A⊿M"JZ-XX":O"QD-ZH":U"X-Ø ":V"Y-Ø ":G"FWJ ":Lbl B:{HRQP}: P"QD-R"：H"XX-CD":R"ZD-R":Q").-1,(.+1,--.Ø":W=H:Z =Ø:Prog"S1":Prog"S4":U=X:V=Y:O=O+H:G=F-9Ø:P=R:"…SR…NEXT…!”:Isz T: Dsz M:Goto B⊿Lbl A：Fixm:{XYSZN}:T=Ø:N"1.JS-XY,2.XY-FS": N=2=>Goto 2⊿S"JS-ZH=":Z"D=":Lbl 1:S＞Z［8T+1+C］=>Isz T:Goto 1⊿Prog "S3":W=S-O: Prog "S1": "JS-X=":X：Pause Ø:"JS-Y=":Y◢Goto AΔLbl 2: X"X…":Y"Y…":I=X:J=Y:Lbl 3: Prog "S3":W=Abs((Y-V)cos(G-9Ø)-(X-U)sin(G-9Ø:W＞H=> Isz T:Goto 3⊿Prog "S2":"FS-ZH…":S=O+W:Pause Ø:"D…":Z ◢Goto A⊿子程序S1A=.1184634425:B=.2393143352:N=.2844444444:K=. Ø4691ØØ77:L=.23Ø7653449:D=Q r W2(P-R)÷2HPR:E=Q r WP-1:Z[1]=G+KE+K2D:Z[2]=G+LE+L2D:Z[3]=G+.5E+.25D:Z[4]=G +(1-L)2D+E(1-L:Z[5]=G +(1-K)2D+E(1-K:F=G+E+D+9Ø：X=U+ZcosF+W(AcosZ[1]+BcosZ[2]+NcosZ[3]+BcosZ[4]+AcosZ[5:Y=V+ZsinF+W(AsinZ[1]+BsinZ[2]+NsinZ[3]+Bsin Z[4]+Asin Z[5：子程序S2W=Abs((Y-V)cos(G-9Ø)-(X-U)sin(G-9Ø:Z=Ø:Lbl 4:Prog "S1":L=G-9Ø+E+D:Z=(J-Y)cosL-(I-X)sinL:AbsZ<.ØØ1=>Goto5:≠>W=W+Z:Goto 4⊿Lbl 5:Z=Ø:Prog "S1":Z=(J-Y)÷sinF:子程序S3H=Z［8T+C+6:O=Z［8T+C+1］-H:U=Z［8T+C+2:V=Z［8T+C+3:G=Z［8T+C+4:P=Z［8T+C+5:R=Z［8T+C+7: Q=Z［8T+C+8:子程序S4Z［8T+C+1］=O+H: Z［8T+C+2］=U: Z［8T+C+3］=V:Z［8T+C+4］=G: Z［8T+C+5］=P:Z［8T+C+6］=H:Z［8T+C+7］=R: Z［8T+C+8］=Q: 结束程序编制说明：1.本程序是为CASIO-fx4850p编制的。
C A S I O F X4850公路计算程序0前言：传统公路测量中，使用的仪器设备和方法都很落后，需带着数学用表、曲线用表、计算盘、计算尺和算盘等一类的工具，完成外业测量工作。
计算器的出现，改变了这一局面。
高速公路建设中，长大曲线比比皆是，传统中对公路中线的测设方法，被极坐标法彻底的否定与取代，但大量的计算工作，只能带着提前计算好的线路逐桩坐标、高程资料，进行外业测量工作，机动性很差，现场查找也不方便。
这些问题都能在CASIO系列可编程计算器上得到很好的解决，对CASIO系列可编程计算器如何使用，直接影响到测量成果的质量和工作效率，本文将对CASIO系列可编程计算器快捷的计算方法进行分析与介绍。
1：以知线外任意点坐标，求对应线路里程在缓和曲线上，要计算任意里程的法线方向及任意宽度的边线坐标，非常简单。
但要计算任意一个已知坐标点，是对应哪一个里程法线方向上的点，就有一些困难。
很难推导一个这样的计算公式。
唯一的方法“渐进”，如果手工计算这可不是一个好方法。
但在有CASIO系列可编程计算器，如：FX-4500的情况下就变的非常简单了。
亦可用于直线和圆曲线的计算。
首先在缓和曲线上任选一点A为起始点，计算该点的坐标和切线方位角，通过坐标反算求起始点A与计算点B的方位角和距离，B点肯定对应A点切线方向上有一个垂足C点，把三点看成一个直角三角形，通过解直角三角形计算A C的距离，当该距离大于某一数值，如0。
001m，A点里程加AC的距离等于C点的里程，回到开始重新进入新一轮的计算，如果AC的距离小于某一规定值，则计算C点的里程与BC的距离即可。
求对应线路里程程序：主程序QLC (已知坐标求里程)Lb1 0：｛LＤE｝：Prog XH：Goto 0子程序：XH (循环)L1 Lb1 1L2 Norm：Prog LYYD：L3 PO1(D-X,E-Y)：W≤0=>W=W+360⊿L4 Z=W-I：A=V×cos Z：L=L+AL5 Abs A≥0.001=>Goto1：≠=>B=V×sinZ：Fix 3：“FXJL=”◢L6 L:Fix3:“DYLC=”◢程序中字母代表D 任意点X坐标，E 任意点Y坐标，DYLC 对应里程，FXJL 中线法线距离。
ZDZBZBJS（匝道中桩、边桩坐标计算）LbI 0：{EG}：A″XA″：B″YA″：C″CA″：D″1/RA″：E″1/RB″：F″DKA″：G″DKB″LbI 1：{HLR}：H″DKI″：L″DL″：R″DR″：H>G＝＞Goto 3⊿P=（E-D）÷Abs（G-F）：Q= Abs（H-F）：I=P×Q：T=D+IJ=C+（I+2D）Q×90÷π◢M=C+（I÷4+2D）×Q×45÷（2π）：N=C+（3I÷4+2D）Q×135÷（2π）K=C+（I÷2+2D）Q×45÷πX=A+Q÷12×（Cos C+4（Cos M+ Cos N）+2 Cos K+ Cos J）◢Y=B+Q÷12×（Sin C+4（Sin M+ Sin N）+2 Sin K+ Sin J）◢U″XL″=X+L Cos（J-90）◢V″YL″=Y+L Sin（J-90）◢W″XR″=X+R Cos（J+90）◢Z″YR″=Y+R Sin（J+90）◢A=X：B=Y：C=J：D=T：F=H：Goto 1LbI 3：A=X：B=Y：D=E：F=G：C=J：Goto 0说明：XA？、YA？——ZH或HZ的X、Y坐标；CA？——切线方位角；1/RA？——起算点的曲率：直线段代0，曲线段代+1/R（右转）或-1/R（左转）；1/RB？——代终点曲率：即±1/R；DKA？——代起算点里程；DKB？——代待求段的终点里程；DKI？——代所求点里程；DL？——代左边桩距中桩的距离；DR？——代右边桩距中桩的距离；——是回车键符号。
（本程序适合于任何线形的计算）（2）ZXZBZJS（直线段中桩、边桩坐标计算）G″KA″A″XA″B″YA″D″XB″E″YB″：POl（D-A，E-B）：J＜0＝＞J=J+360：⊿F=J◢LbI 0：{KSP}：H=K-G：L=Abs H：Rec（L，F）：X=A+I◢Y=B+J◢Rec（S，F-90）：M″XL″=X+I◢N″YL″=Y+J◢Rec（P，F+90）：V″XR″=X+I◢W″YR″=Y+J◢Goto 0说明：KA——起算点里程 XA、YA——起算点坐标 XB、YB——后视点坐标K ——待求点里程 X、Y——待求点中桩坐标 S——左侧桩距中桩距离XL、YL——左侧桩坐标 P——右侧桩距中桩距离 XR、YR——右侧桩坐标（3）DBXMJJS（多边形面积计算）G=0：A=X：B=Y：C=A：D=B：LbI 0：{XY}：X＜0＝＞G=G+AD-BC：A= AbsG/2：A″A=″◢A=0.0015A：A″AA=″◢⊿G=G+XD-YC：C=X：D=Y：Goto 0说明：X、Y——待求点坐标 A——面积（m2） AA——面积（亩）测点坐标输入完毕后在X中输入-1。
Fx4850最精简的万能坐标计算程序(线元法）1、JS-XY(计算主程序)Norm↙Z“KM”：D：N“PJ”:Prog“D”↙P=(S-R)÷T:Q= Z-C:I=PQ↙J=F+90Q(I+2R)÷π:J≥360=>J=J-360: ≠>J<0=>J=J+360⊿↙L=F+45Q(I÷4+2R)÷(2π):○=F+135Q(3I÷4+2R)÷(2π):K=F+45Q(I÷2+2R)÷π↙U=A+Q÷12×(cosF+4(cosL+cos○)+2 cos K+COS J):V =B+Q÷12×(sin F+4(sin L+ sin ○)+2 sin K+ sin J): W=J+N:X=U+D cosW:Y=V+D sin W:Fix 3: “XD=”:X▲“YD=”:Y↙本程序可以计算任何线型，包括直线、缓和曲线、圆曲线、不等半径曲线等，由其是匝道各曲线计算。
注：0为数字（零），○为字母（噢）KM-----计算点桩号（输入）D-------计算边桩距离左负，右正，中0（输入）PJ------路线右偏角，正交边桩为90（输入）XD------计算所得X坐标（计算得出）YD-----计算所得Y坐标（计算得出）FW---所求点桩号的切线方位角（计算得出）2、子程序：D（数据源）Z<803.961=>A=46670.1：B=13269.591:C=700：F=84°50′23″：R=0：S=0：T=103.961：Goto 0: ≠>Z<1085.902=> A=46679.45：B=13373.13:C=803.961：F=83°33′23″：R=0：S=0：T=281.941：Goto 0: ≠>Z<1363.366=> A=46711.091：B=13653.29:C=1085.902:F=83°33′23″：R=-1÷1000：S=-1÷1000：T=277.464：Goto 0:LbI 0:Fixm:A:B:C:F:R:S:T: ↙Z<本段线元的终点桩号=>A=起点X坐标：B=起点Y坐标:C=线元起点桩号：F=起点切线方位角：R=线元起点半径：S=线元终点半径：T=本段线元长度：Goto 0: ≠>Z<本段线元的终点桩号=>A=起点X坐标：B=起点Y坐标:C=线元起点桩号：F=起点切线方位角：R=线元起点半径：S=线元终点半径：T=本段线元长度：Goto 0: ≠>Z<？=>A=？：B=？:C=？：F=？°？′？″：R=？：S=？：T=？：Goto 0:Z<？=>A=？：B=？:C=？：F=？°？′？″：R=？：S=？：T=？：Goto 0:LbI 0:Fixm:A:B:C:F:R:S:T: ↙说明：字程序Z〈后面为每段线元的终点桩号，可以把你标段所有线元段都输入：其程序中的字母：A、B为起点坐标；C为起点桩号；F为起点方位角：R、S为起终点半径（以倒数输入，如R60，1÷60）分左右转如左转为-1÷60,右转为1÷60，直线为0；T为线元长度。
卡西欧fx-4850P万能放线程序主程序：ZBJS↙◢“◢：H◢↙注：↙↙↙↙注：此程序可把全线路要素输入数据库，能执行直线、高程放线正反算。
卵形曲线。
卡西欧fx-4850P直线、各种曲线任意放中线(匝道)1：主程序：TYQXJS“1，KS=＞XY”：“2，XY=＞KS”：N：U“X0”：V“Y0”：o“K0”：G“FWJ”：H“LS”：P“R1”：R“R2”：Q“Y+1；0；Z-1”：C=P-1：D=（P-R）÷（2HPR）：E=180÷π：N=1=＞Goto 1:≠＞Goto 2↙◢“Goto 1↙2：正算子程序：S U B1A=0.1739274226:B=0.3260725774:Z=0.0694318442:L=0.33 00094782:F=1-L:M=1-Z:X=U+W(Acos(G+QEZW(C+ZWD))+Bcos (G+QELW(C+LWD))+Bcos(G+QEFW(C+FWD))+Acos(G+QEMW(C+M WD))):Y=V+W(Asin(G+QEZW(C+ZWD))+Bsin(G+QELW(C+LWD)) +Bsin(G+QEFW(C+FWD))+Asin(G+QEMW(C+MWD))):F=G+QEW(C +WD)+90:X=X+Scos F :Y=Y+SsinF↙3坐标反算子程序 SUB2T=G-90：W=Abs((Y-V)cosT-(X-U)sinT) :S=0:Lbl 0: Prog “SUB1”：L=T+QEW（C+WD）:S=（J-Y）cos L-(I-X)sin L : Abs S＜10-6=＞Goto 1: ≠＞W=W+S:Goto 0↙Lbl 1:S=0:Prog“SUB1”:S=(J-Y)÷sin F↙注：o→代表：欧 0→代表零使用说明：1规定：⑴本程序可执行直线、圆曲线、缓和曲线（完整或非完整型）、复曲线、卵形曲线的任意里程中边桩坐标的正反算。
ZDXLbI 0Prog “FWJ”I”S”▲LbI 1｛ALP｝︰A”J”︰L︰E=X+Lcos (A+J) ︰F=Y+Lsin (A+J) ︰”X”︰E▲“Y”︰F▲P=1=＞Goto 0︰≠=＞Goto 1︰⊿FWJFix 4｛XYMN｝︰X”X0”︰Y”Y 0”︰M”X1”︰N”Y1”︰P ol(M-X,N-Y)FixmJ＜0=>J=J+360︰⊿”A”︰J→DMS▲“S”︰IHYQXLbI 3Fixm﹛K﹜︰Prog “DAT”Prog “QJS”“X”:X▲“Y”:Y▲Goto 3SQXFixmLbl 6Prog “DAT1”Prog “SJS”“H”:X▲Goto 6SDDMDefm 1LbI 0Prog “DAT1”Prog “SJS”Z[1]=XProg “DAT”Prog “QJS”（此处计算坐标必须是隧道中心点坐标）A=X:B=Y:C=D:F=Z[1]｛Z｝︰Z”DM”FixmLbI 1｛VGEH｝︰G”X”︰H”Y”︰E”H”︰V”JING”M=(G-A)cos C+(H-B)sin C︰L=(H-B)cos C-(G-A)sin C︰L=Abs L︰N=E-F︰”K-GAI”︰M▲Prog “DAT2”{W}:W=1=＞Goto 0︰≠=＞Goto 1︰⊿JSKDefm 3LbI 0﹛K﹜︰Prog “DAT”G=Z[3]Prog “QJS”A=X:B=Y:C=D:FixmZ[3]=G:G=Z[1]:H=Z[2]:｛GH｝︰G”X”︰H”Y”Z[1]=G:Z[2]=H:M=(G-A)cos C+(H-B)sin C:”K-GAI”︰M▲K=K+M:”LICH”:K▲L=(H-B)cos C-(G-A)sin C:”BIAN”︰L▲Goto 0QJSS=B²÷24R︰M=C²÷24R︰N=B÷2-B^3÷240R²︰O=C÷2-C^3÷240R²︰T=(R+M-(R+S)cos J)÷sin J+N︰U=(R+S-(R+M)cos J)÷sin J+O︰L=JRπ÷180-B÷2-C÷2+B+C︰H=K-I︰V=L-C︰H≤V=＞W=T︰Z=A+180︰≠=＞W=U︰Z=A+PJ︰⊿E=X+UcosZ︰F=Y+Usin ZH>0=＞H<B=＞Prog”H”︰⊿⊿H≥B=＞H≤V=＞Prog “y”︰⊿⊿H>V=＞H<L=＞H= L-H︰Prog “H”︰H=K-I︰Y=-Y︰D=2A-D+PJ︰⊿⊿H≤0=＞X=H︰Y=0︰D=A︰⊿H≥L＝＞X=L-H︰Y=0︰D=A+PJ︰⊿H≤V=＞O=A︰≠＝＞O=A+PJ +180︰⊿M=E+Xcos O-YsinO︰N=F+Xsin O+Ycos O﹛GQ﹜︰G”B”︰G=0=＞Q= 0︰≠=＞Q”J”︰⊿X=M+Gcos (D+Q) ︰Y=N+Gsin (D+Q) YX=N+Rsin (((2H-B)÷2R)×180÷π) ︰Y=P (S+R(1-cos (90(2H-B)÷πR))) ︰D=A+P (90(2H-B)÷Rπ)HX=H-H^5÷40R²B²+H^9÷3456R^4B^4︰Y=(H^3÷6RB-H^7÷336R^3B^3+H^11÷42240R^5B^5)P︰D=A+P(90H²÷RBπ)SJSW=C-B︰L=WR︰W＜0=＞G=-1︰L= -L︰≠＝＞G=1︰⊿T=L÷2︰L=WR︰E=T²÷2÷R︰H=K-I︰H≤0=＞F= BH︰⊿H>0=＞F=CH︰⊿Abs H<T=＞F=F+(T-A bsH)²÷2÷RG︰⊿X=F+AJDFixmO=0︰V≠1=＞Goto 2︰⊿Pol(N-P,L-Q) ︰J<S=>Goto 3︰⊿J>T=>Goto 3︰⊿o=1︰U=I-R︰U<0=>”R-WAI”︰≠=>”R- NEI”︰⊿U=Abs U︰U▲Goto 3Lbl 2U=R²-(N-P)²︰U<0=> “WUJIE”Goto 3︰⊿U=√U+Q ︰Pol(N-P,U-Q) ︰J<S=>Goto 3 ︰⊿J>T=>Goto 3 ︰⊿o=1 ︰U=L-U ︰U<0=>”H-WAI “︰≠=>”H-NEI”︰⊿U=Abs U ︰U▲Lbl 3DATK<135200＝＞X=1408︰Y=0︰C=280：B=280︰R=6000︰J=23º51º50º︰A=0︰I=131624.36︰P=-1︰Goto 0︰⊿K<137000＝＞X=4475.753︰Y= -1357.128︰A=336º08º10º︰J=20º29º59º︰B=310︰C=310︰R=5000︰I=13 5290.64︰P=1︰Goto0︰⊿K<158950＝＞X=158082.391︰Y= 10000︰A=11º16º41º︰J=11º16º41º︰B=260︰C=260︰R=7000︰I=15726 5.73︰P=-1︰Goto0︰⊿K<166710＝＞X=16096.745︰Y= 10000︰A=0︰J=13º34º14º︰B=280︰C=280︰R=6000︰I=164 591.784︰P=1︰Goto0︰⊿K≥166710＝＞X=14549.512︰Y= 7357.2034︰A=176º15º7.8º︰J=13º42º56.5º︰B=260︰C=260︰R=7000︰I=1671 80.249︰P=1︰Goto0︰⊿Lbl 0DAT1{K}:K<159649=>A=16.692 ︰B=0.003 ︰C=-0.003 ︰R=20000 ︰I=157950 ︰Goto 4 ︰⊿K<161960=>A=8.772︰B=-0.006 ︰C=0︰R=20000 ︰I=161900 ︰G oto 4 ︰⊿K<164558=>A=8.772︰B=0 ︰C=0︰R=20000 ︰I=164500 ︰Goto 4 ︰⊿K<165700=>A=15.152 ︰R=20000︰B=0.0058 ︰C=-0.002︰I=16 5600 ︰Goto 4 ︰⊿K<166750=>A=12.939︰B=-0.002 ︰C=0.003︰R=20000 ︰I=166700 ︰Goto 4 ︰⊿K<168287=>A=17.439︰B=0.003 ︰C=-0.0057︰R=20000 ︰I=16 8200 ︰Goto 4 ︰⊿Lbl 4DAT2Z=330=> Goto A︰⊿Z=320=> Goto B︰⊿Lbl AP=2.27︰Q=0︰R=7.01︰S=0︰T=60︰Prog “JD”:o=1=>Goto Z︰⊿P=2.27︰Q=0︰R=7.01︰S=60︰T=107º53º55º︰Prog “JD” :o=1=>Goto Z︰⊿P=1.217︰Q=3.273︰R=3.57︰S=107º53º55º︰T=157º35º24º︰Prog “JD”:o=1=>Goto Z︰⊿Lbl BP=2.27︰Q=0︰R=7.01︰S=0︰T=60︰Prog“JD”:o=1=>Goto Z︰⊿P=2.27︰Q=0︰R=7.01︰S=60︰T=107º53º55º︰Prog “JD” :o1=>Goto Z︰⊿P=1.217︰Q=3.273︰R=3.57︰S=107º53º55º︰T=157º35º24º︰Prog “JD”:o=1=>Goto Z︰⊿Labl Z程序使用说明此程序适应于CASIO 4850计算器ZDX程序（支导线计算）X0，Y0测站点坐标X1，Y1后视点坐标A 反算方位角S 反算距离J 观测左角L 观测距离X，Y 计算坐标FWJ程序（方位角计算）X0，Y0测站点坐标X1，Y1后视点坐标A 反算方位角S 反算距离HYQX程序（基本型曲线计算）K 计算里程B 计算边距J 右偏角（90为右侧，270为左侧，支持任意角）X，Y 计算坐标SQX程序（竖曲线计算）K 计算里程H 中桩高程SDDM程序（隧道断面计算）K 计算断面所在里程，半径越小，里程越需要准确，7000半径里程偏差十米，误差约1cm，直线无误差。
该程序是本人自己编的适用于casio4850全线坐标计算，发布上来给大家参考，请大家试用指教，如有问题请修改，首先是方位角显示不出来，请高手修改。
一、LXZB（坐标计算主程序）LbI Z：｛P｝A=×××：V=××°××′××″：C=×××：D=×××：P≤直线的终点桩号＝＞Prog“ZX.SUB”：Goto Z（以上各字母含义详直线坐标计算子程序说明）A=×××：V=××°××′××″：C=×××：D=×××：R=×××：I=（1或-1）P≤圆曲线终点桩号＝＞Prog“YQX.SUB”：Goto Z（以上各字母含义详圆曲线坐标计算子程序说明）A=×××：V=××°××′××″：C=×××：D=×××：B=×××：W=××°××′××″：E=×××：F=×××：R=×××：L=×××：I=（1或-1）P≤HY点桩号＝＞Prog“HQ1.SUB”：Goto ZP≤YH点桩号＝＞Prog“HYQ.SUB”：Goto ZP≤HZ点桩号＝＞Prog“HQ2.SUB”：Goto Z（以上各字母含义详（HQ1.SUB、HYQ.SUB、HQ2.SUB）坐标计算子程序说明）该主程序是简单的示例说明，具体的应根据图纸的线型组合，随机应变的做顺序调整，并将其输入Casio4850计算器。
4850线路计算程序1：已知线外任意点坐标，求对应线路里程在缓和曲线上，要计算任意里程的法线方向及任意宽度的边线坐标，非常简单。
但要计算任意一个已知坐标点，是对应哪一个里程法线方向上的点，就有一些困难。
很难推导一个这样的计算公式。
唯一的方法“渐进”，如果手工计算这可不是一个好方法。
但在有CASIO系列可编程计算器，如：FX-4500的情况下就变的非常简单了。
亦可用于直线和圆曲线的计算。
首先在缓和曲线上任选一点A为起始点，计算该点的坐标和切线方位角，通过坐标反算求起始点A与计算点B的方位角和距离，B点肯定对应A点切线方向上有一个垂足C点，把三点看成一个直角三角形，通过解直角三角形计算AC的距离，当该距离大于某一数值，如0。
001m，A点里程加AC的距离等于C点的里程，回到开始重新进入新一轮的计算，如果AC的距离小于某一规定值，则计算C点的里程与BC的距离即可。
求对应线路里程程序：主程序QLC (已知坐标求里程)Lb1 0：｛LＤE｝：Prog XH：Goto 0子程序：XH (循环)L1 Lb1 1L2 Norm：Prog LYYD：L3 PO1(D-X,E-Y)：W≤0=> W=W+360⊿L4 Z=W-I：A=V×cos Z：L=L+AL5 Abs A≥0.001=>Goto 1：≠=>B=V×sinZ：Fix 3：“FXJL=” ◢L6 L:Fix3:“DYLC=” ◢程序中字母代表D 任意点X坐标，E 任意点Y坐标，DYLC 对应里程，FXJL 中线法线距离。
程序中有坐标反算功能。
使用方法：只需输入计算点坐标、和较为接近的桩号。
桩号越接近计算速度越快2：逐桩坐标计算2.1编制方法：线路坐标程序是按照平曲线为单元，直线部分归属在曲线两端的方法，把整段路线分段装进数据库，根据桩号判断采用数据通过共用程序，进行任意点的坐标计算，在坐标转换示意土，第一直线段，是通过方位角和距离直接计算大地坐标，第一缓和曲线和圆曲线段，是先计算任意点切线支距和方位角然后转换大地坐标，第二缓和曲线段和直线段是先计算任意点切线支距和方位角。
fx--4850P圆曲线计算套件程序SJDYQX(实交单圆曲线主程序)Mcl: R:Z:H:V:W←Prog "ZBFS"←K=J:P=I←Prog "ZBFS"←R>800⇒S=20:≠>S=10△R<150⇒S=5△ G=J-K:G<0⇒G=G+180:≠>G=G-180△Q=Abs( G÷2):H=H÷2←Prog "YQXZD"←XJDYQX(虚交单圆曲线主程序)Mcl: R:Z:H:V:W←Prog "ZBFS"←K=J:P=I←Prog "ZBFS"←S=I：L=J：G=J-K:G<0⇒G=G+180:≠>G=G-180△E:F:Pol(E-C,F-D):J<0⇒J=J+360△E=J-L+180:E>360⇒E=E-360△E<0⇒E=E+180:≠>E=E-180△G=G+E:Q=Abs( G÷2):S=Ssin(Abs E)÷sin (180-2Q):"X.GJ":A=A-Scos K ▲"Y.GJ":B=B-Wsin K ▲ H=H÷2:P=P+S：R>800⇒S=20:≠>S=10 △ R<150⇒S=5 △ ←Prog "YQXZD"←YQXZD(圆曲线主点计算)"T":T=Rtan Q ▲ "L":L=2RQπ÷180▲"E":E=R(1÷cos Q-1)▲N=(Int(Z÷S))S:V≠0⇒U=(tan-1(W÷V))÷V △←Lbl 1： G<0⇒F=-90:"ZZJ":2Q▲≠>F=90:"YZJ":2Q ▲ △ "ZY.ZH":Z+P-T▲"QZ.ZH":Z+P-T+L÷2▲ "YZ.ZH":Z+P-T+L▲ "ZD.ZH":Z+P-2T+L+I ▲Z[1]=A+Tcos K:"ZY.X":Z[1]+(w÷2)cos(K-F)▲Z[2]=B+Tsin K:"ZY.Y":Z[2]+(w÷2)sin(K-F)▲ "QZ.X":A+(E+W÷2)cos(K-F+G÷2)▲ "QZ.Y":B+(E+W÷2)sin(K-F+G÷2)▲:Z[3]=A+Tcos J:"YZ.X":Z[3]+(w÷2)cos(J+F) ▲ Z[4]=B+Tsin J:"YZ.Y":Z[4]+(w÷2)sin(J+F) ▲Lbl 2:X=A+(R+E)cos(K-F+G÷2):Y=B+（R+E）sin (K-F+G÷2):M=K+F：H=0⇒Goto 5 △ G>0⇒Z[5]=0:Z[6]=W:≠>Z[5]=W:Z[6]=0 △"X.QZ":A+Pcos K+Hcos (K+90)▲ "Y.QZ":B+Psin K+Hsin (K+90)▲"X.QY":A+Pcos K+Hcos (K-90)▲ "Y.QY":B+Psin K+Hsin (K-90)▲"ZY.XZ":Z[1]+(H+Z[5])cos(K+90) ▲ "ZY.YZ":Z[2]+(H+Z[5])sin (K+90) ▲ "ZY.XY":Z[1]+(H+Z[6])cos(K-90)▲ "ZY.YY":Z[2]+(H+Z[6])sin (K-90) ▲Lbl 3:Z[7]=X+(R+H)cos(K+F+G÷2):Z[8]=Y+(R+H)sin(K+F+G÷2):Z[9]=X+(R-H-W)cos(K+F+G÷2):Z[10]=Y+(R-H-W)sin (K+F+G÷2):G>0⇒"QZ.XZ":Z[7] ▲ "QZ.YZ":Z[8] ▲ "QZ.XY":Z[9] ▲ "QZ.YY":Z[10] ▲:≠>"QZ.XZ":Z[9]▲"QZ.YZ":Z[10] ▲ "QZ.XY":Z[7] ▲ "QZ.YY":Z[8] ▲ △ ←Lbl 4:"YZ.XZ":Z[3]+(H+Z[5])cos(J-90)▲"YZ.YZ":Z[4]+(H+Z[5])sin(J-90)▲"YZ.XY":Z[3]+(H+Z[6])cos(J+90)▲ "YZ.YY":Z[4]+(H+Z[6])sin(J+90)▲"X.ZZ":A+Icos J+Hcos(J-90) ▲ "Y.ZZ":B+Isin J+Hsin (J-90)▲"X.ZY":A+Icos J+Hcos(J+90) ▲ "Y.ZY":B+Isin J+Hsin (J+90)▲Lbl 5：Prog "JKHHJS"←YQXZB(圆曲线坐标计算)Lbl 1: N=N+S:N>Z+L⇒Goto 4 △ ←Lbl 2: "ZH.P":N ▲ G>0⇒P=N-Z:≠>P=Z-N △ ○=180P÷R÷π:"X.P":X+(R-W÷2)cos (M+○)▲"Y.P":Y+(R-W÷2)sin (M+○)▲Lbl 3:H=0⇒Goto 1 △Z[7]=X+(R+H)cos(M+○):Z[8]=Y+(R+H)sin (M+○):Z[9]=X+(R-H-W)cos(M+○):Z[10]=Y+(R-H-W)sin(M+○):G>0⇒"X.ZP": Z[7]▲ "Y.ZP":Z[8] ▲ "X.YP":Z[9] ▲ "Y.YP": Z[10] ▲≠>"X.ZP":Z[9] ▲ "Y.ZP":Z[10] ▲ "X.YP": Z[7] ▲ "Y.YP":Z[8] ▲ △Goto 1 △ ←Lbl 4:N=N-S ←Lbl 5: Prog "JKHHJS" ←JKHHJS(加宽缓和计算含直线段)Z[12]=Z[12]+1:Z[12]=2⇒Goto 5△ C=A:D=B:A=A+Pcos K:B=B+Psin K:K=K+180:W=0⇒Goto 1:≠>"X.HQ":Z[1]=A+(P-T-V)cos K ▲"Y.HQ":Z[2]=B+(P-T-V)sin K ▲ "X.HQZ":Z[1]+Hcos (K-90) ▲"Y.HQZ":Z[2]+Hsin (K-90) ▲ "X.HQY":Z[1]+Hcos(K+90) ▲"Y.HQY":Z[2]+Hsin(K+90)▲ "X.HZ":Z[1]=C+(T+V)cos J ▲ "Y.HZ":Z[2]=D+(T+V)sin J ▲"X.HZZ":Z[1]+Hcos (J-90)▲ "Y.HZZ":Z[2]+Hsin (J-90)▲"X.HZY":Z[1]+Hcos (J+90)▲ "Y.HZY":Z[2]+Hsin (J+90)▲ △←Lbl 1: N=N+S:N≥Z+P-T⇒Goto 3 △ "ZH.P":N ▲ N>Z+P-T-V⇒○=(N-(Z+P-T-V))tan (U(N-(Z+P-T-V))):≠>○=0 △ H≠0⇒Goto 2:≠>"X.P":A+(N-Z)cos K ▲ "Y.P":B+(N-Z)sin K ▲ △ Goto 1 △←Lbl 2: G>0⇒Z[7]=A+(N-Z)cos K:"X.P":Z[7]+(○÷2)cos(K+90) ▲Z[8]=B+(N-Z)sin K:"Y.P":Z[8]+(○÷2)sin(K+90) ▲"X.ZP":Z[7]+Hcos(K-90) ▲ "Y.ZP":Z[8]+Hsin (K-90 )▲"X.YP":Z[7]+(H+○)cos (K+90)▲"Y.YP":Z[8]+(H+○)sin (K+90)▲≠>Z[7]=A+(N-Z)cos K:"X.P": Z[7]+(○÷2)cos(K-90)▲ Z[8]=B+(N-Z)sin K:"Y.P":Z[8]+(○÷2)sin(K-90)▲"X.ZP":Z[7]+(H+○)cos(K-90) ▲ "Y.ZP":Z[8]+(H+○)sin(K-90 )▲ "X.YP":Z[7]+Hcos (K+90)▲"Y.YP":Z[8]+Hsin (K+90)▲△ Goto 1 △← Lbl 3:N=N-S:Z=Z+P-T←Lbl 4: Prog "YQXZB"←Lbl 5: Z=Z+L:A=C:B=D：V≠0⇒U=(tan-1(W÷V))÷V：≠>U=0 △←Lbl 6:N=N+S:N≥Z+I-T⇒Goto 8 △"ZH.P":N ▲N<Z+V⇒○=(Z+V-N)tan (U(Z+V-N)):≠>○=0 △ H≠0⇒Goto 7:≠>"X.P":A+(N-Z+T)cos J ▲"Y.P":B+(N-Z+T)sin J ▲ △ Goto 6 △ ←Lbl 7: G>0⇒Z[7]=A+(N-Z+T)cos J:"X.P":Z[7]+(○÷2)cos(J+90) ▲Z[8]=B+(N-Z+T)sin J :"Y.P":Z[8]+(○÷2)sin (J+90) ▲"X.ZP":Z[7]+Hcos (J-90) ▲ "Y.ZP":Z[8]+Hsin(J-90) ▲ "X.YP":Z[7]+(H+○) cos (J+90) ▲ "Y.YP":Z[8]+(H+○)sin (J+90)▲ ≠>Z[7]=A+(N-Z+T)cos J:"X.P":Z[7]+(○÷2)cos(J-90) ▲ Z[8]=B+(N-Z+T)sin J:"Y.P":Z[8]+(○÷2)sin (J-90)▲ "X.ZP":Z[7]+(H+○)cos(J-90)▲ "Y.ZP":Z[8]+(H+○)sin(J-90)▲"X.YP":Z[7]+Hcos(J+90)▲"Y.YP":Z[8]+Hsin(J+90)▲ △ Goto 6 △←Lbl 8:"END"▲←X.ZZ、Y.ZZ=终点边左（x，y）X.ZY、Y.ZY=终点边右（x，y）ZY.XZ、ZY.YZ=直圆边左（ｘ，ｙ）ZY.XY、ZY.YY=直圆边右（ｘ，ｙ）QZ.XZ、QZ.YZ=曲中边左（x，y）QZ.XY、QZ.YY=曲中边右（x，y）YZ.XZ、YZ.YZ=圆直边左（ｘ，ｙ）YZ.XY、YZ.YY=圆直边右（ｘ，ｙ）X.P、Y.P =所求点中（x，y）XP.Z、YP.Z=所求边左桩点（x，y）XP.Y、YP.Y=所求边右桩点（x，y）单独显示内容：X.GJ、Y.GJ= 公交点（x，y）（在虚交单圆曲线中显示）。
CASIO fx4850万能计算程序[kaixin100sc正式版]卡西欧计算器程序2009-10-13 17:59:51 阅读236 评论0 字号：大中小经本人反复测试与验证,该程序可说是一个绝无仅有的万能程序,对于匝道当然也非常适用,匝道在立交工程中常常用到，由于其采用线型多样，计算复杂，因此一份适用的计算程序是不可缺少的，本文提供一份万能匝道计算程序，该程序简单，容易上手，输入参数少。
一、说明在高速公路和城市建设中常常运用到互通式立交工程，互通式立交由于曲线元的复杂性和多样性，给施工中的测量定位工作带来困难，手工计算已经不能满足要求，时下CASIO系列计算器在土木工程中已经得到普遍运用，为便于对立交工程匝道的准确测量定位，本人根据多年的测量工作经验，编写出一份适合CASIO fx-4850匝道计算程序，该程序界面好，能计算出匝道各种线元的中线坐标和边线坐标。
二、程序清单QXZDSJK (程序运行子程序,可计算n条线路)T=1=>Prog “SJK1”：Prog “QXZD”：Goto 0⊿(回车)T=n=>Prog “SJKn”：Prog “QXZD”⊿(回车)Lbi 0SJK (子程序--如果有多条线路，可改程序名为SJK1...到SJKn){H}：H“DKI”(回车)H≤曲线元终点里程=>A=曲线元起点坐标：B=曲线元终点坐标：C=曲线元起点正切线方位角：D=曲线元起点曲率：E=曲线元终点曲率：F=曲线元起点桩号：G=曲线元终点桩号：Goto 0⊿(回车)H≤曲线元终点里程=>A=曲线元起点坐标：B=曲线元终点坐标：C=曲线元起点正切线方位角：D=曲线元起点曲率：E=曲线元终点曲率：F=曲线元起点桩号：G=曲线元终点桩号：⊿(回车)Lbl 0(说明：一条线路中有多少个线元就要添加多少个进入数据库SJK中）QXZD (程序名)A"XA":B"YA": C"CA": D"RA":E"RB":F"DKA":G"DKB":Deg (回车)H>G=>{H}:Goto 0:≠>H<F=>{H}:Goto0:ΔL=(H-F)÷4:I=180L÷πD:U=90LL(D-E)÷πDE(G-F):T=C+I+U: W=C+2I+4U:K=C+3I+9U:O=C+4I+16U: ｛ZQ｝:Z"BJ(+0+)" :P=O+Q"JIAJIAO(-0+)":X=A+L(cosC+4(cosK+cosT)+2cosW+cosO)÷3+ZcosP:Y=B+L(sinC+4(sinK+sinT)+2sinW+sinO)÷3+Zsin P (回车)"X=":X▲"Y=":Y▲I=0:J=0:Pol（X-M"X0"，Y-N"Y0"）：J<0＝＞J=J+360：Δ“FWJ=”:J←DMS◢“I=”：I : Lbl 0三、变量及符号说明说明XA：YA：CA――――曲线元起点X、Y坐标及起点正切线方位角RA：RB――――――曲线元起点及终点曲率半径(左偏负，右偏正)DKA：DKB：―――曲线元起点和终点桩号DKI――――――――曲线元中待求点桩号Z―――――――――边桩距中线平距X：Y―――――――待求点中边桩桩号的X，Y坐标X0，Y0――――――置镜点的X，Y坐标FWJ――――――――待防样点的方位角L=―――――――――待防样点到置镜点的距离 Q"[JIA(-0+)]"--------------中点切线与两边桩连线的夹角(左负中0右正)四、示例龙门互通立交D匝道线型及相关数据如下(里面卖些关子,就是曲线元的半径没有给足,你会计算吗?)：点号桩号ＸY 方位角起点半径，止点半径QD, K0000 , 2777862.564 ,544381.888 ,295"56"01″1291,1291YH, K0081.721 ,2777900.604 ,544309.576 ,299"33"37″1291,250HY,K 122.038 ,2777921.779 ,544275.286 ,305"04"30″250,250YH, K 277.085,2778043.376 ,544183.130 ,340"36"33″250,10^45GQ,K 317.085 ,2778081.751 ,544171.883 ,345"11"34″10^45,-135HY,K 367.085 ,2778129.138 ,544156.172 ,334"34"57″-135,-135ZD, K 432.211 ,2778179.092 ,544115.375 ,306"56"31″五、计算步骤1、QD－YH曲线元，该段为圆曲线运行程序，按照提示输入曲线元QD纵横坐标及方位角，半径1291,QD和YH点桩号，再输入该曲线元内所求点桩号即可显示该中线坐标和切线方位角，可以输入YH点桩号K0 81.721验证计算出的结果，之后会提示输入“Z”值，即边线宽度(不分左右,均为正值)，"Q"值既中线与边桩连线的夹角(左负右正中0)。
曲线任意里程中边桩坐标正反算(CASIO-4850计算器)专用程序一、程序功能及原理1.功能说明：本程序由一个主程序(TYQXjs)和四个子程——正算子程序(SUB1)、反算子程序(SUB2)等构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。
本修改版程序既可实现正算全线贯通，亦可实现反算全线贯通。
本程序也可以在CASIO fx-4800P计算器运行，但由于容量较小，不建议使用。
2．计算原理：利用Gauss-Legendre 5点通用(可查阅<测绘工程>2006年第二期，由武汉大学花向红所著---道路平面测设计算的通用公式)公式正算线路中边桩坐标、线外测点至曲线元起点和终点的垂距的符号是否相异（即Dca×Dcb<0=>该测点在其线元内）进行判断并利用该线元要素反算中桩里程、支距，最后计算出放样数据。
二、源程序1.主程序(TYQXjs)Lbl 0：Defm 50:"1.SZ => XY,2.XY => SZ"：N：N=1=>Goto 1：≠>Goto 2Δ←Lbl 1：{SZ}：SZS≤Z[18]=>J=1：Prog “DA T1”：Goto 3ΔS≤Z[26]=>J=2：Prog “DA T1”：Goto 3ΔS≤Z[34]=>J=3：Prog “DA T1”：Goto 3ΔS≤Z[42]=>J=4：Prog “DA T1”：Goto 3Δ……………S≤Z[8（N+1）+2]=>J=n：Prog “DAT1”：Goto 3ΔLbl 3：W=Abs(S-O)：Prog "SUB1"："XS=":X=X◢"YS=":Y=Y◢”FWI=”: F=F-M :F→DMS◢Goto 4←┘Lbl 2：{XY}：XY：Z[4]=X：Z[5] =Y：N=0Lbl A：Isz N：N=J: Prog" DA T1": Z[11]=U：Z[12]=V: Z[9]=G:A= Z[9]-M:J=J+1: Prog "DA T1": Z[13]=U: Z[14]=V: Z[10]=G:B= Z[10]-M:Prog "ZX1" ：Z[6]×Z[7]<0=>J=N：Prog DA T1" ：Goto BΔA=Z[9]+M：B=Z[10]+M：Prog "ZX1" ：Z[6]×Z[7]<0=>J=N：Prog "DA T1" ：Goto B：≠> Goto AΔLbl B：Prog "SUB2"："S=":S=O+W◢"Z=":Z=Z◢Goto 2Lbl 4:J=0:I=Pol(X-Z[1],Y-Z[2]):F=J:F<0=>F=F+360Δ”DIST=”: I◢”FW=”:F→DMS◢Goto 12. 正算子程序(SUB1)A=0.1184634425:B=0.2393143352:N=0.2844444444:K=0.0469100770:L=0.2307653449:Z[3]=0.5:X=U+W(Acos (G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Ncos(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))): Y=V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+ Nsin(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)))：M”ANG=”: F=G+QEW(C+WD)+M：X=X+ZcosF：Y=Y+ZsinF3. 反算子程序(SUB2)M”ANG=”:T=G-M：W=Abs((Y-V)cosT-(X-U)sinT)：Z=0：Lbl 0：Prog "SUB1"：L=T+QEW(C+WD)：Z=(Z[5]-Y)cosL-(Z[4]-X)sinL：AbsZ<1E-6=>Goto1：≠>W=W+Z：Goto 0Δ←┘Lbl 1：Z=0：Prog "SUB1"：Z=(Z[5]-Y)÷sinF4. 垂距计算子程序(ZX1)Z[6]=(Z[5]-Z[12])COS A-(Z[4]-Z[11])SIN AZ[7]=(Z[5]-Z[14])COS B-(Z[4]-Z[13])SIN B5.曲线元要素数据库：DA T1Lbl1:J=1=>U=Z[8]:V=Z[9]:O=Z[10]:G=Z[11]:H=Z[12]:P=Z[13]:R=Z[14]:Q=Z[15]ΔJ=2=>U=Z[16]:V=Z[17]:O=Z[18]:G=Z[19]:H=Z[20]:P=Z[21]:R=Z[22]:Q=Z[23]ΔJ=3=> U=Z[24]:V=Z[25]:O=Z[26]:G=Z[27]:H=Z[28]:P=Z[29]:R=Z[30]:Q=Z[31]ΔJ=4=> U=Z[32]:V=Z[33]:O=Z[34]:G=Z[35]:H=Z[36]:P=Z[37]:R=Z[38]:Q=Z[39]Δ……………………………J=N=> U=Z[8N]:V=Z[8N+1]:O=Z[8N+2]:G=Z[8N+3]:H=Z[8N+4]:P=Z[8N+5]:R=Z[8N+6]:Q=Z[8N+7]Δ(注：如有多个曲线元要素继续添加入数据库DAT1中)Lbl 2 ：U"X0"：V"Y0"：O"S0"：G"F0"：H"LS"：P"R0"：R"RN"：Q：C=1÷P：D=(P-R)÷(2HPR)：E=1 80÷π三、使用说明1、规定(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线路往左偏时，Q=-1；当线路往右偏时，Q=1；当线路为直线时，Q=0。
请使用WORD的Times New Roman字体，若使用其它字体，字母O与数字0是区别不出的.切记切记!积分公式计算的线元法CASIO4850计算器程序第十一版1 (主程序)(使用时请运行这个)L1 LbI 0:{ QDV}L2 U<4=> Prog “4”△L3 U=5=> Prog “2”△L4 U=6=> Prog “2”△L5 Int U=7=> D=0:V=0△L6 U=8=>{M}: PSW:D=(B+M“01”P)÷Cos Abs(90-AbsV): Q=W+S√(1-M2)△L7 U>5=> Prog “4” △L8 LbI2:IntU=9=> PS: Prog“F”:A=R+90: Prog“9”:Goto 2△L9 LbI3:U=10=>PSWB: Prog“F”: Pol ( P-W,S-B):J<0=>J=J+360:△A=J: Prog“9”:Goto 3△L10 QDVHEFACZG: Frac Y=0=> Prog “3”: △O: Prog “10”L11 U=2=>T=T+ Int (-5+10 Ran# )÷1000: K=K+ Int (-5+10Ran# )÷1000: Prog “8”: Prog “5”△L12 U=3=> I= -S+2S Ran# :J=√(S2-I2): I= Int I÷1000:J= Int (-J+2J Ran# )÷1000:T+I: Pause 0:K+J: Pause 0:I: Pause 0:J◢√(I2+J2) ◢△L13LbI 1: IntU=7=> {DV}: Prog“10”: S=D:D=B÷2÷Sin Abs V:J=R:U=7.1=>D=B÷2:J=R+V-90: △ Prog “6”: Prog “8”: Prog “5”:D=2D:J=R-180: U=7.1=> J=R+V+90: △ Prog “6”: Prog “8”: Prog “5”: V≤0=>V=V+180: ≠>V=V-180: △D=S:Goto 1△L14 U<5=>Q=Q+B: △Goto 010（子程序）（为共用计算中桩/加转角/用正算计算边桩/反算/显示 的过程）L1 Prog “F”:U=11=>R=0:△J= R+V: Prog “6” : Prog “8” : Prog “5”9 （子程序）（为点垂桩垂距/涵交桩交角的共同运行部分）L1 Prog“7”: Prog“8”: Abs(J-R)<1=>Q=Q-I:≠> Q=Q+I:△I<0.001=>T=P:K=S: Prog “8”: “OK”: Q: Pause 0:I: Pause 0: J=J-R:J>180=>J=-J+180:△J→DMS◢U=9.1=> Prog “11”: △△ Prog “4”7 （子程序）（两点斜直线交点部分）(也可作为两点斜直线交点主程序独立运行)L1E=Tan A:F= Tan R:Y=(EP-S-FT+K) ÷(E-F): Pause 0:O=E(Y-P)+S5 (子程序)( 显示部分)L1T: Pause 0:K: Pause 0:I: Pause 0:J→DMS◢F (子程序) (中桩计算部分)L1 RAD:I= Abs(Q-H):J= Abs(G-H):R=Aπ/180+CI+I2(Z-C)/2/J:T=E+∫(Cos (Aπ/180+CX+X2(Z-C)/2/J),0,I,4):K= F+∫(Sin (Aπ/180+CX+X2(Z-C)/2/J),0,I,4):Q>G =>R=R+π:△DEG:R=180R/πF2 (子程序) (中桩计算部分)(使用时把F2改为F;同时把F改为F1)(这个子程序计算速度快,但在缓曲长到2或3百米且半径小到6或7十米时,计算误差有几厘米)L1 Defm 7L1 I=Q-HL2 J= (Z-C)I÷(G-H)L3 Z[1]=A+90(J÷8+2C)(I÷8)÷πL4 Z[2]=A+90(J÷4+2C)(I÷4)÷πL5 Z[3]=A+90(3J÷8+2C)(3I÷8)÷πL6 Z[4]=A+90(J÷2+2C)(I÷2)÷πL7 Z[5]=A+90(5J÷8+2C)(5I÷8)÷πL8 Z[6]=A+90(6J÷8+2C)(6I÷8)÷πL9 Z[7]=A+90(7J÷8+2C)(7I÷8)÷πL10 R=A+90(J+2C) I÷πL11 T=E+ AbsI(Cos A+4(CosZ[1]+ CosZ[3]+ CosZ[5]+ CosZ[7])+2(CosZ[2]+ CosZ[4]+ CosZ[6])+ CosR)÷24L12 K=F+ AbsI(Sin A+4(SinZ[1]+ SinZ[3]+ SinZ[5]+ SinZ[7])+2(SinZ[2]+ SinZ[4]+ SinZ[6])+ SinR)÷248(子程序) (座标反算部分)(也可作为座标反算主程序独立运行, 方位角按RCL R SHIFT °′″ 显示)L1J=TK: Pol ( T-Y,K-O): Pause 0:J<0=>J=J+360△6(子程序) (边桩和座标正算计算部分)(也可作为座标正算主程序独立运行)L1TK:T=T+ Rec (D,J): Pause 0:K=K+J2(子程序) (推算开挖/填筑桩部分){平原微丘一般推算1-2次,山岭重丘一般4-5次,能达到5厘米以内的精度,看你的熟练程度} L1 {WSPB}：I= -W+ S◢L2 D= Abs IP+B◢L3 {D}:D3(导线点数据库)(数据库中的数字请按你所在公路的数据替换)L1 Y=21=>Y=36640.666:O=65111.666△L2 ……(省略号表示类同上行, 并非真在程序中输入省略号. 注意: 第一个Y为导线点编号, 第二个Y为导线点X座标,O为导线点Y座标, 数据库中的数字请按你所在公路的数据替换. 但注意: 导线点编号必须没字母也没小数;导线点X座标则必须有小数,若刚好没有小数,请加上0.0001并不影响精度)4(平曲线数据库) (数据库中的数字请按你所在公路的数据替换)L1 Q≤444.978=>H=374.978:E=36771.643:F=64929.942:A=151°29′38″:C=0: G=444.978:Z=-600-1: Goto 1△L2 Q≤811.819=>H=444.978: E=36710.800:F=64964.534:A=148°09′06″:C=-600-1: G= 811.819:Z=-600-1: Goto 1△L3 Q≤881.819=> H=811.819:E=36475.596:F=65238.599:A= 113°07′15″:C=-600-1: G=881.819:Z=0: Goto 1△L4 ……(省略号表示类同上行,并非真在程序中输入省略号)L5 LbI 1 (注意:这行始终是平曲线数据库的最后一行, 之前行都是相同的) (为更好的理解这个平曲线数据库的结构,请对照上面数据看下图)11 （子程序）（双圆垂轴隧道断面开挖超欠测量部分）A:实测高程 H:设计高程I:实测宽度 B:超欠值(设棱镜半径0.03)L1 {AHI}:A≤H=>B=110:≠>A≤H+2.75+1.056=>B=√((A-(H+1.905))2+(I+1.813)2) +0.03-4.5 : ≠> B =√ ((A-(H+2.75))2+I2 )+0.03-2.5:△△B◢平曲线数据库的数据图示:(图画得不好,现丑了.图示是对称缓曲,但线元法程序对任意公路线形都适用的)ZH HZ 0+374.978X=36771.643 0Y=64929.942 HY151°29′38″ 0+444.9780 X=36710.800Y=64964.534-600-1 X=36475.596Y=65238.599113°07′15″-600-1使用说明:{本程序是针对CASIO4850计算器的，其它的计算器请作相应修改}{有方框的字为一个字符，不是几个字母，为各种计算状态名或函数名，如RAD与Cos等等，∫为积分符号，L1\ L2\ L3\L4等为程序行号不要在程序中输入, ……省略号表示类同上行,并非真在程序中输入省略号}{如有其它疑问或进行技术交流，请打詹德军的电话015918402837欢迎交流!}运行主程序 1 (按1 FILE进入)U?(判断进入何种计算状态)1为数据库输入原始数据计算中边桩，2为数据库输入计算放线测量记录表，3为数据库输入计算平面位置检查记录表，4为人工输入原始数据计算中边桩，5为人工输入推算开挖/填筑桩状态，6为数据库输入推算开挖/填筑桩状态，7为数据库输入计算斜交斜做涵洞/单跨桥七点放样状态 ，7.1为数据库输入计算斜交正做涵洞七点放样状态 ，(斜交斜做涵洞的伸缩缝处的放样也可用)8为数据库输入计算锥坡放样，9为数据库输入反算线外点垂桩垂距，9.1为数据库输入双圆垂轴隧道 断面开挖 超欠测量，10为数据库输入求线外两点交路线的桩号与夹角，11为计算曲线桥平行布桩的桩/柱放样状态(因为D=0时就是中桩,D≠0时就是边桩,所以不再单独设立中桩或边桩的计算状态,方便程序进一步开发)（U=1）数据库输入原始数据计算中边桩1.U？1(判断进进入何种计算状态)2.Q？所求桩号（计算第二个桩开始Q已经自动加了步长，不用手工输入，但仍然会显示你看，并可修改成别的桩号）3.D？中桩到边桩的距离（求中桩输入0）4.V？中桩到边桩方位角与中桩前进切线方位角间的夹角(注意左转为负值,角度按需要以度分秒输入任意输入)5.Y?置仪点X座标 （也可以输入导线点编号，因为此时使用了导线点数据库）6.O?置仪点Y座标（当X座标输入了导线点编号，则此行不显示，因为数据库已经装入了导线点的X与Y座标）7.=所求桩号X座标 （注:四个结果同时显示于屏幕）=所求桩号Y座标=所求桩号放样距离=所求桩号放样方位角8．B？步长（步长就是间隔长度，当要计算的是大量有相同距离间隔 且 相同中边距或都是中桩 时，步长就是好东西了，计算完第一个桩后所有桩，都只需闭着眼睛连按四次EXE，就是下个桩要计算的所有结果了）9．以后循环第2/3/4/7步（U=2）计算放线测量记录表1．U？2(判断进进入何种计算状态)2．Q？所求桩号（计算第二个桩开始Q已经自动加了步长，不用手工输入，但仍然会显示你看，并可修改成别的桩号）3．D？中桩到边桩的距离（求中桩输入0）4．V？中桩到边桩方位角与中桩前进切线方位角间的夹角(注意左转为负值,角度按需要以度分秒输入任意输入)5．Y?置仪点X座标 （也可以输入导线点编号，因为此时使用了导线点数据库）6．O?置仪点Y座标（当X座标输入了导线点编号，则此行不显示，因为数据库已经装入了导线点的X与Y座标）7．=所求桩号计算X座标 （注:四个结果同时显示于屏幕）=所求桩号计算Y座标=所求桩号计算放样距离=所求桩号计算放样方位角8．=所求桩号复测检查X座标 （注:四个结果同时显示于屏幕）=所求桩号复测检查Y座标=所求桩号复测检查放样距离=所求桩号复测检查放样方位角9．B？步长（步长就是间隔长度）10．以后循环第2/3/4/7/8步（U=3）计算平面位置检查记录表1．U？3(判断进进入何种计算状态)2．Q？所求桩号（计算第二个桩开始Q已经自动加了步长，不用手工输入，但仍然会显示你看，并可修改成别的桩号）3．D？中桩到边桩的距离（求中桩输入0）4．V？中桩到边桩方位角与中桩前进切线方位角间的夹角(注意左转为负值,角度按需要以度分秒输入任意输入)5．Y?置仪点X座标 （也可以输入导线点编号，因为此时使用了导线点数据库）6．O?置仪点Y座标（当X座标输入了导线点编号，则此行不显示，因为数据库已经装入了导线点的X与Y座标）7．=所求桩号设计X座标 （注:四个结果同时显示于屏幕，但后两个数据是用不上的） =所求桩号设计Y座标=所求桩号放样距离=所求桩号放样方位角8．？所允许的最大误差（以毫米为单位，比如高速公路的中线允许的最大误差为20MM，就输入20）9．=所求桩号检测X座标 （注:四个结果同时显示于屏幕）=所求桩号检测Y座标=所求桩号X轴偏差=所求桩号Y轴偏差10．=所求桩号双轴偏差11．B？步长（步长就是间隔长度）12．以后循环第2/3/4/7/9/10步（U=4）人工输入原始数据计算中边桩1．Q？所求桩号（计算第二个桩开始Q已经自动加了步长，不用手工输入，但仍然会显示你看，并可修改成别的桩号）2．D？中桩到边桩的距离（求中桩输入0）3．V？中桩到边桩方位角与中桩前进切线方位角间的夹角(注意左转为负值,角度按需要以度分秒输入任意输入)4．H?起算点桩号5．E?起算点X座标6．F?起算点Y座标7．A?起算点方位角(以度分秒输入)8．C?起算点曲率(即半径分之一，直线为0，注意左转为负值)9．Z?终算点曲率(即半径分之一，直线为0，注意左转为负值)10.G?终算点桩号11.Y?置仪点X座标12.O?置仪点Y座标13.=所求桩号X座标 （注:四个结果同时显示于屏幕）=所求桩号Y座标=所求桩号放样距离=所求桩号放样方位角14.？步长（步长就是间隔长度）15.以后循环第1/2/3/13步（U=5）人工输入原始数据推算开挖/填筑桩状态(平原微丘一般推1-2次，山岭重丘4-5次达到5厘米的精度，要看你的熟练程度)1．W?实测大地高程2．S?变坡点设计高程3．=变坡点设计高程 与 实测大地高程 之高差(正为填方,负为挖方)4．P?填方或挖方的坡比(若1:1.5侧输1.5,无论填挖坡比全为正)5．B?中桩到变坡点宽度6．=推算的中桩到开挖/填筑桩距离7. D?采用的中桩到开挖/填筑桩距离7．Q?所求开挖/填筑桩桩号8．V?中桩到边桩的方位角与中桩的前进切线方位角间的夹角(注意左转为负值,只能输入右边90或左边-90，因为开挖/填筑桩总是垂直的)9．H?起算点桩号10.E?起算点X座标11.F?起算点Y座标12.A?起算点方位角(以度分秒输入)13.C?起算点曲率(即半径分之一，直线为0，注意左转为负值)14.Z?终算点曲率(即半径分之一，直线为0，注意左转为负值)15.G?终算点桩号16. Y?置仪点X座标17. O?置仪点Y座标18.=推算开挖/填筑桩桩号X座标 （注:四个结果同时显示于屏幕）=推算开挖/填筑桩桩号Y座标=推算开挖/填筑桩桩号放样距离=推算开挖/填筑桩桩号放样方位角12．以后循环第1/2/3//4/5/6/7/8/18步（U=6）数据库输入原始数据推算开挖/填筑桩状态(平原微丘一般推1-2次，山岭重丘4-5次达到5厘米的精度，要看你的熟练程度)1．W?实测大地高程2．S?变坡点设计高程3．=变坡点设计高程 与 实测大地高程 之高差(正为填方,负为挖方)4．P?填方或挖方的坡比(若1:1.5侧输1.5,无论填挖坡比全为正)5．B?中桩到变坡点宽度6．=推算的中桩到开挖/填筑桩距离7．D?采用的中桩到开挖/填筑桩距离8．Q?所求开挖/填筑桩桩号9．V?中桩到边桩的方位角与中桩的前进切线方位角间的夹角(注意左转为负值,只能输入右边90或左边-90，因为开挖/填筑桩总是垂直的)10．Y?置仪点X座标 （也可以输入导线点编号，因为此时使用了导线点数据库）11．O?置仪点Y座标（当Y输入了导线点编号，则此行不显示，因为数据库已经装入了导线点的X与Y座标）12．=推算开挖/填筑桩桩号X座标 （注:四个结果同时显示于屏幕）=推算开挖/填筑桩桩号Y座标=推算开挖/填筑桩桩号放样距离=推算开挖/填筑桩桩号放样方位角 12．以后循环第1/2/3//4/5/6/7/8/11步推算开挖/填筑桩的图示(以有二级填方边坡的填筑线为例,其它类推)中桩S变坡点高程中桩至该变坡点距离 B地面线（U=7/7.1）为计算斜交斜做涵洞单跨桥 / 斜交正做涵洞七点放样状态(U=7/7.1并能避免曲线内单跨桥涵进出水口长度不一至)(U=7同时适用于,曲线桥扇形布桩的盖梁放样; 曲线桥扇形布桩的支座放样)1.U? 7/7.1进入计算斜交斜做涵洞单跨桥/ 斜交正做涵洞七点放样状态2.Q? 涵洞桩号3.Y?置仪点X座标 （也可以输入导线点编号，因为此时使用了导线点数据库）4. O?置仪点Y座标（当Y输入了导线点编号，则此行不显示，因为数据库已经装入了导线点的X与Y座标）5.=涵洞中桩X座标 （注:四个结果同时显示于屏幕）=涵洞中桩Y座标=涵洞中桩放样距离=涵洞中桩放样方位角6.D?右侧涵斜长7.V? 右侧转角(进流水轴 与 中桩的前进切线方位角间的夹角,左负右正)8.= 右侧桩X座标 （注:四个结果同时显示于屏幕）=右侧桩Y座标=右侧桩放样距离=右侧桩放样方位角9.B?正的跨径长度(程序自算斜跨径长度)10. =右侧前进侧台身脚X座标 （注:四个结果同时显示于屏幕）=右侧前进侧台身脚Y座标=右侧前进侧台身脚放样距离=右侧前进侧台身脚放样方位角11. =右侧后退侧台身脚X座标 （注:四个结果同时显示于屏幕）=右侧后退侧台身脚Y座标=右侧后退侧台身脚放样距离=右侧后退侧台身脚放样方位角12.D?左侧涵斜长13.V? 左侧转角(出流水轴 与 中桩的前进切线方位角间的夹角,左负右正,程序自动按180度转换了,其实不用输入,只是显示观看)14.= 左侧桩X座标 （注:四个结果同时显示于屏幕）=左侧桩Y座标=左侧桩放样距离=左侧桩放样方位角15. =左侧前进侧台身脚X座标 （注:U=7.1时为后退侧结果，四个结果同时显示于屏幕）=左侧前进侧台身脚Y座标=左侧前进侧台身脚放样距离=左侧前进侧台身脚放样方位角16. =左侧后退侧台身脚X座标 （注: U=7.1时为前进侧结果，四个结果同时显示于屏幕）=左侧后退侧台身脚Y座标=左侧后退侧台身脚放样距离=左侧后退侧台身脚放样方位角17.以后循环第6-8/10-17步(当然这种循环只是重新检查而已,按两次AC/ON退出)计算涵洞单跨桥七点放样图示（注: U=7.1时左侧的15步与16步顺序互换）前进方向10步15步8步 右侧左侧14步 5步11步16步切线 路线（U=8）为数据库输入原始数据计算锥坡放样1．U?8为进入数据库输入原始数据计算锥坡放样2．P？椭圆单侧斜长轴3．S？椭圆单侧斜短轴(当从近桥0点往远桥1点推桩号增加时，S输入负值。
ZHY程序FX5800
辛普森计算器公路测量常用程序修改版一、程序功能
LbI 0:＂X=＂?A: ＂Y=＂?B：＂Z=＂?C：＂1/（R1）＂?D：＂1/（R2）＂?E：＂QT=＂?F: ＂QT=＂?G: ＂V＂?V：＂W＂?W LbI 1：＂L＂?L：＂O＂?O：＂H＂?H：IF L>G：Then Goto 1：IFEnd
（E-D）/（G-F）→P：L-F→Q：P*Q→I
C+（I+2*D）*Q*（90/∏）→J
C+（I/4 +2*D）*Q*（45/（2*∏））→M
C+（3/4*I+2*D）*Q*（135/（2*∏））→N
C+(I÷2+2*D)*Q*(45÷∏) →K
＂X=＂:A+Q÷12*(COS（C）+4*（COS（M）+C0S（N））+2*COS （K）+COS（J）)+O*COS（J+H）→X◢
＂Y=＂:B+Q÷12*(SIN（C）+4*（SIN（M）+SIS（N））+2*SIN （K）+SIN（J）)+O*SIN（J+H）→Y◢
X-V→I：Y-W→J
PoL（I，J）
”J=”:J▲DMS▲
If J <0:Then “J=”:J+360°→J:J▲DMS ▲ if End Goto1
F5800计算器曲线计算JD程序（ZH~HZ）坐标计算放样公式
主程序：
30→Dim Z: "XJD="?B: "YJD="?D: "KJD="?M: "A="?A: "ZJ="?
Q: "R="?R: "LS1="?C: "LS2="?T: "F="?F:prog "A":Lb1 0: "XO="?G: "YO="?H: "LC="?L: "JJ="?Z: "PJ="?S:If L≤M-Z[11]+C:Then M- Z[11]COS(F) →N:D- Z[11]Sin(F) →E:Prog "ZHHY":Prog "FWJ": Goto 0:Else If L≤M-Z[11]+ Z[13]-T:
Then M- Z[11] →K:B- Z[11] COS(F) →N:D- Z[11]Sin(F) →E: Prog "HYYH":Prog "FWJ": Goto 0: Else If L≤M-Z[11]+ Z[13]:Then M- Z[11]+ Z[13] →K:F+AQ+180→V:B-Z[12] COS(V) →N:D- Z[12] Sin(V) →E: Prog "YHHZ":Prog "FWJ": Goto 0: Else "End" ▲
If End: If End: If End
说明：以上主程序需要四个副程序，缺少任何一个副程序不能运行！
"XJD="?: "YJD="?: "KJD="?为曲线交点坐标及桩号
"A="?为曲线左转取-1；右转取1
"ZJ="?为曲线转角值；"R="?为曲线半径
"LS1="?： "LS2="?为前后段缓和曲线长度
"F="?曲线起始方位角； " X0="?、" Y0="?:全站仪测站点坐标
" LC="?为曲线ZH-HZ段内所求任意点桩号；
"JJ="?为切线顺时针夹角；"PJ="?为左右侧所求点距设计线距离；。
