4800全套计算程序.doc

CASIO fx4800P 程序集 公路工程 杨小杰 攀枝花公路建设有限公司 2005.10 整套卡西欧计算器CASIO FX-4800 4500测量程序各位同行：我长年从事公路桥梁测绘工作，收集整理了一整套卡西欧计算器CASIO FX-4800P和卡西欧CASIO FX-4500P测量放线程序，该套程序经过反复修正，具有极高的可靠性。现在免费提供给大家.程序目录：1.1测边大地四边形的条件平差程序1.2测边中点多边形的条件平差程序1.3单一导线的严密平差程序1.4单一结点导线的严密平差程序1.5边、角后方交会的条件平差程序1.6测边交汇点的条件平差程序1.7单一水准路线的平差程序1.8单结点的平差程序1.9多结点水准网的平差程序1.10多边形水准网的平差程序1.11三角高程差的计算程序1.12测距边水平距离的计算程序1.13高斯投影座标正、反解和换带的计算程序1.14平面座标转移的计算程序2.1单边、角后方交会的计算程序2.2双边、角后方交会的计算程序2.3碎部点座标的计算程序2.4内、外分点座标的计算程序2.5房角标的计算程序2.6多次后方测角交会点2.7测边交会点的计算程序2.8导线三维值的计算程序2.9圆心座标和半径的程序2.10测站点、照准点归心改正的计算程序2.11单边、双站后方交会的计算程序3.1线路中心座标的计算程序3.2线路边桩座标的计算程序3.3红线交点座标的计算程序3.4直线与圆曲线交点座标的计算程序3.5垂足座标的计算程序3.6道路不等宽时曲线要素的计算程序3.7缓和曲线敷设时的计算程序3.8圆曲线敷设时的计算程序3.9极座标放线线的计算程序3.10交点座标的计算程序3.11宗地有关项目的计算程序3.12丘地有关项目的计算程序4.1导线X、Y、Z值的近似平差程序4.2无定向导线的计算程序4.3测角前方交会点座标的计算程序4.4测边交会座标的计算程序4.5测角后方交会点座标的计算程序4.6双边、角后方交会点座标的计算程序4.7平面座标转移的计算程序4.8测站点、照准点归心改正的计算程序4.9三角高程高差的计算程序4.10测距边水平距离的计算程序4.11单一导线的严密平差程序4.12边、角后方交会点的严密平差程序4.13测边交会点的严密平差程序5.1图跟导线点座标的计算程序5.2 碎部点座标的计算程序5.3 内、外分点座标的计算程序5.4 角座标的计算程序5.5 圆心座标的和半径的计算程序5.6 交点座标的计算程序5.7 单边、角交会点座标的计算程序5.8 宗地有关项目的计算程序5.9 单边、双站后方交会座标的计算程序6.1 线路边桩座标的计算程序6.2 红线交点座标的计算程序6.3 直线与曲线交点座标的计算程序6.4 道路不等宽时曲线元素的计算程序6.5 垂足座标的计算程序6.6 圆曲线敷设的计算程序6.7 缓各曲线敷设的计算程序6.8 极座标放线的计算程序6.9 单一水准路线的近似平差计算程序6.10丘地有关项目的计算程序 公路中、边桩坐标计算及放样程序 功 能：任意半径曲线中桩、任意交角边桩坐标计算放样、中桩切线方位角、任意点置仪放样、两点间距离计算、两点方位角计算。各子程序可独立执行。适用范围：单圆曲线、基本型缓和曲线、非对称基本型缓和曲线、卵型曲线、回头曲线（转角大于180度）源程序：F 2（2是文件名为） “O=0 X Y Q SO O0 X Y V ”：ORS“LS1”U“LS2”P“ZH-Q”W“Z-G”T“J-G”A“ZH”：S=0=S=E-9：S=S：（或者将S=0=S=E-9：S=S改为S=S+E-9也可以）B=A+S：F“HZ”：E=F-U：G“X-ZH”H“Y-ZH”I“X：HZ”J“Y：HZ”M“XO”N“YO”X“X-HS”Y“Y-HS”：Prog“A”：Prog“D”：Lbi0：Z：Z：D=0=Goto1：Goto2：Lbi2：D：D：Lbi1：ZA=Prog“3”：Prog“0”：Goto0：ZA=ZB=L=Z-A：K=RS：Prog“4”：Prog“5”：Prog“0”：Goto0：ZB=ZE=L=S：K=RS：Prog“4”：Prog“5”：Prog“6”： Prog“0”：Goto0：ZE=ZF=L=F-Z：K=RU：Prog“4”：Prog“9”：Prog“7”：Prog“0”：Goto0：ZF=C=Z-F：Prog“8”：Prog“0”：Goto0F 2（2是文件名）也可用下面这个程序“O=0 X Y Q SO O0 X Y V ”：ORS“LS1”U“LS2”P“ZH-Q”W“Z-G”T“J-G”A“ZH”：B=A+S：F“HZ”：E=F-U：G“X-ZH”H“Y-ZH”I“X：HZ”J“Y：HZ”M“XO”N“YO”X“X-HS”Y“Y-HS”：Prog“A”：Prog“D”：Lbi0：Z：Z：D=0=Goto1：Goto2：Lbi2：D：D：Lbi1：ZA=Prog“3”：Prog“0”：Goto0：ZA=ZB=L=Z-A：K=RS：Prog“4”：Prog“5”：Prog“0”：Goto0：ZB=ZE=L=S：Goto3：ZE=ZF=L=F-Z：K=RU：Prog“4”：Prog“9”：Prog“7”：Prog“0”：Goto0：ZF=C=Z-F：Prog“8”：Prog“0”：Goto0Lbi3：S=0=K=1：Prog“4”：Prog“5”：Prog“6”：Prog“0”：Goto0： S0=K=RS：Prog“4”：Prog“5”：Prog“6”：Prog“0”：Goto0F 3C=A-Z：X=G+Ccos（P+180）：Y=H+Csin（P+180）：V=PF 4Lbi1：X=L-L540K2+L93456K4-L13599040K6+8-L217.80337152E10K10：Y=L36K-L7336K3+L1142240K5-L159676800K7+L193530096640K9-L231.880240947E12K11：C=（X2+Y2）：V=90L2K：X0=Q=tan-1YX：Q=0：F 5Prog“9”：X=G+Ccos（P+Q）：Y=H+Csin（P+Q）： V=P+VF 6L=Z-B：Q=90LR：C=2RsinQ：W0=Q= Q：V= V：Goto1： V= V：Q= -Q：Goto1：Lbi1：X=X+Ccos（V+Q）：Y=Y+Csin（V+Q）：V=V+2QF 7X=I+Ccos（P+W-180-Q）：Y=J+Csin（P+W-180-Q）：V=P+W-VF 8X=I+Ccos（P+W）：Y=J+Csin（P+W）：V=P+WF 9W0=Q：V：Q= -Q：V= -V：F 0X=X+Dcos（V+T）：Y=Y+Dsin（V+T）：O=0=“X=”：X：Pause0：Y=YProg“A”：Prog“D”：O0=“X=”：X：Pause0：Y=YV0=V= V+360 V= VF AXM=Prog“B”：Goto0：X=M=Prog“C”：XM=Prog“B”：Q=Q+180：Lbi0：Q=Q+360：Q360=Q=Q-360Q=QF BQ=tan-1（Y-N）（X-M）F C YN=Q=90Q=270F D C=（X-M）2+（Y-N）2）：“SO=”：Pause0说明：一、程序中字母及符号意义：R：圆曲线半径；LS1（S）：第一缓和曲线长LS2（U）：第二缓和曲线长ZH-Q（P）：直缓（直圆）点切线方位角Z-G（W）弯道转角（左转为负，右转为正）J-G（T）：中桩至右侧某点方向与中桩切线方位角的夹角（大于等于0度且小于等于180度，当正交时为90度）ZH（A）：直缓或直圆点桩号HZ（F）：缓直或圆直点桩号X-ZH（G）：直缓或直圆点X坐标；Y-ZH（H）：直缓或直圆点Y坐标X-HZ（I）：缓直或圆直点X坐标Y-HZ（J）：缓直或圆直点Y坐标XO（M）：测站X坐标； YO（N）：测站X坐标 X-HS：后视点X坐标； Y-HS：后视点Y坐标X、Y：计算或放样点坐标Q（Q）：计算或放样方位角 SO（C）：计算或放样距离Z（Z）：计算点桩号D（D）：边桩距中桩宽度（左为负值，右为正值）V（V）：中桩切线方位角：分数符号Prog：子程序转移命令Pause：显示暂停命令tg-1：反正切函数Goto：无条件转移命令Lbl：对应无条件转移，指定行号。二、输入、计算要点1该程序一次只可输入一个弯道的参数，计算段落为上一弯道终点（HZ或YZ）至下一弯道起点（ZH或ZY）。2计算单圆曲线时LS1、LS2输入时输03当只计算第一缓和曲线及圆曲线，不计算第二缓和曲线时，弯道转角只需输入正或负值（左转为正、右转为负）即可，可不输入准确的角度。4如只计算上一弯道曲线（HZ或YZ）终点至ZH点间的桩号时，除ZH点桩号、方位角和坐标须输入正确外，其它参数可任意输入（但不能有逻辑错误，如HZ点桩号小于了ZH点的桩号的错误情况。）5后视点X、Y坐标在每次重启程序时需重新输入，建议第一次输入计算出后视角度及距离后用笔将计算的后视角度和距离在纸上记录下来，便于中途后视时核对。6用于卵型曲线时坐标计算时，需先计算出卵型曲线的相关参数后再用本程序计算即可。7为保证计算精度，本程序所用切线支距公式（F4子程序中）取至6项，计算非卵型曲线时为节省存储空间可只取前4项，建议至少取3项。8程序启动时，当问O？时，如输入0则计算结果显示坐标X、Y及放样角度Q和放样距离SO；如输入不为0的数时显示坐标X、Y及中桩切线方位角。9当问距中桩宽度D？时，如只放中桩可输入0，程序默认为中桩后以后不再提示输入宽度；如要同时计算中、边桩，在计算中桩时可输入0.0001或E- 4等更小的数即可。10如已知两点坐标只计算其方位角或距离时，直接运行F A或F D子程序即可。11如需将多个曲线参数一次置入程序，减少每次输入的繁锁，可对F 2程序进行修改或另加一主程序即可，其它子程序不用改动。公路路线坐标反算适用范围：单圆曲线、基本型缓和曲线、非对称基本型缓和曲线、卵型曲线、回头曲线（转角大于180度）特点：1已知任意点坐标求对应于中桩任意方向其对应中桩的桩号及到中桩距离，斜交或正交。2无需人为判断即可求得其位于路线的关系。3本程序公用了前面公路中、边桩坐标计算及放样程序中F 3F9子程序,因此节约了空间。源程序一与源程序二任取一个均可源程序一：F E ZUO BIAO FAN SUANZ=0：O=0：D=0：RS“LS1”U“LS2”P“ZH-Q”W“Z-G”T“J-G”A“ZH”：B=A+S：F“HZ”：E=F-U：G“X-ZH”H“Y-ZH”I“X：HZ”J“Y：HZ”M“XO”N“YO”：D：D：Lbi0：Prog“F”： Goto0F FZ=Z+O2：Lbi1：ZA= Prog“3”：Prog“G”： Goto1：ZA=ZB= L=Z-A：K=RS：Prog“4”：Prog“5”：Prog“G”：ZB=ZE= L=S：K=RS：Prog“4”：Prog“5”：Prog“6”：Prog“G”：ZE=ZF= L=F-Z：K=RU：Prog“4”：Prog“9”：Prog“7”：Prog“G”：ZF= C=Z-F：Prog“8”：Prog“G” F G C=（X-M）2+（Y-N）2）：D0=D= -C：Goto0：D0=D=C：Goto0：Lbi0：X=X+Dcos（V+T）：“X=”：Pause0：Y=Y+Dsin（V+T）：“Y=”：Pause0：O=（X-M）2+（Y-N）2）：OE9= D=D：Goto1：OE9=Z=0：O=0：D= -D：Goto1：Lbi1： Abs（X-M）E-4= Abs（Y-N）E-4= X ：“X=”：Pause0： Y ：“Y=” Z ：“Z=”： Pause0： D ：“D=” MND：M“XO”N“YO”D“D”：Z=0： O=0：Abs（X-M）E-4= Abs（Y-N）E-4= Z源程序二：F E ZUO BIAO FAN SUANZ=0：O=0：D=0：RS“LS1”U“LS2”P“ZH-Q”W“Z-G”T“J-G”A“ZH”：B=A+S：F“HZ”：E=F-U：G“X-ZH”H“Y-ZH”I“X：HZ”J“Y：HZ”M“XO”N“YO”：D：D：Lbi0：Z=Z+O2：ZA= Prog“3”：Prog“F”： Goto0：ZA=ZB= L=Z-A：K=RS：Prog“4”：Prog“5”：Prog“F”：Goto0：ZB=ZE= L=S：K=RS：Prog“4”：Prog“5”：Prog“6”：Prog“F”：Goto0：ZE=ZF= L=F-Z：K=RU：Prog“4”：Prog“9”：Prog“7”：Prog“F”：Goto0：ZF= C=Z-F：Prog“8”：Prog“F” ：Goto0F F AbsXE9=X=0：Y=0：Z=0：O=0：D= -D： Goto0：Lbi0： C=（X-M）2+（Y-N）2）：D0=D= -C：Goto1： D=C： Lbi1：X=X+Dcos（V+T）：“X=”：Pause0：Y=Y+Dsin（V+T）：“Y=”：Pause0：O=（X-M）2+（Y-N）2）：Abs（X-M）E-4= Abs（Y-N）E-4= X ：“X=”：Pause0： Y ：“Y=” Z ：“Z=”： Pause0： D ：“D=” MND：M“XO”N“YO”D“D”：Z=0说明：一、程序中字母及符号意义：D：待求点距中桩宽度（左负右正）XO（M）：待求点已知X坐标YO（N）：待求点已知Y坐标Z：待求点对应中桩桩号其它与公路中、边桩坐标计算及放样程序中符号意义相同。二、输入、计算要点1D值的输入：如能判断其位于路线的左或右侧（其值为左负右正），则会缩短程序运行时间；也可不输入由程序自行计算。2参数的输入：如能确定所求点坐标所处的曲线位置时则输入该曲线参数，也会缩短程序运行时间。3假如事先不知道已知点位于哪一曲线上又怎样处理呢？输入你认为已知点可能位于的那个曲线参数，然后看计算结果的桩号是否处于你输入的那个曲线，如在则结果可用，如不在则输入计算结果显示桩号所处的那个曲线参数即可。4程序计算结果精度设定小于0.1mm，为缩短运行时间可将精度重设，方法是在子程序F G中将有E-4的改为E-3（精度设定小于1mm）。5当已知点所处的曲线半径较小时，如R为20、30m时，程序可能无法得出正确结果或出现死循环时，将程序中止进入程序F F中将表达式Z=Z+O2改为Z=Z+O4即可，改后的程序执行所需时间较修改前长，建议一般不改，除非出现上述情况时再作修改。6此程序必须与公路中、边桩坐标计算及放样程序中F 3F 9子程序结合使用，不能独立运行。已知一点A求其对应于中桩的桩号及到中桩的宽席计算原理（试算渐近法）：此分析是基于用于编程而考虑的1. 先假定从路线起点K0+000处为计算起点(如果不是编程序可不从0桩号开始，可根据经验或实际进行估计)，求得0桩号的中桩坐标，再求得其与已知点的距离C1，利用中桩算得C1宽度的边桩坐标点B1，再算B1点与已知点A的距离D1。2. 算得Z1的桩号（为Z0的桩号加上前一个桩号对应边桩B1与已知点A的距离D1的一半，见图），求出其中桩坐标，再算Z1中桩与已知点A的距离C2，再算其C2宽度的边桩B2，计算B2与已知点A的距离D2。3. 算Z2的桩号，同上依此类推。4. 由上我们不难看出，经过多次计算后已知点对应的中桩与对应宽度越来越接近真值。只要中桩Zi对应Ci宽度的边桩坐标与已知点坐标进行比较在规定精度内即认为是正确的可不再计算。5. 对于不同交角，只是在计算边桩坐标时代入不同交角即可。6. 对半径较小的圆曲线或缓和曲线，在计算桩号时可将Di/2中的分母中的2改为3或4，但应该说3就可以了。竖曲线设计标高计算源程序：A“P1”B“P2”RD“JD ZHANG”I“JD-H”：LBI0：Z：Z：W=B-A：T=AbsWR2：E=TW4：C=D-T：F=D+T：W0=K=（Z-