測(cè)量平差課程設(shè)計(jì)報(bào)告

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

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

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

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

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

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

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