PC機測量平差程序

TrueBASIC1.0語言是BASIC語言的創(chuàng)始人JohnG.Kemeny和ThomasE.Kurtz等人于1985年根據(jù)美國國家BASIC標(biāo)準(zhǔn)設(shè)計的一個典型的結(jié)構(gòu)化程序設(shè)計語言中國科學(xué)院軟件研究所引進后，對其進行了漢化漢化后的TrueBASIC可以處理漢字字符1988年1月推出了2.03版，這也是TrueBASIC的最終版本TrueBASIC保留了BASIC語言清晰、易懂、易學(xué)的特點屏棄了非結(jié)構(gòu)化BASIC版本中陳腐有害的內(nèi)容是一個適合于大型數(shù)值計算具有較高和完整性能的結(jié)構(gòu)化程序設(shè)計語言尤其是它的矩陣運算語句MAT非常適合于編制測量平差計算程序TrueBASIC測量平差程序程序段展示了TrueBASIC的矩陣運算語句MAT的魅力平差計算程序變得簡單易讀為使程序易讀，在語句行后輸入中文注釋在注釋文字前加“！”字符實現(xiàn)。TrueBASIC語言編寫的測量計算程序，編譯并匯集成成擴展名為.exe的可執(zhí)行程序?？稍谌我釽indows下運行已知數(shù)據(jù)輸入——從用戶指定的文本格式文件讀入。計算結(jié)果輸出——存入SU文件,CS文件,SV文件。SU文件——存檔，CS文件——CASS展點或上傳到全站儀，SV文件——SV300展點。4.2水準(zhǔn)網(wǎng)條件平差程序(PG4-2)水準(zhǔn)網(wǎng)特點——構(gòu)成閉、附合環(huán)不多，每條水準(zhǔn)路線上的未知點較多，間接平差法——法方程階數(shù)較大，占用PC機內(nèi)存較多。水準(zhǔn)網(wǎng)較適合采用條件平差。PG4-2.exe程序信息設(shè)計1.點編號從1開始編起，先編已知點，后編未知點，未知點的編號順序沒有要求。某四等水準(zhǔn)路線湘1-03-49(1)——國家二等水準(zhǔn)點，網(wǎng)內(nèi)唯一已知點，其余點名為I點的17個點為一級導(dǎo)線點，“龍王橋”和“麻山塘”為新埋設(shè)的兩個四等水準(zhǔn)點，共有19個未知點，點編號如圖所示。

