供应链管理实验报告姓名：学号：班级：***师：***相关问题说明：一、实验性质和教学目的本实验是供应链管理课内安排的上机操作实验。

目的是根据供应链中供应管理和需求管理的实际问题，抽象出相应的数学模型，利用Lingo 优化软件求解模型，通过对求解结果的分析，一方面使学生更好地理解和掌握供应链管理的有关原理和概念，另一方面锻炼学生利用计算机等现代工具分析求解实际问题的动手能力，以达到学以致用的最终目的。

二、实验基本要求要求学生：1. 实验前认真做好理论准备，仔细阅读实验指导书；2. 遵从教师指导，认真完成实验任务，按时按质提交实验报告。

三、主要参考资料1．LINGO软件2. 优化建模与LINDO/LINGO软件，清华大学出版社，20053．运筹学编写组主编，运筹学（第三版），清华大学出版社，19904．《供应链管理：战略、规划与运作》（第3版）（清华管理学系列英文版教材），（美）乔普拉（Chopra，S）,(美)迈因德尔（Meindl，P.）著，清华大学出版社5. 供应链管理（第3版）（工商管理经典译丛），乔普拉等著，陈荣秋等译，中国人民大学出版社实验内容1．Lavare 公司是芝加哥郊区主要的不锈钢水槽制造厂，公司现在正在制定来年需求和供给管理计划。

预计每月分销商的需求如表2所示。

Lavare 公司的产能由工厂雇佣的操作工人数量决定，工人每月工作20天，每天8小时，其他时间的工作算加班，正常工作时间每小时工资15美元，加班费每小时22美元。

每个工人每月的加班时间不得超过20小时。

工厂现雇佣工人数为250名，每个不锈钢水槽的生产需要2小时，单位库存持有成本为每月3美元，单件产品生产成本为40美元。

每单位的销售价格为125美元销售给分销商。

假定没有转包生产。

假定Lavare 公司最初有4000个单位库存，并希望维持年底也有同样多的库存。

市场调查显示，降价1%能够给当月增加20%的销售量，并使得未来两个月销售量的10%提前至当月。

例如，在三月份进行一个百分点的降价促销活动可以使3月的销售量增加3000件（=0.2*15000），并且使得4月的销售量提前1800件（=0.1*18000），5月份的销售量提前2500件（=0.1*25000）至三月。

