首页 > 编程语言 >模糊聚类FCM算法和基于GA遗传优化的FCM聚类算法matlab仿真

模糊聚类FCM算法和基于GA遗传优化的FCM聚类算法matlab仿真

时间:2023-01-03 19:47:17浏览次数:43  
标签:pop 算法 GA 聚类 FCM size

1.算法概述

      在各种糊聚类算法中 ,模糊C -均值聚类算法FCM (FuzzyC MeanClusteringAlgorithm )的应用最为广泛.但在实际的应用中 ,FCM算法却容易陷入局部最优解.因此 ,提出了一种基于遗传算法GA(GeneticAlgorithm )的模糊聚类分析方法 ,它利用了遗传算法随机搜索的特点 ,可以避免陷入局部最优解.实验表明 ,将该遗传算法与FCM算法结合起来进行聚类分析 ,比单一使用遗传算法或单一使用FCM算法进行聚类分析的效果都要好.

 

       针对FCM聚类,GA-FCM聚类采用了GA与FCM(GA-FCM)聚类进行全局搜索,克服了FCM聚类容易陷入局部极小值的问题.GA-FCM聚类是以全局最优个体将GA聚类用一个最优个体,迭代过程中既包括了GA运算.它依据GA的随机性扩大了搜索范围,之后在所找到的个体.

 

       模糊C均值算法是普通C均值算法的改进,普通C均值算法对于数据的划分是硬性的,而FCM则是一种柔性的模糊划分。

       首先,需要一种方法度量两个样本的相似性,这个就是距离。FCM是一种重要的聚类算法,其目标是将n维空间中的数据X = {x_1, ..., x_N}分配到C个聚类中心v_1, ..., v_C。在欧氏距离意义下,数据靠近哪个聚类中心就属于哪个类。

 

      把这种一个样本属于结果的这种相似的程度称为样本的隶属度,一般用表示,表示一个样本相似于不同结果的一个程度指标。

 

        基于此,假定数据集为X,如果把这些数据划分成c类的话,那么对应的就有c个类中心为C,每个样本j属于某一类i的隶属度为,则定义一个FCM目标函数及其约束条件如下所示:

 

 

 

      分别对其中的变量、求导数,求该目标函数的极值。

 

      在最终的结果中,与是相互包含的,在fcm算法的初始阶段既没有,也没有,所以在程序开始时可以任意赋值给和其中一个,迭代计算。

 

FCM类中心的更新法则:

 

 

 

        看式子可以发现也是一个加权平均,类i确定后,首先将所有点到该类的隶属度求和,然后对每个点,隶属度除以这个和就是所占的比重,乘以​就是这个点对于这个类i的贡献值了。   

 

2.仿真效果预览

matlab2022a仿真结果如下:

 

 

 

 

3.MATLAB部分代码预览

 

data=rand(100,2).*10;
max=10;%最大值
min=0;%最小值
cluster_n=5;%所要聚类的个数
pop_size=150;%种群大小
Pc=0.25;%交叉概率
Pm=0.05;%变异概率
max_run=300;%最大运行次数
N=size(data,2);%特征维数
Num_boolean=8;%表示“7.5“所要得二进制个数
G_F=zeros(pop_size,1);%染色体适应度
U=zeros(cluster_n, size(data,1));%初始化隶度矩U
expo=2;%权值m大小
SL=2;
min_impro=0.0001;
 
 
G=round(rand(pop_size,cluster_n*N*Num_boolean));
for h=1:max_run
for i=1:pop_size %遍历整个种群
     center=zeros(cluster_n,N);%初始化单个染色体中心
     center=code(G(i,:),cluster_n,N,max,min,Num_boolean);%(1)对染色体解码,得到实数表示的聚类中心
     [U,obj_fitness(i),obj(i)]=Uconvert(data, cluster_n, expo,center);
 end
 	if i > 1,
		if abs(obj(i) - obj(i-10)) < min_impro, break; end,
	end
 T(h)=h;    
 Fit=obj;
 [Order,Index]=sort(Fit);
 BF=Order(1);
 BFI=BF;
 BG=G(Index(1),:);
 ln=1;
 for i=1:1:SL
     BGG(i,:)=G(Index(ln),:);
     ln=ln+1;
 end 
 Best_value(h)=Order(1);
 Best_chrom=G(Index(1),:);
 
 
 NG=selection(G,obj_fitness,pop_size);
 NG=crossover(NG,Pc,pop_size);
 NG=mutation(NG,Pm,pop_size);
 Rs=4;
 for i=1:1:SL
     NG(pop_size-Rs,:)=BGG(i,:);
     Rs=Rs-1;
 end
 G=NG;
