胡博仁　裴忠民　羅章凱　丁杰

摘要： 鑒于開展衛(wèi)星網(wǎng)絡(luò)局部結(jié)構(gòu)研究是理解網(wǎng)絡(luò)性質(zhì)的重要手段，考慮星間鏈路天線可見性約束等條件，提出了一種基于含時(shí)網(wǎng)絡(luò)數(shù)據(jù)的衛(wèi)星網(wǎng)絡(luò)模體識(shí)別方法，建立了從TLE文件輸入到子結(jié)構(gòu)識(shí)別輸出的模體識(shí)別流程。以GPS衛(wèi)星網(wǎng)絡(luò)的三節(jié)點(diǎn)三邊模體識(shí)別為例，結(jié)果發(fā)現(xiàn)在短時(shí)段內(nèi)衛(wèi)星天線最大掃描范圍與具有特殊結(jié)構(gòu)意義的三角形M4子圖濃度呈正相關(guān)。

關(guān)鍵詞： 衛(wèi)星網(wǎng)絡(luò)；模體識(shí)別；子圖濃度

中圖分類號(hào)： TP393；TP399 文獻(xiàn)標(biāo)識(shí)碼： A

On Motif Counts Method of Satellite Network Based on Temporal Network Data

HU Borena， PEI Zhongmina，LUO Zhangkaia，DING Jieb

（a. Science and Technology on Complex Electronic System Simulation Laboratory;b. Department of Electronics and Optical Engineering， Space Engineering University， Beijing? 101416，China）

Abstract：Conducting research on the local structure of satellite networks is an important means to understand the nature of networks. Considering the visibility constraints of inter-satellite link antennas， a satellite network motif counts method based on temporal network data is proposed， and a motif counts process is established from TLE file input to substructure identification output; taking the three-node three-edge motif counts of GPS satellite network as an example， we found that in a short period of time maximum scanning range of the satellite antenna was positively correlated with the concentration of the triangular M4 subgraph with special structural significance.

Keywords： satellite network; motif counts; subgraph concentration

0 引言

近年來，隨著信息技術(shù)和航天技術(shù)的不斷發(fā)展，在軌衛(wèi)星數(shù)目的增加造成了衛(wèi)星網(wǎng)絡(luò)節(jié)點(diǎn)密集化及網(wǎng)絡(luò)結(jié)構(gòu)復(fù)雜度的提高。從空間資源高效利用和確保網(wǎng)絡(luò)安全運(yùn)行的角度出發(fā)，開展衛(wèi)星網(wǎng)絡(luò)結(jié)構(gòu)特性分析具有重要意義。復(fù)雜網(wǎng)絡(luò)作為網(wǎng)絡(luò)科學(xué)的重要分支，通過研究結(jié)構(gòu)特性和組織方式來探索和理解復(fù)雜系統(tǒng)的結(jié)構(gòu)、行為和功能，為衛(wèi)星網(wǎng)絡(luò)研究提供了有力的方法支持。廣大學(xué)者通過復(fù)雜網(wǎng)絡(luò)的相關(guān)理論對(duì)衛(wèi)星網(wǎng)絡(luò)進(jìn)行了深入的研究。武健等[1]基于復(fù)雜網(wǎng)絡(luò)理論和灰色關(guān)聯(lián)度分析方法，提出了網(wǎng)絡(luò)微觀層面的衛(wèi)星節(jié)點(diǎn)重要度評(píng)估指標(biāo)。朱林等[2]綜合節(jié)點(diǎn)介數(shù)、節(jié)點(diǎn)緊密度和節(jié)點(diǎn)距離對(duì)衛(wèi)星網(wǎng)絡(luò)節(jié)點(diǎn)重要性的貢獻(xiàn)，提出了穩(wěn)態(tài)衛(wèi)星網(wǎng)絡(luò)節(jié)點(diǎn)重要性評(píng)估方法。王瑩[3]提出了衡量衛(wèi)星通信網(wǎng)絡(luò)整體結(jié)構(gòu)性能的衛(wèi)星移動(dòng)通信網(wǎng)約束連通度指標(biāo)，重點(diǎn)在于度量節(jié)點(diǎn)間滿足約束條件的可用路徑。林琪等[4]基于衛(wèi)星網(wǎng)絡(luò)拓?fù)涮卣鳎ǘ确植?、平均最短路徑和網(wǎng)絡(luò)直徑等，提出了衛(wèi)星網(wǎng)絡(luò)綜合效能評(píng)估方法。

