系统工程-第四章-系统动力学仿真.ppt

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,进入,系统工程,（Systems Engineering SE）,第四章 系统仿真及系统动力学方法,1.,系统仿真概述,2.,系统动力学结构模型化原理,3.,基本反馈回路的,DYNAMO,仿真分析,4.,系统动力学模拟步骤,5.,常用软件,教学内容,1.系统仿真概述,2.系统动力学结构模型化原理,3.DYNAMO仿真分析,教学要求,1.熟悉系统仿真的相关概念,2.掌握系统动力学建模的原理、方法及步骤,3.熟练应用STELLA进行系统仿真,教学重点及难点,因果关系图及流图的绘制,DYNAMO,仿真,一.系统仿真及系统动力学概述,（一）概念及作用,1.,基本概念,所谓系统仿真，就是,根据系统分析的目的,，,在分析系统各要素性质及其相互关系的基础上,，建立能,描述系统结构或行为过程的、且具有一定逻辑关系或数量关系的仿真模型,，据此进行试验或定量分析，以获得正确决策所需的各种信息。,2、系统仿真的实质,(1)它是一种对系统问题求数值解的计算技术。,尤其当系统无法通过建立数学模型求解时,，,仿真技术能有效地来处理。,(2)仿真是一种人为的试验手段。它和现实系统,实验的差别在于，仿真实验不是依据实际环,境，而是作为实际系统映象的系统模型以及,相应的,“人造”环境,下进行的。这是仿真的主,要功能。,(3)仿真可以比较真实地描述系统的运行、演变,及其发展过程。,3、系统仿真的作用,(1)仿真的过程也是实验的过程，而且还是系统地收集和积累信息的过程。尤其是对一些复杂的随机问题，应用仿真技术是提供所需信息的唯一令人满意的方法。,(2)对一些难以建立物理模型和数学模型的对象系统，可通过仿真模型来顺利地解决预测、分析和评价等系统问题。,（二）系统仿真方法,系统仿真的基本方法是建立系统的结构模型和量化分析模型，并将其转换为适合在计算机上编程的仿真模型，然后对模型进行仿真实验。,由于连续系统和离散(事件)系统的数学模型有很大差别，所以系统仿真方法基本上分为,两大类,，即,连续系统仿真方法,和,离散系统仿真方法。在以上两类基本方法的基础上，还有一些用于系统(特别是社会经济和管理系统)仿真的特殊而有效的方法，如系统动力学方法、蒙特卡洛法等。,系统动力学方法通过建立系统动力学模型(流图等)、利用DYNAMO仿真语言在计算机上实现对真实系统的仿真实验，从而研究系统结构、功能和行为之间的动态关系。,1、背景,80年代以来,1956年至,60年代初,60年代初至,70年代初,70年代初至,80年代,SD的出现始于1956年，主要应用于工业企业管理，并创立了“Industrial Dynamics”(1959),SD思想和方法的应用范围日益扩大。“Principles of Systems”(1968)，“Urban Dynamics”(1969)的出现.,1972年美国MIT的J.W.Forrester正式提出“Systems Dynamics”。经历了两次严峻的挑战。,SD成为一种重要的系统工程方法论和重要的模型方法。尤其是随着国内外管理界对学习型组织的关注,SD思想和方法的生命力更为强劲。,系统动力学发展,罗马俱乐部：,国际性的未来学研究团体。,1968,年,4,月在意大利经济学家,A.,佩切伊和英国科学家,A.,金倡议下，于罗马成立。宗旨是研究未来的科学技术革命对人类发展的影响，阐明人类面临的主要困难以引起政策制订者和舆论注意。会员限,300,名。现有,100,多名国际上著名的学者和社会活动家为个人会员。每年召开一次大会，并经常召开国际性学术会议。出版了,增长的极限,（即,米都斯报告,）、,重建国际秩序,、,走出浪费的时代,、,人类的目的,、,学无止境,、,第三世界：世界的四分之三,、,关于财富和福利的对话,、,走向未来的道路图,等著作。,The Limits to Growth,简介,：地球是人类目前唯一赖以身存的星球，但是，人类的发展却总是给滋养哺育他的大地带来无尽的折磨和无法修复的毁坏。全球气候变暖、海平面上升、人口的暴涨、土地沙漠化种种迹象表明，人类正在为自己的所作所为付出代价。,增长的极限,曾经是遥远的未来，但今天它们已经广泛存在。崩溃的概念曾经被认为是不可思议的，但今天它已经进入公众的谈论话题,The Limits to Growth,本书可以说是第一次向人们展示了在一个有限的星球上无止境地追求增长所带来的后果。这本震惊了世界并畅销全球的书，在今天，带着30年来新增的数据，再次就人类对气候、水质、鱼类、森林和其他濒危资源的破坏敲晌了警钟,。,第二次挑战,(70,年代初到,80,年代中,),：,Forrester,教授在多方资助之下开始研究美国全国模型，解开了一些在经济方面长期存在、令经济学家困惑不解的疑团。,诸如，,70,年代以来的通货膨胀、失业率和实际利率同时增长等问题。其最有价值的研究成果还在于揭示了美国与西方国家经济长波,(Long Wave),形成的内在奥秘。由于在全国模型与长波理论研究方面取得成就，使系统动力学这一门学科在理论和应用研究两方面都取得了飞跃性进展。从此，系统动力学进入了蓬勃发展时期。,1972年正式定名系统动力学：,“系统动力学是,研究信息反馈系统动态行为的计算机仿真方法,。它有效地,把信息反馈的控制原理与因果关系的逻辑分析结合起来,，面对复杂实际问题，,从研究系统的内部结构入手,，建立系统的仿真模型，并对模型实施各种不同的政策方案，通过计算机仿真展示系统的宏观行为，寻求解决问题的正确途径。”,美,彼得圣吉,（PeterMSenge）,著,，,第五项修炼学习型组织的艺术与实务,，,上海三联书店，1998。,作者简介：1970年从斯坦福大学获工学学士后进入MIT攻读管理硕士学位，在此期间被Forrester教授的SD整体动态搭配的管理新理念所吸引。1978年获得博士学位后，一直和MIT的工作伙伴及企业界人士一道，孜孜不倦地致力于将SD与组织学习、创造原理、认知科学等融合，发展出一种人类梦寐以求的组织蓝图学习型组织。,彼得.圣吉提出“五项修炼”,培养成员的自我超越意识：坚持不懈地实现心中的渴望；,改善心智模式：改变人们对世界的看法；,建立共同愿景：建立组织共同追求的崇高理想；,搞好团队学习：成员间的心灵的相互交融和感悟；,运用系统思考：真正认请自己在系统中的位置和作用。,“系统思考”是“学习型组织”理论的第五部分又是它的核心，为了突出它，该书定名为,第五项修炼,。,2、研究对象及其结构特点,（,1）研究对象社会系统,（2）SD将社会系统当作,非线性（多重）信息反,馈系统,来研究,（3）结构特点,抉择性具有决策环节（人、信息）,自律性具有反馈环节,非线性具有延迟环节,系统动力学模型的特点,(1),多变量,。主要是由SD对象系统的动态特性和复杂性所决定的。SD模型有三种基本变量、五到六种变量。,(2),定性分析与定量分析相结合,。SD模型由结构模型(流图)和数学模型(DYNAMO方程)所组成。,(3),以仿真实验为基本手段和以计算机为工具,。SD作为一种计算机仿真分析方法，是实际系统的“实验室”，可在PD-plus、VENSIM等软件支持下来运行。,(4),可处理高阶次、多回路、非线性的时变复杂系统问题。,认识,问题,界定,系统,要素及其因果关系分析,建立结构模型,建立量化分析模型,仿真,分析,比较与评价,政策,分析,初步分析,规范分析,综合分析,SD工作程序图,4、工作程序,四个基本要素,状态、信息、决策、行动,两个基本变量,水准变量（L）、速率变量（R）,一个基本思想,反馈控制,2、因果关系图和流图,（1）因果关系图（因果反馈回路）,因果箭因果链因果（反馈）回路,利息,(元/年),银行,货币,利率,+,+,(+),正关系,若满足下列条件之一：,A加到B中；A是B的乘积因子；,A变到AA，有B变到BB，即A、B的变化方向相同。,则称A到B具有正因果关系，简称正关系，用“”号标在因果链上。,人口出生率,人口总数,利息,(元/年),银行,货币,利率,+,+,(+),下一年的销售增长,每年的年收入,增长速率,+,+,(+),恶性循环,良性循环,负关系,若满足下列条件之一：,A从B中减去；1/A是B的乘积因子；,A变到AA，有B变到BB，即A、B的变化方向相反。,则称A到B具有负因果关系，简称负关系，用“”号标在因果链上。,人口死亡率,人口总数,反馈,从控制论的观点看，任何一个具有使自身内部保持稳定的系统，都具有某种反馈机制。,反馈（Feedback）:构成系统的某一成分的输出与输入之间的关系，或者说是输出变成了决定系统未来功能的输入,。,生物群落的稳定性,例1：,例2：,例3：,天敌害虫种群系统,天敌数量,（输入）,害虫数量,（输出）,害虫环境系统,初春气温,（输入）,害虫数量,（输出）,害虫作物系统,害虫数量,（输入）,作物长势,（输出）,反馈环,反馈环,反馈环,库存量,订货量,库存差额,期望,库存,+,+,-,（,-,）,负反馈系统实例,库存,发货单,期望库存,定货,途中货物量,收货,一个简单的库存控制系统,：,（）,对服务质量的重视程度,有效的读者意见,读者意见,+,+,-,（,-,）,偏见,负反馈系统实例,人口,分配,生产,消费,流通,消费水平,资源消耗或占用量,资源存量,+,-,+,+,+,+,+,+,+,+,-,(-),(-),2,3.流图绘制程序和方法,明确问题及其构成要素,；,绘制要素间相互作用关系的因果关系图,。注意一定要形成回路,；,确定变量类型（L变量、R变量和A变量）,。将要素转化为变量，是建模的关键一步,。,绘制SD流图,。,在绘制流图时，应特别注意形成正确的回路和用好信息连接线，并注意不要把不同的实物流直连在一起.,a.水准（L）变量是积累变量，可定义在任何时间点；而速率（R)变量只 在一个时段才有意义。,b.决策者最为关注和需要输出的要素一般被处理成L变量。,c.在反馈控制回路中，两个L变量或两个R变量不能直接相连。,d.为降低系统的阶次，应尽可能减少回路中L变量的个数。故在实际系统描述中，辅助（A）变量在数量上一般是较多的。,库存量,每周订货量,库存差额,+,+,-,（-）,期望,库存,库存量,水准变量,每周订货量,速率变量,库存差额,辅助变量,3、举例,库存系统动力学流程图,I,D,目标库存量,实际库存量,库存差,周订货量,R,利息,(元/年),银行,货币,利率,+,+,(+),出生,人口,人口,总量,死亡,人口,（平均）出生率,（平均）死亡率,（,-,）,-,+,+,（+）,组织绩效,组织改善,组织缺陷,（,-,）,-,+,+,3、举例,L1,R1（利息1）,C1(利率),I,R1(订货量),库存量,D,Y(期望库存),（库存差额）,P,R1,R2,(出生人口),(人口总量),(死亡人口),C1（出生率）,C2（死亡率）,组织改善,组织绩 效,组织,缺陷,。,例：,现分析某经营单一商品的零售店的订货策略问题，要求绘制系统动力学流程图。,分析：,由于零售店向顾客销售商品，使零售店的库存量不断减少，为了补充库存，店方要向生产厂家提出订货。接受订货的厂家计划生产该种商品以满足订货要求。这时零售店的库存量又相应增加。系统的边界可以定为由零售店和工厂两部分组成。,1.,确定系统边界,工厂,订货,供应,系统边界,顾客购货,销售,外生变量,或,扰动变量,零售店,2.,系统的组成要素,零售店,零售店的销售量；库存量；订货量,工厂,工厂未供订货量；生产量；生产能力；计划产量,3.,因果关系分析,零售店,销售,零售店,订货,零售店,库存,工厂未,供订货,计划,产量,工厂生,产能力,工厂,生产,商店,工厂,三、基于反馈回路的DYNAMO仿真分析,（一）DYNAMO简介,DYNAmic MOdel 的缩写。它将实际问题构造成具有反馈结构的动态模型，并通过计算机仿真得到该系统随时间变化的动态行为。,DYNAMO模型由两种语句组成：,方程式语句（直接用于仿真计算），,命令语句（用于控制仿真过程、输入输出）,变量的时间标注,DT,DT,JK,KL,前一时刻,现在时刻,下一时刻,J K L,（二）DYNAMO方程,状态（水准）方程,决策（速率）方程,辅助方程,初值方程,常数方程,（1）状态（水准）方程,表述了系统动力学模型中状态积累的过程,L LEVELKLEVEL J DT,（RIN JK,ROUT JK）,用L方程式定义的流位变量必须用初值方程式给定初值,例：,L POPKPOP J DT,（BIRTH JK,DEATH JK）,N POP10000,（2）决策（速率）方程,决策方程（R方程）,描述系统动力学中状态变化速率的方程。基本形式可表示为：,R RATEKLf（状态变量，辅助变量，常量）,速率R的值在K时刻进行计算，而在自K至L的时间间隔内（在DT内），假定保持不变,辅助说明速率变量或简化决策函数的方程。基本形式可表示为：,A AUX.K=g(A.K,L.K,R.JK,C,),时间标注总是K可由当前时刻的其他变量求出,（3）辅助方程（A方程）,（4）初值方程（N方程）,为参数或变量设定初始值，一般形式：,N 变量名称表达式，变量名，数值,仅在仿真过程中第一步运算时使用；左右两边的变量都不加时间标注；由N方程式定义的变量不能直接在重复运行中使用。,（5）常数方程（C方程）,给参数赋值，可以在重复运行中使用,C 变量名常数,DYNAMO中变量名的字符数不超过6个，而且第一个字符必须是字母。,除“*”及“RUN”语句以外，其他语句可以任意安排顺序。DYNAMO语言运行时有一个严格的仿真顺序。例：一阶正反馈回路,年人口,增 加,(+),PR,+,P,PR,P,C1（人口年自然增长率0.02）,。,。,L PK=PJ+DT*PRJK,N P=100,R PRKL=C1*PK,C C1=0.02,P,PR,0,100,2,1,102,2.04,2,104.04,2.0808,p,100,0,一阶正反馈（简单,人口问题）系统输,出特性曲线,人口数,+,例：,L,LC,K,=,LC,J+DT,*,RTJK,A,AJ,K,=,C,2,*,LCK,R,RT,KL=AJ,K,*,C,1,投资,产值,积累,+,+,+,源,工业产值,LC,RT,积累,AJ,积累率C,2,投资率C,1,一级负反馈回路,库存量,库存,差额,订货量,+,(),R1,D,I,+,期望库存Y,1000,Z,（订货调整时间，5）,I,R1,D,Y（6000）,。,。,。,。,L I,K=IJ+DT*R1JK,N I=1000,R R1,KL=D,K/Z,A D,K=Y-IK,C Z=5,C Y=6000,I,D,R1,0,1000,5000,1000,1,2000,4000,800,2,2800,3200,640,I,t,6000,0,一阶负反馈（简单,库存控制）系统输,出特性曲线,1000,4、系统动力学模拟步骤,5、常用软件,小 结,