CN104133988A - 一种基于矩阵乘积的复杂网络构建方法 - Google Patents

Abstract

本发明公开了一种新的基于矩阵乘积的复杂网络构建方法，其主要步骤包括确定生成网络集合、计算生成网络邻接矩阵集合、确定生成网络度分布多项式集合、计算复杂网络的邻接矩阵、计算复杂网络的度分布等。本发明基于多个生成网络邻接矩阵的Kronecker乘积构建复杂网络，通过选取不同的生成网络、不同数量的生成网络或调整生成网络邻接矩阵Kronecker乘积的顺序可以得到不同的生成网络。采用本发明得到的复杂网络不同于通常的随机网络、小世界网络、无标度网络及自相似网络。此外，通过对其度分布多项式表达形式的次数相乘及系数相乘的运算可以从理论上严格计算出此类复杂网络的度分布。

Description

一种基于矩阵乘积的复杂网络构建方法

技术领域

本发明属于电数字数据处理领域，特别适用于特定功能的数据处理方法，具体涉及一种基于一个简单生成超网络邻接矩阵的Tracy-Singh积运算的超网络构建方法。

背景技术

复杂网络研究的深入促进了网络科学的兴起，自然科学和社会科学等许多领域的研究对象均可抽象成复杂网络进行研究。复杂网络模型的构建在社会网络、计算机网络、虚拟社会网络等其它领域的分析研究中占据极为重要的地位。上世纪中叶，ER随机网络模型的提出开创了复杂网络的系统性研究；上世纪末，WS小世界网络模型及BA无标度网络模型的相继提出开辟了复杂网络研究的新纪元。近年来，自相似特性被视为复杂网络的第三个特性而受到了人们越来越大的关注，相关学者提出了LL自相似网络模型。ER随机网络模型、WS小世界网络模型及BA无标度网络模型均是通过对初始基网络添加不同性质的节点或调整节点之间的连接而得到的，LL自相似网络模型则是基于一个生成网络通过其邻接矩阵的迭代运算得到的。现阶段复杂网络的构建主要有如下几种方法：

(1)基于图论的方法

复杂网络的研究发轫于图论，复杂网络的研究继承了图论的研究策略。参考文献[1]“Onthe evolution of random graphs”(Erdos P,Renyi A.Publ.Math.Inst.Hung.Acad.Sci.,[M].1960,5:17–60)采用完全随机的方式处理节点之间的连接，提出了ER随机网络模型——ER(Vertices,Probability)，构造出节点的度分布服从正态分布的随机网络(RandomNetwork)。参考文献[2]“Collective dynamics of'small-world'networks”(Watts D J,Strogatz S H.Nature[J].1998,393:440–442)采用随机重连处理节点之间的连接，提出了WS网络模型——WS(Vertices,Neighbors,Reprobability)，阐述了复杂网络的小世界特征，构造出节点的度分布服从指数分布的小世界网络(‘Small-world’Network)。参考文献[3]“Emergence of Scaling in Random Networks”(Barabasi A L,Albert R.Science[J].1999,286:509-512)采用增长及择优处理节点之间的连接，提出了BA网络模型——BA(InitVertices,InitProbability,AddVertices,AddLink)，阐述了复杂网络的无标度性质，构造出节点的度分布服从幂律分布的无标度网络(Scale-free Network)。小世界特性与无标度特性被誉为复杂网络的两大特性，对复杂网络的许多研究均是基于WS模型或BA模型的，采用这两种方法构建的复杂网络在统计意义上具有小世界或无标度特性。

(2)基于生成网络邻接矩阵的方法

参考文献[4]“一种复杂网络的构建方法”(李天瑞，刘胜久，珠杰，王红军[P],申请号为201410092765.2的中国专利)公开了基于一个生成网络邻接矩阵的Kronecker乘积迭代的生成一系列复杂网络，构造出同时具有自相似及小世界特性的自相似网络模型——LL(InitVertices,InitProbability,IterNum)。其自相似特性源于通过生成网络邻接矩阵的Kronecker乘积迭代产生的分形矩阵形式的邻接矩阵，而其小世界特性源于其直径不超过生成网络直径的两倍。采用此方法构建的自相似网络度分布可以从理论上严格计算得到。

(3)基于超图或超网络的方法

