• <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片在线观看秒播厂| 考比视频在线观看| 日本-黄色视频高清免费观看| 婷婷色综合大香蕉| 国产极品粉嫩免费观看在线| 成人手机av| 男的添女的下面高潮视频| 国产有黄有色有爽视频| 亚洲在久久综合| 国产伦理片在线播放av一区| 国产精品香港三级国产av潘金莲 | 午夜日本视频在线| 国产女主播在线喷水免费视频网站| www.av在线官网国产| 国产精品不卡视频一区二区| 1024香蕉在线观看| 国产不卡av网站在线观看| h视频一区二区三区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 日本猛色少妇xxxxx猛交久久| videossex国产| 婷婷成人精品国产| 国产探花极品一区二区| 亚洲熟女精品中文字幕| av天堂久久9| 五月天丁香电影| av在线app专区| 黄色一级大片看看| 视频区图区小说| 黑人欧美特级aaaaaa片| 97人妻天天添夜夜摸| 国产精品麻豆人妻色哟哟久久| 性高湖久久久久久久久免费观看| 日本av免费视频播放| 久久久精品94久久精品| 国产成人aa在线观看| 国产免费又黄又爽又色| 国产女主播在线喷水免费视频网站| 18禁裸乳无遮挡动漫免费视频| av卡一久久| 欧美亚洲日本最大视频资源| 日本黄色日本黄色录像| 国产欧美日韩综合在线一区二区| 少妇猛男粗大的猛烈进出视频| 成人二区视频| 久久av网站| 亚洲精品视频女| 国产一区二区在线观看av| 中文精品一卡2卡3卡4更新| 大片电影免费在线观看免费| 黑人欧美特级aaaaaa片| 国产高清国产精品国产三级| 久久青草综合色| 在现免费观看毛片| 高清在线视频一区二区三区| 99国产精品免费福利视频| 亚洲精品日韩在线中文字幕| 国产精品人妻久久久影院| 国产在视频线精品| 久久久久久人人人人人| 亚洲男人天堂网一区| 男的添女的下面高潮视频| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲国产最新在线播放| 国产爽快片一区二区三区| 午夜免费男女啪啪视频观看| av网站免费在线观看视频| av又黄又爽大尺度在线免费看| 丰满少妇做爰视频| 伊人久久大香线蕉亚洲五| 国产一区二区三区综合在线观看| 成人手机av| 两个人看的免费小视频| 亚洲欧美日韩另类电影网站| 涩涩av久久男人的天堂| 亚洲精品久久成人aⅴ小说| www.av在线官网国产| 午夜福利影视在线免费观看| 中文字幕制服av| 80岁老熟妇乱子伦牲交| 捣出白浆h1v1| 黄色配什么色好看| 色网站视频免费| 纵有疾风起免费观看全集完整版| 欧美精品一区二区大全| 日韩人妻精品一区2区三区| 亚洲美女黄色视频免费看| 亚洲伊人久久精品综合| 精品久久蜜臀av无| 黑人猛操日本美女一级片| 国产人伦9x9x在线观看 | 国产 精品1| 亚洲精品日本国产第一区| 国产乱来视频区| 一区二区三区激情视频| 久久精品国产鲁丝片午夜精品| 成年av动漫网址| 久久午夜福利片| 熟妇人妻不卡中文字幕| 免费黄色在线免费观看| 日韩 亚洲 欧美在线| 久久久久精品性色| 精品久久久精品久久久| 国产一区二区 视频在线| 精品一区二区免费观看| 国产精品一国产av| 菩萨蛮人人尽说江南好唐韦庄| 亚洲综合色网址| 国产精品女同一区二区软件| 久久久久精品性色| 国产精品麻豆人妻色哟哟久久| 久久亚洲国产成人精品v| 少妇被粗大的猛进出69影院| 熟女av电影| 久久久久久久久久人人人人人人| 天天操日日干夜夜撸| 黄片无遮挡物在线观看| 国产麻豆69| 久久久国产一区二区| 搡老乐熟女国产| 色播在线永久视频| 亚洲欧美日韩另类电影网站| 18禁观看日本| 美女中出高潮动态图| 伦理电影大哥的女人| 母亲3免费完整高清在线观看 | 欧美成人午夜精品| 免费高清在线观看日韩| 多毛熟女@视频| 一级a爱视频在线免费观看| 久久久欧美国产精品| 成人手机av| 黄色毛片三级朝国网站| 精品一区二区三卡| 超色免费av| 亚洲 欧美一区二区三区| 久久久久久久大尺度免费视频| 日韩熟女老妇一区二区性免费视频| 久久这里只有精品19| 人妻人人澡人人爽人人| kizo精华| 精品国产一区二区久久| 国产精品久久久久久精品电影小说| 精品人妻在线不人妻| freevideosex欧美| 青青草视频在线视频观看| 狠狠精品人妻久久久久久综合| 欧美最新免费一区二区三区| 男女边吃奶边做爰视频| 纯流量卡能插随身wifi吗| 国产免费现黄频在线看| √禁漫天堂资源中文www| 麻豆乱淫一区二区| 久久国产精品大桥未久av| 久久久久视频综合| 免费看av在线观看网站| 亚洲成色77777| 亚洲一级一片aⅴ在线观看| 成年人午夜在线观看视频| 少妇被粗大的猛进出69影院| 免费播放大片免费观看视频在线观看| 在线 av 中文字幕| 精品一区二区三卡| 性高湖久久久久久久久免费观看| 99香蕉大伊视频| 各种免费的搞黄视频| 在线观看美女被高潮喷水网站| 国产av国产精品国产| 2021少妇久久久久久久久久久| 国产精品久久久久久精品古装| 久热这里只有精品99| 久久久久久久久久久免费av| 亚洲精品国产av蜜桃| 国产在线一区二区三区精| 亚洲精品一二三| 女性生殖器流出的白浆| 中文欧美无线码| 一级毛片 在线播放| 久久综合国产亚洲精品| 午夜日韩欧美国产| 啦啦啦啦在线视频资源| 一区福利在线观看| 老鸭窝网址在线观看| 日韩人妻精品一区2区三区| 18在线观看网站| 26uuu在线亚洲综合色| 青青草视频在线视频观看| 欧美日韩精品网址| 九九爱精品视频在线观看| 少妇被粗大的猛进出69影院| 久久人妻熟女aⅴ| 国产熟女午夜一区二区三区| 伦理电影免费视频| 在现免费观看毛片| 18禁观看日本| 日本黄色日本黄色录像| 日产精品乱码卡一卡2卡三| 午夜激情久久久久久久| 日韩免费高清中文字幕av| 性色av一级| 亚洲美女搞黄在线观看| 国产探花极品一区二区| 这个男人来自地球电影免费观看 | 国产福利在线免费观看视频| 十八禁高潮呻吟视频| 日本黄色日本黄色录像| 91国产中文字幕| 十分钟在线观看高清视频www| 国产av国产精品国产| 伊人久久大香线蕉亚洲五| 国产在线免费精品| 免费看不卡的av| 国产日韩欧美亚洲二区| 看免费成人av毛片| 大香蕉久久成人网| 午夜免费鲁丝| 久久久久国产精品人妻一区二区| 午夜激情久久久久久久| 日日摸夜夜添夜夜爱| 国产精品99久久99久久久不卡 | 又黄又粗又硬又大视频| 激情视频va一区二区三区| 尾随美女入室| 国产精品 国内视频| 欧美日韩亚洲国产一区二区在线观看 | 亚洲成av片中文字幕在线观看 | 国产欧美日韩综合在线一区二区| 亚洲天堂av无毛| 国产成人精品久久久久久| 欧美亚洲 丝袜 人妻 在线| 一级片免费观看大全| 777米奇影视久久| 在线观看美女被高潮喷水网站| 亚洲精品日本国产第一区| 午夜影院在线不卡| 国产免费现黄频在线看| 免费观看a级毛片全部| 在线观看美女被高潮喷水网站| 亚洲欧美中文字幕日韩二区| 在线观看美女被高潮喷水网站| 亚洲精品国产色婷婷电影| 国产成人aa在线观看| 欧美老熟妇乱子伦牲交| 亚洲精品aⅴ在线观看| 亚洲欧美成人综合另类久久久| 亚洲,一卡二卡三卡| 欧美成人精品欧美一级黄| 捣出白浆h1v1| 午夜日本视频在线| 精品人妻熟女毛片av久久网站| 久久久久久人妻| 久久久久国产网址| 日本wwww免费看| 成人18禁高潮啪啪吃奶动态图| 久久国内精品自在自线图片| 久久婷婷青草| 亚洲av电影在线进入| 亚洲av成人精品一二三区| 91精品国产国语对白视频| 欧美av亚洲av综合av国产av | 精品亚洲成国产av| 熟女电影av网| 精品久久久久久电影网| 免费av中文字幕在线| 欧美 日韩 精品 国产| 久久狼人影院| 午夜福利网站1000一区二区三区| 丝袜人妻中文字幕| 一本色道久久久久久精品综合| 男女国产视频网站| 欧美日韩视频高清一区二区三区二| 亚洲色图综合在线观看| 在线观看美女被高潮喷水网站| 国产人伦9x9x在线观看 | 黄色怎么调成土黄色| 久久久久久久精品精品| 卡戴珊不雅视频在线播放| www.精华液| 亚洲内射少妇av| 免费久久久久久久精品成人欧美视频| 国产高清国产精品国产三级| 国产精品女同一区二区软件| 国产一区二区 视频在线| 最近中文字幕高清免费大全6| 热re99久久国产66热| 校园人妻丝袜中文字幕| 久久人妻熟女aⅴ| 免费高清在线观看视频在线观看| 久久久久久伊人网av| 高清av免费在线| 在线 av 中文字幕| 日韩电影二区| 久久久久久久久久久久大奶| 亚洲av福利一区| 欧美少妇被猛烈插入视频| 国产精品一二三区在线看| 午夜免费男女啪啪视频观看| 国产极品天堂在线| 免费观看性生交大片5| 妹子高潮喷水视频| 久久av网站| 成人国产麻豆网| 久久韩国三级中文字幕| 下体分泌物呈黄色| av视频免费观看在线观看| 亚洲精品aⅴ在线观看| 国产精品免费大片| 九色亚洲精品在线播放| av国产精品久久久久影院| 美女福利国产在线| 啦啦啦啦在线视频资源| 1024香蕉在线观看| 哪个播放器可以免费观看大片| 日日撸夜夜添| 欧美人与性动交α欧美精品济南到 | 欧美黄色片欧美黄色片| 欧美日韩成人在线一区二区| 在线观看国产h片| 一边摸一边做爽爽视频免费| 美女高潮到喷水免费观看| 免费不卡的大黄色大毛片视频在线观看| 狂野欧美激情性bbbbbb| 国产高清不卡午夜福利| 日日撸夜夜添| 国产精品99久久99久久久不卡 | 宅男免费午夜| av在线观看视频网站免费| 国产精品女同一区二区软件| 精品亚洲乱码少妇综合久久| 亚洲av欧美aⅴ国产| 黄色毛片三级朝国网站| 国产日韩欧美亚洲二区| 美女脱内裤让男人舔精品视频| 亚洲欧洲日产国产| 男人爽女人下面视频在线观看| 91成人精品电影| av网站免费在线观看视频| 十分钟在线观看高清视频www| 国产 一区精品| 亚洲精品自拍成人| 日韩三级伦理在线观看| 下体分泌物呈黄色| 国产亚洲av片在线观看秒播厂| 日本免费在线观看一区| 男女边吃奶边做爰视频| 日韩制服丝袜自拍偷拍| 大话2 男鬼变身卡| 色视频在线一区二区三区| 自线自在国产av| 边亲边吃奶的免费视频| 日日摸夜夜添夜夜爱| 亚洲欧美色中文字幕在线| 国产福利在线免费观看视频| 各种免费的搞黄视频| 国产黄色视频一区二区在线观看| 免费在线观看完整版高清| 视频区图区小说| 久久国内精品自在自线图片| 在线精品无人区一区二区三| 少妇人妻 视频| 9191精品国产免费久久| 国产不卡av网站在线观看| 69精品国产乱码久久久| 高清欧美精品videossex| 91成人精品电影| 国产成人精品久久久久久| 日本wwww免费看| xxx大片免费视频| 午夜av观看不卡| xxx大片免费视频| 麻豆精品久久久久久蜜桃| 男女无遮挡免费网站观看| 午夜福利网站1000一区二区三区| 美女福利国产在线| 男人添女人高潮全过程视频| 电影成人av| 一本大道久久a久久精品| 久久这里有精品视频免费| 国产精品国产三级专区第一集| xxxhd国产人妻xxx| 亚洲国产毛片av蜜桃av| 亚洲欧洲国产日韩| 一本大道久久a久久精品| 国产欧美日韩一区二区三区在线| 在线看a的网站| 日本猛色少妇xxxxx猛交久久| 三级国产精品片| 国产欧美亚洲国产| 一级片'在线观看视频| 男人舔女人的私密视频| 国产福利在线免费观看视频| 观看美女的网站| 国产爽快片一区二区三区| 久久人妻熟女aⅴ| www.精华液| 午夜激情av网站| 99国产精品免费福利视频| 国产精品三级大全| 黄色 视频免费看| 成人毛片a级毛片在线播放| 亚洲情色 制服丝袜| 老熟女久久久| 最新的欧美精品一区二区| 久久午夜综合久久蜜桃| 三上悠亚av全集在线观看| 国产欧美日韩综合在线一区二区| 最近手机中文字幕大全| 日韩精品免费视频一区二区三区| 99香蕉大伊视频| 91在线精品国自产拍蜜月| 亚洲美女搞黄在线观看| 久久久亚洲精品成人影院| 天堂8中文在线网| 精品亚洲乱码少妇综合久久| 香蕉国产在线看| 色播在线永久视频| videos熟女内射| 最黄视频免费看| 国产精品久久久久久精品电影小说| 精品国产乱码久久久久久男人| 制服人妻中文乱码| 男人添女人高潮全过程视频| 欧美人与性动交α欧美精品济南到 | 亚洲人成电影观看| 久久综合国产亚洲精品| h视频一区二区三区| 亚洲第一区二区三区不卡| av一本久久久久| av在线app专区| 国产成人精品无人区| 男女国产视频网站| 国产乱人偷精品视频| 人成视频在线观看免费观看| 亚洲av电影在线观看一区二区三区| 在线观看免费视频网站a站| 久久久精品国产亚洲av高清涩受| 欧美成人精品欧美一级黄| 寂寞人妻少妇视频99o| 久久久久久久久免费视频了| 精品人妻在线不人妻| 精品一区二区免费观看| 亚洲精品中文字幕在线视频| 天堂中文最新版在线下载| 校园人妻丝袜中文字幕| 亚洲,欧美精品.| 啦啦啦在线观看免费高清www| 伊人久久国产一区二区| 久久免费观看电影| 亚洲欧洲精品一区二区精品久久久 | 亚洲欧美精品自产自拍| 亚洲 欧美一区二区三区| 两个人免费观看高清视频| 亚洲精品视频女| 91久久精品国产一区二区三区| 亚洲国产精品国产精品| 免费播放大片免费观看视频在线观看| 水蜜桃什么品种好| 男女下面插进去视频免费观看| 亚洲国产av新网站| 亚洲精品久久久久久婷婷小说| 爱豆传媒免费全集在线观看| 亚洲精品av麻豆狂野| 国产精品久久久久久久久免| 少妇的逼水好多| 午夜福利乱码中文字幕| 国产 一区精品| 中文字幕av电影在线播放| 超碰97精品在线观看| 女人久久www免费人成看片| 18禁裸乳无遮挡动漫免费视频| 三级国产精品片| www.av在线官网国产| 国产精品不卡视频一区二区| 欧美日韩综合久久久久久| 色哟哟·www| 老熟女久久久| 99香蕉大伊视频| 亚洲色图综合在线观看| 宅男免费午夜| 亚洲av电影在线观看一区二区三区| 男人添女人高潮全过程视频| 国产极品粉嫩免费观看在线| 亚洲国产欧美日韩在线播放| 免费高清在线观看视频在线观看| 一区二区三区激情视频| 美女大奶头黄色视频| 黄色视频在线播放观看不卡| 精品国产一区二区三区久久久樱花| 制服人妻中文乱码| 午夜福利,免费看| 国产一区二区在线观看av| 国产日韩欧美视频二区| 永久免费av网站大全| 国产亚洲午夜精品一区二区久久| 天天躁夜夜躁狠狠躁躁| 最黄视频免费看| 日韩中文字幕欧美一区二区 | xxx大片免费视频| 999精品在线视频| 男女啪啪激烈高潮av片| 少妇猛男粗大的猛烈进出视频| 久久亚洲国产成人精品v| 午夜久久久在线观看| 亚洲伊人久久精品综合| 国精品久久久久久国模美| 久久国产精品男人的天堂亚洲| 国产黄色免费在线视频| 一级a爱视频在线免费观看| 国产免费现黄频在线看| 丰满迷人的少妇在线观看| 亚洲,欧美精品.| 男女边吃奶边做爰视频| 巨乳人妻的诱惑在线观看| 黄片无遮挡物在线观看| 欧美日韩av久久| 最近中文字幕2019免费版| 9色porny在线观看| 日韩精品免费视频一区二区三区| 免费播放大片免费观看视频在线观看| 美女视频免费永久观看网站| 母亲3免费完整高清在线观看 | 日韩欧美精品免费久久| 老司机影院毛片| 有码 亚洲区| 精品亚洲成国产av| 777米奇影视久久| 天天躁狠狠躁夜夜躁狠狠躁| 各种免费的搞黄视频| 一级黄片播放器| 人成视频在线观看免费观看| 免费日韩欧美在线观看| 成人亚洲精品一区在线观看| 国产成人av激情在线播放| 99久久中文字幕三级久久日本| 精品午夜福利在线看| 9191精品国产免费久久| 午夜福利一区二区在线看| 久久久久精品久久久久真实原创| 日本色播在线视频| 日韩不卡一区二区三区视频在线| 看免费av毛片| 在线观看人妻少妇| 亚洲av中文av极速乱| 成人影院久久| 欧美 日韩 精品 国产| 国产高清不卡午夜福利| 日本wwww免费看| 免费观看性生交大片5| 国产又爽黄色视频| 超碰成人久久| 我要看黄色一级片免费的| 亚洲欧美色中文字幕在线| 人体艺术视频欧美日本| 91久久精品国产一区二区三区| 热99久久久久精品小说推荐| 搡老乐熟女国产| 国产一区二区三区综合在线观看| 丝袜美足系列| 制服丝袜香蕉在线| 狂野欧美激情性bbbbbb| 中文天堂在线官网| 国产日韩欧美视频二区| 男的添女的下面高潮视频| 欧美人与性动交α欧美软件| 久久国产精品大桥未久av| 亚洲 欧美一区二区三区| 国产白丝娇喘喷水9色精品| av在线播放精品| 老司机影院毛片| 美女福利国产在线| 久久国产精品大桥未久av| 日本av手机在线免费观看| av在线老鸭窝| 亚洲综合色惰| 高清在线视频一区二区三区| 中国三级夫妇交换| 丝瓜视频免费看黄片| 国产国语露脸激情在线看| 久久久精品区二区三区| av.在线天堂| 国产男女超爽视频在线观看| 久久久精品94久久精品| 老女人水多毛片| 一区二区av电影网| 亚洲欧美色中文字幕在线| 又大又黄又爽视频免费| 美女视频免费永久观看网站| 久久精品亚洲av国产电影网| 国产成人aa在线观看| 久久精品亚洲av国产电影网| 韩国高清视频一区二区三区| 街头女战士在线观看网站| 人妻系列 视频| 丝袜人妻中文字幕| 麻豆精品久久久久久蜜桃| 日韩一本色道免费dvd| 秋霞在线观看毛片| 最新中文字幕久久久久| 如何舔出高潮| 老鸭窝网址在线观看| 青春草国产在线视频| 日韩人妻精品一区2区三区| 欧美成人午夜精品| 母亲3免费完整高清在线观看 | 亚洲精品在线美女| 电影成人av| 亚洲欧美中文字幕日韩二区| 青草久久国产| 妹子高潮喷水视频| 国产一区二区激情短视频 | 国产av国产精品国产| 免费看不卡的av| 精品第一国产精品|