测量程序6.doc

利用Simpson公式和垂距趋近正反算曲线程序 作者：hanniu_1973 出处：海口 更新时间： 2005年10月13日 一、程序编制原理和功能： 1、原理：利用Simpson公式和点至主点法线垂距逐步迭代趋近单向收敛原理 2、功能：a、正算任意角度、偏距点坐标；b、反算任意角度的中桩里程和偏距。 二、程序主体： 1、主程序：ZFS T1.KDF=XY,2.XY=KD:AQDBZDR1/R0V1/RNCX0DY0EF0:T=1=Goto 1:≠Goto 2△ Lbl 1:{LIF}:LK=:IDP=:FANG=:Prog SIM2:XI=:X=X:Pause 0:YI=:Y=Y▲FW=:W→DMS▲Goto 1 Lbl 2:{XY}:XY:O=X:U=Y:L=A+Abs ((Y-D)cos(E+FANG=-180)-(X-C)sin(E+F-180)) Lbl 3:I=0:Prog SIM2:Z=(U-Y)cos(W+F-180)-(O-X)sin(W+F-180):AbsZ1E-4=Goto 4:≠L=L+Z:Goto 3△ Lbl 4:I=0:Prog SIM2:Z=(U-Y)/sin(W+F) Lbl 5:K=:L:Pause 0:DP=:Z▲Goto 2 2、子程序：SIM2 N=6：T=（V－R）÷（B－A）：H=（L－A）÷N：P=TNH+R：S=9θH÷π：W=E+（P+R）NS G=1： X =C+H÷6×（cosE+cosW+4×∑（cos（E+（（G+θ.5）HT+2R）（G+θ.5）S），G，θ，（N－1））+2∑（cos（E+（THG+2R）GS），G，1，（N－1）））+Icos(W+F)：Y =D+H÷6×（sinE+sinW+4×∑（sin（E+（（G+θ.5）HT+2R）（G+θ.5）S，G，θ，（N－1））+2×∑（sin（E+（THG+2R）GS，G，1,（N－1）））+Isin(W+F) 三、使用说明： QD—起点里程； ZD—终点里程； 1/Rθ、1/RN—起、终点曲率（左负右正） Xθ、Yθ—起点坐标； Fθ—起点方位角； N—可取4～12，视线型而定，亦可改为变量； ANG—线外左右测点连线与路线方向斜交右角； X、Y—待求点中桩坐标 XI、YI—计算点中边桩坐标。 FW—切线方位角 K= —正算时为计算点里程；反算时为测点对应中桩里程 DP—正算时为支距（左负右正）；反算时为测点对应于中线某一角度ANG的支距（若ANG=90则为垂支距，左负右正）。 曲线任意里程中边桩坐标正反算(CASIO fx-4850P计算器)程序(修改版） 作者：佚名 出处：SurM 更新时间： 2005年08月23日 曲线任意里程中边桩坐标正反算(CASIO fx-4850P计算器)程序（修改版） 一、程序功能及原理 1.功能说明：本程序由一个主程序(TYQXjs)和两个子程——正算子程序(SUB1)、反算子程序(SUB2)序构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。另外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲线任意里程中边桩坐标进行正反算。本程序也可以在CASIO fx-4500P计算器及 CASIO fx-4850P计算器上运行。 2．计算原理：利用Gauss-Legendre 5点通用公式计算线路中边桩坐标并计算放样数据。 二、源程序 1.主程序(TYQXjs) Defm 2:1.SZ = XY：2.XY = SZ：N：Prog”DAT1”:UX0：VY0：OS0：GF0：HLS：PR0：RRN：Q：C=1÷P：D=(P-R)÷(2HPR)：E=180÷π：N=1=Goto 1：≠Goto 2Δ← Lbl 1：{SZT}：SZ:T”ANG”：W=Abs(S-O)：Prog SUB1：XS=:X=X◢ YS=:Y=Y◢ Goto 3←┘ Lbl 2：{XY}：XY：I=X：J=Y：Prog SUB2：S=:S=O+W◢ Z=:Z=Z◢ Goto 2 Lbl 3:J=0:I=Pol(X-Z[1],Y-Z[2]):F=J:F0=F=F+360:≠F360=F=F-36ΔΔ ”DIST=”: I◢”FW=”:F→DMS◢ Goto 1 2. 正算子程序(SUB1) A=0.1184634425:B=0.2393143352:N=0.2844444444:K=0.0469100770:L=0.2307653449: