CASIO5800測(cè)量程序

1、CASIO fx-4800P型計(jì)算器用于線路施工曲線中線點(diǎn)坐標(biāo)的計(jì)算程序中鐵十局三建公司工程技術(shù)部 摘要：本文介紹了CASIO fx-4800P型計(jì)算器程序編制用于鐵路、公路曲線線路內(nèi)任意中線點(diǎn)的坐標(biāo)計(jì)算程序及使用方法。本計(jì)算程序具有操作簡(jiǎn)便、計(jì)算快捷、應(yīng)用廣泛等特點(diǎn)、極大地減輕了測(cè)量工作者的內(nèi)業(yè)工作量，對(duì)于測(cè)量工作者有較大的參考和指導(dǎo)作用。關(guān)鍵詞：曲線 線路 施工測(cè)量 計(jì)算程序1概述 過(guò)去，線路中線施工放樣基本依靠經(jīng)緯儀和鋼尺了來(lái)進(jìn)行角度及距離測(cè)量。對(duì)于曲線線路一般的測(cè)量方法是：經(jīng)緯儀置于某一中線點(diǎn)上，采用偏角法撥角再用鋼尺量距來(lái)定出中線點(diǎn)。隨著電子技術(shù)進(jìn)步和經(jīng)濟(jì)發(fā)展，測(cè)量?jī)x器和測(cè)量方法的不

2、斷改進(jìn)，目前，全站儀已廣泛地應(yīng)用于工程施工測(cè)量中，極大的提高了測(cè)量工作效率。但是，在進(jìn)行鐵路、公路工程的曲線線路施工測(cè)設(shè)時(shí)，需要在線路所在區(qū)域建立統(tǒng)一坐標(biāo)系或獨(dú)立坐標(biāo)系，利用坐標(biāo)變換的方法，將整個(gè)曲線的三個(gè)部分（第一緩和曲線、中間圓曲線、第二緩和曲線）統(tǒng)一到同一坐標(biāo)系中。根據(jù)坐標(biāo)系的建立，計(jì)算出整個(gè)曲線內(nèi)任意點(diǎn)的坐標(biāo)，再采用全站儀利用極坐標(biāo)方法進(jìn)行施工放樣。前提是首先利用計(jì)算器計(jì)算出各中線點(diǎn)坐標(biāo)，然后才能進(jìn)行放樣。而普通型計(jì)算器不僅計(jì)算速度慢，且要求計(jì)算者必須正確地記憶很多計(jì)算公式，計(jì)算繁瑣而且容易出錯(cuò)，滿足不了現(xiàn)場(chǎng)測(cè)設(shè)工作的要求。為了能夠快速準(zhǔn)確地為全站儀提供測(cè)設(shè)數(shù)據(jù)，發(fā)揮全站儀快速測(cè)設(shè)的特

3、點(diǎn)，提高測(cè)量工作效率，應(yīng)采用可編程的計(jì)算器，編制計(jì)算程序。本文主要介紹應(yīng)用CASIO fx-4800P型計(jì)算器的計(jì)算程序，供公司測(cè)量同行們參照使用。2計(jì)算程序QXZBJS（文件名：曲線坐標(biāo)計(jì)算fx-4800P） Defm2:R：L：A：N“ZH：X=”：E“ZH：Y=”：F：“FWJ=”：K“ZH：LC=”： P=L2/(24R)-L4/(2688R3):M=L/2-L3/(240R2):T“T”=（R+P）tng(A/2)+MG=RA/180：“S”S=G+LLbiA：C，V：C“CSDLC=”：V“HXPJ=”：D=C-K：DLI=D-D5/（40R2L2）：U=D3/(6RL)-D7/（

4、336R3L3):J=（I2+U2）：Goto1：DGO=90（2D-L）/( R)：I=RsinO+M：U=R（1-cosO）+P：J=（I2+U2）：Goto2：D=S-（C-K）： = D-D5/（40R2L2）：Z2=D3/6RL-D7/（336R3L3）：I=T+（T-Z1）cosA-Z2sinA：U=（T-Z1）sinA+Z2cosA：J=（I2+U2）：Goto3：Lbi1:Q:Q“Z=1；Y=2”：Q=1Q=F-30D2/ （RL）：H=F-90D2/ （RL）：Q=F+30D2/ （RL）：H=F+90D2/（ RL） Goto4：Lbi2:Q:Q“Z=1；Y=2”：Q=1Q

