系统动力学模型SD3研讨.pptVIP

5. 模型与方程的建立(3) 变量与方程 模型参数 方程的初始值 模型实例 5.1 建立SD方程 建立方程是把模型结构“翻译”成数学方程式的过程，把非正规、概念的构思转换成正式的定量的数学表达式——规范模型。 建立方程的目的：在于使模型能用计算机模拟（或得到解析解），以研究模型假设中隐含的动力学特性，并确定解决问题的方法与对策。 5.1.1 状态变量与Level方程 状态变量是随时间而变化的积累量，是物质、能量与信息的储存环节。 如：人口、企业雇员人数、库存、生产能力、银行存款等。 状态变量的输入、输出变化率使积累量增加或减少。 L LEVEL.K=LEVEL.J+DT * (INFLOW.JK- OUTFLOW.JK) LEVEL——状态变量； INFLOW——输入速率； OUTFLOW——输出速率； DT——计算间隔。 如何确定状态变量 状态变量在回路中的作用 一个因果链中的状态变量具有改变系统整体的动力学性质的能力。 状态变量在回路中的作用 状态变量方程小结 状态变量环节能改变随时间变化的输入量的形状，能削弱输入量与输出量之间的联系，使它们多多少少能独立变化，从而使模型可能具备不平衡的动力学性质。 把某变量定义为状态变量？ 1. 积累量：即首先取决于是否能把这个量看作为某种对时间变化的积累过程。 2. 根据积分过程的时间常数来判断： 时间常数很大时→常数项； 时间常数很小时→辅助变量； 一般情况→状态变量。 5.1.2 速率(变化率)方程 速率方程的功用：把影响系统状态的诸因素——来自系统内外的信息、计划与决策，转化成改变系统状态的行动。 构思模型与建立方程时，一个重要的任务便是寻找适当的方程式去描述速率（或变化率）。 典型的变化率方程（构造复杂速率的基本单元）： LEVEL.K*CONST LEVEL.K/LIFE (GOAL.K-LEVEL.K)/ADJTM LEVEL.K*AUX.K与LEVEL.K/AUX.K EFFECT.K+NORM.K（某些因素的影响作用+额定速率） EFFECT.K*NORM.K（额定速率与某个（或几个）因子的乘积） LEVEL.K*CONST (GOAL.K-LEVEL.K)/ADJTM 速率方程式： R RATE.KL= (GOAL.K-LEVEL.K)/ADJTM 速率方程小结 关于RATE方程的两点说明： 1. 避免违背客观事实 例：世界纯出生率NBR R NBR.KL=(NBRN * P.K) * NBRMM.K*NBRFM.K*NBRPM.K*NBRCM.K NBR——出生率 NBRN——正常的出生系数 P—— 人口 四个对出生率的影响影响因素：NBRMM、 NBRFM、 NBRPM、 NBRCM依次是物质条件、食物供应、污染程度与居住条件。 假如发生世界性灾难？（应分别考虑出生率、死亡率） 2. 方程形式要尽可能自然、合理地描述客观因素的影响 5.1.3 辅助变量与方程 在系统动力学模型中，辅助变量表述了系统内部的信息。建模的任务之一是揭示系统内部的机理和对变量之间的关系加以量化。 辅助方程：概念→量化 辅助变量可为常数项、状态变量、速率或其他辅助变量的任一组合。 辅助变量方程的表达类似于速率方程，都是代数运算，且无标准形式。 辅助变量的表函数表示 构思表函数系一种颇具难度的技巧，表函数以图形表示方法在模型中规定变量之间的关系。 例： 辅助变量的表函数表示 A ELBC.K=TABLE(TELBC,LFO.K,0,1,0.1) T TELBC=1.3/1.28/1.25/1.22/1.18/1.1/1/0.7/0.3/0.1/0 建立表函数的原则 建立表函数时大致要考虑：曲线的斜率和形状，一个或一个以上的特殊点和参考曲线。 设置曲线的斜率。使之与其表示的影响的性质吻合，负值斜率代表负反馈，正值斜率代表正反馈。 选择曲线的形状。小心确定在极端条件下和曲线中部的斜率与曲线的值。 尽可能在表函数上把x,y的特殊点标出，如：x,y分别取0和1时，极端条件下的x,y值和某些研究问题所要求的特殊点。 一些参考线如y=x，y=y0，或y=x0，可用来表示理想与希望的条件，作用达饱和的限制或表函数中的其他内容。 实际系统难以测量时，用表函数表示系统的非线性影响与压力是最有效与方便的。但是切勿滥用表函数，在可以用清晰的表达式表示因素的影响时，不要用表函数。 辅助方程小结 量化模型回路中的信息和变量时，要有清晰的概念。构模者可借助