尹 寧 徐桂芝 周 茜

(河北工業(yè)大學(xué)，電磁場與電器可靠性省部共建重點實驗室，天津 300130)

(2012年12月13日收到;2013年1月26日收到修改稿)

1 引言

自然界有時并非如牛頓所說的喜歡簡單化，相反，現(xiàn)實世界中的多數(shù)系統(tǒng)都是由內(nèi)部高度非線性相互作用的動力學(xué)系統(tǒng)構(gòu)成的，具有復(fù)雜的相互關(guān)系和作用機理.復(fù)雜網(wǎng)絡(luò)模型對實際系統(tǒng)進行了合理的抽象與簡化，用節(jié)點來表示動力學(xué)子系統(tǒng)，用邊來表示子系統(tǒng)之間的相互作用，這種近似雖然忽略了子系統(tǒng)自身內(nèi)部的一些細節(jié)，但卻能夠反映出子系統(tǒng)之間的許多重要關(guān)系[1].復(fù)雜網(wǎng)絡(luò)模型可以用來描述人類社會關(guān)系、物種捕食關(guān)系、計算機網(wǎng)絡(luò)聯(lián)接、網(wǎng)絡(luò)病毒傳播、災(zāi)害蔓延、溫度變化、虛擬社區(qū)、交通系統(tǒng)、蛋白質(zhì)相互作用等[2-8].復(fù)雜網(wǎng)絡(luò)的研究已經(jīng)引起了統(tǒng)計物理學(xué)、計算科學(xué)、生態(tài)學(xué)、生命科學(xué)、社會學(xué)、經(jīng)濟學(xué)等各個領(lǐng)域的極大研究興趣，其定性和定量研究成為當今科學(xué)的一大主題.大腦是包含許多不同功能區(qū)域的復(fù)雜動力學(xué)系統(tǒng)，大腦功能的執(zhí)行離不開腦區(qū)之間的相互作用和協(xié)調(diào).因此，在動力學(xué)和復(fù)雜網(wǎng)絡(luò)理論基礎(chǔ)之上，研究跟大腦活動相關(guān)的復(fù)雜腦網(wǎng)絡(luò)是極為必要的.

目前，復(fù)雜腦網(wǎng)絡(luò)可分為反應(yīng)神經(jīng)單元解剖性連接的大腦結(jié)構(gòu)網(wǎng)絡(luò)、描述大腦統(tǒng)計性連接關(guān)系的大腦功能網(wǎng)絡(luò)以及體現(xiàn)大腦信息流向的大腦因效網(wǎng)絡(luò).大腦功能網(wǎng)絡(luò)主要反映空間上分離的不同神經(jīng)元、神經(jīng)元集群或腦區(qū)在時間上的相關(guān)性以及功能活動的統(tǒng)計依賴關(guān)系，是描述大腦動態(tài)活動交互整合的有效手段之一.大腦功能網(wǎng)絡(luò)的研究與功能磁共振成像、腦電圖和腦磁圖等成像技術(shù)密不可分，涉及心理學(xué)、認知神經(jīng)科學(xué)、病理學(xué)等諸多領(lǐng)域.Friston[9]于1994年最早開展了基于功能磁共振成像的腦功能網(wǎng)絡(luò)研究，并給出了功能連接的定義.Stam等[10]基于腦磁圖研究了阿爾茨海默病患者的大腦功能網(wǎng)絡(luò)，結(jié)果發(fā)現(xiàn)患者腦功能網(wǎng)絡(luò)的平均最短路徑長度比正常人顯著增加，“小世界”屬性有所退化.與功能核磁共振成像和腦磁圖相比，腦電圖雖然空間分辨率較低且網(wǎng)絡(luò)節(jié)點較少，但因其具有時間分辨率高、實驗成本低、易于操作等優(yōu)點[11]而被廣泛采用.Smit等[12]利用腦電圖發(fā)現(xiàn)腦功能網(wǎng)絡(luò)具有小世界拓撲結(jié)構(gòu).Li等[13]研究發(fā)現(xiàn)在完成注意任務(wù)時，腦功能網(wǎng)絡(luò)的聚類系數(shù)會高于睜眼無任務(wù)狀態(tài).Micheloyannis等[14]基于腦電圖構(gòu)建了精神分裂癥患者的腦功能網(wǎng)絡(luò)，發(fā)現(xiàn)其在靜息狀態(tài)和工作記憶任務(wù)下的局部聚類系數(shù)與正常人相比都有所下降.Fang等[15]利用腦電圖數(shù)據(jù)建立了酒精成癮癥患者的腦功能網(wǎng)絡(luò)，發(fā)現(xiàn)其復(fù)雜網(wǎng)絡(luò)特征指數(shù)與正常人有明顯不同，患者的信息熵明顯小于正常人.

