卡西欧测量相关程序.docx

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-

2、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” 、 说明 （1）、本程序可计算附和导线和闭合导线的坐标，计算的坐标系经过角度闭合差及坐标增量闭合差分配后的结果，能显示角度闭合差、增量闭合差及导线全长

3、的相对精度； （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坐标的

4、闭合差； K?转换符，当K=0时为计算闭合导线，当K0（任意数）时为计算附和导线。 面积计算（多边形法） 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 、 说明： （1）、本程序适用于所测断面为多边形闭合图形的面积计算。 （2）、折点坐标按顺时针方向输入，得出的面积为正，否则为负，绝对值是一样 的。 、 程序代号注释 A、B计算面积起始点纵横坐标； C、D各转折点的纵横坐标； S代表计算过程

5、中的有关面积； 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断 面间的体积，而后再进行第三

6、断面面积计算，并累计出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

7、: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

8、、 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）、本程序的功能是

9、根据道路施工纵断面图上的设计数据，算出竖曲线上各加 桩点的高程； （2）、本程序适用于由小桩号向大桩号端方向计算，在键入全部已知数据后，当 C出现时，只要键入该点的桩号，高程立即会显示出来。 3、程序代号注释 A起点（或终点）高程 C需计算点的桩号 T竖曲线切线长 D起点（或终点）高程 R竖曲线半径 HC桩号处的高程 B切线交点的高程 Z曲线凹凸判断符Z1时为凸曲线，zJ=J+360J”A1-2=” 操作过程： ZBFSEXE输入X1值EXE输入Y1值EXE输入X2值EXEY2EXE EXEEXE S1-2：计算得出的距离； A1-2：计算得出的角度。（按shift 转换为60进制的角度） 注

