附录
附件1：>> A=[1 3 5;1/3 1 3;1/5 1/3 1] A =
1.0000 3.0000 5.0000 0.3333 1.0000 3.0000 0.2000 0.3333 1.0000
>> [x,y]=eig(A) x =
0.9161 0.9161 0.9161 0.3715 -0.1857 + 0.3217i -0.1857 - 0.3217i 0.1506 -0.0753 - 0.1304i -0.0753 + 0.1304i y =
3.0385 0 0 0 -0.0193 + 0.3415i 0 0 0 -0.0193 - 0.3415i
>> w=x(:,1)/sum(x(:,1)) w =
0.6370 0.2583 0.1047
（?=3.0385，权重向量为w=（0.6370，0.2583，0.1047））
CI???nn?1=(3.0385-3)/(3-1)=0.01925
RI=0.58
CR?CIRI=0.01925/0.58=0.033<0.1,表明A通过了一致性验证。
>> B1=[1 3 5;1/3 1 3;1/5 1/3 1]
16
B1=
1.0000 3.0000 5.0000 0.3333 1.0000 3.0000 0.2000 0.3333 1.0000
>> [x,y]=eig(B1) x =
0.9161 0.9161 0.9161 0.3715 -0.1857 + 0.3217i -0.1857 - 0.3217i 0.1506 -0.0753 - 0.1304i -0.0753 + 0.1304i y =
3.0385 0 0 0 -0.0193 + 0.3415i 0 0 0 -0.0193 - 0.3415i
>> w=x(:,1)/sum(x(:,1)) w =
0.6370 0.2583 0.1047
（?1=3.0385，权重向量w1=（0.6370，0.2583，0.1047））
>> B2=[1 3 7;1/3 1 5;1/7 1/5 1] B2 =
1.0000 3.0000 7.0000 0.3333 1.0000 5.0000 0.1429 0.2000 1.0000
>> [x,y]=eig(B2) x =
0.9140 0.9140 0.9140
17
0.3928 -0.1964 + 0.3402i -0.1964 - 0.3402i 0.1013 -0.0506 - 0.0877i -0.0506 + 0.0877i y =
3.0649 0 0 0 -0.0324 + 0.4448i 0 0 0 -0.0324 - 0.4448i
>> w2=x(:,1)/sum(x(:,1)) w2 =
0.6491 0.2790 0.0719
（?2=3.0649，权重向量w2=（0.6491，0.2790，0.0719））
附件2：
>>
f=0.037*(9/44*0.037+2/11*0.038+3/22*0.039+15/88*0.04+13/88*0.041+7/44*0.042)/((9/44*0.037^2+2/11*0.038^2+3/22*0.039^2+15/88*0.04^2+13/88*0.041^2+7/44*0.042^2)+440*(0.01/1.64)^2)
f =
0.0813 >>
f=0.038*(9/44*0.037+2/11*0.038+3/22*0.039+15/88*0.04+13/88*0.041+7/44*0.042)/((9/44*0.037^2+2/11*0.038^2+3/22*0.039^2+15/88*0.04^2+13/88*0.041^2+7/44*0.042^2)+440*(0.01/1.64)^2)
f =
0.0835 >>
18
f=0.039*(9/44*0.037+2/11*0.038+3/22*0.039+15/88*0.04+13/88*0.041+7/44*0.042)/((9/44*0.037^2+2/11*0.038^2+3/22*0.039^2+15/88*0.04^2+13/88*0.041^2+7/44*0.042^2)+440*(0.01/1.64)^2)
f =
0.0857 >>
f=0.04*(9/44*0.037+2/11*0.038+3/22*0.039+15/88*0.04+13/88*0.041+7/44*0.042)/((9/44*0.037^2+2/11*0.038^2+3/22*0.039^2+15/88*0.04^2+13/88*0.041^2+7/44*0.042^2)+440*(0.01/1.64)^2)
f =
0.0879 >>
f=0.041*(9/44*0.037+2/11*0.038+3/22*0.039+15/88*0.04+13/88*0.041+7/44*0.042)/((9/44*0.037^2+2/11*0.038^2+3/22*0.039^2+15/88*0.04^2+13/88*0.041^2+7/44*0.042^2)+440*(0.01/1.64)^2)
f =
0.0901 >>
f=0.042*(9/44*0.037+2/11*0.038+3/22*0.039+15/88*0.04+13/88*0.041+7/44*0.042)/((9/44*0.037^2+2/11*0.038^2+3/22*0.039^2+15/88*0.04^2+13/88*0.041^2+7/44*0.042^2)+440*(0.01/1.64)^2)
f =
0.0923
?f1，f2，f3,f4，f5，f6???0.0813,0.0835,0.0857,0.0879,0.0901,0.0923?
90，80，60，75，65，70
7.317，6.68，5.142，6.5925，5.8565，6.461=38.327 n=8+7+5+7+6+6=39
19
假设t=1.5,p=60,Pmax=6000，Tmax=176,
?Max=11
附件3：
c=[-39]; A=[39;39;1]; b=[117.3;100;11]; Aeq=[]; beq=[]; vlb=[0;0;0]; vub=[];
[x,fval]=linprog(c,A,b,Aeq,beq,vlb,vub)
x =
2.5641
fval =
-100.0000
>> w=[1 1 3 3 5 5 6 7 1 9 9; 1 1 3 3 5 5 6 7 1 9 9;
1/3 1/3 1 1 5/3 5/3 2 7/3 1/3 3 3; 1/3 1/3 1 1 5/3 5/3 2 7/3 1/3 3 3;
1/5 1/5 3/5 3/5 1 1 6/5 7/5 1/5 9/5 9/5; 1/5 1/5 3/5 3/5 1 1 6/5 7/5 1/5 9/5 9/5; 1/6 1/6 1/2 1/2 5/6 5/6 1 7/6 1/6 3/2 3/2; 1/7 1/7 3/7 3/7 5/7 5/7 6/7 1 1/7 9/7 9/7; 1 1 3 3 5 5 6 7 1 9 9;
1/9 1/9 1/3 1/3 5/9 5/9 2/3 7/9 1/9 1 1; 1/9 1/9 1/3 1/3 5/9 5/9 2/3 7/9 1/9 1 1;] w =
1.0000 1.0000 3.0000 3.0000 5.0000 1.0000 9.0000 9.0000
1.0000 1.0000 3.0000 3.0000 5.0000 1.0000 9.0000 9.0000
0.3333 0.3333 1.0000 1.0000 1.6667 0.3333 3.0000 3.0000
0.3333 0.3333 1.0000 1.0000 1.6667 0.3333 3.0000 3.0000
0.2000 0.2000 0.6000 0.6000
1.0000
20
5.0000 6.0000 5.0000 6.0000 1.6667 2.0000 1.6667 2.0000 1.0000 1.2000 7.0000 7.0000 2.3333 2.3333 1.4000