誤差理論與測量平差課程設(shè)計--水準(zhǔn)網(wǎng)嚴(yán)密平差及精度評定-文庫吧
2025-05-18 07:21 本頁面
【正文】 矩陣 B 和自由項 l，并由它們組成法方程 NBBxW=0（其中 NBB=BTPB， W=BTPl），法方程的個數(shù)等于所選參數(shù)的個數(shù)。 （ 4）解算法方程，求出參數(shù)改正值 x 并計算參數(shù)的平差值 X=X0+x。 （ 5）由誤差方程計算 V，并求出觀測量的平差值。為了檢查平差計算的正確性，將所求的值代入條件方程，看其是否滿足方程。 （ 6）精度評定，計算單位權(quán)中誤差，按照題設(shè)要求列出權(quán)函數(shù)式，再根據(jù)平差參數(shù)的協(xié)方差陣求出協(xié)因數(shù)，最后求出某段高差中誤差，某些點的高程中誤差。 二、程序設(shè)計思想 考慮到在解題過程中一些計算的復(fù)雜性，我們需借助一些技術(shù)將計算簡單化，快捷化，因此在課程設(shè)計過程中，我們把一些 C語言程序設(shè)計引入其中；通過一些簡單、明了的程序及子函數(shù)調(diào)用，我們就可以很方便快捷的求出用 筆算比較繁瑣、費時的矩陣乘積、矩陣的逆（如 BTPB、 BTPl）等運算。 第四部分 程序流程圖 根據(jù)題目列出條件方程并寫成誤差方程的形式 V=Bxl ↓ 5 確定權(quán)陣，根據(jù)誤差方程得到矩陣 B、 l進而寫出 BT ↓ 運用 C 程序語言求出 BTP，進一步得到 NBB=BTPB、 W=BTPl 并求出 NBB1 ↓ 用 C程序求出參數(shù)的改正數(shù) x=NBB1W ↓ 根據(jù) C程序語言求 Bx，進而由 V=Bxl寫出各觀測值的改正數(shù) ↓ 根據(jù) L=L+V 求出各觀測值的平差值 ↓ 檢驗所求各值是否正確，若無誤則往下進行，反之檢查 各步驟查出錯誤并改正 ↓ 由程序計算 VTP進而求出 VTPV，求單位權(quán)中誤差，再根據(jù)權(quán)函數(shù)式、協(xié)因數(shù)傳播定律評定各觀測值及所求高程的精度 第五部分 程序及說明 一、矩陣相乘計算函數(shù) include “” void Matrix(a,b,m,n,k,c) int m,n,k。 double a[],b[],c[]。 { int i,j,l,u。 for(i=0。i=m1。i++) for(j=0。j=k1。j++) { 6 u=i*k+j。c[u]=。 for(l=0。l=n1。l++) c[u]=c[u]+a[i*n+l]*b[l*k+j]。 } return。 } BTP main() { int i,j。 static double a[3][7]=BT； static double c[3][7],b[7][7]=P； Matrixmul(a,b,3,7,7,c)。 printf(“\n”)。 for(i=0。i=2。i++) { for(j=0。j=6。j++) printf(“%\t”,c[i][j]。 printf(“\n”)。 } printf(“\n”)。 return0。 } BTPB，即 NBB main() { int i,j。 static double a[3][7]=BTP； static double c[3][3],b[7][3]=B； Matrixmul(a,b,3,7,3,c)。 printf(“\n”)。 for(i=0。i=2。i++) { for(j=0。j=2。j++) 7 printf(“%\t”,c[i][j]。 printf(“\n”)。 } printf(“\n”)。 return0。 } BTPl，即 W main() { int i,j。 static double a[3][7]=BTP； static double c[3][1],b[7][1]=l； Matrixmul(a,b,3,7,1,c)。 printf(“\n”)。 for
畢業(yè)設(shè)計相關(guān)推薦
鄂ICP備17016276號-1