k-shell分解算法

%function [ret_m]=k_shell(mixedsig)
%邻接矩阵
%  A=load('cdbBA_4000_5_0_.txt');   %里边的数据结构好像是每2个一组，例如 3 4 表示一组，表示点3和节点4有连接。
%  
%  
%  TT=A(:, 1:2)+1;
%  maxtt=max(max(TT));
%  mixedsig=zeros(max(max(TT)));
%  len=length(TT);
%  for i=1:len
%      mixedsig(TT(i,1),TT(i,2))=1;
%      mixedsig(TT(i,2),TT(i,1))=1;
%  end
%  disp('已转化为邻接矩阵！')

[XY]=distri(mixedsig);
mixedsig=[0,1,1,1;1,0,0,1;1,0,0,1;1,1,1,0]
function [ret_m]=kshell(mixedsig)
mixedsig=[0,1,0,1;1,0,0,1;0,0,0,1;1,1,1,0];
mixedsig=[0,1,0,1,0;1,0,0,1,0;0,0,0,1,1;1,1,1,0,0;0,0,1,0,0];

zzz=0;
j=1;   %第j层
t=1;  %
tad=mixedsig;  %邻接矩阵
k=0;
ii=1;
ret_m=[];  %返回值；  i，j  ： 第i层：所有节点
len=length(mixedsig);  %矩阵长度
tt=1;