insysind2.m

用matlab语言编的基于粗糙集理论的属性约简程序

function c=insysind2(a,x)%求不完备系统的等价类

[arow,acol]=size(a);
for k=1:acol
    z(:,k)=x(:,a(1,k));
end
x=z;
[p,q]=size(x);%获得矩阵的行数和列数
for m=1:p
    c{m}=[];
end

for i=1:p%对每行
   %y=1;
    for j=i:p%对每行
       %col=length(setdiff(x(i,:),x(j,:)));
        if(i~=j)
        col = setdiff1(x(i,:),x(j,:));
%         colzhong1 = setdiff1(x(i,:),x(j,:));
%         colzhong2 = setdiff1(x(j,:),x(i,:));
%        
%         col = colzhong2 + colzhong1 - length(setxor(zhong1,zhong2));
        if  (q-col)/q>=0.9%0.99%0.75
       % a(i,y)=j;
      %  y=y+1;
      c{i}=cat(2,[c{i}],[j]);
        c{j}=cat(2,[c{j}],[i]);
        end
        end
    end

     %c{i}=cat(2,[c{i}],a(i,:));
 %    clear a;    
 
end 
for i = 1:p
    c{i}=cat(2,c{i},[i]);
end