4800程序指导

4800计算器程序（缓和曲线任意中、边桩）LBI 0：X：Y：{I}：{J}：A：R：U：N：{D}：{O}：Z：W：W=1＝＞GOTO 1：≠＞GOTO 2←LBI 1：B=90L2÷（ΠRA）：←C=√（（（I-J）-（I-J）5÷（40R2A2））2+（（I-J）3÷（6RA））2）：Q=U+BN÷3：S=U+BN+Z：T=U+BN+Z+180：GOTO 3：←LBI 2：E=90（I-J）÷（ΠR）：C=2R SinE:B=90A÷（ΠR）:Q=U+(B+E)N:S=U+(B+2E)+Z:T=U+(B+2E)N+Z-180:GOTO 3: ←LBI 3：F“XA”=X+C CosQ▲G“YA”=Y+C SinQ▲H“XR”=F+D CosS▲H“YR”=G+D SinS▲H“XL”=F+O CosT▲H“YL”=G+O SinT▲GOTO 0←说明：X Y—-起算点坐标（计算缓和曲线部分输ZH点、计算圆曲线部分输HY点）；I J---为计算桩号，I为大里程，J为小里程，保证（I-J）始终为正；A—-为缓和曲线长；R—-半径；U—-起始方位角（直缓点）；D O—计算边中距；Z—-一般为90度；N—-判别曲线左右偏位，右转输+1、左转输-1；W---当W=1时，程序计算缓和曲线部分，当W=0时，程序计算原曲线部分；F“XA”、 G“YA”---计算桩号的中桩坐标；H“XR”、 H“YR”---计算前进方向（F“XA”、 G“YA”对应）的右边桩坐标，H“XL”、 H“YL”---计算前进方向（F“XA”、G“YA”对应）的左边桩坐标注意：在计算YH至HZ间坐标时，X Y应输入HZ点的坐标值，U为HZ点的方位角-180度，N也相应相反，这时计算的左右边桩与线路前进方向的左右边桩正好相反，计算时要注意取值。

4800程序使用说明书术语：连续曲线数据块、独立曲线数据块、长短链功能：1、放样对应里程中心点2、放样对应里程旋转任意角度外移点（图A点）3、放样对应里程旋转任意角度外移点的外移点（图B点）4、获取放样点坐标及对应里程中心点沿线路前进方向切线方位角程序输入注意事项：1、常用程序文件放在前面，便于调用，一般先建立“CIRC LE”、“XY-AL”、“AL-XY”三文件，文件名千万不能输错。

2、加粗带下划线词组为计算器中的函数或符号，只能从计算器中调出此函数或符号，不能从键盘输入。

符号←表示回车键“EX E”坐标输入通用格式：屏幕显示OPP NAME?（置镜点名称）或CPP NAME?（后视点名称）：1：输入-1：手动输入坐标X,Y2：输入大于1数字：输入用数字所代表的控制点名。

如数据库中没有本数字所代表的点名，则系统显示Syn ERROR in…（文件名）。

曲线放样程序（CIRCLE）运行说明：步骤1：屏幕提示“LOAD?”（输入曲线要素）1.1：输入 1（默认）：输入里程后将自动判断所在曲线并调用（此功能必须保证曲线数据库中有连续曲线数据块，具体详见数据库建立）。

1.2：输入-1：屏幕提示手动输入各曲线要素，ZJ（转角）、R（半径）、L0(缓和曲线长)、ZH KM(直缓里程)、JDX(曲线交点X坐标)、JDY(曲线交点Y坐标)、A0(起始直线边前进方向方位角)、1 R,-1 L(右偏输入1，左偏输入-1)。

1.3：输入大于1的数字：输入用数字代表的曲线名称，比如13，表示放样点位于曲线 13 上（曲线名称在建立数据库是自己命名，但一定要大于1）。

如库中没有此名称，系统显示“Syn ERROR in CI RCLE”。

（此功能必须保证曲线数据库中有独立曲线数据块）步骤2：屏幕提示“OPP NAME?”（输入置镜点名称）详见坐标输入格式。

步骤3：屏幕提示“DK+M?”（输入放样点对应中心里程）如K15+002.35应输入15002.35。

