求:卡西歐計算機(jī)4850 道路測量程序```````

本程序是在yshf及hangniu1973兩位師傅的成果上作的一次改動，修改內(nèi)容主要有一下幾點(diǎn):1、在變量符號上修改，目的是更能與工地實(shí)際結(jié)合；2、把原主程序分為兩個，更能靈活調(diào)用；3、加入放樣程序，做到坐標(biāo)計算到放樣一體化；4、使得整個測量放樣過程更加簡單。附件點(diǎn)擊下載:GAUSSLE坐標(biāo)正反算fx－4850程序源程序1.正算主程序 GSZSI"X0":S"Y0":O"K0":G"F0":H"KN":P"R0":R"RN":Q”Q-Z +Y” :D=P-R÷2H-OPR:KL”L-Z +Y” :M”ANGYJJ”=90:注:此處若不給M賦值，則可計算斜交點(diǎn)J=AbsK-O:Prog"SUB1":”FWJ=”:F=F-M:”X=”:X=X:Pause0: ”Y=”:Y=Y2. 反算主程序 GSFSXY:Z[2]=X:Z[3] =Y:I"X0":S"Y0":O"K0":G"F0":H"KN":P"R0":R"RN":Q”Q-Z +Y” :D=P-R÷2H-OPR:J=AbsY-ScosG-90－X-IsinG-90:L=0:M”MYJJ”=90:Lbl 0:Prog "SUB1":L=Z[3]-YcosG-90+QJ1÷P+JD×180÷π-Z[2]-XsinG-90+QJ1÷P +JD ×180÷π:AbsL<1E-6=>Goto1:≠>J=J+L:Goto 0Δ←┘Lbl 1:L=0:Prog "SUB1":L=Z[3]-Y÷sinF:”K=”:K=O+J:Pause0:”L=”:L=L3. 正算子程序SUB1Defm 4:A=0.1184634425:B=0.2393143352:Z[4]=0.2844444444:C=0.0469100770:E=0.2307653449:Z[1]=0.5:X=I+JAcosG+QCJ1÷P+CJD×180÷π+BcosG+QEJ1÷P+EJD×180÷π+Z[4]cosG+QZ[1]J1÷P+Z[1]JD×180÷π+BcosG+Q1-EJ1÷P+1-EJD×180÷π+AcosG+Q 1-CJ1÷P+1-CJD ×180÷π:Y=S+JAsinG+QCJ1÷P+CJD×180÷π+BsinG+QEJ1÷P+EJD×180÷π+Z[4]sinG+QZ[1]J1÷P+Z[1]JD×180÷π+BsinG+Q1-EJ1÷P+1-EJD×180÷π+AsinG+Q 1-CJ1÷P+1-CJD ×180÷π:F=G+QJ1÷P+JD ×180÷π+M:X=X+LcosF:Y=Y+LsinF4. 曲線元要素數(shù)據(jù)庫:DAT-MK≥O=>K＜H=> I=**:S=**:O=**:G=**:H=**:P=**:R=**:Q=**SS←┘K≥O=>K＜H=> I=**:S=**:O=**:G=**:H=**:P=**:R=**:Q=**SS←┘K≥O=>K＜H=> I=**:S=**:O=**:G=**:H=**:P=**:R=**:Q=**SS←┘K≥O=>K＜H=> I=**:S=**:O=**:G=**:H=**:P=**:R=**:Q=**SS←┘K≥O=>K≤H=> I=**:S=**:O=**:G=**:H=**:P=**:R=**:Q=**SS←┘……………………………K≥O=>K≤H=> I=**:S=**:O=**:G=**:H=**:P=**:R=**:Q=**SS←┘ 注:如有多個曲線元要素繼續(xù)添加入數(shù)據(jù)庫DAT-M中5、M線坐標(biāo)正算組合程序 MG-ZBProg”DAT-M”:Prog”GSZS”6、M線坐標(biāo)計算－放樣組合程序 MG-FYProg”MG-ZB”:Prog”LTKZD”: Prog”FY”7、M線坐標(biāo)反算組合程序 M-GSFBProg”DAT-M”:Prog”GSFS”說明:一、程序功能及原理1.功能說明:本程序由兩個主程序——正算主程序GSZS、反算主程序GSFS和兩個子程——正算子程序SUB1、線元數(shù)據(jù)庫DAT-M構(gòu)成，可以根據(jù)曲線段——直線、圓曲線、緩和曲線完整或非完整型的線元要素起點(diǎn)坐標(biāo)、起點(diǎn)里程、起點(diǎn)切線方位角、終點(diǎn)里程、起點(diǎn)曲率半徑、止點(diǎn)曲率半徑及里程邊距或坐標(biāo)，對該曲線段范圍內(nèi)任意里程中邊樁坐標(biāo)進(jìn)行正反算。本程序可以在CASIO fx-4800P計算器及 CASIO fx-4850P計算器上運(yùn)行。由于加入了數(shù)據(jù)庫DAT-M，可實(shí)現(xiàn)坐標(biāo)正反算的全線貫通。組合程序5可實(shí)現(xiàn)M線的正算貫通，組合程序7可實(shí)現(xiàn)M線的反算貫通，組合程序6可實(shí)現(xiàn)坐標(biāo)計算到放樣一體化。2．計算原理:利用Gauss-Legendre 5點(diǎn)通用公式計算線路中邊樁坐標(biāo)并計算放樣數(shù)據(jù)。利用待求點(diǎn)至線元起點(diǎn)切線作垂線，逐次迭代趨近原理反算里程及邊距。二、使用說明1、規(guī)定1 以道路中線的前進(jìn)方向即里程增大的方向區(qū)分左右；當(dāng)線元往左偏時， Q=-1；當(dāng)線元往右偏時，Q=1；當(dāng)線元為直線時，Q=0。2 當(dāng)所求點(diǎn)位于中線時，L=0；當(dāng)位于中線左側(cè)時，L取負(fù)值；當(dāng)位于中線右側(cè)時，L取正值。3 當(dāng)線元為直線時，其起點(diǎn)、止點(diǎn)的曲率半徑為無窮大，以10的45次代替。4 當(dāng)線元為圓曲線時，無論其起點(diǎn)、止點(diǎn)與什么線元相接，其曲率半徑均等于圓弧的半徑。5 當(dāng)線元為完整緩和曲線時，起點(diǎn)與直線相接時，曲率半徑為無窮大，以10的45次代替；與圓曲線相接時，曲率半徑等于圓曲線的半徑。止點(diǎn)與直線相接時，曲率半徑為無窮大，以10的45次代替；與圓曲線相接時，曲率半徑等于圓曲線的半徑。6 當(dāng)線元為非完整緩和曲線時，起點(diǎn)與直線相接時，曲率半徑等于設(shè)計規(guī)定的值；與圓曲線相接時，曲率半徑等于圓曲線的半徑。止點(diǎn)與直線相接時，曲率半徑等于設(shè)計規(guī)定的值；與圓曲線相接時，曲率半徑等于圓曲線的半徑。7曲線元要素數(shù)據(jù)庫DAT-M可根據(jù)線型不同分為各個線元段輸入到DAT-M中，即分為直線段、緩和曲線、圓曲線等。8正算時可僅輸入里程和邊距及右交角可實(shí)現(xiàn)全線計算，但反算時只能通過首先輸入里程K值讀取數(shù)據(jù)庫DAT-M，計算器自動將里程K所在線元數(shù)據(jù)賦給反算主程序GSFS進(jìn)行試算，試算出的里程和邊距須帶入正算主程序GSZS中計算坐標(biāo)，若坐標(biāo)吻合則反算正確。2、輸入與顯示說明1輸入部分:X0 ？線元起點(diǎn)的X坐標(biāo)Y0 ？線元起點(diǎn)的Y坐標(biāo)K0 ？線元起點(diǎn)里程F0 ？線元起點(diǎn)切線方位角KN ？線元終點(diǎn)里程R0 ？線元起點(diǎn)曲率半徑RN ？線元止點(diǎn)曲率半徑Q ？ 線 元左右偏標(biāo)志左偏Q=-1，右偏Q=1，直線段Q=0K ？ 正算時所求點(diǎn)的里程L ？ 正算時所求點(diǎn)距中線的邊距左側(cè)取負(fù)值，右側(cè)取正值，在中線上取零ANG？正算邊樁時左右邊樁連線與線路中線的右交角 X ？ 反算時所求點(diǎn)的X坐標(biāo)Y ？ 反算時所求點(diǎn)的Y坐標(biāo)M ? 斜交右角線元要素數(shù)據(jù)庫中K≥O=>K＜H=>中的O和H分別為該段線元起點(diǎn)里程和終點(diǎn)里程A、 B、Z[4] 是Gauss-Legendre求積公式中的插值系數(shù)C 、E、Z[1] 是Gauss-Legendre求積公式中的求積節(jié)點(diǎn)2顯示部分: X=××× 正算時，計算得出的所求點(diǎn)的X坐標(biāo)Y=××× 正算時，計算得出的所求點(diǎn)的Y坐標(biāo)K=××× 反算時，計算得出的所求點(diǎn)的里程L=××× 反算時，計算得出的所求點(diǎn)的邊距三、算例某匝道的由五段線元直線+完整緩和曲線+圓曲線+非完整緩和曲線+直線組成，各段線元的要素起點(diǎn)里程S0、起點(diǎn)坐標(biāo)X0 Y0、起點(diǎn)切線方位角F0、線元長度LS、起點(diǎn)曲率半徑R0、止點(diǎn)曲率半徑RN、線元左右偏標(biāo)志Q如下:S0 X0 Y0 F0 LS R0 RN Q500.000 19942.837 28343.561 125 16 31.00 269.256 1E45 1E45 0769.256 19787.340 28563.378 125 16 31.00 37.492 1E45 221.75 -1806.748 19766.566 28594.574 120 25 54.07 112.779 221.75 221.75 -1919.527 19736.072 28701.893 91 17 30.63 80.285 221.75 9579.228 -1999.812 19744.038 28781.659 80 40 50.00 100.000 1E45 1E45 0注:該算例中線元要素Ls為程序修改前須輸入的線元長度，程序修改后改為輸入線元終點(diǎn)里程KN 放樣程序 FYU”XJ”:V”YJ”:W”XH”:Z”YH”:XY:I=W-U:J=Z-V:PolI，J:J<0=>J=J+360Δ“HSJL=”:IO”HF”=J:X=0=>U=U+LCosJ+P”μ”:V=V+LSinJ+P”μ”:”X=”:U:Pause0:”Y=”:V:≠>I=X-U:J=Y-V: PolI，J:J<0=>J=J+360:ΔH”QF”=J:J=H-O:J<0=>J=J+360:Δ“ANG=”:J→DMS 該處不輸Pause0語句亦可實(shí)現(xiàn)“ANG=”和“QSJL=”的同時出現(xiàn)“QSJL=”:I凌鐵大橋控制點(diǎn)數(shù)據(jù)庫 LTKZDA=21:B=22:C=23:D=24:E=25:F=26:G=27:H=28:I=29←┘T “ZJDH”:←┘T=11=>U=****:V=****S←┘T=12=>U=****:V=****S←┘T=13=>U=****:V=****S←┘T=15=>U=****:V=****S←┘T= 4 =>U=****:V=****S←┘T= A=>U=****:V=****S←┘T= B=>U=****:V=****S←┘T= C=>U=****:V=****S←┘……………………………注:如有多個控制點(diǎn)繼續(xù)添加入數(shù)據(jù)庫LTKZD中N “HSDH”:←┘N=11=>W=****:Z=****S←┘N=12=>W=****:Z=****S←┘N=13=>W=****:Z=****S←┘N=15=>W=****:Z=****S←┘N= 4 =>W=****:Z=****S←┘N= A=>W=****:Z=****S←┘N= B=>W=****:Z=****S←┘N= C=>W=****:Z=****S←┘……………………………注:如有多個控制點(diǎn)繼續(xù)添加入數(shù)據(jù)庫LTKZD中組合程式 LTFYProg ”LTKZD”: Prog ”FY” 說明:一、程序功能及原理1、功能說明:FY程序有兩個功能，當(dāng)前視點(diǎn)坐標(biāo)X輸入為零時，程式執(zhí)行轉(zhuǎn)點(diǎn)功能，即首先通過全站儀測定任意轉(zhuǎn)點(diǎn)距置鏡點(diǎn)距離L及置鏡點(diǎn)與轉(zhuǎn)點(diǎn)連線和置鏡點(diǎn)與后視點(diǎn)連線的夾角P”μ”并輸入計算器，即可計算出任意轉(zhuǎn)點(diǎn)坐標(biāo)；另一功能是坐標(biāo)放樣，即輸入放樣點(diǎn)X、Y坐標(biāo)，從而計算出與置鏡點(diǎn)距離 “QSJL=”以及全站儀歸零放樣撥角“ANG=”前視方位角－后視方位角。凌鐵大橋控制點(diǎn)數(shù)據(jù)庫LTKZD功能是與FY程序配合后，可實(shí)現(xiàn)坐標(biāo)放樣時僅輸入置鏡點(diǎn)、后視點(diǎn)的編號，從而減少其數(shù)據(jù)的大量輸入。2、計算原理:利用坐標(biāo)計算和方位角的基本知識以及計算器極坐標(biāo)轉(zhuǎn)換功能鍵POL二、使用說明1、規(guī)定當(dāng)前視點(diǎn)坐標(biāo)X輸入為零時，程式執(zhí)行轉(zhuǎn)點(diǎn)功能2、輸入與顯示說明1輸入部分:XJ ？ 置鏡點(diǎn)X坐標(biāo)YJ ？ 置鏡點(diǎn)Y坐標(biāo)XH ？ 后視點(diǎn)X坐標(biāo)YH ？ 后視點(diǎn)Y坐標(biāo)X ？ 放樣點(diǎn)X坐標(biāo)Y ？ 放樣點(diǎn)Y坐標(biāo)L ？ 實(shí)測轉(zhuǎn)點(diǎn)距置鏡點(diǎn)距離P”μ” ？ 實(shí)測置鏡點(diǎn)與轉(zhuǎn)點(diǎn)連線和置鏡點(diǎn)與后視點(diǎn)連線的夾角T ？ 置鏡點(diǎn)點(diǎn)號ZJDHN ？ 后視點(diǎn)點(diǎn)號HSDH11、12、13、15、4、A、B、C、D、E、F、G、H、I分別為凌鐵大橋控制點(diǎn)點(diǎn)號，計算輸入時按以上數(shù)字或字母輸入計算器即可2顯示部分HSJL=**** 后視距離置鏡點(diǎn)與后視點(diǎn)之間距離ANG= **** 全站儀歸零放樣撥角QSJL=**** 前視距離置鏡點(diǎn)與放樣點(diǎn)之間距離X= **** 轉(zhuǎn)點(diǎn)X坐標(biāo)Y= **** 轉(zhuǎn)點(diǎn)Y坐標(biāo)三、其它說明1、若需要編制另一線路控制點(diǎn)程序，僅需新建另一數(shù)據(jù)庫子程序LTKZD，新建時僅需改變上文中紅色顯示數(shù)據(jù)即可。2、程式中11~15、及A~I均為控制點(diǎn)編號，程式開頭首先給A、B賦值:A=21:B=22等，目的是為了在程式運(yùn)行要求輸入點(diǎn)號時能直接輸入字母A、B、C等。