目前，運(yùn)用復(fù)雜網(wǎng)絡(luò)理論對(duì)衛(wèi)星網(wǎng)絡(luò)進(jìn)行結(jié)構(gòu)分析和研究時(shí)，多從網(wǎng)絡(luò)的整體結(jié)構(gòu)特性出發(fā)，如網(wǎng)絡(luò)連通度、度分布等，或注重節(jié)點(diǎn)或邊的屬性，如節(jié)點(diǎn)重要性評(píng)估，鮮有衛(wèi)星網(wǎng)絡(luò)局部結(jié)構(gòu)相關(guān)的研究。然而局部結(jié)構(gòu)作為網(wǎng)絡(luò)結(jié)構(gòu)分析的中尺度視角，與衛(wèi)星網(wǎng)絡(luò)結(jié)構(gòu)的微觀、宏觀分析角度同樣重要，有時(shí)局部結(jié)構(gòu)特征能更好地揭示網(wǎng)絡(luò)結(jié)構(gòu)與功能的內(nèi)在關(guān)系。模體概念最早于2002年由Milo等[5]提出，定義為網(wǎng)絡(luò)中重復(fù)出現(xiàn)的局部子圖結(jié)構(gòu)。如靜態(tài)三節(jié)點(diǎn)有向圖共有圖1中13種，其中某一種子圖在網(wǎng)絡(luò)中多次出現(xiàn)，且滿足{P，U，D，N}[5]條件，即可稱之為模體。

Martí Rosas-Casals等[6]對(duì)歐洲電網(wǎng)的結(jié)構(gòu)穩(wěn)定性進(jìn)行研究，發(fā)現(xiàn)相對(duì)于分散的去中心化連接模式，四節(jié)點(diǎn)三邊規(guī)模的星型模體數(shù)量增加會(huì)加劇網(wǎng)絡(luò)的脆弱性。Paul Schultz等[7]提出是否存在某種網(wǎng)絡(luò)模體能提高電網(wǎng)結(jié)構(gòu)穩(wěn)定性的問題，發(fā)現(xiàn)彎路模體（Detours motifs）在提高網(wǎng)絡(luò)穩(wěn)定性上具有重要作用，三角形模體即為最簡(jiǎn)單的彎路模體（見圖2）。孫曉偉[8]以引文網(wǎng)絡(luò)、合作網(wǎng)絡(luò)和作者引用網(wǎng)絡(luò)的三階網(wǎng)絡(luò)模體為研究對(duì)象，挖掘論文和作者間引用、合作關(guān)系的演化規(guī)律。

相比電力網(wǎng)絡(luò)和論文引文網(wǎng)絡(luò)等靜止網(wǎng)絡(luò)，衛(wèi)星網(wǎng)絡(luò)具有動(dòng)態(tài)高時(shí)變、網(wǎng)絡(luò)數(shù)據(jù)獲取難等特點(diǎn)，進(jìn)行衛(wèi)星網(wǎng)絡(luò)局部子圖結(jié)構(gòu)研究具有一定的難度。本文旨在為衛(wèi)星網(wǎng)絡(luò)結(jié)構(gòu)分析中引入中尺度的模體角度，提出一種基于含時(shí)網(wǎng)絡(luò)數(shù)據(jù)的衛(wèi)星網(wǎng)絡(luò)模體識(shí)別方法，解決衛(wèi)星網(wǎng)絡(luò)模型建立、網(wǎng)絡(luò)數(shù)據(jù)構(gòu)建和動(dòng)態(tài)模體識(shí)別過程中條件不清晰、數(shù)據(jù)量大等問題。該研究對(duì)于揭示衛(wèi)星網(wǎng)絡(luò)局部結(jié)構(gòu)與整體性質(zhì)之間的關(guān)系有著一定的意義，為衛(wèi)星網(wǎng)絡(luò)局部結(jié)構(gòu)研究提供了方法借鑒和技術(shù)參考。

