亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        低空目標(biāo)被動(dòng)聲識別關(guān)鍵技術(shù)研究*

        2017-04-24 02:30:03
        關(guān)鍵詞:信號

        雷 鳴 喬 柯

        (西安工業(yè)大學(xué)電子信息工程學(xué)院 西安 710021)

        低空目標(biāo)被動(dòng)聲識別關(guān)鍵技術(shù)研究*

        雷 鳴 喬 柯

        (西安工業(yè)大學(xué)電子信息工程學(xué)院 西安 710021)

        低空目標(biāo)被動(dòng)聲探測技術(shù)在現(xiàn)代戰(zhàn)爭中發(fā)揮著越來越重要的作用。針對低空目標(biāo)被動(dòng)聲識別的一些關(guān)鍵技術(shù)進(jìn)行了研究。首先從低空目標(biāo)聲信號與語音信號特征的相似性出發(fā),論證了低空目標(biāo)被動(dòng)聲識別與語音識別可以采用同樣的方法。然后對四旋翼飛行器的聲音信號進(jìn)行了時(shí)頻分析,提出基于MFCC參數(shù)的低空目標(biāo)識別方法。最后論述了低空目標(biāo)聲音識別技術(shù)中的關(guān)鍵技術(shù),重點(diǎn)分析了特征提取中的MFCC技術(shù)以及使用動(dòng)態(tài)時(shí)間規(guī)整(DTW)算法進(jìn)行仿真識別實(shí)驗(yàn),為低空目標(biāo)聲音識別系統(tǒng)的開發(fā)奠定了基礎(chǔ)。

        被動(dòng)聲識別; 端點(diǎn)檢測; 美爾倒譜參數(shù); 動(dòng)態(tài)時(shí)間規(guī)整

        1 引言

        被動(dòng)聲識別也稱為被動(dòng)式聲雷達(dá),與傳統(tǒng)雷達(dá)探測技術(shù)相比,有著抗干擾、低功耗、不易被發(fā)現(xiàn)等優(yōu)點(diǎn),可以彌補(bǔ)雷達(dá)低空探測存在盲區(qū)這一不足。被動(dòng)聲目標(biāo)識別是指聲傳感器接實(shí)時(shí)接收聲目標(biāo)信息,利用包含在聲目標(biāo)信號中特有的特征信息,與典型的聲音信號(如坦克、輪式車輛、直升機(jī)等)建立模式匹配,通過計(jì)算機(jī)自動(dòng)識別聲目標(biāo)。被動(dòng)聲目標(biāo)識別系統(tǒng)可以作為雷達(dá)的輔助系統(tǒng),在未來的戰(zhàn)場上有著非常重要的應(yīng)用前景[1]。

        聲紋識別技術(shù)在過去的幾十年飛速發(fā)展,取得了令人鼓舞的成果。本文在分析低空目標(biāo)聲信號與語音信號相似性的基礎(chǔ)上,借鑒語音信號特征參數(shù)的提取方法,將MFCC參數(shù)用于低空目標(biāo)識別。由于實(shí)時(shí)環(huán)境復(fù)雜,存在噪聲干擾,所以提出MFCC參數(shù)的改進(jìn)算法,提高系統(tǒng)對噪聲的魯棒性和識別成功率。

        2 低空目標(biāo)聲信號與語音信號的相似性

        低空目標(biāo)聲信號和語音信號都是一種媒質(zhì)振動(dòng)形式,兩者在聲學(xué)特征方面具有相似性,體現(xiàn)在以下幾個(gè)方面[2]:

        1) 清音。語音中,由于聲道收縮,使氣流以高速通過這一收縮部分,沒有聲帶振動(dòng)產(chǎn)生。在聲目標(biāo)中,語音清音理解為機(jī)械性噪聲。機(jī)械性噪聲主要由傳動(dòng)系統(tǒng)的軸承部件、曲柄活塞部件、鏈傳動(dòng)部件、液壓系統(tǒng)部件等多種運(yùn)動(dòng)部件的來回重復(fù)運(yùn)動(dòng)所致。

        2) 濁音。語音中,當(dāng)氣流通過聲門時(shí),聲帶的張力使聲帶產(chǎn)生張弛振蕩式振動(dòng),產(chǎn)生準(zhǔn)周期脈沖氣流。聲目標(biāo)中,聲信號主要考慮空氣動(dòng)力性噪聲和排氣噪聲??諝鈩?dòng)力性噪聲主要由發(fā)動(dòng)機(jī)的進(jìn)氣、排氣噪聲以及風(fēng)扇噪聲等組成。

        3) 聲道模型。語音中,把人的語音聲道視為由多個(gè)不同截面積的管子串聯(lián)而成的系統(tǒng)加以研究。聲目標(biāo)中,發(fā)動(dòng)機(jī)罩(無論前置還是后置發(fā)動(dòng)機(jī))至車身底部之間的空間,相當(dāng)于一個(gè)諧振腔或?yàn)V波器。同樣也可以把戰(zhàn)場聲聲道作為非均勻截面管加以研究。

        4) 基音頻率。語音中,發(fā)濁音時(shí),氣流通過聲門時(shí)使聲帶發(fā)生振動(dòng),產(chǎn)生準(zhǔn)周期激勵(lì)脈沖串,這個(gè)脈沖串周期的倒數(shù)就是/基音頻率0。對聲目標(biāo)而言,發(fā)聲濁音時(shí),進(jìn)氣道中的空氣由于激勵(lì)發(fā)生振動(dòng),產(chǎn)生準(zhǔn)周期激勵(lì)脈沖串,這個(gè)脈沖串周期的倒數(shù)就是/基音頻率0。

        5) 共振峰。語音中,人的聲道和鼻道都是非均勻的聲道管,將其視為諧振腔,共振峰就是這個(gè)腔體的諧振頻率。對聲目標(biāo)而言,諧振腔或?yàn)V波器的諧振頻率也可稱為共振峰,因聲目標(biāo)類型的不同而不同。

        戰(zhàn)場聲信號與語音的形成過程以及信號特征的比較結(jié)果表明[3]:

        1) 低空目標(biāo)聲信號和語音信號均可用幅度特征、短時(shí)過零率、聲道共振峰來描述;

        2) 低空目標(biāo)聲信號與語音都是媒質(zhì)振動(dòng)模式。戰(zhàn)場聲音的聲源位置分散,語音聲源則比較集中;

        3) 低空目標(biāo)聲信號特征提取和識別與語音特征提取和識別可以采用同樣的方法。

        3 低空目標(biāo)聲信號時(shí)頻分析

        3.1 時(shí)域特征參數(shù)

        聲音信號的時(shí)域分析方法是最簡單,最直觀的方法。它直接分析聲音信號的時(shí)域波形,提取出語音的特征參數(shù),一般用于最基本的信號處理中,在很多語音信號處理中,都會(huì)用到短時(shí)平均能量,短時(shí)平均過零率,短時(shí)平均幅度,短時(shí)自相關(guān)函數(shù)和短時(shí)平均幅度差函數(shù)等時(shí)域的參數(shù)組合[4~5]。這里簡單介紹一下短時(shí)平均能量,短時(shí)平均過零率。本文研究分析對象為四旋翼飛行器的聲音信號。

        1) 短時(shí)平均能量

        對于語音信號{x(n)}短時(shí)平均能量的定義如下

        (1)

        上式中h(n)=w2(n),En代表第n點(diǎn)的短時(shí)能量,短時(shí)能量En是由語音信號先經(jīng)過平方運(yùn)算之后,再把得到的信號進(jìn)行線性濾波之后輸出。

        h(n)表示線性濾波器的沖擊響應(yīng)。短時(shí)能量的計(jì)算過程如圖1所示。

        圖1 短時(shí)平均能量框圖

        2) 短時(shí)平均過零率

        短時(shí)過零率法:短時(shí)過零率表示一幀聲音信號波形穿過橫軸(零電平)的次數(shù)。對于離散信號序列,如果相鄰的采樣值改變符號則稱為過零。即樣本改變符號的次數(shù),聲音信號x(n)其公式如下

        (2)

        其中,sgn[x(n)]函數(shù)是符號函數(shù)。

        采用雙門限兩級判決的聲音斷點(diǎn)檢測方法,來確定聲音信號的起點(diǎn)和結(jié)束點(diǎn)位置,處理時(shí)間越少,越能抑制噪聲干擾,還可以提高信號的質(zhì)量,這個(gè)時(shí)候的端點(diǎn)檢測最有效。做端點(diǎn)檢測前,設(shè)置閾值,對于能量和過零率都包含一個(gè)低門限和一個(gè)高門限。低門限數(shù)值小,對于信號變化敏感,可能被其他信號超過,高門限數(shù)值大,不容易被超過。

        實(shí)現(xiàn)能量與過零率的端點(diǎn)檢測步驟如下:

        1) 分幀

        在端點(diǎn)檢測的具體運(yùn)行中,首先對聲信號進(jìn)行分幀,在分幀的基礎(chǔ)上求出短時(shí)平均能量和短時(shí)過零率,然后逐幀地依據(jù)閾值進(jìn)行比較和判斷。聲音信號x(n)進(jìn)行分幀處理,每一幀記為si(n),n=1,2,…N,n為離散聲音信號的時(shí)間序列,N為幀長,i表示幀數(shù)。

        2) 計(jì)算每一幀的短時(shí)能量,得到聲音信號的短時(shí)幀能量:

        (3)

        3) 計(jì)算每一幀語音的過零率,得到短時(shí)過零率

        (4)

        4) 第一級判決

        (1)根據(jù)在聲音短時(shí)能量包絡(luò)線上選取一個(gè)較高閾值(門限)T2進(jìn)行一次粗判,就是高于該閾值肯定是聲音信號,聲音信號起止點(diǎn)應(yīng)該位于該閾值與短時(shí)能量包絡(luò)線交點(diǎn)所對應(yīng)的時(shí)間點(diǎn)之外。

        (2)在平均能量上確定一個(gè)較低的閾值(門限)T1,并從第一級判決與包絡(luò)線的兩個(gè)交點(diǎn)分別向左右搜索,分別找到短時(shí)能量包絡(luò)與閾值T1相交的兩個(gè)點(diǎn),于是,這兩個(gè)交點(diǎn)便是雙門限法根據(jù)短時(shí)能量所判定的聲音信號的起止點(diǎn)位置。

        5) 第二級判決

        以短時(shí)平均過零率為準(zhǔn),以短時(shí)能量最終判定的兩個(gè)交點(diǎn)分別向左,向右搜索,找到短時(shí)過零率低于某個(gè)閾值(門限)T3的兩點(diǎn),這便是該聲音信號的起止點(diǎn)。端點(diǎn)檢測如圖2所示。

        圖2 端點(diǎn)檢測圖

        3.2 頻域特征參數(shù)

        Mel頻率倒譜系數(shù)(Mel Frequency Cepstrum Coefficient)的分析是基于人的聽覺機(jī)理,依據(jù)人的聽覺實(shí)驗(yàn)結(jié)果來分析語音的頻譜,期望能獲得好的語音特性。MFCC分析依據(jù)的聽覺機(jī)理有兩個(gè)[6]。

        1) 人的主觀感知頻域的劃定并不是線性的,根據(jù)Stevens和Volkman的工作有下面的公式:

        Fmel=1125log(1+f/700)

        (5)

        Fmel是以美爾(Mel)為單位的感知頻率,f是以Hz為單位的實(shí)際頻率。

        2) 頻率群相應(yīng)于人耳基底膜分成許多很小的部分,每一部分對應(yīng)一個(gè)頻率群,對于同一頻率群的那些頻率的聲音,在大腦中疊加在一起進(jìn)行評價(jià)。按臨界帶的劃分,將語音在頻域上劃分成一系列的頻率群組成的濾波器組,即Mel濾波器組。取漢明窗的Mel濾波器組的響應(yīng)曲線如圖3所示。

        圖3 漢明窗的頻率響應(yīng)曲線

        MFCC特征參數(shù)提取[7~8]原理框圖如圖4所示。

        圖4 MFCC原理框圖

        1) 預(yù)處理

        預(yù)處理包括預(yù)加重、分幀、加窗函數(shù)。

        · 預(yù)加重:聲音信號的頻率響應(yīng)曲線接近于一個(gè)二階低通濾波器,預(yù)加重的目的是為了補(bǔ)償高頻分量的損失,提升高頻分量。預(yù)加重的濾波器常設(shè)為

        H(z)=1-az-1

        (6)

        式中,a為常數(shù)。

        · 分幀處理:由于聲音信號是一個(gè)準(zhǔn)穩(wěn)態(tài)的信號,把它分成較短的幀,在每幀中可將其看作穩(wěn)態(tài)信號,可以用處理穩(wěn)態(tài)信號的方法來處理。同時(shí),為了使一幀與另一幀之間的參數(shù)能較平穩(wěn)的過渡,在相鄰兩幀之間互相有部分重疊。

        · 加窗函數(shù):加窗函數(shù)的目的是減少頻域中的泄露,將對每一幀語音乘以海寧窗。聲音信號x(n)經(jīng)預(yù)處理后為xi(m),其中下標(biāo)i表示分幀后的第i幀。

        2) 快速傅里葉變換

        對每一幀信號進(jìn)行FFT變換,從時(shí)域數(shù)據(jù)轉(zhuǎn)變?yōu)轭l域數(shù)據(jù):

        X(i,k)=FFT[xi(m)]

        (7)

        3) 計(jì)算譜線能量

        對每一幀F(xiàn)FT后的數(shù)據(jù)計(jì)算譜線的能量:

        E(i,k)=[X(i,k)]2

        (8)

        4) 計(jì)算通過MEL濾波器的能量

        把求出的每幀譜線能量譜通過MEL濾波器,計(jì)算在MEL濾波器中的能量。在頻域中相當(dāng)于把每幀的能量譜E(i,k)(其中i表示第i幀,k表示頻域中的第k條譜線)與MEL濾波器的MEL濾波器的頻率響應(yīng)Hm(K)相乘并相加:

        (9)

        5) 計(jì)算DCT倒譜

        (10)

        k=0,1,…,N-1

        (11)

        其中,參數(shù)N是序列x(n)的長度;c(k)是正交因子。

        飛行器X1的聲信號MFCC如圖5所示。

        圖5 飛行器聲信號MFCC參數(shù)圖

        4 對傳統(tǒng)算法的改進(jìn)

        聲音信號特征的Mel倒譜系數(shù)(MFCC)模擬了人耳聽覺感知特性:不同頻率的聲音,會(huì)耳內(nèi)基礎(chǔ)膜的不同位置振動(dòng)。這樣人耳就可以很容易分辨出各種聲音,所以此參數(shù)作為目標(biāo)的識別參數(shù)能更好地反應(yīng)各種聲信號的特性。但標(biāo)準(zhǔn)的MFCC只反映了語音參數(shù)的靜態(tài)特性,而人耳對語音的動(dòng)態(tài)特性更為敏感,本系統(tǒng)特色是引入一階差分(MFCC),二階差分(MFCC)再加上時(shí)域特征信息的幀能量參數(shù)構(gòu)成了(3Q+1)維的特征矢量,更好地消除了語音幀之間的相關(guān)性,更優(yōu)的逼近語音的動(dòng)態(tài)特征,大大提高了識辨率[9]。通過計(jì)算語音動(dòng)態(tài)差分倒譜,k取常數(shù),通常取2:

        (12)

        用以上公式求得的參數(shù)是當(dāng)前幀的前兩幀和后兩幀的線性組合,即是一階差分MFCC差分倒譜參數(shù),把結(jié)果再代入上式就可以得到二階參數(shù)。時(shí)域特征矢量用短時(shí)(幀)平均能量與短時(shí)平均過零率相結(jié)合的兩級雙門限端點(diǎn)檢測法表征,來判斷聲音信號的起始點(diǎn)與終止點(diǎn)。

        5 實(shí)驗(yàn)仿真及分析

        5.1 實(shí)驗(yàn)條件

        實(shí)驗(yàn)中采用實(shí)測噪聲環(huán)境下三種不同型號的四旋翼飛行器1、2、3的聲信號(采樣頻率分別為50kHz、10kHz、10kHz)作為訓(xùn)練、識別數(shù)據(jù)。訓(xùn)練時(shí)分別取不同時(shí)間獲得的聲音信號各10段作為學(xué)習(xí)樣本,每個(gè)樣本3000個(gè)點(diǎn),即樣本分別取得時(shí)間段為:60ms、300ms、300ms。分幀時(shí)幀長為256個(gè)點(diǎn),幀移為80個(gè)點(diǎn);24階的MFCC參數(shù)提取時(shí)數(shù)字濾波器組選取24個(gè),對所有聲音段進(jìn)行預(yù)加重(系數(shù)0.97)加窗(漢明窗)處理,取幀256點(diǎn)幀移10ms分別提取MFCC,MFCC+ΔMFCC,MFCC+ΔMFCC+ΔΔMFCC+Mn特征參數(shù),識別模型采用動(dòng)態(tài)時(shí)間規(guī)整(DTW)方法,利用Matlab進(jìn)行識別仿真[10]。被動(dòng)聲識別框圖如圖6所示。

        圖6 低空目標(biāo)被動(dòng)聲識別框圖

        5.2 實(shí)驗(yàn)結(jié)果分析

        圖7 系統(tǒng)識別率

        由圖7可知:MFCC+ΔMFCC+ΔΔMFCC+Mn識別率為96.63%,MFCC+ΔMFCC識別率為93.45%,MFCC識別率為91.87%,識別率MFCC+ΔMFCC+ΔΔMFCC+Mn明顯最高。

        此方法缺點(diǎn)在于計(jì)算量較大,但隨著計(jì)算機(jī)技術(shù)的飛速發(fā)展,為大計(jì)算量問題的解決提供了較多簡潔的方法,大計(jì)算量已經(jīng)不成問題。以犧牲計(jì)算量和計(jì)算時(shí)間為代價(jià),明顯提高聲音信號的識別率,在聲紋識別系統(tǒng)中已得到了廣泛的應(yīng)用。

        6 結(jié)語

        文中借鑒語音信號的特征提取方法,提出將MFCC參數(shù)應(yīng)用于低空被動(dòng)聲目標(biāo)識別中。考慮到實(shí)時(shí)條件下存在強(qiáng)噪聲干擾,提出一種改進(jìn)的MFCC特征參數(shù)提取方法,通過實(shí)驗(yàn)比較得知:改進(jìn)算法后MFCC+ΔMFCC+ΔΔMFCC+Mn(3Q+1)特征提取的識別率高達(dá)96.6%,識別率得到明顯改善。此法從很大程度上擬合了人耳對語音處理的動(dòng)態(tài)特點(diǎn),而且具有一定的控噪性,在同等條件下大大提高了低空目標(biāo)的聲音信號識別率,從而實(shí)現(xiàn)聲目標(biāo)信號的有效分析,提高了強(qiáng)噪聲條件下聲目標(biāo)識別正確率。其在聲紋識別領(lǐng)域有廣闊的應(yīng)用前景。

        [1] 陳功.戰(zhàn)場被動(dòng)聲目標(biāo)識別關(guān)鍵技術(shù)的研究[D].南京:解放軍理工大學(xué),2007:1-5. CHEN Gong. Research on Key Technologies of passive acoustic target recognition in battlefield[D]. Nanjing: PLA University of Science and Technology,2007:1-5.

        [2] 易克初,田斌,付強(qiáng).語音信號處理[M].北京:國防工業(yè)出版社,2003:21-25. YI Kechu, TIAN Bin, FU Qiang. Processing of speech signal[M]. Beijing: National Defense Industry Press,2003:21-25.

        [3] 朱志松.戰(zhàn)場聲目標(biāo)特征提取研究[J].探測與控制學(xué)報(bào),2006,28(3):9-11. ZHU Zhisong. Study on the feature extraction of acoustic target in battlefield[J]. Journal of Detection and Control,2006,28(3):9-11.

        [4] 夏輝達(dá).基于DSP的戰(zhàn)場聲目標(biāo)識別技術(shù)的研究[D].太原:中北大學(xué)碩士學(xué)位論文,2004:32-36. XIA Huida. Research on the battlefield acoustic target recognition technology based on DSP[D]. Taiyuan: Master Thesis of North Central University,2004:32-36.

        [5] 許可喜.被動(dòng)聲探測若干關(guān)鍵技術(shù)研究[D].南京:南京理工大學(xué)碩士論文,2005:24-52. XU Kexi. Research on several key technologies of passive acoustic detection[D]. Nanjing: Master Thesis of Nanjing University of Science and Technology,2005:24-52.

        [6] 張萬里,劉橋.Mel頻率倒譜系數(shù)提取及其在聲紋識別中的作用[J].貴州大學(xué)學(xué)報(bào),2005,(2):207-230. ZHANG Wanli, LIU Qiao. Mel frequency cepstral coefficients are extracted and the role of voiceprint recognition[J]. Journal of Guizhou University,2005,(2):207-230.

        [7] 呂國云,許學(xué)忠,趙銳.戰(zhàn)場目標(biāo)被動(dòng)噪聲識別技術(shù)[J].探測與控制學(xué)報(bào),2001,23(4):30-32. LV Guoyun, XU Xuezhong, ZHAO Rui. The passive noise identification technology of the battlefield target[J]. Journal of Detection and Control,2001,23(4):30-32.

        [8] 于勝民.多語言語音識別技術(shù)研究[D].北京:中科院自動(dòng)化所博士學(xué)位論文,2005:45-53. YU Shengmin. Multi language speech recognition technology research[D]. Beijing: Chinese Academy of Sciences, Institute of automation, doctoral dissertation,2005:45-53.

        [9] 趙力.語音信號處理[M].北京:機(jī)械工業(yè)出版社,2010:31-56. ZHAO Li. Speech signal processing M[M]. Beijing: Mechanical Industry Press,2010:31-56.

        [10] 楊陽,陳永明.聲紋識別技術(shù)及其應(yīng)用[J].電聲技術(shù),2007,31(2):45-50. YANG Yang, CHEN Yongming. Application of voiceprint recognition technology and[J]. Audio Technology,2007,31(2):45-50.

        Key Technologies of Low Altitude Target Passive Acoustic Recognition

        LEI Ming QIAO Ke

        (Electronic Information Engineering College, Xi’an Technological University, Xi’an 710021)

        Low altitude target passive acoustic detection technology is playing a more and more important role in modern warfare. Some key technologies of passive acoustic recognition for low altitude targets are studied. Firstly, from the similarity of low altitude target sound signal and the characteristics of speech signal, it is proved that the low altitude target passive acoustic recognition and speech recognition can be used in the same way. Then, the time-frequency analysis of the acoustic signals of the four rotor aircraft is carried out, and a low altitude target recognition method based on MFCC parameters is proposed. Finally, the key technologies of low altitude target sound recognition technology is discussed, the MFCC in feature extraction is analyzed emphatically and simulation identification experiment is carried out using dynamic time warping (DTW) algorithm, laid the foundation for the development of low altitude target sound recognition system.

        passive voice recognition, endpoint detection, MFCC, dynamic time warping Class Number TN911

        2016年10月10日,

        2016年11月21日

        雷鳴,男,碩士,副教授,研究方向:測控技術(shù)與通信技術(shù)。喬柯,男,碩士研究生,研究方向:信號處理,通信與信息系統(tǒng)。

        TN911

        10.3969/j.issn.1672-9722.2017.04.011

        猜你喜歡
        信號
        信號
        鴨綠江(2021年35期)2021-04-19 12:24:18
        完形填空二則
        7個(gè)信號,警惕寶寶要感冒
        媽媽寶寶(2019年10期)2019-10-26 02:45:34
        孩子停止長個(gè)的信號
        《鐵道通信信號》訂閱單
        基于FPGA的多功能信號發(fā)生器的設(shè)計(jì)
        電子制作(2018年11期)2018-08-04 03:25:42
        基于Arduino的聯(lián)鎖信號控制接口研究
        《鐵道通信信號》訂閱單
        基于LabVIEW的力加載信號采集與PID控制
        Kisspeptin/GPR54信號通路促使性早熟形成的作用觀察
        又爽又黄又无遮挡网站动态图| 亚洲国产成人无码电影| 国产免费av片在线观看| 色男色女午夜福利影院| 精品 无码 国产观看| 少妇激情一区二区三区视频| 国产一区二区三区亚洲| 日韩AV无码一区二区三区不卡毛片| 亚洲精品中文幕一区二区| 韩国免费一级a一片在线| 女人色毛片女人色毛片18| 国产毛片av最新视频| 琪琪av一区二区三区| 少妇bbwbbw高潮| 国产在线手机视频| 色老头久久综合网老妇女| 久久夜色精品国产亚洲av老牛| 国产精品一区二区黑丝| 少妇性饥渴无码a区免费| 人妻无码一区二区| 日韩精品夜色二区91久久久| 国产区女主播一区在线| 挺进朋友人妻雪白的身体韩国电影| 亚洲一级毛片免费在线观看| 日本视频一区二区二区| 久久精品色福利熟妇丰满人妻91| 人妻激情另类乱人伦人妻| 久久精品无码一区二区乱片子| 青青草视频在线播放81| 欧美大片va欧美在线播放| 午夜一区欧美二区高清三区| 亚洲地区一区二区三区| 激情五月六月婷婷俺来也| 中国精品18videosex性中国| 一本色综合亚洲精品蜜桃冫| 亚洲女同精品久久女同| 国产农村妇女精品一区| 青草视频在线播放| 久久精品国产88久久综合| 丝袜美足在线视频国产在线看| 丰满熟女人妻一区二区三区|