摘要： 鑒于現(xiàn)有超網(wǎng)絡(luò)魯棒性的研究中并未考慮超邊內(nèi)部結(jié)構(gòu)對(duì)超網(wǎng)絡(luò)魯棒性的影響，提出一種能描述超邊內(nèi)部結(jié)構(gòu)與超網(wǎng)絡(luò)魯棒性關(guān)系的容量負(fù)載模型，并通過仿真實(shí)驗(yàn)獲得了超邊內(nèi)節(jié)點(diǎn)在優(yōu)先連接、隨機(jī)連接、全連接3種方式下，k均勻無標(biāo)度超網(wǎng)絡(luò)的魯棒性。通過對(duì)比分析發(fā)現(xiàn)，無標(biāo)度超網(wǎng)絡(luò)的魯棒性與其超邊內(nèi)節(jié)點(diǎn)的連接方式密切相關(guān)，同時(shí)也與超邊內(nèi)節(jié)點(diǎn)的規(guī)模k以及超邊內(nèi)部普通邊的數(shù)量mk有關(guān)。研究結(jié)果表明，超邊內(nèi)部結(jié)構(gòu)對(duì)無標(biāo)度超網(wǎng)絡(luò)整體的魯棒性有較大的影響。

關(guān)鍵詞： 無標(biāo)度超網(wǎng)絡(luò)；超邊內(nèi)部結(jié)構(gòu)；容量負(fù)載模型；級(jí)聯(lián)故障；魯棒性

中圖分類號(hào)： O157.5; N945.15文獻(xiàn)標(biāo)識(shí)碼： A

Influence of Structure Inside Hyperedge on Robustness of Scale-free Hypernetwork

ZHOU Bin， MA Fuxiang GAO Shujie， MA Xiujuan LI Mingjie

（1.College of Computer， Qinghai Normal University， Xining 810016， China;

2.The State Key Laboratory of Tibetan Intelligent Information Processing and Application， Xining 810008， China）

Abstract：In the existing work on the robustness of the hypernetworks， researchers have not considered the effect of internal structure on the robustness of the hypernetworks. Aiming at this problem， this paper proposes a capacity-load model that can describe the relationship between the internal structure and the robustness of the hypernetworks. By simulation experiments， we obtain the robustness of the k-uniform scale-free hypernetwork under three modes： preferential connection， random connection， and completed connection within hyperedges. Analysis of the comparison experiments reveals that the robustness of the scale-free hypernetwork is related to the ways of nodes' connection， the size k of the nodes， and the number of ordinary edges mk within the hyperedges.The results show that the internal structure of the hyperedges has a large impact on the overall robustness of the scale-free hypernetworks.

Keywords： scale-free hypernetwork; internal structure of hyperedge; capacity-load model; cascading failures; robustness

0 引 言

隨著信息技術(shù)的發(fā)展，復(fù)雜網(wǎng)絡(luò)已經(jīng)成為建模各類復(fù)雜系統(tǒng)的有效工具，例如電力網(wǎng)絡(luò)[1]、通信網(wǎng)絡(luò)[2]、交通網(wǎng)絡(luò)[3]和金融網(wǎng)絡(luò)[4]等。建模復(fù)雜系統(tǒng)的最終目標(biāo)是分析系統(tǒng)的性能，并達(dá)到改善系統(tǒng)性能的目的。對(duì)很多現(xiàn)實(shí)復(fù)雜系統(tǒng)來說，魯棒性是其最重要、最基本的系統(tǒng)性能之一。近年來，研究者也依據(jù)復(fù)雜網(wǎng)絡(luò)理論成功研究了各種復(fù)雜系統(tǒng)的魯棒性，獲得了很多重要的結(jié)論，并提出了眾多行之有效的優(yōu)化系統(tǒng)魯棒性的策略[57]。

