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

    自適應(yīng)鄰域構(gòu)造流形學習算法及故障降維診斷*

    2017-01-09 05:49:08張曉濤唐力偉鄧士杰
    振動、測試與診斷 2016年6期
    關(guān)鍵詞:構(gòu)造方法流形降維

    張曉濤, 唐力偉, 王 平, 鄧士杰

    (軍械工程學院火炮工程系 石家莊,050003)

    ?

    自適應(yīng)鄰域構(gòu)造流形學習算法及故障降維診斷*

    張曉濤, 唐力偉, 王 平, 鄧士杰

    (軍械工程學院火炮工程系 石家莊,050003)

    針對流形學習算法中近鄰構(gòu)造問題,提出一種自適應(yīng)鄰域構(gòu)造方法,該方法基于馬氏距離計算樣本間相似系數(shù),由相似系數(shù)均值確定初始近鄰數(shù),根據(jù)樣本高斯核概率密度估計調(diào)整近鄰數(shù),并將自適應(yīng)鄰域構(gòu)造方法用于改進的主成分分析聯(lián)合局部保持投景(principal component analysis-locality preserving projections,簡稱PCA-LPP)流形學習算法中。通過齒輪箱故障類型識別對其特征降維性能進行驗證,結(jié)果表明,自適應(yīng)鄰域PCA-LPP方法比傳統(tǒng)的k近鄰方法及原始無處理的特征識別率都高,可以達到94.67%。

    自適應(yīng); 鄰域; 高斯核概率估計; 流形學習; 故障診斷

    引 言

    流形學習算法是一類非常有效的非線性數(shù)據(jù)降維算法[1],在很多具體的流形學習算法實現(xiàn)中,諸如局部保持投影(locality preserving projection, 簡稱LPP)[2]、拉普拉斯映射(Laplacian eigenmap,簡稱LE)[3]都需要先構(gòu)造數(shù)據(jù)樣本近鄰圖。傳統(tǒng)的近鄰圖構(gòu)造主要是k近鄰和ε近鄰兩種方法[4],其主要缺點在于各個樣本的近鄰數(shù)需要人為經(jīng)驗確定,且所有樣本采用相同的近鄰數(shù),不能針對樣本自身的鄰域分布情況進行數(shù)據(jù)調(diào)整。自適應(yīng)鄰域的思想主要是基于樣本周圍一定范圍數(shù)據(jù)分布的情況,按照預設(shè)規(guī)則自動確定鄰域大小。目前常用的自適應(yīng)鄰域構(gòu)造思想主要基于歐氏距離或者余弦夾角距離計算相似系數(shù),并根據(jù)相似系數(shù)的分布設(shè)定閾值自適應(yīng)求取近鄰數(shù),這樣的方法多用于圖像處理及人臉識別研究中。Yang等[5]基于歐氏距離指數(shù)函數(shù)計算樣本相似系數(shù)。劉鳳連等[6]基于圖像歐氏距離導數(shù)計算樣本相似系數(shù)。黃璞等[7]采用樣本夾角余弦距離作為相似系數(shù)。三種方法在人臉識別研究中取得了良好的分類效果。李城梁等[8]將基于樣本切空間距離的自適應(yīng)鄰域方法應(yīng)用于機械故障信號特征提取,故障分類正確率得到明顯提升。

    自適應(yīng)鄰域構(gòu)造方法能夠有效避免人為選擇近鄰范圍的隨意性,具有更好的樣本局部非線性流形特征表達能力。筆者提出一種自適應(yīng)鄰域構(gòu)造方法,基于馬氏距離相似系數(shù)均值及樣本高斯核概率密度估計調(diào)整自適應(yīng)獲得樣本近鄰數(shù),并將自適應(yīng)鄰域構(gòu)造方法應(yīng)用于PCA-LPP流形學習算法,通過齒輪箱故障信號特征向量降維識別驗證了自適應(yīng)鄰域PCA-LPP方法的有效性。

    1 自適應(yīng)鄰域構(gòu)造方法

    設(shè)高維空間RD中存在數(shù)據(jù)集X=[x1,x2,…,xN],基于近鄰圖構(gòu)造相似矩陣S能夠用于多種流形學習算法,實現(xiàn)數(shù)據(jù)的降維處理。

    傳統(tǒng)的k近鄰構(gòu)造方法,基于樣本間的歐氏距離計算樣本相似性,每個樣本近鄰值相同,相似系數(shù)具體表達如下

    (1)

    k近鄰法存在兩點不足:a.所有樣本采用相同的近鄰數(shù),不能很好地適應(yīng)每個樣本的局部流形結(jié)構(gòu)特征;b.近鄰數(shù)的選擇沒有成熟的算法,常使用的交叉驗證方法存在效率低下等弊端。

    1.1 馬氏距離相似度衡量

    馬氏距離是一種基于數(shù)據(jù)協(xié)方差的距離描述,其定義為

    (2)

    其中:C為數(shù)據(jù)集X的協(xié)方差矩陣。

    馬氏距離對一切非奇異性線性變換具有不變性,不受特征量綱選擇的影響。馬氏距離與歐氏距離的對比如圖1所示,圖中虛線圓為歐氏距離等距線,實線橢圓為馬氏距離等距線。由圖1可以看出,馬氏距離的描述更符合數(shù)據(jù)真實的分布情況,僅當協(xié)方差矩陣C為單位矩陣時,馬氏距離與歐氏距離才會相等[9]。

    圖1 馬氏距離與歐氏距離Fig.1 Mahalanobis distance and Euclidean distance

    基于馬氏距離數(shù)據(jù)分布特征描述方面的優(yōu)勢,提出基于馬氏距離的樣本相似矩陣構(gòu)造方法,相似矩陣元素表達如下

    (3)

    其中:αij=exp(-dij/β);β為所有樣本之間馬氏距離均值的平方。

    1.2 核密度估計調(diào)整

    樣本近鄰構(gòu)造時,空間中密集分布的樣本局部特征往往較為相似。對樣本xi來講,其附近區(qū)域其他樣本出現(xiàn)的概率密度越大,表示與其具有相似局部特征的樣本越多,則xi對應(yīng)的樣本近鄰數(shù)應(yīng)該大一些。據(jù)此提出將高斯核密度估計(Gauss kernel density estimation,簡稱GKDE)[10]用于近鄰圖構(gòu)造,對基于相似系數(shù)均值Mi的自適應(yīng)近鄰結(jié)果進行修正調(diào)整。

    高斯核密度估計是一種非參數(shù)概率密度估計方法,能夠在樣本分布先驗知識未知的情況下,根據(jù)樣本自身信息估計總體分布概率密度。高維空間RD中樣本xi的鄰域概率密度高斯核估計為

    (4)

    (5)

    其中:floor表示數(shù)據(jù)向負無窮大方向取整。

    通過式(5)可知,當xi附近數(shù)據(jù)樣本分布密度大時,該點的近鄰數(shù)自動調(diào)整增大,當其附近數(shù)據(jù)分布稀疏時,近鄰數(shù)自動調(diào)整減小。

    1.3 自適應(yīng)鄰域構(gòu)造流程及特性

    圖2 自適應(yīng)鄰域構(gòu)造流程Fig.2 Flow chart of adaptive neighborhood selection

    自適應(yīng)鄰域構(gòu)造方法最終得到的相似矩陣S一般情況下是非對稱的,主要包含以下幾種情況:

    2) 與第1假設(shè)相同,當樣本xi和xj互相不為近鄰時,相似系數(shù)αij=αji=0,相似矩陣S對稱;

    自適應(yīng)鄰域構(gòu)造方法與傳統(tǒng)k近鄰方法都是為了表達數(shù)據(jù)樣本局部結(jié)構(gòu)的相似性,主要不同為:

    1) k近鄰方法基于歐氏距離計算樣本相似系數(shù),自適應(yīng)鄰域構(gòu)造方法基于馬氏距離計算相似系數(shù),能夠更好地考慮樣本分布特性;

    2) k近鄰方法所有樣本近鄰數(shù)需人為設(shè)置,自適應(yīng)近鄰構(gòu)造方法根據(jù)相似系數(shù)均值Mi確定樣本初始近鄰數(shù),并采用樣本分布概率密度調(diào)整近鄰數(shù);

    3) k近鄰方法所有樣本近鄰數(shù)均相等,自適應(yīng)近鄰構(gòu)造方法每個樣本的近鄰數(shù)不一定相等,且一般都不相等;

    4) k近鄰方法得到的相似矩陣S是對稱的,自適應(yīng)近鄰構(gòu)造方法得到的相似矩陣S一般不對稱。

    1.4 自適應(yīng)鄰域構(gòu)造實例

    采用26個隨機二維數(shù)據(jù)樣本,說明自適應(yīng)鄰域構(gòu)造方法的特點。采用自適應(yīng)鄰域構(gòu)造方法尋找樣本近鄰,以樣本點A(1,4)和點B(8.5,5)的近鄰求解進行說明,并給出k近鄰方法的對比,設(shè)定k=8。表1為自適應(yīng)近鄰求解過程中的參數(shù)及結(jié)果。圖3中4幅圖所示為原始26個樣本點及點A和點B的3種近鄰求解結(jié)果。

    表1 近鄰求解參數(shù)及結(jié)果

    圖3 點A和點B近鄰結(jié)果分布Fig.3 Neighborhood for point A and point B

    2 自適應(yīng)鄰域PCA-LPP流形學習算法

    自適應(yīng)鄰域構(gòu)造方法能夠用于多種基于近鄰圖構(gòu)造的流形學習算法,筆者將其應(yīng)用于一種結(jié)合主元分析與局部保持投影的PCA-LPP改進流形學習算法中,標準局部保持投影算法中相似矩陣S為對稱矩陣,其對稱性對算法的求解帶來很大便利。但自適應(yīng)鄰域構(gòu)造方法得到的相似矩陣S一般為非對稱矩陣,其在PCA-LPP中的求解推導如下。

    PCA-LPP流形學習目的在于將高維空間RD中的數(shù)據(jù)集X=[x1,x2,…,xN],通過投影矩陣W轉(zhuǎn)換投影到低維投影空間Rd(d

    PCA方法尋找數(shù)據(jù)分布方差最大的坐標系,其全局目標函數(shù)為

    (6)

    LPP是一種非線性局部流形保持投影算法,能夠保持投影前后數(shù)據(jù)樣本間的局部結(jié)構(gòu)相似,其局部目標函數(shù)為

    (7)

    其中:L′為拉普拉斯矩陣;Dii和Djj為對角陣,其元素為Dii=∑jsij,Djj=∑isij;D′=Dii+Djj;S′=S+ST;sij為相似系數(shù),構(gòu)成相似矩陣S;W為投影矩陣。

    PCA-LPP流形學習算法兼顧PCA的全局分布方差最大特性及LPP的局部流形保持特性,能夠全面刻畫數(shù)據(jù)的整體和局部特征,其全局目標函數(shù)為

    (8)

    其中:限定條件WTXD′XTW=I可消除隨機尺度因子的影響。

    在LPP局部目標函數(shù)中,自適應(yīng)鄰域構(gòu)造方法得到的相似矩陣S是非對稱矩陣,但局部目標函數(shù)Jl(W)中的矩陣L′是對稱矩陣,其證明如下。

    證明: ? 相似矩陣S

    ∵ Dii=∑jsij,Djj=∑isij均為對角陣

    ∴ Dii和Djj為對稱矩陣

    ∴ D′=Dii+Djj,S′=S+ST為對角對稱矩陣

    ∴ L′=D′-S′為對稱矩陣。

    自適應(yīng)鄰域構(gòu)造PCA-LPP算法的目標函數(shù)式(8)可以通過拉格朗日乘子法[11],將其轉(zhuǎn)化為有約束最大值問題,如式(9)所示

    (9)

    式(9)對W求導并置零可得

    (10)

    化簡式(10)可以得到

    (C-XL′XT)W=λXD′XTW

    (11)

    由式(11)可知,投影矩陣W的求解,實質(zhì)是廣義特征值求解問題,通過求取前d個最大特征值對應(yīng)的特征向量,從而得到投影矩陣W,投影后低維數(shù)據(jù)通過Y=WTX計算。

    由此可以得到自適應(yīng)鄰域PCA-LPP流形學習算法流程,如圖4所示。

    圖4 自適應(yīng)鄰域PCA-LPP算法流程Fig.4 PCA-LPP algorithm based on adaptive neighborhood

    3 齒輪箱故障特征降維分析

    流形學習算法在數(shù)據(jù)可視化以及故障特征降維識別中應(yīng)用廣泛[12-13],下面對自適應(yīng)鄰域PCA-LPP流形學習算法的性能進行驗證。

    采用齒輪箱故障實測聲發(fā)射信號對PCA-LPP的特征降維性能進行驗證。試驗中齒輪箱故障包括正常狀態(tài)、軸承內(nèi)圈故障、外圈故障、內(nèi)外圈復合故障及齒根裂紋故障5種模式,故障軸承為6206,安裝在中間傳動軸,故障齒輪為中間傳動軸大齒輪,一級傳動比為0.5,齒輪箱結(jié)構(gòu)原理及傳感器安裝如圖5所示。試驗中齒輪箱空負載運轉(zhuǎn),其輸入軸轉(zhuǎn)速為1 490 r/min。信號采集設(shè)備為北京軟島DS2-8A型全息聲發(fā)射信號分析儀,設(shè)置采集儀采樣頻率為1 MHz,采集儀硬件濾波參數(shù)為100 kHz~400 kHz帶通濾波,聲發(fā)射傳感器為聲華SR150M型,匹配40 dB前置放大器。每個數(shù)據(jù)樣本長度為1 s,每種故障模式對應(yīng)70個樣本,5種模式共350個樣本。

    對每個故障信號樣本進行db4小波包4層分解,重構(gòu)各子頻帶小波包系數(shù),得到16個子頻帶重構(gòu)分量,各子頻帶寬度為31.25 kHz,計算子頻帶信號能量,得到能量特征向量E=[E1,E2,…,E16],求取16個子頻帶信號的能量熵作為特征向量[8],能量熵的計算方法如下。

    (12)

    采用自適應(yīng)鄰域PCA-LPP方法對所提故障特征向量進行處理,并將降維后特征向量輸入支持向量機進行訓練識別,對比故障模式的分類識別性能,每個故障70個樣本中,40個用于訓練支持向量機分類器,30個用于樣本故障分類識別。

    3.1 自適應(yīng)鄰域PCA-LPP故障特征降維識別分析 齒輪箱故障聲發(fā)射信號的特征向量原始維度為16維,采用自適應(yīng)鄰域PCA-LPP對故障特征向量進行降維處理,降維后特征的維度范圍為3~15維,間隔為1。不同降維數(shù)的特征向量輸入支持向量機的整體故障識別率如圖6所示。由圖6可以看到,原始的16維特征向量經(jīng)降維處理,不同維度的特征向量對應(yīng)的故障辨識率有所不同。辨識率在5維以后出現(xiàn)明顯增高,當降維數(shù)d=8時,整體故障辨識率最高,達到94.67%。其中整體故障辨識率指各單項故障識別率的均值,表示4種故障類型正確識別的總數(shù)占150個測試樣本的比例。

    降維數(shù)d=8時,不同故障的辨識率如表2所示。由表2可知,各類故障模式中,齒根裂紋辨識率最高,因為齒根裂紋與其他幾種故障具有明顯不同,而復合故障辨識率最低,因為復合故障與內(nèi)、外圈故障具有相似特征成分。

    圖6 不同降維數(shù)的故障識別率Fig.6 Fault identification rate of different dimensions

    故障類型齒根裂紋正常狀態(tài)內(nèi)圈故障外圈故障復合故障整體識別率%10096.6793.3393.3390.0094.67

    3.2 不同識別方法性能對比

    對文中的故障特征向量,將k近鄰PCA-LPP方法降維識別結(jié)果、基于歐氏距離的自適應(yīng)PCA-LPP方法降維識別結(jié)果以及原始特征向量不做任何處理的識別結(jié)果與3.1節(jié)的識別結(jié)果進行對比,k近鄰方法的最優(yōu)識別結(jié)果通過交叉驗證獲得。設(shè)置近鄰值k的變換范圍為5~30,間隔為1,降維數(shù)d的范圍為3~15,間隔為1,計算后得到k近鄰方法在k=17,d=7時獲得最高識別率?;跉W氏距離的自適應(yīng)PCA-LPP方法降維識別在降維數(shù)d=8時取得最高識別率。4種方法的識別率對比如表3所示。

    表3 3種情況識別率對比

    從表3可以看到,原始16維特征向量,不經(jīng)任何處理,直接進行支持向量機分類器的訓練和識別,其故障識別率最低,僅為53.33%,而采用流形學習算法降維之后,故障識別率明顯提升。對于同樣的PCA-LPP流形學習算法,不同的鄰域構(gòu)造方法帶來的降維效果差異較大。自適應(yīng)鄰域構(gòu)造方法降維后特征的故障識別率可達94.67%,而k近鄰方法降維后特征的故障識別率為86.67%。不同相似度衡量方法帶來的降維效果也不同,基于馬氏距離的相似度衡量最終故障識別率為94.67%,而基于歐氏距離的相似度衡量最終識別率卻為88.67%。原始特征向量維度高,其中包含故障特征的差異信息,但高維度同時包含過多的冗余信息,因此其故障識別率低。k近鄰方法對所有樣本采用相同的近鄰數(shù),對各樣本局部鄰域的描述不能做到分而劃之,基于歐氏距離的相似度衡量,對數(shù)據(jù)樣本的分布情況不能充分考慮,因此兩種方法特征向量的降維效果不如自適應(yīng)鄰域構(gòu)造方法,故其故障識別率低于自適應(yīng)鄰域PCA-LPP方法。

    4 結(jié)束語

    針對流形學習算法中近鄰圖的構(gòu)造問題,提出一種基于馬氏距離的相似矩陣元素計算。根據(jù)相似矩陣元素均值確定初始鄰域數(shù),并由數(shù)據(jù)樣本鄰域高斯核概率密度估計結(jié)果調(diào)整樣本近鄰數(shù),解決了自適應(yīng)構(gòu)造近鄰圖的問題。將自適應(yīng)鄰域構(gòu)造方法應(yīng)用于兼顧PCA與LPP二者特性的改進流形學習算法,并給出相應(yīng)的理論計算模型,證明了求解拉普拉斯矩陣的對稱性。通過齒輪箱故障聲發(fā)射信號16維能量熵特征的降維識別對自適應(yīng)鄰域PCA-LPP流形學習算法的性能進行驗證,結(jié)果表明故障特征降到8維時識別率最高,可以達到94.67%。對比研究表明,自適應(yīng)鄰域構(gòu)造PCA-LPP方法的降維識別性能優(yōu)于k近鄰PCA-LPP方法以及原始特征向量的識別率。

    [1] 趙振華,郝曉弘. 局部保持鑒別投影及其在人臉識別中的應(yīng)用[J]. 電子與信息學報, 2013,35(2):463-466.

    Zhao Zhenhua, Hao Xiaohong. Linear locality preserving and discriminating projection for face recognition[J]. Journal of Electronics & Information Technology, 2013,35(2): 463-466.(in Chinese)

    [2] He Xiaofei, Niyogi P. Locality preserving projections[C]∥Neural Information Processing Systems16. Cambridge,USA: MIT Press, 2004:153-160.

    [3] 許慶誠,胡建中. 基于改進增量LE的壓縮機故障特征提取方法[J]. 儀器儀表學報,2013,34(4):791-796.

    Xu Qingcheng, Hu Jianzhong. Fault feature extraction method for compressor based on improved incremental Laplacian eigenmap algorithm[J]. Chinese Journal of Scientific Instrument, 2013,34(4):791-796.(in Chinese)

    [4] 何強,蔡洪,韓壯志,等. 基于非線性流形學習的ISAR目標識別研究[J]. 電子學報, 2010,38(3):585-590.

    He Qiang, Cai Hong, Han Zhuangzhi, et al. ISAR target recognition based on non-linear manifold learning[J]. Acta Electronica Sinica,2010,38(3):585-590.(in Chinese)

    [5] Yang Bo, Chen Songcan. Sample-dependent graph construction with application to dimensionality reduction[J]. Neurocomputing, 2010,74(5): 301-314.

    [6] 劉鳳連,汪日偉,程俊,等. 圖像特征提取中領(lǐng)域尺寸和本征維數(shù)的自動選擇算法[J]. 光電子·激光, 2013,24(12): 2416-2420.

    Liu Fenglian, Wang Riwei, Cheng Jun, et al. Auto-selected algorithm of the neighborhood size and intrinsic dimension for image feature extraction[J]. Journal fo Optoelectronics·Laser, 2013,24(12): 2416- 2420.(in Chinese)

    [7] 黃璞,唐振民. 無參數(shù)局部保持投影及人臉識別[J]. 模式識別與人工智能, 2013,26(9): 865-871.

    Huang Pu, Tang Zhenmin. Parameter-free locality preserving projections and face recognition[J]. Pattern Recognition and Artificial Intelligence, 2013,26(9): 865-871.(in Chinese)

    [8] 李城梁,王仲生,姜洪開,等. 自適應(yīng)Hessian LLE在機械故障特征提取中的應(yīng)用[J]. 振動工程學報, 2013,26(5): 758-763.

    Li Chengliang, Wang Zhongsheng, Jiang Hongkai, et al. Adaptive Hessian LLE in mechanical fault feature extraction[J]. Journal of Vibration Engineering,2013, 26(5): 758-763.(in Chinese)

    [9] 駱志高,李旭東,趙俊麗,等. 利用馬氏距離判別法準確實現(xiàn)對裂紋的識別[J]. 振動與沖擊,2013,32(21): 186-188.

    Luo Zhigao, Li Xudong, Zhao Junli, et al. Crack identification with mahalanobis distance discrimination method[J]. Journal of Vibration and Shock,2013,32(21): 186-188.(in Chinese)

    [10]孫即祥. 現(xiàn)代模式識別[M]. 長沙:國防科技大學出版社,2002:136-140.

    [11]李鋒,王家序,楊榮松. 有監(jiān)督不相關(guān)正交局部保持映射故障辨識[J]. 儀器儀表學報, 2013,34(5): 1113-1116.

    Li Feng, Wang Jiaxu, Yang Rongsong. Fault identification method on supervised uncorrelated orthogonal locality preserving projection[J]. Chinese Journal of Scientific Instrument,2013,34(5): 1113-1116.(in Chinese)

    [12]宋濤,湯寶平,李鋒. 基于流形學習和K-最近鄰分類器的旋轉(zhuǎn)機械故障診斷方法[J]. 振動與沖擊,2013,32(5): 149-153.

    Song Tao, Tang Baoping, Li Feng. Fault diagnosis method for rotating machinery based on manifold learning and K-nearest neighbor classifier[J]. Journal of Vibration and Shock, 2013,32(5):149-153.(in Chinese)

    [13] 劉忠寶,潘廣貞,趙文娟. 流形判別分析[J]. 電子與信息學報,2013,35(9):2047-2050.

    Liu Zhongbao, Pan Guangzhen, Zhao Wenjuan. Manifold-based discriminant analysis[J]. Journal of Electronics & Information Technology,2013,35(9): 2047-2050.(in Chinese)

    10.16450/j.cnki.issn.1004-6801.2016.06.027

    *國家自然科學基金資助項目(50775219)

    2014-11-12;

    2015-03-27

    TH165; TN911

    張曉濤,男,1987年5月生,博士生。主要研究方向為機械系統(tǒng)性能檢測與故障診斷。曾發(fā)表《基于SVD與Fast Kurtogram算法的滾動軸承故障診斷》(《振動與沖擊》2014年第33卷第10期)等論文。 E-mail:headic@163.com

    猜你喜歡
    構(gòu)造方法流形降維
    DC-DC變換器分層級構(gòu)造方法
    混動成為降維打擊的實力 東風風神皓極
    車主之友(2022年4期)2022-08-27 00:57:12
    緊流形上的Schr?dinger算子的譜間隙估計
    降維打擊
    海峽姐妹(2019年12期)2020-01-14 03:24:40
    迷向表示分為6個不可約直和的旗流形上不變愛因斯坦度量
    Nearly Kaehler流形S3×S3上的切觸拉格朗日子流形
    《夢溪筆談》“甲子納音”構(gòu)造方法的數(shù)學分析
    幾乎最佳屏蔽二進序列偶構(gòu)造方法
    基于多故障流形的旋轉(zhuǎn)機械故障診斷
    拋物化Navier-Stokes方程的降維仿真模型
    計算物理(2014年1期)2014-03-11 17:00:18
    欧美日韩视频高清一区二区三区二| 久久99精品国语久久久| 免费av观看视频| 亚洲经典国产精华液单| 联通29元200g的流量卡| 又粗又硬又长又爽又黄的视频| 99热这里只有是精品50| 免费高清在线观看视频在线观看| 国产午夜精品一二区理论片| 成人二区视频| 日本色播在线视频| 成人av在线播放网站| 亚洲精品成人久久久久久| 成人一区二区视频在线观看| 欧美另类一区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲av二区三区四区| 亚洲精品第二区| 国产白丝娇喘喷水9色精品| 麻豆精品久久久久久蜜桃| 有码 亚洲区| 亚洲欧美一区二区三区黑人 | 国产精品一区二区三区四区免费观看| 国产人妻一区二区三区在| 日韩强制内射视频| 国产有黄有色有爽视频| 久久午夜福利片| 永久免费av网站大全| 别揉我奶头 嗯啊视频| 亚洲18禁久久av| 国产成人freesex在线| 精品人妻熟女av久视频| 国内揄拍国产精品人妻在线| 国产在线男女| 国产av国产精品国产| 日韩欧美 国产精品| 又爽又黄a免费视频| 国产成人精品婷婷| 亚洲丝袜综合中文字幕| 男女视频在线观看网站免费| 九九在线视频观看精品| 亚洲欧美日韩卡通动漫| 一区二区三区高清视频在线| 国产亚洲一区二区精品| 亚洲国产欧美在线一区| 欧美成人a在线观看| 午夜福利高清视频| 国产高清三级在线| 国产精品女同一区二区软件| 国产爱豆传媒在线观看| 久久久a久久爽久久v久久| 欧美一区二区亚洲| 人妻夜夜爽99麻豆av| 亚洲va在线va天堂va国产| 日日啪夜夜撸| 精品国内亚洲2022精品成人| 99久国产av精品国产电影| 一级二级三级毛片免费看| 亚洲电影在线观看av| 午夜久久久久精精品| 久久韩国三级中文字幕| 国产真实伦视频高清在线观看| 美女被艹到高潮喷水动态| 性色avwww在线观看| 久久精品夜夜夜夜夜久久蜜豆| 欧美97在线视频| 男插女下体视频免费在线播放| 国产黄a三级三级三级人| 大话2 男鬼变身卡| 国产麻豆成人av免费视频| 精品久久久久久久久亚洲| 伦理电影大哥的女人| 黄色欧美视频在线观看| 一区二区三区乱码不卡18| 99热网站在线观看| 男女啪啪激烈高潮av片| 国产 一区 欧美 日韩| 男女边吃奶边做爰视频| 男女啪啪激烈高潮av片| 欧美3d第一页| 黑人高潮一二区| www.av在线官网国产| 精品99又大又爽又粗少妇毛片| 哪个播放器可以免费观看大片| 欧美xxxx性猛交bbbb| 老司机影院毛片| 联通29元200g的流量卡| 免费在线观看成人毛片| 99九九线精品视频在线观看视频| 国产高潮美女av| 成人综合一区亚洲| 亚洲精品第二区| 少妇的逼好多水| 在线观看人妻少妇| 69人妻影院| 免费高清在线观看视频在线观看| 夜夜爽夜夜爽视频| 欧美zozozo另类| 国产黄色免费在线视频| 欧美高清成人免费视频www| 丰满少妇做爰视频| 深夜a级毛片| 亚洲丝袜综合中文字幕| 观看免费一级毛片| 五月天丁香电影| 成人国产麻豆网| 黄色一级大片看看| 超碰97精品在线观看| 一边亲一边摸免费视频| 看十八女毛片水多多多| 久久久久性生活片| 亚洲精品国产成人久久av| 男人舔奶头视频| 日韩电影二区| 亚洲成人久久爱视频| 亚洲精品国产成人久久av| 免费看a级黄色片| 欧美激情国产日韩精品一区| 超碰97精品在线观看| 欧美成人a在线观看| 国产精品久久久久久精品电影| 在线天堂最新版资源| 成人美女网站在线观看视频| 乱系列少妇在线播放| 亚洲av.av天堂| 亚洲av日韩在线播放| 国产成人一区二区在线| 婷婷色av中文字幕| 超碰av人人做人人爽久久| 日韩伦理黄色片| 欧美最新免费一区二区三区| 免费黄网站久久成人精品| 欧美zozozo另类| 国产成人免费观看mmmm| a级一级毛片免费在线观看| 观看美女的网站| av国产免费在线观看| 久久久久久久久久成人| 晚上一个人看的免费电影| 国产精品久久视频播放| 久久99热这里只有精品18| 永久免费av网站大全| 亚洲国产日韩欧美精品在线观看| 精品人妻一区二区三区麻豆| 超碰97精品在线观看| 丝袜美腿在线中文| 99久久九九国产精品国产免费| 丝瓜视频免费看黄片| 亚洲欧美清纯卡通| 久久久精品94久久精品| 精品久久久久久久久亚洲| 久久久精品免费免费高清| 尤物成人国产欧美一区二区三区| 大又大粗又爽又黄少妇毛片口| 亚洲久久久久久中文字幕| 在线观看美女被高潮喷水网站| 中国国产av一级| 国产熟女欧美一区二区| 99热这里只有是精品50| 久久久久网色| 丝袜美腿在线中文| 我的女老师完整版在线观看| 嫩草影院新地址| 淫秽高清视频在线观看| 免费av观看视频| 少妇熟女欧美另类| 亚洲精品色激情综合| 人人妻人人澡人人爽人人夜夜 | 日韩人妻高清精品专区| 国产毛片a区久久久久| 成人综合一区亚洲| 国产精品蜜桃在线观看| 久久精品久久精品一区二区三区| 国产激情偷乱视频一区二区| 国产精品美女特级片免费视频播放器| 日日撸夜夜添| 亚洲欧洲国产日韩| 男女啪啪激烈高潮av片| 欧美精品一区二区大全| 久久久久久久久久久免费av| 久久久久免费精品人妻一区二区| 国产日韩欧美在线精品| 国产综合精华液| 久久这里只有精品中国| 人人妻人人澡欧美一区二区| 少妇高潮的动态图| 免费在线观看成人毛片| 亚洲伊人久久精品综合| 中文欧美无线码| 日韩,欧美,国产一区二区三区| 日韩av在线大香蕉| 日本wwww免费看| 五月伊人婷婷丁香| 国产免费福利视频在线观看| 色综合色国产| 欧美潮喷喷水| 婷婷色av中文字幕| 午夜福利视频1000在线观看| 国国产精品蜜臀av免费| 国产av在哪里看| 黄色日韩在线| 日韩电影二区| 欧美+日韩+精品| 九草在线视频观看| 伦理电影大哥的女人| 国产精品久久久久久久电影| 高清午夜精品一区二区三区| 天天一区二区日本电影三级| 成年av动漫网址| 久久99热这里只频精品6学生| 亚洲一级一片aⅴ在线观看| ponron亚洲| 成人亚洲欧美一区二区av| 老女人水多毛片| 69av精品久久久久久| 三级国产精品欧美在线观看| 国产精品不卡视频一区二区| 亚洲在线自拍视频| 晚上一个人看的免费电影| 97在线视频观看| 国产黄频视频在线观看| 国产伦理片在线播放av一区| 免费观看精品视频网站| 国产成人精品一,二区| 国产美女午夜福利| 久久久精品欧美日韩精品| 欧美日韩精品成人综合77777| 亚洲av成人精品一二三区| 18禁在线播放成人免费| 成人特级av手机在线观看| av一本久久久久| 又粗又硬又长又爽又黄的视频| 国产av码专区亚洲av| 插逼视频在线观看| 中文资源天堂在线| 高清视频免费观看一区二区 | 亚洲精品日韩av片在线观看| 国产精品久久久久久精品电影| 色吧在线观看| 男人舔女人下体高潮全视频| 午夜福利网站1000一区二区三区| 亚洲性久久影院| 久久久久久久久久成人| 99re6热这里在线精品视频| 80岁老熟妇乱子伦牲交| 国产综合懂色| 特大巨黑吊av在线直播| 欧美丝袜亚洲另类| 水蜜桃什么品种好| 在线观看免费高清a一片| 日韩欧美精品免费久久| 免费观看性生交大片5| 国产 亚洲一区二区三区 | 熟妇人妻久久中文字幕3abv| 国产伦在线观看视频一区| 熟妇人妻久久中文字幕3abv| 亚洲经典国产精华液单| 亚洲欧美一区二区三区黑人 | 在线播放无遮挡| 亚洲欧美成人精品一区二区| 日日摸夜夜添夜夜添av毛片| 黄色一级大片看看| 亚洲经典国产精华液单| 色综合色国产| 国产精品一区二区三区四区免费观看| 亚洲成人av在线免费| 国产午夜精品论理片| 欧美潮喷喷水| 可以在线观看毛片的网站| 男女那种视频在线观看| 国产一级毛片在线| 一级av片app| 欧美极品一区二区三区四区| 中文字幕免费在线视频6| 人妻少妇偷人精品九色| 欧美xxxx性猛交bbbb| 国产人妻一区二区三区在| 简卡轻食公司| 久久久久久久大尺度免费视频| 国产亚洲一区二区精品| 欧美日韩亚洲高清精品| 狠狠精品人妻久久久久久综合| 国产精品伦人一区二区| 中文字幕制服av| 国产精品久久久久久久久免| 国产美女午夜福利| 欧美3d第一页| 久久精品夜色国产| 中文字幕免费在线视频6| 精品99又大又爽又粗少妇毛片| 国产久久久一区二区三区| 亚洲不卡免费看| 内地一区二区视频在线| 十八禁国产超污无遮挡网站| 久久久久久久久久久丰满| 精品久久久久久久久久久久久| 别揉我奶头 嗯啊视频| 国产大屁股一区二区在线视频| 91久久精品国产一区二区三区| 高清日韩中文字幕在线| 一个人看的www免费观看视频| 亚洲国产高清在线一区二区三| 网址你懂的国产日韩在线| 国产午夜精品论理片| 日日摸夜夜添夜夜添av毛片| 99久国产av精品| 夜夜爽夜夜爽视频| 丝瓜视频免费看黄片| 联通29元200g的流量卡| 搞女人的毛片| 欧美区成人在线视频| 亚洲真实伦在线观看| xxx大片免费视频| 我要看日韩黄色一级片| 97精品久久久久久久久久精品| 亚洲真实伦在线观看| 亚洲欧美精品专区久久| 国产国拍精品亚洲av在线观看| 亚洲欧美日韩卡通动漫| 午夜爱爱视频在线播放| 99热6这里只有精品| 99热这里只有精品一区| 亚洲丝袜综合中文字幕| 成年女人看的毛片在线观看| 最新中文字幕久久久久| 99久久中文字幕三级久久日本| 一级毛片电影观看| eeuss影院久久| 天堂中文最新版在线下载 | 日韩不卡一区二区三区视频在线| 99热网站在线观看| 一级毛片电影观看| 欧美日韩国产mv在线观看视频 | 国产在线一区二区三区精| 精品国内亚洲2022精品成人| 国内精品一区二区在线观看| 亚洲欧美一区二区三区黑人 | 春色校园在线视频观看| 欧美潮喷喷水| 春色校园在线视频观看| 国产精品国产三级国产专区5o| 久久久久久久大尺度免费视频| 少妇熟女欧美另类| 午夜激情欧美在线| 国产精品国产三级国产专区5o| 国产探花在线观看一区二区| 麻豆久久精品国产亚洲av| 国产成人精品婷婷| 在线观看免费高清a一片| 婷婷色av中文字幕| 午夜免费观看性视频| 全区人妻精品视频| 国内揄拍国产精品人妻在线| 久久久精品94久久精品| 国产淫片久久久久久久久| 午夜精品一区二区三区免费看| 国产黄片美女视频| 午夜日本视频在线| 男女边吃奶边做爰视频| 日韩亚洲欧美综合| 赤兔流量卡办理| 91久久精品电影网| 亚洲欧洲日产国产| 狂野欧美白嫩少妇大欣赏| 免费av观看视频| 大陆偷拍与自拍| 亚洲av男天堂| 亚洲精品国产av成人精品| 老师上课跳d突然被开到最大视频| 精品国产一区二区三区久久久樱花 | 久久精品国产亚洲网站| 赤兔流量卡办理| 汤姆久久久久久久影院中文字幕 | 春色校园在线视频观看| 又黄又爽又刺激的免费视频.| 亚洲精品自拍成人| av在线播放精品| 亚洲人成网站高清观看| 22中文网久久字幕| 天堂√8在线中文| 国产精品一区www在线观看| 久久这里有精品视频免费| 久久久亚洲精品成人影院| 欧美激情在线99| 午夜福利成人在线免费观看| 欧美人与善性xxx| av播播在线观看一区| 亚洲av一区综合| 国产精品伦人一区二区| 国产午夜精品论理片| 只有这里有精品99| 免费大片18禁| 国产精品国产三级国产专区5o| 久久久久久久久中文| 亚洲av中文av极速乱| 国产亚洲最大av| 99热全是精品| 欧美激情国产日韩精品一区| 街头女战士在线观看网站| 美女xxoo啪啪120秒动态图| 亚洲av免费高清在线观看| 午夜精品一区二区三区免费看| 精品国产露脸久久av麻豆 | 激情五月婷婷亚洲| 韩国av在线不卡| 久久久久久久久久成人| 美女大奶头视频| 69人妻影院| 欧美极品一区二区三区四区| www.av在线官网国产| 九九久久精品国产亚洲av麻豆| videossex国产| 国产高清三级在线| 国产精品一区二区性色av| 九色成人免费人妻av| 狂野欧美白嫩少妇大欣赏| 免费少妇av软件| 国产免费福利视频在线观看| 日本猛色少妇xxxxx猛交久久| av女优亚洲男人天堂| 亚洲人成网站在线观看播放| 只有这里有精品99| 偷拍熟女少妇极品色| 国产麻豆成人av免费视频| 国产亚洲精品av在线| 日韩欧美精品v在线| 日本免费a在线| 久热久热在线精品观看| 一个人看的www免费观看视频| 最近最新中文字幕免费大全7| 亚洲精品国产av成人精品| 国产日韩欧美在线精品| 2022亚洲国产成人精品| 久久久久久久久久人人人人人人| 99视频精品全部免费 在线| 国产亚洲午夜精品一区二区久久 | 一区二区三区四区激情视频| 日韩视频在线欧美| 免费高清在线观看视频在线观看| 一个人观看的视频www高清免费观看| 特大巨黑吊av在线直播| 老司机影院毛片| 久久久国产一区二区| 一级黄片播放器| 久久久久免费精品人妻一区二区| 精品人妻视频免费看| 亚洲aⅴ乱码一区二区在线播放| 日本一本二区三区精品| 两个人的视频大全免费| 国产黄色小视频在线观看| 天堂影院成人在线观看| 色网站视频免费| 国产黄a三级三级三级人| 精品99又大又爽又粗少妇毛片| 国产精品一区二区性色av| 搡老妇女老女人老熟妇| 成年av动漫网址| 国产成人精品福利久久| 美女黄网站色视频| 99久国产av精品国产电影| 亚洲无线观看免费| 亚洲精品色激情综合| 日本三级黄在线观看| 久久99蜜桃精品久久| 丰满少妇做爰视频| 国产探花在线观看一区二区| 91aial.com中文字幕在线观看| 大香蕉97超碰在线| 亚洲性久久影院| av又黄又爽大尺度在线免费看| 三级国产精品片| 国产白丝娇喘喷水9色精品| 久久人人爽人人爽人人片va| 麻豆精品久久久久久蜜桃| 国产老妇女一区| a级毛片免费高清观看在线播放| 免费黄网站久久成人精品| 日韩av在线免费看完整版不卡| 大片免费播放器 马上看| 国产午夜精品一二区理论片| 亚洲欧美成人综合另类久久久| 色视频www国产| 日韩人妻高清精品专区| 在线 av 中文字幕| 大陆偷拍与自拍| 亚洲aⅴ乱码一区二区在线播放| 国产成人午夜福利电影在线观看| 99久久精品热视频| 久99久视频精品免费| 精品不卡国产一区二区三区| 国产高清国产精品国产三级 | 国产精品av视频在线免费观看| 不卡视频在线观看欧美| 欧美不卡视频在线免费观看| 一级片'在线观看视频| 晚上一个人看的免费电影| 国内精品美女久久久久久| 欧美性感艳星| 国产亚洲av片在线观看秒播厂 | 97人妻精品一区二区三区麻豆| 日韩亚洲欧美综合| 看免费成人av毛片| 夫妻午夜视频| 最近最新中文字幕大全电影3| 亚洲精品久久久久久婷婷小说| 丝袜美腿在线中文| 亚洲人成网站在线观看播放| 免费av毛片视频| 建设人人有责人人尽责人人享有的 | 国产成人a区在线观看| 亚洲国产精品sss在线观看| 成年人午夜在线观看视频 | ponron亚洲| 亚洲婷婷狠狠爱综合网| 人人妻人人看人人澡| 亚洲精品乱久久久久久| 麻豆av噜噜一区二区三区| 亚洲第一区二区三区不卡| 久久99热这里只频精品6学生| 亚洲天堂国产精品一区在线| 亚洲av国产av综合av卡| 亚洲精品国产av成人精品| 人妻一区二区av| 五月玫瑰六月丁香| 狠狠精品人妻久久久久久综合| 国内少妇人妻偷人精品xxx网站| 69av精品久久久久久| 日韩欧美一区视频在线观看 | 午夜福利在线在线| 最后的刺客免费高清国语| 免费看光身美女| 成人漫画全彩无遮挡| 在线观看人妻少妇| 日日啪夜夜撸| 精品久久久精品久久久| 欧美激情国产日韩精品一区| 床上黄色一级片| 极品少妇高潮喷水抽搐| 日本色播在线视频| 三级毛片av免费| 免费观看在线日韩| 亚洲aⅴ乱码一区二区在线播放| 可以在线观看毛片的网站| 久久久久久久久大av| 国产v大片淫在线免费观看| 天天躁夜夜躁狠狠久久av| 亚洲综合精品二区| 黄色配什么色好看| 激情五月婷婷亚洲| 国产成人a区在线观看| 草草在线视频免费看| 国产视频内射| 欧美精品一区二区大全| 一级毛片我不卡| 精品少妇黑人巨大在线播放| 久久人人爽人人爽人人片va| 91aial.com中文字幕在线观看| 国产人妻一区二区三区在| 黄色配什么色好看| 久久久精品94久久精品| 一级爰片在线观看| 久久久久精品久久久久真实原创| 99热这里只有是精品在线观看| 欧美精品一区二区大全| av卡一久久| 岛国毛片在线播放| 在线观看美女被高潮喷水网站| 大又大粗又爽又黄少妇毛片口| 国产精品久久久久久精品电影| 欧美3d第一页| 成人毛片a级毛片在线播放| 国产精品久久视频播放| 欧美人与善性xxx| 舔av片在线| 波多野结衣巨乳人妻| 久久人人爽人人爽人人片va| 亚洲av在线观看美女高潮| 人人妻人人澡欧美一区二区| 亚洲人与动物交配视频| av在线蜜桃| 亚洲最大成人手机在线| .国产精品久久| 丝袜美腿在线中文| 少妇高潮的动态图| 啦啦啦中文免费视频观看日本| 日韩在线高清观看一区二区三区| 一级毛片 在线播放| 亚洲精品中文字幕在线视频 | 精品99又大又爽又粗少妇毛片| 国产 一区精品| 成人二区视频| av播播在线观看一区| 免费在线观看成人毛片| 亚洲精品aⅴ在线观看| 亚洲精品国产av成人精品| 大陆偷拍与自拍| 好男人视频免费观看在线| 久久久久久久久久久免费av| 欧美一区二区亚洲| 国产欧美日韩精品一区二区| 国产一区二区亚洲精品在线观看| a级一级毛片免费在线观看| 亚洲av成人av| 欧美97在线视频| 亚州av有码| av专区在线播放| av一本久久久久| 免费av观看视频| 亚洲av一区综合| 夜夜爽夜夜爽视频| 成人国产麻豆网| 别揉我奶头 嗯啊视频| 国产av在哪里看| 成人无遮挡网站| 国产男女超爽视频在线观看|