wabc.m

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

function wabc
global pathname net_name c1_datafile
global ed dd g f pn ni si m2 md dir
fit2=fopen(strcat(pathname,net_name,c1_datafile),'wt');
if(fit2==-1)
    msgbox('Error by Opening Output File','Warning','warn');
    return;
end 
fprintf(fit2,'  pn1   pn2   pn3     closure\n');
fprintf(fit2,'========================================================\n');
for i=1:ed+dd-2
    j1=i;
    for j=si(i):si(i)+ni(i)-2
        j2=f(j);
        for k=j+1:si(i)+ni(i)-1
            j3=f(k);
            if(j2>j1&j3>j1)
                [bet2,k2]=beta(j1,j2,j3);
                [bet3,k3]=beta(j2,j3,j1);
                if(k2==1&k3==1)
                    bet1=dir(k)-dir(j);
                    if(bet1<0)
                        bet1=2.*pi+bet1;
                    end
                    if(bet1>pi)
                        bet1=2.*pi-bet1;
                    end
                    ww=(bet1+bet2+bet3-pi)*206265.;
                    fprintf(fit2,'%5d %5d %5d %10.3f\n',pn(j1),pn(j2),pn(j3),ww);
                end
            end
        end
    end
end
fprintf(fit2,'========================================================\n');
fclose(fit2);
open(strcat(pathname,net_name,c1_datafile));
return