fx-4800P缓和曲线和圆曲线坐标正反算程序(正确版)

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计算器测量计算程序抚顺市市政建设(集团)有限公司第一工程处：张和光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注：此程序可循环计算。

曲线计算公式、程序及算例广西建设职业技术学院测量教研室李向民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通用程序清单圆曲线程序：—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δ？--------------------输入交点纵坐标。

线路测设坐标计算程序（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.编制程序过程中已经考虑了零溢出问题。

标准测量曲线坐标计算程序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道路中边桩坐标计算程序程序清单如下：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-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 距中线距离。

一、QXFY 辛甫森公式放样程序1. U“X0”：V“Y0”2. A“XA”：B“Y A”：C“CA”：D“1÷RA”：E“1÷RB”：F“KA”：G“KB”3. Lb1 1：{HLR}：H“KI”：L“JJ”：R“DZ+Y-”4. H>G Goto 15. H<F Goto 16. P=(E-D) ÷Abs(G-F)：Q=Abs(H –F)：Z=P×Q：T=D+Z7. W=C+（Z+2D）Q×90÷π8. M=C+(Z÷4+2D)Q×22.5÷π：N=C+（3Z÷4+2D）Q×67.5÷π：K=C+(Z÷2+2D)Q×45÷π9. X=A+Q（cosC+4(cosM+cosN)+2cosK+cosW）÷12+Rcos(W-L)10. Y=B+Q(sinC+4(sinM+sinN)+2sinK+sinW) ÷12+ Rsin (W-L)11. Pol((X-U)，（Y-V）)J≤0J=J+360J“FWJ”12. I“S”13. Goto 1程序显示说明：须输入参数：X0:置镜点X (对应变量U)Y0:置镜点Y (对应变量V)XA:曲线起算点X (对应变量A)Y A:曲线起算点Y (对应变量B)CA：曲线起算点切线方位角(对应变量C) 1÷RA：1÷半径，即起算点曲率，右偏为正左偏为负(对应变量D)1÷RB：1÷半径，即曲终点曲率，右偏为正左偏为负(对应变量E)KA:起算点里程(对应变量F)KB:曲终点里程(对应变量G)KI:待求点里程(对应变量H)JJ:夹角(与前进方向切线方位角夹角，向右为正)(对应变量L)D“Z+Y－”：偏距，左正右负(对应变量R)计算结果: W: 待求点切线方位角(对应变量W) X:计算点x (对应变量X)Y:计算点y (对应变量Y)FWJ:放样方位角(对应变量J)S:放样距离(对应变量I)K BX 0Y 0本程序依据复化辛甫森公式计算原理改进而成，特点是把曲线按曲率变化点分成若干计算单元单独计算，不论直线（曲率为0）、圆曲线（曲率为1/R ）、卵形曲线，只要按曲率变化点分解弄清变化点曲率半径，右偏左偏曲线由曲率的正负号分别，既可求得该计算单元内任意里程中桩及斜交或正交的边桩。

4800全套计算程序（吐⾎推荐）整套卡西欧计算器CASIO FX-4800 4500测量程序各位同⾏：我长年从事公路桥梁测绘⼯作，收集整理了⼀整套卡西欧计算器CASIOFX-4800P和卡西欧CASIO FX-4500P测量放线程序，该套程序经过反复修正，具有极⾼的可靠性。

