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

    基于雙譜和支持向量機(jī)的小麥碰撞聲分類研究

    2013-07-27 01:25:32張嚴(yán)嚴(yán)
    關(guān)鍵詞:雙譜蟲蛀對(duì)角

    張嚴(yán)嚴(yán),郭 敏

    陜西師范大學(xué) 計(jì)算機(jī)科學(xué)學(xué)院,西安 710062

    基于雙譜和支持向量機(jī)的小麥碰撞聲分類研究

    張嚴(yán)嚴(yán),郭 敏

    陜西師范大學(xué) 計(jì)算機(jī)科學(xué)學(xué)院,西安 710062

    1 引言

    小麥在儲(chǔ)藏期間,易發(fā)生蟲害和霉變,導(dǎo)致小麥質(zhì)量下降,有效地檢測(cè)和防治工作顯得尤為重要。常見的儲(chǔ)糧檢測(cè)方法有糧食探管法、機(jī)器視覺(jué)檢測(cè)法、電導(dǎo)率檢測(cè)法、近紅外反射光譜法、軟X線成像法等[1]。但這些方法存在勞動(dòng)強(qiáng)度大、代價(jià)高等缺陷。因此,具有輕便、簡(jiǎn)單、快捷、廉價(jià)等優(yōu)勢(shì)的聲檢測(cè)技術(shù)日益成為該領(lǐng)域的研究熱點(diǎn)[2]。

    近年來(lái),在聲檢測(cè)領(lǐng)域,碰撞聲技術(shù)因其方便、快捷得到廣泛應(yīng)用。2001年,Pearson等人首次采用聲碰撞技術(shù)實(shí)時(shí)分析開心果碰撞鋼板聲,根據(jù)碰撞聲信號(hào)的不同,分離殼未開和已開的開心果,分類正確率約為97%[3-4]。2007年,Pearson等人利用碰撞聲信號(hào)檢測(cè)小麥蟲蛀粒(Insect Damaged wheat Kernels,IDK),從時(shí)域信號(hào)建模和短時(shí)窗方法等方面處理信號(hào),用神經(jīng)網(wǎng)絡(luò)進(jìn)行特征分類,識(shí)別IDK準(zhǔn)確率為87%[5]。2007年,Ince等人利用碰撞聲信號(hào)檢測(cè)裂殼榛子和小麥IDK,根據(jù)Fisher判別準(zhǔn)則對(duì)自適應(yīng)分段聲信號(hào)提取特征,進(jìn)行分類,對(duì)完好中破損榛子和完好小麥中IDK識(shí)別率均達(dá)90%以上[6]。2008年,Elbatawi等人把碰撞聲分類系統(tǒng)用于土豆塊空心病識(shí)別,根據(jù)完好土豆塊碰撞聲幅度大于空心塊莖、峰值頻率高于空心塊莖的特點(diǎn),提取信號(hào)三部分特征進(jìn)行分類,準(zhǔn)確率約98%[7]。2012年,Khalesi等人采集了兩個(gè)胡桃品種的碰撞聲,采用時(shí)域和基于FFT變換的頻域分析方法提取信號(hào)特征,用多層前饋神經(jīng)網(wǎng)絡(luò)進(jìn)行分類,識(shí)別率達(dá)90%以上[8]。

    本文采用雙譜和支持向量機(jī)(Support Vector Machine,SVM)相結(jié)合的方法,對(duì)小麥完好粒、蟲蛀粒和霉變粒碰撞鋼板聲信號(hào)進(jìn)行識(shí)別研究。首先對(duì)信號(hào)進(jìn)行雙譜估計(jì),提取信號(hào)雙譜峰值和對(duì)角切片譜兩部分特征,再用支持向量機(jī)分類器進(jìn)行分類,得到了較高的識(shí)別率,本研究為麥粒分揀提供了一種新途徑。

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

    2.1 雙譜理論

    高階統(tǒng)計(jì)分析由于在處理非高斯和非線性信號(hào)方面的特殊性質(zhì),逐漸成為一種信號(hào)分析和特征提取的有效技術(shù)[9]。近年來(lái),高階譜在信號(hào)檢測(cè)、參數(shù)估計(jì)、目標(biāo)分類等方面都受到極大關(guān)注,糧食碰撞聲信號(hào)本質(zhì)上是非平穩(wěn)、非高斯信號(hào),因此可利用高階譜進(jìn)行研究分析。高階譜中雙譜的階數(shù)最低,計(jì)算較為簡(jiǎn)單,且包含了高階譜的所有特征,所以本文采用基于雙譜的特征提取方法。

    雙譜可以由信號(hào)的離散傅里葉變換表示:

    式中,X(ωi)為信號(hào)的離散傅里葉變換;ωi為頻率變量;i= 1,2。并且雙譜具有以下性質(zhì):

    本文采用直接法進(jìn)行雙譜估計(jì)。對(duì)于一有限長(zhǎng)度的信號(hào),先估計(jì)其傅里葉值,然后對(duì)該信號(hào)作三重相關(guān)運(yùn)算,即可得到雙譜的估計(jì)。估計(jì)時(shí)采樣頻率設(shè)為 fs,Δ0=fs/L是雙譜區(qū)域內(nèi)沿水平和垂直方向上兩頻率采樣點(diǎn)之間的間隔,采樣點(diǎn)數(shù)為L(zhǎng),步驟如下:

    (1)將數(shù)據(jù)L分成N段,每段M個(gè)數(shù)據(jù),即L=NM,每段減去均值(為便于做FFT算法,可在每段補(bǔ)零)。

    (2)對(duì)每段數(shù)據(jù)作DFT變換,即

    其中,λ=1,2,…,M-1;i=1,2,…,N。

    (3)取M=M1L0,M1為正的奇數(shù),即:M1=2S1+1,適當(dāng)選擇L0使得M為偶數(shù),根據(jù)DFT系數(shù),得到每段數(shù)據(jù)的雙譜:

    2.2 支持向量機(jī)

    支持向量機(jī)是在統(tǒng)計(jì)學(xué)習(xí)理論的基礎(chǔ)上發(fā)展起來(lái)的新一代學(xué)習(xí)算法。其基本思想為:先通過(guò)非線性變換將輸入空間變換到一個(gè)高維空間,即將非線性問(wèn)題轉(zhuǎn)換成高維特征空間中的線性問(wèn)題進(jìn)行分類,然后在這個(gè)新空間中求取最優(yōu)分類面,這種變換是通過(guò)定義適當(dāng)?shù)膬?nèi)積函數(shù)來(lái)實(shí)現(xiàn)的。支持向量機(jī)中所求得的分類函數(shù)在形式上類似于神經(jīng)網(wǎng)絡(luò),輸出是若干中間層節(jié)點(diǎn)的線性組合,每個(gè)中間層節(jié)點(diǎn)對(duì)應(yīng)于輸入樣本與一個(gè)支持向量的內(nèi)積,又被稱為支持向量網(wǎng)絡(luò)。

    支持向量機(jī)關(guān)鍵在于內(nèi)積函數(shù)(核函數(shù))選擇。低維空間向量集映射到高維空間,增加了計(jì)算復(fù)雜度,而核函數(shù)正好巧妙地解決了這個(gè)問(wèn)題,只要選用適當(dāng)?shù)暮撕瘮?shù),就可以得到高維空間的分類函數(shù)。

    3 基于雙譜和支持向量機(jī)的小麥分類識(shí)別

    首先對(duì)信號(hào)進(jìn)行直接雙譜估計(jì),得到雙譜等高線圖,然后從兩方面提取特征:(1)提取雙譜峰值;(2)提取對(duì)角切片譜的譜峰。最后,在支持向量機(jī)分類器中對(duì)樣本進(jìn)行分類,分類過(guò)程分為兩步:樣本訓(xùn)練和樣本測(cè)試。樣本訓(xùn)練用核函數(shù)將樣本映射到高維空間,在該空間中計(jì)算支持向量和最優(yōu)分類面。測(cè)試時(shí),用訓(xùn)練得到的最優(yōu)分類面判別映射到高維空間的測(cè)試樣本,確定測(cè)試樣本分類編號(hào),根據(jù)分類編號(hào)計(jì)算識(shí)別率?;陔p譜和支持向量機(jī)的分類識(shí)別流程如圖1所示。

    圖1 實(shí)驗(yàn)流程圖

    3.1 雙譜分析

    實(shí)驗(yàn)中,取小麥完好粒、蟲害粒和霉變粒各300粒,麥粒掉落在24 cm×11 cm×0.06 cm的鋼板上,掉落高度為50 cm,鋼板的水平傾斜角度為30°,用型號(hào)為舒爾BG 4.1的專業(yè)麥克風(fēng)錄制碰撞聲信號(hào),并傳送至裝有MAYA 44聲卡的計(jì)算機(jī)中,采樣頻率為48 000 Hz,完好粒、IDK和霉變粒各采集100個(gè)聲信號(hào)。由于有效碰撞聲只集中在很小一段樣本點(diǎn)內(nèi),為減少計(jì)算量,提高處理速度,對(duì)信號(hào)雙譜估計(jì)前先進(jìn)行預(yù)處理,提取出有效信號(hào)。

    3.1.1 雙譜估計(jì)

    計(jì)算時(shí),將長(zhǎng)為6 144的樣本點(diǎn)分為128段,經(jīng)估計(jì)得到一個(gè)256×256的雙譜幅度譜Bx( ) ω1,ω2,圖2給出了3類小麥碰撞聲的雙譜等高線及三維圖,不同聲信號(hào)差別顯著。

    雙譜等高線圖中:

    (1)完好粒的幅值分布范圍較廣,幾乎覆蓋低頻及高頻區(qū)域。

    (2)蟲蛀粒幅值在低頻部分沒(méi)有分布,只在高頻區(qū)域疏散分布。

    (3)霉變粒幅值主要集中在低頻部分,圍繞著(0,0)周圍分布。

    雙譜三維圖中:

    (1)完好粒的峰值比較密集,其峰值大小和蟲蛀粒處于一個(gè)數(shù)量級(jí)。

    (2)蟲蛀粒的峰值密度最少,但其峰值是三者中最大的。

    (3)霉變粒的峰值相對(duì)較少,其峰值是三者中最小的。

    這種差別的主要原因與3類麥粒物理結(jié)構(gòu)不同有關(guān)。蟲蛀的麥粒結(jié)構(gòu)會(huì)因?yàn)楹οx的洞穴而被改變,其碰撞鋼板時(shí)會(huì)發(fā)生共振,聲音散列。而發(fā)生霉變的麥粒會(huì)使胚乳多空,麥粒收縮,重量減輕,碰撞鋼板的力度較小,聲音偏低。綜上可見:麥粒的內(nèi)部構(gòu)造發(fā)生變化時(shí),會(huì)造成它們碰撞鋼板的聲音有差別,這種差別在雙譜中表現(xiàn)顯著。

    3.1.2 特征提取

    本文從以下兩方面提取特征:

    (1)提取雙譜峰值

    雙譜的譜峰是反映雙譜結(jié)構(gòu)的一個(gè)重要參數(shù)。在雙譜三維圖中,按1/256的間隔沿Y軸做平行于XZ平面的等間隔截面,可得到雙譜的256個(gè)截面,取各截面最大雙譜值作為特征向量,得到一個(gè)256維的特征向量。圖3為小麥碰撞聲信號(hào)雙譜峰值圖。

    (2)提取對(duì)角切片譜

    對(duì)角切片是指當(dāng)ω1=ω2時(shí),B(ω1,ω2)中對(duì)應(yīng)的雙譜值。雙譜的對(duì)角切片包含有豐富的特征信息,切片譜的譜峰及形態(tài)決定了雙譜的結(jié)構(gòu),同時(shí)也反應(yīng)了信號(hào)的特征。圖4是3類聲信號(hào)的對(duì)角切片圖,差異顯著,因此本文提取了對(duì)角切片譜的256個(gè)譜峰作為特征向量。

    3.2 支持向量機(jī)分類識(shí)別

    3.2.1 最優(yōu)核函數(shù)參數(shù)和懲罰參數(shù)選取

    實(shí)際分類中,支持向量機(jī)性能的優(yōu)劣依賴于核函數(shù)和懲罰參數(shù)的選取。不同核函數(shù)對(duì)應(yīng)于不同算法,會(huì)產(chǎn)生不同的支持向量點(diǎn)和分類面,最終對(duì)分類結(jié)果產(chǎn)生很大影響。

    圖2 3類小麥碰撞聲雙譜等高線圖

    圖3 3類小麥碰撞聲的雙譜峰值圖

    圖4 3類小麥碰撞聲雙譜對(duì)角切片圖

    SVM對(duì)核函數(shù)的要求有兩點(diǎn):(l)必須滿足Mercer條件;(2)在高維空間中可以構(gòu)造一個(gè)支持向量數(shù)相對(duì)較少的最優(yōu)分類界面。目前,較常用的核函數(shù)有3種:徑向基核函數(shù)、多項(xiàng)式核函數(shù)、Sigmoid核函數(shù)。經(jīng)實(shí)驗(yàn)對(duì)比,本文選擇分類正確率較高且計(jì)算量最小的核函數(shù)—重尾徑向基函數(shù)(HTRBF),與普通徑向基函數(shù)相比,它具有指數(shù)級(jí)延遲率,公式為:它包含a、b兩個(gè)參數(shù),默認(rèn) ρ=1,x代表聲音樣本,yi代表第i個(gè)支持向量。對(duì)于HTRBF函數(shù),不同的b值代表不同的分布,但是當(dāng)且僅當(dāng)0≤b≤2的時(shí)候,HTRBF函數(shù)滿足Mercer條件,所以,根據(jù)SVM核函數(shù)的要求,b的取值只可在0到2之間。b的經(jīng)典取值有以下3種:b取0.5時(shí),代表Sublinear分布;b取1.0時(shí),代表Laplacian分布;b取2.0時(shí),代表Gaussian分布。參數(shù)a的選擇不影響Mercer條件,但a取值在0和1之間時(shí),會(huì)決定支持向量數(shù)的衰減速度,從而影響計(jì)算復(fù)雜度。

    表1為a、b取不同值時(shí)對(duì)應(yīng)的識(shí)別正確率,其中,b= 2.0,a=0.25時(shí),識(shí)別率最高,但分類不均勻,完好粒的識(shí)別率相對(duì)較低,因此這對(duì)參數(shù)不是最佳取值;同理b=1.0,a= 0.5時(shí),分類也不均勻;故本文選擇b=2,a=0.125的核函數(shù)作為支持向量的核函數(shù)。

    表1 不同核參數(shù)的SVM識(shí)別正確率(%)

    懲罰參數(shù)是為了解決樣本在高維空中的不完全線性可分問(wèn)題,而合適的參數(shù)C也可以降低誤分率。核函數(shù)參數(shù)固定后,進(jìn)行懲罰參數(shù)C的選擇,表2為C取不同值時(shí)的識(shí)別正確率。表中,C≥4時(shí),識(shí)別結(jié)果均為88%。由于C值越大,算法的時(shí)間復(fù)雜度就越高,因此本文選擇懲罰參數(shù)值C=22=4。

    表2 不同懲罰參數(shù)的SVM識(shí)別正確率(%)

    3.2.2 分類結(jié)果

    選取每類碰撞聲樣本100個(gè),其中50個(gè)用于訓(xùn)練,另外50個(gè)進(jìn)行測(cè)試,識(shí)別結(jié)果如表3所示。3類聲音信號(hào)的正確識(shí)別率均達(dá)到84%以上,且蟲蛀粒的正確率高達(dá)94%。這表明,小麥碰撞聲雙譜特征明顯,支持向量機(jī)分類效果良好,尤其對(duì)蟲蛀粒識(shí)別效果最好。

    表3 3類小麥碰撞聲的識(shí)別結(jié)果

    4 結(jié)束語(yǔ)

    不同麥粒碰撞聲所含信息之間存在差異,本文對(duì)3類小麥碰撞聲進(jìn)行分析,提取了雙譜峰值和對(duì)角切片譜值作為分類依據(jù),用支持向量機(jī)進(jìn)行識(shí)別。得出以下結(jié)論:

    (1)麥粒的內(nèi)部構(gòu)造發(fā)生變化時(shí),會(huì)造成它們碰撞鋼板的聲音有差別,這種差別在雙譜中表現(xiàn)顯著。

    (2)結(jié)合雙譜和支持向量機(jī)方法能夠有效地提取碰撞聲特征,分類效果較好,為小麥粒的識(shí)別研究提供了新方法,同時(shí)也為其他農(nóng)作物的分類研究提供了參考。

    [1]胡麗華,郭敏,張景虎,等.儲(chǔ)糧害蟲檢測(cè)新技術(shù)及應(yīng)用現(xiàn)狀[J].農(nóng)業(yè)工程學(xué)報(bào),2007,23(11):286-290.

    [2]耿森林,尚志遠(yuǎn).儲(chǔ)糧害蟲聲檢測(cè)技術(shù)研究進(jìn)展與展望[J].農(nóng)業(yè)工程學(xué)報(bào),2006,22(4):204-207.

    [3]Pearson T C.Detection of pistachio nuts with closed shells using impactacoustics[J].Applied Engineering in Agriculture,2001,17(2):249-253.

    [4]Pearson T C,Cetin A E,Tewfik A H.Detection of insect damaged wheatkernelsby impactacoustics[C]//ICASSP,2005:649-652.

    [5]Pearson T C,Cetin A E,Tewfik A H,et al.Feasibility of impact-acoustic emissions for detection of damaged wheat kernels[J].Digital signal processing,2007,17:617-633.

    [6]Ince N F,Onaran I,Tewfik A H,et al.Wheat and hazelnut inspection with impact acoustics time-frequency patterns[C]// ASABE Meeting,2007.

    [7]Elbatawi I E.An acoustic impact method to detect hollow heart of potato tubers[J].Biosystems Engineering,2008,100:206-213.

    [8]Khalesi S,Mahmoudi A,Hosainpour A,et al.Detection of walnut varieties using impact acoustics and Artificial Neural Networks(ANNs)[J].Modern Applied Science,2012,6(1):43-49.

    [9]嚴(yán)可國(guó),柳亦兵,徐鴻,等.基于雙譜分析的大型汽輪機(jī)振動(dòng)故障特征提取[J].中國(guó)電機(jī)工程學(xué)報(bào),2007,30(2):98-103.

    ZHANG Yanyan,GUO Min

    College of Computer Science,Shaanxi Normal University,Xi’an 710062,China

    In order to realize the automatic classification of wheat kernels,a new approach that combines the bispectrum and support vector machine is introduced to classify and recognise wheat impact sounds of undamaged kernels,insect damaged kernels and moldy kernels.The impact acoustic signals are processed by bispectrum estimation.Features in bispectrum and diagonal slices spectrum are extracted.Then the features are classified in support vector machine.The recognition accuracy rates in classification of undamaged kernel,insect damaged kernel and moldy kernel are above 84%.The experimental results show that this research has a more comprehensive value in application,and it provides a new method for wheat kernels classification.

    wheat impact acoustic signals;bispectrum estimation;support vector machine

    為實(shí)現(xiàn)小麥顆粒的自動(dòng)分類,采用雙譜和支持向量機(jī)相結(jié)合方法對(duì)小麥完好粒、蟲蛀粒和霉變粒的碰撞聲進(jìn)行分類識(shí)別。對(duì)碰撞聲信號(hào)進(jìn)行雙譜估計(jì),提取信號(hào)雙譜峰值和對(duì)角切片譜兩部分特征,用支持向量機(jī)分類器進(jìn)行分類,對(duì)完好粒、蟲蛀粒和霉變粒3種小麥顆粒識(shí)別正確率均達(dá)84%以上。實(shí)驗(yàn)結(jié)果表明,該研究具有較強(qiáng)的實(shí)際應(yīng)用價(jià)值,為小麥顆粒的分類提供了新的方法和依據(jù)。

    小麥碰撞聲;雙譜估計(jì);支持向量機(jī)

    A

    TP391.42

    10.3778/j.issn.1002-8331.1202-0357

    ZHANG Yanyan,GUO Min.Study on classification of wheat impact acoustic signals based on bispectrum and support vector machine.Computer Engineering and Applications,2013,49(23):203-206.

    國(guó)家自然科學(xué)基金(No.10974130);陜西省教育廳科研計(jì)劃項(xiàng)目(No.11JK0519)。

    張嚴(yán)嚴(yán)(1988—),女,碩士研究生,主要研究方向:信號(hào)處理、模式識(shí)別;郭敏(1964—),通訊作者,女,博士,教授,主要研究領(lǐng)域:數(shù)字信號(hào)處理及模式識(shí)別研究。E-mail:zhangyan-157@163.com

    2012-02-20

    2012-05-08

    1002-8331(2013)23-0203-04

    猜你喜歡
    雙譜蟲蛀對(duì)角
    擬對(duì)角擴(kuò)張Cuntz半群的某些性質(zhì)
    建設(shè)性輿論監(jiān)督:讓公共治理的“大樹”免遭“蟲蛀”
    朱森林養(yǎng)生漫畫
    朱森林養(yǎng)生漫畫
    雙譜圖在語(yǔ)音分析中的應(yīng)用
    蟲蛀川芎揮發(fā)油成分GC-MS研究
    中藥與臨床(2017年3期)2017-10-19 02:42:32
    基于雙譜特征融合的通信輻射源識(shí)別算法
    基于小波包域雙譜的風(fēng)力機(jī)振動(dòng)信號(hào)監(jiān)測(cè)方法研究
    基于倒雙譜的地震子波估計(jì)方法
    非奇異塊α1對(duì)角占優(yōu)矩陣新的實(shí)用簡(jiǎn)捷判據(jù)
    狠狠婷婷综合久久久久久88av| 好男人电影高清在线观看| 王馨瑶露胸无遮挡在线观看| 欧美久久黑人一区二区| 日本vs欧美在线观看视频| 欧美丝袜亚洲另类 | 久久国产精品人妻蜜桃| 色综合欧美亚洲国产小说| 女人久久www免费人成看片| 国产成+人综合+亚洲专区| 视频区图区小说| 久久国产精品人妻蜜桃| 欧美国产精品va在线观看不卡| 女同久久另类99精品国产91| 亚洲精品在线观看二区| 老汉色av国产亚洲站长工具| 麻豆国产av国片精品| 午夜久久久在线观看| netflix在线观看网站| 亚洲熟女毛片儿| 无遮挡黄片免费观看| 精品人妻熟女毛片av久久网站| 免费看十八禁软件| 中文字幕最新亚洲高清| 一二三四在线观看免费中文在| 在线观看一区二区三区激情| 久久国产精品大桥未久av| 别揉我奶头~嗯~啊~动态视频| 国产激情久久老熟女| 咕卡用的链子| tube8黄色片| 精品久久久久久久毛片微露脸| 日韩有码中文字幕| 日日爽夜夜爽网站| 亚洲av成人不卡在线观看播放网| 老司机影院毛片| av福利片在线| 成人精品一区二区免费| 国产在线一区二区三区精| 黄色视频不卡| 99热国产这里只有精品6| 成人18禁高潮啪啪吃奶动态图| 欧美日韩瑟瑟在线播放| 99re在线观看精品视频| 亚洲av欧美aⅴ国产| 99久久综合精品五月天人人| 午夜激情av网站| 亚洲欧美色中文字幕在线| 免费高清在线观看日韩| 狂野欧美激情性xxxx| 欧美乱码精品一区二区三区| 欧美av亚洲av综合av国产av| 满18在线观看网站| 午夜成年电影在线免费观看| 女同久久另类99精品国产91| 老汉色av国产亚洲站长工具| 亚洲人成伊人成综合网2020| 无限看片的www在线观看| 十八禁高潮呻吟视频| videosex国产| 久久这里只有精品19| 日日摸夜夜添夜夜添小说| 亚洲国产欧美日韩在线播放| 亚洲av电影在线进入| 黄色片一级片一级黄色片| 久久 成人 亚洲| 午夜视频精品福利| 天天添夜夜摸| 中出人妻视频一区二区| av免费在线观看网站| 国产淫语在线视频| 国产黄色免费在线视频| 人人妻人人添人人爽欧美一区卜| 天天躁夜夜躁狠狠躁躁| 国产有黄有色有爽视频| 国产男女内射视频| 麻豆乱淫一区二区| 成年人免费黄色播放视频| 后天国语完整版免费观看| 国产在线精品亚洲第一网站| 999久久久国产精品视频| 亚洲精品粉嫩美女一区| 国产亚洲欧美98| 在线观看舔阴道视频| 99国产极品粉嫩在线观看| 90打野战视频偷拍视频| 久久草成人影院| 一进一出抽搐动态| 女人爽到高潮嗷嗷叫在线视频| 亚洲伊人色综图| 免费高清在线观看日韩| 成人三级做爰电影| 中文亚洲av片在线观看爽 | 久久人人97超碰香蕉20202| 美女高潮到喷水免费观看| 99精品久久久久人妻精品| 欧美最黄视频在线播放免费 | 亚洲精品在线观看二区| 老司机影院毛片| 国产亚洲欧美98| 欧美人与性动交α欧美软件| 黄色女人牲交| 精品熟女少妇八av免费久了| 日韩免费高清中文字幕av| 日韩成人在线观看一区二区三区| a级毛片在线看网站| 色尼玛亚洲综合影院| 国产精品一区二区精品视频观看| 人人妻人人爽人人添夜夜欢视频| 人人妻人人澡人人看| 久久精品国产综合久久久| 免费看a级黄色片| 欧美日韩乱码在线| 丰满人妻熟妇乱又伦精品不卡| 亚洲欧美激情在线| 日韩熟女老妇一区二区性免费视频| 国产国语露脸激情在线看| 热re99久久国产66热| 精品人妻1区二区| 亚洲九九香蕉| 男女高潮啪啪啪动态图| 好男人电影高清在线观看| 精品久久久久久久毛片微露脸| 国产区一区二久久| 亚洲专区国产一区二区| 国产一区二区激情短视频| 热99re8久久精品国产| 好看av亚洲va欧美ⅴa在| 国产亚洲精品第一综合不卡| 一边摸一边抽搐一进一小说 | 黑丝袜美女国产一区| 国产区一区二久久| 亚洲欧美激情在线| 成在线人永久免费视频| 国内久久婷婷六月综合欲色啪| 精品国产亚洲在线| 黄色视频,在线免费观看| 狂野欧美激情性xxxx| 老司机在亚洲福利影院| 午夜老司机福利片| 欧美日韩av久久| 亚洲男人天堂网一区| 国产欧美日韩一区二区精品| 高清欧美精品videossex| 久久亚洲精品不卡| 99热国产这里只有精品6| 美国免费a级毛片| 国产一区在线观看成人免费| 制服诱惑二区| 老司机深夜福利视频在线观看| 国产精品香港三级国产av潘金莲| 午夜精品在线福利| 国产精品久久久久成人av| 在线观看免费视频网站a站| a在线观看视频网站| 国精品久久久久久国模美| 91成人精品电影| 69av精品久久久久久| 丰满迷人的少妇在线观看| 女人被躁到高潮嗷嗷叫费观| 亚洲国产精品一区二区三区在线| 亚洲 欧美一区二区三区| 91成人精品电影| 久久亚洲真实| 一区二区三区国产精品乱码| 精品一品国产午夜福利视频| 国产精品秋霞免费鲁丝片| 国产精品久久久久成人av| 国产精品国产av在线观看| 制服诱惑二区| 色老头精品视频在线观看| 91麻豆精品激情在线观看国产 | 免费一级毛片在线播放高清视频 | 91字幕亚洲| 亚洲av第一区精品v没综合| 国内久久婷婷六月综合欲色啪| 欧美亚洲日本最大视频资源| 一夜夜www| 久久久久精品国产欧美久久久| 精品国产超薄肉色丝袜足j| 五月开心婷婷网| 色婷婷av一区二区三区视频| 久久久精品区二区三区| 男女免费视频国产| 91老司机精品| 日本欧美视频一区| 亚洲一卡2卡3卡4卡5卡精品中文| 久久精品亚洲精品国产色婷小说| 中文字幕人妻丝袜一区二区| 成人三级做爰电影| 久久久久久人人人人人| 国产欧美日韩精品亚洲av| 啪啪无遮挡十八禁网站| 国产男女内射视频| 在线观看66精品国产| 国产精品国产av在线观看| 亚洲av欧美aⅴ国产| 首页视频小说图片口味搜索| 91精品国产国语对白视频| 波多野结衣av一区二区av| 十八禁人妻一区二区| 国精品久久久久久国模美| 免费少妇av软件| 嫁个100分男人电影在线观看| 美女福利国产在线| 久热这里只有精品99| 精品熟女少妇八av免费久了| 亚洲五月色婷婷综合| 美女高潮到喷水免费观看| 久久久国产一区二区| 每晚都被弄得嗷嗷叫到高潮| 1024视频免费在线观看| 午夜福利影视在线免费观看| 91精品三级在线观看| 国产成人欧美在线观看 | x7x7x7水蜜桃| 国产亚洲一区二区精品| 最近最新中文字幕大全免费视频| 一区二区三区国产精品乱码| 99精品欧美一区二区三区四区| 国产亚洲欧美98| 一区二区三区激情视频| 51午夜福利影视在线观看| 成年女人毛片免费观看观看9 | 国产又色又爽无遮挡免费看| 丝袜美腿诱惑在线| 免费不卡黄色视频| 免费观看a级毛片全部| 精品久久久久久电影网| 国产激情欧美一区二区| 国产乱人伦免费视频| 免费看a级黄色片| 欧美黄色片欧美黄色片| 亚洲国产精品一区二区三区在线| 极品人妻少妇av视频| 免费观看精品视频网站| 成人特级黄色片久久久久久久| 欧美黑人精品巨大| 在线看a的网站| 国产区一区二久久| 91九色精品人成在线观看| 精品亚洲成a人片在线观看| 丝袜美腿诱惑在线| 亚洲avbb在线观看| 在线观看66精品国产| 99精品欧美一区二区三区四区| 久久香蕉激情| 香蕉久久夜色| 精品国产一区二区三区四区第35| 母亲3免费完整高清在线观看| 操美女的视频在线观看| 亚洲自偷自拍图片 自拍| 99re在线观看精品视频| 午夜福利在线观看吧| 真人做人爱边吃奶动态| 国产又爽黄色视频| 欧美亚洲 丝袜 人妻 在线| 91九色精品人成在线观看| 国产91精品成人一区二区三区| 色婷婷av一区二区三区视频| 九色亚洲精品在线播放| 99在线人妻在线中文字幕 | 亚洲少妇的诱惑av| 人成视频在线观看免费观看| 丰满迷人的少妇在线观看| 好男人电影高清在线观看| 久久中文字幕一级| 高清av免费在线| 久热爱精品视频在线9| 亚洲精品粉嫩美女一区| 国产成人免费观看mmmm| 超碰成人久久| 国产精品久久久人人做人人爽| 欧美乱妇无乱码| 少妇猛男粗大的猛烈进出视频| 国产真人三级小视频在线观看| 国产精品偷伦视频观看了| 国产成人精品久久二区二区免费| 国产精品久久久av美女十八| 欧美日韩黄片免| 午夜福利,免费看| 久久九九热精品免费| 热re99久久国产66热| 亚洲精品在线观看二区| 99国产精品一区二区蜜桃av | 国产一区二区三区视频了| 成人免费观看视频高清| 国产成人免费观看mmmm| 亚洲色图av天堂| 国产91精品成人一区二区三区| 热99re8久久精品国产| 在线观看www视频免费| 69av精品久久久久久| 亚洲一区二区三区不卡视频| 亚洲熟妇熟女久久| 国产高清国产精品国产三级| 成年动漫av网址| 亚洲一区高清亚洲精品| 91成人精品电影| 久久久精品国产亚洲av高清涩受| 午夜精品国产一区二区电影| 两性夫妻黄色片| 少妇 在线观看| 99国产精品一区二区蜜桃av | 丝袜人妻中文字幕| 国产精品98久久久久久宅男小说| 亚洲av第一区精品v没综合| 9191精品国产免费久久| 成人免费观看视频高清| 中文字幕高清在线视频| 国产三级黄色录像| 男男h啪啪无遮挡| 男女床上黄色一级片免费看| 国产精品一区二区免费欧美| 欧美老熟妇乱子伦牲交| 久久ye,这里只有精品| 久久狼人影院| av网站在线播放免费| 精品亚洲成国产av| 国产成人免费观看mmmm| 久久久国产欧美日韩av| 操美女的视频在线观看| 精品亚洲成国产av| 可以免费在线观看a视频的电影网站| 亚洲伊人色综图| 一进一出抽搐动态| 黑人巨大精品欧美一区二区蜜桃| 99热国产这里只有精品6| 身体一侧抽搐| 在线观看日韩欧美| 中文字幕制服av| 免费观看a级毛片全部| 国产99白浆流出| 黄片大片在线免费观看| 丰满饥渴人妻一区二区三| 日韩人妻精品一区2区三区| 又黄又爽又免费观看的视频| 欧美激情高清一区二区三区| 超碰成人久久| 亚洲av熟女| 在线十欧美十亚洲十日本专区| 欧美激情高清一区二区三区| 热99国产精品久久久久久7| 中文字幕人妻丝袜制服| 亚洲午夜理论影院| 久久人妻福利社区极品人妻图片| av不卡在线播放| 狠狠婷婷综合久久久久久88av| 色综合婷婷激情| 岛国毛片在线播放| 啦啦啦视频在线资源免费观看| 亚洲熟妇熟女久久| 亚洲中文日韩欧美视频| 俄罗斯特黄特色一大片| 久久精品亚洲熟妇少妇任你| 一夜夜www| av不卡在线播放| 黄频高清免费视频| 欧美精品啪啪一区二区三区| 久久99一区二区三区| 色婷婷久久久亚洲欧美| 大型av网站在线播放| 免费在线观看视频国产中文字幕亚洲| 男女之事视频高清在线观看| 欧美精品啪啪一区二区三区| av超薄肉色丝袜交足视频| 美女国产高潮福利片在线看| 母亲3免费完整高清在线观看| 成年版毛片免费区| 丝袜人妻中文字幕| 欧美亚洲日本最大视频资源| 女警被强在线播放| 亚洲伊人色综图| 操美女的视频在线观看| 在线观看一区二区三区激情| 国产精品久久电影中文字幕 | 狂野欧美激情性xxxx| 亚洲欧美色中文字幕在线| 久久草成人影院| 国产xxxxx性猛交| 精品欧美一区二区三区在线| 人妻 亚洲 视频| 国产精品欧美亚洲77777| 亚洲精品久久成人aⅴ小说| 欧美中文综合在线视频| 国产精品国产av在线观看| 中文字幕最新亚洲高清| 久久久久精品人妻al黑| 精品久久久精品久久久| 啦啦啦 在线观看视频| 免费一级毛片在线播放高清视频 | 精品人妻熟女毛片av久久网站| xxxhd国产人妻xxx| 午夜福利在线免费观看网站| 超碰成人久久| 国产精品久久视频播放| 黄网站色视频无遮挡免费观看| 两性夫妻黄色片| 悠悠久久av| 99久久国产精品久久久| 色在线成人网| 欧洲精品卡2卡3卡4卡5卡区| 午夜免费观看网址| 高潮久久久久久久久久久不卡| 亚洲熟妇熟女久久| 色播在线永久视频| 亚洲av日韩在线播放| 757午夜福利合集在线观看| 亚洲av日韩在线播放| 欧美日韩精品网址| 黄片大片在线免费观看| 亚洲专区中文字幕在线| 女人精品久久久久毛片| 精品一品国产午夜福利视频| 久久人人爽av亚洲精品天堂| 岛国毛片在线播放| 无限看片的www在线观看| 亚洲一码二码三码区别大吗| 女同久久另类99精品国产91| 99精品久久久久人妻精品| 国产深夜福利视频在线观看| 91麻豆av在线| 丰满的人妻完整版| 人人妻人人添人人爽欧美一区卜| 一个人免费在线观看的高清视频| 国产免费男女视频| 视频在线观看一区二区三区| 成人av一区二区三区在线看| e午夜精品久久久久久久| 国产一区在线观看成人免费| 制服人妻中文乱码| 女性生殖器流出的白浆| 国产成+人综合+亚洲专区| 亚洲黑人精品在线| 日本黄色视频三级网站网址 | 精品一区二区三卡| 最近最新中文字幕大全电影3 | 欧美日韩成人在线一区二区| 亚洲专区字幕在线| 一级a爱片免费观看的视频| 桃红色精品国产亚洲av| 曰老女人黄片| 黄频高清免费视频| 香蕉久久夜色| 亚洲精华国产精华精| 亚洲熟女毛片儿| 日韩有码中文字幕| 天天影视国产精品| 日韩欧美在线二视频 | 亚洲国产欧美网| 国产亚洲精品久久久久5区| 亚洲欧洲精品一区二区精品久久久| 国产精品久久久人人做人人爽| 50天的宝宝边吃奶边哭怎么回事| 成在线人永久免费视频| 曰老女人黄片| 波多野结衣一区麻豆| 精品人妻1区二区| 欧美日韩亚洲国产一区二区在线观看 | 亚洲av成人不卡在线观看播放网| 成人国语在线视频| 国产精品香港三级国产av潘金莲| 人人妻人人添人人爽欧美一区卜| 亚洲九九香蕉| 麻豆成人av在线观看| 一级作爱视频免费观看| 欧美av亚洲av综合av国产av| 制服人妻中文乱码| 精品久久久久久久毛片微露脸| www日本在线高清视频| 亚洲精品自拍成人| 91字幕亚洲| 91麻豆av在线| 高清视频免费观看一区二区| 91九色精品人成在线观看| 久久久久久久久久久久大奶| 午夜视频精品福利| 欧美一级毛片孕妇| 亚洲成av片中文字幕在线观看| 亚洲av片天天在线观看| 国产精品 国内视频| 黄色 视频免费看| 99热网站在线观看| 久久久精品免费免费高清| 我的亚洲天堂| 午夜亚洲福利在线播放| 一区二区三区激情视频| 亚洲成a人片在线一区二区| 国产不卡一卡二| xxx96com| 日韩欧美免费精品| 女同久久另类99精品国产91| 精品国产一区二区久久| 一区二区日韩欧美中文字幕| 日日夜夜操网爽| 18禁美女被吸乳视频| 天堂俺去俺来也www色官网| 性少妇av在线| 午夜免费观看网址| 亚洲欧美一区二区三区久久| a级毛片在线看网站| 国产欧美日韩一区二区三区在线| 久久久国产一区二区| 亚洲av成人av| 亚洲专区字幕在线| 国产精品成人在线| tube8黄色片| xxx96com| 狠狠婷婷综合久久久久久88av| 久久香蕉精品热| 老司机午夜十八禁免费视频| 日韩欧美一区视频在线观看| 窝窝影院91人妻| 亚洲精品中文字幕一二三四区| 在线观看www视频免费| 亚洲免费av在线视频| 另类亚洲欧美激情| 男男h啪啪无遮挡| 老鸭窝网址在线观看| 18禁观看日本| √禁漫天堂资源中文www| 亚洲熟妇熟女久久| 亚洲午夜精品一区,二区,三区| 欧美人与性动交α欧美精品济南到| 精品亚洲成国产av| 精品久久久久久久久久免费视频 | a在线观看视频网站| 十八禁高潮呻吟视频| 中文字幕最新亚洲高清| 国产av精品麻豆| 黄色成人免费大全| 日韩一卡2卡3卡4卡2021年| 一进一出好大好爽视频| 欧美乱码精品一区二区三区| 桃红色精品国产亚洲av| 波多野结衣av一区二区av| 亚洲情色 制服丝袜| 少妇裸体淫交视频免费看高清 | 久久国产乱子伦精品免费另类| 亚洲精品国产色婷婷电影| 免费高清在线观看日韩| 中文亚洲av片在线观看爽 | 1024香蕉在线观看| 亚洲一区二区三区不卡视频| 成年人黄色毛片网站| 欧美精品高潮呻吟av久久| 精品国产乱子伦一区二区三区| 91av网站免费观看| 国产成人欧美| 国产激情久久老熟女| 亚洲av日韩精品久久久久久密| 国产精品1区2区在线观看. | 老司机靠b影院| 欧美成人免费av一区二区三区 | 中文字幕最新亚洲高清| 日韩欧美免费精品| 亚洲国产欧美日韩在线播放| 午夜久久久在线观看| 国产成人精品在线电影| 亚洲中文字幕日韩| 大码成人一级视频| 国产三级黄色录像| 一级黄色大片毛片| 国产精品亚洲av一区麻豆| 在线免费观看的www视频| 成人永久免费在线观看视频| 麻豆国产av国片精品| 成年动漫av网址| 亚洲国产欧美网| 大香蕉久久成人网| 丁香欧美五月| 国产成人免费无遮挡视频| 麻豆av在线久日| 亚洲美女黄片视频| 国产99久久九九免费精品| 69av精品久久久久久| www.自偷自拍.com| 久久热在线av| 国产成人精品久久二区二区91| 91字幕亚洲| 久久影院123| 狠狠婷婷综合久久久久久88av| 飞空精品影院首页| www.熟女人妻精品国产| 亚洲性夜色夜夜综合| 国产亚洲一区二区精品| 成熟少妇高潮喷水视频| 国产精品一区二区在线观看99| 日本黄色视频三级网站网址 | 亚洲精品美女久久久久99蜜臀| 757午夜福利合集在线观看| 黄色片一级片一级黄色片| a级片在线免费高清观看视频| 高清欧美精品videossex| 亚洲美女黄片视频| xxxhd国产人妻xxx| 一级,二级,三级黄色视频| 女人高潮潮喷娇喘18禁视频| 亚洲在线自拍视频| 电影成人av| 亚洲 国产 在线| 9色porny在线观看| 交换朋友夫妻互换小说| 老熟女久久久| 午夜免费鲁丝| 国产精品久久视频播放| 欧美精品高潮呻吟av久久| 美女福利国产在线| 中文欧美无线码| 巨乳人妻的诱惑在线观看| 亚洲成人手机| 欧美日韩视频精品一区| 亚洲精品久久成人aⅴ小说| 国产亚洲精品一区二区www | 色在线成人网|