水準(zhǔn)網(wǎng)條件平差點編號案例1水準(zhǔn)網(wǎng)條件平差點編號案例21.點編號從1開始編起，先編已知點，后編未知點，未知點的編號順序沒有要求。2.水準(zhǔn)測段觀測數(shù)據(jù)的輸入每個水準(zhǔn)測段觀測數(shù)據(jù)的輸入格式為：起訖點號,測段高差(m),路線長(km)每個水準(zhǔn)測段觀測數(shù)據(jù)占一行，6點→5點水準(zhǔn)測段輸入數(shù)據(jù)為：6.05,-8.615,0.95876.05——起訖點號，整數(shù)位為起點號，小數(shù)位為訖點號，根據(jù)總點數(shù)的不同分兩種情形輸入：(1)水準(zhǔn)網(wǎng)總點數(shù)<100，小數(shù)位訖點號為兩位數(shù)，不足兩位數(shù)前面添0湊足兩位數(shù)。如6.05的訖點號為5，而6.5的訖點號為50。(2)水準(zhǔn)網(wǎng)總點數(shù)>100,<1000時，小數(shù)位的訖點號保證三位數(shù)，不足三位數(shù)前面添0湊足。如6.005的訖點號為5，6.05的訖點號為50，而6.5的訖點號為500。水準(zhǔn)測段的高差觀測方向按實際觀測方向定，不要求統(tǒng)一到一個方向，不需要為水準(zhǔn)測段編號，程序根據(jù)測段的起訖點號自動搜索平差計算所需水準(zhǔn)測段數(shù)據(jù)。3.閉合或附合水準(zhǔn)路線條件方程點號列立條件方程時推算高差閉合差的路徑，閉合水準(zhǔn)路線點號的起訖點號應(yīng)相同，附合水準(zhǔn)路線點號的起訖點號不同，應(yīng)為已知點。每條水準(zhǔn)路線點號數(shù)據(jù)占一行，水準(zhǔn)路線點號數(shù)=水準(zhǔn)網(wǎng)多余觀測數(shù)。作用——列條件方程和計算平差后未知點高程，要求第一條水準(zhǔn)路線點號的起點點號應(yīng)為已知點，其余水準(zhǔn)路線點號的起點點號不要求。下圖水準(zhǔn)測段高差觀測數(shù)=22，未知點數(shù)=19，多余觀測數(shù)22-19=3，3個條件的水準(zhǔn)路線點號可以是：1,2,3,4,5,6,16,17,18,11,19,20,15,111,19,20,15,14,13,12,116,7,8,9,10,11,18,17,16,6下圖水準(zhǔn)測段高差觀測數(shù)=12，未知點數(shù)=3，多余觀測數(shù)=12-3=9。9個條件可以是下列4個閉合水準(zhǔn)路線和5個附合水準(zhǔn)路線：1,9,7,17,8,9,72,7,8,25,8,9,51,7,22,8,33,8,44,8,55,9,64.平差值函數(shù)水準(zhǔn)路線點號平差值函數(shù)只能為未知點高程。為計算15,6,17,11四個未知點高程平差值的誤差，應(yīng)給出這些點從網(wǎng)中唯一已知點——1號點推算高程的點號路線：1,151,2,3,4,5,61,2,3,4,5,6,16,171,15,20,19,11為計算7,8,9三個點高程平差值的誤差，給出這些點從網(wǎng)中任一已知點推算高程的點號路線：1,72,86,9水準(zhǔn)網(wǎng)條件平差已知數(shù)據(jù)文件案例標(biāo)題每km全中誤差已知點數(shù),未知點數(shù),測段高差數(shù)點名已知點高程22個測段高差與路線長條件方程路線點號計算高程平差值中誤差推算點號每km高差中數(shù)全中誤差限差二等——2mm三等——6mm四等——10mm圖根——22mm程序PG4-2.exe運行界面水準(zhǔn)網(wǎng)條件平差已知數(shù)據(jù)文件案例4.3三角高程網(wǎng)間接平差程序(PG4-3)1.點編號從1開始編起，先編已知高程點，后編未知點，未知點的編號順序沒有要求。2.測邊數(shù)據(jù)的輸入格式一條三角高程測邊觀測數(shù)據(jù)的輸入格式為：起訖點號,高斯或參考橢球面的邊長(m),往測豎角(dms),儀高(m),覘標(biāo)高(m),返測豎角(dms),儀高(m),覘標(biāo)高(m)若為單向高差觀測，則要求是往測邊方向，返測邊的豎角,儀高,覘標(biāo)高均應(yīng)輸0。測邊的輸入順序可以是任意順序?！袄舷x坡→小院子”測邊的輸入數(shù)據(jù)為：1.05,3353.4068,-2.39463,1.52,3.12,2.43515,1.67,5.726.近似高程推算路線與三角形高差閉合差計算已知數(shù)據(jù)文件尾輸入近似高程推算路線點號。目的——計算未知點高程近似值與三角形高差閉合差。為避免三角形高差閉合差重復(fù)，可編寫多條近似高程推算路線，路線數(shù)由程序自動計算，不需要輸入。下面為兩條近似高程推算路線，要求在第一條路線的前3個點中至少有一個是已知點，且在后面任意一條路線的前3個點中，至少有一個點的點號出現(xiàn)在此前的路線中。1,11,7,-10,-9,-8,4,-6,-2,-3,-5,-11,12,-13,-11,12三角高程網(wǎng)間接平差已知數(shù)據(jù)文件案例數(shù)據(jù)項1——標(biāo)題字符。數(shù)據(jù)項2——測區(qū)平均曲率半徑Rm(m)，x坐標(biāo)加常數(shù)數(shù)據(jù)項3——邊長投影高程面高程，邊長投影面字符為高斯面，Y，否則為參考橢球面。數(shù)據(jù)項4——已知點數(shù)，未知點數(shù)，高差觀測測邊數(shù)。數(shù)據(jù)項5——點名字符。數(shù)據(jù)項6——全部點的概略橫坐標(biāo)值y，邊長為高斯平面長度時才需輸入。數(shù)據(jù)項7——已知點高程。數(shù)據(jù)項8——測邊高差觀測數(shù)據(jù)，每條測邊的往返對向高差觀測數(shù)據(jù)應(yīng)在一行內(nèi)輸完。起訖點號,邊長值(m),往測豎角(ddd.mmss),測站高(m),覘標(biāo)高(m),返測豎角(ddd.mmss),測站高(m),覘標(biāo)高(m)單向觀測測邊，將其編為往測數(shù)據(jù)，返測豎角、測站高、覘標(biāo)高均應(yīng)輸0。數(shù)據(jù)項9——未知點高程推算與三角形高差閉合差計算路線點號。程序PG4-3.exe運行界面4.4光電測距三角高程導(dǎo)線網(wǎng)概算與條件平差程序(PG4-4)1.點編號從1開始編起，先編已知高程點，后編未知點，對未知點編號順序無要求。2.三角高程導(dǎo)線邊的輸入順序與格式