1 基于含時(shí)網(wǎng)絡(luò)數(shù)據(jù)的衛(wèi)星網(wǎng)絡(luò)模體識(shí)別方法

識(shí)別衛(wèi)星網(wǎng)絡(luò)中頻繁出現(xiàn)的局部高階子結(jié)構(gòu)，需考慮網(wǎng)絡(luò)的動(dòng)態(tài)時(shí)變性并自行構(gòu)建衛(wèi)星網(wǎng)絡(luò)數(shù)據(jù)。針對(duì)上述問題，本文提出了一種基于含時(shí)網(wǎng)絡(luò)數(shù)據(jù)的衛(wèi)星網(wǎng)絡(luò)模體識(shí)別方法，流程如圖3所示。借助STK和Matlab等軟件，建立仿真衛(wèi)星網(wǎng)絡(luò)模型，設(shè)置網(wǎng)絡(luò)節(jié)點(diǎn)之間的連接條件；生成包含時(shí)間屬性的衛(wèi)星網(wǎng)絡(luò)數(shù)據(jù)；使用動(dòng)態(tài)模體識(shí)別算法進(jìn)行子圖計(jì)數(shù)并計(jì)算子圖濃度等指標(biāo)，以此來進(jìn)行衛(wèi)星網(wǎng)絡(luò)結(jié)構(gòu)分析。

1.1 建立衛(wèi)星網(wǎng)絡(luò)模型

本文研究的衛(wèi)星網(wǎng)絡(luò)，以衛(wèi)星星座為主體，泛指由存在于近地空間的衛(wèi)星節(jié)點(diǎn)和地面段的地面站節(jié)點(diǎn)組成，節(jié)點(diǎn)之間以星間鏈路和星地鏈路兩種形式建立連接，用于信息傳輸和星座業(yè)務(wù)實(shí)現(xiàn)的網(wǎng)絡(luò)系統(tǒng)，如圖4所示。

下載Celestrak網(wǎng)站的衛(wèi)星星歷（TLE）文件，導(dǎo)入STK，建立近地空間的衛(wèi)星星座模型，或者根據(jù)軌道六參數(shù)，確定衛(wèi)星的軌道和位置。衛(wèi)星網(wǎng)絡(luò)中的地面站節(jié)點(diǎn)設(shè)置與具體衛(wèi)星系統(tǒng)類型有關(guān)，如衛(wèi)星導(dǎo)航系統(tǒng)的地面段包括主控站、注入站和監(jiān)測(cè)站。查閱資料得到站點(diǎn)的經(jīng)度、緯度信息，建立地面站點(diǎn)，或查詢STK自帶地面站庫，直接添加。

衛(wèi)星網(wǎng)絡(luò)中的邊包括星地鏈路和星間鏈路，根據(jù)衛(wèi)星網(wǎng)絡(luò)實(shí)際運(yùn)行情況，結(jié)合簡(jiǎn)化假設(shè)，設(shè)置模型中邊的連接準(zhǔn)則。星間鏈路的建立條件包括[9]：幾何可視條件、天線可視條件和傳輸距離條件。星間鏈路天線可見性的約束條件[9]：

1.2 構(gòu)建衛(wèi)星網(wǎng)絡(luò)數(shù)據(jù)

含時(shí)網(wǎng)絡(luò)[1011]在復(fù)雜網(wǎng)絡(luò)模型的基礎(chǔ)上，加入了時(shí)間維度，用（V，E，t）表示，其中t代表網(wǎng)絡(luò)連邊的發(fā)生時(shí)刻，主要用來刻畫離散時(shí)間內(nèi)網(wǎng)絡(luò)連邊斷續(xù)存在的情形。比如在衛(wèi)星網(wǎng)絡(luò)進(jìn)行數(shù)據(jù)傳輸時(shí)，A和B之間的連邊僅在進(jìn)行數(shù)據(jù)傳輸時(shí)存在，數(shù)據(jù)傳輸結(jié)束后，連邊也隨之消失。只在A、B之間建立一條連邊已無法描述節(jié)點(diǎn)相互作用時(shí)刻變化的特點(diǎn)。

