关于粗糙集的matlab计算

描述了matlab如何计算粗糙集的条件属性的权重的方法

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%读取信息系统文件
file = textread（‘data2.txt‘‘%s‘‘delimiter‘‘\n‘‘whitespace‘‘‘）;
[mn]=size（file）;
for i=1:m
    words=strread（file{i}‘%s‘‘delimiter‘‘ ‘）;
    words=words‘;
    X{i}=words;    
end
X=X‘;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%读取决策系统文件
file = textread（‘decision2.txt‘‘%s‘‘delimiter‘‘\n‘‘whitespace‘‘‘）;
[mn]=size（file）;
for i=1:m
    words=strread（file{i}‘%s‘‘delimiter‘‘ ‘）;
    words=words‘;
    D{i}=words;    
end
D=D‘;

[BAT]=ind_reduct（X）;       %信息系统的约简（基于不可等价关系约简）
ind_AT=ind（X）;     %求信息系统的不可等价关系
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%显示约简信息系统
disp（‘约简后的信息系统为：‘）;
[mn]=size（B）;
for i=1:m
    disp（B{i}）;
end	
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%决策系统的正域约简
[lk]=size（X{1}）;
pos_d=pos（XD）;%求正域
[Breduct_attr]=pos_reduct（XD）;%基于正域启发式函数对决策系统约简

%将在正域规则下约简过的信息系统B连接决策规则D，使B之变成一决策信息系统
[mn]=size（B）;
for i=1:m
    if（~isequal（B{i}[]））
        B{i}{1k+1}=D{i}{1};
    end
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%显示约简决策系统
disp（‘约简后的决策系统为：‘）;
[mn]=size（B）;
for i=1:m
    disp（B{i}）;
end        


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%决策系统分类
ind_D=ind（D）;
[mn]=size（ind_D）;
for i=1:m
    if （~isequal（ind_D{i}[]））
        Low{i}=low（ind_ATind_D{i}）;%求决策规则D在基于不可等价关系ind_AT下的下近似
        Upp{i}=upp（ind_ATind_D{i}）;%求决策规则D在基于不可等价关系ind_AT下的上近似
    end
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%显示确定性分类结果
disp（‘约简后的决策系统其确定性规则有：‘）;
[mn]=size（Low）;
certern=[];
for i=1:n
    certern=union（certernLow{i}）;
end
disp（certern）

disp（‘约简后的决策系统其不确定性规则有：‘）;
[mn]=size（Upp）;
uncertern=[];
for i=1:n
    uncertern=union（uncerternsetdiff（Upp{i}Low{i}））;
end
disp（uncertern）%显示不确定规则

disp（‘其不确定性规则的可信度为：‘）;
if（~isempty（uncertern））
    Cer=cer（uncerternind_AT）;
    disp（Cer）  %显示可行度
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%分别求出信息系统和决策系统的核属性
disp（‘信息系统的核属性为：‘）;
core_ind=core_ind_reduct（X）;
disp（core_ind）
disp（‘决策系统的核属性为：‘）;
core_pos=core_pos_reduct（XD）;
disp（core_pos）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2364  2013-01-14 10:02  cucaoji.m

----------- ---------  ---------- -----  ----

                 2364                    1