• <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)熵的心音信號檢測方法研究
    保持信噪比的相位分解反褶積方法研究
    久热爱精品视频在线9| 美女国产高潮福利片在线看| 久久久久久人人人人人| 中文字幕人妻丝袜一区二区| 久久久国产成人免费| 在线播放国产精品三级| 国产免费男女视频| 2021天堂中文幕一二区在线观 | 精品欧美一区二区三区在线| 亚洲国产日韩欧美精品在线观看 | 99在线视频只有这里精品首页| 高清在线国产一区| 在线观看免费视频日本深夜| 久久婷婷人人爽人人干人人爱| 黄频高清免费视频| 欧美成人免费av一区二区三区| 韩国av一区二区三区四区| 91字幕亚洲| 国产精品亚洲一级av第二区| 国产免费av片在线观看野外av| 性色av乱码一区二区三区2| 99热6这里只有精品| 日日摸夜夜添夜夜添小说| 国产午夜福利久久久久久| 一区福利在线观看| 欧美zozozo另类| 久久精品国产99精品国产亚洲性色| 亚洲av中文字字幕乱码综合 | 性欧美人与动物交配| 日本一本二区三区精品| 成人特级黄色片久久久久久久| 久久精品国产亚洲av香蕉五月| 日日干狠狠操夜夜爽| 欧美黄色淫秽网站| 成人永久免费在线观看视频| 国产一区在线观看成人免费| 亚洲五月天丁香| 国产亚洲欧美精品永久| www日本在线高清视频| 久久人人精品亚洲av| 国产精品爽爽va在线观看网站 | 亚洲 欧美 日韩 在线 免费| 天堂√8在线中文| 成人手机av| 久久国产精品人妻蜜桃| 国产成人欧美| 国产精品自产拍在线观看55亚洲| 日本撒尿小便嘘嘘汇集6| 真人做人爱边吃奶动态| av中文乱码字幕在线| 亚洲人成77777在线视频| 国产精品美女特级片免费视频播放器 | 精品电影一区二区在线| 亚洲狠狠婷婷综合久久图片| 桃色一区二区三区在线观看| 色老头精品视频在线观看| 欧美日本亚洲视频在线播放| 日韩av在线大香蕉| 中文字幕最新亚洲高清| 亚洲人成网站在线播放欧美日韩| 亚洲第一青青草原| 亚洲成av人片免费观看| 免费在线观看日本一区| 久久欧美精品欧美久久欧美| 精品免费久久久久久久清纯| 国产人伦9x9x在线观看| 一进一出抽搐gif免费好疼| 哪里可以看免费的av片| 在线看三级毛片| 99久久99久久久精品蜜桃| www.自偷自拍.com| 俄罗斯特黄特色一大片| 日本在线视频免费播放| 亚洲成人久久性| 国产三级黄色录像| 亚洲成人久久爱视频| 亚洲精品av麻豆狂野| 欧美一区二区精品小视频在线| 亚洲av成人av| 女人被狂操c到高潮| 国产精品亚洲美女久久久| 波多野结衣av一区二区av| 欧美激情高清一区二区三区| 精品久久蜜臀av无| 又紧又爽又黄一区二区| 久久人妻福利社区极品人妻图片| 女人被狂操c到高潮| 天天一区二区日本电影三级| 国产av一区二区精品久久| 亚洲五月天丁香| 欧美大码av| 欧美三级亚洲精品| 亚洲成av片中文字幕在线观看| 两性夫妻黄色片| 国内精品久久久久精免费| 国产成人欧美| 国产精品爽爽va在线观看网站 | 国产高清视频在线播放一区| 精品国产乱码久久久久久男人| 成人国语在线视频| 亚洲五月色婷婷综合| 日本精品一区二区三区蜜桃| 国内精品久久久久久久电影| 曰老女人黄片| 99精品欧美一区二区三区四区| 真人一进一出gif抽搐免费| 精品国内亚洲2022精品成人| 波多野结衣高清无吗| 俄罗斯特黄特色一大片| 午夜成年电影在线免费观看| 身体一侧抽搐| 国内毛片毛片毛片毛片毛片| 国产人伦9x9x在线观看| 嫩草影视91久久| av在线天堂中文字幕| 久久精品国产99精品国产亚洲性色| 成人国产综合亚洲| 亚洲久久久国产精品| 国产麻豆成人av免费视频| 国产亚洲精品综合一区在线观看 | 日本免费a在线| 91麻豆av在线| 999久久久精品免费观看国产| 老汉色av国产亚洲站长工具| 在线av久久热| 欧美精品亚洲一区二区| 国产一区二区在线av高清观看| 精品一区二区三区四区五区乱码| 岛国在线观看网站| 色综合亚洲欧美另类图片| 给我免费播放毛片高清在线观看| 中文资源天堂在线| 免费av毛片视频| 欧美乱妇无乱码| 成人欧美大片| 男女下面进入的视频免费午夜 | 欧美精品啪啪一区二区三区| 一进一出好大好爽视频| aaaaa片日本免费| 色婷婷久久久亚洲欧美| www.999成人在线观看| 午夜激情福利司机影院| 日韩国内少妇激情av| 欧美黑人精品巨大| 操出白浆在线播放| 手机成人av网站| 免费电影在线观看免费观看| 很黄的视频免费| 日本熟妇午夜| 19禁男女啪啪无遮挡网站| 欧美丝袜亚洲另类 | 可以免费在线观看a视频的电影网站| 手机成人av网站| 99热只有精品国产| 久久久久久亚洲精品国产蜜桃av| 一边摸一边做爽爽视频免费| 极品教师在线免费播放| xxx96com| 成人三级黄色视频| 香蕉国产在线看| 欧美日本视频| 一区福利在线观看| 好男人在线观看高清免费视频 | 欧美成人性av电影在线观看| 精品国产乱子伦一区二区三区| 欧美+亚洲+日韩+国产| 亚洲av第一区精品v没综合| 亚洲欧洲精品一区二区精品久久久| 麻豆久久精品国产亚洲av| 99国产综合亚洲精品| 精品久久久久久成人av| 精品一区二区三区视频在线观看免费| 国产亚洲精品久久久久5区| 成人18禁在线播放| av中文乱码字幕在线| av超薄肉色丝袜交足视频| 国产在线精品亚洲第一网站| 日韩 欧美 亚洲 中文字幕| 在线观看免费日韩欧美大片| 国产精品乱码一区二三区的特点| 欧美激情极品国产一区二区三区| 日本三级黄在线观看| 色哟哟哟哟哟哟| 国产黄a三级三级三级人| 久久久久国产一级毛片高清牌| 精品不卡国产一区二区三区| 国产精品av久久久久免费| 美女扒开内裤让男人捅视频| av天堂在线播放| 国产精品 国内视频| 中文字幕久久专区| 国产高清激情床上av| 久久久国产成人免费| 国产熟女xx| 男人舔女人下体高潮全视频| 国产麻豆成人av免费视频| 白带黄色成豆腐渣| 叶爱在线成人免费视频播放| 国内精品久久久久精免费| 亚洲自拍偷在线| 日韩大码丰满熟妇| 超碰成人久久| 国产成年人精品一区二区| 女人被狂操c到高潮| 在线观看66精品国产| 国产av不卡久久| 狂野欧美激情性xxxx| 日本黄色视频三级网站网址| 51午夜福利影视在线观看| 亚洲人成网站高清观看| 国产真人三级小视频在线观看| 好男人在线观看高清免费视频 | 日韩av在线大香蕉| 久久青草综合色| 国产成人系列免费观看| 久久午夜亚洲精品久久| 久热爱精品视频在线9| 欧美在线一区亚洲| 观看免费一级毛片| 亚洲国产看品久久| 久久这里只有精品19| 一二三四在线观看免费中文在| 琪琪午夜伦伦电影理论片6080| 亚洲avbb在线观看| 女人高潮潮喷娇喘18禁视频| 午夜福利视频1000在线观看| 岛国在线观看网站| 黑丝袜美女国产一区| 国产aⅴ精品一区二区三区波| 国产野战对白在线观看| 亚洲人成电影免费在线| 精品少妇一区二区三区视频日本电影| 亚洲一码二码三码区别大吗| 亚洲一区二区三区不卡视频| 亚洲精品国产精品久久久不卡| 亚洲avbb在线观看| 欧美av亚洲av综合av国产av| 国产成人av教育| 国产精品99久久99久久久不卡| 国产精品日韩av在线免费观看| 欧洲精品卡2卡3卡4卡5卡区| 精品久久久久久,| 国产av又大| 手机成人av网站| 亚洲av电影在线进入| 最近最新中文字幕大全电影3 | 18禁黄网站禁片午夜丰满| 欧美性长视频在线观看| a级毛片a级免费在线| 中文字幕最新亚洲高清| 国产主播在线观看一区二区| 亚洲第一欧美日韩一区二区三区| 国产高清视频在线播放一区| 香蕉国产在线看| 欧美不卡视频在线免费观看 | 亚洲欧美激情综合另类| 少妇粗大呻吟视频| 精品久久久久久久末码| 女性被躁到高潮视频| 淫秽高清视频在线观看| 日韩大尺度精品在线看网址| 在线免费观看的www视频| 久久久久国产一级毛片高清牌| 18禁观看日本| 一级作爱视频免费观看| 身体一侧抽搐| 国产亚洲av嫩草精品影院| www日本在线高清视频| 亚洲无线在线观看| 国产99久久九九免费精品| 国内揄拍国产精品人妻在线 | 国产99白浆流出| 亚洲一码二码三码区别大吗| 欧美乱色亚洲激情| 天天添夜夜摸| 麻豆一二三区av精品| 香蕉av资源在线| 老司机午夜福利在线观看视频| 男人舔女人的私密视频| av福利片在线| 国产一区二区三区视频了| 69av精品久久久久久| 黑人操中国人逼视频| 欧美黑人欧美精品刺激| 亚洲天堂国产精品一区在线| 久久精品国产亚洲av香蕉五月| 在线av久久热| 午夜福利18| 日韩有码中文字幕| 老汉色∧v一级毛片| 夜夜爽天天搞| 国产伦在线观看视频一区| 日韩三级视频一区二区三区| 精品国内亚洲2022精品成人| 久久精品影院6| www.熟女人妻精品国产| 男女下面进入的视频免费午夜 | 国产爱豆传媒在线观看 | 伦理电影免费视频| 草草在线视频免费看| 成年女人毛片免费观看观看9| 91字幕亚洲| 一本一本综合久久| 久久精品91蜜桃| cao死你这个sao货| 亚洲国产精品合色在线| 性色av乱码一区二区三区2| 亚洲av片天天在线观看| 在线天堂中文资源库| 国产成人欧美在线观看| 91字幕亚洲| tocl精华| 亚洲一区二区三区色噜噜| 亚洲精华国产精华精| 又黄又粗又硬又大视频| 欧美激情高清一区二区三区| 一本大道久久a久久精品| 午夜福利免费观看在线| 亚洲av中文字字幕乱码综合 | 国产精品亚洲美女久久久| 国产精品日韩av在线免费观看| 我的亚洲天堂| 天堂√8在线中文| 久久久久久久久免费视频了| 久久天堂一区二区三区四区| 正在播放国产对白刺激| 免费电影在线观看免费观看| 正在播放国产对白刺激| 精品国产美女av久久久久小说| 欧美成人免费av一区二区三区| 亚洲成人精品中文字幕电影| www国产在线视频色| 波多野结衣高清无吗| 久久久国产成人免费| 精品一区二区三区视频在线观看免费| 成年人黄色毛片网站| 亚洲三区欧美一区| 午夜精品在线福利| 国产高清视频在线播放一区| 精品熟女少妇八av免费久了| 国内少妇人妻偷人精品xxx网站 | 亚洲黑人精品在线| 777久久人妻少妇嫩草av网站| av视频在线观看入口| 熟妇人妻久久中文字幕3abv| 男女之事视频高清在线观看| 欧美激情久久久久久爽电影| 色播亚洲综合网| 美女国产高潮福利片在线看| a级毛片在线看网站| 色婷婷久久久亚洲欧美| 免费搜索国产男女视频| 国产乱人伦免费视频| 亚洲欧美一区二区三区黑人| 美女国产高潮福利片在线看| 色精品久久人妻99蜜桃| 又紧又爽又黄一区二区| 在线观看免费日韩欧美大片| 久久婷婷成人综合色麻豆| tocl精华| 日韩免费av在线播放| 国产熟女午夜一区二区三区| 午夜福利高清视频| 国产成人一区二区三区免费视频网站| 久久久久亚洲av毛片大全| 国产真实乱freesex| 国产高清视频在线播放一区| 中文字幕精品亚洲无线码一区 | 国产亚洲av嫩草精品影院| 熟女电影av网| tocl精华| 亚洲 国产 在线| 别揉我奶头~嗯~啊~动态视频| 国产99白浆流出| 午夜精品在线福利| 18禁观看日本| 两人在一起打扑克的视频| 亚洲精品在线观看二区| 国产单亲对白刺激| 亚洲片人在线观看| 级片在线观看| 国产精品二区激情视频| 日日夜夜操网爽| 亚洲熟妇中文字幕五十中出| 夜夜躁狠狠躁天天躁| 亚洲男人天堂网一区| 我的亚洲天堂| 欧美性猛交黑人性爽| 国产1区2区3区精品| 91成年电影在线观看| 亚洲国产日韩欧美精品在线观看 | 久久国产精品影院| 大香蕉久久成人网| 老熟妇仑乱视频hdxx| av福利片在线| 人人妻人人澡欧美一区二区| 国产v大片淫在线免费观看| 精品国产一区二区三区四区第35| 欧美成人免费av一区二区三区| 一区福利在线观看| 久久人妻av系列| 亚洲熟妇中文字幕五十中出| 热re99久久国产66热| 亚洲人成77777在线视频| 久久久久免费精品人妻一区二区 | 一级毛片精品| 亚洲人成77777在线视频| 久久久久免费精品人妻一区二区 | 午夜日韩欧美国产| 国产精品精品国产色婷婷| 熟妇人妻久久中文字幕3abv| 国产精品 国内视频| 欧美精品亚洲一区二区| 成年人黄色毛片网站| ponron亚洲| 色精品久久人妻99蜜桃| 日韩精品中文字幕看吧| 国产精品久久视频播放| 欧美三级亚洲精品| 久久精品成人免费网站| 欧美日韩瑟瑟在线播放| 国产精品 国内视频| 久久香蕉精品热| 香蕉国产在线看| 亚洲 欧美一区二区三区| 男女床上黄色一级片免费看| 成人精品一区二区免费| 51午夜福利影视在线观看| 午夜激情av网站| 亚洲五月婷婷丁香| 999久久久国产精品视频| 精品一区二区三区四区五区乱码| 欧美在线黄色| 久久热在线av| 91成年电影在线观看| 中文字幕精品免费在线观看视频| 中亚洲国语对白在线视频| 欧美三级亚洲精品| x7x7x7水蜜桃| 国产又爽黄色视频| 欧美大码av| 啦啦啦 在线观看视频| 久久精品国产综合久久久| а√天堂www在线а√下载| 婷婷亚洲欧美| 欧美乱妇无乱码| 此物有八面人人有两片| 国产亚洲欧美精品永久| 88av欧美| 91在线观看av| 国产成人精品久久二区二区91| 一级a爱视频在线免费观看| 久9热在线精品视频| 日韩免费av在线播放| 国产97色在线日韩免费| av欧美777| 成人18禁高潮啪啪吃奶动态图| 国产午夜精品久久久久久| 国产亚洲精品av在线| 亚洲一码二码三码区别大吗| 动漫黄色视频在线观看| 午夜老司机福利片| 久久国产精品男人的天堂亚洲| svipshipincom国产片| 波多野结衣高清作品| 男男h啪啪无遮挡| 黄频高清免费视频| 亚洲成人免费电影在线观看| 叶爱在线成人免费视频播放| 黄色成人免费大全| 午夜精品久久久久久毛片777| 亚洲av成人不卡在线观看播放网| 久久久久国产一级毛片高清牌| 一进一出好大好爽视频| 一区二区三区国产精品乱码| 久久久久久久精品吃奶| 精品国产超薄肉色丝袜足j| 美女 人体艺术 gogo| 丁香欧美五月| 久久午夜综合久久蜜桃| 亚洲av日韩精品久久久久久密| 欧美国产精品va在线观看不卡| 中国美女看黄片| 国产爱豆传媒在线观看 | 国产精品乱码一区二三区的特点| 在线观看免费日韩欧美大片| 女人高潮潮喷娇喘18禁视频| 黄色丝袜av网址大全| 国产av不卡久久| 熟妇人妻久久中文字幕3abv| 一区二区三区精品91| 日韩欧美国产一区二区入口| 亚洲国产精品久久男人天堂| 亚洲精品美女久久av网站| 无限看片的www在线观看| 男人操女人黄网站| 2021天堂中文幕一二区在线观 | 国产精品九九99| 搞女人的毛片| 大香蕉久久成人网| 午夜影院日韩av| 国产人伦9x9x在线观看| 国产一区二区三区在线臀色熟女| 国产精品1区2区在线观看.| 国产真实乱freesex| 丁香欧美五月| 高潮久久久久久久久久久不卡| 嫩草影院精品99| 黄色成人免费大全| 黄片大片在线免费观看| 2021天堂中文幕一二区在线观 | 视频区欧美日本亚洲| 欧美不卡视频在线免费观看 | 黄色毛片三级朝国网站| av视频在线观看入口| 老汉色∧v一级毛片| 国产单亲对白刺激| 男女床上黄色一级片免费看| 老司机深夜福利视频在线观看| 在线观看免费日韩欧美大片| 久久人人精品亚洲av| 最新美女视频免费是黄的| 嫩草影院精品99| 侵犯人妻中文字幕一二三四区| 美女 人体艺术 gogo| 久久婷婷人人爽人人干人人爱| 99精品在免费线老司机午夜| 叶爱在线成人免费视频播放| 一区二区三区激情视频| 亚洲一区高清亚洲精品| 在线免费观看的www视频| 好看av亚洲va欧美ⅴa在| 麻豆av在线久日| 欧美又色又爽又黄视频| 在线播放国产精品三级| 天天添夜夜摸| 亚洲电影在线观看av| 国产精品亚洲av一区麻豆| 美女高潮喷水抽搐中文字幕| 国产精品亚洲美女久久久| 国产亚洲av嫩草精品影院| 无遮挡黄片免费观看| 人人澡人人妻人| 男女之事视频高清在线观看| 亚洲国产欧美日韩在线播放| 校园春色视频在线观看| 久久久水蜜桃国产精品网| 色婷婷久久久亚洲欧美| 久99久视频精品免费| 国产又色又爽无遮挡免费看| 免费一级毛片在线播放高清视频| av超薄肉色丝袜交足视频| 日韩欧美国产一区二区入口| 制服人妻中文乱码| 久久人妻福利社区极品人妻图片| 欧美日韩黄片免| 在线av久久热| АⅤ资源中文在线天堂| 日韩精品中文字幕看吧| 无人区码免费观看不卡| 香蕉av资源在线| 成人永久免费在线观看视频| 好男人在线观看高清免费视频 | 脱女人内裤的视频| 国产精品久久久久久人妻精品电影| 91字幕亚洲| 正在播放国产对白刺激| 免费观看人在逋| 99久久久亚洲精品蜜臀av| 亚洲avbb在线观看| 一个人免费在线观看的高清视频| 亚洲熟女毛片儿| 中出人妻视频一区二区| 黄色视频不卡| 国产黄片美女视频| 国产视频内射| 亚洲精品色激情综合| 日本熟妇午夜| 一进一出抽搐动态| 精品国产美女av久久久久小说| 国产成+人综合+亚洲专区| 久久午夜亚洲精品久久| av欧美777| 在线观看舔阴道视频| 亚洲成av人片免费观看| 一区福利在线观看| 欧美日本视频| 欧美色视频一区免费| 两个人免费观看高清视频| 女人爽到高潮嗷嗷叫在线视频| 午夜福利视频1000在线观看| 国产又色又爽无遮挡免费看| 久久久久久久久中文| 禁无遮挡网站| 男女下面进入的视频免费午夜 | 久久精品影院6| av免费在线观看网站| 欧美黄色淫秽网站| 黄色视频不卡| 黑人巨大精品欧美一区二区mp4| 久久精品国产99精品国产亚洲性色| 19禁男女啪啪无遮挡网站| 日本成人三级电影网站| 女人被狂操c到高潮| 久久香蕉国产精品| 黑人巨大精品欧美一区二区mp4| 老司机午夜福利在线观看视频| 国产精品日韩av在线免费观看| 亚洲熟女毛片儿| 国产亚洲欧美精品永久| 久久久精品欧美日韩精品| 欧美亚洲日本最大视频资源| 法律面前人人平等表现在哪些方面| 俄罗斯特黄特色一大片| 午夜视频精品福利|