基于Matlab的導(dǎo)線網(wǎng)坐標(biāo)計(jì)算設(shè)計(jì)說明

1、 PAGE33 / NUMPAGES38東華理工大學(xué)長江學(xué)院畢業(yè)設(shè)計(jì) 論文題目：基于Matlab的導(dǎo)線網(wǎng)坐標(biāo)計(jì)算English Title：Traverse Network Coordinate CalculationBased On Matlab學(xué)生：閆曉天 申請學(xué)位門類：工學(xué)學(xué)士 學(xué)號：09351122專業(yè)：測繪工程 系別：測繪工程系指導(dǎo)老師：波職稱：講師二 零 一 三 年 六 月摘 要導(dǎo)線計(jì)算是在所有測量工作中經(jīng)常遇見的問題之一，同時導(dǎo)線計(jì)算的方法也有很多種，本文主要是利用簡單易懂的Matlab對附合導(dǎo)線、閉合導(dǎo)線和支導(dǎo)線進(jìn)行相應(yīng)的平差計(jì)算。文章首先介紹了附合導(dǎo)線、閉合導(dǎo)線、支導(dǎo)線基本概

2、念和計(jì)算方法，其次利用Matlab計(jì)算機(jī)編程語言對三種導(dǎo)線的計(jì)算進(jìn)行編程實(shí)現(xiàn)；最后通過實(shí)例驗(yàn)證，本文利用Matlab編寫的程序正確，通過輸入邊長和角度，能夠快速的得到各控制點(diǎn)的準(zhǔn)確坐標(biāo)。關(guān)鍵詞：Matlab；導(dǎo)線計(jì)算；精度評價；計(jì)算機(jī)編程ABSTRACTTraverse calculated in all measurements often met one of the problems,Also there are many kinds of traverse the way.This article mainly is to use simple Matlab to connecting

3、 traverse、closed traverse、spur traverse to the corresponding adjustment calculation The article first introduces the connecting traverse、closed traverse、spur traverse is basic concept and calculation method.Secondly using Matlab computer programming language was realized by programming calculation o

4、f three traverse.At last,through example validation, In this paper, using the Matlab program, right through the input variable length and Angle, able to quickly get the accurate coordinates of each control point.Key words：Matlab; TraverseCalculated ; Precision Evaluation; Computer Programming目 錄1. T

5、OC o 1-3 h u HYPERLINK l _Toc8103 緒 論 PAGEREF _Toc8103 1 HYPERLINK l _Toc22805 1.1 導(dǎo)線坐標(biāo)計(jì)算常用的方法比較 PAGEREF _Toc22805 1 HYPERLINK l _Toc23673 1.2 基于Matlab的導(dǎo)線計(jì)算背景 PAGEREF _Toc23673 2 HYPERLINK l _Toc10806 2. 導(dǎo)線類型與特點(diǎn) PAGEREF _Toc10806 3 HYPERLINK l _Toc17523 2.1 導(dǎo)線的布設(shè) PAGEREF _Toc17523 3 HYPERLINK l _Toc

6、2509 2.1.1 導(dǎo)線的布設(shè)形式 PAGEREF _Toc2509 3 HYPERLINK l _Toc25147 2.2 導(dǎo)線的觀測 PAGEREF _Toc25147 4 HYPERLINK l _Toc29039 2.2.1 轉(zhuǎn)折角的觀測 PAGEREF _Toc29039 4 HYPERLINK l _Toc1164 2.2.2 導(dǎo)線邊長觀測 PAGEREF _Toc1164 5 HYPERLINK l _Toc10873 2.2.3 三聯(lián)腳架法導(dǎo)線觀測 PAGEREF _Toc10873 5 HYPERLINK l _Toc30360 2.3 導(dǎo)線測量的近似平差計(jì)算 PAGERE

7、F _Toc30360 6 HYPERLINK l _Toc1398 2.3.1 附合導(dǎo)線的計(jì)算 PAGEREF _Toc1398 6 HYPERLINK l _Toc14424 2.3.2 閉合導(dǎo)線的計(jì)算 PAGEREF _Toc14424 7 HYPERLINK l _Toc5555 2.3.3 支導(dǎo)線的計(jì)算 PAGEREF _Toc5555 8 HYPERLINK l _Toc7758 3. 基于Matlab的導(dǎo)線計(jì)算程序設(shè)計(jì) PAGEREF _Toc7758 10 HYPERLINK l _Toc23142 3.1 計(jì)算程序概述 PAGEREF _Toc23142 10 HYPERLI

