5800计算程序.docx

1、5800计算程序1 引言2004年初上市的fx-4850P是fx-4800P的小改进产品，主要改进内容有： 内存容量由4500字节扩充到28500字节； 新增六十进制角度显示命令 与 ； 将测量常用的Pol与Rec函数移至键面； 改进赋值显示命令功能，如在fx-4800P中，语句CSLOPE(m)= (A2 +B2) 的功能是将 (A2+B2)的计算结果赋值给C变量，同时显示字符“SLOPE(m)”与变量C的值，而在fx-4850P中，实现上述功能的语句改为C= (A2+B2):SLOPE(m)=:C ，也即，fx-4850P是将赋值与显示功能分开进行。在卡西欧的fx编程计算器系列产品中，fx

2、-4500PA、fx-4800P与fx-4850P统称为fx-4000机型，设计fx-4000机型的初衷是满足最基本的便携编程计算需求。fx-4000机型使用的简BASIC语言只有 、 、 、Goto、Lbl、Dsz、Isz、Pause、Fixm、等11个语句，而 、 、 、 、 、 只是条件符，缺少逻辑运算语句And、Or、Not。从最早进入我国市场的fx-4500PA算起，fx-4000机型在中国已累计销售了10多年，尤其是2003年夏普公司停产PC-E500后，fx-4000机型销量的增长非常显著，其中，2007年的销量达到了9万台，其中大部分用户分布在与测绘、工程建设企业与学校。在卡西

3、欧的fx编程计算器产品系列中，卡西欧图形编程机fx-7400G、fx-9750G、fx-9860G等使用的是类BASIC语言，其程序语言有?、 、If、Then、Else、IfEnd、Lbl、Goto、Dsz、Isz、 、Locate、Cls、And、Or、Not、For、To、Step、Next、While、WhileEnd、Do、LpWhile、Break、Return、Stop、Getkey，条件符仍为 、 、 、 、 、 。显然，类BASIC语言属于结构化程序语言，比简BASIC语言的功能要丰富且强大得多，尤其是IfIfEnd、ForNext、DoLpWhile、WhileWhileE

4、nd控制结构语句与And、Or、Not逻辑运算语句，对简化编程有非常重要的作用。fx-4000机型在我国测绘与工程建设行业应用的日益普及，其功能缺陷带来的问题也逐渐浮出水面： 内存不是闪存，机内程序需备用电池维护，更换备用电池时如果同时取出了主电池就会丢失机内程序。程序一旦丢失就不可恢复，唯一的办法只有重新输入。fx-4800P的内存只有4500字节，重新输入程序的工作量还可以接受，而fx-4850P有28500字节内存，重新输入丢失程序的工作量是非常惊人的。专著1出版后，我们已收到了很多fx-4850P用户的电子邮件，反映在更换 *广东省高等教育教学改革项目(2006036)，“十一五”国家

5、级规划教材建设项目(703)作者简介：覃辉，1962年11月生，男，壮族，广西南宁，教授，硕士，主要从事测量计算、成图与数字化放样研究。备用电池时丢失了机内程序，询问卡西欧工厂是否可以恢复，并且愿意承担一切高额的费用。用户丢失程序的心情是可以理解的，因为很多程序都是用户在工程建设中经过多年实践工作积累的，工程用户在程序调试通过后普遍很少留有程序的文字记录，丢失程序后，又需要重新编写与调试，还不仅仅是简单地照单输入的问题，更多的是需要重新编写； 没有数据通讯功能，只能以手工按键方式输入程序，不便于用户相互交流程序。针对fx-4000机型应用中的问题，卡西欧于2006年10月推出了换代机型fx-5

6、800P，并于2008年3月在全球范围停产了全部fx-4000机型，国内100多万fx-4000机型用户很快都将面临如何将自己的fx-4000机型程序转换为fx-5800P程序的问题，专著2出版后，笔者也收到了很多询问如何转换fx-4000机型程序的电子邮件。本文以fx-4850P为例，用案例方式介绍了fx-4850P与fx-5800P程序语言的主要差异与程序转换方法。2 fx-5800P的内存与数据通讯fx-5800P的内存容量为28500字节，虽然与fx-4850P的内存容量相同，但fx-5800P是用闪存作为内存，就像普通U盘一样，机内的程序与数据可以脱离电源保存。fx-5800P使用一

