卡西欧4500、4800计算器计算程序

CASIO fx-4800P计算器程序DA DI——N E (公路弯道坐标计算)K“JD”：R：S“LS”：F“PJ——L－R＋”：A“ZH——N”：B“ZH——E”：C“JD——N”：D“JD——E”M“ZH”=K－（R＋S2÷（24R））tan（Abs F÷2）－S÷2＋S^3÷（2 40R2）▲N“HY”=M＋S▲O“QZ”=M＋（πRAbs F÷180＋S）÷2▲Z[1]“YH”=M＋πRAbs F÷180▲Q“HZ”= Z[1]＋S▲T=K－M▲L=Q－M▲E=（R＋S2÷（24R））÷cos（F÷2）－R▲Lb1 0：{G}：G“ZHUANG HAO”：G＜O＝=＞Z=G－M≠=＞Z=Q－G：△{U}：U“I——B OUT＋IN－”Z＜S＋0.005＝=＞H=tan-1（24Z2S2R2－Z^6）÷（48Z^3R^3－6Z^4 SR））：V=Z－Z^5÷（40S2R2）＋Usin H：W=Z^3÷（6SR）－Z ^ 7÷（336Z^3R^3）－Ucos H：≠=＞H=180（Z－S÷2）÷（πR）：V=Rsin H＋S÷2－S^3÷（240R2）＋Usin H：W=R－Rcos H＋S2÷（24R）－Ucos H：△F＜0＝=＞G＜O＝=＞X=V：Y=W：≠=＞X=T＋Tcos F－Vcos F－Ws in Abs F：Y=Tsin Abs F－Vsin Abs F＋Wcos F：△≠=＞G＜O＝=＞X=V：Y=－W：≠=＞X=T＋Tcos F－Vcos F－Wsin Abs F：Y=Vsi n Abs F－Tsin Abs F－Wcos F：△△C=A＝=＞D≥B＝=＞P=90：≠=＞P=270：△≠=＞C＜A＝=＞P= tan-1（（D－B）÷（C－A））＋180：≠=＞P= tan-1（（D－B）÷（C－A））：△△I“COM——N”=A＋√（X2＋Y2）cos（P－tan-1（Y÷X））▲J“COM——E”=B＋√（X2＋Y2）sin（P－tan-1（Y÷X））▲G=G＋20Goto 0程序运行输入：JD？——交点桩号R？——圆曲线半径LS？——缓和曲线长度PJ——L－R＋？——偏角右偏为正ZH——N？——ZH点N坐标ZH——E？——ZH点E坐标JD——N？——JD点N坐标JD——E？——JD点E坐标输出：ZH= ——ZH点桩号HY= ——HY点桩号QZ= ——QZ点桩号YH= ——YH点桩号HZ= ——HZ点桩号T= ——切线长L= ——曲线长E= ——外矢距再输入：ZHUANG HAO？——输入任一点桩号I——B OUT＋IN－？——距中线距离（外侧为正）再输出：COM——N= ——计算出N坐标COM——E= ——计算出E坐标以此循环——再输入——再输出程序中各字母所代表的含义：A——ZH点N坐标B——ZH点E坐标C——JD点N坐标D——JD点E坐标E——外矢距F——偏角（右偏为正）G——弯道内任一点桩号H——弯道内任一点切线角I——计算出的N坐标J——计算出的E坐标K——交点桩号L——曲线长M——ZH点桩号N——HY点桩号O——QZ点桩号P——ZH--JD方位角Q——HZ点桩号R——圆曲线半径S——缓和曲线长度T——切线长U——距中线距离（外侧为正）V——支距W——支距X——弯道内坐标Y——弯道内坐标Z——任一点曲线长度Z[1]——YH点桩号。

