第四章第二节教学教材

第二节系统结构模型化技术

一、系统结构模型化基础

(一)结构分析的概念和意义

任何系统都是由两个以上有机联系、相互作用的要素所组成的，具有特定功能与结构的整体。结构即组成系统诸要素之间相互关联的方式。包括现代企业在内的大规模复杂系统具有要素及其层次众多、结构复杂和社会性突出等特点。在研究和解决这类系统问题时，往往要通过建立系统的结构模型，进行系统的结构分析，以求得对问题全面和本质的认识。

结构模型是定性表示系统构成要素以及它们之间存在着的本质上相互依赖、相互制约和关联情况的模型。结构模型化即建立系统结构模型的过程。该过程注重表现系统要素之间相互作用的性质，是系统认识、准确把握复杂问题，并对问题建立数学模型、进行定量分析的基础。阶层性是大规模复杂系统的基本特性，在结构模型化过程中，对递阶结构的研究是一项重要工作。

结构分析是一个实现系统结构模型化并加以解释的过程。其具体内容包括：对系统目的--功能的认识；系统构成要素的选取；对要素间的联系及其层次关系的分析；系统整体结构的确定及其解释。系统结构模型化是结构分析的基本内容。

结构分析是系统分析的重要内容，是系统优化分析、设计与管理的基础。尤其是在分析与解决社会经济系统问题时，对系统结构的正确认识与描述更具有数学模型和定量分析所无法替代的作用。

(二)系统结构的基本表达方式

系统的要素及其关系形成系统的特定结构。在通常情况下，可采用集合、有向图和矩阵等三种相互对应的方式来表达系统的某种结构。

1、系统结构的集合表达

设系统由n(n≥2)个要素(S1，S2，…，Sn)所组成，其集合为S，则有：

S={S1，S2，…，Sn}

系统的诸多要素有机地联系在一起，并且一般都是以两个要素之间的二元关系为基础的。所谓二元关系是根据系统的性质和研究的目的所约定的一种需要讨论的、存在于系统中的两个要素(Si、Sj)之间的关系Rij(简记为R)。通常有影响关系、因果关系、包含关系、隶属关系以及各种可以比较的关系(如大小、先后、

轻重、优劣等)。二元关系是结构分析中所要讨论的系统构成要素间的基本关系，一般有以下三种情形：

Si与Sj间有某种二元关系R，即SiRSj；

Si与Sj间无某种二元关系R，即Si R Sj；

Si与Sj间的某种二元关系R不明，即Si R%Sj。

在通常情况下，二元关系具有传递性，即：若SiRSj、SjRSk，则有SiRSk(Si、Sj、Sk为系统的任意构成要素)。传递性二元关系反映两个要素的间接联系，可记作R t (t为传递次数)，如何将SiRSk记作SiR2Sk。

有时，对系统的任意构成要素Si和Sj来说，既有SiRSj，又有SjRSi，这种相互关联的二元关系叫强连接关系。具有强连接关系的各要素之间存在替换性。

以系统要素集合S及二元关系的概念为基础，为便于表达所有要素间的关联方式，我们把系统构成要素中满足其种二元关系R的要素Si、Sj的要素对(Si，Sj)的集合，称为S上的二元关系集合，记作Rb，即有：

Rb={(Si,Sj)|Si、Sj∈S,SiRSj,i、j=1,2,…,n}

且在一般情况下，(Si,Sj)和(Sj,Si)表示不同的要素对。

这样，“要素Si和Sj之间是否具有某种二元关系R”，也就等价于“要素对(Si,Sj)是否属于S上的二元关系集合Rb”。

至此，我们就可以用系统的构成要素集合S和在S上确定的某种二元关系集合Rb来共同表示系统的某种基本结构。

例4—1某系统由七个要素(S1、S2、…S7)组成。经过两两判断认为：S2影响S1、S3影响S4、S4影响S5、S7影响S2、S4和S6相互影响。这样，该系统的基本结构可用要素集合S和二元关系集合Rb来表达，其中：

S={S1,S2,S3,S4,S5,S6,S7}