穴位刺激可以等效為大腦神經(jīng)系統(tǒng)的一種外部刺激，對研究大腦復(fù)雜動力學(xué)系統(tǒng)的調(diào)控機理有重要意義，眾多學(xué)者開展了針刺穴位對大腦功能區(qū)域的影響以及相關(guān)腦功能網(wǎng)絡(luò)的研究[16,17]，從而深入分析針灸的作用機理.Li和Luo等[18,19]基于腦電圖構(gòu)建了針刺腦功能網(wǎng)絡(luò)，研究發(fā)現(xiàn)針刺足三里具有提高腦部區(qū)域間信息交流的效果，對腦疾病可能會產(chǎn)生有益的影響.與針刺相比，磁刺激技術(shù)因具有無痛、無創(chuàng)、操作簡便等優(yōu)點，正在作為一種新型的穴位刺激方式應(yīng)用于臨床研究中，在人體神經(jīng)功能調(diào)控、疾病治療與康復(fù)理療等方面具有很好應(yīng)用前景.現(xiàn)有的磁刺激穴位研究[20]主要側(cè)重單通道腦電信號的特性分析，不能深刻揭示磁刺激穴位過程中大腦各個腦區(qū)之間的相互聯(lián)系.本文采用非線性復(fù)雜網(wǎng)絡(luò)理論分析方法，首次基于腦電圖數(shù)據(jù)構(gòu)建了磁刺激穴位復(fù)雜腦功能網(wǎng)絡(luò)，從大腦的整體性角度入手研究磁刺激穴位多通道腦電信號的關(guān)聯(lián)特性和整體特征，并著重研究體表穴位磁刺激對大腦功能網(wǎng)絡(luò)拓撲結(jié)構(gòu)及特征參數(shù)的影響，為更深入的研究體表穴位磁刺激對大腦活動及神經(jīng)功能調(diào)控的影響開辟新的途徑.

2 基于腦電數(shù)據(jù)的復(fù)雜腦功能網(wǎng)絡(luò)的構(gòu)建

2.1 腦電數(shù)據(jù)采集與預(yù)處理

7名健康在校學(xué)生作為被試自愿參加了穴位磁刺激實驗，其中4名男生，3名女生，年齡集中在23—25歲之間，皆為右利手，無任何精神病史.實驗前告知各被試在整個實驗過程中保持安靜閉眼并且清醒的狀態(tài).

選取右側(cè)內(nèi)關(guān)穴作為刺激目標，該穴位與大腦神經(jīng)調(diào)節(jié)有密切關(guān)聯(lián).磁刺激實驗采用英國Magstim公司的Rapid2重復(fù)脈沖磁刺激儀.刺激線圈為8字形線圈，刺激強度為1.76 T，刺激頻率為1 Hz和3 Hz，刺激時間各120 s，刺激模式為重復(fù)刺激模式，各參數(shù)選擇均滿足安全性要求.實驗過程中保持環(huán)境安靜，室溫約為22°C，盡可能避免被試和環(huán)境受干擾.利用美國NeuroScan腦電記錄分析系統(tǒng)同步采集了刺激前、1 Hz刺激、1 Hz刺激后、3 Hz刺激、3 Hz刺激后5種狀態(tài)下的64導(dǎo)頭皮腦電信號，采樣頻率為1000 Hz.將磁刺激實驗過程中采集的腦電信號剔除明顯漂移的數(shù)據(jù)、去除眼電偽跡并進行0.5—40 Hz的帶通濾波等預(yù)處理.分別截取每個被試刺激前、1 Hz刺激、1 Hz刺激后、3 Hz刺激、3 Hz刺激后幾個階段各50 s的腦電數(shù)據(jù)，進行下一步復(fù)雜腦功能網(wǎng)絡(luò)的構(gòu)建與分析.

