intercross1.m

改进GA解决TSP问题

function farm1=intercross1(selection,n,N,p)
c=p*n;
for i=1:floor(c/2)
    mi=ceil(rand*n);
    ni=ceil(rand*n);
    A(i,:)=selection(mi,:);
    B(i,:)=selection(ni,:);
    w1=ceil(rand*N);
    w2=ceil(rand*N);
    mid=zeros(2,N);
    co=1;
    while co
        if w1==w2
           co=1;
           w2=ceil(rand*N);
        else co=0;
        end
    end
    if w1>w2
       for mm=w2:w1
            mid(1,mm)=A(i,mm);
            mid(2,mm)=B(i,mm);
            A(i,mm)=mid(2,mm);
            B(i,mm)=mid(1,mm);
        end
        for mm1=1:w2-1
            for mm4=w2:w1
            for mm2=w2:w1
                if A(i,mm1)==mid(2,mm2)
                    A(i,mm1)=mid(1,mm2);
                end
                if B(i,mm1)==mid(1,mm2)
                   B(i,mm1)=mid(2,mm2);
                end
            end
            end
        end
        for mm3=w1+1:N
            for mm4=w2:w1
            for mm2=w2:w1
                if A(i,mm3)==mid(2,mm2)
                    A(i,mm3)=mid(1,mm2);
                end
                if B(i,mm3)==mid(1,mm2)
                   B(i,mm3)=mid(2,mm2);
                end
            end
            end
        end  
    end
    if w1<w2
        for mm=w1:w2
            mid(1,mm)=A(i,mm);
            mid(2,mm)=B(i,mm);
            A(i,mm)=mid(2,mm);
            B(i,mm)=mid(1,mm);
        end
        for mm1=1:w1-1
            for mm4=w1:w2
            for mm2=w1:w2
                if A(i,mm1)==mid(2,mm2)
                   A(i,mm1)=mid(1,mm2);
                end
                if B(i,mm1)==mid(1,mm2)
                   B(i,mm1)=mid(2,mm2);
                end
            end
            end
        end
        for mm3=w2+1:N
            for mm4=w1:w2
            for mm2=w1:w2
                if A(i,mm3)==mid(2,mm2)
                    A(i,mm3)=mid(1,mm2);
                end
                if B(i,mm3)==mid(1,mm2)
                   B(i,mm3)=mid(2,mm2);
                end
            end
            end
        end
    end
selection(mi,:)=A(i,:);
selection(ni,:)=B(i,:);
farm1=selection;
end