• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    神經(jīng)元網(wǎng)絡(luò)中局部同步引發(fā)的各種效應(yīng)*

    2021-09-17 06:08:32白婧關(guān)富榮唐國寧
    物理學(xué)報(bào) 2021年17期

    白婧 關(guān)富榮 唐國寧

    (廣西師范大學(xué)物理科學(xué)與技術(shù)學(xué)院, 桂林 541004)

    (廣西師范大學(xué), 廣西核物理與核技術(shù)重點(diǎn)實(shí)驗(yàn)室, 桂林 541004)

    在大腦皮層中, 神經(jīng)元大范圍的同步放電可以引發(fā)癲癇, 而癲癇發(fā)作期間可以自發(fā)出現(xiàn)螺旋波, 大量神經(jīng)元的同步放電與螺旋波自發(fā)產(chǎn)生之間的關(guān)系目前仍不清楚.本文通過增加水平長程連接構(gòu)造了具有局域長程耦合區(qū)的二維神經(jīng)元網(wǎng)絡(luò), 采用Morris-Lecar神經(jīng)元模型研究了具有多個(gè)長方形長程耦合區(qū)的神經(jīng)元網(wǎng)絡(luò)中波的傳播, 數(shù)值模擬結(jié)果表明: 傳播方向與長程耦合朝向平行的平面波和靶波經(jīng)過長程耦合區(qū)會導(dǎo)致長程耦合區(qū)內(nèi)的神經(jīng)元同步激發(fā), 這種同步激發(fā)伴隨一部分神經(jīng)元延遲激發(fā), 而另一部分提前激發(fā); 當(dāng)長程耦合區(qū)寬度超過臨界寬度時(shí), 長程耦合區(qū)所有神經(jīng)元延遲激發(fā); 當(dāng)長程耦合區(qū)寬度超過最大導(dǎo)通寬度時(shí), 波將不能通過長程耦合區(qū).當(dāng)適當(dāng)選擇長方形長程耦合區(qū)的尺寸時(shí), 神經(jīng)元同步激發(fā)可使網(wǎng)絡(luò)出現(xiàn)波回傳效應(yīng)和具有波傳播方向的選擇性, 而且這種波傳播方向的選擇性對神經(jīng)元是否處于定態(tài)和耦合強(qiáng)度變化很敏感, 以致高頻平面波列可以部分通過寬度超過最大導(dǎo)通寬度的長程耦合區(qū), 因此可以通過對長程耦合區(qū)內(nèi)的神經(jīng)元施加微擾來控制低頻波是否可以通過一定寬度的長程耦合區(qū).對于適當(dāng)選取的神經(jīng)元網(wǎng)絡(luò)結(jié)構(gòu), 當(dāng)平面波或靶波經(jīng)過長程耦合區(qū)時(shí), 網(wǎng)絡(luò)可自發(fā)出現(xiàn)自維持平面波、螺旋波和靶波等現(xiàn)象.本文對產(chǎn)生這些現(xiàn)象的物理機(jī)制作了分析.

    1 引 言

    非線性波可以在各種反應(yīng)擴(kuò)散系統(tǒng)中自發(fā)產(chǎn)生[1?3], 如心臟中的竇房結(jié)自發(fā)產(chǎn)生靶波, 當(dāng)心肌組織出現(xiàn)疤痕、波傳導(dǎo)阻塞區(qū)和心肌細(xì)胞電性質(zhì)改變時(shí), 心臟中還會自發(fā)產(chǎn)生螺旋波, 這會導(dǎo)致心動(dòng)過速, 如果螺旋波破碎成時(shí)空混沌, 還會導(dǎo)致心室纖維性顫動(dòng)而危及生命[3].Chen等[4]采用B?r-Eiswirth可激發(fā)介質(zhì)模型, 研究了平面波與缺陷相互作用, 觀察到高頻平面波通過缺陷后會自發(fā)產(chǎn)生螺旋波.Zimik和Pandit[5]采用心臟模型研究了平面波與局部異質(zhì)性團(tuán)塊(由心肌細(xì)胞和成纖維細(xì)胞組成)的相互作用, 發(fā)現(xiàn)高頻平面波通過異質(zhì)性團(tuán)塊也會自發(fā)產(chǎn)生螺旋波, 但是這些情況下螺旋波自發(fā)產(chǎn)生機(jī)制與同步無關(guān).

    非線性波(如螺旋波、行波)也可以自發(fā)出現(xiàn)在神經(jīng)元網(wǎng)絡(luò)中[6?8], 如癲癇發(fā)作時(shí)腦神經(jīng)網(wǎng)絡(luò)中可出現(xiàn)螺旋波[6,9], 癲癇發(fā)作的兩個(gè)特征是神經(jīng)元大范圍的同步振蕩和出現(xiàn)特定的節(jié)律[10], 這意味著癲癇發(fā)作時(shí)自發(fā)形成的螺旋波可能與神經(jīng)元的同步有關(guān); 在低血糖雞視網(wǎng)膜中也觀察到螺旋波自發(fā)產(chǎn)生和螺旋波導(dǎo)致的視網(wǎng)膜受損[11]的現(xiàn)象, 由于大腦神經(jīng)網(wǎng)絡(luò)中出現(xiàn)螺旋波往往與某種神經(jīng)疾病有關(guān), 掌握神經(jīng)元網(wǎng)絡(luò)中螺旋波自發(fā)產(chǎn)生的機(jī)制, 有助于有效治療相關(guān)神經(jīng)疾病, 所以螺旋波自發(fā)產(chǎn)生機(jī)制的研究受到科學(xué)工作者的極大關(guān)注, 雖然人們已經(jīng)發(fā)現(xiàn)噪聲[12]、自突觸[13]、興奮-抑制性耦合[14]、排斥耦合[15]等都可以在神經(jīng)元網(wǎng)絡(luò)自發(fā)產(chǎn)生螺旋波, 由于神經(jīng)元網(wǎng)絡(luò)的復(fù)雜性, 神經(jīng)元網(wǎng)絡(luò)中螺旋波的自發(fā)產(chǎn)生機(jī)制仍未完全清楚, 需要進(jìn)一步研究.

    心臟系統(tǒng)和大腦神經(jīng)網(wǎng)絡(luò)系統(tǒng)的結(jié)構(gòu)、工作細(xì)胞的電生理性質(zhì)、細(xì)胞之間耦合方式等都是不相同的, 但是它們有共同特點(diǎn), 都需要細(xì)胞的同步電活動(dòng), 都會出現(xiàn)螺旋波, 如心臟需要心肌細(xì)胞同步激發(fā)收縮來完成泵血功能, 很多大腦的活動(dòng)都要靠同步來實(shí)現(xiàn)[16?18], 如神經(jīng)元間信息的傳遞.研究發(fā)現(xiàn): 時(shí)滯不僅會促進(jìn)神經(jīng)元同步, 還會抑制同步[19];參數(shù)的異質(zhì)性[20]、網(wǎng)絡(luò)的拓?fù)浣Y(jié)構(gòu)[21,22]都會影響神經(jīng)元的同步.神經(jīng)元在適當(dāng)?shù)脑肼曌饔孟乱矔?shí)現(xiàn)同步[23]; 噪聲還會在神經(jīng)元網(wǎng)絡(luò)中誘發(fā)隨機(jī)共振現(xiàn)象[24,25], 部分時(shí)滯也可以在小世界神經(jīng)元網(wǎng)絡(luò)中誘發(fā)隨機(jī)多共振現(xiàn)象[26], 噪聲和延遲耦合共同作用也會使神經(jīng)元網(wǎng)絡(luò)出現(xiàn)隨機(jī)共振[27?29], 但是迄今為止, 同步在神經(jīng)元網(wǎng)絡(luò)自發(fā)出現(xiàn)螺旋波中所扮演的角色仍缺乏研究.

    目前對癲癇研究已經(jīng)取得許多研究成果, 研究結(jié)果表明[30]: 1)電突觸耦合在癲癇發(fā)作時(shí)仍然高水平地存在; 2)癲癇發(fā)作時(shí)會出現(xiàn)網(wǎng)絡(luò)重組現(xiàn)象,例如在局域癲癇中, 在癲癇剛發(fā)作時(shí)會出現(xiàn)一個(gè)大的網(wǎng)絡(luò), 在癲癇傳播中, 這個(gè)大網(wǎng)絡(luò)會破碎成更小的子網(wǎng)絡(luò), 在局部癲癇將要終止時(shí), 這些小的子網(wǎng)絡(luò)又會形成一個(gè)大網(wǎng)絡(luò); 3)在海馬切片中, 局部放電可以將局部網(wǎng)絡(luò)變成不應(yīng)區(qū), 以防止隨后的癲癇活動(dòng)在整個(gè)網(wǎng)絡(luò)中擴(kuò)散.為了透過現(xiàn)象看本質(zhì), 了解神經(jīng)元網(wǎng)絡(luò)中螺旋波的自發(fā)產(chǎn)生機(jī)制, 本文構(gòu)造一個(gè)具有局域長程耦合的二維神經(jīng)元網(wǎng)絡(luò), 它是通過在二維近鄰耦合神經(jīng)網(wǎng)絡(luò)上增加了局域水平長程連接產(chǎn)生, 而且只考慮電耦合, 然后使用Morris-Lecar神經(jīng)元模型[31]研究了平面波和靶波與局域長程耦合形成的神經(jīng)元集團(tuán)的相互作用, 數(shù)值模擬結(jié)果表明波的傳播情況與長程耦合區(qū)的朝向有關(guān).當(dāng)波的傳播方向與長程耦合區(qū)的朝向垂直時(shí), 波的傳播不受影響.當(dāng)波的傳播方向與長程耦合區(qū)的朝向平行時(shí), 波的傳播情況會受長程耦合區(qū)的影響,表現(xiàn)在: 當(dāng)平面波和靶波經(jīng)過局域長程耦合區(qū)時(shí),波與神經(jīng)元集團(tuán)相互作用導(dǎo)致神經(jīng)元集團(tuán)內(nèi)的神經(jīng)元幾乎同步激發(fā), 出現(xiàn)集團(tuán)內(nèi)一些神經(jīng)元比無長程耦合時(shí)提前激發(fā), 另一些神經(jīng)元比無長程耦合時(shí)延遲激發(fā), 當(dāng)長程耦合區(qū)足夠?qū)挄r(shí), 延遲激發(fā)時(shí)間將超過神經(jīng)元細(xì)胞從峰值回到靜息態(tài)的時(shí)間, 產(chǎn)生波回傳效應(yīng); 并且當(dāng)長程耦合區(qū)寬度超過最大導(dǎo)通寬度時(shí), 長程耦合區(qū)將阻礙波的傳播, 這時(shí)長程耦合區(qū)具有波傳播方向的選擇性; 在這兩種行為共同作用下, 當(dāng)平面波和靶波通過不同形狀和尺寸的長程耦合區(qū)時(shí)會自發(fā)產(chǎn)生平面波(或行波)、螺旋波、靶波等現(xiàn)象.這些研究結(jié)果對了解腦神經(jīng)系統(tǒng)中螺旋波產(chǎn)生的機(jī)制有積極意義.此外, 我們還發(fā)現(xiàn)長程耦合區(qū)的波傳播方向選擇性對神經(jīng)元偏離定態(tài)十分敏感, 可以利用這個(gè)性質(zhì)制作生物器件, 控制生物體中信號的傳播.下面先介紹我們的模型, 然后給出模擬結(jié)果, 最后給出結(jié)論.

    2 模 型

    本文使用一個(gè)包含 N ×N=200×200 個(gè)節(jié)點(diǎn)的二維方形神經(jīng)元網(wǎng)絡(luò), Morris-Lecar神經(jīng)元分布在網(wǎng)絡(luò)的每個(gè)節(jié)點(diǎn)上, 每個(gè)神經(jīng)元與其近鄰神經(jīng)元存在電突觸耦合, 在網(wǎng)絡(luò)中存在一塊或多塊長方形長程耦合區(qū), 在長程耦合區(qū)內(nèi)神經(jīng)元之間在水平方向與左、右(即雙邊)8個(gè)鄰居有電突觸耦合, 為了使得所有神經(jīng)元的耦合都是雙向耦合, 以及長程耦合區(qū)域內(nèi)的神經(jīng)元不能與區(qū)域外神經(jīng)元有長程耦合, 因此在該長程耦合區(qū)的邊界附近, 一個(gè)神經(jīng)元與其他位于長程耦合區(qū)的神經(jīng)元耦合的數(shù)量逐漸減少到4個(gè)(假設(shè)長程耦合區(qū)長和寬足夠大), 如圖1所示.在網(wǎng)絡(luò)只有一個(gè)長方形長程耦合區(qū)且其左和右邊界分別在 imin和 imax這一列上時(shí), 網(wǎng)絡(luò)的動(dòng)力學(xué)方程如下[31,32]:

    圖1 神經(jīng)元之間的長程耦合示意圖, 虛線代表近鄰耦合,實(shí)線代表長程耦合, 空心圓代表長程耦合區(qū)外的神經(jīng)元,實(shí)心圓代表長程耦合區(qū)內(nèi)的神經(jīng)元Fig.1.Schematic diagram of long-range coupling between neurons.Dashed and solid lines represent the nearest neighbor and long-range coupling, respectively.Hollow and solid circles represent the neurons outside and inside the longrange coupling region, respectively.

    其中 Vij分別表示網(wǎng)絡(luò)節(jié)點(diǎn) ( i,j) 上神經(jīng)元的膜電位( m V ); I 為背景電流; t 代表時(shí)間( m s ); ωi,j代表恢復(fù)變量; C 代表膜電容; gK, gCa和 gL分別表示鉀、鈣離子電流和漏電電流的最大電導(dǎo)率; VK,VCa和 VL分別表示鉀、鈣和漏電電流的反轉(zhuǎn)電位; ε 為耦合強(qiáng)度; H (ρ) 為 階躍函數(shù), 當(dāng) ρ ≥0 時(shí), H (ρ)=1 ,否則 H (ρ)=0 ; d為長程耦合區(qū)的寬度,d=imax?imin.當(dāng) j 不在長程耦合區(qū)時(shí), (1)式最后一項(xiàng)為0.本 文 各 參 量 取 值 如 下: I =39.7μA/cm2,C=20μF/cm2, gK=8μS/cm2, gCa=4μS/cm2,gL=2μS/cm2, VK=?84mV , VCa=120mV , VL=?60mV ,V1=?1.2mV , V2=18mV , V3=12mV ,V4=17.4mV , ? =0.067.在這些參數(shù)下, 單個(gè)神經(jīng)元是可激發(fā)的, 其定態(tài)為 Vs=?31.17625 , ωs=0.00694.在數(shù)值模擬中, 使用四階龍格庫塔法解微分方程,時(shí)間步長取 ? t=0.01 單位時(shí)間, 采用無流邊界條件, 取神經(jīng)元的定態(tài)為網(wǎng)絡(luò)的初態(tài).

    3 數(shù)值模擬結(jié)果

    首先在神經(jīng)元網(wǎng)絡(luò)中構(gòu)造一塊長度 L =N 個(gè)格點(diǎn)和寬度為 d 個(gè)格點(diǎn)的長方形長程耦合區(qū), 然后在系統(tǒng)左邊界 i ∈[1,10] 范圍內(nèi)通過外部刺激使神經(jīng)元瞬間處于激發(fā)態(tài), 產(chǎn)生一列向右傳播平面波,研究平面波的傳播, 接著在神經(jīng)元網(wǎng)絡(luò)中構(gòu)造一到兩個(gè)長方形長程耦合區(qū), 研究平面波與不同大小的神經(jīng)元集團(tuán)的相互作用, 最后在系統(tǒng)左邊界通過外部刺激產(chǎn)生一個(gè)向右傳播的靶波, 研究靶波與不同大小的神經(jīng)元集團(tuán)的相互作用.

    3.1 平面波的傳播與同步產(chǎn)生的效應(yīng)

    假設(shè)神經(jīng)元網(wǎng)絡(luò)中有一塊長度 L =N 個(gè)格點(diǎn)和寬度為 d 個(gè)格點(diǎn)的長方形長程耦合區(qū), 該長程耦合區(qū)左邊邊界在 imin=20 這一列上, 在平面波向右傳播過程中, 記錄在長程耦合區(qū) j =100 這一行上各格點(diǎn)的激發(fā)時(shí)刻, 這些格點(diǎn)從左到右分別記為l=1,2,···,d , 當(dāng) d =1 時(shí), 表示網(wǎng)絡(luò)中無長程耦合區(qū), 當(dāng) d =2 時(shí), 網(wǎng)絡(luò)中只有次近鄰長程耦合.圖2給出不同耦合強(qiáng)度 ε 和不同寬度 d 下神經(jīng)元的激發(fā)時(shí)刻隨神經(jīng)元位置的變化, 圖2中水平直線所包含的黑點(diǎn)數(shù)就是長程耦合區(qū)的寬度, 斜線對應(yīng)d=1(即無長程耦合), 圖2(b)是等間隔選取寬度 d 做圖.圖2(a)最下方曲線對應(yīng) d =6 , 圖2(b)最下方曲線對應(yīng) d =7.從圖2可以看出, 長程耦合區(qū)內(nèi)的神經(jīng)元激發(fā)時(shí)刻隨 d 增加而延遲, 因?yàn)楫?dāng)波傳播到長程耦合區(qū)左邊界時(shí), 與長程耦合區(qū)左邊界相鄰的神經(jīng)元處于激發(fā)態(tài), 其膜電位大大高于波峰前方與之有耦合的神經(jīng)元的膜電位, 因此這個(gè)神經(jīng)元產(chǎn)生刺激電流, 依次使波前方由近到遠(yuǎn)的神經(jīng)元偏離靜息態(tài), 從而使長程耦合區(qū)內(nèi)所有神經(jīng)元獲得幾乎相近的正刺激電流, 因?yàn)殚L程耦合區(qū)內(nèi)的神經(jīng)元耦合成一個(gè)整體, 顯然 d 越大, 各神經(jīng)元獲得的刺激電流就越小, 需要更長的時(shí)間才能激發(fā), 這就是延遲激發(fā)的原因.我們稱這個(gè)處于激發(fā)態(tài)的神經(jīng)元為產(chǎn)生刺激電流的源, 波前方的神經(jīng)元起到吸收刺激電流的作用, 稱為吸收電流的匯.顯然, 當(dāng)一個(gè)神經(jīng)元獲得的刺激電流小于其臨界值時(shí), 源與匯嚴(yán)重不匹配,神經(jīng)元不能被激發(fā), 表現(xiàn)為波不能通過長程耦合區(qū).

    圖2 不同耦合強(qiáng)度和不同寬度下神經(jīng)元激發(fā)時(shí)刻隨位置 l 的變化, 斜線對應(yīng)無長程耦合情況, 不同水平線上的黑點(diǎn)個(gè)數(shù)對應(yīng)不同寬度的長程耦合區(qū) (a) ε =0.2 ; (b)ε=0.4Fig.2.The firing time point of neurons varies with location l for different coupling strengths and different widths of the long-range coupling region.The oblique line corresponds to the network without long-range coupling, and the number of black dots on different horizontal lines is the width of the corresponding coupling region: (a) ε =0.2 ; (b) ε =0.4.

    從圖2還可以看出以下三點(diǎn).1)無長程耦合區(qū)時(shí), 神經(jīng)元的激發(fā)時(shí)刻隨神經(jīng)元位置的變化呈斜線, 表明神經(jīng)元從左到右依次激發(fā).但是有長程耦合區(qū)時(shí), 神經(jīng)元的激發(fā)時(shí)刻隨神經(jīng)元位置的變化呈水平直線, 而且耦合強(qiáng)度和長程耦合區(qū)寬度越大線越平直, 這個(gè)結(jié)果表明神經(jīng)元幾乎同步激發(fā), 水平線與斜線交點(diǎn)的左側(cè)神經(jīng)元延遲激發(fā), 交點(diǎn)右側(cè)神經(jīng)元提前激發(fā), 隨著 d 的增加, 神經(jīng)元延遲激發(fā)的比率逐漸增加.2)存在一個(gè)長程耦合區(qū)臨界寬度dth, 當(dāng) d ≥dth, 長程耦合區(qū)內(nèi)所有神經(jīng)元都會延遲激發(fā), 耦合強(qiáng)度 ε =0.2,0.4 對應(yīng)的臨界值分別是dth=23,53.可見增大神經(jīng)元之間的耦合強(qiáng)度, 可提高 dth的值.3)存在最大導(dǎo)通寬度 dmax, 當(dāng)長程耦合區(qū)寬度超過其最大導(dǎo)通寬度 dmax時(shí), 波將不能通過長程耦合區(qū), 耦合強(qiáng)度 ε =0.2,0.4 對應(yīng)長程耦合區(qū)最大導(dǎo)通寬度分別是 dmax=26,59.數(shù)值模擬發(fā)現(xiàn), 最大導(dǎo)通寬度敏感依賴神經(jīng)元之間耦合強(qiáng)度變化, 如 ε 從0.2增加到0.2005, dmax從26增加到27.

    圖2給出了波傳播方向與長程耦合朝向平行時(shí)平面波傳播情況, 如果波傳播方向與長程耦合朝向是垂直的, 我們不難得出波可以無延遲通過長程耦合區(qū).為了對長程耦合區(qū)具有對波傳播方向的選擇性有直觀印象, 圖3給出了在水平方向和豎直方向上施加一次刺激產(chǎn)生的平面波通過長程耦合區(qū)時(shí)不同時(shí)刻的斑圖, 斑圖的灰度值與神經(jīng)元的膜電位成正比(在無特別指出下, 后面的斑圖同樣處理), 白色對應(yīng)神經(jīng)元處于激發(fā)態(tài)(即波前), 深黑色區(qū)對應(yīng)神經(jīng)元處于超極化狀態(tài)(波背), 淺黑色區(qū)對應(yīng)神經(jīng)元處于靜息態(tài), 圖3(a)—(c)為水平方向的平面波通過長程耦合區(qū)的斑圖, 也就是波傳播方向與長程耦合朝向平行時(shí)平面波傳播情況, 可以看到只有一列平面波向右傳播時(shí), 在長程耦合區(qū)寬度 d ≥dmax情況下, 平面波不能通過長程耦合區(qū);圖3(d)—(f)為豎直方向的平面波通過長程耦合區(qū)的斑圖, 也就是波傳播方向與長程耦合朝向垂直時(shí)平面波傳播情況, 在長程耦合區(qū)寬度 d ≥dmax情況下, 平面波的傳播沒有受到影響, 與無長程耦合區(qū)相同, 這說明長程耦合區(qū)對波傳播方向有選擇性.

    圖3 在 ε =0.2 的情況下不同時(shí)刻的膜電位斑圖, 白色虛線矩形框內(nèi)為長程耦合區(qū), 其寬度 d =50.上下兩排圖的波傳播方向分別與長程耦合方向平行和垂直 (a), (d) t =74ms ; (b), (e) t =228ms ; (c), (f)t=556msFig.3.Membrane-potential patterns at different time points for ε =0.2.The box marked by a white dotted rectangle shows the long-range coupling region, and its width is equal to 50.The wave propagation directions in the upper and lower row panels are parallel and perpendicular to the long-range coupling direction, respectively: (a), (d) t =74ms ; (b), (e) t =228ms ; (c), (f)t=556ms.

    長程耦合區(qū)除了對波的傳播方向有選擇性外,還具有波回傳效應(yīng), 為了對此有直觀印象, 使用圖2選取的網(wǎng)絡(luò)結(jié)構(gòu)并選取不同的長程耦合區(qū)寬度, 記 錄 j =100 這 一 行 上 從 i =1 列 到 i =80 列 這

    些神經(jīng)元的膜電位隨時(shí)間的變化, 得到膜電位時(shí)空斑圖如圖4所示, 圖4(a)對應(yīng)的長程耦合區(qū)寬度為 d =6 , 可以看出, 長程耦合區(qū)左右兩邊只出現(xiàn)一條傾斜線, 表明平面波從系統(tǒng)左邊界向右傳播, 產(chǎn)生向右的傾斜線, 當(dāng)平面波傳播到長程耦合區(qū)邊界時(shí), 平面波終止在邊界處, 經(jīng)過一段延時(shí)后, 長程耦合區(qū)內(nèi)的神經(jīng)元同時(shí)被激發(fā), 由于長程耦合區(qū)左側(cè)神經(jīng)元還未回到靜息態(tài), 所以同步激發(fā)形成的波只能向右傳播, 那么這樣只出現(xiàn)向右的傾斜線;圖4(b)和圖4(c)對應(yīng)的長程耦合區(qū)寬度分別為d=14, d =21 , 在這兩張圖中可以觀察到: 長程耦合區(qū)的左邊出現(xiàn)兩條傾斜線, 右邊出現(xiàn)一條傾斜線, 長程耦合區(qū)內(nèi)的神經(jīng)元同時(shí)被激發(fā), 一方面產(chǎn)生向右傳播的波, 形成向右傾斜的直線, 另一方面產(chǎn)生向左傳播的波, 形成向左傾斜的直線, 我們把這種現(xiàn)象稱為波回傳效應(yīng).圖4(d)長程耦合區(qū)寬度為 d =27 , 此時(shí)耦合區(qū)寬度已經(jīng)超過最大導(dǎo)通寬度, 波無法通過長程耦合區(qū).可見神經(jīng)元的同步激發(fā)可導(dǎo)致長程耦合區(qū)具有波傳播方向的選擇性及神經(jīng)元網(wǎng)絡(luò)具有波回傳效應(yīng).

    圖4 在 ε =0.2 和不同的長程耦合區(qū)寬度的情況下一行格點(diǎn)的膜電位的時(shí)空斑圖, 兩白色虛線之間的區(qū)域?yàn)殚L程耦合區(qū) (a)d=6 ; (b) d =14 ; (c) d =21 ; (d)d=27Fig.4.Spatiotemporal patterns of membrane potential of a row of grid points for ε =0.2 and different widths of long-range coupling region.The region between the two white dotted lines is the long-range coupling region: (a) d =6 ; (b) d =14 ; (c) d =21 ;(d) d =27.

    通過數(shù)值模擬發(fā)現(xiàn): 長程耦合區(qū)的波傳播方向的選擇性不僅與耦合強(qiáng)度有關(guān), 還與神經(jīng)元偏離定態(tài)的幅度有關(guān), 例如在圖2(a)情況下, 將長程耦合區(qū)寬度取為 d >dmax(d=27,dmax=26) , 神經(jīng)元網(wǎng)絡(luò)的初態(tài)改為 V0=?31.09658 , ω0=0.00696 , 這個(gè)初態(tài)只比定態(tài)值高一點(diǎn), 這時(shí)平面波依然可以通過長程耦合區(qū), 我們把這種現(xiàn)象稱為波傳播方向選擇性對神經(jīng)元偏離定態(tài)幅度的敏感性, 正是這種敏感性, 導(dǎo)致周期大于等于 4 80ms 的低頻平面波不能通過長程耦合區(qū), 因?yàn)樵趫D2(a)情況下延遲時(shí)間超過 2 60ms , 神經(jīng)元?jiǎng)幼麟娢怀掷m(xù)時(shí)間約120ms, 在平面波擾動(dòng)下, 長程耦合區(qū)的神經(jīng)元可以在 4 80ms 內(nèi)回到定態(tài), 使平面波不能通過長程耦合區(qū).周期小于等于 4 70ms 平面波波列可以部分通過長程耦合區(qū), 例如外部刺激產(chǎn)生周期為450ms 的平面波, 第1至第3列波無法通過長程耦合區(qū), 在這幾列波擾動(dòng)下, 第4列波就可以通過長程耦合區(qū), 以后類似現(xiàn)象重復(fù)出現(xiàn), 每4列波通過1列波.如果外部刺激產(chǎn)生周期為 3 00ms 的平面波, 則第1列平面波不能通過長程耦合區(qū), 其產(chǎn)生的擾動(dòng)使長程耦合區(qū)內(nèi)神經(jīng)元偏離了定態(tài), 所以第2列平面波能夠通過長程耦合區(qū).但是當(dāng)?shù)?列平面波達(dá)到長程耦合區(qū)左邊界時(shí), 長程耦合區(qū)內(nèi)的神經(jīng)元?jiǎng)偤锰幱诓粦?yīng)期, 從而阻止了平面波通過,同時(shí)也為第4列平面波通過長程耦合區(qū)創(chuàng)造了條件, 使第4列波可以通過長程耦合區(qū).以后會觀察到第7, 9, 12列平面波可以通過長程耦合區(qū), 波通過長程耦合區(qū)無規(guī)律的原因是, 有些波列遇到回傳波而不能通過長程耦合區(qū).不同周期平面波通過長程耦合區(qū)情況如圖5所示.圖5是在系統(tǒng)左邊界i∈[1,10] 范圍內(nèi)通過周期外部刺激產(chǎn)生周期為 T 、向右傳播平面波得到的結(jié)果, 在無長程耦合區(qū)情況下, 每1列波通過情況用空心圓連線表示, 有長程耦合情況下平面波通過情況用黑方塊連線表示.

    圖5 在 d =27 , ε =0.2 和不同周期T平面波下神經(jīng)元的激發(fā)時(shí)刻隨位置 l 的變化.空心圓對應(yīng)無長程耦合結(jié)果, 實(shí)心方塊對應(yīng)有長程耦合的結(jié)果 (a) T =300ms ; (b) T =450ms ; (c)T=480msFig.5.The firing time point of neurons varies with location l for d =27 , ε =0.2 and different periods of planar wave.The hollow circles show the results obtained without long-range coupling, while the solid squares show the results obtained with long-range coupling: (a) T =300ms ; (b) T =450ms ; (c) T =480ms.

    對于周期更小的平面波, 主要出現(xiàn)平面波遇到長程耦合區(qū)神經(jīng)元處于不應(yīng)態(tài)和平面波遇到回傳波而不能通過長程耦合區(qū).如果增加長程耦合區(qū)的寬度 d , 只要平面波的周期足夠小, 同樣出現(xiàn)部分平面波列可以通過長程耦合區(qū).

    由于長程耦合區(qū)的波傳播方向的選擇性敏感依賴神經(jīng)元偏離定態(tài)的幅度, 我們可以通過對耦合區(qū)內(nèi)的部分神經(jīng)元施加一個(gè)微擾讓原來不能通過長程耦合區(qū)的平面波能通過該區(qū)域, 例如在(1)式右邊第一項(xiàng)方括號內(nèi)加一個(gè)刺激電流 I′=1.4μA/cm2,這個(gè)電流只施加在第20列到第24列這5列格點(diǎn)上.當(dāng)平面波傳播到長程耦合區(qū)左邊界, 第20列的神經(jīng)元偏離定態(tài)達(dá)到 ? 15.5mV 時(shí), 開始加電流刺激, 刺激持續(xù)時(shí)間為 1 0ms , 刺激結(jié)束時(shí)第20列的神經(jīng)元的膜電位只比沒有刺激電流時(shí)增加了約0.17mV, 這時(shí)平面波可以通過長程耦合區(qū).減少刺激神經(jīng)元的數(shù)量, 相應(yīng)要增加刺激電流強(qiáng)度或刺激時(shí)間, 也能使平面波通過長程耦合區(qū).這些結(jié)果表明, 可以利用長程耦合區(qū)的波傳播方向選擇性對神經(jīng)元偏離定態(tài)的幅度, 以及對耦合強(qiáng)度變化的敏感性制作生物控制器件, 通過調(diào)節(jié)耦合強(qiáng)度或?qū)ι窠?jīng)元施加擾動(dòng)控制波的通行.

    3.2 平面波和靶波與不同尺寸的神經(jīng)元集團(tuán)的相互作用

    在二維點(diǎn)陣神經(jīng)元網(wǎng)絡(luò)中構(gòu)造一個(gè)或兩個(gè)不同尺寸的長方形長程耦合區(qū), 在網(wǎng)絡(luò)中存在一個(gè)或兩個(gè)神經(jīng)元集團(tuán), 研究平面波、靶波與神經(jīng)元集團(tuán)的相互作用, 這種作用稱為波-集團(tuán)相互作用, 我們發(fā)現(xiàn)這種相互作用會在神經(jīng)元網(wǎng)絡(luò)中產(chǎn)生各種自維持有序波.為了簡單起見, 下面的研究中, 將耦合強(qiáng)度固定在 ε =0.2.

    3.2.1 平面波與神經(jīng)元集團(tuán)相互作用誘發(fā)平面波

    在二維點(diǎn)陣網(wǎng)絡(luò)構(gòu)造兩塊相同的長方形長程耦合區(qū), 長程耦合區(qū)的長度 L =N 個(gè)格點(diǎn), 寬度為d=14個(gè)格點(diǎn), 長程耦合區(qū)與系統(tǒng)左邊界平行, 兩塊區(qū)相距Δ = 13個(gè)格點(diǎn), 這樣在神經(jīng)元網(wǎng)絡(luò)中存在左右兩個(gè)神經(jīng)元集團(tuán).在系統(tǒng)左邊界 i ∈[1,10] 范圍內(nèi)通過一次刺激產(chǎn)生一列向右傳播的平面波, 當(dāng)平面波與左邊長程耦合區(qū)接觸時(shí), 平面波與左邊神經(jīng)元集團(tuán)相互作用, 由于同步激發(fā), 將產(chǎn)生向左和向右傳播的平面波, 向左傳播的平面波最終運(yùn)動(dòng)出系統(tǒng)的左邊界, 向右傳播的平面波與右邊長程耦合區(qū)的神經(jīng)元集團(tuán)相互作用, 又經(jīng)同步激發(fā)產(chǎn)生向左和向右傳播的平面波, 向右傳播的平面波最終運(yùn)動(dòng)出系統(tǒng)的右邊界, 向左傳播的波與左邊神經(jīng)元集團(tuán)相互作用, 又經(jīng)同步激發(fā)產(chǎn)生向左和向右傳播的平面波, 這種相互激發(fā)作用, 導(dǎo)致自維持平面波出現(xiàn),其產(chǎn)生過程如圖6所示.

    圖6 在 L =N , d =14 和Δ = 13的情況下不同時(shí)刻的膜電位斑圖, 左邊(右邊)兩條白色虛線之間為長程耦合區(qū) (a)t=150ms ; (b) t =292ms ; (c) t =342ms ; (d) t =370ms ; (e) t =428ms ; (f) t =478ms ; (g) t =528ms ; (h) t =570ms.Fig.6.Patterns of the membrane potential at different time points for L =N , d =14 and Δ = 13.The region between two white dotted lines on the left (right) is the long-range coupling region: (a) t =150ms ; (b) t =292ms ; (c) t =342ms ; (d) t =370ms ;(e) t =428ms ; (f) t =478ms ; (g) t =528ms ; (h) t =570ms.

    在上述平面波與集團(tuán)相互作用中, 改變長程耦合區(qū)之間的距離Δ, 原則上不會影響平面波的自發(fā)產(chǎn)生.但是當(dāng)Δ ≤ 10 時(shí), 則不能產(chǎn)生自維持平面波.原因是: 當(dāng)波回傳效應(yīng)產(chǎn)生的波進(jìn)入兩塊長程耦合區(qū)之間的空間后, 波前面的長程耦合區(qū)內(nèi)的神經(jīng)元處于靜息態(tài), 波后面的長程耦合區(qū)內(nèi)的神經(jīng)元處于超極化狀態(tài), 其膜電位更負(fù), 這時(shí)處于激發(fā)態(tài)的神經(jīng)元受到波后方處于超極化神經(jīng)元的拖累, 使其更快復(fù)極化, 導(dǎo)致其無法向波前方神經(jīng)元提供足夠長時(shí)間的電流刺激(即損失部分激發(fā)能), 從而無法使波前方處于靜息態(tài)的神經(jīng)元激發(fā), 結(jié)果兩塊長程耦合區(qū)內(nèi)的神經(jīng)元不能相互激發(fā)產(chǎn)生平面波.

    3.2.2 平面波與神經(jīng)元集團(tuán)相互作用誘發(fā)螺旋波

    本節(jié)研究平面波與小一些的神經(jīng)元集團(tuán)的相互作用.使用類似產(chǎn)生自維持平面波的網(wǎng)絡(luò)結(jié)構(gòu),適當(dāng)選取兩塊相同的長程耦合區(qū), 其長度為 L 、寬度為 d , 兩個(gè)長程耦合區(qū)之間的距離為Δ, 通過平面波與這兩個(gè)神經(jīng)元集團(tuán)相互作用可產(chǎn)生螺旋波,而且是許多小螺旋波.當(dāng)Δ比較大時(shí), 由于波回傳效應(yīng), 產(chǎn)生穩(wěn)定的螺旋波斑圖需要的時(shí)間會更長,當(dāng)Δ比較小時(shí), 可以快速得到穩(wěn)定的螺旋波斑圖,例如取 L =35 個(gè)格點(diǎn), 寬度 d =14 個(gè)格點(diǎn), 兩塊區(qū)相距Δ = 1個(gè)格點(diǎn), 則平面波與這兩個(gè)神經(jīng)元集團(tuán)相互作用會較快自發(fā)產(chǎn)生螺旋波, 其產(chǎn)生過程如圖7所示.

    從圖7可以看出, 由于延遲激發(fā)的原因, 平面波到達(dá)左邊長程耦合區(qū)的左邊界時(shí), 并不會立即導(dǎo)致該耦合區(qū)內(nèi)的神經(jīng)元激發(fā), 平面波只能沿該長程耦合區(qū)的上下邊界向右傳播, 在傳播過程中誘發(fā)耦合區(qū)上下邊界附近的神經(jīng)元激發(fā), 形成向耦合區(qū)中間傳播的波, 當(dāng)平面波即將到達(dá)該長程耦合區(qū)右邊界線時(shí), 該耦合區(qū)內(nèi)其他神經(jīng)元同步激發(fā), 一方面產(chǎn)生向前傳播的波, 導(dǎo)致右邊長程耦合區(qū)內(nèi)神經(jīng)元激發(fā), 另一方面產(chǎn)生向后傳播的行波, 該行波的兩端沿左邊耦合區(qū)邊界運(yùn)動(dòng), 在耦合區(qū)的上下邊界附近形成小螺旋波波頭, 最后逐漸形成穩(wěn)定的兩個(gè)小螺旋波.

    圖7 在 L =35 , d =14 , Δ = 1的情況下不同時(shí)刻的膜電位斑圖, 兩個(gè)白色虛線框內(nèi)為長程耦合區(qū) (a) t =150ms ;(b) t =306ms ; (c) t =318ms ; (d) t =346ms ; (e) t =442ms ; (f) t =490ms ; (g) t =514ms ; (h)t=816msFig.7.Patterns of the membrane potential at different time points for L =35 , d =14 and Δ = 1.The two boxes marked by white dotted rectangle show the long-range coupling regions: (a) t =150ms ; (b) t =306ms ; (c) t =318ms ; (d) t =346ms ;(e) t =442ms ; (f) t =490ms ; (g) t =514ms ; (h) t =816ms.

    如果保持 d =14 , 將Δ改為13, 取長程耦合區(qū)的長度分別為 L =61,101,161 , 得到不同時(shí)刻膜電位斑圖如圖8所示.從圖8可以看出, 適當(dāng)選擇長度, 都可以得到螺旋波, 只是由于 L , Δ比較大, 斑圖長時(shí)間沒有穩(wěn)定, 始終在變化, 形成的螺旋波只有波頭, 無法形成大螺旋波.圖8(a)—(f)是對稱的斑圖, 圖8(g)—(i)隨著時(shí)間的演化, 斑圖上下部分不再對稱, 這主要是由于長程耦合區(qū)上下邊界到系統(tǒng)邊界距離不相等引起的, 如果取 L =162 ,使長程耦合區(qū)上下邊界到系統(tǒng)邊界距離相等, 就可以得到對稱的斑圖, 這個(gè)結(jié)果表明, 當(dāng) L 較大時(shí)會有邊界效應(yīng).

    圖8 在 d =14 , Δ = 13和不同L下不同時(shí)刻的膜電位斑圖, 兩個(gè)白色虛線框內(nèi)為長程耦合區(qū).三個(gè)長程耦合區(qū)的長度分別為L=61,101,161 (a) t =1308ms ; (b) t =2890ms ; (c) t =3718ms ; (d) t =1370ms ; (e) t =2690ms ; (f) t =3710ms ;(g) t =1296ms ; (h) t =2400ms ; (i)t=3750msFig.8.Patterns of the membrane potential at different time points for d =14 , Δ = 13 and different L.The two boxes marked by white dotted rectangle show the long-range coupling regions.The lengths of the three long-range coupling regions are equal to 60,101, 161, respectively: (a) t =1308ms ; (b) t =2890ms ; (c) t =3718ms ; (d) t =1370ms ; (e) t =2690ms ; (f)t=3710ms;(g) t =1296ms ; (h) t =2400ms ; (i) t =3750ms.

    在固定 d =14 的情況下, 能夠產(chǎn)生螺旋波的最小 L 與兩塊長程耦合區(qū)的間隔Δ有關(guān), 當(dāng)Δ = 1時(shí), 只要 L ≥7 , 平面波經(jīng)過長程耦合區(qū)一般都可以產(chǎn)生螺旋波, 其他間隔Δ對應(yīng)的最小L值也接近7.如果網(wǎng)絡(luò)中只有一塊長方形長程耦合區(qū),適當(dāng)選擇L值也可以得到與圖8類似的小螺旋波.

    3.2.3 平面波與神經(jīng)元集團(tuán)相互作用誘發(fā)靶波

    前面我們討論了平面波與兩個(gè)神經(jīng)元集團(tuán)的相互作用, 下面我們討論網(wǎng)絡(luò)中只有一個(gè)長程耦合區(qū).適當(dāng)選取長方形長程耦合區(qū)的長度, 并且讓其寬度足夠?qū)? 則平面波與神經(jīng)元集團(tuán)相互作用可誘發(fā)靶波.這是一種由小螺旋波對為振源的靶波, 例如當(dāng)網(wǎng)絡(luò)中有一塊長度為 L =41 個(gè)格點(diǎn)、寬度為d=21個(gè)格點(diǎn)的長程耦合區(qū)時(shí), 平面波經(jīng)過該區(qū)后,網(wǎng)絡(luò)中會自發(fā)產(chǎn)生靶波, 其形成過程如圖9所示.可以看出, 由于延遲激發(fā)的原因, 平面波不能從耦合區(qū)左邊界進(jìn)入耦合區(qū), 在沿著該耦合區(qū)上下邊界向右傳播過程中, 在豎直方向形成向長程耦合區(qū)中間傳播的平面波, 平面波一端逐漸形成螺旋波波頭沿耦合區(qū)邊界運(yùn)動(dòng)(參見圖9(c)).當(dāng)運(yùn)動(dòng)到耦合區(qū)的右邊界時(shí), 螺旋波波臂在耦合區(qū)右邊界拐角處從波頭附近折斷(參見圖9(d)), 形成小螺旋波對, 最終演化成了靶波.

    圖9 在 ε =0.2 , L =41 和 d =21 的情況下不同時(shí)刻的膜電位斑圖.白色虛線框內(nèi)為長程耦合區(qū) (a) t =292ms ; (b)t=320ms ; (c) t =348ms ; (d) t =640ms ; (e) t =758ms ; (f) t =854ms ; (g) t =890ms ; (h)t=932msFig.9.Patterns of the membrane potential at different time points for ε =0.2 , L =41 and d =21.The box marked by white dotted rectangle shows the long-range coupling region: (a) t =292ms ; (b) t =320ms ; (c) t =348ms ; (d) t =640ms ; (e)t=758ms ; (f) t =854ms ; (g) t =890ms ; (h) t =932ms.

    3.3 靶波與神經(jīng)元集團(tuán)相互作用誘發(fā)有序波

    現(xiàn)在將平面波換成靶波, 仍使用圖6、圖7、圖9的網(wǎng)絡(luò)結(jié)構(gòu)研究靶波與神經(jīng)元集團(tuán)的相互作用.考慮兩種情況.第一種情況, 產(chǎn)生靶波的外部刺激施加在系統(tǒng)左邊界正中間, 覆蓋9個(gè)格點(diǎn), 其坐標(biāo)如下: i ∈[1,3] , j ∈[99,101] , 該靶波源記為S1.第二種情況, 產(chǎn)生靶波的外部刺激施加在系統(tǒng)左邊界的中間偏下方, 覆蓋9個(gè)格點(diǎn), 其坐標(biāo)如下: i ∈[1,3] ,j∈[77,79], 該靶波源記為S2.在三種網(wǎng)絡(luò)結(jié)構(gòu)下,靶波源S1產(chǎn)生的一個(gè)靶波與神經(jīng)元集團(tuán)相互作用,分別得到圖10(a)—(c)所示的結(jié)果, 靶波源S2產(chǎn)生的一個(gè)靶波與神經(jīng)元集團(tuán)相互作用, 分別得到圖10(d)—(f)所示的結(jié)果.從圖10可以看出, 改用靶波與神經(jīng)元集團(tuán)相互作用, 得到與圖6、圖7、圖9類似結(jié)果, 只是在靶波作用下不能形成平面波, 而只能形成行波, 它的波面有一定的彎曲.

    上述數(shù)值模擬結(jié)果表明: 特定結(jié)構(gòu)的神經(jīng)元網(wǎng)絡(luò)可以產(chǎn)生波回傳效應(yīng)和具有波傳播方向的選擇性, 這種網(wǎng)絡(luò)結(jié)構(gòu)極容易產(chǎn)生有序波, 有序波產(chǎn)生的本質(zhì)是波與神經(jīng)元集團(tuán)相互作用結(jié)果必須滿足神經(jīng)元的同步激發(fā)條件, 因?yàn)殚L程耦合將許多神經(jīng)元連接成一個(gè)整體.由于腦神經(jīng)網(wǎng)絡(luò)是具有可朔性的, 在病態(tài)情況下有可能出現(xiàn)本文提出的這種網(wǎng)絡(luò)結(jié)構(gòu), 使網(wǎng)絡(luò)出現(xiàn)螺旋波.

    4 結(jié) 論

    本文采用Morris-Lecar神經(jīng)元模型研究了具有特殊結(jié)構(gòu)的神經(jīng)元網(wǎng)絡(luò)中平面波的傳播和平面波、靶波與神經(jīng)元集團(tuán)的相互作用, 發(fā)現(xiàn)具有特殊結(jié)構(gòu)的長程耦合區(qū)具有波回傳效應(yīng)和波傳播方向的選擇性, 它是由于神經(jīng)元的同步激發(fā)產(chǎn)生的, 此外我們還發(fā)現(xiàn)這種長程耦合區(qū)的波傳播方向選擇性對神經(jīng)元偏離定態(tài)和耦合強(qiáng)度的變化十分敏感,當(dāng)平面波傳播方向與長程耦合方向一致時(shí), 導(dǎo)致足夠高頻率的平面波波列可部分通過寬度超過最大導(dǎo)通寬度的長程耦合區(qū).適當(dāng)選擇網(wǎng)絡(luò)結(jié)構(gòu), 當(dāng)平面波和靶波與具有一定形狀的神經(jīng)元集團(tuán)相互作用時(shí), 可以在網(wǎng)絡(luò)中自發(fā)產(chǎn)生平面波、行波、螺旋波和靶波, 同步在其中扮演了重要角色, 這些結(jié)果有助于我們理解在腦神經(jīng)網(wǎng)絡(luò)中自發(fā)出現(xiàn)的螺旋波和行波.

    av在线播放免费不卡| 黄色片一级片一级黄色片| 中出人妻视频一区二区| 一级毛片高清免费大全| 午夜两性在线视频| 亚洲国产看品久久| 99国产综合亚洲精品| av福利片在线| 久久国产精品人妻蜜桃| 97人妻天天添夜夜摸| 亚洲第一av免费看| 精品久久久久久久毛片微露脸| 波多野结衣一区麻豆| 国产成人啪精品午夜网站| 成人永久免费在线观看视频| 国产深夜福利视频在线观看| 亚洲精品美女久久av网站| 男人的好看免费观看在线视频 | 又黄又爽又免费观看的视频| 国产日韩欧美亚洲二区| 国产精品.久久久| 成人黄色视频免费在线看| 久久人妻av系列| 一级a爱视频在线免费观看| 精品亚洲成国产av| 欧美日韩中文字幕国产精品一区二区三区 | 99精品久久久久人妻精品| 中亚洲国语对白在线视频| 免费高清在线观看日韩| 国产精品九九99| 宅男免费午夜| 一区福利在线观看| 国产又色又爽无遮挡免费看| 一区二区日韩欧美中文字幕| 黄片小视频在线播放| 久久久久精品人妻al黑| 亚洲精品美女久久av网站| 色94色欧美一区二区| 成人亚洲精品一区在线观看| 午夜福利在线免费观看网站| 免费观看人在逋| 久久精品91无色码中文字幕| 久久久国产成人精品二区 | 国产深夜福利视频在线观看| 中文字幕人妻丝袜一区二区| 欧美黑人精品巨大| 亚洲情色 制服丝袜| 美国免费a级毛片| 老汉色∧v一级毛片| 国产一区有黄有色的免费视频| 久久精品国产亚洲av香蕉五月 | 国产精品久久电影中文字幕 | 免费看a级黄色片| 国产一卡二卡三卡精品| 国产免费现黄频在线看| 人成视频在线观看免费观看| 欧美日韩一级在线毛片| 伦理电影免费视频| 黄色毛片三级朝国网站| 曰老女人黄片| 在线观看www视频免费| 精品亚洲成a人片在线观看| 久久精品aⅴ一区二区三区四区| 视频区欧美日本亚洲| 夜夜爽天天搞| av片东京热男人的天堂| 亚洲精品乱久久久久久| 欧美中文综合在线视频| 丰满迷人的少妇在线观看| 国产精品一区二区免费欧美| 亚洲色图综合在线观看| 欧美乱色亚洲激情| 黑人巨大精品欧美一区二区mp4| 悠悠久久av| 亚洲成人免费av在线播放| 十八禁高潮呻吟视频| 国产一区在线观看成人免费| 国产蜜桃级精品一区二区三区 | 久久 成人 亚洲| 国产午夜精品久久久久久| 制服诱惑二区| 免费不卡黄色视频| 黄片播放在线免费| 国产精品1区2区在线观看. | 免费看十八禁软件| 亚洲欧洲精品一区二区精品久久久| 亚洲精品一卡2卡三卡4卡5卡| 久久精品国产亚洲av香蕉五月 | 人人妻,人人澡人人爽秒播| 日韩视频一区二区在线观看| 国产97色在线日韩免费| 欧美性长视频在线观看| 久久狼人影院| 亚洲一区二区三区欧美精品| 黄片播放在线免费| 三上悠亚av全集在线观看| 国产免费现黄频在线看| 国产单亲对白刺激| 欧美激情久久久久久爽电影 | 亚洲成人免费av在线播放| 狠狠狠狠99中文字幕| 校园春色视频在线观看| 亚洲精品一卡2卡三卡4卡5卡| 18禁国产床啪视频网站| 男人操女人黄网站| 成人精品一区二区免费| 久久久国产欧美日韩av| 久久99一区二区三区| 午夜成年电影在线免费观看| 亚洲欧美激情在线| 欧美日本中文国产一区发布| 精品无人区乱码1区二区| 精品久久久精品久久久| 精品视频人人做人人爽| 麻豆av在线久日| 午夜免费鲁丝| 99久久人妻综合| 天堂中文最新版在线下载| 欧美另类亚洲清纯唯美| 亚洲五月色婷婷综合| 亚洲精品美女久久久久99蜜臀| bbb黄色大片| 高清欧美精品videossex| av在线播放免费不卡| 视频在线观看一区二区三区| 国产成人av激情在线播放| 好男人电影高清在线观看| 极品教师在线免费播放| 亚洲国产中文字幕在线视频| 激情视频va一区二区三区| av有码第一页| 亚洲成国产人片在线观看| 欧美日韩av久久| 美女高潮喷水抽搐中文字幕| 欧美日韩福利视频一区二区| 精品久久久精品久久久| 精品一品国产午夜福利视频| 超色免费av| 精品国内亚洲2022精品成人 | 国产熟女午夜一区二区三区| 亚洲国产欧美网| 无遮挡黄片免费观看| 色精品久久人妻99蜜桃| 国产精品电影一区二区三区 | √禁漫天堂资源中文www| 国产亚洲一区二区精品| 视频区欧美日本亚洲| 精品欧美一区二区三区在线| 亚洲欧美激情在线| 色94色欧美一区二区| 精品午夜福利视频在线观看一区| 在线观看免费视频日本深夜| 啦啦啦 在线观看视频| 操美女的视频在线观看| av免费在线观看网站| 高清视频免费观看一区二区| 色婷婷av一区二区三区视频| 777久久人妻少妇嫩草av网站| 91精品三级在线观看| 很黄的视频免费| 成年版毛片免费区| 欧美黑人欧美精品刺激| 精品一区二区三区四区五区乱码| 成年版毛片免费区| 欧美日本中文国产一区发布| 老司机午夜福利在线观看视频| 国产亚洲精品久久久久久毛片 | 岛国在线观看网站| 在线免费观看的www视频| 日本一区二区免费在线视频| 99国产精品一区二区三区| 欧美老熟妇乱子伦牲交| 国产精品秋霞免费鲁丝片| av不卡在线播放| 在线永久观看黄色视频| 别揉我奶头~嗯~啊~动态视频| 欧美日韩国产mv在线观看视频| 国产又爽黄色视频| 嫁个100分男人电影在线观看| 一级毛片高清免费大全| 99re在线观看精品视频| 色精品久久人妻99蜜桃| 啦啦啦在线免费观看视频4| 男人舔女人的私密视频| cao死你这个sao货| 亚洲 欧美一区二区三区| 久久国产乱子伦精品免费另类| 午夜福利在线观看吧| 亚洲色图综合在线观看| 国产真人三级小视频在线观看| 中文欧美无线码| 香蕉久久夜色| 国产91精品成人一区二区三区| 美女 人体艺术 gogo| 成年版毛片免费区| 国产亚洲欧美在线一区二区| 国产男女内射视频| ponron亚洲| 亚洲七黄色美女视频| 黄色怎么调成土黄色| 国产精品久久久av美女十八| 成年人黄色毛片网站| 国产免费男女视频| 91精品三级在线观看| 精品欧美一区二区三区在线| 黄色 视频免费看| 视频区欧美日本亚洲| 在线av久久热| 欧美国产精品va在线观看不卡| 韩国av一区二区三区四区| 两个人看的免费小视频| 动漫黄色视频在线观看| 久久亚洲精品不卡| 51午夜福利影视在线观看| 久久久国产一区二区| 一级a爱视频在线免费观看| 99re6热这里在线精品视频| 老司机午夜十八禁免费视频| 免费在线观看视频国产中文字幕亚洲| 亚洲五月色婷婷综合| 亚洲精品一二三| 久9热在线精品视频| 9191精品国产免费久久| 啦啦啦免费观看视频1| 狠狠狠狠99中文字幕| 午夜久久久在线观看| 黄色成人免费大全| 在线观看一区二区三区激情| 国产av精品麻豆| 日本wwww免费看| 中文字幕高清在线视频| 日本五十路高清| 久久久久久久午夜电影 | 国产精品av久久久久免费| 日韩一卡2卡3卡4卡2021年| 久久国产亚洲av麻豆专区| 每晚都被弄得嗷嗷叫到高潮| 欧美黑人欧美精品刺激| 免费看十八禁软件| 午夜成年电影在线免费观看| 精品少妇一区二区三区视频日本电影| avwww免费| 99热国产这里只有精品6| 97人妻天天添夜夜摸| 亚洲人成77777在线视频| 国产日韩一区二区三区精品不卡| 一级a爱片免费观看的视频| 男人操女人黄网站| 亚洲精品一卡2卡三卡4卡5卡| 天天添夜夜摸| 91老司机精品| 精品一品国产午夜福利视频| 亚洲成a人片在线一区二区| 日本黄色视频三级网站网址 | 日韩中文字幕欧美一区二区| 亚洲一码二码三码区别大吗| 在线观看免费高清a一片| 亚洲精品一二三| 91麻豆av在线| 99国产极品粉嫩在线观看| 大香蕉久久成人网| 下体分泌物呈黄色| 亚洲性夜色夜夜综合| 淫妇啪啪啪对白视频| 亚洲精品久久成人aⅴ小说| 国产精品秋霞免费鲁丝片| 在线观看www视频免费| 天天操日日干夜夜撸| 9191精品国产免费久久| 久久久水蜜桃国产精品网| 欧美国产精品一级二级三级| 91成人精品电影| 一级作爱视频免费观看| 一边摸一边做爽爽视频免费| 久久久久久久精品吃奶| 女人爽到高潮嗷嗷叫在线视频| 精品人妻熟女毛片av久久网站| 夫妻午夜视频| 亚洲国产毛片av蜜桃av| 热99re8久久精品国产| 亚洲中文字幕日韩| 亚洲aⅴ乱码一区二区在线播放 | 国产又爽黄色视频| 99热网站在线观看| 天天躁夜夜躁狠狠躁躁| 国产色视频综合| av免费在线观看网站| 狠狠婷婷综合久久久久久88av| 淫妇啪啪啪对白视频| 两个人免费观看高清视频| 国产成人av教育| 国产视频一区二区在线看| 久久人妻福利社区极品人妻图片| 亚洲在线自拍视频| 黄色女人牲交| 亚洲精品国产一区二区精华液| 女人高潮潮喷娇喘18禁视频| 99热国产这里只有精品6| 亚洲精品在线美女| 18在线观看网站| 久久久久国产精品人妻aⅴ院 | 热99国产精品久久久久久7| 欧美激情极品国产一区二区三区| 两个人看的免费小视频| 热99久久久久精品小说推荐| 国产亚洲欧美精品永久| 91成人精品电影| 久久这里只有精品19| 女人被躁到高潮嗷嗷叫费观| 成年人黄色毛片网站| 久久国产精品影院| 99国产精品一区二区蜜桃av | 亚洲一码二码三码区别大吗| 欧美乱码精品一区二区三区| 亚洲精品美女久久久久99蜜臀| 搡老乐熟女国产| av在线播放免费不卡| 久久精品国产99精品国产亚洲性色 | 超色免费av| 亚洲欧美精品综合一区二区三区| 高清黄色对白视频在线免费看| 国产成+人综合+亚洲专区| 亚洲av电影在线进入| 丁香六月欧美| 国产欧美日韩精品亚洲av| 精品第一国产精品| 国产精品免费一区二区三区在线 | 波多野结衣av一区二区av| 视频区图区小说| 欧美 日韩 精品 国产| 成年人黄色毛片网站| 亚洲中文字幕日韩| 亚洲人成电影观看| 久热爱精品视频在线9| 亚洲黑人精品在线| 手机成人av网站| а√天堂www在线а√下载 | 十八禁人妻一区二区| 欧美日韩av久久| 欧美在线一区亚洲| 亚洲va日本ⅴa欧美va伊人久久| 国产精品自产拍在线观看55亚洲 | 亚洲av欧美aⅴ国产| 岛国在线观看网站| 满18在线观看网站| 欧美 亚洲 国产 日韩一| 欧美激情高清一区二区三区| 亚洲欧美一区二区三区久久| 国产一区二区三区在线臀色熟女 | 日韩欧美国产一区二区入口| 日本wwww免费看| 丝瓜视频免费看黄片| 无人区码免费观看不卡| 老熟女久久久| 日本wwww免费看| 亚洲av成人不卡在线观看播放网| 三级毛片av免费| ponron亚洲| 欧美日韩福利视频一区二区| 天天躁夜夜躁狠狠躁躁| 国产成人一区二区三区免费视频网站| 丝瓜视频免费看黄片| 美女 人体艺术 gogo| 成人国产一区最新在线观看| 视频区图区小说| 国产深夜福利视频在线观看| x7x7x7水蜜桃| 亚洲精品在线观看二区| 99久久精品国产亚洲精品| 国产在线一区二区三区精| 国产精品秋霞免费鲁丝片| 国产精品影院久久| 国产成人免费观看mmmm| 91成人精品电影| 女人被狂操c到高潮| 波多野结衣av一区二区av| 久久久久视频综合| 久久狼人影院| 男女高潮啪啪啪动态图| 精品第一国产精品| 久久精品国产清高在天天线| 久久午夜综合久久蜜桃| 久久精品国产亚洲av香蕉五月 | 亚洲久久久国产精品| 精品久久久久久久久久免费视频 | 亚洲性夜色夜夜综合| 女人爽到高潮嗷嗷叫在线视频| 欧美激情高清一区二区三区| 岛国在线观看网站| 一级,二级,三级黄色视频| 亚洲精品一卡2卡三卡4卡5卡| 中国美女看黄片| 亚洲成人免费电影在线观看| 麻豆国产av国片精品| 亚洲国产欧美日韩在线播放| 精品国产一区二区久久| 男女午夜视频在线观看| 欧美成人免费av一区二区三区 | 在线av久久热| 国产不卡一卡二| 欧美中文综合在线视频| 18禁国产床啪视频网站| 看免费av毛片| videosex国产| 亚洲专区字幕在线| 亚洲一码二码三码区别大吗| 久久国产精品影院| 香蕉国产在线看| 久久久久久人人人人人| 免费在线观看亚洲国产| 一级a爱片免费观看的视频| 国产激情久久老熟女| 亚洲视频免费观看视频| 脱女人内裤的视频| 日韩视频一区二区在线观看| www日本在线高清视频| 亚洲av欧美aⅴ国产| 亚洲精品成人av观看孕妇| 欧美精品一区二区免费开放| 欧美+亚洲+日韩+国产| 亚洲欧美激情综合另类| 国产成人系列免费观看| 高清在线国产一区| 久久久久精品人妻al黑| 热re99久久国产66热| 叶爱在线成人免费视频播放| 成年版毛片免费区| 两性夫妻黄色片| 亚洲专区中文字幕在线| 欧美丝袜亚洲另类 | videos熟女内射| 成年女人毛片免费观看观看9 | 免费不卡黄色视频| 国产精品久久久av美女十八| 国产片内射在线| 99国产精品免费福利视频| 精品久久久久久久久久免费视频 | av有码第一页| 国产亚洲精品一区二区www | 国产成人欧美| 国产国语露脸激情在线看| 成人特级黄色片久久久久久久| 中文欧美无线码| 欧美+亚洲+日韩+国产| 亚洲成国产人片在线观看| 人妻 亚洲 视频| 久久国产精品人妻蜜桃| 成人亚洲精品一区在线观看| 欧美激情久久久久久爽电影 | 51午夜福利影视在线观看| 下体分泌物呈黄色| 制服诱惑二区| 久久精品国产a三级三级三级| 成人国语在线视频| 国产精品国产av在线观看| 看免费av毛片| 午夜亚洲福利在线播放| 国产精品亚洲一级av第二区| 一区在线观看完整版| 女人高潮潮喷娇喘18禁视频| 亚洲精品国产一区二区精华液| 亚洲欧美精品综合一区二区三区| 亚洲五月天丁香| 亚洲专区国产一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲美女黄片视频| 法律面前人人平等表现在哪些方面| 亚洲男人天堂网一区| 啦啦啦在线免费观看视频4| 香蕉丝袜av| 精品福利观看| 欧美色视频一区免费| videos熟女内射| 亚洲国产中文字幕在线视频| 日本黄色日本黄色录像| 亚洲国产欧美网| 欧美在线一区亚洲| 999久久久精品免费观看国产| 国产极品粉嫩免费观看在线| 亚洲国产看品久久| 精品久久久久久久久久免费视频 | 啪啪无遮挡十八禁网站| 亚洲精品国产区一区二| 99热网站在线观看| 精品国产美女av久久久久小说| 国产成人av激情在线播放| 黑丝袜美女国产一区| 又黄又爽又免费观看的视频| 后天国语完整版免费观看| 亚洲伊人色综图| www.精华液| 免费在线观看完整版高清| 一区二区三区国产精品乱码| 亚洲av欧美aⅴ国产| 妹子高潮喷水视频| 国产精品亚洲av一区麻豆| 国产主播在线观看一区二区| 欧美最黄视频在线播放免费 | 超色免费av| 久久人人97超碰香蕉20202| 大片电影免费在线观看免费| 99re6热这里在线精品视频| 国产在视频线精品| 亚洲精华国产精华精| 国产精品综合久久久久久久免费 | 亚洲精品中文字幕一二三四区| 黑人巨大精品欧美一区二区mp4| 波多野结衣一区麻豆| 村上凉子中文字幕在线| 搡老乐熟女国产| 久久久国产精品麻豆| 亚洲午夜精品一区,二区,三区| 大香蕉久久成人网| 色在线成人网| www.熟女人妻精品国产| 亚洲人成77777在线视频| 一边摸一边抽搐一进一出视频| 一区在线观看完整版| 亚洲美女黄片视频| 搡老乐熟女国产| 成在线人永久免费视频| 变态另类成人亚洲欧美熟女 | 欧美精品一区二区免费开放| videos熟女内射| 成人18禁在线播放| 日韩熟女老妇一区二区性免费视频| 在线十欧美十亚洲十日本专区| 如日韩欧美国产精品一区二区三区| 精品人妻1区二区| 成人18禁在线播放| 老司机影院毛片| 欧美日韩精品网址| 亚洲精华国产精华精| 午夜两性在线视频| 91精品三级在线观看| 最近最新中文字幕大全免费视频| 国产成人免费无遮挡视频| 老司机影院毛片| 9191精品国产免费久久| 母亲3免费完整高清在线观看| 久久久久久久国产电影| 性少妇av在线| 亚洲欧美色中文字幕在线| 亚洲综合色网址| 亚洲国产看品久久| 国产亚洲欧美精品永久| av中文乱码字幕在线| e午夜精品久久久久久久| av天堂在线播放| 国产精品亚洲av一区麻豆| 国产蜜桃级精品一区二区三区 | avwww免费| 国产在线精品亚洲第一网站| 人人妻人人澡人人爽人人夜夜| 亚洲自偷自拍图片 自拍| 久久久久国内视频| 欧美日韩黄片免| e午夜精品久久久久久久| 日韩精品免费视频一区二区三区| 国产在线精品亚洲第一网站| 成人免费观看视频高清| 中文字幕人妻熟女乱码| videos熟女内射| av一本久久久久| 女人久久www免费人成看片| 搡老熟女国产l中国老女人| cao死你这个sao货| 老汉色∧v一级毛片| 美女 人体艺术 gogo| 精品一区二区三卡| 在线播放国产精品三级| 男女床上黄色一级片免费看| 最新美女视频免费是黄的| 国产不卡av网站在线观看| svipshipincom国产片| 热99久久久久精品小说推荐| 777米奇影视久久| 欧美精品高潮呻吟av久久| av网站在线播放免费| 午夜福利乱码中文字幕| 欧美大码av| 欧美日韩亚洲综合一区二区三区_| 天堂俺去俺来也www色官网| 久久久国产成人精品二区 | 成年人黄色毛片网站| 巨乳人妻的诱惑在线观看| 国产精品电影一区二区三区 | 欧美 日韩 精品 国产| 美女扒开内裤让男人捅视频| 国产精品.久久久| 中文亚洲av片在线观看爽 | 欧美午夜高清在线| 美国免费a级毛片| 久久午夜亚洲精品久久| 亚洲片人在线观看| 欧美乱码精品一区二区三区| 欧美精品一区二区免费开放| 久久精品熟女亚洲av麻豆精品| 成人精品一区二区免费| 男人的好看免费观看在线视频 | 老鸭窝网址在线观看| 国产在线一区二区三区精| 黄色a级毛片大全视频| 51午夜福利影视在线观看| 欧美+亚洲+日韩+国产| 国产人伦9x9x在线观看| 亚洲精品中文字幕一二三四区| 纯流量卡能插随身wifi吗| 国产精品一区二区在线不卡| 国产伦人伦偷精品视频| 一边摸一边做爽爽视频免费| 黄色片一级片一级黄色片| 亚洲成人免费av在线播放| 国产高清videossex| 久久久国产精品麻豆|