系统动力学模型部分集.doc

第10章系统动力学模型系统动力学模型(System Dynamic)是社会、经济、规划、军事等许多领域进行战略研究的重要工具，如同物理实验室、化学实验室一样，也被称之为战略研究实验室，自从问世以来，可以说是硕果累累。 1 系统动力学概述 2 系统动力学的基础知识 3 系统动力学模型 第1节 系统动力学概述 1.1 概念 系统动力学是一门分析研究复杂反馈系统动态行为的系统科学方法，它是系统科学的一个分支，也是一门沟通自然科学和社会科学领域的横向学科，实质上就是分析研究复杂反馈大系统的计算仿真方法。 系统动力学模型是指以系统动力学的理论与方法为指导，建立用以研究复杂地理系统动态行为的计算机仿真模型体系，其主要含义如下： 1 系统动力学模型的理论基础是系统动力学的理论和方法； 2 系统动力学模型的研究对象是复杂反馈大系统； 3 系统动力学模型的研究内容是社会经济系统发展的战略与决策问题，故称之为计算机仿真法的“战略与策略实验室”； 4 系统动力学模型的研究方法是计算机仿真实验法，但要有计算机仿真语言DYNAMIC的支持，如：PD PLUS，VENSIM等的支持； 5 系统动力学模型的关键任务是建立系统动力学模型体系； 6 系统动力学模型的最终目的是社会经济系统中的战略与策略决策问题计算机仿真实验结果，即坐标图象和二维报表； 系统动力学模型建立的一般步骤是：明确问题，绘制因果关系图，绘制系统动力学模型流图，建立系统动力学模型，仿真实验，检验或修改模型或参数，战略分析与决策。 地理系统也是一个复杂的动态系统，因此，许多地理学者认为应用系统动力学进行地理研究将有极大潜力，并积极开展了区域发展，城市发展，环境规划等方面的推广应用工作，因此，各类地理系统动力学模型即应运而生。 1.2 发展概况 系统动力学是在20世纪50年代末由美国麻省理工学院史隆管理学院教授福雷斯特（JAY.W.FORRESTER）提出来的。目前，风靡全世界，成为社会科学重要实验手段，它已广泛应用于社会经济管理科技和生态灯各个领域。福雷斯特教授及其助手运用系统动力学方法对全球问题，城市发展，企业管理等领域进行了卓有成效的研究，接连发表了工业动力学，城市动力学，世界动力学，增长的极限等著作，引起了世界各国政府和科学家的普遍关注。 在我国关于系统动力学方面的研究始于1980年，后来，陆续做了大量的工作，主要表现如下： 1）人才培养 自从1980年以来，我国非常重视系统动力学人才的培养，主要采用“走出去，请进来”的办法。请进来就是请国外系统动力学专家来华讲学，走出去就是派留学生，如：首批派出去的复旦大学管理学院的王其藩教授等，另外，还多次举办了全国性的讲习班。 2）编译编写专著 组织专家编译了工业动力学，城市动力学等。 编写专著有：王其藩著系统动力学，高级系统动力学；胡玉奎著系统动力学，王洪斌著系统动力学教程，贾仁安著系统动力学教程等。 3）引进专业软件 引进的软件有：MICRO-DYNAMO，DYNAMAP2，DYNAMO，STELLA，PD PLUS等，近几年又引进的最先进实用的VENSIM专业软件。并自行研制了一些专用软件。 4）新设课程 新开设了系统动力学专业课程。在几十所大学的管理系或管理学院以及科研单位的研究生开设了系统动力学课程。 5）组织机构与学术会议 于19 年成立了全国系统动力学委员会。组建了一些专门研究机构和教学机构。开展了许多专项研究工作。建立了国家总体系统动力学模型，省和地区的发展战略研究系统动力学模型，省级能源，环境预测系统动力学模型及科技，工业，农业林业等行业发展战略研究系统动力学模型等。 1986年8月，在上海召开的“全国系统动力学学术研讨会“上，140多名代表提交了95篇有关系统动力学理论和应用研究方面的论文。1987年6月，在上海召开的国际学术会议上我国代表交流了29篇论文，占会议论文数的45%。1988年7月，美国圣迭戈召开了国际学术年会，我国有十名代表参加，交流论文十多篇。1989年7月，在西德斯图加特召开的国际学术年会上，我国学者交流论文14篇，有4人参加会议。 目前，在我国系统动力学已经发展成熟，并正向深入和全面应用延伸，形成了一支强大的研究力量，发展趋势看好，有理由相信，系统动力学必将在我国社会，经济，科技，管理和生态等领域的研究中发挥更大作用。 第2节 系统动力学的基础知识 系统动力学模型建立的基本知识，基本原理主要有：因果关系图，模型流图及模型的组成等。现分别介绍。 2.1 因果关系 1 因果关系 因果关系是指由原因产生某结果的相互关系。从哲学角度讲，原因和结果是揭示客观事物的因果联系的重要哲学概念，它们是客观事物普遍联系和相互作用的表现形式之一。原因是某种事物或现象，是造成某种结果的条件；结果是原因所造成的事物或现象，是在一定阶段上事物发展所达到的目标状态。 通常用箭头线来表示，它有正因果关系和负因果关系两种，如图91。 原因 结果 + 就业机会E 迁入人口数I - 死亡率R 总人口数P 正因果关系：两个变量呈同方向变化趋势，如：E增加，I增加； E减少，I减少。 负因果关系：两个变量呈异方向变化趋势，如：R增加，P减少； R减少，P增加。 2）因果关系环图 因果关系环图是指由两个或两个以上的因果关系连接而成的闭合回路图示。它定性描述了系统中变量之间的因果关系。它有正负因果关系环图两种，如图93，图9-4所示： 正因果关系环图：它会引起系统内部活动加强。 准则：若各因果关系均为正，则该环为正因果关系环； 若各因果关系为负的个数是偶数时，则该环也为正因果关系环。 负因果关系环图：它会引起系统内部活动减弱。 准则：若各因果关系均为负，则该环为负因果关系环； 若因果关系为负的个数是奇数，则该环为负因果关系环。 再如：生态学人口增长因果关系环图，如图95，图9-6 所示： 2.2 系统动力学模型流图 系统动力学模型流图简称SD流图，是指由专用符号组成用以表示因果关系环中各个变量之间相互关系的图示。它能表示出更多系统结构和系统行为的信息，是建立SD模型必不可少的环节，对建立SD模型起着重要作用。其专用符号主要有八个： 1）水平变量水平变量符号是表示水平变量的积累状态的符号，它是SD模型中最主要的变量。它由五部分组成，即：输入速率，输出速率，流线，变量名称及方程代码（L），如图 所示。2）速率变量速率变量符号是表示水平变量变化速率的变量。它能控制水平变量的变化速度，是可控变量。它由三部分组成，即：输入信息变量，变量名称及方程代码（R）。如图 所示。3）辅助变量辅助变量符号是辅助水平变量等的变量。如图 所示。4）外生变量外生变量符号如图 所示。5）表函数表函数符号如图 所示。6）常数常数符号如图 所示。7）流线流线符号又有物质流线，信息流线，资金流线，及订货流线四种：物质流线符号是表示系统中流动着的实体，如图 所示。信息流线符号是表示联接积累与流速的信息通道，如图 所示。资金流线符号是表示资金，存款及货币的流向，如图 所示。订货流线符号是表示订货量与需求量的流向，如图 所示。8）源与沟源符号与沟符号如图 所示。2.3 系统动力学模型系统动力学模型是由六种基本方程和专门的输出语句组成。其六种方程的标志符号分别为：L：水平变量方程； R：速率变量方程；A：辅助变量方程； N：计算初始值方程；C：赋值予常数方程； T：赋值予表函数中Y坐标值。L方程是积累方程；R，A方程是代数运算方程；C，T，N方程是提供参数值方程，并在同一次模拟中其值保持不变。1）L方程L方程是计算水平变量积累值的方程，其一般表示形式为：L 其中， L：水平变量方程代码，表示方程性质。 DT：时间间隔，即时间增量。.J：表示前一刻。 .K：现在时刻。 .L：未来一时刻。 ：过去一时刻人口数。 ：现在时刻人口数。 ：未来一时刻人口数。 ：过去至现在该段时刻的人口出生率。 ：过去至现在该时刻段的人口死亡率。 积累是系统内部流的堆积量，它等于过去一时刻的积累加上积累变动量，即变动增量。积累变动量是时间间隔与输入流速和输出流速之差的乘积。2）R方程R方程是计算单位时间流量的方程，即流速或速率。其一般表示形式为： R R R R 其中，：过去至现在时刻的出生率，单位（人/年）；：过去至现在时刻的死亡率，单位（人/年）； ：现在至未来时刻的出生率；单位（人/年）； ：现在至未来时刻的死亡率，单位（人/年）； ： 出生系数，单位（人/年人）； ： 死亡系数，单位（人/年人）； ：过去时刻人口总数； ：现在时刻人口总数。 3）A方程 A方程是辅助变量方程，用于对辅助变量赋值，其一般表示形式为： A 其中， ：表示现在人口总数。：求和函数，表示求算现在22个年龄组的总和。 4）N方程 N方程是变量初始值方程，表示对变量赋初始值，起一般表示形式为： N 其中，：表示各年龄组人口初始值。：是表函数，表示存储22个年龄组的初始值。 5）T方程 T方程是表函数方程，表示对相应的纵坐标Y赋值。 6）C方程 C方程是常数方程，表示对常数变量赋值。 第3节 系统动力学模型 系统动力学模型应用分析的一般步骤为： 1 明确问题 明确的问题是： 系统的范围：空间范围，如安徽省区域； 时间范围，如1961年 - 2050年； 时间间隔，DT=1年，等等。 解决途径： 计算机仿真实验。 数据资料： 人口总数，出生率，死亡率，自然增长率等。 2 明确目标 人口总数变化趋势； 自然增长率控制目标； 出生率控制目标； 死亡率控制目标等。 3 绘制系统流图 1）因果关系环图 主要变量清单，即列出主要变量的清单，以利于因果关系环流图的绘制。如： 总人口数，出生率，死亡率，出生系数，死亡系数。很容易绘制出下图： 2）SD模型流图 在因果关系环图的基础上可得SD模型流图如图 所示。 4 SD模型的建立 根据上述介绍知识和分析步骤，可得简单的安徽省人口SD模型如下： * POPULAYION SD MODEL OF ANHUI L R R N C C SPEC DT=1/PRINT 1）POP，2）BR，3）DR，PLOT POP，BR，DRPLOT POP 说明： 1）人口数分22个年龄组，即： 1岁，2 4，5 9，10 14，。 。，95 99，100岁数以上计22个年组。 2）增加与迁移人口 在实际应用时要包括：增加人口和迁移人口。 5 模拟或仿真运行 1）进入SD状态 ，即SD专业软件状态，如：PD PLUS，VENSIM等。 2）输入程序 3）运行程序 （1）按初始设计参数，100年总人口的变化情况； （2）调整人口出生系数和人口死亡系数，100年总人口的变化情况； （3）确定各个时期的人口目标，实验取得最佳的BRF，制定合理的人口生育政策。 6 解释分析模拟结果 1）对比分析 按时间段纵向分析，按区域分布横向分析。 2）效果分析 参数的灵敏度，政策的实用性和灵活性。 3）技术分析 SD的理论与方法技术，模型建立技术与技巧，SD专业软件操作技术等。 7 修正再模拟 1）模型结构适合性检验 （1）量纲 量纲必须正确。 （2）极端条件 极端条件下是否仍有意义。 （3）边界条件 2）模型行为适合性检验 （1）结构灵敏度 （2）参数灵敏度 （3）结合灵敏度 3）模型结构与真实系统一致性 4）模型行为与真实系统一致性 现将安徽省人口SD模型列于以下：* POPULATION SD MODEL OF ANHUIC NAGES=22 FOR AGE1=1，NAGES=AG1，AG2，AG3，AG4，AG5，AG6，AG7，AG8AG9，AG10， AG12，AG13，AG14，AG15，AG16，AG17AG18，AG19，AG20，AG21，AG22 FOR AGE2=2，NAGESA L L N POP（AGE1）=IPOP（AGE1）T IPOP(*)=867601,3881562,5772587,7431281,6740139,2488798 4108669，3507741，2741800，2486142，2406997，2031682，1830957 1342043，945653，513026，352871，166180，42022，7165 881,121,29 !29 不详 R ！内积函数，即1620至49岁计7组。T BRF=0.0457795，0.07161650，0.11789292，0.0452570 0.01577322， 0.00636126，0.00148608 R T DRF=0.02216，0.004725，0.001，0.00066，0.00096 0.00156,0.00152，0.00173,0.00213,0.00307,0.00447 0.00741,0.01154,0.01891,0.02899,0.04986,0.07005, 0.0801,0.1751,0.8851,0.9660,0.9752 ！22个年龄组的死亡系数。R ！表示每组每年转出人数。T MT=1，4，5，5，5，5，5，5，5，5，5，5，5，5，5，5，5，55，5，5，1E+30 ！1E+30表示最后一组转出为零。SPEC DT=1/LENGTH=100/PRTPER=2/PLTPER=1/SAVPER=2 PRINT 1）TPOP，2）POP，3）BR，4）DRPLOT TPOP，BR，DRPLOT POP 第4节 系统动力学模型的实习指导4.1 Vensim的基本操作系统动力学模型是由SD专业软件来支撑仿真运行的，例如SD专业软件Micro-Dynamo 、，PD PLUS、Vensim等。但只有Vensim在Windows下运行，而且是非常方便、实用的新版本SD专业软件。故仅介绍Vensim的基本操作。4.1.1 安装与启动(1)安装将装有安装文件和Vensim系统文件一同COPY到硬盘，然后单击安装文件依据提示即可安装。(2)启动安装成功后，快捷方式Vensim图标显示在桌面上，双击该图标即进入Vensim窗口。如图8-2所示。图8-2Vensim主窗口4.1.2窗口组成(1)主菜单File文件管理New Modle建立新模型Open Modle打开模型Save保存模型Print打印模型Check检查模型View观察：对模型的流图进行总体观察Datasets数据集分析：Simulate运行模型Reality Check真实性检验Time Bounds时间控制设置Graph自定义作图：用于由用户选择作图。Build模型建立：实现由其它窗口转换到模型建立窗口Control控制窗口：Variable Selection分析变量选择Time Axis时间轴设置Scalling纵坐标调整Datasets数据集选择Custom Graph自定义作图Output输出窗口：实现从其它跳转到输出窗口(2)主窗口主窗口由三个子窗口组成。Building Window模型建立窗口：用于建立模型Control Window控制窗口：用于调整和控制输出结果，使输出曲线等更加完美Output Window输出窗口：用于结构分析和数据集分析选择三个窗口可通过工具条：BuildingOutputControl实现切换。4.1.3 模型建立模型建立过程是先输入流图，然后输入模型的方程和参数。(1)打开模型建立窗口单击“File”菜单，弹出子菜单；单击“New Modle”或“Open Modle”即可。(2)输入流图绘制流图根据分析的因果图绘制流图如图8-2所示。输入流图采用系统的图示工具条，位于主菜单下的工具行，很容易输入流图。主要工具及使用方法，其功能从左向右依次为：Poiny指针：用