普通图一条边只能连接两个节点，超图中的“边”可包含多个节点。参考文献[5]“一种超网络演化模型构建及特性分析”(胡枫,赵海兴,马秀娟.中国科学:物理学力学天文学[J],2013,43:16-22)构建了一种超网络动态演化模型，从理论上分析了超度分布的特性，并进行了仿真实验，发现随着网络规模的增大，模型出现与已有的增长和优先连接复杂网络一致的结果，即复杂超网络的几种度分布显示出无标度特性。采用此方法构建的复杂网络实际上是另一种形式的无标度网络。

(4)其他方法

除传统的图论、超图及超网络方法外，其他方法也用于复杂网络的构架。参考文献[6]“基于Sierpinski分形垫的确定性复杂网络演化模型研究”(邢长明,刘方爱.物理学报[J].2010,59(3):1608-1614)基于Sierpinski分形垫，通过迭代的方式构造了小世界网络模型S-DSWN和无标度网络模型S-DSFN两个确定性增长的复杂网络模型及一个确定性的统一模型S-DUM。参考文献[7]“多种类型的网络金字塔的研究进展”(方锦清，李永，刘强.复杂系统与复杂性科学.2013.10(2):69-76)总结综述了网络模型复杂性金字塔、高科技网络金字塔及广义Farey树组织的金字塔3种类型的网络金字塔，并分析了这些金字塔的特点和性质。

总体上讲，对复杂网络特性的研究仍是现今复杂网络研究的一大热点，不可否认的是，尽管对随机网络、小世界网络、无标度网络及自相似网络均有较为成熟的理论与方法，大部分研究也与真实复杂网络相符，但仍无法全面反映现实生活中真实复杂网络的各种特点，需要进一步深入研究复杂网络的各项特性。其中，网络模型的构建是重中之重。

发明内容

为了克服现有技术的上述缺点，本发明提供了一种新的基于生成网络邻接矩阵Kronecker乘积的复杂网络构建方法，通过确定生成网络集合、计算生成网络邻接矩阵集合、确定生成网络度分布多项式集合、计算复杂网络的邻接矩阵、计算复杂网络的度分布等步骤，基于多个生成网络邻接矩阵的Kronecker乘积得到复杂网络。通过选取不同的生成网络、不同数量的生成网络或调整生成网络邻接矩阵Kronecker乘积的顺序可以得到不同的生成网络。此外，对其度分布多项式表达形式的次数相乘及系数相乘的运算可以从理论上严格计算出此类复杂网络的度分布。特别地，当只有一种生成网络时，采用该发明构建的复杂网络退化成通常的自相似网络。

一种基于矩阵乘积的复杂网络构建方法，包括如下步骤：

(1)确定生成网络集合U_G＝{G₁，G₂，G₃，…，G_i，…}；

(2)计算生成网络集合U_G中所有网络G_i的邻接矩阵A(G_i)，得到生成网络集合U_G的邻接矩阵集合U_A(G)＝{A(G₁)，A(G₂)，A(G₃)，…，A(G_i)，…}：

在生成网络集合U_G中，对于任一具有n个节点的生成网络G，其邻接矩阵A(G)是n×n的方阵，其中对于方阵中的每一个数据，若节点i与节点j相邻，则有A(G)(i,j)＝1，否则，A(G)(i,j)＝0；若生成网络G的链路数为m，则邻接矩阵A(G)中1的个数也为m，且生成网络的网络密度 $\mathrm{Density}=\frac{2m}{n(n-1)};$

(3)计算生成网络集合U_G中所有网络G_i的度分布，得到生成网络集合U_G的度分布多项式集合U_Poly(G)＝{Poly(G₁)，Poly(G₂)，Poly(G₃)，…，Poly(G_i)，…}：

在生成网络集合U_G中，对于任一具有n个节点的生成网络G，其度分布多项式表达形式Poly(G)可表述如下：

$\mathrm{Poly}\left(G\right)=\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}^{D_i}=\underset{j=1}{\overset{\∞}{\mathrm{\Σ}}}{N}_j{x}^j---\left(1\right)$

其中，n为节点数目，D_i表示第i个节点的度，N_j表示度为j的节点的数目；

