初始化种群的函数：

function pop = initpop(popsize, piecesize)
% 初始化种群
% popsize           input  种群规模
% piecesize         input  工件数量
% pop               output 种群
pop = zeros(popsize, piecesize);
for i = 1:popsize
    pop(i, :) = randperm(piecesize);
end
end

计算适应度值的函数：

function fitness = calfitness(objvalue)
% 计算适应度值
% objvalue      input  目标函数值
% fitness       output 适应度值
fitness = 1 ./ objvalue;
end

制作加工流程矩阵的函数，该矩阵行代表机器，按照工序进行排序（假设工序1有两个机器，则第1行为工序1机器1，第2行为工序1机器2，第3行为工序2机器1……），列代表时间：

function gantt = makegantt(ptr, per, equsize)
% 制作加工流程矩阵
% ptr       input  工件时间记录
% per       input  工件设备记录
% equsize   input  工序设备数量
% gantt     output 加工流程矩阵
finaltime = max(max(ptr));
[piecesize, prosize] = size(per);
cumsumequ = cumsum(equsize);
gantt = zeros(sum(equsize), finaltime);
for pro = 1:prosize
    for i = 1:piecesize
        if pro == 1
            equ = per(i, pro);
        else
            equ = cumsumequ(pro - 1) + per(i, pro);
        end
        starttime = ptr(i, pro*2-1);
        endtime = ptr(i, pro*2);
        gantt(equ, starttime:endtime) = i; 
    end
end
end

生成两个随机数的函数：

function [rl, ru] = makerlru(maxnum)
% 制作两个随机整数， rl < ru
% maxnum        input  最大数
% rl            output 小随机数
% ru            output 大随机数
r1 = randi(maxnum);
r2  =randi(maxnum);
while r2 == r1
    r2 = randi(maxnum);
end
rl = min([r1, r2]);
ru = max([r1, r2]);
end