复杂超图

摘要

提供自然和人类复杂结构的抽象表示是一个具有挑战性的问题。 在考虑系统异构组件的同时考虑到分析的易处理性是一个困难的平衡。 在这里，我介绍复杂超图（chygraph），汇集了超图、多层网络、单纯复形和超结构的概念。 为了说明这种组合结构的适用性，我计算了组件大小统计数据并确定了向巨型组件的过渡。 为此，我介绍了一种矢量化技术，可以解决 chygraph 的多层次性质。 我的结论是，chygraph 是复杂系统的统一表示，可以进行分析洞察。

我简介

图是多粒子系统的结构抽象，顶点代表粒子，边代表粒子之间的成对相互作用Albert and Barabási (2002)。 超图是允许两个或多个顶点之间交互的扩展Ghoshal 等人 (2009)；巴斯克斯（2009）；库蒂尼奥等人 (2020); Civilini 等人 (2021);孙和比安科尼（2021）。 多重网络引入了不同交互类型的层Gómez 等人 (2013)。 单纯复形将连通性扩展到包含层次结构Petri 和 Barrat (2018)；考特尼和比安科尼 (2018)； Iacopini 等人 (2019)。

也有自指结构。 Joslyn 和 Nowak 引入了 ubergraphs Joslyn and Nowak (2017)，这是一种边可以包含其他边的组合结构。 ubergraphs 已用于组织知识数据库中的关系Yadati 等人 (2021)。 在更高的抽象层次上，Baas 使用范畴论 Baas (2019) 中的概念引入了高阶结构，也称为超结构。

在实践中，什么是合适的表示是一个平衡问题。 我们希望超图的灵活性能够超越成对相互作用、多层的可能性以及单纯复形和超结构的分层包含结构。 反过来，我们希望图表能够简单地计算聚合属性，例如巨型组件的大小。

在这里我介绍复杂超图作为一种灵活的组合结构。 足够通用，包括大量以前的工作和新的应用程序。 足够简单，可以进行分析计算。 现将工作安排如下。 昆虫。 II 我定义了复杂的超图并评论了它与其他结构的关系。 昆虫。 III 和 IV 我使用生成函数形式提供了分析方法来表征 chygraph 巨型组件的出现。 我讨论了特殊情况并提供了验证分析结果的数值示例。 我在第 V 节中得出了一些关键结论。

II 关键定义

II.1 动机

复杂的系统被简化以方便分析。 逐渐地，我们删除了一些简化并更接近真实的系统。 例如，科学出版物系统根据所提出的问题由不同的网络结构表示Newman (2001)； Barabási 等人 (2002)。 引文网络表明知识沿着出版物的流动Redner (1998)； Lehmann等人(2003)。 在引文网络中，节点是出版物，引文由定向链接表示。 当专注于合作时，合着网络更适合。 作者-出版物网络可以进一步扩展以明确表示作者和出版物，从而产生二分图Newman (2001)。 相同的作者-出版物关系可以表示为作者超图，其中出版物是关联一位、两位或更多作者 Vasilyeva 等人 (2021) 的超边。 这些网络和超图是简化的，失去了原始系统的某些方面。 我们需要更完整的表示和更丰富的组合结构。 科学出版物包含作者和参考文献。 文档内部结构可以用具有两条边的超图表示，即作者和参考文献列表。 最重要的是，出版物是高阶结构中的一个顶点，其中的构建块是作者和出版物。 通俗地说，是超图的超图，复杂的超图。

II.2 复杂系统

在这里，我在结构意义上使用“复杂”这个术语：由不同的部分组成。 我区分了不可分解为其他部分的部分、原子以及由其他部分（包括原子）组成的复合体。 原子可以具有更精细的结构，但它们已被选为主要构建块。 这些基本知识导致了复杂系统的自洽定义。

定义。复杂系统是原子和络合物的集合，其中络合物是由原子和其他络合物组成的。

后者作为复杂系统的哲学定义。 对于实际应用，我们需要精确的数学结构。 选择取决于复合体的内部结构。

II.3 超图

超图$\mathscr{H}(V,E)$是顶点集$V$和超边集$E$，其中超边是$V$的子集。为了便于说明，我根据上面介绍的复杂系统的措辞重写了这个定义。 超图$\mathscr{H}(A,C)$是一个原子集$A$和一个复数集$C$，其中复数是$A$的子集（超边）。 显然，超图不具备复杂系统的自指属性。 它们的结构很简单。

