999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

電動(dòng)汽車(chē)聲品質(zhì)的評(píng)價(jià)分析及建模

2016-06-12 05:22:51胡騰陸益民合肥工業(yè)大學(xué)合肥230009
汽車(chē)技術(shù) 2016年3期
關(guān)鍵詞:評(píng)價(jià)模型

胡騰陸益民(合肥工業(yè)大學(xué),合肥230009)

?

電動(dòng)汽車(chē)聲品質(zhì)的評(píng)價(jià)分析及建模

胡騰陸益民
(合肥工業(yè)大學(xué),合肥230009)

【摘要】針對(duì)A計(jì)權(quán)聲壓級(jí)不能準(zhǔn)確衡量人對(duì)電動(dòng)汽車(chē)噪聲感受的問(wèn)題,通過(guò)采集兩輛電動(dòng)汽車(chē)不同工況下的噪聲信號(hào),在計(jì)算得到噪聲客觀參量后,結(jié)合成對(duì)比較法的主觀評(píng)價(jià)試驗(yàn),從多方面對(duì)電動(dòng)汽車(chē)噪聲進(jìn)行了評(píng)價(jià)。根據(jù)主觀評(píng)價(jià)結(jié)果與對(duì)各參量的分析結(jié)果,采用回歸方法建立了電動(dòng)汽車(chē)聲品質(zhì)評(píng)價(jià)結(jié)果與客觀參量間的非線性預(yù)測(cè)模型,并驗(yàn)證了模型的有效性和準(zhǔn)確性。

主題詞:電動(dòng)汽車(chē)聲品質(zhì)主觀評(píng)價(jià)回歸分析

1 前言

汽車(chē)聲品質(zhì)是汽車(chē)NVH問(wèn)題研究的重要內(nèi)容,而建立聲音與人主觀感受之間的數(shù)學(xué)模型是聲品質(zhì)研究的重點(diǎn)之一。隨著聲品質(zhì)研究的深入,人們逐漸形成多維研究的共識(shí),強(qiáng)調(diào)在多個(gè)維度對(duì)聲品質(zhì)進(jìn)行建模。近年來(lái),諸多學(xué)者將多元回歸分析方法運(yùn)用于聲品質(zhì)模型建立并取得很好效果[1~3]。但雖然通過(guò)傳統(tǒng)回歸方法建模簡(jiǎn)單方便,將客觀參量進(jìn)行線性組合即可對(duì)聲品質(zhì)進(jìn)行預(yù)測(cè),具有一定的實(shí)際應(yīng)用價(jià)值,然而由于導(dǎo)致聲品質(zhì)主觀感受差異的因素很多,且這些因素與主觀感受往往成非線性關(guān)系[4],因此應(yīng)建立更精確的非線性模型。

本文通過(guò)測(cè)試再計(jì)算得到噪聲的主觀評(píng)價(jià)值和心理聲學(xué)客觀參量,采用回歸方法建立電動(dòng)汽車(chē)聲品質(zhì)評(píng)價(jià)結(jié)果與客觀參量間的非線性預(yù)測(cè)模型,并對(duì)模型進(jìn)行顯著性檢驗(yàn)以及模型間的比較,驗(yàn)證了模型的有效性和準(zhǔn)確度。

2 車(chē)內(nèi)噪聲主觀評(píng)價(jià)

2.1噪聲信號(hào)的獲取

針對(duì)兩款純電動(dòng)汽車(chē)(A車(chē)和B車(chē))進(jìn)行試驗(yàn),采用LMS公司的SCM05數(shù)據(jù)采集系統(tǒng)采集數(shù)據(jù),采樣頻率為441 00 Hz,采樣時(shí)間為20 s。利用B&K生產(chǎn)的4189A型傳聲器采集噪聲信號(hào),并利用Racelogic車(chē)速儀記錄車(chē)速。

由于駕駛員耳旁聲音具有代表性,故選擇采集該處的噪聲信號(hào)用于分析研究,測(cè)點(diǎn)位置如圖1所示。試驗(yàn)時(shí)分別采集兩款車(chē)以車(chē)速為0、20、40、60 km/h和80 km/h勻速行駛時(shí)駕駛員右耳處的噪聲信號(hào),分別用A1、A2、A3、A4、A5和B1、B2、B3、B4、B5表示。