2.2 腦電信號關(guān)聯(lián)特性分析

互信息作為一種雙通道信號非線性分析方法，廣泛應(yīng)用于腦電分析處理領(lǐng)域，可以用來衡量任意兩兩電極測量通道間腦電信號的相互關(guān)聯(lián)關(guān)系.任意雙通道腦電信號可分別用離散隨機變量X和Y表示，它們之間的互信息可以表示為

其中，H(X)，H(Y)分別表示隨機變量X和Y的信息熵，H(XY)表示X和Y之間的聯(lián)合熵.

熵的概念最早起源于物理學(xué)，用于度量一個熱力學(xué)系統(tǒng)的無序程度.Shannon[21]借鑒了熱力學(xué)的概念將熵引入到信息論，把排除了冗余后的平均信息量稱為信息熵，解決了對信息的量化度量問題.任意雙通道腦電信號X和Y的信息熵可以分別表示為

其中，pxi和pyj分別為隨機變量X和Y的概率密度函數(shù).以離散隨機變量X為例，用直方圖法對其概率密度函數(shù)pxi進行估計，即將X分為b個區(qū)域，統(tǒng)計X落在各個區(qū)域的概率.

X和Y之間的聯(lián)合熵可以表示為

其中，pxy為隨機變量X和Y的聯(lián)合概率密度函數(shù).

通過上述方法可以分別計算每個通道腦電信號的信息熵以及任意兩個通道腦電信號之間的聯(lián)合信息熵，從而可以得到各通道腦電信號間的互信息矩陣.為統(tǒng)一度量標準，將所得的互信息進行了歸一化.歸一化互信息的大小反映了通道間信號關(guān)聯(lián)程度的高低，越接近1，表示信號間的關(guān)聯(lián)程度越高，反之表示關(guān)聯(lián)程度越低.

2.3 復(fù)雜腦功能網(wǎng)絡(luò)構(gòu)建

復(fù)雜腦功能網(wǎng)絡(luò)的構(gòu)建關(guān)鍵有三步：一是如何定義網(wǎng)絡(luò)節(jié)點，二是如何選取閾值，三是如何定義網(wǎng)絡(luò)連接邊.

1)節(jié)點選取 將每個電極通道所測量的大腦區(qū)域作為要構(gòu)建的復(fù)雜腦功能網(wǎng)絡(luò)的節(jié)點，通道腦電信號間的歸一化互信息矩陣作為節(jié)點關(guān)聯(lián)關(guān)系矩陣.

2)閾值選取 在得到完整的節(jié)點關(guān)聯(lián)關(guān)系矩陣后，需要確定一個合適的閾值來構(gòu)建二值矩陣.本文在閾值的選取時主要遵循以下原則：避免全連接、保證連接稀疏性且孤立點盡量少.當節(jié)點間的關(guān)系強度大于閾值時，二值矩陣相應(yīng)位置的元素取值為1，反之為0.

3)邊的連接 當二值矩陣元素為1時，在所對應(yīng)節(jié)點之間建立一個連接邊，反之則不建立連接邊.從而可以構(gòu)建出刺激前、1 Hz刺激、1 Hz刺激后、3 Hz刺激、3 Hz刺激后幾種狀態(tài)下的復(fù)雜腦功能網(wǎng)絡(luò).

3 基于復(fù)雜網(wǎng)絡(luò)理論的腦網(wǎng)絡(luò)分析方法

3.1 復(fù)雜腦網(wǎng)絡(luò)的特征參數(shù)

