第四章_系统仿真及系统动力学方法

1、山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院认识认识问题问题界定界定

2、系统系统要素及其因要素及其因果关系分析果关系分析建立结建立结构模型构模型建立数建立数学模型学模型仿真仿真分析分析比较与比较与评价评价政策政策分析分析（流图）（流图）（DYNAMOYDYNAMOY方程）方程） 山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院 流流 实物流实物流信息流信息流 速率变量速率变量 水准变量水准变量 L1 辅助变量辅助变量 A1 。R1R1( )。山东建筑大学管理工程学院山东建筑大

3、学管理工程学院 参数（量）参数（量） （常量）（常量）（初值）（初值） 源与洞源与洞 信息的取出信息的取出 L1。A1。山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院 在绘制流图时，应特别注意形成正确的回路和用好信息在绘制流图时，应特别注意形成正确的回路和用好信息连接线，并注意不要把不同的实物流直连在一起连接线，并注意不要把不同的实物流直连在一起( (参见下例参见下例) )。XFS（分到新房户数）分到新房户数）YMS（已买家俱新房户数）已买家俱新房户数）FFL(FFL(分房数量分房数量) )XSL(家俱数量家俱数量)（未分到新房户数未分到新房户

4、数 ）（未买家俱户数） 错误 山东建筑大学管理工程学院山东建筑大学管理工程学院XFSYMS正确 山东建筑大学管理工程学院山东建筑大学管理工程学院 L1R1（利息1）C1(利率)IR1(订货量)库存量DY(期望库存)（库存差额）PR1R2(出生人口)(人口总量) (死亡人口)C1（出生率）C2（死亡率）组织改善组 织绩 效组织缺陷。 山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院山东建筑大学管理工程学院三、基本反馈回路的三、基本反馈回路的DYNAMO仿真分析仿真分析1 1、基本、基本DYNAMODYNAMO方程方程 DYNAmic MOdel水准方程（水准方程（L L方

5、程）方程） L L1K=L1J+DTL L1K=L1J+DT* *(RIJK-(RIJK-ROJK)ROJK)速率方程（速率方程（R R方程）方程） R R1KL=f ( L1K,A1K,)辅助方程（辅助方程（A方程）方程） A A1K=g(L1K,A2K,R1JK,)赋初值方程（赋初值方程（N方程）方程） N L1=数值数值 或或 L1=L10 L10=数值数值常量方程常量方程 （C C方程）方程） C C1=数值数值 山东建筑大学管理工程学院山东建筑大学管理工程学院2 2、一阶正反馈回路、一阶正反馈回路 年人口增 加人 口数（+）P +PR PR PC1（人口年自然增长率人口年自然增长率0

6、.02）。L P K=PJ+DTJ+DT* *PRJKPRJKN P=100R PRKL=C1KL=C1* *RKRKC C1=0.02 PPR0100211022.042104.042.0808p1000一阶正反馈（简单一阶正反馈（简单人口问题）系统输人口问题）系统输出特性曲线出特性曲线山东建筑大学管理工程学院山东建筑大学管理工程学院3、一级负反馈回路一级负反馈回路 库存量库存差额订货量+()R1DI+期望库存Y1000 Z（订货调整时间，5）I R1D Y（6000）。L IK=IJ+DTK=IJ+DT* *R1JKR1JKN I=1000R R1KL=DK/ZKL=DK/ZA DK=Y-

7、IKK=Y-IKC Z=5C Y=6000 IDR10100050001000120004000800228003200640It10000 一阶负反馈（简单一阶负反馈（简单 库存控制）系统输库存控制）系统输 出特性曲线出特性曲线山东建筑大学管理工程学院山东建筑大学管理工程学院4、简单库存控制系统的扩展简单库存控制系统的扩展 库存量入库量途中存货量订货量库存差额I（）+（）GR2+R1+Z（5）Y（6000） I G100001000R1 R2DW（10）。山东建筑大学管理工程学院山东建筑大学管理工程学院L GK=GJ+DTK=GJ+DT* *(R1KL-R2JK)(R1KL-R2JK)L IK=IJ+DTR2JKK=IJ+DTR2JKR R1KL=D/ZKL=D/ZA D=Y-IKKC Y=6000C W=10,Z=5C I=1000C G=1000060001000t二阶负反馈系统输出特性曲线I山东建筑大学管理工程学