8、NK l _Toc15510 3.1.1 計(jì)算程序特點(diǎn) PAGEREF _Toc15510 10 HYPERLINK l _Toc18960 3.1.2 計(jì)算程序的基本要求 PAGEREF _Toc18960 10 HYPERLINK l _Toc24370 3.1.3 計(jì)算程序的設(shè)計(jì)步驟 PAGEREF _Toc24370 10 HYPERLINK l _Toc6667 3.1.4 坐標(biāo)計(jì)算處理過程 PAGEREF _Toc6667 11 HYPERLINK l _Toc8690 3.2 附合導(dǎo)線計(jì)算程序 PAGEREF _Toc8690 11 HYPERLINK l _Toc11101 3

9、.3 閉合導(dǎo)線計(jì)算程序 PAGEREF _Toc11101 12 HYPERLINK l _Toc3178 3.4 支導(dǎo)線計(jì)算程序 PAGEREF _Toc3178 13 HYPERLINK l _Toc3910 4. 三種導(dǎo)線計(jì)算實(shí)例 PAGEREF _Toc3910 14 HYPERLINK l _Toc31446 4.1 附合導(dǎo)線算例 PAGEREF _Toc31446 14 HYPERLINK l _Toc26920 4.2 閉合導(dǎo)線算例 PAGEREF _Toc26920 18 HYPERLINK l _Toc8924 4.3 支導(dǎo)線算例 PAGEREF _Toc8924 21 HY

10、PERLINK l _Toc24755 結(jié) 論 PAGEREF _Toc24755 25 HYPERLINK l _Toc15236 致 PAGEREF _Toc15236 26 HYPERLINK l _Toc11907 參考文獻(xiàn) PAGEREF _Toc11907 27 HYPERLINK l _Toc18918 附錄1 PAGEREF _Toc18918 28 HYPERLINK l _Toc24202 附錄2 PAGEREF _Toc24202 30 HYPERLINK l _Toc24674 附錄3 PAGEREF _Toc24674 32緒 論導(dǎo)線坐標(biāo)計(jì)算是在所有測量工作中經(jīng)常遇見

11、的問題之一，同時導(dǎo)線計(jì)算的方法也有很多種，主要包括：Excel辦公軟件、CASIO計(jì)算器、VB編程等，每種方法各有優(yōu)缺點(diǎn)，本文主要是利用簡單易懂的Matlab對附和導(dǎo)線、閉合導(dǎo)線和支導(dǎo)線進(jìn)行相應(yīng)的坐標(biāo)平差計(jì)算。文章首先介紹了附合導(dǎo)線、閉合導(dǎo)線、支導(dǎo)線基本概念和計(jì)算方法，其次利用Matlab計(jì)算機(jī)編程語言對三種導(dǎo)線的計(jì)算進(jìn)行編程實(shí)現(xiàn)；最后通過實(shí)例驗(yàn)證，本文利用Matlab編寫的程序正確，通過輸入邊長和角度，能夠快速的得到各控制點(diǎn)的準(zhǔn)確坐標(biāo)。1.1 導(dǎo)線坐標(biāo)計(jì)算常用的方法比較（1）Excel辦公軟件我們知道導(dǎo)線測量的成果計(jì)算雖然有多種方法，但與我們習(xí)慣的導(dǎo)線成果計(jì)算格式不統(tǒng)一，輸入數(shù)據(jù)不直觀，使用

12、起來不習(xí)慣，很難推廣。針對這種情況，我們按照常用的導(dǎo)線成果計(jì)算本的樣式，直接利用excel辦公軟件。EXCEL是常用的辦公軟件，利用其自帶的函數(shù)，和其強(qiáng)大的數(shù)據(jù)處理能力，創(chuàng)建一個EXCEL模板，利用預(yù)先編輯好的公式，輸入外業(yè)測量成果后自動進(jìn)行業(yè)的成果計(jì)算。能夠根據(jù)輸入的外業(yè)測量數(shù)據(jù)自動進(jìn)行坐標(biāo)成果正算；反過來，又能根據(jù)兩點(diǎn)的坐標(biāo)自動反算出兩點(diǎn)之間的距離和方位角，能自動識別測量路線，計(jì)算導(dǎo)線間水平角；也能進(jìn)行簡單的條件分析，對一些誤輸入能自動識別；還能夠有效的保護(hù)公式和表格，防止誤操作的更改，或者蓄意的修改計(jì)算結(jié)果。（2）CASIO計(jì)算器CASIO計(jì)算器是測量人員在野外作業(yè)或礦山井下一線作業(yè)時經(jīng)

