dc.m

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

function c=dc(d,m,t)
[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