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

    小波包分析和FastICA 相結(jié)合對(duì)單通道腦電信號(hào)的去噪研究

    2021-04-08 01:55姚健康熊根良
    現(xiàn)代電子技術(shù) 2021年7期
    關(guān)鍵詞:單通道波包電信號(hào)

    姚健康,熊根良

    (南昌大學(xué) 機(jī)電工程學(xué)院 江西省機(jī)器人與焊接自動(dòng)化重點(diǎn)實(shí)驗(yàn)室,江西 南昌 330031)

    0 引 言

    隨著科學(xué)技術(shù)的發(fā)展,研究人員越來(lái)越重視對(duì)大腦的研究和開發(fā),腦電信號(hào)是大腦內(nèi)部神經(jīng)細(xì)胞在大腦皮層活動(dòng)的綜合體現(xiàn),是一種產(chǎn)生機(jī)理相當(dāng)復(fù)雜的隨機(jī)信號(hào)。通過(guò)電極采集獲取的腦電信號(hào)不但包含了神經(jīng)細(xì)胞的電活動(dòng)特征,還包含了大量的生理與病理信息[1]。但腦電信號(hào)十分微弱,極易受到噪聲的干擾,常見(jiàn)的偽跡有肌電偽跡和眼電偽跡[2]。因此有效地去除偽跡得到純凈的腦電信號(hào),成為了腦電信號(hào)處理分析中的關(guān)鍵,有利于信號(hào)的特征提取和分類。

    腦電信號(hào)去噪的主要方法有:平均偽跡回歸分析[3]、典型相關(guān)分析[4]、小波變換[5-7]、主成分分析[8]等,但這些傳統(tǒng)的去噪方法很難取得理想的效果。平均偽跡回歸分析需要采集獨(dú)立的腦電信號(hào);典型相關(guān)分析中無(wú)法有效解決腦電信號(hào)細(xì)節(jié)反映的問(wèn)題,可能會(huì)導(dǎo)致腦電信號(hào)失真;主成分分析中如果偽跡和腦電信號(hào)不相關(guān),那就不能真正地分離偽跡和腦電信號(hào);小波變換中計(jì)算過(guò)程復(fù)雜而且小波閾值的選定需要先驗(yàn)知識(shí)。

    快速獨(dú)立成分分析是一種有效去除腦電信號(hào)中噪聲的方法[9],能夠改善傳統(tǒng)ICA 算法的收斂速度和分離效果。但是FastICA 算法是以非高斯源信號(hào)為研究對(duì)象,通過(guò)利用分離信號(hào)中各自分量之間的最大獨(dú)立性進(jìn)行對(duì)比,因此FastICA 算法對(duì)獨(dú)立性的要求很高。

    由于不同腦區(qū)之間也存在著關(guān)聯(lián)性,會(huì)導(dǎo)致分離后的成分之間沒(méi)法保持絕對(duì)的獨(dú)立性。因此為更好地解決FastICA 算法對(duì)獨(dú)立性的要求,本文提出基于小波包分解的快速獨(dú)立成分分析對(duì)單通道腦電信號(hào)的去噪算法。首先選擇對(duì)單通道信號(hào)的處理可以避免多通道信號(hào)之間的關(guān)聯(lián)性對(duì)FastICA 算法獨(dú)立性要求的影響,其次利用小波包分解單通道信號(hào)得到頻帶較窄的子帶信號(hào),增強(qiáng)了源信號(hào)的獨(dú)立性且去除了不同腦電信號(hào)之間的統(tǒng)計(jì)相關(guān)性。對(duì)比FastICA 對(duì)多通道信號(hào)的去噪效果,本文提出的方法可以實(shí)現(xiàn)更好的去噪效果。

    1 算法的基本原理和改進(jìn)

    1.1 快速獨(dú)立成分分析

    FastICA 是一種基于負(fù)熵的獨(dú)立成分分析算法,其收斂速度快、分離性能好??赏ㄟ^(guò)觀測(cè)信號(hào)來(lái)預(yù)估個(gè)數(shù)未知的源信號(hào),從而恢復(fù)源信號(hào)。

    觀測(cè)信號(hào)X(t)={x1(t),x2(t),…,xn(t)}是源信號(hào)S(t)={s1(t),s2(t),…,sn(t)}通過(guò)矩陣A混合而成的,即為X=AS,ICA 就是在混合矩陣和源信號(hào)未知時(shí),求解混合矩陣W,從觀測(cè)信號(hào)中分離出源信號(hào)。

    FastICA 算法流程圖如圖1 所示。

    圖1 FastICA 的流程圖

    1.2 小波包分析

    1)在選擇的小波基進(jìn)行多分辨分析中可以將尺度子空間Oi和小波子空間Vi用一個(gè)新的子空間Rni來(lái)替代表示,令其中,Oi+1=Oi+Vi可以采用對(duì)的分解來(lái)代替表示:

    2)在小波包最優(yōu)基的選擇過(guò)程中,對(duì)長(zhǎng)度為N=2t的信號(hào)進(jìn)行分解過(guò)于復(fù)雜,對(duì)出現(xiàn)的每一類情況進(jìn)行逐次處理難以實(shí)現(xiàn)。因此選擇Dmeyer 小波,其具有良好的正交性和對(duì)稱性,可以在一定程度上減少分解和重構(gòu)中的相位失真。

    1.3 電參考和噪聲獨(dú)立成分去除

    1)經(jīng)過(guò)小波包分解得到各子帶信號(hào)作為源信號(hào),由FastICA 變換之后,識(shí)別其中的噪聲獨(dú)立成分并置零。腦電信號(hào)中的電參考成分為周圍設(shè)備對(duì)腦電信號(hào)的影響,電參考信號(hào)對(duì)各觀測(cè)信號(hào)的影響可能是正方向也可能是負(fù)方向的。觀測(cè)信號(hào)和源信號(hào)之間的關(guān)系是X=AS,A為混合矩陣,在FastICA 算法中可知源信號(hào)獨(dú)立成分對(duì)各觀測(cè)信號(hào)的影響與混合矩陣A列相關(guān)。因此通過(guò)分析混合矩陣找出電參考獨(dú)立成分。

    為了區(qū)分出對(duì)觀測(cè)信號(hào)有正影響或者負(fù)影響的獨(dú)立成分,計(jì)算混合矩陣中的每一列和正方向以及負(fù)方向的標(biāo)準(zhǔn)向量的矢量角,最接近正方向和負(fù)方向的兩個(gè)獨(dú)立成分置零,作為電參考噪聲去除。

    2)由于獨(dú)立成分的高斯性差別較大,故采用峰度系數(shù)對(duì)噪聲進(jìn)行識(shí)別。

    峰度系數(shù)對(duì)隨機(jī)變量的分布以及瞬時(shí)特征特別敏感,對(duì)于離散的腦電信號(hào)的計(jì)算公式如下:

    式中:xi為離散腦電信號(hào)值;xˉ是離散腦電信號(hào)的平均值;N為采樣長(zhǎng)度;σt是標(biāo)準(zhǔn)差。

    峰度系數(shù)用于區(qū)分高斯性。當(dāng)K>3 時(shí),表示該獨(dú)立成分具有超高斯性;腦電信號(hào)成分的峰度系數(shù)小于3,將峰度系數(shù)大于3 的獨(dú)立成分作為干擾成分去除。

    2 實(shí) 驗(yàn)

    本實(shí)驗(yàn)所采集的EEG 數(shù)據(jù)采用奧地利g.tec 公司所研制的g.hlamp 腦電接口設(shè)備進(jìn)行采集,采用64 導(dǎo)聯(lián)電極帽記錄腦電信號(hào),采樣頻率為256 Hz。以左右耳垂為參考電極點(diǎn),接地電極則位于前額發(fā)際處。實(shí)驗(yàn)對(duì)象是實(shí)驗(yàn)室一名健康的同學(xué),年齡24 歲,采集實(shí)驗(yàn)對(duì)象靜息狀態(tài)10 min,在此過(guò)程中,實(shí)驗(yàn)對(duì)象伴有眨眼和挑眉若干次,具體實(shí)驗(yàn)環(huán)境如圖2 所示。

    圖2 實(shí)驗(yàn)環(huán)境示意圖

    3 基于小波包分解的快速獨(dú)立成分分析去噪算法

    傳統(tǒng)的FastICA 方法是針對(duì)多通道的腦電信號(hào)進(jìn)行分離,要求源信號(hào)是統(tǒng)計(jì)獨(dú)立的,同時(shí)獨(dú)立成分必須具有非高斯性,但多通道腦電信號(hào)彼此之間是相關(guān)聯(lián)的,并且腦電信號(hào)都是復(fù)雜的非線性結(jié)構(gòu),腦電信號(hào)中不同頻率段自發(fā)性腦電信號(hào)之間也存在著某種關(guān)聯(lián)性,不同導(dǎo)聯(lián)的噪聲信號(hào)之間也存在關(guān)聯(lián)性。這些因素都將導(dǎo)致FastICA 對(duì)多通道腦電信號(hào)進(jìn)行獨(dú)立成分提取結(jié)果并不理想,因?yàn)楠?dú)立成分中不僅包括各種噪聲成分,還包括θ波、α波、β波等不同頻率段的自發(fā)性腦電信號(hào),這些信號(hào)之間并不獨(dú)立。多通道腦電信號(hào)獨(dú)立成分如圖3所示。

    為了改善去噪效果,將小波包分析和FastICA 算法相結(jié)合,對(duì)任一通道的腦電信號(hào)先進(jìn)行小波包分解,得到頻帶較窄的子帶信號(hào),不僅增強(qiáng)了觀測(cè)信號(hào)的非高斯性,而且保證了源信號(hào)具有一定的獨(dú)立性,有效地去除了因?yàn)椴煌X電信號(hào)之間的統(tǒng)計(jì)關(guān)聯(lián)性帶來(lái)的影響,因此用小波包分解實(shí)現(xiàn)子帶信號(hào)FastICA 算法,可以更有效地分離出噪聲獨(dú)立成分。

    圖3 多通道腦電信號(hào)獨(dú)立分量

    依據(jù)源信號(hào)的線性疊加性,將源信號(hào)看作一系列子帶成分的和,如式(10)所示:

    圖4 基于小波包分解的FastICA

    4 實(shí)驗(yàn)分析

    實(shí)驗(yàn)獲得的腦電信號(hào)比較斑駁,因此要對(duì)每一組腦電信號(hào)的數(shù)據(jù)進(jìn)行腦電分段、帶通濾波、去除偽跡、傅里葉變換、功率譜估計(jì)等離線分析。電極所選的腦區(qū)為額葉Fz、FCz、Cz 以及附近的電極,頂葉CPz、Pz 以及附近的電極等15 個(gè)導(dǎo)聯(lián)。

    1)在采集到的多通道腦電信號(hào)中,任選一個(gè)通道的腦電信號(hào),如圖5 所示。

    圖5 單通道原始腦電信號(hào)

    2)對(duì)單通道腦電信號(hào)進(jìn)行4 層小波包分解,得到不同尺度下的高低頻細(xì)節(jié)分量,對(duì)腦電信號(hào)的細(xì)節(jié)做出進(jìn)一步的劃分,并且作為FastICA 算法的觀測(cè)信號(hào),增強(qiáng)了非高斯性以及實(shí)現(xiàn)了初步的相關(guān),將部分噪聲和腦電信號(hào)分離開如圖6 所示。

    3)對(duì)小波包分解之后的分量信號(hào)進(jìn)行快速獨(dú)立成分分析處理,盡可能地分解出同一頻率段下和腦電信號(hào)混合在一起的噪聲分量,最大可能地保證各分量之間的獨(dú)立性。

    4)根據(jù)噪聲成分和電參考識(shí)別方法去除和腦電信號(hào)無(wú)關(guān)的獨(dú)立分量,在圖7 左側(cè)標(biāo)出。

    5)將剩下的獨(dú)立分量經(jīng)過(guò)獨(dú)立成分分析重構(gòu)和小波包重構(gòu)融合成新的觀測(cè)信號(hào),并截取其中某一小段的采樣數(shù)據(jù)便于分析,如圖8 所示。

    5 不同去噪方法的對(duì)比

    用Matlab 對(duì)所獲得的腦電信號(hào)進(jìn)行處理,從而對(duì)比單通道基于小波包分解的FastICA、多通道FastICA 去噪方法的好壞,由于腦電信號(hào)采集的通道較多,因此任選一個(gè)通道腦電信號(hào),從以下三個(gè)方面進(jìn)行比較:進(jìn)行時(shí)域分析;進(jìn)行功率譜密度分析;分別計(jì)算兩種去噪方法的信噪比和均方根誤差。對(duì)比效果如圖9~圖12 所示。

    圖7 單通道腦電信號(hào)獨(dú)立分量

    圖8 去噪之后的腦電信號(hào)

    圖9 去噪前的腦電波形

    圖10 多通道FastICA 去噪后的腦電波形

    為了進(jìn)一步地比較去噪效果,采用信噪比(SNR)和均方根誤差(RMSE)這兩個(gè)指標(biāo)進(jìn)行衡量:式中:si是原始EEG 信號(hào);si為去噪之后的EEG 信號(hào);N為信號(hào)長(zhǎng)度。信號(hào)的信噪比越高,去噪信號(hào)與原始信號(hào)的均方根誤差越小,說(shuō)明信號(hào)去噪的效果越好。分別統(tǒng)計(jì)分析基于小波包分解的FastICA 對(duì)單通道腦電信號(hào)去噪算法和FastICA 對(duì)多通道腦電信號(hào)去噪算法的信噪比以及均方根誤差,結(jié)果如表1 所示。

    表1 兩種方法的SNR 和RMSE 結(jié)果對(duì)比

    圖11 基于小波包分解的FastICA 去噪后腦電波形

    圖12 功率譜密度對(duì)比分析

    圖9~圖12 表示兩種方法在時(shí)域分析和功率譜分析的對(duì)比圖。表1 則代表兩種處理方法的信噪比和均方根誤差。從時(shí)域分析、功率譜密度分析以及信噪比和均方根誤差的角度來(lái)對(duì)比基于小波包的FastICA、多通道FastICA 兩種方法去噪效果的好壞,可知:

    1)在時(shí)域分析圖上,經(jīng)過(guò)本研究去噪之后的腦電信號(hào)幅值低于多通道FastICA 去噪之后的腦電信號(hào),多通道FastICA 雖然對(duì)原始信號(hào)去除了一些噪聲,但是毛刺依舊存在,很雜亂并且多了一些毛刺,說(shuō)明在去噪的過(guò)程中存在一定的失真,去除的噪聲中混有有效成分,而本研究的去噪方法去除了很多毛刺,而且波形較為平滑,腦電信號(hào)中的有效成分也得到了較好的保留。

    2)通過(guò)分析功率譜密度圖,發(fā)現(xiàn)本研究的去噪方法明顯地削弱了噪聲成分,較好地保留了有效成分,而且原本在噪聲信號(hào)掩蓋下的各頻率的腦電特征也顯現(xiàn)出來(lái),功率譜密度圖形呈現(xiàn)出更多的突起,即為腦電信號(hào)的特征。而多通道的FastICA 去噪之后的腦電信號(hào),沒(méi)有明顯的削弱噪聲,還去除了部分有效信號(hào),使腦電信號(hào)出現(xiàn)了失真。

    3)相比多通道FastICA,基于小波包分解的FastICA 更有效地去除了腦電信號(hào)中的噪聲成分,較好地保留了腦電信號(hào),同時(shí)凸顯出更多的腦電特征。

    4)但是僅從時(shí)域分析圖和功率譜密度分析圖去比較兩種方法的電壓幅值以及功率譜值的高低,并不能完全判斷出處理方法的效果,還要觀測(cè)曲線的波動(dòng)穩(wěn)定程度。本研究的去噪方法,曲線出現(xiàn)大大小小的突出,使原本被噪聲信號(hào)掩蓋的腦電信號(hào)凸顯出來(lái)。

    5)一味地降低電壓幅值和功率譜值,只是處理方法強(qiáng)弱的一方面。在去噪聲的過(guò)程中,保證噪聲的有效去除而且保證其中純凈的腦電信號(hào)不受干擾也很重要,因此通過(guò)SNR 和RMSE 對(duì)兩種方法的處理效果進(jìn)行對(duì)比,信號(hào)的SNR 越高,去噪信號(hào)與原始信號(hào)的RMSE 越小說(shuō)明信號(hào)去噪的效果越好。由此可以發(fā)現(xiàn),本研究提出的基于小波包分解的FastICA 對(duì)腦電信號(hào)的去噪效果較為有效。

    6 結(jié) 語(yǔ)

    基于小波包分解的單通道腦電信號(hào)FastICA 與多通道腦電信號(hào)FastICA 相比,其不僅有效地改善了熱刺激中由于長(zhǎng)時(shí)間的采集導(dǎo)致多通道腦電信號(hào)數(shù)據(jù)量大而帶來(lái)的處理上的復(fù)雜程度,同時(shí)也避免了多通道腦電信號(hào)之間的統(tǒng)計(jì)相關(guān)性,以及自發(fā)性腦電信號(hào)之間的關(guān)聯(lián)性對(duì)FastICA 這種新型盲源分離方法的干擾,優(yōu)化了FastICA 對(duì)本實(shí)驗(yàn)的處理效果,而且原本不易去噪的單通道腦電信號(hào)采用該方法之后,去噪的效果很明顯,降低了信號(hào)的電壓幅值,改善了雜亂的波形使其更加平滑,并且提高了信噪比和魯棒性。綜合上述條件表明,本文算法在噪聲和腦電信號(hào)之間的分離、去噪、保真性方面,處理速度更快、效果更好。

    猜你喜歡
    單通道波包電信號(hào)
    基于聯(lián)合聚類分析的單通道腹部心電信號(hào)的胎心率提取
    基于Code Composer Studio3.3完成對(duì)心電信號(hào)的去噪
    基于小波包Tsallis熵和RVM的模擬電路故障診斷
    基于隨機(jī)森林的航天器電信號(hào)多分類識(shí)別方法
    基于擴(kuò)頻碼周期性的單通道直擴(kuò)通信半盲分離抗干擾算法
    基于小波包變換的電力系統(tǒng)諧波分析
    小波包理論與圖像小波包分解
    采用6.25mm×6.25mm×1.8mm LGA封裝的雙通道2.5A、單通道5A超薄微型模塊穩(wěn)壓器
    基于小波包的全信息解調(diào)方法及其應(yīng)用
    基于生物電信號(hào)的駕駛疲勞檢測(cè)方法
    青青草视频在线视频观看| 久久久精品区二区三区| 亚洲第一av免费看| 成人免费观看视频高清| 国产精品 国内视频| 亚洲国产av影院在线观看| 精品第一国产精品| 久久人人爽人人片av| 美女午夜性视频免费| 成人亚洲精品一区在线观看| 国产一区二区三区av在线| 免费av中文字幕在线| 又黄又粗又硬又大视频| 91字幕亚洲| 欧美人与善性xxx| 亚洲七黄色美女视频| 久久精品久久久久久噜噜老黄| 青春草亚洲视频在线观看| 午夜福利影视在线免费观看| tube8黄色片| 性高湖久久久久久久久免费观看| 国产极品粉嫩免费观看在线| 热re99久久国产66热| 免费看av在线观看网站| 国产在线视频一区二区| 日韩伦理黄色片| 久久久久久久久免费视频了| 激情视频va一区二区三区| 久久亚洲精品不卡| 国产成人精品久久二区二区91| 七月丁香在线播放| 王馨瑶露胸无遮挡在线观看| 97在线人人人人妻| 中文字幕人妻丝袜制服| 久久久久国产一级毛片高清牌| 亚洲国产欧美网| 欧美日韩亚洲高清精品| 五月开心婷婷网| 一区在线观看完整版| 欧美大码av| 丝袜脚勾引网站| 国产精品麻豆人妻色哟哟久久| 国产精品一区二区在线观看99| 婷婷色综合大香蕉| 99久久综合免费| 亚洲专区中文字幕在线| 欧美日韩综合久久久久久| 日本一区二区免费在线视频| 亚洲国产欧美网| 日韩制服丝袜自拍偷拍| 丝袜美腿诱惑在线| 精品欧美一区二区三区在线| 精品福利观看| 精品一品国产午夜福利视频| 最新的欧美精品一区二区| √禁漫天堂资源中文www| 2021少妇久久久久久久久久久| av天堂在线播放| 久久久国产精品麻豆| 后天国语完整版免费观看| 丰满迷人的少妇在线观看| 欧美成人午夜精品| 十八禁网站网址无遮挡| 日韩制服骚丝袜av| 欧美国产精品va在线观看不卡| 亚洲精品一区蜜桃| 日本a在线网址| 成人国语在线视频| 成年人午夜在线观看视频| 国产一区有黄有色的免费视频| 欧美日韩亚洲综合一区二区三区_| 精品国产一区二区久久| 国产精品一区二区在线观看99| 国产成人av教育| 一区二区三区四区激情视频| 十八禁高潮呻吟视频| 免费高清在线观看日韩| 欧美黄色淫秽网站| 曰老女人黄片| 亚洲男人天堂网一区| 免费高清在线观看视频在线观看| 亚洲色图综合在线观看| 在线亚洲精品国产二区图片欧美| 久久中文字幕一级| 精品欧美一区二区三区在线| 女性被躁到高潮视频| av在线老鸭窝| 日韩av不卡免费在线播放| 国产精品欧美亚洲77777| 水蜜桃什么品种好| 一二三四社区在线视频社区8| 成年人黄色毛片网站| 色视频在线一区二区三区| 老司机靠b影院| 日韩,欧美,国产一区二区三区| 成人黄色视频免费在线看| 看免费av毛片| 新久久久久国产一级毛片| 叶爱在线成人免费视频播放| 成人国语在线视频| 国产片特级美女逼逼视频| 青春草视频在线免费观看| 国产女主播在线喷水免费视频网站| 欧美日韩亚洲国产一区二区在线观看 | 一本—道久久a久久精品蜜桃钙片| 国产成人免费无遮挡视频| 18禁观看日本| 在线看a的网站| 国产男女内射视频| 建设人人有责人人尽责人人享有的| 国产一卡二卡三卡精品| 亚洲一卡2卡3卡4卡5卡精品中文| 又黄又粗又硬又大视频| 1024香蕉在线观看| 女警被强在线播放| 一本综合久久免费| 青春草视频在线免费观看| 久久天躁狠狠躁夜夜2o2o | 捣出白浆h1v1| 亚洲免费av在线视频| 2021少妇久久久久久久久久久| 日韩 欧美 亚洲 中文字幕| 精品福利永久在线观看| 成人18禁高潮啪啪吃奶动态图| 天堂俺去俺来也www色官网| 亚洲精品美女久久久久99蜜臀 | 天天躁狠狠躁夜夜躁狠狠躁| 色精品久久人妻99蜜桃| 国产日韩一区二区三区精品不卡| 韩国精品一区二区三区| 欧美精品人与动牲交sv欧美| av国产精品久久久久影院| 女人被躁到高潮嗷嗷叫费观| 久久中文字幕一级| 久久影院123| 欧美亚洲日本最大视频资源| 国产熟女欧美一区二区| 制服人妻中文乱码| 亚洲图色成人| 国产精品久久久久久人妻精品电影 | 性色av乱码一区二区三区2| 国产亚洲欧美精品永久| 免费看av在线观看网站| 欧美日韩综合久久久久久| 亚洲综合色网址| 久久国产精品影院| 90打野战视频偷拍视频| 黑人猛操日本美女一级片| 在线观看www视频免费| 国产欧美亚洲国产| 两性夫妻黄色片| 日韩大片免费观看网站| 一区二区日韩欧美中文字幕| 成在线人永久免费视频| 大香蕉久久成人网| 国产精品亚洲av一区麻豆| 亚洲av综合色区一区| 国产视频首页在线观看| 一区二区三区四区激情视频| 美国免费a级毛片| 丝袜在线中文字幕| 欧美人与性动交α欧美软件| 亚洲精品国产av蜜桃| 热99久久久久精品小说推荐| 久久亚洲国产成人精品v| av欧美777| 热99久久久久精品小说推荐| 久久中文字幕一级| 亚洲五月婷婷丁香| 免费不卡黄色视频| 黑人欧美特级aaaaaa片| 9191精品国产免费久久| 自线自在国产av| 亚洲成人国产一区在线观看 | 国产视频首页在线观看| 9热在线视频观看99| 亚洲精品久久成人aⅴ小说| 欧美黑人欧美精品刺激| 性色av乱码一区二区三区2| 欧美人与性动交α欧美精品济南到| 99久久精品国产亚洲精品| xxxhd国产人妻xxx| 久久久久久久国产电影| 欧美日本中文国产一区发布| 亚洲一区二区三区欧美精品| 久久久久视频综合| 高清欧美精品videossex| 亚洲少妇的诱惑av| 欧美激情极品国产一区二区三区| 最新的欧美精品一区二区| 18禁裸乳无遮挡动漫免费视频| 黄色a级毛片大全视频| 丰满饥渴人妻一区二区三| 久久久精品免费免费高清| 国产精品av久久久久免费| 久久久精品国产亚洲av高清涩受| 久久99热这里只频精品6学生| 男人爽女人下面视频在线观看| 国产野战对白在线观看| 天天添夜夜摸| 午夜老司机福利片| 精品少妇内射三级| av天堂在线播放| 久久人妻熟女aⅴ| 男女高潮啪啪啪动态图| 一级黄片播放器| 午夜精品国产一区二区电影| 欧美人与善性xxx| 久久99精品国语久久久| 国产亚洲欧美精品永久| 2021少妇久久久久久久久久久| 黄色 视频免费看| 亚洲精品一二三| 国产成人影院久久av| 又粗又硬又长又爽又黄的视频| 一级毛片我不卡| 老司机靠b影院| 成年人免费黄色播放视频| 美女大奶头黄色视频| 一区二区av电影网| 少妇猛男粗大的猛烈进出视频| 少妇粗大呻吟视频| 亚洲中文av在线| 人人妻人人澡人人看| 午夜日韩欧美国产| 亚洲av在线观看美女高潮| 黄频高清免费视频| 91麻豆av在线| 欧美成人午夜精品| 国产在线免费精品| 日本vs欧美在线观看视频| 高清欧美精品videossex| 日韩欧美一区视频在线观看| 大话2 男鬼变身卡| 亚洲国产精品一区三区| 久久久久久久国产电影| 精品国产一区二区三区久久久樱花| 久久久久国产一级毛片高清牌| 老司机在亚洲福利影院| 亚洲精品一二三| 欧美人与性动交α欧美精品济南到| 亚洲欧美一区二区三区久久| 欧美在线黄色| 亚洲精品自拍成人| 午夜免费鲁丝| 自线自在国产av| 亚洲国产精品一区二区三区在线| 国产高清videossex| 飞空精品影院首页| 午夜福利在线免费观看网站| 亚洲成国产人片在线观看| 久久精品亚洲熟妇少妇任你| 91麻豆精品激情在线观看国产 | 男女无遮挡免费网站观看| a级毛片在线看网站| 久久午夜综合久久蜜桃| 亚洲熟女毛片儿| 黄色一级大片看看| 极品人妻少妇av视频| 精品少妇内射三级| 免费黄频网站在线观看国产| 欧美xxⅹ黑人| 国产精品一区二区在线不卡| 亚洲 国产 在线| 女人久久www免费人成看片| 操美女的视频在线观看| 国产精品99久久99久久久不卡| 日本91视频免费播放| 亚洲成色77777| 国产一区二区在线观看av| 天天躁日日躁夜夜躁夜夜| 在线亚洲精品国产二区图片欧美| 老汉色av国产亚洲站长工具| 亚洲一区二区三区欧美精品| 美女扒开内裤让男人捅视频| 国产在视频线精品| 日本五十路高清| 国产高清不卡午夜福利| 亚洲av成人精品一二三区| 99国产精品99久久久久| 飞空精品影院首页| 黄色视频不卡| 国产黄色免费在线视频| av一本久久久久| 久久久久国产精品人妻一区二区| 国产在线观看jvid| 亚洲欧美精品自产自拍| 亚洲精品乱久久久久久| 免费在线观看视频国产中文字幕亚洲 | 日韩中文字幕视频在线看片| 高清欧美精品videossex| 丝袜喷水一区| 精品国产乱码久久久久久男人| 中文字幕色久视频| 欧美国产精品va在线观看不卡| 精品一区二区三卡| av网站免费在线观看视频| 日韩av免费高清视频| 国产精品99久久99久久久不卡| 亚洲欧美色中文字幕在线| 欧美在线一区亚洲| www.熟女人妻精品国产| www.999成人在线观看| av网站免费在线观看视频| 男女高潮啪啪啪动态图| 亚洲欧美激情在线| 人妻一区二区av| 免费观看a级毛片全部| 久久久久久久久免费视频了| 婷婷色av中文字幕| 亚洲伊人色综图| 国产一区二区三区综合在线观看| 这个男人来自地球电影免费观看| 日本猛色少妇xxxxx猛交久久| 精品第一国产精品| 亚洲国产精品成人久久小说| 亚洲伊人色综图| 精品一品国产午夜福利视频| 母亲3免费完整高清在线观看| 国产精品免费大片| 国产精品99久久99久久久不卡| 久久天堂一区二区三区四区| 1024视频免费在线观看| 黄色a级毛片大全视频| 少妇粗大呻吟视频| 看十八女毛片水多多多| 亚洲国产av新网站| 丰满少妇做爰视频| 精品视频人人做人人爽| 国产黄色免费在线视频| 久久久精品94久久精品| 激情视频va一区二区三区| 国产在线视频一区二区| 国产一区二区在线观看av| 一本久久精品| 亚洲欧美色中文字幕在线| 亚洲av国产av综合av卡| 制服人妻中文乱码| 建设人人有责人人尽责人人享有的| 高清视频免费观看一区二区| 日本午夜av视频| 欧美乱码精品一区二区三区| tube8黄色片| 天天操日日干夜夜撸| 中文字幕最新亚洲高清| 在线观看www视频免费| 国产精品一区二区在线观看99| 国产一区二区在线观看av| 欧美人与性动交α欧美软件| av欧美777| 日韩大码丰满熟妇| 日日摸夜夜添夜夜爱| 国产爽快片一区二区三区| 国产男女超爽视频在线观看| 最新在线观看一区二区三区 | 黑人巨大精品欧美一区二区蜜桃| 午夜福利在线免费观看网站| 久久久久视频综合| 亚洲专区国产一区二区| 亚洲,欧美,日韩| 性少妇av在线| 手机成人av网站| 久热这里只有精品99| 青草久久国产| 啦啦啦中文免费视频观看日本| 精品人妻一区二区三区麻豆| 国产女主播在线喷水免费视频网站| 男女国产视频网站| 欧美 日韩 精品 国产| 亚洲av电影在线进入| 国产亚洲欧美精品永久| 国产精品一二三区在线看| 国产一区二区三区综合在线观看| 各种免费的搞黄视频| 97精品久久久久久久久久精品| 国产成人一区二区在线| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美乱码精品一区二区三区| 少妇被粗大的猛进出69影院| 婷婷色综合大香蕉| 精品久久久精品久久久| 赤兔流量卡办理| 观看av在线不卡| 免费在线观看视频国产中文字幕亚洲 | av国产久精品久网站免费入址| 男人添女人高潮全过程视频| 欧美亚洲 丝袜 人妻 在线| 精品卡一卡二卡四卡免费| 国产精品麻豆人妻色哟哟久久| 美女大奶头黄色视频| 看免费av毛片| 国产野战对白在线观看| 最新的欧美精品一区二区| 日日摸夜夜添夜夜爱| 国产日韩欧美在线精品| 亚洲精品自拍成人| 真人做人爱边吃奶动态| 亚洲国产毛片av蜜桃av| 又大又爽又粗| 丁香六月欧美| 日韩电影二区| 99精品久久久久人妻精品| 一区二区三区乱码不卡18| 午夜福利一区二区在线看| 又大又黄又爽视频免费| 在现免费观看毛片| 老汉色av国产亚洲站长工具| 99热网站在线观看| 亚洲三区欧美一区| 成人亚洲欧美一区二区av| 国产精品免费大片| 18禁观看日本| 亚洲av欧美aⅴ国产| 日日爽夜夜爽网站| 国产主播在线观看一区二区 | 一级a爱视频在线免费观看| 美女高潮到喷水免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 国产日韩欧美视频二区| 成人亚洲欧美一区二区av| 伦理电影免费视频| 国产在线一区二区三区精| 亚洲人成电影观看| 好男人视频免费观看在线| 久热这里只有精品99| 又大又黄又爽视频免费| 中文乱码字字幕精品一区二区三区| 捣出白浆h1v1| 久久久国产一区二区| 亚洲专区中文字幕在线| 亚洲激情五月婷婷啪啪| 亚洲成人手机| 久久久精品免费免费高清| 日韩一卡2卡3卡4卡2021年| 侵犯人妻中文字幕一二三四区| 岛国毛片在线播放| 国产成人啪精品午夜网站| 丝瓜视频免费看黄片| 亚洲一区二区三区欧美精品| 女人爽到高潮嗷嗷叫在线视频| 亚洲第一青青草原| 桃花免费在线播放| 婷婷色av中文字幕| 色精品久久人妻99蜜桃| 男女边摸边吃奶| 可以免费在线观看a视频的电影网站| 欧美日韩av久久| 久久女婷五月综合色啪小说| 欧美精品av麻豆av| 欧美av亚洲av综合av国产av| 99精品久久久久人妻精品| 亚洲三区欧美一区| 精品福利永久在线观看| 9191精品国产免费久久| 久久鲁丝午夜福利片| 国产成人啪精品午夜网站| 国产野战对白在线观看| 婷婷色综合大香蕉| 黄色视频在线播放观看不卡| 成在线人永久免费视频| 50天的宝宝边吃奶边哭怎么回事| 国产免费又黄又爽又色| 国产成人91sexporn| 午夜视频精品福利| 亚洲人成电影观看| 丝袜脚勾引网站| 五月开心婷婷网| 人人澡人人妻人| 另类亚洲欧美激情| 最新在线观看一区二区三区 | 日本黄色日本黄色录像| 伦理电影免费视频| av欧美777| 91麻豆av在线| 国产一区有黄有色的免费视频| 极品人妻少妇av视频| 婷婷色av中文字幕| 久久99精品国语久久久| a 毛片基地| 中文字幕另类日韩欧美亚洲嫩草| 黄色a级毛片大全视频| 久久久精品免费免费高清| 国产精品亚洲av一区麻豆| 老汉色∧v一级毛片| 亚洲精品自拍成人| 老司机亚洲免费影院| 性高湖久久久久久久久免费观看| 国产精品一区二区免费欧美 | 欧美在线黄色| 亚洲av男天堂| 自线自在国产av| 国产精品偷伦视频观看了| 手机成人av网站| 国产伦人伦偷精品视频| 久久 成人 亚洲| 一本一本久久a久久精品综合妖精| 国产亚洲欧美精品永久| av在线app专区| 国产国语露脸激情在线看| 婷婷成人精品国产| 国产黄色免费在线视频| 在线 av 中文字幕| 精品国产乱码久久久久久男人| 国产免费视频播放在线视频| 欧美黄色片欧美黄色片| 男女边吃奶边做爰视频| 搡老岳熟女国产| 黄色片一级片一级黄色片| 黄色 视频免费看| 只有这里有精品99| 国产精品一区二区精品视频观看| 性色av乱码一区二区三区2| netflix在线观看网站| 久久精品国产亚洲av涩爱| 欧美精品av麻豆av| 精品一品国产午夜福利视频| 在线av久久热| 好男人电影高清在线观看| 婷婷色综合大香蕉| 建设人人有责人人尽责人人享有的| 国产成人精品久久二区二区91| 又黄又粗又硬又大视频| 免费高清在线观看日韩| 啦啦啦视频在线资源免费观看| 99re6热这里在线精品视频| 久久综合国产亚洲精品| kizo精华| 国产亚洲午夜精品一区二区久久| 日韩 欧美 亚洲 中文字幕| 日韩大片免费观看网站| 亚洲图色成人| 国产麻豆69| 亚洲 国产 在线| 操出白浆在线播放| 欧美 日韩 精品 国产| 波多野结衣一区麻豆| 亚洲精品在线美女| 日韩大码丰满熟妇| 日韩一区二区三区影片| 欧美日韩成人在线一区二区| 一级黄片播放器| 一级,二级,三级黄色视频| 亚洲欧美一区二区三区国产| 国产女主播在线喷水免费视频网站| 日韩免费高清中文字幕av| 大陆偷拍与自拍| 国产免费一区二区三区四区乱码| 久久精品人人爽人人爽视色| 久久国产亚洲av麻豆专区| 久久鲁丝午夜福利片| av在线播放精品| 国产一区二区激情短视频 | 日日夜夜操网爽| 国产欧美日韩精品亚洲av| 久久久久久久精品精品| 午夜福利,免费看| 久久热在线av| 一本一本久久a久久精品综合妖精| 天堂8中文在线网| 亚洲av美国av| 国产日韩一区二区三区精品不卡| 宅男免费午夜| 久久综合国产亚洲精品| 亚洲第一青青草原| 9色porny在线观看| 中文字幕人妻丝袜一区二区| 91老司机精品| 91精品三级在线观看| 啦啦啦中文免费视频观看日本| 黄色视频不卡| 91精品伊人久久大香线蕉| 国产亚洲欧美在线一区二区| 丝袜美足系列| 人体艺术视频欧美日本| 亚洲精品国产av成人精品| 中文字幕另类日韩欧美亚洲嫩草| 亚洲自偷自拍图片 自拍| 亚洲国产精品一区二区三区在线| 热99国产精品久久久久久7| 国产成人精品无人区| 久热爱精品视频在线9| 亚洲人成电影观看| 在线观看www视频免费| 午夜av观看不卡| 亚洲国产成人一精品久久久| 日韩伦理黄色片| 99久久人妻综合| 99久久精品国产亚洲精品| tube8黄色片| 熟女少妇亚洲综合色aaa.| 精品一品国产午夜福利视频| av一本久久久久| 亚洲久久久国产精品| 亚洲精品av麻豆狂野| 亚洲成色77777| 夜夜骑夜夜射夜夜干| 欧美精品一区二区大全| 国产精品99久久99久久久不卡| 免费看十八禁软件| 精品国产超薄肉色丝袜足j| 精品熟女少妇八av免费久了| 日韩制服丝袜自拍偷拍| 七月丁香在线播放| 国产日韩欧美在线精品| √禁漫天堂资源中文www| 亚洲成人国产一区在线观看 | 一级片免费观看大全| 国产成人欧美| 最新在线观看一区二区三区 | 国产福利在线免费观看视频| 亚洲欧美中文字幕日韩二区| 久久精品成人免费网站| 黄色 视频免费看| 国产一区二区在线观看av| 欧美av亚洲av综合av国产av|