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

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

2、念和計(jì)算方法，其次利用Matlab計(jì)算機(jī)編程語(yǔ)言對(duì)三種導(dǎo)線的計(jì)算進(jìn)行編程實(shí)現(xiàn)；最后通過實(shí)例驗(yàn)證，本文利用Matlab編寫的程序正確，通過輸入邊長(zhǎng)和角度，能夠快速的得到各控制點(diǎn)的準(zhǔn)確坐標(biāo)。關(guān)鍵詞：Matlab；導(dǎo)線計(jì)算；精度評(píng)價(jià)；計(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)線的觀測(cè) PAGEREF _Toc25147 4 HYPERLINK l _Toc29039 2.2.1 轉(zhuǎn)折角的觀測(cè) PAGEREF _Toc29039 4 HYPERLINK l _Toc1164 2.2.2 導(dǎo)線邊長(zhǎng)觀測(cè) PAGEREF _Toc1164 5 HYPERLINK l _Toc10873 2.2.3 三聯(lián)腳架法導(dǎo)線觀測(cè) PAGEREF _Toc10873 5 HYPERLINK l _Toc30360 2.3 導(dǎo)線測(cè)量的近似平差計(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ì)算是在所有測(cè)量工作中經(jīng)常遇見

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

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

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

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

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

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

17、ab語(yǔ)言有豐富的庫(kù)函數(shù)，在進(jìn)行復(fù)雜的數(shù)序運(yùn)算時(shí)可以直接調(diào)用，而且Matlab的庫(kù)函數(shù)同用戶文件在形成上一樣，所以用戶文件也可作為Matlab的庫(kù)函數(shù)來調(diào)用。因而，用戶可以根據(jù)自己的需要方便地建立和擴(kuò)充新的庫(kù)函數(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、線起始于一個(gè)已知控制點(diǎn)而終止于另一個(gè)已知控制點(diǎn)。已知控制點(diǎn)上可以有一條或幾條定向邊與之相連接，也可以沒有定向與之相連接。圖1 附合導(dǎo)線（2） 閉合導(dǎo)線如圖2所示，由一個(gè)已知控制點(diǎn)出發(fā)，最終又回到這一點(diǎn)，形成一個(gè)閉合的多邊形。在閉合導(dǎo)線的已知控制點(diǎn)上至少應(yīng)該有一條定向邊與之相連接。由于閉合導(dǎo)線是一種可靠性極差的控制網(wǎng)圖形，在實(shí)際測(cè)量工作中應(yīng)避免單獨(dú)使用。 圖2 閉合導(dǎo)線（3） 支導(dǎo)線如圖3所示，從一個(gè)已知控制點(diǎn)出發(fā)，既不附合于另一個(gè)已知控制點(diǎn)，也不閉合于原來的起始控制點(diǎn)。由于支導(dǎo)線缺乏檢核條件，故一般只限于地形測(cè)量中的圖根導(dǎo)線中采用。 圖3 支導(dǎo)線2.2 導(dǎo)線的觀測(cè)導(dǎo)線的觀測(cè)包括轉(zhuǎn)折角的觀測(cè)和導(dǎo)

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

20、左角。對(duì)于支導(dǎo)線，應(yīng)分別觀測(cè)導(dǎo)線前進(jìn)方向的左角和右角，以增加檢核條件。當(dāng)觀測(cè)短邊之間的轉(zhuǎn)折角時(shí)，測(cè)站偏心和目標(biāo)偏心對(duì)轉(zhuǎn)折角的影響將十分明顯。因此，應(yīng)對(duì)所用儀器、覘牌和光學(xué)對(duì)中器進(jìn)行嚴(yán)格檢校，并且要特仔細(xì)進(jìn)行對(duì)中和精確照準(zhǔn)。2.2.2 導(dǎo)線邊長(zhǎng)觀測(cè)導(dǎo)線邊長(zhǎng)可采用電磁波測(cè)距儀測(cè)量，也可采用全站儀在測(cè)取導(dǎo)線角的同時(shí)測(cè)取導(dǎo)線邊的邊長(zhǎng)。導(dǎo)線邊長(zhǎng)應(yīng)對(duì)向觀測(cè)，以增加檢核條件。電磁波測(cè)距儀測(cè)量的通常是斜距，還需觀測(cè)豎直角，用以將傾斜距離改變?yōu)樗骄嚯x，必要時(shí)還應(yīng)將其歸算到橢球面上和高斯平面上。2.2.3 三聯(lián)腳架法導(dǎo)線觀測(cè)三聯(lián)腳架法通常使用三個(gè)既能安置全站儀又能安置帶有覘牌的基座和腳架，基座應(yīng)有通用的光學(xué)對(duì)中

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

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

23、和驗(yàn)算，確保觀測(cè)成果無誤并符合限差要求，然后對(duì)邊長(zhǎng)進(jìn)行加常數(shù)改正、乘常數(shù)改正、氣象改正和傾斜改正（改正方法見數(shù)字測(cè)圖原理與方法第五章），對(duì)角度和邊長(zhǎng)進(jìn)行歸心改正（有偏心觀測(cè)時(shí)），以取消系統(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í)測(cè)點(diǎn)的轉(zhuǎn)角依據(jù)公式計(jì)算其角度閉合差用，依次加減，得出最后的方位角，用計(jì)算的方位角=閉合差值。 (3) 將角度閉合差除以測(cè)邊數(shù)分配到各觀測(cè)角中角度閉合差調(diào)整中，觀測(cè)角為左角時(shí)反符號(hào)平均分

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

25、。由于角度觀測(cè)值存在誤差，使得多邊形角和的計(jì)算值不等于其理論值，而產(chǎn)生角度閉合差，即圖6 閉合導(dǎo)線（1）角度閉合差計(jì)算 （2-3） ， n為實(shí)測(cè)角個(gè)數(shù) 應(yīng)進(jìn)行重測(cè)，如果則可進(jìn)行平差計(jì)算。（2）角度閉合差改正 （反號(hào)平均分配） （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);%各觀測(cè)角改正后的角度 %真方位角的計(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)線平差程序：%觀測(cè)數(shù)據(jù)的輸入n=input(請(qǐng)輸入閉合導(dǎo)線的點(diǎn)數(shù));for i=1:n a=num2str(i); string1=(char(請(qǐng)輸入第),a,(條導(dǎo)線的邊長(zhǎng):); bian(i)=input(string1); b=num2str(i); string2=(char(請(qǐng)輸入第),b,(個(gè)轉(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(請(qǐng)輸入起始已知點(diǎn)x的坐標(biāo):);y=input(請(qǐng)輸入起始已知點(diǎn)y的坐標(biāo):);string3=(char(請(qǐng)輸入起始已知方向的起算方位角);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