管理系统仿真

《管理系统仿真》由会员分享，可在线阅读，更多相关《管理系统仿真（35页珍藏版）》请在金锄头文库上搜索。

1、管理系统仿真管理系统仿真 内容简介： 系统仿真（亦称系统模拟）是指通过建立和运行系统 的数学模型，来模仿实际系统的运行状态及其随时间变化的 规律，以实现在计算机上进行试验的全过程。这是近30年来 发展起来的一门新兴技术学科。公共管理的对象通常是社会 、经济、军事等复杂系统，一般都不能通过真实的实验来进 行分析、研究。因此，系统模拟技术就成为十分重要甚至必 不可少的工具。本实验在介绍管理系统模拟的概念以及一般 原理、方法和步骤的基础上，主要介绍四种基本的模拟方法 及其模型，即蒙特卡洛模拟方法、排队模型、系统动力学模 拟、多AGENT系统模拟。通过蒙特卡洛模拟可以具体了解管 理系统模拟的基本原理及

2、方法，排队模型与多AGENT系统体 现了离散事件系统模拟的特点与规律，而系统动力学模拟则 是一种可以广泛应用于公共管理决策及政策分析的连续系统 模拟方法。 教学目标通过本课程及实验的学习，使大家了 解系统模拟的概念、原理、方法和步骤，掌 握蒙特卡洛模拟方法、排队系统模拟方法、 系统动力学模拟方法与多AGENT系统模拟方 法，以期为提高公共管理问题的决策水平服 务。 重点难点v系统模拟的概念v蒙特卡洛模拟方法v排队的概念v单服务台与多服务台模型v系统动力学方法v因果反馈结构v多AGENT系统方法v关键词： 系统模拟 连续系统模拟 离散事件系 统模拟 蒙特卡洛模拟 排队论 单服务 台排队 多服务台

3、排队 系统动力学 因 果关系 因果反馈回路 流图 信息流 物 质流 状态变量 流率变量 多AGENT系统参考文献1谭跃进，陈英武等，系统工程原理。长沙 ：国防科技大学出版社，1999 2（加）唐纳德沃特斯著。张志强等译.管 理科学实务教程。北京：华夏出版社，2000 3黄卫伟编著.管理系统模拟的方法及应用.北 京：中国人民大学出版社，1991第一节 管理系统模拟概述 一、管理系统模拟的概念 二、系统模拟的分类 三、系统模拟的基本步骤 问题描述与定义 建立模拟模型 数据采集 模型的确认 模型的编程实现与验证 模拟试验设计 模型的模拟运行 模拟结果的输出与分析 四、管理系统模拟的应用第二节 蒙特卡

4、洛模拟方法 蒙特卡洛模拟是一种特殊而应用广泛 的计算机模拟方法，它是充分利用计算机 计算能力的随机实验方法。蒙特卡罗模拟的 关键是生成优良的随机数。 对应的实验: 实验一:基于高级编程语言及MATLAB的 随机数的产生 实验二:蒙特卡洛模拟算法实现许多程序语言中都自带生成随机数的方法，如 c 中的 random() 函数，Matlab中的rand()函数等。 但这些生成器生成的随机数效果很不一样，比如c 中的函数生成的随机数性质就比较差，如果用 c ，最好自己再编一个程序。Matlab 中的 rand() 函数，经过了很多优化。可以产生性质很好的随 机数，可以直接利用。第三节 连续系统仿真方法

5、用通用的软件开发工具，开发连续系 统仿真的数字计算机仿真模型，理解连续 系统仿真的数值积分算法（龙格库塔法 等） 对应的实验: 实验三:连续系统仿真 第四节 排队模型 一、排队系统基本概念 二、单服务台排队模型 三、多服务台排队模型 四、排队系统模拟 对应的实验: 实验四:离散系统仿真第四节 系统动力学模拟 一、系统动力学方法 二、因果反馈结构 因果关系 因果反馈回路与反馈系统 反馈系统实列 系统动力学流图 三、系统动力学模拟实例 对应的实验: 实验五:系统动力学模拟第五节 基于多Agent系统的模拟 vAgent概念和特性vAgent模型：BDIvAgent结构vAgent通信语言 对应的实

6、验: 实验六:基于多Agent系统的模拟 实验五 基于系统动力学的复杂系统建模 与仿真工商管理学院 工业工程系赖朝安 博士系统仿真概述系统动力学结构模型化原理基本反馈回路的DYNAMO仿真分析系统仿真及SD方法3.1 系统仿真所谓系统仿真，就是根据系统分析的目的 ，在分析系统各要素性质及其相互关系的基 础上，建立能描述系统结构或行为过程的、 且具有一定逻辑关系或数量关系的仿真模型 ，据此进行试验或定量分析，以获得正确决 策所需的各种信息。 一、概念及作用1.基本概念2、系统仿真的实质3.1 系统仿真(1)它是一种对系统问题求数值解的计算技术 。尤其当系统无法通过建立数学模型求解时， 仿真技术能

7、有效地来处理。 (2)仿真是一种人为的试验手段。它和现实系 统实验的差别在于，仿真实验不是依据实际环 境，而是作为实际系统映象的系统模型以及相 应的“人造”环境下进行的。这是仿真的主要功 能。(3)仿真可以比较真实地描述系统的运行、演 变及其发展过程。 3.1 系统仿真 3、系统仿真的作用(1)仿真的过程也是实验的过程，而且还 是系统地收集和积累信息的过程。尤其是对 一些复杂的随机问题，应用仿真技术是提供 所需信息的唯一令人满意的方法。(2)对一些难以建立物理模型和数学模型 的对象系统，可通过仿真模型来顺利地解决 预测、分析和评价等系统问题。 (3)通过系统仿真，可以把一个复杂系统降 阶成若干