網(wǎng)絡(luò)的動(dòng)態(tài)變化具體表現(xiàn)為邊數(shù)量的增減或節(jié)點(diǎn)數(shù)量的增減，而衛(wèi)星網(wǎng)絡(luò)的節(jié)點(diǎn)數(shù)量通常不發(fā)生變化，邊的數(shù)量隨時(shí)間不斷變化。用三元組（u，v，t）表示衛(wèi)星網(wǎng)絡(luò)中的連邊，u表示邊起點(diǎn)，v表示邊終點(diǎn)，t表示u和v在t時(shí)刻處于連接狀態(tài)。

Matlab與STK互聯(lián)之后，可以調(diào)取衛(wèi)星節(jié)點(diǎn)和地面站節(jié)點(diǎn)的相關(guān)數(shù)據(jù)，包括某時(shí)刻下，衛(wèi)星之間的距離和衛(wèi)星的高程，衛(wèi)星與地面站是否滿足幾何可見性約束等。依據(jù)衛(wèi)星網(wǎng)絡(luò)模型中邊的連接準(zhǔn)則并進(jìn)行條件判斷之后，獲取某一時(shí)間段內(nèi)的網(wǎng)絡(luò)三元組邊數(shù)據(jù)，完成衛(wèi)星網(wǎng)絡(luò)數(shù)據(jù)的構(gòu)建。

1.3 識(shí)別衛(wèi)星網(wǎng)絡(luò)模體

網(wǎng)絡(luò)模體識(shí)別是網(wǎng)絡(luò)模體研究的重點(diǎn)和難點(diǎn)之一。經(jīng)典的模體識(shí)別算法[12]主要是面向靜態(tài)網(wǎng)絡(luò)，生成若干個(gè)與實(shí)證網(wǎng)絡(luò)節(jié)點(diǎn)數(shù)和節(jié)點(diǎn)的度序列相同的隨機(jī)網(wǎng)絡(luò)；在實(shí)證網(wǎng)絡(luò)和隨機(jī)網(wǎng)絡(luò)中搜索某一規(guī)模的子圖，將同構(gòu)的子圖歸為一類；比較每一類子圖在實(shí)證網(wǎng)絡(luò)和隨機(jī)網(wǎng)絡(luò)中的出現(xiàn)次數(shù)以確定其統(tǒng)計(jì)意義，從而確定是否為網(wǎng)絡(luò)模體。常用的靜態(tài)網(wǎng)絡(luò)模體識(shí)別工具包括MFinder[13]，F(xiàn)ANMOD[14]，MODA[15]和NemoMap[16]等。

衛(wèi)星網(wǎng)絡(luò)是具有時(shí)變性的動(dòng)態(tài)網(wǎng)絡(luò)。面向動(dòng)態(tài)網(wǎng)絡(luò)的模體識(shí)別算法包括基于流模型的StreaM[17]和Massive Streaming Data Analytics[18]，SNAP（Stanford Network Analysis Project）框架下的Motifs in temporal networks 動(dòng)態(tài)模體識(shí)別抽樣算法[19]和oDEN算法[20]等。

構(gòu)建衛(wèi)星網(wǎng)絡(luò)數(shù)據(jù)后，選用SNAP框架下面向大規(guī)模網(wǎng)絡(luò)數(shù)據(jù)的抽樣動(dòng)態(tài)三節(jié)點(diǎn)模體識(shí)別算法，該算法適用于二節(jié)點(diǎn)三邊和三節(jié)點(diǎn)三邊的36種有向子圖計(jì)數(shù)，如圖6所示。由于網(wǎng)絡(luò)帶有時(shí)間屬性，引入時(shí)間參數(shù)δ，并定義δ衛(wèi)星網(wǎng)絡(luò)模體：由具有時(shí)間屬性的邊組成的子圖，模體中的有向時(shí)空邊（帶時(shí)間戳）具有先后順序，受時(shí)間段δ約束，即連接關(guān)系均發(fā)生在時(shí)間段δ內(nèi)，用數(shù)學(xué)表達(dá)式可表示為

其中，t1，t2，…為衛(wèi)星網(wǎng)絡(luò)模體中的邊連接時(shí)刻。

本文建立的衛(wèi)星網(wǎng)絡(luò)模型，衛(wèi)星間連接和衛(wèi)星與地面的連接均為無向邊，三節(jié)點(diǎn)三邊無向圖共4種（見圖7），需將算法識(shí)別的有向子圖轉(zhuǎn)換為無向子圖。