(4)从生成网络集合中顺次选取k个生成网络G₍₁₎、G₍₂₎、G₍₃₎、…、G_(k－1)、G_(k)，记为G₍₁₎G₍₂₎G₍₃₎…G_(k－1)G_(k)，允许重复选取，对每个生成网络G_(i)对应的邻接矩阵A(G_(i))按如下方法计算所构建的网络的邻接矩阵A^(l)(G^(l))，其中，l代表运算的次数，A^(l)(G^(l))代表l个生成网络对应的邻接矩阵顺次进行运算后得到的一个新的复杂网络的邻接矩阵：

按照Kronecker乘积的规则进行运算，得到所构建的新的复杂网络的邻接矩阵；矩阵A(a_ij)_m×n及矩阵B(b_ij)_p×q的Kronecker乘积定义如下：

为方便书写，采用此方法顺次选取k个生成网络G₍₁₎、G₍₂₎、G₍₃₎、…、G_(k－1)、G_(k)而得到的复杂网络G^(k)可以记为G^(k)＝G₍₁₎G₍₂₎G₍₃₎…G_(k－1)G_(k)。

(5)按照如下方法计算所构建的新的复杂网络的度分布PolyDD(G^(l))，其中，l代表运算的次数，PolyDD(G^(l))代表l个生成网络对应的邻接矩阵顺次进行运算后得到的新的复杂网络度分布多项式：

$\begin{array}{c}\mathrm{PolyDD}\left(G^{(k+1)}\right)=\mathrm{DD}(\mathrm{Poly}\left(G^{\left(k\right)}\right),\mathrm{Poly}\left(G_{(k+1)}\right))\\ =\mathrm{DD}(\underset{i=1}{\overset{\∞}{\mathrm{\Σ}}}{N}_i^{\left(k\right)}{x}^i,\underset{j=1}{\overset{\∞}{\mathrm{\Σ}}}{N}_{(k+1)j}{x}^j)\\ =\underset{i=1}{\overset{\∞}{\mathrm{\Σ}}}\underset{j=1}{\overset{\∞}{\mathrm{\Σ}}}{N}_i^{\left(k\right)}{N}_{(k+1)j}{x}^{\mathrm{ij}}(i,j=\mathrm{1,2,3},...)\end{array}---\left(3\right)$

(6)重复步骤(4)及步骤(5)，得到指定节点数目、指定链路数目或指定生成网络数目的复杂网络时，终止操作。

与现有技术相比，本发明的积极效果是：

一、区别于经典的复杂网络主要通过添加节点构建复杂网络，通过邻接矩阵构建复杂网络，改进了现有复杂网络的构建方法，计算复杂度低，容易实现。

经典的复杂网络构建方法均是通过对初始基网络添加不同性质的节点或调整节点之间的连接而得到的，从生成网络的邻接矩阵出发构建复杂网络的研究较少。相较于其他复杂网络构建方法，本发明基于多个生成网络邻接矩阵的Kronecker乘积实现复杂网络模型的构建，只涉及到矩阵运算。由于矩阵运算是现有大部分软件工具(如MATLAB等)的基础，本发明提出的复杂网络构建方法借助已有的工具易于实现。

二、区别于自相似网络模型是基于一个生成网络来构建复杂网络，本发明是基于多个不同种类、不同数量、不同顺序的生成网络构建复杂网络，扩展了其应用范围。

现阶段已有的自相似复杂网络模型基于一个确定的生成网络，迭代地通过其邻接矩阵的Kronecker乘积构建复杂网络，只能生成特定节点数目或特定节点数目的复杂网络，应用范围较为狭窄。本发明基于一个包含多个生成网络的集合构建复杂网络，通过对生成网络种类、生成网络数量、生成网络顺序进行调整，并对生成网络对应的邻接矩阵进行Kronecker乘积运算可以得到不同的复杂网络，大大扩展了其应用范围。

三、区别于经典复杂网络模型的度分布是通过统计方法得到的，采用网络的度分布多项式表达形式，可以从理论上严格计算出此类复杂网络的度分布。

复杂网络的小世界特性及无标度特性均是通过统计方法得到的，无标度特性更是直接统计复杂网络的度分布得到的，但对度分布却缺少直观、形象的认识与了解。类似于化学合成中通过小分子化学物质合成大分子化学物质，采用本发明所得到的复杂网络，其度分布可通过类似于多项式乘积的方法从理论上严格计算得到，区别在于通常的多项式乘积是系数相乘次数相加，而本发明所得到的复杂网络，通过系数相乘和次数相乘可以得到复杂网络的度分布。