13、常使用的便攜式計(jì)算器，它與一般計(jì)算器所不同的是有個自編程序步。CASIO系列可編程計(jì)算器，作為、外業(yè)計(jì)算工作的輔助設(shè)備。具有體積小、重量輕、攜帶方便、多行顯示、存貯量大、工作方便等特點(diǎn)，編制好測量程序后，計(jì)算器通過程序計(jì)算，不需要測量人員進(jìn)行逐步計(jì)算，從而消除了人為輸入誤差。而且計(jì)算器在計(jì)算時小數(shù)位數(shù)是自身進(jìn)行取舍，所以它的精度可以得到保證，并且比人工逐步計(jì)算要高。在導(dǎo)線控制測量中，將復(fù)雜的平差計(jì)算過程編成計(jì)算器程序，簡化計(jì)算過程減輕測量人員的業(yè)計(jì)算量。業(yè)計(jì)算時，只需按照提示輸入相應(yīng)觀測值便可得到平差成果，并進(jìn)行必要的精度評定，簡單快捷比專業(yè)的平差軟件更容易掌握。主要缺點(diǎn)是費(fèi)時、工作量大、還易

14、出錯。（3）VB編程采用VB編程進(jìn)行導(dǎo)線計(jì)算,優(yōu)點(diǎn)是使用方便、直觀、計(jì)算速度快。采用計(jì)算器計(jì)算某一導(dǎo)線,可能需要一兩天時間，而采用程序計(jì)算可能只需幾分鐘到十幾分鐘時間。目前市場上,可以買到大型嚴(yán)密平差計(jì)算程序,但對于許多中小型日常的測量計(jì)算，還需要技術(shù)人員自己去開發(fā)；例如：碎部測量、支導(dǎo)線測量、面積計(jì)算、解析交會等。1.2 基于Matlab的導(dǎo)線計(jì)算背景本文依據(jù)Matlab對三種導(dǎo)線計(jì)算進(jìn)行編程設(shè)計(jì)，主要因?yàn)樵撜Z言相對簡單和直接，它主要有如下特點(diǎn)：（1）編程效率高M(jìn)atlab是一種面向科學(xué)與工程計(jì)算的高級語言，允許使用數(shù)學(xué)形式的語言編寫程序，且比VB和C等語言更加接近我們書寫計(jì)算公式的思維方式

15、，用Matlab編寫程序猶如在演算紙上排列出公式與求解問題。因此，Matlab語言也可通俗地稱為演算紙式科學(xué)算法語言。由于它編寫簡單，所以編程效率高，易學(xué)易懂。（2）用戶使用方便 Matlab語言是一種解釋執(zhí)行的語言，它靈活、方便，其調(diào)試程序手段豐富，調(diào)試速度快，需要學(xué)習(xí)時間少。人們用任何一種語言編寫程序一般都要經(jīng)過四個步驟：編輯、編譯、，以與執(zhí)行和調(diào)試。各個步驟之間是順序關(guān)系，編程的過程就是在它們之間做瀑布型的循環(huán)。Matlab語言與其他語言相比，較好的解決了上述問題，把編輯、編譯、和執(zhí)行融為一體。它能在同一畫面上進(jìn)行靈活操作，快速排除輸入程序中的書寫錯誤、語法錯誤以至語義錯誤，從而加快了用

16、戶編寫、修改和調(diào)試程序的速度，可以說在編程和調(diào)試過程中它是一種比VB還要簡單的語言。具體的說，Matlab運(yùn)行時，如直接在命令行輸入Matlab語句（命令），包括調(diào)M文件的語句，每輸入一條語句，就立即對其進(jìn)行處理，完成編譯、和運(yùn)行的全過程。又如，將Matlab源程序編輯為M文件，由于Matlab磁盤文件也是M文件，所以編輯后的源文件就可以直接運(yùn)行，而不需要進(jìn)行編譯和。在運(yùn)行M文件時，如果有錯，計(jì)算機(jī)屏幕上會給出詳細(xì)的出錯信息，用戶經(jīng)修改后再執(zhí)行，直到正確為止。所以可以說，MATLAB語言不僅是一種語言，廣義上講是一種該語言的開發(fā)系統(tǒng)，即語言調(diào)試系統(tǒng)。（3）擴(kuò)充能力強(qiáng)，交互性好高版本的的Matl