卡西欧4800\4850程序一：坐标反算：(ZBFS)①程序步骤：C”X1”:D”Y1”:E”X2”:F”Y2”:Fixm:Pol(E-C,F-D:I”S1-2=”▲J≤Ｏ＝＞J=J+360△J”A1-2=”②操作过程：ZBFS→EXE→输入X1值→EXE→输入Y1值→EXE→输入X2值→EXE→输入Y2值→EXE→EXE→EXE注：①：S1-2:计算得出的距离。

②：A1-2:计算得出的角度。

（按shift°′″转换为60进制的角度）③：此程序可循环计算。

③案例：已知：D1坐标（7811.23，606.136），D2坐标（7805.915，431.910）求解： a D1-D2(A1-2)=268°15′09.56 ″S D1-D2(S1-2)=174.3071二：坐标正算：(ZBZS)①程序步骤：C”X1”:D”Y1”:S”S1-2”:R”A1-2”:Fixm:X”X2”=C+Rec(S,R) ▲Y”Y2”=D+J②操作过程：ZBZS→EXE→输入X1值→EXE→输入Y1值→EXE→输入S1-2距离值→EXE→输入A1-2角度（例268°15′09.56 ″）值→EXE→EXE注：①：X2计算得出坐标，Y2计算得出坐标。

②：此程序可循环计算。

③案例：已知：D1坐标（7811.23，606.136），两点距离：S D1-D2(S1-2)=174.3071方位角：a D1-D2(A1-2)=268°15′09.56 ″求解：D2坐标（7805.915，431.910三：线路曲线计算程序:（XLQXJS）(万能公式)①程序步骤：1.Lbl 0:{E G}：A“XA”：B“YA”：C“CA”：D“1/RA”：E“1/RB”：F“DKA”：G“DKA”2.Lbl 1:{H O R}:H“DKI”：O“DL”：R“DR”：H＞G=﹥Goto 3⊿3.P=(E-D)/Abs(G-F):Q=Abs(H-F):I=P×Q:T=D+I4.J=C+(I+2D)×Q×90/π ▲5.M=C+(I/4+2D)Q×45/(2π):N=C+(3I/4+2D)Q×135/(2π):6.K=C+(I/2+2D)Q×45/π7.X=A+Q(cosC+4(cosM+cosN)+2cosK+cosJ)/12 ▲8.Y=B+Q(sinC+4(sinM+sinN)+2sinK+sinJ)/12 ▲9.U“XL”=X+Ocos(J-90) ▲10.V“YL”=Y+Osin(J-90) ▲11.W“XR”=X+Rcos(J+90) ▲12.Z“YR”=Y+ Rsin(J+90) ▲13. A=X：B=Y：C=J：D=T：F=H：Goto 1Lbl 3: A=X：B=Y：D=E：F=G：C=J：Goto 0程序说明：A- 曲线元起点A的坐标；B- 曲线元起点B的坐标；C- 曲线元起点A的切线坐标方位角；F- 曲线元起点A的里程；G- 曲线元起点B的里程；H- 曲线上待求点i的里程；D- 曲线元起点A的曲率；E- 曲线元终点B的曲率；XL-左边线点位X坐标；YL-左边线点位Y坐标；XR-右边线点位X坐标；YR-右边线点位Y坐标；X- 中线点位纵坐标；Y- 中线点位横坐标；DL-左边线距中线平距；DR-右边线距中线平距；该程序需要输入的数据为：（1）.曲线元起点A的坐标及切线坐标方位角，计算器上用“XA”，“YA”，“CA”显示；（2）.曲线元起点A和B的曲率，计算器上用I÷RA，I÷RB显示（曲线左偏时取“-”）；（3）.曲线元起点A和终点B的里程，计算器上用“DKA”，“DKB”显示；（4）.输入待求点里程和该点距左右的水平距离，计算器上用“DKI”，“DL”，“DR”显示；每算完一个待求点的中线及边线坐标，程序又让输入下一点的“DKI”，“DL”，“DR”当输入的“DKI”大于“DKB”时程序中显示“1÷RB”和“DKB”，此时输入下一个曲线元起点的曲率和里程，然后重复步骤（4），即可计算下一个曲线中线及边线点位坐标。

竖曲线高程的计算公式1.在直线段上的计算公式：H=G±(K-B)*I2.在曲线段上的计算公式：H=G±(K-B)*I+N(K-Z)2/2R式中：K—待求点的桩号Z—直圆点桩号R—曲线半径H—待求点的高程G—变坡点的高程B—变坡点的桩号I —纵坡度N—常数（凸为-1，凹为+1）4500竖曲线计算实例Lbl 0{K}Lbl 1K≤A→H=G1+(K-变坡点1桩号)*I▲≠→Goto 2△Goto 最后一行Lbl 2K≤B→H=G1+(K-变坡点1桩号)*I-(K-A)2/2R ▲≠→Goto 3△Goto 最后一行Lbl 3K≤C→H=GJD1+(K-JD1桩号)*I▲≠→Goto 4△Goto 最后一行Lbl 4K≤D→H=GJD1+(K-JD1桩号)*I+(K-C)2/2R ▲≠→Goto 5△Goto 最后一行以此类推············Lbl 最后一行U=Y-(H-0.7+1.25*0.18+1.1*0.02) ▲V=U+0.02*14.03 ▲Goto04800竖曲线计算实例Lbl 0{k}Lbl 1K≤A→H=G1+(变坡点1桩号-K)*I▲≠→Goto 2△Goto 最后一行Lbl 2K≤B→H=G2+(变坡点2桩号-K)*I-(K-A)2/2R ▲≠→Goto 3△Goto 最后一行Lbl 3K≤C→H=GJD1+(K-JD1桩号)*I▲≠→Goto4△Goto最后一行Lbl 4K≤D→H=GJD1+(K-JD1桩号)*I+(K-C)2/2R ▲≠→Goto5△Goto最后一行以此类推············Lbl 最后一行U=Y-(H-0.7+1.25*0.18+1.1*0.02) ▲V=U+0.02*14.03 ▲Goto 05800竖曲线计算实例“SQXJSCX”↙Lb1 1↙CLs:Fix 3↙“K=”?k◢(计算点里程输入)If k<67549.755 AND K≥66894.3 :Then -0.00052→A : 0→B : 67394.3→S : 67.37→G : 600000→R : Goto 2 : IfEnd↙If k<68708.391 AND K≥67549.755 :Then -0.0048→A : -0.00052→B : 68494.3→S :66.8→G : 100000→R : Goto 2 : IfEnd↙If k<本段竖曲线终点里程AND K≥前一竖曲线终点里程:Then -0.0048→前坡（大里程向）A : -0.00052→后坡（小里程向）B : 68494.3→竖曲线交点里程S :66.8→交点高程G : 100000→曲率半径R : Goto 2 : IfEnd↙………依次类推，计算原始数据完成输入，坡度换算成小数。

4800计算器程序编辑步骤2一、文件名QXFY(曲线放样主程序)W=1:{WQ}:W“LOAD”:W=-1=>A“ZJ”:B“R”:C“L0”:D“ZH KM”:E“JDX”:F“JDY”:G“A0”:H“1 R,-1 L”:⊿←{Q}:Q“OPP NAME”:Q=-1=>U“X0”:V“Y0”:≠>Prog“CON DAT”:U=4.444=>::⊿⊿←Lb1 1:{Z}:Z“DK+M”:Int (Z/1000)>0=>K=Int (Z/1000): ≠>Z=Z+1000K⊿←W≠-1=>Prog “CIR DAT”:A=9=>::⊿⊿←H≠-1=>H≠1=>::⊿⊿←M=Z-D:Prog “CENTER”←{L}:L“OUT L”:L>0=>{R}:R“OUT ANGLE”:X=X+Lcos(P+R):Y=Y+Lsin(P+R) ⊿L<0=>L=AbsL:{R}:R“OUT ANGLE1”:O=P+R:X=X+LcosO:Y=Y+LsinO:{LR}:L“OUT L2”:R“OUT ANGLE2”:L=AbsL:O=O+R +180:X=X+LcosO:Y=Y+LsinO⊿←PoI(X-U,Y-V)←“L=”:I◢( I“L=”◢)J<0=>J=J+360:⊿←“A=”:J→DMS◢(J“A=”◢)O=-1:{O}:O“VIEWXY”:O≠-1=>”X=”:X◢”Y=”:Y◢”P”:P→DMS◢⊿←(X◢Y◢P◢⊿←)Goto 1←二、文件名XY-AL(坐标反算){Q}:Q“OPP NAME” ←Q≠-1=>Prog“CON DAT”:U=4.444=>::⊿A=U:B=V: ≠>A“X0”:B“Y0”: ⊿I=0:J=0←Lb1 1:{Q}:Q“CPP NAME”:Q=-1=>{DC}:C“Xc”:D“Yc”: ≠>Prog“CON DAT”:U=4.444=>::⊿C=U:D=V: ⊿←Po1(C-A,D-B) ←“L=”:I◢(I “L=” ◢)J<0=>J=J+360⊿←“A=”:J→DMS◢(J“A=” ◢)Goto 1←三、文件名CAL XY(导线计算)O=2←Lb1 1:{O}:O“1:X0,Y0-A06个空格2:X0,Y0-XC,YC”:O≠1=>O ≠2=>Goto 1: ⊿⊿←I=0:J=0←{Q}:Q“OPP NAME”:O≠-1 Prog“CON DAT”: U=4.444=>::⊿A=U:B=V: ≠>A“X0”:B“Y0”: ⊿←O=1=>E“A0”: ≠>{Q}:Q“CPP NAME”:Q=-1=> C“Xc”:D“Yc”: ≠>Prog“CON DAT”:U=4.444=>::⊿C=U:D=V: ⊿Po1(C-A,D-B):E=J: ⊿←F=E:X=A:Y=B←Lb1 2:{ML}:M“LEFT ANGLE”:L“DISTANCE”:F=F+M:Rec(L ,F):X=X+I◢Y=Y+J◢F=F-180:Goto 2四、文件名CENTER（曲线放样子程序）S=ABπ/180+C：O=C2/24B：T=（B+O）tan（A/2）+C/2-C^3/240B2：N=-1：M>S-C=>M=S-M：N=1：⊿←M≤0=>X=M:Y=0:P=0: Goto 1: ⊿←M≤C=>X=M-M^5/40B2C2:Y=M^3/6BC-M^7/(336*B^3*C^3):P= 90M2/πBC：Goto 1: ⊿←M≤S-C=>P=180(M-C)/πB+90C/πB:X=BsinP+C/2-C^3/240 B2:Y=B(1-cosP)+O: ⊿Lb1 1:N=1=>P=A-P:I=X:J=Y:X=T+(T-I)cosA-JsinA:Y=(T-I)sinA+JcosA: ⊿←I=X:J=HY:X=E+Tcos(180+G)+IcosG-JsinG:Y=F+Tsin(180+G )+IsinG+JcosG:P=G+HP←五、文件名CIR DAT(曲线要素数据库)A=9：W≠1=> Goto 1⊿←Z≤曲线交点代号=>A=转向角：B=半径：C=缓和曲线长：D=直缓点里程：E=交点X坐标：F=交点Y坐标：G=切线方位角：H=右偏+1，左偏-1：Goto 1⊿←W=编辑顺序同上。

卡西欧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全套计算程序（吐⾎推荐）整套卡西欧计算器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 丘地有关项⽬的计算程序公路中、边桩坐标计算及放样程序功能：任意半径曲线中桩、任意交⾓边桩坐标计算放样、中桩切线⽅位⾓、任意点置仪放样、两点间距离计算、两点⽅位⾓计算。

直线坐标计算：1、N″X0″:K″Y0″:Q:Z″HZ″2、Lb1 0:﹛S﹜:S″JSD″:L=S-Z3、C″X1″=L*COS Q+N◢4、D″Y1″=L*Sin Q+K◢5、﹛BP﹜:A=Q+P6、E″X2″=C+B*COS A◢7、G″Y2″=D+B*Sin A◢8、Goto 0注：X0,Y0表示起算点的纵横坐标；Q表示起算方位角；HZ表示起算点里程；JSD 表示计算点里程；X1，Y1表示计算点坐标；P表示边桩夹角；B表示边桩到计算点的距离；X2，Y2表示边桩坐标。

线路综合计算（由曲线起点至第二条曲线起点）：1、R″R″：I″I″：C″ZXJ″：A″ZX″：B″ZY″：W″ZA″：O″ZH″2、P=I2/24/R-I4/2688/R33、M=I/2-I3/240/R24、T=（R+P）*tan（C/2）+M5、F=I/2/R*180/π6、L=（C-2F）*π/180*R+2I7、G″L0：L=1；R=2″8、G=1=＞N=-1◣9、G=2=＞N=1◣10、Lb1 1：{D}：D″JSD″：Z=D-O11、Pause 012、Z＞L=＞Q=W+N*C◢13、X=A+T*COS W+（T+Z-L）*COS Q◢14、Y=B+T*Sin W+（T+Z-L）*Sin Q◢15、Goto 3:◣16、Z＜I=＞Q=W+N*Z2/2/R/I*180/π:X=Z-Z5/40/R2/I2+Z9/3456/R4/I4:Y=Z3/6/R/I-Z7/336/R3/I3+Z11/42240/R5/I5:Goto 2:◣17、Z＜L-I=＞S=(Z-I)/R*180/π+F:Q=W+N*((Z-I)/R*180/π+F):X=R*Sin S+M：Y=R*（1-COS S）+P：Goto 2：◣18、U=L-Z-（L-Z）5/40/R2/I2+（L-Z）9/3456/R4/I419、V=（L-Z）3/6/R/I-(L-Z)7/336/R3/I3+(L-Z)11/42240/R5/I520、Q=W+N*C-N*(L-Z)2/2/R/I/π*18021、X=T*(1+COS C)-U*COS C-V*Sin C22、Y=T*Sin C-U*Sin C+V*COS C23、Lb1 2:Q＜0=＞Q=Q+360◢24、◣25、Q≥360=＞Q=Q-360◢26、◣27、Q=Q◢28、S=√(X2+Y2)29、K=tan-1(Y/X)30、X=A+S*COS(W+N*K)◢31、Y=B+S*Sin(W+N*K)◢32、Goto 3:◣33、Lb1 334、{JE}35、X″EX″=X+E*COS(Q+J)◢36、Y″EY″=Y+E*Sin(Q+J)◢37、Goto 1注：Ｒ表示曲线半径；Ｉ表示缓和曲线总长；ＺＸＪ表示转向角；ＺＸ、ＺＹ表示直缓点（直圆点）纵横坐标；ＺＡ表示起算方位角；ＺＨ表示直缓（直圆）点里程；Ｌ（１）或Ｒ（２）表示左转或右转；ＪＳＤ表示计算点里程；Ｑ表示计算点切线方位角；Ｘ、Ｙ表示计算点坐标；Ｅ表示边桩到计算点的距离；Ｊ表示边桩与线路的夹角；ＥＸ、ＥＹ表示边桩坐标。

测量实用计算程序（适用于CASIO 4500/4800）一）线路坐标计算程序采用的公式1、曲线要素公式：内移值p= L s2/（24R）切线增量q= L s/2- L s3/（240R2）切线长T=（R+P）tg(A/2)+ q切线角β0= L s/(2R)曲外距E=（R+P）SEC(A/2)-R曲线长L=R(A-2β0)π/180+2 L s2、缓和曲线上任一点与ZH点之间的弦长公式：C=L-L5/（90R2L s2）3、缓和曲线任一点的偏角δ公式：δ=（L/ L s）2δ0δ0=1/3β0（其中L为弧长R为半径Ls为缓和曲线长A为曲线偏角）4、其它均采用座标增量公式：ΔX=LCOSθΔY=Lsinθ然后与JD的座标相加：X=X JD+ΔXY=Y JD+ΔY其中：L为计算点至交点JD的长度θ为JD至计算点的方位角X JD不交点JD的X座标Y JD不交点JD的Y座标圆曲线增设缓和曲线示意图程序功能：可计算直线、缓和曲线、圆曲线上任一点的中桩座标及其法向边桩座标。

A ：曲线偏角 T ：切线距 E ：外矢距 R ：圆曲线半径ZH ：直缓点 HY ：缓圆点 QZ ：曲中点 YH ：圆缓点 HZ ：缓直点JD ZH HZ HY YH AO P q T E β0 R QZ一）单交点对称型曲线坐标计算程序（CASIO 4500/4800）1.Defm 62.A:C”LS=”:D”JD=”:R”R=”3.P=Cx y2/24/R4.Q=C/2-C x y3/240/R x y25.B=90C/R/л6.T=(R+P)tg(absA/2)+Q7.E=(R+P)/cos(A/2)-R8.L=(absA-2B) ЛR/180+2C9.G”ZH”=D-T◢H”HY”=G+C10.z[6]”QZ”=G+L/211.K”YH”=G+L-C12.M”HZ”=G+L13.N”JD X=”:W”JD Y=”:F”FANGWEIJIAO”:J14.A<0→S=-1:≠→S=1⊿15.U=F+A/2+90S16.V=E+R17.B=N+VcosU18.O=W+VsinU19.20.X”X0=”：Y”Y0=”21.Z[1]=X22.Z[2]=Y23.24.Lb1 125.{Z}26.Z≤G→L=D-Z27.V=F+18028.U=F29.Goto 2⊿30.Z≤H→L=Z-G31.V=L-L x y5/90/(RC) x y232.L=30SL x y2/R/л/C33.P=F+18034.Q=F+L35.U=F+3L36.Goto 4⊿37.Z≤K→L=F+A/2+180+90S+180S(Z-I)/R/л38.U=L+90S39.Goto 5⊿40.Z≤M→L=M-Z41.V=L-L x y5/90/(RC) x y242.L=30SL x y2/(RлC)43.P=F+A44.Q=P+180-L45.U=Q-2L-18046.Goto 4⊿47.Z>M→L=Z-M+T48.V=F+A49.U=V50.Goto 2⊿51.Lbl 252.X=N+LcosV◢53.Y=W+LsinV◢54.Goto 6⊿55.Lbl 356.{E}57.P”XL”=X+Ecos(U-90)◢58.Q”YL”=Y+Esin(U-90)◢59.Z[3]=P-Z[1]60.Z[4]=Q-Z[2]61.Pol(Z[3],Z[4])62.J≥0→J”FWJ”=J◢63.≠→J”FWJ”=J+360◢⊿64.Z[5]”L”=I◢65.P”XR”=X+Ecos(U+90)◢66.Q”YR”=Y+Esin(U+90)◢67.Z[3]=P-Z[1]68.Z[4]=Q-Z[2]69.Pol(Z[3],Z[4])70.J≥0→J”FWJ”=J◢71.≠→J”FWJ”=J+360◢⊿72.Z[5]”L”=I◢73.Goto 174.Lbl 475.X=N+TcosP+VcosQ◢76.Y=W+TsinP+VsinQ◢77.Goto 678.Lbl 579.X=B+RcosL◢80.Y=O+RsinL◢81.Goto 682.Lbl 683.J=1→Goto 3⊿84.Goto 1说明：A 曲线偏角（左偏为负数）LS 缓和曲线长度D JD交点里程R 圆曲线半径N JD点X坐标W JD点Y坐标F 起算方位角ZH至JDJ 输入1时计算中桩两侧坐标输入0时只计算中桩坐标E 输入中桩至两侧的距离XL YL 为左侧坐标XR YR 为右侧坐标注：针对具体一个曲线可将已知条件数据直接输入程序中。

4800及4850计算器常用程序ZBZS（坐标正算）E“X1=”:F“Y1=”:G“T1=LbI 0：Fixm：{AS}：G>360=>G=G-360 “T2=”:G DMSE=E+ScosG：“X2=”：F=F+SsinG: “Y2=”X1：起算点的X坐标Y1：起算点的Y坐标T1：起算边指向计算前进方向的坐标方位角A：计算前进方向的左角ZBFS（坐标反算）A“X1=”：B“Y1=”：C“X2=”：D“Y2=”：POL（C-A，D-B）：Fixm：J>0=>W=J:≠“W（1-2）=”“S（1-2）=”：XLYH（线路优化程序）A“XA=”：B“YA=”：C“TA=”：E“1÷RA，Z-，Y+”：P“KA=LbI1: ｛FQ｝:F“1÷RB，Z-，Y+”：Q“KB=LbI2:｛VDZ｝:V“K=”:Z“JIA-JIAO=”:D“WAIYI=(Z-,Y+)”: O=V:O<P=>Goto 3: ≠LbI3: “K-SMALL”LbI4: O>Q=>Goto 5: ≠LbI5: “JISUAN:A-TXY”LbI6: T=C+90+3.1415926×(2E+(F-E)(V-P)÷(Q-P))(V-P):Rad:W=A+∫(cos(C0+0.5×(2E+(F-E)(X-P)÷(Q-P))(X-P), P,V,4)+D×cos(T0+ Z0)Y=B+∫(sin(C0+0.5×(2E+(F-E)(X-P)÷(Q-P))(X-P)),P,V,4)+D×sin(T0+ Z0):LbI 7: O≤Q=>Goto 9:≠LbI 9: W“X=Y“Y=”T“T=”XA：起算点的X坐标YA：起算点的Y坐标TA：起算点的切线方位角1÷RA，Z-，Y+：起算点的曲率（1/R），左偏曲线为正，右偏曲线为负KA：起算点的里程1÷RB，Z-，Y+：终算点的曲率（1/R），左偏曲线为正，右偏曲线为负KB：终算点的里程JIA-JIAO：边桩计算时所要输的与线路切向的夹角，计算点与线路正交时输90度，斜交时则输与线路切向右侧所夹的斜交角度WAIYI=(Z-,Y+)：边桩计算时所要输的计算点到中线点的直线距离（即边桩的外移距离），计算点位于线路左侧时为负，位于线路右侧时为正程序说明：本程序计算线路坐标时要将线路分段（一般分直线段、第一缓和曲线段、圆曲线段、第二缓和曲线段、直线段以此循环）进行计算，例如：缓和曲线加圆曲线计算时，必须将曲线分三段进行计算，先输（起点要素）ZH点的坐标XA、YA，ZH到JD的方位角（即ZH点的切线方位角）TA，ZH 点的曲率1÷RA、里程KA，再输（终点要素）HY点的曲率1÷RB、里程KB，然后输入ZH到HY间任意点的里程便可计算；当所输计算点里程大于HY点里程时，程序显示“JISUAN:A-TXY”表示所要计算点超过第一缓和曲线计算范围，此时程序自动默认HY点的曲线要素为圆曲线计算段的起点要素，用户只需再输入（终点要素）YH点的曲率1÷RB、里程KB，就可计算圆曲线上任意点；当计算点里程大于YH 点里程时，程序显示“JISUAN:A-TXY”表示所要计算点超过圆曲线计算范围，程序自动默认YH点的曲线要素为第二缓和曲线计算段的起点要素，用户只需再输入（终点要素）HZ点的曲率1÷RB、里程KB，就可计算第二缓和曲线上任意点。

主程序F3LbI 1:A“X1”:B“Y1”:C“X2”:D“Y2”N=0:E=C-A:F=D-B ↙ Prog“F1”↙{P}:P≥1 Prog“F2”△↙LbI 6:{XY}:E=X-A:F=Y-B:Prog“F1”:N≠1 L=G ↙Z≥L M=Z-L-180: M=Z-L+180: △↙M >0 M=180-M ▲M=180+Abs M▲△ Goto 6 ↙子程序F1E=0 Goto4 △ I=tan-1Abs(F÷E):S=√ (E2+F2) ▲F≥0 Goto2△E<0 G=180+I▲G = 360-I▲△Goto 5↙LbI 2:E<0 G=180-I ▲G = I ▲△ Goto 5 ↙LbI 4:F>0 G = 90 ▲G =270 ▲△ Goto 5 ↙LbI 5:N=N+1:N=1 Z= G △↙子程序F2LbI 9:{SI}:P=G+180-I:X=A+Scosp▲Y=B+SsinP ▲{o}:o≤0 Goto 9 △ o≤1 A=X:B=Y:G=P+180:G>360 G = G-360:Goto 9 △↙说明:1、主程序F3，已知座标放线；2、P<1主程序进入座标放样计算，P≥1时主程序进入座标计算；3、F3中的S是站点至待放点的距离，M是后视点至待放点的左角。

4、F2中的O表示站点，O>1迁站重新输入（X1，Y1）（X2，Y2）计算座标。

输入转角I时；转角左转I取“+”，右转I取“—”。

5、F2中的S，是表示站点至待求座标点的距离，用测距仪测定。

主程序：PQXLbI O:A“K0”:B“X0”:C“Y0”:D“L0”:E“A0”:F“B0”:G“K1”:N :R:T:Z:L=K-A:S=K-A-D:I=G-K ↙{K} ↙N<1 Prog “P1”:Goto O △↙N=1 Prog “P2”:Goto O △↙Prog “P3”↙Goto O ↙子程序：P1Z≤1 X=B+[L-L5÷(40R2D2)]CosE+[L3cos(E+90)]÷(6RD)▲Y=C+[L-L5÷(40R2D2)]sinE+[ L3sin(E+90)] ÷(6RD) ▲Y=C+[L-L5÷(40R2D2)]sinE+[ L3sin(E-90)] ÷(6RD)▲子程序:P2Z≤1 X=B+[D- D3÷(40R2)]cosE+[(D2÷(6R))]cos(E+90)+[S-S3÷(6 R2 )+ S5 /(120R4 )]Cos[E+(90D) ÷(πR )]+[ S2÷(2R)-S4÷(24R3)+S6÷(720R5)]cos[E+(90D)÷(πR)+90] ▲ Y=C+[D- D3÷(40R2)]sinE+[(D2÷(6R))]sin(E+90)+[S-S3÷(6R2)+S5÷(120R4)]sin[E+(90D)÷(πR)]+[S2÷(2R)-S4÷(24R3)+S6÷(720R5)]sin[E+(90D)÷(πR)+90] ▲X=B+[D-D3÷(40R2)]cosE+[(D2÷(6R)]cos(E-90)+[S-S3÷(6R2) +S5÷(120R4)]cos[E-(90D)÷(πR)]+[S2÷(2R)-S4÷(24R3)+S6÷(720R5)]cos[E-(90D)÷(πR)-90] ▲Y=C+[D-D3÷(40R2)]sinE+[(D2÷(6R)]sin(E-90)+[S-S3÷(6R2) +S5÷(120R4)]sin[E-(90D)÷(πR)]+[S2÷(2R)-S4÷(24R3)+S6÷(720R5)]sin[E-(90D)÷(πR)-90] ▲子程序:P3Z≤1 X=B+TcosE+Tcos(E+F)+[I-I5÷(40R2D2)]cos(E+F+180) +[I3÷(6RD)]cos(E+F+180-90)▲Y=C+TsinE+Tsin(E+F)+[I-I5÷(40R2D2)]sin(E+F+180) +[I3÷(6RD)]sin(E+F+180-90)▲X=B+TcosE+Tcos(E-F)+[I-I5÷(40R2D2)]cos(E-F+180) +[I3÷(6RD)]cos(E-F+180+90)▲Y=C+TsinE+Tsin(E-F)+[I-I5÷(40R2D2)]sin(E-F+180)+[I3÷(6RD)]sin(E-F+180+90)▲说明:1.K0为ZH点桩号;(X0，Y0)为ZH点座标；L0为缓和曲线段长度，A0为直线方位角；B0为转角；K1为HZ点桩号；N为条件参数；R为半径；T为切线长；Z为条件参数；K为待求点桩号。

1.主程序(TYQX js)"1.SZ => XY"："2.XY => SZ"：N：U"X0"：V"Y0"：O"S0"：G"F0"：H"LS"：P"R0"：R"RN"：Q：C=1P：D=(P-R)(2HPR)：E=180π：N=1=>Goto 1：≠>Go to 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◢Z"Z"=Z◢Goto 22. 正算子程序(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(C+F WD))+Acos(G+QEMW(C+MWD)))：Y=V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+QEFW(C+FW D))+Asin(G+QEMW(C+MWD)))：F=G+QEW(C+WD)+90：X=X+ZcosF：Y=Y+ZsinF3. 反算子程序(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二、使用说明1、规定(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。

支导线程序(fx4800P)ZDXProg”006”:A=J:Z”1—2”:X=U:Y=W:LbI0:{Q}:Q”Q=”:{L}:L”L=”:X=X+LCos(A+Q):Y=Y+LSin(A+ Q): X”X=”▲ Y”Y=”▲A=A+Q+180:Z=2=>X=U:Y=W:A=J△Goto0放样程序(fx4800P)FYProg”006”: LbI0:{X}:X”X=”:{Y}:Y”Y=”: Prog”007” :Goto0竖曲线程序(fx4800P)SHU--QXI“I=”:J”J=”:R”R=”:W=I-J:C=Abs(RW)/1000:T=C/2:E=T^2/(2R): C”L=”▲ T”T=”▲ E”E=”▲W<0=>K=1△W>0=>K=-1△D”BIAN—LI=”:Z”BIAN—H=”:A=Z+(-IT)/1000:LbI0:{L}:L”DK=”:L≤D-T=>H=Z+(D-L)(-I/1000)△L>D-T=>H=A+I(L-(D-T))/1000+K(L-(D-T))^2/(2R) △L≥D+T=>H=Z+J(L-D)/1000△H”H=”▲Goto0已知点坐标求对应里程和距离(fx4800P)1、主程序(DYLC)Defm42: AB”N1=”C”N2=”:R:Pro g”000” :Pro g”009”:E”PYL=”:K”Z—Y”: Z”ZH=”:M”ZH—X”:P”ZH—Y”:F:Lbl0:Prog”XH”:Goto02、子程序(XH)Lbl1:Pro g”QX”:{U}:U”YZ—X=”:{W}:W”YZ—Y=”:I=U-X:J=W-Y:Pol(I,J):J<0=>J=J+360△Z[37]=J-Z[34]:Z[38]=cosZ[37]×I:D=D+Z[38]:AbsZ[38]≥0.001=>Goto1:≠> D”DYLC=”△▲Z[35]=sinZ[37] ×I: Z[35] “FXJL=”▲已知3点坐标求圆心和半径(fx4800P)主程序(3DZB--QYX)A”AX=”B”AY=”C”BX=”D”BY=”E”CX=”F”CY=”: X=((F-B)(D^2-B^2+C^2-A^2)+(D-B)(B^2-F^2+A^2-E^2))/(2(C-A)(F-B)-2(E-A)(D-B)):X”X=”▲Y=((E-A)(C^2-A^2+D^2-B^2)+(C-A)(A^2-E^2+B^2-F^2))/(2(D-B)(E-A)-2(F-B)(C-A)):Y”Y=”▲R=√((A-X)^2+(B-Y)^2):R“R=”▲自由设站(fx4800P)(边角后方交会)HFJHLbI0:{A}:A”H1—X=”:{B}:B”H1—Y=”:{C}:C”H2—X=”:{D}:D”H2—Y=”:I=C-A:J=D-B:POL(I,J):J≤0=>J=J+360⊿ I:”S=”▲Z[42]=J: Prog”006”:Z[42]:”FW—Q=”▲{E}:E”L1=”:{F}:F”L2=”:{G}:G”Q=”:H=Cos-1((I2+ E2-F2)/(2IE)):K=Cos-1((I2+ F2-E2)/(2IF)):L=A+E×Cos(J+H) :M=B+E×Sin(J+H):N=C+F×Cos(J+180-K):O=D+F×Sin(J+180-K):X=(L+N)/2:Y=(M+O)/2:X:”X=”▲Y:”Y=”▲Goto0不对称缓和曲线放样数据计算(fx4800P)(带偏移量+中边桩)1 主程序(--QX—Z-B)Defm42:AB”N1=”C”N2=”:R:Pro g”000”: Pro g”009”:E”PYL=”:K”Z—Y”: Z”ZH=”:M”ZH—X”:P”ZH—Y” :F Lbl0:Pro g”QX”: Prog”006”:Prog”007”:Goto02 子程序(QX){D}:Z[1]=D”DK=”-Z:O=Z+Z[2] :N=(M+Z[11]CosF)+Z[12]Cos(F+KA) :S=(P+Z[11]SinF)+Z[12]Sin (F+KA):{L}:L”BIN—L”:{Q}:Q”BIN—Q”:Z[1]≤0=>Prog”001”:≠>Z[1]≤B=>Prog”002”:≠>Z[1]≤(Z[2]-C)=>Prog”003”:≠>Z[1]≤Z[2]=>Prog”004”:≠>Z[1]>Z[2]=>Porg”005”△△△△△ X”X=”▲ Y”Y=”▲Z[34]<0=>Z[34]=Z[34]+360:≠>Z[34]>360=> Z[34]=Z[34]-360△△Z[42]=Z[34]:Pro g”010”:Z[42]”QX—Q”▲(该点切线方位角)3 子程序（000）曲线要素计算Z[4]=B/2-B^3/(240R^2):Z[6]=B^2/(24R)-B^4/(2688R^3):Z[8]=90B/(πR）:Z[5]=C/2-C^3/(240R^2):Z[7]=C^2/(24R)-C^4/(2688R^3):Z[9]=90C/(πR) :Z[11]=(R+Z[6])tng(A/2)+Z[4]-（Z[6]-Z[7]）/sinA:Z[12]=(R+Z[7])tng(A/2)+Z[5]+(Z[6]-Z[7])/sinA:Z[10]=(R+(Z[6]+Z[7])/2)(1/cos(A/2)) -R:Z[2]=(A-Z[8]-Z[9])Rπ/180+B+C: Z[3]=Z[2]-B-C4 子程序（001）第一切线段计算X=M+Z[1]cosF:Y=P+Z[1]sinF:X=X+Lcos(F+Q):Y=Y+Lsin(F+Q): Z[34]=F:5 子程序（002）第一缓和曲线段坐标计算Z[28]=F:Z[27]=E:Z[13]=90B/（πR）：Z[14]=Z[13]/3：Z[15]=90Z[1]^2/(πRB):Z[16]=Z[14](Z[1]/B)^2 :Z[17]=RB:Z[18]=Z[1]-Z[1]^5/(40Z[17]^2)+Z[1]^9/(3456Z[17]^4):Z[19]=Z[1]^3/(6Z[17])-Z[1 ]^7/(336Z[17]^3)+Z[1]^11/(42240Z[17]^5):Z[20]=√(Z[18]^2+Z[19]^2):X=M+Z[20]cos(Z[28]+KZ[16]) :Y=P+Z[20]sin(Z[28]+KZ[16]):K<0=> Z[27]= Z[27]:≠>K≥0=> Z[27]=-1 Z[27]△△Z[23]= Z[27] (Z[1]/B):Z[21]= Z[28]+KZ[15]+90: X=X+Z[23]cosZ[21]:Y=Y+Z[23]sinZ[21]:Z[22]=Z[28]+KZ[15]+Q : X=X+LcosZ[22]:Y=Y+LsinZ[22]: Z[34]=Z[28]+KZ[15]6 子程序（003）圆曲线段坐标计算Z[27]=E:Z[28]=F:Z[13]=90B/(πR):Z[22]=Z[28]+K(Z[13]/3):Z[17]=RB:Z[18]=B-B^5/(40Z[17]^2)+B^9/ (3456Z[17]^4):Z[19]=B^3/(6Z[17])-B^7(336Z[17]^3+B^11/(42240Z[17]^5):Z[20]=√(Z[18]^2+Z[19]^2) :X=M+Z[20]cosZ[22]:Y=P+Z[20]sinZ[22]:Z[14]=Z[1]-B:Z[15]=180Z[14]/( πR):Z[23]= Z[28]+KZ[13] :Z[21]=2Rsin(Z[15]/2):Z[16]=Z[15]/2:X=X+Z[21]cos(Z[23]+KZ[16]):Y=Y+Z[21]sin(Z[23]+KZ[16]):Z[24]= Z[23]+KZ[15]:K<0=> Z[27]= Z[27]:≠>K≥0=> Z[27]=-1 Z[27]△△X=X+ Z[27]cos(Z[24]+90):Y=Y+ Z[27]sin(Z[24]+90):X=X+Lcos(Z[24]+Q):Y=Y+Lsin(Z[24]+Q): Z[34]＝Z[28]+KZ[13]+KZ[15]7 子程序（004）第二缓和曲线段坐标计算Z[15]=O-D:Z[22]=F+KA+180:Z[27]=E:K>0=>Z[27]=Z[27]:≠>K<0=> Z[27]=-1Z[27]△△Z[25]=Z[27]( Z[15]/C):Z[23]=-K:Z[13]=90C/(πR):Z[14]=Z[13]/3:Z[16]=90Z[15]^2/(πRC):Z[17]=Z[14](Z[15]/C)^2 :Z[18]=RC:Z[19]=Z[15]-Z[15]^5/(40Z[18]^2)+Z[15]^9/(3456Z[18]^4):Z[20]=Z[15]^3/(6Z[18])-Z[15]^7 /(336Z[18]^3+Z[15]^11/(42240Z[18]^5):Z[21]=√(Z[19]^2+Z[20]^2):Z[30]=N+Z[21]cos(Z[22]+Z[23]Z[17]): Z[31]=S+Z[21]sin(Z[22]+Z[23]Z[17]):Z[24]=Z[22]+Z[23]Z[16]+90: Z[32]= Z[30]+Z[25]cosZ[24]: Z[33]= Z[31]+Z[25]sinZ[24]:Z[26]=F+KA+Z[23]Z[16]+ Q : X= Z[32]+LcosZ[26]:Y= Z[33]+LsinZ[26]: Z[34]＝Z[22]+（-1K）Z[16]-1808 子程序(005)计算第二切线段坐标Z[13]=F+KA:Z[14]=D-O:X=N+Z[14]cosZ[13]:Y=S+Z[14]sinZ[13]:X=X+Lcos(Z[13]+Q):Y=Y+Lsin(Z[ 13]+Q): Z[34]=Z[13]9 程序(006)计算后视数据{U}:U”ZX=”:{W}:W”ZY=”:{G}:G”HX=”:{H}:H”HY=”:I=G-U:J=H-W:Pol(I,J): I”S=”▲Prog”008”: Z[42]=J: Prog”010”: Z[42]”F0=”▲10 子程序(007)计算放样数据T=J:Pol(X-U, Y-W): Prog”008”: I”FY--L=”▲V=J-T:V<0=>V=V+360:≠>V>360=>V=V-360△△Z[42]=V: Prog”010”: Z[42]”FY--Q=”▲11 子程序(008)角度判断J<0=>J=J+360:≠>J>360=>J=J-360△△12 子程序（009）曲线要素结果显示Z[11]”T1=”▲ Z[12]”T2=”▲ Z[10]”E=”▲ Z[2]”QX--L=”▲ Z[3]”YQX—L=”▲13 子程序(010)角度的度分秒输出Z[40]=60FracZ[42]:Z[41]=60FracZ[40]:Z[42]=IntZ[42]+IntZ[40]÷100+Z[41]÷10000。

卡西欧4800\4850计算器计算程序及说明一：坐标反算：(ZBFS)1 程序步骤：C”X1”:D”Y1”:E”X2”:F”Y2”:Fixm:Pol(E-C,F-D:I”S1-2=”▲J≤Ｏ＝＞J=J+360△J”A1-2=”2 操作过程：ZBFS→EXE→输入X1值→EXE→输入Y1值→EXE→输入X2值→EXE→输入Y2值→EXE→EXE→EXE注：①：S1-2:计算得出的距离。

②：A1-2:计算得出的角度。

（按shift°′″转换为60进制的角度）③：此程序可循环计算。

3 案例：已知：D1坐标（7811.23，606.136），D2坐标（7805.915，431.910）求解： a D1-D2(A1-2)=268°15′09.56 ″S D1-D2(S1-2)=174.3071二：坐标正算：(ZBZS)1 程序步骤：C”X1”:D”Y1”:S”S1-2”:R”A1-2”:Fixm:X”X2”=C+Rec(S,R) ▲Y”Y2”=D+J2 操作过程：ZBZS→EXE→输入X1值→EXE→输入Y1值→EXE→输入S1-2距离值→EXE→输入A1-2角度（例268°15′09.56 ″）值→EXE→EXE注：①：X2计算得出坐标，Y2计算得出坐标。

②：此程序可循环计算。

3 案例：已知：D1坐标（7811.23，606.136），两点距离：S D1-D2(S1-2)=174.3071方位角：a D1-D2(A1-2)=268°15′09.56 ″ 求解：D2坐标（7805.915，431.910三：线路曲线计算程序:（XLQXJS）(万能公式)1 程序步骤：1.Lbl 0:{E G}：A“XA”：B“YA”：C“CA”：D“1/RA”：E“1/RB”：F“DKA”： G“DKA”2.Lbl 1:{H O R}:H“DKI”：O“DL”：R“DR”：H＞G=﹥Goto 3⊿3.P=(E-D)/Abs(G-F):Q=Abs(H-F):I=P×Q:T=D+I4.J=C+(I+2D)×Q×90/π ▲5.M=C+(I/4+2D)Q×45/(2π):N=C+(3I/4+2D)Q×135/(2π):6.K=C+(I/2+2D)Q×45/π7.X=A+Q(cosC+4(cosM+cosN)+2cosK+cosJ)/12 ▲8.Y=B+Q(sinC+4(sinM+sinN)+2sinK+sinJ)/12 ▲9.U“XL”=X+Ocos(J-90) ▲10.V“YL”=Y+Osin(J-90) ▲11.W“XR”=X+Rcos(J+90) ▲12.Z“YR”=Y+ Rsin(J+90) ▲13. A=X：B=Y：C=J：D=T：F=H：Goto 1Lbl 3: A=X：B=Y：D=E：F=G：C=J：Goto 0程序说明：A- 曲线元起点A的坐标；B- 曲线元起点B的坐标；C- 曲线元起点A的切线坐标方位角；F- 曲线元起点A的里程；G- 曲线元起点B的里程；H- 曲线上待求点i的里程；D- 曲线元起点A的曲率；E- 曲线元终点B的曲率；XL-左边线点位X坐标；YL-左边线点位Y坐标；XR-右边线点位X坐标；YR-右边线点位Y坐标；X- 中线点位纵坐标；Y- 中线点位横坐标；DL-左边线距中线平距；DR-右边线距中线平距；该程序需要输入的数据为：（1）.曲线元起点A的坐标及切线坐标方位角，计算器上用“XA”，“YA”，“CA”显示；（2）.曲线元起点A和B的曲率，计算器上用I÷RA，I÷RB 显示（曲线左偏时取“-”）；（3）.曲线元起点A和终点B的里程，计算器上用“DKA”，“DKB”显示；（4）.输入待求点里程和该点距左右的水平距离，计算器上用“DKI”，“DL”，“DR”显示；每算完一个待求点的中线及边线坐标，程序又让输入下一点的“DKI”，“DL”，“DR”当输入的“DKI”大于“DKB”时程序中显示“1÷RB”和“DKB”，此时输入下一个曲线元起点的曲率和里程，然后重复步骤（4），即可计算下一个曲线中线及边线点位坐标。

各位同行,本人领悟yshf的程序后,进行了以下改编,达到以下目标:1.加编数据库,计算中不必逐项输入"线元要素",提高运算速度，避免现场忙中出错2.计算直观,人性化3.正算直接输入里程和边距,反算输入近似里程便可4.增加了“计算点与测站点”的距离和方位角计算语句，方便直接放样5.愿收获与大家共享一.改动后的程序清单:1.主程序(TYQXJS)"1.SZ => XY"："2.XY => SZ"：U"X0"：V"Y0"：O"S0"：G"F0"：H"LS"：P"R0"：R"RN"：Q：D=(P-R)÷(2HPR)：N=1=>Goto 1：≠>Goto 2Δ←┘Lbl 1：{Z}：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=” ◢ΔProg"SJK"←┘Lbl 2：{XY}：XY：I=X：J=Y：Prog "SUB2"：S"S"=O+W◢Z"Z"=Z◢Prog"SJK"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+57.2958QKW(1/P+KWD))+Bcos(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+ZcosF：Y=Y+ZsinF3. 反算子程序(SUB2)T=G-90：W=Abs((Y-V)cosT-(X-U)sinT)：Z=0：Lbl 0：Prog "SUB1"：L=T+Q57.2958W(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二.增设数据库程序（ＳＪＫ）Lbo4：｛S｝：S∠下一线元起点里程≡＞O=本线元起点里程：U=本线元起点X：V=本线元起点Y：G=本线元起算方位角：H=本线元长度：P=起点曲率半径：R=终点曲率半径：Q=0或1、-1：Prog“TYQXJS”：Goto4Δ←┘S∠下一线元起点里程=＞O=。

十一．BI DAO XIAN（闭导线）1．用途：该程序是“闭合导线平差”的计算程序。

2．程序数学模型：3．程序清单：Defm 40：X“X0”：Y“Y0”：F“A0”：N“UNKNOWN POINT”：N=N+1↙I=0：E=0：C=0：Lbl 0 ↙I=I+1：{L，D}：Z[2I-1]=L：Z[2I]=D：E=E+L：C=C+D：I<NGoto 0 Δ B“FW”=E-180（N-2）◢R：A“FR”=RN◢B=-B÷N：I=1：Z[2N+1]=F：Lbl 1 ↙I+I+1◢F=F+Z[2I-1]：F180F=F+180：F=F-180 Δ F360F+F-360 ΔZ[2N+2I-1]“A”=F+（I-1）B◢I<NGoto 1Δ I=0：G=0：H=0：Lbl 2 ↙I+I+1：Z[2I-1]=Z[2I]cosZ[2N+2I-1]：G=G+Z[2I-1]：Z[2N+2I]=Z[2I]sin Z[2N+2I-1]：H=H+Z[2N+2I]：I<N Goto 2 Δ G“FX”◢H“FY”◢I“F”=（G2+H2）◢K=I÷C◢U=-G÷C：V=-H÷C↙I=0：Lbl 3 ↙I+I+1◢Z[2I-1]=Z[2I-1]+U Z[2I]：Z[2N+2I]=Z[2N+2I]+V Z[2I]：X=X+ Z[2I-1]◢Y=Y+ Z[2N+2I]◢I<N-1Goto 3：“END”4．程序说明：X0，Y0，A0——已知的起始点的纵、横坐标及起始边的坐标方位角。

UNKNOWN POINT——待求未知点的个数。

L，D——观测的导线各左转角及平距。

FW，R，FR——导线角度闭合差、角度闭合差的允许误差的系数及允许误差值。

I，A——各条导线边序号及其坐标方位角。

FX，FY——坐标增量ΔX、ΔY的闭合差。

F，K——导线全长闭合差、导线全长相对闭合差。