附图说明

本发明将通过例子并参照附图的方式说明，其中：

图1是采用文献[1]生成的随机网络模型ER(10000,0.5)度分布；

图2是采用文献[2]生成的小世界网络模型WS(10000,100,0.5)度分布；

图3是采用文献[3]生成的无标度网络模型BA(10,0.5,10000,5)度分布；

图4是采用文献[4]生成的自相似网络模型LL(6,0.4,25)度分布；

图5是包含6个生成网络G₁、G₂、G₃、G₄、G₅、G₆的拓扑结构图集合；

图6是包含6个生成网络G₁、G₂、G₃、G₄、G₅、G₆的邻接矩阵集合；

表1是包含6个生成网络G₁、G₂、G₃、G₄、G₅、G₆的度分布多项式集合；

图7是采用本发明随机生成的新的复杂网络G⁽²⁴⁾度分布，其中：

G⁽²⁴⁾＝G₁G₆G₅G₄G₅G₆G₅G₆G₄G₁G₃G₁G₁G₁G₁G₁G₄G₂G₂G₄G₆G₆G₃G₄

图8是采用本发明生成的新的复杂网络G⁽²⁴⁾度分布，其中：

G⁽²⁴⁾＝G₁G₂G₃G₄G₁G₂G₃G₄G₁G₂G₃G₄G₁G₂G₃G₄G₁G₂G₃G₄G₁G₂G₃G₄；

图9是采用本发明生成的新的复杂网络G⁽²⁴⁾度分布，其中：

G⁽²⁴⁾＝G₁G₂G₃G₄G₅G₆G₁G₂G₃G₄G₅G₆G₁G₂G₃G₄G₅G₆G₁G₂G₃G₄G₅G₆；

图10是采用本发明生成的新的复杂网络G⁽²⁴⁾度分布，其中：

G⁽²⁴⁾＝G₆G₆G₆G₆G₆G₆G₆G₆G₆G₆G₆G₆G₆G₆G₆G₆G₆G₆G₆G₆G₆G₆G₆G₆；

具体实施方式

一种基于矩阵乘积的复杂网络构建方法，包括如下步骤：

(1)确定生成网络集合U_G＝{G₁，G₂，G₃，…，G_i，…}；

(2)计算生成网络集合U_G中所有网络G_i的邻接矩阵A(G_i)，得到生成网络集合U_G的邻接矩阵集合U_A(G)＝{A(G₁)，A(G₂)，A(G₃)，…，A(G_i)，…}：

在生成网络集合U_G中，对于任一具有n个节点的生成网络G，其邻接矩阵A(G)是n×n的方阵，其中对于方阵中的每一个数据，若节点i与节点j相邻，则有A(G)(i,j)＝1，否则，A(G)(i,j)＝0；若生成网络G的链路数为m，则邻接矩阵A(G)中1的个数也为m，且生成网络的网络密度 $\mathrm{Density}=\frac{2m}{n(n-1)};$

(3)计算生成网络集合U_G中所有网络G_i的度分布，得到生成网络集合U_G的度分布多项式集合U_Poly(G)＝{Poly(G₁)，Poly(G₂)，Poly(G₃)，…，Poly(G_i)，…}：

在生成网络集合U_G中，对于任一具有n个节点的生成网络G，其度分布多项式表达形式Poly(G)可表述如下：

$\mathrm{Poly}\left(G\right)=\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}^{D_i}=\underset{j=1}{\overset{\∞}{\mathrm{\Σ}}}{N}_j{x}^j---\left(1\right)$

其中，n为节点数目，D_i表示第i个节点的度，N_j表示度为j的节点的数目；

(4)从生成网络集合中顺次选取k个生成网络G₍₁₎、G₍₂₎、G₍₃₎、…、G_(k－1)、G_(k)，记为G₍₁₎G₍₂₎G₍₃₎…G_(k－1)G_(k)，允许重复选取，对每个生成网络G_(i)对应的邻接矩阵A(G_(i))按如下方法计算所构建的网络的邻接矩阵A^(l)(G^(l))，其中，l代表运算的次数，A^(l)(G^(l))代表l个生成网络对应的邻接矩阵顺次进行运算后得到的一个新的复杂网络的邻接矩阵：

