（Matlab）复现《改进NSGA_II算法求解考虑运输约束的柔性作业车间节能调度问题》中的变异操作

昨天的时候想把这篇文章的变异操作和上一篇文章的变异操作发个综合版本，但是在机器的选择问题上卡住了，今天对代码进行再一次调试，最终得到了想要的结果。

文章中的变异操作分为两个部分，一个是工序编码层变异，另外一个是机器编码层的变异。

工序编码层变异：

 代码实现：

for i=1:NIND               
    %取一个个体
    S=Chrom(i,:);
    number=randperm(WNumber);
    pos1=number(1);
    pos2=number(2);
    workpiece1=S(pos1);
    workpiece2=S(pos2);
    B=zeros(PNumber,2);
    C=zeros(PNumber,2);
    a=1;
    b=1;
    for q=1:WNumber
        if S(q)==workpiece1
            B(a,1)=q;
            B(a,2)=S(q+WNumber);
            a=a+1;
        end
        if S(q)==workpiece2
            C(b,1)=q;
            C(b,2)=S(q+WNumber);
            b=b+1;
        end
    end
    
    %% 是否变异
    if MUTR>rand