系统动力学第四章精品文稿.ppt
系统动力学第四章1第1页，本讲稿共26页本章主要内容本章主要内容n4.1 4.1 概述概述n4.2 4.2 变量和方程变量和方程n4.3 4.3 延迟函数和平滑函数延迟函数和平滑函数n4.4 4.4 函数函数n4.5 4.5 输出问题输出问题n4.6 4.6 准确度与运行时间单位的选择准确度与运行时间单位的选择n4.7 Rabbit Limits Model4.7 Rabbit Limits Model2第2页，本讲稿共26页4.1 4.1 概述概述n1 1、Vensim PLE Vensim PLE 软件的特点软件的特点 版本：版本：Vensim PLE 5.9e。（1 1）利用图示化编程建立模型）利用图示化编程建立模型 在在Vensim中中，只只要要在在模模型型建建立立窗窗口口画画出出流流图图，再再通通过过Equations输输入入方方程程和和参参数数，就就可可以以直直接接进进行行模模拟拟了了。（2 2）运行于）运行于WindowsWindows下，数据共享性强，提供丰富的输出信息和下，数据共享性强，提供丰富的输出信息和灵活的输出方式灵活的输出方式Vensim的输出信息是非常丰富且兼容性较强，一般的模拟结果，的输出信息是非常丰富且兼容性较强，一般的模拟结果，除了即时显示外，还提供保存文件和除了即时显示外，还提供保存文件和copy到剪切板。到剪切板。3第3页，本讲稿共26页4.1 4.1 概述概述n1 1、Vensim PLE Vensim PLE 软件的特点软件的特点 n（3 3）模型的多种分析方法）模型的多种分析方法 Vensim提供对于模型的结构分析和数据集分析，其中结构分析提供对于模型的结构分析和数据集分析，其中结构分析包括原因树分析、结果树分析和反馈列表；数据集分析包括数包括原因树分析、结果树分析和反馈列表；数据集分析包括数据表、原因图分析、结果图分析、对比图分析等。据表、原因图分析、结果图分析、对比图分析等。（4 4）真实性检验）真实性检验 对于所研究的系统模型中一些重要变量，依据常识和一些基本原则，可对于所研究的系统模型中一些重要变量，依据常识和一些基本原则，可以预先提出对其正确性的基本要求，这些假设是真实性约束，以预先提出对其正确性的基本要求，这些假设是真实性约束，Vensim 提供相应的检验。提供相应的检验。4第4页，本讲稿共26页4.1 4.1 概述概述n2 2、使用、使用VensimVensim软件处理问题的一般过程软件处理问题的一般过程图图4.1 Vensim4.1 Vensim软件处理问题过程软件处理问题过程5第5页，本讲稿共26页4.2 4.2 变量和方程变量和方程n1 1、积累变量方程、积累变量方程 连续方程形式：连续方程形式：离散化方程：离散化方程：LEVEL.K=LEVEL.J+DT*(INFLOW.JK-OUTFLOW.JK)6第6页，本讲稿共26页4.2 4.2 变量和方程变量和方程n2 2、速率方程、速率方程rateS(t)=g lvS(t),aux(t),exo(t),const(t)图图4.2 4.2 人口自然增长模型人口自然增长模型 Births(t)=Pop(t)*BR，Deaths(t)=Pop(t)/AL7第7页，本讲稿共26页4.2 4.2 变量和方程变量和方程n3 3、辅助方程、辅助方程aux(t)=f lv(t),rate(t),aux*(t),exo(t),const图图4.3 4.3 库存模型库存模型TINV=S*INVCov8第8页，本讲稿共26页4.2 4.2 变量和方程变量和方程n4 4、常量方程、常量方程给常量赋值，如人口系统中：给常量赋值，如人口系统中：births=0.020。n5 5、表函数、表函数有些变量之间不能用函数关系式表示，则用图形方式表示出变量之间的有些变量之间不能用函数关系式表示，则用图形方式表示出变量之间的关系，在关系，在vensim称表函数称表函数(lookup function)；Lookup name(Xmin,Ymin)-(Xmax,Ymax),(X1,Y1),(Xn,Yn)。9第9页，本讲稿共26页4.2 4.2 变量和方程变量和方程n5 5、表函数、表函数建立表函数的一般步骤：建立表函数的一般步骤：（1）确定表函数中自变量和因变量；）确定表函数中自变量和因变量；（2）明确变量的取值变化范围；）明确变量的取值变化范围；（3）考虑曲线的形状和斜率，正斜率表示正反馈，负斜率表示负反馈；）考虑曲线的形状和斜率，正斜率表示正反馈，负斜率表示负反馈；（4）选取合适的端点、驻点、拐点；）选取合适的端点、驻点、拐点；（5）作出图形。）作出图形。10第10页，本讲稿共26页4.3 4.3 延迟函数和平滑函数延迟函数和平滑函数 n1 1、延迟函数、延迟函数 图图4.4 4.4 疾病蔓延模型的基本结构疾病蔓延模型的基本结构 图图4.5 4.5 疾病蔓延模型中的疾病蔓延模型中的DELAY1DELAY1环节环节 11第11页，本讲稿共26页4.3 4.3 延迟函数和平滑函数延迟函数和平滑函数 n1 1、延迟函数、延迟函数 “处潜伏期者处潜伏期者”状态变量的方程有：状态变量的方程有：INC(t0)=INF(t0)*TSSSYMP(t)=INC(t)/TSS在在Vensim中只用一个函数中只用一个函数DELAY就可以表示，即就可以表示，即SYMP(t)=DELAY1INF(t),TSS延迟函数的优缺点：延迟函数的优缺点：优点：方便、简洁；优点：方便、简洁；缺点：隐含了部分状态变量，且不能输出它的结果。缺点：隐含了部分状态变量，且不能输出它的结果。12第12页，本讲稿共26页4.3 4.3 延迟函数和平滑函数延迟函数和平滑函数 n1 1、延迟函数、延迟函数根据延迟环节内部包含的根据延迟环节内部包含的状态变量数状态变量数，可分为一阶延迟函数和三阶，可分为一阶延迟函数和三阶延迟函数。延迟函数。（1 1）一阶延迟函数）一阶延迟函数DELAY1(input,delay time)DELAY11(input,delay time,initial value)13第13页，本讲稿共26页4.3 4.3 延迟函数和平滑函数延迟函数和平滑函数 n1 1、延迟函数、延迟函数 n（2 2）三阶延迟）三阶延迟 图图4.6 4.6 三阶延迟疾病蔓延模型基本结构三阶延迟疾病蔓延模型基本结构14第14页，本讲稿共26页4.3 4.3 延迟函数和平滑函数延迟函数和平滑函数 n1 1、延迟函数、延迟函数 n（2 2）三阶延迟）三阶延迟SYMP(t)=INC3(t)/DLINC3(t0)=IV*DLRT2(t)=INC2(t)/DLINC2(t0)=INC3(t0)RT1(t)=INC1(t)/DLINC1(t0)=INC3(t0)DELAY3(INF,TSS)15第15页，本讲稿共26页4.3 4.3 延迟函数和平滑函数延迟函数和平滑函数 n1 1、延迟函数、延迟函数 n（2 2）三阶延迟）三阶延迟DELAY3(input,delay time)DELAY31(input,delay time,initial value)16第16页，本讲稿共26页4.3 4.3 延迟函数和平滑函数延迟函数和平滑函数 n1 1、延迟函数、延迟函数 n物质延迟的阶次物质延迟的阶次阶次是延迟环节内部包含的状态变量数。阶次是延迟环节内部包含的状态变量数。一阶延迟表现出简单的指数增长特性；三阶延迟开始表现出一阶延迟表现出简单的指数增长特性；三阶延迟开始表现出S形增长特性，而三阶以后延迟也是形增长特性，而三阶以后延迟也是S形增长特性，只是程度上形增长特性，只是程度上有所差异。有所差异。实际使用中，具体用一阶延迟函数还是用三阶延迟函数，根据实际使用中，具体用一阶延迟函数还是用三阶延迟函数，根据实际情况分析！实际情况分析！17第17页，本讲稿共26页4.3 4.3 延迟函数和平滑函数延迟函数和平滑函数 n2 2、平滑函数、平滑函数 n信息平滑函数（信息延迟函数）信息平滑函数（信息延迟函数）由于平滑导致信息延迟，所以信息平滑函数又称信息延迟函数。由于平滑导致信息延迟，所以信息平滑函数又称信息延迟函数。(1)一阶信息延迟函数一阶信息延迟函数SMOOTH1(input,delay time,initial value)SMOOTH(input,delay time)(2)三阶信息延迟函数三阶信息延迟函数SMOOTH31(input,delay time,initial value)SMOOTH3(input,delay time)18第18页，本讲稿共26页4.4 4.4 函数函数 n1 1、数学函数、数学函数 INTEG(X）其中其中 ，n为非零整为非零整数；数；19第19页，本讲稿共26页4.4 4.4 函数函数 n1 1、数学函数、数学函数 20第20页，本讲稿共26页4.4 4.4 函数函数 n2 2、逻辑函数、逻辑函数（1）MAX(A,B)（2）MIN(A,B)（3）条件函数）条件函数IF ELSE THENIF ELSE THEN(cond,ontrue,onfalse)21第21页，本讲稿共26页4.4 4.4 函数函数 n3 3、测试函数、测试函数（1）阶跃函数：阶跃函数：STEP（height,step time）（2）斜坡函数：）斜坡函数：RAMP(slope,start time,end time)（3）脉冲函数）脉冲函数单脉冲函数单脉冲函数PULSE(start,width);周期脉冲函数周期脉冲函数 PULSE TRAIN(start,width,tbetween,end)（4）正弦函数：）正弦函数：SIN(X)（5）正态分布函数）正态分布函数RANDOM NORMAL(min,max,mean,stdev,seed)22第22页，本讲稿共26页4.5 4.5 输出问题输出问题 （1）时间序列的数据输出）时间序列的数据输出提供随时间变化的数据列表；提供随时间变化的数据列表；（2）时间序列的图形输出）时间序列的图形输出提供随时间变化的图形输出；提供随时间变化的图形输出；（3）变量之间的关系图）变量之间的关系图速率变量与其对应的积累变量图形；速率变量与其对应的积累变量图形；（4）静态分析结果）静态分析结果如原因树分析、结果树分析、反馈回路分析等。如原因树分析、结果树分析、反馈回路分析等。23第23页，本讲稿共26页4.6 4.6 准确度与运行时间单位的选择准确度与运行时间单位的选择 n时间步长时间步长(time step)(time step)选择选择图图4.7 4.7 茶水温度冷却系统基本结构茶水温度冷却系统基本结构 24第24页，本讲稿共26页4.6 4.6 准确度与运行时间单位的选择准确度与运行时间单位的选择 n时间步长时间步长(time step)(time step)选择选择表表4.1 4.1 茶水温度变化数值表茶水温度变化数值表 25第25页，本讲稿共26页4.6 4.6 准确度与运行时间单位的选择准确度与运行时间单位的选择 n时间步长时间步长(time step)(time step)选择选择从理论上讲，只要时间步长充分小，使离散计算结果与连续解从理论上讲，只要时间步长充分小，使离散计算结果与连续解的差别小到任意期望值，但是步长太小，一方面计算量增大，的差别小到任意期望值，但是步长太小，一方面计算量增大，另一方面当步长小到一定程度时，离散计算结果与连续解之间另一方面当步长小到一定程度时，离散计算结果与连续解之间的差别没有显著减少，所以步长要去一个合理值。的差别没有显著减少，所以步长要去一个合理值。选择选择time step的经验法则：取（的经验法则：取（0.1-0.5）倍的模型最小常数为步）倍的模型最小常数为步长值。长值。26第26页，本讲稿共26页