2.2評(píng)價(jià)試驗(yàn)方法選擇

聲品質(zhì)主觀評(píng)價(jià)是一種通過(guò)聽(tīng)眾對(duì)聲音的主觀感受,按照一定評(píng)價(jià)要求和測(cè)評(píng)規(guī)則對(duì)聲音進(jìn)行評(píng)價(jià)的方法[5],常用方法包括評(píng)分法、成對(duì)比較法、語(yǔ)義細(xì)分法、排序法等。由于成對(duì)比較法在樣本較少且測(cè)評(píng)者測(cè)評(píng)經(jīng)驗(yàn)相對(duì)不足情況下具有明顯優(yōu)勢(shì)[6],因此根據(jù)此次試驗(yàn)?zāi)康暮蜅l件,選擇成對(duì)比較法對(duì)聲音樣品進(jìn)行主觀評(píng)價(jià)。成對(duì)比較法即將聲音樣品進(jìn)行兩兩對(duì)比,作出甲聲音優(yōu)于乙聲音、劣于乙聲音或兩者難以分辨的判斷,最后將對(duì)比結(jié)果進(jìn)行加權(quán)得到數(shù)值化結(jié)果,從而根據(jù)數(shù)值結(jié)果直觀評(píng)價(jià)聲音的優(yōu)劣。該數(shù)值越大,聲品質(zhì)主觀感受越差。

圖1 駕駛員耳旁測(cè)點(diǎn)位置

2.3測(cè)評(píng)者選擇

根據(jù)本次試驗(yàn)要求以及國(guó)內(nèi)外聲品質(zhì)主觀評(píng)價(jià)試驗(yàn)經(jīng)驗(yàn)[7、8],選擇22名聽(tīng)覺(jué)正常的健康成人作為測(cè)評(píng)者進(jìn)行測(cè)評(píng)試驗(yàn),其中男性16名,女性6名,年齡在22到50歲之間。

2.4主觀評(píng)價(jià)

此次主觀評(píng)價(jià)試驗(yàn)以測(cè)評(píng)者對(duì)聲音的偏好性作為評(píng)價(jià)標(biāo)準(zhǔn),測(cè)評(píng)者根據(jù)自己的理解使用成對(duì)比較法對(duì)聲音進(jìn)行判斷。為了達(dá)到效果,10個(gè)噪聲信號(hào)將組成45組比較對(duì)。在進(jìn)行主觀評(píng)價(jià)時(shí),比較對(duì)將被隨機(jī)播放,測(cè)評(píng)者可對(duì)比較對(duì)進(jìn)行任意重復(fù)播放。

3 試驗(yàn)結(jié)果及分析

3.1測(cè)評(píng)數(shù)據(jù)有效性檢驗(yàn)

對(duì)于成對(duì)比較法的有效性檢驗(yàn),首先要檢驗(yàn)每個(gè)測(cè)評(píng)者測(cè)評(píng)結(jié)果的有效性,排除錯(cuò)誤率較大的測(cè)評(píng)結(jié)果。測(cè)評(píng)者通常會(huì)出現(xiàn)的誤判有相同事件比較誤判、不同回放次序的誤判和三角循環(huán)誤判等3種。由于本文所做主觀評(píng)價(jià)試驗(yàn)對(duì)聲音樣品進(jìn)行隨機(jī)單次比較,并且由于第3種誤判是最常見(jiàn)的錯(cuò)誤[9],故僅考慮三角循環(huán)誤判對(duì)試驗(yàn)數(shù)據(jù)正確性的影響。

設(shè)聲事件為i、j、k,利用Pij、Pjk、Pik分別表示比較對(duì)i-j、j-k、i-k的比較結(jié)果,并將任意兩個(gè)聲事件(i,j)的比較結(jié)果i優(yōu)于j、劣于j以及兩者無(wú)差異的結(jié)果分別賦值Pij=1、Pij=-1以及Pij=0。對(duì)i、j、k 3個(gè)聲事件做成對(duì)比較后,三角循環(huán)誤判時(shí)可能出現(xiàn)的結(jié)果歸納為以下3種情況:

a.當(dāng)Pij>0及Pjk≥0時(shí),卻有Pik≤0;

b.當(dāng)Pij=0時(shí),卻有Pik≠Pjk;

