測量平差程序設(shè)計(jì)-邊角三角網(wǎng)平差程序設(shè)計(jì)

*測量平差程序設(shè)計(jì)課程設(shè)計(jì)任務(wù)書設(shè)計(jì)題目邊角三角網(wǎng)平差程序設(shè)計(jì)設(shè)計(jì)資料:BjBjBjBj邊長觀測值固定誤差(mm)1//*終點(diǎn)邊長(m)Bj832.388Bj770.170Bj685.482Bj041516.951Bj04418.530Bj030Bj03Bj04276.43169Bj03Bj010Bj04Bj0415.48098Bj04Bj17.59294Bj04 11終點(diǎn)起點(diǎn)課程設(shè)計(jì)的任務(wù)和內(nèi)容tlab差程序,要求有簡易的界面，數(shù)據(jù)輸入采用文本輸入，采用間接平差模型完成平/*設(shè)計(jì)階段；制功能設(shè)計(jì)階段設(shè)計(jì)程序要實(shí)現(xiàn)的功能4.流程和界面設(shè)計(jì)階段根據(jù)平差計(jì)算的過程和程序功能，畫出流程圖，設(shè)計(jì)簡易界面實(shí)現(xiàn)數(shù)據(jù)的輸入界面。計(jì)報(bào)告撰寫階段目的意義、/*模型測值數(shù)據(jù)組織成數(shù)據(jù)讀入模塊構(gòu)成模塊調(diào)試調(diào)試調(diào)試綜合考核完成平差計(jì)算模塊。北京建筑工程學(xué)院測量平差計(jì)算自編資料2011/*于學(xué)過的測量平差理論的基礎(chǔ)知識(shí)，在matlab中編寫一個(gè)相對(duì)完整的數(shù)據(jù)平差處理程序，旨在讓同學(xué)們能夠?qū)W(xué)到的平差理論知識(shí)和計(jì)算機(jī)的MATLAB軟件程序編寫此提高每個(gè)人學(xué)生的思考和動(dòng)手能力。處理的matlab的程序代碼編寫與試調(diào)任務(wù)，其中主要包括：數(shù)據(jù)組織與輸入、小組成員各上交一份個(gè)人總結(jié)(設(shè)計(jì)書及總結(jié))。模型：平面控制網(wǎng)中的邊角三角網(wǎng)本次的平差數(shù)據(jù)處理是采用間接平差的方法(數(shù)學(xué)模型)而進(jìn)行的，就三個(gè)待而實(shí)際應(yīng)用是采用所提供的平面控制網(wǎng)中的邊角三角網(wǎng)的平差數(shù)據(jù)模型，分別編寫。通過手動(dòng)地選擇數(shù)據(jù)文件，而后自動(dòng)地對(duì)數(shù)據(jù)文件中涉及的變量進(jìn)行賦值和所獲得的數(shù)據(jù)的數(shù)和該橢圓的自動(dòng)繪制，最后平差成果的輸出。入和變量賦值ymd度ma,mb具體代碼截圖：坐標(biāo)的計(jì)算通過所提供的已知點(diǎn)坐標(biāo)計(jì)算得出的已知方位角，在就方向觀測值，求出其每一次觀測的夾其他待定點(diǎn)所在觀測方向(必須有一個(gè)點(diǎn)是已知點(diǎn))的關(guān)系，求出其他待定點(diǎn)所在觀測方向的方位角(注意方位角是否合理問題)，而后在邊長觀測值中找到對(duì)應(yīng)的邊長，利用極坐標(biāo)法最終求得待定點(diǎn)的近似坐標(biāo)。x0y0.同時(shí)還可以將近似坐標(biāo)以文本的文件輸出，計(jì)算代碼截圖如下：本次平面控制網(wǎng)平差模型中觀測量為方向觀測量和邊長觀測量，未知數(shù)分為待定點(diǎn)坐標(biāo)和定點(diǎn)角兩類未知數(shù)，邊角網(wǎng)誤差方程的形式如下：/*中誤差方程的編寫要注意兩個(gè)問題：系數(shù)的位置確定和是否待定點(diǎn)的判斷。對(duì)于方向觀測值，未知數(shù)不僅包括各待定點(diǎn)坐標(biāo)(6個(gè))，同樣包括各測站的定向角，在未知數(shù)排序中，定向角未知數(shù)在前，待定點(diǎn)坐標(biāo)在后，則未知數(shù)個(gè)數(shù)為nz+2*nd=11，對(duì)于方向觀測值，每一列與定向角對(duì)應(yīng)的系數(shù)均為-1，與坐標(biāo)對(duì)應(yīng)的系數(shù)即每一列有5個(gè)系數(shù)(最多)，其他系數(shù)均為0，及對(duì)應(yīng)的系數(shù)矩陣的數(shù)值的確定(注意正負(fù)號(hào))。而對(duì)于邊長觀測值，每一列包括邊的兩個(gè)點(diǎn)的對(duì)應(yīng)的改正數(shù)(如果兩個(gè)點(diǎn)都是待定點(diǎn)的話)四個(gè)參數(shù)，其他均為0，同樣還要注意系數(shù)的位置問題。由于這里分別對(duì)方向觀測量和邊長觀測量設(shè)置兩個(gè)數(shù)組，分別存儲(chǔ)系數(shù)陣常數(shù)項(xiàng)變量，這里設(shè)置b1b2c1c2四個(gè)數(shù)組。根據(jù)觀測值序列和未知數(shù)序列確定系數(shù)矩陣中各元素的位置和大?。涣?，方向觀測量，等輸出量為誤差方程系數(shù)，常數(shù)項(xiàng)。代碼編寫截圖：和精度評(píng)定在得到誤差方程的系數(shù)之后，根據(jù)方向觀測值和邊長觀測值的初步計(jì)算觀測值權(quán)陣(p1p2)，嚴(yán)格的精確計(jì)算可以采用方差分量估計(jì)方法精確計(jì)算方向和邊長觀測值的權(quán)。在得到誤差方程之后，可以直接利用矩陣進(jìn)行方程組的解算(最小二乘法)，誤差方程分為方向觀測誤差方程，邊長觀測誤差方程，兩者可分別計(jì)算，然后再相加。由于近似坐標(biāo)不精確，在求取參數(shù)程序截圖：差橢圓參數(shù)計(jì)算參考書上的誤差橢圓參數(shù)計(jì)算模型由于誤差橢圓描述了點(diǎn)位精度及其在各個(gè)方向的分布，表示了網(wǎng)中點(diǎn)或點(diǎn)之間的誤差分布情況，常用誤差橢圓對(duì)布網(wǎng)方案做精度分析，誤差橢圓三個(gè)參數(shù)為橢圓的長半軸E，短半軸F和主軸方向代碼截圖：7.控制網(wǎng)圖形繪制和誤差橢圓繪制觀測值，方向觀測值連線，構(gòu)成完整的控制網(wǎng)圖。至于誤差橢圓的繪制則是在控制網(wǎng)圖的基礎(chǔ)上再加上待定點(diǎn)對(duì)應(yīng)的誤差橢圓的繪制就完成其中有一個(gè)特別注意的地方就是參數(shù)的設(shè)置，要根據(jù)實(shí)際情況而定。繪制代碼截圖：8.平面控制網(wǎng)成果輸出以生成txt格式文檔的形式進(jìn)行成果輸出。其內(nèi)容主要包括五個(gè)方面內(nèi)容：平差計(jì)算后的坐標(biāo)值，待定點(diǎn)點(diǎn)位中誤差，誤差橢圓參數(shù)計(jì)算等。/*2.坐標(biāo)平差成果輸出3.點(diǎn)位誤差橢圓參數(shù)4.誤差橢圓的繪制5.方向觀測平差成果：方向觀測值，平差值，改正數(shù)；6.邊長觀測平差成果：邊長觀測值，邊長平差值，改正數(shù)；代碼截圖(其中的幾個(gè))：界面的設(shè)計(jì)平差程序)的基礎(chǔ)上增加了另一個(gè)界面，就是小組信息的界面(包括有小組成員照)，而兩個(gè)界面之間利用一個(gè)總的界面(小組平差程序設(shè)計(jì))進(jìn)行選擇性打開主界面—小組平差程序設(shè)計(jì)代碼書寫與調(diào)試：由于上面已經(jīng)展示過幾乎所有的代碼(根據(jù)每一個(gè)模塊)，真正運(yùn)行整一個(gè)程序的時(shí)候卻是程序。一下列舉了一些出現(xiàn)的主要問題。方位角現(xiàn)方位角的運(yùn)算(如待定點(diǎn)的近似坐標(biāo)的求解過程中出現(xiàn)過)，這就應(yīng)該驚醒一個(gè)重要的細(xì)節(jié)：就是所求的方位角是否合理？怎么樣去判斷呢？就是通過if(elseif)語句去判斷其用這個(gè)語句去判斷，這樣才能保證精密。和提高正確性。單位問題單位問題出現(xiàn)在兩處，就是觀測值的單位權(quán)的求解處，二是誤差方程的系數(shù)求解處。由于這兩者是在后面的法方程求解時(shí)都要用上所以，在此之前要保證兩者的單位相對(duì)統(tǒng)一，這個(gè)統(tǒng)一要分兩處，邊長觀測量方程系數(shù)和方向觀測方程系數(shù)，特別注意的是原始數(shù)據(jù)提供的觀測精度(兩種)的單位，和誤差方程的系數(shù)的單位，而在實(shí)際中，我們組應(yīng)用的是方向定，誤差橢圓計(jì)算等，也要注意單位問題。例參數(shù)的調(diào)整才能更好地在顯示誤差橢圓?；《绒D(zhuǎn)角度緊要的差異，就是弧度轉(zhuǎn)角度的小問題，我們原來采用的是MATLAB系統(tǒng)自帶的弧度轉(zhuǎn)角度3的x,y坐標(biāo)上面，差了好幾米(跟提供的平差結(jié)果數(shù)據(jù)相比)，這樣后面的就影響到誤差/*數(shù)程序之后，我們的數(shù)據(jù)就非常正常的出來的。平差程序設(shè)計(jì)成果展示進(jìn)入界面導(dǎo)入數(shù)據(jù)文件數(shù)據(jù)文件的查詢近似坐標(biāo)的計(jì)算控制網(wǎng)點(diǎn)為分布顯示平差計(jì)算成果平差計(jì)算成果總輸出點(diǎn)位誤差橢圓參數(shù)我們都是學(xué)習(xí)的書本上的理論知識(shí)，這次我們把理論知識(shí)和計(jì)算機(jī)編程結(jié)合在一起，這才是的任務(wù)進(jìn)行了分配。大家各自負(fù)責(zé)一部分的程序代碼編寫內(nèi)容，各自針對(duì)自己所主要面對(duì)的我們的組長是王子強(qiáng)，他和王瑞負(fù)責(zé)的是主要的數(shù)據(jù)計(jì)算部分還有整體修正還有精度評(píng)定部分，這幾點(diǎn)可以說是這次程序代碼中的核心所以需要他們兩個(gè)人分擔(dān)重任，寧子銘負(fù)責(zé)的是誤差橢圓參數(shù)的計(jì)算和圖形的繪制內(nèi)容，繆云龍負(fù)責(zé)的則是數(shù)據(jù)文件的輸入部分的程序編寫，閃永川負(fù)責(zé)的是最后的成果輸出部分，王夢琦則是負(fù)責(zé)最后的界面部分程序代碼的設(shè)合在一起解決問題再繼續(xù)進(jìn)行之后的實(shí)習(xí)。似簡單，但還是給我們帶來了不小的麻煩，需要在matlab中單獨(dú)建立一個(gè)以‘chkdat’命名的‘function’文件。這些看似簡單的問題卻也給我們帶來過不少的煩惱。接下來就是由我們組長王子強(qiáng)和王瑞同學(xué)負(fù)責(zé)的近似坐標(biāo)計(jì)算和相關(guān)精度評(píng)定的代碼缺少nk的定義，而想要修正這一點(diǎn)需要讀懂整段代碼并將nk通過書寫給予定義。位角計(jì)算，未知數(shù)的系數(shù)位置安放問題，還有就是點(diǎn)號(hào)是否是待定點(diǎn)的判斷等等諸多問題。在不斷地學(xué)習(xí)中向其他組成員請教后我們總算輸入了正確的代碼修正了這部分工作。/*出圖也完成之后就是由閃永川負(fù)責(zé)的成果的輸出，這部分不難，我們很順利的完成了，當(dāng)然這也源自閃永川之前的準(zhǔn)備比較充足，讓全組的效率得到了提高。最后的就是，由王夢琦負(fù)責(zé)的簡單界面的設(shè)計(jì)，在這部分我們也遇到了不小的問題，callback后面的屬性值應(yīng)該是txt文件的名稱，但一開始我們卻做錯(cuò)了，后來在請教了別組同學(xué)還有老師之后才知道應(yīng)該選用m文件的名稱，后來修改之后基本解決了問題。練，這次的實(shí)習(xí)對(duì)我們每個(gè)人而言都是難忘的，我們不僅收貨