CASIO fx4800P、fx5800P型计算器用于线路施工曲线中线点坐标的计算程序

1、casio fx-4800p、fx-5800p型计算器用于线路施工曲线中线点坐标的计算程序 摘要：本文介绍了casio fx-4800p 、fx-5800p型计算器程序编制用于铁路、公路曲线线路内任意中线点的坐标计算程序及使用方法。本计算程序具有操作简便、计算快捷、应用广泛等特点、极大地减轻了测量工作者的内业工作量，对于测量工作者有较大的参考和指导作用。关键词：曲线 线路 施工测量 计算程序1概述 过去，线路中线施工放样基本依靠经纬仪和钢尺了来进行角度及距离测量。对于曲线线路一般的测量方法是：经纬仪置于某一中线点上，采用偏角法拨角再用钢尺量距来定出中线点。随着电子技术进步和经济发展，测量仪器和

2、测量方法的不断改进，目前，全站仪已广泛地应用于工程施工测量中，极大的提高了测量工作效率。但是，在进行铁路、公路工程的曲线线路施工测设时，需要在线路所在区域建立统一坐标系或独立坐标系，利用坐标变换的方法，将整个曲线的三个部分（第一缓和曲线、中间圆曲线、第二缓和曲线）统一到同一坐标系中。根据坐标系的建立，计算出整个曲线内任意点的坐标，再采用全站仪利用极坐标方法进行施工放样。前提是首先利用计算器计算出各中线点坐标，然后才能进行放样。而普通型计算器不仅计算速度慢，且要求计算者必须正确地记忆很多计算公式，计算繁琐而且容易出错，满足不了现场测设工作的要求。为了能够快速准确地为全站仪提供测设数据，发挥全站仪

3、快速测设的特点，提高测量工作效率，应采用可编程的计算器，编制计算程序。本文主要介绍应用casio fx-4800p型计算器的计算程序，供公司测量同行们参照使用。2计算程序qxzbjs（文件名：曲线坐标计算fx-4800p） defm2:r：l：a：n“zh：x=”：e“zh：y=”：f：“fwj=”：k“zh：lc=”： p=l2/(24r)-l4/(2688r3):m=l/2-l3/(240r2):t“t”=（r+p）tng(a/2)+mg=ra/180：“s”s=g+llbia：c，v：c“csdlc=”：v“hxpj=”：d=c-k：dli=d-d5/（40r2l2）：u=d3/(6rl

4、)-d7/（336r3l3):j=（i2+u2）：goto1：dgo=90（2d-l）/( r)：i=rsino+m：u=r（1-coso）+p：j=（i2+u2）：goto2：d=s-（c-k）： = d-d5/（40r2l2）：z2=d3/6rl-d7/（336r3l3）：i=t+（t-z1）cosa-z2sina：u=（t-z1）sina+z2cosa：j=（i2+u2）：goto3：lbi1:q:q“z=1；y=2”：q=1q=f-30d2/ （rl）：h=f-90d2/ （rl）：q=f+30d2/ （rl）：h=f+90d2/（ rl） goto4：lbi2:q:q“z=1；y=2

