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

    基于平滑最小方差無失真響應的一致性同步算法研究

    2019-01-21 08:18:10顧廣華齊順愛李小俚
    中國生物醫(yī)學工程學報 2018年6期
    關鍵詞:電信號認知障礙頻段

    顧廣華 崔 冬* 王 娟 齊順愛 李小俚

    1(燕山大學信息科學與工程學院,河北 秦皇島 066004)2(河北省信息傳輸與信號處理重點實驗室,河北 秦皇島 066004)3(北京師范大學認知神經(jīng)科學和學習國家重點實驗室,北京 100875)

    引言

    一致性分析廣泛應用于探討腦電信號的功能連接性上,在知覺、行動和認知等方面的神經(jīng)同步方面起著重要作用[1]。電極對間的一致性值越大,說明兩路信號相互依賴的程度越強,反映大腦相應位點具有較高的功能連接性[2]。

    糖尿病是輕度認知功能障礙(mild cognitive impairment, MCI)的危險因素之一[3]。MCI是介于正常衰老和癡呆的中間狀態(tài),分為遺忘型輕度認知障礙(amnestic mild cognitive impairment, aMCI)和非遺忘型輕度認知障礙(non-amnestic mild cognitive impairment, nMCI),其中aMCI最易發(fā)展成為癡呆,每年約有10%~15% 的患者發(fā)展成癡呆[4],對aMCI盡早進行干預可預防和延緩發(fā)展成癡呆的情況。

    一致性可應用于分析靜息狀態(tài)下大腦功能連接的情況,也可以應用于檢查執(zhí)行認知任務時大腦功能連接的改變。Locatelli對比分析10例癡呆患者和10 例正常老人,發(fā)現(xiàn)癡呆患者在顳-頂區(qū)α頻段一致性下降,顳-頂-枕區(qū)兩半球間θ頻段和β頻段一致性下降,而在額葉和枕部區(qū)域間δ頻段的一致性升高[5]。Babiloni等對大腦內(nèi)部的電極對(F3、Fz、F4、C3、Cz、C4、P3、Pz、P4)進行一致性研究,發(fā)現(xiàn)AD患者在大腦兩半球間一致性下降,而在前額-頂區(qū)一致性升高[6]。Jelles等對15 例AD患者和17 例年齡匹配的正常人進行一致性研究,發(fā)現(xiàn)在α頻段和β頻段患者的一致性低于正常組[7]。Brunovsky發(fā)現(xiàn),癡呆患者的α頻段一致性下降,δ頻段的一致性增加,且這些特征和癡呆嚴重程度有關[8]。Stam發(fā)現(xiàn),AD患者在β頻段一致性下降并且和MMSE得分呈相關性,在α頻段AD患者的一致性值低于正常對照組[9]。Jelic等發(fā)現(xiàn),AD患者在顳-頂區(qū)的一致性值與記憶障礙患者和正常老年人相比下降[10]。

    本研究基于最小方差無失真響應(minimum variance distortionless response, MVDR)一致性算法[11-12],提出結合CB核平滑功率譜的平滑最小方差無失真響應算法(smoothing minimum variance distortionless response, SMVDR)。采用線性模型,在窄帶信號和寬帶信號情況下,對MVDR和SMVDR算法的性能進行仿真分析,包括抗噪性能、信號幅度。仿真結果顯示,SMVDR算法在準確性及抗噪性能方面明顯優(yōu)于MVDR算法。利用SMVDR算法,分析糖尿病患者輕度認知障礙腦電數(shù)據(jù),對aMCI組和對照組在不同頻段、不同通道組間進行獨立樣本t檢驗統(tǒng)計分析,找出有顯著特征的通道組及頻帶,并將一致性結果與認知量表進行相關性分析,探討能夠區(qū)分糖尿病人中具有輕度認知障礙和認知正?;颊叩奶卣?。

    1 方法

    1.1 MVDR一致性算法

    一致性分析可以反映兩個隨機過程在各頻率成分上的同步關系,基于譜分析描述大腦在半球間以及大腦半球內(nèi)不同區(qū)域、不同節(jié)律上電活動的同步情況[13-15]。MVDR算法是一種非參數(shù)譜的估計算法,主要是基于濾波器組的分解,令信號通過一組受限的帶通濾波器,其中每個濾波器都固定一個分析頻率,從而在每個頻帶估計信號的頻譜。

    設x(n)是零均值平穩(wěn)隨機過程,表示為

    x(n)=[x(n)x(n-1)…x(n-L+1)]T

    (1)

    式中,上標T表示向量的轉置。

    令x(n)是K個濾波器的輸入,每個濾波器的長度為L[11-12],有

    (2)

    濾波器gk的輸出yk(n)的功率為

    (3)

    式中,E{·}為數(shù)學期望,上標H表示向量或矩陣的轉置共軛,協(xié)方差矩陣Rxx為

    Rxx=E{x(n)x(n)H}

    (4)

    (5)

    在頻率ωk處x(n)的頻譜為

    (6)

    對于MVDR譜估計,每個濾波器gk的設計相當于最小化下面的函數(shù),有

    (7)

    則有

    (8)

    因此,MVDR譜為

    (9)

    對于兩個零均值平穩(wěn)隨機信號x1(n)和x2(n),兩個濾波器可表示為

    (10)

    在頻率ωk處的自譜Sxixi(ωk) (i=1,2)和交叉譜Sx1x2(ωk)為

    (i,j=1,2)

    (11)

    式中,上標*表示復共軛。

    最后,兩個信號x1(n)和x2(n)之間的MVDR一致性函數(shù)定義為

    (12)

    1.2 SMVDR一致性算法

    本研究提出新的計算一致性的算法,通過CB核平滑的最小方差無失真地響應算法MVDR頻譜。CB核廣泛應用在時頻分析中,克服了高斯核的缺點,是一種緊致支撐核[16]。CB核通過控制核的參數(shù),可以有效地抑制交叉項,提高時頻分辨率,恢復由于截斷引起的信息丟失,同時提高算法的準確性[17]。通過計算CB核和瞬時自相關函數(shù)的卷積來產(chǎn)生廣義瞬時自相關函數(shù),可以減少有效信息的損失,使能量集中,提高準確性。CB核的定義[18]如下:

    (13)

    式中:N為采樣點,n=1,2,…,N;M代表窗長,通常取奇數(shù)值,m=1,2,…,M。

    核的寬度通過參數(shù)M/B控制,峰值通過參數(shù)C調整,以上參數(shù)使CB核具有較高的分辨率和足夠的交叉相抑制。

    對于計算SMVDR一致性,首先把通過MVDR方法計算出來的自譜的平方值和交叉頻譜的平方值作為信號,再計算這3種平方譜的瞬時頻率自相關函數(shù),有

    (14)

    式中,m是頻率延遲。

    給定該CB核的窗長M,計算核與瞬時頻率自相關函數(shù)的卷積,即進行核平滑操作,定義廣義瞬時頻率自相關函數(shù)如下:

    (15)

    式中,?表示卷積操作。

    令頻率延遲m為某一常數(shù)a時,得到新的SMVDR一致性計算方法為

    這種改進的計算一致性的SMVDR方法 也是關于頻率的函數(shù),并且取值范圍同樣在0~1之間,可以評估兩信號間的相關性程度,適用于雙通道腦電信號的同步分析。

    1.3 算法仿真分析

    選取線性仿真模型,包含兩個正弦信號x1(n)和x2(n),表達形式[11]如下:

    (17)

    式中:fi是正弦信號的頻率,其中i=0,1,…,Nf-1;fs為采樣頻率;φi是隨機相位;r1(n)和r2(n)是兩個獨立的均值為0、方差為1的高斯白噪聲。

    為獲得窄帶信號,設定5個一致性頻率:f0=10,f1=15,f2=25,f3=30,f4=35 Hz。理論上,在這5個頻率點處一致性值應該接近1,而在其他頻率處一致性值應該接近0。

    濾波器個數(shù)K設為1 500,濾波器長度L設為150。這種參數(shù)設置可以保證一致性估計值較為準確,并且計算復雜度不會太高[19]。CB核中參數(shù)的選取,C=1.1,M/B=1.2,可以保證分辨率較高,并且能較好地抑制交叉項[20]。平滑操作中,M=3,a=1。

    1.3.1抗噪性能

    為了比較兩種一致性算法的抗噪性能,在線性模型產(chǎn)生的窄帶信號中加高斯白噪聲,其信噪比SNR分別為-10、-5、5和20 dB。噪聲對MVDR和SMVDR算法的影響如圖1所示??梢钥闯?,隨著信噪比的增加,兩種一致性算法性能更好,并且在5個一致性頻率處的峰值也隨著信噪比的增大而增大,如圖1 (b)~(d) 所示。SMVDR算法在一致性頻率處具有更窄的峰值,而在其他頻率處更接近0,因此該算法在抗噪性能方面要優(yōu)于MVDR算法。

    圖1 噪聲對窄帶信號一致性算法的影響。(a) SNR=-10 dB;(b) SNR=-5 dB;(c) SNR=5 dB;(d) SNR=20 dBFig.1 Effects of the noise impact on coherence methods in narrow-band signals. (a) SNR=-10 dB; (b) SNR=-5 dB; (c) SNR=5 dB; (d) SNR=20 dB

    1.3.2信號幅度

    為了觀察信號不同頻率間相對功率對一致性算法的影響,將最初選定的5個一致性頻率分成兩組,設定兩個信號在前3個頻率點f0=10、f1=15和f2=25 Hz處以及后兩個頻率點f3=30、f4=35Hz處的信號幅度之比分別為10∶1和5∶1,一致性對比的仿真圖形如圖2所示。

    如圖2(a) 所示,當兩組信號的幅度對比是10∶1時,后兩個頻率點f3=30 Hz和f4=35 Hz處的一致性值不明顯;如圖2(b) 所示,當信號幅度增大時,后兩個頻率點處的一致性值有所提高??梢?,兩種一致性算法均受信號幅度的影響。因此,在對實際腦電信號進行數(shù)據(jù)分析時,為減小信號不同頻段振幅的影響,先對實際的腦電信號進行分頻段的濾波處理,再對每一頻段的信號用SMVDR算法進行一致性計算,使分析的結果更準確。

    圖2 信號幅度對一致性算法的影響。(a) 幅度對比10∶1;(b) 幅度對比5∶1Fig.2 Effects of the amplitude on coherence methods. (a) Amplitude contrast is 10∶1; (b) Amplitude contrast is 5∶1

    1.4 腦電信號一致性分析

    1.4.1糖尿病患者資料

    本次研究的實驗對象選取在人民解放軍火箭軍總醫(yī)院就診的31 名糖尿病患者,均符合aMCI的納入及排除標準。所有受試者均被告知研究內(nèi)容并同意參加實驗,經(jīng)北京師范大學倫理委員會批準。

    記錄患者的姓名、年齡、性別、受教育程度、MOCA得分情況。依據(jù)MOCA得分對患者進行分組:MOCA分數(shù)<26為aMCI組,18 例;MOCA分數(shù)≥26為對照組,13例。運用SPSS統(tǒng)計軟件,對兩組患者的基本人口特征和神經(jīng)心理學量表進行獨立樣本t檢驗,結果用“均值±標準差”的形式表示,如表1所示??梢钥闯觯瑑山M患者在人口統(tǒng)計特征中不存在顯著差異,即年齡、受教育程度方面是匹配的,并且MOCA分數(shù)達到顯著水平。

    表1人口學特征和神經(jīng)心理學檢查統(tǒng)計結果

    Tab.1Theresultsofdemographiccharacteristicsandneuropsychologicalexamination

    因素aMCI組對照組P值年齡/歲70.44±7.2671.92±6.700.568受教育程度13.06±4.5613.38±2.730.818MOCA21.94±1.2626.92±1.190.000???

    1.4.2腦電信號采集與預處理

    采用EGI的GES300 型128 通道腦電記錄儀,試者處于閉眼靜息狀態(tài)進行靜息EEG數(shù)據(jù)采集,雙側乳突為參考電極,采樣頻率500 Hz,時長5 min,通過0~200 Hz帶通濾波,阻抗<20 kΩ。采用基于小波增強的獨立分量分析算法[21]對信號進行預處理,去除眼電等偽跡。

    1.4.3腦電信號數(shù)據(jù)分析

    本研究選擇19 路電極并劃分在5 個大腦區(qū)域,分別為前額(frontal, F)、頂區(qū)(parietal, P)、枕區(qū)(occipital, O)、左顳(left temporal, LT)、右顳(right temporal, RT),如圖3 所示。

    圖3 腦電信號19電極安置區(qū)域Fig.3 The location of nineteen EEG electrodes

    為了研究和比較aMCI組和對照組在大腦不同節(jié)律、不同腦區(qū)的一致性特征,將所選的雙通道電極對分為3個組,分別是大腦兩側電極對、大腦一側短距離電極對和大腦一側長距離電極對,具體選取情況如表2所示。

    表2腦區(qū)間電極對的選取情況

    Tab.2Theselectionofelectrodepairsbetweenbrainregions

    區(qū)域分組通道組大腦兩側Fp1Fp2,F3F4,C3C4,P3P4,O1O2,F7F8,T3T4,T5T6,F7T4,T3T6,F8T5,P4O1,F7T6,F8T3,Fp1F4,F8C3,PzT4,F4O1,F8O1,T3O2,F8Pz,C4F3,FzPz,F4T5大腦一側短距離P3O1,P4O2,F3C3,F4C4,C4P4,C3P3,F7T3,F8T4,T3T5,T4T6,F8F4大腦一側長距離F3T3,F4T4,F3T5,F4T6,F3P3,F4P4,F3O1,F4O2,T3O1,T4O2,T3P3,T4P4,T6C4,T6P6,T4C4,F8P4,Fp1O1,F8C4

    在真實腦電數(shù)據(jù)分析中,算法參數(shù)設置與仿真相同,即濾波器長度L和濾波器個數(shù)K分別取150和1 500,CB核中參數(shù)M和a分別取3 和1。實際腦電信號噪聲干擾較多,且頻率成分在各頻段分布不均勻,SMVDR算法具有較好的抗噪性能,但估計出的一致性值會受到各頻率振幅的影響。因此,在用SMVDR算法進行腦電信號數(shù)據(jù)分析前,首先對實際腦電信號進行濾波。本研究利用Gabor小波變換,將EEG信號分為Delta(1~4 Hz)、Theta(4~8 Hz),alpha(8~13 Hz)和beta(13~30 Hz)4個頻段。在每個頻段上,各電極對之間的一致性值通過下式進行計算,有

    (18)

    式中,U和L分別為指定頻段的上限頻率和下限頻率。

    數(shù)據(jù)長度截取為5 s一段,計算出每對電極的一致性值,剔除超出均值±3倍標準方差的異常值,求得的均值即為每對信號的一致性值,可用于后續(xù)的統(tǒng)計分析。

    利用SPSS 20.0統(tǒng)計軟件對數(shù)據(jù)進行統(tǒng)計分析,兩組間的比較采用獨立樣本t檢驗,P<0.05認為差異有統(tǒng)計學意義。采用皮爾森線性相關分析所有患者一致性值與MOCA測試得分的相關性,利用Bonferroni方法對顯著相關值進行了校正,以排除由于概率的原因出現(xiàn)的顯著相關,探究腦電信號一致性值與認知功能的內(nèi)在聯(lián)系。

    2 結果

    2.1 aMCI組與對照組統(tǒng)計分析

    統(tǒng)計分析aMCI組和對照組在各個頻段中有顯著特征的通道組的情況,結果如圖4所示??梢钥闯?,delta頻段aMCI組的一致性值要普遍低于對照組;theta頻段aMCI組在大腦兩側的右額-左枕(F4O1)以及右顳-左枕(F8O1)通道組間的一致性值顯著高于對照組,且在F4O1處具有較強的顯著性;alpha頻段aMCI組的一致性值都要普遍低于對照組,且在T4O2處具有較強的顯著性;beta頻段aMCI組的一致性值都要普遍高于對照組,且通道組都集中在大腦兩側,如T5T6、F8T5、F8O1和F4O1。

    圖4 各頻段aMCI組與對照組的一致性統(tǒng)計分析結果(每列上為腦地形圖,下為一致性對比分析;腦地形圖中,?表示aMCI組一致性值小于對照組,表示aMCI組一致性值大于對照組;一致性對比分析圖中,*表示P<0.05,**表示P<0.01)。(a) delta頻段;(b) theta頻段;(c) alpha頻段; (d)beta頻段Fig.4 The coherence statistic analysis results of aMCI group and normal control group in each frequency band(The top line of each row is brain map, the bottom line is the graph of coherence comparison analysis between the two groups; in brain map, ? donates that coherence values of aMCI is lower than that of NC, ★ donates that coherence values of aMCI is lager than that of NC; in coherence comparison analysis graph, * donates P<0.05, ** donates P<0.01). (a) Delta frequency band; (b)Theta frequency band; (c) Alpha frequency band; (d) Beta frequency band

    2.2 一致性與認知功能的相關性分析

    腦電信號一致性值與MOCA得分的相關性分析結果如圖5所示??梢钥闯?,在F7F8和T4T6的δ頻段一致性均與MOCA得分存在顯著的正相關;在θ頻段,F(xiàn)4O1和F8O1通道組的一致性均與MOCA得分存在顯著的負相關;在F7F8、F3O1、T4O2、F7T3、F8F4和FP1O1的α頻段,一致性均與MOCA得分存在顯著的正相關;在左右顳間T5T6、F8T5,右額-左枕F4O1和右顳-左枕F8O1的β頻段,一致性均與MOCA得分呈顯著的負相關。

    3 討論

    本研究提出基于CB核平滑的SMVDR一致性算法。在算法的仿真分析中,利用線性模型從抗噪性能、信號振幅影響進行兩種方法的對比。從仿真結果中可以看出,SMVDR算法比MVDR算法具有較好的抗噪性及準確性。

    利用SMVDR一致性算法分析糖尿病患者輕度認知障礙的腦電信號,發(fā)現(xiàn)在有顯著特征的通道組中,在δ頻段aMCI組的一致性值要普遍低于對照組,尤其是在大腦兩側的左右顳區(qū)(F7F8)δ頻段,與Babiloni研究結果相同[6]。近幾年對老年癡呆患者進行一致性研究,發(fā)現(xiàn)其中的某些頻帶內(nèi)一致性值會增大,而在某些頻帶中一致性值會減小。

    在θ頻段,aMCI組的一致性值要高于對照組。Locatelli等的研究結果也顯示,癡呆患者在大腦兩側θ頻段一致性值增大[5]。θ頻段aMCI的一致性值增大,可能與皮層下結構的皮層傳入神經(jīng)有關[21],也可能與補償機制有關,因為認知障礙損傷患者在記憶能力以及認知功能等方面都會減弱,使得在q頻段要對其進行補償[22]。這些原因可能導致某些通道組在θ頻段的一致性值升高。

    α頻段aMCI組的一致性值都要普遍低于對照組。從近些年針對老年癡呆患者腦電信號的一致性研究結果來看,Locatelli、Knott以及Stam等得出的結論表明,癡呆患者的腦電信號一致性值要低于正常對照組,尤其是在α頻段[5,9,23]。其中,Locatelli提出癡呆患者在大腦皮層的前額-枕區(qū)以及右顳-枕區(qū)的一致性值要低于對照組,本研究F3O1和T4O2通道組的研究結果與其相符[5]。有研究顯示,2型糖尿病認知障礙患者在耦合強度方面,大腦皮層的前額-左枕以及右額-右顳的α頻段一致性值要低于對照組[24],本研究Fp1O1和F8F4通道組與其研究結果相符。α頻段的一致性情況可以揭示皮質區(qū)域之間的功能連接水平,糖尿病患者中aMCI組的一致性值降低可能與不同腦區(qū)大腦皮質之間的連接能力受阻有關[21]。

    圖5 不同頻段神經(jīng)心理學量表(MOCA)與一致性值線性分析結果(每個子圖中第一行為相關性分析腦地形圖,其余為一致性值與MOCA測試得分散點圖;腦地形圖中,?表示一致性值與MOCA得分呈顯著性正相關,表示一致性值與MOCA得分呈顯著性負相關 )。(a)delta頻段;(b)theta頻段; (c)alpha頻段;(d)beta頻段Fig.5 The results of correlation analysis results of neuropsychological scale and coherence values in each frequency band(The first line of each subgraph is brain map of correlation analysis, the others are the scatter graphs of correlation between coherence values and MOCA scores; in brain map, ? donates a significant positive correlation between coherence values and MOCA scores, donates a significant negative correlation between coherence values and MOCA scores). (a) Delta frequency band;(b) Theta frequency band;(c) Alpha frequency band;(d)Beta frequency band

    β頻段aMCI組的一致性值要高于對照組。有研究表明,糖尿病認知障礙患者的腦電信號在耦合強度方面,位于大腦區(qū)域左右顳間的aMCI組一致性值要高于對照組,與本研究中左右顳間的T5T6和F8T5通道組與糖尿病認知障礙患者的腦電信號在耦合強度方面的研究結果相符合[25]。

    本研究主要對新算法SMVDR方法和原有方法做仿真分析,并將其應用于實際腦電信號的分析中。仿真實驗和實際數(shù)據(jù)的分析分別得出了一些結論,但研究中尚存在一些局限性:樣本數(shù)量較少,只有31人;通道數(shù)目較少,只選取了19路腦電信號等。在以后的工作中,將對此做進一步分析。

    4 結論

    本研究提出了一種新的計算雙通道腦電信號的一致性算法——SMVDR方法,該方法具有較好的抗噪性能及準確性。通過對糖尿病患者輕度認知障礙的腦電信號一致性的統(tǒng)計分析,發(fā)現(xiàn)這些患者在δ頻段的左右顳間一致性下降,在θ頻段的前額-枕區(qū)一致性增加,在α頻段的右顳-枕區(qū)和前額-右顳區(qū)域一致性下降,在β頻段的左右顳間一致性增加。通過一致性值與MOCA得分的相關性分析,發(fā)現(xiàn)在特征通道下δ和α頻段的一致性與MOCA得分存在顯著的正相關,θ和β頻段的一致性與MOCA分數(shù)呈負相關。SMVDR方法以及相應節(jié)律下的不同腦區(qū)之間的一致性特征,可以輔助醫(yī)生盡早地診斷糖尿病輕度認知障礙,對預測輕度認知障礙并及時進行干預具有重要意義。

    猜你喜歡
    電信號認知障礙頻段
    防跌倒,警惕認知障礙性疾病
    基于聯(lián)合聚類分析的單通道腹部心電信號的胎心率提取
    關愛父母,關注老年認知障礙癥
    中老年保健(2021年6期)2021-08-24 06:54:06
    gPhone重力儀的面波頻段響應實測研究
    地震研究(2021年1期)2021-04-13 01:04:56
    慢性心力衰竭與認知障礙的相關性
    基于Code Composer Studio3.3完成對心電信號的去噪
    科技傳播(2019年24期)2019-06-15 09:29:28
    基于隨機森林的航天器電信號多分類識別方法
    推擠的5GHz頻段
    CHIP新電腦(2016年3期)2016-03-10 14:07:52
    TD—LTE在D頻段和F頻段的覆蓋能力差異
    中國新通信(2015年1期)2015-05-30 10:30:46
    老年抑郁癥與輕度認知障礙的相關性研究
    亚洲黑人精品在线| 日韩免费av在线播放| 国产精品久久久人人做人人爽| av在线天堂中文字幕| 麻豆国产av国片精品| www.熟女人妻精品国产| 欧美日韩精品网址| 男女下面进入的视频免费午夜| 免费高清视频大片| 给我免费播放毛片高清在线观看| 最新美女视频免费是黄的| 亚洲av五月六月丁香网| 亚洲真实伦在线观看| 老汉色av国产亚洲站长工具| 午夜精品在线福利| 黄片大片在线免费观看| 精品久久久久久久久久免费视频| 亚洲国产高清在线一区二区三| 久久精品夜夜夜夜夜久久蜜豆| 一个人免费在线观看电影| 亚洲人与动物交配视频| 中文字幕久久专区| 伊人久久大香线蕉亚洲五| 天堂√8在线中文| 嫁个100分男人电影在线观看| 日本在线视频免费播放| 在线a可以看的网站| 在线观看一区二区三区| 日本 欧美在线| 99久国产av精品| netflix在线观看网站| 亚洲专区中文字幕在线| 亚洲av电影在线进入| 国产精品美女特级片免费视频播放器| 免费一级毛片在线播放高清视频| 此物有八面人人有两片| 深爱激情五月婷婷| 夜夜夜夜夜久久久久| 免费av观看视频| 嫩草影视91久久| 国产三级黄色录像| 狂野欧美白嫩少妇大欣赏| 窝窝影院91人妻| av在线蜜桃| 日本精品一区二区三区蜜桃| 美女cb高潮喷水在线观看| 欧美大码av| 精品不卡国产一区二区三区| 国产精品一及| 国产欧美日韩一区二区精品| 真人做人爱边吃奶动态| 午夜福利在线在线| 久久久精品欧美日韩精品| 美女 人体艺术 gogo| 免费一级毛片在线播放高清视频| 午夜影院日韩av| 俄罗斯特黄特色一大片| 免费看日本二区| 久久精品国产亚洲av涩爱 | eeuss影院久久| 搡老妇女老女人老熟妇| 内地一区二区视频在线| 熟女电影av网| 老汉色av国产亚洲站长工具| 欧美色欧美亚洲另类二区| 色综合婷婷激情| 男人舔奶头视频| 琪琪午夜伦伦电影理论片6080| 国产真实乱freesex| 啦啦啦韩国在线观看视频| 国内少妇人妻偷人精品xxx网站| 深夜精品福利| 一进一出抽搐动态| 日韩国内少妇激情av| 国产真人三级小视频在线观看| 观看美女的网站| 老司机午夜福利在线观看视频| 伊人久久精品亚洲午夜| 日本五十路高清| 欧美一区二区亚洲| 国产久久久一区二区三区| 亚洲av五月六月丁香网| 亚洲国产欧洲综合997久久,| 国产精华一区二区三区| 天堂网av新在线| ponron亚洲| 久久精品综合一区二区三区| 亚洲精品456在线播放app | 一级毛片高清免费大全| 亚洲真实伦在线观看| 黄色视频,在线免费观看| xxxwww97欧美| 亚洲va日本ⅴa欧美va伊人久久| 丁香欧美五月| 18美女黄网站色大片免费观看| 中文字幕人妻熟人妻熟丝袜美 | 午夜精品久久久久久毛片777| 国产精品电影一区二区三区| 亚洲乱码一区二区免费版| 在线观看美女被高潮喷水网站 | 日本一本二区三区精品| 欧美又色又爽又黄视频| 亚洲国产精品成人综合色| 内地一区二区视频在线| 我的老师免费观看完整版| 有码 亚洲区| x7x7x7水蜜桃| 成人高潮视频无遮挡免费网站| 国产一区二区在线观看日韩 | 久久精品国产亚洲av涩爱 | 母亲3免费完整高清在线观看| 日韩欧美在线二视频| 国产精品亚洲一级av第二区| 久久久久久久午夜电影| 久久草成人影院| 欧美日韩精品网址| 九色成人免费人妻av| 午夜激情欧美在线| 成人国产综合亚洲| 亚洲精品美女久久久久99蜜臀| 别揉我奶头~嗯~啊~动态视频| 国产av麻豆久久久久久久| 人人妻,人人澡人人爽秒播| 少妇裸体淫交视频免费看高清| 午夜福利成人在线免费观看| 国产午夜福利久久久久久| 久久香蕉国产精品| 蜜桃亚洲精品一区二区三区| 亚洲国产高清在线一区二区三| 有码 亚洲区| 成年免费大片在线观看| 99久久九九国产精品国产免费| 高清日韩中文字幕在线| 欧美日韩一级在线毛片| 亚洲最大成人中文| 黄片小视频在线播放| 亚洲国产色片| 99国产精品一区二区三区| 亚洲片人在线观看| 国产淫片久久久久久久久 | 一个人看的www免费观看视频| 国产私拍福利视频在线观看| 精品无人区乱码1区二区| 欧美黄色淫秽网站| h日本视频在线播放| 99热只有精品国产| 97超级碰碰碰精品色视频在线观看| 国产免费男女视频| 亚洲,欧美精品.| 一进一出抽搐动态| 老司机在亚洲福利影院| 亚洲精品国产精品久久久不卡| xxxwww97欧美| 九九在线视频观看精品| 国产一区二区三区在线臀色熟女| 日韩欧美在线二视频| 亚洲一区二区三区色噜噜| 淫秽高清视频在线观看| 精品久久久久久成人av| 女同久久另类99精品国产91| 舔av片在线| 精品99又大又爽又粗少妇毛片 | 老司机在亚洲福利影院| 亚洲精品粉嫩美女一区| 欧美另类亚洲清纯唯美| 一进一出抽搐gif免费好疼| 午夜亚洲福利在线播放| 天堂网av新在线| 国产三级在线视频| 国产亚洲av嫩草精品影院| 看片在线看免费视频| 欧美av亚洲av综合av国产av| 男人和女人高潮做爰伦理| 小蜜桃在线观看免费完整版高清| 热99在线观看视频| 中文字幕精品亚洲无线码一区| 国产免费男女视频| 99精品欧美一区二区三区四区| 久久久久性生活片| 日本成人三级电影网站| 嫩草影院精品99| 欧洲精品卡2卡3卡4卡5卡区| 人人妻,人人澡人人爽秒播| 母亲3免费完整高清在线观看| 国产探花在线观看一区二区| 日韩 欧美 亚洲 中文字幕| 免费观看的影片在线观看| 午夜日韩欧美国产| 一进一出好大好爽视频| 精品无人区乱码1区二区| 黄色日韩在线| av在线天堂中文字幕| 尤物成人国产欧美一区二区三区| 好男人电影高清在线观看| 女人十人毛片免费观看3o分钟| 日韩大尺度精品在线看网址| 久久亚洲精品不卡| 亚洲成人久久爱视频| 久久久久国产精品人妻aⅴ院| 亚洲熟妇中文字幕五十中出| 老司机午夜福利在线观看视频| 欧美三级亚洲精品| 狠狠狠狠99中文字幕| 亚洲在线观看片| 成人三级黄色视频| 可以在线观看的亚洲视频| 婷婷亚洲欧美| 最近视频中文字幕2019在线8| 亚洲不卡免费看| 国产乱人视频| 一本久久中文字幕| 久久精品国产亚洲av涩爱 | 1024手机看黄色片| 日韩精品中文字幕看吧| 久久精品影院6| 午夜福利成人在线免费观看| 老熟妇仑乱视频hdxx| а√天堂www在线а√下载| 久久久精品大字幕| 久久精品夜夜夜夜夜久久蜜豆| 日韩有码中文字幕| 国产精品98久久久久久宅男小说| 超碰av人人做人人爽久久 | eeuss影院久久| 国产真实乱freesex| 激情在线观看视频在线高清| 精品久久久久久久毛片微露脸| 高清在线国产一区| 女人高潮潮喷娇喘18禁视频| 国产真实乱freesex| 免费无遮挡裸体视频| 国产精品三级大全| 日本一本二区三区精品| eeuss影院久久| 黄色片一级片一级黄色片| 亚洲精品一卡2卡三卡4卡5卡| 欧美性猛交黑人性爽| 99久久成人亚洲精品观看| 精品人妻1区二区| 嫁个100分男人电影在线观看| 成人永久免费在线观看视频| 琪琪午夜伦伦电影理论片6080| 尤物成人国产欧美一区二区三区| 亚洲精品一区av在线观看| 亚洲最大成人手机在线| 欧美激情久久久久久爽电影| 嫩草影院精品99| 日本一本二区三区精品| 九九久久精品国产亚洲av麻豆| 中文字幕精品亚洲无线码一区| 中文亚洲av片在线观看爽| 男女之事视频高清在线观看| 欧美高清成人免费视频www| 国产精品一及| 欧洲精品卡2卡3卡4卡5卡区| 91av网一区二区| 操出白浆在线播放| 亚洲精品在线美女| 99riav亚洲国产免费| 日本熟妇午夜| 舔av片在线| 欧美一级a爱片免费观看看| 一区福利在线观看| 国产成人欧美在线观看| 午夜福利在线在线| 一进一出抽搐动态| 日韩欧美 国产精品| 一个人观看的视频www高清免费观看| 亚洲av中文字字幕乱码综合| 亚洲18禁久久av| 亚洲人成网站在线播| 国产三级中文精品| 精品熟女少妇八av免费久了| 精品一区二区三区人妻视频| 欧美黑人巨大hd| 亚洲成av人片在线播放无| 国产高清视频在线观看网站| 亚洲av不卡在线观看| 伊人久久大香线蕉亚洲五| 免费观看的影片在线观看| 手机成人av网站| 欧美性猛交╳xxx乱大交人| 国产成人a区在线观看| 久久精品国产自在天天线| 在线看三级毛片| 两个人的视频大全免费| 欧美成狂野欧美在线观看| 波野结衣二区三区在线 | 99久久精品一区二区三区| 国产精品,欧美在线| 99久久综合精品五月天人人| 国产午夜精品久久久久久一区二区三区 | 美女cb高潮喷水在线观看| 婷婷亚洲欧美| 亚洲av第一区精品v没综合| 麻豆一二三区av精品| 色吧在线观看| 99久久成人亚洲精品观看| 五月伊人婷婷丁香| 亚洲在线观看片| 丰满的人妻完整版| 国产精品 欧美亚洲| 亚洲成人精品中文字幕电影| 久久久久性生活片| 午夜福利在线在线| 欧美日韩福利视频一区二区| av在线蜜桃| 99riav亚洲国产免费| 成人国产综合亚洲| 国产高清有码在线观看视频| 欧美精品啪啪一区二区三区| 欧美日韩乱码在线| 欧美绝顶高潮抽搐喷水| 午夜亚洲福利在线播放| 搡老熟女国产l中国老女人| 亚洲精品粉嫩美女一区| 丁香六月欧美| 国产一区二区亚洲精品在线观看| 国产伦精品一区二区三区四那| 国产精品一及| 真实男女啪啪啪动态图| 美女高潮喷水抽搐中文字幕| 黑人欧美特级aaaaaa片| 日本精品一区二区三区蜜桃| 亚洲在线观看片| 国产精品,欧美在线| 日本在线视频免费播放| 丰满的人妻完整版| 欧美大码av| 久久婷婷人人爽人人干人人爱| 人妻久久中文字幕网| 欧美日本亚洲视频在线播放| 欧美丝袜亚洲另类 | 亚洲国产高清在线一区二区三| 老鸭窝网址在线观看| 在线十欧美十亚洲十日本专区| 中文字幕熟女人妻在线| 在线播放无遮挡| 国产亚洲欧美在线一区二区| 一区二区三区国产精品乱码| 琪琪午夜伦伦电影理论片6080| 免费在线观看日本一区| 亚洲熟妇熟女久久| 国产精品电影一区二区三区| 久久久久久久午夜电影| 国产亚洲精品综合一区在线观看| 亚洲av日韩精品久久久久久密| 成人欧美大片| 午夜福利免费观看在线| 真人做人爱边吃奶动态| 亚洲av电影在线进入| 少妇的逼水好多| 久久久久久久午夜电影| 9191精品国产免费久久| 免费看美女性在线毛片视频| 亚洲va日本ⅴa欧美va伊人久久| 香蕉丝袜av| 99久久99久久久精品蜜桃| 淫秽高清视频在线观看| 好看av亚洲va欧美ⅴa在| 亚洲性夜色夜夜综合| 免费av观看视频| 精品久久久久久久毛片微露脸| 中文字幕av在线有码专区| 亚洲精品粉嫩美女一区| 色综合婷婷激情| 噜噜噜噜噜久久久久久91| 久久久久久久久中文| 国产三级黄色录像| 在线播放无遮挡| 国产精品自产拍在线观看55亚洲| 中文资源天堂在线| 国产伦人伦偷精品视频| 岛国在线观看网站| 久久久国产精品麻豆| 天堂影院成人在线观看| 欧美av亚洲av综合av国产av| 一进一出抽搐动态| 午夜两性在线视频| 久久久久久大精品| 午夜久久久久精精品| www.色视频.com| 国产在视频线在精品| 他把我摸到了高潮在线观看| 国产aⅴ精品一区二区三区波| 日韩欧美一区二区三区在线观看| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 欧美黄色片欧美黄色片| 精品国内亚洲2022精品成人| 欧美日韩中文字幕国产精品一区二区三区| 日韩欧美在线乱码| 深爱激情五月婷婷| 激情在线观看视频在线高清| 性色avwww在线观看| 欧美三级亚洲精品| 国产精品美女特级片免费视频播放器| av中文乱码字幕在线| 最新在线观看一区二区三区| 国产真人三级小视频在线观看| 亚洲成人久久爱视频| 日韩中文字幕欧美一区二区| 69av精品久久久久久| 18+在线观看网站| 亚洲av免费在线观看| 国产精品久久电影中文字幕| 国产精品影院久久| 午夜免费激情av| 国内精品一区二区在线观看| 成人特级av手机在线观看| 在线国产一区二区在线| 51午夜福利影视在线观看| 少妇熟女aⅴ在线视频| 啪啪无遮挡十八禁网站| 国产视频一区二区在线看| 一本久久中文字幕| ponron亚洲| 免费在线观看成人毛片| 亚洲国产精品999在线| 久久草成人影院| 一区福利在线观看| 制服人妻中文乱码| 天堂√8在线中文| 波多野结衣巨乳人妻| 欧美性感艳星| 国产毛片a区久久久久| 亚洲精华国产精华精| 在线天堂最新版资源| 久久九九热精品免费| ponron亚洲| 国产日本99.免费观看| 亚洲精品日韩av片在线观看 | 亚洲精品一区av在线观看| 狂野欧美激情性xxxx| 亚洲片人在线观看| 淫妇啪啪啪对白视频| 国产亚洲欧美在线一区二区| 欧美绝顶高潮抽搐喷水| 母亲3免费完整高清在线观看| 久久香蕉精品热| 99精品久久久久人妻精品| 亚洲精品亚洲一区二区| 91麻豆av在线| 丰满人妻一区二区三区视频av | 宅男免费午夜| 国产精品野战在线观看| 搡老妇女老女人老熟妇| 女同久久另类99精品国产91| 亚洲人成电影免费在线| 无遮挡黄片免费观看| 人妻久久中文字幕网| 午夜免费观看网址| 丁香六月欧美| 亚洲av第一区精品v没综合| av天堂中文字幕网| 午夜精品久久久久久毛片777| 精品人妻1区二区| 欧美日韩亚洲国产一区二区在线观看| 岛国视频午夜一区免费看| 成人欧美大片| 51午夜福利影视在线观看| av福利片在线观看| 在线观看日韩欧美| 美女大奶头视频| 亚洲avbb在线观看| 一二三四社区在线视频社区8| 精品久久久久久久人妻蜜臀av| 亚洲无线在线观看| 99国产精品一区二区三区| 国模一区二区三区四区视频| 久久久久久人人人人人| 亚洲第一欧美日韩一区二区三区| 午夜免费观看网址| 一级毛片高清免费大全| 丝袜美腿在线中文| 人妻夜夜爽99麻豆av| e午夜精品久久久久久久| 搡女人真爽免费视频火全软件 | 久久久成人免费电影| 精品国内亚洲2022精品成人| 一个人免费在线观看电影| 日本成人三级电影网站| 日韩亚洲欧美综合| 色在线成人网| 午夜久久久久精精品| 国产一区二区亚洲精品在线观看| 亚洲国产精品合色在线| 国产精品嫩草影院av在线观看 | 亚洲成av人片在线播放无| 91久久精品国产一区二区成人 | 午夜福利在线观看吧| 每晚都被弄得嗷嗷叫到高潮| 性欧美人与动物交配| 久久久成人免费电影| 丁香六月欧美| 亚洲欧美日韩卡通动漫| 免费在线观看影片大全网站| 国产单亲对白刺激| 亚洲av免费在线观看| 亚洲成av人片在线播放无| 男女视频在线观看网站免费| 久久久久久九九精品二区国产| 欧美黄色淫秽网站| 99热只有精品国产| 搡老熟女国产l中国老女人| 久久天躁狠狠躁夜夜2o2o| 国内少妇人妻偷人精品xxx网站| 国产v大片淫在线免费观看| 国产精品免费一区二区三区在线| 亚洲国产欧美网| 亚洲国产精品成人综合色| 我的老师免费观看完整版| 少妇的丰满在线观看| 高清日韩中文字幕在线| 亚洲国产日韩欧美精品在线观看 | 悠悠久久av| 国产精品久久久久久久久免 | 欧美性感艳星| 午夜两性在线视频| 99久久精品国产亚洲精品| 免费av毛片视频| 午夜免费男女啪啪视频观看 | 不卡一级毛片| 久久草成人影院| 日本精品一区二区三区蜜桃| 欧美绝顶高潮抽搐喷水| 男人和女人高潮做爰伦理| 精品久久久久久,| 男女那种视频在线观看| 制服丝袜大香蕉在线| 精品国产亚洲在线| 国产单亲对白刺激| 一二三四社区在线视频社区8| 伊人久久精品亚洲午夜| 麻豆国产av国片精品| www日本黄色视频网| 亚洲在线观看片| 亚洲av一区综合| 熟女电影av网| 久久久国产成人精品二区| 美女黄网站色视频| 又黄又粗又硬又大视频| 亚洲成av人片免费观看| 精华霜和精华液先用哪个| 天堂网av新在线| 国产高清有码在线观看视频| 老熟妇仑乱视频hdxx| 亚洲国产中文字幕在线视频| 亚洲无线在线观看| 国产视频一区二区在线看| 中文字幕精品亚洲无线码一区| 黄色女人牲交| 欧美乱码精品一区二区三区| 高清毛片免费观看视频网站| 午夜福利18| 夜夜躁狠狠躁天天躁| 69人妻影院| 老熟妇仑乱视频hdxx| 国产久久久一区二区三区| 中文字幕人妻丝袜一区二区| 99久久精品一区二区三区| 亚洲精品成人久久久久久| 免费观看人在逋| 免费看光身美女| 欧美zozozo另类| 日本成人三级电影网站| 桃红色精品国产亚洲av| 国内精品美女久久久久久| 亚洲成人久久爱视频| 国产av在哪里看| 人妻久久中文字幕网| 国产91精品成人一区二区三区| 啦啦啦观看免费观看视频高清| 欧美一区二区国产精品久久精品| 日韩欧美在线乱码| 国产高清有码在线观看视频| av福利片在线观看| 久久久色成人| 99国产极品粉嫩在线观看| 亚洲欧美精品综合久久99| 久久久成人免费电影| 怎么达到女性高潮| 最近视频中文字幕2019在线8| 亚洲最大成人中文| 性色avwww在线观看| 日韩欧美 国产精品| 国产视频内射| 色噜噜av男人的天堂激情| 国产 一区 欧美 日韩| 天堂√8在线中文| 国产亚洲精品综合一区在线观看| 最近最新中文字幕大全免费视频| 在线观看66精品国产| 精品欧美国产一区二区三| 国产爱豆传媒在线观看| 日本黄色片子视频| 搡女人真爽免费视频火全软件 | 99久久九九国产精品国产免费| 在线免费观看不下载黄p国产 | 搡老岳熟女国产| 麻豆一二三区av精品| 嫁个100分男人电影在线观看| 五月伊人婷婷丁香| 人人妻人人看人人澡| 欧美成狂野欧美在线观看| 久久亚洲精品不卡| 国产伦在线观看视频一区| 亚洲国产色片| 亚洲激情在线av| 亚洲最大成人中文| 99热6这里只有精品| 中文字幕人妻丝袜一区二区| 国语自产精品视频在线第100页| 色视频www国产| 国产一区二区激情短视频| 欧美日韩综合久久久久久 |