【正文】 ultfile){ double pvv=ca_V()。 // 殘差計算 ofstream outfile(resultfile,ios::out)。//以輸出方式打開文件，若文件不存在，創(chuàng)建文件 //輸出原始觀測數(shù)據(jù) outfileendl觀測總數(shù):gcz_zs 總點數(shù):szd_zs。 outfile 已知點數(shù):yz_szd_zsendlendl。 outfile單位權(quán)水準路線長:s0(km)endl。 outfileendl===================== 已知高程 =====================endl。//輸出原始觀測數(shù)據(jù)已知點點號、高程 for(int i=0。i=yz_szd_zs1。i++) { outfile dm[i]。 outfilesetiosflags(ios::fixed)。 outfilesetw(10)setprecision(4)szd_gc[i]endl。 } outfileendlendl===================== 高差觀測值=====================endlendl。//輸出原始觀測數(shù)據(jù)高程觀測值與路線長 outfile起始點名 終點點名 高差觀測值(m) 兩點間距離(km)endl。 for(int i=0。i=gcz_zs1。i++) { outfile dm[qsd_dh[i]]setw(9)dm[zd_dh[i]]。 outfilesetiosflags(ios::fixed)。 outfilesetw(16)setprecision(4)gcz[i]。 outfilesetiosflags(ios::fixed)。 outfilesetw(16)setprecision(4)s0/P[i]endl。 } zwc=sqrt(pvv/(gcz_zs(szd_zsyz_szd_zs)))。//計算單位權(quán)中誤差 outfileendl===================== 單位權(quán)中誤差=====================endl。//輸出單位權(quán)中誤差 outfileendlσ0=zwc(mm)endl。 outfileendl===================== 高程平差值及其精度 =====================endlendl。//輸出高程平差值及其精度 outfile點名 近似高程(m) 改正數(shù)(mm) 高程平差值(m) 中誤差(mm)endl。 for(int i=0。 iszd_zs。 i++) { outfilesetw(2)dm[i]。 double dx=dX[i]。 double qii=BTPB[ij(i,i)]。 outfilesetiosflags(ios::fixed)。//用定點格式顯示浮點數(shù) outfilesetw(12)setprecision(4)szd_gc[i]dx。 outfilesetiosflags(ios::fixed)。 outfilesetw(15)setprecision(4)dx*1000。 outfilesetiosflags(ios::fixed)。 outfilesetw(13)setprecision(4)szd_gc[i]。 outfilesetiosflags(ios::fixed)。 outfilesetw(16)setprecision(4)sqrt(qii)*zwcendl。 }//輸出觀測值平差值及其精度 outfileendlendl===================== 觀測值平差值及其精度 =====================endlendl。 outfile起點 終點 觀測高差(m) ｖ(mm) 高差平差值(m) 觀測權(quán) 中誤差(mm)endl。 for(int i=0。i=gcz_zs1。i++) { int k1=qsd_dh[i]。 int k2=zd_dh[i]。 double qii=BTPB[ij(k1,k1)]。 double qjj= BTPB[ij(k2,k2)] 。 double qij=BTPB[ij(k1,k2)]。 double ml=sqrt(qii+*qij)*zwc。 (2)。 outfiledm[k1]。 (7)。 outfiledm[k2]。 outfilesetiosflags(ios::fixed)。 outfilesetw(14)setprecision(4)gcz[i]。 outfilesetiosflags(ios::fixed)。 outfilesetw(13)setprecision(4)V[i]。 outfilesetiosflags(ios::fixed)。 outfilesetw(11)setprecision(4)gcz[i]+V[i]。 outfilesetiosflags(ios::fixed)。 outfilesetw(12)setprecision(4)P[i]。 outfilesetiosflags(ios::fixed)。 outfilesetw(11)setprecision(4)mlendl。 } ()。}///////////////////////////////////////////////////////////// // 水準網(wǎng)最小二乘平差void SZWPC::zxecpc(char *resultfile){ ca_H0()。 //近似高程計算 ca_BTPB()。 // 組成法方程 ca_dX()。 // 高程平差值計算}///////////////////////////////////////////int main(){ char *datafile,a[100]=\0,*resultfile,b[100]=\0。 cout原始數(shù)據(jù)文件存儲的路徑及文件名字:。 cina。//原始數(shù)據(jù)文件存儲地址指針 datafile=a。 cout平差后文件存儲的路徑及文件名字:。 cinb。//平差后文件存儲地址指針 resultfile=b。 coutendlendl水準網(wǎng)經(jīng)典間接平差endlendl。 SZWPC new_net。//定義新的對象 (datafile)。//輸入原始數(shù)據(jù) (resultfile)。//最小二乘平差計算 (resultfile)。//輸出平差結(jié)果 return 0。}23