线路中边桩任意里程坐标正反算改编程序CASIO4850.doc

线路中边桩任意里程坐标正反算改编程序CASIO4850作者：白黑文章来源：本站原创点击数： 742更新时间：2006-11-8含数据库、测站放样计算（第五次修改-正式版）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+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 = XY 2. 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.72590 S=700 Z=0 计算得 XS=19827.33592 YS=28506.83837 S=700 Z= 5 计算得 XS=19823.25398 YS=28503.95084 S=780 Z=-5 计算得 XS=19785.25749 YS=28575.02270 S=780 Z=0 计算得 XS=19781.15561 YS=28572.16358 S=780 Z= 5 计算得 XS=19777.05373 YS=28569.30446 S=870 Z=-5 计算得 XS=19747.53609 YS=28654.13091 S=870 Z=0 计算得 XS=19742.68648 YS=28652.91379 S=870 Z= 5 计算得 XS=19737.83688 YS=28651.69668 S=940 Z=-5.123 计算得 XS=19741. 59118 YS=28722.05802 S=940 Z=0 计算得 XS=19736.47687 YS=28722.35642 S=940 Z= 3.009 计算得 XS=19733.47298 YS=28722.53168 2、 反算 X=19831.418 Y=28509.726 计算得 S=699.9999974 Z= -5 .00018164 X=19827.336 Y=28506.838 计算得 S=699.9996493 Z= 0.000145136 X=19823.25398 Y=28503.95084 计算得 S=699.9999985 Z= 5.000003137 X=19785.25749 Y=28575.02270 计算得 S=780.0000035 Z= -5 .000001663 X=19781.15561 Y=28572.16358 计算得 S=780.0000025 Z=- 0.000002979 X=19777.05373 Y=28569.30446 计算得 S=780.0000016 Z= 4.99999578 X=19747.536 Y=28654.131 计算得 S=870.0001137 Z= -4.99941049 X=19742.686 Y=28652.914 计算得 S=870.0003175 Z=- 0.00041814 X=19737.837 Y=28651.697 计算得 S=870.0002748 Z= 4.999808656 X=19741.5912 Y=28722.0580 计算得 S=939.9999786 Z= -5.123024937 X=19736.4769 Y=28722.3564 计算得 S=939.9999862 Z=- 0.000027710 X=19733.4730 Y=28722.5317 计算得 S=940.0000238 Z= 3.00898694casio4850全线坐标计算程序该程序是本人自己编的适用于casio4850全线坐标计算，发布上来给大家参考，请大家试用指教，如有问题请修改，首先是方位角显示不出来，请高手修改。一、LXZB（坐标计算主程序）LbIZ：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)