原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
%参数的输入
p二input (Please put into p
)
m二input (How many juries?
)
s二input (How many school s?
)
num二input(How many papers?
,)
1 二ones (s, 1);
%1数据从txt文件中读取
l=load(，d: \L txt);
tic;
%打开计时器
t二zeros (s+1, 1);
i0=l;
%初值变量的设定
%统计各个数据段的数据之和
t(2)=l(l);
for i18=3:s+1
t(il8)=t(il8-l)+l (118-1)
?
end
zh二zeros (p,1);
n二0;
for i21=l:s
n二n+1 (i21);
%定义存放不平衡度的矩阵%定义存放加权不平衡度的矩阵
%定义存放不平衡度的矩阵
%定义存放加权不平衡度的矩阵
%设定方差较大的初始不均衡度
%设定方差较大的初始加权不均衡度
varl=zeros(10000,1);
var3(1)=10000;
varl(1)=10000;
c=zeros (m, s);
w=l;
%用组合数求解总的组合数量
for i19=m：-l:m-p+l
w二il9*w/(il9~m+p);
end
%?产生总的组合数
for i1=1:m-p+l
for i2=il+l:m-p+2
for i3二i2+l:m-p+3
for i4=i3+l:m-p+4
a(i0, l:p) = [il, i2, 13 i4];
iO二i0+l;
end end
%循环一万次，得出较好的结果
for i25=2:10000
i0=l;
zz=zcros (m, s);
per=zeros (m, 1);
r=l;
var2=zeros (s,1);
for i5=l:s
for i6=l:1 (i5)
k=ceil (unifrnd (1, w)) ; %随机产生 1—w的整数
if (i5m+l)
wh订e((all (a(k, :)~=i5))二二0) %保证评委不评自己所在学校试卷
k=ceil(unifrnd (1,w));
end
b(r, 1 :p)=a (k, 1 :p);
r=r+l;
else
b(r, l:p)=a(k, l:p);
r=r+l; %满足条件则赋值给b
end
%统计每所学校分配给各个评委的试卷数 for i10=1:s
for il2=t(il0)+l:t(il0+l)
for il1=1:p
n=b(il2, ill);
zz (n, il0)=l+zz (n, ilO);;
end
end
end
%求均衡度和加权不均衡度 avel=num*p/m;
for il2二l:m
for il3=l:s
per (i 12)=zz(订2, il3)+per(il2);
end
end
for il4=l:m varl (i25)=varl (i25) + (per(il4)-ave1)*(per(il4)-avel); %fAt
y22i
end for i15=1:s
if(il5〈m+l)
ave2 (il5)=l (il5)/(m-1);
else
ave2 (il5)=l(il5)/m;
end
end
for il6=l:s
for i17=1:m
var2(il6)=var2 (il6) + (zz (il7, il6)-ave2 (il6))*(zz (il7, il6)-av e2(il6));
end
var2(il6)=var2(il6)*1(il6)/num;
end
for i18=1:s
var3(i25)=var3(i25)+var2(i18);
%取不均衡度与加权不均衡度都较小的方案
if ((varl (i25) varl (i25T))(var3 (i25) var3(i25-l)))
o=varl(i25);
u=var3(i25);
c=zz;
end
end
o
u
x=l:m;
bar (x, sum(c，))
grid on
tl二toe
实例运算结果
o 二
1786 %不均衡度
u 二
144.0735 %加权不均衡度
tl 二
266.0034
%总共运行时间为266. 0034s
%每个评委在给各个学校打分的试卷数（该模型中将没有评委的学校
放在后面
C =
Columns 1 through 16
0
1
0
2
1
5
3
2
1
0
0
3
0
1
2
0
1
0
4
3
3
2
2
3
2
2
0
2
1
1
0
1
1
1
0
4
2
1
2
1
0
3
3
0
3
2
2
2
0
0
2
0
1
4
3
0
4
2
1
2
2
0
0
0
0
1
3
3
0
文档评论（0）