c.當(dāng)Pik≤0時(shí),卻有Pik≥0。

采用計(jì)權(quán)一致性系數(shù)來(lái)判斷主觀評(píng)價(jià)數(shù)據(jù)的有效性[10],計(jì)算式如下:

統(tǒng)計(jì)得到各測(cè)評(píng)者所得結(jié)果的計(jì)權(quán)一致性系數(shù)如表1所列。

表1各測(cè)評(píng)者所得結(jié)果的計(jì)權(quán)一致性系數(shù)

為了保證評(píng)價(jià)結(jié)果的可靠性,將排除錯(cuò)誤率較大測(cè)評(píng)者的測(cè)評(píng)結(jié)果,數(shù)據(jù)剔除一般遵循兩個(gè)原則:允許剔除10%~20%的測(cè)評(píng)結(jié)果;計(jì)權(quán)一致性系數(shù)盡量在0.7以上。由表1可知,測(cè)評(píng)結(jié)果中所有測(cè)評(píng)者計(jì)權(quán)一致性系數(shù)都在0.87以上,誤判率極低,表明測(cè)評(píng)者很少出現(xiàn)誤判現(xiàn)象,試驗(yàn)數(shù)據(jù)能正確反映出各測(cè)評(píng)者的真實(shí)感受。

在確定22組測(cè)評(píng)結(jié)果作為有效數(shù)據(jù)后,對(duì)各組數(shù)據(jù)進(jìn)行加權(quán)計(jì)分,得到各測(cè)評(píng)者對(duì)各聲音樣品的評(píng)分結(jié)果。本文采用spss統(tǒng)計(jì)軟件對(duì)評(píng)分結(jié)果進(jìn)行分析,分別計(jì)算22組評(píng)分結(jié)果兩兩之間的相關(guān)系數(shù),再取算數(shù)平均得到測(cè)評(píng)者平均相關(guān)性系數(shù),如表2所列。

由表2可知,Test9、Test10兩組測(cè)評(píng)結(jié)果與其它測(cè)評(píng)結(jié)果相關(guān)性較小,可認(rèn)為是偶然事件,故將其排除。對(duì)剩余20組測(cè)評(píng)結(jié)果進(jìn)行處理,取各組針對(duì)某一工況打分之和作為最終結(jié)果,該結(jié)果即代表大眾對(duì)噪聲作出的正確評(píng)價(jià)。

3.2心理聲學(xué)客觀參量

心理聲學(xué)是研究大腦如何處理進(jìn)入人耳聲音事件的學(xué)科,運(yùn)用于噪聲主觀評(píng)價(jià)可定量反映聽(tīng)覺(jué)感受的差別,消除個(gè)體影響。常用的心理聲學(xué)客觀參量有Zwick?er響度、尖銳度、粗糙度、波動(dòng)度及AI指數(shù)(語(yǔ)音清晰度)等[11]。

表2測(cè)評(píng)者平均相關(guān)性系數(shù)

Zwicker響度更適合描述強(qiáng)度相等、頻率不同聲音的主觀感受;尖銳度是首次特征頻帶分布隨指定響度級(jí)變化關(guān)系的加權(quán),反映高頻成分與低頻成分的比率關(guān)系;波動(dòng)度用來(lái)描述小于20 Hz的低頻噪聲周期性變化特征,粗糙度用來(lái)描述大于20 Hz的高頻噪聲周期性變化特征,兩者都反映聲音大小隨時(shí)間的變化。而AI指數(shù)用來(lái)描述在噪聲環(huán)境下說(shuō)話(huà)的清晰程度,其中,0%表示說(shuō)話(huà)受到很大的干擾,完全聽(tīng)不清;100%表示說(shuō)話(huà)無(wú)干擾,完全能夠聽(tīng)得清楚。SCM05數(shù)據(jù)采集系統(tǒng)能夠采集聲音信號(hào)并對(duì)以上心理聲學(xué)客觀參量進(jìn)行計(jì)算,表3為根據(jù)聲音樣品計(jì)算的客觀參量以及主觀評(píng)分結(jié)果。

表3聲音樣品客觀參量及主觀評(píng)分結(jié)果

3.3試驗(yàn)結(jié)果分析

