附和导线和闭合导线的坐标的计算程序.doc

附和导线和闭合导线的坐标计算 1、 源程序 F1 A1 L1 Defm 4N-2 L2 N:A:B:Pol(C-A,D-B):WW=W+360?T=W L3 K=0=M=T+180:E=C:F=D:GOTO 0: E:F:Pol(G-E,H-F):WW=W+360?M=W L4 Lbl 0:L=0:U=0:I=0:R=2:Z1=T L5 Lbl 1:J:ZR+360: ?R=N+1=GOTO 2: R=R+1:GOTO 1 L6 Lbl 2:P”JB”=(ZN+1-M) ?Q”JL”=40N?R=2 L7 Lb1 3:S:ZN+R=S:L=L+S? L8 Z2N-1+R=Rec(S,(ZR-P(R-1)/N):U=U+V L9 Z3N-2+R=W:I=I+W:N=R=GOT 4: R=R+1:GOTO 3 L10 Lbl 4:P=U+C-E?Q=I+D-F? L11 G”1:M”=L/Pol(P,Q) ?R=2 L12 Lbl 5:X”XI”=C+Z2N-1+R-PZN+R/L?Y”YI”=D+Z3N-2+R-QZN+R/L? L13 R=N=GOTO 6: R=R+1:C=X:D=Y:GOTO 5 L14 Lbl 6:”END” 2、 说明 （1）、本程序可计算附和导线和闭合导线的坐标，计算的坐标系经过角度闭合差及坐标增量闭合差分配后的结果，能显示角度闭合差、增量闭合差及导线全长的相对精度； （2）、输入的观测角为导线的左角。 3、程序代号注释 N?导线观测角的折角数； A、B?导线起始点所后视的已知点的坐标x,y； C、D?导线起始点（即设站点）的坐标x,y； E、F?导线终点（已知点）的坐标x,y； G、H?在导线终点设站观测前视已知点的坐标x,y； T?起始站后视至起始点的方位角； M?终点站至前视已知点的方位角； J?观测的左角值； JB?角度闭合差； JL?允许的角度闭合差，程序中是以40n计算的，如和要求的不一致，可改一下L6语句中的有关部分。 S?所测导线的边长； L?边长的累计数； U?x的累计数； I?y的累计数； P?x坐标的闭合差； Q?y坐标的闭合差； K?转换符，当K=0时为计算闭合导线，当K0（任意数）时为计算附和导线。 铁路曲线坐标计算程序 ( for CASIO fx-4500P ) 一、主程序：M（M是文件名，下同) L1 A”CX”B”CY”C”HX”D”HY”E”ZHX”F”ZHY”Q”HZX”L”HZY”G”ZHDK”H”ZHJD” I”A”J”L0”K”L1”R L2 Pol(C-A,D-B):M=W L3 Lbl 2:W=90:NOWZ:Z”DK”N”Z/!/Y(1/2/3)”:N=2=Goto 1O”M”W”A” L4 Lbl 1:P=Z-G:S=E:T=F:U=H:V=H:Fixm L5 P0=Prog 3:=PJ=Prog 1:=PJ+K=Prog 2:V=90J/R+180- 1R-1 (P-J)PJ+K=Prog 4:V=H+V L6 PJ+K=P=2J+K-P:S=Q:T=L:P0=Prog 1:X=-X:U=H+I:Prog 4: V=U-V: = U=H+I:V=U:P=-P:Prog 3 L7 N2=P=O:U=V-W:N=3=P=-PS=X:T=Y:Prog 3 L8 Pol(X-A,Y-B):W=W-M:WW=W+360 L9 W:”=”V:”S=”X:”X=”Y:”Y=”Goto 2 二、子程序：1 L1 V=90P2(RJ)-1 L2 X=P-Pxy5(40R2J2)-1 L3 Y=PV/540 三、子程序：2 L1 Y=180-1R-1(P-0.5J) L2 X=RsinY+J/2- Jxy3/240/R2 L3 Y=J/24/R+R-RcosY 四、子程序：3 L1 X=S+PcosU L2 Y=T+PsinU 五、子程序：4 L1 IY=-Y:V=-V L2 S=S+XcosU-YsinU L3 Y=T+XsinU+YcosU L4 X=S 说明：1、该程序适用于计算器 CASIO fx-4500PA。 2、程序符号定义说明：测站点坐标CX,CY；后视点坐标HX,HY；直缓点坐标 ZHX,ZHY；缓直 点坐标HZX,HZY；直缓点里程ZHDK；经过直缓点和交点的直线的方位角ZHJD； 曲线偏 角（曲线左偏为负，右偏为正）A；缓和曲线长L0；圆曲线长L1;圆曲线半径R； 计算 点里程DK；计算曲线的左边点或者右边点）Z/!/Y(1/2/3)，曲线的右边输入1，左 边输 入3，线路上输入2；左边的点与当前里程点的连线的距离和当前里程点切线的逆 时针 方向的角度A（如果上一项不选择2的话）。 3、计算结果显示的数据是：极角、极距和坐标。 4、该程序适用于极坐标法测量曲线时的数据（极角和极距）计算，以及曲线坐标 计算。 面积计算（多边形法） 1、源程序 F1 A2 L1 N:P=A:Q=B:S=0:I=2 L2 Lbl 0:C,D:F=(A+C)(B-D):S=S+F L3 A=C:B=D:I=I+1 L4 IN=GOTO 0 L5 F=(C+P)(D-Q):S=S+F:S”W”=S/2 3、 说明： （1）、本程序适用于所测断面为多边形闭合图形的面积计算。 （2）、折点坐标按顺时针方向输入，得出的面积为正，否则为负，绝对值是一样 的。 4、 程序代号注释 A、B计算面积起始点纵横坐标； C、D各转折点的纵横坐标； S代表计算过程中的有关面积； S“W”为图形最后需要的计算面积。 N多边形的折点个数。 体积计算 1、源程序 F1 A3 L1 J=0:H=0:WG L2 Lbl 0:NAB:NAB:P=A:Q=B:S=0:I=1 L3 Lbl 1:CD:S=S+(A+C)(B-D)/2:A=C:B=D:I=I+1 L4 IGOTO 1S=S+(C+P)(D+Q)/2 L5 J1=GOTO 2: L=G-H:V=(R+S+（R*S)*L/3W=W+V L6 Lbl 2:R=S:H=G:J=1:G:G:GOTO 0 2、说明 （1）、程序可自动计算每一断面面积，当进行到第二个断面时就会显示出12断 面间的体积，而后再进行第三断面面积计算，并累计出13断面之间的体 积。，直到最后得出需算断面间的总体积。 （2）、坐标输入时，应按顺时针方向逐个输入折点坐标，这样得出的面积为正 值，一个桩号折点输入完后，程序自动进入下一桩号的输入状态。 3、程序代号注释 G断面桩号； A、B断面起算折点的坐标； C、D断面上其他折点坐标； S断面面积； L断面间距； V本断面与前一断面之间计算出的体积。 NG桩号断面上的折点个数； W本断面之前所有体积之和。 在任意控制点上测定直线上的任意位置与高程数据计算 1、源程序 F1 A4 L1 A”X”:B”Y”:E”Z0”:C”X1”:D”Y1”:F”Z1”:G”H0”:I L2 Lbl 0:LJV:LJV L3 K= tanF:T=tan(E+L) L4 X”XP”=(B-D-AT+KC)/(K-T) Y”YP”=B+TX-AT L5 S=（X-A）2+（Y-B）2） L6 H=G+StanJ+I-VGOTO 0 2、说明 （1）、本程序功能：在已知断面上某一点的坐标及断面方向的方位角后，可直接 将仪器架设在邻近控制点上，为测设各个断面上的点提供数据，不需要一定要将 仪器架设在断面桩上测设断面，而所测的断面点均有坐标和高程。 （2）、基本原理：如图所示，A、B为已知控制点，其方位角为Z0，断1-断2为断 面线控制桩，其方位角Z1可以算出，仪器架设在A点，后视B点，转角L1、L2、 L3。，这时A1、A2，。的方位角也为已知，根据解析 几何，两方位直线相交，可解出1#，2#。的坐标及其与A的距离。 3、程序代号注释 A、B测站点坐标 C、D断面桩点坐标 L观测断面点的水平角 J观测断面点的竖直角 V觇标高 S测站至断面测点的距离 H断面测点的高程 I仪器高 Z0测站至后视点的方位角 Z1断面线的方位角 竖曲线计算 1、源程序 F1 A5 L1 BADTRZ L2 Lbl 0:C:C L3 Z1=H=A+(B-A)/T*(C-D)-(C-D)2/(2R) H=A+(B-A)/T*(C-D)+ (C-D) 2/(2R) L4 GOTO 0 2、说明 （1）、本程序的功能是根据道路施工纵断面图上的设计数据，算出竖曲线上各加 桩点的高程； （2）、本程序适用于由小桩号向大桩号端方向计算，在键入全部已知数据后，当 C出现时，只要键入该点的桩号，高程立即会显示出来。 3、程序代号注释 A起点（或终点）高程 C需计算点的桩号 T竖曲线切线长 D起点（或终点）高程 R竖曲线半径 HC桩号处的高程 B切线交点的高程 Z曲线凹凸判断符Z1时为凸曲线，zH=A+(B-A)/T*(C-D)-(C-D)2/(2R) H=A+(B-A)/T*(C-D)+ (C-D) 2/(2R) L4 GOTO 0 2、说明 （1）、本程序的功能是根据道路施工纵断面图上的设计数据，算出竖曲线上各加 桩点的高程； （2）、本程序适用于由小桩号向大桩号端方向计算，在键入全部已知数据后，当 C出现时，只要键入该点的桩号，高程立即会显示出来。 3、程序代号注释 A起点（或终点）高程 C需计算点的桩号 T竖曲线切线长 D起点（或终点）高程 R竖曲线半径 HC桩号处的高程 B切线交点的高程 Z曲线凹凸判断符Z1时为凸曲线，z0时为凹曲线。 两点测角前方交会坐标计算 1、源程序 F1 A6 L1 ABCDEF L2 X“XP”=（A/tanF+C/tanE-B+D）/(1/tanE+1/tanF) L3 Y“YP”=（B/tanF+D/tanE-C+A）/(1/tanE+1/tanF) 说明： E1#点的观测角 F2#点的观测角 1#、2#点的编号时应注意：面向交会点P的左侧定为1#点，右侧定为2#点。 线路中、边桩测量放样程序 (Ver 3.2) F1 XLCS （主程序，步数385） L1 Norm:Deg:U=OA0:Prog 1:Q=U:CX-JD:DY-JD:U=AA:R+L- :Prog 1: B=Abs U:R:SL0:EK-ZH L2 Fix 3:M=.5S-Sxy3/240R2:P=S2/24R:TT=(R+P)tan.5B+MFL= RB/180 +SF=F+E:Norm L3 KZJ:XY1:K1L:V=LK:U=0:Prog 4: G=X:H=Y:G X:HY KJ:KHS:XY1:K=1I:X=IX:Y=JY:V=JK:U=0:Prog 4 Fix 3 L4 VD0=Pol(X-G,Y-HFix 4:N=W:W0W=W+360Prog 2:W A0=W Norm:U=0:U:UAB:Prog 1:K=U L5 Lb1 0:U=0:U:UCS:XY1:U=0UZ:V=ZK: UBZ:R+L- :A0U =-UProg 4: XY:X:YProg 5: Prog 6:Goto 0 F2 1 （十进制六十进制子程序，步数19） L1 U=Int U+Frac U/.6+Frac 100U/90 F3 2 （六十进制十进制子程序，步数22） L1 60Frac W:W=Int W+.01Int Ans+.006Frac Ans F4 3 （缓和曲线上任意点坐标计算子程序，步数57） L1 Y=RS:X=V-Vxy5/40Y2:Y=Vxy3(1-Vxy4/56Y2)/6Y+URec(1,90V2/ Y:X=X-UW F5 4 （中桩、边桩坐标计算子程序，步数195） L1 VFX=T-UsinB+Rec(V-F+T,B:Y=W+UcosB:Goto 0 L2 VF-SV=F-V:Prog 3: Rec(1,B:U=X:X=T+TV-XV-YW:Y=TW- UW+YV: Goto 0 L3 VE+SY=P+R-Rec(R-U,180(V-E-.5S)/R:X=M+W:Goto 0 L4 VEV=V-E:Prog 3: X=V-E:Y=U L5 Lbl 0: A0Y=-YRec(1,Q:U=X:X=C-TV+XV-YW:Y=D- TW+UW+YV F6 5 （测设数据输出子程序，步数76） L1 Pol(X-G,Y-H:W=W-N:W0W=W+360W=K+W:W360 W= W-360 Fix 4:Prog 2:WAC=WFix 3:VDC=VNorm F7 6 （测设时移桩数据计算子程序，步数92） L1 Lbl 0:Norm:L=0:W=0 L2 Lbl 1:U=0:U:USC:U0W=W+U:L=L+1:Goto 1W0 U:UV: Prog 1:Fix 3:WSD=W/Lsin UWMOVE=V-WGoto 0 【使用说明】 1、 程序说明 (1) 本程序特点：可置镜任意点放样任意点（中桩点、边桩点、导线 点）；人性化设计，提示信息全面，并充分考虑工程术语和习惯；专门设计的角 度输入输出方式；绝对优秀的存储器分配方案。 (2) 本程序适用于放射螺旋线作为缓和曲线且前、后缓和曲线长度相同的 单曲线放样，可放样线路中桩和边桩，按里程和坐标放样均可。 (3) 角度输入（出）方式：度.分秒，如 56319 输入（出）为 56.0319。注意：分、秒必须为