【正文】 長觀測類的動(dòng)態(tài)數(shù)組進(jìn)行遍歷①角度觀測部分：設(shè)置變量j，k，h保存角度觀測對(duì)象的成員變量_j,_k,_h。 ④各導(dǎo)線點(diǎn)近似坐標(biāo)都已求解完畢，break跳出循環(huán)。 ③起算邊jk坐標(biāo)方位角由坐標(biāo)反算求得，與角度觀測值相加得到j(luò)，h邊近似坐標(biāo)方位角，j，h邊邊長通過點(diǎn)的編號(hào)在邊長觀測數(shù)組_dxLSet中索引得到。 //觀測角度數(shù)組 ②使用一個(gè)while循環(huán)對(duì)角度觀測數(shù)組_dxASet進(jìn)行多輪的遍歷。 vectorCJk _dxLSet。3． 關(guān)鍵問題描述 未知點(diǎn)近似坐標(biāo)計(jì)算 根據(jù)導(dǎo)線網(wǎng)布設(shè)原理，未知點(diǎn)近似坐標(biāo)通過兩起算點(diǎn)構(gòu)成起算邊方位角，以及觀測邊長和觀測角從而解得坐標(biāo)增量△x和△y，通過起算點(diǎn)坐標(biāo)求得未知點(diǎn)近似坐標(biāo)。使我們初步具備編寫測繪軟件常用算法的能力以及開發(fā)中小規(guī)模測繪專業(yè)軟件的能力。 } DFS((*iter1).first,lpvec,visited)。 iter3++) { if((*iter1).first+1==(*iter3)._index) { break。 //+1，從lpvec中找出該點(diǎn) for(iter3=()。 iter1++) { if(visited[(*iter1).first]==false) { visited[(*iter1).first]=true。 for(iter1=NodeSet[cloc].()。}void CMyGraph::DFS(int cloc,vectorCLevelPointamp。 i++) visited[i]=false。 for(int i=0。 lpvec){ bool* visited。 j++) { if (g[i][j] != 0) { NodeSet[i].(mapint, double::value_type(j, g[i][j]))。 i++) //鄰接點(diǎn)和高差值的輸入 { for (int j = 0。 } for (i = 0。 = i。 i numPoints。}。 lpvec)。 lpvec,bool* visited)。 virtual ~CMyGraph()。 // 與當(dāng)前結(jié)點(diǎn)相鄰的結(jié)點(diǎn)集合，及其與相鄰結(jié)點(diǎn)之間路徑高差值}。struct SGNode { int key。 } } return result。 j _row。 i _row。 } } } CMatrix result(_row, _row)。 l _row * 2。 double tem = (k, i)。 k _row。 k++) (i, k, (i, k) / mk)。 } } // 交換行完成 // 消去 for (int k = 0。 (i, k, (kindex, k))。 k _row * 2。 kindex = j。 j _row。 kindex = 1。 i _row。 // mk為列主元的值，kindex用于存儲(chǔ)列主元所在的位置（行） double aiitem。 } } } // 用高斯若當(dāng)消去法求解平面方程的系數(shù)，只要對(duì)NTW進(jìn)行消去即可求得系數(shù) long kindex。 else { if (i == (j _row)) (i, j, )。 j _row * 2。 i _row。 CMatrix NXW(_row, _row * 2)。 } } return m。 j _row。 i _column。}// 轉(zhuǎn)置CMatrix CMatrix::transpose(){ CMatrix m(_column, _row)。 j++) (i, j, _A[i][j])。 i++) { for (int j = 0。 for (int i = 0。 } } return result。 j _column。 i _row。 s){ CMatrix result(_row, _column)。 } } return result。 k++) // 這里之前出現(xiàn)了錯(cuò)誤 value += getValue(i, k) * (k, j)。 for (int k = 0。 j ()。 i _row。 CMatrix result(_row, ())。}// 乘法運(yùn)算CMatrix CMatrix::operator*(const CMatrixamp。 j++) _A[i][j] = (i, j)。 i++) { for (int j = 0。 for (int i = 0。amp。 CMatrix::operator=(const CMatrixamp。 } } return *this。 j _column。 i _row。 _column == ())。 m){ assert(_row == () amp。}// 自加運(yùn)算符CMatrixamp。 j++) { (i, j, _A[i][j] (i, j))。 i++) { for (int j = 0。 for (int i = 0。 _column == ())。 m){ assert(_row == () amp。 } } return result。 j _column。 i _row。 CMatrix result(_row, _column)。amp。}// 加法運(yùn)算符CMatrix CMatrix::operator+(const CMatrixamp。 else return true。 j _column。 i _row。 m){ if ((_row != ()) || (_column != ())) return true。 return true。 j _column。 i _row。 m){ if ((_row != ()) || (_column != ())) return false。 } } return true。 ++j) { if (i == j) _A[i][j] = 。 ++i) { for (int j = 0。 for (int i = 0。 } } return *this。 j _column。 i _row。 i++) _A[i] = new double[_column]。 for (int i = 0。 _column = 。 delete [] _A。 i _row。m == this) return *this。 CMatrix::operator=(const CMatrixamp。 j++) { _A[i][j] = [i][j]。 i++) { for (int j = 0。 for (i = 0。 i _row。 _A = new double *[_row]。 m){ _row = 。 j++) _A[i][j] = 。 i++) for (int j = 0。 // Initialize for (i = 0。 i row。 _A = new double *[row]。 }}// 構(gòu)造函數(shù)一CMatrix::CMatrix(int row, int column){ _row = row。 i++) delete[] _A[i]。：矩陣類cpp文件include include CMatrix::CMatrix():_A(NULL), _row(1), _column(1){}CMatrix::~CMatrix(){ if (NULL != _A) { for (int i = 0。 // 行 int _column。 // 列主元素法private: // 成員變量 double** _A。 }public: // 矩陣轉(zhuǎn)置 CMatrix transpose()。 outfendl。 j。 i。outf,CMatrix amp。 } } friend ostream amp。 j++) { cout_A[i][j] 。 i++) { for(int j=0。 //輸出到屏幕 void Output() { for(int i=0。 } int getColumn() const { return _column。 } int getRow() const { return _row。 } double getValue(int row, int column) const { return _A[row][column]。 operator[](int i){return _A[i]。 s)。 m)。 // 單目運(yùn)算符 CMatrix operator()。 operator=(const CMatrixamp。 m)。 // 自加/減運(yùn)算符 CMatrixamp。 CMatrix operator(const CMatrixamp。 // 加/減運(yùn)算符 CMatrix operator+(const CMatrixamp。 bool operator!=(const CMatrixamp。 // 比較運(yùn)算符 bool operator==(const CMatrixamp。 operator=(const CMatrixamp。public: // 默認(rèn)析構(gòu)函數(shù) virtual ~CMatrix()。 // 復(fù)制構(gòu)造函數(shù) CMatrix(const CMatrixamp。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()。：構(gòu)造函數(shù)+析構(gòu)函數(shù)。s){return in。 operator(istream amp。 // 高程值改正數(shù)（初始化為0） bool _isKnown。 // 水準(zhǔn)點(diǎn)編號(hào) double _eleValue。 virtual ~CLevelPoint()。：構(gòu)造函數(shù)+析構(gòu)函數(shù)。s){return in。 operator(istream amp。 // 起始點(diǎn)編號(hào) long _endPoint。 // 觀測值 double _weight。 virtual ~CElvDif()。iter++) { (*iter)。 for(iter=()。 }}void CElevationNet::InitCP(vectorCLevelPoint amp。iter!=()。ss) //高差觀測動(dòng)態(tài)數(shù)組初始化{ vectorCElvDif::iterator iter。 ()。 inumPoints。 int i,j。 ouf改正數(shù)V:endlVendl。 CMatrix X=(())*W。 ouf未知點(diǎn)高程協(xié)因數(shù)陣QXX:endlQXXendl。 CMatrix N=(())*P*B。 CMatrix L=LMatrix()。 CMatrix B=BMatrix()。 Jsgc()。 return P。 iter4++,i++) { (i,i,(iter4_weight))。 for(iter4=()。 vectorCElvDif::iterator iter4。 return L。 i++。 iter3++) { l=(*iter3)._value_lpVec[(*iter3)._endPoint1]._eleValue+_lpVec[(*iter3)._startPoint1]._eleValue。 for(iter3=()。 int i=0。}CMatrix CElevationNet::LMatrix() //求解自由項(xiàng)L陣{ CMatrix L(numElvDif,1)。 } ()。 if((*iter2)._endPointnumKnPoint) (i,(*iter2)._endPointnumKnPoint1,1)。 iter2!=()。 int