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

    小波包節(jié)點分段閾值降噪在水聲監(jiān)聽中的應(yīng)用?

    2019-12-04 07:41:54英123力123志123初士博123劉茂科123
    應(yīng)用聲學(xué) 2019年6期
    關(guān)鍵詞:波包水聲小波

    趙 杰 楊 英123 惠 力123 王 志123 初士博123 劉茂科123

    (1齊魯工業(yè)大學(xué)(山東省科學(xué)院) 濟南 250353)

    (2山東省科學(xué)院海洋儀器儀表研究所 山東省海洋環(huán)境監(jiān)測技術(shù)重點實驗室 青島 266061)

    (3國家海洋監(jiān)測設(shè)備工程技術(shù)研究中心 青島 266061)

    0 引言

    目前,水聲監(jiān)聽在海洋生物調(diào)查研究、艦艇監(jiān)聽監(jiān)測、深海勘探等方面需求量具大且應(yīng)用前景十分廣闊,其監(jiān)聽信號主要為寬帶中低頻信號和窄帶高頻信號,如魚類的摩擦聲頻率一般在100 Hz~8 kHz,魚鰾發(fā)聲多數(shù)在75~100 Hz;艦船噪聲譜主要集中在1 kHz以下;水聲場測量、海洋地震測量、淺地層剖面測量等主要采用寬帶中低頻進行測量,頻率范圍在幾百Hz到十幾kHz;側(cè)掃聲吶及地貌掃描聲吶等采用窄帶高頻信號,頻率范圍在十幾kHz到幾百kHz。在實際的水聲探測、接收和識別過程中,海洋環(huán)境噪聲信息會弱化相關(guān)聲目標信號監(jiān)聽的準確度[1-2],嚴重影響目標的后續(xù)識別與處理,特別是寬帶中低頻目標信號更容易摻雜大量復(fù)雜的噪聲信息。因此,如何將目標信號與噪聲有效分離與提取成為水聲監(jiān)聽研究的重中之重。

    傳統(tǒng)的數(shù)字信號處理及去噪方法多數(shù)是以傅里葉變換為基礎(chǔ)的,如司新新等[3]利用短時傅里葉變換進行時變和非平穩(wěn)信號的處理和分析,姚家駿等[4]在地震波特征量分析中采用短時傅里葉變換,曲麗榮[5]選用合適的短時傅里葉窗口函數(shù)對信號頻率特性分析。以上研究均以傅里葉變換為基礎(chǔ)的傳統(tǒng)數(shù)字信號去噪法屬于全局變換,無法同時表征時頻雙域特性,只能反映信號整體特性,時域無法局部化,難以檢測局部突變信號,去噪局限性較大。小波變換去噪方法具有良好的時頻局部性、多分辨率、去相關(guān)性、快速算法和選基靈活等特點,與傅里葉變換去噪方法相比具有明顯的優(yōu)越性,成為國內(nèi)外許多學(xué)者專家應(yīng)用研究的去噪方法。Kalpana等[6]對水聲通信中提高信噪比的去噪技術(shù)進行研究,利用Gabor小波變換來實現(xiàn)去噪提高信噪比,并與香農(nóng)小波進行了對比研究。Hill等[7]使用雙數(shù)復(fù)小波系數(shù)法實現(xiàn)圖像噪聲處理研究。普利達大學(xué)研究人員Ghribi等[8]提出用小波變換分解的雙線性結(jié)構(gòu),通過雙通道對稱方法實現(xiàn)小波自適應(yīng)去噪聲。周建等[9]利用一種新的小波閾值去噪方法對振動信號進行去噪分析。張振鳳等[10]利用改進閾值隨實際分解層數(shù)而變化的方法進行去噪,通過信噪比和均方差的比較,該方法優(yōu)于傳統(tǒng)固定閾值、stein無偏似然估計閾值、極大極小值的小波去噪方法。曹建華等[11]利用小波分析法對振動信號進行去噪處理,對去噪后信號進行功率譜密度分析。余本富等[12]基于自適應(yīng)分解層數(shù)和閾值的小波去噪算法降低電能質(zhì)量擾動信號中的噪聲。李戰(zhàn)明等[13]對小波分析中4種去噪方法分析比較,通過對比分析采用改進閾值函數(shù)提升小波變化閾值法用于心電信號去噪分析。以上小波變換的研究方法在閾值和閾值函數(shù)選取等方面具有較高的參考價值,但其主要針對頻率變化緩慢和頻帶較窄的非平穩(wěn)振動、心電等生物醫(yī)學(xué)、圖像處理等方面信號,不能同時兼顧低頻和高頻目標,無法實現(xiàn)寬頻帶信號分頻段聚焦處理,因此,對頻段變換范圍較大的水聲監(jiān)聽信號消噪及提取處理適應(yīng)性較差。

    水聲信號屬于非平穩(wěn)、非線性時變信號,受噪聲及混響等干擾的影響,探測獲得的水聲信號往往會存在很多突變的奇異點。小波變換具有多分辨和時頻局部分析的能力,與短時傅里葉去噪法相比,能更好地去噪,小波和小波包降噪處理的核心是閾值和閾值函數(shù)的選取,常用的小波硬閾值和軟閾值函數(shù)達到了一定的去噪效果,但硬閾值處理后的小波系數(shù)在正負閾值處不連續(xù),軟閾值可有效避免間斷連續(xù)性較好,但小波系數(shù)較大時,系數(shù)絕對值與系數(shù)之間存在偏差,易造成高頻信息損失,重構(gòu)信號易產(chǎn)生誤差。小波包分解系數(shù)按照不同的閾值規(guī)則能更有效去除噪聲,然而不同頻段的小波包分解系數(shù)對目標信號和噪聲信號反映不同,用同一閾值和閾值規(guī)則對整個頻段進行處理無法達到較為理想的效果。因此,本文采用小波包相對能量法進行最優(yōu)分解層判別,考慮各個節(jié)點頻段有用信息的比例,按照比例大小確定閾值準則及函數(shù),進一步提高分辨率,優(yōu)化高頻信號分析,去噪提取結(jié)果更加精細優(yōu)化。利用小波包節(jié)點系數(shù)分段閾值處理的方法對監(jiān)聽獲得的100 Hz~50 kHz范圍內(nèi)水下聲音數(shù)據(jù)降噪處理分析,對不同頻段范圍內(nèi)的水聲信號進行提取,噪聲信號分離,結(jié)果表明:該方法既可以對1 kHz上下范圍內(nèi)的較強信號進行降噪提取,也可對較高頻率相對較弱的信號實現(xiàn)降噪分離。

    1 小波包算法基本原理

    小波包算法為最佳子帶樹結(jié)構(gòu)(Optimal subband tree structuring),包括分解算法和重構(gòu)算法。用子帶樹來表示小波包分解,利用多次迭代的小波轉(zhuǎn)換對信號的低頻和高頻部分進行多層次劃分,根據(jù)信號特征,自適應(yīng)選擇頻帶范圍,分析信號細節(jié)部分。小波包克服了小波分析的高頻分辨率低的缺點,對信號在全頻帶范圍內(nèi)進行正交分解,提高時-頻分辨率[14]。子帶樹結(jié)構(gòu)圖如圖1所示。

    圖1 子帶樹結(jié)構(gòu)Fig.1 Subband tree structure

    將尺度空間Vj和小波子空間Wj統(tǒng)一用U(n,j)新的子空間來表示:

    定義子空間U(n,j)為函數(shù)Un(t)的閉包空間,Un(t)是函數(shù)U2n(t)的閉包空間,令Un(t)滿足雙尺度方程:

    式(2)系數(shù)為多分辨率分析中的濾波器系數(shù),滿足正交關(guān)系,g(k)=(-1)kh(1-k)。當(dāng)n=0,

    在多分辨率分析中,U0(t)和U1(t)分別退化為尺度函數(shù)φ(t)和小波基函數(shù)ψ(t):

    小波包變換的分解與合成中,存在由多種小波基組成的基底庫,因此需要根據(jù)不同場合、不同信號并且在綜合考慮小波基的正交性、緊支性、對稱性等特性的基礎(chǔ)上選擇合適的小波基。從小波庫中尋找使代價函數(shù)最小的小波包基,選擇關(guān)鍵是代價函數(shù)[15]。

    式(5)中,s代表信號,si表示信號在小波包基上的分解系數(shù)。求出使代價函數(shù)最小的小波包序列,從而求出最優(yōu)基。代價函數(shù)稱為熵或平均信息量,代價函數(shù)自變量是信號在某個基向量組上分解系數(shù),函數(shù)實部因變量即為代價。采用二叉樹自下向上搜索最優(yōu)小波包基,小波包基的選取是后續(xù)數(shù)據(jù)處理的關(guān)鍵。

    1.1 小波包相對能量確定最優(yōu)分解層數(shù)

    小波包分解需要獲得一個最優(yōu)的分解層數(shù)。分解層數(shù)直接影響信號的去噪效果,分解層數(shù)過多對各層小波系數(shù)閾值處理過程中造成信息量丟失嚴重、信噪比下降、運算量增加。分解層數(shù)過少,信噪比提高不多,去噪效果不明顯[16]。原始信號和噪聲信號的小波系數(shù)在不同層次有不同特性,原始信號的小波系數(shù)會隨著分解層次增加而增大,噪聲信號小波系數(shù)減小。

    通過最優(yōu)小波包中節(jié)點能量判斷的方法確定小波包分解層數(shù)。根據(jù)最優(yōu)小波基節(jié)點能量的大小,判斷該節(jié)點有用信號能量大小。能量越高,表示該頻段的有效成分越多,反之,噪聲成分越多。據(jù)此,根據(jù)能量值判斷有效成分,確定最優(yōu)分解層數(shù)。

    如圖1所示,U(i,j)表示第j層的第i個節(jié)點,每個節(jié)點系數(shù)表示某一頻段內(nèi)的信號特征,其中j=0,1,···,n,i=0,1,···,2j-1。能量計算公式為[17]

    式(6)中,s(i,j)為U(i,j)節(jié)點的小波包系數(shù)。由此可得出每層分解后子帶的相對能量公式為

    其中,k(i,j)為節(jié)點i在第j層的相對能量,代表有效成分的含量,對每層的節(jié)點相對能量k(i,j)由大到小進行排序,k(i,j)的最大值記為C(1,j),最小值記為C(n,j),即C(1,j)≥C(2,j)≥···≥C(n,j),其中j代表層數(shù)。首先對信號第一層分解相對能量進行初判:C(1,1)≥0.9、C(2,1)<0.1信號能量頻段相對集中的情況,這種情況下對每層的k(2,j)進行對比,當(dāng)?shù)趈層的C(2,j)大于第j-1層C(2,j-1)時,將j-1作為最優(yōu)分解層數(shù)。

    其他情況下,對信號逐層分解,將C(1,j)首次小于0.2的層數(shù)j進行記錄,如果j>1,最優(yōu)分層數(shù)為m=(j-1);如果j=1,最優(yōu)分解層數(shù)為1。

    1.2 小波包分段閾值處理法

    水聲信號的頻率范圍比較廣,不能簡單地將低頻段作為有用信息、高頻段作為噪聲信號處理,按小波相對能量法實現(xiàn)分段閾值降噪處理,該方法與單一閾值處理相比充分考慮信號頻率段與噪聲分布情況。針對不同的分解系數(shù)選取最合適的閾值函數(shù)和閾值th,最大程度去除噪聲,保留有效成分[18]。目前,常用的閾值估計方法主要有以下幾種:

    (1)固定閾值估計法(sqtwolog),在正態(tài)高斯噪聲模型下,針對多維獨立正態(tài)變量分布得出的結(jié)論。

    (2)無偏似然閾值估計法(rigrsure),小波分解系數(shù)的平方記為λ,從小到大的順序排列,即λ1≤λ2≤···≤λn。設(shè)風(fēng)險向量R,其元素為ri,

    風(fēng)險向量R的最小值rmin作為風(fēng)險值,對應(yīng)求出對應(yīng)的λmin,則閾值為

    (3)啟發(fā)式閾值法(heursure),是前兩種閾值的綜合,設(shè)S為N個小波系數(shù)的平方和。令

    (4)極大極小閾值法(minimaxi),也是固定閾值的一種,產(chǎn)生最小均方誤差的極值。

    以上各式中σ=MAD/0.6745為噪聲標準方差,MAD為各高頻子帶系數(shù)的中值;N為信號長度。

    如1.1節(jié)所述,若C(1,1)≥ 0.9、C(2,1)<0.1信號能量頻段相對集中,只對C(1,1)≥0.9進行小波去噪。其他情況下,信號最優(yōu)分解層m上,節(jié)點相對能量k(i,j)≥0.2的節(jié)點系數(shù)采用無偏似然估計閾值估計方法;k(i,j)≥0.1、k(i,j)<0.2的節(jié)點系數(shù)采用極大極小閾值法;k(i,j)<0.1的節(jié)點系數(shù)采用統(tǒng)一閾值法選擇。多閾值準則選取表如1所示。

    表1 多閾值準則選取表Table1 Multiple threshold criteria selection table

    目前常用的閾值函數(shù)為硬閾值函數(shù)和軟閾值函數(shù)。硬閾值將絕對值不大于閾值的元素置零,存在不連續(xù)現(xiàn)象,數(shù)學(xué)表達式為

    軟閾值函數(shù)在式(13)的基礎(chǔ)上實現(xiàn)連續(xù)點收縮為零,有效避免間斷,連續(xù)性較好,但小波系數(shù)較大時,系數(shù)絕對值與系數(shù)之間存在偏差,易造成高頻信息損失,重構(gòu)信號易產(chǎn)生誤差,其表達式為

    考慮到硬閾值和軟閾值函數(shù)處理信號細節(jié)成分的局限性,利用文獻[19]中的新型閾值函數(shù),公式如下:

    式(15)中,x表示小波包分解后的節(jié)點系數(shù),y為閾值函數(shù)處理后的小波節(jié)點系數(shù),th為選定的閾值。式(14)和式(15)中sign(x)為符號函數(shù),當(dāng)x>0,sign(x)=1;當(dāng)x=0,sign(x)=0;當(dāng)x<0,sign(x)=-1。閾值函數(shù)處理結(jié)果對比如圖2、圖3所示。

    圖2 軟、硬閾值函數(shù)結(jié)果對比Fig.2 Comparison of soft and hard threshold function results

    圖3 新、軟閾值函數(shù)結(jié)果對比Fig.3 Comparison of new and soft threshold function results

    2 監(jiān)聽水聲信號去噪分析

    監(jiān)聽水聲信號來源于阿拉斯加費班克大學(xué)Erin Pettit和華盛頓大學(xué)JeffNystuen等在阿拉斯加州冰灣利用自主水聽器(被動水聽器(PAL))獲得的被動水聲記錄[20],該海域的水深為160 m,水聽器的布放水深為90 m。所有數(shù)據(jù)采樣頻率為100 kHz,從100 Hz~50 kHz共分為64個頻段,其中100 Hz~3 kHz的范圍內(nèi)的頻率分辨率為200 Hz,3 kHz~50 kHz范圍的為分辨率為1 kHz。讀取2009年某一時間段內(nèi)的原始音頻文件,其時頻特性標度圖如圖4所示。

    由圖4的時-頻域分布標度圖可以看出原始水聲監(jiān)聽信號頻率隨時間變化的強弱分布,本文的主要工作是對不同頻段范圍、不同強弱的信號進行消噪并且提取。將最原始水聲監(jiān)聽信號按db6進行小波包分解計算,計算每層節(jié)點相對能量并由大到小排序(每層節(jié)點僅排前四位),如表2所示。根據(jù)每層節(jié)點相對能量的大小按照1.1節(jié)判別方法獲得最優(yōu)分解層數(shù)為5層。按小波包5層將水聲信號分解,進行最優(yōu)樹結(jié)構(gòu)和節(jié)點系數(shù)頻率大小排序,分別如圖5和圖6所示。的方式實現(xiàn)消噪分離,處理后時-頻域標度分布如圖7所示。依此頻段來進行本文方法與全局單一閾值對比,將小波包硬閾值、小波包軟閾值在5層分解的基礎(chǔ)上對信號進行消噪處理,獲得的時-頻域信號分布圖,分別如圖8和圖9所示。由圖7、圖8和圖9三幅圖片對比可以看出,本文方法具有良好的優(yōu)越性。

    表2 節(jié)點相對能量排序表Table2 Nodes relative energy ranking table

    圖4 監(jiān)聽水聲信號時-頻分布標度圖Fig.4 Time-frequency distribution scale diagram for monitoring underwater acoustic signals

    圖5 最優(yōu)樹結(jié)構(gòu)Fig.5 Optimal tree structure

    圖6 節(jié)點系數(shù)頻率排序圖Fig.6 Node coefficient frequency sorting diagram

    圖7 分段閾值降噪信號時-頻分布標度圖Fig.7 Time-frequency distribution scale diagram of segmental threshold denoising signal

    圖6中橫軸為時間,左邊縱軸為頻率大小順序,右邊縱軸為小波包分解第五層節(jié)點系數(shù)31-62按照頻率順序由下往上進行排序。由圖6可看出監(jiān)聽水聲信號頻段主要分布在31、32、34三個節(jié)點上,頻率范圍在0~5 kHz。通過本文方法,首先對頻段在節(jié)點31處監(jiān)聽水聲信號,頻率范圍大約在0.1 kHz~1.6 kHz,按照節(jié)點系數(shù)分頻段閾值處理

    圖8 硬閾值降噪信號時-頻分布標度圖Fig.8 Time-frequency distribution scale diagram of hard threshold denoising signal

    圖9 軟閾值降噪信號時-頻分布標度圖Fig.9 Time-frequency distribution scale diagram of soft threshold denoising signal

    圖7為利用本文方法的消噪提取結(jié)果,信號提取結(jié)果較為理想,噪聲基本消除。圖8中為小波包硬閾值5層分解處理后的監(jiān)聽水聲信號時-頻域標度分布圖,與圖4原始信號相比,信號降噪較為明顯,但依然存在噪聲。圖9為小波包軟閾值降噪處理,與圖8相比,噪聲頻率信號進一步減少,但降噪不夠徹底,仍存在小部分噪聲。小波包分段閾值降噪分離的關(guān)鍵在于在未知真實有效的水聲信號前提下,節(jié)點系數(shù)按頻率排序,節(jié)點系數(shù)頻率由低到高、信號強度由強到弱的排序為31、32、34、33、37,其頻率范圍大約在0.1 kHz~1.6 kHz、1.7 kHz~3.3 kHz、 3.4 kHz~5 kHz、 5.1 kHz~6.7 kHz、6.8 kHz~8.4 kHz,按照本文方法對剩余4個節(jié)點頻率信號進行提取,其消噪提取結(jié)果較為理想,結(jié)果如圖10~13所示。由圖12、圖13可以看出頻率范圍相對較高、信號強度較弱的節(jié)點33、節(jié)點37的信號消噪提取亦可通過本文方法實現(xiàn)。

    由于監(jiān)聽的真實水聲信號未知,參照文獻[21]平滑度指標、信噪比增益和峰值信噪比三種方法對監(jiān)聽水聲信號進行去噪評價。三種方法的定義及公式表示如下:

    (1)平滑度指標是指去噪后的差分數(shù)的方差根和原始信號的差分數(shù)的方差根之比,記為r,

    其中,f′(n)為去噪后函數(shù),f(n)為原始信號。

    (2)信噪比增益是指去噪后的信噪比和去噪前的信噪比的比值,記為GSNR,

    其中,SNRdn為去噪后信噪比,SNRn為原始信號信噪比。

    (3)峰值信噪比信號最大可能功率和影響它的表示精度的破壞性噪聲功率的比值,記為PSNR,

    其中,f′(n)為去噪后函數(shù),f(n)為原始信號。

    圖10 節(jié)點32頻段降噪信號時-頻分布標度圖Fig.10 Node 32 time-frequency distribution scale diagram of segmental threshold denoising signal

    圖11 節(jié)點34頻段降噪信號時-頻分布標度圖Fig.11 Node 34 time-frequency distribution scale diagram of segmental threshold denoising signal

    圖12 節(jié)點33頻段降噪信號時-頻分布標度圖Fig.12 Node 33 time-frequency distribution scale diagram of segmental threshold denoising signal

    圖13 節(jié)點37降噪信號時-頻分布標度圖Fig.13 Node 37 time-frequency distribution scale diagram of segmental threshold denoising signal

    平滑度指標數(shù)值越小,信噪比增益和峰值信噪比越大,去噪效果越好。表3是采用上述三種方法分別對小波包硬閾值,小波包軟閾值和小波包分段閾值的降噪處理結(jié)果進行去噪評價,由表3可以看出本文方法是可靠的。

    表3 評價結(jié)果Table3 The evaluation results

    3 結(jié)論

    通過對水下監(jiān)聽信號的去噪分析可以得出本文方法的關(guān)鍵在于(1)按照每層小波包節(jié)點系數(shù)相對能量的變化確定小波包最優(yōu)分解層數(shù);(2)按相對能量法實現(xiàn)分頻段閾值降噪處理。根據(jù)節(jié)點相對能量大小,判斷該節(jié)點有用信號能量大小。能量越高,表示該頻段的有效成分越多,反之,噪聲成分越多。針對節(jié)點系數(shù)相對能量大小選取最合適的閾值函數(shù)和閾值th,最大程度去除噪聲,保留有效成分。本文方法去噪效果優(yōu)于傳統(tǒng)小波、小波包全局單一閾值處理等方法,可在信號噪聲的處理中進行廣泛的推廣和應(yīng)用,特別在100 Hz~50 kHz頻率范圍內(nèi)的水聲監(jiān)聽信號降噪處理中具有重要參考價值。

    猜你喜歡
    波包水聲小波
    構(gòu)造Daubechies小波的一些注記
    基于MATLAB的小波降噪研究
    電子制作(2019年13期)2020-01-14 03:15:32
    基于小波包Tsallis熵和RVM的模擬電路故障診斷
    基于改進的G-SVS LMS 與冗余提升小波的滾動軸承故障診斷
    認知水聲通信系統(tǒng)中OFDM技術(shù)的應(yīng)用
    電子制作(2017年22期)2017-02-02 07:10:34
    新型多功能水聲應(yīng)答器電子系統(tǒng)設(shè)計
    電子制作(2017年19期)2017-02-02 07:08:28
    FRFT在水聲信道時延頻移聯(lián)合估計中的應(yīng)用
    基于小波包變換的電力系統(tǒng)諧波分析
    小波包理論與圖像小波包分解
    基于壓縮感知的水聲數(shù)據(jù)壓縮與重構(gòu)技術(shù)
    狂野欧美激情性xxxx在线观看| 如何舔出高潮| 国产麻豆成人av免费视频| 午夜a级毛片| 在线天堂最新版资源| 少妇的逼水好多| 欧美成人免费av一区二区三区| 欧美成人a在线观看| 禁无遮挡网站| 久久6这里有精品| 免费看a级黄色片| 久久国产乱子免费精品| 久久精品人妻少妇| 国产午夜精品论理片| 久久久久精品国产欧美久久久| 亚洲七黄色美女视频| 精品人妻偷拍中文字幕| 中出人妻视频一区二区| 狂野欧美白嫩少妇大欣赏| 色综合亚洲欧美另类图片| 天天躁夜夜躁狠狠久久av| 国产精品国产高清国产av| 婷婷亚洲欧美| 亚洲人成网站在线观看播放| 中文字幕精品亚洲无线码一区| 日本免费a在线| 久久久午夜欧美精品| 精品久久久久久久人妻蜜臀av| 国产精品一二三区在线看| 91午夜精品亚洲一区二区三区| 国产 一区精品| 国产 一区精品| 免费观看在线日韩| 午夜福利在线观看免费完整高清在 | 免费在线观看影片大全网站| 天堂动漫精品| 3wmmmm亚洲av在线观看| 亚洲成人久久性| 午夜福利高清视频| 在线a可以看的网站| 国产高清不卡午夜福利| 国产亚洲精品久久久久久毛片| 男人舔奶头视频| 99久久九九国产精品国产免费| 啦啦啦观看免费观看视频高清| 国产aⅴ精品一区二区三区波| 女的被弄到高潮叫床怎么办| 欧美高清性xxxxhd video| 淫妇啪啪啪对白视频| 国内精品久久久久精免费| 国产毛片a区久久久久| 此物有八面人人有两片| 亚洲成人精品中文字幕电影| 99视频精品全部免费 在线| 变态另类成人亚洲欧美熟女| 黄色日韩在线| a级一级毛片免费在线观看| 午夜精品国产一区二区电影 | 一级黄色大片毛片| 两性午夜刺激爽爽歪歪视频在线观看| 国产视频一区二区在线看| av在线蜜桃| 亚洲国产精品久久男人天堂| 日韩欧美国产在线观看| 男女那种视频在线观看| 亚洲最大成人中文| 身体一侧抽搐| 国产亚洲精品久久久久久毛片| 国产白丝娇喘喷水9色精品| 波野结衣二区三区在线| 1000部很黄的大片| 91午夜精品亚洲一区二区三区| 久久久午夜欧美精品| 国产男靠女视频免费网站| 欧美三级亚洲精品| 身体一侧抽搐| 人人妻,人人澡人人爽秒播| 国产男靠女视频免费网站| 男女下面进入的视频免费午夜| 久久婷婷人人爽人人干人人爱| 日韩欧美一区二区三区在线观看| 久久久久久伊人网av| 舔av片在线| 一级毛片电影观看 | 男女视频在线观看网站免费| av视频在线观看入口| 在线观看av片永久免费下载| 免费av观看视频| 日韩成人伦理影院| 一个人看的www免费观看视频| 亚洲性久久影院| 亚洲无线观看免费| 国产免费男女视频| 别揉我奶头 嗯啊视频| 久久精品国产鲁丝片午夜精品| 秋霞在线观看毛片| 久久草成人影院| 日韩欧美国产在线观看| 亚洲国产精品国产精品| 欧美成人精品欧美一级黄| 亚洲第一电影网av| 日韩欧美国产在线观看| 日本一本二区三区精品| 精品久久久久久成人av| 成熟少妇高潮喷水视频| 久久久国产成人免费| 婷婷六月久久综合丁香| 一级a爱片免费观看的视频| 淫妇啪啪啪对白视频| 午夜激情福利司机影院| 亚洲欧美日韩无卡精品| 村上凉子中文字幕在线| h日本视频在线播放| 91狼人影院| 午夜老司机福利剧场| 嫩草影院新地址| 免费搜索国产男女视频| 欧美成人一区二区免费高清观看| 色尼玛亚洲综合影院| 日韩,欧美,国产一区二区三区 | 1000部很黄的大片| 亚州av有码| 最近手机中文字幕大全| 小蜜桃在线观看免费完整版高清| 亚洲美女黄片视频| 国产精品久久久久久av不卡| 免费观看精品视频网站| av天堂在线播放| 亚洲av一区综合| 成人鲁丝片一二三区免费| 免费看光身美女| 久久精品国产自在天天线| 免费人成视频x8x8入口观看| 大又大粗又爽又黄少妇毛片口| 91久久精品电影网| 欧美性猛交╳xxx乱大交人| 免费看av在线观看网站| h日本视频在线播放| 亚洲av美国av| 久久精品国产亚洲网站| 亚洲经典国产精华液单| 亚洲精品色激情综合| 一级毛片我不卡| 一卡2卡三卡四卡精品乱码亚洲| 精品福利观看| 久久精品影院6| 午夜激情欧美在线| 中文字幕av在线有码专区| 欧美又色又爽又黄视频| 一区二区三区免费毛片| 看片在线看免费视频| 男女啪啪激烈高潮av片| 91在线精品国自产拍蜜月| 晚上一个人看的免费电影| 久久久久久国产a免费观看| 免费在线观看成人毛片| 伦理电影大哥的女人| 一进一出好大好爽视频| 少妇人妻精品综合一区二区 | 亚洲色图av天堂| 香蕉av资源在线| 99久国产av精品国产电影| 日韩一区二区视频免费看| 久久热精品热| 日韩av在线大香蕉| 大又大粗又爽又黄少妇毛片口| 此物有八面人人有两片| 18禁在线无遮挡免费观看视频 | 国产精品精品国产色婷婷| 一级av片app| 免费av不卡在线播放| 美女黄网站色视频| 最新在线观看一区二区三区| 欧美三级亚洲精品| 精品免费久久久久久久清纯| 精品久久久噜噜| 蜜臀久久99精品久久宅男| 亚洲av第一区精品v没综合| 干丝袜人妻中文字幕| 三级经典国产精品| 久久久精品94久久精品| 国产日本99.免费观看| 不卡视频在线观看欧美| 亚洲av成人av| 性欧美人与动物交配| 黄色配什么色好看| 99热只有精品国产| 观看美女的网站| 久久久久久久久久久丰满| 亚洲中文字幕一区二区三区有码在线看| 久久精品综合一区二区三区| 国产精品永久免费网站| 人人妻,人人澡人人爽秒播| 国产精品久久久久久av不卡| 亚洲三级黄色毛片| 久久精品国产自在天天线| 97热精品久久久久久| 国产综合懂色| 国产男靠女视频免费网站| 女人十人毛片免费观看3o分钟| 国产 一区 欧美 日韩| 久久精品人妻少妇| 一级毛片我不卡| 亚洲欧美成人精品一区二区| 99久久精品热视频| 久久久久性生活片| 精品熟女少妇av免费看| 欧美在线一区亚洲| 国产精品人妻久久久久久| 国产伦精品一区二区三区视频9| 最近手机中文字幕大全| 久久久久国产精品人妻aⅴ院| 国产精品一及| 最新在线观看一区二区三区| 成人av一区二区三区在线看| 成人无遮挡网站| 久久久欧美国产精品| 国产高清激情床上av| 春色校园在线视频观看| 亚洲性久久影院| 国产私拍福利视频在线观看| 亚洲欧美精品综合久久99| 不卡一级毛片| 丝袜美腿在线中文| 国产精品不卡视频一区二区| 天堂动漫精品| 女的被弄到高潮叫床怎么办| 中文字幕av成人在线电影| 一级毛片久久久久久久久女| 欧美最黄视频在线播放免费| www.色视频.com| 久久精品夜色国产| 欧美三级亚洲精品| 国产成人a∨麻豆精品| 啦啦啦韩国在线观看视频| 久久国产乱子免费精品| 一进一出抽搐动态| 国产高清不卡午夜福利| 狠狠狠狠99中文字幕| 成人一区二区视频在线观看| 久久精品国产亚洲av香蕉五月| 最近手机中文字幕大全| 在线播放无遮挡| 国产精品av视频在线免费观看| 成人av在线播放网站| 大香蕉久久网| 黄色视频,在线免费观看| 99久久成人亚洲精品观看| 亚洲精品一区av在线观看| 91在线精品国自产拍蜜月| 99精品在免费线老司机午夜| 亚洲综合色惰| 色av中文字幕| 国产三级在线视频| 精品人妻熟女av久视频| 日本黄色视频三级网站网址| 免费搜索国产男女视频| 久久精品国产亚洲av涩爱 | 久久久色成人| 亚洲专区国产一区二区| 狠狠狠狠99中文字幕| 日日摸夜夜添夜夜添小说| 国产伦精品一区二区三区视频9| 非洲黑人性xxxx精品又粗又长| .国产精品久久| 噜噜噜噜噜久久久久久91| 国产综合懂色| 日韩欧美一区二区三区在线观看| 如何舔出高潮| 亚洲欧美日韩无卡精品| 麻豆乱淫一区二区| 国内少妇人妻偷人精品xxx网站| 免费人成视频x8x8入口观看| 不卡视频在线观看欧美| 欧美性猛交╳xxx乱大交人| 成人美女网站在线观看视频| 日韩制服骚丝袜av| 99久国产av精品| 中文字幕av成人在线电影| 插阴视频在线观看视频| 在线观看66精品国产| 欧美一区二区国产精品久久精品| 亚洲av成人精品一区久久| 91av网一区二区| 1024手机看黄色片| 亚洲国产高清在线一区二区三| 日本爱情动作片www.在线观看 | 色哟哟·www| 好男人在线观看高清免费视频| 国产成人freesex在线 | 黄色配什么色好看| 校园人妻丝袜中文字幕| 亚洲无线观看免费| 九九热线精品视视频播放| 中文在线观看免费www的网站| 搡老岳熟女国产| 老司机影院成人| 亚洲精品一卡2卡三卡4卡5卡| 我要看日韩黄色一级片| a级一级毛片免费在线观看| 99热网站在线观看| 别揉我奶头 嗯啊视频| 久久人人精品亚洲av| 国产成人aa在线观看| 插逼视频在线观看| 国产成人a∨麻豆精品| 久久久久久久久久成人| 激情 狠狠 欧美| 精品国内亚洲2022精品成人| 亚洲av免费在线观看| 直男gayav资源| 国产 一区精品| 51国产日韩欧美| 欧美日本亚洲视频在线播放| 少妇人妻精品综合一区二区 | 国产精品福利在线免费观看| 麻豆乱淫一区二区| 蜜臀久久99精品久久宅男| 精品人妻熟女av久视频| 国产爱豆传媒在线观看| 观看美女的网站| 亚洲国产精品国产精品| 永久网站在线| 国产精品三级大全| 午夜福利视频1000在线观看| 最近在线观看免费完整版| 亚洲av电影不卡..在线观看| 亚洲国产精品合色在线| 男女做爰动态图高潮gif福利片| 一级毛片我不卡| 天堂网av新在线| 一级黄色大片毛片| 99热只有精品国产| 一级黄片播放器| 日本熟妇午夜| 久久精品国产自在天天线| 午夜老司机福利剧场| av专区在线播放| 最好的美女福利视频网| 卡戴珊不雅视频在线播放| 亚洲久久久久久中文字幕| 国产熟女欧美一区二区| 国产黄a三级三级三级人| 麻豆精品久久久久久蜜桃| 成人永久免费在线观看视频| 亚洲精品成人久久久久久| 成人av在线播放网站| 寂寞人妻少妇视频99o| 一边摸一边抽搐一进一小说| 日本成人三级电影网站| 成人无遮挡网站| 国产精品1区2区在线观看.| 69人妻影院| 亚洲性久久影院| 黄色视频,在线免费观看| 在线国产一区二区在线| 最近视频中文字幕2019在线8| 啦啦啦观看免费观看视频高清| 亚洲av美国av| 啦啦啦观看免费观看视频高清| 亚洲av美国av| 大又大粗又爽又黄少妇毛片口| 深夜精品福利| 天堂网av新在线| 丰满人妻一区二区三区视频av| 欧美色欧美亚洲另类二区| 插逼视频在线观看| 在现免费观看毛片| 国产一级毛片七仙女欲春2| 18禁裸乳无遮挡免费网站照片| 99久久九九国产精品国产免费| 校园人妻丝袜中文字幕| 国产精品国产高清国产av| 99国产精品一区二区蜜桃av| 老熟妇乱子伦视频在线观看| 亚洲七黄色美女视频| 国产精品乱码一区二三区的特点| 天美传媒精品一区二区| 不卡视频在线观看欧美| 一区二区三区四区激情视频 | 深夜a级毛片| 女人被狂操c到高潮| 色在线成人网| 极品教师在线视频| 亚洲欧美成人综合另类久久久 | 免费看日本二区| 一级a爱片免费观看的视频| 99久国产av精品| 欧美bdsm另类| 欧美国产日韩亚洲一区| 国产成人91sexporn| 天堂影院成人在线观看| 天堂动漫精品| 少妇的逼水好多| 九九热线精品视视频播放| 日韩在线高清观看一区二区三区| 一夜夜www| 日韩,欧美,国产一区二区三区 | 18禁在线播放成人免费| 日韩精品中文字幕看吧| 18禁在线无遮挡免费观看视频 | 99热6这里只有精品| 欧美一区二区精品小视频在线| 国产精品一区www在线观看| 午夜福利高清视频| 在线看三级毛片| 日本与韩国留学比较| 男人狂女人下面高潮的视频| 国产精品av视频在线免费观看| 级片在线观看| 少妇的逼好多水| 日日摸夜夜添夜夜添小说| 淫秽高清视频在线观看| 欧美人与善性xxx| 久久6这里有精品| 国产精品国产高清国产av| 亚洲高清免费不卡视频| 国产熟女欧美一区二区| 久99久视频精品免费| 日韩人妻高清精品专区| 在线免费观看不下载黄p国产| 麻豆久久精品国产亚洲av| av女优亚洲男人天堂| 亚洲激情五月婷婷啪啪| 亚洲一区高清亚洲精品| 男女下面进入的视频免费午夜| 国产av一区在线观看免费| 亚洲精品色激情综合| av专区在线播放| 舔av片在线| 不卡视频在线观看欧美| 在线观看午夜福利视频| 国产极品精品免费视频能看的| 乱系列少妇在线播放| 久久精品人妻少妇| 别揉我奶头~嗯~啊~动态视频| 欧美区成人在线视频| 亚洲欧美日韩卡通动漫| 中文字幕av在线有码专区| 99久久中文字幕三级久久日本| 一卡2卡三卡四卡精品乱码亚洲| 国内少妇人妻偷人精品xxx网站| 亚洲欧美清纯卡通| 免费无遮挡裸体视频| 欧美性感艳星| 久久久a久久爽久久v久久| 久久精品久久久久久噜噜老黄 | 日本与韩国留学比较| 国产精品一区www在线观看| 精品不卡国产一区二区三区| 一夜夜www| 亚洲国产日韩欧美精品在线观看| 国产成年人精品一区二区| 九九爱精品视频在线观看| 国产精品无大码| 51国产日韩欧美| 麻豆久久精品国产亚洲av| 日韩欧美精品免费久久| 一级毛片aaaaaa免费看小| 久久久久久伊人网av| 久久草成人影院| 色在线成人网| 精品无人区乱码1区二区| 寂寞人妻少妇视频99o| 午夜免费激情av| 免费观看在线日韩| 欧美xxxx性猛交bbbb| 最好的美女福利视频网| 久久精品国产鲁丝片午夜精品| 欧美日韩精品成人综合77777| 免费人成视频x8x8入口观看| 国产精品99久久久久久久久| 综合色av麻豆| 亚洲欧美日韩高清在线视频| 久久精品国产亚洲av香蕉五月| 日本五十路高清| 最近在线观看免费完整版| 亚洲精品一卡2卡三卡4卡5卡| 精品免费久久久久久久清纯| 久久亚洲国产成人精品v| 国产一区二区亚洲精品在线观看| 久久久久久久亚洲中文字幕| 婷婷色综合大香蕉| 成熟少妇高潮喷水视频| 丰满人妻一区二区三区视频av| 赤兔流量卡办理| 深爱激情五月婷婷| 免费av不卡在线播放| 麻豆成人午夜福利视频| 国产精品福利在线免费观看| 精品久久久久久久久久久久久| 亚洲av成人av| 国产免费男女视频| 国模一区二区三区四区视频| 日韩一本色道免费dvd| 免费看av在线观看网站| 丰满的人妻完整版| 免费人成在线观看视频色| a级毛片a级免费在线| 久久鲁丝午夜福利片| 久久久久国内视频| 97碰自拍视频| 国产精品久久久久久久久免| 天堂动漫精品| 亚洲成av人片在线播放无| 少妇猛男粗大的猛烈进出视频 | 亚洲自拍偷在线| 国产av一区在线观看免费| 久久久久久国产a免费观看| 亚洲最大成人中文| 国产亚洲av嫩草精品影院| 国产爱豆传媒在线观看| 国产91av在线免费观看| 久久久久久久久久黄片| 国产亚洲精品av在线| 午夜亚洲福利在线播放| 中文亚洲av片在线观看爽| av女优亚洲男人天堂| 亚洲精品成人久久久久久| 老师上课跳d突然被开到最大视频| 欧美性猛交黑人性爽| av在线亚洲专区| 俺也久久电影网| 97碰自拍视频| 日本在线视频免费播放| 99久久精品热视频| 久久人妻av系列| 色综合亚洲欧美另类图片| av黄色大香蕉| 菩萨蛮人人尽说江南好唐韦庄 | 午夜福利在线观看免费完整高清在 | 日本成人三级电影网站| 毛片一级片免费看久久久久| 精品福利观看| 一级毛片我不卡| 国产视频一区二区在线看| 国产成人影院久久av| av.在线天堂| 97热精品久久久久久| 国产老妇女一区| 亚洲最大成人手机在线| 国产精品亚洲一级av第二区| 久久天躁狠狠躁夜夜2o2o| 国产毛片a区久久久久| 最近中文字幕高清免费大全6| 亚洲中文字幕日韩| 桃色一区二区三区在线观看| 国产精品av视频在线免费观看| 赤兔流量卡办理| 少妇猛男粗大的猛烈进出视频 | 国产男人的电影天堂91| 欧美激情国产日韩精品一区| 99久国产av精品国产电影| 日韩欧美精品v在线| 十八禁国产超污无遮挡网站| 乱人视频在线观看| 在线观看66精品国产| 久久精品国产亚洲网站| 99热只有精品国产| 欧美激情国产日韩精品一区| 欧美色视频一区免费| 久久国产乱子免费精品| 免费人成在线观看视频色| 久久6这里有精品| 精品人妻偷拍中文字幕| 蜜桃久久精品国产亚洲av| 久久久久久大精品| 免费av不卡在线播放| 在线免费观看不下载黄p国产| 免费av毛片视频| 久久久久久久久中文| 六月丁香七月| av天堂在线播放| av在线老鸭窝| 欧美一区二区亚洲| 全区人妻精品视频| 久久久国产成人精品二区| 亚洲av美国av| 久久久久国产网址| 简卡轻食公司| 亚洲欧美成人综合另类久久久 | 天堂网av新在线| 91在线精品国自产拍蜜月| 国产高潮美女av| 亚洲国产精品国产精品| 听说在线观看完整版免费高清| 中文亚洲av片在线观看爽| 中出人妻视频一区二区| 国产精品三级大全| 欧美区成人在线视频| 国产日本99.免费观看| 亚洲18禁久久av| 欧美激情久久久久久爽电影| 在线免费观看不下载黄p国产| 成人三级黄色视频| 精品久久久久久成人av| 高清日韩中文字幕在线| 青春草视频在线免费观看| 色哟哟·www| 欧美人与善性xxx| 两个人视频免费观看高清| 亚洲人成网站在线播放欧美日韩| 国产欧美日韩一区二区精品| 又粗又爽又猛毛片免费看| 国产视频内射| 久久国产乱子免费精品| 尤物成人国产欧美一区二区三区| 深夜精品福利| av专区在线播放| 69人妻影院| 日本成人三级电影网站| 能在线免费观看的黄片| 国产午夜精品久久久久久一区二区三区 | 亚洲一区高清亚洲精品| 精品国产三级普通话版|