《系统动力学方法》由会员分享,可在线阅读,更多相关《系统动力学方法(41页珍藏版)》请在金锄头文库上搜索。
1、系统动力学方法小组成员:牛头Dynamo 只需要理解即可,无需去学习系统动力学的由来与发展v系统动力学(Systems Dynamics, SD)是美国麻省理工学院(MIT)J.W.弗雷斯特(J. W. Forrester)教授最早提出的一种对社会经济问题进行系统分析的方法论和定性与定量相结合的分析方法。目的在于综合控制论、信息论和决策论的成果,以电子计算机为工具,分析研究信息反馈系统的结构和行为。 v1972年正式提出“Systems Dynamics” 系统动力学的由来与发展vSD的出现始于20世纪50年代后期,当时,主要应用于工商企业管理,处理诸如生产与雇员情况的波动、企业的供销、生产与
2、库存、股票与市场增长的不稳定性等问题,并创立“Industrial Dynamics” (1959)。系统动力学的由来与发展v在整个60年代,动力学思想与方法的应用范围日益扩大,几乎遍及各类系统,深入到各种领域。作为方法论基础,出现了“Principles of Systems(1968)”, 总结美国城市兴衰问题的理论与应用研究成果的“Urban Dynamics(1969)”和著名的“World Dynamics (1971)”等也是J.W.弗雷斯特等人的重要成就。 系统动力学的由来与发展v70年代以来,SD经历两次严峻的挑战并走向世界,进入蓬勃发展时期。 v第一次挑战(70年代初到70年
3、代中):SD与罗马俱乐部一起闻名于世,走向世界,其主要标志是两个世界模型(WORLD,)的研制与分析 v第二次挑战(70年代初到80年代中):对美国全国SD模型的研制和对美国与整个西方国家经济长波(Long Wave)问题的研究 系统动力学的由来与发展v近年来SD正在成为一种新的系统工程方法论和重要的模型方法,渗透到许多领域,在国土规划、区域开发、环境治理和企业战略研究等方面,正显示出它的重要作用。v尤其是随着国内外管理界对学习型组织的关注,SD思想和方法的生命力更为强劲。v但目前应更加注重SD的方法论意义,并注意其定量分析手段的应用场合及条件。SD的研究对象及其特点 vSD的研究对象主要是社
4、会(经济)系统。 v(1)社会系统中存在着决策环节v社会系统的行为总是经过采集信息,并按照某个政策进行信息加工处理作出决策后出现的,决定是一个经过多次比较、反复选择、优化的过程。v对于大规模复杂的社会系统来说,其决策环节所需要的信息的信息量是十分庞大的。其中既有看得见、摸得着的实体,又有看不见、摸不到的价值、伦理、道德观念及个人、团体的偏见等因素。SD的研究对象及其特点v(2)社会系统具有自律性 v自律性就是自己作主进行决策,自己管理、控制、约束自身行为的能力和特性。 工程系统是由于导入反馈机构而具有自律性的;社会系统因其内部固有的“反馈机构”而具有自律性。因此,研究社会系统的结构,首先(也是
5、最重要的)就在于认识和发现社会系统中所存在着的由因果关系形成的反馈机构。SD的研究对象及其特点v(3)社会系统的非线性v指社会现象中原因和结果之间的极端非线性关系:原因和结果在时、空上的分离性、出现事件的意外性、难以直观性等。v高度非线性是由于社会问题的原因和结果相互作用的多样性、复杂性造成的。具体来说,一方面是由于社会问题的原因和结果在时间、空间上的滞后,另一方面是由于社会系统具有多重反馈结构。SD的研究对象及其特点vSDSD将社会系统当作将社会系统当作非线性(多重)信息反非线性(多重)信息反馈系统馈系统来研究来研究vSDSD将社会经济问题模型化,对社会经济现将社会经济问题模型化,对社会经济
6、现象进行预测、对社会系统结构和行为进行象进行预测、对社会系统结构和行为进行分析,为企业、地区、国家、国际制定发分析,为企业、地区、国家、国际制定发展战略、进行决策,提供有用的信息。展战略、进行决策,提供有用的信息。 SD的模型特点v(1)多变量。这主要是由SD对象系统的动态特性和复杂性所决定的。v(2)定性分析与定量分析相结合。SD模型由结构模型(流图)和数学模型(DYNAMO方程)所组成 v(3)以仿真实验为基本手段和以计算机为工具。 v(4)可处理高阶次、多回路、非线性的时变复杂系统问题。 SD的一般工作过程 v明确系统建模目的-认识问题v确定系统边界-假定:封闭社会系统,系统行为来自系统
7、内部v因果关系分析-系统内部要素之间的因果关系分析,用以明确要素之间的因果关系,并用表示因果关系的反馈回路来描述v建模-a.流程图根据反馈回路,用专门设计的描述各种变量的符号绘制而成;b.结构方程式,用专门的DYNAMO语言建立,用以进行定量分析和计算机仿真实验。SD的一般工作过程v计算机仿真实验v结果分析v模型修正v比较与评价v政策分析认识认识问题问题界定界定系统系统要素及其因要素及其因果关系分析果关系分析建立结建立结构模型构模型建立数建立数学模型学模型仿真仿真分析分析比较与比较与评价评价政策政策分析分析(流图)(流图)(DYNAMOY方程)方程)SDSD结构模型化原理结构模型化原理 决策决
8、策信息信息行动行动系统系统状态状态速率变量速率变量水准变量水准变量信息信息 流流(行动)(行动)(Rate)(Level)四个基本要素四个基本要素状态、信息、决策、行动状态、信息、决策、行动两个基本变量两个基本变量水准变量(水准变量(L L)、)、速率变量(速率变量(R R)一个基本思想一个基本思想反馈控制反馈控制因果关系分析因果关系分析是构成系统动力学模型的基础,是社会系统内部关系的真实写照。因果关系通过因果关系图来描述l因果箭因果关系也称影响关系,可通过连接因果要素的“因果箭”表示。因果关系的极性(正,负)表示因果的变化方向是否一致l因果链根据因果关系的传递性;因果链的极性:因果链的极性符
9、号与因果箭极性乘积符号相同AB+ABAB+CDABCDE-+-反馈回路原因与结果之间的相互作用在自然现象中非常普遍,从而形成因果关系的反馈回路反馈回路的基本特征:原因和结果具有相对性反馈回路的极性:正反馈回路自我强化(自我弱化)负反馈回路自我调节(内部稳定)系统动力学认为:系统性质和行为完全取决于系统中存在的反馈回路;系统的结构主要就是系统中反馈回路的结构。 库存量库存量订货量订货量库存差额库存差额期望期望库存库存+-( -)利息利息(元元/年年)银行货币利率利率+(+)多重反馈回路在复杂的社会系统中存在的两个或两个以上的反馈回路,称做多重反馈回路。这些反馈回路可能存在相互促进或制约的作用,可
10、能在不同时期表现出不同的主导程度。出生出生人口人口人口人口总量总量死亡死亡人口人口(平均)死亡率(平均)死亡率(-)-+(+)(平均)出生率(平均)出生率流图符号流图符号 1 1 流流 实物流实物流信息流信息流2 2 速率变量速率变量 3 3 水准变量水准变量 L14 4 辅助变量辅助变量A1R1R1()。5 5 参数(量)参数(量) (常量)(常量)(初值)(初值)6 6 源与洞源与洞 7 7 信息的取出信息的取出L1。A1。8 滞后或延迟滞后或延迟(Delay) 由于信息和物质运动需要一定的时间,于是就带来原由于信息和物质运动需要一定的时间,于是就带来原因和结果、输入和输出、发送和接收等之
11、间的时差,因和结果、输入和输出、发送和接收等之间的时差,并有物流和信息流滞后之分。并有物流和信息流滞后之分。 输入率输出率DELAY1,3。延迟时间SMOOTH平滑变量。DLINF3延迟变量。A,L或R平滑时间A,L或R延迟时间流图绘制程序和方法流图绘制程序和方法v 明确问题及其构成要素;明确问题及其构成要素;v 绘制要素相互作用的因果关系图。一定要形成回路。绘制要素相互作用的因果关系图。一定要形成回路。v 确定变量类型(确定变量类型(L L变量、变量、R R变量和变量和A A变量)。应考虑以变量)。应考虑以下几个具体原则:下几个具体原则:va. a. 水准(水准(L L)变量是积累变量,可定
12、义在任何时点;而变量是积累变量,可定义在任何时点;而速率(速率(R)R)变量只在一个时段才有意义。变量只在一个时段才有意义。vb. b. 决策者最关注和需输出的要素一般被处理成决策者最关注和需输出的要素一般被处理成L L变量。变量。 vc. c. 在反馈回路中两个在反馈回路中两个L L变量或两个变量或两个R R变量不能直接相连变量不能直接相连 。vd. d. 为降低系统的阶次,应尽可能减少回路中为降低系统的阶次,应尽可能减少回路中L L变量的个变量的个数。因此,辅助(数。因此,辅助(A A)变量在数量上一般是较多的。变量在数量上一般是较多的。v 绘制绘制SDSD流图。流图。流图举例流图举例L1
13、R1(利息1)C1(利率)IR1(订货量)库存量DY(期望库存)(库存差额)PR1R2(出生人口)(人口总量) (死亡人口)C1(出生率)C2(死亡率)组织改善组 织绩效组织缺陷。基本基本DYNAMODYNAMO方程方程JKKLJKLDTDTSD的对象系统是随时间连续变化的动态系统。在DYNAMO方程中变量一般带有时间标号 ,如右图J 表示过去时刻K 表示现在时刻L 表示将来时刻JK表示由过去时刻到现在时刻的时间间隔KL表示由现在时刻到未来时刻的时间间隔DT仿真步长变量名:字符数不超过6个,第一个为字母保留字:L 、R、A、S、N、C、T、PRINT、PIOT、RUN、NOISE、NOTE基本
14、基本DYNAMODYNAMO方程方程v水准方程(水准方程(L L方程)方程) L L1.K=L1.J+DT*(RI.JK-L L1.K=L1.J+DT*(RI.JK-RO.JK)RO.JK)v速率方程(速率方程(R R方程)方程) R R1.KL=f ( L1.K,A1.K,)R R1.KL=f ( L1.K,A1.K,)v辅助方程(辅助方程(A A方程)方程) A A A1.K=g(L1.K,A2.K,R1.JK,)A1.K=g(L1.K,A2.K,R1.JK,)v赋初值方程(赋初值方程(N N方程)方程) N L1=N L1=数值数值 或或 N L1=L10N L1=L10 C L10=L
15、10=数值数值v常量方程常量方程 (C C方程)方程) C C1=C C1=数值数值典型反馈回路及其仿真计算 v1、一阶正反馈回路(以人口的增加机理为例) 年人口增加人口数(+)P+PRPRPC1(0.02)典型反馈回路及其仿真计算vL P.K=PJ+DT*PR.JKvN P=100vR PR.KL=C1*R.KvC C1=0.02PPR0100211022.042104.042.0808p1000简单人口问题简单人口问题典型反馈回路及其仿真计算v2、一阶负反馈回路(以简单库存系统为例) 库存量库存差额订货量+()R1DI+期望库存Y1000Z(订货调整时间,5)IR1DY(6000)典型反馈
16、回路及其仿真计算v数学模型及仿真计算L I.K=I.J+DT*R1.JK N I=IO C IO=1000 R R1.KL=D.K/Z A D.K=Y-I.K C Z=5 C Y=6000 IDR10100050001000120004000800228003200640334402560512.It10000一阶负反馈回路一阶负反馈回路典型反馈回路及其仿真计算v3、二阶负反馈回路库存量入库量途中存货量订货量库存差额I()+()GR2+R1+Z(5)Y(6000)IG100001000R1R2DW(10)简单库存控制的数学模型(01) arrival rate=in-transit/arriv
17、al timeUnits: *undefined*(02) arrival time=10Units: *undefined*(03) expect inventory=6000Units: *undefined*(04) FINAL TIME = 100Units: DayThe final time for the simulation.(05) in-transit= INTEG ( +order rate-arrival rate, 10000)Units: *undefined*(06) INITIAL TIME = 0Units: DayThe initial time for t
18、he simulation.(07) inventory= INTEG (arrival rate, 1000)Units: *undefined*(08) order rate= shortage/order timeUnits: *undefined*(09) order time= 5Units: *undefined*(10) SAVEPER = TIME STEP Units: Day 0,?The frequency with which output is stored.(11) shortage= expect inventory-inventoryUnits: *undefined*(12) TIME STEP = 1Units: Day 0,? The time step for the simulation.基本介绍结束谢谢大家!