卡西欧测量相关程序.docx
《卡西欧测量相关程序.docx》由会员分享，可在线阅读，更多相关《卡西欧测量相关程序.docx（53页珍藏版）》请在冰豆网上搜索。
卡西欧测量相关程序
附闭合导线计算
１、源程序
F1A1
L1Defm4N-2
L2N:
A:
B:
Pol（C-A,D-B）:
W<0=>W=W+360?
T=W
L3K=0=>M=T+180:
E=C:
F=D:
GOTO0:
≠>E:
F:
Pol（G-E,H-F）:
W<0=>W=W+360?
M=W
L4Lbl0:
L=0:
U=0:
I=0:
R=2:
Z[1]=T
L5Lbl1:
{J}:
Z[R]+360:
?
R=N+1=>GOTO2:
≠>R=R+1:
GOTO1
L6Lbl2:
P”JB”=（Z[N+1]-M）?
Q”JL”=40√N?
R=2
L7Lb13:
{S}:
Z[N+R]=S:
L=L+S?
L8Z[2N-1+R]=Rec（S,（Z[R]-P（R-1）/N））:
U=U+V
L9Z[3N-2+R]=W:
I=I+W:
N=R=>GOT4:
≠>R=R+1:
GOTO3
L10Lbl4:
P=U+C-E?
Q=I+D-F?
L11G”1:
M”=L/Pol（P,Q）?
R=2
L12Lbl5:
X”XI”=C+Z[2N-1+R]-PZ[N+R]/L?
Y”YI”=D+Z[3N-2+R]-QZ[N+R]/L?
L13R=N=>GOTO6:
≠>R=R+1:
C=X:
D=Y:
GOTO5
L14Lbl6:
”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、源程序
F1A2
L1N:
P=A:
Q=B:
S=0:
I=2
L2Lbl0:
{C,D}:
F=（A+C）（B-D）:
S=S+F
L3A=C:
B=D:
I=I+1
L4I≦N=>GOTO0⊿
L5F=（C+P）（D-Q）:
S=S+F:
S”W”=S/2◢
３、说明：
（1）、本程序适用于所测断面为多边形闭合图形的面积计算。
（2）、折点坐标按顺时针方向输入，得出的面积为正，否则为负，绝对值是一样
的。
４、程序代号注释
A、B—计算面积起始点纵横坐标；
C、D—各转折点的纵横坐标；
S—代表计算过程中的有关面积；
S“W”—为图形最后需要的计算面积。
N—多边形的折点个数。
体积计算
1、源程序
F1A3
L1J=0:
H=0:
WG
L2Lbl0:
{NAB}:
NAB:
P=A:
Q=B:
S=0:
I=1
L3Lbl1:
{CD}:
S=S+（A+C）（B-D）/2:
A=C:
B=D:
I=I+1
L4IGOTO1⊿S=S+（C+P）（D+Q）/2◢
L5J≠1=>GOTO2:
≠>L=G-H:
V=（R+S+√（R*S））*L/3◢⊿W=W+V◢
L6Lbl2:
R=S:
H=G:
J=1:
{G}:
G:
GOTO0
2、说明
（1）、程序可自动计算每一断面面积，当进行到第二个断面时就会显示出1~2断
面间的体积，而后再进行第三断面面积计算，并累计出1~3断面之间的体
积。
。
。
。
。
。
，直到最后得出需算断面间的总体积。
（2）、坐标输入时，应按顺时针方向逐个输入折点坐标，这样得出的面积为正
值，一个桩号折点输入完后，程序自动进入下一桩号的输入状态。
3、程序代号注释
G—断面桩号；
A、B—断面起算折点的坐标；
C、D—断面上其他折点坐标；
S—断面面积；
L—断面间距；
V—本断面与前一断面之间计算出的体积。
N—G桩号断面上的折点个数；
W—本断面之前所有体积之和。
在任意控制点上测定直线上的任意位置与高程数据计算
1、源程序
F1A4
L1A”X”:
B”Y”:
E”Z0”:
C”X1”:
D”Y1”:
F”Z1”:
G”H0”:
I
L2Lbl0:
{LJV}:
LJV
L3K=tanF:
T=tan（E+L）
L4X”XP”=（B-D-AT+KC）/（K-T）◢Y”YP”=B+TX-AT◢
L5S=√（（X-A）2+（Y-B）2）
L6H=G+StanJ+I-V◢GOTO0
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、源程序
F1A5
L1BADTRZ
L2Lbl0:
{C}:
C
L3Z≧1=>H=A+（B-A）/T*（C-D）-（C-D）2/（2R）◢≠>H=A+（B-A）/T*（C-D）+（C-D）
2/（2R）◢
L4GOTO0
2、说明
（1）、本程序的功能是根据道路施工纵断面图上的设计数据，算出竖曲线上各加
桩点的高程；
（2）、本程序适用于由小桩号向大桩号端方向计算，在键入全部已知数据后，当
C出现时，只要键入该点的桩号，高程立即会显示出来。
3、程序代号注释
A—起点（或终点）高程
C—需计算点的桩号
T—竖曲线切线长
D—起点（或终点）高程
R—竖曲线半径
H—C桩号处的高程
B—切线交点的高程
Z—曲线凹凸判断符Z≧1时为凸曲线，z<0时为凹曲线。
两点测角前方交会坐标计算
1、源程序
F1A6
L1ABCDEF
L2X“XP”=（A/tanF+C/tanE-B+D）/（1/tanE+1/tanF）◢
L3Y“YP”=（B/tanF+D/tanE-C+A）/（1/tanE+1/tanF）◢
说明：
E—1#点的观测角
F—2#点的观测角
1#、2#点的编号时应注意：
面向交会点P的左侧定为1#点，右侧定为2#点。
坐标反算
（forCASIOfx-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进制的角度）
注：
此程序可循环计算。
坐标反算
（forCASIOfx-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进制的角度）
注：
此程序可循环计算。
缓和曲线曲线要素
（forCASIOfx-4800P）
程序步骤：
B=0°1718.87′﹡L/R◢
X=L-L＾3/40/R2◢
Y=L2/6/R◢
P=Y-R（1-cosB◢
Q=X-RsinB◢
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/tanB
操作过程：
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值
注：
此程序可循环计算。
铁路曲线坐标计算程序
（forCASIOfx-4500P）
一、主程序：
M（M是文件名，下同）
L1A”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
L2Pol（C-A,D-B）:
M=W
L3Lbl2:
W=90:
{NOWZ}:
Z”DK”N”Z/!
/Y（1/2/3）”:
N=2=>Goto1⊿O”M”W”A°”
L4Lbl1:
P=Z-G:
S=E:
T=F:
U=H:
V=H:
Fixm
L5P≤0=>Prog3:
≠=>P≤J=>Prog1:
≠=>P≤J+K=>Prog2:
V=90J/π/R+180π-
1R-1
（P-J）⊿⊿P≤J+K=>Prog4:
V=H+V⊿⊿
L6P>J+K=>P=2J+K-P:
S=Q:
T=L:
P>0=>Prog1:
X=-X:
U=H+I:
Prog4:
V=U-V:
≠=>
U=H+I:
V=U:
P=-P:
Prog3⊿⊿
L7N≠2=>P=U=V-W:
N=3=>P=-P⊿S=X:
T=Y:
Prog3⊿
L8Pol（X-A,Y-B）:
W=W-M:
W<0=>W=W+360⊿
L9W:
”°°°=”◢V:
”S=”◢X:
”X=”◢Y:
”Y=”◢Goto2
二、子程序：
1
L1V=90P2（πRJ）-1
L2X=P-Pxy5（40R2J2）-1
L3Y=PVπ/540
三、子程序：
2
L1Y=180π-1R-1（P-0.5J）
L2X=RsinY+J/2-Jxy3/240/R2
L3Y=J/24/R+R-RcosY
四、子程序：
3
L1X=S+PcosU
L2Y=T+PsinU
五、子程序：
4
L1I<0=>Y=-Y:
V=-V⊿
L2S=S+XcosU-YsinU
L3Y=T+XsinU+YcosU
L4X=S
说明：
1、该程序适用于计算器CASIOfx-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"BJ"T"QX"Q"QZL"
H"QZH"C"CP"E"HP"
:
A=Q-T:
B=Q+T?
Lbl0:
｛I｝:
I"SQ"=H-C（Q-I）?
≠=>I>B⇒G=H+E（I-Q）?
≠=>S=I-Q:
I:
P=SC:
Goto1:
≠=>D=B-I:
:
P=SE:
Goto1
Lbl1:
C>0⇒Goto2?
≠=>Goto3?
Lbl2:
E>0⇒Goto4?
≠=>Goto6?
Lbl3:
E>0⇒Goto7?
≠=>Goto5?
Lbl4:
C>E⇒Goto6?
≠=>Goto7?
Lbl5:
C>E⇒Goto7?
≠=>Goto6?
Lbl6:
G"GH"=H+P-D2÷2
÷R?
Goto0?
Lbl7:
G"GH"=H+P+D2÷2
÷R?
Goto0?
BJ半径QX切线长QZL交点里程
QZH交点高程CP前坡HP后坡
SQ所求里程
曲线任意里程中边桩坐标正反算（CASIOfx-4800P计算器）程序
一、程序功能
本程序由一个主程序（TYQXJS）和两个子程?
?
正算子程序（SUB1）、反算子程序（
SUB2）序构成，可以根据曲线段?
?
直线、圆曲线、缓和曲线（完整或非完整型）的线
元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲
率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。
另
外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲
线任意里程中边桩坐标进行正反算。
本程序也可以在CASIOfx-4500P计算器及CASIOfx-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=>Goto1：
≠>Goto2Δ←┘
Lbl1：
{SZ}：
SZ：
W=Abs（S-O）：
Prog"SUB1"：
X"XS"=X?
Y"YS"=Y?
F"FS"=F-90?
Goto1←┘
Lbl2：
{XY}：
XY：
I=X：
J=Y：
Prog"SUB2"：
S"S"=O+W?
Z"Z"=Z?
Goto2
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：
Lbl0：
Prog"SUB1"：
L=T+QEW（C+
WD）：
Z=（J-Y）cosL-（I-X）sinL：
AbsZ<1E-6=>Goto1：
≠>W=W+Z：
Goto0Δ←┘
Lbl1：
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、起点坐标X0Y0、起点切线方位角F0、线元长度
LS、起点曲率半径R0、止点曲率半径RN、线元左右偏标志Q）如下：
S0X0Y0F0LSR0RNQ
500.00019942.83728343.5611251631.00269.2561E451E450
769.25619787.34028563.3781251631.0037.4921E45221.75-1
806.74819766.56628594.5741202554.07112.779221.7
升级会员 