poly_close.m

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

function poly_close(e,d,sid,g,f,dir,si,ni,ed,dd,pn,x0,y0)
%闭合导线闭合差计算
global pathname net_name c4_datafile
fit2=fopen(strcat(pathname,net_name,c4_datafile),'wt');
[p1,p2,s,dir12,dir21]=orientcollect(e,d,sid,g,f,dir,si,ni,ed,dd,pn);%整理平面网观测数据
cpp=ccpoints(p1',p2',ed+dd,length(p1),s);%按余支搜索法进行闭合环搜索
fprintf(fit2,'     闭合导线条件\n');
fprintf(fit2,'环号     闭合环点号       X坐标闭合差    Y坐标闭合差    角度闭合差    长度闭合差   相对闭合差\n');
sitefinal2=[];
for i=1:size(cpp,2)
	site=cpp{i};
	[closediff,err]=orientlead(4,site,si,ni,pn,e,d,sid,g,f,dir,x0,y0);%计算平面网闭合导线的闭合差
	for k=1:length(site)
       site(k)=pn(site(k));
	end
    fprintf(fit2,'---------------------------------------------\n');
	fprintf(fit2,' %5d  %5d%5d%5d%5d%5d%5d%5d',i,site);
	fprintf(fit2,' %10.3f%10.3f%10.3f%10.3f   %10.0f\n',closediff(1:5));
end
fprintf(fit2,'---------------------------------------------\n');
fclose(fit2);
open(strcat(pathname,net_name,c4_datafile));
return