17、ab語言有豐富的庫函數(shù)，在進(jìn)行復(fù)雜的數(shù)序運(yùn)算時可以直接調(diào)用，而且Matlab的庫函數(shù)同用戶文件在形成上一樣，所以用戶文件也可作為Matlab的庫函數(shù)來調(diào)用。因而，用戶可以根據(jù)自己的需要方便地建立和擴(kuò)充新的庫函數(shù)，以便提高M(jìn)atlab的使用效率和擴(kuò)充它的功能1。2. 導(dǎo)線類型與特點(diǎn)2.1 導(dǎo)線的布設(shè)2.1.1 導(dǎo)線的布設(shè)形式導(dǎo)線可被布設(shè)為單一導(dǎo)線和導(dǎo)線網(wǎng)。兩條以上導(dǎo)線的匯聚點(diǎn)，稱為導(dǎo)線的結(jié)點(diǎn)。單一導(dǎo)線與導(dǎo)線網(wǎng)的區(qū)別，在于導(dǎo)線網(wǎng)具有結(jié)點(diǎn)，而單一導(dǎo)線則不具有結(jié)點(diǎn)。按照不同情況和要求，單一導(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)

18、線起始于一個已知控制點(diǎn)而終止于另一個已知控制點(diǎn)。已知控制點(diǎn)上可以有一條或幾條定向邊與之相連接，也可以沒有定向與之相連接。圖1 附合導(dǎo)線（2） 閉合導(dǎo)線如圖2所示，由一個已知控制點(diǎn)出發(fā)，最終又回到這一點(diǎn)，形成一個閉合的多邊形。在閉合導(dǎo)線的已知控制點(diǎn)上至少應(yīng)該有一條定向邊與之相連接。由于閉合導(dǎo)線是一種可靠性極差的控制網(wǎng)圖形，在實(shí)際測量工作中應(yīng)避免單獨(dú)使用。 圖2 閉合導(dǎo)線（3） 支導(dǎo)線如圖3所示，從一個已知控制點(diǎn)出發(fā)，既不附合于另一個已知控制點(diǎn)，也不閉合于原來的起始控制點(diǎn)。由于支導(dǎo)線缺乏檢核條件，故一般只限于地形測量中的圖根導(dǎo)線中采用。 圖3 支導(dǎo)線2.2 導(dǎo)線的觀測導(dǎo)線的觀測包括轉(zhuǎn)折角的觀測和導(dǎo)

19、線邊的觀測以與導(dǎo)線點(diǎn)高程的觀測。2.2.1 轉(zhuǎn)折角的觀測轉(zhuǎn)折角的觀測一般采用測回法進(jìn)行。當(dāng)導(dǎo)線點(diǎn)上應(yīng)觀測的方向數(shù)多于2個時，應(yīng)采用方向觀測法進(jìn)行。各測回間應(yīng)按規(guī)定進(jìn)行水平度盤配置。各等級導(dǎo)線測量水平角觀測的技術(shù)要求見表1： 表1 各等級導(dǎo)線測量水平角觀測的技術(shù)要求在進(jìn)行國家等級導(dǎo)線轉(zhuǎn)折角觀測時，應(yīng)以奇數(shù)測回和偶數(shù)測回分別觀測導(dǎo)線前進(jìn)方向的左角和右角；左角和右角分別取中數(shù)后，再計(jì)算圓周角閉合差，值對于三、四等導(dǎo)線應(yīng)分別不超過3.5和5.0。在進(jìn)行一、二級和三級導(dǎo)線轉(zhuǎn)折角觀測時，一般應(yīng)觀測導(dǎo)線前進(jìn)方向的左角。對于閉合導(dǎo)線，若按逆時針方向進(jìn)行觀測，則觀測的導(dǎo)線角既是閉合多邊形的角，又是導(dǎo)線前進(jìn)方向的