现在免费提供给⼤家.程序⽬录：1.1 测边⼤地四边形的条件平差程序1.2 测边中点多边形的条件平差程序1.3 单⼀导线的严密平差程序1.4 单⼀结点导线的严密平差程序1.5 边、⾓后⽅交会的条件平差程序1.6 测边交汇点的条件平差程序1.7 单⼀⽔准路线的平差程序1.8 单结点的平差程序1.9 多结点⽔准⽹的平差程序1.10 多边形⽔准⽹的平差程序1.11 三⾓⾼程差的计算程序1.12 测距边⽔平距离的计算程序1.13 ⾼斯投影座标正、反解和换带的计算程序1.14 平⾯座标转移的计算程序2.1 单边、⾓后⽅交会的计算程序2.2 双边、⾓后⽅交会的计算程序2.3 碎部点座标的计算程序2.4 内、外分点座标的计算程序2.5 房⾓标的计算程序2.6 多次后⽅测⾓交会点2.7 测边交会点的计算程序2.8 导线三维值的计算程序2.9 圆⼼座标和半径的程序2.10 测站点、照准点归⼼改正的计算程序2.11 单边、双站后⽅交会的计算程序3.1 线路中⼼座标的计算程序3.3 红线交点座标的计算程序3.4 直线与圆曲线交点座标的计算程序3.5 垂⾜座标的计算程序3.6 道路不等宽时曲线要素的计算程序3.7 缓和曲线敷设时的计算程序3.8 圆曲线敷设时的计算程序3.9 极座标放线线的计算程序3.10 交点座标的计算程序3.11 宗地有关项⽬的计算程序3.12 丘地有关项⽬的计算程序4.1 导线X、Y、Z值的近似平差程序4.2 ⽆定向导线的计算程序4.3 测⾓前⽅交会点座标的计算程序4.4 测边交会座标的计算程序4.5 测⾓后⽅交会点座标的计算程序4.6 双边、⾓后⽅交会点座标的计算程序4.7 平⾯座标转移的计算程序4.8 测站点、照准点归⼼改正的计算程序4.9 三⾓⾼程⾼差的计算程序4.10 测距边⽔平距离的计算程序4.11 单⼀导线的严密平差程序4.12 边、⾓后⽅交会点的严密平差程序4.13 测边交会点的严密平差程序5.1 图跟导线点座标的计算程序5.2 碎部点座标的计算程序5.3 内、外分点座标的计算程序5.4 ⾓座标的计算程序5.5 圆⼼座标的和半径的计算程序5.6 交点座标的计算程序5.7 单边、⾓交会点座标的计算程序5.8 宗地有关项⽬的计算程序5.9 单边、双站后⽅交会座标的计算程序6.1 线路边桩座标的计算程序6.3 直线与曲线交点座标的计算程序6.4 道路不等宽时曲线元素的计算程序6.5 垂⾜座标的计算程序6.6 圆曲线敷设的计算程序6.7 缓各曲线敷设的计算程序6.8 极座标放线的计算程序6.9 单⼀⽔准路线的近似平差计算程序6.10 丘地有关项⽬的计算程序公路中、边桩坐标计算及放样程序功能：任意半径曲线中桩、任意交⾓边桩坐标计算放样、中桩切线⽅位⾓、任意点置仪放样、两点间距离计算、两点⽅位⾓计算。

