CASIO fx系列编程计算器测量程序集编写: 王劲松广东交通职业技术学院2004年9月第一部分：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——竖盘盘左读数LX——测站点至碎部点的平距；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”◢ J“A12”3． 程序说明：X1，Y1——点1的平面坐标；X2，Y2——点2的平面坐标； D12——点1至点2间平距；A12——点1至点2坐标方位角由于编程时用到了极坐标函数POL（），故在程序运行中，会出现I=？XXX（数字）及J=？XXX（数字）这是函数POL计算出的两个量，不能改变，直接回车默认即可三．ZHENG SUAN （正算）1．用途：该程序是“坐标正算公式”程序，用于根据点1的平面坐标（X1，Y1），及点1至点2的平距D12、坐标方位角A12，来计算点2的平面坐标（X2，Y2）2． 程序清单：B“X1”：C“Y1”：D“D12”：A“A12” ↙X“X2”=B+DcosA ◢ Y“Y2”=C+DsinA3． 程序说明：X1，Y1——点1的平面坐标；X2，Y2——点2的平面坐标； D12——点1至点2间平距；A12——点1至点2坐标方位角；X2，Y2——点2的平面坐标另外在输入A12时，可直接输入“度分秒”四．YUAN ZHU （圆主）1．用途：该程序是“单圆曲线主点测设元素及里程计算”程序。
2．程序清单：A“TURNING ANGLE”：R：B“JD” ↙T=Rtan(A÷2)◢ L=RAπ÷180◢ E=R÷cos(A÷2)-R◢ D=2T-L◢ C“ZY”=B-T◢ F“YZ”=C+L◢G“QZ”=F-L÷23．程序说明：TURNING ANGLE——交点JD处的转角；R——圆曲线半径；JD——交点T——切线长；L——圆曲线长；E——外矢距；D——切曲差；ZY——ZY点里程；YZ——YZ点里程；QZ——QZ点里程五．YUAN XIANG （圆详）1．用途：该程序是“单圆曲线切线支距法及长弦偏角法详细测设”程序2．程序清单：A“ZY OR YZ”：R：Lbl 1 ↙{L}：B=180Abs(L-A) ÷R÷π：X=RsinB◢Y=R-RcosB◢ D“ANGLE”=B÷2◢C“DISTANCE”=2RsinD◢Goto 13．程序说明：ZY OR YZ——ZY点或YZ点里程（桩号）；R——圆曲线半径； L——某桩的里程（桩号）；X——切线支距法的X值；Y——切线支距法的Y值；ANGLE——偏角法的偏角值，单位为“度”若要查看其“度分秒”，可在程序运行显示时先后按“SHIFT”、“度分秒”这2个键；DISTANCE——长弦偏角法的弦长。
六．HUAN ZHU （缓主）1．用途：该程序是“完整对称带缓和曲线的圆曲线”的主点测设元素及里程计算程序2．程序清单：A“TURNING ANGLE”：R：S“LS”：H“JD” ↙B=180S÷（2 Rπ）↙P=S2÷24÷R-S^4÷(2688 R ^3) ↙Q=S÷2-S^3÷(240R2) ↙T=(R+P)tan(A÷2)+Q◢ M=R(A-2B)π÷180↙L=M+2S◢ E=(R+P)÷cos(A÷2)-R◢ D=2T-L ◢ N“ZH”=H-T◢ O“HY”=N+S◢ U“QZ”=O+M÷2◢ V“YH”=O+M◢ W“HZ”=V+S3．程序说明：TURNING ANGLE——交点JD处的转角；R——圆曲线半径；LS——缓和曲线长；JD——交点里程（桩号）T——切线长；L——曲线长；E——外矢距；D——切曲差；ZH——ZH点里程；HY——HY点里程；QZ——QZ点里程；YH——YH点里程；HZ——HZ点里程七．HUAN QIE （缓切）1．用途：该程序是“完整对称带缓和曲线的圆曲线”的切线支距法详细测设坐标计算程序2．程序清单：A“ZH”：R：S“LS”：Lbl 1 ↙{L，B}↙：↙Lbl 2↙C=Abs(L-A)：D=RS：X=C-C^5÷40D2+C^9÷3456D^4-C^13÷599040D^6+C^17÷17542600D^8◢Y=C^3÷6D-C^7÷336D^3+C^11÷42240D^5-C^15÷9676800D^7+C^19÷3530097000D^9◢ Goto 1↙Lbl 3↙E=180(Abs(L-A)-S)÷R÷π+180S÷(2πR)：P=S2÷24÷R-S^4÷2688÷R^3：Q=S÷2-S^3÷240÷R2↙X=RsinE+Q◢Y=R-RcosE+P◢Goto 1↙3．程序说明：ZH——ZH点桩号（里程）；R——圆曲线半径；LS——缓和曲线长；L——待测设桩的桩号（里程）；B——当待测设中桩位于缓和曲线段，则输入“1”，当待测设中桩位于圆曲线段，则输入“1”以外的数值。
X——切线支距法的X值；Y——切线支距法的Y值八．HUAN PIAN （缓偏）1．用途：该程序是“完整对称带缓和曲线的圆曲线”的偏角法详细测设计算程序2．程序清单：A“ZH”：R：S“LS” ↙B=30S÷R÷π：F“B0”=2B◢Lbl 1 ↙{L}↙E=Abs(L-A)：D“ANGLE”=（E÷S）2B◢C“DISTANCE”=E-E^5÷(90R2S2) ◢Goto 13．程序说明：ZH——ZH点桩号（里程）；R——圆曲线半径；LS——缓和曲线长；B0——在以后要测设圆曲线段上的中桩时，架仪于HY（或YH）点，瞄准ZH（或HZ）点，拨角B0为圆曲线的切线方向，找到切线方向后，即按单圆曲线偏角法去测设即可； L——待测设桩的桩号（里程）；ANGLE——偏角法的偏角值，在程序显示中其单位为“度”若要查看其“度分秒”，可在显示时先后按“SHIFT”、“度分秒”这2个键；DISTANCE——长弦偏角法的弦长九．ZHUAN HUAN （转换）1．用途：该程序是“两平面坐标系间坐标转换”的计算程序2．程序清单：C“X0”：E“Y0”：D“ANGLE”：F“SIGN” ↙Lbl 0↙ {A，B}↙F1A=A：B=-BΔ X=C+AcosD-BsinD◢Y=E+BcosD+AsinD◢Goto 03．程序说明：X0，Y0——施工坐标系（A-O’-B）的原点O’在统一坐标系（x-o-y）中的坐标。
ANGLE——为统一坐标系的x轴顺时针旋转至施工坐标系的A轴的角值SIGN——为符号函数，若输入“1”时，则表明x-o-y为左手系，且A-O’-B也为左手系；若输入“1”之外值，则表明x-o-y为左手系，而A-O’-B为右手系A，B——某点在施工坐标系中的纵、横坐标X，Y——该点在相应统一坐标系中的纵、横坐标十．FU DAO XIAN（附导线）1．用途：该程序是“附合导线平差”的计算程序2．程序清单：Defm 40：X“X0”：Y“Y0”：F“A0”：S“X1”：E“Y1”：A“A1”：C“UNKNOWN POINT”↙N=0：M=0：Lbl 0 ↙N=N+1：{L，D}：G=180：F=F+L：FGF=F+G：F=F-G Δ F2GF+F-2G Δ Z[2N-1]=F：Z[2N]=D：M=M+D：N
X1，Y1，A1——已知的结束点的纵、横坐标及结束边的坐标方位角UNKNOWN POINT——待求未知点的个数L，D——观测的导线各左转角及平距FW，R，FR——导线角度闭合差、角度闭合差的允许误差的系数及允许误差值N，A——各条导线边序号及其坐标方位角FX，FY——坐标增量ΔX、ΔY的闭合差F，K——导线全长闭合差、导线全长相对闭合差N，X，Y——各导线点的序号及其纵横坐标注：（1）本程序是按10个未知点来设计的，若未知点数超过10个时，按“4倍点数”来改变源程序中的“Defm 40”2）显示的角度允许误差FR的单位为“秒”3）在程序显示中角度单位为“度”若要查看其“度分秒”，可在显示时先后按“SHIFT”、“度分秒”这2个键十一．BI DAO XIAN（闭导线）1．用途：该程序是“闭合导线平差”的计算程序2．程序清单：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
UNKNOWN POINT——待求未知点的个数L，D——观测的导线各左转角及平距FW，R，FR——导线角度闭合差、角度闭合差的允许误差的系数及允许误差值I，A——各条导线边序号及其坐标方位角FX，FY——坐标增量ΔX、ΔY的闭合差F，K——导线全长闭合差、导线全长相对闭合差I，X，Y——各导线点的序号及其纵横坐标注：（1）本程序是按10个未知点来设计的，若未知点数超过10个时，按“4倍点数”来改变源程序中的“Defm 40”2）显示的角度允许误差FR的单位为“秒”3）在程序显示中角度单位为“度”若要查看其“度分秒”，可在显示时先后按“SHIFT”、“度分秒”这2个键十二．SHUI ZHUN（水准）1．用途：该程序是“单一水准路线的平差”的计算程序2．程序清单：A“CODE”：Defm 20：G=0：M=0：I=0：C“UNKNOWN POINT”↙Lbl 0：I=I+1：{H}：Z[2I]=H：G=G+H ↙A=1{L}：Z[2I-1]=L：M=M+L：{N}：Z[2I-1]=N：M=M+N Δ I
UNKNOWN POINT——待求高程点的个数L——水准点间的路线长N——水准点间的测站数H——水准点间的实测高差H0，H1——已知的起始点和结束点的高程F，FR——高差闭合差、高差允许闭合差R——高差允许闭合差系数等外水准：12（按测站数）或40（按公里数）；四等水准：6（按测站数）或20（按公里数）I，H——第I段及第I个待求水准点的高程注：（1）本程序是按10个未知点来设计的，若未知点数超过10个时，按“2倍点数”来改变源程序中的“Defm 20”2）当水准路线是闭合水准路线时，只需将H0，H1输入相同的值十三．FANG PIN GAO（方平高）1．用途：该程序是“方格网法平整场地，填挖方量基本平衡时的设计高程”的计算程序2．程序清单：N：X“INPUT ANGLE NUMBER”：I=0：C=0：Lbl 0↙I=I+1：{H}：C=C+H：I
INPUT ANGLE NUMBER——输入角点的个数INPUT EDGE NUMBER——输入边点的个数INPUT CORNER NUMBER——输入拐点的个数INPUT MIDDLE NUMBER——输入中点的个数H——各点的实测高程H0——填挖方量基本平衡时的设计高程十四．FANG TU FANG（方土方）1．用途：该程序是“方格网法平整场地时，填、挖方量”的计算程序2．程序清单：D：E“H0”：X“INPUT ANGLE NUMBER”：G=0：J=0：A=0：Lbl 1↙J=J+1：{H}：F“DH”=H-E◢ F>0G=G+F÷4：A=A+ F÷4 Δ J0L=L+M÷2：B=B+ M÷2 Δ K0P=P+3T÷4：C=C+ 3T÷4 ΔQ0R=R+U：I=I+U Δ S
H0——平整场地时的设计高程INPUT ANGLE NUMBER——输入角点的个数INPUT EDGE NUMBER——输入边点的个数INPUT CORNER NUMBER——输入拐点的个数INPUT MIDDLE NUMBER——输入中点的个数H——各点的实测高程DH——填挖高度（各点的实测高程与设计高程之差），正为挖，负为填VW ，VT——挖方、填方十五．JI ZHUO BIAO （极坐标放样程序）1．用途：该程序是用于根据测站点坐标(x1,y1)和后视点坐标(x2,y2)，计算在测站点后视后视点，用极坐标法放样任一待定点(x,y)所需正拨的水平角及边长2． 程序清单：G“X1”：B“Y1”：C“X2”：H“Y2”：E=C-G：F=H-B ↙Pol（E，F）：↙ L=J：Lbl 1↙{X，Y}：M=X-G：N=Y-B↙Pol（M，N）：D=I◢↙ A=J-L：A◢ Goto 13．程序说明：X1，Y1——测站点的坐标X2，Y2——后视点的坐标X，Y——待放样点的坐标A——测站点后视后视点，旋转至待放样点所需正拨的水平角角值D——测站点至待放样点的水平距离由于编程时用到了极坐标函数POL（），故在程序运行中，放样第一个点时，会出现J=？XXX（数字）及I=？XXX（数字）这是函数POL计算出的两个量，不能改变，直接回车默认即可。
放样以后几个点时，不会出现此现象第二部分：CASIO fx-4500P 编程计算器施工测量程序名代码变 量 与 说明2曲线要素：P=I2/24R-I4/2688R3M=I/2-I3/240R2T=(R+P)Tan(A/2)+ML=A/180×πR+IB=I/R×90/πA,R,I,m,p,B,T,L 本系统的变量是统一的，代表特定含义，在各段代码中通用，不宜更改 A:曲线转角R:曲线半径I:缓和曲线长T:切线长L:曲线长m,p,B:曲线要素 3直线相交（四点）点:NSOTPUQVPol(Q-P,V-U)Q=U-PtanWU=WPol(O-N,T-S)N=T-OtanWO=W▲ U▲A=U-O:Pol(cosA,sinA):A=AbsW▲ K=W/A▲ J=(Q-N)/(tanO-tanU) ▲ D=JtanO+N▲O,U,J,D连续输入四点坐标 两条直线方位转角交点坐标O,U:第一、第二条直线的方位（计算后）J,D:交点坐标 （计算后）31交点方位:JDOUA=U-OPol(cosA,sinA)A=AbsAK=W/AJ,D,O,U 输入交点坐标与切线方位转角转向K:曲线转向,右转为正4任意点坐标：Lbl 2:{S}S≤0→G=S:H=0:C=0:Goto 4△ S≥L→G=L-S:H=0:C=0: Goto 3△ S→Prog 5: Goto 4△ S→C=B+180(S-I)/RπG=RsinC+M:H=R-RcosC+P: Goto 4≠>S=L-S:Prog 5Lbl 3X=J-KhsinU+cosU(T-G) ▲ Y=D+KhcosU+sinU(T-G) ▲ Z=U-KC:Goto 5Lbl 4X=J-KhsinO-cosO(T-G) ▲ Y=D+KhcosO-sinO(T-G) ▲ Z=O+KCLbl 5:{N}:G=X+KNsinZ▲ H=Y-KNcosZ▲ Goto 2S,X,Y,G,H,N 输入到曲线起点的弧长坐标坐标方向距离（外侧为正）法线点坐标S:距曲线起点的弧长,S可用到曲线起点的里程差得到N:法向距离(曲线外侧为正)X,Y:中线点坐标Z:切线方位 G,H:法向点坐标 41程序连接:Prog 31:Prog 2:Prog 4 48曲线长反算:Z=O+(90+A/2)KS=(R+P)/cos(A/2)G=J+ScosZH=D+SsinZLbl 8:{XY}Pol(X-G,Y-H)Z=K(W-O)+90S=RZ×/180+I/2▲Goto 8X,Y,SOAKRPJD为曲线变量 法向点坐标曲线长由法向点坐标反算成对应曲线长49E=G+G5/40Q2-G9/3456Q4F=H-E3/6Q+E7/336Q3E2/2Q×F/(1-EF/Q)+E 缓和曲线段的反算5缓和曲线:Q=RIG=S-S5/40Q2+S9/3456Q2Q2H=S3/6Q-S7/336Q3+S11/42240Q5C=S2/RI×90/πG,H,C6放样: Defm 2Z[1] ▲ Z[1]=Ans:Z[2] ▲ Z[2]=AnsEFPol(Z[1]-E,Z[2]-F)Q=WLbl 0:{XY}Pol(X-E,Y-F) ▲W-QAns<0→Ans+360△”B=”▲Goto 0 输入后视点坐标输入测站坐标EF 放样距离V放样顺拨角度(度)Z[1],Z[2]:后视点坐标E,F: 测站点坐标Q:后视方位X,Y: 放样点坐标8单导线:KLbL 1:{AS}Prog A:SB=B+180-KA▲X=X+ScosB▲ Y=Y+SsinB▲ Goto 1A,S,B,X,Y输入左角右角输入角度距离 前视方位前视点坐标连续计算单导线坐标，角度按度.分秒输入左角K=-1,右角K=181支导线:ZSZ=Z+QX=E+ScosZ▲ Y=F+SsinZ▲Z,S,Q,X,Y前视角度(度),距离前视方位前视坐标放样时计算前视点的坐标,以便反算A角度转换:A=IntA+Int(100FracA)/60+Frac(100A)/36A输入度.分秒，转换成度13。