5、=F-tng-1(U/I)：H=F-O：Q= F+tng-1(U/I)：H=F+OGoto4：Lbi3:Q:Q“Z=1；Y=2”：Q=1Q=F-tng-1(U/I):H=F-(A-90(S-(C-K)2/ （RL）): Q= F+tng-1(U/I):H=F+(A-90(S-(C-K)2/ （RL）): Goto4： Lbi4:B=90+H:H0H“QXFWJ”=B+360 H360H“QXFWJ”=H-360 H“QXFWJ”=HGoto5：Lbi5: X“CSD:X”=JcosQ+N+VcosB Y“CSD:Y”=JsinQ+E+VsinBGotoA 3程序說(shuō)明3.1 輸入已知變量 R圓曲

6、線半徑，顯示R？ L緩和曲線長(zhǎng)，顯示L？ A曲線轉(zhuǎn)向角，顯示A？ E直緩點(diǎn)縱坐標(biāo)，顯示ZH：X=？ N直緩點(diǎn)橫坐標(biāo)，顯示ZH：Y=？ F第一切線方位角，即ZH至JD的方位角，顯示FWJ=？ K直緩點(diǎn)里程，顯示ZH：LC=？ 3.2 計(jì)算待求量 T切線長(zhǎng)度，顯示T= S曲線全長(zhǎng)，顯示S= Z3外矢距，即JD到QZ的距離，顯示E0= 3.3 輸入待求變量 K輸入待求（測(cè)設(shè)）點(diǎn)的里程，顯示LC=？ V橫向偏距，即測(cè)設(shè)點(diǎn)左、右側(cè)外移距的偏移量，若為中線點(diǎn)輸入0；右側(cè)輸入“+”值，左側(cè)輸入“-”值。3.4 計(jì)算結(jié)果和說(shuō)明X待求（測(cè)設(shè)）點(diǎn)的縱坐標(biāo)，顯示X= Y待求（測(cè)設(shè)）點(diǎn)的橫坐標(biāo)，顯示Y= 該程序既可對(duì)

7、曲線內(nèi)中線任意點(diǎn)的坐標(biāo)進(jìn)行計(jì)算，又可以對(duì)中線外移樁的坐標(biāo)進(jìn)行計(jì)算。計(jì)算原理如下：先建立整個(gè)曲線的坐標(biāo)系統(tǒng)，確定曲線起點(diǎn)（ZH）的坐標(biāo)，以ZH到JD的第一切線方位角為計(jì)算方位角。 在程序計(jì)算過(guò)程中，當(dāng)屏幕顯示Z=1；Y=2？的提示時(shí)可根據(jù)曲線的轉(zhuǎn)向，從下列兩種情況中選擇一種予以回答。若曲線為左轉(zhuǎn)，則輸入數(shù)值1；若曲線為右轉(zhuǎn)，則輸入數(shù)值2。 本程序式輸入時(shí)，凡遇到、兩種符號(hào)時(shí)，不應(yīng)將其分開，應(yīng)視為一個(gè)整體的符號(hào)。 由于變量的增加，程序式中使用到Z1、Z2、Z3三個(gè)計(jì)算量，此時(shí)應(yīng)將計(jì)算器的內(nèi)存變量擴(kuò)大，另行增加36個(gè)變量即可。在計(jì)算器上的操作過(guò)程為：按SHIFT鍵再按Defm鍵，然后按數(shù)字36即可。

8、3.5 本程序在實(shí)際測(cè)量工作中的應(yīng)用 本程序不僅可以計(jì)算曲線線路上任意中線點(diǎn)的坐標(biāo)，亦可以計(jì)算線路左、右側(cè)的法線偏距點(diǎn)的坐標(biāo)，用于曲線中樁和邊樁放樣。 可以計(jì)算曲線線路中心點(diǎn)的切線方位角和法線方位角，該方位角可以用于橋梁墩、臺(tái)（亦為曲線上橋梁墩、臺(tái)截面縱、橫軸線的方位角）定位。其切線方位角在計(jì)算器上顯示為“QXFWJ= ”。 本程序具有多方面的計(jì)算功能，可根據(jù)儀器狀況和測(cè)設(shè)對(duì)象以及工作需要而靈活應(yīng)用，可在測(cè)設(shè)現(xiàn)場(chǎng)計(jì)算，減輕內(nèi)業(yè)的計(jì)算時(shí)間，提高工作效率。4計(jì)算舉例 結(jié)合某線路的曲線資料，計(jì)算曲線線路上中樁或邊樁的坐標(biāo)。4.1已知要素 R=1500， L=150， AZ=30° 主點(diǎn)里程

9、：ZHDK8+150.000 HYDK8+300.000 YHDK8+935.398 HZDK9+085.398 4.2 建立平面控制坐標(biāo)系 建立獨(dú)立坐標(biāo)系統(tǒng)，若以ZH為坐標(biāo)原點(diǎn)，以ZHJD的切線為X軸的正方向，則ZH點(diǎn)的坐標(biāo)可設(shè)為：X=8150.000，Y=1000.000；切線的初始方位角F“FWJ= ”可設(shè)為0°0000?，F(xiàn)以計(jì)算圓曲線內(nèi)DK8+400的中線點(diǎn)或左側(cè)15米邊樁點(diǎn)的坐標(biāo)為例。4.3 計(jì)算操作步驟 1輸入曲線半徑R=1500按EXE鍵 2輸入緩和曲線長(zhǎng)L=150按EXE鍵 3輸入曲線轉(zhuǎn)向角A=30°按EXE鍵 4輸入ZH點(diǎn)X坐標(biāo)ZH：X=8150.000按

10、EXE鍵 5. 輸入ZH點(diǎn)Y坐標(biāo)ZH：Y=1000.000按EXE鍵 6輸入切線方位角FWJ=0°0000按EXE鍵 7輸入ZH點(diǎn)里程ZH：LC=8150按EXE鍵 顯示T=477.0849919（切線長(zhǎng)）按EXE鍵；顯示S=935.398164（曲線長(zhǎng)）；按EXE鍵顯示E0=53.56126046（外矢距）；按EXE鍵8輸入測(cè)設(shè)點(diǎn)里程LC=8400；按EXE鍵9輸入中線點(diǎn)或左側(cè)邊樁的距離，如計(jì)算中樁點(diǎn)坐標(biāo)輸入PJU=0，如需計(jì)算左側(cè)15米邊樁的坐標(biāo)則輸入PJU=-15米；按EXE鍵10屏幕顯示Z=1；Y=2？提示，本例為左偏曲線故輸入1，按EXE鍵,顯示FXFWJ=83.31549

11、239，按SHIFT鍵再按°，即DK8+400中樁點(diǎn)的法線方位角為83°1855.77按EXE鍵11. 顯示CSD：X=8399.5970289，按EXE鍵； 顯示CSD：Y=989.1782961，按EXE鍵。若需計(jì)算DK8+400左側(cè)15米邊樁點(diǎn)的坐標(biāo)，按EXE鍵，屏幕提示LC=8400；再按EXE鍵，提示PJU=？，輸入-15米，再按EXE鍵，顯示Z=1；Y=2？輸入數(shù)字1按EXE鍵，顯示QXFWJ=83.31549239，按EXE鍵；顯示CSD：X=8397.850997，按EXE鍵；顯示CSD：Y=974.2802637，按EXE鍵。12若需要再計(jì)算其他線路中線點(diǎn)

12、或邊樁點(diǎn)的坐標(biāo)，重復(fù)811步驟操作，不必再輸入各已知數(shù)據(jù)，僅輸入需要計(jì)算點(diǎn)的里程或橫向偏距及Z=1；Y=2？的參數(shù)即可。ZBFS(坐標(biāo)正算)A”X1”:B”Y1”:C”X2”:D”Y2”:Fixm：Pol（C-A,D-B）:I”S=” J<0J=J+360J”FWJ” 測(cè)站點(diǎn)到后視點(diǎn)的距離及方位角X1: Y1為測(cè)站點(diǎn)坐標(biāo)X2: Y2為后視點(diǎn)坐標(biāo)S=為測(cè)站點(diǎn)到后視點(diǎn)的距離FWJ為測(cè)站點(diǎn)到后視點(diǎn)的方位角ZBZS(坐標(biāo)正算)A”X1”:S:F:D”FWJ”: C”X2”=A+SCOS(F) D”Y2”=B+SSIN(F) X1: Y1為測(cè)站點(diǎn)坐標(biāo)X2: Y2為后視點(diǎn)坐標(biāo)S=為測(cè)站點(diǎn)到后視點(diǎn)的距

13、離FWJ為測(cè)站點(diǎn)到后視點(diǎn)的方位角JDZB（文件名：曲線坐標(biāo)計(jì)算CASIO fx-5800P）“JDX1”?A:“JDY1”?B: “JDX2”?C:“JDY2”?D:Pol(C-A,D-B):R：L：A：JF：“ZH：LC=”?K：L2/24R-L4/2688R3P:L/2-L3/240R2M：（R+P）*tng(A/2)+MT: “T=”:TRA/180G： G+LS: “S=”:SC+TCOS(F+180) N: “ZH.X=”:ND+TSIN(F+180) E: “ZH.Y=”:E LbiA：C，V：C“DK?+”：V“PJ”：C-KD：If DL:Then D-D5/（40R2L2)

14、I: D3/6RL-D7/（336R3L3）U:（I2+U2）J：Goto1：IfEnd：IF DG:Then O=90（2D-L）/ （R）： RsinO+MI：R（1-cosO）+PU：（I2+U2）J：Goto2：Else S-（C-K）D:D-D5/（40R2L2）Z1： D3/6RL-D7/（336R3L3）Z2：T+（T-Z1）cosA-Z2sinAI：（T-Z1）sinA+Z2cosAU：（I2+U2）J：Goto3：IfEndLbi1:“Z=1；Y=2”？X：If X=1 Then F-30D2/ （RL）Q：F-90D2/ （RL）H：Else F+30D2/ （RL）Q：F

15、+90D2/（ RL）H:IfEnd: Goto4：Lbi2: X“Z=1；Y=2”：If X=1 Then F-tng-1(U/I) Q： F-OH：Else F+tng-1(U/I) Q：F+OH: IfEnd:Goto4：Lbi3: X“Z=1；Y=2”：If X=1 Then F-tng-1(U/I) Q:F-(A-90(S-(C-K)2/ （RL）) H: Else F+tng-1(U/I) Q: F+(A-90(S-(C-K)2/ （RL）) H: IfEnd: Goto4： Lbi4: 90+HB: If H0 Then B+360H: IfEnd: If H360 Then H

16、-360H : Else HH: IfEnd:Goto5：Lbi5:“X=”:JcosQ+N+VcosBX “Y=”:JsinQ+E+VsinB“FWJ=”:HGotoA 3程序說(shuō)明3.1 輸入已知變量“JDX1” 交點(diǎn)1X坐標(biāo)“JDY1” 交點(diǎn)1Y坐標(biāo) “JDX2” 交點(diǎn)2X坐標(biāo)“JDY2” 交點(diǎn)2Y坐標(biāo) R圓曲線半徑，顯示R？ L緩和曲線長(zhǎng)，顯示L？ A曲線轉(zhuǎn)向角，顯示A？ E直緩點(diǎn)縱坐標(biāo)，顯示ZH：X= N直緩點(diǎn)橫坐標(biāo)，顯示ZH：Y= F第一切線方位角，即ZH至JD的方位角 K直緩點(diǎn)里程，顯示ZH：LC=？ 3.2 計(jì)算待求量 T切線長(zhǎng)度，顯示T= S曲線全長(zhǎng)，顯示S= 3.3 輸入待求

17、變量 K輸入待求（測(cè)設(shè)）點(diǎn)的里程，顯示DK?+？ V橫向偏距，即測(cè)設(shè)點(diǎn)左、右側(cè)外移距的偏移量，若為中線點(diǎn)輸入0；右側(cè)輸入“+”值，左側(cè)輸入“-”值。3.4 計(jì)算結(jié)果和說(shuō)明X待求（測(cè)設(shè)）點(diǎn)的縱坐標(biāo)，顯示X= Y待求（測(cè)設(shè)）點(diǎn)的橫坐標(biāo)，顯示Y= 該程序既可對(duì)曲線內(nèi)中線任意點(diǎn)的坐標(biāo)進(jìn)行計(jì)算，又可以對(duì)中線外移樁的坐標(biāo)進(jìn)行計(jì)算。計(jì)算原理如下：先建立整個(gè)曲線的坐標(biāo)系統(tǒng)，確定曲線起點(diǎn)（ZH）的坐標(biāo)，以ZH到JD的第一切線方位角為計(jì)算方位角。 在程序計(jì)算過(guò)程中，當(dāng)屏幕顯示Z=1；Y=2？的提示時(shí)可根據(jù)曲線的轉(zhuǎn)向，從下列兩種情況中選擇一種予以回答。若曲線為左轉(zhuǎn)，則輸入數(shù)值1；若曲線為右轉(zhuǎn)，則輸入數(shù)值2(非1也對(duì)

18、)。3.5 本程序在實(shí)際測(cè)量工作中的應(yīng)用 本程序不僅可以計(jì)算曲線線路上任意中線點(diǎn)的坐標(biāo)，亦可以計(jì)算線路左、右側(cè)的法線偏距點(diǎn)的坐標(biāo)，用于曲線中樁和邊樁放樣。 可以計(jì)算曲線線路中心點(diǎn)的切線方位角和法線方位角，該方位角可以用于橋梁墩、臺(tái)（亦為曲線上橋梁墩、臺(tái)截面縱、橫軸線的方位角）定位。其切線方位角在計(jì)算器上顯示為“FWJ= ”。 本程序具有多方面的計(jì)算功能，可根據(jù)儀器狀況和測(cè)設(shè)對(duì)象以及工作需要而靈活應(yīng)用，可在測(cè)設(shè)現(xiàn)場(chǎng)計(jì)算，減輕內(nèi)業(yè)的計(jì)算時(shí)間，提高工作效率。4計(jì)算舉例 結(jié)合某線路的曲線資料，計(jì)算曲線線路上中樁或邊樁的坐標(biāo)。4.1已知要素 “JDX1” 交點(diǎn)1X坐標(biāo)“JDY1” 交點(diǎn)1Y坐標(biāo) “JDX

19、2” 交點(diǎn)2X坐標(biāo)“JDY2” 交點(diǎn)2Y坐標(biāo) R圓曲線半徑，顯示R？ L緩和曲線長(zhǎng)，顯示L？ A曲線轉(zhuǎn)向角，顯示A？ E直緩點(diǎn)縱坐標(biāo)，顯示ZH：X= N直緩點(diǎn)橫坐標(biāo)，顯示ZH：Y= F第一切線方位角，即ZH至JD的方位角 K直緩點(diǎn)里程，顯示ZH：LC=？ “JDX1” =2920660.136“JDY1” =475778.751 “JDX2” =2919624.99“JDY2” =477022.978R=4500， L=340， AZ=26°2423”(左轉(zhuǎn)) 主點(diǎn)里程：ZH.LC=DK348+692.554 4.3 計(jì)算操作步驟 1輸入交點(diǎn)1的X坐標(biāo) “JDX1” =2920660

20、.136按EXE鍵 2輸入交點(diǎn)1的Y坐標(biāo) “JDY1” =475778.751按EXE鍵3輸入交點(diǎn)2的X坐標(biāo) “JDX2” =2919624.990按EXE鍵 4輸入交點(diǎn)2的Y坐標(biāo) “JDY2” =477022.978按EXE鍵5. 輸入曲線半徑R=4500按EXE鍵 6輸入緩和曲線長(zhǎng)L=340按EXE鍵 7輸入曲線轉(zhuǎn)向角A=26°2423”按EXE鍵 8輸入交點(diǎn)2的ZH點(diǎn)里程ZH：LC=348692.554按EXE鍵 顯示交點(diǎn)2的T=1225.972（切線長(zhǎng)）按EXE鍵；顯示交點(diǎn)2的S=2413.952（曲線長(zhǎng)）；按EXE鍵E直緩點(diǎn)縱坐標(biāo)，顯示交點(diǎn)2的ZH：X=2920409.07

21、4 N直緩點(diǎn)橫坐標(biāo)，顯示交點(diǎn)2的ZH：Y=476080.52279輸入交點(diǎn)2的測(cè)設(shè)點(diǎn)里程DK?+?=349032.554；按EXE鍵10輸入交點(diǎn)2中線點(diǎn)或左側(cè)邊樁的距離，如計(jì)算中樁點(diǎn)坐標(biāo)輸入PJU=0，如需計(jì)算左側(cè)15米邊樁的坐標(biāo)則輸入PJU=-15米；按EXE鍵11屏幕顯示Z=1；Y=2？提示，本例為左偏曲線故輸入1，按EXE鍵12. 顯示X=2920194.946，按EXE鍵； 顯示Y=476344.5952，按EXE鍵。顯示QXFWJ=127°3540.43”， 即DK349+032.554中樁點(diǎn)的切線方位角為127°3540.43”按EXE鍵若需計(jì)算DK349+03

22、2.554左側(cè)15米邊樁點(diǎn)的坐標(biāo)，按EXE鍵，屏幕提示DK?+?= 349032.554；再按EXE鍵，提示PJU=？，輸入-15米，再按EXE鍵，顯示Z=1；Y=2？輸入數(shù)字1按EXE鍵， 顯示X=2920206.831，按EXE鍵；顯示476353.746，按EXE鍵。顯示QXFWJ=127°3540.43”，按EXE鍵；13若需要再計(jì)算其他線路中線點(diǎn)或邊樁點(diǎn)的坐標(biāo)，重復(fù)811步驟操作，不必再輸入各已知數(shù)據(jù)，僅輸入需要計(jì)算點(diǎn)的里程或橫向偏距及Z=1；Y=2？的參數(shù)即可。豎曲線程序fx-4800p一、程序名：SQX二、程序:RDNCI1：FI2：T=0.5Abs(F-C)RE=T2

23、/(2R)VTU=1,AO=-1LbI1：L：LD=>Goto2L-DT=>H=M+ (L-D)F -V(T-L +D)2/(2R)：Goto3：H=M+ (L-D)F：Goto3H=M （L-D）F：Goto3LbI2：D-LT=>H=M-（D-L）C -V（T-D+ L）2/(2R)：Goto3：H=M-（D-L）C：Goto3LbI3：HGoto1三、程序執(zhí)行1、CI1 前進(jìn)方向輸入變坡點(diǎn)后段坡率2、FI2 輸入變坡點(diǎn)前段坡率以路線前進(jìn)方向3、R 輸入豎曲線半徑。計(jì)算出T豎曲線切線長(zhǎng)，E豎曲線外矢距。4、L 輸入所求點(diǎn)的樁號(hào)（里程）5、D 輸入變坡點(diǎn)樁號(hào)（里程）6、M

24、輸入變坡點(diǎn)高程。7、H 以路線計(jì)算出所求點(diǎn)高程8、VTU=1,AO=-1 凸曲線時(shí)輸入1，凹曲線時(shí)輸入-1四、說(shuō)明1、計(jì)算凸曲線，所求點(diǎn)樁號(hào)（里程）不小于后一個(gè)凹曲線臨界（終點(diǎn)）樁號(hào)，不大于前一個(gè)凹曲線臨界（起點(diǎn)）樁號(hào)。2、計(jì)算凹曲線，所求點(diǎn)樁號(hào)（里程）不小于后一個(gè)凸曲線臨界（終點(diǎn)）樁號(hào)，不大于前一個(gè)凸曲線臨界（起點(diǎn)）樁號(hào)。曲線任意里程中邊樁坐標(biāo)正反算(CASIO fx-4800P計(jì)算器)程序    一、程序功能    本程序由一個(gè)主程序(TYQXJS)和兩個(gè)子程正算子程序(SUB1)、反算子程序(SUB2)序構(gòu)成，可以根據(jù)曲線段直線

25、、圓曲線、緩和曲線（完整或非完整型）的線元要素（起點(diǎn)坐標(biāo)、起點(diǎn)里程、起點(diǎn)切線方位角、線元長(zhǎng)度、起點(diǎn)曲率半徑、止點(diǎn)曲率半徑）及里程邊距或坐標(biāo)，對(duì)該曲線段范圍內(nèi)任意里程中邊樁坐標(biāo)進(jìn)行正反算。另外也可以將本程序中核心算法部分的兩個(gè)子程序移植到其它相關(guān)的程序中，用于對(duì)曲線任意里程中邊樁坐標(biāo)進(jìn)行正反算。本程序也可以在CASIO fx-4500P計(jì)算器及 CASIO fx-4850P計(jì)算器上運(yùn)行。    特別申明：(1). 適用于弧長(zhǎng)小于2倍半徑的各種線元坐標(biāo)正反算，精度優(yōu)于1mm；        

26、;      (2). 在引用該核心計(jì)算部分時(shí)，請(qǐng)注明來(lái)源。    二、源程序1.主程序(TYQXJS)"1.SZ => XY"："2.XY => SZ"：N：U"X0"：V"Y0"：O"S0"：G"F0"：H"LS"：P"R0"：R"RN"：Q：C=1÷P：D=(P-R)÷(2HPR)：E=180

27、7;：N=1=>Goto 1：>Goto 2Lbl 1：SZ：SZ：W=Abs(S-O)：Prog "SUB1"：X"XS"=XY"YS"=YF"FS"=F-90Goto 1Lbl 2：XY：XY：I=X：J=Y：Prog "SUB2"：S"S"=O+WZ"Z"=ZGoto 22. 正算子程序(SUB1)A=0.1739274226：B=0.3260725774：K=0.0694318442：L=0.3300094782：F=1-L：M=1-K

28、：X=U+W(Acos(G+QEKW(C+KWD)+Bcos(G+QELW(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

29、)cosL-(I-X)sinL：AbsZ<1E-6=>Goto1：>W=W+Z：Goto 0Lbl 1：Z=0：Prog "SUB1"：Z=(J-Y)÷sinF三、使用說(shuō)明1、規(guī)定   (1) 以道路中線的前進(jìn)方向（即里程增大的方向）區(qū)分左右；當(dāng)線元往左偏時(shí)，Q=-1；當(dāng)線元往右偏時(shí)，Q=1；當(dāng)線元為直線時(shí)，Q=0。   (2) 當(dāng)所求點(diǎn)位于中線時(shí)，Z=0；當(dāng)位于中線左鍘時(shí)，Z取負(fù)值；當(dāng)位于中線中線右側(cè)時(shí)，Z取正值。   (3) 當(dāng)線元為直線時(shí)，其起點(diǎn)、止點(diǎn)的曲率半徑為無(wú)窮大，以10的4

30、5次代替。   (4) 當(dāng)線元為圓曲線時(shí)，無(wú)論其起點(diǎn)、止點(diǎn)與什么線元相接，其曲率半徑均等于圓弧的半徑。   (5) 當(dāng)線元為完整緩和曲線時(shí)，起點(diǎn)與直線相接時(shí)，曲率半徑為無(wú)窮大，以10的45次代替；與圓曲線相接時(shí)，曲率半徑等于圓曲線的半徑。止點(diǎn)與直線相接時(shí)，曲率半徑為無(wú)窮大，以10的45次代替；與圓曲線相接時(shí)，曲率半徑等于圓曲線的半徑。   (6) 當(dāng)線元為非完整緩和曲線時(shí)，起點(diǎn)與直線相接時(shí)，曲率半徑等于設(shè)計(jì)規(guī)定的值；與圓曲線相接時(shí)，曲率半徑等于圓曲線的半徑。止點(diǎn)與直線相接時(shí)，曲率半徑等于設(shè)計(jì)規(guī)定的值；與圓曲線相接時(shí)，曲率半徑等于圓曲線

31、的半徑。2、輸入與顯示說(shuō)明   輸入部分：   1. SZ => XY   2. XY = > SZ    N ?  選擇計(jì)算方式，輸入1表示進(jìn)行由里程、邊距計(jì)算坐標(biāo) ；輸入2表示由坐標(biāo)反算里程和邊距。   X0 ？線元起點(diǎn)的X坐標(biāo)    Y0 ？線元起點(diǎn)的Y坐標(biāo)    S0 ？線元起點(diǎn)里程   F0 ？線元起點(diǎn)切線方位角   LS ？線元長(zhǎng)度   R0 ？線元起點(diǎn)曲率半徑&

32、#160;  RN ？線元止點(diǎn)曲率半徑   Q ？  線 元左右偏標(biāo)志(左偏Q=-1，右偏Q=1，直線段Q=0）   S ？ 正算時(shí)所求點(diǎn)的里程   Z ？正算時(shí)所求點(diǎn)距中線的邊距(左側(cè)取負(fù)，值右側(cè)取正值，在中線上取零)   X ？反算時(shí)所求點(diǎn)的X坐標(biāo)   Y ？反算時(shí)所求點(diǎn)的Y坐標(biāo)   顯示部分：   XS=×××    正算時(shí)，計(jì)算得出的所求點(diǎn)的X坐標(biāo)  

33、; YS=×××    正算時(shí)，計(jì)算得出的所求點(diǎn)的Y坐標(biāo)   FS=×××    正算時(shí)，所求點(diǎn)對(duì)應(yīng)的中線點(diǎn)的切線方位角   S=×××    反算時(shí)，計(jì)算得出的所求點(diǎn)的里程   Z=×××    反算時(shí)，計(jì)算得出的所求點(diǎn)的邊距高程源公式|編程計(jì)算器卡西歐4800 2008-10-11 01:02:

34、42  作者：LEE.CI  來(lái)源：網(wǎng)易博客  瀏覽次數(shù)：165  文字大?。骸敬蟆俊局小俊拘　?超高及高程（源程序）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+TN<0 D=O+(SN)F: SN<0 D=O+(SN)FU(SN+T)2÷2÷V: SNT D=O+(SN)GU(TS+N)2&#

35、247;2÷V: D=O+(SN)G D“ H=” L4: R< E= : M= Goto 1: R< E= : M= Goto 1:  L5: Lbl 1 : M L C=L: Goto 2: C=M  L6: Lbl 2 : P X >0 Goto 3 : Goto 4  L7: Lbl 4 : S Z I=0.02 : Goto 5: S<Z+C I=0.02(0.02EPX)÷C×(SZ) : Goto 5: S<HC I=EPX : Goto 5: S<H I=EPX+(0.02EPX)&

36、#247;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: S<Z+C I=0.02+(EPX0.02)÷(CZ1)×(SZZ1) : Goto 5: S<HC I=EPX : Goto 5: S<HZ1 I=EPX(EPX0.02)÷(CZ1)×(SH+C) : Goto 5: I=0.02 : Goto 5L9: Lbl 5 : B : W“HB”=DIBL10:

37、Goto 6說(shuō)明： R:平曲線半徑 L“LS”:緩和曲線長(zhǎng) P：偏轉(zhuǎn)，Z=1,Y=1X“XL”：線路（Z,Y）左線1,右線1 Z“ZH”: 直緩點(diǎn)里程H“HZ” ：緩直點(diǎn)里程 V“RS”：豎曲線半徑 F“I 1” ：前坡堵（帶符號(hào)）G“I2”： 后坡度值（帶符號(hào)） O“BPH”：變坡點(diǎn)高程 N“BP”： 變坡點(diǎn)里程T： 切線長(zhǎng) E：最大超高值 C ：漸變段長(zhǎng)度（LC） I： 橫坡值 M：最大漸變長(zhǎng)度（Lmax） B ：距中心距離 W“HB”：邊部高程 字四、QXZB(曲線坐標(biāo)) A”ZH(ZH)=:B”X(ZH)=:C”Y(ZH)= ”X(HZ)=:E”Y(HZ)=:F”FWJ-1=:G”FW

38、J-2=:H”R=:I”L1=:J”L=:K”L2=:L”ZJ(Z=1,Y=-1)= O=180/pi (pi為圓周率，電腦中找不到此符號(hào)以此代替) LbI A MN M”ZHUANG HAO= P=M-A MM>(A+I+J+K)?GOTO A N”ZHUANG JU(Z+,Y-)= M>(A+I+J)?GOTO D M>(A+I)?GOTO C LbI B(第一緩和曲線) X=B+(P-(P5÷40÷H2÷I2)COSF+(P3÷6÷H÷I-P7÷336÷H3÷I3)LSINF Y=

39、C+(P-(P5÷40÷H2÷I2)SINF-(P3÷6÷H÷I-P7÷336÷H3÷I3)LCOSF X=X+NCOS(F-(P2÷2÷H÷I)OL-90) Y=Y+NSIN(F-(P2÷2÷H÷I)OL-90) GOTO A LbI C(圓曲線) X=B+(HSIN(P-I)÷H+I÷2÷H)O)+(I÷2-I3÷240÷H2)COSF-(H(1-COS(P-I)÷H+I&

40、#247;2÷H)O)+ I2÷24÷H)(-1)LSINF Y=C+(HSIN(P-I)÷H+I÷2÷H)O)+(I÷2-I3÷240÷H2)SINF+(H(1-COS(P-I)÷H+I÷2÷H)O)+I2÷24÷H)(-1)LCOSF X=X+NCOS(F-(P-I)÷H+I÷2÷H)OL-90) Y=Y+NSIN(F-(P-I)÷H+I÷2÷H)OL-90) GOTO A LbI D(第二緩

41、和曲線) Q=A+I+J+K-M X=D-(Q-Q5÷40÷H2÷K2)COSG+(Q3÷6÷H÷K-Q7÷336÷H3÷K3)LSING Y=E-(Q-Q5÷40÷H2÷K2)SING-(Q3÷6÷H÷K-Q7÷336÷H3÷K3)LCOSG X=X+NCOS(F-(I÷2÷H+J÷H+K÷2÷H-Q2÷2÷H÷K)OL-90) Y=Y

42、+NSIN(F-(I÷2÷H+J÷H+K÷2÷H-Q2÷2÷H÷K)OL-90) GOTO A 變量說(shuō)明： 以上變量依次為： ZH(ZH)=直緩點(diǎn)樁號(hào)，無(wú)緩和曲線則為直圓點(diǎn)； X(ZH)=直緩點(diǎn)X坐標(biāo)；無(wú)緩和曲線則為直圓點(diǎn)； Y(ZH)=直緩點(diǎn)Y坐標(biāo)，無(wú)緩和曲線則為直圓點(diǎn)； X(HZ)=緩直點(diǎn)X坐標(biāo)，無(wú)緩和曲線則為圓直點(diǎn)； Y(HZ)=緩直點(diǎn)Y坐標(biāo)，無(wú)緩和曲線則為圓直點(diǎn)； FWJ-1=第一切線方位角； FWJ-2=第二切線方位角； R=半徑； L1=第一緩和曲線長(zhǎng)； L=圓曲線長(zhǎng)； L2=第二緩和曲線長(zhǎng)； ZJ(Z

43、=1,Y=-1)=轉(zhuǎn)角形式(左=1，右=-1) ZHUANG HAO=所求點(diǎn)樁號(hào) ZHUANG JU(Z+,Y-)=邊樁樁距，左邊樁為正，右為負(fù)，中樁0 我已經(jīng)用過(guò)了，沒發(fā)現(xiàn)什么問(wèn)題，請(qǐng)大家提意見。 1、極坐標(biāo)法放樣Prog：FYLb1 0：A“X0”：B“Y0”：I=0：J=0：Pol（C“XA”-A），（D“YA”-B）：J<0=>G“FW- OA”=J+360L“L0”=IGoto 1：> G“FW O-A”=JL“L0”=ILb1 1：EQ：E“Xi”：Q“Yi”：Pol（E-A），（Q-B）：J<0=>J=J+360：Goto 2：> Goto 2

44、Lb1 2： F“FW-OB”=JL=I0=F-G：O<0=>O“BJ”=O+360Goto 3：> O “BJ” Lb1 3：P=O-180Goto 1注：a、輸入：（X0、Y0）、（XA、YA）測(cè)站點(diǎn)坐標(biāo)、后視點(diǎn)坐標(biāo) Xi、Yi 放樣點(diǎn)坐標(biāo)b、輸出：FW-OA測(cè)站至后視邊方位角、L0后視邊長(zhǎng) FW-OB測(cè)站至放樣點(diǎn)方位角、L放樣邊長(zhǎng)BJ后視邊置零，放樣點(diǎn)順時(shí)針撥角P偏角（+為右偏、-為左偏）本值用于計(jì)算路線偏角2、公路豎曲線高程計(jì)算程序Prog:SQXLbl A:A“+（-）i1”:B“+（-）i2” W=（B-A）÷100：R：T=Abs（RW）÷2

45、：L=T*2：E=T2÷（2R）：K“JD K+”：G“JD H”：C=K-T：D=K+T：Lbl 0：J“Ki+”：J<0=>Goto 1：> Goto 2Lb1 1：“Out QX1”：H=G-(K-J)A÷100Goto 5Lb1 2：J>D=>Goto 4 W<0=>F=-1W>0=>F=1J>K=>Goto 3H=G-（K-J）A÷100+F(J-C)2÷（2R）Goto 5Lb1 3：H=G+（J-K）B÷100+F（D-J）2÷（2R）Goto 5Lb1

46、4：“OUT QX2”：H=G+（J-K）B÷100Goto 5Lb1 5：M“DHi”：H=H+M注：a、公式：L=|R(i2-i1)| 、T=L÷2、E=T2÷（2R）、h=l2÷(2R)b、功能：已知前后坡度%、豎曲線半徑，計(jì)算各樁高程。c、輸入：i1、i2前、后坡度%（上坡+、下坡-） JD K+、 JD H、R變坡點(diǎn)樁號(hào)、高程、豎曲線半徑 Ki+計(jì)算樁號(hào) DHi計(jì)算點(diǎn)與設(shè)計(jì)線高差（高+、低-）d、輸出：H待計(jì)算點(diǎn)高程3、直線坐標(biāo)計(jì)算Prog：ZXLb1 0：M“K0”：A“X1”：B“Y1”：C“X2”：D“Y2”：I=0：J=0：Pol（C-

47、A），（D-B）：L=IJ<0=>F”FW”=360+JGoto 1：>F”FW”=J Lb1 1：KZQ:K“Ki”：Z“Lb”：Q=90：Q“ANG”：L=K-M：Rec（L，F(xiàn)）：X=A+I：Y=B+J：Rec（Z，F(xiàn)+Q）：O=I：P=J：X=X+OY=Y+P Goto 1注：a、公式：X=X0+L*cos(Fw)、Y=Y0+L*sin(Fw)b、功能：已知直線起終點(diǎn)坐標(biāo)，計(jì)算中間樁中邊樁坐標(biāo)c、輸入： K0、（X1、Y1）、（X2、Y2）起點(diǎn)樁號(hào)、起點(diǎn)坐標(biāo)、終點(diǎn)坐標(biāo)Ki、Lb、Ang 待求樁號(hào)、邊樁距（左-、右+）、邊樁與中線右夾角（90）d、輸出：L、FW直線長(zhǎng)度

48、、起點(diǎn)至終點(diǎn)方位角 X、Y計(jì)算點(diǎn)坐標(biāo)a多邊形面積計(jì)算（一）Prog：AreaLb1 0：S=0：L=0：E“X0”：F“Y0”：A=E：B=FLb1 1：CD：C“Xi： XiX0，Yi=Y0 END”：D“Yi”：M=（C-E）2+（D-F）2）：S=S+（AD-CB）：L=L+（C-A）2+（D-B）2）：A=C：B=D：S=0Goto 1：Lb1 3：M0=Goto 1：S“Area”=Abs S÷2L注：a、公式：Area =（XiYi+1-Xi+1Yi）÷2b、功能：已知多邊形各轉(zhuǎn)點(diǎn)坐標(biāo)，計(jì)算多邊形面積。c、輸入：X0、Y0起點(diǎn)坐標(biāo) Xi、Yi各點(diǎn)坐標(biāo)（按順序輸入

49、、最后要閉合即最后一點(diǎn)要輸入起點(diǎn)坐標(biāo)）d、輸出：Area面積 L周長(zhǎng)4、對(duì)稱標(biāo)準(zhǔn)曲線中邊樁坐標(biāo)計(jì)算Prog：QXLb1 A: M“JDK+”:A“PJ=>R（+），L（-）”:R：W“LS”: A0=>G=1:Goto E: >G=-1:Lb1 E:Defm 10:A=Abs A:B=90W÷R÷：P=W2÷（6R）-W4÷（336R3）-R（1-cosB）:Q=W-W3÷(40R2)-RsinB：T=（R+P）tan（A÷2）+Q：L=RA÷180+W:E=(R+P)÷cos(A÷2)

50、-R:N=M-T:Z5=N+W:Z=N+L÷2:C=N+L-W:D=N+L:A=AG:U:U“ZD-2,XY-1”=2=>TELN“ZH”Z5“HY”Z“QZ”C“YH”D“HZ”Goto A：> Goto B Lbl B：F“FW”：V“JD-X”：H“JD-Y”：I=0：J=0：Rec（T，F(xiàn)+180）：Z1“ZH X”=I+V：Z2“ZH Y”=J+H：Rec（T，F(xiàn)+A）：Z3“HZ X”=I+V：Z4“HZ Y”=J+H：L=0Lbl 0：E=90：KSE：K“JS K+”：S“LB”：E“ANG”：S=0=>U=1：Goto C：>U=S÷

51、;（Abs S）Lbl C：S=Abs S：W=0=> Goto D：>K<N=> Goto 1：>K<Z5 => Goto 2：>K<D=> Goto 4：> Goto 5Lbl D：K<N=>Goto 1：>K<D=>Goto 3：> Goto 5Lbl 1：Rec（M-K，F(xiàn)+180）：X=V+I：Y=J+H：Rec（S，F(xiàn)+180-（180-E）U）：X“Xi”=X+IY“Yi”=Y+JS=SU：Goto 7Lbl 2：Z=K-N：O=90Z2÷R÷W÷：

52、X=Z-Z5÷40÷R2÷W2+Z9÷3456÷R4÷W4：Z=Z3÷6÷R÷W-Z7÷336÷R3÷W3+Z11÷42240÷R5÷W5：Lbl 6：Rec（X，F(xiàn)）：X=Z1+I：Y=Z2+J：Rec（Z，F(xiàn)+90G）：X=X+I：Y=Y+J：Rec（S，F(xiàn)+OG+EU）：X“Xi”=X+IY“Yi”=Y+JS=SU：Goto 7Lbl 3：Z=K-Z5：O=（K-N-W÷2）*180÷÷R：X=Q+Rsin（B

53、+180Z÷÷R）：Goto 6Lbl 4：Z=D-K：O=90Z2÷R÷W÷：X=Z-Z5÷40÷R2÷W2+Z9÷3456÷R4÷W4：Z=Z3÷6÷R÷W-Z7÷336÷R3÷W3+Z11÷42240÷R5÷W5：Rec（X，F(xiàn)+A+180）：X=Z3+I：Y=Z4+J：Rec（Z，F(xiàn)+A+180-90G）：X=X+I：Y=Y+J：Rec（S，F(xiàn)+A+180-OG-（180-E）U）：X“X

54、i”=X+IY“Yi”=Y+JS=SU：Goto 7Lbl 5：Rec（K-D+T），（F+A）：X=V+I：Y=H+J：Rec（S，F(xiàn)+A+EU）：X“Xi”=X+IY“Yi“=Y+J S=SU： Goto 7Lbl 7：Goto 0注：輸入：JDK+、PJ、R、LS、交點(diǎn)樁號(hào)、路線偏角（右+、左-）、半徑、緩和曲線長(zhǎng)“ZD-2,XY-1”要顯示主點(diǎn)要素輸入2，否則進(jìn)入坐標(biāo)計(jì)算FW、（JD X、JD Y）起始邊方位角（即ZHJD方位角）、交點(diǎn)坐標(biāo)JS K+、LB、Ang待計(jì)算樁號(hào)、邊樁距（右+、左-）、邊樁與路中線右夾角（默認(rèn)90度）輸出：T、E、L、ZH、HY、YH、HZ、曲線要素、 Xi

55、,Yi待求點(diǎn)坐標(biāo)5、公路任何曲線通用中邊樁坐標(biāo)正反算（復(fù)化辛卜生公式也稱積木法）Prog：HQXDefm 4V“V=1=>XY： 2=>Ki”Lbl 0：KLWLbl 4：Q“one=0：1、2” =0=> Goto 1Q=1=>Prog“ZD”（線元要素?cái)?shù)據(jù)文件名） Q=3=>Prog“要素表文件名” （注：本句是為了全線及多條線路計(jì)算用，如僅單個(gè)可省略本句）Lbl 1：A"K1"：B"K2"：C"FW0"：D"R1"：E"R2"：F"X0"：

56、G"Y0"D0=>I=1÷D：=>I=DE0=>J=1÷E：=>J=EAbsD+AbsE=0=>M=2：=>M=16V=2=>L=0：W=90K“Ki”：L“L-BZ”：W"ANG"：N=0：Z1=0：Z2=0：Z3=0：Z4=0 進(jìn)入坐標(biāo)迭代計(jì)算Lbl 2：N=N+1：H=2(K-A)÷M：R=NH÷2+A：R=C+180÷×(I+(J-I)÷2÷(B-A)×(R-A)×(R-A)：Int(N÷2)=N

57、÷2=>Z1=Z1+cosR：Z2=Z2+sinR：=>Z3=Z3+cosR：Z4=Z4+ sinRN=M=>Goto 3：>Goto 2Lbl 3：X=F+H÷6×(cosC+4Z3+2Z1-cosR)+Lcos(R+W)：Y=G+H÷6×(sinC+4Z4+2Z2-sinR)+Lsin(R+W)：V=2=>Goto 6X"Xi="Y"Yi="R"FW"=R-360Intg(R÷360) Q=0=>KLW：Goto 1 ：> Got

58、o 0Lbl 6： Pol(T"XF"-X,U"YF"-Y) 進(jìn)入樁號(hào)求算O=Icos(J-R)：AbsO0.0001=>O"DL"K=K+O O"LP"=Isin(J-R)TU：Goto 6：=>K=K+O：L=0：Goto 4線元要素?cái)?shù)據(jù)文件格式：K>起點(diǎn)樁號(hào)=> A"K1"=？：B"K2"=？：C"FW0"=？：D"R1"=？：E"R2"=？：F"X0"=？：G"Y0"=？：K> B=>A=B：B=？（下曲線要素）注：a、功能：利用公路通用復(fù)化辛卜生積分公式進(jìn)行中邊樁坐