基于matlab的導(dǎo)線網(wǎng)坐標(biāo)計算畢業(yè)設(shè)計

東華理工大學(xué)長江學(xué)院畢業(yè)設(shè)計論文題目：基于Matlab的導(dǎo)線網(wǎng)坐標(biāo)計算EnglishTitle：TraverseNetworkCoordinateCalculationBasedOnMatlab學(xué)生姓名：申請學(xué)位門類：工學(xué)學(xué)士專業(yè)：測繪工程系別：測繪工程系東華理工大學(xué)長江學(xué)院畢業(yè)設(shè)計摘要摘要導(dǎo)線計算是在所有測量工作中經(jīng)常遇見的問題之一，同時導(dǎo)線計算的方法也有很多種，本文主要是利用簡單易懂的Matlab對附合導(dǎo)線、閉合導(dǎo)線和支導(dǎo)線進行相應(yīng)的平差計算。文章首先介紹了附合導(dǎo)線、閉合導(dǎo)線、支導(dǎo)線基本概念和計算方法，其次利用Matlab計算機編程語言對三種導(dǎo)線的計算進行編程實現(xiàn)；最后通過實例驗證，本文利用Matlab編寫的程序正確，通過輸入邊長和角度，能夠快速的得到各控制點的準(zhǔn)確坐標(biāo)。關(guān)鍵詞：Matlab；導(dǎo)線計算；精度評價；計算機編程東華理工大學(xué)長江學(xué)院畢業(yè)設(shè)計ABSTRACTABSTRACTTraversecalculatedinallmeasurementsoftenmetoneoftheproblems,Alsotherearemanykindsoftraversetheway.ThisarticlemainlyistousesimpleMatlabtoconnectingtraverse、closedtraverse、spurtraversetothecorrespondingadjustmentcalculationThearticlefirstintroducestheconnectingtraverse、closedtraverse、spurtraverseisbasicconceptandcalculationmethod.SecondlyusingMatlabcomputerprogramminglanguagewasrealizedbyprogrammingcalculationofthreetraverse.Atlast,throughexamplevalidation,Inthispaper,usingtheMatlabprogram,rightthroughtheinputvariablelengthandAngle,abletoquicklygettheaccuratecoordinatesofeachcontrolpoint.Keywords：Matlab;TraverseCalculated;PrecisionEvaluation;ComputerProgramming東華理工大學(xué)長江學(xué)院畢業(yè)設(shè)計目錄目錄1.TOC\o"1-3"\h\u8103緒論 緒論導(dǎo)線坐標(biāo)計算是在所有測量工作中經(jīng)常遇見的問題之一，同時導(dǎo)線計算的方法也有很多種，主要包括：Excel辦公軟件、CASIO計算器、VB編程等，每種方法各有優(yōu)缺點，本文主要是利用簡單易懂的Matlab對附和導(dǎo)線、閉合導(dǎo)線和支導(dǎo)線進行相應(yīng)的坐標(biāo)平差計算。文章首先介紹了附合導(dǎo)線、閉合導(dǎo)線、支導(dǎo)線基本概念和計算方法，其次利用Matlab計算機編程語言對三種導(dǎo)線的計算進行編程實現(xiàn)；最后通過實例驗證，本文利用Matlab編寫的程序正確，通過輸入邊長和角度，能夠快速的得到各控制點的準(zhǔn)確坐標(biāo)。1.1導(dǎo)線坐標(biāo)計算常用的方法比較（1）Excel辦公軟件我們知道導(dǎo)線測量的成果計算雖然有多種方法，但與我們習(xí)慣的導(dǎo)線成果計算格式不統(tǒng)一，輸入數(shù)據(jù)不直觀，使用起來不習(xí)慣，很難推廣。針對這種情況，我們按照常用的導(dǎo)線成果計算本的樣式，直接利用excel辦公軟件。EXCEL是常用的辦公軟件，利用其自帶的函數(shù)，和其強大的數(shù)據(jù)處理能力，創(chuàng)建一個EXCEL模板，利用預(yù)先編輯好的公式，輸入外業(yè)測量成果后自動進行內(nèi)業(yè)的成果計算。能夠根據(jù)輸入的外業(yè)測量數(shù)據(jù)自動進行坐標(biāo)成果正算；反過來，又能根據(jù)兩點的坐標(biāo)自動反算出兩點之間的距離和方位角，能自動識別測量路線，計算導(dǎo)線間水平角；也能進行簡單的條件分析，對一些誤輸入能自動識別；還能夠有效的保護公式和表格，防止誤操作的更改，或者蓄意的修改計算結(jié)果。（2）CASIO計算器CASIO計算器是測量人員在野外作業(yè)或礦山井下一線作業(yè)時經(jīng)常使用的便攜式計算器，它與一般計算器所不同的是有個自編程序步。CASIO系列可編程計算器，作為內(nèi)、外業(yè)計算工作的輔助設(shè)備。具有體積小、重量輕、攜帶方便、多行顯示、存貯量大、工作方便等特點，編制好測量程序后，計算器通過程序計算，不需要測量人員進行逐步計算，從而消除了人為輸入誤差。而且計算器在計算時小數(shù)位數(shù)是自身進行取舍，所以它的精度可以得到保證，并且比人工逐步計算要高。在導(dǎo)線控制測量中，將復(fù)雜的平差計算過程編成計算器程序，簡化計算過程減輕測量人員的內(nèi)業(yè)計算量。內(nèi)業(yè)計算時，只需按照提示輸入相應(yīng)觀測值便可得到平差成果，并進行必要的精度評定，簡單快捷比專業(yè)的平差軟件更容易掌握。主要缺點是費時、工作量大、還易出錯。（3）VB編程采用VB編程進行導(dǎo)線計算,優(yōu)點是使用方便、直觀、計算速度快。采用計算器計算某一導(dǎo)線,可能需要一兩天時間，而采用程序計算可能只需幾分鐘到十幾分鐘時間。目前市場上,可以買到大型嚴(yán)密平差計算程序,但對于許多中小型日常的測量計算，還需要技術(shù)人員自己去開發(fā)；例如：碎部測量、支導(dǎo)線測量、面積計算、解析交會等。1.2基于Matlab的導(dǎo)線計算背景本文依據(jù)Matlab對三種導(dǎo)線計算進行編程設(shè)計，主要因為該語言相對簡單和直接，它主要有如下特點：（1）編程效率高Matlab是一種面向科學(xué)與工程計算的高級語言，允許使用數(shù)學(xué)形式的語言編寫程序，且比VB和C等語言更加接近我們書寫計算公式的思維方式，用Matlab編寫程序猶如在演算紙上排列出公式與求解問題。因此，Matlab語言也可通俗地稱為演算紙式科學(xué)算法語言。由于它編寫簡單，所以編程效率高，易學(xué)易懂。（2）用戶使用方便Matlab語言是一種解釋執(zhí)行的語言，它靈活、方便，其調(diào)試程序手段豐富，調(diào)試速度快，需要學(xué)習(xí)時間少。人們用任何一種語言編寫程序一般都要經(jīng)過四個步驟：編輯、編譯、鏈接，以及執(zhí)行和調(diào)試。各個步驟之間是順序關(guān)系，編程的過程就是在它們之間做瀑布型的循環(huán)。Matlab語言與其他語言相比，較好的解決了上述問題，把編輯、編譯、鏈接和執(zhí)行融為一體。它能在同一畫面上進行靈活操作，快速排除輸入程序中的書寫錯誤、語法錯誤以至語義錯誤，從而加快了用戶編寫、修改和調(diào)試程序的速度，可以說在編程和調(diào)試過程中它是一種比VB還要簡單的語言。具體的說，Matlab運行時，如直接在命令行輸入Matlab語句（命令），包括調(diào)M文件的語句，每輸入一條語句，就立即對其進行處理，完成編譯、鏈接和運行的全過程。又如，將Matlab源程序編輯為M文件，由于Matlab磁盤文件也是M文件，所以編輯后的源文件就可以直接運行，而不需要進行編譯和鏈接。在運行M文件時，如果有錯，計算機屏幕上會給出詳細(xì)的出錯信息，用戶經(jīng)修改后再執(zhí)行，直到正確為止。所以可以說，MATLAB語言不僅是一種語言，廣義上講是一種該語言的開發(fā)系統(tǒng)，即語言調(diào)試系統(tǒng)。（3）擴充能力強，交互性好高版本的的Matlab語言有豐富的庫函數(shù)，在進行復(fù)雜的數(shù)序運算時可以直接調(diào)用，而且Matlab的庫函數(shù)同用戶文件在形成上一樣，所以用戶文件也可作為Matlab的庫函數(shù)來調(diào)用。因而，用戶可以根據(jù)自己的需要方便地建立和擴充新的庫函數(shù)，以便提高Matlab的使用效率和擴充它的功能[1]。東華理工大學(xué)長江學(xué)院畢業(yè)設(shè)計導(dǎo)線類型及特點2.導(dǎo)線類型及特點2.1導(dǎo)線的布設(shè)2.1.1導(dǎo)線的布設(shè)形式導(dǎo)線可被布設(shè)為單一導(dǎo)線和導(dǎo)線網(wǎng)。兩條以上導(dǎo)線的匯聚點，稱為導(dǎo)線的結(jié)點。單一導(dǎo)線與導(dǎo)線網(wǎng)的區(qū)別，在于導(dǎo)線網(wǎng)具有結(jié)點，而單一導(dǎo)線則不具有結(jié)點。按照不同情況和要求，單一導(dǎo)線可被布設(shè)為附和導(dǎo)線、閉合導(dǎo)線和支導(dǎo)線。導(dǎo)線網(wǎng)可被布設(shè)為自由導(dǎo)線網(wǎng)和附合導(dǎo)線網(wǎng)。（1）附合導(dǎo)線如下圖1所示，導(dǎo)線起始于一個已知控制點而終止于另一個已知控制點。已知控制點上可以有一條或幾條定向邊與之相連接，也可以沒有定向與之相連接。圖1附合導(dǎo)線（2）閉合導(dǎo)線如圖2所示，由一個已知控制點出發(fā)，最終又回到這一點，形成一個閉合的多邊形。在閉合導(dǎo)線的已知控制點上至少應(yīng)該有一條定向邊與之相連接。由于閉合導(dǎo)線是一種可靠性極差的控制網(wǎng)圖形，在實際測量工作中應(yīng)避免單獨使用。圖2閉合導(dǎo)線（3）支導(dǎo)線如圖3所示，從一個已知控制點出發(fā)，既不附合于另一個已知控制點，也不閉合于原來的起始控制點。由于支導(dǎo)線缺乏檢核條件，故一般只限于地形測量中的圖根導(dǎo)線中采用。圖3支導(dǎo)線2.2導(dǎo)線的觀測導(dǎo)線的觀測包括轉(zhuǎn)折角的觀測和導(dǎo)線邊的觀測以及導(dǎo)線點高程的觀測。2.2.1轉(zhuǎn)折角的觀測轉(zhuǎn)折角的觀測一般采用測回法進行。當(dāng)導(dǎo)線點上應(yīng)觀測的方向數(shù)多于2個時，應(yīng)采用方向觀測法進行。各測回間應(yīng)按規(guī)定進行水平度盤配置。各等級導(dǎo)線測量水平角觀測的技術(shù)要求見表1：表1各等級導(dǎo)線測量水平角觀測的技術(shù)要求在進行國家等級導(dǎo)線轉(zhuǎn)折角觀測時，應(yīng)以奇數(shù)測回和偶數(shù)測回分別觀測導(dǎo)線前進方向的左角和右角；左角和右角分別取中數(shù)后，再計算圓周角閉合差，值對于三、四等導(dǎo)線應(yīng)分別不超過±3.5″和±5.0″。在進行一、二級和三級導(dǎo)線轉(zhuǎn)折角觀測時，一般應(yīng)觀測導(dǎo)線前進方向的左角。對于閉合導(dǎo)線，若按逆時針方向進行觀測，則觀測的導(dǎo)線角既是閉合多邊形的內(nèi)角，又是導(dǎo)線前進方向的左角。對于支導(dǎo)線，應(yīng)分別觀測導(dǎo)線前進方向的左角和右角，以增加檢核條件。當(dāng)觀測短邊之間的轉(zhuǎn)折角時，測站偏心和目標(biāo)偏心對轉(zhuǎn)折角的影響將十分明顯。因此，應(yīng)對所用儀器、覘牌和光學(xué)對中器進行嚴(yán)格檢校，并且要特仔細(xì)進行對中和精確照準(zhǔn)。2.2.2導(dǎo)線邊長觀測導(dǎo)線邊長可采用電磁波測距儀測量，也可采用全站儀在測取導(dǎo)線角的同時測取導(dǎo)線邊的邊長。導(dǎo)線邊長應(yīng)對向觀測，以增加檢核條件。電磁波測距儀測量的通常是斜距，還需觀測豎直角，用以將傾斜距離改變?yōu)樗骄嚯x，必要時還應(yīng)將其歸算到橢球面上和高斯平面上。2.2.3三聯(lián)腳架法導(dǎo)線觀測三聯(lián)腳架法通常使用三個既能安置全站儀又能安置帶有覘牌的基座和腳架，基座應(yīng)有通用的光學(xué)對中器。如圖4所示，將全站儀安置在測站i的基座中，帶有覘牌的反射棱鏡安置在后視點i-1和前視點i+1的基座中，進行導(dǎo)線測量。遷站時，導(dǎo)線點i和i+1的腳架和基座不動，只取下全站儀和帶有覘牌的反射棱鏡，在導(dǎo)線點i+1上安置全站儀，在導(dǎo)線點i的基座上安置帶有覘牌的反射棱鏡，并將導(dǎo)線點i-1上的腳架遷至導(dǎo)線點i+2處并予以安置，這樣直到測完整條導(dǎo)線為止。在觀測者精心安置儀器的情況下，三聯(lián)腳架法可以減弱儀器和目標(biāo)對中誤差對測角和測距的影響，從而提高導(dǎo)線的觀測精度，減少了坐標(biāo)傳遞誤差。i-1ii+1i+2圖4三聯(lián)腳架法導(dǎo)線觀測在城市或工業(yè)區(qū)進行導(dǎo)線測量時，可在夜間進行作業(yè)，以避免白天作業(yè)時行人、車輛的干擾，夜間作業(yè)，空氣穩(wěn)定、儀器振動小，并可避免太陽暴曬，從而可提高觀測成果的精度。2.3導(dǎo)線測量的近似平差計算導(dǎo)線測量的目的是獲得各導(dǎo)線點的平面直角坐標(biāo)計算的起始數(shù)據(jù)是已知點坐標(biāo)、已知坐標(biāo)方位角，觀測數(shù)據(jù)為觀測角值和觀測邊長。通常情況下，導(dǎo)線平差應(yīng)進行嚴(yán)密平差，但對于二級及其以下等級的圖根導(dǎo)線允許對以單一導(dǎo)線、單結(jié)點導(dǎo)線網(wǎng)采用近似平差方法進行計算。導(dǎo)線近似平差的基本思路是將角度誤差和邊長誤差分別進行平差計算，先進行角度閉合差的分配，在此基礎(chǔ)在進行坐標(biāo)閉合差的分配，通過調(diào)整坐標(biāo)閉合差，以達(dá)到角度的剩余誤差和邊長誤差的目的。在進行導(dǎo)線測量平差計算之前，首先要按照規(guī)范要求對外業(yè)觀測成果進行檢查和驗算，確保觀測成果無誤并符合限差要求，然后對邊長進行加常數(shù)改正、乘常數(shù)改正、氣象改正和傾斜改正（改正方法見《數(shù)字測圖原理與方法》第五章），對角度和邊長進行歸心改正（有偏心觀測時），以取消系統(tǒng)誤差的影響[2]。2.3.1附合導(dǎo)線的計算如下圖5為附合導(dǎo)線，其坐標(biāo)平差計算步驟如下：圖5附合導(dǎo)線(1)方位角計算（2-1）判斷或＜0，以便判斷所在相應(yīng)象限，加上或減去180度，推算出每條導(dǎo)線邊的坐標(biāo)方位角。(2)由地面實測點的轉(zhuǎn)角依據(jù)公式計算其角度閉合差用，依次加減，得出最后的方位角，用計算的方位角=閉合差值。(3)將角度閉合差除以測邊數(shù)分配到各觀測角中角度閉合差調(diào)整中，觀測角為左角時反符號平均分配到各觀測角中,觀測角為右角時，則按閉合差同符號分配到測角，如有小數(shù),按長邊少分,短邊多分原則。(4)用調(diào)整后的觀測角計算方位角(5)坐標(biāo)增量計算,（2-2）(6)測站測出各邊的距離，用調(diào)整后的方位角進行計算各點計算坐標(biāo)(7)計算導(dǎo)線點坐標(biāo)閉合差為:=計算終點X坐標(biāo)-設(shè)計終點X坐標(biāo)=計算終點Y坐標(biāo)-設(shè)計終點Y坐標(biāo)將閉合邊長差按路線邊長與總邊長的比例反向分配到坐標(biāo)增量后計算各點坐標(biāo)。2.3.2閉合導(dǎo)線的計算如圖6所示為閉合導(dǎo)線，閉合導(dǎo)線的精度評定與具有兩個連接角的附和導(dǎo)線精度評定相同，可以采用角度閉合差和導(dǎo)線全長相對閉合差來評定，閉合導(dǎo)線的可靠性較差，在實際測量中避免單獨使用。由于角度觀測值存在誤差，使得多邊形內(nèi)角和的計算值不等于其理論值，而產(chǎn)生角度閉合差，即圖6閉合導(dǎo)線（1）角度閉合差計算（2-3），n為實測內(nèi)角個數(shù)應(yīng)進行重測，如果則可進行平差計算。（2）角度閉合差改正（反號平均分配）（2-4）（3）計算改正后的角度（2-5）（4）推算方位角公式為：（2-6）判斷△x>或<0，△y>或<0以便判斷所在相應(yīng)象限，加上或減去180度，推算出每條導(dǎo)線邊的坐標(biāo)方位角。(5)計算坐標(biāo)增量,（2-7）（6）坐標(biāo)增量閉合差,（2-8）坐標(biāo)增量閉合差的分配坐標(biāo)改正數(shù)的計算：（2-9）計算改正后的坐標(biāo)增量：（2-10）計算坐標(biāo)（2-11）2.3.3支導(dǎo)線的計算以下圖7為例，支導(dǎo)線計算步驟如下：圖7支導(dǎo)線（1）坐標(biāo)方位角推算設(shè)直線MA的坐標(biāo)方位角已知為，按方位角推算公式：計算導(dǎo)線邊的坐標(biāo)方位角。（2）坐標(biāo)增量的計算由各邊的坐標(biāo)方位角和邊長，按公式、，計算坐標(biāo)增量。導(dǎo)線點坐標(biāo)計算按公式、，計算各導(dǎo)線點坐標(biāo)[3]。東華理工大學(xué)長江學(xué)院畢業(yè)設(shè)計基于Matlab的導(dǎo)線計算程序設(shè)計3.基于Matlab的導(dǎo)線計算程序設(shè)計3.1計算程序概述3.1.1計算程序特點相對于手工計算，坐標(biāo)計算程序計算的主要特點是計算速度快、精度高、數(shù)據(jù)處理自動化，從而把人從繁重的計算工作中解放出來。從程序設(shè)計的角度看，程序設(shè)計與平差計算相對獨立。在平差手工計算時，我們總是面對需要計算的具體問題，所以其數(shù)據(jù)是特定的，計算過程由人實時控制；在計算機程序計算中，在程序設(shè)計時數(shù)據(jù)是抽象的，必須考慮到實際計算中問題的多樣性，以及數(shù)據(jù)計算過程的自動化，所以在程序設(shè)計時必須考慮需要處理的所有問題的普遍性和規(guī)律性。另外，相對于手工計算，在程序計算時，選擇平差方法的依據(jù)不同。在手工計算時，我們通常希望盡量降低計算工作量。當(dāng)必要觀測數(shù)t大于多余觀測數(shù)r時，我們可以選擇條件平差；當(dāng)必要觀測數(shù)t小于多余觀測數(shù)r時，我們可以選擇間接平差，這樣，可以降低平差計算量。然而，在計算機程序計算時，由于計算機計算的快速高效性，我們不是很關(guān)心計算量的問題，而把主要精力集中于方法實現(xiàn)的現(xiàn)實性方面，也就是要求該方法具有較強的規(guī)律性，便于程序設(shè)計的技術(shù)實現(xiàn)。在坐標(biāo)計算程序設(shè)計中，使用間接平差，對于一般控制網(wǎng)，誤差方程形式統(tǒng)一、規(guī)律性強、便于程序設(shè)計；而使用條件平差，誤差方程形式多樣、規(guī)律性差，不利于程序設(shè)計。所以，在本課程中，我們主要使用間接平差方法進行程序設(shè)計?？傊覀冊谶x擇數(shù)學(xué)模型的時候，一定要考慮算法同計算機程序設(shè)計的特點相統(tǒng)一。3.1.2計算程序的基本要求坐標(biāo)計算程序設(shè)計與其它程序設(shè)計相同，應(yīng)當(dāng)滿足一定的要求：（1）程序邏輯結(jié)構(gòu)簡單，清晰易讀，符合結(jié)構(gòu)化程序設(shè)計要求，便于擴展；（2）運算速度快，占用內(nèi)存小，內(nèi)外存之間的交換不宜過于頻繁；（3）數(shù)學(xué)模型及計算方法正確、先進，計算結(jié)果精度高；（4）適應(yīng)性強，便于移植，充分考慮各種可能形式，盡量滿足不同要求與需要。3.1.3計算程序的設(shè)計步驟(1)結(jié)構(gòu)總體設(shè)計；(2)數(shù)據(jù)結(jié)構(gòu)設(shè)計；(3)確定軟件各組成部分的算法及數(shù)據(jù)組織；(4)選定某種表達(dá)式來描述各種算法；(5)程序編寫；(6)程序調(diào)試。3.1.4坐標(biāo)計算處理過程（1）數(shù)據(jù)輸入；（2）角度閉合差計算及分配；（3）方位角推算；（4）坐標(biāo)增量計算；（5）坐標(biāo)增量閉合差分配；（6）坐標(biāo)計算。3.2附合導(dǎo)線計算程序附合導(dǎo)線包括一個連接角的附合導(dǎo)線、兩個連接角的附合導(dǎo)線和無連接角的附合導(dǎo)線等，依據(jù)本文介紹的附合導(dǎo)線計算方法和基本公式，依據(jù)matlab語言對附合導(dǎo)線進行編程,主要的程序代碼如下（具體程序代碼見附錄1）：%觀測數(shù)據(jù)輸入n=input('請輸入導(dǎo)線邊數(shù)');string1=(char(('請輸入第')',a,('條導(dǎo)線邊長')'))';string2=(char(('請輸入第')',s,('個轉(zhuǎn)折角')'))';%起算數(shù)據(jù)輸入%起算方位角的計算%角度閉合差計算fb=fwjAB-fwjCD+sum(jiao1)-n*180;%各觀測角改正后的角度jiao2=jiao1+(fb/6)*(-1);%真方位角的計算fwjz(k)=fwjAB+jiao2(k);fwjz(k)=fwjz(k-1)+jiao2(k)-180;%坐標(biāo)閉合差計算xz(i)=bian(i)*cosd(fwjz(i));yz(i)=bian(i)*sind(fwjz(i));fx=xA+sum(xz)-xC;fy=yA+sum(yz)-yC;%坐標(biāo)閉合差分配vx(i)=(-1)*fx/sum(bian)*bian(i);vy(i)=(-1)*fy/sum(bian)*bian(i);%坐標(biāo)計算x(1)=xA+xz(1)+vx(1)y(1)=yA+yz(1)+vy(1);x(i)=x(i-1)+xz(i)+vx(i);y(i)=y(i-1)+yz(i)+vy(i);3.3閉合導(dǎo)線計算程序閉合導(dǎo)線就是已知一條邊，測量若干個邊長和夾角后又閉合到已知邊的導(dǎo)線測量方法。通過計算平差后，可計算得到經(jīng)過的未知點的平面坐標(biāo)。依據(jù)本章介紹的閉合導(dǎo)線計算方法和基本公式，依據(jù)matlab語言對閉合導(dǎo)線進行編程，主要的程序代碼如下（具體程序代碼見附錄2）：%觀測數(shù)據(jù)的輸入n=input('請輸入閉合導(dǎo)線的點數(shù)');string1=(char(('請輸入第')',a,('條導(dǎo)線的邊長:')'))';string2=(char(('請輸入第')',b,('個轉(zhuǎn)折角:')'))';%起算數(shù)據(jù)的輸入string3=(char('請輸入起始已知方向的起算方位角'));%角度閉合差的計算fb=sum(jiao2)-(n-2)*180;%真方位角計算fwjz(1)=fwj2+jiao2(1)+v(1)-180;%坐標(biāo)閉合差計算xzq=cosd(fwj2)*bian(1);yzq=sind(fwj2)*bian(1);xz(i)=cosd(fwjz1(i))*bian(i+1);yz(i)=sind(fwjz1(i))*bian(i+1);fx=xzq+sum(xz);fy=yzq+sum(yz);fs=sqrt(fx^2+fy^2);k=fs/sum(bian);%坐標(biāo)計算X(1)=x+xzq+vx(1);Y(1)=y+yzq+vy(1);X(i)=X(i-1)+xz(i-1)+vx(i);Y(i)=Y(i-1)+yz(i-1)+vy(i);3.4支導(dǎo)線計算程序由已知控制點出發(fā)，不附合、不閉合于任何已知點的導(dǎo)線。依據(jù)本章介紹的支導(dǎo)線計算方法和基本公式，依據(jù)matlab語言對支導(dǎo)線進行編程，主要的程序代碼如下（具體程序代碼見附錄3）：%觀測數(shù)據(jù)輸入n=input('pleasethe待定點個數(shù):');s1=char('pleaseinputthe邊長ofa:');s2=char('pleaseinputthe轉(zhuǎn)折角ofa:');jiao2(i)=degree1+minute1/60+second1/3600;%已知數(shù)據(jù)輸入%方位角推算fwj(j)=f0+jiao2(j)-180;fwj(j)=fwj(j-1)+jiao2(j)-180;%坐標(biāo)增量計算xzl(k)=S(k)*cosd(fwj(k));yzl(k)=S(k)*sind(fwj(k));%坐標(biāo)推算x(l)=x0+xzl(l);y(l)=y0+yzl(l);x(l)=x(l-1)+xzl(l);y(l)=y(l-1)+yzl(l);東華理工大學(xué)長江學(xué)院畢業(yè)設(shè)計三種導(dǎo)線坐標(biāo)平差實例4.三種導(dǎo)線計算實例4.1附合導(dǎo)線算例利用現(xiàn)實中獲得的一組附合導(dǎo)線數(shù)據(jù)，輸入到所設(shè)計的matlab平差程序進行實現(xiàn)，如下圖8所示：圖8附合導(dǎo)線如上圖8所示，已知兩控制點平面坐標(biāo)分別為：B（2507.69,1215.63）、C（2166.72,1757.29），其中坐標(biāo)方位角=237°59′30″，各觀測角數(shù)據(jù)如下：∠1=99°01′00″，∠2=167°45′36″，∠3=123°11′24″，∠4=189°20′36″，∠5=179°59′18″，∠C=129°27′24″，其導(dǎo)線各邊長分別為：S12=225.85mS23=139.03m,S34=172.57m,S45=100.07m,S5C=102.48m.依據(jù)本文程序設(shè)計對上述數(shù)據(jù)進行平差操作，輸出如下成果截圖：基于Matlab的附合導(dǎo)線程序輸入界面如下圖9所示：圖9程序設(shè)計界面(2)結(jié)果輸出界面如下圖10所示：圖10結(jié)果輸出界面（3）根據(jù)程序輸出的平差結(jié)果，整理得下表2：表2具有兩個連接角的附合導(dǎo)線計算上表為相應(yīng)實例基于Matlab程序設(shè)計的輸出結(jié)果，誤差均在限差之內(nèi)。4.2閉合導(dǎo)線算例利用一組由全站儀觀測的閉合導(dǎo)線的數(shù)據(jù)，對matlab平差程序進行實現(xiàn)，如下圖11所示：圖11閉合導(dǎo)線已知1點坐標(biāo)為（10.00,100.00），S12=50.340，S23=91.515，S34=36.863，S45=51.860，S51=41.999。α12=0°00′00″。各觀測內(nèi)角分別為：∠1=89°58′53″，∠2=91°9′26″，∠3=89°19′25″，∠4=105°37′40″，∠5=163°53′47″。依據(jù)本文程序設(shè)計對上述數(shù)據(jù)進行平差操作，輸出如下成果：（1）基于Matlab的閉合導(dǎo)線程序輸入界面如圖12所示：圖12程序輸入界面（2）結(jié)果輸出界面如下圖13所示：圖13結(jié)果輸出界面（3）根據(jù)程序輸出的平差結(jié)果，整理得下表3：表3閉合導(dǎo)線的計算上表為相應(yīng)實例基于Matlab程序設(shè)計的輸出結(jié)果，誤差均在限差之內(nèi)。4.3支導(dǎo)線算例利用一組支導(dǎo)線的數(shù)據(jù)對matlab程序進行實現(xiàn)，實例如下圖14：圖14支導(dǎo)線如上圖所示，A為一已知點其坐標(biāo)為（2507.69，1215.63），已知坐標(biāo)方位角為237°59′30″，SA-1=225.85，S1-2=139.03，S2-3=172.57，轉(zhuǎn)折角分別為：099°01′00″，167°45′36″，123°11′24″。依據(jù)本文程序設(shè)計對上述數(shù)據(jù)進行平差操作，輸出如下成果：(1）基于Matlab的支導(dǎo)線程序輸入界面,如下圖15所示：圖15程序輸入界面（3）結(jié)果輸出界面如下圖16所示：圖16結(jié)果輸出界面（3）根據(jù)程序輸出的平差結(jié)果，整理得下表4：表4支導(dǎo)線的計算上表為相應(yīng)實例基于Matlab程序設(shè)計的輸出結(jié)果，誤差均在限差之內(nèi)。東華理工大學(xué)長江學(xué)院畢業(yè)設(shè)計結(jié)論結(jié)論導(dǎo)線坐標(biāo)計算是在所有測量工作中經(jīng)常遇見的問題之一，本文在比較了Excel辦公軟件、CASIO計算器、VB編程等導(dǎo)線計算方法的優(yōu)缺點后，決定利用簡單易懂的Matlab對附和導(dǎo)線、閉合導(dǎo)線和支導(dǎo)線進行相應(yīng)的坐標(biāo)平差計算。文章首先介紹了幾種導(dǎo)線坐標(biāo)計算方法的優(yōu)缺點，也敘述了Matlab的特點，其次依據(jù)《數(shù)字測圖原理與方法》介紹了附合導(dǎo)線、閉合導(dǎo)線、支導(dǎo)線基本概念和計算方法，詳細(xì)描述了導(dǎo)線布設(shè)和導(dǎo)線觀測，導(dǎo)線觀測包括轉(zhuǎn)折角和邊長觀測，也敘述了三聯(lián)腳架法的相應(yīng)內(nèi)容。接下來利用Matlab計算機編程語言對三種導(dǎo)線的計算進行編程實現(xiàn)，主要思路為：數(shù)據(jù)輸入、角度閉合差計算及分配、坐標(biāo)閉合差計算及分配、坐標(biāo)計算；最后通過三種導(dǎo)線實例驗證，本文利用Matlab編寫的程序正確，通過輸入邊長和角度，能夠快速的得到各控制點的準(zhǔn)確坐標(biāo)。致謝四年大學(xué)生活一晃而過，回首走過的歲月，留下了很多美好的回憶，心中倍感充實，當(dāng)我寫完這篇畢業(yè)論文的時候，有一種如釋重負(fù)的感覺，感慨良多。首先誠摯的感謝我的論文指導(dǎo)老師劉波老師。他在忙碌的教學(xué)工作中擠出時間來審查、修改我的論文。還有教過我的所有老師們，你們嚴(yán)謹(jǐn)細(xì)致、一絲不茍的作風(fēng)一直是我工作、學(xué)習(xí)中的榜樣；他們循循善誘的教導(dǎo)和不拘一格的思路給予我無盡的啟迪。感謝四年中陪伴在我身邊的同學(xué)、朋友，感謝你們?yōu)槲姨岢龅挠幸娴慕ㄗh和意見，有了你們的支持、鼓勵和幫助，我才能充實的度過了四年的學(xué)習(xí)生活。本論文是在導(dǎo)師劉波老師的悉心指導(dǎo)下完成的。老師淵博的專業(yè)知識，嚴(yán)謹(jǐn)?shù)闹螌W(xué)態(tài)度，精益求精的工作作風(fēng)，誨人不倦的高尚師德，嚴(yán)以律己、寬以待人的崇高風(fēng)范，樸實無華、平易近人的人格魅力對我影響深遠(yuǎn)。不僅使我樹立了遠(yuǎn)大的學(xué)術(shù)目標(biāo)、掌握了基本的研究方法，還使我明白了對待一件事要有認(rèn)真、負(fù)責(zé)的態(tài)度。本論文從選題到完成，每一步都是在指導(dǎo)的悉心指導(dǎo)下完成的，傾注了指導(dǎo)老師大量的心血。在此，謹(jǐn)向劉波老師表示崇高的敬意和衷心的感謝！東華理工大學(xué)長江學(xué)院畢業(yè)設(shè)計參考文獻參考文獻[1]姚連壁、周小平編著.《基于MATLAB的控制網(wǎng)平差程序設(shè)計》，上海：同濟大學(xué)出版社，2006年，1-32.[2]著.武漢大學(xué)出版社，2004年，181-190.[3]武漢大學(xué)測繪學(xué)院測量平差學(xué)科組編著.《誤差理論與測量平差基礎(chǔ)》，武漢：武漢大學(xué)出版社，2010年，1-5.[4]高寧、高彩云.MATLAB在測繪領(lǐng)域中的應(yīng)用[J]，平頂山工學(xué)院學(xué)報，2008年01期.[5]魯鐵定、周世健、張立亭、官云蘭.測量教學(xué)中MATLAB軟件的應(yīng)用[J]，地礦測繪;2004年01期.[6]崔利、武文波.測繪領(lǐng)域中MATLAB的應(yīng)用[J]，遼寧工程技術(shù)大學(xué)學(xué)報，2005年S1期.[7]陳以新.Matlab的幾則程序設(shè)計經(jīng)驗[J]，計算機應(yīng)用研究，1999年09期.[8]孫慶華、熊偉、張憲柱、王磊.基于Excel和MATLA導(dǎo)線網(wǎng)平差[J]，北京測繪，2009年03期，1-32.[9]白征東.Matlab在測量平差教學(xué)中的應(yīng)用[J]，測繪通報，2009年11期,1-23.[10]高彩云、高寧、孫敬克.MATLAB在測量平差課程教學(xué)中的嘗試[J]，科技信息，2009年05期，1-12.[11]尹暉主編.《測繪工程專業(yè)英語》，武漢：武漢大學(xué)出版社，2005年1月第1版，1-43.[12]Ingenieurvermessung2004,InternationalcourseonEngineeringSurveying[C].ETHZuerich,2004.[13]KahmenH.HybridMeasurementRobotsinEngineeringSurveys[J].ICongres,Washington,東華理工大學(xué)長江學(xué)院畢業(yè)設(shè)計附錄附錄1基于Matlab的附合導(dǎo)線平差程序：%觀測數(shù)據(jù)輸入n=input('請輸入導(dǎo)線邊數(shù)');fori=1:na=num2str(i);string1=(char(('請輸入第')',a,('條導(dǎo)線邊長')'))';bian(i)=input(string1);endfork=1:(1+n)s=num2str(k);string2=(char(('請輸入第')',s,('個轉(zhuǎn)折角')'))';jiao=input(string2,'s');B=jiao;du=str2num(B(1:3));fen=str2num(B(4:5));miao=str2num(B(6:7));jiao1(k)=du+fen/60+miao/3600;end%起算數(shù)據(jù)輸入xA=input('請輸入已知點A的x坐標(biāo)');yA=input('請輸入已知點A的y坐標(biāo)');xB=input('請輸入已知點B的x坐標(biāo)');yB=input('請輸入已知點B的y坐標(biāo)');%起算方位角的計算Rab=atand((abs(yB-yA))/(abs(xB-xA)));%AB的起算方位角if(yB-yA)>0&(xB-xA)>0fwjAB=Rab;endif(yB-yA)>0&(xB-xA)<0fwjAB=180-Rab;endif(yB-yA)<0&(xB-xA)<0fwjAB=180+Rab;endif(yB-yA)<0&(xB-xA)>0fwjAB=360-Rab;endRcd=atand((abs(yD-yC))/(abs(xD-xC)));%CD的起算方位角if(yD-yC)>0&(xD-xC)>0fwjCD=Rcd;endif(yD-yC)>0&(xD-xC)<0fwjCD=180-Rcd;endif(yD-yC)<0&(xD-xC)<0fwjCD=180+Rcd;endif(yD-yC)<0&(xD-xC)>0fwjCD=360-Rcd;end%角度閉合差計算fb=fwjAB-fwjCD+sum(jiao1)-n*180;jiao2=jiao1+(fb/6)*(-1);%各觀測角改正后的角度%真方位角的計算fork=1:(n+1)ifk==1fwjz(k)=fwjAB+jiao2(k);endifk>1fwjz(k)=fwjz(k-1)+jiao2(k)-180;endend%坐標(biāo)閉合差計算fori=1:nxz(i)=bian(i)*cosd(fwjz(i));fori=1:nyz(i)=bian(i)*sind(fwjz(i));endendfx=xA+sum(xz)-xC;fy=yA+sum(yz)-yC;%坐標(biāo)閉合差分配fori=1:nvx(i)=(-1)*fx/sum(bian)*bian(i);endfori=1:nvy(i)=(-1)*fy/sum(bian)*bian(i);end%坐標(biāo)計算fori=1:nifi==1x(1)=xA+xz(1)+vx(1)y(1)=yA+yz(1)+vy(1);endifi>1x(i)=x(i-1)+xz(i)+vx(i);y(i)=y(i-1)+yz(i)+vy(i);end附錄2基于Matlab的閉合導(dǎo)線平差程序：%觀測數(shù)據(jù)的輸入n=input('請輸入閉合導(dǎo)線的點數(shù)');fori=1:na=num2str(i);string1=(char(('請輸入第')',a,('條導(dǎo)線的邊長:')'))';bian(i)=input(string1);b=num2str(i);string2=(char(('請輸入第')',b,('個轉(zhuǎn)折角:')'))';jiao1=input(string2,'s');c=jiao1;du=str2num(c(1:3));fen=str2num(c(4:5));miao=str2num(c(6:7));jiao2(i)=du+fen/60+miao/3600;end%起算數(shù)據(jù)的輸入x=input('請輸入起始已知點x的坐標(biāo):');y=input('請輸入起始已知點y的坐標(biāo):');string3=(char('請輸入起始已知方向的起算方位角'));fwj1=input(string3,'s');c=fwj1;du=str2num(c(1:3));fen=str2num(c(4:5));miao=str2num(c(6:7));fwj2=du+fen/60+miao/3600;%角度閉合差的計算fb=sum(jiao2)-(n-2)*180;fori=1:nv(i)=(-1)*fb/n;end%真方位角計算fori=1:nifi==1fwjz(1)=fwj2+jiao2(1)+v(1)-180;endifi>1fwjz(i)=fwjz(i-1)+jiao2(i)+v(i)-180;endiffwjz(i)>0fwjz1(i)=fwjz(i);endiffwjz(i)<0fwjz1(i)=fwjz(i)+360;end%坐標(biāo)閉合差計算xzq=cosd(fwj2)*bian(1);yzq=sind(fwj2)*bian(1);fori=1:(n-1)xz(i)=cosd(fwjz1(i))*bian(i+1);yz(i)=sind(fwjz1(i))*bian(i+1);endfx=xzq+sum(xz);fy=yzq+sum(yz);fs=sqrt(fx^2+fy^2);k=fs/sum(bian);%坐標(biāo)計算fori=1:nvx=(-1)*fx/sum(bian)*bian;vy=(-1)*fy/sum(bian)*bian;endfori=1:nifi==1X(1)=x+xzq+vx(1);Y(1)=y+yzq+vy(1);endifi>1X(i)=X(i-1)+xz(i-1)+vx(i);Y(i)=Y(i-1)+yz(i-1)+vy(i);end附錄3基于Matlab的支導(dǎo)線平差程序：%觀測數(shù)據(jù)輸入n=input('pleasethe待定點個數(shù):');fori=1:na=num2str(i);s1=char('pleaseinputthe邊長ofa:');S(i)=input(s1);s2=char('pleaseinputthe轉(zhuǎn)折角ofa:');jiao1=input(s2,'s');degree1=str2num(jiao1(1:3));minute1=str2num(jiao1(4:5));second1=str2num(jiao1(6:7));jiao2(i)=degree1+minute1/60+second1/3600;end%已知數(shù)據(jù)輸入x0=input('pleaseinputtheknownpointofx0:');y0=input('pleaseinputtheknownpointofy0:');s3=char('pleaseinputtheknownazimuth:');fwj0=input(s3,'s');degree2=str2num(fwj0(1:3));minute2=str2num(fwj0(4:5));second2=str2num(fwj0(6:7));f0=degree2+minute2/60+second2/3600%方位角推算forj=1:nifj==1fwj(j)=f0+jiao2(j)-180;endifj>1fwj(j)=fwj(j-1)+jiao2(j)-180;endend%坐標(biāo)增量計算fork=1:nxzl(k)=S(k)*cosd(fwj(k));yzl(k)=S(k)*sind(fwj(k));end%坐標(biāo)推算forl=1:nifl==1x(l)=x0+xzl(l);y(l)=y0+yzl(l);endifl>1x(l)=x(l-1)+xzl(l);y(l)=y(l-1)+yzl(l);endend基于C8051F單片機直流電動機反饋控制系統(tǒng)的設(shè)計與研究基于單片機的嵌入式Web服務(wù)器的研究MOTOROLA單片機MC68HC（8）05PV8/A內(nèi)嵌EEPROM的工藝和制程方法及對良率的影響研究基于模糊控制的電阻釬焊單片機溫度控制系統(tǒng)的研制基于MCS-51系列單片機的通用控制模塊的研究基于單片機實現(xiàn)的供暖系統(tǒng)最佳啟停自校正（STR）調(diào)節(jié)器單片機控制的二級倒立擺系統(tǒng)的研究基于增強型51系列單片機的TCP/IP協(xié)議棧的實現(xiàn)基于單片機的蓄電池自動監(jiān)測系統(tǒng)基于32位嵌入式單片機系統(tǒng)的圖像采集與處理技術(shù)的研究基于單片機的作物營養(yǎng)診斷專家系統(tǒng)的研究基于單片機的交流伺服電機運動控制系統(tǒng)研究與開發(fā)基于單片機的泵管內(nèi)壁硬度測試儀的研制基于單片機的自動找平控制系統(tǒng)研究基于C8051F040單片機的嵌入式系統(tǒng)開發(fā)基于單片機的液壓動力系統(tǒng)狀態(tài)監(jiān)測儀開發(fā)模糊Smith智能控制方法的研究及其單片機實現(xiàn)一種基于單片機的軸快流CO〈,2〉激光器的手持控制面板的研制基于雙單片機沖床數(shù)控系統(tǒng)的研究基于CYGNAL單片機的在線間歇式濁度儀的研制基于單片機的噴油泵試驗臺控制器的研制基于單片機的軟起動器的研究和設(shè)計基于單片機控制的高速快走絲電火花線切割機床短循環(huán)走絲方式研究基于單片機的機電產(chǎn)品控制系統(tǒng)開發(fā)基于PIC單片機的智能手機充電器基于單片機的實時內(nèi)核設(shè)計及其應(yīng)用研究基于單片機的遠(yuǎn)程抄表系統(tǒng)的設(shè)計與研究基于單片機的煙氣二氧化硫濃度檢測儀的研制基于微型光譜儀的單片機系統(tǒng)單片機系統(tǒng)軟件構(gòu)件開發(fā)的技術(shù)研究基于單片機的液體點滴速度自動檢測儀的研制基于單片機系統(tǒng)的多功能溫度測量儀的研制基于PIC單片機的電能采集終端的設(shè)計和應(yīng)用基于單片機的光纖光柵解調(diào)儀的研制氣壓式線性摩擦焊機單片機控制系統(tǒng)的研制基于單片機的數(shù)字磁通門傳感器基于單片機的旋轉(zhuǎn)變壓器-數(shù)字轉(zhuǎn)換器的研究基于單片機的光纖Bragg光柵解調(diào)系統(tǒng)的研究單片機控制的便攜式多功能乳腺治療儀的研制基于C8051F020單片機的多生理信號檢測儀基于單片機的電機運動控制系統(tǒng)設(shè)計Pico專用單片機核的可測性設(shè)計研究基于MCS-51單片機的熱量計基于雙單片機的智能遙測微型氣象站MCS-51單片機構(gòu)建機器人的實踐研究基于單片機的輪軌力檢測基于單片機的GPS定位儀的研究與實現(xiàn)基于單片機的電液伺服控制系統(tǒng)用于單片機系統(tǒng)的MMC卡文件系統(tǒng)研制基于單片機的時控和計數(shù)系統(tǒng)性能優(yōu)化的研究基于單片機和CPLD的粗光柵位移測量系統(tǒng)研究單片機控制的后備式方波UPS提升高職學(xué)生單片機應(yīng)用能力的探究基于單片機控制的自動低頻減載裝置研究基于單片機控制的水下焊接電源的研究基于單片機的多通道數(shù)據(jù)采集系統(tǒng)基于uPSD3234單片機的氚表面污染測量儀的研制基于單片機的紅外測油儀的研究96系列單片機仿真器研究與設(shè)計基于單片機的單晶金剛石刀具刃磨設(shè)備的數(shù)控改造基于單片機的溫度智能控制系統(tǒng)的設(shè)計與實現(xiàn)基于MSP430單片機的電梯門機控制器的研制基于單片機的氣體測漏儀的研究基于三菱M16C/6N系列單片機的CAN/USB協(xié)議轉(zhuǎn)換器基于單片機和DSP的變壓器油色譜在線監(jiān)測技術(shù)研究基于單片機的膛壁溫度報警系統(tǒng)設(shè)計基于AVR單片機的低壓無功補償控制器的設(shè)計基于單片機船舶電力推進電機監(jiān)測系統(tǒng)基于單片機網(wǎng)絡(luò)的振動信號的采集系統(tǒng)基于單片機的大容量數(shù)據(jù)存儲技術(shù)的應(yīng)用研究基于單片機的疊圖機研究與教學(xué)方法實踐基于單片機嵌入式Web服務(wù)器技術(shù)的研究及實現(xiàn)基于AT89S52單片機的通用數(shù)據(jù)采集系統(tǒng)基于單片機的多道脈沖幅度分析儀研究機器人旋轉(zhuǎn)電弧傳感角焊縫跟蹤單片機控制系統(tǒng)基于單片機的控制系統(tǒng)在PLC虛擬教學(xué)實驗中的應(yīng)用研究基于單片機系統(tǒng)的網(wǎng)絡(luò)通信研究與應(yīng)用基于PIC16F877單片機的莫爾斯碼自動譯碼系統(tǒng)設(shè)計與研究基于單片機的模糊控制器在工業(yè)電阻爐上的應(yīng)用研究HYPERLINK"/det