第五章系统仿真及系统动力学系统仿真概述系统动力学结构模型化原理基本反馈回路的DYNAMO仿真分析一.系统仿真概述一.系统仿真概述一.系统仿真概述一.系统仿真概述（二）系统仿真方法一.系统仿真概述（二）系统仿真方法一.系统仿真概述（二）系统仿真方法一.系统仿真概述（二）系统仿真方法一.系统仿真概述（二）系统仿真方法一.系统仿真概述（二）系统仿真方法一.系统仿真概述（二）系统仿真方法一.系统仿真概述（二）系统仿真方法二.系统动力学结构模型化原理2、研究对象及其结构特点3、工作程序（二）SD结构模型化原理（二）SD结构模型化原理2、因果关系图和流图库存量（2）流图符号⑤参数（量）明确问题及其构成要素；确定系统的边界；绘制要素间相互作用关系的因果关系图。注意一定要形成回路；确定变量类型（L变量、R变量和A变量）。将要素转化为变量，是建模的关键一步。在此，应考虑以下几个具体原则：a.水准（L）变量是积累变量，可定义在任何时点；而速率（R)变量只在一个时段才有意义。b.决策者最为关注和需要输出的要素一般被处理成L变量。c.在反馈控制回路中，两个L变量或两个R变量不能直接相连。d.为降低系统的阶次，应尽可能减少回路中L变量的个数。故在实际系统描述中，辅助（A）变量在数量上一般是较多的。绘制SD流图。注意：在绘制流图时，应特别注意形成正确的回路和用好信息连接线，并注意不要把不同的实物流直连在一起(参见下例)。举例P举例举例：商店库存问题举例：商店库存问题三、基本反馈回路的DYNAMO仿真分析水准方程（流位方程，L方程）LL1·K=L1·J+DT*(RI·JK-RO·JK)注：①方程式“=”左端第一项的变量名称必须和右端第一项名称相同；Eg:I.K=DT*(XIN.JK-Xou.JK)+I.J×②L方程等号右边第二项至少包含一个速率变量（流率变量）③凡是用L方程式定义的变量名称，只能用N方程式赋初值，不设定将自动取零。Eg:NP=10√CP=10×④DT只能出现在L方程，不能出现在其它方程。速率方程（流率方程,R方程）RR1·KL=f(L1·K,A1·K,…)注意：①R方程无固定表达式，其右边是根据实际情况写，但不允许含有时间间隔DT(速率的值在DT内不变);Eg:ROR.KL=(DI-I.K)/AT②R方程没有速率变量；③速率的值是对正在进行计算的时刻K以后立即出现的KL间隔来说的。辅助方程（A方程）AA1·K=g(L1·K,A2·K,R1·JK,…)注意：①没有统一的标准格式；②时间表示总是K，无前一时刻（.J）的值，一旦出现过去时刻辅助变量，则被认为是环；Eg:AP.K=P.J+D.K×RPSR.KL=(IDR.K-IAR.K)/ATAIDR.K=(AIR)(RSR.K)√③可由现在时刻的其他变量（A,L,R）求出.赋初值方程（N方程）NL1=数值或L1=L10L10=数值注意：①初始条件方程仅在开始仿真时计算一次；②N方程式左、右两侧的变量都不加注时间符号常量方程（C方程）CC1=数值2、一阶正反馈回路3、一级负反馈回路4、二阶负反馈回路（简单库存控制系统的扩展）使用蒙特卡洛法对下面的报童问题进行仿真