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

    低信噪比下高精度通信信號(hào)載頻盲估計(jì)算法

    2017-04-20 07:56:46張亞豐趙成林
    軟件 2017年3期
    關(guān)鍵詞:高斯復(fù)雜度信噪比

    張亞豐,王 翔,趙成林

    (1.北京郵電大學(xué) 信息與通信工程學(xué)院,北京 100876;2.北京簡(jiǎn)易科技有限公司,北京 100088)

    低信噪比下高精度通信信號(hào)載頻盲估計(jì)算法

    張亞豐1,王 翔2,趙成林1

    (1.北京郵電大學(xué) 信息與通信工程學(xué)院,北京 100876;2.北京簡(jiǎn)易科技有限公司,北京 100088)

    未知制式的通信信號(hào)的載頻估計(jì)是臺(tái)站監(jiān)測(cè)與頻譜管理中的一個(gè)重要問(wèn)題,針對(duì)現(xiàn)有載頻估計(jì)算法的估計(jì)精度和復(fù)雜度不能兼得的不足,在分析信號(hào)功率譜特性的基礎(chǔ)上,結(jié)合線性調(diào)制信號(hào)的功率譜的對(duì)稱性,在頻譜重心法的基礎(chǔ)上提出了以一種改進(jìn)的估計(jì)算法,首先對(duì)信號(hào)的功率譜進(jìn)行平滑和高斯擬合的預(yù)處理,然后通過(guò)迭代搜索的方法估計(jì)載頻大小,該算法無(wú)需知道信號(hào)的制式以及其他信息即可實(shí)現(xiàn)載頻的盲估計(jì),仿真結(jié)果表明,在高斯信道環(huán)境下,該算法可以實(shí)現(xiàn)對(duì)載頻的快速估計(jì),并且其在低信噪比下的估計(jì)較為精確,在臺(tái)站監(jiān)測(cè)方面有著良好的應(yīng)用前景。

    載頻估計(jì);譜重心法;高斯擬合;迭代搜索

    0 引言

    缺乏先驗(yàn)知識(shí)的通信信號(hào)載波頻率盲估計(jì)是臺(tái)站以及頻譜監(jiān)測(cè)領(lǐng)域[1-2]的一個(gè)非常重要的部分,載波頻率估計(jì)的精確度對(duì)監(jiān)測(cè)或截獲信號(hào)后續(xù)的調(diào)制方式識(shí)別以及解調(diào)有著很大的影響。為了監(jiān)測(cè)合法臺(tái)站的工作情況以及頻點(diǎn)使用情況,并且需要適應(yīng)惡劣信道狀況,因此需要一種抗干擾比較強(qiáng)的高精度載頻估計(jì)算法。目前針對(duì)載頻估計(jì)方法已經(jīng)有了很多研究,大致可以歸為時(shí)域與頻域兩類。文獻(xiàn)[3]根據(jù)數(shù)字調(diào)制信號(hào)在每個(gè)載波周期會(huì)兩次經(jīng)過(guò)零點(diǎn)的特性,統(tǒng)計(jì)了接收到的信號(hào)在一段時(shí)間內(nèi)經(jīng)過(guò)零點(diǎn)的次數(shù)來(lái)估計(jì)載頻,該方法復(fù)雜度低,但是在低信噪比下的估計(jì)精確度迅速下降。文獻(xiàn)[4]采用了一種相位差分法,在高信噪比下對(duì)相位調(diào)制信號(hào)的估計(jì)性能良好,但是對(duì)非相位調(diào)制信號(hào)估計(jì)性能受限,且由于是時(shí)域估計(jì)方法,對(duì)噪聲比較敏感。文獻(xiàn)[5]對(duì)信號(hào)進(jìn)行了M倍頻,頻譜搬移到了M*fc處,然后搜索頻譜最大值,就可以得到載頻的估計(jì)值,該方法需要知道信號(hào)調(diào)制方式等先驗(yàn)知識(shí),而且在倍頻的時(shí)候,噪聲也相應(yīng)的被同樣的處理,噪聲影響被放大,估計(jì)結(jié)果精度降低。文獻(xiàn)[6]針對(duì)時(shí)域算法對(duì)噪聲的敏感性,提出了一種經(jīng)典的載頻估計(jì)方法,頻率重心法,通過(guò)計(jì)算頻譜的重心位置來(lái)估計(jì)信號(hào)載頻,方法簡(jiǎn)單,但是在估計(jì)性能同樣在低信噪比下受限。因?yàn)槌S靡苿?dòng)通信信號(hào)都具有循環(huán)平穩(wěn)的特性,文獻(xiàn)[7,8]提出了基于循環(huán)譜和循環(huán)累積量的載頻估計(jì)算法,這兩種方法可以抑制噪聲,尤其高階循環(huán)累積量理論上可以完全消除噪聲,相對(duì)前幾種方法估計(jì)精度有了很大提高,但是這兩種方法運(yùn)算復(fù)雜度也成倍的增加。

    本文在研究了移動(dòng)通信信號(hào)功率譜特性和文獻(xiàn)中所提到的算法的特性,針對(duì)目前估計(jì)算法低復(fù)雜度與高準(zhǔn)確性不可兼得的不足,在復(fù)雜度較低的頻譜中心法基礎(chǔ)上,提出了一種改進(jìn)的低信噪比下低復(fù)雜度載頻估計(jì)算法,該方法首先利用了直接法求取信號(hào)的功率譜,并對(duì)其平滑和進(jìn)行高斯擬合得到平滑的功率譜曲線,然后利用功率譜對(duì)稱性,對(duì)功率譜曲線左右快速上升和下降的區(qū)間進(jìn)行求和逼近,進(jìn)而求得載頻估計(jì)值,最后針對(duì)2G、3G、4G移動(dòng)通信信號(hào)的載頻估計(jì)進(jìn)行了MonteCarlo仿真實(shí)驗(yàn)。

    1 信號(hào)模型的建立和功率譜估計(jì)

    1.1 信號(hào)模型

    此處三種制式的信號(hào)可以處理為一個(gè)接收信號(hào)模型,設(shè)接收到的信號(hào)有以下的形式:

    其中,s(t)為經(jīng)過(guò)載頻調(diào)制后的移動(dòng)通信信號(hào),I (t)和Q(t)分別為基帶調(diào)制信號(hào)的同相分量和正交分量,具體表達(dá)形式如式(3)與(4)所示,n(t)為零均值的高斯白噪聲, fc為載波頻率。

    其中cn為發(fā)送符號(hào),承載了需要傳輸?shù)男畔?,T為符號(hào)周期,g(t)為脈沖成型。

    1.2 功率譜估計(jì)

    在實(shí)際的監(jiān)測(cè)以及非協(xié)作通信系統(tǒng)中,接收到的通信信號(hào)是具有某種統(tǒng)計(jì)特性的隨機(jī)信號(hào),并且由于需要對(duì)信號(hào)進(jìn)行快速的分析,因此,本文選用最簡(jiǎn)單的譜估計(jì)方法--周期圖法,對(duì)該信號(hào)進(jìn)行譜分析,功率譜Q(f)的表達(dá)式,如式(5)所示:

    其中XN(f)為采樣得到的信號(hào)x(n)的傅里葉變換,周期圖法雖然是漸進(jìn)無(wú)偏估計(jì),但是因?yàn)橄鄬?duì)其他平譜分析方法,周期圖法的分辨率較高,因此得到的功率譜不平滑,且變化幅度較大。所以,針對(duì)周期圖法的一些缺點(diǎn),采用對(duì)譜線加窗進(jìn)行中值平滑的方法來(lái)功率譜進(jìn)行第一步處理,經(jīng)過(guò)平滑后的功率譜線仍然會(huì)有很多毛刺,因此針對(duì)平滑后的功率譜線采用高斯曲線擬合對(duì)譜線進(jìn)行第二步處理。

    2 改進(jìn)的載頻盲估計(jì)算法

    2.1 高斯非線性擬合[9]

    經(jīng)過(guò)平滑的功率譜細(xì)節(jié)處還是有很多毛刺,因此在此基礎(chǔ)上對(duì)功率譜做一個(gè)簡(jiǎn)單的高斯多項(xiàng)式擬合會(huì)提升功率譜平滑效果,以便后續(xù)的處理。一般對(duì)數(shù)據(jù)進(jìn)行擬合的方法是:對(duì)給定的仿真數(shù)據(jù)(xi,yi),(i=0,1,2,3,···,m),在選取的函數(shù)類φ中,通過(guò)若干次的迭代,求得 f(x)∈φ,使 f(x)與yi的誤差平方和最小[10],即

    求得的 f(x)即為實(shí)際數(shù)據(jù)的擬合函數(shù),求取f(x)的過(guò)程為曲線擬合。在曲線擬合中,可以選擇不同的基函數(shù)進(jìn)行擬合,本文選用高斯函數(shù)進(jìn)行曲線擬合。

    高斯擬合主要使用形如式(8)的函數(shù)對(duì)數(shù)據(jù)集進(jìn)行擬合,高斯擬合和多項(xiàng)式擬合類似,也是運(yùn)用最小二乘法進(jìn)行逼近,不同的是函數(shù)不一樣:

    k為階數(shù), f(x)即為k個(gè)高斯函數(shù)的組合。

    求取高斯函數(shù)的主要原理如下:

    設(shè)現(xiàn)在有一組需要擬合的數(shù)據(jù)為(xi,yi)(i=0, 1,2,3,···,m),可用式(9)高斯函數(shù)描述,即

    式(9)中需要估計(jì)的參數(shù)有三個(gè),為a,b,c,分別代表的物理意義為高斯曲線的峰高、半寬度信息和峰位置,將(9)兩邊取自然對(duì)數(shù)得:

    則式(10)化為二次多項(xiàng)式擬合函數(shù)

    考慮全部數(shù)據(jù)和量測(cè)誤差,用矩陣形式表示如下:

    記為Z=XB,可求得矩陣B的廣義最小二乘解為B=(XTX)-1XTZ ,由此可以求得參數(shù)a,b,c,進(jìn)而求得(8),本文選用高斯基函數(shù)為4個(gè)。

    2.2 算法原理

    用WCDMA信號(hào)做事例分析,信號(hào)模型如第一節(jié)所示,碼率為3.84 Mcps,圖1是WCDMA信號(hào)未加噪聲的功率譜梯度曲線,其他兩種信號(hào)類似,從圖中可以看出,功率譜上升和下降區(qū)間比較窄,即譜線上升或者下降比較迅速,因此受噪聲影響比較小。

    圖1 功率譜梯度曲線

    經(jīng)過(guò)線性調(diào)制的信號(hào)的功率譜與信號(hào)載頻有式(12)的關(guān)系,即功率譜上升段與下降段是關(guān)于載頻對(duì)稱的,因此利用功率譜對(duì)稱性提出了在譜重心法估計(jì)載頻的基礎(chǔ)上的修正算法。

    算法估計(jì)流程的詳細(xì)描述如下:

    Step1:先對(duì)信號(hào)進(jìn)行求取功率譜,本文采取周期圖法,由公式(1)對(duì)x(n)求取FFT,可以的到對(duì)應(yīng)的頻譜,然后根據(jù)公式(14)可以得到功率譜

    Step2:對(duì)功率譜進(jìn)行平滑和高斯擬合,然后頻譜中心法公式(3),可以得到 fd:

    Step3:根據(jù)Step2中的 fd或者根據(jù)搜到的最大功率譜值Maxi,向兩邊搜索3dB帶寬的大致位置,可以得到信號(hào)的3db帶寬band=Band_3dB(Hz )

    Step4:在 fd左右兩側(cè)選取相同長(zhǎng)度大小的頻帶,即需要保證,fL1和 fL2以及 fR1和 fR2的確定,由上一步頻率重心法得到的 fd以及下面的公式求取,公式如下:Step5:根據(jù)兩頻帶區(qū)間的功率譜值之和進(jìn)行算法迭代,迭代的邏輯控制偽代碼如下:

    停止,否則繼續(xù)進(jìn)行step5,經(jīng)過(guò)若干次迭代以后,可以得到載頻的最終估計(jì)結(jié)果是

    注:

    1.在步驟 4中,c1取常數(shù)(試驗(yàn)中取 0.6),distanc是常量,在算法初始化以后,在迭代過(guò)程中保持不變。

    2.上面在步驟 5中,step取值在5~10之間,取值越小越精確,取值越大迭代速度越快。

    3.Threshold由如下公式確定

    c根據(jù)實(shí)際情況取1~3之間的常數(shù)。

    3 仿真結(jié)果結(jié)果與分析

    3.1 仿真結(jié)果

    本文分別對(duì)三種制式的移動(dòng)通信信號(hào)在不同信噪比進(jìn)行了仿真和分析,仿真結(jié)果主要包括兩部分內(nèi)容,第一部分主要論證了算法的收斂性;第二部分主要評(píng)估算法對(duì)信號(hào)載頻的估計(jì)性能,其中采用了兩中估計(jì)方案:1)沒(méi)有經(jīng)過(guò)優(yōu)化的頻譜重心法,2)本文提出的迭代算法。在不知道信號(hào)先驗(yàn)條件下,以及較低信噪比的情況下,可以明顯看出,算法2)的估計(jì)性能比算法1)的性能更好。

    3.1.1 算法的收斂性能的仿真

    在不同信噪比下,分別對(duì)三種制式的通信信號(hào)進(jìn)行100次蒙特卡洛仿真實(shí)驗(yàn),得到的算法平均迭代次數(shù)的仿真圖如圖2所示。

    圖2 不同信噪比下算法的迭代次數(shù)

    從圖2我們可以看出算法速度收斂非???,并且隨著信噪比的提升,算法的迭代次數(shù)在減少,因?yàn)樵陔S著信噪比的升高,譜重心法的估計(jì)性能也在提升,因此文中所提在此基礎(chǔ)上的迭代修正次數(shù)也會(huì)減少,運(yùn)算量也會(huì)隨之降低。

    3.1.2 算法對(duì)信號(hào)載頻的估計(jì)性能的仿真

    這一小節(jié)主要對(duì)算法性能進(jìn)行了仿真分析,因?yàn)槭菍?duì)通信信號(hào)的載頻進(jìn)行估計(jì),因此選用GSM、WCDMA、LTE三種制式的信號(hào),信號(hào)背景噪聲為高斯白噪聲,即信號(hào)傳播信道為AWGN信道,對(duì)背景噪聲大小的描述用SNR(信噪比)來(lái)進(jìn)行描述,采用平均相對(duì)誤差(式(21))來(lái)評(píng)價(jià)載頻估計(jì)的精確度。

    在不同信噪比下分別進(jìn)行100次蒙特卡洛仿真實(shí)驗(yàn),得到兩種算法的估計(jì)平均相對(duì)誤差曲線分別如圖3、圖4和圖5所示:

    圖3 GSM信號(hào)估計(jì)誤差均值曲線

    圖4 WCDMA信號(hào)估計(jì)誤差均值曲線

    圖5 LTE信號(hào)估計(jì)誤差均值曲線

    分析上面三幅圖的結(jié)果可以看出,總體上,所提算法相對(duì)譜重心法的估計(jì)性能表現(xiàn)良好;隨著信噪比的降低,頻率重心法對(duì)信號(hào)載頻的估計(jì)精確度下降較快,在0 dB到10 dB區(qū)間估計(jì)結(jié)果的精確度急劇下降,文中所提改進(jìn)算法對(duì)載頻估計(jì)精度提升明顯,尤其當(dāng)信噪比在0 dB以下的時(shí)候,所提算法的估計(jì)精確度仍然比較高,較原來(lái)算法估計(jì)精度有所提升,證明算法的抗噪性能良好。通過(guò)上面的實(shí)驗(yàn)仿真可以看出文中所提算法對(duì)通信信號(hào)的載頻估計(jì)準(zhǔn)確度較高而且非常有效。

    4 結(jié)束與展望

    本文根據(jù)信號(hào)監(jiān)測(cè)系統(tǒng)的需要,針對(duì)三種制式的通信信號(hào),提出了一種在缺乏先驗(yàn)知識(shí)的情況下,準(zhǔn)確快速估計(jì)信號(hào)的載頻的方法。利用平滑和高斯擬合對(duì)信號(hào)進(jìn)行預(yù)處理,結(jié)合信號(hào)功率譜的對(duì)稱性對(duì)載頻估計(jì)結(jié)果進(jìn)行迭代修正,使算法在低信噪比下也有良好的性能,能夠滿足監(jiān)測(cè)系統(tǒng)在全盲狀態(tài)下對(duì)信號(hào)載頻估計(jì)的需要。在以后的研究中,還會(huì)對(duì)步長(zhǎng)和門限動(dòng)態(tài)調(diào)整進(jìn)行研究,更好的適應(yīng)信號(hào)監(jiān)測(cè)系統(tǒng)。

    [1]孫健, 張錦南.頻譜資源利用評(píng)價(jià)體系研究[J].軟件, 2016(2): 17-21.

    [2]孫大飛, 宋鐵成, 吳名, 等.頻譜感知能量檢測(cè)器的性能分析[J].新型工業(yè)化, 2012(10).

    [3]Cai H.Fast Frequency Measurement Algorithm Based On Zero Crossing Method[C]// Computer Engineering and Technology (ICCET), 2010 2nd International Conference on.2010: V4-606-V4-608.

    [4]Wei J.Modification of frequency estimation algorithms for sinusoidal signals based on phase difference of overlap FFT[C]// International Conference on Communications, Circuits and Systems.IEEE, 2008: 927-929.

    [5]Hill D A, Bodie J B.Experimental carrier detection of BPSK and QPSK direct sequence spread spectrum signals[C]// Military Communications Conference, 1995.MILCOM '95, Conference Record, IEEE.1995: 362-367 vol.1.

    [6]Hui L, Dai B Q, Wei L.A Pitch Detection Algorithm Based on AMDF and ACF[C]// IEEE International Conference on Acoustics.2006: I-I.

    [7]黃春琳, 柳征, 姜文利, 等.基于循環(huán)譜包絡(luò)的擴(kuò)譜直序信號(hào)的碼片時(shí)寬、載頻、幅度估計(jì)[J].電子學(xué)報(bào), 2002(09): 1353-1356.

    [8]鄭文秀, 趙國(guó)慶, 羅明.基于高階循環(huán)累積量的OFDM子載波盲估計(jì)[J].電子與信息學(xué)報(bào), 2008(02): 346-349.

    [9]唐沖, 惠輝輝.基于Matlab的高斯曲線擬合求解[J].計(jì)算機(jī)與數(shù)字工程, 2013, 41(8): 1262-1263.

    [10]劉明騫, 李兵兵, 王婧舒, 等.低信噪比下低復(fù)雜度的OFDM信號(hào)帶寬盲估計(jì)方法[J].江蘇大學(xué)學(xué)報(bào)(自然科學(xué)版), 2013, 34(01): 76-80.

    High-Precision Carrier Frequency Blind Estimation Algorithm for Communication Signals with Low SNR

    ZHANG Ya-feng1, WANG Xiang2, ZHAO Cheng-lin1
    (1.Beijing University of Posts and Telecommunications, Department of Communication and Information System, Beijing 100876; 2.Beijing jianea technology Co., Ltd., Beijng 100088)

    The carrier frequency estimation of the communication signal of unknown standard is an important problem in the station monitoring and spectrum management.Considering the shortcomings that existing carrier frequency estimation algorithms do not have both high accuracy and low complexity, so aiming at the symmetry of the power spectral of the linear modulated signals, an improved estimation algorithm based on the spectral center of gravity method is presented.First, the power spectrum of the signal is smoothed and pretreated by Gaussian curve fitting.Then, the carrier frequency can be estimated accurately by the proposesd method of iterative search.The algorithm can realize blind estimation of carrier frequency without knowing the signal format and other information.The simulation results indicate that the algorithm can realize fast estimation of carrier frequency in Gaussian channel environment, as well as more accurate in low SNR.In summary, it has a promising application prospect in the station monitoring.

    Carrier frequency estimation; Frequency mediacy method; Gaussian fitting; Iterative search

    TN911

    A

    10.3969/j.issn.1003-6970.2017.03.017

    張亞豐(1993-),男,碩士研究生,主要研究方向?yàn)闊o(wú)線信號(hào)處理,參數(shù)估計(jì);王翔(1969-),男,主要研究方向?yàn)轭l譜監(jiān)測(cè),網(wǎng)絡(luò)通信。

    趙成林(1964-),男,教授,主要研究方向?yàn)闊o(wú)線信號(hào)處理與信號(hào)檢測(cè)。

    本文著錄格式:張亞豐,王翔,趙成林.低信噪比下高精度通信信號(hào)載頻盲估計(jì)算法[J].軟件,2017,38(3):82-86

    猜你喜歡
    高斯復(fù)雜度信噪比
    小高斯的大發(fā)現(xiàn)
    基于深度學(xué)習(xí)的無(wú)人機(jī)數(shù)據(jù)鏈信噪比估計(jì)算法
    天才數(shù)學(xué)家——高斯
    一種低復(fù)雜度的慣性/GNSS矢量深組合方法
    低信噪比下LFMCW信號(hào)調(diào)頻參數(shù)估計(jì)
    低信噪比下基于Hough變換的前視陣列SAR稀疏三維成像
    求圖上廣探樹(shù)的時(shí)間復(fù)雜度
    某雷達(dá)導(dǎo)51 頭中心控制軟件圈復(fù)雜度分析與改進(jìn)
    保持信噪比的相位分解反褶積方法研究
    出口技術(shù)復(fù)雜度研究回顧與評(píng)述
    91精品三级在线观看| 中亚洲国语对白在线视频| 亚洲欧美精品综合一区二区三区| 国产区一区二久久| 肉色欧美久久久久久久蜜桃| 日本五十路高清| 亚洲va日本ⅴa欧美va伊人久久 | 菩萨蛮人人尽说江南好唐韦庄| 一本一本久久a久久精品综合妖精| 亚洲精品中文字幕在线视频| 久久人人爽av亚洲精品天堂| 日本一区二区免费在线视频| 天天添夜夜摸| 亚洲 国产 在线| 国精品久久久久久国模美| 亚洲性夜色夜夜综合| 99re6热这里在线精品视频| 天天影视国产精品| 久久久久久久大尺度免费视频| 午夜日韩欧美国产| 日韩中文字幕视频在线看片| 欧美日韩福利视频一区二区| 99国产综合亚洲精品| 精品一品国产午夜福利视频| 飞空精品影院首页| 在线十欧美十亚洲十日本专区| 欧美成狂野欧美在线观看| 欧美午夜高清在线| 欧美+亚洲+日韩+国产| 黄色视频在线播放观看不卡| 丰满迷人的少妇在线观看| 成人手机av| 美女高潮到喷水免费观看| 亚洲少妇的诱惑av| 97精品久久久久久久久久精品| 国产成人a∨麻豆精品| 美女扒开内裤让男人捅视频| 一区二区三区精品91| 亚洲av欧美aⅴ国产| 久久女婷五月综合色啪小说| 一区福利在线观看| 男女午夜视频在线观看| 日韩中文字幕欧美一区二区| 两个人看的免费小视频| 免费观看人在逋| 少妇 在线观看| 1024视频免费在线观看| 亚洲欧美激情在线| 女人精品久久久久毛片| 免费少妇av软件| 国产精品 欧美亚洲| 精品人妻熟女毛片av久久网站| 99国产精品一区二区三区| 别揉我奶头~嗯~啊~动态视频 | 岛国毛片在线播放| 国产精品熟女久久久久浪| 在线看a的网站| 青春草视频在线免费观看| 一二三四在线观看免费中文在| 我要看黄色一级片免费的| av电影中文网址| 美女主播在线视频| 无遮挡黄片免费观看| 一区在线观看完整版| 2018国产大陆天天弄谢| 久久99一区二区三区| 国产精品 国内视频| 美女大奶头黄色视频| 久久精品成人免费网站| 男女免费视频国产| 久久久国产精品麻豆| 国产在线免费精品| 日韩熟女老妇一区二区性免费视频| 午夜福利在线免费观看网站| 亚洲第一青青草原| 午夜免费观看性视频| 他把我摸到了高潮在线观看 | 日韩一区二区三区影片| 久热爱精品视频在线9| 一本大道久久a久久精品| av电影中文网址| 不卡av一区二区三区| 亚洲国产欧美网| 91精品伊人久久大香线蕉| 大片免费播放器 马上看| 欧美精品av麻豆av| 久久久久久亚洲精品国产蜜桃av| 精品国产一区二区久久| 久久精品亚洲av国产电影网| 亚洲久久久国产精品| 亚洲三区欧美一区| www.熟女人妻精品国产| 久久天躁狠狠躁夜夜2o2o| avwww免费| 亚洲精品成人av观看孕妇| 国产人伦9x9x在线观看| 在线观看免费日韩欧美大片| 国产欧美日韩一区二区三 | 国产免费一区二区三区四区乱码| 大香蕉久久成人网| 19禁男女啪啪无遮挡网站| 久久久久久人人人人人| 天堂8中文在线网| 成人手机av| 黄色 视频免费看| av又黄又爽大尺度在线免费看| 国产欧美日韩综合在线一区二区| 国产成人影院久久av| 久久久久久久久久久久大奶| 日韩欧美免费精品| 女性生殖器流出的白浆| 亚洲第一青青草原| 热99久久久久精品小说推荐| 久久亚洲精品不卡| 欧美在线黄色| 亚洲国产看品久久| 国产主播在线观看一区二区| 免费黄频网站在线观看国产| 男人添女人高潮全过程视频| 国产日韩欧美亚洲二区| 日韩欧美国产一区二区入口| 亚洲精品中文字幕在线视频| 一边摸一边做爽爽视频免费| 激情视频va一区二区三区| 久久久精品区二区三区| 亚洲精品国产av蜜桃| 999久久久国产精品视频| 欧美+亚洲+日韩+国产| 十八禁高潮呻吟视频| 亚洲精品国产av蜜桃| 亚洲欧美色中文字幕在线| 12—13女人毛片做爰片一| 美女视频免费永久观看网站| 天天躁夜夜躁狠狠躁躁| 精品亚洲成国产av| 新久久久久国产一级毛片| 精品视频人人做人人爽| 国产主播在线观看一区二区| 欧美激情久久久久久爽电影 | 亚洲 国产 在线| 亚洲国产欧美日韩在线播放| 老司机亚洲免费影院| 国产精品久久久久久精品电影小说| 欧美人与性动交α欧美精品济南到| 成人影院久久| 热99国产精品久久久久久7| 99国产极品粉嫩在线观看| 欧美激情极品国产一区二区三区| 久久久久国内视频| 亚洲精品日韩在线中文字幕| 久久久国产成人免费| 夜夜夜夜夜久久久久| 色综合欧美亚洲国产小说| 亚洲国产欧美在线一区| 国产片内射在线| 天堂8中文在线网| 在线观看免费高清a一片| 天天躁日日躁夜夜躁夜夜| 在线十欧美十亚洲十日本专区| 亚洲国产精品成人久久小说| 亚洲国产精品成人久久小说| 国产在线一区二区三区精| 亚洲国产毛片av蜜桃av| av在线老鸭窝| av欧美777| 国产av一区二区精品久久| 欧美激情 高清一区二区三区| 久久久欧美国产精品| 热99久久久久精品小说推荐| 国产一区二区三区综合在线观看| 国产精品九九99| 国产成人精品在线电影| 免费女性裸体啪啪无遮挡网站| 欧美在线一区亚洲| 18在线观看网站| 国产成人一区二区三区免费视频网站| 亚洲伊人久久精品综合| 日韩中文字幕视频在线看片| 亚洲欧洲日产国产| 国产高清国产精品国产三级| 精品国产乱子伦一区二区三区 | 欧美97在线视频| 黄片小视频在线播放| 19禁男女啪啪无遮挡网站| 美女脱内裤让男人舔精品视频| 天天操日日干夜夜撸| 欧美日韩成人在线一区二区| 免费在线观看黄色视频的| 欧美黄色淫秽网站| 午夜久久久在线观看| 啦啦啦中文免费视频观看日本| 真人做人爱边吃奶动态| 久久精品亚洲熟妇少妇任你| 下体分泌物呈黄色| 丝袜在线中文字幕| 亚洲国产毛片av蜜桃av| 国产免费av片在线观看野外av| 国产一区二区三区在线臀色熟女 | 他把我摸到了高潮在线观看 | 亚洲欧美色中文字幕在线| 一区二区日韩欧美中文字幕| 捣出白浆h1v1| 男女之事视频高清在线观看| 十八禁人妻一区二区| 中文字幕人妻熟女乱码| 国产精品亚洲av一区麻豆| 伦理电影免费视频| 亚洲色图综合在线观看| 亚洲一区中文字幕在线| 国产高清视频在线播放一区 | 超色免费av| 最黄视频免费看| 18禁裸乳无遮挡动漫免费视频| 久久久久精品国产欧美久久久 | 人人妻,人人澡人人爽秒播| 捣出白浆h1v1| 免费在线观看完整版高清| 国产精品影院久久| 男人爽女人下面视频在线观看| 亚洲色图综合在线观看| 亚洲欧美成人综合另类久久久| 香蕉丝袜av| 少妇被粗大的猛进出69影院| 波多野结衣av一区二区av| 波多野结衣av一区二区av| 在线观看免费视频网站a站| 老熟妇仑乱视频hdxx| www日本在线高清视频| 亚洲欧洲精品一区二区精品久久久| 免费观看a级毛片全部| 国产有黄有色有爽视频| 亚洲第一青青草原| 高清在线国产一区| 亚洲精品国产一区二区精华液| 青春草亚洲视频在线观看| 国精品久久久久久国模美| 色婷婷久久久亚洲欧美| 国产精品一二三区在线看| tube8黄色片| 十八禁网站网址无遮挡| 性色av乱码一区二区三区2| 亚洲精品自拍成人| 亚洲国产av影院在线观看| 我的亚洲天堂| 精品福利观看| av在线app专区| 久久久国产精品麻豆| 日本av免费视频播放| 18禁裸乳无遮挡动漫免费视频| 黄片播放在线免费| 国产精品1区2区在线观看. | 精品亚洲成a人片在线观看| 日韩电影二区| 久久久欧美国产精品| 欧美日韩中文字幕国产精品一区二区三区 | 久久 成人 亚洲| 下体分泌物呈黄色| 涩涩av久久男人的天堂| 国产一卡二卡三卡精品| 久久精品人人爽人人爽视色| 最新的欧美精品一区二区| 精品国产一区二区久久| 成人影院久久| 王馨瑶露胸无遮挡在线观看| 国产视频一区二区在线看| 男女午夜视频在线观看| 精品乱码久久久久久99久播| 啦啦啦 在线观看视频| 精品福利永久在线观看| 午夜福利在线观看吧| 国产成人欧美在线观看 | 欧美国产精品一级二级三级| 中文字幕人妻熟女乱码| 久久性视频一级片| 啪啪无遮挡十八禁网站| 国产av精品麻豆| 亚洲国产日韩一区二区| 欧美老熟妇乱子伦牲交| 日韩欧美一区视频在线观看| 亚洲伊人色综图| 在线亚洲精品国产二区图片欧美| 一二三四在线观看免费中文在| 97人妻天天添夜夜摸| 不卡av一区二区三区| 99久久综合免费| 色视频在线一区二区三区| 精品一区在线观看国产| 久久久精品94久久精品| 婷婷丁香在线五月| 成人国产av品久久久| 涩涩av久久男人的天堂| 国内毛片毛片毛片毛片毛片| 日韩制服骚丝袜av| 91精品伊人久久大香线蕉| 丝瓜视频免费看黄片| 91精品三级在线观看| 久久久久久久久久久久大奶| 最新在线观看一区二区三区| 咕卡用的链子| www.自偷自拍.com| 捣出白浆h1v1| 国产极品粉嫩免费观看在线| 久久久久久久精品精品| 建设人人有责人人尽责人人享有的| 少妇精品久久久久久久| 91成年电影在线观看| 亚洲黑人精品在线| 欧美日韩亚洲国产一区二区在线观看 | 黄色视频不卡| 免费高清在线观看视频在线观看| 黑人欧美特级aaaaaa片| 午夜福利影视在线免费观看| 成人国产一区最新在线观看| h视频一区二区三区| 免费观看人在逋| 777米奇影视久久| 久久人人爽av亚洲精品天堂| 一本大道久久a久久精品| 超碰成人久久| 国产一区有黄有色的免费视频| 最近中文字幕2019免费版| 少妇的丰满在线观看| 大片免费播放器 马上看| 亚洲精品一区蜜桃| 一个人免费在线观看的高清视频 | 久久热在线av| 韩国高清视频一区二区三区| 人妻一区二区av| 中文欧美无线码| 午夜免费鲁丝| 婷婷成人精品国产| 亚洲成人国产一区在线观看| 久久久久国内视频| 亚洲精品久久成人aⅴ小说| 18禁国产床啪视频网站| 男女下面插进去视频免费观看| 亚洲国产av影院在线观看| 91大片在线观看| 国产亚洲av高清不卡| 国产精品国产av在线观看| 视频在线观看一区二区三区| 精品欧美一区二区三区在线| 亚洲成人手机| 国产亚洲精品第一综合不卡| 男女之事视频高清在线观看| 婷婷丁香在线五月| 欧美 亚洲 国产 日韩一| 巨乳人妻的诱惑在线观看| 51午夜福利影视在线观看| 可以免费在线观看a视频的电影网站| 免费在线观看日本一区| 蜜桃国产av成人99| 大香蕉久久网| 亚洲国产av新网站| 亚洲va日本ⅴa欧美va伊人久久 | 大陆偷拍与自拍| 一进一出抽搐动态| 国产1区2区3区精品| 欧美日韩亚洲高清精品| 久久久水蜜桃国产精品网| 麻豆国产av国片精品| 看免费av毛片| 男人舔女人的私密视频| 国产免费福利视频在线观看| 777米奇影视久久| 天堂俺去俺来也www色官网| 日本精品一区二区三区蜜桃| 波多野结衣一区麻豆| 国产在视频线精品| 亚洲成人国产一区在线观看| tocl精华| 丁香六月欧美| 亚洲国产看品久久| 久久久久国内视频| 窝窝影院91人妻| 日本五十路高清| 亚洲激情五月婷婷啪啪| 中文字幕av电影在线播放| 涩涩av久久男人的天堂| 欧美精品一区二区大全| 久久午夜综合久久蜜桃| 久久久久久久大尺度免费视频| 男女床上黄色一级片免费看| 久久天躁狠狠躁夜夜2o2o| 如日韩欧美国产精品一区二区三区| 亚洲精品一二三| 国产成人精品久久二区二区免费| 亚洲一码二码三码区别大吗| 午夜福利视频在线观看免费| 日本欧美视频一区| 中文字幕另类日韩欧美亚洲嫩草| 午夜久久久在线观看| 超碰成人久久| 亚洲男人天堂网一区| 一区二区av电影网| 免费日韩欧美在线观看| 欧美老熟妇乱子伦牲交| 一区二区av电影网| 9热在线视频观看99| 久久亚洲国产成人精品v| 后天国语完整版免费观看| 看免费av毛片| av天堂在线播放| 亚洲国产欧美一区二区综合| 9色porny在线观看| 18禁观看日本| 51午夜福利影视在线观看| 亚洲国产精品一区三区| 十八禁高潮呻吟视频| 亚洲精华国产精华精| 我的亚洲天堂| 亚洲伊人色综图| 国产免费福利视频在线观看| 国产精品一区二区免费欧美 | 最近最新免费中文字幕在线| 黄色毛片三级朝国网站| 精品亚洲成国产av| 这个男人来自地球电影免费观看| 真人做人爱边吃奶动态| 欧美97在线视频| 亚洲色图 男人天堂 中文字幕| 黑人操中国人逼视频| 亚洲第一av免费看| 免费高清在线观看视频在线观看| 9191精品国产免费久久| 热99久久久久精品小说推荐| 最黄视频免费看| 欧美国产精品va在线观看不卡| 精品乱码久久久久久99久播| e午夜精品久久久久久久| 亚洲五月色婷婷综合| 天天添夜夜摸| 国产一区二区 视频在线| 精品亚洲乱码少妇综合久久| 亚洲精品日韩在线中文字幕| 9191精品国产免费久久| 久久天堂一区二区三区四区| 免费不卡黄色视频| 亚洲熟女毛片儿| 国产免费现黄频在线看| 国产精品免费大片| 欧美+亚洲+日韩+国产| cao死你这个sao货| 真人做人爱边吃奶动态| 美女中出高潮动态图| 久久性视频一级片| 黄色毛片三级朝国网站| 精品国产一区二区三区四区第35| 黑人巨大精品欧美一区二区蜜桃| 激情视频va一区二区三区| 精品人妻在线不人妻| 亚洲欧美精品自产自拍| 最新在线观看一区二区三区| 国产免费视频播放在线视频| 亚洲 国产 在线| 国产人伦9x9x在线观看| 国产视频一区二区在线看| 免费观看av网站的网址| 我的亚洲天堂| 精品视频人人做人人爽| 精品一区二区三区四区五区乱码| 欧美成人午夜精品| 99国产精品一区二区蜜桃av | 国产欧美日韩一区二区三 | 久久青草综合色| 日韩电影二区| 免费不卡黄色视频| 欧美日韩黄片免| 麻豆av在线久日| 亚洲精品国产av成人精品| 捣出白浆h1v1| 制服人妻中文乱码| 午夜福利视频在线观看免费| av福利片在线| 一二三四在线观看免费中文在| 黑丝袜美女国产一区| 国产99久久九九免费精品| 久久毛片免费看一区二区三区| 各种免费的搞黄视频| 在线亚洲精品国产二区图片欧美| 电影成人av| 丰满少妇做爰视频| 亚洲精品美女久久久久99蜜臀| 首页视频小说图片口味搜索| 老汉色∧v一级毛片| 久久性视频一级片| 亚洲精品中文字幕在线视频| 制服诱惑二区| 亚洲精品在线美女| 国产一区二区三区av在线| 亚洲精品在线美女| 一级片免费观看大全| 精品久久久精品久久久| 捣出白浆h1v1| 91老司机精品| 国产在线观看jvid| 一区二区av电影网| 日韩 亚洲 欧美在线| 免费女性裸体啪啪无遮挡网站| 精品一区二区三区av网在线观看 | 手机成人av网站| 日韩熟女老妇一区二区性免费视频| 久久影院123| 欧美变态另类bdsm刘玥| 日本撒尿小便嘘嘘汇集6| 亚洲精品国产区一区二| cao死你这个sao货| 搡老熟女国产l中国老女人| 日本精品一区二区三区蜜桃| 日本猛色少妇xxxxx猛交久久| 波多野结衣av一区二区av| 免费在线观看黄色视频的| a级毛片黄视频| 俄罗斯特黄特色一大片| av在线播放精品| 国产91精品成人一区二区三区 | 日韩熟女老妇一区二区性免费视频| 欧美黄色淫秽网站| 免费日韩欧美在线观看| h视频一区二区三区| 成年人黄色毛片网站| 中文字幕高清在线视频| 国产成人系列免费观看| 啪啪无遮挡十八禁网站| 国产精品一区二区免费欧美 | 午夜福利在线免费观看网站| 成年女人毛片免费观看观看9 | 日日夜夜操网爽| 午夜福利视频在线观看免费| 国产成人精品无人区| 久久影院123| 男人舔女人的私密视频| 国内毛片毛片毛片毛片毛片| www.精华液| 午夜精品国产一区二区电影| 丝袜美腿诱惑在线| 亚洲欧洲日产国产| 亚洲精品中文字幕在线视频| 日韩欧美免费精品| 正在播放国产对白刺激| 伊人久久大香线蕉亚洲五| 成年美女黄网站色视频大全免费| 欧美中文综合在线视频| 999久久久精品免费观看国产| 美女脱内裤让男人舔精品视频| 午夜久久久在线观看| 欧美日韩黄片免| 免费观看av网站的网址| svipshipincom国产片| 巨乳人妻的诱惑在线观看| 国产一区二区在线观看av| 亚洲国产成人一精品久久久| 香蕉国产在线看| 欧美日韩亚洲高清精品| 亚洲熟女精品中文字幕| 国产精品1区2区在线观看. | 欧美国产精品一级二级三级| 日韩视频在线欧美| 2018国产大陆天天弄谢| 国产xxxxx性猛交| 一本色道久久久久久精品综合| 国产亚洲精品第一综合不卡| 日韩制服骚丝袜av| 新久久久久国产一级毛片| av片东京热男人的天堂| 波多野结衣av一区二区av| 日本vs欧美在线观看视频| 亚洲国产成人一精品久久久| 国产一区二区三区av在线| 亚洲欧美日韩另类电影网站| 美女大奶头黄色视频| 国产精品1区2区在线观看. | 成人18禁高潮啪啪吃奶动态图| 人妻 亚洲 视频| 欧美另类一区| 两个人免费观看高清视频| 成人影院久久| 欧美黄色片欧美黄色片| 高潮久久久久久久久久久不卡| 美女高潮喷水抽搐中文字幕| 正在播放国产对白刺激| kizo精华| 国产高清国产精品国产三级| a级片在线免费高清观看视频| 国产人伦9x9x在线观看| 免费在线观看完整版高清| 精品第一国产精品| 精品少妇黑人巨大在线播放| 午夜免费鲁丝| 在线天堂中文资源库| 免费在线观看日本一区| 欧美日韩中文字幕国产精品一区二区三区 | 18禁裸乳无遮挡动漫免费视频| 久久久国产欧美日韩av| 久久久久视频综合| 美女中出高潮动态图| 亚洲欧洲日产国产| 1024香蕉在线观看| 黄色a级毛片大全视频| 女人久久www免费人成看片| 在线观看舔阴道视频| 久久人妻熟女aⅴ| 老司机亚洲免费影院| av网站免费在线观看视频| 九色亚洲精品在线播放| 亚洲男人天堂网一区| 制服人妻中文乱码| 99精品欧美一区二区三区四区| 国产日韩一区二区三区精品不卡| 十八禁网站网址无遮挡| 亚洲精品av麻豆狂野| 捣出白浆h1v1| 午夜日韩欧美国产| 欧美日韩亚洲综合一区二区三区_| 欧美性长视频在线观看|