muoshao.asv

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

function qe=muoshao(z,m,b,t)
m=[0 1 1 0 0 1 0 1 0 0 1 0 0]
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一个具有两个连通系的信息矩阵
b=[0 0 0 0 1 0 0 0 0 1 0 1 0];
t=[0 1 1 0 1 1 0 1 0 1 1 1 0]
d=qiu_d(z)
c=dc(d,m,t)
[c1,c2]=size(c)
qb1=[]
qe1=[]
tb=[]
muoshao1=[]
muoshao2=[]
[m1,m2]=find(m==1)
[m3,m4]=size(m2)
w=1
for i=1:m4
    qb1(1)=m2(i)
    tb(1)=m2(i)
    while(length(qb1)~=0)
        a=qb1(1)
        qb1(1)=[]
        for j=1:c2
            if(c(a,j)==1)
                qb1(length(qb1)+1)=j
                tb(length(tb)+1)=j
             end
                if(c(a,:)==0)
                qe1(length(qe1)+1)=a
                break
            end
        end
    end
    muoshao1=[m2(i) 0 qe1 0 tb]   %电源点,末稍点,树节点队列分别用0隔开
    qe1=[]
    qb1=[]
    tb=[]
    [b1,b2]=size(muoshao1)
        for s=1:b2
            muoshao2(w,s)=muoshao1(1,s)
        end
   w=w+1
end
muoshao2
end