CASIO fx4800P 编程计算器测量程序集

1、CASIO fx系列编程计算器测量程序集编写: 王劲松广东交通职业技术学院2004年9月CASIO fx-4800P 编程计算器测量程序集一SHI JU CE TU （视距测图）1 用途：该程序用于“经纬仪视距法碎部测量”时，计算测站点至碎部点间的平距及碎部点的高程。2 程序数学模型：3 程序清单：C“H0”：I：A“UP”：B“DOWN”：V“MIDDLE”：L：E=90-L D=100Abs（A-B）（cosE）2 H=（C+I）+DtanE-V4 程序说明：H0测站点高程；I测站仪器高；UP塔尺上丝读数；DOWN塔尺下丝读数；MIDDLE塔尺中丝读数；L经纬仪竖盘读数盘左读数；D测站点至

2、碎部点间平距；H碎部点的高程。5该程序若在fx-3950中的程序清单为：？M：？A：？B：？C：？D：100（A-B）（cos(90-D)）2X M+Xtan(90-D)-CY说明：M将测站高程H0+仪器高I；A上丝读数；B下丝读数；C中丝读数；D竖盘盘左读数L。X测站点至碎部点的平距；Y碎部点的高程。二FAN SUAN （反算）1用途：该程序是“坐标反算公式”程序，用于根据两点1、2的平面坐标（X1，Y1）、（X2，Y2），计算两点间的距离D12及坐标方位角A12。2 程序数学模型：3 程序清单：G“X1”：B“Y1”：C“X2”：H“Y2”：E=C-G：F=H-B Pol（E，F）：I“D

3、12” J“A12”4 程序说明：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 程序数学模型：3 程序清单：B“X1”：C“Y1”：D“D12”：A“A12” X“X2”=B+DcosA

4、Y“Y2”=C+DsinA4 程序说明：X1，Y1点1的平面坐标；X2，Y2点2的平面坐标； D12点1至点2间平距；A12点1至点2坐标方位角；X2，Y2点2的平面坐标。另外在输入A12时，可直接输入“度分秒”。四YUAN ZHU （圆主）1用途：该程序是“单圆曲线主点测设元素及里程计算”程序。2程序数学模型：3程序清单：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+LG“QZ”=F-L÷24程序说明：TURNING

5、ANGLE交点JD处的转角；R圆曲线半径；JD交点。T切线长；L圆曲线长；E外矢距；D切曲差；ZYZY点里程；YZYZ点里程；QZQZ点里程。五YUAN XIANG （圆详）1用途：该程序是“单圆曲线切线支距法及长弦偏角法详细测设”程序。2程序数学模型：3程序清单：A“ZY OR YZ”：R：Lbl 1 L：B=180Abs(L-A) ÷R÷：X=RsinBY=R-RcosB D“ANGLE”=B÷2C“DISTANCE”=2RsinDGoto 14程序说明：ZY OR YZZY点或YZ点里程（桩号）；R圆曲线半径； L某桩的里程（桩号）；X切线支距法的X值；Y切

6、线支距法的Y值；ANGLE偏角法的偏角值，单位为“度”。若要查看其“度分秒”，可在程序运行显示时先后按“SHIFT”、“度分秒”这2个键；DISTANCE长弦偏角法的弦长。六HUAN ZHU （缓主）1用途：该程序是“完整对称带缓和曲线的圆曲线”的主点测设元素及里程计算程序。2程序数学模型：3程序清单：A“TURNING ANGLE”：R：S“LS”：H“JD” B=180S÷（2 R）P=S2÷24÷R-S4÷(2688 R 3) Q=S÷2-S3÷(240R2) T=(R+P)tan(A÷2)+Q M=R(A-2B)&#

7、247;180L=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+S4程序说明：TURNING ANGLE交点JD处的转角；R圆曲线半径；LS缓和曲线长；JD交点里程（桩号）。T切线长；L曲线长；E外矢距；D切曲差；ZHZH点里程；HYHY点里程；QZQZ点里程；YHYH点里程；HZHZ点里程。七HUAN QIE （缓切）1用途：该程序是“完整对称带缓和曲线的圆曲线”的切线支距法详细测设坐标计算程序。2程序数学模型：3程序清单：A“ZH”：R：S“LS