II.4 超级图

为了使超图具有自指性，Joslyn 和 Nowak 引入了超级图 Joslyn and Nowak (2017)，这是一种超边可以包含其他超边的组合结构。 这里我用上面介绍的复杂系统的措辞重新定义它们。

定义。 ubergraph $𝒰(A,C)$ 是一组原子 $A$ 和一组复合体 $C$，其中复合体是 $A\cup C$。

图1： 一个 chygraph 示例。 该 chygraph 是科学出版物系统的表示，其中原子代表作者、复合出版物和包含引用。

II.5 复杂超图

为了添加更多内部结构，我选择超图。 这就引出了复杂超图的定义。

定义。复杂超图（chygraph） $\chi (A,C)$是一组原子$A$和一组复合体$C$，其中复合体是具有顶点集的超图在$A\cup C$中。

让我们用一些例子来解释这个定义。 图 $𝒢(V,E)$ 由 chygraph $\chi (V,E)$ 表示，其中复合体是边。 具有分层图 $G_l(V,E_l)$、$l=1,\mathrm{\dots },L$ 的多重图 Gómez 等人 (2013) 由 chygraph $\chi (V,{\cup }_{l=1}^L{E}_l)$ 加上一些分区表示下面讨论的结构。 科学出版物系统由 chygraph $\chi (A,\{{\mathscr{H}}_i(A_i\cup R_i,\{A_i,R_i\})\})$ 表示，其中原子是作者，复合体是出版物，出版物由具有两条边的超图表示($A_i$ 代表作者，$R_i$ 用于参考），索引 $i$ 贯穿所有出版物。 两条边没有重叠，因此代表科学出版物的复合体具有两个内部组件（图1)。

我们可以将 chygraph 视为具有两种边缘类型的超级图：与超图相关的复合体中的边和复合体。 该复合体充当一个特殊的边缘，我们将其指定为包含的基本单元。 反过来，超级图是 chygraph，其中所有复杂内超图都包含一条边。 显然，在这种情况下，复杂内超图只有一个其大小等于超边基数的分量。 因此，在用超边基数替换组件大小后，以下部分中得出的所有结果对于超级图均有效。

II.6 Cygraphs 属性

图/超图的许多属性都被延续到了 chygraph 中。 为了区别与复杂超图结构相关的度量，我将使用希腊字母来命名 chygraph 级别的数量。 给定一个复杂的 $j$，令 $\mathscr{H}(V_j,E_j)$ 为其关联的超图。 令 $i\in A\cup C$ 为原子或复合物，其中 $i=1,\mathrm{\dots },n$ 和 $n=|A\cup C|$。 我将使用与 $i\in V_j$ 等效的符号 $i\in C_j$。 chy-邻接矩阵 $\alpha$ 是具有矩阵元素的 $n\times n$ 矩阵

(1)

其中$i,j\in A\cup C$。 关联的顶点 chy 度

$${\kappa }_i=\sum _j{\alpha }_{ij}.$$ (2)

需要额外的指标来表征潜在的多层结构。 Chygraph 可能包含多个级别的包含内容。 在较低的层次上，我们有原子，其精细结构为空或未指定。 在上面一层我们有复合体，它们的内部结构被指定为包含原子和/或其他复合体作为顶点的超图。 在某些系统中，定义更高级别的内含物是有意义的。 例如，当按地点描述人口时，我们会谈论社区、城市、国家、大陆和世界。 这种包含层次结构导致了 chygraph 长度的定义。

定义：chygraph长度$L(\chi )$。 令 $\chi (A,C=\{{\mathscr{H}}_i(V_i,E_i)\})$ 为 chygraph。 令 $\mathrm{\Pi }={\mathrm{\Pi }}_1\cup \mathrm{\dots }\cup {\mathrm{\Pi }}_l$ 为 $C$ 的一个分区，它是非拦截的(${\mathrm{\Pi }}_i\cap {\mathrm{\Pi }}_j=\mathrm{\varnothing }$ 代表 $i\ne j$)、分层（如果 $C_i\in {\mathrm{\Pi }}_j$ 然后 $V_i\subset A\cup ({\cup }_{k\le j}{\mathrm{\Pi }}_k)$)并且同一分区内的复合体具有相似的统计特性。 chygraph 长度由 $L(\chi )$ 表示，是所有此类分区中的最大值 $l$。

通过统计属性区分分区允许指定多类型结构。 多重图和超图就是这种情况。 在这种情况下，两个层可以具有相同的顶点集，但是根据层的不同，复合体可以具有不同的统计属性。

