第三章DYNAMOL模拟语言分解.ppt

第三章 DYNAMO模拟语言 DYNAMO，取名来自Dynamic Models的缩写，是一种计算机模拟语言系统。 用DYNAMO写成的反馈系统模型经计算机进行模拟，可得到随时间连续变化的系统图像。建模者可据此分析系统的结构、功能及行为。在此基础上，进行政策分析与决策，实现预测和择优的目的。 1 DYNAMO中的时间下标 例如： 库存系统的方程可用DYNAMO表示为： INV.K=INV.J+DT*(ORRE.JK-SH.JK) 式中： INV.K——库存现有量； INV.J——DT前的库存量（J时刻的库存量）； ORRE——JK间隔内的入库量； SH——KL间隔内的发货量。 假如某仓库每月入库量为80件，发货量为100件， 5个月前的库存量为1200件，求现在的库存量？ INV（现在）＝INV（过去）＋（时间间隔）×（纯速率） ＝1200＋5×（80－100） ＝1200＋5×（－20） ＝1200－100 ＝1100件 2 DYNAMO的有关规定与规则 2．1 变量名字符的规定 变量名字符不超过6个。变量名的第一字符必须是字母，其后则可为字母或数字。 如：LEV、LEV1、LEV2、X1、X2等 2. 2 代数运算符的表示 SD通用的代数运算符： 加法：= 减法：－ 乘法：* 除法：/ 代数运算的顺序是：先乘方、开方，再乘、除，最后加、减。 括号内的加、减优于括号外的运算。 同一层次的运算符，按先左后右的原则进行运算。 2．3 方程的列数 一个方程中不能有空格，每行不能超过72列。一行不够可另起一行，但第一列必须以字符“X”作为标记。（P66） 2 ． 4 变量与常量 DYNAMO中的所以数量可分为大两类： （1）变量：其值是可变的，有时间下标。 （2）常量：其值在一次模拟过程中不变，没有时间下标。 例如：SALE.K ——变化的销售量 SALE——固定的销售量 3 DYNAMO的方程式 3.1 状态（State，Level）变量方程 在SD中，凡是能对输入和输出变量进行积累的变量称为状态变量，计算状态变量的方程称为状态变量方程。 状态变量方程在 DYNAMO模型中，以字母L为标志写在第一列。 例如1： L LEVEL.K=LEVEL.J+DT*(INFLOW.JK-OUTFLOW.JK) 式中： LEVEL——状态变量； INFLOW——输入速率（变化率）； OUTFLOW——输出速率（变化率）； DT——计算间隔（从J时刻到K时刻）。 例如2： L POP.K=POP.J+DT*(BIRHS.JK-DEATHS.JK) 式中： POP——人口（人）； BIRTHS——出生率（人/年）； DEATHS——死亡率（人/年）。 3．2 速率（Rate）方程 在状态变量方程中代表输入与输出的变量称为速率，它由速率方程求出。在DYNAMO中，速率方程以字母R为标志。与状态变量方程不同的是，速率方程没有标准格式，只能根据系统的具体情况书写。 例如1： R BIRTHS.KL=BRF*POP.K 式中： BIRTHS——出生率（人/年）； BRF——出生率系数（1/年）； POP——人口（人）。 例如： R DEATHS.KL=POP.K/AVLIFE 式中： DEATHS——死亡率（人/年）； POP——人口（人）； AVLIFE——平均寿命（年）。 速率方程的特性： （1）速率方程无一定的格式； （2）速率的值在DT时间内式不变的。速率方程式在K时刻进行计算，而在自K至L的时间间隔（DT）中保持不变； （3）速率的时间下标为KL。 3．3 辅助（Auxiliary）方程 在建立速率方程之前，若未先做好某些代数计算，把速率方程中必需的信息仔细加以考虑，那么将遇到很大的困难。这些附加的代数计算，在DYNAMO中称为辅助方程，方程中的变量则称为辅助变量。 在DYNAMO语言中，辅助方程要以字母A为标志写在方程的前面。 辅助方程的特性： （1）辅助方程没有统一的标准格式； （2）辅助变量的时间下标为K； （3）辅助变量可由现在时刻的其他变量：状态变量、速率变量或其他辅助变量求出。 3．4 表函数（Table Function）