【正文】 分復(fù)雜，它是許多地質(zhì)采礦因素綜合影響的結(jié)果。沿走向主斷面布設(shè)的觀測線稱為走向觀測線，沿傾斜主斷面布設(shè)的觀測線稱為傾斜觀測線。 getchar()。 for(int i=0。i++) { fprintf(fp,\n%3d ,i+1)。 else a0[i]= ai0/a00。 for(int k=0。 } } return q。 //觀測值數(shù)組 double *Height。//擬穩(wěn)平差函數(shù) }。 im_Pnumber。 if((fp=fopen(datafile,r))==NULL) { MyBreak(\n 數(shù)據(jù)文件打不開!)。 fscanf(fp,%lf,amp。 //讀取高程起點名 StartP[i]=GetStationNumber(buffer)。}12. Printdata原始數(shù)據(jù)寫至結(jié)果文件函數(shù)此函數(shù)的功能就是將網(wǎng)的已知信息反饋到結(jié)果文件中。各點的高程值用Height數(shù)組保存，近似高程計算計算之前。 jj++。 for(int i=0。 ATPA[ij(i,j)]=Pk。以觀測值序號k為循環(huán)變量，按觀測值循環(huán)，按下式計算各觀測值的殘差，（k=0,1,2,…,n1）式中k為觀測值編號；是觀測高差；是觀測值的平差改正數(shù)，是稱為殘差；表示高差兩端點的編號（即點號）；分別表示觀測高差起點和終點的高程平差值，即平差中未知數(shù)。 i++) { fprintf(resultfp,\n %5s ,Pname[i])。 fprintf(resultfp,% %%%%, L[i],V[i],L[i]+V[i],P[i],ml)。 fprintf(resultfp,\n%s\n,\n 擬穩(wěn)點：)。i++) { for(int j=0。 PrintM2(resultfp,ATPA,m_Pnumber,4,%15e ,Qx)。 fclose(fp)?！皵M穩(wěn)平差”的基本思想是，考慮到監(jiān)測網(wǎng)中的點，處于不同的地質(zhì)構(gòu)造和地球物理環(huán)境，隨著時間的延伸，都可能發(fā)生變動，但是總存在相對變化小的，即相對穩(wěn)定的點。本程序包含的C++庫函數(shù)頭文件如下：include include include include include include include ，把（）中的函數(shù)源代碼按照順序全部復(fù)制粘貼到這個源代碼文件中即可。余學(xué)祥老師能及時給予我建議，指出不足之處，幫助我更好地組織論文結(jié)構(gòu)，不斷充實論文的內(nèi)容?！啊敝?，其中擬穩(wěn)點選取按照（）中的選擇方法和判斷原則來選取，因為擬穩(wěn)點的選取是不固定的，故本程序在此不易討論。文件中包含擬穩(wěn)點數(shù)和全部擬穩(wěn)點的點名，平差程序可以從這個文件中讀取到平差所需的擬穩(wě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)點的總個數(shù)，第2行以下是擬穩(wěn)點的點名，點名之間用空格隔開。 =fp。 m_mu= sqrt(pvv/(m_Lnumber(m_Pnumber1)))。 } fclose(fp)。 } IsStable=new int[m_Pnumber]。 double qjj= ATPA[ij(k2,k2)] 。 函數(shù)的源代碼如下：void CLevelingAdjust:: PrintResult(){ fprintf(resultfp,\n\n ==== 高程平差值及其精度 ====\n)。 } dX[i]=xi。 ATPL[i]=Pk*Lk。 exit(0)。 //高差起點號 int k2=EndP[i]。i=m_Lnumber1。L[i],amp。Height[c])。 EndP=new int [m_Lnumber]。 //Pname數(shù)組已經(jīng)存滿，且沒有待查點名 }11. Inputdata原始數(shù)據(jù)輸入函數(shù)該函數(shù)有一個形參datafile,定義為string型地址變量，其內(nèi)容是數(shù)據(jù)文件名稱字符串的地址，文件名包括文件的全路徑及擴(kuò)張名，如文件在程序所在的文件目錄下的路徑為“”。 delete []ATPA。 //法方程組成函數(shù) void ca_dX()。 //文件指針,輸出計算結(jié)果 int *StartP。 for(int s=0。ir。 } for(int i=1。函數(shù)將此方程組的系數(shù)矩陣和常數(shù)項向量輸出至指定文件，函數(shù)原型如下：void PrintEquation(FILE* fp, double A[], double b[], int n, int t, char *fmt, char* title)fp——文件指針A——方程組系數(shù)矩陣，數(shù)組長度為b——方程組常數(shù)項向量，數(shù)組長度為n。 } fprintf(fp,\n)。 vsprintf(buffer, fmt, argptr)。3. 按步站形式的不同分為：（1）網(wǎng)狀觀測站：在產(chǎn)狀復(fù)雜的煤層或在建筑物密集的地區(qū)開采時，可考慮多布設(shè)一下測點，組成網(wǎng)格狀觀測站。一般是當(dāng)回采工作面自開切眼開始向前推進(jìn)的距離相當(dāng)于1/4~1/2H0（H0 為平均采深）時，開采影響即波及到地表，引起地表下沉。它不僅是由各種不同性質(zhì)的巖層組成，而且還由于各種地質(zhì)作用（如褶皺、斷層、開裂、火成巖、侵入，陷落柱等）而產(chǎn)生了大量的不連續(xù)面。這些采空區(qū)對國民經(jīng)濟(jì)的發(fā)展有著不可小覷的阻礙。當(dāng)采深和采厚的比值較小（一般小于30）或具有較大的地質(zhì)構(gòu)造時，地表的移動和變形在空間和時間上將是不連續(xù)的，移動和變形的分布沒有嚴(yán)格的規(guī)律性，地表可能出現(xiàn)較大的裂縫或坍塌坑。所謂觀測站，是指在開采影響范圍內(nèi)的地表、巖層內(nèi)部或其它研究對象上，按一定要求設(shè)置的一系列互相聯(lián)系的觀測點。擬穩(wěn)平差將全部高程點分為擬穩(wěn)點和非擬穩(wěn)點兩部分，僅在擬穩(wěn)點參數(shù)的平方和最小的約束下求法方程式的解。 int j=0。j=i。 }}從誤差方程到權(quán)矩陣，到最后的求其平差值和改正數(shù)，都離不開矩陣的求逆運算，函數(shù)原型：bool inverse(double a[],int n)a——函數(shù)調(diào)用前為待求逆矩陣的元素；n——矩陣的階數(shù)；函數(shù)返回值——若計算成功返回true；若計算失敗返回false。in。 N[ij(i,j)]+=nij。 virtual ~CLevelingAdjust()。 //殘差數(shù)組 double m_mu。}析構(gòu)函數(shù)是在對象撤銷時調(diào)用的，利用析構(gòu)函數(shù)刪除對象中由指針成員所指向的動態(tài)分配的存儲空間。 im_Pnumber。m_knPnumber)。i++) { // GetStationNumber函數(shù)根據(jù)Pname[i]是否為NULL // 確定Pname[i]是否為點名地址 Pname[i] = NULL。 } fscanf(fp,%s,buffer)。 fprintf(resultfp,\n\n 已知高程:\n)。im_Pnumber。 } } if(jj==(m_Pnumberm_knPnumber))break。 i++) ATPL[i]=。 exit(0)。i++) { int k1=StartP[i]。 fprintf(resultfp,\n No. 起 點 終 點 觀測高差 ｖ)。 fclose(resultfp)。 int k=GetStationNumber(name)。 } ca_dX()。 printf(\n\n擬穩(wěn)平差計算\n)。下面的數(shù)據(jù)文件內(nèi)容就是本程序的原始數(shù)據(jù)文件部分內(nèi)容的格式。假設(shè)檢驗用于擬穩(wěn)點選擇的具體步驟如下：1)由兩期擬穩(wěn)平差結(jié)果對監(jiān)測網(wǎng)中的各點進(jìn)行單個檢驗，確定出變形點和擬穩(wěn)點。由mian主函數(shù)的建立的文件路徑可知，程序運行的結(jié)果保存在新建工程的文件目錄下，保存在””文件中，””文件運行的結(jié)果和顯示器上顯示的結(jié)果分別如圖（52）和圖（53） 圖52 ””文件運行的結(jié)果 圖53 顯示器顯示結(jié)果結(jié)論自周江文教授將擬穩(wěn)平差的基本原理提出來以后，許多測量工作者將其運用到實踐中，讓擬穩(wěn)平差在沉降監(jiān)測網(wǎng)的數(shù)據(jù)處理的應(yīng)用中，取得了重大的進(jìn)步，具有重要的意義。本文的擬穩(wěn)平差程序設(shè)計借鑒了大量的圖書資料，秉承著“創(chuàng)新為主、借鑒為輔”的原則，成功的實現(xiàn)了從程序源代碼編寫到實現(xiàn)程序運行。三條判別原則：1)擬穩(wěn)平差中所選取的點必須是相對穩(wěn)定點：2)區(qū)分出的變形點應(yīng)處在不穩(wěn)定地帶，應(yīng)符合該區(qū)的形變特征；3)擬穩(wěn)點所構(gòu)成的區(qū)域必須是相對穩(wěn)定區(qū)域。驗前單位權(quán)中誤差即觀測高差的每千米高差中誤差，以米為單位。 printf(結(jié)果文件： %s\n,resultfile)。im_Pnumber。 fprintf(resultfp,\n出現(xiàn)了數(shù)據(jù)文件中沒有的點名：%s,name)。 } fscanf(fp,%d,amp。 for(i=0。 V[i]=Height[k2]Height[k1]L[i]。 im_Pnumber。 km_Lnumber。 for(i=0。 int jj=0。i=m_knPnumber1。 if(EndP[i]0) { fprintf(resultfp,\n數(shù)據(jù)文件出錯：)。 //臨時數(shù)組，保存從文件中讀到的點名 // 讀取已知高程數(shù)據(jù) for( i=0。 Height=new double [m_Pnumber]。 } else { //待查點名是一