muoshao1.m

用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')
qe1