LX、RX：左侧或右侧点X坐标
LY、RY：左侧或右侧点Y坐标
I：置镜点至左或右点距离
J：置镜点至左或右点方位角
CONT：输入0计算下一里程；输入1则从本程序第一步开始
QX1计算范围：
DK214+832.674~DK216+053.229
QX2计算范围：
DK215+707.196~DK217+253.771(反算)
LbI 5 ↙
H=1=〉Y= -1×Y△ U=B+X cosT-YsinT：V=C+XsinT+YcosT：↙
X=U▲
Y=V▲
PoL（（X-M，Y-N）↙
I ▲
J＜0=〉J=J+360△↙
J ▲
{W}：W“NEXT”：W=1=〉Goto 1△↙
LHale Waihona Puke I 6 ↙ {O}：O“O=”：O=1=〉Goto 7△↙
QX3计算范围：
DK216+623.229~DK219+480.262
QX4计算范围：
DK217+923.771~DK222+114.707(反算)
J ▲
F=T ▲
Goto 6 △↙
LbI 3 ↙
L＞A=〉Goto 4 △ X=L-L5÷40÷（RA）2：Y=L3÷（6RA）-L7÷336÷（RA）3：Z=60×L2÷（πRA）+tan-1（Y÷X）：H=1=〉Z= -1×Z △F=T+Z▲
Goto 5 △↙
LbI 4 ↙
Z=180（Abs（E-D）-A）÷（πR）+90A÷（πR）：X=R sinZ+A÷2-A3÷（240R2）：Y=R（1-cos Z）+A2÷（24R）-A4÷（2688R3）：H=1=〉Z= -1×Z △ F=T+Z▲

竖曲线高程的计算公式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↙………依次类推，计算原始数据完成输入，坡度换算成小数。

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 4800计算器程序使用说明一、程序列表：1。

PQX：主程序，计算路线的中边桩坐标。

(1)"L0"--- 输入测段的近似桩号，来调取合适的平曲线要素数据，输入0则默认上次输的平曲线要素。

输入-1可以手工输入交点桩号JL，交点X坐标JX，交点Y坐标JY，I0为前一交点与本交点的直线方位角，J0为本交点的偏角(左负右正)，R为圆曲线半径，LS为缓和曲线长度，L1为上一交点的HZ(或YZ)的桩号，L2为下一交点的ZH(或ZY)桩号。

(2)"L"---输入计算点的桩号(3)"JJ"---输入与路线前进方向的右交角（正交为90）(4)"YC"---输入与中桩的距离(左侧为负,右侧为正)注:若JJ和YC中任一个输入为0,则可返回上一级,重新输入桩号计算.若YC输入为-1,可以输入实测点的坐标X和Y来反算与中桩宽度YC,桩号差值DL和该点的近似桩号L.(5)"X,Y"数值1数值2显示计算点的坐标值X,Y(6)"A,D"ddd°mm′ss″数值4显示测站点与放样点的方位角及平距(7)"TMP"---进入自由运算状态,若计算结果为0则退出该状态.该状态会破坏变量I的内容.(8)"YC"---输入-1进入反算桩号模式,程序要求输x,y坐标,若x,y坐标其中任一输入为0则退出该状态.2.LYC程序:两个作用:一是根据实测坐标反算桩号L和与中桩宽度YC, 二是根据输入的三维坐标x,y,z和输入的设计标高SG以及"YC0","HC0","M"来进行边坡放样.(1)"L0"--- 输入测段的近似桩号，(2)"XF","YF","ZF"----输入测点的三维坐标x,y,z若ZF输为0则只反算桩号,若不为0则进行边桩放样.(3)"DL"---反算桩号时显示中间结果桩号差值DL,若DL值过大,可以重新输入测段的近似桩号来计算.(4) "L,YC"---显示桩号L和与中桩距YC(左负右正),存在变量L和O中.(5) "SG:"---输入路基边缘的设计高程(请分清左右),存在变量Z[10]中.(6) TW---显示填挖高度值(负值为挖方高度,正值为填方高度),存在变量C中.(7) "YC0"-路基横断面图上最后一个变坡点与中桩的水平宽度(左右均为正值),存在变量P中.(8)"HC0"---路基横断面图上最后一个变坡点与路基边缘设计标高的高差,挖方平台为正值,填方平台为负值,存在变量Q中.(9) "M="---路基横断面图上最后一级边坡的坡比,存在变量T中.(10) "BL"---显示在实测标高位置的路基设计宽度,存在变量U中.(11)"DB"---显示设计宽度和实测宽度的差值,存在变量J中.若DB为负值,则应往中桩位置移动(宽度减小),若DB为正值,则应往路线外侧移动(宽度增大).3、ZZ子程序---计算中桩坐标以及前进方位角的子程序。

CASIO 4500全站仪测设公路中桩和边桩程序作者：lq52搜集人气：22 来源：网络【字体：大中小】【打印】【关闭】用全站仪测设公路中桩和边桩程序[转帖]由本人从网络是整理而得。

