曲线任意里程中边桩坐标正反算修改程序5800.docx

1、曲线任意里程中边桩坐标正反算修改程序5800曲线任意里程中边桩坐标正反算修改程序本程序由一个主程序(ZB)和四个子程正算子程序(SUB1)、反算子程序(SUB2)、放样子程序（SUB3）、数据源子程序(SUB0)序构成。此程序可以根据曲线段直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。也可以在CASIO fx-4500P计算器及 CASIO fx-4800P计算器上运行。1.主程序(ZB)1.SZ = XY：2.XY = SZ：N：N=2=Goto2：

2、Goto 1Lbl 1：SZN：S：Prog SUB0：Z：N=90：W=Abs(S-O)：Prog SUB1：XS：XYS： YProg SUB3Goto 1Lbl 2：XYS：XYS：I=X：J=Y：Prog SUB2：S：S=O+WZ： ZGoto 2注：如将N=90改为N，可计算斜交。2. 正算子程序(SUB1)A=0.1739274226：B=0.3260725774：K=0.0694318442：L=0.3300094782：F=1-L：M=1-K：X=U+W(Acos(G+Q(180)KW(1P)+KW(P-R) (2HPR)+Bcos(G+Q(180)LW(1P)+LW(P-R

3、) (2HPR)+Bcos(G+Q(180)FW(1P)+FW(P-R) (2HPR)+Acos(G+Q(180)MW(1P)+MW(P-R) (2HPR)：Y=V+W(Asin(G+Q(180)KW(1P)+KW(P-R) (2HPR)+Bsin(G+Q(180)LW(1P)+LW(P-R) (2HPR)+Bsin(G+Q(180)FW(1P)+FW(P-R) (2HPR)+Asin(G+Q(180)MW(1P)+MW(P-R) (2HPR)：F=G+Q(180)W(1P)+W(P-R) (2HPR)：X=X+Zcos（F+N）：Y=Y+Zsin（F+N）注：如应用Defm 命令用Zi将部分

4、字母代替出来，用以扩充变量，以实现二次运行同一程序时输入值不变。3. 反算子程序(SUB2)Lbl 0：Prog SUB0：T =G-90：W=Abs(Y-V)cosT-(X-U)sinT)：WH=S=S+H：Goto 0Z=0：Lbl 1：N=90：Prog SUB1：L=T+Q(180)W(1P)+W(P-R) (2HPR)：Z=(J-Y)cosL-(I-X)sinL：AbsZGoto 2：W=W+Z：Goto 1Lbl 2：Z=0：Prog SUB1：Z=(J-Y)sin（F+90）注：SUB1、SUB2中(1P) 、(P-R) (2HPR)、 (180)为原程序中的C、D、E。4. 放

5、样子程序(SUB3)CDET：CXJ：D：YJ：EXH：T：YH：X：Y：I=0：J=0：Pol（E-C，T-D）：J J=J+360HSFW：J-DMSN=J：HSJL：IPol（X-C，Y-D）：J J=J+360QSFW：J-DMSN=J-N：N N=N+360QSZJ：N-DMSQSJL：I5. 数据源子程序(SUB0)：对 OS0：UX0：VY0： GF0：HLS：PR0：RRN：Q 进行赋值。6数据库（SUB4）S500= S O =500:U =19942.837:V = 28343.561:G=125。16。 31.00。：H=269.256：P 1e45: R=1e45:Q=

6、0S769.256= S O =769.256:U =19787.340:V = 28563.378:G=125。16。 31.00。：H=37.492：P 1e45: R=221.75:Q=-1S806.748= S O =806.748:U =19766.566:V = 28594.574:G=120。25。 54.07。：H=112.779：P 221.75: R=221.75:Q=-1S919.527= S O =919.527:U =19736.072:V = 28701.893:G=91。17。 30.63。：H=80.285：P 221.75: R=9579.228:Q=-1S9

7、99.812= S O =999.812:U =19744.038:V = 28781.659:G=80。40。 50.00。：H=100.000：P 1e45: R=1e45:Q=0三、使用说明1、规定 (1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时， Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。 (2) 当所求点位于中线时，Z=0；当位于中线左侧时，Z取负值；当位于中线中线右侧时，Z取正值。 (3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。 (4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半

8、径。 (5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径 (6) 当线元为非完整缓和曲线时，起点（止点）与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。2、输入与显示说明 输入部分： 1. SZ = XY 2. XY = SZ N ? 选择计算方式，只有输入2时进行坐标反算，其余均为正算。 S0 ？线元起点里程X0 ？线元起点的X坐标 Y0 ？线元起点的Y坐标 LS ？线元长度F0 ？线元起点切线方位角 R0 ？线元起点曲率半径 RN ？线元止点曲率半径 Q ？ 线元左右偏标志(左偏

9、Q=-1，右偏Q=1，直线段Q=0） S ？ 正算时所求点的里程 Z ？正算时所求点距中线的边距(左侧取负，值右侧取正值，在中线上取零) XJ？正算时置镜点X坐标YJ？正算时置镜点Y坐标XH？正算时后视点X坐标YH？正算时后视点Y坐标X ？反算时所求点的X坐标 Y ？反算时所求点的Y坐标 S ？反算时给定一初始里程（必须起点里程小于计算点里程，无法判断时输起始里程，越接近计算越快）。显示部分： XS= 正算时，计算得出的所求点的X坐标 YS= 正算时，计算得出的所求点的Y坐标 HSFW 正算时，置镜点到后视点的后视方位角HSJL 正算时，置镜点到后视点的后视距离QSFW 正算时，置镜点到前视点

10、的前视方位角QSZJ 正算时，视线自后视到前视的顺时针转角QSJL 正算时，置镜点到前视点的前视距离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

11、 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 计算得 X

12、S=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=8

13、70 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.0001

14、45136 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.

15、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.00898694

16、3、放样例：以起点为置镜点XJ=19942.837， YJ=28343.561 后视点XH=19948.837，YH=28343.561 S=700 Z=-5 计算得 XS=19831.41785 YS=28509.72590 HSFW=0。00。00。 HSJL=6 QSFW=123。50。35.45。 QSZJ=123。50。35.45。 QSJL=200.0624902 S=700 Z=0 计算得 XS=19827.33592 YS=28506.83837HSFW=0。00。00。 HSJL=6 QSFW=125。16。31。 QSZJ=125。16。31。 QSJL=200 S=700

17、 Z= 5 计算得 XS=19823.25398 YS=28503.95084HSFW=0。00。00。 HSJL=6 QSFW=126。42。26.55。 QSZJ=126。42。26.55。 QSJL=200.0624902注意事项：1、程序使用前状态：（1） 角度必须60进制，屏幕显示“D”（2） 屏幕设置Norm 1，设置方式Function=6DSP/CLR=3Norm=1设定2、当使用数据库运算时，计算的里程必须在最小里程和最大里程范围之内，否则无法显示。例如在运行程序6中的数据库时所输入的里程S XY：2.XY = SZ：N：UX0：VY0：OS0：GF0：HLS：PR0：RRN

18、：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)+Aco

19、s(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、规定

20、 (1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。 (2) 当所求点位于中线时，Z=0；当位于中线左铡时，Z取负值；当位于中线中线右侧时，Z取正值。 (3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。 (4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。 (5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时

21、，曲率半径等于圆曲线的半径。 (6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。2、输入与显示说明 输入部分： 1. SZ = XY 2. XY = SZ N ? 选择计算方式，输入1表示进行由里程、边距计算坐标 ；输入2表示由坐标反算 里程和边距。 X0 ？线元起点的X坐标 Y0 ？线元起点的Y坐标 S0 ？线元起点里程 F0 ？线元起点切线方位角 LS ？线元长度 R0 ？线元起点曲率半径 RN ？线元止点曲率半径 Q ？ 线 元左右偏

22、标志(左偏Q=-1，右偏Q=1，直线段Q=0） S ？ 正算时所求点的里程 Z ？正算时所求点距中线的边距(左侧取负，值右侧取正值，在中线上取零) X ？反算时所求点的X坐标 Y ？反算时所求点的Y坐标 显示部分： XS= 正算时，计算得出的所求点的X坐标 YS= 正算时，计算得出的所求点的Y坐标 FS= 正算时，所求点对应的中线点的切线方位角 S= 反算时，计算得出的所求点的里程 Z= 反算时，计算得出的所求点的边距验证过的数据, 缓和曲线左偏 S0=87418.227X0=4581068.056 Y0=484543.0267 F0=3472550.78 LS=340 R0=1045RN=3

23、000Q=-1 (已验证过:正算正确反算正确)第一缓和曲线段桩号:正算K87+460左(-)13ms?87460z?-13 XS=4581105.985 YS=484521.2385K87+460中线上0m s?87460 z?0 XS=4581108.825 YS=484533.9245K87+460右(+)13m s?87460 z?13XS=4581111.665 YS=484546.6105反算X=4581105.985Y=484521.2385 S=87459.99964Z=-13.00004001X=4581108.825Y=484533.9245 S=87459.99965 Z=

24、-0.00005237499X=4581111.665Y=484546.6105 S=87459.99984Z=12.99997507K87+600左(-)13m XS=4581242.22YS=484489.8717K87+600中线上0m XS=4581245.254YS=484502.5126K87+600右(+)13m XS=4581248.289 YS=484515.1536反算X=4581242.22 Y=484489.8717 S=87599.99981 Z=-13.00001607X=4581245.254Y=484502.5126 S=87599.99982 Z=-0.000

25、06918891X=4581248.289Y=484515.1536 S=87600.00043 Z=13.00012293K87+740左(-)13mXS=4581377.383YS=484455.1799K87+740中线上0mXS=4581380.852 YS=484467.7085K87+740右(+)13mXS=4581384.321YS=484480.2371反算X=4581377.383 Y=484455.1799S=87740.00013 Z=-12.99998575X=4581380.852Y=484467.7085S=87740.00013Z=-0.00000001556X

26、=4581384.321 Y=484480.2371 S=87739.9999Z=12.99999725K87+758.227左(-)13m XS=4581394.861YS=484450.2853K87+758.227中线上0m XS=4581398.404YS=484462.7931K87+758.227右(+)13m XS=4581401.947YS=484475.3009反算X=4581394.861Y=484450.2853 S=87758.22744 Z=-12.99983181X=4581398.404Y=484462.7931 S=87758.22741 Z=0.00011225994X=4581401.947Y=484475.3009 S=87758.22721 Z=13.00000741姓名：wendi168506-1-20 21:36评：+5 分 S0=87758.227X0=4581398.409 Y0=484462.791 F0=344113 LS=2603.31 R0=3000 RN=3000