adj.m

利用ＭＡＴＬＡＢ实现了对测量中的平面网

global net ed dd sd dd1 pn x0 y0 m1 m2 m3 ms pp e d sid md g f dir ni si ma s t az 
global aa bb cc rt rr tt 
global a ql pa3 qls w c x y uw0
global pathname net_name s_datafile a_datafile;
s_datafile='.txt'; b_datafile='_bak.txt';  x0y0_datafile='_xy.out';
a_datafile='_adj.out';d_datafile='.dxf';
rddat1;
calcux0y0;
fit2=fopen(strcat(pathname,net_name,a_datafile),'wt');
if(fit2==-1)
    msgbox('Input File or Path is not correct','Warning','warn');
    return;
end          
k=1;
while(k)
	m=m1+m2+m3;
	obnorm;
	c=invsqr(c,2*dd);
	[uw0,k]=adjxy(fit2);
end
ellipse(uw0,fit2)
n=2*dd;
sum=n*(n+1)/2.0;
n1=2*(ed+dd);
sum1=n1*(n1+1)/2.0;
for i=sum1:-1:sum1-sum+1
    c(i)=c(i-(sum1-sum));
end
for i=1:2*ed
     di=(i-1)*(n1-i/2.0);
     for j=i:n1
          if(j==i)
              c(di+j)=0.00000001;
          else
              c(di+j)=0.0;
          end
     end
end
if(m1>0)
    adjs(uw0,fit2);
end
if(m2>0)
    adjd(uw0,fit2);
end
if(m3>0)
    adja(uw0,fit2);
end
fclose(fit2);
open(strcat(pathname,net_name,a_datafile));
return