20、左角。對于支導(dǎo)線，應(yīng)分別觀測導(dǎo)線前進(jìn)方向的左角和右角，以增加檢核條件。當(dāng)觀測短邊之間的轉(zhuǎn)折角時，測站偏心和目標(biāo)偏心對轉(zhuǎn)折角的影響將十分明顯。因此，應(yīng)對所用儀器、覘牌和光學(xué)對中器進(jìn)行嚴(yán)格檢校，并且要特仔細(xì)進(jìn)行對中和精確照準(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é)對中

21、器。如圖4所示，將全站儀安置在測站i的基座中，帶有覘牌的反射棱鏡安置在后視點(diǎn)i-1和前視點(diǎn)i+1的基座中，進(jìn)行導(dǎo)線測量。遷站時，導(dǎo)線點(diǎn)i和i+1的腳架和基座不動，只取下全站儀和帶有覘牌的反射棱鏡，在導(dǎo)線點(diǎn)i+1上安置全站儀，在導(dǎo)線點(diǎn)i的基座上安置帶有覘牌的反射棱鏡，并將導(dǎo)線點(diǎn)i-1上的腳架遷至導(dǎo)線點(diǎn)i+2處并予以安置，這樣直到測完整條導(dǎo)線為止。在觀測者精心安置儀器的情況下，三聯(lián)腳架法可以減弱儀器和目標(biāo)對中誤差對測角和測距的影響，從而提高導(dǎo)線的觀測精度，減少了坐標(biāo)傳遞誤差。 i-1 i i+1 i+2 圖4 三聯(lián)腳架法導(dǎo)線觀測在城市或工業(yè)區(qū)進(jìn)行導(dǎo)線測量時，可在夜間進(jìn)行作業(yè)，以避免白天作業(yè)時行人、

22、車輛的干擾，夜間作業(yè)，空氣穩(wěn)定、儀器振動小，并可避免太陽暴曬，從而可提高觀測成果的精度。2.3 導(dǎo)線測量的近似平差計(jì)算導(dǎo)線測量的目的是獲得各導(dǎo)線點(diǎn)的平面直角坐標(biāo)計(jì)算的起始數(shù)據(jù)是已知點(diǎn)坐標(biāo)、已知坐標(biāo)方位角，觀測數(shù)據(jù)為觀測角值和觀測邊長。通常情況下，導(dǎo)線平差應(yīng)進(jìn)行嚴(yán)密平差，但對于二級與其以下等級的圖根導(dǎo)線允許對以單一導(dǎo)線、單結(jié)點(diǎn)導(dǎo)線網(wǎng)采用近似平差方法進(jìn)行計(jì)算。導(dǎo)線近似平差的基本思路是將角度誤差和邊長誤差分別進(jìn)行平差計(jì)算，先進(jìn)行角度閉合差的分配，在此基礎(chǔ)在進(jìn)行坐標(biāo)閉合差的分配，通過調(diào)整坐標(biāo)閉合差，以達(dá)到角度的剩余誤差和邊長誤差的目的。在進(jìn)行導(dǎo)線測量平差計(jì)算之前，首先要按照規(guī)要求對外業(yè)觀測成果進(jìn)行檢查

23、和驗(yàn)算，確保觀測成果無誤并符合限差要求，然后對邊長進(jìn)行加常數(shù)改正、乘常數(shù)改正、氣象改正和傾斜改正（改正方法見數(shù)字測圖原理與方法第五章），對角度和邊長進(jìn)行歸心改正（有偏心觀測時），以取消系統(tǒng)誤差的影響2。2.3.1 附合導(dǎo)線的計(jì)算如下圖5為附合導(dǎo)線，其坐標(biāo)平差計(jì)算步驟如下：圖5 附合導(dǎo)線(1) 方位角計(jì)算 （2-1）判斷或0， 以便判斷所在相應(yīng)象限，加上或減去180度，推算出每條導(dǎo)線邊的坐標(biāo)方位角。(2) 由地面實(shí)測點(diǎn)的轉(zhuǎn)角依據(jù)公式計(jì)算其角度閉合差用，依次加減，得出最后的方位角，用計(jì)算的方位角=閉合差值。 (3) 將角度閉合差除以測邊數(shù)分配到各觀測角中角度閉合差調(diào)整中，觀測角為左角時反符號平均分

