卡西欧测量相关程序.doc

附闭合导线计算
１、 源程序
F1 A1
L1 Defm 4N-2
L2 N:A:B:Pol(C-A,D-B):W<0=>W=W+360?T=W
L3 K=0=>M=T+180:E=C:F=D:GOTO 0: ≠>E:F:Pol(G-E,H-F):W<0=>W=W+360?M=W
L4 Lbl 0:L=0:U=0:I=0:R=2:Z[1]=T
L5 Lbl 1:{J}:Z[R]+360: ?R=N+1=>GOTO 2: ≠>R=R+1:GOTO 1
L6 Lbl 2:P”JB”=(Z[N+1]-M) ?Q”JL”=40√N?R=2
L7 Lb1 3:{S}:Z[N+R]=S:L=L+S?
L8 Z[2N-1+R]=Rec(S,(Z[R]-P(R-1)/N)):U=U+V
L9 Z[3N-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+Z[2N-1+R]-PZ[N+R]/L?Y”YI”=D+Z[3N-2+R]-QZ[N+R]/L?
L13 R=N=>GOTO 6: ≠>R=R+1:C=X:D=Y:GOTO 5
L14 Lbl 6:”END”

２、 说明
（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?允许的角度闭合差，程序中是以40√n计算的，如和要求的不一致，可改一下L6语句中的有关部分。
S?所测导线的边长；
L?边长的累计数；
U?△x的累计数；
I?△y的累计数；
P?x坐标的闭合差；
Q?y坐标的闭合差；
K?转换符，当K=0时为计算闭合导线，当K≠0（任意数）时为计算附和导线。

面积计算（多边形法）
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 I≦N=>GOTO 0⊿
L5 F=(C+P)(D-Q):S=S+F:S”W”=S/2◢

３、 说明：
（1）、本程序适用于所测断面为多边形闭合图形的面积计算。
（2）、折点坐标按顺时针方向输入，得出的面积为正，否则为负，绝对值是一样
的。
４、 程序代号注释
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 1⊿S=S+(C+P)(D+Q)/2◢
L5 J≠1=>GOTO 2: ≠>L=G-H:V=(R+S+√（R*S))*L/3◢⊿W=W+V◢
L6 Lbl 2:R=S:H=G:J=1:{G}:G:GOTO 0
2、说明
（1）、程序可自动计算每一断面面积，当进行到第二个断面时就会显示出1~2断
面间的体积，而后再进行第三断面面积计算，并累计出1~3断面之间的体
积。。。。。。，直到最后得出需算断面间的总体积。
（2）、坐标输入时，应按顺时针方向逐个输入折点坐标，这样得出的面积为正
值，一个桩号折点输入完后，程序自动进入下一桩号的输入状态。
3、程序代号注释
G—断面桩号；
A、B—断面起算折点的坐标；
C、D—断面上其他折点坐标；
S—断面面积；
L—断面间距；
V—本断面与前一断面之间计算出的体积。
N—G桩号断面上的折点个数；
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-V◢GOTO 0

2、说明
（1）、本程序功能：在已知断面上某一点的坐标及断面方向的方位角后，可直接
将仪器架设在邻近控制点上，为测设各个断面上的点提供数据，不需要一定要将
仪器架设在断面桩上测设断面，而所测的断面点均有坐标和高程。
（2）、基本原理：如图所示，A、B为已知控制点，其方位角为Z0，断1-断2为断
面线控制桩，其方位角Z1可以算出，仪器架设在A点，后视B点，转角L1、L2、
L3。。。。。。，这时A~1、A~2，。。。。。。的方位角也为已知，根据解析
几何，两方位直线相交，可解出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 Z≧1=>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—竖曲线半径
H—C桩号处的高程
B—切线交点的高程

Z—曲线凹凸判断符Z≧1时为凸曲线，z<0时为凹曲线。

两点测角前方交会坐标计算
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) ◢

说明：
E—1#点的观测角
F—2#点的观测角
1#、2#点的编号时应注意：面向交会点P的左侧定为1#点，右侧定为2#点。
坐标反算
( for CASIO fx-4800P )
程序步骤：

