Lingo求解线性规划问题Matlab与Excel数据交互.docx

1、Lingo求解线性规划问题Matlab与Excel数据交互云南财经大学实验报告 系 （院）： 专 业： 班 级： 学 号： 姓 名： 课程名称： 实验时间： 指导教师： 云南财经大学教务处制实验名称实验七、用LINGO求解线性规划问题、Matlab与Excel数据交互实验目的(1)用Lingo软件读入数据求解“六个产地八个销地”问题，学习用file命令读入*.dat、*.txt数据，ole命令读入*.xls文件的数据进行计算（可进行数据交互）。(2)编写MATLAB程序，用xlsread和xlswrite读入和写出数据，了解MATLAB与EXCEL的数据交互及相关排序计算程序的编写实验内容（算

2、法、程序、步骤和方法）实验内容（算法、程序、步骤和方法）实验内容（算法、程序、步骤和方法）实验内容（算法、程序、步骤和方法）一、 Lingo数据的读入和计算1、读入txt数据：建立1_2.txt文件：WH1 WH2 WH3 WH4 WH5 WH6 V1 V2 V3 V4 V5 V6 V7 V8 60 55 51 43 41 52 35 37 22 32 41 32 43 38 6 2 6 7 4 2 5 94 9 5 3 8 5 8 25 2 1 9 7 4 3 37 6 7 3 9 2 7 12 3 9 5 7 2 6 55 5 2 2 8 1 4 3在Lingo窗口输入：sets:ware