采用平均度、平均聚類系數(shù)、平均路徑長度、全局效率等復(fù)雜網(wǎng)絡(luò)測度方法對所建立的大腦功能網(wǎng)絡(luò)進行評價.通過對比研究刺激前、刺激中和刺激后的腦功能網(wǎng)絡(luò)拓撲結(jié)構(gòu)的異同，嘗試探索大腦在外界電磁刺激下的功能拓撲動態(tài)變化過程和功能活動機理.

1)平均度 度是復(fù)雜網(wǎng)絡(luò)節(jié)點的屬性中最簡單但也是最重要的性質(zhì).一個節(jié)點的度定義為與它相連的節(jié)點的數(shù)目.從直觀上看，一個節(jié)點的度越大，意味著該節(jié)點重要性越強[22].網(wǎng)絡(luò)的平均度可以通過下式計算：

其中，N表示網(wǎng)絡(luò)中的節(jié)點個數(shù)，Di表示網(wǎng)絡(luò)中第i個節(jié)點的度.

2)平均聚類系數(shù) 聚類系數(shù)是衡量網(wǎng)絡(luò)內(nèi)部聚集程度的重要參數(shù).復(fù)雜腦功能網(wǎng)絡(luò)的聚類系數(shù)反映了大腦皮層的功能性分化機理，即局部腦功能區(qū)域內(nèi)的神經(jīng)元之間連接緊密程度[22].節(jié)點i的聚類系數(shù)可以表示為

其中，Ei表示與節(jié)點i相連的ki個節(jié)點之間實際存在的邊數(shù)，Mi表示與節(jié)點i相連的ki個節(jié)點之間最多可能存在的邊數(shù).對網(wǎng)絡(luò)中所有節(jié)點的聚類系數(shù)取平均值可以得到整個網(wǎng)絡(luò)的平均聚類系數(shù)，即

3)平均路徑長度 平均路徑長度指網(wǎng)絡(luò)中任意兩個節(jié)點間最短距離的平均值，是研究復(fù)雜特性的關(guān)鍵參數(shù)[22].不考慮節(jié)點到其自身的距離對平均路徑長度的影響，腦功能網(wǎng)絡(luò)中任意節(jié)點間的平均路徑長度可由公式表示為

其中，N表示網(wǎng)絡(luò)中節(jié)點的個數(shù)，dij表示連接兩個節(jié)點i和 j的最短路徑上邊的個數(shù)，也稱為兩個節(jié)點之間的距離.

4)全局效率 全局效率是衡量網(wǎng)絡(luò)中信息傳遞快慢的一個綜合指標[23]，它可以用來評價腦電信息是如何在復(fù)雜腦功能網(wǎng)絡(luò)中快速傳遞的，具體可以表示為

其中，N表示網(wǎng)絡(luò)中節(jié)點的個數(shù)，dij表示連接兩個節(jié)點i和 j的最短路徑.

從定義上來看，平均路徑長度和全局效率都是衡量信息在腦功能網(wǎng)絡(luò)中的傳遞快慢的指標，主要用來刻畫大腦在特定狀態(tài)下能量及物質(zhì)消耗的是否“經(jīng)濟”.若腦功能網(wǎng)絡(luò)的平均路徑長度較短，即說明其全局效率較高，信息在該網(wǎng)絡(luò)中傳遞速度也較快.

3.2 復(fù)雜腦網(wǎng)絡(luò)的“小世界”屬性

由于大腦的復(fù)雜性，規(guī)則網(wǎng)絡(luò)和隨機網(wǎng)絡(luò)并不能很好地反映出腦網(wǎng)絡(luò)的特性.小世界網(wǎng)絡(luò)是介于規(guī)則網(wǎng)絡(luò)和隨機網(wǎng)絡(luò)之間的一種網(wǎng)絡(luò)模型，具有較短的平均路徑和較高的聚類系數(shù).小世界網(wǎng)絡(luò)概念的提出，為量化研究復(fù)雜腦網(wǎng)絡(luò)的性質(zhì)和腦功能性疾病提供了新手段.Watts和Strogatz建議將隨機網(wǎng)絡(luò)作為基準來量化網(wǎng)絡(luò)的“小世界”屬性[24].

