--PAGE1-fx-5800Pfx-5800P矩陣編程計(jì)算器《工程測(cè)量與公路測(cè)設(shè)程序加油站操作手冊(cè)》fx-5800Pfx-4500PA、fx-4800Pfx-4850P中運(yùn)行！fx-5800Pfx-5800P20個(gè)常用工程測(cè)量與公路測(cè)量程序、程序名、子程序名及其功能列于表1。表1《工程測(cè)量與公路測(cè)設(shè)程序加油站》程序列表主程序名 需要調(diào)用的子程序名 功能說(shuō)明PR1計(jì)算一個(gè)起點(diǎn)至任意個(gè)端點(diǎn)的邊長(zhǎng)與方位角PR2極坐標(biāo)法放樣元素計(jì)算PR3經(jīng)緯儀視距法測(cè)圖計(jì)算記錄PR4單一閉附合圖根水準(zhǔn)路線近似平差計(jì)算PR5單一閉附合圖根導(dǎo)線近似平差計(jì)算PR6圖根支導(dǎo)線坐標(biāo)計(jì)算PR7測(cè)角前方交會(huì)坐標(biāo)計(jì)算PR8ANGLE測(cè)角后方交會(huì)坐標(biāo)計(jì)算PR9DIST-BEAR測(cè)邊后方交會(huì)點(diǎn)坐標(biāo)計(jì)算PR10COEFFICIENT,ZS,FS54北京系與80西安系高斯投影正算PR11COEFFICIENT,ZS,FS,DM-S54北京系與80西安系高斯投影反算PR12COEFFICIENT,ZS,FS54北京系與80西安系高斯投影換帶計(jì)算PR13X[K]-GC單交點(diǎn)單圓曲線偏角法中樁放樣要素計(jì)算PR14X[K]-XY單交點(diǎn)單圓曲線切線支距法中樁放樣要素計(jì)算PR15X[K]-GCXY單交點(diǎn)單圓曲線中樁坐標(biāo)計(jì)算PR16Z[10]-XY,X[K]-HGC單交點(diǎn)對(duì)稱基本形曲線偏角法中樁放樣要素計(jì)算PR17Z[10]-XY,X[K]-HXY單交點(diǎn)對(duì)稱基本形曲線切線支距法中樁放樣要素計(jì)算PR18Z[10]-XY,X[K]-HOXY單交點(diǎn)對(duì)稱基本形曲線中樁坐標(biāo)計(jì)算PR19X[K]-GQXY雙圓復(fù)曲線中樁測(cè)量坐標(biāo)計(jì)算PR20X[K]-XH豎曲線的計(jì)算表1中共有20個(gè)主程序，15個(gè)子程序，占用內(nèi)存21520字節(jié)，占計(jì)算器總內(nèi)存容量28500字節(jié)的75.5%。主程序名都是以字符“PR”開(kāi)頭加數(shù)字1~20命名的，其余字符開(kāi)頭的程序均為子程序，用戶只能運(yùn)行PR1~PR20主程序，不能直接運(yùn)行子程序。圖1兩臺(tái)fx-5800P間數(shù)據(jù)通訊卡西歐(上海)貿(mào)易有限公司向中國(guó)各地的卡西歐代理商提供一臺(tái)內(nèi)含20個(gè)程序的母機(jī)，如果您的計(jì)算器是在卡西歐代理商處購(gòu)買(mǎi)的，請(qǐng)您及時(shí)向代理商申請(qǐng)將母機(jī)內(nèi)的20個(gè)程序通過(guò)數(shù)據(jù)通訊的方式傳輸?shù)侥?gòu)買(mǎi)機(jī)器的內(nèi)存中。操作方法為：①使用SB-62數(shù)據(jù)線(需要另外購(gòu)買(mǎi))連接兩臺(tái)fx-5800P的3Pin數(shù)據(jù)通訊口，見(jiàn)圖1所示；②在您購(gòu)買(mǎi)的fx-5800P上按鍵打開(kāi)機(jī)器電源，按 (Receive)鍵使的機(jī)器處于接收數(shù)據(jù)狀態(tài)，屏幕顯示“eceiving，表示機(jī)器正等待接收母機(jī)發(fā)送來(lái)的數(shù)據(jù)；③在代理商的母機(jī)上按鍵打開(kāi)機(jī)器電源，按 鍵啟動(dòng)代理商的母機(jī)發(fā)送數(shù)據(jù)，屏幕顯示“…，表示正在發(fā)送數(shù)據(jù)。傳輸完全部20主程序及15個(gè)子程序大約需要時(shí)間56秒，程序傳輸完畢后，兩臺(tái)機(jī)器的屏幕都顯示“Complete2所示。圖2發(fā)送母機(jī)的全部數(shù)據(jù)到用戶機(jī)器的操作過(guò)程成后即可拔出數(shù)據(jù)線，按本操作手冊(cè)使用程序。請(qǐng)您將程序使用中的問(wèn)題與建議及時(shí)發(fā)送到qh-506@163.com，以便改進(jìn)我們的工作，更好地為您服務(wù)，卡西歐(上海)貿(mào)易有限公司感謝您的大力支持！1、程序PR1，計(jì)算一個(gè)起點(diǎn)至任意個(gè)端點(diǎn)的邊長(zhǎng)與方位角"X,YDISTAndBEARINGPR1"顯示程序標(biāo)題Deg:ClrStat:FreqOff:Fix3基本設(shè)置提示輸入0點(diǎn)的坐標(biāo)0K 變量計(jì)數(shù)清零Lbl0K+1K:Norm1:"n=":K 計(jì)數(shù)變量計(jì)數(shù)Fix3:"Xn(m)=,(0END)"?C提示輸入n點(diǎn)的x坐標(biāo)，輸入0結(jié)束程序計(jì)算WhileC0"Yn(m)="?D提示輸入n點(diǎn)的y坐標(biāo)Pol(C-A,D-B):Cls 調(diào)用極坐標(biāo)函數(shù)并清除屏幕顯示IfJ<0:ThenJ+360F:ElseJF:IfEnd"DISTO0n(m)=":I 顯示水平距離"BEAR0n(DMS)=":FDMS 以六十進(jìn)制度顯示計(jì)算出的方位角IListX[K]:FListY[K]存儲(chǔ)邊長(zhǎng)與方位角到統(tǒng)計(jì)串列Goto0:WhileEnd 提示重復(fù)輸入端點(diǎn)的坐標(biāo)"PR1END"x0時(shí)，停止程序運(yùn)行。1-1PR101，2，3，4點(diǎn)的邊長(zhǎng)與方位角，表中的紅色數(shù)字為計(jì)算結(jié)果。表1-1使用程序PR1計(jì)算邊長(zhǎng)與方位角案例點(diǎn)點(diǎn)號(hào)On(°′″)(m)DOn起訖點(diǎn)號(hào)y(m)x(m)O3885.6343114.47114281.7393592.881O→1621.108502235.623356.6683419.507O→2610.616150146.0933373.3972385.189O→3891.2012345458.8943968.1033005.750O→4136.4603071054.11完成計(jì)算后起點(diǎn)至端點(diǎn)的邊長(zhǎng)與方位角分別存儲(chǔ)在統(tǒng)計(jì)串列ListY中請(qǐng)按(REG)鍵進(jìn)入雙變量統(tǒng)計(jì)與回歸模式查看在ListY中的方位角是以十進(jìn)制度顯示的若要查看六十進(jìn)制的角度值，應(yīng)將光標(biāo)移動(dòng)到ListY的某個(gè)單元，按 鍵。2、程序PR2，極坐標(biāo)法放樣元素計(jì)算"POLARSETTING-OUTPR2"顯示程序標(biāo)題Deg:ClrStat:FreqOff:Fix3基本設(shè)置提示輸入測(cè)站點(diǎn)的坐標(biāo)提示輸入0方向點(diǎn)的坐標(biāo)Pol(C-A,D-B):Cls 調(diào)用極坐標(biāo)函數(shù)并清除屏幕顯示IfJ<0:ThenJ+360E:ElseJE:IfEnd判斷并計(jì)算方位角"DISTO0(m)=":I 顯示測(cè)站→0方向的水平距離"BEARO0(DMS)=":EDMS 以六十進(jìn)制度顯示測(cè)站→0方向的方位角1K:IListX[K]:EListY[K]存儲(chǔ)0方向的邊長(zhǎng)與方位角到統(tǒng)計(jì)串列Lbl0K+1K:Norm1:"n=":K 計(jì)數(shù)變量計(jì)數(shù)Fix3:"Xn(m)=,(0END)"?U提示輸入放樣點(diǎn)的x坐標(biāo)輸入0結(jié)束程序計(jì)算WhileU0"Yn=(m)"?V 提示輸入放樣點(diǎn)的y坐標(biāo) 調(diào)用極坐標(biāo)函數(shù)并清除屏幕顯示IfJ<0:ThenJ+360F:ElseJF:IfEnd判斷并計(jì)算方位角H計(jì)算0方向與n方向的水平夾角IfH<0:ThenH+360H:IfEnd判斷并計(jì)算水平夾角"DISTOn(m)=":I 顯示測(cè)站→n方向的水平距離"BEAROn(DMS)=":FDMS 以六十進(jìn)制度顯示測(cè)站→n方向的方位角"ANGLE0-O-n(DMS)=":HDMS 以六十進(jìn)制度顯示0方向與n方向的水平IListX[K]:FListY[K]存儲(chǔ)O→n方向的邊長(zhǎng)與方位角到統(tǒng)計(jì)串列HListFreq[K]存儲(chǔ)O→n方向的放樣水平角到統(tǒng)計(jì)串列Goto0:WhileEnd 提示重復(fù)輸入端點(diǎn)的坐標(biāo)"PR2END"程序說(shuō)明：程序先計(jì)算測(cè)站至零方向的邊長(zhǎng)與方位角，再分別計(jì)算測(cè)站至任意端點(diǎn)的邊長(zhǎng)與方位角，并由此計(jì)算出零方向與任意端點(diǎn)方向的水平夾角。x0時(shí)，停止程序運(yùn)行。2-1PR2O01，2，3，4點(diǎn)的邊長(zhǎng)、方位角與水平夾角，表中的紅色數(shù)字為計(jì)算結(jié)果。表2-1使用程序PR2計(jì)算極坐標(biāo)法放樣數(shù)據(jù)點(diǎn)號(hào)x(m)y(m)起訖點(diǎn)號(hào)DOn(m)On(°′″)0On(°′″)測(cè)站O4854.2002760.2790點(diǎn)4826.1072844.794O→089.0621082312.9814842.6162727.851O→134.4352502031.731415718.7524819.5522718.204O→254.5052303144.991220832.0134838.8452672.077O→389.5292600727.861514414.8944861.9092681.723O→478.9332753616.951671303.97ListX[1]，ListY[1]2ListX，ListY，List中，請(qǐng)按(REG)ListYList中的六十進(jìn)制的角值，應(yīng)將光標(biāo)移動(dòng)到其中的某個(gè)單元，按鍵。3、程序PR3，經(jīng)緯儀視距法測(cè)圖計(jì)算記錄"STADIAMAPPINGPR3"顯示程序標(biāo)題Deg:FreqOn:Fix3 基本設(shè)置提示輸入測(cè)站點(diǎn)高程 提示輸入測(cè)站儀器高提示輸入經(jīng)緯儀豎盤(pán)指標(biāo)差I(lǐng)fN=0:ThenClrStat:IfEnd清除統(tǒng)計(jì)串列ListX,ListFreqLbl0:Norm1"UPPER(mm)="?A:"LOWER(mm)="?B提示輸入碎部點(diǎn)觀測(cè)上、下絲讀數(shù)Fix3:"VER(Deg)="?V提示輸入碎部點(diǎn)方向豎盤(pán)讀數(shù)"HOR(Deg)="?L 提示輸入碎部點(diǎn)方向水平盤(pán)讀0.1Abs(B-A)cos(90-V+X)2D計(jì)算測(cè)站至碎部點(diǎn)的水平距離H+Dtan(90-V+X)+I-(A+B)÷2000G計(jì)算碎部點(diǎn)的高程"DI(m)=":D 顯示水平距離"HI(m)=":G 顯示水平距離Norm1 設(shè)置正常顯示格式"STOList,YES(1),NO(Else)."?Y是否將碎部點(diǎn)展點(diǎn)數(shù)據(jù)存儲(chǔ)到串列IfY=1:ThenN+1N:LListX[N]DListY[N]:GListFreq[N]"ListNUMBER=":N 顯示已存儲(chǔ)到串列中的碎部點(diǎn)數(shù)據(jù)個(gè)數(shù)Else"ListNUMBER=":N 顯示已存儲(chǔ)到串列中的碎部點(diǎn)數(shù)據(jù)個(gè)數(shù)"GOON(1),Stop(Else)."?Z按 鍵繼續(xù)，按其余鍵停止程序運(yùn)行IfZ=1:ThenGoto0:IfEnd判斷是否繼續(xù)運(yùn)行程序"PR3END"ListListYListN。程序運(yùn)行前，如要清除串列單元的數(shù)據(jù)，應(yīng)在執(zhí)COMP0NN清零。已知數(shù)據(jù)相同，則可分別按鍵使用上次輸入值。每觀測(cè)一個(gè)碎部點(diǎn)都要求用戶依次輸入上絲讀數(shù)、下絲讀數(shù)、豎盤(pán)讀數(shù)與水平盤(pán)讀數(shù)四個(gè)觀測(cè)數(shù)據(jù)，程序依此計(jì)算出測(cè)站至碎部點(diǎn)的水平距離與碎部點(diǎn)高程，水平盤(pán)讀數(shù)不參與計(jì)算。STOList,YES(1),NO(Else).?，表示是否將該碎部點(diǎn)的水平盤(pán)讀數(shù)、水平距離與高程值存儲(chǔ)到統(tǒng)計(jì)串列中，按鍵為存儲(chǔ)，按其余鍵后按鍵為不存儲(chǔ)。GOON(1),Stop(Else).?鍵為繼續(xù)，按其余鍵后按鍵為停止程序運(yùn)行。停止程序運(yùn)行后，可以按鍵進(jìn)入REG模式查看存儲(chǔ)在統(tǒng)計(jì)串列中的碎部點(diǎn)展點(diǎn)數(shù)據(jù)，--PAGE6-ListXYFreq存儲(chǔ)的是碎部點(diǎn)的高程值。3-1103-1PR3記錄計(jì)算經(jīng)緯儀視距法測(cè)圖數(shù)據(jù)=5.553m，測(cè)站儀器高=1.42m，豎盤(pán)指標(biāo)差=-0°15′序上絲讀數(shù)(mm)下絲讀數(shù)(mm)豎盤(pán)讀數(shù)(°′)水平盤(pán)讀數(shù)(°′)水平距離(m)高程(m)15001485910084598.4533.8322800169890223584589.7904.758316001844892901624.3995.36541700205589293593635.4995.26151100175390193585665.2944.9016170025108958166.0580.9994.56271400231289211780191.1965.75481200204189291781284.0985.74491300174089301783743.9995.64510130013429102187134.1985.5584、程序PR4，單一閉附合圖根水準(zhǔn)路線近似平差計(jì)算"CLOSEDOrCONNECTINGMAPPINGLEVELPR4"顯示程序標(biāo)題"LEVELTYPE,PLATE(1)="?P輸入水準(zhǔn)路線類(lèi)型數(shù)字，1代表平坦，其余數(shù)代表山地H(m)="?A輸入起點(diǎn)已知高程"ENDH(m)="?B輸入終點(diǎn)已知高程閉合水準(zhǔn)路線時(shí)為起點(diǎn)已知高程 基本設(shè)置"HEIGHTNUM="?N輸入測(cè)段高差數(shù)NDimZ 定義額外數(shù)組變量用于存儲(chǔ)平差后點(diǎn)的高程0X:0Y累加高差,路線長(zhǎng)變量清零For1IToN"n=":I 顯示測(cè)段號(hào)"h(m)="?ListX[I]輸入測(cè)段高差"L(km)Orn="?ListY[I]輸入測(cè)段路線長(zhǎng)或測(cè)站數(shù)ListX[I]+XX:ListY[I]+YY累加路線長(zhǎng)與高差A(yù)+X-BF 計(jì)算高差閉合差I(lǐng)fP=1:Then0.04 (Y)W:Else0.012 (Y)W:IfEnd圖根水準(zhǔn)高差閉合差限差"HCLOSEERROR(MM)=":1000F 顯示高差閉合差Σ(L)kmOrΣ(n)="Y 顯示總路線長(zhǎng)或總測(cè)站數(shù)IfAbs(F)>W:Then"CLOSEERROROVER.":IfEnd顯示高差閉合差超限-F÷YV 計(jì)算每km或每站高差改正數(shù)For1IToNListX[I]+VListY[I]ListFreq[I]計(jì)算改正后的測(cè)段高差I(lǐng)fI=1:ThenA+ListFreq[I]Z[I]:ElseZ[I-1]+ListFreq[I]Z[I]:IfEnd"n=":I 顯示測(cè)段號(hào)"hADJUST(m)=":ListFreq[I] 顯示改正后的測(cè)段高"HnADJUST(m)=":Z[I] 顯示改正后點(diǎn)的高程"CLOSETEST(mm)=":1000(Z[I]-B) 顯示高程檢核結(jié)果"PR4END"程序說(shuō)明程序先提示用戶輸入水準(zhǔn)路線類(lèi)型按 鍵為平坦路線要求其后輸入以km為單位的水準(zhǔn)測(cè)段路線長(zhǎng)按其余數(shù)字鍵與 鍵為山地水準(zhǔn)路線要求其后輸入水準(zhǔn)測(cè)段測(cè)站數(shù)。提示用戶輸入起點(diǎn)已知高程與終點(diǎn)已知高程，當(dāng)為閉合水準(zhǔn)路線時(shí)，輸入的終點(diǎn)高程應(yīng)等于起點(diǎn)高程。提示用戶輸入水準(zhǔn)測(cè)段數(shù)，按從起點(diǎn)到終點(diǎn)方向的順序分別輸入測(cè)段高差與路線長(zhǎng)(或測(cè)站數(shù))，并分別顯示以mm為單位的高差閉合差、以km為單位的總路線長(zhǎng)(或總測(cè)站數(shù))。當(dāng)閉合差超限時(shí)，程序給出提示"CLOSEERROROVER."，并繼續(xù)計(jì)算。程序按從起點(diǎn)到終點(diǎn)方向的順序分別顯示測(cè)段高差平差值、未知點(diǎn)高程值及檢核計(jì)算結(jié)果。4-1BM-ABM-B為已知高程的水準(zhǔn)點(diǎn)，圖中箭頭表示水準(zhǔn)測(cè)量前進(jìn)方向，路線上方的數(shù)字為測(cè)得的兩點(diǎn)間的高差(m位(km為單位123點(diǎn)的高程。圖4-1附合水準(zhǔn)路線略圖表4-1使用程序PR4計(jì)算圖根水準(zhǔn)測(cè)量的成果點(diǎn)名路線長(zhǎng)Li(km)觀測(cè)高差hi(m)改正數(shù)Vi(m)改正后高差h?(m)i高程H(m)BM-A1.6+2.331-0.0082.32345.286147.6092.1+2.813-0.0112.802250.4111.7-2.244-0.008-2.252348.1592.0+1.430-0.010+1.420BM-B49.5797.4-0.0375、程序PR5，單一閉附合圖根導(dǎo)線近似平差計(jì)算"CLOSEDOrCONNECTINGMAPPINGPR5"顯示程序標(biāo)題Norm1"CLOSED(0)OrCONNECTING(1)="?Z導(dǎo)線類(lèi)型，1為閉合導(dǎo)線，0為附合導(dǎo)線"UNKNOWNPOINTn="?P未知導(dǎo)線點(diǎn)數(shù)Deg:Fix3 設(shè)置十進(jìn)制度為單位輸入已知點(diǎn)A的坐標(biāo)輸入已知點(diǎn)B的坐標(biāo)計(jì)算A→B方向的方位角IfJ<0:ThenJ+360R:ElseJR:IfEnd"DISTAB(m)=":I A→B方向的水平距離"BEARAB(DMS)=":R A→B方向的方位角IfZ=1:Then"XC(m)="?E:"YC(m)="?F若為附合導(dǎo)線，則繼續(xù)輸入已知點(diǎn)C的坐標(biāo)輸入已知點(diǎn)D的坐標(biāo)計(jì)算C→D方向的方位角IfJ<0:ThenJ+360S:ElseJS:IfEnd"DISTCD(m)=":I C→D方向的水平距離--PAGE10-"BEARCD(DMS)=":S C→D方向的方位角ElseIfR>180:ThenR-180S:ElseR+180S:IfEnd若為閉合導(dǎo)線，則計(jì)算A→B方向的反方位角"BEARBA(DMS)=":S B→A方向的方位角 計(jì)算閉合導(dǎo)線的反方位角ClrStat:reqOn 清除統(tǒng)計(jì)串列ListXLit,Liste打開(kāi)頻度串列0M累加邊長(zhǎng)和變量清零1IP+1Norm1:"n=":I 顯示當(dāng)前輸入的觀測(cè)數(shù)據(jù)計(jì)數(shù)"ANGLE(Deg)="?ListX[I]輸入水平角觀測(cè)值"DIST(m)="?ListY[I]輸入水平距離觀測(cè)值ListY[I]+MM累加水平距離和IfI=1:ThenR+ListX[I]L:ElseListFreq[I-1]+ListX[I]L:IfEnd推算導(dǎo)線邊方位角IfL>180:ThenL-180L:ElseL+180L:IfEndIfL>360:ThenL-360L:IfEnd判斷方位角是否大于360LListFreq[I]存儲(chǔ)導(dǎo)線邊方位角"LASTANGLE(Deg)="?ListX[P+2]輸入最后一個(gè)水平角ListFreq[P+1]+ListX[P+2]LIfL>180:ThenL-180L:ElseL+180L:IfEndIfL>360:ThenL-360L:IfEnd判斷方位角是否大于360LListFreq[P+2]存儲(chǔ)最后一個(gè)方位角3600(L-S)U 以秒為單位的方位角閉合差Fix1:"ANGLECLOSEERROR(S)=":U 顯示方位角閉合差60 (P+2)W 以秒為單位的方位角閉合差限差I(lǐng)fAbs(U)>W:Then"ANGLECLOSEERROROVER.":IfEnd顯示角度閉合差超限-U÷(P+2)V 計(jì)算角度改正數(shù)2DimZ:0Z[1]:0Z[2]定義額外變量數(shù)組用于存儲(chǔ)坐標(biāo)增量累加和1IP+1 分配角度改正數(shù)并計(jì)算導(dǎo)線邊的坐標(biāo)增量ListX[I]+V÷3600ListX[I]計(jì)算并存儲(chǔ)改正后的角度IfI=1:ThenR+ListX[I]L:ElseListFreq[I-1]+ListX[I]L:IfEnd推算導(dǎo)線邊方位角IfL>180:ThenL-180L:ElseL+180L:IfEndIfL>360:ThenL-360L:IfEnd判斷方位角是否大于360LListFreq[I]存儲(chǔ)角度改正后的導(dǎo)線邊方位角ListY[I]cos(L)X:ListY[I]sin(L)Y計(jì)算導(dǎo)線邊坐標(biāo)增量X+Z[1]Z[1]:Y+Z[2]Z[2]累加導(dǎo)線邊坐標(biāo)增量ListFreq[P+1]+ListX[P+2]+V÷3600L計(jì)算最后一條導(dǎo)線邊的方位角IfL>180:ThenL-180L:ElseL+180L:IfEndIfL>360:ThenL-360L:IfEnd判斷方位角是否大于360LListFreq[P+2]存儲(chǔ)最后一條導(dǎo)線邊的方位角3600(L-S)U 以秒為單位的方位角閉合差檢核計(jì)算"CHECKANGLECLOSEERROR(S)=":U 顯示方位角閉合差檢核結(jié)果IfZ=0:ThenZ[1]U:Z[2]V 計(jì)算閉合導(dǎo)線的坐標(biāo)增量閉合ElseC+Z[1]-EU:D+Z[2]-F計(jì)算附合導(dǎo)線的坐標(biāo)增量閉合M÷(U2+V2)K 計(jì)算導(dǎo)線全長(zhǎng)相對(duì)閉合差FixX(m)=":U X坐標(biāo)增量閉合差Y(m)=":V Y坐標(biāo)增量閉合差--PAGE8-Fix0:"RELATCLOSEERROR=":K 顯示導(dǎo)線全長(zhǎng)相對(duì)閉合差-U÷MV 計(jì)算坐標(biāo)增量閉合差每米改正數(shù)IfK<2000:Then"RELATCLOSEERROROVER.":IfEnd顯示全長(zhǎng)相對(duì)閉合差超限2(P+1)DimZ 定義額外變量數(shù)組用于存儲(chǔ)未知點(diǎn)的坐標(biāo)1IP+1 分配坐標(biāo)閉合差并計(jì)算未知點(diǎn)的坐ListY[I]cos(ListFreq[I])X:ListY[I]sin(ListFreq[I])Y計(jì)算導(dǎo)線邊坐標(biāo)增量X+UListY[I]X:Y+VListY[I]Y計(jì)算改正后的導(dǎo)線邊坐標(biāo)增量IfI=1:ThenC+XZ[2I-1]:D+YZ[2I]ElseZ[2(I-1)-1]+XZ[2I-1]:Z[2(I-1)]+YZ[2I]:IfEndNorm1:"POINTn=":I 顯示未知點(diǎn)號(hào)Fix3:"XP(m)=":Z[2I-1] X坐標(biāo)"YP(m)=":Z[2I] Y坐標(biāo)IfZ=0:ThenZ[2(P+1)-1]-CX:Z[2(P+1)]-DY計(jì)算閉合導(dǎo)線的坐標(biāo)檢核結(jié)果ElseZ[2(P+1)-1]-EX:Z[2(P+1)]-F計(jì)算附合導(dǎo)線的坐標(biāo)檢核結(jié)果"CHECKX(m)=":X 顯示X坐標(biāo)檢核計(jì)算結(jié)果"CHECKY(m)=":Y 顯示Y坐標(biāo)檢核計(jì)算結(jié)果ClrStat:reqOff 清除統(tǒng)計(jì)串列ListXLit,Liste關(guān)閉頻度串1IP分配坐標(biāo)閉合差并計(jì)算未知點(diǎn)的坐標(biāo)Z[2I-1]ListX[I]:Z[2I]ListY[I]將未知點(diǎn)的坐標(biāo)存儲(chǔ)到統(tǒng)計(jì)串列ListX,ListY中"PR5END"程序說(shuō)明：程序先提示用戶輸入導(dǎo)線類(lèi)型與未知點(diǎn)總數(shù)。當(dāng)提示導(dǎo)線類(lèi)型時(shí)，按鍵為A，BA→BB→A的方位角；A，B，C，DB，C→D的方位角。PP+2A→B→1→2→3→……→B→AA→B→1→2→3→……→C→D輸入正數(shù)；水平角位于方位角推算路線右邊時(shí)，角度應(yīng)輸入負(fù)數(shù)。合差后的檢核結(jié)果、X，Y結(jié)果等。A，B，C，D，E，F(xiàn)，G，HListXListYListXListX中，YListYListFreq。程序運(yùn)行完成后，用戶可以按(REG)鍵進(jìn)入雙變量統(tǒng)計(jì)回歸模式查看未知點(diǎn)的坐標(biāo)成果。[5-1]5-16角，其中的兩個(gè)為右角；5PR55-1。表5-1使用程序PR5計(jì)算圖根附合導(dǎo)線的成果角度閉合差-36″點(diǎn)號(hào)x?(m)y?(m)X坐標(biāo)增量閉合差-0.149m12299.8241303.798Y坐標(biāo)增量閉合差0.149m22186.2821383.969全長(zhǎng)相對(duì)閉合差1/352032192.4501556.39642179.7401655.637--PAGE9-圖5-1鋼尺量距圖根附合導(dǎo)線[5-2]5-25角；4PR55-2。表5-2使用程序PR5計(jì)算圖根閉合導(dǎo)線的成果角度閉合差50″點(diǎn)號(hào)x?(m)y?(m)X坐標(biāo)增量閉合差0.057m1445.195301.329Y坐標(biāo)增量閉合差-0.070m2493.090365.640全長(zhǎng)相對(duì)閉合差1/43363569.667261.442圖5-2光電測(cè)距圖根閉合導(dǎo)線6、程序PR6，圖根支導(dǎo)線坐標(biāo)計(jì)算"MAPPINGOPENPR6"顯示程序標(biāo)題Norm1:"OPENPOINTNUM="?P支導(dǎo)線點(diǎn)數(shù)Fix3:Deg:ClrStat:FreqOn基本設(shè)置輸入已知點(diǎn)A的坐標(biāo)輸入已知點(diǎn)B的坐標(biāo)計(jì)算A→B方向的方位角IfJ<0:ThenJ+360R:ElseJR:IfEnd"DISTAB(m)=":I A→B方向的水平距離"BEARAB(DMS)=":R A→B方向的方位角2PDimZ 定義額外數(shù)組變量用于存儲(chǔ)觀測(cè)水平角與水平距離For1IToPNorm1:"POINTn=":I 顯示支導(dǎo)線點(diǎn)號(hào)"ANGLE(Deg)="?Z[2I-1] 輸入水平角觀測(cè)值--PAGE10-"DIST(m)="?Z[2I] 輸入水平距離觀測(cè)值IfI=1:ThenR+Z[2I-1]L:ElseListFreq[I-1]+Z[2I-1]L:IfEnd推算導(dǎo)線邊方位角IfL>180:ThenL-180L:ElseL+180L:IfEndIfL>360:ThenL-360L:IfEnd判斷方位角是否大于360LListFreq[I]存儲(chǔ)導(dǎo)線邊方位角Z[2I]cos(L)X:Z[2I]sin(L)Y計(jì)算導(dǎo)線邊的坐標(biāo)增量IfI=1:ThenC+XListX[I]:D+YListY[I]ElseListX[I-1]+XListX[I]:ListY[I-1]+YListY[I]:IfEnd推算支導(dǎo)線點(diǎn)的坐標(biāo)Fix3:"XP(m)=":ListX[I] X坐標(biāo)"YP(m)=":ListY[I] Y坐標(biāo)"PR6END"A→B的邊長(zhǎng)與方位角。XListX中，YListYListFreq中。程序運(yùn)行完成后，用戶可以按(REG)鍵進(jìn)入雙變量統(tǒng)計(jì)回歸模式查看支導(dǎo)線點(diǎn)的坐標(biāo)與方位角成果。[6-1]5-1B兩個(gè)已知點(diǎn)，五個(gè)未知點(diǎn)，觀測(cè)了5個(gè)水平角，其中的2個(gè)為右角；5條水平距離，運(yùn)行程序PR6的計(jì)算結(jié)果列于表6-1。7、程序PR7，測(cè)角前方交會(huì)坐標(biāo)計(jì)算7-1A、BP觀測(cè)水平角,和檢查角CAP定待定點(diǎn)P的坐標(biāo)。為保證交會(huì)定點(diǎn)的精度，在選定P點(diǎn)時(shí)，使交會(huì)角 位于30°~120°之間，最好近于90°。

表6-1使用程序PR6計(jì)算圖根支導(dǎo)線的成果點(diǎn)號(hào)x?(m)y?(m)12299.7811303.85022186.2161384.05532192.3641556.51642179.6461655.7795 2166.641 1757.435xAcotxBcot(yByA)P cotcot ycoty

cot(x

x)

(7-1)A B A BP cotcot 應(yīng)用(7-1)A、B、P成的旋轉(zhuǎn)方向?yàn)槟鏁r(shí)針?lè)较蚯遗c實(shí)際情況相符。按(7-1)計(jì)算出的P點(diǎn)坐標(biāo)，可以通過(guò)反算坐標(biāo)方位角求出檢查角，與其觀測(cè)值之差為：一般要求應(yīng)小于給定的限差。

