%賦初值
value1=180;
f11=1;
f21=1;
f31=1;
f41=1;
%通過循環求解
for i=1:1000;   %循環次數設為1000次
%四項成績的和
v0=[34.6	35	34.4	34.2	35.1	35	34.6	34.1	34.5	34.4];   
%每個運動員在每項比賽中的最差成績（最悲觀）
C0=[8.4	9.3	8.4	8.1	8.4	9.4	9.5	8.4	8.4	9                   
    8.4	8.4	8.1	8.7	9	8.7	8.4	8.8	8.4	8.1
    9.1	8.4	8.4	9	8.3	8.5	8.3	8.7	8.4	8.2
    8.7	8.9	9.5	8.4	9.4	8.4	8.4	8.2	9.3	9.1];
C1=C0;            %把C0賦值給C1，以保護C0的數據
v1=v0;            %把v0賦值給v1，以保護v0的數據

%通過“隨機數法”，不重復地取出C0和v0中對應的四列
f1=fix(10*rand(1,1))+1;    %從1~10十個數中隨機抽取一個f1
C1(:,f1)=[];               %清除C1的第f1列
val1=v1(1,f1);             %將v1的第f1項數值賦給val1
v1(:,f1)=[];               %清除v1的第f1項
f2=fix(9*rand(1,1))+1;     %從1~9九個數中隨機抽取一個f2
C1(:,f2)=[];               %清除C1的第f2列
val2=v1(1,f2);             %將v1的第f2項數值賦給val2
v1(:,f2)=[];               %清除v1的第f2項
f3=fix(8*rand(1,1))+1;     %從1~8八個數中隨機抽取一個f3
C1(:,f3)=[];               %清除C1的第f3列
val3=v1(1,f3);             %將v1的第f3項數值賦給val3
v1(:,f3)=[];               %清除v1的第f3項
f4=fix(7*rand(1,1))+1;     %從1~7七個數中隨機抽取一個f4
C1(:,f4)=[];               %清除C1的第f4列
val4=v1(1,f4);             %將v1的第f4項數值賦給val4
v1(:,f4)=[];               %清除v1的第f4項

C=-[C1(1,:),C1(2,:),C1(3,:),C1(4,:)];    %0-1規劃中目標函數的系數陣
a1=[ones(1,6),zeros(1,18)];
a2=[zeros(1,6),ones(1,6),zeros(1,12)];
a3=[zeros(1,12),ones(1,6),zeros(1,6)];
a4=[zeros(1,18),ones(1,6)];
b1=[1,zeros(1,5),1,zeros(1,5),1,zeros(1,5),1,zeros(1,5)];
b2=[0,1,zeros(1,4),0,1,zeros(1,4),0,1,zeros(1,4),0,1,zeros(1,4)];
b3=[0,0,1,zeros(1,3),0,0,1,zeros(1,3),0,0,1,zeros(1,3),0,0,1,zeros(1,3)];
b4=[0,0,0,1,0,0,0,0,0,1,0,0,0,0,0,1,0,0,0,0,0,1,0,0];
b5=[zeros(1,4),1,0,zeros(1,4),1,0,zeros(1,4),1,0,zeros(1,4),1,0];
b6=[zeros(1,5),1,zeros(1,5),1,zeros(1,5),1,zeros(1,5),1];
A=[b1;b2;b3;b4;b5;b6];                  %0-1規劃中不等號限制條件的系數陣
b=[3;3;3;3;3;3];                        %0-1規劃中不等號限制條件的右端向量
Aeq=[a1;a2;a3;a4];                      %0-1規劃中等號限制條件的系數陣
beq=[2;2;2;2];                          %0-1規劃中等號限制條件的右端向量
[x,val]=bintprog(C,A,b,Aeq,beq);        %使用matlab自帶的0-1規劃函數求解
value2=-val+val1+val2+val3+val4;        %在本次循環中求得的總成績值
%通過一個小的if循環將本次循環所得的總成績值與初始總成績值對比，
...取大的作為總成績初值賦給value1存儲，并將其對應的“取四個人”
...時產生的隨機數賦值給f11,f22,f33,f44作為存儲值
if value2>value1;
   value1=value2;
    f11=f1;
    f21=f2;
    f31=f3;
    f41=f4;
end
end
[value1,f11,f21,f31,f41]  %通過循環得到1000次選擇中的最大值及其方案