如果所研究的網(wǎng)絡(luò)相對于相同規(guī)模(具有相同節(jié)點數(shù)和網(wǎng)絡(luò)平均度)的隨機網(wǎng)絡(luò)而言具有較大的聚類系數(shù)和近似的最短路徑長度，即滿足

則證明該網(wǎng)絡(luò)屬于“小世界”網(wǎng)絡(luò)范疇.其中，Creal和Lreal分別代表所構(gòu)建的真實網(wǎng)絡(luò)的平均聚類系數(shù)和平均路徑長度，Crand和Lrand分別代表隨機網(wǎng)絡(luò)的平均聚類系數(shù)和平均路徑長度.

兩個度量指標可以統(tǒng)一成一個“小世界”屬性綜合指標σ=γ/λ.當σ＞1時，網(wǎng)絡(luò)具有“小世界”屬性，且σ越大說明該網(wǎng)絡(luò)的“小世界”屬性越強[25,26].

4 結(jié)果與分析

4.1 復(fù)雜腦功能網(wǎng)絡(luò)構(gòu)建結(jié)果

閾值T的選取對大腦功能網(wǎng)絡(luò)的構(gòu)建有重要的影響，探討復(fù)雜腦網(wǎng)絡(luò)特征隨閾值的變化情況是非常必要的，作為基準態(tài)的刺激前狀態(tài)下的腦功能網(wǎng)絡(luò)特征參數(shù)隨閾值的變化曲線如圖1所示.

圖1 腦功能網(wǎng)絡(luò)特征參數(shù)隨閾值的變化曲線 (a)平均度;(b)平均聚類系數(shù);(c)平均路徑長度;(d)全局效率

由于當閾值T＜0.2時，磁刺激穴位各個狀態(tài)下的腦功能網(wǎng)絡(luò)近似為全連接狀態(tài)且特征參數(shù)基本趨于一致，不具有比較意義;而當閾值T＞0.5時，大腦功能網(wǎng)絡(luò)的孤立點大量增加，也失去了相應(yīng)的研究意義.因此，圖中給出了從閾值T=0.2逐步遞增到T=0.5過程中的特征參數(shù)變化結(jié)果.從圖中可以看出，腦功能網(wǎng)絡(luò)的平均度D、平均聚類系數(shù)C、全局效率E均隨著閾值T的逐步增大而減小，而平均路徑長度L隨閾值T的逐步增大而增大.

考慮到閾值對腦功能網(wǎng)絡(luò)特征參數(shù)的重要影響，本文在閾值選取過程中主要兼顧考慮了以下因素：盡可能保證腦功能網(wǎng)絡(luò)連接的稀疏性，盡量減少孤立點個數(shù).以被試2為例，其刺激前、1 Hz刺激、3 Hz刺激三個狀態(tài)下的腦功能網(wǎng)絡(luò)如圖2所示，其中閾值T=0.4.從圖中可以看出，1 Hz刺激和3 Hz刺激與刺激前狀態(tài)相比，大腦功能網(wǎng)絡(luò)的功能連接數(shù)有所增加，且3 Hz刺激時變化更為顯著.該閾值的選取保留了幾種狀態(tài)的差異性，同時腦功能網(wǎng)絡(luò)的連接較為稀疏且孤立點很少，滿足閾值選取要求.

4.2 磁刺激穴位對復(fù)雜腦網(wǎng)絡(luò)特征參數(shù)的影響

