CASIOfx4500P计算器程序.docx
《CASIOfx4500P计算器程序.docx》由会员分享,可在线阅读,更多相关《CASIOfx4500P计算器程序.docx(26页珍藏版)》请在冰豆网上搜索。
CASIOfx4500P计算器程序
CASIOfx4500P计算器
实用程序
XCAJCJ编著
二○○一年八月
说 明
本程序稿由于时间仓促及编者水平有限,缺点、错误以及疏漏与不足之处在所难免,恳请大家批评指正。
程序步骤计算:
1个字符的文件名算3步(每增加一个字符多算一步),行号算1步,每个字符算1步。
CASIOfx-4500P共有1103步。
XCAJCJ
二○○一年八月
目 录
程序1——角度判断1
程序2——角度换算①1
程序3——角度换算②1
程序4——角度换算③1
程序5——角度换算④1
程序6——视距测量2
程序7——坐标反算2
程序8——坐标方位角计算2
程序9——坐标正算3
程序10——坐标方位角交会计算3
程序11——垂距计算4
程序12——视差法测距计算5
程序13——剖面交点计算6
程序14——前方交绘计算6
程序15——附合导线近似平差7
程序16——渠道中线里程计算9
程序17——多边形面积计算11
程序18——闭合导线近似平差13
程序19——直线交点坐标计算14
程序20——极坐标法放样计算15
附录1:
CASIOfx-4500P错误标示一览表16
附录2:
常用数学公式17
附录3:
常用测量公式17
附录4:
单一附合导线内业平差计算20
附录5:
勘察分院CAD制图《试行》规定21
程序1——角度判断
文件名
A
角 度 判 断
29步
L1
B≥360B=B-360△B<0B=B+360
说明:
①程序A为子程序,对≥360°或<0°的角值进行判断处理。
程序2——角度换算①
文件名
B
角 度 换 算 ①
23步
L1
B=IntB+FracB/.6+Frac100B/90
程序3——角度换算②
文件名
C
角 度 换 算 ②
22步
L1
B=IntB+.6FracB-.004Frac60B
程序4——角度换算③
文件名
D
角 度 换 算 ③
25步
L1
F=B:
B=I:
ProgB:
I=B:
B=F:
ProgB
程序5——角度换算④
文件名
E
角 度 换 算 ④
25步
L1
F=B:
B=I:
ProgB:
I=B:
B=F:
ProgB
说明:
①程序B、C、D、E均为子程序,简化测量计算中角度的输入及显示。
注:
子程序应先输入CAISOfx-4500P中,以备其它程序调用。
程序6——视距测量
文件名
1
视 距 测 量
41步
L1
Lbl0:
{SCB}:
HSC:
ProgB:
S(sinB)2▲Ans/tanB+H-C▲ProgC:
Goto0
说明:
①程序用以计算视距测量中测站至照准点的平距及高程;
②程序中:
H=测站高程+仪高(程序只要求输入该值一次);S=视距、C=切尺、B=竖盘读数;
③算例:
按Prog1调用程序显示H?
输入测站高程+仪高445.79+1.41,按EXE键、显示S?
输入视距39,按EXE键、显示C?
输入切尺2.9,按EXE键、显示B?
输入竖盘读数90.12,按EXE键、显示(计算结果)平距38.9995248,按EXE键、显示高程444.1638654,按EXE键、显示S?
、输入视距……
④竖盘读数90.12是将度分秒值(90°12′00″)以小数的形式输入;其它程序的角度值输入及显示均如此。
程序7——坐标反算
文件名
2
视 标 反 算
33步
L1
A"X0"C"Y0":
Pol(X-A,Y-C▲B=W:
PorgA:
ProgC
说明:
①程序用以计算两点间的距离及坐标方位角;
②程序中:
A=起点纵坐标X0,C=起点横坐标Y0,X=终点纵坐标,Y=终点横坐标;
③算例:
按Prog2调用程序显示X0?
输入起点纵坐标667.132,按EXE键、显示Y0?
输入起点横坐标993.505,按EXE键、显示X?
、输入终点纵坐标395.372,按EXE键、显示Y?
输入终点横坐标914.821,按EXE键、显示两点间距离282.9216666,按EXE键、显示起点至终点的方位角196.0851267。
程序8——坐标方位角计算
文件名
3
坐 标 方 位 角 计 算
26步
L1
ProgD:
B=B+I-180:
ProgA:
D=B:
ProgC
说明:
①程序用以推算导线边的坐标方位角(左角公式);
②程序中:
B=已知方位角、I=转折角(左角);
③语句B=D是为其它程序调用本程序而预设的,此处无具体用途;
④算例:
按Prog4调用程序显示B?
输入已知方位角120.10214按EXE键、显示I?
输入转折角160.01217,按EXE键、显示所求方位角100.11431;
⑤若转折角为右角,将语句B=B+I-180改为B=B-I+180即可。
程序9——坐标正算
文件名
4
坐 标 正 算
22步
L1
ProgB:
X=Rec(L,B)+X▲Y=Y+W
说明:
①程序根据已知坐标、边长及其坐标方位角计算待定点坐标;
②程序中:
B=坐标方位角、L=边长、X=已知点纵坐标、Y=已知点横坐标;
③算例:
按Prog3调用程序显示B?
输入坐标方位角100.11431,按EXE键、显示L?
输入边长151.213,按EXE键、显示X?
输入已知点纵坐标1000,按EXE键、显示待定点纵坐标973.2346789,按EXE键、显示Y?
输入已知点横坐标2000,按EXE键、显示待定点横坐标2148.825364。
程序10——坐标方位角交会计算
文件名
5
坐 标 正 算
100步
L1
Lbl0:
{I}:
ProgD:
S=B+I:
A"X0"C"YO"XYN:
B=N:
ProgB:
E=(AtanS-X
tanB+Y-C)/(tanS-tanB▲C+tanS(E-A▲Pol(A-E,C-Ans▲B=B-W:
ProgA:
ProgC▲B=F:
Goto0
说明:
①程序利用已知方位角,计算测站上任意方向线与轴线相交的距离及交点坐标和该交点作为测站后与轴线的夹角;
②如下图,程序中:
B=已知方位角、I=观测角、A(X0)=测站纵坐标、C(Y0)=测站横坐标、X=轴线上任意点纵坐标、Y=轴线上任意点横坐标、N=轴线方位角;
③算例:
按Prog5调用程序显示B?
输入已知方位角150.01252,按EXE键、显示I?
输入观测角110.30,按EXE键、显示X0?
输入测站纵坐标100,按EXE键、显示Y0?
输入测站横坐标100,按EXE键、显示X?
输入轴线方向上任意点纵坐标80,按EXE键、显示Y?
输入任意点横坐标150,按EXE键、显示N?
输入轴线方位角210.12374,按EXE键、显示交点纵坐标111.3973848,按EXE键、显示交点横坐标168.2813366,按EXE键、显示测点至交点的距离69.22601608,按EXE键、显示该交点作为测站后与轴线的夹角09.41122,按EXE键、显示I?
输入观测角……
④程序所用正切公式如图所示:
已知两方向、两点坐标,求其交点P的坐标;
程序11——垂距计算
文件名
6
垂 距 计 算
27步
L1
ProgB:
D=B:
Prog2▲B=W-D:
VcosB▲VsinB
说明:
①程序用以计算勘探点偏离轴线的距离;
②如图,程序中:
B=轴线坐标方位角、A(X0)=轴线控制点的纵坐标、C(Y0)=轴线控制点的横坐标、X=勘探点纵坐标、Y=勘探点横坐标;
③算例:
按Prog6调用程序显示B?
输入轴线坐标方位角301.09505,按EXE键、显示X0?
输入轴线上控制点的纵坐标138.579,按EXE键、显示Y0?
输入轴线上控制点的横坐标1321.474,按EXE键、显示X?
输入勘探点的纵坐标341.871,按EXE键、显示Y?
输入勘探点的横坐标986.408,按EXE键、显示控制点至勘探点的距离391.914355(斜距),按EXE键、显示控制点至勘探点的坐标方位角301.1446034,按EXE键、显示控制点至勘探点的距离(平距)391.9139527,按EXE键、显示勘探点偏离轴线的距离0.561530147(垂距);
④如果勘探点偏离轴线的距离为正值,表示该点在轴线右侧;若为负值,则该点在轴线左侧。
程序12——视差法测距计算
文件名
7
视 差 法 测 距 计 算
30步
L1
ProgD:
KTH:
1/tan(B/2▲H=AnstanI+H+K-T
说明:
①程序用以计算2M横基尺视差法测距的边长及高程;
②程序中:
B=视差角、I=垂直角、K=仪高、T=觇标高、H=测站高程;
③算例:
按Prog7调用程序显示B?
输入视差角2.33176,按EXE键、显示I?
输入垂直角1.45386,按EXE键、显示K?
输入仪高1.432,按EXE键、显示T?
输入觇标高1.405,按EXE键、显示H?
输入测站高程443.425,按EXE键、显示测站至测点的距离44.84444364,按EXE键、显示测点高程444.8305217。
程序13——剖面交点计算
文件名
8
剖 面 交 点 计 算
55步
L1
ProgB:
L:
D=B:
B=J:
ProgD:
N=D-B:
M=I-D:
Abs(LsinM/sin(M+N▲Abs(L
sinN/sin(M+N
说明:
①程序8用以计算剖面交点距离;
②输入坐标方位角时不用考虑方向性;
③算例:
如图,按Prog8调用程序显示B?
输入坝1的坐标方位角96.0303或276.0303,按EXE键、显示L?
输入坝1交坝Ⅱ围Ⅰ的距离21.9,按EXE键、显示J?
输入坝Ⅱ的坐标方位角301.0950或121.0950,按EXE键、显示I?
输入围Ⅰ的坐标方位角307.1202或127.1202,按EXE键、显示坝Ⅱ交围Ⅰ及坝1的距离107.7198781,按EXE键、显示围Ⅰ交坝1坝Ⅱ的距离88.38012441。
程序14——前方交绘计算
文件名
9
前 方 交 绘 计 算
65步
L1
A"X0"C"Y0":
ProgD:
M=1/tanB:
N=1/tanI:
(AN+XM-C+Y)/(M+N▲(CN+
YM+A-X)/(M+N
说明:
①程序9用以计算前方交会法中待定点的坐标;
②程序采用余切公式计算坐标,三角形点号应按逆时针编排。
若是顺时针排列,观测角应以负值输入;
③算例,如图表所示:
逆 时 针
顺 时 针
Xa
2235.6
Ya
845.0
Xa
6735.7
Ya
10738.5
Xb
6735.7
Yb
10738.5
Xb
2235.6
Yb
845.0
∠cab
48.5221
∠cba
22.4045
∠cab
-22.4045
∠cba
-48.5221
Xc
6467.705
Yc
2112.188
Xc
6467.705
Yc
2112.188
④交会角∠C在30°~50°之间时交会误差最小。
程序15——附合导线近似平差
文件名
10
附 合 导 线 近 似 平 差
315步
L1
C=0:
A>9Defm2(A-9△Lbll:
C=C+1:
{I}:
Prog3:
A[C+8]=D:
C≠A
Goto1△{B}:
ProgB:
B=D-B:
D=B/A:
ProgC▲C=0
L2
Lbl2:
C=C+1:
A[C+8]=A[C+8]-DC:
B=A[C+8:
ProgC▲C≠AGoto2△C=
8:
D=0:
E=0:
F=0
L3
Lbl3:
C=C+1:
{B}:
D=D+B:
E=E+BcosA[C:
F=F+BsinA[C:
A[A+C]=B:
C≠
A+7Goto3△D▲{BI}:
BGHI:
E=E+B-H▲F=F+G-I▲D/√(E2+F2▲E=E/D
:
F=F/D
L4
C=8:
Lbl4:
C=C+1:
B=A[A+C]cosA[C]+B-EA[A+C▲G=A[A+C]sinA[C]
+G-FA[A+C▲C≠A+7Goto4△Defm0
说明:
①程序10用以单一附合导线的近似平差计算;
②程序中:
A=测站数、B=起算边坐标方位角、附合边坐标方位角、导线边边长、起点纵坐标、I=观测角(左角)、终点横坐标、G=起点横坐标、H=终点纵坐标;
③程序设计可以计算9个测站,大于是9站时只要计算器剩余有步骤数,程序会自动扩充储存器增加储存变量;
④算例如表:
按Prog10调用程序显示A?
输入测站数5,按EXE键、显示B?
输入起算边的坐标方位角169.5359,按EXE键、显示I?
输入B点的观测角93.3642,按EXE键、显示I?
输入C点的观测角179.56,依次输入各观测角…,按EXE键、显示B?
输入附合边的坐标方位角250.06,按EXE键、显示导线的角度闭合差-0.003699999(-37″),按EXE键、显示平差后B点至C点的坐标方位角83.30484,按EXE键依次显示各导线边平差后的坐标方位角…,按EXE键、显示B?
输入B点至C点的边长126.11,按EXE键、显示B?
依次输入各导线边的边长…,按EXE键、显示各边长之和508.96,按EXE键、显示B?
输入起点纵坐标801.53,按EXE键、显示G?
输入起点横坐标724.84,按EXE键、显示H?
输入终点纵坐标959.29,按EXE键、显示I?
输入终点横坐标1158.41,按EXE键、显示纵坐标闭合差0.093449934,按EXE键、显示横坐标闭合差-0.02045631,按EXE键、显示导线全长相对闭合差5320.360577(K=1/5320),按EXE键、显示C点纵坐标815.754,按EXE键、显示C点横坐标850.148,按EXE键依次显示D、E、F、G点的纵横坐标。
点号
观测角
坐标方位角
边长
纵坐标X
横坐标Y
A
169 53 59.0
B
93 36 42
801.530
724.840
83 30 48.0
126.11
C
179 56 00
815.754
850.148
83 26 55.8
126.92
D
180 03 18
830.211
976.244
83 30 21.2
129.65
E
121 28 36
844.850
1105.067
24 59 04.0
126.28
F
45 06 48
959.290
1158.410
250 06 00.0
G
备
注
N=5 fβ=-37″ ΣD=508.96 fx=+0.093 fy=-0.020f=√(fx2+fy2) K=1/(Σd/f) k=1/5350
角度闭合差容许值:
(图根导线
)
导线全长相对闭合差:
K≤1/2000
程序16——渠道中线里程计算
文件名
11
渠 道 中 线 里 程 计 算
164步
L1
K=0:
N=0:
A"X0"C"Y0":
Lbl1:
{XY}:
K=K+1:
Pol(X-A,Y-C▲B=W:
Porg
A:
P=B:
PorgC▲A=X:
C=Y
L2
ProgF:
D=Ans:
K=1J=D△K≥2B=N-P:
B≥180B=B-360△B≤-
180B=B+360△G=B:
J▲ProgC▲ProgG:
B=P+3G/4:
ProgC▲B=P+G/4:
Prog
C▲△N=P:
Goto1
文件名
F
取 位 子 程 序
39步
L1
U=Frac10V/10:
10U≥.5V+.1-U△10U<.5V-U△
文件名
G
曲 线 要 素 计 算 子 程 序
101步
L1
{R}:
V=Rtan(AbsG/2:
ProgF:
T=Ans▲V=AbsGRπ/180:
ProgF:
W=5(Frac
V-U
L2
FracW≠0V=V+.05△FracW=0V=V-.05△ProgF:
L=Ans▲V=R/cos
(AbsG/2)-R:
ProgF▲V=J-T▲V+L/2▲V+L▲J=Ans+D-T
说明:
①程序11用以渠道中线测量的里程及其曲线要素的计算
②程序中:
A=渠线起点IP0的纵坐标、C=渠线起点IP0的横坐标、X=IP1…IPn的纵坐标、Y=IP1…IPn的横坐标、R=半径;
③程序计算时会自动取舍弧长L的小数位为偶数,其它数据取位后均无余留尾数,用以消除累加误差;
④程序设计只保留1位小数,若要保留2位小数,应修改取位子程序中的语句U=Frac10V/10为U=Frac100V/100;修改取位子程序中的语句:
10U≥.5V+.1-U△10U<.5V-U△为:
100U≥.5V+.01-U△100U<.5V-U△;修改曲线要素计算子程序中的语句W=5(FracV-U为W=5(Frac10V-10U;修改语句FracW≠0V=V+.05△FracW=0V=V-.05△为FracW≠0V=V+.005△FracW=0V=V-.005△
⑤若要保留3位小数,应修改取位子程序中的语句U=Frac10V/10为U=Frac1000V/1000;修改取位子程序中的语句10U≥.5V+.1-U△10U<.5V-U△为1000U≥.5V+.001-U△1000U<.5V-U△;修改曲线要素计算子程序中的语句W=5(FracV-U为W=5(Frac100V-100U;修改语句FracW≠0V=V+.05△FracW=0V=V-.05△为FracW≠0V=V+.0005△FracW=0V=V-.0005△;
⑥算例见下表,按Prog11调用程序显示X0?
输入IP0的纵坐标,按EXE键、显示Y0?
输入IP0的横坐标,按EXE键、显示X?
输入IP1的纵坐标,按EXE键、显示Y?
输入IP1的横坐标,按EXE键、显示IP0至P1的坐标方位角,按EXE键、显示IP0至IP1的距离,按EXE键、显示X?
输入IP2的纵坐标,按EXE键、输入IP2的横坐标,按EXE键、显示IP1至IP2的坐标方位角,按EXE键、显示IP1至IP2的距离,按EXE键、显示IP1的里程J,按EXE键、显示IP1的外偏角α,按EXE键、显示R?
输入半径,按EXE键、显示切线长T,按EXE键、显示弧长L,按EXE键、显示外矢距E,按EXE键、显示直圆点(ZY)里程,按EXE键、显示曲中点(QZ)里程,按EXE键、显示圆直点(YZ)里程,按EXE键、显示直圆点(ZY)至曲中点(QZ)的坐标方位角,按EXE键、显示曲中点(QZ)至圆直点(YZ)的坐标方位角,按EXE键、显示X?
依次输入IP3…IPn的纵横坐标……;
⑦计算结束后按AC、RCL、M+键将显示终点IP桩的里程。
⑧如果外偏角α显示为正值,并且纵坐标X是递增的,则该角位于前进方向的左侧;反之纵坐标X是递减的,遇该角位于前进方向的右侧。
如果外偏角α显示为负值,并且纵坐标X是递减的,则该角位于前进方向的左侧;反之纵坐标X是递增的,则该角位于前进方向的右侧。
IP
桩
点
号
X
Y
坐标方位角
边长
IP桩
里程
JD
外偏角
半
径
R
切
线
长
T
弧
长
L
外
矢
距
E
0
54110
22903
158 52 09.0
396.672
1
53740
23046
396.7
+47 53 44.4
60
26.6
50.2
5.7
110 58 24.6
128.515
2
53694
23166
522.2
-29 33 14.6
60
15.8
31.0
2.1
140 31 39.1
572.591
3
53252
23530
1094.2
-36 29 22.2
60
19.8
38.2
3.2
177 01 11.4
461.624
4
52791
23554
1554.4
IP
桩
点
号
直圆点
ZY里程
曲中点
QZ里程
圆直点
YZ里程
直圆点~曲中点
坐标方位角
曲中点~圆直点
坐标方位角
0
1
370.1
395.2
420.3
146 53 42.9
122 56 50.7
2
506.4
521.9
537.4
118 21 43.2
133 08 20.5
3
1074.4
1093.5
1112.6
149 39 02.2
167 53 48.3
4
计算公式:
T=R×tg(α/2) L=π×α×R E=R/(cos(α/2))-R校核:
q=2T-L JD=QZ+q/2
程序17——多边形面积计算
文件名
12
多 边 形 面 积 计 算
100步
L1
A:
XY:
M=X:
N=Y:
C=0:
S=0:
D=0:
B=0:
Lbl1:
C=C+1:
S=S+3(-BX+DY)/
4000:
D=X:
B=Y:
C=A→X=M:
Y=N:
Goto1△C≠A+1→{XY}:
Goto1△S
说明:
①程序12采用多边形顶点坐标计算多边形的面积;
②程序中A=多边形边数(或顶点数)、X=多边形各顶点的纵坐标、Y=多边形各顶点的横坐标;
③计算时如果是按顺时针方向输入多边形各顶点的纵横坐标,则面积显示为正值,反之则为负
④面积是以亩为单位,亩与平方米的换算常数是:
2000/3,即1亩=2000/3M2、1M2=3/2000亩;
⑤算例:
按Prog12调用程序显示A?
输入多边形的边数(或顶点数)4,按EXE键、显示X?
输入顶点1的纵坐标10,按EXE键、显示Y?
输入顶点1的横坐标10,按EXE键、显示X?
输入顶点2的纵坐标40,按键、显示Y?
输入顶点2的横坐标10,按EXE键、显示X?
输入顶点3的纵坐标40,按EXE键、显示Y?
输入顶点3的横坐标60,按EXE键、显示X?
输入顶点4的纵坐标10,按EXE键、显示Y?
输入顶点4的横坐标60,按EXE键、显示四边形的面积2.25(亩)。
⑥坐标法计算多边形面积的公式如下:
(1)
(2)
其中:
i=1、2、3、4…n
公式
(1)和
(2)中:
S=多边形面积、n=图形界址点个数、
Xi=界址纵坐标、Yi=界址点横坐标;当i=1时,i-1=n;i=n时,i+1=1
程序18——闭合导线近似平差
文件名
13
闭 合 导 线 近 似 平 差
376步
L1
C=-1:
E=0:
A