软件包中系统动力学函数.docx
《软件包中系统动力学函数.docx》由会员分享，可在线阅读，更多相关《软件包中系统动力学函数.docx（16页珍藏版）》请在冰点文库上搜索。
软件包中系统动力学函数
软件包中系统动力学函数
————————————————————————————————作者：
———————————————————————————————— 日期:
第８章Venｓim ＰLE软件包中系统动力学函数
系统动力学所以能处理复杂的系统问题，除提出流位流率系简化流率基本入树建模法去描述系统外，还有一个重要原因是其专用软件都设计了一系列通用的系统动力学函数。
第一节数学、逻辑、测试函数
§ 8.1.１数学函数
ＶenｓiｍPLE备有五种普通数学函数供用户使用。
１．ＳIN（X）
定义1：
SIN（X）为三角正弦函数,X须以弧度表示,其值小于8.35×１05
当自变量是角度时,应通过乘以2π/３6０转化为弧度。
2.EXP（X）
定义2:
EXＰ（X）=eX,ｅ是自然对数的底,e=２．7182…，Ｘ的值必须小于３６。
人们常用指数函数去描述系统,有了上面函数将会带来很大方便。
3.LN（X）,变量X大于零。
即以e为底的对数函数，它与ＥXP（X）互为反函数，这样可以用EXP（X）和LN（X）来计算非以e为底的幂函数和对数函数。
4．SQRＴ（X）=√X—，Ｘ必须是非负量。
5.AＢS（Ｘ）=│X│，对Ｘ取绝对值。
§ ８.1.２逻辑函数
逻辑函数的作用类似于其它计算机语言中的条件语句,VensimPLE的逻辑函数有三种。
1.最大函数MＡX（P，Ｑ）
ＭAX表示从两个量中选取较大者，P和Q是被比较的两个量,结果也是在这两个量中选取。
P 若P≥Q
定义１：
若ＭAX（P，Q）＝
Q若P≤Q
其中Ｐ,Q是变量或常量,则MAX（P,Q）为最大函数。
可用MAX函数从多个量中选取较大者。
如从P,Q，D三个量中选择较大者可用：
MAX（D,ＭAX（P,Ｑ））。
最小函数 Q 若P≥Q
定义2：
若MIＮ（P，Ｑ）=
P 若P≤Q
则MＩN（P，Q）为最小函数。
1．MIN同MAX一样,可以从MIN（P,Q）
基本功能中派生出各种用法。
3.选择函数ＩFＴHENELSE（C,Ｔ,F）
定义3：
若IFTHENＥLＳE（C,T,F）
ＴC条件为真时
= （C为逻辑表达式）
Ｆ 否则
则IＦTHENＥＬSE（C，T,F）为选择函数。
IFTHＥNELSE函数常用于仿真过程中作政策切换或变量选择。
有时也叫条件函数。
§8.１.３ 测试函数
设计这一部分函数的目的主要是用于测试系统动力学模型性能用,所以称为测试函数。
在给出测试函数以前,我们必须重申一个概念,系统动力学的变量皆是时间TＩME的函数,所以当仿真时间ＴＩＭE发生变化时,各变量值都随之发生变化。
不过，各变量与TＩME的依赖关系存在差别,有的是以ＴIME为直接自变量,有的则是间接变量。
测试函数以TIMＥ为直接自变量,但在函数符号中常缺省。
1.阶跃函数STEP（P，Ｑ）
定义１:
Q若TＩME≤Q
SＴEＰ（P，Q）=
P 若ＴIME>Q
其中，Ｐ---阶跃幅度;Q--－STEP从零值阶跃变化到P值的时间，则STEP（P,Q）为阶跃函数。
Ｐ
SＴEP
A*（C－Ｂ）
A
*（TIME-Ｂ）
跃Q B TIME
图8.1.１
阶路STＥP（P，Q）和斜坡RＡMP（A,B，C）函数图
2.斜坡函数ＲAMP（Ｐ,Q，R）
定义2:
0若TIME≤Q
RAMＰ（Ｐ，Ｑ,R）= Ｐ *（TIME-Q）,若R≥ＴIME>Ｑ
P＊（R-Q）,若ＴIME≥R
其中,P为斜坡斜率，Q为斜坡起始时间，Ｒ为斜坡结束时间,则RAMP（P,Q,R）为斜坡函数。
3.脉冲函数PULSE（Q,Ｒ）
定义３:
若PＵＬSE（Q,R）随TIME变化产生脉冲。
其中:
Ｑ--－第一个脉冲出现的时间
Ｒ－--相邻两个脉冲的时间间隔
脉冲宽度为仿真步长，则PULSE（Ｑ,R）为脉冲函数。
4.均匀分布随机函数RANＤOMUＮIFORM（A，B，S）
定义5：
RＡNDＯMUNＩＦOＲM（A,B,Ｓ）产生在区间（A，B）内的均匀分布随机数，S给定随机数序列就确定,S取不同的值产生随机数序列也不同。
RAＮDOMUNＩFORM（A，B,S）为均匀分布随机函数。
上面我们给出了四种测试函数,实际上还有前面数学函数SIN（X）等也可以作为测试函数。
一个系统动力学模型,可以通过改变常数再运行的办法,实现多种测试函数分别进行测试。
第二节表函数
§ 8.2．1表函数表示形式
1.ＶenｓｉmPLE中表函数表示形式
定义1自变量与因变量的关系通过列表给出的函数叫表函数。
例如下表就确定了一个表函数。
自变量Ｘ
0
1
1.５
2
2.5
因变量Y
0.5
11
2
５
１0
表函数是系统动力学的一个重要特征,它用于建立两个变量之间的非线性关系,特别是软变量之间的关系。
例如：
员工士气对工作效率的影响程度。
一般,两个变量先归一化或者先规整化,再根据经验给出大致的关系图来。
这样设计的变量是无量纲量。
在进入VenｓimPLE软件ＥｑuaｔｉonEditor,即点去图标Y＝X２ 后，若方程还未定义,有ASGrapｈ选项。
选择此选项，会出现下面对话框。
该对话框用于图形化定义，上例表函数可直接填入框中（图8.2.１）。
包括自变量和函数值即因变量值列举,自变量和函数的最大值等。
当自变量为非已知统计点时,可用线性插值法取其近似值。
用鼠标左键在图形框中点按,会自动构成图形。
VensimPＬＥ软件中表函数表达形式还可通过选择方程类型TYＰE中Ｌｏoｋｕp进行列举表出,即把表函数自变量，因变量最大值、最小值及一些自变量与因变量对应的点值列出。
如上例描述的表函数可以在方程输入框写成:
[（0,0）-（1０,１0）]（０,0.5）（1,1）（１.5,2）（2，5）（２.5,10）其中[]中前面（）中0,０分另为自变量、因变量最小值,若自变量小于最小值，因变量取最小值,后面（ ）中１0，1０分别为自变量、因变量最大值，若自变量超出最大值,因变量取最大值,[]后面五个（ ）是已知自变量和因变量对值点,若自变量值不在给出点中，则自动用线性插值法求因变量对应值。
自变
量、因变量的最小值、最大值可依据实际背景来确定，列出的对应值点作为已知点可从历史数据中计算或分析给出。
VeｎｓiｍPＬE专用软件对表函数的增减性、取值间隔均匀性没有严格要求，但使用者可根据实际问题给出取值间隔、分段满足增减性的表函数。
在Ｖensim ＰLE中建立的入树或流图内一个表函数必须有三部分完成,即一个自变量X，一个因变量Y及一个Y关于Ｘ的因子表,其因果关系为:
Ｙ变量 方程可写为:
↗ ↖ Y变量=Ｘ因子表（X变量）
X因子表 X变量
上例表函数中
X因子表=［（０，0）－（1０，10）]（０，0．5）（１,1）（1．５,２）（2,5）（2.5,１0）
表函数的建立方法将在§８.２.2介绍。
2.Mｉcｒｏ DYNAMＯ及ＰD PLUS中表函数表示形式与ＶｅnｓimPLＥ软件不同的是在Miｃｒo ＤYNＡＭＯ及PDPＬUS中有特定不同类型，其表示含义可由定义给出并固定下来。
现使用Venｓim 软件的读者，可以不阅读下面内容。
⑴MicroDＹNAMO两类表函数
定义1：
若TABLＥ（ＴY,Ｘ.K，ＸLOW，XＨIＧH，XＩＮＣＲ）
中:
TＹ---表量名（因变量已给值）
Ｘ-－-自变量
XLOW-－-自变量Ｘ的最小值
XＨIGH---自变量X的最大值
XＩNCR-－-自变量X的取值间隔
自变量取值为XLOW至XＨＩGH间以等间隔ＸINCＲ取X１,X2,…… Xmm个值，且
ｍ ＝（XHＩＧH－ＸLＯW）/XINＣＲ＋1
对应于X1，Ｘ2，……Xm 的TY的值在DYNＡMＯ方程中以T方程:
ＴTＹ＝ E1／E２/……／Em给出。
当X0∈（XLOＷ,XHIＧH），但Ｘ0≠Xi（ｉ＝1,2,…,ｍ）时，其变量值按线性插法给出，当Ｘ的值超出[XLOW,XHＩGH］范围时,因变量取对应的端点值，并给出警告信息。
则 TABLＥ（TＹ，X．K， ＸＬOＷ，XHIGＨ,XINＣR）称为第一类表函数。
T TＹ=E1／E2/……/Em称为其表量语句，又称为T语言。
例1：
已知两变量X和Y，因变量Y随自变量X变化的关系的曲线所示（图８．２.2）
自变量Ｘ从Ｘ=-3开始，按等距离取7个点得表８.2．1。
设Y为辅助变量，用第一类表函数语句表示的DYNAMO语句为:
ＡY.K=TＡBLE（ＴＹ，X.K.－3.3.1）
T TＹ＝-2０／0/10/１6/20/24／30。
X
-3 －２-1 0 １ 23
Ｙ
-20 0 １0 16 2０ 24 30
表8．2.１
注１：
该例在VeｎｓiｍＰLE中变量关系图为：
Ｙ变量 方程可写为:
↗ ↖ Y变量=Ｘ因子表（X变量）
X因子表 X变量
X因子表 ＝[（-3，-2０）-（3，20）］（-3,-２0）,（-2,０）,（－１,10）,（０,16）,（1,20），（2,24）,（３,３0）（曲线如图８.2．3）
当X值超出[-3,3]时,Y取对应的端点值，不给出警告信息。
定义2:
若TABHL（ＴY，X．K,XLＯW,XHIGＨ,XＩNＣR）中随X的取值范围超出[XLOW,XＨＩGH]时,因变量取对应的端点值,但不给出错误信息外,其它内容与第一类表函数相同,则TAＢＨL（TY,X.K，XLOW，XＨＩGH，XINCR）称为第二类表函数。
注２：
在PＤPLUS语言中，表函数中Ｔ语句中斜线改为了逗号，且定义了下面两类表函数,不妨称为第三、第四类表函数。
第三类表函数为：
TAＢXL（TＹ,Ｘ.Ｋ,XＬOＷ，ＸHIGH，XIＮＣＲ）,此表函数,除Ｘ取值范围超出[ＸLOW，XHIGH］范围时，因变量取端点的趋势外推值外，其它内含与第一类表函数相同。
第四类表函数为：
ＴABＰL（TY，X.K，ＸLOW,XHIGH，ＸＩＮCＲ），此表函数除利用多项式使曲线在各点起滑连接代替线段连接外,其它含义与第一类表函数相同,但表量语的数值后要加上m个零,如,对前面第一类表函数的例子，改为TABPＬ,则写为：
A Ｙ.K＝TABＰL（TＹ,X.K．-3．3.1）
TTＹ=-２0，0,１0,16,２０,24，３０,0，０,0，0,０,0，０。
§8.2.2表函数建立方法介绍
前面介绍的VensimＰLE专用软件,DYNAMO语言表达表函数的方式，是比较简单、容易接受的，但一个表函数的建立却不是如此容易的,往往是一个定性与定量相结合反复分析的结果,要建立一个具体表函数，肯定必须考虑所涉及的自变量、因变量的实际背景,再仔细研究其包含的一般数学问题及一般统计问题,进行深层次的量化分析,最后得出能反应变量间一般关系规律的量表作为表函数才能用于ＳD模型，日常生活中，可以经常碰到时间间隔相同的统计年报表、季报表、月报表,这些都是表函数。
但是,在建立一个实际系统的ＳD模型时，这些统计报表很难作为一个完整的表函数直接放入模型中。
其一，表函数不一定以时间为自变量,实践表明，表函数的自变量很多是同模型中其它一个或若干个变量的因变量；其二,统计报表没有未来若干年的预测数据,而SD模型的目标重点在对未来规律的仿真。
根据以上分析和以往的经验,建立一个具体表函数,必须涉及下面基本步骤:
一、确定变量变化范围及取值间隔：
首先根据实际背景,初步统计或估计因变量,自变量变化范围,即最小值、最大值,再依据获取数据的难易程度及灵敏度、精确度的要求来确定变量点间的取值间隔，由于VensimPLE中自变量间隔不一定要求均匀,在定取值间隔时可依实际背景把变量范围看作一个阶段或若干阶段来定,在不同阶段，取值间隔可不一样，目的是必须准确反应变量间变化规律。
二、确定函数的变化趋势
根据变量间因果关系极性来确定函数的增减规律,对整个变量范围要进行分析,可能某阶段呈递增态,另一个阶段呈现递减态,有的阶段不明显，根据变化幅度大小可以调整各阶段取值间隔点间隔及点密度。
三、找出特殊点与特殊线
针对一个实际系统,建立SＤ模型，调用的表函数往往涉及到一些特殊取值点,如极值点、参照点、临界点,若有这些起点，在用VｅｎsimＰLE时，最好能直接给出在表函数的表达形式中,比如一个表函数的因变量Y是模型中其它变量乘积因子，最好找出Ｙ＝1对应的自变量X,把（X，1）直接放入表函数表达式中;如Y是模型中其它变量的和式，最好找出Ｙ=０对应的自变量X,把（Ｘ,0）直接放入表函数表达式中;如Y=Y0 是表函数单调性改变的点,最好找出相应自变量X0,把（X0，Y0）直接放入表函数表达式中,等等。
虽然表函数一般情况下是非线性的,但不排除在某阶段呈线性,若在某阶段呈线性,可在表函数表达式直接给出此阶段的两端点，（X1,Ｙ1），（X2，Y2）。
（此两点代表确定的那条直线）。
特殊线是指Y=X,Ｙ＝Y0或Yˊ=X0，的一些线段,反映了表函数的某种特征,若有，最好能直接给出。
对于如何确定特殊点、特殊线，必须据实际背景,采用各种系统工程方法来定,可参考后面表函数实例。
四、确定斜率
这里讲的斜率主要是指非特殊点的变化情况，在一个表函数中除找出特殊点外,更重要的是非特殊点的变化规律,在VｅnsimＰLE中，表函数表达式内必须列出反应变化规律的若干非特殊点,这些点通过分析历史数据、预测数据得到,直接关系到表函数是否有效，其分析方法可根据实际背景借用各种相关理论及技术手段。
为了让读者更直观的了解表函数的建立过程,下面给出二个现有SD模型的表函数进行简单分析。
例1:
如图８.2.4为一个城市建设简化模型中剩余可建面积对其事业单位新建面积的影响因子ELBC关于已占有土地比ＬＦO的表函数，即ＥLBC=f（ＬFO）,两个变量均为无量纲,此函数无法用基本初等函数表示,但可采用定性定量相结合的方法建立表函数。
第一步：
确立自变量LFO的变化范围及间距。
假设建模研究时,已建土地比是0．1,又不考虑土地余留,故LFO变化范围为[0.1,1］。
对于[0.１，1]区间如何等份?
根据整个模型的精度要求及日常实际，分成9个等份。
间隔取0.1为宜。
因若在仿真中LFO需取两位小数,也可由表函数的线性取值办法解决（因一般统计土地比只保留到两位小数）。
第二步:
确定函数增减性。
当LFO∈[0.1，0.4］时,ELBC随LＦＯ递增
当LFＯ∈[0.４，1]时,EＬBC随ＬFO递减
前式表明年新增长面积随城市建筑的增多而增加。
这是因为,在城市发展的早期阶段,大量土地有待开发，且已有企事业单位的建成会为更多企事业单位的建立创造了有利的条件。
如砖厂促进了建筑加快；路面的铺砌使材料运输更方便、更迅速；水、气、电企事业单位建立能提供生活基本保障。
另外,经实地分析，当已占土地比不超过４0%时，建地挑选也有更大余地。
出于上述定性分析得前式成立。
根据定性分析,与其它城市建立的历史事实横向比较,已占面积LFO超4０%以后，由于各种主要类型企事业单位已基本建立，市场供给网已基本建立,好的建筑环境为数不多,建筑土地资金费相应增加，这些都制约着开发地区的年建筑面积的新增,则有后式成立。
第三步:
确定特殊点。
1、由历史统计数据,当土地面积比LFO=0.1时，年新增建筑面积比为60％，则
1.7×ＥLＢC＝0．6，
ELBC=0.86
得点（０.1，０.８6）
2、有企业建设年新建面积方程分析，当LFＯ＝0.２时,ＥLBC＝1
3、确定极大值点LＦＯ=0．4时，ELBC的值，这是一个预测值,确定这种值,系统动力学本身未提供有效的方法，一般建模者常借用其它预测方法来帮助解决,如借用特尔菲方法,趋势外推，时间序列法,回归分析法、GM模型等方法。
最简单的是专家咨询法。
通过定性分析的ELBC最大值为１.1。
得特殊点（0．4，1.1）。
由实际情况,显然还有特殊点（1，0）。
第四步:
确定斜率
也就是确定非特殊点对应的ＥLBC值。
这些数据来自两部分，一部分是历史数据,另一部分是预测数据。
据系统分析综合两部分结果,参考有关资料得到了图8.2.5的表函数。
例2：
如图8．2.5为我们建立的珠海市宏观经济SD模型中建城区绿地面积第三产业影响因子关于第三产业指数的表函数。
这里第三产业指数是第三产业增加值的函数,根据珠海市的实际情况,珠海市１9９7年的建城区绿地面积占有率作为标准年,即当第三产业增加值为珠海市1997年实际值时，第三产业指数为１，这样得到了表函数的特殊点（1，1），根据珠海市第三产业的发展规模和速度，第三产业指数的取值范围定为0.5到１0,根据珠海市过去的绿地面积占有率和城市的发展理念，对比1９９7年的情况，建城区绿地面积第三产业影响因子的取值范围定为0.9到２，珠海市是一个高度重视环境绿化的城市，近几年其绿地占有率逐年上升，可以预知随着第三产业的发展，绿地面积占有率还会提高，最后达到一个较稳定的数据，通过以上分析和专家咨询，综合珠海市的发展规划得到了上面图８．2.5表函数。
该表函数可以随时根据珠海的发展进行调整，适当的时候可纳入珠海市的宏观调控计划模拟中。
此SD模型在珠海市已使用了两年，得到了各方面的肯定，表函数的功能得到了充分发挥。
第三节延迟函数
§８.3.1 VenｓimPLE中延迟函数表示及使用方法
一、定延迟函数的概念
义１ 量变化需要经过一段时间的滞后才能得到响应，这种现象称为延迟。
刻划延迟现象的函数称为延迟函数。
延迟是系统动力学中一个重要概念，因为在系统中存在大量延迟现象，例如培训的学员要经过一段时间才能发挥作用;投资要经过一段时间才能成为新的增生产能力;人得病,有潜伏期;污染物排放到江河之中,要经过扩散才能使江河发生污染等。
另外,延迟函数的构造丰富了系统动力学理论。
二、延迟函数的分类
发生的物流流线上的延迟称为物流延迟;发生在信息流线上的延迟称为 信息延迟。
根据以上概念，原则上所有的物流和信息流，在其流线上都会出现延迟,但我们在建模时,应抓住主要延迟进行设计，才能使复杂与精确性得到统一。
由物流和信息流的不同,延迟函数分为物流指数延迟函数和信息延迟函数,在VeｎsiｍPLE中其函数名有固定的表示形式，分别为ＤELAＹ1、DＥLAY2、DＥLＡＹ３和SMOＯTH、SMOOTH3等。
三、使用方法
延迟函数在VｅnsｉｍPＬE软件中直接给出，其函数在用到时可直接调用。
图8.３．１是仓存Ｉｎvent.dmｌ模型的简化流率基本入树，其中变量接到的订单ORＤRCV是关于延迟时间ＤＥL的延迟函数变量,其方程为：
接到的订单ORDRCV=ＤELＡＹ3（订单率ORＤRS,延迟时间DEL）
流位方程为:
仓存INV=ＩNTEG（接到的订单OＲＤRＣV-货运率SＨIP，期望的仓存DSIＮV）
在此该延迟变量可以不出现在入树或流图中，而可以直接放在流位方程中,这样在入树或流图中保证了流率对流位的直接作用对应关系,此时流位方程可直接写为:
仓存INＶ=INTEG（DＥLAＹ3（订单率ORＤＲS，延迟时间DEL）-货运率SHＩP，期望的仓存DＳINＶ）,其运行结果和前面一样。
下面是仓存iｎvｅnt.ｄml模型的所有方程，其中TEST是测试变量由测试函数组成,上机对TＥSＴ1、TEST2、ＴEST3、ＴＥST4中的某个赋值1，其余的值仍为零，观察运行结果可以通过测试函数了解该模型变量的实增、稳态的下降、振动和随机扰动,这样能帮助我们弄清楚模型的反馈结构及其动态行为之间的联系（有兴趣的读者可以上机试试）。
（01）ﻩ标准货运NSHIP=
1０0
ﻩUｎｉｔs：
货运单位/周
ﻩ
（02）ﻩ测试输入量TEST=
ﻩﻩＳTＥP（10，2 ）*TEＳT1+TEST2*RAMP（２０,２, 20）+TＥST３*PＵLSE（2,200）+TEST４*RANDOM ＵNＩＦＯＲM
（-5,５,０.25）
ﻩＵnｉts:
*＊undｅfinｅd**
（０3）仓存调整时间IAT=
ﻩﻩ２
ﻩＵｎiｔs:
周
ﻩ
（04）ﻩ仓存调整IＮVADJ=
ﻩﻩ（期望的仓存ＤSINV-仓存INV）/（仓存调整时间ＩAT）
ﻩUnits:
货运单位/周
ﻩ
（05）ﻩ仓存INV=INTEG（
ﻩﻩ接到的订单ORDRCＶ-货运率SＨIP,
ﻩ期望的仓存DSIＮＶ）
Unｉｔs:
货运单位
（06）ﻩ订单率ＯＲDRS=
平均货运率AVSHIP+仓存调整INVＡDＪ
Unｉtｓ：
货运单位／周
ﻩ
（07）货运率SHＩＰ=
ﻩ标准货运NSHIＰ+测试输入量TEST
ﻩUｎits:
货运单位／周
（08）ﻩ接到的订单ＯRDRCV＝
ﻩDEＬＡY3（订单率ORDＲS， 延迟时间DEL）
Units:
货运单位/周
ﻩ
（０9）平均货运率AＶＳHＩP=
ﻩＳMＯOTH（货运率ＳＨＩP， 信息延迟时间TAS ）
Ｕniｔs:
货运单位/周
ﻩ
（１0）期望的仓存DSINV=
ﻩ3*标准货运ＮSHIＰ
Ｕniｔs:
货运单位
（11）ﻩ信息延迟时间TＡS=
ﻩ２
ﻩUnits：
周
（12）ＦINALＴIME ＝2５
ﻩＵnits:
Wｅek
Ｔhefinalｔimｅｆortｈeｓiｍulatiｏn.
（13）ﻩINITIALTIＭE＝ 0
ﻩUｎｉts：
Wｅeｋ
Theiｎiｔialtｉmｅｆortｈｅｓｉmulation．
（14）延迟时间DEＬ=
3
Units:
周
ﻩ
（1５）SAVEPER =0.5
Units：
Weeｋ
The fｒequencｙ withwhich outputisstored．
（1６）TEST1＝
ﻩ０
ﻩUｎiｔｓ：
*＊ｕｎdefined**
（1７）ﻩTESＴ２=
0
Ｕnits:
＊*unｄeｆined＊*
ﻩ
（1８）ﻩTEST3=
ﻩ0
Unｉts:
＊*ｕndefinｅｄ**
ﻩ
（19）ﻩTEＳT4=
ﻩ0
ﻩUnｉｔｓ：
**uｎｄefｉned**
ﻩ
（20）ﻩＴIME ＳＴEＰ = 0.2５
Units：
Week
升级会员 