yuandian.m

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

function qe=yuandian(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);
qb2=[];
qe2=[];
for i=1:c1
    if (c(i,:)==0)
        qe2(length(qe2)+1)=i;
    end
end
a=length(qe2);
fprintf('请输入末稍点:\n')
v=input('');
sign=1;
while sign
    [b1,b2]=find(qe2==v);
    if (length(b2)==0)
        fprintf('输入的节点不是末稍点!\n请重新输入:\n')
        v=input('')
    else
        break
    end
end
q=[];
q(1)=v;
while (length(q)~=0)
       v=q(1);
       q(1)=[];
       for i=1:c1
           if (c(i,v)==1)
               q(length(q)+1)=i;
           elseif(c(:,v)==0)
              qb2(length(qb2)+1)=v;
           break 
          end
       end
end
fprintf('该末稍点的源点是:\n')
qb2