MATLAB測量平差程序?qū)嵙?xí)報告

1、測量平差程序設(shè)計課程設(shè)計任務(wù)書專業(yè)班級：_ _ _指導(dǎo)教師：_ _小組成員： 目 錄設(shè)計題目1設(shè)計資料:1一、課程設(shè)計的目的2二、課程設(shè)計的任務(wù)和內(nèi)容2三、課程設(shè)計階段2四、組織方式進度安排3五、考核與成績評定4六、參考文獻：4七、實習(xí)報告：5設(shè)計題目 邊角三角網(wǎng)平差程序設(shè)計設(shè)計資料:表1 已知點坐標(biāo)點號X坐標(biāo)Y坐標(biāo)備注Bj011149.5525529.757已知點Bj021351.0826907.547已知點Bj031466.8006231.200待定點Bj041548.2005403.600待定點表2 觀測值精度邊長觀測值固定誤差(mm)1比例誤差(ppm)1方向觀測值()1表3邊長觀測值

2、起點終點邊長(m)Bj03Bj04832.388Bj01Bj03770.170Bj03Bj02685.482Bj02Bj011392.450Bj02Bj041516.951Bj01Bj04418.530表4 方向觀測值起點終點方向觀測值(°)起點終點方向觀測值(°)Bj01Bj030Bj03Bj010Bj01Bj0215.57244Bj03Bj0429.57251Bj01Bj04276.43169Bj03Bj02213.5653Bj02Bj010Bj04Bj030Bj02Bj0415.48098Bj04Bj021.48093Bj02Bj0317.59294Bj04Bj016

3、6.45528一、課程設(shè)計的目的學(xué)生在學(xué)習(xí)完誤差理論與測量平差基礎(chǔ)、測量平差程序設(shè)計基礎(chǔ)等課程的基礎(chǔ)上，設(shè)計一個完整的測量數(shù)據(jù)處理程序，培養(yǎng)學(xué)生綜合應(yīng)用量數(shù)據(jù)處理與計算機應(yīng)用能力，培養(yǎng)學(xué)生主動學(xué)習(xí)，創(chuàng)新設(shè)計能力。二、課程設(shè)計的任務(wù)和內(nèi)容1.課程設(shè)計任務(wù)：在兩周的時間內(nèi)應(yīng)用者Matlab程序設(shè)計語言編制一個完整的邊角網(wǎng)嚴(yán)密平差程序,要求有簡易的界面，數(shù)據(jù)輸入采用文本輸入，采用間接平差模型完成平差的基本計算，能夠畫出控制網(wǎng)圖，輸出基本的計算結(jié)果，并根據(jù)設(shè)計過程完成設(shè)計報告。2. 程序設(shè)計主要內(nèi)容包括：1. 系統(tǒng)功能設(shè)計2. 界面設(shè)計3. 流程設(shè)計4. 代碼書寫5. 程序調(diào)試三、課程設(shè)計階段1. 準(zhǔn)

4、備階段研究設(shè)計任務(wù)書，分析設(shè)計題目，熟悉原始數(shù)據(jù)，明確設(shè)計內(nèi)容和要求；制定課程設(shè)計計劃和進度。2. 熟悉算法模型閱讀誤差理論與測量平差基礎(chǔ)教材，掌握平面控制網(wǎng)數(shù)據(jù)處理的數(shù)學(xué)模型，這里主要是指方向觀測量、角度觀測量、邊長觀測量的觀測方程和誤差方程的構(gòu)成，研究平面觀測數(shù)據(jù)的組織方法，設(shè)計Matlab算法，實現(xiàn)計算的自動表達。3. 功能設(shè)計階段 設(shè)計程序要實現(xiàn)的功能平差程序的基本功能包括數(shù)據(jù)的輸入，平差計算，精度評定、成果輸出等；4.流程和界面設(shè)計階段1. 根據(jù)平差計算的過程和程序功能，畫出流程圖，設(shè)計簡易界面實現(xiàn)數(shù)據(jù)的輸入和平差計算和成果輸出。在此基礎(chǔ)上，根據(jù)功能要求，設(shè)計簡便的界面。5代碼書寫和

