CASIO4500坐标计算程序.docx

1、CASIO4500坐标计算程序以文本方式查看主题-合铁工校(-=专业交流=(-CASIO4500计算程序(-CASIO4500坐标计算程序 HUANHEQUXIANKJDDLS：B=D2/24R：M=D/2-DB/10RT=M+tan.5A(R+B)L=RA/180+DE=(R+B)/cos.5A-RZZH=K-THHY=Z+DQQZ= Z+L/2JYH=Z+L-DOHZJ+Dprog 1NN=1V:P:prog 3:L=W:FixmLb1 2:S:SKX:SJGoto4SHGoto3U=S-Z:E=U-UXY5/40R2D2:F=UXY3/6RD:Goto5 Lb1 3:U=90(2S-2H

2、+D)/ R:E=RsinU+M:F=R-RcosU+B:Goto5 Lb1 4:F=O-S:D0I=30F2/RDV=F-FXY3/90R2:N=1U=360-A:F=180-A+I: U=A:F-180+A-IE=T+TcosU+VcosF:F=TsinU+VsinF:Goto1 Lb1 5:N=1F=-FLb1 1:X=E:Y=F:P=1prog2: prog4: prog3:W=W-L:W0W=W+360V: SW: RGoto2 X=C+EcosL-FsinLY=G+EsinL+FcosLCX0GY0:Pol(X-C,Y-G):W0W=W+360W D0Q90U2/RDFixm:I=

3、A-3I:N=1Q=-Q:U=-U:I=-ISHF=L+QSJF=L+I: F=L+UV=1: E:EB:FR=F+E:Lb1 6: I:ISL:X+IcosFY+IsinFV2V=V+1：Goto6ZHIXIAN EX0FY0AR0KCZ:Lb1 1: S=SKX:D=S-K:X=E+DcosAY=F+DsinAV=1:Lb1 2: BI:H=A+B-180:ISL:X+IcosHY+IsinHV2V=V+1:Goto2Goto1计算要素：JD交点里程 LS缓和曲线长R 圆曲线半径 A 线路转角T 切线长 L 圆弧总长度E 外矢距 N 曲线方向，左偏取“1”，右偏取“0”P 取“1”X0.Y

4、0ZH点坐标X.Y交点坐标KX待求点里程B 与中线夹角 SL边距，左“+”；右“-”-作者：admin-发布时间：2005-8-15 9:43:10-CASIO4800超高计算程序 超高及高程（源程序）L1: R:L“LS”:Z“ZH”:H“HZ”:V“RS”:F“I1”:G“I2”:O“BPH”:N“BP”:T=Vabs(GF)2:X“X(R+1,L-1)”:P“P(R+,L-1)”L2: Lbl 6 : SL3: S+TN0 D=O+(SN)F: SN0 D=O+(SN)FU(SN+T)22V: SNT D=O+(SN)GU(TS+N)22V: D=O+(SN)G D“ H=” L4: R

5、 E= : M= Goto 1: R0 Goto 3 : Goto 4L7: Lbl 4 : S Z I=0.02 : Goto 5: SZ+C I=0.02(0.02EPX)C(SZ) : Goto 5: SHC I=EPX : Goto 5: SH I=EPX+(0.02EPX)C(SH+C) : Goto 5: I=0.02 : Goto 5L8: Lbl 3 : Z1=0.04(EPX+0.02)C : S Z + Z1 I=0.02 : Goto 5: SZ+C I=0.02+(EPX0.02)(CZ1)(SZZ1) : Goto 5: SHC I=EPX : Goto 5: S X

6、Y：2.XY = SZ：N：UX0：VY0：OS0：GF0：HLS：PR0：RRN：Q：C=1P：D=(P-R)(2HPR)：E=180：N=1=Goto 1：Goto 2Lbl 1：SZ：SZ：W=Abs(S-O)：Prog SUB1：XXS=XYYS=YGoto 1Lbl 2：XY：XY：I=X：J=Y：Prog SUB2：SS=O+WZZ=ZGoto 22. 正算子程序(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+Q

7、ELW(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+ZsinF3. 反算子程序(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 0Lbl 1：Z=

8、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次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相

9、接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。 (6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。2、输入与显示说明 输入部分：1. SZ = XY2. XY = SZ N ?选择计算方式，输入1表示进行由里程、边距计算坐标 ；输入2表示由坐标反算里程和边距。 X0 ？线元起点的X坐标 Y0 ？线元起点的Y坐标 S0 ？线元起点里程 F0 ？线元起点切线方位角 LS ？线元长度 R0 ？线元起点曲率半

10、径 RN ？线元止点曲率半径 Q ？线 元左右偏标志(左偏Q=-1，右偏Q=1，直线段Q=0）S ？ 正算时所求点的里程Z ？正算时所求点距中线的边距(左侧取负，值右侧取正值，在中线上取零)X ？反算时所求点的X坐标Y ？反算时所求点的Y坐标 显示部分： XS=正算时，计算得出的所求点的X坐标 YS=正算时，计算得出的所求点的Y坐标 S=反算时，计算得出的所求点的里程 Z=反算时，计算得出的所求点的边距四、算例某匝道的由五段线元（直线+完整缓和曲线+圆曲线+非完整缓和曲线+直线）组成，各段线元的要素（起点里程S0、起点坐标X0 Y0、起点切线方位角F0、线元长度LS、起点曲率半径R0、止点曲率

11、半径RN、线 元左右偏标志Q）如下：S0 X0 Y0 F0LS R0 RN Q500.00019942.83725643.561 125 16 31.00269.256 1E45 1E450769.25619787.34028563.378 125 16 31.0037.492 1E45 221.75 -1806.74819766.56628594.574 120 25 54.07 112.779221.75 221.75 -1919.52719736.07228701.893 91 17 30.63 80.285221.75 9579.228 -1999.81219744.03828781.

12、659 80 40 50.00 100.0001E45 1E45 0 1、正算（注意：略去计算方式及线元要素输入，请自行根据所求点所在的线元输入线元要素）S=700Z=-5计算得XS=19831.41785 YS=25809.72590S=700Z=0 计算得XS=19827.33592 YS=25806.83837S=700Z= 5计算得XS=19823.25398 YS=25803.95084S=780Z=-5计算得XS=19785.25749 YS=28575.02270S=780Z=0 计算得XS=19781.15561 YS=28572.16358S=780Z= 5计算得XS=197

13、77.05373 YS=28569.30446 S=870Z=-5计算得XS=19747.53609 YS=28654.13091S=870Z=0 计算得XS=19742.68648 YS=28652.91379S=870Z= 5计算得XS=19737.83688 YS=28651.69668 S=940Z=-5.123计算得XS=19741. 59118YS=28722.05802S=940Z=0计算得XS=19736.47687 YS=28722.35642S=940Z= 3.009计算得XS=19733.47298 YS=28722.53168 2、 反算 X=19831.418 Y=2

14、5809.726 计算得S=699.9999974Z= -5 .00018164 X=19827.336 Y=25806.838 计算得S=699.9996493Z= 0.000145136 X=19823.25398 Y=25803.95084计算得S=699.9999985Z= 5.000003137 X=19785.25749 Y=28575.02270计算得S=780.0000035Z= -5 .000001663 X=19781.15561 Y=28572.16358 计算得S=780.0000025Z=- 0.000002979 X=19777.05373 Y=28569.3044

15、6 计算得S=780.0000016Z= 4.99999578 X=19747.536 Y=28654.131 计算得S=870.0001137Z= -4.99941049 X=19742.686 Y=28652.914 计算得S=870.0003175Z=- 0.00041814 X=19737.837 Y=28651.697 计算得S=870.0002748Z= 4.999808656X=19741.5912 Y=28722.0580 计算得S=939.9999786Z= -5.123024937X=19736.4769 Y=28722.3564 计算得S=939.9999862Z=- 0

16、.000027710X=19733.4730 Y=28722.5317 计算得S=940.0000238Z= 3.00898694-作者：花下独眠-发布时间：2006-2-15 20:42:58-看到斑竹为大家提供这么好的东西，非常感动，本人毕业后从事一段时间测量工作，有几个实用的Casio Fx-4500PA计算器程序，现场使用比较方便，也在此献丑！绝对好用，特别是曲线任意坐标计算！Casio Fx-4500PA计算器程序库一、整机概况Fx-4500PA计算器是日本Casio公司的产品，它除了具备函数型计算器的全部功能之外，还具备如下特殊功能： 1、二进制和八进制字组显示； 2、二进制、八进

17、制、十进制和十六进制的转换和加、减、乘、除运算。 3、进行f（x）函数的输入与积分的计算； 4、进行标准差和回归计算。 Fx-4500PA计算器具有同时显示运算表达式和计算结果的双行式显示的特点。该计算器具有记存1103个程序步的记忆储存装置，计算器储存了程序之后，只要您输入数据，再按一下EXE键，计算器就会用数据进行程序执行,对于随数据而变化的反复计算的运用十分方便。 Fx-4500PA计算器具备条件指令转移和非条件指令转移的能力，还具备逻辑判断能力，因此它可以进行测量工作中比较复杂的计算，当然其它专业的计算工作也同样可以解决。 Fx-4500PA计算器操作方法比较简单、易学，程序编排灵活、

18、易懂，可称之谓各专业技术人员的物美、价廉、得心应手的计算工具。 每一台Fx-4500PA计算器都具备一本说明书，该计算器各部键的功能以及操作方法，只要认真按照说明书各章节学习，很快就能掌握和运用，现将我在测量实践中应用的部分计算程序介绍给同行以便共同探讨。二、程序库 = 1 = 坐标反算程序（已知两坐标求算其平距和方位角） 程序名：ZBFS Norm：A”X1”：B”Y1”：pol(C”X2” - A，D”Y2” - B：”S=”WW=W+360W：”A=” 注意：本程序在CASIO4500上调试运行通过。如果是在4800上运行请将源程序中的V，W分别用I，J替代即可。 运行后界面如下： 输入

19、项目: X1？请输入点1的X坐标，按EXE键输入下一个量。 Y1？请输入点1的Y坐标。 X2？请输入点2的X坐标。 Y2？请输入点2的Y坐标。 输出项目: S为两点之间距离，按EXE键输出下一个量。 A为1点到2点的方位角。按EXE键循环。 = 2 = 坐标正算程序（已知一坐标及其平距和方位角求算另一点坐标） 程序名：ZBZS C”X1”D”Y1”S”S=”R”A”：X”X2”=C+Rec(S，RY”Y2”=D+W 注意：本程序在CASIO4500上调试运行通过。如果是在4800上运行请将源程序中的V，W分别用I，J替代即可。 运行后界面如下： 输入项目 X1？ 请输入点1的X坐标，按EXE键

20、输入下一个量。 Y1？ 请输入点1的Y坐标。 S=？ 为两点之间距离。 A？为1点到2点的方位角。 输出项目 X2输出点2的X坐标。 Y2输入点2的Y坐标。按EXE键进入下一个循环。 = 3 = 竖曲线计算程序 程序名：SQX LBl 1：DEIHRQPJ：W=DEIHRQPJ：T=(I-H)R/2：N=T2/2R LBl 2：A：Z=A LBl 3：AP=Goto 4：AE=Goto 6：AQ=Goto 5：X=D+(A-E)I Goto 2 LBl 4：X= D+(A-E)HGoto 2 LBl 5：M=(Q-A)2/2RJ：X=D+(A-E)I+M Goto 2 LBl 6：M=(P-A

21、)2/2RJ：X=D+(A-E)H+M Goto 2 注意：本程序在CASIO4500上调试运行通过。如果是在4800上运行请将源程序中的V，W分别用I，J替代即可。 字母含义： D：变坡点高程 E：变坡点里程J：曲线凸凹（凸为-1、凹为+1） L：开始坡度 H：结束坡度 A：待求点标高R：竖曲线半径 Q：起点里程 P：终点里程 = 4 = 计算曲线上任一点放样极坐标 F1：QXJS L1Abs WFixmC=7=Fixm:Goto FFix 4:Z=E7AA0/Abs AN9:A=Abs A:Q=R+SLsci2/24R L2T=S/2-Sxy3/240R2+Anstan .5AL=RA/1

22、80 +SHZH=JJDT-TVHY=H+SVQZ=H+.5LVYH=H+L-SLHZ=H+LDxJKyJ BGraOXEUYEGXmMXm L3LbI F:FN:FH=X=F-J:ZN:FL-S=I= L-S=I= F-H:FProg Q:Rec(R-ZN,90(2I-S)/R: X=W-Qtan.5A:Q-VI=L-F:PoI(T-I,ZN:FProg Q:PoI(-X,AnsX=Rec(Ans,A+W L4WX=D+Rec(PoI(X,ZAns,W+B:Y=K+W: EnD OK!:PoI(G-O,M-U:I=W:PoI(X-O,Y-U:W-I: AnsAns+360I=Ans:60Frac Ans:IDeg=IntI+.01Int Ans+6mFrac AnsVP=VGoto F 主程序结束 F2：Q子程序 L1Y=2RS:W=180I2/Y:X=I-Ixy5/10Y2-T-ZNsin W:Ixy3/3Y-Ixy7/42Yxy3+