按照上述閾值選取原則，分別構(gòu)建了所有被試在磁刺激穴位不同狀態(tài)下的腦功能網(wǎng)絡(luò).從直觀上觀察，被試在磁刺激穴位過程中大腦的功能網(wǎng)絡(luò)的連接數(shù)量發(fā)生了不同程度的變化.為了更加深入且定量的研究磁刺激對大腦功能網(wǎng)絡(luò)產(chǎn)生的影響，分別針對不同狀態(tài)腦功能網(wǎng)絡(luò)的平均度D、平均聚類系數(shù)C、平均路徑長度L和全局效率E四個參數(shù)進行了對比分析，結(jié)果如圖3所示.從圖中可以看出，1 Hz刺激與刺激前狀態(tài)相比，大部分被試腦功能網(wǎng)絡(luò)的平均度D、平均聚類系數(shù)C和全局效率E有所升高，平均路徑長度L有所降低.3 Hz刺激與刺激前和1 Hz刺激狀態(tài)相比，所有被試腦功能網(wǎng)絡(luò)的平均度、平均聚類系數(shù)和全局效率均呈現(xiàn)增長趨勢，平均路徑長度減小，且經(jīng)過配對樣本t檢驗，結(jié)果有顯著性差異(p＜0.05).

復(fù)雜腦功能網(wǎng)絡(luò)幾個特征參數(shù)并不是孤立存在的，它們之間存在著內(nèi)在的聯(lián)系.平均聚類系數(shù)和平均路徑長度隨平均度的變化曲線如圖4所示.從圖中可以看出，平均聚類系數(shù)隨著平均度的增加呈現(xiàn)明顯的上升趨勢，而平均路徑長度則隨著平均度的增加而遞減，說明更高的邊密度導(dǎo)致更多的節(jié)點聚集到了一起，從而使節(jié)點間的距離縮短.

圖2 大腦功能網(wǎng)絡(luò)圖 (a)刺激前;(b)1 Hz刺激;(c)3 Hz刺激

圖3 不同磁刺激狀態(tài)下大腦功能網(wǎng)絡(luò)特征參數(shù)對比圖 (a)平均度;(b)平均聚類系數(shù);(c)平均路徑長度;(d)全局效率

圖4 平均聚類系數(shù)和平均路徑長度隨平均度的變化曲線

結(jié)合幾個特征參數(shù)之間的內(nèi)在聯(lián)系，以及復(fù)雜腦功能網(wǎng)絡(luò)的聚類性和結(jié)構(gòu)化分析結(jié)果表明，適當頻率的內(nèi)關(guān)穴位磁刺激會使復(fù)雜腦功能網(wǎng)絡(luò)節(jié)點間的連接數(shù)增加，網(wǎng)絡(luò)的平均度隨之增加;相鄰節(jié)點之間的聯(lián)系變得更加緊密，即平均聚類系數(shù)增加;節(jié)點之間的距離隨之縮短，于是網(wǎng)絡(luò)的平均路徑長度減小而全局效率提高.這一系列變化的結(jié)果使得復(fù)雜腦功能網(wǎng)絡(luò)區(qū)域間的內(nèi)部連通性變強，對信息交流效果產(chǎn)生一定的改善作用.

4.3 磁刺激穴位復(fù)雜腦網(wǎng)絡(luò)的小世界屬性

所構(gòu)建的腦功能網(wǎng)絡(luò)與具有相同規(guī)模(相同節(jié)點和平均度)的隨機網(wǎng)絡(luò)的拓撲結(jié)構(gòu)對比結(jié)果如圖5所示.從圖中可以直觀看出二者在網(wǎng)絡(luò)連接結(jié)構(gòu)上具有明顯的差異性.

圖6給出了磁刺激為3 Hz時大腦功能網(wǎng)絡(luò)的特征參數(shù)與相應(yīng)具有相同規(guī)模隨機網(wǎng)絡(luò)的比較結(jié)果.其中，C3和L3分別代表了3 Hz刺激腦功能網(wǎng)絡(luò)的平均聚類系數(shù)及平均路徑長度.從結(jié)果可以看出，3 Hz刺激時的腦功能網(wǎng)絡(luò)的平均聚類系數(shù)遠遠大于相應(yīng)的隨機網(wǎng)絡(luò)，而平均路徑長度與隨機網(wǎng)絡(luò)較為接近.

