cloudmodel.m

云模型的matlab实现~具体实例~~~~~~~~~~~~~~~~~~~~~~~~~~

U(1)=0;U(2)=0;U(3)=0;U(4)=0;
for i=1:1000
    e=0.165+0.01*randn(1);
    x(1)=7.81+e*randn(1);
    e=(9-x(1)+0.5)/3+0.01*randn(1);
    if x(1)<9
        y(1)=exp(-(x(1)-9)^2/(2*e*e));%指标1对于评语“优”的录属度
    else
        y(1)=1; 
    end
    e=0.066+0.01*randn(1);
    x(2)=9.12+e*randn(1);
    e=(9-x(2)+0.5)/3+0.01*randn(1);
    if x(2)<9
        y(2)=exp(-(x(2)-9)^2/(2*e*e));%指标2对于评语“优”的录属度
    else
        y(2)=1;
    end
    e=0.165+0.01*randn(1);
    x(3)=7.725+e*randn(1);
    e=(9-x(3)+0.5)/3+0.01*randn(1);
    if x(3)<9
        y(3)=exp(-(x(3)-9)^2/(2*e*e));%指标3对于评语“优”的录属度
    else
        y(3)=1;
    end
    e=0.066+0.01*randn(1);
    x(4)=8.22+e*randn(1);
    e=(9-x(4)+0.5)/3+0.01*randn(1);
    if x(4)<9
        y(4)=exp(-(x(4)-9)^2/(2*e*e));%指标4对于评语“优”的录属度
    else
        y(4)=1;
    end
    e=0.099+0.01*randn(1);
    x(5)=9.18+e*randn(1);
    e=(9-x(5)+0.5)/3+0.01*randn(1);
    if x(5)<9
        y(5)=exp(-(x(5)-9)^2/(2*e*e));%指标5对于评语“优”的录属度
    else
        y(5)=1;
    end
    SC(1)=0.3*y(1)+0.2*y(2)+0.3*y(3)+0.1*y(4)+0.1*y(5);%所有指标综合相对于评语“优”的录属度
    e=0.165+0.01*randn(1);
    x(1)=7.81+e*randn(1);
    if x(1)>9
        e=(x(1)-9+0.5)/3+0.01*randn(1);
        y(1)=exp(-(x(1)-9)^2/(2*e*e));%指标1对于评语“良”的录属度
    elseif x(1)<8
        e=(8-x(1)+0.5)/3+0.01*randn(1);
        y(1)=exp(-(x(1)-8)^2/(2*e*e));
    else
        y(1)=1;
    end      
    e=0.066+0.01*randn(1);
    x(2)=9.12+e*randn(1);
    if x(2)>9
        e=(x(2)-9+0.5)/3+0.01*randn(1);
        y(2)=exp(-(x(2)-9)^2/(2*e*e));%指标2对于评语“良”的录属度
    elseif x(2)<8
        e=(8-x(2)+0.5)/3+0.01*randn(1);
        y(2)=exp(-(x(2)-8)^2/(2*e*e));
    else
        y(2)=1;
    end
    e=0.165+0.01*randn(1);
    x(3)=7.725+e*randn(1);
    if x(3)>9
        e=(x(3)-9+0.5)/3+0.01*randn(1);
        y(3)=exp(-(x(3)-9)^2/(2*e*e));%指标3对于评语“良”的录属度
    elseif x(3)<8
        e=(8-x(3)+0.5)/3+0.01*randn(1);
        y(3)=exp(-(x(3)-8)^2/(2*e*e));
    else
        y(3)=1;
    end
    e=0.066+0.01*randn(1);
    x(4)=8.22+e*randn(1);
    if x(4)>9
        e=(x(4)-9+0.5)/3+0.01*randn(1);
        y(4)=exp(-(x(4)-9)^2/(2*e*e));%指标4对于评语“良”的录属度
    elseif x(4)<8
        e=(8-x(4)+0.5)/3+0.01*randn(1);
        y(4)=exp(-(x(4)-8)^2/(2*e*e));
    else
        y(4)=1;
    end
    e=0.099+0.01*randn(1);
    x(5)=9.18+e*randn(1);
    if x(5)>9
        e=(x(5)-9+0.5)/3+0.01*randn(1);
        y(5)=exp(-(x(5)-9)^2/(2*e*e));%指标5对于评语“良”的录属度
    elseif x(5)<8
        e=(8-x(5)+0.5)/3+0.01*randn(1);
        y(5)=exp(-(x(5)-8)^2/(2*e*e));
    else
        y(5)=1;
    end
    SC(2)=0.3*y(1)+0.2*y(2)+0.3*y(3)+0.1*y(4)+0.1*y(5);%所有指标综合相对于评语“良”的录属度
    e=0.165+0.01*randn(1);
    x(1)=7.81+e*randn(1);
    if x(1)<6
        e=(6-x(1)+0.5)/3+0.01*randn(1);
        y(1)=exp(-(x(1)-6)^2/(2*e*e));%指标1对于评语“合格”的录属度
    elseif x(1)>8
        e=(x(1)-8+0.5)/3+0.01*randn(1);
        y(1)=exp(-(x(1)-8)^2/(2*e*e));
    else
        y(1)=1;
    end
    e=0.066+0.01*randn(1);
    x(2)=9.12+e*randn(1);
    if x(2)<6
        e=(6-x(2)+0.5)/3+0.01*randn(1);
        y(2)=exp(-(x(2)-6)^2/(2*e*e));%指标2对于评语“合格”的录属度
    elseif x(2)>8
        e=(x(2)-8+0.5)/3+0.01*randn(1);
        y(2)=exp(-(x(2)-8)^2/(2*e*e));
    else
        y(2)=1;
    end
    e=0.165+0.01*randn(1);
    x(3)=7.725+e*randn(1);
    if x(3)<6
        e=(6-x(3)+0.5)/3+0.01*randn(1);
        y(3)=exp(-(x(3)-6)^2/(2*e*e));%指标3对于评语“合格”的录属度
    elseif x(3)>8
        e=(x(3)-8+0.5)/3+0.01*randn(1);
        y(3)=exp(-(x(3)-8)^2/(2*e*e));
    else
        y(3)=1;
    end
    e=0.066+0.01*randn(1);
    x(4)=8.22+e*randn(1);
    if x(4)<6
        e=(6-x(4)+0.5)/3+0.01*randn(1);
        y(4)=exp(-(x(4)-6)^2/(2*e*e));%指标4对于评语“合格”的录属度
    elseif x(4)>8
        e=(x(4)-8+0.5)/3+0.01*randn(1);
        y(4)=exp(-(x(4)-8)^2/(2*e*e));
    else
        y(4)=1;
    end
    e=0.099+0.01*randn(1);
    x(5)=9.18+e*randn(1);
    if x(5)<6
        e=(6-x(5)+0.5)/3+0.01*randn(1);
        y(5)=exp(-(x(5)-6)^2/(2*e*e));%指标5对于评语“合格”的录属度
    elseif x(5)>8
        e=(x(5)-8+0.5)/3+0.01*randn(1);
        y(5)=exp(-(x(5)-8)^2/(2*e*e));
    else
        y(5)=1;
    end
    SC(3)=0.3*y(1)+0.2*y(2)+0.3*y(3)+0.1*y(4)+0.1*y(5);%所有指标综合相对于评语“合格”的录属度
    e=0.165+0.01*randn(1);
    x(1)=7.81+e*randn(1);
    e=(x(1)-6+0.5)/3+0.01*randn(1);
    if x(1)>6
        y(1)=exp(-(x(1)-6)^2/(2*e*e));%指标1对于评语“不合格”的录属度
    else
        y(1)=1;
    end
    e=0.066+0.01*randn(1);
    x(2)=9.12+e*randn(1);
    e=(x(2)-6+0.5)/3+0.01*randn(1);
    if x(2)>6
        y(2)=exp(-(x(2)-6)^2/(2*e*e));%指标2对于评语“不合格”的录属度
    else
        y(2)=1;
    end
    e=0.165+0.01*randn(1);
    x(3)=7.725+e*randn(1);
    e=(x(3)-6+0.5)/3+0.01*randn(1);
    if x(3)>6
        y(3)=exp(-(x(3)-6)^2/(2*e*e));%指标3对于评语“不合格”的录属度
    else
        y(3)=1;
    end
    e=0.066+0.01*randn(1);
    x(4)=8.22+e*randn(1);
    e=(x(4)-6+0.5)/3+0.01*randn(1);
    if x(4)>6
        y(4)=exp(-(x(4)-6)^2/(2*e*e));%指标4对于评语“不合格”的录属度
    else
        y(4)=1;
    end
    e=0.099+0.01*randn(1);
    x(5)=9.18+e*randn(1);
    e=(x(5)-6+0.5)/3+0.01*randn(1);
    if x(5)>6
        y(5)=exp(-(x(5)-6)^2/(2*e*e));%指标5对于评语“不合格”的录属度
    else
        y(5)=1;
    end
    SC(4)=0.3*y(1)+0.2*y(2)+0.3*y(3)+0.1*y(4)+0.1*y(5);%所有指标综合相对于评语“不合格”的录属度
    U(1)=U(1)+SC(1);
    U(2)=U(2)+SC(2);
    U(3)=U(3)+SC(3);
    U(4)=U(4)+SC(4);
end
x=[1:4];
y=[U(1)/1000 U(3)/1000 U(3)/1000 U(4)/1000];
bar(x,y);