由表3可知,隨車(chē)速的增加,兩款車(chē)的AI指數(shù)值減小,響度值、A聲壓級(jí)、評(píng)價(jià)得分呈增長(zhǎng)趨勢(shì),說(shuō)明車(chē)速越高,測(cè)評(píng)者對(duì)車(chē)內(nèi)噪聲主觀感受越差;在同一車(chē)速下,大多數(shù)情況下A車(chē)的A聲壓級(jí)、響度均大于B車(chē),A車(chē)AI指數(shù)小于B車(chē),B車(chē)主觀評(píng)價(jià)值小于A車(chē),即B車(chē)給人主觀感受優(yōu)于A車(chē)。

兩車(chē)的尖銳度、粗糙度、波動(dòng)度則隨車(chē)速不同表現(xiàn)出不同的變化趨勢(shì)。

4 聲品質(zhì)模型建立

雖然各心理聲學(xué)參數(shù)可在不同方面描述聲音特點(diǎn),比較聲音樣本優(yōu)劣,并可通過(guò)主觀評(píng)價(jià)試驗(yàn)將主觀感受數(shù)值化而直觀地顯示聲音的優(yōu)劣,但進(jìn)行主觀評(píng)價(jià)試驗(yàn)費(fèi)時(shí)費(fèi)力。為節(jié)省試驗(yàn)時(shí)間和成本,可通過(guò)各心理聲學(xué)參數(shù)與主觀感受之間的關(guān)系建立回歸模型,在遇到同類(lèi)事件時(shí),通過(guò)模型對(duì)評(píng)價(jià)值進(jìn)行預(yù)測(cè)。

回歸分析方法是常用的建模方法。首先選取常用的客觀參量作為解釋變量,并采用逐步線性回歸方法對(duì)聲品質(zhì)評(píng)價(jià)值作回歸分析。F檢驗(yàn)用于檢驗(yàn)被解釋變量與所有解釋變量之間線性關(guān)系的顯著性,在此通過(guò)F檢驗(yàn)對(duì)不顯著變量進(jìn)行剔除,最終得到如下回歸模型(模型1):

式中,Y為電動(dòng)汽車(chē)聲品質(zhì)主觀評(píng)價(jià)值;A為AI指數(shù)值;D為噪聲A計(jì)權(quán)聲壓值。

實(shí)際中客觀參量與主觀感受之間關(guān)系往往是非線性的。為了建立非線性模型,利用了一種改進(jìn)的方法,是顯著的、有意義的。即通過(guò)各客觀參量對(duì)測(cè)評(píng)得分值進(jìn)行曲線估計(jì),探究與主觀評(píng)價(jià)結(jié)果關(guān)系密切的客觀參量以及其變化形式,采用的曲線類(lèi)型包括一元線性、二次函數(shù)、三次函數(shù)、對(duì)數(shù)函數(shù)、指數(shù)函數(shù)等。部分曲線估計(jì)結(jié)果如圖2所示。

由圖2可看出,主觀評(píng)價(jià)值可分別與AI指數(shù)、響度、聲壓級(jí)等構(gòu)成多個(gè)一元多次方程。由于主觀評(píng)價(jià)值與多重因素有關(guān),因此判斷模型與某些客觀參量及其變化形式相關(guān),故嘗試將各客觀參量取二次方、三次方、對(duì)數(shù)等變化后作為解釋變量,將主觀評(píng)價(jià)值作為被解釋變量,再次通過(guò)逐步回歸方法做回歸分析,得到表4所列結(jié)果。表4中非標(biāo)準(zhǔn)化系數(shù)即所求的回歸方程系數(shù),由此得到以下兩個(gè)模型(模型2和模型3):

式中,L為噪聲響度值;S為噪聲尖銳度。

由表4可知,兩個(gè)模型的擬合優(yōu)度R2都大于0.969,說(shuō)明兩個(gè)模型都具有很好的擬合度;T檢驗(yàn)置信概率(sig)都小于0.05,即被解釋變量與每個(gè)解釋變量之間線性關(guān)系都是顯著的。因此,可認(rèn)為兩種模型在統(tǒng)計(jì)學(xué)上是顯著的、有意義的。

圖2 部分曲線估計(jì)結(jié)果

表4回歸方程系數(shù)檢驗(yàn)結(jié)果

