【正文】 y(l)=y(l1)+yzl(l)。y(l)=y0+yzl(l)。yzl(k)=S(k)*sind(fwj(k))。endif j1fwj(j)=fwj(j1)+jiao2(j)180。second2=str2num(fwj0(6:7))。degree2=str2num(fwj0(1:3))。s39。)。s3=char(39。please input the known point of y0:39。)。end%已知數(shù)據(jù)輸入x0=input(39。second1=str2num(jiao1(6:7))。degree1=str2num(jiao1(1:3))。s39。)。s2=char(39。)。s1=char(39。)。 end 附錄3 基于Matlab的支導(dǎo)線平差程序：%觀測數(shù)據(jù)輸入n=input(39。 end if i1 X(i)=X(i1)+xz(i1)+vx(i)。endfor i=1:n if i==1 X(1)=x+xzq+vx(1)。%坐標(biāo)計(jì)算for i=1:nvx=(1)*fx/sum(bian)*bian。fs=sqrt(fx^2+fy^2)。endfx=xzq+sum(xz)。for i=1:(n1)xz(i)=cosd(fwjz1(i))*bian(i+1)。 end%坐標(biāo)閉合差計(jì)算xzq=cosd(fwj2)*bian(1)。 end if fwjz(i)0 fwjz1(i)=fwjz(i)。end%真方位角計(jì)算for i=1:n if i==1 fwjz(1)=fwj2+jiao2(1)+v(1)180。%角度閉合差的計(jì)算fb=sum(jiao2)(n2)*180。miao=str2num(c(6:7))。du=str2num(c(1:3))。)。fwj1=input(string3,39。請輸入起始已知方向的起算方位角39。)。y=input(39。請輸入起始已知點(diǎn)x的坐標(biāo):39。 jiao2(i)=du+fen/60+miao/3600。 fen=str2num(c(4:5))。 c=jiao1。s39。))39。個(gè)轉(zhuǎn)折角:39。)39。 string2=(char((39。 bian(i)=input(string1)。)39。,a,(39。請輸入第39。for i=1:n a=num2str(i)。請輸入閉合導(dǎo)線的點(diǎn)數(shù)39。 y(i)=y(i1)+yz(i)+vy(i)。 end %坐標(biāo)計(jì)算 for i=1:n if i==1 x(1)=xA+xz(1)+vx(1) y(1)=yA+yz(1)+vy(1)。 %坐標(biāo)閉合差分配 for i=1:n vx(i)=(1)*fx/sum(bian)*bian(i)。end end fx=xA+sum(xz)xC。 end end %坐標(biāo)閉合差計(jì)算 for i=1:n xz(i)=bian(i)*cosd(fwjz(i))。%各觀測角改正后的角度 %真方位角的計(jì)算 for k=1:(n+1) if k==1 fwjz(k)=fwjAB+jiao2(k)。end%角度閉合差計(jì)算 fb=fwjABfwjCD+sum(jiao1)n*180。endif (yDyC)0amp。endif (yDyC)0amp。end if (yDyC)0amp。 %CD的起算方位角if (yDyC)0amp。(xBxA)0fwjAB=360Rab。(xBxA)0fwjAB=180+Rab。(xBxA)0fwjAB=180Rab。(xBxA)0fwjAB=Rab。%起算方位角的計(jì)算Rab=atand((abs(yByA))/(abs(xBxA)))。請輸入已知點(diǎn)B的y坐標(biāo)39。)。xB=input(39。請輸入已知點(diǎn)A的y坐標(biāo)39。)。end%起算數(shù)據(jù)輸入xA=input(39。 miao=str2num(B(6:7))。 du=str2num(B(1:3))。)。 jiao=input(string2,39。)39。,s,(39。請輸入第39。endfor k=1:(1+n) s=num2str(k)。))39。條導(dǎo)線邊長39。)39。 string1=(char((39。)。2004年01期.[6] 崔利、[J]，遼寧工程技術(shù)大學(xué)學(xué)報(bào)，2005年S1期.[7] [J]，計(jì)算機(jī)應(yīng)用研究，1999年09期.[8] 孫慶華、熊偉、張憲柱、[J]，北京測繪，2009年03期，132.[9] [J]，測繪通報(bào)，2009年11期,123.[10] 高彩云、高寧、[J]，科技信息，2009年05期，112.[11] 尹暉主編.《測繪工程專業(yè)英語》，武漢：武漢大學(xué)出版社，2005年1月第1版，143.[12] Ingenieurvermessung 2004,14th International course on Engineering Surveying[C].ETH Zuerich,2004.[13] Kahmen Measurement Robots in Engineering Surveys[J]. I SPRS Congres,Washington,1992.東華理工大學(xué)長江學(xué)院畢業(yè)設(shè)計(jì) 附錄 附錄1 基于Matlab的附合導(dǎo)線平差程序：%觀測數(shù)據(jù)輸入n=input(39。本論文從選題到完成，每一步都是在指導(dǎo)的悉心指導(dǎo)下完成的，傾注了指導(dǎo)老師大量的心血。老師淵博的專業(yè)知識，嚴(yán)謹(jǐn)?shù)闹螌W(xué)態(tài)度，精益求精的工作作風(fēng)，誨人不倦的高尚師德，嚴(yán)以律己、寬以待人的崇高風(fēng)范，樸實(shí)無華、平易近人的人格魅力對我影響深遠(yuǎn)。感謝四年中陪伴在我身邊的同學(xué)、朋友，感謝你們?yōu)槲姨岢龅挠幸娴慕ㄗh和意見，有了你們的支持、鼓勵(lì)和幫助，我才能充實(shí)的度過了四年的學(xué)習(xí)生活。他在忙碌的教學(xué)工作中擠出時(shí)間來審查、修改我的論文。致 謝四年大學(xué)生活一晃而過，回首走過的歲月，留下了很多美好的回憶，心中倍感充實(shí)，當(dāng)我寫完這篇畢業(yè)論文的時(shí)候，有一種如釋重負(fù)的感覺，感慨良多。 文章首先介紹了幾種導(dǎo)線坐標(biāo)計(jì)算方法的優(yōu)缺點(diǎn)，也敘述了Matlab的特點(diǎn)，其次依據(jù)《數(shù)字測圖原理與方法》介紹了附合導(dǎo)線、閉合導(dǎo)線、支導(dǎo)線基本概念和計(jì)算方法，詳細(xì)描述了導(dǎo)線布設(shè)和導(dǎo)線觀測，導(dǎo)線觀測包括轉(zhuǎn)折角和邊長觀測，也敘述了三聯(lián)腳架法的相應(yīng)內(nèi)容。依據(jù)本文程序設(shè)計(jì)對上述數(shù)據(jù)進(jìn)行平差操作，輸出如下成果：(1）基于Matlab的支導(dǎo)線程序輸入界面,如下圖15所示：圖15 程序輸入界面（3）結(jié)果輸出界面如下圖16所示：圖16 結(jié)果輸出界面（3）根據(jù)程序輸出的平差結(jié)果，整理得下表4：表4 支導(dǎo)線的計(jì)算上表為相應(yīng)實(shí)例基于Matlab程序設(shè)計(jì)的輸出結(jié)果，誤差均在限差之內(nèi)。45′36″，123176。59′30″，SA1=，S12=，S23=，轉(zhuǎn)折角分別為：099176。依據(jù)本文程序設(shè)計(jì)對上述數(shù)據(jù)進(jìn)行平差操作，輸出如下成果：（1）基于Matlab的閉合導(dǎo)線程序輸入界面如圖12所示：圖12 程序輸入界面（2）結(jié)果輸出界面如下圖13所示：圖13 結(jié)果輸出界面（3）根據(jù)程序輸出的平差結(jié)果，整理得下表3：表3 閉合導(dǎo)線的計(jì)算上表為相應(yīng)實(shí)例基于Matlab程序設(shè)計(jì)的輸出結(jié)果，誤差均在限差之內(nèi)。37′40″，∠5=163176。9′26″，∠3=89176。各觀測內(nèi)角分別為：∠1=89176。α12=0176。27′24″，其導(dǎo)線各邊長分別為：S12==,S34=,S45=,S5C=.依據(jù)本文程序設(shè)計(jì)對上述數(shù)據(jù)進(jìn)行平差