[理学]系统动力学课程PPT共五章全.ppt

系统动力学：组织学习实验室 课程内容 第一章 系统动力学的孕育、诞生、成长 第二章 系统动力学的问题、思想、方法 第三章 系统动力学建模基础 第四章 DYNAMO：系统思考的计算机实现 第五章 系统动力学模型的经典案例 第六章 系统动力学模型：政策实验室 第七章 系统动力学：组织学习实验室 （实际案例分析） 第八章 思考与展望 第一章 系统动力学的孕育、诞生、成长 1. 定义 系统动力学 (System Dynamics,简称 SD) 是一门分析研究信息反馈(IF) 系统问题的学科 2. 学科性质 SD是一门沟通自然科学与社会科学等的交叉学科 3. 学科地位 纵向：SD是系统科学的一个分支 横向：SD是控制论、信息论、一般系统论加上计算机模拟技术的综合 案例二: 污染净化模型 问题: 分析环境系统中污染净化的动力学机制. (1) 线性污染模型 系统分析: 考虑污染量、污染发生、污染净化与污染净化时间等因素 非线性污染净化模型DYNAMO程序 * POLUTION MODEL L POL.K=POL.J+DT*(POLGR.JK-POLAR.JK) N POL=0 R POLGR.KL=CONST C CONST=3 R POLAR.KL=POL.K/PAT.K A POLR.K=POL.K/POLS C POLS=1 A PAT.K=TABLE(TPAT,POLR.K,O,80,10) T TPAT=0.6/2.5/5/8/11.5/15.5/20/31/50 NAR=POLGR-POLAR 非线性污染净化模型的初值依赖效应 案例三: 库存模型 问题: 分析库存系统进货与发货变化的力学机制 理想的库存系统的两个基本的实际要求： (1)应具有应付各种市场波动的能力； (2)期望库存量应当保持在一个相对稳定的范围内. 系统分析: 考虑库存量、定货率、发货率、期望库存、定货系数等因数。 （1）简单库存模型 因果关系图 系统动力学流图 先不考虑捕捞因素（E=0），系统方程为dx/dt=ax(1-x/L)-bxydy/dt=-cy+dxy可以看出系统中有一个非平凡平衡点 =c/d, =a/b[1-c/(dL)]再考虑捕捞因素(E0)系统方程为 dx/dt=ax(1x/L)-bxy-Ex  dy/dt=-cy+dxy-Ey 此时系统中有一个非平凡平衡点为 (E)=(c+E)/d,  (E)=(a-E)/b{1-a(c+E)/[dL(a-E)]} 模型的动态特性分析不失一般性，取E=0，将系统在其稳定平衡点附近展开并保留其线性部分，略去高次项, 得 式中 方程组的特征方程为 λ2+Eλ+A*D=0 它的两个特征根为 讨论： 1）当IAT小于4倍的WFAT时，系统的两个特征值为一对共轭复数，且它们的实部均小于零，即Reλ10，Reλ20。系统的平衡点为稳定焦点，意味着系统在调节过程中将会出现衰减振荡的特性。 t DINV INV 2) 当IAT小于4倍的WFAT时，系统的两个特征值为一对共轭复数，且它们的实部均小于零，即Reλ10，Reλ20。系统的平衡点为稳定焦点，意味着系统在调节过程中将会出现衰减振荡的特性。 t DINV INV 案例5：捕食与被捕食模型 基本模型 问题：渔业悖论。 伏尔泰拉 （捕食者与被捕食者）模型 模型假设： 1）捕食者是控制被捕食者的唯一因素，被捕食者具有充分的食物，当没有捕食者时，被捕食者以指数规律增长 。 （dx/dt=ax, 当y=0时， a0） 2）被捕食者是捕食者唯一的食物来源，当没有被捕食者时，捕食者以指数规律衰亡。（dy/dt=-cy, 当x=0时) 3) 捕食者与被捕食者共存时发生捕食现象；捕食者因被捕食而减少；被捕食者因捕食而增加。 即有： dx/dt=ax-bxy dy/dt=-cy+dxy a,b,c,d0 捕食者-被捕食者模型因果关系图 捕食者-被捕食者模型系统动力学流图 先不考虑捕捞因素（E=0），系统方程为 dx/dt=ax-bxy dy/dt=-cy+dxy a,b,c,d0 可以看出系统中有一个非平凡平衡点 =c/d， =a/b 再考虑捕捞因素(E0)系统方程为 dx/dt=ax-bxy-Ex=(a-E)x-bxy dy/dt=-cy+dxy-Ey =-(c+E)y+dxy 此时系统中有一个非平凡平衡点