(1)建立一个M文件,导入数据进入矩阵
filename=‘data.xls’;
A=xlsread(filename);
(2)通过矩阵记录均值和标准差
u=A(:,1);
s=A(:,2);
e=zeros(16,1) ;//负责纪录夺冠次数
(3)将标准差和均值进行整体比较,将大的数据进行记录(在模型中相当于晋级的队伍)
If
normrnd(u(c(j,m)),s(c(j,m)))>normrnd(u(c(j+1,m)),s(c(j+1,m)))
c(k,m+1)=c(j,m);
else
c(k,m+1)=c(j+1,m);
(3)根据3的公式,将此比较方式进行循环,得出每两个相邻队伍获胜的一方,再进行轮次的循环,记为m,从1到4,(16支队伍四轮比赛)
for m=1:4
k=1;
for j= 1:2:(2^(5-m)-1)(m为第几轮
(3)循环
end
k=k+1;
end
(4) 将(3)的冠军求法循环一万次,就可以得出误差比较小的答案 将每次冠军记录在矩阵中
e(c(1,5))=e(c(1,5))+1;