最新CASIOfx4500P计算器程序.docx

1、最新CASIOfx4500P计算器程序CASIOfx4500P计算器程序CASIO fx4500P计算器实 用 程 序XCAJCJ 编著二一年八月说明本程序稿由于时间仓促及编者水平有限，缺点、错误以及疏漏与不足之处在所难免，恳请大家批评指正。程序步骤计算：1个字符的文件名算3步(每增加一个字符多算一步)，行号算1步，每个字符算1步。CASIO fx-4500P共有1103步。XCAJCJ二一年八月目录程序1角度判断 1程序2角度换算 1程序3角度换算 1程序4角度换算 1程序5角度换算 1程序6视距测量 2程序7坐标反算 2程序8坐标方位角计算 2程序9坐标正算 3程序10坐标方位角交会计算

2、3程序11垂距计算 4程序12视差法测距计算 5程序13剖面交点计算 6程序14前方交绘计算 6程序15附合导线近似平差 7程序16渠道中线里程计算 9程序17多边形面积计算 11程序18闭合导线近似平差 13程序19直线交点坐标计算 14程序20极坐标法放样计算 15附录1：CASIO fx-4500P错误标示一览表 16附录2：常用数学公式 17附录3：常用测量公式 17附录4：单一附合导线内业平差计算 20程序1角度判断文 件 名A角度判断29步L1B360B=B360B0B=B+360说明：程序A为子程序，对360或0的角值进行判断处理。 程序2角度换算文 件 名B角度换算23步L1B

3、=Int B+Frac B/.6+Frac100B/90程序3角度换算文 件 名C角度换算22步L1B=Int B+.6Frac B.004Frac 60B程序4角度换算文 件 名D角度换算25步L1F=B:B=I:Prog B:I=B:B=F:Prog B程序5角度换算文 件 名E角度换算25步L1F=B:B=I:Prog B:I=B:B=F:Prog B说明：程序B、C、D、E均为子程序，简化测量计算中角度的输入及显示。注：子程序应先输入CAISOfx-4500P中，以备其它程序调用。程序6视距测量文 件 名1视距测量41步L1Lbl 0 :SCB:HSC:Prog B:S(sin B)2

4、Ans /tan B+HCProg C: Goto 0说明：程序用以计算视距测量中测站至照准点的平距及高程；程序中：H=测站高程+仪高(程序只要求输入该值一次)；S=视距、C=切尺、B=竖盘读数；算例：按Prog 1调用程序显示H？输入测站高程+仪高445.79+1.41，按EXE键、显示S？输入视距39，按EXE键、显示 C？输入切尺2.9，按EXE键、显示B？输入竖盘读数90.12，按EXE键、显示(计算结果)平距38.9995248，按EXE键、显示高程444.1638654，按EXE键、显示S？、输入视距竖盘读数90.12是将度分秒值(901200)以小数的形式输入；其它程序的角度值输

