属性约简(MATLAB算法)有实例有讲解资料讲解.docxVIP

属性约简（MATLAB算 法）有实例有讲解 粗糙集属性约简 题目：SARS诊断 注：这是我学粗糙集时老师给的作业 U 干 咳 呼 吸 困 难 血 液 检 测 高 烧 38 度 X 射 线 浓 痰 血 细 胞 多 寒 战 肌 肉 酸 痛 乏 力 胸 膜 痛 头 痛 非 典 1 1 1 1 1 0 0 0 0 1 1 0 1 1 2 0 0 0 0 0 0 0 0 0 0 0 0 0 3 1 0 1〔 0 ] 0 : 0 ] 0 0 0 1 0 0 0 4 0 0 0 1 1 : 1 1 0 1 0 1 1 0 5 1 0 0 1 1 1 1 1 0 1 1 0 0 6 ] 0 1 0 ] 1」 1 J 1 1 1 1 0 0 1 0 7 1 0 0 0 1 1 1 0 0 1 1 1 0 8 1 1 1 11 0 0 0 0 ] 1 1 0 1 1 9 1 0 1 ] 1」 1 : 0 : 0 0 1 1 0 1 1 10 1 1 1 1 0 0 0 0 1 1 0 1 1 11 1 0 1 1 1 : 0 0 0 1 1 0 1 1 12 1 0 1 1 1 0 0 0 1 1 0 1 1」 属性约简的顺序如下：求正域、生成未经处理的区分矩阵、对区分矩阵进行化 简、求核、对已经处理过的区分矩阵进行属性约简。约简后的决策表有 26行, 所有12个属性都是正域中的属性，核为空 程序： % mai n.m tic; a=[ 1,1,1,1,0,0,0,0,1,1,0,1; 0,0,0,0,0,0,0,0,0,0,0,0; 1,0,1,0,0,0,0,0,0,1,0,0; 0,0,0,1,1,1,1,0,1,0,1,1; 1,0,0,1,1,1,1,1,0,1,1,0; 0,1,0,1,1,1,1,1,1,0,0,1; 1,0,0,0,1,1,1,0,0,1,1,1; 1,1,1,1,0,0,0,0,1,1,0,1; 1,0,1,1,1,0,0,0,1,1,0,1; 1,1,1,1,0,0,0,0,1,1,0,1; 1,0,1,1,1,0,0,0,1,1,0,1; 1,0,1,1,1,0,0,0,1,1,0,1 ]； d=[1;0;0;0;0;0;0;1;1;1;1;1]; pos=posCD(a,d); dismat=dismatrix(a,d,pos); dism=disbe(dismat); core=cor(dism); [red,row]=redu(dism); time=toc % dismatrix.m %生成未经处理的区分矩阵dismat fun cti on dismat=dismatrix(a,d,pos) [m, n]=size(a); p=1; in dex1=0;i ndex2=0;i ndex=0; dis=-1*o nes(m*(m-1)/2, n); for i=1:m for j=i+1:m if (isxbel on gtopos(i,pos )~ isxbel on gtopos(j,pos))... |(~isxbel on gtopos(i,pos )isxbel on gtopos(j,pos))... |(isxbel on gtopos(i,pos )isxbel on gtopos(j,pos )~ isxybel on gto in dD(i,j,d)) in dex2=1; end if in dex2==0 con ti nue; end for k=1: n if a(i,k)~=a(j,k) dis(p,k)=1; in dex1=1; else dis(p,k)=0; end end if in dex1==1 p=p+1;i ndex=1; end in dex 1=0; in dex2=0; end end if p=m*(m-1)/2 if in dex==0 dismat=[]; return; end if dis(p,1)==-1 P=P-1； end else p=m*(m-1)/2; end dismat=dis(1:p,:); % redu.m %对已经处理过的区分矩阵进行知识约简 fun cti on [red,row]=redu(dism) [m, n]=size(dism); red=[]; row=0; if m=0 return; end for i=1: n if dism(1,i)~=0 row=row+1; end end red(1:row,:)=zeros(row, n); j=1; for i=1:row while dism(1,j)==0 j=j+1; end red(i,j)=1; j=j+1; end temp=[];temp