• SIMULINK是MATLAB提供的实现动态系统建模和仿真的一 个软件包，它适用于连续系统和离散系统，也适用于线性系 统和非线性系统。它采用系统模块直观的描述系统典型环节， 因此可十分方便的建立系统模型而不需花较多时间编程。此 软件有两个明显的功能：仿真和连接，亦即可以利用鼠标器 在模型窗口上画出所需的控制系统模型，然后利用该软件提 供的功能来对系统直接进行仿真。很明显，这种做法使得一 个很复杂系统的输入变得相当容易。因此，本文选用 SIMULINK作为仿真语言。
• 2. (R, S)策略 该策略和（Q, R）策略一样，都是连续
性检查类型的策略，也就是要随时检查库存 状态，当发现库存降低到订货点水平R时， 开始订货，订货后使最大库存保持不变，即 为常量 S，若发出订单时库存量为I，则其订 货量即为（S-I）。该策略和（Q, R）策略的 不同之处在于其订货量是按实际库存而定， 因而订货量是可变的。
基于啤酒游戏的供应链牛鞭效 应分析与研究
案例：啤酒游戏
• 一、供应链牛鞭效应研究的目的及意义
随着市场全球化和竞争的加剧，企业之间 的竞争己变成供应链之间的竞争。所以，如何 提高整条供应链的增值能力，增强供应链的竞 争力，成为各节点企业共同的目标，但是“牛 鞭效应”的存在削弱了供应链的增值能力和竞 争力。
Dt
应链模型。假设以下的事件序列:在每个周期t， 零售商①接收下游订单,②满足下游需求,③收到 上游来货,④察看库存水平,⑤向上游发出订单。 在零售商发出订单和接收到货物之间有一个固定 的时间间隔，因此在一个周期t末发出的订单将在 周期t+L的开始接收到货物。特别是，提前期L包 含一段时间的订单传送和TP时间周期的商品生产和 发送。这里假设订单传送为1个周期，TP=3倍的时 间周期。
L t
t
DD D D D z z Fz
t t1
t
t1
1
1 1
t1
1 11
1
指数平滑算法的传递函数如下：
F
z
1
z 9 8 1
S 当 Ot =0.9时， t
已知 O' =
St
D
11
z1TP2
-库存量，库存量记为I，净
库存量记为 。k＝0 O z z z zz '
O
1
TP
D
1 1
1
O
T 1
订货量（订货单位）
Step Response
时间（订货周期）（sec） 图4-5 订货量曲线
图4-5表明，在某一周顾客需求突增后， 零售商的订货量将会产生较长时间（即大 于10周）的波动，直至40周后才变得平稳。
为了看出供应链牛鞭效应从下游到上游传递 的现象、验证牛鞭效应的存在 ，我们将前面获取 的模型串联，并在起始端加入一组在[-1，l]区间上 的随机信号，由仿真结果曲线可知初始需求的波动 被逐级放大。
②连续性检查的固定订货点、最大库存策略， 即（R, S）策略；
③周期性检查策略，即（t, S）策略；
④综合库存策略，即（t, R, S）策略。
• 1. (Q, R)策略 该策略的基本思想是：对库存进行连续
性检查，当库存降低到订货点水平R时，即 发出一个订货，每次的订货量保持不变，都 为固定值Q。该策略适用于需求量大、缺货 费用较高、需求波动性很大的情形。
• 讨论：有哪些库存订货策略？
订货点法
• 订货点法又称订购点法. 订货点法指的是：对于某
种物料或产品，由于生产或销售的原因而逐渐减
少，当库存量降低到某一预先设定的点时，即开
始发出订货单(采购单或加工单) 来补充库存，直
至库存量降低到安全库存时，发出的订单所定购 的物料(产品)刚好到达仓库，补充前一时期的消耗， 此一订货的数值点，即称为订货点。
图4-6 供应链牛鞭效应的传递过程
初始 需求 曲线
供应链 第一层 的响应
订货量（订货单位）
需求量（需求单位）
时间（订货周期）（sec）
图4-7
时间（订货周期）（sec）
图4-8
供应链 第二层 的响应
供应链 第三层 的响应
订货量（订货单位）
订货量（订货单位）
时间（订货周期）（sec）
时间（订货周期）（sec） 图4-9
• 二、国内外研究现状
目前对供应链中牛鞭效应的研究国外学者起步 较早，大多数的研究侧重于证实牛鞭效应的存在， 分析可能产生的原因，提出一些降低和缓解牛鞭 效应影响的方法。
到目前为止，牛鞭效应的的研究工作主要集 中于三个方面:第一，牛鞭效应的存在性;第二， 牛鞭效应的定量化研究;第三，对牛鞭效应的抑制 及其稳定性优化。但目前的研究尚处于游离阶段， 绝大多数的文献往往是对上述三个方面的其中一 个或两个方面进行研究。
• 1.2决策规则
• 在任何的最大库存订货策略下，订货决策如
下:
Ot=St-库存量 • Ot是周期t末所做出的订货决策， • St是周期t内的最大库存， • 库存量等于净库存量加上马上到来的货(或者
叫WIP)，净库存量等于仓库存货减去未交付 订货的积压。
L L
S t Dt k t
L
•
这 计里 值，Dt
该策略不设订货点，只设固定检查周期和最大
库存量。该策略适用于一些不很重要的、或使用量
不大的物资。
• 4. (t, R, S)策略 该策略是策略（t,S）和策略（R,S）的综合。
这种补给策略有一个固定的检查周期t、最大库存量 S、固定订货点水平R。当经 过一定的检查周期t后，
若库存低于订货点，则发出订货，否则，不订货。 订货量的大小等于最大库存量减去检查时的库存量。
s 0.1178
转换到S域为:
L
t
• 2传递函数模型PID的Simulink控制仿真 下图是PID控制器的控制结构框图：
比例
r(t) + e(t)
-
+ u(t)
积分
被控对象
c(t)
+ +
微分
图 5-1 PID控制结构图
加入到传递函数模型中去。得到PID控制 下的供应链模型，如下所示：
• 讨论：如何调节PID的系数？
1
P
1
D
O z z z zz I
' O
1 O
T
1
P
1 1
1
O 1
1 D
1
其中
No Image
延迟因子来确保事件序列能够与实际情
况相符。 1
1 z
1
表示延迟
1t
个采样周期，
对应的时间函数为
S T z zz O
D 2
I
t
P
1 1 1
O 1 D 1 1
• 库存量等于净库存量加上马上到来 的货(或者叫WIP)
• 3. (t, S)策略
该策略是每隔一定时期检查一次库存，并发出 一次订货，把现有库存补充到最大库存水平S，如 果检查时库存量为I，则订货量为S-I。经过固定的 检查期t，发出订货，这时，库存量为I1，订货量为 (S-I1)。经过一定的时间（LT），库存补充(S-I1)， 库存到达A点。再经过一个固定的检查时期 t，又 发出一次订货，订货量为(S-I2)，经过一定的时间 （LT-订货提前期，可以为随机变量），库存有达 到新的高度B。如此周期性检查库存，不断补给。
O
z 1
1
D T P 2 D
1
T z
1
z T z T O
T z T T z D
1 1
1
a
2
P 1
Hale Waihona Puke 1 1 1 1a
2
P 1
1
a
z 1 T
P 2 T
a z 1
z
z T
P T
a 3 T
P T
2
a
T a z 1 z
T a z 1 z
• 当经过固定的检查 时期到达A点时，此时库存已降 低到订货点水平线R之下，因而应发出一次订货，订 货量等于最大库存量S与当时的库存量I1的差（SI1）。经过一定的订货 提前期后在B点订货到达， 库存补充到C点，在第二个检查期到来时，此时库存 位置在D，比订货点水平位置线高，无须订货。第三 个检查期到来时，库存点在E， 等于订货点，又发 出一次订货，订货量为(S-I3),如此，周期进行下去， 实现周期性库存补给。
• 比例（P）调节作用：是按比例反应系统的偏差， 系统一旦出现了偏差，比例调节立即产生调节作 用用以减少 偏差。比例作用大，可以加快调节，
减少误差，但是过大的比例，使系统的稳定性下 降，甚至造成系统的不稳定。
• 积分（I）调节作用：是使系统消除稳态误差，提 高无差度。因为有误差，积分调节就进行，直至 无差，积分调节停止，积分调节输出一常值。积 分作用的强弱取决与积分时间常数Ti，Ti越小， 积分作用就越强。反之Ti大则积分作用弱，加入 积分调节可使系统稳定性下降，动态响应变慢。 积分作用常与另两种调节规律结合，组成PI调节 器或PID调节器。
• 三、本案例研究内容
本案例主要是对供应链牛鞭效应进行分析 研究，应用Matlab/Simulink仿真软件对啤酒分 销模型的三级供应链中牛鞭效应现象进行了仿真 设计，分析了不同平滑指数对牛鞭效应的影响， 基于控制理论提出减弱牛鞭效应的方法和措施， 重点阐述了PID控制，进而验证策略的可行性。
• 四、供应链牛鞭效应建模
• 首先，我们必须查出放大的来源。这种从下游到上 游的库存放大可能来源于订货策略，预测方法，提 前期等方面，然而，这些都与库存直接相关，因此 我们从简单的库存模型开始研究。
• 1传递函数模型
• 传递函数模型是零初始条件下线性系统输出量的拉 普拉斯变换与输入量的拉普拉斯变换之比。