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

    運動相關(guān)腦電信號的運動意圖預(yù)測方法研究

    2018-05-23 00:45:48,,
    計算機(jī)測量與控制 2018年5期
    關(guān)鍵詞:電信號識別率特征向量

    ,,

    (中國電子科技集團(tuán)公司第二十七研究所,鄭州 450047 )

    0 引言

    運動意圖的神經(jīng)解碼已被廣泛研究并應(yīng)用于大腦運動皮質(zhì)(M1)的神經(jīng)元活動,其神經(jīng)元放電模式編碼肢體運動的方向信息的方式可用于目標(biāo)動作預(yù)測[1]。研究表明,在運動規(guī)劃和執(zhí)行過程中,后頂葉皮層區(qū)(posterior parietal cortex, PPC) 的運動相關(guān)電位扮演著重要的作用。

    運動相關(guān)電位是一種慢性皮質(zhì)電位,近年來運動意圖檢測將其被廣泛用于驅(qū)動腦-機(jī)接口系統(tǒng)外部設(shè)備。文獻(xiàn)[2]通過研究腦電圖的時間動態(tài)特性,探討了運動意圖的基本原理。文獻(xiàn)[3]通過對九名受試者EEG中運動相關(guān)電位的實時檢測,可提高BCI系統(tǒng)的精度和實時性。文獻(xiàn)[4-6]對健康受試者和中風(fēng)患者的一系列研究中,已經(jīng)確定可以通過運動相關(guān)電位檢測到運動意圖,其具有足夠短的延遲以通過觸發(fā)外周神經(jīng)引起皮質(zhì)快速興奮。在人類腦電信號解析方面,運動方向解碼雖然取得了一定的研究成果及突破。對于頭皮腦電信號(electroencephalographic, EEG)的方向識別準(zhǔn)確率比較低,一般在65%以下[7-8],文獻(xiàn)[9]利用共空間模式和功率譜密度算法識別左右腳運動意圖監(jiān)測,并獲得72.6%和72.2%的分類識別率。如何提取并優(yōu)化提取EEG信號中的有效成分,用于運動方向預(yù)測以及如何利用運動準(zhǔn)備前(即運動規(guī)劃階段)的EEG信號提高運動方向預(yù)測的正確率是目前國內(nèi)外腦電信號研究和運動方向預(yù)測的瓶頸之一。

    想象運動和實際運動之間具有類似的神經(jīng)機(jī)制,并且能夠通過想象運動相關(guān)的腦電信號識別參與想象運動的肢體類型[10],且運動前的神經(jīng)活動或者伴隨運動的神經(jīng)活動編碼了運動的方向、速度和其他信息[11]。事件相關(guān)電位(event-related potentials, ERPs)在一定程度上反映感覺、知覺或認(rèn)知活動動態(tài)加工的過程,其中運動相關(guān)電位(movement-related potentials, MRPs)作為ERPs的一種,可以反映運動準(zhǔn)備、運動過程和運動結(jié)束的動態(tài)加工過程[12];mu/beta節(jié)律的事件相關(guān)同步/去同步(event-related (de)synchronization, ERS/ERD),用來反映運動前、運動過程中和運動結(jié)束后感覺運動區(qū)能量的變化關(guān)系,其對側(cè)化特征也可以反映參與運動或想象運動的肢體類型[13]。

    為研究運動相關(guān)腦電如何表征運動意圖預(yù)測,本研究通過左、右手兩任務(wù)運動想象實驗范式,采集運動想象任務(wù)開始前運動準(zhǔn)備階段的腦電信號,以及運動想象任務(wù)執(zhí)行時的腦電信號數(shù)據(jù),通過聯(lián)合MRPs與mu/beta節(jié)律的ERD/ERS特征,采用WPD-CSP特征提取方法,利用支持向量機(jī)(support vector machine, SVM)進(jìn)行運動意圖預(yù)測解碼,根據(jù)運動意圖預(yù)測準(zhǔn)確率驗證提取運動意圖特征的有效性。

    1 實驗

    11名健康受試者參與實驗,其中8男3女,平均年齡26.4歲,碩士以上學(xué)歷,工科背景。所有受試者皆為右利手,視覺與聽覺均正常。除了受試者S1、S5、S6和S8參與過腦-機(jī)器接口EEG數(shù)據(jù)采集外,其他均沒有EEG和BCI實驗的經(jīng)驗,也沒有受試者有已知的感覺運動疾病或心理病史。

    1.1 范式設(shè)計

    實驗開始前,受試者按自己舒適的習(xí)慣方式坐在椅子上,處于放松狀態(tài),準(zhǔn)備開始兩種任務(wù)的運動想象實驗。每種任務(wù)包含30個trial,每個trial包含3個階段:運動準(zhǔn)備階段、運動規(guī)劃及運動執(zhí)行階段。

    圖1 完整過程的刺激范式示意圖

    一次完整的實驗過程如圖1所示,運動準(zhǔn)備階段(起始時刻t=0 ms):刺激范式顯示“+”,提示受試者集中注意力,準(zhǔn)備實驗;運動規(guī)劃階段(t=1 000 ms):單次運動想象任務(wù)實驗中,左、右兩個任務(wù)(“┤”、“├”) 隨機(jī)出現(xiàn)一次,持續(xù)時間為1 000 ms,提示受試者準(zhǔn)備進(jìn)行運動想象任務(wù),受試者身體保持靜止?fàn)顟B(tài);運動執(zhí)行階段(t=2 000 ms):范式出現(xiàn)“←”或“→”,提示受試者開始想象運動,持續(xù)時間為3 000 ms。每組實驗連續(xù)采集,實驗期間加入短暫的休息提示,受試者可適當(dāng)休息。

    1.2 數(shù)據(jù)采集

    本研究的數(shù)據(jù)采集設(shè)備為本研究機(jī)構(gòu)自主研制的16導(dǎo)EEG放大器[14],0~100 Hz帶寬,1 000 Hz采樣頻率,24位A/D轉(zhuǎn)換器。

    運動準(zhǔn)備過程會誘發(fā)大腦額葉皮質(zhì)區(qū)和初級運動皮層區(qū)的電位變化,不同的受試者,關(guān)聯(lián)通道會存在差異[15]。在本研究中采用國際標(biāo)準(zhǔn)的10-20系統(tǒng),主要采集運動功能區(qū)的9個電極FC3、FCz、FC4、C3、Cz、C4、CP3、CPz 和CP4,如圖2所示。

    圖2 電極位置分布圖

    2 信號處理及分類

    算法和實驗在Windows XP系統(tǒng)Matlab R2012b環(huán)境下實現(xiàn),配置環(huán)境為:Intel(R) Core(TM) i7-5500U CPU @2.4 GHz。

    本研究的信號處理及分類算法設(shè)計由四部分組成:數(shù)據(jù)采集、預(yù)處理、特征提取和分類識別。算法設(shè)計流程圖如圖3所示。

    圖3 算法設(shè)計流程圖

    2.1 數(shù)據(jù)預(yù)處理

    數(shù)據(jù)預(yù)處理包括:降采樣、陷波、去基線漂移、基線校準(zhǔn)、去眼電。具體如表1所示。

    表1 腦電信號預(yù)處理算法

    2.2 特征提取及分類識別

    2.2.1 特征提取

    本研究提出一種聯(lián)合頻域和空域特征的多維特征提取算法,首先對待分析腦電信號進(jìn)行小波包分解[16](wavelet packet decomposition,WPD),再重構(gòu);然后采用共空間模式 (common spatial patterns, CSP)進(jìn)行空域濾波特征提取[17],實現(xiàn)時-頻-空域的多維聯(lián)合。WPD是一種高效的時頻分析快速算法,通過多尺度小波分解,構(gòu)建特征頻段的小波包分解系數(shù)特征向量;CSP提取不同電極所采集的運動想象信號的高維空域特征進(jìn)行對角化,分別最大限度的最大化和最小化兩類方差。聯(lián)合離散小波分解和共空間模式不僅能有效利用WPD的時頻特征,而且能彌補(bǔ)CSP缺乏頻域信息等缺陷。

    WPD不僅將頻帶進(jìn)行多層次劃分,而且對高頻部分也進(jìn)一步的分解,具有更好的時頻特性,提高了運動想象EEG信號的分析精度。以db4小波對每個導(dǎo)聯(lián)的信號進(jìn)行6層分解,得到S(6,0)-S(6,63)共64個子帶信號,其中S(6,0)對應(yīng)的子帶頻率為[0 7.8]Hz,在EEG信號中,有關(guān)方向信息主要集中在低頻和高頻部分,由于EEG信號通過顱骨傳輸?shù)筋^皮時出現(xiàn)高頻的劇烈衰減,主要分析低頻段數(shù)據(jù),故選取S(6,0)子帶信號。

    經(jīng)WPD重構(gòu)后得到一個具有頻率特征的矩陣EH×L,其中H是通道數(shù),L是采樣點數(shù),歸一化后的腦電數(shù)據(jù)的協(xié)方差矩陣為:

    (1)

    使用C1和C2分別表示兩種狀態(tài)下的協(xié)方差矩陣,則合成的協(xié)方差矩陣Ct:

    Ct=C1+C2

    (2)

    又可以將Ct寫成Ct=UtλtUt,為矩陣Ct的特征向量,λt是相對應(yīng)的特征值,通過該變化,特征值按照降序進(jìn)行了排列,與其對應(yīng)的特征向量也重新進(jìn)行了排列,經(jīng)過白化變化得到:

    (3)

    則C1和C2的協(xié)方差矩陣變換為:

    S1=PC1PT,S2=PC2PT

    (4)

    S1與S2有相同的特征向量,也就是如果S1=Bλ1BT,S2=Bλ2BT,并且λ1+λ2=I,S1與S2有相同的特征向量B,I是單位矩陣。由于兩類矩陣特征值之和為I,則S1的最大特征值所對應(yīng)的特征向量使S2有最小的特征值。

    將白化后的EEG信號投影在特征向量B的最小和最大列特征向量上,能夠得到最佳的分類特征。投影矩陣為:

    W=BTP

    (5)

    腦電信號最終可轉(zhuǎn)換為:

    Z=WE

    (6)

    CSP投影矩陣用于訓(xùn)練分類器,只有少數(shù)信號m可以有效區(qū)分各類別。能夠最大區(qū)分的信號Zp通常為λ1和λ2的最大值的對角化矩陣,其為Z的第一和最后m行。特征向量由下列公式獲得:

    (7)

    在本研究中m=3。

    2.2.2 分類識別

    對于單次實驗運動意圖腦電信號,具有高度的非線性和非平穩(wěn)性,可供學(xué)習(xí)的樣本較少且與電極組合后樣本特征向量維度較高,降維則會損失有用的信息。SVM[18]利用非線性核函數(shù)將輸入空間的特征向量變換到一個高維空間,使樣本線性可分,在高維空間構(gòu)造線性判別函數(shù),建立一個超平面作為決策曲面,使得正例和反例之間的隔離邊緣被最大化,把樣本空間映射到高維空間,把低維空間的線性不可分問題映射為高維空間的線性可分問題。

    研究中采用徑向基函數(shù)內(nèi)核:K(xi,xj)=exp(-g(xi-xj))2,作為SVM的核函數(shù),核函數(shù)系數(shù)g>0,懲罰因子C則決定了SVM的收斂速度及推廣能力,是影響SVM性能的另一個重要參數(shù)。在SVM參數(shù)選取時采用網(wǎng)格法,首先設(shè)定C和g取值范圍,然后進(jìn)行網(wǎng)格劃分,遍歷網(wǎng)格內(nèi)的所有取值點,計算單次遍歷的C和g參數(shù)下訓(xùn)練集的分類正確率,同時利用五折交叉驗證方法保證結(jié)果有效性,依次計算各次遍歷點選出分類正確率最高時對應(yīng)的C和g的值即可。

    圖4 參數(shù)C和g選取等高線

    3 結(jié)果和討論

    3.1 導(dǎo)聯(lián)通道優(yōu)化

    采集多名受試者EEG信號,疊加平均EEG中的ERP成分,繪制出一個Trail的腦地形圖,分析時序關(guān)系和各個時刻激活的相關(guān)區(qū)域,確定與運動準(zhǔn)備相關(guān)的通道。其中受試者S8的腦地形圖序列如圖5所示,由圖5所示,運動的起始時刻為2 000 ms,從運動起始前的1 000 ms開始繪制腦地形圖,每隔100 ms畫一幅,觀察大腦的激活腦區(qū),在1 000~1 100 ms期間,大腦處于靜息的狀態(tài),從1 400 ms開始,可觀察到大腦激活腦區(qū)開始有明顯變化。在運動開始時刻前的220~280 ms,內(nèi)側(cè)前額及額中央部觀察到ERP的差異,這些活動反映運動意圖的預(yù)測過程。實驗結(jié)果表明,本研究選取的FC3、FCz、FC4、C3、Cz、C4、CP3、CPz 和CP4通道所覆蓋的后頂葉皮層區(qū)能夠表征有關(guān)運動方向的信息。

    如圖5所示,運動開始前至少500 ms(約1 500 ms)出現(xiàn)緩慢負(fù)向變化的運動準(zhǔn)備電位可能反映了運動規(guī)劃,運動想象開始后出現(xiàn)由負(fù)電位向正電位變化可能反映了對運動的精細(xì)控制,結(jié)果進(jìn)一步證實了MRPs可以表征運動準(zhǔn)備、運動規(guī)劃和運動執(zhí)行的腦神經(jīng)機(jī)制。

    圖5 不同時刻腦區(qū)激活的序列圖

    3.2 MRPs及mu和beta節(jié)律特征分析

    針對右手想象運動腦電信號,分別計算C3和C4通道的mu/beta節(jié)律的功率譜變化, S8受試者的C3和C4通道腦電信號功率譜見圖6。

    圖6 mu/beta節(jié)律的功率譜

    由圖6結(jié)果可知,對于C3通道,右手想象時mu和beta節(jié)律的功率譜要比左手小,而對于C4電通道結(jié)果正好相反,這就說明ERD也具有對側(cè)半球優(yōu)勢。實驗結(jié)果發(fā)現(xiàn),在動作發(fā)生前10 Hz頻率附近的mu/beta節(jié)律會有能量遞減變化。

    為了進(jìn)一步細(xì)化研究運動想象信號中表征運動意圖的頻段,將S8受試者的C3和C4導(dǎo)聯(lián)的MRPs信號進(jìn)行了1~10 Hz的時-頻變換,時-頻圖譜如圖7所示。

    圖7 C3通道和C4通道MRPs時-頻圖譜

    圖7的時-頻圖譜研究表明,運動規(guī)劃和執(zhí)行階段EEG信號中的運動意圖信息主要集中在低頻段,通過對11名受試者腦電信號時頻圖譜的比較分析,發(fā)現(xiàn)在1~7 Hz的信號中包含運動方向的信息。

    3.3 分類結(jié)果

    對左、右手想象數(shù)據(jù)進(jìn)行預(yù)處理,采用200 ms的滑動時間窗進(jìn)行截取信號,對兩類數(shù)據(jù)按照公式(1)到(7)求協(xié)方差,求白化矩陣和空間源成分,尋找特征值最大和最小的空間源成分對信號進(jìn)行空間濾波,提取低頻分量MRPs與mu和beta節(jié)律的有效的特征。本研究中每位受試者分別進(jìn)行十組實驗,采用SVM算法進(jìn)行分類識別,各次結(jié)果均為十折交叉檢驗后的結(jié)果,分類識別率結(jié)果如圖8和表2所示。

    圖8 六名受試者分類識別率結(jié)果

    表2 11名受試者平均分類識別結(jié)果

    由圖8結(jié)果可知,受試者S3對運動想象不敏感,分類識別率結(jié)果接近隨機(jī),出現(xiàn)運動想象“盲”現(xiàn)象,其余受試者由于個體差異性,分類識別率略有差異但基本穩(wěn)定。根據(jù)圖8和表2結(jié)果可知,無腦-機(jī)接口實驗經(jīng)驗受試者分類識別率在70%左右,四名有運動想象經(jīng)驗的受試者平均分類識別率為74%左右,最高可達(dá)到85%。本實驗結(jié)果與文獻(xiàn)[19]中Lakany和 Conway的兩任務(wù)分類識別率相當(dāng),高于文獻(xiàn)[7-8]中的識別率。本文所提算法獲得分類結(jié)果高于文獻(xiàn)[20]中采用ICA和SVM算法所得結(jié)果,驗證了本實驗的可行性和有效性,同時也為下一步確定運動意圖、解析動作狀態(tài)提供參考。

    4 結(jié)論

    本研究以手部運動想象腦電信號為研究對象,以融合MRPs與ERS/ERD特征為途徑,設(shè)計了包含時頻域與空域的特征提取算法,實現(xiàn)了手部運動意圖預(yù)測,且識別正確率有了顯著提高,經(jīng)過運動想象訓(xùn)練的受試者最高分類識別率達(dá)到85%。初步驗證了手部運動想象腦電信號能夠解碼運動過相關(guān)電位信息,且集中在10 Hz以下的低頻信號成分中,為下一步的運動精細(xì)控制提供技術(shù)支持。

    然而仍有一些研究需要進(jìn)一步明確:1)能否設(shè)計自適應(yīng)的算法進(jìn)一步提升運動預(yù)測的準(zhǔn)確率;2)結(jié)合運動準(zhǔn)備電位以及運動執(zhí)行電位是否可以提高識別準(zhǔn)確率;3)腦電信號能否解碼想象運動的運動速度、力的變化率和目標(biāo)力大小、運動軌跡等或參數(shù)。

    參考文獻(xiàn):

    [1]Hochberg L R, Serruya M D, Friehs G M, et al. Neuronal ensemble control of prosthetic devices by a human with tetraplegia[J]. Nature, 2006, 442(7099): 164-171.

    [2]Wairagkar M, Hayashi Y, Nasuto S. Movement Intention Detection from Autocorrelation of EEG for BCI[M]. Brain Informatics and Health. Springer International Publishing, 2015: 212-221.

    [3]Xu R, Jiang N, Lin C, et al. Enhanced Low-Latency Detection of Motor Intention From EEG for Closed-Loop Brain-Computer Interface Applications [J]. IEEE Trans Biomed Eng, 2014, 61 (2) :288-296.

    [4]Niazi I K, Jiang N, Tiberghien O, et al. Detection of movement intention from single-trial movement-related cortical potentials[J]. Journal of Neural Engineering , 2011 , 8 (6) :066009.

    [5]Niazi I K, Mrachacz-Kersting N, Jiang N, et al. Peripheral electrical stimulation triggered by self-paced detection of motor intention enhances motor evoked potentials[J]. IEEE Trans Neural Syst Rehabil Eng, 2012, 20 (4) :595-604.

    [6]Mrachacz-Kersting N, Jiang N, Dremstrup K, et al. A Novel Brain-Computer Interface for Chronic Stroke Patients[M]. Brain-Computer Interface Research. Springer Berlin Heidelberg, 2014: 51-61.

    [7] Waldert S, Preissl H, Demandt E, et al. Hand movement direction decoded from MEG and EEG [J]. Journal of Neuroscience, 2008, 28(4):1000-1008.

    [8]伏云發(fā), 徐保磊, 裴立力,等. 基于腦電運動速度想象的單次識別研究[J]. 中國生物醫(yī)學(xué)工程學(xué)報, 2011,30(4): 555-561.

    [9]Kee Y J, Won D O, Lee S W. Classification of left and right foot movement intention based on steady-state somatosensory evoked potentials[A].International Winter Conference on Brain-Computer Interface[C]. 2017:106-108.

    [10]Pfurtscheller G, Neuper C. Motor imagery and direct brain-computer communication[J]. Proceedings of the IEEE, 2002, 89(7):1123-1134.

    [11]Gu Y, Dremstrup K, Farina D. Single-trial discrimination of type and speed of wrist movements from EEG recordings [J]. Clinical Neurophysiology, 2009, 120(8): 1596-1600.

    [12]OF do Nascimento, KD Nielsen, M Voigt. Movement related parameters modulate cortical activity during imaginary isometric plantar-flexions[J]. Experimental Brain Research, 2006 , 171 (1) :78-90.

    [13]Neuper C, Wortz M, Pfurtscheller G. ERD/ERS patterns reflecting sensorimotor activation and deactivation [J]. Progress in Brain Research , 2006 , 159 (1) :211-222.

    [14]孫廣金. 基于ADS1299的新型腦電采集系統(tǒng)設(shè)計[J]. 科技信息,2014 (1) :73-74.

    [15]Takahashi J, Yasunaga D. Reported visual imagery and apparent motion: an event-related potential study[J]. Neuroreport, 2012, 23(15):904-910.

    [16]Lin E B, Shen X P. Wavelet analysis of EEG signals. Aerospace and Electronics[A].IEEE National Conference - NAECON[C]. 2011, 105-110.

    [17]孟 明,朱俊青,佘青山,等. 多類運動想象腦電信號的兩級特征提取方法[J]. 自動化學(xué)報, 2016 , 42 (12) :1915-1922.

    [18] Müller KR, Krauledat M, Dornhege G, et al. Machine learning techniques for brain-computer interfaces [J]. Biomedical Engineering , 2004 :11-22.

    [19]Sporns O. The human connectome: a complex network [J]. Ann N Y Acad Sci , 2011 , 1224 (1) :109-125.

    [20]He Y, Chen ZJ, Evans AC. Small-world anatomical networks in the human brain revealed by cortical thickness from MRI [J]. Cerebral Cortex, 2007, 17 (10) :2407-2419.

    猜你喜歡
    電信號識別率特征向量
    二年制職教本科線性代數(shù)課程的幾何化教學(xué)設(shè)計——以特征值和特征向量為例
    克羅內(nèi)克積的特征向量
    基于聯(lián)合聚類分析的單通道腹部心電信號的胎心率提取
    基于類圖像處理與向量化的大數(shù)據(jù)腳本攻擊智能檢測
    基于真耳分析的助聽器配戴者言語可懂度指數(shù)與言語識別率的關(guān)系
    基于Code Composer Studio3.3完成對心電信號的去噪
    科技傳播(2019年24期)2019-06-15 09:29:28
    提升高速公路MTC二次抓拍車牌識別率方案研究
    一類特殊矩陣特征向量的求法
    基于隨機(jī)森林的航天器電信號多分類識別方法
    EXCEL表格計算判斷矩陣近似特征向量在AHP法檢驗上的應(yīng)用
    色播亚洲综合网| 波多野结衣高清无吗| 亚洲美女黄片视频| 乱人视频在线观看| 久久久久性生活片| 99热这里只有是精品50| 小说图片视频综合网站| 亚洲成人中文字幕在线播放| 18禁在线播放成人免费| 国产视频内射| 色综合婷婷激情| av在线蜜桃| 男人舔女人下体高潮全视频| 日韩强制内射视频| 国产中年淑女户外野战色| 美女被艹到高潮喷水动态| 全区人妻精品视频| 深夜精品福利| 美女cb高潮喷水在线观看| 极品教师在线视频| 亚洲欧美日韩东京热| 午夜福利在线在线| 有码 亚洲区| 给我免费播放毛片高清在线观看| 级片在线观看| 黄色欧美视频在线观看| 综合色av麻豆| 午夜福利在线观看免费完整高清在 | 中文在线观看免费www的网站| 99热精品在线国产| 亚洲,欧美,日韩| h日本视频在线播放| 亚洲电影在线观看av| 亚洲精品亚洲一区二区| 亚洲自拍偷在线| 九九久久精品国产亚洲av麻豆| 国产亚洲91精品色在线| 俄罗斯特黄特色一大片| 嫩草影院入口| 人妻少妇偷人精品九色| 色av中文字幕| av在线亚洲专区| 国产色爽女视频免费观看| 日韩大尺度精品在线看网址| 两个人的视频大全免费| 综合色av麻豆| 99久国产av精品| 亚洲,欧美,日韩| 久久欧美精品欧美久久欧美| 黄色一级大片看看| 少妇熟女aⅴ在线视频| 亚洲成人免费电影在线观看| 日韩中文字幕欧美一区二区| 欧美区成人在线视频| 精品福利观看| av.在线天堂| 成人国产一区最新在线观看| 精品人妻视频免费看| 高清日韩中文字幕在线| 久久久精品欧美日韩精品| 国产成人影院久久av| 麻豆精品久久久久久蜜桃| 最近最新中文字幕大全电影3| 午夜视频国产福利| 欧美激情久久久久久爽电影| ponron亚洲| 黄片wwwwww| 亚洲18禁久久av| 69人妻影院| 男人的好看免费观看在线视频| 日本三级黄在线观看| 看十八女毛片水多多多| 天堂动漫精品| 亚洲成人久久爱视频| 成人欧美大片| 日本与韩国留学比较| 国产一区二区三区在线臀色熟女| 亚洲精品456在线播放app | 欧美成人a在线观看| 九九爱精品视频在线观看| 亚洲专区国产一区二区| 99久久精品热视频| 简卡轻食公司| 韩国av一区二区三区四区| 国产av在哪里看| 听说在线观看完整版免费高清| 日本撒尿小便嘘嘘汇集6| 在线播放无遮挡| 人人妻人人澡欧美一区二区| 色播亚洲综合网| 尤物成人国产欧美一区二区三区| 老熟妇仑乱视频hdxx| 床上黄色一级片| 国产中年淑女户外野战色| av在线亚洲专区| 日韩中字成人| 亚洲欧美清纯卡通| 欧美一级a爱片免费观看看| 色哟哟哟哟哟哟| 变态另类丝袜制服| 欧美丝袜亚洲另类 | 欧美xxxx黑人xx丫x性爽| 成人永久免费在线观看视频| 97热精品久久久久久| 九九爱精品视频在线观看| 九九在线视频观看精品| 亚洲国产欧洲综合997久久,| xxxwww97欧美| 女同久久另类99精品国产91| 久久精品影院6| 亚洲成人久久爱视频| 韩国av一区二区三区四区| 亚洲精品成人久久久久久| 女人十人毛片免费观看3o分钟| 欧洲精品卡2卡3卡4卡5卡区| 两个人的视频大全免费| 狠狠狠狠99中文字幕| 亚洲人成网站在线播| 欧美最新免费一区二区三区| 有码 亚洲区| 桃色一区二区三区在线观看| 久久香蕉精品热| 国产高清视频在线观看网站| 蜜桃亚洲精品一区二区三区| av天堂中文字幕网| 亚洲国产欧美人成| 欧美成人免费av一区二区三区| 免费无遮挡裸体视频| 日韩欧美一区二区三区在线观看| 免费观看的影片在线观看| 成人综合一区亚洲| 国产 一区精品| 亚洲国产精品成人综合色| 国产久久久一区二区三区| 久久香蕉精品热| 麻豆国产97在线/欧美| 很黄的视频免费| 成人永久免费在线观看视频| 亚洲18禁久久av| 午夜激情福利司机影院| 精品福利观看| 亚洲精品国产成人久久av| 性欧美人与动物交配| 国产精品国产高清国产av| 免费av毛片视频| 在现免费观看毛片| 两性午夜刺激爽爽歪歪视频在线观看| 在线免费观看的www视频| 午夜精品久久久久久毛片777| 国产激情偷乱视频一区二区| 日韩精品有码人妻一区| 美女cb高潮喷水在线观看| 免费黄网站久久成人精品| 国产精品人妻久久久影院| 免费在线观看影片大全网站| 神马国产精品三级电影在线观看| 免费看a级黄色片| 99久久精品热视频| 久久久成人免费电影| 国产精品日韩av在线免费观看| 久久欧美精品欧美久久欧美| 国产精品98久久久久久宅男小说| 国产毛片a区久久久久| 亚洲一区二区三区色噜噜| 3wmmmm亚洲av在线观看| 亚洲精品粉嫩美女一区| 欧美精品啪啪一区二区三区| 国产伦一二天堂av在线观看| 国产极品精品免费视频能看的| 日日撸夜夜添| 欧美成人a在线观看| 夜夜夜夜夜久久久久| 日日摸夜夜添夜夜添av毛片 | 久久亚洲精品不卡| 国产高潮美女av| 亚洲av五月六月丁香网| 伊人久久精品亚洲午夜| 可以在线观看毛片的网站| 精品日产1卡2卡| 真实男女啪啪啪动态图| 亚洲av二区三区四区| 12—13女人毛片做爰片一| 国产一区二区在线观看日韩| 床上黄色一级片| 给我免费播放毛片高清在线观看| 国产在视频线在精品| av天堂中文字幕网| 国产伦精品一区二区三区视频9| 91久久精品国产一区二区成人| 美女高潮的动态| 午夜免费男女啪啪视频观看 | 五月伊人婷婷丁香| 国内毛片毛片毛片毛片毛片| 51国产日韩欧美| 国产av麻豆久久久久久久| 欧美一区二区精品小视频在线| АⅤ资源中文在线天堂| 欧美成人免费av一区二区三区| 亚洲中文字幕一区二区三区有码在线看| 淫秽高清视频在线观看| 国产在线精品亚洲第一网站| 少妇丰满av| 麻豆国产av国片精品| 色5月婷婷丁香| 久久久成人免费电影| 狂野欧美白嫩少妇大欣赏| 网址你懂的国产日韩在线| 床上黄色一级片| 欧美另类亚洲清纯唯美| 亚洲精品亚洲一区二区| 男人和女人高潮做爰伦理| 国产成人一区二区在线| 亚洲熟妇中文字幕五十中出| 国产在线男女| 99国产极品粉嫩在线观看| 一边摸一边抽搐一进一小说| 日本熟妇午夜| 欧美性猛交黑人性爽| 精品久久久久久久久久久久久| 午夜免费成人在线视频| 一本一本综合久久| 一级黄色大片毛片| 99久久中文字幕三级久久日本| 亚洲欧美日韩卡通动漫| 久久久久国内视频| 国产一区二区三区在线臀色熟女| 少妇的逼好多水| 国产国拍精品亚洲av在线观看| 欧美成人a在线观看| 最近视频中文字幕2019在线8| 波多野结衣高清作品| 色播亚洲综合网| 成人午夜高清在线视频| 超碰av人人做人人爽久久| 美女xxoo啪啪120秒动态图| а√天堂www在线а√下载| 观看美女的网站| 舔av片在线| 岛国在线免费视频观看| 少妇猛男粗大的猛烈进出视频 | 国产精品爽爽va在线观看网站| 国产美女午夜福利| 欧美最黄视频在线播放免费| 又紧又爽又黄一区二区| 日本精品一区二区三区蜜桃| 亚洲国产色片| 又黄又爽又刺激的免费视频.| 伊人久久精品亚洲午夜| 人人妻人人澡欧美一区二区| 99在线人妻在线中文字幕| 一区二区三区四区激情视频 | netflix在线观看网站| 亚洲人成网站高清观看| 黄色欧美视频在线观看| 黄色欧美视频在线观看| 成人美女网站在线观看视频| 窝窝影院91人妻| 国产激情偷乱视频一区二区| 亚洲av成人精品一区久久| 久久人妻av系列| 在线观看66精品国产| 久久久久久久久中文| 国产精品久久久久久久电影| 国产男人的电影天堂91| 韩国av在线不卡| 国产精品亚洲美女久久久| 亚洲美女视频黄频| 久久久色成人| 中出人妻视频一区二区| 中文资源天堂在线| 欧美日韩黄片免| 丰满人妻一区二区三区视频av| 天天躁日日操中文字幕| 成人亚洲精品av一区二区| 亚洲av中文字字幕乱码综合| 日本五十路高清| 亚洲av一区综合| 女的被弄到高潮叫床怎么办 | 欧美成人性av电影在线观看| 亚洲久久久久久中文字幕| 欧美日韩综合久久久久久 | 乱系列少妇在线播放| 日韩在线高清观看一区二区三区 | 久久6这里有精品| 亚洲国产精品久久男人天堂| 人妻夜夜爽99麻豆av| 欧美一区二区亚洲| 午夜视频国产福利| 国产一区二区激情短视频| 最近最新中文字幕大全电影3| 午夜福利欧美成人| 女人被狂操c到高潮| 两个人的视频大全免费| 性色avwww在线观看| 少妇高潮的动态图| 日本爱情动作片www.在线观看 | 午夜福利18| 亚洲人与动物交配视频| 日韩人妻高清精品专区| 最近在线观看免费完整版| av天堂中文字幕网| 午夜福利在线观看吧| 亚洲最大成人中文| 日日干狠狠操夜夜爽| 久久久久久久久中文| 亚洲熟妇熟女久久| 国产精品一区二区免费欧美| 国内揄拍国产精品人妻在线| 内地一区二区视频在线| 3wmmmm亚洲av在线观看| 人妻久久中文字幕网| 婷婷精品国产亚洲av在线| 国产精品嫩草影院av在线观看 | 搡老妇女老女人老熟妇| 老司机深夜福利视频在线观看| 琪琪午夜伦伦电影理论片6080| 精品久久久噜噜| 日韩欧美一区二区三区在线观看| 直男gayav资源| 久久精品国产自在天天线| 韩国av在线不卡| 在线观看av片永久免费下载| 亚洲欧美日韩高清专用| 免费在线观看影片大全网站| 成人av在线播放网站| 不卡一级毛片| 麻豆成人午夜福利视频| 日本黄色视频三级网站网址| 性插视频无遮挡在线免费观看| 成人高潮视频无遮挡免费网站| 露出奶头的视频| 联通29元200g的流量卡| 国产一区二区激情短视频| 日韩大尺度精品在线看网址| 国产高清激情床上av| 日韩亚洲欧美综合| 亚洲av不卡在线观看| 久久九九热精品免费| 少妇的逼水好多| 久久人人精品亚洲av| 丝袜美腿在线中文| 大又大粗又爽又黄少妇毛片口| 国产午夜精品久久久久久一区二区三区 | 国产高清三级在线| 国产免费一级a男人的天堂| 91久久精品国产一区二区三区| 人妻久久中文字幕网| 成年女人看的毛片在线观看| 国产精品自产拍在线观看55亚洲| 老司机福利观看| 999久久久精品免费观看国产| 我要看日韩黄色一级片| av黄色大香蕉| 午夜福利在线在线| 亚洲va在线va天堂va国产| 免费av观看视频| 国产免费av片在线观看野外av| 日韩 亚洲 欧美在线| 别揉我奶头~嗯~啊~动态视频| 国产精品不卡视频一区二区| 国产精品久久久久久精品电影| 很黄的视频免费| 国产毛片a区久久久久| 久久久久久国产a免费观看| 99久久成人亚洲精品观看| 一a级毛片在线观看| 极品教师在线视频| 国产主播在线观看一区二区| 日韩欧美国产在线观看| 国产91精品成人一区二区三区| 欧美高清性xxxxhd video| 热99在线观看视频| 久久精品国产鲁丝片午夜精品 | 精品一区二区三区av网在线观看| 国产精品一及| 精品久久久久久久久久久久久| 亚洲综合色惰| 午夜精品一区二区三区免费看| 波多野结衣高清作品| 国产探花在线观看一区二区| 成人特级av手机在线观看| a级一级毛片免费在线观看| 久久久久久久午夜电影| 91狼人影院| av在线天堂中文字幕| 国产精品人妻久久久久久| 久久精品综合一区二区三区| 免费av观看视频| 俄罗斯特黄特色一大片| 欧美精品啪啪一区二区三区| 色精品久久人妻99蜜桃| 午夜福利在线观看免费完整高清在 | 久久久久久久久大av| 精品久久久久久久久久久久久| 我的女老师完整版在线观看| 久久久久久久亚洲中文字幕| 天堂√8在线中文| 日日摸夜夜添夜夜添小说| 国产精品三级大全| 亚洲性久久影院| 我要搜黄色片| 婷婷精品国产亚洲av| www.www免费av| 永久网站在线| 国产中年淑女户外野战色| 最好的美女福利视频网| 亚洲av电影不卡..在线观看| 最新中文字幕久久久久| 最近最新中文字幕大全电影3| 国产高清视频在线观看网站| 十八禁国产超污无遮挡网站| av国产免费在线观看| 久久天躁狠狠躁夜夜2o2o| 久久久久九九精品影院| 少妇高潮的动态图| 少妇裸体淫交视频免费看高清| 亚洲中文日韩欧美视频| 国产精品,欧美在线| 日本黄色片子视频| 午夜福利在线观看免费完整高清在 | 久久久久久久亚洲中文字幕| 黄色女人牲交| 国产精品伦人一区二区| 精品人妻1区二区| 国产淫片久久久久久久久| 日本熟妇午夜| 亚洲精品粉嫩美女一区| 蜜桃久久精品国产亚洲av| 久久久久久久久中文| 色综合亚洲欧美另类图片| 一级黄片播放器| 国产精品久久久久久久电影| 日本 欧美在线| 老熟妇乱子伦视频在线观看| 91在线精品国自产拍蜜月| 日本熟妇午夜| 久久久久久久精品吃奶| 亚洲国产精品合色在线| 99九九线精品视频在线观看视频| 日韩在线高清观看一区二区三区 | 最后的刺客免费高清国语| 国产精品人妻久久久久久| 国产av不卡久久| 精品午夜福利视频在线观看一区| 亚洲 国产 在线| 中国美女看黄片| 特级一级黄色大片| 亚洲国产欧美人成| 欧美一级a爱片免费观看看| 一进一出好大好爽视频| 欧美极品一区二区三区四区| 日本五十路高清| 一进一出好大好爽视频| 免费看日本二区| 日韩人妻高清精品专区| 精品欧美国产一区二区三| 国产熟女欧美一区二区| 午夜老司机福利剧场| 悠悠久久av| 色视频www国产| 日本 av在线| 亚洲在线观看片| 亚洲最大成人av| 亚洲无线观看免费| 九九热线精品视视频播放| 欧美成人免费av一区二区三区| 国产在视频线在精品| 人人妻人人澡欧美一区二区| 别揉我奶头~嗯~啊~动态视频| 国产高潮美女av| 无人区码免费观看不卡| 人人妻,人人澡人人爽秒播| 成人av在线播放网站| 成人毛片a级毛片在线播放| 一进一出抽搐动态| 人妻夜夜爽99麻豆av| 亚洲人与动物交配视频| 日韩精品青青久久久久久| 熟女人妻精品中文字幕| 久久热精品热| 麻豆成人午夜福利视频| 不卡一级毛片| 亚洲最大成人手机在线| 国产一区二区三区视频了| 别揉我奶头~嗯~啊~动态视频| 12—13女人毛片做爰片一| 99精品久久久久人妻精品| 91久久精品电影网| 国产精品一区www在线观看 | 欧美zozozo另类| 国产高清不卡午夜福利| 久久午夜福利片| 亚洲专区国产一区二区| 性色avwww在线观看| 九九在线视频观看精品| 国产亚洲精品久久久com| 亚洲人成网站在线播放欧美日韩| www.www免费av| 国产精品久久久久久久电影| 国产精华一区二区三区| 亚洲电影在线观看av| 亚洲欧美清纯卡通| av.在线天堂| 国产精品嫩草影院av在线观看 | 久久久久久国产a免费观看| 亚洲成av人片在线播放无| 男女之事视频高清在线观看| 一个人看视频在线观看www免费| 舔av片在线| 熟妇人妻久久中文字幕3abv| 18禁在线播放成人免费| 国产精品电影一区二区三区| 哪里可以看免费的av片| 99久久成人亚洲精品观看| 精品国产三级普通话版| 欧美3d第一页| av在线观看视频网站免费| 欧美xxxx性猛交bbbb| 国内少妇人妻偷人精品xxx网站| 国产欧美日韩一区二区精品| 熟女电影av网| 一边摸一边抽搐一进一小说| 特级一级黄色大片| 长腿黑丝高跟| 天堂影院成人在线观看| 国产精品,欧美在线| 日韩大尺度精品在线看网址| 最近中文字幕高清免费大全6 | 亚洲久久久久久中文字幕| 国产黄a三级三级三级人| 日本熟妇午夜| 99精品久久久久人妻精品| 少妇被粗大猛烈的视频| 国产精品电影一区二区三区| 日韩av在线大香蕉| 久久精品人妻少妇| 国产精品久久电影中文字幕| 色综合站精品国产| 三级男女做爰猛烈吃奶摸视频| www.www免费av| 久久精品91蜜桃| 午夜免费激情av| 九九在线视频观看精品| 男人的好看免费观看在线视频| 国产精品亚洲一级av第二区| 午夜福利视频1000在线观看| 中文亚洲av片在线观看爽| 午夜福利在线观看吧| 国产精品久久久久久久电影| 男人和女人高潮做爰伦理| 97人妻精品一区二区三区麻豆| 国产色婷婷99| 91在线精品国自产拍蜜月| 悠悠久久av| 亚洲男人的天堂狠狠| 有码 亚洲区| 国语自产精品视频在线第100页| 成人性生交大片免费视频hd| 婷婷六月久久综合丁香| 波多野结衣巨乳人妻| 欧美成人一区二区免费高清观看| 免费看a级黄色片| 99热这里只有是精品50| 国产高清有码在线观看视频| 1000部很黄的大片| 欧洲精品卡2卡3卡4卡5卡区| 噜噜噜噜噜久久久久久91| 在线a可以看的网站| 日本黄色片子视频| 亚洲国产欧洲综合997久久,| 色综合色国产| 一级黄片播放器| 97超视频在线观看视频| 欧美日韩瑟瑟在线播放| 国产成人av教育| 国产私拍福利视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看| av女优亚洲男人天堂| 男女做爰动态图高潮gif福利片| 日韩在线高清观看一区二区三区 | 日韩精品中文字幕看吧| 午夜爱爱视频在线播放| 久久99热这里只有精品18| 国产午夜精品论理片| 禁无遮挡网站| 一级a爱片免费观看的视频| 黄色欧美视频在线观看| 性插视频无遮挡在线免费观看| av中文乱码字幕在线| 亚洲欧美日韩高清在线视频| 久久6这里有精品| 久久欧美精品欧美久久欧美| 午夜久久久久精精品| 日本爱情动作片www.在线观看 | or卡值多少钱| 国产午夜福利久久久久久| 亚洲人成网站在线播| 舔av片在线| 三级国产精品欧美在线观看| 久久精品国产亚洲av天美| 亚洲狠狠婷婷综合久久图片| 国内精品美女久久久久久| 亚洲美女搞黄在线观看 | 国产单亲对白刺激| 成年女人永久免费观看视频| 在线天堂最新版资源| 国产精品伦人一区二区| 国国产精品蜜臀av免费| 国产老妇女一区| 如何舔出高潮| or卡值多少钱| av在线蜜桃| av黄色大香蕉| 岛国在线免费视频观看| 国产av在哪里看|