《计算器工程程序大全》由会员分享，可在线阅读，更多相关《计算器工程程序大全（34页珍藏版）》请在金锄头文库上搜索。
1、SOUTH南方测绘 CASIO工程简明程序集卡西欧（CASIO）工程简明程序散集地 址：杭州市文三路111-2号 电 话:0571-88061065 13515713388 13587304710(绍兴办事处)传 真: 0571-88057475 E-mail：编 者：吴 东目录 前言 02第一章 专业篇 1-1 大地测量 03 1-2 工程放样 061-3 圆与曲线 101-4 平差计算 18第二章 生活篇2-1 不翻日历便知某年某月某日是星期几182-2 连加（减）程序 192-3 量算面积程序 19第三章 经验篇3-1 Fx-4500PA和Fx-4800P之比较 213-2 编程经验初探
2、 22第四章 信息篇4-1 测绘网站 224-2 卡西欧新品CASIO Fx-4850P 234-3 CASM的产业化 244-4 我的测绘理念 28前 言 电子技术的发展带动了测绘技术的进步，繁杂的测量工作也由此变得轻松起来。卡西欧（CASIO）公司一直致力于计算器的研发，对我们测界人士来说，理推Fx-4500PA和Fx-4800P程序型计算器。它们携带方便，存贮量大，实为测绘工作者不可缺少的辅助工具。笔者从事测绘一线工作多年，和4500PA和4800P结下了深情厚谊，也从中摸索出很多所谓的诀窍，今整理出来，供同行参考斧正，同时也作为2004年“南方公司新产品展示会暨新技术发布会绍兴”的一份
3、献礼。本册子分四大块：专业篇，生活篇，经验篇，信息篇。从章节安排上，突出新颖；内容安排上，力求精简。由于篇幅关系，所有程序略去了数学模型，按程序内容、字符说明、功能介绍、示例示图等进行阐述。部分略去的程序恕没列上，感兴趣的读者可与笔者本人联系。书中所有程序笔者释为至宝，已在多次生产实践中得到应验，读者可放心使用。笔者致力于程序的优化和实用，读者可在具体使用中体会到这些。限于业务水平，纰漏之处难勉，恳请同行多多批评赐教。 吴 东 二00四年六月于中国杭州 第一章 专业篇 1-1 大地测量1.1高斯正形投影正算程序1.1.1 高斯正形投影正算程序（4500PA）：主程序：XYL1 P=206265
4、: Lbl 0:BLJZ:BLJ”L0”Z:U=B:Prog2:B=U:Z=5=Prog 5:Prog 4L2 U=L:Prog 2:L=U:U=J:Prog 2:J=U:T=tanB:G=F.5cosB:N=A/(1E(sinB)2).5:K=3600(LJ)/PL3 X=Q+NK2sin(2B)/4+Nsin(2B)(cosB)2(5T2+9G2+4G4)K4/48+N(cosB)4(6158T2+T4)K5sin(2B)/1440L4 Y=NkcosB+N(1T2+G2)(KcosB)3/6+N(518T2+T4+14G258(GT)2) (KcosB)5/120:XYGoto 0子程序
5、：2U=IntU+5FracU/3+Frac100U/90子程序：4Q=111134.8611B(32005.7799sinB+133.9238(sinB)3+.6973(sinB)5+.0039(sinB)7)cosB:A=6378245:E=.006693:F=.006738子程序：5Q=111133.00478B(32009.8575sinB+133.9602(sinB)3+.6976(sinB)5+.0039(sinB)7)cosB:A=6378140:E=.006694:F=.006739字符说明：1B-纬度，Z-调节数（1975年椭球(5)和克氏椭球(4)的选择），L-经度，J-中
6、央子午线经度2、4500PA在本程序运行中速度很慢。3、本程序约446STEPS。4、为了达到计算精度，小数点位数要设置到9位。示例：B:29355.5817L:1065159.5438（克氏椭球） Fom-1序号字符输 入结 果序号字符输 入结 果1B29350558174Z42L106515954385X3275611.1873J1056Y180859.668说明：由于4500PA的精度不够，计算结果跟实际有点小差异。 本例的实际计算结果：X：3275611.188Y: 180859.8681.1.2高斯正形投影反算程序（4500PA）：主程序：LBL1 Lbl 0:XYK:XYK：X=X
7、/106:K=5=A=6378140:E=.006694:F=.00674:D=6399596.652:Prog 5:A=6378140:E=.006693:F=.006738:D=6399698.90178:Prog 4L2 T=tanC:G=FcosC:N=A/(1E(sinC)2):V=D/N:H=Y/N: B=C90V2T/(H2H4(5+3T2+G29(GT)2)/12+H6(61+90T2+45T4)/360)L3 L=180/(cosC)(HH3(1+2T2+G2)/6+H5(5+28T2+24T4+6G2+8(GT)2)/120):BLGoto 0子程序：4X3000=Prog
8、 44:C=9.04353301X.0000005X2.00075311X3.00000084X4.00000426X5.0000001X6子程序：44C=27.11115372+9.02468257(X3).0057974(X3)2.00043532(X3)3+.00004857(X3)4+.00000216(X3)5.00000019(X3)6子程序：5X3000=Prog 55:C=9.04369066X.0000005X2.00075326X3.00000084X4+.00000426X5.0000001X6子程序：55C=27.11162289+9.02483658(X3).0057
9、9851(X3)2.0004354(X3)3+4.85810 5(X3)4+.00000216(X3)5说明：1、 XY-纵横坐标，K-调节数（1975年椭球(5)和克氏椭球(4)的选择）,指根号.2、 显示结果中，B为纬度，L为经差。后者要加上中央经度才是坐标的经度。3、 本程序在经差3时，精度较高。4、 由于计算器本身精度的限制，计算结果与实际有差异，但完全可满足GPS测量的需要。5、 本程序约637STEPS。示例：X：3275611.188,Y: 180859.868（克氏椭球） Fom-2序号字符输 入结 果序号字符输 入结 果1X 3275611.1884 B 29.5816736
10、2(293454.03)2Y 180859.8685 L1. 866511677(15159.44)3K 4 说明：由于4500PA的精度不够，计算结果跟实际有点小差异。 本例的实际计算结果：B：29355.58L: 15159.5438 1-2 工程放样1.1工程放点计算(4500PA)主程序名:FWL1 ABCD：Pol(AC,BD):W=W+180L=VZ=W:Lbl 0:JR:M=A:N=B:O=J:P=R:Prog 1:F=WZ:F0=W=F:W=360+FWGoto 0子程序：1L1 MNOP：Pol（MO，NP）：W=180+WL=V说明：1、 AB指测站坐标，CD指定向坐标，J
11、R指放点坐标，W=W+180指放点方位，L=V指放点平距，W指放点方向和定向方向间的夹角。2、 子程序1可专门用来算任意两点间的方位和平距。示例示图：【Fig-1】 Fom-3序号字符输入结果序号字符输入结果1A5007J4802B5008R7003C5809W=180+W957106（95-42-38）4D60010L=V2009985W=180+W513402（51-20-25）11W443704（44-22-13）6L=V1280621.1.1点的加密计算(4500PA)主程序名：JML1 GIUZ：Pol(GU,IZ):W=180+WVA=WL2 Lbl 0:WUZN:WU“UT”Z“
12、ZS”:Prog 9:M=U:X=G+Zcos(W+M) Y=I+Zsin(W+M) N:N0=G=X:I=Y:W=180+W+M:Goto 0子程序：9U=IntU+5FracU/3+Frac100U/90说明：1、 GIUZ指任意方向上两点的坐标。2、 放点置方位时，W可输为零，U为加密方向上的方位值；放点置零时，W可直接按过，U为加密方向相对于起始方向的角度。A始终是起始方向的方位角值，Z为所加密点和起始点的平距。3、 N不为零说明此方向上在连续单向加密；N为零说明是非连续单向加密。4、 加密点不外乎1）本方向上的点 2）垂直方向上的点 3）非上两种情况的点示例示图：【Fig-2】本例为W直接按过，并有一个连续加密的情况。 Fom-4序号字符输入结果序号字符输入结果1G5009ZS1502I50010X55958573U58011
