【正文】 即使借鑒的過程，也是自己理解、學習的過程，通過本次實驗，我切實的體會到了理論聯(lián)系實際需要做的東西很多。 } }實驗結果:實驗體會：做完本實驗收獲頗多。 } for(i=0。 double m。i++) [i][i]=avg/[i].weight。j。//水準路線的總長度 avg=sum/。 double avg。i++) {//未知點的近似高程矩陣 [i][0]=[+i].eleValue。[[i].endPoint1].isKnown==0) { [i][[i].]=1。amp。i++) {//如果起始點為已知點并且結束點為未知點 則系數(shù)設置如下 if([[i].startPoint1].isKnown==1amp。j()。 X,LJH_shzwNet A)//友元函數(shù),求誤差方程得系數(shù)矩陣 { int i,j。amp。i++) {//如果觀測段的起始點等于水準點的編號 并且結束點得高程不為0則該該起始點的高程等于結束點高程減去高差值（觀測值） if(thisedVec[i].startPoint==thisgczhi[i].bhamp。 } } for(j=thisnumKnPoint。thisgczhi[thisedVec[i].endPoint1].eleValue==0) { thisgczhi[thisedVec[i].endPoint1].eleValue=thisgczhi[thisedVec[i].startPoint1].eleValue+thisedVec[i].value。 for(i=0。 for(i=0。 } } void LJH_shzwNet::output() { int i。 thisedVec[i].value=k。 for(i=0。 gczhi[i].isKnown=1。ib。 gczhi[i].dv=0。 thisnumKnPoint=b。 double m,n,k。 //水準網(wǎng)中已知點的數(shù)目 LJH_gaocha edVec[maxnum]。 L, LJH_shzwNet A)。 //求取系數(shù)矩陣和未知點高程矩陣 friend void quanzhen(LJH_CMatrixamp。} //修改改正數(shù) double getdv(int i){return gczhi[i1].dv。} //返回未知點數(shù) double getgcValue(int i){return thisgczhi[i1].eleValue。////求每個點得近似高程 int getgczs(){return numgaocha。 // 構造函數(shù) ~LJH_shzwNet(){}。主要代碼：class LJH_shzwNet { public: // 成員函數(shù) LJH_shzwNet(){numgaocha=0。實驗五、面向對象的水準網(wǎng)平差程序設計與實現(xiàn)實驗目的：在矩陣實現(xiàn)的基礎上，實現(xiàn)對水準網(wǎng)平差。i。 break。 for(i=0。i++) { temp=fabs([i]Y[i])。 //epsilon += fabs(Y[i][i])。 for (j=0。 k++。i。 m) //GaussSeidel高斯賽德爾迭代法{ int i,j,k=0。i。i。 break。j++) { if(i!=j) s += [i][j]*[j]。i。 Y = new double []。 double epsilon,s,*Y。i++) coutx[i+1]=[i] 。 [i]= Y[i]/[i][i]。i=0。j=i。i++) Y[i]=0。 } double *Y。 } coutb矩陣：endl。i++) { for (j=0。 } coutendl。i。t。i++) { for (j=i+1。i。j。i++) { L[i] = new double []。 b = new double []。 coutendl。 } cout用列主元素法求得的方程組的解：endl。i) //自下往上逐步回代求得真解 { for(j=。amp。 for(t=i。j++) { swap([i][j],[k][j])。j。 for (i=0。同時掌握這些函數(shù)的實現(xiàn)。實驗四、基于直接法（列主元素法）的線性方程組求解實驗目的：基于直接法（列主元素法）與迭代法（Jacobi迭代與GaussSeidel迭代法），和收斂速度做個比較。比如在實現(xiàn)矩陣的轉置和矩陣的求逆這方面，涉及遞歸調(diào)用。 coutendl。 ()。 ()。 array6=array1*array2。 cout兩矩陣求和：endl。 cout請輸入矩陣array2相關信息：。 cout請輸入矩陣array1相關信息：。j++) { [i][j]*=xishu3。i。j++) { [m][j]=[m][j]xishu1*[i][j]。mi。}} ////上面為求上三角陣 } ///////////////////////////化對角陣 for(i=this_column1。 for(j=0。 } double xishu。i。i。i++) for(j=0。//創(chuàng)建一個和當前方陣階數(shù)相同的單位矩陣 =。 E1=*this。 m=i。 double max=fabs(_A[k][k])。 _A[k][i]=_A[l][i]。 double j。j++) { [i][j]/=n。 for(i=0。 int n。j_column。 int i,j。 } return sum。 double sum=0。 } if(_row==1amp。 } if(n!=j) l++。n_column。 for(m=0。 } LJH_CMatrix LJH_CMatrix::yuzishi(int i,int j)//求矩陣的余子式 { LJH_CMatrix temp。i++) for(j=0。 =this_row。j++) { cout_A[i][j] 。 for(i=0。i++) for(j=0。 cin_column。 void LJH_CMatrix::input()//輸入函數(shù) { int i,j。 }public:// 成員變量 double** _A。 } // 設置行、列的值 void setRow(const int row) { _row = row。 change(int k,int l)。//求矩陣的行列式 LJH_CMatrix bansui()。//矩陣輸入 void output()。 // 構造函數(shù)一 LJH_CMatrix(const LJH_CMatrixamp。同時加深對矩陣的理論的理解，并用計算機程序算法來描述矩陣生成及運算。 closeedge[p].begvex=VerticesList[n]。 closeedge[n].lowcost=maxWeight。 } for (int m=1。jnumVertices。 cinu。 } } else Edge1[a][b]=Edge[b][a]=maxWeight。 k++。j++) { if (Edge1[i][j]min) { a=i。 for (int i=0。i++) for(int j=0。mnumVertices。}//////kruskal函數(shù)的實現(xiàn)//////void LJH_Graphmtx::kruskal(){ int a,b,k=0。//訪問頂點wvisited[w]=1。//定義一個輔助隊列(loc)。int loc=(v)。//去圖中的定點個數(shù) bool *visited=new bool[n]。//找v的第一個鄰接頂點w while(w!=1)//若鄰接頂點w存在 { if(visited[w]==0) DFS(G,w,visited)。}void DFS(LJH_Graphmtx amp。i++)//初始化輔助數(shù)組visited { visited[i]=0。 v)//從頂點v出發(fā)，對圖G進行深度優(yōu)先遍歷的主要過程{ int i,loc,n=()。}//圖的深度優(yōu)先搜索函數(shù)////////void DFS(LJH_Graphmtx amp。j++) { if(!visit[j]){if(closeedge[min].lowcostcloseedge[j].lowcost) { min=j。}int LJH_Graphmtx::mini()//求圖中所有邊的最小權值，并返回{ static int i。wmaxWeight) { e1=thisgetValue(i)。jn。 cout所有邊的信息為：endl。 n=thisNumberOfVertices()。 }} return true。//查頂點號 k=thisgetVertexPos(e2)。 } i=0。 for(i=0。 int weight。 return true。Edge[v1][v2]0amp。v21amp。}bool LJH_Graphmtx::removeEdge(int v1,int v2)//在圖中刪去邊(v1,v2){if(v11amp。//頂點個數(shù)減1 for(j=0。 for(i=0。inumVertices。v=numVertices)return false。 return true。v2numVerticesamp。v1numVerticesamp。//頂點表滿，不插入 VerticesList[numVertices++]=vertex。amp。amp。amp。j++) Edge[i][j]=(i==j)? 0 : maxWeight。for(i=0。//創(chuàng)建頂點表數(shù)組 Edge=(int * *)new int *[maxVert