与高阶网络 De Domenico 等人 (2013) 的情况一样，我们可以将聚类系数和熵等其他指标扩展到 chygraph。 在这里，我将注意力限制在以下部分所需的指标上。

III渗流理论一：无层内夹杂物

尽管其复杂性，chygraph组合结构仍适合分析处理。 我专注于渗透，这是图论的一个核心问题。

生成函数技术已被用来解决图和超图中的渗流问题 Callaway 等人 (2000); Ghoshal 等人 (2009);库蒂尼奥等人 (2020);孙和比安科尼（2021）。 它遵循一个简单的方法：将组件尺寸分布的生成函数表示为自身的递归函数，并由其他相关分布的生成函数进行调制。 在糜烂图中，关键量是组件大小 ${\sigma }^l$ 和多余组件大小 ${\overline{\sigma }}^{ml}$，当组件分别从层 $l$ 随机取样或从另一层 $m$ 取样时。 这里我们需要精确地了解我们如何从另一层到达。

我们可以通过两种不同的方式得出参考复合体。 来自它包含的复合体（来自下面）或来自它所包含的复合体（来自上面）。 如果没有层内连接，则符号 ${\overline{\sigma }}^{ml}$ 不需要任何进一步的规范。 当时，很明显我们是从下面到达的。 相反，当$m>l$时，我们是从上面到达的。 在本节中，我假设没有层内连接。 层内连接的情况将是下一节的主题。

组件大小取决于 chy 度 ${\widehat{\kappa }}_l$ 的联合分布、复杂超图组件大小 ${\check{s}}^l$ 及其多余等价 ${\check{\overline{\kappa }}}_l^m$ 和 ${\check{\overline{s}}}_l^m$ 从层 $m$ 到达时。符号${\check{s}}_l=(s_{l0},s_{l1}\mathrm{\dots },s_{ll-1}$)表示$l$层的复合体内的组件由来自$A,{\mathrm{\Pi }}_1,\mathrm{\dots },{\mathrm{\Pi }}_{l-1}$的顶点组成。 反过来，${\widehat{\kappa }}_l=({\kappa }_{ll+1},\mathrm{\dots },{\kappa }_{lL})$表示$l$层的顶点的chy-度被分解为${\mathrm{\Pi }}_{l+1},\mathrm{\dots },{\mathrm{\Pi }}_L$层的顶点的chy-度。 ${\widehat{\kappa }}_l$、${\check{\overline{\kappa }}}_{ml}$、${\check{s}}^l$、${\check{\overline{s}}}^{ml}$、${\sigma }^l$和${\overline{\sigma }}^{ml}$的概率生成函数分别用${\mathrm{\Phi }}^l({\widehat{x}}_l)$、${\mathrm{\Psi }}^{ml}({\widehat{x}}_l)$、$G^l({\check{x}}_l)$、$U^{ml}({\check{x}}_l)$、${\mathrm{\Gamma }}^l(x)$和${\mathrm{{\rm Y}}}^{ml}(x)$表示。 由于它们是概率分布的生成函数，因此在 $x=1$ 处求值时它们都等于 1，并且它们的一阶导数等于相应的期望值。

III.1 平均组件尺寸

chygraph 的定义被转化为组件尺寸生成函数的一组自洽方程。 作为指导，我们到达一个参考复合体，然后从复合体中我们导航 chygraph。 遵循 chygraph 邻接矩阵，引导我们找到包含参考复合体的复合体，或者遍历参考复合体超图中包含的复合体。 更确切地说，

${\mathrm{\Gamma }}^l(x)$ $=$ $x{\mathrm{\Phi }}^l[{\mathrm{{\rm Y}}}^{ll+1}(x),\mathrm{\dots },{\mathrm{{\rm Y}}}^{lL}(x)]$ (3)

$\times$ $G^l[{\mathrm{{\rm Y}}}^{l0}(x),\mathrm{\dots },{\mathrm{{\rm Y}}}^{ll-1}(x)],$

${\mathrm{{\rm Y}}}^{ml}(x)$ $=$ $x{\mathrm{\Psi }}^{ml}[{\mathrm{{\rm Y}}}^{ll+1}(x),\mathrm{\dots },{\mathrm{{\rm Y}}}^{lL}(x)]$ (4)

