【正文】 public: int _j。 operator(istream amp。：觀測(cè)邊類頭文件include fstreamusing namespace std。 double dy。 outf編號(hào)為i的點(diǎn)平差值坐標(biāo)： setiosflags(ios::fixed)setprecision(3)(*iterp).x+(*iterp).dx/100 (*iterp).y+(*iterp).dy/100endl。 i++。 CMatrix W =()*P*L。 //權(quán)計(jì)算，單位為秒平方/厘米平方 } return P。 int i=0。 j=(*iterl).startPoint1。 double fwjhc=afjhafjk。 } else if( dxjh=0 amp。 double dyjh=_dxPSet[h1].y_dxPSet[j1].y。 } else if( dxjk0 ) { afjk=afjkt+PI。 double afjk。 //角度自由項(xiàng) int i=0。 B[i][2*(knumKnPoint)]=dxjk/sjk0。 if(j=(numKnPoint) amp。 for(iterl=()。 } if(j(numKnPoint)amp。 } if(j=numKnPoint amp。 k=numKnPoint) { B[i][2*(jnumKnPoint1)]=RO*(dyjk)/sjk0。 double sjk0=dxjk*dxjk+dyjk*dyjk。 itera++,i++) { int j,k,h。 _dxPSet[k1]._isKnown=true。 iterl!=()。 _dxPSet[h1].y=_dxPSet[j1].y + ( (*iterl)._lvalue * sin(afjh) )。 for(iterl=()。 dyjh0) { afjh=2*PI+afjht。 double dyjh=_dxPSet[h1].y_dxPSet[j1].y。 } else if( dxjk0 ) { afjk=afjkt+PI。 _dxPSet[h1]._isKnown)) { double afjk。 if((_dxPSet[j1]._isKnown==true amp。 } return inf。 i++) { CJk p。 } for(i=0。ss) //重載操作符{ inf。 //權(quán)陣P CMatrix LMatrix()。 //總點(diǎn)數(shù) int numKnPoint。（列主元素法和高斯順序消元）。②邊長(zhǎng)觀測(cè)部分：設(shè)置變量j，k保存相應(yīng)端點(diǎn)編號(hào)，對(duì)兩點(diǎn)是否為起算點(diǎn)進(jìn)行討論通過(guò)下式確定誤差方程系數(shù)： 系數(shù)輸入到矩陣相應(yīng)位置。若當(dāng)前遍歷的角度觀測(cè)類對(duì)象相應(yīng)三個(gè)點(diǎn)編號(hào)的成員變量j，k，h其中j，k點(diǎn)已知（或近似坐標(biāo)已求解）且h點(diǎn)未知，則對(duì)h點(diǎn)進(jìn)行近似坐標(biāo)求解。 } }} 設(shè)計(jì)二：導(dǎo)線網(wǎng)平差程序的設(shè)計(jì)與實(shí)現(xiàn)一．設(shè)計(jì)目的立足于《數(shù)據(jù)結(jié)構(gòu)與測(cè)繪軟件開(kāi)發(fā)》這一課程的課堂教學(xué)及實(shí)驗(yàn)課程設(shè)計(jì)，為著實(shí)提高我們基于計(jì)算機(jī)輔助的方式解決工程實(shí)際問(wèn)題的動(dòng)手能力，通過(guò)本次實(shí)習(xí)，一方面我們深入了解了課堂所學(xué)知識(shí)，另一方面，通過(guò)實(shí)踐掌握測(cè)繪行業(yè)軟件設(shè)計(jì)與開(kāi)發(fā)的基本方法，深刻掌握矩陣運(yùn)算，曲線/曲面擬合的數(shù)值算法，掌握不同類型的典型測(cè)繪軟件的設(shè)計(jì)方法。 iter1!=NodeSet[cloc].()。 visited=new bool[zds]。 (sg)。 void DFS(int cloc,int zds,vectorCLevelPointamp。}：圖類頭文件include vectorpragma warning(disable: 4786)include mappragma warning(disable: 4786)include iostreampragma warning(disable: 4786)include using namespace std。 l++) { (k, l, (k, l) tem * (i, l))。 k _row * 2。 j++) { if (fabs((j, i)) fabs(mk)) // 需要考慮絕對(duì)值得情形 { mk = (j, i)。 double mk。 // 設(shè)置NXW的值 for (int i = 0。 for (int i = 0。}// 單目運(yùn)算符：取負(fù)CMatrix CMatrix::operator(){ CMatrix m(_row, _column)。}CMatrix CMatrix::operator*(const doubleamp。 i++) { for (int j = 0。 j _column。}// 自減運(yùn)算符CMatrixamp。amp。 i _row。 j++) { (i, j, _A[i][j] + (i, j))。 m){ assert(_row == () amp。 for (int i = 0。 for (int i = 0。 i _row。 for (i = 0。 i++) delete [] _A[i]。 j _column。 _column = 。 i++) _A[i] = new double[column]。 i _row。 } return outf。 operator(ostream amp。 } // 設(shè)置單位陣 bool setUnitMatrix()。 double * amp。 CMatrixamp。 m)。 m)。in,CLevelPoint amp。：水準(zhǔn)點(diǎn)類頭文件class CLevelPoint {public: CLevelPoint()。 // 權(quán)重 long _startPoint。ss) //水準(zhǔn)點(diǎn)動(dòng)態(tài)數(shù)組初始化{ vectorCLevelPoint::iterator iter。 i++,j++) _lpVec[i]._dv=(j,0)。 CMatrix W=(())*P*L。 oufendl計(jì)算結(jié)果：endl。 int i=0。 iter3!=()。 i++。 delete []g。 int k2=(*iter)._endPoint1。inumPoints。in,CElevationNet amp。 //權(quán)陣 void PCfunc()。} // 控制網(wǎng)中已知點(diǎn)的數(shù)目 void InitCD(vectorCElvDif amp。 未知點(diǎn)相應(yīng)的高程中誤差由上式求得。 2）自由項(xiàng)列向量l求解定義L矩陣，行數(shù)為高差觀測(cè)值個(gè)數(shù)，列數(shù)為1。③在進(jìn)行當(dāng)前節(jié)點(diǎn)的遍歷時(shí)，在圖的初始化中，利用STL中的容器定義的動(dòng)態(tài)數(shù)組來(lái)存放鄰接點(diǎn)和與鄰接點(diǎn)之間的高差值，因此若下一個(gè)節(jié)點(diǎn)的_iskown==0,在當(dāng)前節(jié)點(diǎn)高程（或近似高程）已知的情況下，利用高差值，求出下一個(gè)即將遍歷節(jié)點(diǎn)的近似高程，保存在水準(zhǔn)點(diǎn)對(duì)象當(dāng)中。擬采用簡(jiǎn)潔平差法求各個(gè)待定點(diǎn)的平差高程。二．設(shè)計(jì)內(nèi)容 如圖所示水準(zhǔn)網(wǎng)，A和B是已知高程的水準(zhǔn)點(diǎn)，并設(shè)這些點(diǎn)的已知高程無(wú)誤差。②圖中先關(guān)信息從文件中通過(guò)鄰接矩陣讀入完畢。設(shè)當(dāng)前遍歷第k個(gè)觀測(cè)值，起點(diǎn)，末點(diǎn)編號(hào)分別為i，j。通過(guò)法方程解得未知點(diǎn)高程改正數(shù)，與高程近似值相加得到高程平差值。 int GetnumElvDif(){return numElvDif。 //系數(shù)矩陣 void Jsgc()。 // 高程值數(shù)組 vectorCLevelPoint _lpVec。inumPoints。 for(iter=()。 for(i=0。 iter2!=()。 int i=0。 return L。 return P。 CMatrix N=(())*P*B。 int i,j。iter!=()。 virtual ~CElvDif()。s){return in。 // 高程值改正數(shù)（初始化為0） bool _isKnown。if !defined(AFX_MATRIX_H__D5670A71_87F4_49BA_9D6E_CCA34B98AD9B__INCLUDED_)define AFX_MATRIX_H__D5670A71_87F4_49BA_9D6E_CCA34B98AD9B__INCLUDED_if _MSC_VER 1000pragma onceendif // _MSC_VER 1000class CMatrix {public: CMatrix()。 // 比較運(yùn)算符 bool operator==(const CMatrixamp。 // 自加/減運(yùn)算符 CMatrixamp。 m)。 } int getRow() const { return _row。 j++) { cout_A[i][j] 。 j。 // 行 int _column。 _A = new double *[row]。 j++) _A[i][j] = 。 for (i = 0。m == this) return *this。 for (int i = 0。 } } return *this。 } } return true。 return true。 else return true。 i _row。 _column == ())。}// 自加運(yùn)算符CMatrixamp。 j _column。 for (int i = 0。 CMatrix result(_row, ())。 k++) // 這里之前出現(xiàn)了錯(cuò)誤 value += getValue(i, k) * (k, j)。 j _column。 j++) (i, j, _A[i][j])。 } } return m。 else { if (i == (j _row)) (i, j, )。 kindex = 1。 (i, k, (kindex, k))。 double tem = (k, i)。 j _row。 virtual ~CMyGraph()。 i numPoints。 j++) { if (g[i][j] != 0) { NodeSet[i].(mapint, double::value_type(j, g[i][j]))。}void CMyGraph::DFS(int cloc,vectorCLevelPointamp。 iter3++) { if((*iter1).first+1==(*iter3)._index) { break。 vectorCJk _dxLSet。 1）誤差方程系數(shù)矩陣B的建立 對(duì)角度觀測(cè)類對(duì)象的動(dòng)態(tài)數(shù)組和邊長(zhǎng)觀測(cè)類的動(dòng)態(tài)數(shù)組進(jìn)行遍歷①角度觀測(cè)部分：設(shè)置變量j，k，h保存角度觀測(cè)對(duì)象的成員變量_j,_k,_h。 法方程列立與解算 根據(jù)間接平差原理，列立法方程如下： ①B矩陣由平差計(jì)算類成員函數(shù)BMatrix（）實(shí)現(xiàn)； ②L陣由平差計(jì)算類成員函數(shù)LMatrix（）實(shí)現(xiàn)； ③P陣由邊角權(quán)運(yùn)算獲得； 矩陣運(yùn)算的功能通過(guò)設(shè)計(jì)矩陣類CMatrix實(shí)現(xiàn)。class CDxNet {public: CDxNet()。 //觀測(cè)邊長(zhǎng)public: void Jszb()。 //觀測(cè)角度數(shù)組friend istream amp。 infp。 (p)。 i++) { CJkh p。 while (1) { if(i=numPointnumKnPoint) break。amp。 if(dxjk=0 amp。 if (dyjk=0) { afjk=*PIafjkt。 } else if( dxjh0 ) { afjh=afjht+PI。 double avalue=(floor(c)+floor((cfloor(c))*100)/60+(cfloor(c)floor((cfloor(c))*100)/100)*100/36)*PI/