24、配到各觀測角中,觀測角為右角時，則按閉合差同符號分配到測角，如有小數(shù),按長邊少分,短邊多分原則。(4) 用調(diào)整后的觀測角計(jì)算方位角(5) 坐標(biāo)增量計(jì)算, （2-2） (6) 測站測出各邊的距離，用調(diào)整后的方位角進(jìn)行計(jì)算各點(diǎn)計(jì)算坐標(biāo)(7) 計(jì)算導(dǎo)線點(diǎn)坐標(biāo) 閉合差為: =計(jì)算終點(diǎn)X坐標(biāo)-設(shè)計(jì)終點(diǎn)X坐標(biāo)=計(jì)算終點(diǎn)Y坐標(biāo)-設(shè)計(jì)終點(diǎn)Y坐標(biāo) 將閉合邊長差按路線邊長與總邊長的比例反向分配到坐標(biāo)增量后計(jì)算各點(diǎn)坐標(biāo)。2.3.2 閉合導(dǎo)線的計(jì)算如圖6所示為閉合導(dǎo)線，閉合導(dǎo)線的精度評定與具有兩個連接角的附和導(dǎo)線精度評定一樣，可以采用角度閉合差和導(dǎo)線全長相對閉合差來評定，閉合導(dǎo)線的可靠性較差，在實(shí)際測量中避免單獨(dú)使用

25、。由于角度觀測值存在誤差，使得多邊形角和的計(jì)算值不等于其理論值，而產(chǎn)生角度閉合差，即圖6 閉合導(dǎo)線（1）角度閉合差計(jì)算 （2-3） ， n為實(shí)測角個數(shù) 應(yīng)進(jìn)行重測，如果則可進(jìn)行平差計(jì)算。（2）角度閉合差改正 （反號平均分配） （2-4）（3）計(jì)算改正后的角度 （2-5）（4）推算方位角 公式為： （2-6） 判斷x或或0&(xB-xA)0fwjAB=Rab;end if (yB-yA)0&(xB-xA)0fwjAB=180-Rab;endif (yB-yA)0&(xB-xA)0fwjAB=180+Rab;endif (yB-yA)0fwjAB=360-Rab;endRcd=atand(abs(

26、yD-yC)/(abs(xD-xC); %CD的起算方位角if (yD-yC)0&(xD-xC)0fwjCD=Rcd;end if (yD-yC)0&(xD-xC)0fwjCD=180-Rcd;endif (yD-yC)0&(xD-xC)0fwjCD=180+Rcd;endif (yD-yC)0fwjCD=360-Rcd;end%角度閉合差計(jì)算 fb=fwjAB-fwjCD+sum(jiao1)-n*180; jiao2=jiao1+(fb/6)*(-1);%各觀測角改正后的角度 %真方位角的計(jì)算 for k=1:(n+1) if k=1 fwjz(k)=fwjAB+jiao2(k); end

27、 if k1 fwjz(k)=fwjz(k-1)+jiao2(k)-180; end end %坐標(biāo)閉合差計(jì)算 for i=1:n xz(i)=bian(i)*cosd(fwjz(i); for i=1:n yz(i)=bian(i)*sind(fwjz(i);end end fx=xA+sum(xz)-xC; fy=yA+sum(yz)-yC; %坐標(biāo)閉合差分配 for i=1:n vx(i)=(-1)*fx/sum(bian)*bian(i); end for i=1:n vy(i)=(-1)*fy/sum(bian)*bian(i); end %坐標(biāo)計(jì)算 for i=1:n if i=1

28、 x(1)=xA+xz(1)+vx(1) y(1)=yA+yz(1)+vy(1); end if i1 x(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)線的點(diǎn)數(shù));for i=1:n a=num2str(i); string1=(char(請輸入第),a,(條導(dǎo)線的邊長:); bian(i)=input(string1); b=num2str(i); string2=(char(請輸入第),b,(個轉(zhuǎn)折角:); jiao1=input(string2

29、,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(請輸入起始已知點(diǎn)x的坐標(biāo):);y=input(請輸入起始已知點(diǎn)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;%角度閉合差的計(jì)算fb=sum(jiao2)-(n-2)*180;for i=1:nv(i)=(-1)*fb/n;end%真方位角計(jì)算for i=1:n if i=1 fwjz(1)=fwj2+jiao2(1)+v(1)-180; end i