dc.m

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

function c=dc(d,m,t)
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]
t=[0 1 1 0 1 1 0 1 0 1 1 1 0]
d=qiu_d(z)
[b1,b2]=size(d)
c=zeros(b1,b2)
qs=[]
[i1 j]=find(m==1)
[j1 j2]=size(j)
for i=1:j2
    qs(i)=j(i)
end
while(length(qs)~=0)
    a=qs(1)
    qs(1)=[]
    for i2=1:b2
        if (d(a,i2)==1&t(i2)==0)
            d(a,i2)=1,d(i2,a)=0
        elseif(d(a,i2)==1&t(i2)==1&d(i2,a)==1)
            d(a,i2)=1,d(i2,a)=0
            qs=[i2 qs]
        end
    end
    
end
c=d