5、調(diào)試階段按照計算流程圖和界面設(shè)計，根據(jù)方向觀測值，邊長觀測值的誤差方程的組成，設(shè)計Matlab算法，實現(xiàn)誤差方程的自動構(gòu)成，分階段書寫代碼，調(diào)試實現(xiàn)各個階段的功能。6.設(shè)計報告撰寫階段設(shè)計報告是對整個設(shè)計過程進行綜合總結(jié)提高，內(nèi)容包括課設(shè)的目的意義、程序設(shè)計的內(nèi)容、算法設(shè)計、設(shè)計心得等根據(jù)設(shè)計過程和對測量數(shù)據(jù)處理以及程序設(shè)計的理解進行獨立撰寫。四、組織方式進度安排以小組為單位，每小組5-6人，分工合作共同完成程序設(shè)計任務(wù)，時間兩周，進度安排如下：第一周周一 熟悉設(shè)計書內(nèi)容，制定計劃 熟悉平差模型周二 功能設(shè)計 流程設(shè)計 繪制流程圖 觀測值數(shù)據(jù)組織周三 熟悉觀測方程 誤差方程的構(gòu)成，列出計算公式

6、，設(shè)計Matlab算法 周四 界面設(shè)計 按照分工進行分模塊代碼設(shè)計 完成數(shù)據(jù)讀入模塊 周五 分模塊進行調(diào)試 完成 誤差方程構(gòu)成模塊 第二周周一 代碼書寫 調(diào)試 完成平差計算模塊周二 代碼書寫 調(diào)試 完成輸出模塊周三 代碼書寫 調(diào)試 完成整體調(diào)試 周四 撰寫設(shè)計報告 周五 綜合考核五、考核與成績評定1. 考核內(nèi)容根據(jù)考勤、課程設(shè)計報告、程序結(jié)構(gòu)、程序代碼進行綜合評定。2. 成績評定方法平時成績20%，設(shè)計報告30%；程序代碼50%（程序20%，運行正確30%），總評成績分為優(yōu)、良、中、差四個等級。設(shè)計報告內(nèi)容完整，圖表清晰，代碼書寫規(guī)范，計算結(jié)果正確評為優(yōu)秀。六、參考文獻：1. 武漢大學(xué)測繪學(xué)院

7、 誤差理論與測量平差基礎(chǔ) 20092. 北京建筑工程學(xué)院 測量平差計算 自編資料 20113. 同濟大學(xué) Matlab 控制網(wǎng)測量平差 2006七、實習(xí)報告：目的與任務(wù)：目的：基于學(xué)過的測量平差理論的基礎(chǔ)知識，在matlab中編寫一個相對完整的數(shù)據(jù)平差處理程序，旨在讓同學(xué)們能夠?qū)W(xué)到的平差理論知識和計算機的MATLAB軟件程序編寫有機地結(jié)合在起來，讓同學(xué)在實踐過程中更加深刻地理解理論知識并嘗試去應(yīng)用，由此提高每個人學(xué)生的思考和動手能力。任務(wù)：主要任務(wù)就是在實習(xí)的兩周時間內(nèi)，分組完成一套完整的關(guān)于邊角三角網(wǎng)的平差數(shù)據(jù)處理的的程序代碼編寫與試調(diào)任務(wù)，其中主要包括：數(shù)據(jù)組織與輸入、平差計算、精度評定

