谱聚类算法MATLAB

function [IDX,C,sumd,D,U,perm1,A100] = spectrul(A, num_clusters) %最后的聚类程序 A100=[]; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% n = size(A, 1); %[U11, S11, V11] = svd(A); %svds = S11(1:num_clusters, 1:num_clusters); %s11 = diag(S11); %scatter(1:length(s11), s11) %hold on %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% n=size(A,1); sum1=0; sum2=0; sum31=0; sum32=0; S1=[]; S2=[]; for i=1:n-1 for j=i+1:n for k=1:n sum1=sum1+(A(i,k)+A(j,k)); if A(i,k)>0 && A(j,k)>0 sum2=sum2+(A(i,k)+A(j,k)); end end S1(i,j)=exp(sum2/sum1); sum1=0; sum2=0; end end for i=2:n for j=1:i-1 S1(i,j)=S1(j,i); end end for i=1:n S1(i,i)=1; end S1 sums1=0; sums2=0; Ss1=[]; Ss2=[]; Ss3=[]; for j=1:n-1 for k=j+1:n for i=1:n sums1=sums1+(A(i,j)+A(i,k)); if A(i,j)>0 && A(i,k)>0 sums2=sums2+(A(i,j)+A(i,k)); end end Ss1(j,k)=exp(sums2/sums1); sums1=0; sums2=0; end end for i=2:n for j=1:i-1 Ss1(i,j)=Ss1(j,i); end end for i=1:n Ss1(i,i)=1; end Ss1 for i=1:n for j=1:n S0(i,j)=(S1(i,j)+Ss1(i,j))/2; end end S0 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Do laplacian, L = D^(-1/2) * S * D^(-1/2) % disp('Doing Laplacian...'); D1 = sum(S0, 2) + (1e-10); D1 = sqrt(1./D1); % D^(-1/2) D1 = spdiags(D1, 0, n, n); L = D1 * S0 * D1; [U111, S111, V111] = svd(L); svds1 = S111(1:num_clusters, 1:num_clusters); s111 = diag(S111); scatter(1:length(s111), s111) hold on % % Do eigendecomposition, if L = % D^(-1/2) * S * D(-1/2) : set 'LM' (Largest Magnitude), or % I - D^(-1/2) * S * D(-1/2): set 'SM' (Smallest Magnitude). % disp('Performing eigendecomposition...'); OPTS.disp = 0; [V, val] = eigs(L, num_clusters, 'LM', OPTS); % % Do k-means % disp('Performing kmeans...'); % Normalize each row to be of unit length sq_sum = sqrt(sum(V.*V, 2)) + 1e-20; U = V ./ repmat(sq_sum, 1, num_clusters); [IDX,C,sumd,D] = kmeans(U,num_clusters,'distance','cosine'); U; %算法树 figure Y = pdist(U,'cosine'); Z = linkage(Y,'average'); [H,T,perm] = dendrogram(Z,0); perm; perm1=perm'%算法树的次序 set(H,'LineWidth',2) hold on % % Calculate and show time statistics % V_sequ=perm; vlength=length(V_sequ); A100=[]; for xx=1:vlength %column exchange A_M(:,xx)=A(:,V_sequ(xx)); end for xx=1:vlength % row exchange A100(xx,:)=A_M(V_sequ(xx),:); end A100 disp('Finished!');