4800计算器常用测量计算程序

fx—4800PA计算器在工程测量中的应用编写测量坐标计算程序考虑到施工的便捷性以及现有工具，我们采用fx—4800PA计算器进行编程，这样既可以保证测量的准确又可使施工测量便捷，编程如下：1: F 坐标通用公式Z=0=>GOTO 1:≠>GOTO 2LbI 1K＝F：L=Q-A:X=M+LcosF◢Y=N+LsinF◢GOTO 0LBI 2S:W=0=>GOTO 3:≠>GOTO 6LBI 3L=Abs(Q-A):I=L-L^5/40(RS)2:J=L^3/6RS-L^7/336(RS)^3:T=0=>GOTO 4:≠>GOTO5LBI 4K＝F-90GL^2/RSπ：X=M+GIcosF+JsinF◢Y=N+GIsinF-JcosF◢GOTO 0LBI 5K＝F+90GL^2/RSπ：X=M+GIcosF-JsinF◢Y=N+GIsinF+JcosF◢GOTO 0LBI 6L=Q-A:D=180L-90S:I=Rsin(D/Rπ)+S/2-S^3/240R^2:J=R(1-cos(D/R π))+S^2/24R:P=0=>GOTO 7:≠>GOTO 8LBI 7K＝F-180（L-S）/Rπ：X=M+GIcosF+JsinF◢Y=N+GIsinF-JcosF◢LBI 8K＝F+180（L-S）/Rπ：X=M+GIcosF-JsinF◢Y=N+GIsinF+JcosF◢GOTO 0LBI 0U“ZX”＝X+Bcos(K-90) ◢(90度可以变为任意角度)V“ZY”＝Y+Bsin(K-90) ◢(90度可以变为任意角度)E“YX”＝X+Ccos(K+90) ◢(90度可以变为任意角度)O“YX”＝Y+Csin(K+90) ◢(90度可以变为任意角度)说明：Q—待测点里程 A—起点里程X--待测点X坐标：Y--待测点Y坐标M—起点X坐标： N—起点Y坐标：L—直线段待测点至起点距离 F—起点方位角S—缓和曲线总长 R—圆曲线半径 K--待测点切线方位角B—左边桩点至中心距离 C—左边桩点至中心距离U“ZX”—左边桩X坐标： V“ZY”—左边桩Y坐标：E“YX”—右边桩X坐标： O“YX”—右边桩Y坐标：Z=0—直线段坐标计算Z≠0，W＝0，T=0,G=1—左转角ZH—HY和右转角HZZ≠0，W＝0，T=0,G=-1—左转角YH—HZZ≠0，W＝0，T≠0,G=1—右转角ZH—HY和左转角HZZ≠0，W＝0，T≠0,G=-1—右转角YH—HZZ≠0，W≠0，P=0,G=1—左转角HY—YH或ZY—YZZ≠0，W≠0，P≠0,G=1—右转角HY—YH或ZY—YZ该程序能够计算出曲线，对称及不对称缓和曲线任意桩号的中桩坐标、征地线边桩坐标及中桩前进方向切线顺时针旋转任意角度至中桩任意距离的点的坐标。

FX4800P常用测量程序集1、由一个已知点坐标和一个已知方位角计算支导线坐标（程序名：ZBJS）程序清单：X：Y：T：Lb1 0：Fixm：{AS}：T=T+A-180：T〈0 =〉T=T+360⊿T〉360 =〉T =T-360⊿T◢X=X+Scos T◢Y=Y+Ssin T◢Goto 0程序运行说明：运行“ZBJS”，每运行一步按EXE键一次，显示X？时输入已知点X坐标，显示Y？时输入已知点Y坐标，显示T？时输入已知点由后视到前视的已知坐标方位角，运行一次本程序只需输入一次已知数据。

显示A？时输入第一个左转角值，显示S？时输入测站到前视的水平距离。

接下来计算器显示的T是测站到前视的坐标方位角，按+ 可将显示的十进制角度换算为度分秒。

如此继续可计算出所有支导线未知坐标。

2、由两个已知点坐标计算支导线坐标（程序名：ZBJS2）程序清单：M“X0”：N“Y0”：D“XH”：E“YH” ↲Pol(M-D，N-E)：Fixm：T=J：X=M：Y=N：T〈0=〉T=T+360⊿T“TH-J=” ◢I“S=”◢Lb1 0:{AS}: T=T+A-180：T〈0 =〉T=T+360⊿T〉360 =〉T =T-360⊿T◢X=X+Scos T◢Y=Y+Ssin T◢Goto 0程序运行说明：运行后，显示X0，Y0时输入已知置镜点坐标，显示XH，YH时输入已知后视点坐标，其它同ZBJS的输入与显示。

3、由一个已知点坐标和一个已知方位角进行放射型导线坐标计算（程序名FSXZB JS）程序清单：C“X0”：D“Y0”：E“T0”：Fixm：Lb1 1:{SA}↲T=E+A-180：T〈0=〉T=T+360⊿T〉360 =〉T=T-360⊿T◢X=C+Scos T◢Y=D+Ssin T◢Goto 1程序运行说明：变量输入参照ZBJS程序。

4、由两个已知点坐标进行放射型导线坐标计算（程序名FSXZBJS2）M“X0”：N“Y0”：D“XH”：E“YH” ↲Pol(M-D，N-E)：Fixm：F=J：X=M：Y=N：F〈0=〉F=F+360⊿F“TH-J=” ◢I“S=” ◢Lb1 1:{SA}↲T=F+A-180：T〈0=〉T=T+360⊿T〉360 =〉T=T-360⊿T◢X=X+Scos T◢Y=Y+Ssin T◢Goto 15、坐标计算优化程序（可进行连续式（LX）、放射式（FS）已知一点（1-D）或者两点（2-D）进行坐标计算）：V“1=>LX，≠>FS”：W“1 =>1-D≠>2-D”：W=1=> Goto 1:≠> Goto2↲Lb1 1：F“X0=”：G“Y0=”：H“T H-0=”：Goto 3↲Lb1 2:F“X0=”：G“Y0=”：M“XH=”：N“YH=”:Pol(F-M，G-N)：Fixm: J〈0=〉J=J+ 360⊿H=J：H“TH-0=”◢I“S-H=”◢Goto 3↲Lb1 3: X=F：Y=G：T=H ：V=1=> Goto 4:≠> Goto 5↲Lb1 4: {AS}: A:S:T=T+A-180：T〈0 =〉T=T+360⊿T〉360 =〉T =T-360⊿T◢X=X+Scos T◢Y=Y+Ssin T◢Goto 4Lb1 5：{AS}:A:ST=H+A-180：T〈0=〉T=T+360⊿T〉360 =〉T=T-360⊿T◢X=F+Scos T◢Y=G+Ssin T◢Goto 56、已知置镜点坐标和置镜点到后视点坐标方位角求放样角度和放样距离的放样程序（程序名FY1）程序清单：F“X0”：G“Y0”：H“T0-H”：Fixm：Lb1 1:{XY}↲Pol(X-F，Y-G) ↲J〈0 =〉J=J+360⊿A=J-H↲A〈0 =〉A=A+360⊿A“A=” ◢I“S=” ◢Goto 1程序运行说明：变量输入参照坐标计算，A=为放样顺时针角度，S=为放样距离。

