测量天地

设为首页  
加入收藏

CSAIO计算器全自动任意交角中边桩坐标计算程序
时间:2017-04-06 17:04:42,点击:0
一、平曲线测量简单说明:

  单个程序,单独运行,专为CASIO4500而改编,亦可用于4800,需按输入习惯稍作调整。此程序能用于计算对称完整曲线(直—缓—圆—缓—直)、单圆曲线(直—圆—直)、两缓和曲线不相等的复曲线(直—缓1—圆—缓2—直)和无圆曲线的两缓和曲线相连的曲线(直—缓—缓—直)。
程序名:【H.CURVE】 L01:LbI 0………………………………………………………………………………………将该处行标设为0
L02:Norm………………………………………………………………………………………将计算器设置为普通模式
L03:Deg………………………………………………………………………………………将计算器设置为度分秒格式
L04:{QW}:RL“LSci1”M“LSci2”F“JD”E“JDX”N“JDY”B“PJ”A“AZIMUTH”C“R.ANG”K“STAT”Q“PEG NO.”W“WIDTH”…………………………………………将变量赋予Q、W等变量并显示
L05:QGoto0………………………………………………………………如果桩号超出上一交点的HZ点则转到执行0
L06:C=90=>GotoX
L07:C“ANG=”………………………………………………………………L5~L6如交角不是正交则显示交角值
L08:LbIX……………………………………………………………………………………将该处行标设为X
L09:G=AbsBπR/180-(L+M)/2……………………………………………………计算圆曲线长
L10:S=M:Z=1………………………………………………………………………………将M(Ls2)值赋予S以计算p2和q2
L11:LbIY……………………………………………………………………………………将该处行标设为Y
L12:P=S2/24R-Sxy4/2688Rxy3………………………………………………………计算内移值
L13:O=S/2-Sxy3/240R2…………………………………………………………………计算切线增值
L14:Z=1=>Z=0:S=L:U=P:V=O:GotoY…将算得的p2和q2分别存入U、V,将L(Ls1)值赋予S以计算p1和q1
L15:T=(R+U-(R+P)cosAbsB)/sinAbsB+O…………………………………………计算切线长T(T1)
L16:S=(R+P-(R+U)cosAbsB)/sinAbsB+V…………………………………………计算切线长S(T2)
L17:D=F-T+0.001………………………………………………………………………………计算ZH点桩号
L18:Q>D+M+L+G=>Goto0
L19:Q≤D=>Goto1 L18~L22判断要求的桩号所在
L20:Q≤D+L=>J=Q-D:U=L:Goto2 的位置是在直线段还是在第一
L21:Q≤D+L+G=>J=1:Goto3 缓和段或圆曲线段或第二缓和
L22:Q≤D+M+L+G =>J=Q-D-M-L-G:U=M:Goto2 段并作出转换
L23:LbI1…………………………………………………………………………………将该处行标设为1
L24:X=E+(F-Q)cos(A+180)………………………………………………………计算直线段的中桩X坐标
L25:Y=N+(F-Q)sin(A+180) ……………………………………………………计算直线段的中桩Y坐标
L26:S=A+C………………………………………………………………………………直线段边桩的方位角 L27:Goto7……………………………………………………………………………… 转换到7
L28:LbI2…………………………………………………………………………………将该处行标设为2
L29:H=J-Jxy5/40R2U2+Jxy9/3456Rxy4Uxy4………………………………计算缓和段横支距
L30:I=Jxy3/6RU-Jxy7/336Rxy3Uxy3+Jxy11/42240Rxy5Uxy5………计算缓和段纵支距
L31:J=Q-D=>Goto4………………………………………………………………若为第一缓和段则转到4
L32:Goto5………………………………………………………………………………若为第二缓和段则转到5
L33:LbI3…………………………………………………………………………………将该处行标设为3
L34:U=90(2Q-2D-L)/πR…………………………………………………为所求点到ZH(HZ)点的曲线长所对应的圆心角
L35:H=RsinU+O………………………………………………………………计算圆曲线段到HZ点横支距
L36:I=R(1-cosU)+P…………………………………………………………计算圆曲线段到HZ点纵支距
L37:LbI4…………………………………………………………………………将该处行标设为4
L38:B<0=>I=-I………………………………………………………………曲线左偏则纵支距为负 L39:X=E+Tcos(A+180)+HcosA+Icos(A+90) 计算第一缓和段和圆曲
L40:Y=N+Tsin(A+180) +HsinA+Isin(A+90) 线的中桩X、Y坐标
L41:Goto6………………………………………………………………………转到6去计算边桩
L42:LbI5………………………………………………………………………将该处行标设为5,计算第二缓和段中桩X、Y坐标
L43:B>0=>I=-I…………………………………………………………………与L36相反,曲线右偏则纵支距为负
L44:X=E+Scos(A+B)+Hcos(A+B)+Icos(A+B+90) ………………计算第二缓和段的中桩X坐标
L45:Y=N+Ssin(A+B)+Hsin(A+B)+Isin(A+B+90) ………………计算第二缓和段的中桩Y坐标
L46:LbI6……………………………………………………………………………将该处行标设为6
L47:U=R……………………………………………………………………………将R用U代替,以免构成错误循环
L48:B<0=>U=-R………………………………………………………………曲线左偏则半径为负
L49:J=Q-D=>S=C+A+90J2/πUL……………………………………计算第一缓和段切线方位角
L50:J=Q-D-L-M-G=>S=C+A+B-90J2/πUM………………………计算圆曲线段切线方位角
L51:J=1=>S=C+A+180(Q-D-L/2)/πU……………………………计算第二缓和段切线方位角
L52:LbI7…………………………………………………………………………将该处行标设为6
L53:Fix3…………………………………………………………………………将小数点位设为3位
L54:X=X+WcosS…………………………………………………………计算边桩X坐标
L55:Y=Y+WsinS…………………………………………………………计算边桩X坐标
L56:Goto0………………………………………………………………………转到0。如果输出值只为坐标时此行才写
L57:D=#####.###:G=#####.###…………………………… 写入测站点X、Y坐标
L58:Prog AZIMUTH…………………………………………………………运行子程序“AZIMUTH”
L59:Goto0………………………………………………………………………转到0