7、节7号电池供电，没有备用电池，只要用户不手工删除程序，机内程序永远不会丢失，确保了程序的安全。由于采用闪存作内存，fx-5800P还设计有现场保护与恢复功能，在任意模式或状态下，甚至在运行程序的中途，都可以按 键关机，机器自动保存当前模式与状态，再次按 键开机时，机器自动恢复最近一次关机的模式、状态与数据。fx-5800P能实现两台机器相互数据通讯。用SB-62数据线连接两台fx-5800P的3Pin音频通讯口，在接收机上按 (LINK) (Receive)键，使接收计算器处于接受数据状态，在发送计算器上，按 (LINK) (Transmit) (All) 键，发送计算器开始发送机内的全部程序

8、、自定义公式及设置内容。数据传输完成后，发送计算器与接收计算器同时显示“Complete!”，分别在发送计算器与接收计算器上按 键退出数据传输状态。虽然fx-5800P没有卡西欧图形编程机fx-7400G、fx-9750G、fx-9860G与PC机双向数据通讯的功能，但相比fx-4850P还是前进了一步，至少给用户相互交流程序带来了很大的便利。3 程序语言fx-4850P使用简BASIC语言，fx-5800P使用类BASIC语言，与卡西欧图形编程机fx-7400G、fx-9750G、fx-9860G的程序语言完全相同，用户学会了fx-5800P的使用，就很容易触类旁通地学习图形编程机的使用。

9、按 (PROG)键调出程序命令菜单，共3页，内容见图1所示。图1 fx-4850P的程序命令菜单 按 (PROG)键调出程序命令菜单，共5页，内容见图2所示。图2 程序输入与编辑状态下的程序命令菜单其中，代表fx-4850P，代表fx-5800P，下同。(1) 相同功能程序命令的转换1) 赋值语句案例：C= (A2+B2)； (A2+B2) C。2) 变量输入语句案例：A:ASIDE(m)=； SIDE(m)=?A或SIDE(m)=? A。3) 简单条件语句案例：J 0 J=J+360 ； J 0 J+360 J或If J 0:Then J+360 J:IfEnd。4) 复杂条件语句案例：J

10、0 F=J+360: F=J ；If J 0:Then J+360 F:Else J F:IfEnd。5) 逻辑运算案例：A 0 B 0 C= (A+B) ； A 0 And B 0 (A+B) C6) 定义扩充变量案例：Defm n，n=12850，n只能是数字，不能是变量； n DimZ，n=12372，n可以是数字，也可以是变量。7) 计数转移命令：Dsz，Isz；Dsz，Isz，功能完全相同。8) 清屏命令：Cls；Cls，功能完全相同。9) 标记与转移命令：Lbl n，Goto n；Lbl n，Goto n，功能完全相同。10) 暂停命令：Pause n；Locate与Getkey配

11、套使用，功能更强。11) 子程序返回调用主程序命令：无；Return。12) 终止程序命令：无；Stop。(2) fx-5800P的控制结构命令1) ForToStepNext句法1：For To Next控制变量的取值从始值开始，步长为1，重复执行语句块，直至终值为止。一般始值小于终值，如始值大于终值，则不执行语句块而是直接执行Next后的语句。句法2：For ToStep Next它与ForToNext句法的功能相同，唯一区别是增加了步长语句Step。2) DoLpWhile句法：Do LpWhile 先执行语句块，然后测试条件，条件为真时重复执行语句块，否则执行后的语句。无论条件是否为真

12、，语句块至少被执行一次。3) WhileWhileEnd句法：While WhileEnd 先测试条件，当条件为真时执行语句块，然后再测试条件，条件为真时重复执行语句块，直到条件为假时执行WhileEnd后的语句结束循环。熟悉PC机程序语言的读者都知道，控制结构命令是编写大型复杂程序必须具备的，fx-4850P没有控制结构命令，编写大型复杂程序时，只能使用Dsz、Isz或条件语句加Lbl、Goto语句代替，这样编写的程序逻辑较乱，不易读懂。4) Break中断For、Do、While控制结构命令，Break一般位于控制结构命令的条件命令中。4 基本计算1) 省略乘号的运算级别省略乘号的运算级别

