【正文】 cassertpragma warning(disable: 4786)include pragma warning(disable: 4786)include iostreamusing namespace std。 m)。 m)。 m)。 m)。 CMatrixamp。 // 取負(fù)數(shù) // 乘法運(yùn)算符 CMatrix operator*(const CMatrixamp。 double * amp。 } // 設(shè)置行、列的值 void setRow(const int row) { _row = row。 } // 設(shè)置單位陣 bool setUnitMatrix()。 j_column。 operator(ostream amp。 i++) { for(int j=0。 } return outf。 int _row。 i _row。 _column = column。 i++) _A[i] = new double[column]。 j column。 _column = 。 i++) _A[i] = new double[_column]。 j _column。 m){ if (amp。 i++) delete [] _A[i]。 _A = new double *[_row]。 for (i = 0。 j++) { _A[i][j] = [i][j]。 i _row。 else _A[i][j] = 。 for (int i = 0。 j++) if (fabs(_A[i][j] (i, j)) 10e6) return false。 for (int i = 0。 j++) { if (fabs(_A[i][j] (i, j)) 10e6) continue。 m){ assert(_row == () amp。 for (int i = 0。 j++) { (i, j, _A[i][j] + (i, j))。amp。 i _row。 } } return result。amp。 i++) { for (int j = 0。}// 自減運(yùn)算符CMatrixamp。 _column == ())。 j _column。 m){ assert(_column == ())。 i++) { for (int j = 0。 k getColumn()。}CMatrix CMatrix::operator*(const doubleamp。 ++i) { for (int j = 0。}// 單目運(yùn)算符：取負(fù)CMatrix CMatrix::operator(){ CMatrix m(_row, _column)。 j _column。 for (int i = 0。 j++) { (i, j, getValue(j, i))。 // 設(shè)置NXW的值 for (int i = 0。 j++) { if (j _row) (i, j, getValue(i, j))。 double mk。 i++) { // 選取列主元 mk = (i, i)。 j++) { if (fabs((j, i)) fabs(mk)) // 需要考慮絕對(duì)值得情形 { mk = (j, i)。 k++) { aiitem = (i, k)。 k _row * 2。 k++) { if (k == i) continue。 l++) { (k, l, (k, l) tem * (i, l))。 i++) { for (int j = 0。}：圖類頭文件include vectorpragma warning(disable: 4786)include mappragma warning(disable: 4786)include iostreampragma warning(disable: 4786)include using namespace std。class CMyGraph {public: CMyGraph(int numpoints,double **g)。 void DFS(int cloc,int zds,vectorCLevelPointamp。：圖類cpp文件include pragma warning(disable: 4786)include pragma warning(disable: 4786)CMyGraph::CMyGraph(int numPoints,double **g){ for (int i = 0。 (sg)。 j numPoints。 visited=new bool[zds]。 DFS(cloc, lpvec, visited)。 iter1!=NodeSet[cloc].()。 iter3!=()。 } }} 設(shè)計(jì)二：導(dǎo)線網(wǎng)平差程序的設(shè)計(jì)與實(shí)現(xiàn)一．設(shè)計(jì)目的立足于《數(shù)據(jù)結(jié)構(gòu)與測(cè)繪軟件開發(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ì)與開發(fā)的基本方法，深刻掌握矩陣運(yùn)算，曲線/曲面擬合的數(shù)值算法，掌握不同類型的典型測(cè)繪軟件的設(shè)計(jì)方法。 ①設(shè)計(jì)觀測(cè)角度類和觀測(cè)邊長(zhǎng)類，使用stl中的容器創(chuàng)建相應(yīng)動(dòng)態(tài)數(shù)組來(lái)保存導(dǎo)線網(wǎng)觀測(cè)量的相關(guān)信息。若當(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ù)矩陣列立和自由項(xiàng)列向量的求解。②邊長(zhǎng)觀測(cè)部分：設(shè)置變量j，k保存相應(yīng)端點(diǎn)編號(hào)，對(duì)兩點(diǎn)是否為起算點(diǎn)進(jìn)行討論通過(guò)下式確定誤差方程系數(shù)： 系數(shù)輸入到矩陣相應(yīng)位置。邊角網(wǎng)的定權(quán)涉及到邊角權(quán)的運(yùn)算，一般以測(cè)角中誤差作為單位權(quán)中誤差，測(cè)距中誤差由： 確定得到邊長(zhǎng)觀測(cè)中誤差，相應(yīng)觀測(cè)邊的權(quán)按權(quán)的定義求解： 單位為mm2/″（）函數(shù)實(shí)現(xiàn)。（列主元素法和高斯順序消元）。設(shè)計(jì)流程圖： 附錄（代碼）：:平差計(jì)算類頭文件include vectorinclude include include include using namespace std。 //總點(diǎn)數(shù) int numKnPoint。 //觀測(cè)角數(shù) int numLine。 //權(quán)陣P CMatrix LMatrix()。 //觀測(cè)邊長(zhǎng)數(shù)組 vectorCJkh _dxASet。ss) //重載操作符{ inf。 i++) { CDxPoint p。 } for(i=0。 =false。 i++) { CJk p。 i。 } return inf。 itera=()。 if((_dxPSet[j1]._isKnown==true amp。!(_dxPSet[j1]._isKnown amp。 _dxPSet[h1]._isKnown)) { double afjk。 if(abs(dxjk) = abs(dyjk)) { double afjkt=atan(dyjk/dxjk)。 } else if( dxjk0 ) { afjk=afjkt+PI。 } } else { double afjkt=atan(dxjk/dyjk)。 double dyjh=_dxPSet[h1].y_dxPSet[j1].y。 dyjh=0) { afjh=afjht。 dyjh0) { afjh=2*PI+afjht。 } } } if(_dxPSet[k1]._isKnown==true) { double c=(*itera)._avalue。 for(iterl=()。 h==(*iterl).endPoint) || (h==(*iterl).startPoint amp。 _dxPSet[h1].y=_dxPSet[j1].y + ( (*iterl)._lvalue * sin(afjh) )。 double avalue=(floor(c)+floor((cfloor(c))*100)/60+(cfloor(c)floor((cfloor(c))*100)/100)*100/36)*PI/180。 iterl!=()。amp。 _dxPSet[k1]._isKnown=true。 //定義B矩陣 vectorCJkh::iterator itera。 itera++,i++) { int j,k,h。 double dxjk=_dxPSet[k1].x _dxPSet[j1].x。 double sjk0=dxjk*dxjk+dyjk*dyjk。 knumKnPoint) { B[i][2*(knumKnPoint1)]=RO*(dyjk)/sjk0。 k=numKnPoint) { B[i][2*(jnumKnPoint1)]=RO*(dyjk)/sjk0。 knumKnPoint) { B[i][2*(knumKnPoint1)]=RO*(dyjk)/sjk0。 } if(j=numKnPoint amp。 } if(jnumKnPoint amp。 } if(j(numKnPoint)amp。 B[i][2*(hnumKnPoint1)]=RO*(dyjh)/sjh0。 for(iterl=()。 k=(*iterl).endPoint1。 if(j=(numKnPoint) amp。 } if(j(numKnPoint) amp。 B[i][2*(knumKnPoint)]=dxjk/sjk0。 k=(numKnPoint)) { B[i][2*(jnumKnPoint)]=dxjk/sjk0。 //角度自由項(xiàng) int i=0。 itera++,i++) { int j,k,h。 double afjk。 if(abs(dxjk) = abs(dyjk)) { afjkt=atan(dyjk/dxjk)。 } else if( dxjk0 ) { afjk=afjkt+PI。 } } else { afjkt=atan(dxjk/dyjk)。 double dyjh=_dxPSet[h1].y_dxPSet[j1].y。 if(dxjh=0 amp。 } else if( dxjh=0 amp。 if (dyjh=0) { afjh=*PIafjht。 double fwjhc=afjhafjk。 //邊長(zhǎng)自由項(xiàng) vectorCJk::iterator iterl。 j=(*iterl).startPoint1。 double sjk0=sqrt(dxjk*dxjk+dyjk*dyjk)。 int i=0。 } vectorCJk::iterator iterl。 //權(quán)計(jì)算，單位為秒平方/厘米平方 } return P。 CMatrix L = LMatrix()。 CMatrix W =()*P*L。 for(iterp=()。 i++。 i=1。 outf編號(hào)為i的點(diǎn)平差值坐標(biāo)： setiosflags(ios::fixed)setprecision(3)(*iterp).x+(*iterp).dx/100 (*iterp).y+(*iterp).dy/100endl。public: double x。 double dy。inf,CDxPoint amp。：觀測(cè)邊類頭文件include fstreamusing namespace std。 int endPoint。 operator(istream amp。：觀測(cè)邊類cpp文件析構(gòu)函數(shù)+構(gòu)造函數(shù)。public: int _j