matlab模糊聚类程序

1、实用标准文档3.数据标准化1) 数据矩阵设论域 U=x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12 为被分类的对象，每个 对 象 又 由 指 标 Y=y1,y2,y3,y4,y5,y6,y7,y8,y9 表 示 其 性 状 即xi=xi1,xi2,xi3,xi4,xi5,xi6,xi7,xi8,xi9,xi10,xi11,xi12 (i=1,2, ,12) 于是得到原是数据矩阵75250134212178219238437832959653720541326531331362123121814178691127810436943147232536111211246

2、16A=1013253528652413894286788203291695158724930483714632791126928969792949932754642417231149187951218382) 数据标准化将模糊矩阵的每一个数据压缩到0,1 上，采用平移 .极差变换进行数据标准B(i,k)=A(i,k)- min A(i,k)1 i nm1ai nx A(i,k)- m1 iinn A(i,k)1 i n 1 i n(k=1,2, ,m)运用 matlab 编程由函数 F_jisjbzh.m 【见附录 3.4 】的标准化矩阵是附录 3.4function X=F_JISjBzh

3、(cs,X)%模糊聚类分析数据标准化变换文案大全实用标准文档%X 原始数据矩阵； cs=0 ，不变换； cs=1 ，标准差变换 %cs=2, 极差变换 if(cs=0) return ;endn,m=size(X);% 获得矩阵的行列数 if(cs=1) % 平移极差变换for(k=1:m) xk=0;for(i=1:n) xk=xk+X(i,k);end xk=xk/n;sk=0;for(i=1:n) sk=sk+(X(i,k)-xk)2;end sk=sqrt(sk/n);for(i=1:n) X(i,k)=(X(i,k)-xk)/sk;end endelse % 平移 * 极差变换for

4、(k=1:m) xmin=X(1,k);xmax=X(1,k);for(i=1:n)if(xminX(i,k) xmin=X(i,k);end if(xmaxX(i,k) xmax=X(i,k);end endfor(i=1:n) X(i,k)=(X(i,k)-xmin)/(xmax-xmin);end endend文案大全实用标准文档0000000.03190.028600.01560.13950.04840.18390.08650.01470.4043 0 2 08.6 0.24310.23750.27910.45970.68970.35580.27940.57450.28570.1667

5、0.14370.09300.23390.35630.20190.32350.12770.42860.08330.53440.74420.88710.83911.00000.51471.00000.80000.31250.05000.23260.27420.06900.11540.14710.25530.08570.09720.29380.31400.41130.54020.82690.75000.43621.00000.63890.06560.01160.04030.03450.07690.0147000.18750.50620.53490.45160.77010.47120.42650.51

6、060.97141.00001.00001.00001.00001.00000.85581.00000.84040.74290.32640.26870.25580.41940.67820.23080.23530.24470.22860.32640.03440.02330.056500.00960.01470.191500.0417B=第二步：标定(建立模糊相似矩阵)对标定我们运用了直接欧几里得距离法：rij 1 c d(x i ,x j )其中 c为任意选区的参数，它使得 0rij 1, d(x i ,x j )表示 xi与xj的距离，d(xi,xj )xjk)2运用 matlab 软件编写

7、F_jir.m 函数【见附录 3.5 】，取 cs=8 ，的模糊相似矩阵附录 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)-X(j,k)2;endd=sqrt(d);%直接海明距离法elseif(cs=9)for(k=1:m)d=d+abs(X(i,k)-X(j,k);

8、end%直接切比雪夫距离法elsefor(k=1:m)if(dabs(X(i,k)-X(j,k)d=abs(X(i,k)-X(j,k); endend文案大全实用标准文档endif(Cd)C=d;endendendC=1/(1+C);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)2; end d=sqrt(d);%直接海明距离法elseif(cs=9)for(k=1:m) d=d+abs(X(i,k)-X(j,k); end%直接切比雪夫距离法文案大全实用标准文档elsefor(k=1:m)if(da

9、bs(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.46340.67310.76381.00000.81400.66940.79950.81500.81401.00000.53490.37150.46340.66940.53491.0000R0.86680.89730.77360.85340.48630.49300.56080.69610.67050.71040.93830.87460.69070.81050

10、.39280.46020.54900.68120.62040.69050.27450.35410.56180.44490.78630.71510.78660.89070.84910.59980.94990.89720.70160.80630.40013）聚类（求动态聚类图）0.8668 0.4930 0.9383 0.4602 0.2745 00.8973 04.65 6 0 08. 5 4 09.08 7 0.3541 0.7866 0.8970.77360.69610.69070.68120.561800.85340.67050.81050.62040.444900.48630.71040

11、.39280.69050.786301.00000.58010.87550.54940.388100.58011.00000.52160.80260.619900.87550.52161.00000.49590.297900.54940.80260.49591.00000.621400.38810.61990.29790.62141.000000.79910.67830.74460.68520.51610.89720.50910.93000.48020.30020 传递闭包法根据标定所得的模糊矩阵，只是一个模糊相似矩阵 R，不一定具有传递性，文案大全实用标准文档即 R 不一定是模糊等价矩阵，还

12、需要对其改造成模糊等价矩阵R ,根据定理，用二次方法求传递闭包 t(R)，t(R)就是所求模糊等价矩阵 R,即： t(R) =R ,再让 由大变到小，就可形成动态聚类图。通过 matlab 软件编的函数 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 ;endfor

13、(i=1:n) R(i,i)=1;% 修正错误for(j=i+1:n)if(R(i,j)1) R(i,j)=1;endR(i,j)=round(10000*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;endendImd(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

14、) 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(Imd(i)=Imd(x) js=js+1;Sz(js)=j;end;end flsz(x)=flsz(x)+1;文案大全实用标准文档endendendfor(i=1:k-1)for(j=i+1:k)if(flsz(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:

15、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;endif(pd) if(js=0) y=0;endfor(j=1:n) if(R(i,j)=Imd(x) js=js+1;Sz(js)=j;end;endflsz(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);end文案大全实用标准文档

16、for(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;end;endif(pd) if(js=0) y=0;endfor(j=1:n) if(R(i,j)=Imd(x) js=js+1;Sz(js)=j;end;endflsz(x)=fl

17、sz(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);end文案大全实用标准文档endendF_dtjltx=figure(name, 动态聚类图 ,color,w);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)

18、,y+Gd/2,int2str(xhsz(i);line(lx+i*Kd,lx+i*Kd,y,y-Gd);linesz(i)=lx+i*Kd;endtext(lx*1.5+i*Kd,y+Gd/2, 分类数 );y=y-Gd;for(x=1:fl)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+

19、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,i);js2=0;pd=0;for(j=1:flsz(x-1)js2=js2+flqksz(x-1,j);if(js2=js1) pd=1;break;endend if(j=js0)line(linesz(js0),linesz(j),y,y);end line(linesz(js0)+linesz(j)/2,(linesz(js0

20、)+linesz(j)/2,y,y-Gd);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图六：动态聚类图文案大全实用标准文档根据动态聚类图，选定不同的 的值，将就可以得到不同的分类附录 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)-X(j,k)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)if(dabs(X(i,k)-X(j,k)d=abs(X(i,k)-