5800计算器程序.docx

1、5800计算器程序曲線任意里程中邊樁坐標正反算5800pTYQXJS(主程序名)LbI 0:4Dimz:Cls1.SZXY2.XYSZ?N:ClsXO?U: YO?V: SO?O:FO?G: LS?H: RO?P:RN?R:?Q1PC：(P-R)(2HPR) D：180EN=2=Goto 2LbI 1:Cls:?SS=0=Goto 0(注:當S輸入0時,程序會由LbI 0開始重新運行)?Z:Abs(S-O)WProgSUB1:ClsXS=:XYS=:YFS=:F-90 DMSGoto 1LbI 2:Cls:?XX=0=Goto 0(注:當X輸入0時,程序由會LbI 0開始重新運行)?Y:XI:

2、YJProgSUB2:O+WS:ClsS=:SZ=:ZGoto 2注:紅色部可以不要,亦可以要,它的作用是將S或 X輸入0時,程序會重新由頭開始運行,(即是可以計算下段線元)SUB1(子程序名)0.A:0.B0.K:0.L1-LF:1-KMG+QEKW(C+KWD)Z1G+QELW(C+LWD)Z2G+QEFW(C+FWD)Z3G+QEMW(C+MWD)Z4Acos(Z1)X(剩號可省略)X+Bcos(Z2)XX+Bcos(Z3)XX+Acos(Z4)XU+WXXAsin(Z1)YY+Bsin(Z2)YY+Bsin(Z3)YY+Asin(Z4)YV+WYYG+QEW(C+WD)+90FX+Zc

3、os(F)XY+Zsin(F)YSUB2(子程序名)G-90T(Y-V)cos(T)-(X-U)sin(T)WAbs(W)W:0ZLbI 4:ProgSUB1T+QEW(C+WD)L(J-Y)cos(L)-(I-X)sin(L)ZIf Abs(Z)1E-6:Then 0Z :ProgSUB1(J-Y)sin(F)ZElse W+ZW:Goto 4:IfEnd注:1E-6即是10的負6次方,亦即是等於0.輸入時按SHIFT-log-(-)-6亦可直接將If Abs(Z)1E-6寫成If Abs(Z)0.注意:本程序是根據yhhf編寫的4800P改寫而成,以下是源程序及算例大家有問題亦可在以下網址

4、查看及向yshf發問曲线任意里程中边桩坐标正反算(CASIO fx-4800P计算器)程序一、程序功能本程序由一个主程序(TYQXJS)和两个子程正算子程序(SUB1)、反算子程序(SUB2)序构成，可以根据曲线段直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。另外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲线任意里程中边桩坐标进行正反算。本程序也可以在CASIO fx-4500P计算器及 CASIO fx-4850P计算器上

5、运行。 特别申明：(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

6、.：B=0.：K=0.：L=0.：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+

7、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) 当线元为完