供大家使用！！(原程序是这样的，我未加整理，见谅)程序说明：本程序用于用全站仪测设中桩和边桩，采用坐标变换的方法编制而成，在使用fx-4500计算器时，可直接输入以下程序便可使用；而在使用fx-4800计算器时，见fx-4800程序清单。

注意：因程序的0(零)和O易混淆,所以字母O用较大的全角符号，零仍用标准字。

程序输入时应连续输入，无空格。

而且输入时文件名需和给出的一样，否则会出现调用文件错误。

变量说明： <程序运算符定义>X 输入摆站点的x坐标 C 输入第一个曲线角桩x坐标 I 输入中间曲线角桩的x坐标Y 输入摆站点的y坐标 D 输入第一个曲线角桩y坐标 J 输入中间曲线角桩的y坐标E 最后一个曲线角桩x坐标 R 中间曲线半径 L 中间曲线缓和曲线长F 最后一个曲线角桩y坐标 K 第一个K为中间角桩桩号 S 中间角桩曲线总长，检验用K 第二个K为要测设的桩号 B 测中桩时输入0，测边桩时输入宽度（左输负值，右输正值）程序清单注释文件:F1 文件名：1L1L2L3Defm 8XYCDIJEFRLK:Fixm:Z[2]=0:Prog BZ[5]=0.5L-L Xy 3/240R2:Z[3]=Z[5]+(R+L2/24R)Abs tan T┛2:Z=K-Ans:S=R Abs Tπ/180+L:S◢H=Z+Ans:G=RL:Prog 2文件:F2 文件名：2L1Lbl 0:Fixm:{K}L2Lbl 6:{B}:Norm:K(-1) Xy Z[2]>(H-L)(-1) Xy Z[2]＝>Goto 4⊿Frac Z[2]┛2=0=>Z[6]=Z:≠>Z[6]=H⊿A=90T/Abs TL3B<0＝>N=0:≠>N=1⊿B=Abs BL4０=Abs(K-Z[6]:K(-1) Xy Z[2]>Z[6](-1) Xy Z[2]＝>Goto 2⊿Q=Z[3]+０:Z[4]=B L5B≠0＝>Prog C⊿Prog D:P=0L6Prog E:Goto 0L7Lbl 2:０≥L＝>Goto 3⊿L8Q=Z[3]-０+０ Xy 5/40G2:Z[4]=０ Xy 3/6G-０ Xy 7/336G Xy 3:Prog D:P=BL9B≠0＝>Prog F⊿Prog E:Goto 0L10Lbl 3:Q=Z[3]-Z[5]:Z[4]=R+L2/24R:Prog DL11Z[8]=U+A(2(０-L)/πR-1+L/πR):B≠0＝>Prog G:≠>P=R⊿Prog E:Goto 0L12Lbl 4:M=C:０=D:C=E:D=F:E=M:F=０:Z[2]=Z[2]+1:Prog B:Goto 6文件:F3 文件名：0L1Pol(M,０):W<0＝>W=W+360文件:F4 文件名：BL1L2M=E-I:０=F-J:Prog 0:T=W:M=I-C:０=J-D:Prog 0:U=W:T=T-WT<-180＝>T=T+360⊿T>180＝>T=T-360文件:F5 文件名：CL1 A=-90(-1) Xy (N+Z[2])文件:F6 文件名：DL1 Z[7]=I-QcosU+Z[4]cos(U+A):Z[1]=J-QsinU+Z[4]sin(U+A)文件:F7 文件名：EL1L2M=Z[7]+PcosZ[8]-X:０=Z[1]+PsinZ[8]-Y:Prog 0W◢Fix 3:V◢Norm文件:F8 文件名：FL1Z[8]=U+A+A０2/πG+90+A(-1) Xy (N+Z[2])文件:F9 文件名：GL1P=R+AB(-1) Xy (N+Z[2])/90此处L2句加入： M+X◢０+Y◢可显示测点坐标值曲线上任意点中边桩坐标计算程序(CASIO fx-4800P可全线计算)作者：冯第胜人气：14 来源：网络【字体：大中小】【打印】【关闭】主:PQXLbi 0：prog“TAD”Goto 1↙Lbi 1:Z[1]“ZH-X”=M+Tcos（F+180）:Z[2]“ZH-Y”=N+Tsin（F+180）：V[1]“HZ-X”=M+Tcos（F+Z）：V[2]“HZ-Y”=N+Tsin（F+Z）：H=0=>Goto 2:≠>D=H:prog“QXY”:A“HY-X”= Z[1]+I cos(F+SZ÷AbsZ) ：C“HY-Y”= Z[2]+I sin(F+SZ÷AbsZ):Q“YH-X”=V[1]+I cos(F+Z+180-SZ÷AbsZ ) ：E“YH-Y”= V[2]+ Isin(F+Z+180-SZ÷AbsZ )↙Lbi 2:K≤G=>D=Abs(G-K):O=M:P=N:I=D+T:J=F+180:W=F:Goto 4⊿↙H=0→Goto 3⊿↙K≤G+H=>D=Abs(K-G):O=Z[1]:P=Z[2]:prog“QXY” :J=F+SZ÷AbsZ:W=F+3SZ÷AbsZ: Goto 4⊿↙Lbi 3:K≤G+L-H=>D=Abs(K-G-H):H=0=>H=1:O=Z[1]:P=Z[2]:≠>H=H:O=A:P=C:prog“QXY”H=1=>Q=0⊿I=D:J=F+(Q+E)Z÷AbsZ:W=F+(Q+2E)Z÷AbsZ: Goto 4⊿↙K≤V=>D=Abs(V-K):O=V[1]:P=V[2]:prog“QXY”:F=F+Z+180: :J=F-SZ÷AbsZ:W=F-3SZ÷AbsZ:U=180+U: Goto 4⊿↙Lbi 4:X=O+IcosJ+Bcos(W+U)▲Y=P+IsinJ+Bsin(W+U)▲I=0：J=0：I=POL（X-x0，Y-y0：J<0→J=J+360△J“FWJ”=J▲I“r”=I▲（x0、y0为置仪点坐标）子:QXYA=D-D^5/40R²H²+ D^9/3456R^4H^4- D^13/599040R^6H^6+ D^17/17542600R^8H^8+1E-99:C=D³/6RH-D^7/336R³H³+D^11/42240R^5H^5-D^15/976800R^7H^7+ D^19/3530097000R^9H^9: S=tan¯1(C÷A):I=√(A²+C²):Q=90H÷лR:E=90D÷лR:D=2RsinE子:TADLbi 0:{KBU}:KBU↙K≤xxx=>M=xxx:N=xxx: F=xxx: Z=xxx: R=xxx: H=xxx: T=xxx: L=xxx: G=xxx: V=xxx:≠>K≤xxx=>M=xxx:N=xxx: F=xxx: Z=xxx: R=xxx: H=xxx: T=xxx: L=xxx: G=xxx: V=xxx:≠>K≤xxx=>M=xxx……………⊿注:K？输入计算点桩号;B？输入距中桩距离,左－,右＋;U？输入与前进方向偏角计算中桩坐标时B、U输入0,M－交点X坐标;N－交点Y坐标;F－路线前进方向方位角;Z－转角,左转时连同“-”号输入; R－半径;H-缓和曲线长,无缓和曲线时为0;T－切线长；L－曲线总长;G－直缓点桩号;V－缓直的桩号在Lbi1行处,把各主点后的“:”改为“▲”，如Z[1]“ZH-X”=M+Tcos（F+180）:改成Z[1]“ZH-X”=M+Tcos（F+180）▲可直接显示主点桩号坐标,在子程“TAD”内可输入全线角桩资料,K≤本转角缓直点桩号时输入本转角曲线要素。

QL（坐标反算）Lbi 0:C“X0”:D“Y0”:{XY}:X“X2”Y“Y2”:(X-C)=O=>Goto 1:≠Lbi 1:W“Q”=0◢Goto 0Lbi 2:Prog“FWJ”Fix 4W“Q”=IntW+0.01Int(60FracW)+0.006Frac(60FracW) ◢NormT“L”= ((X-C)2+(Y-D)2)◢Goto 0说明：（已知两点坐标，反算方位角及两点间的距离）先按FILE键，光标选择QL，按EXE运行。

输入：X0? 测站点（或起算点）X轴坐标，按EXEY0? 测站点（或起算点）Y轴坐标，按EXEX2? 后视点（或方向点）X轴坐标，按EXEY2? 后视点（或方向点）Y轴坐标，按EXE显示： Q 方位角（如：125.0325即表示125。

03’25”）,按EXE L 两点间距离再按EXE进行循环运算，重复输入X2？、Y2?否则按AC键退出。

XY(坐标正算)Lbi 0:C“X0”:D“Y0”:{WM}:W“Q”M“L”:X=C+McosW◢Y=D+MsinW◢Goto 0说明：（已知一点坐标、方位角、距离，求另一点坐标）先按FILE键，光标选择XY，按EXE运行。

输入：X0? 测站点（或起算点）X轴坐标，按EXE Y0? 测站点（或起算点）Y轴坐标，按EXE Q? 已知方位角，按EXEL? 两点间距离，按EXE显示：X 所求点X轴坐标，按EXEY 所求点Y轴坐标，按EXE再按EXE进行循环运算，重复输入Q? 、L?否则按AC键退出。

ZX（直线坐标计算）Lbi 0:G“QD”：A“QDX”: B“QDY”:Q“Q”:{L,K}:Prog“ZZ”:Goto 0说明：（计算直线段任意里程桩号坐标）先按FILE键，光标选择XY，按EXE运行。

输入：QD? 直线段起点里程，按EXEQDX? 起点X轴坐标，按EXEQDY? 起点Y轴坐标，按EXEQ? 已知方位角，按EXEL? 所求点里程，按EXEK? 所求点距中线的宽度（左负右正），按EXE显示：X 所求点X轴坐标，按EXEY 所求点Y轴坐标，再按EXE进行循环运算，重复输入L? 、K?否则按AC键退出。

/
N
+
计算曲线长。G、H加E-21的目的是防止“0”溢出，增强程序的通用性。
(
O
+
U
)
/
2
:
″
Q
X
C
=
″
◢
G
=
R
O
+
EXP
(-)
2
1
:
H
=
R
U
+
EXP
(-)
2
1
L3
V(I)
=
O
:
W(J)
=
L
-
O
-
U
:
K
=
U
调用子程序Z2，计算曲线坐标。
:
Prog
Z
2
L4
Y
=
Y
Z
/
Sin
P
:
X
-
Y
Cos
P
:
″
计算第一切线长L1、第二切线长L2。
(
O
+
W(J)
)
/
2
R
L4
E
=
K
-
(
U
xy
5
-
(
U
-
K
)
xy
求第二缓和曲线坐标增量计算元素E和F。
5
)
/
4
0
H2
:
F
=
(
U
xy
3
-
(
U
-
K
)
xy
3
)
/
6
H
-
(
U
xy

卡西欧4500、4800计算器计算程序卡西欧4500/4800计算器程序任意里程桩号计算程序1.主程序(TYQXjs)"1.SZ => XY"："2.XY => SZ"：N：U"X0"：V"Y0"：O"S0"：G"F0"：H"LS"：P"R0"：R"RN"：Q：C =1÷P：D=(P-R)÷(2HPR)：E=180÷π：N=1=>Goto 1：≠>Goto 2Δ←┘Lbl 1：{SZ}：SZ：W=Abs(S-O)：Prog "SUB1"：X"XS"=X◢Y"YS"=Y◢F”FS”=F-90◢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))+Bc os(G+QEFW(C+FWD))+Acos(G+QE M W(C+M WD))):Y=V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD) )+Bsin(G+QEFW(C+FWD)) +Asin(G+QE M W(C+M WD))) :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。

