系统动力学与动态系统描述-方程.docVIP

系统动力学与动态系统描述 李旭 教授 复旦大学管理学院 从库存系统开始认识SD方程 ? 对右图的库存系统考虑： – 库存是如何变化的？ – 如何进行订货决策？ – 如何用数学方法描述？ ? 库存变化规律： ? 销售量的描述： ? 决策过程描述： ? 辅助计算描述： 方程及其理解 ? SD方程的概念： – SD方程是在流图基础上对系统要素之间的关系定量描述的一组数学 关系式； – SD方程是从一组已知的初始状态开始确定下一组状态的递推关系式； – SD方程中要有一个恰当的时间间隔，以完成方程的递推； ? SD方程的理解： – SD方程的实质是微分方程组，由于规模和非线性等原因不能求得解 析解，所以只能求其数值解。即差分化处理后仿真； – 按照上述规则递推就可以得到各个变量随时间变化的曲线。即系统的 变化过程。 SD方程的种类 ? 水平方程 (L) ? 速率方程 (R) ? 辅助方程 (A) ? 常量方程 (C) ? 初值方程 (N) SD方程中的时间描述 ? 为了完成递推计算，需要首先明确三个基本时间 参数：时点、区间、差分步长。 ? 时间参数的描述： – K：现在时刻； – J：前一个时刻； – L：下一个时刻； – JK：时刻J和K之间的区间； – KL：时刻K和L之间的区间； – DT：差分步长。 水平方程（L） ? 反映系统状态随时间的变化，是变化对时间的积 累。因此具有固定的形式： ? SD中采用差分方程的形式： L 方程的理解 ? 对水平方程的理解： – 水平方程是一个一阶差分方程，具有固定的表现形式； – 水平方程是一个有记忆的量，方程中一定有其前一时刻 的状态值； – 水平方程是将决策变成行动，即将速率变量转换成水平 量的方程，因此方程中一定含有速率量； – 水平方程是变化对时间的积累，因此方程中一定含有 DT，并且 DT 只能出现在水平方程中。 速率方程（R） 观测状态 偏差 行动→R 目标状态 ? 方程原理： – 系统变化的自然规律。例如，人口的死亡。 – 人们控制系统的主观愿望。例如，订货决策。 ? 一般形式： R 方程的理解 ? 速率方程的实质是自然规律或决策策略，由这些 规律或决策策略改变系统的状态； ? 速率方程最终是水平变量和常量的函数，但为了 更好地描述决策过程或表达清楚，速率方程中经 常包括辅助变量； ? 速率方程中不出现具有积分意义的差分步长DT。 常见 R 方程的形式 ? R = Level×Const – 例如： ? R = Level/Life – 例如： 常见 R 方程的形式 ? R = (Goal-Level)/Adjustment – 例如： ? R = Normal×Effect – 例如： EIR：监禁率； PC：监狱容量； ASL：平均服刑期 CR：容量比； ECTS：容量比对服刑期的影响 ATS：实际服刑时间 常见 R 方程的形式 ? R = Normal + Effect – 例如： 被捕食者 死亡数量 被捕食者自 ＝ ＋ 然死亡数量 捕食者吃 掉的数量 辅助方程（A） ? 从理论上看： – 水平方程 完全确定和计 – 速率方程 算系统的状态 – 常量、初始条件 ? 问题： – 计算式过于冗长、复杂； – 不便于描述自然规律或决策过程； – 不便于利用中间结果分析问题。 辅助方程（A） ? 辅助方程及其作用： – A 方程是计算 R 方程的子方程，即用于计算辅助变量的取值； – A 方程描述自然规律或决策过程； – A 方程一般具有实际意义。 ? 辅助方程的建立： – 建立辅助方程一般采用“跟踪”法，按逻辑顺序计算； – 变量之间的运算规则，可以根据实际意义确定； – 量纲分析是建立辅助方程的重要技巧； – 辅助方程之间不能出现“环”。 辅助方程之间的“环” ? 环：辅助变量“闭合”引用： B A C D ? 用水平变量解开“环”： B A C D 人口问题(1) 人口问题(1) 人口问题(2) 人口问题(2) 库存问题(1) 库存问题(1) 库存问题(2) 库存问题(2) 传染病问题(1) 传染病问题(1) 传染病问题(2) 捕食者与被捕食者问题 ? X：被捕食者数量； ? BX：被捕食者出生数量； ? DX：被捕食者死亡数量； ? Y：捕食者数量； ? BY：捕食者出生数量； ? DY：捕食者死亡数量； ? BXY：被捕食者数量对捕食者成活率的影响； ? DYX：单位捕食者、单位时间内，作为食物，对被捕食者 的需求。 捕食者与被捕食者问题 捕食者与被捕食者问题 方程的计算 J K L t J K L ? 计算顺序如下： 时间向前推进一个DT 计算K时刻 计算K时刻 计算KL区间 的水平量 的辅助变量 的速率量 方程计算举例 ? 0时刻(取 DT＝1)： – L