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

        ?

        基于子帶譜熵法和PSO-GA-SVM的汽車?guó)Q笛識(shí)別

        2021-10-20 09:28:44余凌浩陸鐵文曾毓敏
        關(guān)鍵詞:鳴笛子帶譜線

        余凌浩,陸鐵文,李 晨,曾毓敏,袁 芳

        (1.南京師范大學(xué)計(jì)算機(jī)與電子信息學(xué)院/人工智能學(xué)院,江蘇 南京 210023)(2.杭州愛(ài)華智能科技有限公司,浙江 杭州 311121)

        近年來(lái),交通噪聲的關(guān)注度日益提升. 鳴笛噪聲是交通噪聲的主要組成成分,汽車在禁鳴區(qū)違章鳴笛,會(huì)嚴(yán)重影響他人的生活、學(xué)習(xí)與工作. 為了遏制違章鳴笛,鳴笛抓拍系統(tǒng)應(yīng)運(yùn)而生. 鳴笛抓拍系統(tǒng)會(huì)實(shí)時(shí)檢測(cè)當(dāng)前的聲壓級(jí),當(dāng)聲壓級(jí)超過(guò)一定閾值時(shí),系統(tǒng)通過(guò)遠(yuǎn)場(chǎng)波束形成技術(shù)定位聲源的位置并進(jìn)行音視頻記錄[1]. 然而交通環(huán)境中聲源情況非常復(fù)雜,有很多非鳴笛聲也能觸發(fā)鳴笛抓拍系統(tǒng)的閾值,只依靠聲壓級(jí)觸發(fā)門限會(huì)出現(xiàn)較多誤判的情況. 這也是目前鳴笛抓拍系統(tǒng)應(yīng)用較少的主要原因. 因此,對(duì)鳴笛抓拍系統(tǒng)采集的樣本進(jìn)行鳴笛與非鳴笛的分類是非常有必要的.

        在汽車?guó)Q笛識(shí)別方面,目前已有許多學(xué)者進(jìn)行了研究[2-8],較為常用的方法主要分為以下兩類:(1)閾值判斷,如文獻(xiàn)[8]中使用的短時(shí)能量和短時(shí)過(guò)零率法等. 該方法將特征值高于或低于一定門限的聲音樣本直接進(jìn)行分類,其分類效果主要取決于特征和閾值的選取. (2)機(jī)器學(xué)習(xí),如文獻(xiàn)[2]中使用的卷積神經(jīng)網(wǎng)絡(luò)、文獻(xiàn)[6]中使用的隱馬爾科夫模型等. 該方法首先對(duì)樣本提取特征,通過(guò)訓(xùn)練一定數(shù)量的樣本得到一個(gè)模型,再通過(guò)采集的樣本是否與模型相匹配來(lái)進(jìn)行分類.

        進(jìn)行閾值判斷時(shí),閾值的設(shè)置會(huì)直接影響分類結(jié)果的準(zhǔn)確性,然而交通環(huán)境中聲源情況非常復(fù)雜,會(huì)給閾值設(shè)置帶來(lái)較大阻礙. 機(jī)器學(xué)習(xí)雖然準(zhǔn)確率較高,但對(duì)模型有比較高的要求. 如果實(shí)測(cè)的鳴笛樣本中含有較多噪聲,也會(huì)對(duì)分類準(zhǔn)確率產(chǎn)生很大影響. 因此,我們需要結(jié)合實(shí)際交通環(huán)境中的鳴笛噪聲,研究具有較高準(zhǔn)確率的鳴笛識(shí)別算法.

        由于鳴笛聲具有豐富的諧波特征,其子帶譜熵值相對(duì)較低. 相比較于傳統(tǒng)的短時(shí)能量和短時(shí)過(guò)零率特征,子帶譜熵具有更高的區(qū)分度[9-10],閾值的設(shè)置具有普適性. 機(jī)器學(xué)習(xí)主要分為傳統(tǒng)機(jī)器學(xué)習(xí)和深度學(xué)習(xí)兩大類. 深度學(xué)習(xí)非常適合處理大數(shù)據(jù)下的特征,但是在較小的數(shù)據(jù)量下,深度學(xué)習(xí)并沒(méi)有表現(xiàn)出很好的效果,并且由于其計(jì)算量過(guò)大,很難滿足實(shí)時(shí)性的需求. 而支持向量機(jī)(support vector machine,SVM)作為傳統(tǒng)機(jī)器學(xué)習(xí)的方法,在進(jìn)行小樣本集分類時(shí)效果較好[11].

        SVM參數(shù)優(yōu)化的傳統(tǒng)方法是網(wǎng)格搜索法,該方法通過(guò)遍歷參數(shù)的所有組合來(lái)選取最優(yōu)結(jié)果,運(yùn)行速度較慢. 我們使用粒子群優(yōu)化(particle swarm optimization,PSO)和遺傳算法(genetic algorithm,GA)的融合(PSO-GA)來(lái)優(yōu)化SVM的參數(shù)[12-13],在準(zhǔn)確率沒(méi)有降低的情況下,提高了選取最優(yōu)參數(shù)的速度.

        本文結(jié)合了閾值判斷和機(jī)器學(xué)習(xí)的方法,提出了一種基于子帶譜熵法和PSO-GA-SVM的汽車?guó)Q笛識(shí)別算法. 仿真結(jié)果表明,該算法具有較好的魯棒性. 在對(duì)實(shí)際采集樣本的鳴笛識(shí)別中,本文算法也取得了較高的準(zhǔn)確率.

        1 算法介紹

        1.1 算法總流程

        本文首先使用子帶譜熵法對(duì)鳴笛抓拍系統(tǒng)采集的聲音樣本進(jìn)行初判,將子帶譜熵高于閾值的樣本直接判定為非鳴笛樣本. 為了使SVM訓(xùn)練出較好的模型,在子帶譜熵法中再結(jié)合子帶譜線能量門限,對(duì)初判為鳴笛的樣本中的疑似鳴笛部分進(jìn)行分割. 然后使用SVM對(duì)分割結(jié)果進(jìn)行進(jìn)一步分類,并使用PSO-GA來(lái)優(yōu)化SVM的參數(shù). 算法總流程如圖1所示.

        圖1 算法總流程Fig.1 Overall algorithm process

        1.2 預(yù)處理

        文中預(yù)處理的主要流程為:帶通濾波、預(yù)加重、分幀、加窗.

        一般汽車?guó)Q笛的頻率分布在1 500~6 000 Hz. 為了保留汽車?guó)Q笛的主要特征并去除噪聲的干擾,對(duì)聲音樣本進(jìn)行1 500~6 000 Hz的帶通濾波.

        預(yù)加重是一個(gè)一階高通濾波器,用一階FIR濾波器表示為:

        s′(n)=s(n)-as(n-1).

        (1)

        式中,a為一常數(shù),s(n)為原始信號(hào),s′(n)為預(yù)加重后的信號(hào). 預(yù)加重可以對(duì)信號(hào)進(jìn)行高頻提升,使頻譜變得更平坦,減少頻譜的動(dòng)態(tài)范圍. 一般來(lái)說(shuō),0.9

        分幀可以將一段較長(zhǎng)聲音信號(hào)分成多個(gè)小段,使每小段聲音信號(hào)可以作為平穩(wěn)信號(hào)來(lái)處理. 為了使每一幀信號(hào)之間平滑地過(guò)渡,兩個(gè)相鄰幀之間需要有一些重疊區(qū)域. 本文中幀長(zhǎng)選為50 ms,幀移長(zhǎng)度選為幀長(zhǎng)的1/3.

        加窗可以避免因信號(hào)兩端不連續(xù)而導(dǎo)致的頻譜泄漏. 本文選用漢明窗作為窗函數(shù)w(n),若共有N個(gè)采樣點(diǎn),則第n個(gè)采樣點(diǎn)的窗函數(shù)公式為:

        (2)

        1.3 子帶譜熵法

        1.3.1 算法介紹

        子帶譜熵的思想是將每一幀的頻譜分成若干個(gè)子帶,這樣就消除了每一條譜線受噪聲影響的問(wèn)題. 設(shè)每個(gè)子帶由4條譜線組成,共有Nb個(gè)子帶. 第i幀中的第n個(gè)子帶的譜線能量E(n,i)為:

        (3)

        式中,Yi(k)為傅里葉變換后的譜線能量. 子帶譜線能量的概率分布p(n,i)為:

        (4)

        引入一個(gè)正常量C,得到新的子帶譜線能量的概率分布式:

        (5)

        文獻(xiàn)[9]中推導(dǎo)證明,在噪聲環(huán)境下,引入正常量C后,分割的準(zhǔn)確度能得到提高. 本文中C取0.1. 子帶譜熵Hb(i)為:

        (6)

        1.3.2 算法流程

        子帶譜熵法的主要流程如圖2所示. 在信號(hào)的子帶譜線能量包絡(luò)線上選取一個(gè)較高的閾值eth2,包絡(luò)線與eth2相交于左右兩個(gè)點(diǎn)A和B. 認(rèn)定高于eth2的部分一定是鳴笛,而鳴笛的起止點(diǎn)應(yīng)該在AB段之外. 再在信號(hào)的子帶譜線能量包絡(luò)線上選取一個(gè)較低的閾值eth1,從A點(diǎn)向左、從B點(diǎn)向右搜索,分別找到包絡(luò)線與eth1相交的左右兩個(gè)點(diǎn)C和D,如圖3(a)所示.CD段即為基于子帶譜線能量門限分割的鳴笛部分.

        圖2 子帶譜熵法流程Fig.2 Sub-band spectral entropy method process

        圖3 雙門限分割示意圖Fig.3 Schematic diagram of double threshold segmentation

        在對(duì)應(yīng)的子帶譜熵包絡(luò)線上,從C點(diǎn)向右、從D點(diǎn)向左進(jìn)行更加細(xì)致的分割,找到子帶譜熵低于閾值eth3的左右兩個(gè)點(diǎn)E和F,如圖3(b)所示.EF段即為基于子帶譜線能量和子帶譜熵雙門限分割的鳴笛部分. 圖3(c)為聲音樣本的頻譜圖以及分割結(jié)果,左邊的豎線為分割的起始點(diǎn),右邊的豎線為分割的停止點(diǎn).

        由于汽車?guó)Q笛聲具有豐富的諧波特征,其子帶譜熵值相對(duì)較低,子帶譜熵特征的區(qū)分度更高. 為了驗(yàn)證子帶譜熵法的有效性,我們使用文獻(xiàn)[8]的短時(shí)能量和短時(shí)過(guò)零率法與本文的子帶譜熵法進(jìn)行了對(duì)比實(shí)驗(yàn). 實(shí)驗(yàn)結(jié)果如圖4和圖5所示.

        圖4的樣本中沒(méi)有鳴笛聲,正確的結(jié)果應(yīng)沒(méi)有分割線,如圖4(a)所示. 文獻(xiàn)[8]的方法產(chǎn)生了分割線,而本文的方法沒(méi)有產(chǎn)生分割線,并將其分類為非鳴笛樣本. 圖5的樣本中含有鳴笛和噪聲,正確的分割結(jié)果如圖5(a)所示. 文獻(xiàn)[8]的方法不僅分割出了鳴笛部分,還分割出了多個(gè)噪聲片段. 而本文的方法則準(zhǔn)確地分割出了鳴笛聲,沒(méi)有分割出噪聲部分.

        圖4 對(duì)一段非鳴笛聲音樣本的分割Fig.4 Segmentation of a non-whistle sound samples

        圖5 對(duì)一段含有鳴笛和噪聲的聲音樣本的分割Fig.5 Segmentation of a sound samples with car whistle and noise

        使用子帶譜熵法可以將子帶譜熵高于閾值的樣本直接判為非鳴笛樣本,如圖4(c)所示. 然而在實(shí)際交通環(huán)境中也存在一些非鳴笛聲(如汽車轉(zhuǎn)向提示聲、公交車報(bào)站聲、警笛聲等),由于具有一定的諧波特征,產(chǎn)生了分割結(jié)果,會(huì)被子帶譜熵法誤判為鳴笛聲,如圖6所示. 因此,我們需要對(duì)子帶譜熵法的分割結(jié)果進(jìn)行進(jìn)一步的分類,來(lái)消除此類誤判.

        圖6 對(duì)部分非鳴笛樣本的分割Fig.6 Segmentation of some non-whistle sound samples

        1.4 PSO-GA-SVM

        1.4.1 算法介紹

        SVM是一種二分類模型. 設(shè)樣本集為(xi,yi),其中xi為樣本特征,yi為類別符號(hào). SVM的目標(biāo)是尋找一個(gè)分類面ωTx+b=0將樣本正確地分成兩類,且使‖ω2‖最小. 在樣本集線性不可分時(shí),需要引入懲罰因子c和核函數(shù)K(xi,xj)來(lái)進(jìn)行分析. 核函數(shù)可以將低維向量?jī)?nèi)積的結(jié)果轉(zhuǎn)化為高維向量?jī)?nèi)積的結(jié)果,但不需要在高維中進(jìn)行計(jì)算. 本文選用的核函數(shù)為徑向基核函數(shù),

        K(xi,xj)=e-g‖xi-xj‖2.

        (7)

        懲罰因子c和核函數(shù)參數(shù)g對(duì)SVM的分類性能有很大的影響. 本文使用PSO-GA來(lái)優(yōu)化c和g參數(shù).

        PSO將粒子i在搜索空間中的飛行速度表示為向量vi,位置表示為向量xi. 在每一次迭代中,計(jì)算每個(gè)粒子的適應(yīng)度,通過(guò)個(gè)體極值pbest和全局極值gbest來(lái)更新每一個(gè)粒子的速度和位置,如下兩式所示:

        vi=K[vi+φ1r(pbest-xi)+φ2r(gbest-xi)],

        (8)

        xi=xi+vi.

        (9)

        式(8)中,r為均勻分布在(0,1)之間的隨機(jī)數(shù),K由下式表示[14]:

        (10)

        本文中φ取4.1,則K取0.730.

        GA將搜索空間中每一個(gè)可能的解編碼為染色體,隨機(jī)選擇一組染色體來(lái)形成初始群體. 在每一次迭代中,計(jì)算每個(gè)粒子的適應(yīng)度,按照適應(yīng)度大小的順序來(lái)選擇個(gè)體,進(jìn)行交叉、變異等操作. 反復(fù)執(zhí)行遺傳操作,直到滿足終止條件為止.

        1.4.2 算法流程

        PSO-GA-SVM主要流程如圖7所示,其主要步驟如下:

        圖7 PSO-GA-SVM流程Fig.7 PSO-GA-SVM process

        步驟1 對(duì)訓(xùn)練集進(jìn)行k折交叉驗(yàn)證. 本文中k取5.

        步驟2 將搜索空間中每一個(gè)c和g可能的解編碼為染色體,隨機(jī)選擇N個(gè)染色體來(lái)形成初始種群,并初始化粒子速度. 計(jì)算每一個(gè)粒子的適應(yīng)度,求出pbest和gbest. 本文中N取20.

        步驟3 按照式(8)和(9)更新每一個(gè)粒子的速度和位置.

        步驟4 按照適應(yīng)度選出N-1個(gè)粒子,對(duì)它們進(jìn)行交叉、變異操作,再與上一代的最優(yōu)粒子組合成新的N個(gè)粒子. 計(jì)算每一個(gè)粒子的適應(yīng)度,更新pbest和gbest.

        步驟5 轉(zhuǎn)步驟3,直到滿足一定迭代次數(shù)為止. 本文中迭代次數(shù)取100.

        步驟6 輸出最優(yōu)c和g.

        步驟7 將最優(yōu)c和g代入SVM的c和g參數(shù)中,對(duì)測(cè)試集進(jìn)行分類.

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

        本文使用CPU為Intel i7-8750H、內(nèi)存為8G的計(jì)算機(jī)作為硬件平臺(tái). 本文中所有音頻的格式為wav格式,位深度為32位. 所有程序均在Python 3.7.0下編寫(xiě). 每次實(shí)驗(yàn)從數(shù)據(jù)集中隨機(jī)選取80%的樣本作為訓(xùn)練集,剩下20%的樣本作為測(cè)試集,計(jì)算100次實(shí)驗(yàn)的平均準(zhǔn)確率. 本文選用Mel頻率倒譜系數(shù)作為分類器輸入的聲音特征[15],其中Mel濾波器的個(gè)數(shù)選為96個(gè),離散余弦變換后取第2~65維,再對(duì)每一幀信號(hào)求平均,最終形成一個(gè)64維的向量.

        2.1 仿真實(shí)驗(yàn)

        本節(jié)中對(duì)不同底噪、不同信噪比的聲音樣本進(jìn)行仿真實(shí)驗(yàn). 鳴笛聲取自杭州愛(ài)華智能科技有限公司提供的純凈鳴笛聲樣本,包含奧迪A3、大眾高爾夫、吉利等13種不同車型的鳴笛聲,共120個(gè). 非鳴笛聲取自UrbanSound8K數(shù)據(jù)集中的城市聲音樣本,包含空調(diào)外機(jī)聲、小孩玩耍聲、警笛聲等8種類型的聲音,共211個(gè). 本節(jié)中音頻的時(shí)長(zhǎng)為1 s,采樣頻率為48 kHz. 為了模擬鳴笛抓拍系統(tǒng)的采集方式,每個(gè)聲音樣本中只有0.25 s存在待識(shí)別的聲音,并且隨機(jī)設(shè)置該聲音在音頻中的起始位置. 為了驗(yàn)證本文算法的魯棒性,在上述聲音樣本的基礎(chǔ)上,分別添加信噪比為20 dB、10 dB、5 dB、0 dB的白噪聲和粉噪聲作為底噪進(jìn)行仿真實(shí)驗(yàn). 白噪聲和粉噪聲取自NOISE92數(shù)據(jù)集. 仿真實(shí)驗(yàn)結(jié)果如表1所示.

        分析表1可知,在信噪比較高的情況下,本文算法對(duì)仿真樣本的識(shí)別準(zhǔn)確率可以達(dá)到99%以上. 隨著信噪比不斷降低,本文算法的識(shí)別效果雖然也有所下降,但依然維持在相對(duì)較高的準(zhǔn)確率上. 在極差的信噪比環(huán)境下(信噪比為0 dB時(shí)),本文算法依然能夠有75%以上的準(zhǔn)確率. 說(shuō)明本文算法具有較好的魯棒性.

        2.2 結(jié)果分析

        本節(jié)中對(duì)杭州愛(ài)華智能科技有限公司的鳴笛抓拍系統(tǒng)實(shí)際采集的1 067個(gè)聲音樣本進(jìn)行鳴笛識(shí)別,其中包括683個(gè)不含有鳴笛聲的樣本,包含卡車發(fā)動(dòng)機(jī)聲、汽車轉(zhuǎn)向提示聲、公交車報(bào)站聲、警笛聲等超過(guò)10種城市環(huán)境音. 包括384個(gè)含有鳴笛聲的樣本,包含各種車型的鳴笛聲,背景噪聲中也存在多種上述城市環(huán)境音. 本節(jié)中音頻的時(shí)長(zhǎng)為3.8 s,采樣頻率為24 kHz. 我們選擇文獻(xiàn)[2]算法、文獻(xiàn)[7]算法與本文算法進(jìn)行了對(duì)比實(shí)驗(yàn). 實(shí)驗(yàn)結(jié)果如表2所示.

        從表2可見(jiàn)出,本文算法相比文獻(xiàn)[2]算法準(zhǔn)確率提高了21.3%,相比文獻(xiàn)[7]算法準(zhǔn)確率提高了14.4%. 文獻(xiàn)[2]和文獻(xiàn)[7]只使用深度學(xué)習(xí)的方法對(duì)聲音樣本進(jìn)行分類,而本文首先使用子帶譜熵法對(duì)聲音樣本進(jìn)行初判和分割,再使用分類器對(duì)分割結(jié)果進(jìn)行進(jìn)一步分類. 同時(shí),文獻(xiàn)[2]和文獻(xiàn)[7]中深度學(xué)習(xí)的方法在樣本數(shù)量較少時(shí)分類效果并不理想,而本文使用的PSO-GA-SVM更加適合小樣本集的分類. 因此,本文算法在汽車?guó)Q笛識(shí)別上具有更高的準(zhǔn)確率.

        表2 不同算法的平均準(zhǔn)確率Table 2 Average accuracy with different method

        我們使用子帶譜熵法結(jié)合文獻(xiàn)[7]使用的反向傳播神經(jīng)網(wǎng)絡(luò)(backpropagation neural network,BPNN)和本文使用的PSO-GA-SVM進(jìn)一步進(jìn)行了對(duì)比試驗(yàn). 實(shí)驗(yàn)結(jié)果如表3所示.

        表3 使用子帶譜熵法前后的平均準(zhǔn)確率Table 3 Average accuracy with and without sub-bandspectral entropy method

        子帶譜熵法對(duì)非鳴笛樣本的誤判率在1%以內(nèi),但是對(duì)鳴笛樣本的誤判率在34%以上,不能滿足鳴笛抓拍系統(tǒng)對(duì)準(zhǔn)確率的要求,因此需要對(duì)子帶譜熵法初判為鳴笛的樣本作進(jìn)一步的分類. 由表3可知,使用子帶譜熵法前,BPNN和PSO-GA-SVM的準(zhǔn)確率并不是非常理想,而使用子帶譜熵法后,兩種分類器的準(zhǔn)確率均有了較大提高. 這是因?yàn)樽訋ёV熵法可以分割出疑似鳴笛的部分,大幅減少了噪聲對(duì)分類器模型造成的影響. 同時(shí)從表3中還可以看出,本文使用的PSO-GA-SVM相比較于BPNN具有更高的準(zhǔn)確率,更加適合小樣本集的分類.

        此外,由于本文使用的聲音特征維數(shù)較低,僅有64維,提取特征以及分類的速度較快,在本文使用的計(jì)算機(jī)上處理一個(gè)樣本僅需約0.3 s,可以很好地應(yīng)用到實(shí)時(shí)檢測(cè)中.

        3 結(jié)論

        針對(duì)現(xiàn)有的汽車?guó)Q笛識(shí)別算法易受噪聲影響、計(jì)算量較大的問(wèn)題,本文提出如下算法:首先利用閾值判斷對(duì)聲音樣本進(jìn)行初判,使用子帶譜熵法分割出樣本中疑似的鳴笛部分;然后利用機(jī)器學(xué)習(xí)對(duì)疑似鳴笛樣本作進(jìn)一步分類,使用PSO-GA-SVM訓(xùn)練模型,最終分類出所有的鳴笛與非鳴笛樣本. 實(shí)驗(yàn)表明,本文算法準(zhǔn)確率較高且計(jì)算量較小,可以滿足鳴笛抓拍系統(tǒng)實(shí)時(shí)檢測(cè)的要求.

        猜你喜歡
        鳴笛子帶譜線
        一種基于奇偶判斷WPT的多音干擾抑制方法*
        基于HITRAN光譜數(shù)據(jù)庫(kù)的合并譜線測(cè)溫仿真研究
        朋友圈
        子帶編碼在圖像壓縮編碼中的應(yīng)用
        電子制作(2019年22期)2020-01-14 03:16:24
        鐵合金光譜譜線分離實(shí)驗(yàn)研究
        鍶原子光鐘鐘躍遷譜線探測(cè)中的程序控制
        基于虛擬孔徑擴(kuò)展的子帶信息融合寬帶DOA估計(jì)
        藥芯焊絲GMAW電弧光譜的研究
        加拿大司機(jī)很少鳴笛
        基于子帶模式的AMC技術(shù)算法分析
        中文有码人妻字幕在线 | 99国产精品丝袜久久久久| 激情五月婷婷六月俺也去| 日本不卡一区二区三区在线视频 | 亚洲av成人无码网站…| Y111111国产精品久久久| 男女啪啪免费视频网址| 男女男精品视频网站免费看| 无码人妻久久一区二区三区免费| 欧美国产亚洲日韩在线二区| 日本一区二区三区四区在线看| 白白色发布会在线观看免费| 亚洲热线99精品视频| 免费啪啪视频一区| 日韩av一区在线播放| 国产成人一区二区三区乱| 18禁真人抽搐一进一出在线| 国产a级午夜毛片| 久久综合老鸭窝色综合久久| 免费a级毛片在线播放| 国产精品三级在线观看无码| 亚洲福利天堂网福利在线观看| 青青草是针对华人绿色超碰| 亚洲日韩精品无码av海量| 日本乱子人伦在线视频| 无码一区二区三区人| 国内免费自拍9偷1拍| 一本大道无码人妻精品专区| 久久网视频中文字幕综合| 在线观看黄片在线播放视频| 白嫩人妻少妇偷人精品| 国产97在线 | 中文| 中文字幕人成人乱码亚洲| 永久免费观看的黄网站在线| 国产成人无码精品久久久露脸 | 国产精品99久久久久久宅男| 久久精品国产亚洲片| 日本中文一区二区在线| 一个人看的视频www免费| 国产老妇伦国产熟女老妇高清| 精品一区二区三区a桃蜜|