同時(shí),為了對(duì)比模型的精度,比較了各模型預(yù)測(cè)評(píng)價(jià)值與試驗(yàn)值的差異,如圖3所示。由圖3可看出,模型3誤差不超過(guò)6.2%,而模型1和模型2最大誤差達(dá)到30%以上,因此認(rèn)為模型3建立的多元非線性模型比模型1和模型2精準(zhǔn)度高。實(shí)際應(yīng)用中,根據(jù)響度值、A聲壓級(jí)及尖銳度即可對(duì)主觀評(píng)價(jià)值進(jìn)行較為準(zhǔn)確的估計(jì),大大節(jié)省了主觀評(píng)價(jià)試驗(yàn)可能花費(fèi)的人力和物力。

圖3 評(píng)價(jià)值與試驗(yàn)值相對(duì)誤差對(duì)比結(jié)果

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

a.對(duì)兩款電動(dòng)汽車(chē)共10個(gè)工況下的車(chē)內(nèi)噪聲進(jìn)行了測(cè)試,得到噪聲客觀參量,從多方面評(píng)價(jià)了汽車(chē)噪聲特點(diǎn)。

b.通過(guò)成對(duì)比較法進(jìn)行了主觀評(píng)價(jià)試驗(yàn),結(jié)果表明,數(shù)值化后的主觀評(píng)價(jià)結(jié)果能夠充分代表人的主觀感受。

c.通過(guò)曲線估計(jì)和線性回歸分析,對(duì)電動(dòng)汽車(chē)聲品質(zhì)測(cè)評(píng)結(jié)果進(jìn)行建模,建立了聲品質(zhì)非線性模型,模型顯示主觀感受量值與客觀參量之間作非線性建模誤差更小,主觀感受主要受響度值、A計(jì)權(quán)聲壓級(jí)、尖銳度等影響。

參考文獻(xiàn)

1陳雙籍,陳瑞石.基于心理聲學(xué)參數(shù)的車(chē)內(nèi)聲品質(zhì)偏好性評(píng)價(jià).噪聲與振動(dòng)控制,2005,3(13):45~47.

2姜吉光,王登峰,蘇麗俐,等.車(chē)內(nèi)噪聲品質(zhì)偏好性主客觀評(píng)價(jià)及相關(guān)性分析.汽車(chē)技術(shù),2012(8):6~10.

3孟祥德,張俊紅,李立順,等.基于回歸分析的車(chē)用柴油機(jī)聲品質(zhì)預(yù)測(cè)技術(shù).內(nèi)燃機(jī)學(xué)報(bào),2011,29(6):524~527.

4徐中明,張芳,周小林.汽車(chē)?yán)嚷暺焚|(zhì)主觀評(píng)價(jià)與分析.汽車(chē)工程,2013;35(2):188~192.

5Otto N,Amman S,Eaton C,et al.Guidelines for jury evalua?tions of automotive sounds.SAE 1999-01-1822,1999.

6Fry J,Jennings P.Jury evaluation of sound quality using neural networks.RASC,2002:301~306.

7Hempel T,Chouard N.Evaluation of interior car sound with a new specific semantic differential design.J.Acous?tic.Am.1999;105(2):1280~1285.

8申秀敏,左曙光,王勇,等.轎車(chē)聲品質(zhì)主觀評(píng)價(jià)研究.噪聲與振動(dòng)控制,2013(1):96~101.

9Stephen Baker,Paul Jennings,Garry Dunne,Roger Wil?liams,Improving the effectiveness of paired comparison tests for automotive sound quality.11thInternational Con?gress and Sound and Vibration,2004.3029~3036.

10毛東興,俞悟周,王佐民.聲品質(zhì)成對(duì)比較主觀評(píng)價(jià)的數(shù)

據(jù)檢驗(yàn)及判據(jù).聲學(xué)學(xué)報(bào),2005(9):46~47.

11毛東興.車(chē)內(nèi)聲品質(zhì)主觀評(píng)價(jià)與分析方法的研究:[學(xué)位論文].上海:同濟(jì)大學(xué),2003.

(責(zé)任編輯文楫)

修改稿收到日期為2015年12月1日。

Evaluation and Modeling of Sound Quality for Electric Vehicle

Hu Teng,Lu Yimin
(Hefei university of technology,Hefei 230009)

