1、附闭合导线计算 、 源程序 F1 A1 L1 Defm 4N-2 L2 N:A:B:Pol(C-A,D-B):WW=W+360T=W L3 K=0=M=T+180:E=C:F=D:GOTO 0: E:F:Pol(G-E,H-F):WW=W+360M=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”=40NR=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
2、 L9 Z3N-2+R=W:I=I+W:N=R= GOTO 4: R=R+1:GOTO 3 L10 Lbl 4:P=U+C-EQ=I+D-F L11 G”1:M”=L/Pol(P,Q) R=2 L12 Lbl 5:X”XI”=C+Z2N-1+R-PZN+R/LY”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）、本程序可计算附和导线和闭合导线的坐标，计算的坐标系经过角度闭合差及坐标增量闭合差分配后的结果，能显示角度闭合差、增量闭合差及导线全长的相对精度； （2）、输入的观测角为
3、导线的左角。 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?
4、转换符，当 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”:
7、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
8、 可以算出，仪器架设在 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
9、-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 时为凸曲线，zJ=J+360J”A1-2=” 操作过程： ZBFSEXE输入 X1 值EXE输
10、入 Y1 值EXE输入 X2 值EXEY2EXE EXEEXE S1-2：计算得出的距离； A1-2：计算得出的角度。 （按 shift 转换为 60 进制的角度） 注：此程序可循环计算。坐标反算 ( 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 进制的
11、角度） 注：此程序可循环计算。缓和曲线曲线要素 ( for CASIO fx-4800P ) 程序步骤： B=01718.87L/R X=L-L3/40/ R2 Y=L2/6/ R P=Y-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 角（例 1253
12、123.25）EXE得 T 值EXE得 E 值EXE 得 Z 值 EXE得 J 值 EXE 得 D 值 注：此程序可循环计算。铁路曲线坐标计算程序 ( 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:
13、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=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
14、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；缓直 点坐标 HZ
15、X,HZY；直缓点里程 ZHDK；经过直缓点和交点的直线的方位角 ZHJD； 曲线偏 角（曲线左偏为负，右偏为正）A；缓和曲线长 L0；圆曲线长 L1;圆曲线半径 R； 计算 点里程 DK；计算曲线的左边点或者右边点）Z/!/Y(1/2/3)，曲线的右边输入 1，左 边输 入 3，线路上输入 2；左边的点与当前里程点的连线的距离和当前里程点切线的逆 时针 方向的角度 A（如果上一项不选择 2 的话） 。 3、计算结果显示的数据是：极角、极距和坐标。 4、该程序适用于极坐标法测量曲线时的数据（极角和极距）计算，以及曲线坐标 计算。 程序设计：杨序邦 2002 年 7 月 重庆直线段边坡超欠挖检查
16、及开口线放样程序 程序名：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(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 为边坡
17、底线的起点，显示为: 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； 中间计算测量点至边坡底线的垂足点，保存变量为 X、Y、Z。 O 为边坡的高程超欠挖值，显示为 Hcw 或 Hqw，分别表示超挖或欠挖，即测点高程设计高程减实测高程； N 为边坡的距离超欠挖值，显示为 Scw 或 Sqw，分别表示超挖或欠挖，按高差计算的设计距离减实测点到垂足点距离。圆弧段边坡超欠挖检查及开口线放样程序
18、 程序名：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”=S S0=M“Sqw”=SVK“Hcw”=V V0=K“Hqw”=V“END” 说明： A、 B 为圆心坐标，显示为 YX？YY？，依次输入圆心点的 X、Y； X、Y、Z 为测点坐标，显示为 XP？YP？ZP ？依次输入测点的 X、Y，Z； R 为到边坡底线的距离（半径） ； H 为边坡底线的高程（马道或平台的高程） ； I 不用管它，直接按“EXE”，显示下一个输入；
19、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 - T : B = Q + T ? Lbl 0 : I : I “ S Q “ I B G = H + E ( I - Q ) ? = S = I - Q : I D =
20、 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 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 交点里
21、程 QZH 交点高程 CP 前坡 HP 后坡 SQ 所求里程 曲线任意里程中边桩坐标正反算(CASIO fx-4800P 计算器) 程序 一、程序功能 本程序由一个主程序(TYQXJS)和两个子程?正算子程序(SUB1)、反算子程序( SUB2)序构成，可以根据曲线段?直线、圆曲线、缓和曲线（完整或非完整型）的线 元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲 率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。另 外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲 线任意里程中边桩坐标进行正反算。本程序也可以在 CASI
22、O fx-4500P 计算器及 CASIO fx-4850P 计算器上运行。 特别申明：(1). 适用于弧长小于 2 倍半径的各种线元坐标正反算，精度优 于 1mm； (2). 在引用该核心计算部分时，请注明来源。 二、源程序 1.主程序(TYQXJS) “1.SZ = XY“：“2.XY = SZ“：N ：U“X0“ ：V“Y0“ ：O“S0“：G“F0“：H“LS“ ：P“R0“：R“ RN“：Q ：C=1P ：D=(P-R)(2HPR)：E=180：N=1=Goto 1：Goto 2 Lbl 1： SZ：SZ：W=Abs(S-O)：Prog “SUB1“：X“XS“=X? Y“YS“=Y
23、? F“FS“=F-90? Goto 1 Lbl 2： XY：XY：I=X：J=Y：Prog “SUB2“：S“S“=O+W? Z“Z“=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(G+QEMW(C+MWD)：Y=V+W(Asin(G+QEKW(C+KWD)+Bsin(G+ QELW(C+LWD)+Bsin
24、(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 三、使用说明 1、规定 (1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时， Q=-1；当线元往右偏时
25、，Q=1；当线元为直线时，Q=0。 (2) 当所求点位于中线时，Z=0；当位于中线左铡时，Z 取负值；当位于中线中线右 侧时，Z 取正值。 (3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以 10 的 45 次代替。 (4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆 弧的半径。 (5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以 10 的 45 次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半 径为无穷大，以 10 的 45 次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。 (6) 当线元为非完整缓和曲线时，起
26、点与直线相接时，曲率半径等于设计规定的 值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等 于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。 2、输入与显示说明 输入部分： 1. SZ = XY 2. XY = SZ N ? 选择计算方式，输入 1 表示进行由里程、边距计算坐标 ；输入 2 表示由坐标反算 里程和边距。 X0 ？线元起点的 X 坐标 Y0 ？线元起点的 Y 坐标 S0 ？线元起点里程 F0 ？线元起点切线方位角 LS ？线元长度 R0 ？线元起点曲率半径 RN ？线元止点曲率半径 Q ？ 线 元左右偏标志(左偏 Q=-1，右偏 Q=1，直线段
27、Q=0） S ？ 正算时所求点的里程 Z ？正算时所求点距中线的边距(左侧取负，值右侧取正值，在中线上取零) X ？反算时所求点的 X 坐标 Y ？反算时所求点的 Y 坐标 显示部分： XS= 正算时，计算得出的所求点的 X 坐标 YS= 正算时，计算得出的所求点的 Y 坐标 FS= 正算时，所求点对应的中线点的切线方位角 S= 反算时，计算得出的所求点的里程 Z= 反算时，计算得出的所求点的边距 四、算例 某匝道的由五段线元（直线+完整缓和曲线+圆曲线+ 非完整缓和曲线+ 直线）组 成，各段线元的要素（起点里程 S0、起点坐标 X0 Y0、起点切线方位角 F0、线元长度 LS、起点曲率半径
28、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.75 221.75 -1 919.527 19736.072 28701.893 91 17 30.63 80.285 221.75 9579.
29、228 -1 999.812 19744.038 28781.659 80 40 50.00 100.000 1E45 1E45 0 1、正算 （注意：略去计算方式及线元要素输入，请自行根据所求点所在的线元输入线元 要素） S=700 Z=-5 计算得 XS=19831.41785 YS=28509.72590 S=700 Z=0 计算得 XS=19827.33592 YS=28506.83837 S=700 Z= 5 计算得 XS=19823.25398 YS=28503.95084 S=780 Z=-5 计算得 XS=19785.25749 YS=28575.02270 S=780 Z=0
30、 计算得 XS=19781.15561 YS=28572.16358 S=780 Z= 5 计算得 XS=19777.05373 YS=28569.30446 S=870 Z=-5 计算得 XS=19747.53609 YS=28654.13091 S=870 Z=0 计算得 XS=19742.68648 YS=28652.91379 S=870 Z= 5 计算得 XS=19737.83688 YS=28651.69668 S=940 Z=-5.123 计算得 XS=19741. 59118 YS=28722.05802 S=940 Z=0 计算得 XS=19736.47687 YS=2872
31、2.35642 S=940 Z= 3.009 计算得 XS=19733.47298 YS=28722.53168 2、 反算 X=19831.418 Y=28509.726 计算得 S=699.9999974 Z= -5 .00018164 X=19827.336 Y=28506.838 计算得 S=699.9996493 Z= 0.000145136 X=19823.25398 Y=28503.95084 计算得 S=699.9999985 Z= 5.000003137 X=19785.25749 Y=28575.02270 计算得 S=780.0000035 Z= -5 .00000166
32、3 X=19781.15561 Y=28572.16358 计算得 S=780.0000025 Z=- 0.000002979 X=19777.05373 Y=28569.30446 计算得 S=780.0000016 Z= 4.99999578 X=19747.536 Y=28654.131 计算得 S=870.0001137 Z= -4.99941049 X=19742.686 Y=28652.914 计算得 S=870.0003175 Z=- 0.00041814 X=19737.837 Y=28651.697 计算得 S=870.0002748 Z= 4.999808656 X=197
33、41.5912 Y=28722.0580 计算得 S=939.9999786 Z= -5.123024937 X=19736.4769 Y=28722.3564 计算得 S=939.9999862 Z=- 0.000027710 X=19733.4730 Y=28722.5317 计算得 S=940.0000238 Z= 3.00898694任意里程中边桩坐标正反算（CASIO fx-4800P 计算器）程序 一、源程序(SZ-XY-SZ) “1.SZ = XY“：“2. XY = SZ“：N ：U“X0“：V“Y0“：O“S0“：G“F0“：Q：Q=0=R=1E-45： RJ=0：M=Abs
34、Q：K=1-M：N=1=Goto 1：Goto 2 Lbl 1：SZ：SZ ：L=S-O：B=90QLR：F=G+B：E=F+B+90 ：C=LK+2MRsinQB：X“XS“ = U+Rec(C,F)+ZcosE? Y“YS“=V+J+ZsinE? Goto 1 Lbl 2： W=G+90Q：E=U+Rec(R,W)：F=V+J Lbl 3： XY：XY ：D=Pol(X-E,Y-F):T=J：JT=J+360B=T-W-180M：P=Abs B：S“S“=O+ DKcosP+ MRP180? Z“Z“=DKsinB+MQ(R-D)? Goto 3 二、使用说明 1.由里程、中边桩计算坐标
35、启动程序后按以下提示输入数据： N ？1 （N=1 ，由里程中边桩计算坐标； N1，由坐标反算里程及中边桩) X0 ？142.353 (直线或圆曲线起点的 X 坐标） Y0 ？368.166 (直线或圆曲线起点的 Y 坐标） S0 ？155.252 (直线或圆曲线起点的里程） F0 ？605403(直线或圆曲线起点的切线方位角） Q ？-1 (当线元为直线时 Q=0；当线元为圆曲线且左偏时 Q=-1；当线元为圆曲线且 右偏时 Q=1) R ？ 100 (圆曲线的半径) S ？ 194.016 (中线里程) Z ？ -5.251 (中线左侧 5.251m 的边桩；在中线上 Z=0，在左侧 0)
36、显示计算结果 XS=172.7467767 YS=395.2316807 2.由坐标反算里程、中边桩 启动程序后按以下提示输入数据： N ？2 （N=1 ，由里程中边桩计算坐标； N1，由坐标反算里程及中边桩) X0 ？142.353 (直线或圆曲线起点的 X 坐标） Y0 ？368.166 (直线或圆曲线起点的 Y 坐标） S0 ？155.252 (直线或圆曲线起点的里程） F0 ？605403(直线或圆曲线起点的切线方位角） Q ？-1 (当线元为直线时 Q=0；当线元为圆曲线且左偏时 Q=-1；当线元为圆曲线且 右偏时 Q=1) R ？ 100 (圆曲线的半径) X ？ 172.7467
37、767 (待求里程及中边桩的点的 X 坐标) Y ？ 395.2316807 (待求里程及中边桩的点的 Y 坐标) 显示计算结果 S=194.016 Z=-5.251 三、注意事项 1.本程序主要是针对市政道路工程中的直线和圆曲线而编制的，未考虑缓和曲线等 曲线的计算； 2.程序中“R=1E-45”的意思为 10 的-45 次方，输入程序时切莫输错。 3.程序中设置的计算顺序是由小里程往大里程方向（即里程增大的方向） ，曲线左偏、右偏也是以里程增大的方向来区分的。公路三维坐标控制系统高程控制部分(for 4850) 主程序 DH SZ:SZ SS1N=1:Prog“HDAT”:GOTO 1 S
38、S2N=2:Prog“HDAT”:GOTO 1 SS3N=3:Prog“HDAT”:GOTO 1 Lbl 1:S I Q= F+G(S- E S I Q= F+G(S- E)+(-1)J(S-I)2/(2 M R= Q+11.75 P-0.03:”DH=”: Q= Q ? Z=0 Q= Q Z0Q= R H:H”H=”: T= Q-H “H=”: T= T ? ABKL: T=Abs T: U= T / 8 : T 8W= TA+14+L T 8 W=8A+B(T-8)+14+2Int U+L T 16W=8A+8B+K(T-16)+14+2 Int U+L “DW”: W= W ? “W”:
39、 D= W-Abs Z? 竖曲线数据库 HDAT N=1E= xxxx.xxx: F= xxxx.xxx: G= xxxx.xxx: J= xxxx.xxx: M= xxxx.xxx: O= xxxx.xxxN=2E= xxxx.xxx: F= xxxx.xxx: G= xxxx.xxx: J= xxxx.xxx: M= xxxx.xxx: O= xxxx.xxxN=3E= xxxx.xxx: F= xxxx.xxx: G= xxxx.xxx: J= xxxx.xxx: M= xxxx.xxx: O= xxxx.xxx I= E- O 说明： (1)变量赋值说明 E 变坡点桩号 F 变坡点高程
40、 G 变坡点前段坡比 I 该段竖曲线起点桩号即高差改正值为零处桩号 J 竖曲线系数 1 为凸 2 为凹 M 竖曲线半径 O 竖曲线切线长 P 路面横坡（已在坐标反算程序中赋值） Q 路线设计高程，当 Z=0 时为中桩设计高程（或中分带边缘设计高程） ，当 Z0 时为路基边缘设计高程 R 路基边缘设计高程 T：“H” 高差 V 边坡分级，第 n 级坡 W 开口、坡脚设计宽度（根据高差及坡比求得）可计算任意高度处的设计宽度 D 宽差（当其为负? 向中桩方向移相应宽度，当其为正则向外移相应宽度） S、Z 所测点对应桩号及距中桩处的距离 A、B、K 横断面一、二、三级坡比 L 附加宽度（以便根据横断面
41、设计实际情况调整，象坠落台、护坡道、路堑处边沟设置宽度的不固定值设置及填方段为使路基满足设计的压实要求超填的宽度） H 实测高程 (2)编制思路 该程序编者根据工地的实际情况编制在此仅供一个编制思路，用户需根据实际情况进行修改。 注意：程序中有下划线处数据说明（单位：cm） 11.75 为中央分隔带边缘距硬路肩宽度，该段横坡为一定值（超高段需加入超高计算子程序，另外在设计宽度计算时加入加宽子程序即可） 003 为土路肩内外侧高差，因土路肩与行车道及硬路肩横坡不同，故分开计算。 （用户需根据实际情况修改） 2 为平台设置宽度，第 n 级坡的平台设置总宽度即为 2n. （用户需根据实际情况修改）
42、8 每级边坡设置高度. 14 路面设计宽度（路线中线处距土路肩边缘设计宽度） 该程序为公路三维控制系统中的一部分，即高程控制部分。该程序可与 hanniu_1973 老师的曲线任意里程中边桩坐标正反算(CASIO fx-4850P)程序（第三次修改版） 及大部分坐标正反算程序联合使用（修改变量名后） 。求两直线交点(CASIO fx?4800P 计算器)程序 ;直线(x1, y1)-(x2, y2) 与直线 (x3, y3)-(x3, y4)的交点 Lbl 1： ABCDEFGH：A“X1“：B“Y1“：C“X2“：D“Y2“ ：E“X3“：F“Y3“：G“X4“： H“Y4“： N=C-A：
43、O=D-B ：P=G-E：Q=H-F ：K=QN-OP： U=NB-OA：V=PF-QE： X“XP“=(PU-NV)K? Y“YP=“(QU-OV)K? Goto 1CASIO fx-4800P 编程计算器测量程序集 一SHI JU CE TU （视距测图） 1 用途：该程序用于“经纬仪视距法碎部测量 ”时，计算测站点至碎部点间的平距及碎部点的高程。 2 程序清单： C“H0”：I：A“UP”：B“DOWN”：V“MIDDLE” ：L：E=90-L D=100Abs（A-B） （cosE ）2 H=（C+I）+DtanE-V 3 程序说明： H0测站点高程；I测站仪器高； UP塔尺上丝读数；
44、DOWN塔尺下丝读数；MIDDLE塔尺中丝读数；L经纬仪竖盘读数盘左读数；D测站点至碎部点间平距；H碎部点的高程。 4该程序若在 fx-3950 中的程序清单为： ？M：？A：？B：？C：？D：100（A-B ） （ cos(90-D)）2X M+Xtan(90-D)-CY 说明：M将测站高程 H0+仪器高 I；A 上丝读数；B下丝读数；C中丝读数；D竖盘盘左读数 L。X 测站点至碎部点的平距；Y碎部点的高程。 二FAN SUAN （反算） 1用途：该程序是“坐标反算公式”程序，用于根据两点 1、2 的平面坐标（X1，Y1 ） 、（X2， Y2） ，计算两点间的距离 D12 及坐标方位角 A1
45、2。 2 程序清单： G“X1”：B“Y1”：C“X2”：H“Y2”：E=C-G：F=H-B Pol（ E，F ）：I“D12” J“A12” 3 程序说明： X1，Y1点 1 的平面坐标；X2，Y2点 2 的平面坐标； D12点 1 至点 2 间平距；A12点 1 至点 2 坐标方位角。 由于编程时用到了极坐标函数 POL（） ，故在程序运行中，会出现 I=？XXX（数字）及J=？XXX（数字）这是函数 POL 计算出的两个量，不能改变，直接回车默认即可。 三ZHENG SUAN （正算） 1用途：该程序是“坐标正算公式”程序，用于根据点 1 的平面坐标（X1，Y1） ，及点 1 至点 2
46、的平距 D12、坐标方位角 A12，来计算点 2 的平面坐标（X2，Y2 ） 。 2 程序清单： B“X1”：C“Y1”：D“D12”：A“A12” X“X2”=B+DcosA Y“Y2”=C+DsinA 3 程序说明： X1，Y1点 1 的平面坐标；X2，Y2点 2 的平面坐标； D12点 1 至点 2 间平距；A12点 1 至点 2 坐标方位角；X2 ，Y2点 2 的平面坐标。另外在输入 A12 时，可直接输入“度分秒”。 四YUAN ZHU （圆主） 1用途：该程序是“单圆曲线主点测设元素及里程计算”程序。 2程序清单： A“TURNING ANGLE”：R：B“JD” T=Rtan(A
47、2) L=RA180 E=Rcos(A2)-R D=2T-L C“ZY”=B-T F“YZ”=C+L G“QZ”=F-L2 3程序说明： TURNING ANGLE交点 JD 处的转角；R圆曲线半径；JD交点。 T切线长；L圆曲线长；E外矢距；D切曲差；ZYZY 点里程；YZYZ 点里程；QZQZ 点里程。CASIO FX-4800P 计算器在导线测量中的应用 一、 概述 导线测量是建立国家平面控制网的方法之一，也可用于工程建 设、城市建设的平面控制测量中。在水电、矿山、道路建设中也经常遇到。随着可编程计算器（FX-4800P）在导线测量中的应用。我们根据导线的计算方法，利用 FX-4800P 可编程功能适时计算出导线点的坐标。在速度、精度、准确性方面都能满足要求。 二、 附和导线的计算 附和导线就是起始于一个已知控制点，而终止于另一个已知控 制点的导线。 计算公式： 1、 角度闭合差及其分配 f=-(BM-NA)-N*180 V=- f/N 2、 坐标闭合差及其分配 fx=XB-XB fy=YB-YB Vxi=-(Si/S)*fx Vyi=-(Si/S)*fy 3、 计算待定点坐标 Xi=Xi+Vxi Yi=Yi+Vyi 4、 FX-4800P 计算程序 其中 N 为测站数 fix 0:N:fix 3:Defm 6