end
center
b=sort(Best_value);
Pbest=b(1)
 
plot(T,Best_value);
plotflag=2;
M=2;
[U2,P,Dist,Cluster_Res,Obj_Fcn,iter]=fuzzycm(data,cluster_n,plotflag,M)
A_019

 

  

 

标签:pop,算法,GA,聚类,FCM,size
From: https://www.cnblogs.com/51matlab/p/17023187.html

相关文章

  • jsjiami.v6算法分析
    这是一个例子:constjsjiami={v6:(str)=>{letresult='';for(leti=0;i<str.length;i++){result+=String.fromCharCode(str.charCodeAt(i)......
  • 大数的高精度算法
    高精度加法:#include<bits/stdc++.h>usingnamespacestd;stringadd(string&a,string&b){stringc;intt=0;if(a.length()<b.length())returnadd......
  • KMP算法
    KMP算法参考:如何更好地理解和掌握KMP算法?模式串匹配算法,在一个主串(文本串s)中查找子串(模式串p)第一次匹配的位置算法两个关键操作根据模式串建立next数组根据next......
  • 算法竞赛进阶指南 0x43 线段树
    文章目录​​线段树简介​​​​线段树的简单代码实现​​​​建树代码​​​​修改操作​​​​查询操作​​​​线段树的查询操作的时间复杂度分析:​​​​[AcWing245.你......
  • 美颜sdk磨皮算法与人脸皮肤识别技术
    在之前的文章中小编曾提起过,美颜sdk以及其它美颜工具的核心技术都是人脸关键点识别,只有先识别人脸关键点才能进行后续的美颜操作。今天小编要讲的美颜sdk磨皮算法同样不例外......
  • 小程序 SHA1加密算法使用
    创建一个js文件,或写入util.js中//SHA1加密functionencodeUTF8(s){vari,r=[],c,x;for(i=0;i<s.length;i++)if((c=s.charCodeAt(i))<0x80)r.pu......
  • 排序算法
    选择排序选择排序(Selectionsort)是一种简单直观的排序算法。它的工作原理如下。首先在未排序序列中找到最小(大)元素,存放到排序序列的起始位置,然后,再从剩余未排序元素中继续......
  • 每日算法之把二叉树打印成多行
    JZ78把二叉树打印成多行题目给定一个节点数为n二叉树,要求从上到下按层打印二叉树的val值,同一层结点从左至右输出,每一层输出一行,将输出的结果存放到一个二维数组中......
  • 工程师应该学点算法——图论2
    为什么QQ要给女朋友推送前女友?这还是从图的算法说起。前篇->​​图论1​​图的遍历在图的遍历中我们一定要掌握两种最基础的算法:深度优先和广度优先。深度优先遍历(DFS)这......
  • 算法刷题 Day 6 | 242.有效的字母异位词 349. 两个数组的交集 202. 快乐数 1. 两数之
    哈希表理论基础建议:大家要了解哈希表的内部实现原理,哈希函数,哈希碰撞,以及常见哈希表的区别,数组,set和map。什么时候想到用哈希法,当我们遇到了要快速判断一个元素是......