對(duì)無向衛(wèi)星網(wǎng)絡(luò)數(shù)據(jù)進(jìn)行模體識(shí)別時(shí)，有向圖中有帶環(huán)的子圖如M1，2，會(huì)出現(xiàn)重復(fù)計(jì)數(shù)的現(xiàn)象，使得計(jì)數(shù)結(jié)果多于實(shí)際無向圖，故不考慮。 SNAP框架下的temporal motifs 算法為提高識(shí)別速度運(yùn)用了抽樣思想，使對(duì)應(yīng)同一種無向子圖的有向子圖計(jì)數(shù)結(jié)果存在細(xì)微差別，故取有向圖計(jì)數(shù)結(jié)果的均值為無向圖計(jì)數(shù)結(jié)果，具體計(jì)算見公式（4）：

4種三節(jié)點(diǎn)三邊無向圖中，M1、M2、M3均為星型子圖，指單節(jié)點(diǎn)為中心節(jié)點(diǎn)，其他節(jié)點(diǎn)直接與中心節(jié)點(diǎn)相連構(gòu)成的子圖，M4為三角形彎路子圖，有利于網(wǎng)絡(luò)結(jié)構(gòu)的穩(wěn)定，重點(diǎn)關(guān)注子圖M4的出現(xiàn)次數(shù)和子圖濃度。子圖濃度是指相同實(shí)驗(yàn)條件下，同等子圖規(guī)模的某種子圖所占的比例，體現(xiàn)了網(wǎng)絡(luò)中同等規(guī)模不同連接模式子圖的分布情況，計(jì)算如式（5）所示。

其中，Ck為某子圖規(guī)模下第k個(gè)子圖結(jié)構(gòu)的子圖濃度，Mk為第k個(gè)子圖的出現(xiàn)次數(shù)，該子圖規(guī)模下，共有N個(gè)異構(gòu)子圖結(jié)構(gòu)。

2 實(shí)例分析

以GPS為例，建立由空間在軌衛(wèi)星和地面監(jiān)測(cè)站組成的衛(wèi)星網(wǎng)絡(luò)模型。GPS由空間段、地面段和用戶終端組成，本文不考慮用戶終端?？臻g段由30顆中軌道衛(wèi)星組成。地面段包括主控站、監(jiān)測(cè)站及注入站，主控站[21]主要是收集和處理監(jiān)測(cè)站的觀測(cè)數(shù)據(jù)；監(jiān)測(cè)站利用復(fù)雜的GPS接收機(jī)跟蹤從監(jiān)測(cè)站上空經(jīng)過的GPS衛(wèi)星，收集導(dǎo)航信號(hào)、范圍測(cè)量數(shù)據(jù)和大氣數(shù)據(jù)等；注入站在衛(wèi)星離開其作用范圍之前進(jìn)行指令等信息注入。綜上，為簡(jiǎn)化模型，故本文只考慮GPS衛(wèi)星與監(jiān)測(cè)站之間的信息傳輸，建立由空間段衛(wèi)星節(jié)點(diǎn)和監(jiān)測(cè)站節(jié)點(diǎn)組成的衛(wèi)星網(wǎng)絡(luò)，且不考慮地面節(jié)點(diǎn)間的連接。在STK軟件中建立的16個(gè)GPS監(jiān)測(cè)站點(diǎn)如圖8所示。

使用Celestrak網(wǎng)站的兩行根數(shù)（TLE）文件，構(gòu)建GPS空間段。TLE文件的時(shí)間為2021年12月21日，導(dǎo)入STK軟件，衛(wèi)星節(jié)點(diǎn)數(shù)量為30，地面站節(jié)點(diǎn)和衛(wèi)星節(jié)點(diǎn)數(shù)目總和為46，3D示意如圖9所示。

