《线路中边桩任意里程坐标正反算改》由会员分享，可在线阅读，更多相关《线路中边桩任意里程坐标正反算改（5页珍藏版）》请在金锄头文库上搜索。
1、线路中边桩任意里程坐标正反算改 含数据库、测站放样计算（第五次修改-正式版）1.加编数据库及计算总调度程序,计算中不必逐项输入“线元要素“,提高运算速度，避免现场忙中出错2.可加入多条线路的数据库,内业输入数据,外业一目了然3.计算直观,人性化4.正算直接输入里程和边距,反算输入近似里程便可5.增加了“计算点与测站点”的距离和方位角计算语句，方便直接放样6.愿收获与大家共享,同时也希望大家提出心得和宝贵建议7.核心计算程序摘自“yshf”一.改动后的程序清单，增设数据库程序,可输入 n 条线路的数据库1.(QXZDJS 计算总调度程序)M=1=Prog “SJK1“:Prog “SUBSJK“
2、.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=
3、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”运行，程序
4、提示输入里程“S”？时，正算直接输入待求点里程，反算输入所求点“近似“里程3.如果在正算中途想进行反算，重新给”N“？赋值“2“便可4.程序中“XC、YC”为测站坐标，“I”、“FWJ”为放样距离和方位角5.程序中程序中“QX-FWJ=”为计算点在中线上的正切线方位角，此方位角在新线路输入数据库时非常重要，我为计算点在中线上的正切线方位角，此方位角在新线路输入数据库时非常重要，我们在每输入完一段曲线元要素后，立即退出返回到们在每输入完一段曲线元要素后，立即退出返回到“计算总调度程序计算总调度程序(SUBSJK)”中，然后输入该段曲线元中，然后输入该段曲线元的终点里程，其计算出来的的终点里程，其
5、计算出来的 X，Y 中桩坐标应该与设计图纸中桩坐标应该与设计图纸(当然是经过复查无误的图纸当然是经过复查无误的图纸)上的该点中桩坐上的该点中桩坐标一致，否则一定是哪个地方数据输入或者程序有误，这时候我们一定要把该计算点在中线上的正切线方标一致，否则一定是哪个地方数据输入或者程序有误，这时候我们一定要把该计算点在中线上的正切线方位角位角“QX-FWJ=”记录下来为下一段曲线元要素数据输入服务，因为它就是下一段曲线元的起点正切线方位记录下来为下一段曲线元要素数据输入服务，因为它就是下一段曲线元的起点正切线方位角，嘿嘿！角，嘿嘿！6.其它输入与“yshf“原程序一致-附件:(原 yshf 程序)曲线
6、任意里程中边桩坐标正反算(CASIO fx-4800P 计算器)程序一、程序功能本程序由一个主程序(TYQXJS)和两个子程正算子程序(SUB1)、反算子程序(SUB2)序构成，可以根据曲线段直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。另外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲线任意里程中边桩坐标进行正反算。本程序也可以在 CASIO fx-4500P 计算器及 CASIO fx-4850P 计算器上运行。特别申明：
7、(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“
8、=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. 反算子
9、程序(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、率半径为无穷大，以 10 的 45 次代替。(4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。(5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以 10 的 45 次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以 10 的 45 次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。(6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。2、输入与显示说明输
11、入部分：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 坐
12、标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.3
13、40 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=198
14、31.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
