untitled3.m

来自「这是一个模糊聚类代码（matlab）供大家学习」· M 代码 · 共 22 行

%程序 L10_2.m: 模糊聚类%
%模糊 C-平均(FCM)算法%
%N：输入样本数，dimen：样本维度%
%X：输入样本(N*dimen)%
%cluster：群数， m：指数型权重(1<r)%
%tolerance：误差临界值%
%U：分割矩阵(大小为cluster*N)%
clear; figure;
X=rand(200,2);
[N,dimen]=size(X);
plot(X(:,1),X(:,2),'*'), axis([-0.1 1.1 -0.1 1.1])
hold on
%步骤一：设定一些初值及初始分割矩阵%
cluster=4; r=2; tolerance=1e-5; cycle=500;
U=rand(cluster,N);
ss=sum(U);
U_new=U./ss(ones(cluster,1),:);    %初始模糊分割矩阵
for L=1:cycle
%步骤二：计算模糊聚类中心%
 U=U_new;
 Uexp=U.^r;
 C=Uexp*X./( (ones(dimen,1)*sum(Uexp'))' );    %新群心的值