(7-2)

圖7-1前方交會(huì)INTERSECTIONPR7"顯示程序標(biāo)題Fix3:Deg 基本設(shè)置輸入已知點(diǎn)A的坐標(biāo)輸入已知點(diǎn)B的坐標(biāo)Norm1:"CPOINTYES(1)OrNO(0)"?Z是否有檢查點(diǎn)Fix3IfZ=1:Then"XC(m)="?E:"YC(m)="?F:IfEnd輸入已知點(diǎn)C的坐標(biāo)"∠A(Deg)="?O 輸入 角--PAGE11-"∠B(Deg)="?P 角IfZ=1:Then∠C(Deg)="?Q:IfEnd C點(diǎn)坐標(biāo)時(shí)輸入角1÷tan(O)M:1÷tan(P)N計(jì)算角度的余切(AN+CM+D-B)÷(M+N)X使用式(7-2)計(jì)算交會(huì)點(diǎn)的坐標(biāo)(BN+DM+A-C)÷(M+N)Y"XP(m)=":X X坐標(biāo)"YP(m)=":Y Y坐標(biāo)IfZ=1:ThenPol(X-A,Y-B):ClsIfJ<0:ThenJ+360R:ElseJR:IfEnd計(jì)算A→P方向的方位角IfJ<0:ThenJ+360S:ElseJS:IfEnd計(jì)算A→C方向的方位角V<0:ThenV+360計(jì)算∠CAP"CAL-CHECKANGLE(DMS)=":V 顯示計(jì)算出的檢查角Q-VW 計(jì)算檢查角觀測(cè)值與計(jì)算值之差"CHECKANGLEERROR(DMS)=":W 顯示檢查角觀測(cè)值與計(jì)算值之差Else"NOCHECKANGLE.":IfEnd顯示沒(méi)有觀測(cè)檢查角"PR7END"CC點(diǎn)的坐標(biāo)。觀測(cè)數(shù)據(jù)輸入時(shí)要求必須輸入 ，角若觀測(cè)了檢查角則應(yīng)輸入檢查角否則不需要輸入檢查角。P反算出的檢查角及其與觀測(cè)角的差。PX，Y7-17-1的計(jì)算結(jié)果，紅色數(shù)據(jù)為程序計(jì)算結(jié)果。點(diǎn)名x坐標(biāo)(m) y點(diǎn)名x坐標(biāo)(m) y坐標(biāo)(m) 角名水平角A781.2271323.38553°33′42″B769.8621428.98950°34′55″C866.2941306.11754°03′08″P843.5341380.64154°03′20″差0°00′11.92″8、程序PR8，測(cè)角后方交會(huì)坐標(biāo)計(jì)算8-1P,和檢查角PABCPPPP點(diǎn)位于危險(xiǎn)圓上。ABCABC的三個(gè)內(nèi)角∠AB，按下列公式計(jì)算出系數(shù)。圖8-1測(cè)角后方交會(huì)--PAGE12-P 1 tantanAA cotAcot tantanAP 1