8、”：Lbl 1 L，B：Lbl 2C=Abs(L-A)：D=RS：X=C-C5÷40D2+C9÷3456D4-C13÷599040D6+C17÷17542600D8Y=C3÷6D-C7÷336D3+C11÷42240D5-C15÷9676800D7+C19÷3530097000D9 Goto 1Lbl 3E=180(Abs(L-A)-S)÷R÷+180S÷(2R)：P=S2÷24÷R-S4÷2688÷R3：Q=S÷2-S3&#

9、247;240÷R2X=RsinE+QY=R-RcosE+PGoto 14程序说明：ZHZH点桩号（里程）；R圆曲线半径；LS缓和曲线长；L待测设桩的桩号（里程）；B当待测设中桩位于缓和曲线段，则输入“1”，当待测设中桩位于圆曲线段，则输入“1”以外的数值。X切线支距法的X值；Y切线支距法的Y值。八HUAN PIAN （缓偏）1用途：该程序是“完整对称带缓和曲线的圆曲线”的偏角法详细测设计算程序。2程序数学模型：3程序清单：A“ZH”：R：S“LS” B=30S÷R÷：F“B0”=2BLbl 1 LE=Abs(L-A)：D“ANGLE”=（E÷S）2BC

10、“DISTANCE”=E-E5÷(90R2S2) Goto 14程序说明：ZHZH点桩号（里程）；R圆曲线半径；LS缓和曲线长；B0在以后要测设圆曲线段上的中桩时，架仪于HY（或YH）点，瞄准ZH（或HZ）点，拨角B0为圆曲线的切线方向，找到切线方向后，即按单圆曲线偏角法去测设即可； L待测设桩的桩号（里程）；ANGLE偏角法的偏角值，在程序显示中其单位为“度”。若要查看其“度分秒”，可在显示时先后按“SHIFT”、“度分秒”这2个键；DISTANCE长弦偏角法的弦长。九ZHUAN HUAN （转换）1用途：该程序是“两平面坐标系间坐标转换”的计算程序。2程序数学模型：3程序清单：C

11、“X0”：E“Y0”：D“ANGLE”：F“SIGN” Lbl 0 A，BF1A=A：B=-B X=C+AcosD-BsinDY=E+BcosD+AsinDGoto 04程序说明：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用途：该

12、程序是“附合导线平差”的计算程序。2程序数学模型：3程序清单：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 Z2N-1=F：Z2N=D：M=M+D：N<C+1Goto 0 L：F=F+L：FGF=F+G：F=F-G F2GF+F-2G B“FW”=F-A R：H“FR”=R（C+2） B=-B÷（C+2）N=0：P=0：Q=0：Lbl 1 N=N+1Z2N-1“A”=Z2N-1+NBI=

13、Z2NcosZ2N-1：P=P+I：Z2C+2N+1=I：J= Z2NsinZ2N-1：Q=Q+J：Z2C+2N+2=J：N<C+1Goto 1 F“FX”=P-S+X G“FY”=Q-E+Y H“F”=（F2+G2） K=H÷MN=0：V=-F÷M：W=-G÷M：Lbl 2 N=N+1Z2C+2N+1= Z2C+2N+1+VZ2N ：Z2C+2N+2= Z2C+2N+2+WZ2N：X=X+ Z2C+2N+1 Y=Y+Z2C+2N+2 N<CGoto 2：“END”4程序说明：X0，Y0，A0已知的起始点的纵、横坐标及起始边的坐标方位角。X1，Y1，A

14、1已知的结束点的纵、横坐标及结束边的坐标方位角。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”、“度分

