matlab模糊聚类程序

1、实用标准文档3.数据标准化(1)数据矩阵设论域URx!,”,。,%,外,%,%,%;,”?为被分类的对象，每个Xi标 Y =yY2y,y§y6生表，刀x 4ix6 x6i x,nx8ix9ix(i=M2,x i752501342、121782192384378329596537205413265313313621231218141786911278104 36943147232536111211246161013253528652413894286788203291695158724930483714632791126928969792949932754642417231149187

2、95121838 >象又由指=Xi X2i X原是数据矩阵A=，其州状即1)2)于是得到(2)数据标准化将模糊矩阵的每一个数据压缩到0,1上，采用平移.极差变换进行数据标准A(i,k)-B(i,k尸rpjn A(i,k)m1ax A(i，k)-rpiinA(i,k)(k=1,2,m)运用matlab编程由函数F_jisjbzh.m【见附录3.4 的标准化矩阵是附录3.4function X=F_JISjBzh(cs,X)峨糊聚类分最据标准化变换%X®始数据矩阵；cs=0,不变换；cs=1,标准差变换%cs=2极差变换if(cs=0) return ;endn,m=size(X)

3、;%获得矩阵的行列数if(cs=1) % 平移极差变换for(k=1:m) xk=0;for(i=1:n) xk=xk+X(i,k);endxk=xk/n;sk=0;for(i=1:n) sk=sk+(X(i,k)-xk)A2;endsk=sqrt(sk/n);for(i=1:n) X(i,k)=(X(i,k)-xk)/sk;end endelse %平移*极差变换for(k=1:m) xmin=X(1,k);xmax=X(1,k);for(i=1:n)if(xmin>X(i,k) xmin=X(i,k);endif(xmax<X(i,k) xmax=X(i,k);endendfo

4、r(i=1:n) X(i,k)=(X(i,k)-xmin)/(xmax-xmin);endendend0000000.03190.028600.01560.13950.04840.18390.08650.01470.4043 02a6 0.24310.23750.27910.45970.68970.35580.27940.57450.28570.16>670.14370.09300.23390.35630.20190.32350.12770.42860.081330.53440.74420.88710.83911.00000.51471.00000.80000.3125B=0.05000

5、.23260.27420.06900.11540.14710.25530.08570.09)720.29380.31400.41130.54020.82690.75000.43621.00000.63；890.06560.01160.04030.03450.07690.0147000.18750.50620.53490.45160.77010.47120.42650.51060.97141.00)001.00001.00001.00001.00000.85581.00000.84040.74290.32640.26870.25580.41940.67820.23080.23530.24470.

6、22860.3264、0.03440.02330.056500.00960.01470.191500.0417第二步：标定（建立模糊相似矩阵） 对标定我们运用了直接欧几里得距离法:rj =1 -c d(xi,xj)其中c为任意选区的参数，它使得0& rj W1, d（xi,xj）表示x与xj的距离,d（xi,xj） =（。- xjk）2 k 1运用matlab软件编写F_jir.m 函数【见附录3.5，取cs=8,的模糊相 似矩阵附录3.5 :（仅附录了一段用到的程序）function R=F_jir（cs,X）%cs=8直接欧几里得距离法%cs=9直接海明距离法（绝对值减数法）%cs

7、=10直接切比雪夫距离法文案大全实用标准文档elseif(cs<=10)C=0;for(i=1:n)for(j=i+1:n)d=0;% 直接欧几里得距离法if(cs=8)for(k=1:m)d=d+(X(i,k)-X(j,k)A2;endd=sqrt(d);% 直接海明距离法elseif(cs=9)for(k=1:m)d=d+abs(X(i,k)-X(j,k);end% 直接切比雪夫距离法elsefor(k=1:m)if(d<abs(X(i,k)-X(j,k) d=abs(X(i,k)-X(j,k);endendendif(C<d)C=d;endendendC=1/(1+C)

8、;for(i=1:n)for(j=1:n)d=0;% 直接欧几里得距离法if(cs=8)for(k=1:m)d=d+(X(i,k)-X(j,k)A2;endd=sqrt(d);% 直接海明距离法elseif(cs=9)for(k=1:m)d=d+abs(X(i,k)-X(j,k);end%直接切比雪夫距离法elsefor(k=1:m)if(d<abs(X(i,k)-X(j,k) d=abs(X(i,k)-X(j,k);endendendR(i,j)=1-C*d;endend1.00000.85960.67310.79950.37150.85961.00000.76380.81500.46

9、340.67310.76381.00000.81400.66940.79950.81500.81401.00000.53490.37150.46340.66940.53491.0000R =0.86680.89730.77360.85340.48630.49300.56080.69610.67050.71040.93830.87460.69070.81050.39280.46020.54900.68120.62040.69050.27450.35410.56180.44490.78630.71510.78660.89070.84910.5998、0.94990.89720.70160.8063

10、0.4001(3)聚类(求动态聚类图)<1>专递闭包法(3)聚类(求动态聚类图)<1>专递闭包法0.86680.49300.93830.46020.2745(0.8973465608549087 0.3541 0.7866 0.89,0.77360.69610.69070.68120.5618 (0.85340.67050.81050.62040.4449 (0.48630.71040.39280.69050.7863 (1.00000.58010.87550.54940.3881 (0.58011.00000.52160.80260.6199 (0.87550.521

11、61.00000.49590.2979 (0.54940.80260.49591.00000.6214 (0.38810.61990.29790.62141.0000 (0.79910.67830.74460.68520.51610.89720.50910.93000.48020.3002 (根据标定所得的模糊矩阵，只是一个模糊相似矩阵 R,不一定具有传递性,即R不一定是模糊等价矩阵，还需要对其改造成模糊等价矩阵R',根据定理,用二次方法求传递闭包t (R), t (R)就是所求模糊等价矩阵R',即：t (R) =R , 再让儿由大变到小，就可形成动态聚类图。通过matlab软

12、件编的函数F_JIDtjl.m 【见附录3.6，得到动态聚类图或 者直接运用matlab软件编的函数F_Jlfx.m【见附录3.7，运行F_Jlfx(2,8,A) 得动态聚类图是:附录3.6 :function F_JIDtjl(R)%定义函数%真糊聚类分析动态聚类%R 真糊相似矩阵m,n=size(R);% 获得矩阵的行列数 if(m=n|m=0) return ;end for(i=1:n) R(i,i)=1;%修正错误for(j=i+1:n)if(R(i,j)<0) R(i,j)=0;elseif(R(i,j)>1) R(i,j)=1;endR(i,j)=round(1000

13、0*R(i,j)/10000;%保留四位小数R(j,i)=R(i,j);endendjs0=0;while(1)% 求传递闭包R1=Max_Min(R,R);% 【见附录3.6.1 】js0=js0+1;if(R1=R) break;else R=R1;end endImd(1)=1;k=1;for(i=1:n) for(j=i+1:n) pd=1;% 找出所有不相同的元素for(x=1:k)if(R(i,j)=Imd(x) pd=0;break;end;endif(pd) k=k+1;Imd(k)=R(i,j);endend;endfor(i=1:k-1) for(j=i+1:k)if(Im

14、d(i)<Imd(j)% 从大到小排序x=Imd(j);Imd(j)=Imd(i);Imd(i)=x;end;end;endfor(x=1:k) 哪Imd (x)分类，分类数为flsz(x),临时用Sz记录元素序 号js=0;flsz(x)=0;for(i=1:n) pd=1;for(y=1:js) if(Sz(y)=i) pd=0;break;end;endif(pd)for(j=1:n)if(R(i,j)>=Imd(x) js=js+1;Sz(js)=j;end;endflsz(x)=flsz(x)+1;endendendfor(i=1:k-1)for(j=i+1:k)if(f

15、lsz(j)=flsz(i) flsz(j)=0;end;end;endfl=0;% 排除相同的分类for(i=1:k) if(flsz(i) fl=fl+1;Imd(fl)=Imd(i);end;endfor(i=1:n) xhsz(i)=i;endfor(x=1:fl)% 获得分类情况：对元素分类进行排序js=0;flsz(x)=0;for(i=1:n) pd=1;for(y=1:js) if(Sz(y)=i) pd=0;break;end;end文案大全实用标准文档if(pd) if(js=0) y=0;endfor(j=1:n)js=js+1;Sz(js)=j;end;endflsz(

16、x)=flsz(x)+1;Sz0(flsz(x)=js-y;endendjs0=0;for(i=1:flsz(x)for(j=1:Sz0(i) Sz1(j)=Sz(js0+j);endfor(j=1:n) for(y=1:Sz0(i)if(xhsz(j)=Sz1(y)js0=js0+1;Sz(js0)=xhsz(j);end;end;endendfor(i=1:n) xhsz(i)=Sz(i);endendfor(x=1:fl)% 获得分类情况：每一子类的元素个数js=0;flsz(x)=0;for(i=1:n) pd=1;for(y=1:js) if(Sz(y)=i) pd=0;break;

17、end;endif(pd) if(js=0) y=0;endfor(j=1:n)js=js+1;Sz(js)=j;end;endflsz(x)=flsz(x)+1;Sz0(flsz(x)=js-y;endendjs0=1;for(i=1:flsz(x) y=1;for(j=1:flsz(x)if(Sz(y)=xhsz(js0)flqksz(x,i)=Sz0(j);js0=js0+Sz0(j);break;endy=y+Sz0(j);endendendF_dtjltx=figure('name',' 动态聚类图','color','w

18、9;);axis('off');Kd=30;Gd=40;y=fl*Gd+Gd;lx=80;text(24,y+Gd/2,' 入')；for(i=1:n)text(lx-5+i*Kd-0.4*Kd*(xhsz(i)>9),y+Gd/2,int2str(xhsz(i);line(lx+i*Kd,lx+i*Kd,y,y-Gd);if(R(i,j)>=Imd(x)if(R(i,j)>=Imd(x)linesz(i)=lx+i*Kd;end text(lx*1.5+i*Kd,y+Gd/2,' 分类数 ');y=y-Gd;for(x=1:f

19、l)text(8,y-Gd/2,num2str(Imd(x);js0=1;js1=0;if(x=1)for(i=1:flsz(x)js1=flqksz(x,i)-1;if(js1)line(linesz(js0),linesz(js0+js1),y,y);endline(linesz(js0+js1)+linesz(js0)/2,(linesz(js0+js1)+linesz(js0 )/2,y,y-Gd);linesz(i)=(linesz(js0+js1)+linesz(js0)/2;js0=js0+js1+1;endelse for(i=1:flsz(x)js1=js1+flqksz(x

20、,i);js2=0;pd=0;for(j=1:flsz(x-1)js2=js2+flqksz(x-1,j);if(js2=js1) pd=1;break;end endif(j=js0)line(linesz(js0),linesz(j),y,y);endline(linesz(js0)+linesz(j)/2,(linesz(js0)+linesz(j)/2,y,y-G d);linesz(i)=(linesz(js0)+linesz(j)/2;js0=j+1;end;endtext(2*lx+n*Kd,y-Gd/3,int2str(flsz(x);y=y-Gd;end图六：动态聚类图文案大

21、全实用标准文档10.94990.93830阴73 0.3S72 0.8907 0.85340.84910.30260.78630.7104 口 6961根据动态聚类图，选定不同的 人的值，将就可以得到不同的分类附录3.5 :(仅附录了一段用到的程序)function R=F_jir(cs,X)%cs=8直接欧几里得距离法%cs=9直接海明距离法(绝对值减数法)%cs=10直接切比雪夫距离法elseif(cs<=10)C=0;for(i=1:n)for(j=i+1:n)d=0;%直接欧几里得距离法if(cs=8)for(k=1:m)d=d+(X(i,k)-XGk)，2;endd=sqrt(d);% 直接海明距离法elseif(cs=9)for(k=1:m)d=d+abs(X(i,k)-X(j,k); end% 直接切比雪夫距离法elsefor(k=1:m)i