模糊c均值聚類 FCM算法的MATLAB代碼 我做畢業論文時需要模糊C-均值聚類，找了好長時間才找到這個，分享給大家：FCM算法的兩種迭代形式的MATLAB代碼寫于下,也許有的同學會用得著:m文件1/7:function U,P,Dist,Cluster_Res,Obj_Fcn,iter=fuzzycm(Data,C,plotflag,M,epsm)% 模糊 C 均值聚類 FCM: 從隨機初始化劃分矩陣開始迭代% U,P,Dist,Cluster_Res,Obj_Fcn,iter = fuzzycm(Data,C,plotflag,M,epsm)% 輸入:% Data: NS 型矩陣,聚類的原始數據,即一組有限的觀測樣本集,% Data 的每一行為一個觀測樣本的特征矢量,S 為特征矢量% 的維數,N 為樣本點的個數% C: 聚類數,1CN% plotflag: 聚類結果 2D/3D 繪圖標記,0 表示不繪圖,為缺省值 % M: 加權指數,缺省值為 2% epsm: FCM 算法的迭代停止閾值,缺省值為 1.0e-6% 輸出:% U: CN 型矩陣,FCM 的劃分矩陣% P: CS 型矩陣,FCM 的聚類中心,每一行對應一個聚類原型% Dist: CN 型矩陣,FCM 各聚類中心到各樣本點的距離,聚類中% 心 i 到樣本點 j 的距離為 Dist(i,j)% Cluster_Res: 聚類結果,共 C 行,每一行對應一類% Obj_Fcn: 目標函數值% iter: FCM 算法迭代次數% See also: fuzzydist maxrowf fcmplotif nargin5 epsm=1.0e-6; endif nargin4 M=2;endif nargin4 | plotflag Obj_Fcn(iter)=sum(sum(Um.*Dist.2); end % FCM 算法迭代停止條件 if norm(U-U0,Inf) 3 res = maxrowf(U); for c = 1:C v = find(res=c); Cluster_Res(c,1:length(v)=v; endend% 繪圖if plotflag fcmplot(Data,U,P,Obj_Fcn);endm文件2/7:function U,P,Dist,Cluster_Res,Obj_Fcn,iter=fuzzycm2(Data,P0,plotflag,M,epsm)% 模糊 C 均值聚類 FCM: 從指定初始聚類中心開始迭代% U,P,Dist,Cluster_Res,Obj_Fcn,iter = fuzzycm2(Data,P0,plotflag,M,epsm)% 輸入: Data,plotflag,M,epsm: 見 fuzzycm.m% P0: 初始聚類中心% 輸出: U,P,Dist,Cluster_Res,Obj_Fcn,iter: 見 fuzzycm.m % See also: fuzzycmif nargin5 epsm=1.0e-6; endif nargin4 M=2;endif nargin4 | plotflag Obj_Fcn(iter)=sum(sum(Um.*Dist.2); end % FCM 算法迭代停止條件 if norm(P-P0,Inf) 3 res = maxrowf(U); for c = 1:C v = find(res=c); Cluster_Res(c,1:length(v)=v; endend% 繪圖if plotflag fcmplot(Data,U,P,Obj_Fcn);endm文件3/7:function fcmplot(Data,U,P,Obj_Fcn)% FCM 結果繪圖函數% See also: fuzzycm maxrowf ellipseC,S = size(P); res = maxrowf(U);str = po*x+dv2 figure(2),plot3(P(:,1),P(:,2),P(:,3),rs),hold on for i=1:C v=Data(find(res=i),:); plot3(v(:,1),v(:,2),v(:,3),str(rem(i,12)+1) ellipse(max(v(:,1)-min(v(:,1), . max(v(:,2)-min(v(:,2), . max(v(:,1)+min(v(:,1), . max(v(:,2)+min(v(:,2)/2, . r:,(max(v(:,3)+min(v(:,3)/2) end grid on,title(3D 聚類結果圖,fontsize,8),hold offendm文件4/7:function D=fuzzydist(A,B)% 模糊聚類分析: 樣本間的距離% D = fuzzydist(A,B)D=norm(A-B);m文件5/7:function mr=maxrowf(U,c)% 求矩陣 U 每列第 c 大元素所在行,c 的缺省值為 1% 調用格式: mr = maxrowf(U,c)% See also: addrif nargin2 c=1;endN=size(U,2);mr(1,N)=0;for j=1:N aj=addr(U(:,j),descend); mr(j)=aj(c);endm文件6/7:function ellipse(a,b,center,style,c_3d)% 繪制一個橢圓% 調用: ellipse(a,b,center,style,c_3d)% 輸入:% a: 橢圓的軸長(平行于 x 軸)% b: 橢圓的軸長(平行于 y 軸)% center: 橢圓的中心 x0,y0,缺省值為 0,0% style: 繪制的線型和顏色,缺省值為實線藍色% c_3d: 橢圓的中心在 3D 空間中的 z 軸坐標,可缺省if nargin4 style=b;endif nargin4 plot3(x,y,ones(1,360)*c_3d,style)else plot(x,y,style)endm文件7/7:function f = addr(a,strsort)% 返回向量升序或降序排列后各分量在原始向量中的索引% 函數調用:f = addr(a,strsort)% strsort: ascend or descend% default is ascend% - example -% addr( 4 5 1 2 ) returns ans:% 3 4 1 2 if nargin=1 strsort=ascend;endsa=sort(a); ca=a;la=length(a);f(la)=0;for i=1:la f(i)=find(ca=sa(i),1); ca(f(i)=NaN;endif strcmp(strsort,descend) f=fliplr(f);end幾天前我還在這里發帖求助，可是很幸運在其他地方找到了，在這里和大家分享一下！function center, U, obj_fcn = FCMClust(data, cluster_n, options) % FCMClust.m 采用模糊C均值對數據集data聚為cluster_n類% % 用法： % 1.center,U,obj_fcn = FCMClust(Data,N_cluster,options); % 2.center,U,obj_fcn = FCMClust(Data,N_cluster); % % 輸入： % data - nxm矩陣,表示n個樣本,每個樣本具有m的維特征值 % N_cluster - 標量,表示聚合中心數目,即類別數 % options - 4x1矩陣，其中 % options(1):隸屬度矩陣U的指數，1 (缺省值: 2.0) % options(2):最大迭代次數 (缺省值: 100) % options(3):隸屬度最小變化量,迭代終止條件 (缺省值: 1e-5) % options(4):每次迭代是否輸出信息標志 (缺省值: 1) % 輸出： % center - 聚類中心 % U - 隸屬度矩陣 % obj_fcn - 目標函數值 % Example: % data = rand(100,2); % center,U,obj_fcn = FCMClust(data,2); % plot(data(:,1), data(:,2),o); % hold on; % maxU = max(U); % index1 = find(U(1,:) = maxU); % index2 = find(U(2,:) = maxU); % line(data(index1,1),data(index1,2),marker,*,color,g); % line(data(index2,1),data(index2,2),marker,*,color,r); % plot(center(1 2,1),center(1 2,2),*,color,k) % hold off; if nargin = 2 & nargin = 3, %判斷輸入參數個數只能是2個或3個 error(Too many or too few input arguments!); end data_n = size(data, 1); % 求出data的第一維(rows)數,即樣本個數 in_n = size(data, 2); % 求出data的第二維(columns)數，即特征值長度 % 默認操作參數 default_options = 2; % 隸屬度矩陣U的指數 100; % 最大迭代次數 1e-5; % 隸屬度最小變化量,迭代終止條件 1; % 每次迭代是否輸出信息標志if nargin = 2, options = default_options; else %分析有options做參數時候的情況 % 如果輸入參數個數是二那么就調用默認的option; if length(options) 4, %如果用戶給的opition數少于4個那么其他用默認值; tmp = default_options; tmp(1:length(options) = options; options = tmp; end % 返回options中是數的值為0(如NaN),不是數時為1 nan_index = find(isnan(options