MATLAB调用LINGO程序并交换数据

思路

第一次写文章完全不熟练…在做比赛时想要用列生成解决混合整数规划问题，一开始在python上写，约束的种类非常多，debug时非常痛苦，稍微修改下约束工作量都很大。做了大概一周后还是想用lingo。lingo书写约束和修改约束非常舒适，而且能顺便测试有没有写错的地方，很节省时间。
主要思路是：
1、用lingo完成模型建立和求解，并将数据输出到excel表格中，方便运行中读取和后续模型计算使用；
2、用matlab调用lingo脚本，借助excel获得lingo的计算结果，处理计算结果后写入excel，用于后续迭代。
需要解决的两个关键问题：
1、按要求调用lingo脚本；
2、计算结果顺利传递。

一、lingo脚本调用

这一部分是为了在MATLAB中通过命令的形式运行lingo代码。

主要方法：
1、将lingo模型(.lg4文件)变成lingo命令脚本文件(.ltf文件)格式。
可以在lingo窗口的File-New中找到。变量、约束不变，文件结尾加一行 GO 即可（不加封号）。开头处加一行 SET TERSEO 1 可以避免在后面matlab窗口中输出冗长的结果。新加这类语句都不需要加封号。还有很多其他命令可以参考谢金星老师的11580405《优化建模与LINDOLINGO软件》这本书。

2、.bat批处理文件
根据lingo安装的位置和命令脚本文件的位置写批处理文件。用Notepad可以在新建文件处找到。具体写法参考Matlab调用Lingo脚本文件的方法中的方法二。

3、获得执行命令
例如上述文章中dos(‘C:\Users\10096\Desktop\test.bat’)，可直接添加在matlab代码中运行

二、lingo接收数据与输出数据

lingo中自带文件输出函数。这里我选择了EXCEL作为数据交换的媒介，但matlab操作EXCEL比较慢，导致效率不高。lingo也可以向txt文件输出数据，也许可以提高速度，有试过的兄弟麻烦留个言分享下。

lingo自带函数@ole用于和EXCEL交换数据，这一部分有很多资料。主要的方法就是在EXCEL表格中先定义和变量数量对应的输出位置的名称，选中输出位置的表格后，右键定义名称。之后再@ole函数中写好文件位置和定义的名称就行。

我在调试时发现@ole函数的文件路径写错了，但由于文件再lingo程序运行时被打开，并没有影响lingo书写和读取数据。lingo输出结果时需要EXCEL表格在内存中，而且不会自动保存。这一点在后面写matlab代码时要特别注意。

三、matlab代码

这里遇到最大的问题是lingo读取数据后不自动关闭保存文件，这导致lingo程序运行后，EXCEL文件被lingo锁定，matlab自带的xlswrite和xlsread都不能正常使用。我这里的解决办法就是先在matlab中打开excel表格，运行lingo程序后，在matlab中保存数据并关闭表格。这样操作后xlswrite和xlsread就不会出问题了。

有关在matlab中操作excel的函数参考matlab和Excel的交互 非xlsread和xlswrite（1-Excel基础）。

主要代码结构：

try
    Excel=actxGetRunningServer('Excel.Application');%如果Excel 服务器已经打开，返回其句柄
catch
    Excel=actxserver('Excel.Application');%如果Excel服务器没有打开，则创建一个Excel服务器，并返回句柄
end

Workbook = Excel.Workbooks.Open('C:\Users\LENVOLE1\Desktop\matlab_lingo\lingo_excel.xlsx');%打开excel表格
dos('C:\Users\LENVOLE1\Desktop\matlab_lingo\test.bat');%运行lingo程序
invoke(Workbook,'save');%保存lingo运行结果
invoke(Workbook,'close');%关闭excel表格
X = xlsread('C:\Users\LENVOLE1\Desktop\matlab_lingo\lingo_excel.xlsx', 'XXout', 'A2:EG2');%此时读取数据正常
xlswrite('C:\Users\LENVOLE1\Desktop\matlab_lingo\lingo_excel.xlsx',AA,'A','B2'); 

Excel.Quit; % 关闭 Excel
Excel.delete; % 删除对象

目前这个解决办法看起来挺憨憨的…有兄弟有好方法麻烦留言嗷！

最后

现在看这个问题也不是很麻烦，但是当时自己试着做时真的很搞心态，希望能节省大家一点时间。祝一切顺利！