曲线任意里程全线中边桩坐标正反算及放样(CASIO fx-4800P计算器)程序档

合集下载

FX-4800P计算器在测量放样中的应用

FX-4800P计算器在测量放样中的应用摘要:本文介绍了CASIO FX-4800P计算器在测量放样中的应用.关键词: 对称、不对称曲线,已知坐标、反算对应中桩里程、坐标、边距，起坡线、开挖线放样。

一、概述测量放样工作是工程建设中的一项重要工作。

哪里有建设工程施工，哪里就必须有测量人员去进行测量放样工作，任何工程如果没有测量人员首先在现场放样，那么这项工程其他工序是很难开展工作的。

测量放样工作好比建设者的眼睛，给他们指明了工作的目标，因此在工程建设中的重要性是不可忽视的。

放样工作及时性牵涉到工程的进展，如有延误，就会使施工计划受到影响，工程的质量要求测量放样点线的正确性达到百分之一百的保证率。

在这种情况下，为了使测量放样人员更好更快做好放样工作，随着可编程计算器（FX-4800P）在测量放样中的应用。

我们根据放样的内容，利用FX-4800P可编程功能适时计算出各放样内容点的坐标，特别是公路的缓和曲线，起坡线、开挖线，反算中桩里程、坐标、边距放样，在外面施工放样中用FX-4800P编程实用率很高。

二、对称、不对称曲线计算程序在各类建筑工程中，常常会遇到设计和各种曲线形建筑物，要求在现场进行测量放样，尤其在道路的施工中，频繁出现各种曲线，担任放样的工作人员必须掌握各类曲线的计算，能采用相应的测量放样方法，特别高速公路和立交桥匝道各种曲线特别多，这个程序就可以满足这种放样的要求。

