FCM聚类matlab代码

下面是一个简单的FCM聚类的MATLAB代码：

% 生成随机数据
data = rand(100,2);

% 设置聚类数目
c = 3;

% 设置模糊因子m
m = 2;

% 初始化隶属度矩阵U
U = rand(size(data,1),c);

% 迭代次数
max_iter = 100;

% 迭代过程
for i = 1:max_iter
    % 计算聚类中心
    centroids = U'*data./sum(U)';
    
    % 计算距离矩阵
    dist = pdist2(data,centroids);
    
    % 计算隶属度矩阵
    U_new = zeros(size(U));
    for j = 1:c
        U_new(:,j) = 1./sum((dist(:,j)./dist).^2,2);
    end
    
    % 更新隶属度矩阵
    U = U_new.^m;
    
    % 判断是否收敛
    if norm(U-U_new) < 1e-6
        break;
    end
end

% 绘制聚类结果
[~,labels] = max(U,[],2);
scatter(data(:,1),data(:,2),[],labels,'filled');

其中，`data`是要聚类的数据，`c`是聚类数目，`m`是模糊因子，`max_iter`是迭代次数。在每次迭代中，先计算聚类中心，再计算距离矩阵和隶属度矩阵，并根据隶属度矩阵更新隶属度矩阵。最后根据隶属度矩阵确定每个数据点的聚类标签，并绘制聚类结果。