一條三角高程導(dǎo)線邊觀測數(shù)據(jù)的輸入格式為：起訖點號,往測斜距(m),往測豎角(dms),儀高(m),覘標(biāo)高(m),返測斜距(m),返測豎角(dms),儀高(m),覘標(biāo)高(m)“I45→高頂山”三角高程導(dǎo)線邊的輸入數(shù)據(jù)為：8.01,1316.6287,3.1318,1.605,1.535,1316.6374,-3.1434,1.55,1.365單向觀測邊，應(yīng)是往測邊，返測邊的“斜距,豎角,儀高,覘標(biāo)高”均應(yīng)輸0。

三角高程觀測邊數(shù)=57，未知點數(shù)也即必要觀測數(shù)=46，多余觀測數(shù)也即條件方程數(shù)=11?？闪谐?個閉合三角高程路線,6個附合三角高程路線。點號為：3,13,14,15,16,17,4,35,34,33,32,31,30,33,30,31,32,33,34,35,4,18,19,20,21,5,39,38,37,36,2,11,12,32,36,37,38,39,5,22,23,24,43,42,41,40,22,40,41,42,43,24,50,49,48,47,46,45,44,246,47,48,49,50,24,6,25,26,53,52,51,461,8,9,10,22,11,12,33,13,14,15,16,17,44,18,19,20,21,55,22,23,24,66,25,26,27,28,29,73.平差值函數(shù)水準(zhǔn)路線點號平差值函數(shù)只能為未知點高程。為計算24,26,46三個未知點高程平差值的誤差，給出從網(wǎng)中任意已知點推算這些點高程的點號路線：6,246,25,262,44,45,46光電測距三角高程導(dǎo)線網(wǎng)平差已知數(shù)據(jù)文件數(shù)據(jù)項1——標(biāo)題字符。數(shù)據(jù)項2——測區(qū)平均曲率半徑Rm,如<6000000,取6371000。數(shù)據(jù)項3——三角高程測量等級數(shù)字,4——四等,5——五等。數(shù)據(jù)項4——已知點數(shù),未知點數(shù),三角高程導(dǎo)線測邊數(shù)n。數(shù)據(jù)項5——點名字符。數(shù)據(jù)項6——已知點高程。數(shù)據(jù)項7——三角高程導(dǎo)線邊觀測數(shù)據(jù),格式為：起訖點號,往測斜距(m),往測豎角(dms),儀高(m),覘標(biāo)高(m),返測斜距(m),返測豎角(dms),儀高(m),覘標(biāo)高(m)

