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

    參數(shù)優(yōu)化GBDT算法的密封繼電器材質(zhì)識別

    2021-02-15 03:58:14呂冰澤王國濤汪國強
    黑龍江大學自然科學學報 2021年6期
    關(guān)鍵詞:微粒材質(zhì)分類

    呂冰澤, 王國濤,2, 汪國強, 李 碩

    (1.黑龍江大學 電子工程學院, 哈爾濱 150080;2.哈爾濱工業(yè)大學 電器與電子可靠性研究所, 哈爾濱 150001)

    0 引 言

    密封繼電器作為一種常用的航天元器件,在航天設(shè)備中起著重要的作用。因此,密封繼電器的可靠性會直接影響到航天設(shè)備的可靠性[1]。在諸多問題中對密封繼電器可靠性影響最大的就是內(nèi)部多余物問題。繼電器在制作中容易將不易發(fā)現(xiàn)的鋁粒、焊錫和熱熔膠等多余物微?;烊朊芊馇惑w內(nèi)部,當設(shè)備處于失重或者超重的環(huán)境下,會對腔體內(nèi)部的多余物微粒產(chǎn)生激勵,使其產(chǎn)生無規(guī)則運動,微??赡芤苿又两M件連接處或觸點間[2]。當金屬多余物處于電路觸點時易產(chǎn)生電弧破壞內(nèi)部結(jié)構(gòu),而非金屬多余物卡在組件連接處會使得繼電器無法正常運行導致系統(tǒng)失效[3]。

    美國是最早關(guān)注多余物問題和研究微粒碰撞檢測法的國家,國內(nèi)對于多余物的材質(zhì)識別的研究主要集中于哈爾濱工業(yè)大學[4]。文獻[2]首次提出基于時域和頻域特征,實現(xiàn)了對金屬和非金屬多余物材質(zhì)特征識別。文獻[5]提出基于學習矢量量化網(wǎng)絡(luò)(Learning vector quantization, LVQ)和特征級數(shù)據(jù)融合的多余物材質(zhì)識別方法,對繼電器中錫粒、銅絲、橡膠和導線皮多余物進行了材質(zhì)識別,但是其研究的材質(zhì)識別分類器參數(shù)依靠經(jīng)驗確定,材質(zhì)識別模型復雜且泛化能力較差,最終并未應(yīng)用于實際。文獻[6]提出了基于改進型梅爾頻率倒譜系數(shù)(Mel frequency cepstral coefficients, MFCC)的多余物材質(zhì)特征識別方法,建立了多余物材質(zhì)概率統(tǒng)計模型,選用多余物材質(zhì)敏感頻域特征,實現(xiàn)了對導線皮、芯片、鋁和錫粒的有效識別。依據(jù)Mel濾波器設(shè)計原理和多余物信號的頻譜特征,設(shè)計了一種基于敏感頻域的Mel濾波器,提出適用于多余物信號的改進型MFCC特征提取方法,提升了多余物信號的有用信息,但是該方法缺乏實現(xiàn)設(shè)備,并沒有與多余物檢測設(shè)備結(jié)合起來,仍然未應(yīng)用于實際。

    綜上所述,由于目前在對多余物信號進行研究時一般使用傳統(tǒng)信號分析方法或者采用如KNN和支持向量機等經(jīng)典的機器學習算法[7-8],盡管可以對多余物材質(zhì)進行識別分類,但準確率還有待提高。所以,本文提出一種采用參數(shù)優(yōu)化的梯度GBDT算法對密封繼電器中的多余物進行金屬和非金屬微粒材質(zhì)分類。并且在原有特征的基礎(chǔ)上,采用小波變換的方法提取了多余物信號能量特征。為了提高金屬和非金屬材質(zhì)識別精度,加快模型訓練過程中參數(shù)優(yōu)化的效率,引入了Python中的分布式異步超參數(shù)優(yōu)化模塊Hyperopt,使用樹形評估器Parzen對基于默認參數(shù)的GBDT算法進行參數(shù)調(diào)優(yōu)。實驗表明,基于貝葉斯優(yōu)化的GBDT分類算法可以有效提高密封繼電器多余物材質(zhì)識別的準確率。

    1 多余物特征提取

    1.1 傳統(tǒng)特征提取

    受實驗環(huán)境等因素的影響,在進行顆粒碰撞噪聲檢測(Particle impact noise detection,PIND)實驗時,會采集到一部分噪聲信號隨機分布在多余物信號中,需要先去除數(shù)據(jù)中的噪聲[9]。常見的多余物信號特征一般從頻域和時域兩個部分進行提取,在頻域上目前采用短時傅里葉算法來進行多余物特征提取。主要將信號通過窗函數(shù)處理使其由不平穩(wěn)信號加以平穩(wěn),再經(jīng)過傅里葉變換得到局部的頻譜特征值。在時域上,可以直觀地表現(xiàn)出信號的根本特征,本文選取了脈沖面積、脈沖左右對稱度、脈沖上下對稱度、脈沖持續(xù)時間、脈沖上升占比、能量密度、脈沖占比、波峰系數(shù)、面積占比、頻譜質(zhì)心、均方頻率、方差、過零點率和均方根差14個特征。

    1.2 小波變換特征提取

    1.2.1 小波變換

    由于短時傅里葉變換需要對窗函數(shù)進行設(shè)置,但是在實際應(yīng)用中需要對窗函數(shù)的長度只能設(shè)置常見值或者進行多次嘗試,當長度較短時,會導致截取的局部頻率不夠精準;當長度過長時,會出現(xiàn)范圍過廣不能對實時變化情況作出描述的情況,所以需要一種新的方法能精確地將信號的頻譜特征加以表達。小波變換的選擇就是能對這個問題進行解決,小波變換通過對原始信號進行層級分解,在不同的層級根據(jù)所分解的信號進行適當?shù)拇昂瘮?shù)調(diào)整,這樣能夠?qū)π盘柕念l譜信息更好地進行表達[10]。本文在原信號的基礎(chǔ)上,采用小波變換對新的特征加以提取。

    在小波分析中,可以任意選擇小波基函數(shù)。本文對采集信號低頻部分連續(xù)進行3層分解,得到關(guān)于原始信號小波變換的能量特征,小波變換是一種具有低時間分辨率和高頻譜分辨率的快速算法[11]。

    1.2.2 能量特征提取

    小波變換主要通過小波函數(shù)來分析非平穩(wěn)信號,在材質(zhì)識別中,提取金屬和非金屬微粒所包含的隱藏特征。因此,通過數(shù)值分析,計算出小波變換系數(shù)的能量值。

    一維信號X(t)的小波分解[10]為:

    X(t)=Gj+∑Zj

    (1)

    式中:G為低頻分量的近似信號;Z為高頻分量的細節(jié)信號。

    一維信號X(t)與尺度函數(shù)V通過卷積運算可得到低頻的近似系數(shù)cGj(t)。一維信號X(t)與小波函數(shù)通過卷積運算可得到高頻的細節(jié)系數(shù)cGDj(l)。其近似系數(shù)和細節(jié)系數(shù)為:

    (2)

    小波能量系數(shù)為:

    (3)

    式中l(wèi)= 1,2,…,N/2j-1,N是信號長度。

    選擇第3層近似系數(shù)的能量和各層細節(jié)系數(shù)的能量作為信號特征,然后構(gòu)建小波能量系數(shù)特征向量模型為:

    [EGj,EZ1,EZ2,···,EZj]

    (4)

    選擇小波系數(shù)中的雙正交(Biorthogona)小波,分別對金屬和非金屬微粒進行小波3層分解,使用式(2)分別得到1個小波近似系數(shù)和3個細節(jié)系數(shù),使用式(3)提取近似系數(shù)的能量和各層系數(shù)的能量,作為多余物微粒特征。

    2 算法原理

    2.1 GBDT算法

    梯度增強決策樹(GBDT)是在增強算法(Boosting)的基礎(chǔ)上,采用梯度增強算法通過依次迭代生成決策樹,高效地解決二分類問題的方法。梯度提升的訓練目的是為了減少與上一次訓練之間的殘差,通過不斷地在梯度方向上訓練新的模型,通過逐漸降低殘差逼近目標函數(shù),最后達到誤差最小化[11-12]。在GBDT的訓練中,假設(shè)上一輪的迭代得到了強學習器是fm-1(X),損失函數(shù)是L(fm-1(X)),本輪的目標是通過訓練學習找到一個弱學習器hm(X),讓損失函數(shù)L(y,fm(X))=L(fm-1(X))+hm(X))最小,在這里一般定義 GBDT分類的損失函數(shù)為:

    L(y,fm(Xi))=-(yilogρi+(1-yi)log(1-ρi))

    (5)

    GBDT算法流程主要分為以下步驟:

    步驟1:對弱學習器

    (6)

    進行初始化,并確定初始化函數(shù)式(7),利用損失函數(shù)式(5)計算出最佳初始值。

    (7)

    步驟2:對于輸入樣本m=1,2,…,M,建立出M棵決策樹,建立過程如下:

    (a)對每個決策樣本i=1,2,…,M,計算其負梯度,利用式(8)確定每一輪的殘差值。

    (8)

    (b)將(a)中得到的殘差作為訓練數(shù)據(jù),再次對其訓練學習,得到新的決策樹fm(X),其對應(yīng)的葉子節(jié)點區(qū)域為Rjm,j=1,2,…,J,其中J為決策樹m的葉子節(jié)點個數(shù)。

    (c)通過函數(shù)(9)對葉子節(jié)點j=1,2,…,J計算最佳相似度,即擬合值。

    (9)

    該步驟會對葉子節(jié)點的殘差值進行變換,γjm表示第j個葉子節(jié)點殘差轉(zhuǎn)換值,yi表示第j個葉子節(jié)點上樣本xi的類別概率觀測值,fm-1(xi)表示第j個葉子節(jié)點的樣本xi在上一棵樹上的預(yù)測值,γ是常數(shù)(第j個節(jié)點的殘差轉(zhuǎn)換后的值)。對于每個節(jié)點都有一個最優(yōu)殘差轉(zhuǎn)換值γjm是第j個葉子節(jié)點的觀測值yi與fm-1(xi)之間的最小誤差[13]。通過損失函數(shù)對γ微分即可找到每個節(jié)點的最優(yōu)殘差轉(zhuǎn)換值。

    (10)

    (d)對前一輪訓練得到的強學習器進行更新,得到函數(shù)為:

    (11)

    步驟3:通過不斷的迭代學習,得到強學習器為:

    (12)

    在迭代過程中,如果損失誤差小于設(shè)定的閾值,未達到給定的迭代次數(shù),則返回式(12)繼續(xù)學習。然后算出分類的概率值與設(shè)定值對比,并將結(jié)果標記為相應(yīng)的所屬類別。GBDT算法訓練過程如圖1所示,其中T表示輸入數(shù)據(jù)集,F(xiàn)(X)表示輸出梯度提升模型。

    圖1 GBDT算法訓練過程

    2.2 Hyperopt參數(shù)調(diào)優(yōu)

    超參數(shù)優(yōu)化對機器學習算法的必要性已在許多文獻的實驗研究中都得到了證明。大多數(shù)的機器學習算法,如深度神經(jīng)網(wǎng)絡(luò)和梯度提升模型等都需要借助優(yōu)化算法將模型性能調(diào)整至最佳。通常情況下有兩種類型的參數(shù)調(diào)試方法,網(wǎng)格搜索(Grid search)和隨機搜索(Random search)。網(wǎng)格搜索速度慢但是適用于需要隨整個參數(shù)空間進行搜索的情況;隨機搜索速度很快,但是容易遺漏一些重要信息。網(wǎng)格搜索調(diào)參方法在搜索空間方面效果更好,但速度較慢,而隨機搜索很快,但可能會錯過搜索區(qū)域中的重要信息。而Hyperopt使用貝葉斯優(yōu)化的形式進行參數(shù)調(diào)整,可以自適應(yīng)的為給定模型獲得最佳參數(shù),在大范圍內(nèi)優(yōu)化具有數(shù)百個參數(shù)的模型[14]。綜合考慮下,本文使用Python中的Hyperopt工具執(zhí)行貝葉斯優(yōu)化,從而對所建立的模型進行調(diào)參。

    使用Hyperopt主要采用四個步驟:

    步驟1:確定目標函數(shù)。對于所要評估的模型首先需要將其最小化,由于Hyperopt將目標函數(shù)進行“黑箱”處理,所以只是對輸入和輸出值進行處理,使得輸出為最小化的實值,找到損失最小的輸入。

    步驟2:對目標模型中的超參數(shù)劃分不同參數(shù)值的空間范圍。對于貝葉斯優(yōu)化中的每個超參數(shù)來說,空間閾值的概率密度是呈離散分布的。由于不同訓練集之間的最佳模型空間范圍很難確定,因此在參數(shù)優(yōu)化時,采用多個閾值點來設(shè)定每個參數(shù)之間的邊界值。

    步驟3:構(gòu)造代理模型并選擇樹結(jié)構(gòu)估計算法(Tree-structured parzen estimator, TPE)評估超參數(shù)值。

    步驟4:將評估后的損失函數(shù)和最優(yōu)參數(shù)進行輸出。

    由于GBDT算法采用Boosting框架構(gòu)成,其弱分類器為CART回歸樹組合而成,所以此時首先對其弱分類器進行參數(shù)調(diào)優(yōu),對回歸樹進行修剪,在決策樹的訓練學習過程中,有些節(jié)點分類是意義不大的,而這些節(jié)點的存在不僅增加復雜度,而且減少了分類器的正確率。決策樹模型有很多超參數(shù),但不是所有參數(shù)都需要調(diào)試,主要需要調(diào)試樹的最大深度、最大葉子節(jié)點數(shù)和內(nèi)部節(jié)點再劃分所需最小樣本數(shù)這幾個超參數(shù)。

    3 材質(zhì)識別模型構(gòu)建

    基于GBDT算法和Hyperopt模塊的原理介紹構(gòu)建多余物材質(zhì)分類模型,實現(xiàn)對多余物微粒金屬和非金屬的材質(zhì)識別。由于GBDT算法以回歸樹作為基本分類器,因此,弱分類器一般選擇具有偏差較大的回歸樹模型。在迭代訓練生成決策樹時,通過不斷降低偏差來提高弱學習器集成的強學習器的模型精度。同時,預(yù)測模型的因變量是連續(xù)的數(shù)值型變量,因此,選擇平方損失函數(shù)作為模型的損失函數(shù)。

    輸入:T={(x11,…,x17,y1),…,(x11,…,x17,y1)},損失函數(shù)為:

    (13)

    初始化弱分類器:

    對式(13)中的平方損失函數(shù)直接求導,令導數(shù)等于零,得到:

    (14)

    式中:c的取值為所有訓練樣本標簽值的均值。對于m=1,2,…,M,對樣本i=1,2,…,n,根據(jù)式(8)計算損失函數(shù)的負梯度值,利用數(shù)據(jù)集(xi,rmi)擬合下一輪基礎(chǔ)模型。得到第m棵決策樹的葉子節(jié)點,并計算各葉子節(jié)點的最佳擬合值。

    (15)

    式(15)表示cmj的取值為第m棵樹的第j個結(jié)點中殘差的均值。

    更新式(16),得:

    (16)

    得到新一輪函數(shù)值,結(jié)合前(m-1)輪的基礎(chǔ)模型,加入學習率防止模型過擬合,得到梯度提升模型函數(shù)為:

    FM(x)=FM-1(x)+vfm(x)

    (17)

    式中v為學習率。

    為了提高識別準確率,采用默認參數(shù)建立基于GBDT多余物金屬非金屬分類模型后,使用Hyperopt對其進行參數(shù)優(yōu)化提高模型整體效率[15]。為克服數(shù)據(jù)集自身造成的缺陷,還將采用k折交叉檢驗的方法。通過將整體數(shù)據(jù)集分為k份,在對模型進行訓練時,將其中K-1份進行訓練,然后將剩下的一份數(shù)據(jù)集作為數(shù)據(jù)集進行測試,對上述操作重分k次,每次采用不同份的數(shù)據(jù)集進行實驗,最終將k次檢驗的結(jié)果求均值輸出作為最終結(jié)果。

    4 實驗分析

    4.1 多余物材質(zhì)數(shù)據(jù)集

    選擇來自某型號密封繼電器采集到的一組數(shù)據(jù)集來評估該材質(zhì)識別模型的有效性,采集到的主要數(shù)據(jù)集信息如表1所示,其中包括了樣本總數(shù),金屬微粒數(shù)量和非金屬微粒數(shù)量,以及數(shù)據(jù)集選擇的特征總數(shù)。

    表1 數(shù)據(jù)集的主要信息

    為了識別試驗效果,對處理后的數(shù)據(jù)集采用參數(shù)優(yōu)化的GBDT算法進行分類驗證。在進行實驗時,70%的數(shù)據(jù)用于訓練GBDT模型;30%的數(shù)據(jù)用來預(yù)測分類模型,評估模型的識別效果。

    通常情況下,信號的特征信息越豐富,對于模型識別精度提高就越明顯。在原有脈沖面積、脈沖左右對稱度、脈沖上下對稱度、脈沖持續(xù)時間、脈沖上升占比、能量密度、脈沖占比、波峰系數(shù)、面積占比、頻譜質(zhì)心、均方頻率、方差、過零點率和均方根差14個特征的基礎(chǔ)上,對原始信號進行小波變換得到信號頻帶分析,將計算得到的每個頻帶能量作為新特征。與本文原有特征相結(jié)合組成新的樣本數(shù)據(jù)集進行了分析,提取出的部分特征如表2所示。

    表2 小波變換部分特征集

    4.2 評價指標

    對于機器學習分類算法,存在樣本過擬合問題。為了更好地對模型性能加以反映,在保證了準確率不變的基礎(chǔ)上使用混淆矩陣評估模型調(diào)參的效果,如表3所示。

    表3 二分類混淆矩陣

    在對該模型進行評估時,經(jīng)過各個評價指標對比,最后選擇查準率(Precison)、召回率(Recall)和F1-score三個值作為評價指標。

    “查準率”是用來表示預(yù)測集中被預(yù)測正確的個數(shù),定義為:

    (18)

    召回率主要是查找有多少樣本能被找回公式為:

    (19)

    F1-score是綜合評價查準率和查全率的評價指標,為查準率和查全率的調(diào)和平均數(shù)[14]。計算方法為:

    (20)

    式中,β=1時,F(xiàn)1-score值一般情況下β=1。F1-score值能很好地反映二分類的分類效果。

    4.3 實驗驗證

    為了驗證加入小波變換能量特征后的整體特征集的識別效果,使用原有特征組合進行機器學習模型識別,結(jié)果如表4所示。

    表4 原特征集測試結(jié)果

    此時,模型的精確度為0.90,使用加入小波變換后的特征集采用同樣的模型進行金屬和非金屬識別,得到的測試結(jié)果如表5所示。

    表5 新特征集測試結(jié)果

    加入新特征后模型精確度為0.914,準確率提高1.4%。此外,由表4和表5中的判別指標對比可以發(fā)現(xiàn),對非金屬的檢測效果有所提升,證明加入小波變換的能量特征可以提高多余物材質(zhì)分類精度。

    本實驗采用GBDT算法進行密封繼電器中金屬和非金屬多余物材質(zhì)識別,由于GBDT算法是基于Boosting框架下的集成算法,由弱分類器回歸樹(CART)組成,主要選擇貝葉斯優(yōu)化法對弱分類器中的最大深度(max_depth)和最大特征數(shù)(max_features)進行單棵樹的優(yōu)化。在Boosting框架下,n_estimators表示整個模型弱分類個數(shù),在實際參數(shù)優(yōu)化過程中,一般與learning_rate,即每個弱分類器的權(quán)重縮減系數(shù)共同進行調(diào)優(yōu)。

    在訓練集上將learning_rate設(shè)置為0.5,并采用5重交叉驗證對參數(shù)n_estimators進行尋優(yōu),結(jié)果為800。再用同樣的方法對max_depth和max_features進行超參數(shù)的選擇,參數(shù)值均為5。

    在進行參數(shù)選擇后,基于前面的模型構(gòu)建,利用數(shù)據(jù)集分別在優(yōu)化前后的GBDT模型上進行測試。此外,為了分析GBDT模型的分類效果,分別以KNN分類模型和SVM支持向量機模型測試結(jié)果進行對比。采用默認參數(shù)進行測試得到準確率為0.902,采用Hyperopt對其進行參數(shù)優(yōu)化,得到整體模型參數(shù)優(yōu)化數(shù)值,如表6所示。

    表6 模型參數(shù)優(yōu)化

    參數(shù)優(yōu)化后與優(yōu)化前模型對比測試結(jié)果如表7所示。

    表7 默認參數(shù)測試結(jié)果

    經(jīng)過參數(shù)優(yōu)化后,模型整體準確率可以達到93.2%,比默認參數(shù)模型準確率提升了3.1%,此外,其他評估指標也有一定的提高。將優(yōu)化后的模型與目前常用的KNN分類模型和SVM支持向量機模型進行對比,分類準確率對比如表8所示。

    表8 不同模型的分類準確率

    由預(yù)測結(jié)果可以看出,參數(shù)優(yōu)化后的GBDT算法的分類準確率更高,準確率為0.932。分類結(jié)果優(yōu)于KNN、支持向量機模型和基于默認參數(shù)的GBDT分類模型。因此,參數(shù)優(yōu)化后的GBDT模型為分類結(jié)果最佳的選擇。由表9可以看出,與默認參數(shù)的GBDT算法相比,使用參數(shù)優(yōu)化的GBDT算法對金屬和非金屬微粒識別時,各項指標都有提升,這相比傳統(tǒng)的檢測方法的精確度也有一定程度的提升,尤其對非金屬識別的準確度提升較高。雖然使用SVM支持向量機算法的二分類各項指標也能很好地反映金屬微粒和非金屬微粒分類效果,但由于泛化性較低,因此,使用參數(shù)優(yōu)化的GBDT算法能夠有效地提升金屬微粒和非金屬微粒的準確度。

    表9 評估指標值

    5 結(jié) 論

    采用貝葉斯參數(shù)優(yōu)化的GBDT算法對航天密封繼電器內(nèi)部的金屬和非金屬微粒進行材質(zhì)識別。在原有特征的基礎(chǔ)上,利用小波變換提取了能量系數(shù)作為特征值。實驗驗證,加入新特征后模型識別準確率提升了1.4%。根據(jù)GBDT算法流程建立了材質(zhì)識別的分類器模型,并使用Hyperopt模塊的TPE評估器對多余物微粒識別模型進行超參數(shù)調(diào)優(yōu)。實驗結(jié)果表明,調(diào)參后的GBDT算法具有更好的分類效果,分類準確率提升到93%,提升了3.1%。同時,3個評價指標相對于KNN分類模型和SVM支持向量機模型也有一定的提升,其中精確度、召回率和F1-score值均提升超過了3%。此方案可以提高密封繼電器金屬和非金屬微粒的識別準確率。通過對多余物的材質(zhì)信息更準確的識別,還可以改善多余物的生產(chǎn)環(huán)節(jié),進一步推動相關(guān)材質(zhì)的改進,提高整個航天系統(tǒng)的可靠性。

    猜你喜歡
    微粒材質(zhì)分類
    塑料微粒的旅程
    塑料微粒的旅程
    塑料微粒的旅程
    分類算一算
    分類討論求坐標
    衣服材質(zhì)大揭秘
    數(shù)據(jù)分析中的分類討論
    教你一招:數(shù)的分類
    致今天的你,致年輕的你
    外套之材質(zhì)對比戰(zhàn)
    Coco薇(2015年10期)2015-10-19 12:17:50
    成人av在线播放网站| 天堂中文最新版在线下载 | 亚洲av国产av综合av卡| 高清av免费在线| 欧美性感艳星| 久久人人爽人人片av| 女人被狂操c到高潮| 亚洲国产色片| 免费看美女性在线毛片视频| 亚洲四区av| 亚洲av在线观看美女高潮| 一个人看视频在线观看www免费| 国内少妇人妻偷人精品xxx网站| 免费不卡的大黄色大毛片视频在线观看 | 欧美日本视频| 亚洲av中文av极速乱| 精品人妻视频免费看| 毛片女人毛片| 亚洲国产精品专区欧美| 日本色播在线视频| 国产真实伦视频高清在线观看| 少妇丰满av| av免费在线看不卡| 日韩欧美三级三区| 国产黄a三级三级三级人| 夫妻午夜视频| 国产av国产精品国产| 日韩欧美精品免费久久| 男女视频在线观看网站免费| 午夜精品在线福利| 99久国产av精品| 可以在线观看毛片的网站| 老司机影院毛片| 久久久久久九九精品二区国产| 最新中文字幕久久久久| 国产综合精华液| 啦啦啦啦在线视频资源| 亚洲色图av天堂| 国产一区二区在线观看日韩| 日本爱情动作片www.在线观看| 成年免费大片在线观看| 日韩国内少妇激情av| 亚洲欧美精品专区久久| 一夜夜www| 国产视频内射| 中文在线观看免费www的网站| 日韩视频在线欧美| 亚洲精品乱码久久久久久按摩| 少妇的逼好多水| 精品欧美国产一区二区三| 丝袜喷水一区| 精品国产露脸久久av麻豆 | 日本av手机在线免费观看| 欧美xxxx黑人xx丫x性爽| 高清毛片免费看| 99久久九九国产精品国产免费| 精华霜和精华液先用哪个| 免费看美女性在线毛片视频| 国产一级毛片在线| 成人欧美大片| 高清午夜精品一区二区三区| 一级毛片 在线播放| 欧美xxxx黑人xx丫x性爽| 欧美丝袜亚洲另类| 少妇熟女aⅴ在线视频| 91av网一区二区| 91在线精品国自产拍蜜月| 久热久热在线精品观看| 人妻制服诱惑在线中文字幕| 日本午夜av视频| 两个人视频免费观看高清| 亚洲av电影不卡..在线观看| 女人久久www免费人成看片| 日韩伦理黄色片| 一二三四中文在线观看免费高清| 欧美3d第一页| 午夜视频国产福利| 一级毛片黄色毛片免费观看视频| 国产午夜精品久久久久久一区二区三区| 日韩欧美三级三区| 国产精品美女特级片免费视频播放器| 联通29元200g的流量卡| 久久精品国产亚洲av涩爱| 午夜福利在线在线| 免费观看无遮挡的男女| 韩国av在线不卡| 可以在线观看毛片的网站| 三级经典国产精品| 爱豆传媒免费全集在线观看| 国产片特级美女逼逼视频| 成人一区二区视频在线观看| 男人舔奶头视频| 狂野欧美激情性xxxx在线观看| 少妇熟女欧美另类| www.av在线官网国产| 又大又黄又爽视频免费| 中文乱码字字幕精品一区二区三区 | 国产毛片a区久久久久| 少妇被粗大猛烈的视频| 国产亚洲精品久久久com| 国产精品久久久久久久久免| 成人欧美大片| 亚洲熟女精品中文字幕| 18禁裸乳无遮挡免费网站照片| 国产一区有黄有色的免费视频 | 男女视频在线观看网站免费| 最近最新中文字幕大全电影3| 国产一级毛片七仙女欲春2| 国产成人福利小说| 极品教师在线视频| 免费观看av网站的网址| 偷拍熟女少妇极品色| 久久草成人影院| 欧美不卡视频在线免费观看| 亚洲国产精品专区欧美| 亚洲高清免费不卡视频| 六月丁香七月| 欧美丝袜亚洲另类| 国产熟女欧美一区二区| 精品国产三级普通话版| 中国国产av一级| 免费看日本二区| 亚洲色图av天堂| 美女黄网站色视频| 欧美丝袜亚洲另类| 亚洲色图av天堂| 97超视频在线观看视频| 女人被狂操c到高潮| h日本视频在线播放| 99久久中文字幕三级久久日本| 午夜精品在线福利| 国产午夜精品一二区理论片| 国产精品无大码| 天美传媒精品一区二区| 老司机影院成人| 国产午夜精品一二区理论片| 亚洲精品成人久久久久久| 成人无遮挡网站| 成年女人在线观看亚洲视频 | av在线观看视频网站免费| 99九九线精品视频在线观看视频| 精品99又大又爽又粗少妇毛片| 色网站视频免费| xxx大片免费视频| 少妇熟女欧美另类| 国产有黄有色有爽视频| 听说在线观看完整版免费高清| 99久国产av精品| 99视频精品全部免费 在线| 中文精品一卡2卡3卡4更新| 最近的中文字幕免费完整| 99re6热这里在线精品视频| 精品国产三级普通话版| 日韩亚洲欧美综合| 亚洲国产日韩欧美精品在线观看| 国内揄拍国产精品人妻在线| 久久国内精品自在自线图片| 看黄色毛片网站| 久久精品国产自在天天线| 久久久久九九精品影院| 高清视频免费观看一区二区 | 亚洲丝袜综合中文字幕| 亚洲乱码一区二区免费版| 91久久精品国产一区二区成人| 午夜久久久久精精品| 国内少妇人妻偷人精品xxx网站| 日本免费a在线| 直男gayav资源| 国产成人精品久久久久久| 国产亚洲91精品色在线| 国产一级毛片在线| 中文天堂在线官网| 男女视频在线观看网站免费| 啦啦啦韩国在线观看视频| 黄片无遮挡物在线观看| 一区二区三区四区激情视频| 亚洲av不卡在线观看| 伦理电影大哥的女人| xxx大片免费视频| 麻豆av噜噜一区二区三区| 久久久久久久久久久免费av| 国产国拍精品亚洲av在线观看| 少妇人妻精品综合一区二区| 白带黄色成豆腐渣| 亚洲精品久久久久久婷婷小说| 视频中文字幕在线观看| 免费观看精品视频网站| 国产不卡一卡二| 男人和女人高潮做爰伦理| 国产欧美日韩精品一区二区| 亚洲激情五月婷婷啪啪| 午夜福利在线观看吧| 欧美97在线视频| 亚洲欧美日韩无卡精品| av女优亚洲男人天堂| 内射极品少妇av片p| 国产精品一及| 嘟嘟电影网在线观看| 你懂的网址亚洲精品在线观看| 成人特级av手机在线观看| 狂野欧美白嫩少妇大欣赏| 特级一级黄色大片| 久久久精品94久久精品| 日日啪夜夜爽| 97人妻精品一区二区三区麻豆| av国产久精品久网站免费入址| 久久久久久久久中文| 九九久久精品国产亚洲av麻豆| 国产成人aa在线观看| 亚洲欧美日韩东京热| 国产亚洲91精品色在线| av黄色大香蕉| 日韩伦理黄色片| 97热精品久久久久久| videos熟女内射| 欧美日韩综合久久久久久| 亚洲av成人av| 国产在线一区二区三区精| 边亲边吃奶的免费视频| 久久久久久久午夜电影| 男女边吃奶边做爰视频| 日本爱情动作片www.在线观看| 干丝袜人妻中文字幕| 你懂的网址亚洲精品在线观看| 欧美一区二区亚洲| 日韩欧美一区视频在线观看 | 精品国产露脸久久av麻豆 | 精品国内亚洲2022精品成人| 国产不卡一卡二| 亚洲成人久久爱视频| 三级经典国产精品| 亚洲精品久久午夜乱码| 最近最新中文字幕大全电影3| 日本黄大片高清| 精品一区二区三区视频在线| 久久久久久久久中文| 亚洲综合色惰| 男人狂女人下面高潮的视频| 国产成人精品久久久久久| 两个人视频免费观看高清| 成人欧美大片| 国产精品人妻久久久影院| 男人狂女人下面高潮的视频| 国内精品美女久久久久久| 女人十人毛片免费观看3o分钟| a级一级毛片免费在线观看| 久久久成人免费电影| 国产精品国产三级国产专区5o| 肉色欧美久久久久久久蜜桃 | 国产免费视频播放在线视频 | 日韩成人av中文字幕在线观看| 亚洲人成网站高清观看| 九九久久精品国产亚洲av麻豆| 精品不卡国产一区二区三区| 久久久久久伊人网av| 高清欧美精品videossex| 免费观看精品视频网站| freevideosex欧美| 91久久精品国产一区二区成人| 乱码一卡2卡4卡精品| 人人妻人人澡人人爽人人夜夜 | 一个人免费在线观看电影| 一级av片app| av天堂中文字幕网| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲av成人av| 精品欧美国产一区二区三| 可以在线观看毛片的网站| 黄片wwwwww| 国产在视频线精品| 国产伦精品一区二区三区四那| 成人毛片60女人毛片免费| 一级片'在线观看视频| 老司机影院成人| 小蜜桃在线观看免费完整版高清| 亚洲av国产av综合av卡| 亚洲av电影在线观看一区二区三区 | 91久久精品国产一区二区三区| 欧美xxxx性猛交bbbb| 日韩在线高清观看一区二区三区| 2022亚洲国产成人精品| 欧美zozozo另类| 直男gayav资源| 国产老妇女一区| 国产精品国产三级国产av玫瑰| 亚洲国产色片| 国产一区二区三区综合在线观看 | 国产精品一二三区在线看| 午夜福利在线观看免费完整高清在| 亚洲国产精品成人综合色| 欧美高清成人免费视频www| 亚洲精品色激情综合| 蜜桃久久精品国产亚洲av| 久久99精品国语久久久| 午夜激情福利司机影院| 久久久久久九九精品二区国产| 亚洲综合精品二区| 成人亚洲欧美一区二区av| 国产av国产精品国产| kizo精华| 久久99热这里只有精品18| 久99久视频精品免费| 国产一区二区三区综合在线观看 | 九色成人免费人妻av| 国产 一区 欧美 日韩| 亚洲欧洲国产日韩| 中文字幕制服av| 成人毛片60女人毛片免费| 69av精品久久久久久| 亚洲欧美成人精品一区二区| 国产成人福利小说| 欧美潮喷喷水| 最近最新中文字幕免费大全7| 免费黄网站久久成人精品| 2021天堂中文幕一二区在线观| av女优亚洲男人天堂| 欧美日韩一区二区视频在线观看视频在线 | 丰满乱子伦码专区| 国产伦精品一区二区三区四那| av又黄又爽大尺度在线免费看| 高清在线视频一区二区三区| 草草在线视频免费看| 亚洲精品456在线播放app| 午夜福利在线在线| 青春草国产在线视频| 精品一区二区三区人妻视频| 亚洲va在线va天堂va国产| 国产精品.久久久| 亚洲内射少妇av| 亚洲国产成人一精品久久久| 亚洲欧美精品专区久久| 日韩一区二区三区影片| 日韩精品青青久久久久久| 在线 av 中文字幕| 熟女电影av网| videossex国产| 一个人看的www免费观看视频| 国产av国产精品国产| 熟妇人妻久久中文字幕3abv| a级毛片免费高清观看在线播放| 成人一区二区视频在线观看| 五月玫瑰六月丁香| 啦啦啦韩国在线观看视频| 国产精品人妻久久久影院| 在线观看人妻少妇| 亚洲18禁久久av| 乱码一卡2卡4卡精品| 日本wwww免费看| 99热6这里只有精品| 欧美xxⅹ黑人| 在现免费观看毛片| 精品一区二区三区人妻视频| 中文乱码字字幕精品一区二区三区 | 色吧在线观看| 国产乱来视频区| av在线蜜桃| 国模一区二区三区四区视频| 国产老妇女一区| 日韩欧美精品免费久久| 国产熟女欧美一区二区| 日本欧美国产在线视频| 身体一侧抽搐| 乱系列少妇在线播放| av网站免费在线观看视频 | 九九在线视频观看精品| 午夜免费激情av| 色网站视频免费| 日韩欧美精品免费久久| 亚洲欧美一区二区三区黑人 | 深夜a级毛片| 超碰97精品在线观看| 亚洲乱码一区二区免费版| 春色校园在线视频观看| 黄色一级大片看看| 亚洲国产精品成人综合色| 白带黄色成豆腐渣| 午夜老司机福利剧场| 人体艺术视频欧美日本| 99久国产av精品国产电影| 亚洲电影在线观看av| 97热精品久久久久久| 国产成人freesex在线| 久久久精品94久久精品| 99久久精品一区二区三区| 国产黄a三级三级三级人| 美女被艹到高潮喷水动态| 天堂√8在线中文| 精品国内亚洲2022精品成人| 国国产精品蜜臀av免费| 亚洲无线观看免费| 女人十人毛片免费观看3o分钟| 九草在线视频观看| 日韩av在线大香蕉| av天堂中文字幕网| 男人和女人高潮做爰伦理| 欧美精品一区二区大全| 亚洲精品第二区| 亚洲av中文av极速乱| a级毛色黄片| 免费看a级黄色片| 女人被狂操c到高潮| 日本午夜av视频| 亚洲久久久久久中文字幕| 亚洲国产精品专区欧美| 免费黄频网站在线观看国产| 天堂√8在线中文| 80岁老熟妇乱子伦牲交| 偷拍熟女少妇极品色| 国产一区二区三区综合在线观看 | 蜜桃久久精品国产亚洲av| 亚洲国产av新网站| 欧美一级a爱片免费观看看| 九色成人免费人妻av| 小蜜桃在线观看免费完整版高清| 免费看美女性在线毛片视频| 精品酒店卫生间| 欧美成人一区二区免费高清观看| 你懂的网址亚洲精品在线观看| 日本欧美国产在线视频| 亚洲精品视频女| 色视频www国产| 五月伊人婷婷丁香| 天堂影院成人在线观看| 国产黄a三级三级三级人| 欧美高清成人免费视频www| 国产综合精华液| 九草在线视频观看| 最近手机中文字幕大全| 欧美精品一区二区大全| 久久精品久久久久久久性| 美女大奶头视频| 久久久久久久国产电影| 国产精品一及| 日韩欧美一区视频在线观看 | 99视频精品全部免费 在线| 在现免费观看毛片| 亚洲精品久久久久久婷婷小说| 成人二区视频| 美女国产视频在线观看| 夜夜爽夜夜爽视频| 国产视频首页在线观看| 大香蕉久久网| 婷婷色综合www| 国内少妇人妻偷人精品xxx网站| 国产精品国产三级国产专区5o| 哪个播放器可以免费观看大片| 别揉我奶头 嗯啊视频| 国产成年人精品一区二区| 免费观看精品视频网站| 亚州av有码| 亚洲图色成人| 高清av免费在线| 尾随美女入室| 18禁在线无遮挡免费观看视频| 你懂的网址亚洲精品在线观看| 人妻少妇偷人精品九色| 国产真实伦视频高清在线观看| 日韩强制内射视频| a级一级毛片免费在线观看| 一级毛片久久久久久久久女| 夫妻性生交免费视频一级片| 成人性生交大片免费视频hd| 自拍偷自拍亚洲精品老妇| 亚洲人成网站在线播| 久久久久久久久久黄片| 国产午夜精品一二区理论片| 国产av国产精品国产| av福利片在线观看| 好男人在线观看高清免费视频| 中文天堂在线官网| 麻豆成人午夜福利视频| 国产精品综合久久久久久久免费| 日本欧美国产在线视频| 免费看a级黄色片| 熟妇人妻不卡中文字幕| 五月玫瑰六月丁香| 国产午夜精品论理片| 亚洲高清免费不卡视频| 人妻夜夜爽99麻豆av| 国产精品一二三区在线看| 日韩不卡一区二区三区视频在线| 联通29元200g的流量卡| 视频中文字幕在线观看| 亚洲成人av在线免费| 中文乱码字字幕精品一区二区三区 | 禁无遮挡网站| 久久久午夜欧美精品| 蜜桃久久精品国产亚洲av| 亚洲精品乱久久久久久| 在线播放无遮挡| 一个人观看的视频www高清免费观看| 免费观看在线日韩| 国产亚洲一区二区精品| 国内揄拍国产精品人妻在线| 国产亚洲av片在线观看秒播厂 | 日本黄色片子视频| 亚洲av.av天堂| 午夜免费观看性视频| 日本免费a在线| av在线亚洲专区| 国产综合精华液| 国产午夜精品论理片| 听说在线观看完整版免费高清| 91狼人影院| 精品久久久久久久末码| 男女边摸边吃奶| 18禁动态无遮挡网站| 国产精品人妻久久久影院| 国产伦理片在线播放av一区| 日韩欧美精品免费久久| 国产精品国产三级专区第一集| 欧美高清性xxxxhd video| 亚洲内射少妇av| 日韩,欧美,国产一区二区三区| 一本久久精品| 久久99精品国语久久久| 中文在线观看免费www的网站| 国产亚洲最大av| 伦理电影大哥的女人| av女优亚洲男人天堂| 美女主播在线视频| 精品久久久噜噜| 天美传媒精品一区二区| 寂寞人妻少妇视频99o| 欧美丝袜亚洲另类| av免费在线看不卡| 91狼人影院| 午夜免费激情av| 日本与韩国留学比较| 国产综合懂色| 乱人视频在线观看| 全区人妻精品视频| 亚洲一区高清亚洲精品| 久久久久久久久久黄片| 亚洲在久久综合| 婷婷色综合www| 国产午夜精品久久久久久一区二区三区| videos熟女内射| 成人漫画全彩无遮挡| 最近中文字幕2019免费版| 人妻系列 视频| 七月丁香在线播放| 91狼人影院| 毛片女人毛片| 亚洲精品成人久久久久久| 国产探花在线观看一区二区| 免费观看无遮挡的男女| 天天躁夜夜躁狠狠久久av| 欧美成人午夜免费资源| 国产欧美另类精品又又久久亚洲欧美| 久久这里有精品视频免费| 亚洲av一区综合| 国产精品久久久久久久电影| 婷婷色av中文字幕| 天堂√8在线中文| 国产精品一区二区三区四区免费观看| 水蜜桃什么品种好| 国产伦精品一区二区三区视频9| 国产伦精品一区二区三区四那| 毛片一级片免费看久久久久| 日韩av不卡免费在线播放| 亚洲久久久久久中文字幕| 男女国产视频网站| 熟女电影av网| 嘟嘟电影网在线观看| 久久久国产一区二区| 国产精品.久久久| 亚洲精品第二区| 免费电影在线观看免费观看| 国产伦理片在线播放av一区| 国产成人a∨麻豆精品| 国产成人精品久久久久久| 国产成人aa在线观看| 免费高清在线观看视频在线观看| 国产精品伦人一区二区| 看非洲黑人一级黄片| 国产乱人偷精品视频| 美女黄网站色视频| 国产淫语在线视频| 久久久久久国产a免费观看| 国产精品人妻久久久久久| 一区二区三区四区激情视频| .国产精品久久| 久久国内精品自在自线图片| 一区二区三区高清视频在线| 非洲黑人性xxxx精品又粗又长| 日韩欧美精品v在线| 如何舔出高潮| www.色视频.com| 久久精品熟女亚洲av麻豆精品 | 国产亚洲一区二区精品| 国产不卡一卡二| 又黄又爽又刺激的免费视频.| 97在线视频观看| 日韩,欧美,国产一区二区三区| 能在线免费看毛片的网站| 我的老师免费观看完整版| 99热这里只有是精品在线观看| 免费播放大片免费观看视频在线观看| 男女那种视频在线观看| 淫秽高清视频在线观看| 黄色一级大片看看| 大香蕉97超碰在线| 91午夜精品亚洲一区二区三区| 晚上一个人看的免费电影| 美女黄网站色视频| 看黄色毛片网站| 丝袜美腿在线中文| 午夜日本视频在线| 国产一区有黄有色的免费视频 | 伦精品一区二区三区| 成人性生交大片免费视频hd| 亚洲av男天堂| 不卡视频在线观看欧美| 午夜激情欧美在线|