adjs.m

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

function adjs(uw,fit2)    %边长观测精度
global ed dd sd pn m1 e d sid x y
fprintf(fit2,'adjusted sides \n');
fprintf(fit2,'   i       pn(e)    pn(d)       side        vs(cm)        side+vs        ms\n');
for i=1:m1
	dx=x(d(i))-x(e(i));
	dy=y(d(i))-y(e(i));
	if(e(i)<d(i))
        [maa,mss]=trel(uw,e(i),d(i));
	else
        [maa,mss]=trel(uw,d(i),e(i));
	end
	ss=dx*dx+dy*dy;
	ss=sqrt(ss);
	vs=(ss-sid(i))*100.0;
	sid1=sid(i)+vs/100.0;
	fprintf(fit2,' %3d   %8d %8d %15.4f %10.4f %15.4f %8.2f\n',i,pn(e(i)),pn(d(i)),sid(i),vs,sid1,mss);
end
return