【導(dǎo)讀】東華理工大學(xué)長(zhǎng)江學(xué)院。東華理工大學(xué)長(zhǎng)江學(xué)院畢業(yè)設(shè)計(jì)摘要。導(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)線基本概念和計(jì)。算方法，其次利用Matlab計(jì)算機(jī)編程語言對(duì)三種導(dǎo)線的計(jì)算進(jìn)行編程實(shí)現(xiàn)；最后。通過實(shí)例驗(yàn)證，本文利用Matlab編寫的程序正確，通過輸入邊長(zhǎng)和角度，能夠快。速的得到各控制點(diǎn)的準(zhǔn)確坐標(biāo)。東華理工大學(xué)長(zhǎng)江學(xué)院畢業(yè)設(shè)計(jì)ABSTRACT

　　

【正文】 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。請(qǐng)輸入導(dǎo)線邊數(shù) 39。)。 for i=1:n a=num2str(i)。 string1=(char((39。請(qǐng)輸入第 39。)39。,a,(39。條導(dǎo)線邊長(zhǎng) 39。)39。))39。 bian(i)=input(string1)。 end for k=1:(1+n) s=num2str(k)。 string2=(char((39。請(qǐng)輸入第 39。)39。,s,(39。個(gè)轉(zhuǎn)折角 39。)39。))39。 jiao=input(string2,39。s39。)。 B=jiao。 du=str2num(B(1:3))。 fen=str2num(B(4:5))。 miao=str2num(B(6:7))。 jiao1(k)=du+fen/60+miao/3600。 end %起算數(shù)據(jù)輸入 xA=input(39。請(qǐng) 輸入已知點(diǎn) A的 x坐標(biāo) 39。)。 yA=input(39。請(qǐng)輸入已知點(diǎn) A的 y坐標(biāo) 39。)。 xB=input(39。請(qǐng)輸入已知點(diǎn) B的 x坐標(biāo) 39。)。 yB=input(39。請(qǐng)輸入已知點(diǎn) B的 y坐標(biāo) 39。)。 %起算方位角的計(jì)算 Rab=atand((abs(yByA))/(abs(xBxA)))。 %AB 的起算方位角 if (yByA)0amp。(xBxA)0 fwjAB=Rab。 end if (yByA)0amp。(xBxA)0 fwjAB=180Rab。 end if (yByA)0amp。(xBxA)0 fwjAB=180+Rab。 end 東華理工大學(xué)長(zhǎng)江學(xué)院畢業(yè)設(shè)計(jì) 附錄 29 if (yByA)0amp。(xBxA)0 fwjAB=360Rab。 end Rcd=atand((abs(yDyC))/(abs(xDxC)))。 %CD 的起算方位角 if (yDyC)0amp。(xDxC)0 fwjCD=Rcd。 end if (yDyC)0amp。(xDxC)0 fwjCD=180Rcd。 end if (yDyC)0amp。(xDxC)0 fwjCD=180+Rcd。 end if (yDyC)0amp。(xDxC)0 fwjCD=360Rcd。 end %角度閉合差計(jì)算 fb=fwjABfwjCD+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 if k1 fwjz(k)=fwjz(k1)+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。 東華理工大學(xué)長(zhǎng)江學(xué)院畢業(yè)設(shè)計(jì) 附錄 30 %坐標(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 x(1)=xA+xz(1)+vx(1) y(1)=yA+yz(1)+vy(1)。 end if i1 x(i)=x(i1)+xz(i)+vx(i)。 y(i)=y(i1)+yz(i)+vy(i)。 end 附錄 2 基于 Matlab 的閉合導(dǎo)線平差程序： %觀測(cè)數(shù)據(jù)的輸入 n=input(39。請(qǐng)輸入閉合導(dǎo)線的點(diǎn)數(shù) 39。)。 for i=1:n a=num2str(i)。 string1=(char((39。請(qǐng)輸入第 39。)39。,a,(39。條導(dǎo)線的邊長(zhǎng) :39。)39。))39。 bian(i)=input(string1)。 b=num2str(i)。 string2=(char((39。請(qǐng)輸入第 39。)39。,b,(39。個(gè)轉(zhuǎn)折角 :39。)39。))39。 jiao1=input(string2,39。s39。)。 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(39。請(qǐng)輸入起始已知點(diǎn) x 的坐標(biāo) :39。)。 東華理工大學(xué)長(zhǎng)江學(xué)院畢業(yè)設(shè)計(jì) 附錄 31 y=input(39。請(qǐng)輸入起始已知點(diǎn) y 的坐標(biāo) :39。)。 string3=(char(39。請(qǐng)輸入起始已知方向的起算方位角 39。))。 fwj1=input(string3,39。s39。)。 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)(n2)*180。 for i=1:n v(i)=(1)*fb/n。 end %真方位角計(jì)算 for i=1:n if i==1 fwjz(1)=fwj2+jiao2(1)+v(1)180。 end if i1 fwjz(i)=fwjz(i1)+jiao2(i)+v(i)180。 end if fwjz(i)0 fwjz1(i)=fwjz(i)。 end if fwjz(i)0 fwjz1(i)=fwjz(i)+360。 end %坐標(biāo)閉合差計(jì)算 xzq=cosd(fwj2)*bian(1)。 yzq=sind(fwj2)*bian(1)。 for i=1:(n1) xz(i)=cosd(fwjz1(i))*bian(i+1)。 yz(i)=sind(fwjz1(i))*bian(i+1)。 end fx=xzq+sum(xz)。 fy=yzq+sum(yz)。 fs=sqrt(fx^2+fy^2)。 k=fs/sum(bian)。 東華理工大學(xué)長(zhǎng)江學(xué)院畢業(yè)設(shè)計(jì) 附錄 32 %坐標(biāo)計(jì)算 for i=1:n vx=(1)*fx/sum(bian)*bian。 vy=(1)*fy/sum(bian)*bian。 end for i=1:n if i==1 X(1)=x+xzq+vx(1)。 Y(1)=y+yzq+vy(1)。 end if i1 X(i)=X(i1)+xz(i1)+vx(i)。 Y(i)=Y(i1)+yz(i1)+vy(i)。 end 附錄 3 基于 Matlab 的支導(dǎo)線平差程序： %觀測(cè)數(shù)據(jù)輸入 n=input(39。please the 待定點(diǎn)個(gè)數(shù) :39。)。 for i=1:n a=num2str(i)。 s1=char(39。please input the 邊長(zhǎng) of a:39。)。 S(i)=input(s1)。 s2=char(39。please input the 轉(zhuǎn)折角 of a:39。)。 jiao1=input(s2,39。s39。)。 degree1=str2num(jiao1(1:3))。 minute1=str2num(jiao1(4:5))。 second1=str2num(jiao1(6:7))。 jiao2(i)=degree1+minute1/60+second1/3600。 end %已知數(shù)據(jù)輸入 x0=input(39。please input the known point of x0:39。)。 y0=input(39。please input the known point of y0:39。)。 s3=char(39。please input the known azimuth:39。)。 fwj0=input(s3,39。s39。)。 degree2=str2num(fwj0(1:3))。 minute2=str2num(fwj0(4:5))。 second2=str2num(fwj0(6:7))。 f0=degree2+minute2/60+second2/3600 %方位角推算 東華理工大學(xué)長(zhǎng)江學(xué)院畢業(yè)設(shè)計(jì) 附錄 33 for j=1:n if j==1 fwj(j)=f0+jiao2(j)180。 end if j1 fwj(j)=fwj(j1)+jiao2(j)180。 end end %坐標(biāo)增量計(jì)算 for k=1:n xzl(k)=S(k)*cosd(fwj(k))。 yzl(k)=S(k)*sind(fwj(k))。 end %坐標(biāo)推算 for l=1:n if l==1 x(l)=x0+xzl(l)。 y(l)=y0+yzl(l)。 end if l1 x(l)=x(l1)+xzl(l)。 y(l)=y(l1)+yzl(l)。 end end