光電測距三角高程導(dǎo)線網(wǎng)平差已知數(shù)據(jù)文件數(shù)據(jù)項8——閉附合路線點號，每環(huán)線或附合路線點號占一行，行數(shù)應(yīng)等于網(wǎng)的多余觀測數(shù)。數(shù)據(jù)項9——需評定未知點高程中誤差的平差值函數(shù)路線點號，從已知點號開始，沿網(wǎng)中測邊經(jīng)歷的點號至需要評定高程精度的未知點止，圖中最后三行為評定24，26，46點高程精度。光電測距三角高程導(dǎo)線網(wǎng)平差已知數(shù)據(jù)文件程序PG4-4.exe運行界面4.5三角網(wǎng)、邊角網(wǎng)概算與按方向附條件間接平差程序(PG4-5)1.點編號從1開始編起，先編已知高程點，后編未知點，對未知點編號順序無要求。2.近似坐標(biāo)計算與坐標(biāo)推算路線點號保證A,B,P三點構(gòu)成的旋轉(zhuǎn)方向為逆時針方向,實際情況相符，保證三角型閉合差計算。三種情形討論。(1)有兩個以上已知點且已知點間的連線為三角網(wǎng)的觀測邊“橋東飯店”,“南山”為兩個已知點，已知點連線為三角網(wǎng)的一條方向觀測邊，坐標(biāo)推算路線的頭兩個點為這兩個已知點，沿坐標(biāo)推算路線計算一次即可得到全網(wǎng)未知點的真坐標(biāo)。1,2,-9,-10,-3,-4,2,-5,-6,-7,-8,-9(2)線形鎖——坐標(biāo)推算路線的頭一個點為已知點，箭頭所示方向的坐標(biāo)推算路線點號為：1,6,7,5,3,-4,-2(3)一個已知點，一個或多個已知邊及已知方位角已知點“觀音山”，有兩條已知基線邊和兩個已知方位角。圖中箭頭所示方向的兩條坐標(biāo)推算路線點號分別為：1,2,4,-3,-5,-6,7,11,12,131,4,8,-7,9,-10,-11數(shù)據(jù)項1——標(biāo)題字符。數(shù)據(jù)項2——等級三角網(wǎng)與邊角網(wǎng)測角中誤差允許值1→二等,1.8→三等,2.5→四等,5→一級,10→二級。數(shù)據(jù)項3——是否計算邊長的橢球與高斯投影改正?偏心改正?Y或y，計算邊長的橢球與高斯投影改正，Y或y，計算偏心改正。數(shù)據(jù)項4——測區(qū)平均曲率半徑，橫坐標(biāo)加常數(shù)。數(shù)據(jù)項5——高程異常，邊長投影高程面的高程H0。只有當(dāng)Y時才需要。數(shù)據(jù)項6——已知點數(shù)，未知點數(shù)，方向觀測數(shù)，邊長數(shù)，已知方位角數(shù)。沒有測距邊或沒有已知方位角時，應(yīng)輸0。

數(shù)據(jù)項7——點名字符，可在一行內(nèi)輸完，也可分行輸入。分行輸入時，中間行應(yīng)以“,”結(jié)束。數(shù)據(jù)項8——點的概略高程值，應(yīng)在一行內(nèi)全部輸完。只有當(dāng)Y時，才需要輸入。數(shù)據(jù)項9——已知點平面坐標(biāo)。數(shù)據(jù)項10——已知方位角數(shù)據(jù)，格式為：起訖點號，方位角。只有當(dāng)方位角數(shù)≠0時才需要輸入。

數(shù)據(jù)項11——按點號順序的方向觀測值與覘點號每行一個點的全部方向觀測值和覘點號，共13個點。格式：方向觀測值，覘點號。數(shù)據(jù)項12——全網(wǎng)測站歸心元素和照準(zhǔn)點歸心元素，13個點。只有當(dāng)eYT$=“Y”或“y”時才需要輸入。數(shù)據(jù)項13——已知基線邊或測距邊數(shù)據(jù)，每條邊長的數(shù)據(jù)占一行，格式為:起訖點號，邊長值,權(quán)。當(dāng)邊長為已知基線邊時，其權(quán)可輸2000。數(shù)據(jù)項14——近似坐標(biāo)推算路線點號每條坐標(biāo)推算路線點號數(shù)據(jù)占用一行。程序PG4-5.exe運行界面4.6導(dǎo)線網(wǎng)概算與按方向與邊長間接平差程序(PG4-6)1.點編號從1開始編起，先編已知坐標(biāo)點，后編未知坐標(biāo)點，未知坐標(biāo)點的編號順序沒有要求。

