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

    基于經驗模分解的陀螺信號消噪

    2011-01-04 07:57:14隋立芬
    測繪學報 2011年6期
    關鍵詞:陀螺小波方差

    甘 雨,隋立芬

    信息工程大學測繪學院,河南鄭州450052

    基于經驗模分解的陀螺信號消噪

    甘 雨,隋立芬

    信息工程大學測繪學院,河南鄭州450052

    陀螺隨機漂移是影響慣性導航精度的重要因素。小波消噪方法對異常噪聲效果不明顯,且對小波基和分解尺度等因素依賴性較強。提出陀螺信號經驗模分解(EMD)消噪方法,將信號進行經驗模分解得到一個本征模態(tài)函數(IMF)組,先基于2σ準則處理異常噪聲IMF分量,再利用相關系數確定高頻噪聲IMF分量個數,將噪聲分量去除以實現(xiàn)陀螺信號消噪。詳細對比小波方法與EMD方法,利用交疊式Allan方差分析兩者的消噪效果,通過慣導算例進一步驗證EMD方法的實效性。結果表明,相比小波方法,EMD消噪法能剔除異常噪聲,可以更有效地抑制陀螺漂移。

    陀螺隨機漂移;小波;經驗模分解;消噪

    1 引 言

    陀螺隨機漂移是影響慣性導航精度的重要因素。抑制陀螺漂移的方法主要有兩種[1]:① 建立漂移模型,使用Kalman濾波等方法進行補償;②對陀螺輸出信號進行消噪處理。由于隨機漂移往往表現(xiàn)為弱非線性、非平穩(wěn)、慢時變,且易受到外部環(huán)境等多種不確定因素的影響[2],無法建立其準確的系統(tǒng)模型,故需要采用陀螺信號消噪的方法。

    目前對陀螺信號消噪主要采用小波方法。小波具有優(yōu)良的多分辨率分析特性,小波消噪不需要系統(tǒng)的誤差模型,因此被廣泛用于陀螺信號的消噪處理中[1-4]。然而,小波分解雖然也能實現(xiàn)對非平穩(wěn)信號的濾波,但其實質是帶通濾波器,限制了濾波的精確性,且小波變換中小波基一經選定,整個信號分析過程中就只能使用這一個小波基,即小波變換是非適應性的[5]。有研究表明[6]:對于如陀螺信號一類的非平穩(wěn)信號,當異常噪聲淹沒了有用信號時,采用小波消噪效果也不甚理想,對這類信號的消噪目前還未發(fā)現(xiàn)較好的方法。

    文獻[7]提出一種分析非平穩(wěn)、非線性信號的自適應分解方法:經驗模分解方法(empirical mode decomposition,EMD),它將復雜的信號分解成若干個按頻率高低排列的本征模態(tài)函數(intrinsic mode function,IMF),每個IMF是一個單分量信號。該方法與小波分析的區(qū)別在于它不需要事先選定基函數,而是根據信號本身的特性自適應地產生合適的模態(tài)函數,這些模態(tài)函數能很好地反映信號在任何時間局部的頻率特征。該方法已經用于機械振動信號分析、SAR影像濾波、氣象以及GPS信號處理等領域[8-11]。

    將EMD方法引入到陀螺信號的處理中,給出陀螺EMD消噪的方法,按照2σ準則剔除振幅偏大的異常噪聲,利用相關系數確定隨機高頻噪聲IMF分量的個數。與小波消噪依賴小波基、分解尺度、閾值估計方法不同,EMD消噪過程完全依賴陀螺信號本身特性。用交疊式Allan方差對比分析本文方法與小波消噪方法,通過慣性導航算例進一步驗證方法的實效性。

    2 EMD消噪方法

    2.1 EMD基本原理

    EMD分解方法認為任何待分解信號都由一組固有振動模式構成,并據此將信號分解為若干本征模態(tài)函數IMF的和。這些本征模態(tài)函數既可以是線性的,也可以是非線性的;既可以是平穩(wěn)的,也可以是非平穩(wěn)的[11]。分解得到的IMF分量滿足:①零點數目與極值點數目相同或至多相差1;②函數由局部極大值點構成的包絡線和由局部極小值構成的包絡線的均值為零。分解過程通過一個稱為“篩選”的步驟來完成。“篩選”過程可以表示為[7]:

    (1)分別由原始信號x(t)的極小值點與極大值點用三次樣條插值得到x(t)的上下包絡線,計算上下包絡的均值m1,進而計算x(t)和m1的差值h1=x(t)-m1,判斷h1是否滿足IMF的兩個條件。若滿足,則h1為x(t)的第一個分量imf1。

    (2)若不滿足,則將h1作為新的信號繼續(xù)步驟(1),得到h11,判斷h11是否滿足條件。若不滿足則繼續(xù)(1)的步驟,直到重復k次后h1k滿足IMF的條件,則有imf1=h1k,求出原始信號與imf1的差值r1=x(t)-imf1。

    (3)將r1作為新的“原始”信號重復上述步驟,直到提取出第2個,第3個,直至第n個IMF分量imf2,imf3,…,imfn,有r2=r1-imf2,…,rn=rn-1-imfn,當IMF分量imfn或余項rn小于預先設定的值,或者余項rn已經成為單調函數時,整個篩選分解過程結束,本文以余項為單調函數作為分解的終止條件。

    經過上述步驟后,原始信號x(t)可分解為n個IMF分量和1個余項的和

    從上述的經驗模態(tài)分解方法可以看出,越早分解出來的IMF頻率越高[7],第一個分解出來的代表原信號的最高頻率成份,各個IMF的頻率幾乎是按2的負冪次方的形式遞減[6,12]。對于混有隨機噪聲的信號,其先分解出的IMF分量通常對應于信號的高頻噪聲[8]。若IMF組去除了先分解的幾個IMF,把其余的IMF組合起來形成一個信號,可以削弱信號的噪聲。可以看出,用EMD對陀螺信號消噪的關鍵在于從信號分解出的IMF組中辨識出噪聲IMF分量。

    2.2 處理異常噪聲

    由于外界環(huán)境的異常變化,陀螺信號會在局部時段內受到振幅較大的異常噪聲的干擾,導致信號輸出失真,對陀螺信號的處理首先要消除這類噪聲的影響,而小波對這種噪聲的消噪效果不佳。經驗模分解將信號按照其固有振動模式分成IMF組[13],而異常噪聲的頻率和振動模式一般既不同于隨機高頻噪聲,也不同于真實信號,對信號進行分解之后,異常噪聲一般被分解到其中一個或幾個IMF之中,只要能將相應的IMF進行識別并加以處理,就能減弱異常噪聲的影響。

    利用誤差理論中制定極限誤差的2σ準則,可以對陀螺信號的異常噪聲IMF分量進行識別。設原始信號為x(t),求得其采樣標準差為σ,計算IMF各個分量imfi的最大振幅Ai,對各分量按如下判斷準則進行處理

    異常噪聲作用的時間相對比較短,σ通常反映陀螺信號的整體觀測精度,若某個imfi的最大振幅大于2σ,則可認為它是異常噪聲分量。IMF分量的頻率各不相同,且如前面所述按2的負冪次方遞減,因而異常噪聲通常會被分解到少數IMF分量且是中高頻分量,非常便于進行識別和剔除。

    判別和處理異常噪聲是陀螺EMD消噪方法的第一步,除非陀螺的觀測環(huán)境非常好,否則都應該在對信號進行EMD分解后進行異常噪聲的辨識和剔除。

    圖1給出某次試驗中陀螺(標稱陀螺漂移為1°/h)靜基座下X軸8 000個歷元的輸出信號,在1 500~2 000、4 100~4 300、5 600~6 900歷元處受到異常噪聲的干擾。對該信號進行EMD分解,得到10個IMF分量,圖2顯示為前6個IMF以及最終的余項r10。可直觀看出,1 500~2 000、4 100~4 300歷元的異常噪聲主要分解在imf3內,5 600~6 900歷元的分量分解到imf2和imf3內。

    計算陀螺信號標準差為σ=34.02,則限差2σ=68.04,各IMF分量的最大振幅分別為A1=10.66,A2=216.38,A3=87.51,A4=39.85,A5=42.96,A6=12.28。按照式(2),imf2和 imf3的振幅超限,應該將這兩個分量置零,這與圖形所示的結果一致。

    圖1 陀螺原始信號Fig.1 Original gyro signal

    圖2 陀螺信號EMD分解的IMF分量Fig.2 IMFs from EMD of gyro signal

    2.3 確定噪聲IMF分量個數

    剔除異常噪聲后,需要消除陀螺信號中的高頻噪聲。設前m級IMF分量(已經過異常噪聲的處理)為高頻噪聲分量,則消噪過程可表示為

    若選取的m偏小,濾波后的信號x′(t)仍會受到噪聲干擾,若選取的m偏大,可能會把有用的信號濾去。為避免此問題,有文獻借鑒小波閾值消噪的思想提出對每個IMF分量也進行閾值處理[14]。然而,每個IMF分量都是具有特定物理意義的信號成分,反映原始信號某一方面的內在特征[7,13],利用閾值消噪的方法可能會破壞信號特征的完整性,且閾值消噪法存在確定閾值估計方法的難點,其實用性有待進一步研究。EMD消噪的關鍵還是在于確定噪聲分量的個數m。

    長度為N的信號序列x(t)和y(t)的采樣相關系數可定義為

    相關系數反映兩個信號之間的相似程度或依賴程度。設前m個IMF分量之和為若這m個分量都是噪聲IMF分量,不含原信號的有用成分,則按照式(3)濾波后的信號x′m(t)與原信號x(t)之間應該具有較強的依賴程度。因此可考慮由相關系數確定噪聲分量的個數,計算

    從m=1開始計算,當ρxx′m>c(c為常量,可取經驗值0.75~0.8)時,m=m+1,計算下一個ρxx′m,直到某個ρxx′m<c,停止計算,說明第m個分量imfm為信號的有用成分,不是噪聲IMF,取M=m-1作為噪聲分量個數。最終的EMD消噪結果為

    陀螺信號的EMD消噪方法歸納為:

    (1)對陀螺信號進行EMD分解,得到IMF分量;

    (2)按2σ準則處理異常噪聲IMF;

    (3)根據相關系數確定噪聲IMF個數,原信號減去噪聲IMF完成消噪過程。

    以上的EMD過程無需事先的參數設置,只與陀螺信號本身有關。陀螺信號的EMD分解是按照信號的固有振動模式進行的,采樣標準差σ和相關系數ρxx′m都由陀螺信號及其分解的IMF分量來計算,整個消噪過程不受外在因素影響。需要注意的是,常量c的取值及2σ準則的確定是EMD消噪中的重要環(huán)節(jié),如果選取不合適可能會影響最終的消噪結果,因此實用中需要結合實際的數據類型進行分析確定。

    確定圖2中噪聲分量個數,從m=4開始,相關系數均小于0.7,則取噪聲分量個數為3,EMD消噪的結果如圖3所示。

    圖3 EMD消噪信號Fig.3 EMD de-noised signal

    3 消噪方法對比

    3.1 消噪信號直接比較

    為對比小波消噪法與本文EMD消噪法,對圖1中的陀螺原始信號進行小波閾值消噪,采用軟閾值函數,閾值估計方法為SUREShrink閾值[3,15]。分別取不同小波基與分解尺度進行消噪試驗,消噪的效果總體一致,圖4、圖5分別為db8小波9尺度消噪結果和Haar小波7尺度消噪結果。

    圖4 db8小波消噪信號Fig.4 db8wavelet de-noised signal

    圖5 Haar小波消噪信號Fig.5 Haar wavelet de-noised signal

    比較圖1、圖3、圖4和圖5,可以看出:

    (1)EMD利用異常噪聲與有用信號和高頻噪聲不同的內在振動模態(tài)將異常噪聲分離,因此EMD消噪法既能消去隨機高頻噪聲,又能剔除異常噪聲;

    (2)小波消噪方法能有效地去除高頻噪聲,但是不能抵制異常噪聲的干擾;

    (3)小波分析的本質決定了它不能分離出異常噪聲,通過選取不同的小波基和分解尺度也無法克服這一缺陷。

    3.2 交疊式Allan方差對比分析

    采用交疊式Allan方差進一步對比兩種方法的消噪效果。Allan方差適合于分析非平穩(wěn)的隨機信號,是IEEE推薦的陀螺隨機誤差分析方法[16],其中的交疊式Allan方差在相同的置信水平下比普通Allan方差分析方法具有更大的置信區(qū)間,是對Allan方差的改進[17],其具體定義見文獻[17—18]。

    信號的Allan方差可表示為σ2(τ)(τ=nτ0,τ0為采樣間隔),則Allan標準差為σ(τ)~τ雙對數曲線圖可以描述陀螺的各種隨機誤差成分,不同的成分具有不同的斜率特性,這些成分包括量化噪聲,角度隨機游走,零偏不穩(wěn)定性,速率隨機游走和速率斜坡,可由Allan方差擬合得到,具體計算方法和它們在σ(τ)~τ圖中的斜率特性可參考文獻[18]。

    分別求取陀螺原始信號、小波消噪信號、EMD消噪信號的交疊式Allan標準差σ(τ),σ(τ)~τ雙對數曲線如圖6所示,其中db8和Haar小波的結果接近,只顯示db8小波結果。根據計算的Allan方差擬合得到的各種誤差成分的值如表1所示。

    圖6 σ(τ)~τ曲線對比Fig.6 Comparison ofσ(τ)~τcurve

    表1 陀螺誤差成分對比Tab.1 Comparison of error components of gyro

    分析圖6和表1的結果可知:

    (1)小波對陀螺信號具有一定的消噪作用,但是由于受到異常噪聲的干擾,小波消噪后的陀螺信號中仍然含有較大的誤差成分;

    (2)EMD消噪方法有效地削弱了陀螺各種誤差成分;

    (3)在τ=0.2~0.3s處,EMD及小波消噪后信號的Allan標準差與原始信號相比都沒有明顯改善,是因為受到殘余低頻有色噪聲的影響。一般可通過建立時間序列模型的方法削弱陀螺信號的有色噪聲成分。

    4 慣導解算與分析

    所用數據為一組靜態(tài)慣性測量單元(inertial measurement unit,IMU)數據,標稱陀螺漂移和加速度計偏置分別為1°/h和10-4g,IMU采樣頻率100Hz,取歷元282 930.83~283 232.82s(GPST,對應的周數為1 467)的數據進行試驗。利用282 930.83~283 112.83s的數據進行Kalman濾波精對準得到初始姿態(tài)角,初始航向角、俯仰角和翻滾角分別為8.916°、-0.013°、1.094°,分別采取3種方案對283 112.83~283 232.82s的IMU數據進行慣導解算,得到各歷元的速度值。各方案的加速度計數據均為原始輸出,但陀螺數據的處理方法不同。由于靜止狀態(tài)下慣導的速度真值實際上是“零”,因此,計算的速度值其實就是速度誤差。采用如下3種方案進行解算:

    方案1 使用IMU輸出的各軸陀螺原始數據解算;

    方案2 使用db8小波消噪后的各軸陀螺數據進行解算;

    方案3 使用本文EMD方法消噪后的各軸陀螺數據進行解算。

    3種方案的速度誤差見圖7~圖9,東向和北向的誤差結果類似,這里給出的是東向的結果,RMS的比較如表2所示。

    圖7 方案1東向速度誤差Fig.7 Veerror of scheme 1

    圖8 方案2東向速度誤差Fig.8 Veerror of scheme 2

    圖9 方案3東向速度誤差Fig.9 Veerror of scheme 3

    表2 3種方案RMS比較Tab.2 Comparison of RMS for three schemes(m/s)

    分析解算結果,可知:

    (1)受到加速度計誤差和殘余陀螺漂移的影響,單獨慣導解算的誤差隨著歷元數增加而不斷積累,3種方案的誤差均呈現(xiàn)擴大趨勢,需要引入其他的導航信息(如衛(wèi)星導航)才能予以修正。

    (2)方案2計算的速度誤差明顯小于方案1的結果,說明小波具有削弱陀螺漂移、減小慣導誤差的能力。

    (3)用EMD方法不僅削弱了陀螺漂移中隨機高頻噪聲的影響,而且消除了異常噪聲的干擾,相比小波又進一步減小了慣導誤差。

    5 結束語

    小波消噪方法能在一定程度上抑制陀螺漂移,但是它有一系列固有缺陷:不能消除異常噪聲;分解精度受測不準原理影響;對小波基、分解尺度和閾值估計方法等依賴性太大,需要繁瑣的調試才能達到好的效果。EMD消噪方法根據陀螺信號本身特性進行分解得到IMF,無需任何事先的參數設置,消噪過程中所需的閾值及相關系數由原始陀螺信號及其IMF分量計算得到,因此EMD方法的消噪過程只與陀螺信號本身有關。利用EMD方法對陀螺信號消噪,既能夠抵制異常噪聲干擾,又能削弱高頻噪聲,減少陀螺各項誤差成分,提高慣導解算的精度。

    [1] HUO Ju,WANG Shijing,YANG Ming,et al.Noise Processing of FOG Signal Based on Wavelet Threshold-value[J].Journal of Chinese Inertial Technology,2008,16(3):343-347.(霍炬,王石靜,楊明,等.基于小波變換閾值法處理光纖陀螺信號噪聲[J].中國慣性技術學報,2008,16(3):343-347.)

    [2] WAN Yanhui,QING Yongyuan.Application of Wavelet Analysis in Gyro Signal Filtering[J].Piezoelectrics and Acoustooptics,2005,27(4):455-457.(萬彥輝,秦永元.小波分析在陀螺信號濾波中的研究[J].壓電與聲光,2005,27(4):455-457.)

    [3] WU Fumei,YANG Yuanxi.GPS/INS Integrated Navigation by Adaptive Filtering Based on Wavelet Threshold De-noising[J].Acta Geodaetica et Cartographica Sinica,2007,36(2):124-128.(吳富梅,楊元喜.基于小波閾值消噪自適應濾波的GPS/INS組合導航[J].測繪學報,2007,36(2):124-128.)

    [4] TANG Wei,LI Shixin,LIU Luyuan,et al.Select of Wavelet Basis in Gyro Signal Processing[J].Journal of Chinese Inertial Technology,2002,10(5):28-30.(湯巍,李士心,劉魯源,等.關于陀螺信號處理中小波基選取的研究[J].中國慣性技術學報,2002,10(5):28-30.)

    [5] LIU Bin,JIANG Jinshui,YU Weikai.EMD De-noising Method and Its Application in the Processing for Rolling Mill Signals[J].Acta Metrologica Sinica,2009,30(1):73-77.(劉彬,蔣金水,于偉凱.EMD相關度消噪及其在軋機信號處理中的應用[J].計量學報,2009,30(1):73-77.)

    [6] JIANG Li,LI Changyun.A Study of Wavelet Threshold Filtering Based on Empirical Mode Decomposition[J].Signal Processing,2005,21(6):659-662.(江力,李長云.基于經驗模分解的小波閾值濾波方法研究[J].信號處理,2005,21(6):659-662.)

    [7] HUANG N E,SHEN Z,LONG S R,et al.The Empirical Mode Decomposition and the Hilbert Spectrum for Nonlinear and Nonstationary Time Series Analysis[C]∥Proceedings of Royal Society.London:[s.n.],1998:903-993.

    [8] DAI Wujiao,DING Xiaoli,ZHU Jianjun,et al.EMD Filter Method and Its Application in GPS Multipath[J].Acta Geodaetica et Cartographica Sinica,2006,35(11):321-327.(戴吾蛟,丁曉利,朱建軍,等.基于經驗模式分解的濾波消噪法及其在GPS多路徑效應中的應用[J].測繪學報,2006,35(11):321-327.)

    [9] KATHLEEN T.Stratospheric and Tropospheric Signals Extracted Using the Empirical Mode Decomposition Method[D].Washington:University of Washington,2003.

    [10] WANG Jian,GAO Jingxiang,WANG Jinling.GPS Baseline Solution Based on Empirical Mode Decomposition[J].Acta Geodaetica et Cartographica Sinica,2008,37(1):10-14.(王堅,高井祥,王金嶺.基于經驗模態(tài)分解的GPS基線解算模型[J].測繪學報,2008,37(1):10-14.)

    [11] CHEN Jun,XU Youli.Structural Damage Identification Based on Response of Forced Vibration[J].Journal of Vibration,Measurement and Diagnosis,2005,25(2):101-104.(陳雋,徐幼麟.經驗模分解在信號趨勢項提取中的應用[J].振動、測試與診斷,2005,25(2):101-104.)

    [12] FLANDRIN P,RILLING G,GONCALVES P.Empirical Mode Decomposition as a Filter Bank[J].IEEE Signal Processing Letters,2004,11(2):112-114.

    [13] HUANG N E,WU Zhaohua.A Review on Hilbert-Huang Transform:Method and Its Applications to Geophysical Studies[J].Reviews of Geophysics,2008,46:1-23.

    [14] BOUDRAA A O,CEXUS J C,SAIDI Z.EMD-based Signal Noise Reduction[J].International Journal of Signal Processing,2004,1:33-36.

    [15] DONOHO D L,JOHNSTONE I M.Adapting to Unknown Smoothness via Wavelet Shrinkage[J].Journal of the American Statistical Association,1995,90(12):1200-1224.

    [16] IEEE STD 952-1997.IEEE Standard Specification Format Guide and Test Procedure for Single-Axis Interferometric Fiber Optic Gyros[S].[S.l.]:IEEE Standard Board,1997.

    [17] IEEE 1139.Definitions of Physical Quantities for Fundamental Frequency and Time Metrology-Random Instabilities[S].[S.l.]:IEEE Standard Board,2008.

    [18] LI Xiaoying,HU Min,ZHANG Peng,et al.Applying Overlapping Allan Variance Theory to Better Stochastic Modeling of Microgyro[J].Journal of Northwestern Polytechnical University,2007,25(2):225-229.(李曉瑩,胡敏,張鵬,等.交疊式Allan方差在微機械陀螺隨機誤差辨識中的應用[J].西北工業(yè)大學學報,2007,25(2):225-229.)

    De-noising Method for Gyro Signal Based on EMD

    GAN Yu,SUI Lifen
    Institute of Surveying and Mapping,Information Engineering University,Zhengzhou 450052,China

    Gyro random drift is a remarkable factor that can affect the precision of inertial navigation system(INS).Wavelet de-noising method is poor in coping with exceptional noise,and it depends greatly on the selection of wavelet base and decomposition scale.Empirical mode decomposition(EMD)de-noising method for gyro signal is presented.The signal is decomposed into an intrinsic mode function(IMF)group.Based on this group,IMFs of exceptional noise are first disposed by2σcriterion and then the number of IMFs of high frequency noise is determined by correlation coefficient.The de-noising process is finally done by removing the noisy IMFs.Detailed comparison between EMD method and wavelet method is given.Overlapping Allan variance is used to analyze the effect of the two methods,and the applicable ability of EMD method is tested through an INS calculation.It is shown that EMD method outperforms wavelet method in removing exceptional noise and is more efficient in weakening random drift.

    gyro random drift;wavelet;empirical mode decomposition;de-noising

    GAN Yu(1988—),male,postgraduate,majors in dynamic geodetic data processing.

    1001-1595(2011)06-0745-06

    P228

    A

    國家自然科學基金(40974010);信息工程大學測繪學院碩士學位論文創(chuàng)新與創(chuàng)優(yōu)基金(S201101)

    宋啟凡)

    2010-12-06

    2010-12-30

    甘雨(1988—),男,碩士生,主要從事動態(tài)大地測量數據處理研究。

    E-mail:ganyu099@163.com

    猜你喜歡
    陀螺小波方差
    方差怎么算
    構造Daubechies小波的一些注記
    科技風(2021年19期)2021-09-07 14:04:29
    概率與統(tǒng)計(2)——離散型隨機變量的期望與方差
    基于MATLAB的小波降噪研究
    電子制作(2019年13期)2020-01-14 03:15:32
    做個紙陀螺
    計算方差用哪個公式
    玩陀螺
    學生天地(2019年6期)2019-03-07 01:10:46
    陀螺轉轉轉
    軍事文摘(2018年24期)2018-12-26 00:58:18
    我最喜歡的陀螺
    快樂語文(2018年36期)2018-03-12 00:56:02
    基于改進的G-SVS LMS 與冗余提升小波的滾動軸承故障診斷
    免费av中文字幕在线| 国产免费视频播放在线视频| 97精品久久久久久久久久精品| 老熟女久久久| 老汉色av国产亚洲站长工具| 伊人亚洲综合成人网| 最新在线观看一区二区三区 | 一二三四在线观看免费中文在| 成年av动漫网址| 欧美老熟妇乱子伦牲交| 婷婷色麻豆天堂久久| 久久久国产一区二区| 国语对白做爰xxxⅹ性视频网站| 建设人人有责人人尽责人人享有的| 国产精品久久久av美女十八| 大型av网站在线播放| 成人国产av品久久久| 精品一区二区三区四区五区乱码 | 赤兔流量卡办理| 久久精品人人爽人人爽视色| 日本午夜av视频| 丰满迷人的少妇在线观看| 中文乱码字字幕精品一区二区三区| 亚洲少妇的诱惑av| 精品国产乱码久久久久久男人| 日韩制服骚丝袜av| 国语对白做爰xxxⅹ性视频网站| 乱人伦中国视频| 亚洲av日韩精品久久久久久密 | 久久久久国产精品人妻一区二区| 蜜桃在线观看..| 电影成人av| 十分钟在线观看高清视频www| 丰满少妇做爰视频| 国产又爽黄色视频| 日韩av不卡免费在线播放| 日韩av免费高清视频| 免费日韩欧美在线观看| 国产淫语在线视频| 亚洲激情五月婷婷啪啪| 亚洲av电影在线观看一区二区三区| 777久久人妻少妇嫩草av网站| 美女午夜性视频免费| 免费女性裸体啪啪无遮挡网站| 最新在线观看一区二区三区 | 黄色一级大片看看| 日本wwww免费看| av国产精品久久久久影院| 丰满饥渴人妻一区二区三| 日韩电影二区| 久久久久久久久久久久大奶| 亚洲国产精品一区三区| 久久久久国产精品人妻一区二区| av视频免费观看在线观看| 伊人亚洲综合成人网| 国产亚洲午夜精品一区二区久久| 韩国高清视频一区二区三区| 1024香蕉在线观看| 日本av免费视频播放| 日韩一区二区三区影片| 一级黄片播放器| a级片在线免费高清观看视频| 汤姆久久久久久久影院中文字幕| av国产精品久久久久影院| 国产欧美日韩一区二区三区在线| 中文字幕另类日韩欧美亚洲嫩草| av网站在线播放免费| 曰老女人黄片| 午夜福利,免费看| 久久精品亚洲av国产电影网| 国产免费视频播放在线视频| 亚洲欧美色中文字幕在线| 性色av乱码一区二区三区2| 国产精品人妻久久久影院| 亚洲自偷自拍图片 自拍| 欧美精品高潮呻吟av久久| 久久人妻熟女aⅴ| 777米奇影视久久| 亚洲一区中文字幕在线| 久久亚洲国产成人精品v| 日本av免费视频播放| av在线播放精品| 纯流量卡能插随身wifi吗| 岛国毛片在线播放| 亚洲欧美成人综合另类久久久| 丝袜美足系列| 91字幕亚洲| 欧美在线黄色| 久久精品亚洲av国产电影网| 免费久久久久久久精品成人欧美视频| 亚洲国产精品成人久久小说| 少妇被粗大的猛进出69影院| 欧美人与性动交α欧美精品济南到| 久久青草综合色| 国产爽快片一区二区三区| 亚洲人成77777在线视频| 久久久久久亚洲精品国产蜜桃av| 欧美日韩国产mv在线观看视频| 涩涩av久久男人的天堂| 免费看十八禁软件| 久久精品久久精品一区二区三区| 久久久国产欧美日韩av| 麻豆av在线久日| 亚洲国产中文字幕在线视频| 亚洲精品国产区一区二| 国产男女超爽视频在线观看| 亚洲国产中文字幕在线视频| 超碰97精品在线观看| 九草在线视频观看| 国产成人91sexporn| 国产精品熟女久久久久浪| 久久狼人影院| 又大又爽又粗| 夫妻性生交免费视频一级片| 桃花免费在线播放| 午夜两性在线视频| 国产一区二区三区综合在线观看| 久久久久久亚洲精品国产蜜桃av| 麻豆乱淫一区二区| 免费高清在线观看日韩| 欧美日韩亚洲国产一区二区在线观看 | 亚洲欧美日韩另类电影网站| 黄片小视频在线播放| 一区在线观看完整版| 亚洲,一卡二卡三卡| av线在线观看网站| 日韩av在线免费看完整版不卡| 青春草亚洲视频在线观看| 亚洲成人免费电影在线观看 | 日本欧美国产在线视频| 亚洲男人天堂网一区| 最黄视频免费看| 女人爽到高潮嗷嗷叫在线视频| 国产真人三级小视频在线观看| 1024视频免费在线观看| 亚洲中文日韩欧美视频| 国产在线观看jvid| 后天国语完整版免费观看| 真人做人爱边吃奶动态| 一区福利在线观看| 久久热在线av| 狠狠精品人妻久久久久久综合| 大陆偷拍与自拍| 婷婷丁香在线五月| 飞空精品影院首页| 亚洲九九香蕉| 国产精品一区二区免费欧美 | 国产男女超爽视频在线观看| av天堂久久9| 欧美 日韩 精品 国产| 人体艺术视频欧美日本| 美女扒开内裤让男人捅视频| 热re99久久精品国产66热6| 国产精品一区二区在线观看99| 爱豆传媒免费全集在线观看| 激情视频va一区二区三区| 久久精品人人爽人人爽视色| 久久久精品94久久精品| 成人国产一区最新在线观看 | 天堂俺去俺来也www色官网| 啦啦啦 在线观看视频| 国产精品亚洲av一区麻豆| 亚洲精品一卡2卡三卡4卡5卡 | 男女之事视频高清在线观看 | 欧美在线黄色| 久久精品久久久久久噜噜老黄| 最近手机中文字幕大全| 深夜精品福利| 18禁裸乳无遮挡动漫免费视频| 亚洲欧美日韩另类电影网站| 免费高清在线观看日韩| www.熟女人妻精品国产| 伊人久久大香线蕉亚洲五| 最近中文字幕2019免费版| 又粗又硬又长又爽又黄的视频| 国产视频一区二区在线看| 大香蕉久久网| 一级毛片电影观看| 精品亚洲成国产av| 亚洲美女黄色视频免费看| 91精品伊人久久大香线蕉| 国产精品香港三级国产av潘金莲 | 新久久久久国产一级毛片| 欧美精品一区二区大全| 99九九在线精品视频| av视频免费观看在线观看| 国产欧美日韩综合在线一区二区| 欧美精品亚洲一区二区| 精品久久久久久久毛片微露脸 | 国产免费视频播放在线视频| 日韩 欧美 亚洲 中文字幕| 日本av免费视频播放| 亚洲欧美清纯卡通| 日本黄色日本黄色录像| 日本猛色少妇xxxxx猛交久久| 少妇粗大呻吟视频| 男女之事视频高清在线观看 | 人妻 亚洲 视频| 免费观看a级毛片全部| 国语对白做爰xxxⅹ性视频网站| 国产又色又爽无遮挡免| 少妇的丰满在线观看| 麻豆乱淫一区二区| 国产成人精品无人区| 久久人人爽av亚洲精品天堂| 午夜免费观看性视频| 国产极品粉嫩免费观看在线| a 毛片基地| 国产在线免费精品| 黑人欧美特级aaaaaa片| 在线观看一区二区三区激情| 少妇人妻 视频| 狂野欧美激情性bbbbbb| 天天添夜夜摸| 精品一区二区三区av网在线观看 | 热99国产精品久久久久久7| 亚洲久久久国产精品| 成年女人毛片免费观看观看9 | 色精品久久人妻99蜜桃| 一本综合久久免费| 看免费成人av毛片| www.熟女人妻精品国产| 国产免费又黄又爽又色| 99热网站在线观看| 久久久久久久大尺度免费视频| 亚洲欧美一区二区三区黑人| 王馨瑶露胸无遮挡在线观看| 亚洲成av片中文字幕在线观看| 狂野欧美激情性xxxx| 久久久久久免费高清国产稀缺| 国产亚洲精品久久久久5区| 我的亚洲天堂| 国产日韩欧美在线精品| 日本vs欧美在线观看视频| 亚洲成人免费电影在线观看 | 性色av一级| 国产亚洲精品久久久久5区| 人人妻人人添人人爽欧美一区卜| 亚洲人成网站在线观看播放| av又黄又爽大尺度在线免费看| 婷婷成人精品国产| 亚洲精品自拍成人| 国产老妇伦熟女老妇高清| 欧美激情极品国产一区二区三区| 黑人欧美特级aaaaaa片| 日韩精品免费视频一区二区三区| 欧美日韩福利视频一区二区| 丰满饥渴人妻一区二区三| 亚洲国产精品一区二区三区在线| 色播在线永久视频| 啦啦啦 在线观看视频| 欧美少妇被猛烈插入视频| 国产野战对白在线观看| 18禁黄网站禁片午夜丰满| 午夜激情av网站| 精品少妇久久久久久888优播| 在线观看免费视频网站a站| 国产成人av教育| 一级a爱视频在线免费观看| 亚洲人成77777在线视频| 国产一区二区 视频在线| 一个人免费看片子| 一区二区av电影网| 在线观看www视频免费| 国产女主播在线喷水免费视频网站| 免费一级毛片在线播放高清视频 | 国产在线一区二区三区精| 嫁个100分男人电影在线观看 | 欧美亚洲 丝袜 人妻 在线| 一本久久精品| 在线精品无人区一区二区三| 久久午夜综合久久蜜桃| 美女中出高潮动态图| 国产一区二区在线观看av| 亚洲精品中文字幕在线视频| 日韩人妻精品一区2区三区| 无限看片的www在线观看| 亚洲国产日韩一区二区| 久久ye,这里只有精品| 国产野战对白在线观看| 精品久久久久久电影网| 在线亚洲精品国产二区图片欧美| 一级,二级,三级黄色视频| 欧美乱码精品一区二区三区| 免费高清在线观看日韩| 大片免费播放器 马上看| 精品少妇黑人巨大在线播放| 中国国产av一级| 亚洲第一av免费看| 亚洲国产av影院在线观看| 精品国产超薄肉色丝袜足j| 欧美精品人与动牲交sv欧美| 国产精品久久久久久精品电影小说| 天天添夜夜摸| av国产久精品久网站免费入址| 国产一区二区三区综合在线观看| 色视频在线一区二区三区| 成人午夜精彩视频在线观看| 久久久国产一区二区| 午夜福利视频在线观看免费| 日韩制服骚丝袜av| 中文字幕亚洲精品专区| 如日韩欧美国产精品一区二区三区| a 毛片基地| 黄频高清免费视频| 一本大道久久a久久精品| 日韩熟女老妇一区二区性免费视频| 天堂8中文在线网| 欧美变态另类bdsm刘玥| 精品一区在线观看国产| 午夜免费观看性视频| 大码成人一级视频| 国产精品久久久av美女十八| 各种免费的搞黄视频| 成年av动漫网址| 一区在线观看完整版| 激情视频va一区二区三区| 色精品久久人妻99蜜桃| 少妇 在线观看| 亚洲黑人精品在线| 黄色视频在线播放观看不卡| 美女大奶头黄色视频| 无限看片的www在线观看| 免费高清在线观看日韩| 午夜91福利影院| 成人免费观看视频高清| 中文字幕高清在线视频| 国产一区二区三区综合在线观看| 亚洲色图综合在线观看| 99国产精品免费福利视频| 黑人猛操日本美女一级片| 日韩视频在线欧美| 2018国产大陆天天弄谢| 精品视频人人做人人爽| 国产亚洲午夜精品一区二区久久| 久久精品aⅴ一区二区三区四区| 国产精品人妻久久久影院| 日本91视频免费播放| 一区二区av电影网| 精品人妻一区二区三区麻豆| 日本vs欧美在线观看视频| 亚洲av成人不卡在线观看播放网 | 国产成人a∨麻豆精品| 在线天堂中文资源库| 91老司机精品| 国产在线视频一区二区| xxxhd国产人妻xxx| 黄色一级大片看看| 丝袜美足系列| 叶爱在线成人免费视频播放| 国产在视频线精品| 少妇裸体淫交视频免费看高清 | 午夜激情av网站| 人成视频在线观看免费观看| 午夜福利影视在线免费观看| 51午夜福利影视在线观看| 夜夜骑夜夜射夜夜干| videosex国产| 美女福利国产在线| 最近中文字幕2019免费版| 久久精品国产综合久久久| 波野结衣二区三区在线| 久久 成人 亚洲| 岛国毛片在线播放| 99久久99久久久精品蜜桃| 国产欧美日韩精品亚洲av| 久久久久久久久免费视频了| 国产在线视频一区二区| 精品高清国产在线一区| 久久人妻福利社区极品人妻图片 | 最新的欧美精品一区二区| 久热爱精品视频在线9| 亚洲激情五月婷婷啪啪| 亚洲五月色婷婷综合| 婷婷色麻豆天堂久久| 国产亚洲av高清不卡| 久久久亚洲精品成人影院| 亚洲欧美一区二区三区国产| 这个男人来自地球电影免费观看| 在线观看www视频免费| 欧美精品一区二区免费开放| 国产欧美日韩一区二区三区在线| 精品久久久久久电影网| 国产精品麻豆人妻色哟哟久久| 首页视频小说图片口味搜索 | 各种免费的搞黄视频| 亚洲国产毛片av蜜桃av| 国产97色在线日韩免费| 欧美在线黄色| 国产高清videossex| 宅男免费午夜| 亚洲色图综合在线观看| 在线精品无人区一区二区三| 国精品久久久久久国模美| 岛国毛片在线播放| av电影中文网址| 国产免费现黄频在线看| 亚洲 欧美一区二区三区| 国产在线观看jvid| 极品人妻少妇av视频| 亚洲 国产 在线| 又大又黄又爽视频免费| 国产精品国产三级国产专区5o| 成年美女黄网站色视频大全免费| 国产精品一区二区精品视频观看| 国产色视频综合| 亚洲av综合色区一区| 观看av在线不卡| 亚洲 国产 在线| 操美女的视频在线观看| 美女午夜性视频免费| 国产亚洲午夜精品一区二区久久| 人妻一区二区av| 啦啦啦中文免费视频观看日本| 女人久久www免费人成看片| 色播在线永久视频| 大话2 男鬼变身卡| 精品福利观看| 久久性视频一级片| 18禁观看日本| 桃花免费在线播放| 一区二区日韩欧美中文字幕| 欧美日韩视频高清一区二区三区二| 在线天堂中文资源库| 亚洲国产欧美网| 久久中文字幕一级| 国产精品99久久99久久久不卡| 久久久亚洲精品成人影院| 国产亚洲一区二区精品| 午夜老司机福利片| 日韩视频在线欧美| 久久影院123| 在线精品无人区一区二区三| 亚洲av片天天在线观看| 国产成人精品无人区| 国产视频一区二区在线看| 亚洲天堂av无毛| 亚洲图色成人| 成人三级做爰电影| 天天影视国产精品| 一级毛片电影观看| 啦啦啦 在线观看视频| 一级黄片播放器| 欧美精品一区二区免费开放| 亚洲av成人不卡在线观看播放网 | 久久久久网色| 国产成人影院久久av| 丰满饥渴人妻一区二区三| 日韩伦理黄色片| 午夜福利视频在线观看免费| 男女高潮啪啪啪动态图| 亚洲国产日韩一区二区| 我要看黄色一级片免费的| av网站免费在线观看视频| 久久精品久久久久久久性| bbb黄色大片| 国产精品麻豆人妻色哟哟久久| 国产成人啪精品午夜网站| 亚洲av美国av| 99九九在线精品视频| 高清黄色对白视频在线免费看| 秋霞在线观看毛片| 亚洲av成人精品一二三区| 99久久人妻综合| 又粗又硬又长又爽又黄的视频| 宅男免费午夜| 午夜日韩欧美国产| 国产精品久久久久久人妻精品电影 | 国产成人啪精品午夜网站| 人妻一区二区av| 国产成人精品久久二区二区免费| 高清黄色对白视频在线免费看| 午夜两性在线视频| 亚洲国产av新网站| 亚洲国产看品久久| 欧美日韩视频高清一区二区三区二| 亚洲成色77777| 超碰成人久久| 日韩大码丰满熟妇| 人人妻,人人澡人人爽秒播 | 欧美av亚洲av综合av国产av| av天堂久久9| 五月开心婷婷网| 在现免费观看毛片| 嫁个100分男人电影在线观看 | 国产亚洲午夜精品一区二区久久| 国产一区二区三区综合在线观看| 又黄又粗又硬又大视频| 最新的欧美精品一区二区| 91字幕亚洲| av网站在线播放免费| 欧美人与善性xxx| 天天躁日日躁夜夜躁夜夜| 肉色欧美久久久久久久蜜桃| 亚洲欧美日韩另类电影网站| 国产麻豆69| 精品一品国产午夜福利视频| 高潮久久久久久久久久久不卡| 成在线人永久免费视频| 亚洲精品第二区| 午夜老司机福利片| 两个人免费观看高清视频| 日本av手机在线免费观看| 中文字幕另类日韩欧美亚洲嫩草| 国产精品国产三级国产专区5o| 日韩制服骚丝袜av| 久久毛片免费看一区二区三区| 亚洲,一卡二卡三卡| 亚洲欧美精品自产自拍| 妹子高潮喷水视频| 一本—道久久a久久精品蜜桃钙片| 少妇被粗大的猛进出69影院| 国产成人a∨麻豆精品| 亚洲黑人精品在线| 亚洲专区中文字幕在线| 亚洲欧洲日产国产| 视频区欧美日本亚洲| 下体分泌物呈黄色| 日韩伦理黄色片| 精品视频人人做人人爽| 啦啦啦视频在线资源免费观看| 国产成人精品久久二区二区91| 性色av一级| 日韩大片免费观看网站| 亚洲男人天堂网一区| 一本久久精品| 成人手机av| 在线观看www视频免费| 中文欧美无线码| 美女视频免费永久观看网站| 欧美日韩一级在线毛片| 亚洲 国产 在线| 日韩一卡2卡3卡4卡2021年| 欧美+亚洲+日韩+国产| 国产麻豆69| 国产xxxxx性猛交| 久久精品久久精品一区二区三区| 大型av网站在线播放| 在线亚洲精品国产二区图片欧美| 日本猛色少妇xxxxx猛交久久| 九草在线视频观看| 水蜜桃什么品种好| 日本一区二区免费在线视频| 日韩电影二区| 男男h啪啪无遮挡| 1024视频免费在线观看| 新久久久久国产一级毛片| 久久国产精品人妻蜜桃| 叶爱在线成人免费视频播放| 妹子高潮喷水视频| 自拍欧美九色日韩亚洲蝌蚪91| 日韩,欧美,国产一区二区三区| a级毛片黄视频| bbb黄色大片| 午夜免费男女啪啪视频观看| 精品久久久精品久久久| 欧美精品一区二区免费开放| 国产xxxxx性猛交| 啦啦啦在线免费观看视频4| 高清黄色对白视频在线免费看| 午夜视频精品福利| 欧美另类一区| 99精国产麻豆久久婷婷| 妹子高潮喷水视频| 2018国产大陆天天弄谢| 久久久国产精品麻豆| 国产亚洲av高清不卡| 成人国语在线视频| 久久鲁丝午夜福利片| 一边摸一边抽搐一进一出视频| 在线观看免费午夜福利视频| 午夜免费观看性视频| 9191精品国产免费久久| 日韩大码丰满熟妇| 午夜影院在线不卡| 免费一级毛片在线播放高清视频 | 大话2 男鬼变身卡| 狠狠精品人妻久久久久久综合| 男女国产视频网站| 少妇裸体淫交视频免费看高清 | 精品一品国产午夜福利视频| 国产成人精品久久二区二区免费| 嫁个100分男人电影在线观看 | 国产一区二区 视频在线| 日本a在线网址| 午夜精品国产一区二区电影| 亚洲欧洲精品一区二区精品久久久| 一级毛片黄色毛片免费观看视频| 每晚都被弄得嗷嗷叫到高潮| √禁漫天堂资源中文www| 天天添夜夜摸| 老汉色∧v一级毛片| 大码成人一级视频| 国产免费现黄频在线看| 激情视频va一区二区三区| av网站在线播放免费| 精品久久久久久久毛片微露脸 | 免费黄频网站在线观看国产| 亚洲色图综合在线观看| 亚洲av国产av综合av卡| 国产在线观看jvid| 男男h啪啪无遮挡| 久久性视频一级片| 一区福利在线观看| 精品久久久精品久久久| 亚洲一卡2卡3卡4卡5卡精品中文| 国产av一区二区精品久久| 夜夜骑夜夜射夜夜干| 国产欧美日韩综合在线一区二区| 国产一卡二卡三卡精品| 在线 av 中文字幕| 老司机午夜十八禁免费视频| 一二三四社区在线视频社区8| 9热在线视频观看99| 欧美黄色淫秽网站|