本文将结合官方文档和其他相关介绍,对Matlab + Gurobi的使用做一个完全零基础的入门介绍(我也是小白,这也是对自己学习的一个记录)
一、Gurobi优化器快速入门指导 - Windows
1.1从求解一个简单模型开始 - 了解Gurobi命令行
该问题是一个非常简单的铸币厂生产硬币问题,一共有五种类型的硬币,价值分别对应0.01$, 0.05$, 0.1$, 0.25$, 和1$,每种硬币由不同的原材料组成,一共四种,使用的原材料分量如表1-1所示(原问题详见官网:Solving a Simple Model - The Gurobi Command Line):
Penny | Nickel | Dime | Quarter | Dollar | |
Copper(Cu) | 0.06g | 3.8g | 2.1g | 5.2g | 7.2g |
Nickel(Ni) | 1.2g | 0.2g | 0.5g | 0.2g | |
Zinc(Zi) | 2.4g | 0.5g | |||
Manganses(Mn) | 0.3g |
表1-1
不同材料的剩余量分别为Cu=1000, Ni=50, Zi=50, Mn=50,决策变量不同硬币的生产数量(注意这里的变量需为正数),求如何生产能让产出的总价值最大。
在具体的LP格式的Gurobi文档中,需要注意以下四点:
1)不同部分的顺序。需要按照目标函数,约束条件,变量范围,整数要求等顺序进行;
2)分离符号。不同的符号如+-号,数字,变量之间必须用空格或新的一行间隔开;
3)变量安排。约束条件中所有的变量必须放在约束的左手边;
4)变量默认范围。如果没有特殊说明,变量的下界为0,上界无穷大。
在cmd中运行如下代码,可以得到铸币厂的最大产值为113.45$。
gurobi_cl D:\Software\Gurobiv9\win64\examples\data\coins
用如下命令可以得到包含决策变量最优解的文件:
gurobi_cl ResultFile=coins.sol D:\Software\Gurobiv9\win64\examples\data\coins.lp
若没有指定目录,便会用户文件夹下生成结果文件。
1.2gurobi交互层(桌面快捷方式)
1.2.1读取和优化模型
m=read('D:\Software\Gurobiv9\win64\examples\data\coins')
可以直接读取一个写好的模型,并将其赋予变量m中。
m.read('D:\Software\Gurobiv9\win64\examples\*\coin*')
也可以使用通配符*,但注意末尾的*必须存在。如果在电脑上安装了gzip,bzip2,7zip等,可以在文件名后面加上合适的后缀来读写压缩文件。
m.optimize()
对调用优化方法对模型m进行优化。
1.2.2检查结果
当一个模型被求解后,可以用以下命令查看模型的最优解
m.printAttr(&#