但隨著人類社會(huì)的不斷發(fā)展，與人類生活息息相關(guān)的各類系統(tǒng)也越來越復(fù)雜。國內(nèi)外多位學(xué)者[810]提出很多現(xiàn)有的系統(tǒng)已不能單純地抽象為任意兩個(gè)節(jié)點(diǎn)之間的關(guān)系，而是多個(gè)節(jié)點(diǎn)間存在更加復(fù)雜的聯(lián)系。因此，基于圖的復(fù)雜網(wǎng)絡(luò)理論已不能完全表示日益復(fù)雜的系統(tǒng)結(jié)構(gòu)[1113]，也無法深刻分析這類復(fù)雜系統(tǒng)的性能。超網(wǎng)絡(luò)理論的出現(xiàn)，為研究這類復(fù)雜系統(tǒng)帶來了新的思路和研究方法。超網(wǎng)絡(luò)中的超邊可以更好地表示多個(gè)節(jié)點(diǎn)間同時(shí)存在的某種復(fù)雜關(guān)系，因此近年來，超網(wǎng)絡(luò)被用于建模許多真實(shí)復(fù)雜系統(tǒng)，且表現(xiàn)出了比復(fù)雜網(wǎng)絡(luò)更好的優(yōu)越性[1416]。雖然超網(wǎng)絡(luò)的建模研究已日趨成熟，但由于超網(wǎng)絡(luò)結(jié)構(gòu)的復(fù)雜性，目前對(duì)超網(wǎng)絡(luò)動(dòng)力學(xué)性質(zhì)的研究還處在起步階段。同時(shí)，已有的研究結(jié)果也表明，超網(wǎng)絡(luò)確實(shí)表現(xiàn)出了和普通復(fù)雜網(wǎng)絡(luò)不同的動(dòng)力學(xué)性質(zhì)。例如，Wu等[17]將復(fù)雜系統(tǒng)的同步首次引入到超網(wǎng)絡(luò)中，發(fā)現(xiàn)超網(wǎng)絡(luò)相較于普通網(wǎng)絡(luò)可以更恰當(dāng)?shù)孛枋龆喾矫娴鸟詈蠀f(xié)作關(guān)系，更有利于同步問題的研究。Gong等[18]構(gòu)建了在線社交超網(wǎng)絡(luò)，解決了普通社交網(wǎng)絡(luò)不能準(zhǔn)確描述社交網(wǎng)絡(luò)中群聚特性的問題，并對(duì)在線社交超網(wǎng)絡(luò)中信息全局傳播的動(dòng)態(tài)過程進(jìn)行了仿真分析，發(fā)現(xiàn)同復(fù)雜網(wǎng)絡(luò)結(jié)構(gòu)相比，基于超圖結(jié)構(gòu)的超網(wǎng)絡(luò)中信息全局傳播速度更快，波及范圍更廣。Ma等[19]通過對(duì)比研究發(fā)現(xiàn)，雖然超網(wǎng)絡(luò)面對(duì)外部攻擊表現(xiàn)出了和普通網(wǎng)絡(luò)類似的性質(zhì)，但超網(wǎng)絡(luò)的級(jí)聯(lián)故障過程卻比普通的復(fù)雜網(wǎng)絡(luò)緩慢，面對(duì)同樣的外部擾動(dòng)，表現(xiàn)出了比普通網(wǎng)絡(luò)更強(qiáng)的魯棒性。Ma等[20]同樣根據(jù)超網(wǎng)絡(luò)的結(jié)構(gòu)提出了一類可以描述超網(wǎng)絡(luò)級(jí)聯(lián)故障過程的級(jí)聯(lián)故障模型，通過該模型的仿真分析，也獲得了類似的結(jié)論。此外，Chen等[21]提出了基于超邊擴(kuò)散的容量負(fù)載級(jí)聯(lián)故障模型，將模型分別應(yīng)用在隨機(jī)超網(wǎng)絡(luò)和小世界超網(wǎng)絡(luò)上，對(duì)比發(fā)現(xiàn)兩類超網(wǎng)絡(luò)相較于普通網(wǎng)絡(luò)具有更強(qiáng)的魯棒性，且隨機(jī)超網(wǎng)絡(luò)比小世界超網(wǎng)絡(luò)具有更強(qiáng)的魯棒性。Luo等[22]發(fā)現(xiàn)超網(wǎng)絡(luò)相較于普通網(wǎng)絡(luò)更適合用于表示復(fù)雜的城市公交網(wǎng)絡(luò)，節(jié)點(diǎn)為公交站點(diǎn)，超邊為線路，并分析了該超網(wǎng)絡(luò)的魯棒性，發(fā)現(xiàn)公交超網(wǎng)絡(luò)對(duì)隨機(jī)攻擊魯棒，蓄意攻擊脆弱。但上述研究結(jié)果大多都默認(rèn)超網(wǎng)絡(luò)超邊內(nèi)節(jié)點(diǎn)的連接方式為全連接，然而現(xiàn)實(shí)超網(wǎng)絡(luò)中節(jié)點(diǎn)內(nèi)部的連接方式往往都是非全連接的，例如在微信群超網(wǎng)絡(luò)，交通超網(wǎng)絡(luò)等。因此，以往的研究忽略了超邊的內(nèi)部結(jié)構(gòu)對(duì)超網(wǎng)絡(luò)魯棒性的影響，也即忽略了微觀結(jié)構(gòu)對(duì)系統(tǒng)整體性能產(chǎn)生的作用。例如，在集成電路開發(fā)中，通過超網(wǎng)絡(luò)的思想去設(shè)計(jì)集成電路，研究集成功能塊內(nèi)部電子元件間的連接關(guān)系，可以進(jìn)一步優(yōu)化集成電路的魯棒性；再比如，在快遞服務(wù)網(wǎng)的布局當(dāng)中，通過布局快遞網(wǎng)內(nèi)部的交互關(guān)系，從而對(duì)快遞服務(wù)網(wǎng)的結(jié)構(gòu)進(jìn)行優(yōu)化，增強(qiáng)其魯棒性。因此，研究超邊內(nèi)部結(jié)構(gòu)對(duì)超網(wǎng)絡(luò)魯棒性的影響，將有利于更好地發(fā)現(xiàn)超網(wǎng)絡(luò)結(jié)構(gòu)與超網(wǎng)絡(luò)魯棒性之間的關(guān)系，同時(shí)能更全面地認(rèn)識(shí)和理解超網(wǎng)絡(luò)魯棒性的影響因素，從而能提出更好的優(yōu)化策略，提高超網(wǎng)絡(luò)抵抗各類攻擊的能力，為進(jìn)一步增強(qiáng)復(fù)雜系統(tǒng)的魯棒性提供可靠的理論支撐。

為此，本文以Barabasi-Albert （BA）無標(biāo)度超網(wǎng)絡(luò)為基礎(chǔ)，考慮了超邊內(nèi)節(jié)點(diǎn)的連接方式，構(gòu)建了三類k均勻無標(biāo)度超網(wǎng)絡(luò)，并基于Motter等提出的容量負(fù)載模型[23]的思想提出了一種新的容量負(fù)載模型，研究了超邊內(nèi)部結(jié)構(gòu)對(duì)超網(wǎng)絡(luò)魯棒性的影響。

1 超網(wǎng)絡(luò)與超網(wǎng)絡(luò)模型

超網(wǎng)絡(luò)是以超圖[2425]作為拓?fù)浣Y(jié)構(gòu)對(duì)復(fù)雜系統(tǒng)進(jìn)行抽象后形成的網(wǎng)絡(luò)，超網(wǎng)絡(luò)通過超邊有效表達(dá)了多個(gè)節(jié)點(diǎn)同時(shí)滿足的某種關(guān)系，比普通網(wǎng)絡(luò)更加直觀。例如，作者合作超網(wǎng)絡(luò)表達(dá)了多個(gè)作者合作同一篇論文的情況[15]；交通超網(wǎng)絡(luò)表達(dá)了多個(gè)交通站點(diǎn)被同一條交通線路經(jīng)過的情況[22]。

超網(wǎng)絡(luò)中節(jié)點(diǎn)vi的超度定義為包含節(jié)點(diǎn)vi的超邊的個(gè)數(shù)，記為dH（vi）[26]。而節(jié)點(diǎn)vi的節(jié)點(diǎn)度與普通網(wǎng)絡(luò)類似，依然定義為節(jié)點(diǎn)vi所關(guān)聯(lián)的普通邊的個(gè)數(shù)，記為d（vi）[26]。節(jié)點(diǎn)vi在單獨(dú)一條超邊ei內(nèi)的普通度記為kei（vi）。

