• <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ù)最近鄰查詢算法
    一般非齊次非線性擴散方程的等價變換和高維不變子空間
    欧美中文综合在线视频| 国产成人精品无人区| 人人妻,人人澡人人爽秒播| 一级毛片高清免费大全| 日本撒尿小便嘘嘘汇集6| 亚洲精品粉嫩美女一区| 亚洲精品在线观看二区| 欧美激情久久久久久爽电影 | 成熟少妇高潮喷水视频| 俄罗斯特黄特色一大片| 日韩欧美国产一区二区入口| 精品高清国产在线一区| aaaaa片日本免费| 欧美乱妇无乱码| 丰满的人妻完整版| 777久久人妻少妇嫩草av网站| 一级片'在线观看视频| 一级,二级,三级黄色视频| 欧美激情 高清一区二区三区| 日韩中文字幕欧美一区二区| 久久午夜综合久久蜜桃| 大香蕉久久成人网| 国产不卡av网站在线观看| 精品国产亚洲在线| 国产人伦9x9x在线观看| 最新的欧美精品一区二区| 精品乱码久久久久久99久播| 一进一出抽搐gif免费好疼 | 久久久国产一区二区| 久久精品成人免费网站| 午夜福利,免费看| 婷婷成人精品国产| 一级毛片高清免费大全| 国产精品亚洲av一区麻豆| 男女下面插进去视频免费观看| 国产精品香港三级国产av潘金莲| 国产精品免费视频内射| 国产精品一区二区免费欧美| 欧美人与性动交α欧美软件| 老熟妇乱子伦视频在线观看| 男人操女人黄网站| 久久久久精品人妻al黑| 午夜免费成人在线视频| 亚洲国产精品合色在线| 国产高清videossex| 亚洲精华国产精华精| 美女高潮喷水抽搐中文字幕| 欧美日韩亚洲高清精品| 天天躁日日躁夜夜躁夜夜| 欧美另类亚洲清纯唯美| 久久性视频一级片| 久久久久久久午夜电影 | 在线观看66精品国产| 黄色怎么调成土黄色| av片东京热男人的天堂| 狠狠婷婷综合久久久久久88av| tube8黄色片| 午夜两性在线视频| 久久人人爽av亚洲精品天堂| 精品无人区乱码1区二区| 香蕉丝袜av| 国产av一区二区精品久久| 久久久久久免费高清国产稀缺| 精品国内亚洲2022精品成人 | 国产高清国产精品国产三级| 久久性视频一级片| 久久香蕉国产精品| 天天影视国产精品| 建设人人有责人人尽责人人享有的| 怎么达到女性高潮| av片东京热男人的天堂| 亚洲精品自拍成人| 黄色视频不卡| 国产一区二区激情短视频| 国产成人精品在线电影| 天堂中文最新版在线下载| 亚洲一区二区三区欧美精品| 国产精品久久久久久人妻精品电影| 国产精品久久视频播放| 国产成人av教育| 伦理电影免费视频| 天天躁日日躁夜夜躁夜夜| 天堂俺去俺来也www色官网| 亚洲成人免费av在线播放| 美女午夜性视频免费| 捣出白浆h1v1| 精品一品国产午夜福利视频| av在线播放免费不卡| 好看av亚洲va欧美ⅴa在| 亚洲 欧美一区二区三区| 国产精品久久久久久精品古装| 亚洲人成电影观看| 午夜福利在线免费观看网站| 国产免费av片在线观看野外av| av中文乱码字幕在线| 成人特级黄色片久久久久久久| 国产在视频线精品| www日本在线高清视频| 99精品在免费线老司机午夜| 亚洲精品国产色婷婷电影| 午夜免费鲁丝| 水蜜桃什么品种好| 午夜免费成人在线视频| 嫩草影视91久久| 久久国产精品影院| 亚洲一区二区三区不卡视频| 中文字幕高清在线视频| 正在播放国产对白刺激| 熟女少妇亚洲综合色aaa.| 国产成人免费观看mmmm| 国产欧美日韩综合在线一区二区| 免费在线观看亚洲国产| 精品国产国语对白av| 侵犯人妻中文字幕一二三四区| 国产xxxxx性猛交| 视频区欧美日本亚洲| 交换朋友夫妻互换小说| 免费不卡黄色视频| 亚洲精品久久午夜乱码| 亚洲专区国产一区二区| 99国产极品粉嫩在线观看| 午夜福利视频在线观看免费| 亚洲av成人av| 捣出白浆h1v1| 老熟妇仑乱视频hdxx| 搡老岳熟女国产| 俄罗斯特黄特色一大片| 伦理电影免费视频| 宅男免费午夜| 久热爱精品视频在线9| 乱人伦中国视频| 久久中文看片网| 视频在线观看一区二区三区| 国产男女内射视频| 一级黄色大片毛片| 91精品三级在线观看| 国产精品成人在线| 久久 成人 亚洲| 午夜视频精品福利| 国产亚洲精品久久久久5区| 精品久久久久久电影网| 激情在线观看视频在线高清 | 婷婷成人精品国产| 国产精品av久久久久免费| 日本vs欧美在线观看视频| 一本大道久久a久久精品| 欧美在线黄色| 黄片大片在线免费观看| av国产精品久久久久影院| 亚洲 国产 在线| 俄罗斯特黄特色一大片| 丰满迷人的少妇在线观看| 999久久久精品免费观看国产| 12—13女人毛片做爰片一| 精品国产超薄肉色丝袜足j| 精品卡一卡二卡四卡免费| 少妇的丰满在线观看| 国产精华一区二区三区| a在线观看视频网站| 91国产中文字幕| 十八禁人妻一区二区| 桃红色精品国产亚洲av| 热re99久久国产66热| 一a级毛片在线观看| 嫩草影视91久久| 久久精品国产a三级三级三级| 欧美黄色片欧美黄色片| 人人妻人人澡人人看| 色在线成人网| 51午夜福利影视在线观看| 国产精品免费大片| 老熟妇仑乱视频hdxx| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲精品久久午夜乱码| 亚洲欧美色中文字幕在线| 最近最新免费中文字幕在线| 亚洲成人手机| 欧美午夜高清在线| 久久久久久亚洲精品国产蜜桃av| 久久香蕉精品热| 欧美激情极品国产一区二区三区| 欧美成人免费av一区二区三区 | 99精国产麻豆久久婷婷| 正在播放国产对白刺激| av线在线观看网站| 这个男人来自地球电影免费观看| 亚洲一区二区三区不卡视频| 亚洲,欧美精品.| 日本一区二区免费在线视频| 欧美日韩精品网址| videosex国产| www.熟女人妻精品国产| 欧美 日韩 精品 国产| 黄色视频,在线免费观看| 欧美日韩亚洲国产一区二区在线观看 | 亚洲国产看品久久| 精品熟女少妇八av免费久了| 成年版毛片免费区| 午夜福利免费观看在线| 国产高清videossex| 成人特级黄色片久久久久久久| 午夜福利,免费看| 国产精品久久视频播放| 免费不卡黄色视频| 国产精品美女特级片免费视频播放器 | 精品人妻熟女毛片av久久网站| 一二三四社区在线视频社区8| 中文字幕高清在线视频| 一区在线观看完整版| 午夜久久久在线观看| 搡老熟女国产l中国老女人| 丰满饥渴人妻一区二区三| 一级a爱片免费观看的视频| 亚洲精品国产精品久久久不卡| 人成视频在线观看免费观看| 日日爽夜夜爽网站| 欧美黑人精品巨大| 美女 人体艺术 gogo| 色婷婷久久久亚洲欧美| 日韩熟女老妇一区二区性免费视频| 成年女人毛片免费观看观看9 | 国产成人精品在线电影| 久久久久久久久免费视频了| 午夜福利欧美成人| 三级毛片av免费| 一二三四在线观看免费中文在| 国产欧美亚洲国产| 正在播放国产对白刺激| 99re在线观看精品视频| 人人妻人人澡人人看| 黄色成人免费大全| 国产有黄有色有爽视频| 日日夜夜操网爽| 国产主播在线观看一区二区| 99久久精品国产亚洲精品| 久久国产精品影院| 天堂中文最新版在线下载| 欧美久久黑人一区二区| 久久精品国产亚洲av香蕉五月 | 国产真人三级小视频在线观看| 成人手机av| 9色porny在线观看| 韩国av一区二区三区四区| 欧美日韩福利视频一区二区| 黄片播放在线免费| 精品国产超薄肉色丝袜足j| 国产精品九九99| 黄片大片在线免费观看| 国产精品免费大片| 自线自在国产av| av有码第一页| 久久久久久免费高清国产稀缺| 欧美日韩亚洲高清精品| 精品国产超薄肉色丝袜足j| 不卡一级毛片| 精品亚洲成a人片在线观看| 亚洲男人天堂网一区| 丰满的人妻完整版| 国产麻豆69| 一区在线观看完整版| 自拍欧美九色日韩亚洲蝌蚪91| 黑人猛操日本美女一级片| 热99国产精品久久久久久7| 手机成人av网站| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲久久久国产精品| a级毛片在线看网站| 欧美成人免费av一区二区三区 | 成人特级黄色片久久久久久久| 久久中文看片网| 国产黄色免费在线视频| 成熟少妇高潮喷水视频| 午夜福利影视在线免费观看| 国产欧美日韩一区二区三区在线| 热99久久久久精品小说推荐| 在线天堂中文资源库| 制服诱惑二区| 国产亚洲av高清不卡| 黑丝袜美女国产一区| 午夜免费鲁丝| 国产91精品成人一区二区三区| 色婷婷av一区二区三区视频| 好男人电影高清在线观看| 一边摸一边抽搐一进一小说 | 69av精品久久久久久| 黄色毛片三级朝国网站| 欧美黄色片欧美黄色片| 久久ye,这里只有精品| 久久国产乱子伦精品免费另类| 久久精品国产亚洲av高清一级| 两个人看的免费小视频| 热99久久久久精品小说推荐| 夫妻午夜视频| 亚洲国产精品sss在线观看 | 色尼玛亚洲综合影院| 新久久久久国产一级毛片| 久久性视频一级片| 9热在线视频观看99| videosex国产| 夫妻午夜视频| 欧美日韩av久久| 国产成人啪精品午夜网站| 国精品久久久久久国模美| 亚洲熟女毛片儿| 免费在线观看完整版高清| 男人操女人黄网站| 咕卡用的链子| av视频免费观看在线观看| avwww免费| 久久中文字幕人妻熟女| 国产成人欧美在线观看 | 久久精品国产亚洲av香蕉五月 | 变态另类成人亚洲欧美熟女 | 国产在视频线精品| 国产xxxxx性猛交| 黄频高清免费视频| 热re99久久国产66热| 欧美日本中文国产一区发布| 国产99久久九九免费精品| 日韩欧美一区二区三区在线观看 | 亚洲 国产 在线| 成人黄色视频免费在线看| 亚洲成人免费电影在线观看| 国产精品偷伦视频观看了| 黄色成人免费大全| 亚洲av欧美aⅴ国产| 韩国av一区二区三区四区| 久久久久久久午夜电影 | 自线自在国产av| 久久香蕉激情| 丝袜在线中文字幕| 麻豆国产av国片精品| 动漫黄色视频在线观看| 如日韩欧美国产精品一区二区三区| 亚洲专区国产一区二区| 高清黄色对白视频在线免费看| 午夜日韩欧美国产| 国产av精品麻豆| 91精品国产国语对白视频| 十八禁网站免费在线| 欧美久久黑人一区二区| 亚洲精品一二三| 大片电影免费在线观看免费| 国产主播在线观看一区二区| 国内久久婷婷六月综合欲色啪| 在线观看一区二区三区激情| 成年女人毛片免费观看观看9 | 黄色片一级片一级黄色片| 身体一侧抽搐| 操出白浆在线播放| 老熟女久久久| 日本撒尿小便嘘嘘汇集6| 久久人人97超碰香蕉20202| 精品熟女少妇八av免费久了| 亚洲熟女毛片儿| 欧美日韩乱码在线| 老司机午夜福利在线观看视频| 亚洲色图av天堂| 大香蕉久久网| 女人精品久久久久毛片| 国产视频一区二区在线看| 深夜精品福利| 欧美成人午夜精品| 久久国产精品影院| 国产三级黄色录像| 女人精品久久久久毛片| 美女国产高潮福利片在线看| 一区二区日韩欧美中文字幕| 别揉我奶头~嗯~啊~动态视频| 岛国在线观看网站| 亚洲成人国产一区在线观看| 亚洲欧美激情综合另类| 国产精品一区二区精品视频观看| 国产欧美日韩精品亚洲av| 在线av久久热| 中文亚洲av片在线观看爽 | 一级毛片精品| 精品人妻1区二区| 最新的欧美精品一区二区| 久久国产亚洲av麻豆专区| 国产不卡一卡二| 欧美日韩视频精品一区| 久久久久国产精品人妻aⅴ院 | 精品福利观看| 欧美亚洲 丝袜 人妻 在线| 日韩中文字幕欧美一区二区| 久久香蕉激情| 婷婷精品国产亚洲av在线 | 国产不卡av网站在线观看| 国产精品一区二区在线不卡| 午夜福利,免费看| 人人妻人人添人人爽欧美一区卜| 亚洲熟妇熟女久久| 中文字幕高清在线视频| 中文字幕人妻丝袜制服| 久久久久国产精品人妻aⅴ院 | ponron亚洲| 人人妻人人澡人人看| av天堂在线播放| 色婷婷av一区二区三区视频| 美女午夜性视频免费| 三级毛片av免费| 女同久久另类99精品国产91| a级片在线免费高清观看视频| 国产成人av教育| 欧美乱色亚洲激情| 国产av一区二区精品久久| 欧美人与性动交α欧美精品济南到| 无限看片的www在线观看| 亚洲av美国av| 亚洲欧美激情综合另类| 国产精品亚洲一级av第二区| 国产精品一区二区在线观看99| 精品国产美女av久久久久小说| 精品乱码久久久久久99久播| 最新在线观看一区二区三区| 国产99白浆流出| 国产亚洲精品久久久久5区| 日日爽夜夜爽网站| 久久久精品免费免费高清| 国产极品粉嫩免费观看在线| 老熟妇乱子伦视频在线观看| 亚洲av日韩精品久久久久久密| 国产日韩一区二区三区精品不卡| 国产高清激情床上av| 法律面前人人平等表现在哪些方面| 一进一出好大好爽视频| 国产激情欧美一区二区| 午夜福利在线免费观看网站| 757午夜福利合集在线观看| 久久精品91无色码中文字幕| 妹子高潮喷水视频| 视频区图区小说| 天天躁狠狠躁夜夜躁狠狠躁| 国内久久婷婷六月综合欲色啪| 国产一卡二卡三卡精品| 黄色a级毛片大全视频| 亚洲精品美女久久av网站| 啦啦啦视频在线资源免费观看| 欧美精品高潮呻吟av久久| 女人精品久久久久毛片| 国产主播在线观看一区二区| 成熟少妇高潮喷水视频| 天天影视国产精品| 多毛熟女@视频| 欧美老熟妇乱子伦牲交| 一区二区三区国产精品乱码| 在线看a的网站| 国产亚洲av高清不卡| 一级作爱视频免费观看| 色在线成人网| 久久久久久人人人人人| 精品国产一区二区久久| 精品一区二区三卡| 黑人欧美特级aaaaaa片| av有码第一页| 757午夜福利合集在线观看| 这个男人来自地球电影免费观看| 午夜精品国产一区二区电影| 精品欧美一区二区三区在线| 精品国产美女av久久久久小说| 免费在线观看日本一区| 中文字幕另类日韩欧美亚洲嫩草| 9热在线视频观看99| 国产精品成人在线| 国产男靠女视频免费网站| 免费久久久久久久精品成人欧美视频| 丰满的人妻完整版| 正在播放国产对白刺激| 久久精品国产亚洲av香蕉五月 | a级毛片在线看网站| 久久精品国产综合久久久| 精品人妻熟女毛片av久久网站| 国产精品一区二区在线不卡| 久久狼人影院| 十分钟在线观看高清视频www| www.自偷自拍.com| 69精品国产乱码久久久| 成年人午夜在线观看视频| 成在线人永久免费视频| 久久 成人 亚洲| 亚洲熟妇熟女久久| 亚洲在线自拍视频| 18禁国产床啪视频网站| 精品午夜福利视频在线观看一区| 99国产极品粉嫩在线观看| 国产成人av教育| 日韩免费av在线播放| 操出白浆在线播放| 精品电影一区二区在线| 亚洲熟女精品中文字幕| 妹子高潮喷水视频| 两性夫妻黄色片| 午夜两性在线视频| 大香蕉久久成人网| 日本黄色视频三级网站网址 | 精品免费久久久久久久清纯 | 在线永久观看黄色视频| 久久久国产成人精品二区 | 欧美激情 高清一区二区三区| 黄色a级毛片大全视频| 国产亚洲欧美在线一区二区| 精品卡一卡二卡四卡免费| 51午夜福利影视在线观看| av免费在线观看网站| 成人影院久久| 搡老熟女国产l中国老女人| 91av网站免费观看| 精品福利观看| 久久九九热精品免费| 91字幕亚洲| aaaaa片日本免费| 在线视频色国产色| 国产亚洲欧美精品永久| 亚洲精品久久成人aⅴ小说| av天堂久久9| 日韩三级视频一区二区三区| 男男h啪啪无遮挡| 99在线人妻在线中文字幕 | 午夜福利欧美成人| 欧美日韩亚洲国产一区二区在线观看 | av中文乱码字幕在线| 国产精品秋霞免费鲁丝片| 女人高潮潮喷娇喘18禁视频| 美女 人体艺术 gogo| 十八禁高潮呻吟视频| 成人国语在线视频| 丁香欧美五月| 亚洲第一欧美日韩一区二区三区| 久久国产乱子伦精品免费另类| 成年女人毛片免费观看观看9 | 超色免费av| 亚洲精品国产区一区二| av福利片在线| 日韩欧美国产一区二区入口| 夜夜夜夜夜久久久久| 91成人精品电影| 99久久精品国产亚洲精品| 精品亚洲成国产av| av网站在线播放免费| 亚洲成av片中文字幕在线观看| 老司机深夜福利视频在线观看| 成年人黄色毛片网站| 婷婷精品国产亚洲av在线 | 午夜福利在线免费观看网站| 在线观看免费视频日本深夜| 国产97色在线日韩免费| 99精品欧美一区二区三区四区| 亚洲一区二区三区不卡视频| 午夜精品在线福利| 精品一区二区三卡| 99国产精品99久久久久| 母亲3免费完整高清在线观看| 好男人电影高清在线观看| 高清在线国产一区| 国产精品av久久久久免费| 热re99久久精品国产66热6| 久久久国产成人精品二区 | 国产精品美女特级片免费视频播放器 | 欧美激情 高清一区二区三区| 亚洲欧洲精品一区二区精品久久久| 亚洲精品一二三| 最近最新中文字幕大全电影3 | 女人久久www免费人成看片| 老司机午夜福利在线观看视频| 亚洲美女黄片视频| 夫妻午夜视频| 成年人免费黄色播放视频| 一本综合久久免费| 日韩人妻精品一区2区三区| 中文亚洲av片在线观看爽 | 999精品在线视频| 久久这里只有精品19| 国产淫语在线视频| xxx96com| 视频区欧美日本亚洲| 亚洲va日本ⅴa欧美va伊人久久| 亚洲一区二区三区欧美精品| 妹子高潮喷水视频| 国产成人av教育| 久久国产精品男人的天堂亚洲| 久久久久精品国产欧美久久久| 久久精品国产清高在天天线| 久久久久久久久免费视频了| 国产精品.久久久| 欧美乱色亚洲激情| 高清毛片免费观看视频网站 | 757午夜福利合集在线观看| 狠狠狠狠99中文字幕| 亚洲中文字幕日韩| 精品熟女少妇八av免费久了| 91成人精品电影| 国产精品永久免费网站| 久久久久精品人妻al黑| 亚洲av成人一区二区三| 一夜夜www| 国产精品1区2区在线观看. | 黑人操中国人逼视频| 精品免费久久久久久久清纯 | 国产色视频综合| 自线自在国产av| 美女扒开内裤让男人捅视频| 美女高潮喷水抽搐中文字幕| 亚洲av成人一区二区三| 国产成人精品无人区| 国产成人欧美| 免费不卡黄色视频| 在线av久久热| 99re在线观看精品视频| 欧美久久黑人一区二区| 国内毛片毛片毛片毛片毛片| 下体分泌物呈黄色| 久99久视频精品免费| 亚洲国产看品久久| 50天的宝宝边吃奶边哭怎么回事|