FX5800P全线任意里程中边桩坐标正反算 曲线任意里程中边桩坐标正反算.doc

fx5800p全线任意里程中边桩坐标正反算（加了数据库）本程序非本人原创，结合上传者长处，进行改编，愿收获和大家共享，不足之处多多指教。希望过路神仙把竖曲线加进去带算隧道超欠挖，谢谢曲線任意里程中邊樁坐標正反算5800p本程序由内蒙文哥和师父阿玖哥共同改编完成主程序zbjslbl 6：1.sz xy：2.xy sz：？s：？n：prog“sub0”1pc：(p-r)(2hpr) d：180e:n=1 =goto 1: goto 2lbl 1:?z: abs(s-o)w: progsub1:xs=:xys=:yfs=:f-90 dms“xc” ?c： “yc” ? e：0i：0j：pol（x-c，y-e）:“i=”：iifj0: then “j=”: j+360dmsifendifj 0：then “j=” ：jdmsifend：goto 6lbl 2?x:?y:xi:yj:progsub2:o+ws:cls:s=:sz=:zgoto 6正算子程序 sub10.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(g+qelw(c+lwd)+bsin(g+qefw(c+fwd)+asin(g+qemw(c+mwd) y：g+qew(c+wd)+90f： x+zcos（f）x： y+zsin（f）y反算子程序sub2 g-90t(y-v)cost-(x-u)sin(t) wabs(w)w:0zlbl6：prog sub1t+qew(c+wd) l： (j-y)cos（l）-(i-x)sin（l）zif abs（z）1e-6：then0z:prog sub1:(j-y)sin（f）z:else w+zw：goto6:ifend数据库子程序 sub0goto 1(线元可输入多条，分离式可在前多加一位，匝道一样。例：左幅为k129+500,右幅输线元参数里程为1129+500，其他不变，前面1为任意数字，计算机便于区分)lbl1：ifs线元终点里程：theno(线元起点里程) ： u(线元起点x坐标)：v(线元起点y坐标)：g(线元起点计算方位角)：p(线元起点半径)：r(线元止点半径)：h(线元长度)：q(线元左、右偏标志，左偏-1，右偏1，直线为0)：return：ifendifs xy2. xy = sz s ?选择里程 ,n?，z?左或右偏距。选择计算方式，1.表示进行由里程、边距计算坐标 ；输入2表示由坐标反算里程和边距。显示部分： xs=正算时，计算得出的所求点的x坐标 ys=正算时，计算得出的所求点的y坐标 fs=正算时，所求点对应的中线点的切线方位角xc=测站点x值yc= 测站点y值 i=xxx 测站到放样点距离j=xxx测站到放样点方位角 s=反算时，计算得出的所求点的里程 z=反算时，计算得出的所求点的边距三算例请到/celiang/article_view.asp?id=1101&author=tung 有算例。注:1.本程序根据空间里程序改编而成，并非本人原创。2.正算跨线元计算正确，反算时跨线元计算第一次结果错误，反应速度比一般输出结果慢，给出结果s=xxx为输入坐标的线元概略里程，在进行计算，显示结果*近似里程*，里程及偏距误差0.001m.，可更改数据（1e-6）改为（1e-3），速度稍快一点。3本程序可进行偏角法放样。4.输入每行数据看数据输完毕没有，输完不要忘记按（exe）5.使用者需细心输入，有问题在以下网址查找原因。/celiang/article_view.asp?id=380&author=yshf61e-6即是10的負6次方,亦即是等於0.0000017.希望高人指点把此程序改为数据库形式，并加入高程，希望哪位神仙路过指点一下，把程序改为坐标、高程一体的这程序就完美了，程序中不足之处请指点。曲線任意里程中邊樁坐標正反算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:yjprogsub2:o+ws:clss=:sz=:zgoto 2注:紅色部可以不要,亦可以要,它的作用是將s或x輸入0時,程序會重新由頭開始運行,(即是可以計算下段線元)sub1(子程序名)0.1739274226a:0.3260725774b0.0694318442k:0.3300094782l1-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+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)1e-6:then 0z :progsub1(j-y)sin(f)zelse w+zw:goto 4:ifend注:1e-6即是10的負6次方,亦即是等於0.000001輸入時按shift-log-(-)-6亦可直接將if abs(z)1e-6寫成if abs(z)0.000001注意:本程序是根據yhhf編寫的4800p改寫而成,以下是源程序及算例大家有問題亦可在以下網址查看及向yshf發問/celiang/article_view.asp?id=380&author=yshf曲线任意里程中边桩坐标正反算(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=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)+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 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次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。 (6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。2、输入与显示说明 输入部分：1. sz = xy2. xy = sz n ?选择计算方式，输入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坐标 ys=正算时，计算得出的所求点的y坐标 fs=正算时，所求点对应的中线点的切线方位角 s=反算时，计算得出的所求点的里程 z=反算时，计算得出的所求点的边距四、算例某匝道的由五段线元（直线+完整缓和曲线+圆曲线+非完整缓和曲线+直线）组成，各段线元的要素（起点里程s0、起点坐标x0 y0、起点切线方位角f0、线元长度ls、起点曲率半径r0、止点曲率半径rn、线 元左右偏标志q）如下： s0 x0y0f0ls r0rn q500.00019942.83728343.561 125 16 31.00 269.2561e45 1e45 0769.25619787.34028563.378 125 16 31.00 37.492 1e45 221.75 -1806.74819766.56628594.574 120 25 54.07 112.779221.75 221.75 -1919.52719736.07228701.893 91 17 30.6380.285 221.75 9579.228 -1999.81219744.03828781.659 80 40 50.00100.0001e45 1e450 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=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.47687 ys=28722.35642s=940z= 3.009计算得xs=19733.47298 ys=28722.53168 2、 反算 x=19831.418 y=28509.726计算得s=699.9999974z= -5 .00018164 x=19827.336 y=28506.838计算得s=699.9996493z= 0.000145136 x=19823.25398 y=28503.95084计算得s=699.9999985z= 5.000003137 x=19785.25749 y=28575.02270计算得s=780.0000035z= -5 .000001663 x=19781.15561 y=28572.16358计算得s=780.0000025z=- 0.000002979 x=19777.05373 y=28569.30446计算得s=780.0000016z= 4.99999578 x=19747.536y=28654.131 计算得s=870.0001137z= -4.99941049 x=19742.686y=28652.914 计算得s=870.0003175z=- 0.00041814 x=19737.837y=28651.697 计算得s=870.0002748z= 4.999808656 x=19741.5912 y=28722.0580计算得s=939.9999786z= -5.123024937 x=19736.4769 y=28722.3564计算得s=939.9999862z=- 0.000027710 x=19733.4730 y=28722.5317计算得s=940.0000238z= 3.00898694casio5800直线段中边桩坐标计算（且能计算斜交） lbi0/k(d0)?d:x(d0)?u:y(d0)?v:?f:x(00)?m:y(00)?n:0-i:1-j:1-g/?k:w(l)?a:w?b:n?ck-d-lu+lcosf-xv+lsinf-yf+c-lx(l)=:x-acos(l)-i黑三角y(l)=:y-asin(l)-j黑三角pol(x-acos(l)-m),(y-asin(l)-n)d(l)=:i黑三角jj+360-jf(l)=:j黑三角dms黑三角x(z)=:x黑三角y(z)=:y黑三角pol(x-m),(y-n)d(z)=:i黑三角jj+360-jf(z)=:j黑三角dms黑三角x(r)=:x+bcos(l)-i黑三角y=:y+bsin(l)-j黑三角pol(x+bcos(l)-m,(y+bsin(l)-n)d=:i黑三角jj+360-jf=:j黑三角dms黑三角goto0说明：d-线路起点桩号u-起点x坐标v-起点y坐标f-方位角度m-置仪点xn-置仪点yk-所求桩号a-左侧宽b-右侧宽c-与线路顺时针夹角90或任意角度casio5800计算器测量计算程序 上casio5800计算器测量计算程序 上上月做这个东西的时候没仔细检查，有好几处输错了的地方，今天把它修改过来。简要介绍：1.新版程序把线元法和交点法已经集成在一个模块中了，用户只需修改jd程序和za程序中的数据部分即可，其余不需作任何的改动。2.因为每条路高程计算不尽相同，且比较复杂，现在可利用pc机excel计算好打印成表格带到工地上使用，所以本版程序未对线路高程序进行专门的编程计算，而是利用统计计算模式中来输入桩号（第一列x）及左、右高程（第二、三列y，freq），这种输入数据的方式最为直观，易发现错误，也易修改，输入完毕后运行s程序对数据按桩号进行排序，在程序中通过调用gg程序来进行内插计算，sg=-1得左标高，sg=1得右标高（若sg输入0，则可进行一般的线性内插计算）。3.在jd程序和xy程序中，先将一个计算单元的数据置入矩阵f中（1行8列或1行9列），这样程序可读性极好。4相比原casio4850程序操作习惯，作了一点小小的改动，测站坐标存在z10，n中，x坐标原存在m中容易被误操作修改，而设计标高存在m中，这样易于修改，因为casio5800没有in，out功能，很不方便。4.程序利用z2变量值来判断是采用交点法还是线元法模型计算，z20为线元法，否则为交点法。一、pqx程序：计算中边桩坐标及近似的桩号反算，在运行模式直接调用。z10s:”xo”?s:sz10:”yo”?n:prog“au”lbl2:?l:prog“z”:prog“e”:1no:90slbl4:”jj”?s:”yc”?o:so=0=goto2原来lbl后没有标号的。o=-1=goto6“x，y”:r+ocos(z+s)xu+osin(z+s)yprog“d”:goto4lbl6:z7x:z8y:”xf”?x:”yf”?y:xy=0=goto4xz7:yz8:pol(x-r，y-u+1p):z+s-jj:”yc，dl，l”:icos(j)oisin(j)il+igoto6二、p程序：在程序中提供一个自由运算的模式。lbl1:”tmp”?i:ifi0:then“rst”:igoto1:ifend二、lyc程序：进行桩号反算及边坡放样，在运行模式直接调用。prog“au”lbl1:z7x:z8y:z6s:”xf”?x:xz7:”yf”?y:yz8:”zf”?s:sz6lbl2:prog“z”:y=u=y+1pypol(x-r，y-u):j-zj:isin(j)o:icos(j)iifabs(i)0.1:thenprog“e”:”l，yc”:l+ilogoto3:ifendifz90:thenpol(z9-so，i):jz9180i:ifend”dl”:il+il:goto2lbl3:z6s:ifs=0:thengoto1:ifendmzlbl4:”sg”?z:zm:ifabs(z)=1:thenprog“gg”:yz:ifx=1:thenxz:ifend:zm:ifendifz=0:thenz11c:z12d:z13e:z14f:”la”?c:cz11:”ha”?d:dz12:”lb”?e:ez13:”hb”?f:fz14:d+(e-c)-1(f-d)(l-c)z:zm:ifendlbl5:”sg”:z:”tw”:z-sc:z3p:z4q:z5t:”yc0”?p:pz3:”hc0”?q:qz4:”m”?t:tz5:oabs(o)i:”bl”:i(p+tc(c+q)abs(c)u”db”:iu-iojansi:prog“p”:goto1四、i程序：通过手工方式输要素z2=0=goto2z1r:z2u:”jl”?k:”jx”?g:”jy”?h:”ja”?v:”jj”?w:”r”?a:”ls”?b:”l1”?r:”l2”?u:rz1:uz2:returnlbl2:”l0”?k:”x0”?g:”y0”?h:”a0”?v”:”l1”?f:”l-0z=-1y=1”?w:w=0=return”r0”?a:”r1”?b五、d程序：根据测站点的坐标及放样点的坐标计算水平距离及方位角y=n=y+1nypol(x-z10，y-n):jj+360j”a，d”:jdmsi六、e程序：桩号出界判断z2=0=goto2iflz2:then“l2”:lifendiflz1:then“l1”:lifendgoto3lbl2:iflk:then“f:then“l1”:l-fifendlbl3七、gg程序：用统计数据中桩号，左右高程数据进行内插计算。freqoff:nc:freqon:iflminx:then“lmaxx:then“lmaxx”:c-1i:goto8:ifendfor1itoc-1ifllistxi+1:thenbreak:ifend:nextlbl8:(l-listx)(listxi+1-listx)j:listy+j(listyi+1-listy)x:listfreq+j(listfreqi+1-listfreq)y八、s程序:对统计数据列表中数据按桩号进行排序。freqoff:nc:freqon“waitn=”:locate10，1，cifc2thenreturn:ifendfor1dtoc:di:ford+1etoc:iflistxereturnprog“gg”:”x，y”:xygoto1十一、gzw程序：构造物坐标放样计算，在运行模式直接调用。z10s:”xo”?s:sz10:”yo”?n：0tlbl1;”x0”?r:”y0”?u:”a0”?z:ifz=-1:then“a0，dd”:jziifendifz=0:thenrx:uy:prog“d”:goto1:ifend在后增加一句，解决了程序再次运行时，交角的初值被置为测站点的坐标的bug.o-zs:”j0”?s:z+so:”0-xx，yyx，y”:”1-x，yxx，yy”?t:ift0:thengoto3:ifendlbl2:”xx”?p:”yy”?q:ifpq=0:thengoto1:ifend”x，y”:r+pcos(z)+qcos(o)xu+psin(z)+qsin(o)yprog“d”:goto2lbl3:”xf”?x:”yf”?y:ifxy=0:thengoto1:ifend”xx，yy”:(x-r)sin(o)-(y-u)cos(o)sin(s)p(y-u)cos(z)-(x-r)sin(z)sin(s)qgoto3十二、y程序：线元法计算匝道坐标l-kc:iflf:thenf-kc:ifendifw=0orcf:thenr+rec(l-f，z)r:u+ju:0z9:ifend十三、z程序：交点法计算中桩坐标ifz2=0:thenprog“y”:return:ifendb2a24p:0.1b(5-pa)q:q+(a+p)tan(0.5abs(w)t:k-tc:c+bd:d+aabs(w)180f:f-bewabs(w)s:0z9:iflc:thenl-ko:vz:goto2:ifendiflf:thenl-f+to:v+wz:goto2:ifendifle:thenf-lo:prog”x”:v+wo:o-srz:g+(t-i)cos(o)-sjsin(o)r:h+(t-i)sin(o)+sjcos(o)u:goto3:ifendifld:then(l-d+0.5b)a)ro:q+asin(o)i:p+a(1-cos(o)j:az9:v+soz:goto4:ifendl-co:prog“x”:v+srzlbl4:g+(i-t)cos(v)-sjsin(v)r:h+(i-t)sin(v)+sjcos(v)u:goto3lbl2:g+ocos(z)r:h+osin(z)ulbl3十四、au程序:注：l0-输测段桩号近似值，”au”子程序根据桩号近似值自动调用曲线要素（包括线元法和交点法要素），z2=0则调用线元法数据，否则为交点法数据（l0=0-默认为原曲线要素l0=-1-手工输入曲线要素）cls:”l0=0-conti.”:”l0=-1-input”:ifz2=0:thenlocate13，2，”(xy)”:elselocate13，2，”(jd)”:ifendz2s:”0xye.jd”?s:sz2:”l0”?l:ifl=0:thenreturn:ifendifl=-1:thenprog“i”:return:ifendifz2=0:thenprog“za”:elseprog“jd”:ifendmatf1，1k:matf1，2g:matf1，3h:matf1，4v:matf1，5w:matf1，6a:matf1，7b:ifz2=0:thenmatf1，8f:elsematf1，8z11:matf1，9z12:ifend十五、jd程序：交点法矩阵每一行数据依次为交点桩号，交点x坐标，交点y坐标，起始边方位角，偏角，圆曲线半径，缓和曲线长度，计算单元起点，计算单元终点ifl24900:then“lminx”:lstop:ifendifl30801.27:then“lmaxx”:lstop:ifendifl26227.359:then25676.236，142225.377，504270.949，62459.6，-705010.6，770，300.007，24900，26227.359matf:return:ifendifl27318.276:then27026.458，142873.016，502918.028，2953449.，-153118.5，1500，180.003，26227.359，27318.276matf:return:ifendifl28522.182:then28148.818，143069.516，501810.234，2800330.5，414502.1，780，206.532，27318.276，28522.182matf:return:ifendifl29246.112:then28895.742，143666，501341，3214832.6，-223308.8，1250，215.74，28538.355，29246.112matf:return:ifendifl30801.27:then30431.525，144420，499995，2991523.8，492952.7，720，160.598，29246.112，30801.27matf:return:ifend十六、za程序注2：线元法矩阵每一行数据依次为起点桩号，起点x坐标，起点y坐标，起点方位角，左右转标志，起点曲率半径，终点曲率半径，终点桩号ifl0:then“lminx”:lstop:ifendifl1479.85:then“lmaxx”:lstop:ifendifl1000:thengoto2:ifendifl553.49:then“lmaxx”:lstop:ifendifl82.699:then0，142872.513，502902.111，2900300，-1，1509.5，1509.5，82.699matf:return:ifendifl111.519:then82.699，142898.724，502823.687，2865440，-1，1509.5，1t，111.519matf:return:ifendifl155.022:then111.519，142906.932，502796.061，2862151，-1，1t，1000，155.022matf:return:ifendifl269.54:then155.022，142918.885，502754.233，2850704，-1，1000，1000，269.54matf:return:ifendifl359.475:then269.54，142942.363，502642.211，2783323，-1，1000，1t，359.475matf:return:ifendifl434.421:then359.475，142953.072，502552.924，2755848，1，1t，1200，434.421matf:return:ifendifl463.597:then434.421，142961.655，502478.474，2774609，1，1200，1200，463.597matf:return:ifendifl501.14:then463.597，142965.95，502449.617，2790944，1，1200，1t，501.14matf:return:ifendifl553.488:then501.14，142972.314，502412.618，2800331，0，1t，1t，553.488matf:return:ifendlbl2:ifl1169.959:then1000，142957.323，502388.342，1025507，0，1t，1t，1169.959matf:return:ifendifl1215.026:then1169.959，142919.325，502554，1025507，-1，1t，1500，1215.026matf:return:ifendifl1261.836:then1215.026，142909.471，502597.975，1020328，-1，1500，1500，1261.836matf:return:ifendifl1283.436:then1261.836，142900.408，502643.898，1001611，-1，1500，1t，1283.436matf:return:ifendifl1307.44:then1283.436，142896.66，502665.17，995126，-1，1t，700，1307.44matf:return:ifendifl1346.219:then1307.44，142892.686，502688.842，985230，-1，700，700，1346.219matf:return:ifendifl1363.61:then1346.219，142887.767，502727.303，954203，-1，700，1t，1363