fx4800P曲线任意里程中边桩坐标正反算程序
卡西欧计算器坐标计算程序4800
CASIO fx-4800P计算器程序DA DI——N E (公路弯道坐标计算)K“JD”：R：S“LS”：F“PJ——L－R＋”：A“ZH——N”：B“ZH——E”：C“JD——N”：D“JD——E”M“ZH”=K－（R＋S2÷（24R））tan（Abs F÷2）－S÷2＋S^3÷（2 40R2）▲N“HY”=M＋S▲O“QZ”=M＋（πRAbs F÷180＋S）÷2▲Z[1]“YH”=M＋πRAbs F÷180▲Q“HZ”= Z[1]＋S▲T=K－M▲L=Q－M▲E=（R＋S2÷（24R））÷cos（F÷2）－R▲Lb1 0：{G}：G“ZHUANG HAO”：G＜O＝=＞Z=G－M≠=＞Z=Q－G：△{U}：U“I——B OUT＋IN－”Z＜S＋0.005＝=＞H=tan-1（24Z2S2R2－Z^6）÷（48Z^3R^3－6Z^4 SR））：V=Z－Z^5÷（40S2R2）＋Usin H：W=Z^3÷（6SR）－Z ^ 7÷（336Z^3R^3）－Ucos H：≠=＞H=180（Z－S÷2）÷（πR）：V=Rsin H＋S÷2－S^3÷（240R2）＋Usin H：W=R－Rcos H＋S2÷（24R）－Ucos H：△F＜0＝=＞G＜O＝=＞X=V：Y=W：≠=＞X=T＋Tcos F－Vcos F－Ws in Abs F：Y=Tsin Abs F－Vsin Abs F＋Wcos F：△≠=＞G＜O＝=＞X=V：Y=－W：≠=＞X=T＋Tcos F－Vcos F－Wsin Abs F：Y=Vsi n Abs F－Tsin Abs F－Wcos F：△△C=A＝=＞D≥B＝=＞P=90：≠=＞P=270：△≠=＞C＜A＝=＞P= tan-1（（D－B）÷（C－A））＋180：≠=＞P= tan-1（（D－B）÷（C－A））：△△I“COM——N”=A＋√（X2＋Y2）cos（P－tan-1（Y÷X））▲J“COM——E”=B＋√（X2＋Y2）sin（P－tan-1（Y÷X））▲G=G＋20Goto 0程序运行输入：JD？——交点桩号R？——圆曲线半径LS？——缓和曲线长度PJ——L－R＋？——偏角右偏为正ZH——N？——ZH点N坐标ZH——E？——ZH点E坐标JD——N？——JD点N坐标JD——E？——JD点E坐标输出：ZH= ——ZH点桩号HY= ——HY点桩号QZ= ——QZ点桩号YH= ——YH点桩号HZ= ——HZ点桩号T= ——切线长L= ——曲线长E= ——外矢距再输入：ZHUANG HAO？——输入任一点桩号I——B OUT＋IN－？——距中线距离（外侧为正）再输出：COM——N= ——计算出N坐标COM——E= ——计算出E坐标以此循环——再输入——再输出程序中各字母所代表的含义：A——ZH点N坐标B——ZH点E坐标C——JD点N坐标D——JD点E坐标E——外矢距F——偏角（右偏为正）G——弯道内任一点桩号H——弯道内任一点切线角I——计算出的N坐标J——计算出的E坐标K——交点桩号L——曲线长M——ZH点桩号N——HY点桩号O——QZ点桩号P——ZH--JD方位角Q——HZ点桩号R——圆曲线半径S——缓和曲线长度T——切线长U——距中线距离（外侧为正）V——支距W——支距X——弯道内坐标Y——弯道内坐标Z——任一点曲线长度Z[1]——YH点桩号。
4800计算器程序
CASIO4800坐标计算程序CASIO4800缓和曲线、圆曲线和直线段上的任意中、边桩坐标计算程序1、该程序适用于计算器 CASVO fx-4800P，可计算与线路中心成任意夹角的缓和曲线、圆曲线、直线段中、边桩坐标及切线方位角。
2、A？输入转角：左转为负，右转为正3、R？输入圆曲线半径4、LS？输入缓和曲线长度5、JD？输入交点桩号6、J？输入0程序计算中桩，输入1程序计算边桩7、N？ E？输入交点X、Y坐标8、FW 待求点切线方位角9、当计算交角为60°时，则公式改为（U+60），（U-120）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◢6 W=(R+P)/cos(A/2)-R◢7 L=((AbsA)-2B)兀R/180+2C◢8 G=D-T◢9 H=G+C◢10 I=G+L/2◢11 K=G+L-C◢12 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◢31 Goto 4⊿32 Z≤K=>L=F+A/2+90S+180+180(Z-I)S/R/兀33 U“FW”=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◢64 Goto 665 Lbl 666 J=1=> Goto 3⊿67 Goto 1。
曲线(含直线)任意里程中边桩坐标正反算
曲线（含直线）任意里程中边桩坐标正反算（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计算器曲线计算程序
线路测设坐标计算程序（Casio fx-4800p）在铁路、公路工程勘测和施工测量中，相对复杂的线路测设工作占了很大比重，且传统线路测设公式已不能满足现代化曲线测设的需要。
为此，笔者根据工程测量的特点，严格定义了坐标系，推导了线路测设的通用计算公式。
该公式不但可以计算单交点平曲线坐标，同时还可以计算复曲线坐标，如不完整缓和曲线等。
根据该公式，笔者编写了Casio fx-4800p 计算器上的单交点平曲线中桩、边桩坐标计算程序。
该程序可以判断已知数据输入的正确性，并可限制输入的测设点里程。
在程序编制过程中，已考虑了零溢出问题。
一．使用说明1.起动名为“ZXZB”的计算程序。
2.输入已知数据：X(QD)—起点纵坐标(EXE)Y(QD)—起点横坐标(EXE)X(JD)—交点纵坐标(EXE)Y(JD)—交点横坐标(EXE)S(JD)—交点里程(EXE)X(ZD)—终点纵坐标(EXE)Y(ZD)—终点横坐标(EXE)RY —圆曲线半径(EXE)L1 —前缓和曲线长度(EXE)RS —前缓和曲线起始半径（0表示无穷大）(EXE)L2 —后缓和曲线长度(EXE)RE —后缓和曲线结束半径（0表示无穷大）(EXE)S(HZ1) —起点缓直点里程(EXE)S(ZH1) —交点直缓点里程(EXE)S(HZ2) —交点缓直点里程(EXE)S(ZH2) —终点直缓点里程(EXE)S? —测设点里程(EXE)X=*** —显示测设点纵坐标(EXE)Y=*** —显示测设点横坐标(EXE)I? —边桩距（左负右正）(EXE)XB=*** —显示边桩纵坐标(EXE)YB=*** —显示边桩横坐标(EXE)I? —输入0结束边桩计算，继续下一测设点计算。
二．计算举例1．已知数据：X(QD)=2175.199 Y(QD)=1015.467X(JD)=2000 Y(JD)=2000S(JD)=6000X(ZD)=2621.088 Y(ZD)=2783.741RY=600L1=220 RS=5800L2=180 RE=3300S(HZ1)=5200 S(ZH1)=5631.047S(HZ2)=6311.047 S(ZH2)=67002.求：各测设点中桩坐标及左边桩10M、右边桩20M坐标3.计算结果三．源程序清单ZXZB1. Defm 15：Rad：Norm↙2. A＂X(QD)＂:B＂Y(QD)＂:C＂X(JD)＂:D＂Y(JD)＂:Z＂S(JD)＂:E＂X(ZD)＂:F＂Y(ZD)＂:N＂RY＂↙3. Q＂L1＂:R＂RS＂:V＂L2＂:W＂RE＂:G＂S(HZ1)＂:H＂S(ZH1)＂:K＂S(HZ2)＂:L＂S(ZH2)＂:Fixm↙4. O＝H＋Q:P＝K－V:M＝Pol (C－A,D－B):T＝J:Z[3]＝J: Z[1]＝C＋Rec (Z－H,J＋π):Z[2]＝D＋J↙5. M＝Pol(E－C,F－D):Z [12]＝J:T＝J－T:Abs T＞π＝＞T＞π＝＞T＝T－2π: ≠＞T＝T＋2π△△↙6. T＞0＝＞Z[13]＝1:≠＞Z[13]＝－1△↙7. R＝0＝＞Z[14]＝0:≠＞Z[14]＝QN÷(R－N)△↙8. W＝0＝＞Z[15]＝0:≠＞Z[15]＝VN÷(W－N)△↙9.Q＝0＝＞I＝0:≠＞I＝(2QZ[14]＋Q²) ÷2÷N÷( Z[14]＋Q)△↙10.V＝0＝＞J＝0:≠＞J＝(2V(Z[15]＋V)－V²)÷2÷N÷( Z[15]＋V)△↙11.Z[6]＝Z[3]＋I Z[13]: Z[9]＝Z[12]－JZ[13] ↙12.T＝Abs T－I－J:J＝TN:Abs(P－O－J)＞0.05＝＞＂ERROR＂▲Goto 1△↙13.M＝Q:Prog ＂ZHD＂:Z[4]＝X:Z[5]＝Y:M＝J:Porg ＂HYD＂:Z[7]＝X:Z[8]＝Y:M＝V:Prog ＂YHD＂: Z[10]=X:Z[11]＝Y↙14.Prog ＂ZXZS＂↙15.Lbl 1:0:＂＂:Pause 0ZXZS1.Rad:Norm:Fixm2.Lbl 2:{S}:S＜G＝＞＂RANGE OUT＂▲Goto 5△↙3.S＞L＝＞＂RANGE OUT＂▲Goto 5△↙4.S＜H＝＞M＝Z[3]:X＝Z[1]＋Rec(S－H,M):Y＝Z[2]＋J:Goto 3△↙5.S＜O＝＞M＝S－H:Prog ＂ZHD＂:Goto 3△↙6.S＜P＝＞M＝S－O:Prog ＂HYD＂:Goto 3△↙7.S＜K＝＞M＝S－P:Prog ＂YHD＂:Goto 3△↙8.M＝Z[12]:X＝Z[10]＋Rec(S－K,M):Y＝Z[11]+J↙9.Lbl 3:I＂X＂＝Int (1000X＋0.5)÷1000▲J＂Y＂＝Int(1000Y＋0.5)÷1000▲10.Lb1 4:I＝0:{I}:I＞0＝＞Goto 2△↙11.T＝X＋Rec(I,M+π÷2):U＝Y＋J:T＂XB＂＝Int (1000T+0.5)÷1000▲U＂YB＂＝Int(1000U＋0.5)÷1000▲Goto 4↙12.Lbl 5:0:＂＂:Pause 0ZHD1.Fixm:M＜0.001＝＞T＝0:U＝0:M＝Z[3]:Goto 6△↙2.I＝Z[14]:X＝N(I＋Q):Y＝Z[13]:T＝M－(20I²＋15IM＋3M²)M^3÷120÷X²:T＝T＋((1008I+1680M)I^3＋1080I²M²＋(315I＋35M)M^3)M^5÷120960÷X^4↙3.U＝((70I+84M)I²＋(35I+5M)M²)M^4÷1680÷X^3:U＝(3I＋M)M²÷6÷X－U:U＝U＋((3696I＋7920M)I^4＋(6930I＋3080M)I²M²＋(693I＋63M)M^4)M^6÷2661120÷X^5:U＝UY↙4.M＝(2IM+M²)÷2÷X:I＝Z[3]:M＝I＋MY↙5.Lbl 6:X＝Z[1]＋Tcos I－Usin I:Y＝Z[2]＋Tsin I＋Ucos IHYD1.Fixm:M＝M÷N:I＝Z[6]:Y＝Z[13]:T＝Nsin M:U＝(N－Ncos M)Y:M＝I＋MY2.X＝Z[4]＋Tcos I－Usin I:Y＝Z[5]＋Tsin I＋Ucos IYHD1.Fixm:M＜0.001＝＞T＝0:U＝0:M＝Z[9]:Goto 7△↙2.I＝Z[15] ＋V:X＝NI:Y＝Z[13]:T＝M－(20I²－15IM＋3M²)M^3÷120÷X²:T＝T＋((1008I－1680M)I^3＋1080I²M²－(315I－35M)M^3)M^5÷120960÷X^4↙3.U＝((70I－84M)I²＋(35I－5M)M²)M^4÷1680÷X^3:U＝(3I－M)M²÷6÷X－U:U＝U＋((3696I－7920M)I^4＋(6930I－3080M)I²M²＋(693I－63M)M^4)M^6÷2661120÷X^5:U＝UY↙4.M＝(2IM－M²)÷2÷X:I＝Z[9]:M＝I＋MY↙5.Lbl 7:X＝Z[7]＋Tcos I－Usin I:Y＝Z[8]＋Tsin I＋UCos I注: 1.由于本计算器程序容量为4500步，本程序只能计算单交点平曲线；2.当已知数据输入存在错误时，程序能够自动判断错误后退出；3.当输入的测设点里程超出范围时，计算自动终止；4.编制程序过程中已经考虑了零溢出问题。
fx-4800计算程序
标准测量曲线坐标计算程序1.文件名MAILLbi 0{U,D}（U为桩号，D为所求点到中点的距离左为负，右为正）U<12111.892(ZH点)=>PROG“JD1”：≠>U《12000=>“JD2”⊿⊿（有多少语句就有多少⊿）X◢Y◢GOTO 02.数据记录文件“JD1”Z=12000（直缓点桩号）H=13000（缓直点桩号）L=100（缓和曲线长度）T=100（切线长）R=2000（圆曲线半径）A=175.2315（第一方位角）B=150.2000（第二方位角）I=+ -1（左为负，右为正）P=（交点X坐标）Q=（交点Y坐标）V=（缓圆点X坐标）W=（缓圆点Y坐标）PROG“SE”3.“SE”判断J=Z+LG=H-LU〈Z=> PROG“LINE”：≠>U〈J=> PROG“Q1”：≠>U〈G=> PROG“CIR”：≠>U〈H=> PROG “Q2”⊿⊿⊿⊿（有多少语句就有多少⊿）4. “LINE”直线S=U-ZC=AX=P+（S+T）*COS（A+180）+D*COS（A+90）Y=Q+（S+T）* SIN（A+180）+D*SIN（A+90）5. 第一缓和曲线“Q1”S=U-ZC=A+I*S2/2/L/R*180/PIM=S-S5/40/R2/L2+S9/3456/（RL）4N=S3/6/L/R-S7/336/L3/R3+S11/42240/（RL）5X=P+T*COS（A+180）+M*COSA+N*COS（A+I*90）+D*COS（C+90）Y=Q+T*SIN （A+180）+M*SIN A+N* SIN（A+I*90）+D* SIN（C+90）6.圆曲线“CIR”S=U-JC=A+I*（L/2/R+S/R）*180/PIX=V+2*R*SIN（S/2/R*180/PI）*COS（A+I*（（L+S）/2/R）*180/PI）+D*COS （C+90）Y=W+2*R*SIN（S/2/R*180/PI）*SIN（A+I*（（L+S）/2/R）*180/PI）+D*SIN （C+90）7.第二缓和曲线“Q2”S=H-UC=B- I*S2/2/L/R*180/PIM=S-S5/40/R2/L2+S9/3456/（RL）4N=S3/6/L/R-S7/336/L3/R3+S11/42240/（RL）5X=P+M*COS（B+180）+T*COSB+N*COS（B+I*90）+D*COS（C+90）Y=Q+M*SIN （B+180）+T*SIN B+N* SIN（B+I*90）+D* SIN（C+90）匝道程序1.主程序（TYQXJS）“1.SZ=XY”：“2. XY =SZ”：N：U“X0”：V“Y0”：O“S0”：G“F0”：H“LS”：P“RO”：R“RN”：Q：C=1/P：D=（P-R）/（2HPR）：E=180/PI：N=1=>GOTO 1：≠>GOTO 2⊿Lbi 1：{SZ}：SZ：W=AbS（S-O）：PROG“SUB1”：X“XS”=X⊿Y“YS”=Y⊿GOTO 1Lbi 2：{XY}。
CASIO fx-4800P程序
1、方位角计算程序:【FWJ】 (2)2、转换坐标与高程程序:【0】 (2)3、坐标计算主程序:【XY】 (2)4、高程计算主程序:【H】 (5)5、坐标计算子程序:【XY-1】 (5)6、高程计算子程序:【H-1】 (6)7、水平距离归算到高斯平面计算程序:【JLGF】 (8)8、任意多边形面积计算程序:【MJ】 (8)9、四边形平差计算程序:【DDSBX】 (8)10、坐标换带计算程序:【ZBFD】 (11)11、两直线交点计算程序:【L-XY-L】 (12)12、高斯投影坐标计算程序:【GSBL-XY】 (12)13、竖曲线计算程序:【SQX】 (13)14、复曲线计算程序:【FQX】 (13)15、匝道计算子程序:【FQX】 (14)1、方位角计算程序:【FWJ】置镜点A＂X【Z】＂:B＂Y【Z】＂:后视点C＂X【H】＂:D＂Y【H】＂:前视点E＂X【Q】＂:F＂Y【Q】＂↵M=C-A:N=D-B:V=E-A:W=F-B:L=√（M2+N2）:S=√（V2+W2）:I=AbsTan-1（N÷M）:Q= AbsTan-1（W÷V）↵M＞0＝＞Goto 1:≠＝＞Goto 2↵Lbl 1↵N＞0＝＞O=I:Goto 3:≠＝＞O=360-I:Goto 3↵Lbl 2↵N＞0＝＞O=180-I:Goto 3:≠＝＞O=180+I:Goto 3↵Lbl 3↵V＞0＝＞Goto 4:≠＝＞Goto 5↵Lbl 4↵W＞0＝＞H=Q:Goto 6:≠＝＞H=360-Q:Goto 6↵Lbl 5↵W＞0＝＞H=180-Q:Goto 6:≠＝＞H=180+Q:Goto 6↵Lbl 6↵O＞H＝＞X=O-H:Goto 9:≠＝＞X=O-H+360:Goto 9↵Lbl 9↵＂F【Z-H】=＂:O ◢＂L【Z-H】=＂:L ◢＂F【Z-Q】=＂:H ◢＂L【Z-Q】=＂:S ◢＂f【H-Z-Q】=＂:360-X（左角）◢＂f【Q-Z-H】=＂: X（右角）计算示例:2、转换坐标与高程程序:【0】计算选择:1.坐标计算;2.高程计算。
坐标正算与反算、施工与大地坐标换算程序(4800P)
Goto 1
算例：1点，X=10.000，Y=20.000
2点，X=500.000，Y=600.000
3点，X=700.000，Y=400.000
专业的3S站
DB子程序：V=V+1:B=Z[V]+T:Prog“A”:D=Z[V+N]:Prog“XY”
DEG子程序：P=Frac B*100:B=INTB+INTP/60+Frac P/36
三、 支导线的计算
支导线就是从一个已知点出发，既不符合到另一个已知点也不回到原来的点上。支导线没有检核条件，不易发现错误，过不宜采用。
X？ 输入5点X坐标
Y？ 输入5点Y坐标
S=218500 （1.2.3.4.5点围成的面积）
X？ 输入6点X坐标
Y？ 输入6点Y坐标
S=230500 （1.2.3.4.5.6点围成的面积）
.........
.........
Lbl 3：｛XY｝：X“X=”：Y“Y=”：A=（X－C）cos R+(Y－D) sin R:“A=” ：Pause 0:B=（Y－D）cos R－(X－C) sin R:B“B=” ▲ Goto 0 △
Lbl 4：｛AB｝：A“A=”：B“B=”：X=C+A cos R－B sin R：“X=”：Pause 0:Y=D+A sin R+B cos R：Y“Y=” ▲ Goto 0 △
坐标正算与反算、施工与大地坐标换算程序（4800P）
Deg：Fix4：J=0：C“X=”：D“Y=”：R“R=“：
单条曲线坐标计算程序(新) 卡西欧4800计算器
单条曲线坐标计算程序(casio fx-4800P)(只限于对称型综合曲线使用)包含1个主程序和2个子程序1、主程序: “QU-XIAN”O“JDLC”：A“JDX”：B“JDY”：R：G“L0”：Z“ZUAN-JIAO”：H“FANG-XIANG”：V“Z-1，Y+1”LbL 1{S,M}G>0=>Prog“ZH-Y-HZ”:≠=>Prog“YY”⊿U<360=>U:≠=>U=U-360⊿E=X+cos(U-90)×MF=Y+sin(U-90)×MK=X+cos(U+90)×MP=Y+sin(U+90)×MX◢Y◢E“ZX”◢F“ZY”◢K“YX”◢P“YY”◢U“HUANG-XIANG”◢Goto 12、子程序：“ZH-Y-HZ”T=(R+G2÷24R)tg(Z÷2)+G÷2－G3÷240R2L=R(Z－G÷R×180÷π)×π÷180+2GC=O－TD=C+GI=C+L－GN=C+LE=A+Tcos(H+180)F=B+Tsin(H+180)K=E+cos(H+VG2÷6÷R÷G×180÷π)×(G－G5÷90÷(RG)2)P=F+sin(H+VG2÷6÷R÷G×180÷π)×(G－G5÷90÷(RG)2)Q=A+cos(H+VZ)×TW=B+sin(H+VZ)×TS<c=>L=C－S:X=A+cos(H+180)×(T+L):Y=B+sin(H+180)×(T+L):U=H⊿S>C=>L=S－C:J=H+VL2÷6÷R÷G×180÷π:L=L－L^5÷90÷(RG)2:X=E+cosJ ×L: Y=F+sinJ×L:U=H+3V(S－C)2÷6÷R÷G×180÷π⊿S>D=>L=S－D:J=H+VG÷R×90÷π+90VL÷π÷R:L=2R×sin(90L÷R÷π):X=K+cosJ×L:Y=P+sinJ×L:U=H+VG÷R×90÷π+V×180(S－D)÷π÷R⊿S>I=>L=N－S:J=H+VZ－VL2÷6÷R÷G×180÷π:L=L－L^5÷90÷(RG)2:X=Q+cosJ×(－L): Y=W+sinJ×(－L):U=H+VZ－3V(N－S)2÷6÷R÷G×180÷π⊿S>N=>L=S－N:J=H+VZ:X=A+cosJ×(T+L):Y=B+sinJ×(T+L):U=J⊿3、子程序：“YY”T=Rtg(Z÷2)L=RZ×π÷180C= O－TN=C+LE=A+Tcos(H+180)F=B+Tsin(H+180)S<C=>L=C－S:X=A+cos(H+180)×(T+L):Y=B+sin(H+180)×(T+L):U=H⊿S>C=>L=S-C:J=H+90VL÷πR:X=E+cosJ×(2Rsin(90L÷πR)):Y=F+sinJ×(2Rsin(90L÷πR)):U=H+180VL÷πR⊿S>N=>L=S-N:J=H+VZ:X=A+cosJ×(T+L):Y=B+sinJ×(T+L):U=J⊿使用说明：JDLC--------------交点里程JDX---------------交点X坐标JDY---------------交点Y坐标R-----------------圆曲线半径L0----------------缓和曲线长ZUAN-JIAO---------线路转角FANG-XIANG--------ZH（ZY）方向角V“Z-1，Y+1”-----系数，左转输“-1”，右转输“+1”S-----------------待求点桩号M-----------------边桩距离显示计算结果：X、Y--------------- 待求点中桩坐标ZX、ZY--------------待求点左边桩坐标YX、YY--------------待求点右边桩坐标HUANG-XIANG---------待求点方向角匝道上回旋线元任意点坐标计算程序(casio fx-4800P)文件名：ZA DAO HUANC ：O ：LbL 1{S,J,W}H=1=>M=0△R=0=>G=0:≠=>G=C ÷R △ I=G+S-K H=2=>M=C ÷RD:G=C ÷R:I=G-S+K △ V=1=>L=M+S-K △ V=-1=>L=M-S+K △ T=O-90VM 2÷C ÷π：A=L-M-(L^5-M^5)÷40C 2+(L^9-M^9)÷3456C^4-(L^13-M^13)÷599040C^6：B=(L^3-M^3)÷6C-(L^7-M^7)÷336C^3+(L^11-M^11)÷42240C^5-(L^15-M^15) ÷9676800C^7： U=X+V ACosT-BSinT ◢Z=Y+V ASinT+BCosT ◢P=O+Abs(G 2-I 2) ÷2C×180÷π×DP<360=>P ◢≠=>P=P-360◢△E=U+JCos(P-90) ◢F=Z+JSin(P-90) ◢N=U+WCos(P+90) ◢Q=Z+WSin(P+90) ◢Goto 1一、 使用说明C ——为回旋线曲率半径变化率，C=A 2=RL 0O ——为（ZH 、YH ）方向角S ——回旋上待求点里程K ——回旋线起算点（ZH 、YH ）里程X 、Y ——起算点（ZH 、YH ）的横、纵坐标U 、Z ——待求点横、纵坐标P ——待求点方向角J 、W ——待求点左、右边桩距离E 、F ——待算点左边桩横、纵坐标N 、Q ——待算点右边桩横、纵坐标D ——线路转向判断值，左转为“-1”，右转为“+1”二、 注意事项H ——缓和曲线类型判断值，第一缓和曲线为“1”，第二缓和曲线为“2”（起算点半径为无穷大为第一缓和曲线，起算点为YH 点时为第二缓和曲线）V ——曲率半径大小比较值，曲率半径ρ由大变小为“+1”，反之为“-1”，线路左转时，曲率半径ρ要带负号进行比较R ——起算点（ZH 、YH ）的半径值，不需输入正负号，在ZH 时，半径为∞（无穷大），此时半径输入0 。
卡西欧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Δ←┘(第三线元数据要素)。
卡西欧4800坐标正反算通用程序(终极篇)
卡西欧4800坐标正反算通用程序(终极篇)1. 坐标正算主程序(命名为ZBZS)第1行：Lbl 0:{W,T,B}:W”K=”:T”LEN=”:B第2行：Prog “A”第3行：X=D+Tcos(F+B)◢第4行：Y=E+Tsin(F+B)◢第5行：F=F◢第6行：Goto 0K——计算点的里程LEN——计算点到中桩的距离（左负右正）B——取前右夹角为正2. 坐标反算桩号和偏距主程序(命名为ZBFS)第1行：{U,V,K}:U”X1”:V”Y1”: W”K1”第2行：I=0:J=0第3行：Lbl 0:Prog “A”第4行：Pol(U-D,V-E):S=Icos(F-J):W=W+S第5行：Abs(S)>0.0001=>Goto 0△第6行：T=Isin(J-F)第7行：K=W◢第8行：T”LEN”=T◢X1——取样点的X坐标Y1——取样点的Y坐标K1——输入时为计算起始点(在线路内即可)，输出时为反算点的桩号T——偏距（左负右正）3. 计算坐标子程序(命名为XYF)为了简洁，本程序由数据库直接调用，上述中的正反算主程序不直接调用此程序第1行：S=W-A：I=(Q-P)÷L第2行：D=D+∫(cos(F+X(2P+XI)×90÷π),0,S,4)第3行：E=E+∫(sin(F+X(2P+XI)×90÷π),0,S,4)第4行：F=F+S(2P+SI)×90÷π4. 数据库（命名为A）第1行：W≤175.191=> GoTo 2△第2行：A=175.191:D=428513.730:E=557954.037:F=92°26′40″:P=0:Q=1/240:L=70.417: W≤A+L =>GoTo 1△第3行：A=245.607: D=428507.298:E=558024.092:F=100°50′59.4″:P=1/240:Q=1/240: L=72.915: W≤A+L =>Goto 1△第4行：A=318.522: D=428482.988:E=558092.538:F=118°15′25.2″:P=1/240:Q=0:L=55.104: W≤A+L =>Goto 1△第5行：A=373.627:D=428453.283:E=558138.912:F=124°50′4.5″0:P=0:Q=-1/180:L=67.222:W≤A+L=>Goto 1△Goto 2第6行：Lbl 1:Prog “XYF”: Goto 3第7行：Lbl 2:D=0:E=0:F=0第8行：Lbl 3A——曲线段起点的里程D——曲线段起点的x坐标E——曲线段起点的y坐标F——曲线段起点的坐标方位角P——曲线段起点的曲率(左负右正)Q——曲线段终点的曲率(左负右正)L——曲线段长度（尽量使用长度，为计算断链方便）说明：（1）正算主程序可以计算一般边桩的坐标，如要计算类似涵洞端墙的坐标需增加两个变量，具体方法参考本程序集中的第1篇辛普生公式的坐标计算通用程序（2）程序规定，左偏曲线曲率（半径倒数）输入负值，右偏曲线曲率输入正值，直线上点曲率输入0。
任意曲线坐标放样正反算程序(4800,4850,5800均可以用)
任意曲线坐标放样正反算程序（fx－4800p）正算ZSLbl 0↙｛K,L,Q,S｝:K:L:Q:Prog“FXL”:Prog“SUB”↙“X=”:X=X:Pause0:Y=Y◢K=S+K↙Goto0↙反算FSLbl 0:｛PT｝:P“X0”:T“Y0”↙Lbl 1:｛K｝:L=0:Q=0:Prog“FXL”: Prog“SUB”↙U=I-90↙V=（T-Y）CosU-（P-X）SinU↙AbsV≥0.001=>K=K+V:Goto1:≠>“K0=”:K=K:Pause0:W“L0”=（T-Y）CosI-（P-X）SinI◢Goto0↙子程序SUBNorm↙Deg↙U=（E-D）÷Abs（G-F）:V=Abs（K-F）:W=UV↙J=C+90V（W+2D）÷π:J<0=>J=J+360⊿J≥360=>J=J-360⊿↙M=C+45V（W÷4+2D）÷（2π）:N=C+135V（3W÷4+2D）÷（2π）:O=C+45V（W÷2+2D）÷π↙I=J+Q↙Fix3↙X=A+V÷12×（CosC+4（CosM+CosN）+2CosO+CosJ）+L CosI↙Y=B+V÷12×（SinC+4（SinM+SinN）+2SinO+SinJ）+L SinI↙数据库FXLK≥0=>K≤1211.753=>A=9764.547:B=3344.148:C=142，6，26，:D=0:E=0:F=0:G=1211.753⊿⊿↙K≥1211.753=>K≤1406.455=>A=8808.279:B=4088.39:C=142，6，26，:D=1÷700:E=1÷700:F=1211.753:G=1406.455⊿⊿↙使用说明：K－里程桩号; A,B－每段线元的起点X,Y坐标; C－每段线元的起点的切线方位角; F－每段线元的起点桩号;G－每段线元的终点桩号;D－每段线元的起点半径曲率;E－每段线元的终点半径曲率; （左转取-1÷半径曲率,右转取1÷半径曲率,直线为下一个桩号（例: 赋值20确认,即上个桩号加20）D,E取值说明：当线元为直线段,不论其起,止与何线元相接, D,E都取0.当线元为圆曲线,不论其起,止与何线元相接, D,E都取(右偏)1/(左偏)-1÷该圆曲线的半径.当线元为缓和段, 起,止与直线相接, D,E都取0. 与圆曲线相接, D,E都取(右偏)1/(左偏)-1÷该圆曲线的半径. 与缓和段相接, D,E都取0.任意点的切线方位角C说明：编数据库时,若两曲线相接,例：先编ZH至HY段,运行算出HY点的切线方位角,(即求HY点的坐标,然后终止程序调出C值),然后再以计算出的C编HY至YH段.依此类推.反算桩号时需多次确认所输入的X,Y坐标,程序运算至桩号误差小于允许值时,会同时显示结果:K0-桩号,L0-距离.坐标反算程序:ZBFSC“X1”：D“Y1”：E“X2”：F“Y2”：Fixm：Pol(E－C,F－D：I“S=”◢J≤0=>J=J+360⊿J“J=”↙坐标正算程序ZBZSC“X1”：D“Y1”：L：O“J”：Fixm：X=C+Rec(L,O◢Y=D+J↙竖曲线SQXJ“BPD”：B“BPH”：C“I1”：D“I2”：RW=D－C◢A=W÷AbsW：T=Abs（0.5RW）◢L=2T↙E=T2÷2R◢X“QD”=J－T◢Y“ZD”=J+T◢Lbl1：P“K”：H=B+C（P－J）+A（P－J+T）2÷2R◢W“CONT”：W=1=>Goto1：≠>Goto2：⊿↙Lbl2↙说明:BPD－变坡点桩号BPH－变坡点高程I1－前一个坡度（上“+”,下“-”）I2－后一个坡度（上“+”,下“-”）R－半径E－外距QD－该竖曲线起点桩号ZD－该竖曲线终点桩号K－要求的桩号H－要求桩号的高程任意曲线坐标放样正反算程序（fx－4850p）正算ZSLbl 0↙｛K,L,Q,S｝:K:L:Q:Prog“FXL”:Prog“SUB”↙“X=”:X=X:Pause0:“Y=”:Y=Y◢K=S+K↙Goto0↙反算FSLbl 0:｛PT｝:P“X0”:T“Y0”↙Lbl 1:｛K｝:L=0:Q=0:Prog“FXL”: Prog“SUB”↙U=I-90↙V=（T-Y）CosU-（P-X）SinU↙AbsV≥0.001=>“K0=”:K=K+V:Goto1:≠>“K0=”:K=K:Pause0:W=（T-Y）CosI-（P-X）SinI ↙“L0=”:W=W◢Goto0↙子程序SUBNorm↙Deg↙U=（E-D）÷Abs（G-F）:V=Abs（K-F）:W=UV↙J=C+90V（W+2D）÷π:J<0=>J=J+360⊿J≥360=>J=J-360⊿↙M=C+45V（W÷4+2D）÷（2π）:N=C+135V（3W÷4+2D）÷（2π）:O=C+45V（W÷2+2D）÷π↙I=J+Q↙Fix3↙X=A+V÷12×（CosC+4（CosM+CosN）+2CosO+CosJ）+L CosI↙Y=B+V÷12×（SinC+4（SinM+SinN）+2SinO+SinJ）+L SinI↙数据库FXLK≥0=>K≤1211.753=>A=9764.547:B=3344.148:C=142，6，26，:D=0:E=0:F=0:G=1211.753⊿⊿↙K≥1211.753=>K≤1406.455=>A=8808.279:B=4088.39:C=142，6，26，:D=1÷700:E=1÷700:F=1211.753:G=1406.455⊿⊿↙使用说明：K－里程桩号; A,B－每段线元的起点X,Y坐标; C－每段线元的起点的切线方位角;D－每段线元的起点半径曲率;E－每段线元的终点半径曲率; （左转取-1÷半径曲率,右转取1÷半径曲率,直线为0. ）;L－距中桩距离（左负,右正,中桩为0）; Q－与线路前进方向的右交角（正交为90）; S－赋值加上一个桩号放下一个桩号（例: 赋值20确认,即上个桩号加20）D,E取值说明：当线元为直线段,不论其起,止与何线元相接, D,E都取0.当线元为圆曲线,不论其起,止与何线元相接, D,E都取(右偏)1/(左偏)-1÷该圆曲线的半径.当线元为缓和段, 起,止与直线相接, D,E都取0. 与圆曲线相接, D,E都取(右偏)1/(左偏)-1÷该圆曲线的半径. 与缓和段相接, D,E都取0.任意点的切线方位角C说明：编数据库时,若两曲线相接,例：先编ZH至HY段,运行算出HY点的切线方位角,(即求HY点的坐标,然后终止程序调出C值),然后再以计算出的C编HY至YH段.依此类推.反算桩号时需多次确认所输入的X,Y坐标,程序运算至桩号误差小于允许值时,会同时显示结果:K0-桩号,L0-距离.任意曲线坐标放样正反算程序（fx－4800p）任意曲线坐标放样正反算程序（fx－4850p）。
公路曲线中桩边桩坐标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 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次方，输入程序时切莫输错。
卡西欧4800P计算器测量程序集共10页
福建林业职业技术学院工程系土木工程室高国兴CASIO fx-4500PA计算中桩、边桩坐标程序直线、圆曲线中桩及左、右边桩计算程序fx-4800p=”：C“Xθ=”：D“Yθ=”：E“Wθ=”：N=Ans：R“R=”：L“Lμo=”：K“D=”：M“JμoQ=9θL÷（πR）：S=2RsinQ ：O =C＋Scos（N-Q）◢P=D+Ssin（N-Q）◢W=E-2Q◢X=O+Kcos（N-2Q-M）◢Y=P+Ksin（N-2Q-M）◢注：Xθ=起算点坐标Yθ=起算点坐标Wθ=起算点方位角R=曲线半径，线路前进方向，左+，右－直线R为无穷大=置镜点到待求点距离LμoD左、右边距离=前进方向左（+90），右（－90）JμoO=中桩坐标P=中桩坐标W=待求点方位角X=左、右边桩坐标Y=左、右边桩坐标缓和曲线中桩及左、右边桩计算程序Fx-4800pC“Xθ=”：D“Yθ=”：E“Wθ=”：N=Ans：R“R=”：F“Lθ=”：=”：K“D=”：M“Jμo=”←L“LμoU=L－L^5÷（40R2F2）：V=L^3÷（6RF）－L^7÷（336R^3F^3）：Q=tan－1（V÷U）：S=√（U2+V2）←B=9θL2÷（πRF）：O=C+Scos（N－Q）◢P=D+Ssin（N－Q）◢W=E－B◢X=O+Kcos（N－B－M）◢Y=P+Ksin（N－B－M）◢注：Xθ=起算点坐标Yθ=起算点坐标Wθ=起算点方位角R=曲线半径，线路前进方向，左+，右－Lθ=缓和曲线长=置镜点到待求点距离LμoD左、右待求距离=前进方向左（+90），右（－90）JμoO=中桩坐标P=中桩坐标W=待求点方位角X=左、右边桩坐标Y=左、右边桩坐标放样程序fx-4800pRLI“XZ”：W：“YZ”：G“XH”：H“YH”←B=H－W ←C=G－I ←Prog“P”：U=A←K=√（B2＋C2）←C=X“X－”－I←B=Y“Y－”－W←Prog“P”←L“L”=√（B2＋C2）◢V=A－U：V<θ==>V=V＋36θ⊿V “V=”◢PA=tan－1（B÷C）←B≥θ==> C≥θ==> A=A：≠=> C<θ==> A=A＋18θ⊿⊿≠=> B<θ==> C≥θ==> A=A＋36θ: ≠=>C<θ==>A= A＋18θ⊿⊿⊿←XZ、YZ—置镜点X、Y坐标；XH、YH—后视点坐标；X-、Y-—前视点坐标；输出：R=水平角、L=置镜点至前视点水平距离。
fx4800编程计算器程序二
6. 圆曲线段 文件名：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/π,(左转-，右转 +)。
3.缓和曲线曲线要素： 程序步骤： B=0°1718.87′﹡L/R◢ X=L-L＾3/40/ R2◢ Y=L2/6/ R◢ P=Y-R(1-cos B◢ Q= X-Rsin B◢ T=(R+P)tan(A/2)+Q◢ E=(R+P)(cos(A/2))-1-R◢ Z= R(A-2B)л/180+2L◢ J=2T-Z◢ D=X-Y/tan B 操作过程： HQQXYS→EXE→输入 L 值（即缓和曲线总长）→EXE→输入 R 值（即圆曲线半径） →EXE→得β角度→EXE→得 Xh 值→EXE→得 Yh 值→EXE→得 P 值→EXE→得 Q 值 →EXE→输入 A 角（例 125°31′23.25″）→EXE→得 T 值→EXE→得 E 值→EXE →得 Z 值→EXE→得 J 值→EXE→得 D 值 注：此程序可循环计算。
卡西欧计算器fx4800P
福建林业职业技术学院工程系土木工程室高国兴CASIO fx-4500PA计算中桩、边桩坐标程序直线、圆曲线中桩及左、右边桩计算程序fx-4800pC“Xθ=”：D“Yθ=”：E“Wθ=”：N=Ans：R“R=”：L“Lμo=”：K“D=”：M“Jμo=”：Q=9θL÷（πR）：S=2RsinQ ：O =C＋Scos（N-Q）◢P=D+Ssin（N-Q）◢W=E-2Q◢X=O+Kcos（N-2Q-M）◢Y=P+Ksin（N-2Q-M）◢注：Xθ=起算点坐标Yθ=起算点坐标Wθ=起算点方位角R=曲线半径，线路前进方向，左+，右－直线R为无穷大Lμo=置镜点到待求点距离D左、右边距离Jμo=前进方向左（+90），右（－90）O=中桩坐标P=中桩坐标W=待求点方位角X=左、右边桩坐标Y=左、右边桩坐标缓和曲线中桩及左、右边桩计算程序Fx-4800pC“Xθ=”：D“Yθ=”：E“Wθ=”：N=Ans：R“R=”：F“Lθ=”：L“Lμo=”：K“D=”：M“Jμo=”←U=L－L^5÷（40R2F2）：V=L^3÷（6RF）－L^7÷（336R^3F^3）：Q=tan－1如有你有帮助，请购买下载，谢谢！（V÷U）：S=√（U2+V2）←B=9θL2÷（πRF）：O=C+Scos（N－Q）◢P=D+Ssin（N－Q）◢W=E－B◢X=O+Kcos（N－B－M）◢Y=P+Ksin（N－B－M）◢注：Xθ=起算点坐标Yθ=起算点坐标Wθ=起算点方位角R=曲线半径，线路前进方向，左+，右－Lθ=缓和曲线长Lμo=置镜点到待求点距离D左、右待求距离Jμo=前进方向左（+90），右（－90）O=中桩坐标P=中桩坐标W=待求点方位角X=左、右边桩坐标Y=左、右边桩坐标放样程序fx-4800pRLI“XZ”：W：“YZ”：G“XH”：H“YH”←B=H－W ←C=G－I ←Prog“P”：U=A←K=√（B2＋C2）←。
卡西欧fx-4800P编程
一、本人编写整段曲线坐标程序一则。
可作为参考，从ZH-HY-YH-HZ可一次性计算出中心或切点方向或法线方法等坐标，当只有圆曲线时，L0输入0，只计算单圆周曲线。
示意图无法粘贴。
若有问题请回复。
本程序经调试计算准确。
缓和曲线、圆曲线计算程序A“JDX”：B“JDY”：C“R”：D“L0”：J“HR”：F“T”：G“QD”：H“ZD”：I“Z，Y，-，+”：T“E”：I>0=>Z=1：≠=>Z= -1:△LB1 1：{K，L，M}：K“ZHUANG”：L“L1”：M“L2”：K<G+D=>N“Li”=K-G：Y=1：E=J：GOTO 2：≠=>K≤H-D=> GOTO 3：≠=>N=H-K：E=J+I-180：Y=-1: GOTO 2：△△Lb1 2：O“X0”=A-COS E×F：P“Y0”=B-SIN E×F：Q=N-N^5÷（40×C2×D2）：R=N^3÷（6×C×D）：S=TAN-1（R÷Q）×YZ：U=Y Z N2×90÷（ДCD）+E：K ≤G+D=>U=U-90：≠=>U=U+90:△V=U+90：W“X1”=O+COS（E+S）×√（Q2+R2）+COS U×L+COS V×M▲X“Y1”=P+SIN（E+S）×√（Q2+R2）+SIN U×L+SIN V×M▲GOTO 1：Lb1 3：O=J+Z（0.5×Abs I+90）：P“X0”=Cos O×（C+T）+A：Q“Y0”=Sin O×（C+T）+B：R=O-180-Z（（H-G）-2D）×90÷Д÷C：S=（K-（G+D））×180÷Д÷C：U=R+ZS：W“X1”=P+COS U×（C+ZL）+COS（U+Z×90）×M▲X“Y1”=Q+SIN U×（C+ZL）+ SIN（U+Z×90）×M▲GOTO 1说明：JDX：交点X坐标 JDY：交点Y坐标R：半径 L0：缓和曲线长 HR：前进方位角T：切线长 QD：起点桩号 ZD：终点桩号Z，Y，-，+：偏角，左为负，右为正E：外矢距ZHUANG：待求点桩号L1：待求点法线方向距离，左侧为正，右侧为负，求中桩坐标时为0L2：待求点切线方向距离，左侧为正，右侧为负，求中桩坐标时为0本程序用于计算缓和曲线圆曲线缓和曲线缓和曲线为对称型缓和曲线，若无缓和曲线，计算单圆曲线时L0输入0卡西欧4800计算全站仪坐标程序！这个是我月个朋友编写的“卡西欧4800计算全站仪坐标程序”一共是5个程序。
CASIO fx4800P曲线任意里程中边桩坐标正反算程序
Y0 ？线元起点的Y坐标
S0 ？线元起点里程
F0 ？线元起点切线方位角
LS ？线元长度
R0 ？线元起点曲率半径
RN ？线元止点曲率半径
Q ？ 线 元左右偏标志(左偏Q=-1，右偏Q=1，直线段Q=0）
S=700 Z= 5 计算得 XS=19823.25398 YS=28503.95084
S=780 Z=-5 计算得 XS=19785.25749 YS=28575.02270
S=780 Z=0 计算得 XS=19781.15561 YS=28572.16358
X=19742.686 Y=28652.914 计算得 S=870.0003175 Z=- 0.00041814
X=19737.837 Y=28651.697 计算得 S=870.0002748 Z= 4.999808656
X=19741.5912 Y=28722.0580 计算得 S=939.9999786 Z= -5.123024937 80.0000016 Z= 4.99999578
X=19747.536 Y=28654.131 计算得 S=870.0001137 Z= -4.99941049
919.527 19736.072 28701.893 91 17 30.63 80.285 221.75 9579.228 -1
999.812 19744.038 28781.659 80 40 50.00 100.000 1E45 1E45 0
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+
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
999.812 19744.038 28781.659 80 40 50.00 100.000 1E45 1E45 0
1、正算
（注意：略去计算方式及线元要素输入，请自行根据所求点所在的线元输入线元
要素）
S=700 Z=-5 计算得 XS=19831.41785 YS=28509.72590
X ？反算时所求点的X坐标
Y ？反算时所求点的Y坐标
显示部分：
XS=××× 正算时，计算得出的所求点的X坐标
YS=××× 正算时，计算得出的所求点的Y坐标
S=××× 反算时，计算得出的所求点的里程
SUB2)序构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线
元要素（止点曲
率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。另
外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲
X=19737.837 Y=28651.697 计算得 S=870.0002748 Z= 4.999808656
X=19741.5912 Y=28722.0580 计算得 S=939.9999786 Z= -5.123024937
X=19823.25398 Y=28503.95084 计算得 S=699.9999985 Z= 5.000003137
X=19785.25749 Y=28575.02270 计算得 S=780.0000035 Z= -5 .000001663
(6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的
值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等
于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。
2、输入与显示说明
输入部分：
1. SZ => XY
806.748 19766.566 28594.574 120 25 54.07 112.779 221.75 221.75 -1
919.527 19736.072 28701.893 91 17 30.63 80.285 221.75 9579.228 -1
三、使用说明
1、规定
(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，
Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。
(2) 当所求点位于中线时，Z=0；当位于中线左铡时，Z取负值；当位于中线中线右
侧时，Z取正值。
(3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。
S=940 Z=0 计算得 XS=19736.47687 YS=28722.35642
S=940 Z= 3.009 计算得 XS=19733.47298 YS=28722.53168
(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
2、 反算
X=19831.418 Y=28509.726 计算得 S=699.9999974 Z= -5 .00018164
X=19827.336 Y=28506.838 计算得 S=699.9996493 Z= 0.000145136
S=870 Z=0 计算得 XS=19742.68648 YS=28652.91379
S=870 Z= 5 计算得 XS=19737.83688 YS=28651.69668
S=940 Z=-5.123 计算得 XS=19741. 59118 YS=28722.05802
S0 X0 Y0 F0 LS R0 RN Q
500.000 19942.837 28343.561 125 16 31.00 269.256 1E45 1E45 0
769.256 19787.340 28563.378 125 16 31.00 37.492 1E45 221.75 -1
X=19781.15561 Y=28572.16358 计算得 S=780.0000025 Z=- 0.000002979
X=19777.05373 Y=28569.30446 计算得 S=780.0000016 Z= 4.99999578
X=19736.4769 Y=28722.3564 计算得 S=939.9999862 Z=- 0.000027710
X=19733.4730 Y=28722.5317 计算得 S=940.0000238 Z= 3.00898694
线任意里程中边桩坐标进行正反算。本程序也可以在CASIO fx-4500P计算器及 CASIO fx-4850P计算器上运行。
二、源程序
1.主程序(TYQXjs)
"1.SZ => XY"："2.XY => SZ"：N：U"X0"：V"Y0"：O"S0"：G"F0"：H"LS"：P"R0"：R"
2. XY = > SZ
N ? 选择计算方式，输入1表示进行由里程、边距计算坐标 ；输入2表示由坐标反算
里程和边距。
X0 ？线元起点的X坐标
Y0 ？线元起点的Y坐标
S0 ？线元起点里程
F0 ？线元起点切线方位角
Z=××× 反算时，计算得出的所求点的边距
四、算例
某匝道的由五段线元（直线+完整缓和曲线+圆曲线+非完整缓和曲线+直线）组
成，各段线元的要素（起点里程S0、起点坐标X0 Y0、起点切线方位角F0、线元长度
LS、起点曲率半径R0、止点曲率半径RN、线 元左右偏标志Q）如下：
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
CASIO fx4800P曲线任意里程中边桩坐标正反算程序
作者：lq52搜集 人气： 来源：网络 【字体：大 中 小】【打印】【关闭】
CASIO fx4800P曲线任意里程中边桩坐标正反算程序
曲线任意里程中边桩坐标正反算(CASIO fx-4800P计算器)程序
一、程序功能
本程序由一个主程序(TYQXjs)和两个子程——正算子程序(SUB1)、反算子程序(
(4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆
弧的半径。
(5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45
次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半
径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。
S=700 Z=0 计算得 XS=19827.33592 YS=28506.83837
S=700 Z= 5 计算得 XS=19823.25398 YS=28503.95084
S=780 Z=-5 计算得 XS=19785.25749 YS=28575.02270
S=780 Z=0 计算得 XS=19781.15561 YS=28572.16358
S=780 Z= 5 计算得 XS=19777.05373 YS=28569.30446
S=870 Z=-5 计算得 XS=19747.53609 YS=28654.13091
Z"Z"=Z◢
Goto 2
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
X=19747.536 Y=28654.131 计算得 S=870.0001137 Z= -4.99941049
X=19742.686 Y=28652.914 计算得 S=870.0003175 Z=- 0.00041814
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◢
Goto 1←┘
Lbl 2：{XY}：XY：I=X：J=Y：Prog "SUB2"：S"S"=O+W◢
LS ？线元长度
R0 ？线元起点曲率半径
RN ？线元止点曲率半径
Q ？ 线 元左右偏标志(左偏Q=-1，右偏Q=1，直线段Q=0）