$\times$ $U^{ml}[{\mathrm{{\rm Y}}}^{l0}(x),\mathrm{\dots },{\mathrm{{\rm Y}}}^{ll-1}(x)],$

其中$l=0,\mathrm{\dots },L(\chi )$。 请注意，这些方程右侧的第一个 $x$ 对应于参考复合体，包含 ${\mathrm{\Phi }}^l[\mathrm{\cdots }]$ 和 ${\mathrm{\Psi }}^{ml}[\mathrm{\cdots }]$ 的项使用 chygraph 邻接矩阵进行导航（包括参考复合体的复合体）以及包含 $G^l[\mathrm{\cdots }]$ 和 $U^{ml}[\mathrm{\cdots }]$ 的术语，用于在参考复合体超图中导航复合体。

III.2 平均组件尺寸

平均多余元件尺寸${⟨\overline{\sigma }⟩}^{ml}={\dot{\mathrm{{\rm Y}}}}^{ml}(1)$可以通过等式计算。 (4)，导致

$${⟨\overline{\sigma }⟩}^{ml}=1+\sum _{k=l+1}^L{⟨\overline{\kappa }⟩}_{lk}^m{⟨\overline{\sigma }⟩}^{lk}+\sum _{k=0}^{l-1}{⟨\overline{s}⟩}_{lk}^m{⟨\overline{\sigma }⟩}^{lk},$$ (5)

其中$l,m=0,\mathrm{\dots },L$。 请注意，仅当 $k=m$ 时才出现 ${⟨\overline{\kappa }⟩}_{lk}^m\ne {⟨\kappa ⟩}_{lk}$ 和 ${⟨\overline{s}⟩}_{lk}^m\ne {⟨s⟩}_{lk}$。 当我们从层$m$进入层$l$，然后返回层$m$时。现在是矢量化技巧。

矩阵方程(5)可以通过矩阵方程Horn等人(1994)的矢量化方法求解。 向量化运算符 $\mathrm{vec}X$ 通过堆叠 $X$ 的列，将 $(M,N)$ 矩阵转换为 $M\times N$ 列向量。 例如，

$$\mathrm{vec}{X}^{(2,2)}=\left[\begin{array}{c}{X}^{00}\\ X^{10}\\ X^{01}\\ X^{11}.\end{array}\right]$$ (6)

为了处理具有四个索引的张量，我概括了向量化运算符。 作用于 $X_{(O,P)}^{(M,N)}$ 张量的矢量化运算符通过沿列堆叠上部索引和沿行堆叠下部索引来创建 $(M\times N,O\times P)$ 矩阵。 例如，

(7)

应用矢量化方程。 (5) 转换为

$$\mathrm{vec}\{A\}\mathrm{vec}\{⟨\overline{\sigma }⟩\}=\mathrm{vec}\{B\},$$ (8)

在哪里

$$A_{nk}^{ml}={\delta }_{mn}{\delta }_{lk}-{⟨\overline{\kappa }⟩}_{nk}^m{\mathrm{\Theta }}_{k-n}{\delta }_{nl}-{⟨\overline{s}⟩}_{nk}^m{\mathrm{\Theta }}_{n-k}{\delta }_{nl},$$ (9)

$$B^{ml}=1,$$ (10)

${\mathrm{\Theta }}_i$ 是整数 Heaviside 阶跃函数（如果 $i>0$ 则为 ${\mathrm{\Theta }}_i=1$，否则为 0）。 线性方程组 (8) 将克拉默规则作为形式解

$$\mathrm{vec}{\{⟨\overline{\sigma }⟩\}}_i=\frac{det(\mathrm{vec}{\{A\}}_i)}{det(\mathrm{vec}\{A\})},$$ (11)

其中 $\mathrm{vec}{\{A\}}_i$ 派生自 $\mathrm{vec}\{A\}$，通过用 $\mathrm{vec}\{B\}$ 替换第 $i$ 列。 如果 $A$ 不是单数，则该解决方案有效。 当$det(A)\to 0^{+}$时，平均成分尺寸发散，系统实现渗透。 因此，chygraph临界渗透条件由下式给出

$$det(\mathrm{vec}\{A\})=0.$$ (12)

III.3 巨型组件

平均组件尺寸 (8) 的公式在 $det(\mathrm{vec}\{A\})>0$ 条件下有效。 下面我将证明 $det(\mathrm{vec}\{A\})>0$ 对应于亚临界相。 令 $P^l$ 为随机选择的层 $l$ 中的顶点不属于巨型组件的概率，令 $Q^{ml}$ 为层中的顶点不属于巨型组件的概率从$m$层的复合体中选择的$l$不属于巨型组件。 这些概率满足自洽方程

