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

    抗沖激干擾的稀疏懲罰約束遺漏最小均方算法*

    2016-11-12 06:53:53晏國杰
    電訊技術(shù) 2016年10期
    關(guān)鍵詞:遺漏均方高斯

    晏國杰,林 云

    (重慶郵電大學(xué) 移動通信技術(shù)重慶市重點實驗室,重慶 400065)

    ?

    抗沖激干擾的稀疏懲罰約束遺漏最小均方算法*

    晏國杰**,林云

    (重慶郵電大學(xué) 移動通信技術(shù)重慶市重點實驗室,重慶 400065)

    當(dāng)被識別系統(tǒng)是稀疏系統(tǒng)時,傳統(tǒng)的遺漏最小均方(LLMS)自適應(yīng)算法收斂性能較差,特別在非高斯噪聲環(huán)境中,該算法性能進(jìn)一步惡化甚至算法不平穩(wěn)收斂。為了解決因信道的稀疏性使算法收斂變慢的問題,對LLMS算法的代價函數(shù)分別利用加權(quán)1-norm和加權(quán)零吸引兩種稀疏懲罰項進(jìn)行改進(jìn);為了優(yōu)化算法的抗沖激干擾的性能,利用符號函數(shù)對已改進(jìn)的算法迭代式作進(jìn)一步改進(jìn)。同時,將提出的兩個算法運用于非高斯噪聲環(huán)境下的稀疏系統(tǒng)識別,仿真結(jié)果顯示提出的算法性能優(yōu)于現(xiàn)存的同類稀疏算法。

    稀疏系統(tǒng)識別;自適應(yīng)算法;沖激干擾;收斂性

    1 引 言

    最小均方(Least-Mean Square,LMS)自適應(yīng)算法[1]由于具有結(jié)構(gòu)簡單和計算復(fù)雜度低的特點[2],已經(jīng)被廣泛應(yīng)用在各工程領(lǐng)域,如控制工程、通信工程、雷達(dá)信號處理、干擾消除、有噪聲源控制(Active Noise Control,ANC)、生物醫(yī)學(xué)工程、波束成形、語音回聲信道等[3-4]。而當(dāng)輸入信號是特征分散較大的時,其傳統(tǒng)的LMS自適應(yīng)算法收斂較慢。為了解決這個問題,Guitlin等利用輸入信號特征分散較大這一特點在文獻(xiàn)[5]中提出了改進(jìn)的LMS算法,即遺漏的最小均方(Leaky Least-Mean Square,LLMS)自適應(yīng)算法。在輸入信號的特征分散度較大情況下,LLMS算法的確解決了LMS算法收斂慢的缺陷。

    在實際的工程應(yīng)用中遇到的系統(tǒng)并不總是分散的(稠密的),而有可能是稀疏系統(tǒng)(大部分抽頭系數(shù)為零或其絕對值小到可以忽略的,只有很少部分活躍抽頭的系數(shù)較大)。例如:在稀疏的網(wǎng)絡(luò)回聲消除系統(tǒng)中,由于時延效應(yīng)的影響,在總時隙為64~128 ms的響應(yīng)中,其活躍的時隙大概只有8~12 ms[6]。當(dāng)系統(tǒng)是稀疏的時候,LLMS算法收斂性較差。為了解決這個問題,應(yīng)該研究稀疏的LLMS算法。稀疏的LLMS算法目前研究相對較少,現(xiàn)有的稀疏的LLMS算法主要有Zero-Attracting Leaky Least-Mean-Square(ZA-LLMS)[7]和Weighted Zero-Attracting Leaky-LMS(WZA-LLMS)[8]。在噪聲為高斯干擾的情況下,以上兩種算法雖然增加了計算抽頭系數(shù)矢量的一范數(shù)的計算復(fù)雜度,但是在一定程度上解決了當(dāng)系統(tǒng)是稀疏的時候LLMS算法收斂慢的問題。在實際的工程應(yīng)用中,非高斯沖激干擾(包括沖激干擾、沖激干擾和高斯干擾的混合)也是經(jīng)常遇到的。在非高斯干擾的環(huán)境中,以上存在的稀疏類的LLMS算法的性能嚴(yán)重惡化,甚至使算法不能收斂。

    本文為了解決稀疏抗非高斯沖激干擾的問題,提出了兩個改進(jìn)的稀疏LLMS算法,即加權(quán)零吸引的符號遺漏最小均方算法(Reweighted Zero-Attracting Sign Least Mean Square,RZASLLMS)和基于加權(quán)一范數(shù)的符號遺漏最小均方算法(Reweighted1-norm Sign Least Mean Square,RL1-SLLMS)。首先,為了利用系統(tǒng)稀疏性這一先驗知識,在傳統(tǒng)的LLMS算法的代價函數(shù)的右邊分別加上稀疏約束懲罰函數(shù)項(即加權(quán)的零吸引項和加權(quán)的抽頭系數(shù)矢量的一范數(shù)項),從而使傳統(tǒng)的LLMS算法在識別稀疏系統(tǒng)時的收斂速度得到提高。其次,由于符號函數(shù)能防止誤差因沖激干擾突然劇增而具有抗沖激噪聲的性能,被廣泛應(yīng)用于各種自適應(yīng)算法中[9-10]。受符號函數(shù)這一特性的啟發(fā),對上一步改進(jìn)得到的兩個迭代方程做進(jìn)一步修正:分別將兩個迭代方程中的誤差項用自變量為誤差的符號函數(shù)代替,以此降低因誤差突增而帶來算法性能惡化的程度。經(jīng)過以上兩步改進(jìn)之后就可以分別得到提出的兩個算法的迭代方程了,并對提出的算法給出相應(yīng)的計算復(fù)雜度分析。同時,將提出的兩個算法運用于稀疏系統(tǒng)識別的Matlab仿真并且與現(xiàn)有的遺漏類算法對比。由于WZA-LLMS的算法性能優(yōu)于ZA-LLMS的性能,實驗只將對WZA-LLMS算法、傳統(tǒng)的LLMS算法與提出的算法進(jìn)行對比。

    2 傳統(tǒng)的遺漏最小均方自適應(yīng)算法

    2.1系統(tǒng)模型

    首先定義一個線性系統(tǒng)識別模型:

    d(n)=xT(n)wo+v(n)。

    (1)

    式中:d(n)表示期望信號;列矢量x(n)=[x(n),x(n-1),…,x(n-L+1)]T為輸入信號;wo為待估計系統(tǒng)的抽頭系數(shù)列矢量;v(n)是背景加性噪聲;L是濾波器長度。濾波器的輸出可以表示為y(n)=xT(n)w(n),估計誤差可表示為

    e(n)=d(n)-xT(n)w(n) 。

    (2)

    式中:w(n)=[w0,w2,…,wL-1]T是wo在第n次迭代的估計列矢量。

    2.2遺漏最小均方(LLMS)自適應(yīng)算法

    傳統(tǒng)的LLMS算法的代價函數(shù)定義如下[11-13]:

    (3)

    式中:γ>0為遺漏因子。

    由負(fù)梯度方向可推導(dǎo)其迭代式如下:

    (1-μγ)w(n)+μe(n)x(n)。

    (4)

    式中:2>μ>0為算法迭代步長。

    3 符號最小均方自適應(yīng)算法SLMS

    SLMS算法作為一種抗沖激干擾的算法,其表達(dá)式為[14]

    w(n+1)=w(n)+μsgn(e(n))x(n)。

    (5)

    式中:sgn(·)表示符號函數(shù)。

    3.1提出的RZASLLMS算法

    源于Reweighted Zero-Attracting LMS(RZA-LMS)[15]算法的思想,為了使傳統(tǒng)的LLMS算法具有稀疏識別性,在傳統(tǒng)遺漏的最小均方算法的代價函數(shù)的基礎(chǔ)上增加稀疏約束懲罰函數(shù)項,即

    JRZA(n)=J(n)+λRZAJRZA1(n)=

    (6)

    式中:λRZA>0是一個為了調(diào)節(jié)估計誤差與稀疏控制項的權(quán)衡因子。再由負(fù)梯度方向可得

    (7)

    為了減小因沖激噪聲而使誤差發(fā)生突變給算法收斂性帶來的不利影響,利用符號函數(shù)對式(7)進(jìn)一步改進(jìn)可得RZASLLMS算法的表達(dá)式:

    w(n+1)=(1-μγ)w(n)+

    (8)

    式中:ρ=μλRZA/ε′,ε=1/ε′,一般ε=20[14]。

    3.2提出的RL1-SLLMS算法

    由文獻(xiàn)[14]知加權(quán)的一范數(shù)在估計稀疏信道方面有較好的性能,因此,為了增加LLMS算法的稀疏識別性,可以用加權(quán)的一范數(shù)改進(jìn)LLMS算法的代價函數(shù)為

    λRL1‖f(n)w(n)‖1。

    (9)

    式中:λRL1>0與λRZA一樣是權(quán)衡因子;符號‖·‖1表示向量的一范數(shù);而f(n)被定義[14]為

    (10)

    式中:1>δRL1>0,例如δRL1=0.05,因此,[f(n)]i>0滿足i=0,1,…,L-1。

    由式(9)負(fù)梯度方向迭代有

    (1-μγ)w(n)+μe(n)x(n)-

    (11)

    為了使算法具有抗沖激干擾性,對式(11)進(jìn)一步改進(jìn)如下:

    w(n+1)=(1-μγ)w(n)+μsgn(e(n))x(n)-

    (12)

    式中:δRL1>0防止算法初始迭代時除零。

    4 提出算法的計算復(fù)雜度分析

    表1是對LLMS、WZALLMS、RZASLLMS和RL1-SLLMS算法在一次迭代中計算復(fù)雜度的對比分析,其中L表示濾波器長度。在統(tǒng)計運算次數(shù)時,將加法和減法的計算次數(shù)歸納在一起,乘法和除法的計算次數(shù)歸納在一起,符號函數(shù)的計算次數(shù)單獨歸納為一項。從表中可以看出傳統(tǒng)的LLMS算法計算復(fù)雜度是最小的,WZALLMS、RZASLLMS和RL1-SLLMS算法由于都增加了稀疏約束懲罰項,所以都在LLMS算法的基礎(chǔ)上增加了2L次加法運算和2L次乘法運算。WZALLMS、RZASLLMS和RL1-SLLMS算法的計算復(fù)雜度保持基本相同,其原因是這3個算法都是在LLMS算法的代價函數(shù)上增加類似的稀疏約束項進(jìn)行處理的。而對于提出的RZASLLMS和RL1-SLLMS算法之所以比WZALLMS算法多一次符號函數(shù)的運算,是因為提出的算法為了抗沖激干擾對測量誤差使用了符號函數(shù)。

    表1 LLMS、WZALLMS、RZASLLMS和RL1-SLLMS算法在一次迭代中的計算復(fù)雜度的分析

    總之,提出的RZASLLMS和RL1-SLLMS算法雖然與傳統(tǒng)的LLMS算法相比計算復(fù)雜度相對有些增加,但是與WZALLMS算法相比每一次迭代只增加了一次符號函數(shù)的計算。

    5 仿真結(jié)果

    下面對提出的算法做仿真實驗分析。由于大部分真實環(huán)境的噪聲是高斯干擾,而對于另一些情況,如電子設(shè)備中的人為噪聲、自然現(xiàn)象中的照明高峰等可以描述為非高斯干擾模型[16]。所以實驗分為兩個部分:第一部分主要驗證提出的算法在一般的高斯干擾環(huán)境下的稀疏自適應(yīng)系統(tǒng)識別有良好的均方偏差性能,在被識別系統(tǒng)稀疏程度不同情況下,與現(xiàn)有算法對比;第二部分主要是將噪聲換成沖激噪聲,與現(xiàn)有算法的對比,同時進(jìn)一步在語音信道中進(jìn)行驗證。在整個仿真實驗過程中,性能學(xué)習(xí)追蹤的歸一化均方誤差(Normalized Mean Square Deviation,NMSD)表達(dá)式為

    所有的仿真實驗的輸入信號都是零均值的隨機序列通過x(n)=0.9x(n-1)+b(n)產(chǎn)生的,b(n)是零均值、方差為1的高斯噪聲。為了避免實驗結(jié)果的偶然性,所有實驗結(jié)果都是計算機重復(fù)100次獨立仿真得到的平均結(jié)果。

    5.1在高斯干擾下的性能

    在輸入信號中加入10 dB的零均值的高斯噪聲。輸入信號通過一個抽頭系數(shù)為32個的稀疏濾波器,其中2個抽頭系數(shù)為1,其他的均為0,并定義其稀疏度為2/32。仿真結(jié)果如圖1所示。所有算法的步長取μ=0.002,遺漏因子取λ=0.000 1[8],λRZA=0.02,ε=20,δRL1=0.05,λRL1=0.000 05[14]。從圖1中可以看出,提出的RZASLLMS和RL1-SLLMS算法歸一化均方偏差曲線比傳統(tǒng)的LLMS算法低,而與WZALLMS算法的歸一化學(xué)習(xí)曲線相差不大。當(dāng)?shù)? 000次時,系統(tǒng)由wo變到-wo,提出的算法和WZALLMS算法依然保持良好的健壯性。因此,本文提出的兩個改進(jìn)算法優(yōu)于LLMS算法,和已存在的WZALLMS有相同的稀疏識別性。

    圖1 在高斯干擾環(huán)境下4種算法的學(xué)習(xí)曲線比較

    為了進(jìn)一步證實本文提出的算法具有稀疏識別性,在環(huán)境噪聲與上面實驗保持一致的情況下,比較在不同稀疏度(2/32,4/32,6/32,8/32,10/32)算法的性能。如圖2,各算的平穩(wěn)時的NMSD的取值均是在保持相同的收斂速度的前提。從圖2中可以看到WZALLMS和提出的RZASLLMS、RL1-LLMS算法的歸一化的均方偏差隨著稀疏度的增加逐漸向LLMS算法靠攏,這是因為隨著稀疏度增加,系統(tǒng)越來越分散,稀疏的自適應(yīng)算法WZALLMS、 RZASLLMS、RL1-LLMS的收斂性就越來越接近LLMS算法了。除此之外,LLMS算法的曲線基本保持水平,這是因為其步長一直沒有改變。由于符號算法自身有收斂慢的特點[9-10],所以正如圖1和圖2所示,提出的算法的均方偏差稍微高于WZALLMS算法。雖然RZASLLMS和RL1-SLLMS算法與WZALLMS算法相比犧牲了少于1 dB的均方偏差,但是與之前的LLMS算法相比具有良好的收斂性。同時,正如接下來的實驗要分析的,提出的算法擁有以上算法不具有的抗沖激性能。

    圖2 不同稀疏度下的各算法的平穩(wěn)時的NMSD學(xué)習(xí)曲線

    5.2非高斯干擾下的性能

    圖3 在非高斯干擾下各算法的學(xué)習(xí)曲線對比

    為了進(jìn)一步驗證提出的算法在稀疏系統(tǒng)識別的實用性,下面將實驗信道換成如圖4所示的隨機產(chǎn)生稀疏語音回聲信道。利用這個信道做兩次仿真,一次仿真的噪聲是高斯的,一次為非高斯干擾,其結(jié)果分別如圖5和圖6所示。

    圖4 隨機產(chǎn)生抽頭數(shù)為128個的稀疏語音信道

    圖5 高斯干擾環(huán)境下的語音信道識別

    圖6 非高斯干擾環(huán)境下的語音信道識別

    從圖5可以知道,在稀疏度較高的語音信道中提出的RZASLLMS、RL1-LLMS算法比現(xiàn)有的LLMS和 WZALLMS算法收斂更快。同時,從圖6更可以看出提出的算法在非高斯干擾下的語音信道識別中有更強的健壯性。

    6 結(jié)束語

    本文針對傳統(tǒng)的LLMS算法在識別稀疏系統(tǒng)時收斂速度慢的問題和不能抗沖激干擾的問題,首先利用兩種不同的一范數(shù)加權(quán)因子對傳統(tǒng)的LLMS算法的代價函數(shù)進(jìn)行改進(jìn),使傳統(tǒng)的LLMS算法具有稀疏識別性;然后又利用符號算法的思想對算法的迭代方程作進(jìn)一步改進(jìn),使提出的算法具有抗沖激干擾的性能。由仿真結(jié)果得知,提出的兩個算法的收斂性不僅在高斯干擾環(huán)境下明顯優(yōu)于LLMS算法,而且在噪聲為非高斯沖激干擾環(huán)境中也明顯優(yōu)于LLMS算法和WZALLMS算法的收斂性能。因此,提出的RZASLLMS、RL1-SLLMS算法雖然犧牲了很小的計算復(fù)雜度,但是解決了LLMS算法不能識別稀疏系統(tǒng)的問題和WZALLMS算法不能抗沖激干擾的問題,具有很好的工程意義。下一步的研究工作是將本文提出的算法運用到?jīng)_激干擾環(huán)境下稀疏度時變的自適應(yīng)系統(tǒng)識別問題中。

    [1]HAYKIN S.Adaptive filter theory[J].Englewood Clilfs,NJ:Prentice Hall Co.,1991:231-238.

    [2]葉挺,朱賽男.新變步長自適應(yīng)算法及其時延估計性能[J].電訊技術(shù),2014,54(6):775-779.

    YE Ting,ZHU Sainan.A new variable step-size LMS adaptive filtering algorithm and its performance in time delay estimation[J].Telecommunication Engineering,2014,54(6):775-779.(in Chinese)

    [3]FARHANG-BOROUJENY B.Adaptive filters:theory and applications[M].New York:John Wiley & Sons,2013:9-20.

    [4]KADAGE M R,PATIL M S N.A new LMS for beamforming in mobile communication[J].International Research Journal of Engineering and Technology,2015,2(5):1-5.

    [5]GITLIN R D,MEADORS H C,WEINSTEIN S B. The tap-leakage algorithm:an algorithm for the stable operation of a digitally implemented,fractionally spaced adaptive equalizer[J].Bell System Technical Journal,1982,61(8):1817-1839.

    [6]DAS B K,CHAKRABORTY M. Sparse adaptive filtering by an adaptive convex combination of the LMS and the ZA-LMS algorithms[J].IEEE Transactions on Circuits and Systems I:Regular Papers,2014,61(5):1499-1507.

    [7]SALMAN M S. Sparse leaky-LMS algorithm for system identification and its convergence analysis[J].International Journal of Adaptive Control and Signal Processing,2014,28(10):1065-1072.

    [8]SALMAN M S,JAHROMI M N S,HOCANIN A,et al.A weighted zero-attracting leaky-LMS algorithm[C]//Proceedings of 2012 20th International Conference on Software,Telecommunications and Computer Networks(SoftCOM).Split,Croatia:IEEE,2012:1-3.

    [9]PALEOLOGU C,BENESTY J,CIOCHINA S. Sparse adaptive filters for echo cancellation[J].Synthesis Lectures on Speech and Audio Processing,2010,6(1):34-36.

    [10]SEO J H,JUNG S M,PARK P G. Proportionate sign subband adpative filtering algorithm for network echo cancellers[J].Journal of Advances in Computer Networks,2015,3(2):99-105.

    [11]MAYYAS K A,ABOULNASR T. Leaky lms:a detailed analysis[C]//Proceedings of 1995 IEEE International Symposium on Circuits and Systems.Seatlle,Washington,USA:IEEE,1995:1255-1258.

    [12]KAMENETSKY M,WIDROW B.A variable leaky LMS adaptive algorithm[C]//Proceedings of the Thirty-Eighth Asilomar Conference on Signals,Systems and Computers.Pacific Grove,CA:IEEE,2004:125-128.

    [13]MAYYAS K,ABOULNASR T. Leaky LMS algorithm:MSE analysis for Gaussian data[J].IEEE Transactions on Signal Processing,1997,45(4):927-934.

    [14]GUAN G,XU L,MA W,et al.Robust adaptive sparse channel estimation in the presence of impulsive noises[C]//Proceedings of 2015 IEEE International Conference on Digital Signal Processing(DSP).Singapore:IEEE,2015:628-632.

    [15]CHEN Y,GU Y,HERO III A O. Sparse LMS for system identification[C]//Proceedings of 2009 IEEE International Conference on Acoustics,Speech and Signal Processing.Taipei:IEEE,2009:3125-3128.

    [16]WU Z,PENG S,MA W,et al.Minimum error entropy algorithms with sparsity penalty constraints[J].Entropy,2015,17(5):3419-3437.

    [17]SHI L,LIN Y.Convex combination of adaptive filters under the maximum correntropy criterion in impulsive interference[J].IEEE Signal Processing Letters,2014,21(11):1385-1388.

    晏國杰(1989—),男,貴州銅仁人,現(xiàn)為碩士研究生,主要研究方向為稀疏自適應(yīng)信號處理;

    YAN Guojie was born in Tongren,Guizhou Province,in 1989. He is now a graduate student. His research concerns sparse signal processing.

    Email:guojieyancqupt@foxmail.com

    林云(1968—),男,四川南充人,現(xiàn)為副教授、碩士生導(dǎo)師,主要研究方向為壓縮感知技術(shù)及稀疏信號處理。

    LIN Yun was born in Nanchong,Sichuan Province,in 1968. He is now an associate professor and also the instructor of graduate students.His research concerns compressing sense and sparse signal processing.

    Email:lycqupt@sina.com

    Sparse Penalty Constraint Leaky Least Mean Square Algorithms against Impulsive Interference

    YAN Guojie,LIN Yun

    (Chongqing Key Laboratory of Mobile Communications Technology,Chongqing University of Posts and Telecommunications,Chongqing 400065,China)

    The leaky least mean square(LLMS) adaptive filtering algorithm converges slowly when the identified system is sparse. Especially when the noise is non-Gaussian impulsive interference,the performance of LLMS algorithm deteriorates severely.To solve the problem that the convergent rate becomes slower because the system is sparse,the cost function of the conventional LLMS algorithm is improved by the two penalty functions,the reweighted zero-attracting(the log-sum penalty) and reweighted1-norm(RL1),respectively. To address the problem of the impulsive interference,the iterative functions are improved by introducing the sign function. Simultaneously,the simulations are made for the proposed algorithms to prove to be better performances compared with existing leaky-style algorithms in the case of impulsive interference.

    sparse system identification;adaptive filtering algorithm;impulsive interference;convergence

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

    2016-01-15;

    2016-04-18Received date:2016-01-15;Revised date:2016-04-18

    TN912.3

    A

    1001-893X(2016)10-1153-06

    引用格式:晏國杰,林云.抗沖激干擾的稀疏懲罰約束遺漏最小均方算法[J].電訊技術(shù),2016,56(10):1153-1158.[YAN Guojie,LIN Yun.Sparse penalty constraint leaky least mean square algorithms against impulsive interference[J].Telecommunication Engineering,2016,56(10):1153-1158.]

    **通信作者:guojieyancqupt@foxmail.comCorresponding author:guojieyancqupt@foxmail.com

    猜你喜歡
    遺漏均方高斯
    小高斯的大發(fā)現(xiàn)
    來自動物星球的挑戰(zhàn)(二)小五狼遺漏的線索
    一類隨機積分微分方程的均方漸近概周期解
    遺漏的光陰
    鴨綠江(2021年17期)2021-11-11 13:03:41
    Beidou, le système de navigation par satellite compatible et interopérable
    天才數(shù)學(xué)家——高斯
    應(yīng)用品管圈降低腹腔鏡抗反流手術(shù)術(shù)前準(zhǔn)備遺漏率的實踐
    基于抗差最小均方估計的輸電線路參數(shù)辨識
    基于隨機牽制控制的復(fù)雜網(wǎng)絡(luò)均方簇同步
    有限域上高斯正規(guī)基的一個注記
    91麻豆精品激情在线观看国产 | 99久久综合免费| 老司机在亚洲福利影院| 精品国产一区二区三区四区第35| 狠狠婷婷综合久久久久久88av| 男人操女人黄网站| 日日夜夜操网爽| 亚洲欧美精品综合一区二区三区| 久久精品熟女亚洲av麻豆精品| 亚洲九九香蕉| 亚洲欧美成人综合另类久久久| 国产视频首页在线观看| 国产精品一二三区在线看| 久热这里只有精品99| 国产免费现黄频在线看| 性少妇av在线| 又黄又粗又硬又大视频| 久久国产精品男人的天堂亚洲| 精品久久久久久久毛片微露脸 | 人人妻人人爽人人添夜夜欢视频| 欧美精品一区二区免费开放| 欧美少妇被猛烈插入视频| 国产精品熟女久久久久浪| 女人高潮潮喷娇喘18禁视频| 狠狠婷婷综合久久久久久88av| 操美女的视频在线观看| 超碰成人久久| 亚洲国产精品国产精品| 免费女性裸体啪啪无遮挡网站| 亚洲av成人精品一二三区| 亚洲 欧美一区二区三区| 亚洲精品在线美女| 黄色视频不卡| 麻豆国产av国片精品| 国产精品一区二区免费欧美 | 91麻豆精品激情在线观看国产 | 在线观看免费视频网站a站| 国产av一区二区精品久久| 老汉色av国产亚洲站长工具| 国产免费现黄频在线看| 电影成人av| 亚洲精品日韩在线中文字幕| 日本色播在线视频| 18禁黄网站禁片午夜丰满| av片东京热男人的天堂| 成人三级做爰电影| 亚洲欧洲精品一区二区精品久久久| 在线观看免费高清a一片| 9热在线视频观看99| 亚洲欧美色中文字幕在线| 你懂的网址亚洲精品在线观看| e午夜精品久久久久久久| 王馨瑶露胸无遮挡在线观看| 国产在线观看jvid| 免费女性裸体啪啪无遮挡网站| 亚洲成国产人片在线观看| 精品卡一卡二卡四卡免费| 日韩制服丝袜自拍偷拍| 中文字幕另类日韩欧美亚洲嫩草| 夫妻午夜视频| 在线观看www视频免费| 亚洲精品一二三| 中国美女看黄片| 免费av中文字幕在线| 下体分泌物呈黄色| 亚洲,欧美精品.| 超色免费av| 亚洲情色 制服丝袜| 一区二区三区激情视频| 国产在线免费精品| 欧美黄色淫秽网站| 丝袜脚勾引网站| 美女中出高潮动态图| 99国产精品一区二区蜜桃av | 欧美人与善性xxx| 午夜av观看不卡| 我的亚洲天堂| www.自偷自拍.com| 男女之事视频高清在线观看 | 自线自在国产av| 老司机影院毛片| av片东京热男人的天堂| 97人妻天天添夜夜摸| 久热这里只有精品99| 视频在线观看一区二区三区| 日本av手机在线免费观看| 在线观看www视频免费| 一二三四社区在线视频社区8| 亚洲熟女精品中文字幕| 精品卡一卡二卡四卡免费| 91成人精品电影| 中文字幕另类日韩欧美亚洲嫩草| 亚洲欧洲精品一区二区精品久久久| 日韩中文字幕欧美一区二区 | 亚洲欧美清纯卡通| 午夜影院在线不卡| 青草久久国产| 在线观看免费高清a一片| 国产日韩欧美亚洲二区| 国产精品 国内视频| 丰满人妻熟妇乱又伦精品不卡| 一级黄色大片毛片| 看免费成人av毛片| 精品一区二区三区四区五区乱码 | 大香蕉久久网| 各种免费的搞黄视频| 老司机亚洲免费影院| 一区二区日韩欧美中文字幕| 久9热在线精品视频| 丝袜美足系列| 国产淫语在线视频| 国产精品偷伦视频观看了| 搡老乐熟女国产| 中文乱码字字幕精品一区二区三区| 免费在线观看影片大全网站 | 99精国产麻豆久久婷婷| 亚洲av欧美aⅴ国产| 黑人猛操日本美女一级片| 一边摸一边做爽爽视频免费| 日韩,欧美,国产一区二区三区| 精品国产乱码久久久久久小说| 欧美日韩国产mv在线观看视频| 搡老岳熟女国产| 男女边摸边吃奶| 日本av手机在线免费观看| 久久人人97超碰香蕉20202| 99国产综合亚洲精品| 91麻豆av在线| 国产成人欧美| 精品免费久久久久久久清纯 | 精品人妻熟女毛片av久久网站| 欧美日韩福利视频一区二区| 欧美国产精品va在线观看不卡| 日韩电影二区| 99香蕉大伊视频| 亚洲国产欧美网| 欧美97在线视频| 亚洲精品自拍成人| 美女福利国产在线| 欧美黑人精品巨大| 日本一区二区免费在线视频| 高清欧美精品videossex| 午夜免费观看性视频| 国产精品偷伦视频观看了| 真人做人爱边吃奶动态| 各种免费的搞黄视频| 亚洲欧美精品综合一区二区三区| 色94色欧美一区二区| 国产精品人妻久久久影院| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲九九香蕉| 国产精品一区二区在线不卡| 亚洲国产欧美日韩在线播放| 精品国产一区二区三区久久久樱花| av片东京热男人的天堂| 亚洲专区国产一区二区| 国产成人精品在线电影| av一本久久久久| 飞空精品影院首页| av网站免费在线观看视频| 黑人欧美特级aaaaaa片| 免费少妇av软件| 亚洲精品久久久久久婷婷小说| 亚洲精品av麻豆狂野| 丝袜喷水一区| 日本色播在线视频| 9热在线视频观看99| 亚洲精品美女久久久久99蜜臀 | 亚洲国产毛片av蜜桃av| 欧美日韩视频精品一区| www日本在线高清视频| 国产欧美日韩一区二区三 | 狂野欧美激情性xxxx| 久热爱精品视频在线9| 国产精品九九99| 日韩制服骚丝袜av| 黄色毛片三级朝国网站| 又紧又爽又黄一区二区| 久久精品国产a三级三级三级| 一区二区日韩欧美中文字幕| 亚洲男人天堂网一区| 成年美女黄网站色视频大全免费| 久久精品亚洲av国产电影网| 少妇的丰满在线观看| 亚洲色图综合在线观看| 亚洲精品第二区| av天堂在线播放| 99re6热这里在线精品视频| 久久久亚洲精品成人影院| 最新在线观看一区二区三区 | 久久人人爽人人片av| 一区福利在线观看| 在线看a的网站| 十分钟在线观看高清视频www| 亚洲国产精品成人久久小说| 久久精品人人爽人人爽视色| 亚洲情色 制服丝袜| 国产日韩欧美亚洲二区| 久9热在线精品视频| 99国产精品一区二区三区| 高清视频免费观看一区二区| 欧美日韩亚洲国产一区二区在线观看 | 十八禁高潮呻吟视频| 美女大奶头黄色视频| 国产99久久九九免费精品| 午夜影院在线不卡| 国产成人av教育| 黄色视频在线播放观看不卡| 五月开心婷婷网| 国产男人的电影天堂91| 久久人人97超碰香蕉20202| 电影成人av| 美女福利国产在线| 久久天堂一区二区三区四区| 自拍欧美九色日韩亚洲蝌蚪91| 国产成人欧美| 亚洲精品久久久久久婷婷小说| 午夜91福利影院| 国产精品久久久人人做人人爽| 女警被强在线播放| 国产精品一二三区在线看| 丝袜美腿诱惑在线| 少妇猛男粗大的猛烈进出视频| 日韩精品免费视频一区二区三区| 又大又黄又爽视频免费| 韩国高清视频一区二区三区| 交换朋友夫妻互换小说| 亚洲欧美一区二区三区久久| 深夜精品福利| 久久久久精品国产欧美久久久 | 亚洲精品久久成人aⅴ小说| 最黄视频免费看| 九色亚洲精品在线播放| 90打野战视频偷拍视频| 亚洲欧美日韩高清在线视频 | 国产深夜福利视频在线观看| 18在线观看网站| 女人爽到高潮嗷嗷叫在线视频| 久久亚洲国产成人精品v| 女人精品久久久久毛片| 成年av动漫网址| 又大又爽又粗| 欧美日韩精品网址| 国产男女超爽视频在线观看| 亚洲精品av麻豆狂野| 考比视频在线观看| 18禁裸乳无遮挡动漫免费视频| 热re99久久精品国产66热6| 满18在线观看网站| 在线观看免费高清a一片| 午夜免费男女啪啪视频观看| 美女扒开内裤让男人捅视频| tube8黄色片| 久9热在线精品视频| 水蜜桃什么品种好| 大话2 男鬼变身卡| av电影中文网址| 久久精品国产亚洲av涩爱| 午夜激情久久久久久久| 在线观看人妻少妇| 青草久久国产| 国产日韩欧美视频二区| 亚洲国产日韩一区二区| 操美女的视频在线观看| 天天影视国产精品| 亚洲精品日韩在线中文字幕| 黄片播放在线免费| 每晚都被弄得嗷嗷叫到高潮| 波野结衣二区三区在线| 国产男人的电影天堂91| 黄色视频在线播放观看不卡| 免费在线观看完整版高清| 99精品久久久久人妻精品| 国产免费一区二区三区四区乱码| 午夜福利一区二区在线看| 七月丁香在线播放| 欧美日韩精品网址| 精品人妻熟女毛片av久久网站| 麻豆国产av国片精品| 精品高清国产在线一区| 国产亚洲av片在线观看秒播厂| 亚洲久久久国产精品| 只有这里有精品99| 国产精品 欧美亚洲| 国产精品久久久久成人av| 欧美+亚洲+日韩+国产| 亚洲欧美一区二区三区黑人| 国产成人欧美在线观看 | 超碰成人久久| 精品视频人人做人人爽| 亚洲一码二码三码区别大吗| 国产片内射在线| 韩国高清视频一区二区三区| 最近最新中文字幕大全免费视频 | 亚洲av欧美aⅴ国产| 一本久久精品| 大陆偷拍与自拍| 一级毛片女人18水好多 | 精品卡一卡二卡四卡免费| 亚洲精品美女久久久久99蜜臀 | 90打野战视频偷拍视频| 亚洲中文日韩欧美视频| 成人免费观看视频高清| 日日夜夜操网爽| 欧美黄色片欧美黄色片| 9色porny在线观看| 欧美黑人欧美精品刺激| 9色porny在线观看| 男女床上黄色一级片免费看| 黑人欧美特级aaaaaa片| 人人妻人人爽人人添夜夜欢视频| 久久综合国产亚洲精品| 欧美人与善性xxx| 久久久久精品国产欧美久久久 | 日本黄色日本黄色录像| 欧美久久黑人一区二区| 亚洲第一青青草原| 在线观看国产h片| 两个人免费观看高清视频| 久久久久视频综合| 欧美精品高潮呻吟av久久| 色婷婷av一区二区三区视频| 亚洲专区中文字幕在线| 操出白浆在线播放| 天天躁夜夜躁狠狠躁躁| 国产精品秋霞免费鲁丝片| 婷婷成人精品国产| 亚洲一卡2卡3卡4卡5卡精品中文| 美女视频免费永久观看网站| 亚洲天堂av无毛| 啦啦啦啦在线视频资源| 18在线观看网站| 亚洲欧美精品综合一区二区三区| 麻豆av在线久日| a 毛片基地| 永久免费av网站大全| 高清欧美精品videossex| 大陆偷拍与自拍| 午夜日韩欧美国产| 叶爱在线成人免费视频播放| 91九色精品人成在线观看| 成人免费观看视频高清| 王馨瑶露胸无遮挡在线观看| 成人亚洲精品一区在线观看| netflix在线观看网站| 大香蕉久久网| 亚洲成人手机| 夫妻性生交免费视频一级片| 中文字幕高清在线视频| 午夜福利在线免费观看网站| 亚洲国产最新在线播放| 国产伦人伦偷精品视频| 人妻人人澡人人爽人人| 国产成人精品久久久久久| 亚洲精品国产av成人精品| 国产精品av久久久久免费| 国产成人a∨麻豆精品| 欧美激情极品国产一区二区三区| 国产成人a∨麻豆精品| 五月开心婷婷网| 日日爽夜夜爽网站| 亚洲人成77777在线视频| 91精品伊人久久大香线蕉| 国产在线免费精品| av国产久精品久网站免费入址| 日韩熟女老妇一区二区性免费视频| 啦啦啦视频在线资源免费观看| 岛国毛片在线播放| 大片免费播放器 马上看| 日韩熟女老妇一区二区性免费视频| 亚洲色图综合在线观看| 亚洲av电影在线进入| 精品少妇内射三级| 中文字幕另类日韩欧美亚洲嫩草| 在线 av 中文字幕| 亚洲精品自拍成人| 妹子高潮喷水视频| 美女主播在线视频| 每晚都被弄得嗷嗷叫到高潮| 国产麻豆69| 国产深夜福利视频在线观看| 一本—道久久a久久精品蜜桃钙片| 人人妻,人人澡人人爽秒播 | www.自偷自拍.com| 成人国语在线视频| 欧美97在线视频| 一边亲一边摸免费视频| 国产欧美日韩一区二区三 | 亚洲av片天天在线观看| 狠狠婷婷综合久久久久久88av| 超碰成人久久| 国产精品 国内视频| 一二三四在线观看免费中文在| 日韩一本色道免费dvd| 一级黄片播放器| 欧美日韩成人在线一区二区| 女人被躁到高潮嗷嗷叫费观| 色婷婷av一区二区三区视频| 美女午夜性视频免费| 久久久国产一区二区| 一区二区日韩欧美中文字幕| 9色porny在线观看| 成人免费观看视频高清| 欧美人与性动交α欧美软件| 欧美精品一区二区大全| 亚洲精品成人av观看孕妇| 嫁个100分男人电影在线观看 | 成人国语在线视频| 午夜免费观看性视频| 国产精品一区二区在线不卡| 看免费成人av毛片| 咕卡用的链子| 满18在线观看网站| 午夜两性在线视频| 久久精品久久久久久噜噜老黄| 日韩制服丝袜自拍偷拍| 日韩中文字幕视频在线看片| 赤兔流量卡办理| 久久人人爽av亚洲精品天堂| 看十八女毛片水多多多| 无限看片的www在线观看| 在线看a的网站| 五月天丁香电影| 国产精品av久久久久免费| 精品少妇内射三级| 亚洲图色成人| av在线播放精品| 亚洲国产精品999| √禁漫天堂资源中文www| 国产精品麻豆人妻色哟哟久久| 国产成人一区二区三区免费视频网站 | 中文字幕av电影在线播放| 在线精品无人区一区二区三| av天堂久久9| 欧美精品高潮呻吟av久久| 少妇人妻久久综合中文| 免费在线观看影片大全网站 | 午夜日韩欧美国产| 欧美黄色片欧美黄色片| 中文精品一卡2卡3卡4更新| 欧美 日韩 精品 国产| 天天操日日干夜夜撸| 一区二区三区乱码不卡18| 黄色一级大片看看| xxxhd国产人妻xxx| 精品少妇黑人巨大在线播放| 亚洲国产精品一区三区| 日韩 欧美 亚洲 中文字幕| 一本—道久久a久久精品蜜桃钙片| 纯流量卡能插随身wifi吗| 又大又爽又粗| 一边摸一边抽搐一进一出视频| 国产亚洲欧美在线一区二区| 成在线人永久免费视频| 久久国产精品人妻蜜桃| 国产xxxxx性猛交| 欧美日韩一级在线毛片| 妹子高潮喷水视频| 国产成人一区二区三区免费视频网站 | av有码第一页| 国产精品熟女久久久久浪| 亚洲精品自拍成人| 18禁观看日本| 中文字幕色久视频| kizo精华| 99国产精品99久久久久| 男女下面插进去视频免费观看| 午夜两性在线视频| 日韩制服骚丝袜av| 精品人妻在线不人妻| 高潮久久久久久久久久久不卡| 国产男女超爽视频在线观看| 欧美精品亚洲一区二区| 69精品国产乱码久久久| 国产野战对白在线观看| 男女国产视频网站| 国产精品二区激情视频| 午夜91福利影院| 久久久久久久久久久久大奶| 亚洲精品中文字幕在线视频| 国产成人免费观看mmmm| 热99久久久久精品小说推荐| 91精品伊人久久大香线蕉| 在线看a的网站| 丰满饥渴人妻一区二区三| 日韩制服骚丝袜av| 午夜福利一区二区在线看| 日日夜夜操网爽| 天堂中文最新版在线下载| 国产日韩欧美视频二区| 老司机午夜十八禁免费视频| 成人午夜精彩视频在线观看| 国产1区2区3区精品| 高清黄色对白视频在线免费看| 激情视频va一区二区三区| 少妇人妻久久综合中文| 国产精品国产av在线观看| 国产免费视频播放在线视频| 久久精品国产亚洲av高清一级| 日韩中文字幕欧美一区二区 | 午夜福利一区二区在线看| 亚洲色图综合在线观看| 国产福利在线免费观看视频| 国产成人欧美| 婷婷丁香在线五月| 亚洲,欧美,日韩| 一本综合久久免费| 亚洲人成电影观看| videos熟女内射| av一本久久久久| 国产免费现黄频在线看| 一级a爱视频在线免费观看| 女警被强在线播放| 人妻 亚洲 视频| 国产1区2区3区精品| 男女下面插进去视频免费观看| 欧美亚洲日本最大视频资源| 久久久久国产一级毛片高清牌| 人人妻人人爽人人添夜夜欢视频| 亚洲精品国产av成人精品| 欧美亚洲 丝袜 人妻 在线| 亚洲伊人色综图| 国产精品欧美亚洲77777| 国产亚洲欧美在线一区二区| 97人妻天天添夜夜摸| 日本一区二区免费在线视频| 日韩欧美一区视频在线观看| 免费一级毛片在线播放高清视频 | 丝袜美腿诱惑在线| 久久久精品区二区三区| 各种免费的搞黄视频| 国产精品二区激情视频| 国产精品熟女久久久久浪| 欧美成狂野欧美在线观看| 在线观看人妻少妇| 国产成人av激情在线播放| 成年美女黄网站色视频大全免费| 婷婷成人精品国产| 午夜激情久久久久久久| 美女高潮到喷水免费观看| 97人妻天天添夜夜摸| 久久久精品国产亚洲av高清涩受| 热99久久久久精品小说推荐| 91九色精品人成在线观看| 一区二区三区精品91| 精品亚洲成国产av| 午夜免费鲁丝| 永久免费av网站大全| 午夜精品国产一区二区电影| 美女扒开内裤让男人捅视频| 欧美激情极品国产一区二区三区| 水蜜桃什么品种好| 老司机午夜十八禁免费视频| 亚洲,欧美精品.| 国产精品久久久久成人av| 国产一区二区三区av在线| 国产又爽黄色视频| 黄色毛片三级朝国网站| 男女边摸边吃奶| 新久久久久国产一级毛片| 久久人人爽av亚洲精品天堂| 大话2 男鬼变身卡| 狠狠精品人妻久久久久久综合| 欧美日韩亚洲国产一区二区在线观看 | 男人操女人黄网站| 日韩电影二区| 大片免费播放器 马上看| 免费高清在线观看视频在线观看| 少妇精品久久久久久久| av在线播放精品| 电影成人av| 亚洲国产精品国产精品| 久久人人爽av亚洲精品天堂| 国产成人精品无人区| 999久久久国产精品视频| 国产精品三级大全| 欧美日韩综合久久久久久| 亚洲熟女毛片儿| 久久久久视频综合| 欧美激情极品国产一区二区三区| 一级a爱视频在线免费观看| av片东京热男人的天堂| 99国产精品免费福利视频| 亚洲一区中文字幕在线| 人人妻人人澡人人看| 91麻豆精品激情在线观看国产 | 精品国产一区二区三区四区第35| 99国产精品99久久久久| 亚洲五月婷婷丁香| 男人添女人高潮全过程视频| 大片免费播放器 马上看| 久久性视频一级片| avwww免费| 亚洲av片天天在线观看| 两性夫妻黄色片| 久久国产精品男人的天堂亚洲| 在线观看免费高清a一片| 欧美亚洲 丝袜 人妻 在线| 巨乳人妻的诱惑在线观看| 热re99久久国产66热| 丝袜人妻中文字幕| 国产又爽黄色视频| 亚洲成人手机| 成年人黄色毛片网站| 久久亚洲精品不卡| 一本综合久久免费| 免费在线观看黄色视频的| 欧美老熟妇乱子伦牲交| 亚洲精品日本国产第一区| 亚洲国产av影院在线观看| 99re6热这里在线精品视频| 国产精品一区二区精品视频观看| 久久久久久亚洲精品国产蜜桃av| 精品少妇黑人巨大在线播放| 久热爱精品视频在线9|