poly_close.asv

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

function poly_close(e,d,sid,g,f,dir,si,ni,ed,dd,pn,x0,y0)
%闭合导线闭合差计算
global pathname net_name c4_datafile
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 aa bb cc rt rr tt
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,'===============================================\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);%计算平面网闭合导线的闭合差
    fprintf(fit2,' %5d%5d   %5d%5d%5d%5d%5d%5d',i,site);
    fprintf(fit2,' %10.3f %10.3f %10.3f%10.3f   %10.3f\n',closediff(1:5));
end
fprintf(fit2,'===============================================\n');
fclose(fit2);
open(strcat(pathname,net_name,c4_datafile));
return