各输入值含义: R-圆曲线半径;LSci1-第一缓和曲线长度;LSci2-第一缓和曲线长度;JD-交点桩号;JDX-交点X坐标;JDY-交点Y坐标;PJ-偏角,左负右正输入;AZIMUTH-起始边方位角;R.ANG-(右角)边桩与中桩切线方向右交角;STAT-起点桩号;PEG NO.-所求点桩号;WIDTH-(边中距)边桩到中桩的距离

二、方位角程序
简单说明:可单独运行(将Z值输入为一个非0的数),亦可配合平曲线程序计算放样点到测站点的平距和方位角。
程序名:【AZIMUTH】
L01:LbI0…………………………………………………………………………………………将该处行标设为0
L02:Norm………………………………………………………………………………………将计算器设为普通模式
L03:Deg…………………………………………………………………………………………将计算器设为角度状态。
L04:Z=0=>Goto2…………………………………………………………………………输入测站点(station)X、Y坐标,转到2
L05:V=0:W=0……………………………………………………………………………………将I,J的值分别置零
L06:LbI1…………………………………………………………………………………………将该处行标设为1
L07:{XYDG}……………………………………………………………………………………将XYDG设为变量
L08:D“XA”G“YA”X“XB”Y“YB” ………………………………………………………将XYDG注释
L09:LbI2……………………………………………………………………………………………将该处行标设为2
L10:POL(X-D,Y-G) …………………………………………………………………………计算两点间的角度和距离
L11:W<0=>W=360+W………………………………………………………………………如果算得的角度<0则将其加个360度
L12:Fix6……………………………………………………………………………………………将小数点设这6位
L13:W“AZIMUTH A-B”=Int W+.01Int(60FracW)+.006Frac(60FracW)…显示所算得的方位角值
L14:Fix3……………………………………………………………………………………………将小数点设这3位
L15:I“H.DIST=” ……………………………………………………………………………显示平距值
L16:Z≠0=>Goto0……………………………………………………………………………判断Z值并构成循环

怎么大写都变成了小写??还有幂的形式,例如: h=j-jxy5/40r2u2+jxy9/3456rxy4uxy4=>H=J-J^5/40R^2U^2+J^9/3456R^4U^4 下同,输入时请千万注意。

由本站论坛版主 半边书生 所写,作者不允许以任何形式发表于杂志与报纸上。
打印】【关闭
            本站的部分源程序是由站长由网络收集整理的,如有侵权,请告之,我会第一时间删除相关内容。
因时间原因,源码不可能每一个都进行了测试,所以不能保证源码全是正确的,提供源码只是提供一份思路,一个参考,方便写出专属于您自己的程序
Copyright 2003-2025 测量天地 (SurveySky.Com) All Rights Reserved.        
         备案许可证:新ICP备12001392号-1         | 关于我们 |  联系我们 |网站留言