测量实用计算程序（适用于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 为右侧坐标注：针对具体一个曲线可将已知条件数据直接输入程序中。

CASIO4800系列计算器线路线路正反算程序经本人苦心钻研，奋战多日，终于编写出了速度快，精度高，功能全的线路坐标正反算程序，欢迎试用并提出宝贵意见。

部分引用请注明出处。

功能简介及特点：1、选用高斯-勒让德公式作计算内核，保证精度，模块化设计，便于扩充功能。

2、线元数据可自动从数据库调用，也可手工输入。

3、可管理多条线路,如里程不在线路或线元范围，将警告里程偏大、偏小。

4、边桩计算设计为导线式递推方式，可用于由一个中桩推出结构物所有角点坐标。

5、反算实现了智能化操作，只需输入线路号（或手工输线元资料）、坐标，不需近似里程，即可自动从起点向后开始试算出里程、位置，如对算出里程、位置表示怀疑，还可以让计算器从终点起再向前试算下一个可能的位置（匝道、回头曲线同一坐标可能会有一个以上结果）。

6、程序代码规范简洁，便于阅读、理解。

程序清单:ZFS///正反算主程序Norm:Lbl 1:U"1 ZS 2 FS"=1=>Prog "ZS":≠>U=2=>Prog"FS":≠>Goto 1ZS///正算子程序{T}:Lbl 1:{K}:T"1 B"=1=>Prog"ZZ":Prog"WY":Goto 1△Prog"ZZ":Prog"ZB":Goto 1FS///反算子程序Lbl 1:{VW}:V"X"W"Y":U=2:Lbl 2:Prog "ZZ":Pol(V-X,W-Y:J=J-O:I=Rec(I,J:K=K+I:AbsI<0.001=>"KJ":K:Pause1:J◢{Z}:Z"0 NEXT"=0=>U=0:Goto 2△Goto 1△Goto 2GS///高斯法中桩子程序S=K-L:G=(P-R)÷2PQR:A=.1739274226:B=.5-A:C=.0694318442:D=.3300094782:E=1-C:F=1-D: I=4:Lbl 1:B[I]=O+180π-1SB[I](P-1+SGB[I]:Dsz I:Goto 1:O=O+180π-1S(P-1+SG：X=X+S(AcosC+BcosD+AcosE+BcosF:Y=Y+S(AsinC+BsinD+AsinE+BsinFZZ///中桩计算子程序Prog"XL":Prog"GS"WY///外移点计算子程序Lbl 1:I=0:J=90:{IJ}:I"L"≠0=>Z=J"<":J=O+J:O=O+Z:X=X+Rec(I,J:Y=Y+J:Prog"ZB":Goto 1ZB///坐标显示子程序"XY":X:Pause 1:Y◢YC///里程异常处理子程序U=2=>K=LΔU=0=>K=MΔU=1:K<L=>"<<!"◢H=1ΔK>M=>">>!"◢H=1Δ---------------------------以上为程序运算部分,以下为数据库部分-------------------------------------- XL///线路数据库选择子程序Lbl 1:H=0:N"0 SD":N=0=>Prog"0"△N=1=>Prog"1"△N=2=>Prog"2"△...有几条线路仿上行格式输几行H=1=>{NLXYOPQRK}:Goto 10///手工输入子程序L"K0"XYO"A0"Q"LS"P"R0"R"RN":M=L+Q:Prog"YC"1///线路一数据库子程序L=线路起点里程:M=线路终点里程:Prog"YC":H=1=>Goto EΔX=起点X坐标:Y=起点Y坐标:O=起点方位角:Q=线元长:P=起点半径:R=终点半径:K<L+Q=>Goto EΔL=123:X=123:Y=123:O=123:Q=123:P=123:R=123:K≤L+Q=>Goto EΔ...按里程由小到大，每条线元输一行。