卡西欧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），即可计算下一个曲线中线及边线点位坐标。

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÷（240R2）▲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^4SR））：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－Wsin 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=Vsin 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点桩号。

CASIO fx－4500P与CASIO fx－4800P计算器部分按键及程序语句不同之处比较【坐标反算程序】文件名：ZBFS（调用子程序文件名：Z.FS）L1 Lbl 0L2 {XY}:X″X1″:Y″Y1″:U″X0″:Z″Y0″L3 Prog ″Z.FS″L4 ″DIST=″:D:Pause 0L5 D=0=>Goto 0:◣A″ANGLE=″◢Goto 0使用说明屏幕提示执行操作屏幕显示计算结果X1? 输入视点坐标XDIST= 置镜点与视点的距离Y1? 输入视点坐标YX0? 输入置镜点坐标XANGLE= 置镜点→视点方位角Y0? 输入置镜点坐标Y注意fx-4500P型计算器将L3行改成：Prog Z.FS【坐标正算程序】文件名：ZBZSL1 Lbl 0L2 {DA}:D″DIST″:A″ANGLE″:U″X0″:Z″Y0″L3 ″X=″:X=U+Dcos A:Pause 0L4 Y=Z+Dsin A◢Goto 0使用说明屏幕提示执行操作屏幕显示计算结果DIST? 输入置镜点与视点的距离X= 视点坐标X ANGLE? 输入置镜点→视点方位角X0? 输入置镜点坐标XY= 视点坐标Y Y0? 输入置镜点坐标Y 【转角放样程序】文件名：ZJFY（调用子程序文件名：Z.FS）L1 Lbl 0L2 {STUZ}:X=S″HSX″:Y=T″HSY″:U″X0″:Z″Y0″L3 Prog ″Z.FS″L4 D=0=>Goto 0:◣J=AL5 Lbl 1L6 {XY}:X″QSX″:Y″QSY″L7 Prog ″Z.FS″L8 ″DIST=″:D:Pause 0L9 D=0=>Goto 1:◣A≥J=>A″R-ANGLE″=A-J◢≠=>A″R-ANGLE″=A-J+360◢◣Goto 1使用说明屏幕提示执行操作屏幕显示计算结果HSX? 输入后视点坐标XDIST= 置镜点与前视点距离HSY? 输入后视点坐标YX0? 输入置镜点坐标XY0? 输入置镜点坐标YR-ANGLE=置镜点→后视点到置镜点→前视点的转角QSX? 输入前视点坐标XQSY? 输入前视点坐标Y注意fx-4500P型计算器将L3、L7行改成：Prog Z.FS【单交点基本型（对称、非对称）平曲线要素、中桩、边桩坐标计算程序】文件名：PQXJS（调用子程序文件名：PQXZB）L1 R:H″LS1″:H=0=>I=0:◣I″LS2″:Z:D″JD″:S=Abs Z L2 N=H2÷24R-H^4÷2384R^3L3 P=I2÷24R-I^4÷2384R^3L4 K=H÷2-H^3÷240R2L5 M=I÷2-I^3÷240R2L6 B=H÷2R:C=I÷2RL7 T″T1″=(R+P-(R+N)cos S)÷sin S+K◢U″T2″=(R+N-(R+P)cos S)÷sin S+M◢L=H+I+(Sπ÷180-B-C)R◢H=I=>E=(R+N)÷cos (S÷2)-R◢J=2T-L◢◣N″LY″=L-H-I◢B″ZH″=D-T◢C″HY″=B+H◢E″YH″=C+L-H-I◢F″HZ″=E+I◢G″QZ″=C+N÷2◢Prog ″PQXZB″使用说明屏幕提示执行操作屏幕显示计算结果R? 输入圆曲线半径T1= 小里程段切线长LS1?输入小里程段缓和曲线长（不设缓和曲线时，输入0）T2= 大里程段切线长LS2?输入大里程段缓和曲线长（不设缓和曲线时，不提示）L= 平曲线总长Z? 输入转向角(左转负、右转正) E= 外距(T1=T2时显示)JD? 输入交点桩号J= 校正数(T1=T2时显示)LY= 圆曲线长ZH= 直缓点桩号不设缓和曲线时，ZH=HY，为直圆点桩号HY= 缓圆点桩号YH= 圆缓点桩号不设缓和曲线时，YH=HZ，为圆直点桩号HZ= 缓直点桩号QZ= 曲中点桩号XJD? 输入交点的X坐标XZH= 直缓点X坐标（不设缓和曲线时，为直圆点X坐标）YJD? 输入交点的Y坐标YZH= 直缓点Y坐标（不设缓和曲线时，为直圆点Y坐标）F1? 输入小里程段导线的方位角XHZ= 缓直点X坐标（不设缓和曲线时，为圆直点X坐标）YHZ= 缓直点Y坐标（不设缓和曲线时，为圆直点Y坐标）XHY= 缓圆点X坐标不设缓和曲线时，不显示YHY= 缓圆点Y坐标XYH= 圆缓点X坐标YYH= 圆缓点Y坐标DQD? 输入待求点桩号FWJ= 待求点切线方位角X= 待求点中桩X坐标Y= 待求点中桩Y坐标ZB? 输入左幅宽度XZ= 待求点左幅边桩X坐标YZ= 待求点左幅边桩Y坐标YB? 输入右幅宽度XY= 待求点右幅边桩X坐标YY= 待求点右幅边桩Y坐标注意fx-4500P型计算器将L7行最后一句改成：Prog PQXZB；子程序中Prog ″PQX.F″改成Prog PQX.F，Prog ″PQX.BZ″改成Prog PQX.BZ【单交点基本型（对称、非对称）平曲线中桩、边桩坐标计算子程序】文件名：PQXZB（可单独执行，调用子程序PQX.F，PQX.BZ）L1 R:H″LS1″:H=0=>I=0:◣I″LS2″:Z:{XY}:X″XJD″:Y″YJD″:T″T1″:H=0=>U=T:◣U″T2″:B″ZH″:H=0=>C=B:◣C″HY″:E″YH″:H=0=>F=E:◣F″HZ″:{A}:A″F1″:Q=A+Z:Z>0=>S=1:≠=>S=-1:◣ J″XZH″=X+Tcos (A+180)◢K″YZH″=Y+Tsin (A+180)◢P″XHZ″=X+Ucos Q◢O″YHZ″=Y+Usin Q◢H≠0=>V=C-B-(C-B)^5÷40R2H2:W=(C-B)^3÷6RH-(C-B)^7÷336R^3H^3:M″XHY″=J+Vcos A+Wcos (A+90S)◢N″YHY″=K+Vsin A+Wsin (A+90S)◢V=F-E-(F-E)^5÷40R2I2:W=(F-E)^3÷6RI-(F-E)^7÷336R^3I^3:Z ″XYH″=P+Vcos (Q+180)+Wcos (Q+90S)◢Z ″YYH″=O+Vsin (Q+180)+Wsin (Q+90S)◢◣Lbl 0L2 {Z}:Z″DQD″L3 Z<B=>Goto 0:◣Z>F=>Goto 0:◣H=0=>G=A+180(Z-B)S÷πR:Prog ″PQX.F″:M=J:N=K:Goto 1:◣Z≤C=>G=A+90(Z-B)2S÷πRH:Prog ″PQX.F″:V=Z-B-(Z-B)^5÷40R2H2:W=(Z-B)^3÷6RH-(Z-B)^7÷336R^3H^3:X=J+Vcos A+Wcos (A+90S)◢Y=K+Vsin A+Wsin (A+90S)◢Prog ″PQX.BZ″:Goto 0:◣Lbl 1L4 Z≤E=>G=A+90(C-B)2S÷πRH+180(Z-C)S÷πR:Prog ″PQX.F″:V=90H÷πR:W=90(Z-C)÷πR:X=M+2Rsin Wcos (A+SV+SW)◢Y=N+2Rsin Wsin (A+SV+SW)◢Prog ″PQX.BZ″:Goto 0:◣Z>E=>G=Q-90(F-Z)2S÷πRI:Prog ″PQX.F″:V=F-Z-(F-Z)^5÷40R2I2:W=(F-Z)^3÷6RI-(F-Z)^7÷336R^3I^3:X=P+Vcos (Q+180)+Wcos (Q+90S)◢Y=O+Vsin (Q+180)+Wsin (Q+90S)◢Prog ″PQX.BZ″:◣Goto 0【平曲线切线方位角计算结果处理子程序】文件名：PQX.FL1 G>0=>G<360=>G″FWJ=″◢≠>G″FWJ″=G-360◢◣≠>G″FWJ″=G+360◢◣【平曲线边桩坐标计算子程序】文件名：PQX.BZL1 {L}:V=G-90:W″XZ″=X+L″ZB″cos V◢W″YZ″=Y+Lsin V◢{L}:V=G+90:W″XY″=X+L″YB″cos V◢W″YY″=Y+Lsin V◢【契形体体积计算程序】文件名：XXTTJL1 V=H÷6×(A″A1″×B″B1″+C″A2″×D″B2″+(A+C)(B+D)) 【竖曲线要素及高程计算程序】文件名：SQXJSL1 M″BZH″:G″BGC″:A″I1″:B″I2″:RL2 W=B-A◢L=RAbs W◢T=L÷2◢E=T2÷2R◢F″QZH″=M-T◢J″ZZH″=M+T◢Lbl 0L3 {Z}:Z″ZH″L4 Z≤F=>H=G-(M-Z)A◢Goto 0:◣Z≥J=>H=G+(Z-M)B◢Goto 0:◣Y=(Z-F)2÷2RL5 H=G-(M-Z)AL6 W≥0=>H=H+Y◢≠=>H=H-Y◢◣Goto 0 使用说明屏幕提示执行操作屏幕显示计算结果BZH? 输入变坡点桩号W= 前后纵坡坡度差BGC? 输入变坡顶点高程L= 竖曲线长度I1? 输入小里程段纵坡(前进方向上坡为+，下坡为-) T= 竖曲线切线长I2? 输入大里程段纵坡(前进方向上坡为+，下坡为-) E= 竖曲线外距R? 输入竖曲线半径QZH= 竖曲线起点桩号ZH? 输入计算点的桩号ZZH= 竖曲线终点桩号H= 计算点的中桩高程【无竖曲线段纵坡高程连续计算程序】文件名：ZPGC L1 F″ZH0″:G″GC0″:IL2 Lbl 0:{N}:N=0=>Goto 0:◣Z=FL3 Lbl 1L4 ″ZH=″:Z=Z+N:Pause 0L5 ″H=″:H=G+(Z-F)I:Pause 2L6 Goto 1使用说明屏幕提示执行操作屏幕显示计算结果ZH0? 输入已知桩号ZH= 计算点的桩号GC0? 输入已知桩号的高程I? 输入纵坡(前进方向下坡为负，上坡为正)H= 计算点的高程N? 输入桩号增量【余弦定理的坐标运用】文件名：YXDL-ZBL1 A″AX″:B″AY″:C″BX″:D″BY″:E″CX″:F″CY″L2 G″AB″=√((C-A)2+(D-B)2)◢H″BC″=√((E-C)2+(F-D)2)◢I″CA″=√((E-A)2+(F-B)2)◢J″<A″=cos-1 ((G2+I2-H2)÷2GI)◢J″360°-<A″=360-J◢J″<B″=cos-1 ((G2+H2-I2)÷2GH)◢J″360°-<B″=360-J◢J″<C″=cos-1 ((H2+I2-G2)÷2HI)◢J″360°-<C″=360-J使用说明屏幕提示执行操作屏幕显示计算结果AX? 输入A点坐标X AB= AB距离AY? 输入A点坐标Y BC= BC距离BX? 输入B点坐标X CA= CA距离BY? 输入B点坐标Y <A= <A角度CX? 输入C点坐标X <B= <B角度CY? 输入C点坐标Y <C= <C角度【施工坐标与大地坐标相互转换程序】文件名：ZBZH（调用子程序文件名：Z.FS，Z.ZH）L1 Lbl 0L2 {XYFG}:X″AX1″:Y″AY1″:U=F″OX1″:Z=G″OY1″L3 Prog ″Z.FS″L4 D=0=>Goto 0:◣E=AL5 Lbl 1L6 {XYHI}:X″AX2″:Y″AY2″:U=H″OX2″:Z=I″OY2″L7 Prog ″Z.FS″L8 D=0=>Goto 1:◣E≥A=>E″ANGLE1->2″=E-A◢≠=>E″ANGLE1->2″=E-A+360◢◣Prog ″Z.ZH″使用说明屏幕提示执行操作屏幕显示计算结果AX1? 输入A点在坐标系1中的坐标XANGLE1->2= 坐标系1到坐标系2的转角AY1? 输入A点在坐标系1中的坐标Y OX1? 输入O点在坐标系1中的坐标X OY1? 输入O点在坐标系1中的坐标Y AX2? 输入A点在坐标系2中的坐标X AY2? 输入A点在坐标系2中的坐标YOX2? 输入O点在坐标系2中的坐标XX2= 计算点在坐标系2中的坐标XOY2? 输入O点在坐标系2中的坐标YX1? 输入计算点在坐标系1中的坐标XY2= 计算点在坐标系2中的坐标YY1? 输入计算点在坐标系1中的坐标Y注意fx-4500P型计算器将L3、L7行改成：Prog Z.FSL8行最后一句改成：Prog Z.ZH【坐标反算子程序】文件名：Z.FSL1 B=X-U:C=Y-ZL2 D=√(B2+C2)L3 D=0=>Goto 0:◣C≥0=>A=cos-1(B÷D):≠=>A=360-cos-1(B÷D):◣Lbl 0【坐标转换子程序】文件名：Z.ZH（可单独执行，调用子程序：Z.FS）L1 U=F″OX1″:Z=G″0Y1″:H″0X2″:I″OY2″:E″ANGLE1->2″L2 Lbl 0L3 {XY}:X″X1″:Y″Y1″L4 Prog ″Z.FS″L5 D=0=>Goto 0:◣A=A-EL6 ″X2=″:X=H+Dcos A:Pause 0L7 Y″Y2″=I+Dsin A◢Goto 0使用说明屏幕提示执行操作屏幕显示计算结果OX1? 输入O点在坐标系1中的坐标XX2=计算点在坐标系2中的坐标XOY1? 输入O点在坐标系1中的坐标YOX2? 输入O点在坐标系2中的坐标XY2=计算点在坐标系2中的坐标YOY2? 输入O点在坐标系2中的坐标YANGLE1->2? 输入坐标系1到坐标系2的转角X1? 输入计算点在坐标系1中的坐标XY1? 输入计算点在坐标系1中的坐标Y注意fx-4500P型计算器将L4行改成：Prog Z.FS【仁宗海交通桥拱圈悬链线坐标计算程序】文件名：XLXL1 M=2.514:F=12.195:L=30.484L2 Lbl 0L3 {X}:J=30.484-XL4 P=-0.00011837194X2+0.0072169X◢A=Jln (M+√(M2-1))÷LL5 D=tan-1(FAsinh A÷J(M-1))◢Y″YZ″=F-F(cosh A-1)÷(M-1)+2363.821+P◢Z″YS″=Y+0.7÷cosD◢{H}:U=Z-H◢W″YX″=Y-0.7÷cosD◢Q=W-H◢Goto 0使用说明屏幕提示执行操作屏幕显示计算结果X? 输入拱圈X坐标(以拱圈左下角为坐标原点)P= 计算点的设计预拱度YZ= 中轴设计高程（已加了设计预拱度P）YS= 上缘设计高程（已加了设计预拱度P）H? 输入实测高程U=上缘设计高程与实测高程的高差：负数表明实测高程高于设计高程正数表明实测高程低于设计高程YX= 下缘设计高程（已加了设计预拱度P）Q=下缘设计高程与实测高程的高差：负数表明实测高程高于设计高程正数表明实测高程低于设计高程注意fx-4500P型计算器将L5行中sinh A改成(e A-e-A)÷2；将coshA改成(e A+e-A)÷2【高斯坐标转换程序】文件名：GSZHL1 A″H0″:Lbl 1:{BCH}:E=B″H1″+C″H2″:F″HD0″=(1+(A-.5E)÷6378245)H″HD″◢Goto 1【高斯坐标转换程序1】文件名：GSZH1L1 A″H0″:R=6378245:Lbl 1:{BCH}:E=B″H1″+C″H2″:F″HD0″=(R+A)H″HD″÷(R+.5E)◢Goto 1使用说明屏幕提示执行操作屏幕显示计算结果H0? 输入投影面高程HD0=高程投影后的平距H1?、H2? 输入该边两端点的高程HD? 输入平距【40米T型梁梁底模反拱度高程控制程序】文件名：FGD-Q（使用全站仪）L1 Lbl 0L2 {XQ}L3 P=.003X-.04+(20-X)2÷10000◢Z=P+2717.90358◢B=Z-Q◢Goto 0使用说明屏幕提示执行操作屏幕显示计算结果X? 输入X坐标（以0＃台梁端为(0,0)点，输入范围0～40）P= 计算点与0,0点的设计高差Z= 计算点的设计高程Q? 输入实测高程B= 设计高程－实测高程（正值表示该点偏低，负值表示该点偏高）注意fx-4500P型计算器将L3行中“÷”改为“/”，0＃台梁端高程为2717.90358【40米T型梁梁底模反拱度高程控制程序】文件名：FGD-S（使用水准仪）L1 Lbl 0L2 {XQ}L3 P=.003X-.04+(20-X)2÷10000◢Z=P+2717.90358◢B=Z-S+Q◢Goto 0使用说明屏幕提示执行操作屏幕显示计算结果X?输入X坐标（以0＃台梁端为(0,0)点，输入范围0～40）P= 计算点与0,0点的设计高差Z= 计算点的设计高程S?输入视线高程（单位：m，只输入一次）B=设计高程－实测高程（正值表示该点偏低，负值表示该点偏高）Q? 输入前视读数（单位：m）注意fx-4500P型计算器将L3行中“÷”改为“/”，0＃台梁端高程为2717.90358。

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-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≤下段直缓点桩号……△△△ (注：有多少弯道就加多少个△)。