tantanB

(8-1)B cotBcot

tantanBP 1

tantanCC cotCcot

tantanC則P點(diǎn)坐標(biāo)的計(jì)算公式為：PxAxBxCPPAPB

P CC

(8-2)yAyByCP PPP A B C 主程序：PR8"ANGLERESECTIONPR8"顯示程序標(biāo)題Fix3:Deg:ClrStat:FreqOn基本設(shè)置"XA(m)="?A:AListX[1]輸入已知點(diǎn)A的X坐標(biāo)"YA(m)="?B:BListY[1]輸入已知點(diǎn)A的Y坐標(biāo)"XB(m)="?C:CListX[2]輸入已知點(diǎn)B的X坐標(biāo)"YB(m)="?D:DListY[2]輸入已知點(diǎn)B的Y坐標(biāo)"XC(m)="?E:EListX[3]輸入已知點(diǎn)C的X坐標(biāo)"YC(m)="?F:FListY[3]輸入已知點(diǎn)C的Y坐標(biāo)Norm1:"DPOINTYES(1)OrNO(0)"?Z是否有檢查點(diǎn)Fix3IfZ=1:Then"XD(m)="?G:GListX[4]輸入已知點(diǎn)D的X坐標(biāo)"YD(m)="?H:HListY[4]:IfEnd輸入已知點(diǎn)D的Y坐標(biāo)14DimZ 定義額外數(shù)組變量"ANGLEA(Deg)="?Z[11]:Z[11]ListFreq[1]輸入角A"ANGLEB(Deg)="?Z[12]:Z[12]ListFreq[2]輸入角B"ANGLEC(Deg)="?Z[13]:Z[13]ListFreq[3]輸入角CIfZ=1:Then"CHECKANGLED(Deg)="?Z[14]:Z[14]ListFreq[4]:IfEnd有D點(diǎn)坐標(biāo)輸入角O:2P:3Q為調(diào)子程序計(jì)算∠A準(zhǔn)備Prog"ANGLE"調(diào)子程序計(jì)算∠A"∠A(DMS)=":V 顯示∠A1÷(tan(V)-1-tan(ListFreq[1])-1)Z[1] O:3P:1Q為調(diào)子程序計(jì)算∠B準(zhǔn)備Prog"ANGLE"調(diào)子程序計(jì)算∠B"∠B(DMS)=":V 顯示∠B1÷(tan(V)-1-tan(ListFreq[2])-1)Z[2] O:1P:2Q為調(diào)子程序計(jì)算∠C準(zhǔn)備Prog"ANGLE"調(diào)子程序計(jì)算∠C"∠C(DMS)=":V 顯示∠C1÷(tan(V)-1-tan(ListFreq[3])-1)Z[3] (Z[1]ListX[1]+Z[2]ListX[2]+Z[3]ListX[3])÷(Z[1]+Z[2]+Z[3])X 計(jì)算交會(huì)點(diǎn)坐標(biāo)(Z[1]ListY[1]+Z[2]ListY[2]+Z[3]ListY[3])÷(Z[1]+Z[2]+Z[3])Y--PAGE14-"XP(m)=":X 顯示交會(huì)點(diǎn)坐標(biāo)"YP(m)=":YIfZ=1:ThenXListX[5]:YListY[5]如果觀測(cè)了檢查角5O:3P:4Q為調(diào)子程序計(jì)算檢查角準(zhǔn)備Prog"ANGLE"調(diào)子程序計(jì)算檢查角"∠CHECK(DMS)=":V 顯示計(jì)算的檢查角ListK 計(jì)算檢查角觀測(cè)值與計(jì)算值之差"CHECKANGLEERROR(DMS)=":K 顯示檢查角差Else"NOCHECKANGLE.":IfEnd顯示沒(méi)有觀測(cè)檢查角"PR8END"ANGLEPol(ListX[P]-ListX[O],ListY[P]-ListY[O]):ClsIfJ<0:ThenJ+360R:ElseJR:IfEnd計(jì)算O→P方向的方位角Pol(ListX[Q]-ListX[O],ListY[Q]-ListY[O]):ClsIfJ<0:ThenJ+360S:ElseJS:IfEnd計(jì)算O→Q方向的方位角V 計(jì)算方位角差I(lǐng)fV<0:ThenV+360 返回主程序程序說(shuō)明：程序顯示完計(jì)算標(biāo)題后，先提示用戶輸入A，B，C點(diǎn)的已知坐標(biāo)，然后提示是否觀測(cè)了檢查角，若觀測(cè)了檢查角，則提示輸入D點(diǎn)的坐標(biāo)，否則不提示輸入D點(diǎn)的坐標(biāo)。觀測(cè)數(shù)據(jù)輸入時(shí)要求必須輸入 ，，角若觀測(cè)了檢查角則應(yīng)輸入檢查角否則不需要輸入檢查角。完成已知與觀測(cè)數(shù)據(jù)輸入后，屏幕顯示計(jì)算出的A，B，C及后方交會(huì)點(diǎn)P的坐標(biāo)，若有檢查角，還將顯示根據(jù)坐標(biāo)反算出的檢查角及其與觀測(cè)角的差。PX，Y8-18-1的計(jì)算結(jié)果，紅色數(shù)據(jù)為程序計(jì)算結(jié)果。表8-1使用程序PR8計(jì)算測(cè)角后方交會(huì)點(diǎn)的成果點(diǎn)名x坐標(biāo)(m)y坐標(biāo)(m)角名水平角頂角計(jì)算值A(chǔ)2502.9325044.00978°12′17″∠A41°01′26.93″B3229.4165185.604161°05′58″∠B87°47′23.94″C3325.0064569.476120°41′45″∠C51°11′09.13″D2724.7014416.40071°10′24″P2896.7634906.24171°10′23.77″差0°00′0.23″9、程序PR9，測(cè)邊后方交會(huì)點(diǎn)坐標(biāo)計(jì)算8-1RARBRCPABCabc。由余弦定理可以求出下列三個(gè)偏角R2c2R2

