模糊数学MATLAB实验

《模糊数学MATLAB实验》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、模糊数学MATLAB实验例1：程序：1）a=[0.40.70;10.80.5],b=[10.7;0.40.6;00.3]m=size(a,1);n=size(b,2);fori=1:mforj=1:nc(i,j)=max(min([a(i,:);b(:,j)']));endendc2）a=[0.40.70;10.80.5],b=[10.7;0.40.6;00.3]m=size(b,1);n=size(a,2);fori=1:mforj=1:nc(i,j)=max(min([b(i,:);a(:,j)']));endendc3)a*b4)b*a例2现有五个等级的茶叶样品，待识别茶叶B。反映茶叶质

2、量的因素有六项指标，构成论域U，其中U={（条索），（色泽），（净度），（汤色），（色泽），（滋味）}设五个等级的样品对6项指标的数值为：A1=(0.5,0.4,0.3,0.6,0.5,0.4)A2=(0.3,0.2,0.2,0.1,0.2,0.2)A3=(0.2,0.2,0.2,0.1,0.1,0.2)A4=(0,0.1,0.2,0.1,0.1,0.1)A5=(0,0.1,0.1,0.1,0.1,0.1)待识别茶叶的各项指标值为B=(0.4,0.2,0.1,0.4,0.5,0.6)确定B的属类。程序a=[0.50.40.30.60.50.4;0.30.20.20.10.20.2;0.20.2

3、0.20.10.10.2;00.10.20.10.10.1;00.10.10.10.10.1];b=[0.40.20.10.40.50.6];fori=1:5x=[a(i,:);b];t(i)=min([max(min(x))1-min(max(x))]);endt例3：用MATLAB编程计算这个模糊关系矩阵的截矩阵程序A=[10.50.20;0.510.10.3;0.20.110.8;00.30.81];fori=1:4forj=1:4ifA(i,j)>=0.8A(i,j)=1;elseA(i,j)=0;endendendA例4.求解的MATLAB程序如下：i）求模糊相似矩阵的MATLAB程

4、序a=[276.2324.5158.6412.5292.8258.4334.1303.2292.9243.2159.7331.2;251.5287.3349.5297.4227.8453.6321.5451.0466.2307.5421.1455.1;192.7433.2289.9366.3466.2239.1357.4219.7245.7411.1357.0353.2;246.2232.4243.7372.5460.4158.9298.7314.5256.6327.0296.5423.0;291.7311.0502.4254.0245.6324.8401.0266.5251.3289.925

5、5.4362.1;466.5158.9223.5425.1251.4321.0315.4317.4246.2277.5304.2410.7;258.6327.4432.1403.9256.6282.9389.7413.2466.5199.3282.1387.6;453.4365.5357.6258.1278.8467.2355.2228.5453.6315.6456.3407.2;158.2271.0410.2344.2250.0360.7376.4179.4159.2342.4331.2377.7;324.8406.5235.7288.8192.6284.9290.5343.7283.428

6、1.2243.7411.1];mu=mean(a),sigma=std(a)fori=1:12forj=1:12r(i,j)=exp(-(mu(j)-mu(i))^2/(sigma(i)+sigma(j))^2);endendrsavedata1raii）矩阵合成的MATLAB函数functionrhat=hecheng(r);n=length(r);fori=1:nforj=1:nrhat(i,j)=max(min([r(i,:);r(:,j)']));endendiii）求模糊等价矩阵和聚类的程序loaddata1r1=hecheng(r)r2=hecheng(r1)r3=hecheng(

7、r2)bh=zeros(12);bh(find(r2>0.998))=1iv)计算表6的程序编写计算误差平方和的函数如下：functionerr=wucha(a,t);b=a;b(:,t)=[];mu1=mean(a,2);mu2=mean(b,2);err=sum((mu1-mu2).^2);计算28个方案的主程序如下：loaddata1ind1=[1,5];ind2=[2:3,6,8:11];