系统动力学案例分析.ppt

例 经营单一商品的零售店订货问题 1、系统边界 零售店 工厂 订货 供应 顾客购货 销售 2、因果关系图和反馈回路 零售店 订货 工厂未 供订货 工厂 产量 计划 产量 工厂生 产能力 零售店 库存 零售店 销售 + + + + + + - - - R R L L A A R 3、流图 判断变量类型的其他原则 同一个反馈回路中至少有一个水准变量和一个速率变量，同时为了降低系统阶次，一般一个反馈回路中只有一个水准变量。 在因果关系图中确定变量类型时，一般先找L变量，再找R变量，最后找A变量。 一个L变量至少有一个输入R与之相连，至多有两个R变量与之相连，且这两个R变量一个为输入R，一个为输出R. 流经A变量的流一般为信息流，因此如果因果关系图中流入和流出某变量的流均为信息流，则该变量一定为A变量。 四、举例-根据下面因果关系图绘制流图 四、举例 L1 R1（利息1） C1(利率) I R1(订货量) 库存量 D Y(期望库存) （库存差额） P R1 R2 (出生人口) (人口总量) (死亡人口) C1（出生率） C2（死亡率） 组织改善 组织绩 效 组织 缺陷 。 3 DYNAMO语言 DYNAMO来源 DYNAMO，取名来自 Dynamic Models(动态模型)的混合缩写。顾名思义，DYNAMO命名的涵意在于建立真实系统的模型，藉助计算机进行系统结构、功能与动态行为的模拟。 DYNAMO和系统动力学的关系，可追溯到50年代系统动力学发展的初期。DYNAMO的前身称SIMPLE (Simulation of Industrial Management Problems with Lots of Equation) 一、DYNAMO描述动力学系统的基本原理 以库存系统为例，为简单起见，考虑输入、输出速率为常数的情况。 假定每月发货与入库各为100与80件，则库存INV每月减少20件，其动态行为是线性的，以图形表示就是随时间变化的直线。 可用数学式表达： INV现在＝INV过去十(时间间隔)x(纯速率) 若库存量在5个月前为l200件，则： INV现在＝1200件+(5月)*(80件/月-100件/月) ＝1200+5*(-20) ＝1200-100 ＝1100(件) 一、DYNAMO描述动力学系统的基本原理 在DYNAMO中，给变量带上时间下标以区别在时间上的先后。 英文字母K表示现在，J表示刚刚过去的那一时刻，L表示紧随当前的未来的那一时刻。DT表示J与K或K与L之间的时间长度。 DT DT J K L JK KL 一、DYNAMO描述动力学系统的基本原理 库存方程可用DYNAMO表示如下： 　　L INV.K = INV.J + DT*(ORRE.JK - SH.JK) 式中： INV.K –––– 库存现有量； INV.J –––– DT前的库存量； DT –––– 计算的时间间隔； ORRE –––– 在JK间隔内收到的订货量 SH –––– 在JK间隔内的发货量。 二、DYNAMO方程水准方程 水准变量：对输入和输出变量(或其中之一)进行积累的变量。 在DYNAMO中计算水准变量的方程称为水准变量方程。 在DYNAMO模型中，以L为标志写在第一列。 例如： L LEVEL.K=LEVEL.J+DT*(INFLOW.JK-OUTFLOW.JK) 式中： LEVEL –––– 水平(状态)变量； INFLOW –––– 输入速率(变化率)； OUTFLOW –––– 输出速率(变化率)； DT –––– 计算间隔(从J时刻到K时刻)。 速率变量与速率方程 在水准变量方程中代表输入与输出的变量称为速率，它由速率方程求出。 DYNAMO中，速率方程以R为标志。 例如： R BIRTHS.KL = BRF*POP.K 式中： BIRTHS —— 出生率(人/年)； BRF —— 出生率系数(1/年)； POP —— 人口(人)。 速率的值在DT时间内是不变的。进一步说，速率方程是在K时刻进行计算，而在自K至L的时间间隔(DT)中保持不变。速率的时间下标为KL。 辅助变量与方程 辅助方程以字母A为标志 例如： A DISC.K = ROOM – TEA.K R CHNG.KL= CONST*DISC.K 式中： DISC —— 茶水与室温度差(℃)； ROOM —— 室温度(℃)； TEA —— 茶水温度(℃)； CHNG —— 茶水的温度变化率(℃/分)“