Lb1 0:{DEG}:ABCDEFG:Lb1 1:{H}:Q=Abs(H-F:I=Q(E-D)÷Abs(G-F:R=45QI÷32л:Z=45QD÷2л:U=C+64R+8Z◢X=C+R+Z:Y=C+9R+3Z:L=C+25R+5Z:T=C+49R+7Z:K=C+4R+2Z:O=C+16R+4Z: P=C+36R+6Z:M=A+Q(cosC+4(cosX+cosY+cosL+cosT)+2(cosK+cosO+cosP)+cos U)÷24:N=B+Q(sinC+4(sinX+ sinY+ sinL+ sinT)+2( sinK+ sinO+ sinP)+ sinU)÷24:Lb1 2:{S}:X=M-SsinU◢Y=N+ScosU◢S≠0 =>Goto 2 Lb1 3:{PLK}:PLK:M=X+LcosU÷sinP+Kcos(U+P)÷sinP◢Z=Y+LsinU÷sinP+Ksin(U+P)÷sinP◢L≠0 =>Goto 3 K≠0=>Goto 3H≠G =>Goto 1 A=X:B=Y:F=G:C=U:Goto 01、字符说明：A、B为起点x、y坐标；C为起点方位角；D、E为起、终点曲率1/R(左偏角输入“－”，右偏角输入“+”值)；F、G为起、终点里程桩号；H为待求桩号；S为待求点距中桩边距(左－、右+、中0)；P桥墩轴线(从左右)与中桩切线方向的夹角；L 为桥涵上各点至桥中轴线的垂直距离（前+、后－）；K为桥涵上各点至中桩切线方向的垂直距离（左－、右+）；U为待求点的计算切线方位角(U0时U=U－360)；X、Y为各待求点计算坐标；M、Z为桥涵上各待求点计算坐标；2、文字说明：①、只可直线、缓和曲线、圆曲线、卵型曲线分线段依次计算，至本段终点后自动进入下段计算，然后输入下一段对应曲率；②、求桥涵时(弯桥直做)先计算全桥的中心桩号坐标,后输入P、L、K值；③、求桥涵时(弯桥弯做)计算每一墩号中桩坐标,后输入P、L、K值；④、将L、K值输0即可退出桥涵状态返回路基放样状态(P不可为0)；3、曲率说明：①、直线段D=0，E=0；②、独立圆曲线段D=E=1/R(Ls1=Ls2=0)；③、对称曲线ZH~HY段D=0,E=1/R;HY~YH段D=E=1/R;YH~HZ段D=1/R,E=0(Ls1=Ls2≠0)；④、不对称曲线【两段圆曲线共用缓和曲线】ZH~HY段D=0,E=1/R1;HY~YH段D=1/R1,E=1/R2;HY~YH段D=E=1/R2；(Ls1≠Ls2≠0)⑤、不对称(卵型)曲线【两段圆曲线共用缓和点】ZH~HY段D=0,E=1/R1;HY~YH段D=1/R1,E=1/R2;YH~YH段D=E=1/R2；(Ls1≠0, Ls2=0)三、起坡线、开挖线放样计算程序在公路开工前，测量人员须先进行测量放出起坡线和开挖边线，由于施工现场是多变的，测量工作人员必须配合开挖工作人员的现场开挖，随时测出开挖边坡线，这个程序就能很好地完成这些工作须要。

曲线(含直线)任意里程中边桩坐标正反算

曲线（含直线）任意里程中边桩坐标正反算（CASIO fx-4800p&fx-4850）J-SQX(竖曲线数据输入)“J-SQX”：{NHUVMQP}：Z[9]=N“SJD”：Z[10]=H“JDZ”：Z[11]=U“I1”：Z[12]=V“I2”：Z[13]=M“R”：Z[14]=Q“QD”：Z[15]=P“ZD”：“TO J-JSMS”J-PQX(平曲线数据输入,自动切换到J-JSMS)A“JD”B“JDX”C“JDY”F“FJ”O“ZJ：Z-，Y+”RE“LS1”K“LS2”：E＜1=＞E=1E-9⊿K＜1=＞K=1E-9⊿Z[1]=EE÷(24R)-E^4÷(2688RRR)：Z[2]=E÷2-EEE÷(240RR)：X=(EE-KK)÷(24R)÷sinAbsO：Z[3]“T1”=(R+Z[1])tan(AbsO÷2)+Z[2]-X▲Z[4]“T2”=(R+KK÷(24R)-K^4÷(2688RRR))tan(AbsO÷2)+K÷2-KKK÷(240RR)+X▲L=AbsOπR÷180+(E+K)÷2▲J=tan-1((R+Z[1])÷(Z[3]-Z[2]))：X“E”=(R+Z[1])÷sinJ-R▲X=A-Z[3]：Y=X+E：E＜1=＞X“ZY”▲≠=＞X“ZH”▲Y“HY”▲⊿Y“QZ”=X+(L-K-E)÷2+E▲Y=X+L-K：X=X+L：K＜1=＞X“YZ”▲≠=＞Y“YH”▲X“HZ”▲⊿Prog“J-JSMS”J-JSMS(放样模式主程序)“1-ZS,2-FS,4-DMFY” Lb1 0：{Z}：Z“MS”≤1=＞Goto 1：≠=＞Z“MS”≤4=＞Goto 2⊿⊿Lb1 1：{PDW}：PD“BZ”W“BJ”：Prog“JP”：X“X=”▲Y“Y=”▲Goto 0⊿Lb1 2：{XYW}：XYW“BJ”： Prog“JF”：P“P=”▲D“BZ=”▲Z=3=＞Prog“JS”：Prog“DMFY”⊿(运行竖曲线高程计算程序)Z=4=＞Z[26] “H”=12.417+(P-75360)*5.2/1000: Prog“DMFY”⊿（运行单面坡比高程计算,语句中12.417为起点桩号DK75+360的高程,可以根据实际情况进行调整；5.2/1000为单面上坡率,上坡输正值，下坡输负值）Goto 0JS(竖曲线计算主程序)P＜Z[14] =＞Prog“J-SQX”⊿P＞Z[15] =＞ Prog“J-SQX”⊿N=Z[9]：U=Z[11]：V=Z[12]：H=Z[13]：G=H Abs(U-V)÷200：P＜N=＞I=U：M=N-G：M＞P=＞M=P⊿≠=＞I=V：M=N+G：M＜P=＞M=P⊿⊿J=(P-M)2÷(2H)：U-V＞0=＞J=-J⊿H=Z[10]+(P-N)×I÷100+JZ[26] “H”=HJP(平曲线正算子程序)FixmLb1 1：J=F：X=B-Z[3]cosF：Y=C-Z[3]sinF：G=EP≤A-Z[3] =＞I=A-Z[3]-P：M=-I：N=0：H=F+W：Goto5≠=＞P≤A-Z[3]+E=＞I=P-A+Z[3]：H=90II÷(REπ)：O＜0=＞H=-H⊿H=H+W+F： Goto 3≠=＞P≤A-Z[3]+L-K=＞I=P-A+Z[3]-E： Goto 4：≠=＞Goto 2⊿⊿Lb1 2：X=B+Z[4]cos（F+O）：Y=C+Z[4]sin（F+O）：J=F+O+180：G=KP≤A-Z[3]+L=＞I= A-Z[3]+L-P：H=90II÷(REπ)：O＞0=＞H=-H⊿H=H+J+W+180： Goto 3≠=＞I=P-A+Z[3]-L：M=-I：N=0：H=J+W+180： Goto 5Lb1 3：M=I-I^5÷(40RRGG)：N=III÷(6GR)-I^7÷(336RRRGGG)： Goto 5Lb1 4：H=(E+2I)×90÷(πR)：M=RsinH+Z[2]：N=R(1-cosH)+Z[1]：O＜0=＞H=-H⊿H=J+W+H： Goto 5Lb1 5：P≤A-Z[3]+L-K=＞O＜0=＞N=-N⊿≠＞O＞0=＞N=-N⊿⊿Goto 6Lb1 6：X=X+Mcos J-Nsin J+Dcos H：Y=Y+Msin J+N cosJ+Dsin HJF(平曲线反算子程序)FixmU=X：V=Y：D=0：J=F-W：P=A+(Y-C)cos J-(X-B)sin JLb1 1：Prog“JP”：J=H-180：I=(V-Y)cos J-(U-X)sin J：Abs I＜1E-4=＞Goto 2：≠=＞P=P+I： Goto 1⊿Lb1 2：D=(V-Y)÷sin HDMFY隧道开挖断面(超欠挖情况)计算程序{HIJMNG}：Z[20]= H“C”：Z[21]=I“A”：Z[22]= J“H1”：Z[23]= G“GC”：Z[24]= M“R1”：Z[25]= N“R2”D≤2.3=＞D“ZD”=Abs(D-2.3)▲≠=＞D＞2.3=＞D“YD”= D-2.3▲⊿⊿(C为圆心O1与O2的高差值,A为圆心O1与O2的宽度差值,H1为圆心O1至设计高程的高度,GC为实测高程,R1为上部第1个半径,R2为下部第2个半径)(显示值ZD为至隧道中心线左边的宽度,YD为至隧道中心线右边的宽度,2.3为设计线路距隧道中心线的宽度,可以实际情况进行改变调整)Z[26]“H”▲(显示所求桩号点设计高程值)Z[27]“O1” =Z[26]+ Z[22](第1个圆心高程)Z[28]“O2”=Z[26]+ Z[22]+ Z[20] (第2个圆心高程)Z[29]= Z[24]×Sin 60▲(此数据为上半弧60度范围的宽度值，可根据实际情况调整；数值可不显示) D>Z[29] =＞Goto 1⊿Z[30]=√（（Z[24]）2-D2）Z[31]= √（Abs(（Z[24]）2-(Z[23]-Z[27])2)）Z[32] “HGD”=Z[27]+ Z[30]- Z[23] ▲Z[33] “SKD”=Z[31]-D ▲Prog“J-JSMS”Lb1 1Z[23]<Z[27] =＞Goto 2⊿Z[31]= √(（Z[25]）2-(Z[23]-Z[28])2)Z[33] “SKD”=Z[21]+ Z[31]-D ▲Prog“J-JSMS”Lb1 2Z[31]=Z[25]-（Z[27]-Z[23]）*0.35/(0.78+Z[20]+Z[22])(0.35为下边墙往边墙底内缩的距离，0.78为设计高程至边墙底部的高度)Z[33] “SKD”= Z[31]-D ▲Prog“J-JSMS”一、程序简介1、本套程序共有2个主程序，5个子程序。

4800程序

卡西欧Fx-4800P计算器测量程序一、坐标放样源程序ZBFYU“Xc”:M“Yc”:O“XH”:Q“YH”:I=O-U:J=Q-M: Pol(I，J):J＜0＝＞J=J+360:≠＞J=J:⊿T“HF” =J◢I“HL”◢X=0=＞W“XQ” =U+Lcos(T+P“a”)◢S“YQ” =M+Lsin(T+P)◢≠＞I=X-U:J=Y-M: Pol(I，J):J＜0＝＞J=J+360:≠＞J=J⊿I“QL=”◢H“QF” =J◢J=H-T:J＜0=＞J=J+360◢≠＞J=J使用说明:XH—后视点X坐标 YH—后视点Y坐标XC—置镜点X坐标 YC—置镜点Y坐标HF—后视方位角 QF—前视方位角HL—后视距离 QL—前视距离X、Y—前视点坐标 J—仪器后视归零后拨角当X=0时运行坐标反算QL—测设点到置镜点的距离а—测设点与后视点夹角XQ—测设点的X坐标 YQ—测设点的Y坐标二、复合辛卜生公式计算曲线元坐标源程序XBS—XY(已知里程求坐标)A“QD”:Z“ZD”:R“R0”:V“RN”:D“X0”:E“Y0”:F“F0”:K:L“B”:P“а”:J:Prog “SIM1”:W＜0＝＞W=W+360:≠＞W=W⊿W“FWI=”◢X=X◢Y=Y◢XBS—L、B(已知坐标求里程和边距)A“QD”:Z“ZD”:R“R0”:V“RN”:D“X0”:E“Y0”:F“F0”:K”K”(输入大致里程)I=X:B=Y:K=A+Abs((Y-E)cos(F+P“a”-180)+(X-D)sin(F+P-180):Lbl 1:L=0:Prog“SIM1”:C=(B-Y)cos(W+P-180)- (I-X)sin(W+P-180): AbsC＜10∧-4＝＞Goto2: ≠＞K=K+C∶Goto1:⊿Lbl2:L=0:Prog“SIM1”:C=(B-Y)÷sin(W+P):K=K◢C“B=”◢子程序SIM1T=(J÷V-J÷R)÷(Z-A)H=(K-A)÷NS=90H÷πW=F+((TNH+J÷R)+J÷R)NS: G=1:X=D+H÷6×(cosF+cosW+4×Σ(cos(F+((G+0.5)HT+2J÷R)(G+0.5)S),G,0,(N-1))+2Σ(cos(F+(THG+2J÷R)GS,G,1, (N-1)))+Lcos(W+P):Y=E+H÷6×(sinF+sinW+4×Σ(sin(F+((G+0.5)HT+2J÷R)(G+0.5)S,G,0,(N-1))+2×Σ(sin(F+(THG+2J÷R)GS,G,1,(N-1)))+Lsin(W+P)使用说明:QD—起点里程 ZD—终点里程R0—起点半径 RN—终点半径J—曲线转向 (左转为-1 直线为0 右转为+1)X0、Y0—线元起点坐标 F0—起点方位角K—待求点里程 N—计算精度(可取4～12的整数)а—边桩与中线右交角 B—边桩到中桩的距离，左“-”右“+”X、Y —待求点坐标 K—待求点里程三、数据库的编制以每个线元为一段先付参数值，例K≥542.759＝＞K≤686.148＝＞A”起点里程”=542.7590:Z终点里程=163.5:R起点半径=10^45: V终点半径=10^45:D起点坐标=22236.147:E起点坐标=32145.669:F起点切线方位角=81°22′42″:J＝0:⊿⊿K≥987.324＝＞K＜10234.021＝＞A=……J=－1⊿⊿……Prog“XBS—XY”四、纵断高程计算Lbi 0R:T:D”K(JD)”:H”H(JD)”:U”I1”:V”I2”:｛K｝V＞U=＞G=1:≠＞V＜U=＞G=-1⊿⊿K≤D-T=＞L=D-KP=H-LU/100≠＞K≤D=＞L=K-D+TP=H-U(D-K)/100+GL2/2R≠＞K≤D+T=＞L=D+T-KP=H+V(K-D)/100+GL2/2R≠＞K＞D+T=＞L=K-DP=H+LV/100⊿⊿⊿⊿P”H=”▲Goto 0H:高程、 K：变坡点里程、H：变坡点高程、I1、I2前后坡度五、中轴旋转、超高计算: (1)IR,(2)IL线路左转输2，W（Z）：中分带宽度，K：桩号，H（Z）：中桩设计高程，W（L）左侧宽度，W（R）：右侧宽度Lbi 1G”(1)IR,(2)IL”:D”K(ZH)”:S”LS1”:H”LS2”:Q”LY”:W”W(Z)”:I”I1”:J”IB”｛KVAE｝:K:V”h(Z)”:A”W(L)”:E”W(R)”G=2=＞G=-1⊿｛KPAB｝:K:P”H(Z)”:A”W(L)”:B”W(R)”Prog”XZ”:U”H(L)=”▲V”H(R)=”▲Goto11、子程序：XZG=1=＞U=A：V=B：≠＞G=-1=＞U=B：V=A⊿⊿K≤D=＞R=P-I（U-W/2）/100：E=P-I（V-W/2）/100：≠＞K≤D+S=＞Prog”ZXA”: ≠＞K≤D+S+Q=＞R=P+J(U-W/2)/100:E=P-J(V-W/2)/100: ≠＞K≤D+S+Q+H=＞Prog”ZXB”: ≠＞K＞D+S+Q+H=＞R=P-I(U-W/2)/100:E=P-I(V-W/2)/100⊿⊿⊿⊿⊿G=1=＞U=R:V=E: ≠＞G=-1=＞U=E:V=R⊿⊿2、子程序： ZXAL=K-DR=P-I(U-W/2)/100+L(U-W/2)(I+J)/100SL≤2IS/(I+J) =＞E=P-I(V-W/2)/100: ≠＞L＞2IS/(I+J) =＞E=P-LJ(V-W/2)/100S ⊿⊿3、子程序： ZXBL=D+S+Q+H-KR=P-I(U-W/2)/100+L(U-W/2)(I+J)/100HL≤2IH/(I+J) =＞E=P-I(V-W/2)/100: ≠＞L＞2IH/(I+J) =＞E=P-IJ(V-W/2)/100H ⊿⊿六：边轴旋转超高计算：L：平曲线长度，I1正常横坡，IB超高横坡。

CASIO-4800完整缓和曲线任意中边桩坐标计算(绝对正确实用)

回
Lbl3:Q”X(Z1)”=C+T*COS(A+180) ▲V”Y(Z1)” =D+T*SIN(A+180) ▲K=A: GOTO D
Lbl 4: Z≤Z【2】=>GOTO 5: ≠=>Z＞Z【2】=>GOTO 7
Lbl 5: F=(Z- Z【1】)-( Z- Z【1】)^5/(40R^2Y^2)
I=(Z- Z【1】)^3/(6RY)- (Z- Z【1】)^7/(336R^3Y^3):M=A+30*J*( Z- Z【1】)^2/(∏RY)
Lbl A: Z﹤Z【4】=>GOTO B: ≠=>Z≥Z【4】=> Q”X(Z2)”=C+T*COS B▲V”Y(Z2)” =D+T*SIN B▲K=B: GOTO D
Lbl B: F=(Z【4】- Z)-( Z【4】- Z)^5/(40R^2Y^2) :I=( Z【4】- Z)^3/(6RY): M=B-30*J*( Z【4】- Z)^2/(∏RY)
Lbl 6: Q”X(H1)”=C+T*COS(A+180)+ ABS(F^2+I^2)^1/2*COS M▲V”Y(H1)” =D+T*SIN(A+180) + ABS(F^2+I^2)^1/2*SIN M▲K=A+90*J( Z- Z【1】)^2/(∏RY): GOTO D
Lbl 7: Z≤Z【3】=>GOTO 8: ≠=>Z＞Z【3】=>GOTO A
CASIO-4800完整缓和曲线任意中边桩坐标计算(绝对正确实用)
DEFM 4
Lbl 0:{Z,S,O,U}
A:B:C”JDX”:D”JDY”:E”JD”:R:Y”LS”:T”LC”:L:J

CASIOf4800可编程计算器曲线放样程序

曲线计算公式、程序及算例广西建设职业技术学院测量教研室李向民2007年4月一、圆曲线中桩与边桩坐标计算公式1．根据已知的中线偏角、半径和交点桩号计算切线长T 、曲线长L 及外距E ，以及曲线主点的桩号切线长 2tan αR T = 曲线长180παR L =外距 )12(sec2cos-=-=ααR R R E切曲差（超距） D = 2T －L 式中α以度为单位。

交点的桩号已由中线丈量得到，根据交点的桩号和曲线测设元素，可计算出各主点的桩号，ZY 桩号 =JD 桩号－TQZ 桩号 = ZY 桩号＋2LYZ 桩号 = QZ 桩号＋2L为了避免计算中的错误，可用下式进行计算检核：JD 桩号 = YZ 桩号－T ＋D2．根据已知的两个交点的坐标，计算曲线主点的坐标和细部桩点的坐标（1）计算直圆点（ZY ）坐标根据JD 1和JD 2的坐标(x 1 , y 1)、(x 2 , y 2)，用坐标反算公式计算第一条切线的方位角α2-1和ZY 点坐标（x ZY , y ZY ），212112arctanx x y y --=-α122122sin cos --+=+=ααT y y T x x ZY ZY（2）计算圆心坐标因ZY 点至圆心方向与切线方向垂直，其方位角为αZY-O =α2-1 ± 90°（左偏角时用“+”，右偏角时用“－”）则圆心坐标（x o , y o ）为ozy ZY o o zy ZY o R y y R x x --+=+=ααsin cos（3）计算圆心至各细部点（和主点）的方位角设ZY 点至曲线上某细部里程桩点的弧长为l i ，其所对应的圆心角i β按下式计算得到：πβ180⋅=R l i i 则圆心至各细部点的方位角αi 为αi =（αZY-O ＋180°）±i β （左偏角时用“－”，右偏角时用“+”）（4）计算各细部点（和主点）的坐标根据圆心至细部点的方位角和半径，可计算细部点坐标io i i o i R y y R x x ααsin cos +=+=（5）计算各桩点边桩坐标设左边桩与中桩之间的垂直距离为D 左，右边桩与中桩之间的垂直距离为D 右，则左边桩的坐标为：io i i o i D R y y D R x x ααsin )(cos )(左左左左±+=±+=（左偏角时用“－”，右偏角时用“+”）右边桩的坐标为：io i i o i D R y y D R x x ααsin )(cos )(右右右右±+=±+=（左偏角时用“+”，右偏角时用“－”）注：可用全站仪按极坐标法测设线路主点和细部中桩，以及边桩；二、圆曲线计算程序（CASIO fx-4800P ）“YQS ”文件（主程序，计算圆曲线、直线的中桩坐标、边桩坐标）{}{}{}[][][][][][][][]{}{}2▲)90sin()(""▲)90cos()("":1:10:)sin )((:)cos )((:,211▲sin )(""▲cos )("")(180)(,11180sin cos 90sin cos 12""4:212""3:""2:180""1:)2tan(:,31:3:20""Pr :"":"":"":"",,,,2211Goto H H C B AbsH N Y V C B AbsH M X U C C H B K I P N B K I D M HK H K Lb Goto H H Q HZ R N Y V Q HZ R M X U OZ C Q R L O T I K L H K Lb J C J R F N J R E M Z B J B T P F B T D E Z Z YZ Z Z Z QZ Z T I ZY Z RA L Z A R T ZA Z A Lb Goto Goto R R OB og I Y P X D Y Y X X I P D Y X -=++=++=-=≠⇒=⇒≥-+=-+=-=++=++=⨯+=÷=--=+=+=+=⨯-=+=+=+=÷+=-=÷=÷=≠⇒⇒=ππ“OB ”文件（子程序，计算两点间的坐标方位角）BB Lb Goto W B W B P Y Lb Goto W B Goto D X D X P Y W Lb Goto B B P Y Lb Goto Goto D X =+=≠⇒=⇒>-+=≠⇒⇒>--÷-==≠⇒=⇒>-≠⇒⇒=----:414:360tan :tan 0:314:180tan :30:)()(:214:270:900:11:2:10111注释:R 为圆曲线半径（当算直线时输入0） I 为2JD 桩号1X ,1Y 是1JD 的XY 坐标2X ,2Y 是2JD 的XY 坐标A 为偏角Z ：右偏角时Z=1，左偏角Z=－1 K 为待放点桩号H 为边桩距：当为路线左边桩时符号为正，右边桩时符号为负，求中桩时输入0。

CASIO fx-4800p编程计算路线中桩及边桩坐标

豳３
断面成果圈
Ｋ３＋３２５０２．１９４３８５０８３９０７８２９８ห้องสมุดไป่ตู้６０８３５０７１５３１．０３３２９８６６．９。１０．３１．２９８．８。１８．１
Ｋ３＋３７５ｏ２．ｌ９５８０５０８４４２３９２９８６８．９５０７４１８２．０３１２９８６７．１．１１．０１．６９８８０．１９．９
维普资讯
・
３２・５
第３２卷第ｌ４期２００６年７月
山西建筑
ＳＨ￣ＮＸＩＡＲＣＨＩＴＥＣＴＵＲＥ
Ｖ０．２Ｎｏ１１３．４
Ｊ１２０．０６ｕ
・
计算机技术与应用・
中图分类号：７ＴＵ１文献标识码：Ａ
引言
在路基施工前，需要恢复路基中桩、放样路基边桩以便征地
方便施工测量的同时也加快了测量速度。在计算路基边桩的过程中，关键是计算路线横断面方向方位角。对于缓和曲线段可根据
拆迁工作进行，但是设计单位一般只提供路基中桩坐标，先测缓和曲线上任意一点的切线方向与ｘ轴的夹角为该点到缓和曲线若起点的长度的平方与回旋线参数平方的２倍的比值 … 。而横断面１量出路基中线后再放样路基边线比较费时费力，别是需要加密特中桩时。若预先计算出路基边桩坐标，再采用全站仪来直接放样方向方位角可根据直缓点到交点的方位角和缓和曲线角求得。圆路基边桩以完成规划红线放样及征地拆迁工作则比较方便。在施工实践中用ＣＳＯ一８０ＡＩｆ４０ｐ计算器编制程序计算路基边桩。ｘ

卡西欧Casiofx-4800P道路中边桩坐标计算程序

卡西欧Casiofx-4800P道路中边桩坐标计算程序程序清单如下：Casio fx-4800P道路中边桩坐标计算程序程序目的：依平曲线要素计算直线、圆曲线、缓和曲线的任意中桩、左、右桩坐标。

程序说明：K0：起始桩号X0：起始X坐标Y0：起始Y坐标ALF：起始方位角R：半径LS：缓和曲线长N：曲线左转N=1，右转N=2K：待求桩号LL、LR：左、右桩距离Q：左、右桩与中线斜交角求得XZ、YZ、XL、YL、XR、YR分别为中桩、左、右桩坐标。

一、直线段文件名：ZX (COMP)程式：L”K0”:O”X0”:P”Y0”:W”ALF”:Lbl 0:{K}:X”XZ”=O+(K-L)cosW◢Y”YZ”=P+(K-L)sinW◢{B}:S”XL”=X-B”LL”cos(W+Q) ◢T”YL”=Y-Bsin(W+Q) ◢{C}:U”XR”=X+C”LR”cos(W+Q) ◢V”YR”=Y+Csin(W+Q) ◢Goto 0注：在程序执行过程中，赋给的要素变数的值被固定不变，可对变数（K、LL、LR）赋予不同值，迅速求得所需坐标。

二、圆曲线段文件名：YQX (COMP)程式：L”K0”:O”X0”:P”Y0”:W”ALF”: Lbl 1:{K}:J=(-1)^N (K-L)÷R×180÷π:D=2Rsin((-1)^N J÷2):X”XZ”=O+Dcos(W+J÷2) ◢Y”YZ”=P+Dsin(W+J÷2) ◢{B}:S”XL”=X-B”LL”cos(W+J+Q) ◢T”YL”=Y-Bsin(W+J+Q) ◢{C}:U”XR”=X+C”LR”cos(W+J+Q) ◢V”YR”=Y+Csin(W+J+Q) ◢Goto 1注：若没有直接HY点方位角，则ALFHY=ALFZH±Ls/2/R×180/π,(左转-，右转+)。

三、缓和曲线文件名：HHQX (COMP)程式：L”K0”:O”X0”:P”Y0”:W”ALF”: M”LS”:Lbl 3:{K}:I=(-1)^N×(K-L)^2÷M÷R÷6×180÷π:D=(K-L)-(K-L)^5÷90÷(RM)^2:X”XZ”=O+Dcos(W+I)◢Y”YZ”=P+Dsin(W+I) ◢{B}:S”XL”=X-B”LL”cos(W+3I+Q) ◢T”YL”=Y-Bsin(W+3I+Q) ◢{C}:U”XR”=X+C”LR”cos(W+3I+Q) ◢V”YR”=Y+Csin(W+3I+Q) ◢Goto 3注：1、坐标计算方法是根据偏角法原理；2、缓和曲线（ZH~HY或YH~HZ）以ZH（或HZ）为起始点；3、平曲线左转（ZH~HY段N=1，YH~HZ段N=2），曲线右转（ZH~HY段N=2，YH~HZ段N=1）。

卡西欧计算器坐标的正反算

可以算任意斜交涵洞轴线的坐标，增加T为斜交角度，规定T为涵轴右侧方向与“线路前进方向切线”之间的夹角，当涵轴与线路正交时，T=90,其他操作与原程序一样；1. 正算子程序(SUB1) [color=Red]A=0.26：B=0.74：K=0.02：L=0.82：F=1-L：M=1-K：X=U+W(Acos(G+57.2958QKW(1/P+KWD))+Bcos(G+57.2958QLW(1/P+LW D))+Bcos(G+57.2958QFW(1/P+FWD))+Acos(G+57.2958QMW(1/P+MWD)))：Y=V+W(Asin(G+57.2958QKW(1/ P+KWD))+Bsin(G+57.2958QLW(1/P+LWD))+Bsin(G+57.2958QFW(1/P+FWD))+Asin(G+57.2958QMW (1/P+MWD)))：F=G+57.2958QW(1/P+WD)+90：X=X+Zcos(F-90+T)：Y=Y+Zsin(F-90+T)2. 反算子程序(SUB2)W=Abs((Y-V)cos(G-90)-(X-U)sin(G-90))：Z=0：Lbl 0：Prog "SUB1"：L=(G-90)+5 7.2958QW(1/P+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二.增设数据库程序（ＳＪＫ主程序）Lb1 4："1.SZ => XY"："2.XY => SZ"：｛NS｝：S∠下一线元起点里程＝＞O =本线元起点里程：U=本线元起点X：V=本线元起点Y：G=本线元起算方位角：H =本线元长度：P=起点曲率半径：R=终点曲率半径：Q=0或1、-1：Prog“TYQXJS”：Goto0Δ←┘(第一线元数据要素)S∠下一线元起点里程＝＞O=本线元起点里程：U=本线元起点X：V=本线元起点Y：G=本线元起算方位角：H=本线元长度：P=起点曲率半径：R=终点曲率半径：Q=0或1、-1：Goto0Δ←┘(第二线元数据要素)S∠下一线元起点里程＝＞O=本线元起点里程：U=本线元起点X：V=本线元起点Y：G=本线元起算方位角：H=本线元长度：P=起点曲率半径：R=终点曲率半径：Q=0或1、-1：Goto0Δ←┘(第三线元数据要素)。

CASIO-fx4800P型计算器连续计算线路高程、坐标Word文档

CASIO-fx4800P型计算器连续计算线路高程、坐标及放样程序【简述】为了适应现场快速、准确、灵活的放样要求，作者根据线路施工测量特点，运用计算器的编程功能，使整条线路的施工测量计算数据能连续快速的计算，只要把整条线路各个交点处的计算要素输入子程序中，你便拥有整条线路的“数据库”。

在忙碌的施工现场，你不必携带大量的施工图和资料来查找平面、纵断面的设计要素，更不必担心设计要素的输入错误，从而使得计算和放样轻松、方便、快捷。

【计算范围】坐标部分包括：各等级公路和高速路的直线段、圆曲线段、加对称缓和曲线的圆曲线段中、边桩坐标及放样计算。

高程部分包括：各等级公路和高速路的直线段、竖曲线段、缓和超高段中、边桩高程及放样计算。

一、坐标及放样程序【起算数据】：交点里程桩号、交点坐标、前直线方位角、交点转角、圆曲线半径、缓和曲线长（一）主程序：XYFY2Fix4M"Xc=":P"Yc=" ＇如须放样输入测站坐标值,否则输入0Lbl0{HS}:H"K0+0":S"B=(-B,0,B)":Prog"A" ＇变量输入和声明S=0=>E=0:≠>E=90 ＇求边桩坐标时的偏角值判定L=πRN÷180+V ＇含有缓曲线的曲线总长T=(V÷2-V^3÷(240R^2))+(R+(V^2÷(24R)-V^4÷(2688R^3)))tan(N÷2)＇切线长A=Q-T:B=A+V:D=A+L:C=D-V ＇A-ZH(ZY)点桩号,B-HY点桩号,C-YH点桩号,D-HZ(YZ)点桩号Rec(T,F+180) ＇求ZH(ZY)点的坐标增量Z[1]=W+I:Z[2]=K+J ＇ZH(ZY)点坐标值Rec(T,F+GN) ＇求HZ(YZ)点的坐标增量Z[3]=W+I:Z[4]=K+J ＇HZ(YZ)点的坐标值Lbl1U=S÷(Abs(S+10E-19)):S=AbsS ＇左右偏角每件判定H<A=>Goto2: ≠> ＇第一直线上坐标计算判定H<B=>Goto3: ≠> ＇第一缓和曲线上坐标计算判定H<C=>Goto4: ≠> ＇净圆曲线上坐标计算判定H<D=>Goto5: ≠> Goto6⊿⊿⊿⊿＇第二缓和曲线和第二直线上坐标计算判定Lbl2Rec(Q-H,F+180) ＇第一直线上中桩坐标增量计算X=W+I:Y=K+J ＇中桩坐标值计算Rec(S,F+180-(180-E)U) ＇第一直线上与中桩同一断面的边桩增量计算X=X+I:Y=Y+J ＇边桩坐标值计算S=SU ＇边长条件判定Goto8Lbl3＇进入第一缓曲线段计算Z=H-A ＇待求点至ZH点的距离O=90Z^2÷(πRV) ＇所求点缓曲线对应的圆心角X=Z-Z^5÷(40R^2V^2)+Z^9÷(3456R^4V^4) ＇缓曲线上ZH点至待求的支距X坐标Z=Z^3÷(6RV)-Z^7÷(336R^3V^3)+Z^11÷(42240R^5V^5) ＇缓曲线上ZH点至待求点的支距Y坐标Lbl7Rec(X,F)X=Z[1]+I:Y=Z[2]+JRec(Z,F+90G) ＇把中桩支距坐标转换成测量统一坐标X=X+I:Y=Y+JRec(S,F+O G+EU) ＇曲线边桩坐标增量计算X=X+I:Y=Y+J ＇坐标计算结果S=SU ＇边长条件Goto8Lbl4＇进入主圆计算部分Z=H-A-V÷2 ＇曲线上待求点到HY点的里程O=180Z÷(Rπ) ＇待求点方位角所对应的圆心角X=Rsin O+(V÷2-V^3÷(240R^2))Z=R(1-cos O)+(V^2÷(24R)-V^4÷(2688R^3))＇主圆部分待求点支距坐标XY Goto7Lbl5＇进入第二缓和曲线计算Z=D-H ＇曲线上待求点到HZ点的里程O=90Z^2÷(πRV) ＇待求点方位角所对应的圆心角X=Z-Z^5÷(40Z R^2V^2)+^9÷(3456R^4V^4)Z=Z^3÷(6RV)-Z^7÷(336R^3V^3)+Z^11÷(42240R^5V^5) ＇待求点支距坐标XY Rec(X,F+GN+180)X=Z[3]+I:Y=Z[4]+J ＇把中桩支距坐标转换成测量统一坐标Rec(Z,F+GN+180-90G)X=X+I:Y=Y+JRec(S,F+GN+180-O G-(180-E)U) ＇曲线边桩坐标增量计算X=X+I:Y=Y+J ＇坐标计算结果S=SUGoto8Lbl6＇进入第二直线段计算部分Rec(H-D+T,F+GN) ＇中桩坐标增量计算X=W+I:Y=K+J ＇中桩坐标计算结果Rec(S,F+GN+EU) ＇边桩坐标增量计算X=X+I:Y=Y+J ＇边桩坐标计算结果S=SUGoto8Lbl8X=X◢＇显示坐标计算结果Y=Y◢M≠0=>Goto9: ≠>Goto0 ⊿Lbl9 ＇进入放样部分计算Pol(X-M,Y-P) ＇增量计算J<0=>J=J+360: ≠>J=J ⊿J"A"=Int J+0.01Int(60Frac J)+0.006Frac(60Frac J) ◢＇方位角计算结果（此显示值为度分秒格式）I"L"=I◢＇极距计算结果Goto0(二)子程序:AH≤ZY2(ZH2)=>Q=JD1:W=X1:K=Y1:R=R1:F=F0:N=N1:V=Ls1:G=1(-1):=>H≤ZY3(ZH3)=>Q=JD2:W=X2:K=Y2:R=R2:F=F1:N=N2:V=Ls2:G=1(-1):=>H≤ZY4(ZH4)=>Q=JD3:W=X3:K=Y3:R=R3:F=F2:N=N3:V=Ls3:G=1(-1):=>H≤ZY5(ZH5)=>Q=JD4:W=X4:K=Y4:R=R4:F=F3:N=N4:V=Ls4:G=1(-1):=>……………………………H≤ZY n+1(ZH n+1)=>Q=JD n:W=X n:K=Y n:R=R n:F=F n-1:N=N n:V=Ls n:G=1(-1)注:H—待求点桩号ZY n+1(ZH n+1)—后曲线起点桩号 Q—本曲线交点桩号 W—本曲线交点X坐标K—本曲线交点Y坐标 R—本曲线内圆曲线半径 F—前直线方位角 N—本交点转角V—本曲线内缓和曲线长 G—线路转向符,左转角-1,右转角为+1二、高程及放样程序【起算数据】：变坡点桩号A、变坡点高程B、前后纵坡C D、竖曲半径R、横坡度F、最大超高横坡度G、ZH点桩号、YH点桩号、HZ点桩号、施工层厚度N、视线高W （一）主程序：GCFYW"SXG="：N＇常量输入（视线高程、施工层厚度）Lbl 1｛LО｝:L"K0+0"：О"B=" :Prog"B"＇变量输入及声明K=C-D ＇坡度差及竖曲线偏角K>0=>J=-1:≠>J=1⊿＇凹凸竖曲线判断条件T=Abs(RK)÷2 ＇竖曲线切线长,如需显示则在后面加上◢H=A-T ＇竖曲线起点ZY桩号,如需显示则在后面加上◢X=A+T ＇竖曲线终点YZ桩号,如需显示则在后面加上◢M=L-A ＇待求点至变坡点的距离L≤H=>P=B+MC-N-W: ≠> ＇前纵坡直线段中桩高程计算L≤A=>P=B+MC+J(T-AbsM)2÷(2R)-N-W: ≠>＇前纵坡竖曲线中桩高程计算L≤X=>P=B+MD+J(T-AbsM)2÷(2R)-N-W: ≠>＇后纵坡竖曲线中桩高程计算L>X=>P=B+MD-N-W ⊿⊿⊿⊿⊿＇后纵坡直线段中桩高程计算E=0=>Goto 2 ＇直线部分中、边桩计算判断条件≠>L≤E=> Goto 2 ＇直线部分中、边桩计算判断条件≠>L>V=> Goto 2 ＇直线部分中、边桩计算判断条件≠>L>E=> Goto 3 ＇前缓曲超高部分中、边桩计算判断条件≠> L>U=> Goto 4⊿⊿⊿⊿＇后缓曲超高部分中、边桩计算判断条件Lbl 2O=0=>Z“H”=AbsP◢＇显示直线部分中、边桩高程计算结果≠>Y“HB”=Abs(P-OF)◢不需放样显示高程,否则显示放样水准尺计数⊿Goto 1Lbl 3S=（L-E）（F+G）÷（V-U）-F＇前缓曲线超高横坡度计算S≤F=> Goto 5：≠>S≤G=> Goto 6 ＇超高起点(终点)至零界面、零界面至全超高及全超高≠> Goto 7 ⊿⊿部分的计算判断条件Lbl 4S=（V-L）（F+G）÷（V-U）-F＇后缓曲线超高横坡度计算S≤F=> Goto 5：≠>S≤G=> Goto 6 ＇超高起点(终点)至超高变坡零界面、零界面至全≠> Goto 7 ⊿⊿超高、全超高部分计算判断条件Lbl 5O=0=>Z“H1”=AbsP◢＇超高起点(终点)至超高变坡零界面中、边桩高程≠>Y“H-D1”=Abs(P-OF)◢计算结果显示Q“H-G1”=Abs(P+OS)◢⊿Goto 1Lbl 6O=0=>Z“H2”=AbsP◢＇超高变坡零界面至全超高中、边桩高程≠>Y“H-D2”=Abs(P-OS)◢计算结果显示Q“H-G2”=Abs(P+OS)◢⊿Goto 1Lbl 7O=0=>Z“H3”=AbsP◢＇全超高中、边桩高程计算结果显示≠>Y“H-D3”=Abs(P-OG)◢Q“H-G3”=Abs(P+OG)◢⊿Goto 1注：W“SXG”-视线高N-施工层厚度L"K0+0"-待求点桩号О"B="-中桩至边桩距离（二）子程序：BL≤ZY2=>A=JD1:B=H1:C=I0:D=I1: R=R1:F=I h:E=ZH1:U=YH1:V=HZ1:G=I max1 :≠>L≤ZY3=>A=JD2:B=H2:C=I1:D=I2: R=R2:F=I h:E=ZH2:U=YH2:V=HZ2:G=I max2 :≠>L≤ZY4=>A=JD3:B=H3:C=I2:D=I3: R=R3:F=I h:E=ZH3:U=YH3:V=HZ3:G=I max3 :≠>L≤ZY4=>A=JD3:B=H3:C=I2:D=I3: R=R3:F=I h:E=ZH3:U=YH3:V=HZ3:G=I max3 :≠> ……………………L≤ZY n+1=>A=JD n:B=H n:C=I n-1:D=I n: R=R n:F=I h:E=ZH n:U=YH n:V=HZ n:G=I maxn注:ZY N+1:-下一个竖曲线起点桩号（如：变坡点2的计算范围为前竖曲线终点YZ1至后竖曲线起点ZY3，则在程序中“L≤ZY3”就应输入ZY3的里程桩。

公路曲线中桩边桩坐标fx4800P计算程序

公路曲线中桩边桩坐标fx4800P计算程序在快速发展的公路建设事业中,使用GPS坐标控制系统的高等级公路日渐增多,并且,随着高精度的测量设备的广泛采用,其他等级公路也开始普遍应用坐标控制。

使用坐标控制的公路项目,有着线型流畅、行车速度快而舒适、环境优美等特点,这些都来源于高质量的设计与施工,卡西欧fx-4800p计算器由于其轻便、精度高、易于操作,能够解决设计、施工当中很多的实际困难,因此,倍受公路工程技术人员的青睐。

本人通过在高速公路多年的施工经验,编写了许多公路测量实用程序,如常用的中桩边桩计算程序,匝道坐标计算程序,施工放样程序,换带程序等等。

下面把公路中常用的中桩边桩计算程序介绍给大家,希望本程序能够对公路测量同行有所帮助.。

1 程序清单:LbI 0:Defm 6↙″CURVE MODE″:{M}:M″0:T1=T2,1:T1≠T2″↙M=0=>Goto 7:≠>M=1=>Goto 8:≠>Goto 0⊿⊿↙LbI 7:{ABCDEFGR}↙A″ZHD″:B″HZD″:C″XJD″:D″YJD″:E″ELPHA″:F″PJ:Z(-),Y(+)″:G″L0″:R↙T=(R+G2÷(24R)-G ÷(2688R ))tan(AbsF÷2)+0.5G-G ÷(240R2) ↙U=G:V=T: Goto 1↙LbI 8:{ABCDEFGUTVR}↙A″ZHD″:B″HZD″:C″XJD″:D″YJD″:E″ELPHA″:F″PJ:Z(-),Y(+)″:G″L01″:U″L0 2″:T″T1″:V″T2″: R↙LbI 1:{H}:H″DQD″↙H≥B=>M=1: ≠>H>B-U=>M=2: ≠>H≥A+G=>M=3: ≠>H≥A=>M=4: ≠>HM=5⊿⊿⊿⊿⊿↙M≤2=>L=Abs(H-B): ≠>L=Abs(H-A)⊿↙M=1=>Z[1]=L: ≠>M=2=>Z[1]=√((L-L ÷(40R2U2)+L ÷(3456R U ))2+(L ÷(6RU)- L ÷(336R U )+ L÷(42240R U ))2): ≠>M=3=>Z[1]=√((G-G ÷(40R2)+G ÷(3456R ))2+(G2÷(6R)-G ÷(336R )+G ÷(42240R ))2):≠>M=4=>Z[1]=√((L-L ÷(40R2G2)+L÷(3456R G ))2+ (L ÷(6RG)- L ÷(336R G )+ L÷(42240R G ))2): ≠>M=5=>Z[1]=L⊿⊿⊿⊿⊿↙M≤2=>Z[2]=C+Vcos(E+F): Z[3]=D+Vsin(E+F): ≠>Z[2]=C+Tcos(E+180): Z[3]=D+Tsin(E+180) ⊿↙M=1=>Z[4]=E+F:≠>M=2=>Z[4]=E+F+180:≠>M=5=>Z[4]=E+ 180:≠>Z[4]=E ⊿⊿⊿↙N=F÷AbsF↙M=1=>I=0:≠>M=2=>I=-N(30L2÷(ΠRU)):≠>M=3=>I=N(30G÷(ΠR)):≠>M=4=>I=N(30L2÷(ΠR G)):≠>M=5=>I=0⊿⊿⊿⊿⊿↙Z[4]=Z[4]+I↙Z[5]=Z[2]+Z[1]cosZ[4]↙Z[6]=Z[3]+Z[1]sinZ[4]↙M≠3=>Goto 3⊿↙LbI 2:P=90L÷(ΠR):Q=90(L-G)÷(ΠR)↙Z[5]= Z[5]+2RsinQcos(E+NP)↙Z[6]= Z[6]+2RsinQsin(E+NP)↙LbI 3↙X″X(ZX)″=Z[5]◢Y″Y(ZX)″=Z[6]◢LbI 4↙M=3=>I=90N(2L-G) ÷(ΠR): ≠>I=3I⊿↙M=1=>P=E+F: ≠>M=2=>P=E+F+180: ≠>M=5=>P=E+180: ≠>P=E⊿⊿⊿↙P=P+I-90:Q=P+180↙M=2=>X=P-90: ≠>M=5=>X=P+180: ≠>X=Q-90⊿⊿↙XX=X+360: ≠>X>360=>X=X-360⊿⊿↙X″QXFWJ(Z-Y)″=X◢{SW}:S″SZ ″: W″SY ″:AbsS+AbsW=0=>Goto 1⊿↙M=2=>Goto 5: ≠>M=5=>Goto 5: ≠>Goto 6⊿⊿↙LbI 5:S=-S:W=-W↙LbI 6↙X″X(Z)″=Z[5]+ScosP◢Y″Y(Z)″=Z[6]+SsinP◢X″X(Y)″=Z[5]+WcosQ◢Y″Y(Y)″=Z[6]+WsinQ◢Goto 1↙2 程序参数说明:开始运行程序时选择曲线类型(CURVE MODE),0为等切线曲线,1为不等长切线曲线。

CASIOfx-4800P 路线计算(偏角法对称型)

CASIO fx-4800P 路线放样计算程序(偏角法对称型计算)程序说明：1、本程序是计算对称型圆曲线中、边桩坐标。

2、程序所占内存2400步左右。

3、本程序单个交点输入。

4、程序运行计算时有五大桩号提示。

L1 DegL2 Defm 19L3 Fixm 程序菜单1-设站，2-数据输入，3-路线坐标计算，4-设新点L4 Lbl 0L5 {A}:A”SZ-1,INPUT.SJ-2,LXFY-3,XD-4”L6 A=1=> Goto 1L7 A=2=> Goto 2L8 A=3=> Goto 3L9 A=4=> Goto BL10 Goto 0 输入站点与后视点坐标L11 Lbl 1L12 {BCDE}:B”ZDX”:C”ZDY”:D”HDX”:E:”HDY”L13 Pol(D-B,E-C)L14 Z[1]”HD”=I◢L15 J>0=>Z[2]”HA”=J◢≠>Z[2]”HA”=360+J◢L13 Goto 0 交点要素数据输入L14 Lbl 2L15 {FGHIKLMNRP}:F”JD1X”:G”JD2Y”:H”JDX”:K”JDY”:L”JD2X”:M”JD2Y”:N”JDLC”:R”JDR”:P”JDLS1”L16 Pol(H-F,K-G)L17 J>0=>Z[3]”1=>0”=J◢≠>Z[3]”1=>0”=360+J◢L18 Pol(L-H,M-K)L19 J>0=>Z[4]”0=>2”=J◢≠>Z[4]”0=>2”=360+J◢Z[3]+180>360=>Z[18]=Z[3]-180 计算交点转角和判断是左角还是右角Z[3]+180<360=>Z[18]=Z[3]+180Z[18]<Z[4]=>Z[19]=Abs(Z[18]+360-Z[4]])Z[18]>Z[4]=>Z[19]=Abs(Z[18]-Z[4])Z[19]>180=>U=1:O=Z[19]-180Z[19]<180=>U=-1:O=180-Z[19]L20 Z[5]=90P/R/πL21 Z[6]=(R+P^2/24/R)Tan(O/2)+(P/2-P^3/240/R^2) 切线长L22 Z[7]= πON/180+P 曲线长L23 Z[8]= πR(O-2Z[5])/180 圆曲线长L24 Z[13]=N-Z[6] ZH点L25 Z[14]=Z[13]+P HY点L26 Z[15]=Z[14]+Z[8]/2 QZ点L27 Z[16]=Z[14]+Z[8] YH点L28 Z[17]=Z[16]+P HZ点L29 Goto 0 路线计算输入里程桩号L30 Lbl 3L31 {S}:S”INPUT LC”L32 S<Z[13]+10=>S>Z[13]-10=>Goto C ZH点L33 P≠0=>S<Z[14]+10=>S>Z[14]-10=>Goto D YH点L34 S<Z[15]+10=>S>Z[15]-10=>Goto E QZ点L35 P≠0=>S<Z[16]+10=>S>Z[16]-10=>Goto F HY点L36 S<Z[17]+10=>S>Z[17]-10=>Goto G HZ点L37 S>Z[13]=>S≤Z[14]=>T=S-Z[13]:Goto 4 ZH点至HY点段L38 S>Z[14]=>S<Z[16]=>T=S-Z[14]:Goto 6 HY点至YH点段L39 S≥Z[16]=>S<Z[17]=>T=Z[17]-S: Goto 4 YH点至ZH点段L40 S≤Z[13]=>Goto 8 ZH点前直线段L41 S≥Z[17]=>Goto 9 HZ点后直线段L42 Goto 0 计算缓和段长偏角L43 Lbl 4L44 P≠0=>Z[9]=30T^2/π/R/PL45 P≠0=>Z[10]=√((T-T^5/(40R^2P^2))^2+(T^3/6RP)^2)L46 P≠0=>S≥Z[13]=>S≤Z[14]=>Goto 5 ZH点至HY点段L47 P≠0=>S≥Z[16]=>S≤Z[17]=>Goto 7 YH点至HZ点段L48 S>Z[14]=>S<Z[16]=>T=S-Z[14]:Goto 6 HY点至YH点段L49 S≤Z[13]=>Goto 8L50 S≥Z[17]=>Goto 9L48 Goto 0 计算ZH点至YH点段中线坐标L51 Lbl 5L52 U=1=>V=Z[3]-Z[9]:Z[11]”H1”=Z[3]-3Z[9] ◢U=-1=>V=Z[3]+Z[9]:Z[11]”H1”=Z[3]+3Z[9] ◢L53 Rec(Z[10],V)L54 X=H+Z[6]COS(180+Z[3])+IL55 Y=K+Z[6]SIN(180+Z[3])+JL56 Goto A 计算圆曲线段中线坐标L57 Lbl 6L58 Z[9]=90T/π/RL59 U=1=>V=Z[3]-Z[5]-Z[9]:Z[12]=Z[3]-Z[5]/3:Z[11]”YQ”=Z[3]-Z[5]-2Z[9] ◢U=-1=>V=Z[3]+Z[5]+Z[9]:Z[12]=Z[3]+Z[5]/3:Z[11]”YQ”=Z[3]+Z[5]+2Z[9] ◢L60 Rec(ABS(2RSinZ[9]),V)L61 X=H+Z[6]Cos(180+Z[3])+ √((P-P^3/40/R^2)^2+(P^2/6/R)^2)CosZ[12]+IL62 Y=K+Z[6]Sin(180+Z[3])+ √((P-P^3/40/R^2)^2+(P^2/6/R)^2)SinZ[12]+JL63 Goto A 计算YH点至HZ点段坐标L64 Lbl 7L65 U=1=>V=Z[4]+180+Z[9]:Z[11]”H2”=Z[4]+3Z[9] ◢U=-1=>V=Z[4]+180-Z[9]:Z[11]”H2”=Z[4]-3Z[9] ◢L66 Rec(Z[10],V)L67 X=H+Z[6]CosZ[4]+IL68 Y=K+Z[6]SinZ[4]+JL69 Goto A 计算ZH点前直线坐标L70 Lbl 8L71 Z[11]=Z[3]L72 Rec(N-S,Z[3])L73 X=H-IL74 Y=K-JL75 Goto A 计算HZ点后直线坐标L76 Lbl 9L77 Z[11]=Z[4]L78 Rec(S-N+2Z[6]-Z[7],Z[4])L79 X=H+IL80 Y=K+jL81 Goto A 计算边桩坐标与站点到计算点的平距和方位角L82 Lbl AL83 {W}:W”INPUT FX”L84 X”X”=X+WCos(90+Z[11]) ◢L85 Y”Y”=Y+WSin(90+Z[11]) ◢L86 Pol(X-B,Y-C)L87 I”HD”=I◢L88 J>0=>J”HA”=J◢≠>J”HA”=360+J◢L89 {A}:A”Cont.FY-3,EXT-4”L90 A=3=>Goto 3L91 A=4=>Goto 0 计算新点的坐标L92 Lbl bL93 {O}:O”XD.HD”:{Z}:Z”XD.HA”L94 Rec(O,Z)L95 X”XP”=B+I◢L96 Y”YP”=C+J◢L97 Goto 0 ZH点L98 Lbl CL99 P=0=>Z[13]”ZYD”=Z[13] ◢P≠0=>Z[13]”ZHD”=Z[13] ◢L100 {A}:A”WuDZ-1,NO-3”L101 A=1=>S=Z[13]: T=S-Z[13]:Goto 4L102 A=3⇒S<Z[13]⇒Goto 8≠> T=S-Z[13]:Goto 4L103 Lbl D YH点L104 Z[14]”HYD”=Z[14] ◢L105 {A}:A”WuDZ-1,NO-3”L106 A=1=>S=Z[14]:T=S-Z[13]:Goto 4L107 A=3=> T=S-Z[13]:Goto 4L108 Lbl E QZ点L109 Z[15]”QZD”=Z[15] ◢L110 {A}:A”WuDZ-1,NO-3”L111 A=1=>S=Z[15]:T=S-Z[14]:Goto 6L112 A=3=> T=S-Z[14]:Goto 6L113 Lbl F YH点L114 Z[16]”YHD”=Z[16] ◢L115 {A}:A”WuDZ-1,NO-3”L116 A=1=>S=Z[16]:T=Z[17]-S:Goto 4L117 A=3=> T=Z[17]-S:Goto 4L118 Lbl G HZ点L119 P=0=>Z[17]”YZD”=Z[17] ◢P≠0=>Z[17]”HZD”=Z[17] ◢L120 {A}:A”WuDZ-1,NO-3”L121 A=1=>S=Z[17]:T=Z[17]-S:Goto 4L122 A=3=>T=Z[17]-S:Goto 4。

Casio fx4800任意曲线程序

Casio fx-4800(4850)计算器程序任意曲线坐标计算程序(HHQX)LbI 0 ：｛QD｝←┘G=I=>Goto 1ΔG=-1=>Goto 2△G=0=>Goto 4 △←┘LbI 1 ：L=C÷R ：M=L+Q-K ：F=A+90(M2-L2)÷（πC）：F◢T=A-90L2/πC：Goto 3△←┘LbI 2 ：L=C÷R ：M=L+K-Q ：:T=A+90L2÷πC+180 ：F=A-90(M2-L2)÷πC:F◢Goto 3△←┘LbI 3 ：N=(M-L)-(M^5-L^5)÷40C2+(M^9-L^9)÷3456C^4-(M^13-L^13)÷599040C^6：E=(M^3-L^3)/6C-(M^7-L^7)÷336C^3+(M^11-L^11)÷42240C^5-(M^15-L^15)÷9676800 C^7 ：X=I+NcosT-GEsinT+Dcos(F-90)◢Y=J+NsinT+GEcosT+Dsin(F-90)◢Goto 0←┘LbI 4 ：L=Q-K ：N=Rsin(180L÷πR)：E=R(1-cos(180L÷πR))：F=A+180L÷πR：F◢X=I+NcosA-EsinA+Dcos(F-90)◢Y=J+NsinA+EcosA+Dsin(F-90)◢Goto 0说明:当G=1时,计算第一缓和曲线或由大圆到小圆的卵型曲线；G=-1时,计算第二缓和曲线或由小圆到大圆的卵型曲线；G=0时,计算直线或圆曲线。

R为起算点半径,线路左转时R、C取负值，右转时R、C取正值，计算第一缓和曲线时R取10^45。

计算第二缓和曲线时R取圆曲线半径，计算直线时R取10^45.此计算思路是沿路线的前进方向依次进行分段计算。

卡西欧4800P计算器大地坐标放样程序

计算点K
LX
置仪点
置方位角预
后视点
eh曲线外距p输入1时计算测点的左边点输入0时计算中点
主程序 WTO
(卡西欧 4800P 计算器用坐标放样程序)
L“L1”E“L2”A“ZH”：Z：F：U“JX”V“JY”M“X0”N“Y0”：J=180/πR：B“HZ”=A+Z/J+L/2+E/2▲ W“-1，1”：Prog“K”↙ LIb0：｛K｝：｛H｝：｛P｝：K≤A ＞Goto1△K≤A+L ＞Goto2△K＜B-E ＞Goto3△ K≤B ＞Goto4△Goto5△↙ LIb1：S=K-A：X=U-TcosF+ScosF+Hcos（F+90P）▲Y=V-TsinF+SsinF+Hsin（F+90P）▲Goto6↙ LIb2：S=K-A：Prog“A”：Goto6↙ LIb3：S=K-A-L：Prog“B”：Goto6↙ LIb4：S=B-K：Prog“C”：Goto6↙ LIb5：S=K-B：X=U+(I+S)cos(F+WZ)+Hcos(F+WZ+90P)▲ Y=V+(I+S)sin(F+WZ)+Hsin(F+WZ+90P)▲Goto6↙ LIb6：X≠M ＞Prog“G”↙ Goto0 子程序G：
L1－第一缓和曲线总长；L2－第二缓和曲线总长；Z－曲线总转角；F－ZH点至JD点的方位角：
-1 2 缓点桩号；HZ－缓直点桩号；JX－交点的横坐标(N)；JY－交点的纵坐标(E)； X0－置仪点的横坐标(N)；ZY－置仪点的纵坐标(E)；R－曲线半径：HZ－缓直点桩号； “-1，1?"—曲线左转时输如“-1”，曲线右转时输入“+1”；Eh——曲线外距 K－测点的桩号；H－测点边桩与中线之间的距离； P－输入“-1”时计算测点的左边点，输入“0”时计算中点；输入“+1”时计算右边点 X－测点的点的横坐标(N)；Y－测点的纵坐标(E)； LX－测距；C－置仪点与测点间的方位角(即拨角)。注：应用此程序时须预置置仪点至后视点间的方位角。

Casio4800平曲线坐标计算教学资料

C a s i o4800平曲线坐标计算Casio fx-4800P平曲线坐标计算程序湖南新光工程有限公司三板溪项目部李甲生程序目的：依平曲线要素计算直线、圆曲线、缓和曲线的任意中桩、左、右桩坐标。

程序说明：K0：起始桩号X0：起始X坐标Y0：起始Y坐标ALF：起始方位角R：半径LS：缓和曲线长N：曲线左转N=1，右转N=2 K：待求桩号LL、LR：左、右桩距离Q：左、右桩与中线斜交角求得XZ、YZ、XL、YL、XR、YR分别为中桩、左、右桩坐标。

一、直线段文件名：ZX (COMP)程式：L”K0”:O”X0”:P”Y0”:W”ALF”:Lbl 0:{K}:X”XZ”=O+(K-L)cosW◢Y”YZ”=P+(K-L)sinW◢{B}:S”XL”=X-B”LL”cos(W+Q)◢T”YL”=Y-Bsin(W+Q) ◢{C}:U”XR”=X+C”LR”cos(W+Q)◢V”YR”=Y+Csin(W+Q)◢Goto 0注：在程序执行过程中，赋给的要素变数的值被固定不变，可对变数（K、LL、LR）赋予不同值，迅速求得所需坐标。

二、圆曲线段文件名：YQX (COMP)程式：L”K0”:O”X0”:P”Y0”:W”ALF”:Lbl 1:{K}:J=(-1)^N (K-L)÷R×180÷π:D=2Rsin((-1)^N J÷2):X”XZ”=O+Dcos(W+J÷2)◢Y”YZ”=P+Dsin(W+J÷2)◢{B}:S”XL”=X-B”LL”cos(W+J+Q)◢T”YL”=Y-Bsin(W+J+Q) ◢{C}:U”XR”=X+C”LR”cos(W+J+Q)◢V”YR”=Y+Csin(W+J+Q) ◢Goto 1注：若没有直接HY点方位角，则ALFHY=ALFZH±Ls/2/R×180/π,(左转-，右转+)。

三、缓和曲线文件名：HHQX (COMP)程式：L”K0”:O”X0”:P”Y0”:W”ALF”:M”LS”:Lbl3:{K}:I=(-1)^N×(K-L)^2÷M÷R÷6×180÷π:D=(K-L)-(K-L)^5÷90÷(RM)^2:X”XZ”=O+Dcos(W+I)◢Y”YZ”=P+Dsin(W+I)◢{B}:S”XL”=X-B”LL”cos(W+3I+Q)◢T”YL”=Y-Bsin(W+3I+Q) ◢{C}:U”XR”=X+C”LR”cos(W+3I+Q)◢V”YR”=Y+Csin(W+3I+Q)◢Goto 3注：1、坐标计算方法是根据偏角法原理；2、缓和曲线（ZH~HY或YH~HZ）以ZH（或HZ）为起始点；3、平曲线左转（ZH~HY段N=1，YH~HZ段N=2），曲线右转（ZH~HY段N=2，YH~HZ段N=1）。

CASIO-fx4800P型计算器连续计算线路高程、坐标

【计算范围】坐标部分包括：各等级公路和高速路的直线段、圆曲线段、加对称缓和曲线的圆曲线段中、边桩坐标及放样计算。

高程部分包括：各等级公路和高速路的直线段、竖曲线段、缓和超高段中、边桩高程及放样计算。

一、坐标及放样程序【起算数据】：交点里程桩号、交点坐标、前直线方位角、交点转角、圆曲线半径、缓和曲线长（一）主程序：XYFY2Fix4M"Xc=":P"Yc=" ＇如须放样输入测站坐标值,否则输入0Lbl0{HS}:H"K0+0":S"B=(-B,0,B)":Prog"A" ＇变量输入和声明S=0=>E=0:≠>E=90 ＇求边桩坐标时的偏角值判定L=πRN÷180+V ＇含有缓曲线的曲线总长T=(V÷2-V^3÷(240R^2))+(R+(V^2÷(24R)-V^4÷(2688R^3)))tan(N÷2)  ＇切线长A=Q-T:B=A+V:D=A+L:C=D-V ＇A-ZH(ZY)点桩号,B-HY点桩号,C-YH点桩号,D-HZ(YZ)点桩号Rec(T,F+180) ＇求ZH(ZY)点的坐标增量Z[1]=W+I:Z[2]=K+J ＇ZH(ZY)点坐标值Rec(T,F+GN) ＇求HZ(YZ)点的坐标增量Z[3]=W+I:Z[4]=K+J ＇HZ(YZ)点的坐标值Lbl1U=S÷(Abs(S+10E-19)):S=AbsS ＇左右偏角每件判定H<A=>Goto2: ≠> ＇第一直线上坐标计算判定H<B=>Goto3: ≠> ＇第一缓和曲线上坐标计算判定H<C=>Goto4: ≠> ＇净圆曲线上坐标计算判定H<D=>Goto5: ≠> Goto6⊿⊿⊿⊿ ＇第二缓和曲线和第二直线上坐标计算判定Lbl2Rec(Q-H,F+180) ＇第一直线上中桩坐标增量计算X=W+I:Y=K+J ＇中桩坐标值计算Rec(S,F+180-(180-E)U) ＇第一直线上与中桩同一断面的边桩增量计算X=X+I:Y=Y+J ＇边桩坐标值计算S=SU ＇边长条件判定Goto8Lbl3 ＇进入第一缓曲线段计算Z=H-A ＇待求点至ZH点的距离O=90Z^2÷(πRV) ＇所求点缓曲线对应的圆心角X=Z-Z^5÷(40R^2V^2)+Z^9÷(3456R^4V^4) ＇缓曲线上ZH点至待求的支距X坐标Z=Z^3÷(6RV)-Z^7÷(336R^3V^3)+Z^11÷(42240R^5V^5) ＇缓曲线上ZH点至待求点的支距Y坐标Lbl7Rec(X,F)X=Z[1]+I:Y=Z[2]+JRec(Z,F+90G) ＇把中桩支距坐标转换成测量统一坐标X=X+I:Y=Y+JRec(S,F+O G+EU) ＇曲线边桩坐标增量计算X=X+I:Y=Y+J ＇坐标计算结果S=SU ＇边长条件Goto8Lbl4 ＇进入主圆计算部分Z=H-A-V÷2 ＇曲线上待求点到HY点的里程O=180Z÷(Rπ) ＇待求点方位角所对应的圆心角X=Rsin O+(V÷2-V^3÷(240R^2))Z=R(1-cos O)+(V^2÷(24R)-V^4÷(2688R^3)) ＇主圆部分待求点支距坐标XY Goto7Lbl5 ＇进入第二缓和曲线计算Z=D-H ＇曲线上待求点到HZ点的里程O=90Z^2÷(πRV) ＇待求点方位角所对应的圆心角X=Z-Z^5÷(40Z R^2V^2)+^9÷(3456R^4V^4)Z=Z^3÷(6RV)-Z^7÷(336R^3V^3)+Z^11÷(42240R^5V^5) ＇待求点支距坐标XY Rec(X,F+GN+180)X=Z[3]+I:Y=Z[4]+J ＇把中桩支距坐标转换成测量统一坐标Rec(Z,F+GN+180-90G)X=X+I:Y=Y+JRec(S,F+GN+180-O G-(180-E)U) ＇曲线边桩坐标增量计算X=X+I:Y=Y+J ＇坐标计算结果S=SUGoto8Lbl6 ＇进入第二直线段计算部分Rec(H-D+T,F+GN) ＇中桩坐标增量计算X=W+I:Y=K+J ＇中桩坐标计算结果Rec(S,F+GN+EU) ＇边桩坐标增量计算X=X+I:Y=Y+J ＇边桩坐标计算结果S=SUGoto8Lbl8X=X◢ ＇显示坐标计算结果Y=Y◢M≠0=>Goto9: ≠>Goto0 ⊿Lbl9 ＇进入放样部分计算Pol(X-M,Y-P) ＇增量计算J<0=>J=J+360: ≠>J=J ⊿J"A"=Int J+0.01Int(60Frac J)+0.006Frac(60Frac J) ◢ ＇方位角计算结果（此显示值为度分秒格式）I"L"=I◢ ＇极距计算结果Goto0∙(二)子程序:AH≤ZY2(ZH2)=>Q=JD1:W=X1:K=Y1:R=R1:F=F0:N=N1:V=Ls1:G=1(-1):=>H≤ZY3(ZH3)=>Q=JD2:W=X2:K=Y2:R=R2:F=F1:N=N2:V=Ls2:G=1(-1):=>H≤ZY4(ZH4)=>Q=JD3:W=X3:K=Y3:R=R3:F=F2:N=N3:V=Ls3:G=1(-1):=>H≤ZY5(ZH5)=>Q=JD4:W=X4:K=Y4:R=R4:F=F3:N=N4:V=Ls4:G=1(-1):=>……………………………H≤ZY n+1(ZH n+1)=>Q=JD n:W=X n:K=Y n:R=R n:F=F n-1:N=N n:V=Ls n:G=1(-1)注:H—待求点桩号ZY n+1(ZH n+1)—后曲线起点桩号 Q—本曲线交点桩号 W—本曲线交点X坐标K—本曲线交点Y坐标 R—本曲线内圆曲线半径 F—前直线方位角 N—本交点转角V—本曲线内缓和曲线长 G—线路转向符,左转角-1,右转角为+1二、高程及放样程序【起算数据】：变坡点桩号A、变坡点高程B、前后纵坡C D、竖曲半径R、横坡度F、最大超高横坡度G、ZH点桩号、YH点桩号、HZ点桩号、施工层厚度N、视线高W （一）主程序：GCFYW"SXG="：N ＇常量输入（视线高程、施工层厚度）Lbl 1｛LО｝:L"K0+0"：О"B=" :Prog"B" ＇变量输入及声明K=C-D ＇坡度差及竖曲线偏角K>0=>J=-1:≠>J=1⊿ ＇凹凸竖曲线判断条件T=Abs(RK)÷2 ＇竖曲线切线长,如需显示则在后面加上◢H=A-T ＇竖曲线起点ZY桩号,如需显示则在后面加上◢X=A+T ＇竖曲线终点YZ桩号,如需显示则在后面加上◢M=L-A ＇待求点至变坡点的距离L≤H=>P=B+MC-N-W: ≠> ＇前纵坡直线段中桩高程计算L≤A=>P=B+MC+J(T-AbsM)2÷(2R)-N-W: ≠>  ＇前纵坡竖曲线中桩高程计算L≤X=>P=B+MD+J(T-AbsM)2÷(2R)-N-W: ≠> ＇后纵坡竖曲线中桩高程计算L>X=>P=B+MD-N-W⊿⊿⊿⊿⊿ ＇后纵坡直线段中桩高程计算E=0=>Goto 2 ＇直线部分中、边桩计算判断条件≠>L≤E=> Goto 2 ＇直线部分中、边桩计算判断条件≠>L>V=> Goto 2 ＇直线部分中、边桩计算判断条件≠>L>E=> Goto 3 ＇前缓曲超高部分中、边桩计算判断条件≠> L>U=> Goto 4⊿⊿⊿⊿ ＇后缓曲超高部分中、边桩计算判断条件Lbl 2O=0=>Z“H”=AbsP◢＇显示直线部分中、边桩高程计算结果≠>Y“HB”=Abs(P-OF)◢不需放样显示高程,否则显示放样水准尺计数⊿Goto 1Lbl 3S=（L-E）（F+G）÷（V-U）-F  ＇前缓曲线超高横坡度计算S≤F=> Goto 5：≠>S≤G=> Goto 6 ＇超高起点(终点)至零界面、零界面至全超高及全超高≠> Goto 7 ⊿⊿ 部分的计算判断条件Lbl 4S=（V-L）（F+G）÷（V-U）-F  ＇后缓曲线超高横坡度计算S≤F=> Goto 5：≠>S≤G=> Goto 6 ＇超高起点(终点)至超高变坡零界面、零界面至全≠> Goto 7 ⊿⊿ 超高、全超高部分计算判断条件Lbl 5O=0=>Z“H1”=AbsP◢＇超高起点(终点)至超高变坡零界面中、边桩高程≠>Y“H-D1”=Abs(P-OF)◢计算结果显示Q“H-G1”=Abs(P+OS)◢⊿Goto 1Lbl 6O=0=>Z“H2”=AbsP◢＇超高变坡零界面至全超高中、边桩高程≠>Y“H-D2”=Abs(P-OS)◢计算结果显示Q“H-G2”=Abs(P+OS)◢⊿Goto 1Lbl 7O=0=>Z“H3”=AbsP◢＇全超高中、边桩高程计算结果显示≠>Y“H-D3”=Abs(P-OG)◢Q“H-G3”=Abs(P+OG)◢⊿Goto 1∙∙注：W“SXG”-视线高N-施工层厚度L"K0+0"-待求点桩号О"B="-中桩至边桩距离（二）子程序：BL≤ZY2=>A=JD1:B=H1:C=I0:D=I1: R=R1:F=I h:E=ZH1:U=YH1:V=HZ1:G=I max1 :≠>L≤ZY3=>A=JD2:B=H2:C=I1:D=I2: R=R2:F=I h:E=ZH2:U=YH2:V=HZ2:G=I max2 :≠>L≤ZY4=>A=JD3:B=H3:C=I2:D=I3: R=R3:F=I h:E=ZH3:U=YH3:V=HZ3:G=I max3 :≠>L≤ZY4=>A=JD3:B=H3:C=I2:D=I3: R=R3:F=I h:E=ZH3:U=YH3:V=HZ3:G=I max3 :≠> ……………………L≤ZY n+1=>A=JD n:B=H n:C=I n-1:D=I n: R=R n:F=I h:E=ZH n:U=YH n:V=HZ n:G=I maxn注:ZY N+1:-下一个竖曲线起点桩号（如：变坡点2的计算范围为前竖曲线终点YZ1至后竖曲线起点ZY3，则在程序中“L≤ZY3”就应输入ZY3的里程桩。

巧用卡西欧FX—4800 P计算器快速精确测放竖曲线

论文申报表巧用卡西欧FX—4800 P计算器快速精确测放竖曲线李肇霞【摘要】巧妙利用CASIO FX—4800P计算器本身的内设函数，修改曲线计算公式和方法，从而求出道路纵断面上任意点的中桩高程，并由此推出横断面上任意点的高程，结合水准仪可精确测定线路内任意点的高程或与设计高程之差。

【关键词】自动；快速；精确；任意测定；1前言近年来在高等级公路铁路建设中对道路的质量要求越来越高，尤其是对线路纵断面的要求。

设计院为了提高道路的行车质量，不断优化线路纵断面，在许多高速公路中竖曲线占路线总长的比例都在80%以上，竖曲线的数量已远远超过平面曲线。

在实际施工中如何快速准确地确定线路的高程是测量放样工作中的一大难题。

竖曲线上各点高程的计算程序设计，要在PC—1500袖珍计算机或计算机上完成程序设计比较容易，但要在FX—4800P计算器上完成线路各点的高程的计算是比较繁琐的。

FX—4800P计算器变量只有26个，各种函数也受到限制。

但是考虑到FX—4800P计算器使用比较普及，价格低，为了充分发挥现有设备的作用，本文介绍了一套简单实用的新程序，供同行们参考使用。

2 fx-4800P的特点从所周知，CASIO fx -4500PA编程函数计算器早在20世纪80年代就已经在我国工程界得到了广泛的应用，fx -4800P是比fx -4500PA更高一个级别的产品，它采用硬翻盖设计，除了具有fx -4500PA的全部功能以外，还主要改进了下列特性：2.1 具有16列×4行显示屏（fx -4500PA只有12列×2行），可以显示更多的信息，这对于编写与修改程序特别方便。

2.2 在程序区域中最多可以存储4500字节的程序与数据（fx -4500PA只能存储1103字节）。

2.3 求方程的数值解。

例如可以求方程2X2+6X-9＝0的数值解。

2.4 复数计算功能，例如可以进行下列复数表达式的运算√(3+4i)=1.27201965+0.786151378i2.5 一、二次微分计算，例如_d_(x3+4x2+x-6)∣x=3 = 52dx_d2_(x3+4x2+x-6)∣x=3 = 26dx22.6 求和计算，例如6∑（K2-3K+5）＝55K＝22.7 fx -4800P是使用15位的尾数进行内部计算，但小数位最多只能显示9位数字。

CASIO fx-4800P程序编制(里程中边桩计算坐标及反算)

一、源程序(SZ-XY-SZ)"1.SZ => XY"："2. XY => SZ"：NU"X0"V"Y0"O"S0"G"F0"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：X"XS" = U+Rec(C,F)+ZcosE◢Y"YS"=V+J+ZsinE◢Goto 1←┙Lbl 2：W=G+90Q：E=U+Rec(R,W)：F=V+J←┙Lbl 3：{XY}：XY：D=Pol(X-E,Y-F):T=J：J<0=>T=J+360△B=T-W-180M：P=Abs B：S"S"=O+ DKcosP+ MπRP÷180◢Z"Z"=DKsinB+MQ(R-D)◢Goto 3←┙二、使用说明1.由里程、中边桩计算坐标启动程序后按以下提示输入数据：N ？1 （N=1，由里程中边桩计算坐标；N≠1，由坐标反算里程及中边桩)X0 ？142.353 (直线或圆曲线起点的X坐标）Y0 ？368.166 (直线或圆曲线起点的Y坐标）S0 ？155.252 (直线或圆曲线起点的里程）F0 ？60°54′03″(直线或圆曲线起点的切线方位角）Q ？-1 (当线元为直线时Q=0；当线元为圆曲线且左偏时Q=-1；当线元为圆曲线且右偏时Q=1)R ？100 (圆曲线的半径)S ？194.016 (中线里程)Z ？-5.251 (中线左侧5.251m的边桩；在中线上Z=0，在左侧<0，在右侧Z>0)显示计算结果XS=172.7467767YS=395.23168072.由坐标反算里程、中边桩启动程序后按以下提示输入数据：N ？2 （N=1，由里程中边桩计算坐标；N≠1，由坐标反算里程及中边桩)X0 ？142.353 (直线或圆曲线起点的X坐标）Y0 ？368.166 (直线或圆曲线起点的Y坐标）S0 ？155.252 (直线或圆曲线起点的里程）F0 ？60°54′03″(直线或圆曲线起点的切线方位角）Q ？-1 (当线元为直线时Q=0；当线元为圆曲线且左偏时Q=-1；当线元为圆曲线且右偏时Q=1)R ？100 (圆曲线的半径)X ？172.7467767 (待求里程及中边桩的点的X坐标)Y ？395.2316807 (待求里程及中边桩的点的Y坐标)显示计算结果S=194.016Z=-5.251三、注意事项1.本程序主要是针对市政道路工程中的直线和圆曲线而编制的，未考虑缓和曲线等曲线的计算；2.程序中“R=1E-45”的意思为10的-45次方，输入程序时切莫输错。

CASIO4800综合曲线坐标计算程序

[正确] CASIO fx-4800P综合曲线坐标计算程序说明：该程序适用于计算器 CASIO fx-4800P，可计算线路中心的缓和曲线、圆曲线、直线段，中、边桩坐标及切线方位角。

1、A? 输入转角:左转为负，右转为正2、R? 输入圆曲线半径3、LS? 输入缓和曲线长度4、JD(DK)? 输入交点里程桩号5、X(JD)? 输入本交点X坐标6、Y(JD)? 输入本交点Y坐标7、FWJ? 输入待求点切线方位角9、J? 输入0程序计算中桩，输入1程序计算边桩10、Z? 输入里程桩号1 A：R：C“LS”：D“JD(DK)”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◢6 W=(R+P)/cos(A/2)-R◢7 L=((AbsA)-2B)兀R/180+2C◢8 G“ZH”=D-T◢9 H“HY”=G+C◢10 I“QZ”=G+L/2◢11 K“YH”=G+L-C◢12 M“HZ”=G+L◢13 N＂X(JD)＂：E＂Y(JD)＂：F＂FWJ＂: J14 A＜0=>S=-1:≠=>S=1⊿ (提示：0为数字“0”)15 U=F+A/2+90S16 V=W+R17 B=N+VcosU18 O=E+VsinU (提示：O为字母“O”)19 Lbl 120 {Z}21 Z≤G=>L=T+G-Z22 V=F+18023 U=F25 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“FWJ”=F+3L◢31 Goto 4⊿32 Z≤K=>L=F+A/2+90S+180+180(Z-I)S/R/兀33 U“FWJ”=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+A◢41 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}52 P“XL”=X+Wcos(U-90) ◢53 Q“YL”=Y+Wsin(U-90) ◢54 P“XR”=X+Wcos(U+90) ◢55 Q“YR”=Y+Wsin (U+90) ◢56 Goto 157 Lbl 458 X=N+TcosP+VcosQ◢59 Y=E+TsinP+VsinQ◢60 Goto 661 Lbl 562 X=B+R cosL◢63 Y=O+RsinL◢ (提示：O为字母“O”)65 Lbl 666 J=1=> Goto 3⊿67 Goto 1注:1、◢为输出指令，若在后面加上，即可显示前面的计算结果输出在屏幕上。