為了更好地發(fā)現(xiàn)超網(wǎng)絡(luò)中超邊的內(nèi)部結(jié)構(gòu)對(duì)超網(wǎng)絡(luò)魯棒性的影響，本文參考了Hu等[27]提出的k均勻無標(biāo)度超網(wǎng)絡(luò)的演化思想，構(gòu)建了三類超邊內(nèi)節(jié)點(diǎn)不同連接方式的無標(biāo)度超網(wǎng)絡(luò)模型，如表1所示。三類超網(wǎng)絡(luò)模型的具體構(gòu)建過程如下所述。

步驟1 分別構(gòu)建節(jié)點(diǎn)總數(shù)為N，超邊總數(shù)為M的k均勻無標(biāo)度超網(wǎng)絡(luò)；

步驟2 斷開超邊內(nèi)節(jié)點(diǎn)間的原有連接。當(dāng)超邊內(nèi)的節(jié)點(diǎn)之間采用優(yōu)先連接的方式時(shí)，轉(zhuǎn)1）執(zhí)行；采用隨機(jī)連接方式時(shí)，轉(zhuǎn)2）執(zhí)行；采用全連接方式時(shí)，轉(zhuǎn)3）執(zhí)行；

1）超邊內(nèi)部節(jié)點(diǎn)采用優(yōu)先連接方式：

（1）初始化：設(shè)Vei表示超邊ei內(nèi)的節(jié)點(diǎn)集合，超邊ei（1≤i≤M）內(nèi)的k（1≤k≤N）個(gè)節(jié)點(diǎn)構(gòu)成超邊ei的導(dǎo)出子圖G（ei）中的孤立節(jié)點(diǎn)；選定3個(gè)節(jié)點(diǎn)進(jìn)行全連接使其成為超邊內(nèi)初始的連通分支C0；設(shè)Γt表示t時(shí)刻超邊內(nèi)連通分支Ct的節(jié)點(diǎn)集合；

（2）增長：t+1時(shí)刻隨機(jī)選擇連通分支Ct外的一個(gè)節(jié)點(diǎn)vk（vk∈Vei，vkCt），增加到Ct中，構(gòu)成新的連通分支Ct+1；vk加入連通分支時(shí)，以（3）中的方式選擇連通分支Ct內(nèi)的一個(gè)節(jié)點(diǎn)vi（vi∈Vei，vi∈Ct）；重復(fù)（2），不斷增大連通分支，直至其包含超邊ei內(nèi)的所有節(jié)點(diǎn)；

（3）優(yōu)先連接：設(shè)kei，Ct（vi）表示超邊ei內(nèi)的節(jié)點(diǎn)vi在連通分支Ct中的普通度。連通分支Ct外的節(jié)點(diǎn)vk選擇連通分支內(nèi)的δm個(gè)節(jié)點(diǎn)并用普通邊進(jìn)行連接；連通分支內(nèi)的節(jié)點(diǎn)vi以優(yōu)先連接概率Πi被選擇。

Πi=kei（vi）∑Γtj=1kei（vj），vi，vj∈Ct（1）

2）超邊內(nèi)部節(jié)點(diǎn)采用隨機(jī)連接方式：

（1）初始化：使超邊ei（1≤i≤M）內(nèi)的k（1≤k≤N）個(gè)節(jié)點(diǎn)成為超邊ei的導(dǎo)出子圖G（ei）中的孤立節(jié)點(diǎn)，然后給定節(jié)點(diǎn)重連邊的概率p∈［0，1］；

（2）超邊內(nèi)節(jié)點(diǎn)以如下的方式隨機(jī)重連：a.生成一個(gè)隨機(jī)數(shù)r；b.隨機(jī)選擇同一超邊內(nèi)的一對(duì)節(jié)點(diǎn)；c.如果r<p，判斷 b 中所選擇的節(jié)點(diǎn)對(duì)之間是否有連邊，若無連邊，則用一條普通邊連接，否則不做任何操作；

（3）重復(fù)步驟（2），直至生成p×k（k-1）/2條普通邊。

3）超邊內(nèi)部節(jié)點(diǎn)全連接方式：

（1）初始化：使超邊ei（1≤i≤M）內(nèi)的k（1≤k≤N）個(gè)節(jié)點(diǎn)成為超邊ei的導(dǎo)出子圖G（ei）中的孤立節(jié)點(diǎn)；

（2）超邊內(nèi)節(jié)點(diǎn)重連：使k個(gè)節(jié)點(diǎn)中每對(duì)不同的節(jié)點(diǎn)之間都恰有一條邊相連，生成包含k個(gè)節(jié)點(diǎn)、k（k-1）/2條邊的全連通子圖；

步驟3 重復(fù)步驟2，直至遍歷超網(wǎng)絡(luò)中的所有超邊。

2 基于超邊內(nèi)部結(jié)構(gòu)的超網(wǎng)絡(luò)容量負(fù)載級(jí)聯(lián)故障模型

以BA_C超網(wǎng)絡(luò)為例，當(dāng)超網(wǎng)絡(luò)中一個(gè)節(jié)點(diǎn)發(fā)生崩潰后，會(huì)將自身的負(fù)載分配到其未失效的鄰居節(jié)點(diǎn)，整個(gè)超網(wǎng)絡(luò)的級(jí)聯(lián)故障進(jìn)程如圖1所示。圖1a中節(jié)點(diǎn)v4（填充狀態(tài)為失效節(jié)點(diǎn)，否則為正常節(jié)點(diǎn)）是初始時(shí)刻受到攻擊而失效的節(jié)點(diǎn)，其余為正常節(jié)點(diǎn)；圖1b中顯示初始攻擊節(jié)點(diǎn)v4失效后，向其鄰居節(jié)點(diǎn)分配負(fù)載（箭頭表示為失效節(jié)點(diǎn)負(fù)載傳播的方向）；圖1c中顯示第一輪負(fù)載分配結(jié)束后，節(jié)點(diǎn)的失效情況。此時(shí)超邊e e2和e3中的所有節(jié)點(diǎn)均已失效，則這三條超邊失效（虛線表示超邊失效，實(shí)線表示超邊未失效），負(fù)載通過節(jié)點(diǎn)v7和v8繼續(xù)向未失效的超邊e4傳播；圖1d表示所有超邊均已失效，級(jí)聯(lián)故障的擴(kuò)散到此結(jié)束。

