系统动力学与供应链建模 课件

系统动力学

及

DYNAMO语言系统动力学

及

DYNAMO语言1什么是计算机模拟？计算机模拟是数值分析方法的一种。它用计算机程序直接建立真实系统的模型，并且通过计算机的计算了解系统随时间变化的行为或系统的特性。什么是计算机模拟？计算机模拟是数值分析方法的一种。它用计算机2为什么要进行计算机模拟？数学方法用解析式子反应变量之间的精确关系。数学模型方法例如运筹学的方法，以及概率统计的方法，对研究复杂系统问题，最优化问题以及各种决策问题都起到了巨大的作用。但是在寻求数学表达式及求解的时候，都会遇到很多问题。这些问题有的可以设法解决，有些根本解决不了。而计算机模拟采用了一种全新的思想，它充分利用计算机的优势，只凭经验数据，直接模仿客观现象，不仅利用数据关系，还利用逻辑关系描述复杂的现象。它可以利用程序把难以用数学式子表示的事件、活动、进程都模仿下来。为什么要进行计算机模拟？数学方法用解析式子反应变量之间的精确3计算机模拟语言的分类离散模拟语言如GPSS连续模拟语言如DYNAMO动态性，是指系统所包含的量具有随时间而变化的特征。比如，企业雇用员工的变动、股票市场上股票价格和交易额的波动、城市中税收和生活标准的变化、甚至糖尿病的血糖指标的变化，这些都是动态问题，可以用变量随时间变化的图形来表示。这里所说的动态性，不是随机的不稳定的动态性，而是可以预期的，有一定规律的动态性。计算机模拟语言的分类离散模拟语言如GPSS4系统动力学概述

系统动力学一种以反馈控制理论为基础，以计算机仿真技术为手段，通常用以研究复杂的社会经济系统的定量方法。适用于处理长期性和周期性的问题。如自然界的生态平衡、社会问题中的经济危机等都呈现周期性规律并需通过较长的历史阶段来观察。DYNAMO语言的名字是由dynamic（动态）和model两词的词头合并而成，表明了其预期的用途：模拟真实世界系统，使得它们随时间变化的动态行为能用计算机来跟踪、模拟。系统动力学把现实生活中的复杂系统映射成系统动力学流图，DYNAMO语言则把系统流图模型送入计算机并计算出数字结果。系统动力学概述系统动力学一种以反馈控制理论为基础，以计算机5基本思想系统动力学的基本思想是充分认识系统中的反馈和延迟，并按一定的规则从因果关系图逐步的建立系统动力学流式图的结构模式。基本思想系统动力学的基本思想是充分认识系统中的反馈和延迟，6反馈“反馈”是指信息的传送和返回。“反馈”一词的重点是在“返回”上。反馈的概念是普遍存在的。以取暖系统产生热量温暖房间为例，屋内一个和它相连的探测器将室温的信息返回给取暖系统，以此来控制系统的开关，因此也控制了屋内的温度。室温探测器是反馈装置，它和炉子、管道、抽风机一起组成了一个反馈系统。反馈“反馈”是指信息的传送和返回。“反馈”一词的重点是在“7负反馈室温高，则热风量应减小，可在室温对热风调节影响的箭头上加一个负号。反之，热风量大，则室温增加，可在热风调节对室温影响的箭头上加一个正号。从整体上看，室温影响热风量，热风量又影响了室温。从室温回到了室温，这就是一个反馈关系。另一方面，这些互相影响是相互制约的。因为温度高，则热风量减小，使室温降低。反之，室温低，则增大热风量，使室温升高。这种关系称为负反馈。图中用一个带负号的环来表示，这个环称为负反馈环，此处，负反馈环的目的是使室温接近恒定的温度。热风调节室温_负反馈室温高，则热风量应减小，可在室温对热风调节影响的箭头8正反馈相反，正反馈环总是加大环内的偏差或扰动，它具有不平衡、不断增长的特性。例如在人口系统中，人口数增加了，每年所出生的人就增加，这就使人口数按指数规律很快的增长下去。这样，从“人口数”到“每年出生的人”又返回到“人口数”之间就存在一个正反馈。增强而不是抵消环中某个元素的变化是所有正反馈环的共同特征。人口每年出生人口++正反馈相反，正反馈环总是加大环内的偏差或扰动，它具有不平衡9延迟系统动力学的另一个基本点是延迟。延迟也是普遍存在的，物理学认为“在宏观的自然界中不存在突变”。一般来说，原因并非立即就能产生结果，往往某个原因经过了一段时间才能作用产生了效果。比如在调节室温的过程中，空调机马力加大以后，室温不会马上就上升。要经过一定的热传导的过程，热量不断积累，才会使室温上升。热风调节热量积累室温室温增加延迟系统动力学的另一个基本点是延迟。延迟也是普遍存在的，物10延迟+反馈把这个因果关系图转换成流图。室温是一个存量，室温增加是对应的速率，它们用实物流相连。热风调节使热量积累，所以热量积累是一个存量，热风调节是使其增加的对应速率，它们也用实物流相连。热量积累使室温增加，这是存量对速率的影响，它是一个信息链，用温升时间常数辅助表示这个关系。热量积累延迟+反馈把这个因果关系图转换成流图。室温是一个存量，室11系统动力学流图系统动力学流图12系统动力学流图控制作用从系统获取信息由信息确定如何控制系统动力学流图控制作用从系统由信息13基本概念

（1）存量（LEVEL），图中的长方形，象一个水池，可称为水平、存量、积累量、流位，它是系统的状态，也就是系统的某个指标值。流图中用水池的水位高低来模仿系统状态值的大小。（2）速率（RATE），图中的阀门形，可称为决策函数、速度、速率、流率等。它控制着存量的变化。流图用阀门模拟控制水位的高低的机制。（3）实线称为实物流。实线连接“阀门”和“水池”，模仿控制的通路。它贯穿存量和速率。（4）虚线称为信息链，模仿信息传递的过程。它指向速率，表示根据什么信息控制速率。但不能从这条曲线看出具体控制方式的细节；基本概念（1）存量（LEVEL）14基本概念

