1、粗糙集 -屬性約簡 -matlab 程序Data2為條件屬性，decision2為決策屬性%my_test 函數(shù)實現(xiàn)clc;% 讀取信息系統(tǒng)文件file = textread('data2.txt','%s','delimiter','n','whitespace',''); %讀取文件信息，每一行為一個胞元 m,n=size(file); % 胞元的大小for i=1:mwords=strread(filei,'%s','delimiter',' '

2、);% 讀取每個胞元中字符，即分解胞元為新的胞元 words=words'% 轉(zhuǎn)置Xi=words;endX=X'% B,num,AT=my_reduct(X);%信息系統(tǒng)的約簡ind_AT=ind(X); %信息系統(tǒng)的不可等價關(guān)系% 顯示約簡信息系統(tǒng)disp('約簡后的條件系統(tǒng)為：)；m,n=size(B);for i=1:mdisp(Bi)；end% 讀取決策系統(tǒng)文件file = textread('decision2.txt','%s','delimiter','n','whitespace&

3、#39;,'')；m,n=size(file)；for i=1:mwords=strread(filei,'%s','delimiter',' ')；words=words'；Di=words；endD=D'；% 決策系統(tǒng)的正 域約簡X_D=X；l,k=size(X_D1)；pos_d=pos(X_D,D)；% 正域fori=1:m%正域有問題% if(ismember(num(i),pos_d)Bi=''% 若約簡后的信息系統(tǒng) Bi 不在正域中則刪除該行end% 因為相同的條件得到的決策不一樣，e

4、nd%將在正域規(guī)則下約簡過的信息系統(tǒng)B 連接決策系統(tǒng) Dm,n=size(B);for i=1:mif(isequal(Bi,'')Bi1,k+1=Di1;end end顯示約簡決策% 系統(tǒng)disp('約簡后的決策系統(tǒng)為：)；m,n=size(B); for i=1:mdisp(Bi)；end%my_reduct 函數(shù)實現(xiàn)function C,num,reduct_attr=my_reduct(X)%y 為約簡后的 cell 數(shù)組， reduct_attr 為可約去的屬性%X 為行向量(元素為胞元)clc；% 約簡 m,n=size(X)；p,k=size(X1)；in

5、d_AT=ind(X)；% 尋找不可等價關(guān)系reduct_attr=； % 可約去的的屬性num=zeros(m,1);%約簡后的信息對應(yīng)的個體for i=1:kB=delete_AT(X,i);if(isequal(ind_A T,ind(B) %若 IND(AT-a=IND(A T) reduct_attr=union(reduct_attr,i);%則尋找到可約去的屬性X=B;endend% 剔除重復(fù)的行k=1;for i=1:mif(isequal(ind_A Ti,)%返回約簡后的信息系統(tǒng)C_i=ind_A Ti,1(1); num(k)=i;Ck,1=XC_i;k=k+1;ende

6、nd%ind 函數(shù)實現(xiàn)function yy=ind(X) % 尋找不可分辨關(guān)系 m,n=size(X);k=1;ind_A T=cell(m,1);for i=1:mfor j=(i+1):m %潛在問題，如 i=m 是終止循環(huán)，此時若最后一行不為空的話，將漏掃 if(isequal(Xi,'')%若 Xi 不為空ind_A Tk=union(ind_A Tk,i); %不可等價關(guān)系賦初值 if(isequal(Xi,Xj)Xj=''%若 Xi=Xj, 則刪除 Xjind_ATk=union(ind_A Tk,j); % 尋找不可等價關(guān)系 endendendk

7、=k+1;end%delete_AT 函數(shù)的源代碼function y=delete_AT(X,ATi) % 刪除 X 中第 i 列的屬性值 m,n=size(X);l,k=size(X1);for i=1:mXiA Ti=''endy=X;%pos 函數(shù)實現(xiàn)function pos_d=pos(X,D)% 求決策系統(tǒng)的正域函數(shù) %X 為條件屬性， D 為決策屬性 ind_D=ind(D); % 求決策屬性 D 的不可等價關(guān)系 m,n=size(ind_D);ind_X=ind(X); % 求信息系統(tǒng)屬性 X 的不可等價關(guān)系 low=; % 存儲正域個體的編號for i=1:mfor j=1:m if(isequal(ind_Xi,)&&isequal(ind_Dj,