8、子系统以便于分析。 (4)通过系统仿真，能启发新的思想或产生 新的策略，还能暴露出原系统中隐藏着的 一些问题，以便及时解决。3.1 系统仿真3.1 系统仿真 二、系统仿真方法系统仿真的基本方法是建立系统的结构 模型和量化分析模型，并将其转换为适合 在计算机上编程的仿真模型，然后对模型 进行仿真实验。由于连续系统和离散(事件)系统的数学 模型有很大差别，所以系统仿真方法基本 上分为两大类，即连续系统仿真方法和离 散系统仿真方法。3.1 系统仿真在以上两类基本方法的基础上，还有 一些用于系统(特别是社会经济和管理系统) 仿真的特殊而有效的方法，如系统动力学 方法、蒙特卡洛法等。系统动力学方法通过建

9、立系统动力学模 型(流图等)、利用DYNAMO仿真语言在计 算机上实现对真实系统的仿真实验，从而 研究系统结构、功能和行为之间的动态关 系。3.2 系统动力学原理1、由来与发展 Systems Dynamics, SD/ J.W. Forrester(MIT)Industridl Dynamics (ID), 1959Principles of Systems, 1968Urban Dynamics (UD), 1969World Dynamics (WD), 1971SD, 19723.2 系统动力学原理2、研究对象及其结构特点 （1）研究对象社会系统（2）结构特点抉择性具有决策环节（人、信

10、息）自律性具有反馈环节非线性具有延迟环节（3）SD将社会系统当作非线性(多重)信息反 馈系统来研究3.2 系统动力学原理认识 问题界定 系统要素及其因 果关系分析建立结 构模型建立数 学模型仿真 分析比较与 评价政策 分析（流图）（DYNAMOY方程） 3、工作程序3.2 系统动力学原理4、系统动力学模型（1）常用要素流速率水平变量源与汇参数（2）流图符号 流 实物流信息流 速率变量 水准变量 L1 辅助变量 A1 。R1 R1( )。3.2 系统动力学原理 明确问题及其构成要素； 绘制要素间相互作用关系的因果关系 图。注意一定要形成回路； 确定变量类型（L变量、R变量和A变量 ）。将要素转化

11、为变量，是建模的关键一步。 在此，应考虑以下几个具体原则：（3）流图绘制程序和方法3.2 系统动力学原理a. 水准（L）变量是积累变量，可定义在任 何时点；而速率（R)变量只在一个时段才有 意义。b. 决策者最为关注和需要输出的要素一般 被处理成L变量。 c. 在反馈控制回路中，两个L变量或两个 R变量不能直接相连 。d. 为降低系统的阶次，应尽可能减少回路 中L变量的个数。故在实际系统描述中，辅 助（A）变量在数量上一般是较多的。 绘制SD流图。3.2 系统动力学原理5、举例L1R1（利息1）C1(利率)IR1(订货量)库存量DY(期望库存)（库存差额）PR1R2(出生人口)(人口总量) (

12、死亡人口)C1（出生率）C2（死亡率）组织改善组织 绩 效组织 缺陷。3.2 系统动力学原理1、基本DYNAMO方程( DYNAmic Model)水准方程（L方程） L L1K=L1J+DT*(RIJK- ROJK)速率方程（R方程） R R1KL=f ( L1K,A1K,)辅助方程（A方程） A A1K=g(L1K,A2K, R1JK, )赋初值方程（N方程） N L1=数值 或 L1=L10 L10=数值常量方程 （C方程） C C1=数值3.3 基本反馈回路的DYNAMO仿真分析2、一阶正反馈回路 年人口 增 加人 口 数（+）P +PRPRPC1（人口年自然增长率0.02）。L PK

13、=PJ+DT*PRJKN P=100R PRKL=C1*RKC C1=0.02PPR 01002 11022.04 2104.042.0808 p100 0 一阶正反馈（简单 人口问题）系统输 出特性曲线3.3 基本反馈回路的DYNAMO仿真分析3、一级负反馈回路 库存量库存 差额订货 量+()R1DI +期望库存Y1000Z （订货调整时间，5）IR1DY（6000）。 。 。L IK=IJ+DT*R1JKN I=1000R R1KL=DK/ZA DK=Y-IKC Z=5C Y=6000IDR10100050001000120004000800228003200640It10000一阶负反馈

14、（简单库存控制）系统输出特性曲线4、简单库存控制系统的扩展 库存量入库量途中存货量订货量库存 差额I （）+（）GR2+R1+Z（5）Y（6000）I G100001000R1 R2DW（10）。L GK=GJ+DT*(R1KL- R2JK)L IK=IJ+DTR2JKR R1KL=D/ZA D=Y-IKC Y=6000C W=10,Z=5C I=1000C G=1000060001000t二阶负反馈系统输出特性曲线I交 （到 ）货 率库存量测试函数 正常销售 （发货）率。销售（发货）率T3MAXT2库存差额 期望 库存SMOOTH平均销 售（发 货）率订货率Step Ramp Pulse Sin NoiseDELAY。 。