以往研究者提出的超網(wǎng)絡(luò)級(jí)聯(lián)故障模型，存在考慮指標(biāo)過于單一的問題，以Chen等[21]提出的基于超邊擴(kuò)散的超網(wǎng)絡(luò)級(jí)聯(lián)故障模型為例，作者將超網(wǎng)絡(luò)轉(zhuǎn)換成線圖進(jìn)行研究，默認(rèn)超邊內(nèi)部節(jié)點(diǎn)全連接，當(dāng)一條超邊受到攻擊失效時(shí)，負(fù)載會(huì)傳播到與其相鄰的未失效超邊上，并沒有考慮負(fù)載在超邊內(nèi)部通過節(jié)點(diǎn)傳播的情況。為更好地描述超邊內(nèi)部結(jié)構(gòu)對(duì)超網(wǎng)絡(luò)魯棒性的影響，本文提出了一種新的容量負(fù)載模型用于分析超邊內(nèi)部結(jié)構(gòu)與超網(wǎng)絡(luò)魯棒性之間的關(guān)系。

在一個(gè)有N個(gè)節(jié)點(diǎn)的超網(wǎng)絡(luò)中，節(jié)點(diǎn)vi的初始負(fù)載與該節(jié)點(diǎn)的超度dH（vi）和節(jié)點(diǎn)度d（vi）有關(guān)，其初始負(fù)載Lvi（0）被定義為

Lvi（0）=α（dH（vi）+d（vi））β i=，…，N，α≥ β≥1（2）

為了控制節(jié)點(diǎn)vi的初始負(fù)載，設(shè)α為負(fù)載參數(shù)，β為可調(diào)參數(shù)。由式（2）可知節(jié)點(diǎn)初始負(fù)載與其節(jié)點(diǎn)超度和節(jié)點(diǎn)度成正比。

設(shè)節(jié)點(diǎn)vi與節(jié)點(diǎn)vj在某條超邊內(nèi)通過普通邊連接。在t時(shí)刻，當(dāng)節(jié)點(diǎn)vi因故障失效時(shí)，首先其負(fù)載Lvi將按照失效節(jié)點(diǎn)vi的超度dH（vi）等分，并將等分后的負(fù)載分配給與其關(guān)聯(lián)的所有未失效超邊，則與vi關(guān)聯(lián)的所有未失效超邊接收到的負(fù)載都為Lvi/dH（vi）。然后將某條超邊ei接收到的負(fù)載按照節(jié)點(diǎn)vi在超邊ei中的節(jié)點(diǎn)度kei（vi）等分，則節(jié)點(diǎn)vj接收到的負(fù)載如式（3）所示。

ΔLvi→vjt=Lvit×1dθHvi×1kθeivi（3）

其中，θ為擾動(dòng)參數(shù)，且0≤θ≤1。由式（3）可以看出，t時(shí)刻超邊ei內(nèi)的節(jié)點(diǎn)vj接收到的額外載荷ΔLvi→vj（t）與t時(shí)刻失效節(jié)點(diǎn)vi的超度、節(jié)點(diǎn)vi在超邊ei內(nèi)的普通度以及節(jié)點(diǎn)vi的初始負(fù)載有關(guān)。若t時(shí)刻節(jié)點(diǎn)vj未失效，則其負(fù)載為節(jié)點(diǎn)vj在t-1時(shí)刻的負(fù)載加上從失效節(jié)點(diǎn)vi重分配到的負(fù)載，即：

Lvj（t）=Lvj（t-1）+ΔLvi→vj（t）（4）

在現(xiàn)實(shí)超網(wǎng)絡(luò)中，容量是節(jié)點(diǎn)或者超邊所能處理負(fù)載的最大值，與節(jié)點(diǎn)的初始負(fù)載成正比。設(shè)節(jié)點(diǎn)vj的容量為Cvj，如式（5）所示。

Cvj=（1+T）Lvj T≥0（5）

其中，T為容量參數(shù)，T的值越大，節(jié)點(diǎn)的容量就越大，抵御故障的能力就越強(qiáng)，但抵御成本會(huì)增加。臨界閾值TC是避免超網(wǎng)絡(luò)發(fā)生全局崩潰的最小容量值。當(dāng)T>TC時(shí)，整個(gè)超網(wǎng)絡(luò)不會(huì)出現(xiàn)全局崩潰。當(dāng)T<TC時(shí)，整個(gè)超網(wǎng)絡(luò)會(huì)發(fā)生全局崩潰。因此，T的臨界閾值TC是衡量超網(wǎng)絡(luò)魯棒性的重要指標(biāo)。顯然，TC越小，表明超網(wǎng)絡(luò)越魯棒。

如果節(jié)點(diǎn)vj在獲得額外載荷之后發(fā)生失效，則其應(yīng)滿足不等式：

Lvj（t-1）+ΔLvi→vj（t）>Cvj（6）

如果式（6）成立，則節(jié)點(diǎn)vj將過載失效，當(dāng)vj的負(fù)載被重新分配后，可能會(huì)導(dǎo)致其它的節(jié)點(diǎn)發(fā)生失效。結(jié)合式（3），式（6）可以表示為

Lvj（t-1）+Lvi（t）×1dθH（vi）×1kθei（vi）>Cvj（7）

為了衡量超網(wǎng)絡(luò)的魯棒性，初始攻擊節(jié)點(diǎn)vi并使其失效，然后對(duì)它的負(fù)載進(jìn)行重新分配。對(duì)于負(fù)載重新分配后的其它節(jié)點(diǎn)，若滿足式（7），則該節(jié)點(diǎn)失效。當(dāng)超邊內(nèi)節(jié)點(diǎn)全失效后，則此條超邊發(fā)生失效。當(dāng)超網(wǎng)絡(luò)中的失效節(jié)點(diǎn)數(shù)達(dá)到一個(gè)穩(wěn)定狀態(tài)或者所有節(jié)點(diǎn)全部失效（全局崩潰）后，統(tǒng)計(jì)超網(wǎng)絡(luò)中的失效超邊數(shù)FM（0≤FM≤M），計(jì)算超邊失效比例fM如式（8）所示：