R2a2R2

R2b2R2A B cos1A B， cos1B C， cos1C AA B 2RAc

2RBa

2RCb 則A→P，B→P，C→P方向的方位角為APABA，BPBCB，CPCAC式中的“±”應(yīng)按圖9-1的四種情形區(qū)別對(duì)待，其中圖9-1(a)為全取“－”號(hào)，圖9-1(b)依次9-c9-1(依次分別取，B，CP點(diǎn)坐標(biāo)的公式為xPxARAcosAP xxBRBcosBP

xxCCcosCPyyRsin ，y

Rsin

，yyRsin P A

AP P B

BP P C

CP圖9-1測(cè)邊后方交會(huì)的方位角推算的四種情形PP主程序：PR9"SIDERESECTIONPR9"顯示程序標(biāo)題Fix3:Deg:ClrStat:FreqOn基本設(shè)置"XA(m)="?A:AListX[1]輸入已知點(diǎn)A的X坐標(biāo)"YA(m)="?B:BListY[1]輸入已知點(diǎn)A的Y坐標(biāo)"XB(m)="?C:CListX[2]輸入已知點(diǎn)B的X坐標(biāo)"YB(m)="?D:DListY[2]輸入已知點(diǎn)B的Y坐標(biāo)"XC(m)="?E:EListX[3]輸入已知點(diǎn)C的X坐標(biāo)"YC(m)="?F:FListY[3]輸入已知點(diǎn)C的Y坐標(biāo)26DimZ 定義額外數(shù)組變量"SIDEPA(m)="?ListFreq[1] RA"SIDEPB(m)="?ListFreq[2] RB"SIDEPC(m)="?ListFreq[3] RCO:2P為調(diào)子程序計(jì)算邊長(zhǎng)與方位角賦起訖點(diǎn)號(hào)ProgA→B的邊長(zhǎng)與方位角SZ[1]:RZ[2] A→B的邊長(zhǎng)c與方位角O:3P為調(diào)子程序計(jì)算邊長(zhǎng)與方位角賦起訖點(diǎn)號(hào)ProgB→Ca與方位角SZ[3]:RZ[4] B→Ca與方位角O:1P為調(diào)子程序計(jì)算邊長(zhǎng)與方位角賦起訖點(diǎn)號(hào)ProgC→A的邊長(zhǎng)b與方位角SZ[5]:RZ[6] C→A的邊長(zhǎng)b與方位角cos-1((ListFreq[1]2+Z[1]2-ListFreq[2]2)÷(2ListFreq[1]Z[1]))Z[7]計(jì)算Acos-1((ListFreq[2]2+Z[3]2-ListFreq[3]2)÷(2ListFreq[2]Z[3]))Z[8]計(jì)算Bcos-1((ListFreq[3]2+Z[5]2-ListFreq[1]2)÷(2ListFreq[3]Z[5]))Z[9]計(jì)算CO:2P:3Q 為調(diào)子程序計(jì)算∠A準(zhǔn)備ProgANGLE" 調(diào)子程序計(jì)算∠AVZ[17] 存儲(chǔ)∠AO:3P:1Q 為調(diào)子程序計(jì)算∠B準(zhǔn)備Prog"ANGLE"調(diào)子程序計(jì)算∠BVZ[18]存儲(chǔ)∠BO:1P:2Q為調(diào)子程序計(jì)算∠C準(zhǔn)備Prog"ANGLE"調(diào)子程序計(jì)算∠CVZ[19]存儲(chǔ)∠CZ[1]sin(Z[18])Z[11]計(jì)算A點(diǎn)至BC邊長(zhǎng)的垂直距離Z[3]sin(Z[19])Z[12]計(jì)算B點(diǎn)至AC邊長(zhǎng)的垂直距離Z[5]sin(Z[17])Z[13]計(jì)算C點(diǎn)至AB邊長(zhǎng)的垂直距離IfZ[7]<Z[17]AndZ[8]<Z[18]AndZ[9]<Z[19]P點(diǎn)位于△ABC內(nèi)ThenZ[2]-Z[7]Z[14]:Z[4]-Z[8]Z[15]:Z[6]-Z[9]Z[16]:IfEndIfListFreq[1]>Z[11]AndZ[9]>Z[19]P點(diǎn)位于BC邊外ThenZ[2]-Z[7]Z[14]:Z[4]+Z[8]Z[15]:Z[6]-Z[9]Z[16]:IfEndIfListFreq[2]>Z[12]AndZ[7]>Z[17]P點(diǎn)位于AC邊外ThenZ[2]-Z[7]Z[14]:Z[4]-Z[8]Z[15]:Z[6]+Z[9]Z[16]:IfEndIfListFreq[3]>Z[13]AndZ[8]>Z[18]P點(diǎn)位于AB邊外ThenZ[2]+Z[7]Z[14]:Z[4]-Z[8]Z[15]:Z[6]-Z[9]Z[16]:IfEndListX[1]+ListFreq[1]cos(Z[14])Z[21]由A點(diǎn)計(jì)算P點(diǎn)的坐標(biāo)ListY[1]+ListFreq[1]sin(Z[14])Z[22]ListX[2]+ListFreq[2]cos(Z[15])Z[23]由B點(diǎn)計(jì)算P點(diǎn)的坐標(biāo)ListY[2]+ListFreq[2]sin(Z[15])Z[24]ListX[3]+ListFreq[3]cos(Z[16])Z[25]由C點(diǎn)計(jì)算P點(diǎn)的坐標(biāo)ListY[3]+ListFreq[3]sin(Z[16])Z[26](Z[21]+Z[23]+Z[25])÷3X計(jì)算P點(diǎn)坐標(biāo)的平均值(Z[22]+Z[24]+Z[26])÷3Y計(jì)算P點(diǎn)坐標(biāo)的平均值"P(m)=":X P點(diǎn)坐標(biāo)三組坐標(biāo)的平均值"P(m)=":Y"P點(diǎn)第一組坐標(biāo)與平均值的差"P-Y1(m)=":Y-Z[22]"P點(diǎn)第二組坐標(biāo)與平均值的差"P-Y2(m)=":Y-Z[24]"P點(diǎn)第三組坐標(biāo)與平均值的差"P-Y3(m)=":Y-Z[25]"PR9END"子程序：DIST-BEARPol(ListX[P]-ListX[O],ListY[P]-ListY[O]):ClsIS 保存O→P方向的邊長(zhǎng)IfJ<0:ThenJ+360R:ElseJR:IfEnd計(jì)算O→P方向的方位角 返回主程序程序說(shuō)明：A，B，CPA，B，CRARBRC。A，B，CPPXYPY中，X，YU，V中。表9-1為圖8-1的計(jì)算結(jié)果，紅色數(shù)據(jù)為程序計(jì)算結(jié)果。表9-1使用程序PR9計(jì)算圖9-1四種情形的測(cè)邊后方交會(huì)點(diǎn)的成果情形點(diǎn)名x(m)y(m)邊名平距(m)xp(m)yp(m)xmax(m)ymax(m)圖A2502.9325044.009RA417.2322896.7634906.2420.0000.0019-1B3229.4165185.604RB434.3980.0000.000(a)C3325.0064569.476RC544.7970.000-0.001圖A2408.0266132.987RA480.0932787.9386426.5110.0000.0009-1B3134.5106274.581RB378.4110.0000.000(b)C3230.1005658.453RC886.2390.0000.000圖A2408.0267255.919RA996.6173395.0757118.1510.0000.0009-1B3134.5107397.513RB382.0170.0000.000(c)C3230.1006781.385RC375.0040.0000.000圖A2408.0268477.046RA535.8782648.3887998.0980.0000.0009-1B3134.5108618.641RB788.2820.0000.000(d)C3230.1008002.513RC581.7290.0000.00010、程序PR10，54北京系與80西安系高斯投影正算af310-1af計(jì)算而來(lái)。表10-1參考橢球的初始元素值坐標(biāo)系名1954北京坐標(biāo)系1980西安坐標(biāo)系 WGS-84坐標(biāo)系(GPS用)參考橢球名克拉索夫斯基橢球IUGG—1975橢球 IUGG—1979橢球a63782456378140 6378137f1:298.31:298.257 1:298.257223563主程序：PR10"GAUSSPROJECTZSPR10"顯示程序標(biāo)題Norm1:Deg:ClrStat:FreqOn:Norm1基本設(shè)置Lbl0"COORDINATESYSTEM(54)Or(80)="?W坐標(biāo)系數(shù)字，只能輸入數(shù)值54或80IfW54:ThenIfW80:Then"NOTBE54Or80,REINPUT.":Goto0:IfEnd:IfEnd10DimZProg"COEFFICIENT""CENTERL0(Deg)="?I定義額外數(shù)組變量輸入中央子午線經(jīng)度"X+CONSTANT(m)="?Z[6]輸入X坐標(biāo)加常數(shù)"Y+CONSTANT(m)="?Z[7]輸入Y坐標(biāo)加常數(shù)0KLbl1K+1K點(diǎn)計(jì)數(shù)變量清零點(diǎn)計(jì)數(shù)變量計(jì)數(shù)"POINTL(Deg)="?J輸入點(diǎn)的經(jīng)度"POINTB(Deg)="?B輸入點(diǎn)的緯度JListX[2K-1]存點(diǎn)的經(jīng)度到統(tǒng)計(jì)串列BListY[2K-1]存點(diǎn)的緯度到統(tǒng)計(jì)串列J-IL計(jì)算點(diǎn)的經(jīng)度差Prog"ZS"調(diào)正算子程序X+Z[6]ListX[2K]X坐標(biāo)到統(tǒng)計(jì)串列Y+Z[7]ListY[2K]Y坐標(biāo)到統(tǒng)計(jì)串列RListFreq[2K]存子午線收斂角到統(tǒng)計(jì)串列"n=":K顯示點(diǎn)號(hào)Fix3:"X(m)=":ListX[2K]顯示高斯X坐標(biāo)"Y(m)=":ListY[2K]顯示高斯Y坐標(biāo)"r(DMS)=":R顯示子午線收斂角Prog"FS"調(diào)反算子程序檢核計(jì)算"GAUSSPROJECTFSCHECK.提示高斯投影反算檢核結(jié)果Norm1L(S)=":3600(ListX[2K-1]-I-L) 顯示反算檢核的經(jīng)度差B(S)=":3600(ListY[2K-1]-B) 顯示反算檢核的緯度差"COMPNextPOINTYES(1)OrNO(0)="?Z[8]IfZ[8]=1:ThenGoto1:IfEnd"PR10END"計(jì)算子午圈弧長(zhǎng)系數(shù)子程序：COEFFICIENTW=54W=80Z[1]~Z[5]。IfW=54:Then6378245A:298.3-1F54北京坐標(biāo)系橢球參數(shù)Else6378140A:298.257-1F:IfEnd80西安坐標(biāo)系橢球參數(shù)(2F-F2)E第一偏心率(1÷(1-F)2-1)Z 第二偏心率A(1-F)N 橢球短半徑A2÷NC 極點(diǎn)子午曲率半徑N2÷AD赤道子午曲率半徑D(1+34E2+4564E^(4)+175256E^(6)+1102516384E^(8))p÷180Z[1]D(34E2+4564E^(4)+175256E^(6)+1102516384E^(8))Z[2]D(1532E^(4)+175384E^(6)+36758192E^(8))Z[3]D(3596E^(6)+7352048E^(8))Z[4]D×3151024E^(8)Z[5]計(jì)算子午圈弧長(zhǎng)公式中的5個(gè)系數(shù)值 返回主程序高斯投影正算子程序：ZS入口參數(shù)為子午圈弧長(zhǎng)系數(shù)值Z[1]~Z[5]系數(shù)值、經(jīng)度差L與緯度B，輸出結(jié)果為x,y坐標(biāo)。Sin(B)P:cos(B)Q:Zcos(B)U:tan(B)T緯度的函數(shù)Z[1]B-(Z[2]P+Z[3]P^(3)+Z[4]P^(5)+Z[5]P^(7))QX子午圈弧長(zhǎng)(1+U2)N 卯酉圈曲率半徑QLp÷180M緯度與經(jīng)度差的函數(shù)X+NT(0.5M2+124(5-T2+9U2+4U^(4))M^(4)+1720(61-58T2+T^(4)+270U2-330U2T2)M^(6))Xx坐標(biāo)N(M+16(1-T2+U2)M^(3)+1120(5-18T2+T^(4)+14U2-58U2T2)M^(5))Yy坐標(biāo)T÷p×(180M+60(1+3U2+2U^(4))M^(3)+12(2-T2)M^(5))R子午線收斂角r 返回主程序高斯投影反算子程序：FS入口參數(shù)子午圈弧長(zhǎng)系數(shù)值Z[1]~Z[5]系數(shù)值、x，y坐標(biāo)，輸出計(jì)算結(jié)果為經(jīng)度差L與緯度B。X÷Z[1]B 底點(diǎn)緯度初始值迭代計(jì)算底點(diǎn)緯度BH:sin(H)P:cos(H)Q底點(diǎn)緯度的函數(shù)-(Z[2]P+Z[3]P^(3)+Z[4]P^(5)+Z[5]P^(7))QG(X-G)÷Z[1]BAbsB-H)O 底點(diǎn)緯度迭代誤差LpWhileO>1×10-12 底點(diǎn)緯度迭代誤差Zcos(B)U:tan(B)T:cos(B)Q 底點(diǎn)緯度的函數(shù)(1+U2)V:YV÷CNB-(1+U2)T÷p×(90N2-7.5(5+3T2+U2-9U2T2)N^(4)+0.25(61+90T2+45T^(4))N^(6))B 緯度--PAGE18-(pQ)-1(180N-30(1+2T2+U2)N^(3)+1.5(5+28T2+24T^(4))N^(5))L 經(jīng)度 返回主程序程序說(shuō)明：主程序需要調(diào)用COEFFICIENT，ZS，F(xiàn)S三個(gè)子程序，其中反算子程序的作用為檢核正算計(jì)算結(jié)果的正確性。5480IUGG1975X，Y0響應(yīng)。XYfx-5800P×10-6″，緯度差<-5×10-7″。計(jì)算案例數(shù)據(jù)見(jiàn)表10-2所示。PtPOINT)r10-2PR10進(jìn)行高斯投影正算成果案例1954北京坐標(biāo)系，中央子午線經(jīng)度L0=111°,x坐標(biāo)加常數(shù)=0m,y坐標(biāo)加常數(shù)=500000m序 L(°′″) B(°′″) X(m) Y(m) r(°′″) 經(jīng)度檢核L(″)

緯度檢核差B(″)1 0513.685 223849.526 2507005.772 714567.798 04814.28 2.88×10-6 -5.9292×10-72 0455.002 223154.967 2494239.989 714212.261 04753.17 2.861136×10-6 -5.886×10-73 0320.216 223839.704 2506658.355 711330.656 04730.22 2.683584×10-6 -5.5872×10-74 0242.902 223803.524 2505530.213 710280.051 04714.64 2.62224×10-6 -5.4756×10-75 0229.798 223654.359 2503396.508 709934.887 04707.32 2.602764×10-6 -5.436×10-76 0402.832 223652.280 2503369.195 712593.951 04743.06 2.759904×10-6 -5.7168×10-77 0343.973 223809.061 2505724.695 712022.373 04738.36 2.724912×10-6 -5.6556×10-78 0542.003 223800.707 2505514.694 715397.885 04823.55 2.932452×10-6 -6.0228×10-7完成計(jì)算后，用戶輸入點(diǎn)的經(jīng)度、緯度分別存儲(chǔ)在ListX與ListY的奇數(shù)行，如1行、3行、5行等……，正算結(jié)果高斯X，Y坐標(biāo)與子午線收斂角r分別存儲(chǔ)在ListX，ListY，List的偶數(shù)行如2行3行5行等……用戶可以按 (REG)鍵進(jìn)入雙變量統(tǒng)計(jì)與回歸模式查看11、程序PR11，54北京系與80西安系高斯投影反算主程序：PR11"GAUSSPROJECTFSPR11"顯示程序標(biāo)題基本設(shè)置Lbl0"COORDINATESYSTEM(54)Or(80)="?W坐標(biāo)系數(shù)字，只能輸入數(shù)值54或80IfW54:ThenIfW80:Then"NOTBE54Or80,REINPUT.":Goto0:IfEnd:IfEnd14DimZ 定義額外數(shù)組變量Prog"COEFFICIENT""CENTERL0(Deg)="?I 輸入中央子午線經(jīng)度"X+CONSTANT(m)="?Z[6]輸入X坐標(biāo)加常數(shù)"Y+CONSTANT(m)="?Z[7]輸入Y坐標(biāo)加常數(shù)0K 點(diǎn)計(jì)數(shù)變量清零Lbl1K+1K 點(diǎn)計(jì)數(shù)變量計(jì)數(shù)"POINTX(m)="?X輸入點(diǎn)的X坐標(biāo)"POINTY(m)="?Y輸入點(diǎn)的Y坐標(biāo)ListX[2K-1]存點(diǎn)的X坐標(biāo)到統(tǒng)計(jì)串列ListY[2K-1] Y坐標(biāo)到統(tǒng)計(jì)串列X-Z[6]X XX坐標(biāo)--PAGE20-Y YY坐標(biāo)Prog"FS"調(diào)反算子程序L+IListX[2K]存點(diǎn)的經(jīng)度到統(tǒng)計(jì)串列BListY[2K]存點(diǎn)的緯度到統(tǒng)計(jì)串列Norm1:"NUMBER=":K 顯示點(diǎn)號(hào)ListX[2K]JProg"DM-S":Fix4 調(diào)子程序分離經(jīng)度的度分與秒值"L(DM)=":Z[13] 顯示點(diǎn)的經(jīng)度到度分"+L(S)=":Z[14] 顯示點(diǎn)的經(jīng)度到秒ListY[2K]JProg"DM-S"調(diào)子程序分離緯度的度分與秒值"B(DM)=":Z[13] 顯示點(diǎn)的緯度到度分"+B(S)=":Z[14] 顯示點(diǎn)的緯度到秒Prog"ZS"調(diào)正算子程序檢核RListFreq[2K]存子午線收斂角到統(tǒng)計(jì)串列RJProg"DM-S"調(diào)子程序分離經(jīng)度的度分與秒值"r(DM)=":Z[13] 顯示點(diǎn)的子午線收斂角到度分"+r(S)=":Z[14] 顯示點(diǎn)的子午線收斂角到秒"GAUSSPROJECTZSCHECK."提示高斯投影正算檢核結(jié)果Norm1 設(shè)置顯示格式X(m)=":ListX[2K-1]-Z[6]-X 顯示反算檢核的X坐標(biāo)Y(m)=":List顯示反算檢核的X坐標(biāo)"COMPNextPOINTYES(1)OrNO(0)="?Z[8]IfZ[8]=1:ThenGoto1:IfEnd"PR11END"分離角度的度分與秒值子程序：DM-S入口參數(shù)為角度變量J，輸出結(jié)果的度分存儲(chǔ)在Z[13]中，秒值存儲(chǔ)在Z[14]中。Int(J)Z[11]取出經(jīng)度值的度數(shù)Int(60(J-Z[11]))Z[12]取出經(jīng)度值的分?jǐn)?shù)Z[11]+Z[12]÷60Z[13]取出經(jīng)度值的度分值3600(J-Z[13])Z[14]取出經(jīng)度值的秒數(shù)PR10使用的同名子程序完全相同，其中正算子程序的作用為檢核反算計(jì)算結(jié)果的正確性。0.00015480IUGG1975X，Y0響應(yīng)。LBrL，BX，YmXY坐標(biāo)差值。fx-5800P的數(shù)值計(jì)算經(jīng)度非常高，一般情況下，正算X坐標(biāo)差<2×105Y坐標(biāo)差<×10511每完成一個(gè)點(diǎn)的換帶計(jì)算后屏幕“PtPOINT)r按鍵為繼續(xù)下一個(gè)點(diǎn)的反算，按 鍵為結(jié)束程序，屏幕顯示“。進(jìn)行高斯投影反算成果案例1954北京坐標(biāo)系，中央子午線經(jīng)度L0=111°,x坐標(biāo)加常數(shù)=0m,y坐標(biāo)加常數(shù)=500000m序 X(m) Y(m)

