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

    基于變點分析思想的高維機械噪聲數(shù)據(jù)本征維估計方法

    2016-05-04 18:40:44梁勝杰張志華高紹忠胡俊波李大偉
    船舶力學(xué) 2016年11期
    關(guān)鍵詞:變點本征高維

    梁勝杰,張志華,高紹忠,胡俊波,李大偉

    (1.91697部隊,山東 青島 266405;2.海軍工程大學(xué) 科研部,武漢 430033;3.91550部隊,遼寧 大連116023)

    基于變點分析思想的高維機械噪聲數(shù)據(jù)本征維估計方法

    梁勝杰1,張志華2,高紹忠1,胡俊波2,李大偉3

    (1.91697部隊,山東 青島 266405;2.海軍工程大學(xué) 科研部,武漢 430033;3.91550部隊,遼寧 大連116023)

    利用變點分析思想,針對高維數(shù)據(jù)協(xié)方差矩陣的特征值曲線特點,提出一種本征維估計方法:基于均值估計的特征值變點法,并應(yīng)用于某雙層圓柱殼體的機械噪聲數(shù)據(jù)處理中。通過與其他現(xiàn)有常用方法的比較分析,驗證了基于均值估計特征值變點法的可行性。結(jié)果顯示:基于均值估計的特征值變點法理論簡單、計算簡便,計算結(jié)果客觀、合理,對某雙層圓柱殼體機械噪聲數(shù)據(jù)的處理效果較好。

    高維數(shù)據(jù);本征維;變點分析;判別分析;機械噪聲數(shù)據(jù)

    0 引 言

    智能計算技術(shù)的發(fā)展使得工程上待處理的數(shù)據(jù)對象變得越來越復(fù)雜,數(shù)據(jù)維數(shù)越來越高,如遙感圖像數(shù)據(jù)、機械振動信號、雷達目標信號、地震資料數(shù)據(jù)等等。高維數(shù)據(jù)的特征空間中通常含有許多冗余特征,這些特征增加了學(xué)習(xí)、訓(xùn)練的時間及空間復(fù)雜度,降低了分析處理的精度。因此,在對高維數(shù)據(jù)進行分析前,應(yīng)進行降維預(yù)處理[1]。

    降維是指從原始數(shù)據(jù)較多的特征量中利用數(shù)學(xué)統(tǒng)計方法提取出具有代表意義的典型特征變量,以達到簡化數(shù)據(jù)的目的。常用方法如主成分分析法、因子分析法等。為敘述方便,我們將提取出來的典型特征量稱為本征特征量(如主成分分析法中的主元、因子分析法中的因子等),典型特征量的個數(shù)稱為本征維。降維過程中,一個重要步驟就是確定本征維,如主成分分析方法中的主元數(shù)目[2]。錯誤的本征維(低估或高估)將導(dǎo)致數(shù)據(jù)分析結(jié)果的錯誤,如低估本征維將導(dǎo)致有用信息的損失及造成對計算結(jié)果本質(zhì)的曲解;高估本征維將保留具有較少重要載荷的特征量,這些特征量很難被解釋其意義[3]。

    目前本征維估計方法主要分為兩類:基于特征值的估計方法和基于分形思想的估計方法?;诜中嗡枷氲墓烙嫹椒梢暂^準確地估計出數(shù)據(jù)的本征維[4],但為了獲得D維數(shù)據(jù)準確的本征維,要求觀測樣本量滿足的條件,樣本量過大,大大增加了試驗難度;基于特征值的估計方法原理簡單、計算方便,在工程上得到廣泛應(yīng)用,常用方法如以“特征值大于1”為標準的方法(Eigenvalue-greaterthan-one rule,K1法)[5]、平行分析法(Parallel Analysis,PA法)[6],碎石圖檢驗法(Scree Test,SCREE法)[7],最小平均偏相關(guān)法(Minimum Average Partial,MAP法)[8-9]和累計貢獻率法(Proportion of total variation,PTV法)[10]等。因此,本文將基于特征值研究高維數(shù)據(jù)的本征維估計方法。

    基于特征值的估計方法均依據(jù)原始數(shù)據(jù)協(xié)方差矩陣的特征值進行計算,其原理不同,估計結(jié)果也不盡相同,各有特點:K1法過于主觀,人們無法解釋為什么認為特征值1.01對應(yīng)的主元是重要的,而認為特征值0.99對應(yīng)的主元不重要[11];PA法是一種蒙特卡洛模擬法,計算中將花費大量時間,不適宜處理大樣本和高維數(shù)據(jù)集;SCREE法比K1法更加主觀和模棱兩可,尤其是在“肘點”不明顯或多于兩個“肘點”情況下[12];MAP法相比上述幾種方法來說能夠得到較精確的結(jié)果,但要求主元中含有較多的原始變量[3],不適宜處理維數(shù)較低的數(shù)據(jù)集;PTV法在確定累計貢獻率標準時同樣具有較強的主觀性,如可選50%[13]、80%[14]或95%[15]等。由于上述方法的種種缺點,導(dǎo)致目前比較普遍使用的統(tǒng)計程序也未能為使用者提供最精確的本征維估計方法[16]。

    鑒于上述問題,本文將針對特征值曲線簡單、直觀的優(yōu)點,利用統(tǒng)計學(xué)中成熟的變點分析思想,提出一種明確、客觀的高維數(shù)據(jù)本征維估計方法:基于均值估計的特征值變點法(Eigenvalue Chang-point Method based mean estimating,meanECP),并應(yīng)用于某雙層圓柱殼體機械噪聲數(shù)據(jù)的預(yù)處理之中。與現(xiàn)有方法相比,本文方法效果較好。

    1 特征值變點法

    變點分析研究始于20世紀50年代。自誕生以來,變點問題一直是統(tǒng)計中的一個熱門話題。變點技術(shù)早期主要應(yīng)用于工業(yè)自動控制領(lǐng)域,隨著計算機技術(shù)的發(fā)展,現(xiàn)在已在醫(yī)學(xué)、金融、水文研究、交通流研究等[17-21]領(lǐng)域得到廣泛應(yīng)用。

    1.1 特征值變點定義

    變點是指“模型中的某個或某些量起突然變化之點”[22]。這種突然變化往往反映事物的某種質(zhì)的變化。對工程中的高維數(shù)據(jù)來說,在利用降維方法降維過程中,首先需要計算其協(xié)方差矩陣的特征值,并將這些特征值按從大至小的次序排列。則在高維數(shù)據(jù)協(xié)方差矩陣的特征值曲線中某個特征值起突然變化之點就是變點,就稱為特征值變點。

    利用特征值變點可以將高維數(shù)據(jù)協(xié)方差矩陣特征值對應(yīng)的特征量按典型和非典型分開,得到的典型特征量即認為是原始數(shù)據(jù)經(jīng)降維后的本征維,典型特征量的個數(shù)即原始數(shù)據(jù)的本征維。特征值變點分析的主要任務(wù)是對這種突變點給予客觀的劃分標準。事實上,SCREE法利用觀測所確定的“肘點”也是特征值曲線的一個變點,而此變點是通過觀察碎石圖得到,受主觀影響較大。本文將利用變點分析的統(tǒng)計知識給出特征值曲線變點的客觀數(shù)值解。

    1.2 數(shù)學(xué)模型

    變點分析包括兩類模型:離散型和連續(xù)型。已知高維數(shù)據(jù)協(xié)方差矩陣的特征值是按從大到小排列的離散值,本文將建立離散型特征值變點分析模型。設(shè)λ1,λ2,…,λn為來自母體的特征值抽樣,且

    λ1,λ2,…,λn之間相互獨立,若

    現(xiàn)在成熟的變點分析方法有:最小二乘法、極大似然法、Bayes法和局部比較法等[23]。本文將在上述方法的基礎(chǔ)上,利用變點分析思想提出一種用于估計高維數(shù)據(jù)本征維的特征值變點法。

    1.3 均值估計法

    對于相互獨立的λ1,λ2,…,λn特征值抽樣,設(shè)

    其中:a1≠a2,隨機誤差ξ1,ξ2,…,ξn相互獨立,期望為0,方差為σ2,則稱k為特征值變點。

    在特征值曲線上取第i個特征值點,位置坐標為(i,λi),記為b點。其前后兩個特征值點c點:(i+ 1,λi+1)和a點記線段與夾角為B。

    由 λi-1≥λi≥λi+1可知如圖1中(Ⅰ)、(Ⅱ)兩種情況所示。

    圖1 相鄰三個特征值點夾角示意圖Fig.1 Sketch map of included angle between the three contiguous eigenvalue points

    當(dāng)B=π時,點a、b、c在同一直線上,此時無變點,三點對應(yīng)的特征值均值相同。如圖1(Ⅲ)所示;隨著夾角B逐漸趨向于當(dāng)λi-1=λi≠λi+1或λi-1≠λi=λi+1時,B達到最小值,此時變點最明顯,變點位置為i。由上述分析可知:當(dāng)無變點時,有即點a、c的縱坐標均值估計與點b的橫坐標相等;當(dāng)有變點時即點a、c的縱坐標均值估計與點b的橫坐標存在一定誤差,且當(dāng)時誤差最大,此時變點最明顯。因此,本文將利用均值估計誤差Δλi(1≤i≤n)對變點位置進行估計:

    除變點估計技術(shù)外,變點分析方法還應(yīng)包括變點檢驗技術(shù)[22]。由于本文處理的對象是特征值曲線,且主要是通過估計一個較小的本征維達到降低觀測數(shù)據(jù)維數(shù)的目的,故認為該特征值曲線存在變點,而不需要做變點檢驗。

    1.4 估計步驟

    設(shè)高維數(shù)據(jù)Y是一n×N型矩陣,即

    其中:n為維數(shù),N為樣本個數(shù)。為簡化敘述,通常假定Y是經(jīng)過中心化后樣本矩陣。則meanECP法估計高維數(shù)據(jù)本征維的步驟如下:

    (1)標準化高維數(shù)據(jù)矩陣:

    其中:Q是Y的協(xié)方差矩陣平方根的逆矩陣;

    (3)計算

    其中:i=2,…,n-1;

    2 方法應(yīng)用及比較分析

    將meanECP法用于某雙層圓柱殼體機械噪聲數(shù)據(jù)的本征維估計中,并與其它常用方法進行比較。

    2.1 數(shù)據(jù)來源

    某雙層圓柱殼體的機械振動噪聲數(shù)據(jù)主要來自三個設(shè)備:海水泵、電機和激振器。在殼體的不同位置布置19個加速度傳感器(每個傳感器稱為一通道),設(shè)定八個不同工況,如表1所示[24]。

    設(shè)定采樣時間8 s,采樣頻率2 048 Hz,則每一工況在每一通道上都可測得16 384個樣本數(shù)據(jù)。

    由于同一個通道下的采樣點較多,且特征單一,不易分析數(shù)據(jù)的規(guī)律特點,因此需要對樣本點進行預(yù)處理:

    (1)選定信號處理中常用的10個特征量[24]:波形因數(shù),峰值因數(shù),脈沖因數(shù),振動加速度總級,期望,方差,標準差,三階累量,四階累量,有效值;

    (2)以2 048個相鄰樣本點為一個樣本段,并以70個樣本點為步長,將樣本段逐步向后平移,可在16 384個樣本點內(nèi)得到約205個樣本段;

    (3)每一工況下隨機選擇100個樣本段,在每個樣本段內(nèi)計算選定的十個特征量。

    本文是為了驗證基于變點分析思想估計本征維的可行性,及比較與現(xiàn)有方法的效果好壞,故只選取10個常用的特征量和隨機抽取了100個樣本段作為舉例,實際上可以選擇更多的特征量和樣本段。

    預(yù)處理后,在每一工況下的每一通道上可得到一個10*100噪聲數(shù)據(jù)。其中,10是維數(shù),100是樣本量。則在八種工況下,每一條通道上可得到一個10*800數(shù)據(jù)。

    2.2 計算結(jié)果

    對2.1節(jié)的試驗數(shù)據(jù),分別使用K1法、PA法、SCREE法、MAP法、PTV法及meanECP法進行處理,得到噪聲數(shù)據(jù)在19條通道上的本征維,及不同方法的程序運行時間(程序運行平臺為Matlab R2009a),結(jié)果如表2和表3所示。

    表2 本征維估計結(jié)果Tab.2 Latent characteristic variables dimensionality

    表3 不同估計方法運算時間(單位:秒)Tab.3 Computing time of different methods(Unit:s)

    續(xù)表3

    表3中最后一行“時間”表示不同方法在19條通道上程序執(zhí)行時間的平均值。為直觀比較不同方法效果,圖2~7以第7通道為例,顯示了6種方法的本征維估計結(jié)果。其中PA法以蒙特卡洛模擬的特征值均值線和0.95百分位線為標準[25],兩個標準估計的結(jié)果相同;MAP法中原始MAP指的是Velicer于1976年提出的以偏相關(guān)系數(shù)平方的最小平均值為標準,改進MAP是指Velicer于2000年提出的以偏相關(guān)系數(shù)四次方的最小平均值為標準;PTV法中以0.95為標準[15]。

    圖2 第七通道上K1法估計結(jié)果圖Fig.2 Evaluating result of K1 method in seventh sensor

    圖3 第七通道上PA法估計結(jié)果圖Fig.3 Evaluating result of PA method in seventh sensor

    圖4 第七通道上SCREE法估計結(jié)果圖Fig.4 Evaluating result of SCREE method in seventh sensor

    圖5 第七通道上MAP法估計結(jié)果圖Fig.5 Evaluating result of MAP method in seventh sensor

    圖6 第七通道上PTV法估計結(jié)果圖Fig.6 Evaluating result of PTV method in seventh sensor

    圖7 第七通道上meanECP法估計結(jié)果圖Fig.7 Evaluating result of meanECP method in seventh sensor

    2.3 比較分析

    表2中結(jié)果顯示,MAP方法估計的本征維較其他方法大,且在不同通道上的變化波動也較大;同時,第18、19通道上估計的本征維為0,說明MAP方法對這兩通道不適用;由K1法、PA法、SCREE法、PTV法的計算結(jié)果可知,在所有19條通道上的某雙層圓柱殼體機械噪聲數(shù)據(jù)的本征維在2與3范圍內(nèi),由meanECP法估計的本征維同樣也在2、3兩位置浮動(只在第18通道上估計結(jié)果為4),從而驗證了meanECP的可行性。

    由表3可知,PA方法的程序運算時間明顯高于其他方法,說明該法不適宜處理大樣本或高維數(shù)的觀測數(shù)據(jù)。同時,由于SCREE法與PTV法受主觀影響較大。因此,對于某雙層圓柱殼體機械噪聲數(shù)據(jù)來說K1法與meanECP法較合適。

    為進一步比較K1法與ECP法的效果,下面將利用判別分析方法確定噪聲數(shù)據(jù)在兩種方法估計的本征維下對每一條通道上來自不同工況樣本的判別分析結(jié)果?;舅枷胧牵合薅骋煌ǖ?,對該通道上來自八個不同工況的100個觀測樣本組成的10*800數(shù)據(jù)矩陣進行主成分分析,并分別利用K1法和meanECP法估計本征維,然后根據(jù)估計的本征維對來自不同工況的樣本利用判別分析進行分類,依據(jù)分類正確率的大小比較估計效果:分類正確率越高,說明估計的本征維越合理,估計方法效果越好。19條通道上的判別分析結(jié)果如表4所示(程序運行平臺:SPSS 17.0)。

    表4 K1法與meanECP法的判別分析結(jié)果比較Tab.4 Results of discrimnant analysis between the K1 method and the meanECP method

    表4中,除第17、18通道外,其余通道上均顯示meanECP法的估計效果比K1法的估計效果好。

    3 結(jié) 論

    基于均值估計的特征值變點法是在現(xiàn)有常用本征維估計方法基礎(chǔ)上,充分利用變點分析思想提出來的。通過對現(xiàn)有變點位置估計技術(shù)的改進,基于均值估計提出了一種針對特征值曲線特點的變點估計方法。通過在某雙層圓柱殼體機械噪聲數(shù)據(jù)中的實際應(yīng)用,驗證了基于均值估計的特征值變點法(meanECP法)在高維數(shù)據(jù)本征維估計中的可行性。最后通過分析,與現(xiàn)有常用的本征維估計方法進行了比較。結(jié)果顯示:meanECP法理論簡單、計算簡便,計算結(jié)果客觀、合理,對某雙層圓柱殼體機械噪聲數(shù)據(jù)具有較好的處理效果。

    參 考 文 獻:

    [1]van der Maaten L J P,Postma E O,van den Herik H J.Dimensionality reduction:a comparative review[J].IEEE Transactions on Pattern Analysis and Machine Intelligence,2009(10):1-35.

    [2]Rubén D L,Pedro V.Determining the number of factors to retain in EFA:An easy-to-use computer program for carrying out Parallel Analysis[J].Practical Assessment,Research&Evaluation,2007,12(2):1-11.

    [3]Zwick W R,Velicer W F.Comparison of five rules for determining the number of components to retain[J].Psychological Bulletin,1986,99(3):432-442.

    [4]Camastra F,Vinciarelli A.Estimating the intrinsic dimension of data with a fractal-based method[J].IEEE Trans.Pattern Anal Mach Intell,2002,24(10):1404-1407.

    [5]Kaiser H F.The application of electronic computers to factor analysis[J].Educational and Psychological Measurement, 1960,20:141-151.

    [6]Horn J L.A rationale and test for the number of factors in factor analysis[J].Psychometrika,1965,30(2):179-185.

    [7]Cattell R B.The scree test for the number of factors[J].Multivariate Behavioral Research,1966(1):245-276.

    [8]Velicer W F.Determining the number of components from the matrix of partial correlations[J].Psychometrika,1976,41 (3):321-327.

    [9]Velicer W F,Eaton C A,Fava J L.Construct explication through factor or component analysis:A review and evaluation of alternative procedures for determining the number of factors or components[M].Problems and solutions in human assessment:Honoring Douglas N.Jackson at seventy,Boston:Kluwer Academic Publishers,2000:41-71.

    [10]Jolliffe I T.Principal component analysis(second edition)[M].New York:Springer-Verlag,2002.

    [11]Fabrigar L R,Wegener D T,Maccallum R C,et al.Evaluating the use of exploratory factor analysis in psychological research[J].Psychological Methods,1999(3):272-299.

    [12]Hayton J C,Allen D G,Scarpello V.Factor retention decisions in exploratory factor analysis:a tutorial on parallel analysis[J].Organizational Research Methods,2004,7:191-205.

    [13]Streiner D L.Figuring out factors:The use and misuse of factor analysis[J].Canadian Journal of Psychiatry,1994,39: 135-140.

    [14]Floyd F J,Widaman K F.Factor analysis in the development and refinement of clinical assessment instruments[J].Psychological Assessment,1995,7:286-299.

    [15]Jackson D A.Stopping rules in principal components analysis:a comparison of heuristical and statistical approaches[J]. Ecology,1993,74(8):2204-2214.

    [16]Costello A B,Osborne J W.Best practices in exploratory factor analysis:four recommendations for getting the most from your analysis[J].Practical Assessment Research&Evaluation,2005,10(7):1-9.

    [17]Andreou E,Ghysels E.Structural breaks in financial time series[M].Handbook of Financial Time Series,Berlin Heidelberg:Springer-Verlag,2009:6,839-870.

    [18]關(guān)文韜,金百鎖,繆柏其.原油價格與黃金價格的變點分析[J].中國科學(xué)技術(shù)大學(xué)學(xué)報,2014,44(6):502-507. Guan Wentao,Jin Baisuo,Miao Baiqi.Change-point analysis of crude oil and gold prices[J].Journal of University of Science and Technology of China,2014,44(6):502-507.

    [19]熊立華,于坤霞,董磊華等.水文時間序列變點分析的可靠性檢驗[J].武漢大學(xué)學(xué)報(工學(xué)版),2011,44(2):137-141. Xiong Lihua,Yu Kunxia,Dong Leihua,et al.Reliability test for detecting change point of hydrological time series[J]. Engineering Journal of Wuhan University,2011,44(2):137-141.

    [20]危潤初,肖長來,方 樟.黑龍江建三江地區(qū)地下水動態(tài)趨勢突變點分析[J].吉林大學(xué)學(xué)報:地球科學(xué)版,2016,46 (1):202-210. Wei Runchu,Xiao Changlai,Fang Zhang.Trends mutation nodes of groundwater dynamic in jiangsanjiang area of heilongjiang province[J].Journal of Jilin University(Earth Science Edition),2016,46(1):202-210.

    [21]Céline L,Fran?ois R.Detection and localization of change-points in high-dimensional network traffic data[J].The Annals of Applied Statistics,2009,3(2):637-662.

    [22]陳希孺.變點統(tǒng)計分析簡介(I)[J].數(shù)理統(tǒng)計與管理,1991,10(1):55-58. Chen Xiru.The overview of statistical analysis and case study of change points(I)[J].Journal of Applied Statistics and Management,1991,10(1):55-58.

    [23]Csorgo M,Horvath L.Limit theorems in change-point analysis[M].London:John Wiley&Sons,1997.

    [24]張志華,梁勝杰,尹曰建等.基于SVM的雙層圓柱殼體機械振動噪聲數(shù)據(jù)特征提取方法比較[J].船舶力學(xué),2015,19 (1-2):215-220. Zhang Zhihua,Liang Shengjie,Yin Yuejian,et al.Comparison of feature extraction methods on mechanical vibrating noise of ribbed cylindrical double-shells based SVM[J].Journal of Ship Mechanics,2015,19(1-2):215-220.

    [25]Raiche G,Riopel M,Blais J G.Non Graphical Solutions for the Cattel’s Scree Test[C]//International Meeting of the Psychometric Society(IMPS 2006).HEC,Montreal,2006.

    A method of latent characteristic variables dimensionality evaluating of High-dimensional mechanical noise data based Change-Point analysis

    LIANG Sheng-jie1,ZHANG Zhi-hua2,GAO Shao-zhong1,HU Jun-bo2,LI Da-wei3
    (1.91697 PLA troops,Qingdao 266405,China;2.Dept.of Scientific Research,Naval University of Engineering, Wuhan 430033,China;3.91550 PLA troops,Dalian 116023,China)

    According to eigenvalue-curve’s characteristic of covariance matrix from a High-dimensional data,a new method of latent characteristic variables dimensionality evaluating is proposed based change-point analysis technology,which called Eigenvalue Change-point method based mean estimating,and then apply it to deal with the mechanical noise data of ribbed cylindrical double-shells.Through the analysis comparing with other methods,the proposed method was validated in its feasibility.It is also verified that Eigenvalue Change-point method based mean estimating is simply in theory and convenient in computing;its compute result is impersonal and reasonable;meanwhile,it has attained a preferable result when dealing with the mechanical noise data of a ribbed cylindrical double-shells.

    High-dimensional data;latent characteristic variables dimensionality;Change-Point Analysis; Discriminant Analysis;mechanical noise data

    TB53 O212.4

    A

    10.3969/j.issn.1007-7294.2016.07.015

    1007-7294(2016)11-1485-09

    2016-07-20

    國家自然科學(xué)基金(50775218)

    梁勝杰(1981-),男,博士,E-mail:lsj990308@163.com;張志華(1965-),男,教授,博士生導(dǎo)師;高紹忠(1963-),男,高級工程師;胡俊波(1982-),男,博士;李大偉(1985-),男,博士。

    猜你喜歡
    變點本征高維
    基于本征正交分解的水平軸風(fēng)力機非定常尾跡特性分析
    回歸模型參數(shù)的變點檢測方法研究
    正態(tài)分布序列均值變點檢測的貝葉斯方法
    基于二元分割的多變點估計
    KP和mKP可積系列的平方本征對稱和Miura變換
    獨立二項分布序列變點的識別方法
    一種改進的GP-CLIQUE自適應(yīng)高維子空間聚類算法
    本征平方函數(shù)在變指數(shù)Herz及Herz-Hardy空間上的有界性
    基于加權(quán)自學(xué)習(xí)散列的高維數(shù)據(jù)最近鄰查詢算法
    一般非齊次非線性擴散方程的等價變換和高維不變子空間
    欧美黄色淫秽网站| 99精品在免费线老司机午夜| 色老头精品视频在线观看| 丝瓜视频免费看黄片| 在线永久观看黄色视频| 人成视频在线观看免费观看| 成人国语在线视频| 老司机深夜福利视频在线观看| 丝瓜视频免费看黄片| 欧美激情极品国产一区二区三区| a级片在线免费高清观看视频| 国产av精品麻豆| 99热国产这里只有精品6| 国产精品一区二区免费欧美| 亚洲精品一卡2卡三卡4卡5卡| 极品少妇高潮喷水抽搐| 黄色视频,在线免费观看| 国产欧美日韩综合在线一区二区| 中文字幕高清在线视频| 免费在线观看影片大全网站| 国产午夜精品久久久久久| 两性夫妻黄色片| 校园春色视频在线观看| 欧美精品啪啪一区二区三区| 天堂俺去俺来也www色官网| 成年动漫av网址| 国产精品国产av在线观看| 国产有黄有色有爽视频| 亚洲第一欧美日韩一区二区三区| 老司机深夜福利视频在线观看| 精品亚洲成国产av| 狠狠婷婷综合久久久久久88av| x7x7x7水蜜桃| 一本综合久久免费| tube8黄色片| 99精国产麻豆久久婷婷| 十分钟在线观看高清视频www| 色综合婷婷激情| 三上悠亚av全集在线观看| 91精品三级在线观看| 亚洲精品自拍成人| 人人澡人人妻人| 日本黄色视频三级网站网址 | 性少妇av在线| 欧美乱色亚洲激情| 国产真人三级小视频在线观看| 日韩欧美免费精品| 国产激情久久老熟女| 欧美日韩瑟瑟在线播放| 国产免费av片在线观看野外av| 母亲3免费完整高清在线观看| 婷婷丁香在线五月| 欧美精品亚洲一区二区| 免费在线观看完整版高清| 在线观看午夜福利视频| 亚洲精品久久午夜乱码| 亚洲精品在线美女| 一本一本久久a久久精品综合妖精| 超色免费av| 制服人妻中文乱码| 亚洲国产欧美日韩在线播放| 99国产精品免费福利视频| 两性夫妻黄色片| 美女扒开内裤让男人捅视频| 亚洲国产中文字幕在线视频| 久久国产精品男人的天堂亚洲| 麻豆乱淫一区二区| 97人妻天天添夜夜摸| 午夜两性在线视频| 少妇粗大呻吟视频| 在线观看www视频免费| 亚洲 国产 在线| 亚洲国产看品久久| 高清毛片免费观看视频网站 | svipshipincom国产片| 水蜜桃什么品种好| 久久天堂一区二区三区四区| 欧美亚洲 丝袜 人妻 在线| 国产高清国产精品国产三级| 亚洲成人免费电影在线观看| 女人久久www免费人成看片| 国产成人av教育| 热re99久久国产66热| 久久精品国产a三级三级三级| 色94色欧美一区二区| 亚洲人成伊人成综合网2020| 757午夜福利合集在线观看| av欧美777| 国产精品 国内视频| 丰满人妻熟妇乱又伦精品不卡| 在线观看免费视频日本深夜| 深夜精品福利| 亚洲欧美激情综合另类| 无限看片的www在线观看| 亚洲国产毛片av蜜桃av| 久久久久国内视频| 一级毛片精品| 国产精品久久久av美女十八| 在线观看免费日韩欧美大片| 色综合婷婷激情| 波多野结衣av一区二区av| 中文字幕另类日韩欧美亚洲嫩草| 亚洲欧洲精品一区二区精品久久久| 在线观看免费视频日本深夜| 成人亚洲精品一区在线观看| 国产成人欧美在线观看 | 青草久久国产| 美国免费a级毛片| 天天影视国产精品| 欧美成狂野欧美在线观看| 亚洲色图综合在线观看| 亚洲少妇的诱惑av| 日韩人妻精品一区2区三区| 久久ye,这里只有精品| 亚洲精品国产区一区二| 少妇裸体淫交视频免费看高清 | 亚洲欧美日韩高清在线视频| 国产亚洲精品第一综合不卡| 久久天躁狠狠躁夜夜2o2o| 久久国产精品人妻蜜桃| 成人av一区二区三区在线看| 性色av乱码一区二区三区2| 欧美另类亚洲清纯唯美| 飞空精品影院首页| 黄频高清免费视频| 久久热在线av| 亚洲美女黄片视频| 又大又爽又粗| 亚洲av成人一区二区三| 国产精品国产av在线观看| 国产av又大| 美女福利国产在线| 80岁老熟妇乱子伦牲交| 久久久久国内视频| 亚洲熟女精品中文字幕| 男女午夜视频在线观看| 纯流量卡能插随身wifi吗| 女人被躁到高潮嗷嗷叫费观| av中文乱码字幕在线| 女人高潮潮喷娇喘18禁视频| 亚洲av成人不卡在线观看播放网| 久久久久久人人人人人| 精品国产超薄肉色丝袜足j| 在线视频色国产色| 一a级毛片在线观看| 久久久久久人人人人人| 大型黄色视频在线免费观看| 999久久久国产精品视频| 国产一区有黄有色的免费视频| 999久久久国产精品视频| 国产不卡一卡二| 亚洲,欧美精品.| 国产欧美日韩精品亚洲av| 一区福利在线观看| 巨乳人妻的诱惑在线观看| 久久精品国产清高在天天线| 亚洲第一青青草原| 色94色欧美一区二区| 十八禁人妻一区二区| 久久久国产一区二区| 亚洲欧美一区二区三区久久| 欧美日韩黄片免| 美女午夜性视频免费| 高清视频免费观看一区二区| 人妻 亚洲 视频| 日韩欧美一区视频在线观看| 亚洲欧美色中文字幕在线| a级片在线免费高清观看视频| 成年人免费黄色播放视频| 黑丝袜美女国产一区| 亚洲精品自拍成人| 欧美一级毛片孕妇| 午夜日韩欧美国产| 国产亚洲一区二区精品| 日韩 欧美 亚洲 中文字幕| 午夜免费观看网址| 一本一本久久a久久精品综合妖精| 制服诱惑二区| 法律面前人人平等表现在哪些方面| 色综合婷婷激情| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品美女久久久久99蜜臀| 欧美精品人与动牲交sv欧美| 国产三级黄色录像| 少妇的丰满在线观看| 亚洲午夜理论影院| 大型黄色视频在线免费观看| 女性生殖器流出的白浆| 国产欧美日韩综合在线一区二区| 搡老熟女国产l中国老女人| 巨乳人妻的诱惑在线观看| 男女免费视频国产| 国产激情久久老熟女| 91在线观看av| 三上悠亚av全集在线观看| 视频区欧美日本亚洲| 久久久久国产一级毛片高清牌| 欧美性长视频在线观看| 精品人妻在线不人妻| 欧美日韩成人在线一区二区| 欧美最黄视频在线播放免费 | 99国产精品一区二区蜜桃av | 涩涩av久久男人的天堂| 亚洲七黄色美女视频| 久久影院123| 欧美亚洲日本最大视频资源| 亚洲成av片中文字幕在线观看| 中文字幕人妻熟女乱码| 丰满迷人的少妇在线观看| 777久久人妻少妇嫩草av网站| 人妻一区二区av| 午夜福利在线观看吧| 正在播放国产对白刺激| 老司机亚洲免费影院| 老熟女久久久| 99re6热这里在线精品视频| 亚洲精品粉嫩美女一区| 国产精品秋霞免费鲁丝片| 亚洲欧美一区二区三区久久| 国产一区二区三区视频了| 亚洲精品在线观看二区| 十八禁网站免费在线| 老司机午夜十八禁免费视频| 亚洲三区欧美一区| 99精国产麻豆久久婷婷| 亚洲在线自拍视频| 久久精品aⅴ一区二区三区四区| 男人操女人黄网站| 侵犯人妻中文字幕一二三四区| 久久久久精品人妻al黑| 国产成人精品久久二区二区91| 伊人久久大香线蕉亚洲五| 午夜福利,免费看| 亚洲五月天丁香| 欧美日韩亚洲国产一区二区在线观看 | 视频区欧美日本亚洲| 色精品久久人妻99蜜桃| 成年人黄色毛片网站| 一级片'在线观看视频| 一本大道久久a久久精品| 亚洲av成人不卡在线观看播放网| 99在线人妻在线中文字幕 | 日韩欧美在线二视频 | 国产欧美亚洲国产| 18禁裸乳无遮挡免费网站照片 | av片东京热男人的天堂| 黄片播放在线免费| 成人精品一区二区免费| 少妇 在线观看| 久久精品国产a三级三级三级| 男女床上黄色一级片免费看| 一级片免费观看大全| 国产精品亚洲av一区麻豆| 久久精品国产综合久久久| www.自偷自拍.com| 亚洲精品久久成人aⅴ小说| 老司机影院毛片| 日本一区二区免费在线视频| 日日摸夜夜添夜夜添小说| 身体一侧抽搐| 国产不卡av网站在线观看| 性色av乱码一区二区三区2| 最新在线观看一区二区三区| 成人国语在线视频| 亚洲精品av麻豆狂野| 国产无遮挡羞羞视频在线观看| 无遮挡黄片免费观看| 精品人妻熟女毛片av久久网站| 亚洲一码二码三码区别大吗| 麻豆国产av国片精品| 日韩欧美在线二视频 | 国产主播在线观看一区二区| 国产av又大| 久久精品国产综合久久久| 乱人伦中国视频| 久久中文字幕人妻熟女| 午夜精品久久久久久毛片777| 国产男女内射视频| 国产精品国产高清国产av | 亚洲精华国产精华精| 狂野欧美激情性xxxx| 老司机亚洲免费影院| 亚洲精品国产区一区二| 亚洲精品久久成人aⅴ小说| 国产精品av久久久久免费| 国产亚洲欧美98| 亚洲精品av麻豆狂野| 一区二区三区激情视频| 国产av一区二区精品久久| 亚洲精品粉嫩美女一区| 一级,二级,三级黄色视频| 久久久久国产一级毛片高清牌| 国产成+人综合+亚洲专区| 久久精品国产亚洲av香蕉五月 | 亚洲 欧美一区二区三区| 国产亚洲av高清不卡| 午夜免费观看网址| 久久精品91无色码中文字幕| 国产激情久久老熟女| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美黄色片欧美黄色片| 亚洲aⅴ乱码一区二区在线播放 | 欧美亚洲日本最大视频资源| 亚洲第一av免费看| 久久亚洲真实| 欧美成人免费av一区二区三区 | 一区二区三区激情视频| 国产一区二区三区综合在线观看| 亚洲五月天丁香| 国产又爽黄色视频| 欧美日韩亚洲综合一区二区三区_| 女人精品久久久久毛片| 人人澡人人妻人| 日韩中文字幕欧美一区二区| 亚洲国产看品久久| 1024视频免费在线观看| √禁漫天堂资源中文www| 国产视频一区二区在线看| 在线免费观看的www视频| 又大又爽又粗| 最新美女视频免费是黄的| 不卡av一区二区三区| 日韩精品免费视频一区二区三区| 91av网站免费观看| 正在播放国产对白刺激| 美女视频免费永久观看网站| 亚洲熟妇中文字幕五十中出 | 俄罗斯特黄特色一大片| 精品久久久久久电影网| 不卡av一区二区三区| 在线永久观看黄色视频| 亚洲一区高清亚洲精品| 日本五十路高清| 1024视频免费在线观看| 淫妇啪啪啪对白视频| 12—13女人毛片做爰片一| 亚洲成人国产一区在线观看| 亚洲av欧美aⅴ国产| 国产区一区二久久| 真人做人爱边吃奶动态| 老熟妇仑乱视频hdxx| 捣出白浆h1v1| 亚洲欧美激情综合另类| 9191精品国产免费久久| 悠悠久久av| 欧美激情 高清一区二区三区| 久9热在线精品视频| 一级毛片精品| 啦啦啦 在线观看视频| xxx96com| 国产成人系列免费观看| 性色av乱码一区二区三区2| 久久午夜综合久久蜜桃| 精品一区二区三区视频在线观看免费 | 一级,二级,三级黄色视频| 高清在线国产一区| 国产视频一区二区在线看| 在线免费观看的www视频| 777米奇影视久久| 免费黄频网站在线观看国产| 国产精品偷伦视频观看了| 韩国av一区二区三区四区| 亚洲国产看品久久| 精品福利永久在线观看| 午夜久久久在线观看| 国产精品亚洲一级av第二区| 免费在线观看完整版高清| 99riav亚洲国产免费| 一级毛片精品| 在线视频色国产色| 国产精品 国内视频| 亚洲精品av麻豆狂野| 老熟妇仑乱视频hdxx| 日韩三级视频一区二区三区| 女警被强在线播放| 香蕉久久夜色| 精品乱码久久久久久99久播| 我的亚洲天堂| 亚洲熟妇中文字幕五十中出 | 欧美日韩一级在线毛片| 国产激情欧美一区二区| 午夜影院日韩av| 久久久久久久久久久久大奶| xxxhd国产人妻xxx| 免费黄频网站在线观看国产| 黄频高清免费视频| 亚洲综合色网址| 在线观看www视频免费| 欧美日韩亚洲综合一区二区三区_| 真人做人爱边吃奶动态| 叶爱在线成人免费视频播放| 人妻丰满熟妇av一区二区三区 | 久久久国产精品麻豆| 免费不卡黄色视频| 老司机深夜福利视频在线观看| 国产xxxxx性猛交| 免费观看a级毛片全部| av一本久久久久| 精品一区二区三区视频在线观看免费 | 人人妻人人澡人人看| 亚洲av片天天在线观看| 国产一区有黄有色的免费视频| 一级,二级,三级黄色视频| 欧美最黄视频在线播放免费 | 国产99久久九九免费精品| 精品一区二区三卡| 国产色视频综合| 一级毛片精品| 成年动漫av网址| 日韩有码中文字幕| 精品久久久久久久久久免费视频 | 亚洲色图综合在线观看| 巨乳人妻的诱惑在线观看| 别揉我奶头~嗯~啊~动态视频| 热99国产精品久久久久久7| 国产欧美日韩精品亚洲av| 人人妻人人澡人人看| 免费在线观看视频国产中文字幕亚洲| av视频免费观看在线观看| 国产欧美日韩综合在线一区二区| 免费日韩欧美在线观看| 亚洲中文av在线| 亚洲avbb在线观看| 飞空精品影院首页| 天堂中文最新版在线下载| 久久久久视频综合| 亚洲熟女毛片儿| 国产一区二区三区在线臀色熟女 | bbb黄色大片| 亚洲色图av天堂| 五月开心婷婷网| 三级毛片av免费| 无限看片的www在线观看| 国产蜜桃级精品一区二区三区 | 国产精品亚洲一级av第二区| 欧美日韩亚洲综合一区二区三区_| 动漫黄色视频在线观看| 高清黄色对白视频在线免费看| xxxhd国产人妻xxx| 人妻 亚洲 视频| 三级毛片av免费| 欧美日韩av久久| 精品乱码久久久久久99久播| 国产亚洲欧美98| 国产精品乱码一区二三区的特点 | 国产高清videossex| 老司机靠b影院| 免费人成视频x8x8入口观看| 一区二区三区国产精品乱码| 日韩欧美一区视频在线观看| a在线观看视频网站| 久久热在线av| 嫩草影视91久久| 欧美色视频一区免费| 久久香蕉精品热| 亚洲人成77777在线视频| 黑人猛操日本美女一级片| 99国产精品一区二区三区| 国产日韩欧美亚洲二区| 在线观看午夜福利视频| 亚洲专区字幕在线| 日本精品一区二区三区蜜桃| 国产成人免费观看mmmm| 身体一侧抽搐| 超色免费av| 欧美人与性动交α欧美软件| 国产精品九九99| 国产在线精品亚洲第一网站| 国产一区在线观看成人免费| 欧美人与性动交α欧美精品济南到| av免费在线观看网站| 亚洲性夜色夜夜综合| 50天的宝宝边吃奶边哭怎么回事| 天天躁狠狠躁夜夜躁狠狠躁| 一级片免费观看大全| 一本大道久久a久久精品| 视频区欧美日本亚洲| 在线播放国产精品三级| 免费黄频网站在线观看国产| 99香蕉大伊视频| 国产精品1区2区在线观看. | 在线观看免费午夜福利视频| 国产亚洲一区二区精品| 在线观看午夜福利视频| 欧美国产精品va在线观看不卡| 亚洲久久久国产精品| 法律面前人人平等表现在哪些方面| 国产精品久久久av美女十八| 人人妻人人澡人人看| 国产视频一区二区在线看| 成在线人永久免费视频| av国产精品久久久久影院| 午夜福利乱码中文字幕| e午夜精品久久久久久久| 黄色毛片三级朝国网站| 婷婷精品国产亚洲av在线 | 啦啦啦在线免费观看视频4| 国产免费av片在线观看野外av| 欧美一级毛片孕妇| 久久人人97超碰香蕉20202| 亚洲一卡2卡3卡4卡5卡精品中文| av线在线观看网站| 亚洲一区高清亚洲精品| 欧美av亚洲av综合av国产av| 啦啦啦 在线观看视频| 一级作爱视频免费观看| 亚洲熟女毛片儿| 亚洲五月婷婷丁香| 欧美日韩黄片免| 波多野结衣一区麻豆| 麻豆成人av在线观看| 1024视频免费在线观看| 亚洲欧美日韩另类电影网站| 久久天躁狠狠躁夜夜2o2o| 变态另类成人亚洲欧美熟女 | 国产视频一区二区在线看| 美女 人体艺术 gogo| 精品视频人人做人人爽| 啦啦啦在线免费观看视频4| 在线看a的网站| 美女视频免费永久观看网站| 99久久综合精品五月天人人| 亚洲精品美女久久久久99蜜臀| 亚洲自偷自拍图片 自拍| 91国产中文字幕| 久久久国产欧美日韩av| 久久人人爽av亚洲精品天堂| 中文字幕另类日韩欧美亚洲嫩草| 精品国产亚洲在线| 可以免费在线观看a视频的电影网站| 人妻一区二区av| 国产免费男女视频| 国产av又大| 日本欧美视频一区| 国产精品电影一区二区三区 | 免费在线观看影片大全网站| 成人三级做爰电影| 国产精品九九99| 国产乱人伦免费视频| 国产欧美日韩综合在线一区二区| 一区二区日韩欧美中文字幕| 在线观看免费高清a一片| 国产精品 国内视频| 日韩欧美国产一区二区入口| www.熟女人妻精品国产| 成在线人永久免费视频| 一区二区三区激情视频| 成年女人毛片免费观看观看9 | 国产高清videossex| 天天添夜夜摸| 动漫黄色视频在线观看| 777米奇影视久久| 国产欧美亚洲国产| 欧美精品亚洲一区二区| 很黄的视频免费| 黄片小视频在线播放| 国产精品久久久久成人av| 久久性视频一级片| 欧美在线黄色| 99国产精品99久久久久| 亚洲精品自拍成人| 亚洲成人免费av在线播放| 久久热在线av| 久久精品国产亚洲av高清一级| 午夜视频精品福利| 国产精品欧美亚洲77777| 亚洲人成77777在线视频| 亚洲成人免费电影在线观看| 可以免费在线观看a视频的电影网站| 久久久久精品人妻al黑| 亚洲色图av天堂| 亚洲欧美激情综合另类| 美女高潮喷水抽搐中文字幕| 欧美日韩瑟瑟在线播放| 国产成人av教育| 欧美老熟妇乱子伦牲交| 免费日韩欧美在线观看| 18禁黄网站禁片午夜丰满| 午夜免费观看网址| 99国产精品一区二区三区| 在线观看免费午夜福利视频| 他把我摸到了高潮在线观看| 中文字幕精品免费在线观看视频| 午夜精品久久久久久毛片777| 久热这里只有精品99| 在线国产一区二区在线| 国产色视频综合| 国产淫语在线视频| 王馨瑶露胸无遮挡在线观看| 久久九九热精品免费| 午夜福利欧美成人| 99久久99久久久精品蜜桃| 两性夫妻黄色片| 午夜福利乱码中文字幕| 香蕉久久夜色| 在线观看日韩欧美| 国产av一区二区精品久久| 日韩免费高清中文字幕av| 一级作爱视频免费观看| 热99久久久久精品小说推荐| 中文字幕精品免费在线观看视频| 在线观看免费午夜福利视频| 国产不卡av网站在线观看| av不卡在线播放| 热99re8久久精品国产| 宅男免费午夜| a级毛片黄视频| 两人在一起打扑克的视频| 天堂√8在线中文| 黄色成人免费大全| 亚洲成a人片在线一区二区| 国产一区二区三区综合在线观看| 久久国产精品影院| 大码成人一级视频| 午夜福利在线免费观看网站|