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

    基于LE和IIN算法的正弦信號頻率估計(jì)算法

    2017-09-18 01:01:13譚釧章李宏偉樊昌周
    探測與控制學(xué)報(bào) 2017年4期
    關(guān)鍵詞:運(yùn)算量正弦插值

    譚釧章,李宏偉,樊昌周,耿 耿

    (1.空軍工程大學(xué)信息與導(dǎo)航學(xué)院,陜西 西安 710077;2.解放軍95999部隊(duì),北京 100078)

    基于LE和IIN算法的正弦信號頻率估計(jì)算法

    譚釧章1,李宏偉1,樊昌周1,耿 耿2

    (1.空軍工程大學(xué)信息與導(dǎo)航學(xué)院,陜西西安710077;2.解放軍95999部隊(duì),北京100078)

    針對信號頻率位于兩個相鄰離散頻率點(diǎn)的中心區(qū)域時,迭代插值(Iterative Interpolation Based Algorithm, IIN)算法估計(jì)誤差較大,而當(dāng)信號頻率位于離散頻率點(diǎn)附近時,線性方程(Linear Equation, LE)算法估計(jì)誤差較大的問題,提出了一種基于LE和IIN算法的正弦信號頻率估計(jì)算法。該算法運(yùn)用頻譜搬移的思想,首先利用LE算法進(jìn)行頻率粗估計(jì),然后將原信號往離散頻率點(diǎn)附近頻移,再利用IIN算法進(jìn)行二次迭代。Monte-Carlo仿真結(jié)果表明,該算法頻率估計(jì)的均方根誤差在全頻段內(nèi)逼近克拉美-羅界(Cramer-Rao Lower Bound,CRLB),精度和穩(wěn)定性皆優(yōu)于R-IIN算法、M-LE算法和IIN算法,并且運(yùn)算量小于IIN算法(二次迭代)和M-LE算法。

    頻率估計(jì);頻譜搬移;插值;克拉美-羅界

    0 引言

    低信噪比條件下正弦信號的頻率估計(jì)一直是信號處理領(lǐng)域的經(jīng)典課題,在雷達(dá)和通信等領(lǐng)域有廣泛的應(yīng)用。國內(nèi)外許多學(xué)者對此有深入的研究。文獻(xiàn)[1]給出了加性復(fù)高斯白噪聲中正弦信號頻率的最大似然(Maximum Likelihood, ML)估計(jì),該算法的估計(jì)誤差逼近克拉美-羅界(CRLB),為最優(yōu)估計(jì),它需要在頻率域進(jìn)行一維搜索,計(jì)算量太大?;陔x散傅里葉變換(Discrete Fourier transform, DFT)的頻率估計(jì)方法[2]可以利用快速傅里葉變換(Fast Fourier Transform, FFT)實(shí)現(xiàn),因而在工程上得到了廣泛的應(yīng)用。對截?cái)嗟男盘栠M(jìn)行抽樣離散,然后作FFT變換后,會帶來頻譜泄漏與柵欄效應(yīng),通過對信號加窗可以有效減少頻譜泄漏的問題[3-6]。為減小柵欄效應(yīng)所帶來的估計(jì)誤差,文獻(xiàn)[7]提出的Rife算法利用信號頻譜幅度最大的兩根譜線幅度的比值進(jìn)行插值來估計(jì)頻率,高信噪比環(huán)境下能有效提高估計(jì)精度,在低信噪比條件下容易出現(xiàn)插值方向錯誤[8]。Quinn算法[9-10]利用FFT譜線的相位對相位噪聲有180°的抵抗容限,用幅值最大的三根譜線的復(fù)比值的實(shí)部進(jìn)行插值,有效地避免了Rife算法插值方向錯誤的問題。但在低信噪比條件下時,對某些頻率點(diǎn)的估計(jì)仍然波動明顯。針對這一問題,Serge提出線性方程(LE)[11-12]算法。該方法對離散正弦信號建立了線性方程模型,較Quinn算法有更強(qiáng)的抗噪性能,在信號頻率位于兩個離散頻率點(diǎn)中心區(qū)域時精度較高,在信號頻率點(diǎn)位于離散頻率點(diǎn)附近時估計(jì)精度會下降。文獻(xiàn)[13]在Quinn算法的基礎(chǔ)上,利用LE算法來判斷插值方向,提出了基于頻移修正的線性方程頻率估計(jì)算法,即M-LE算法,相比Quinn算法在低噪比條件下的估計(jì)精度和穩(wěn)定性有進(jìn)一步提高。但該算法受Quinn算法在信號頻率越靠近離散頻率點(diǎn)時估計(jì)誤差會變大的限制,在全頻段內(nèi)估計(jì)性能不穩(wěn)定,出現(xiàn)波動。Aboutanios提出一種迭代插值(IIN)算法[14-15],通過兩次以上的迭代使得MSE逼近CRLB。該算法在信號頻率位于量化頻率點(diǎn)附近時精度較高,否則要經(jīng)過多次迭代來減小估計(jì)誤差,而多次迭代會產(chǎn)生巨大的運(yùn)算量。文獻(xiàn)[16]結(jié)合Rife算法在信號頻率位于兩個相鄰量化頻率點(diǎn)的中心區(qū)域時精度較高和IIN算法在信號頻率位于量化頻率點(diǎn)附近時精度較高的優(yōu)點(diǎn),提出了R-IIN算法,在高信噪比條件下有良好的估計(jì)性能。但是Rife算法在低信噪比條件下易出現(xiàn)插值方向錯誤,因此R-IIN算法在低信噪比條件下估計(jì)誤差較大,在低信噪比條件下進(jìn)行頻率估計(jì)時精度下降。

    本文針對LE算法和IIN算法的估計(jì)誤差依信號頻率而變化且波動較大的問題,結(jié)合兩種算法的優(yōu)點(diǎn),提出了基于LE和IIN算法的正弦信號頻率估計(jì)算法。

    1 LE算法和IIN算法原理

    1.1 信號模型

    加性高斯白噪聲污染的復(fù)正弦信號模型可表示為:

    x(n)=Aej(2πfn+φ)+w(n),n=0,1,…,N-1

    (1)

    其中,A、f、φ分別為信號的振幅、歸一化頻率和初始相位,N為信號序列長度。w(n)為均值為零,方差為σ2的復(fù)高斯白噪聲序列。對x(n)作FFT,得到其傅里葉變換:

    (2)

    其中,W(k)為w(n)的傅里葉變換。其頻譜圖分辨率為Δf-fs/N,fs為采樣頻率。設(shè)k=k0時X(k)幅度最大,當(dāng)信號頻率f正好為Δf的整數(shù)倍時,真實(shí)頻率與最大譜線位置重合,估計(jì)值不會出現(xiàn)偏差。在實(shí)際中出現(xiàn)這種情況的概率很小,更常見的是信號頻率位于最大譜線位置與次大譜線位置之間,不失一般性,信號的真實(shí)頻率可表示為:

    f=(k0+δ)×Δf;δ∈[-0.5,0.5]

    (3)

    因此,插值類算法可以分兩步進(jìn)行。第一步,粗估計(jì):對信號序列作FFT得到最大幅度譜線對應(yīng)位置k0(整數(shù));第二步,精估計(jì):利用插值算法對偏差δ(分?jǐn)?shù))進(jìn)行精確估計(jì),最后根據(jù)式(3)得出頻率估計(jì)值。

    1.2LE算法原理

    在無噪聲污染環(huán)境下時,復(fù)正弦信號模型由式(1)可進(jìn)一步表示為[11-12]:

    s(n)=Aλm

    (4)

    其中,λ=ejω,A為信號幅度。 對s(n)作DFT得到:

    (5)

    式(5)中,WN=e-j2π/N,p=A(1-λN)。等式兩邊同時乘以分母,式(5)變形為:

    (6)

    在噪聲環(huán)境下,式(1)和式(2)可由式(4)、式(5)表示為:

    x(n)=s(n)+ω(n)

    (7)

    X(k)=S(k)+ψ(k)

    (8)

    為區(qū)分W(k)與WN,將W(k)用ψ(k)代替(如式(8)中所示)。式(7)、式(8)直接代入式(6)得到:

    (9)

    (10)

    (11)

    (12)

    (13)

    1.3IIN算法原理

    1)迭代一:

    (14)

    (15)

    2)迭代二:

    (16)

    (17)

    3)

    (18)

    IIN算法是把信號的頻譜往靠近|δ|=0的方向進(jìn)行頻譜搬移,即往離散頻率點(diǎn)附近頻移。當(dāng)信號頻率t位于離散頻率點(diǎn)附近時,X0.5與X-0.5的實(shí)部值大小十分接近,同時值較大,受噪聲影響較小,此時IIN算法的估計(jì)精度較高。當(dāng)信號頻率f位于兩個離散頻率點(diǎn)的中心區(qū)域時,X0.5與X-0.5兩者之中遠(yuǎn)離Xk0的實(shí)部值較小,受噪聲干擾較大,算法估計(jì)的精度下降,此時需要多次迭代才能逼近CRLB。

    2 基于LE和IIN算法的正弦信號頻率估計(jì)算法

    由以上分析可知:IIN算法在信號頻率位于離散頻率點(diǎn)附近時估計(jì)精度較高,位于兩離散頻率點(diǎn)中心區(qū)域時,受噪聲影響估計(jì)精度下降。而LE算法在低信噪比條件下具有良好的抗噪聲性能,在全頻段內(nèi)性能穩(wěn)定,且在兩個離散頻率點(diǎn)中心區(qū)域時具有較高的估計(jì)精度[11]。

    因此,對原始信號進(jìn)行離散化處理之后,首先用LE算法來判斷插值方向,確定插值方向的準(zhǔn)則為:

    (19)

    綜上所述,基于LE和IIN算法的正弦信號頻率估計(jì)算法實(shí)現(xiàn)步驟如下:

    1)對信號作FFT,得到最大譜線位置;

    3 實(shí)驗(yàn)仿真和計(jì)算量分析

    為了驗(yàn)證基于LE和IIN算法的正弦信號頻率估計(jì)算法的穩(wěn)定性與估計(jì)精度,利用Monte-Carlo模擬實(shí)驗(yàn)的方法對算法進(jìn)行仿真,并與M-LE算法[13]、IIN算法[14]和 R-IIN算法[16]進(jìn)行對比。仿真實(shí)驗(yàn)分為兩個部分:第一部分是各算法在信號頻率位于兩個離散頻率點(diǎn)之間不同位置時的估計(jì)誤差的對比,實(shí)驗(yàn)結(jié)果采用均方誤差(Mean Square Error, MSE)與克拉美-羅界的比值,即MSE/CRLB來度量;第二部分分析了各算法在不同信噪比條件下的估計(jì)誤差,實(shí)驗(yàn)結(jié)果為MSE。定義信噪比SNR=A2/σ2,在正弦信號的幅度和初始相位都未知的情況下,頻率估計(jì)的CRLB為[1]:

    (22)

    3.1 實(shí)驗(yàn)仿真與結(jié)果分析

    仿真條件1): 采樣頻率fs=1 MHz,采樣點(diǎn)數(shù)N=256,信號頻率f=Δf(20+δ),其中δ∈[-0.5,0.5],區(qū)間間隔為0.05。分別在SNR為10 dB、和-6 dB的條件下,對每個頻點(diǎn)處做1 000次Monte-Carlo仿真,得到其性能曲線分別如下圖 1 (a)和圖1 (b)所示。

    圖1 不同信號頻率處各算法MSE/CRLB值Fig.1 MSE/CRLBof each algorithm

    由圖1仿真結(jié)果可知,在信噪比較高(SNR=10 dB)的情況下,由于R-IIN和M-LE分別是在Rife算法和Quinn算法的基礎(chǔ)之上的改進(jìn),受到Rife算法和Quinn算法在全頻段內(nèi)不穩(wěn)定的限制,部分頻率點(diǎn)的誤差較大,基于LE和IIN算法的正弦信號頻率估計(jì)算法的性能曲線幾乎與IIN算法的性能曲線重合,算法在全頻段內(nèi)的穩(wěn)定性和整體的估計(jì)精度都明顯優(yōu)于R-IIN和M-LE算法;在信噪比較低(SNR=-6 dB)的情況下,R-IIN算法的估計(jì)誤差變得比較大,達(dá)到了2.5倍以上CRLB。M-LE算法的估計(jì)誤差波動明顯,MSE整體在1.4倍CRLB以下,IIN算法在|δ|→0時MSE逐步變大,整體在1.4倍CRLB以下。而基于LE和IIN算法的正弦信號頻率估計(jì)算法在全頻段內(nèi)性能都比較穩(wěn)定,均方誤差(MSE)在1.2倍CRLB以下,優(yōu)于R-IIN算法、M-LE算法和IIN算法。

    仿真條件2): 采樣頻率fs=1 MHz,采樣點(diǎn)數(shù)N=256,SNR∈[-8,2],區(qū)間間隔為0.5,分別在信號頻率f=Δf(20+0.05)和f=Δf(20+0.45)時,在每個信噪比條件下做1 000次Monte-Carlo仿真實(shí)驗(yàn),得出MSE性能曲線如下圖2(a)和圖2 (b)所示。

    圖2 不同信噪比條件下各算法估計(jì)性能Fig.2 The comparison of MSE in different SNR

    由圖2仿真結(jié)果可知M-LE算法在f=Δf(20+0.05)頻點(diǎn)處信噪比門限高于另外三種算法,信噪比達(dá)到-7 dB以下時,算法失效。R-IIN算法的MSE明顯高于另外三種算法。當(dāng)f=Δf(20+0.45)時,在信噪比達(dá)到-6 dB的條件下,R-IIN算法波動較大,估計(jì)精度下降。而基于LE和IIN算法的正弦信號頻率估計(jì)算法性能始終表現(xiàn)穩(wěn)定。當(dāng)|δ|→0時,其MSE大小與IIN算法相當(dāng) ,當(dāng)|δ|→0.5時,其MSE小于IIN算法,更加逼近CRLB。

    3.2 計(jì)算量分析

    作一次N點(diǎn)的FFT需要(N/2)log2N次復(fù)數(shù)乘法和Nlog2N次復(fù)數(shù)加法。此外,基于LE和IIN算法的正弦信號頻率估計(jì)算法在判斷插值方向時要增加兩次復(fù)數(shù)乘法和兩次復(fù)數(shù)加法。進(jìn)行二次迭代計(jì)算X±0.5時增加2N次復(fù)數(shù)乘法2(N-1)次復(fù)數(shù)加法。各算法運(yùn)算量比較如表1所示。

    表1 各算法運(yùn)算量分析

    由表1分析可知,基于LE和IIN算法的正弦信號頻率估計(jì)算法的運(yùn)算量小于M-LE算法和IIN算法(二次迭代),比R-IIN算法只多兩次復(fù)數(shù)乘法和兩次復(fù)數(shù)加法。因此,基于LE和IIN算法的正弦信號頻率估計(jì)算法在不顯著增加運(yùn)算量的情況下,提高了估計(jì)精度和穩(wěn)定性,整體計(jì)算量較小,便于硬件實(shí)現(xiàn)。

    4 結(jié)論

    本文分析了LE算法和IIN算法的特點(diǎn),并結(jié)合兩種算法提出了基于LE和IIN算法的正弦信號頻率估計(jì)算法。該算法既利用了LE算法在信號頻率位于兩離散頻率點(diǎn)中心區(qū)域時精度較高和抗噪性能良好的優(yōu)點(diǎn),避免了插值方向錯誤的問題,又發(fā)揮了IIN算法在信號頻率位于離散頻率點(diǎn)附近時精度較高的優(yōu)點(diǎn)。實(shí)驗(yàn)仿真結(jié)果表明,基于LE和IIN算法的正弦信號頻率估計(jì)算法在全頻段內(nèi)性能的穩(wěn)定性和精度整體優(yōu)于R-IIN算法、M-LE算法和IIN算法,在保證運(yùn)算量較小的同時有效地提高了估計(jì)精度和穩(wěn)定性,適用于信號頻率的實(shí)時估計(jì)且精度要求較高的場合。

    [1]Rife D C, Boorstyn R R. Single-tone parameter estimation from discrete-time observations[J]. IEEE Transactions on Information Theory, 1974, 20(5):591-598.

    [2]Palmer L C. Coarse frequency estimation using the discrete fourier transform[J]. IEEE Transactions on Information Theory, 1974, IT-20(1):104-109.

    [3]Duda K. DFT interpolation algorithm for Kaiser-Bessel and Dolph-Chebyshev windows[J]. IEEE Transactions on Instrumentation and Measurement, 2011, 60(3): 784-790.

    [4]Belega D, Petri D. Frequency estimation by two-or three-point interpolated Fourier algorithms based on cosine windows[J]. Signal Processing, 2015, 117: 115-125.

    [5]Candan ?. Fine resolution frequency estimation from three DFT samples: Case of windowed data[J]. Signal Processing, 2015, 114: 245-250.

    [6]Werner K J, Germain F G. Sinusoidal parameter estimation using quadratic interpolation around power-scaled magnitude spectrum peaks[J]. Applied Sciences, 2016, 6(10): 306.

    [7]Rife D C, Vincent G A. Use of the discrete fourier transform in the measurement of frequencies and levels of tones[J]. Bell Labs Technical Journal, 1970, 49(2):197-228.

    [8]齊國清. 幾種基于FFT的頻率估計(jì)方法精度分析[J]. 振動工程學(xué)報(bào), 2006, 19(1):86-92.

    [9]Quinn B G. Estimating frequency by interpolation using Fourier coefficients[J]. Signal Processing IEEE Transactions on, 1994, 42(5):1264-1268.

    [10]Quinn B G. Estimation of frequency, amplitude, and phase from the DFT of a time series[J]. IEEE Transactions on Signal Processing, 1997, 45(3):814-817.

    [11]Provencher S. Estimation of complex single-tone parameters in the DFT domain[J]. IEEE Transactions on Signal Processing, 2010, 58(7):3879-3883.

    [12]Provencher S. Parameters estimation of complex multitone signal in the DFT domain[J]. IEEE Transactions on Signal Processing, 2011, 59(7): 3001-3012.

    [13]李夏, 郭英, 張坤峰,等. 基于頻移修正的線性方程頻率估計(jì)算法[J]. 探測與控制學(xué)報(bào), 2015,37(4):66-71.

    [14]Aboutanios E, Mulgrew B. Iterative frequency estimation by interpolation on Fourier coefficients[J]. IEEE Transactions on Signal Processing, 2005, 53(4):1237-1242.

    [15]Aboutanios E. Estimating the parameters of sinusoids and decaying sinusoids in noise[J]. IEEE Instrumentation & Measurement Magazine, 2011, 2(14): 8-14.

    [16]周龍健, 羅景青, 房明星. 基于IIN算法和Rife算法的正弦波頻率估計(jì)算法[J]. 數(shù)據(jù)采集與處理, 2013, 28(6):839-842.

    FrequencyEstimatorofSinusoidSignalBasedonLEandIINAlgorithm

    TAN Chuanzhang1, LI Hongwei1,F(xiàn)AN Changzhou1,GENG Geng2

    (1. Information and Navigation College, Air Force Engineering University, Xi’an 710077, China; 2. Unit 95999 of PLA, Beijing 100078, China)

    Aiming at the problem that when the frequency locates in the central region of two adjacent discrete frequency points, the variance of iterative interpolation (IIN) based estimator becomes larger and the estimation precision of linear equation (LE) algorithm decrease when the frequency is near to the discrete frequency points, a new algorithm which based on LE and IIN algorithm was proposed in this paper. Firstly, a coarse frequency was estimated by the LE algorithm. Then move the original signal to the discrete frequency point via frequency shift modification. Finally, the second iteration of IIN algorithm was used to get a fine signal frequency estimate. Monte-Carlo simulation showed that the root mean square error (RMSE) approached to Cramer-Rao lower bound(CRLB) throughout whole frequency range. The performance of the proposed algorithm was better than R-IIN algorithm, M-LE algorithm and IIN algorithm both in accuracy and stability, and the computation was smaller than IIN algorithm(two iterations) and M-LE algorithm.

    frequency estimation; frequency shift modification; interpolation; Cramer-Rao lower bound

    2017-01-12

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

    :譚釧章(1993—),男,湖南湘潭人,碩士研究生,研究方向:雷達(dá)信號處理。E-mail:jzscion@163.com。

    TN911.6

    :A

    :1008-1194(2017)04-0119-05

    猜你喜歡
    運(yùn)算量正弦插值
    例說正弦定理的七大應(yīng)用
    正弦、余弦定理的應(yīng)用
    用平面幾何知識解平面解析幾何題
    基于Sinc插值與相關(guān)譜的縱橫波速度比掃描方法
    “美”在二倍角正弦公式中的應(yīng)用
    減少運(yùn)算量的途徑
    讓拋物線動起來吧,為運(yùn)算量“瘦身”
    一種改進(jìn)FFT多譜線插值諧波分析方法
    基于四項(xiàng)最低旁瓣Nuttall窗的插值FFT諧波分析
    基于VSG的正弦鎖定技術(shù)研究
    午夜免费成人在线视频| 国产成人精品无人区| 久9热在线精品视频| 色综合欧美亚洲国产小说| 精品一区二区三区av网在线观看| 丰满饥渴人妻一区二区三| 国产xxxxx性猛交| 国产三级黄色录像| 国产极品粉嫩免费观看在线| 黄色a级毛片大全视频| 亚洲av五月六月丁香网| 国产av在哪里看| 国产精品九九99| 大香蕉久久成人网| 丰满饥渴人妻一区二区三| 黑人巨大精品欧美一区二区蜜桃| 日本黄色日本黄色录像| 中文欧美无线码| 国产精品久久电影中文字幕| 无限看片的www在线观看| 黄网站色视频无遮挡免费观看| 久久狼人影院| 一夜夜www| 99久久精品国产亚洲精品| 亚洲 欧美一区二区三区| 丰满饥渴人妻一区二区三| 久久精品aⅴ一区二区三区四区| 一夜夜www| 女人精品久久久久毛片| 一级毛片精品| 99久久人妻综合| 久久久水蜜桃国产精品网| 国产精品久久视频播放| 久久午夜综合久久蜜桃| 999精品在线视频| 亚洲,欧美精品.| 啪啪无遮挡十八禁网站| 级片在线观看| 久久国产精品影院| 人人妻人人爽人人添夜夜欢视频| av福利片在线| 亚洲国产欧美网| 亚洲全国av大片| 啦啦啦 在线观看视频| 青草久久国产| 亚洲色图综合在线观看| 久久欧美精品欧美久久欧美| 香蕉丝袜av| 桃色一区二区三区在线观看| 好看av亚洲va欧美ⅴa在| 日韩欧美国产一区二区入口| 五月开心婷婷网| 制服诱惑二区| 精品一区二区三区av网在线观看| 欧美日韩亚洲高清精品| 淫秽高清视频在线观看| 久久这里只有精品19| 久久 成人 亚洲| ponron亚洲| 啦啦啦 在线观看视频| 波多野结衣一区麻豆| 我的亚洲天堂| 久久人人97超碰香蕉20202| 久久久水蜜桃国产精品网| 久久久国产成人免费| 人人妻人人添人人爽欧美一区卜| 欧美日韩亚洲高清精品| 女性生殖器流出的白浆| 桃色一区二区三区在线观看| 中文字幕精品免费在线观看视频| 一区福利在线观看| 大香蕉久久成人网| 热re99久久国产66热| 亚洲专区中文字幕在线| 成人免费观看视频高清| 欧美中文日本在线观看视频| 国产亚洲精品一区二区www| 亚洲成人免费av在线播放| 琪琪午夜伦伦电影理论片6080| 中文字幕人妻熟女乱码| 午夜福利影视在线免费观看| 精品一品国产午夜福利视频| 高清毛片免费观看视频网站 | 悠悠久久av| 国产成年人精品一区二区 | 日韩精品青青久久久久久| 国产午夜精品久久久久久| 欧美日韩视频精品一区| 一级毛片女人18水好多| 欧美一区二区精品小视频在线| 国产黄a三级三级三级人| 亚洲av熟女| 亚洲精品国产一区二区精华液| 成年版毛片免费区| 欧美午夜高清在线| 别揉我奶头~嗯~啊~动态视频| 亚洲 欧美一区二区三区| 天天躁夜夜躁狠狠躁躁| 亚洲精品中文字幕一二三四区| 狂野欧美激情性xxxx| 岛国在线观看网站| 亚洲男人天堂网一区| 在线观看免费日韩欧美大片| 精品高清国产在线一区| 亚洲第一av免费看| 国产精华一区二区三区| 丰满迷人的少妇在线观看| 亚洲精品国产区一区二| 一个人免费在线观看的高清视频| 亚洲免费av在线视频| 国产亚洲精品一区二区www| 在线观看一区二区三区激情| 国产精品久久久久久人妻精品电影| 精品久久久久久成人av| 国产无遮挡羞羞视频在线观看| www.999成人在线观看| 国产欧美日韩一区二区精品| 91九色精品人成在线观看| 日韩精品青青久久久久久| 99精国产麻豆久久婷婷| 午夜久久久在线观看| 变态另类成人亚洲欧美熟女 | 国产成人影院久久av| 久久人妻av系列| 999精品在线视频| 国产精品国产av在线观看| 国产一区二区三区在线臀色熟女 | 国产深夜福利视频在线观看| 欧美激情极品国产一区二区三区| 成人三级做爰电影| 老司机靠b影院| 成人亚洲精品一区在线观看| 欧美丝袜亚洲另类 | 在线视频色国产色| 国产精品 国内视频| 中文字幕最新亚洲高清| 波多野结衣高清无吗| 丝袜美腿诱惑在线| 正在播放国产对白刺激| 99久久人妻综合| 日日干狠狠操夜夜爽| √禁漫天堂资源中文www| 女同久久另类99精品国产91| 一个人免费在线观看的高清视频| 亚洲avbb在线观看| 免费不卡黄色视频| 国产精品98久久久久久宅男小说| 国产精品日韩av在线免费观看 | 亚洲午夜理论影院| 国产xxxxx性猛交| 在线观看午夜福利视频| 国产熟女xx| 99re在线观看精品视频| 99riav亚洲国产免费| 欧美最黄视频在线播放免费 | 欧美人与性动交α欧美精品济南到| 看免费av毛片| 国产又爽黄色视频| 欧美日韩视频精品一区| 国产精品国产高清国产av| 亚洲,欧美精品.| 国产成人av教育| 国产精品乱码一区二三区的特点 | 久久精品成人免费网站| 黑人操中国人逼视频| 国产真人三级小视频在线观看| 久久精品91无色码中文字幕| 91精品国产国语对白视频| 自线自在国产av| 亚洲av熟女| 国产精品一区二区免费欧美| 亚洲欧美日韩无卡精品| 成年人黄色毛片网站| 亚洲中文日韩欧美视频| 黄色女人牲交| 日本vs欧美在线观看视频| 色老头精品视频在线观看| 岛国视频午夜一区免费看| 国产精品秋霞免费鲁丝片| 国产亚洲精品综合一区在线观看 | 人人妻人人添人人爽欧美一区卜| 免费少妇av软件| 中国美女看黄片| 可以在线观看毛片的网站| 国产亚洲精品久久久久久毛片| 搡老岳熟女国产| 亚洲av熟女| 亚洲欧美激情综合另类| 亚洲人成网站在线播放欧美日韩| 精品免费久久久久久久清纯| 久久国产精品人妻蜜桃| 国产一区在线观看成人免费| 中文字幕人妻熟女乱码| 精品少妇一区二区三区视频日本电影| 免费观看人在逋| 级片在线观看| 成人三级黄色视频| 日韩有码中文字幕| 亚洲国产精品合色在线| videosex国产| 成人18禁高潮啪啪吃奶动态图| 久久精品91蜜桃| 男女床上黄色一级片免费看| 精品久久久久久久毛片微露脸| 精品久久久久久,| 国产三级黄色录像| 男人操女人黄网站| 久久国产精品人妻蜜桃| 久久久久久免费高清国产稀缺| 两个人免费观看高清视频| 国产av精品麻豆| 亚洲成a人片在线一区二区| 成年版毛片免费区| 国产有黄有色有爽视频| 脱女人内裤的视频| 精品一区二区三卡| 美女 人体艺术 gogo| 岛国在线观看网站| 国产一区二区三区在线臀色熟女 | 亚洲精品在线观看二区| 一级,二级,三级黄色视频| 窝窝影院91人妻| 亚洲国产精品sss在线观看 | 亚洲欧美激情在线| 黄频高清免费视频| 麻豆一二三区av精品| 精品一品国产午夜福利视频| 国产97色在线日韩免费| 亚洲成人久久性| 亚洲一码二码三码区别大吗| 99久久99久久久精品蜜桃| 亚洲激情在线av| 亚洲人成伊人成综合网2020| 色哟哟哟哟哟哟| 亚洲七黄色美女视频| 黄色怎么调成土黄色| 午夜福利欧美成人| 日本wwww免费看| 精品福利观看| 国产亚洲精品久久久久5区| 人妻丰满熟妇av一区二区三区| 国产亚洲精品第一综合不卡| 超色免费av| 欧美人与性动交α欧美软件| 免费高清在线观看日韩| 可以在线观看毛片的网站| 两个人免费观看高清视频| 99久久精品国产亚洲精品| 亚洲精品在线美女| 日韩大码丰满熟妇| 一边摸一边抽搐一进一出视频| 色综合站精品国产| 国产精品电影一区二区三区| 欧美黄色片欧美黄色片| 中文字幕色久视频| 91成人精品电影| 免费看十八禁软件| 亚洲色图av天堂| av在线天堂中文字幕 | 他把我摸到了高潮在线观看| www国产在线视频色| 日本免费一区二区三区高清不卡 | 天堂影院成人在线观看| 国产真人三级小视频在线观看| 三级毛片av免费| 久久精品91蜜桃| 国产成人精品无人区| 亚洲国产欧美网| 88av欧美| 在线播放国产精品三级| 一区二区日韩欧美中文字幕| 欧美日韩亚洲国产一区二区在线观看| 亚洲成人免费av在线播放| 欧美日韩一级在线毛片| 在线观看www视频免费| 欧美日韩乱码在线| 久久国产亚洲av麻豆专区| 午夜福利一区二区在线看| 首页视频小说图片口味搜索| 新久久久久国产一级毛片| 欧美日韩视频精品一区| 久久久久久免费高清国产稀缺| 国产精品久久久久久人妻精品电影| 91在线观看av| 久久久国产精品麻豆| xxx96com| 久久中文看片网| av网站在线播放免费| 男女午夜视频在线观看| 日本vs欧美在线观看视频| 老汉色∧v一级毛片| 中亚洲国语对白在线视频| 超碰成人久久| 国产伦人伦偷精品视频| 精品久久久精品久久久| 一个人免费在线观看的高清视频| 一区二区三区国产精品乱码| 女性被躁到高潮视频| 成年人免费黄色播放视频| 最新美女视频免费是黄的| 在线观看免费视频网站a站| 国产区一区二久久| 男人舔女人的私密视频| 亚洲人成电影观看| 亚洲欧美一区二区三区久久| 精品人妻1区二区| 日韩av在线大香蕉| 91老司机精品| 99re在线观看精品视频| www国产在线视频色| 18禁观看日本| 免费看十八禁软件| 黄色成人免费大全| 丝袜美腿诱惑在线| 亚洲一区二区三区不卡视频| 午夜视频精品福利| 成在线人永久免费视频| 午夜老司机福利片| 色综合站精品国产| av超薄肉色丝袜交足视频| 中出人妻视频一区二区| 真人做人爱边吃奶动态| 亚洲专区中文字幕在线| 久久精品人人爽人人爽视色| av超薄肉色丝袜交足视频| 正在播放国产对白刺激| 日韩欧美国产一区二区入口| 欧美成狂野欧美在线观看| 亚洲欧美日韩另类电影网站| 亚洲国产精品一区二区三区在线| 欧美日韩瑟瑟在线播放| 精品久久久久久久毛片微露脸| 国产成人免费无遮挡视频| 国产xxxxx性猛交| 亚洲一区高清亚洲精品| 日韩免费av在线播放| 成年人黄色毛片网站| 亚洲专区字幕在线| 精品免费久久久久久久清纯| 可以在线观看毛片的网站| 99国产精品免费福利视频| 成熟少妇高潮喷水视频| 在线观看一区二区三区激情| 十分钟在线观看高清视频www| 久久精品91蜜桃| 欧美日韩乱码在线| 欧美激情久久久久久爽电影 | 俄罗斯特黄特色一大片| 精品一区二区三区av网在线观看| 日本黄色日本黄色录像| 欧美日韩福利视频一区二区| 国产精品亚洲一级av第二区| 国产一区二区三区在线臀色熟女 | 精品国产国语对白av| 亚洲精品在线观看二区| 午夜亚洲福利在线播放| 成年女人毛片免费观看观看9| 色老头精品视频在线观看| 国产免费av片在线观看野外av| 99riav亚洲国产免费| 男女高潮啪啪啪动态图| 一级a爱视频在线免费观看| 亚洲成人免费av在线播放| 久久中文字幕人妻熟女| 精品无人区乱码1区二区| 色综合婷婷激情| 精品国产一区二区三区四区第35| 亚洲国产精品一区二区三区在线| 成熟少妇高潮喷水视频| 妹子高潮喷水视频| 色综合站精品国产| 国产单亲对白刺激| 亚洲精品在线美女| 国产午夜精品久久久久久| 女人被狂操c到高潮| 亚洲精品国产区一区二| 国产黄a三级三级三级人| 国产精华一区二区三区| 精品第一国产精品| 999久久久国产精品视频| 久久人妻福利社区极品人妻图片| 成熟少妇高潮喷水视频| 一级黄色大片毛片| 露出奶头的视频| 中文字幕人妻熟女乱码| 在线国产一区二区在线| 纯流量卡能插随身wifi吗| 99在线人妻在线中文字幕| 亚洲精品国产一区二区精华液| 午夜a级毛片| 久久久精品国产亚洲av高清涩受| 日韩免费av在线播放| 在线观看免费日韩欧美大片| 亚洲精品国产区一区二| 青草久久国产| 久久人人97超碰香蕉20202| 中国美女看黄片| 成人特级黄色片久久久久久久| 男女下面插进去视频免费观看| 欧美乱码精品一区二区三区| 国产成人精品在线电影| 一级毛片高清免费大全| 91成人精品电影| 长腿黑丝高跟| 成人黄色视频免费在线看| 丰满人妻熟妇乱又伦精品不卡| 亚洲国产欧美网| 99re在线观看精品视频| 精品午夜福利视频在线观看一区| 亚洲片人在线观看| 美女扒开内裤让男人捅视频| 国产精品久久久人人做人人爽| 久久草成人影院| 精品国产一区二区三区四区第35| 中文字幕色久视频| 婷婷六月久久综合丁香| 法律面前人人平等表现在哪些方面| 亚洲片人在线观看| √禁漫天堂资源中文www| 久久精品国产清高在天天线| 国产高清videossex| 黄片小视频在线播放| www.精华液| 另类亚洲欧美激情| 多毛熟女@视频| 在线免费观看的www视频| 在线av久久热| 男女下面进入的视频免费午夜 | 久久久国产欧美日韩av| 成人亚洲精品av一区二区 | www.精华液| 亚洲欧美日韩无卡精品| 男女高潮啪啪啪动态图| 男女下面进入的视频免费午夜 | 99精品在免费线老司机午夜| 国产精品久久电影中文字幕| 变态另类成人亚洲欧美熟女 | 国产精品免费视频内射| 国产精品二区激情视频| 两人在一起打扑克的视频| 日本五十路高清| 91国产中文字幕| 热re99久久精品国产66热6| 最近最新中文字幕大全免费视频| 精品免费久久久久久久清纯| 国产单亲对白刺激| 欧美成人午夜精品| 国产免费现黄频在线看| 人人澡人人妻人| 午夜a级毛片| 99久久国产精品久久久| 男人舔女人的私密视频| 大型av网站在线播放| 免费日韩欧美在线观看| 久久精品影院6| 精品欧美一区二区三区在线| 身体一侧抽搐| 啪啪无遮挡十八禁网站| 欧美日韩瑟瑟在线播放| www.999成人在线观看| 国产精品美女特级片免费视频播放器 | 最好的美女福利视频网| 俄罗斯特黄特色一大片| 在线播放国产精品三级| 亚洲熟女毛片儿| 色综合站精品国产| 老司机福利观看| 色哟哟哟哟哟哟| 日日干狠狠操夜夜爽| 丝袜美足系列| 国产高清激情床上av| 欧美日韩一级在线毛片| 99国产精品一区二区三区| 国产精品二区激情视频| 少妇的丰满在线观看| 妹子高潮喷水视频| 久久人妻熟女aⅴ| 最新美女视频免费是黄的| 成人三级黄色视频| 两性夫妻黄色片| 久久精品亚洲熟妇少妇任你| 亚洲人成77777在线视频| 男女高潮啪啪啪动态图| 亚洲av成人不卡在线观看播放网| 久久精品91无色码中文字幕| 在线十欧美十亚洲十日本专区| 久久精品影院6| 成人亚洲精品av一区二区 | 88av欧美| 成人影院久久| 久久国产乱子伦精品免费另类| 夜夜夜夜夜久久久久| 久久久久亚洲av毛片大全| 中文字幕色久视频| 国产真人三级小视频在线观看| 欧美中文日本在线观看视频| 精品人妻1区二区| 在线观看免费日韩欧美大片| 久久精品亚洲精品国产色婷小说| 日韩视频一区二区在线观看| 久久人妻av系列| 国产高清视频在线播放一区| 69av精品久久久久久| 欧美日韩精品网址| 麻豆av在线久日| 亚洲一区二区三区色噜噜 | 日日摸夜夜添夜夜添小说| 男女之事视频高清在线观看| 不卡av一区二区三区| 狠狠狠狠99中文字幕| 亚洲中文字幕日韩| 欧美精品一区二区免费开放| 亚洲美女黄片视频| 欧美另类亚洲清纯唯美| 久久狼人影院| 亚洲色图综合在线观看| 久久香蕉国产精品| 亚洲av五月六月丁香网| 免费看十八禁软件| 亚洲伊人色综图| 免费在线观看亚洲国产| 精品久久久久久,| 国产精品久久电影中文字幕| 久热这里只有精品99| 一级毛片高清免费大全| 91av网站免费观看| 动漫黄色视频在线观看| 中文字幕人妻熟女乱码| 久久亚洲精品不卡| 亚洲情色 制服丝袜| 91成年电影在线观看| 亚洲自偷自拍图片 自拍| 在线观看午夜福利视频| 不卡一级毛片| 久久天躁狠狠躁夜夜2o2o| xxxhd国产人妻xxx| 久久天堂一区二区三区四区| 国产成人欧美在线观看| 三上悠亚av全集在线观看| 最近最新免费中文字幕在线| 国产精品99久久99久久久不卡| 国产精品自产拍在线观看55亚洲| 一级毛片女人18水好多| ponron亚洲| 交换朋友夫妻互换小说| 国产高清videossex| 免费在线观看视频国产中文字幕亚洲| 欧美午夜高清在线| 后天国语完整版免费观看| 手机成人av网站| 搡老岳熟女国产| www国产在线视频色| 三上悠亚av全集在线观看| 狠狠狠狠99中文字幕| 黄色毛片三级朝国网站| 亚洲熟妇熟女久久| 亚洲在线自拍视频| 日韩欧美一区视频在线观看| 麻豆国产av国片精品| 黄色a级毛片大全视频| 久久天堂一区二区三区四区| 久久久久久亚洲精品国产蜜桃av| 亚洲精品一卡2卡三卡4卡5卡| 五月开心婷婷网| 国产1区2区3区精品| а√天堂www在线а√下载| 国产成人免费无遮挡视频| 免费观看人在逋| 不卡一级毛片| 高清毛片免费观看视频网站 | 成人18禁在线播放| 国产av在哪里看| 国产日韩一区二区三区精品不卡| 国产主播在线观看一区二区| 最近最新中文字幕大全电影3 | 精品欧美一区二区三区在线| 在线永久观看黄色视频| 亚洲精品av麻豆狂野| 国产成人精品在线电影| 日韩中文字幕欧美一区二区| 久久久国产成人免费| 少妇粗大呻吟视频| 精品久久久久久成人av| 亚洲熟妇中文字幕五十中出 | 久久中文字幕一级| 亚洲精品中文字幕在线视频| 97超级碰碰碰精品色视频在线观看| 人人妻人人爽人人添夜夜欢视频| 99久久99久久久精品蜜桃| 国产色视频综合| 人人妻,人人澡人人爽秒播| 少妇的丰满在线观看| 久久精品国产亚洲av香蕉五月| 9色porny在线观看| 日本欧美视频一区| 99久久精品国产亚洲精品| 欧美在线黄色| 在线国产一区二区在线| 亚洲精品中文字幕一二三四区| 91字幕亚洲| 免费不卡黄色视频| 日本黄色日本黄色录像| 午夜福利影视在线免费观看| 99香蕉大伊视频| 在线观看一区二区三区激情| 麻豆一二三区av精品| 久久99一区二区三区| www.自偷自拍.com| 久久久久国产一级毛片高清牌| 亚洲国产精品999在线| 日本黄色视频三级网站网址| 99精品久久久久人妻精品| 亚洲视频免费观看视频| 国产精品 欧美亚洲| 欧美激情 高清一区二区三区|