adja.m

来自「利用ＭＡＴＬＡＢ实现了对测量中的平面网」· M 代码 · 共 25 行

function adja(uw,fit2)     %方位角精度评定
global m3 x y s t az
lo=206264.8062470964;
fprintf(fit2,'adjusted azimath \n  i    pn(s)     pn(t)        az       va(")        az+va       ma');
for i=1:m3
	dx=x(t(i))-x(s(i));
	dy=y(t(i))-y(s(i));
	a0=alfa(dx,dy);
	if(t(i)<s(i))
		[maa,mss]=trel(uw,t(i),s(i));
	else
		[maa,mss]=trel(uw,s(i),t(i));
	end
	if((a0-az(i))>pi)
		va=(a0-az(i)-2.0*pi)*lo;
	else
		va=(a0-az(i))*lo;
	end
	az1=az(i)+va/lo;
	az1=wg(az1);
	az(i)=wg(az(i));
	fprintf(fit2,'      %3d %8d %8d %13.5f %8.2f %12.5f %8.2f',i,pn(s(i)),pn(t(i)),az(i),va,az1,maa);
end
return