求，关于聚类分析方面的问题？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

　大家好！我现在做的是聚类分析，高手老师们能不能给点matlab关于聚类的原代码？多谢！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：聚类分析 MATLAB atlab matla 大家好 聚类分析

这个能弄到吗？

function y=kMeansCluster(m,k,isRand)

if nargin<3,        isRand=0;   end
if nargin<2,        k=1;        end   
[maxRow, maxCol]=size(m)
if maxRow<=k,     
     y=[m, 1:maxRow]
else  
  % initial value of centroid
    if isRand,
       p = randperm(size(m,1));      % random initialization
        for i=1:k
            c(i,:)=m(p(i),:)
       end
    else
        for i=1:k
           c(i,:)=m(i,:)        % sequential initialization
     end
    end
     temp=zeros(maxRow,1);   % initialize as zero vector
     while 1,
        d=DistMatrix(m,c);  % calculate objcets-centroid distances
        [z,g]=min(d,[],2);  % find group matrix g
        if g==temp,
            break;          % stop the iteration
        else
            temp=g;         % copy group matrix to temporary variable
        end
        for i=1:k
            f=find(g==i);
            if f            % only compute centroid if f is not empty
                c(i,:)=mean(m(find(g==i),:),1);
            end
        end
end
     y=[m,g];
end

function d=DistMatrix(A,B)
                     
             [hA,wA]=size(A);
             [hB,wB]=size(B);
             if wA ~= wB, error(' second dimension of A and B must be the same'); end
             for k=1:wA
                 C{k}= repmat(A(:,k),1,hB);
                 D{k}= repmat(B(:,k),1,hA);
             end
             S=zeros(hA,hB);
             for k=1:wA
                 S=S+(C{k}-D{k}').^2;
             end
             d=sqrt(S);

这是k—Means聚类算法

我会用NTSYS和SPSS做，方便的话发到QQ，我帮你处理吧

感觉你是求别人想你提问题的！

谢谢大家的帮助！！