《模糊聚类分析及matlab程序实.pdf》由会员分享,可在线阅读,更多相关《模糊聚类分析及matlab程序实.pdf(3页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、模糊聚类分析及 matlab 程序实现采用模糊数学语言对按一定的要求进行描述和分类的数学方法称为模糊聚类分析。聚类分析主要经过标定和聚类两步骤。【1】1 标定(建立模糊相似矩阵)城市居民食品零售价格,第t 时刻第 i 种食品的零售价记为),(tix。相似矩阵 R的构建方法:NTV法设时间序列),(jiA表示食品 i 在时间 t 的价格,其中 i=1,242;t=1,239。mkjkikmkjkikxxxxjiR11),max(1),((其中 i,j,k=1,242,m=39)42*42),(jiRR2 聚类2.1 计算 R的传递闭包:对模糊相似矩阵R,依次用平方法计算,2R,4R,t2R,当第
2、一次出现kkkRRR*时,则称kR为传递闭包。【1】2.2 开始聚类:【2】(1)令 T=1,2,3 42,取)1(xiT,令 X、Q为空集;(2)令0j;(3)若),(jxiR且Xxj,则令 jXX,jQQ;(4)1jj;(5)若nj,返回(1);(6)若 Q为空集,怎输出聚类x,X-TT;(7))1(xiQ,xiQQ,返回(2)。设置不同的置信水平值,就可以得到不同的分类。Matlab 程序实现:A=data;N M=size(A);for i=1:N for j=1:N R(i,j)=abs(1-sum(abs(A(i,:)-A(j,:)/sum(max(A(i,:);A(j,:);en
3、d endfor j=1:42 for i=1:42 y(i,j)=0;for k=1:42 mn(k)=min(R(i,k),R(k,j);end y(i,j)=max(mn);end end numda=1 0.9 0.95 0.85 0.8 0.75 0.55 0.7 0.655 0.65 0.6 0.55 0.5 0.45 0.454 0.4 0.45 0.3 0.35 0.255 0.25 0.2 0.15 0.1;for i=1:42 TT(i)=i;end for i=1:length(numda)disp(当分类系数是);disp(numda(i);a=numda(i);T=TT;disp(分类为);while 1 if isempty(T)xi=T(1);end X=;Q=;while 1 for j=1:42 if(y(xi,j)=a)&isempty(intersect(X,j)X=union(X,j);Q(length(Q)+1)=j;end end if isempty(Q)disp(X);break else xi=Q(1);Q(1)=;end end T=setdiff(T,X);if isempty(T)break end end end