• <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ù)研究
    国产一区二区三区av在线 | a级毛片免费高清观看在线播放| 丰满乱子伦码专区| 亚洲精品久久国产高清桃花| 一级黄片播放器| 精品99又大又爽又粗少妇毛片| 成人美女网站在线观看视频| 在线a可以看的网站| 久久99精品国语久久久| АⅤ资源中文在线天堂| 春色校园在线视频观看| 99久久人妻综合| 特级一级黄色大片| 国产乱人视频| 亚洲国产欧美人成| 午夜精品在线福利| 精品国内亚洲2022精品成人| 两个人的视频大全免费| 亚洲av成人av| 国产精品久久久久久精品电影小说 | 99在线人妻在线中文字幕| 18禁在线无遮挡免费观看视频| 久久精品国产亚洲网站| 看片在线看免费视频| 亚洲电影在线观看av| 热99在线观看视频| 美女内射精品一级片tv| 亚洲不卡免费看| 午夜视频国产福利| 亚洲av中文av极速乱| 国产亚洲91精品色在线| 99riav亚洲国产免费| 如何舔出高潮| 老熟妇乱子伦视频在线观看| avwww免费| 淫秽高清视频在线观看| 51国产日韩欧美| 亚洲欧美日韩高清在线视频| 黄色日韩在线| 久久这里只有精品中国| 亚洲国产精品成人综合色| 桃色一区二区三区在线观看| 国产伦理片在线播放av一区 | 18禁裸乳无遮挡免费网站照片| 91精品国产九色| 亚洲成av人片在线播放无| 悠悠久久av| 免费av毛片视频| 色播亚洲综合网| 国产精品伦人一区二区| 搞女人的毛片| 3wmmmm亚洲av在线观看| 亚洲av第一区精品v没综合| 99热6这里只有精品| 联通29元200g的流量卡| 国产精品综合久久久久久久免费| 五月玫瑰六月丁香| 午夜爱爱视频在线播放| 亚洲精华国产精华液的使用体验 | 中出人妻视频一区二区| 中国美白少妇内射xxxbb| 久久久久久九九精品二区国产| 久久中文看片网| 亚洲精品影视一区二区三区av| 色视频www国产| 看黄色毛片网站| 日韩一本色道免费dvd| 天堂√8在线中文| 国产亚洲av嫩草精品影院| 午夜福利成人在线免费观看| 黄色日韩在线| 国产一区亚洲一区在线观看| 能在线免费观看的黄片| 日韩在线高清观看一区二区三区| 大型黄色视频在线免费观看| 久久久久久国产a免费观看| 国产精品嫩草影院av在线观看| 观看免费一级毛片| 人人妻人人澡人人爽人人夜夜 | 女同久久另类99精品国产91| 日韩三级伦理在线观看| 欧美xxxx性猛交bbbb| 看片在线看免费视频| 91午夜精品亚洲一区二区三区| 久久人人爽人人爽人人片va| 最近2019中文字幕mv第一页| 久久精品国产鲁丝片午夜精品| 国产色婷婷99| 亚洲成人久久爱视频| 在现免费观看毛片| 午夜a级毛片| 日日撸夜夜添| 欧美性猛交黑人性爽| 青青草视频在线视频观看| 亚洲精品影视一区二区三区av| 亚洲国产精品成人久久小说 | 日本免费一区二区三区高清不卡| 国内精品美女久久久久久| 免费看a级黄色片| 26uuu在线亚洲综合色| 国产黄a三级三级三级人| 综合色丁香网| 亚洲婷婷狠狠爱综合网| 亚洲成人久久爱视频| 别揉我奶头 嗯啊视频| 看黄色毛片网站| 国产亚洲精品久久久久久毛片| 免费看a级黄色片| 久久精品夜夜夜夜夜久久蜜豆| 九色成人免费人妻av| 一级毛片久久久久久久久女| 免费人成视频x8x8入口观看| 国模一区二区三区四区视频| 亚洲国产精品成人久久小说 | 国产黄片美女视频| 欧美日韩综合久久久久久| 国产日韩欧美在线精品| 蜜桃久久精品国产亚洲av| 美女高潮的动态| 久99久视频精品免费| 最后的刺客免费高清国语| 69人妻影院| 好男人视频免费观看在线| 日韩av在线大香蕉| 国产视频首页在线观看| 国产国拍精品亚洲av在线观看| 亚洲av第一区精品v没综合| 可以在线观看的亚洲视频| 午夜视频国产福利| 欧美日韩精品成人综合77777| 一个人观看的视频www高清免费观看| 国产精品蜜桃在线观看 | 97在线视频观看| 午夜激情欧美在线| 熟妇人妻久久中文字幕3abv| 亚洲欧洲国产日韩| 国产亚洲av嫩草精品影院| 亚洲三级黄色毛片| 九九在线视频观看精品| 日韩欧美三级三区| 春色校园在线视频观看| 国产精品久久电影中文字幕| 91久久精品电影网| 一级二级三级毛片免费看| 成年免费大片在线观看| 在线天堂最新版资源| 爱豆传媒免费全集在线观看| 夜夜夜夜夜久久久久| 免费看光身美女| 久久精品夜夜夜夜夜久久蜜豆| 性色avwww在线观看| 国内精品久久久久精免费| 边亲边吃奶的免费视频| 欧美日韩乱码在线| 99久久成人亚洲精品观看| 日本黄色视频三级网站网址| a级毛片免费高清观看在线播放| 又爽又黄无遮挡网站| 国产伦一二天堂av在线观看| 久久久久久伊人网av| 国内精品宾馆在线| 亚洲精华国产精华液的使用体验 | 丝袜喷水一区| 久久精品人妻少妇| 午夜亚洲福利在线播放| 国产色爽女视频免费观看| 22中文网久久字幕| 亚洲av熟女| 一卡2卡三卡四卡精品乱码亚洲| 国产亚洲精品久久久久久毛片| 变态另类丝袜制服| 99久国产av精品| 国产白丝娇喘喷水9色精品| 亚洲综合色惰| 一进一出抽搐动态| 日韩成人伦理影院| 丰满的人妻完整版| 午夜精品在线福利| 亚洲最大成人中文| 亚洲熟妇中文字幕五十中出| 亚洲一区高清亚洲精品| 欧美区成人在线视频| 亚洲自偷自拍三级| 久久草成人影院| 在线观看66精品国产| 亚洲自偷自拍三级| 婷婷色av中文字幕| 国产成人精品久久久久久| 淫秽高清视频在线观看| 国产视频内射| av.在线天堂| 久久久精品94久久精品| 亚洲精品乱码久久久v下载方式| 在线观看免费视频日本深夜| 国产黄a三级三级三级人| 亚洲一区高清亚洲精品| 两个人的视频大全免费| 精品久久久噜噜| 成人毛片60女人毛片免费| 久久综合国产亚洲精品| 国产 一区 欧美 日韩| 日韩人妻高清精品专区| 国产精品久久久久久精品电影| 亚洲熟妇中文字幕五十中出| 日韩一区二区视频免费看| 国产免费男女视频| 在线a可以看的网站| 黄色日韩在线| av免费观看日本| 两个人视频免费观看高清| 国产人妻一区二区三区在| 一本精品99久久精品77| 亚洲美女搞黄在线观看| 成年免费大片在线观看| 亚洲第一电影网av| 高清在线视频一区二区三区 | 日本免费a在线| 最后的刺客免费高清国语| 亚洲欧美中文字幕日韩二区| 成人av在线播放网站| 人妻系列 视频| 国产激情偷乱视频一区二区| 日本与韩国留学比较| 日韩av在线大香蕉| 国内精品久久久久精免费| 村上凉子中文字幕在线| 美女cb高潮喷水在线观看| 欧美变态另类bdsm刘玥| 久久这里只有精品中国| a级毛色黄片| 日韩,欧美,国产一区二区三区 | 69人妻影院| 国内精品久久久久精免费| 国产午夜精品久久久久久一区二区三区| av在线播放精品| 国产综合懂色| 亚洲自拍偷在线| 国产成人影院久久av| 久久这里有精品视频免费| 成人三级黄色视频| 能在线免费看毛片的网站| 亚洲丝袜综合中文字幕| 欧美bdsm另类| 日韩一区二区三区影片| 九九在线视频观看精品| 亚洲七黄色美女视频| 欧美一级a爱片免费观看看| av在线观看视频网站免费| 欧美在线一区亚洲| 赤兔流量卡办理| eeuss影院久久| 精品欧美国产一区二区三| 国产精品99久久久久久久久| 哪里可以看免费的av片| 又爽又黄a免费视频| 午夜福利高清视频| 成人无遮挡网站| 国产午夜精品一二区理论片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | kizo精华| 久久久成人免费电影| 午夜激情欧美在线| 国产极品精品免费视频能看的| 日本免费a在线| 久久人妻av系列| 午夜福利在线在线| 禁无遮挡网站| 婷婷精品国产亚洲av| 国产精品三级大全| av在线老鸭窝| 搡女人真爽免费视频火全软件| 国产精品久久久久久精品电影小说 | 大香蕉久久网| 日韩强制内射视频| 国产 一区 欧美 日韩| 蜜臀久久99精品久久宅男| 欧美xxxx黑人xx丫x性爽| 亚洲人成网站在线观看播放| 国产精品人妻久久久影院| 女人被狂操c到高潮| 天堂中文最新版在线下载 | 99久国产av精品| 国产免费一级a男人的天堂| 夜夜爽天天搞| 国产精品.久久久| 免费大片18禁| 欧美日本视频| 少妇的逼好多水| 老女人水多毛片| 亚洲精品成人久久久久久| 最近2019中文字幕mv第一页| 久久热精品热| 成人性生交大片免费视频hd| 永久网站在线| 国模一区二区三区四区视频| 黄色配什么色好看| 一夜夜www| 欧美变态另类bdsm刘玥| 久久久久久久午夜电影| 成人午夜精彩视频在线观看| 国国产精品蜜臀av免费| 联通29元200g的流量卡| 国产av不卡久久| 99久久成人亚洲精品观看| 18禁黄网站禁片免费观看直播| 韩国av在线不卡| 1000部很黄的大片| 国产片特级美女逼逼视频| 三级男女做爰猛烈吃奶摸视频| 午夜亚洲福利在线播放| 亚洲在久久综合| 舔av片在线| 成人一区二区视频在线观看| 国产成人精品婷婷| 国产成人a区在线观看| 婷婷六月久久综合丁香| 精品久久久久久成人av| 国产高清激情床上av| 国产单亲对白刺激| 久久99蜜桃精品久久| 国产成人aa在线观看| 精品熟女少妇av免费看| 男插女下体视频免费在线播放| 亚洲经典国产精华液单| 男女视频在线观看网站免费| 亚洲av二区三区四区| 久久久久久久午夜电影| 高清在线视频一区二区三区 | 午夜福利成人在线免费观看| 亚洲欧美精品自产自拍| 亚洲高清免费不卡视频| 中文资源天堂在线| 国产一区二区亚洲精品在线观看| 高清毛片免费看| 国产精品人妻久久久久久| 麻豆国产97在线/欧美| 国产亚洲欧美98| 亚洲人成网站在线播放欧美日韩| 色播亚洲综合网| 国内精品久久久久精免费| 亚洲精品自拍成人| 亚洲国产精品成人久久小说 | 久久久久久久久久成人| 久久久精品欧美日韩精品| 内地一区二区视频在线| 精品一区二区免费观看| 亚洲va在线va天堂va国产| 乱码一卡2卡4卡精品| 欧美精品一区二区大全| 日本一本二区三区精品| 久久久欧美国产精品| 此物有八面人人有两片| 久久婷婷人人爽人人干人人爱| 国产精华一区二区三区| 国产久久久一区二区三区| 国产伦在线观看视频一区| 亚洲国产精品sss在线观看| 男女视频在线观看网站免费| 热99re8久久精品国产| 久久久久久久久中文| 免费看a级黄色片| 在线观看美女被高潮喷水网站| 麻豆精品久久久久久蜜桃| 老熟妇乱子伦视频在线观看| 亚洲成人中文字幕在线播放| 一级毛片电影观看 | 黄色欧美视频在线观看| av国产免费在线观看| 如何舔出高潮| 国产精品久久久久久亚洲av鲁大| 91久久精品国产一区二区成人| 日韩欧美在线乱码| 国产精品免费一区二区三区在线| 如何舔出高潮| 在线观看一区二区三区| 美女 人体艺术 gogo| 中国美白少妇内射xxxbb| 青春草亚洲视频在线观看| 尤物成人国产欧美一区二区三区| av.在线天堂| 亚洲最大成人av| 在现免费观看毛片| 亚洲天堂国产精品一区在线| 婷婷精品国产亚洲av| 国产精品av视频在线免费观看| 国产精品一二三区在线看| 搞女人的毛片| 3wmmmm亚洲av在线观看| 国产中年淑女户外野战色| 久久久久久久久久成人| 久久精品人妻少妇| 精品久久久久久久久av| 午夜激情福利司机影院| 99久久成人亚洲精品观看| 九草在线视频观看| 亚洲真实伦在线观看| 亚洲无线在线观看| 在线免费观看不下载黄p国产| 欧美一级a爱片免费观看看| 高清午夜精品一区二区三区 | 日日撸夜夜添| 欧美激情在线99| 色视频www国产| 熟妇人妻久久中文字幕3abv| 波多野结衣高清作品| 精品久久久久久成人av| 丰满的人妻完整版| 久久精品国产99精品国产亚洲性色| 伦精品一区二区三区| 国产黄片视频在线免费观看| 亚洲人成网站在线观看播放| 一个人免费在线观看电影| 成人综合一区亚洲| 菩萨蛮人人尽说江南好唐韦庄 | 欧美性猛交黑人性爽| 国产成人午夜福利电影在线观看| 国产精品,欧美在线| 国产v大片淫在线免费观看| avwww免费| 午夜激情福利司机影院| 亚洲成人中文字幕在线播放| 99热6这里只有精品| 寂寞人妻少妇视频99o| 国产精品美女特级片免费视频播放器| 国产极品精品免费视频能看的| 亚洲欧美中文字幕日韩二区| 国产亚洲精品久久久久久毛片| av在线播放精品| 夫妻性生交免费视频一级片| 校园人妻丝袜中文字幕| 日韩成人伦理影院| 欧美在线一区亚洲| 日日摸夜夜添夜夜添av毛片| 亚洲av一区综合| 国产精品电影一区二区三区| 国产精品,欧美在线| 99热精品在线国产| 丰满的人妻完整版| 99在线视频只有这里精品首页| 亚洲av中文av极速乱| 好男人在线观看高清免费视频| ponron亚洲| 精品久久久久久久末码| 成年版毛片免费区| 两个人视频免费观看高清| 波多野结衣巨乳人妻| 在线免费观看的www视频| 亚洲av电影不卡..在线观看| 久久婷婷人人爽人人干人人爱| 亚洲精品456在线播放app| 国产极品精品免费视频能看的| 久久久久网色| 日本成人三级电影网站| 国产成人精品婷婷| 国产伦一二天堂av在线观看| 久久久久久大精品| 亚洲最大成人中文| 18禁在线播放成人免费| 久久久久久国产a免费观看| 国产精品永久免费网站| 99热只有精品国产| 国产亚洲av片在线观看秒播厂 | 午夜精品国产一区二区电影 | 欧美又色又爽又黄视频| 亚洲一区高清亚洲精品| 亚洲av中文av极速乱| 特大巨黑吊av在线直播| 日本色播在线视频| 真实男女啪啪啪动态图| 日韩欧美一区二区三区在线观看| 可以在线观看毛片的网站| 国产精品免费一区二区三区在线| 中文字幕av在线有码专区| 国产极品天堂在线| 欧美另类亚洲清纯唯美| 免费看a级黄色片| 欧美高清性xxxxhd video| 2022亚洲国产成人精品| 非洲黑人性xxxx精品又粗又长| 久久国产乱子免费精品| 国产单亲对白刺激| 91麻豆精品激情在线观看国产| 亚洲国产精品久久男人天堂| 国产成人精品婷婷| 久久久久久久亚洲中文字幕| 偷拍熟女少妇极品色| 国产精品野战在线观看| 成人漫画全彩无遮挡| 国产精品一区二区在线观看99 | 小说图片视频综合网站| 中国美女看黄片| 国产精品一区二区性色av| 99久久精品国产国产毛片| 2022亚洲国产成人精品| 51国产日韩欧美| 在线观看一区二区三区| 精品99又大又爽又粗少妇毛片| 一进一出抽搐gif免费好疼| 婷婷色av中文字幕| 亚洲av二区三区四区| 精品久久久久久久久av| 成人欧美大片| 亚洲aⅴ乱码一区二区在线播放| 欧美+日韩+精品| 亚洲va在线va天堂va国产| 日韩,欧美,国产一区二区三区 | 99热这里只有是精品50| 欧美一区二区亚洲| 国产亚洲精品av在线| 国产精品久久视频播放| 日日摸夜夜添夜夜添av毛片| 美女xxoo啪啪120秒动态图| 精品少妇黑人巨大在线播放 | 成年免费大片在线观看| 一进一出抽搐gif免费好疼| 国产v大片淫在线免费观看| 久久国产乱子免费精品| 欧美xxxx性猛交bbbb| 天天躁夜夜躁狠狠久久av| 国产片特级美女逼逼视频| 亚洲精品日韩av片在线观看| 成人二区视频| 日韩欧美三级三区| 精品国内亚洲2022精品成人| 亚洲美女视频黄频| 级片在线观看| 亚洲精品久久久久久婷婷小说 | 欧美+亚洲+日韩+国产| 噜噜噜噜噜久久久久久91| 精品久久久久久久人妻蜜臀av| 久久99热这里只有精品18| 51国产日韩欧美| 久久精品91蜜桃| 国产人妻一区二区三区在| 成人特级黄色片久久久久久久| 男人和女人高潮做爰伦理| 久久久欧美国产精品| АⅤ资源中文在线天堂| 亚洲av电影不卡..在线观看| 最新中文字幕久久久久| 国产一区亚洲一区在线观看| 欧美xxxx性猛交bbbb| 亚洲欧美日韩高清在线视频| 亚洲精品日韩在线中文字幕 | 人妻系列 视频| 国内精品美女久久久久久| 18禁裸乳无遮挡免费网站照片| 日本熟妇午夜| 国产精品三级大全| 精品日产1卡2卡| 哪个播放器可以免费观看大片| 麻豆成人午夜福利视频| 中文亚洲av片在线观看爽| 麻豆国产97在线/欧美| 老司机福利观看| 91久久精品电影网| 免费人成视频x8x8入口观看| 真实男女啪啪啪动态图| 亚洲国产欧美在线一区| 午夜精品国产一区二区电影 | 日韩高清综合在线| 看免费成人av毛片| 久久精品国产自在天天线| 国产片特级美女逼逼视频| 国产高清不卡午夜福利| 国产亚洲精品久久久com| 天天一区二区日本电影三级| 欧美在线一区亚洲| 国产伦在线观看视频一区| 国产黄片视频在线免费观看| 国产伦在线观看视频一区| 欧美激情在线99| 国产av一区在线观看免费| 日韩欧美一区二区三区在线观看| 又黄又爽又刺激的免费视频.| 你懂的网址亚洲精品在线观看 | av在线亚洲专区| 男女下面进入的视频免费午夜| 国产中年淑女户外野战色| 成人三级黄色视频| 毛片女人毛片| 日韩一本色道免费dvd| 亚洲国产色片| av在线观看视频网站免费| 午夜爱爱视频在线播放| 国产精品无大码| 精品人妻熟女av久视频| 2021天堂中文幕一二区在线观| 免费搜索国产男女视频| 女同久久另类99精品国产91| 69人妻影院| 天堂网av新在线| 少妇丰满av| 亚洲av不卡在线观看| 九九久久精品国产亚洲av麻豆| 一本久久精品| h日本视频在线播放| 免费观看的影片在线观看| 亚洲经典国产精华液单| 欧美性猛交黑人性爽| 日韩欧美国产在线观看| 国产精品无大码| 老司机福利观看| 日本-黄色视频高清免费观看| 淫秽高清视频在线观看| 能在线免费观看的黄片| 日本与韩国留学比较| 卡戴珊不雅视频在线播放| 亚洲成人av在线免费| 免费观看人在逋| 欧洲精品卡2卡3卡4卡5卡区| av又黄又爽大尺度在线免费看 | 黄色日韩在线| 国产精品三级大全| 精品日产1卡2卡|