poleterm.m

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

function [cpp2,mr,tr]=poleterm(sd,e,d,sid,g,f,dir,si,ni,pn,md)   %应用余枝搜索法自动搜索所有极条件   
global pathname net_name c2_datafile
fit2=fopen(strcat(pathname,net_name,c2_datafile),'wt');
cpp2=[];
n=0;
tr=[];
mr=[];
for i=1:sd
    num=0;
    cpp=[];temp4=[];
    temp4=find(i==f);
    xx=[]; %存放外围点
    for k=temp4
        xx=[xx g(k)];
    end
    num=lengeh(xx);
    if num>2
        numb=0;
        d1=[];d2=[];ds=[];
        for i0=1:num-1
            for j0=i0+1:num
                [dir1,kdir1]=finddir(xx(i0),xx(j0),g,f,dir,si,ni); %是否存在1-2方向观测
                [dir2,kdir2]=finddir(xx(j0),xx(i0),g,f,dir,si,ni); %是否存在2-1方向观测
                if kdir1+kdir2==2
                    numb=numb+1;
                    d1(numb)=xx(i0);
                    d2(numb)=xx(j0);
                    ds(numb)=1000;
                end
            end
        end
        dd1=zero(1,length(d1));
        dd2=zero(1,length(d1));
        for y=1:num %对外围点重新从1-num编号
            for y1=1:length(d1)
                if d1(y1)==xx(y)
                    dd1(y1)=y;
                end
                if d2(y1)==xx(y)
                    dd2(y1)=y;
                end
            end
        end
        cpp=ccpoints(dd1',dd2',num,length(dd1),ds); %寻找外围点闭合环
    end
    if isempty(cpp)~=1
        for t=1:size(cpp,2)
            n=n+1;
            site=cpp{t};
            cpp2{n}=xx(site);
            mr=[mr i];
        end
        tr=[tr num];
    end
end