15、秒”这2个键。十一BI DAO XIAN（闭导线）1用途：该程序是“闭合导线平差”的计算程序。2程序数学模型：3程序清单：Defm 40：X“X0”：Y“Y0”：F“A0”：N“UNKNOWN POINT”：N=N+1I=0：E=0：C=0：Lbl 0 I=I+1：L，D：Z2I-1=L：Z2I=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：Z2N+1=F：Lbl 1 I+I+1 F=F+Z2I-1：F180F=F+180：F=F-180 F360F+F-360 Z2N+2I-1“A”=F+（I-1

16、）B I<N Goto 1 I=0：G=0：H=0：Lbl 2 I+I+1：Z2I-1=Z2IcosZ2N+2I-1：G=G+Z2I-1：Z2N+2I=Z2Isin Z2N+2I-1：H=H+Z2N+2I：I<N Goto 2 G“FX” H“FY” I“F”=（G2+H2） K=I÷C U=-G÷C：V=-H÷CI=0：Lbl 3 I+I+1 Z2I-1=Z2I-1+U Z2I：Z2N+2I=Z2N+2I+V Z2I：X=X+ Z2I-1 Y=Y+ Z2N+2I I<N-1Goto 3：“END”4程序说明：X0，Y0，A0已知的起始点的纵、横

17、坐标及起始边的坐标方位角。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个键。十二SH

18、UI ZHUN（水准）1用途：该程序是“单一水准路线的平差”的计算程序。2程序数学模型：3程序清单：A“CODE”：Defm 20：G=0：M=0：I=0：C“UNKNOWN POINT”Lbl 0：I=I+1：H：Z2I=H：G=G+H A=1L：Z2I-1=L：M=M+L：N：Z2I-1=N：M=M+N I<C+1Goto 0 B“H0”：D“H1”：F=B+G-D J“FR”=RM÷1000 AbsF<J“CHAO XIAN CHA”：V=-F÷M：I=0：H=B：Lbl 1：I=I+1 H=H+Z2I+VZ2I-1 I<CGoto 1：“END”4

19、程序说明：CODE地形编码，当是平地，统计公里数时输入“1”；当是山地，统计测站数时输入非“1”。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输入相

20、同的值。十三FANG PIN GAO（方平高）1用途：该程序是“方格网法平整场地，填挖方量基本平衡时的设计高程”的计算程序。2程序数学模型：3程序清单：N：X“INPUT ANGLE NUMBER”：I=0：C=0：Lbl 0I=I+1：H：C=C+H：I<X Goto 0 M“INPUT EDGE NUMBER”：J=0：D=0：Lbl 1J“I”=J+1：H：D=D+2H：J<M Goto 1 P“INPUT CORNER NUMBER”：K=0：E=0：Lbl 2K“I”=K+1：H：E=E+3H：K<P Goto 2 Q“INPUT MIDDLE NUMBER”：L=

21、0：F=0：Lbl 3L“I”=L+1：H：F=F+4H：L<Q Goto 3 G“H0”=（C+D+E+F）÷4÷N4程序说明：N方格的个数。INPUT ANGLE NUMBER输入角点的个数。INPUT EDGE NUMBER输入边点的个数。INPUT CORNER NUMBER输入拐点的个数。INPUT MIDDLE NUMBER输入中点的个数。H各点的实测高程。H0填挖方量基本平衡时的设计高程。十四FANG TU FANG（方土方）1用途：该程序是“方格网法平整场地时，填、挖方量”的计算程序。2程序数学模型：3程序清单：D：E“H0”：X“INPUT ANGLE NUMBER”：G=0：J=0：A=0：Lbl 1J=J+1：H：F“DH”=H-E F>0G=G+F÷4：A=A+ F÷4 J<X Goto 1 Y“INPUT EDGE NUMBER”：L=0：K=0：B=0：Lbl 2：K=K+1：H：M“DH”=H-E M>0L=L+M÷2：B=B+ M÷2 K<Y Goto 2 Z“INPUT CORNER NUMBER”：P=0：Q=0：C=0：Lbl 3Q=Q+1：H：T“DH”=H-E T>0P=P+3T÷4：C=C+ 3T÷4