对卡西欧4800计算器编程

对卡西欧 4800 计算器在编程中的浅见作者: luzhmu 发布日期: 2005-10-08 查看数: 793 出自: 在计算能力飞速发展的今天，各种计算工具应运而生代替了以往我们大量、繁琐的分析计算工作，简化了计算步骤节省了计算时间。在种类繁多的计算工具中适合测量工作的计算器也大有所在。比如有一款 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=NW 2G:I=T-W+W540V2+W 36V: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(0.5A):M=E(N-1):W=NL2(24N-24):Q=M2-4W:Q0“R=”:0.5(M+Q)程序说明:本程序适用于外距控制单交点平曲线PJ? 输入路线转角E? 输入外距LS? 输入缓和曲线长2006-03-28 17:52:15大中小程序名:LZXJD(两直线交点)A“X1“:B“Y1“:C“X2“:D“Y2“:E“X3“:F“Y3“:G“X4“:H“Y4“:N=C-A:O=D-B:P=G-E:Q=H-F:K=QN-OP:U=NB-OA:V=PF-QE:“X=“:(PU-NV)/K“Y=“:(QU-OV)/K程序说明:X1? 输入第一点 X坐标Y1? 输入第一点 Y坐标X2? 输入第二点 X坐标Y2? 输入第二点 Y坐标X3? 输入第三点 X坐标Y3? 输入第三点 Y坐标X4? 输入第四点 X坐标Y4? 输入第四点 Y坐casio4850 全线坐标计算程序 一、LXZB（坐标计算主程序）LbI Z： PA=： V=：C=：D=：P直线的终点桩号Prog“ZX.SUB”：Goto Z（以上各字母含义详直线坐标计算子程序说明） A=： V=：C=：D=：R=： I=（1 或-1）P圆曲线终点桩号Prog“YQX.SUB”：Goto Z（以上各字母含义详圆曲线坐标计算子程序说明）A=： V=：C=：D=：B=： W=：E=：F=：R=： L=：I=（1 或-1）PHY 点桩号Prog“HQ1.SUB” ：Goto ZPYH 点桩号Prog“HYQ.SUB”：Goto ZPHZ 点桩号Prog“HQ2.SUB”：Goto Z（以上各字母含义详（HQ1.SUB、HYQ.SUB、HQ2.SUB）坐标计算子程序说明）该主程序是简单的示例说明，具体的应根据图纸的线型组合，随机应变的做顺序调整，并将其输入 Casio4850 计算器。二、ZX.SUB（直线坐标计算子程序）G=P-A：T=V“X=”：X=C+GcosT “Y=”：Y=D+GsinT Q：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+JcosV Q: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 （599040 R6L6）J=(G3(6RL)-G7(336R3L3)+G11(42240R5L5)IT=V+90G2(RL)I“X=”:X=C+HcosV-JsinV “Y=”:Y=D+HsinV+JcosV Q: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+JcosV Q: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 （599040 R6L6）J=-(G3(6RL)-G7(336R3L3)+G11(42240R5L5)IT=W-90G2(RL)I“X=”:X=E-HcosW+JsinW “Y=”:Y=F-HsinW-JcosW Q: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-O cos(T+90) “YZ=”：N=Y-O sin(T+90) “XY=”：K=X+S cos(T+90) “YY=”：U=Y+S sin(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 “SUBSJK“2.(SJK1 数据库程序)“1.SZ = XY“：“2.XY = SZ“：N,S：N：S“DKI“S本线元终点里程U= 本线元起点 X 坐标：V=本线元起点 Y 坐标：G=本线元起点正切线方位角：P=本线元起点曲率半径：R=本线元终点曲率半径：O=本线元起点里程：H= 本线元终点里程： Q=本线元左直右偏向(0 或 1、-1):Goto0S下一线元终点里程= O=。： Goto0LB1 03.(SUBSJK 运算主程序)UVGPROHQ：U“XA“：V“YA“ ：G“FA“：P“RA“：R“RB“：O“DKA“：H“DKB“：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=“：X“YS=“：Y“QX-FWJ=“： F=F-90DMSI=0：J=0 ：Pol（X-C“XC“，Y-E“YC“）：JGoto1：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：U“X0“：V“Y0“：O“S0“：G“F0“：H“LS“：P“R0“：R“RN“：Q：C=1P：D=(P-R)(2HPR)：E=180：N=1=Goto 1：Goto 2Lbl 1：SZ：SZ：W=Abs(S-O) ：Prog “SUB1“：X“XS“=XY“YS“=YF“FS“=F-90Goto 1Lbl 2：XY：XY ：I=X：J=Y：Prog “SUB2“：S“S“=O+WZ“Z“=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+QEW(C+WD)+90：X=X+ZcosF ： Y=Y+ZsinF3. 反算子程序 (SUB2)T=G-90：W=Abs(Y-V)cosT-(X-U)sinT)：Z=0 ：Lbl 0：Prog “SUB1“：L=T+QEW(C+WD)：Z=(J-Y)cosL-(I-X)sinL：AbsZGoto1：W=W+Z ：Goto 0Lbl 1：Z=0 ：Prog “SUB1“：Z=(J-Y)sinF三、使用说明1、规定(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。(2) 当所求点位于中线时，Z=0；当位于中线左铡时，Z 取负值；当位于中线中线右侧时，Z 取正值。(3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以 10 的 45 次代替。(4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。(5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以 10 的 45 次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以 10 的 45 次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。(6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。2、输入与显示说明输入部分：1. SZ = XY2. XY = SZ N ? 选择计算方式，输入 1 表示进行由里程、边距计算坐标 ；输入 2 表示由坐标反算里程和边距。X0 ？线元起点的 X 坐标 Y0 ？线元起点的 Y 坐标 S0 ？线元起点里程F0 ？线元起点切线方位角LS ？线元长度R0 ？线元起点曲率半径RN ？线元止点曲率半径Q ？ 线 元左右偏标志(左偏 Q=-1，右偏 Q=1，直线段 Q=0）S ？ 正算时所求点的里程Z ？正算时所求点距中线的边距(左侧取负，值右侧取正值，在中线上取零 )X ？反算时所求点的 X 坐标Y ？反算时所求点的 Y 坐标显示部分：XS= 正算时，计算得出的所求点的 X 坐标YS= 正算时，计算得出的所求点的 Y 坐标FS= 正算时，所求点对应的中线点的切线方位角S= 反算时，计算得出的所求点的里程Z= 反算时，计算得出的所求点的边距四、算例某匝道的由五段线元（直线+完整缓和曲线+ 圆曲线 +非完整缓和曲线+ 直线）组成，各段线元的要素（起点里程 S0、起点坐标 X0 Y0、起点切线方位角 F0、线元长度LS、起点曲率半径 R0、止点曲率半径 RN、线 元左右偏标志 Q）如下：S0 X0 Y0 F0 LS R0 RN Q500.000 19942.837 28343.561 125 16 31.00 269.256 1E45 1E45 0769.256 19787.340 28563.378 125 16 31.00 37.492 1E45 221.75 -1806.748 19766.566 28594.574 120 25 54.07 112.779 221.75 221.75 -1919.527 19736.072 28701.893 91 17 30.63 80.285 221.75 9579.228 -1999.812 19744.038 28781.659 80 40 50.00 100.000 1E45 1E45 0 1、正算（注意：略去计算方式及线元要素输入，请自行根据所求点所在的线元输入线元要素）S=700 Z=-5 计算得 XS=19831.41785 YS=28509.72590S=700 Z=0 计算得 XS=19827.33592 YS=28506.83837S=700 Z= 5 计算得 XS=19823.25398 YS=28503.95084S=780 Z=-5 计算得 XS=19785.25749 YS=28575.02270S=780 Z=0 计算得 XS=19781.15561 YS=28572.16358S=780 Z= 5 计算得 XS=19777.05373 YS=28569.30446S=870 Z=-5 计算得 XS=19747.53609 YS=28654.13091S=870 Z=0 计算得 XS=19742.68648 YS=28652.91379S=870 Z= 5 计算得 XS=19737.83688 YS=28651.69668S=940 Z=-5.123 计算得 XS=19741. 59118 YS=28722.05802S=940 Z=0 计算得 XS=19736.47687 YS=28722.35642S=940 Z= 3.009 计算得 XS=19733.47298 YS=28722.531682、 反算X=19831.418 Y=28509.726 计算得 S=699.9999974 Z= -5 .00018164X=19827.336 Y=28506.838 计算得 S=699.9996493 Z= 0.000145136X=19823.25398 Y=28503.95084 计算得 S=699.9999985 Z= 5.000003137X=19785.25749 Y=28575.02270 计算得 S=780.0000035 Z= -5 .000001663X=19781.15561 Y=28572.16358 计算得 S=780.0000025 Z=- 0.000002979X=19777.05373 Y=28569.30446 计算得 S=780.0000016 Z= 4.99999578X=19747.536 Y=28654.131 计算得 S=870.0001137 Z= -4.99941049X=19742.686 Y=28652.914 计算得 S=870.0003175 Z=- 0.00041814X=19737.837 Y=28651.697 计算得 S=870.0002748 Z= 4.999808656X=19741.5912 Y=28722.0580 计算得 S=939.9999786 Z= -5.123024937X=19736.4769 Y=28722.3564 计算得 S=939.9999862 Z=- 0.000027710X=19733.4730 Y=28722.5317 计算得 S=940.0000238 Z= 3.00898694常用 CASIO4800 计算程序集FX4800P 常用测量程序集1、由一个已知点坐标和一个已知方位角计算支导线坐标（程序名：ZBJS）程序清单：X：Y：T ：Lb1 0：Fixm：AS：T=T+A-180：T0 =T=T+360 T360 =T=T-360T X=X+Scos TY=Y+Ssin TGoto 0程序运行说明：运行“ZBJS”，每运行一步按 EXE 键一次，显示 X？时输入已知点 X 坐标，显示 Y？时输入已知点 Y 坐标，显示 T？时输入已知点由后视到前视的已知坐 标方位角，运行一次本程序只需输入一次已知数据。显示 A？时输入第一个左转角值，显示S？时输入测站到前视的水平距离。接下来计算器显示的 T 是测站到前视的坐标方位角，按 + 可将显示的十进制角度换算为度分秒。如此继续可计算出所有支导线未知坐标。2、由两个已知点坐标计算支导线坐标（程序名：ZBJS2）程序清单：M“X0”：N“Y0”：D“XH”：E“YH” Pol(M-D，N-E)：Fixm：T=J：X=M ：Y=N：T0=T=T+360T“TH-J=” I“S=”Lb1 0:AS: T=T+A-180：T0 =T=T+360 T360 =T =T-360TX=X+Scos TY=Y+Ssin TGoto 0程序运行说明：运行后，显示 X0，Y0 时输入已知置镜点坐标，显示 XH，YH 时输入已知后视点坐标，其它同 ZBJS 的输入与显示。3、由一个已知点坐标和一个已知方位角进行放射型导线坐标计算（程序名 FSXZBJS）程序清单：C“X0”：D“Y0”：E“T0” ：Fixm ：Lb1 1:SAT=E+A-180：T 0=T=T+360 T360 =T=T-360TX=C+Scos TY=D+Ssin TGoto 1程序运行说明：变量输入参照 ZBJS 程序。4、由两个已知点坐标进行放射型导线坐标计算（程序名 FSXZBJS2）M“X0”：N“Y0”：D“XH”：E“YH