• <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
    老年抑郁癥與輕度認知障礙的相關性研究
    纯流量卡能插随身wifi吗| 一级,二级,三级黄色视频| 人人妻,人人澡人人爽秒播| 亚洲精品成人av观看孕妇| 亚洲成人免费电影在线观看| 欧美日韩黄片免| 在线观看免费视频日本深夜| 亚洲午夜精品一区,二区,三区| 日韩免费高清中文字幕av| 咕卡用的链子| 一级毛片女人18水好多| 国产1区2区3区精品| 欧美久久黑人一区二区| 成人影院久久| 国产精品永久免费网站| 亚洲成a人片在线一区二区| 国产区一区二久久| 亚洲五月婷婷丁香| 欧美乱妇无乱码| 超色免费av| 国产成人精品久久二区二区91| 黑丝袜美女国产一区| 黄色视频,在线免费观看| 久久久久精品人妻al黑| 午夜福利乱码中文字幕| 日韩精品免费视频一区二区三区| 黑人巨大精品欧美一区二区蜜桃| 欧美日韩黄片免| 国产精品一区二区精品视频观看| 亚洲一区中文字幕在线| 亚洲午夜理论影院| 99精品在免费线老司机午夜| 日日摸夜夜添夜夜添小说| 老司机福利观看| 两性夫妻黄色片| 欧美激情高清一区二区三区| 色94色欧美一区二区| 极品教师在线免费播放| 成人免费观看视频高清| 啦啦啦 在线观看视频| 成人免费观看视频高清| 丝袜美腿诱惑在线| 久久久久久久久久久久大奶| 一二三四社区在线视频社区8| 国产欧美日韩一区二区三区在线| 9热在线视频观看99| 两性午夜刺激爽爽歪歪视频在线观看 | 久久性视频一级片| 丝袜人妻中文字幕| 交换朋友夫妻互换小说| 免费观看人在逋| 精品久久久久久久毛片微露脸| 国产野战对白在线观看| 欧美日韩av久久| 美女福利国产在线| 国产单亲对白刺激| 久99久视频精品免费| 男人的好看免费观看在线视频 | 国产1区2区3区精品| 啦啦啦在线免费观看视频4| 欧美久久黑人一区二区| 国产人伦9x9x在线观看| www.熟女人妻精品国产| 欧美日韩一级在线毛片| 黄片大片在线免费观看| 国产伦人伦偷精品视频| 午夜免费观看网址| 三上悠亚av全集在线观看| 欧美日韩av久久| 色婷婷久久久亚洲欧美| 国产成+人综合+亚洲专区| 精品卡一卡二卡四卡免费| 啦啦啦在线免费观看视频4| 免费观看精品视频网站| 国产99久久九九免费精品| 老司机福利观看| 中亚洲国语对白在线视频| 国内毛片毛片毛片毛片毛片| www.999成人在线观看| 欧美av亚洲av综合av国产av| 午夜激情av网站| 女人爽到高潮嗷嗷叫在线视频| 美女福利国产在线| 午夜日韩欧美国产| 国产免费男女视频| 欧美成人免费av一区二区三区 | 欧美成人免费av一区二区三区 | av一本久久久久| 国产精品免费大片| 制服人妻中文乱码| av国产精品久久久久影院| 夫妻午夜视频| 亚洲综合色网址| 午夜成年电影在线免费观看| 18禁裸乳无遮挡动漫免费视频| 免费不卡黄色视频| 国产欧美日韩综合在线一区二区| 亚洲人成伊人成综合网2020| 国产精华一区二区三区| av网站免费在线观看视频| 99精国产麻豆久久婷婷| 丰满的人妻完整版| 99热只有精品国产| 天堂√8在线中文| av有码第一页| 男女下面插进去视频免费观看| videosex国产| 男女之事视频高清在线观看| 一区在线观看完整版| 可以免费在线观看a视频的电影网站| 亚洲国产中文字幕在线视频| netflix在线观看网站| 日本五十路高清| 中文字幕人妻丝袜制服| 国内久久婷婷六月综合欲色啪| 建设人人有责人人尽责人人享有的| 日本精品一区二区三区蜜桃| 国产亚洲欧美98| 欧美精品啪啪一区二区三区| 18禁国产床啪视频网站| 日本wwww免费看| 久久午夜综合久久蜜桃| 国产精品乱码一区二三区的特点 | 一二三四在线观看免费中文在| 丝袜在线中文字幕| 国产午夜精品久久久久久| 五月开心婷婷网| 欧美日韩中文字幕国产精品一区二区三区 | 国产亚洲精品久久久久久毛片 | 亚洲av成人av| 成人亚洲精品一区在线观看| 成人免费观看视频高清| 老司机影院毛片| 高清黄色对白视频在线免费看| 久久久久国内视频| 淫妇啪啪啪对白视频| 大型黄色视频在线免费观看| 老鸭窝网址在线观看| av天堂久久9| 女人被躁到高潮嗷嗷叫费观| 丝袜美腿诱惑在线| xxx96com| 黑人巨大精品欧美一区二区蜜桃| 夜夜爽天天搞| 黄片播放在线免费| 村上凉子中文字幕在线| 日韩大码丰满熟妇| 女人爽到高潮嗷嗷叫在线视频| 午夜福利欧美成人| 精品熟女少妇八av免费久了| 美女国产高潮福利片在线看| 天堂√8在线中文| 亚洲免费av在线视频| a级毛片黄视频| 亚洲成人国产一区在线观看| 18禁黄网站禁片午夜丰满| 午夜福利影视在线免费观看| 国产精品免费一区二区三区在线 | 人人妻人人添人人爽欧美一区卜| 9191精品国产免费久久| 一本大道久久a久久精品| 激情视频va一区二区三区| 女人被躁到高潮嗷嗷叫费观| 国产在线观看jvid| 一本一本久久a久久精品综合妖精| 香蕉丝袜av| 亚洲人成电影观看| 精品无人区乱码1区二区| xxxhd国产人妻xxx| 老司机在亚洲福利影院| 麻豆成人av在线观看| 国产亚洲一区二区精品| 欧美+亚洲+日韩+国产| 天天躁狠狠躁夜夜躁狠狠躁| 一级a爱片免费观看的视频| 99热只有精品国产| 久久中文字幕人妻熟女| 视频在线观看一区二区三区| 一进一出抽搐动态| 国产麻豆69| 高清av免费在线| av天堂久久9| 热re99久久国产66热| 亚洲自偷自拍图片 自拍| 曰老女人黄片| 午夜激情av网站| 黄色 视频免费看| 日韩精品免费视频一区二区三区| 精品乱码久久久久久99久播| 亚洲精品美女久久av网站| 岛国毛片在线播放| 免费女性裸体啪啪无遮挡网站| 亚洲精品中文字幕一二三四区| 久久久国产成人免费| 国产高清videossex| 久久国产亚洲av麻豆专区| 黄色成人免费大全| 亚洲精品在线观看二区| 深夜精品福利| 高清黄色对白视频在线免费看| 精品久久久久久电影网| av国产精品久久久久影院| 亚洲精品中文字幕在线视频| 国产av一区二区精品久久| 黄色视频,在线免费观看| 后天国语完整版免费观看| 久久久久精品国产欧美久久久| 两性午夜刺激爽爽歪歪视频在线观看 | 岛国毛片在线播放| 国产精品成人在线| 中文字幕最新亚洲高清| 9热在线视频观看99| 男人舔女人的私密视频| 国产精品免费一区二区三区在线 | 法律面前人人平等表现在哪些方面| 中文字幕精品免费在线观看视频| 在线视频色国产色| 国产成人精品在线电影| 后天国语完整版免费观看| √禁漫天堂资源中文www| 在线国产一区二区在线| 成人av一区二区三区在线看| 精品乱码久久久久久99久播| 亚洲五月婷婷丁香| 国内久久婷婷六月综合欲色啪| 女人爽到高潮嗷嗷叫在线视频| 国产精品99久久99久久久不卡| 少妇粗大呻吟视频| 国产真人三级小视频在线观看| 亚洲精品国产精品久久久不卡| 欧美精品一区二区免费开放| 搡老乐熟女国产| 黑人巨大精品欧美一区二区mp4| 一区二区日韩欧美中文字幕| 免费在线观看黄色视频的| 欧美国产精品va在线观看不卡| 一二三四社区在线视频社区8| 777米奇影视久久| 欧美 日韩 精品 国产| 国产成人精品无人区| 日日爽夜夜爽网站| 另类亚洲欧美激情| 高清av免费在线| www.自偷自拍.com| 一区在线观看完整版| 亚洲一区中文字幕在线| 女性生殖器流出的白浆| 搡老乐熟女国产| 波多野结衣av一区二区av| 中文欧美无线码| 超碰成人久久| 久久中文字幕人妻熟女| 美女福利国产在线| 亚洲精品美女久久av网站| 国产精品免费大片| 丁香欧美五月| 极品人妻少妇av视频| 精品国产超薄肉色丝袜足j| 中文字幕人妻丝袜制服| 91老司机精品| 自线自在国产av| 欧美日韩乱码在线| 久久这里只有精品19| 两个人看的免费小视频| 一二三四在线观看免费中文在| 悠悠久久av| 色94色欧美一区二区| 狂野欧美激情性xxxx| 又黄又粗又硬又大视频| av天堂在线播放| 欧美日韩乱码在线| 极品人妻少妇av视频| 国产成人av激情在线播放| 老司机在亚洲福利影院| 色在线成人网| 国产欧美亚洲国产| 最新的欧美精品一区二区| 国内久久婷婷六月综合欲色啪| 欧美乱色亚洲激情| av线在线观看网站| 亚洲精品美女久久av网站| 国产精品秋霞免费鲁丝片| 在线国产一区二区在线| 亚洲精品在线观看二区| 精品国产美女av久久久久小说| 夜夜爽天天搞| 国产人伦9x9x在线观看| 又紧又爽又黄一区二区| 亚洲精品一二三| 成人黄色视频免费在线看| 久久国产亚洲av麻豆专区| videos熟女内射| а√天堂www在线а√下载 | 亚洲五月色婷婷综合| 国产高清激情床上av| 最近最新中文字幕大全免费视频| 欧美一级毛片孕妇| 九色亚洲精品在线播放| 日韩欧美三级三区| 久99久视频精品免费| 丝袜美足系列| 丝瓜视频免费看黄片| 国产男女超爽视频在线观看| aaaaa片日本免费| 高清欧美精品videossex| 亚洲片人在线观看| 最新的欧美精品一区二区| 最近最新中文字幕大全电影3 | 国产精品免费大片| 午夜福利影视在线免费观看| 精品一区二区三卡| 大型av网站在线播放| e午夜精品久久久久久久| 亚洲欧洲精品一区二区精品久久久| 欧美乱码精品一区二区三区| 午夜视频精品福利| 亚洲熟妇熟女久久| 国产深夜福利视频在线观看| 精品人妻在线不人妻| 天堂中文最新版在线下载| 99热网站在线观看| 王馨瑶露胸无遮挡在线观看| 国产精品影院久久| 90打野战视频偷拍视频| 三级毛片av免费| 黑人操中国人逼视频| 18禁观看日本| 久久久久精品人妻al黑| 成人黄色视频免费在线看| 久久 成人 亚洲| 免费在线观看影片大全网站| 99在线人妻在线中文字幕 | 亚洲自偷自拍图片 自拍| 在线永久观看黄色视频| 一级片免费观看大全| 国产人伦9x9x在线观看| 亚洲精品中文字幕在线视频| 人人妻人人爽人人添夜夜欢视频| 桃红色精品国产亚洲av| 伦理电影免费视频| 中文字幕制服av| 女性生殖器流出的白浆| 老司机午夜福利在线观看视频| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲,欧美精品.| 欧美日韩视频精品一区| 9色porny在线观看| 美女视频免费永久观看网站| 90打野战视频偷拍视频| 正在播放国产对白刺激| 亚洲精品国产区一区二| 深夜精品福利| 黄片大片在线免费观看| 欧美精品av麻豆av| 777久久人妻少妇嫩草av网站| 久久国产乱子伦精品免费另类| 国产伦人伦偷精品视频| 久久久国产一区二区| 在线十欧美十亚洲十日本专区| 一级a爱片免费观看的视频| 中文字幕另类日韩欧美亚洲嫩草| 老鸭窝网址在线观看| 高清av免费在线| 又黄又爽又免费观看的视频| 日韩精品免费视频一区二区三区| 欧美日本中文国产一区发布| 天天操日日干夜夜撸| 老司机福利观看| 日韩欧美在线二视频 | 亚洲精品国产精品久久久不卡| 首页视频小说图片口味搜索| 久久国产精品影院| 日韩三级视频一区二区三区| 亚洲人成伊人成综合网2020| 国产乱人伦免费视频| 捣出白浆h1v1| 国产亚洲欧美精品永久| 一进一出抽搐gif免费好疼 | 国产成人精品在线电影| 黑人猛操日本美女一级片| www.自偷自拍.com| 国产精品 国内视频| 真人做人爱边吃奶动态| 国产野战对白在线观看| 国产精品 国内视频| 欧美大码av| 老汉色∧v一级毛片| 国产亚洲欧美在线一区二区| 操出白浆在线播放| 热re99久久国产66热| 精品人妻1区二区| 香蕉丝袜av| 丰满人妻熟妇乱又伦精品不卡| 国产无遮挡羞羞视频在线观看| 国产精品一区二区在线观看99| 国产精品秋霞免费鲁丝片| 看片在线看免费视频| 亚洲专区中文字幕在线| 99在线人妻在线中文字幕 | 美国免费a级毛片| videos熟女内射| 久久精品国产亚洲av高清一级| 欧美黄色淫秽网站| 激情在线观看视频在线高清 | 欧美人与性动交α欧美精品济南到| 精品少妇久久久久久888优播| 自线自在国产av| 51午夜福利影视在线观看| 成人18禁在线播放| 欧美亚洲日本最大视频资源| 视频区欧美日本亚洲| 国产单亲对白刺激| 一边摸一边做爽爽视频免费| 岛国毛片在线播放| 久久精品国产99精品国产亚洲性色 | 我的亚洲天堂| 久久久久久免费高清国产稀缺| 一二三四社区在线视频社区8| 波多野结衣一区麻豆| 脱女人内裤的视频| 精品亚洲成国产av| 国产成人精品久久二区二区91| 老司机福利观看| 色婷婷久久久亚洲欧美| 久久亚洲精品不卡| а√天堂www在线а√下载 | 国产成人欧美| 一本一本久久a久久精品综合妖精| 在线观看午夜福利视频| 国产人伦9x9x在线观看| 不卡一级毛片| 久久人妻福利社区极品人妻图片| 欧美日韩亚洲国产一区二区在线观看 | 50天的宝宝边吃奶边哭怎么回事| 日韩免费av在线播放| 亚洲国产精品一区二区三区在线| 制服人妻中文乱码| 国产精品影院久久| 电影成人av| 这个男人来自地球电影免费观看| 国产真人三级小视频在线观看| 最新的欧美精品一区二区| 淫妇啪啪啪对白视频| 啦啦啦免费观看视频1| tocl精华| 国产免费男女视频| 亚洲,欧美精品.| 亚洲欧美激情综合另类| 精品福利观看| 国产免费现黄频在线看| av网站在线播放免费| 热re99久久精品国产66热6| 亚洲av电影在线进入| 多毛熟女@视频| 久久国产亚洲av麻豆专区| 午夜福利乱码中文字幕| 极品人妻少妇av视频| 身体一侧抽搐| 人成视频在线观看免费观看| 在线国产一区二区在线| 99精品欧美一区二区三区四区| 亚洲免费av在线视频| 国产精品免费视频内射| 国产欧美日韩精品亚洲av| 久久久国产成人免费| 色老头精品视频在线观看| 男女床上黄色一级片免费看| 夜夜夜夜夜久久久久| 中文亚洲av片在线观看爽 | 日本撒尿小便嘘嘘汇集6| 变态另类成人亚洲欧美熟女 | 色综合欧美亚洲国产小说| 黄色成人免费大全| 成年人免费黄色播放视频| 国产成人免费观看mmmm| 一二三四在线观看免费中文在| 午夜福利免费观看在线| 国产激情久久老熟女| 男人舔女人的私密视频| 色精品久久人妻99蜜桃| 色综合婷婷激情| 热99久久久久精品小说推荐| 美女 人体艺术 gogo| 亚洲色图av天堂| 国产欧美亚洲国产| av网站在线播放免费| 每晚都被弄得嗷嗷叫到高潮| 亚洲黑人精品在线| 国产激情欧美一区二区| 国产淫语在线视频| 精品无人区乱码1区二区| 久久久精品区二区三区| 69av精品久久久久久| 精品国产国语对白av| 啦啦啦 在线观看视频| 亚洲avbb在线观看| 国产深夜福利视频在线观看| 天天操日日干夜夜撸| 丁香六月欧美| 久久精品国产亚洲av香蕉五月 | 色综合婷婷激情| 欧美激情 高清一区二区三区| 十分钟在线观看高清视频www| 麻豆成人av在线观看| 日韩欧美免费精品| svipshipincom国产片| 日韩欧美国产一区二区入口| 18在线观看网站| 国产精品久久久久成人av| 中文字幕人妻丝袜一区二区| 老司机午夜福利在线观看视频| 国产亚洲一区二区精品| 国内毛片毛片毛片毛片毛片| 国精品久久久久久国模美| 久久亚洲精品不卡| 两人在一起打扑克的视频| 亚洲精品国产区一区二| 国产精品 国内视频| 日韩免费高清中文字幕av| 十八禁网站免费在线| svipshipincom国产片| 国产野战对白在线观看| 精品熟女少妇八av免费久了| 侵犯人妻中文字幕一二三四区| 久久中文字幕一级| 黄色毛片三级朝国网站| 他把我摸到了高潮在线观看| 精品欧美一区二区三区在线| 国产欧美日韩一区二区三区在线| 亚洲精品自拍成人| 欧美不卡视频在线免费观看 | 日本撒尿小便嘘嘘汇集6| 国产精品1区2区在线观看. | 老鸭窝网址在线观看| 中文字幕高清在线视频| 深夜精品福利| 亚洲精品国产一区二区精华液| 午夜精品国产一区二区电影| 久久天躁狠狠躁夜夜2o2o| 搡老熟女国产l中国老女人| 日本黄色日本黄色录像| 91麻豆av在线| 亚洲精品在线观看二区| 美女 人体艺术 gogo| 午夜福利视频在线观看免费| aaaaa片日本免费| 国产日韩欧美亚洲二区| 久久精品国产99精品国产亚洲性色 | 国产亚洲精品久久久久久毛片 | 久久精品国产综合久久久| 久久久久精品人妻al黑| 国产日韩一区二区三区精品不卡| 国产淫语在线视频| 亚洲成国产人片在线观看| 亚洲 国产 在线| 两性夫妻黄色片| 建设人人有责人人尽责人人享有的| 精品福利永久在线观看| 国产在视频线精品| 99精国产麻豆久久婷婷| 国产精品二区激情视频| 欧美日韩视频精品一区| 日本黄色视频三级网站网址 | 亚洲欧美一区二区三区久久| 美女高潮到喷水免费观看| 久久青草综合色| 国产免费av片在线观看野外av| 狂野欧美激情性xxxx| 久久久久久久国产电影| 青草久久国产| 欧美激情久久久久久爽电影 | 欧美大码av| 黄色女人牲交| 欧美中文综合在线视频| 波多野结衣av一区二区av| 女人精品久久久久毛片| 天天躁夜夜躁狠狠躁躁| 波多野结衣一区麻豆| 亚洲专区国产一区二区| 一级,二级,三级黄色视频| 亚洲午夜精品一区,二区,三区| 俄罗斯特黄特色一大片| 热re99久久国产66热| 女人久久www免费人成看片| 亚洲在线自拍视频| 国产蜜桃级精品一区二区三区 | 女同久久另类99精品国产91| 国产精品久久久久久人妻精品电影| 欧美激情久久久久久爽电影 | 国精品久久久久久国模美| 亚洲精品美女久久久久99蜜臀| 黄片播放在线免费| 99精品欧美一区二区三区四区| xxx96com| 99精品欧美一区二区三区四区| 国产精品一区二区免费欧美| 宅男免费午夜| 国产aⅴ精品一区二区三区波| 精品国产美女av久久久久小说| 国产高清videossex| 国产欧美日韩一区二区三| 国产片内射在线| 18禁国产床啪视频网站| 精品久久久久久久久久免费视频 | 亚洲国产精品合色在线| 久久人妻熟女aⅴ| 欧美人与性动交α欧美软件| 亚洲av成人av| 亚洲在线自拍视频| 99精品欧美一区二区三区四区| 国产高清videossex| 精品国产美女av久久久久小说| 亚洲精品国产色婷婷电影| 热99久久久久精品小说推荐| 色尼玛亚洲综合影院|