poletermerr.m

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

function[err,maxdiff]=poletermerr(cpp2,mr,tr,g,f,dir,si,ni,pn,md)  
%计算极条件闭合差和限差=====================================================
n=size(cpp2,2);
err=[];
maxdiff=[];
for i=1:n
    site=cpp2{i};
    m=size(site,2);
    w=1;
    w1=0;
    for j=1:m-1
        [bet1,k1]=beta(mr(i),site(j),site(j+1));
        [bet2,k2]=beta(site(j),site(j+1),mr(i);
        w=w*sin(bet1)/sin(bet2);
        w1=w1+(cot(bet1))^2+(cot(bet2))^2;
    end
    err(i)=w-1;
    maxdiff(i)=2*sqrt(2*w1)*md;
end
err=err*206265;

fprintf(fit2,'      极条件\n');
fprintf(fit2,'      中心点    环点      闭合差      限差\n');
fprintf(fit2,'=====================================================\n');
for i=1:size(cpp2,2)
    site=cpp2{i};
    if abs(maxdiff(i)-abs(err(i))>=0
        fprntf(fit2,'  %5d %5d  %5d %5d %5d %5d %5d %5d',i,mr(i),site);
        fprintf(fit2,' %10.3f %10.3f\n',err(i),maxdiff(i));
    else
        fprntf(fit2,'* * %5d %5d  %5d %5d %5d %5d %5d %5d',i,mr(i),site);
        fprintf(fit2,' %10.3f %10.3f\n',err(i),maxdiff(i));
    end
end
fprintf(fit2,'=====================================================\n');
fclose(fit2);
open(strcat(pathname,net_name,c2_datafile));
return