mcm.m

mcm-美国2007年数学建模大赛的题目

function f=MCM(N)
%0 seat_interference
%2 aisle_interference
seat_sum1=0; aisle_sum1=0; all_sum1=0;
for n=1:N
a=rand(26,7);  %generate the random matrix
for i=1:26
	a(i,4)=1;      %the 4th col is the aisle
end
%tercher liu's method%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%1
for i=1:3
    for j=2
        if a(i,j)>a(i,j+1)
            a(i,j)=0;
        end
    end
end
for i=1:3
    for j=5
        if a(i,j)<a(i,j+1)
            a(i,j+1)=0;
        end
    end
end
%2
for i=7:25
    if a(i,1)>a(i+1,1)
        a(i,1)=2;
    end
end
%3
for i=4:6
    for j=1:2
        if a(i,j)>a(i,j+1)
            a(i,j)=0;
        end 
    end 
end
for i=4:6
    for j=5:6
        if a(i,j)<a(i,j+1)
            a(i,j+1)=0;
        end
    end
end
%4
for i=10:25
    if a(i,2)>a(i+1,2)
        a(i,2)=2;
    end
end
%5
for i=7:9
    for j=2
        if a(i,j)>a(i,j+1)
            a(i,j)=0;
        end
    end
end
for i=7:9
    for j=5:6
        if a(i,j)<a(i,j+1)
            a(i,j+1)=0;
        end
    end
end
%6
for i=13:25
    if a(i,3)>a(i+1,3)
        a(i,3)=2;
    end
end
%7
for i=10:12
    if a(i,3)>a(i+1,3)
        a(i,3)=2;
    end
end
for i=10:12
    for j=5:6
        if a(i,j)<a(i,j+1)
            a(i,j+1)=0;
        end
    end
end
%8
for i=15:25
    if a(i,7)>a(i+1,7)
        a(i,7)=2;
    end
end
%9
for i=13:14
    for j=5:6
        if a(i,j)<a(i,j+1)
            a(i,j+1)=0;
        end
    end
end
%10
for i=17:25
    if a(i,6)>a(i+1,6)
        a(i,6)=2;
    end
end
%11
for i=15:16
    for j=5
        if a(i,j)<a(i,j+1)
            a(i,j+1)=0;
        end
    end
end
%12
for i=19:25
    if a(i,5)>a(i+1,5)
        a(i,5)=2;
    end
end
%13
for i=17
    if a(i,5)>a(i+1,5)
        a(i,5)=2;
    end
end
seat_interference1=0;
aisle_interference1=0;
for i=1:26
    for j=1:7
        if a(i,j)==0
            seat_interference1=seat_interference1+1;
        end
    end
end
for i=1:26
    for j=1:7
        if a(i,j)==2
            aisle_interference1=aisle_interference1+1;
        end
    end
end
aisle_interference1=0;
seat_interference1;
all_interference1=aisle_interference1+seat_interference1;

               
seat_sum1=seat_sum1+seat_interference1 ;
aisle_sum1=aisle_sum1+aisle_interference1 ;
all_sum1=all_sum1+all_interference1 ;


end
seat_average1=seat_sum1/N;
aisle_average1=aisle_sum1/N;
all_average1=all_sum1/N;


X=[seat_average1 aisle_average1 all_average1]