CASIO4500坐标计算程序首先，我们需要明确CASIO4500是一款什么样的计算器。

如果CASIO4500是一款用于数学和几何运算的计算器，那么下面是一个用于坐标计算的简单程序。

```pythondef distance(x1, y1, x2, y2):#计算两点之间的距离return ((x2 - x1)**2 + (y2 - y1)**2)**0.5def midpoint(x1, y1, x2, y2):#计算两点连线的中点坐标return (x1 + x2)/2, (y1 + y2)/2def slope(x1, y1, x2, y2):#计算两点连线的斜率return (y2 - y1) / (x2 - x1)def equation_of_line(x1, y1, x2, y2):#计算两点连线的直线方程m = slope(x1, y1, x2, y2)b=y1-m*x1return f"y = {m}x + {b}"def main(:print("CASIO 4500坐标计算程序")print("-----------------------")print("请选择要进行的操作:")print("1. 计算两点之间的距离")print("2. 计算两点连线的中点坐标")print("3. 计算两点连线的斜率")print("4. 计算两点连线的直线方程")print("-----------------------")choice = int(input()if choice == 1:x1 = float(input("请输入第一个点的x坐标: ")) y1 = float(input("请输入第一个点的y坐标: ")) x2 = float(input("请输入第二个点的x坐标: ")) y2 = float(input("请输入第二个点的y坐标: ")) result = distance(x1, y1, x2, y2)print("两点之间的距离为:", result)elif choice == 2:x1 = float(input("请输入第一个点的x坐标: "))y1 = float(input("请输入第一个点的y坐标: ")) x2 = float(input("请输入第二个点的x坐标: ")) y2 = float(input("请输入第二个点的y坐标: ")) result = midpoint(x1, y1, x2, y2)print("两点连线的中点坐标为:", result)elif choice == 3:x1 = float(input("请输入第一个点的x坐标: ")) y1 = float(input("请输入第一个点的y坐标: ")) x2 = float(input("请输入第二个点的x坐标: ")) y2 = float(input("请输入第二个点的y坐标: ")) result = slope(x1, y1, x2, y2)print("两点连线的斜率为:", result)elif choice == 4:x1 = float(input("请输入第一个点的x坐标: ")) y1 = float(input("请输入第一个点的y坐标: ")) x2 = float(input("请输入第二个点的x坐标: ")) y2 = float(input("请输入第二个点的y坐标: ")) result = equation_of_line(x1, y1, x2, y2) print("两点连线的直线方程为:", result)else:print("无效的选择!")if __name__ == "__main__":main```这个程序提供了四种操作：计算两点之间的距离、计算两点连线的中点坐标、计算两点连线的斜率和计算两点连线的直线方程。