設(shè)定節(jié)點(diǎn)間的連邊條件。某時(shí)刻下，若衛(wèi)星間的距離lAB 滿足星間鏈路的幾何可見性約束及天線可見性約束，則認(rèn)為衛(wèi)星節(jié)點(diǎn)間建立了雙向連接，不考慮傳輸距離等因素的影響；某時(shí)刻下，若地面站節(jié)點(diǎn)與衛(wèi)星節(jié)點(diǎn)滿足可見性約束，則認(rèn)為建立了從衛(wèi)星到地面站的雙向連接，不考慮時(shí)延和誤碼等因素。在星間鏈路建立時(shí)，天線可見性約束（見式（1））中的衛(wèi)星天線最大掃描范圍αmax起關(guān)鍵作用，對(duì)網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)具有較大影響[22]。結(jié)合實(shí)際衛(wèi)星網(wǎng)絡(luò)運(yùn)行情況并考慮數(shù)據(jù)量規(guī)模等問題，設(shè)置實(shí)驗(yàn)時(shí)長(zhǎng)為86 400 s，時(shí)間步長(zhǎng)為1 s，天線最大掃描范圍分別為30°，45°和60°，三種情況下的GPS衛(wèi)星網(wǎng)絡(luò)邊數(shù)據(jù)，均以txt文件格式輸出，每一行皆為三元組（u，v，t），代表一條邊，輸出數(shù)據(jù)的情況如表1所示。

若忽視衛(wèi)星與地面站之間的連接（簡(jiǎn)稱為第1種情況），GPS衛(wèi)星網(wǎng)絡(luò)將是一個(gè)獨(dú)立自主運(yùn)行的系統(tǒng)?？紤]衛(wèi)星網(wǎng)絡(luò)信息傳輸?shù)臅r(shí)效性，時(shí)間參數(shù)δ取值為1s、2s、3s（三邊連接時(shí)間點(diǎn)的最大間隔不超過δs），三節(jié)點(diǎn)三邊無向圖子圖濃度統(tǒng)計(jì)結(jié)果如圖10、圖11和圖12所示，子圖濃度計(jì)算式如式（5）所示。

既考慮星間鏈路，又考慮星地鏈路（簡(jiǎn)稱為第2種情況）時(shí)的結(jié)果如圖11、圖12所示。

由圖10、圖11和圖12可得：1）圖10、圖11和圖12所處情況下，天線最大掃描范圍都與M4子圖濃度呈正相關(guān)，αmax越大，衛(wèi)星節(jié)點(diǎn)間的連接越多，全連通結(jié)構(gòu)子圖的出現(xiàn)概率越大，M4子圖濃度越高，越有利于網(wǎng)絡(luò)的穩(wěn)定。 2）αmax=30°和60°時(shí)，第2種情況的M4子圖濃度均高于第1種情況，說明在αmax值較小時(shí)，可以借助星地鏈路建立 “衛(wèi)星－地面站－衛(wèi)星”的三節(jié)點(diǎn)全連通結(jié)構(gòu)；在αmax值較大時(shí)，衛(wèi)星節(jié)點(diǎn)之間的聯(lián)系緊密，第1種情況下的網(wǎng)絡(luò)具有較高的連通度，增加星地鏈路之后，地面站同時(shí)與多個(gè)衛(wèi)星建立連接，增加了三節(jié)點(diǎn)全連通結(jié)構(gòu)的出現(xiàn)次數(shù)，提高了M4子圖濃度。3）按照三條邊的發(fā)生順序，M1和M2的邊路徑發(fā)生了一次轉(zhuǎn)換，而M3發(fā)生了兩次轉(zhuǎn)換，發(fā)現(xiàn)在第1種情況和第2種情況的不同αmax下，M1和M2總是具有相同的子圖濃度。在本文的衛(wèi)星網(wǎng)絡(luò)模型中，受節(jié)點(diǎn)連接關(guān)系特性的影響，節(jié)點(diǎn)之間的連接在某一時(shí)間段內(nèi)持續(xù)存在，使得邊路徑發(fā)生一次轉(zhuǎn)換的M1和M2同時(shí)出現(xiàn)并具有相同的出現(xiàn)次數(shù)。

3 結(jié)語

