【正文】 y(l)=y(l1)+yzl(l)。 y(l)=y0+yzl(l)。 yzl(k)=S(k)*sind(fwj(k))。 end if j1 fwj(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)線平差程序： %觀測(cè)數(shù)據(jù)輸入 n=input(39。 end if i1 X(i)=X(i1)+xz(i1)+vx(i)。 end for i=1:n if i==1 X(1)=x+xzq+vx(1)。 東華理工大學(xué)長(zhǎng)江學(xué)院畢業(yè)設(shè)計(jì) 附錄 32 %坐標(biāo)計(jì)算 for i=1:n vx=(1)*fx/sum(bian)*bian。 fs=sqrt(fx^2+fy^2)。 end fx=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。請(qǐng)輸入起始已知方向的起算方位角 39。)。 東華理工大學(xué)長(zhǎng)江學(xué)院畢業(yè)設(shè)計(jì) 附錄 31 y=input(39。請(qǐng)輸入起始已知點(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。請(qǐng)輸入第 39。 for i=1:n a=num2str(i)。請(qǐng)輸入閉合導(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)。 東華理工大學(xué)長(zhǎng)江學(xué)院畢業(yè)設(shè)計(jì) 附錄 30 %坐標(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))。%各觀測(cè)角改正后的角度 %真方位角的計(jì)算 for k=1:(n+1) if k==1 fwjz(k)=fwjAB+jiao2(k)。 end %角度閉合差計(jì)算 fb=fwjABfwjCD+sum(jiao1)n*180。 end if (yDyC)0amp。 end if (yDyC)0amp。 end if (yDyC)0amp。 %CD 的起算方位角 if (yDyC)0amp。(xBxA)0 fwjAB=360Rab。(xBxA)0 fwjAB=180+Rab。(xBxA)0 fwjAB=180Rab。(xBxA)0 fwjAB=Rab。 %起算方位角的計(jì)算 Rab=atand((abs(yByA))/(abs(xBxA)))。請(qǐng)輸入已知點(diǎn) B的 y坐標(biāo) 39。)。 xB=input(39。請(qǐng)輸入已知點(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。請(qǐng)輸入第 39。 end for k=1:(1+n) s=num2str(k)。))39。條導(dǎo)線邊長(zhǎng) 39。)39。 string1=(char((39。)。2020 年 01期 . [6] 崔利、武文波 .測(cè)繪領(lǐng)域中 MATLAB的應(yīng)用 [J]，遼寧工程技術(shù)大學(xué)學(xué)報(bào)， 2020年 S1 期 . [7] 陳以新 .Matlab的幾則程序設(shè)計(jì)經(jīng)驗(yàn) [J]，計(jì)算機(jī)應(yīng)用研究， 1999年 09期 . [8] 孫慶華、熊偉、張憲柱、王磊 .基于 Excel 和 MATLA 導(dǎo)線網(wǎng)平差 [J]，北京測(cè)繪， 2020年 03 期， 132. [9] 白征東 .Matlab在測(cè)量平差教學(xué)中的應(yīng)用 [J]，測(cè)繪通報(bào)， 2020年 11期 ,123. [10] 高彩云、高寧、孫敬克 .MATLAB 在測(cè)量平差課程教學(xué)中的嘗試 [J]，科技信息， 2020 年 05 期， 112. [11] 尹暉主編 .《測(cè)繪工程專業(yè)英語(yǔ)》，武漢：武漢大學(xué)出版社， 2020 年 1 月第 1 版， 143. [12] Ingenieurvermessung 2020,14th International course on Engineering Surveying[C].ETH Zuerich,2020. [13] Kahmen Measurement Robots in Engineering Surveys[J]. I SPRS Congres,Washington,1992. 東華理工大學(xué)長(zhǎng)江學(xué)院畢業(yè)設(shè)計(jì) 附錄 28 附錄 1 基于 Matlab 的附合導(dǎo)線平差程序： %觀測(cè)數(shù)據(jù)輸入 n=input(39。本論文從選題到完成，每一步都是在指導(dǎo)的悉心指導(dǎo)下完成的，傾注了指導(dǎo)老師大量的心血。老師淵博的專業(yè)知 識(shí)，嚴(yán)謹(jǐn)?shù)闹螌W(xué)態(tài)度，精益求精的工作作風(fēng)，誨人不倦的高尚師德，嚴(yán)以律己、寬以待人的崇高風(fēng)范，樸實(shí)無(wú)華、平易近人的人格魅力對(duì)我影響深遠(yuǎn)。 感謝四年中陪伴在我身邊的同學(xué)、朋友，感謝你們?yōu)槲姨岢龅挠幸娴慕ㄗh和意見(jiàn)，有了你們的支持、鼓勵(lì)和幫助，我才能充實(shí)的度過(guò)了四年的學(xué)習(xí)生活。他在忙碌的教學(xué)工作中擠出時(shí)間來(lái)審查、修改我的論文。 致 謝 東華理工大學(xué)長(zhǎng)江學(xué)院畢業(yè)設(shè)計(jì) 結(jié)論 26 四年大學(xué)生活一晃而過(guò)，回首走過(guò)的歲月，留下 了很多美好的回憶，心中倍感充實(shí)，當(dāng)我寫完這篇畢業(yè)論文的時(shí)候，有一種如釋重負(fù)的感覺(jué)，感慨良多。 文章首先介紹了幾種導(dǎo)線坐標(biāo)計(jì)算方法的優(yōu)缺點(diǎn)，也敘述了 Matlab 的特點(diǎn) ，其次依據(jù)《數(shù)字測(cè)圖原理與方法》介紹了附合導(dǎo)線、閉合導(dǎo)線、支導(dǎo)線基本概念和計(jì)算方法，詳細(xì)描述了導(dǎo)線布設(shè)和導(dǎo)線觀測(cè)，導(dǎo)線觀測(cè)包括轉(zhuǎn)折角和邊長(zhǎng)觀測(cè)，也敘述了三聯(lián)腳架法的相應(yīng)內(nèi)容?！洹澹?/m /m /m /m /m∑= 5 3 7 . 4 5M237 59 30A （P 1 ） 99 01 00 157 00 30 1 167 45 36 144 46 06 2 123 11 24 87 57 30 3 ∑ 389 58 00 上表為相應(yīng)實(shí)例基于 Matlab 程序設(shè)計(jì)的輸出結(jié)果，誤差均在限差之內(nèi)。 依據(jù)本文程序設(shè)計(jì)對(duì)上述數(shù)據(jù)進(jìn)行平差操作，輸出如下成果： (1）基于 Matlab 的支導(dǎo)線程序輸入界面 ,如下圖 15所示： 圖 15 程序輸入界面 東華理工大學(xué)長(zhǎng)江學(xué)院畢業(yè)設(shè)計(jì) 三種導(dǎo)線坐標(biāo)平差實(shí)例 23 （ 3）結(jié)果輸出界面如下圖 16 所示： 圖 16 結(jié)果 輸出界面 東華理工大學(xué)長(zhǎng)江學(xué)院畢業(yè)設(shè)計(jì) 三種導(dǎo)線坐標(biāo)平差實(shí)例 24 （ 3）根據(jù)程序輸出的平差結(jié)果，整理得下表 4： 表 4 支導(dǎo)線的計(jì)算 點(diǎn)名 觀測(cè)角 坐標(biāo)方位角 邊長(zhǎng)（S ） △X △Y x y（176。 45′ 36″， 123176。 59′ 30″ ， SA1=， S12=， S23=，轉(zhuǎn)折角分別為：099176?！洹澹?/m /m /m /m /m∑= 2 7 2 . 5 8fx=,fy=,fs=,k=512 +1091 09 263 +1089 19 25 6 +64 000 00 00 7 +9 271 09 36 12 +16 180 29 11 +10105 37 40 106 07 01 7 +95 +10163 53 471 +989 58 53∑ 539 59 1190 00 58 6 +7 上表為相應(yīng)實(shí)例基于 Matlab 程序設(shè)計(jì)的輸出結(jié)果，誤差均在限差之內(nèi)。 依據(jù)本文程序設(shè)計(jì)對(duì)上述數(shù)據(jù)進(jìn)行平差操作，輸出如下成果： 東華理工大學(xué)長(zhǎng)江學(xué)院畢業(yè)設(shè)計(jì) 三種導(dǎo)線坐標(biāo)平差實(shí)例 19 （ 1）基于 Matlab的閉合導(dǎo)線程序輸入界面如圖 12所示： 圖 12 程序輸入界面 東華理工大學(xué)長(zhǎng)江學(xué)院畢業(yè)設(shè)計(jì) 三種導(dǎo)線坐標(biāo)平差實(shí)例 20 （ 2）結(jié)果輸出界面如下圖 13 所示： 圖 13 結(jié)果輸出界面 東華理工大學(xué)長(zhǎng)江學(xué)院畢業(yè)設(shè)計(jì)