赞
踩
首先从YLMIP获取YALMIP
,然后从ShareApps获取CPLEX
。
注意:CPLEX
与matlab
存在兼容性问题,因此需要在IBM产品兼容性检测平台确定matlab
版本对应的CPLEX
版本。
YALMIP
压缩包下载之后,解压至可用路径,然后在matlab
中把解压文件所在路径进行路径添加。YALMIP
安装完毕,可以通过yalmiptest
命令进行检测。
CPLEX
CPLEX
为软件,进行正常的软件安装操作即可。把安装目录下的\CPLEX_Studio\cplex\matlab
文件夹路径进行路径添加。CPLEX
安装完毕,同样可以通过yalmiptest
命令进行检测。
上述操作完成后,键入yalmiptest
,我们可以发现CPLEX
对应版本处于found
状态。
这里主要参考知乎专栏文章Matlab+yalmip+cplex安装教程及算例。
m i n f ( x ) = 2 x 1 + 3 x 2 s . t . x 1 + x 2 ≥ 350 x 1 ≥ 100 2 x 1 + x 2 ≤ 600 x 1 , x 2 ≥ 0 min \quad f(x)=2x_1+3x_2 \\ s.t. \quad x_1+x_2 \geq 350\\ x_1 \geq 100 \\ 2x_1 + x_2 \leq 600 \\ x_1,x_2 \geq 0 minf(x)=2x1+3x2s.t.x1+x2≥350x1≥1002x1+x2≤600x1,x2≥0
%定义变量 x=sdpvar(2,1); %目标函数 obj=2*x(1)+3*x(2); %约束条件 constraint=[]; constraint=[constraint,x(1)+x(2)>350]; constraint=[constraint,x(1)>100]; constraint=[constraint,2*x(1)+x(2)<600]; constraint=[constraint,x(2)>0]; %求解 ops = sdpsettings('solver','cplex','verbose',1); ops.cplex.display='on'; ops.cplex.timelimit=600; ops.cplex.mip.tolerances.mipgap=0.001; % 诊断求解可行性 disp('开始求解') diagnostics=optimize(constraint,obj,ops); if diagnostics.problem==0 disp('Solver thinks it is feasible') elseif diagnostics.problem == 1 disp('Solver thinks it is infeasible') pause(); else disp('Timeout, Display the current optimal solution') end
sdpvar:实数变量
intvar:整数变量
binvar:0-1变量
value:变量或表达式的值
ops = sdpsettings('solver','cplex','verbose',1):设置求解方法为调用 Cplex
CPXPARAM_TimeLimit 600
CPXPARAM_MIP_Tolerances_MIPGap 0.001
Tried aggregator 1 time.
LP Presolve eliminated 2 rows and 0 columns.
Reduced LP has 2 rows, 2 columns, and 4 nonzeros.
Presolve time = 0.17 sec. (0.00 ticks)
Iteration log . . .
Iteration: 1 Dual objective = 700.000000
Solver thinks it is feasible
Copyright © 2003-2013 www.wpsshop.cn 版权所有,并保留所有权利。