CASIOFX4850公路计算程序.doc

0前言：传统公路测量中，使用的仪器设备和方法都很落后，需带着数学用表、曲线用表、计算盘、计算尺和算盘等一类的工具，完成外业测量工作。计算器的出现，改变了这一局面。高速公路建设中，长大曲线比比皆是，传统中对公路中线的测设方法，被极坐标法彻底的否定与取代，但大量的计算工作，只能带着提前计算好的线路逐桩坐标、高程资料，进行外业测量工作，机动性很差，现场查找也不方便。这些问题都能在CASIO系列可编程计算器上得到很好的解决，对CASIO系列可编程计算器如何使用，直接影响到测量成果的质量和工作效率，本文将对CASIO系列可编程计算器快捷的计算方法进行分析与介绍。 1：以知线外任意点坐标，求对应线路里程在缓和曲线上，要计算任意里程的法线方向及任意宽度的边线坐标，非常简单。但要计算任意一个已知坐标点，是对应哪一个里程法线方向上的点，就有一些困难。很难推导一个这样的计算公式。唯一的方法“渐进”，如果手工计算这可不是一个好方法。但在有CASIO系列可编程计算器，如：FX-4500的情况下就变的非常简单了。亦可用于直线和圆曲线的计算。首先在缓和曲线上任选一点A为起始点，计算该点的坐标和切线方位角，通过坐标反算求起始点A与计算点B的方位角和距离，B点肯定对应A点切线方向上有一个垂足C点，把三点看成一个直角三角形，通过解直角三角形计算AC的距离，当该距离大于某一数值，如0。001m，A点里程加AC的距离等于C点的里程，回到开始重新进入新一轮的计算，如果AC的距离小于某一规定值，则计算C点的里程与BC的距离即可。求对应线路里程程序：主程序QLC(已知坐标求里程)Lb10：LE：ProgXH：Goto0子程序：XH(循环)L1Lb11L2Norm：ProgLYYD：L3PO1(D-X,E-Y)：W0=W=W+360L4Z=W-I：A=VcosZ：L=L+AL5AbsA0.001=Goto1：=B=VsinZ：Fix3：“FXJL=”L6L:Fix3:“DYLC=”程序中字母代表D任意点X坐标，E任意点Y坐标，DYLC对应里程，FXJL中线法线距离。程序中有坐标反算功能。使用方法：只需输入计算点坐标、和较为接近的桩号。桩号越接近计算速度越快2：逐桩坐标计算2.1编制方法：线路坐标程序是按照平曲线为单元，直线部分归属在曲线两端的方法，把整段路线分段装进数据库，根据桩号判断采用数据通过共用程序，进行任意点的坐标计算，在坐标转换示意土，第一直线段，是通过方位角和距离直接计算大地坐标，第一缓和曲线和圆曲线段，是先计算任意点切线支距和方位角然后转换大地坐标，第二缓和曲线段和直线段是先计算任意点切线支距和方位角。然后转换为ZH坐标系的坐标，通过ZH坐标系的坐标再转换为大地坐标。2.2使用方法2.2.1准备工作：室内把已知曲线条件，装进数据库，曲线划分界线、判断条件装进子程序LYYD(路由引导)。2.2.2现场使用：根据计算机提示输入相关数据即可。提示情况如下：K公里桩号如312，启动程序出现一次。L细部里程桩号如518.如采用渐进只出现一次，否则逐桩输入。过千米桩时需输入1000确认。O渐进长度，如20米一点,取O=20，公里桩号也自动渐进。否则O=0，启动程序出现一次。Y断链条件，执行输0不执行默认，不输入【】以内的程序，Y不出现。E边线角度，法线为90度，分正负值，输E=0此后则不在出现计算边线的过程。D边线点至中线点的距离VW输出的边线1的大地坐标XY输出的边线2和中线的大地坐标,2.3逐桩坐标计算程序主程序：XLZB(线路坐标)L1Lbl0：L1000=P=P+1：L=L-1000L2O=0=progFJJL=L+O：progLYYD：progXSZB：E0=progBXGoto0子程序：FJJ(非渐进)L：L=L子程序：LYYD(路由引导)N=（P“K”+L/1000）1000：【X=0：Y=6】N*.*=prog1：progPQX：progZJN*.*=prog2：progPQX：progYJN*.*=prog3：progPQX：progZJN*.*=prog4：progPQX：progYJN*.*=prog5：progPQX：progYJN*progZB子程序：PQX（单圆曲线及带缓和曲线的平曲线）L1B=Z+Q-S：H=Z+Q：A=S/2-S3/240R2：T=A+（R+S2/24R）tan（F/2）L2NZ=V=N-Z：W=0：I=0：Goto1L3U=N-Z：N(Z+S)=V=U-U5/40R2S2：W=U7/336R3S3-U3/6RS：I=90U2/RS：Goto1【L*Y：Y=0=U=U+X】L4NB=I=90（2U-S）/R：V=RsinI+A：W=RcosI-R-S2/24R：Goto1L5NH=U=H-N：C=U-U5/40R2S2：G=U3/6RS-U7/336R3S3：V=（T-C）cosF-GsinF+T：W=（C-T）sinF-GcosF：I=F-90U2/RS：Goto1L6NH=V=T+（T+N-H）cosF：W=（H-T-N）sinF：I=F：Goto1L7Lbl1子程序：YJ（右角）W=-W：I=K+I子程序：ZJ（左角）I=K-I子程序：ZB（坐标）X=J+VcosK-WsinK:Y=M+VsinK+WcosK子程序：XSZB（显示坐标）O0=L=L：pause5X=XY=Y子程序：BX（边线）DE：I=I+E：V=X+DcosIW=Y+DsinIDE：I=I+E：X=V+DcosIY=W+DsinI2.4数据库：(每一组曲线占用一个子程序)1K=*：F=*：R=*：J=*：M=*：Z=*：Q=*：S=*：22K=*：F=*：R=*：J=*：M=*：Z=*：Q=*：S=*：【N*.*=X=*.*】2.5注解：程序XLZB：线路坐标，它是计算逐桩坐标的主程序。程序：FJJ(非渐进)O0时，只需输入起始点桩号如计算为每20米一点时,取O=20，此后则自动渐进，公里桩号也自动渐进。起始桩号应输入第一个计算点桩号减渐进长度。如K36+700输入36+680即可。否则取O=0。每一个点均需输入细部点桩号。当公里桩号发生变化时，如：计算K25+910-K26+110，每20米一点。K25+990完了便是K26+010，此时无需重新输入公里桩号，只需输入一个大于等于1000的桩号，此后则按正常方法输入。例：K25+970K25+990K25+1000K26+010K26+050程序LYYD：路由引导，段数根据曲线数量确定增减。L是曲线间的分界点桩号。用选ZH点或ZY点以前的桩号。用选HZ点或YZ点以后的桩号。为了做为QLC(已知坐标求里程)的子程序，故于XLZB：(线路坐标)分为两个程序。否则可和二为一。P显示K，为公里桩号，为输入方便，可省略公里桩号中的相同之处，如K315+200-K395+800，输入公里桩号时，可省略百位的3，只输十位和个位的15-95即可。L为细部桩号，如+660.318，程序PQX：（单圆曲线及带缓和曲线的平曲线）第一行，计算曲线要素，第二行，计算第一直线段任意点坐标，第三行，计算第一缓和曲线上任意点坐标，L*断链；建议不采用第四行，计算圆曲线上任意点坐标。第五行，计算第二缓和曲线上任意点坐标，第六行，计算第二直线段任意点坐标。程序J：曲线偏角为右角时，进入该程序W=-W，转换偏角F为左角，执行程序后，ZH坐标系统的坐标（V，W）转换成大地坐标（X，Y），I为曲线上任意点到ZH坐标系统中X轴的夹角，K+I是该点切线沿线路前进方向的方位角。程序J：曲线偏角为左角时。进入该程序，其它意义同上。程序：是坐标转换程序。计算线路坐标时不显示，以程序XSZB：显示计算结果。程序XSZB：O=0时显示坐标，O0时显示桩号和坐标，为了做为QLC(已知坐标求里程)的子程序，故于ZB：(坐标转换)分为两个程序。否则可以取消。程序：是求线路外任意点的坐标（V，W）其中E为夹角，有正负之分，顺时针为正，逆时针为负，线路法线为正、负90度。D为线路中线点到计算点间的距离，(X，Y)线路外第二任意点的坐标，如斜交桥、涵的坐标计算，（V，W）为涵口边墙或桥台坐标，(X，Y)为八字墙端部坐标。输E=0计算边线的过程此后则不出现。，程序：数据库程序，用数字、表示，根据曲线数量确定增减，其中Q：曲线总长；F：偏角；R：半径；S：缓和曲线长，在单圆曲线中输0；Z：ZH点里程；（J，M）ZH点大地坐标；K：ZH点至JD点的起始方位角。N*=X=*计算点桩号和断链长度，(下文详述)程序中部分字符以标出，有些字符在不同位置意义不同，循环使用，不宜标出。只要把需输入和输出的字符搞对即可。2.6单圆曲线的平曲线是通过ZY点坐标计算圆心的坐标，通过圆心的坐标计算曲线上任意点坐标（X，Y），B：中线到边线的距离分正、负值。除S转向角为左角输+1转向角为右角输-1以外。数据库和PQX（单圆曲线及带缓和曲线的平曲线）基本相同，它可代替除数据库和路径引导程序外的所有子程序，未编第二直线段部分，因为本曲线第二直线段部分也就是下一个曲线的第一直线段部分。未考虑与其它程序的配合和对断链的处理。可根据所管工程线型情况选用该程序。主程序：XLZB(线路坐标)Lbl0：L：L=L+O：N=（P“K”+L/1000）1000：N*.*=prog1：progPQXY：Goto0子程序：PQXY(单圆曲线的平曲线)L1U=N-Z：NZ=X=J+UcosKY=M+UsinKX=X+Bcos(K+90)Y=Y+Bsin(K+90)Goto1L2=V=J+Rcos(K+90S)：W=M+Rsin(K+90S)：E=180U/RL3I=K-90S：I0=I=I+360I=I+SEL4B：X=V+(R+B)cosIY=W+(R+B)sinIGoto1L5Lbl13：逐桩高程计算3.1编制方法：纵断高程程序是按照竖曲线为单元，同坡部分归属在曲线两端，把整段路线分段装进数据库，根据桩号判断采用数据通过共用程序，进行任意点的高程计算，3.2使用方法3.2.1准备工作：室内把已知曲线条件，装进数据库，曲线划分界线、判断条件装进子程序LJYD(路径引导)。3.2.2现场使用：根据计算机提示输入相关数据即可。提示与输入情况如下：提示K、L、O、Y同逐桩坐标计算程序，Z输出高程3.3：逐桩高程计算程序主程序ZDGC(纵断高程)Lbl0：L：L=L+C“O”：N=（K+L/1000）1000：progLJYD：Goto0子程序：LJYD(路径引导)L1【V=0：X=0：Y=9】L2N*.*=progA：N*.*=progB：N*.*=progC：N*.*=progD：N*.*=progE：NProgSQX子程序SQX(竖曲线)L1【Y：Y=0=V=X】M=A-T【-V】：W=A+T【+V】：U=Abs(A-N)【-V】：NM=Z=H-JUGoto1NA=Z=H-JU+F(N-M)2/2RGoto1NW=Z=H+IU+F(W-N)2/2RGoto1NW=Z=H+IULb113.4数据库：(每一组曲线占用一个子程序)A：R=：T=：A=：H=：J=-：I=-：F=1：【N*=X=-*】B：R=：T=：A=：H=：J=-：I=-：F=-13.5注解：主程序ZDGC(纵断高程)为了于三维坐标段落法隧道断面测量程序配合，于LJYD(路径引导)一分为二，否则可合二为一。程序：LJYD(路径引导)是路径引导程序，段数根据曲线数量确定增减。L是曲线间的分界点桩号。用选曲线起点以前的桩号。用选曲线终点以后的桩号。程序SQX(竖曲线)第一段计算曲线起点以前的高程，第二段计算曲线起点以后的高程，第三段计算曲线终点以前的高程，第四段计算曲线终点以后的高程程序：数据库程序用字母A、B、C表示，根据曲线数量确定增减。程序中字母代表R表示竖曲线半径，T表示切线长，A表示变坡点里程，H表示变坡点高程，F=-1表示凸曲线，F=1表示凹曲线，J表示前一竖曲线坡度，下坡为负，上坡为正。I表示后一竖曲线坡度，下坡为负，上坡为正X断链长度，分正、负值，。4.对断链的处理方法：高速公路中坐标法控制线路的平面位置，断链较少。设计上以考虑到施工计算方便的问题。平曲线内一般不会出现断链，尽可能也不设在竖曲线内，一般会将断链推到直线同坡段。直线部分归属和划分，应考虑到断链，依断链桩号为划分界线。有时能躲的开平曲线但躲不开竖曲线，程序ZDGC(纵断高程)以考虑到竖曲线内出现断链的情况，处理方法是：数据库中赋值，当计算点大于或小于某一桩号时，修正计算点到切点和变坡点的曲线长度，程序如下：N*.*=X=*.*，N：为计算点桩号，*为断链点桩号，X：为断链距离。分正、负值。无断链时，程序自动赋值X=0，其中：N*=X=*。变坡点以前出现长链，如：K*+530=K*+480N*+530=X=50变坡点以前出现短链，如：K*+480=K*+530N*+530=X=-50变坡点以后出现长链，如：K*+630=K*+580N*+580=X=-50变坡点以后出现短链，如：K*+580=K*+630N*+580=X=50当遇短链如：K*+480=K*+530，+480至+530之间没有距离，输Y=0即可，当遇长链如：K*+530=K*+480，+480至+530之间有二倍的距离，有两个完全一样的里程，输Y=0只算了后一个+480至+530，若计算前一个+480至+530，取Y为任意值。偶遇平曲线内出现断链，要有就在圆曲线上，肯定不会在缓和曲线上。处理方法类似竖曲线。程序中【】符号并非计算机运算符，没有断链时，【】符号内的程序不输入计算机。线路中断链不多的情况下，为了提高运算速度，建议不采用【】以内的程序处理断链。对断链进行单独的处理。5坐标反算主程序：ZBFSL1Lb10：DE：Norm：PO1(D-X,E-Y)：Fix3：“S=”L2W0=W=W+360IntW+Int(fracW60)/100+frac(fracW60)0.006：Fix4：“AV=”Goto0程序中字母代表D任意点X坐标E任意点Y坐标AV输出角度S输出距离坐标反算输出角度小数点后四位为分和秒，如：168.3639为168度36分39秒。6结语公路施工测量工作，全站仪完全满足了极坐标法放样的硬件要求，CASIO系列可编程计算器完善了全站仪在公路测量中的软件不足之处，珠联璧合。使的极坐标法在公路测量中得到了良好应用。极坐标法放样和可编程计算器改变了施工测量中的放样模式，解决了很多过去不好解决的问题，对可编程计算器如何使用，直接影响到测量成果的质量和工作效率，对可编程计算器充分利用，公路外业测量工作不需要再带线路逐桩坐标、高程资料，只带一台CASIO系列可编程计算器即可。外业测量工作中，只需输入里程，即可提供线路任意点坐标、高程。不但方便而且及时准确。相当于把线路平面和纵断面装进了计算机。并解决了缓和曲线段以知线外任意点坐标，求对应线路里程的难点问题，亦可用于直线和圆曲线的计算。建议的程序排列顺序，主程序：1：线路坐标2：纵断高程、3：求里程4：坐标反算5隧道断面6-10预留空位或其它程序，子程序：11-17线路坐标的各子程序，18-19纵断高程的两个子程序。19以后为线路坐标和纵断高程的数据库。建议使用CASIO系列FX-4800或FX-4850大容量机型，FX-4500装不下常用的全部程序，即便装一部分或单个程序，数据库也无足够的空间装载线路数据.正如一句话的描写，4500反应不快，4800即将淘汰，4850大容量风行时代。CASIOFX-4850有28K字节的容量，以上所有程序不过1200字节。所以上百公里的平、纵断面仅需一台CASIOFX-4850即可。CASIO fx4850万能坐标计算程序一、程序内容程序：McW“1.js 2.SZ”:W=1=Goto 1 W=2= O“KOU LING”:O123456=O=0“OUT” Goto 4 O=0:Defm 83:Z81=1:Goto 0Lbi 0ABCDEFGU:A“X0”:B“Y0”:C“F0:D“R0”:E“RN:F“D0”:G“LS”:U“G” ZZ81=A:ZZ81+10=B:ZZ81+20=C:ZZ81+30=1D:ZZ81+40=1E:ZZ81+50=F: ZZ81+60=F+G: ZZ81+70=UA=0= Dsz Z81:Goto 1Isz Z81:Goto 0Lbi 1HS:H“D:S“Z”:Z82=1:Goto 2Lbi 2Z81Z82 =”OUT” :Goto5HZZ82+60=A=ZZ82:B=ZZ82+10:ZZ82+20=C:D=ZZ82+30: E=ZZ82+40: F=ZZ81+50: G=ZZ81+60: U=ZZ81+70: Goto3:Isz Z82:Goto 2Lbi 3PU（E-D）Abs（G-F）：QAbs（H-F）：IPQ：JC+90 Q（I+2UD）/： JJ=J+360MC+45 Q（I4+2UD）2：NC+135Q（3I4+2UD）2：KC+45Q（I2+ 2UD） S=0=Z=0：Goto4S0=Z1：Goto4Z1：Goto4Lbi 4XA+Q(Cos C+4（Cos M+Cos N）+2Cos K+Cos J) 12 +ZSCos（J+90Z) YB+Q(Sin C+4（Sin M+Sin N）+2Sin K+Sin J) 12 +ZSSin（J+90Z) S=0=“X”:X:Pause 0: “Y” :Y Goto 1S“XL”:X:Pause 0: “YL”:Y Goto 1S0=“XR”:X:Pause 0: “YR”:Y Goto 1 Lbi5 二、 变量及说明X0：Y0：F0曲线元起点X、Y坐标及起点正切线方位角R0：RN曲线元起点及终点半径D0：LS：Q曲线元起点桩号、路线长度及线路左右偏标志（左=-1,直线=0,右=1）D曲线元中待求点桩号Z计算边桩距中线平距，左边输入负值，右边输入正值，中桩输入0X：YXL:YLXR:YR 待求点中边桩桩号的X，Y坐标扩充变量：ZZ81: ZZ81+10: ZZ81+20: ZZ81+30: ZZ81+40: ZZ81+50: ZZ81+60: ZZ81+70: 分别为各线元X0:Y0:F0;R0:RN:D0:LS:Q (1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。 (2) 当所求点位于中线时，Z=0，坐标显示X Y；当位于中线左侧时，Z取负值,坐标显示XL YL,；当位于中线右侧时，Z取正值,坐标显示XR YR。 (3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。 (4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。 (5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。 (6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。本程序可以根据曲线段直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距，对该曲线段范围内任意里程中边桩坐标进行正算。本程序是对扩充变量的应用，实现了真正意义上的的全线贯通。程序分为两部分：1为计算2为设置。首先对内存变量进行扩充，增加了83个扩充变量，为Z1Z83，为防止误操作在输入正确口令123456时方可进行设置。本程序设计为十段线元要素，如有需要可以根据内存容量进行增加。在执行程序中将各个线元要素按规律输入到扩充内存变量中，在坐标计算时首先判断其在哪一线元内，并把其线元要素调出进行计算。在进行计算时直接输入桩号、宽度就可以进行全线坐标计算。 本程序也可在4800上运行，只需将输出部分进行修改。四、 示例郑石高速路面9标主线线元要素及相关数据如下：X0 Y0 F0 R0 RN D0 LS Q744450.244 510454.431 22041 31.8 7000 7000 116075.889 3386.293 1742505.351 507722.676 24824 33.6 1045 1045 119462.182 1980.825 0741776.461 505880.832 24824 33.6 8500 8500 121443.007 1517.999 1741346.546 504427.086 25838 30.1 1045 1045 122961.007 1985.977 0桩号 宽度 X Y 117000 0 743791.308 509807.481118000 -20 743164.306 509027.467119000 20 742707.550 508137.933119462.182 0 742505.349 507722.673120000 -15.253 742293.266 507228.206121000 15.260 741953.665 506287.142121443.007 0 741776.461 505880.833122000 -10.20 741578.902 505359.701122961.007 50.355 741395.915 504417.169CASIO fx - 4850 最豪华的三坐标计算 程序本程序可同时录入两路的全线原始平曲线要素和坚曲线要素，计算任意桩号点中桩，法向：斜向桩三个坐标值，并且有多断链处理功能由于本程序使用阵列变量编写，所以只要你的内存足够不管你线路有多少交点，都能一次录入（并且只需录入一次）永久使用，不用每次输入参数缓和曲线可以输入不对称值技术支持QQ:285507730主程序 文件名:MAILLbl Arad:NormA=Z1:A:AZK1 YK-1 IN0A=0= Prog IN0:Goto AZ1=AA=Z11:A:AJJ:Z11=ALbl 1NormG=Z9+Z14:G:GDKProg DLG=-1=Goto 2A=Z10:A:AJL:Z10=AGZ9=Prog COREGGoto 2Rec(Z10,W+/2+Z11/180)Fix 3X,Y,Z=X+J:Pause 0Y+I:Pause 0ZGoto 1Lbl 2Norm:DegEXITA=1:A:AYES 1 NO -1AGoto A子程序 文件名：CORELbl AZ1=1=Goto 1Z1=-1=Goto 2Lbl 1D=20:E=0:F=Z3-2:V=Z7:I=0:J=0G=-2=F=Z3: Prog IN1:Goto BG=-3=D=20+5Z3+5Z4:E=0:F=Z5: Prog IN2: Goto BG=-4=D=20+5Z3+5Z4+4Z5+4Z6:E=0:F=Z12: Prog IN3: Goto BZ9=GProg ZBD=20+5Z3+5Z4:E=0:F=Z5-2Prog GCGoto BLbl 2D=20+5Z3:E=0:F=Z4-2:V=Z8:I=0:J=0G=-2=F=Z4: Prog IN1: Goto BG=-3=D=20+5Z3+5Z4+4Z5:E=0:F=Z6: Prog IN2: Goto BG=-4=D=20+5Z3+5Z4+4Z5+4Z6+2Z12:E=0:F=Z13: Prog IN3: Goto BZ9=GProg ZBD=20+5Z3+5Z4+4Z5:E=0:F=Z6-2Prog GCLbl B子程序文件名:IN0A=Z3:A:AZKPQXDS:Z3=AA=Z4:A:AYKPQXDS:Z3=AA=Z5:A:AZKSQXDS:Z5=AA=Z6:A:AYKSQXDS:Z6=AA=Z7:A:AZKJDDK:Z7=AA=Z8:A:AYKJDDK:Z8=AA=Z12:A:AZKDLS:Z12=AA=Z13:A:AYKDLS:Z13=APLASE RUN:Dfem20+5Z3+5Z4+4Z5+4Z6+2Z12+2Z13 子程序文件名:IN1Lbl AC=5E+DE:Pause 0A=ZC+1:A:AJDN:ZC+1=AE:Pause 0A=ZC+2:A:AJDE:ZC+2=AE:Pause 0A=ZC+3:A:AJDR:ZC+3=AE:Pause 0A=ZC+4:A:AJDA1:ZC+4=AE:Pause 0A=ZC+5:A:AJDA2:ZC+5=AE=E+1Dsz FGoto A子程序文件名:IN2Lbl AClsC=4E+DE:Pause 0A=ZC+1:A:AJDDK:ZC+1=AE:Pause 0A=ZC+2:A:AJDH:ZC+2=AE:Pause 0A=ZC+3:A:AJDR:ZC+3=AE:Pause 0A=ZC+4:A:AJDI:ZC+4=AE=E+1Dsz FGoto A子程序文件名:IN3Lbl AClsC=2E+DE:Pause 0A=ZC+1:A:ADLDK1:ZC+1=AE:Pause 0A=ZC+2:A:ADLDK2:ZC+2=AE=E+1Dsz FGoto A子程序文件名:ZB坐标Lbl AE=E+1C=5E+DPol(zC+2-ZC-3,ZC+1-ZC-4)JJ=J+2H=JPol(zC+7-ZC+2,ZC+6-ZC+1)JJ=J+2K=J-HR=zC+3:L=zC+4:M=zC+5N=L/2RO=M/2RP=(Abs K-N-O)RX=L-L3/40R2Y=L2/6R-L4/336R3Q=Y-R(1-CosN)S=X-RSinNT=(R+Q)Tan Abs(K/2)+SX=M-M3/40R2Y=M2/6R-M4/336R3Q=Y-R(1-CosO)S=X-RSinOU=(R+Q)Tan Abs(K/2)+SGV-T=Goto 1Rec(G-V,H)X=J+zC+1Y=I+zC+2W=HGoto BLbl 1GV-T+L+P+M=Goto4GV-T+L= Goto 2Prog ZHHYGoto BLbl 2GV-T+L+P= Goto 3Prog HYYHGoto BLbl 3GV-T+L+P+M= Goto 4Prog YHHZGoto BLbl 4Pol(zC+7-ZC+2,ZC+6-ZC+1)V=V-T+L+P+M-U+IDsz FGoto AGV= Goto EQ=G-VW=JRec(Q,J)X=ZC+6+JY=ZC+7+IGoto BLbl EX=-1Lbl B子程序文件名:ZHHYRec(T,H)N=ZC+