- 1、本文档共6页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
运筹学上机实践报告(LINGO软件)
Southwest university of science and technology
实验报告
LINGO软件在线性规划中的运用
学院名称环境与资源学院专业名称采矿工程学生姓名学号指导教师陈星明 教授
二〇一五年十一月
实验 LINGO软件在线性规划中的运用
实验目的
掌握LINGO软件求解线性规划问题的基本步骤,了解LINGO软件解决线性规划问题的基本原理,熟悉常用的线性规划计算代码,理解线性规划问题的迭代关系。
实验仪器、设备或软件
电脑,LINGO软件
实验内容
1.LINGO软件求解线性规划问题的基本原理;
2.编写并调试LINGO软件求解线性规划问题的计算代码;
实验步骤
1.使用LINGO计算并求解线性规划问题;
2.写出实验报告,并浅谈学习心得体会(线性规划的基本求解思路与方法及求解过程中出现的问题及解决方法)。
实验过程
有一艘货轮,分为前、中、后三个舱位,它们的容积与允许载重量如下表所示。现有三种商品待运,已知有关数据列于下表中。又为了航运安全,要求前、中、后舱在实际载重量上大体保持各舱最大允许载重量的比例关系。具体要求前、后舱分别与中舱之间的载重量比例偏差不超过15%,前、后舱之间不超过10%。问货轮应装载A、B、C各多少件,运费收入为最大?试建立这个问题的线性规划模型。
前 舱中 舱后 舱最大允许载重量(t)
容积(m3)2000
40003000
54001500
1500
商品数量(件)每件体积(m3/件)每件重量(t/件)运价(元/件)A
B
C600
1000
80010
5
78
6
51000
700
600
首先分析问题,建立数学模型:
确定决策变量
假设i=1,2,3分别代表商品A、B、C,8用j=1,2,3分别代表前、中、后舱,设决策变量xij为装于j舱位的第i种商品的数量(件)。
确定目标函数
商品A的件数为:
商品B的件数为:
商品A的件数为:
为使运费最高,目标函数为:
确定约束条件
前、中、后舱位载重限制为:
前、中、后舱位体积限制为:
A、B、C三种商品数量的限制条件:
各舱最大允许载重量的比例关系构成的约束条件:
且决策变量要求非负,即xij≥0,i=1,2,3;j=1,2,3。
综上所述,此问题的线性规划数学模型为:
xij≥0,i=1,2,3;j=1,2,3。
把数学模型编写成代码写入LINGO程序
编入如下代码:
!设前舱运A为x11,运B:x12,运C:x13;
!设中舱运A为x21,运B:x22,运C:x23;
!设后舱运A为x31,运B:x32,运C:x33;!单位:件;
!目标函数;
max=1000*(x11+x12+x13)+700*(x21+x22+x23)+600*(x31+x32+x33);
!数量约束;
x11+x12+x13=600;
x21+x22+x23=1000;
x31+x32+x33=800;
!容量约束;
x11*10+x21*5+x31*7=4000;
x12*10+x22*5+x32*7=5400;
x13*10+x23*5+x33*7=1500;
!重量约束;
x11*8+x21*6+x31*5=2000;
x12*8+x22*6+x32*5=3000;
x13*8+x23*6+x33*5=1500;
!平衡约束;
(2/3)*(1-0.15)=(x11*8+x21*6+x31*5)/(x12*8+x22*6+x32*5);
(1/2)*(1-0.15)=(x13*8+x23*6+x33*5)/(x12*8+x22*6+x32*5);
(4/3)*(1-0.10)=(x11*8+x21*6+x31*5)/(x13*8+x23*6+x33*5);
(x11*8+x21*6+x31*5)/(x12*8+x22*6+x32*5)=(2/3)*(1+0.15);
(x13*8+x23*6+x33*5)/(x12*8+x22*6+x32*5)=(1/2)*(1+0.15);
(x11*8+x21*6+x31*5)/(x13*8+x23*6+x33*5)=(4/3)*(1+0.10);
!整数约束;
@gin(x11);@gin(x12);@gin(x13);
@gin(x21);@gin(x22);@gin(x23);
@gin(x31);@gin(x32);@gin(x33);
如下图所示:
然后运行得出结果:
Local optimal solution found.
Objective value: 801000.0
Extended solve
文档评论(0)