level_close2.m

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

function level_close2     %余枝搜索法
global pathname net_name s_datafile b_datafile a1_datafile;;
global ed dd sd pn gd h0 k1 k2 h1 dh s  k11 k12;
fit2=fopen(strcat(pathname,net_name,'_close.txt'),'wt');
fprintf(fit2,'                水准路线闭合差\n');
fprintf(fit2,'=================================================================\n');
fprintf(fit2,'环号     闭合环点号      闭合环长度    闭合差    允许闭合差\n');
cpp=ccpoints(k1,k2,sd,gd,s);
m=size(cpp,2);n=[];
for i=1:m
n=cpp{i};
j=size(n,2);
no=0;s0=0;
  for jj=1:j-1
      k=size(k1,1);
          for kk=1:k
        if (n(jj)==k1(kk))& (n(jj+1)==k2(kk))
            no=no+h1(kk);
            s0=s0+s(kk);
        end
       if (n(jj)==k2(kk))& (n(jj+1)==k1(kk)) 
            no=no-h1(kk);
            s0=s0+s(kk);
        else
            continue;
        end
    end
end
clo(i)=no;           %计算每个环的闭合差   
ss(i)=s0;                 %闭合环的长度
sx(i)=0.02*sqrt(ss(i));       %计算每个环的允许闭合差
     fprintf(fit2,'%3d  ',i);
     fprintf(fit2,'%3d %3d %3d %3d ',n);   
     fprintf(fit2,' %10.3f %10.3f %10.3f\n',ss(i),clo(i),sx(i));
 end
 fprintf(fit2,'============================================================\n');
fclose(fit2);
open(strcat(pathname,net_name,'_close.txt'));
end