第3讲 连续模拟原理.pptVIP

第3讲 连续模拟原理 第3讲 连续模拟原理 3.1 连续模拟基础 3.2 系统动力学 3.3 连续系统模型的建模示例 3.1 连续模拟基础 3.1.1 术语 1）状态变量 两类状态变量：水平变量（Level Variables）和速率变量（Rates Variables）。 2）模拟时间 3）状态方程 3.1.2 连续模拟模型 1）微分方程模型 a. 常微分方程：称为集中参数系统 b. 偏微分方程：称为分布参数系统，在管理系统中，主要表现为具有多个自变量的系统。 2）差分方程模型 3.1.3 连续模拟的过程 1）系统结构的建立 就是系统动力学中的流图绘制过程。 2）系统模型的建立 即在流图的基础上建立一组反映系统行为特征的数学模型，它描述了各项状态变量与主要自变量－模拟时间的关系。 3）离散化转换 包括单步法、多步法和预测校正方法，单步法包括欧拉法和龙格-库塔法。 其中，称为步长（时间增量），fn为y(t)在时刻时的导数值。因已知初始值，由欧拉公式可以推出，然后求出，依次类推。 欧拉法（Eular） 欧拉法又称折线法，其近似计算原理如图1.3.2所示。 欧拉几何图形表示 欧拉递推公式： 龙格-库塔法（Runge-Kutta） 德国Runge和Kutta两人先后提出：用函数值f的线性组合来代替f的高阶导数项，则可避免计算高阶导数，又可提高数值积分的精度。其方法如下。 下一时刻的变化量并不是取前一时刻的变化率与步长之积，而是取及两时刻的斜率平均值与步长相乘，即： Yn+1=Yn+T(fn+fn+1)/2 所以精度要比欧拉法高一些。 二阶龙格-库塔法几何图形表示 4）模拟模型的建立 即将数学模型转换为能在计算机上运行的计算机程序。一般采用专用的模拟语言：Anylogic 、Vensim 。 5）模拟实验与分析 模拟结果的分析形式： a. 以时间为横轴，状态变量为纵轴。 b. 以状态变量为横、纵轴，分析系统状态变化的相轨线。 3.2 系统动力学 3.2.1 基本原理 1）物流与信息流 2）延迟 3）因果关系 4）反馈环 5）结构 3.2.2 系统动力学图解 （1）流位（LEVEL） （2）流率（RATE） （3）流（Flow） （4）辅助变量（AUXILIARY） （5）常量（CONSTANT） （6）信息的取出 （7）图外的变量 （8）函数 （9）外生变量 （10）延迟环节（DELAY） （11）源（Source）和漏（Sink） 3.2.3 系统结构的建立 1）划定系统边界 划出系统边界的一条准则是把系统中的反馈回路考虑成闭合的回路。应力图把哪些与建模目的关系密切、重要的量都划入边界。 2）找出系统的因果关系环 在工程系统中讨论的几乎都是负反馈环，但管理系统中除了负反馈环之外，还存在着指数增长过程的正反馈环，例如，人口和物质财富等的增长。 变量之间的复杂关系常常导致反馈环极性的变化。复杂系统中不仅有正反馈环，而且还有负反馈环，常常会出现反馈极性的转移。 通常一个复杂系统的因果关系链可如下图所示。 反馈极性转移 3）找出反馈环中的流位与流率 系统结构的每一个反馈环中至少包含两种基本的变量，即流位与流率。 在任一反馈环中流位变量与流率变量总是同时出现，交错出现的。 4）决定流率的结构 流率方程表达了系统状态的控制策略，描述了系统流位动态变化的内在规律。要很好的决定流率的结构必须摸清系统运行的机制。管理系统一般都受到人为的控制，哪些流率适合人为控制，需要去思索去研究。 例如，在人口增长系统中，控制人口总数量，要控制人口出生率而不是死亡率；在商店系统中，控制商品的存货量，要控制商品的进货率而不是销售率。 3.3 连续系统模型的建模示例 3.3.1 基于微分方程组的建模 现以一个湖的生态系统为例，说明基于微分方程组的系统模型的建模。 该系统模型包含有以下状态变量：（1）植物生成数量XP；（2）吞食植物的虫类生成数量XH；（3）食虫植物的生成数量XC；（4）太阳能量供应量XS；（5）湖底的有机物沉淀量XO；（6）扩散于周围环境的能量XE。 代表XP的增加率，代表XH的增加率，代表XC的增加率。 太阳能量供应量XS、湖底的有机物沉淀量XO和扩散于周围环境的能量XE为辅助变量。可得该系