5、入及显示均如此。程序7坐标反算文 件 名2视标反算33步L1AX0CY0:Pol(XA,YCB=W:Porg A:Prog C说明：程序用以计算两点间的距离及坐标方位角；程序中：A=起点纵坐标X0，C=起点横坐标Y0，X=终点纵坐标，Y=终点横坐标；算例：按Prog 2调用程序显示X0？输入起点纵坐标667.132，按EXE键、显示Y0？输入起点横坐标993.505，按EXE键、显示X？、输入终点纵坐标395.372，按EXE键、显示Y？输入终点横坐标914.821，按EXE键、显示两点间距离282.9216666，按EXE键、显示起点至终点的方位角196.0851267。程序8坐标方位角计算

6、文 件 名3坐标方位角计算26步L1Prog D:B=B+I180:Prog A:D=B:Prog C说明：程序用以推算导线边的坐标方位角(左角公式)；程序中：B=已知方位角、I=转折角(左角)；语句B=D是为其它程序调用本程序而预设的，此处无具体用途；算例：按Prog 4调用程序显示B？输入已知方位角120.10214按EXE键、显示I？输入转折角160.01217，按EXE键、显示所求方位角100.11431；若转折角为右角，将语句B=B+I180改为B=BI+180即可。程序9坐标正算文 件 名4坐标正算22步L1Prog B:X=Rec(L, B)+XY=Y+W说明：程序根据已知坐标、

7、边长及其坐标方位角计算待定点坐标；程序中：B=坐标方位角、L=边长、X=已知点纵坐标、Y=已知点横坐标；算例：按Prog3调用程序显示B？输入坐标方位角100.11431，按EXE键、显示L？输入边长151.213，按EXE键、显示X？输入已知点纵坐标1000，按EXE键、显示待定点纵坐标973.2346789，按EXE键、显示Y？输入已知点横坐标2000，按EXE键、显示待定点横坐标2148.825364。程序10坐标方位角交会计算文 件 名5坐标正算100步L1Lbl 0: I :Prog D:S=B+I:AX0CYOXYN:B=N:Prog B:E=(AtanSXtanB +YC)/(t

8、anStanBC+tanS(EAPol(AE,CAnsB=BW:Prog A:Prog CB=F:Goto 0说明：程序利用已知方位角，计算测站上任意方向线与轴线相交的距离及交点坐标和该交点作为测站后与轴线的夹角；如下图，程序中：B=已知方位角、I=观测角、A(X0)=测站纵坐标、C(Y0)=测站横坐标、X=轴线上任意点纵坐标、Y=轴线上任意点横坐标、N=轴线方位角；算例：按Prog5调用程序显示B？输入已知方位角150.01252，按EXE键、显示I？输入观测角110.30，按EXE键、显示X0？输入测站纵坐标100，按EXE键、显示Y0？输入测站横坐标100，按EXE键、显示X？输入轴线方

9、向上任意点纵坐标80，按EXE键、显示Y？输入任意点横坐标150，按EXE键、显示N？输入轴线方位角210.12374，按EXE键、显示交点纵坐标111.3973848，按EXE键、显示交点横坐标168.2813366，按EXE键、显示测点至交点的距离69.22601608，按EXE键、显示该交点作为测站后与轴线的夹角09.41122，按EXE键、显示I？输入观测角程序所用正切公式如图所示：已知两方向、两点坐标，求其交点P的坐标；程序11垂距计算文 件 名6垂距计算27步L1Prog B:D=B:Prog 2B=WD:Vcos BVsin B说明：程序用以计算勘探点偏离轴线的距离；如图，程序中

10、：B=轴线坐标方位角、A(X0)=轴线控制点的纵坐标、C(Y0)=轴线控制点的横坐标、X=勘探点纵坐标、Y=勘探点横坐标；算例：按Prog 6调用程序显示B？输入轴线坐标方位角301.09505，按EXE键、显示X0？输入轴线上控制点的纵坐标138.579，按EXE键、显示Y0？输入轴线上控制点的横坐标1321.474，按EXE键、显示X？输入勘探点的纵坐标341.871，按EXE键、显示Y？输入勘探点的横坐标986.408，按EXE键、显示控制点至勘探点的距离391.914355(斜距)，按EXE键、显示控制点至勘探点的坐标方位角301.1446034，按EXE键、显示控制点至勘探点的距离(

11、平距)391.9139527，按EXE键、显示勘探点偏离轴线的距离0.561530147(垂距)；如果勘探点偏离轴线的距离为正值，表示该点在轴线右侧；若为负值，则该点在轴线左侧。程序12视差法测距计算文 件 名7视差法测距计算30步L1Prog D:KTH:1/tan(B/2H=Ans tan IHKT说明：程序用以计算2M横基尺视差法测距的边长及高程；程序中：B=视差角、I=垂直角、K=仪高、T=觇标高、H=测站高程；算例：按Prog 7调用程序显示B？输入视差角2.33176，按EXE键、显示I？输入垂直角1.45386，按EXE键、显示K？输入仪高1.432，按EXE键、显示T？输入觇标

12、高1.405，按EXE键、显示H？输入测站高程443.425，按EXE键、显示测站至测点的距离44.84444364，按EXE键、显示测点高程444.8305217。程序13剖面交点计算文 件 名8剖面交点计算55步L1Prog B:L:D=B:B=J:Prog D:N=DB:M=ID:Abs(LsinM/sin(M+NAbs(Lsin N/sin(M+N说明：程序8用以计算剖面交点距离；输入坐标方位角时不用考虑方向性；算例：如图,按Prog 8调用程序显示B？输入坝1的坐标方位角96.0303或276.0303，按EXE键、显示L？输入坝1交坝围的距离21.9，按EXE键、显示J？输入坝的坐

13、标方位角301.0950或121.0950，按EXE键、显示I？输入围的坐标方位角307.1202或127.1202，按EXE键、显示坝交围及坝1的距离107.7198781，按EXE键、显示围交坝1坝的距离88.38012441。程序14前方交绘计算文 件 名9前方交绘计算65步L1AX0CY0:Prog D:M=1/tan B:N=1/tan I:(AN+XMC+Y)/(M+N(CN+Y M+AX)/(M+N说明：程序9用以计算前方交会法中待定点的坐标；程序采用余切公式计算坐标，三角形点号应按逆时针编排。若是顺时针排列，观测角应以负值输入；算例，如图表所示：逆时针顺时针Xa2235.6Ya

14、845.0Xa6735.7Ya10738.5Xb6735.7Yb10738.5Xb2235.6Yb845.0cab48.5221cba22.4045cab22.4045cba48.5221Xc6467.705Yc2112.188Xc6467.705Yc2112.188交会角C在3050之间时交会误差最小。程序15附合导线近似平差文 件 名10附合导线近似平差315步L1C=0:A9Defm 2(A9Lbl l:C=C+1:I:Prog 3:AC+8=D:CAGoto 1B:Prog B:B=DB:D=B/A:Prog CC=0L2Lbl 2:C=C+1:AC+8=AC+8DC:B=AC+8:P

15、rog CCAGoto 2C=8:D=0:E=0:F=0L3Lbl 3:C=C+1:B:D=D+B:E=E+Bcos AC:F=F+Bsin AC:AA+C =B:CA +7Goto 3DBI:BGHI:E=E+BHF=F+GID/(E2+F2E=E/D:F=F /DL4C=8:Lbl 4:C=C+1:B=AA+Ccos AC+BEAA+CG=AA+C sin AC+GFAA+CCA+7Goto4Defm 0说明：程序10用以单一附合导线的近似平差计算；程序中：A=测站数、B=起算边坐标方位角、附合边坐标方位角、导线边边长、起点纵坐标、I=观测角(左角)、终点横坐标、G=起点横坐标、H=终点纵

16、坐标；程序设计可以计算9个测站，大于是9站时只要计算器剩余有步骤数，程序会自动扩充储存器增加储存变量；算例如表：按Prog 10调用程序显示A？输入测站数5，按EXE键、显示B？输入起算边的坐标方位角169.5359，按EXE键、显示I？输入B点的观测角93.3642，按EXE键、显示I？输入C点的观测角179.56，依次输入各观测角，按EXE键、显示B？输入附合边的坐标方位角250.06，按EXE键、显示导线的角度闭合差0.003699999(37)，按EXE键、显示平差后B点至C点的坐标方位角83.30484，按EXE键依次显示各导线边平差后的坐标方位角，按EXE键、显示B？输入B点至C点

17、的边长126.11，按EXE键、显示B？依次输入各导线边的边长，按EXE键、显示各边长之和508.96，按EXE键、显示B？输入起点纵坐标801.53，按EXE键、显示G？输入起点横坐标724.84，按EXE键、显示H？输入终点纵坐标959.29，按EXE键、显示I？输入终点横坐标1158.41，按EXE键、显示纵坐标闭合差0.093449934，按EXE键、显示横坐标闭合差0.02045631，按EXE键、显示导线全长相对闭合差5320.360577(K=1/5320)，按EXE键、显示C点纵坐标815.754，按EXE键、显示C点横坐标850. 148，按EXE键依次显示D、E、F、G点的

18、纵横坐标。点号观测角坐标方位角边长纵坐标X横坐标YA1695359.0B 933642801.530724.840 833048.0126.11C1795600815.754850.148 832655.8126.92D1800318830.211976.244 833021.2129.65E1212836844.8501105.067 245904.0126.28F 450648959.2901158.4102500600.0G备注N=5f=37D=508.96fx=+0.093fy=0.020 f=(fx2+fy2)K=1/(d/f)k=1/5350角度闭合差容许值： (图根导线)导线全长

19、相对闭合差：K1/2000程序16渠道中线里程计算文 件 名11渠道中线里程计算164步L1K=0:N=0:AX0CY0:Lbl 1:XY:K=K+1:Pol(XA,YCB= W:Porg A:P=B:Porg CA=X:C=YL2Prog F:D=Ans:K=1J=DK2B=NP:B180B=B360B180B=B+360G=B:JProg CProg G:B=P+3G/4:Prog CB=P+G/4:Prog CN =P:Goto 1文 件 名F取位子程序39步L1U=Frac 10V/10:10U.5V+.1U10U.5VU文 件 名G曲线要素计算子程序101步L1R:V=R tan(A

20、bs G/2:Prog F:T=AnsV=Abs GR/180:Prog F:W =5(Frac VUL2Frac W0V=V+.05Frac W=0V=V.05Prog F:L=AnsV=R/cos(Abs G/2)R:Prog FV=JTV+L/2V+LJ=Ans+DT说明：程序11用以渠道中线测量的里程及其曲线要素的计算程序中：A=渠线起点IP0的纵坐标、C=渠线起点IP0的横坐标、X=IP1IPn的纵坐标、Y=IP1IPn的横坐标、R=半径；程序计算时会自动取舍弧长L的小数位为偶数，其它数据取位后均无余留尾数，用以消除累加误差；程序设计只保留1位小数，若要保留2位小数，应修改取位子程序

21、中的语句U=Frac 10V/10为U=Frac 100V/100；修改取位子程序中的语句:10U.5V+.1U10U.5VU为:100U.5V+.01U100U.5VU；修改曲线要素计算子程序中的语句W=5(Frac VU为W=5(Frac 10V10U；修改语句Frac W0V=V+.05Frac W=0V=V.05为Frac W0V=V+.005Frac W=0V=V.005若要保留3位小数，应修改取位子程序中的语句U=Frac 10V/10为U=Frac 1000V/1000；修改取位子程序中的语句10U.5V+.1U10U.5VU为1000U.5V+.001U1000U.5VU；修改

22、曲线要素计算子程序中的语句W=5(Frac VU为W=5(Frac 100V100U；修改语句Frac W0V=V+.05Frac W=0V=V.05为Frac W0V=V+.0005Frac W=0V=V.0005；算例见下表，按Prog 11调用程序显示X0？输入IP0的纵坐标，按EXE键、显示Y0？输入IP0的横坐标，按EXE键、显示X？输入IP1的纵坐标，按EXE键、显示Y？输入IP1的横坐标，按EXE键、显示IP0至P1的坐标方位角，按EXE键、显示IP0至IP1的距离，按EXE键、显示X？输入IP2的纵坐标，按EXE键、输入IP2的横坐标，按EXE键、显示IP1至IP2的坐标方位角

23、，按EXE键、显示IP1至IP2的距离，按EXE键、显示IP1的里程J，按EXE键、显示IP1的外偏角，按EXE键、显示R？输入半径，按EXE键、显示切线长T，按EXE键、显示弧长L，按EXE键、显示外矢距E，按EXE键、显示直圆点(ZY)里程，按EXE键、显示曲中点(QZ)里程，按EXE键、显示圆直点(YZ)里程，按EXE键、显示直圆点(ZY)至曲中点(QZ)的坐标方位角，按EXE键、显示曲中点(QZ)至圆直点(YZ)的坐标方位角，按EXE键、显示X？依次输入IP3IPn的纵横坐标；计算结束后按AC、RCL、M+键将显示终点IP桩的里程。如果外偏角显示为正值，并且纵坐标X是递增的，则该角位于

24、前进方向的左侧；反之纵坐标X是递减的，遇该角位于前进方向的右侧。如果外偏角显示为负值，并且纵坐标X是递减的，则该角位于前进方向的左侧；反之纵坐标X是递增的，则该角位于前进方向的右侧。IP桩点号XY坐标方位角边长IP桩里程JD外偏角半径R切线长T弧长L外矢距E054110229031585209.0396.67215374023046 396.7475344.46026.650.25.71105824.6128.51525369423166522.2293314.66015.831.02.11403139.1572.591353252235301094.2362922.26019.838.23.

25、21770111.4461.624452791235541554.4IP桩点号直圆点ZY里程曲中点QZ里程圆直点YZ里程直圆点曲中点坐标方位角曲中点圆直点坐标方位角01 370.1 395.2 420.31465342.91225650.72 506.4 521.9 537.41182143.21330820.531074.41093.51112.61493902.21675348.34计算公式：T=Rtg(/2) L=R E=R/(cos(/2)R 校核：q=2TLJD=QZ+q/2 程序17多边形面积计算文 件 名12多边形面积计算100步L1A:XY:M=X:N=Y:C=0:S=0:D=

26、0:B=0:Lbl 1:C=C+1:S=S+3(BX+ DY)/4000:D=X:B=Y:C=AX=M:Y=N:Goto 1CA+1XY:Goto 1S说明：程序12采用多边形顶点坐标计算多边形的面积；程序中A=多边形边数(或顶点数)、X=多边形各顶点的纵坐标、Y=多边形各顶点的横坐标；计算时如果是按顺时针方向输入多边形各顶点的纵横坐标，则面积显示为正值，反之则为负面积是以亩为单位，亩与平方米的换算常数是：2000/3，即 1亩=2000/3M2、1M2=3/2000亩；算例：按Prog 12调用程序显示A？输入多边形的边数(或顶点数)4，按EXE键、显示X？输入顶点1的纵坐标10，按EXE键、显示Y？输入顶点1的横坐标10，按EXE键、显示X？输入顶点2的纵坐标40，按键、显示Y？输入顶点2的横坐标10，按EXE键、显示X？输入顶点3的纵坐标40，按EXE键、显示Y？输入顶点3的横坐标 60，按EXE键、显示X？输入顶点4的纵坐标10，按EXE键、显示Y？输入顶点4的横坐标60，按EXE键、显示四边形的面积2.25(亩)。坐标法计算多边形面积的公式如下：（1）（2）其中：i=1、2、3、4n公式（1）和（2）中：S=多边形面积、n=图形界址点个数、Xi=界址纵坐标、Yi=界址点横坐标；当i=1时，i1=n；i=n时，i+1=1程序18闭合导线近似平差文 件 名13