wa2__对卡西欧4800计算器编程

对卡西欧4800计算器在编程中的浅见作者: luzhmu发布日期: 2005-10-08 查看数: 在计算能力飞速发展的今天，各种计算工具应运而生代替了以往我们大量、繁琐的分析计算工作，简化了计算步骤节省了计算时间。在种类繁多的计算工具中适合测量工作的计算器也大有所在。比如有一款SHAPPC-E500计算器，它相当于是基于BASIC语言的微型电脑，功能也相当齐全能编写较复杂的程序、内存也相当大有32KB、64 KB、128 KB、256 KB、几种型号，同时具有输入输出通讯设备。但价格不菲，现在市场价格大约为2000元到3000元，个人经济承受能力有限。通过几年的实际测量工作我个人为认为真正适合测量人员较为专业的计算器应首推基于科学计算器具有简单程序语言的casiox计算器。它有很多优点：一、使用简单，最大的一个优点就是只要具有初中水平会运用数学公式对函数有一般的了解的人，能基本理解casiox的内装函数即可进行简单的编程。非常适合测量初学者和数学程序爱好者学习，也适合专业的测量人员的使用和能力提高。不像其它的专业一定要具有较高的专业技术水平才能进行电脑编程。二、成本低廉，价格在400元左右。三、携带方便，体积很小可随时放在口袋里随拿随用。casio系列较好的编程型号有casio（以下简称）、casio（以下简称），（好象近来还推出了一款casio4850）前者较内存小，只有个字节，能应付一些较为简单的公式计算和科学计算，但由于内存有限，对一些较复杂或子程序过多的程序就力不从心了，不能出色的完成测量任务。就比有较大的改进，内存达到的个字节，而且显示屏是的几倍大，能更准确的显示数据，内装函数字符一目了然。且具有简单的人机对话功能，出现了菜单子菜单。还在增加了啊佛加德罗常数、万有引力长常数、详见操作说明书。有很多很实用的功能如他有一个公式解答功能，其原理是用牛顿法解方程。举个例子：有一方程式：a=2b-c求 当a=2、c=5时的b值。将该方程式存入公式存贮器中，：先按按解答键“SOLVE” a输入2 、c输入5，再按解答键“SOLVE”计算器就会显示：b=3.5。此项功能被称为自动解答功能。同时它也是非常实用的，在实际工作中通常要有一个经常使用的小公式，可借助它来完成。使用的程序语言可以算做简单的BASIC语言，有的命令如GOTO（转移到）、PAUSE（暂停）就与BASIC语言的一模一样。现在的程序语言来说说。 其主要命令有： 1、= 条件转移成立符号，其用法相当于BASIC中的IFTHEN（假设语句相当于假如然后，IF相当于条件THEN相当于结果）语句2、条件转移不成立符号，其用法相当于BASIC中的IFELSE语句通常二者连用，相当于BASIC中的IFTHENELSE语句（它的英语形式一般为if ab then cd else if ba the）3、 条件转移结束符号，与=和配合使用，放在条件语句最后面。4、 LbI标记命令。用于将一段语句作转换标记。后可接字母、数字、符号，但不能超过两个字节，如不能用10的数字作行标，否则会出现出错信息。5、Goto （条件）转移命令。前面可加条件语句，与BASIC中的GOTO作用相同。通常与LbI一起用，如果所转移的行号无效，则会显示：GO ERROR（详见说明书）出错信息 6、 Dsz减量循环命令。可减少未知数的数量。7、 Isz增量循环命令。8、 Pause暂停命令。后可接09之间的整数n，可使某一数据显示n/2秒钟，然后继续运行下面的程序。它被认为是一个语句。9、 Fixm变量锁定命令。该命令能使其所有变量值（AZ）均当成定数处理。当程序运行时，将不需要输入变数（“”内的变数除外），而是将存贮器中原有的数值来完成计算。10、 变量输入命令。只程序在循环使用时经常发生改变的数字，如里程、和宽度。它的输入方式可以使很多种如ABA,BA B都可以。注意“”和“”必需成对输入。否则会出现Syn ERROR（详见说明书）出错信息。11、=、I=C-K:F=0:Goto 1KU+G=W=K-U:Lb1 2:F=NW2G:I=T-W+W540V2+W36V:Goto 1KU+L-G=W=K-U-G:F=N(2W+G):I=T-RsinF-Q+(P-RcosF):Goto 1KU+L=W=U+L-K:D=E+180:H=-H:S=-S:Goto 2Lb1 1:“X=”:A-cos(D-HArgI)AbsI+Scos(D+HF+Z)“Y=”:B-sin(D-HArgI)AbsI+Ssin(D+HF+Z)数据库:SJK(全线交点要素)KHZ点桩号=A=本交点X坐标:B=本交点Y坐标:C=本交点桩号:D=前直线方位角:E=后直线方位角:R=曲线半径:G=缓和曲线长:Goto 0Lbl 0:Prog“ZBJS”程序说明:K? 输入桩号PJU? 输入边桩至中桩的斜距PJIAO? 输入边桩方向与路线前进方向的右夹角(如程序用于4800，显示部分需要更改才能正常：“X=”：更改为X“X”=； “Y=”：更改为Y“Y”=程序名：ZDZH(主点桩号计算)Q=C“LS”/2-C3/240R2:P=C2/24R-C4/2384R3:T=(R+P)tan(AbsJ“PJ”/2)+Q:L=RAbsJ/180+C:B=Z“JDZH”-T:“ZH=”:B“HY=”:B+C“QZ=”:B+L/2“YH=”:B+L-C“HZ=”:B+L程序说明:LS? 输入缓和曲线长R? 输入半径PJ? 偏角JDZH? 交点桩号程序名：WJPQXA “PJ”:E:L“LS”:N=1cos():M=E(N-1):W=NL2(24N-24):Q=M2-4W:Q0“R=”:0.5(M+Q)程序说明:本程序适用于外距控制单交点平曲线PJ? 输入路线转角E? 输入外距LS? 输入缓和曲线长2006-03-2817:52:15大中小程序名:LZXJD(两直线交点)AX1:BY1:CX2:DY2:EX3:FY3:GX4:HY4:N=C-A:O=D-B:P=G-E:Q=H-F:K=QN-OP:U=NB-OA:V=PF-QE:X=:(PU-NV)/KY=:(QU-OV)/K程序说明:X1? 输入第一点X坐标Y1? 输入第一点Y坐标X2? 输入第二点X坐标Y2? 输入第二点Y坐标X3? 输入第三点X坐标Y3? 输入第三点Y坐标X4? 输入第四点X坐标Y4? 输入第四点Y坐casio4850全线坐标计算程序一、LXZB（坐标计算主程序）LbIZ：PA=：V=：C=：D=：P直线的终点桩号Prog“ZX.SUB”：GotoZ（以上各字母含义详直线坐标计算子程序说明） A=：V=：C=：D=：R=：I=（1或-1）P圆曲线终点桩号Prog“YQX.SUB”：GotoZ（以上各字母含义详圆曲线坐标计算子程序说明）A=：V=：C=：D=：B=：W=：E=：F=：R=：L=：I=（1或-1）PHY点桩号Prog“HQ1.SUB”：GotoZPYH点桩号Prog“HYQ.SUB”：GotoZPHZ点桩号Prog“HQ2.SUB”：GotoZ（以上各字母含义详（HQ1.SUB、HYQ.SUB、HQ2.SUB）坐标计算子程序说明）该主程序是简单的示例说明，具体的应根据图纸的线型组合，随机应变的做顺序调整，并将其输入Casio4850计算器。二、ZX.SUB（直线坐标计算子程序）G=P-A：T=V“X=”：X=C+GcosT“Y=”：Y=D+GsinTQ：Q=1=Prog“BZ.SUB”说明：G：所求桩号至起算点的距离P：所求点的桩号A：起算点的桩号T：所求桩号的方位角V：起始方位角C、D：起算点的X、Y坐标X、Y：显示所求桩号的坐标Q：参数，为1时计算边标坐标，为0时不计算三、YQX.SUB（圆曲线坐标计算子程序）G=P-AH=Rsin(GR)(180)J=(R(1-cos(GR)(180)IT=V+(GR)(180)I“X=”:X=C+HcosV-JsinV“Y=”:Y=D+HsinV+JcosVQ:Q=1=Prog“BZ.SUB”说明：G：所求桩号至起算点的距离P：所求点的桩号A：起算点的桩号(指ZY点)R：圆曲线半径T：所求桩号的方位角V：起始方位角(指ZY点的方位角)C、D：起算点（ZY点）的X、Y坐标X、Y：显示所求桩号的坐标I：参数，（左偏为1，右偏为1）Q：参数，为1时计算边标坐标，为0时不计算四、HQ1.SUB（第一缓和曲线坐标计算子程序）G=P-AH=G-G5(40R2L2)+G9（3456R4L4）-G13（599040R6L6）J=(G3(6RL)-G7(336R3L3)+G11(42240R5L5)IT=V+90G2(RL)I“X=”:X=C+HcosV-JsinV“Y=”:Y=D+HsinV+JcosVQ:Q=1=Prog“BZ.SUB”说明：G：所求桩号至起算点的距离P：所求点的桩号A：起算点的桩号(指ZH点)R：曲线半径L：缓和曲线长LsT：所求桩号的方位角V：起始方位角(指ZH点的方位角)C、D：起算点（ZH点）的X、Y坐标X、Y：显示所求桩号的坐标I：参数，（左偏为1，右偏为1）Q：参数，为1时计算边标坐标，为0时不计算五、HYQ.SUB（带缓和曲线的圆曲线坐标计算子程序）G=P-AH=Rsin(L2R)180+(G-L)R180)+L2-L3(240R2)J=(R(1-cos(L2R)180+(G-L)R180)+L2(24R)-L4(2688R3)IT=V+90L(R)I+(180(G-L)(R)I“X=”:X=C+HcosV-JsinV“Y=”:Y=D+HsinV+JcosVQ:Q=1=Prog“BZ.SUB”说明：G：所求桩号至起算点的距离P：所求点的桩号A：起算点的桩号(指ZH点)R：曲线半径L：缓和曲线长LsT：所求桩号的方位角V：起始方位角(指路线前进方向ZH点的方位角)C、D：起算点（ZH点）的X、Y坐标X、Y：显示所求桩号的坐标I：参数，（左偏为1，右偏为1）Q：参数，为1时计算边标坐标，为0时不计算六、HQ2.SUB（第二缓和曲线坐标计算子程序）G=B-PH=G-G5(40R2L2)+G9（3456R4L4）-G13（599040R6L6）J=-(G3(6RL)-G7(336R3L3)+G11(42240R5L5)IT=W-90G2(RL)I“X=”:X=E-HcosW+JsinW“Y=”:Y=F-HsinW-JcosWQ:Q=1=Prog“BZ.SUB”说明：G：所求桩号至起算点的距离P：所求点的桩号B：起算点的桩号(指HZ点桩号)R：曲线半径L：缓和曲线长LsT：所求桩号的方位角W：起始方位角(指路线前进方向HZ点的方位角)E、F：起算点（HZ点）的X、Y坐标X、Y：显示所求桩号的坐标I：参数，（左偏为1，右偏为1）Q：参数，为1时计算边标坐标，为0时不计算七、BZ.SUB（边桩坐标计算子程序）O，S：O“ZBK”：S“YBK”：“XZ=”：M=X-Ocos(T+90)“YZ=”：N=Y-Osin(T+90)“XY=”：K=X+Scos(T+90)“YY=”：U=Y+Ssin(T+90)说明：ZBK：中桩至左边桩的距离YBK：中桩至右边桩的距离XZ、YZ：显示所计算左边桩的X、Y坐标XY、YY：显示所计算右边桩的X、Y坐标线路中边桩任意里程坐标正反算改编程序CASIO4850 Posted by 忘不了糊涂2006-10-5 23:45:00含数据库、测站放样计算（第五次修改-正式版）1.加编数据库及计算总调度程序,计算中不必逐项输入线元要素,提高运算速度，避免现场忙中出错2.可加入多条线路的数据库,内业输入数据,外业一目了然3.计算直观,人性化4.正算直接输入里程和边距,反算输入近似里程便可5.增加了“计算点与测站点”的距离和方位角计算语句，方便直接放样6.愿收获与大家共享,同时也希望大家提出心得和宝贵建议7.核心计算程序摘自“yshf” 一.改动后的程序清单，增设数据库程序,可输入n条线路的数据库 1.(QXZDJS 计算总调度程序)M=1=Prog SJK1:Prog SUBSJK.M=n=Prog SJKn:Prog SUBSJK2.(SJK1 数据库程序)1.SZ = XY：2.XY = SZ：N,S：N：SDKIS本线元终点里程U=本线元起点X坐标：V=本线元起点Y坐标：G=本线元起点正切线方位角：P=本线元起点曲率半径：R=本线元终点曲率半径：O=本线元起点里程：H=本线元终点里程：Q=本线元左直右偏向(0或1、-1):Goto0S下一线元终点里程=O=。： Goto0LB1 0 3.(SUBSJK 运算主程序)UVGPROHQ：UXA：VYA：GFA：PRA：RRB：ODKA：HDKB：QDeg:SH=Goto AD=(P-R)2PR(H-O)：N=1=Goto 1：Goto 2Lbl 1：Z：Z：W=Abs(S-O)：Prog SUB1：XS=：XYS=：YQX-FWJ=：F=F-90DMSI=0：J=0：Pol（X-CXC，Y-EYC）：J0J=J+360：“FWJ=”:JDMS“I=”：IGoto ALbl 2：XY：I=X：J=Y：Prog SUB2：S=:S=O+WZ=:ZLbl A 4. (SUB1 正算子程序) A=0.1739274226：B=0.3260725774：K=0.0694318442：L=0.3300094782：X=U+W(Acos(G+180QKW(1P+KWD)+Bcos(G+180QLW(1P+LWD)+Bcos(G+180(1-L)QW(1P+(1-L)WD)+Acos(G+180(1-K)QW(1P+(1-K)WD)：Y=V+W(Asin(G+180QKW(1P+KWD)+Bsin(G+180QLW(1P+LWD)+Bsin(G+180(1-L)QW(1P+(1-L)WD)+Asin(G+180(1-K)QW(1P+(1-K)WD)：F=G+180QW(1P+WD)+90：X=X+ZcosF：Y=Y+ZsinF 5.(SUB2 反算子程序)T=G-90：W=Abs(Y-V)cosT-(X-U)sinT)：Z=0：Lbl 0：Prog SUB1：L=T+180QW(1P+WD)：Z=(J-Y)cosL-(I-X)sinL：AbsZGoto1：W=W+Z：Goto 0Lbl 1：Z=0：Prog SUB1：Z=(J-Y)sinF二. 使用说明： 1.把所有相关的”线元要素“依次输入”SJK“，如果有多条线路，可把“SJK”分成SJK1.SJKn等n条线路数据库，然后对应“计算总调度程序”中的M1.Mn备用 2.运算时直接调用“SJK”运行，程序提示输入里程“S”？时，正算直接输入待求点里程，反算输入所求点“近似“里程 3.如果在正算中途想进行反算，重新给”N“？赋值2便可 4.程序中“XC、YC”为测站坐标，“I”、“FWJ”为放样距离和方位角5.程序中“QX-FWJ=”为计算点在中线上的正切线方位角，此方位角在新线路输入数据库时非常重要，我们在每输入完一段曲线元要素后，立即退出返回到“计算总调度程序(SUBSJK)”中，然后输入该段曲线元的终点里程，其计算出来的X，Y中桩坐标应该与设计图纸(当然是经过复查无误的图纸)上的该点中桩坐标一致，否则一定是哪个地方数据输入或者程序有误，这时候我们一定要把该计算点在中线上的正切线方位角“QX-FWJ=”记录下来为下一段曲线元要素数据输入服务，因为它就是下一段曲线元的起点正切线方位角，嘿嘿！ 6.其它输入与yshf原程序一致- 附件:(原yshf程序)曲线任意里程中边桩坐标正反算(CASIO fx-4800P计算器)程序 一、程序功能 本程序由一个主程序(TYQXJS)和两个子程正算子程序(SUB1)、反算子程序(SUB2)序构成，可以根据曲线段直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。另外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲线任意里程中边桩坐标进行正反算。本程序也可以在CASIO fx-4500P计算器及 CASIO fx-4850P计算器上运行。 特别申明：(1). 适用于弧长小于2倍半径的各种线元坐标正反算，精度优于1mm； (2). 在引用该核心计算部分时，请注明来源。 二、源程序1.主程序(TYQXJS)1.SZ = XY：2.XY = SZ：N：UX0：VY0：OS0：GF0：HLS：PR0：RRN：Q：C=1P：D=(P-R)(2HPR)：E=180：N=1=Goto 1：Goto 2Lbl 1：SZ：SZ：W=Abs(S-O)：Prog SUB1：XXS=XYYS=YFFS=F-90Goto 1Lbl 2：XY：XY：I=X：J=Y：Prog SUB2：SS=O+WZZ=ZGoto 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)+Bcos(G+QEFW(C+FWD)+Acos(G+QEMW(C+MWD)：Y=V+W(Asin(G+QEKW(C+KWD)+Bsin(G+QELW(C+LWD)+Bsin(G+QEFW(C+FWD)+Asin(G+QEMW(C+MWD)：F=G+