2.近似坐標(biāo)計算與坐標(biāo)推算路線四種類型：(1)附合導(dǎo)線——坐標(biāo)推算路線最前面兩個點和最尾部兩個點為已知坐標(biāo)點或其坐標(biāo)已經(jīng)在該坐標(biāo)推算路線前的坐標(biāo)推算路線中計算出。1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,1,25,6,30,31,32,33,34,35,26,27(3)單邊定向?qū)Ь€只能計算出坐標(biāo)增量閉合差，無方位角閉合差，1,2,4,5,6,3(4)無定向?qū)Ь€這類導(dǎo)線兩端的兩個已知點1,2號點不通視，先在假定坐標(biāo)系中推算未知點的假定坐標(biāo)，再將其變換到真坐標(biāo)系。1,3,4,5,2數(shù)據(jù)項1——標(biāo)題字符。數(shù)據(jù)項2——等級導(dǎo)線網(wǎng)測角中誤差允許值1.5→三等,2.5→四等,5→一級,8——二級,12→三級。測距儀標(biāo)稱測距精度，5mm+5ppm→5,5數(shù)據(jù)項3——是(y)否計算邊長的橢球與高斯投影改正，導(dǎo)線邊長是(y)否已為水平距離。數(shù)據(jù)項4——測區(qū)平均曲率半徑，橫坐標(biāo)的加常數(shù)。數(shù)據(jù)項5——高程異常，邊長投影高程面的高程。只有當(dāng)project$=“Y”時才需要。數(shù)據(jù)項6——已知點數(shù)，未知點數(shù)，方向觀測數(shù)，邊長數(shù)。湖南望城火電廠大眾坡方案一級導(dǎo)線平面網(wǎng)嚴(yán)密平差計算資料5,5,5Y,N6366335,385000000,02,51,116,57數(shù)據(jù)項7——點名字符。數(shù)據(jù)項8——點概略高程值，當(dāng)project$=“Y”

