【正文】 如下：double Calculate_q(double Q[],double F[],int Fin[],int n){ double q=。 for(int k=0。 } delete []a0。 else a0[i]= ai0/a00。kn。i++) { fprintf(fp,\n%3d ,i+1)。 j%t==0)fprintf(fp,\n )。 for(int i=0。i++) { if(i%t==0) { j++。 getchar()。：應(yīng)用程序都需要通過計算機界面動態(tài)地顯示一些提示性信息，例如，正定矩陣求逆計算遇到被求逆的矩陣是非正定矩陣，程序就無法繼續(xù)運行，應(yīng)該通過用戶界面顯示“矩陣求逆失敗”的提示信息，然后終止程序運行。沿走向主斷面布設(shè)的觀測線稱為走向觀測線，沿傾斜主斷面布設(shè)的觀測線稱為傾斜觀測線。主要研究地表移動和變形的規(guī)律；（2）巖層內(nèi)部觀測站：測點一般布設(shè)在井下巷道或巖層內(nèi)部的鉆孔中，用于研究巖層內(nèi)部的移動和變形規(guī)律；（3）專門觀測站：為了某一個特定的目的所設(shè)立的觀測站，如建筑物觀測站、鐵路觀測站、邊坡移動觀測站等。巖層與地表移動的過程十分復(fù)雜，它是許多地質(zhì)采礦因素綜合影響的結(jié)果。在開采影響波及到地表以后，受采動影響的地表從原有標(biāo)高向下沉降，從而在采空區(qū)上方地表形成一個比采空區(qū)面積還大的沉陷區(qū)域。沉降監(jiān)測網(wǎng)是為了監(jiān)測“地表移動”的表面上指定位置的點的沉降變化而建立的一種水準(zhǔn)網(wǎng)。本文正是基于這一理論的基礎(chǔ)上，研究擬穩(wěn)平差在沉降觀測網(wǎng)中的應(yīng)用以及程序設(shè)計。沉降監(jiān)測網(wǎng)擬穩(wěn)平差及程序設(shè)計畢業(yè)論文目錄ABSTRACT 31沉降監(jiān)測網(wǎng)概述 6 沉降監(jiān)測網(wǎng)的形成 62地下開采引起的地表移動和破壞 6 6 73沉降監(jiān)測網(wǎng)的數(shù)據(jù)來源 7 7 84擬穩(wěn)平差程序概述 9 9 10 11 30 31 31 31 32 33 33 34 35緒論中國的經(jīng)濟發(fā)展隨著改革開放而騰飛，而經(jīng)濟的發(fā)展離不開對各種能源的需求。 1沉降監(jiān)測網(wǎng)概述 沉降監(jiān)測網(wǎng)的形成有用礦物被采出以后，開采區(qū)域周圍的巖體的原始應(yīng)力平衡狀態(tài)受到破壞，應(yīng)力重新分布，到達新的平衡。通過實地多次觀測這些指定位置的點的高程和平面坐標(biāo)，確定這些點的空間位置坐標(biāo)，進行數(shù)據(jù)統(tǒng)計分析，觀察這些點的沉降情況、沉降速度、沉降變化等，判斷并預(yù)測這些點位在空間中的變化趨勢。這種地表沉陷區(qū)稱為地表移動盆地，或稱下沉盆地。認(rèn)識巖層與地表移動這一復(fù)雜過程，目前的主要方法是實地觀測。2. 普通觀測站和短期觀測站普通觀測站的觀測時間較長（一般在一年以上），它是在地表移動的開始到移動結(jié)束的整個過程中定期進行觀測，主要用于研究地表移動和變形的規(guī)律。 圖31觀測站的布置形式示意圖本程序處理的原始數(shù)據(jù)包括各個點的平面坐標(biāo)和各個點之間的觀測高差，其中點位的平面坐標(biāo)用于計算水準(zhǔn)路線的長度，用以定權(quán)，原始數(shù)據(jù)以excel表格形式如下圖（32）和（33）所示 圖32 點位的平面坐標(biāo) 圖33 各點的觀測高差4擬穩(wěn)平差程序概述自由網(wǎng)擬穩(wěn)平差公式設(shè)誤差方程式為 （式41）式中：為n維自由項向量；V是n維殘差向量；A為系數(shù)矩陣；X是t維參數(shù)向量，這里X是高程點近似值的平差改正數(shù)。函數(shù)源代碼如下：void MyBreak(char* fmt, ...){ char buffer[256]。endif // VC_EXTRALEAN}平差程序經(jīng)常會處理對稱矩陣，例如觀測值的權(quán)矩陣、權(quán)逆陣、法方程系數(shù)陣等都是對稱矩陣，為了節(jié)省存儲空間和避免重復(fù)計算，采用僅存下三角矩陣（含主對角線元素）的存儲方案存儲對稱矩陣，亦即將對稱矩陣的下三角元素按順序存到數(shù)組中。 if(IsLabel)fprintf(fp,\n%3d ,j)。in。 fprintf(fp,fmt,M[index++])。 for(int j=0。k++) { double a00=a[0]。 for(int j=1。 return true。kn。 for(int k=0。}9. CLevelingAdjust類設(shè)計C語言中的類是指對具有共同屬性的一類事物的抽象描述，共同屬性被描述為類中的數(shù)據(jù)成員，共同行為被描述為類中的成員函數(shù)。 //已知點數(shù) int m_StablePnumber。 //高程值數(shù)組 double *P。 //輸出原始數(shù)據(jù)函數(shù) int GetStationNumber(char *name)。構(gòu)造函數(shù)的作用是為類對象中的數(shù)據(jù)成員賦初值，這里指的是為m_Pnumber（總點數(shù)）和m_Lnumber（觀測值）總數(shù)賦初值為0。 delete []P。i++) if(Pname[i]!=NULL)delete[](Pname[i])。 strcpy(Pname[i], name)。 exit(0)。 ATPA=new double [m_Pnumber*(m_Pnumber+1)/2]。m_Sigma)。i++) { fscanf(fp,%s,buffer)。 if(StartP[i]0) { fprintf(resultfp,\n數(shù)據(jù)文件出錯：)。 fclose(resultfp)。在結(jié)果文件中顯示網(wǎng)額信息，例如；網(wǎng)的總點數(shù)，總觀測值，已知點個數(shù)。 fprintf(resultfp,% ,Height[i])。當(dāng)某點的近似高程計算出來之后，近似高程計算就此結(jié)束。ii++) { for(i=0。 } else if(Height[k1] amp。i++) { if(Height[i]) printf(\n%s,Pname[i])。 it*(t+1)/2。 int j=EndP[k]。 }}15. ca_dX高程平差值計算函數(shù)該函數(shù)首先調(diào)用對稱正定矩陣求逆函數(shù)inverse，計算法方程系數(shù)矩陣的逆陣，然后以總點數(shù)為循環(huán)界逐點計算各點高程的改正數(shù)，將高程改正數(shù)保存dX數(shù)組中，同時將高程改正數(shù)累加到數(shù)組Height的對應(yīng)單元中。 for(int j=0。第k次循環(huán)中所要進行的工作包括：1) 獲得高差的起點號i和終點j；2) 獲得起點和終點的高程值（已經(jīng)是平差值）；3) 計算殘差，存在V[k]中；4) 將累加到[pvv]中。 } return(pvv)。 double dx=dX[i]。i++) { int k1=StartP[i]。 }}18. Quasi_Stable擬穩(wěn)平差函數(shù)擬穩(wěn)平差的法方程系數(shù)矩陣為（），可先調(diào)用ca_ATPA函數(shù)計算，再向數(shù)組累加。 if(m_StablePnumberm_Pnumber || m_StablePnumber1) { MyBreak(擬穩(wěn)點數(shù)錯誤!)。 for( i=0。 exit(0)。j=i。j=i。 }主函數(shù)的功能就是調(diào)用所有其它非主函數(shù)，并最終在主函數(shù)的調(diào)用下，執(zhí)行這個程序。 if (fp==NULL) { MyBreak(創(chuàng)建結(jié)果文件失??！)。}水準(zhǔn)網(wǎng)平差所需的已知數(shù)據(jù)和觀測數(shù)據(jù)全部從數(shù)據(jù)文件中讀取，所以將水準(zhǔn)網(wǎng)原始數(shù)據(jù)文件的內(nèi)容分為三類：網(wǎng)的概況、已知數(shù)據(jù)、觀測數(shù)據(jù)。（3）第4行至第10行為觀測高差，每行的內(nèi)容為高差起始點點名、高差終點點名、高差觀測值和高差觀測路線長度，高差以米為單位，路線長度以千米為單位。我們稱之為“擬穩(wěn)點”。所以我將第一觀測的數(shù)據(jù)與第二次觀測的數(shù)據(jù)、第三次觀測的數(shù)據(jù)做了相對比較，即用第一次的觀測數(shù)據(jù)減去第二次觀測的數(shù)據(jù)、第三次的觀測數(shù)據(jù)，從中挑選出高差變化較小的點。然后單機菜單欄下面的“組建”，從其下拉列表中單擊“編譯”。在我進行畢業(yè)論文的收集資料，初稿，最后定稿的各個階段都給予了我極大的幫助，尤其是當(dāng)我遇到一些困難請教老師時，他都會不予余力的給予我最大的幫助以及指引。余學(xué)祥老師一直是我非常敬佩和尊重的教授，不僅是因為他在學(xué)術(shù)上的成就，更由于他對教學(xué)及科研的態(tài)度?！敖M建”，從其下拉菜單中單擊“執(zhí)行”即可。7) 通過表格最后三列的高差之差數(shù)據(jù)，可以較容易的判斷和選取高差之差變化相對較小的點，最后結(jié)合點位的平面分布圖，選擇擬穩(wěn)點8) 表格中沒有數(shù)據(jù)的單元格，表明這三次的測站不是同一測站，因為只有在同一測站的情況下，才能判斷點位的高差變化，通過點位的平面坐標(biāo)，使用Excel表格，計算出各個測站之間的距離，計算結(jié)果如下圖（51）： 圖51 點位坐標(biāo)與測站對應(yīng)圖表格中最后一列為測站間的平面距離，將表格中最后四列信息復(fù)制粘貼到“”文件中，放到第四行即可觀測數(shù)據(jù)的高差和距離的輸入工作。擬穩(wěn)平差的關(guān)鍵是選取符合實際的“擬穩(wěn)點”，即確定陣。為了進行擬穩(wěn)平差，首先要知道哪些點是擬穩(wěn)點，哪些點不是擬穩(wěn)點。已知數(shù)據(jù)包括：已知點名、已知高程值。 } CLevelingAdjust level。 char *resultfile =。 } double pvv = ca_V()。amp。 fprintf(resultfp,\n%10s ,Pname[k])。i++) { char name[100]。 exit(0)。函數(shù)的源代碼如下：void CLevelingAdjust:: Quasi_Stable(char *file) //file是擬穩(wěn)點名文件{ //打開文件，獲取擬穩(wěn)點名 FILE *fp=fopen(file,r)。 double qii=ATPA[ij(k1,k1)]。 fprintf(resultfp,%%%%, Height[i]dx,dx,Height[i],sqrt(qii)*m_mu)。輸出的內(nèi)容包括高程點名、高程平差值及其中誤差，高差端點點名