【正文】 算的思路第一步，計算位置點標(biāo)志。i++) { fprintf(resultfp,\n%5d%8s%8s,i+1,Pname[StartP[i]],Pname[EndP[i]])。 for(i=0。 fprintf(resultfp,% ,Height[i])。i=m_knPnumber1。 fprintf(resultfp,\n\n 已知高程:\n)。 fprintf(resultfp,\n 觀測值總數(shù): %d 總點數(shù): %d 已知點數(shù)：%d \n, m_Lnumber, m_Pnumber,m_knPnumber)。在結(jié)果文件中顯示網(wǎng)額信息，例如；網(wǎng)的總點數(shù)，總觀測值，已知點個數(shù)。 } fclose(fp)。P[i])。 } fscanf(fp,%lf%lf,amp。 fclose(resultfp)。 if(EndP[i]0) { fprintf(resultfp,\n數(shù)據(jù)文件出錯：)。 } fscanf(fp,%s,buffer)。 fclose(resultfp)。 if(StartP[i]0) { fprintf(resultfp,\n數(shù)據(jù)文件出錯：)。i++) { fscanf(fp,%s,buffer)。 } // 讀取觀測數(shù)據(jù) for(i=0。 fscanf(fp,%lf,amp。i++) { fscanf(fp,%s,buffer)。 //臨時數(shù)組，保存從文件中讀到的點名 // 讀取已知高程數(shù)據(jù) for( i=0。i++) { // GetStationNumber函數(shù)根據(jù)Pname[i]是否為NULL // 確定Pname[i]是否為點名地址 Pname[i] = NULL。 for(int i=0。m_Sigma)。 P=new double [m_Lnumber]。 L=new double [m_Lnumber]。 StartP=new int [m_Lnumber]。 ATPA=new double [m_Pnumber*(m_Pnumber+1)/2]。 Height=new double [m_Pnumber]。m_knPnumber)。m_Lnumber,amp。 exit(0)。函數(shù)院代碼如下：void CLevelingAdjust:: Inputdata(char *datafile){ FILE *fp。為了從數(shù)據(jù)文件讀取數(shù)據(jù)，函數(shù)內(nèi)首先聲明了FILE(文件指針定義的關(guān)鍵字)型變量fp，并調(diào)用fopen函數(shù)打開數(shù)據(jù)文件，將文件的指針賦給fp。 } } return 1。 strcpy(Pname[i], name)。 } else { //待查點名是一個新的點名，將新點名的地址放到Pname數(shù)組中 int len = strlen(name)。 im_Pnumber。 }}10. GetStationNumber點名獲取函數(shù)工作流程如下：（1）將待查點名name與Pname數(shù)組中已經(jīng)保存的點名利用系統(tǒng)內(nèi)置函數(shù)strcmp逐一比較，檢查Pname中是否已經(jīng)保存有name這個點名，如果Pname中存在name這個點名，就返回到name在Pname中的下標(biāo)。i++) if(Pname[i]!=NULL)delete[](Pname[i])。 for(int i=0。 delete []ATPL。 } if(m_Pnumber0) { delete []Height。 delete []P。 delete []EndP。}析構(gòu)函數(shù)是在對象撤銷時調(diào)用的，利用析構(gòu)函數(shù)刪除對象中由指針成員所指向的動態(tài)分配的存儲空間。 m_Pnumber=0。構(gòu)造函數(shù)的作用是為類對象中的數(shù)據(jù)成員賦初值，這里指的是為m_Pnumber（總點數(shù)）和m_Lnumber（觀測值）總數(shù)賦初值為0。 //殘差計算函數(shù) void Quasi_Stable( char *file)。 //平差值計算函數(shù) void PrintResult()。 //近似高程計算函數(shù) void ca_ATPA()。 //輸出原始數(shù)據(jù)函數(shù) int GetStationNumber(char *name)。 //是否為擬穩(wěn)點號 void Inputdata(char *datafile)。 //殘差數(shù)組 double m_mu。 //法方程系數(shù)矩陣與自由項 double *dX。 //高程值數(shù)組 double *P。 //點名地址數(shù)組 double *L。 //高差起點號 int *EndP。 //[pvv] FILE *resultfp。 //已知點數(shù) int m_StablePnumber。 //高差總數(shù) int m_Pnumber。 virtual ~CLevelingAdjust()。程序中具有相同屬性的數(shù)據(jù)一般用數(shù)組來保存，例如，高程值、觀測高差、點名等等。}9. CLevelingAdjust類設(shè)計C語言中的類是指對具有共同屬性的一類事物的抽象描述，共同屬性被描述為類中的數(shù)據(jù)成員，共同行為被描述為類中的成員函數(shù)。 q+=Q[ij(i,j)]*F[k]*F[s]。sn。k++) { int i=Fin[k]。 for(int k=0。函數(shù)原型：double Calculate_q(double Q[],double F[],int Fin[],int n)Q——權(quán)逆陣數(shù)組，僅存下三角矩陣，數(shù)組長度等于；F——系數(shù)數(shù)組，數(shù)組長度為t；t——權(quán)逆陣的階數(shù)；返回值——返回式（）計算的結(jié)果。 N[ij(i,j)]+=nij。sn。kn。j++) { double nij=。i++) for(int j=0。函數(shù)源代碼如下：void Calculate_BQBT(double B[],double Q[],int r,int n,double N[]){ for(int i=0。 return true。 } a[n*(n+1)/21]=。in。j++) { a[(i1)*i/2+j1]=a[i*(i+1)/2+j]+ai0*a0[j]。 for(int j=1。 if(i=nk1)a0[i]= ai0/a00。in。 return false。k++) { double a00=a[0]。 for(int k=0。 }}從誤差方程到權(quán)矩陣，到最后的求其平差值和改正數(shù)，都離不開矩陣的求逆運算，函數(shù)原型：bool inverse(double a[],int n)a——函數(shù)調(diào)用前為待求逆矩陣的元素；n——矩陣的階數(shù)；函數(shù)返回值——若計算成功返回true；若計算失敗返回false。j++) { fprintf(fp,fmt,A[i*t+j])。 for(int j=0。in。n——方程的個數(shù)t——方程未知數(shù)個數(shù)fmt——輸出格式title——標(biāo)題字符串地址，默認值為空字符串函數(shù)源代碼如下：void PrintEquation(FILE* fp, double A[], double b[], int n, int t, char *fmt, char* title){ if(title)fprintf(fp,\n %s: ,title)。 }設(shè)有線性方程組，A是矩陣，b是n維向量。 fprintf(fp,fmt,M[index++])。amp。j=i。 else fprintf(fp,\n )。in。 int index=0。}4. PrintM2對稱矩陣輸出函數(shù)設(shè)有階對稱矩陣，程序中僅存矩陣下三角元素，矩陣元素為雙精度型，函數(shù)PrintM2將矩陣輸出至指定文件，函數(shù)原型如下：void PrintM2(FILE* fp, double M[], int n, int t,char *fmt, char* title,bool IsLabel)fp——文件指針M——待輸出的數(shù)組n——矩陣的階數(shù)t——格式控制變量fmt——輸出格式title——標(biāo)題字符串地址，缺省值空IsLabel——值為true時每行前添加行號，等于false時不加行號，默認值為true。 } fprintf(fp,fmt,A[i])。 if(IsLabel)fprintf(fp,\n%3d ,j)。isize。 int j=0。函數(shù)源代碼如下：int ij(int i,int j){ return (i=j)? i*(i+1)/2+j :j*(j+1)/2+i。endif // VC_EXTRALEAN}平差程序經(jīng)常會處理對稱矩陣，例如觀測值的權(quán)矩陣、權(quán)逆陣、法方程系數(shù)陣等都是對稱矩陣，為了節(jié)省存儲空間和避免重復(fù)計算，采用僅存下三角矩陣（含主對角線元素）的存儲方案存儲對稱矩陣，亦即將對稱矩陣的下三角元素按順序存到數(shù)組中。else printf(buffer)。 va_end(argptr)。 va_start(argptr, fmt)。函數(shù)源代碼如下：void MyBreak(char* fmt, ...){ char buffer[256]。設(shè)G、G2都是t維向量， 其中（i=0,1,…,t1）為 i號點位擬穩(wěn)點 （式45） i號點為非擬穩(wěn)點 擬穩(wěn)平差參數(shù)解及其權(quán)逆陣的公式為 （式46）