（5）圆称为辅助变量，辅助表示流速变动的规律。（6）小圆加斜线称为常数，是系统中重要的参数。（7）不规则的闭曲线称为源或汇。源指实物的来源，汇指实物的去向，好比水的源泉和去向。它是系统之外的元素。（8）函数变量、表变量。（9）延迟。基本概念15系统动力学基本概念流图中区分实物和信息。存量是实物的积累；速率是实物随时间变化的规律；实物流是实物集散的路径；源和汇是实物的来源或去向。用信息来控制速率变动的规律，信息来自系统外部变量或内部存量；信息流是从信息源到速率的路径。系统动力学基本概念流图中区分实物和信息。存量是实物的积累；速16Part2.Dynamo语言介绍语法时间下标综述Part2.Dynamo语言介绍语法17热风调节模型L热量积累.K=热量积累.J+DT*(热风调节.JK-室温增加.JK)R热风调节.KL=温差.K*调节系数A温差.K=给定温度-室温.KC调节系数=0.1C给定温度=100L室温.K=室温.J+DT*室温增加.JKR室温增加.KL=热量积累.K/温升时间C温升时间=15N室温=80N热量积累=0SPECDT=1,LENGTH=80,PRTPER=5,PLTPER=1PRINT热风调节,热量积累,室温增加,室温PLOT热风调节=1,热量积累=2,室温增加=3,室温=4,温差=5热风调节模型L热量积累.K=热量积累.J+DT*(热风调节18语法元素语句类别标识：L、R、A、T、N、C、SPEC、PRINT、PLOT变量名时间下标：J、K、L、JK、KL等号（=）圆点(.)函数SIN、COS、NOISE、TABLE等。运算符：+、-、*、/语法元素语句类别标识：L、R、A、T、N、C、SPEC、PR19控制语句说明语句

SPECDT=0.5，LENGTH=50，PRTPER=2，PLTPER=3制表语句PRINT制图语句PLOT控制语句说明语句20函数——延迟宏函数一阶物流延迟DELAY1（IN，DEL)一阶信息延迟SMOOTH（IN，DEL）三阶物流延迟DELAY3（IN，DEL)三阶信息延迟DLNF3（IN，DEL）函数——延迟宏函数一阶物流延迟DELAY1（IN，DEL21物流延迟假设某工厂每周向河里倒一次农药残渣，两天农药就自然的被吸收净化了。在河水中农药的存量—农药吸收速度+农药流入河水的速度+由于对农药的吸收是花时间的，所以它在水中的含量不是立即减少，而是逐渐减少，因此吸收速度是农药倒入速度的延迟。物流延迟假设某工厂每周向河里倒一次农药残渣，两天农药就自然的22一阶物流延迟LSLV.K=SLV.J+DT*(IN.JK-OUT.JK)ROUT.KL=SLV.K/DELNSLV=IN*DEL一阶物流延迟LSLV.K=SLV.J+DT*(IN.JK-23一阶物流延迟如果无需了解其中动态变化的细节，前面的模型可简化为：ROUT.KL=DELAY1(IN.JK,DEL)SLVDELAY1一阶物流延迟如果无需了解其中动态变化的细节，前面的模型可简化24一阶信息延迟用过去一段时间内的销售情况作出订货的决策数，便用到了信息延迟，它所对应的数学方法是指数平均。A(N)=(4/5)*A(N-1)+(1/5)*S(N-1)前一天的销量所占的比重是4/5，前两天的销量占的比重是1/5*4/5,前三天的销量占的比重是1/5*4/5*4/5，……。A(N)=A(N-1)+(1/5)*[S(N-1)-A(N-1)]如果以天为单位，令DT＝1，得到：A(N)=A(N-1)+DT*[S(N-1)-A(N-1)]/5一阶信息延迟用过去一段时间内的销售情况作出订货的决策数，便用25一阶信息延迟LA.K=A.J+DT*ADT.JKRADT.KL=GAP.K/ATCAT=5AGAP.K=S.K-A.KNA=S（S）送入的信息S经过了四个元素的相互作用变成了输出信息A，A在时间上比S落后。因为S增大以后，它与原来A的差变大，因此控制ADT增大，使实物流加速向存量积累，存量就逐渐等于SGAPSA一阶信息延迟LA.K=A.J+DT*ADT.JK（S）26一阶信息延迟A.K=SMOOTH(S.K,AT)SMOOTH的宏LSMOOTH.K=SMOOTH.J+DT*(S.K-SMOOTH.J)/ATNSMOOTH=S（S）ATASMOOTH一阶信息延迟A.K=SMOOTH(S.K,AT)（S）AT27三阶物流延迟三阶物流延迟28三阶信息延迟（IN）三阶信息延迟（IN）29三阶延迟的响应曲线三阶物流延迟及信息延迟的响应曲线都是一样的，呈S型，如图：三阶延迟的响应曲线三阶物流延迟及信息延迟的响应曲线都是一样的30标准函数EXP(P)LONG(P)SQRT(P)SIN(P)COS(Y)标准函数EXP(P)31表函数自变量等间距变化；如果自变量没有取道表格上的数值，采用线性插值的方法去函数值。AY.K=TABLE(TY,X.K,0,8,2)TTY=7,4,3,9,4X02468Y74394表函数自变量等间距变化；X02468Y7439432逻辑函数1、跳跃函数CLIP(P,Q,R,S)2、开关函数SWITCH(P,Q,R)3、最大值函数MAX(A,B)4、最小值函数MIN(A,B)逻辑函数1、跳跃函数CLIP(P,Q,R,S)33时间控制函数1、阶跃函数STEP(P,Q)2、脉冲函数PULSE(P,Q,R)QPPQR时间控制函数1、阶跃函数STEP(P,Q)QPPQR34时间控制函数3、斜坡函数RAMP(P,Q)4、取样函数SAMPLE(P,Q,R)斜率为pQRppp时间控制函数3、斜坡函数RAMP(P,Q)斜率为pQRpp35随机变量均匀分布随机变量NOISE（）表示产生一个从-0.5到0.5之间均匀分布的随机变量。正态分布随机变量NORMRN（MEAN，STDV）随机变量均匀分布随机变量NOISE（）36生态系统的例子——简化版在一个封闭的海岛上，有植物、食草动物和食肉动物构成一个简单的生态系统。粗略的分析三者之间的关系，得到右面的因果关系图。生态系统的例子——简化版在一个封闭的海岛上，有植物、食草动物37系统动力学与供应链建模ppt课件38存量速率方程L植物.K＝植物.J＋DT*（植物出生率.JK－植物死亡率.JK）R植物出生率.KL=植物.K*K1R植物死亡率.KL=食草动物.K*M1L食草动物.K＝食草动物.J＋DT*（食草动物出生率.JK－食草动物死亡率.JK）R食草动物出生率.KL=食草动物.K*K2R食草动物死亡率.KL=食肉动物.K*M2L食肉动物.K＝食肉动物.J＋DT*（食肉动物出生率.JK－食肉动物死亡率.JK）R食肉动物出生率.KL=食肉动物.K*K3R食肉动物死亡率.KL=食肉动物.K*M3存量速率方程L植物.K＝植物.J＋DT*（植物出生率.JK39常数和初始值方程N植物＝10,000N食草动物＝1,000N食肉动物＝100CK1＝3CK2＝2CK3＝1CM1=15CM2=20CM3=0.5常数和初始值方程N植物＝10,00040生态系统的例子——复杂版植物出生率受到气候影响和其他不可预知的因素的影响，因此，植物出生率不是一个常数。在此我们引入辅助变量。生态系统的例子——复杂版植物出生率受到气候影响和其他不可预知41植物出生率的方程：A植物出生率.K=平均出生率常数+气候影响.K+随机噪音.KC平均出生率常数＝3A气候影响.K=SIN（6.28*TIME.K/10）A随机噪音.K=0.1*NOISE（）植物出生率的方程：42对食肉动物：A食肉动物寿命.K=平均寿命＋瘟疫影响.KA瘟疫影响.K=0.5*COS(6.28*TIME.K/20)C平均寿命＝2对食肉动物：43植物和食草动物的关系植物作为食草动物的食物，影响食草动物的出生率。X51015Y434食草动物占有的植物数食草动物出生率51015134植物和食草动物的关系植物作为食草动物的食物，影响食草动物的44食草动物和食肉动物的关系食肉动物出生率510150.522.5食肉动物占有的食草动物数食草动物和食肉动物的关系食肉动物出生率510150.522.45系统动力学与供应链建模ppt课件46食草动物出生率的方程：A食草动物出生率.K=Y.KAY.K=TABLE（TY,食草动物所占的植物数.K，0，15，5）TTY＝0，1，3，4A食草动物所占有的植物数.K=植物.K/食草动物.K食肉动物出生率的方程：A食肉动物出生率.K=Z.KAY.K=TABLE（TZ,食肉动物所占的食草动物数.K，0，15，5）TTZ＝0，0.5，2，2.5A食肉动物所占有的植物数.K=食草动物.K/食肉动物.K食草动物出生率的方程：47系统动力学与供应链建模ppt课件48L植物.K=植物.J+DT*(植物出生.JK-植物死亡.JK)R植物出生.KL=植物.K*植物出生率.KA植物出生率.K=平均植物出生率+气候影响.K+噪音.KA气候影响.K=SIN(6.28*TIME.K/10)A噪音.K=0.1*NOISE()R植物死亡.KL=食草动物.K*植物死亡率L食草动物.K=食草动物.J+DT*(食草动物出生.JK-食草动物死亡.JK)R食草动物出生.KL=食草动物.K*食草动物出生率.KA食草动物出生率.K=TABLE(TY,食草动物占有植物数.K,0,15,5)TTY=0,1,3,4A食草动物占有植物数.K=植物.K/食草动物.KR食草动物死亡.KL=食肉动物.K*食草动物死亡率L食肉动物.K=食肉动物.J+DT*(食肉动物出生.JK-食肉动物死亡.JK)R食肉动物出生.KL=食肉动物.K*食肉动物出生率.KA食肉动物出生率.K=TABLE(TZ,食肉动物占有食草动物.K,0,15,5)TTZ=0,0.5,2,2.5A食肉动物占有食草动物.K=食草动物.K/食肉动物.KR食肉动物死亡.KL=食肉动物.K*食肉动物死亡率.KA食肉动物死亡率.K=1/寿命.KA寿命.K=平均寿命+瘟疫影响.KA瘟疫影响.K=0.5*COS(6.28*TIME.K/20)C平均植物出生率=3C植物死亡率=0.5C食草动物死亡率=1.2C平均寿命=2N植物=10000N食草动物=1000N食肉动物=100L植物.K=植物.J+DT*(植物出生.JK-植物死亡.J49定货购物商店模拟商店经理希望了解商店中定货单、库存变化的动态，例如，当外界对商店货物需求发生了一个突然增加后，定货单如何变化，库存如何变化。定货商店业务处理流程：顾客定货后，由店员处理这些定单。按照定货单到商店的仓库取货，并将货发送给顾客。发货结束后，这笔业务就算完成了。为了保持商店的仓库总有货物，商店就必需去定货，由这些定货来补充仓库中的库存。定货购物商店模拟商店经理希望了解商店中定货单、库存变化的动态50第一步：确定基本框架存量：积存定单、库存速率：顾客定单到达速度、发货速度、商店定货速度＋－＋－＋－＋＋－第一步：确定基本框架存量：积存定单、库存＋－＋－＋－＋＋－51两条实物流：顾客定单流及商店货物流定货延迟两条实物流：顾客定单流及商店货物流定货延迟52第二步：确定速率、补充细节顾客定货速率是一个外部不可控因素，模型将其假设为一个阶跃变化是所引起的系统状态的动态变化。R顾客定货.KL=正常值＋STEP(100,5)C正常值＝1000发货速率等于积存定单/满足顾客定货时间。而满足顾客定货时间与实际库存成反比，与平均的顾客定货成正比。平均的顾客定货是顾客定货的一阶信息延迟。A平均的顾客定货.K=SMOOTH(顾客定货.JK,平均时间）C平均时间＝8满足顾客定货时间＝a/(实际库存/平均的顾客定货),比例系数a可用表函数来确定：第二步：确定速率、补充细节顾客定货速率是一个外部不可控因素，53R发货.KL=积存定单.K/满足顾客定货时间.KA满足顾客定货时间.K＝TABLE(TY,实际库存.K/平均的顾客定货.K,4,12,4)TTY=1.5,1,0.751.510.75满足顾客定货时间4812R发货.KL=积存定单.K/满足顾客定货时间.K1.51054商店定货速率是商店的定货策略，此商店经理的策略是将定货速率分成两部分。一部分是与前8周内的平均顾客定货一致，另一部分是调节实际库存到一个期望库存。调节时间是4周。期望库存是希望维持8周的平均顾客定货量。R商店定货.KL=平均顾客定货.K+(期望库存.K-实际库存.K)/调节时间A期望库存.K=平均顾客定货.K*期望库存周数C期望库存周数＝8C调节时间＝4商店定货速率是商店的定货策略，此商店经理的策略是将定货速率分55模型的初值可以用间接的方式来定。为了使模型从一个稳定的初始状态开始：N实际库存=期望库存N积存订单=顾客定货*满足顾客定货时间模型的初值可以用间接的方式来定。为了使模型从一个稳定的初始状56系统动力学与供应链建模ppt课件57L积存定单.K=积存定单.J+DT*(顾客定货.JK-发货.JK)N积存定单=满足顾客定货时间*顾客定货R发货.KL=积存定单.K/满足顾客定货时间.KA满足顾客定货时间.K=TABLE(TDFR,X.K,4,12,4)AX.K=实际库存.K/平均顾客定货.KTTDFR=1.5,1,0.75A平均顾客定货.K=SMOOTH(顾客定货.JK,平均时间)C平均时间=8A期望库存.K=库存周数*平均顾客定货.KC库存周数=8R定货延迟.KL=DELAY3(商店定货.JK,定货延迟时间)R商店定货.KL=平均顾客定货.K+(期望库存.K-实际库存.K)/调节时间C调节时间=4C定货延迟时间=6L实际库存.K=实际库存.J+DT*(定货延迟.JK-发货.JK)N实际库存=期望库存R顾客定货.KL=正常值+STEP(修正值,5)C正常值=1000C修正值=100L积存定单.K=积存定单.J+DT*(顾客定货.JK-发货58VENSIM应用实例——牛鞭效应宝洁公司（P&G）在研究“尿不湿”的市场需求时发现，该产品的零售数量相当稳定，波动性不大，但在考察分销中心的订货情况时却吃惊地发现其订单的变动程度比零售数量的波动大得多，而分销中心是将批发商的订货需求量汇总后进行订货的。通过进一步研究后发现，零售商往往根据对历史和现实销售情况的预测，确定一个较客观的订货量，但为了能应付客户需求增加的变化，他们通常会将预测订货量进行一定的放大后向批发商订货，而批发商也出于同样的考虑，会在其订货量的基础上再进行一定的放大后向分销中心订货——就这样，虽然顾客需求量并没有大的波动，但经过零售商、批发商和分销中心的订货放大后，订货量便一级一级地被放大了。供应链的信息流从末端（最终客户）向源端（原始生产商）传递时，需求信息的波动会越来越大，这种信息扭曲的放大作用在图形上很像一条甩起来的牛鞭，因此被形象地称为牛鞭效应（BullwhipEffect）。

工厂分销商批发商零售商客户VENSIM应用实例——牛鞭效应宝洁公司（P&G）在研究“尿59系统动力学与供应链建模ppt课件60流图的说明：三个存量：DistributorInventory（分销商库存）WholesalerInventory（批发商库存）RetailInventory（零售商库存）四个速率：FactoryShipments（工厂发货速度）DistributorShipments（分销商发货速度）WholesalerShipments（批发商发货速度）RetailSales（零售速度）流图的说明：三个存量：61四个常量：InventoryCoverageTime：各个企业应该持有多少周的库存ShippingDelay：表示下游企业从发出订单到接收上游企业的订货需要多长时间SaleAverageTime：给出了各个企业在进行订单预测时从多少周以前的销售情况考察起InventoryAdjustTime：表示企业调整库存所需要的时间两个开关变量：VMISwitch：是否使用VMI策略RandomOrders：零售情况是否随机变化四个常量：62三个三阶物流延迟：FactoryShipments＝DELAY3(DistributorOrders,ShippingDelay)DistributorShipments＝DELAY3(WholesalerOrders,ShippingDelay)WholesalerShipments＝DELAY3(RetailOrders,ShippingDelay)三个一阶信息延迟：DistributorSalesForecast＝SMOOTH(DistributorShipments,SaleAverageTime)WholesalerSalesForecast＝SMOOTH(WholesalerShipments,SaleAverageTime)RetailSalesForecast＝SMOOTH(RetailSales,SaleAverageTime)三个三阶物流延迟：63VMI策略VMI即供货商管理库存（VendorManagementInventory）VMI是由供应商来为客户管理库存，为他们制定库存策略和补货计划，根据客户的销售信息和库存水平为客户进行补货的一种库存管理策略和管理模式，它是供应链上成员间达成紧密义务伙伴关系后的一种结果，既是一种有效的供应链管理优化方法，也是供应链上企业联盟的一种库存管理策略。

VMI策略VMI即供货商管理库存（VendorManage64VMI策略与传统模式的比较传统模式：当一个分销商的库存水平低于安全库存量时，他需要向生产商发出采购定单，分销商自己控制了库存策略、计划、补货时间和数量等业务，由生产商为其供货。

生产商分销商基于预测的定单货物送达VMI策略与传统模式的比较传统模式：当一个分销商的库存水平低65VMI策略与传统模式的比较VMI模式：生产商从分销商处接收电子数据，这些数据代表了分销商销售和库存的真实信息（如POS和库存水平的信息等），然后生产商通过处理和分析这些信息得知分销商仓库里每一种货物的库存情况和市场需求，就可以根据它们为分销商制定和维护库存计划，在这种模式下，定单是由生产商生成的，而不是由分销商完成。

生产商分销商实时的客户需求信息按实际需求发货图7VMI模式VMI策略与传统模式的比较VMI模式：生产商从分销商处接收电66两种订单策略DistributorOrdersIFTHENELSE(VMISwitch=0,(DistributorSalesForecast+(DesiredDistributorInventory-DistributorInventory)/InventoryAdjustTime),(RetailSales+(DesiredRetailInventory*3-RetailInventory-WholesalerInventory-DistributorInventory)/InventoryAdjustTime))WholesalerOrdersIFTHENELSE(VMISwitch=0,(WholesalerSalesForecast+(DesiredWholesalerInventory-WholesalerInventory)/InventoryAdjustTime),(RetailSales+(DesiredRetailInventory*2-RetailInventory-WholesalerInventory)/InventoryAdjustTime))RetailOrders两种订单策略DistributorOrders67模拟结果（一）顾客的需求发生一次变化，传统模式订货模拟结果（一）顾客的需求发生一次变化，传统模式订货68模拟结果（二）顾客的需求随机变化，传统模式订货模拟结果（二）顾客的需求随机变化，传统模式订货69模拟结果（三）顾客的需求发生一次变化，VMI模式订货模拟结果（三）顾客的需求发生一次变化，VMI模式订货70模拟结果（二）顾客的需求随机变化，传统模式订货模拟结果（二）顾客的需求随机变化，传统模式订货71Part3.建模过程思路步骤动态性问题流图概念化Part3.建模过程思路72系统动力学建模研究步骤策略的执行对一个系统的认识问题的定义策略分析计算机模拟系统的概念化模型的建立系统动力学建模研究步骤策略的执行对一个系统的认识问题的定义73流图与系统概念化

无论是对于初学者还是经验丰富的建模者，系统概念化阶段都是系统动力学中最困难的阶段之一。概念化阶段的两大准则是：（1）明确建模目的；（2）研究问题，而不是研究整个系统。这两点是建模取得成功的先决条件。明确了研究的问题和建模目的，就好像佩带了一副滤色镜，可以滤去不必要的细节，集中精力研究反馈系统中的主要方面。它们使系统动力学的应用者和DYNAMO的建模者，不必为了研究某个事物而去考虑所有的事物。－流图与系统概念化无论是对于初学者还是经验丰富的建模者，系74结果分析系统动力学研究所得出的策略建议，不仅来自模型的最终计算结果，而且还来自模拟过程中通过各阶段迭代所得到的认识。系统动力学研究最终应该能够提出切实可行的策略建议。因为模型只是达到目的的工具，最终的目的是提高对现实世界的认识，增加对客观规律的理解。结果分析系统动力学研究所得出的策略建议，不仅来自模型的最终计75Part6.建模软件功能介绍使用说明Part6.建模软件功能介绍76图形DYNAMO图形DYNAMO建模软件是由中国人民大学信息学院经济科学实验室研发制作的利用DYNAMO语言建模的工具。它在传统DYNAMO语言的基础上，充分运用图形技术，为连续模拟语言的实际应用提供了全新的更为良好的环境和人机界面。对于各类模拟项目具有特别的意义。图形DYNAMO软件有系统动力学流图和DYNAMO语言的传统方程两种输入方式，可以根据用户的各种需要对模型进行处理，将模型、模型的编译结果、计算结果存储于数据库中，并可以按照用户的要求显示输出。图形DYNAMO图形DYNAMO建模软件是由中国人民大学信息77传统方程输入方式在传统方程模式下文本输入所要运行的传统方程软件根据DYNAMO语言的规则进行编译，形成动方程和静方程编译成功后由用户进行变量的挑选若挑选完毕则可进行计算最后根据用户的要求显示计算结果，可以图形显示，也可以文本显示传统方程输入方式在传统方程模式下文本输入所要运行的传统方程78系统动力学流图输入方式在流图模式下在画布上画出流图检查流图并形成动方程和静方程由用户进行对动方程和静方程的补充由用户进行变量挑选补充成功且已挑选变量，则可以进行计算最后根据用户的要求显示计算结果，可以图形显示，也可以文本显示系统动力学流图输入方式79创建模型创建模型80编辑方程编辑方程81模型编译模型编译82结果计算结果计算83绘制流图绘制流图84补充方程补充方程85变量选择变量选择86Dynamo软件使用综述传统方程、流图两种输入方式；模型建立方程编辑/流图绘制模型编译静方程补充/动方程补充变量选择/结果显示。Dynamo软件使用综述传统方程、流图两种输入方式；87VENSIMVENSIM是一个建模工具,可以建立动态系统的概念化的,文档化的,仿真,分析和优化模型.PLE(个人学习版)和PLEplus是VENSIM的缩减版,主要用来简单化学习动态系统.提供了一种简单富有弹性的方法从常规的循环或储存过程和流程图建立模型.VENSIMVENSIM是一个建模工具,可以建立动态系统的概88系统动力学

及

DYNAMO语言系统动力学

及

DYNAMO语言89什么是计算机模拟？计算机模拟是数值分析方法的一种。它用计算机程序直接建立真实系统的模型，并且通过计算机的计算了解系统随时间变化的行为或系统的特性。什么是计算机模拟？计算机模拟是数值分析方法的一种。它用计算机90为什么要进行计算机模拟？数学方法用解析式子反应变量之间的精确关系。数学模型方法例如运筹学的方法，以及概率统计的方法，对研究复杂系统问题，最优化问题以及各种决策问题都起到了巨大的作用。但是在寻求数学表达式及求解的时候，都会遇到很多问题。这些问题有的可以设法解决，有些根本解决不了。而计算机模拟采用了一种全新的思想，它充分利用计算机的优势，只凭经验数据，直接模仿客观现象，不仅利用数据关系，还利用逻辑关系描述复杂的现象。它可以利用程序把难以用数学式子表示的事件、活动、进程都模仿下来。为什么要进行计算机模拟？数学方法用解析式子反应变量之间的精确91计算机模拟语言的分类离散模拟语言如GPSS连续模拟语言如DYNAMO动态性，是指系统所包含的量具有随时间而变化的特征。比如，企业雇用员工的变动、股票市场上股票价格和交易额的波动、城市中税收和生活标准的变化、甚至糖尿病的血糖指标的变化，这些都是动态问题，可以用变量随时间变化的图形来表示。这里所说的动态性，不是随机的不稳定的动态性，而是可以预期的，有一定规律的动态性。计算机模拟语言的分类离散模拟语言如GPSS92系统动力学概述

系统动力学一种以反馈控制理论为基础，以计算机仿真技术为手段，通常用以研究复杂的社会经济系统的定量方法。适用于处理长期性和周期性的问题。如自然界的生态平衡、社会问题中的经济危机等都呈现周期性规律并需通过较长的历史阶段来观察。DYNAMO语言的名字是由dynamic（动态）和model两词的词头合并而成，表明了其预期的用途：模拟真实世界系统，使得它们随时间变化的动态行为能用计算机来跟踪、模拟。系统动力学把现实生活中的复杂系统映射成系统动力学流图，DYNAMO语言则把系统流图模型送入计算机并计算出数字结果。系统动力学概述系统动力学一种以反馈控制理论为基础，以计算机93基本思想系统动力学的基本思想是充分认识系统中的反馈和延迟，并按一定的规则从因果关系图逐步的建立系统动力学流式图的结构模式。基本思想系统动力学的基本思想是充分认识系统中的反馈和延迟，94反馈“反馈”是指信息的传送和返回。“反馈”一词的重点是在“返回”上。反馈的概念是普遍存在的。以取暖系统产生热量温暖房间为例，屋内一个和它相连的探测器将室温的信息返回给取暖系统，以此来控制系统的开关，因此也控制了屋内的温度。室温探测器是反馈装置，它和炉子、管道、抽风机一起组成了一个反馈系统。反馈“反馈”是指信息的传送和返回。“反馈”一词的重点是在“95负反馈室温高，则热风量应减小，可在室温对热风调节影响的箭头上加一个负号。反之，热风量大，则室温增加，可在热风调节对室温影响的箭头上加一个正号。从整体上看，室温影响热风量，热风量又影响了室温。从室温回到了室温，这就是一个反馈关系。另一方面，这些互相影响是相互制约的。因为温度高，则热风量减小，使室温降低。反之，室温低，则增大热风量，使室温升高。这种关系称为负反馈。图中用一个带负号的环来表示，这个环称为负反馈环，此处，负反馈环的目的是使室温接近恒定的温度。热风调节室温_负反馈室温高，则热风量应减小，可在室温对热风调节影响的箭头96正反馈相反，正反馈环总是加大环内的偏差或扰动，它具有不平衡、不断增长的特性。例如在人口系统中，人口数增加了，每年所出生的人就增加，这就使人口数按指数规律很快的增长下去。这样，从“人口数”到“每年出生的人”又返回到“人口数”之间就存在一个正反馈。增强而不是抵消环中某个元素的变化是所有正反馈环的共同特征。人口每年出生人口++正反馈相反，正反馈环总是加大环内的偏差或扰动，它具有不平衡97延迟系统动力学的另一个基本点是延迟。延迟也是普遍存在的，物理学认为“在宏观的自然界中不存在突变”。一般来说，原因并非立即就能产生结果，往往某个原因经过了一段时间才能作用产生了效果。比如在调节室温的过程中，空调机马力加大以后，室温不会马上就上升。要经过一定的热传导的过程，热量不断积累，才会使室温上升。热风调节热量积累室温室温增加延迟系统动力学的另一个基本点是延迟。延迟也是普遍存在的，物98延迟+反馈把这个因果关系图转换成流图。室温是一个存量，室温增加是对应的速率，它们用实物流相连。热风调节使热量积累，所以热量积累是一个存量，热风调节是使其增加的对应速率，它们也用实物流相连。热量积累使室温增加，这是存量对速率的影响，它是一个信息链，用温升时间常数辅助表示这个关系。热量积累延迟+反馈把这个因果关系图转换成流图。室温是一个存量，室99系统动力学流图系统动力学流图100系统动力学流图控制作用从系统获取信息由信息确定如何控制系统动力学流图控制作用从系统由信息101基本概念

（1）存量（LEVEL），图中的长方形，象一个水池，可称为水平、存量、积累量、流位，它是系统的状态，也就是系统的某个指标值。流图中用水池的水位高低来模仿系统状态值的大小。（2）速率（RATE），图中的阀门形，可称为决策函数、速度、速率、流率等。它控制着存量的变化。流图用阀门模拟控制水位的高低的机制。（3）实线称为实物流。实线连接“阀门”和“水池”，模仿控制的通路。它贯穿存量和速率。（4）虚线称为信息链，模仿信息传递的过程。它指向速率，表示根据什么信息控制速率。但不能从这条曲线看出具体控制方式的细节；基本概念（1）存量（LEVEL）102基本概念

（5）圆称为辅助变量，辅助表示流速变动的规律。（6）小圆加斜线称为常数，是系统中重要的参数。（7）不规则的闭曲线称为源或汇。源指实物的来源，汇指实物的去向，好比水的源泉和去向。它是系统之外的元素。（8）函数变量、表变量。（9）延迟。基本概念103系统动力学基本概念流图中区分实物和信息。存量是实物的积累；速率是实物随时间变化的规律；实物流是实物集散的路径；源和汇是实物的来源或去向。用信息来控制速率变动的规律，信息来自系统外部变量或内部存量；信息流是从信息源到速率的路径。系统动力学基本概念流图中区分实物和信息。存量是实物的积累；速104Part2.Dynamo语言介绍语法时间下标综述Part2.Dynamo语言介绍语法105热风调节模型L热量积累.K=热量积累.J+DT*(热风调节.JK-室温增加.JK)R热风调节.KL=温差.K*调节系数A温差.K=给定温度-室温.KC调节系数=0.1C给定温度=100L室温.K=室温.J+DT*室温增加.JKR室温增加.KL=热量积累.K/温升时间C温升时间=15N室温=80N热量积累=0SPECDT=1,LENGTH=80,PRTPER=5,PLTPER=1PRINT热风调节,热量积累,室温增加,室温PLOT热风调节=1,热量积累=2,室温增加=3,室温=4,温差=5热风调节模型L热量积累.K=热量积累.J+DT*(热风调节106语法元素语句类别标识：L、R、A、T、N、C、SPEC、PRINT、PLOT变量名时间下标：J、K、L、JK、KL等号（=）圆点(.)函数SIN、COS、NOISE、TABLE等。运算符：+、-、*、/语法元素语句类别标识：L、R、A、T、N、C、SPEC、PR107控制语句说明语句

SPECDT=0.5，LENGTH=50，PRTPER=2，PLTPER=3制表语句PRINT制图语句PLOT控制语句说明语句108函数——延迟宏函数一阶物流延迟DELAY1（IN，DEL)一阶信息延迟SMOOTH（IN，DEL）三阶物流延迟DELAY3（IN，DEL)三阶信息延迟DLNF3（IN，DEL）函数——延迟宏函数一阶物流延迟DELAY1（IN，DEL109物流延迟假设某工厂每周向河里倒一次农药残渣，两天农药就自然的被吸收净化了。在河水中农药的存量—农药吸收速度+农药流入河水的速度+由于对农药的吸收是花时间的，所以它在水中的含量不是立即减少，而是逐渐减少，因此吸收速度是农药倒入速度的延迟。物流延迟假设某工厂每周向河里倒一次农药残渣，两天农药就自然的110一阶物流延迟LSLV.K=SLV.J+DT*(IN.JK-OUT.JK)ROUT.KL=SLV.K/DELNSLV=IN*DEL一阶物流延迟LSLV.K=SLV.J+DT*(IN.JK-111一阶物流延迟如果无需了解其中动态变化的细节，前面的模型可简化为：ROUT.KL=DELAY1(IN.JK,DEL)SLVDELAY1一阶物流延迟如果无需了解其中动态变化的细节，前面的模型可简化112一阶信息延迟用过去一段时间内的销售情况作出订货的决策数，便用到了信息延迟，它所对应的数学方法是指数平均。A(N)=(4/5)*A(N-1)+(1/5)*S(N-1)前一天的销量所占的比重是4/5，前两天的销量占的比重是1/5*4/5,前三天的销量占的比重是1/5*4/5*4/5，……。A(N)=A(N-1)+(1/5)*[S(N-1)-A(N-1)]如果以天为单位，令DT＝1，得到：A(N)=A(N-1)+DT*[S(N-1)-A(N-1)]/5一阶信息延迟用过去一段时间内的销售情况作出订货的决策数，便用113一阶信息延迟LA.K=A.J+DT*ADT.JKRADT.KL=GAP.K/ATCAT=5AGAP.K=S.K-A.KNA=S（S）送入的信息S经过了四个元素的相互作用变成了输出信息A，A在时间上比S落后。因为S增大以后，它与原来A的差变大，因此控制ADT增大，使实物流加速向存量积累，存量就逐渐等于SGAPSA一阶信息延迟LA.K=A.J+DT*ADT.JK（S）114一阶信息延迟A.K=SMOOTH(S.K,AT)SMOOTH的宏LSMOOTH.K=SMOOTH.J+DT*(S.K-SMOOTH.J)/ATNSMOOTH=S（S）ATASMOOTH一阶信息延迟A.K=SMOOTH(S.K,AT)（S）AT115三阶物流延迟三阶物流延迟116三阶信息延迟（IN）三阶信息延迟（IN）117三阶延迟的响应曲线三阶物流延迟及信息延迟的响应曲线都是一样的，呈S型，如图：三阶延迟的响应曲线三阶物流延迟及信息延迟的响应曲线都是一样的118标准函数EXP(P)LONG(P)SQRT(P)SIN(P)COS(Y)标准函数EXP(P)119表函数自变量等间距变化；如果自变量没有取道表格上的数值，采用线性插值的方法去函数值。AY.K=TABLE(TY,X.K,0,8,2)TTY=7,4,3,9,4X02468Y74394表函数自变量等间距变化；X02468Y74394120逻辑函数1、跳跃函数CLIP(P,Q,R,S)2、开关函数SWITCH(P,Q,R)3、最大值函数MAX(A,B)4、最小值函数MIN(A,B)逻辑函数1、跳跃函数CLIP(P,Q,R,S)121时间控制函数1、阶跃函数STEP(P,Q)2、脉冲函数PULSE(P,Q,R)QPPQR时间控制函数1、阶跃函数STEP(P,Q)QPPQR122时间控制函数3、斜坡函数RAMP(P,Q)4、取样函数SAMPLE(P,Q,R)斜率为pQRppp时间控制函数3、斜坡函数RAMP(P,Q)斜率为pQRpp123随机变量均匀分布随机变量NOISE（）表示产生一个从-0.5到0.5之间均匀分布的随机变量。正态分布随机变量NORMRN（MEAN，STDV）随机变量均匀分布随机变量NOISE（）124生态系统的例子——简化版在一个封闭的海岛上，有植物、食草动物和食肉动物构成一个简单的生态系统。粗略的分析三者之间的关系，得到右面的因果关系图。生态系统的例子——简化版在一个封闭的海岛上，有植物、食草动物125系统动力学与供应链建模ppt课件126存量速率方程L植物.K＝植物.J＋DT*（植物出生率.JK－植物死亡率.JK）R植物出生率.KL=植物.K*K1R植物死亡率.KL=食草动物.K*M1L食草动物.K＝食草动物.J＋DT*（食草动物出生率.JK－食草动物死亡率.JK）R食草动物出生率.KL=食草动物.K*K2R食草动物死亡率.KL=食肉动物.K*M2L食肉动物.K＝食肉动物.J＋DT*（食肉动物出生率.JK－食肉动物死亡率.JK）R食肉动物出生率.KL=食肉动物.K*K3R食肉动物死亡率.KL=食肉动物.K*M3存量速率方程L植物.K＝植物.J＋DT*（植物出生率.JK127常数和初始值方程N植物＝10,000N食草动物＝1,000N食肉动物＝100CK1＝3CK2＝2CK3＝1CM1=15CM2=20CM3=0.5常数和初始值方程N植物＝10,000128生态系统的例子——复杂版植物出生率受到气候影响和其他不可预知的因素的影响，因此，植物出生率不是一个常数。在此我们引入辅助变量。生态系统的例子——复杂版植物出生率受到气候影响和其他不可预知129植物出生率的方程：A植物出生率.K=平均出生率常数+气候影响.K+随机噪音.KC平均出生率常数＝3A气候影响.K=SIN（6.28*TIME.K/10）A随机噪音.K=0.1*NOISE（）植物出生率的方程：130对食肉动物：A食肉动物寿命.K=平均寿命＋瘟疫影响.KA瘟疫影响.K=0.5*COS(6.28*TIME.K/20)C平均寿命＝2对食肉动物：131植物和食草动物的关系植物作为食草动物的食物，影响食草动物的出生率。X51015Y434食草动物占有的植物数食草动物出生率51015134植物和食草动物的关系植物作为食草动物的食物，影响食草动物的132食草动物和食肉动物的关系食肉动物出生率510150.522.5食肉动物占有的食草动物数食草动物和食肉动物的关系食肉动物出生率510150.522.133系统动力学与供应链建模ppt课件134食草动物出生率的方程：A食草动物出生率.K=Y.KAY.K=TABLE（TY,食草动物所占的植物数.K，0，15，5）TTY＝0，1，3，4A食草动物所占有的植物数.K=植物.K/食草动物.K食肉动物出生率的方程：A食肉动物出生率.K=Z.KAY.K=TABLE（TZ,食肉动物所占的食草动物数.K，0，15，5）TTZ＝0，0.5，2，2.5A食肉动物所占有的植物数.K=食草动物.K/食肉动物.K食草动物出生率的方程：135系统动力学与供应链建模ppt课件136L植物.K=植物.J+DT*(植物出生.JK-植物死亡.JK)R植物出生.KL=植物.K*植物出生率.KA植物出生率.K=平均植物出生率+气候影响.K+噪音.KA气候影响.K=SIN(6.28*TIME.K/10)A噪音.K=0.1*NOISE()R植物死亡.KL=食草动物.K*植物死亡率L食草动物.K=食草动物.J+DT*(食草动物出生.JK-食草动物死亡.JK)R食草动物出生.KL=食草动物.K*食草动物出生率.KA食草动物出生率.K=TABLE(TY,食草动物占有植物数.K,0,15,5)TTY=0,1,3,4A食草动物占有植物数.K=植物.K/食草动物.KR食草动物死亡.KL=食肉动物.K*食草动物死亡率L食肉动物.K=食肉动物.J+DT*(食肉动物出生.JK-食肉动物死亡.JK)R食肉动物出生.KL=食肉动物.K*食肉动物出生率.KA食肉动物出生率.K=TABLE(TZ,食肉动物占有食草动物.K,0,15,5)TTZ=0,0.5,2,2.5A食肉动物占有食草动物.K=食草动物.K/食肉动物.KR食肉动物死亡.KL=食肉动物.K*食肉动物死亡率.KA食肉动物死亡率.K=1/寿命.KA寿命.K=平均寿命+瘟疫影响.KA瘟疫影响.K=0.5*COS(6.28*TIME.K/20)C平均植物出生率=3C植物死亡率=0.5C食草动物死亡率=1.2C平均寿命=2N植物=10000N食草动物=1000N食肉动物=100L植物.K=植物.J+DT*(植物出生.JK-植物死亡.J137定货购物商店模拟商店经理希望了解商店中定货单、库存变化的动态，例如，当外界对商店货物需求发生了一个突然增加后，定货单如何变化，库存如何变化。定货商店业务处理流程：顾客定货后，由店员处理这些定单。按照定货单到商店的仓库取货，并将货发送给顾客。发货结束后，这笔业务就算完成了。为了保持商店的仓库总有货物，商店就必需去定货，由这些定货来补充仓库中的库存。定货购物商店模拟商店经理希望了解商店中定货单、库存变化的动态138第一步：确定基本框架存量：积存定单、库存速率：顾客定单到达速度、发货速度、商店定货速度＋－＋－＋－＋＋－第一步：确定基本框架存量：积存定单、库存＋－＋－＋－＋＋－139两条实物流：顾客定单流及商店货物流定货延迟两条实物流：顾客定单流及商店货物流定货延迟140第二步：确定速率、补充细节顾客定货速率是一个外部不可控因素，模型将其假设为一个阶跃变化是所引起的系统状态的动态变化。R顾客定货.KL=正常值＋STEP(100,5)C正常值＝1000发货速率等于积存定单/满足顾客定货时间。而满足顾客定货时间与实际库存成反比，与平均的顾客定货成正比。平均的顾客定货是顾客定货的一阶信息延迟。A平均的顾客定货.K=SMOOTH(顾客定货.JK,平均时间）C平均时间＝8满足顾客定货时间＝a/(实际库存/平均的顾客定货),比例系数a可用表函数来确定：第二步：确定速率、补充细节顾客定货速率是一个外部不可控因素，141R发货.KL=积存定单.K/满足顾客定货时间.KA满足顾客定货时间.K＝TABLE(TY,实际库存.K/平均的顾客定货.K,4,12,4)TTY=1.5,1,0.751.510.75满足顾客定货时间4812R发货.KL=积存定单.K/满足顾客定货时间.K1.510142商店定货速率是商店的定货策略，此商店经理的策略是将定货速率分成两部分。一部分是与前8周内的平均顾客定货一致，另一部分是调节实际库存到一个期望库存。调节时间是4周。期望库存是希望维持8周的平均顾客定货量。R商店定货.KL=平均顾客定货.K+(期望库存.K-实际库存.K)/调节时间A期望库存.K=平均顾客定货.K*期望库存周数C期望库存周数＝8C调节时间＝4商店定货速率是商店的定货策略，此商店经理的策略是将定货速率分143模型的初值可以用间接的方式来定。为了使模型从一个稳定的初始状态开始：N实际库存=期望库存N积存订单=顾客定货*满足顾客定货时间模型的初值可以用间接的方式来定。为了使模型从一个稳定的初始状144系统动力学与供应链建模ppt课件145L积存定单.K=积存定单.J+DT*(顾客定货.JK-发货.JK)N积存定单=满足顾客定货时间*顾客定货R发货.KL=积存定单.K/满足顾客定货时间.KA满足顾客定货时间.K=TABLE(TDFR,X.K,4,12,4)AX.K=实际库存.K/平均顾客定货.KTTDFR=1.5,1,0.75A平均顾客定货.K=SMOOTH(顾客定货.JK,平均时间)C平均时间=8A期望库存.K=库存周数*平均顾客定货.KC库存周数=8R定货延迟.KL=DELAY3(商店定货.JK,定货延迟时间)R商店定货.KL=平均顾客定货.K+(期望库存.K-实际库存.K)/调节时间C调节时间=4C定货延迟时间=6L实际库存.K=实际库存.J+DT*(定货延迟.JK-发货.JK)N实际库存=期望库存R顾客定货.KL=正常值+STEP(修正值,5)C正常值=1000C修正值=100L积存定单.K=积存定单.J+DT*(顾客定货.JK-发货146VENSIM应用实例——牛鞭效应宝洁公司（P&G）在研究“尿不湿”的市场需求时发现，该产品的零售数量相当稳定，波动性不大，但在考察分销中心的订货情况时却吃惊地发现其订单的变动程度比零售数量的波动大得多，而分销中心是将批发商的订货需求量汇总后进行订货的。通过进一步研究后发现，零售商往往根据对历史和现实销售情况的预测，确定一个较客观的订货量，但为了能应付客户需求增加的变化，他们通常会将预测订货量进行一定的放大后向批发商订货，而批发商也出于同样的考虑，会在其订货量的基础上再进行一定的放大后向分销中心订货——就这样，虽然顾客需求量并没有大的波动，但经过零售商、批发商和分销中心的订货放大后，订货量便一级一级地被放大了。供应链的信息流从末端（最终客户）向源端（原始生产商）传递时，需求信息的波动会越来越大，这种信息扭曲的放大作用在图形上很像一条甩起来的牛鞭，因此被形象地称为牛鞭效应（BullwhipEffect）。

工厂分销商批发商零售商客户VENSIM应用实例——牛鞭效应宝洁公司（P&G）在研究“尿147系统动力学与供应链建模ppt课件148流图的说明：三个存量：DistributorInventory（分销商库存）WholesalerInventory（批发商库存）RetailInventory（零售商库存）四个速率：FactoryShipments（工厂发货速度）DistributorShipments（分销商发货速度）WholesalerShipments（批发商发货速度）RetailSales（零售速度）流图的说明：三个存量：149四个常量：InventoryCoverageTime：各个企业应该持有多少周的库存ShippingDelay：表示下游企业从发出订单到接收上游企业的订货需要多长时间SaleAverageTime：给出了各个企业在进行订单预测时从多少周以前的销售情况考察起InventoryAdjustTime：表示企业调整库存所需要的时间两个开关变量：VMISwitch：是否使用VMI策略RandomOrders：零售情况是否随机变化四个常量：150三个三阶物流延迟：FactoryShipments＝DELAY3(DistributorOrders,ShippingDelay)DistributorShipments＝DELAY3(WholesalerOrders,ShippingDelay)WholesalerShipments＝DELAY3(RetailOrders,ShippingDelay)三个一阶信息延迟：DistributorSalesForecast＝SMOOTH(DistributorShipments,SaleAverageTime)WholesalerSalesForecast＝SMOOTH(WholesalerShipments,SaleAverageTime)RetailSalesForecast＝SMOOTH(RetailSales,SaleAverageTime)三个三阶物流延迟：151VMI策略VMI即供货商管理库存（VendorManagementInventory）VMI是由供应商来为客户管理库存，为他们制定库存策略和补货计划，根据客户的销售信息和库存水平为客户进行补货的一种库存管理策略和管理模式，它是供应链上成员间达成紧密义务伙伴关系后的一种结果，既是一种有效的供应链管理优化方法，也是供应链上企业联盟的一种库存管理策略。

VMI策略VMI即供货商管理库存（VendorManage152VMI策略与传统模式的比较传统模式：当一个分销商的库存水平低于安全库存量时，他需要向生产商发出采购定单，分