Lingo求解线性规划案例1——生产计划问题
凯鲁嘎吉 - 博客园
http://www.cnblogs.com/kailugaji/
说明：
Lingo版本：
某工厂明年根据合同，每个季度末向销售公司提供产品，有关信息如下表。若当季生产的产品过多，季末有积余，则一个季度每积压一吨产品需支付存贮费O.2万元。现该厂考虑明年的最佳生产方案，使该厂在完成合同的情况下，全年的生产费用最低。试建立模型。
|
季度j |
生产能力aj(吨) |
生产成本dj (万元／吨) |
需求量bj(吨) |
|
1
|
30
|
15．O
|
20
|
|
2
|
40
|
14.O
|
20
|
|
3
|
20
|
15．3
|
30
|
|
4
|
10
|
14．8
|
10 |
解：现在我们对本问题定义三种不同形式的决策变量，从而从不同的途径来构建模型。
(1)设工厂第j季度生产产品xj吨。
首先，考虑约束条件：
第一季度末工厂需交货20吨；故应有x1≥20；
第一季度末交货后积余(x1-20)吨；
第二季度末工厂需交货20吨，故应有x1-20+ x2≥20；类似地，应有x1+ x2-40+ x3≥30；
第四季度末供货后工厂不能积压产品，故应有x1+ x2+ x3-70+ x4=10；
又考虑到工厂每个季度的生产能力，故应有0≤xj≤aj。
其次，考虑目标函数：
第一季度工厂的生产费用为15.0x1，
第二季度工厂的费用包括生产费用14 x2及积压产品的存贮费0.2(x1-20)；
类似地，第三季度费用为15.3x3+0.2(x1+ x2-40)，
第四季度费用为14.8x4+0.2(x1+ x2+ x3-70)。
工厂一年的费用即为这四个季度费用之和。
整理后，得下列线性规划模型：
Min f =15.6 x1+14.4 x2+15.5 x3+14.8 x4-26
s.t．x1+ x2 ≥40
x1+ x2+ x3 ≥70
x1+ x2+ x3+ x4=80
20≤x1≤30
0≤x2≤40
0≤x3≤20
0≤x4≤10
Lingo程序：
min=15.6*x1+14.4*x2+15.5*x3+14.8*x4-26;
x1>=20;
x1<=30;
x1+x2>=40;
x2<=40;
x1+x2+x3>=70;
x3<=20;
x1+x2+x3+x4=80;
x4<=10;
end
结果为：
Global optimal solution found.
Objective value: 1165.000
Infeasibilities: 0.000000
Total solver iterations: 0 Variable Value Reduced Cost
X1 20.00000 0.000000
X2 40.00000 0.000000
X3 10.00000 0.000000
X4 10.00000 0.000000 Row Slack or Surplus Dual Price
1 1165.000 -1.000000
2 0.000000 -0.1000000
3 10.00000 0.000000
4 20.00000 0.000000
5 0.000000 1.100000
6 0.000000 0.000000
7 10.00000 0.000000
8 0.000000 -15.50000
9 0.000000 0.7000000
(2)设第j季度工厂生产的产品为xj吨，第j季度初存贮的产品为yj吨(显然，y1＝0)。
因为每季度初的存贮量为上季度存贮量、生产量之和与上季度的需求量之差，又考虑到第四季度末存贮量为零，故有；
x1-20=y2，
y2+x2-20=y3，
y3+x3-30=y4，
y4+x4=10；
同时，每季度的生产量不能超过生产能力：xj≤aj；而工厂四个季度的总费用由每季的生产费用与存贮费用组成，于是得线性规划：
minf=15.Ox1+O.2y2+14x2+O.2y3+15.3x3+O.2y4+14.8x4
s．t． x1-y2=20
y2+x2-y3＝20
y3+x3-y4=30
y4+x4=10
0≤x1≤30 0≤x2≤40
0≤x3≤20 0≤x4≤10
0≤yj j=2,3,4
Lingo程序：
min=15*x1+14*x2+15.3*x3+14.8*x4+0.2*y2+0.2*y3+0.2*y4;
x1-y2=20;
x2+y2-y3=20;
y3+x3-y4=30;
y4+x4=10;
x1<30;
x2<40;
x3<20;
x4<10;
end
结果为：
Global optimal solution found.
Objective value: 1165.000
Infeasibilities: 0.000000
Total solver iterations: 0 Variable Value Reduced Cost
X1 20.00000 0.000000
X2 40.00000 0.000000
X3 10.00000 0.000000
X4 10.00000 0.000000
Y2 0.000000 0.1000000
Y3 20.00000 0.000000
Y4 0.000000 0.7000000 Row Slack or Surplus Dual Price
1 1165.000 -1.000000
2 0.000000 -15.00000
3 0.000000 -15.10000
4 0.000000 -15.30000
5 0.000000 -14.80000
6 10.00000 0.000000
7 0.000000 1.100000
8 10.00000 0.000000
9 0.000000 0.000000
(3)设第i季度生产而用于第j季度末交货的产品数量为xij吨。
根据合同要求，必须有：
x11=20，
x12+x22=20，
x13+x23+x33=30，
x14+x24+x34+x44=10。
又每季度生产而用于当季和以后各季交货的产品数不可能超过该季度工厂的生产能力，故应有。
X11+x12+x13+x14≤30，
x22+x23+x24≤40，
x33+x34≤20，
x44≤10。
第i季度生产的用于第j季度交货的每吨产品的费用cij=dj+0.2(j-i)，于是，有线性规划模型。
minf=15.0x11+15.2x12+15.4x13+15.6x14+14x22+14.2x23+14.4x24+15.3x33+15.5x34+14.8x44
s.t. x11=20
x12+x22=20
x13+x23+x33=30
x14+x24+x34+x44=10
x11+x12+x13+x14≤30
x22+x23+x24≤40
x33+x34≤20
x44≤10
xij≥0， i=1，…，4;j=1，…，4，j≥i。
Lingo程序为：
min=15*x11+15.2*x12+15.4*x13+15.6*x14+14*x22+14.2*x23+14.4*x24+15.3*x33+15.5*x34+14.8*x44;
x11=20;
x12+x22=20;
x13+x23+x33=30;
x14+x24+x34+x44=10;
x11+x12+x13+x14<30;
x22+x23+x24<40;
x33+x34<20;
x44<10;
end
结果为：
Global optimal solution found.
Objective value: 1165.000
Infeasibilities: 0.000000
Total solver iterations: 4 Variable Value Reduced Cost
X11 20.00000 0.000000
X12 0.000000 0.1000000
X13 0.000000 0.1000000
X14 0.000000 0.1000000
X22 20.00000 0.000000
X23 20.00000 0.000000
X24 0.000000 0.000000
X33 10.00000 0.000000
X34 0.000000 0.000000
X44 10.00000 0.000000 Row Slack or Surplus Dual Price
1 1165.000 -1.000000
2 0.000000 -15.00000
3 0.000000 -15.10000
4 0.000000 -15.30000
5 0.000000 -15.50000
6 10.00000 0.000000
7 0.000000 1.100000
8 10.00000 0.000000
9 0.000000 0.7000000
Lingo求解线性规划案例1——生产计划问题的更多相关文章
- Lingo求解线性规划案例4——下料问题
凯鲁嘎吉 - 博客园 http://www.cnblogs.com/kailugaji/ 造纸厂接到定单,所需卷纸的宽度和长度如表 卷纸的宽度 长度 5 7 9 10000 30000 20000 工 ...
- Lingo求解线性规划案例3——混料问题
凯鲁嘎吉 - 博客园 http://www.cnblogs.com/kailugaji/ 某糖果厂用原料A.B和C按不向比率混合加工而成甲.乙.丙三种糖果(假设混合加工中不损耗原料).原料A.B.C ...
- Lingo求解线性规划案例2——多阶段投资问题
凯鲁嘎吉 - 博客园 http://www.cnblogs.com/kailugaji/ 某公司现有资金30万元可用于投资,5年内有下列方案可供采纳: 1号方案:在年初投资1元,2年后可收回1. ...
- 用Lingo求解线性规划问题
第一步:输入目标条件和约束条件.每行以分号隔开.然后点击工具栏上的Solve按钮,或Lingo菜单下的Solve子菜单. 第二步:检查report中的结果. 默认情况下,Lingo不进行灵敏度分析. ...
- 易普优APS应用案例：线束行业生产计划排产
一.线束行业生产现状 (1)产品种类以及标准繁多,生产计划难协调 线束行业的生产,虽然原材料不多,但线束产品却多达几万种.一般线束企业,虽然不是每个月都生产数万种产品,但每月生产的产品品种在300种以 ...
- Python求解线性规划——PuLP使用教程
简洁是智慧的灵魂,冗长是肤浅的藻饰.--莎士比亚<哈姆雷特> 1 PuLP 库的安装 如果您使用的是 Anaconda[1] 的话(事实上我也更推荐这样做),需要先激活你想要安装的虚拟环境 ...
- 什么是APS高级计划排程（生产计划排产）系统主要功能模块有哪些？
什么是APS高级计划排程(生产计划排产)系统? APS高级计划排程(高级计划排产)系统主要解决“在有限产能条件下,交期产能精确预测.工序生产与物料供应最优详细计划”的问题.APS高级计划排程(高级计划 ...
- APS中生产计划排程模块的基本原理
高级计划系统(APS)作为ERP和MES的补充,用于协调物流.开发瓶颈资源和保证交货日期. APS包括需求和供应计划.运输和生产计划排程等各种供应链计划模块,本文主要介绍APS中生产计划排程模块的基本 ...
- 为什么众多软件厂商无法提供APS高级计划排程系统？工厂目前生产计划是怎么排产的？
一.行业现状如想了解一下目前现状,去考察一下上了ERP的企业,会发现一个有趣的现象该企业无论ERP软件搞得如何如火如荼,似乎都与生产调度人员无关. 车间里或者生产线上的生产作业计划.生产过程的调度和管 ...
随机推荐
- UVC 驱动调用过程与驱动框架的简单分析
内核:Linux-3.4.2 驱动:drivers\media\video\uvc\uvc_driver.c UVC 驱动整体调用流程: /* 打开设备描述符 */ 1. open: uvc_v4l2 ...
- iOS面试准备之思维导图
以思维导图的方式对iOS常见的面试题知识点进行梳理复习,文章xmind点这下载,文章图片太大查看不了也点这下载 你可以在公众号 五分钟学算法 获取数据结构与算法相关的内容,准备算法面试 公众号回复 g ...
- DLL加载，设置相对路径
DLL加载,设置相对路径 1. 加载dll方法之一:(./ 代表当前目录,../ 代表上层目录)包含头文件的相对路径(当前路径为源代码路径,路径 “../../” 当前项目文件夹上级目录),链接lib ...
- vb.net 使用ip查詢(Host Name)(WorkGroup Name)(MAC Address)-運用cmd及nbtstat命令
Sub nbtstat(ByVal ip As String) Dim strRst, strRst1, strRst2, strRst3 As String Dim n1, n2, n3 As In ...
- 软件测试工程师这样面试，拿到offer的几率是80%
面试难还是不难?取决于面试者的底蕴(气场+技能).心态和认知及沟通技巧.面试其实可以理解为一场聊天和谈判,在这过程中有心理.思想上的碰撞和博弈.其实你只需要搞清楚一个逻辑:“面试官为什么会这样问?他希 ...
- Netty实战三之Netty的组件和设计
有关Netty,我们可以从两个视角来讨论Netty:类库的视角以及框架的视角,对于使用Netty编写高效的.可重用的和可维护的代码来说,两者缺一不可. Netty解决了两个响应的关注领域,可以大致标志 ...
- 浅谈spring中AOP以及spring中AOP的注解方式
AOP(Aspect Oriented Programming):AOP的专业术语是"面向切面编程" 什么是面向切面编程,我的理解就是:在不修改源代码的情况下增强功能.好了,下面在 ...
- Yii2基本概念之——事件(Event)
说起事件(event),我们可是一点都不陌生.现实生活当中的事件无处不在,比如你发了一条微博,触发了一条事件,导致关注你的人收到了一条消息,看到你发的内容:比如你通过支付宝买东西,付了款,触发一个事件 ...
- hihoCoder编程练习赛67
题目1 : 序列 时间限制:20000ms 单点时限:1000ms 内存限制:256MB 描述 给定两个正整数 n, P,求满足以下两个条件的长度为 n 的序列 ai 个数: 1. 1 ≤ ai ≤ ...
- img图像标签和超链接标签a
图像标签语法:<img src="" alt="".../> img属性:src="" 显示图像的URLalt="& ...