k_means聚类

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

学习第一篇

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：means mean ans

wenjian=xlsread('wenjian.xlsx');
[m1,n1]=size(wenjian);
%初始中心点
a1=3;
a2=5;
b1=8;
b2=6;
d1=zeros(m1,1);
d2=zeros(m1,1);
biaozhun1=10;%初始运行条件
biaozhun2=10;
while biaozhun1>0.00000001 & biaozhun2>0.0000001%设置终止条件
for i=1:m1
    d1(i)=sqrt((wenjian(i,1)-a1)^2+(wenjian(i,2)-a2)^2);
end
for i=1:m1
    d2(i)=sqrt((wenjian(i,1)-b1)^2+(wenjian(i,2)-b2)^2);
end
%保存分类的数据
for i=1:m1
    if d1(i)-d2(i)<=0
    s1(i)=i;
    else
    s2(i)=i;
    end
end
s1(s1==0)=[];
s2(s2==0)=[];

n1=length(s1);%计算新的聚类中心
n2=length(s2);
sum11=sum(wenjian(s1,1));
sum12=sum(wenjian(s1,2));
mean11=sum11/n1;
mean12=sum12/n1;
c1=[mean11,mean12];

sum21=sum(wenjian(s2,1));
sum22=sum(wenjian(s2,2));
mean21=sum21/n2;
mean22=sum22/n2;
c2=[mean21,mean22];
biaozhun1=sqrt((mean11-a1)^2+(mean12-a2)^2);
biaozhun2=sqrt((mean21-b1)^2+(mean22-b2)^2);

a1=mean11;%更新聚类中心点
a2=mean12;
b1=mean21;
b2=mean22;

end
center1=[a1,a2]
center2=[b1,b2]

x=wenjian(:,1);
y=wenjian(:,2);
plot(x,y,'o',a1,a2,'b*',b1,b2,'r*')
%小结 距离公式必须加上根号，否则偏差特别大

1. yuan=xlsread('yuan.xlsx',1,'a2:c1001');
   
2. [m,n]=size(yuan);
   
3. x=sum(yuan(:,2));
   
4. y=sum(yuan(:,3));
   
5. xbar1=x/m;
   
6. ybar1=y/m;
   
7. temp=zeros(m,1);
   
8. d1=zeros(m,1);
   
9. d2=zeros(m,1);
   
10. s1=zeros(m,1);
    
11. s2=zeros(m,1);
    
12. %% 设置初始值
    
13. for i=1:m
    
14.     temp(i)=sqrt((yuan(i,2)-xbar1)^2+(yuan(i,3)-ybar1)^2);
    
15. end
    
16. c1=min(temp);
    
17. c2=max(temp);
    
18. biaozhun1=10;%初始运行条件
    
19. biaozhun2=10;

复制代码