Rb={(S2,S1),(S3,S4),(S4,S5),(S7,S2),(S4,S6),(S6,S4)}

4—4 例4—1有向图〖TS)〗

有向图(D)由节点和连接各节点的有向弧(箭线)组成，可用来表达系统的结构。具体方法是：用节点表示系统的各构成要素，用有向弧表示要素之间的二元关系。从节点i(Si)到j(Sj)的最小(少)的有向弧数称为D 中节点间通路长度(路长)，也即要素Si 与Sj 间二元关系的传递次数。在有向图中，从某节点出发，沿着有向弧通过其它某些节点各一次可回到该节点时，在D 中形成回路。呈强连接关系的要素节点间具有双向回路。

表达例4—1给出的系统要素及其二元关系的有向图如图4—4所示。其中S3到S5、S3到S6和S7到S1的路长均为2。另外，S4和S6间具有强连接关系，S4和S6相互到达，在其间形成双向回路。

3、系统结构的矩阵表达

(1)邻接矩阵

邻接矩阵(A)是表示系统要素间基本二元关系或直接联系情况的方阵。若A=(aij)n ×n ，则其定义式为： aij= 1,SiRSj 或(Si,Sj)∈Rb(Si 对Sj 有某种二元关系)

0,Si R Sj 或(Si,Sj)∈Rb(Si 对Sj 没有某种二元关系)

有了表达系统结构的集合(S,R b)或有向图(D)，就可很容易地将A 写出，反之亦然。与例4—1和图4—4对应的邻接矩阵如下：

S1 S2 S3 S4 S5 S6 S7

A= 00000001

0000000

0010000000110000000000010000100000??????????????????????

很明显，A 中“1”的个数与例4—1中Rb 所包含的要素对数目和图4—4中有向弧的条数相等，均为6。在邻接矩阵中，若有一列(如第j 列)元素全为0，则Sj 是系统的输入要素，如图4—4中的S3和S7；若有一行(如第i 行)元素全为0，则Si 是系统的输出要素，如图4—4中的S1和S5。

(2)可达矩阵

S1 S2 S3 S4 S5 S6 S7

若在要素Si 和Sj 间存在着某种传递性二元关系，或在有向图上存在着由节点i 至j 的有向通路时，称Si 是可以到达Sj 的，或者说Sj 是Si 可以到达的。所谓可达矩阵(M)，就是表示系统要素之间任意次传递性二元关系或有向图上两个节点之间通过任意长的路径可以到达情况的方阵。若M=(mij)n ×n ，且在无回路条件下的最大路长或传递次数为r ，即有0≤t ≤r ，则可达矩阵的定义式为： mij= 1，SiR t Sj (存在着i 至j 的路长最大为r 的通路)

0,Si t R Sj (不存在i 至j 的通路)

当t=1时，表示基本的二元关系，M 即为A ；当t=0时，表示Si 自身到达，或SiRSi ，也称反射性二元关系；当t ≥2时，表示传递性二元关系。

矩阵A 和M 的元素均为“1”或“0”，是n ×n 阶0—1矩阵，且符合布尔代数的运算规则，即：0+0=0,0+1=1,1+0=1,1+1=1,0×0=0,0×1=0,1×0=0,1×1=1。通过对邻接矩阵A 的运算，可求出系统要素的可达矩阵M ，其计算公式为：

M=(A+I)r (4—1)

其中I 为与A 同阶次的单位矩阵(即其主对角线元素全为“1”，其余元素为“0”)，反映要素自身到达；最大传递次数(路长)r 根据下式确定：

(A+I)≠(A+I)2≠(A+I)3≠…≠(A+I)r-1≠(A+I)r =(A+I)r+1=…=(A+I)n

(4—2) 以与例4—1和图4—4对应的邻接矩阵为例有：

S1 S2 S3 S4 S5 S6 S7

A+I= 1000000110000100110000

00111000001000

0010100100001?????????????????????? 其中主对角线上的“1”表示诸要素通过零步(自身)到达情况(单位矩阵I)，其余“1”表示要素间通过一步(直接)到达情况(邻接矩阵A)。

