【正文】 f0=degree2+minute2/60+second2/3600%方位角推算for j=1:nif j==1fwj(j)=f0+jiao2(j)180。jiao2(i)=degree1+minute1/60+second1/3600。please the 待定點個數(shù):39。 end if fwjz(i)0 fwjz1(i)=fwjz(i)+360。string3=(char(39。)39。)。(xDxC)0fwjCD=360Rcd。 %AB的起算方位角if (yByA)0amp。 fen=str2num(B(4:5))。)39。還有教過我的所有老師們，你們嚴謹細致、一絲不茍的作風一直是我工作、學習中的榜樣；他們循循善誘的教導和不拘一格的思路給予我無盡的啟迪。19′25″，∠4=105176。 x(l)=x(l1)+xzl(l)。依據(jù)本章介紹的支導線計算方法和基本公式，依據(jù)matlab語言對支導線進行編程，主要的程序代碼如下（具體程序代碼見附錄3）： %觀測數(shù)據(jù)輸入 n=input(39。))。 string1=(char((39。 %真方位角的計算 fwjz(k)=fwjAB+jiao2(k)。 string1=(char((39。 東華理工大學長江學院畢業(yè)設計 基于Matlab的導線計算程序設計3. 基于Matlab的導線計算程序設計 計算程序概述 計算程序特點相對于手工計算，坐標計算程序計算的主要特點是計算速度快、精度高、數(shù)據(jù)處理自動化，從而把人從繁重的計算工作中解放出來。在觀測者精心安置儀器的情況下，三聯(lián)腳架法可以減弱儀器和目標對中誤差對測角和測距的影響，從而提高導線的觀測精度，減少了坐標傳遞誤差。當導線點上應觀測的方向數(shù)多于2個時，應采用方向觀測法進行。（3）擴充能力強，交互性好高版本的的Matlab語言有豐富的庫函數(shù)，在進行復雜的數(shù)序運算時可以直接調(diào)用，而且Matlab的庫函數(shù)同用戶文件在形成上一樣，所以用戶文件也可作為Matlab的庫函數(shù)來調(diào)用。主要缺點是費時、工作量大、還易出錯。文章首先介紹了附合導線、閉合導線、支導線基本概念和計算方法，其次利用Matlab計算機編程語言對三種導線的計算進行編程實現(xiàn)；最后通過實例驗證，本文利用Matlab編寫的程序正確，通過輸入邊長和角度，能夠快速的得到各控制點的準確坐標。而且計算器在計算時小數(shù)位數(shù)是自身進行取舍，所以它的精度可以得到保證，并且比人工逐步計算要高。又如，將Matlab源程序編輯為M文件，由于Matlab磁盤文件也是M文件，所以編輯后的源文件就可以直接運行，而不需要進行編譯和鏈接。由于支導線缺乏檢核條件，故一般只限于地形測量中的圖根導線中采用。 三聯(lián)腳架法導線觀測三聯(lián)腳架法通常使用三個既能安置全站儀又能安置帶有覘牌的基座和腳架，基座應有通用的光學對中器。 (5) 計算坐標增量 , （27）（6）坐標增量閉合差 , （28）（7） 坐標增量閉合差的分配坐標改正數(shù)的計算： （29）計算改正后的坐標增量： （210）（8） 計算坐標 （211） 支導線的計算以下圖7為例，支導線計算步驟如下： 圖7 支導線（1）坐標方位角推算設直線MA的坐標方位角已知為，按方位角推算公式：計算導線邊的坐標方位角。 附合導線計算程序附合導線包括一個連接角的附合導線、兩個連接角的附合導線和無連接角的附合導線等，依據(jù)本文介紹的附合導線計算方法和基本公式，依據(jù)matlab語言對附合導線進行編程,主要的程序代碼如下（具體程序代碼見附錄1）： %觀測數(shù)據(jù)輸入 n=input(39。))39。依據(jù)本章介紹的閉合導線計算方法和基本公式，依據(jù)matlab語言對閉合導線進行編程，主要的程序代碼如下（具體程序代碼見附錄2）： %觀測數(shù)據(jù)的輸入 n=input(39。))39。 X(i)=X(i1)+xz(i1)+vx(i)。 yzl(k)=S(k)*sind(fwj(k))。各觀測內(nèi)角分別為：∠1=89176。致 謝四年大學生活一晃而過，回首走過的歲月，留下了很多美好的回憶，心中倍感充實，當我寫完這篇畢業(yè)論文的時候，有一種如釋重負的感覺，感慨良多。)39。)。請輸入已知點B的y坐標39。endif (yDyC)0amp。 y(i)=y(i1)+yz(i)+vy(i)。)39。y=input(39。end%真方位角計算for i=1:n if i==1 fwjz(1)=fwj2+jiao2(1)+v(1)180。 end if i1 X(i)=X(i1)+xz(i1)+vx(i)。degree1=str2num(jiao1(1:3))。degree2=str2num(fwj0(1:3))。yzl(k)=S(k)*sind(fwj(k))。)。s1=char(39。for i=1:(n1)xz(i)=cosd(fwjz1(i))*bian(i+1)。fwj1=input(string3,39。s39。請輸入第39。%各觀測角改正后的角度 %真方位角的計算 for k=1:(n+1) if k==1 fwjz(k)=fwjAB+jiao2(k)。(xBxA)0fwjAB=180Rab。end%起算數(shù)據(jù)輸入xA=input(39。endfor k=1:(1+n) s=num2str(k)。老師淵博的專業(yè)知識，嚴謹?shù)闹螌W態(tài)度，精益求精的工作作風，誨人不倦的高尚師德，嚴以律己、寬以待人的崇高風范，樸實無華、平易近人的人格魅力對我影響深遠。依據(jù)本文程序設計對上述數(shù)據(jù)進行平差操作，輸出如下成果：（1）基于Matlab的閉合導線程序輸入界面如圖12所示：圖12 程序輸入界面（2）結(jié)果輸出界面如下圖13所示：圖13 結(jié)果輸出界面（3）根據(jù)程序輸出的平差結(jié)果，整理得下表3：表3 閉合導線的計算上表為相應實例基于Matlab程序設計的輸出結(jié)果，誤差均在限差之內(nèi)。59′30″，各觀測角數(shù)據(jù)如下：∠1 =99176。 s1=char(39。 %坐標閉合差計算 xzq=cosd(fwj2)*bian(1)。,a,(39。 yz(i)=bian(i)*sind(fwjz(i))。,a,(39。另外，相對于手工計算，在程序計算時，選擇平差方法的依據(jù)不同。通常情況下，導線平差應進行嚴密平差，但對于二級及其以下等級的圖根導線允許對以單一導線、單結(jié)點導線網(wǎng)采用近似平差方法進行計算?！搴?77。兩條以上導線的匯聚點，稱為導線的結(jié)點。目前市場上,可以買到大型嚴密平差計算程序,但對于許多中小型日常的測量計算，還需要技術(shù)人員自己去開發(fā)；例如：碎部測量、支導線測量、面積計算、解析交會等。 Precision Evaluation。（2）CASIO計算器CASIO計算器是測量人員在野外作業(yè)或礦山井下一線作業(yè)時經(jīng)常使用的便攜式計算器，它與一般計算器所不同的是有個自編程序步。Matlab語言與其他語言相比，較好的解決了上述問題，把編輯、編譯、鏈接和執(zhí)行融為一體。在閉合導線的已知控制點上至少應該有一條定向邊與之相連接。 導線邊長觀測導線邊長可采用電磁波測距儀測量，也可采用全站儀在測取導線角的同時測取導線邊的邊長。 閉合導線的計算如圖6所示為閉合導線，閉合導線的精度評定與具有兩個連接角的附和導線精度評定相同，可以采用角度閉合差和導線