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

    基于活動性檢測動態(tài)估計噪聲的心音降噪算法

    2024-01-22 06:04:40許春冬辛鵬麗應冬文李海兵
    計算機工程與設計 2024年1期
    關鍵詞:心音于小波信噪比

    許春冬,辛鵬麗,閔 源+,應冬文,2,周 靜,李海兵

    (1.江西理工大學 信息工程學院,江西 贛州 341000;2.中國科學院大學 電子電氣與通信工程學院,北京 100049)

    0 引 言

    心音信號蘊含著大量心臟和周邊血管功能狀態(tài)的生理與病理信息[1,2],是心血管疾病臨床診斷的重要依據(jù)之一[3]。然而,在進行心音信號的采集與處理時,由于心音信號能量較弱[4],致使其極易受到來自外部環(huán)境噪聲、采集設備自噪聲及人體其它生理器官雜音等噪聲的干擾[1,2,5]。這些噪聲的存在往往會導致心音信號的質(zhì)量降低,進而對后續(xù)的分析與診斷結果造成嚴重的干擾[1,6]。

    近十幾年來,國內(nèi)外學者提出了許多心音降噪的方法,主要可分為3類:基于小波的心音降噪方法、基于經(jīng)驗模態(tài)分解的心音降噪方法、基于噪聲功率估計的心音降噪方法[7-10]。雖然上述算法均取得了一定的降噪效果,但仍存在一定的不足,例如基于小波的心音降噪算法[7,8]在處理異常心音時易造成心音丟失;基于經(jīng)驗模態(tài)分解的降噪方法[9,10]容易使降噪后的心音信號出現(xiàn)失真或低、中頻噪聲殘留?;谠肼晞討B(tài)估計的方法[11]存在低頻噪聲殘留和部分心音失真的問題。

    針對上述心音降噪算法所存在的不足,本文提出了一種基于心音活動性檢測動態(tài)估計噪聲的心音降噪算法(后文中均稱為HSAD-IMCRA算法)。該方法通過設計的HSAD方法來判斷心音幀是否為基礎心音,并對基礎心音和非基礎心音中的噪聲采用不同的方法進行估計,從而克服了傳統(tǒng)噪聲估計的不足,能夠在保證心音信號不失真的情況下,實現(xiàn)最大程度的噪聲抑制。

    1 算法框架

    提出算法的流程如圖1所示。首先,該算法將時域含噪心音信號經(jīng)短時傅里葉變換(short-time Fourier transform,STFT)變換到頻域,并根據(jù)心音信號的結構特點[12],設計適合的HSAD方法,判別基礎心音(fundamental heart sound,F(xiàn)HS)與非基礎心音(non-FHS)成分;其次,分別采用IMCRA和遞歸平滑的方法進行噪聲估計;然后,通過估計的噪聲功率譜與含噪心音的前后幀信息來估計非因果先驗信噪比(signal-to-noise ratio,SNR),求取干凈心音估計器的增益函數(shù)GLSA[13];最后,結合原始相位信息,通過逆STFT(ISTFT)重構時域降噪心音。其中,采用非因果SNR估計方式代替因果SNR的估計方式,更準確地估計了FHS出現(xiàn)和截止邊緣部分及一些異常突變心音信號的先驗SNR[14,15],在一定程度上減少了FHS的失真;采用HSAD輔助噪聲估計,能更準確地跟蹤噪聲功率的變化情況,從而減少了噪聲的殘留。

    圖1 心音降噪算法流程

    2 算法原理

    2.1 MMSE-LSA算法原理

    由Ephraim等[16]提出的最小均方誤差對數(shù)幅度譜(minimum mean square error of log-spectral amplitude,MMSE-LSA)估計算法通過對聲學信號進行最優(yōu)幅度譜估計[17]來恢復干凈信號,相較于一些經(jīng)典的降噪算法,能更有效降低信號失真[18],提高信號的可分析性,故適用于心音信號的降噪處理。假設心音與噪聲不相關[13],設s(n) 為干凈心音,d(n) 為噪聲,則含噪心音y(n) 可表示為

    y(n)=s(n)+d(n)

    (1)

    對y(n) 進行短時傅里葉變換[12,13],可得

    Y(k,l)=S(k,l)+D(k,l)

    (2)

    將式(2)表示為極坐標形式[13]

    R(k,l)ejφ(k,l)=A(k,l)ejφ(k,l)+N(k,l)ej?(k,l)

    (3)

    其中,l和k分別代表幀數(shù)和頻點數(shù)的索引,Y(k,l)、S(k,l)和D(k,l) 分別是y(n)、s(n) 和d(n) 進行STFT后的結果,R(k,l) 和φ(k,l) 分別是Y(k,l) 的頻譜幅度和相位;A(k,l) 和φ(k,l) 分別是S(k,l) 的頻譜幅度和相位;N(k,l) 和?(k,l) 分別是D(k,l) 的頻譜幅度和相位。

    (4)

    (5)

    (6)

    其中,ζ1(k,l)=λS(k,l)/λD(k,l) 為先驗信噪比,γ1(k,l)=|Y(k,l)|2/λD(k,l) 為后驗信噪比,λS(k,l) 為心音功率譜,λD(k,l) 為噪聲功率譜。

    2.2 非因果先驗信噪比估計

    在增益函數(shù)GLSA(k,l) 的求解過程中,用非因果信噪比[15]代替因果信噪比,能更準確估計存在突變心音信號或脈沖干擾噪聲時[13]的噪聲功率。在非因果信噪比的估計中,利用第0至l+M幀含噪心音來估計第l幀的先驗信噪比[15],而不是局限于使用前一幀信息。

    (7)

    (8)

    (9)

    2.3 噪聲估計

    2.3.1 IMCRA噪聲估計算法

    (10)

    基于最小均方誤差(minimum mean square error,MMSE)的準則下[13],噪聲功率譜的估計如式(11)所示[13,20]

    (11)

    (12)

    α2=α+(1-α)·P(k,l)

    (13)

    式(13)中,心音存在條件概率P(k,l), 可通過貝葉斯公式[22]計算得到,如式(14)和式(15)所示

    (14)

    (15)

    第一次迭代過程中,根據(jù)判決準則對是否存在基礎心音進行粗略估計[19],如式(16)所示

    (16)

    其中,γ0和ζ0為常數(shù)閾值參數(shù),基礎心音存在時I(k,l) 為1,基礎心音不存在時I(k,l) 為0。γmin(k,l) 和ζ2(k,l) 的定義如式(17)和式(18)所示

    (17)

    (18)

    (19)

    (20)

    (21)

    (22)

    (23)

    (24)

    (25)

    (26)

    (27)

    其中,γ3為常數(shù)閾值。

    2.3.2 基于HSAD-IMCRA的噪聲功率估計

    傳統(tǒng)的基于MMSE-LSA的降噪算法采用IMCRA方法對信號噪聲功率譜進行估計,但IMCRA噪聲估計對基礎心音開始邊緣部位的先/后驗信噪比參數(shù)估計往往偏低,易造成降噪后的心音信號出現(xiàn)基礎心音成分失真。此外,若心動周期內(nèi)存在人體其它器官引起的額外噪聲,也易造成先/后驗信噪比偏高,使得噪聲殘留。而心音信號主要是由第一心音(S1)、收縮期(Sys)、第二心音(S2)、舒張期(Dia)這4個狀態(tài)組成[9],且具備一定的周期性[7,9],可對應圖2所示。其中,Sys和Dia一般不包含基礎心音信號成分,表現(xiàn)為噪聲特性,故可根據(jù)遞歸平滑法來估計并更新噪聲功率譜[20],獲得更為準確的先/后驗信噪比估計結果。而S1和S2持續(xù)時間段是由噪聲和基礎心音信號混疊構成,故采用IMCRA對噪聲進行有效追蹤,動態(tài)估計噪聲功率譜并更新[22]。需要額外指出的是,IMCRA噪聲估計過程中,采用2.2節(jié)所述非因果信噪比計算方式替換因果信噪比的計算方式,減小心音失真情況的出現(xiàn)。

    圖2 心音信號時域波形

    故本文算法中采用HSAD-IMCRA的算法來估計當前含噪心音幀的噪聲功率譜,能進一步降低估計的干凈心音的失真和噪聲殘留。所提算法噪聲功率譜估計流程如圖3所示。

    圖3 HSAD-IMCRA算法噪聲功率譜估計流程

    活動性檢測是聲學信號處理領域里較為常用的處理手段,且類別繁多[13,20]。為了在有效檢測心音活動性的前提下進一步考慮低的算法復雜度,本文采用對數(shù)頻譜距離(logarithmic amplitude spectrum distance,LASD)來設計適合心音信號的HSAD方法,用于判別基礎心音幀和非基礎心音幀[21]?;谇皫讕瑸榉腔A心音幀的假設,可通過式(28)將LASD計算如下

    (28)

    (29)

    T(l)=τD(l-1)+(τ-1)D(l),τ>1&T(l)>D(1)

    (30)

    (31)

    其中,ρ(l-1)=1表示前一幀為基礎心音幀,ρ(l-1)=0表示前一幀為非基礎心音幀,NIMCRA(k,l) 為IMCRA算法估計的噪聲幅度譜,τ為自適應閾值平滑參數(shù),αn為平滑因子[20]。

    此外,根據(jù)式(31)的原理,噪聲功率譜的平滑估計λD(k,l) 可表示如下

    (32)

    其中,λIMCRA(k,l) 為IMCRA估計的噪聲功率譜,在獲得估計的噪聲功率譜后,即可計算先/后驗信噪比,并按式(5)獲得增益函數(shù)GLSA,從而重構出干凈心音信號。實驗中所涉及的參數(shù)設置見表1。

    表1 參數(shù)的取值

    3 實驗仿真與結果分析

    3.1 實驗數(shù)據(jù)

    實驗仿真所用數(shù)據(jù)集由公開數(shù)據(jù)集(https://github.com/yaseen21khan/Classification-of-Heart-Sound-Si-gnal-Using-Multiple-Features-)[23]和自采集心音數(shù)據(jù)集[24]兩部分組成,共計1030條。參考文獻[10]指出,高斯白噪聲是心音信號中的主要噪聲成分,故本文主要以去除高斯白噪聲為主。在實驗仿真中,從構建的數(shù)據(jù)集中抽取出80條相對干凈的正常心音、320條相對干凈的異常心音,并將高斯白噪聲分別以0 dB、5 dB、10 dB、15 dB、20 dB的SNR疊加到抽取的心音信號中[10],構建含噪心音數(shù)據(jù)集,且所有心音信號均降采樣至2000 Hz。

    3.2 HSAD-IMCRA噪聲估計性能分析

    為了驗證HSAD-IMCRA算法的噪聲估計性能,將其與IMCRA噪聲估計性能進行了對比。選用的評測指標為對數(shù)似然比(log likelihood ratio,LLR)和加權譜斜率(weighted spectral slope,WSS)[13,20]。其中,LLR是對頻譜距離的度量,取值范圍在0~2之間,估計的干凈心音與原始干凈心音間的LLR越小,即噪聲估計越準確[13];WSS是對頻帶譜斜率之間加權差距的度量,WSS越小,則估計的干凈心音與原始干凈心音越接近,即噪聲估計越準確[20]。噪聲估計算法相關指標對比結果如圖4所示。

    圖4 噪聲估計算法相關指標結果對比

    圖4(a)和圖4(b)分別是基于IMCRA算法和基于HSAD-IMCRA算法估計的干凈心音與真實的干凈心音間的LLR與WSS的對比結果。顯然,基于HSAD-IMCRA算法的WSS和LLR評測結果均低于IMCRA算法,這表明基于HSAD-IMCRA算法的噪聲估計性能要優(yōu)于IMCRA算法。

    3.3 各類算法降噪結果的時頻圖對比

    提出算法與基于小波閾值的降噪算法[8]、基于OMLSA-小波的降噪算法[11]進行實驗對比。其中,基于小波閾值降噪算法中采用的小波基為coif-5,分解層數(shù)為5層,閾值函數(shù)為自適應非線性閾值;基于OMLSA-小波降噪算法中,先通過OMLSA對含噪心音降噪,再利用小波降噪算法對OMLSA降噪后的心音繼續(xù)降噪處理,其中小波降噪算法中的相關參數(shù)參考文獻[11]。算法降噪性能的評測指標為輸出SNR和均方根誤差(root mean square error,RMSE)[2]。

    采用3種降噪算法對公開數(shù)據(jù)集和自采集的心音信號進行降噪,來驗證各類算法的降噪效果。其中,公開數(shù)據(jù)集中的數(shù)據(jù)相對比較干凈,故對該類數(shù)據(jù)集中的心音進行降噪之前,將高斯白噪聲以10 dB的SNR分別疊加到所選的心音中,來構建含噪心音。而部分自采集數(shù)據(jù)集中本身含有大量的噪聲,故不對其進行加噪處理。各類算法對不同類型心音的降噪結果如圖5所示。

    圖5 3種算法對不同類型心音的降噪結果

    圖5(a)、圖5(b)和圖5(c)分別是上述3種算法對正常/異常心音降噪后的時域波形圖和頻譜圖的對比結果。其中,圖5(a)和圖5(b)分別使用的是公開數(shù)據(jù)集中的正常心音和異常心音,而圖5(c)使用的是自采集心音數(shù)據(jù)集的正常心音。為了便于對比分析,圖中用橢圓進行標記。從圖5(a)和圖5(c)的時域波形圖中可以看出,基于OMLSA-小波算法因?qū)⒉糠诌吘売杏眯盘柡图捌湮⑷醯碾s音當成噪聲信號進行估計造成后期估計的基礎心音成分有所丟失,而基于小波閾值算法和本文算法降噪后的波形與原始心音波形差別不大;然而,從圖5(a)和圖5(c)的頻譜圖中可以看出,基于小波閾值算法能通過頻帶去除高頻和低頻噪聲使非基礎心音部分殘留的噪聲較少,但基礎心音部分卻因為丟失的高頻和低頻分量中包含部分有用心音成分而造成了FHS的失真。另外從圖5(b)的時域波形圖和頻譜圖中也可以看出,基于小波閾值算法和基于OMLSA-小波算法在對噪聲去除的同時將噪聲和基礎心音信號混疊構成的異常心音信號成分一同丟失,導致出現(xiàn)較為明顯的心音失真,且基于OMLSA-小波算法降噪后的心音在非基礎心音部分仍存在部分中頻噪聲。相比上述兩種算法,本文算法可使心音信號在減小FHS失真的前提下,更為有效地抑制噪聲,重構出一個更為平滑自然的頻譜,更加符合后續(xù)分類處理的需求。綜上時頻圖的對比分析來看,提出算法較基于小波閾值的算法和基于OMLSA-小波的算法更優(yōu)。

    3.4 各類算法降噪結果的指標評測

    此外,本文對3種算法的降噪性能進行了客觀的對比,評測所用心音包括正常心音(N)與異常心音。其中,異常心音為二尖瓣反流(mitral regurgitation,MR)心音、二尖瓣脫垂(mitral valve prolapse,MVP)心音、二尖瓣狹窄(mitral stenosis,MS)心音和主動脈狹窄(aortic stenosis,AS)心音4類[23]。3種算法對不同信噪比的正常與異常心音降噪后輸出SNR的統(tǒng)計平均值見表2。輸出SNR值越大,表示算法的降噪性能越好[2,6,7]。

    表2 對不同信噪比的含噪心音降噪后的SNR

    從表2中可以看出,基于OMLSA-小波算法和基于小波閾值算法要明顯低于提出算法的輸出SNR;原因是前兩種算法降噪后的心音信號殘留了較多的中頻噪聲,且存在FHS失真,尤其是異常心音更為明顯。而提出算法對SNR為0 dB的MVP所得SNR值要低于基于OMLSA-小波算法,這是源于提出算法中的IMCRA為了防止噪聲過估計,使部分噪聲殘留,而MVP相對其它類型的心音殘留的噪聲較多,但整體上對多種類型的心音有著較好的降噪效果,且能在保證心音信號不失真的前提下,對噪聲最大程度的抑制。

    鑒于SNR的部分不良表現(xiàn),故進一步通過RMSE來檢驗3種算法在降噪中對正常心音與異常心音的失真與降噪平衡度。RMSE越小,表示降噪后的心音越接近于真實干凈心音,即降噪性能越好[2,6,7],反之降噪性能越差。從數(shù)據(jù)中抽選N、MVP、MS和MR這4種類型的心音,利用3種算法分別對不同信噪比下的4類心音進行降噪處理,其RMSE對比結果[2]如圖6所示。

    圖6 不同信噪比的4類心音降噪后的RMSE對比曲線

    從圖6中可以看出,隨著噪聲水平的降低,提出算法的RMSE也逐漸減小。圖6(a)顯示,SNR為20 dB時,提出算法降噪后的RMSE為0.0124,說明降噪后的信號與原始信號之間相似度較高,即降噪效果較好;SNR為0 dB時,提出算法降噪后的RMSE為0.0399,降噪后的信號與原始信號之間相似度降低,降噪性能下降;這是因為低SNR下,為了確保FHS無失真引起了部分噪聲殘留。從圖6(a)、圖6(b)、圖6(c)和圖6(d)可以看出,除圖6(c)中因防止噪聲過估計造成MVP信號內(nèi)殘留部分噪聲,從而使提出算法在0 dB~5 dB間的RMSE值略高于基于OMLSA-小波的算法之外,其余情況下,提出算法對N、MR、MVP和MS這4種類型心音進行降噪所得RMSE值均低于對比算法,說明提出算法降噪后的信號更接近于原始干凈心音信號,能較好地保留FHS成分,更為適用與不同噪聲水平下的多種類型心音降噪。

    4 結束語

    本文根據(jù)心音信號的基本特征,設計HSAD方法判決當前心音幀是否為基礎心音幀,并分別采用遞歸平滑和IMCRA動態(tài)對非基礎心音和基礎心音的噪聲功率進行估計與更新。該方法一方面結合HSAD,來減小收縮期和舒張期雜音對噪聲估計的影響,彌補IMCRA的不足,減小雜音和噪聲的殘留;另一方面采用非因果信噪比替代因果信噪比,避免了由非基礎心音向基礎心音過渡時先/后驗信噪比偏小的問題,減小了基礎心音成分的失真。

    猜你喜歡
    心音于小波信噪比
    基于深度學習的無人機數(shù)據(jù)鏈信噪比估計算法
    基于小波去噪的稱重雨量數(shù)據(jù)分析
    低信噪比下LFMCW信號調(diào)頻參數(shù)估計
    電子測試(2018年11期)2018-06-26 05:56:02
    一種新的基于小波基的時變信道估計
    鐵道學報(2018年5期)2018-06-21 06:21:08
    基于雙閾值的心音快速分段算法及其應用研究
    基于小波和Hu 矩的颮線雷達回波識別
    低信噪比下基于Hough變換的前視陣列SAR稀疏三維成像
    雷達學報(2017年3期)2018-01-19 02:01:27
    雙聲道心音能量熵比的提取與識別研究
    基于香農(nóng)熵的心音信號檢測方法研究
    保持信噪比的相位分解反褶積方法研究
    亚洲经典国产精华液单| 九色亚洲精品在线播放| 狠狠精品人妻久久久久久综合| 午夜福利乱码中文字幕| 亚洲经典国产精华液单| 熟女电影av网| 永久网站在线| 成人综合一区亚洲| 亚洲精品第二区| 久久久久久久精品精品| 国产成人精品婷婷| 全区人妻精品视频| 考比视频在线观看| 久久久国产欧美日韩av| 亚洲成av片中文字幕在线观看 | 人人妻人人爽人人添夜夜欢视频| 成年人免费黄色播放视频| 久久久精品免费免费高清| 亚洲,一卡二卡三卡| 国产av码专区亚洲av| 国产又爽黄色视频| 国产精品99久久99久久久不卡 | 亚洲成国产人片在线观看| 久久久久久人妻| 男人舔女人的私密视频| 青春草视频在线免费观看| 国产欧美亚洲国产| 国产熟女欧美一区二区| 国产乱人偷精品视频| 咕卡用的链子| 一级片'在线观看视频| 色婷婷av一区二区三区视频| 国产一区二区三区av在线| 中文字幕人妻丝袜制服| 免费播放大片免费观看视频在线观看| 免费久久久久久久精品成人欧美视频 | 日韩av在线免费看完整版不卡| 亚洲精品美女久久av网站| 国产精品蜜桃在线观看| 午夜精品国产一区二区电影| 日韩一区二区三区影片| a级毛色黄片| 99久久精品国产国产毛片| 精品亚洲乱码少妇综合久久| 欧美xxⅹ黑人| 少妇人妻 视频| 五月玫瑰六月丁香| 久久久久久久久久久久大奶| 久久鲁丝午夜福利片| 少妇高潮的动态图| 亚洲经典国产精华液单| 成年人免费黄色播放视频| 捣出白浆h1v1| 亚洲经典国产精华液单| 女性生殖器流出的白浆| 日韩 亚洲 欧美在线| 国产精品久久久久久av不卡| 高清不卡的av网站| 精品99又大又爽又粗少妇毛片| 人人妻人人爽人人添夜夜欢视频| 九九在线视频观看精品| 国产毛片在线视频| 日韩制服丝袜自拍偷拍| 午夜精品国产一区二区电影| a级毛色黄片| 少妇猛男粗大的猛烈进出视频| av免费在线看不卡| 精品卡一卡二卡四卡免费| 国语对白做爰xxxⅹ性视频网站| 日本猛色少妇xxxxx猛交久久| 国产熟女午夜一区二区三区| 国产爽快片一区二区三区| 国产精品国产av在线观看| 国产精品 国内视频| 亚洲熟女精品中文字幕| 精品一区二区三区四区五区乱码 | 亚洲情色 制服丝袜| 中文字幕最新亚洲高清| 制服人妻中文乱码| 亚洲欧美清纯卡通| 成年女人在线观看亚洲视频| 午夜av观看不卡| 国产国拍精品亚洲av在线观看| 精品一区在线观看国产| 搡女人真爽免费视频火全软件| 亚洲av欧美aⅴ国产| 在线看a的网站| 寂寞人妻少妇视频99o| 成年美女黄网站色视频大全免费| 国产白丝娇喘喷水9色精品| 国产精品一国产av| 大片免费播放器 马上看| 精品人妻一区二区三区麻豆| 人妻一区二区av| 少妇高潮的动态图| 伊人久久国产一区二区| 成年av动漫网址| 80岁老熟妇乱子伦牲交| 熟女av电影| av视频免费观看在线观看| 日韩一区二区三区影片| 麻豆精品久久久久久蜜桃| 国产精品久久久久久久久免| 男人添女人高潮全过程视频| 18禁观看日本| 久久久久久人人人人人| 啦啦啦啦在线视频资源| 9热在线视频观看99| 午夜福利视频在线观看免费| 看免费成人av毛片| 日本与韩国留学比较| 纵有疾风起免费观看全集完整版| 色94色欧美一区二区| 91精品三级在线观看| 成年人免费黄色播放视频| 十八禁网站网址无遮挡| 亚洲国产欧美日韩在线播放| 下体分泌物呈黄色| av天堂久久9| 久久久久久久久久成人| 久久国产精品大桥未久av| 91午夜精品亚洲一区二区三区| 国产男女超爽视频在线观看| 日本免费在线观看一区| 日韩人妻精品一区2区三区| 日韩av不卡免费在线播放| 亚洲熟女精品中文字幕| 制服诱惑二区| 香蕉丝袜av| 啦啦啦在线观看免费高清www| 在线 av 中文字幕| 巨乳人妻的诱惑在线观看| 久久狼人影院| kizo精华| 午夜视频国产福利| 欧美+日韩+精品| 欧美激情极品国产一区二区三区 | 日本av免费视频播放| 久久99精品国语久久久| 国精品久久久久久国模美| 久久精品国产亚洲av天美| 一级毛片黄色毛片免费观看视频| 99热网站在线观看| 欧美日韩视频高清一区二区三区二| 少妇人妻 视频| 日韩制服骚丝袜av| 男女国产视频网站| 校园人妻丝袜中文字幕| 啦啦啦中文免费视频观看日本| 看免费av毛片| 国产精品嫩草影院av在线观看| 中文字幕av电影在线播放| 久久久久久久精品精品| 高清不卡的av网站| av卡一久久| 在线观看免费日韩欧美大片| 亚洲av.av天堂| 亚洲精品久久午夜乱码| 99久国产av精品国产电影| 久久久久久人妻| 看十八女毛片水多多多| 婷婷色综合大香蕉| 成人影院久久| a级毛片黄视频| 秋霞伦理黄片| 欧美精品一区二区大全| 多毛熟女@视频| 韩国av在线不卡| 欧美老熟妇乱子伦牲交| 久久99一区二区三区| 国产男人的电影天堂91| 中文乱码字字幕精品一区二区三区| 夜夜爽夜夜爽视频| 免费高清在线观看日韩| 中文天堂在线官网| 蜜臀久久99精品久久宅男| 国产亚洲最大av| 亚洲婷婷狠狠爱综合网| 天堂8中文在线网| 麻豆精品久久久久久蜜桃| 国产精品99久久99久久久不卡 | 少妇熟女欧美另类| 91成人精品电影| 少妇的逼水好多| 久久久久久久久久人人人人人人| 免费女性裸体啪啪无遮挡网站| 午夜日本视频在线| av在线播放精品| 国产精品久久久av美女十八| 成年美女黄网站色视频大全免费| 国产成人免费观看mmmm| 中文欧美无线码| 免费观看av网站的网址| 日韩人妻精品一区2区三区| 美女大奶头黄色视频| 99香蕉大伊视频| 中国三级夫妇交换| 欧美性感艳星| 国语对白做爰xxxⅹ性视频网站| 婷婷色av中文字幕| 天堂俺去俺来也www色官网| 啦啦啦啦在线视频资源| 精品人妻一区二区三区麻豆| 99久久人妻综合| 又黄又爽又刺激的免费视频.| 大香蕉久久网| 狠狠婷婷综合久久久久久88av| 亚洲欧洲国产日韩| 国产亚洲最大av| 搡老乐熟女国产| 99久久综合免费| 热re99久久国产66热| 观看美女的网站| 国产成人免费观看mmmm| 国国产精品蜜臀av免费| 精品一区二区三区视频在线| 亚洲精品国产av蜜桃| 欧美老熟妇乱子伦牲交| 成人毛片60女人毛片免费| 成人亚洲欧美一区二区av| 国产国语露脸激情在线看| 90打野战视频偷拍视频| av免费观看日本| 人人妻人人澡人人看| 制服人妻中文乱码| 王馨瑶露胸无遮挡在线观看| 在线观看www视频免费| 日本欧美视频一区| 久久久国产精品麻豆| av.在线天堂| 国产黄色免费在线视频| av福利片在线| av天堂久久9| 美女xxoo啪啪120秒动态图| 日日摸夜夜添夜夜爱| 久久久久久久大尺度免费视频| tube8黄色片| 宅男免费午夜| 亚洲一码二码三码区别大吗| 18在线观看网站| 少妇人妻 视频| 51国产日韩欧美| 亚洲经典国产精华液单| 成人国产av品久久久| 国产成人午夜福利电影在线观看| 日韩成人伦理影院| 色94色欧美一区二区| 日本-黄色视频高清免费观看| 精品人妻偷拍中文字幕| 99热国产这里只有精品6| 曰老女人黄片| 99视频精品全部免费 在线| 超色免费av| av.在线天堂| 这个男人来自地球电影免费观看 | 啦啦啦视频在线资源免费观看| 熟女人妻精品中文字幕| 中文字幕免费在线视频6| 蜜桃在线观看..| 只有这里有精品99| 大码成人一级视频| av在线播放精品| 成人二区视频| 欧美97在线视频| 国产毛片在线视频| 国产黄色免费在线视频| 成人国语在线视频| 黄色毛片三级朝国网站| 中文字幕免费在线视频6| 国产不卡av网站在线观看| 天天操日日干夜夜撸| 亚洲国产av影院在线观看| av国产精品久久久久影院| 中文字幕精品免费在线观看视频 | 亚洲三级黄色毛片| 深夜精品福利| 中文欧美无线码| 看免费av毛片| 蜜桃国产av成人99| 三上悠亚av全集在线观看| 午夜视频国产福利| 国产欧美亚洲国产| 观看av在线不卡| 成人漫画全彩无遮挡| 亚洲四区av| 久久精品国产a三级三级三级| 久久久久国产网址| 波多野结衣一区麻豆| 看免费av毛片| 亚洲,欧美,日韩| 国产亚洲午夜精品一区二区久久| 亚洲伊人色综图| 久久久久久久久久人人人人人人| 色网站视频免费| 满18在线观看网站| 大码成人一级视频| 欧美精品av麻豆av| 男女国产视频网站| 欧美精品亚洲一区二区| 新久久久久国产一级毛片| 夫妻性生交免费视频一级片| av片东京热男人的天堂| 成年av动漫网址| 精品人妻偷拍中文字幕| 男女国产视频网站| 97精品久久久久久久久久精品| 国产精品人妻久久久久久| 90打野战视频偷拍视频| 日韩大片免费观看网站| 久久 成人 亚洲| 国产黄色视频一区二区在线观看| 亚洲丝袜综合中文字幕| 午夜福利在线观看免费完整高清在| 精品国产一区二区久久| 免费观看无遮挡的男女| 人妻少妇偷人精品九色| 亚洲精品一区蜜桃| 久久国产精品大桥未久av| 国产精品偷伦视频观看了| 国产精品 国内视频| 永久网站在线| 亚洲av国产av综合av卡| 欧美日韩视频精品一区| 自拍欧美九色日韩亚洲蝌蚪91| 日韩三级伦理在线观看| 青青草视频在线视频观看| 男女高潮啪啪啪动态图| 国产精品女同一区二区软件| 校园人妻丝袜中文字幕| 色婷婷久久久亚洲欧美| 亚洲av日韩在线播放| 大码成人一级视频| 亚洲情色 制服丝袜| 国产精品99久久99久久久不卡 | 国产精品蜜桃在线观看| 女性生殖器流出的白浆| 国产一区二区在线观看日韩| 免费观看无遮挡的男女| 九九爱精品视频在线观看| 夫妻性生交免费视频一级片| 国产国语露脸激情在线看| 9191精品国产免费久久| 亚洲av在线观看美女高潮| 久久亚洲国产成人精品v| av网站免费在线观看视频| 97在线视频观看| 日韩一区二区三区影片| 狠狠精品人妻久久久久久综合| av在线老鸭窝| 男的添女的下面高潮视频| 99久久中文字幕三级久久日本| 大香蕉97超碰在线| 亚洲成国产人片在线观看| 精品人妻一区二区三区麻豆| 久久久久久久久久人人人人人人| 极品人妻少妇av视频| 男女边吃奶边做爰视频| 一级片免费观看大全| 天天操日日干夜夜撸| 一区二区av电影网| 亚洲av欧美aⅴ国产| 亚洲内射少妇av| 韩国高清视频一区二区三区| av在线老鸭窝| 在线看a的网站| 国产熟女午夜一区二区三区| 亚洲精品乱久久久久久| 亚洲av电影在线观看一区二区三区| 一级黄片播放器| 欧美日韩亚洲高清精品| 人妻人人澡人人爽人人| 狠狠精品人妻久久久久久综合| 少妇被粗大的猛进出69影院 | 国产黄频视频在线观看| 国产精品人妻久久久影院| 国产精品无大码| 免费大片18禁| 久热这里只有精品99| 狠狠婷婷综合久久久久久88av| 国产精品欧美亚洲77777| 最近中文字幕高清免费大全6| 老司机亚洲免费影院| av在线app专区| 成年人免费黄色播放视频| 高清视频免费观看一区二区| 国产精品国产三级专区第一集| 久久久久久伊人网av| 亚洲精品国产av成人精品| 亚洲人与动物交配视频| 成人亚洲精品一区在线观看| 成人无遮挡网站| av有码第一页| 国产视频首页在线观看| 丝瓜视频免费看黄片| 国产成人午夜福利电影在线观看| 久久鲁丝午夜福利片| 大话2 男鬼变身卡| 欧美精品国产亚洲| 一本大道久久a久久精品| 久久久久国产精品人妻一区二区| 丝袜人妻中文字幕| 亚洲欧美日韩另类电影网站| 日韩人妻精品一区2区三区| 日韩 亚洲 欧美在线| 成年女人在线观看亚洲视频| 午夜日本视频在线| 五月伊人婷婷丁香| 亚洲内射少妇av| 97精品久久久久久久久久精品| 亚洲丝袜综合中文字幕| 国产永久视频网站| 最近2019中文字幕mv第一页| 亚洲欧美精品自产自拍| 狂野欧美激情性bbbbbb| 亚洲av中文av极速乱| 黄片播放在线免费| 国产亚洲一区二区精品| 久久久a久久爽久久v久久| 国产一区二区激情短视频 | 国产精品偷伦视频观看了| 国产精品国产三级国产专区5o| 香蕉丝袜av| 91精品国产国语对白视频| a级毛色黄片| 久久久久精品人妻al黑| 蜜臀久久99精品久久宅男| 国产综合精华液| 久久 成人 亚洲| 最近2019中文字幕mv第一页| 久久久久国产网址| 丰满迷人的少妇在线观看| 久久av网站| 爱豆传媒免费全集在线观看| 亚洲丝袜综合中文字幕| 精品久久蜜臀av无| 99热全是精品| 伦理电影大哥的女人| 国产黄色免费在线视频| www.色视频.com| 亚洲四区av| 在线观看免费视频网站a站| 亚洲精品乱码久久久久久按摩| 日日爽夜夜爽网站| 精品一区二区三卡| 国产精品一区二区在线不卡| 黄色一级大片看看| 亚洲av成人精品一二三区| 国产亚洲一区二区精品| 国产成人精品一,二区| 免费高清在线观看日韩| 日韩精品有码人妻一区| 高清不卡的av网站| 精品亚洲成国产av| av在线观看视频网站免费| 国产69精品久久久久777片| 亚洲av日韩在线播放| 亚洲人与动物交配视频| 女性生殖器流出的白浆| 一区二区av电影网| 久久久精品94久久精品| 黑人巨大精品欧美一区二区蜜桃 | 国产一区二区三区av在线| 免费观看a级毛片全部| 99国产精品免费福利视频| 久久99热这里只频精品6学生| 爱豆传媒免费全集在线观看| 国内精品宾馆在线| 老熟女久久久| 汤姆久久久久久久影院中文字幕| 亚洲内射少妇av| 麻豆乱淫一区二区| 另类亚洲欧美激情| 日产精品乱码卡一卡2卡三| 国产有黄有色有爽视频| 亚洲精品国产av成人精品| 日韩欧美一区视频在线观看| 国产极品天堂在线| 欧美日韩视频高清一区二区三区二| 国产精品久久久久成人av| 久久久久久久久久成人| 亚洲美女视频黄频| 内地一区二区视频在线| 久久国产精品大桥未久av| 波多野结衣一区麻豆| 人妻系列 视频| 亚洲经典国产精华液单| 久久99热这里只频精品6学生| 久久久久久人妻| 国产免费现黄频在线看| 男女高潮啪啪啪动态图| 男人爽女人下面视频在线观看| 亚洲经典国产精华液单| 欧美精品av麻豆av| 赤兔流量卡办理| 蜜臀久久99精品久久宅男| 亚洲欧美日韩卡通动漫| 九草在线视频观看| 国产欧美日韩综合在线一区二区| 亚洲精品成人av观看孕妇| 国产成人精品福利久久| 久久精品久久精品一区二区三区| 免费人成在线观看视频色| 久久人人爽人人爽人人片va| 欧美丝袜亚洲另类| 亚洲精品乱码久久久久久按摩| 亚洲 欧美一区二区三区| 欧美激情国产日韩精品一区| 99久久中文字幕三级久久日本| 极品少妇高潮喷水抽搐| 免费观看av网站的网址| 日韩大片免费观看网站| 精品人妻一区二区三区麻豆| 汤姆久久久久久久影院中文字幕| 欧美成人午夜精品| 两性夫妻黄色片 | 青青草视频在线视频观看| 97超碰精品成人国产| 精品国产一区二区久久| 久久久久久久国产电影| 国产女主播在线喷水免费视频网站| 国产精品嫩草影院av在线观看| 一个人免费看片子| 日韩人妻精品一区2区三区| 亚洲国产av新网站| 亚洲国产av新网站| 在线免费观看不下载黄p国产| 久久久久久人妻| 午夜福利网站1000一区二区三区| av不卡在线播放| 午夜免费观看性视频| 一级毛片 在线播放| 丝袜美足系列| 日韩精品免费视频一区二区三区 | 激情视频va一区二区三区| 成年动漫av网址| 久久久久精品久久久久真实原创| 99久久综合免费| 久久久久国产网址| 精品久久久久久电影网| 咕卡用的链子| 久久久久久久久久久久大奶| 晚上一个人看的免费电影| 成人国产麻豆网| 成年女人在线观看亚洲视频| 一区二区三区精品91| 在线 av 中文字幕| 永久网站在线| av在线播放精品| 日本黄色日本黄色录像| 视频在线观看一区二区三区| 在线 av 中文字幕| 亚洲久久久国产精品| √禁漫天堂资源中文www| 热re99久久精品国产66热6| 美女大奶头黄色视频| 久久人人爽av亚洲精品天堂| av有码第一页| 久久久精品94久久精品| 老司机影院毛片| 青春草国产在线视频| 99视频精品全部免费 在线| 一区二区三区乱码不卡18| 国产成人欧美| 亚洲,欧美精品.| 亚洲国产精品专区欧美| 中文字幕制服av| 男的添女的下面高潮视频| www日本在线高清视频| 热99久久久久精品小说推荐| 又大又黄又爽视频免费| 三上悠亚av全集在线观看| 波野结衣二区三区在线| 好男人视频免费观看在线| 蜜桃在线观看..| 9色porny在线观看| 伦精品一区二区三区| 桃花免费在线播放| av黄色大香蕉| 亚洲激情五月婷婷啪啪| 在线 av 中文字幕| 久久久久国产网址| 日本与韩国留学比较| 人妻系列 视频| 久久鲁丝午夜福利片| 精品少妇久久久久久888优播| 久久久久久久亚洲中文字幕| 久久久国产一区二区| av视频免费观看在线观看| 亚洲精品aⅴ在线观看| 国产免费一区二区三区四区乱码| 18禁国产床啪视频网站| 日韩伦理黄色片| 久久久精品区二区三区| 两性夫妻黄色片 | 老司机影院成人| 久久久久久久久久久免费av| 啦啦啦视频在线资源免费观看| 在线观看免费视频网站a站| 黑人猛操日本美女一级片| 国产精品国产av在线观看| www.熟女人妻精品国产 | 99国产精品免费福利视频| 久久99蜜桃精品久久| 香蕉丝袜av| 精品久久久久久电影网| 国产精品99久久99久久久不卡 | 黑人猛操日本美女一级片| 在线免费观看不下载黄p国产| 亚洲丝袜综合中文字幕| 国产精品一国产av| 国产精品久久久久成人av| 免费看光身美女| 成人影院久久| 少妇精品久久久久久久| 精品久久久久久电影网|