【正文】 double sjh0=dxjh*dxjh+dyjh*dyjh。 if(j=numKnPoint amp。amp。 knumKnPoint) { B[i][2*(knumKnPoint1)]=RO*(dyjk)/sjk0。 B[i][2*(knumKnPoint1)+1]=RO*(dxjk)/sjk0。 } if(jnumKnPoint amp。amp。 k=numKnPoint) { B[i][2*(jnumKnPoint1)]=RO*(dyjk)/sjk0。 B[i][2*(jnumKnPoint1)+1]=RO*(dxjk)/sjk0。 } if(jnumKnPoint amp。amp。 knumKnPoint) { B[i][2*(knumKnPoint1)]=RO*(dyjk)/sjk0。 B[i][2*(knumKnPoint1)+1]=RO*(dxjk)/sjk0。 B[i][2*(jnumKnPoint1)]=RO*(dyjk)/sjk0。 B[i][2*(jnumKnPoint1)+1]=RO*(dxjk)/sjk0。 } if(j=numKnPoint amp。amp。 hnumKnPoint) { B[i][2*(hnumKnPoint1)]=RO*(dyjh)/sjh0。 B[i][2*(hnumKnPoint1)+1]=RO*(dxjh)/sjh0。 } if(jnumKnPoint amp。amp。 h=numKnPoint) { B[i][(jnumKnPoint1)*2]=B[i][(jnumKnPoint1)*2]+RO*dyjh/sjh0。 B[i][(jnumKnPoint1)*2+1]=B[i][(jnumKnPoint1)*2+1]RO*dxjh/sjh0。 } if(j(numKnPoint)amp。amp。 hnumKnPoint) { B[i][2*(jnumKnPoint1)]=B[i][2*(jnumKnPoint1)]+RO*dyjh/sjh0。 B[i][2*(jnumKnPoint1)+1]=B[i][2*(jnumKnPoint1)+1]RO*dxjh/sjh0。 B[i][2*(hnumKnPoint1)]=RO*(dyjh)/sjh0。 B[i][2*(hnumKnPoint1)+1]=RO*(dxjh)/sjh0。 } } //角度系數(shù)計(jì)算完畢 vectorCJk::iterator iterl。 //邊長(zhǎng)系數(shù)計(jì)算開(kāi)始 i=numAngel。 for(iterl=()。 iterl!=()。 iterl++,i++) { int j,k。 j=(*iterl).startPoint1。 k=(*iterl).endPoint1。 double dxjk=_dxPSet[k1].x _dxPSet[j1].x。 double dyjk=_dxPSet[k1].y _dxPSet[j1].y。 double sjk0=sqrt(dxjk*dxjk+dyjk*dyjk)。 if(j=(numKnPoint) amp。amp。 k(numKnPoint)) { B[i][2*(knumKnPoint)]=dxjk/sjk0。 B[i][2*(knumKnPoint)+1] =dyjk/sjk0。 } if(j(numKnPoint) amp。amp。 k(numKnPoint)) { B[i][2*(jnumKnPoint)]=dxjk/sjk0。 B[i][2*(jnumKnPoint)+1] =dyjk/sjk0。 B[i][2*(knumKnPoint)]=dxjk/sjk0。 B[i][2*(knumKnPoint)+1] =dyjk/sjk0。 } if(j(numKnPoint) amp。amp。 k=(numKnPoint)) { B[i][2*(jnumKnPoint)]=dxjk/sjk0。 B[i][2*(jnumKnPoint)+1] =dyjk/sjk0。 } } //邊長(zhǎng)系數(shù)計(jì)算結(jié)束 return B。}CMatrix CDxNet::LMatrix(){ CMatrix L(numAngel+numLine,1)。 //角度自由項(xiàng) int i=0。 vectorCJkh::iterator itera。 for(itera=()。 itera != ()。 itera++,i++) { int j,k,h。 j=(*itera)._j。 k=(*itera)._k。 h=(*itera)._h。 double afjk。 double afjkt。 double dxjk=_dxPSet[k1].x_dxPSet[j1].x。 double dyjk=_dxPSet[k1].y_dxPSet[j1].y。 if(abs(dxjk) = abs(dyjk)) { afjkt=atan(dyjk/dxjk)。 if(dxjk=0 amp。amp。 dyjk=0) { afjk=afjkt。 } else if( dxjk0 ) { afjk=afjkt+PI。 } else if( dxjk=0 amp。amp。 dyjk0) { afjk=2*PI+afjkt。 } } else { afjkt=atan(dxjk/dyjk)。 if (dyjk=0) { afjk=*PIafjkt。 } else { afjk=*PIafjkt。 } } double dxjh=_dxPSet[h1].x_dxPSet[j1].x。 double dyjh=_dxPSet[h1].y_dxPSet[j1].y。 double afjh。 double afjht。 if(abs(dxjh) = abs(dyjh)) { afjht=atan(dyjh/dxjh)。 if(dxjh=0 amp。amp。 dyjh=0) { afjh=afjht。 } else if( dxjh0 ) { afjh=afjht+PI。 } else if( dxjh=0 amp。amp。 dyjh0) { afjh=2*PI+afjht。 } } else { afjht=atan(dxjh/dyjh)。 if (dyjh=0) { afjh=*PIafjht。 } else { afjh=*PIafjht。 } } double c=(*itera)._avalue。 double avalue=(floor(c)+floor((cfloor(c))*100)/60+(cfloor(c)floor((cfloor(c))*100)/100)*100/36)*PI/180。 double fwjhc=afjhafjk。 if (fwjhc0) {fwjhc=2*PI+fwjhc。} L[i][0]=(avaluefwjhc)*206265。 //角度自由項(xiàng)計(jì)算完成，單位為秒 } i=numAngel。 //邊長(zhǎng)自由項(xiàng) vectorCJk::iterator iterl。 for(iterl=()。 iterl != ()。 iterl++,i++) { int j,k。 j=(*iterl).startPoint1。 k=(*iterl).endPoint1。 double dxjk=_dxPSet[k1].x_dxPSet[j1].x。 double dyjk=_dxPSet[k1].y_dxPSet[j1].y。 double sjk0=sqrt(dxjk*dxjk+dyjk*dyjk)。 L[i][0]=((*iterl)._lvaluesjk0)*100 。 //邊長(zhǎng)自由項(xiàng)計(jì)算結(jié)束，單位為cm } return L。}CMatrix CDxNet::PMatrix(){ CMatrix P(numAngel+numLine,numAngel+numLine)。 int i=0。 for(i=0。 inumAngel。 i++) { P[i][i]=1。 } vectorCJk::iterator iterl。 for(iterl=()。 iterl!=()。 iterl++) { P[i][i]=(sigmabeta*sigmabeta/(*iterl)._lvalue)*100。 //權(quán)計(jì)算，單位為秒平方/厘米平方 } return P。}void CDxNet::PCFunc(){ ofstream outf(C:\\Users\\liyu\\Desktop\\,ios::out | ios::app)。 Jszb()。 CMatrix B = BMatrix()。 CMatrix L = LMatrix()。 CMatrix P = PMatrix()。 CMatrix N = ()*P*B。 CMatrix QXX=()。 CMatrix W =()*P*L。 CMatrix X = ()*W。 int i=0。 vectorCDxPoint::iterator iterp。 for(iterp=()。 iterp!=()。 iterp++) { if((*iterp)._valuenumKnPoint) { (*iterp).dx=X[i][0]。 (*iterp).dy=X[i+1][0]。 i++。 } } outfendlendl誤差方程系數(shù)矩陣B：endlB。 outf權(quán)陣：endlPendl。 outfQXX:endlQXXendl。 i=1。 for(iterp=()。 iterp!=()。 iterp++,i++) { outf編號(hào)為i的點(diǎn)近似坐標(biāo)： setiosflags(ios::fixed)setprecision(3)(*iterp).x (*iterp).y 。 outf編號(hào)為i的點(diǎn)平差值坐標(biāo)： setiosflags(ios::fixed)setprecision(3)(*iterp).x+(*iterp).dx/100 (*iterp).y+(*iterp).dy/100endl。 }}:導(dǎo)線點(diǎn)類頭文件include fstreamusing namespace std。class CDxPoint {public: CDxPoint()。 virtual ~CDxPoint()。public: double x。 double y。 int _value。 double dx。 double dy。 bool _isKnown。friend istream amp。 operator(istream amp。inf,CDxPoint amp。ss){ return inf。}}。:導(dǎo)線點(diǎn)類cpp文件析構(gòu)函數(shù)+構(gòu)造函數(shù)。：觀測(cè)邊類頭文件include fstreamusing namespace std。class CJk {public: CJk()。 virtual ~CJk()。public: int startPoint。 int endPoint。 double _lvalue。 double weight。friend istream amp。 operator(istream amp。inf,CJk amp。ss){ return inf。}}。：觀測(cè)邊類cpp文件析構(gòu)函數(shù)+構(gòu)造函數(shù)。:觀測(cè)角度類頭文件include fstreamusing namespace std。class CJkh {public: CJkh()。 virtual ~CJkh()。public: int _j。 //角度定點(diǎn)號(hào) int _k。 int _h。 //從jk沿所測(cè)角度旋轉(zhuǎn)至jh為順時(shí)針