1、22、實驗結(jié)果題目：數(shù)據(jù)挖掘?qū)W院：電子工程學院專業(yè)：智能科學與技術(shù)學生姓名：*學號：02115*k-means實驗報告waveformwaveform1、算法描述從數(shù)據(jù)集XN中任意選取k個賦給初始的聚類中心c1,nn-11c2,ck；對數(shù)據(jù)集中的每個樣本點x計算其與各個聚類中心&的歐氏距離并獲取其類別標號：label(i)=argmin|x-c.|2,i=1,，N,j=1,，kij按下式重新計算k個聚類中心；Vxc=s：label(s)=j=Nj重復步驟2和步驟3，直到達到最大迭代次數(shù)為止22、實驗結(jié)果匚QimmandWinch訓-hX慈中宙3Columns1through1D-Ck1D0S0

2、L26710.52410.317Q0/91101.35232.D7dP2.217d2.15772.5199O.fllSO0.43651.G3D61.542D2.2561.54594,D5223.3257302523.8155-0.1837(X137fl0.0787-0.DD02D.0694D.16010.6.271(k43376361.4M&oll23E：S11through202.97773.04343.04783.27143.14072.53491.79931.3905085880.63fiB32132.D2S11.El61.1122Duf21QO.0130.1M1-Ou22S10.162

3、10.D7E72.42263.130G3.761d.B15L5.36CiSi.36613.52262.2376L”75130.8695Colvin210.0682-0.0915-0.02510.6367XT1、算法描述同上；回FileEditViewInsertToolsDesktopWindowHelp回FileEditViewInsertToolsDesktopWindowHelpcia篷既紗晏凰日匡I原團像換度團像二蓿化圖像代碼：k_means:%K_means%函數(shù)說明%輸入:%sample樣本集；%k聚類數(shù)目；%輸出：%y類標（從0開始）%cnew聚類中心%n迭代次數(shù)function

4、ycnewn=k_means(sample,k)NV=size(sample);%N為樣本的個數(shù)K為樣本的維數(shù)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(cold=cnew)cold=cnew;n=n+1;%記錄迭代次數(shù)%對樣本進行重新分類fori=1:Nforj=1:kif(V=1)dist(1,j)=abs(sample(i,:)-cold(j,:);elsedist(1,j)=norm

5、(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;cnew(j,=cnew(j,+sample(I,;endendendforj=1:kcnew(j,=cnew(j,/flag(j,1);endendk_means_waveform:clear;clc;%數(shù)據(jù)讀入%data=load(G:西電2014大三下大

6、作業(yè)DataMiningk_meanswaveform.data);NK=size(data);%數(shù)據(jù)集的數(shù)目data0=zeros(1,K);data1=zeros(1,K);data2=zeros(1,K);fori=1:Nif(data(i,K)=0)data0=cat(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,

7、:);label=sample(:,K);%樣本的正確類標sample=sample(:,1:K-1);%樣本集k=3;%聚類中心的數(shù)目%K_means%ycnewn=k_means(sample,k);%正確率統(tǒng)計%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)|(y(i,1)=1)&la

8、bel(i,1)=2)|(y(i,1)=2)&label(i,1)=1)sum(1,3)=sum(1,3)+1;endendfori=1:Nif(y(i,1)=0)&(label(i,1)=1)|(y(i,1)=1)&label(i,1)=0)|(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)|(y(i,1)=1)&label(i,1)=2)|(y(i,1)=2)&label(i,1)=0)sum(1,5)=sum(1,5)+1;endendfori=1:Nif(y(i,1

9、)=0)&(label(i,1)=2)|(y(i,1)=1)&label(i,1)=0)|(y(i,1)=2)&label(i,1)=1)sum(1,6)=sum(1,6)+1;endendsum=sum/N;creatrate=max(sum);disp（循環(huán)次數(shù):）;disp(n);disp(聚類中心為:);disp(cnew);disp(正確率為:);disp(creatrate);k_means_picture:clear;clc;%數(shù)據(jù)讀入%口=imread(G:西電2014大三下大作業(yè)DataMiningk_meanslena.jpg);I2=rgb2gray(I1);%轉(zhuǎn)化為灰度

10、圖像I=im2double(I2);numv=size(I);sample=reshape(I,v*num,1);%樣本集k=2;%聚類中心的數(shù)目%K_means%ycnewn=k_means(sample,k);%v%I3=sample;if(cnew(1,1)=cnew(2,1)F0=255;F1=0;elseF0=0;F1=255;endfori=1:num*vif(y(i,1)=0)I3(i,1)=F0;elseI3(i,1)=F1;endendI3=reshape(I3,num,v);figure(1)subplot(1,3,1);imshow(I1);title(原圖像);subp

11、lot(1,3,2);imshow(I2);title(灰度圖像);subplot(1,3,3);imshow(I3);title(二值化圖像);k_medoids實驗報告waveformwaveform1、算法描述(1)隨機選擇k個對象作為初始的代表對象;(2)repeat(3)指派每個剩余的對象給離它最近的代表對象所代表的簇；隨意地選擇一個非代表對象Orandom;計算用Orandom代替Oj的總代價S；如果S0.001)%將所有樣本分配到最近的代表點fori=1:Nforj=1:kdist(1,j)=norm(sample(i,:)-med(j,:);endtemp(i,1)=min(d

12、ist);fors=1:kif(dist(1,s)=temp(i,1)temp(i,2)=s;endendendy=temp(:,2);Eold=sum(temp(:,1);%隨機的選擇一個非代表點，生成新的代表點集合index=index+1;med_temp=med;E=zeros(1,k);forj=1:kmed_temp(j,:)=sample(index,:);%將所有樣本分配到最近的代表點fori=1:Nfort=1:kdist(1,t)=norm(sample(i,:)-med_temp(t,:);endtemp(i,1)=min(dist);fors=1:kif(dist(1,

13、s)=temp(i,1)temp(i,2)=s;6uiuiiai6uiuiiai刃旳亦為丄三枕0乙曲厘：9.)peo|=eiep%那11人心%：DpJeap廠斗puapuaP|O3=M9U3as|a：(:/xapui)a|dujes=(:/isaqo)pauj(p|O3Mau3)ipuapua：l=isaqo(Mau3=(n)3)!|：l=lJO(3)uiiu=Mau3pua：(-,:)diuai)iuns=(PT)3puapuapuak_medoidswaveform-+noise.data);NK=size(data);%數(shù)據(jù)集的數(shù)目data0=zeros(1,K);data1=zeros

14、(1,K);data2=zeros(1,K);fori=1:Nif(data(i,K)=0)data0=cat(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,:);Iabel=sample(:,K);%Nu%pAOyEAaesample=sample(:,1:K-1);%NU士場一k=3;%聚類中心的數(shù)目%k_medoid

15、s%ymed=k_medoids(sample,k);%正確率統(tǒng)計圧%sum=zeros(1,6);NV=size(sample);fori=1:Ny(i,1)=y(i,1)-1;endfori=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)|(y(i,1)=1)&label(i,1)=2)|(y(i,1)=2)&label(i,1)=1)sum(1,3)=su

