【正文】 )。 }int NumberOfVertices()//返回當(dāng)前頂點(diǎn)數(shù) { return numVertices。//中間權(quán)值}。實(shí)驗(yàn)內(nèi)容： 。 coutendl。 cout（以作為分支結(jié)尾，例如：A B C ）endl。 if (rchild!=NULL) thisrchildPreOrder()。 else thislchild= newChild0。 char ch。//中序遍歷 void PostOrder()。中國(guó)礦業(yè)大學(xué)測(cè)繪軟件設(shè)計(jì)與實(shí)現(xiàn)實(shí)驗(yàn)報(bào)告學(xué)　　號(hào)：　 　　 姓　　名：　 　　 班　　級(jí)：　 指導(dǎo)教師：　　王永波　　實(shí)驗(yàn)一 二叉樹(shù)的構(gòu)建及其遍歷算法的實(shí)現(xiàn)實(shí)驗(yàn)?zāi)康模和瓿啥鏄?shù)的構(gòu)建以及二叉樹(shù)的遍歷等，加深對(duì)樹(shù)以及二叉樹(shù)的遍歷相關(guān)知識(shí)的理解。//后序遍歷 private: T data。 cinch。 //右子樹(shù) C_LJH_BinTree *newChild1 = new C_LJH_BinTreeT()。 } //中序遍歷 template class T void C_LJH_BinTreeT::InOrder() { if (lchild) lchildInOrder()。 C_LJH_BinTreechar m_tree。 cout后序遍歷為：endl。 。class LJH_Graphmtx//圖的鄰接矩陣定義{ public: LJH_Graphmtx(int sz=DefaultVertices)。 } int NumberOfEdges()//返回當(dāng)前邊數(shù) { return numEdges。//取頂點(diǎn)v的第一個(gè)鄰接頂點(diǎn) int getNextNeighbor(int v,int w)。 return 1。//圖中當(dāng)前邊數(shù) int numVertices。 numEdges=0。imaxVertices。Edge[v][i]maxWeight)return i。Edge[v][i]maxWeight) return i。amp。 } else return false。i++)//減去與v相關(guān)聯(lián)的邊數(shù) if(Edge[i][v]0amp。jnumVertices。amp。 } else return false。in。 if(j==1||k==1) cout邊兩端點(diǎn)信息輸入有誤，請(qǐng)重新輸入！endl。 m=thisNumberOfEdges()。j++) { w=thisgetWeight(i,j)。 int min=0。 G,int v,bool visited[])。 } loc=(v)。//若w未被訪問(wèn)，遞歸訪問(wèn)頂點(diǎn)w w=(v,w)。//取頂點(diǎn)v的位置號(hào)cout(loc)endl。//標(biāo)記w已經(jīng)被訪問(wèn)(w)。m++) visit[m]=m。inumVertices。 for (int n=0。 int i=thisgetVertexPos(u)。mnumVertices。} } } }}實(shí)驗(yàn)結(jié)果：實(shí)驗(yàn)體會(huì)：經(jīng)過(guò)這次實(shí)驗(yàn)讓我更深刻的理解了類的創(chuàng)建及相互間調(diào)用，能夠?qū)ΧS數(shù)組的動(dòng)態(tài)開(kāi)辟空間和釋放空間有了更深刻的理解，對(duì)圖的遍歷及構(gòu)建最小生成樹(shù)也有了深刻的體會(huì)。 m)。//求矩陣的伴隨矩陣 LJH_CMatrix inverse()。 } int getRow() const { return _row。 cout請(qǐng)輸入矩陣的行數(shù)：endl。j_column。 } coutendl。j_column。m_row。 } if(m!=i) k++。 for(i=0。 for(i=0。 n=thishanglieshi()。 } return temp。 _A[l][i]=j。 } } return m。 =。i++) [i][i]=1。 for(m=i+1。i0。 [m][j]=[m][j]xishu1*[i][j]。 [i][j]*=xishu3。 ()。 ()。 coutendl。通過(guò)構(gòu)造一個(gè)跟需求矩陣同階的單位矩陣，然后通過(guò)行變換把該矩陣變成單位陣，原來(lái)的單位陣也做相應(yīng)變化，便得到所求的逆矩陣。主要代碼：void Matrix::ColPivot(Matrixamp。j++) { if(fabs([k][i])fabs([j][i])) k=j。t+1。j=i+1。}void Matrix::DirectLU(Matrixamp。 } for(i=0。i++) //得到b矩陣 { b[i]=[i][]。t++) { [j][t] = L[j][i]*[i][t]。 } coutU矩陣:endl。 for (i=0。 for (i=0。i) //利用公式Ux=Y，求解得到x { for(j=。 coutendl。 while (true) { epsilon=0。 } Y[i] = ([i][]s)/[i][i]。i++) [i]= Y[i]。 double epsilon,s,*Y,max,temp。 for (i=0。 } max=fabs([0]Y[0])。i。i++) [i]= Y[i]。numPoints=0。} //返回觀測(cè)值數(shù)目 int getzds(){return numPoints。} //返回改正數(shù) friend void xishu(LJH_CMatrixamp。 //求取觀測(cè)值L矩陣和l陣 public: // 成員變量 int numgaocha。 cout請(qǐng)輸入水準(zhǔn)網(wǎng)中水準(zhǔn)點(diǎn)總數(shù)和已知點(diǎn)數(shù)：endl。 gczhi[i].eleValue=0。 } cout請(qǐng)輸入觀測(cè)值的個(gè)數(shù)：endl。 thisedVec[i].weight=n。ithisnumgaocha。 }else if(thisgczhi[thisedVec[i].startPoint1].eleValue==0amp。amp。 for(i=0。amp。 [i][[i].]=1。 for(i=0。j++) [i][j]=0。 for(i=0。首先通過(guò)本實(shí)驗(yàn)，我發(fā)現(xiàn)雖然自己以前平差課學(xué)的很好，但真正對(duì)應(yīng)到編程實(shí)現(xiàn)時(shí)還會(huì)存在很多問(wèn)題。通過(guò)查資料，以及參考往年水準(zhǔn)網(wǎng)平差設(shè)計(jì)的程序，漸漸的理清了思路。i++) { [i][0]=[i].value。i。i++) sum+=[i].weight。i。 }//如果起始點(diǎn)為未知點(diǎn)并且結(jié)束點(diǎn)為已知點(diǎn) 則系數(shù)設(shè)置如下 else if([[i].startPoint1].isKnown==0amp。i++) for(j=0。 }//如果觀測(cè)段的結(jié)束點(diǎn)等于水準(zhǔn)點(diǎn)的編號(hào) 并且起始點(diǎn)點(diǎn)得高程不為0則該該結(jié)束點(diǎn)的高程等于起始點(diǎn)高程加上高差值（觀測(cè)值） else if(thisedVec[i].endPoint==thisgczhi[i].bhamp。thisgczhi[thisedVec[i].endPoint1].eleValue!=0) { thisgczhi[thisedVec[i].startPoint1].eleValue=thisgczhi[thisedVec[i].endPoint1].eleValuethisedVec[i].value。 } void LJH_shzwNet::jsgc()//求每個(gè)點(diǎn)得近似高程 { int i,j。 thisedVec[i].endPoint=b。 thisnumgaocha=c。 } for(i=0。 thisnumPoints=a。 // 水準(zhǔn)網(wǎng)中點(diǎn)的總數(shù)目 int numKnPoint。 X,LJH_shzwNet A)。} //返回已知點(diǎn)數(shù) int getwzds(){return numPointsnumKnPoint。}。明白了，只要迭代矩陣的譜半徑小于1，迭代就是收斂的，與右邊所取得的近似值無(wú)關(guān)。 coutendl。i。i++) { s=0。 for(i=0。 for(i=0。 } if (epsilon1e6) { cout迭代次數(shù)為：kendl。 for (i=0。 m) //Jacobi雅克比迭代法{ int i,j,k=0。j) Y[i] = [j]*[i][j]。i++) { for(j=0。i++) { coutb[i]endl。i。 for(i=0。i。i++) for (j=0。 double **L,*b。 [i]= [i][]/[i][i]。 } } } if ([][]==0amp。j+1。 double l。經(jīng)過(guò)這次實(shí)驗(yàn)，獲益匪淺。 ()。 array7=array1array2。 ()。 }}}void main(){ LJH_CMatrix array1,array2,array3,array4,array5,array6,array7。 for(i=0。 for(m=0。m++)//第m行 { xishu=[m][i]/[i][i]。 for(i=0。i。 LJH_CMatrix E1。 } int LJH_CMatrix::max_cloumn(int k)//求矩陣第k列中從第k個(gè)元素之后絕對(duì)值最大的行數(shù) { int m=k。 LJH_CMatrix::change(int k,int l)//交換矩陣的第k行和第l行 { int i。 int i,j。i++) for(j=0。i++) { sum+=pow(1,i)*_A[0][i]*thisyuzishi(0,i).hanglieshi()。 } double LJH_CMatrix::hanglieshi()//求矩陣的行列式 { if(_row!=_column) {cerr此矩陣無(wú)行列式endl。 for(n=0