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

    運用分頻段希爾伯特黃變換進行多分量信號的頻散分析?

    2012-07-01 18:03:39蔣禮
    電訊技術 2012年4期
    關鍵詞:希爾伯特時頻頻段

    蔣禮

    (1.中國地質(zhì)大學(武漢)地球物理與空間信息學院,武漢430074;2.華北水利水電學院數(shù)學與信息科學學院,鄭州450011)

    運用分頻段希爾伯特黃變換進行多分量信號的頻散分析?

    蔣禮1,2

    (1.中國地質(zhì)大學(武漢)地球物理與空間信息學院,武漢430074;2.華北水利水電學院數(shù)學與信息科學學院,鄭州450011)

    將分頻段希爾伯特黃變換應用于多分量信號的頻散分析中。首先,利用帶通濾波器和經(jīng)驗模態(tài)分解相結(jié)合,成功實現(xiàn)了經(jīng)驗模態(tài)頻率分解,并準確提取了經(jīng)驗頻率模態(tài)函數(shù);然后,使用該方法準確地獲取了多分量含噪信號的時頻能量譜和時頻相位譜;最后,基于同步相差和異步相差算法,精確繪制了原信號的相速度頻散分析曲線。數(shù)值試驗表明,該算法擁有較高的時頻分辨能力和良好的抗噪性能,對于復雜且信噪比較低的信號,能獲得比傳統(tǒng)希爾伯特黃變換更準確的頻散分析結(jié)果。

    希爾伯特黃變換;經(jīng)驗模態(tài)頻率分解;經(jīng)驗頻率模態(tài)函數(shù);瞬時頻率;頻散曲線

    1 引言

    基于瞬時頻率的希爾伯特黃變換(Hilbert-Huang Transform)是一種進行時頻分析的有效方法[1],因為該方法可以準確刻畫信號的瞬時頻率和瞬時相位,所以它成為進行多分量信號頻散分析的有效工具[2-3]。但由于傳統(tǒng)希爾伯特黃變換的核心是基于時間序列的經(jīng)驗模態(tài)分解,該方法在擁有較高的時間分辨能力的同時往往頻率分辨能力較差[4]。雖然對于普通的多分量信號,使用該方法能夠取得令人滿意的結(jié)果,但是對于所含頻率分量較多的情況,使用該方法往往不能令人滿意。

    本文對傳統(tǒng)的希爾伯特黃變換進行改進,運用經(jīng)驗模態(tài)頻率分解理論[5]代替?zhèn)鹘y(tǒng)的經(jīng)驗模態(tài)分解,從而提出分頻段希爾伯特黃變換。該方法的頻率分辨能力較傳統(tǒng)的希爾伯特黃變換有明顯提高。將該方法用于進行多分量信號的頻散分析,可以得到令人滿意的結(jié)果。

    2 原理簡介

    傳統(tǒng)希爾伯特黃變換的核心為經(jīng)驗模態(tài)分解(EMD)和固有模態(tài)函數(shù)(IMF)。由于經(jīng)驗模態(tài)分解進行的是一種時域濾波,所以時間分辨率通常是令人滿意的;但由于分解后的固有模態(tài)函數(shù)個數(shù)有限,故頻率分辨率往往較低。對于頻率成分較為豐富的信號直接進行希爾伯特黃變換通常得不到理想的結(jié)果。

    2.1 經(jīng)驗模態(tài)頻率分解

    分頻段希爾伯特黃變換可以有效解決上述問題。分頻段希爾伯特黃變換的核心為經(jīng)驗模態(tài)頻率分解(EMFD)[5],它是由北京大學程乾生教授于2005

    年提出的一種時頻分解理論,初始目的是為了給希爾伯特黃變換提供一組正交基。但在實際應用的過程中,經(jīng)此方法分解所得的經(jīng)驗頻率模態(tài)函數(shù)(EFMF)不僅正交,而且也具有比固有模態(tài)函數(shù)更高的頻率分辨能力。

    假設某一信號s(t)存在K個分量,通過經(jīng)驗模態(tài)分解后可得

    式中,每個ci(t)稱為一個固有模態(tài)函數(shù)。設每個固有模態(tài)函數(shù)的頻譜為Ci(f),由于希爾伯特黃變換的頻域分辨率較差,雖然Ci(f)和Ci+1(f)的主頻并不相同,但兩者在頻譜卻有不少重疊的部分[6-7](這也正是傳統(tǒng)的固有模態(tài)函數(shù)不是正交基的原因)。假設兩者的主頻分別位于fi的兩側(cè),則設定一個理想帶通濾波器:Hi(f)=1,fi≤f≤fi-1,其余頻段Hi(f)=0,將該濾波器作用于原信號,可得

    濾波所得信號的時間序列di(t)稱為一個經(jīng)驗頻率模態(tài)函數(shù)。上述的整個過程稱為一個經(jīng)驗模態(tài)頻率分解過程。對該多分量信號通過完整的經(jīng)驗模態(tài)頻率分解可得

    此時不僅由于分解所得的di(t)頻譜部分沒有重疊,從而構(gòu)成一組正交基;而且因為分解過程中K個帶通濾波器的作用,所以可以準確無誤地提取每個頻率分量,故該方案具有極高的頻率分辨能力。

    2.2 時頻分析過程

    經(jīng)驗模態(tài)頻率分解的數(shù)學原理雖然十分簡單,但是按照上述的推導流程來實現(xiàn)整個分解過程將會存在下列難點:

    (1)對于一個未知的信號,由于無法預知其頻率分量,從而對于濾波器的數(shù)目及參數(shù)難以確定;

    (2)雖然可以在進行經(jīng)驗模態(tài)頻率分解前先運用經(jīng)驗模態(tài)分解來確定濾波器的數(shù)目及參數(shù),但這樣會浪費分解時間,降低分解效率。

    為了高效實現(xiàn)經(jīng)驗模態(tài)頻率分解過程,可首先在整個頻段上設置一系列的頻點,將相鄰兩個頻點組成帶通濾波器;然后將待分析的信號通過每個帶通濾波器可以得到一系列的頻段信號;最后對每個分頻段信號進行經(jīng)驗模態(tài)分解,可以得到每個頻段的固有模態(tài)函數(shù)。這時每個頻段的固有模態(tài)函數(shù)實際上就是經(jīng)驗頻率模態(tài)函數(shù)(每個經(jīng)驗頻率模態(tài)函數(shù)的頻譜不會有重疊的部分),而對于分頻段信號的經(jīng)驗模態(tài)分解過程實際上就是經(jīng)驗模態(tài)頻率分解過程。

    完整的分頻段希爾伯特黃變換實現(xiàn)流程如圖1所示。與數(shù)學推導流程的不同之處在于,數(shù)學推導是先進行經(jīng)驗模態(tài)分解,而后再制定帶通濾波器濾波;而此處對于未知信號是先做帶通濾波再做經(jīng)驗模態(tài)分解,但兩者最終所得結(jié)果完全一樣(K個固有模態(tài)函數(shù)即為原信號的K個頻率分量,而且頻譜部分不會重疊完全正交)。

    圖1 分頻段希爾伯特黃變換算法流程Fig.1 The algorithm flow of sub-band Hilbert-Huang transform

    2.3 頻散分析過程

    對原始多分量信號進行分頻段希爾伯特黃變換可以得到時頻能量譜和時頻相位譜[3];對上述兩譜進行分析,可以獲得任意分量信號的瞬時頻率、瞬時相位、瞬時能量、持續(xù)時間等信息[3]。根據(jù)已獲知的信息運用同步相差和異步相差理論可以求取任意頻率分量的相速度值[3]。完整的頻散分析過程如圖2所示。

    圖2 運用分頻段希爾伯特黃變換進行頻散分析流程Fig.2 The algorithm flow of dispersion analysis by sub-band Hilbert-Huang transform

    3 仿真實例

    為了驗證分頻段希爾伯特黃變換的正確性以及較傳統(tǒng)希爾伯特黃變換的優(yōu)越性,此處設定一個較為復雜的多分量信號模型。

    3.1 模型建立

    假設瑞利波源由10個分量信號疊加而成,初始信號持續(xù)時間為0.2 s;各分量信號的群速度和相速度隨著頻率的增加而減小,且群速度小于相速度;信號能量集中在50 Hz、60 Hz,隨著頻率的增加和減少,能量依次遞減。各分量信號的具體參數(shù)詳見表1。

    表1 瑞利波源分量參數(shù)列表Table 1 Rayleigh wave source component parameter list

    兩個檢波器分別距離震源12 m和18m,檢波時間為1.023 s,采樣時間間隔為1ms。為了檢驗該算法的抗噪性能,檢波信號除了添加10%的高斯白噪聲模擬檢波干擾外,再整體添加5%的高斯白噪聲作為環(huán)境噪聲。12m處獲取的信號如圖3所示,18m處獲取的信號如圖4所示。

    圖3 信號1Fig.3 Signal1

    圖4 信號2Fig.4 Signal2

    3.2 譜分析

    3.2.1 希爾伯特黃變換

    運用傳統(tǒng)希爾伯特黃變換所得譜分析結(jié)果如圖5所示。圖5(a)和(b)分別為12m和18m處信號的時間-頻率-能量譜;而圖5(c)和(d)則分別為12m和18m處信號的時間-頻率-相位譜。從上述各圖可以發(fā)現(xiàn),無論是時頻能量譜還是時頻相位譜都不能準確刻畫出10個頻率分量的完整信息。雖然由于希爾伯特黃變換良好的時間分辨能力,對比12m處和18m處從高到低的各頻率分量能量譜可以清晰地看出頻散效應;但是由于不能將10個頻率分量的時頻信息完整地反映出來,故而無法求取每個分量的準確相速度值。

    圖5 傳統(tǒng)希爾伯特黃變換的譜分析圖Fig.5 Spectrum analysis diagram of Hilbert-Huang transform

    圖6 分頻段希爾伯特黃變換的譜分析圖Fig.6 Spectrum analysis diagram of the sub-band Hilbert-Huang transform

    3.2.2 分頻段希爾伯特黃變換

    分頻段帶通濾波器必須預設兩個參數(shù):有效頻段和濾波器帶寬。有效頻段反映的是波源信號的物理特性,而帶寬則反映了頻率分辨能力。理論情況下,無窮大的有效頻段和無限小的帶寬可以滿足任何多分量信號的分解;但在實際應用中,有效頻段設定越小,帶寬設定越大,則運算效率越高。此處有效頻段設為5~105 Hz,以10 Hz為帶寬進行窄帶濾波。

    運用分頻段希爾伯特黃變換所得譜分析結(jié)果如圖6所示。從上述各圖可以發(fā)現(xiàn),無論是在時頻能量譜還是時頻相位譜中,10個頻率分量都顯示得十分清晰,而且可以看出各分量的瞬時頻率非常接近于真實的物理頻率,延遲時間、持續(xù)時間和能量分布也基本正確。

    從圖6(a)和(b)中可以看出12m處及18m處的各分量信號均通過0.2 s;18m處的信號各分量均通過0.25 s。為了分析方便,只測量這兩個時刻的瞬時頻率和瞬時相位,所測量的數(shù)據(jù)如表2所示(表2中的平均頻率為某一分量所有時刻瞬時頻率的加權(quán)平均[3])。從表2中可以看出,無論是瞬時頻率還是平均頻率都非常接近接近于該分量信號的真實物理頻率。

    表2 實測時間-頻率-相位Table 2 Measured values of time-frequency-phase

    3.3 計算結(jié)果分析

    對表2所測得的數(shù)據(jù)進行相速度計算,所得結(jié)果及誤差如表3所示。從計算結(jié)果來看,所得相速度較為理想。第三分量至第十分量基本接近于相應模型分量的真實速度;第一分量和第二分量雖然誤差較大,但仍在可接受的范圍之內(nèi)。

    表3 相速度及誤差計算結(jié)果Table 3 Results of phase velocity and error

    上述兩個方法出現(xiàn)的誤差主要產(chǎn)生于分母部分的相差和頻率的比值項Δφ/ω0[3]。理論上Δφ/ω0的誤差會隨著相位解纏時頻率的升高而逐漸減小。

    3.4 頻散曲線

    根據(jù)表3的計算結(jié)果可以繪制4條頻散曲線。若在位于24m的測點處增加一個檢波器記錄,則可以繪制12條頻散曲線。其中8條是道間距為6m的頻散曲線,4條為道間距為12m的頻散曲線。所有的相速度頻散曲線如圖7所示。

    圖7 相速度頻散曲線Fig.7 Phase velocity dispersion curves

    從圖7中可以看出,頻散曲線正確反應了頻散模型的相速度變化規(guī)律。隨著頻率的升高,頻散曲線的收斂性越來越好;而且道間距為12m的頻散曲線誤差要整體小于道間距為6m的頻散曲線誤差。

    4 結(jié)論

    經(jīng)驗模態(tài)頻率分解在保留經(jīng)驗模態(tài)分解較高的時間分辨能力同時,又提高了頻率分辨能力,以其為基礎的分頻段希爾伯特黃變換具有極為準確的時頻-相位定位能力。將該方法用于進行多分量信號的頻散分析,不僅計算結(jié)果準確,而且具有較強的抗噪能力。

    使用該方法時有兩個細節(jié)需要引起注意。

    (1)進行分頻段處理時選擇合適的頻窗,不同窗型及窗口的大小會影響最終的結(jié)果。因為進行窄帶濾波必然會降低希爾伯特譜的時間分辨率。此時應注意,頻窗不應取得太小,窗型不應選擇劇烈跳變類型。

    (2)如果對源信號各分量的頻率成分較為熟悉,則可采用頻段交叉的頻窗進行窄帶濾波,即在不減小窗口大小的情況下使窗口頻段交叉(頻窗的主頻部分不可交叉)。這樣可以在相同有效頻段的基礎上增加帶通濾波器的數(shù)量,即能夠在不降低時間分辨率的同時進一步提高頻率分辨率。

    運用譜分析方法計算多分量信號的頻散特性曲線關鍵在于兩點:準確的時頻分析和精確的速度計算。從本文和文獻[3]的對比可以看出,S-HHT與傳統(tǒng)的HHT相比雖然增加了計算量,但提高了頻域分辨能力及抗噪性能。當待分析信號較為復雜或信噪比較差時,使用S-HHT能夠獲得更準確的時頻譜。而同步相差和異步相差算法,在進行相速度計算時精度相似,都是進行頻散分析的優(yōu)秀算法。

    [1]Huang N E,Shen Zheng,Long SR.The empiricalmode decomposition and the Hilbertspectrum for nonlinear and nonstationary time series analysis[J].Proceedings of the Royal Society of London,Series A-Mathematical Physical and Engineering Sciences,1998,454(1971):903-995.

    [2]Chen Chau-Huei,Li Cheng-Pling,Teng Ta-Liang.Surface-Wave Dispersion Measurements Using Hilbert-Huang Transform[J].Terrestrial Atmospheric and Oceanic Sciences,2002,13(2):171-184.

    [3]蔣禮.運用希爾伯特黃變換進行多分量信號的頻散分析[J].電訊技術,2011,51(7):60-66. JIANG Li.Frequency Dispersion Analysis of Multi-component Signalusing Hilbert-Huang Transform[J].Telecommunication Engineering,2011,51(7):60-66.(in Chinese)

    [4]Patrick F,GabrielR,Paulo G.EmpiricalMode Decomposition as a Filter Bank[J].IEEE Signal Processing Letters,2004,11(2):112-114.

    [5]程乾生,武連文.時間序列的經(jīng)驗模態(tài)頻率分解EMFD[J].數(shù)學的實踐與認識,2005,36(5):151-153. CHEN Qian-sheng,WU Lian-wen.The EmpiricalMode Frequency Decomposition for Time Series Analysis[J].Mathematics in Practice and Theory,2005,36(5):151-153.(in Chinese)

    [6]Huang N E,Shen SSP.Hilbert-Huang Transform and Its Application,Interdisciplinary Mathematical Sciences vol 5[M].London:World Scientific Public Co Information,2005:57-74.

    [7]WU Zhaohua,Huang N E.A Study of The Characteristics of White Noise Using The Empirical Mode Decomposition Method[J].Proceedings of the Royal Society London,Series A-Mathematical Physical and Engineering Sciences,2004,460(2046):1597-1611.

    Analysis of Frequency Dispersion of M ulti-com ponent Signal using Sub-band Hilbert-Huang Transform

    JIANG Li1,2
    (1.Institute of Geophysics and Geomatics,China University of Geosciences(Wuhan),Wuhan 430074,China;2.College of Mathematics and Information Science,North China University ofWater Resources and Electric Power,Zhengzhou 450011,China)

    In this paper,the sub-band Hilbert-Huang transform(S-HHT)is used to analyse the frequency dispersion ofmulti-component signal.Firstly,the band-pass filter and the empiricalmode decomposition are combined to achieve the empiricalmode frequency decomposition(EMFD)successfully,and accurately obtain the empirical frequencymode functions(EFMF).Secondly,the time frequency energy spectrums and time frequency phase spectrums from multi-componentsignalwith noise are precisely drawn by themethod.Finally,by using the synchronous phase difference and asynchronous phase difference arithmetic the phase velocity dispersion curves are drawn accurately.The numerical results show that the S-HHTarithmetic hashigh time frequency resolution and good resistance to noise.For complex and low Signal-to-Noise Ratio(SNR)signal,the SHHT can get themore accurate results of dispersion analysis than the HHT.

    Hilbert-Huang transform;empirical mode frequency decomposition;empirical frequency mode function;instant frequency;dispersion curves

    The National Natural Science Foundation of China(No.40974079)

    the B.S.degree and the M.S.degree in 2001 and 2004,respectively.He is now a lecturer and currently working toward the Ph.D.degree.His research concerns geoscience-signal analysis and processing

    1001-893X(2012)04-0472-06

    2011-12-21;

    2012-03-09

    國家自然科學基金資助項目(40974079)

    TN911

    A

    10.3969/j.issn.1001-893x.2012.04.010

    蔣禮(1979—),男,江蘇江陰人,2001年獲學士學位,2004年獲碩士學位,現(xiàn)為博士研究生,講師,主要從事地學信號的分析及處理。

    Email:jlncwu@126.com

    JIANG Li was born in Jiangyin,Jiangsu Province,in 1979. He

    猜你喜歡
    希爾伯特時頻頻段
    一個真值函項偶然邏輯的希爾伯特演算系統(tǒng)
    邏輯學研究(2021年3期)2021-09-29 06:54:34
    gPhone重力儀的面波頻段響應實測研究
    地震研究(2021年1期)2021-04-13 01:04:56
    下一個程序是睡覺——數(shù)學家希爾伯特的故事
    基于希爾伯特-黃變換和小波變換的500kV變電站諧振數(shù)據(jù)對比分析
    電測與儀表(2016年7期)2016-04-12 00:22:14
    推擠的5GHz頻段
    CHIP新電腦(2016年3期)2016-03-10 14:07:52
    基于希爾伯特- 黃變換的去噪法在外測數(shù)據(jù)處理中的應用
    TD—LTE在D頻段和F頻段的覆蓋能力差異
    中國新通信(2015年1期)2015-05-30 10:30:46
    基于時頻分析的逆合成孔徑雷達成像技術
    對采樣數(shù)據(jù)序列進行時頻分解法的改進
    雙線性時頻分布交叉項提取及損傷識別應用
    久久久久网色| 你懂的网址亚洲精品在线观看| 国产精品欧美亚洲77777| 丝袜喷水一区| 26uuu在线亚洲综合色| 99热全是精品| 黄片无遮挡物在线观看| 国产黄片视频在线免费观看| 国产成人freesex在线| 久久97久久精品| 日日摸夜夜添夜夜添av毛片| 国产精品免费大片| 赤兔流量卡办理| 青春草亚洲视频在线观看| 成人午夜精彩视频在线观看| 欧美高清性xxxxhd video| 老熟女久久久| 亚洲av综合色区一区| 国产色爽女视频免费观看| 国产欧美另类精品又又久久亚洲欧美| 精品亚洲乱码少妇综合久久| 国产精品一及| 99热这里只有是精品50| 又粗又硬又长又爽又黄的视频| h视频一区二区三区| 丰满乱子伦码专区| 中文资源天堂在线| 日韩成人av中文字幕在线观看| 国产探花极品一区二区| 日本爱情动作片www.在线观看| 十八禁网站网址无遮挡 | 亚洲精品日韩av片在线观看| 国产欧美日韩一区二区三区在线 | 最后的刺客免费高清国语| 亚洲国产高清在线一区二区三| 中文乱码字字幕精品一区二区三区| 国产永久视频网站| 免费黄色在线免费观看| 蜜臀久久99精品久久宅男| 香蕉精品网在线| 国产 一区精品| 亚洲不卡免费看| xxx大片免费视频| 国产在线男女| 亚洲精品久久午夜乱码| h日本视频在线播放| 国产亚洲精品久久久com| 交换朋友夫妻互换小说| 亚洲欧美一区二区三区国产| 黄片无遮挡物在线观看| 99久久人妻综合| 欧美精品国产亚洲| 麻豆成人av视频| 嫩草影院入口| 国产女主播在线喷水免费视频网站| 国产中年淑女户外野战色| 精品人妻视频免费看| 亚洲精品一区蜜桃| 免费久久久久久久精品成人欧美视频 | 国产成人午夜福利电影在线观看| 在线观看美女被高潮喷水网站| 人人妻人人澡人人爽人人夜夜| 国产精品嫩草影院av在线观看| 亚洲av国产av综合av卡| 狠狠精品人妻久久久久久综合| 18禁裸乳无遮挡动漫免费视频| 欧美高清成人免费视频www| 国产精品久久久久久av不卡| 久久久久国产精品人妻一区二区| 老熟女久久久| 中文资源天堂在线| 水蜜桃什么品种好| 一区在线观看完整版| 国产高清三级在线| 精品99又大又爽又粗少妇毛片| 国产成人a区在线观看| 久热久热在线精品观看| 在现免费观看毛片| av在线播放精品| 一区二区三区免费毛片| 蜜桃亚洲精品一区二区三区| 男女无遮挡免费网站观看| 欧美变态另类bdsm刘玥| 国模一区二区三区四区视频| 国产精品99久久99久久久不卡 | 国产真实伦视频高清在线观看| a级毛片免费高清观看在线播放| 少妇的逼好多水| 最近中文字幕高清免费大全6| 一边亲一边摸免费视频| 男人爽女人下面视频在线观看| 成人二区视频| 午夜视频国产福利| 久久午夜福利片| 91午夜精品亚洲一区二区三区| 久久精品国产a三级三级三级| 777米奇影视久久| 日本黄色日本黄色录像| 国产成人aa在线观看| 国产伦精品一区二区三区视频9| 国产淫语在线视频| 天天躁日日操中文字幕| freevideosex欧美| 国产一区有黄有色的免费视频| 天美传媒精品一区二区| 成人毛片60女人毛片免费| 欧美高清性xxxxhd video| 中国国产av一级| 最近中文字幕2019免费版| 91精品国产九色| 亚洲欧洲日产国产| 国产有黄有色有爽视频| 欧美成人精品欧美一级黄| 色哟哟·www| 你懂的网址亚洲精品在线观看| 中文字幕亚洲精品专区| 久久久久人妻精品一区果冻| 日本-黄色视频高清免费观看| 久热久热在线精品观看| 国产高清国产精品国产三级 | 日本av免费视频播放| 久久久精品免费免费高清| 丰满少妇做爰视频| 久久鲁丝午夜福利片| 日日摸夜夜添夜夜爱| 中文字幕av成人在线电影| 国产高清三级在线| 国产免费一区二区三区四区乱码| 亚洲美女黄色视频免费看| 国产在线男女| 亚洲国产精品成人久久小说| 我的女老师完整版在线观看| 大香蕉久久网| 久久精品久久久久久久性| 国产亚洲91精品色在线| 国产熟女欧美一区二区| 亚洲美女搞黄在线观看| 最近中文字幕高清免费大全6| 2022亚洲国产成人精品| 久久婷婷青草| 国产精品嫩草影院av在线观看| 久久久成人免费电影| 亚洲最大成人中文| 色视频www国产| 国产精品女同一区二区软件| 国产永久视频网站| 久久婷婷青草| 乱系列少妇在线播放| 国产成人精品婷婷| 亚洲电影在线观看av| 日日啪夜夜爽| 99久久综合免费| 日韩不卡一区二区三区视频在线| 亚洲最大成人中文| 欧美成人精品欧美一级黄| 日本爱情动作片www.在线观看| 亚洲欧美精品专区久久| kizo精华| 日本黄色片子视频| 搡老乐熟女国产| 免费人妻精品一区二区三区视频| 亚洲人成网站在线观看播放| 国产黄色视频一区二区在线观看| 久久人人爽人人爽人人片va| 大话2 男鬼变身卡| 男人添女人高潮全过程视频| 国产精品国产三级国产专区5o| 最近最新中文字幕大全电影3| 色综合色国产| 黄色一级大片看看| 网址你懂的国产日韩在线| 波野结衣二区三区在线| 夫妻性生交免费视频一级片| 日本-黄色视频高清免费观看| 在线观看免费视频网站a站| 午夜激情福利司机影院| a级毛片免费高清观看在线播放| 久久99蜜桃精品久久| 国产在线视频一区二区| 成年av动漫网址| 性色avwww在线观看| 国产欧美另类精品又又久久亚洲欧美| 国产精品国产三级国产专区5o| 卡戴珊不雅视频在线播放| 午夜精品国产一区二区电影| 少妇精品久久久久久久| 精品一区二区三卡| av视频免费观看在线观看| 18+在线观看网站| 人妻一区二区av| 伦理电影大哥的女人| 三级经典国产精品| 成人美女网站在线观看视频| 免费在线观看成人毛片| 久久精品久久精品一区二区三区| 少妇人妻久久综合中文| 大话2 男鬼变身卡| 日韩成人伦理影院| 国产久久久一区二区三区| 春色校园在线视频观看| 国产国拍精品亚洲av在线观看| 男女啪啪激烈高潮av片| 国产探花极品一区二区| 又爽又黄a免费视频| 久久久色成人| 超碰av人人做人人爽久久| 麻豆精品久久久久久蜜桃| 国产精品人妻久久久久久| 91aial.com中文字幕在线观看| 九九在线视频观看精品| 午夜视频国产福利| 午夜精品国产一区二区电影| 美女cb高潮喷水在线观看| 亚洲精品国产色婷婷电影| 少妇丰满av| 男女免费视频国产| 亚洲,欧美,日韩| 新久久久久国产一级毛片| 久久久久国产精品人妻一区二区| 国产高清有码在线观看视频| 免费黄色在线免费观看| 国产成人a∨麻豆精品| 亚洲欧洲日产国产| av国产免费在线观看| 亚洲电影在线观看av| 亚洲精品一二三| 伊人久久国产一区二区| 日韩 亚洲 欧美在线| 国产av国产精品国产| 国产男人的电影天堂91| 国产亚洲5aaaaa淫片| 夫妻性生交免费视频一级片| 99久久精品一区二区三区| 亚洲精品456在线播放app| 国产亚洲精品久久久com| 偷拍熟女少妇极品色| 日韩大片免费观看网站| av播播在线观看一区| 国产av码专区亚洲av| 亚洲欧美成人精品一区二区| 亚洲av在线观看美女高潮| 最近的中文字幕免费完整| 18禁在线无遮挡免费观看视频| 欧美另类一区| 九九久久精品国产亚洲av麻豆| 91午夜精品亚洲一区二区三区| 免费av不卡在线播放| 亚洲精华国产精华液的使用体验| 97在线人人人人妻| 亚洲av综合色区一区| 国产片特级美女逼逼视频| 日日摸夜夜添夜夜爱| 伦理电影免费视频| 国产精品国产三级国产av玫瑰| 亚洲欧美精品专区久久| 寂寞人妻少妇视频99o| 黄色一级大片看看| 另类亚洲欧美激情| tube8黄色片| 纯流量卡能插随身wifi吗| a 毛片基地| 欧美xxxx黑人xx丫x性爽| 青青草视频在线视频观看| 国产男人的电影天堂91| 精品国产三级普通话版| 91午夜精品亚洲一区二区三区| 欧美区成人在线视频| 久久精品国产亚洲av涩爱| 大又大粗又爽又黄少妇毛片口| 免费黄网站久久成人精品| 99热全是精品| 国产片特级美女逼逼视频| 欧美xxxx黑人xx丫x性爽| 黄色怎么调成土黄色| 久久久精品94久久精品| 一本色道久久久久久精品综合| 精华霜和精华液先用哪个| 汤姆久久久久久久影院中文字幕| 永久网站在线| 亚洲成色77777| 国产在线视频一区二区| 国产精品人妻久久久久久| 亚洲三级黄色毛片| 日日啪夜夜爽| 国产欧美日韩一区二区三区在线 | 日韩欧美一区视频在线观看 | 国产黄色视频一区二区在线观看| 亚洲欧美清纯卡通| 少妇猛男粗大的猛烈进出视频| 偷拍熟女少妇极品色| 男女国产视频网站| 少妇猛男粗大的猛烈进出视频| 国产亚洲av片在线观看秒播厂| 九草在线视频观看| 久久久午夜欧美精品| 爱豆传媒免费全集在线观看| 久久青草综合色| 少妇 在线观看| 免费少妇av软件| 夜夜爽夜夜爽视频| 免费少妇av软件| 夜夜爽夜夜爽视频| 男人添女人高潮全过程视频| 黄色欧美视频在线观看| 国产成人精品一,二区| 少妇精品久久久久久久| 国产高清三级在线| 久久久a久久爽久久v久久| 精品亚洲成国产av| 亚洲国产精品国产精品| 精品人妻一区二区三区麻豆| 亚洲不卡免费看| 国产成人一区二区在线| 99热网站在线观看| 国产人妻一区二区三区在| 麻豆成人av视频| 亚洲av中文av极速乱| 亚洲国产欧美在线一区| 能在线免费看毛片的网站| 特大巨黑吊av在线直播| a级一级毛片免费在线观看| 国产av精品麻豆| 一本—道久久a久久精品蜜桃钙片| 国产精品嫩草影院av在线观看| 简卡轻食公司| 国产精品一二三区在线看| 99热6这里只有精品| 一级黄片播放器| 男人狂女人下面高潮的视频| 欧美日韩亚洲高清精品| 久久久久久久大尺度免费视频| 国产爽快片一区二区三区| 亚洲图色成人| 涩涩av久久男人的天堂| 日本欧美视频一区| 成人午夜精彩视频在线观看| 亚洲精品一二三| 人妻制服诱惑在线中文字幕| 成人免费观看视频高清| 人体艺术视频欧美日本| 久久鲁丝午夜福利片| 丝袜喷水一区| 成年女人在线观看亚洲视频| 18+在线观看网站| 亚洲久久久国产精品| 国产黄片美女视频| 欧美xxⅹ黑人| 91午夜精品亚洲一区二区三区| 只有这里有精品99| 少妇被粗大猛烈的视频| 七月丁香在线播放| 精品久久久久久久久av| 青春草亚洲视频在线观看| 国精品久久久久久国模美| 一区二区三区免费毛片| 在线 av 中文字幕| 秋霞在线观看毛片| 女的被弄到高潮叫床怎么办| 精品亚洲乱码少妇综合久久| 国产免费又黄又爽又色| 欧美日本视频| 人妻制服诱惑在线中文字幕| 国产av码专区亚洲av| 午夜福利高清视频| 一区在线观看完整版| 高清日韩中文字幕在线| 免费av中文字幕在线| 亚洲,一卡二卡三卡| 美女cb高潮喷水在线观看| 欧美xxxx性猛交bbbb| 中文字幕亚洲精品专区| 亚洲人成网站高清观看| 免费看不卡的av| 女人久久www免费人成看片| 肉色欧美久久久久久久蜜桃| 久久热精品热| 看非洲黑人一级黄片| 国产亚洲91精品色在线| 天美传媒精品一区二区| 国产在视频线精品| 麻豆乱淫一区二区| 新久久久久国产一级毛片| 久久久欧美国产精品| 成人亚洲精品一区在线观看 | 自拍偷自拍亚洲精品老妇| 国产伦精品一区二区三区视频9| 国产亚洲欧美精品永久| 精品少妇久久久久久888优播| 国产在视频线精品| 性色av一级| 国产成人精品一,二区| 亚洲在久久综合| 亚洲欧美一区二区三区黑人 | 国产视频首页在线观看| 国产精品久久久久久av不卡| 国产成人aa在线观看| 免费黄网站久久成人精品| 人妻 亚洲 视频| 日本黄色日本黄色录像| 肉色欧美久久久久久久蜜桃| 九草在线视频观看| 欧美日韩一区二区视频在线观看视频在线| 国产在线视频一区二区| 亚洲精品aⅴ在线观看| av免费观看日本| 日本猛色少妇xxxxx猛交久久| 久久国产亚洲av麻豆专区| 51国产日韩欧美| av黄色大香蕉| 久久av网站| 永久免费av网站大全| 成人国产麻豆网| 青春草亚洲视频在线观看| 国产亚洲精品久久久com| 少妇熟女欧美另类| 大话2 男鬼变身卡| 大片免费播放器 马上看| 纯流量卡能插随身wifi吗| 国产在线男女| 夜夜看夜夜爽夜夜摸| 观看av在线不卡| 精品一区在线观看国产| 日韩免费高清中文字幕av| 嘟嘟电影网在线观看| 赤兔流量卡办理| 亚洲国产高清在线一区二区三| 一区在线观看完整版| 久久久精品94久久精品| 夜夜骑夜夜射夜夜干| 国产爽快片一区二区三区| 熟女av电影| 高清黄色对白视频在线免费看 | 男女免费视频国产| 亚洲国产av新网站| 日韩成人av中文字幕在线观看| 精品人妻熟女av久视频| 欧美最新免费一区二区三区| 国产精品女同一区二区软件| 久久精品国产鲁丝片午夜精品| 人妻夜夜爽99麻豆av| 国产精品精品国产色婷婷| 中文字幕亚洲精品专区| 国产熟女欧美一区二区| 精品一区在线观看国产| 精品亚洲成国产av| 久久国产精品大桥未久av | 一级av片app| 精品酒店卫生间| 九色成人免费人妻av| 婷婷色av中文字幕| 美女脱内裤让男人舔精品视频| 国产精品三级大全| 如何舔出高潮| 国产成人精品久久久久久| 偷拍熟女少妇极品色| 一级爰片在线观看| 色吧在线观看| 内地一区二区视频在线| 伦理电影大哥的女人| 日韩精品有码人妻一区| 久久6这里有精品| 如何舔出高潮| 欧美激情极品国产一区二区三区 | 九草在线视频观看| 婷婷色综合大香蕉| 久久97久久精品| 中文字幕亚洲精品专区| www.av在线官网国产| 亚洲成人中文字幕在线播放| 日韩成人av中文字幕在线观看| 成年免费大片在线观看| 久久av网站| 国产精品无大码| 日韩成人av中文字幕在线观看| 啦啦啦啦在线视频资源| 欧美日本视频| 久久99蜜桃精品久久| 少妇裸体淫交视频免费看高清| 干丝袜人妻中文字幕| 99国产精品免费福利视频| 中文欧美无线码| 久久青草综合色| 国产欧美日韩一区二区三区在线 | www.色视频.com| 春色校园在线视频观看| 国产日韩欧美在线精品| 成人毛片60女人毛片免费| 有码 亚洲区| 国产精品国产三级国产专区5o| 中文资源天堂在线| 亚洲欧美日韩无卡精品| 啦啦啦中文免费视频观看日本| 亚洲国产成人一精品久久久| 三级经典国产精品| 99久久中文字幕三级久久日本| 国产亚洲精品久久久com| 男人添女人高潮全过程视频| 午夜福利高清视频| 舔av片在线| 蜜桃亚洲精品一区二区三区| 女人久久www免费人成看片| www.色视频.com| 热99国产精品久久久久久7| 日韩一区二区视频免费看| 欧美亚洲 丝袜 人妻 在线| 视频中文字幕在线观看| 欧美高清性xxxxhd video| 九九在线视频观看精品| 中文精品一卡2卡3卡4更新| 狂野欧美白嫩少妇大欣赏| 欧美人与善性xxx| 一区二区av电影网| 国内精品宾馆在线| 国产成人午夜福利电影在线观看| 欧美日韩在线观看h| 18禁在线无遮挡免费观看视频| 一级爰片在线观看| 男人舔奶头视频| 国产高清有码在线观看视频| 高清在线视频一区二区三区| 国产精品一区二区性色av| 噜噜噜噜噜久久久久久91| 欧美日韩亚洲高清精品| 国产熟女欧美一区二区| 男女边吃奶边做爰视频| 亚洲精品456在线播放app| 久久国产精品大桥未久av | 街头女战士在线观看网站| av免费在线看不卡| 黄色视频在线播放观看不卡| 欧美+日韩+精品| 精品久久久久久电影网| 在线免费十八禁| 精品少妇黑人巨大在线播放| 少妇被粗大猛烈的视频| 热99国产精品久久久久久7| 亚洲欧美成人精品一区二区| 最新中文字幕久久久久| 女性被躁到高潮视频| 亚洲无线观看免费| 久久久久久久久大av| 少妇 在线观看| 久久女婷五月综合色啪小说| 国产精品免费大片| 久久久久久久大尺度免费视频| 成年美女黄网站色视频大全免费 | 欧美成人精品欧美一级黄| 亚洲电影在线观看av| 久久综合国产亚洲精品| 国产大屁股一区二区在线视频| 永久网站在线| av国产精品久久久久影院| 亚洲第一av免费看| 日本午夜av视频| 国产淫片久久久久久久久| 久久鲁丝午夜福利片| 2018国产大陆天天弄谢| 纯流量卡能插随身wifi吗| 边亲边吃奶的免费视频| 一二三四中文在线观看免费高清| 一级毛片aaaaaa免费看小| 国产一区二区在线观看日韩| 人妻少妇偷人精品九色| 亚洲经典国产精华液单| 国产亚洲最大av| av线在线观看网站| 欧美亚洲 丝袜 人妻 在线| 国语对白做爰xxxⅹ性视频网站| 日韩成人伦理影院| av又黄又爽大尺度在线免费看| 国产v大片淫在线免费观看| 亚洲精品aⅴ在线观看| 欧美成人a在线观看| 日韩成人伦理影院| 精品午夜福利在线看| 国产爽快片一区二区三区| 亚洲av国产av综合av卡| 美女cb高潮喷水在线观看| 免费看光身美女| a级毛片免费高清观看在线播放| 97在线人人人人妻| 成人黄色视频免费在线看| 丝袜脚勾引网站| 欧美区成人在线视频| 最近中文字幕高清免费大全6| 99久久精品一区二区三区| 国产成人免费观看mmmm| 久热久热在线精品观看| 这个男人来自地球电影免费观看 | 久久精品久久久久久噜噜老黄| 少妇人妻精品综合一区二区| 国产熟女欧美一区二区| 免费观看无遮挡的男女| 亚洲欧美成人精品一区二区| 777米奇影视久久| 尤物成人国产欧美一区二区三区| 成人黄色视频免费在线看| 国产中年淑女户外野战色| 国产精品一区二区性色av| 成人无遮挡网站| 精品酒店卫生间| 日韩成人av中文字幕在线观看| 久久毛片免费看一区二区三区| 777米奇影视久久| 中文在线观看免费www的网站| 最黄视频免费看| 亚洲av综合色区一区| 亚洲精品日韩在线中文字幕| 亚洲国产色片| 街头女战士在线观看网站| 亚洲精品乱码久久久v下载方式| 免费久久久久久久精品成人欧美视频 | av又黄又爽大尺度在线免费看| 99热这里只有是精品50| 九色成人免费人妻av|