【Abstract】Using A-weighted SPL cannot measure the feel of human being on electric vehicle noise accurately.To eliminate this shortcoming,noise signals of two EVs under different conditions are collected.After the noise objective parameters are calculated,and combined with subjective evaluation test of paired comparison method,noise evaluation of EV is made from many aspects.Based on results of the subjective evaluation and analysis of various parameters,a nonlinear forecast model between the evaluation results and the objective parameters of the sound quality is established with regression analysis,and effectiveness and accuracy of this model are verified.

Key words:EV,Sound quality,Subjective evaluation,Regression analysis

中圖分類(lèi)號(hào):U467.7+93

文獻(xiàn)標(biāo)識(shí)碼:A

文章編號(hào):1000-3703(2016)03-0026-04

猜你喜歡
評(píng)價(jià)模型
一半模型
SBR改性瀝青的穩(wěn)定性評(píng)價(jià)
石油瀝青(2021年4期)2021-10-14 08:50:44
中藥治療室性早搏系統(tǒng)評(píng)價(jià)再評(píng)價(jià)
重要模型『一線三等角』
重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
3D打印中的模型分割與打包
FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
基于Moodle的學(xué)習(xí)評(píng)價(jià)
關(guān)于項(xiàng)目后評(píng)價(jià)中“專(zhuān)項(xiàng)”后評(píng)價(jià)的探討
保加利亞轉(zhuǎn)軌20年評(píng)價(jià)
主站蜘蛛池模板: 国产成人精品免费视频大全五级| 亚洲天堂日本| 亚洲免费福利视频| 久久黄色视频影| 午夜日b视频| 91蝌蚪视频在线观看| 国产精品嫩草影院视频| 亚洲精品第五页| 白浆免费视频国产精品视频 | 一级看片免费视频| 中文字幕日韩视频欧美一区| 色噜噜在线观看| 国产网站免费观看| av天堂最新版在线| 亚洲天堂视频在线观看| 少妇极品熟妇人妻专区视频| 免费国产高清精品一区在线| 国产一级毛片高清完整视频版| 久久精品电影| 为你提供最新久久精品久久综合| 91午夜福利在线观看| 久久亚洲高清国产| 综合久久五月天| 一级毛片在线播放免费观看| 久青草国产高清在线视频| 亚洲天堂区| 亚洲性视频网站| 久久亚洲中文字幕精品一区 | 香蕉精品在线| 亚洲三级网站| 日本人妻一区二区三区不卡影院 | 日本一区二区不卡视频| 特级做a爰片毛片免费69| 亚洲无码一区在线观看| 爽爽影院十八禁在线观看| 精品人妻无码中字系列| 婷婷久久综合九色综合88| 日日碰狠狠添天天爽| 999精品视频在线| 国产乱子伦视频在线播放| 国产精品欧美在线观看| 久久婷婷色综合老司机| 日韩av在线直播| 国产午夜精品鲁丝片| 中文国产成人精品久久| 日韩123欧美字幕| 国产精品人人做人人爽人人添| 国产91视频观看| 久久亚洲国产最新网站| 天天躁夜夜躁狠狠躁躁88| 不卡无码h在线观看| 亚洲 欧美 偷自乱 图片| 欧美日韩在线国产| 看av免费毛片手机播放| 天天婬欲婬香婬色婬视频播放| 成人福利在线视频免费观看| 亚洲国语自产一区第二页| jizz国产视频| 99久久国产综合精品2020| 热久久这里是精品6免费观看| 99热在线只有精品| 日韩欧美在线观看| 凹凸国产分类在线观看| 91精品最新国内在线播放| 首页亚洲国产丝袜长腿综合| 成年片色大黄全免费网站久久| www精品久久| 亚洲精品高清视频| 99热精品久久| 激情六月丁香婷婷四房播| 亚洲精品第五页| 呦视频在线一区二区三区| 99re在线免费视频| 香蕉eeww99国产精选播放| 免费A级毛片无码免费视频| 91青青在线视频| 91精品情国产情侣高潮对白蜜| 亚洲中文无码h在线观看| 国产区福利小视频在线观看尤物| 18禁影院亚洲专区| 中文无码精品A∨在线观看不卡 | 天天躁夜夜躁狠狠躁躁88|