求模糊相似矩阵的 MATLAB 程序

i）求模糊相似矩阵的MA TLAB 程序

a=[276.2 324.5 158.6 412.5 292.8 258.4 334.1 303.2 292.9 243.2 159.7 331.2

251.5 287.3 349.5 297.4 227.8 453.6 321.5 451.0 466.2 307.5 421.1 455.1

192.7 433.2 289.9 366.3 466.2 239.1 357.4 219.7 245.7 411.1 357.0 353.2

246.2 232.4 243.7 372.5 460.4 158.9 298.7 314.5 256.6 327.0 296.5 423.0

291.7 311.0 502.4 254.0 245.6 324.8 401.0 266.5 251.3 289.9 255.4 362.1

466.5 158.9 223.5 425.1 251.4 321.0 315.4 317.4 246.2 277.5 304.2 410.7

258.6 327.4 432.1 403.9 256.6 282.9 389.7 413.2 466.5 199.3 282.1 387.6

453.4 365.5 357.6 258.1 278.8 467.2 355.2 228.5 453.6 315.6 456.3 407.2

158.2 271.0 410.2 344.2 250.0 360.7 376.4 179.4 159.2 342.4 331.2 377.7

324.8 406.5 235.7 288.8 192.6 284.9 290.5 343.7 283.4 281.2 243.7 411.1];

mu=mean(a),sigma=std(a)

for i=1:12

for j=1:12

r(i,j)=exp(-(mu(j)-mu(i))^2/(sigma(i)+sigma(j))^2);

end

end

r

save data1 r a

ii）矩阵合成的MA TLAB函数

function rhat=hecheng(r);

n=length(r);

for i=1:n

for j=1:n

rhat(i,j)=max(min([r(i,:);r(:,j)']));

end

end

iii）求模糊等价矩阵和聚类的程序

load data1

r1=hecheng(r)

r2=hecheng(r1)

r3=hecheng(r2)

bh=zeros(12);

bh(find(r2>0.998))=1

iv)计算表6的程序

编写计算误差平方和的函数如下：

function err=wucha(a,t);

b=a;b(:,t)=[];

mu1=mean(a,2);mu2=mean(b,2);

err=sum((mu1-mu2).^2);

模糊算子程序：max（a^b）

function ab=synt(a,b)

m=size(a,1);n=size(b,2);

for i=1:m

for j=1:n

ab(i,j)=max(min([a(i,:);b(:,j)']));

end

end