资源描述：
《西电数据挖掘大作业k-means和k-medoids》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、题目：数据挖掘学院：电子工程学院专业：智能科学与技术学生姓名：**学号：02115***k-means实验报告一、waveform数据1、算法描述1.从数据集{Xn}n-1N中任意选取k个赋给初始的聚类中心c1,c2,…,ck；2.对数据集中的每个样本点xi，计算其与各个聚类中心cj的欧氏距离并获取其类别标号：labeli=argmin
2、
3、xi-cj
4、
5、2,i=1,…,N,j=1,…,k3.按下式重新计算k个聚类中心；cj=s:labels=jxjNj,j=1,2,…k重复步骤2和步骤3，直到达到最大迭代次数为止2、实验结果一、图像处理1、算法描述同上；2、实验结果代码：k_means:%
6、%%%%%%%%K_means%%%%%%%%%%%%%%%%%函数说明%%%%%%%%%输入:%sample——样本集；%k——聚类数目；%输出：%y——类标(从0开始)%cnew——聚类中心%n——迭代次数function[ycnewn]=k_means(sample,k)[NV]=size(sample);%N为样本的个数K为样本的维数y=zeros(N,1);%记录样本类标dist=zeros(1,k);rand_num=randperm(N);cnew=(sample(rand_num(1,1:k),:));%随机初始化聚类中心cold=zeros(k,V);n=0;while(
7、cold~=cnew)cold=cnew;n=n+1;%记录迭代次数%对样本进行重新分类fori=1:Nforj=1:kif(V==1)dist(1,j)=abs(sample(i,:)-cold(j,:));elsedist(1,j)=norm(sample(i,:)-cold(j,:));endendfors=1:kif(dist(1,s)==min(dist))y(i,1)=s-1;endendend%更新聚类中心cnew=zeros(k,V);flag=zeros(k,1);fori=1:Nforj=1:kif(y(I,1)==j-1)flag(j,1)=flag(j,1)+1;c
8、new(j,J=cnew(j,J+sample(I,J;endendendforj=1:kcnew(j,J=cnew(j,J/flag(j,1);endendk_means_waveform:clear;clc;%%%%%%%%%数据读入%%%%%%%data=load('G:西电2014大三下大作业DataMiningk_meanswaveform.data');[NK]=size(data);%数据集的数目data0=zeros(1,K);data1=zeros(1,K);data2=zeros(1,K);fori=1:Nif(data(i,K)==0)data0=cat(
9、1,data(i,:),data0);elseif(data(i,K)==1)data1=cat(1,data(i,:),data1);elsedata2=cat(1,data(i,:),data2);endendsample=cat(1,data0(1:100,:),data1(1:100,:),data2(1:100,:));label=sample(:,K);%样本的正确类标sample=sample(:,1:K-1);%样本集k=3;%聚类中心的数目%%%%%%%%%K_means%%%%%%%%[ycnewn]=k_means(sample,k);%%%%%%%%%%正确率统计%
10、%%%%%%sum=zeros(1,6);[NV]=size(sample);fori=1:Nif(y(i,1)==label(i,1))sum(1,1)=sum(1,1)+1;endendfori=1:Nif((y(i,1)+label(i,1))==2)sum(1,2)=sum(1,2)+1;endendfori=1:Nif(((y(i,1)==0)&&(label(i,1)==0))
11、
12、((y(i,1)==1)&&label(i,1)==2)
13、
14、((y(i,1)==2)&&label(i,1)==1))sum(1,3)=sum(1,3)+1;endendfori=1:Nif(((y(
15、i,1)==0)&&(label(i,1)==1))
16、
17、((y(i,1)==1)&&label(i,1)==0)
18、
19、((y(i,1)==2)&&label(i,1)==2))sum(1,4)=sum(1,4)+1;endendfori=1:Nif(((y(i,1)==0)&&(label(i,1)==1))
20、
21、((y(i,1)==1)&&label(i,1)==2)
22、
23、((y(i,1)==2)&&label(i,1)==
