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);
% 定义变量   bin0 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

  • 19
    点赞
  • 141
    收藏
    觉得还不错? 一键收藏
  • 36
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 36
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值