(a) 假定没有促销活动，这一年的最优生产计划是什么？年利润是多少？这项计划的成本是多少（b ）在4月还是7月进行促销活动更好？各能增加多少利润？（c ）如果水槽的售价从125美元涨至250美元，最佳的促销时间是否需要调整？为什么？解：（a ）决策变量：H t = 第t 月新雇佣的人数(t = 1,..,12) L t =第 t 月解聘的人数(t = 1,..,12) W t =第t 月的工人的人数 (t = 1,..,12) O t = 第t 月加班的小时数 (t = 1,..,12)I t = 第t 月末持有的库存量，单位千 (t = 1,..,12) C t =第t 月外包的生产量，单位千(t = 1,..,12) P t =第t 月的生产量，单位千(t = 1,..,12) 已知参数D t =第t 月的需求量(t = 1,..,12) 数学模型：Minimize 121212121111240022340tttt t i i i W O IP ====+++∑∑∑∑Subject to: 库存约束：1,1,..,12t t t t I P I D t -+-==0124000I I ==加班时间限制200,1,...,12t t O W t -≤=生产能力限制21600,1,...12t t t P O W t --≤=劳动力限制250,0,...,12t W t ==（a ）：最优生产计划：成本计算源代码： s e t s :j i e d u a n /1 2 3 4 5 6 7 8 9 10 11 12/; P e r i o d s (j i e d u a n ): D , W , O , I , P ; e n d s e t s !目标函数;m i n=@s u m(p e r i o d s(t)|t#G E#1:2400*W(t)+22*O(t)+3*I(t)+40*P(t));!需求约束;@f o r(p e r i o d s(t)|t#G T#1:I(t-1)+P(t)=D(t)+I(t));I0+P(1)=D(1)+I(1);I0=4000;I(12)=4000;@f o r(p e r i o d s(t)|t#G T#1:W(t)=W(t-1));W(1)=W0;W0=250;@f o r(p e r i o d s(t)|t#G E#1:P(t)<=80*W(t)+O(t)/2);@f o r(p e r i o d s(t)|t#G E#1:O(t)<=20*W(t));!这里是数据;d a t a:D=10000,11000,15000,18000,25000,26000,30000,29000,21000,18000,14000,11000;e n d d a t ae n dLINGO计算运行后得到成本为16820000美元利润计算源代码：s e t s:j i e d u a n/123456789101112/;P e r i o d s(j i e d u a n):D,W,O,I,P;e n d s e t s!目标函数;m i n=@s u m(p e r i o d s(t)|t#G E#1:125*D(t)-2400*W(t)-22*O(t)-3*I(t)-40*P(t));!需求约束;@f o r(p e r i o d s(t)|t#G T#1:I(t-1)+P(t)=D(t)+I(t));I0+P(1)=D(1)+I(1);I0=4000;I(12)=4000;@f o r(p e r i o d s(t)|t#G T#1:W(t)=W(t-1));W(1)=W0;W0=250;@f o r(p e r i o d s(t)|t#G E#1:P(t)<=80*W(t)+O(t)/2);@f o r(p e r i o d s(t)|t#G E#1:O(t)<=20*W(t));!这里是数据;d a t a:D=10000,11000,15000,18000,25000,26000,30000,29000,21000,18000,14000,11000;e n d d a t ae n dL I N G O运行结果得利润=11680000美元四月份促销：成本计算代码：s e t s:j i e d u a n/123456789101112/;P e r i o d s(j i e d u a n):D,W,O,I,P;e n d s e t s!目标函数;m i n=@s u m(p e r i o d s(t)|t#G E#1:2400*W(t)+22*O(t)+3*I(t)+40*P(t));!需求约束;@f o r(p e r i o d s(t)|t#G T#1:I(t-1)+P(t)=D(t)+I(t));I0+P(1)=D(1)+I(1);I0=4000;I(12)=4000;@f o r(p e r i o d s(t)|t#G T#1:W(t)=W(t-1));W(1)=W0;W0=250;@f o r(p e r i o d s(t)|t#G E#1:P(t)<=80*W(t)+O(t)/2);@f o r(p e r i o d s(t)|t#G E#1:O(t)<=20*W(t));!这里是数据;d a t a:D=10000,11000,15000,26700,22500,23400,30000,29000,21000,18000,14000,11000;e n d d a t ae n d利润计算代码：s e t s:j i e d u a n/123456789101112/;P e r i o d s(j i e d u a n):D,W,O,I,P;e n d s e t s!目标函数;max=@sum(periods(t)| t #GE# 1:125*D(t) )-1.25*D(4)-@sum(periods(t)| t #GE# 1: 2400*W(t) +22*O(t) + 3*I(t) + 40*P(t));!需求约束;@f o r(p e r i o d s(t)|t#G T#1:I(t-1)+P(t)=D(t)+I(t));I0+P(1)=D(1)+I(1);I0=4000;I(12)=4000;@f o r(p e r i o d s(t)|t#G T#1:W(t)=W(t-1));W(1)=W0;W0=250;@f o r(p e r i o d s(t)|t#G E#1:P(t)<=80*W(t)+O(t)/2);@f o r(p e r i o d s(t)|t#G E#1:O(t)<=20*W(t));!这里是数据;D=10000,11000,15000,26700,22500,23400,30000,29000,21000,18000,14000,11000;e n d d a t ae n d11857220美七月份促销：成本计算代码：s e t s:j i e d u a n/123456789101112/;P e r i o d s(j i e d u a n):D,W,O,I,P;e n d s e t s!目标函数;m i n=@s u m(p e r i o d s(t)|t#G E#1:2400*W(t)+22*O(t)+3*I(t)+40*P(t));!需求约束;@f o r(p e r i o d s(t)|t#G T#1:I(t-1)+P(t)=D(t)+I(t));I0+P(1)=D(1)+I(1);I0=4000;I(12)=4000;@f o r(p e r i o d s(t)|t#G T#1:W(t)=W(t-1));W(1)=W0;W0=250;@f o r(p e r i o d s(t)|t#G E#1:P(t)<=80*W(t)+O(t)/2);@f o r(p e r i o d s(t)|t#G E#1:O(t)<=20*W(t));!这里是数据;D=10000,11000,15000,18000,25000,26000,41000,26100,18900,18000,14000,11000;e n d d a t ae n d利润计算代码：s e t s:j i e d u a n/123456789101112/;P e r i o d s(j i e d u a n):D,W,O,I,P;e n d s e t s!目标函数;max=@sum(periods(t)| t #GE# 1:125*D(t)-2.5*D(7))-@sum(periods(t)| t #GE# 1: 2400*W(t) +22*O(t) + 3*I(t) + 40*P(t));!需求约束;@f o r(p e r i o d s(t)|t#G T#1:I(t-1)+P(t)=D(t)+I(t));I0+P(1)=D(1)+I(1);I0=4000;I(12)=4000;@f o r(p e r i o d s(t)|t#G T#1:W(t)=W(t-1));W(1)=W0;W0=250;@f o r(p e r i o d s(t)|t#G E#1:P(t)<=80*W(t)+O(t)/2);@f o r(p e r i o d s(t)|t#G E#1:O(t)<=20*W(t));!这里是数据;d a t a:D=10000,11000,15000,18000,25000,26000,41000,26100,18900,18000,14000,11000;e n d d a t ae n dLINGO运行结果显示7月份成=17367300 利润=11831450综上结果，四月利润增加177200美元，七月增加151450美元。