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

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

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

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

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

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

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

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

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

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

10、坐標(biāo)計(jì)算得出的已知方位角，在就方向觀測值，求出其每一次觀測的夾角，而后找出每個(gè)測站所擁有的已知點(diǎn)（一個(gè)觀測方向的兩個(gè)點(diǎn)）確定其所在觀測的行數(shù)與其他待定點(diǎn)所在觀測方向（必須有一個(gè)點(diǎn)是已知點(diǎn)）的關(guān)系，求出其他待定點(diǎn)所在觀測方向的方位角（注意方位角是否合理問題），而后在邊長觀測值中找到對應(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)誤差方程的形式如下：方向觀測量：誤差方程：（系數(shù)為：b1,c1）式中常數(shù)項(xiàng)：變長觀測量

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

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

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

14、繪制控制網(wǎng)點(diǎn)位分布圖，點(diǎn)位標(biāo)注，邊長觀測量和方向觀測量的繪制等。基本思路是首先根據(jù)坐標(biāo)繪制控制網(wǎng)散點(diǎn)圖，并進(jìn)行標(biāo)注，然后根據(jù)距離觀測值，方向觀測值連線，構(gòu)成完整的控制網(wǎng)圖。至于誤差橢圓的繪制則是在控制網(wǎng)圖的基礎(chǔ)上再加上待定點(diǎn)對應(yīng)的誤差橢圓的繪制就完成了。（利用Plot命令）其中有一個(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ì)算等。信息具體內(nèi)容如下：1.平差計(jì)算成果成果總輸出；2.坐標(biāo)平差成果輸出3.點(diǎn)位誤差橢圓參數(shù)4.誤差橢圓的繪制

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

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

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

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

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

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

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

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