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

    基于非線性混沌理論的非平穩(wěn)信號的比較分析*

    2014-09-17 06:01:38魏春雨楊威王梓卉敏韓清鵬
    動力學(xué)與控制學(xué)報 2014年2期
    關(guān)鍵詞:振動誤差預(yù)測

    魏春雨 楊威 王梓卉敏 韓清鵬

    (遼寧科技大學(xué)機(jī)械工程與自動化學(xué)院,鞍山 114051)

    基于非線性混沌理論的非平穩(wěn)信號的比較分析*

    魏春雨 楊威 王梓卉敏 韓清鵬?

    (遼寧科技大學(xué)機(jī)械工程與自動化學(xué)院,鞍山 114051)

    本文研究采用基于非線性混沌理論的兩種非線性參數(shù)估計方法(代替數(shù)據(jù)法和Lyapunov指數(shù)估計法)對非平穩(wěn)信號進(jìn)行分析.首先對上述兩種非線性方法的具體算法進(jìn)行介紹,然后對兩組本質(zhì)不同的非平穩(wěn)振動信號進(jìn)行對比分析.這兩組信號是通過測試具有不同非線性約束邊界條件的薄壁構(gòu)件獲得.分析結(jié)果表明,在時域波形上直觀相似的非平穩(wěn)信號,用上述非線性混沌分析的方法可以有效地加以定量區(qū)分.

    非線性混沌理論, 非平穩(wěn)信號, 代替數(shù)據(jù)法, Lyapunov指數(shù)

    引言

    人們通常考慮振動信號是否具有平穩(wěn)或非平穩(wěn)的性質(zhì),而線性或非線性特性則是對其源系統(tǒng)而言[1].在工程中,大多數(shù)系統(tǒng)所測得的信號是非平穩(wěn)的,而系統(tǒng)是非線性的.區(qū)分非平穩(wěn)信號在現(xiàn)代信號處理中,非平穩(wěn)信號的處理是最為活躍、發(fā)展最為迅速的方向之一[2],在通訊、雷達(dá)、自動控制、模式識別、機(jī)械振動和生物醫(yī)學(xué)工程等領(lǐng)域廣泛應(yīng)用.通過分析信號的時變特征,構(gòu)造合適的時頻分布并進(jìn)行恰當(dāng)?shù)奶幚?,達(dá)到不同的信號處理目的.近年來,人們采用多種時頻分析方法對非線性和非平穩(wěn)信號進(jìn)行分析,包括小波分析、Wagner-Ville分布等[3-6].此外,還有許多非線性分析的方法也加以應(yīng)用[7-9].但是,從非線性角度,基于混沌理論的一些方法,如Lyapunov指數(shù)和分?jǐn)?shù)維,對于具有非線性特性的非平穩(wěn)信號分析更為有效[10,11].另外,對信號進(jìn)行混沌辨識,并將非線性信號從隨機(jī)噪聲中區(qū)分出來往往是比較困難的.這是因為混沌信號和隨機(jī)噪聲信號通常都具有相似的寬頻特征.目前區(qū)分混沌和噪聲的常用方法主要有兩類:維數(shù)估計方法和非線性預(yù)測法.第一類方法是基于隨機(jī)只能存在無限維吸引子的原理,而有限維吸引子意味著混沌,如在 Grassberger-Procaccia算法中所指出的那樣[10].第二類方法則是根據(jù)時間序列的非線性預(yù)測原理、混沌具有不同于隨機(jī)數(shù)據(jù)的短期預(yù)測性的性質(zhì)來區(qū)分混沌和噪聲信號.也就是對于一個混沌時間序列,它的預(yù)測值和實際值之間的相關(guān)系數(shù)會隨預(yù)測時間的增加而減弱;而對于隨機(jī)時間序列,這種相關(guān)性不會因預(yù)測時間的增加而變化.

    本文針對由金屬薄壁構(gòu)件實測得到的具有不同非線性約束邊界條件、不同性質(zhì)的非平穩(wěn)振動信號,采用兩種不同的非線性預(yù)測方法進(jìn)行對比分析,即代替數(shù)據(jù)法和Lyapunov指數(shù)方法.其中代替數(shù)據(jù)法主要參考了文獻(xiàn)[12].這種方法可以通過對比混沌信號和噪聲信號(包括白噪聲和有色噪聲)的預(yù)測誤差分布及其所對應(yīng)的代替數(shù)據(jù)集的預(yù)測誤差分布的差異來實現(xiàn).Lyapunov指數(shù)是某過程在其相平面內(nèi)相鄰軌線的平均發(fā)散速率的量化定義.正的Lyapunov指數(shù)(一個或多個)是進(jìn)行混沌辨識的重要指標(biāo).時間序列的Lyapunov指數(shù)的估算方法有多種,見文獻(xiàn)[13-17].

    1 兩種非線性預(yù)測方法的算法原理

    1.1 基于代替數(shù)據(jù)法的混沌辨識

    (1)信號的代替數(shù)據(jù)集的生成

    信號所對應(yīng)的代替數(shù)據(jù)集可以由基于時間序列的 Gauss隨機(jī)過程假設(shè)得到.本文通過對非Gauss過程有效的非線性直方圖變換方法來實現(xiàn).對于一個原始信號,一般生成包含128個不同時間

    接下來,通過將上述復(fù)數(shù)乘以eiφ使其相位角在每個頻率上隨機(jī)化,其中φ是歸一化的、在[0,2π]區(qū)間內(nèi)變化的隨機(jī)量,得到新的X'(k).對其進(jìn)行逆傅立葉變換,可以得到Gauss型代替數(shù)據(jù)序列,也就是得到具有原始信號相同的幅值分布形式的代替數(shù)據(jù)序列x'(n),如下式序列的代替數(shù)據(jù)集.

    進(jìn)行直方圖變換時,首先生成一個與原給定時間序列長度相同的Gauss隨機(jī)數(shù)集合,然后對這個Gauss數(shù)據(jù)集的順序進(jìn)行重排.新的時間序列應(yīng)是具有Gauss概率密度函數(shù)分布的、與原序列相對應(yīng)的非線性尺度變換的結(jié)果.對于原始信號是Gauss隨機(jī)的情況,變換后的序列也具有Gauss分布.第一步是,將時間序列x(n)進(jìn)行傅立葉變換,得到

    由于上式的逆傅立葉變換是實數(shù),其相位角具有對稱性,即φ(k)=-φ(N-k).

    (2)統(tǒng)計量的計算

    由于低維混沌意味著其系統(tǒng)在短期內(nèi)可以視為是確定性的,而隨機(jī)過程與此不同.可以將預(yù)測誤差ε取為統(tǒng)計量.首先利用狀態(tài)變量x將時間序列x(n)進(jìn)行相空間重構(gòu),采用時間滯后的嵌入法,即

    其中嵌入維數(shù)d應(yīng)該滿足d≥2D+1,D是系統(tǒng)的真實吸引子維數(shù),在實際計算過程中往往需要依靠經(jīng)驗選?。⒁獾竭@里的n小于數(shù)據(jù)長度N,即1≤n≤N.時間滯后點τ的選擇有時具有一定隨意性,在這里簡單地取作1.

    下面將時間序列的數(shù)據(jù)點集分成長度相等(Nf=Nt)的擬合集和檢驗集兩部分.在擬合集中,尋找與當(dāng)前點歐幾里德距離最為相鄰的k個狀態(tài)點.這k個當(dāng)前時刻為m的狀態(tài)點與對應(yīng)的下一時刻m+1的狀態(tài)點組成如下k個點對:

    采用如(5)式所示的預(yù)測公式對這k個點對進(jìn)行擬合

    其中,a和b是擬合系數(shù),ˉx是x的預(yù)測點.

    接下來,計算檢驗集Nt中的所有點對應(yīng)的預(yù)測誤差,將得到N/2個預(yù)測誤差.預(yù)測誤差的定義為由(6)式得到的下一時刻(m+1)的預(yù)測值和實際值的差值.

    進(jìn)行混沌識別的統(tǒng)計量取為上述全部預(yù)測誤差的平均絕對誤差(MAE).

    (3)假設(shè)檢驗

    對于生成的所有代替數(shù)據(jù)集,根據(jù)如下公式(7)來計算顯著度:

    其中,QD是由原始信號時間序列計算得到的統(tǒng)計量值(MAE值),us和σs分別是由生成的128個代替數(shù)據(jù)序列計算得到的統(tǒng)計量值的均值和方差.

    計算得到的χ值可以用于分析原始信號數(shù)據(jù)和代替數(shù)據(jù)的差異.如果χ是一個較小的數(shù),這意味著原始信號和它的所有代替數(shù)據(jù)集具有相同的性質(zhì),因此隨機(jī)假設(shè)可以接受,也就是原始信號是隨機(jī)的.相反,如果χ值較大,可以認(rèn)為代替數(shù)據(jù)序列與原始信號有較大的差別,拒絕隨機(jī)假設(shè).

    更進(jìn)一步,為了辨識原始數(shù)據(jù)序列是隨機(jī)的還是混沌的,定義如式(8)所示的置信判據(jù).拒絕隨機(jī)假設(shè)的最大概率也就是相應(yīng)的顯著度P定義為

    其中erf()為數(shù)據(jù)序列所有元素的誤差函數(shù),其定義如式(9)所示.

    根據(jù)經(jīng)驗,如果計算得到的概率P值小于0.05,原始信號數(shù)據(jù)將顯著地不同于它的代替數(shù)據(jù)集,這時可以拒絕隨機(jī)假設(shè),認(rèn)為原始信號在95%置信度下是混沌的.如果P值大于0.05,則認(rèn)為原始信號是隨機(jī)的(95%置信度).當(dāng)然,這里的臨界值0.05可以根據(jù)實際情況的不同而不同.

    1.2 Lyapunov指數(shù)的估計算法

    Lyapunov指數(shù)的估計算法也是基于非線性預(yù)測理論的.進(jìn)行相空間重構(gòu)后,考慮兩條具有不同初始點的相鄰軌線L1和L2,它們的初始點分別是x0和z0.這兩個初始點之間的距離d0=|z0-x0|.經(jīng)過時間Δt后,x0和z0將沿著各自的軌線到達(dá)x1和y1,這時新的距離d1=|y1-x1|.在x1和y1之間選擇一個新點z1,并設(shè)d0=|z1-x1|.點x1和z1分別位于L1和L3之上.再過時間Δt,以x1和z1為起始點的軌線L1和L3將到達(dá)它們對應(yīng)的新點x2和y2.這樣經(jīng)過p次重復(fù),將得到di=|yi-xi|(i=1,2,…,p).Lyapunov指數(shù)的計算公式如式(10)所示.

    按數(shù)值大小進(jìn)行重排,得到該時間序列的Lyapunov指數(shù)如下

    因為p是一個較大的整數(shù),由此得到的Lyapunov指數(shù)是每個相點在其軌線上以指數(shù)形式發(fā)散的統(tǒng)計平均值.在上述Lyapunov指數(shù)中,一個或多個Lyapunov指數(shù)可能都是正的.根據(jù)非線性理論,正的Lyapunov指數(shù)意味著該時間序列(信號)是混沌的.

    2 實例

    對某薄壁構(gòu)件進(jìn)行寬頻隨機(jī)激振試驗.外激勵為500Hz范圍內(nèi)的標(biāo)準(zhǔn)白噪聲,對構(gòu)件的振動響應(yīng)進(jìn)行測試.由于是兩種性質(zhì)不同的非線性邊界條件,得到的構(gòu)件振動響應(yīng)呈現(xiàn)出明顯的非平穩(wěn)振動信號特征,其中應(yīng)該蘊(yùn)含著不同的特征.

    圖1所示為測試得到的兩組振動響應(yīng)時間信號,分別記為x1和x2.它們所對應(yīng)的功率譜密度(PSD)如圖2所示.從圖(2)只能看出,這兩個信號都具有0-500Hz有限頻帶分布的特征,形狀上僅有稍微區(qū)別.因此,僅從圖1和圖2人們很難區(qū)分這兩個信號.

    首先從非線性定性分析的角度對這兩組信號進(jìn)行比較.可以利用連續(xù)峰值來繪制這兩個信號的偽Poincare映射圖.設(shè)是一個時間序列的峰值集合,則它的偽Poincare映射圖是指對應(yīng)繪出的點圖.如果考慮采樣誤差,信號的周期點會在偽Poincare映射圖表現(xiàn)為一個較小的區(qū)域.而偽Poincare映射圖中出現(xiàn)分散的點區(qū)域時,根據(jù)非線性混沌理論,表明存在不規(guī)則或奇怪吸引子.對于這兩個信號x1和x2,它們的偽Poincare映射圖如圖3所示.從圖3(a)和(b)可以看出,信號x1和x2的偽Poincare映射圖在形貌上還是有區(qū)別的.

    圖1 非平穩(wěn)時間序列信號x1和x2Fig.1 Non - stationary time series of x1and x2

    圖2 x1和x2的功率譜密度Fig.2 Power spectral densities of x1and x2

    根據(jù)前文所介紹的非線性混沌分析理論,用代替數(shù)據(jù)法計算得到的關(guān)于x1和x2的概率P值見表1.如1.1節(jié)所言,如果P值小于0.05,則拒絕隨機(jī)假設(shè),信號的混沌特性得到辨識(95%的置信度).因此信號x1和x2都可以認(rèn)為是混沌時間序列.另外,從表1中的其它特征值來看,如QD,us和χ,也可以對這兩組信號加以定量地對比區(qū)分.

    圖3 x1和x2偽Poincare映射圖Fig.3 Pseudo Poincare mapping portraits of x1and x2

    計算得到的x1和x2的Lyapunov指數(shù)如表2所示.在估算這兩組信號的Lyapunov指數(shù)的過程中,相空間重構(gòu)的嵌入維數(shù)設(shè)為5,滯后點數(shù)設(shè)為7.從表2可以看出,x1信號有2個正的Lyapunov指數(shù),即 1.1785、0.1846,信號x2有 3 個正的 Lyapunov指數(shù),即 0.9674、0.1857 和 0.0184.所以,這兩組信號可以視為是混沌的,但具有不同的混沌階數(shù).

    表1 代替數(shù)據(jù)法得到的x1和x2的特征參數(shù)數(shù)值Table 1 Characteristic parameter values of x1and x2by surrogate data method

    表2 估算得到的x1和x2的Lyapunov指數(shù)Table 2 Lyapunov index of estimated x1and x2

    3 結(jié)論

    本文利用基于混沌理論的非線性預(yù)測方法對非平穩(wěn)信號進(jìn)行了比較分析.從定量分析研究的角度,所采用的代替數(shù)據(jù)法和Lyapunov指數(shù)估計方法對區(qū)分具有不同非線性性質(zhì)的非平穩(wěn)信號是有效的.

    對于兩個典型的非平穩(wěn)振動信號x1和x2,代替數(shù)據(jù)法所得到的特征參數(shù)的數(shù)值是互不相同的,估算得到的最大Lyapunov指數(shù)也是正的.根據(jù)代替數(shù)據(jù)法中的概率值的大小可以看出,x1比x2具有更明顯的混沌特征,但是x1有2個正的Lyapunov指數(shù),而x2有3個正的Lyapunov指數(shù).

    1 Kantz H,Schreiber T.Nonlinear time series analysis.Cambridge:Cambridge University Press,1999

    2 張賢達(dá),保錚.非平穩(wěn)信號分析與處理.北京:國防工業(yè)出版社,1998(Zhang X D,Bao Z.Nonstationary signal analysis and processing.Beijing:National Defence Industry Press,1998(in Chinese))

    3 He Zhengjia,Meng Qingfeng,Zhao Juyuan,et al.Timefrequency(scale)analysis and diagnosis for nonstationary dynamic signal of machinery.International Journal of Plant Engineering and Management,1996,1(1):40 ~47

    4 Meltzer G S,Ivanov Y Y.Fault detection in gear drives with non-stationary rotational speed- Part I:the time-frequency approach.Mechanical Systems and Signal Processing,2003,17(5):1033 ~1047

    5 Wang Wenyi,Harrap M J.Condition monitoring of rolling element bearings by using cone kernel time-frequency distribution.In:Processings of the SPIE Conference on Measurement Technology and Intelligent Instrument.Bellingham,WA,1993:290~298

    6 Katvonik V,Stankovic L.Instantaneous frequency estimation using the Wigner distribution with varying and datadriven widow length.IEEE Transactions on Signal Processing,1998,46(9):2351 ~2325

    7 Casdagli M.Nonlinear prediction of chaotic time series.Physica D,1989,35(3):335~356

    8 Engbert R.Testing for nonlinearity:the role of surrogate data.Chaos,Solitons&Fractals,2002,13(1):79 ~84

    9 Banbrook M,Ushaw G,McLaughlin S.Lyapunov exponents from a time series:a noise-robust extraction algorithm.Chaos,Solitons&Fractals,1996,7(7):973 ~976

    10 Theiler J,Eubank S,Longtin A,et al.Testing for nonlin-earity in time series:the method of surrogate data.Physica D:Nonlinear Phenomena,1992,58(1-4):77 ~94

    11 Yonemoto K,Yanagawa T.Estimating the Lyapunov exponent from chaotic time series with dynamic noise.Statistical Methodology,2007,4(4):461 ~480

    12 Kennel M,Isabelle S.Method to distinguish possible chaos from coloured noise and to determine embedding parameters.Physical Review A,1992,46(6):3111 ~3118

    13 Kantz H.A robust method to estimate the maximal Lyapunov exponent of a time series.Physics Letters A,1994,185(1):77~87

    14 Wolf A,Swift J V,Swinney H L,et al.Determining Lyapunov exponents from a time series.Physica D:Nonlinear Phenomena,1985,16(3):285 ~317

    15 Rosenstein M T,Collins J J,Luca C J.A practical method for calculating largest Lyapunov exponents from small data sets.Physica D:Nonlinear Phenomena,1993,65(1-2):117~134

    16 周良強(qiáng),陳芳啟,葛鵬.一類大展弦比機(jī)翼非線性模型的穩(wěn)定性與分岔.動力學(xué)與控制學(xué)報,2012,10(3):222 ~227(Zhou L Q,Chen Q F,Ge P.Stability and bifurcation for a nonlinear model of high-aspect-ratio wings.Journal of Dynamics and Control,2012,10(3):222~27(in Chinese))

    17 呂書鋒,胡宇達(dá).正交各向異性疊層板的非線性主共振分析.動力學(xué)與控制學(xué)報,2009,7(1):35~38(Lu S F,Hu Y D.Nonlinear principal resonance of orthotropic laminated plates.Journal of Dynamics and Control,2009,7(1):35~38(in Chinese))

    *The project supported by the National Natural Science Foundation of China(10972192)

    ? Corresponding author E-mail:han1011@163.com

    COMPARATIVE ANALYSES OF NON-STATIONARY SIGNALS BASED ON NONLINEAR CHAOTIC THEORIES*

    Wei Chunyu Yan WeiWang Zihuimin Han Qingpeng?
    (College of Mechanical Engineering,Liaoning Science and Technology University,Anshan114051,China)

    In the paper,two nonlinear estimation methods based on nonlinear chaotic theory,surrogate data method and Lyapunov exponents,are used to distinguish the difference of non-stationary signals.After brief introduction of the corresponding algorithms,two typical different non-stationary signals measured from a thin-plate structure with different nonlinear constraining boundaries are taken to compare by using the above two methods respectively.The obtained results demonstrate that the apparently similar signals are distinguished effectively in quantitative way with applying above nonlinear chaotic analyses.

    nonlinear chaotic theory, non-stationary signals, surrogate data method, Lyapunov exponents

    9 October 2012,

    1 July 2013.

    10.6052/1672-6553-2014-014

    2012-10-09 收到第 1 稿,2013-07-01 收到修改稿.

    *國家自然科學(xué)基金資助項目(10972192)

    E-mail:han1011@163.com

    猜你喜歡
    振動誤差預(yù)測
    振動的思考
    無可預(yù)測
    黃河之聲(2022年10期)2022-09-27 13:59:46
    選修2-2期中考試預(yù)測卷(A卷)
    選修2-2期中考試預(yù)測卷(B卷)
    角接觸球軸承接觸角誤差控制
    哈爾濱軸承(2020年2期)2020-11-06 09:22:26
    振動與頻率
    Beidou, le système de navigation par satellite compatible et interopérable
    壓力容器制造誤差探究
    中立型Emden-Fowler微分方程的振動性
    不必預(yù)測未來,只需把握現(xiàn)在
    午夜福利在线观看吧| 91精品三级在线观看| 成人永久免费在线观看视频 | 欧美日韩成人在线一区二区| 日本五十路高清| av天堂久久9| 亚洲av国产av综合av卡| 不卡一级毛片| 国产片内射在线| 久久久久久久精品吃奶| 成人国产一区最新在线观看| 搡老乐熟女国产| 麻豆av在线久日| 狠狠精品人妻久久久久久综合| 久久影院123| 十八禁网站免费在线| 黄色 视频免费看| 国产av国产精品国产| 女警被强在线播放| 女警被强在线播放| 精品亚洲成国产av| 丝袜美足系列| 日韩有码中文字幕| 最近最新中文字幕大全免费视频| 欧美变态另类bdsm刘玥| 又黄又粗又硬又大视频| 国产欧美日韩一区二区精品| 1024香蕉在线观看| 少妇粗大呻吟视频| 热re99久久国产66热| 日本一区二区免费在线视频| 视频区欧美日本亚洲| 亚洲美女黄片视频| 一本大道久久a久久精品| 操出白浆在线播放| 1024香蕉在线观看| 久久久久久久国产电影| 欧美久久黑人一区二区| 日本a在线网址| 国产不卡一卡二| 热re99久久国产66热| 成人国产av品久久久| 一个人免费在线观看的高清视频| 久久久久久人人人人人| 精品高清国产在线一区| 三级毛片av免费| 日韩欧美免费精品| h视频一区二区三区| 悠悠久久av| 精品国内亚洲2022精品成人 | 欧美日韩黄片免| 18禁黄网站禁片午夜丰满| 欧美黄色片欧美黄色片| 亚洲av片天天在线观看| 日韩一卡2卡3卡4卡2021年| 性色av乱码一区二区三区2| 大型黄色视频在线免费观看| 不卡av一区二区三区| 成人18禁在线播放| 日本av免费视频播放| 国内毛片毛片毛片毛片毛片| 免费少妇av软件| 精品一区二区三区视频在线观看免费 | 久久精品国产综合久久久| 女人精品久久久久毛片| cao死你这个sao货| 精品少妇黑人巨大在线播放| 99国产极品粉嫩在线观看| 午夜两性在线视频| 免费一级毛片在线播放高清视频 | 欧美黑人欧美精品刺激| 999久久久国产精品视频| 老熟妇乱子伦视频在线观看| 啦啦啦中文免费视频观看日本| 色婷婷av一区二区三区视频| 夜夜爽天天搞| 一进一出好大好爽视频| 亚洲欧美日韩另类电影网站| 99九九在线精品视频| 高清在线国产一区| 国产成人一区二区三区免费视频网站| 亚洲av电影在线进入| 亚洲视频免费观看视频| 国产成人欧美| 两性午夜刺激爽爽歪歪视频在线观看 | 老司机福利观看| 亚洲av欧美aⅴ国产| 精品人妻1区二区| 九色亚洲精品在线播放| 亚洲精品国产色婷婷电影| 久9热在线精品视频| 色在线成人网| 热re99久久国产66热| av一本久久久久| 一边摸一边抽搐一进一小说 | 少妇裸体淫交视频免费看高清 | 国产精品 国内视频| 日本欧美视频一区| 亚洲中文字幕日韩| 热re99久久国产66热| 狂野欧美激情性xxxx| 亚洲五月婷婷丁香| 国产免费现黄频在线看| kizo精华| 国产日韩欧美亚洲二区| 国产日韩欧美视频二区| 欧美日韩国产mv在线观看视频| 久久久国产一区二区| 最黄视频免费看| 老鸭窝网址在线观看| 老汉色av国产亚洲站长工具| 国产日韩欧美亚洲二区| 成人特级黄色片久久久久久久 | 久久久久网色| 精品国产国语对白av| www.熟女人妻精品国产| 亚洲人成伊人成综合网2020| 夜夜爽天天搞| 午夜福利免费观看在线| 精品亚洲成国产av| 欧美成人午夜精品| 久久精品国产综合久久久| 成人国产av品久久久| 别揉我奶头~嗯~啊~动态视频| 性高湖久久久久久久久免费观看| 人妻 亚洲 视频| www.自偷自拍.com| 国产精品1区2区在线观看. | 亚洲国产成人一精品久久久| 天天躁日日躁夜夜躁夜夜| 久久人妻av系列| 女人精品久久久久毛片| 欧美成狂野欧美在线观看| 亚洲熟妇熟女久久| 日韩 欧美 亚洲 中文字幕| 欧美 日韩 精品 国产| 在线十欧美十亚洲十日本专区| 亚洲精品国产色婷婷电影| 日韩大片免费观看网站| 最黄视频免费看| 色精品久久人妻99蜜桃| 亚洲中文av在线| 蜜桃在线观看..| 午夜福利视频在线观看免费| 欧美人与性动交α欧美软件| 国产精品一区二区免费欧美| 一本—道久久a久久精品蜜桃钙片| 亚洲一卡2卡3卡4卡5卡精品中文| 国产男女超爽视频在线观看| 午夜福利,免费看| 久久精品熟女亚洲av麻豆精品| 女人高潮潮喷娇喘18禁视频| 国产无遮挡羞羞视频在线观看| 亚洲欧美日韩另类电影网站| 国产精品免费视频内射| 美国免费a级毛片| 久久青草综合色| 亚洲精品自拍成人| 欧美精品人与动牲交sv欧美| 午夜成年电影在线免费观看| 欧美在线黄色| av国产精品久久久久影院| 国产成人一区二区三区免费视频网站| 桃红色精品国产亚洲av| av超薄肉色丝袜交足视频| 国产精品自产拍在线观看55亚洲 | 男人舔女人的私密视频| 精品久久久久久电影网| 69av精品久久久久久 | 日韩大片免费观看网站| 一区在线观看完整版| 国产一区二区在线观看av| 精品久久久久久电影网| 中文字幕人妻丝袜制服| 日韩熟女老妇一区二区性免费视频| 人妻 亚洲 视频| 在线观看免费视频网站a站| 一边摸一边抽搐一进一出视频| 熟女少妇亚洲综合色aaa.| 国产精品亚洲av一区麻豆| 99精品久久久久人妻精品| 中亚洲国语对白在线视频| 性色av乱码一区二区三区2| videosex国产| 久久久久久久精品吃奶| cao死你这个sao货| 满18在线观看网站| 丝袜美腿诱惑在线| 国产精品一区二区在线观看99| 亚洲,欧美精品.| 日韩欧美一区二区三区在线观看 | 69av精品久久久久久 | 国产色视频综合| 这个男人来自地球电影免费观看| 国产有黄有色有爽视频| 国产男女超爽视频在线观看| 天天添夜夜摸| 91成年电影在线观看| 99国产精品免费福利视频| 午夜精品国产一区二区电影| 真人做人爱边吃奶动态| 在线av久久热| 视频区图区小说| 女警被强在线播放| 国产熟女午夜一区二区三区| 嫁个100分男人电影在线观看| 男女下面插进去视频免费观看| 精品国产一区二区三区久久久樱花| 成年人免费黄色播放视频| 久久精品熟女亚洲av麻豆精品| 一个人免费看片子| 免费看a级黄色片| 亚洲av成人不卡在线观看播放网| 999精品在线视频| 亚洲av成人不卡在线观看播放网| 日本五十路高清| 日韩 欧美 亚洲 中文字幕| 日本欧美视频一区| 男人舔女人的私密视频| 亚洲成人免费av在线播放| tube8黄色片| 精品一区二区三卡| 91麻豆精品激情在线观看国产 | 免费在线观看完整版高清| kizo精华| 性高湖久久久久久久久免费观看| 岛国毛片在线播放| 成年人免费黄色播放视频| 亚洲av日韩精品久久久久久密| 国产欧美日韩一区二区三| 亚洲成av片中文字幕在线观看| 日韩三级视频一区二区三区| 国产一区二区激情短视频| 久久人妻福利社区极品人妻图片| 狠狠婷婷综合久久久久久88av| 久久99一区二区三区| 亚洲精华国产精华精| 国产精品成人在线| 国产欧美日韩综合在线一区二区| aaaaa片日本免费| 成人国产av品久久久| 黑人欧美特级aaaaaa片| 淫妇啪啪啪对白视频| 亚洲中文日韩欧美视频| 亚洲一区中文字幕在线| 国产视频一区二区在线看| 国产一区二区三区视频了| 免费看a级黄色片| 亚洲精品久久成人aⅴ小说| 99国产精品免费福利视频| 岛国毛片在线播放| 午夜福利在线免费观看网站| 天堂中文最新版在线下载| 91av网站免费观看| 亚洲精品自拍成人| 久久性视频一级片| 亚洲国产av影院在线观看| 久久久国产一区二区| 真人做人爱边吃奶动态| 国产精品欧美亚洲77777| 国产精品一区二区免费欧美| 亚洲av成人一区二区三| 香蕉久久夜色| 国产精品亚洲一级av第二区| 精品一区二区三卡| 久久中文字幕人妻熟女| 夜夜爽天天搞| 亚洲一卡2卡3卡4卡5卡精品中文| 一本—道久久a久久精品蜜桃钙片| 美女扒开内裤让男人捅视频| 91老司机精品| 欧美日韩国产mv在线观看视频| 亚洲成人手机| 亚洲欧美精品综合一区二区三区| 亚洲三区欧美一区| 91字幕亚洲| 亚洲成人免费av在线播放| 国产免费av片在线观看野外av| 国产精品.久久久| 狠狠婷婷综合久久久久久88av| 高清毛片免费观看视频网站 | 日韩有码中文字幕| 99热网站在线观看| 国产精品欧美亚洲77777| 在线观看免费高清a一片| 亚洲一区中文字幕在线| 中文亚洲av片在线观看爽 | 国产精品电影一区二区三区 | 欧美黄色淫秽网站| 法律面前人人平等表现在哪些方面| 中文字幕最新亚洲高清| 我的亚洲天堂| 夜夜骑夜夜射夜夜干| 免费日韩欧美在线观看| www日本在线高清视频| 亚洲中文字幕日韩| 伊人久久大香线蕉亚洲五| 天天躁狠狠躁夜夜躁狠狠躁| av在线播放免费不卡| 欧美变态另类bdsm刘玥| 捣出白浆h1v1| 精品国产乱码久久久久久男人| 超色免费av| 女人精品久久久久毛片| 在线播放国产精品三级| 青青草视频在线视频观看| av电影中文网址| 亚洲精品成人av观看孕妇| 久久精品国产a三级三级三级| 久久精品aⅴ一区二区三区四区| 大香蕉久久网| 亚洲av欧美aⅴ国产| 热99re8久久精品国产| 夫妻午夜视频| 久久国产亚洲av麻豆专区| 午夜激情av网站| 久久久久久人人人人人| 国产成人系列免费观看| 男女床上黄色一级片免费看| 人人妻人人澡人人爽人人夜夜| 欧美精品一区二区大全| 精品熟女少妇八av免费久了| 免费不卡黄色视频| 黄片播放在线免费| 国产成+人综合+亚洲专区| 国产亚洲av高清不卡| 免费看a级黄色片| 18禁国产床啪视频网站| 成人亚洲精品一区在线观看| 俄罗斯特黄特色一大片| 如日韩欧美国产精品一区二区三区| 久久九九热精品免费| 亚洲精品乱久久久久久| 亚洲成人手机| 在线观看免费视频日本深夜| 国产高清视频在线播放一区| 亚洲精品在线观看二区| 国产真人三级小视频在线观看| 欧美老熟妇乱子伦牲交| 国产精品二区激情视频| 欧美激情高清一区二区三区| 国产成人影院久久av| 久久久精品国产亚洲av高清涩受| 狂野欧美激情性xxxx| 国产日韩一区二区三区精品不卡| 99久久99久久久精品蜜桃| 亚洲欧美激情在线| 18禁观看日本| 国产淫语在线视频| 女性生殖器流出的白浆| 又大又爽又粗| 亚洲第一欧美日韩一区二区三区 | 国产成+人综合+亚洲专区| 黑人巨大精品欧美一区二区mp4| 黑人巨大精品欧美一区二区蜜桃| 少妇被粗大的猛进出69影院| 国内毛片毛片毛片毛片毛片| 中文欧美无线码| 国产男女内射视频| 国产1区2区3区精品| 少妇 在线观看| 女人高潮潮喷娇喘18禁视频| 日本精品一区二区三区蜜桃| 午夜精品国产一区二区电影| 在线观看免费视频日本深夜| 97在线人人人人妻| 久久中文看片网| 亚洲人成77777在线视频| 国产成人影院久久av| 午夜免费鲁丝| 久久精品国产a三级三级三级| 精品一品国产午夜福利视频| 国产男女超爽视频在线观看| 精品福利永久在线观看| 两人在一起打扑克的视频| 日日夜夜操网爽| 亚洲欧美激情在线| 无遮挡黄片免费观看| 国产一区二区三区综合在线观看| 天堂中文最新版在线下载| 如日韩欧美国产精品一区二区三区| 色综合欧美亚洲国产小说| 亚洲午夜理论影院| 老汉色∧v一级毛片| 国产免费av片在线观看野外av| 99国产综合亚洲精品| 欧美日韩福利视频一区二区| 伦理电影免费视频| 两性午夜刺激爽爽歪歪视频在线观看 | 久久狼人影院| 成人亚洲精品一区在线观看| 视频在线观看一区二区三区| 少妇的丰满在线观看| 精品高清国产在线一区| 欧美日韩黄片免| 日本五十路高清| 国产淫语在线视频| 亚洲伊人色综图| 亚洲精品国产色婷婷电影| 丝袜喷水一区| 岛国在线观看网站| 一区二区三区乱码不卡18| 一本综合久久免费| 黄色视频在线播放观看不卡| 亚洲av成人不卡在线观看播放网| 亚洲av美国av| 久久久久国内视频| 久久久精品94久久精品| 不卡av一区二区三区| 午夜福利视频在线观看免费| av有码第一页| av福利片在线| 日日爽夜夜爽网站| 国产无遮挡羞羞视频在线观看| 久久天堂一区二区三区四区| 亚洲av电影在线进入| 国产成人av教育| 丁香六月欧美| 久久久久久久大尺度免费视频| 日本黄色日本黄色录像| 99久久国产精品久久久| 亚洲欧美日韩高清在线视频 | 久久精品人人爽人人爽视色| 18禁黄网站禁片午夜丰满| 亚洲精品av麻豆狂野| 亚洲国产欧美日韩在线播放| 一本一本久久a久久精品综合妖精| 汤姆久久久久久久影院中文字幕| 高清av免费在线| 欧美人与性动交α欧美软件| 丝袜喷水一区| 男女之事视频高清在线观看| 大香蕉久久成人网| 亚洲第一欧美日韩一区二区三区 | 久久久久网色| 亚洲第一av免费看| 69av精品久久久久久 | 久久久久久人人人人人| 国产欧美日韩精品亚洲av| 国产欧美日韩综合在线一区二区| 亚洲午夜理论影院| 少妇精品久久久久久久| 久久中文字幕人妻熟女| 色综合欧美亚洲国产小说| 黄色成人免费大全| 99re6热这里在线精品视频| 欧美日韩福利视频一区二区| 99精品在免费线老司机午夜| 日韩制服丝袜自拍偷拍| 国产福利在线免费观看视频| 脱女人内裤的视频| av不卡在线播放| 三上悠亚av全集在线观看| 国产精品二区激情视频| 中文字幕人妻熟女乱码| 黄片播放在线免费| 在线观看www视频免费| 国产成人欧美在线观看 | 国产精品av久久久久免费| 一进一出抽搐动态| 久久久久精品国产欧美久久久| 国产精品一区二区精品视频观看| 午夜福利乱码中文字幕| 一进一出抽搐动态| 香蕉久久夜色| 男人舔女人的私密视频| 19禁男女啪啪无遮挡网站| 午夜老司机福利片| 亚洲av成人一区二区三| 欧美 日韩 精品 国产| 国产精品免费一区二区三区在线 | 国产精品久久久久久精品古装| 99在线人妻在线中文字幕 | 亚洲国产精品一区二区三区在线| 菩萨蛮人人尽说江南好唐韦庄| 一本大道久久a久久精品| 国产不卡一卡二| 99精国产麻豆久久婷婷| 成人三级做爰电影| 亚洲专区国产一区二区| 深夜精品福利| 99热网站在线观看| 欧美日韩精品网址| 国产不卡av网站在线观看| 夜夜夜夜夜久久久久| 国产熟女午夜一区二区三区| 久久久久精品国产欧美久久久| 黄网站色视频无遮挡免费观看| 亚洲精品美女久久av网站| 女人精品久久久久毛片| 欧美性长视频在线观看| 久久久欧美国产精品| 免费观看a级毛片全部| 咕卡用的链子| 亚洲 国产 在线| 国产有黄有色有爽视频| 国产成人免费观看mmmm| 色老头精品视频在线观看| 国产av一区二区精品久久| 一区二区三区精品91| 精品国产国语对白av| 久久精品国产亚洲av香蕉五月 | 性色av乱码一区二区三区2| av片东京热男人的天堂| 久热爱精品视频在线9| 99在线人妻在线中文字幕 | 1024视频免费在线观看| 国产av国产精品国产| 涩涩av久久男人的天堂| 女性生殖器流出的白浆| 国产精品久久电影中文字幕 | 人人妻人人澡人人看| 国产免费视频播放在线视频| 波多野结衣av一区二区av| 汤姆久久久久久久影院中文字幕| 亚洲av成人一区二区三| 免费少妇av软件| 国产成人影院久久av| 亚洲色图综合在线观看| 在线观看一区二区三区激情| 一边摸一边抽搐一进一小说 | 在线av久久热| 国产精品98久久久久久宅男小说| 亚洲,欧美精品.| 十分钟在线观看高清视频www| 亚洲专区国产一区二区| 久久久久精品人妻al黑| √禁漫天堂资源中文www| 999精品在线视频| 狠狠狠狠99中文字幕| 国产在线精品亚洲第一网站| 97在线人人人人妻| 视频在线观看一区二区三区| 午夜免费鲁丝| 老熟妇仑乱视频hdxx| 久久久久久久精品吃奶| 亚洲欧美日韩高清在线视频 | 国产熟女午夜一区二区三区| 精品久久久精品久久久| 国产av一区二区精品久久| 亚洲人成电影免费在线| 下体分泌物呈黄色| 久久精品人人爽人人爽视色| 国产精品免费大片| 精品亚洲成a人片在线观看| 超色免费av| 久久精品亚洲熟妇少妇任你| 国产精品二区激情视频| 多毛熟女@视频| 亚洲伊人色综图| 欧美成狂野欧美在线观看| 热99国产精品久久久久久7| 日本精品一区二区三区蜜桃| 亚洲第一欧美日韩一区二区三区 | 国产淫语在线视频| 欧美日韩福利视频一区二区| av电影中文网址| 国产精品二区激情视频| 亚洲国产中文字幕在线视频| 亚洲国产成人一精品久久久| 中亚洲国语对白在线视频| 天天躁日日躁夜夜躁夜夜| 深夜精品福利| av天堂久久9| 国产在线观看jvid| 黄色片一级片一级黄色片| 国产亚洲av高清不卡| 极品教师在线免费播放| 看免费av毛片| 一进一出好大好爽视频| 日韩中文字幕视频在线看片| 久久人妻福利社区极品人妻图片| 亚洲 欧美一区二区三区| kizo精华| 欧美中文综合在线视频| 国产一区有黄有色的免费视频| 日本一区二区免费在线视频| 9热在线视频观看99| 人人妻人人澡人人爽人人夜夜| 高清在线国产一区| 99久久人妻综合| 亚洲熟女精品中文字幕| 超碰成人久久| 亚洲精品中文字幕一二三四区 | 一区福利在线观看| 欧美日韩亚洲国产一区二区在线观看 | 国产在线免费精品| 黑人猛操日本美女一级片| www.自偷自拍.com| 欧美亚洲 丝袜 人妻 在线| 国产精品99久久99久久久不卡| 成人免费观看视频高清| 午夜久久久在线观看| 波多野结衣av一区二区av| 99久久国产精品久久久| 飞空精品影院首页| 在线观看免费日韩欧美大片| 黄色成人免费大全| 纵有疾风起免费观看全集完整版| 91九色精品人成在线观看| 国产欧美日韩精品亚洲av| 国产成人av激情在线播放| 欧美午夜高清在线| 黑人猛操日本美女一级片| 精品一品国产午夜福利视频| 可以免费在线观看a视频的电影网站| 国产精品影院久久| 视频区图区小说| 亚洲av成人不卡在线观看播放网| 一区二区三区精品91| 一级黄色大片毛片| 中文字幕人妻丝袜制服| 一区二区三区乱码不卡18| 老司机影院毛片| 精品少妇内射三级|