O —X N Y N计算L1 ZX？输入置镜点X 坐标X “ ZX ” Y “ ZY ” ZKA “ ZH ” B ZY？输入置镜点丫坐标“ XZH ” S“ YZH ” O“ EJ”：R“ R Z？输入置镜点高程(L- , R+ ) ” = 9 — Goto E / I “ LS ”K？输入仪器高度L : Lbl E : Q= 9 : J “ HX ” U ZH？输入直缓点里程Goto 6 ： Lbl 9 XZH？输入直缓点X 坐标L2 YZH？输入直缓点丫坐标G : “ — E ” 丄Lbl 1 : {N} : N “ —L” ：EJ？输入直缓点方位角W=Abs ( A-AbsN ： N= 9 — Goto 9 R(L— , R +)？输入圆曲线半径,/ E=O: F= 9 : P= 9 : R 工9 — Q=90/ 左偏为负, 右偏为正n R : Goto2 / Rec ( W, 18 9 +O : LS？输入直缓点至缓圆点长度Goto4 / Lbl 2 L？输入完整曲线长度,即两倍L3 LS 加圆曲线长I < W — E=Q ( W-I : Goto 3 / HX？输入后视点X 坐标P=QW2/I Rec( W-Wx y5/9 9 R2I 2 , HY？输入后视点丫坐标P/3+O : Goto4 / Lbl 3 —E 显示后视方位角L4 —L？输入计算里程W > L-I — Goto 1 / F=2RsinE : E = 显示视方位角P=E+O+IQ FcosP+Rec ( I － Ix y3/9 S = 显示前视距离9 R2, IQ/3+O Lbl 4 M？输入1 进入①中桩计算模L5 式, 输入9 进入②边桩计算模式 , Q=B+Ans J=Q F=S+W+FsinP 输入竖直角度进入③测距仪计算U=Ans : N > 9 — Goto 6 / Lbl 5 模式L6{D} : D= 9 — Goto 1 / D=1 — Goto 8 ① 中桩计算模式/ H=P+E — 9 9 AbsD/D : J=Q+Rec XD = 显示计算里程中桩X 坐标( Abs D , P+E-9 9 AbsD/D YD = 显示计算里程中桩Y 坐标U=F+W Lbl 6 —L？输入下一计算里程L7 ② 边桩计算模式Pol ( J-X , U-Y : W v 9 — W+36 D？输入边桩距离,左侧输正值, 9 :工—W / Q=9 —G=Ans : Goto 右侧输负值9 / “ E=” 丄V: “ S=” 丄Lbl 7 E = 显示前视方位角L8 s= 显示前视距离{CMTH} : M= 9 — Goto 5 / M =1 M？输入1 显示边桩坐标—Goto 8 / C Abs sinM : “ — D= ”—L？输入下一计算里程丄V— Ans : “ / S=” 丄Z= 9 — Goto 7 ③测距仪计算模式/ T= 9 — Goto 7 / C？输入斜距K-T+H+Z+CcosM : “ —H ” 丄Goto —D = 显示平距5 / Lb1 8 / S= 显示距离仪器？mL9 说明“ 9 ” 为零 , O 为英文字J: “XD= ” 丄U: “ YD = ” 丄Goto 1C “ X 9 = ”：D “ Y 9 = ”：E “ W 9 = ”： N=Ans : R “ R= ”： L “ L 卩。

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

福建林业职业技术学院工程系土木工程室高国兴CASIO fx-4500PA计算中桩、边桩坐标程序直线、圆曲线中桩及左、右边桩计算程序fx-4800p=”：C“Xθ=”：D“Yθ=”：E“Wθ=”：N=Ans：R“R=”：L“Lμo=”：K“D=”：M“JμoQ=9θL÷（πR）：S=2RsinQ ：O =C＋Scos（N-Q）◢P=D+Ssin（N-Q）◢W=E-2Q◢X=O+Kcos（N-2Q-M）◢Y=P+Ksin（N-2Q-M）◢注：Xθ=起算点坐标Yθ=起算点坐标Wθ=起算点方位角R=曲线半径，线路前进方向，左+，右－直线R为无穷大=置镜点到待求点距离LμoD左、右边距离=前进方向左（+90），右（－90）JμoO=中桩坐标P=中桩坐标W=待求点方位角X=左、右边桩坐标Y=左、右边桩坐标缓和曲线中桩及左、右边桩计算程序Fx-4800pC“Xθ=”：D“Yθ=”：E“Wθ=”：N=Ans：R“R=”：F“Lθ=”：=”：K“D=”：M“Jμo=”←L“LμoU=L－L^5÷（40R2F2）：V=L^3÷（6RF）－L^7÷（336R^3F^3）：Q=tan－1（V÷U）：S=√（U2+V2）←B=9θL2÷（πRF）：O=C+Scos（N－Q）◢P=D+Ssin（N－Q）◢W=E－B◢X=O+Kcos（N－B－M）◢Y=P+Ksin（N－B－M）◢注：Xθ=起算点坐标Yθ=起算点坐标Wθ=起算点方位角R=曲线半径，线路前进方向，左+，右－Lθ=缓和曲线长=置镜点到待求点距离LμoD左、右待求距离=前进方向左（+90），右（－90）JμoO=中桩坐标P=中桩坐标W=待求点方位角X=左、右边桩坐标Y=左、右边桩坐标放样程序fx-4800pRLI“XZ”：W：“YZ”：G“XH”：H“YH”←B=H－W ←C=G－I ←Prog“P”：U=A←K=√（B2＋C2）←C=X“X－”－I←B=Y“Y－”－W←Prog“P”←L“L”=√（B2＋C2）◢V=A－U：V<θ==>V=V＋36θ⊿V “V=”◢PA=tan－1（B÷C）←B≥θ==> C≥θ==> A=A：≠=> C<θ==> A=A＋18θ⊿⊿≠=> B<θ==> C≥θ==> A=A＋36θ: ≠=>C<θ==>A= A＋18θ⊿⊿⊿←XZ、YZ—置镜点X、Y坐标；XH、YH—后视点坐标；X-、Y-—前视点坐标；输出：R=水平角、L=置镜点至前视点水平距离。

Fx—4800P计算器公路工程测量程序汇编蒋立日期：1998年12月《fx4800p》由坐标计算多边形面积程序程序清单程序名：ALbi0：{NAB}：N“N=”：S=0：I=1：I“N”=I◢A“X=”：B“Y=”：C=A：D=B ↙Lbi1：{XY}：I=I+1：I“N”=I◢X“X=”：Y“Y=”↙Lbi2：S=S+0.5 (AY-BX) ：A=X ：B=Y ↙Lbi3：I<N＝＞Goto1: ≠＞S=S+0.5(DX-CY) ↙Lbi4：S“S”=AbsS◢↙Lbi5：Goto0 ↙使用说明：按FILE ■Ａ→EXE ↙显示Ｎ＝？问多边形个数，如Ｎ＝３Ｎ＝１提示输入第１个角点的Ｘ、Ｙ坐标Ｘ＝？１↙Ｙ＝？１↙Ｎ＝２提示输入第２个角点的Ｘ、Ｙ坐标Ｘ＝？４↙Ｙ＝？５↙Ｎ＝３提示↙输入第３个角点的Ｘ、Ｙ坐标Ｘ＝？１↙Ｙ＝？５↙Ｓ＝６多边形（三角形）面积的等于６《fx4800p》计算直线、园曲线、缓和曲线上的中桩坐标、法线方向上边线坐标程序程序清单程序名：BLbi0 ：Defm12 ：{ABCDER}:A“JDK”：B“X0”：C“Y0”：D“PJ0”：E“FWJ：JD-ZH”：R“R”：H“Ls”↙Lbi1 ：F=28.6479*H/R：P=H2/（24R：Q=H/2-H∧3/（240R2：L=R（D-2F）*π/180+2H◢↙Lbi3 ：Z[12]=A：Z[11]=E：G=（R+P）*tan（D/2）+Q：V=Z[12]-G： W=V+L：V“ZH”=V◢W“HZ”=W◢↙Lbi4：Rec（G，Z[11]）：X=B+I：Y=C+J：X“ZHX”=X◢Y“ZHY”=Y◢Z[1]=X：Z[2]=Y↙Lbi5：{N，T}：N“LP：1，RP：0”：T“S：0，S”↙Lbi6:N=1＝＞Z=Z[11]-D+180：≠＞N=0＝＞Z=[11]+D-180≨≨↙Lbi7：Z＜0＝＞Z=Z+360：≨Rec（G，Z）：X=B+I：Y=C+J：X“HZX”=X◢Y“HZY”=Y◢Z[3]= X：Z[4]=Y ↙Lbi9：Z[11]-180＜0＝＞Z[11]=Z[11]+180：≠＞Z[11]=Z[11]-180≨↙LbiA：{M}：M“PK”↙LbiB：M≤L/2+V＝＞GotoC：≨GotoG ↙LbiC：M≤V＝＞Rec（V-M，Z[11]）：X=Z[1]-I◢ Y=Z[2]-J◢ GotoI：≠＞M≤H+V＝＞U=M-V：Prog“Z2”：≠＞O=F+180（M-V-H）/（πR）：Prog“Z3”≨↙LbiD：N=1＝＞S=Z[11]-tan-1（Y/X）：≠＞N=0＝＞S=Z[11]+tan-1（Y/X）≨≨↙LbiE：prog“Z4”：X=Z[1]+I◢Y=Z[2]+J◢↙LbiI：T=0＝＞GotoA：≨GotoN ↙LbiN：{T}：T“S”：M≤V＝＞Rec（T，Z[11]-90）：Prog“Z6”：GotoA：≠＞M ≤H+V＝＞G=28.6479*U2/（RH）：GotoO：≠＞GotoP ≨↙LbiO：N=1＝＞Z[12]=Z[11]-G-90：≠＞N=0＝＞Z[12]=Z[11]+G-90≨≨Z[12] ＜0＝＞Z[12]=Z[12]+360：≨Rec（T，Z[12]）：Prog“Z6”：GotoA ↙LbiP：Prog“Z7”：Prog“Z6”：GotoA ↙LbiG：N=1＝＞S=Z[11]-D：≠＞S=Z[11]+D≨M≥W＝＞Rec（M-W，S）：X=Z[3]+I◢Y=Z[4]+J◢GotoV：≠＞M≤W-H＝＞O=F+180（W-H-M）/（πR）：Prog“Z3”：≠＞U=W-M：Prog“Z2”≨↙LbiH：N=1＝＞S=Z[11]-D+180+tan-1（Y/X）：≠＞N=0＝＞S=Z[11]+D+180- tan-1（Y/X）≨≨↙LbiK：Prog“Z4”：X=Z[3]+I ◢Y=Z[4]+J◢↙LbiV：T=0＝＞GotoA：≨{T}：T“S”：M≥W＝＞Rec（T，S-90）：Prog“Z6”：GotoA：≠＞M≤W-H＝＞GotoS：≨GotoT ↙LbiS：Prog“Z7”：Prog“Z6”：GotoA ↙LbiT：G=28.6479*U2/（RH）↙LbiU：N=1＝＞Z[12]=Z[11]-D+G-90：≠＞N=0＝＞Z[12]=Z[11]+D-G-90≨≨Z[12]＜0＝＞Z[12]=Z[12]+360：≨Rec（T，Z[12]）：Prog“Z6”：GotoA ↙子程序程序名：Z2程序清单：X=U-U∧5/（40R2H2）+U∧9/（3456R∧4H∧4）-U∧13/（599040R∧6H∧6）：Y=U∧3/（6RH）-U∧7/（336R∧3H∧3）+U∧11/（422400 R∧5H∧5）↙程序名：Z3程序清单：X=RsinO+Q：Y=R（1-cosO）+P ↙程序名：Z4程序清单：Rec（√—（X2+Y2），S）↙程序名：Z6程序清单：Z[5]=X+I：Z[6]=Y+J：Z[7]=2X-Z[5]：Z[8]=2Y-Z[6]：Z[5]“XL”=Z[5]◢ Z[6]“YL”=Z[6]◢ Z[7]“XR”=Z[7]◢ Z[8]“YR”=Z[8]◢↙程序名：Z7程序清单：N=1＝＞Z[12]=Z[11]-90-D/2：≠＞N=0＝＞Z[12]=Z[11]+90+D/2≨≨Z[12] ＜0＝＞Z[12]=Z[12]+360：≨Rec（（P+R）/cos（D/2），Z[12]）：Z[9]=B+I：Z[10]=C+J：N=1＝＞Pol（Z[9]-X，Z[10]-Y）：≠＞N=0＝＞Pol（X-Z[9]，Y-Z[10]）≨≨J ＜0＝＞J=J+360：≨Rec（T，J）↙B程序使用说明启动程序：prog“B”↙显示输入说明例如JDK ？交点桩号10035.89 ↙X0 ？交点X坐标26714.542 ↙Y0 ？交点Y坐标145904.039 ↙PJ°？交点偏角(度、分、秒) 49°27′32″↙FWJ°：JD—ZH ？交点至ZH点切线方位角204°29′01″↙R ？园曲线半径710 ↙Ls ？缓和曲线长度，（当Ls=0时，为园曲线）125 ↙﹎﹎﹎﹎﹎﹎﹎﹎计算结果显示﹎﹎﹎﹎﹎﹎﹎﹎L=737.886 曲线总长度↙ZH=9645.978 ZH点桩号↙HZ=10383.865 HZ点桩号↙I ？中间显示↙J ？中间显示↙ZHX=26359.691 ZH点的X坐标↙ZHY=145742.447 ZH点的Y坐标↙LP:1,RP:0 ？[LP—线路左偏为1,RP—线路右偏为0 ] 0 ↙S：0，S ？0—只计算线路中线坐标，S—不等于0的任何数，计算中线坐标和边坐标，边坐标是指对应中桩法线方向的坐标，10 ↙HZX=26822.392 HZ点的X坐标↙ZHY=146278.738 HZ点的Y坐标↙PK ？待求点的桩号9770.98 ↙X=26471.845 待求点中桩X坐标↙Y=145797.548 待求点中桩Y坐标↙S= ？线路中点至待求点法线方向的距离15 ↙XL=26479.238 线路左15m的X坐标↙YL=145784.496 线路左15m的Y坐标↙XR=26464.453 线路右15m的X坐标↙YR=145810.600 线路右15m的Y坐标↙PK ？另一个待求点的桩号《fx4800p》计算高速公路匝道中心线、边线坐标程序程序清单程序名：CLbi0 ：Defm12 ：{ABCDFGKRP}: B“BP”：C“XB”：D“YB”：F“XP”：G“YP”：K “PJ°”:R“R”：A“A”：P“R1”：Pol（C-F，D-G）：E=J+360 ↙Lbi1：{N，T，M}：N“LP：1，RP：0”：T“S：0，S”：M“L＝＞R：1，R＝＞L：0”：A=0＝＞GotoJ：≨Goto2 ↙Lbi2：U=A2/P：H= A2/R：P=9E99＝＞Z[3]=C：Z[4]=D：Goto3：≨U= A2/P：H= A2/R：Goto3 ↙Lbi3：M=1＝＞Z[1]=B-U：Z[2]=Z[1]+H：Z[1]“ZH”=Z[1]◢Z[2]“HY”=Z[2] ◢Goto4：≠＞M=0＝＞Z[1]=B+U：Z[2]=Z[1]-H：Z[1]“HZ”=Z[1]◢Z[2]“YH”=Z[2] ◢≨≨↙Lbi4:P=9E99＝＞Goto5：≨Prog“Z2”：W=28.6479*U2/（RH：S=X-Y/tanW：Rec（Y/sinW，E）：X=C+I：Y=D+J：N=1＝＞E=E+W：≠＞N=0＝＞E-W≨≨Rec （S，E）：Z[3]=X+I：Z[4]=Y+J ↙Lbi5：{Q}：Q“PK”↙Lbi6：M=1＝＞Goto7：≠＞M=0＝＞Goto8≨≨↙Lbi7：Q≦Z[2] ＝＞Goto9：≠＞GotoD≨↙Lbi8:Q≧Z[2] ＝＞Goto9：≠＞GotoD≨↙Lbi9：M=1＝＞U=Q-Z[1]：≠＞M=0＝＞U=Z[1]-Q≨≨Prog“Z2”↙LbiA：N=1＝＞V=E+180-tan-1（Y/X）：≠＞N=0＝＞V= E+180+tan-1（Y/X）≨≨Rec （√—（X2+Y2），V）：X=Z[3]+I：Y=Z[4]+J：X“XP”=X ◢Y“YP”=Y ◢↙LbiB：T=0＝＞Goto5：≨GotoC↙LbiC：{T}：T“S”：W=28.6479*U2/（RH：N=1＝＞W=E+180-W-90：≠＞N=0＝＞W=E+180+W-90≨≨Rec（T，W）：Prog“Z5”：Goto5↙LbiD：U=H：Prog“Z2”：W=28.6479*H/R：N=1＝＞O=E+180-tan-1（Y/X）：≠＞N=0＝＞O=E+180+tan-1（Y/X）≨≨Rec（√—（X2+Y2），O）：Z[9] =Z[3]+I：Z[10]=Z[4]+J ↙LbiE：N=1＝＞V=E+180-W-90：≠＞N=0＝＞V=E+180+W+90≨≨Rec（R，V）：Z[11]=Z[9]+I：Z[12]=Z[10]+J ↙LbiF：M=1＝＞L=Q-Z[2]：：≠＞M=0＝＞L=Z[2]-Q≨≨Prog“Z8”↙LbiH：{T}：T“S”：T=0＝＞Goto5：≨Prog“Z9”：Goto5↙LbiJ：Pol（F-C，G-D）：W=J+360 ↙LbiK：{Q}：Q“PK”↙LbiM：N=1＝＞V=W-90：≠＞N=0＝＞V=W+90≨≨Rec（R，V）：Z[11]=C+I：Z[12]=D+J ↙LbiN：L=Abs（Q-B）：Prog“Z8”↙LbiO：{T}：T“S”：T=0＝＞GotoK：≨Prog“Z9”：GotoK↙子程序程序名：Z5Lbi0：Z[5]=X+I：Z[6]=Y+J：Z[7]=2X-Z[5]：Z[8]=2Y-Z[6] ↙Lbi1：M=1＝＞Z[5]“XL”=Z[5]◢Z[6]“YL”=Z[6]◢Z[7]“XR”=Z[7]◢Z[8]“YR”=Z[8]◢≠＞M=0＝＞Z[7]“XL”=Z[7]◢Z[8]“YL”=Z[8]◢Z[5]“XR”=Z[5]◢Z[6]“YR”=Z[6]◢≨≨↙程序名：Z8Lbi0：N=1＝＞Z=V+180-L*180/（πR）：≠＞N=0＝＞Z= V+180+L*180/（πR）（R，Z）：X=Z[11] +I：Y=Z[12] +J ：X“XP“=X ◢Y“YP“=Y ◢↙≨≨Rec程序名：Z9N=1＝＞Z=Z-180：≠＞N=0＝＞Z=Z≨≨Rec（T，Z）：Prog“Z5”↙C匝道程序使用说明FUNCTION 6 5 （变量清零）↙Prog“C”显示输入例如BP ？起点桩号754.493 ↙XB ？起点X坐标26541.168 ↙YB ？起点Y坐标145907.887 ↙XP ？第一段曲线X坐标26552.825 ↙YP ？第一段曲线Y坐标145935.665 ↙PJ°？园曲线偏角118.41°↙R ？园曲线偏半径80 ↙A？缓和曲线参数(园曲线A=0) 60 ↙R1 ？缓和曲线起点半径(当R无穷大时,R1=9EXP99) 9EXP99 ↙LP:1,RP:0 ？(左偏为1,右偏为0) 1 ↙S：0，S ？{ 计算边坐标为S（S为不为零的数），不计算为0} 10 ↙L＝＞R：1，R＝＞L：0？（计算顺序，从左往右为1，从右往左为0）0 ↙﹎﹎﹎﹎﹎﹎﹎﹎计算结果显示﹎﹎﹎﹎﹎﹎﹎﹎HZ=754.493 缓直点桩号↙(ZH) (直缓点桩号)YH=709.493 园缓点桩号↙PK ？输入待求点桩号544.161 ↙I ？中间显示↙XP=26698.387 待求点X坐标↙YP=145928.101 待求点Y坐标↙S ？从中点至法线方向的距离10 ↙XL=26707.612 线路前进方向左侧X坐标↙YL=145931.962 线路前进方向左侧Y坐标↙XR=26689.162 线路前进方向右侧X坐标↙YL=145924.241 线路前进方向右侧Y坐标↙PK ？另一个待求点的桩号《fx4800p》纵曲线计算程序（包括直线、曲线、直线）程序名：D程序清单：Lbi0：{ABCDERNT}：A“PVI”：B“HI”：R“R”：C“I1”：D“I2”↙Lbi1：N“LP：1，RP：0”↙Lbi2：N=1＝＞T=R（D-C）/200：≠＞N=0＝＞T=R（C-D）/200≨≨E=A-T：F=A+T：E“PVC”=E◢F“PVT”=F◢↙Lbi3：{G}：G“PK”：G≦A＝＞Goto4：≨Goto5 ↙Lbi4：G≦E＝＞H=B-C（A-G）/100：H“HP”=H◢Goto3：≠＞Prog“H1”：N=1＝＞H=B+I-J：H“HP”=H◢Goto3：≠＞N=0＝＞H=B-I-J：H“HP”=H◢Goto3≨≨↙Lbi5：G≧F＝＞H=B+D（G-A）/100：H“HP”=H◢Goto3：≠＞Prog“H2”：N=1＝＞H=B+I+J：H“HP”=H◢Goto3：≠＞N=0＝＞H=B-I+J：H“HP”=H◢Goto3≨≨↙子程序程序名：H1I=（G-E）2/（2R：J=C（A-G）/100 ↙程序名：H2I=（F-G）2/（2R：J=D（G-A）/100 ↙D程序使用说明Prog“D”显示输入说明例如PVI ？变坡点桩号1040 ↙HI ？变坡点高程36.778 ↙R ？园曲线半径25285.340 ↙I1 ？变坡点左侧坡度值(百分数) 1.018 ↙I2 ？变坡点右侧坡度值(百分数) -1.260 ↙LP:1,RP:0 ？逆时针左偏为1, 顺时针右偏为0 0 ↙﹎﹎﹎﹎﹎﹎﹎﹎计算结果显示﹎﹎﹎﹎﹎﹎﹎﹎PVC=752.00 纵曲线起点桩号↙PVT=1328.00 纵曲线终点桩号↙PK= ？待求点桩号↙HP=35.165 待求点P点的设计高程↙PK= ？输入待求点桩号↙。

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

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

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

如数据库中没有本数字所代表的点名，则系统显示SynERROR 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 CIRCLE”。

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

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

4800P计算器计算任意曲线点位坐标程序A“DK0”B“DK1”R“R0”V“R1”C“X0”D“Y0”E“W0”：N：LbI0：{LFJ}：L“DK”F“BL”J“JU”：T=（1÷R－1÷V）÷(B－A):H=(L－A)÷N:P=TNH－1÷R:S=90H÷π:W=E+(P－1÷R)NS▲G=1:X=C+H÷6×(cosE+cosW+4×∑(cos(E+((G+0.5)HT －2÷R)(G+0.5)S,G,0,(N－1))+2×∑(cos(E+（THG－2÷R)GS,G,1,(N－1)))+Fcos(W－J)▲Y=D+H÷6×(sinE+sinW+4×∑(sin(E+((G+0.5)HT－2÷R)(G+0.5)S,G,0,(N－1))+2×∑(sin(E+（THG－2÷R)GS,G,1,(N－1)))+Fsin(W－J)▲goto0:△注：A“DK0”——曲线起点里程B“DK1”——曲线终点里程R“R0”——曲线起点半径，曲线左偏为正、右偏为负。

（如为直线，则为无穷大）V“R1”——曲线终点半径，曲线左偏为正、右偏为负。

（如为直线，则为无穷大）C“X0”D“Y0”E“W0”——起点坐标、切线方位角N——求和分段数，一般可取2，如为了保证较高精度，则N可取3、4、5。

L“DK”——待求点处的里程F“BL”——距中桩距离J“JU”——横距偏角（与切线前进方向左偏为正、右偏为负）结果：W——待求点处的切线方位角、X——、Y——。

Program[EDIT]ULbI0：{ONZHVQ}：OVQ：Prog“Q”：N≤2→{XYSD}:SDXY:Prog“P”：W=J△LbI1：{KE}：K=0→Goto0△K＞H-L→X=T:Y=0:I=H-K:K＞H→AbsA＜180→I=I-T△Prog“L”：Goto3△Prog“H”：J=AbsA-J：C=AbSA：X=T-X：B=T：U=0：Prog“X”：Goto2△K＜Z+L→I=K-Z:K＜Z→X=I:Y=E:J=0:Goto3△Prog“H”：Goto2△I=K-Z-L：J=（90L+180I）/πR：A＞0→E=-E△I=R+E:X=M+IsinJ:Y=R+P-IcosJ：LbI2：A＜0→J=-J:Y=-Y△LbI3:O≠0→C=O:B=V:U=Q:Prog“X”△N=3→X▲Y▲Goto1△N=2→Prog“G”：Goto1Q{RLA}：M=L/2-L∧3/（240R2）▲P=L2/（24R）-L∧4/（2688R∧3）▲T=M+（R+P）tan（AbSA/2）▲F=πRAbsA/180+L▲XB=B+XcosC-YsinC：U=U+XsinC+YcosC：X=B：Y=ULPoI（-I，E）：C=A+J：Prog“Z”：J=APPoI(X-S,Y-D):Prog“T”HA＞0→E=-E△J=90I2/(πRL):X=I-I∧5/(40R2L2)+EsinJ:Y=I∧3/(6RL)-I∧7/(336R∧3L∧3)-Eco sJTJ≤0→J=360+J△J＞360→J=J-360GProg“P”：J=J-W：Prog“T”：J▲I▲ZX=X+Rec（I，C）：Y=Y+JXYLbI0：{SDXY}：SDXY：Prog“P”：W=J：{BA}：BA：E=S+Rec （B，W+A）：F=D+J：E▲F▲Goto0IJLbI0：{SDXY}：SDXY：Prog“P”：W=J：{EF}：EF：PoI（E-S，F-D）：Prog“T”：G=J：J=G-W：Prog“T”：J▲I▲Goto0为了分辩“零”和“O”，可变换字体U—曲线放样与坐标计算程序O—第一切线方位角V—ZH点X坐标Q—ZH点Y坐标L—缓和曲线长R—半径A—曲线偏角（左偏为负、右偏为正）N—N等于2时为曲线放样、等于3时为计算坐标N等于3时计算坐标K—计算点里程H—HZ点里程Z—ZH点里程E—偏距（左偏为负、右偏为正）结果：X=Y=N等于2时放样S—置镜点X坐标D—置镜点Y坐标X—后视点X坐标Y—后视点Y坐标J不管K—放样点里程H—HZ点里程Z—ZH点里程E—偏距（左偏为负、右偏为正）结果：J—顺拨角I—前视距IJ—放样程序S—置镜点X坐标D—置镜点Y坐标X—后视点X坐标Y—后视点Y坐标E—前视点X坐标F—前视点Y坐标结果：J—顺拨角I—前视距XY—计算前视点坐标程序S—置镜点X坐标D—置镜点Y坐标X—后视点X坐标Y—后视点Y坐标B—前视距A—顺拨角结果：E—前视点X坐标F—前视点Y坐标。

卡西欧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——导线全长闭合差、导线全长相对闭合差。