$$P^l={\mathrm{\Phi }}^l[Q^{ll+1},\mathrm{\dots },Q^{lL}]G^l[Q^{0l},\mathrm{\dots },Q^{l-1l}],$$ (13)

$$Q^{ml}={\mathrm{\Psi }}^{ml}[Q^{ll+1},\mathrm{\dots },Q^{lL}]U^{ml}[Q^{0l},\mathrm{\dots },Q^{l-1l}].$$ (14)

该方程组没有显式解析解。 可以通过逐次逼近找到解决方案，其中将迭代 $t$ 插入右侧后，将左侧解释为 $t+1$ 迭代。 特别是，在没有巨大组件的情况下，方程。 (13) 和 (14) 承认解$P^l=Q^{ml}=1$。 让我们假设 $Q^{lm}=1-x^{ml}$，其中 $x^{ml}\to 0$。 将等式中 $x^{ml}$ 中的项保持为一阶。 (14) 得出递归近似方程

$$\mathrm{vec}\{x\}(t+1)=(I-\mathrm{vec}\{A\})\mathrm{vec}\{x\}(t).$$ (15)

当且仅当 $\mathrm{\Lambda }(\mathrm{vec}\{A\})>0$ 时，线性映射 (15) 收敛于 $\mathrm{vec}\{x\}=0$，其中 $\mathrm{\Lambda }(\mathrm{vec}\{A\})$ 是 $\mathrm{vec}\{A\}$。 因此$\mathrm{\Lambda }(\mathrm{vec}\{A\})$是巨型组件存在的控制参数。 在亚临界（超临界）相 $\mathrm{\Lambda }>0$ () 中，不存在（存在）巨型组件，且 $P^l=1$ () 。 渗流转变发生在临界条件$\mathrm{\Lambda }(\mathrm{vec}\{A\})=0$，这相当于等式： (12)。

三.4 $L=1$

当只有一个复杂层($L(\chi )=1$)时，则式（1）中的$\mathrm{vec}\{A\}$。 (9) 减少为

(16)

具有行列式

$$det(\mathrm{vec}\{A\})=1-A_{01}^{10}{A}_{10}^{01}=1-{\overline{k}}_{01}^1{\overline{s}}_{10}^0.$$ (17)

