誤差理論與測(cè)量平差課程設(shè)計(jì)--水準(zhǔn)網(wǎng)嚴(yán)密平差及精度評(píng)定

上傳人：奔*** IP屬地：河北更新時(shí)間：2024-03-01 格式：doc 頁(yè)數(shù)：17 大?。?0.50KB 人氣指數(shù)：12 舉報(bào) 版權(quán)申訴

2023年全國(guó)碩士研究生考試考研英語(yǔ)一試題真題（含答案詳解+作文范文）_第1頁(yè)

已閱讀1頁(yè)，還剩16頁(yè)未讀，繼續(xù)免費(fèi)閱讀

版權(quán)說(shuō)明：本文檔由用戶提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權(quán)，請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡(jiǎn)介

1、　　測(cè)量平差程序設(shè)計(jì)　　課　　程　　設(shè)　　計(jì)

2、　　報(bào)　　告　　課題：水準(zhǔn)網(wǎng)嚴(yán)密平差及精度評(píng)定　　目錄　　一、目錄 ----------------------------1</

3、p>　　二、序言 ---------------------------- 2　　三、設(shè)計(jì)思路 ------------------------ 3　　四、程序流程圖 ---------------------- 4　　五、程序及說(shuō)明 ---------------------- 5&

4、lt;/p>　　六、計(jì)算結(jié)果 -----------------------12　　七、總結(jié) --------------------------- 15　　第二部分序言　　1、課程設(shè)計(jì)的性質(zhì)、目的和任務(wù)

5、;　　誤差理論與測(cè)量平差是一門(mén)理論與實(shí)踐并重的課程，其課程設(shè)計(jì)是測(cè)量數(shù)據(jù)處理理論學(xué)習(xí)的一個(gè)重要的實(shí)踐環(huán)節(jié)，它是在我們學(xué)習(xí)了專(zhuān)業(yè)基礎(chǔ)課“誤差理論與測(cè)量平差基礎(chǔ)”課程后進(jìn)行的一門(mén)實(shí)踐課程。其目的是增強(qiáng)我們對(duì)誤差理論與測(cè)量平差基礎(chǔ)理論的理解，牢固掌握測(cè)量平差的基本原理和基本公式，熟悉測(cè)量數(shù)據(jù)處理的基本技能和計(jì)算方法，靈活準(zhǔn)確地應(yīng)用于解決各類(lèi)數(shù)據(jù)處理的實(shí)際問(wèn)題，并能用所學(xué)的計(jì)算機(jī)理論知識(shí)，編制簡(jiǎn)單的計(jì)算程序。<p&

