Cplex的安装和使用实例

最近有文章复现需要用到Cplex，所以捣鼓了一下，看其他博客感觉不解渴，所以索性在这里记录一些key point

1）下载与安装：Cplex的下载可以在官方网站下载，学生是免费下载使用的，非学生可以在网上找一找，有可用的网盘下载方案；

安装过程比较简单，个人是把应用程序安装好之后再在matlab中添加的路径，参考链接：https://blog.csdn.net/robert_chen1988/article/details/44464369

使用的过程中最好把Yalmip也安装好，因为有一些必要的函数在Cplex里面没有。

做整数规划的过程中，需要让Yalmip能识别到Cplex，使用yalmiptest可以看到Cplex是否被yalmip看到，如果没有，则需要根据自身安装的Cplex版本号对yalmip目录下的solvers/definesolvers.m文件做更改，将里面的某一个版本号改为你的版本号，例如将12.8.0都改为12.9.0，对应的cplexlink1280也全改为cplexlink1290。这样保证两个套件相互都能引用。

2）功能测试：

参考链接：https://www.cnblogs.com/trux/p/7701230.html

根据上述链接里的例子运行了一下，例子稍微看一下，不用特别较真，运行没有问题即可开始写自己的程序。

上述例子运行过程基本都采用了cplex的solver来解，yalmip建模问题下面总结一下，与cplex可以做一个对比。

通常来说，yalmip中一个线性规划的写法包含三个部分：

定义变量：intvar (整型)，spdvar（实型），binvar（0/1型）等，例如：x=intvar(3,1);

定义目标函数和限制条件：为f赋值，例如：f=-[1,2,3]，g=f*x为目标函数(x1+2*x2+3*x3)，F=set(限定条件)+set(限定条件)+...+set(限定条件)；

设置Cplex为求解工具求解：sdpsettings('verbose','1','solver','cplex'); optimize(F,-f,ops);double(f);double(x);

3）解决问题：

由于Cplex是IBM搞的一个比较成熟的软件，相关的类写得已经很丰富了，建立对象后按照相应的方法格式去写就行了。所以Cplex写问题与前面描述的yalmip的写法有一些区别：Cplex的方法可以参考IBM的类描述https://www.ibm.com/support/knowledgecenter/SSSA5P_12.7.0/ilog.odms.cplex.help/refmatlabcplex/html/classCplex.html

问题：

Maximize     x1 + 2x2 + 3x3 + x4

Subject to

                   -x1 + x2 + x3 +10x4<=20

                    x1 - 3x2+ x3            <=30

                             x2         - 3x4   =  0

Bounds

                      0<=x1<=40, x2>=0, x3>=0, 2<=x4<=3, 且x4为整数

程序：

%% 初始化Cplex对象
cplex = Cplex('mipex1');
cplex.Model.sense = 'maximize';

%% 问题描述
cplex.addCols([1; 2; 3; 1], [], [0; 0; 0; 2], [40; inf; inf; 3],'CCCI');
cplex.addRows(-inf, [-1  1 1   10], 20);
cplex.addRows(-inf, [ 1 -3 1    0], 30);
cplex.addRows(   0, [ 0  1 0 -3.5],  0);

%% 问题求解
cplex.solve();

%% 写出结果
fprintf ('\nSolution status = %s \n', cplex.Solution.statusstring);
fprintf ('Solution value = %f \n', cplex.Solution.objval);
disp ('Values =');
disp (cplex.Solution.x);
disp ('Slacks =');
disp (cplex.Model.rhs - cplex.Solution.ax);

cplex下的方法可以参考上述链接，addcols方法中给出了目的函数、上下限以及变量的类型，解的结果如下：

Slacks表示的是将求得的x代入到原限制条件中，求得的值与条件中右侧值之间的差。