fx-4800P缓和曲线和圆曲线坐标正反算程序(正确版)
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。
卡西欧计算器程序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注：此程序可循环计算。
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。
4800P计算器圆曲线、缓和曲线通用程序清单
4800P通用程序清单圆曲线程序：—YQX（拼音）O“Xδ”：H“Yδ”：G“а0”：M“μв”：P“YP,ZP＝+1,-1”：R“R”： L“L C”： C“L h”: Q“ZYD.ZH”：LbI 1：{KD}：S＝K－Q：S＜0=＞S＝0:⊿ S＞L=＞S＝L:⊿V＝S÷R×180÷π：N＝(M－G)÷2：N＜0=＞N＝-N:⊿N＞90=＞N＝180-N:⊿U＝C÷2R×180÷π：C＝0=＞E＝R÷CosN:≠＞E＝(R＋C2÷24R)÷CosN:⊿ P＝1=＞F＝M＋(180－2N)÷2：T＝G＋270：W＝T＋U＋V：Z＝R－D:≠＞F＝M－(180－2N)÷2：T＝G＋90：W＝T－U－V:Z＝R＋D:⊿ Fix 3：X“X”＝O＋E×CosF＋Z×CosW◢Y“Y”＝H＋E×SinF＋Z×SinW◢Goto 1缓和曲线程序：—HHQX（拼音）O“Xδ”：H“Yδ”： G“а0”: Z“Bp g”：P“P=1,-1,2,-2”： R“R”： C“L h”： Q“QSD.ZH”：LbI 0：{KD}：S＝K－Q：J＝Abs P：J＝1=＞Goto 1:≠＞J＝2=＞Goto 2:≠＞P“P=1，-1，2，-2；Ran#”◢⊿ Goto 0：LbI 1：S＜0=＞S＝0:⊿ S＞C=＞S＝C:⊿ Goto 3：LbI 2：S≥0=＞S＝0:≠＞S＝Abs S:⊿ S＞C=＞S＝C:⊿ Goto 3：LbI 3：V＝S2÷(2×R×C)×180÷π： E＝S∧3÷(6×R×C)：F＝S－S∧5÷(40×R2×C2)：Goto 4：LbI 4：P＝-1=＞Goto 5:≠＞P＝1=＞Goto 6:≠＞P＝-2=＞Goto 7:≠＞P＝2=＞Goto 8:⊿LbI 5：U＝G＋270：W＝G－V＋90：Goto 9：LbI 6：U＝G＋90：W＝G＋V＋90：F＝-F：Goto 9：LbI 7：U＝G＋270：W＝G－V＋270：Goto 9：LbI 8：U＝G＋90：W＝G＋V＋270：F＝-F：Goto 9：LbI 9：L＝C÷2－C∧3÷(240×R2)：Z=AbsZ：T＝(R＋C2÷(24×R))×Tan(Z÷2)：M＝O＋(L＋T)×Cos(G＋180):N＝H＋(L＋T)×Sin(G＋180): Fix 3：X“X”＝M＋E×CosU－F×SinU＋D×CosW◢Y“Y”＝N＋E×SinU＋F×CosU＋D×SinW◢Goto 0：程序注译圆曲线程序---YQX（拼音）Xδ？--------------------输入交点纵坐标。
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计算器程序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点桩号。
卡西欧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）。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
缓和曲线和圆曲线坐标正反算程序主程序“TYQXJS”
Lb1 0↙→（EXE）
{NUVOGHPRQ}:“1.SZ=>XY”: “2.xy=>SZ”:
N:U“QDX”:V“QDY”:O“QDLC”:G“FWJ”:
H“LS”:P“RO”:R“RN”:Q“ZP=-1,YP=+1,ZZ=0”
:C=1÷P:D=(P-R)÷(2HPR):E=180÷π:N=1=>Goto 1:≠=>Goto 2◣↙Lb1 1:{SZ}:SZ:W= Abs(S-O): Prog“1”:
X“XS”=X◢
Y“YS”=Y◢
F“FS”=F-90◢
Goto 3↙
Lb1 2:{XY}:XY:I=X:J=Y: Prog“2”:S“S”=O+W◢
Z“Z”=Z◢
Goto 3↙
Lb1 3↙
{DE}:E“QX-JJ,Z-1,Y+1”:D“BZ-JL”↙
F=F+E↙
X=X+D Cos F◢
Y=Y+D Sin F◢
{DE}:D“BZ-JJ”:E“JJ,Z-1,Y+1”↙
F=F+E↙
X=X+D Cos F◢
Y=Y+D Sin F◢
Goto 0
子程序1：“1”
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+Q EFW(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+QE W(C+WD)+90:X=X+Z Cos F:Y=Y+Z Sin F
子程序2: “2”
T=G-90:W=Abs((Y-V)Cos T-(X-U)Sin T):Z=0:
Lb1 0:Prog“1”:L=T+QEW(C+WD):Z=(J-Y)Cos L-(I-X)Sin L:Abs
Z<1E-6=>Goto 1: ≠=>W=W+Z:Goto 0◣↙(E为：4800P键盘的EXE键) Lb1 1:Z=0Prog“1”:Z=(J-Y)÷Sin F◣
注：第一缓和曲线起点半径输入无穷大（10 45），终点输入圆曲线半径；第二缓和曲线起点半径输入圆曲线半径，终点半径输入无穷大（10 45）；圆曲线输入给出的起点和终点半径；直线段则都输入无穷大（10 45）。
其中起点切线方位角用此程序计算。
◣(代替空心)
坐标正算
程序“ZBZS”
Lb1 0↙
X“HSX”:Y“HSY”:U“CZX”:V“CZY”↙
X-U≥0=>Goto 1: ≠=>Goto 2◣↙
Y-V≥0=>F“FWJ”=tan-1((Y-V)÷(X-U)): ≠=>F“FWJ”=360+tan-1((Y-V)÷(X-U))◢
F“FWJ”◢
◣↙
Goto 3◣↙
Lb1 2↙
F“FWJ”=180+ tan-1((Y-V)÷(X-U))◢
Goto 3◣↙
Lb1 3↙
L“QSJL”:O“JJ”↙
N“QSX”=U+L Cos（F+O）◢
E“QSY”=V+L Sin（F+O）◢
Goto 0↙
↙
坐标反算
程序“ZBFS”
Lb1 0↙
U“CZX”:V“CZY”:X“HSX”: Y“HSY”↙
X-U≥0=>A“CZFWJ”=tan-1（（Y-V）÷(X-U)）◢
≠=>A“CZFWJ”=tan-1（（Y-V）÷(X-U)）+180◢
Goto 1◣↙
N“QSX”：E“QSY”↙
N-U≥0=>B“QSFWJ”=tan-1（（E-V）÷(N-U)）◢
≠=>B“QSFWJ”=tan-1（（E-V）÷(N-U)）+180◢
Goto 2◣↙
Lb1 2↙
O“JJ”=B“F2”-A“F1”◢
L“HSJL”=√（（Y-V）2+（X-U）2◢
L“QSJL”=√（（E-V）2+（N-U）2◢
Goto 1↙
竖曲线
程序“SQX”
A“JDK”：B“HD”：C“R”：M“I1”:D=M÷100:N“I2”：E=N÷100:H=C ÷2×Abs（D-E）：D﹥E=>F=-1：≠=> F=1：◣↙
Lb1 1：{G}：G“ZHUANG”:G﹤A-H=>I=-D：Goto 2：≠=>G≤A =>I=-D：Goto 3：≠=>G<A+H=>I=E: Goto 3：≠=> I=E: Goto 2：◣◣◣↙Lb1 2：J“HD”=B+I×Abs（G-A）◢
Goto 1：↙
Lb1 3：J“HD”=B+ I×Abs（G-A）+ F×(H-Abs(G-A))2÷(2C)◢Goto 1：↙。