所有被試不同狀態(tài)下的復(fù)雜腦功能網(wǎng)絡(luò)“小世界”屬性綜合參數(shù)σ如表1所示.由結(jié)果可以看出，所有被試腦功能網(wǎng)絡(luò)的綜合參數(shù)σ均大于1，可以說明所構(gòu)建的磁刺激穴位大腦功能網(wǎng)絡(luò)均具有“小世界”屬性，即在保證有較短的平均路徑的同時又具有較高的聚類系數(shù).從結(jié)果還可以發(fā)現(xiàn)，3 Hz刺激時的綜合參數(shù)σ均高于刺激前和1 Hz刺激時的σ，表明3 Hz穴位磁刺激時復(fù)雜腦功能網(wǎng)絡(luò)的“小世界”屬性有所增強，使得網(wǎng)絡(luò)中的信息傳遞更加高效.

圖5 腦功能網(wǎng)絡(luò)與相同規(guī)模隨機網(wǎng)絡(luò)的結(jié)構(gòu)對比 (a)腦功能網(wǎng)絡(luò);(b)隨機網(wǎng)絡(luò)

圖6 大腦功能網(wǎng)絡(luò)與相同規(guī)模隨機網(wǎng)絡(luò)特征參數(shù)對比 (a)平均聚類系數(shù);(b)平均路徑長度

表1 不同狀態(tài)下的腦功能網(wǎng)絡(luò)“小世界”屬性綜合參數(shù)σ

5 結(jié)論

本文基于非線性和復(fù)雜網(wǎng)絡(luò)分析方法，從整體性角度研究了磁刺激內(nèi)關(guān)穴對大腦功能連接產(chǎn)生的影響.利用互信息方法對多通道腦電信號進行了非線性時域關(guān)聯(lián)特性分析，并在此基礎(chǔ)上構(gòu)建了磁刺激穴位不同狀態(tài)下的大腦功能網(wǎng)絡(luò).利用復(fù)雜網(wǎng)絡(luò)測度方法對靜息態(tài)以及不同刺激態(tài)下的大腦功能網(wǎng)絡(luò)結(jié)構(gòu)分布、特征參數(shù)、小世界屬性進行了深入研究與分析，探索了體表穴位磁刺激時大腦功能網(wǎng)絡(luò)拓撲結(jié)構(gòu)的動態(tài)變化.結(jié)果發(fā)現(xiàn)，3 Hz磁刺激與刺激前相比，復(fù)雜腦功能網(wǎng)絡(luò)的平均度、平均聚類系數(shù)和全局效率均有顯著升高，平均路徑長度顯著降低.同時，3 Hz刺激時腦功能網(wǎng)絡(luò)的“小世界”屬性也有所增強，信息在大腦區(qū)域間的傳遞更加高效.

本文首次開展了體表穴位磁刺激復(fù)雜腦功能網(wǎng)絡(luò)的分析和研究，將各個腦區(qū)作為整體來探索磁刺激穴位對大腦產(chǎn)生的影響，為磁刺激穴位的研究開拓了新思路和新方法.研究結(jié)果表明，適當頻率的內(nèi)關(guān)穴位磁刺激使得復(fù)雜腦功能網(wǎng)絡(luò)的節(jié)點連接有所增加，相鄰節(jié)點之間的聯(lián)系更加緊密，節(jié)點之間的距離縮短，網(wǎng)絡(luò)區(qū)域間的連通性增強，信息交流效果有所改善，對一些大腦功能性疾病可能會產(chǎn)生有益的影響.磁刺激穴位復(fù)雜腦功能網(wǎng)絡(luò)研究為腦功能性疾病的康復(fù)理療提供了新的思路，值得繼續(xù)開展大量實驗和理論研究進行深入探討.

[1]Zhang F F，Zheng Z G 2012 J.ShanghaiUniv.Sci.Tech.34 138(in Chinese)[張方風，鄭志剛2012上海理工大學(xué)學(xué)報34 138]

[2]Kaiser-Bunbury C N，Muff S，Memmott J，Mller C B，Ca flisch A 2010 Ecol.Lett.13 442