按照Kronecker乘积的规则进行运算，得到所构建的新的复杂网络的邻接矩阵；矩阵A(a_ij)_m×n及矩阵B(b_ij)_p×q的Kronecker乘积定义如下：

为方便书写，采用此方法顺次选取k个生成网络G₍₁₎、G₍₂₎、G₍₃₎、…、G_(k－1)、G_(k)而得到的复杂网络G^(k)可以记为G^(k)＝G₍₁₎G₍₂₎G₍₃₎…G_(k－1)G_(k)。

(5)按照如下方法计算所构建的新的复杂网络的度分布PolyDD(G^(l))，其中，l代表运算的次数，PolyDD(G^(l))代表l个生成网络对应的邻接矩阵顺次进行运算后得到的新的复杂网络度分布多项式：

$\begin{array}{c}\mathrm{PolyDD}\left(G^{(k+1)}\right)=\mathrm{DD}(\mathrm{Poly}\left(G^{\left(k\right)}\right),\mathrm{Poly}\left(G_{(k+1)}\right))\\ =\mathrm{DD}(\underset{i=1}{\overset{\∞}{\mathrm{\Σ}}}{N}_i^{\left(k\right)}{x}^i,\underset{j=1}{\overset{\∞}{\mathrm{\Σ}}}{N}_{(k+1)j}{x}^j)\\ =\underset{i=1}{\overset{\∞}{\mathrm{\Σ}}}\underset{j=1}{\overset{\∞}{\mathrm{\Σ}}}{N}_i^{\left(k\right)}{N}_{(k+1)j}{x}^{\mathrm{ij}}(i,j=\mathrm{1,2,3},...)\end{array}---\left(3\right)$

(6)重复步骤(4)及步骤(5)，得到指定节点数目、指定链路数目或指定生成网络数目的复杂网络时，终止操作。

对于顺次选取k个生成网络G₍₁₎、G₍₂₎、G₍₃₎、…、G_(k－1)、G_(k)而得到的复杂网络G^(k)，得到的复杂网络的节点数 $n\left(k\right)=\underset{i=1}{\overset{k}{\mathrm{\Π}}}n\left(G_{\left(i\right)}\right),$ 链路数 $n\left(k\right)=2^{k-1}\underset{i=1}{\overset{k}{\mathrm{\Π}}}m\left(G_{\left(i\right)}\right),$ 网络密度其中n(G_(i))表示生成网络的节点数，m(G_(i))生成网络的链路数。

采用本发明提出的新的复杂网络构建方法，基于一个包含多个生成网络的集合，通过顺次选取的一系列生成网络的邻接矩阵的Kronecker乘积运算构建复杂网络。通过调整生成网络种类、生成网络数量及生成网络顺序可以得到不同的复杂网络，区别于经典的随机网络模型、小世界网络模型及无标度网络模型，也不同于基于一个生成网络而构建的自相似网络模型。通过对生成网络度分布多项式表达形式的次数相乘及系数相乘的运算可以从理论上严格计算出此类复杂网络的度分布。在极端情况下，当只有一种生成网络，即生成网络集合基数为1时，采用本发明得到的复杂网络退化为通常的自相似网络模型。

仿真实验

为了验证本发明基于一个包含多个生成网络的集合，并通过顺次选取的生成网络对应的邻接矩阵的Kronecker乘积而得到的复杂网络构建方法的有效性，进行了仿真实验。

确定一个包含6个生成网络的集合U_G，其拓扑结构图如图5所示，按从上至下，从左到右的顺序依次为G₁、G₂、G₃、G₄、G₅、G₆；其对应的邻接矩阵集合U_A(G)如图6所示，按从上至下，从左到右的顺序依次为A(G₁)、A(G₂)、A(G₃)、A(G₄)、A(G₅)、A(G₆)，其对应的度分布多项式表达形式集合U_Poly(G)如表1所示。

对图G₁而言，图G₁共有4个节点，其中，第1个节点的度为3，第2个节点的度为2，第3个节点的度为2，第4个节点的度为1，故其度分布多项式表达形式为：Poly(G)＝x+2x²+x³。

对图G₂而言，图G₂共有4个节点，其中，第1个节点的度为3，第2个节点的度为3，第3个节点的度为2，第4个节点的度为2，故其度分布多项式表达形式为：Poly(G)＝2x²+2x³。

