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

    一種改進(jìn)的基于EMD分解的超寬帶信號消噪算法?

    2012-07-01 18:04:23王海梁熊華鋼吳慶劉成
    電訊技術(shù) 2012年4期
    關(guān)鍵詞:超寬帶階數(shù)信噪比

    王海梁,熊華鋼,吳慶,劉成

    (北京航空航天大學(xué)電子信息工程學(xué)院,北京100191)

    一種改進(jìn)的基于EMD分解的超寬帶信號消噪算法?

    王海梁,熊華鋼,吳慶,劉成

    (北京航空航天大學(xué)電子信息工程學(xué)院,北京100191)

    針對低信噪比超寬帶信號的消噪問題,提出一種改進(jìn)的基于經(jīng)驗(yàn)?zāi)J椒纸猓‥MD)的消噪算法。該算法首先對含噪信號進(jìn)行EMD分解,得到多個(gè)固有模態(tài)函數(shù)(IMF)分量,然后選取高階IMF重構(gòu)原信號,達(dá)到消噪的目的。針對對UWB信號的IMF重構(gòu)過程中階數(shù)閾值難以確定的問題,通過數(shù)值仿真的方法,得到信號分量和噪聲分量在不同階IMF上的能量分布特性;在對所得特性進(jìn)行分析的基礎(chǔ)上,設(shè)計(jì)了一種數(shù)據(jù)自適應(yīng)的階數(shù)閾值選取算法,解決了EMD消噪中的階數(shù)閾值選取問題。仿真結(jié)果表明,EMD消噪算法能夠在較低信噪比下提供平均10 dB的信噪比增益,可以有效地對超寬帶信號進(jìn)行消噪。

    超寬帶;信號消噪;經(jīng)驗(yàn)?zāi)J椒纸?/p>

    1 引言

    脈沖超寬帶(Ultra-wideband,UWB)無線電技術(shù)具有高隱蔽性、抗多徑效應(yīng)、較高的距離分辨力等優(yōu)良特性,在通信、定位和雷達(dá)等領(lǐng)域都有應(yīng)用,是當(dāng)前無線電技術(shù)的一個(gè)研究熱點(diǎn)[1-2]。受發(fā)射功率的限制,脈沖UWB信號通常信噪比較低,因此,對UWB信號的消噪是UWB信號處理中的一個(gè)重要課題。

    處理信號消噪問題的傳統(tǒng)方法是利用Fourier變換或者小波變換消噪。由于UWB信號具有高帶寬、低功率譜的特性,利用Fourier變換進(jìn)行消噪的方法對UWB信號難以取得較好的消噪效果;而小波消噪需要預(yù)先選擇好基函數(shù),對于不同種類的UWB脈沖波形,很難找到一種通用的基函數(shù)選擇方法,目前,尚未有針對UWB信號的有效消噪方法的相關(guān)文獻(xiàn)見諸報(bào)道。

    經(jīng)驗(yàn)?zāi)B(tài)分解(Empirical Mode Decomposition,EMD)是近年來提出的一種時(shí)頻處理算法[3],該算法無需預(yù)先指定基函數(shù),特別適用于對非平穩(wěn)、非線性信號的處理。文獻(xiàn)[4-7]提出對EMD分解得到的固有模態(tài)函數(shù)(Intrinsic Mode Function,IMF)進(jìn)行篩選,選取部分IMF重構(gòu)原信號可以達(dá)到消噪的目的。

    本文將IMF重構(gòu)消噪的思想應(yīng)用于超寬帶信號消噪問題的解決中。針對重構(gòu)過程中IMF的階數(shù)閾值選取問題,設(shè)計(jì)了一種自適應(yīng)的閾值選取算法。該選取算法以超寬帶信號的各階IMF能量分布特性為基礎(chǔ),從數(shù)據(jù)中直接提取閾值信息,無需對信號和噪聲做特殊假設(shè)。IMF重構(gòu)方法和自適應(yīng)的閾值選取方法相結(jié)合,構(gòu)成了基于EMD分解的超寬帶信號消噪方法。通過對人工合成信號和實(shí)際采集信號的仿真實(shí)驗(yàn),驗(yàn)證了該消噪方法的有效性。

    2 超寬帶信號消噪問題

    2.1 脈沖超寬帶信號

    脈沖超寬帶信號由一系列時(shí)間上不連續(xù)的窄脈沖構(gòu)成,脈沖的寬度通常在納秒或者百皮秒的量級,根據(jù)使用目的的不同,脈沖的幅度或者位置可以被編碼,形成移位的窄脈沖序列。脈沖超寬帶信號一般具有較大的帶寬和較低的功率譜密度,受到發(fā)射功率的的限制,一般接收到的超寬帶信號具有較低的信噪比。圖1給出了文獻(xiàn)[1]中描述的跳時(shí)超寬帶信號原始脈沖串和附加了噪聲的含噪信號波形圖。

    2.2 超寬帶信號的消噪

    已知含噪信號y(t)=x(t)+n(t),原始信號為x(t),n(t)是高斯白噪聲。信號消噪即求對原始信號的估計(jì)信號^x(t),使得均方誤差

    盡可能小。

    為衡量消噪的效果,定義信號輸入信噪比Sin和消噪后的輸出信噪比Sout分別為

    那么可以用信噪比增益

    來描述消噪的性能。

    3 EMD分解

    EMD方法利用一個(gè)迭代的篩分過程將一個(gè)給定的信號s(t)分解成L個(gè)IMF分量hi(t)和一個(gè)殘余分量rL(t),每個(gè)IMF代表了信號s(t)在某個(gè)時(shí)間尺度或者頻帶上的細(xì)節(jié)。每個(gè)IMF要滿足兩個(gè)條件[3]:一是IMF中的極值點(diǎn)和過零點(diǎn)的數(shù)目最多相差1,二是上包絡(luò)線和下包絡(luò)線的局部均值為0。

    篩分過程按照從高頻到低頻的順序形成不同階的IMF,最先得到的IMF(低階IMF)具有最高的頻率,最后得到的IMF(高階)具有最低的頻率。EMD的主要優(yōu)點(diǎn)是篩分過程是數(shù)據(jù)驅(qū)動(dòng)的,不需要預(yù)先構(gòu)造基函數(shù)。

    由文獻(xiàn)[3],IMF的篩分過程是一個(gè)迭代的過程,每次迭代篩分出一個(gè)IMF分量,對于第j次迭代,rj-1(t)為當(dāng)前殘余分量,并且當(dāng)j=1時(shí),有r0(t)=s(t)。迭代過程如下:

    (1)令hj,i(t)←rj-1(t),i←1;(2)求hj,i(t)的極大值和極小值,并利用3次樣條插值求得極大值確定的上包絡(luò)線和極小值確定的下包絡(luò)線;

    (3)求上下包絡(luò)線的均值mj,i(t);

    (4)更新hj,i+1(t)←hj,i(t)-mj,i(t),i←i+1;

    (5)判斷hj,i(t)是否滿足IMF的兩個(gè)條件,若滿足,則hj(t)=hj,i(t)為第j階IMF分量,若不滿足,重復(fù)步驟(2)~(4);

    (6)更新殘余值rj=rj-1(t)-hj(t),j←j+1,如果rj(t)中的極值點(diǎn)少于兩個(gè),終止迭代,否則,從第1步開始繼續(xù)迭代過程。

    篩分過程結(jié)束后,原信號s(t)可以表示成分解得到的各階IMF hj(t)和殘余值rL(t)的組合:

    式中,hj(t)代表不同頻段的信號分量,rL(t)代表信號整體的變化趨勢。

    4 基于EMD的超寬帶信號消噪

    4.1 IMF重構(gòu)消噪

    EMD過程將信號分解成一系列代表不同特征頻率的IMF分量,低階IMF代表信號的高頻部分,高階IMF對應(yīng)信號的低頻部分。由文獻(xiàn)[5,8]可知,對于白噪聲,EMD分解后噪聲的能量將主要集中于高頻IMF分量中。當(dāng)對過采樣的含噪信號進(jìn)行EMD分解時(shí),由文獻(xiàn)[6],可以認(rèn)為噪聲的主要能量分布在高頻IMF中(前幾階IMF),信號的主要能量分布在低頻IMF中(后幾階IMF),并隨著階數(shù)的降低而遞減。因此,可以認(rèn)為存在某個(gè)階數(shù)閾值th,對于任何階數(shù)大于等于該閾值的IMF,其中信號部分占據(jù)主要成分;階數(shù)小于該閾值的IMF中,噪聲占據(jù)主要成分。

    通過濾除階數(shù)小于該閾值的IMF,利用高階IMF和趨勢分量重構(gòu)原信號,可以達(dá)到消噪的目的,重構(gòu)信號可以表示為

    能否選擇合理的階數(shù)閾值是有效實(shí)現(xiàn)該重構(gòu)算法的關(guān)鍵。在實(shí)際應(yīng)用中,由于帶寬和采樣率的變化,使得對于不同信號的分解,有效部分對應(yīng)的階數(shù)不固定,無法預(yù)先確定。下面通過對超寬帶含噪信號和噪聲的EMD分解結(jié)果的對比分析,給出一種閾值選取算法。

    4.2 超寬帶信號IMF重構(gòu)的閾值選取

    目前,EMD方法還沒有比較完備的解析解釋,大多數(shù)針對EMD的研究結(jié)果都是在數(shù)值仿真的基礎(chǔ)上通過對實(shí)驗(yàn)結(jié)果進(jìn)行總結(jié)得出的。本文也采取計(jì)算機(jī)數(shù)值仿真的方式研究超寬帶信號EMD分解后不同階IMF的能量分布特性,從中找出選取階數(shù)閾值的方法。

    仿真采用跳時(shí)UWB信號,采用高斯0階脈沖作為基本脈沖波形,脈沖寬度為1 ns,脈沖的重復(fù)周期為10 ns,信號總長度為100 ns,對該信號疊加高斯白噪聲,形成信噪比為-5 dB的含噪信號。仿真在20 GHz和50 GHz兩種采樣率下進(jìn)行,在每種采樣率下各運(yùn)行了10 000次仿真,每次仿真生成含噪信號和噪聲,并分別進(jìn)行EMD分解,求各階IMF的能量,最后對10 000次結(jié)果求平均值,得到了含噪信號和噪聲的各階IMF能量分布圖,如圖2和圖3所示。

    圖2 采樣率20 GHz時(shí)噪聲和含噪信號的各階IMF能量對比(對數(shù)坐標(biāo))Fig.2 IMF energy of noise and noised signal at 20 GHz sample frequency(semi-log axis)

    圖3 采樣率50 GHz時(shí)噪聲和含噪信號的各階IMF能量對比(對數(shù)坐標(biāo))Fig.3 IMF energy of noise and noised signal at 50 GHz sample frequency(semi-log axis)

    從圖中可得,對于噪聲的EMD分解,各階IMF的能量隨著階數(shù)的增加而逐漸減弱;在對數(shù)坐標(biāo)軸中,各階IMF的能量近似構(gòu)成一條下降的直線,說明各階IMF的能量近似構(gòu)成等比數(shù)列,這與文獻(xiàn)[8]中對白噪聲IMF的研究結(jié)果是一致的。

    對于含噪信號,其前幾階IMF和噪聲的IMF曲線基本重合,當(dāng)階數(shù)到達(dá)某個(gè)值時(shí)(20 GHz采樣率時(shí)為3,50 GHz采樣率時(shí)為4),兩條曲線的差值出現(xiàn)躍變,含噪信號該階IMF的能量遠(yuǎn)大于對應(yīng)的噪聲IMF的能量,說明在該階IMF中,信號能量開始占據(jù)主要成分,這個(gè)階數(shù)即需要求的階數(shù)閾值。

    令含噪信號第i階IMF的能量密度為Es,i,噪聲第i階IMF的能量密度En,i, 兩者的比值為

    那么,當(dāng)βi大于某個(gè)常數(shù)λ(一般可采用1.5)時(shí),可以認(rèn)為該IMF中信號能量占據(jù)主要成分。于是,可以選擇滿足βi>λ的最小i值為階數(shù)閾值。

    在實(shí)際應(yīng)用中,無法準(zhǔn)確獲得噪聲的IMF,也無法獲得En,i的準(zhǔn)確值。在過采樣的條件下,可以假設(shè)含噪信號的前兩階IMF和噪聲前兩階IMF近似相等,利用噪聲IMF能量的等比數(shù)列特性可以獲得對En,i的 估計(jì):

    采用這種閾值選取方法,可以直接由獲得的信號中估計(jì)閾值,不需要對噪聲分布的先驗(yàn)知識,因此這是一種自適應(yīng)的閾值選取方法,具有較好的靈活性和實(shí)用性。

    4.3 EMD消噪算法描述

    由上面分析,超寬帶信號的EMD消噪算法可以表述如下:

    (1)對含噪信號進(jìn)行EMD分解,得到IMF集合{hi(t),i=1,2,…,L};

    5 仿真與分析

    為了考察EMD消噪算法對超寬帶信號的消噪效果,本文進(jìn)行了一系列的仿真實(shí)驗(yàn)。

    實(shí)驗(yàn)一采用計(jì)算機(jī)合成的、被白噪聲污染的跳時(shí)UWB信號做為待處理的含噪信號,用0階高斯脈沖做為基本脈沖,脈沖的重復(fù)周期取10 ns,脈寬1 ns,信號長度取3個(gè)脈沖重復(fù)周期。分別在不同的信噪比條件下對EMD消噪算法的性能進(jìn)行了仿真。作為對比,同時(shí)采用小波消噪方法對相同的信號進(jìn)行了仿真。

    表1給出了不同信噪比條件下,EMD消噪和小波消噪算法的信噪比增益??梢钥吹?,使用EMD消噪算法,在不同信噪比下,平均可以獲得10 dB左右的信噪比增益;在相同信噪比條件下,可以比小波消噪算法獲得3 dB以上的信噪比增益。

    表1 消噪后信噪比增益Table 1 SNR gain after denosing

    圖5給出了信噪比為-5 dB時(shí)原始信號、含噪信號和消噪后的信號對比圖。從圖中可以看到,相比小波消噪,EMD消噪算法可以消除信號波形中由噪聲引起的毛刺,更有效地突出原始窄脈沖的輪廓和位置。

    圖4 消噪算法流程圖Fig.4 Flow chart of proposed denosing algorithm

    圖5 -5 dB輸入信噪比信號的消噪結(jié)果Fig.5 Denoising resultswhen SNR=-5 dB

    實(shí)驗(yàn)二采用了一段實(shí)際采集的含噪超寬帶脈沖信號,該信號來自于某超寬帶穿墻雷達(dá)[9]的脈沖回波,雷達(dá)發(fā)射脈沖為類高斯脈沖,脈寬1.4 ns,采樣率為20 GHz?;夭ㄖ泻卸鄠€(gè)從不同反射體反射回來的脈沖,因?yàn)榘l(fā)射功率受限,并且回波信號穿墻后受到較大衰減,使得接收信號的信噪比較低。圖6給出了含噪回波信號和采用小波消噪和EMD消噪后的信號對比圖。從圖中可以看到,采用EMD消噪算法,將附加在脈沖回波中的大量噪聲毛刺濾除,消噪后的信號與原信號相比,脈沖輪廓更加清晰,相鄰脈沖之間的界限更明顯。在原信號中,由于信噪比較低,一些低能量的脈沖回波被埋藏在噪聲水平線下,無法分辨,而EMD消噪后,這些低能量脈沖可以被清晰地分辨出來。EMD消噪算法增強(qiáng)了UWB雷達(dá)對回波脈沖的分辨能力。

    圖6 UWB穿墻雷達(dá)回波信號消噪效果Fig.6 Denoising results of UWB radar signal

    6 結(jié)論

    本文提出了一種利用EMD分解來解決超寬帶信號消噪問題的算法。該算法通過部分選取EMD分解得到的高階IMF分量重構(gòu)原始信號,從而達(dá)到消噪的目的;該算法針對各階IMF的能量分布特征,給出一種自適應(yīng)的階數(shù)閾值選擇方法,解決了IMF重構(gòu)中的閾值選擇問題。仿真結(jié)果表明,該算法在低信噪比下具有較好的消噪效果,可以有效地應(yīng)用在低信噪比超寬帶通信和雷達(dá)信號處理領(lǐng)域,具有較好的實(shí)用前景。如何降低對信號采樣率的要求和如何獲得高效的硬件實(shí)現(xiàn)可以作為該問題的進(jìn)一步研究方向。

    [1]Win M,Scholtz R.Ultra-wide bandwidth time hopping spread-spectrum impulse radio for wirelessmultiple access communications[J].IEEE Transactions on Communications,2000,48(4):679-691.

    [2]Di Benedetto,Guerino Giancola.超寬帶無線電基礎(chǔ)[M].北京:電子工業(yè)出版社,2005. Di Benedetto,Guerino Giancola.Understanding Ultra Wide Band Radio Fundamentals[M].Beijing:Publishing House of Electronics Industry,2005.(in Chinese)

    [3]Huang N E,Shen Zheng,Long SR,etal.The empiricalmode decomposition and the hilbert spectrum for nonlinear and non -stationary time series analysis[J].Proceedings of the Royal Society A,1998,454(1971):903-995.

    [4]Flandrin P,Goncalves P,Rilling G.Detrending and denoisingwith empiricalmode decompositions[C]//Proceedings of the 2004 European Signal Processing Conference.Vienna,Austria:IEEE,2004:1581-1584.

    [5]Flandrin P,Goncalves P,Rilling G.Empiricalmode decomposition as a filter bank[J].IEEE Signal Processing Letters,2004,11(2):112-114.

    [6]Boudraa A,Cexus JC.EMD-based signal filtering[J]. IEEE Transactions on Intsrument and Measurement,2007,56(6):2196-2202.

    [7]Elgamel SA,Soraghan JJ.Using EMD-FrFT filtering tomitigate very high power interference in chirp tracking radars[J]. IEEESignal Processing Letters,2011,18(4):263-266.

    [8]Wu Z,Huang NE.A study of the characteristics ofwhite noise using the empiricalmode decompositionmethod[J].Proceedings of the Royal Society A,2004,460(2046):1597-1611.

    [9]Wu Qing,Zheng Kuangyu,Li Qiao,et al.Design and implementation of sub-GHz transmitterfor ultra-wideband throughwall radar[C]//Proceedings of 2010 IEEE International Conference on Ultra-wideband.Nanjing:IEEE,2010:1-4.

    WANG Hai-liang was born in Shenyang,Liaoning Province,in 1983.He is currently working toward the Ph.D.degree.His research interests include UWB communicationsand signalprocessing.

    Email:gabbar@163.com,gabbar@ee.buaa.edu.cn

    熊華鋼(1961—),男,1998年獲博士學(xué)位,現(xiàn)為教授、博士生導(dǎo)師,主要研究方向?yàn)楹娇针娮泳W(wǎng)絡(luò)、超寬帶通信;

    吳慶(1984—),男,博士研究生,主要研究方向?yàn)槌瑢拵ㄐ藕托盘柼幚恚?/p>

    WUQingwas born in 1984.He is currentlyworking toward the Ph.D.degree.His research interests include UWB communications and signal processing.

    劉成(1985—),男,博士研究生,主要研究方向?yàn)楹娇针娮泳W(wǎng)絡(luò)、超寬帶通信。

    LIUCheng was born in 1985.He is currently working toward the Ph.D.degree.His research interests include avionics networks and UWB communications.

    A Novel Denoising Algorithm for UWB Signals Based on Em pirical M ode Decom position

    WANGHai-liang,XIONGHua-gang,WU Qing,LIU Cheng
    (School of Electronics and Information Engineering,Beijing University of Aeronautics and Astronautics,Beijing 100191,China)

    To solve the problem of denoising ultra-wideband(UWB)signals,a novel denosing algorithm based on empiricalmode decomposition(EMD)is proposed.The algorithm first uses EMD to decompose the noised signal into intrinsicmode functions(IMF),then uses high order IMFs to reconstruct the original signal.In order to get the proper threshold in IMF reconstructing,the energy distribution property of each IMF is examined by numeric simulation.According to the distribution property,a data-adaptive algorithm on choosing the threshold is proposed.Simulation results show that the proposed denoising algorithm can improve the signal signal-tonoise ratio(SNR)by 10 dB under low SNR conditions.

    ultra-wideband;signal denoising;empiricalmode decomposition

    -gangwas born in 1961.He

    the Ph.D. degree in 1998.He is now a professor and also the Ph.D.supervisor.His research interests include avionics networks and UWB communications.

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

    TN911.6

    A

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

    王海梁(1983—),男,遼寧沈陽人,博士研究生,主要研究方向?yàn)槌瑢拵ㄐ藕托盘柼幚恚?/p>

    1001-893X(2012)04-0461-05

    2011-12-20;

    2012-02-20

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

    猜你喜歡
    超寬帶階數(shù)信噪比
    基于多模諧振器的超寬帶濾波器設(shè)計(jì)
    關(guān)于無窮小階數(shù)的幾點(diǎn)注記
    確定有限級數(shù)解的階數(shù)上界的一種n階展開方法
    基于深度學(xué)習(xí)的無人機(jī)數(shù)據(jù)鏈信噪比估計(jì)算法
    超寬帶雷達(dá)建筑物結(jié)構(gòu)稀疏成像
    低信噪比下LFMCW信號調(diào)頻參數(shù)估計(jì)
    電子測試(2018年11期)2018-06-26 05:56:02
    低信噪比下基于Hough變換的前視陣列SAR稀疏三維成像
    一種基于射頻開關(guān)切換的偽隨機(jī)編碼超寬帶多發(fā)多收雷達(dá)設(shè)計(jì)
    一種脊波導(dǎo)超寬帶濾波器
    電子制作(2016年11期)2016-11-07 08:43:34
    保持信噪比的相位分解反褶積方法研究
    久久毛片免费看一区二区三区| 搡老熟女国产l中国老女人| 99久久99久久久精品蜜桃| 97精品久久久久久久久久精品| 久久精品国产a三级三级三级| 他把我摸到了高潮在线观看 | 视频区欧美日本亚洲| 五月开心婷婷网| 婷婷色av中文字幕| 美女主播在线视频| 制服诱惑二区| 国产精品久久久人人做人人爽| 日本黄色日本黄色录像| 日韩视频一区二区在线观看| 欧美日韩成人在线一区二区| 高清黄色对白视频在线免费看| 欧美激情高清一区二区三区| 亚洲精品久久成人aⅴ小说| xxxhd国产人妻xxx| 欧美精品av麻豆av| 精品一区在线观看国产| 熟女少妇亚洲综合色aaa.| 日本av免费视频播放| 99精品久久久久人妻精品| 最新在线观看一区二区三区| 亚洲欧洲日产国产| 日本av手机在线免费观看| 一级毛片精品| 国产精品久久久久成人av| 国产精品1区2区在线观看. | 亚洲精品美女久久av网站| 成在线人永久免费视频| 中国美女看黄片| 一级片'在线观看视频| 国产精品久久久久成人av| 丝袜喷水一区| 电影成人av| 桃花免费在线播放| 国产片内射在线| 久久国产精品男人的天堂亚洲| 亚洲欧美激情在线| 制服人妻中文乱码| 国产99久久九九免费精品| 不卡av一区二区三区| 97在线人人人人妻| 一个人免费在线观看的高清视频 | 91九色精品人成在线观看| 伊人久久大香线蕉亚洲五| 19禁男女啪啪无遮挡网站| 麻豆国产av国片精品| 青青草视频在线视频观看| 国产精品国产三级国产专区5o| 精品人妻在线不人妻| 电影成人av| 精品一区在线观看国产| 亚洲国产欧美日韩在线播放| 啦啦啦 在线观看视频| 两性夫妻黄色片| 精品一区二区三区av网在线观看 | 十八禁高潮呻吟视频| 国产成人精品在线电影| 亚洲伊人久久精品综合| 精品国产国语对白av| 自线自在国产av| 国产精品九九99| 99香蕉大伊视频| 亚洲色图综合在线观看| 老熟女久久久| 男女免费视频国产| 中文字幕最新亚洲高清| 又紧又爽又黄一区二区| 日韩,欧美,国产一区二区三区| 欧美人与性动交α欧美精品济南到| 黄色a级毛片大全视频| 国产无遮挡羞羞视频在线观看| 美女高潮喷水抽搐中文字幕| 极品人妻少妇av视频| av电影中文网址| 亚洲成人免费av在线播放| 少妇裸体淫交视频免费看高清 | 亚洲国产欧美日韩在线播放| 99国产综合亚洲精品| 首页视频小说图片口味搜索| 色婷婷久久久亚洲欧美| 蜜桃在线观看..| 一区二区三区精品91| 久久综合国产亚洲精品| av天堂在线播放| 精品国产一区二区三区久久久樱花| 成人国语在线视频| 久久国产精品影院| 国产精品九九99| 极品人妻少妇av视频| 夜夜骑夜夜射夜夜干| 午夜影院在线不卡| 亚洲人成电影观看| 中文字幕av电影在线播放| 久久久久久久国产电影| 免费观看av网站的网址| 成人18禁高潮啪啪吃奶动态图| 国产无遮挡羞羞视频在线观看| 在线天堂中文资源库| 国产极品粉嫩免费观看在线| 中国国产av一级| 久久久久网色| 久久中文看片网| 亚洲色图 男人天堂 中文字幕| 久久国产亚洲av麻豆专区| av天堂在线播放| 日本a在线网址| 午夜成年电影在线免费观看| 纯流量卡能插随身wifi吗| 亚洲精品粉嫩美女一区| 少妇的丰满在线观看| 男人爽女人下面视频在线观看| 国产又色又爽无遮挡免| 精品国产一区二区久久| 国产精品1区2区在线观看. | 精品人妻1区二区| 欧美xxⅹ黑人| 精品一区二区三区av网在线观看 | av欧美777| 亚洲精品久久久久久婷婷小说| 人人妻人人添人人爽欧美一区卜| av在线老鸭窝| 亚洲美女黄色视频免费看| 国产免费视频播放在线视频| 日韩电影二区| 亚洲人成电影免费在线| 下体分泌物呈黄色| 亚洲精品国产精品久久久不卡| 亚洲国产av影院在线观看| 欧美性长视频在线观看| kizo精华| 午夜精品久久久久久毛片777| 精品国产国语对白av| 精品国产一区二区三区四区第35| 成人手机av| 亚洲成人手机| 欧美午夜高清在线| bbb黄色大片| 一本综合久久免费| 99精品欧美一区二区三区四区| 最黄视频免费看| 丰满少妇做爰视频| 一区二区日韩欧美中文字幕| 亚洲九九香蕉| 国产真人三级小视频在线观看| 丰满饥渴人妻一区二区三| 国产精品一二三区在线看| 丝瓜视频免费看黄片| 欧美黄色片欧美黄色片| 亚洲av欧美aⅴ国产| 精品久久蜜臀av无| 热re99久久精品国产66热6| 国产三级黄色录像| 欧美另类一区| av线在线观看网站| 一级,二级,三级黄色视频| 色视频在线一区二区三区| 99国产极品粉嫩在线观看| 大型av网站在线播放| 久久人人爽av亚洲精品天堂| 国产日韩欧美亚洲二区| 欧美精品一区二区免费开放| 欧美人与性动交α欧美软件| 一本一本久久a久久精品综合妖精| 午夜老司机福利片| 97精品久久久久久久久久精品| www.av在线官网国产| 最新在线观看一区二区三区| 不卡一级毛片| 纵有疾风起免费观看全集完整版| 国产免费视频播放在线视频| kizo精华| 可以免费在线观看a视频的电影网站| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲国产av影院在线观看| av不卡在线播放| 亚洲成人国产一区在线观看| 久久久久视频综合| 在线观看www视频免费| 国产一区二区三区在线臀色熟女 | 捣出白浆h1v1| 国产亚洲午夜精品一区二区久久| 纯流量卡能插随身wifi吗| 久久精品国产亚洲av高清一级| 两个人看的免费小视频| 日韩电影二区| 99久久精品国产亚洲精品| 97人妻天天添夜夜摸| 黄色视频不卡| 国产片内射在线| 乱人伦中国视频| 伊人亚洲综合成人网| 啦啦啦中文免费视频观看日本| 国产亚洲av片在线观看秒播厂| 亚洲九九香蕉| 亚洲国产av新网站| 人人妻人人添人人爽欧美一区卜| 天堂8中文在线网| 青春草视频在线免费观看| 国产免费福利视频在线观看| 中文字幕最新亚洲高清| 热re99久久精品国产66热6| 男女床上黄色一级片免费看| 久久久精品国产亚洲av高清涩受| 我要看黄色一级片免费的| 亚洲少妇的诱惑av| 大型av网站在线播放| 老司机福利观看| 亚洲自偷自拍图片 自拍| tube8黄色片| 汤姆久久久久久久影院中文字幕| 美女福利国产在线| 男女无遮挡免费网站观看| 香蕉国产在线看| 50天的宝宝边吃奶边哭怎么回事| 老司机靠b影院| 丝袜美腿诱惑在线| 自拍欧美九色日韩亚洲蝌蚪91| 97精品久久久久久久久久精品| 欧美激情 高清一区二区三区| 两性夫妻黄色片| 在线观看免费高清a一片| 国产日韩欧美在线精品| 国产福利在线免费观看视频| 免费av中文字幕在线| 婷婷成人精品国产| 天天添夜夜摸| 国产精品亚洲av一区麻豆| 欧美 亚洲 国产 日韩一| 免费黄频网站在线观看国产| 狠狠婷婷综合久久久久久88av| 正在播放国产对白刺激| 黄色 视频免费看| 久久久水蜜桃国产精品网| 最黄视频免费看| 国产高清国产精品国产三级| 午夜两性在线视频| 免费av中文字幕在线| 色老头精品视频在线观看| 搡老岳熟女国产| av超薄肉色丝袜交足视频| 欧美中文综合在线视频| 男人爽女人下面视频在线观看| 久久亚洲精品不卡| 午夜福利在线观看吧| 国产日韩一区二区三区精品不卡| 激情视频va一区二区三区| 十八禁网站免费在线| 午夜老司机福利片| 一本大道久久a久久精品| 一个人免费看片子| 国产成人精品无人区| 亚洲国产中文字幕在线视频| 日韩 亚洲 欧美在线| 在线永久观看黄色视频| 久久精品成人免费网站| 超碰成人久久| 最近中文字幕2019免费版| 久久人人爽av亚洲精品天堂| 1024香蕉在线观看| 2018国产大陆天天弄谢| 人人妻,人人澡人人爽秒播| 免费看十八禁软件| 夜夜骑夜夜射夜夜干| 啦啦啦 在线观看视频| 久久久久精品国产欧美久久久 | 国产免费福利视频在线观看| 午夜福利影视在线免费观看| 久久国产精品大桥未久av| 黑人欧美特级aaaaaa片| 久久久欧美国产精品| 在线av久久热| 99久久人妻综合| 大香蕉久久成人网| 国产精品99久久99久久久不卡| 丝袜美腿诱惑在线| 午夜91福利影院| 十八禁网站网址无遮挡| 久久久久久久久免费视频了| 国产免费av片在线观看野外av| 国产一区二区在线观看av| 老司机靠b影院| 久久久久国内视频| 少妇 在线观看| 亚洲欧美一区二区三区黑人| 午夜福利影视在线免费观看| 国产av精品麻豆| 亚洲精品av麻豆狂野| 热99久久久久精品小说推荐| 免费久久久久久久精品成人欧美视频| 国产精品麻豆人妻色哟哟久久| 国产国语露脸激情在线看| 亚洲熟女毛片儿| 国产成人精品久久二区二区91| 精品国产国语对白av| 在线天堂中文资源库| 丰满迷人的少妇在线观看| 欧美黑人精品巨大| 熟女少妇亚洲综合色aaa.| 丝袜喷水一区| 一本—道久久a久久精品蜜桃钙片| 国产高清视频在线播放一区 | 精品人妻1区二区| 中文字幕人妻熟女乱码| 少妇粗大呻吟视频| 欧美人与性动交α欧美精品济南到| 一边摸一边做爽爽视频免费| 最新的欧美精品一区二区| 岛国在线观看网站| 亚洲欧美一区二区三区黑人| 国产精品熟女久久久久浪| 日韩视频一区二区在线观看| 亚洲欧美一区二区三区久久| 黄色视频,在线免费观看| 一级毛片精品| 国产亚洲欧美在线一区二区| 国产成人av激情在线播放| 亚洲精品一二三| 多毛熟女@视频| 久久国产精品男人的天堂亚洲| 操美女的视频在线观看| 国产亚洲精品第一综合不卡| 在线av久久热| 国产成人系列免费观看| 亚洲精品乱久久久久久| 搡老岳熟女国产| 最新在线观看一区二区三区| 久久久国产欧美日韩av| 可以免费在线观看a视频的电影网站| 精品亚洲乱码少妇综合久久| 日韩一卡2卡3卡4卡2021年| 男人舔女人的私密视频| 国产一区二区三区av在线| 日韩一区二区三区影片| 91国产中文字幕| 国产视频一区二区在线看| 老鸭窝网址在线观看| 色婷婷久久久亚洲欧美| 国产成人精品久久二区二区91| 国产精品久久久久成人av| 叶爱在线成人免费视频播放| 国产成+人综合+亚洲专区| 亚洲av电影在线进入| 国产av国产精品国产| 精品人妻熟女毛片av久久网站| 9色porny在线观看| 精品人妻熟女毛片av久久网站| 淫妇啪啪啪对白视频 | 午夜激情久久久久久久| 日韩精品免费视频一区二区三区| 国产精品 欧美亚洲| 在线观看免费视频网站a站| 久久久国产一区二区| tube8黄色片| 首页视频小说图片口味搜索| 日本欧美视频一区| 国产成人啪精品午夜网站| 啦啦啦免费观看视频1| av又黄又爽大尺度在线免费看| 波多野结衣av一区二区av| 欧美黄色淫秽网站| 波多野结衣av一区二区av| 首页视频小说图片口味搜索| 日韩视频在线欧美| 黑人猛操日本美女一级片| 久久久精品免费免费高清| 亚洲欧洲日产国产| 亚洲精品一二三| 亚洲国产精品一区三区| 捣出白浆h1v1| 亚洲成人手机| 国内毛片毛片毛片毛片毛片| 老司机影院成人| 黄频高清免费视频| 一边摸一边抽搐一进一出视频| 亚洲一码二码三码区别大吗| 久热爱精品视频在线9| 母亲3免费完整高清在线观看| 精品福利观看| 最近最新免费中文字幕在线| 久久久精品区二区三区| 最黄视频免费看| 男人舔女人的私密视频| 国产区一区二久久| 男人添女人高潮全过程视频| 国产男人的电影天堂91| 国产av一区二区精品久久| 欧美少妇被猛烈插入视频| 精品人妻在线不人妻| 男女国产视频网站| 精品乱码久久久久久99久播| 窝窝影院91人妻| 国产精品九九99| 99久久综合免费| 国产一区二区三区在线臀色熟女 | 成人亚洲精品一区在线观看| 人人澡人人妻人| 一本一本久久a久久精品综合妖精| 欧美中文综合在线视频| 国产亚洲av片在线观看秒播厂| 精品国产一区二区久久| 欧美精品高潮呻吟av久久| 午夜激情久久久久久久| 1024香蕉在线观看| 亚洲色图综合在线观看| 大陆偷拍与自拍| 欧美人与性动交α欧美软件| 久久久久久久久久久久大奶| 少妇 在线观看| 一区二区三区精品91| 国产一级毛片在线| 精品一区二区三区四区五区乱码| 国产日韩欧美在线精品| www.熟女人妻精品国产| 亚洲av美国av| av国产精品久久久久影院| 国产精品九九99| 男男h啪啪无遮挡| 亚洲国产av影院在线观看| 精品亚洲成国产av| 日本猛色少妇xxxxx猛交久久| 他把我摸到了高潮在线观看 | 亚洲精品中文字幕在线视频| 日韩有码中文字幕| www.自偷自拍.com| 自拍欧美九色日韩亚洲蝌蚪91| 精品人妻在线不人妻| 亚洲国产精品一区三区| 中文字幕人妻丝袜一区二区| 久久国产亚洲av麻豆专区| 最近中文字幕2019免费版| 欧美国产精品va在线观看不卡| 日韩人妻精品一区2区三区| 久久人人爽人人片av| 制服诱惑二区| 女人高潮潮喷娇喘18禁视频| 色综合欧美亚洲国产小说| 午夜福利一区二区在线看| 黄片播放在线免费| 日韩熟女老妇一区二区性免费视频| 搡老熟女国产l中国老女人| 国产伦人伦偷精品视频| 国产亚洲av高清不卡| 69精品国产乱码久久久| av天堂在线播放| 国产国语露脸激情在线看| 国产无遮挡羞羞视频在线观看| 电影成人av| 精品卡一卡二卡四卡免费| 国产欧美亚洲国产| 国产精品 欧美亚洲| 青春草视频在线免费观看| 夜夜夜夜夜久久久久| 国产高清视频在线播放一区 | 中文精品一卡2卡3卡4更新| 韩国高清视频一区二区三区| 国产成人精品久久二区二区91| 美国免费a级毛片| 91成年电影在线观看| 9191精品国产免费久久| 大型av网站在线播放| 国产精品国产三级国产专区5o| 欧美日韩一级在线毛片| 午夜福利视频在线观看免费| 亚洲欧洲精品一区二区精品久久久| 亚洲成av片中文字幕在线观看| 国产成人精品久久二区二区免费| 国产成人精品无人区| 韩国高清视频一区二区三区| 色婷婷av一区二区三区视频| 久久精品国产亚洲av香蕉五月 | 欧美另类一区| 国产三级黄色录像| 热99国产精品久久久久久7| 丝袜人妻中文字幕| 在线亚洲精品国产二区图片欧美| 亚洲第一欧美日韩一区二区三区 | 老司机影院毛片| 国产精品久久久久成人av| 国产黄色免费在线视频| 韩国精品一区二区三区| www.精华液| 国产一区二区三区综合在线观看| 我要看黄色一级片免费的| 久久久久国内视频| av线在线观看网站| 午夜成年电影在线免费观看| 岛国毛片在线播放| 男人爽女人下面视频在线观看| 久久九九热精品免费| 女人高潮潮喷娇喘18禁视频| 国产片内射在线| 日韩熟女老妇一区二区性免费视频| www日本在线高清视频| 大香蕉久久成人网| 精品亚洲成国产av| 亚洲国产精品一区三区| 精品人妻熟女毛片av久久网站| 看免费av毛片| 亚洲精品国产精品久久久不卡| 久久九九热精品免费| 亚洲欧美日韩高清在线视频 | 免费少妇av软件| 午夜激情久久久久久久| 日韩精品免费视频一区二区三区| 亚洲第一欧美日韩一区二区三区 | 又紧又爽又黄一区二区| 在线看a的网站| 夜夜骑夜夜射夜夜干| 丝袜美腿诱惑在线| 欧美人与性动交α欧美软件| 久久人妻熟女aⅴ| 两性午夜刺激爽爽歪歪视频在线观看 | 日韩免费高清中文字幕av| 香蕉丝袜av| 亚洲欧美清纯卡通| 国产亚洲精品久久久久5区| 中文字幕精品免费在线观看视频| 99国产综合亚洲精品| 亚洲自偷自拍图片 自拍| 久久精品国产亚洲av高清一级| 欧美 亚洲 国产 日韩一| 色精品久久人妻99蜜桃| 亚洲专区中文字幕在线| 国产精品欧美亚洲77777| 美女高潮到喷水免费观看| 成年人午夜在线观看视频| 国产av又大| 国产在线观看jvid| 黄色a级毛片大全视频| 叶爱在线成人免费视频播放| 国产精品一二三区在线看| 国产区一区二久久| 亚洲av美国av| 国产精品熟女久久久久浪| 亚洲天堂av无毛| 欧美97在线视频| 日日夜夜操网爽| 国产精品久久久久久人妻精品电影 | 午夜精品国产一区二区电影| 日本猛色少妇xxxxx猛交久久| av有码第一页| 亚洲专区国产一区二区| 性少妇av在线| 亚洲,欧美精品.| 9191精品国产免费久久| 免费人妻精品一区二区三区视频| 新久久久久国产一级毛片| 成人18禁高潮啪啪吃奶动态图| 久久人妻福利社区极品人妻图片| 999久久久精品免费观看国产| 国产成人欧美在线观看 | av视频免费观看在线观看| 亚洲国产欧美网| 免费高清在线观看视频在线观看| 丝袜喷水一区| 脱女人内裤的视频| 电影成人av| 精品亚洲乱码少妇综合久久| 亚洲国产精品一区三区| 亚洲欧美精品自产自拍| 亚洲精品国产av蜜桃| 男女高潮啪啪啪动态图| 亚洲人成77777在线视频| 亚洲中文字幕日韩| 日本wwww免费看| 久久中文看片网| 99国产精品一区二区三区| 69精品国产乱码久久久| 日韩欧美一区视频在线观看| 精品少妇内射三级| 十八禁人妻一区二区| 精品亚洲成国产av| 亚洲精品久久午夜乱码| 亚洲一码二码三码区别大吗| 嫩草影视91久久| 成人18禁高潮啪啪吃奶动态图| 精品免费久久久久久久清纯 | 一区二区三区激情视频| 日韩中文字幕视频在线看片| 国产在线视频一区二区| 亚洲精品一卡2卡三卡4卡5卡 | 国产欧美日韩一区二区三区在线| 欧美少妇被猛烈插入视频| 丝袜美腿诱惑在线| 欧美 日韩 精品 国产| 欧美另类一区| 少妇裸体淫交视频免费看高清 | 国产精品偷伦视频观看了| 天天影视国产精品| 午夜福利视频精品| 免费高清在线观看视频在线观看| 亚洲第一av免费看| 啦啦啦在线免费观看视频4| 久久精品成人免费网站| 国产成人系列免费观看| 真人做人爱边吃奶动态| 超色免费av| 亚洲第一av免费看| 精品少妇内射三级| 亚洲美女黄色视频免费看| 精品一品国产午夜福利视频| 天天躁日日躁夜夜躁夜夜| 免费观看a级毛片全部| 另类亚洲欧美激情| 我要看黄色一级片免费的| 18禁黄网站禁片午夜丰满| 国产免费视频播放在线视频| 一级,二级,三级黄色视频| 欧美av亚洲av综合av国产av| 少妇猛男粗大的猛烈进出视频| 一区二区三区精品91|