時，才輸入。數(shù)據(jù)項9——已知點的平面坐標(biāo)。高頂山,高沖,I49,I48,I47,I46,I17,I18,I19,I20,I21,I5,I4,I3,I2,I1,I30,I42,I41,I40,I39,I38,I37,I29,I28,I10,I43,I44,I45,I22,I23,I24,I25,I26,I27,I52,I51,I50,I6,I7,I8,I9,I31,I32,I33,I34,I35,I36,I53,I16,I15,I14,I13156.96,96.22,69.86,58.26,44.37,32.10,35.91,42.796,56.91,62.95,69.31,39.608,30.995,31.442,36.207,35.982,36.201,36.24,35.54,35.363,50.25,69.13,71.559,73.06,81.91,80.465,50.93,55.33,82.78,33.34,36.05,36.98,72.73,66.38,90.11,51.98,56.21,77.19,67.799,60.50,68.03,74.65,35.721,36.05,37.588,30.45,58.10,71.67,46.72,61.78,52.037,54.706,74.313130334.4515,38387901.6828,3134621.62,38384764.47數(shù)據(jù)項10——導(dǎo)線點方向觀測值與覘點號每行輸一個點的全部方向觀測值和覘點號，每個方向觀測的輸入格式為：方向觀測值，覘點號。0,2,104.4701,290,3,40.0841,10,4,254.0213,20,5,84.5819,30,6,178.4755,40,7,29.1807,30,201.0704,50,8,198.5540,60,9,159.3435,70,10,31.2620,49,155.0625,50,185.121,80,11,194.5230,90,12,185.181,100,13,148.0445,39,187.5427,110,14,227.0234,120,15,167.2806,130,16,88.2648,140,15,189.2322,170,16,181.5724,18,271.0652,430,17,276.2339,190,18,158.1622,200,19,181.3314,210,20,179.5442,220,21,175.5938,230,24,131.0145,48,173.0024,220,23,207.4050,250,24,249.2348,260,25,45.0006,27,129.3753,35,187.1409,38,252.0405,420,28,142.5029,260,29,185.4840,270,1,145.1711,280,31,144.4501,60,30,189.1903,320,31,161.4545,330,32,77.3405,53,95.2619,340,33,264.3524,350,26,253.1406,340,37,143.2806,490,38,199.4501,360,26,201.0819,370,12,144.2736,400,39,220.2800,410,40,169.1126,420,41,174.2731,260,17,153.1808,440,43,214.4150,450,44,221.3906,460,45,82.4244,470,46,153.3130,480,47,195.0608,230,36,163.2038,90,51,284.3656,90,52,175.4900,500,51,270.3850,530,52,171.3514,33數(shù)據(jù)項11——測距邊數(shù)據(jù)，每條邊長數(shù)據(jù)占一行，格式分為兩種情況，hor$=“Y”時為：起訖點號,水平邊長值；hor$≠“Y”時為：起訖點號,往測斜距,往測豎角,返測斜距,返測豎角。2.03,1117.3677,-1.2048,1117.3626,1.20143.04,1049.2690,-0.3850,1049.2608,0.37404.05,692.6772,-1.0902,692.6733,1.08325.06,589.5769,-1.1136,589.5736,1.11106.07,684.4627,0.1830,684.4597,-0.18597.08,361.3508,1.0506,361.3500,-1.05388.09,703.9727,1.0834,703.9697,-1.09139.10,617.9871,0.3250,617.9890,-0.331910.11,647.5212,0.3424,647.5207,-0.345511.12,728.5811,-2.2050,728.5770,2.202712.13,721.5765,-0.4122,721.5711,0.404413.14,618.3475,0.0234,618.3465,-0.030614.15,565.5008,0.2833,565.4979,-0.290215.16,578.9502,-0.0220,578.9322,0.021816.17,303.9428,0.0240,303.9429,-0.031517.18,550.3615,0.0026,550.3548,-0.002618.19,387.4466,-0.1648,387.4434,0.161819.20,695.8640,0.0402,695.8654,-0.043420.21,577.9320,1.2844,577.9376,-1.291621.22,405.6676,2.3930,405.6734,-2.395922.23,296.5693,0.2800,296.5823,-0.290023.24,599.3939,0.0826,599.4045,-0.090124.25,1051.0746,0.2804,1051.0805,-0.284025.26,894.8822,-0.0516,894.8884,0.044226.27,538.9529,-3.0832,538.9199,3.070627.28,395.6487,0.3650,395.6519,-0.400028.29,921.3101,1.4229,921.3148,-1.440229.01,1316.6623,3.1318,1316.6675,-3.14346.30,775.1006,0.0552,775.1015,-0.060130.31,444.2571,0.2034,444.2594,-0.211031.32,337.7719,0.0904,337.7780,-0.094432.33,663.8379,3.0512,663.8437,-3.054033.34,777.6143,-0.2826,777.6172,0.280034.35,424.4006,3.1248,424.3883,-3.115635.26,733.0975,-0.4527,733.0913,0.445426.38,474.9799,-0.2450,474.9835,0.241838.37,425.9490,-2.4940,425.9425,2.482437.36,414.8751,-0.3658,414.8828,0.335036.49,327.0196,-0.5740,327.0175,0.540449.09,401.5972,1.2548,401.5951,-1.270012.39,468.8005,3.2543,468.7998,-3.261639.40,337.4099,-1.1346,337.4098,1.131340.41,492.9877,0.5146,492.9919,-0.521841.42,332.7452,1.0757,332.7492,-1.083042.26,577.8875,0.3434,577.8828,-0.345617.43,342.8990,-0.0430,342.9008,0.040043.44,537.7473,0.0150,537.7537,-0.022244.45,337.9014,0.1500,337.9051,-0.153445.46,682.3638,-0.3634,682.3653,0.361246.47,533.5230,2.5812,533.5229,-2.584247.48,660.3458,1.1038,660.3516,-1.110348.23,326.5478,-0.0156,326.5476,0.01269.50,609.0899,0.2707,609.0946,-0.275250.51,619.8298,-0.5430,619.8299,0.535751.52,541.2319,0.1625,541.2295,-0.170452.53,368.8258,3.0112,368.8379,-3.021053.33,491.9328,-0.1128,491.9330,0.1104數(shù)據(jù)項12——近似坐標(biāo)推算路線點