模糊C均值聚类程序文件

./实验二模糊C均值聚类实验目的：学会使用MATLAB软件进行模糊C均值聚类,学会如何进行迭代并观察迭代过程。实验学时：4学时实验容：1、认真阅读guide.doc文件,通过给出的英文的例子学习进行C均值聚类的具体步骤。2、在学习完所给的例子后进行实际操作。根据表格提供的固定资本和人力资本等进行聚类分布。进一步熟悉和掌握熟悉FUZZYCLUSTERING.实验日期：20XX4月24日实验过程：查看所给数据表格〔如下,由经济学理论知,GDP的产出状况是由固定资本的投入和人力资源的投入决定的。因此,实际上我们只需要选取固定资本和人力资源这两组数据进行处理就行了。地区〔DistrictGDP固定资本<FixedCapital>人力资本<manpowerCapital>地区〔DistrictGDP固定资本<FixedCapital>人力资本<manpowerCapital>2.05765.94891.36000.78031.38201.00001.83284.03081.39900.60390.91710.66600.82361.62230.86901.36812.72981.02600.54401.03370.7960广西0.46600.83420.54600.65031.10990.93100.68590.81270.62001.20012.03140.38500.56500.81270.62000.75530.31141.02200.51181.04100.56300.93441.42081.28100.28560.57560.29903.06748.03961.74800.48401.01660.46601.29332.24501.08800.52751.35880.52401.46293.10380.89400.50401.03070.57400.51990.81120.61400.41650.85440.45901.23652.05231.15500.57541.5080.55000.51980.74940.78500.53001.50360.51801.04391.65340.94700.79812.02260.91100.59031.92100.6530平均0.89141.89360.85732、通过学习guide中的例,将所给的defcm.m程序进行重新编辑。其具体程序如下：function[NCentres,M]=defcm<Centres,q>Tiles=[5.94891.360014.03081.399012.03140.385011.42081.281018.03961.748012.24501.088013.10380.894012.05231.155011.65340.947012.72981.026011.62230.8690-11.03370.7960-11.10990.9310-10.31141.0220-10.81120.6140-10.74940.7850-11.92100.6530-11.38201.0000-10.91710.6660-10.83420.5460-10.81270.6200-10.81270.6200-11.04100.5630-10.57560.2990-11.01660.4660-11.35880.5240-11.03070.5740-10.85440.4590-11.5080.5500-11.50360.5180-12.02260.9110-1];%将固定资本和人力资本的数据按GDP的平均值进行分类,大于平均值的分为一类,记为1,小于平均值的分为一类,记为-1Tiles<:,1>=log<Tiles<:,1>>;Tiles<:,2>=log<Tiles<:,2>>;clf;holdoff;plot<Tiles<1:16,1>,Tiles<1:16,2>,'ob'>;axis<[-1.52.5-1.52.5]>;xlabel<'固定资本'>;ylabel<'人力资本'>;title<'Tilesdata:o=wholetiles,*=crackedtiles,x=centres'>;holdon;plot<Tiles<17:31,1>,Tiles<17:31,2>,'*b'>;plot<Centres<:,1>,Centres<:,2>,'xr'>;Object=Tiles<:,[12]>';Centres=Centres';k=length<Centres<:,1>>;D=dist<Object,Centres>;M=member<D,q>;ClustSums=sum<M>;[NoOfObs,NoOfCtrs]=size<M>;z=M./<ones<NoOfObs,1>*ClustSums>;NCentres=Centres';fori=1:NoOfCtrs,w=z<:,i>*ones<1,2>;NCentres<i,:>=sum<Object'.*w>;end;3、将上述程序保存为defcm.m文件,然后在命令窗口中输入初始聚类中心Cts=[0.50.5;0.80.8],即在这里面我们取了两个聚类中心点,c1=<0.5,0.5>、c2=<0.8,0.8>,并设q=2。需要注意的是,q越接近1收敛速度越快,但若q=1则为硬C均值聚〔HCM。然后不断调用函数defcm,[Cts,M]=defcm<Cts,q>。即在Command窗口中输入：Cts=[0.50.5;0.80.8]；q=2；[Cts,M]=defcm<Cts,q>。本次实验一共迭代了13次,M矩阵才不再变化,达到了最佳聚类中心。迭代的仿真图形如下所示：〔迭代时只要将Command窗口中Cts中的值换成前一次迭代输出的值即可第1次迭代仿真的图形第6次迭代仿真的图像第13次的迭代图形：下面为最终的聚类中心和M矩阵。Cts=0.0733-0.44950.7166-0.0713M=0.26830.73170.20940.79060.54220.45780.29450.70550.30960.69040.03820.96180.12410.87590.05700.94300.11870.88130.07800.92200.18220.81780.90660.09340.72260.27740.66880.33120.92680.07320.85650.14350.27870.72130.37590.62410.96500.03500.92440.07560.92740.07260.92740.07260.97700.02300.75010.24990.90480.09520.84260.15740.98180.01820.88650.11350.73170.26830.73870.26130.00120.99884、由于不能确保FCM收敛于一个最优解,所以我们可以采用不同的初始聚类中心启动该算法,以便用来比较结果。在Command窗口中输入：Cts=[0.51.0;1.52.0]；q=2；[Cts,M]=defcm<Cts,q>在这次实验一共迭代了15次,M矩阵才不再变化,达到了最佳聚类中心。在这里只给出第15次,Command窗口中显示的聚类中心值Cts和M矩阵：通过最终的迭代图可以看出：当初始聚类中心取值不同时,最后的聚类中心点的值也有细微的差别。FCM是否收敛于一个最优解,不仅与初始聚类中心的选择有关,还与加权指数q有关。由模糊集的知识我们知道,加权指数q取得过大,聚类效果会很差；q取得过小则算法会接近与HCM聚类算法。所以下面的仿真中,其他数据不变,取：q=1.1,聚类中心为Cts=[0.50.5;0.80.8]。经过MATLAB仿真后知道,此时只需迭代8次就可达到了最佳聚类中心。下面只给出第八次迭代的相关数据：实验小结：通过这次实验对模糊C均值聚类〔FCM有了近一步的了解,并且结合例子能够进行简单的聚类中心迭代,寻找最佳聚类中心。另外,通过q