5800公路测量程序(修改)

..CASIO-5800P计算器以其轻巧灵活、便于携带、功能强大、经济实用的特点,深受广大测量爱好者的喜爱。尤其在公路施工测量方面,全站仪与CASIO-5800P计算器的组合是目前公路施工测量应用最广泛的工具。为使其功能得到更好的应用,特编辑整理此程序,供广大公路测量爱好者参考。不妥之处敬请指教！编者：管东利2009-3-15简要说明本程序具有以下特点：简洁实用,结构合理,功能齐全。由项目程序、主程序和子程序构成本程序系统。有功能强大的路线坐标计算程序,同时配备了其它常用的应用程序。如：后方交会、面积计算、点到直线的距离及两直线交点坐标、距离方位反算等。路线坐标计算：可以将整个项目的多条路线的多条曲线要素一次编入子程序<P>；竖曲线要素编入子程序<G>,横坡度编入子程序〔I。做到一次编辑,项目程序<如：YL2>多次运行,直到本项目结束,避免每次计算都必须输入曲线要素的重复劳动,也减少了计算出错机会,极大地提高工作效率。这一点也是本程序的精华所在。可对每条曲线进行坐标计算：即可使用曲线要素法计算,也可使用主点桩坐标法进行计算。其中主点桩坐标法可以用来计算任何曲线,尤其适用于对付互通立交。与设计文件紧密结合。几乎所有已知曲线要素都可以在设计文件中直接找到,避免了重新计算的麻烦。变量名力求与设计、习惯或数学公式相匹配,易于联想记忆。如：半径-R、直缓点-ZH、水平角-HR、里程桩号-LX等。功能齐全的子程序,使编制新的应用程序简便易行。CASIO-5800P公路施工测量应用程序<GDL-2009>管东利YL2<程序名：永蓝2段>Lb11:Prog"1":Prog"P":Prog"2":Goto1输入：LX:计算点桩号〔对应变量V;S:计算点距该中桩距离〔对应变量S;W:计算点与中桩连线与路线前进方向的夹角〔对应变量W计算并显示：计算点坐标X,Y说明:本程序为项目程序,先将本项目的平曲线要素编辑到子程序P,一次输入,长期使用。YLGC<永蓝高程>Lb11:Prog"1":Prog"G":Prog"I":"H=":H+I<S-1→H◢Goto1输入：LX:计算点桩号〔对应变量V;S:计算点距该中桩距离〔对应变量S;W:计算点与中桩连线与路线前进方向的夹角〔对应变量W计算并显示：计算点对应的路面设计高程H说明:本程序为项目程序,先将本项目的竖曲线要素分别编辑到子程序G和子程序I中,一次输入,长期使用。<注：S要大于或等于零；W为±900；1为高程设计线到路中线的距离,一般为中央分隔带宽度的一半,要与设计图纸相符。>ZX<直线段>Prog"8":Lb11:Prog"1":Prog"Z":Prog"2":Goto1已知：LA:第一点桩号〔对应变量K;X1:第一点X坐标〔对应变量A;Y1:第一点Y坐标〔对应变量B;LN:第二点桩号〔对应变量N;X2:第二点X坐标〔对应变量C;Y2:第二点Y坐标〔对应变量D输入：LX:计算点桩号〔对应变量V;S:计算点距该中桩距离〔对应变量S;W:计算点与中桩连线与路线前进方向的夹角〔对应变量W计算并显示：计算点坐标X,Y说明：对应变量,在编辑子程序"P"和"G"时用的到,要对照输入相应的值。YQX<圆曲线>Prog"8":?R:Prog"7":Lb11:Prog"1":Prog"Y":Prog"2":Goto1已知：LA:第一点桩号〔对应变量K;X1:第一点X坐标〔对应变量A;Y1:第一点Y坐标〔对应变量B;LN:第二点桩号〔对应变量N;X2:第二点X坐标〔对应变量C;Y2:第二点Y坐标〔对应变量D；R:圆曲线半径〔对应变量R;Q:曲线转向<Y=1右转Z=-1左转>〔对应变量Q输入：LX:计算点桩号;S:计算点距中桩距离;W:中桩与计算点连线与路线前进方向夹角计算并显示：计算点坐标X,Y示意图如下：HHQX<缓和曲线>"ZH"？K："HZ"？N："XD"？C："YD"?D："XH"?A："YH"?B:"R"?R:"T"?T:"LS"?U:Prog"7":Lb11:Prog"1":Prog"H":Prog"2":Goto1已知：ZH:直缓点桩号〔对应变量K;HZ:缓直点桩号〔对应变量N;XD:交点X坐标〔对应变量C;YD:交点Y坐标〔对应变量D;XH:后方交点X坐标〔对应变量A;YH:后方交点Y坐标〔对应变量B;R:圆曲线半径〔对应变量R;T:切线长〔对应变量T;Ls:缓和曲线长<对应变量U;Q:Y=1右转Z=-1左转〔对应变量Q输入：LX:计算点桩号〔对应变量V;S:计算点距该中桩距离〔对应变量S;W:计算点与中桩连线与路线前进方向的夹角〔对应变量W计算结果：计算点坐标X,Y说明:1.适用于标准的缓和曲线,既第一,第二缓和曲线长度相等。2.计算范围为：含两端直线段在内的整条曲线。FHQX<复合曲线>Prog"8":"R1"?H:"R2"?R:Prog"7":Lb11:Prog"1":Prog"F":Prog"2":Goto1已知：LA:第一点桩号〔对应变量K;X1:第一点X坐标〔对应变量A;Y1:第一点Y坐标〔对应变量B;LN:第二点桩号〔对应变量N;X2:第二点X坐标〔对应变量C;Y2:第二点Y坐标〔对应变量D；Q:曲线转向<Y=1右转Z=-1左转>〔对应变量Q;R1:第一点曲线半径〔桩号LA点>〔对应变量H;R2:第二点曲线半径<桩号LN点>〔对应变量R〔如：缓圆点或圆缓点的曲线半径等于圆曲线半径R；直缓点或缓直点的曲线半径为+∞,可输入1G<千兆>；其他项同HHQX说明：1、LA和LN在缓和曲线段上即可,不必是起点或终点；对应的曲率半径计算要正确。2、缓和曲线曲率半径计算公式：R=A2/L〔A为缓和曲线参数；L为计算点至直缓点之间的缓和曲线长度；3、LA可以大于LNSQX<竖曲线>"K1"?A:"K2"?B:"LD"?C:"HD"?D"R"?R:Lb11:"LX"?V:Pros"S":"H="H→H◢Goto1已知：K1:前一段路线纵坡〔对应变量AK2:后一段路线纵坡〔对应变量BLD:竖曲线交点桩号〔对应变量CHD：交点高程〔对应变量DR:竖曲线半径〔对应变量R输入：LX计算点桩号计算:LX桩号的路线设计高程HLXSW<桩号反算>Lb10:"LC"？V:"X"？X："Y"？Y:X→Z[1]:Y→Z[2]:DO:Prog"P":Z[1]→E:Z[2]→G:0=<X-A>2+<Y-B>2=>E+.0001→E:Pol<E-X,G-Y:Icos<J-M>→K:V+K→V:Lpwhile.005≤Abs<K:"LX=":V◢"S=":I→S◢J-M→W:180＜Abs<W=>W-360Int<W÷360+.5W÷Abs<W→W:"W=":WDMS◢Z[1]→X:Z[2]→Y:Goto0已知：猜测桩号LC,已知点P坐标〔X,Y。计算：已知点P〔X,Y相对于路线的位置。即对应的桩号LX,距中桩距离S,位于路线的哪一侧〔W=90在右侧；W=-90在左侧。说明：1、本程序在路线曲线要素已编入子程序"P",即项目程序正常运行时方可运行。2、输入程序时注意区别：字母O和数字0、G和1G<千兆>、K和1k<千>、M和1M<兆>。HDZB<涵洞坐标>Lb11:"LX"？V:"W"？W："LZ"？E："LY"？G:"B/2"？B：Sin<W>=0=>Goto1:V→Z[1]:W→Z[2]:E→Z[3]:G→Z[4]:B→Z[5]:Prog"P":X→A:Y→B:Rec<Z[4],M+W:X+I→C:Y+J→D:0→K:X→A:Z[4]→N:-Z[3]→V:-Z[2]→W:Abs<Z[5]÷Sin<W→S:Prog"Z":Rec<S,M+W:I→E:J→G:"X1=":X+E◢"Y1=":Y+G◢"X4=":X-E◢"Y4=":Y-G◢0→V:Prog"Z":"X2=":X+E◢"Y2=":Y+G◢"X5=":X-E◢"Y5=":Y-G◢Z[4]→V:Prog"Z":"X3=":X+E◢"Y3=":Y+G◢"X6=":X-E◢"Y6=":Y-G◢Z[1]→V:Z[2]→W:Z[3]→E:Z[4]→G:Z[5]→B:Goto1输入：LX：涵中心桩号〔对应变量VW：涵交角〔对应变量WLZ：左涵长〔对应变量ELY：右涵长〔对应变量GB/2：涵正宽的一半〔对应变量B计算:涵放样点坐标〔X1,Y1-<X6,Y6>说明：本程序在路线曲线要素已编入子程序"P",即项目程序正常运行时方可运行。HCJS<涵长计算>Lb11:"LX"？U:"W"？W："B/2"？B："Z"？Z："H+D"？P：B→Z[1]:B→S:Sin<W>=0=>Goto1:Abs<S÷Sin<W→E:U+Ecos<W→V:Prog"G":Prog"I":H+I<S-1-.75>-.03→H:If8＞H-Z-P:Then<S+1.5<H-Z-P>>÷Sin<W→S:Else<S+1.5<H-Z-P>+.25<H-Z-P-8>+2>÷Sin<W→S:IfEnd:"L=":S◢Z[1]→B:Goto1输入：LX：涵中心桩号〔对应变量VW：涵交角〔对应变量WB/2：计算一侧的路面宽度〔对应变量BZ：涵底标高〔对应变量ZH+D：涵底至盖板顶高度〔对应变量P计算:对应W值一侧的涵长L〔对应变量S<计算另一侧涵长,将W输入对应的值,一般为W+180>说明：本程序在竖曲线要素已编入子程序"G",横坡度要素编入子程序"I"时方可运行,内容也要根据边坡设计图纸做适当修改。TTBX<填土边线>"m1"?E:"m2"?F:"m3"?G:"B1"?O:"B2"?P:E→Z[3]:F→Z[4]:G→Z[5]:O→Z[6]:P→Z[7]:Lb11:"LX"？V:"W"？W："B/2"？B："Z"？Z：B→Z[1]:Z→Z[2]:B→S:1≠Abs<Sin<W=>Goto1:Prog"G":Prog"I":H+I<S-1-.75>-.03→H:H≤Z=>Goto1:IfH≤Z+8:ThenS+Z[3]<H-Z>→S:Goto2:IfEnd:IfH≤Z+20:ThenS+8Z[3]+Z[4]<H-Z-8>+O→S:Goto2:IfEnd:IfH≤Z+32:ThenS+8Z[3]+12Z[4]+Z[5]<H-Z-20>+O+Z[7]→S:Goto2:IfEnd:H＞Z+32=>Goto1:Lb12:"TH=":H-Z◢"S=":S◢Prog"P":Prog"2":Z[1]→B:Z[2]→Z:Z[3]→E:Z[4]→F:Z[5]→G:Z[6]→O:Z[7]→P:Goto1输入：一级、二级、三级坡比分别为m1、m2、m3；一级、二级平台宽度分别为B1、B2LX：路中心桩号〔对应变量VW：在路的哪一侧,交角±90〔对应变量WB/2：半幅路面宽度〔对应变量BZ：实测地面标高〔对应变量Z计算:计算显示填土高度"TH=";填土边线到中桩距离"S=";对应W值一侧填土边线坐标<X,Y>说明：本程序在在路线曲线要素已编入子程序"P",竖曲线要素已编入子程序"G",横坡度要素编入子程序"I"时方可运行,内容也要根据边坡设计图纸做适当修改。<本次边坡设计见下图,距中1m为高程设计线,土路肩宽度0.75m.暂时编辑三级坡以下的。>HJ<后方交汇>Prog"3":"H"?H:"IH"?K:"RH"?R:"HR1"?O:"HD1"?L:"HR2"?P:"HD2"?S:"VD"?T:Prog"6":I→Q:P-O→W:Rec<S,W:Pol<I-L,J:Q÷I→Q:Rec<QS,W:0→U:0→V:QL→E:0→G:I→M:J→N:Prog"B":"Z=":H-T-K+R→Z已知：两已知点坐标〔X1,Y1;〔X2,Y2,H：任意一点高程；实测：IH：仪器高；RH：棱镜高；HR1：第一点水平角读数；HD1：第一点水平距离读数；HR2：第二点水平角读数；HD2：第二点水平距离读数；VD：垂直高差读数<对应H>计算并显示：所求测站点坐标<X,Y,Z>FWFS<方位反算>"X0"?A:"Y0"?B:Lb11:"X"?X:"Y"?Y:Pol<X-A,Y-B:"L=":I→L◢"HR=":J+180<1-AbsJ÷J◢DMS◢Goto1已知：两点坐标P0<X0,Y0>；P<X,Y>计算:P0点到P点的距离L和方向角HR〔显示度、分、秒。MJ<面积>Prog"3":Prog"6":0→S:I→L:0→E:Lb11:"X"？X："Y"？Y:.5<AD+BX+CY-DX-AY-CB→K:"S=":S+K→S◢√<<X-C>2+<Y-D>2→G:√<<X-A>2+<Y-B>2→H:"L=":L-E+G+H→L◢X→C:Y→D:H→E:Goto1已知:多边形各角点坐标,按顺时针或按逆时针顺序输入各角点坐标,计算闭合多边形面积S〔为累计值及周长L。示意图如下:JD<交点>Prog"3":Prog"6":Lb11:"X3"？E："Y3"？G:Pol<E-A,G-B:"S="：Isin<J-F→S◢"X4"？M："Y4"？N：Pol<M-E,N-G:J-F→W:IfInt<W÷180>≠W÷180:Then"XD=":<CEB-CMB+AMD-AED+CGM-AGM+AEN-CEN>÷<EB-MB+MD-ED+CG-AG+AN-CN→X◢"YD=":<DGA-DNA+BNC-BGC+DNE-BNE+BMG-DMG>÷<GA-NA+NC-GC+DE-BE+BM-DM→Y◢IfEnd:Goto1已知：第一条直线上两点坐标分别为P1〔X1,Y1、P2〔X2,Y2;第二条直线上两点坐标分别为P3〔X3,Y3、P4〔X4,Y4。计算并显示：P3〔X3,Y3到直线〔P1P2的垂直距离S;以及两直线交点坐标<XD,YD>。示意图如下：子程序1"Lx"？V："S"？S："W"？W：Return2IfS≥0:ThenRec<S,M+W:ElseRec<-S,M+W-180:IfEnd:"X=":X+I→X◢"Y=":Y+J→Y◢：Return3"X1"?A:"Y1"?B:"X2"?C:"Y2"?D：Return40→Z:K>N=>180→Z:cos<Z→P：Return5A+I→X:B+J→Y：Return6Pol<C-A,D-B:J→F：Return7"Y=1,Z=-1"?Q：Return8"LA"?K:"X1"?A:"Y1"?B:"LN"?N:"X2"?C:"Y2"?D:Return9N→K:C→A:B→D:ReturnXRU÷L→O:L-L^<3>÷<40O2>+L^<5>÷<3456O^<4>>-L^<7>÷<599040O^<6→X:L2÷<6O>-L^<4>÷<336O^<3>>+L^<6>÷<42240O^<5>>-L^<8>÷<9676800O^<7→Y：ReturnHU=0=>1μ→U:Prog"4":Prog"6":P<V-K→L:0→H:If<PV>>.5P<K+N:ThenP<N-V→L:180→H:F+180Q<N-K-PU>÷<πR>+H→F:IfEnd:Rec<T,F:C-I→E:D-J→G:IfL≦0:ThenRec<-L,F+180:F+Z+H→M:1M→O:Goto2:IfEnd:IfL≦U:ThenProg"X":Pol<X,Y:Rec<I,F+PQJcos<H:F+90PQLcos<H>÷<πO>+Z+H→M:Goto2:IfEnd:45<2L-U>÷<πR→J:45U÷<πR→I:F+2PQJcos<H>+Z+H→M:U→L:Prog"X":Rec<2Rsin<J-I,J+I:X+I→X:Y+J→Y:Pol<X,Y:Rec<I,F+PQJcos<H:Lb12:E+I→X:G+J→Y：ReturnS1→E:A>B=>-1→E:.5RAbs<A-B→T:A→M:V>C=>B→M:IfT≤Abs<V-C:ThenD+M<V-C→H:ElseV-C+T→G:D-AT+Gtan<tan-1<A>+90EG÷<πR→H：IfEnd:ReturnFR<N-K>÷<H-R→M:N-K+M→U:M→L:Prog"X":X→E:Y→G:U→L:Prog"X":Prog"6":I→Z:Pol<X-E,Y-G:Z÷I→Z:F-QJ→F:V=K=>K+1μ→V:V-K+M→L:Prog"X":Pol<X-E,Y-G:Rec<ZI,F+QJ:F+90QL÷<πO→M：Prog"5"：ReturnZProg"4":Prog"6":F+Z→M:IfV≥K:ThenRec<V-K,M:ElseRec<K-V,M+180:IfEnd:Prog"5":1G→O:1→Q：ReturnYProg"6":Prog"4":F-90Q<N-K>÷<πR>+Z→E:90Q<V-K>÷<πR→G:IfQRsin<G>≥0:ThenRec<2QRsin<G,E+G:ElseRec<-2QRsin<G,E+G-180:IfEnd:E+2G→M:Prog"5":R→O：ReturnBProg"6":I→Q:Pol<M-E,N-G:Q÷I→Q:F-J→F:Pol<U-E,V-G:Rec<QI,F+J:"X=":A+I→X◢"Y=":B+J→Y◢：Return匝道"1-4：WA-WB：5-8：EA-ED：9WN：10WS：11EN：12ES"?T：Lb10：Porg"1":IfT=1:ThenPorg"P"：Porg"2"：Goto0:IfEnd:IfT=2:ThenPorg"WA"：Porg"2"：Goto0:IfEnd:IfT=3:Porg"WB"：Porg"2"：Goto0:IfEndP〔示例：永蓝主线K99+000-K106+012第二缓和曲线终点桩号→N:IfV≤N:Then第二缓和曲线起点桩号→K:起点x坐标→A:起点y