超图 $\mathscr{H}(V,E)$ 映射到一层 chygraph，其中复合体是超图边缘：$\chi (V,\{\mathscr{H}(V,e_l),e_l\in E\})$。 在这种情况下，多余的组件大小${\overline{s}}_{10}^0$是多余的超边基数$⟨\overline{c}⟩$，多余的原子度${\overline{k}}_{01}^1$是顶点的多余度$⟨\overline{k}⟩$ 。 代入等式。 (17) 得到超图的临界条件：$⟨\overline{c}⟩⟨\overline{k}⟩=1$，与 Coutinho et al Coutinho 等人 (2020 ）。 此外，图是基数超过 1 的超图，因此，临界条件降低为 $⟨\overline{k}⟩=1$，如 Molloy 和 Reed Molloy 和 Reed (1998) 以及 Callaway et al Callaway 等人 (2000).

三.5 $L=2$

当有两个复杂层($L(\chi )=2$)时，则方程中的$\mathrm{vec}\{A\}$。 (9) 减少为

(18)

具有行列式

(19)

我们可以扩展行列式，但这会导致复杂的代数表达式。 对于大多数实际目的，最好使用方程式中的显式矩阵形式。 (19)，包括数值计算。 当我们考虑具有附加约束的组合结构时，会得到进一步的简化。

III.6分层包含

在许多系统中，包容性遵循层次结构。 地理缩放，例如 Okabe 和 Sadahiro (1996)。 在 chygraph 中构建分层包含很简单：原子包含在第 1 层复合体中，第 1 层复合体包含在第 2 层复合体中，……层 $L-1$ 复合体包含在层 $L$ 复合体中。 在这种情况下，如果 $|m-l|\ne 1$ 或 $|n-k|\ne 1$，则 $A_{nk}^{ml}$ 的张量元素为零。对于 $L=2$ 的情况，此约束导致消除等式中的第二行和第五行以及第二行和第五列(19)，导致

(20)

展开后一个行列式并代入张量 $A$ 的显式形式（方程 (9)），我们得到

$det(\mathrm{vec}\{A\})$ $=$ $(1-{⟨\overline{\kappa }⟩}_{01}^1{⟨\overline{s}⟩}_{10}^0)(1-{⟨\overline{\kappa }⟩}_{12}^2{⟨\overline{s}⟩}_{21}^1)$ (21)

$-$ ${⟨\overline{\kappa }⟩}_{12}^0{⟨\overline{s}⟩}_{10}^2{⟨\overline{\kappa }⟩}_{01}^1{⟨\overline{s}⟩}_{21}^1$

在括号内，我们单独列出了每一层的临界条件。 最后一项表示包括两层在内的交互。 后者是层次包容的真正复杂性。 前两个括号内的项可以为正，这意味着没有标准的 1 层渗透，并且 $det(\mathrm{vec}\{A\})$ 可能由于最后一个交互项而变为 0。

III.7 多重超图

多重超图是另一种多层结构Sun and Bianconi (2021)。 多重超图是一组具有相同顶点集的超图$\{{\mathscr{H}}_l(V,E_l),l=1,\mathrm{\dots },L\}$。 请注意，当超图 ${\mathscr{H}}_l$ 的统计属性不同时，多重超图在统计上不等于 $\mathscr{H}(V,{\cup }_{l=1}^L{E}_l)$。 多重超图可以映射到 chygraph $\chi (V,{\cup }_{l=1}^L{E}_l,\mathrm{\Pi }=E_1\mathrm{\dots },E_L)$，其中所有边都由复合体表示，并且复合体根据它们源自的超图进行分区。 对于情况$L=2$，由于没有将复合体包含在复合体中，导致消除了(19)的第四行和第六行以及第四列和第六列，从而导致

(22)

展开后一个行列式并代入张量 $A$ 的显式形式（方程 (9)），我们得到

$det(\mathrm{vec}\{A\})$ $=$ $(1-{\overline{k}}_{01}^1{\overline{s}}_{10}^0)(1-{\overline{k}}_{02}^2{\overline{s}}_{20}^0)$ (23)

$-$ ${\overline{k}}_{01}^2{\overline{s}}_{10}^0{\overline{k}}_{02}^1{\overline{s}}_{20}^0.$

在前两个括号内，有每个超图层单独的贡献。 最后一项表示两个超图之间通过顶点的相互作用。 后者是多重超图的真正复杂性。 前两个括号内的项可以是正数，这意味着没有标准的超图渗透，并且 $det(\mathrm{vec}\{A\})$ 可能由于最后一个交互项而变为 0。

图2： 无层内夹杂物。 对于具有多重结构的随机 chygraph，对 (A) 巨型成分分数和 (B) 平均成分大小作为 $-det(\mathrm{vec}\{A\})$ 的函数进行数值估计：一层 $n_0$ 原子，两层配合物，其中$n_1$和$n_2$配合物以及$m_i$随机选择的原子包含在层$i$中随机选择的配合物中， $i=1,2$。 模型参数为$n_0=n_1=n_2={10}^6$，图例中指示的${\theta }_1=m_1^2/(n_0{n}_1)$的两个值和区间$[0.01,4]$中的${\theta }_2=m_2^2/(n_0{n}_2)$。 取 10 个随机 chygraph 的平均值。

III.8数值示例

为了测试分析结果，我将使用具有多重结构的随机 chygraph 模型：一层 $n_0$ 原子，两层具有 $n_1$ 和 $n_2$ 复合物的复合物$m_i$将随机选择的原子包含到层$i$、$i=1,2$中随机选择的络合物中。 由于内含物是随机的，所以乳糜度和成分大小均具有平均值为 $k_{0i}=m_i/n_i$ 和 $s_{i0}=m_i/n_0$、$i=1,2$ 的泊松分布。 对于泊松分布，超额平均值与平均值一致，因此 ${\overline{k}}_{0i}^j=m_i/n_i$ 和 ${\overline{s}}_{i0}^j=m_i/n_0$、$i=1,2$。 将这些值代入方程。 (23) 我们到达临界条件

$$det(\mathrm{vec}\{A\})=(1-{\theta }_1)(1-{\theta }_2)-{\theta }_1{\theta }_2$$ (24)

其中${\theta }_i=m_i^2/n_i{n}_0$,$i=1,2$。

为了测试这个表达式，我计算了巨大成分分数和不包括巨大成分的平均成分大小。 为此，我将 chygraph 投影到一个网络中，其中节点代表原子和复合物，而链接代表包含关系。 图2报告了巨型成分分数和平均成分尺寸的数值估计，作为$-det(\mathrm{vec}\{A\})$的函数，使用方程1计算。 (24)。 在$-det(\mathrm{vec}\{A\})=0$处存在相变，出现巨型组分以及排除巨型组分的平均组分尺寸的最大值。 该协议证明了多重 chygraph 分析结果的有效性。

事实证明，这个例子也是对层次包含的验证。 考虑另一个具有以下包含结构的随机 chygraph 模型：一层 $n_0$ 原子、用 $n_1$ 复合物标记为 1 的层以及随机选择的原子的 $m_1$ 包含物进入第 1 层随机选择的复合体，以及标记为 2 的层，其中带有 $n_2$ 复合体和 $m_2$ 包含从第 1 层随机选择的复合体到第 2 层随机选择的复合体。 人们可以检查方程 1 中具有分层包含的 chygraph 的临界条件。 (21) 简化为等式： (24)。 此外，将 chygraph 投影到网络中，其中节点代表原子和复合物，而链接代表包含关系，我们获得与上述多重 chygraph 相同的网络。 请注意，此等式适用于 $L=2$。 对于 $L>2$ 来说，具有多重或分层包含结构的随机 chygraph 之间不存在等价性。 对于多路复用结构，网络投影有一层连接到 $L-1$ 层。 对于分层包含，网络投影具有 $L$ 有序层，仅在相邻层之间有链接。

IV渗流理论二：有层内包裹体

我们可以从它包含的复合体（从下面）或它包含的位置（从上面）得出参考复合体。 当没有层内连接时，索引对 $ml$ 足以区分我们如何到达参考复合体。 如果我们来自下面。 如果$m>l$我们是从上面来的。 然而，当存在层内包含时，两个索引不足以指定我们如何得出案例 $m=l$ 的复合体。 我们需要一个额外的索引。 在生成函数形式主义的背景下，这意味着当从原子或另一个复合体到达复合体时，我们需要两个分量大小的生成函数。 当从下面到达时，我将用 ${\mathrm{{\rm Y}}}_{-}^{ml}$ 表示它们；当从上面到达时，我将用 ${\mathrm{{\rm Y}}}_{+}^{ml}$ 表示它们。 有了这个区别，方程。 (3) 和 (4) 被重写为

${\mathrm{\Gamma }}^l(x)$ $=$ $x{\mathrm{\Phi }}^l[{\mathrm{{\rm Y}}}_{-}^{ll}(x),\mathrm{\dots },{\mathrm{{\rm Y}}}_{-}^{lL}(x)]$ (25)

$\times$ $G^l[{\mathrm{{\rm Y}}}_{+}^{l0}(x),\mathrm{\dots },{\mathrm{{\rm Y}}}_{+}^{ll}(x)],$

${\mathrm{{\rm Y}}}_i^{ml}(x)$ $=$ $x{\mathrm{\Psi }}_i^{ml}[{\mathrm{{\rm Y}}}_{-}^{ll}(x),\mathrm{\dots },{\mathrm{{\rm Y}}}_{-}^{lL}(x)]$ (26)

$\times$ $U_i^{ml}[{\mathrm{{\rm Y}}}_{+}^{l0}(x),\mathrm{\dots },{\mathrm{{\rm Y}}}_{+}^{ll}(x)],$

其中 $l=0,\mathrm{\dots },L(\chi )$ 和

(27)

(28)

注意我隐含地假设乳糜度和复杂内超图分量是独立的。

平均多余元件尺寸${⟨\overline{\sigma }⟩}_i^{ml}={\dot{\mathrm{{\rm Y}}}}_i^{ml}(1)$、$i=-,+$可以根据等式计算。 (26)，导致

$${⟨\overline{\sigma }⟩}_i^{ml}=1+\sum _{k=l}^L{⟨\overline{\kappa }⟩}_{ilk}^m{⟨\overline{\sigma }⟩}_{-}^{lk}+\sum _{k=0}^l{⟨\overline{s}⟩}_{ilk}^m{⟨\overline{\sigma }⟩}_{+}^{lk},$$ (29)

在哪里

(30)

(31)

我们可以应用矢量化技巧来变换方程： (29) 转化为标准线性方程组。 为此，我们需要应用矢量化技巧两次。 一次像之前一样对索引 $(m,l,k)$ 进行，另一次对索引 $i$ 进行。 在第一次矢量化之后(29) 可以重写为

(32)

第二次向量化后为

$${\mathrm{vec}}^2\{A\}{\mathrm{vec}}^2\{\overline{\sigma }\}={\mathrm{vec}}^2\{B\},$$ (33)

在哪里

${(A_{--})}_{nk}^{ml}$ $=$ ${\delta }_{mn}{\delta }_{lk}-{⟨\kappa ⟩}_{nk}{\mathrm{\Theta }}_{k-n+1}{\delta }_{nl},$

${(A_{-+})}_{nk}^{ml}$ $=$ $-{⟨\overline{s}⟩}_{nk}^m{\mathrm{\Theta }}_{n-k+1}{\delta }_{nl},$

${(A_{+-})}_{nk}^{ml}$ $=$ $-{⟨\overline{\kappa }⟩}_{nk}^m{\mathrm{\Theta }}_{k-n+1}{\delta }_{nl},$

${(A_{++})}_{nk}^{ml}$ $=$ ${\delta }_{mn}{\delta }_{lk}-{⟨s⟩}_{nk}{\mathrm{\Theta }}_{n-k+1}{\delta }_{nl},$ (34)

$$B_{+}^{ml}=B_{-}^{ml}=1,$$ (35)

最后，临界条件为

$$det({\mathrm{vec}}^2\{A\})=0.$$ (36)

IV.1 $L=1$

当等式中只有一个复杂层($L(\chi )=1$)时： (34) 我们得到

(37)

人们可以继续扩大行列式。 然而，我还没有发现任何可以得出紧凑代数方程的项分组。 因此，我将直接使用方程式。 (37)。

图3： 具有层内夹杂物。 对于具有一层 $n_0$ 原子、一层 $n_1$ 复合物、随机选择的原子在随机选择的复合物中的夹杂 $m_0$ 以及随机选择的复合物在随机选择的复合物中的夹杂 $m_1$ 的随机糜晶图，作为 $-det(\mathrm{vec}\{\mathrm{vec}\{A\}\})$ 函数的 (A) 巨型分量分数和 (B) 平均分量大小的数值估计。 模型参数为 $n_0=n_1={10}^6$、图例中所示的 ${\theta }_0=m_0/n_0$ 和 $[0.1,2]$ 区间内的 ${\theta }_2=m_1/n_1$ 的两个值。 取 10 个随机 chygraph 的平均值。

IV.2数值示例

测试方程。 (37) 我考虑一个随机 chygraph，其中包含一层 $n_0$ 原子、一层 $n_1$ 复合物、$m_0$ 内含物随机选择的原子进入随机选择的复合体，并且$m_1$将随机选择的复合体包含到随机选择的复合体中。 该模型的可能实现如图1所示。

由于内含物是随机的，所以乳糜度和成分大小都具有平均值为 ${⟨\kappa ⟩}_{01}=m_0/n_0$、${⟨\kappa ⟩}_{11}=m_1/n_1$、${⟨s⟩}_{10}=m_0/n_1$ 和 ${⟨s⟩}_{11}=m_1/n_1$ 的泊松分布。 对于泊松分布，超额平均值与平均值一致，因此 ${⟨\overline{\kappa }⟩}_{01}^1=m_0/n_0$、${⟨\overline{\kappa }⟩}_{11}^1=m_1/n_1$、${⟨\overline{s}⟩}_{10}^0=m_0/n_1$ 和 ${⟨\overline{s}⟩}_{11}^1=m_1/n_1$。 将这些值代入方程。 (37) 我们得到

(38)

其中 $r=n_0/n_1$ 和 ${\theta }_i=m_i/n_i$、$i=0,1$。

图 3 报告了巨型成分分数和平均成分大小的数值估计，作为 $-det({\mathrm{vec}}^2\{A\})=0$ 的函数，使用方程 1 计算。 (38)。 在$-det({\mathrm{vec}}^2\{A\})=0$处存在相变，出现巨型组分以及排除巨型组分的平均组分尺寸的最大值。 该协议证明了具有层内连接的 chygraph 分析结果的有效性。

V结论

总之，chygraph 是一种表示复杂系统的通用组合结构。 它们允许编码不同类型的结构异质性和层次结构。 关键成分是 chygraph 的分形性质：复合体由原子和其他复合体组成，并且它也可以是其他复合体的一部分。 我已经使用矢量化计算了 chygraphs 的组件大小统计数据。 未来的工作需要将这种形式主义扩展到其他问题。