基于matlab的導(dǎo)線網(wǎng)坐標(biāo)計算畢業(yè)論文-資料下載頁
2025-06-27 17:59本頁面
【正文】 ((abs(yDyC))/(abs(xDxC)))。 %CD的起算方位角if (yDyC)0amp。(xDxC)0fwjCD=Rcd。end if (yDyC)0amp。(xDxC)0fwjCD=180Rcd。endif (yDyC)0amp。(xDxC)0fwjCD=180+Rcd。endif (yDyC)0amp。(xDxC)0fwjCD=360Rcd。end%角度閉合差計算 fb=fwjABfwjCD+sum(jiao1)n*180。 jiao2=jiao1+(fb/6)*(1)。%各觀測角改正后的角度 %真方位角的計算 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)閉合差計算 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)計算 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)線平差程序：%觀測數(shù)據(jù)的輸入n=input(39。請輸入閉合導(dǎo)線的點數(shù)39。)。for i=1:n a=num2str(i)。 string1=(char((39。請輸入第39。)39。,a,(39。條導(dǎo)線的邊長:39。)39。))39。 bian(i)=input(string1)。 b=num2str(i)。 string2=(char((39。請輸入第39。)39。,b,(39。個轉(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。請輸入起始已知點x的坐標(biāo):39。)。y=input(39。請輸入起始已知點y的坐標(biāo):39。)。string3=(char(39。請輸入起始已知方向的起算方位角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。%角度閉合差的計算fb=sum(jiao2)(n2)*180。for i=1:nv(i)=(1)*fb/n。end%真方位角計算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)閉合差計算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)。endfx=xzq+sum(xz)。fy=yzq+sum(yz)。fs=sqrt(fx^2+fy^2)。k=fs/sum(bian)。%坐標(biāo)計算for i=1:nvx=(1)*fx/sum(bian)*bian。vy=(1)*fy/sum(bian)*bian。endfor 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)線平差程序：%觀測數(shù)據(jù)輸入n=input(39。please the 待定點個數(shù):39。)。for i=1:na=num2str(i)。s1=char(39。please input the 邊長 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%方位角推算for j=1:nif j==1fwj(j)=f0+jiao2(j)180。endif j1fwj(j)=fwj(j1)+jiao2(j)180。endend%坐標(biāo)增量計算for k=1:nxzl(k)=S(k)*cosd(fwj(k))。yzl(k)=S(k)*sind(fwj(k))。end%坐標(biāo)推算for l=1:nif l==1x(l)=x0+xzl(l)。y(l)=y0+yzl(l)。endif l1x(l)=x(l1)+xzl(l)。y(l)=y(l1)+yzl(l)。endend
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1