8、，成果可視化、平差成果輸出等。然后小組上交一份實習(xí)總結(jié)和小組成員各上交一份個人總結(jié)（設(shè)計書及總結(jié)）。算法與模型：算法：間接平差方法；模型：平面控制網(wǎng)中的邊角三角網(wǎng)本次的平差數(shù)據(jù)處理是采用間接平差的方法（數(shù)學(xué)模型）而進行的，就三個待定點xy坐標(biāo)改正數(shù)及涉及到的五個觀測站的方位角改正數(shù)為未知數(shù)（11個），進行誤差方程的編寫和平差計算，精度評定等。而實際應(yīng)用是采用所提供的平面控制網(wǎng)中的邊角三角網(wǎng)的平差數(shù)據(jù)模型，分別是就著邊長觀測量和方向角觀測量的觀測數(shù)據(jù)進行平差程序的編寫。程序功能設(shè)計：通過手動地選擇數(shù)據(jù)文件，而后自動地對數(shù)據(jù)文件中涉及的變量進行賦值和所獲得的數(shù)據(jù)的輸出，而后通過實際情況分析得出需要

9、進行平差計算和精度評定等的必要步驟，如本次數(shù)據(jù)平差的必要步驟就是待定點近似坐標(biāo)的計算，緊接著就是誤差方程的編寫，包括系數(shù)矩陣和觀測權(quán)陣的求得，而后就是利用最小二乘法的平差計算和精度評定，并且計算誤差橢圓的參數(shù)和該橢圓的自動繪制，最后平差成果的輸出。分模塊講解：數(shù)據(jù)組織輸入和變量賦值已知點坐標(biāo)信息：包括已知點數(shù)n0，待定點數(shù)nd，點號名pn和已知點坐標(biāo)xy0 觀測值精度信息：方向觀測值精度md，邊長觀測值精度ma ,mb邊長觀測值信息：包括起點bf1, 終點bf2 ,邊長觀測值L2方向觀測值：包括測站點df1 ,照準(zhǔn)點df2, 方向觀測值L1 具體代碼截圖：待定點近似坐標(biāo)的計算通過所提供的已知點

10、坐標(biāo)計算得出的已知方位角，在就方向觀測值，求出其每一次觀測的夾角，而后找出每個測站所擁有的已知點（一個觀測方向的兩個點）確定其所在觀測的行數(shù)與其他待定點所在觀測方向（必須有一個點是已知點）的關(guān)系，求出其他待定點所在觀測方向的方位角（注意方位角是否合理問題），而后在邊長觀測值中找到對應(yīng)的邊長，利用極坐標(biāo)法最終求得待定點的近似坐標(biāo)。x0y0. 同時還可以將近似坐標(biāo)以文本的文件輸出，計算代碼截圖如下：誤差方程式組成本次平面控制網(wǎng)平差模型中觀測量為方向觀測量和邊長觀測量，未知數(shù)分為待定點坐標(biāo)和定點角兩類未知數(shù)，邊角網(wǎng)誤差方程的形式如下：方向觀測量：誤差方程：（系數(shù)為：b1,c1）式中常數(shù)項：變長觀測量

11、：誤差方程：（系數(shù)為：b2,c2）其中誤差方程的編寫要注意兩個問題：系數(shù)的位置確定和是否待定點的判斷。對于方向觀測值，未知數(shù)不僅包括各待定點坐標(biāo)（6個），同樣包括各測站的定向角，在未知數(shù)排序中，定向角未知數(shù)在前，待定點坐標(biāo)在后，則未知數(shù)個數(shù)為nz+2*nd=11，對于方向觀測值，每一列與定向角對應(yīng)的系數(shù)均為-1，與坐標(biāo)對應(yīng)的系數(shù) 即每一列有5個系數(shù)（最多），其他系數(shù)均為0，及對應(yīng)的系數(shù)矩陣的數(shù)值的確定（注意正負號）。而對于邊長觀測值，每一列包括邊的兩個點的對應(yīng)的改正數(shù)（如果兩個點都是待定點的話）四個參數(shù)，其他均為0，同樣還要注意系數(shù)的位置問題。由于這里分別對方向觀測量和邊長觀測量設(shè)置兩個數(shù)組，