[3]Ouyang M，F(xiàn)eiQ，Yu M H 2008 Acta Phys.Sin.57 6763(in Chinese)[歐陽敏，費奇，余明暉2008物理學(xué)報57 6763]

[4]Zhou L，Gong Z Q，ZhiR，F(xiàn)eng G L 2008 Acta Phys.Sin.57 7380(in Chinese)[周磊，龔志強，支蓉，封國林2008物理學(xué)報57 7380]

[5]Zhang L，Liu Y 2008 Acta Phys.Sin.57 5419(in Chinese)[張立，劉云2008物理學(xué)報57 5419]

[6]Zheng X，Chen J P，ShaoJ L，Bie L D 2012 Acta Phys.Sin.61 190510(in Chinese)[鄭嘯，陳建平，邵佳麗，別立東2012物理學(xué)報 61 190510]

[7]Liu H K，Zhou T 2007 Acta Phys.Sin.56 106(in Chinese)[劉宏鯤，周濤2007物理學(xué)報56 106]

[8]ColomboE，Alcalay M，PelicciP G 2011 Oncogene 30 2595

[9]Friston K J 1994 Hum.Brain Mapp.2 56

[10]Stam C J，de Haan W，Daffertshofer A，Jones B F，Manshanden I，van Cappellen van Walsum A M，Montez T，Verbunt J P A，de Munck J C，van Dijk B W，Berendse H W，Scheltens P 2009 Brain 132 213

[11]LiL，Jin Z L，LiB 2011 Acta Phys.Sin.60 048703(in Chinese)[李凌，金貞蘭，李斌2011物理學(xué)報60 048703]

[12]Smit D J A，Stam C J，Posthuma D，Boomsma D I，de Geus E J C 2008 Hum.Brain Mapp.29 1368

[13]LiL，Jin Z L，LiB 2011 Chin.Phys.B 20 038701

[14]Micheloyannis S，Pachou E，Stam C J，Breakspear M，Bitsios P，Vourkas M，ErimakiS，Zervakis M 2006 Schizophr.Res.87 60

[15]Fang X L，Jiang Z L 2007 Acta Phys.Sin.56 7330(in Chinese)[方小玲，姜宗來2007物理學(xué)報56 7330]

[16]Bian R J，Wang J，Han C X，Deng B，WeiX L，Che Y Q 2011 Acta Phys.Sin.60 118701(in Chinese)[邊洪瑞，王江，韓春曉，鄧斌，魏熙樂，車艷秋2011物理學(xué)報60 118701]

[17]Han C X，Wang J，Che Y Q，Deng B，GuoYi，GuoY M，Liu Y Y 2010 Acta Phys.Sin.59 5881(in Chinese)[韓春曉，王江，車艷秋，鄧斌，郭義，郭永明，劉陽陽2010物理學(xué)報59 5881]

[18]LiN，Wang J，Deng B，WeiX L，Che Y Q，Jia C H，GuoY，Wang C 2011 Acup.Res.36 278(in Chinese)[李諾，王江，鄧斌，魏熙樂，車艷秋，賈晨輝，郭義，王超2011針刺研究36 278]

[19]LuoX L，Wang J，Han C X，Deng B，WeiX L，Bian H R 2012 Chin.Phys.B 21 028701

[20]Xu G Z，Zhang X，Yang Q X，Yang S，Yan W L 2010 Trans.Chin.Electrotech.Soc.25 34(in Chinese)[徐桂芝，張秀，楊慶新，楊碩，顏威利2010電工技術(shù)學(xué)報25 34]

[21]Shannon C E 1949 BellSyst.Tech.J.28 656

[22]Albert R，BarabasiA L 2002 Rev.Mod.Phys.74 47

[23]Latora V，MarchioriM 2001 Phys.Rev.Lett.87 198701

[24]Watts D J，Strogatz S H 1998 Nature 393 440

[25]Humphries M D，Gurney K，Prescott T J 2006 Proc.R.Soc.Lond.B 273 503

[26]SuiD N，ZhaoQ B，Tang Y Y 2010 Chin.Phys.Lett.27 018702