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

    頻域匹配濾波-FFT聯(lián)合捕獲算法*

    2018-01-02 17:25:52歐春湘吳智杰
    現(xiàn)代防御技術(shù) 2017年6期
    關(guān)鍵詞:偽碼頻域時(shí)域

    歐春湘,吳智杰

    (北京遙感設(shè)備研究所,北京 100854)

    0 引言

    由于電磁波的空間傳輸延時(shí),以及收發(fā)機(jī)間相對(duì)運(yùn)動(dòng)所導(dǎo)致的頻率多普勒效應(yīng),無線電信號(hào)接收時(shí)需要進(jìn)行定時(shí)同步和載波同步。捕獲是指在時(shí)間和頻率上分別設(shè)置固定的步進(jìn)對(duì)碼相位和多普勒頻率進(jìn)行二維搜索相關(guān)峰值[1-3]。對(duì)于直擴(kuò)序列的快速捕獲主要有2種技術(shù)途徑。一種是在時(shí)域上進(jìn)行偽碼和載波的串行搜索、時(shí)域匹配濾波器法[2-5];另一種是在頻域上利用快速傅立葉變換(fast fourier transform,F(xiàn)FT)和快速傅立葉逆變換(inverse fast fourier transform,IFFT)計(jì)算多普勒頻移估計(jì)和碼相位估計(jì)。傳統(tǒng)的捕獲方法主要有滑動(dòng)相關(guān)法、時(shí)域匹配濾波法、頻域匹配濾波法、平均相關(guān)捕獲算法、折疊匹配濾波法和PMF-FFT(partial matched filtering)部分匹配濾波法等[6-10]。本文提出的頻域匹配濾波-FFT聯(lián)合捕獲法不僅實(shí)現(xiàn)簡(jiǎn)單且大幅度地提高了多普勒頻率的估計(jì)精度。

    1 頻域匹配濾波

    1.1 基本原理

    假設(shè)偽碼長度為N,接收到的偽碼序列為c(n),本地復(fù)制偽碼序列為c′(n),相關(guān)值為corr(n),則有

    (1)

    對(duì)相關(guān)值corr(n)進(jìn)行離散傅里葉變換:

    (2)

    式中:C(k)和C′(k)分別為c(n)和c′(n)的離散傅里葉變換;C′*(k)代表復(fù)數(shù)C′(k)的共軛。

    C(k)=FFT(c(n)),

    (3)

    C′(k)=FFT(c′(n)).

    (4)

    式(2)表明2個(gè)碼序列c(n)和c′(n)在時(shí)域做相關(guān)運(yùn)算,在頻域內(nèi)相當(dāng)于一個(gè)序列的離散傅里葉變換C(k)與另一個(gè)序列的離散傅里葉變換的共軛值C′(k)相乘。因此CORR(k)的離散傅里葉反變換正好是2個(gè)碼序列在各個(gè)碼相位處的相關(guān)值。

    corr(n)=IFFT(CORR(k)).

    (5)

    1.2 實(shí)現(xiàn)流程

    捕獲過程需要完成載波多普勒和偽碼相位的二維掃描[11],如圖1所示,fbin和tbin分別為頻率搜索步進(jìn)和碼相位搜索步進(jìn)。當(dāng)在第i行第j列得到相關(guān)峰值且超過門限時(shí),則得到與接收信號(hào)一致的載波多普勒fd(i)和碼相位值p(j),完成了碼并行捕獲。如何按照一定的策略快速地遍歷所有的載波多普勒頻點(diǎn)和碼相位是捕獲的關(guān)鍵[12-14]。頻率匹配濾波法就是在每一次掃描一個(gè)載波多普勒頻率,在該頻點(diǎn)按照式(2)~(4)完成偽碼的快速相關(guān),實(shí)現(xiàn)流程圖如圖2所示。

    圖1 捕獲二維搜索圖Fig.1 Acquisition searching map in two dimension

    圖2 頻譜匹配濾波原理圖Fig.2 Construction of frequency domain matching filter

    ui(t)為接收信號(hào),載波NCO產(chǎn)生查找地址,通過查找正余弦查找表得到本地同相載波uoc(t)和正交載波uos(t),分別表示為

    (6)

    (7)

    (8)

    式中:a和a′分別為接收信號(hào)和本地載波的幅值;d(t)為調(diào)制的數(shù)據(jù);c(t)為偽碼序列;ωi和ωo分別為接收到的載波頻率和本地載波頻率;φ為接收信號(hào)的初相位。

    正交下變頻后同相支路再經(jīng)過低通濾波器濾除高頻分量后得到

    i(t)=ui(t)uoc(t)=2aa′d(t)c(t)·

    cos(ωit+φ)cos(ωot)=

    aa′d(t)c(t){cos[(ωi+ωo)t+φ]+

    aa′d(t)c(t)cos[(ωi-ωo)t+φ].

    (9)

    正交下變頻后正交支路再經(jīng)過低通濾波器濾除高頻分量后得到

    q(t)=ui(t)uos(t)=2aa′d(t)c(t)·

    cos(ωit+φ)sin(ωot)=

    aa′d(t)c(t){sin[(ωi+ωo)t+φ]-

    -aa′d(t)c(t)sin[(ωi-ωo)t+φ].

    (10)

    I和Q 2個(gè)支路合并成i(t)+j·q(t)后做FFT運(yùn)算,再與偽碼FFT共軛相乘,做FFT逆變換得到各個(gè)偽碼相位上的相關(guān)值,相關(guān)運(yùn)算后得到

    uout(t)=aa′d(t)R(τ)e-j(Δωt+φ),

    (11)

    式中:R(τ)為接收到的偽碼c(t)與本地偽碼c′(t)的相關(guān)值;τ為c(t)與c′(t)的時(shí)間差;Δω=ωi-ωo為接收到的載波和本地載波的頻率差。

    對(duì)uout(t)取模進(jìn)行幅值判斷,若大于門限值則表示捕獲成功,否則繼續(xù)下一個(gè)頻點(diǎn)重復(fù)上述過程。因?yàn)閭未a的強(qiáng)相關(guān)性和弱互相關(guān)性,當(dāng)τ不為零時(shí),R(τ)將嚴(yán)重衰減使得幅值減小;當(dāng)Δω不為零時(shí),由指數(shù)函數(shù)的衰減性可知,幅值將成指數(shù)衰減[15-16]。因此由式(10)可知,只有當(dāng)載波和碼相位都嚴(yán)格對(duì)齊時(shí)才能得到峰值,捕獲到峰值的同時(shí)得到了載波多普勒值和偽碼相位值。

    2 聯(lián)合捕獲法

    聯(lián)合捕獲法是基于頻域匹配濾波法的一種捕獲方法,連續(xù)k個(gè)碼周期做相關(guān)運(yùn)算,對(duì)同相位的k

    個(gè)相關(guān)值做二次FFT運(yùn)算,具體的實(shí)現(xiàn)方法如圖3所示。

    第1步,利用FFT運(yùn)算求解輸入信號(hào)的頻域表示,F(xiàn)FT的頻移特性即頻域移相等效時(shí)域相乘的特性如式(12)所示,因此在各支路上對(duì)N1點(diǎn)輸出進(jìn)行2N2+1次移相,等效于在多普勒頻率范圍[-Fd,F(xiàn)d](kHz)內(nèi)對(duì)輸入信號(hào)進(jìn)行2N2+1路并行下變頻。

    若FFT[x(n)]=X(k),則

    (12)

    相關(guān)積分時(shí)間Tcoh為一個(gè)碼周期時(shí)間,將一個(gè)碼周期劃分成N1個(gè)搜索單元,

    (13)

    式中:fs為采樣頻率;Rc為碼速率;L為碼周期。

    圖3 頻域匹配濾波-FFT聯(lián)合捕獲法實(shí)現(xiàn)流程Fig.3 Construction of frequency domain matching filter and FFT union acquisition

    (14)

    將多普勒頻率范圍[-Fd,F(xiàn)d](kHz)分成2N2+1個(gè)搜索單元,每個(gè)搜索頻率記為fdi,本地載波記為CLi(t),即

    (15)

    fdi=i·fbin,

    i=-N2,-(N2-1),…,0,…,N2,

    (16)

    CLi(t)=cos(2πfdit)-i·sin(2πfdit),

    i=-N2,-(N2-1),…,0,…,N2,

    (17)

    式中:[·]表示取整數(shù)。

    假設(shè)輸入的基帶信號(hào)為S1(t),S1(t)與CLi(t)正交下變頻得到基帶復(fù)信號(hào)S2i(t)。

    i=-N2,-(N2-1),…,0,…,N2,

    (18)

    i=-N2,-(N2-1),…,0,…,N2.

    (19)

    第2步,2N2+1路下變頻信號(hào)S2i(t)與原位輸出的本地偽碼的FFT運(yùn)算結(jié)果共軛相乘,并對(duì)輸出結(jié)果進(jìn)行IFFT運(yùn)算,由式(2)可知,該步驟相當(dāng)于求得各通道上的一個(gè)碼周期長度的時(shí)域相關(guān)積分,得到N1(2N2+1)個(gè)相關(guān)積分值。至此,完成了單個(gè)偽碼周期內(nèi)的掃頻搜索,頻率搜索步進(jìn)為fbin。

    2N2+1通道的第k通道零中頻基帶信號(hào)S2k(t)與本地偽碼pn(t)相關(guān)累加后得到N1個(gè)相關(guān)值S3k(l),采用MSK調(diào)制。

    (20)

    S3k(l)=S2k(t)?pn(t+lTc)=

    i(t)?pn(t+lTc)+j·q(t)?pn(t+lTc)=

    I(l)+j·Q(l),

    l=0,1,…,N1-1,

    (21)

    式中:Tc=TL/N1,TL為一個(gè)碼周期時(shí)間。

    I(l)=i(t)?pn(t+lTc)=

    l=0,1,…,N1-1,

    (22)

    Q(l)=q(t)?pn(t++lTc)=

    l=0,1,…,N1-1,

    (23)

    fe=fd-fdi,

    (24)

    S3k(l)=I(l)+j·Q(l)=

    l=0,1,…,N1-1,

    (25)

    式中:t1為相關(guān)起始時(shí)刻;ωe=2πfe;R(·)表示自相關(guān)函數(shù)。

    第3步,在連續(xù)K個(gè)偽碼周期上做上述相關(guān)運(yùn)算,每個(gè)通道得到KN1個(gè)相關(guān)值,將KN1值存儲(chǔ)在一個(gè)K×N1的矩陣X1里。對(duì)每一列K個(gè)點(diǎn)進(jìn)行FFT運(yùn)算,得到K×N1的矩陣X2,等效為對(duì)各相位上的相關(guān)結(jié)果按基帶速率進(jìn)行K次相干累加,理論捕獲精度提高為fbin/K。

    矩陣X2的第p列為S3k(p),F(xiàn)FT峰值為v,則

    v=K·S3k(p)=

    (26)

    τp+pTc≈0,

    (27)

    fe=fd-fdi≈0.

    (28)

    式(27),(28)表示只有當(dāng)碼相位差小于碼相位搜索步徑,且本地多普勒頻率與接收信號(hào)的多普勒頻率基本保持一致時(shí)才會(huì)出現(xiàn)相關(guān)峰值。

    第4步,在總樣本為(2N2+1)N1K的基礎(chǔ)上,計(jì)算相關(guān)值的最大值和均值,并得到峰均比。至此,完成了第1輪K個(gè)偽碼周期的相干累積搜索。

    第5步,當(dāng)連續(xù)2次判決輸出均過門限Vt,2次捕獲結(jié)果碼片差在碼相位搜索步進(jìn)以內(nèi),且頻率差在fbin/K以內(nèi),則判決成功,啟動(dòng)后續(xù)信號(hào)處理模塊,否則,捕獲模塊繼續(xù)搜索。

    3 仿真與分析

    與理論分析部分的假設(shè)相對(duì)應(yīng),設(shè)置基本的數(shù)值和仿真環(huán)境為:采樣頻率fs=4 MHz,碼速率Rc=2 MHz,碼周期L=128,多普勒頻偏范圍±100 kHz,聯(lián)合捕獲的頻率搜索步進(jìn)為15.625 kHz,碼片搜索步進(jìn)為0.5個(gè)碼片,相關(guān)時(shí)間為一個(gè)碼周期,在連續(xù)K=20個(gè)偽碼周期上做相關(guān)運(yùn)算,每個(gè)碼相位上得到20個(gè)相關(guān)值,再補(bǔ)10個(gè)0做32點(diǎn)FFT運(yùn)算。FFT運(yùn)算提高了相關(guān)峰均值,多普勒頻率估計(jì)精度提高為15.625 kHz/32=488.28 Hz,設(shè)置以下2個(gè)仿真條件分別進(jìn)行驗(yàn)證。

    設(shè)置載波多普勒頻率為66.7 kHz,碼相位差為155個(gè)半碼片,不存在跳頻同步誤差。圖4為一個(gè)偽碼周期的相關(guān)峰均值,圖5為聯(lián)合捕獲最后一級(jí)32點(diǎn)FFT得到的相關(guān)峰均值,對(duì)比可知,聯(lián)合捕獲得到的相關(guān)峰均值明顯高于直接相關(guān)累加的相關(guān)峰均值,提高了檢測(cè)概率。

    圖4 頻域匹配濾波相關(guān)峰值Fig.4 Correlation result of frequency domain matching filter

    圖5 聯(lián)合捕獲相關(guān)峰值Fig.5 Correlation result of union acquisition

    設(shè)置中頻信噪比為-15 dB,多普勒頻率在±100 kHz范圍內(nèi),通過仿真得到聯(lián)合捕獲法和頻率匹配濾波法的載波多普勒頻率估計(jì)誤差的變化情況,如圖6所示。由圖6可知,頻率匹配濾波法殘余的多普勒頻偏精度為8 kHz,聯(lián)合捕獲殘余的為1 kHz,在捕獲精度上明顯優(yōu)于傳統(tǒng)頻域匹配濾波法。

    圖6 多普勒頻率估計(jì)誤差Fig.6 Estimated error of Doppler frequency

    為了得到在±100 kHz多普勒頻偏范圍內(nèi),捕獲概率隨信噪比的變化,設(shè)置多普勒頻偏分別為-100,0,100 kHz,得到3條捕獲概率隨中頻輸入信噪比變化的曲線如圖7所示。由圖7可知,隨著信噪比的增加捕獲概率不斷增加,由于位于頻率搜索單元的不同頻率位置,0和±100 kHz存在2 dB的差別,多普勒頻偏為0時(shí)在-16 dB捕獲概率接近1,多普勒頻偏為±100 kHz時(shí)在-14 dB捕獲概率接近1。

    圖7 不同信噪比下的捕獲概率Fig.7 Acquisition probability in different SNR

    4 結(jié)束語

    頻域匹配濾波-FFT聯(lián)合捕獲法在頻域匹配濾波的基礎(chǔ)上對(duì)連續(xù)K個(gè)相同碼相位的相關(guān)結(jié)果進(jìn)行二次FFT運(yùn)算,相當(dāng)于K次非相干累加。經(jīng)過實(shí)驗(yàn)仿真可知,該方法提高了多普勒頻率的估計(jì)精度,在信噪比大于-14 dB時(shí)捕獲概率能達(dá)到0.988,多普勒頻率的捕獲范圍能夠達(dá)到±100 kHz。由此可見,頻域匹配濾波-FFT聯(lián)合捕獲法在估計(jì)精度和捕獲概率上優(yōu)于傳統(tǒng)的頻域匹配濾波法。

    [1] 王江浩,李署堅(jiān).一種改進(jìn)的偽碼捕獲方法[J].電子測(cè)量技術(shù),2008,31(12):10-14.

    WANG Jiang-hao,LI Shu-jian.Improved Method of Acquisition of Pseudo[J].Electronic Measurement Technology,2008,31(12):10-14.

    [2] 展放,張君昌.數(shù)字匹配濾波器在直序擴(kuò)頻快速捕獲中的應(yīng)用[J].航空電子技術(shù),2003,34(3):15-18.

    ZHAN Fang,ZHANG Jun-chang.Application of DFM to DS/SS Rapid Acquisition[J].Avionics Technology,2003,34(3):15-18.

    [3] 金俊坤,吳嗣亮,李菊.基于基2-FFT的偽碼快速捕獲實(shí)現(xiàn)新算法[J].系統(tǒng)工程與電子技術(shù),2005,27(11):1957-1960.

    JIN Jun-kun,WU Si-liang,LI Ju.Implementation of A New Fast PN Code-Acquisition Using Radix-2 FFT[J].Systems Engineering and Electronics,2005,27(11):1957-1960.

    [4] 胡建波,楊莘元,盧滿宏.一種基于FFT的高動(dòng)態(tài)擴(kuò)頻信號(hào)的快速捕獲方法[J].遙測(cè)技術(shù),2004,25(6):19-24.

    HU Jian-bo,YANG Xin-yuan,LU Man-hong.An FFT Based Method for Fast Acquisition of High Dynamic DS Spread Spectrum Signals[J].Telemetry & Telecontrol,2004,25(6):19-24.

    [5] 國世超,陳安樂,韓方劍.基于部分匹配濾波的偽碼捕獲方法研究[J].電子信息對(duì)抗技術(shù),2010,25(2):20-23.

    GUO Shi-chao,CHEN An-le,HAN Fang-jian.Pseudo-Noise Code Acquisition Method Based on Partial Matched Filter[J].Electronic Information Warfare Technology,2010,25(2):20-23.

    [6] 黃振,陸建華,楊士中.基于DMF直擴(kuò)系統(tǒng)捕獲性能研究[J].電路與系統(tǒng)學(xué)報(bào),2002,7(1):92-95.

    HUANG Zhen,LU Jian-hua,YANG Shi-zhong.Perfor-mance of Acquisition in a Matched-Ffilter for DSSS[J].Journal of Circuits and Systems,2002,7(1):92-95.

    [7] 黃磊,張其善,寇艷紅.基于差分碼和塊處理方法的GPS信號(hào)快速捕獲算法[J].通信學(xué)報(bào),2006,27(8):14-19.

    HUANG Lei,ZHANG Qi-shan,KOU Yan-hong.GPS Signal Fast Acquisition Algorithm Using Differential Code and Block Processing Techniques[J].Journal on Communications,2006,27(8):14-19.

    [8] LIN D,TSUI J.Comparison of Acquisition Methods for Software GPS Receiver[J].Ion GPS,2002,91(5):2385-2390.

    [9] LEE Yong-hwan,TANTARATANA S.Sequential Acqui-sition of PN Sequences for DS/ SS Communications:Design and Performance[J].IEEE Journal on Select Areas in Communications,1992,10(4):750- 759.

    [10] TANIGAWA M,HARA T,F(xiàn)UKUI M,et al.Fast Acquisition and Low Correlation Multi-Component Code[J].IEEE Transactions on Communications,2004,4(1):2454-2457.

    [11] FENG G.Block Processing Techniques for the Global Positioning System[D].Ohio:Ohio University,2003.

    [12] ZHU Z.Averaging Correlation for Weak Global Positioning System Signal Processing[D].Ohio:Ohio University,2002.

    [13] STARZYK J A,ZHU Z.Averaging Correlation for C/A code Acquisition and Tracking in Frequency Domain[J].MWSCS Conf,2001,12(15):24-28.

    [14] IUNTA G,BENEDETTO F.Spread-Spectrum Code Acquisition in The Presence of Cell Correlation[J].IEEE Transactions on Communications,2007,55 (2):257-261.

    [15] JARI H,IINATTI J.On the Threshold Setting Principles in Code Acquisition of DS-SS Signals[J].IEEE Journal on Selected Areas in Communications,2000,8(1):62-72.

    [16] SIMON M K,OMURA J K,Scholtz R A,et al.Spread Spectrum Communication Handbook[M].Beijing:Post & Telecommunications Press,2002.

    猜你喜歡
    偽碼頻域時(shí)域
    非協(xié)作多用戶短碼直擴(kuò)信號(hào)偽碼估計(jì)
    基于時(shí)域信號(hào)的三電平逆變器復(fù)合故障診斷
    頻域稀疏毫米波人體安檢成像處理和快速成像稀疏陣列設(shè)計(jì)
    基于極大似然準(zhǔn)則與滾動(dòng)時(shí)域估計(jì)的自適應(yīng)UKF算法
    基于改進(jìn)Radon-Wigner變換的目標(biāo)和拖曳式誘餌頻域分離
    基于時(shí)域逆濾波的寬帶脈沖聲生成技術(shù)
    一種基于頻域的QPSK窄帶干擾抑制算法
    偽碼體制脈沖串引信信號(hào)參數(shù)提取方法
    基于頻域伸縮的改進(jìn)DFT算法
    基于時(shí)域波形特征的輸電線雷擊識(shí)別
    免费观看在线日韩| 欧美日本视频| 色综合婷婷激情| 亚洲精品一卡2卡三卡4卡5卡| 国产免费一级a男人的天堂| www.色视频.com| 九色成人免费人妻av| 成人特级av手机在线观看| 午夜免费激情av| 亚洲精品456在线播放app | 婷婷亚洲欧美| 久久久久久九九精品二区国产| 黄色女人牲交| 俄罗斯特黄特色一大片| 18禁黄网站禁片午夜丰满| 高清在线国产一区| 国产精品亚洲美女久久久| 色哟哟哟哟哟哟| 国产午夜精品久久久久久一区二区三区 | 丝袜美腿在线中文| 白带黄色成豆腐渣| a级一级毛片免费在线观看| 国产精品久久久久久久久免| 亚洲美女搞黄在线观看 | 国产精品三级大全| 免费观看人在逋| 国产人妻一区二区三区在| 久久午夜福利片| 一个人看的www免费观看视频| 麻豆精品久久久久久蜜桃| 亚洲精华国产精华精| 91久久精品国产一区二区三区| 亚洲精品粉嫩美女一区| 露出奶头的视频| 国产欧美日韩精品亚洲av| 色播亚洲综合网| 国产美女午夜福利| 亚洲在线观看片| 99久久九九国产精品国产免费| 日韩精品青青久久久久久| 久久久久久伊人网av| 狠狠狠狠99中文字幕| 午夜爱爱视频在线播放| 免费观看精品视频网站| 国产极品精品免费视频能看的| 黄片wwwwww| 成年女人永久免费观看视频| 久久九九热精品免费| 成人亚洲精品av一区二区| 最近最新中文字幕大全电影3| 欧美激情久久久久久爽电影| 少妇人妻一区二区三区视频| 欧美一区二区亚洲| 最近在线观看免费完整版| 能在线免费观看的黄片| 国产探花在线观看一区二区| 国产精品综合久久久久久久免费| 美女高潮喷水抽搐中文字幕| 禁无遮挡网站| 国产高清三级在线| 日韩中文字幕欧美一区二区| 丰满人妻一区二区三区视频av| 国产精品三级大全| 亚洲成人久久爱视频| 99国产精品一区二区蜜桃av| 日日夜夜操网爽| 在线免费观看不下载黄p国产 | 很黄的视频免费| 国产人妻一区二区三区在| 1024手机看黄色片| av女优亚洲男人天堂| 91麻豆av在线| 看免费成人av毛片| 一区福利在线观看| 亚洲黑人精品在线| 日本在线视频免费播放| 永久网站在线| 少妇的逼好多水| 嫁个100分男人电影在线观看| 欧洲精品卡2卡3卡4卡5卡区| 变态另类成人亚洲欧美熟女| 美女黄网站色视频| 99久久精品国产国产毛片| 别揉我奶头 嗯啊视频| 伊人久久精品亚洲午夜| 美女免费视频网站| 欧美黑人巨大hd| 免费黄网站久久成人精品| 一个人看视频在线观看www免费| 又粗又爽又猛毛片免费看| 国产精品久久久久久精品电影| 国产高清有码在线观看视频| 最新在线观看一区二区三区| 亚洲午夜理论影院| 波多野结衣高清无吗| 波多野结衣高清作品| 我要搜黄色片| 他把我摸到了高潮在线观看| 成人亚洲精品av一区二区| 午夜激情福利司机影院| 久久午夜福利片| 国模一区二区三区四区视频| 少妇猛男粗大的猛烈进出视频 | 国产伦精品一区二区三区四那| 国产在视频线在精品| 十八禁国产超污无遮挡网站| 国产亚洲精品久久久久久毛片| 香蕉av资源在线| 免费人成视频x8x8入口观看| 日本在线视频免费播放| 国产高清不卡午夜福利| 色综合亚洲欧美另类图片| 精品久久久久久久末码| av在线老鸭窝| 悠悠久久av| 国产av不卡久久| 亚洲四区av| 国产高清三级在线| 日本黄大片高清| 国产高清视频在线观看网站| 99精品久久久久人妻精品| 午夜精品一区二区三区免费看| 成人永久免费在线观看视频| 免费一级毛片在线播放高清视频| 亚洲最大成人av| 男女下面进入的视频免费午夜| 亚洲自拍偷在线| avwww免费| 丰满人妻一区二区三区视频av| 我要搜黄色片| 午夜精品在线福利| 春色校园在线视频观看| 一夜夜www| 校园人妻丝袜中文字幕| 久久久久国内视频| 日本撒尿小便嘘嘘汇集6| 亚洲美女搞黄在线观看 | 久久久久久久精品吃奶| 日本一二三区视频观看| 亚洲四区av| 亚洲成人久久爱视频| 久久精品国产亚洲av天美| 久久精品国产鲁丝片午夜精品 | 亚洲va日本ⅴa欧美va伊人久久| 久久久久久九九精品二区国产| 亚洲自拍偷在线| 丰满的人妻完整版| 国产精品爽爽va在线观看网站| 一区二区三区免费毛片| 美女高潮的动态| 亚洲av成人精品一区久久| 亚洲在线观看片| 日韩欧美国产在线观看| 久久久久免费精品人妻一区二区| 在线天堂最新版资源| 免费看a级黄色片| 国产一区二区三区在线臀色熟女| 久久久久久久精品吃奶| 成人无遮挡网站| 简卡轻食公司| 精品久久久久久久久久免费视频| 欧美激情国产日韩精品一区| 最好的美女福利视频网| 淫秽高清视频在线观看| 久久亚洲真实| 成人鲁丝片一二三区免费| 国产综合懂色| 欧美激情国产日韩精品一区| 亚洲一区二区三区色噜噜| 99久久久亚洲精品蜜臀av| 国产一区二区三区av在线 | 国产精品一区二区三区四区免费观看 | 变态另类丝袜制服| 日韩欧美三级三区| 欧美一区二区国产精品久久精品| www.www免费av| 国产高清三级在线| 日韩精品青青久久久久久| 18禁在线播放成人免费| 欧美高清性xxxxhd video| 国产精品人妻久久久久久| 成人三级黄色视频| 欧美成人一区二区免费高清观看| 十八禁国产超污无遮挡网站| 亚洲精品在线观看二区| 99久久精品国产国产毛片| 欧美潮喷喷水| 亚洲精品日韩av片在线观看| 久久亚洲真实| 国产亚洲精品av在线| eeuss影院久久| 色精品久久人妻99蜜桃| 午夜免费男女啪啪视频观看 | 听说在线观看完整版免费高清| 国产精华一区二区三区| 精品久久久久久久人妻蜜臀av| 亚洲精品一区av在线观看| 搡老妇女老女人老熟妇| 亚洲av电影不卡..在线观看| 日韩中文字幕欧美一区二区| 国产精品女同一区二区软件 | 亚州av有码| 欧美色视频一区免费| 日韩欧美在线乱码| 男女下面进入的视频免费午夜| 久久国内精品自在自线图片| 欧美日韩亚洲国产一区二区在线观看| 美女xxoo啪啪120秒动态图| 97人妻精品一区二区三区麻豆| 免费在线观看影片大全网站| 久久欧美精品欧美久久欧美| 日韩欧美一区二区三区在线观看| 九九热线精品视视频播放| 午夜老司机福利剧场| 女的被弄到高潮叫床怎么办 | 午夜免费成人在线视频| 日本免费a在线| 国内久久婷婷六月综合欲色啪| 永久网站在线| 国产高清有码在线观看视频| 少妇人妻精品综合一区二区 | 亚洲精华国产精华液的使用体验 | 国产 一区 欧美 日韩| 中文字幕熟女人妻在线| 国产一区二区三区在线臀色熟女| 国产精品三级大全| 国产精品爽爽va在线观看网站| 婷婷精品国产亚洲av| 亚洲va日本ⅴa欧美va伊人久久| 国产综合懂色| 亚洲av中文字字幕乱码综合| 波野结衣二区三区在线| 国产精品久久久久久亚洲av鲁大| 色在线成人网| 国产精品免费一区二区三区在线| 老熟妇乱子伦视频在线观看| av在线老鸭窝| netflix在线观看网站| 日韩中文字幕欧美一区二区| 欧美日韩精品成人综合77777| 色综合站精品国产| 日韩高清综合在线| 欧美日韩亚洲国产一区二区在线观看| 国产亚洲精品av在线| av.在线天堂| 成人国产麻豆网| 俄罗斯特黄特色一大片| 欧美黑人巨大hd| 国产真实乱freesex| 欧美国产日韩亚洲一区| 亚洲性久久影院| 在线观看舔阴道视频| 一进一出抽搐gif免费好疼| 窝窝影院91人妻| av福利片在线观看| 91午夜精品亚洲一区二区三区 | 亚洲性夜色夜夜综合| 日韩 亚洲 欧美在线| 国产在视频线在精品| 午夜福利在线观看免费完整高清在 | 不卡一级毛片| 99久久中文字幕三级久久日本| 亚洲五月天丁香| 欧美成人免费av一区二区三区| 搡老熟女国产l中国老女人| 在线观看av片永久免费下载| 波多野结衣高清作品| 久久九九热精品免费| 一本精品99久久精品77| 免费看光身美女| 国产一区二区三区av在线 | 日日摸夜夜添夜夜添小说| 人妻夜夜爽99麻豆av| 免费不卡的大黄色大毛片视频在线观看 | 国产高清有码在线观看视频| 亚洲一区高清亚洲精品| 亚洲美女搞黄在线观看 | 在现免费观看毛片| 午夜激情欧美在线| 日日夜夜操网爽| 欧美日韩综合久久久久久 | 亚洲乱码一区二区免费版| 国产蜜桃级精品一区二区三区| 国产精品久久久久久av不卡| 桃色一区二区三区在线观看| 亚洲av中文字字幕乱码综合| 免费在线观看影片大全网站| 精品午夜福利在线看| 免费在线观看成人毛片| 日韩在线高清观看一区二区三区 | 日本黄色视频三级网站网址| 88av欧美| 亚洲国产色片| 真人做人爱边吃奶动态| 在线观看av片永久免费下载| 九九热线精品视视频播放| 国产亚洲av嫩草精品影院| 老女人水多毛片| 夜夜看夜夜爽夜夜摸| 精品乱码久久久久久99久播| 男插女下体视频免费在线播放| 亚洲图色成人| 变态另类成人亚洲欧美熟女| 91麻豆精品激情在线观看国产| 嫁个100分男人电影在线观看| 精品久久久久久久久久免费视频| 国产极品精品免费视频能看的| 亚洲av成人av| 乱系列少妇在线播放| 中文字幕av在线有码专区| 22中文网久久字幕| 我要看日韩黄色一级片| 少妇人妻精品综合一区二区 | 欧美精品啪啪一区二区三区| 大型黄色视频在线免费观看| 国产精品野战在线观看| 91在线观看av| 国产极品精品免费视频能看的| 亚洲美女黄片视频| 精品人妻偷拍中文字幕| 欧美日韩中文字幕国产精品一区二区三区| 国产精品美女特级片免费视频播放器| 真人做人爱边吃奶动态| av天堂在线播放| 美女xxoo啪啪120秒动态图| 91麻豆av在线| 一个人看的www免费观看视频| 日本 av在线| 欧美xxxx黑人xx丫x性爽| 极品教师在线视频| 国产黄片美女视频| 国产免费av片在线观看野外av| 国模一区二区三区四区视频| 久久久久久久亚洲中文字幕| 伊人久久精品亚洲午夜| 欧美另类亚洲清纯唯美| 欧美高清成人免费视频www| 欧美+亚洲+日韩+国产| 在线观看美女被高潮喷水网站| 日韩 亚洲 欧美在线| 99九九线精品视频在线观看视频| 日韩亚洲欧美综合| 91麻豆av在线| 欧美丝袜亚洲另类 | 国产精品亚洲美女久久久| 日日啪夜夜撸| 日本爱情动作片www.在线观看 | 亚洲国产精品成人综合色| 色av中文字幕| 久久香蕉精品热| 有码 亚洲区| or卡值多少钱| 美女cb高潮喷水在线观看| 大又大粗又爽又黄少妇毛片口| 久久精品人妻少妇| 午夜亚洲福利在线播放| 欧洲精品卡2卡3卡4卡5卡区| 亚洲av五月六月丁香网| 亚洲av.av天堂| 国产亚洲av嫩草精品影院| 中文字幕av在线有码专区| 国产精品无大码| 91久久精品国产一区二区三区| 狂野欧美白嫩少妇大欣赏| 亚洲性久久影院| 久久99热这里只有精品18| 一本一本综合久久| 看黄色毛片网站| 一级a爱片免费观看的视频| 99久久无色码亚洲精品果冻| 老司机福利观看| 全区人妻精品视频| 99在线人妻在线中文字幕| 亚洲国产色片| 国产精品亚洲美女久久久| 精品久久久久久,| 国模一区二区三区四区视频| 午夜福利成人在线免费观看| 欧美日韩综合久久久久久 | 亚洲18禁久久av| 国产中年淑女户外野战色| 亚洲国产精品合色在线| 国产私拍福利视频在线观看| 成人av在线播放网站| 国产伦一二天堂av在线观看| 久久精品夜夜夜夜夜久久蜜豆| 亚洲aⅴ乱码一区二区在线播放| 欧美丝袜亚洲另类 | 亚洲电影在线观看av| 九色成人免费人妻av| 一级av片app| 在线天堂最新版资源| 久久精品综合一区二区三区| 亚洲av不卡在线观看| 欧美日韩综合久久久久久 | 国产三级中文精品| 久久久精品欧美日韩精品| 欧美成人免费av一区二区三区| 亚洲人成网站在线播放欧美日韩| 亚洲av.av天堂| 在线免费十八禁| 欧美人与善性xxx| 国内揄拍国产精品人妻在线| 一夜夜www| 午夜久久久久精精品| 91久久精品电影网| 91麻豆精品激情在线观看国产| 精品人妻熟女av久视频| 日韩欧美在线乱码| 日本 av在线| 国产精品美女特级片免费视频播放器| 亚洲最大成人中文| 国产精品免费一区二区三区在线| 亚洲欧美清纯卡通| 欧美日韩亚洲国产一区二区在线观看| 国产亚洲精品av在线| 久久人人精品亚洲av| 18禁裸乳无遮挡免费网站照片| 看黄色毛片网站| 午夜a级毛片| 可以在线观看的亚洲视频| 日本熟妇午夜| 亚洲熟妇熟女久久| 人妻少妇偷人精品九色| 欧美激情久久久久久爽电影| 免费高清视频大片| 成人性生交大片免费视频hd| 成人欧美大片| 又黄又爽又免费观看的视频| 亚洲精品粉嫩美女一区| 国产精品美女特级片免费视频播放器| 91在线观看av| 热99在线观看视频| 一进一出抽搐动态| av福利片在线观看| 嫩草影院精品99| 国产精品永久免费网站| 欧美又色又爽又黄视频| 嫩草影院入口| 成人精品一区二区免费| 亚洲成人精品中文字幕电影| 国产男人的电影天堂91| 精品人妻熟女av久视频| 国产精品一区二区三区四区免费观看 | 美女cb高潮喷水在线观看| 欧美日韩国产亚洲二区| 亚洲性夜色夜夜综合| netflix在线观看网站| 国产私拍福利视频在线观看| 国产大屁股一区二区在线视频| 亚洲成人精品中文字幕电影| 亚洲人成网站在线播| 男女下面进入的视频免费午夜| 国产 一区 欧美 日韩| 美女高潮喷水抽搐中文字幕| 久99久视频精品免费| 99久久九九国产精品国产免费| 亚洲欧美日韩卡通动漫| 日本成人三级电影网站| 最后的刺客免费高清国语| 日本成人三级电影网站| 亚洲五月天丁香| 国产伦在线观看视频一区| av中文乱码字幕在线| 午夜免费成人在线视频| 露出奶头的视频| 亚洲久久久久久中文字幕| 精品人妻偷拍中文字幕| 欧美xxxx黑人xx丫x性爽| 亚洲无线观看免费| 久久久久久久亚洲中文字幕| 亚洲乱码一区二区免费版| 亚洲国产精品成人综合色| 成人美女网站在线观看视频| 亚洲成人精品中文字幕电影| 免费大片18禁| 久久国产精品人妻蜜桃| 亚洲av中文av极速乱 | 在线观看美女被高潮喷水网站| 午夜日韩欧美国产| 91久久精品国产一区二区三区| 亚洲成av人片在线播放无| 偷拍熟女少妇极品色| 日韩国内少妇激情av| 亚洲七黄色美女视频| 国产毛片a区久久久久| 亚洲精品影视一区二区三区av| 国产毛片a区久久久久| 国模一区二区三区四区视频| 精品一区二区三区视频在线观看免费| 好男人在线观看高清免费视频| 成人毛片a级毛片在线播放| 麻豆国产av国片精品| 一进一出抽搐gif免费好疼| 国产中年淑女户外野战色| 国产精品无大码| 免费看av在线观看网站| 美女免费视频网站| 十八禁国产超污无遮挡网站| 99久久久亚洲精品蜜臀av| 成人性生交大片免费视频hd| 久久精品国产亚洲av涩爱 | 亚洲国产高清在线一区二区三| 国产欧美日韩一区二区精品| 亚洲 国产 在线| 国产亚洲欧美98| av国产免费在线观看| 精品久久久久久成人av| 国产91精品成人一区二区三区| 99热这里只有精品一区| 欧美精品国产亚洲| 18+在线观看网站| 日本a在线网址| 在线观看免费视频日本深夜| 桃红色精品国产亚洲av| 一进一出抽搐gif免费好疼| 麻豆精品久久久久久蜜桃| 日日啪夜夜撸| 中文字幕免费在线视频6| 中文字幕av成人在线电影| 午夜免费激情av| 国产午夜福利久久久久久| 午夜精品久久久久久毛片777| av天堂在线播放| 97热精品久久久久久| 国产极品精品免费视频能看的| 国产淫片久久久久久久久| 日韩精品青青久久久久久| 成人性生交大片免费视频hd| 亚洲自拍偷在线| 最近最新中文字幕大全电影3| 两性午夜刺激爽爽歪歪视频在线观看| 综合色av麻豆| 久久久久国产精品人妻aⅴ院| 久9热在线精品视频| 精品久久久久久久久亚洲 | 亚洲成人久久性| 三级毛片av免费| 国产蜜桃级精品一区二区三区| 精品久久久久久久久久免费视频| 在线观看美女被高潮喷水网站| 老熟妇乱子伦视频在线观看| 国产美女午夜福利| 国产一区二区三区av在线 | 网址你懂的国产日韩在线| 最近最新免费中文字幕在线| 久久久久国内视频| 国产综合懂色| 男女下面进入的视频免费午夜| 午夜爱爱视频在线播放| 九九在线视频观看精品| 性色avwww在线观看| 白带黄色成豆腐渣| 免费看美女性在线毛片视频| 日韩一区二区视频免费看| 久久中文看片网| 美女 人体艺术 gogo| 日本成人三级电影网站| 高清毛片免费观看视频网站| 看十八女毛片水多多多| 国产午夜精品久久久久久一区二区三区 | 十八禁网站免费在线| 真人做人爱边吃奶动态| 欧美一区二区国产精品久久精品| 国产av一区在线观看免费| 久久久久久伊人网av| 直男gayav资源| 久久精品久久久久久噜噜老黄 | 国产 一区 欧美 日韩| 久久99热这里只有精品18| 99热这里只有是精品50| 国产av不卡久久| 真实男女啪啪啪动态图| 国产男人的电影天堂91| 窝窝影院91人妻| 精品久久国产蜜桃| 免费一级毛片在线播放高清视频| 亚洲自偷自拍三级| 久久精品国产鲁丝片午夜精品 | 毛片女人毛片| 国产精品久久久久久久电影| 最近最新免费中文字幕在线| 看十八女毛片水多多多| 99久久无色码亚洲精品果冻| 不卡视频在线观看欧美| 日本免费a在线| 国产精品一区二区三区四区免费观看 | 淫秽高清视频在线观看| 亚洲人成网站在线播放欧美日韩| 三级毛片av免费| 日韩欧美三级三区| 亚洲无线观看免费| av女优亚洲男人天堂| 亚洲成人中文字幕在线播放| 真人做人爱边吃奶动态| 欧美zozozo另类| 亚洲 国产 在线| 内地一区二区视频在线| av在线天堂中文字幕| 国产精品一区二区性色av| 欧美丝袜亚洲另类 | 日韩大尺度精品在线看网址| 国产高清有码在线观看视频| 亚洲第一区二区三区不卡| 国产亚洲欧美98| 18禁黄网站禁片免费观看直播| 亚洲男人的天堂狠狠| 啦啦啦观看免费观看视频高清| 久久精品夜夜夜夜夜久久蜜豆| 人妻久久中文字幕网| 日韩强制内射视频| 精品久久久久久久久久久久久| 亚洲美女视频黄频| 直男gayav资源| 国产蜜桃级精品一区二区三区|