CPLEX 在matlab中实现的一个例子-RCPSP

CPLEX 在matlab中实现的一个例子-RCPSP

很多人运行不了或者不方便，于是我将去年写的代码发到网盘中，大家可以下载看看，我刚刚运行了一下，代码没问题。贴上链接和提取码
链接：https://pan.baidu.com/s/1AQtgE5RIBq8HeMjuxekJCw 提取码：3g71
本人去年想采用使用cplex求解数学模型，本想在博客上找一个例子仿照练习，无奈找遍全网都没有一个合适的，幸好有师兄们的代码。

matlab调用CPLEX的方法网上有很多，在这里我就不介绍了。直接上例子。这个代码以RCPSP为例，熟悉RCPSP的人肯定只要，如果不熟悉也没关系，我附上数学模型。

资源受限项目调度问题描述：

附上代码

% CPLEX  matlab上实现
% a cplex code of RCPSP on matlab
function solution=cplex_function()
% 定义的全部变量，储存项目信息
global Case 
% 假设的时间最大值
H=sum(Case.dur);
% 定义变量   bin是0 1 值  full代表是满的矩阵，非稀疏矩阵
x=binvar(Case.job_count,H,'full');
% 资源上限，是常数
r=Case.Con_R;
% 定义中间变量   int 是整数值，可以是任意整数
C=intvar(1);
% 目标函数
Objective=min(C); 

 % 约束
Constraints=[];
% 作业时间约束
% dur 代表作业持续时间
x2=sum(x,2);
for i=1:Case.job_count    
    Constraints=[Constraints,x2(i)==Case.dur(i)];
end
for j=1:Case.job_count 
    for i=1:H
        Constraints=[Constraints,x(j,i)*i<=C];
    end
end
    
% prec代表作业间约束 
% 紧前紧后约束
for j=2:Case.job_count
    Pjj=[];
    for jj=1:j
        if(jj~=j && Case.prec(jj,j)==1)
            Pjj=[Pjj jj];
        end
    end
    for index=1:length(Pjj) %对于任意一个j的紧前任务
        for d=1:H
            last=0;
            for dd=1:d-1
                last=last+x(Pjj(index),dd);
            end
            Constraints=[Constraints,Case.dur(Pjj(index))*x(j,d)<=last];
        end
    end
end
% 作业工期约束
% for j=1:Case.job_count
%     for d=1:Cases.Times
%          Constraints=[Constraints,d*x(j,d)<=Paras.T]; 
%     end
% end
% 作业不可中断
for j=1:Case.job_count
    for d=1:H-1
        dur_instant=0;
        for q=d+2:H
            dur_instant=dur_instant+x(j,q);
        end
        Constraints=[Constraints,Case.dur(j)*x(j,d)-Case.dur(j)*x(j,d+1)+ dur_instant<=Case.dur(j)];
    end
end
% 资源约束
% r代表给定资源数
for k=1:4
    for d=1:H
        total_alloc=0;
        for j=1:Case.job_count
            total_alloc=total_alloc+x(j,d)*Case.res(j,k);
        end
         Constraints=[Constraints,total_alloc<=r(k)];
    end
end
Constraints=[Constraints,0<C<H];    %目标函数C的约束
% 约束构建好了接下来是求解过程，所有的都一样
options = sdpsettings('solver','cplex','showprogress',1,'verbose',2);
% 这是最大求解时间
options.cplex.MaxTime=36000;  
options.cplex.NodeDisplayInterval=1;
sol=optimize(Constraints,Objective,options);
if strcmp(sol.info,'Successfully solved (CPLEX-IBM)')==true
    solution.Objective=value(Objective);
    display(solution.Objective);
    solution.feasibility='feasible'
    solution.info=sol.info;
    solution.x=value(C);     
    for j=1:Case.job_count   
         for t=1:H
             if(double(x(j,t))==1)
                 solution.st(j)=t;
                 break;
             end
         end
     end
             
elseif strcmp(sol.info,'Maximum iterations or time limit exceeded (CPLEX-IBM)')==true
    solution.Objective=value(Objective);
    display(solution.Objective);
    solution.feasibility='unknown';
    solution.info=sol.info;
else
    solution.feasibility='infeasible';
    solution.solvertime=sol.solvertime;
    solution.x=value(C);
end

如果这个代码看着舒服，可点击下面网页，下载代码，建议直接复制，下载后中文解释都乱码了
[1]:https://github.com/XXUZEXIN/CPLEX.git