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

    基于LPZ譜分析法的非高斯脈動(dòng)風(fēng)壓模擬

    2017-09-19 09:53:18蔣磊李春祥鄧瑩
    關(guān)鍵詞:偏度峰度譜分析

    蔣磊,李春祥,鄧瑩

    (上海大學(xué)土木工程系,上海200444)

    基于LPZ譜分析法的非高斯脈動(dòng)風(fēng)壓模擬

    蔣磊,李春祥,鄧瑩

    (上海大學(xué)土木工程系,上海200444)

    提出一種基于線性預(yù)測和Z轉(zhuǎn)換(linear prediction and Z-transform,LPZ)譜分析法的非高斯脈動(dòng)風(fēng)壓模擬算法.運(yùn)用Johnson變換系統(tǒng)實(shí)現(xiàn)高斯隨機(jī)過程到非高斯白噪聲的轉(zhuǎn)換;再使用LPZ對(duì)其進(jìn)行數(shù)字濾波,進(jìn)而得到所需的非高斯脈動(dòng)風(fēng)壓.采用基于LPZ譜分析法對(duì)單變量非高斯隨機(jī)信號(hào)和非高斯脈動(dòng)風(fēng)壓進(jìn)行了數(shù)值模擬.通過對(duì)模擬的非高斯信號(hào)和脈動(dòng)風(fēng)壓的統(tǒng)計(jì)參數(shù)(峰度和偏度)與目標(biāo)統(tǒng)計(jì)參數(shù),以及模擬的功率譜與目標(biāo)功率譜進(jìn)行比較,驗(yàn)證了基于LPZ譜分析法的非高斯脈動(dòng)風(fēng)壓模擬算法的有效性.

    非高斯脈動(dòng)風(fēng)壓模擬;線性預(yù)測和Z轉(zhuǎn)換譜分析法;Johnson變換系統(tǒng);數(shù)字濾波

    非高斯過程的模擬方法可以分為兩類:第一類是根據(jù)任意指定的高階特征統(tǒng)計(jì)參數(shù)(例如斜度與峰度)及其目標(biāo)功率譜密度(power spectral density,PSD)函數(shù)來實(shí)現(xiàn);第二類是通過任意指定的邊緣概率分布函數(shù)(probability density function,PDF)及其目標(biāo)功率譜密度(PSD)函數(shù)來實(shí)現(xiàn).Gurley等[1]、Seong等[2]和Suresh等[3]等在20世紀(jì)90年代就開始了第一類非高斯過程的模擬研究.Seong等[4]采用指數(shù)峰值模型進(jìn)行了單變量和多變量非高斯風(fēng)壓時(shí)程的模擬.Suresh等[3,5]基于快速傅里葉變換(fast Fourier transform,FFT)技術(shù),采用參數(shù)較少的指數(shù)峰值模型模擬了一維單變量非高斯風(fēng)壓時(shí)程,并用于大跨低矮屋蓋的風(fēng)振分析.第一類模擬方法過于繁瑣,建立的轉(zhuǎn)換關(guān)系往往為非線性,且在建立過程中并沒有考慮PSD,因而轉(zhuǎn)換后的PSD可能出現(xiàn)負(fù)值.對(duì)于第二類非高斯模擬,主要有KLE(Karhunen-Loeve)法[6]和Hermite譜修正法[7-9]等.Yang等[9]通過Hermite PDF模型模擬出了多變量非高斯風(fēng)壓.Ma等[10]提出了基于增強(qiáng)現(xiàn)實(shí)(augmented reality,AR)模型與Johnson逆變換的風(fēng)壓模擬方法.第二類模擬通常需要采用迭代算法.然而,隨著迭代次數(shù)的增加,模擬數(shù)值不一定會(huì)逐漸逼近目標(biāo)值,且當(dāng)?shù)揭欢ù螖?shù)時(shí),反而可能出現(xiàn)發(fā)散現(xiàn)象.針對(duì)第一類模擬算法,本工作提出一種基于線性預(yù)測和Z轉(zhuǎn)換(linear prediction and Z-transform,LPZ)譜分析法的非高斯脈動(dòng)風(fēng)壓非迭代模擬算法,簡化了已有算法,提高了運(yùn)算效率與模擬精度.

    1 LPZ譜分析法

    根據(jù)隨機(jī)過程理論,表征隨機(jī)變量數(shù)學(xué)特征的統(tǒng)計(jì)參數(shù)如下:均值(一階)μ=w=表∫示隨機(jī)過程的擺動(dòng)中心,其中p(w)表示隨機(jī)過程w的概率密度函數(shù);方差(二階)表示隨機(jī)過程在某時(shí)刻對(duì)于其均值的偏離程度;偏度(三階)是定量描述獨(dú)立隨機(jī)過程與高斯分布偏離程度的基本量;峰表示分布曲線尖峭或平坦的程度.

    高斯過程的偏度為0而峰度為3.偏度不為0或峰度不為3的隨機(jī)過程即為非高斯過程.偏度反映了以其平均值為中心的分布的不對(duì)稱程度,正偏度表示不對(duì)稱部分的分布更趨向正值,負(fù)偏度表示不對(duì)稱部分的分布更趨向負(fù)值.峰度能夠表明時(shí)域信號(hào)峰值的概率大小,描述了概率密度函數(shù)中間部分的尖銳度或平滑度以及尾部的窄度或?qū)挾?當(dāng)峰度K>3時(shí),表示分布呈高峰度,具有相對(duì)多的高峰值和低谷值,此時(shí)的非高斯過程稱為超高斯隨機(jī)過程(softening non-Gaussian processes).當(dāng)K<3時(shí),表示分布呈低峰度,具有相對(duì)較少的高峰值和低谷值,此時(shí)的非高斯過程稱為亞高斯隨機(jī)過程(hardening non-Gaussian processes).

    把離散的隨機(jī)脈動(dòng)風(fēng)壓序列{wr|r=0,1,2,···,N-1}作為研究的隨機(jī)過程,則表征該隨機(jī)脈動(dòng)風(fēng)壓的數(shù)學(xué)統(tǒng)計(jì)參數(shù)分別為

    LPZ譜分析法是FFT的一種替代和改進(jìn)[11].對(duì)于一個(gè)有限時(shí)間內(nèi)無限數(shù)據(jù)序列的離散信號(hào),運(yùn)用LPZ譜分析法可以避免因FFT的局限性和數(shù)據(jù)截?cái)嗨斐傻娜毕莺驼`差.盡管FFT被普遍認(rèn)為是信號(hào)從時(shí)域轉(zhuǎn)換到頻域的最終機(jī)制,但LPZ譜分析法具有更為優(yōu)良的頻率分辨率和靈敏度.此外,LPZ譜分析法可降噪且不損失分辨率,這是FFT算法無法達(dá)到的.正因?yàn)閮?yōu)良靈敏度和高分辨率,LPZ譜分析法替代FFT可獲得一些不容易獲得的頻譜信息,提高了模擬計(jì)算的準(zhǔn)確性.

    線性預(yù)測(linear prediction,LP)理論表明,每個(gè)復(fù)值數(shù)據(jù)點(diǎn)yn可由如下數(shù)值點(diǎn)yn+m的線性組合來求得其近似結(jié)果:

    式中,am為未知線性系數(shù),M為濾波長度,N為采樣點(diǎn)總數(shù).

    由式(1)可以推導(dǎo)出無限的數(shù)據(jù)序列以及由Z轉(zhuǎn)換定義的相應(yīng)譜函數(shù)S(z):

    式中,z=exp(-iωτ),1/τ為采樣率,ω/2π為頻率值.式(2)可以通過式(1)進(jìn)行簡化,得到LPZ譜分析法公式:

    同時(shí),式(1)的PSD函數(shù)可表示為

    2 Johnson變換系統(tǒng)

    Johnson變換系統(tǒng)可將高斯過程轉(zhuǎn)變?yōu)閾碛胁煌群头宥鹊姆歉咚剐蛄?基于統(tǒng)計(jì)參數(shù)算法,Johnson變換系統(tǒng)產(chǎn)生的頻率曲線可用于計(jì)算已知四階統(tǒng)計(jì)參數(shù)的隨機(jī)分布.在這一系統(tǒng)中,主要有3種不同形式的頻率曲線:無界系統(tǒng)(unbound system),即有界系統(tǒng)(bounded system),即為高斯隨機(jī)序列,η′′為給定偏度與峰度的非高斯隨機(jī)序列,γ,δ,ξ和λ為常數(shù),可以通過擬合前四階參數(shù)得到.Hill等[12]給出了求解這些常數(shù)的公式.

    3 基于LPZ譜分析法的非高斯脈動(dòng)風(fēng)壓模擬新算法

    本工作采用LPZ數(shù)字濾波技術(shù)[13]和Johnson變換系統(tǒng),建立基于LPZ譜分析法的非高斯脈動(dòng)風(fēng)壓模擬新算法.本算法屬于第一類模擬,即根據(jù)任意指定的高階特征統(tǒng)計(jì)參數(shù)(如斜度與峰度)及其目標(biāo)功率譜密度(PSD)函數(shù)來實(shí)現(xiàn)非高斯隨機(jī)過程的模擬.

    首先假設(shè){ηr}為一組高斯隨機(jī)序列,則這組高斯序列可通過有限脈沖響應(yīng)(finite impulse response,FIR)濾波器表示為

    對(duì)式(5)的兩邊分別進(jìn)行LPZ轉(zhuǎn)換,得到

    式中,Ak為高斯序列{ηr}的LPZ轉(zhuǎn)換形式,Wk為wr的LPZ轉(zhuǎn)換形式,Hk為系統(tǒng)的轉(zhuǎn)換方程(或頻率特性),具體等式關(guān)系為

    其中z=exp(-iωτ),r,k=0,1,2,···,N-1.Hk的逆LPZ轉(zhuǎn)換形式,即為濾波因子h:

    這二者之間的PSD函數(shù)存在如下轉(zhuǎn)化關(guān)系:

    式中,Sk為所求的非高斯隨機(jī)序列PSD函數(shù),而Sη為非高斯白噪聲隨機(jī)序列的PSD函數(shù),為常數(shù).隨機(jī)生成的高斯序列{ηr},通過Johnson變換系統(tǒng)轉(zhuǎn)化成擁有一定偏度SSkη和峰度Kη的非高斯白噪聲輸入序列.最后輸出序列中的偏度SSk和峰度K與輸入序列中SSkη和Kη的關(guān)系為

    綜上所述,為模擬指定偏度SSk、峰度K和PSD函數(shù)SSk的非高斯脈動(dòng)風(fēng)壓,整個(gè)模擬過程如下.

    步驟1已知Sk,并令Sη=C(C為常數(shù)),使用式(9)求解濾波方程

    步驟2使用式(8)求解濾波系數(shù)hr.

    步驟3已知非高斯序列的峰度K和偏度SSk,使用式(10)和(11)求解原輸入序列,即求解非高斯白噪聲風(fēng)壓時(shí)間序列η的偏度SSkη和峰度Kη.

    步驟4通過Johnson變換系統(tǒng),將隨機(jī)生成的高斯序列{ηr}轉(zhuǎn)換為具有指定偏度SSkη和峰度Kη的非高斯白噪聲隨機(jī)序列η.

    步驟5已知η和hr,通過方程(5)對(duì)步驟4中求出的非高斯白噪聲隨機(jī)風(fēng)壓η進(jìn)行降噪處理,求出擁有指定偏度SSk、峰度K和PSD函數(shù)SSk的非高斯隨機(jī)風(fēng)壓.

    4 數(shù)值驗(yàn)證

    4.1 白噪聲非高斯信號(hào)模擬

    首先,對(duì)一維單變量白噪聲非高斯信號(hào)進(jìn)行模擬,以驗(yàn)證上述方法的有效性.選圖1中的目標(biāo)PSD曲線為所需模擬的非高斯信號(hào)目標(biāo)功率譜密度函數(shù)Sk,該頻譜為200~500 Hz上的平直譜,均方根值σ約為5.96 g.由于fu=2 000 Hz,根據(jù)采樣定理確定fs= 5 120 Hz,N=4 096,Δf=1.25 Hz.模擬的目標(biāo)PSD函數(shù)為Sk,目標(biāo)峰度K=7.0,偏度SSk=-1.10.假設(shè)Sη=1,由式(10)得Hk=使用式(11)和(12),求得非高斯白噪聲隨機(jī)序列η的指定偏度SSkη=-2.33,峰度Kη=17.53.

    圖1 濾波輸出后的非高斯白噪聲隨機(jī)序列PSD函數(shù)F ig.1 PSD function of the fi ltered non-Gaussian white noise randomsequences

    通過Johnson變換系統(tǒng),使隨機(jī)生成的高斯序列{ηr}轉(zhuǎn)換為具有指定偏度SSkη和峰度Kη的非高斯白噪聲隨機(jī)序列η,如圖2所示.圖3為非高斯白噪聲隨機(jī)序列η的PSD函數(shù).

    圖2 Johnson變換系統(tǒng)產(chǎn)生的非高斯白噪聲隨機(jī)序列F ig.2 Non-Gaussian white noise randomsequences generated by Johnson translator system

    經(jīng)過FIR濾波器的降噪過濾,輸出的非高斯脈動(dòng)風(fēng)壓如圖4所示,其偏度和峰度分別為-1.06與7.08,與目標(biāo)偏度SSk=-1.10和峰度K=7.0相差很小.圖1為濾波輸出后的非高斯脈動(dòng)風(fēng)壓的PSD函數(shù)與目標(biāo)PSD函數(shù)的對(duì)比,發(fā)現(xiàn)二者吻合很好.值得指出的是,該信號(hào)模擬過程中的功率譜屬于分貝頻率坐標(biāo)設(shè)置下的形式,與線性頻率坐標(biāo)下的PSD不為負(fù)數(shù)并不沖突.因此,本工作中非高斯脈動(dòng)風(fēng)壓模擬將使用線性坐標(biāo).

    圖3 Johnson變換系統(tǒng)產(chǎn)生的非高斯白噪聲隨機(jī)序列PSD函數(shù)F ig.3 PSD function of the non-Gaussian white noise randomsequences generated by Johnson translator system

    圖4 濾波輸出后的非高斯白噪聲隨機(jī)序列Fig.4 Filtered non-Gaussian white noise randomsequences

    4.2 非高斯脈動(dòng)風(fēng)壓模擬與算法對(duì)比

    下面進(jìn)行一維單變量非高斯脈動(dòng)風(fēng)壓的模擬.使用李錦華等[14]推導(dǎo)出的脈動(dòng)風(fēng)壓功率譜方程模擬符合目標(biāo)脈動(dòng)風(fēng)壓功率譜、偏度和峰度的非高斯風(fēng)壓.Davenport水平脈動(dòng)風(fēng)速功率譜為

    表1 相關(guān)模擬參數(shù)Tab le 1 Related simu lation parameters

    圖5為通過Johnson變換系統(tǒng)轉(zhuǎn)換得到的非高斯白噪聲風(fēng)壓.圖6為非高斯白噪聲隨機(jī)風(fēng)壓的PSD函數(shù).經(jīng)過FIR濾波器的降噪過濾,輸出的非高斯脈動(dòng)風(fēng)壓如圖7所示,其偏度和峰度分別為-2.49與4.77,與目標(biāo)偏度SSk=-2.4和峰度K=4.8非常接近.圖8為濾波輸出后的非高斯脈動(dòng)風(fēng)壓的PSD函數(shù),與目標(biāo)PSD函數(shù)吻合良好.

    圖5 Johnson變換系統(tǒng)產(chǎn)生的非高斯白噪聲脈動(dòng)風(fēng)壓Fig.5 Non-Gaussian white noise fluctuating w ind pressure generated by Johnson translator system

    圖6 Johnson變換系統(tǒng)產(chǎn)生的非高斯白噪聲脈動(dòng)風(fēng)壓PSD函數(shù)Fig.6 PSD function of the non-Gaussian white noise fluctuating w ind pressure generated by Johnson translator system

    5 結(jié)束語

    本工作提出了一種基于LPZ譜分析法的非高斯脈動(dòng)風(fēng)壓模擬新算法,即運(yùn)用Johnson變換系統(tǒng)實(shí)行高斯隨機(jī)過程到非高斯白噪聲的轉(zhuǎn)換,再使用LPZ譜分析實(shí)現(xiàn)了非高斯白噪聲譜到目標(biāo)譜的轉(zhuǎn)化.模擬譜與目標(biāo)譜基本一致,偏度、峰度值與目標(biāo)值近似,故而所模擬的非高斯脈動(dòng)風(fēng)壓與目標(biāo)風(fēng)壓趨于一致.

    圖7 濾波輸出后的非高斯白噪聲脈動(dòng)風(fēng)壓Fig.7 Filtered non-Gaussian white noise fluctuating w ind pressure

    圖8 模擬非高斯白噪聲脈動(dòng)風(fēng)壓PSD函數(shù)與目標(biāo)PSD函數(shù)的比較Fig.8 Comparisons of the PSD function of the simulated and the target non-Gaussian white noise fluctuating w ind pressure

    通過MATLAB編程,使用本算法分別對(duì)非高斯隨機(jī)信號(hào)和非高斯脈動(dòng)風(fēng)壓進(jìn)行了模擬.通過比較模擬非高斯脈動(dòng)風(fēng)壓的統(tǒng)計(jì)參數(shù)(峰度與偏度)與目標(biāo)統(tǒng)計(jì)參數(shù),以及模擬功率譜與目標(biāo)功率譜,驗(yàn)證了本算法的有效性.

    [1]G URLEY K,K AREEMA.Analysis interpretation modeling and simu lation of unsteady w ind and pressure data[J].JournalofW ind Engineering and Industrial Aerodynamics,1997,69:657-669. [2]S EONG S H,P ETERKAJ A.Computer simulation of non-Gaussian mu ltiple w ind press time series[J].Journal ofW ind Engineering and Industrial Aerodynamics,1997,72(1):95-105.

    [3]S URESH K K,S TATHOpOULOS T.Synthesis of non-Gaussian w ind pressure time series on low building roofs[J].Engineering Structures,1999,21(12):1086-1100.

    [4]S EONG S H,P ETERKAJ A.D igital generation of surface-pressure fluctuations with spiky features[J].Journal of W ind Engineering and Industrial Aerodynamics,1998,73(2):181-192.

    [5]S URESH K K,S TATHOpOULOS T.Computer simulation of fluctuating w ind pressures on low building roofs[J].Journal ofW ind Engineering and Industrial Aerodynamics,1997,69(1):485-495.

    [6]V OˇRCHOVSK′Y M.Simulation of simply cross correlated randomfields by series expansion methods[J].Structural Safety,2008,30(4):337-373.

    [7]Y ANG Q,T IAN Y.Amodel of probability density function of non-Gaussian w ind pressure with multiple samples[J].Journal of W ind Engineering and Industrial Aerodynamics,2015,140: 67-78.

    [8]H UANG G,L UO Y,G URLEY K R,et al.Revisiting moment-based characterization for w ind pressures[J].Journal ofW ind Engineering and Industrial Aerodynamics,2016,151:158-168.

    [9]Y ANG L,G URLEY K R.Effi cient stationary mu ltivariate non-Gaussian simulation based on a Hermite PDF model[J].Probabilistic Engineering Mechanics,2015,42(4):31-41.

    [10]MAX L,X U F Y,K AREEMA,et al.Estimation of surface pressure extremes:hybrid data and simulation-based approach[J].Journal of Engineering Mechanics,2016,DOI: 10.1061/(ASCE)EM.1943-7889.0001127.

    [11]T ANG J,N ORRIS J R.LPZ spectral analysis using linear prediction and the Z transform[J]. Journal of Chemical Physics,1986,84(9):5210-5211.

    [12]H ILL ID,H ILL R,H OLDER R L.Fitting Johnson curves by moments[J].Applied Statistics, 1976,25:149-176.

    [13]H U Y Z,T ONDER K.Simulation of 3-D randomrough surface by 2-D digital fi lter and Fourier analysis[J].International Journal of Machine Tools Manufacturing,1992,32(1/2):83-90.

    [14]李錦華,李春祥.土木工程隨機(jī)風(fēng)場數(shù)值模擬研究的進(jìn)展[J].振動(dòng)與沖擊,2008,19(9):116-125, 187.

    Simu lation of non-Gaussian fl uctuating w ind pressu re based on LPZ spectral analysis

    JIANG Lei,LIChunxiang,DENG Ying
    (Department of Civil Engineering,Shanghai University,Shanghai 200444,China)

    This paper presents an algorithmbased on linear prediction and Z-transform(LPZ)spectralanalysis to simulate non-Gaussian fluctuating w ind pressure.The Gaussian process is converted into a non-Gaussian white noise process by Johnson translator system. The non-Gaussian white noise process is then fi ltered with LPZ spectral analysis to obtain fluctuating w ind pressure.The univariate non-Gaussian randomsignals and non-Gaussian fluctuating w ind pressure are simulated using the algorithm.By comparing the statistical parameters including skewness,kurtosis and power spectral density of the non-Gaussian white noise process and fl uctuating w ind pressure with their target values,it is confi rmed that the algorithmbased on LZP spectral analysis can eff ectively simulate non-Gaussian fluctuating w ind pressure.

    non-Gaussian fl uctuating w ind pressure simu lation;linear prediction and Z-transform(LZP)spectral analysis;Johnson translator system;digital fi lter

    TU 311

    A

    1007-2861(2017)04-0600-09

    DO I:10.12066/j.issn.1007-2861.1852

    2016-09-30

    國家自然科學(xué)基金資助項(xiàng)目(51378304)

    李春祥(1964—),男,教授,博士生導(dǎo)師,博士,研究方向?yàn)榻Y(jié)構(gòu)振動(dòng)控制、結(jié)構(gòu)風(fēng)工程、結(jié)構(gòu)抗震.

    E-mail:Li-chunxiang@vip.sina.com

    猜你喜歡
    偏度峰度譜分析
    納譜分析技術(shù)(蘇州)有限公司
    色譜(2022年5期)2022-04-28 02:49:10
    擴(kuò)散峰度成像技術(shù)檢測急性期癲癇大鼠模型的成像改變
    對(duì)稱分布的矩刻畫
    磁共振擴(kuò)散峰度成像在肝臟病變中的研究進(jìn)展
    基于自動(dòng)反相校正和峰度值比較的探地雷達(dá)回波信號(hào)去噪方法
    基于偏度的滾動(dòng)軸承聲信號(hào)故障分析方法
    Cr12MoV冷作模具鋼滲鉻層界面能譜分析
    Rotenberg模型中一類遷移算子的譜分析
    沉香GC-MS指紋圖譜分析
    中成藥(2016年8期)2016-05-17 06:08:26
    考慮偏度特征的動(dòng)態(tài)多響應(yīng)穩(wěn)健參數(shù)設(shè)計(jì)與優(yōu)化
    熟女人妻精品中文字幕| 又粗又爽又猛毛片免费看| 内射极品少妇av片p| 亚洲人成网站在线播放欧美日韩| 色哟哟哟哟哟哟| 婷婷精品国产亚洲av在线| 色综合亚洲欧美另类图片| 亚洲内射少妇av| 亚洲内射少妇av| 一级毛片久久久久久久久女| 男插女下体视频免费在线播放| 免费在线观看成人毛片| 亚洲精品乱码久久久v下载方式| 亚洲男人的天堂狠狠| 日韩欧美在线二视频| 国产久久久一区二区三区| 精品久久久久久久人妻蜜臀av| 在线播放国产精品三级| www.熟女人妻精品国产| 国产成人av教育| 99riav亚洲国产免费| 午夜福利在线观看吧| 最近最新免费中文字幕在线| 国产成人影院久久av| 国产亚洲av嫩草精品影院| 亚洲五月婷婷丁香| 99在线视频只有这里精品首页| 禁无遮挡网站| 国产中年淑女户外野战色| 中文字幕精品亚洲无线码一区| 色尼玛亚洲综合影院| 成人永久免费在线观看视频| 午夜福利在线在线| 99国产综合亚洲精品| 在现免费观看毛片| 免费大片18禁| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 亚洲在线自拍视频| 亚洲欧美精品综合久久99| 成熟少妇高潮喷水视频| 亚洲国产精品久久男人天堂| 久久久色成人| 欧美成狂野欧美在线观看| 午夜老司机福利剧场| 亚洲精华国产精华精| 午夜免费激情av| 别揉我奶头~嗯~啊~动态视频| 狂野欧美白嫩少妇大欣赏| 亚洲一区高清亚洲精品| 国产精品伦人一区二区| 欧美日韩亚洲国产一区二区在线观看| 蜜桃久久精品国产亚洲av| 九色国产91popny在线| 色5月婷婷丁香| 97碰自拍视频| 哪里可以看免费的av片| 午夜免费激情av| 一本精品99久久精品77| 99久久久亚洲精品蜜臀av| 男女床上黄色一级片免费看| www.www免费av| 尤物成人国产欧美一区二区三区| 日本黄大片高清| 五月玫瑰六月丁香| 一进一出抽搐动态| 熟妇人妻久久中文字幕3abv| 在线观看一区二区三区| 俺也久久电影网| 久久人妻av系列| 久久国产精品影院| 免费高清视频大片| 国产美女午夜福利| 亚洲国产日韩欧美精品在线观看| 国产精品影院久久| 一区二区三区四区激情视频 | 成熟少妇高潮喷水视频| 久久热精品热| 听说在线观看完整版免费高清| 好看av亚洲va欧美ⅴa在| 激情在线观看视频在线高清| 国产一区二区三区视频了| 国产精品影院久久| 久久这里只有精品中国| 亚洲av免费在线观看| 天堂动漫精品| 男女之事视频高清在线观看| 久久久久九九精品影院| 两人在一起打扑克的视频| 乱码一卡2卡4卡精品| 亚洲最大成人中文| 亚洲,欧美精品.| 午夜福利在线观看吧| 级片在线观看| 午夜福利在线观看免费完整高清在 | 欧美一区二区国产精品久久精品| 两人在一起打扑克的视频| 亚洲无线观看免费| 无遮挡黄片免费观看| 性色av乱码一区二区三区2| 久久精品国产99精品国产亚洲性色| 特级一级黄色大片| 精品久久久久久久久亚洲 | 蜜桃久久精品国产亚洲av| 黄色日韩在线| 国产精品三级大全| 热99re8久久精品国产| 国产精品人妻久久久久久| 国产av在哪里看| 免费观看人在逋| 蜜桃亚洲精品一区二区三区| 亚洲av二区三区四区| 亚洲综合色惰| 国产人妻一区二区三区在| av专区在线播放| ponron亚洲| 淫秽高清视频在线观看| 欧美另类亚洲清纯唯美| 一进一出抽搐gif免费好疼| 国产一区二区在线av高清观看| 中文字幕精品亚洲无线码一区| 亚洲午夜理论影院| 男女下面进入的视频免费午夜| 久久久久国内视频| 九色国产91popny在线| 国产精品综合久久久久久久免费| 免费看美女性在线毛片视频| 国产在视频线在精品| 三级国产精品欧美在线观看| 国产亚洲欧美在线一区二区| 亚洲人成伊人成综合网2020| 国产精品98久久久久久宅男小说| 美女高潮喷水抽搐中文字幕| 麻豆成人午夜福利视频| 成人高潮视频无遮挡免费网站| 男人舔女人下体高潮全视频| 精品欧美国产一区二区三| 午夜视频国产福利| 久久人人精品亚洲av| 精品免费久久久久久久清纯| 啦啦啦观看免费观看视频高清| 久久久久久大精品| 精品午夜福利视频在线观看一区| 国产乱人伦免费视频| 国产欧美日韩一区二区三| 在现免费观看毛片| 国产精品乱码一区二三区的特点| 亚洲激情在线av| 99精品久久久久人妻精品| 少妇的逼水好多| 免费看a级黄色片| 久久亚洲真实| 日日夜夜操网爽| 国产精品久久久久久亚洲av鲁大| 毛片一级片免费看久久久久 | 亚洲中文日韩欧美视频| 国内揄拍国产精品人妻在线| 欧美xxxx黑人xx丫x性爽| 夜夜躁狠狠躁天天躁| 午夜福利在线观看免费完整高清在 | 欧美日韩亚洲国产一区二区在线观看| av在线蜜桃| 淫秽高清视频在线观看| 国产精品久久电影中文字幕| 老熟妇仑乱视频hdxx| 俄罗斯特黄特色一大片| 级片在线观看| 亚洲乱码一区二区免费版| 最近最新免费中文字幕在线| 午夜免费男女啪啪视频观看 | 大型黄色视频在线免费观看| 丝袜美腿在线中文| 乱码一卡2卡4卡精品| 国内少妇人妻偷人精品xxx网站| 国产伦人伦偷精品视频| 看免费av毛片| 97人妻精品一区二区三区麻豆| 一个人观看的视频www高清免费观看| 乱码一卡2卡4卡精品| 在线观看av片永久免费下载| 男女视频在线观看网站免费| 99热精品在线国产| 色综合亚洲欧美另类图片| 欧美乱妇无乱码| 精品久久久久久久人妻蜜臀av| 亚洲自拍偷在线| 黄色视频,在线免费观看| 如何舔出高潮| 麻豆久久精品国产亚洲av| 极品教师在线免费播放| 天堂网av新在线| 大型黄色视频在线免费观看| 免费人成视频x8x8入口观看| 国产探花在线观看一区二区| 国产不卡一卡二| av女优亚洲男人天堂| 久久久久免费精品人妻一区二区| 精华霜和精华液先用哪个| 美女cb高潮喷水在线观看| 国产在线男女| 亚洲综合色惰| 欧美成人免费av一区二区三区| 久久香蕉精品热| 久久久精品欧美日韩精品| 日韩欧美在线二视频| 成人毛片a级毛片在线播放| 一个人免费在线观看电影| 中文在线观看免费www的网站| av在线观看视频网站免费| 欧美黑人巨大hd| 欧美潮喷喷水| 精品一区二区免费观看| 午夜激情欧美在线| 女同久久另类99精品国产91| 琪琪午夜伦伦电影理论片6080| 亚洲第一欧美日韩一区二区三区| 波多野结衣高清无吗| 91久久精品国产一区二区成人| 夜夜躁狠狠躁天天躁| 丁香欧美五月| 国产美女午夜福利| 国产视频一区二区在线看| 色哟哟·www| 国产成人欧美在线观看| 日韩精品青青久久久久久| 日韩av在线大香蕉| 日日摸夜夜添夜夜添av毛片 | 国产成人影院久久av| 免费av观看视频| 婷婷六月久久综合丁香| 女人被狂操c到高潮| 精品一区二区三区视频在线观看免费| 91在线精品国自产拍蜜月| 久久久精品大字幕| 黄色女人牲交| 成人无遮挡网站| 免费观看精品视频网站| 国产伦精品一区二区三区视频9| 国产淫片久久久久久久久 | 国产精品爽爽va在线观看网站| 欧洲精品卡2卡3卡4卡5卡区| 大型黄色视频在线免费观看| 久久九九热精品免费| 国产熟女xx| 亚洲av成人av| 欧美国产日韩亚洲一区| 99国产极品粉嫩在线观看| 日韩高清综合在线| 国产精品一区二区三区四区免费观看 | 两个人的视频大全免费| 欧美xxxx黑人xx丫x性爽| 中文资源天堂在线| 国产精品三级大全| 99热这里只有是精品50| 少妇人妻一区二区三区视频| 午夜激情福利司机影院| 国产一级毛片七仙女欲春2| 最近中文字幕高清免费大全6 | 内射极品少妇av片p| 欧美精品国产亚洲| 日韩欧美国产在线观看| 午夜a级毛片| 在线免费观看不下载黄p国产 | 国产成年人精品一区二区| 91麻豆av在线| 无遮挡黄片免费观看| 精品欧美国产一区二区三| 国语自产精品视频在线第100页| 丰满乱子伦码专区| 国产色婷婷99| 日本精品一区二区三区蜜桃| 尤物成人国产欧美一区二区三区| 久久精品夜夜夜夜夜久久蜜豆| 亚洲成人久久爱视频| 麻豆成人av在线观看| 美女xxoo啪啪120秒动态图 | 日日干狠狠操夜夜爽| av专区在线播放| 丰满人妻一区二区三区视频av| 美女高潮的动态| 亚洲中文字幕一区二区三区有码在线看| 丝袜美腿在线中文| 伊人久久精品亚洲午夜| 亚洲无线在线观看| 亚洲三级黄色毛片| 欧美黄色片欧美黄色片| 国产色爽女视频免费观看| eeuss影院久久| 黄色配什么色好看| 一二三四社区在线视频社区8| 亚洲久久久久久中文字幕| 午夜免费成人在线视频| 一进一出抽搐动态| 两个人的视频大全免费| 亚洲最大成人手机在线| 精品一区二区三区视频在线观看免费| 成人无遮挡网站| 久久久久国产精品人妻aⅴ院| 琪琪午夜伦伦电影理论片6080| 久久人妻av系列| 国产高清激情床上av| 午夜老司机福利剧场| 亚洲无线观看免费| 欧美性感艳星| 赤兔流量卡办理| 美女黄网站色视频| 亚洲男人的天堂狠狠| 午夜影院日韩av| 美女 人体艺术 gogo| 欧美日韩黄片免| 在线观看午夜福利视频| 亚洲经典国产精华液单 | 亚洲av二区三区四区| 制服丝袜大香蕉在线| 啦啦啦观看免费观看视频高清| 免费电影在线观看免费观看| 欧美国产日韩亚洲一区| 搡女人真爽免费视频火全软件 | 一夜夜www| 国产伦精品一区二区三区四那| 亚洲18禁久久av| 色尼玛亚洲综合影院| 欧美bdsm另类| 特大巨黑吊av在线直播| 日本与韩国留学比较| 在线免费观看不下载黄p国产 | 国产三级黄色录像| 国产色婷婷99| 免费高清视频大片| 欧美潮喷喷水| 国产精品三级大全| 看黄色毛片网站| 精品午夜福利在线看| 国产在视频线在精品| 欧美bdsm另类| 99热这里只有精品一区| 一区福利在线观看| 亚洲avbb在线观看| 99久久精品热视频| 黄色丝袜av网址大全| 亚洲美女视频黄频| 国产精品永久免费网站| 九九久久精品国产亚洲av麻豆| 亚洲欧美精品综合久久99| 亚洲美女黄片视频| 首页视频小说图片口味搜索| 国产男靠女视频免费网站| 午夜激情福利司机影院| 国产精品影院久久| 97超级碰碰碰精品色视频在线观看| 在线观看av片永久免费下载| 国产麻豆成人av免费视频| ponron亚洲| 最近在线观看免费完整版| 亚洲国产高清在线一区二区三| 国产精品久久视频播放| 精品不卡国产一区二区三区| 夜夜爽天天搞| 黄色一级大片看看| 久久亚洲真实| 在线观看66精品国产| 欧美成人一区二区免费高清观看| av福利片在线观看| 十八禁人妻一区二区| 黄色日韩在线| 日韩高清综合在线| 熟女电影av网| 午夜福利视频1000在线观看| 国产精品一区二区性色av| 亚洲成av人片免费观看| 每晚都被弄得嗷嗷叫到高潮| 午夜激情欧美在线| 久久国产乱子伦精品免费另类| 首页视频小说图片口味搜索| 久久精品国产自在天天线| 婷婷丁香在线五月| 男女做爰动态图高潮gif福利片| 特大巨黑吊av在线直播| 久久精品夜夜夜夜夜久久蜜豆| 亚洲美女搞黄在线观看 | 人人妻人人澡欧美一区二区| 女人被狂操c到高潮| 丁香欧美五月| 精品99又大又爽又粗少妇毛片 | 免费搜索国产男女视频| 三级男女做爰猛烈吃奶摸视频| 一边摸一边抽搐一进一小说| 欧美日韩综合久久久久久 | 欧美色欧美亚洲另类二区| 中文字幕av成人在线电影| 国产成人aa在线观看| 久久欧美精品欧美久久欧美| 免费在线观看成人毛片| 2021天堂中文幕一二区在线观| 欧美色欧美亚洲另类二区| 人人妻人人澡欧美一区二区| 国产精品伦人一区二区| 国内揄拍国产精品人妻在线| 久久国产精品人妻蜜桃| 免费高清视频大片| 亚洲成av人片在线播放无| a级毛片免费高清观看在线播放| 免费av观看视频| 久久久精品欧美日韩精品| 午夜免费男女啪啪视频观看 | 亚洲18禁久久av| 高潮久久久久久久久久久不卡| 亚洲av免费在线观看| 欧美最黄视频在线播放免费| 综合色av麻豆| 亚洲va日本ⅴa欧美va伊人久久| 在线观看美女被高潮喷水网站 | a级毛片a级免费在线| 午夜视频国产福利| 搡老妇女老女人老熟妇| 又黄又爽又免费观看的视频| 中国美女看黄片| 熟女电影av网| 亚洲成人久久爱视频| 国产三级中文精品| 精品午夜福利视频在线观看一区| 午夜福利欧美成人| 国内久久婷婷六月综合欲色啪| 一进一出好大好爽视频| 日本一本二区三区精品| 国产精品三级大全| 国产免费男女视频| 日本 欧美在线| 99视频精品全部免费 在线| 青草久久国产| 中文字幕av成人在线电影| 97超级碰碰碰精品色视频在线观看| 亚洲乱码一区二区免费版| 老司机福利观看| 亚洲综合色惰| 99热精品在线国产| 超碰av人人做人人爽久久| 亚洲最大成人中文| www.www免费av| 国产一区二区三区视频了| 三级毛片av免费| 国产男靠女视频免费网站| 特大巨黑吊av在线直播| 欧美午夜高清在线| 脱女人内裤的视频| 听说在线观看完整版免费高清| 久久久久国产精品人妻aⅴ院| 亚洲精品影视一区二区三区av| 九九在线视频观看精品| 别揉我奶头~嗯~啊~动态视频| 国产av在哪里看| 99久国产av精品| 精品国产三级普通话版| 亚洲片人在线观看| 久久伊人香网站| 免费av毛片视频| 国产大屁股一区二区在线视频| 久久精品影院6| 色综合婷婷激情| 欧美成人性av电影在线观看| 女同久久另类99精品国产91| 18禁黄网站禁片午夜丰满| 国产av在哪里看| 欧美成人免费av一区二区三区| 久久精品国产清高在天天线| 免费在线观看日本一区| 制服丝袜大香蕉在线| 亚洲色图av天堂| 亚洲av免费高清在线观看| 91久久精品电影网| 美女高潮的动态| avwww免费| 成人美女网站在线观看视频| 琪琪午夜伦伦电影理论片6080| 内射极品少妇av片p| 亚洲精华国产精华精| 1000部很黄的大片| 日本黄色片子视频| 香蕉av资源在线| 又黄又爽又刺激的免费视频.| 亚洲无线在线观看| 国产毛片a区久久久久| 一本久久中文字幕| 一进一出抽搐gif免费好疼| 国产精品亚洲av一区麻豆| 午夜福利欧美成人| 日韩大尺度精品在线看网址| 桃红色精品国产亚洲av| 午夜福利在线在线| 直男gayav资源| 麻豆国产97在线/欧美| 欧美绝顶高潮抽搐喷水| a在线观看视频网站| 日韩人妻高清精品专区| 欧美中文日本在线观看视频| 精品乱码久久久久久99久播| 波多野结衣巨乳人妻| 久久久久久久精品吃奶| 午夜两性在线视频| x7x7x7水蜜桃| 国内揄拍国产精品人妻在线| 日日夜夜操网爽| 少妇裸体淫交视频免费看高清| 国产精品1区2区在线观看.| 久久午夜亚洲精品久久| 亚洲乱码一区二区免费版| 日本在线视频免费播放| a级毛片免费高清观看在线播放| 天天一区二区日本电影三级| 老鸭窝网址在线观看| 成年女人看的毛片在线观看| 欧美色欧美亚洲另类二区| 国产主播在线观看一区二区| 窝窝影院91人妻| 一夜夜www| 国产午夜福利久久久久久| 在线看三级毛片| 免费av毛片视频| 真人做人爱边吃奶动态| 欧美成人a在线观看| 日韩欧美 国产精品| 老司机福利观看| 国产精品一区二区三区四区久久| 免费看光身美女| 午夜老司机福利剧场| 一进一出抽搐gif免费好疼| 国产精品伦人一区二区| 欧美精品国产亚洲| 一个人看视频在线观看www免费| 久久性视频一级片| 村上凉子中文字幕在线| 中文在线观看免费www的网站| 成人特级av手机在线观看| 成人国产综合亚洲| 别揉我奶头~嗯~啊~动态视频| 精品久久久久久久久亚洲 | 99久久久亚洲精品蜜臀av| 久久精品国产99精品国产亚洲性色| 国产精品一及| 韩国av一区二区三区四区| 18美女黄网站色大片免费观看| 男人舔奶头视频| 国产久久久一区二区三区| 成人特级av手机在线观看| 老熟妇乱子伦视频在线观看| 18禁黄网站禁片午夜丰满| 男人舔女人下体高潮全视频| 亚洲av中文字字幕乱码综合| 免费搜索国产男女视频| 此物有八面人人有两片| 亚洲美女搞黄在线观看 | 日本精品一区二区三区蜜桃| 欧美+日韩+精品| 久久久久九九精品影院| 少妇人妻一区二区三区视频| 中文字幕久久专区| 久久6这里有精品| 老熟妇仑乱视频hdxx| 欧美日韩亚洲国产一区二区在线观看| 国产三级黄色录像| 婷婷精品国产亚洲av| 精品无人区乱码1区二区| 色综合亚洲欧美另类图片| 中文资源天堂在线| 日本一二三区视频观看| 亚洲精品日韩av片在线观看| 国产在线精品亚洲第一网站| 男女之事视频高清在线观看| av在线蜜桃| 免费黄网站久久成人精品 | 亚洲国产精品合色在线| 99精品久久久久人妻精品| 亚洲美女视频黄频| 人人妻人人看人人澡| 桃红色精品国产亚洲av| 亚洲成人中文字幕在线播放| 99热精品在线国产| 最近最新免费中文字幕在线| 久久欧美精品欧美久久欧美| 九色国产91popny在线| 少妇裸体淫交视频免费看高清| 国产一级毛片七仙女欲春2| 国产精品久久久久久久久免 | 国产大屁股一区二区在线视频| 听说在线观看完整版免费高清| 午夜福利在线观看免费完整高清在 | 又粗又爽又猛毛片免费看| 久久精品国产亚洲av天美| 极品教师在线免费播放| 亚洲精品在线美女| 久久草成人影院| 国产av不卡久久| av天堂中文字幕网| 国产日本99.免费观看| 在线免费观看的www视频| 一本综合久久免费| 成人av在线播放网站| 一a级毛片在线观看| 99久久精品一区二区三区| 最近最新免费中文字幕在线| 午夜福利18| 免费观看精品视频网站| 啪啪无遮挡十八禁网站| 最好的美女福利视频网| 久久久精品大字幕| 亚洲精华国产精华精| 日本五十路高清| 国产探花在线观看一区二区| 嫩草影院新地址| 精品人妻偷拍中文字幕| 久久伊人香网站| 在线观看美女被高潮喷水网站 | 午夜免费激情av| 免费看a级黄色片| 美女被艹到高潮喷水动态| 欧美xxxx黑人xx丫x性爽|