• <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片永久免费下载| 女生性感内裤真人,穿戴方法视频| 99久久无色码亚洲精品果冻| 校园春色视频在线观看| 搞女人的毛片| 免费看a级黄色片| 深夜精品福利| 亚洲国产色片| 大又大粗又爽又黄少妇毛片口| 亚洲色图av天堂| 国产精品久久久久久亚洲av鲁大| 亚洲精品国产成人久久av| 午夜福利视频1000在线观看| 国产在视频线在精品| 精品福利观看| avwww免费| 国产黄色小视频在线观看| 久久人人爽人人片av| av专区在线播放| 永久网站在线| 免费不卡的大黄色大毛片视频在线观看 | 国产一区亚洲一区在线观看| 最近视频中文字幕2019在线8| av在线播放精品| 男女视频在线观看网站免费| 99久久成人亚洲精品观看| 国产毛片a区久久久久| 欧美日韩在线观看h| 中文字幕熟女人妻在线| 欧美色欧美亚洲另类二区| 黄色欧美视频在线观看| 天天躁夜夜躁狠狠久久av| 国产欧美日韩精品亚洲av| 最近手机中文字幕大全| 国产成人a区在线观看| 欧美日韩综合久久久久久| 最新中文字幕久久久久| 97热精品久久久久久| 中文字幕免费在线视频6| 国产探花极品一区二区| 一级a爱片免费观看的视频| 日韩欧美免费精品| 日本黄色片子视频| 黄色一级大片看看| 亚洲av.av天堂| 国产精品久久视频播放| 成人美女网站在线观看视频| 日韩精品青青久久久久久| 色哟哟·www| 国产片特级美女逼逼视频| 午夜精品一区二区三区免费看| 全区人妻精品视频| 在线免费观看的www视频| 一级毛片我不卡| 久久精品国产鲁丝片午夜精品| 老熟妇乱子伦视频在线观看| 99热网站在线观看| 简卡轻食公司| 搡老熟女国产l中国老女人| 亚洲精品一区av在线观看| 国产一区二区三区在线臀色熟女| 国产中年淑女户外野战色| 精品久久久久久久人妻蜜臀av| 老熟妇仑乱视频hdxx| 亚洲精品粉嫩美女一区| 国产成人aa在线观看| 成人欧美大片| 婷婷精品国产亚洲av在线| 噜噜噜噜噜久久久久久91| 最近的中文字幕免费完整| 天堂av国产一区二区熟女人妻| 国产高清视频在线播放一区| 听说在线观看完整版免费高清| 一级黄片播放器| 日韩在线高清观看一区二区三区| 少妇被粗大猛烈的视频| 99久国产av精品国产电影| 亚洲在线自拍视频| 悠悠久久av| 精品久久久久久久久av| 少妇熟女欧美另类| 精品福利观看| 国产精品一区二区三区四区久久| 午夜激情福利司机影院| 日韩,欧美,国产一区二区三区 | 国内少妇人妻偷人精品xxx网站| 天堂√8在线中文| 人妻少妇偷人精品九色| 久久久久久久久久成人| 久久久成人免费电影| 成人无遮挡网站| 日韩欧美国产在线观看| 午夜福利高清视频| 校园人妻丝袜中文字幕| 亚洲不卡免费看| 国产高清不卡午夜福利| 国产av不卡久久| 国产精品美女特级片免费视频播放器| 国产免费男女视频| 草草在线视频免费看| 午夜福利在线观看吧| 中文在线观看免费www的网站| 全区人妻精品视频| 听说在线观看完整版免费高清| 国产单亲对白刺激| 成人高潮视频无遮挡免费网站| 男人舔女人下体高潮全视频| 精品久久久久久久末码| 欧美潮喷喷水| 波多野结衣高清无吗| 99九九线精品视频在线观看视频| 亚洲国产精品久久男人天堂| 国产高潮美女av| 国产精品三级大全| 女的被弄到高潮叫床怎么办| 你懂的网址亚洲精品在线观看 | 搡女人真爽免费视频火全软件 | 亚洲中文日韩欧美视频| 此物有八面人人有两片| 亚洲国产精品成人综合色| 午夜久久久久精精品| 午夜福利在线观看免费完整高清在 | 九色成人免费人妻av| 国产黄色视频一区二区在线观看 | 啦啦啦观看免费观看视频高清| 久久久久久九九精品二区国产| 最近视频中文字幕2019在线8| 如何舔出高潮| 欧美最新免费一区二区三区| 99热网站在线观看| 亚洲成a人片在线一区二区| 国产精品乱码一区二三区的特点| 亚州av有码| 熟女人妻精品中文字幕| 免费无遮挡裸体视频| 91麻豆精品激情在线观看国产| 亚洲av免费高清在线观看| 中文字幕av在线有码专区| 色在线成人网| 高清毛片免费观看视频网站| 99热精品在线国产| 国产亚洲精品久久久久久毛片| av中文乱码字幕在线| 国产高清激情床上av| 亚洲av一区综合| 桃色一区二区三区在线观看| 亚洲国产精品久久男人天堂| 亚洲中文日韩欧美视频| 欧美激情久久久久久爽电影| 日本 av在线| 春色校园在线视频观看| 熟妇人妻久久中文字幕3abv| 亚洲一级一片aⅴ在线观看| 精品一区二区三区视频在线| 熟妇人妻久久中文字幕3abv| 深夜精品福利| 麻豆精品久久久久久蜜桃| 人妻制服诱惑在线中文字幕| 哪里可以看免费的av片| 久久99热6这里只有精品| 精品福利观看| 亚洲最大成人av| 亚洲人成网站在线播| 国产精品久久久久久久电影| 国产精品女同一区二区软件| 又黄又爽又刺激的免费视频.| av.在线天堂| 亚洲无线在线观看| 99久久成人亚洲精品观看| 欧美一区二区精品小视频在线| 一个人免费在线观看电影| 亚洲国产精品国产精品| 国产精品久久久久久亚洲av鲁大| 亚洲婷婷狠狠爱综合网| 日日摸夜夜添夜夜添小说| 亚洲成人精品中文字幕电影| 1024手机看黄色片| 看免费成人av毛片| 长腿黑丝高跟| 99久国产av精品| 高清午夜精品一区二区三区 | 乱系列少妇在线播放| 欧美激情国产日韩精品一区| 日韩三级伦理在线观看| 欧美成人精品欧美一级黄| 成人特级黄色片久久久久久久| 狂野欧美激情性xxxx在线观看| 九九爱精品视频在线观看| 91久久精品电影网| 一边摸一边抽搐一进一小说| 色5月婷婷丁香| 高清日韩中文字幕在线| 亚洲精品456在线播放app| 熟女电影av网| 成人漫画全彩无遮挡| 欧美另类亚洲清纯唯美| 精品一区二区三区人妻视频| 观看美女的网站| 久久精品国产99精品国产亚洲性色| 久久亚洲国产成人精品v| 久久精品国产亚洲av涩爱 | 日本熟妇午夜| 欧美+亚洲+日韩+国产| 免费高清视频大片| 欧美成人一区二区免费高清观看| 亚洲av不卡在线观看| 搞女人的毛片| 亚洲乱码一区二区免费版| 色哟哟·www| 国产美女午夜福利| 国产精品一区二区免费欧美| 免费av不卡在线播放| 精品国产三级普通话版| 亚洲五月天丁香| 色尼玛亚洲综合影院| 日韩,欧美,国产一区二区三区 | 九九久久精品国产亚洲av麻豆| 乱系列少妇在线播放| 亚洲色图av天堂| 波多野结衣高清作品| 亚洲欧美日韩高清专用| 最近最新中文字幕大全电影3| 欧美区成人在线视频| 一进一出好大好爽视频| 亚洲精品国产av成人精品 | 香蕉av资源在线| 最近在线观看免费完整版| 国产探花在线观看一区二区| 22中文网久久字幕| 狂野欧美激情性xxxx在线观看| 国产视频内射| 国产精品99久久久久久久久| 菩萨蛮人人尽说江南好唐韦庄 | 噜噜噜噜噜久久久久久91| 天堂影院成人在线观看| 2021天堂中文幕一二区在线观| 日本一二三区视频观看| 国语自产精品视频在线第100页| 亚洲国产欧美人成| 午夜精品国产一区二区电影 | 国内揄拍国产精品人妻在线| 久久99热6这里只有精品| 亚洲第一区二区三区不卡| av在线播放精品| 亚洲av中文字字幕乱码综合| 黑人高潮一二区| 干丝袜人妻中文字幕| 免费在线观看影片大全网站| 99riav亚洲国产免费| 搡老熟女国产l中国老女人| 精品午夜福利在线看| 色av中文字幕| 亚洲欧美日韩无卡精品| 久久精品人妻少妇| 国内精品宾馆在线| 国产精品1区2区在线观看.| 美女xxoo啪啪120秒动态图| 国产在线精品亚洲第一网站| 久久久成人免费电影| 啦啦啦啦在线视频资源| 亚洲激情五月婷婷啪啪| 麻豆成人午夜福利视频| 国产精品久久久久久久电影| 欧美成人一区二区免费高清观看| 成年版毛片免费区| 欧美性猛交╳xxx乱大交人| 国产成人a区在线观看| 男女下面进入的视频免费午夜| 秋霞在线观看毛片| 寂寞人妻少妇视频99o| 91久久精品国产一区二区成人| 午夜精品一区二区三区免费看| 成人国产麻豆网| 亚洲成a人片在线一区二区| 人妻夜夜爽99麻豆av| 美女内射精品一级片tv| 亚洲国产精品合色在线| 精品免费久久久久久久清纯| 赤兔流量卡办理| 日韩欧美三级三区| 国产乱人偷精品视频| 欧美日韩在线观看h| 能在线免费观看的黄片| 免费在线观看影片大全网站| 亚洲最大成人手机在线| 你懂的网址亚洲精品在线观看 | 看十八女毛片水多多多| 热99re8久久精品国产| 国产毛片a区久久久久| 国产成人福利小说| 成年av动漫网址| 12—13女人毛片做爰片一| 一区二区三区四区激情视频 | 国产 一区 欧美 日韩| 村上凉子中文字幕在线| 天天躁日日操中文字幕| av专区在线播放| 亚洲欧美成人综合另类久久久 | 美女大奶头视频| 最近手机中文字幕大全| 亚洲精华国产精华液的使用体验 | 2021天堂中文幕一二区在线观| 亚洲五月天丁香| 麻豆一二三区av精品| 日日啪夜夜撸| 日日摸夜夜添夜夜爱| 欧美高清成人免费视频www| 国产欧美日韩精品亚洲av| 亚洲天堂国产精品一区在线| 热99在线观看视频| 欧美另类亚洲清纯唯美| 日韩精品有码人妻一区| 一卡2卡三卡四卡精品乱码亚洲| 亚洲成人精品中文字幕电影| 女人被狂操c到高潮| 久久热精品热| 在线免费观看不下载黄p国产| 1000部很黄的大片| 亚洲欧美清纯卡通| 插逼视频在线观看| 最后的刺客免费高清国语| 国内精品一区二区在线观看| 偷拍熟女少妇极品色| 久久人妻av系列| 国产男靠女视频免费网站| 精品久久久久久久末码| 日日摸夜夜添夜夜添av毛片| 精品乱码久久久久久99久播| 亚洲精品456在线播放app| 91久久精品电影网| 日本a在线网址| 亚洲精品亚洲一区二区| 亚洲欧美清纯卡通| 日本-黄色视频高清免费观看| 欧美成人a在线观看| 人妻制服诱惑在线中文字幕| 亚洲aⅴ乱码一区二区在线播放| 国产伦在线观看视频一区| 中文字幕熟女人妻在线| 91狼人影院| 国产精品爽爽va在线观看网站| 中文字幕久久专区| 长腿黑丝高跟| a级一级毛片免费在线观看| АⅤ资源中文在线天堂| 欧美潮喷喷水| 91久久精品国产一区二区三区| 免费观看精品视频网站| 变态另类丝袜制服| 少妇的逼水好多| 有码 亚洲区| 国产精品电影一区二区三区| 欧美一区二区国产精品久久精品| 免费av观看视频| 亚洲无线在线观看| av在线亚洲专区| 俄罗斯特黄特色一大片| 国产精品99久久久久久久久| 亚洲av熟女| 欧美激情久久久久久爽电影| av卡一久久| 人人妻人人澡人人爽人人夜夜 | 啦啦啦韩国在线观看视频| 99热只有精品国产| 久久久精品大字幕| 有码 亚洲区| 日日摸夜夜添夜夜添小说| 亚洲精品粉嫩美女一区| 免费av观看视频| or卡值多少钱| 欧美绝顶高潮抽搐喷水| 一区二区三区四区激情视频 | 男女做爰动态图高潮gif福利片| 色av中文字幕| 国产蜜桃级精品一区二区三区| 欧美国产日韩亚洲一区| 国产私拍福利视频在线观看| 久久久久久久亚洲中文字幕| 五月玫瑰六月丁香| 国产探花在线观看一区二区| a级毛片a级免费在线| a级一级毛片免费在线观看| 欧美一区二区国产精品久久精品| 中国美白少妇内射xxxbb| 欧美一区二区亚洲| 日韩,欧美,国产一区二区三区 | av在线天堂中文字幕| 色5月婷婷丁香| 观看美女的网站| 天天一区二区日本电影三级| 美女内射精品一级片tv| 成年女人看的毛片在线观看| 亚洲成人久久爱视频| 22中文网久久字幕| 麻豆av噜噜一区二区三区| 国产精品福利在线免费观看| 国产精品久久久久久久久免| 久久精品国产亚洲av香蕉五月| 尤物成人国产欧美一区二区三区| 亚洲av美国av| 可以在线观看毛片的网站| 丰满乱子伦码专区| 国产成人一区二区在线| 卡戴珊不雅视频在线播放| .国产精品久久| 高清毛片免费观看视频网站| 一进一出抽搐gif免费好疼| 狂野欧美白嫩少妇大欣赏| 亚洲真实伦在线观看| 久久精品国产99精品国产亚洲性色| 97人妻精品一区二区三区麻豆| 校园春色视频在线观看| 99久久中文字幕三级久久日本| 国产又黄又爽又无遮挡在线| 国产欧美日韩精品亚洲av| 真人做人爱边吃奶动态| 日韩成人av中文字幕在线观看 | 久久午夜亚洲精品久久| 久久久久久九九精品二区国产| 韩国av在线不卡| 少妇人妻一区二区三区视频| 蜜桃久久精品国产亚洲av| 99热网站在线观看| 啦啦啦韩国在线观看视频| 午夜免费激情av| 国产成人精品久久久久久| 精品一区二区三区视频在线| 亚洲真实伦在线观看| 国产精品久久久久久久久免| 99热网站在线观看| av专区在线播放| 亚洲性夜色夜夜综合| 欧美+日韩+精品| 色播亚洲综合网| 熟女电影av网| 亚洲成人久久爱视频| 2021天堂中文幕一二区在线观| 特大巨黑吊av在线直播| 国产精品久久久久久精品电影| 亚洲在线自拍视频| 久久久精品欧美日韩精品| 日本熟妇午夜| 亚洲欧美清纯卡通| 波多野结衣高清无吗| 蜜桃亚洲精品一区二区三区| 久久精品国产清高在天天线| 国产激情偷乱视频一区二区| 国产精品女同一区二区软件| 免费av不卡在线播放| 亚洲av免费高清在线观看| 嫩草影视91久久| 又爽又黄无遮挡网站| 久久这里只有精品中国| 淫妇啪啪啪对白视频| 麻豆av噜噜一区二区三区| 国产av在哪里看| 亚洲av不卡在线观看| 国产私拍福利视频在线观看| 波多野结衣巨乳人妻| 一夜夜www| 波野结衣二区三区在线| 亚洲欧美日韩高清专用| 日本与韩国留学比较| 免费黄网站久久成人精品| 日韩三级伦理在线观看| 国产精品不卡视频一区二区| 欧美在线一区亚洲| 亚洲电影在线观看av| 一进一出抽搐动态| 黄片wwwwww| 国产高清三级在线| 黄色视频,在线免费观看| 老熟妇乱子伦视频在线观看| 国产欧美日韩一区二区精品| 2021天堂中文幕一二区在线观| 欧洲精品卡2卡3卡4卡5卡区| 床上黄色一级片| av在线亚洲专区| 亚洲婷婷狠狠爱综合网| 中文字幕av成人在线电影| 99久国产av精品| 亚洲一级一片aⅴ在线观看| 真人做人爱边吃奶动态| 日韩av在线大香蕉| 国产女主播在线喷水免费视频网站 | 欧美区成人在线视频| 乱系列少妇在线播放| 18禁黄网站禁片免费观看直播| 国产综合懂色| 亚洲无线在线观看| h日本视频在线播放| 精品一区二区三区视频在线| 97人妻精品一区二区三区麻豆| 国产又黄又爽又无遮挡在线| or卡值多少钱| 欧美最新免费一区二区三区| 亚洲人与动物交配视频| 久久精品夜夜夜夜夜久久蜜豆| 成人综合一区亚洲| 我要搜黄色片| 亚洲电影在线观看av| 日本撒尿小便嘘嘘汇集6| 午夜a级毛片| 啦啦啦韩国在线观看视频| 日本爱情动作片www.在线观看 | 欧美色欧美亚洲另类二区| 搡女人真爽免费视频火全软件 | 国产私拍福利视频在线观看| 日韩在线高清观看一区二区三区| 午夜福利在线观看吧| 国内少妇人妻偷人精品xxx网站| 特大巨黑吊av在线直播| 久久人人爽人人爽人人片va| 久久久久国产网址| 麻豆成人午夜福利视频| 看非洲黑人一级黄片| 免费观看精品视频网站| 两个人视频免费观看高清| 国内精品美女久久久久久| 久久久国产成人精品二区| 亚洲一级一片aⅴ在线观看| 日韩成人伦理影院| 午夜免费激情av| 蜜臀久久99精品久久宅男| 久久精品国产亚洲av天美| 国产国拍精品亚洲av在线观看| 久久久久久久久久成人| 久久午夜亚洲精品久久| 91av网一区二区| 国产av一区在线观看免费| 国产一区二区激情短视频| 看免费成人av毛片| 成人性生交大片免费视频hd| 日韩精品有码人妻一区| 成人精品一区二区免费| 一区二区三区免费毛片| 国产av在哪里看| 国产又黄又爽又无遮挡在线| 日韩欧美免费精品| 成年女人看的毛片在线观看| 老司机午夜福利在线观看视频| 精品午夜福利在线看| 不卡视频在线观看欧美| av女优亚洲男人天堂| 色噜噜av男人的天堂激情| 美女被艹到高潮喷水动态| 色播亚洲综合网| 大型黄色视频在线免费观看| 国产精品永久免费网站| 91久久精品国产一区二区三区| 一级黄色大片毛片| 秋霞在线观看毛片| 国产精品,欧美在线| 亚洲精华国产精华液的使用体验 | 亚洲av成人精品一区久久| 一个人免费在线观看电影| 国产精品一区二区性色av| 亚洲av不卡在线观看| 男女边吃奶边做爰视频| av天堂中文字幕网| 夜夜夜夜夜久久久久| 国产一区亚洲一区在线观看| 成人高潮视频无遮挡免费网站| 别揉我奶头~嗯~啊~动态视频| 一区福利在线观看| 麻豆一二三区av精品| 97热精品久久久久久| 日本与韩国留学比较| 午夜精品一区二区三区免费看| 女人十人毛片免费观看3o分钟| 成人特级黄色片久久久久久久| 成年女人毛片免费观看观看9| 精品午夜福利在线看| 少妇高潮的动态图| 久久人人爽人人片av| 97超视频在线观看视频| 国产免费男女视频| 简卡轻食公司| 免费av毛片视频| 亚洲无线观看免费| 最近在线观看免费完整版| 亚洲欧美成人综合另类久久久 | 成人一区二区视频在线观看| 六月丁香七月| 国产av在哪里看| 久久久精品大字幕| 久久韩国三级中文字幕| 午夜福利在线在线| 国产 一区 欧美 日韩| 亚洲七黄色美女视频| 日韩一本色道免费dvd| 中文在线观看免费www的网站| 男女之事视频高清在线观看| 悠悠久久av| 亚洲国产精品国产精品| 一个人免费在线观看电影| 国内少妇人妻偷人精品xxx网站| av黄色大香蕉| 日本a在线网址| 亚洲激情五月婷婷啪啪| 十八禁国产超污无遮挡网站| 亚洲七黄色美女视频| 嫩草影视91久久| 久久久久久久亚洲中文字幕| 免费搜索国产男女视频| 免费看a级黄色片| 联通29元200g的流量卡| videossex国产| 国产精品无大码| 亚洲av中文字字幕乱码综合| 春色校园在线视频观看| 午夜a级毛片| 国产伦一二天堂av在线观看| 能在线免费观看的黄片| 亚洲一级一片aⅴ在线观看|