8、整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。 (6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。2、输入与显示说明 输入部分：1. SZ = XY2. XY = SZ N ?选择计算方式，输入1表示进行由里程、边距计算坐标 ；输入2表示由坐标反算里程和边距。 X0 ？线元起点的

9、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=反算时，计算得出的所求点的边距四、算例某匝道的由五段线

10、元（直线+完整缓和曲线+圆曲线+非完整缓和曲线+直线）组成，各段线元的要素（起点里程S0、起点坐标X0 Y0、起点切线方位角F0、线元长度LS、起点曲率半径R0、止点曲率半径RN、线 元左右偏标志Q）如下：S0 X0 Y0 F0LS R0 RN Q500.00019942.83728343.561 125 16 31.00269.256 1E45 1E450769.25619787.34028563.378 125 16 31.0037.492 1E45 221.75 -1806.74819766.56628594.574 120 25 54.07 112.779221.75 221.75 -

11、1919.52719736.07228701.893 91 17 30.63 80.285221.75 9579.228 -1999.81219744.03828781.659 80 40 50.00 100.0001E45 1E45 0 1、正算（注意：略去计算方式及线元要素输入，请自行根据所求点所在的线元输入线元要素）S=700Z=-5计算得XS=19831.41785 YS=28509.72590S=700Z=0 计算得XS=19827.33592 YS=28506.83837S=700Z= 5计算得XS=19823.25398 YS=28503.95084S=780Z=-5计算得XS=

12、19785.25749 YS=28575.02270S=780Z=0 计算得XS=19781.15561 YS=28572.16358S=780Z= 5计算得XS=19777.05373 YS=28569.30446 S=870Z=-5计算得XS=19747.53609 YS=28654.13091S=870Z=0 计算得XS=19742.68648 YS=28652.91379S=870Z= 5计算得XS=19737.83688 YS=28651.69668 S=940Z=-5.123计算得XS=19741. 59118YS=28722.05802S=940Z=0计算得XS=19736.47

13、687 YS=28722.35642S=940Z= 3.009计算得XS=19733.47298 YS=28722.53168 2、 反算 X=19831.418 Y=28509.726 计算得S=699.Z= -5 . X=19827.336 Y=28506.838 计算得S=699.Z= 0. X=19823.25398 Y=28503.95084计算得S=699.Z= 5. X=19785.25749 Y=28575.02270计算得S=780.Z= -5 . X=19781.15561 Y=28572.16358 计算得S=780.Z=- 0. X=19777.05373 Y=2856

14、9.30446 计算得S=780.Z= 4. X=19747.536 Y=28654.131 计算得S=870.Z= -4. X=19742.686 Y=28652.914 计算得S=870.Z=- 0. X=19737.837 Y=28651.697 计算得S=870.Z= 4.X=19741.5912 Y=28722.0580 计算得S=939.Z= -5.X=19736.4769 Y=28722.3564 计算得S=939.Z=- 0.X=19733.4730 Y=28722.5317 计算得S=940.Z= 3.谢谢发贴：+24 分 【测量空间】0查看楼主其它文章Goto 2LbI 1

15、:Cls:?SS=0=Goto 0?Z:Abs(S-O)WProgSUB1:ClsXS=:XYS=:YFS=:F-90 DMS (90与DMS有空格?)Goto 1LbI 2:Cls:?XX=0=Goto 0?Y:XI:YJProgSUB2:O+WS:ClsS=:SZ=:ZGoto 2SUB1(子程序名)0.A:0.B0.K:0.L1-LF:1-KMG+QEKW(C+KWD)Z1G+QELW(C+LWD)Z2G+QEFW(C+FWD)Z3G+QEMW(C+MWD)Z4Acos(Z1)XX+Bcos(Z2)XX+Bcos(Z3)XX+Acos(Z4)XU+WXXAsin(Z1)YY+Bsin(Z

16、2)YY+Bsin(Z3)YY+Asin(Z4)YV+WYYG+QEW(C+WD)+90FX+Zcos(F)XY+Zsin(F)YSUB2(子程序名)G-90T(Y-V)cos(T)-(X-U)sin(T)WAbs(W)W:0ZLbI 4:ProgSUB1T+QEW(C+WD)L(J-Y)cos(L)-(I-X)sin(L)ZIf Abs(Z)10（-6）:Then 0Z ProgSUB1(Z与ProgSUB1之间的空格?)(J-Y)sin(F)ZElse W+ZW:Goto 4:IfEnd6 楼：tung 楼主 09-2-8 15:34 回所有程序都沒有空格,一直輸入就可以,轉行別忘記按EX

17、E鍵,本人輸入后,計算所得答案与yshf提供的答案一樣=-是一 內置符號.輸入按FUNCTION 3 REPLAY(向下鍵二下) 3大家有問題亦可在以下網址查看及向yshf發問7 楼：it8752009-2-15 20:18 回评：+6 分很好,很简洁,辛苦啦,非常感谢8 楼：hong09-2-19 22:49 回评：+10 分楼主，哪些地方要转行的呢?9 楼：hong09-2-20 0:19 回评：+10 分楼主，怎么我到Q就错了呢？1. SZ = XY2. XY = SZ N ?选择计算方式，输入1表示进行由里程、边距计算坐标 ；输入2表示由坐标反算里程和边距。 X0 ？线元起点的X坐标

18、Y0 ？线元起点的Y坐标 S0 ？线元起点里程 F0 ？线元起点切线方位角 LS ？线元长度 R0 ？线元起点曲率半径 RN ？线元止点曲率半径到了 Q?就错了10 楼：鸿09-2-20 0:25 回评：+10 分楼主，你能用电话帮我下吗？我真的不知道错在哪，你发条信息来，我给你回打电话。谢谢了11 楼：geergril09-2-20 18:41 回评：+10 分输入变量运算正确,但强烈建议楼主将其改编为数据库型的最好,期待楼主出新版本!12 楼：tung 楼主 09-2-20 19:50 回10楼:鸿建议再检查清楚13 楼：it8752009-2-20 23:24 回楼主,请问一下,为什么我

19、反算的时侯桩号会对而距离不对?还有就是该语句If Abs(Z)10（-6）:Then 0Z ProgSUB1后面的0Z：ProgSUB1是否要加冒号？14 楼：11109-2-22 9:29 回评：+4 分楼主，我为什么正算行反算就不行了 那里有问题请你仔细检查检查15 楼：11109-2-22 9:30 回评：+4 分楼主，我为什么正算行反算就不行了 那里有问题请你仔细检查检查16 楼：tung 楼主 09-2-22 18:55 回13楼:it87520要加冒号的,多谢提示,已更正If Abs(Z)10（-6）:Then 0Z :ProgSUB117 楼：tung 楼主 09-2-22 20

20、:37 回14楼:111本人初次输入后亦有问题,经多次检查后才能使用运行程序时程序如显示Syntax ERROR,这表示程序中欠缺了某种符号,例如冒号.变量(英文字母)等这时你要按AC键或者EXIT键看下光标停在什么位置,并检查光标附近句子是缺少了什么符号18 楼：hanhan09-2-22 19:55 回评：+8 分请问楼主，我按照你上面程序编写还是不能进行计算我编写完成后，输入第一个命令就是选择正算还是反算（我那个计算器上面没有出现N这个符号，不知道是什么原因）19 楼：tung 楼主 09-2-22 20:35 回18楼:hanhan没错,N这个变量是不会显示的,应该显示如下:1.SZX

21、Y2.XYSZ?执行SZXY时-输入1 EXE执行XYSZ时-输入2 EXE20 楼：11109-2-23 9:13 回评：+10 分对了 谢谢21 楼：hanhan09-2-23 9:37 回评：+8 分请问楼主，整个程序我照你上面编写，执行还是不成功哦不知道是什么原因哦（执行程序之后出现：GOERROR22 楼：hanhan09-2-23 10:37 回评：+8 分请问楼主，程序里面0和o如何区分23 楼：geergirl09-2-23 15:42 回这个程序手输线元参数后计算结果是正确的,但不是数据库,感觉不好24 楼：hanhan09-2-23 17:47 回听不懂，数据库程序跟这个有

22、区别吗？可以讲解一下吗？25 楼：tung 楼主 09-2-23 20:17 回LbI 0X=0=Goto 0 以上的是零(比较扁瘦)XO?U: YO?V: SO?O:FO?G: LS?H: RO?P:RN?R:?Q以上的是英文字母O(比较园肥)26 楼：geergirl09-2-23 20:33 回评：+6 分输入变量运算正确,但强烈建议楼主将其改编为数据库型的最好能全线连续计算,期待楼主出新版本!27 楼：包西09-2-23 21:47 回楼主你好，可以给我说说在那里可以找到这几个固定值吗？谢谢了，0.A:0.B0.K:0.L 28 楼：vbvbvbvb09-2-23 22:32 回楼主：这个5800程序，怎样加入数据库？4850已解决！29 楼：liugang09-2-25 18:0 回这个程序4800版本我一直在用（已经2年了），现在出来了5800版本，真高兴啊，只是给楼主提个问题：应用中发现反算程序计算里程和偏移值相差太大，只有把反算程序运行两次差值才可以达到精度，请问楼主是否注意到，能否改进30 楼：tung 楼主 09-2-25 19:59 回27.28.28搂:这个程序不是我编写的,是改编自本论坛yshf的4800p大家有問題亦可在以下網址查看及向yshf發問