5800公路測量程序(修改)

..CASIO-5800P計算器以其輕巧靈活、便于攜帶、功能強大、經濟實用的特點,深受廣大測量愛好者的喜愛。尤其在公路施工測量方面,全站儀與CASIO-5800P計算器的組合是目前公路施工測量應用最廣泛的工具。為使其功能得到更好的應用,特編輯整理此程序,供廣大公路測量愛好者參考。不妥之處敬請指教！編者：管東利2009-3-15簡要說明本程序具有以下特點：簡潔實用,結構合理,功能齊全。由項目程序、主程序和子程序構成本程序系統(tǒng)。有功能強大的路線坐標計算程序,同時配備了其它常用的應用程序。如：后方交會、面積計算、點到直線的距離及兩直線交點坐標、距離方位反算等。路線坐標計算：可以將整個項目的多條路線的多條曲線要素一次編入子程序<P>；豎曲線要素編入子程序<G>,橫坡度編入子程序〔I。做到一次編輯,項目程序<如：YL2>多次運行,直到本項目結束,避免每次計算都必須輸入曲線要素的重復勞動,也減少了計算出錯機會,極大地提高工作效率。這一點也是本程序的精華所在。可對每條曲線進行坐標計算：即可使用曲線要素法計算,也可使用主點樁坐標法進行計算。其中主點樁坐標法可以用來計算任何曲線,尤其適用于對付互通立交。與設計文件緊密結合。幾乎所有已知曲線要素都可以在設計文件中直接找到,避免了重新計算的麻煩。變量名力求與設計、習慣或數學公式相匹配,易于聯(lián)想記憶。如：半徑-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、輸入程序時注意區(qū)別：字母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