4800万能程序[精校版本]

1、07-1-11楼主：MACHAO作者：马超 本文已被浏览30727次CASIO万能坐标程序MC2008.01.29程序特点：真正的全线贯通坐标计算，在曲线元要素输入时仅需要输入第一段全部曲线元要素，后面曲线元要素除起点半径、终点半径、曲线长、转向需输入外其他要素均从前一曲线按辛普森8等分计算得出；辛普森公式任意等分，满足所有精度要求；全线曲线元数据一次性程序化输入，扩充变量数据库，无需修改程序内容；傻瓜操作，适用初级用户。一、程序：MCW“1.JS 2.SZ”:W=1=Z2=0:Goto 1 W=2= O“KOULING”:O123456=O=0: “OUT” Goto 4O=0: Z1=0:

2、Goto 0Lbi 0”N0.”：Z1+1 Z1=0=ABCREFGU:A“X0”:B“Y0”:C“F0”:R“R0”:E“RN”:F“D0”:G“LS”:U“G” : ZZ18+3=A:ZZ18+4=B:ZZ18+5=C:ZZ18+6= R-1:ZZ18+7= E-1:ZZ18+8=F: ZZ18+9=F+G: ZZ18+10=U: “NEXT” Isz Z1: Goto 0Z1=1=D=Z9:Z=0:Z2=0:GOTO 2D=Z(Z1-1)8+9:Z=0:Z2=Z1-1:GOTO 2Lbi A: ZZ18+3=X:ZZ18+4=Y:ZZ18+5=J: ZZ18+8=D: REGU:R“Ro

3、”：E“RN”: G“LS”:U“G”: ZZ18+6=R-1 :ZZ18+7=E-1: ZZ18+9=D+G: ZZ18+10=U: “NEXT” Isz Z1: Goto 0Lbi 1DZ:D：Z：Z2=0:Goto 2Lbi 2Z2Z1=GoTo 4DZZ28+9=A=ZZ28+3:B=ZZ28+4: C =ZZ28+5:R=ZZ28+6: E=ZZ28+7: F=ZZ28+8: G=ZZ28+9: U=ZZ28+10: Goto3Isz Z2:Goto 2Lbi 3W=2 =N=8：N=5PU（E-R）Abs（G-F）：QAbs（D-F）N：S=90Q： JC+(NPQ+2UR)NS：

4、L=1XA+Q6(Cos C+Cos J +4（Cos (C+(L+0.5)PQ+2UR)(L+.5)S),L,0,(N-1)）+2（Cos (C+(LPQ+2UR)LS）,L,1,(N-1)）+ZCos（J+90) YB+Q6(Sin C+Sin J +4（Sin (C+(L+0.5)PQ+2UR)(L+.5)S),L,0,(N-1)）+2（Sin (C+(LPQ+2UR)LS）,L,1,(N-1)）+ZSin（J+90):W=2=GOTO AZ=0=“X”:X:Pause 0: “Y” :Y Goto 1Z“XL”:X:Pause 0: “YL”:Y Goto 1 fx4850 Z0=“X

5、R”:X:Pause 0: “YR”:Y Goto 1 Z=0= X “X” Y “Y” Goto 1Z X “XL” Y “YL” Goto 1 fx4800 Z0= X “XR” Y “YR” Goto 1 Lbi 4二、说明a、编制说明本程序是运用复化辛普生公式根据曲线段直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距，对该曲线段范围内任意里程中边桩坐标进行计算，以及对卡西欧扩充变量的灵活应用，实现了真正意义上的的全线贯通及曲线要素输入程序化（在不修改程序内容的情况下可通过运行程序输入任意多段曲线元要

6、素）。通过对N=？ 进行修改，可对辛普森公式进行任意等分进行运算。（注：N为不小于2的整数，N越大精度越高，计算速度越慢；N越小精度越低，运算速度越快，一般曲线取N=4就能满足精度要求，在能满足精度的情况下尽量N取小值，已获得最佳运算速度，不要盲目的追求精度）。b、程序操作说明程序分为两部分：1.JS为计算，2.SZ为设置。1、首先用Shift+Defm键对计算器内存变量进行扩充，扩充变量数为8X+2（X为曲线元段数，变量数视内存情况尽量大些），运行程序，选2进行曲线要素设置，为防止误操作在正确输入口令123456时方可进行设置，否则显示OUT跳到程序尾（LBI4）。在执行程序中将第一段曲线元

