load data41
ind1 = [1,5];ind2 = [2,3,6,8,9,10,11];ind3 = [4,7];so = [0];
for i =1:length(ind1)
for j = 1:length(ind3)
for k = 1:length(ind2)
t = [ind1(i),ind3(j),ind2(k)];
err = wucha(a,t);
so =[so,[t,err]];
end
end
end
tm = find(so(:,4)==min(so(:,4)));
shanchu = so(tm,1:3);
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
save data41 r a%保存数据r以及a在data41中
load data41
r1 = hecheng(r);
r2 = hecheng(r1);
r3 = hecheng(r2);
bh = zeros(12);
bh(r2>0.998) = 1;
function rhat = hecheng(r)
n = length(r);
rhat = zeros(n);
for i = 1:n
for j = 1:n
rhat(i,j) = max(min([r(i,:);r(:,j)']))
end
end
function err = wucha(a,t)
b = a;b(:,t) = [];
mu1 = mean(a,2);mu2 = mean(b,2);
err = sum((mu1-mu2).^2);
自己学习留下的记录(错误地方请斧正)