【正文】 /返回已知點(diǎn)數(shù) int getwzds(){return numPointsnumKnPoint。} //返回未知點(diǎn)數(shù) double getgcValue(int i){return thisgczhi[i1].eleValue。} //獲得高程值 long getgcbh(int i){return gczhi[i1].bh。} //返回高程點(diǎn)編號(hào) void setgczhi(int i, double dv){gczhi[i1].eleValue+=dv。} //改正高程值 void setdv(int i,double dv){thisgczhi[i1].dv=dv。} //修改改正數(shù) double getdv(int i){return gczhi[i1].dv。} //返回改正數(shù) friend void xishu(LJH_CMatrixamp。 B, LJH_CMatrixamp。 X,LJH_shzwNet A)。 //求取系數(shù)矩陣和未知點(diǎn)高程矩陣 friend void quanzhen(LJH_CMatrixamp。 Q, LJH_shzwNet A)。 //求取權(quán)陣 friend void l_zhen( LJH_CMatrixamp。 l ,LJH_CMatrixamp。 L, LJH_shzwNet A)。 //求取觀測(cè)值L矩陣和l陣 public: // 成員變量 int numgaocha。 // 高差總數(shù) int numPoints。 // 水準(zhǔn)網(wǎng)中點(diǎn)的總數(shù)目 int numKnPoint。 //水準(zhǔn)網(wǎng)中已知點(diǎn)的數(shù)目 LJH_gaocha edVec[maxnum]。 // 觀測(cè)值數(shù)組 LJH_shzpoint gczhi[maxnum]。 // 水準(zhǔn)點(diǎn)數(shù)組 }。 void LJH_shzwNet::input() { int a,b,c,i。 double m,n,k。 cout請(qǐng)輸入水準(zhǔn)網(wǎng)中水準(zhǔn)點(diǎn)總數(shù)和已知點(diǎn)數(shù)：endl。 cinab。 thisnumPoints=a。 thisnumKnPoint=b。 for(i=0。ia。i++)//對(duì)水準(zhǔn)網(wǎng)中的點(diǎn)進(jìn)行編號(hào)，并對(duì)水準(zhǔn)點(diǎn)進(jìn)行初始化 { gczhi[i].bh=i+1。 gczhi[i].dv=0。 gczhi[i].eleValue=0。 gczhi[i].isKnown=0。 } for(i=0。ib。i++)//為水準(zhǔn)網(wǎng)中的已知點(diǎn)輸入高程 { cout請(qǐng)輸入第i+1個(gè)已知點(diǎn)的高程值：endl。 cinm。 gczhi[i].eleValue=m。 gczhi[i].isKnown=1。 } cout請(qǐng)輸入觀測(cè)值的個(gè)數(shù)：endl。 cinc。 thisnumgaocha=c。 for(i=0。ithisnumgaocha。i++)//輸入觀測(cè)值信息 { cout請(qǐng)輸入第i+1段觀測(cè)段的高差值（m）、長(zhǎng)度（km）、起始編號(hào)和終點(diǎn)編號(hào)endl。 cinknab。 thisedVec[i].value=k。 thisedVec[i].weight=n。 thisedVec[i].startPoint=a。 thisedVec[i].endPoint=b。 } } void LJH_shzwNet::output() { int i。 for(i=0。ithisnumPoints。i++)//輸出水準(zhǔn)點(diǎn)的信息 cout編號(hào)為thisgczhi[i].bh的高程值為：thisgczhi[i].eleValueendl。 for(i=0。ithisnumgaocha。i++)//輸出高差的信息 cout觀測(cè)段i+1平差后的值為：thisedVec[i].valueendl。 } void LJH_shzwNet::jsgc()//求每個(gè)點(diǎn)得近似高程 { int i,j。 for(i=0。ithisnumgaocha。i++) {//如果高差觀測(cè)段的起始點(diǎn)高程值不為0，結(jié)束點(diǎn)高程為0 則：結(jié)束點(diǎn)高程=等于起始點(diǎn)高程+該測(cè)段的高差值 (0是默認(rèn)值，代表還沒(méi)參與計(jì)算) //如果高差觀測(cè)段的起始點(diǎn)高程值為0，結(jié)束點(diǎn)高程為不0 則：起始點(diǎn)高程=等于結(jié)束點(diǎn)高程該測(cè)段的高差值 if(thisgczhi[thisedVec[i].startPoint1].eleValue!=0amp。amp。thisgczhi[thisedVec[i].endPoint1].eleValue==0) { thisgczhi[thisedVec[i].endPoint1].eleValue=thisgczhi[thisedVec[i].startPoint1].eleValue+thisedVec[i].value。 }else if(thisgczhi[thisedVec[i].startPoint1].eleValue==0amp。amp。thisgczhi[thisedVec[i].endPoint1].eleValue!=0) { thisgczhi[thisedVec[i].startPoint1].eleValue=thisgczhi[thisedVec[i].endPoint1].eleValuethisedVec[i].value。 } } for(j=thisnumKnPoint。jthisnumPoints。j++) {//如果水準(zhǔn)點(diǎn)的高程等于0則執(zhí)行下面語(yǔ)句 if(thisgczhi[j].eleValue==0) { for(i=0。ithisnumgaocha。i++) {//如果觀測(cè)段的起始點(diǎn)等于水準(zhǔn)點(diǎn)的編號(hào) 并且結(jié)束點(diǎn)得高程不為0則該該起始點(diǎn)的高程等于結(jié)束點(diǎn)高程減去高差值（觀測(cè)值） if(thisedVec[i].startPoint==thisgczhi[i].bhamp。amp。thisgczhi[thisedVec[i].endPoint1].eleValue!=0) { thisgczhi[i].eleValue=thisgczhi[thisedVec[i].endPoint1].eleValuethisedVec[i].value。 }//如果觀測(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。amp。thisgczhi[thisedVec[i].startPoint1].eleValue!=0) { thisgczhi[i].eleValue=thisgczhi[thisedVec[i].startPoint1].eleValue+thisedVec[i].value。 } } } } } void xishu(LJH_CMatrix amp。 B, LJH_CMatrix amp。 X,LJH_shzwNet A)//友元函數(shù),求誤差方程得系數(shù)矩陣 { int i,j。 for(i=0。i()。i++) for(j=0。j()。j++) { [i][j]=0。 } for(i=0。i()。i++) {//如果起始點(diǎn)為已知點(diǎn)并且結(jié)束點(diǎn)為未知點(diǎn) 則系數(shù)設(shè)置如下 if([[i].startPoint1].isKnown==1amp。amp。[[i].endPoint1].isKnown==0) { [i][[i].()1]=1。 }//如果起始點(diǎn)為未知點(diǎn)并且結(jié)束點(diǎn)為已知點(diǎn) 則系數(shù)設(shè)置如下 else if([[i].startPoint1].isKnown==0amp。amp。[[i].endPoint1].isKnown==1) { [i][[i].()1]=1。 }//如果起始點(diǎn)和結(jié)束點(diǎn)均為未知點(diǎn) 則系數(shù)設(shè)置如下 else if([[i].startPoint1].isKnown==0amp。amp。[[i].endPoint1].isKnown==0) { [i][[i].]=1。 [i][[i].]=1。 } } for(i=0。i。i++) {//未知點(diǎn)的近似高程矩陣 [i][0]=[+i].eleValue。 } } void quanzhen(LJH_CMatrix amp。 P,LJH_shzwNet A)//求權(quán)陣,與距離（單位km）的長(zhǎng)度成反比 { int i,j。 double sum=0。 double avg。 for(i=0。i。i++) sum+=[i].weight。//水準(zhǔn)路線(xiàn)的總長(zhǎng)度 avg=sum/。//水準(zhǔn)路線(xiàn)平均長(zhǎng)度 for(i=0。i。i++) { for(j=0。j。j++) [i][j]=0。 } for(i=0。i。i++) [i][i]=avg/[i].weight。 } void l_zhen(LJH_CMatrix amp。 l,LJH_CMatrix amp。L,LJH_shzwNet A)//求取觀測(cè)值L矩陣和l陣 { int i。 double m。 for(i=0。i()。i++) { [i][0]=[i].value。 } for(i=0。i。i++) {//（起始點(diǎn)高程+高差結(jié)束點(diǎn)高程）*1000作為l陣 m=([[i].startPoint1].eleValue+[i].[[i].endPoint1].eleValue)*1000。 [i][0]=m。 } }實(shí)驗(yàn)結(jié)果:實(shí)驗(yàn)體會(huì)：做完本實(shí)驗(yàn)收獲頗多。首先通過(guò)本實(shí)驗(yàn)，我發(fā)現(xiàn)雖然自己以前平差課學(xué)的很好，但真正對(duì)應(yīng)到編程實(shí)現(xiàn)時(shí)還會(huì)存在很多問(wèn)題。從開(kāi)始的水準(zhǔn)網(wǎng)平差中水準(zhǔn)點(diǎn)類(lèi)和觀測(cè)邊類(lèi)的設(shè)計(jì)開(kāi)始，其中包括水準(zhǔn)點(diǎn)編號(hào)、怎樣對(duì)應(yīng)到每段的高差等，都遇到了不少問(wèn)題。通過(guò)查資料，以及參考往年水準(zhǔn)網(wǎng)平差設(shè)計(jì)的程序，漸漸的理清了思路。即使借鑒的過(guò)程，也是自己理解、學(xué)習(xí)的過(guò)程，通過(guò)本次實(shí)驗(yàn)，我切實(shí)的體會(huì)到了理論聯(lián)系實(shí)際需要做的東西很多。真正完成本次實(shí)驗(yàn)后，我也收獲了很多。