一、源程序(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次方，输入程序时切莫输错。

线路中边桩任意里程坐标正、反算改编程序（CASIO fx-5800P）1. 正算子程序(SUB1)A=0.26：B=0.74：K=0.02：L=0.3300094782：F=1-L：M=1-K：X=U+W(Acos(G+57.2958QKW(1/P+KWD))+B cos(G+57.2958QLW(1/P+LWD))+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)+57.2958QW(1/P+WD)：Z=(J-Y)cosL-(I-X)sinL：AbsZ<1 E-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＜75685.656 ＝＞O=75172.813：U=10257.769：V=615287.070：G=325-46-1.5：H =512.843：P=2600：R=2600：Q=1：Goto0Δ←┘(JD65)S＜75884.442＝＞O=75685.656：U=10707.383：V=615042.113：G= 337-04-06.6：H=198.786：P=2600：R=1E45：Q=0：Goto0Δ←┘(直线)S＜76044.442＝＞O=75884.442：U=10890.459：V=614964.660：G=337-04-06.6：H=160：P=1E45：R=1050：Q=-1：Goto0Δ←┘(JD66第一缓和段)S＜76392.111＝＞O=76044.442：U=11036.146：V=614898.615：G= 332-42-11.19：H=347.669：P=1050：R=1050：Q=-1：Goto0Δ←┘(J D66圆曲线段)S＜76552.111＝＞O=76392.111：U=11313.328：V=614691.387：G= 313-43-54.11：H=160：P=1050：R=1E45：Q=-1：Goto0Δ←┘(JD66第二缓和段)S＜78271.813＝＞O=76552.111：U=11417.893：V=614570.338：G= 309-21-58.7：H=1719.702：P=1E45：R=1E45：Q=0：Goto0Δ←┘(直线段)S＜78431.813＝＞O=78271.813：U=12508.659：V=613240.824：G= 309-21-58.7：H=160：P=1E45：R=1050：Q=-1：Goto0Δ←┘(JD67第一缓和段)S＜78810.831＝＞O=78431.813：U=12606.944：V=613114.623：G= 305-00-3.29：H=379.018：P=1050：R=1050：Q=-1：Goto0Δ←┘(JD 67圆曲线段)S＜78970.831＝＞O=78810.831：U=12764.225：V=612772.038：G= 284-19-7.91：H=160：P=1050：R=1E45：Q=-1：Goto0Δ←┘(JD67第二缓和段)S＜79284.977＝＞O=78970.831：U=12795.865：V=612615.240：G= 279-57-12.5：H=314.146：P=1E45：R=1E45：Q=0：Goto0Δ←┘(直线段)S＜79444.977＝＞O=79284.977：U=12850.165：V=612305.822：G= 279-57-12.5：H=160：P=1E45：R=750：Q= 1：Goto0Δ←┘(JD68第一缓和段)S＜79607.747＝＞O=79444.977：U=12883.388：V=612149.392：G= 286-03-54.08：H=162.77：P=750：R=750：Q= 1：Goto0Δ←┘(JD6 8圆曲线段)S＜79767.747＝＞O=79607.747：U=12944.985：V=611999.072：G= 298-29-59.22：H=160：P=750：R=1E45：Q= 1：Goto0Δ←┘(JD68第二缓和段)S＜79927.747＝＞O=79767.747：U=13031.084：V=611864.310：G= 304-36-40.8：H=160：P=1E45：R=750：Q=-1：Goto0Δ←┘(JD69第一缓和段)S＜80152.430＝＞O=79927.747：U=13117.183：V=611729.547：G= 298-29-59.22：H=224.683：P=750：R=750：Q=-1：Goto0Δ←┘(JD6 9圆曲线段)S＜80312.430＝＞O=80152.430：U=13193.440：V=611519.093：G= 281-20-6.98：H=160：P=750：R=1E45：Q=-1：Goto0Δ←┘(JD69第二缓和段)S＜80472.430＝＞O=80312.430：U=13213.651：V=611360.456：G= 275-13-25.4：H=160：P=1E45：R=920：Q= 1：Goto0Δ←┘(JD70A 第一缓和段)S＜80764.750＝＞O=80472.430：U=13232.840：V=611201.665：G= 280-12-21.47：H=292.320：P=920：R=920：Q= 1：Goto0Δ←┘(JD7 0A圆曲线段)S＜80924.750＝＞O=80764.750：U=13329.149：V=610926.967：G= 301-56-13.33：H=160：P=920：R=1E45：Q= 1：Goto0Δ←┘(JD70 A第二缓和段)S＜＝＞O=75018.789：U=10119.920：V=615344.580：G=325-46-1. 5：H=591.742：P=3000：R=3000：Q= 1：Goto0Δ←┘(Z JD2圆曲线段)S＜75868.894＝＞O=75610.531：U=10638.706：V=615061.939：G= 337-04-6.6：H=258.363：P=3000：R=1E45：Q=0：Goto0Δ←┘(直线段)S＜76028.894＝＞O=75868.894：U=10876.651：V=614961.272：G= 337-04-6.6：H=160：P=1E45：R=900：Q=-1：Goto0Δ←┘(Z JD3第一缓和段)S＜76304.04＝＞O=76028.894：U=11022.044：V=614894.617：G= 331-58-31.95：H=275.146：P=900：R=900：Q=-1：Goto0Δ←┘(Z J D3圆曲线段)S＜76464.04＝＞O=76304.04：U=11241.555：V=614730.506：G= 314-27-33.35：H=160：P=900：R=1E45：Q=-1：Goto0Δ←┘(Z JD3第二缓和段)S＜78325.868＝＞O=76464.04：U=11346.622：V=614609.911：G= 309-21-58.7：H=1861.828：P=1E45：R=1E45：Q=0：Goto0Δ←┘(直线段)S＜78485.868＝＞O=78325.868：U=12527.535：V=613170.519：G= 309-21-58.7：H=160：P=1E45：R=1000：Q=-1：Goto0Δ←┘(Z JD4第一缓和段)S＜78839.218＝＞O=78485.868：U=12625.657：V=613044.196：G= 304-46-57.52：H=353.350：P=1000：R=1000：Q=-1：Goto0Δ←┘(Z JD4圆曲线段)S＜78999.218＝＞O=78839.218：U=12772.319：V=612724.740：G= 284-32-13.68：H=160：P=1000：R=1E45：Q=-1：Goto0Δ←┘(Z JD4第二缓和段)S＜80472.391＝＞O=80312.391：U=13213.651：V=611360.455：G= 275-13-25.4：H=160：P=1E45：R=936.13：Q= 1：Goto0Δ←┘(Z JD 1A第一缓和段)S＜80707.393＝＞O=80472.391：U=13232.720：V=611201.648：G= 280-07-12.42：H=235.002：P=936.13：R=936.13：Q= 1：Goto0Δ←┘(Z JD1A圆曲线段)S＜80867.397＝＞O=80707.393：U=13302.398：V=610977.859：G= 297-27-38.28：H=160：P=936.13：R=1E45：Q= 1：Goto0Δ←┘(Z J D1A第二缓和段)LB1 0：D=(P-R)÷(2HPR)：N=1=>Goto 1：≠>Goto 2Δ←┘Lbl 1：{TZ}：Z：W=Abs(S-O)：Prog "SUB1"：X"XS"=X◢Y"YS"=Y◢F"FS"=F-90◢C“XC”：E“YC”：I=0：J=0：Pol（X-C，Y-E）：I“I=” ◢J∠0＝＞J“J=”+360◢≠＞J“J=” ◢ΔGoto4←┘Lbl 2：{XY}：XY：I=X：J=Y：Prog "SUB2"：S"S"=O+W◢Z"Z"=Z◢G o t o4。

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 FX-4800P施工放样程序二零零六年四月六日坐标正算程序1 X“X0”:Y“Y0”:S:T2 N=X+ScosT:E=Y+SsinT3 N=N▲E=E▲4 输入数据说明:X0─已知点X坐标；Y0─已知点Y坐标；S─已知点到待定点之间的距离；T─已知点到待定点之间的方位角5 计算结果说明：N─待定点X坐标；E─待定点Y坐标6 检核数据已知：X O=500.00m；Y O=500.00m；S OA=150.00m；T OA=40°25′10″结果：N A=614.1977m；E A=597.2567m7 简图：O(原点)坐标反算程序1Lb1 0:{NE}:X“X0”:Y“Y0”:pol((N-X),(E-Y):J≤0 => J “＜T”=J+360▲≠=>J“＜T”=J▲△I“＜S”=I▲G o t o 02输入数据说明: （“=>”为条件转移符号、“≠=>”为无条件转移符号）X0─O点X坐标；Y0─O点Y坐标；N─A点X坐标；E─A点Y坐标3计算结果说明：S─O点到A点之间的距离；T─O点到A点之间的方位角4检核数据已知：X O=500.00m；Y O=500.00m； N A=614.1977m；E A=597.2567m结果： S OA=150.00m；T OA=40°25′10″5 简图：(已知点) O(已知点)大地坐标转施工坐标程序1 Lb1 1:A“X0”:B“Y0”:T:{NE}:X=（N-A）×cosT+（E-B）×s inT:Y=（E-B）×cosT-（N-A）×s inT2 X=X▲Y=Y▲G o t o 13 输入数据说明:X0、Y0─施工坐标原点在大地中的坐标；T─施工坐标系相对于大地坐标系的夹角；N─大地坐标系中的X坐标；E─大地坐标系中的Y坐标4计算结果说明：X─施工坐标系中的X坐标；Y─施工坐标系中的Y坐标5 检核数据已知： X O=316.689m；Y O=501.123m；T=-8°13′36″N=755.500m；E=740.800m结果： X=400.000m；Y=300.000m6 简图及数据：详见中等专业学校教材《工程测量学》第一章第3～6页图1-1坐标系统换算示意图；表1-1 测图坐标换算施工坐标已知数据表；表1-2 坐标换算与校核计算表。

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的里程桩。

[正确] 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、◢为输出指令，若在后面加上，即可显示前面的计算结果输出在屏幕上。

福建林业职业技术学院工程系土木工程室高国兴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）←。

一、本人编写整段曲线坐标程序一则。

可作为参考，从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个程序。