誤差理論與測(cè)量平差課程設(shè)計(jì)報(bào)告(參考版)
2025-06-26 07:53本頁(yè)面
【正文】 }23。//最小二乘平差計(jì)算 (resultfile)。//定義新的對(duì)象 (datafile)。 coutendlendl水準(zhǔn)網(wǎng)經(jīng)典間接平差endlendl。 cinb。//原始數(shù)據(jù)文件存儲(chǔ)地址指針 datafile=a。 cout原始數(shù)據(jù)文件存儲(chǔ)的路徑及文件名字:。 // 組成法方程 ca_dX()。}///////////////////////////////////////////////////////////// // 水準(zhǔn)網(wǎng)最小二乘平差void SZWPC::zxecpc(char *resultfile){ ca_H0()。 outfilesetw(11)setprecision(4)mlendl。 outfilesetw(12)setprecision(4)P[i]。 outfilesetw(11)setprecision(4)gcz[i]+V[i]。 outfilesetw(13)setprecision(4)V[i]。 outfilesetw(14)setprecision(4)gcz[i]。 outfiledm[k2]。 outfiledm[k1]。 double ml=sqrt(qii+*qij)*zwc。 double qjj= BTPB[ij(k2,k2)] 。 int k2=zd_dh[i]。i=gcz_zs1。 outfile起點(diǎn) 終點(diǎn) 觀測(cè)高差(m) ｖ(mm) 高差平差值(m) 觀測(cè)權(quán) 中誤差(mm)endl。 outfilesetw(16)setprecision(4)sqrt(qii)*zwcendl。 outfilesetw(13)setprecision(4)szd_gc[i]。 outfilesetw(15)setprecision(4)dx*1000。//用定點(diǎn)格式顯示浮點(diǎn)數(shù) outfilesetw(12)setprecision(4)szd_gc[i]dx。 double qii=BTPB[ij(i,i)]。 i++) { outfilesetw(2)dm[i]。 for(int i=0。 outfileendl===================== 高程平差值及其精度 =====================endlendl。//計(jì)算單位權(quán)中誤差 outfileendl===================== 單位權(quán)中誤差=====================endl。 outfilesetw(16)setprecision(4)s0/P[i]endl。 outfilesetw(16)setprecision(4)gcz[i]。i++) { outfile dm[qsd_dh[i]]setw(9)dm[zd_dh[i]]。 for(int i=0。 } outfileendlendl===================== 高差觀測(cè)值=====================endlendl。 outfilesetiosflags(ios::fixed)。i=yz_szd_zs1。 outfileendl===================== 已知高程 =====================endl。 outfile 已知點(diǎn)數(shù):yz_szd_zsendlendl。 // 殘差計(jì)算 ofstream outfile(resultfile,ios::out)。 } return(pvv)。 V[i]=(szd_gc[k2]szd_gc[k1]gcz[i])*1000。i++) { int k1=qsd_dh[i]。 for(int i=0。 szd_gc[i]+=xi。 j++) { xi+=BTPB[ij(i,j)]*BTPL[j]。 for(int j=0。 iszd_zs。//矩陣為奇異矩陣，無法求逆 exit(0)。i++) BTPB[ij(i,i)]=。 }}/////////////////////////// 高程平差值計(jì)算void SZWPC::ca_dX(){ for(int i=0。//獲得法方程系數(shù)矩陣 BTPB[ij(j,j)]+=Pk。//獲得法方程自由項(xiàng) BTPL[j]+=Pk*lk。//獲取權(quán)值 double lk=gcz[k](szd_gc[j]szd_gc[i])。//獲取點(diǎn)號(hào) int j=zd_dh[k]。 kgcz_zs。 i++) BTPL[i]=。//賦初值 for(int i=0。 iszd_zs*(szd_zs+1)/2。//所有的近似高程計(jì)算完成，退出 }}//////////////////////////////// 組成法方程void SZWPC::ca_BTPB(){ int t=szd_zs。//計(jì)算近似高程 k++。amp。//計(jì)算近似高程 k++。amp。 //高差起點(diǎn)號(hào) int k2=zd_dh[i]。igcz_zs。j++) { int k=0。i++)szd_gc[i]=。 //dm數(shù)組已經(jīng)存滿，且沒有待查點(diǎn)名}//////////高程近似值計(jì)算void SZWPC::ca_H0(){ for(int i=yz_szd_zs。//待查點(diǎn)名是一個(gè)新的點(diǎn)名，將新點(diǎn)名的地址放到dm數(shù)組中 return i。//判斷點(diǎn)名長(zhǎng)度 dm[i]=new char[len+1]。i++) { if(dm[i]!=NULL) { if(strcmp(name,dm[i])==0)return i。}////////// 點(diǎn)名轉(zhuǎn)點(diǎn)號(hào)，返回點(diǎn)名對(duì)應(yīng)的點(diǎn)號(hào)int SZWPC::dm_dh(char *name){ for(int i=0。 //讀取高差值與路線長(zhǎng)度 P[i]=s0/P[i]。//讀取高程終點(diǎn) zd_dh[i]=dm_dh(buffer)。 //讀取高程起點(diǎn)名 qsd_dh[i]=dm_dh(buffer)。igcz_zs。 infileszd_gc[i]。i++)// 讀取已知高程數(shù)據(jù) { infilebuffer。 //臨時(shí)數(shù)組，保存從文件中讀到的點(diǎn)名 for(int i=0。i++)
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1