建立衛(wèi)星網(wǎng)絡(luò)模型、構(gòu)建衛(wèi)星網(wǎng)絡(luò)數(shù)據(jù)和識(shí)別衛(wèi)星網(wǎng)絡(luò)模體，提出了衛(wèi)星網(wǎng)絡(luò)模體識(shí)別的一般方法，建立了從TLE文件輸入到子結(jié)構(gòu)識(shí)別輸出的衛(wèi)星網(wǎng)絡(luò)模體識(shí)別流程；以GPS衛(wèi)星網(wǎng)絡(luò)的三節(jié)點(diǎn)模體識(shí)別為例進(jìn)行了實(shí)證分析，結(jié)果表明：此方法基本可以實(shí)現(xiàn)衛(wèi)星網(wǎng)絡(luò)三節(jié)點(diǎn)無向圖模體識(shí)別，子圖濃度分布受模型相關(guān)參數(shù)值的影響較大，短時(shí)段內(nèi)天線最大掃描范圍與M4子圖濃度呈正相關(guān)。

本文提出的衛(wèi)星網(wǎng)絡(luò)模體識(shí)別方法，包括如下不足：建立的衛(wèi)星網(wǎng)絡(luò)模型較為簡(jiǎn)單，特別是網(wǎng)絡(luò)在建立星間鏈路和星地鏈路時(shí)，考慮的相關(guān)約束和條件較少，假設(shè)較多；受動(dòng)態(tài)網(wǎng)絡(luò)模體識(shí)別算法的影響，僅能識(shí)別三節(jié)點(diǎn)模體。下一步可以分析網(wǎng)絡(luò)遭受攻擊和損害情況下M4子圖濃度的變化情況，進(jìn)行衛(wèi)星網(wǎng)絡(luò)節(jié)點(diǎn)重要性評(píng)估和魯棒性評(píng)估等。

參考文獻(xiàn)：

[1]武健，劉新學(xué)，舒健生，等. 基于復(fù)雜網(wǎng)絡(luò)的衛(wèi)星重要度評(píng)估[J]. 火力與指揮控制， 2014， 39（5）： 60－63.

WU J， LIU X X， SHU J S， et al. Satellite significance assessment based on complex network[J]. Fire Power and Command Control， 2014， 39（5）： 60－63.

[2]朱林，方勝良，胡卿，等.衛(wèi)星時(shí)變拓?fù)渚W(wǎng)絡(luò)節(jié)點(diǎn)重要度評(píng)估方法[J].系統(tǒng)工程與電子技術(shù)，2017，39（6）：1274－1279.

ZHU L， FANG S L， HU Q， et al. Evaluation method of node significance of satellite time-varying topology network[J]. System Engineering and ElectronicTechnology，2017，39（6）：1274－1279.

[3]王瑩.衛(wèi)星移動(dòng)通信網(wǎng)若干理論和技術(shù)研究[D].武漢：華中科技大學(xué)，2008.

WANG Y.Research on theoretical and technical of satellite mobile communication network[D]. Wuhan： Huazhong University of Science and Technology，2008.

[4]林琪，李智.基于拓?fù)涮卣鞯男l(wèi)星網(wǎng)絡(luò)效能評(píng)估[J].中南大學(xué)學(xué)報(bào)（自然科學(xué)版），2013，44（S2）：368－371.

LIN Q， LI Z. Evaluation of Satellite network performance based on topological characteristics[J]. Journal of Central South University（Natural Science Edition），2013，44（S－2）：368－371.

[5]MILO R， SHEN-ORR S， ITZKOVITZ N，et al. Network motifs：simple building blocks of complex networks[J]. Science， 298（5594）：824－827.

[6]CASALS M R， Corominas-Murtra B. Assessing european power grid reliability by means of topological measures[J]. OAI， 2009，121：527－537.

[7]PAUL S， JOBST H， KURTHS J，et al. Detours around basin stability in power networks[J].New Journal of Physics，2014，16（12）：125001.

[8]孫曉偉. 基于網(wǎng)絡(luò)模體的科學(xué)學(xué)數(shù)據(jù)研究[D].成都：電子科技大學(xué)，2019.

SUN X W. Scientific data research based on networkmotif[D]. Chengdu： University of Electronic Science and Technology of China，2019.

[9]高賀. 北斗導(dǎo)航系統(tǒng)星間鏈路分配方法研究[D].湖南：湖南大學(xué)，2018.

GAO H. Research on inter-satellite link allocation method of BDS[D]. Hunan： Hunan University，2018.

