(40)住房建设 - (49)住房折旧 -

Uses: (34)土地占用因子 - (40)住房建设 - (49)住房折旧 -

(51)住户数与住房比率 -

(51)

住户数与住房比率 = 人口总量/(家庭人口数*住房总量)

Units: 户/房 (31)人口总量 - (50)住房总量 - (08)家庭人口数 -

Uses: (42)住房可得性因子 - (46)住房吸引力因子 -

SAVEPER = 1 Units: 年

The frequency with which output is stored.

TIME STEP = 1 Units: 年

The time step for the simulation.

(52)

(53)

6.3模型模拟

先按原模型数据模拟一次，再改变住房土地因子表和企业土地因子表模

拟一次，得到数据文件city.vdf和city1.vdf。由运行参数比较获得以下：

Comparing CITY and CITY1

******Constant differences between CITY and CITY1****** FINAL TIME - has changed in value 50 CITY 100 CITY1

-----------------------------------------------------------

******Lookup differences between CITY and CITY1******

企业土地系数表 - has changed in value

CITY CITY1 X | Y X | Y 0 | 1 -0.005155 | 0.7254 0.1 | 1.15 0.1289 | 0.838 0.2 | 1.3 0.2294 | 0.9437 0.3 | 1.4 0.3093 | 1.049 0.4 | 1.45 0.4175 | 1.19 0.5 | 1.4 0.4948 | 1.282 0.6 | 1.3 0.567 | 1.394 0.7 | 0.9 0.6211 | 1.493 0.8 | 0.5 0.6985 | 1.585 0.9 | 0.25 0.7448 | 1.62 1 | 0 0.7938 | 1.655

17

0.8376|1.634 0.8686|1.535 0.8969|1.373 0.9046|1.169 0.9072|1.155 1|0

住房土地因子表 - has changed in value

CITY CITY1 X | Y X | Y 0 | 0.4 0 | 0.4 0.1 | 0.7 0.1 | 0.7 0.2 | 1 0.2 | 1 0.3 | 1.25 0.2655 | 1.282 0.4 | 1.45 0.3247 | 1.493 0.5 | 1.5 0.3814 | 1.606 0.6 | 1.5 0.3814 | 1.599 0.7 | 1.4 0.4562 | 1.606 0.8 | 1 0.5541 | 1.415 0.9 | 0.5 0.634 | 1 1 | 0 0.7088 | 0.6549 0.8299|0.3451 1|0

6.3数据集分析

下图分别是住房总量随时间的变化图，它的直接决定变量随时间的变化图,及它所直接决定的变量随时间的变化图。 CITYCITY1住房总量200,000100,0000025Time (年)50

18

CITYCITY1住房总量200,000100,0000住房建设4,0002,0000住房折旧2,0001,0000012.525Time (年)37.550住房初始总量CITY: 14,000 CITYCITY1住房总量200,000100,0000土地占用因子21住房建设4,0002,0000住房折旧2,0001,0000住户数与住房比率1.50012.50CITY1: 14,000

25Time (年)37.550

第7章 真实性检验

证实所建立的模型的合理性和有效性，有好多方法。而Vensim所提供的真

实性检验(Reality Check)是一种非常有效的方法和工具。它也是Vensim的一个主要特征。

19

7.1真实性检验的概念

模型是对于现实存在的系统的简化和抽象。我们可以由建立的模型方便的得到系统的行为模式。但如何确定这些行为模式正确性？常用的一种方法是对照现实的情况，仔细分析运行得到的各个变量时间序列数据及其相互关系。在系统复杂、变量多情况下，这是很困难的。一种简单有效的方法是首先确定模型的正确性。如果模型不违反基本的常识和规则，模型本身的合理性得到了检验，则系统的行为模式的真实性就有了保证。真实性检验由这一思想出发，提供了一种简便和有效的方法和工具。

根据要研究的系统的基本情况，我们一般可以提出模型应该遵守的一些基本假定。这些基本假定是保证模型正确的必要条件，构成了对于模型正确性的一系列约束。真实性检验就是在模型建立后，通过标准的约束描述语言，具体地描述这些约束，并通过运行模型，自动的检验系统的行为模式是否违反这些约束。若违反，说明模型尚不完善，则要调整模型的结构和参数，直到符合这些约束条件。真实性检验的描述部分完全独立于模型本身，不会改变模型结构和行为模式。

7.2模型真实性检验的定义

1.变量的命名

检验语句是对于假设的描述，是由实际系统出发对于模型的行为模式的断言。它一般包含的变量是“判断句”。有别于一般的变量用名词命名，它一般用表示判断的短语命名。例如：no worker no production, rain means flooding。 2. 输入测试变量(Test Input) 用来代替一个等式。它只能用作约束的前提条件，可以使约束变得简单易读。在约束表达式中，在其条件部分出现以一个逻辑变量形式出现。其定义形式为： name :TEST INPUT: variable = expression

其中，name表示变量名，:TEST INPUT: 是变量类型标识，而后面等式是定义的对象。等式右边可以是允许的任意表达式。等式左右边的构成变量均来自模型。

例如：best condition :TEST INPUT: temperature = 3*initial temperature 3. 时间演进函数(Time Transition) 该函数的自变量是时间，而函数值是随时间变化的变量值。其一般表达式为： variable = TIME TRANSITION (X，… ，Y ，….)

其中X是时间，而函数值variable按一定的对应规则取相应的Y。在Y的列举中，可以有“-”，表示保持模型的值，而“*”表示采用模型的计算值。实际上，Y是我们对于模型行为模式随时间变化的一种假定的表示。 主要形式有以下几例： profit = TIME TRANSITION (0,30,80,100,20,50,50,15)

表示profit随时间变化的折线图。时间由0、30、80到100，profit由20斜线上升到50，保持50不变，然后斜线下降到15。 profit = TIME TRANSITION (30,80,100,-,50,*)

表示profit在时间30以前采用模型正常值A；在时间区间30-80，profit采用以(30,A)为起点，以(80,50)为终点直线上的值；在时间区间80-100，将采用模型的计算值。 profit = TIME TRANSITION (30,80,100,*,50,*)

表示在时间区间30-80，profit采用模型的计算值，在时间为80时值为50。

20