niantongxi.asv

用MATLAB编的配电网自动化的一些程序

function niantongxi(z,m)
z=[6 10;10 13;13 8;10 1;1 5;5 11;5 9;9 2;3 12;12 7;12 4]   %p12图2-2一个具有两个连通系的信息矩阵
m=[0 1 1 0 0 1 0 1 0 0 1 0 0]
d=qiu_d(z)
[m1,m2]=size(m)
vs=zeros(m1,m2)
qc=[]
n=0
[m3,mm]=find(m==1)
[m5,m6]=size(mm)
for i=1:m6
    qc(1)=mm(i)
    n=n+1
    [e,f]=find(vs==0)
    [g,h]=size(f)
    if (h==0)                  %当vs的每个元素均不为0时结束循环
        break
    end
    while(length(qc)~=0)
        a=qc(1)
        qc(1)=[]
        vs(a)=n
         for j=1:m2
             if(d(a,j)==1&vs(j)~=n)
                    qc(length(qc)+1)=j
             end
         end
     end
end
aa=max(vs)
bb=length(vs)
for i=1:aa
    gg=1
    for j=1:bb
        if (vs(j)==i)
            s{i}(gg)=j
            gg=gg+1
        end
    end
end
for i=1:aa
    fprintf('此网')
    disp(s{i})
end