12、分別存儲系數(shù)陣常數(shù)項變量，這里設(shè)置b1 b2 c1 c2四個數(shù)組。根據(jù)觀測值序列和未知數(shù)序列確定系數(shù)矩陣中各元素的位置和大??；這里將誤差方程構(gòu)成設(shè)計成一個function函數(shù)文件，輸入量為已知點信息，近似坐標(biāo)，邊長觀測量，方向觀測量，等輸出量為誤差方程系數(shù)，常數(shù)項。代碼編寫截圖：5.平差計算和精度評定本次平差的特殊就是結(jié)合兩個誤差方程的數(shù)據(jù)進行同時平差，具體步驟在得到誤差方程的系數(shù)之后，根據(jù)方向觀測值和邊長觀測值的初步計算觀測值權(quán)陣（p1 p2），嚴(yán)格的精確計算可以采用方差分量估計方法精確計算方向和邊長觀測值的權(quán)。在得到誤差方程之后，可以直接利用矩陣進行方程組的解算(最小二乘法)，誤差方程分為

13、方向觀測誤差方程，邊長觀測誤差方程，兩者可分別計算，然后再相加。由于近似坐標(biāo)不精確，在求取參數(shù)時，程序設(shè)置了迭代計算并確定終止條件，在改正數(shù)滿足條件后，即可進行平差值計算和精度統(tǒng)計。程序截圖：6.誤差橢圓參數(shù)計算參考書上的誤差橢圓參數(shù)計算模型由于誤差橢圓描述了點位精度及其在各個方向的分布，表示了網(wǎng)中點或點之間的誤差分布情況，常用誤差橢圓對布網(wǎng)方案做精度分析，誤差橢圓三個參數(shù)為橢圓的長半軸E，短半軸F和主軸方向 本次參數(shù)計算是通過單位權(quán)中誤差，協(xié)因素陣，待定點的協(xié)方差陣等上一步平差得出的結(jié)果或推算結(jié)果，來進行誤差橢圓確定。代碼截圖：7.控制網(wǎng)圖形繪制和誤差橢圓繪制邊角三角網(wǎng)控制網(wǎng)圖形繪制主要考慮

14、繪制控制網(wǎng)點位分布圖，點位標(biāo)注，邊長觀測量和方向觀測量的繪制等。基本思路是首先根據(jù)坐標(biāo)繪制控制網(wǎng)散點圖，并進行標(biāo)注，然后根據(jù)距離觀測值，方向觀測值連線，構(gòu)成完整的控制網(wǎng)圖。至于誤差橢圓的繪制則是在控制網(wǎng)圖的基礎(chǔ)上再加上待定點對應(yīng)的誤差橢圓的繪制就完成了。（利用Plot命令）其中有一個特別注意的地方就是參數(shù)的設(shè)置，要根據(jù)實際情況而定。繪制代碼截圖：8.平面控制網(wǎng)成果輸出以生成txt格式文檔的形式進行成果輸出。其內(nèi)容主要包括五個方面內(nèi)容：平差計算后的坐標(biāo)值，待定點點位中誤差，誤差橢圓參數(shù)計算等。信息具體內(nèi)容如下：1.平差計算成果成果總輸出；2.坐標(biāo)平差成果輸出3.點位誤差橢圓參數(shù)4.誤差橢圓的繪制

15、5.方向觀測平差成果：方向觀測值，平差值，改正數(shù)；6.邊長觀測平差成果：邊長觀測值，邊長平差值，改正數(shù)；代碼截圖（其中的幾個）：平差主界面的設(shè)計 此界面為本組實習(xí)成果的最大創(chuàng)新之一，就是我們在原本要求的平差程序界面（邊角三角網(wǎng)平差程序）的基礎(chǔ)上增加了另一個界面，就是小組信息的界面（包括有小組成員照），而兩個界面之間利用一個總的界面（小組平差程序設(shè)計）進行選擇性打開代碼截圖：主界面小組平差程序設(shè)計總的設(shè)計流程圖：輸入已知點信息和觀測信息數(shù)據(jù)編輯閉合差檢驗合格計算近似坐標(biāo)計算坐標(biāo)改正值及驗后單位權(quán)方差組成法方程并求逆組成邊長觀測誤差方程組成方向觀測方程判斷參數(shù)改正數(shù)精度評定輸出計算結(jié)果NOYES代

