muoshao1.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)
qe1=[]
tb=[]
qb1=[]
[m1,m2]=find(m==1)
[m3,m4]=size(m2)
 fprintf('请输入一个源点:\n')
 mm=input('')
sign=1
while sign
        [q1,q2]=find(m2==mm)
        if length(q2)==0
        fprintf('对不起您输入的不是源点,请重新输入:\n')
        mm=input('')
        else
        break
        end
end
qb1(1)=mm
while(length(qb1)~=0)
        a=qb1(1)
        qb1(1)=[]
        for j=1:c2
            if(c(a,j)==1)
                qb1=[qb1 j]
                tb=[tb j]
             end
                if(c(a,:)==0)
                qe1=[qe1 a]
                break
            end
        end
end
fprintf('此源点的末稍点是:\n')