CASIO4500坐标计算程序F1: AL1 Lbl 9 {IBGJF}:K=AR兀/180+I“L0”◢ P=I 2/24R-I 4/2384R3◢T=(R+P)tan(A/2)+I/2-I 3/240R 2◢ B“A0”: G“XJD”: J“YJD”:F“Y=1:Z=-1”:Prog XL2 Lbl 2 : {SDH}:L=S“KC”- Z“KZH”: D“DM:1；FW:2”L3 L≤0=>ProgZZ: ≠>L≤I=>ProgHQ:≠>L≤K-I=>ProgYQ:≠>L＜K=>ProgYH:≠>ProgHZ⊿⊿⊿⊿X◢Y◢L4 M=X:U=Y:ProgX:D=1=>ProgR: ≠>D=2=>ProgD⊿⊿Goto2F2: XE“X0”: O“Y0”: M“X1”: U“Y1”: V=Pol(M-E,U-O) ◢W＜0=>N=W+360◢≠>N=W◢F3: HQX=L-L5/40 R2 I2 +I9/3456R4/I4：:Y=L3/60RI-I7/336R3/I3:Q=90L2/兀RI: Prog AF4: YQQ=180(L-I)/兀R+90I/兀R:X=RsinQ+I/2-I3/240R2:Y=R(1-cosQ)+P:ProgAF5: YHL=Abs(K-L):Q=90L2/兀RI:Q=A-Q:N=L-L5/40R2I2:V=L3/6RI:X=T+(T-N)cosA-VsinA:Y=(T-N)sinA+VcosA:Prog A.F6: AV=Pol(T-X,Y):W=B-180-FW:X=G+VcosW:Y=J+V sinW: Q=B+FQF7: ZZL=Abs L+T : X=G-LcosB:Y=J-LsinB:Q=BF8: HZL=L-K+T:Q=B+FA:X=G+LcosQ:Y=J+LsinQF9 R{CD}: C“Q”: D“SD”: Q=Q+C:M“XD”=X+DcosQ◢U“YD”=Y+DsinQ◢ ProgXF10: DL1: {H}:A“WA”: H=H-N: H<0=>U=N-Q-90: ≠>U=N-Q+90⊿H=Abs H: M=V/sin(180-H-U): D=MsinU◢C=MsinU◢C=MsinH◢◢L2: Goto 9 ◢F11: ZHL1: A“XZ”: B“YZ”: C“Q0”L2 Lb1 5: {DK}: D“A0”:K“B0”:F=Pol(D,K):W<0=>G=W+360: ≠>G=W⊿M“XP”=A+Fcos(G+C) ◢U“YP”=B+Fsin(G+C) ◢L3 Prog X : Goto 5F12: SQL1: I“I2”: D“I1”: N“L2”: E“H2”L2: W=I-D: L=RW: T=L/2◢ V=N-T: U=N+T: X=K-V: Y=X2/2R: W>0=>G=E-(N-K)D:H=G+Y◢≠>W<0=>X=K-V:G=E-(N-K)I: H=G-Y◢说明:1、运行F1程序时：A ——转角值 R ——曲线半径L 0 ——缓和曲线长 K ——曲线总长P ——曲线移动量 T ——切线长A0 ——起始段直线方位角 XJD、YJD ——交点坐标Y=1: Z=-1 (转角值为左时输入-1,为右时输入1)X0、Y0 ——置镜点坐标 X1、Y1 ——放样点坐标V ——放样距离 N ——放样方位角KC ——所求点的里程 KZH ——曲线起点里程DM：1 ；FW：2 ——若求曲线上任意点断面上的坐标输入1X、Y ——所求点坐标 Q、S——断面点与线路的夹角和距离V、N ——置镜点到所求点的距离和方位角2、运行F2程序时：X0、Y0 ——测站点坐标 X1、Y1 ——任意点坐标V、N ——所求两点间的距离和方位角3、运行F9程序时：Q ——须输入0（程序内部冲突）SD ——所求点到已知点距离 X、Y ——已知点坐标XD、YD ——所求点坐标4、运行F11程序时：XZ、YZ ——已知点坐标（作为基点）Q0 ——线路方位角 A0 ——所求点到已知点的纵向距离B0 ——所求点到已知点的横向距离XP、YP ——所求点坐标值5、运行F12程序时：I2 ——竖曲线后面坡比 H2 ——竖曲线交点高程I1 ——竖曲线前面坡比 L2 ——竖曲线中心里程R ——竖曲线半径 K——所求点里程 H——所求点高程。