16、碼書寫與調(diào)試：由于上面已經(jīng)展示過幾乎所有的代碼（根據(jù)每一個模塊），真正運行整一個程序的時候卻是按照步驟一步步地一個模塊地運行，但是問題一定會出現(xiàn)的，而且還不是很少，關(guān)鍵就是怎么樣通過一步步地發(fā)現(xiàn)問題所在，并根據(jù)實際情況去試調(diào)程序，合理地按照平差理論去運行程序。一下列舉了一些出現(xiàn)的主要問題。方位角 由于所提供的數(shù)據(jù)時邊角三角網(wǎng)的數(shù)據(jù)，這樣就會涉及到方位角問題，而且還不止一次的出現(xiàn)方位角的運算（如待定點的近似坐標(biāo)的求解過程中出現(xiàn)過），這就應(yīng)該驚醒一個重要的細節(jié)：就是所求的方位角是否合理？怎么樣去判斷呢？就是通過if（elseif）語句去判斷其與0和360度或者2*pi 的關(guān)系，在去改正方位角。原則

17、上一旦出現(xiàn)方位角的求解出都是要用這個語句去判斷，這樣才能保證精密。和提高正確性。單位問題單位問題出現(xiàn)在兩處，就是觀測值的單位權(quán)的求解處，二是誤差方程的系數(shù)求解處。由于這兩者是在后面的法方程求解時都要用上所以，在此之前要保證兩者的單位相對統(tǒng)一，這個統(tǒng)一要分兩處，邊長觀測量方程系數(shù)和方向觀測方程系數(shù)，特別注意的是原始數(shù)據(jù)提供的觀測精度（兩種）的單位，和誤差方程的系數(shù)的單位，而在實際中，我們組應(yīng)用的是方向觀測用秒為單位，而變成觀測以米為單位，這樣出來的結(jié)果最后的平差數(shù)據(jù)處理，如精度評定，誤差橢圓計算等，也要注意單位問題。比例參數(shù)根據(jù)實際求出的誤差橢圓參數(shù)（E F fi）的大小和所繪制的三角網(wǎng)的大小，

18、這兩者之間通過比例參數(shù)的調(diào)整才能更好地在顯示誤差橢圓。弧度轉(zhuǎn)角度對于弧度轉(zhuǎn)角度這個小問題，我不得不提，這是一個關(guān)鍵問題。為什么？因為我們組把其他問題都解決了后， 組的坐標(biāo)改正數(shù)，誤差橢圓參數(shù)等都很不正常的顯示在電腦屏幕上，一次次地打擊我們的信心。最后，我在仔細地對照書本，細心查詢，觀察異同發(fā)現(xiàn)，一個最不關(guān)緊要的差異，就是弧度轉(zhuǎn)角度的小問題，我們原來采用的是MATLAB系統(tǒng)自帶的弧度轉(zhuǎn)角度的函數(shù)，結(jié)果我們發(fā)現(xiàn)這個問題直接就是導(dǎo)致待定點的近似坐標(biāo)差異特別大，尤其在待定點3的x,y坐標(biāo)上面，差了好幾米（跟提供的平差結(jié)果數(shù)據(jù)相比），這樣后面的就影響到誤差方程系數(shù)的不正常，和平差計算得出的改正數(shù)，單位全