对图G₃而言，图G₃共有5个节点，其中，第1个节点的度为3，第2个节点的度为1，第3个节点的度为2，第4个节点的度为2，第5个节点的度为2，故其度分布多项式表达形式为：Poly(G)＝x+3x²+x³。

对图G₄而言，图G₄共有5个节点，其中，第1个节点的度为2，第2个节点的度为3，第3个节点的度为3，第4个节点的度为1，第5个节点的度为1，故其度分布多项式表达形式为：Poly(G)＝2x+x²+2x³。

对图G₅而言，图G₅共有6个节点，其中，第1个节点的度为3，第2个节点的度为1，第3个节点的度为2，第4个节点的度为2，第5个节点的度为3，第6个节点的度为1，故其度分布多项式表达形式为：Poly(G)＝2x+2x²+2x³。

对图G₆而言，图G₆共有6个节点，其中，第1个节点的度为3，第2个节点的度为2，第3个节点的度为2，第4个节点的度为2，第5个节点的度为2，第6个节点的度为1，故其度分布多项式表达形式为：Poly(G)＝x+4x²+x³。

表1度分布多项式表达形式集合U_Poly(G)

实验一

在生成网络集合U_G中随机选取24个生成网络，允许重复，采用本发明构建得到的复杂网络记为G⁽²⁴⁾，其中：

G⁽²⁴⁾＝G₁G₆G₅G₄G₅G₆G₅G₆G₄G₁G₃G₁G₁G₁G₁G₁G₄G₂G₂G₄G₆G₆G₃G₄；

其度分布如图7所示。

实验二

在生成网络集合U_G中循环6遍顺次选取生成网络G₁、G₂、G₃、G₄，采用本发明构建得到的复杂网络记为G⁽²⁴⁾，其中：

G⁽²⁴⁾＝G₁G₂G₃G₄G₁G₂G₃G₄G₁G₂G₃G₄G₁G₂G₃G₄G₁G₂G₃G₄G₁G₂G₃G₄；

其度分布如图8所示。

实验三

在生成网络集合U_G中循环4遍顺次选取生成网络G₁、G₂、G₃、G₄、G₅、G₆，采用本发明构建得到的复杂网络记为G⁽²⁴⁾，其中：

G⁽²⁴⁾＝G₁G₂G₃G₄G₅G₆G₁G₂G₃G₄G₅G₆G₁G₂G₃G₄G₅G₆G₁G₂G₃G₄G₅G₆；

其度分布如图9所示。

实验四

在生成网络集合U_G中循环24遍顺次选取生成网络G₆，采用本发明构建得到的复杂网络记为G⁽²⁴⁾，其中：

G⁽²⁴⁾＝G₆G₆G₆G₆G₆G₆G₆G₆G₆G₆G₆G₆G₆G₆G₆G₆G₆G₆G₆G₆G₆G₆G₆G₆；

其度分布如图10所示。

从图7、图8、图9、图10可以看出，采用发明得到的复杂网络度分布不同于经典的ER随机网络模型(如图1等)、WS小世界网络模型(如图2等)、BA无标度网络模型(如图3等)，而是与通常的自相似网络模型相似(如图4等)，主要原因在于与基于一个生成网络而构建复杂网络的自相似网络模型类似，均是基于邻接矩阵的Kronecker乘积而得到复杂网络的。从图7、图8、图9、图10还可以看出，采用本发明得到的复杂网络，通过调整生成网络种类、生成网络数量及生成网络顺序可以得到不同的复杂网络。对比图4及图10，可以看出，在极端情况下，当只有一种生成网络，即生成网络集合基数为1时，采用本发明得到的复杂网络退化为通常的自相似网络模型。

Claims (3)

1.一种基于矩阵乘积的复杂网络构建方法，其特征在于：包括如下步骤：

(1)确定生成网络集合U_G＝{G₁，G₂，G₃，…，G_i，…}；

(2)计算生成网络集合U_G中所有网络G_i的邻接矩阵A(G_i)，得到生成网络集合U_G的邻接矩阵集合U_A(G)＝{A(G₁)，A(G₂)，A(G₃)，…，A(G_i)，…}：

