卡西欧4800P曲线任意里程中边桩坐标正反算

卡西欧4800P曲线任意里程中边桩坐标正反算曲线任意里程中边桩坐标正反算(CASIO fx-4800P计算器)程序曲线任意里程中边桩坐标正反算(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=1 - P :D=(P-R) - (2HPR) : E=180 -n : N=1=Goto 1:工 Goto 2 Lbl 1 ：SZ ：SZ：W=Abs(S-O)：Prog SUB1 ：XXS=XJYY S=Y 丄FFS二F-90丄Goto 1 JLbl 2 : XY : XY: I=X : J=Y : Prog SUB2 : SS=O+W 丄ZZ=Z 丄Goto 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: ProgSUB1 : L=T+QEW(C+WD) : Z=(J-Y)cosL-(I-X)sinL: AbsZGoto1 :工 W=W+Z : Goto 0 XY2. XY = SZN ? 选择计算方式，输入 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 xx 正算时，计算得出的所求点的X坐标YS=x xx 正算时，计算得出的所求点的丫坐标FS=x xx 正算时，所求点对应的中线点的切线方位角S=x xx 反算时，计算得出的所求点的里程Z=x xx 反算时，计算得出的所求点的边距四、算例 某匝道的由五段线元（直线+ 完整缓和曲线+ 圆曲线+ 非完整缓和曲线+直线）组成，各段线元的要素（起点里程SO、起点坐标X0 Y0、起点切线方位角 F0 、线元长度 LS 、起点曲率半径 R0 、止点曲率半径 RN、 线 元左右偏标志 Q ）如下：S0 X0 Y0F0LS R0 RN Q500.000 19942.837 28343.561 12516 .00 269.2561E45 1E45 0769.25619787.34028563.3781251631.0037.4921E45 221.75 -1806.74819766.56628594.5741202554.07112.779221.75 221.75 -1919.52719736.07228701.893911730.6380.285221.75 9579.228-1999.81219744.03828781.659804050.00100.0001E45 1E45 01 、正算注意：略去计算方式及线元要素输入，请自行根据所求点所在的线元输入线元要素）S=700 Z=-5得 XS=19831.41785YS=28509.72590S=700 Z=0得 XS=19827.33592YS=28506.83837S=700 Z=5 计得 XS=19823.25398YS=28503.95084S=780 Z=-5算算算算得 XS=19785.25749YS=28575.02270S=780 Z=0计算得 XS=19781.15561YS=28572.16358S=780 Z=5计算得 XS=19777.05373YS=28569.30446S=870 Z=-5计算得 XS=19747.53609YS=28654.13091S=870 Z=0计算得 XS=19742.68648YS=28652.91379S=870 Z=5计算得 XS=19737.83688YS=28651.69668S=940 Z=-5.123计算得XS=19741.59118 YS=28722.05802S=940 Z=0计算得 XS=19736.47687YS=28722.35642S=940 Z=3.009计算得 XS=19733.47298YS=28722.531682 、 反算X=19831.418 Y=28509.726计算得 S=699.9999974 Z= -5 .00018164X=19827.336Y=28506.838计算得 S=699.9996493 Z= 0.000145136X=19823.25398Y=28503.95084计算得 S=699.9999985 Z= 5.000003137X=19785.25749Y=28575.02270计算得 S=780.0000035 Z= -5 .000001663X=19781.15561Y=28572.16358计算得 S=780.0000025 Z=- 0.000002979X=19777.05373Y=28569.30446计算得 S=780.0000016 Z= 4.99999578X=19747.536Y=28654.131计算得 S=870.0001137 Z= -4.99941049X=19742.686Y=28652.914计算得 S=870.0003175 Z=- 0.00041814X=19737.837Y=28651.697计算得 S=870.0002748 Z= 4.999808656X=19741.5912Y=28722.0580计算得 S=939.9999786 Z= -5.123024937X=19736.4769Y=28722.3564计算得 S=939.9999862 Z=- 0.000027710X=19733.4730Y=28722.5317计算得 S=940.0000238 Z= 3.00898694