16、m(1,3)+1;endendfori=1:Nif(y(i,1)=0)&(label(i,1)=1)|(y(i,1)=1)&label(i,1)=0)|(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)|(y(i,1)=1)&label(i,1)=2)|(y(i,1)=2)&label(i,1)=0)sum(1,5)=sum(1,5)+1;endendfori=1:Nif(y(i,1)=0)&(label(i,1)=2)|(y(i,1)=1)&label(i,1)=0)|(y

17、(i,1)=2)&label(i,1)=1)sum(1,6)=sum(1,6)+1;endendsum=sum/N;creatrate=max(sum);disp(-$0A2ooo;a。)disp(med)jdisp(oyLAEia。)disp(creafrat:e)jkmedoidspicr+ure-c-earnp%Ey$WAEe%loll.lnread(-G,0(ffiol4汁川-Hv?f-TTrf=XDaarMiningklmedoidsx-ena.jpg-)DuiljIlHimnoiseuo-gaussiarTOD)!12229目(11)旅左護對畫鄭Ililn2doub-e(12)j-FluinvTsizessainp-eHreshape(LV*num、1)j%NbakHz%$0A2oooAHAEyA、%Kmeans%-ylnedTkllnedoids(salnp-pk)j%畫香逵%I3usalnppif(med(l、lvHITIe