clc;
clear;
aa=[];
t=1;
aa(:,:,t)=xlsread('Topologyidentification.xlsx',1);%网络结构
AS = aa(:,:,t);
B = size(AS,1);
for i = 1:B
for j = 1:B
if (aa(i,j,t)==1)&&(i<j)
AS(i,j)=0;
AS(j,i)=0;
AS(j,j)=0;
end
end
end
BS = [];%孤岛节点
n = 0;%孤岛数
for i=1:B
if AS(i,i)==1
n=n+1;
c=1;
BS(n,c)=i;
for j=i+1:B
if aa(i,j,t)==1
c=c+1;
BS(n,c)=j;
end
end
for j=i+1:B
if (AS(j,j,t)==0)&&(aa(i,j,t)==0)
for k=i+1:B
if (aa(k,j,t)==1)&&(aa(i,k,t)==1)
c=c+1;
BS(n,c)=j;
end
end
end
end
end
end