在生成网络集合U_G中，对于任一具有n个节点的生成网络G，其邻接矩阵A(G)是n×n的方阵，其中对于方阵中的每一个数据，若节点i与节点j相邻，则有A(G)(i,j)＝1，否则，A(G)(i,j)＝0；若生成网络G的链路数为m，则邻接矩阵A(G)中1的个数也为m，且生成网络的网络密度 $\mathrm{Density}=\frac{2m}{n(n-1)};$

(3)计算生成网络集合U_G中所有网络G_i的度分布，得到生成网络集合U_G的度分布多项式集合U_Poly(G)＝{Poly(G₁)，Poly(G₂)，Poly(G₃)，…，Poly(G_i)，…}：

在生成网络集合U_G中，对于任一具有n个节点的生成网络G，其度分布多项式表达形式Poly(G)可表述如下：

$\mathrm{Poly}\left(G\right)=\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{x}^{D_i}=\underset{j=1}{\overset{\∞}{\mathrm{\Σ}}}{N}_j{x}^j;$

其中，n为节点数目，D_i表示第i个节点的度，N_j表示度为j的节点的数目；

(4)从生成网络集合中顺次选取k个生成网络G₍₁₎、G₍₂₎、G₍₃₎、…、G_(k－1)、G_(k)，记为G₍₁₎G₍₂₎G₍₃₎…G_(k－1)G_(k)，允许重复选取，对每个生成网络G_(i)对应的邻接矩阵A(G_(i))按如下方法计算所构建的网络的邻接矩阵A^(l)(G^(l))，其中，l代表运算的次数，A^(l)(G^(l))代表l个生成网络对应的邻接矩阵顺次进行运算后得到的一个新的复杂网络的邻接矩阵：

按照Kronecker乘积的规则进行运算，得到所构建的新的复杂网络的邻接矩阵；矩阵A(a_ij)_m×n及矩阵B(b_ij)_p×q的Kronecker乘积定义如下：

采用此方法顺次选取k个生成网络G₍₁₎、G₍₂₎、G₍₃₎、…、G_(k－1)、G_(k)得到新的复杂网络G^(k)；

(5)按照如下方法计算所构建的新的复杂网络G^(k)的度分布PolyDD(G^(l))，其中，l代表运算的次数，PolyDD(G^(l))代表l个生成网络对应的邻接矩阵顺次进行运算后得到的新的复杂网络G^(k)度分布多项式：

$\begin{array}{c}\mathrm{PolyDD}\left(G^{(k+1)}\right)=\mathrm{DD}(\mathrm{Poly}\left(G^{\left(k\right)}\right),\mathrm{Poly}\left(G_{(k+1)}\right))\\ =\mathrm{DD}(\underset{i=1}{\overset{\∞}{\mathrm{\Σ}}}{N}_i^{\left(k\right)}{x}^i,\underset{j=1}{\overset{\∞}{\mathrm{\Σ}}}{N}_{(k+1)j}{x}^j)\\ =\underset{i=1}{\overset{\∞}{\mathrm{\Σ}}}\underset{j=1}{\overset{\∞}{\mathrm{\Σ}}}{N}_i^{\left(k\right)}{N}_{(k+1)j}{x}^{\mathrm{ij}}(i,j=\mathrm{1,2,3},...)\end{array};$

(6)重复步骤(4)及步骤(5)，得到指定节点数目、指定链路数目或指定生成网络数目的复杂网络时，终止操作。

2.根据权利要求1所述的一种基于矩阵乘积的复杂网络构建方法，其特征在于：确定生成网络集合U_G时，选择节点数目n小于等于10且节点间连接较少的简单网络作为生成网络。

3.根据权利要求1所述的一种基于矩阵乘积的复杂网络构建方法，其特征在于：对于顺次选取k个生成网络G₍₁₎、G₍₂₎、G₍₃₎、…、G_(k－1)、G_(k)而得到的复杂网络G^(k)，得到的复杂网络的节点数 $n\left(k\right)=\underset{i=1}{\overset{k}{\mathrm{\Π}}}n\left(G_{\left(i\right)}\right),$ 链路数 $n\left(k\right)=2^{k-1}\underset{i=1}{\overset{k}{\mathrm{\Π}}}m\left(G_{\left(i\right)}\right),$ 网络密度其中n(G_(i))表示生成网络的节点数，m(G_(i))生成网络的链路数。