[10] 辜姣，郭龍，江健，等.多層網(wǎng)絡(luò)和含時(shí)網(wǎng)絡(luò)的相關(guān)問題研究[J].復(fù)雜系統(tǒng)與復(fù)雜性科學(xué)，2016，13（1）：58－63，67.

GU J， GUO L， JIANG J， et al. Research on related problems of multilayer networks andtemporal networks[J]. Complex Systems and Complexity Science，2016，13（1）：58－63，67.

[11] HOLME P， SARAMKI J. Temporal networks[J]. Physics Reports， 2012， 519（3）：97－125.

[12] 覃桂敏.復(fù)雜網(wǎng)絡(luò)模式挖掘算法研究[D]. 西安：西安電子科技大學(xué)， 2013.

QIN G M. Research on complex network pattern mining algorithm[D]. Xi′an： Xidian University， 2013.

[13] KASHTAN N， ITZKOVITZ S， MILO R， et al. Efficient sampling algorithm for estimating subgraph concentrations and detecting network motifs.[J]. Bioinformatics，2004，20（11）：1746－1758.

[14] WERNICKE S， RASCHE F. FANMOD： a tool for fast network motif detection[J]. Bioinformatics，2006，22（9）：1152－1153.

[15] SAEED O， FALK S， MASOUDI-NEJAD， et al. MODA： an efficient algorithm for network motif discovery in biological networks[J]. Genes & Genetic Systems，2009，84（5）：385－395.

[16] TIEN H， SOMADINA M， KIM W， et al. NemoMap improved motifcentric network motif discovery algorithm[J]. Advances in Science， Technology and Engineering Systems，2018，3（5）：186－199.

[17] SCHILLER B， JAGER S， HAMACHER K， et al. Strea m-a stream-based algorithm for counting motifs in dynamic graphs[C]// International Conference on Algorithms for Computational Biology. Mexico： Springer， Cham， 2015.

[18] EDIGER D， JIANG K， RIEDY J， et al. Massive streaming data analytics： a case study with clustering coefficients[C]// IEEE International Symposium on Parallel & Distributed Processing. Atlanta， Georgia，? USA： IEEE， 2010.

[19] ASHWIN P， AUSTIN R? B， JURE L. Motifs in temporal-networks[J].CoRR，2016，abs/1612.09259.

[20] SARPE I， VANDIN F. odeN： simultaneous approximation of multiple motif counts in large temporal networks[C]. Proceedings of the 30th ACM International Conference on Information & Knowledge Management， Virtual Event. Queensland， Australia： 2021： 1568－1577.

[21] 劉天雄.GPS全球定位系統(tǒng)由幾部分組成？[J].衛(wèi)星與網(wǎng)絡(luò)，2012（4）：56－62.

LIU T X. How many parts does GPS system consist of？ [J]. Satellite and Network，2012（4）：56－62.

[22] 李朝瑞，孟新. 星座仿真中天線掃描范圍對(duì)系統(tǒng)的影響分析[C]//中國(guó)空間科學(xué)學(xué)會(huì)空間探測(cè)專業(yè)委員會(huì)第十九次學(xué)術(shù)會(huì)議論文集（下冊(cè)）.寧波， 2006，2：433－437.

LI C R， MENG X. Analysis of the influence of antenna scanning range on the system in constellation simulation[C]//Proceedings of the 19th Academic Conference of the Space Exploration Professional Committee of the Chinese Society of Space Sciences? 2. Ningbo， 2006， 2： 433－437.

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

收稿日期： 2022－04－07；修回日期：2022－06－20

基金項(xiàng)目： 復(fù)雜電子系統(tǒng)仿真重點(diǎn)實(shí)驗(yàn)室基礎(chǔ)研究項(xiàng)目（DXZT-JC-ZZ-2020-001）

第一作者： 胡博仁（1999－），男，湖南寧鄉(xiāng)人，碩士研究生，主要研究方向?yàn)閺?fù)雜網(wǎng)絡(luò)、系統(tǒng)科學(xué)。

通信作者： 裴忠民（1976－），男，山東濟(jì)寧人，博士，副研究員，主要研究方向?yàn)橛?jì)算機(jī)科學(xué)與技術(shù)、系統(tǒng)科學(xué)。