測量平差課程設(shè)計報告

1、東南大學(xué)交通學(xué)院測量平差課程設(shè)計報告設(shè)計題目：專業(yè)：測繪工程專業(yè) 班級：學(xué)號：姓名：指導(dǎo)老師：日期：目錄1.課程設(shè)計目的22.課程設(shè)計任務(wù)23. 課程設(shè)計重點以及基本要求24. 課程設(shè)計具體要求 25. 課程設(shè)計案例及分析36.課程設(shè)計展示成果107.課程設(shè)計源代碼 128.課程設(shè)計心得體會 17 1、課程設(shè)計目的誤差理論與測量平差是一門理論與實踐并重的課程，其課程設(shè)計是測量數(shù)據(jù)處理理論學(xué)習(xí)的一個重要的實踐環(huán)節(jié)，它是在我們學(xué)習(xí)了專業(yè)基礎(chǔ)課“誤差理論與測量平差基礎(chǔ)”課程后進(jìn)行的一門實踐課程。其目的是增強(qiáng)我們對誤差理論與測量平差基礎(chǔ)理論的理解，牢固掌握測量平差的基本原理和基本公式，熟悉測量數(shù)據(jù)處理

2、的基本技能和計算方法，靈活準(zhǔn)確地應(yīng)用于解決各類數(shù)據(jù)處理的實際問題，并能用所學(xué)的計算機(jī)理論知識，編制簡單的計算程序。2、課程設(shè)計的任務(wù)（1）該課的課程設(shè)計安排在理論學(xué)習(xí)結(jié)束之后進(jìn)行的，主要是平面控制網(wǎng)和高程控制網(wǎng)嚴(yán)密平差。（2）通過課程設(shè)計，培養(yǎng)學(xué)生運(yùn)用本課程基本理論知識和技能，分析和解決本課程范圍內(nèi)的實際工程問題的能力，加深對課程理論的理解與應(yīng)用。(3）在指導(dǎo)老師的指導(dǎo)下，要求每個學(xué)生獨(dú)立完成本課程設(shè)計的全部內(nèi)容。3、課程設(shè)計重點以及基本要求課程設(shè)計要求每一個學(xué)生必須遵守課程設(shè)計的具體項目的要求，獨(dú)立完成設(shè)計內(nèi)容，并上交設(shè)計報告。在學(xué)習(xí)知識、培養(yǎng)能力的過程中，樹立嚴(yán)謹(jǐn)、求實、勤奮、進(jìn)取的良好學(xué)

3、風(fēng)。課程設(shè)計前學(xué)生應(yīng)認(rèn)真復(fù)習(xí)教材有關(guān)內(nèi)容和測量平差課程設(shè)計指導(dǎo)書，務(wù)必弄清基本概念和本次課程設(shè)計的目的、要求及應(yīng)注意的事項，以保證保質(zhì)保量的按時完成設(shè)計任務(wù)。本次課程設(shè)計重點是培養(yǎng)我們正確應(yīng)用公式、綜合分析和解決問題的能力，以及計算機(jī)編程能力。另外它要求我們完成12個綜合性的結(jié)合生產(chǎn)實踐的題目，如目前生產(chǎn)實踐中經(jīng)常用到的水準(zhǔn)網(wǎng)嚴(yán)密平差及精度評定，邊角網(wǎng)（導(dǎo)線）嚴(yán)密平差及精度評定等。4、課程設(shè)計具體設(shè)計項目內(nèi)容總體思路：現(xiàn)有等級水準(zhǔn)網(wǎng)的全部觀測數(shù)據(jù)及網(wǎng)型、起算數(shù)據(jù)。要求對該水準(zhǔn)網(wǎng)，分別用條件、間接兩種方法進(jìn)行嚴(yán)密平差，并進(jìn)行平差模型的正確性檢驗。水準(zhǔn)網(wǎng)的條件平差：列條件平差值方程、改正數(shù)條件方程

4、、法方程；利用自編計算程序解算基礎(chǔ)方程，求出觀測值的平差值、待定點的高程平差值；評定觀測值平差值的精度和高程平差值的精度。進(jìn)行平差模型正確性的假設(shè)檢驗。水準(zhǔn)網(wǎng)的間接平差：列觀測值平差值方程、誤差方程、法方程；利用自編計算程序解算基礎(chǔ)方程，求出觀測值的平差值、待定點的高程平差值；評定觀測值平差值的精度和高程平差值的精度。進(jìn)行平差模型正確性的假設(shè)檢驗。平面控制網(wǎng)（導(dǎo)線網(wǎng)）嚴(yán)密平差及精度評定總體思路：現(xiàn)有等級導(dǎo)線網(wǎng)的全部觀測數(shù)據(jù)及網(wǎng)型、起算數(shù)據(jù)。要求對該導(dǎo)線網(wǎng)，用間接進(jìn)行嚴(yán)密平差，并進(jìn)行平差模型的正確性檢驗。報告的編寫對手工解算控制網(wǎng)進(jìn)行程序驗證，編寫課程設(shè)計報告。報告應(yīng)包括起算數(shù)據(jù)、控制網(wǎng)圖形、平

5、差結(jié)果、精度指標(biāo)，點位誤差橢圓圖等。成果應(yīng)以表格形式給出，封面統(tǒng)一格式見附錄。5. 課程設(shè)計案例及分析如圖所示水準(zhǔn)網(wǎng)，a、b兩點為高程已知，各觀測高差及路線長度如表1。已知數(shù)據(jù) 表1高差觀測值(m)對應(yīng)線路長度(km)已知點高程（m）h1 = 1.359h2 = 2.009h3 = 0.363h4 =-0.640h5 = 0.657h6 = 1.000h7 = 1.6501122112h1= 35.000h2= 36.000要求：按條件以及間接平差法分別求：（1） 待定點高程平差值；（2） 待定點高程中誤差；（3） p2和p3點之間平差后高差值的中誤差；（4） 平差模型正確性檢驗（四等水準(zhǔn)測量

6、每公里高差觀測中誤差為5毫米）。5.1 水準(zhǔn)網(wǎng)條件平差求平差值、列條件方程：由題意可知：n=7，t=4，r=n-t=3.觀測方程為線性化得條件方程：其中系數(shù)陣a=，w=組建法方程：令1km的觀測高差為單位權(quán)觀測，即，又。即，且法方程為，其中：，由此可得法方程：解算法方程由矩陣運(yùn)算程序可得.4 計算改正數(shù)。由矩陣運(yùn)算程序可得：計算平差值，得：又可得：精度評定，其中，又由矩陣運(yùn)算程序可得又得：平差模型的正確性檢驗原假設(shè)和備選假設(shè)為又計算統(tǒng)計量以自由度=4，=0.05查得的分布表：可見，在內(nèi)，該平差模型正確，平差結(jié)果可用。5.2水準(zhǔn)網(wǎng)間接平差平差值 列誤差方程選取、三點的高程、為參數(shù)，

7、由題意知可列出七個平差值方程 且代入可得誤差方程 可得系數(shù)陣組成法方程取2km的觀測高差為單位權(quán)觀測，即有權(quán)陣，其中由矩陣運(yùn)算程序可得：,由可得法方程解算法方程并由矩陣運(yùn)算程序可得:計算改正數(shù)計算平差值，得：又可得：精度評定,又且故有：平差模型的正確性檢驗原假設(shè)和備選假設(shè)為又計算統(tǒng)計量以自由度=4，=0.05查得的分布表：可見，在內(nèi)，該平差模型正確，平差結(jié)果可用。6.課程設(shè)計展示成果檢驗過程部分截圖如下7.課程設(shè)計源代碼using system;usingsystem.collections.generic;usingsystem.componentmodel;usingsystem.data

8、;usingsystem.drawing;usingsystem.linq;usingsystem.text;usingsystem.threading.tasks;usingsystem.windows.forms;publicpartialclassform1 : form public form1() initializecomponent(); privatevoid label12_click(object sender, eventargs e) privatevoid button3_click(object sender, eventargs e) application.ex

9、it(); privatevoid button2_click(object sender, eventargs e) foreach (control ctrl in controls)/或為groupbox1.controls/panel1.controls if (ctrl istextbox)ctrl.text = ; privatevoid textbox1_textchanged(object sender, eventargs e) privatevoid textbox2_textchanged(object sender, eventargs e) privatevoid t

10、extbox3_textchanged(object sender, eventargs e) privatevoid textbox4_textchanged(object sender, eventargs e) privatevoid textbox5_textchanged(object sender, eventargs e) privatevoid textbox6_textchanged(object sender, eventargs e) privatevoid textbox7_textchanged(object sender, eventargs e) privatev

11、oid label10_click(object sender, eventargs e) privatevoid textbox8_textchanged(object sender, eventargs e) privatevoid textbox9_textchanged(object sender, eventargs e) privatevoid textbox10_textchanged(object sender, eventargs e) privatevoid button1_click(object sender, eventargs e) double ax = doub

12、le.parse(textbox1.text);double ay = double.parse(textbox2.text);doublebx = double.parse(textbox3.text);double by = double.parse(textbox4.text);double a1 = double.parse(textbox5.text);double a2 = double.parse(textbox6.text);double a3 = double.parse(textbox7.text);double b1 = double.parse(textbox8.tex

13、t);double b2 = double.parse(textbox9.text);double b3 = double.parse(textbox10.text);double m, n, px, py; m = (a1 + a2 / 60 +a3 / 3600) * (math.pi / 180); n = (a1 + a2 / 60 + a3 / 3600) * (math.pi / 180);px = (ax * 1 / (math.tan(n) + bx * 1 / (math.tan(m) - (ay - by) / (1 / (math.tan(n) + 1 / (math.t

14、an(m);using system;usingsystem.collections.generic;usingsystem.componentmodel;usingsystem.data;usingsystem.drawing;usingsystem.linq;usingsystem.text;usingsystem.threading.tasks;usingsystem.windows.forms;public partial class form1 : form public form1() initializecomponent(); private void label12_clic

15、k(object sender, eventargs e) private void button3_click(object sender, eventargs e) application.exit(); private void button2_click(object sender, eventargs e) foreach (control ctrl in controls)/或為groupbox1.controls/panel1.controls if (ctrl is textbox)ctrl.text = ; private void textbox1_textchanged(

16、object sender, eventargs e) private void textbox2_textchanged(object sender, eventargs e) private void textbox3_textchanged(object sender, eventargs e) private void textbox4_textchanged(object sender, eventargs e) private void textbox5_textchanged(object sender, eventargs e) private void textbox6_te

17、xtchanged(object sender, eventargs e) private void textbox7_textchanged(object sender, eventargs e) private void label10_click(object sender, eventargs e) private void textbox8_textchanged(object sender, eventargs e) private void textbox9_textchanged(object sender, eventargs e) private void textbox1

18、0_textchanged(object sender, eventargs e) private void button1_click(object sender, eventargs e) double ax = double.parse(textbox1.text);double ay = double.parse(textbox2.text);doublebx = double.parse(textbox3.text);double by = double.parse(textbox4.text);double a1 = double.parse(textbox5.text);doub

19、le a2 = double.parse(textbox6.text);double a3 = double.parse(textbox7.text);double b1 = double.parse(textbox8.text);double b2 = double.parse(textbox9.text);double b3 = double.parse(textbox10.text);double m, n, px, py; m = (a1 + a2 / 60 +a3 / 3600) * (math.pi / 180); n = (a1 + a2 / 60 + a3 / 3600) *

20、(math.pi / 180);px = (ax * 1 / (math.tan(n) + bx * 1 / (math.tan(m) - (ay - by) / (1 / (math.tan(n) + 1 / (math.tan(m);py = (ay * 1 / (math.tan(n) + by * 1 / (math.tan(m) + (ax - bx) / (1 / (math.tan(n) + 1 / (math.tan(m); label15.text = + px; label16.text = + py; private void label2_click(object se

21、nder, eventargs e) private void label3_click(object sender, eventargs e) private void label1_click(object sender, eventargs e) private void label17_click(object sender, eventargs e) private void label15_click(object sender, eventargs e) private void label16_click(object sender, eventargs e) py = (ay

22、 * 1 / (math.tan(n) + by * 1 / (math.tan(m) + (ax - bx) / (1 / (math.tan(n) + 1 / (math.tan(m); label15.text = + px; label16.text = + py; privatevoid label2_click(object sender, eventargs e) privatevoid label3_click(object sender, eventargs e) privatevoid label1_click(object sender, eventargs e) privatevoid label17_click(object sender, eventargs e) privatevoid label1