“第1章 绪论”CASIOfx-4800P计算程序
卡西欧fx-4800P计算器程序集(小刚)
计算器对应程序名称：{ZDZB}（转点坐标）A“XZ”B“YZ”C“XH”D“YH”：{RL}I=1：J=1：Pol（（C－A），（D－B））：F=R＋JRec（L，F）X“X”=I＋A◢Y“Y”=J＋B◢注：用于已知两点坐标，及待测点与测站点的连线与已知边（两已知点的连线）的水平夹角和待测点至测站点的水平距离求待测点坐标。
XZ、YZ—测站点坐标；XH、YH—后视点坐标；R—水平角（顺时针为正，逆时针为负）；L—待测点至测站点水平距离。
计算器对应程序名称：{FWJ}(方位角)X“X1”：Y “Y1”：U“X2”：V“Y2”G=U-X:H=V-YPol(G,H):I◢J<0==>J=+360⊿J“J0=”注：用于已知两点坐标，求两点之间的方位角和距离。
X1、Y1—前视（后视）点坐标；X2、Y2—测站点坐标；L—测站点至前视（后视）点距离；J0—方位角(按SHIFT键再按°′″键转换为60制角度)。
任意多边形面积计算程序fx-4800p（三）计算器对应程序名称：{MJ}Lbl 0:M“N”：N“E”：E=M：F=N：S=0Lbl 1:{XY}:S=S-(XN-YM)/2:X=E==>Y=F==>“㎡”:S ◢Goto O ⊿⊿M=X:N=Y:Goto 1注：用于计算任意多边形的面积。
N-起始点X坐标（闭合点X坐标）E-起始点Y坐标（闭合点Y坐标）X-前进点X坐标Y-前进点Y坐标S-多边形面积（取绝对值）应顺序依次输入点坐标，输完所有点坐标后应再次输入起点坐标，这样就形成了多边形的闭合，程序将自动计算出多边形的面积。
极坐标测设建筑物（放样）计算程序fx-4800p（四）计算器对应程序名称：{ZBCS}E:F:K:LG=K-E:H=L-FPOL(G,H):J<0==>B=J+360：≠==>B=J⊿Lbl 0{XY}:X=0==>G o t o E⊿G=X-E:H=Y-FPol(G,H):I◢J<0==>J=J+360⊿C=J-B:C<0==>C=C+360⊿C◢Goto 0Lbl E注：用于计算已知方向边AB（已知点A、B两点的连线）与A（测站点）点至任意放样点i方向边（测站点与待放样点的连线）的水平夹角Bi和水平距离。
卡西欧计算器坐标计算程序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点桩号。
卡西欧计算器程序fx-4800p
Casio Fx-4800P计算器测量计算程序抚顺市市政建设(集团)有限公司第一工程处：张和光Casio Fx-4800P一、整机概况Fx-4800P计算器是日本Casio公司的产品，它除了具备函数型计算器的全部功能之外，还具备如下特殊功能：1、二进制和八进制字组显示；2、二进制、八进制、十进制和十六进制的转换和加、减、乘、除运算。
3、进行f（x）函数的输入与积分的计算；4、进行标准差和回归计算。
Fx-4800P计算器具有同时显示运算表达式和计算结果的四行式显示的特点。
该计算器具有记存4500个程序步的记忆储存装置，计算器储存了程序之后，只要您输入数据，再按一下EXE键，计算器就会用数据进行程序执行,对于随数据而变化的反复计算的运用十分方便。
Fx-4800P计算器具备条件指令转移和非条件指令转移的能力，还具备逻辑判断能力，因此它可以进行测量工作中比较复杂的计算，当然其它专业的计算工作也同样可以解决。
（1）Fx-4800P计算器操作方法比较简单、易学，程序编排灵活、易懂，可称之谓各专业技术人员的物美、价廉、得心应手的计算工具。
每一台Fx-4800P计算器都具备一本说明书，该计算器各部键的功能以及操作方法，只要认真按照说明书各章节学习，很快就能掌握和运用，现将我在生产实践中应用的部分计算程序介绍给测绘界的同行以便共同探讨。
（2）二、程序库1．坐标反算：程序步骤：C”X1”:D”Y1”:E”X2”:F”Y2”:Fixm:Pol(E-C ,F-D:I”S1-2=”◢J≤O=>J=J+360△J”A1-2=”操作过程：ZBFS→EXE→输入X1值→EXE→输入Y1值→EXE→输入X2值→EXE→Y2→EXE→EXE→EXES1-2：计算得出的距离；A1-2：计算得出的角度。
（按shift°′″转换为60进制的角度）注：此程序可循环计算。
（3）2．坐标正算：程序步骤：C”X1”:D”Y1”:S”S1-2”:R”A1-2”:Fixm:X”X 2”=C+Rec(S,R◢Y” Y2”=D+J操作过程：ZBZS→EXE→输入X1值→EXE→输入Y1值→EXE→输入S1-2距离→EXE→输入A1-2角度（例125°31′23.25″）→EXE→EXE注：此程序可循环计算。
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.编制程序过程中已经考虑了零溢出问题。
FX4800P测量程序集
CASIO fx 系列编程计算器测量程序集编写: 王劲松第一部分：CASIO fx-4800P 编程计算器测量程序集一．SHI JU CE TU （视距测图）1． 用途：该程序用于“经纬仪视距法碎部测量”时，计算测站点至碎部点间的平距及碎部点的高程。
2． 程序清单：C “H0”：I ：A “UP ”：B “DOWN ”：V “MIDDLE ”：L ：E=90-L ↙D=100Abs （A-B ）（cosE ）2◢ H=（C+I ）+DtanE-V3． 程序说明：H0——测站点高程；I ——测站仪器高；UP ——塔尺上丝读数；DOWN ——塔尺下丝读数；MIDDLE ——塔尺中丝读数；L ——经纬仪竖盘读数盘左读数；D ——测站点至碎部点间平距；H ——碎部点的高程。
5．该程序若在fx-3950中的程序清单为：？→M ：？→A ：？→B ：？→C ：？→D ：100（A-B ）（cos(90-D)）2→X ◢ M+Xtan(90-D)-C →Y说明：M ——将测站高程H0+仪器高I ；A ——上丝读数；B ——下丝读数；C ——中丝读数；D ——竖盘盘左读数L 。
X ——测站点至碎部点的平距；Y ——碎部点的高程。
二．FAN SUAN （反算）1．用途：该程序是“坐标反算公式”程序，用于根据两点1、2的平面坐标（X1，Y1）、（X2，Y2），计算两点间的距离D12及坐标方位角A12。
2． 程序清单：G “X1”：B “Y1”：C “X2”：H “Y2”：E=C-G ：F=H-B ↙Pol （E ，F ）：I “D12”◢ ∆+=⇒<3600J J J J “A12”3． 程序说明：X1，Y1——点1的平面坐标；X2，Y2——点2的平面坐标； D12——点1至点2间平距；A12——点1至点2坐标方位角。
由于编程时用到了极坐标函数POL （），故在程序运行中，会出现I=？XXX （数字）及J=？XXX （数字）这是函数POL计算出的两个量，不能改变，直接回车默认即可。
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.高程计算。
CASiO fx-4800计算器程序
CASiO fx-4800计算器程序1、CE ZHANFix3:Prog"X":U=X:V=Y:Prog"X":X=X:Y=Y:Fixm:Pol(X-U,Y-V:Pause0:W=J:J<0=>J=J+360△IntJ+.01(.6Fr ac 60J+Int(60 FracJ:Fix 4:J"D,F"= Ans2、QU XIANFix3:L"K+"D"B÷2":Fixm:AbsD<3E 3=>Prog"Q":≠>X=L:Y=D △X:Pause0:Y"N,E="◢Pol ( X-U ，Y-V:Pause0: J <0=>J=J+360△IntJ+.01( .6 Frac 60J+Int(60 FracJ:Fix 4:J"D,FWJ"=Ans3、X{NXY}:N=控制点编号=>X=控制点X 坐标：Y=控制点Y 坐标△N=下个控制点编号=>X=下个控制点X 坐标：Y=下个控制点Y 坐标△……X:Pause0:Y4、QProg"SJ":Prog"P":C=L-H:C>.5B=>K=B-C:≠>K=C △K>0=>K≤S=>A=RS:X=K -K 5÷40A 2:Y=K 3÷6A -K 7÷336A 3:≠>A=180(K-S)÷πR+M:Y=R(1-cosA)+W:X=RsinA+Z ［1△Pol(X,Y:≠>I=K:J=0△E-F>0=>C<.5B=>X=N+Rec(I,J+F:Y=J+Q:Prog"J":J=Ans+F:≠>J=-J:Prog"K":C=B-C:Prog"J":J=E-Ans:△≠>C<.5B=>J=-J:X=N+Rec(I,J+F:Y=J+Q:Prog"J":J=F- Ans: ≠>Prog"K":C=B-C:Prog"J":J=Ans+E: △△X=X+Rec(D,J+90:Y+J5、KX=P+Rec(I,J+180+E:Y=J+O6、JC≤0=>0: C≤S=>90C 2÷πA: ≠>A △△△7、PE=F-Z:G=H+B:W=S 2÷24R: Z ［1］=S-S 3÷40R 2-Rsin(180S÷2πR:M =180S÷2πR ：Rec （T,F-180：N=I+X:Q=J+Y: Rec （T,E:P=I:+X:O =Y+J8、SJL ≤缓直点桩号=>B=曲线长度:F=计算方位角:H=直缓点桩号:T=切线长度：X=交点X 坐标：Y=交点Y 坐标： S=缓和曲线长度:R=曲线半径：Z=转角（左转输正，右转输负）：≠>L≤下段直缓点桩号=>B=曲线长度：F=计算方位角:H=直缓点桩号：T=切线长度：X=交点X 坐标：Y=交点Y 坐标： S=缓和曲线长度:R=曲线半径：Z=转角值: ≠>L≤下段直缓点桩号……△△△ (注：有多少弯道就加多少个△)。
4800程序
CASIO（fx-4800p）程序单源公式：水平距：S= √（X［ZD］－X［QD］）+（Y［ZD］－Y［QD］）方位角：A=arctg（Y［ZD］－Y［QD］）/（X［ZD］－X［QD］）程序：■FANG WEI JIAO ［文件名］Lbl 0{M}M″X［ZD］″{P}P″Y［ZD］″N″X［QD］″Q″Y［QD］″Pol（（M－N），（P－Q））▲J＜0=＞A″TA″=J+360▲≠＞A″TA″=J▲△Goto 0说明：输入后视点X［ZD］、Y［ZD］坐标；输入测站点X［QD］、Y［QD］坐标；显示（Disp）水平距（r =…）显示（Disp）方位角（θ=…）若为负值则加360度；显示（Disp）方位角（TA=…）二、文件名：直线上任意点坐标计算源公式：横坐标：X=X［QD］+ScosA纵坐标：Y=Y［QD］+SsinA程序：ZX ［文件名］■Lbl 1B″X［QD］″C″Y［QD］″E″QD ZHUANG HAO″{D}D″ZHUANG HAO″X=B+（D－E）cosA▲Y=C+（D－E）sinA▲J=A▲Goto 1说明：输入起算点X［QD］、Y［QD］坐标；输入起算点桩号；（QD ZHUANG HAO）输入所求点桩号；（ZHUANG HAO）输入直线方位角；（A）显示（Disp）所求点横坐标（X）；显示（Disp）所求点纵坐标（Y）；显示（Disp）所求点切线方位角（J）。
三、文件名：（左偏）圆曲线任意点坐标计算源公式：横坐标：X=X［ZY］+2Rsin（90L/（πR）cos（A－90L/（πR）））纵坐标：Y=Y［ZY］+2Rsin（90L/（πR）sin（A－90L/（πR）））程序：YL ［文件名］■Lbl 2B″X［ZY］″C″Y［ZY］″{D}D″ZHUANG HAO″{E}E″ZHI YUAN DIAN″F=90（D－E）÷（πR）X=B+2RsinFcos（A－F）▲Y=C+2RsinFsin（A－F）▲J=A－180（D－E）÷（πR）▲Goto 2说明：输入直圆点X［ZY］、Y［ZY］坐标；输入所求点桩号；（ZHUANG HAO）输入直圆点桩号；（ZHI YUAN DIAN）输入圆曲线半径；（R）输入直圆点~交点的方位角（A）显示（Disp）所求点横坐标（X）；显示（Disp）所求点纵坐标（Y）；显示（Disp）所求点切线方位角（J）。
卡西欧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）。
CASIO-fx4800程序
CASIO-fx4800程序本人自网络上进行整理，谢谢大家支持！！！！CASIO-fx4800p 1 KWX Lbl1:{ABC}:Z[5]=A:Z[6]=B:Z[7]=C:Lbl2:Prog“ZBJX”:Abs(Z[5]-X)<0.001ÞG oto 4∆ Pol(Z[5]-X,Z[6]-Y):K=K+Icos(J-V):U=U-Isin(J-V):Goto2:Lbl4:Prog“SQX”:A=Z[7]-H:A$0ÞL=AbsU -1.5A-13:¹ÞL=AbsU-1.0A-15∆ Ky A+0.75y Uy Ay Goto 1¿ 2 ZBJX Fixm:{K}:Lbl 0:{K}:K:K# $ÞW:C: D:A:B:R:S∆ P=S2/24R-S4/2688R3:Q=S/2-S3/240R2:T=(R+P)tg(AbsA/2)+Q:L=pRAbsA/180:E=W-T:F= E+L¿ K#EÞL=W-K:X=C+Lcos(B-180):Y=D+Lsin(B-180):V=B:Goto 1∆ K#E+SÞL=K-E:X=L-L5/40R2S2: Y=L3-6RS-L7/336R3S3:Pol(T-X,Y):X=C+Icos(B-180-AbsAJ/A):Y=D+Isin(B-180- AbsAJ/A):V=B+90L2Ab sA/ApRS:Goto 1∆ K#FÞL=K-E-S:I=(S/2+L)×180/pR:X=T-RsinI-Q:Y=R(1-cosI)+P:V=IAbsA/A+B:Pol(X, Y):X=C+Icos(B-180-JAbsA/A):Y=D+ Isin(B-180-JAbsA/A):Goto 1∆ K#F+SÞL=F+S-K:X=T-L+L5/40R2S 2:Y=L3/6RS-L7/336R3S3:Pol(X,Y):X=C+Icos(B+A+JAbsA/A):Y=D+ Isin(B+A+JAbsA/A):V=B+A-90L2A bsA/ApRS:Goto 1∆ K>F+SÞL=K+T-S-F:X=C+Lcos(A+B):Y=D+Lsin(A+B):V=A+B:Goto 1∆ Lbl 1:X=X+ Ucos(V-090):Y=Y+Usin(V-090)¿ O=V¿ 3 SQX Fixm:{K}:K#$ÞG:C:R:E:F∆ W=F-E:T=Abs(RW/2):K<G-TÞH=C-E(G-K)∆ K$G-TÞH=C-E(G-K)+(K+T-G)2AbsW/2RW∆ K$GÞH=C-F(G-K)+(K-G-T)2AbsW/2RW∆K$G+TÞH=C-F(G-K)∆¿ 4 HGHXP {KU}:K:U=U:Prog“ZBJX”:Prog:“SQX”:Xy Yy H¿ 备注：1(KWX)中，A、B、C为某点的坐标X，Y，Z；K 桩号；L 当L=0时处于开挖线位置；A 高程；A+0.75 测点高程+路面厚度；U 距中线距离。
CASIO4800P计算器坐标计算程序
坐标放样程序Defm 60：A”JX”:B”JY”:Q:R:L”L1”:M”L2”:T”T1”: C”T2”:E”ZH”:D=E+L:O”HZ”:I=O-M:F:U :Fix3：回车Lbi 0:{K}:K≤E=〉G0t0 1：≠=〉G0t0 2⊿回车Lbi1：{VJZ}:Z[1]=A-T CosF:Z[2]=B-TSinF: Z[3]=E-K: Z[4]=F+VJ: Z[20]= Z[1]- Z[3] CosF:Z[21]=Z[2]-Z[3]SinF：X=Z[20]+ZCos （Z[4]）：Y=Z[21] +ZSin（Z[4]）：X“X=”▲Y“Y=”▲{N}：N＜1=〉Prog“FMH”:≠=> Goto 0回车Lbi 2：K≤D=> Goto 3: ≠=> Goto 4⊿回车Lbi3:{VJZ}:Z[5]=A-TCosF:Z[6]=B-TSinF: H=K-E: Z[7]=H-（(H＾5)÷（40R＾2L＾2)）-（(H＾13)÷（599040R＾6L＾6)）+（(H＾9)÷（3456R＾4L＾4)）:W=F+U(90H＾2÷（RL ∏）)+VJ:Z[22]=Z[5]+Z[7] ÷Cos(（30H＾2）÷（RL∏)） *Cos(F+U(30H^2)÷(RL∏)): Z[23]= Z[6]+ Z[7] ÷Cos((30H^2)÷(RL∏)) *Sin(F+U(30H^2)÷(RL∏)):回车X= Z[22]+Z CosW:Y=Z[23]+Z SinW: 回车X“X=”▲Y“Y=”▲Prog“AA”:回车{N}：N＜1=〉Prog“FMH”:≠=> Goto 0回车Lbi 4:K≤I=> Goto 5:≠=> Goto 6⊿回车Lbi5:{VJZ}:Z[8]=A+TCos(F+180):Z[9]=B+ TSin(F+180):Z[10]=L-(L^3)÷(40R^2)-(L^7)÷(599040R^6)+(L^5)÷(3456R^4): Z[11]=(30L)÷(R∏): Z[12]=F+U(30L)÷(R∏):P=K-D: Z[13]=F+90U（L+P）÷(R∏): Z[14]=(90P)÷(R∏):W=F+U(90L)÷(R∏)+(180UP)÷(R ∏)+VJ: 回车Z[24]= Z[8]+ Z[10]÷Cos(Z[11]) *Cos(Z[12])+2RSin(Z[14]) *Cos(Z[13]): Z[25]= Z[9]+ Z[10] ÷Cos(Z[11]) *Sin(Z[12])+2RSin(Z[14]) Sin(Z[13]):X=Z[24]+Z CosW:Y= Z[25]+Z SinW: X“X=”▲Y“Y=”▲Prog“BB”:{N}：回车N〈1=>Prog“FMH”:≠=> Goto 0⊿回车Lbi 6:K≤O=> Goto 7: ≠=> Goto 8⊿回车Lbi 7:{VJZ}:P=F+Q:Z[15]=A+C CosP: Z[16]=B+C SinP:H=O-K: Z[17]=H-(H∧5)÷(40R∧2M∧2)-(H∧13)÷(599040R∧6M∧6)+(H∧9)÷(3456R∧4M∧4): Z[18]=(30H∧2) ÷(RM ∏): Z[19]=P+180-U(30H^2)÷(RM ∏):W=P+180-U((90H∧2) ÷(RM∏))-VJ:Z[26]= Z[15]+ Z[17] ÷Cos (Z[18]) *Cos(Z[19]): Z[27]= Z[16]+ Z[17] ÷Cos( Z[18])* Sin (Z[19]):X= Z[26]+ Z CosW:Y= Z[27]+Z SinW:X“X=”▲Y“Y=”▲Prog“CC”:回车{N}：N〈1=>Prog“FMH”:≠=> Goto 0回车Lbi 8:K≤O=> Goto7: ≠=> Goto9⊿回车Lbi9:{VJZ}:P=F+Q:W=P+VJ:Z[28]=A+(C+K-O) CosP: Z[29]=B+(C+K-O) SinP:X= Z[28]+Z CosW:Y= Z[29]+Z SinW:X“X=”▲Y“X=”▲Prog“DD”:{N}：N〈1=>Prog“FM H”:≠=> Goto 0⊿2、文件名“ FMH”Lbi1:{S}:S<0=> Goto2: ≠=> Goto3⊿回车Lbi2{J}:Z[30]=X- Z[20]: Z[31]=Y- Z[21]: Z[32]= Z[31] ÷Z: Z[33]= Z[30] ÷Z:V= Z[20]+JZ[33]-SZ[32]:K= Z[21]+ JZ[32]+S Z[33]:V“X=”▲K“Y=”▲Goto13、文件名“AA”Z[20]= Z[22]: Z[21]= Z[23]4、文件名“BB”Z[20]= Z[24]: Z[21]= Z[25]5、文件名“CC”Z[20]= Z[26]: Z[21]= Z[27]6、文件名“DD”Z[20]= Z[28]: Z[21]= Z[29]道路前进方向操作：A:B:坐标:Q转角（左负，右正）：R：半径：L1：第一缓和曲线长度（L）：L2:第二缓和曲线长度(M):T1:第一起线长度（T）:T2: 第二起线长度(C):ZH:直缓点里程(E):HZ:缓直点里程(O):F：方位角：U：转向系数（左转负，U=-1,右正，U=1）:K:待求点里程：V:方向系数（左取－1，右取1）：J：边桩与道路中线的夹角:Z:边距：在左边上时S取正J取0左边下时S取负J 为距中宽度，在右边上时S取负J取距宽度中，在右边下时S取正J取0。
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次方，输入程序时切莫输错。
卡西欧计算器程序fx-4800p
Casio Fx-4800P计算器测量计算程序抚顺市市政建设(集团)有限公司第一工程处：张和光Casio Fx-4800P一、整机概况Fx-4800P计算器是日本Casio公司的产品，它除了具备函数型计算器的全部功能之外，还具备如下特殊功能：1、二进制和八进制字组显示；2、二进制、八进制、十进制和十六进制的转换和加、减、乘、除运算。
3、进行f（x）函数的输入与积分的计算；4、进行标准差和回归计算。
Fx-4800P计算器具有同时显示运算表达式和计算结果的四行式显示的特点。
该计算器具有记存4500个程序步的记忆储存装置，计算器储存了程序之后，只要您输入数据，再按一下EXE键，计算器就会用数据进行程序执行,对于随数据而变化的反复计算的运用十分方便。
Fx-4800P计算器具备条件指令转移和非条件指令转移的能力，还具备逻辑判断能力，因此它可以进行测量工作中比较复杂的计算，当然其它专业的计算工作也同样可以解决。
（1）Fx-4800P计算器操作方法比较简单、易学，程序编排灵活、易懂，可称之谓各专业技术人员的物美、价廉、得心应手的计算工具。
每一台Fx-4800P计算器都具备一本说明书，该计算器各部键的功能以及操作方法，只要认真按照说明书各章节学习，很快就能掌握和运用，现将我在生产实践中应用的部分计算程序介绍给测绘界的同行以便共同探讨。
（2）二、程序库1．坐标反算：程序步骤：C”X1”:D”Y1”:E”X2”:F”Y2”:Fixm:Pol(E-C ,F-D:I”S1-2=”◢J≤O=>J=J+360△J”A1-2=”操作过程：ZBFS→EXE→输入X1值→EXE→输入Y1值→EXE→输入X2值→EXE→Y2→EXE→EXE→EXES1-2：计算得出的距离；A1-2：计算得出的角度。
（按shift°′″转换为60进制的角度）注：此程序可循环计算。
（3）2．坐标正算：程序步骤：C”X1”:D”Y1”:S”S1-2”:R”A1-2”:Fixm:X”X 2”=C+Rec(S,R◢Y” Y2”=D+J操作过程：ZBZS→EXE→输入X1值→EXE→输入Y1值→EXE→输入S1-2距离→EXE→输入A1-2角度（例125°31′23.25″）→EXE→EXE注：此程序可循环计算。
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）。
卡西欧4800公路测量程序
卡西欧4800公路测量程序第一章：基本符号1、X0——始点横坐标2、Y0——始点纵坐标3、X1——终点横坐标4、Y1——终点纵坐标5、D——始点到终点的水平距离6、F——始点到终点的方位角7、X(P1)——P1点的横坐标8、Y(P1)——P1点的纵坐标9、X(P2)——P2点的横坐标10、Y(P2)——P2点的纵坐标11、D(OP1)——置仪点到P1的水平距离12、D(OP2)——置仪点到P2的水平距离13、F(OP2)——置仪点到P2的方位角14、X(00)——置仪点的横坐标15、Y(00)——置仪点的纵坐标16、JD，JD(AB)——选择输入1或者2（1）JD——实交（2）JD(AB)——虚交17、IR、IL——选择输入1或者2（1）IR——路线右转（2）IL——路线左转19、(1) E-R; (2) T-R——选择输入(1)或(2)（1）E-R——继续E长控制半径R（2）T-R——T长控制半径R19、(1) E-R; (2) END——选择输入(1)或(2)（1）E-R——继续E长控制半径R测曲中（2）END——结束20、(1)S1，(2)S2，(3)R2——选择输入(1)或(2)或(3)(反复曲线S型)（1）S1——进入第一条曲线测量（2）S2——进入第二条曲线测量（3）R2——这是由第二条曲线测设后，进入第一条曲线测，在进入第二条曲线直接测量曲中21、(1)R1，(2)R2，(3)END——选择输入(1)或(2)或(3)(反复曲线S型)（1）R1——如对曲线一、二的测设皆不满意，可输入1，进入第一条曲线的测设，结束后输入3，进入第二条曲线的测设（2）R2——(2)R2:如对曲线二的邪恶设不满意，可输入2，再测曲中（3）END——结束22、X(QZ)——曲中的横坐标23、Y(QZ)——曲中的纵坐标24、D(QZ)——置仪点到曲中的水平距离25、F(QZ)——置仪点到曲中的方位角26、K——输入桩号27、W(L)——左侧边桩距离中桩的水平距离28、W(R)——右侧边桩距离中桩的水平距离29、X(L)——左侧边的桩横坐标30、Y(L)——左侧边的桩纵坐标31、H(L)——左侧边的高程32、D(L)——置仪点到左侧边桩的水平距离33、F(L)——置仪点到左侧边桩的方位角34、X(Z)——中桩的横坐标35、Y(Z)——中桩的纵坐标36、H(Z)——中桩的高程37、D(Z)——置仪点到中桩的水平距离38、F(Z)——置仪点到中桩的方位角39、X(R)——右侧边桩的横坐标40、Y(R)——右侧边桩的纵坐标41、H(R)——右侧边桩的高程42、D(R)——置仪点到右侧边桩的水平距离43、F(R)——置仪点到右侧边桩的方位角44、F1——第一条线路的方位角45、F2——第二条线路的方位角46、JA——虚交JDA转角47、JB——虚交JDB转角48、ａ——路线转角49、R——平、竖曲线半径50、R1——第一天曲线半径51、R2——第二条曲线半径52、L S——缓和曲线长度53、L S1——第一缓和曲线长度54、L S2——第二缓和曲线长度55、T——平、竖曲线切线长度56、T1——第一平曲线切线长度57、T2——第二平曲线切线长度62、L——平、竖曲线长度63、L1——第一平曲线长度64、L2——第二平曲线长度65、E——平、竖曲线外距66、E1——第一平曲线外距67、E2——第二平曲线外距68、LY——圆曲线长度69、L Y1——第一圆曲线长度70、L Y2——第二圆曲线长度72、K(JD) ——交点桩号73、AB——虚交JDA—JDB的水平距离74、ZH——直缓点桩号75、HY——缓圆点桩号76、QZ——曲中点桩号77、YH——圆缓点桩号78、HZ——缓直点桩号79、GQ——公切点桩号80、N°——右角顺时针第二章：路线测量与道路施工放样计算程序第一节：已知两点坐标求距离和方位角计算程序【例题1】如图1-1，已知A点坐标A(100,100)， XB点坐标B(-100,-100)。
卡西欧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个程序。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
“第1章 绪论”CASIOfx-4800P 计算程序§1.1 参考椭球元素计算测量中，高程的基准面是大地水准面，地表某点距离大地水准面的垂直距离称为该点的高程H ；大地坐标系的基准是参考椭球，用大地经度L 和大地纬度B 表示点在参考椭球表面的位置，用大地方位角A 表示球面上边长的方向。
定义参考椭球大小的元素值为：a ——长半轴b ——短半轴f ——扁率 aba f −=e ——第一偏心率 2222a b a e −=e ′ ——第二偏心率 2222b b a e −=c ——极曲率半径 ba c 2=我国采用过的两个参考椭球元素值以及GPS 测量使用的参考椭球元素值列于表1-1。
表1-1 参考椭球元素值坐标系名称 1954北京坐标系 1980西安坐标系 WGS-84坐标系(GPS 用)参考椭球名称 克拉索夫斯基椭球 IUGG —1975椭球 IUGG —1979椭球a 6378245 6378140 6378137b6356863.01877 6356755.28815 6356752.31424066 f 1:298.3 1:298.257 1:298.2572235632e 0.006693421622966 0.006694384999588 0.006694379990132e ′ 0.006738525414684 0.006739501819473 0.00673949674227 c6399698.90178 6399596.65199 6399593.62601(1) 数学模型通过参考椭球面上任一点的法线，可以作无数个法截面，法截面与参考椭球面的交线称为法截线，球面上不同纬度及方向法截弧线的曲率半径是不相同的。
1) 子午圈曲率半径M通过参考椭球旋转轴南北极的法截面(也称子午面)与参考椭球面相截形成的闭合圈称为子午圈，其曲率半径M 的计算公式为：3Vc M =(1-1) 式中21η+=V(1-2) B e cos ′=η(1-3)2) 卯酉圈曲率半径N与该点子午面相垂直的法截面与参考椭球面相截形成的闭合圈称为卯酉圈，其曲率半径N 的计算公式为：Vc N =(1-4)3) 平均曲率半径R2V c MN R == (1-5)由于V >1，所以有M R N >>。
4) 任意法截弧的曲率半径A R方位角为A 方向法截弧的曲率半径A R 的计算公式为ANR A 22cos 1η+=(1-6)当A =0时，M R A =；当A =90时，N R A =。
(2) 程序与案例分别在1954北京坐标系、1980西安坐标系和WGS-84坐标系中计算参考椭球表面上任意点的M 、N 、R 及A R 。
1) 变量对照表数学模型变量 fx-4800P 变量 单位 注释c Cm 极曲率半径 2e ′ E第二偏心率 B B°′″大地纬度 A A°′″大地方位角 η U中间参数 V V中间参数 M Mm 子午圈曲率半径 N Nm 卯酉圈曲率半径 R Rm 平均曲率半径 A R P m 任意方向法截弧曲率半径 2) 程序程序名：MNR Lbl 0 {Z}Z=54C=6399698.90178: E=0.006738525414684 : Z=80C=6399596.65199:E=0.006739501819473 : Z=84C=6399593.62601: E=0.00673949674227 : Goto 0:{BA}U=EcosB: V=(1+U 2)M=C ÷V^3N=C ÷VR=C ÷V 2P=N ÷(1+(UcosA)2)3) 操作步骤设=B 22°32′50″，=A 153°19′18″，分别计算其在1954北京坐标系、1980西安坐标系及WGS-84坐标系中的M 、N 、R 及A R 值。
将上述程序以MNR 的文件名输入计算器后，按键 5 2及选择程序MNR ，按键，屏幕提示及操作步骤如下：步骤 显示 按键 注释1 Z? ××× 54 大地坐标系2 B? ××× 223250 大地纬度3 M=6344916.774子午圈曲率半径4 N=6381385.838 卯酉圈曲率半径5 R=6363125.179 平均曲率半径6 T? ××× 1531918大地方位角7 P=6352234.773 方位角方向曲率半径8 Z ？54 80 9 B ？22.54722222 10 M=6344807.518 11 N=6381281.239 12 R=6363018.244 13 T ？153.3216667 14 P=6352126.447 15 Z? 80 84 16 B ？22.54722222 17 M=6344804.56 18 N=6381278.236 19 R=6363015.264 20 T ？153.3216667 21 P=6352123.479§1.2边长投影计算使用电磁波测距仪或全站仪观测得到地面某直线的水平距离D 需要经过两次投影才能在高斯平面坐标系中进行计算。
第一次投影是从地面投影到参考椭球面上，第二次投影是从参考椭球面投影到高斯平面上。
(1) 数学模型设投影到参考椭球面上的长度为0S ′，则计算公式为 D R h H R h H D S n g m ng m ⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧+−+−=′220)( (1-7) 式中，m H 为边长D 高程，g h 为边长所在地区大地水准面对于参考椭球面的高度，n R 为边长中点方向上参考椭球面法截弧的曲率半径。
设参考椭球面上的长度0S ′投影到高斯平面上的长度为0S ，则计算公式为 0222200242S R y R y S S m mm ′⎭⎬⎫⎩⎨⎧∆++′= (1-8) 式中，m y 为直线两端点近似y 坐标的平均值，y ∆为直线两端点近似y 坐标的增量，m R 为参考椭球在直线中点的平均曲率半径。
当边长不长时，可以取边长中点方向法截弧曲率半径A R 代替n R 和m R 。
测区内，可以认为g h 为一常数，如果未知，工程计算中可以取g h =0。
(2)、程序与案例分别在1954北京坐标系和1980西安坐标系中计算边长投影到高斯平面上的距离0S 。
1) 变量对照表数学模型变量fx-4800P 变量单位 注释 cCm 极曲率半径 2e ′ E 第二偏心率 B B °′″边长中点的大地纬度 AA°′″边长方向方位角ηU 中间参数 V V 中间参数A R R m 边长中点方向法截弧曲率半径D D m 地面水平距离观测值 m HH m 地面水平距离的高程g hG m 大地水准面对于参考椭球面的高度1y Y m 边长一端的高斯横坐标 2y W m 边长另一端的高斯横坐标 m y K m 边长两端高斯横坐标平均值 y ∆ L m 边长两端高斯横坐标差值 0S ′ J m 投影到椭球面的边长值 0SSm投影到高斯平面的边长值2) 程序程序名：DS Lbl 0 {Z}Z=54C=6399698.90178: E=0.006738525414684: Z=80C=6399596.65199: E=0.006739501819473:Goto 0 {BA}U=EcosB: V=(1+U 2):R= C ÷V(1+(UcosA)2){DHGYW}J=D-((H+G)÷R-(H+G)2÷R 2)D K=(Y+W)÷2: L=Y-WS=J+(K 2÷2÷R 2+L 2÷24÷R 2)J 3) 操作步骤设=B 22°33′04″，=A 107°41′02″，1y =210154m ，2y =213577m ，D =3591.293m ，m H =35.4m ，g h =4.3m ，分别计算其在1954北京坐标系和1980西安坐标系中的高斯平面距离值。
将上述程序以DS 的文件名输入计算器后，按键 5 2及选择程序DS ，按键，屏幕提示及操作步骤如下：步骤 显示 按键 注释1 Z? ××× 54 大地坐标系2 B? ××× 223304 大地纬度3 A? ××× 1074102大地经度4 R=6378004.388 平均曲率半径5 D? ××× 3591.293 地面水平距离观测值6 H ？××× 35.4 地面水平距离的高程7 G ？××× 4.3 大地水准面相对参考椭球面的高度 8 J=3591.270646投影到椭球面的边长值9 Y? ××× 210154 边长一端的高斯横坐标10 W? ××× 213577 边长另一端的高斯横坐标 11 S=3593.252075投影到高斯平面的边长值12 Z? 54 8013 B? 22.5511111114 A? 107.683888915 R=6377899.35416 D? 3591.29317 H ？35.418 G ？4.319 J=3591.27064620 Y? 21015421 W? 21357722 S=3593.25214§1.3高斯投影正算与子午线收敛角的计算将大地经纬度为L 、B 的点投影到高斯平面上并求出其坐标x 、y 称高斯投影正算，令γ为子午线收敛角。
(1) 数学模型⎥⎦⎤⎢⎣⎡+−+++−++=64244222)5861(7201)495(24121m t t m t m Nt X x ηη (1-9)⎥⎦⎤⎢⎣⎡−++−++−+=522242322)5814185(1201)1(61m t t t m t m N y ηηη (1-10)[]523420)2(12)231(60180m t m m t−++++=ηηπγ (1-11)0γ表示以十进制度为单位的子午线收敛角。
式中B t tan =(1-12) 180cos πBlm =(1-13)N 为卯酉圈曲率半径，计算公式见(1-4)，η的计算公式见(1-3)，X 为子午线弧长。
对于克拉索夫斯基椭球：BB B B B B X cos )sin 0039.0 sin 6973.0sin 9238.133sin 7799.32005(8611.1111347530+++−= (1-14)对于IUGG1975椭球：BB B B B B X cos )sin 0039.0 sin 6976.0sin 9602.133sin 8575.32009(0047.1111337530+++−= (1-15)(2)、程序与案例分别在1954北京坐标系和1980西安坐标系中进行高斯投影正算。