6、gt;　　2、誤差理論與測(cè)量平差課程和其它課程的聯(lián)系和分工　　這次課程設(shè)計(jì)中所用的數(shù)學(xué)模型和計(jì)算方法是我們?cè)谡`差理論與測(cè)量平差課程中所學(xué)的內(nèi)容，所使用的C程序語(yǔ)言使我們?cè)谟?jì)算機(jī)基礎(chǔ)課程中所學(xué)知識(shí)。誤差理論與測(cè)量平差課程設(shè)計(jì)是測(cè)量平差和計(jì)算機(jī)程序設(shè)計(jì)等課程的綜合實(shí)踐與應(yīng)用，同時(shí)也為我們今后步入工作崗位打下了一定基礎(chǔ)。　　3、課程設(shè)計(jì)重點(diǎn)及內(nèi)容</p

7、>　　本次課程設(shè)計(jì)重點(diǎn)是培養(yǎng)我們正確應(yīng)用公式、綜合分析和解決問(wèn)題的能力，以及計(jì)算機(jī)編程能力。另外它要求我們完成1－2個(gè)綜合性的結(jié)合生產(chǎn)實(shí)踐的題目。如目前生產(chǎn)實(shí)踐中經(jīng)常用到的水準(zhǔn)網(wǎng)嚴(yán)密平差及精度評(píng)定，邊角網(wǎng)（導(dǎo)線）嚴(yán)密平差及精度評(píng)定等。此次我所選的課程設(shè)計(jì)課題是水準(zhǔn)網(wǎng)嚴(yán)密平差及精度評(píng)定，其具體內(nèi)容如下：　　根據(jù)題目要求，正確應(yīng)用平差模型列出觀測(cè)值條件方程、誤差方程和法方

8、程；解算法方程，得出平差后的平差值及各待定點(diǎn)的高程平差值；評(píng)定各平差值的精度和各高程平差值的精度。　　具體算例為：　　如圖所示水準(zhǔn)網(wǎng)，有2個(gè)已知點(diǎn)，3個(gè)未知點(diǎn)，7個(gè)測(cè)段。各已知數(shù)據(jù)及觀測(cè)值見(jiàn)下表　　已知點(diǎn)高程H1=5.016m ， H2=6.016m （2）高差觀測(cè)值(

9、m)　?。?）求各待定點(diǎn)的高程；3-4點(diǎn)的高差中誤差；3號(hào)點(diǎn)、4號(hào)點(diǎn)的高程中誤差。　　第三部分設(shè)計(jì)思路　　一、解題步驟　?。?）此次設(shè)計(jì)我所采用的模型為間接平差模型，根據(jù)已知條件我們可知觀測(cè)總數(shù)

10、n=7，必要觀測(cè)數(shù)t=3（則多余觀測(cè)數(shù)r=n-t=4），因此我需先選定三個(gè)參數(shù)，即3、4、5點(diǎn)的最或然高程X3、X4、X5（X=X0+x，X30=6.375、X40=7.025、X50=6.611；其中X0為參數(shù)的近似值，x為其改正值）為參數(shù)。　　（2）列出條件方程，即將每一個(gè)觀測(cè)量的平差值分別表達(dá)成所選參數(shù)的函數(shù)，H1+h1=X3、H1+h2=X4、H2+h3=X3、H2+h4=X4、X3+h

11、5=X4、X3+h6=X5、X5+h7=H2；整理后得出誤差方程，v1=x3、v2=x4、v3=x3-4、v4=x4-3、v5=-x3+x4-7、v6=-x3+x5-2、v7=-x5，即v=Bx-l的形式。　?。?）定權(quán)，令每千米的觀測(cè)高差為單位權(quán)觀測(cè)，即Pi=1/Si，從而可寫(xiě)出權(quán)陣P；根據(jù)誤差方程式又可得其系數(shù)矩陣B和自由項(xiàng)l，并由它們組成法方程N(yùn)BBx-W=0（其中NBB=BTPB，W=B

12、TPl），法方程的個(gè)數(shù)等于所選參數(shù)的個(gè)數(shù)。　?。?）解算法方程，求出參數(shù)改正值x并計(jì)算參數(shù)的平差值X=X0+x。　?。?）由誤差方程計(jì)算V，并求出觀測(cè)量的平差值。為了檢查平差計(jì)算的正確性，將所求的值代入條件方程，看其是否滿足方程。　?。?）精度評(píng)定，計(jì)算單位權(quán)中誤差，按照題設(shè)要求列出權(quán)函數(shù)式，再根據(jù)平差參數(shù)的協(xié)方

13、差陣求出協(xié)因數(shù)，最后求出某段高差中誤差，某些點(diǎn)的高程中誤差。　　二、程序設(shè)計(jì)思想　　考慮到在解題過(guò)程中一些計(jì)算的復(fù)雜性，我們需借助一些技術(shù)將計(jì)算簡(jiǎn)單化，快捷化，因此在課程設(shè)計(jì)過(guò)程中，我們把一些C語(yǔ)言程序設(shè)計(jì)引入其中；通過(guò)一些簡(jiǎn)單、明了的程序及子函數(shù)調(diào)用，我們就可以很方便快捷的求出用筆算比較繁瑣、費(fèi)時(shí)的矩陣乘積、矩陣的

14、逆（如BTPB、BTPl）等運(yùn)算。　　第四部分程序流程圖　　↓　　↓　　↓　　↓</b&g

15、t;　　↓　　↓　　↓　　第五部分程序及說(shuō)明　　一、矩陣相乘計(jì)算函數(shù)

16、#include “stdio.h”　　void Matrix(a,b,m,n,k,c)　　int m,n,k;　　double a[],b[],c[];　　{　　int i,j,l,u;&l

17、t;/p>　　for(i=0;i<=m-1;i++)　　for(j=0;j<=k-1;j++)　　{　　u=i*k+j;c[u]=0.0;　　for(l=0;l<=n-1;l++)<

18、/p>　　c[u]=c[u]+a[i*n+l]*b[l*k+j];　　}　　return;　　}　　1.計(jì)算BT

19、P　　main()　　{　　int i,j;　　static double a[3][7]=BT；

20、;　　static double c[3][7],b[7][7]=P；　　Matrixmul(a,b,3,7,7,c);　　printf(“\n”);　　for(i=0;i<=2;i++)　　{<p

21、>　　for(j=0;j<=6;j++)　　printf(“%8.4f\t”,c[i][j];　　printf(“\n”);　　}　　printf(“\n”);

22、　return0;　　}　　2.計(jì)算BTPB，即NBB　　main()　　{

23、　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”);<

24、;p>　　for(i=0;i<=2;i++)　　{　　for(j=0;j<=2;j++)　　printf(“%8.4f\t”,c[i][j];　　printf(“\n”);&

25、lt;b>　　}　　printf(“\n”);　　return0;　　}　　3.計(jì)算BTPl，即W　　main

26、()　　{　　int i,j;　　static double a[3][7]=BTP；　　static double c[3][1],b[7][1]=l；

27、　　Matrixmul(a,b,3,7,1,c);　　printf(“\n”);　　for(i=0;i<=2;i++)　　{　　for(j=0;j<=0;j++)

28、;　　printf(“%8.4f\t”,c[i][j];　　printf(“\n”);　　}　　printf(“\n”);　　return0;<b&

29、gt;　　}　　二、矩陣的逆計(jì)算函數(shù)（求NBB-1）　　#include "stdio.h" 　　#define M 3 　　void main() 　　{ </b&

30、gt;　　float MAT[M][2*M]; 　　float MAT1[M][M]; 　　float t; 　　int i,j,k,l; 　　/***********************

31、************************/ 　　/*對(duì)矩陣進(jìn)行初始化*/ 　　for(i=0;i<M;i++) 　　for(j=0;j<2*M;j++) 　　MAT1[j]='\0'; 　　/*

32、對(duì)MAT1矩陣賦初值 */ 　　for(i=0;i<M;i++) 　　for (j=0;j<M;j++) 　　scanf("%f",&MAT1[j]); 　　/*打印目標(biāo)矩陣?*/ 　　pr

33、intf("原矩陣為：\n"); 　　for (i=0;i<M;i++) 　　{ 　　for (j=0;j<M;j++) 　　printf("%13.7f",MAT1[j]);

34、　　printf("\n"); 　　} /********************************************/ 　　/*對(duì)MAT1矩陣進(jìn)行擴(kuò)展,MAT1矩陣添加單位陣，由M*M變成2M*2M矩陣 */ 　　for(i=0;i<

35、;M;i++) 　　for(j=0;j<2*M;j++) 　　if (j<M) MAT[j]=MAT1[j]; 　　else if (j==M+i) MAT[j]=1; 　　else MAT[j]=0; 　　/*對(duì)M矩

36、陣進(jìn)行變換，使得前半部分矩陣成為單位陣，則 */ 　　/*后半部分矩陣即為所求矩陣逆陣 */ 　　for(i=0;i<M;i++) 　　{ 　　/*對(duì)第i行進(jìn)行歸一化 */ 　　fo

37、r (j=0;j<2*M;j++) 　　for(k=i+1;k<M;k++) 　　MAT[j]=MAT[j]+MAT[k][j]; 　　t=MAT; 　　for(j=i;j<2*M;j++)

38、　　MAT[j]=MAT[j]/t; 　　/*對(duì)矩陣進(jìn)行行變換，使得第i 列只有一個(gè)元素不為零，且為1*/ 　　for(k=0;k<M;k++) 　　if(k!=i) 　　{ </b&g

39、t;　　t=MAT[k]; 　　for (l=i;l<2*M;l++) 　　MAT[k][l]=MAT[k][l]-MAT[l]*t; 　　} 　　} </b&g

40、t;　　/*將后半部分矩陣即所求矩陣逆陣存入MAT2矩陣。*/ 　　for(i=0;i<M;i++) 　　{ 　　for(j=0;j<M;j++) 　　MAT1[j]=MAT[j+

41、M]; 　　printf("\n"); 　　} 　　/*********************************************/ 　　/*輸出所求的逆陣*/ <p&g

42、t;　　printf("逆陣為：\n"); 　　for(i=0;i<M;i++) 　　{ 　　for(j=0;j<M;j++) 　　printf("%8.4f",MAT1[j]);

43、　　printf("\n"); 　　} 　　}　　4.求NBB-1W，即改正數(shù)x　　main()

44、　　{　　int i,j;　　static double a[3][3]=NBB-1；　　static double c[3][1],b[3][1]=W；

45、　　Matrixmul(a,b,3,3,1,c);　　printf(“\n”);　　for(i=0;i<=2;i++)　　{　　for(j=0;j<=0;j++)　　printf(

46、“%8.4f\t”,c[i][j];　　printf(“\n”);　　}　　printf(“\n”);　　return0;　　}<

47、/b>　　5.計(jì)算Bx　　main()　　{　　int i,j;　　st

48、atic double a[7][3]=B；　　static double c[7][1],b[3][1]=x；　　Matrixmul(a,b,7,3,1,c);　　printf(“\n”);　　for(i=0;i<=6;i++)<p&

49、gt;　　{　　for(j=0;j<=0;j++)　　printf(“%8.4f\t”,c[i][j];　　printf(“\n”);　　}　　pr

50、intf(“\n”);　　return0;　　}　　6.計(jì)算VTP　　main()<p&g

51、t;　　{　　int i,j;　　static double a[1][7]=VT；　　static double c[1][7],b[7][7]=P；　　Matrixmul(a,b,1,7,7,

52、c);　　printf(“\n”);　　for(i=0;i<=0;i++)　　{　　for(j=0;j<=6;j++)　　printf(“%8.4f\t”,c[i][j];<

53、/p>　　printf(“\n”);　　}　　printf(“\n”);　　return0;　　}<

54、p>　　7.計(jì)算VTPV　　main()　　{　　int i,j;　　static double a[1][7]=

55、VTP；　　static double c[1][1],b[7][1]=V；　　Matrixmul(a,b,1,7,1,c);　　printf(“\n”);　　for(i=0;i<=0;i++)　　{<

56、;/b>　　for(j=0;j<=0;j++)　　printf(“%8.4f\t”,c[i][j];　　printf(“\n”);　　}　　printf(“\n”);</p&

57、gt;　　return0;　　注：程序中有下劃線部分在C語(yǔ)言環(huán)境中運(yùn)行時(shí)，需根據(jù)已知條件及所求結(jié)果進(jìn)行替換！　　第六部分計(jì)算結(jié)果　　根據(jù)條件方程及定權(quán)原則寫(xiě)出B、l、P及BT

58、;　　B={{1.0，0.0，0.0},　　{0.0，1.0，0.0}，　　{1.0，0.0，0.0}，　　{0.0，1.0，0.0}，　　{-1.0，1.0，0.0}，　　{-1.0，0.0，1.0}，&

59、lt;p>　　{0.0，0.0，-1.0}}　　l={{0.0},　　{0.0}，　　{4.0}，　　{3.0}，

60、　　{7.0}，　　{2.0}，　　{0.0}}　　P={{0.9091,0.0,0.0,0.0,0.0,0.0,0.0},<p

61、>　　{0.0,0.5882,0.0,0.0,0.0,0.0,0.0}，　　{0.0,0.0,0.4348,0.0,0.0,0.0,0.0}，　　{0.0,0.0,0.0,0.3704,0.0,0.0,0.0}，　　{0.0,0.0,0.0,0.0,0.4167,0.0,0.0}，&

62、lt;p>　　{0.0,0.0,0.0,0.0,0.0,0.7143,0.0}，　　{0.0,0.0,0.0,0.0,0.0,0.0,0.3846}}　　BT={{1.0,0.0,1.0,0.0,-1.0,-1.0,0.0},　　{0.0,1.0,0.0,1.0,1.0,0.0,0.0}，

63、;　　{0.0,0.0,0.0,0.0,0.0,1.0,-1.0}}　　一、在矩陣相乘計(jì)算函數(shù)的程序前提下，進(jìn)行以下子程序的調(diào)用　　1.替換第1個(gè)程序中的BT 、P并運(yùn)行程序得到BTP 　　BTP={{0.9091,0.0,0.4348,0.0,-0.4167,-0.7143,0.0},&l

64、t;/p>　　{0.0,0.5882,0.0,0.374,0.4167,0.0,0.0}，　　{0.0,0.0,0.0,0.0,0.0,0.7143,-0.3846}}　　2.替換第2個(gè)程序中的BTP、B并運(yùn)行程序得到BTPB，即NBB　　NBB={{2.4748，-0.4167，-0

65、.7143},　　{-0.4167，1.3753，0.0}，　　{-0.7143，0.0，1.0989}}　　3.替換第3個(gè)程序中的BTP、l并運(yùn)行程序得到BTPl，即W　　W={{-2.6063},　　{4.0

66、281}，　　{1.4286}}　　二、在矩陣的逆計(jì)算函數(shù)程序中進(jìn)行以下操作　　運(yùn)行程序，按照提示及以上運(yùn)算得到的矩陣NBB輸入其元素，運(yùn)行的結(jié)果即為NBB-1 NBB-1={{0.5307，0.1608，0.3450},<p

67、>　　{0.1608，0.7758，0.1045}，　　{0.3450，0.1045，1.1342}}　　三、再次在矩陣相乘計(jì)算函數(shù)的程序前提下，進(jìn)行以下子程序的調(diào)用　　1.替換第4個(gè)程序中的NBB-1、W并運(yùn)行程序得到NBB-1W，即所選參數(shù)的改正數(shù)x　　x={{

68、-0.2426},　　{2.8552}，　　{1.1421}}　　2.替換第5個(gè)程序中的B、x并運(yùn)行程序得到Bx　　Bx={{-0.2426},<b

69、>　　{2.8552}，　　{-0.2464}，　　{2.8552}，　　{3.0978}，　　{1.3847}，

70、;　　{-1.1421}}　　3.根據(jù)V=Bx-l求出各觀測(cè)值的改正數(shù)V，并寫(xiě)出VT，然后替換第6個(gè)程序中的VT、P并運(yùn)行程序得到VTP　　V={{-0.2426},　　{2.8552}，　　{-4.2426}，

71、　　{-0.1448}，　　{-3.9022}，　　{-0.6153}，　　{-1.1421}}　　VT=｛{-0.2426,2.8552,-4.2426,-0.1448,-3.9022,-0.6153,-1.1421}｝<

72、/p>　　VTP={{-0.2205,1.6794,-1.8447,-0.0536,1.6260,-0.4395,-0.4393}}　　4.替換第7個(gè)程序中的VTP、V并運(yùn)行程序得到VTPV　　VTPV=19.7997　　四、求出各個(gè)觀測(cè)值平差值并按要求平定精度

73、　　X3=6.3748 m X4=7.0279 m X5=6.6122 m　　h1=1.3588m h2=2.0119m h3=0.3588m h4=1.0119m h5=0.6531m h6=0.2374m h7=-0.5961m　　根據(jù)公式可求得單位權(quán)中誤差為 2.225mm<p&g

74、t;　　h34= X3- X4 Q34=[1 -1 0] NBB-1[1 -1 0]T=0.9849　　H3=X3 Q34=[1 0 0] NBB-1[1 0 0]T=0.5307　　H4= X4 Q34=[0 1 0] NBB-1[0 1 0]T=0.7758　　3、4點(diǎn)高差中誤差為 2.208mm<

75、;/p>　　3號(hào)點(diǎn)高程中誤差為 1.621mm　　4號(hào)點(diǎn)高程中誤差為 1.96mm　　第七部分總結(jié)　　通過(guò)這次誤差理論與測(cè)量平差的課程設(shè)計(jì)，我又對(duì)整本書(shū)有了一個(gè)更深的理解。其實(shí)課程設(shè)計(jì)就是將我們所學(xué)的理論知識(shí)應(yīng)用于實(shí)踐的過(guò)程，在這一過(guò)程中，進(jìn)一步

76、掌握測(cè)量平差的基本原理和基本公式，并熟悉測(cè)量數(shù)據(jù)處理的基本技能和計(jì)算方法。　　或許我們已對(duì)《誤差理論與測(cè)量平差》這本書(shū)的理論知識(shí)有了一定了解，但將它應(yīng)用于實(shí)踐依然是我們的一個(gè)難點(diǎn)，尤其是將這門(mén)課程與計(jì)算機(jī)程序完美地結(jié)合。這便要求我們?cè)谠械慕忸}思路中加入C語(yǔ)言程序，并讓它來(lái)幫助我們解決矩陣的復(fù)雜運(yùn)算。既然用到了程序，我們就必須保證其運(yùn)算的簡(jiǎn)潔性、正確性，尤其是在編寫(xiě)過(guò)程中要認(rèn)真檢查，為程序順利運(yùn)

77、行打下基礎(chǔ)。另外在各個(gè)子程序調(diào)用過(guò)程中，我們要充分考慮其順序性并反復(fù)調(diào)試，以便得到理想結(jié)果。　　盡管在這次課程設(shè)計(jì)中遇到了很多困難，但我卻得到了不少收獲，并培養(yǎng)了自己正確應(yīng)用公式、綜合分析和解決問(wèn)題的能力，同時(shí)也為今后步入社會(huì)打下了一定的基礎(chǔ)。另外，我們還要學(xué)會(huì)綜合利用自身所學(xué)的知識(shí)，并將它們聯(lián)系起來(lái)幫助自己有效地解決實(shí)際中的問(wèn)題。　　總之，在這次課程設(shè)

眾賞文庫(kù)> 全部分類(lèi)> 畢業(yè)設(shè)計(jì)

溫馨提示

1. 本站所有資源如無(wú)特殊說(shuō)明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
3. 本站RAR壓縮包中若帶圖紙，網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽，若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
5. 眾賞文庫(kù)僅提供信息存儲(chǔ)空間，僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理，對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯，并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容，請(qǐng)與我們聯(lián)系，我們立即糾正。
7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。

誤差理論與測(cè)量平差課程設(shè)計(jì)--水準(zhǔn)網(wǎng)嚴(yán)密平差及精度評(píng)定

文檔簡(jiǎn)介

溫馨提示

最新文檔

評(píng)論

誤差理論與測(cè)量平差課程設(shè)計(jì)--水準(zhǔn)網(wǎng)嚴(yán)密平差及精度評(píng)定

文檔簡(jiǎn)介

溫馨提示

最新文檔

評(píng)論

免費(fèi)下載