matlab 实现模糊C均值聚类

1. 原理

模糊c均值算法步骤:

  1. 设定聚类数目c和加权指数b:

  2. 初始化各个聚类中心m

  1. 重复下面的运算,直到各个样本的隶属度值稳定:用当前的聚类中心根据下式计算隶属度函数:

用当前的隶属度函数按下式更新计算各类聚类中心:

当模糊C均值算法收敛时,就得到了各类的聚类中心和各个样本对于各类的隶属度值,从而完成了模糊聚类划分

Matlab中模糊c均值聚类函数(b=2):

Center,U,obj_fcn\]=fcm(data,cluster_n) data:要聚类的数据集合,每一行为一个样本; cluster_n:聚类数; Center:最终的聚类中心矩阵,每一行为聚类中心的坐标值; U:最终的模糊分区矩阵; obj_fcn:在迭代过程中的目标函数值。 ## 2. 过程 ### 2.1 数据集合并及可视化 ![](https://file.jishuzhan.net/article/1766375306697904130/b5b553864d385a0234325407c0343002.webp) ![](https://file.jishuzhan.net/article/1766375306697904130/0b0afe0fe36e9dcdcdca5dca27a75f79.webp) 将MALE.txt和FEMALE.txt文件合并成一个data.txt文件,在matlab中导入data文件并绘制样本点分布: ![](https://file.jishuzhan.net/article/1766375306697904130/e954d1ddadde77958077447c30247fa1.webp) 图1 初始样本点分布图 ### 2.2 模糊c均值聚类 调用MATLAB中模糊c均值聚类函数fcm实现聚类,初始设置聚类簇的个数为2,求出每个样本的最大隶属度,然后找出分别属于两类样本的点并绘制图像: ![](https://file.jishuzhan.net/article/1766375306697904130/25c3a7b11b09e597a8f2edda589c1e79.webp) 程序如上图所示,运行后聚类结果如下: ![](https://file.jishuzhan.net/article/1766375306697904130/6066cc914845adb81b39b23bc375ae75.webp) 图2 迭代过程中目标函数值 可以看出目标函数值在迭代中是不断减小的 ![](https://file.jishuzhan.net/article/1766375306697904130/9581ba2766a8571dd7348f99f61f8f57.webp) 图3 聚类结果 ![](https://file.jishuzhan.net/article/1766375306697904130/75774f01d3dfb1ea6fcdab3f3d949f31.webp) 图4 簇为2聚类结果样本点分布 为了能够寻找到更合理的聚类数目,接下来分别进行3类、4类、5类聚类 ![](https://file.jishuzhan.net/article/1766375306697904130/ceb962cdf4eef8538fd759ef6b1bff52.webp) 图5 簇个数为3 ![](https://file.jishuzhan.net/article/1766375306697904130/6a3d1f884da35e663ce6190126f9a2dd.webp) 图6 簇个数为4 ![](https://file.jishuzhan.net/article/1766375306697904130/260175bef6b24330dfe8cbe026c89abf.webp) 图7 簇个数为5 ## 3. 结果与分析 在上述实验过程中,我们成功实现了模糊c均值聚类,因为数据集是由两类样本合成的所以初始聚类簇的个数设置为2,结果如图4所示,然后从图中可以发现聚类效果似乎并非最佳,于是后续又增大了簇的个数进行多次实验,如图5到图7所示,通过观察可知聚类个数为4时或许会更加合理。在实验迭代的过程中可以了解到这种聚类的方式主要是通过优化目标函数然后得到每个样本点对所有类中心的隶属度,从而决定样本点的类属以达到自动对样本数据进行聚类的目的。同时我也通过改变聚类个数的方式去寻找更加合理的类别数。但是实验中还缺少对于聚类优劣的描述,可以通过一些聚类指标来对实验进行进一步优化。

相关推荐
程思扬4 天前
利用JSONCrack与cpolar提升数据可视化及跨团队协作效率
网络·人工智能·经验分享·docker·信息可视化·容器·架构
月阳羊4 天前
【硬件-笔试面试题-95】硬件/电子工程师,笔试面试题(知识点:RC电路中的时间常数)
java·经验分享·单片机·嵌入式硬件·面试
tyatyatya4 天前
MATLAB中进行视觉检测入门教程
开发语言·matlab·视觉检测
智者知已应修善业4 天前
【矩阵找最大小所在位置】2022-11-13
c语言·c++·经验分享·笔记·算法·矩阵
SccTsAxR4 天前
[C语言]常见排序算法①
c语言·开发语言·经验分享·笔记·其他·排序算法
智者知已应修善业4 天前
【51单片机单按键控制2个LED循环闪烁】2022-12-7
c语言·经验分享·笔记·嵌入式硬件·51单片机
艾莉丝努力练剑4 天前
【C++】类和对象(下):初始化列表、类型转换、Static、友元、内部类、匿名对象/有名对象、优化
linux·运维·c++·经验分享
2zcode5 天前
基于Matlab可见光通信系统中OOK调制的误码率性能建模与分析
算法·matlab·php
听情歌落俗5 天前
MATLAB3-1变量-台大郭彦甫
开发语言·笔记·算法·matlab·矩阵