【正文】 esultfp,\n 單位權(quán)中誤差: m=%lf,m_mu)。 PrintResult()。 PrintM2(resultfp,ATPA,m_Pnumber,4,%15e ,Qx)。 }主函數(shù)的功能就是調(diào)用所有其它非主函數(shù)，并最終在主函數(shù)的調(diào)用下，執(zhí)行這個(gè)程序。函數(shù)的代碼如下：void main(){ char *datafile =。 char *resultfile =。 char *QuasiStablefile =。 printf(\n\n擬穩(wěn)平差計(jì)算\n)。 printf(數(shù)據(jù)文件： %s\n,datafile)。 printf(結(jié)果文件： %s\n,resultfile)。 FILE *fp=fopen(resultfile,w)。 if (fp==NULL) { MyBreak(創(chuàng)建結(jié)果文件失敗！)。 return。 } CLevelingAdjust level。 =fp。 (datafile)。 //()。 (QuasiStablefile)。 fclose(fp)。}水準(zhǔn)網(wǎng)平差所需的已知數(shù)據(jù)和觀測(cè)數(shù)據(jù)全部從數(shù)據(jù)文件中讀取，所以將水準(zhǔn)網(wǎng)原始數(shù)據(jù)文件的內(nèi)容分為三類：網(wǎng)的概況、已知數(shù)據(jù)、觀測(cè)數(shù)據(jù)。網(wǎng)的概況包括：總點(diǎn)數(shù)、已知點(diǎn)數(shù)、觀測(cè)值總數(shù)、驗(yàn)前單位權(quán)中誤差。已知數(shù)據(jù)包括：已知點(diǎn)名、已知高程值。觀測(cè)數(shù)據(jù)包括：高差起點(diǎn)點(diǎn)名、高差終點(diǎn)點(diǎn)名、觀測(cè)高差、路線長(zhǎng)度。下面的數(shù)據(jù)文件內(nèi)容就是本程序的原始數(shù)據(jù)文件部分內(nèi)容的格式。97 98 2 CL06 CL02 CL01 CL02 CL02 CL03 CL03 ML001 ML001 ML002 ML002 ML003 ML003 ML004 ML004 ML005 格式說明：（1）第一行稱為網(wǎng)的概況信息，分別是觀測(cè)值總數(shù)、總點(diǎn)數(shù)、已知點(diǎn)總數(shù)、驗(yàn)前單位權(quán)中誤差。驗(yàn)前單位權(quán)中誤差即觀測(cè)高差的每千米高差中誤差，以米為單位。（2）第第3行是已知點(diǎn)名及點(diǎn)名對(duì)應(yīng)的已知高程值，高程值以米為單位。（3）第4行至第10行為觀測(cè)高差，每行的內(nèi)容為高差起始點(diǎn)點(diǎn)名、高差終點(diǎn)點(diǎn)名、高差觀測(cè)值和高差觀測(cè)路線長(zhǎng)度，高差以米為單位，路線長(zhǎng)度以千米為單位。（4）網(wǎng)中的點(diǎn)名在程序中作為字符串處理，點(diǎn)名可以包含漢字、字母和數(shù)字等，字母區(qū)分大小寫，每一個(gè)點(diǎn)只能有唯一的點(diǎn)名，點(diǎn)名中間不能有空格和控制字符。為了進(jìn)行擬穩(wěn)平差，首先要知道哪些點(diǎn)是擬穩(wěn)點(diǎn)，哪些點(diǎn)不是擬穩(wěn)點(diǎn)。，文件中包含擬穩(wěn)點(diǎn)數(shù)和全部擬穩(wěn)點(diǎn)的點(diǎn)名，平差程序可以從這個(gè)文件中讀取到平差所需的擬穩(wěn)點(diǎn)名，file文件的格式如下：60CL01 CL02 CL03 ML001 ML002 ML003 ML004 ML005 ML006 ML007 ML008 ML009 ML010 ML011 ML012 ML013 ML016 ML017 ML018 ML021 ML022格式說明：第1行為擬穩(wěn)點(diǎn)的總個(gè)數(shù)，第2行以下是擬穩(wěn)點(diǎn)的點(diǎn)名，點(diǎn)名之間用空格隔開。[1]監(jiān)測(cè)網(wǎng)的目的是精確提供未知量的變化?，F(xiàn)有的平差方法或者固定若干未知量作強(qiáng)制符合，或者作自由網(wǎng)平差，全部未知量(包括穩(wěn)定的和不穩(wěn)定的)擬合于它們的起始值，改正數(shù)范數(shù)極小。前一種方法不免使觀測(cè)受到歪曲，損害精度；后一種方法即偽逆解法，其缺點(diǎn)是未知量沒有穩(wěn)定的基準(zhǔn)，因而不能提供所要求的變化…。“擬穩(wěn)平差”的基本思想是，考慮到監(jiān)測(cè)網(wǎng)中的點(diǎn)，處于不同的地質(zhì)構(gòu)造和地球物理環(huán)境，隨著時(shí)間的延伸，都可能發(fā)生變動(dòng)，但是總存在相對(duì)變化小的，即相對(duì)穩(wěn)定的點(diǎn)。我們稱之為“擬穩(wěn)點(diǎn)”。把它們作為基準(zhǔn)，附加合理的限制與擬穩(wěn)點(diǎn)有關(guān)的未知量范數(shù)極小的約束條件，原來(lái)的秩虧自由網(wǎng)可以得到準(zhǔn)確的確定解，這樣既有明確的測(cè)量根據(jù)，可以為監(jiān)測(cè)網(wǎng)提供比較前后變化的基準(zhǔn)面；又沒有強(qiáng)制附合，可保持網(wǎng)形不被歪曲，得到所需要的精確信息。擬穩(wěn)平差的關(guān)鍵是選取符合實(shí)際的“擬穩(wěn)點(diǎn)”，即確定陣。目前選擇擬穩(wěn)點(diǎn)的方法主要有假設(shè)檢驗(yàn)、糊聚類、名次法等技術(shù)途徑。假設(shè)檢驗(yàn)用于擬穩(wěn)點(diǎn)選擇的具體步驟如下：1)由兩期擬穩(wěn)平差結(jié)果對(duì)監(jiān)測(cè)網(wǎng)中的各點(diǎn)進(jìn)行單個(gè)檢驗(yàn)，確定出變形點(diǎn)和擬穩(wěn)點(diǎn)。2)對(duì)擬穩(wěn)點(diǎn)所在區(qū)域進(jìn)行整體檢驗(yàn)，然后由下面三條原則對(duì)所得結(jié)果進(jìn)行判別。三條判別原則：1)擬穩(wěn)平差中所選取的點(diǎn)必須是相對(duì)穩(wěn)定點(diǎn)：2)區(qū)分出的變形點(diǎn)應(yīng)處在不穩(wěn)定地帶，應(yīng)符合該區(qū)的形變特征；3)擬穩(wěn)點(diǎn)所構(gòu)成的區(qū)域必須是相對(duì)穩(wěn)定區(qū)域。因?yàn)閿M穩(wěn)點(diǎn)是選取相對(duì)穩(wěn)定的點(diǎn)，即擬穩(wěn)點(diǎn)的相對(duì)高差變化相對(duì)于其它點(diǎn)而言，高差的變化相對(duì)較小。所以我將第一觀測(cè)的數(shù)據(jù)與第二次觀測(cè)的數(shù)據(jù)、第三次觀測(cè)的數(shù)據(jù)做了相對(duì)比較，即用第一次的觀測(cè)數(shù)據(jù)減去第二次觀測(cè)的數(shù)據(jù)、第三次的觀測(cè)數(shù)據(jù)，從中挑選出高差變化較小的點(diǎn)。本次的平面坐標(biāo)數(shù)據(jù)展繪到AutoCAD中，呈十字形狀，如圖（41） 圖41 點(diǎn)位分布平面圖利用excel2003將分別將第一次觀測(cè)數(shù)據(jù)分別與第二次觀測(cè)數(shù)據(jù)和第三次觀測(cè)數(shù)據(jù)相減，運(yùn)行結(jié)果如下圖數(shù)據(jù)處理說明：1) 第一列的粉紅色表示：第一次觀測(cè)值與第二次觀測(cè)值具有相同的測(cè)站2) 第一列的黃色表示：第一次觀測(cè)值與第三觀測(cè)值具有相同的測(cè)站，而且在第三次觀測(cè)值中也用黃色標(biāo)在了對(duì)應(yīng)的測(cè)站位置3) 第七列的青色表示：這三次的觀測(cè)值具有相同的測(cè)站4) 第四列的紅色表示：第二次的觀測(cè)值與第三次的觀測(cè)值具有相同的測(cè)站5) 表格后面的1與1與2與3分別表示：第一次觀測(cè)數(shù)據(jù)與第二次觀測(cè)數(shù)據(jù)高差之差、第一次觀測(cè)數(shù)據(jù)與第三觀測(cè)數(shù)據(jù)高差之差、第二次觀測(cè)數(shù)據(jù)與第三次觀測(cè)數(shù)據(jù)高差之差6) 表中沒有用顏色標(biāo)記的行，表明該行的三次觀測(cè)數(shù)據(jù)的測(cè)站不一致，這點(diǎn)可以從表中直接觀察到。7) 通過表格最后三列的高差之差數(shù)據(jù)，可以較容易的判斷和選取高差之差變化相對(duì)較小的點(diǎn)，最后結(jié)合點(diǎn)位的平面分布圖，選擇擬穩(wěn)點(diǎn)8) 表格中沒有數(shù)據(jù)的單元格，表明這三次的測(cè)站不是同一測(cè)站，因?yàn)橹挥性谕粶y(cè)站的情況下，才能判斷點(diǎn)位的高差變化，通過點(diǎn)位的平面坐標(biāo)，使用Excel表格，計(jì)算出各個(gè)測(cè)站之間的距離，計(jì)算結(jié)果如下圖（51）： 圖51 點(diǎn)位坐標(biāo)與測(cè)站對(duì)應(yīng)圖表格中最后一列為測(cè)站間的平面距離，將表格中最后四列信息復(fù)制粘貼到“”文件中，放到第四行即可觀測(cè)數(shù)據(jù)的高差和距離的輸入工作?！啊敝?，其中擬穩(wěn)點(diǎn)選取按照（）中的選擇方法和判斷原則來(lái)選取，因?yàn)閿M穩(wěn)點(diǎn)的選取是不固定的，故本程序在此不易討論。 Visual C++，在“文件”下拉菜單中選擇“新建”，在“工程”下拉列表中選擇“Win32 console Application”，輸入工程名稱及保存的位置，將（）中的“”和“”兩個(gè)數(shù)據(jù)文件復(fù)制到該新建工程保存的目錄下?！靶陆ā?，在“文件”下拉列表中選擇“C++ Source File”，輸入源文件的名稱即可。，以“include文件名；”的形式添加C++庫(kù)函數(shù)頭文件。本程序包含的C++庫(kù)函數(shù)頭文件如下：include include include include include include include ，把（）中的函數(shù)源代碼按照順序全部復(fù)制粘貼到這個(gè)源代碼文件中即可。然后單機(jī)菜單欄下面的“組建”，從其下拉列表中單擊“編譯”。若編譯過程中沒有發(fā)現(xiàn)任何錯(cuò)誤，則表示編譯成功，此時(shí)在狀態(tài)輸出窗口顯示出“0 error(s),0 warning(s)”信息，若編譯過程中發(fā)現(xiàn)，則將在狀態(tài)輸出窗口顯示出發(fā)生錯(cuò)誤的文件、所在的行號(hào)和出錯(cuò)原因，將這些信息復(fù)制，粘貼到百度搜索網(wǎng)站中，即可找到錯(cuò)誤的原因和解決的方法?！敖M建”，從其下拉菜單中單擊“執(zhí)行”即可。程序運(yùn)行結(jié)束前，將在顯示器自動(dòng)顯示出“Press any key to continue”提示信息，此時(shí)按下鍵盤上任意鍵，即可關(guān)閉該窗口，重新回到VC++操作界面窗口。由mian主函數(shù)的建立的文件路徑可知，程序運(yùn)行的結(jié)果保存在新建工程的文件目錄下，保存在””文件中，””文件運(yùn)行的結(jié)果和顯示器上顯示的結(jié)果分別如圖（52）和圖（53） 圖52 ””文件運(yùn)行的結(jié)果 圖53 顯示器顯示結(jié)果結(jié)論自周江文教授將擬穩(wěn)平差的基本原理提出來(lái)以后，許多測(cè)量工作者將其運(yùn)用到實(shí)踐中，讓擬穩(wěn)平差在沉降監(jiān)測(cè)網(wǎng)的數(shù)據(jù)處理的應(yīng)用中，取得了重大的進(jìn)步，具有重要的意義。打破了國(guó)外原有經(jīng)典平差的局面，是平差的方法提升到了一個(gè)新的層次。本文的擬穩(wěn)平差程序設(shè)計(jì)借鑒了大量的圖書資料，秉承著“創(chuàng)新為主、借鑒為輔”的原則，成功的實(shí)現(xiàn)了從程序源代碼編寫到實(shí)現(xiàn)程序運(yùn)行。該套編程具有代碼簡(jiǎn)單、功能齊全、語(yǔ)言簡(jiǎn)練等特點(diǎn)參考文獻(xiàn)[1] 萬(wàn)斐, 陳艷艷. 擬穩(wěn)平差在測(cè)量數(shù)據(jù)處理中的應(yīng)用 [J]. 地理空間信息, 2008, 6(7982.[2]何國(guó)清,[M].徐州:中國(guó)礦業(yè)大學(xué)出版社,1991:2632. [2]陶本藻、邱衛(wèi)寧、[M].武漢:武漢大學(xué)出版社，2003:102142[3]胡振琪、王家貴、[M].北京:煤礦工業(yè)出版社,2008:305307.[1][M].北京:國(guó)防工業(yè)出版社,2008:154194.[1]++基礎(chǔ)教程[M].北京:清華大學(xué)出版社,2004:242243.謝辭感謝我的導(dǎo)師余學(xué)祥老師，他為人和藹，對(duì)學(xué)生認(rèn)真負(fù)責(zé)，對(duì)工作一絲不茍。在我進(jìn)行畢業(yè)論文的收集資料，初稿，最后定稿的各個(gè)階段都給予了我極大的幫助，尤其是當(dāng)我遇到一些困難請(qǐng)教老師時(shí)，他都會(huì)不予余力的給予我最大的幫助以及指引。這些都讓我受益匪淺。余學(xué)祥老師一直是我非常敬佩和尊重的教授，不僅是因?yàn)樗趯W(xué)術(shù)上的成就，更由于他對(duì)教學(xué)及科研的態(tài)度。余學(xué)祥老師能及時(shí)給予我建議，指出不足之處，幫助我更好地組織論文結(jié)構(gòu)，不斷充實(shí)論文的內(nèi)容。這些都是我能順利完成論文的必要條件最后，以我的誠(chéng)摯之心，致敬余老師。33