【正文】 矩陣為：endl。 coutendl。但我收獲很多。比如在實(shí)現(xiàn)矩陣的轉(zhuǎn)置和矩陣的求逆這方面，涉及遞歸調(diào)用。此外，這個(gè)實(shí)驗(yàn)還讓我們復(fù)習(xí)了以前所學(xué)的函數(shù)重載，運(yùn)算符重載等知識。實(shí)驗(yàn)四、基于直接法（列主元素法）的線性方程組求解實(shí)驗(yàn)?zāi)康模夯谥苯臃ǎ兄髟胤ǎ┡c迭代法（Jacobi迭代與GaussSeidel迭代法），和收斂速度做個(gè)比較。實(shí)驗(yàn)內(nèi)容：完成基于直接法（列主元素法）與迭代法（Jacobi迭代與GaussSeidel迭代法）的線性方程組求解并且對各種迭代的收斂條件，和收斂速度做個(gè)比較。同時(shí)掌握這些函數(shù)的實(shí)現(xiàn)。 m) //列主元素法{ int i,j,k,t。 for (i=0。i++) { k=i。j。 } if(k!=i) { for (j=0。j++) { swap([i][j],[k][j])。j。 for(t=i。t++) { [j][t] = l*[i][t]。amp。 } for (i=。i) //自下往上逐步回代求得真解 { for(j=。j) [i][] = [j]*[i][j]。 } cout用列主元素法求得的方程組的解：endl。i。 coutendl。 m) //LU分解法{ int i,j,/*k,*/t。 b = new double []。 //對LU分解中的L矩陣的內(nèi)存分配、初始化 for (i=0。i++) { L[i] = new double []。i。j。 else L[i][j]=0。i。 } for (i=0。i++) { for (j=i+1。j++) //將系數(shù)矩陣化為上三角矩陣得到U矩陣 { L[j][i]=[j][i]/[i][i]。t。 } } } coutL矩陣:endl。i。j。 } coutendl。 for(i=0。i++) { for (j=0。j++) { cout[i][j] 。 } coutb矩陣：endl。i。 } double *Y。 for(i=0。i++) Y[i]=0。i。j=i。 Y[i]= b[i]/L[i][i]。i=0。j=i+1。 [i]= Y[i]/[i][i]。 for(i=0。i++) coutx[i+1]=[i] 。}void Matrix::Jacobi(Matrixamp。 double epsilon,s,*Y。i。 Y = new double []。 k++。i。 for (j=0。j++) { if(i!=j) s += [i][j]*[j]。 epsilon += fabs(Y[i][i])。 break。 break。i。 } cout用Jacobi雅克比迭代法求得的方程組的解：endl。i。 coutendl。 m) //GaussSeidel高斯賽德爾迭代法{ int i,j,k=0。 Y = new double []。i。 Y[i]=0。 k++。i。 for (j=0。j++) { if(j!=i) s += [i][j]*Y[j]。 //epsilon += fabs(Y[i][i])。 for (i=0。i++) { temp=fabs([i]Y[i])。 } if (max1e6) { cout迭代次數(shù)為：kendl。 for(i=0。i++) coutx[i+1]=Y[i] 。 break。 break。i。 }實(shí)驗(yàn)結(jié)果:實(shí)驗(yàn)體會(huì)：通過本次實(shí)驗(yàn)，我學(xué)會(huì)了基于直接法（列主元素法）與迭代法（Jacobi迭代與GaussSeidel迭代法），和收斂速度做了詳盡比較。實(shí)驗(yàn)五、面向?qū)ο蟮乃疁?zhǔn)網(wǎng)平差程序設(shè)計(jì)與實(shí)現(xiàn)實(shí)驗(yàn)?zāi)康模涸诰仃噷?shí)現(xiàn)的基礎(chǔ)上，實(shí)現(xiàn)對水準(zhǔn)網(wǎng)平差。實(shí)驗(yàn)內(nèi)容：在深刻理解平差原理的基礎(chǔ)上，在原有矩陣實(shí)現(xiàn)的基礎(chǔ)上，考慮到水準(zhǔn)網(wǎng)數(shù)據(jù)的組織、水準(zhǔn)網(wǎng)數(shù)據(jù)編號、水準(zhǔn)網(wǎng)類的設(shè)計(jì)等問題基礎(chǔ)上，實(shí)現(xiàn)水準(zhǔn)網(wǎng)平差程序。主要代碼：class LJH_shzwNet { public: // 成員函數(shù) LJH_shzwNet(){numgaocha=0。numKnPoint=0。 // 構(gòu)造函數(shù) ~LJH_shzwNet(){}。//水準(zhǔn)網(wǎng)輸入函數(shù) void output()。////求每個(gè)點(diǎn)得近似高程 int getgczs(){return numgaocha。} //返回總點(diǎn)數(shù) int getyzds(){return numKnPoint。} //返回未知點(diǎn)數(shù) double getgcValue(int i){return thisgczhi[i1].eleValue。} //返回高程點(diǎn)編號 void setgczhi(int i, double dv){gczhi[i1].eleValue+=dv。} //修改改正數(shù) double getdv(int i){return gczhi[i1].dv。 B, LJH_CMatrixamp。 //求取系數(shù)矩陣和未知點(diǎn)高程矩陣 friend void quanzhen(LJH_CMatrixamp。 //求取權(quán)陣 friend void l_zhen( LJH_CMatrixamp。 L, LJH_shzwNet A)。 // 高差總數(shù) int numPoints。 //水準(zhǔn)網(wǎng)中已知點(diǎn)的數(shù)目 LJH_gaocha edVec[maxnum]。 // 水準(zhǔn)點(diǎn)數(shù)組 }。 double m,n,k。 cinab。 thisnumKnPoint=b。ia。 gczhi[i].dv=0。 gczhi[i].isKnown=0。ib。 cinm。 gczhi[i].isKnown=1。 cinc。 for(i=0。i++)//輸入觀測值信息 { cout請輸入第i+1段觀測段的高差值（m）、長度（km）、起始編號和終點(diǎn)編號endl。 thisedVec[i].value=k。 thisedVec[i].startPoint=a。 } } void LJH_shzwNet::output() { int i。ithisnumPoints。 for(i=0。i++)//輸出高差的信息 cout觀測段i+1平差后的值為：thisedVec[i].valueendl。 for(i=0。i++) {//如果高差觀測段的起始點(diǎn)高程值不為0，結(jié)束點(diǎn)高程為0 則：結(jié)束點(diǎn)高程=等于起始點(diǎn)高程+該測段的高差值 (0是默認(rèn)值，代表還沒參與計(jì)算) //如果高差觀測段的起始點(diǎn)高程值為0，結(jié)束點(diǎn)高程為不0 則：起始點(diǎn)高程=等于結(jié)束點(diǎn)高程該測段的高差值 if(thisgczhi[thisedVec[i].startPoint1].eleValue!=0amp。thisgczhi[thisedVec[i].endPoint1].eleValue==0) { thisgczhi[thisedVec[i].endPoint1].eleValue=thisgczhi[thisedVec[i].startPoint1].eleValue+thisedVec[i].value。amp。 } } for(j=thisnumKnPoint。j++) {//如果水準(zhǔn)點(diǎn)的高程等于0則執(zhí)行下面語句 if(thisgczhi[j].eleValue==0) { for(i=0。i++) {//如果觀測段的起始點(diǎn)等于水準(zhǔn)點(diǎn)的編號 并且結(jié)束點(diǎn)得高程不為0則該該起始點(diǎn)的高程等于結(jié)束點(diǎn)高程減去高差值（觀測值） if(thisedVec[i].startPoint==thisgczhi[i].bhamp。thisgczhi[thisedVec[i].endPoint1].eleValue!=0) { thisgczhi[i].eleValue=thisgczhi[thisedVec[i].endPoint1].eleValuethisedVec[i].value。amp。 } } } } } void xishu(LJH_CMatrix amp。 X,LJH_shzwNet A)//友元函數(shù),求誤差方程得系數(shù)矩陣 { int i,j。i()。j()。 } for(i=0。i++) {//如果起始點(diǎn)為已知點(diǎn)并且結(jié)束點(diǎn)為未知點(diǎn) 則系數(shù)設(shè)置如下 if([[i].startPoint1].isKnown==1amp。[[i].endPoint1].isKnown==0) { [i][[i].()1]=1。amp。 }//如果起始點(diǎn)和結(jié)束點(diǎn)均為未知點(diǎn) 則系數(shù)設(shè)置如下 else if([[i].startPoint1].isKnown==0amp。[[i].endPoint1].isKnown==0) { [i][[i].]=1。 } } for(i=0。i++) {//未知點(diǎn)的近似高程矩陣 [i][0]=[+i].eleValue。 P,LJH_shzwNet A)//求權(quán)陣,與距離（單位km）的長度成反比 { int i,j。 double avg。i。//水準(zhǔn)路線的總長度 avg=sum/。i。j。 } for(i=0。i++) [i][i]=avg/[i].weight。 l,LJH_CMatrix amp。 double m。i()。 } for(i=0。i++) {//（起始點(diǎn)高程+高差結(jié)束點(diǎn)高程）*1000作為l陣 m=([[i].startPoint1].eleValue+[i].[[i].endPoint1].eleValue)*1000。 } }實(shí)驗(yàn)結(jié)果:實(shí)驗(yàn)體會(huì)：做完本實(shí)驗(yàn)收獲頗多。從開始的水準(zhǔn)網(wǎng)平差中水準(zhǔn)點(diǎn)類和觀測邊類的設(shè)計(jì)開始，其中包括水準(zhǔn)點(diǎn)編號、怎樣對應(yīng)到每段的高差等，都遇到了不少問題。即使借鑒的過程，也是自己理解、學(xué)習(xí)的過程，通過本次實(shí)驗(yàn)，我切實(shí)的體會(huì)到了理論聯(lián)系實(shí)際需要做的東西很多。