系统工程6系统仿真及SD方法.ppt
系统工程,第六讲 系统仿真及SD方法,经济管理学院工业工程系,系统仿真概述系统动力学结构模型化原理基本反馈回路的DYNAMO仿真分析,系统仿真及SD方法,6.1 系统仿真,所谓系统仿真，就是根据系统分析的目的，在分析系统各要素性质及其相互关系的基础上，建立能描述系统结构或行为过程的、且具有一定逻辑关系或数量关系的仿真模型，据此进行试验或定量分析，以获得正确决策所需的各种信息。,一、概念及作用,1.基本概念,2、系统仿真的实质,6.1 系统仿真,(1)它是一种对系统问题求数值解的计算技术。尤其当系统无法通过建立数学模型求解时，仿真技术能有效地来处理。(2)仿真是一种人为的试验手段。它和现实系统实验的差别在于，仿真实验不是依据实际环境，而是作为实际系统映象的系统模型以及相应的“人造”环境下进行的。这是仿真的主要功能。(3)仿真可以比较真实地描述系统的运行、演变及其发展过程。,6.1 系统仿真,3、系统仿真的作用,(1)仿真的过程也是实验的过程，而且还是系统地收集和积累信息的过程。尤其是对一些复杂的随机问题，应用仿真技术是提供所需信息的唯一令人满意的方法。(2)对一些难以建立物理模型和数学模型的对象系统，可通过仿真模型来顺利地解决预测、分析和评价等系统问题。,(3)通过系统仿真，可以把一个复杂系统降阶成若干子系统以便于分析。(4)通过系统仿真，能启发新的思想或产生新的策略，还能暴露出原系统中隐藏着的一些问题，以便及时解决。,6.1 系统仿真,6.1 系统仿真,二、系统仿真方法,系统仿真的基本方法是建立系统的结构模型和量化分析模型，并将其转换为适合在计算机上编程的仿真模型，然后对模型进行仿真实验。由于连续系统和离散(事件)系统的数学模型有很大差别，所以系统仿真方法基本上分为两大类，即连续系统仿真方法和离散系统仿真方法。,6.1 系统仿真,在以上两类基本方法的基础上，还有一些用于系统(特别是社会经济和管理系统)仿真的特殊而有效的方法，如系统动力学方法、蒙特卡洛法等。系统动力学方法通过建立系统动力学模型(流图等)、利用DYNAMO仿真语言在计算机上实现对真实系统的仿真实验，从而研究系统结构、功能和行为之间的动态关系。,6.2 系统动力学原理,1、由来与发展,Systems Dynamics,SD/J.W.Forrester(MIT)Industridl Dynamics(ID),1959Principles of Systems,1968Urban Dynamics(UD),1969World Dynamics(WD),1971SD,1972,6.2 系统动力学原理,2、研究对象及其结构特点,（1）研究对象社会系统（2）结构特点 抉择性具有决策环节（人、信息）自律性具有反馈环节 非线性具有延迟环节（3）SD将社会系统当作非线性(多重)信息反馈系统来研究,6.2 系统动力学原理,（流图）（DYNAMOY方程）,3、工作程序,6.2 系统动力学原理,4、系统动力学模型,（1）常用要素流速率水平变量源与汇参数,（2）流图符号,流,实物流,信息流,速率变量,水准变量,L1,辅助变量,A1,。,R1,R1,(),。,6.2 系统动力学原理,明确问题及其构成要素；绘制要素间相互作用关系的因果关系图。注意一定要形成回路；确定变量类型（L变量、R变量和A变量）。将要素转化为变量，是建模的关键一步。在此，应考虑以下几个具体原则：,（3）流图绘制程序和方法,6.2 系统动力学原理,a.水准（L）变量是积累变量，可定义在任何时点；而速率（R)变量只在一个时段才有意义。b.决策者最为关注和需要输出的要素一般被处理成L变量。c.在反馈控制回路中，两个L变量或两个R变量不能直接相连。d.为降低系统的阶次，应尽可能减少回路中L变量的个数。故在实际系统描述中，辅助（A）变量在数量上一般是较多的。绘制SD流图。,6.2 系统动力学原理,5、举例,6.2 系统动力学原理,1、基本DYNAMO方程(DYNAmic Model),水准方程（L方程）L L1K=L1J+DT*(RIJK-ROJK)速率方程（R方程）R R1KL=f(L1K,A1K,)辅助方程（A方程）A A1K=g(L1K,A2K,R1JK,)赋初值方程（N方程）N L1=数值 或 L1=L10 L10=数值常量方程（C方程）C C1=数值,6.3 基本反馈回路的DYNAMO仿真分析,2、一阶正反馈回路,。,。,L PK=PJ+DT*PRJKN P=100R PRKL=C1*RKC C1=0.02,6.3 基本反馈回路的DYNAMO仿真分析,3、一级负反馈回路,。,。,。,。,L IK=IJ+DT*R1JKN I=1000R R1KL=DK/ZA DK=Y-IKC Z=5C Y=6000,4、简单库存控制系统的扩展,。,。,。,。,。,。,L GK=GJ+DT*(R1KL-R2JK)L IK=IJ+DTR2JKR R1KL=D/ZA D=Y-IKC Y=6000C W=10,Z=5C I=1000C G=10000,DELAY,。,。,。,。,。,。,。,。,。,。,。,