fM=FMM，0≤fM≤1（8）

其中，M為超網(wǎng)絡(luò)的超邊總數(shù)，由式（8）可知，fM越大說明超網(wǎng)絡(luò)中失效超邊數(shù)越多，即超網(wǎng)絡(luò)的魯棒性越差。

3 仿真分析

在分析k均勻無標(biāo)度超網(wǎng)絡(luò)超邊內(nèi)部結(jié)構(gòu)對(duì)超網(wǎng)絡(luò)魯棒性的影響時(shí)，本文以k均勻無標(biāo)度超網(wǎng)絡(luò)為底層網(wǎng)絡(luò)構(gòu)建了三類超網(wǎng)絡(luò)：BA_P超網(wǎng)絡(luò)、BA_R超網(wǎng)絡(luò)和BA_C超網(wǎng)絡(luò)，對(duì)三類無標(biāo)度超網(wǎng)絡(luò)在蓄意攻擊和隨機(jī)攻擊兩種策略下的級(jí)聯(lián)故障進(jìn)程進(jìn)行仿真模擬，并記錄了相關(guān)實(shí)驗(yàn)數(shù)據(jù)。仿真實(shí)驗(yàn)中的隨機(jī)攻擊方式為隨機(jī)選取超網(wǎng)絡(luò)中的一個(gè)節(jié)點(diǎn)進(jìn)行攻擊；而蓄意攻擊則是選取超網(wǎng)絡(luò)中最重要的節(jié)點(diǎn)進(jìn)行攻擊，因?yàn)楸疚难芯康氖浅W(wǎng)絡(luò)，所以在選擇節(jié)點(diǎn)重點(diǎn)性指標(biāo)時(shí)，選取了能體現(xiàn)超邊內(nèi)與超邊間結(jié)構(gòu)共同作用的指標(biāo)，即選擇超網(wǎng)絡(luò)中滿足節(jié)點(diǎn)度值與節(jié)點(diǎn)超度之和最大的節(jié)點(diǎn)進(jìn)行攻擊。

3.1 參數(shù)敏感性分析

每類超網(wǎng)絡(luò)的規(guī)模與參數(shù)k相關(guān)，因此為了進(jìn)行不同規(guī)模下的仿真分析，本文分別采用k為20， 40和60三種規(guī)模的網(wǎng)絡(luò)。本文研究的重點(diǎn)在于發(fā)現(xiàn)超邊內(nèi)部結(jié)構(gòu)對(duì)超網(wǎng)絡(luò)魯棒性的影響，為了消除其他不確定因素對(duì)魯棒性的影響，這里本文做了參數(shù)敏感性分析，如圖2所示。從圖2a中可以看出，三類超網(wǎng)絡(luò)對(duì)參數(shù)α的敏感性較弱，不論α的值如何增大，三類超網(wǎng)絡(luò)的魯棒性臨界閾值TC始終保持在一個(gè)相對(duì)穩(wěn)定的狀態(tài)。α值為10時(shí)，三類超網(wǎng)絡(luò)魯棒性最優(yōu)；由圖2b中可知，三類超網(wǎng)絡(luò)對(duì)參數(shù)β的敏感性大致相同，且魯棒性臨界閾值TC與參數(shù)β成正比例關(guān)系，也即魯棒性會(huì)隨著參數(shù)β的增大而變差。β值為1時(shí)，三類超網(wǎng)絡(luò)魯棒性最優(yōu)；由圖2c可以發(fā)現(xiàn)，三類超網(wǎng)絡(luò)對(duì)參數(shù)θ的敏感性較強(qiáng)，其中參數(shù)θ對(duì)BA_P超網(wǎng)絡(luò)的影響尤為明顯。三類超網(wǎng)絡(luò)的魯棒性臨界閾值TC與參數(shù)θ成反比例關(guān)系，即魯棒性隨著參數(shù)θ的增大而變強(qiáng)。θ值為1時(shí)，三類超網(wǎng)絡(luò)魯棒性最優(yōu)。

因此在仿真分析當(dāng)中使用控制變量法，選取使得三類超網(wǎng)絡(luò)魯棒性達(dá)到最優(yōu)時(shí)的參數(shù)值，即負(fù)載參數(shù)α取值為10，參數(shù)β取值為1，擾動(dòng)參數(shù)θ取值為1，為了保證結(jié)果的有效性與真實(shí)性，實(shí)驗(yàn)結(jié)果取50次結(jié)果的平均值。

3.2 三類k均勻BA無標(biāo)度超網(wǎng)絡(luò)魯棒性分析

圖3顯示了在不同超網(wǎng)絡(luò)規(guī)模下，三類k均勻無標(biāo)度超網(wǎng)絡(luò)在受到攻擊時(shí)，其超邊失效比例fM隨容量參數(shù)T的變化情況。圖中結(jié)果表明：在蓄意攻擊和隨機(jī)攻擊兩種策略下，超邊失效比例隨容量參數(shù)值的增大呈遞減趨勢(shì)，并在一定的容量參數(shù)下達(dá)到全局崩潰的臨界閾值。當(dāng)T≤TC時(shí)，三類k均勻無標(biāo)度超網(wǎng)絡(luò)均處于全局崩潰的狀態(tài)，當(dāng)T>TC時(shí)，故障規(guī)模開始減小并最終達(dá)到全局不失效的狀態(tài)。

當(dāng)k=20時(shí)，k均勻BA_P超網(wǎng)絡(luò)，在蓄意攻擊方式下，當(dāng)容量參數(shù)T=0.185，超網(wǎng)絡(luò)的崩潰比例fM=1，即此時(shí)超網(wǎng)絡(luò)發(fā)生全局崩潰；當(dāng)容量參數(shù)T>0.185時(shí)，超網(wǎng)絡(luò)的崩潰比例fM<1，因此超網(wǎng)絡(luò)在受到蓄意攻擊時(shí)其臨界閾值TC=0.185；同理分析可得，k均勻BA_P超網(wǎng)絡(luò)在受到隨機(jī)攻擊時(shí)其臨界閾值TC=0.058。由于TC越小，魯棒性越強(qiáng)，結(jié)合表2的結(jié)果可以明顯看出，k均勻BA_P超網(wǎng)絡(luò)在隨機(jī)攻擊下魯棒，在蓄意攻擊下脆弱。同理分析可得，k均勻BA_R超網(wǎng)絡(luò)在蓄意攻擊下魯棒，在隨機(jī)攻擊下脆弱；k均勻BA_C超網(wǎng)絡(luò)在蓄意攻擊下魯棒，在隨機(jī)攻擊下脆弱。