X檢核差

Y檢核差L(L(°′″)B(°′″)r(°′″)X(m)Y(m)12507005.772714567.79791130513.685223849.52604814.28-1.709×10-58.2528×10-522494239.989714212.26111130455.002223154.96704753.17-1.697×10-58.2052×10-532506658.355711330.65631130320.216223839.70404730.22-1.615×10-57.6898×10-542505530.213710280.05141130242.902223803.52404714.64-1.583×10-57.5144×10-552503396.508709934.88691130229.798223654.35904707.32-1.572×10-57.4594×10-562503369.195712593.95111130402.832223652.28004743.06-1.649×10-57.9105×10-572505724.695712022.37251130343.973223809.06104738.36-1.632×10-57.8088×10-582505514.694715397.88511130542.003223800.70704823.55-1.734×10-58.4038×10-5完成計(jì)算后用戶輸入的點(diǎn)的坐標(biāo)分別存儲(chǔ)在ListX與ListY的奇數(shù)行如1行3行、5行等反算結(jié)果大地經(jīng)度緯度B與子午線收斂角r分別存儲(chǔ)在ListFreq的偶數(shù)行如2行3行5行等……用戶可以按 (REG)鍵進(jìn)入雙變量統(tǒng)計(jì)與回歸模式查看。12、程序PR12，54北京系與80西安系高斯投影換帶計(jì)算主程序：PR12"GAUSSPROJECTHDPR12"顯示程序標(biāo)題基本設(shè)置Lbl0"COORDINATESYSTEM(54)Or(80)="?W坐標(biāo)系數(shù)字，只能輸入數(shù)值54或80IfW54:ThenIfW80:Then"NOTBE54Or80,REINPUT.":Goto0:IfEnd:IfEnd10DimZ 定義額外數(shù)組變量Prog"COEFFICIENT""CENTERL01(Deg)="?I輸入源坐標(biāo)中央子午線經(jīng)度"CENTERL02(Deg)="?J輸入目標(biāo)坐標(biāo)中央子午線經(jīng)度"X1+CONSTANT(m)="?Z[6]輸入源坐標(biāo)的X坐標(biāo)加常數(shù)"Y1+CONSTANT(m)="?Z[7]輸入源坐標(biāo)的Y坐標(biāo)加常數(shù)"X2+CONSTANT(m)="?Z[8]輸入目標(biāo)坐標(biāo)的X坐標(biāo)加常數(shù)"Y2+CONSTANT(m)="?Z[9]輸入目標(biāo)坐標(biāo)的Y坐標(biāo)加常數(shù)0K 點(diǎn)計(jì)數(shù)變量清零Lbl1K+1K 點(diǎn)計(jì)數(shù)變量計(jì)數(shù)"POINTX(m)="?X輸入點(diǎn)的X坐標(biāo)"POINTY(m)="?Y輸入點(diǎn)的Y坐標(biāo)ListX[2K-1]存點(diǎn)的X坐標(biāo)到統(tǒng)計(jì)串列ListY[2K-1]存點(diǎn)的Y坐標(biāo)到統(tǒng)計(jì)串列X-Z[6]X去掉X坐標(biāo)加常數(shù)后的X坐標(biāo)Y去掉Y坐標(biāo)加常數(shù)后的Y坐標(biāo)Prog"FS"調(diào)反算子程序L+I-JL 計(jì)算點(diǎn)在目標(biāo)坐標(biāo)系的經(jīng)度差Prog"ZS"調(diào)正算子程序X+Z[8]ListX[2K]存加常數(shù)后的高斯X坐標(biāo)到統(tǒng)計(jì)串列Y+Z[9]ListY[2K]存加常數(shù)后的高斯Y坐標(biāo)到統(tǒng)計(jì)串列RListFreq[2K]存子午線收斂角到統(tǒng)計(jì)串列Norm1:"NUMBER=":K 顯示點(diǎn)號(hào)"X(m)=":ListX[2K] X坐標(biāo)"Y(m)=":ListY[2K] Y坐標(biāo)"r(DMS)=":R 顯示子午線收斂角"GAUSSPROJECTHDCHECK.提示高斯投影正算檢核結(jié)果Prog"FS"調(diào)反算子程序L+J-IL 計(jì)算點(diǎn)在源坐標(biāo)系的經(jīng)度差Prog"ZS"調(diào)正算子程序RListFreq[2K-1]存源坐標(biāo)系子午線收斂角到統(tǒng)計(jì)串列X(m)=":X+Z[6]-ListX[2K-1] 顯示源坐標(biāo)檢核計(jì)算結(jié)Y(m)=":Y+Z[7]-ListY[2K-1] 顯示源坐標(biāo)檢核計(jì)算結(jié)"COMPNextPOINTYES(1)OrNO(0)="?Z[10]IfZ[10]=1:ThenGoto1:IfEnd"PR12END"PR10使用的同X1，Y1→L，BB→X2，Y2，r2X2，Y2→L，BL，B→X1，Y1，r1。5480IUGG1975X，Y0響應(yīng)；目標(biāo)坐標(biāo)系中央子午線經(jīng)度，X，Y0響應(yīng)。由計(jì)算出源坐標(biāo)系的X1，Y1mXY12-1所示。每完成一個(gè)點(diǎn)的換帶計(jì)算后屏幕“PtPOINT)r按鍵為繼續(xù)下一個(gè)點(diǎn)的換帶計(jì)算，按 鍵為結(jié)束程序，屏幕顯示“。表12-1使用程序PR12進(jìn)行高斯投影換帶成果案例源坐標(biāo)：1980西安坐標(biāo)系，3.38帶，中央子午線經(jīng)度L0=114°,x坐標(biāo)加常數(shù)=0m,y坐標(biāo)加常數(shù)=500000mr(°′″)X(m)Y(m)12505729.753406105.9442506945.934714504.163404813.451.937×10-5-8..0677r(°′″)X(m)Y(m)12505729.753406105.9442506945.934714504.163404813.451.937×10-5-8..0677×10-522492979.498405493.9182494180.107714148.597104752.341.924×10-5-8.0074×10-532491915.738390106.5012492807.288698779.514804424.731.512×10-5-5.4332×10-542493333.329422112.6712494867.775730765.527005135.952.51×10-5-1.1459×10-452480319.093404845.4842481503.502713753.578004730.901.911×10-5-7.9421×10-562475928.8893907