7、要素按规律输入到扩充内存变量中（A“X0”:B“Y0”:C“F0”:R“R0”:E“RN”:F“D0”:G“LS”:U“G”），以后曲线仅需输入R“R0”、E“RN” :G“LS”:U“G”，其他参数自动计算得出。在显示NEXT时为提示是否输入下一曲线，按EXE继续输入，如不需再输入则退出程序。2、选1.JS进行坐标运算，D输入桩号，Z输入左右距离（负为左，正为右，0为中）程序首先自动判断其在哪一线元内，并把其线元要素调出进行计算，即可对全线进行坐标计算。c、变量说明X0：Y0：F0曲线元起点X、Y坐标及起点正切线方位角R0：RN曲线元起点及终点半径D0：LS：Q 曲线元起点桩号、路线长度及线

8、路左右偏标志（左=-1,直线=0,右=1）D 曲线元中待求点桩号Z 计算边桩距中线平距，左边输入负值，右边输入正值，中桩输入0X：Y（中）XL:YL （左）XR:YR （右）待求点的X，Y坐标扩充变量：ZZ18+3: ZZ18+4: ZZ18+4: ZZ18+6: ZZ18+7: ZZ18+8: ZZ18+9: ZZ18+10: 分别为各线元X0:Y0:F0;R0:RN:D0:LS:G,(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。(2) 当所求点位于中线时，Z=0，坐标显示XY；当位于中线左侧时，Z取负值,

9、坐标显示XLYL,；当位于中线右侧时，Z取正值,坐标显示XRYR。(3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次方代替。(4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。(5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次方代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次方代替；与圆曲线相接时，曲率半径等于圆曲线的半径。(6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，

10、曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。/color(7)当两段缓和曲线相连时，相连处要素均输为10的45次方（可以看作缓直+直缓）。B本程序可在48004850上运行，只需将输出部分进行修改。（48504800）/B郑石高速路面NO.9标余官营互通E匝道序号XoYoFoRoRND0LSG1744383.112510406.45643。3400.0700070000170-1中点桩号中点X中点Y终点桩号终点X终点Y85744445.055510464.662170744507.700 510522.1122744507.700 510522.112 42。1030.

11、87000500170116.0711中点桩号中点X中点Y终点桩号终点X终点Y228.036744550.193510561.637286.071744590.136510603.7133744590.136510603.71349。1802.3500500286.071299.2821中点桩号中点X中点Y终点桩号终点X终点Y435.712744669.415510729.969585.353744707.945510873.9874744707.945510873.98783。3544.8500无穷大585.353601中点桩号中点X中点Y终点桩号终点X终点Y615.353744710.54

12、5510903.872645.353744712.247510933.8232008年1月27日更新，加入了坐标反算功能，坐标正算加入斜角计算功能。程序特点：真正的全线贯通坐标正反计算！在曲线元要素输入时仅需要输入第一段全部曲线元要素，后面曲线元要素除起点半径、终点半径、曲线长、转向需输入外其他要素均从前一曲线按辛普森8等分计算得出，解决了主线坐标计算无法获得第二段及其以后曲线元起点参数的问题；辛普森公式任意等分，满足所有精度要求；全线曲线元数据一次性程序化输入，参数存储采用扩充变量数据库，无需修改程序内容；多功能采用单程序编程，苊馄捣钡饔米映绦颍岣咴怂闶涠龋簧倒喜僮鳎乇鹗视贸跹摺?一、程序：

13、MCW“1.ZS 2.FS 3.SZ”:W=1=Z2=0:V=0：Goto 1 W=2= Goto 4W=3= O“KOULING”:O123456=O=0: “OUT” Goto DO=0: V=0:Z1=0:Goto 0Lbi 0”N0.”：Z1+1 Z1=0=ABCREFGU:A“X0”:B“Y0”:C“F0”:R“R0”:E“RN”:F“D0”:G“LS”:U“G” : ZZ18+3=A:ZZ18+4=B:ZZ18+5=C:ZZ18+6= R-1:ZZ18+7= E-1:ZZ18+8=F: ZZ18+9=F+G: ZZ18+10=U: “NEXT” Isz Z1: Goto 0Z1=1

14、=D=Z9:Z=0:Z2=0:GOTO 2D=Z(Z1-1)8+9:Z=0:Z2=Z1-1:GOTO 2Lbi AZZ18+3=X:ZZ18+4=Y:ZZ18+5=J: ZZ18+8=D: REGU:R“Ro”：E“RN”: G“LS”:U“G”: ZZ18+6=R-1 :ZZ18+7=E-1: ZZ18+9=D+G: ZZ18+10=U: “NEXT” Isz Z1: Goto 0Lbi 1DZO:D：Z：O“RJ”：Z2=0:Goto 2Lbi 2V1=Z2Z1 =GoToDDZZ28+9= A=ZZ28+3:B=ZZ28+4: C =ZZ28+5:R=ZZ28+6: E=ZZ28+7:

15、F=ZZ28+8: G=ZZ28+9: U=ZZ28+10: Goto3Isz Z2:Goto 2Lbi 3W=3 =N=8：N=4PU（E-R）Abs（G-F）：QAbs（D-F）N：S=90Q： JC+(NPQ+2UR)NS：L=1XA+Q6(Cos C+Cos J +4（Cos (C+(L+0.5)PQ+2UR)(L+.5)S),L,0,(N-1)）+2（Cos (C+(LPQ+2UR)LS),L,1,(N-1)+ZCos（J+ O） YB+Q6(Sin C+Sin J +4（Sin (C+(L+0.5)PQ+2UR)(L+.5)S),L,0,(N-1)）+2（Sin (C+(LPQ+2

16、UR)LS),L,1,(N-1)+ZSin（J+ O）:V=1=Goto6V=2=Goto9V=3= GOTO CW=3=GOTO AZ=0=“X”:X:Pause 0: “Y” :Y Goto 1Z“XL”:X:Pause 0: “YL”:Y Goto 1 fx4850 Z0=“XR”:X:Pause 0: “YR”:Y Goto 1 Z=0= X “X” Y “Y” Goto 1Z X “XL” Y “YL” Goto 1 fx4800 Z0= X “XR” Y “YR” Goto 1 Lbi4 MH ：M“X” ：H“Y” ：Z2=0：GOTO 5Lbi5V=1：D= ZZ28+9：Z=

17、0 ：O=90：GOTO 2Lbi6 K=（H -B）Cos（C-90）-（M-A）Sin（C-90）（H -Y）Cos（J-90） -（M-X）Sin（J-90）：K0= Goto 7Isz Z2：Goto5Lbi7 D=F+Abs(H -B)Cos(C-90)-(M-A)Sin（C-90）)：DG= Isz Z2: Goto5Goto 8Lbi8V=2 ：GOTO 3Lbi9 K=(H -Y)Cos(J-90)-(M-X)Sin(J-90)：AbsKGoto BD=D+K ：GOTO 8LbiB V=3 ：Z=0 ：Goto 3LbiCZ=(H-Y) Sin(J+90)：“D”:D:Pau

18、se 0: “Z” Z(Z=(H-Y) Sin(J+90)：D“D” “Z” Z)4800输出GOTO 4LbiD二、说明 a、编制说明 本程序是运用复化辛普生公式根据曲线段直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距，对该曲线段范围内任意里程中边桩坐标进行计算，以及对卡西欧扩充变量的灵活应用，实现了真正意义上的的全线贯通及曲线要素输入程序化（在不修改程序内容的情况下可通过运行程序输入任意多段曲线元要素）。通过对N=？ 进行修改，可对辛普森公式进行任意等分进行运算。（注：N为不小于2的整数，N越大精度越

19、高，计算速度越慢；N越小精度越低，运算速度越快，一般曲线取N=4就能满足精度要求，在能满足精度的情况下尽量N取小值，已获得最佳运算速度，不要盲目的追求精度）b、程序操作说明程序分为两部分：1.JS为计算，2.SZ为设置。1、首先用Shift+Defm键对计算器内存变量进行扩充，扩充变量数为8X+2（X为曲线元段数，变量数视内存情况尽量大些），运行程序，选2进行曲线要素设置，为防止误操作在正确输入口令123456时方可进行设置，否则显示OUT跳到程序尾（LBI4）。在执行程序中将第一段曲线元要素按规律输入到扩充内存变量中（A“X0”:B“Y0”:C“F0”:R“R0”:E“RN”:F“D0”:G

20、“LS”:U“G”），以后曲线仅需输入R“R0”、E“RN” :G“LS”:U“G”，其他参数自动计算得出。在显示NEXT时为提示是否输入下一曲线，按EXE继续输入，如不需再输入则退出程序。2、选1.JS进行坐标运算，D输入桩号，Z输入左右距离（负为左，正为右，0为中）程序首先自动判断其在哪一线元内，并把其线元要素调出进行计算，即可对全线进行坐标计算。c、变量说明X0：Y0：F0曲线元起点X、Y坐标及起点正切线方位角R0：RN曲线元起点及终点半径D0：LS：Q曲线元起点桩号、路线长度及线路左右偏标志（左=-1,直线=0,右=1）D曲线元中待求点桩号Z计算边桩距中线平距，左边输入负值，右边输入正

21、值，中桩输入0RJ ：坐标正算右夹角X：Y（中）XL:YL （左）XR:YR （右）待求点的X，Y坐标扩充变量：ZZ18+3: ZZ18+4: ZZ18+4: ZZ18+6: ZZ18+7: ZZ18+8: ZZ18+9: ZZ18+10: 分别为各线元X0:Y0:F0;R0:RN:D0:LS:G,(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。(2) 当所求点位于中线时，Z=0，坐标显示XY；当位于中线左侧时，Z取负值,坐标显示XLYL,；当位于中线右侧时，Z取正值,坐标显示XRYR。(3) 当线元为直线时，

22、其起点、止点的曲率半径为无穷大，以10的45次方代替。（对于特别长的直线应适当提高10的次方数，以防止以大半径曲线误算，建议以10的90次方计算）(4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。 (5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次方代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次方代替；与圆曲线相接时，曲率半径等于圆曲线的半径。 (6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时， 本程序可在48004850上运行，只

23、需将输出部分进行修改。（48504800）郑石高速路面NO.9标余官营互通E匝道序号XoYoFoRoRND0LSG1744383.112510406.45643。3400.0700070000170-1中点桩号中点X中点Y终点桩号终点X终点Y85744445.055510464.662170744507.700 510522.1122744507.700 510522.112 42。1030.87000500170116.0711中点桩号中点X中点Y终点桩号终点X终点Y228.036744550.193510561.637286.071744590.136510603.7133744590.1

24、36510603.71349。1802.3500500286.071299.2821中点桩号中点X中点Y终点桩号终点X终点Y435.712744669.415510729.969585.353744707.945510873.9874744707.945510873.98783。3544.8500无穷大585.353601中点桩号中点X中点Y终点桩号终点X终点Y615.353744710.545510903.872645.353744712.247510933.823计算时间测试（秒），计算精度满足1mm时测试每段曲线元的终点4等分正算： 2.3 2.4 2.5 2.6反算： 9.612.416.517.05等分正算： 2.5 2.7 2.8 2.9反算：11.813.818.719.1程序特点：全线标高数据一次性程序化输入，扩充变量数据库，无需修改程序内容；全线贯通计算，标高计算时将全线数据分为4等分，首先判断待求点在哪个区间，再循环选择，提高运算时间；傻瓜操作，适用初级用户。一、程序：MCW“1.BG2.SZ ”：W=1= Got