此外從圖3中可以明顯看出，隨著三類k均勻無標(biāo)度超網(wǎng)絡(luò)規(guī)模的增大，其臨界閾值TC均呈逐漸變小的趨勢(shì)。這是因?yàn)殡S著節(jié)點(diǎn)規(guī)模的增大，當(dāng)超網(wǎng)絡(luò)中的節(jié)點(diǎn)發(fā)生故障后，失效節(jié)點(diǎn)負(fù)載重分配時(shí)就會(huì)有更多的鄰居節(jié)點(diǎn)分擔(dān)其負(fù)載，使超網(wǎng)絡(luò)更不容易崩潰。因此，k均勻無標(biāo)度超網(wǎng)絡(luò)的規(guī)模越大，超網(wǎng)絡(luò)在受到攻擊時(shí)越魯棒。

以往的研究結(jié)果表明[1921]，當(dāng)忽略超邊內(nèi)節(jié)點(diǎn)的連接方式時(shí)，k均勻無標(biāo)度超網(wǎng)絡(luò)在遭受外部攻擊時(shí)與普通的無標(biāo)度網(wǎng)絡(luò)類似，即表現(xiàn)出對(duì)隨機(jī)攻擊魯棒，對(duì)蓄意攻擊脆弱的特性。但本文分析發(fā)現(xiàn)，當(dāng)k均勻無標(biāo)度超網(wǎng)絡(luò)超邊內(nèi)部節(jié)點(diǎn)的連接方式不同時(shí)，k均勻無標(biāo)度超網(wǎng)絡(luò)表現(xiàn)出了不同的魯棒性，即當(dāng)超網(wǎng)絡(luò)超邊內(nèi)節(jié)點(diǎn)采用優(yōu)先連接方式時(shí)，超網(wǎng)絡(luò)表現(xiàn)出了與已有研究結(jié)果一致的魯棒性，而當(dāng)超邊內(nèi)的節(jié)點(diǎn)采用其他兩種連接方式時(shí)，超網(wǎng)絡(luò)表現(xiàn)出了相反的魯棒性，但在兩種攻擊方式下的魯棒性差異不大，且超邊內(nèi)節(jié)點(diǎn)的規(guī)模對(duì)其有重要的影響，即超邊內(nèi)節(jié)點(diǎn)越多，則超邊內(nèi)節(jié)點(diǎn)的連接方式產(chǎn)生的影響也越大。同時(shí)，本文對(duì)超邊內(nèi)節(jié)點(diǎn)連接方式對(duì)魯棒性產(chǎn)生影響的原因進(jìn)行了具體的分析。這里k均勻無標(biāo)度超網(wǎng)絡(luò)的各超邊之間的連接方式為優(yōu)先連接。

1）當(dāng)超邊內(nèi)部節(jié)點(diǎn)之間的連接方式為優(yōu)先連接時(shí)，每條超邊內(nèi)部的節(jié)點(diǎn)度分布與節(jié)點(diǎn)超度分布均為冪律分布，如圖4和圖5所示。這表明這種連接方式加強(qiáng)了超網(wǎng)絡(luò)的無標(biāo)度特性，因此BA_P超網(wǎng)絡(luò)在受到攻擊時(shí)與無標(biāo)度超網(wǎng)絡(luò)表現(xiàn)出一致的特性，即對(duì)隨機(jī)攻擊魯棒，對(duì)蓄意攻擊脆弱；

2）當(dāng)超邊內(nèi)部節(jié)點(diǎn)之間的連接方式為隨機(jī)連接時(shí)，每條超邊內(nèi)部的節(jié)點(diǎn)度分布為泊松分布，又因節(jié)點(diǎn)的超度分布仍為冪律分布，所以BA_R超網(wǎng)絡(luò)整體的節(jié)點(diǎn)度分布在超度分布的影響下就會(huì)呈現(xiàn)偏態(tài)的泊松分布，如圖6所示。因此超邊內(nèi)節(jié)點(diǎn)的隨機(jī)連接方式改變了超網(wǎng)絡(luò)整體的無標(biāo)度特性，所以BA_R超網(wǎng)絡(luò)在受到攻擊時(shí)呈現(xiàn)出與ER隨機(jī)網(wǎng)絡(luò)相同的特性，即對(duì)隨機(jī)攻擊脆弱，對(duì)蓄意攻擊魯棒；

3）當(dāng)超邊內(nèi)部節(jié)點(diǎn)之間采用全連接的方式時(shí)，每條超邊內(nèi)部為一個(gè)規(guī)則網(wǎng)絡(luò)，因?yàn)楣?jié)點(diǎn)超度分布為冪律分布，如圖4所示，則當(dāng)一個(gè)超度較大的節(jié)點(diǎn)遭受到攻擊時(shí)，介于本文所規(guī)定的負(fù)載傳播方式，如第2部分的式（3）所示，失效節(jié)點(diǎn)在進(jìn)行負(fù)載傳播的過程中，當(dāng)超度大的節(jié)點(diǎn)vi其負(fù)載傳播時(shí)，分配到其相鄰未失效節(jié)點(diǎn)vj所在超邊上的負(fù)載就會(huì)相對(duì)更少，為Lvit/d θHvi。又因超邊內(nèi)部節(jié)點(diǎn)采用全連接的方式，那么失效節(jié)點(diǎn)vi的每個(gè)未失效的鄰居節(jié)點(diǎn)vj就會(huì)分到相對(duì)更少的負(fù)載，即Lvit/d θHvikθeivi。所以此時(shí)超網(wǎng)絡(luò)會(huì)表現(xiàn)出對(duì)蓄意攻擊魯棒，對(duì)隨機(jī)攻擊脆弱的特性。此外，Ma等[20]研究的無標(biāo)度超網(wǎng)絡(luò)雖然也采用了超邊內(nèi)全連接的方式，但其研究的k均勻無標(biāo)度超網(wǎng)絡(luò)中k=3，也即超邊內(nèi)節(jié)點(diǎn)的規(guī)模僅為3。因此，在Ma等的研究[20]中，超邊內(nèi)節(jié)點(diǎn)的連接方式對(duì)超網(wǎng)絡(luò)魯棒性產(chǎn)生的影響非常小，不足以改變超網(wǎng)絡(luò)整體的魯棒性。本文研究發(fā)現(xiàn)，當(dāng)k均勻無標(biāo)度超網(wǎng)絡(luò)超邊內(nèi)節(jié)點(diǎn)間全連接，且k≥8時(shí)，會(huì)對(duì)k均勻無標(biāo)度超網(wǎng)絡(luò)的魯棒性產(chǎn)生重要的影響，甚至改變其原有的魯棒性，如圖7所示。

