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

    兩房室神經(jīng)元模型的同步簇放電的動(dòng)力學(xué)分析*

    2016-05-19 07:47:14孟盼陸啟韶趙勇董健衛(wèi)
    關(guān)鍵詞:胞體樹突房室

    孟盼陸啟韶趙勇董健衛(wèi)

    (1.廣東藥科大學(xué)基礎(chǔ)學(xué)院,廣州 510006)(2.北京航空航天大學(xué)航空科學(xué)與工程學(xué)院,北京 100191)(3.河南理工大學(xué)數(shù)學(xué)與信息科學(xué)學(xué)院,焦作 454000)

    兩房室神經(jīng)元模型的同步簇放電的動(dòng)力學(xué)分析*

    孟盼1?陸啟韶2趙勇3董健衛(wèi)1

    (1.廣東藥科大學(xué)基礎(chǔ)學(xué)院,廣州 510006)(2.北京航空航天大學(xué)航空科學(xué)與工程學(xué)院,北京 100191)(3.河南理工大學(xué)數(shù)學(xué)與信息科學(xué)學(xué)院,焦作 454000)

    研究一個(gè)具有電流反饋的兩房室錐體神經(jīng)元的簇放電動(dòng)力學(xué)問題.首先,利用數(shù)值模擬研究了該神經(jīng)元的樹突和胞體兩個(gè)房室之間的峰相位同步行為.其次,基于快慢動(dòng)力學(xué)分析方法,給出了快子系統(tǒng)的雙參數(shù)分岔曲線圖,結(jié)果表明快子系統(tǒng)產(chǎn)生了3個(gè)余維-2分岔點(diǎn),它們分別為fold-Hopf分岔點(diǎn)ZH、Cusp分岔點(diǎn)CP和Bogdanov-Takens分岔點(diǎn)BT.最后,根據(jù)這些余維-2分岔點(diǎn)確定系統(tǒng)存在著豐富的簇放電模式,如“fold/fold”型,“Hopf/Hopf”型,“Hopf/homoclinic”型以“fold/homoclinic”型簇振蕩等,并分析研究了簇放電模式的產(chǎn)生以及它們之間相互轉(zhuǎn)遷的動(dòng)力學(xué)機(jī)理.本文的研究結(jié)果進(jìn)一步表明了神經(jīng)元的幾何形態(tài)特性對于神經(jīng)元放電模式多樣性的重要影響.

    簇放電, 同步, 電流反饋, 快慢動(dòng)力學(xué)分析, 分岔

    引言

    神經(jīng)元作為神經(jīng)系統(tǒng)的基本功能單位,通過放電活動(dòng)感受外界刺激,從而實(shí)現(xiàn)神經(jīng)信息的傳遞[1].不同的放電模式是對不同外界刺激的不同反映[2-3].研究表明,簇放電(bursting)是神經(jīng)元最重要的放電模式之一,在生物信息編碼中起著關(guān)鍵作用.簇振蕩比單個(gè)峰放電(spiking)可以攜帶更多的信息并且更加可靠,因此理解簇振蕩活動(dòng)的產(chǎn)生和轉(zhuǎn)遷機(jī)制具有重要意義[4].簇振蕩表現(xiàn)為慢變量引起的靜息態(tài)(silent phase)與快變量引起的重復(fù)放電狀態(tài)(active phase)之間的相互轉(zhuǎn)遷.利用快慢動(dòng)力學(xué)分析來考察簇放電的靜息態(tài)和放電狀態(tài)之間的分岔,并依此給出簇模式的拓?fù)漕愋停抢斫獯卣袷幃a(chǎn)生機(jī)理的主要方法[5].

    為了更好的探究神經(jīng)元的各部分在信息傳遞中的作用,神經(jīng)元的形態(tài)結(jié)構(gòu)多利用房室模型來刻畫,如錐體神經(jīng)元[6],蒲肯野神經(jīng)元[7]和中間神經(jīng)元[8]等.房室模型是將空間上連續(xù)的單個(gè)神經(jīng)元分成離散的多個(gè)房室,各個(gè)房室相互連接后整體描述神經(jīng)元的時(shí)空變化.和單房室神經(jīng)元模型相比,多房室模型同樣可以產(chǎn)生豐富的放電模式,并且能夠更好的體現(xiàn)神經(jīng)元的形態(tài)結(jié)構(gòu)特征[6-8].利用房室模型來研究神經(jīng)元的電活動(dòng)行為,一直是計(jì)算神經(jīng)領(lǐng)域的熱點(diǎn).本文主要研究具有電流反饋的兩房室錐體神經(jīng)元模型,通過數(shù)值模擬和理論分析,討論了神經(jīng)元通過形態(tài)參數(shù)變化時(shí)展現(xiàn)的放電節(jié)律模式以及各種模式之間的相互轉(zhuǎn)遷機(jī)理.

    1 模型描述

    考慮錐體神經(jīng)元,它由樹突和軸突兩個(gè)房室構(gòu)成.模型描述如下[9]:

    其中門控函數(shù)為:

    (1)式和(3)式分別為胞體膜電壓VS和樹突膜電壓VD的動(dòng)力學(xué)方程,其中胞體房室的電流包含了胞體的突觸輸入電流、樹突流向胞體的內(nèi)部電流、Na+離子電流、K+離子電流和泄漏電流;樹突房室的電流包含了樹突的突觸輸入電流、樹突流向胞體的內(nèi)部電流和泄漏電流.ENa和EK分別表示Na+離子通道和K+離子通道的平衡電位,ESL和EDL分別為胞體和樹突膜上泄漏電流的平衡電位;gNa、gK、gSL和gDL分別代表各通道的最大電導(dǎo);gc表示兩房室連接電導(dǎo);C表示膜電容;IS和ID分別為胞體和樹突的突觸輸入電流;p和1-p分別表示胞體和樹突所占的面積比例.p是錐體神經(jīng)元的一個(gè)重要幾何形態(tài)參數(shù).(2)式描述了K+離子通道門開通概率w的變化規(guī)律.作為外界電流刺激的影響,(1)式中的胞體突觸輸入電流IS由一個(gè)線性反饋的動(dòng)力學(xué)方程(4)支配,這里ε是時(shí)間尺度因子,取為很小的正數(shù)以保證IS的慢變性.方程(1)-(4)稱為具有電流反饋的兩房室神經(jīng)元模型.由于ε?1,因此方程(1)-(3)可視為快子系統(tǒng),其中取IS作為慢變量,p作為控制參數(shù).在這篇文章中,其它參數(shù)取值為:ε=0.01,gNa=50,gK=15,gSL=2,gDL=2,gc=2,ENa=10,EK=-90,ESL=-70,EDL=-70,Φ=0.3,C=2,v0=-22,v1=1.2,v2=18,v3=10,v4=6.

    為了便于以后的研究,先考慮在上述參數(shù)條件下兩個(gè)房室的同步動(dòng)力學(xué)行為.引入峰相位函數(shù)φS(t)[10],定義為:

    其中ti是第i次峰放電發(fā)放的時(shí)間.因此每經(jīng)過一次峰放電,峰相位就增加2π.基于相位函數(shù)的定義,引入頻率函數(shù)ω(t)=˙φ(t)和平均峰頻率函數(shù)Ω(t)=<˙φ(t)>,這里<·>表示對時(shí)間取平均,因此有

    如果兩個(gè)房室的平均峰頻率差ΔΩ=ΩS-ΩD為0,那么兩個(gè)房室的放電頻率一致(在時(shí)間平均的意義上),從而實(shí)現(xiàn)放電相位的同步[11].

    2 神經(jīng)元形態(tài)參數(shù)引起的兩室間簇放電的同步

    通過數(shù)值計(jì)算得到ΔΩ和形態(tài)參數(shù)p的關(guān)系圖(見圖1).由圖1可以看出在0.18<p<1.0的參數(shù)范圍內(nèi)有平均峰頻率差ΔΩ=0,即兩個(gè)房室總是處于峰相位同步狀態(tài),下面的圖3~6的數(shù)值結(jié)果也證實(shí)這一點(diǎn).因此在下面的神經(jīng)元簇放電模式研究中,我們可以著重考察其中一個(gè)房室(如胞體)的放電動(dòng)力學(xué)行為,例如只依據(jù)VS的分岔圖就可以進(jìn)行快慢動(dòng)力學(xué)分析.

    圖1 平均峰頻率差ΔΩ關(guān)于p的變化圖Fig.1 Relationship of themean spike frequency differenceΔΩ and parameter p

    3 簇放電模式及轉(zhuǎn)遷

    現(xiàn)在考慮前述的兩房室錐體神經(jīng)元模型的復(fù)雜放電模式,圖2給出一些分岔結(jié)果(具體的分析過程從略).圖2(a)是全系統(tǒng)(1)-(4)關(guān)于形態(tài)參數(shù)p的分岔圖.可以看出,平衡點(diǎn)的分岔曲線是一條直線,在HB處出現(xiàn)supHopf分岔,系統(tǒng)的靜息態(tài)消失,使得穩(wěn)定平衡點(diǎn)變?yōu)椴环€(wěn)定,同時(shí)產(chǎn)生穩(wěn)定極限環(huán),系統(tǒng)開始振蕩.圖2(b)是雙參數(shù)(IS,p)-平面的分岔曲線圖.其中余維-1分岔曲線包括Hopf分岔曲線(h),平衡點(diǎn)的鞍結(jié)分岔曲線(f1,f2),同宿軌分岔曲線(hc).余維-2分岔點(diǎn)分別為fold-hopf分岔點(diǎn)ZH(Hopf分岔曲線h和鞍結(jié)分岔曲線f2的交點(diǎn))、Bogdanov-Takens分岔點(diǎn)BT(fold分岔曲線f1和同宿軌分岔曲線hc相切的點(diǎn))以及Cusp分岔點(diǎn)CP(兩條鞍結(jié)分岔曲線f1和f2處的交點(diǎn)).左下角是點(diǎn)BT附近的局部放大圖.右上角是點(diǎn)ZH附近的局部放大圖,右下角是點(diǎn)CP附近的放大圖.

    圖2 (a)全系統(tǒng)(1)-(4)關(guān)于形態(tài)參數(shù)p的分岔圖.(b)(IS,p)-參數(shù)平面上快子系統(tǒng)(1)-(3)的分岔曲線.Fig.2 (a)One-parameter bifurcation diagram of the whole system with the variation of parameter p.(b)Bifurcation curves in(IS,p)-parametric plane for the fast subsystem(1)-(3)

    為了進(jìn)一步說明問題,下面結(jié)合圖2來考慮形態(tài)參數(shù)p的變化對系統(tǒng)(1)-(4)的影響.通過取參數(shù)p不同的值來給出整個(gè)系統(tǒng)的某些典型的放電模式.隨著參數(shù)p的變化,系統(tǒng)呈現(xiàn)不同形式的簇放電模式,如圖3~6所示,其中黑線和紅線分別表示胞體和樹突的放電時(shí)間序列圖.可以看到這兩個(gè)房室的簇放電行為一直處于峰相位同步狀態(tài).接下來利用快慢動(dòng)力學(xué)分析,解釋這些簇模式之間轉(zhuǎn)遷的動(dòng)力學(xué)機(jī)理.

    首先,當(dāng)形態(tài)參數(shù)p=0.19時(shí),系統(tǒng)產(chǎn)生簇放電模式,如圖3(a)所示.快慢動(dòng)力學(xué)分析如圖3(a)所示.快子系統(tǒng)的分岔曲線是一條S形曲線,上支由穩(wěn)定焦點(diǎn)(實(shí)線)和不穩(wěn)定焦點(diǎn)(虛線)組成;中支是由鞍點(diǎn)(虛線)組成;下支是由穩(wěn)定結(jié)點(diǎn)(實(shí)線)組成.由圖2(b)可以看出,現(xiàn)在在慢變參數(shù)區(qū)域內(nèi)的快子系統(tǒng)的分岔曲線從左到右僅有平衡點(diǎn)的鞍結(jié)分岔曲線f1和f2,相應(yīng)的快子系統(tǒng)的分岔點(diǎn)分別記為F1和F2.全系統(tǒng)的軌線也疊加在快子系統(tǒng)的分岔圖上.隨著參數(shù)IS的增加,相應(yīng)于快子系統(tǒng)分岔曲線下支穩(wěn)定結(jié)點(diǎn)的簇放電的靜息態(tài),經(jīng)鞍結(jié)分岔F2轉(zhuǎn)遷到上支的放電狀態(tài).隨著IS的減小,放電狀態(tài)經(jīng)由鞍結(jié)分岔F1轉(zhuǎn)遷到下支的靜息態(tài).因此簇放電類型是“fold/fold”型簇放電[5].

    圖3 當(dāng)p=0.19時(shí),“fold/fold”型簇放電.(a)兩個(gè)房室的膜電位VS(黑線),VD(紅線)的時(shí)間歷程圖.(b)快慢動(dòng)力學(xué)分岔分析.Fig.3 Bursting of“fold/fold”type when p=0.19.(a)Time history ofmembrane potentials VSin the black line and VDin the red one.(b)Fast-slow dynamics of the fast subsystem.

    圖4 當(dāng)p=0.35時(shí),“Hopf/Hopf”型簇放電.(a)兩個(gè)房室的膜電位VS(黑線),VD(紅線)的時(shí)間歷程圖.(b)快慢動(dòng)力學(xué)分岔分析.Fig.4 Bursting of“fold/fold”type when p=0.35.(a)Time history ofmembrane potentials VSin the black line and VDin the red one.(b)Fast-slow dynamics of the fast subsystem

    其次,當(dāng)p=0.35時(shí),系統(tǒng)可產(chǎn)生如圖4(a)的簇放電模式.運(yùn)用圖4(b)的快慢動(dòng)力學(xué)分析,來研究此放電的動(dòng)力學(xué)性質(zhì)及其類型.通過雙參數(shù)(IS,p)-平面分岔圖2(b)可以看出,在兩個(gè)鞍結(jié)分岔點(diǎn)F1和F2之間,產(chǎn)生了兩個(gè)supHopf分岔點(diǎn),在快子系統(tǒng)的分岔圖上分別記為HB1和HB2.這樣,S形分岔曲線的上支的穩(wěn)定焦點(diǎn)(實(shí)線)經(jīng)由點(diǎn)HB1處的Hopf分岔而失穩(wěn)(虛線),同時(shí)產(chǎn)生穩(wěn)定極限環(huán).隨著參數(shù)IS的增加,該穩(wěn)定極限環(huán)經(jīng)由點(diǎn)HB2處的Hopf分岔消失而轉(zhuǎn)變?yōu)榉€(wěn)定焦點(diǎn).現(xiàn)在結(jié)合全系統(tǒng)的軌線來討論簇放電的產(chǎn)生過程.隨著慢變量IS的增加,由于supHopf分岔點(diǎn)HB2處的慢變效應(yīng)[11],下支的靜息狀態(tài)轉(zhuǎn)遷到上支的放電狀態(tài);隨著慢變量IS的減小,放電狀態(tài)收斂于經(jīng)由supHopf分岔點(diǎn)HB1產(chǎn)生的穩(wěn)定焦點(diǎn),再經(jīng)由鞍結(jié)分岔F1轉(zhuǎn)遷到下支的靜息狀態(tài).這種簇振蕩模式稱為“Hopf/Hopf”型簇放電[5].綜上所述,當(dāng)參數(shù)p的取值小于余維-2分岔點(diǎn)BT的縱坐標(biāo)pBT值時(shí),系統(tǒng)會(huì)產(chǎn)生“fold/fold”型或“Hopf/Hopf”型簇放電.

    接下來,當(dāng)形態(tài)參數(shù)p的取值位于兩個(gè)余維-2分岔點(diǎn)Bogdanov-Takens分岔點(diǎn)BT和fold-hopf分岔點(diǎn)ZH的縱坐標(biāo)值之間,例如p=0.42時(shí),兩個(gè)房室的簇放電模式如圖5(a)所示,其相應(yīng)的分岔分析見圖5(b).和圖4(b)相比,現(xiàn)在快子系統(tǒng)僅有一個(gè)SupHopf分岔點(diǎn)HB.神經(jīng)元的靜息態(tài)經(jīng)由SupHopf分岔而轉(zhuǎn)遷到放電態(tài),放電態(tài)又經(jīng)極限環(huán)的同宿軌分岔HC轉(zhuǎn)變?yōu)殪o息態(tài),從而形成了一個(gè)“Hopf/Homoclinic”型簇放電[5].

    圖5 當(dāng)p=0.42時(shí),“Hopf/Homoclinic”型簇放電.(a)兩個(gè)房室的膜電位VS(黑線)、VD(紅線)的時(shí)間歷程圖.(b)快慢動(dòng)力學(xué)分岔分析.Fig.5 Bursting of“Hopf/Homoclinic”type when p=0.42.(a)Time history ofmembrane potentials VSin the black line and VDin the red one,(b)Fast-slow dynamics of the fast subsystem.

    圖6 當(dāng)p=0.9時(shí),“fold/homoclinic”型簇放電.(a)兩個(gè)房室的膜電位VS(黑線)、VD(紅線)的時(shí)間歷程圖.(b)快慢動(dòng)力學(xué)分岔分析.Fig.6 Bursting of“fold/Homoclinic”type when p=0.9,(a)Time history ofmembrane potentials VSin the black line and VDin the red one.(b)Fast-slow dynamics of the fast subsystem.

    最后,當(dāng)形態(tài)參數(shù)p的取值大于余維-2分岔點(diǎn)fold-hopf分岔點(diǎn)ZH的縱坐標(biāo)pZH值,例如p=0.9時(shí),系統(tǒng)可產(chǎn)生如圖6(a)的簇放電模式,相應(yīng)的快慢動(dòng)力學(xué)分析如圖6(b)所示.隨著慢變量的增加,快子系統(tǒng)依次經(jīng)過鞍結(jié)分岔F1、同宿軌分岔HC、鞍結(jié)分岔F2、supHopf分岔HB.值得注意的是,現(xiàn)在鞍結(jié)分岔點(diǎn)F2和supHopf分岔點(diǎn)HB的先后順序和圖5(b)有所不同.因此,系統(tǒng)軌線的靜息態(tài)由平衡點(diǎn)的鞍結(jié)分岔(F2)轉(zhuǎn)變?yōu)榉烹姞顟B(tài),然后放電狀態(tài)沿著由HB產(chǎn)生的穩(wěn)定極限環(huán)運(yùn)動(dòng),最后消失于極限環(huán)的同宿軌分岔(HC).這種簇振蕩模式稱為“fold/homoclinic”型簇放電[5].

    4 結(jié)論

    通過數(shù)值模擬和快慢動(dòng)力學(xué)分析,研究了具有電流反饋的兩房室錐體神經(jīng)元的簇放電模式.結(jié)果表明,胞體和樹突這兩個(gè)房室在一定條件下可以處于峰相位同步狀態(tài):神經(jīng)元可以表現(xiàn)出豐富的放電模式,如“fold/fold”型、“Hopf/Hopf”型、“Hopf/Homoclinic”型和“fold/homoclinic”型簇放電模式;特別地,神經(jīng)元的幾何形態(tài)參數(shù)對神經(jīng)元放電節(jié)律的產(chǎn)生和轉(zhuǎn)遷有著重要的影響.還注意到,如果不考慮電流反饋的話,此兩房室神經(jīng)元模型并不存在這些豐富的簇放電模式,因此本文的研究對進(jìn)一步揭示外界激勵(lì)對神經(jīng)元放電活動(dòng)的重要作用也有一定意義.

    1 汪云九.神經(jīng)信息學(xué)-神經(jīng)系統(tǒng)的理論和模型.北京,高等教育出版社,2006(Wang Y J.Neuro-informationtics:theories and models of neural systems.Beijing:Higher Education Press,2006(in Chinese))

    2 Fan Y S,Chay T R.Generation of periodic and chaotic bursting in an excitable cellmodel.Biological Gybernetics,1994,71(5):417~431

    3 Gu H G,Pan B B.Identification of neural firing patterns,frequency and temporal coding mechanisms in individual aortic baroreceptors.Frontiers in Computational Neuroscience,2015,9:108

    4 Lisman JE.Bursts as a unit of neural information:making unreliable synapses reliable.Trends in Neuroscience,1997,20(1):38~43

    5 Rinzel J.Bursting oscillations in an excitable membranemodel.In:Ordinary and Partial Differential Equations,Lecture Notes in Mathematics.edited by Sleeman B D,Jarvis R J,New York:Springer,1995:304~316

    6 馬丹,劉深泉,汪雷.兩房室神經(jīng)元模型的分岔現(xiàn)象.北京生物醫(yī)學(xué)工程,2011,30(6):567~573(Ma d,Liu S Q,Wang L.Bifurcation phenomenon in a two-compartment neuron model.Beijing Biomedical Engineering,2011,30(6):567~573(in Chinese))

    7 張煦晨,劉深泉,任會(huì)霞.浦肯野神經(jīng)元回路中的同步性研究.動(dòng)力學(xué)與控制學(xué)報(bào),2015,13(5):377~381(Zhang X C,Liu SQ,Ren H X.The research on synchronization in purinje neural circuit.Journal of Dynamics and Control,2015,13(5):377~381(in Chinese))

    8 江小芳,劉深泉,張煦晨.中等多棘神經(jīng)元的多房室模型分析.生物物理學(xué)報(bào),2013,29(7):527~542(Jiang X F,Liu SQ,Zhang X C.Multi-compartmentmodelanalysis of the medium spiny neuron.Biphysics Reports,2013,29(7):527~542(in Chinese))

    9 伊國勝,王江,魏熙樂,鄧斌.電場調(diào)制神經(jīng)元放電動(dòng)力學(xué)的生物物理機(jī)制.中國科學(xué):物理學(xué)·力學(xué)·天文學(xué),2013,42(2):210~216.(Yin G S,Wang J,Wei X L,Deng B.The mechanism of electrical field modulates the neuron firing pattern dynamics.Science China Physics,Mechanics&Astronomy,2013,42(2):210~216(in Chinese))

    10王海俠,陸啟韶.對稱電耦合連接方式下的多尺度同步行為和轉(zhuǎn)遷過程.動(dòng)力學(xué)與控制學(xué)報(bào),2011,9(3):214~218(Wang H X,Lu Q S.Multi-time-scale synchronization and transitions with symmetrical and electrical coupling schemes.Journal of Dynamics and Control,2011,9(3):214~218(in Chinese))

    11 Baer S M,Erneux T,Rinzel J.Slow passage through a Hopf bifurcation:Delay,memory effects,and resonance. SIAM Journal on Applied Mathematics,1989,49:55~71

    DYNAM ICS ANALYSISOF SYNCHRONOUS BURSTING IN TWO-COMPARTMENT NEURON MODEL*

    Meng Pan1?Lu Qishao2Zhao Yong3Dong Jianwei1
    (1.School of Basic Courses,Guangdong Pharmaceutical University,Guangzhou 510006,China)(2.School ofAeronautical Science and Engineering,Beihang University,Beijing 100083,China)(3.School of Mathematics and Information Science,Henan Polytechnic University,Jiaozuo 454000,China)

    A two-compartment neuronmodelwith outside current feedback is considered in order to study the dynamics of bursting patterns.Firstly,on the base of numerical simulation,the behavior of phase synchronization for the spikes in the soma and dendrite compartments is investigated.Secondly,by using fast-slow dynamics analysis,two-parameter bifurcation analysis of the fast subsystem is performed in the two-parameter plane.It is shown that three condimension-2 bifurcation points are generated in the fast subsystem,i.e.,fold-Hopf bifurcation ZH,cusp bifurcation CP and Bogdanov-Takens bifurcation BT.Finally,these condimension-2 bifurcation points determine that different bursting oscillations occurs in the system,such as“fold/fold”type,“Hopf/Hopf”type,“Hopf/homoclinic”type and“fold/homoclinic”type.Furthermore,the generation and transition mechanisms of different bursting patterns are discussed.The results offer significant data for understanding the important role of the geometric parameter in the variety of neuronal firing activities

    bursting, synchronization, current-feedback, fast-slow analysis, bifurcation

    10.6052/1672-6553-2016-023

    2016-01-31收到第1稿,2016-03-26收到修改稿.

    *國家自然科學(xué)基金資助項(xiàng)目(11402057,11502073,11272034,11371046)

    ?通訊作者E-mail:mengpan200e@163.com

    Received 31 January 2016,revised 26 March 2016.

    *The project supported by the National Natural Science Foundation of China(11402057,11502073,11272024 and 11371046)

    ?Corresponding author E-mail:mengpan200e@163.com

    猜你喜歡
    胞體樹突房室
    防沖液壓支架點(diǎn)陣吸能結(jié)構(gòu)設(shè)計(jì)及其性能分析
    科學(xué)家揭示大腦連接的真實(shí)結(jié)構(gòu) 解決了樹突棘保存難題
    海外星云(2021年6期)2021-10-14 07:20:40
    房室交接區(qū)期前收縮致復(fù)雜心電圖表現(xiàn)1 例
    基于正多胞體空間擴(kuò)展濾波的時(shí)變參數(shù)系統(tǒng)辨識方法
    復(fù)雜開曲面的魯棒布爾運(yùn)算
    房室阻滯表現(xiàn)多變的臨床心電圖分析
    siRNA干預(yù)樹突狀細(xì)胞CD40表達(dá)對大鼠炎癥性腸病的治療作用
    經(jīng)食管心臟電生理檢測房室交界區(qū)前傳功能
    樹突狀細(xì)胞疫苗抗腫瘤免疫研究進(jìn)展
    徽章樣真皮樹突細(xì)胞錯(cuò)構(gòu)瘤三例
    成人鲁丝片一二三区免费| 亚洲最大成人手机在线| 亚洲成人av在线免费| 女人被狂操c到高潮| 久久影院123| 边亲边吃奶的免费视频| 天堂网av新在线| 美女内射精品一级片tv| 国产极品天堂在线| 97在线人人人人妻| 国产综合精华液| 亚洲成人中文字幕在线播放| 十八禁网站网址无遮挡 | 黄色视频在线播放观看不卡| 国产成人精品婷婷| 日韩一本色道免费dvd| 成人欧美大片| 99久国产av精品国产电影| 日韩欧美精品v在线| 99re6热这里在线精品视频| 国产精品久久久久久精品电影小说 | 97热精品久久久久久| 中文欧美无线码| 青春草视频在线免费观看| 欧美日韩视频高清一区二区三区二| 午夜福利在线观看免费完整高清在| 精品亚洲乱码少妇综合久久| 一级毛片久久久久久久久女| 91久久精品电影网| 丝袜喷水一区| 欧美性感艳星| 欧美精品一区二区大全| 夜夜爽夜夜爽视频| 一本一本综合久久| 日韩欧美精品v在线| 国产成人午夜福利电影在线观看| 国模一区二区三区四区视频| 草草在线视频免费看| 色5月婷婷丁香| 人妻系列 视频| 久久国产乱子免费精品| 熟女av电影| 性色av一级| 欧美另类一区| 亚洲成人av在线免费| 欧美成人午夜免费资源| 一级毛片我不卡| 最新中文字幕久久久久| 男女无遮挡免费网站观看| 毛片一级片免费看久久久久| 最新中文字幕久久久久| 少妇高潮的动态图| 晚上一个人看的免费电影| 精品少妇久久久久久888优播| 亚洲av欧美aⅴ国产| 香蕉精品网在线| 亚洲欧洲日产国产| 99视频精品全部免费 在线| 日本三级黄在线观看| 69人妻影院| 国产极品天堂在线| 国产精品人妻久久久久久| 亚洲av免费高清在线观看| 免费不卡的大黄色大毛片视频在线观看| 亚洲第一区二区三区不卡| 午夜爱爱视频在线播放| 久久99热这里只有精品18| 大香蕉久久网| 国产精品久久久久久久久免| 综合色av麻豆| 1000部很黄的大片| 国产成人精品一,二区| av国产精品久久久久影院| 男女边摸边吃奶| 日本一本二区三区精品| 久久久久网色| 成人毛片a级毛片在线播放| 欧美丝袜亚洲另类| 婷婷色综合大香蕉| 午夜精品一区二区三区免费看| 精品久久久久久久末码| 亚洲精品视频女| 特级一级黄色大片| 精品人妻熟女av久视频| 日韩三级伦理在线观看| 三级国产精品欧美在线观看| 一级爰片在线观看| 亚洲国产精品成人综合色| 久久热精品热| 欧美潮喷喷水| 男人爽女人下面视频在线观看| 国产精品一及| 国产高清不卡午夜福利| 欧美成人午夜免费资源| 大片免费播放器 马上看| 欧美区成人在线视频| 大又大粗又爽又黄少妇毛片口| 乱系列少妇在线播放| 在线观看免费高清a一片| 久久精品久久久久久噜噜老黄| 超碰97精品在线观看| 久久人人爽人人爽人人片va| 色综合色国产| 熟妇人妻不卡中文字幕| 欧美成人a在线观看| 91午夜精品亚洲一区二区三区| 天堂网av新在线| 精品久久久久久久末码| 久久久久精品久久久久真实原创| 午夜视频国产福利| 久久国内精品自在自线图片| 久久99蜜桃精品久久| 成年av动漫网址| 午夜免费鲁丝| 菩萨蛮人人尽说江南好唐韦庄| 精品久久久噜噜| 91精品伊人久久大香线蕉| 国产一区二区在线观看日韩| 午夜福利在线观看免费完整高清在| a级毛色黄片| 神马国产精品三级电影在线观看| 美女国产视频在线观看| 最近2019中文字幕mv第一页| 18禁裸乳无遮挡免费网站照片| 日韩中字成人| 一个人看的www免费观看视频| 成人午夜精彩视频在线观看| 99久久人妻综合| 嫩草影院入口| 夫妻午夜视频| 3wmmmm亚洲av在线观看| 国内精品美女久久久久久| 久久久久国产网址| 日韩在线高清观看一区二区三区| 国国产精品蜜臀av免费| 久久久久精品久久久久真实原创| 精品国产乱码久久久久久小说| 纵有疾风起免费观看全集完整版| 亚洲精品国产av蜜桃| 午夜福利网站1000一区二区三区| 亚洲丝袜综合中文字幕| 久久99热6这里只有精品| 王馨瑶露胸无遮挡在线观看| 久久精品国产自在天天线| 夫妻性生交免费视频一级片| 亚洲三级黄色毛片| 亚洲精品第二区| 国产男人的电影天堂91| 亚洲精品aⅴ在线观看| 尾随美女入室| 国产日韩欧美亚洲二区| 国产 精品1| 亚洲色图综合在线观看| 国产黄a三级三级三级人| 欧美日韩在线观看h| 制服丝袜香蕉在线| 成人综合一区亚洲| 久久久久性生活片| 天天躁日日操中文字幕| 综合色av麻豆| 亚洲天堂av无毛| 欧美激情久久久久久爽电影| 久久6这里有精品| 国产有黄有色有爽视频| 免费看日本二区| 女人十人毛片免费观看3o分钟| 尾随美女入室| 真实男女啪啪啪动态图| 久久精品久久精品一区二区三区| 少妇人妻 视频| av在线蜜桃| 国产色婷婷99| 欧美xxxx黑人xx丫x性爽| 男的添女的下面高潮视频| 国产一区二区三区综合在线观看 | 色哟哟·www| 久久久久久久久久久免费av| 国产色爽女视频免费观看| 精品人妻熟女av久视频| 中国国产av一级| 大码成人一级视频| 日本与韩国留学比较| 一级毛片黄色毛片免费观看视频| 如何舔出高潮| 丰满少妇做爰视频| 18禁在线播放成人免费| 国产探花极品一区二区| 如何舔出高潮| av女优亚洲男人天堂| 女人十人毛片免费观看3o分钟| 天堂中文最新版在线下载 | 男的添女的下面高潮视频| 熟妇人妻不卡中文字幕| 国产免费一区二区三区四区乱码| 18禁在线播放成人免费| 日本猛色少妇xxxxx猛交久久| 亚洲怡红院男人天堂| 亚洲丝袜综合中文字幕| 亚洲精品成人av观看孕妇| 婷婷色综合大香蕉| 大香蕉久久网| 久久久精品免费免费高清| 亚洲不卡免费看| 亚洲国产精品成人久久小说| 久久综合国产亚洲精品| 国产精品秋霞免费鲁丝片| 99久久九九国产精品国产免费| 在现免费观看毛片| 天天躁日日操中文字幕| 美女主播在线视频| 午夜福利高清视频| 三级经典国产精品| 亚洲,欧美,日韩| 国产精品蜜桃在线观看| 插阴视频在线观看视频| 有码 亚洲区| 视频区图区小说| 肉色欧美久久久久久久蜜桃 | 亚洲av男天堂| 亚洲一区二区三区欧美精品 | 国产女主播在线喷水免费视频网站| 国产毛片a区久久久久| 秋霞伦理黄片| 高清在线视频一区二区三区| 99热这里只有精品一区| 水蜜桃什么品种好| 国产高清有码在线观看视频| 欧美国产精品一级二级三级 | 欧美老熟妇乱子伦牲交| 精品亚洲乱码少妇综合久久| 水蜜桃什么品种好| 一区二区三区乱码不卡18| 久久精品久久精品一区二区三区| 国产乱人视频| 黄色一级大片看看| 国产成人一区二区在线| 欧美高清性xxxxhd video| 男女啪啪激烈高潮av片| www.av在线官网国产| 尤物成人国产欧美一区二区三区| 亚洲欧美清纯卡通| 亚洲av男天堂| 国产精品99久久久久久久久| 深夜a级毛片| 1000部很黄的大片| 日韩成人伦理影院| 啦啦啦啦在线视频资源| 搡老乐熟女国产| 国产男人的电影天堂91| www.av在线官网国产| 日韩视频在线欧美| 欧美xxxx性猛交bbbb| 男男h啪啪无遮挡| 国产一区二区亚洲精品在线观看| 美女主播在线视频| 中文资源天堂在线| 国产在线男女| 蜜桃久久精品国产亚洲av| 国产av国产精品国产| 亚洲国产欧美人成| 亚州av有码| 波野结衣二区三区在线| 精品少妇黑人巨大在线播放| 国产伦理片在线播放av一区| 极品少妇高潮喷水抽搐| 高清毛片免费看| 丰满乱子伦码专区| 久久久午夜欧美精品| 国产老妇女一区| 亚洲成色77777| 搞女人的毛片| 18+在线观看网站| 又黄又爽又刺激的免费视频.| 久久精品久久久久久久性| 欧美日韩国产mv在线观看视频 | 久久精品国产鲁丝片午夜精品| 国模一区二区三区四区视频| 七月丁香在线播放| 欧美日本视频| 欧美变态另类bdsm刘玥| 建设人人有责人人尽责人人享有的 | 女人久久www免费人成看片| 中文字幕久久专区| 日本一本二区三区精品| 国产日韩欧美亚洲二区| 午夜老司机福利剧场| 久久影院123| 黑人高潮一二区| 精品久久国产蜜桃| 日日撸夜夜添| 精品一区二区免费观看| 久久女婷五月综合色啪小说 | 女人久久www免费人成看片| 赤兔流量卡办理| 91精品伊人久久大香线蕉| 国产av不卡久久| av在线蜜桃| 亚洲成人精品中文字幕电影| 交换朋友夫妻互换小说| 26uuu在线亚洲综合色| 女人被狂操c到高潮| 校园人妻丝袜中文字幕| 国产高潮美女av| 国产午夜福利久久久久久| 免费观看在线日韩| 亚洲欧美日韩东京热| 国产毛片在线视频| 如何舔出高潮| 欧美成人一区二区免费高清观看| videossex国产| 亚洲天堂av无毛| 汤姆久久久久久久影院中文字幕| 精品酒店卫生间| 亚洲最大成人中文| 熟女av电影| 久久久久网色| 日本与韩国留学比较| 久久99蜜桃精品久久| 精品久久久久久电影网| 自拍偷自拍亚洲精品老妇| freevideosex欧美| 亚洲精品日韩在线中文字幕| 午夜日本视频在线| www.色视频.com| 午夜免费观看性视频| 91久久精品国产一区二区成人| 国产精品三级大全| 欧美精品国产亚洲| 国产极品天堂在线| 国产黄色免费在线视频| 精品人妻视频免费看| 国产成人精品福利久久| 人妻一区二区av| 美女视频免费永久观看网站| 久久99精品国语久久久| 人妻制服诱惑在线中文字幕| 99精国产麻豆久久婷婷| 欧美日本视频| 久久精品国产亚洲av涩爱| 黄色一级大片看看| 国产精品一区二区三区四区免费观看| 国产白丝娇喘喷水9色精品| 黄色配什么色好看| 国模一区二区三区四区视频| 亚洲自拍偷在线| 亚洲精品日韩av片在线观看| 黄色一级大片看看| 人人妻人人澡人人爽人人夜夜| 国产精品人妻久久久久久| 成年av动漫网址| 久久精品久久久久久久性| 国产中年淑女户外野战色| 日本欧美国产在线视频| 六月丁香七月| 午夜老司机福利剧场| 乱系列少妇在线播放| 男人狂女人下面高潮的视频| 久久久久国产精品人妻一区二区| 亚洲精品影视一区二区三区av| 精品久久久久久久末码| 看十八女毛片水多多多| 少妇高潮的动态图| 国产白丝娇喘喷水9色精品| 三级经典国产精品| 特级一级黄色大片| 欧美+日韩+精品| 国产免费又黄又爽又色| 国产午夜精品一二区理论片| 国产探花极品一区二区| 欧美最新免费一区二区三区| 内射极品少妇av片p| 日本熟妇午夜| 免费观看在线日韩| 国产黄频视频在线观看| 五月伊人婷婷丁香| 在线观看三级黄色| 激情五月婷婷亚洲| 在线观看免费高清a一片| av卡一久久| 亚洲欧美清纯卡通| 丝袜美腿在线中文| 一级毛片久久久久久久久女| 噜噜噜噜噜久久久久久91| 2021天堂中文幕一二区在线观| 欧美精品人与动牲交sv欧美| 国产精品秋霞免费鲁丝片| 欧美少妇被猛烈插入视频| 美女内射精品一级片tv| 18禁裸乳无遮挡动漫免费视频 | 男人狂女人下面高潮的视频| 最近最新中文字幕免费大全7| 在线观看三级黄色| 欧美国产精品一级二级三级 | 人人妻人人爽人人添夜夜欢视频 | 老女人水多毛片| 女人久久www免费人成看片| 欧美老熟妇乱子伦牲交| 免费看光身美女| 国产精品一二三区在线看| 女人被狂操c到高潮| 久久久欧美国产精品| 18禁裸乳无遮挡免费网站照片| 哪个播放器可以免费观看大片| 91久久精品电影网| av免费观看日本| 大片电影免费在线观看免费| 国产精品人妻久久久影院| 国产黄色免费在线视频| 黄片无遮挡物在线观看| 国产午夜精品久久久久久一区二区三区| 天堂俺去俺来也www色官网| av在线老鸭窝| 天天躁夜夜躁狠狠久久av| 美女xxoo啪啪120秒动态图| 亚洲va在线va天堂va国产| 男女无遮挡免费网站观看| 九九在线视频观看精品| 自拍偷自拍亚洲精品老妇| 国产黄片视频在线免费观看| 51国产日韩欧美| 免费看日本二区| 一级爰片在线观看| 国产精品人妻久久久久久| 日本免费在线观看一区| 国产黄色免费在线视频| 丰满人妻一区二区三区视频av| 欧美性猛交╳xxx乱大交人| 国产精品熟女久久久久浪| 精品人妻一区二区三区麻豆| 日韩国内少妇激情av| freevideosex欧美| 啦啦啦在线观看免费高清www| 欧美激情在线99| 国产黄色视频一区二区在线观看| 久久久久久久久久久丰满| 亚洲无线观看免费| 又大又黄又爽视频免费| 麻豆国产97在线/欧美| 天天躁日日操中文字幕| 国产精品久久久久久精品古装| 国产午夜精品一二区理论片| 色网站视频免费| 搡女人真爽免费视频火全软件| 国产精品人妻久久久影院| 久久韩国三级中文字幕| 中文乱码字字幕精品一区二区三区| 亚洲不卡免费看| 一个人看视频在线观看www免费| 男女边摸边吃奶| 一个人看视频在线观看www免费| 免费观看在线日韩| 精品一区二区免费观看| 亚洲av一区综合| videossex国产| 在线亚洲精品国产二区图片欧美 | 午夜精品一区二区三区免费看| 欧美+日韩+精品| 亚洲精品视频女| 精品午夜福利在线看| 国产成人a∨麻豆精品| 亚洲成人一二三区av| 极品教师在线视频| 人人妻人人爽人人添夜夜欢视频 | 高清欧美精品videossex| 内地一区二区视频在线| 久久久久久久亚洲中文字幕| 亚洲三级黄色毛片| 国产老妇女一区| 亚洲欧美成人综合另类久久久| 久久久精品欧美日韩精品| 少妇人妻精品综合一区二区| 欧美 日韩 精品 国产| 成人午夜精彩视频在线观看| 在线a可以看的网站| 亚洲欧美精品自产自拍| 国产真实伦视频高清在线观看| a级毛片免费高清观看在线播放| 在线免费十八禁| 好男人视频免费观看在线| 精品久久久久久久末码| 一级黄片播放器| 日韩一区二区视频免费看| 夜夜看夜夜爽夜夜摸| 一级毛片 在线播放| 婷婷色综合www| 97在线人人人人妻| 国产熟女欧美一区二区| 精品国产三级普通话版| 18禁在线播放成人免费| 青春草国产在线视频| 尾随美女入室| 麻豆久久精品国产亚洲av| 亚洲一级一片aⅴ在线观看| 一个人观看的视频www高清免费观看| 亚洲国产成人一精品久久久| 91狼人影院| 国产黄色免费在线视频| 男女国产视频网站| 黄色视频在线播放观看不卡| 成年免费大片在线观看| 欧美高清成人免费视频www| 国产免费一级a男人的天堂| 精品国产三级普通话版| 欧美成人a在线观看| 亚洲国产成人一精品久久久| av专区在线播放| 永久免费av网站大全| 日韩不卡一区二区三区视频在线| 亚洲av一区综合| 成人鲁丝片一二三区免费| 色5月婷婷丁香| 亚洲精品久久久久久婷婷小说| 大片电影免费在线观看免费| 成人鲁丝片一二三区免费| 搞女人的毛片| 亚洲欧美成人综合另类久久久| 国语对白做爰xxxⅹ性视频网站| 国产精品国产av在线观看| 黑人高潮一二区| 欧美高清性xxxxhd video| 中文欧美无线码| 亚洲欧美中文字幕日韩二区| 国产探花极品一区二区| 久久久国产一区二区| 久久久a久久爽久久v久久| 91aial.com中文字幕在线观看| 最近2019中文字幕mv第一页| 国产熟女欧美一区二区| av卡一久久| 国产精品99久久久久久久久| 一级毛片电影观看| 国产一区有黄有色的免费视频| 中文欧美无线码| 国产美女午夜福利| 亚洲av日韩在线播放| 天天躁日日操中文字幕| 99re6热这里在线精品视频| 听说在线观看完整版免费高清| 午夜免费观看性视频| 国产精品国产三级国产av玫瑰| 国产精品蜜桃在线观看| 欧美日韩在线观看h| 国产精品一区www在线观看| h日本视频在线播放| 日韩成人av中文字幕在线观看| 97超碰精品成人国产| 国产一区二区亚洲精品在线观看| 国产在线一区二区三区精| 成人黄色视频免费在线看| 日韩一区二区三区影片| 国产有黄有色有爽视频| 欧美激情在线99| 免费看不卡的av| 国产精品三级大全| 成人综合一区亚洲| 亚洲一级一片aⅴ在线观看| 直男gayav资源| 亚洲美女视频黄频| 国产爽快片一区二区三区| 精华霜和精华液先用哪个| 麻豆精品久久久久久蜜桃| 亚洲av欧美aⅴ国产| 2021少妇久久久久久久久久久| 你懂的网址亚洲精品在线观看| 国模一区二区三区四区视频| 亚洲在久久综合| 日本av手机在线免费观看| 男女无遮挡免费网站观看| 亚洲av福利一区| 日本一本二区三区精品| 亚洲av免费在线观看| 涩涩av久久男人的天堂| 久久99热这里只有精品18| 丝瓜视频免费看黄片| 国产精品福利在线免费观看| 久久精品夜色国产| 亚洲精品亚洲一区二区| 亚洲精品第二区| 国产极品天堂在线| 国产中年淑女户外野战色| 免费高清在线观看视频在线观看| 麻豆精品久久久久久蜜桃| 免费看a级黄色片| 超碰97精品在线观看| 亚洲天堂av无毛| 中文字幕免费在线视频6| 欧美日韩精品成人综合77777| 免费看日本二区| 久久精品国产自在天天线| 九九爱精品视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 别揉我奶头 嗯啊视频| 亚洲美女视频黄频| 亚洲精品一区蜜桃| 色吧在线观看| 日本一本二区三区精品| 国产亚洲91精品色在线| 一级毛片aaaaaa免费看小| 午夜激情久久久久久久| 联通29元200g的流量卡| 一级毛片电影观看| 日韩成人伦理影院| 深爱激情五月婷婷| 伊人久久精品亚洲午夜| 欧美zozozo另类| 亚洲av二区三区四区| 精品久久久久久久末码| 免费观看在线日韩| 简卡轻食公司| 汤姆久久久久久久影院中文字幕| 99热这里只有是精品在线观看| 日本与韩国留学比较| 狠狠精品人妻久久久久久综合| 国产精品av视频在线免费观看| 九九爱精品视频在线观看|