13、高于，运算，如执行A=12:B=2:C=3:ABC的结果为2。省略乘号的运算级别等于，运算，如执行12 A:2 B:3 C:ABC的结果为18。2) 带括号函数：sin，cos，Log等均为不带括号函数，sin 30的结果为0.5；sin(，cos(，Log(等均为带括号函数，sin(30)的结果为0.5。3) 清除存储器命令：Scl清除统计存储器，Mcl清除存储器；ClrStat清除统计串列List X、List Y、List Freq，ClrMemory清除存储器，ClrMat清除矩阵存储器，ClrVar清除公式变量。5 矩阵计算fx-4850P没有矩阵功能，fx-5800P的矩阵功能是从

14、图形机fx-9750G与fx-9860G移植过来的，可以在COMP模式或程序中使用矩阵功能。fx-5800P内置了Mat AMat F六个矩阵，最多可以定义10行10列矩阵，可以对矩阵进行加、减、乘、平方、绝对值与转置运算，还可以对方阵求行列式与逆。将矩阵功能应用到测量平差的组成与解算法方程是很方便的。例如，设水准网条件方程为 ，高差观测值的协因数阵为，当系数阵 已输入到矩阵Mat A，闭合差阵 已输入到矩阵Mat B，协因数阵 已输入到矩阵Mat C。由条件平差原理可知，法方程联系数的解为 ，可以在COMP模式执行矩阵表达式(Mat AMat CTrn(Mat A)-1Mat B Mat E

15、求出；观测值改正数的计算公式为 ，可以执行矩阵表达式Mat CTrn(Mat A)Mat E Mat B求出。6 统计计算(1) 样本数据的输入与编辑 用 键输入样本数据，不能查看已输入的样本数据。单变量统计(SD模式)按 或 键删除已输入的样本数据；双变量统计(LR模式)按 或 键删除已输入的样本数据。 单变量统计(SD)模式，在统计串列List X中输入样本数据，List Freq中输入样本数据频率；双变量统计与回归(REG)模式，在统计串列List X与List Y中输入样本数据，List Freq中输入样本数据频率，可对已输入的样本数据与频率进行实时编辑。(2) 回归方程 实际只能进行

16、线性回归( )，对数回归( )、指数回归( )与乘方回归( )需要输入样本变量 或 的函数符。例如，输入对数回归方程样本数据的格式为： ；对指数回归方程两边取自然对数为 ，样本数据输入格式为： ；对乘方回归方程两边取自然对数为 ，样本数据输入格式为： 。 有线性回归(Line)、二次回归(Quad)、对数回归(Log)、e指数回归(eExp)、ab指数回归(abExp)、乘方回归(Power)、逆性回归(Inv)共七种回归类型，无论哪种回归类型，都是将 数据输入到List X串列， 数据输入到List Y串列，频率输入到List Freq串列，回归类型是在事后按 (RESULT) (REG)键

17、，在图3的回归菜单中，按 键选择需要的回归类型计算。图3 fx-5800P回归菜单这种事后选择回归类型的方法也可用于程序中，在程序输入模式，按 (STAT) (REG)键调出与图5类似的回归命令菜单，按 键分别为输入回归命令：LinearReg、QuadReg、LogReg、ExpReg、abExpReg、PowerReg、InvReg。7 fx-4850P与fx-5800P程序转换案例如果只要将fx-4850P程序简单地转换为fx-5800P程序，则只需要转换变量输入语句、赋值语句、条件语句即可。表1程序A-12的功能是根据用户输入的任意直线两端点的平面坐标计算直线的平距与方位角，将其转换为

18、fx-5800P的程序A-13列于表中。表1 坐标反算程序A-12转换案例fx-4850P程序A-12 fx-5800P程序A-13 语句差异说明Deg:Fix 3 Deg:Fix 3 AXO(m)=:BYO(m)= XO(m)=?A:YO(m)=?B ?为fx-5800P的变量输入命令I=0:J=0 fx-5800P只提示输入?后的变量Lbl 0:CD:CX(m)=,0 END Lbl 0:X(m)=,0 END?C fx-5800P用?命令代替了命令C 0 Goto 1 C 0 Goto 1 fx-5800P没有 命令DY(m)= Y(m)=?D ?为fx-5800P的变量输入命令Pol(

19、C-A,D-B):Cls Pol(C-A,D-B):Cls 相同J0 J=J+360 J0 J+360 J 为fx-5800P的赋值命令DIST(m)=:I DIST(m)=:I 相同BEAR(DMS)=:J BEAR(DMS)=:J 基本相同Goto 0 Goto 0 相同Lbl 1:A-12 ENDLbl 1:A-12 END相同由表1可知，简单转换是很容易实现的。28 程序浏览:FileName:RESET 初始化程序Norm 1:50C:12345JRESET PW?I:I=J=500DimZ 为数据库增加额外变量500个，在SET、SETPFDYS程序根据实际再增减变量50ZC+22

20、PASSWORDS?I:IZC+39 重设要素保护密码Cls:StopFileName:DATLOCK 要素保护密码确定认程序Cls:Norm 1:50C PASSWORDS?J:Cls:JZC+39=Stop 2008-08-07修改FileName:SHELL(外壳程序)50C ;在扩充变量预留前50个给别的程序用。如不够就适量加大。RESET，SHELL，SET，SETPFDYS，这几个程序中C值必需一致ZC+35I:STATION-N?I:IZC+35 设置测站N坐标ZC+36I:STATION-E?I:IZC+36 设置测站E坐标ZC+44I:STATION-Z?I:IZC+44设置

21、测站Z坐标ZC+45I:STATION-HI?I:IZC+45 设置仪高ZC+41I:GC-DH?I:IZC+41 放样点高差常数Norm 1ZC+23N当前分段要素N坐标ZC+24E当前分段要素E坐标ZC+25M当前分段要素起点桩号ZC+26H当前分段要素起点方位角（正北）单位：弧度ZC+27A当前分段要素起点曲率有左偏负右偏正（注意不半径）ZC+28R当前分段要素终点曲率有左偏负右偏正（注意不半径）ZC+29L当前分段要素长度ZC+3WZC+32I:DAT1 2 3?I:IZC+32 平曲线要素数库类型选择1为内置式，2文件式，3实时输入(查看当前要素值) 0IZC+4I:1 2 3?I:

22、IZC+4 选择本程序模式默认为坐标正算，1坐标反算，2横向边仰坡放样，3，纵向边仰坡放样(隧道进出口用到)Lbi 1ZC+4II=1=ProgTURNZHI=2=ProgHXBYPI=3=ProgZXBYPLbi ADeg:Norm 1:ClsZC+1G:ZH?G:GZC+1 ;G桩求桩号ZC+2B:JL?B:BZC+2 ;B横向距离，左正右负Lbi BProg ZBJS ;坐标正算Prog GCJS ;高程计算Fix 3:Cls ;设置三位小数ZH:Locate 4,1,G ;第一行显示桩号X:Locate 3,2,X ;第二行显示X(N)坐标Y:Locate 3,3,Y；第三行显示Y(E

23、)坐标Prog GCJS ;高程计算Z:Locate 2,4,Z+ZC+41 ；第四行显示Z坐标Locate 10,4,ZC+2 ；第四行显示横向距离0IDo:I+1I:I=1000=Goto 1:LpWhile Not(Getkey=57 Or Getkey=27) ;锁定键盘，并在几十秒后自动返回Getkey=57= Goto 1Z+ZC+41+ZC+49-ZC+44-ZC+45ZC+48Pol(ZC+35-X,ZC+36-Y)Cls:DL:Locate 3,1,I:Locate 12,1,ZC+45DH:Locate 3,2,ZC+48:Locate 12,2,ZC+49FWJJJ+36

24、0J:JDMS 2008-08-07 修改JZC+46IZC+47Goto 1FileName：ZBJS 坐标计算程序ProgREADDATRadG-MQIF AR=0 :Then If A=R :Then 1J:Else 3J :IfEnd:Else If A=R :Then 2J:Else 3J:IfEnd:IfEndIf J=1 :Then HF:N+QCos(H)X:E+QSin(H)Y:IfEnd 直线段直接计算If J=2 :Then H+QRF:Rec(Abs(2Sin(Abs(0.5QR)R),H+QR2):N+IX:E+JY:IfEnd 圆弧段直接计算If J=3 :Then

25、 Goto 5 :IfEndX+BCos(F+W)XY+BSin(F+W)YDeg:ReturnLbi 5 用五点通用坐标计算计算缓和段0.5(R-A)LKAQIKQ2J0.0469100770P:H+IP+JP2U 0.2307653449P:H+IP+JP2V0.5P :H+IP+JP2D 0.7692346551P:H+IP+JP2F 0.9530899230P:H+IP+JP2T 0.1184634425I0.2393143352J0.2844444444ON+Q(ICos(U)+JCos(V)+OCos(D)+JCos(F)+ICos(T)XE+Q(ISin(U)+JSin(V)+O

26、Sin(D)+JSin(F)+ISin(T)YH+AQ+KQ2FX+BCos(F+W)XY+BSin(F+W)YDeg:ReturnFileName:TURNZH(坐标反算) Norm 1ZC+37I:DQD-N?I:IZC+37 输入待求点N坐标I=-1=ProgINFWJDLDHZC+38I:DQD-E?I:IZC+38 输入待求点E坐标0B：M+L2G：ProgZBJSLbi S:RadZC+37-XI:ZC+38-YJIf I=0 And J=0 :Then Goto A:IfEndPol(I,J)Lbi A:Rec(I,J-F):JZC+2G+IG:IF Abs(I)0.0001 :

27、Then ProgZBJS:Goto S GZC+1FileName：INFWJDLDH实测坐标（用方位角，距离，高差）输入程序Deg:Norm 1:ClsZC+46J:DQD-FWJ?J:JZC+46ZC+47I:DQD-DL?I:IZC+47ZC+48K:DQD-DZ?K:KZC+48ZC+49P:RHT?P:PZC+49Rec(I,J)ZC+35+IZC+37ZC+36+JZC+38ZC+44+ZC+45+K-PZC+43FileName:HXBYP 横向边仰坡放样程序Cls:Norm 1ZC+40I:QPD-DL?I:IZC+40 输入起坡点与中桩距离常数ZC+41I:QPD-DZ?I

28、:IZC+41 输入起坡点与中桩高差常数ZC+42I:i?I:IZC+42 输入边仰坡坡度，左仰坡(路堑)为正，右仰坡(路堑)为负，左边坡为负，右边坡为正Cls:ZC+43I:DQD-Z?I:IZC+43 输入实测高程ProgTURNZHProg GCJS(ZC+43-Z-ZC+41)ZC+42+ZC+40BFix 3:B-ZC+2 显示与上一次测点的偏差BZC+2FileName:ZXBYP 纵向仰坡放样程序(隧道进口使用)Cls:Norm 1:ZC+33I:QPD-ZH?I:IZC+33 输入起坡点桩号ZC+34I:QPD-Z?I:IZC+34 输入起坡点高程ZC+42I:i?I:IZC+42 仰坡时(隧道进出口仰坡)，进口为正，出口为负)Cls:ZC+43I:DQD-Z?I:IZC+43 输入实测高程ProgTURNZHZC+42(ZC+43-ZC+34)+ZC+33GFix 3:G-ZC+1 显示与上一次测点的偏差GZC+1FileName:SET 设置程序1Lbi S2Norm 1350C450ZC+2250K6ZC+5I:1 Or 2?I:IZC+5 ;选择1线或者2线7180ZC+3J:XZJJ?J:J180ZC+3;斜桩夹角8If I=1:Then ZC+8ZC+6:ZC+9