另外，結(jié)合表2中的數(shù)據(jù)分析得出，k均勻無標(biāo)度超網(wǎng)絡(luò)在面對(duì)蓄意攻擊時(shí)，BA_C超網(wǎng)絡(luò)的魯棒性最好，其次是BA_R超網(wǎng)絡(luò)，最差的是BA_P超網(wǎng)絡(luò)；分析其原因?yàn)椋簁均勻無標(biāo)度超網(wǎng)絡(luò)在面對(duì)蓄意攻擊時(shí)的魯棒性與超邊內(nèi)部的普通邊數(shù)量mk有關(guān)。超邊內(nèi)部不同結(jié)構(gòu)的普通邊數(shù)量如圖8所示。因?yàn)樵谛钜夤魰r(shí)，攻擊的是超邊內(nèi)部超度占比與節(jié)點(diǎn)度占比之和最大的節(jié)點(diǎn)，普通邊越多，就會(huì)有更多的鄰居節(jié)點(diǎn)在負(fù)載傳播的過程中分擔(dān)其負(fù)載。所以在蓄意攻擊策略下，超邊內(nèi)部普通邊的數(shù)量mk越大，k均勻無標(biāo)度超網(wǎng)絡(luò)的魯棒性也越好。而k均勻無標(biāo)度超網(wǎng)絡(luò)在面對(duì)隨機(jī)攻擊時(shí)，因隨機(jī)攻擊具有一定的隨機(jī)性，又因?yàn)闊o標(biāo)度超網(wǎng)絡(luò)本身的無標(biāo)度特性，從而減弱了超邊內(nèi)部普通邊的數(shù)量對(duì)魯棒性的影響。

4 結(jié)論

為了突破現(xiàn)有超網(wǎng)絡(luò)結(jié)構(gòu)魯棒性研究的局限性，探究超邊內(nèi)部結(jié)構(gòu)對(duì)超網(wǎng)絡(luò)魯棒性的影響，本文提出三類超邊內(nèi)部結(jié)構(gòu)不同的k均勻無標(biāo)度超網(wǎng)絡(luò)模型，并提出一種新的容量負(fù)載模型分析了超邊內(nèi)部結(jié)構(gòu)對(duì)超網(wǎng)絡(luò)整體魯棒性的影響。得到結(jié)論：當(dāng)超邊內(nèi)節(jié)點(diǎn)的連接方式與超網(wǎng)絡(luò)各個(gè)超邊間的連接方式一致時(shí)，k均勻無標(biāo)度超網(wǎng)絡(luò)對(duì)隨機(jī)攻擊魯棒，對(duì)蓄意攻擊脆弱；當(dāng)超邊內(nèi)節(jié)點(diǎn)的連接方式與超網(wǎng)絡(luò)各個(gè)超邊間的連接方式不一致時(shí)，k均勻無標(biāo)度超網(wǎng)絡(luò)對(duì)隨機(jī)攻擊脆弱，對(duì)蓄意攻擊魯棒；當(dāng)超邊內(nèi)節(jié)點(diǎn)的規(guī)模k較大時(shí)，其連接方式對(duì)k均勻無標(biāo)度超網(wǎng)絡(luò)魯棒性的影響也較大；當(dāng)超邊內(nèi)普通邊數(shù)量mk越大，k均勻無標(biāo)度超網(wǎng)絡(luò)魯棒性越好。

本文的研究結(jié)果對(duì)進(jìn)一步理解微觀結(jié)構(gòu)對(duì)系統(tǒng)整體性能的影響有重要的意義，同時(shí)也為進(jìn)一步優(yōu)化超網(wǎng)絡(luò)的結(jié)構(gòu)魯棒性提供了重要的參考。

參考文獻(xiàn)：

[1]REN H P， GAO Y， HUO L， et al. Frequency stability in modern power network from complex network viewpoint[J]. Physica A， 2020， 545（3）： 123558.

[2]WANG T， CHENG H， WANG X. A link addition method based on uniformity of node degree in interdependent power grids and communication networks[J]. Physica A， 2020， 560（5）： 125112.

[3]CAI W X， LIANG F F， WANG Y C， et al. An innovative approach for constructing a shipping index based on dynamic weighted complex networks[J]. Physica A， 202 578（5）： 126101.

[4]毛昌梅， 韓景倜， 劉舉勝. 基于復(fù)雜網(wǎng)絡(luò)的銀行波動(dòng)溢出效應(yīng)研究[J]. 復(fù)雜系統(tǒng)與復(fù)雜性科學(xué)， 2020， 17（2）： 1121.

MAO C M， HAN J T， LIU J S. Volatility spillover effect of chinese listed commercial banks based on complex network[J]. Complex Systems and Complexity Science， 2020， 17（2）： 1121.

[5]王哲， 李建華， 康東， 等. 復(fù)雜網(wǎng)絡(luò)魯棒性增強(qiáng)策略研究綜述[J]. 復(fù)雜系統(tǒng)與復(fù)雜性科學(xué)， 2020， 17（3）： 126+46.

WANG Z， LI J H， KANG D， et al. Review on strategies enhancing the robustness of complex network[J]. Complex Systems and Complexity Science， 2020， 17（3）： 126，46.

