实验报告：模糊统计与模糊分布.docVIP

[姓名][学号]

实验一模糊统计与模糊分布

实验目的：

会使用数学软件MATLAB进行模糊统计，掌握得到模糊分布(隶属函数)近似解析表达式的方法.

实验学时：4学时

实验内容：

⑴根据模糊集的数据统计区间(例如随机产生“青年人”的年龄区间)，绘制模糊分布图.

⑵根据上述数据，导出该模糊集隶属函数的近似解析表达式，说明你的方法的可行性(在同一坐标系中绘制模糊分布图与模糊集隶属函数的图形，给出残差平方和，用F分布检验).

实验日期：200x年x月x日

操作步骤：

⑴随机产生模糊集的数据统计区间，论域[a,b]，绘制模糊分布图.

n=20;a=12;b=40;A=rand(2,n);

fori=1:1:2*n

A(i)=a+A(i)*(b-a);%取整数A(i)=a+round(A(i)*(b-a));

end

fori=1:1:n

ifA(1,i)A(2,i)%调整区间的上下限

xxxxx=A(1,i);

A(1,i)=A(2,i);

A(2,i)=xxxxx;

end

end

A%显示模糊集A的数据

A1=b;A2=a;

fori=1:1:n

if(A(1,i)A1)A1=A(1,i);%A的最小值

end

if(A(2,i)A2)A2=A(2,i);%A的最大值

end

end

x=A1:(A2-A1)/100:A2;

Ax=[];

fori=1:1:length(x)

Ax(i)=0;

forj=1:n

if(x(i)=A(1,j)x(i)=A(2,j))Ax(i)=Ax(i)+1;

end

end

Ax(i)=Ax(i)/n;

end

plot(x,Ax)

结果：

⑵模糊集A隶属函数的近似解析表达式，在同一坐标系中绘制模糊分布图与模糊集隶属函数的图形，并给出残差平方和，用F分布检验.

以下仅给出多项式回归并用F分布检验的方法，供参考

x0=37.0:0.5:43.0;

y0=[3.403.003.002.272.101.831.531.701.801.902.352.542.90];

x=x0;%转置

y=y0;%转置

m=2;n=length(x);

xs=polyfit(x,y,m)%方法1多项式系数行向量降幂排列

yy=0;jxj=[];

fori=1:n

jxj(i)=xs(m+1);

fork=1:m

jxj(i)=jxj(i)+xs(k)*x(i)^(m-k+1);%解析解

end

yy=yy+y(i);

end

jxj

plot(x,jxj,x,y,-+)

X(n,m+1)=1;%定义n行m+1列矩阵X

X(:,1)=1;%X的第1列全为1

fori=1:m

X(:,i+1)=x.^i;%X的第i+1列=x的i次幂

end

b=inv(X*X)*X*y%方法2多项式系数列向量升幂排列

yy=yy/n;%计算y的均值

Syy=0;

fori=1:n

Syy=Syy+(y(i)-yy)*(y(i)-yy);%计算总和

end

Q=(y-X*b)*(y-X*b)%计算残差平方和

Ftjl=((Syy-Q)/m)/(Q/(n-m-1))%计算F统计量

Ffbz=Ftjl/10;

whilefcdf(Ffbz,m-1,n-m-1)0.99

Ffbz=Ffbz+0.1;%查表F分布值

end

Ffbz