19、中誤差，點位誤差和誤差橢圓參數(shù)的異常。而這樣我不得不自己再編兩個弧度和角度相互轉(zhuǎn)換的函數(shù)文件。于是重新采用轉(zhuǎn)換函數(shù)程序之后，我們的數(shù)據(jù)就非常正常的出來的。平差程序設(shè)計成果展示進入界面導(dǎo)入數(shù)據(jù)文件數(shù)據(jù)文件的查詢查詢結(jié)果：近似坐標(biāo)的計算計算結(jié)果：控制網(wǎng)點為分布顯示平差計算成果平差計算成果總輸出點位誤差橢圓參數(shù)實習(xí)總結(jié)：兩周的實習(xí)歷經(jīng)百般波折磨難，終于結(jié)束了，這兩周的時間里，大家經(jīng)歷的是前所未有的實習(xí)體驗，讓我們不單單收貨了專業(yè)能力，更是給我們的人生上了重大的一課，在課堂上我們都是學(xué)習(xí)的書本上的理論知識，這次我們把理論知識和計算機編程結(jié)合在一起，這才是真正學(xué)習(xí)到了未來生活工作中所實用的東西，提高了我

20、們真正面對測量平差工作還有matlab軟件時的操作能力。我們在實習(xí)中發(fā)現(xiàn)自己的問題并及時的解決，讓自己在以后的學(xué)習(xí)工作中不犯同樣錯誤。在剛拿到實習(xí)任務(wù)的第一天我們就進行了分組，分完組后，為了能夠更快更有效率的完成此次實習(xí)，也為了能讓組內(nèi)的每個成員參與進來，首先，我們在組內(nèi)對各自的任務(wù)進行了分配。大家各自負責(zé)一部分的程序代碼編寫內(nèi)容，各自針對自己所主要面對的內(nèi)容查書并向其他各組詢問討論，然后根據(jù)書上的代碼提示，調(diào)通自己負責(zé)的這段程序代碼，每個人把自己的部分完成后，大家在一起編程查錯，經(jīng)過反復(fù)修改弄出整個的成果輸出，還有簡單的界面設(shè)計。我們的組長是王子強，他和王瑞負責(zé)的是主要的數(shù)據(jù)計算部分還有整體

21、修正還有精度評定部分，這幾點可以說是這次程序代碼中的核心所以需要他們兩個人分擔(dān)重任，寧子銘負責(zé)的是誤差橢圓參數(shù)的計算和圖形的繪制內(nèi)容，繆云龍負責(zé)的則是數(shù)據(jù)文件的輸入部分的程序編寫，閃永川負責(zé)的是最后的成果輸出部分，王夢琦則是負責(zé)最后的界面部分程序代碼的設(shè)計，這些分工只是各自主要負責(zé)的項目，但是因為大家能力都有限，所以出現(xiàn)問題后還是集合在一起解決問題再繼續(xù)進行之后的實習(xí)。一開始，先是將已知點信息錄入，這塊代碼是由我們組的繆云龍主要負責(zé)的，這部分的代碼看似簡單，但還是給我們帶來了不小的麻煩，需要在matlab中單獨建立一個以chkdat命名的function文件。這些看似簡單的問題卻也給我們帶來過

22、不少的煩惱。接下來就是由我們組長王子強和王瑞同學(xué)負責(zé)的近似坐標(biāo)計算和相關(guān)精度評定的代碼編寫，這部分非常的復(fù)雜，它包含了方位角計算，編寫誤差方式的程序代碼等或大或小的問題。遇到的首要問題就是我們看不懂書上的近似坐標(biāo)計算的代碼，而書上的代碼又存在問題，缺少nk的定義，而想要修正這一點需要讀懂整段代碼并將nk通過書寫給予定義。還有就是誤差方程的系數(shù)矩陣的編寫程序。這里有好幾個因素要考慮：方向觀測邊的方位角計算，未知數(shù)的系數(shù)位置安放問題，還有就是點號是否是待定點的判斷等等諸多問題。在不斷地學(xué)習(xí)中向其他組成員請教后我們總算輸入了正確的代碼修正了這部分工作。算完上面這些計算之后，就是由寧子銘負責(zé)的誤差橢圓參數(shù)的計算和圖形的繪制了，這部分的代碼編與之前的計算相比就顯得簡單了一點，但由于我們只能參考教材，而