C”X1”:D”Y1”:E”X2”:F”Y2”:Fixm:Pol(E-C,F-D:I”S1-2=”◢

J≤O=>J=J+360△J”A1-2=”

操作过程：

ZBFS→EXE→输入X1值→EXE→输入Y1值→EXE→输入X2值→EXE→Y2→EXE→
EXE→EXE

S1-2：计算得出的距离；

A1-2：计算得出的角度。（按shift°′″ 转换为60进制的角度）

注：此程序可循环计算。

坐标反算
( for CASIO fx-4800P )
程序步骤：

C”X1”:D”Y1”:E”X2”:F”Y2”:Fixm:Pol(E-C,F-D:IS1-2=”◢

J≤O=>J=J+360△J”A1-2=”

操作过程：

ZBFS→EXE→输入X1值→EXE→输入Y1值→EXE→输入X2值→EXE→Y2→EXE→
EXE→EXE

S1-2：计算得出的距离；

A1-2：计算得出的角度。（按shift°′″ 转换为60进制的角度）

注：此程序可循环计算。
缓和曲线曲线要素
( for CASIO fx-4800P )

程序步骤：

B=0°1718.87′﹡L/R◢

X=L-L＾3/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

操作过程：

HQQXYS→EXE→输入L值（即缓和曲线总长）→EXE→输入R值（即圆曲线半径）
→EXE→得β角度→EXE→得Xh值→EXE→得Yh值→EXE→得P值→EXE→得Q值→
EXE→输入A角（例125°31′23.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 1⊿O”M”W”A°”

L4 Lbl 1:P=Z-G:S=E:T=F:U=H:V=H:Fixm

L5 P≤0=>Prog 3:≠=>P≤J=>Prog 1:≠=>P≤J+K=>Prog 2:V=90J/π/R+180π-
1R-1

(P-J)⊿⊿P≤J+K=>Prog 4:V=H+V⊿⊿

L6 P>J+K=>P=2J+K-P:S=Q:T=L:P>0=>Prog 1:X=-X:U=H+I:Prog 4: V=U-V:
≠=>

U=H+I:V=U:P=-P:Prog 3⊿⊿

L7 N≠2=>P=U=V-W:N=3=>P=-P⊿S=X:T=Y:Prog 3⊿

L8 Pol(X-A,Y-B):W=W-M:W<0=>W=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 I<0=>Y=-Y:V=-V⊿

L2 S=S+XcosUYsinU

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、该程序适用于极坐标法测量曲线时的数据（极角和极距）计算，以及曲线坐标
计算。

程序设计：杨序邦 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：J≤0=>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：
V<0=>O“Hcw”=V◢⊿
V≥0=>O“Hqw”=V◢⊿
U=W*P-S：U<0=>N“Scw”=U◢⊿U≥0=>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；
中间计算测量点至边坡底线的垂足点，保存变量为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：S<0=>M“Scw”=S◢⊿S≥0=>M“Sqw”=S◢⊿V<0=>K“Hcw”=V◢⊿
V≥0=>K“Hqw”=V◢⊿“END”
说明：
A、 B为圆心坐标，显示为YX？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 - T : B = Q + T ? 　 　 　
Lbl 0 : ｛ I ｝ : I " S Q " < A ⇒ G
= H - C ( Q - I ) ? 　 　 　 　 　 　
≠=> I > B ⇒ G = H + E ( I - Q ) ?
≠=> S = I - Q : I < Q ⇒ D = I - A
: P = S C : Goto 1 : ≠=> 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 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)序构成，可以根据曲线段??直线、圆曲线、缓和曲线（完整或非完整型）的线
元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲
率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。另
外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲
线任意里程中边桩坐标进行正反算。本程序也可以在CASIO 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=1÷P：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?
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(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：AbsZ<1E-6=>Goto1：≠>W=W+Z：Goto 0Δ←┘
Lbl 1：Z=0：Prog "SUB1"：Z=(J-Y)÷sinF

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

2、输入与显示说明
输入部分：
1. SZ => XY
2. XY = > SZ
N ? 选择计算方式，输入1表示进行由里程、边距计算坐标 ；输入2表示由坐标反算
里程和边距。
X0 ？线元起点的X坐标
Y0 ？线元起点的Y坐标
S0 ？线元起点里程
F0 ？线元起点切线方位角
LS ？线元长度
R0 ？线元起点曲率半径
RN ？线元止点曲率半径
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 28563378 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.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 计算得 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=28722.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 .000001663
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=19741.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：≠>
R△J=0：M=AbsQ：K=1-M：N=1=>Goto 1：≠>Goto 2△←┙
Lbl 1：{SZ}：SZ：L=S-O：B=90QL÷π÷R：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：J<0=>T=J+360△B=T-W-180M：P=Abs B：S"S"=O+
DKcosP+ MπRP÷180?
Z"Z"=DKsinB+MQ(R-D)?
Goto 3←┙
二、使用说明
1.由里程、中边桩计算坐标
启动程序后按以下提示输入数据：
N ？1 （N=1，由里程中边桩计算坐标；N≠1，由坐标反算里程及中边桩)
X0 ？142.353 (直线或圆曲线起点的X坐标）
Y0 ？368.166 (直线或圆曲线起点的Y坐标）
S0 ？155.252 (直线或圆曲线起点的里程）
F0 ？60°54′03″(直线或圆曲线起点的切线方位角）
Q ？-1 (当线元为直线时Q=0；当线元为圆曲线且左偏时Q=-1；当线元为圆曲线且
右偏时Q=1)
R ？ 100 (圆曲线的半径)
S ？ 194.016 (中线里程)
Z ？ -5.251 (中线左侧5.251m的边桩；在中线上Z=0，在左侧<0，在右侧Z>0)
显示计算结果
XS=172.7467767
YS=395.2316807
2.由坐标反算里程、中边桩
启动程序后按以下提示输入数据：
N ？2 （N=1，由里程中边桩计算坐标；N≠1，由坐标反算里程及中边桩)
X0 ？142.353 (直线或圆曲线起点的X坐标）
Y0 ？368.166 (直线或圆曲线起点的Y坐标）
S0 ？155.252 (直线或圆曲线起点的里程）
F0 ？60°54′03″(直线或圆曲线起点的切线方位角）
Q ？-1 (当线元为直线时Q=0；当线元为圆曲线且左偏时Q=-1；当线元为圆曲线且
右偏时Q=1)
R ？ 100 (圆曲线的半径)
X ？ 172.7467767 (待求里程及中边桩的点的X坐标)
Y ？ 395.2316807 (待求里程及中边桩的点的Y坐标)
显示计算结果
S=194.016
Z=-5.251
三、注意事项
1.本程序主要是针对市政道路工程中的直线和圆曲线而编制的，未考虑缓和曲线等
曲线的计算；
2.程序中“R=1E-45”的意思为10的-45次方，输入程序时切莫输错。
3.程序中设置的计算顺序是由小里程往大里程方向（即里程增大的方向），曲线左偏、右偏也是以里程增大的方向来区分的。
公路三维坐标控制系统高程控制部分(for 4850)

主程序DH

{SZ}:SZ
S≤S1＝＞N=1:Prog“HDAT”:GOTO 1△
S≤S2＝＞N=2:Prog“HDAT”:GOTO 1△
S≤S3＝＞N=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 △
Z≠0＝＞Q= R △
{H}:H”H=”: T= Q-H
“△H=”: T T ?
{ABKL}: T=Abs T: U= T / 8 : T ≤8＝＞W= TA+14+L△
T ＞8＝＞W=8A+B(T-8)+14+2Int U+L△
T ＞16＝＞W=8A+8B+K(T-16)+14+2 Int U+L△
“DW”: W= W ?
“△W”: D= W-Abs Z?

竖曲线数据库HDAT
N=1＝＞E= xxxx.xxx: F= xxxx.xxx: G= xxxx.xxx: J= xxxx.xxx: M= xxxx.xxx:
O= xxxx.xxx△N=2＝＞E= xxxx.xxx: F= xxxx.xxx: G= xxxxxxx: J= xxxx.xxx: M= xxxx.xxx: O= xxxx.xxx△N=3＝＞E= xxxx.xxx: F= xxxx.xxx: G= xxxx.xxx: J= xxxx.xxx: M= xxxx.xxx: O= xxxx.xxx△
……
I= E- O

说明：
(1)变量赋值说明

E变坡点桩号
F变坡点高程
G变坡点前段坡比
I该段竖曲线起点桩号即高差改正值为零处桩号
J竖曲线系数 1为凸2为凹
M竖曲线半径
O 竖曲线切线长
P路面横坡（已在坐标反算程序中赋值）
Q路线设计高程，当Z=0时为中桩设计高程（或中分带边缘设计高程），当Z≠0
时为路基边缘设计高程
R 路基边缘设计高程
T：“△H”高差
V边坡分级，第n级坡
W开口、坡脚设计宽度（根据高差及坡比求得）可计算任意高度处的设计宽度
D宽差（当其为负?向中桩方向移相应宽度，当其为正则向外移相应宽度）
S、Z所测点对应桩号及距中桩处的距离
A、B、K横断面一、二、三级坡比
L附加宽度（以便根据横断面设计实际情况调整，象坠落台、护坡道、路堑处边沟设置宽度的不固定值设置及填方段为使路基满足设计的压实要求超填的宽度）
H实测高程
(2)编制思路
该程序编者根据工地的实际情况编制在此仅供一个编制思路，用户需根据实际情况进行修改。
注意：程序中有下划线处数据说明（单位：cm）
11.75为中央分隔带边缘距硬路肩宽度，该段横坡为一定值（超高段需加入超高计算子程序，另外在设计宽度计算时加入加宽子程序即可）
0．03为土路肩内外侧高差，因土路肩与行车道及硬路肩横坡不同，故分开计算。（用户需根据实际情况修改）
2 为平台设置宽度，第n级坡的平台设置总宽度即为2n. （用户需根据实际情况修改）
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：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 1
CASIO 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——塔尺上丝读数；DOWN——塔尺下丝读数；MIDDLE——塔尺中丝读数；L——经纬仪竖盘读数盘左读数；D——测站点至碎部点间平距；