7 参考文献
参考文献:
[1] 白其峥,数学建模案例分析[M].北京:海洋出版社,1999.
[2] 阮沈勇,王永利,桑群芳,MATLAB程序设计[M].北京:电子工业出版社,2004. [3] 杨纶标,模糊数学原理及应用[M].武汉:华南理工大学出版社,1998. [4] 姜启源,谢金星,叶俊,数学模型(第三版),北京:高等教育出版社,2003.
[5] 张莉,陈利琼,彭云飞,研究生录取多级综合测评数学模型[J].贵州大学学报(自然科学版),2010(2). [6] 钱颂迪,运筹学(修订版)[M].北京:清华大学出版社,1999. [7] 薛毅,数学建模基础[M].北京:北京工业大学出版社,2004.
8 附录
8.1 程序算法 A=1; B=0.9126; C=0.8; D=0.5245;
A0=[A B B B B;B A A B B;A B A C A;B B A B A;B B B B A;A A A C B;B B B A C;B A B A B;B A A C B;A A B B B];
A1=[A B A B A;A C A B B;B A C D C;B A B B A;B A B C A;B B A B D;A B C B B;B A A C A;B B A B B;D B A C C;D B C A B;A B A C A;B C B A D;D C A B C;A B C B B];
A2=[B A A B A;A B A C B;B A A C B;A B D C B;B A B C B;B D A B C;A B C B D;B A A C B;B B A B B;D B A C C;B C B A B;A B B A A;B A B C B;B B A A A;A B C B B];
A3=[A B B B B;A B A B A;B A D C B;A C B B B;B A B C B;A D A B A;A B C B B;B B A C B;A B C B B;C B A C D;D B B A B;A B B C B;B D C A C;D C A B C;A B B B A];
A4=[A B B B A;A B A B A;B A B B B;B B B B B;B B C B B;B C C B C;A C B D C;B B A C B;A C A B B;C A B B D;D B C B C;A B B B B;A B B A A;B B A B B;A C B B B];
A5=[A B A A A;B B A B A;B A B B A;B B B C B;A B B B A;C A B B C;A C B D C;B B A B A;B C A B B;B C A B D;C C A B C;A B B A B;A B B C D;B C B C C;A B B B B];
A6=[A B B B B;A B A B A;A B B B B;B A C D C;A A B B A;B C A B B;B A B C D;C B A B B;B C B B C;C A B C D;C B B A B;B B B B B;B B C B B;A C A B B;B B C A B];
A7=[B C A A B;A A B B A;A B C A C;B B B B A;A B C B B;A B B B A;A D B B B;B C B C C;A C A B B;C A B C D;D B A B C;A B B B B;B C A B C;C A B B D;B A A A A];
A8=[A A B B A;B B A A A;A B C B B;B C A B B;A C C B C;C A B B D;A B B C B;B C C B C;A B B B B;C B B C C;C B B B B;B A B B B;B B B A A;C B B C A;B B B A B]; rij=A1+A2+A3+A4+A5+A6+A7+A8; Bjk=rij/8; Bkj=Bjk'; Rj=sum(Bkj)/5
a=[416 410 405 397 392 389 385 382 380 378 377 372 360 358 356]; Aj=(a-min(a))/(max(a)-min(a)); n=0.4 m=0.6
Cj=n*Aj+m*Rj; b=sort(Cj);
[b,in]=sort(Cj);
11
d=fliplr(in) A=5; B=4; C=3; D=2;
F(:,:,1)=A1; F(:,:,2)=A2; F(:,:,3)=A3; F(:,:,4)=A4; F(:,:,5)=A5; F(:,:,6)=A6; F(:,:,7)=A7; F(:,:,8)=A8; for k=1:8
for i=1:10
E(i,:,k)=F(d(i),:,k); end end
for k=1:8
for i=1:10 for j=1:10
H(i,j,k)=sum(((E(i,:,k)-A0(j,:)+4))')/5; end end end
G=zeros(10,10); for k=1:8
G=G+H(:,:,k); end R=G/8
oo=2.4944; pp=0.8413; ff=0.1787; hh=0.6523; for i=1:10
for j=1:10 if R(i,j)<=4
RR(i,j)=(1+oo*(R(i,j)-pp)^(-2))^(-1); else
RR(i,j)=ff*log(R(i,j))+hh; end end end Sij=RR
I=[15 8 2 2;32 4 1 1;36 9 3 3;12 2 3 3;17 1 1 1;21 3 0 6;18 2 1 1;10 2 0 1;32 6 2 4;17 5 2 2]
12
I0=zeros(10,5); for j=1:4
for i=1:10
II(i,j)=(I(i,j)-min(I(:,j)))/(max(I(:,j))-min(I(:,j))); end end
K=sum(II')/4
P1=[2 3 1 4 3 3 2 1 4 3;3 1 2 3 2 4 4 3 1 4] P3=P1'
P0=[1 1 1 2 2 3 3 3 4 4] for k=1:2
for i=1:10 for j=1:10
P2(i,j,k)=P3(i,k)-P0(1,j); if (P2(i,j,k)==0)&&(k==1) PP(i,j,k)=1;
else if (P2(i,j,k)==0)&&(k==2) PP(i,j,k)=0.6309; else
PP(i,j,k)=0; end end end end end
P3=PP(:,:,1)+PP(:,:,2) for i=1:10
Tij(i,:)=P3(i,:).*K; end
CIj=sqrt(Sij.*Tij) Cij=CIJ'
8.2 LINGO 计算指派模型
13