Bchengxu.zip_The Program_recordyba_费

function [links,H]=fenlei(jianhuidu) H=zeros; for k=1:length(jianhuidu) temp=find(jianhuidu{k}==1); hang=0; j=0; for i=2:length(temp) j=j+1; if temp(i)-temp(i-1)~=1 hang=hang+1; if mod(j,2) H(k,hang)=i+1-(j+1)/2; else H(k,hang)=i+1-j/2; end j=0; end end end for i=1:length(H(1,:)) for j=1:length(H(:,1)) if H(j,i)==0 H(j,i)=nan; end end end %% for i=1:length(jianhuidu) for j=1:length(jianhuidu) if sum(ismember(H(j,:),H(i,:)))==length(find(H(j,:)>=0)) temp(i,j)=1; else temp(i,j)=0; end end end [a,b]=find(temp==1); links=[a';b']; list=zeros; leishu=0; tem={}; list={}; for i=1:209 tem{i}=[(find(temp(i,:)==1)),(find(temp(:,i)==1))]; end list{1}=tem{1}; for i=1:209 flag=1; for j=1:209 if intersect(list{i},tem{j}) list{i}=union(list{i},tem{j}); flag=0; end end if flag i=i+1; end end end % leishu=0; % list=zeros; % for i=1:length(jianhuidu) % if ~any(list==i) % leishu=leishu+1; % list(leishu,1)=i; % list(leishu,2)=0.5; % end % for j=1:length(jianhuidu) % if i~=j % if sum(ismember(H(j,:),H(i,:)))==length(find(H(j,:)>=0)) % if ~any(list==j) % [a,b]=find(list==i); % [c,d]=find(list(a,:)==0.5); % list(c,d)=j; % list(c,d+1)=0.5; % end % end % if sum(ismember(H(i,:),H(j,:)))==length(find(H(i,:)>=0)) % if~any(list==j) % [a,b]=find(list==i); % [c,d]=find(list(a,:)==0.5); % list(c,d)=j; % list(c,d+1)=0.5; % end % end % end % end % % if flag % % leishu=leishu+1; % % list(leishu,1)=i; % % list(leishu,2)=0.5; % % end % end % % % % % % if ~any(list==i) % % if ~any(list==j) % % leishu=leishu+1; % % list(leishu,1)=j; % % list(leishu,2)=i; % % list(leishu,3)=0.5; % % else [a,b]=find(list==j); % % list(find(list(a,:))==0.5)=i; % % list(a,b+1)=0.5; % % end % % end % for i=1:length(list(1,:)) % for j=1:length(list(:,1)) % if list(j,i)==0.5 % list(j,i)=0; % end % end % end % jianyan=[1:209]; % flag=sum(sum(ismember(jianyan,list))); % end