10、：此程序可循环计算。坐标反算 ( for CASIO fx-4800P ) 程序步骤： C”X1”:D”Y1”:E”X2”:F”Y2”:Fixm:Pol(E-C,F-D:I”S1-2=” JO=J=J+360J”A1-2=” 操作过程： ZBFSEXE输入X1值EXE输入Y1值EXE输入X2值EXEY2EXE EXEEXE S1-2：计算得出的距离； A1-2：计算得出的角度。（按shift 转换为60进制的角度） 注：此程序可循环计算。缓和曲线曲线要素 ( for CASIO fx-4800P ) 程序步骤： B=01718.87L/R X=L-L3/40/ R2 Y=L2/6/ R P=Y

11、-R(1-cos B Q= X-Rsin B T=(R+P)tan(A/2)+Q E=(R+P)(cos(A/2)-1-R Z= R(A-2B)/180+2L J=2T-Z D=X-Y/tan B 操作过程： HQQXYSEXE输入L值（即缓和曲线总长）EXE输入R值（即圆曲线半径） EXE得角度EXE得Xh值EXE得Yh值EXE得P值EXE得Q值 EXE输入A角（例1253123.25）EXE得T值EXE得E值EXE得Z值 EXE得J值EXE得D值 注：此程序可循环计算。铁路曲线坐标计算程序 ( for CASIO fx-4500P ) 一、主程序：M（M是文件名，下同) L1 A”CX”B

12、”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=-

13、X:U=H+I:Prog 4: V=U-V: = U=H+I:V=U:P=-P:Prog 3 L7 N2=P=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+P

14、cosU 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，左 边

15、输 入3，线路上输入2；左边的点与当前里程点的连线的距离和当前里程点切线的逆 时针 方向的角度A（如果上一项不选择2的话）。 3、计算结果显示的数据是：极角、极距和坐标。 4、该程序适用于极坐标法测量曲线时的数据（极角和极距）计算，以及曲线坐标 计算。 程序设计：杨序邦 2002年7月 重庆直线段边坡超欠挖检查及开口线放样程序 程序名：BPFY 程序： A“QX”：B“QY”：C“QZ”：E“ZX”：F“ZY”：G“ZZ”：P“PB”： K“XP”：L“YP”：M“ZP”：Fixm：Pol(E-A，F-B：J0=J=J+360 D=Abs(K-A)*CosJ+(L-B)*SinJ)：S=Abs

16、(L-F)*CosJ-(K-E)*SinJ)： X=A+D*CosJ：Y=B+D*SinJ：Z=(G-C)/I*D+C：W=AbS(M-Z)：V=S/P-W： VO“Hcw”=V V0=O“Hqw”=V U=W*P-S：UN“Scw”=UU0=N“Sqw”=U“END” 说明： A、B、C为边坡底线的起点，显示为: QX?QY?QZ? 依次输入地线起点的X,Y,Z； E、F、G为边坡底线的终点，显示为: ZX?ZY?ZZ? 依次输入地线终点的X,Y,Z； P为坡比，无正负条件； K、L、M为测量点，显示为: XP?YP?PZ? 依次输入测量点的X,Y,Z； 中间计算测量点至边坡底线的垂足点，保

17、存变量为X、Y、Z。 O为边坡的高程超欠挖值，显示为Hcw或Hqw，分别表示超挖或欠挖，即测点高程设计高程减实测高程； N为边坡的距离超欠挖值，显示为Scw或Sqw，分别表示超挖或欠挖，按高差计算的设计距离减实测点到垂足点距离。圆弧段边坡超欠挖检查及开口线放样程序 程序名：YFFY 程序： A“YX”：B“YY”：X“XP”：Y“YP”：Z“ZP”： Fixm：Pol(A-X，B-Y)：T=Abs(I-R)：W=Abs(H-Z)： S=W*P-T：V=T/P-W：SM“Scw”=SS0=M“Sqw”=SVK“Hcw”=V V0=K“Hqw”=V“END” 说明： A、 B为圆心坐标，显示为YX

18、？YY？，依次输入圆心点的X、Y； X、Y、Z为测点坐标，显示为XP？YP？ZP？依次输入测点的X、Y，Z； R为到边坡底线的距离（半径）； H为边坡底线的高程（马道或平台的高程）； I不用管它，直接按“EXE”，显示下一个输入； P为坡比，无须输入正负号； M为计算的边坡距离超欠值，显示为Scw或Sqw，分别表示超挖或欠挖，距离超欠值为：设计距离减测量距离； K为计算的边坡高程超欠值，显示为Hcw或Hqw，分别表示超挖或欠挖，高程超欠值为：设计高程减测量高程。CASIO 计算器竖曲线计算程式 R B J T Q X Q Q Z L H Q Z H C C P E H P : A = Q -

19、T : B = Q + T ? Lbl 0 : I : I S Q I B ⇒ G = H + E ( I - Q ) ? = S = I - Q : I D = B - I : : P = S E : Goto 1 Lbl 1 : C 0 ⇒ Goto 2 ? = Goto 3 ? Lbl 2 : E 0 ⇒ Goto 4 ? = Goto 6 ? Lbl 3 : E 0 ⇒ Goto 7 ? = Goto 5 ? Lbl 4 : C E ⇒ Goto 6 ? = Goto 7 ? Lbl 5 : C E ⇒ Goto

20、7 ? = Goto 6 ? Lbl 6 : G G H = H + P - D2 2 R ? Goto 0 ? Lbl 7 : G G H = H + P + D2 2 R ? Goto 0 ? BJ半径 QX切线长 QZL交点里程 QZH交点高程 CP前坡 HP后坡 SQ所求里程 曲线任意里程中边桩坐标正反算(CASIO fx-4800P计算器)程序 一、程序功能 本程序由一个主程序(TYQXJS)和两个子程?正算子程序(SUB1)、反算子程序( SUB2)序构成，可以根据曲线段?直线、圆曲线、缓和曲线（完整或非完整型）的线 元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半

21、径、止点曲 率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。另 外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲 线任意里程中边桩坐标进行正反算。本程序也可以在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：R RN：Q：C=1P：D=(P-R)(2H

22、PR)：E=180：N=1=Goto 1：Goto 2 Lbl 1：SZ：SZ：W=Abs(S-O)：Prog SUB1：XXS=X? YYS=Y? FFS=F-90? Goto 1 Lbl 2：XY：XY：I=X：J=Y：Prog SUB2：SS=O+W? ZZ=Z? Goto 2 2. 正算子程序(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(

23、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+ZsinF 3. 反算子程序(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 0 Lbl 1：Z=0：Prog SUB1：Z=(J-Y)sinF 三、使用说

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

25、次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。 (6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的 值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等 于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。 2、输入与显示说明 输入部分： 1. SZ = XY 2. XY = SZ N ? 选择计算方式，输入1表示进行由里程、边距计算坐标 ；输入2表示由坐标反算 里程和边距。 X0 ？线元起点的X坐标 Y0 ？线元起点的Y坐标 S0 ？线元起点里程 F0 ？线元起点切线方位角 LS ？线元长度 R0 ？线元起点曲率半径 RN ？线元止点曲

26、率半径 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 Q 500.000 19942.837 28343.561 125 16 31.00 269.256 1E45 1E45 0 769.256 19787.340 28563.378 125 16 31.00 37.492 1E45 221.75 -1 806.748 19766.566 28594.574 120 25 54.07 112.779 221.7