yj2.m

余弦距离进行聚类分析

% q 26*54
load q
load t

b(:,1)=0*q(:,1);
b(:,2)=0.2*q(:,37);
b(:,3)=0.2*q(:,39);
b(:,4)=0.2*q(:,40);
b(:,5)=0.2*q(:,41);
b(:,6)=0.2*q(:,42);

[i,c,s,d]=kmeans(b,3,'Distance','cosine');

m(:,1)=0.01*c(:,1);
m(:,2)=0.01*c(:,2);
m(:,3)=0.01*c(:,3);
m(:,4)=0.01*c(:,4);
m(:,5)=0.01*c(:,5);
m(:,6)=0.01*c(:,6);




for ii=1:3
xx(ii)=odist(t,m(ii,:));
end
[temp,idx]=min(xx);
ii=1;
for jj=1:length(i)
    if i(jj)==idx
        tt(ii,:)=b(jj,:);
        dist(ii)=odist(t,tt(ii,:));
        ii=ii+1;
        
    end
end

for ii=1:3
   [temp,resultidx(ii)]=min(dist);
   dist(resultidx(ii))=9999999;
   result(ii,:)=q(resultidx(ii),:);
end