• <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)化
    亚洲精品久久久久久婷婷小说| 日日啪夜夜爽| 中文资源天堂在线| 欧美精品一区二区大全| 日韩不卡一区二区三区视频在线| 国产精品一区二区在线不卡| 99热6这里只有精品| 97在线人人人人妻| 国产成人精品福利久久| 欧美日韩在线观看h| 大片电影免费在线观看免费| 久久国内精品自在自线图片| 涩涩av久久男人的天堂| 亚洲国产高清在线一区二区三| 成人黄色视频免费在线看| 多毛熟女@视频| 亚洲aⅴ乱码一区二区在线播放| 五月天丁香电影| 蜜桃久久精品国产亚洲av| 18+在线观看网站| 精品国产露脸久久av麻豆| 91aial.com中文字幕在线观看| 九色成人免费人妻av| 夜夜骑夜夜射夜夜干| 婷婷色av中文字幕| 18禁动态无遮挡网站| 黄片无遮挡物在线观看| 永久免费av网站大全| 成年美女黄网站色视频大全免费 | av免费在线看不卡| 亚洲天堂av无毛| 亚洲美女视频黄频| 黑人高潮一二区| 日韩中文字幕视频在线看片 | 韩国高清视频一区二区三区| 成人美女网站在线观看视频| 欧美xxⅹ黑人| 自拍欧美九色日韩亚洲蝌蚪91 | 久久久久久久大尺度免费视频| 成人午夜精彩视频在线观看| 欧美xxⅹ黑人| 日日撸夜夜添| 插阴视频在线观看视频| 免费黄频网站在线观看国产| 中文天堂在线官网| 狂野欧美激情性xxxx在线观看| 最后的刺客免费高清国语| 亚洲av综合色区一区| 久久av网站| 91午夜精品亚洲一区二区三区| 国产精品av视频在线免费观看| 97超视频在线观看视频| 五月玫瑰六月丁香| 国产精品欧美亚洲77777| 亚洲四区av| 99热这里只有是精品50| 国产免费一级a男人的天堂| 国产精品秋霞免费鲁丝片| 亚洲精品久久午夜乱码| 国产69精品久久久久777片| 大又大粗又爽又黄少妇毛片口| av.在线天堂| 日韩大片免费观看网站| 国产黄色视频一区二区在线观看| 亚洲精品一二三| 久久热精品热| 国产精品精品国产色婷婷| 国产精品一二三区在线看| 国产一区有黄有色的免费视频| 在线观看人妻少妇| 亚洲精品视频女| 久久人妻熟女aⅴ| 免费人妻精品一区二区三区视频| 日本色播在线视频| 久久青草综合色| 国产高清国产精品国产三级 | 黄色配什么色好看| 亚洲综合色惰| 高清欧美精品videossex| 国产乱人偷精品视频| 网址你懂的国产日韩在线| 美女主播在线视频| 国产免费一区二区三区四区乱码| 少妇精品久久久久久久| 亚洲精品aⅴ在线观看| 欧美精品国产亚洲| 熟女人妻精品中文字幕| 精品一品国产午夜福利视频| 中文天堂在线官网| av一本久久久久| 国产精品国产三级国产专区5o| 国产伦精品一区二区三区四那| 日韩av在线免费看完整版不卡| 久久ye,这里只有精品| 99热6这里只有精品| 国产毛片在线视频| 亚洲欧洲国产日韩| 亚洲欧美清纯卡通| 久久久精品免费免费高清| 中国三级夫妇交换| 春色校园在线视频观看| av专区在线播放| 亚洲精品亚洲一区二区| 久久99热这里只频精品6学生| 一级黄片播放器| 熟女av电影| 亚洲熟女精品中文字幕| 美女内射精品一级片tv| 中文字幕久久专区| 日韩人妻高清精品专区| 视频中文字幕在线观看| 国产91av在线免费观看| 亚洲经典国产精华液单| 男女下面进入的视频免费午夜| 超碰av人人做人人爽久久| 久久99热6这里只有精品| 一个人看视频在线观看www免费| 五月伊人婷婷丁香| 美女内射精品一级片tv| 欧美变态另类bdsm刘玥| 免费播放大片免费观看视频在线观看| 一区二区三区免费毛片| 久久精品熟女亚洲av麻豆精品| 日韩人妻高清精品专区| 亚洲aⅴ乱码一区二区在线播放| 又爽又黄a免费视频| 国产熟女欧美一区二区| 久久久精品94久久精品| 中文字幕制服av| 国产精品一区二区在线观看99| 国产成人免费观看mmmm| 国产美女午夜福利| 国产白丝娇喘喷水9色精品| 我的女老师完整版在线观看| 中文乱码字字幕精品一区二区三区| 日韩电影二区| 国产 一区精品| 啦啦啦视频在线资源免费观看| 国产日韩欧美亚洲二区| 国产亚洲5aaaaa淫片| 亚洲精品久久午夜乱码| 国产一区有黄有色的免费视频| 亚洲精品视频女| 天美传媒精品一区二区| 乱码一卡2卡4卡精品| 男女免费视频国产| 亚洲av电影在线观看一区二区三区| 亚洲va在线va天堂va国产| 一本久久精品| 高清午夜精品一区二区三区| 美女脱内裤让男人舔精品视频| 少妇被粗大猛烈的视频| 国产男人的电影天堂91| 久久韩国三级中文字幕| 性色av一级| 美女中出高潮动态图| 人妻制服诱惑在线中文字幕| 永久免费av网站大全| 777米奇影视久久| 精品久久久精品久久久| 日本黄色片子视频| 欧美97在线视频| 五月开心婷婷网| 日韩av不卡免费在线播放| 蜜臀久久99精品久久宅男| 视频区图区小说| 精品久久久久久久末码| 精品一品国产午夜福利视频| 亚洲美女视频黄频| 国产人妻一区二区三区在| 熟女av电影| 国产伦精品一区二区三区视频9| 亚洲精品国产av成人精品| 成年女人在线观看亚洲视频| 亚洲av成人精品一二三区| 热99国产精品久久久久久7| 免费人成在线观看视频色| 亚洲欧美精品专区久久| 丰满少妇做爰视频| 国产又色又爽无遮挡免| 亚洲伊人久久精品综合| 狂野欧美白嫩少妇大欣赏| 国产成人一区二区在线| 高清日韩中文字幕在线| 狂野欧美激情性bbbbbb| 一个人看的www免费观看视频| 色吧在线观看| 干丝袜人妻中文字幕| 亚洲,欧美,日韩| 伦理电影大哥的女人| 日韩 亚洲 欧美在线| 麻豆精品久久久久久蜜桃| 日韩一区二区视频免费看| 自拍偷自拍亚洲精品老妇| 91精品伊人久久大香线蕉| 精品久久久噜噜| 日产精品乱码卡一卡2卡三| 大片免费播放器 马上看| 久久精品熟女亚洲av麻豆精品| 亚洲成人中文字幕在线播放| 国产片特级美女逼逼视频| 人妻一区二区av| 美女国产视频在线观看| 亚洲成人一二三区av| 欧美日韩综合久久久久久| 国产精品爽爽va在线观看网站| 欧美高清性xxxxhd video| 寂寞人妻少妇视频99o| 亚洲av电影在线观看一区二区三区| 纯流量卡能插随身wifi吗| 欧美三级亚洲精品| 国产成人a∨麻豆精品| 又黄又爽又刺激的免费视频.| 国产爱豆传媒在线观看| 国模一区二区三区四区视频| 亚洲精华国产精华液的使用体验| 丰满人妻一区二区三区视频av| 久久久久性生活片| 久久亚洲国产成人精品v| 边亲边吃奶的免费视频| 国产精品国产三级国产av玫瑰| 欧美zozozo另类| av在线老鸭窝| 少妇被粗大猛烈的视频| 久久久久性生活片| 欧美精品一区二区大全| 欧美日韩国产mv在线观看视频 | 国产亚洲91精品色在线| 欧美成人精品欧美一级黄| 免费观看av网站的网址| 精品少妇久久久久久888优播| 国产成人a区在线观看| 亚洲欧美精品专区久久| 亚洲高清免费不卡视频| 中文字幕人妻熟人妻熟丝袜美| 十八禁网站网址无遮挡 | 最后的刺客免费高清国语| 一区二区三区乱码不卡18| 涩涩av久久男人的天堂| 久久av网站| 欧美亚洲 丝袜 人妻 在线| 精品午夜福利在线看| 插阴视频在线观看视频| 久久久欧美国产精品| 亚洲综合色惰| 午夜老司机福利剧场| 精品亚洲乱码少妇综合久久| 麻豆成人av视频| 亚洲精品,欧美精品| 一区二区三区四区激情视频| 国产又色又爽无遮挡免| 国产成人精品婷婷| 免费黄网站久久成人精品| 国产永久视频网站| 97热精品久久久久久| h视频一区二区三区| 精品久久久久久久久亚洲| 高清欧美精品videossex| 人妻制服诱惑在线中文字幕| 国产精品精品国产色婷婷| 极品少妇高潮喷水抽搐| 一级毛片黄色毛片免费观看视频| 精品久久久久久电影网| 成人亚洲精品一区在线观看 | 国产精品偷伦视频观看了| 成人国产av品久久久| 人人妻人人爽人人添夜夜欢视频 | 男女无遮挡免费网站观看| 大陆偷拍与自拍| 国产亚洲最大av| 直男gayav资源| 亚洲精品自拍成人| 热99国产精品久久久久久7| 亚洲国产日韩一区二区| 国产亚洲午夜精品一区二区久久| 国产日韩欧美亚洲二区| 国产真实伦视频高清在线观看| 国语对白做爰xxxⅹ性视频网站| 亚洲无线观看免费| 国产亚洲午夜精品一区二区久久| 精品亚洲成a人片在线观看 | 国产亚洲一区二区精品| 国产高清三级在线| 男女边摸边吃奶| 一区二区三区乱码不卡18| 亚洲性久久影院| 男女边吃奶边做爰视频| 黄色欧美视频在线观看| 久久久成人免费电影| 亚洲欧美日韩另类电影网站 | 免费大片18禁| av视频免费观看在线观看| 亚洲综合色惰| 老熟女久久久| 男女下面进入的视频免费午夜| 日韩一区二区三区影片| 制服丝袜香蕉在线| 国语对白做爰xxxⅹ性视频网站| 亚洲在久久综合| 国产男女内射视频| 欧美一级a爱片免费观看看| 一区在线观看完整版| 日韩欧美精品免费久久| 天堂8中文在线网| 少妇的逼好多水| 啦啦啦视频在线资源免费观看| 啦啦啦啦在线视频资源| 久久av网站| 亚洲天堂av无毛| 国产老妇伦熟女老妇高清| 黑丝袜美女国产一区| 国产男人的电影天堂91| 国产精品麻豆人妻色哟哟久久| 欧美三级亚洲精品| 大片免费播放器 马上看| 日日啪夜夜爽| 亚洲精品一二三| av免费观看日本| 久久青草综合色| 国产成人精品一,二区| 日本av免费视频播放| 能在线免费看毛片的网站| 在线观看人妻少妇| 日韩制服骚丝袜av| 日韩一区二区三区影片| 午夜福利高清视频| 99久久综合免费| 亚洲欧美日韩东京热| 我的女老师完整版在线观看| 91精品国产九色| 黑丝袜美女国产一区| 久久久精品94久久精品| 女的被弄到高潮叫床怎么办| av天堂中文字幕网| 免费不卡的大黄色大毛片视频在线观看| 97精品久久久久久久久久精品| 亚洲av日韩在线播放| 免费大片18禁| 另类亚洲欧美激情| www.av在线官网国产| 97超视频在线观看视频| 一二三四中文在线观看免费高清| 女性被躁到高潮视频| 亚洲欧洲日产国产| 丰满少妇做爰视频| 精品久久久精品久久久| 香蕉精品网在线| 欧美人与善性xxx| 日韩强制内射视频| 精品亚洲成国产av| 亚洲伊人久久精品综合| 简卡轻食公司| 18+在线观看网站| 午夜免费观看性视频| 国产在视频线精品| 热re99久久精品国产66热6| 欧美丝袜亚洲另类| 亚洲aⅴ乱码一区二区在线播放| 韩国高清视频一区二区三区| 尾随美女入室| 亚洲av中文av极速乱| 国产一区亚洲一区在线观看| 校园人妻丝袜中文字幕| 精品一区二区免费观看| 卡戴珊不雅视频在线播放| 成人亚洲欧美一区二区av| 亚洲成色77777| 99精国产麻豆久久婷婷| 免费观看av网站的网址| 久久久久久九九精品二区国产| 国产伦理片在线播放av一区| 嘟嘟电影网在线观看| 久久精品国产鲁丝片午夜精品| 久久久久久久大尺度免费视频| 国产欧美日韩一区二区三区在线 | av免费在线看不卡| 国产极品天堂在线| 观看av在线不卡| 色婷婷av一区二区三区视频| 26uuu在线亚洲综合色| 一级片'在线观看视频| 亚洲欧美精品自产自拍| 亚洲国产毛片av蜜桃av| 黑丝袜美女国产一区| 搡女人真爽免费视频火全软件| 久久国内精品自在自线图片| 欧美丝袜亚洲另类| 欧美xxxx性猛交bbbb| 一本一本综合久久| 日韩 亚洲 欧美在线| 天天躁日日操中文字幕| 国产精品熟女久久久久浪| av在线播放精品| 久久人人爽人人片av| 日本欧美国产在线视频| 国产 一区 欧美 日韩| 秋霞伦理黄片| 亚洲无线观看免费| 身体一侧抽搐| 亚洲精品久久久久久婷婷小说| 三级国产精品欧美在线观看| 国产极品天堂在线| 欧美国产精品一级二级三级 | 日本猛色少妇xxxxx猛交久久| 91精品国产国语对白视频| 精品久久久精品久久久| 王馨瑶露胸无遮挡在线观看| 女性被躁到高潮视频| 香蕉精品网在线| 国产大屁股一区二区在线视频| 久久久精品免费免费高清| 国产男女超爽视频在线观看| 内射极品少妇av片p| 蜜臀久久99精品久久宅男| 婷婷色综合大香蕉| 日韩大片免费观看网站| 我要看日韩黄色一级片| 国模一区二区三区四区视频| 麻豆成人午夜福利视频| 亚洲精品国产av蜜桃| 午夜免费男女啪啪视频观看| 大陆偷拍与自拍| 国产精品福利在线免费观看| 精品久久国产蜜桃| 青青草视频在线视频观看| 免费观看的影片在线观看| 男人和女人高潮做爰伦理| 一级av片app| 在线观看一区二区三区激情| 日本-黄色视频高清免费观看| 久久久久久九九精品二区国产| 免费高清在线观看视频在线观看| 亚洲国产av新网站| 九九久久精品国产亚洲av麻豆| 亚洲欧洲国产日韩| 网址你懂的国产日韩在线| 夜夜看夜夜爽夜夜摸| 高清午夜精品一区二区三区| 国产免费一级a男人的天堂| 国产黄频视频在线观看| 亚洲美女黄色视频免费看| 搡女人真爽免费视频火全软件| 国产精品不卡视频一区二区| 国产成人精品一,二区| 夫妻午夜视频| 纵有疾风起免费观看全集完整版| 久久久久久久久久成人| 丝袜喷水一区| 国产精品.久久久| 人人妻人人爽人人添夜夜欢视频 | 欧美3d第一页| 麻豆乱淫一区二区| 中文字幕久久专区| 美女脱内裤让男人舔精品视频| 少妇人妻一区二区三区视频| 一级片'在线观看视频| 一级毛片我不卡| 国产精品一二三区在线看| 高清视频免费观看一区二区| 欧美日韩在线观看h| 一区在线观看完整版| 丰满乱子伦码专区| 日韩免费高清中文字幕av| 精品99又大又爽又粗少妇毛片| 精品亚洲成国产av| av在线app专区| 日日摸夜夜添夜夜添av毛片| 精品国产一区二区三区久久久樱花 | 一本久久精品| 亚洲va在线va天堂va国产| 国产爽快片一区二区三区| 国产伦在线观看视频一区| 网址你懂的国产日韩在线| 国产精品国产三级专区第一集| 亚洲成人av在线免费| 日韩av在线免费看完整版不卡| 国产伦精品一区二区三区视频9| 久久久久久久久大av| 最近最新中文字幕免费大全7| 免费久久久久久久精品成人欧美视频 | 精品亚洲乱码少妇综合久久| 91精品一卡2卡3卡4卡| 亚洲精品aⅴ在线观看| 欧美日韩综合久久久久久| 黄色一级大片看看| 高清黄色对白视频在线免费看 | 亚洲精品成人av观看孕妇| 精品人妻熟女av久视频| 亚洲综合精品二区| 男人狂女人下面高潮的视频| 日本黄色片子视频| 精品亚洲乱码少妇综合久久| 国产精品不卡视频一区二区| 这个男人来自地球电影免费观看 | av线在线观看网站| 精品少妇黑人巨大在线播放| 女人久久www免费人成看片| 免费观看在线日韩| 91狼人影院| 久久亚洲国产成人精品v| 天天躁夜夜躁狠狠久久av| 2018国产大陆天天弄谢| 亚洲精品国产av成人精品| 久久精品国产a三级三级三级| 99re6热这里在线精品视频| 久久精品国产自在天天线| 成人18禁高潮啪啪吃奶动态图 | 久久久久久人妻| 香蕉精品网在线| 80岁老熟妇乱子伦牲交| 国产精品一区二区性色av| 一本—道久久a久久精品蜜桃钙片| 少妇人妻一区二区三区视频| 在线观看美女被高潮喷水网站| 国产 一区 欧美 日韩| 自拍偷自拍亚洲精品老妇| 三级国产精品欧美在线观看| 亚洲精品亚洲一区二区| 99久久精品热视频| 最近中文字幕高清免费大全6| 精品视频人人做人人爽| 夜夜看夜夜爽夜夜摸| 亚州av有码| 成人漫画全彩无遮挡| 久久午夜福利片| 高清在线视频一区二区三区| 在线天堂最新版资源| 久久精品久久久久久噜噜老黄| 欧美日韩国产mv在线观看视频 | 涩涩av久久男人的天堂| 成人黄色视频免费在线看| 久久6这里有精品| 又爽又黄a免费视频| 制服丝袜香蕉在线| 国国产精品蜜臀av免费| 国产老妇伦熟女老妇高清| 午夜免费观看性视频| 久久99蜜桃精品久久| 久久久久久久国产电影| 日韩 亚洲 欧美在线| 亚洲,欧美,日韩| 精品国产乱码久久久久久小说| 一级毛片我不卡| 春色校园在线视频观看| 亚洲国产成人一精品久久久| 熟女电影av网| 十分钟在线观看高清视频www | 日韩大片免费观看网站| 国产精品一区www在线观看| 亚洲,欧美,日韩| 又大又黄又爽视频免费| 能在线免费看毛片的网站| 国产亚洲精品久久久com| 一级爰片在线观看| av线在线观看网站| 人妻 亚洲 视频| 成人毛片60女人毛片免费| 超碰av人人做人人爽久久| 亚洲精品一区蜜桃| 精品人妻偷拍中文字幕| 午夜免费鲁丝| a级毛色黄片| 国产精品.久久久| 哪个播放器可以免费观看大片| 97超视频在线观看视频| 女性被躁到高潮视频| 又爽又黄a免费视频| 国产 一区精品| 亚洲欧美清纯卡通| 精品视频人人做人人爽| 国产成人freesex在线| 亚洲av二区三区四区| 国产黄频视频在线观看| 亚洲国产精品专区欧美| 精品人妻熟女av久视频| 国产成人freesex在线| 男人舔奶头视频| 夫妻午夜视频| 久久国产精品男人的天堂亚洲 | 国产精品熟女久久久久浪| kizo精华| 国产白丝娇喘喷水9色精品| 少妇的逼好多水| videossex国产| 在线播放无遮挡| 久久人妻熟女aⅴ| 少妇人妻久久综合中文| 一本—道久久a久久精品蜜桃钙片| 一级片'在线观看视频| 一级毛片电影观看| 欧美日韩亚洲高清精品| 黄色一级大片看看| 欧美另类一区| 韩国高清视频一区二区三区| www.av在线官网国产| 十分钟在线观看高清视频www | av国产久精品久网站免费入址| 精品一品国产午夜福利视频| 黄片wwwwww| 午夜福利影视在线免费观看| 婷婷色av中文字幕| 亚洲真实伦在线观看| 女性被躁到高潮视频| h日本视频在线播放| av又黄又爽大尺度在线免费看| 少妇丰满av| 亚洲精品国产色婷婷电影| 大香蕉97超碰在线| 日韩视频在线欧美| 国产 一区精品| 成年人午夜在线观看视频| 国产成人精品福利久久| 搡女人真爽免费视频火全软件| 韩国av在线不卡| av在线蜜桃|