S1 S2 S3 S4 S5 S6 S7

S1 S2 S3 S4 S5 S6 S7

(A+I)2=A 2+A+I= ()()()()100000011000000011`00

00111000001000

00110100001????????I I ????????I ??

??I ????

其中带圆圈的“1”表示要素间通过两步(间接)到达情况(矩阵A 2)。按照前述布尔代数的运算规则，在原式(A+I)2的展开中利用了A+A=A 的关系。

进一步计算发现：(A+I)3=(A+I)2。由(4—2)式即有r=2。

这样，根据(4—1)式，与例4—1和图4—4对应的可达矩阵为：

S1 S2 S3 S4 S5 S6 S7

M=(A+I)2= 1000000110000000111100

00111000001000

0011101100001??????????????????????

(3)其它矩阵

在邻接矩阵和可达矩阵的基础上，还有其它表达系统结构并有助于实现系统结构模型化的矩阵形式，如缩减矩阵、骨架矩阵等。

①缩减矩阵

根据强连接要素的可替换性，在已有的可达矩阵M 中，将具有强连接关系的一组要素看作一个要素，保留其中的某个代表要素，删除掉其余要素及其在M 中的行和列，即得到该可达矩阵M 的缩减矩阵M ′。如原例可达矩阵的缩减矩阵为：

S1 S2 S3 S4 S5 S6 S7 S1 S2 S3 S4 S5 S6 S7

S1 S2 S3 S4 S5 S7

M ′= 1000001

100000

011100001100000101

10001??????????????????????

②骨架矩阵

对于给定系统，A 的可达矩阵M 是唯一的，但实现某一可达矩阵M 的邻接矩阵A 可以具有多个。我们把实现某一可达矩阵M 、具有最小二元关系个数(“1”元素最少)的邻接矩阵叫M 的最小实现二元关系矩阵，或称之为骨架矩阵，记作A ′。

系统结构的三种基本表达方式相互对应，各有特色。用集合来表达系统结构概念清楚，在各种表达方式中处于基础地位；有向图形式较为直观、易于理解；矩阵形式便于通过逻辑运算，用数学方法对系统结构进行分析处理。以它们为基础和工具，通过采用各种技术，可实现复杂系统结构的模型化。

(三)常用系统结构模型化技术

系统结构模型化技术是以各种创造性技术为基础的系统整体结构的决定技术。它们通过探寻系统构成要素、定义要素间关联的意义、给出要素间以二元关系为基础的具体关系，并且将其整理成图、矩阵等较为直观、易于理解和便于处理的形式，逐步建立起复杂系统的结构模型。常用的系统结构模型化技术有：关联树法、解释结构模型化技术、系统动力学结构模型化技术等，其中解释结构模型化(ISM)技术是最基本和最具特色的系统结构模型化技术。

ISM 技术是美国J ·N ·沃菲尔德教授于1973年作为分析复杂的社会经济系统结构问题的一种方法而开发的。其基本思想是：通过各种创造性技术，提取问题的构成要素，利用有向图、矩阵等工具和计算机技术，对要素及其相互关系等信息进行处理，最后用文字加以解释说明，明确问题的层次和整体结构，提高对问题的认识和理解程度。该技术由于具有不需高深的数学知识、模型直观且有启发性、可吸收各种有关人员参加等特点，因而广泛适用于认识和处理各类社会经济S1 S2 S3 S4 S5 S7

系统的问题。ISM的基本工作原理如图4—5所示。

图4—5 ISM工作原理图

由图4—5可知，实施ISM技术，首先是提出问题，组建ISM实施小组；接着采用集体创造性技术，搜集和初步整理问题的构成要素，并设定某种必须考虑的二元关系(如因果关系)，经小组成员及与其他有关人员的讨论，形成对问题初步认识的意识(构思)模型。在此基础上，实现意识模型的具体化、规范化、系统化和结构模型化，即进一步明确定义各要素，通过人机对话，判断各要素之间的二元关系情况(即SiRSj?)，形成某种形式的“信息库”；根据要素间关系的传递性，通过对邻接矩阵的计算或逻辑推断，得到可达矩阵；将可达矩阵进行分解、缩约和简化处理，得到反映系统递阶结构的骨架矩阵，据此绘制要素间多级递阶有向图，形成递阶结构模型；通过对要素的解释说明，建立起反映系统问题某种二元关系的解释结构模型。最后，将解释结构模型与人们已有的意识模型进行比较，如不相符合，一方面可对有关要素及其二元关系和解释结构模型的建立进行修正；更重要的是，人们通过对解释结构模型的研究和学习，可对原有的意识模型有所启发和修正。经过反馈、比较、修正、学习，最终得到一个令人满意、具有启发性和指导意义的结构分析结果。

通过对可达矩阵的处理，建立系统问题的递结构模型，这是ISM技术的核心

内容。根据问题规模和分析条件，可在掌握基本原理及其规范方法的基础上，采用多种手段、选择不同方法来完成此项工作。

二、建立递阶结构模型的规范方法

建立反映系统问题要素间层次关系的递阶结构模型，可在可达矩阵M 的基础上进行，且一般要经过区域划分、级位划分、骨架矩阵提取和多级递阶有向图绘制等四个阶段。这是建立递阶结构模型的基本方法。现以例4—1所示问题为例说明。与图4—4对应的可达矩阵(其中将Si 简记为i)为：

1 2 3 4 5 6 7

M= 1000000110000000111100

00111000001000

0011101100001??????????????????????

1、区域划分

区域划分即将系统的构成要素集合S ，分割成关于给定二元关系R 的相互独立的区域的过程。

为此，需要首先以可达矩阵M 为基础，划分与要素Si(i=1,2,…,n)相关联的系统要素的类型，并找出在整个系统(所有要素集合S)中有明显特征的要素。有关要素集合的定义如下：

(1)可达集R(Si)

系统要素Si 的可达集是在可达矩阵或有向图中由Si 可到达的诸要素所构成的集合，记为R(Si)。其定义式为：

R(Si)={Sj|Sj ∈S,m ij =1,j=1,2,…,n} i=1,2,…,n

如在给出的可达矩阵中有：R(S1)={S1},R(S2)={S1,S2}，R(S3)={S3,S4,S5,S6},R(S4)=R(S6)={S4,S5,S6},R(S5)={S5},R(S7)={S1,S2,S7}。

(2)先行集A(Si)

1 2 3 4 5 6 7

系统要素Si 的先行集是在可达矩阵或有向图中可到达Si 的诸系统要素所构成的集合，记为A(Si)。其定义式为：

A(Si)={Sj|Sj ∈S,m ji =1,j=1,2,…,n}i=1,2,…,n

如在给出的可达矩阵中有：A(S1)={S1,S2,S7},A(S2)={(S2,S7),A(S3)}={S3},A(S4)=A(S6)={S3,S4,S6},A(S

5)={S3,S4,S5,S6},A(S7)={S7}。

(3)共同集C(Si)

系统要素Si 的共同集是Si 在可达集和先行集的共同部分，即交集，记为C(S i)。其定义式为：

C(Ｓi)={Sj|Ｓj ∈S,m ij =1, m ij =1,j=1,2,,n} i=1,2，… ,n 如：

C(S1)={S1},C(S2)={S2},C(S3)={S3},C(S4)=C(S6)={S4,S6},C(S5)={S5},C(S7)={S7}。

系统要素Si 的可达集R(Si)、先行集A(Si)、共同集C(Si)之间的关系如图4—6所示。

(4)起始集B(S)和终止集E(S)

图4—6 可达集、先行集、共同集关系示意图

系统要素集合S 的起始集是在S 中只影响(到达)其他要素而不受其他要素影响(不被其他要素到达)的要素所构成的集合，记为B(S)。B(S)中的要素在有向图中只有箭线流出，而无箭线流入，是系统的输入要素。其定义式为：B(S)={Si|Si ∈S,C(Si)=A(Si),i=1,2,…,n}。

如在与图4—4所对应的可达矩阵中，B(S)={S3,S7}。 R(Si) A(Sj) C(Sj) Sj