MATLAB实现FSP调度的最大完工时间求解

编写FSP问题求总完工时间函数

函数形式：function [Cmax] = FSP (n,m)，其中n表示工件个数，m表示机器个数。

函数功能：在函数中随机生成加工时间矩阵T（推荐取值为1到9），随机生成一个调度解，然后在该加工时间矩阵下，求取该调度解的总完工时间，以Cmax输出。

function [ Cmax ] = FSP( n,m )
%FSP  求n个工件在m台机器的完工时间
%   输入工件数n和机器数m，随机生成（n,m）n个工件在m台机器上的加工矩阵，然后随机生成调度解并求出调度解下的完工时间
T=randi(9,n,m);  %生成1~9加工时间矩阵
fprintf('随机生成%d个工件在%d台机器的加工时间矩阵:\n',n,m);
disp(T);  
A=randperm(n);   %随机生成n个工件加工调度解，存储到A
disp('加工顺序为：')
disp(A);
for i=1:n     %将加工时间按调度解排序赋值
    for j=1:m
        S(i,j)=T(A(1,i),j); %将每个工件在机器上的加工时间按调度解排序
    end
end
disp('按加工顺序调整后的加工矩阵：')
disp(S);
C(1,1)=S(1,1); %C为完工时间矩阵，初始赋值为S（1,1）
for j=2:m
    C(1,j)=C(1,j-1)+S(1,j); %计算第1个工件在1~m台机器上的加工完工时间矩阵
end
for i=2:n
    C(i,1)=C(i