[6]WANG S L， L W Z， ZHANG J H， et al. Method of power network critical nodes identification and robustness enhancement based on a cooperative framework[J]. Reliability Engineering & System Safety， 202 207：107313.

[7]GAO J， BULDYREV S V， HAVLIN S， et al. Robustness of a network formed by n interdependent networks with a one-to-one correspondence of dependent nodes[J]. Physical Review E， 2012， 85（6）： 066134.

[8]陳關(guān)榮. 探索復(fù)雜網(wǎng)絡(luò)的高階拓?fù)浼捌鋺?yīng)用[R]. 北京： 中國指揮與控制學(xué)會(huì)， 2021.

CHEN G R. Exploring higher-order topologies of complex networks and applications[R]. Beijing： Chinese Institute of Command and Control.

[9]FEDERICO B， GIULIA C， IACOPO I， et al. Networks beyond pairwise interactions： structure and dynamics[J]. Physics Reports， 2020， 874： 192.

[10] SINAN G A， CLIFF J， CARLOS O M， et al. Hypernetwork science via high-order hypergraph walks[J]. EPJ Data Science， 2020， 9（1）： 519535.

[11] WANG J W， RONG LL， DENG Q H， et al. Evolving hypernetwork model[J]. The European Physical Journal B， 2010， 77（4）： 493498.

[12] ESTRADA E， RODRGUEZ-VELZQUEZ J A. Subgraph centrality and clustering in complex hyper-networks[J]. Physica A， 2006， 364（3）： 581594.

[13] 索琪， 郭進(jìn)利. 基于超圖的超網(wǎng)絡(luò)：結(jié)構(gòu)及演化機(jī)制[J]. 系統(tǒng)工程理論與實(shí)踐， 2017， 37（3）： 720734.

SUO Q， GUO J L. The structure and dynamics of hypernetworks[J]. Systems Engineering-Theory & Practice， 2017， 37（3）： 720734.

[14] 李甍娜， 郭進(jìn)利， 卞聞， 等. 網(wǎng)絡(luò)視角下的唐詩[J]. 復(fù)雜系統(tǒng)與復(fù)雜性科學(xué)， 2017， 14（4）： 6671.

LI M N， GUO J L， BIAN W， et al. Tang poetry from the perspective of network[J]. Complex Systems and Complexity Science， 2017， 14（4）： 6671.

[15] 胡楓， 趙海興， 何佳倍， 等. 基于超圖結(jié)構(gòu)的科研合作網(wǎng)絡(luò)演化模型[J]. 物理學(xué)報(bào)， 2013， 62（19）： 547554.

HU F， ZHAO H X， HE J B， et al. An evolving model for hypergraph-structure-based scientific collaboration networks[J]. Acta Physica Sinica， 2013， 62（19）： 547554.

[16] 胡楓， 劉猛， 趙靜， 等. 蛋白復(fù)合物超網(wǎng)絡(luò)特性分析及應(yīng)用[J]. 復(fù)雜系統(tǒng)與復(fù)雜性科學(xué)， 2018， 15（4）： 3138.

HU F， LIU M， ZHAO J， et al. Analysisand application of the topological properties of protein complex hypernetworks[J]. Complex Systems and Complexity Science， 2018， 15（4）： 3138.

[17] Wu Z Y， Duan J Q， Fu X C. Synchronization of an evolving complex hyper-network[J]. Applied Mathematical Modelling， 2014， 38（11/12）： 29612968.

[18] 鞏云超， 李發(fā)旭， 周麗娜， 等. 在線社交超網(wǎng)絡(luò)的信息全局傳播模型[J]. 電子科技大學(xué)學(xué)報(bào)， 202 50（3）： 437445.

GONG Y C， LI F X， ZHOU L N， et al. Global Dissemination of Information Based on Online Social Hypernetwork[J]. Journal of University of Electronic Science and Technology of China， 202 50（3）： 437445.

[19] 馬秀娟， 趙海興， 胡楓. 基于超圖的超網(wǎng)絡(luò)相繼故障分析[J]. 物理學(xué)報(bào)， 2016， 65（8）： 374383.

MA X J， ZHAO H X， HU F. Cascading failure analysis in hyper-network based on the hypergraph[J]. Acta Physica Sinica， 2016， 65（8）： 374383.

[20] MA X J， MA F X， YIN J， et al. Cascading failures of k uniform hyper-network based on the hyper adjacent matrix[J]. Physica A， 2018， 510： 281289.

[21] CHEN Y， MA X J， MA F X， et al. The capacity load model of k-uniform hyper-network based on equal load distribution[J]. Journal of Physics： Conference Series， 202 1828（1）： 012060.

[22] 羅海秀， 趙海興， 肖玉芝， 等. 基于超圖的公交超網(wǎng)絡(luò)拓?fù)涮匦约棒敯粜苑治鯷J]. 西南大學(xué)學(xué)報(bào)（自然科學(xué)版）， 202 43（10）： 181191.

LUO H X， ZHAO H X， XIAO Y Z. A hypergraph-based analysis of the topology and robustness of bus hypernetworks[J]. Journal of Southwest University （Natural Science Edition）， 202 43（10）： 181191.

[23] MOTTER A E， CHENG L Y. Cascade-based attacks on complex networks[J]. Physical Review E， 2002， 66（6）： 065102.

[24] BERGE C. Graphs and Hpergraphs[M]. New York： Elsevier， 1973.

[25] BERGE C， STERBOUL F. Equipartite colorings in graphs and hypergraphs[J]. Journal of Combinatorial Theory， Series B， 1977， 22（2）： 97113.

[26] BRETTO A. Hypergraph Theory[M]. Heidelberg： Springer， 2013.

[27] 胡楓， 趙海興， 馬秀娟. 一種超網(wǎng)絡(luò)演化模型構(gòu)建及特性分析[J]. 中國科學(xué)： 物理學(xué) 力學(xué) 天文學(xué)， 2013， 43（1）： 1622.

HU F， ZHAO H X， MA X J. An evolving hypernetwork model and its properties[J]. Scientia Sinica（Physica， Mechanica& Astronomica）， 2013， 43（1）： 1622.

（責(zé)任編輯 李 進(jìn)）