整套卡西欧计算器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.2 线路边桩座标的计算程序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.2 红线交点座标的计算程序6.3 直线与曲线交点座标的计算程序6.4 道路不等宽时曲线元素的计算程序6.5 垂足座标的计算程序6.6 圆曲线敷设的计算程序6.7 缓各曲线敷设的计算程序6.8 极座标放线的计算程序6.9 单一水准路线的近似平差计算程序6.10 丘地有关项目的计算程序公路中、边桩坐标计算及放样程序功能：任意半径曲线中桩、任意交角边桩坐标计算放样、中桩切线方位角、任意点置仪放样、两点间距离计算、两点方位角计算。

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

十一．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——导线全长闭合差、导线全长相对闭合差。

文件名“ZH”:中线程序RS”LS”Z”ZH”FPXY:Lbl 1:{CD K}:CDK:Fixm:C=Z⇒I=0:J=1:W=F:Goto 4⊿M=C －Z:L=Abs M:M<0⇒N=F+180:≠⇒N=F⊿L>S⇒Goto 2⊿V=L－L^5÷(40R2S2)+L^9÷(3456R^4S^4):O=L^3÷(6RS)－L^7÷(336R^3S^3):U=90L2÷(∏RS):Goto 3:Lbl 2:W=180(L－S÷2)÷R÷∏:V=S÷2－S^3÷240÷R2＋Rsin W:O=S2÷24÷R－S^4÷(2384R^3)+R－Rcos W:U=90(2L－S)÷R÷∏:S=0⇒U=W⊿Lbl 3:Pol (V,O):P>0⇒E=N+J:W=N+U:≠⇒E=N－J:W=N－U⊿J=M÷L:Lbi 4:A=X+I cos E+DJcos (W+K) ◢ B=Y+I sin E+DJsin (W+K) ◢Prog “1”:Goto 1文件名：“1”T”ZX”G”ZY”AB:Pol (A－T,B－G):Fixm:J<0⇒J=360+J⊿J◢I◢注：文件名“ZH” (适用于casio4800)输入：R：半径；Ls：缓和曲线长（单圆曲线Ls=0;直线Ls=1E20;R≠0）；ZH：直缓点（直圆点）、（圆直点）缓直点桩号；F：方位角（顺公路走向）；P：判断左偏右偏（左=－1；右=1）；X，Y：直缓点，直圆点，圆直点，缓直点的横纵坐标；C：所求点桩号；D：至中桩距离（左为负，右为正）；K：偏角（正交=90）。

显示：A、B：计算点的X、Y坐标；J、I：计算点到置仪点的方位角、距离。

文件名“1”ZX，ZY：测站点的横纵坐标；A，B：后视点的横纵坐标。