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

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