3、house/file(G:study数学建模LINGO10LINGO101_2.txt)/:capacity;vendors/file(G:study数学建模LINGO10LINGO101_2.txt)/:demand;links(warehouse,vendors):cost,volume;endsetsmin=sum(links:cost*volume);for(vendors(J): sum(warehouse(I):volume(I,J)=demand(J);for(warehouse(I): sum(vendors(J):volume(I,J)=capacity(I);data:ca

4、pacity=file(G:study数学建模LINGO10LINGO101_2.txt);demand=file(G:study数学建模LINGO10LINGO101_2.txt);cost=file(G:study数学建模LINGO10LINGO101_2.txt);enddataend运行2、读入*.dat数据:在Lingo窗口输入：sets:warehouse/file(G:study数学建模LINGO10LINGO101_2.dat)/:capacity;vendors/file(G:study数学建模LINGO10LINGO101_2.dat)/:demand;links(ware

5、house,vendors):cost,volume;endsetsmin=sum(links:cost*volume);for(vendors(J): sum(warehouse(I):volume(I,J)=demand(J);for(warehouse(I): sum(vendors(J):volume(I,J)=capacity(I);data:capacity=file(G:study数学建模LINGO10LINGO101_2.dat);demand=file(G:study数学建模LINGO10LINGO101_2.dat);cost=file(G:study数学建模LINGO10

6、LINGO101_2.dat);enddataend运行3、读入*.Xls数据建立1_2.xls文件：选中A1：F1，插入姓名，定义为”warehouse”选中A2：H2，插入姓名，定义为”vendors”选中A3：F3，插入姓名，定义为”capacity”选中A4：H4，插入姓名，定义为”demand”选中A2：H2，插入姓名，定义为”vendors”选中B6：I11，插入姓名，定义为”cost”在Lingo窗口输入sets:warehouse/ole(D:study数学建模2013LINGO10LINGO101_2.xls)/:capacity;vendors/ole(D:study数学建

7、模2013LINGO10LINGO101_2.xls)/:demand;links(warehouse,vendors):cost,volume;endsetsmin=sum(links:cost*volume);for(vendors(J): sum(warehouse(I):volume(I,J)=demand(J);for(warehouse(I): sum(vendors(J):volume(I,J)=capacity(I);data:capacity=ole(D:study数学建模2013LINGO10LINGO101_2.xls);demand=ole(D:study数学建模201

8、3LINGO10LINGO101_2.xls);cost=ole(D:study数学建模2013LINGO10LINGO101_2.xls);enddataend运行二、MATLAB与EXCEL数据的交互，xlsread、xlsread的使用建立Excel，输入姓名，学号，成绩在Matlab命令窗口输入：number string score=xlsread(cjb.xls); row_length=length(score(:,1);temp=score(2:row_length,:);A=(number(:,2)+number(:,3)+number(:,4)/3for i=1:5 tem

9、p1(i,:)=temp(i,:) A(i);endscore1=score(1,:) average;data_sorted=sortrows(temp1,6);data_out=score1;data_sorted;success=xlswrite(result.xls,data_out);if success=1 disp(成功!)else disp(失败!)end结论（结果）结论（结果）结论（结果）结论（结果）一、 Lingo数据的读入和计算Lingo读入*.dat,*.txt,*.xls数据计算结果都为：Global optimal solution found. Objective

10、 value: 664.0000 Total solver iterations: 15 Variable Value Reduced Cost CAPACITY( WH1) 60.00000 0.000000 CAPACITY( WH2) 55.00000 0.000000 CAPACITY( WH3) 51.00000 0.000000 CAPACITY( WH4) 43.00000 0.000000 CAPACITY( WH5) 41.00000 0.000000 CAPACITY( WH6) 52.00000 0.000000 DEMAND( V1) 35.00000 0.000000

11、 DEMAND( V2) 37.00000 0.000000 DEMAND( V3) 22.00000 0.000000 DEMAND( V4) 32.00000 0.000000 DEMAND( V5) 41.00000 0.000000 DEMAND( V6) 32.00000 0.000000 DEMAND( V7) 43.00000 0.000000 DEMAND( V8) 38.00000 0.000000 COST( WH1, V1) 6.000000 0.000000 COST( WH1, V2) 2.000000 0.000000 COST( WH1, V3) 6.000000

12、 0.000000 COST( WH1, V4) 7.000000 0.000000 COST( WH1, V5) 4.000000 0.000000 COST( WH1, V6) 2.000000 0.000000 COST( WH1, V7) 5.000000 0.000000 COST( WH1, V8) 9.000000 0.000000 COST( WH2, V1) 4.000000 0.000000 COST( WH2, V2) 9.000000 0.000000 COST( WH2, V3) 5.000000 0.000000 COST( WH2, V4) 3.000000 0.

13、000000 COST( WH2, V5) 8.000000 0.000000 COST( WH2, V6) 5.000000 0.000000 COST( WH2, V7) 8.000000 0.000000 COST( WH2, V8) 2.000000 0.000000 COST( WH3, V1) 5.000000 0.000000 COST( WH3, V2) 2.000000 0.000000 COST( WH3, V3) 1.000000 0.000000 COST( WH3, V4) 9.000000 0.000000 COST( WH3, V5) 7.000000 0.000

14、000 COST( WH3, V6) 4.000000 0.000000 COST( WH3, V7) 3.000000 0.000000 COST( WH3, V8) 3.000000 0.000000 COST( WH4, V1) 7.000000 0.000000 COST( WH4, V2) 6.000000 0.000000 COST( WH4, V3) 7.000000 0.000000 COST( WH4, V4) 3.000000 0.000000 COST( WH4, V5) 9.000000 0.000000 COST( WH4, V6) 2.000000 0.000000

15、 COST( WH4, V7) 7.000000 0.000000 COST( WH4, V8) 1.000000 0.000000 COST( WH5, V1) 2.000000 0.000000 COST( WH5, V2) 3.000000 0.000000 COST( WH5, V3) 9.000000 0.000000 COST( WH5, V4) 5.000000 0.000000 COST( WH5, V5) 7.000000 0.000000 COST( WH5, V6) 2.000000 0.000000 COST( WH5, V7) 6.000000 0.000000 CO

16、ST( WH5, V8) 5.000000 0.000000 COST( WH6, V1) 5.000000 0.000000 COST( WH6, V2) 5.000000 0.000000 COST( WH6, V3) 2.000000 0.000000 COST( WH6, V4) 2.000000 0.000000 COST( WH6, V5) 8.000000 0.000000 COST( WH6, V6) 1.000000 0.000000 COST( WH6, V7) 4.000000 0.000000 COST( WH6, V8) 3.000000 0.000000 VOLUM

17、E( WH1, V1) 0.000000 5.000000 VOLUME( WH1, V2) 19.00000 0.000000 VOLUME( WH1, V3) 0.000000 5.000000 VOLUME( WH1, V4) 0.000000 7.000000 VOLUME( WH1, V5) 41.00000 0.000000 VOLUME( WH1, V6) 0.000000 2.000000 VOLUME( WH1, V7) 0.000000 2.000000 VOLUME( WH1, V8) 0.000000 10.00000 VOLUME( WH2, V1) 1.000000

18、 0.000000 VOLUME( WH2, V2) 0.000000 4.000000 VOLUME( WH2, V3) 0.000000 1.000000 VOLUME( WH2, V4) 32.00000 0.000000 VOLUME( WH2, V5) 0.000000 1.000000 VOLUME( WH2, V6) 0.000000 2.000000 VOLUME( WH2, V7) 0.000000 2.000000 VOLUME( WH2, V8) 0.000000 0.000000 VOLUME( WH3, V1) 0.000000 4.000000 VOLUME( WH

19、3, V2) 11.00000 0.000000 VOLUME( WH3, V3) 0.000000 0.000000 VOLUME( WH3, V4) 0.000000 9.000000 VOLUME( WH3, V5) 0.000000 3.000000 VOLUME( WH3, V6) 0.000000 4.000000 VOLUME( WH3, V7) 40.00000 0.000000 VOLUME( WH3, V8) 0.000000 4.000000 VOLUME( WH4, V1) 0.000000 4.000000 VOLUME( WH4, V2) 0.000000 2.00

20、0000 VOLUME( WH4, V3) 0.000000 4.000000 VOLUME( WH4, V4) 0.000000 1.000000 VOLUME( WH4, V5) 0.000000 3.000000 VOLUME( WH4, V6) 5.000000 0.000000 VOLUME( WH4, V7) 0.000000 2.000000 VOLUME( WH4, V8) 38.00000 0.000000 VOLUME( WH5, V1) 34.00000 0.000000 VOLUME( WH5, V2) 7.000000 0.000000 VOLUME( WH5, V3

21、) 0.000000 7.000000 VOLUME( WH5, V4) 0.000000 4.000000 VOLUME( WH5, V5) 0.000000 2.000000 VOLUME( WH5, V6) 0.000000 1.000000 VOLUME( WH5, V7) 0.000000 2.000000 VOLUME( WH5, V8) 0.000000 5.000000 VOLUME( WH6, V1) 0.000000 3.000000 VOLUME( WH6, V2) 0.000000 2.000000 VOLUME( WH6, V3) 22.00000 0.000000

22、VOLUME( WH6, V4) 0.000000 1.000000 VOLUME( WH6, V5) 0.000000 3.000000 VOLUME( WH6, V6) 27.00000 0.000000 VOLUME( WH6, V7) 3.000000 0.000000 VOLUME( WH6, V8) 0.000000 3.000000 Row Slack or Surplus Dual Price 1 664.0000 -1.000000 2 0.000000 -4.000000 3 0.000000 -5.000000 4 0.000000 -4.000000 5 0.00000

23、0 -3.000000 6 0.000000 -7.000000 7 0.000000 -3.000000 8 0.000000 -6.000000 9 0.000000 -2.000000 10 0.000000 3.000000 11 22.00000 0.000000 12 0.000000 3.000000 13 0.000000 1.000000 14 0.000000 2.000000 15 0.000000 2.000000二、MATLAB与EXCEL数据交互A = 87.0000 81.0000 81.3333 87.6667 85.0000成功!备注或说明（成功或失败的原因、实验后的心得体会）当把成绩表中的姓名、学号两列对调时，则无法输出平均成绩，所以一直得不到正确结果，在定义A时，number(:,n)应对4、5、6列求平均，但在输入n值时，应该是number的2、3、4列，在程序编写时也出现了错误，所以一直没能得出结果，以后应加强对程序语句的具体细节的理解。指导教师评分（包括对实验的预习、操作和结果的综合评分）：指导教师总评：签名： 日期：