5、”：q=1q=f-tng-1(u/i)：h=f-o：q= f+tng-1(u/i)：h=f+ogoto4：lbi3:q:q“z=1；y=2”：q=1q=f-tng-1(u/i):h=f-(a-90(s-(c-k)2/ （rl）): q= f+tng-1(u/i):h=f+(a-90(s-(c-k)2/ （rl）): goto4： lbi4:b=90+h:h0h“qxfwj”=b+360 h360h“qxfwj”=h-360 h“qxfwj”=hgoto5：lbi5: x“csd:x”=jcosq+n+vcosb y“csd:y”=jsinq+e+vsinbgotoa 3程序说明3.1 输入已知

6、变量 r圆曲线半径，显示r？ l缓和曲线长，显示l？ a曲线转向角，显示a？ e直缓点纵坐标，显示zh：x=？ n直缓点横坐标，显示zh：y=？ f第一切线方位角，即zh至jd的方位角，显示fwj=？ k直缓点里程，显示zh：lc=？ 3.2 计算待求量 t切线长度，显示t= s曲线全长，显示s= z3外矢距，即jd到qz的距离，显示e0= 3.3 输入待求变量 k输入待求（测设）点的里程，显示lc=？ v横向偏距，即测设点左、右侧外移距的偏移量，若为中线点输入0；右侧输入“+”值，左侧输入“-”值。3.4 计算结果和说明x待求（测设）点的纵坐标，显示x= y待求（测设）点的横坐标，显示y=

7、该程序既可对曲线内中线任意点的坐标进行计算，又可以对中线外移桩的坐标进行计算。计算原理如下：先建立整个曲线的坐标系统，确定曲线起点（zh）的坐标，以zh到jd的第一切线方位角为计算方位角。 在程序计算过程中，当屏幕显示z=1；y=2？的提示时可根据曲线的转向，从下列两种情况中选择一种予以回答。若曲线为左转，则输入数值1；若曲线为右转，则输入数值2。 本程序式输入时，凡遇到、两种符号时，不应将其分开，应视为一个整体的符号。 由于变量的增加，程序式中使用到z1、z2、z3三个计算量，此时应将计算器的内存变量扩大，另行增加36个变量即可。在计算器上的操作过程为：按shift键再按defm键，然后按数

8、字36即可。3.5 本程序在实际测量工作中的应用 本程序不仅可以计算曲线线路上任意中线点的坐标，亦可以计算线路左、右侧的法线偏距点的坐标，用于曲线中桩和边桩放样。 可以计算曲线线路中心点的切线方位角和法线方位角，该方位角可以用于桥梁墩、台（亦为曲线上桥梁墩、台截面纵、横轴线的方位角）定位。其切线方位角在计算器上显示为“qxfwj= ”。 本程序具有多方面的计算功能，可根据仪器状况和测设对象以及工作需要而灵活应用，可在测设现场计算，减轻内业的计算时间，提高工作效率。4计算举例 结合某线路的曲线资料，计算曲线线路上中桩或边桩的坐标。4.1已知要素 r=1500， l=150， az=30 主点里程

9、：zhdk8+150.000 hydk8+300.000 yhdk8+935.398 hzdk9+085.398 4.2 建立平面控制坐标系 建立独立坐标系统，若以zh为坐标原点，以zhjd的切线为x轴的正方向，则zh点的坐标可设为：x=8150.000，y=1000.000；切线的初始方位角f“fwj= ”可设为00000。现以计算圆曲线内dk8+400的中线点或左侧15米边桩点的坐标为例。4.3 计算操作步骤 1输入曲线半径r=1500按exe键 2输入缓和曲线长l=150按exe键 3输入曲线转向角a=30按exe键 4输入zh点x坐标zh：x=8150.000按exe键 5. 输入zh

10、点y坐标zh：y=1000.000按exe键 6输入切线方位角fwj=00000按exe键 7输入zh点里程zh：lc=8150按exe键 显示t=477.0849919（切线长）按exe键；显示s=935.398164（曲线长）；按exe键显示e0=53.56126046（外矢距）；按exe键8输入测设点里程lc=8400；按exe键9输入中线点或左侧边桩的距离，如计算中桩点坐标输入pju=0，如需计算左侧15米边桩的坐标则输入pju=-15米；按exe键10屏幕显示z=1；y=2？提示，本例为左偏曲线故输入1，按exe键,显示fxfwj=83.31549239，按shift键再按，即dk8

11、+400中桩点的法线方位角为831855.77按exe键11. 显示csd：x=8399.5970289，按exe键； 显示csd：y=989.1782961，按exe键。若需计算dk8+400左侧15米边桩点的坐标，按exe键，屏幕提示lc=8400；再按exe键，提示pju=？，输入-15米，再按exe键，显示z=1；y=2？输入数字1按exe键，显示qxfwj=83.31549239，按exe键；显示csd：x=8397.850997，按exe键；显示csd：y=974.2802637，按exe键。12若需要再计算其他线路中线点或边桩点的坐标，重复811步骤操作，不必再输入各已知数据，仅

12、输入需要计算点的里程或横向偏距及z=1；y=2？的参数即可。zbfs(坐标正算)a”x1”:b”y1”:c”x2”:d”y2”:fixm：pol（c-a,d-b）:i”s=” jgoto2l-dt=h=m+ (l-d)f +（f-c）abs(f-c) (t-l +d)2(2r)：goto3：h=m+ (l-d)f：goto3lbi2：d-lt=h=m-（d-l）c+（f-c） abs(f-c)（t-d+ l）2(2r)：goto3：h=m-（d-l）c：goto3lbi3：hgoto1三、程序执行1、ci1 前进方向输入变坡点后段坡率2、fi2 输入变坡点前段坡率以路线前进方向3、r 输入竖曲

13、线半径。计算出t竖曲线切线长，e竖曲线外矢距。4、l 输入所求点的桩号（里程）5、dbp.lc 输入变坡点桩号（里程）6、mbp.gc 输入变坡点高程。7、h 以路线计算出所求点高程四、说明1、计算凸曲线，所求点桩号（里程）不小于后一个凹曲线临界（终点）桩号，不大于前一个凹曲线临界（起点）桩号。2、计算凹曲线，所求点桩号（里程）不小于后一个凸曲线临界（终点）桩号，不大于前一个凸曲线临界（起点）桩号。曲线任意里程中边桩坐标正反算(casio fx-4800p计算器)程序 一、程序功能 本程序由一个主程序(tyqxjs)和两个子程正算子程序(sub1)、反算子程序(sub2)序构成，可以根据曲线段

14、直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。另外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲线任意里程中边桩坐标进行正反算。本程序也可以在casio fx-4500p计算器及 casio fx-4850p计算器上运行。特别申明：(1). 适用于弧长小于2倍半径的各种线元坐标正反算，精度优于1mm； (2). 在引用该核心计算部分时，请注明来源。二、源程序1.主程序(tyqxjs)1.sz = xy：2.xy = sz：n：

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

16、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+qew(c+wd)：z=(j-y)cosl-(i-x)sinl：abszgoto1：w=w+z：goto 0lbl 1：z=0：prog

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

18、，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。 (6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。2、输入与显示说明 输入部分： 1. sz = xy 2. xy = sz n ? 选择计算方式，输入1表示进行由里程、边距计算坐标 ；输入2表示由坐标反算里程和边距。 x0 ？线元起点的x坐标 y0 ？线元起点的y坐标 s0 ？线元起点里程 f0 ？线元起点切线方位角 ls ？线元长度 r0 ？线元起点曲率

19、半径 rn ？线元止点曲率半径 q ？ 线 元左右偏标志(左偏q=-1，右偏q=1，直线段q=0） s ？ 正算时所求点的里程 z ？正算时所求点距中线的边距(左侧取负，值右侧取正值，在中线上取零) x ？反算时所求点的x坐标 y ？反算时所求点的y坐标 显示部分： xs= 正算时，计算得出的所求点的x坐标 ys= 正算时，计算得出的所求点的y坐标 fs= 正算时，所求点对应的中线点的切线方位角 s= 反算时，计算得出的所求点的里程 z= 反算时，计算得出的所求点的边距四、算例 某匝道的由五段线元（直线+完整缓和曲线+圆曲线+非完整缓和曲线+直线）组成，各段线元的要素（起点里程s0、起点坐标x

20、0 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 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.

21、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 计算得 xs=19827.33592 ys=28506.83837 s=700 z= 5 计算得 xs=19823.25398 ys=28503.95084 s=780 z=-5 计算得 xs=19785.25749 ys=2

22、8575.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=870 z= 5 计算得 xs=19737.83688 ys=28651.69668 s=940 z=-5.123 计算得 xs=19741. 59118 ys=28722.05802 s=940 z=0 计算得 xs

23、=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.000145136 x=19823.25398 y=28503.95084 计算得 s=699.9999985 z= 5.000003137 x=19785.25749 y=28575.02270 计算得 s=780.000

24、0035 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.99941049 x=19742.686 y=28652.914 计算得 s=870.0003175 z=- 0.00041814 x=19737.837 y=28651.697 计算得 s=870.0002748

25、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（dmcl）断面测量（fx -4800）abcdekifgj：a：b：c”r1” ：d”r2” ：e”r3” ：k”h1” ：i：f”h0” ：g”s0” ：j”lc”lbl 1：hz：h：z”pj” ：x=f+(j-g)i+k：y=x+c/2

26、：w=x+a：u=w-dsin（175355）：q=w-(d-e) sin（175355）：hx+c=goto 3hy=p=(c2-(h-x)2):goto 2hu=p=b+(d2-(h-w)2):goto 2: p=b+(d-e) cos（175355）+(e2-(h-q)2) :goto 2lbl 2：“z0=“:z-pgoto 1lbl 3：“h=“:h-(x+c)goto 1jdzb（文件名：曲线坐标计算casio fx-5800p）“jdx1”?a:“jdy1”?b: “jdx2”?c:“jdy2”?d:pol(c-a,d-b):r：l：a：jf：“zh：lc=”?k：l2/24r-

27、l4/2688r3p:l/2-l3/240r2m：（r+p）*tng(a/2)+mt: “t=”:tra/180g： g+ls: “s=”:sc+tcos(f+180) n: “zh.x=”:nd+tsin(f+180) e: “zh.y=”:e lbia：c，v：c“dk?+”：v“pj”：c-kd：if dl:then d-d5/（40r2l2) i: d3/6rl-d7/（336r3l3）u:（i2+u2）j：goto1：ifend：if dg:then o=90（2d-l）/ （r）： rsino+mi：r（1-coso）+pu：（i2+u2）j：goto2：else s-（c-k）d

28、:d-d5/（40r2l2）z1： d3/6rl-d7/（336r3l3）z2：t+（t-z1）cosa-z2sinai：（t-z1）sina+z2cosau：（i2+u2）j：goto3：ifendlbi1:“z=1；y=2”？x：if x=1 then f-30d2/ （rl）q：f-90d2/ （rl）h：else f+30d2/ （rl）q：f+90d2/（ rl）h:ifend: goto4：lbi2: x“z=1；y=2”：if x=1 then f-tng-1(u/i) q： f-oh：else f+tng-1(u/i) q：f+oh: ifend:goto4：lbi3: x“z

29、=1；y=2”：if x=1 then f-tng-1(u/i) q:f-(a-90(s-(c-k)2/ （rl）) h: else f+tng-1(u/i) q: f+(a-90(s-(c-k)2/ （rl）) h: ifend: goto4： lbi4: 90+hb: if h0 then b+360h: ifend: if h360 then h-360h : else hh: ifend:goto5：lbi5:“x=”:jcosq+n+vcosbx “y=”:jsinq+e+vsinb“fwj=”:hgotoa 3程序说明3.1 输入已知变量“jdx1” 交点1x坐标“jdy1” 交点

30、1y坐标 “jdx2” 交点2x坐标“jdy2” 交点2y坐标 r圆曲线半径，显示r？ l缓和曲线长，显示l？ a曲线转向角，显示a？ e直缓点纵坐标，显示zh：x= n直缓点横坐标，显示zh：y= f第一切线方位角，即zh至jd的方位角 k直缓点里程，显示zh：lc=？ 3.2 计算待求量 t切线长度，显示t= s曲线全长，显示s= 3.3 输入待求变量 k输入待求（测设）点的里程，显示dk?+？ v横向偏距，即测设点左、右侧外移距的偏移量，若为中线点输入0；右侧输入“+”值，左侧输入“-”值。3.4 计算结果和说明x待求（测设）点的纵坐标，显示x= y待求（测设）点的横坐标，显示y= 该程

31、序既可对曲线内中线任意点的坐标进行计算，又可以对中线外移桩的坐标进行计算。计算原理如下：先建立整个曲线的坐标系统，确定曲线起点（zh）的坐标，以zh到jd的第一切线方位角为计算方位角。 在程序计算过程中，当屏幕显示z=1；y=2？的提示时可根据曲线的转向，从下列两种情况中选择一种予以回答。若曲线为左转，则输入数值1；若曲线为右转，则输入数值2(非1也对)。3.5 本程序在实际测量工作中的应用 本程序不仅可以计算曲线线路上任意中线点的坐标，亦可以计算线路左、右侧的法线偏距点的坐标，用于曲线中桩和边桩放样。 可以计算曲线线路中心点的切线方位角和法线方位角，该方位角可以用于桥梁墩、台（亦为曲线上桥梁

32、墩、台截面纵、横轴线的方位角）定位。其切线方位角在计算器上显示为“fwj= ”。 本程序具有多方面的计算功能，可根据仪器状况和测设对象以及工作需要而灵活应用，可在测设现场计算，减轻内业的计算时间，提高工作效率。4计算举例 结合某线路的曲线资料，计算曲线线路上中桩或边桩的坐标。4.1已知要素 “jdx1” 交点1x坐标“jdy1” 交点1y坐标 “jdx2” 交点2x坐标“jdy2” 交点2y坐标 r圆曲线半径，显示r？ l缓和曲线长，显示l？ a曲线转向角，显示a？ e直缓点纵坐标，显示zh：x= n直缓点横坐标，显示zh：y= f第一切线方位角，即zh至jd的方位角 k直缓点里程，显示zh：

33、lc=？ “jdx1” =2920660.136“jdy1” =475778.751 “jdx2” =2919624.99“jdy2” =477022.978r=4500， l=340， az=262423”(左转) 主点里程：zh.lc=dk348+692.554 4.3 计算操作步骤 1输入交点1的x坐标 “jdx1” =2920660.136按exe键 2输入交点1的y坐标 “jdy1” =475778.751按exe键3输入交点2的x坐标 “jdx2” =2919624.990按exe键 4输入交点2的y坐标 “jdy2” =477022.978按exe键5. 输入曲线半径r=4500

34、按exe键 6输入缓和曲线长l=340按exe键 7输入曲线转向角a=262423”按exe键 8输入交点2的zh点里程zh：lc=348692.554按exe键 显示交点2的t=1225.972（切线长）按exe键；显示交点2的s=2413.952（曲线长）；按exe键e直缓点纵坐标，显示交点2的zh：x=2920409.074 n直缓点横坐标，显示交点2的zh：y=476080.52279输入交点2的测设点里程dk?+?=349032.554；按exe键10输入交点2中线点或左侧边桩的距离，如计算中桩点坐标输入pju=0，如需计算左侧15米边桩的坐标则输入pju=-15米；按exe键11屏

35、幕显示z=1；y=2？提示，本例为左偏曲线故输入1，按exe键12. 显示x=2920194.946，按exe键； 显示y=476344.5952，按exe键。显示qxfwj=1273540.43”， 即dk349+032.554中桩点的切线方位角为1273540.43”按exe键若需计算dk349+032.554左侧15米边桩点的坐标，按exe键，屏幕提示dk?+?= 349032.554；再按exe键，提示pju=？，输入-15米，再按exe键，显示z=1；y=2？输入数字1按exe键， 显示x=2920206.831，按exe键；显示476353.746，按exe键。显示qxfwj=12

36、73540.43”，按exe键；13若需要再计算其他线路中线点或边桩点的坐标，重复811步骤操作，不必再输入各已知数据，仅输入需要计算点的里程或横向偏距及z=1；y=2？的参数即可。竖曲线程序fx-5800p一、程序名：sqx二、程序:i1?c：i2? f：r：bp.lc? d：bp.gc? m：: 0.5abs(f-c)rtt2/(2r) elbi1：“dk？+”? l：if ld：then goto2：ifend：if l-dt：then m+ (l-d)f +（f-c）abs(f-c) (t-l +d)2(2r) h：goto3：else m+ (l-d)fh：goto3：else m+

37、(l-d) f h：goto3：ifend：lbi2：if d-lt：then m-（d-l）c+（f-c） abs(f-c)（t-d+ l）2(2r) h：goto3：else m-（d-l）ch：goto3：ifend：lbi3：”h=” ：hgoto1三、程序执行1、ci1 前进方向输入变坡点后段坡率2、fi2 输入变坡点前段坡率以路线前进方向3、r 输入竖曲线半径。计算出t竖曲线切线长，e竖曲线外矢距。4、l 输入所求点的桩号（里程）5、dbp.lc 输入变坡点桩号（里程）6、mbp.gc 输入变坡点高程。7、h 以路线计算出所求点高程四、说明1、计算凸曲线，所求点桩号（里程）不小于后

38、一个凹曲线临界（终点）桩号，不大于前一个凹曲线临界（起点）桩号。2、计算凹曲线，所求点桩号（里程）不小于后一个凸曲线临界（终点）桩号，不大于前一个凸曲线临界（起点）桩号。曲线任意里程中边桩坐标正反算(casio fx-5800p计算器)程序 一、程序功能 本程序由一个主程序(tyqxjs)和两个子程正算子程序(sub1)、反算子程序(sub2)序构成，可以根据曲线段直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。另外也可以将本程序中核心算法部分的两个子程

39、序移植到其它相关的程序中，用于对曲线任意里程中边桩坐标进行正反算。本程序也可以在casio fx-4500p、casio fx-4800p计算器及 casio fx-48500p计算器上运行。 特别申明：(1). 适用于弧长小于2倍半径的各种线元坐标正反算，精度优于1mm； (2). 在引用该核心计算部分时，请注明来源。 二、源程序1.主程序(tyqxjs)1.sz = xy：2.xy = sz：？n： x0？u： y0？v： s0？o： f0？g： ls？h： r0？p： rn？r：？q： 1pc： (p-r)(2hpr) d： 180e：if n=1：then goto 1：else go

40、to 2：ifend lbl 1：? s：?z： abs(s-o) w：if wh or sh or w0 ：then stop ：ifend： o+ws： s=：sz：zf-90f： fs=：fgoto 22. 正算子程序(sub1)0.1739274226a： 0.3260725774b： 0.0694318442k： 0.3300094782l： 1-lf： 1-km： u+w(acos(g+qekw(c+kwd)+bcos(g+qelw(c+lwd)+bcos(g+qefw(c+fwd)+acos(g+qemw(c+mwd) x：v+w(asin(g+qekw(c+kwd)+bsin(

41、g+qelw(c+lwd)+bsin(g+qefw(c+fwd)+asin(g+qemw(c+mwd) y： g+qew(c+wd)+90f： x+zcosfx： y+zsinfy3. 反算子程序(sub2)g-90t： abs(y-v)cost-(x-u)sint) w： 0z：lbl 0：prog sub1： t+qew(c+wd) l： (j-y)cosl-(i-x)sinlz：if absz xy 2. xy = sz n ? 选择计算方式，输入1表示进行由里程、边距计算坐标 ；输入2表示由坐标反算里程和边距。 x0 ？线元起点的x坐标 y0 ？线元起点的y坐标 s0 ？线元起点里程

42、f0 ？线元起点切线方位角 ls ？线元长度 r0 ？线元起点曲率半径 rn ？线元止点曲率半径 q ？ 线 元左右偏标志(左偏q=-1，右偏q=1，直线段q=0） s ？ 正算时所求点的里程 z ？正算时所求点距中线的边距(左侧取负，值右侧取正值，在中线上取零) x ？反算时所求点的x坐标 y ？反算时所求点的y坐标 显示部分： xs= 正算时，计算得出的所求点的x坐标 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.340 28563.378 125 16 31.00 37.492 1e45 221.75 -1806.748 19766.566 28594.574 120 25 54.07 112.779