• <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蜜桃av| 亚洲精品国产区一区二| 国产精品亚洲一级av第二区| 777久久人妻少妇嫩草av网站| 精品久久久精品久久久| 一区二区三区国产精品乱码| 淫妇啪啪啪对白视频| 成人影院久久| 黄色成人免费大全| 国产精品国产av在线观看| 欧美在线一区亚洲| 极品人妻少妇av视频| 丁香六月欧美| 国产深夜福利视频在线观看| 国产欧美日韩精品亚洲av| 18禁裸乳无遮挡免费网站照片 | 免费在线观看完整版高清| 交换朋友夫妻互换小说| 狂野欧美激情性xxxx| 精品国产国语对白av| 身体一侧抽搐| 99在线人妻在线中文字幕| 久久性视频一级片| 两个人看的免费小视频| 老司机靠b影院| 欧美黄色片欧美黄色片| 欧美日韩福利视频一区二区| 丝袜美足系列| 婷婷精品国产亚洲av在线| 国产一区二区激情短视频| 成人黄色视频免费在线看| 丰满迷人的少妇在线观看| 亚洲五月色婷婷综合| 婷婷六月久久综合丁香| 在线观看免费日韩欧美大片| 怎么达到女性高潮| 国产精品综合久久久久久久免费 | 亚洲精品国产精品久久久不卡| 十八禁人妻一区二区| 桃色一区二区三区在线观看| 国产高清videossex| 午夜福利影视在线免费观看| 国产精品国产av在线观看| 久久久久久久精品吃奶| 高清av免费在线| 国产精品久久久久成人av| 91成人精品电影| 国产精品自产拍在线观看55亚洲| 老司机午夜福利在线观看视频| 国产熟女xx| 午夜免费成人在线视频| 露出奶头的视频| 中文字幕人妻熟女乱码| 亚洲午夜理论影院| 亚洲自拍偷在线| 99精国产麻豆久久婷婷| 十分钟在线观看高清视频www| 国产高清国产精品国产三级| 麻豆久久精品国产亚洲av | 国产1区2区3区精品| 亚洲成人久久性| 成年人免费黄色播放视频| 久久久久九九精品影院| 亚洲成国产人片在线观看| 中亚洲国语对白在线视频| 黄色女人牲交| 亚洲人成网站在线播放欧美日韩| 亚洲av片天天在线观看| 777久久人妻少妇嫩草av网站| 亚洲精品成人av观看孕妇| 一级黄色大片毛片| 国产视频一区二区在线看| 黄色 视频免费看| 99久久人妻综合| 91字幕亚洲| 国产精品久久久av美女十八| 欧美激情极品国产一区二区三区| 99热国产这里只有精品6| 国产1区2区3区精品| 好男人电影高清在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 精品福利观看| 亚洲成av片中文字幕在线观看| 欧美成狂野欧美在线观看| 国产在线精品亚洲第一网站| 老司机靠b影院| 久久久久久久午夜电影 | 精品国产一区二区三区四区第35| av超薄肉色丝袜交足视频| 日韩欧美在线二视频| 国产精品一区二区精品视频观看| 天堂动漫精品| 国产欧美日韩综合在线一区二区| 黑人操中国人逼视频| 男女午夜视频在线观看| 亚洲精华国产精华精| 黄色丝袜av网址大全| 久久精品91蜜桃| 人人妻人人爽人人添夜夜欢视频| 免费在线观看亚洲国产| 黄色成人免费大全| 亚洲精品美女久久av网站| 国产国语露脸激情在线看| 99久久99久久久精品蜜桃| 欧美日本亚洲视频在线播放| 亚洲一码二码三码区别大吗| 国产aⅴ精品一区二区三区波| 超碰97精品在线观看| 日本黄色日本黄色录像| 桃色一区二区三区在线观看| 亚洲午夜精品一区,二区,三区| 国产又色又爽无遮挡免费看| 久久九九热精品免费| 一级毛片精品| 在线视频色国产色| 十八禁人妻一区二区| 免费在线观看日本一区| 欧美+亚洲+日韩+国产| 久久中文字幕人妻熟女| 老司机福利观看| 在线观看66精品国产| 日韩欧美三级三区| 在线观看免费午夜福利视频| 亚洲一区中文字幕在线| 免费在线观看完整版高清| 日韩有码中文字幕| www日本在线高清视频| 久热这里只有精品99| 桃红色精品国产亚洲av| 免费一级毛片在线播放高清视频 | 少妇裸体淫交视频免费看高清 | 波多野结衣av一区二区av| 国产成人精品无人区| 伦理电影免费视频| 欧美日韩亚洲高清精品| 在线看a的网站| 久久久久精品国产欧美久久久| 丰满饥渴人妻一区二区三| 亚洲熟妇熟女久久| 亚洲avbb在线观看| 一进一出抽搐gif免费好疼 | 丝袜在线中文字幕| 久久99一区二区三区| 国产片内射在线| 老汉色∧v一级毛片| 成人亚洲精品一区在线观看| 欧洲精品卡2卡3卡4卡5卡区| 日韩免费av在线播放| 中文字幕另类日韩欧美亚洲嫩草| 久久影院123| 国产欧美日韩一区二区三区在线| 久久亚洲真实| 在线观看舔阴道视频| 日韩 欧美 亚洲 中文字幕| 午夜福利在线观看吧| 窝窝影院91人妻| 国产亚洲精品第一综合不卡| 超色免费av| 麻豆一二三区av精品| 国产在线观看jvid| 别揉我奶头~嗯~啊~动态视频| 搡老熟女国产l中国老女人| 男人舔女人下体高潮全视频| 别揉我奶头~嗯~啊~动态视频| 国产伦一二天堂av在线观看| 老汉色av国产亚洲站长工具| 高清av免费在线| 久久香蕉精品热| 91麻豆精品激情在线观看国产 | 亚洲成人免费av在线播放| 亚洲三区欧美一区| 黑人操中国人逼视频| 午夜两性在线视频| 免费在线观看完整版高清| 一级,二级,三级黄色视频| 老司机福利观看| 搡老熟女国产l中国老女人| 男女午夜视频在线观看| 精品国产一区二区三区四区第35| 叶爱在线成人免费视频播放| 在线观看免费午夜福利视频| 久久热在线av| 久久香蕉精品热| 精品人妻在线不人妻| 一夜夜www| 99久久国产精品久久久| 久久久久久久久免费视频了| 国产精品98久久久久久宅男小说| 色在线成人网| 国产精品爽爽va在线观看网站 | 亚洲专区国产一区二区| 一a级毛片在线观看| 久久久久久久午夜电影 | 久久久国产精品麻豆| 十八禁网站免费在线| 一边摸一边抽搐一进一小说| 人人妻人人澡人人看| 国产日韩一区二区三区精品不卡| 9热在线视频观看99| 欧美av亚洲av综合av国产av| 久久久国产精品麻豆| 波多野结衣一区麻豆| 级片在线观看| 91在线观看av| 日韩免费av在线播放| 亚洲熟妇熟女久久| 国产色视频综合| 成人18禁高潮啪啪吃奶动态图| 在线观看66精品国产| 在线永久观看黄色视频| 夜夜爽天天搞| 精品国产国语对白av| 亚洲精品一二三| 国产精品秋霞免费鲁丝片| 少妇 在线观看| 午夜福利免费观看在线| 国产人伦9x9x在线观看| 国产xxxxx性猛交| 高清毛片免费观看视频网站 | 久久中文看片网| 亚洲人成伊人成综合网2020| 国产一区二区三区视频了| 少妇裸体淫交视频免费看高清 | 亚洲视频免费观看视频| av片东京热男人的天堂| 悠悠久久av| 亚洲国产精品一区二区三区在线| 在线国产一区二区在线| 99在线人妻在线中文字幕| 精品久久久久久成人av| 涩涩av久久男人的天堂| 黄色片一级片一级黄色片| 欧美日韩av久久| 国产精品日韩av在线免费观看 | 亚洲精品一卡2卡三卡4卡5卡| 淫秽高清视频在线观看| 国产精品电影一区二区三区| av在线播放免费不卡| 国产精品美女特级片免费视频播放器 | 女同久久另类99精品国产91| 免费女性裸体啪啪无遮挡网站| av天堂久久9| 亚洲成人免费av在线播放| 国产精华一区二区三区| 九色亚洲精品在线播放| 男女之事视频高清在线观看| 久久久久久久久中文| 欧美老熟妇乱子伦牲交| 熟女少妇亚洲综合色aaa.| 亚洲精品国产区一区二| 亚洲人成电影观看| 亚洲色图av天堂| 久久精品人人爽人人爽视色| 99精品在免费线老司机午夜| 日韩欧美一区视频在线观看| 精品欧美一区二区三区在线| 色在线成人网| 色婷婷av一区二区三区视频| 三上悠亚av全集在线观看| 午夜福利在线观看吧| 国产成+人综合+亚洲专区| 中文字幕人妻丝袜一区二区| 99久久综合精品五月天人人| 日韩欧美国产一区二区入口| 免费在线观看黄色视频的| 国产亚洲av高清不卡| 男女之事视频高清在线观看| 亚洲中文日韩欧美视频| 国产午夜精品久久久久久| 欧美中文日本在线观看视频| 9热在线视频观看99| 亚洲性夜色夜夜综合| 50天的宝宝边吃奶边哭怎么回事| 他把我摸到了高潮在线观看| 国产亚洲欧美精品永久| 国产在线观看jvid| 午夜老司机福利片| 女人被躁到高潮嗷嗷叫费观| 欧美色视频一区免费| 水蜜桃什么品种好| 啪啪无遮挡十八禁网站| 日韩精品中文字幕看吧| 黄频高清免费视频| 一级毛片高清免费大全| 亚洲精品粉嫩美女一区| 看免费av毛片| 日本精品一区二区三区蜜桃| 国产大屁股一区二区在线视频| 精华霜和精华液先用哪个| 国产麻豆成人av免费视频| 国产午夜福利久久久久久| 性插视频无遮挡在线免费观看| 欧美+日韩+精品| 亚洲内射少妇av| 一本精品99久久精品77| 亚洲av.av天堂| 校园春色视频在线观看| 日韩欧美精品v在线| 国产麻豆成人av免费视频| 不卡一级毛片| 国产亚洲精品久久久com| 国产乱人伦免费视频| 免费看美女性在线毛片视频| 亚洲国产日韩欧美精品在线观看| 精品久久国产蜜桃| 国产精品野战在线观看| 亚洲av免费在线观看| 蜜桃久久精品国产亚洲av| 美女大奶头视频| 精品无人区乱码1区二区| 两个人的视频大全免费| 一本一本综合久久| 国产一区二区三区视频了| 午夜老司机福利剧场| 久久久久国产精品人妻aⅴ院| 亚洲av电影在线进入| 黄色视频,在线免费观看| 啦啦啦韩国在线观看视频| 久久精品国产清高在天天线| www.熟女人妻精品国产| 久久这里只有精品中国| 欧美国产日韩亚洲一区| 亚洲综合色惰| av黄色大香蕉| 国产欧美日韩一区二区三| 色av中文字幕| 男人舔女人下体高潮全视频| 欧美高清性xxxxhd video| 深夜a级毛片| 最新在线观看一区二区三区| 国产av在哪里看| 少妇人妻精品综合一区二区 | 性插视频无遮挡在线免费观看| 亚洲欧美日韩东京热| 亚洲av一区综合| 中文字幕av在线有码专区| АⅤ资源中文在线天堂| av专区在线播放| 国产欧美日韩精品亚洲av| 99热这里只有是精品50| 欧美黑人欧美精品刺激| 看免费av毛片| 国产亚洲欧美在线一区二区| 国产欧美日韩精品一区二区| 九九热线精品视视频播放| 免费人成在线观看视频色| www.999成人在线观看| 观看免费一级毛片| 波多野结衣高清作品| 国产成人a区在线观看| 香蕉av资源在线| 中文字幕人妻熟人妻熟丝袜美| 久久婷婷人人爽人人干人人爱| 免费黄网站久久成人精品 | 看免费av毛片| 亚洲av日韩精品久久久久久密| 国产精品野战在线观看| 一个人免费在线观看的高清视频| 久久久久免费精品人妻一区二区| 一个人免费在线观看电影| 欧美激情国产日韩精品一区| 成人午夜高清在线视频| 黄色日韩在线| 在线观看午夜福利视频| 亚洲精品在线美女| 久9热在线精品视频| 免费在线观看影片大全网站| 亚洲人与动物交配视频| 最好的美女福利视频网| 首页视频小说图片口味搜索| 91av网一区二区| 国产v大片淫在线免费观看| 亚洲成人中文字幕在线播放| 免费黄网站久久成人精品 | 日日摸夜夜添夜夜添小说| 久久人人爽人人爽人人片va | 免费电影在线观看免费观看| 人妻丰满熟妇av一区二区三区| 99热只有精品国产| 成人欧美大片| 90打野战视频偷拍视频| 欧美激情国产日韩精品一区| 丝袜美腿在线中文| .国产精品久久| 97热精品久久久久久| 九九在线视频观看精品| 久9热在线精品视频| 欧美一区二区精品小视频在线| 亚洲成人久久爱视频| 97碰自拍视频| 夜夜夜夜夜久久久久| 人人妻人人澡欧美一区二区| 亚洲精品456在线播放app | 亚洲国产精品成人综合色| 国产黄片美女视频| 禁无遮挡网站| 天天一区二区日本电影三级| 综合色av麻豆| 狂野欧美白嫩少妇大欣赏| 日韩 亚洲 欧美在线| 日日摸夜夜添夜夜添小说| 日韩欧美精品v在线| 亚洲精华国产精华精| a级一级毛片免费在线观看| 日本精品一区二区三区蜜桃| 成人亚洲精品av一区二区| 五月伊人婷婷丁香| 午夜免费激情av| 少妇熟女aⅴ在线视频| 嫩草影院精品99| 国产伦精品一区二区三区四那| 草草在线视频免费看| 欧美成人免费av一区二区三区| 日韩人妻高清精品专区| 97热精品久久久久久| 国语自产精品视频在线第100页| 国产精品自产拍在线观看55亚洲| a级一级毛片免费在线观看| 国产私拍福利视频在线观看| 亚洲avbb在线观看| 亚洲欧美日韩无卡精品| 午夜激情欧美在线| 悠悠久久av| 久久九九热精品免费| 三级男女做爰猛烈吃奶摸视频| 在线观看美女被高潮喷水网站 | 国产白丝娇喘喷水9色精品| 美女xxoo啪啪120秒动态图 | 久久久久久国产a免费观看| 精品久久久久久久久久免费视频| 久久精品国产亚洲av香蕉五月| 免费看日本二区| 波多野结衣高清无吗| 亚洲av免费高清在线观看| 在线观看免费视频日本深夜| 久久久久久久精品吃奶| 搡老熟女国产l中国老女人| 九九热线精品视视频播放| 脱女人内裤的视频| 男女做爰动态图高潮gif福利片| 成年女人毛片免费观看观看9| 国产乱人视频| 十八禁人妻一区二区| 亚洲色图av天堂| 最后的刺客免费高清国语| 亚洲第一欧美日韩一区二区三区| av在线蜜桃| a在线观看视频网站| 黄色配什么色好看| 黄色视频,在线免费观看| 窝窝影院91人妻| www.色视频.com| 精品久久久久久成人av| 亚洲国产精品sss在线观看| 一级黄片播放器| 悠悠久久av| 亚洲av免费高清在线观看| 国产视频一区二区在线看| 日本 欧美在线| 成人国产综合亚洲| 搡老妇女老女人老熟妇| 亚洲中文日韩欧美视频| 九色成人免费人妻av| 又黄又爽又刺激的免费视频.| 日本 av在线| av黄色大香蕉| 色综合站精品国产| 99久久九九国产精品国产免费| 国模一区二区三区四区视频| 我的老师免费观看完整版| 99久久精品国产亚洲精品| 午夜福利欧美成人| 男女那种视频在线观看| 极品教师在线免费播放| 变态另类成人亚洲欧美熟女| 麻豆成人午夜福利视频| 亚洲美女黄片视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 看十八女毛片水多多多| 1024手机看黄色片| 成人三级黄色视频| 国产aⅴ精品一区二区三区波| 色综合站精品国产| 看十八女毛片水多多多| 欧美日韩福利视频一区二区| 亚洲精品成人久久久久久| 99久久成人亚洲精品观看| 精品久久久久久久久av| 99国产精品一区二区三区| 一本久久中文字幕| 国产男靠女视频免费网站| 直男gayav资源| 亚洲精品粉嫩美女一区| 欧美激情久久久久久爽电影| 1024手机看黄色片| 亚洲精品成人久久久久久| 日本一本二区三区精品| 国内毛片毛片毛片毛片毛片| 亚洲成人免费电影在线观看| 国产视频内射| 高清日韩中文字幕在线| 亚洲中文日韩欧美视频| 嫩草影院入口| 亚洲精品亚洲一区二区| 日本三级黄在线观看| 色哟哟·www| 日韩精品中文字幕看吧| 别揉我奶头~嗯~啊~动态视频| 亚洲国产精品999在线| 九色国产91popny在线| 国产一区二区在线av高清观看| 国产精品99久久久久久久久| 性色avwww在线观看| 91字幕亚洲| 欧美高清性xxxxhd video| 亚洲内射少妇av| 黄色配什么色好看| 成人性生交大片免费视频hd| 婷婷色综合大香蕉| 高潮久久久久久久久久久不卡| 欧美日韩乱码在线| 精品国产亚洲在线| 美女高潮的动态| 在线免费观看的www视频| 在线观看舔阴道视频| 国产成人影院久久av| 亚洲中文字幕一区二区三区有码在线看| 麻豆av噜噜一区二区三区| www.色视频.com| 99国产精品一区二区蜜桃av| 国产精品不卡视频一区二区 | 亚洲一区二区三区色噜噜| 国产单亲对白刺激| 国产主播在线观看一区二区| 中出人妻视频一区二区| 亚洲av成人不卡在线观看播放网| 午夜影院日韩av| 久久九九热精品免费| 亚洲第一区二区三区不卡| 夜夜爽天天搞| 亚洲在线观看片| 久久这里只有精品中国| 亚洲avbb在线观看| www日本黄色视频网| 中文字幕免费在线视频6| 日韩欧美一区二区三区在线观看| 国产在视频线在精品| 国产精品久久久久久精品电影| 男女那种视频在线观看| 精品乱码久久久久久99久播| 久久久久久久亚洲中文字幕 | 一本一本综合久久| 国产精品综合久久久久久久免费| 非洲黑人性xxxx精品又粗又长| 最好的美女福利视频网| 国产v大片淫在线免费观看| 99久久精品热视频| 久久人妻av系列| 国产av不卡久久| 亚洲午夜理论影院| 日日干狠狠操夜夜爽| av福利片在线观看| 桃红色精品国产亚洲av| 国产高清视频在线播放一区| 99热6这里只有精品| 亚洲成人中文字幕在线播放| 脱女人内裤的视频| 日本 欧美在线| av在线观看视频网站免费| 久久久久久久亚洲中文字幕 | 国产精品一区二区性色av| 欧美极品一区二区三区四区| 亚洲无线在线观看| 小说图片视频综合网站| 久久99热6这里只有精品| 久久久久久久久久黄片| 国产成+人综合+亚洲专区| 久久久国产成人免费| 在线观看一区二区三区| 日日夜夜操网爽| 一区二区三区四区激情视频 | 女人被狂操c到高潮| 国产av不卡久久| 老司机深夜福利视频在线观看| 国产精品嫩草影院av在线观看 | 国产精品三级大全| 精品不卡国产一区二区三区| 一级av片app| 少妇的逼好多水| 麻豆成人av在线观看| 国产精品伦人一区二区| 变态另类丝袜制服| 女同久久另类99精品国产91| 久久草成人影院| 色综合站精品国产| 亚洲电影在线观看av| 一卡2卡三卡四卡精品乱码亚洲| 淫秽高清视频在线观看| 999久久久精品免费观看国产| 精品午夜福利视频在线观看一区| 久久人人爽人人爽人人片va | 女人十人毛片免费观看3o分钟| 亚洲精品亚洲一区二区| 十八禁国产超污无遮挡网站| 99国产精品一区二区蜜桃av| 国产av在哪里看| 脱女人内裤的视频| 日韩欧美精品免费久久 | 国产精品永久免费网站| 制服丝袜大香蕉在线| 91久久精品电影网| 老司机福利观看| 国产一区二区亚洲精品在线观看|