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

    一種循環(huán)迭代的寬帶MIMO雷達(dá)正交稀疏頻譜波形設(shè)計(jì)方法

    2017-10-14 02:55:40李玉翔任修坤鄭娜娥
    電子與信息學(xué)報(bào) 2017年4期
    關(guān)鍵詞:阻帶陷波旁瓣

    李玉翔 任修坤 孫 揚(yáng) 鄭娜娥

    ?

    一種循環(huán)迭代的寬帶MIMO雷達(dá)正交稀疏頻譜波形設(shè)計(jì)方法

    李玉翔*任修坤 孫 揚(yáng) 鄭娜娥

    (信息工程大學(xué)導(dǎo)航與空天目標(biāo)工程學(xué)院 鄭州 450001)

    甚高頻(VHF)和超高頻(UHF)寬帶MIMO雷達(dá)在實(shí)現(xiàn)目標(biāo)高分辨與成像方面表現(xiàn)出巨大潛力,但也面臨著工作頻段擁塞和電磁信號(hào)干擾問(wèn)題。針對(duì)這一問(wèn)題,該文提出一種基于循環(huán)迭代的寬帶MIMO雷達(dá)正交稀疏頻譜波形設(shè)計(jì)方法。首先,將期望頻譜矩陣作為待優(yōu)化的輔助變量,綜合考慮最小化發(fā)射波形頻譜與期望頻譜之間的均方誤差以及積分旁瓣電平為目標(biāo)函數(shù),以發(fā)射波形恒模和期望頻譜幅度滿足上、下界約束為條件建立優(yōu)化模型;然后,在循環(huán)迭代的算法框架下,為提高運(yùn)算效率,利用快速傅里葉變換和譜分解技術(shù)實(shí)現(xiàn)模型的求解。仿真結(jié)果表明,優(yōu)化波形較現(xiàn)有方法具有更好的稀疏頻譜特性,同時(shí)還具有低自相關(guān)旁瓣和互相關(guān)性。

    寬帶MIMO雷達(dá);稀疏頻譜;正交波形;循環(huán)迭代

    1 引言

    UHF和VHF波段寬帶MIMO成像雷達(dá)[1]對(duì)樹(shù)林和地表有很強(qiáng)的穿透能力,能夠探測(cè)隱蔽目標(biāo),具有分辨率高、攜有的目標(biāo)信息量大的特點(diǎn),符合現(xiàn)代雷達(dá)高分辨與成像等技術(shù)的發(fā)展趨勢(shì),受到越來(lái)越廣泛的關(guān)注。然而,隨著雷達(dá)和通信技術(shù)的迅速發(fā)展,戰(zhàn)場(chǎng)的電磁環(huán)境變得錯(cuò)綜復(fù)雜,工作于VHF以及UHF擁塞頻段的寬帶MIMO雷達(dá)不僅受到通訊、導(dǎo)航等系統(tǒng)保留頻帶的限制[4],也面臨被其它輻射源信號(hào)干擾的問(wèn)題。為了避開(kāi)電磁干擾頻段,如何設(shè)計(jì)能量在頻域上稀疏分布的信號(hào),即稀疏頻譜波形[7](Sparse Frequency Waveform, SFW),同時(shí)滿足系統(tǒng)的寬頻帶要求,成為工作于擁塞頻段的寬帶MIMO雷達(dá)亟待解決的問(wèn)題。所謂稀疏信號(hào),是指能量在頻域上稀疏分布的信號(hào),使其避開(kāi)干擾較大的頻帶,盡量將信號(hào)能量分布在干擾較小的頻帶上。

    針對(duì)傳統(tǒng)雷達(dá)的頻譜干擾抑制問(wèn)題,Lindenfeld[8]采用“零陷”思想來(lái)設(shè)計(jì)稀疏頻譜波形,在干擾頻帶形成阻帶(陷波),然后根據(jù)優(yōu)化后的波形設(shè)計(jì)時(shí)域失配濾波器來(lái)降低脈沖壓縮后的距離旁瓣。但是該方法在接收濾波器約束主瓣為常數(shù)時(shí)存在失配損失。He等人[9]基于迭代的思想提出了阻帶循環(huán)算法SCAN(Stopband Cyclic Algorithm New),實(shí)現(xiàn)了SFW波形設(shè)計(jì)。Aubry等人[10]在設(shè)計(jì)稀疏頻譜波形時(shí)考慮最大化接收端信干噪比,建立非凸優(yōu)化模型,然后采用半正定松弛技術(shù)進(jìn)行求解。但是上述方法只能用于單發(fā)射的稀疏頻譜波形設(shè)計(jì),無(wú)法應(yīng)用到MIMO雷達(dá)系統(tǒng)中。在文獻(xiàn)[9]的基礎(chǔ)上,文獻(xiàn)[11]在波形設(shè)計(jì)時(shí)綜合考慮功率譜密度(Power Spectrum Density, PSD)和自相關(guān)函數(shù)性能建立優(yōu)化模型,然后采用共軛梯度法進(jìn)行求解。該方法的優(yōu)化結(jié)果依賴(lài)算法初始值,且求解共軛梯度的計(jì)算復(fù)雜度也較高。文獻(xiàn)[12]在恒模約束條件下考慮了組網(wǎng)雷達(dá)系統(tǒng)的頻譜限制問(wèn)題,利用交替迭代[13]的思想設(shè)計(jì)了多站發(fā)射稀疏頻譜波形。該方法的缺點(diǎn)是優(yōu)化波形的PSD在通帶和陷波內(nèi)都存在較大的波動(dòng)。

    針對(duì)上述問(wèn)題,本文摒棄傳統(tǒng)SFW設(shè)計(jì)方法將期望功率譜密度作為固定變量(通常根據(jù)頻譜分配情況設(shè)定)的思路,提出了一種基于循環(huán)迭代的寬帶MIMO雷達(dá)正交稀疏頻譜波形設(shè)計(jì)新方法。首先,將期望頻譜作為待優(yōu)化的輔助變量之一,綜合考慮最小化發(fā)射波形頻譜與期望頻譜之間的均方誤差以及積分旁瓣電平為目標(biāo)函數(shù),以各陣元發(fā)射波形恒模以及期望頻譜幅度上、下界約束為條件建立優(yōu)化模型,從而使MIMO雷達(dá)各陣元的發(fā)射波形具有稀疏頻譜特性,同時(shí)滿足低自相關(guān)旁瓣和低互相關(guān);然后在循環(huán)迭代的算法框架下,利用快速傅里葉變換(Fast Fourier Transform, FFT)和譜分解技術(shù)對(duì)優(yōu)化模型進(jìn)行求解。仿真結(jié)果驗(yàn)證了本文方法的有效性。

    2 信號(hào)模型

    3 正交稀疏頻譜波形設(shè)計(jì)

    3.1目標(biāo)函數(shù)的建立

    傳統(tǒng)頻譜逼近方法優(yōu)化得到波形的PSD存在“震蕩”現(xiàn)象,即在期望PSD附近存在較大波動(dòng)。因此,本文考慮在設(shè)計(jì)稀疏頻譜時(shí)不將期望頻譜矩陣作為固定變量,而是將期望頻譜作為待優(yōu)化的輔助變量,通過(guò)設(shè)置期望頻譜幅度的上界和下界,在求解時(shí)使波形頻譜滿足該約束條件。假設(shè)期望頻譜幅度的上、下界向量分別為,以最小均方誤差準(zhǔn)則使MIMO雷達(dá)在恒模條件下逼近期望頻譜矩陣,可用式(4)的優(yōu)化問(wèn)題來(lái)表示。

    此外,為了在接收端獲得較低的距離旁瓣,MIMO雷達(dá)通常需要發(fā)射正交信號(hào),即發(fā)射波形具有低自相關(guān)和互相關(guān)旁瓣。發(fā)射信號(hào)和之間的相關(guān)函數(shù)定義為

    (5)

    (7)

    則式(6)的目標(biāo)函數(shù)可轉(zhuǎn)化為

    (8)

    為了使MIMO雷達(dá)發(fā)射波形同時(shí)具有稀疏頻譜特性和正交性能,本文建立聯(lián)合優(yōu)化模型為

    3.2 循環(huán)迭代求解

    (12)

    (13)

    (15)

    根據(jù)矩陣跡的性質(zhì),進(jìn)一步推導(dǎo)可以得到

    綜上所述,本文基于循環(huán)迭代的正交稀疏頻譜波形設(shè)計(jì)算法的具體流程如表1所示。

    表1正交稀疏頻譜波形設(shè)計(jì)算法流程

    步驟1 設(shè)定期望頻譜幅度的上界和下界,初始化迭代次數(shù),設(shè)置最大迭代次數(shù)為和權(quán)重因子。采用長(zhǎng)度為的QPSK信號(hào)作為初始信號(hào)記為,以隨機(jī)產(chǎn)生的方式初始化輔助變量; 步驟2 根據(jù)波形矩陣和輔助變量,利用式(13)計(jì)算; 步驟3 根據(jù)波形矩陣和期望頻譜矩陣,由式(15)得到; 步驟4 根據(jù)期望頻譜矩陣和輔助變量,利用式(18)求解; 步驟5 判斷臨時(shí)變量( )是否滿足終止條件或者迭代次數(shù)溢出。若條件成立終止迭代,否則,執(zhí)行并跳轉(zhuǎn)至步驟2。

    4 實(shí)驗(yàn)仿真

    實(shí)驗(yàn)1 為了驗(yàn)證優(yōu)化波形的干擾抑制能力和正交特性,假設(shè)超寬帶合成孔徑MIMO雷達(dá)工作于150~550 MHz。當(dāng)前頻段內(nèi)無(wú)干擾的可用信道分別為[194,206][238,274][294,314][370,470][490,530] MHz,稀疏頻譜波形需要在其余頻段內(nèi)形成阻帶。假設(shè)波形時(shí)寬為,以MHz時(shí)鐘采樣可得發(fā)射波形碼長(zhǎng), FFT點(diǎn)數(shù)取,權(quán)重因子取。期望頻譜下界設(shè)置為充分小,取值為dB。在無(wú)干擾頻帶內(nèi)期望頻譜上界設(shè)置為充分大,取值為dB,干擾頻帶設(shè)置為dB,偏置量取dB。圖1給出了文獻(xiàn)[11],文獻(xiàn)[12]以及本文方法得到的各個(gè)陣元功率譜密度曲線。圖2以MIMO雷達(dá)第1個(gè)陣元為參考,給出了3種方法得到波形在持續(xù)時(shí)間為的功率分布情況??梢钥闯觯墨I(xiàn)[11]波形的PSD平均陷波深度約為dB,且在通帶和陷波內(nèi)存在較大的波動(dòng),文獻(xiàn)[12]波形的PSD在通帶和陷波內(nèi)的波動(dòng)幅度較文獻(xiàn)[11]小,性能優(yōu)于文獻(xiàn)[11]。本文方法的PSD被嚴(yán)格控制在dB之下,且功率譜密度較為平坦,波形陷波能力優(yōu)于文獻(xiàn)[11],文獻(xiàn)[12],因此具有更好的電磁干擾抑制能力。

    圖3給出了文獻(xiàn)[11],文獻(xiàn)[12]以及本文所優(yōu)化波形的ACF和CCF曲線,表2對(duì)應(yīng)給出了各陣元的自相關(guān)峰值旁瓣和互相關(guān)峰值旁瓣。文獻(xiàn)[11]得到陣元1~3的平均為dB,平均為dB,文獻(xiàn)[12]得到陣元1~3的平均為dB,平均為dB,本文方法陣元1~3波形的平均為dB,平均為dB。實(shí)驗(yàn)結(jié)果可知所提方法得到波形的自相關(guān)旁瓣電平更低,同時(shí)具有更低的互相關(guān),可以保證在雷達(dá)接收端獲得具有更好的脈沖壓縮性能。

    實(shí)驗(yàn)2 為了分析陷波個(gè)數(shù)對(duì)算法性能的影響,假設(shè)150~550 MHz工作帶寬內(nèi)的干擾信道為[370, 450] MHz,稀疏頻譜波形只需在該頻段形成單阻帶。圖4給出了本文優(yōu)化波形的功率譜密度曲線,可以看出PSD形成更深的陷波,平均陷波深度約為dB。圖5給出了本文算法在單阻帶條件的ACF和CCF曲線,此時(shí)MIMO雷達(dá)3個(gè)陣元的分別為dB,dB,dB,分別為dB,dB,dB,波形的PSD和ACF性能均優(yōu)于實(shí)驗(yàn)1。實(shí)驗(yàn)結(jié)果可知,較少的阻帶個(gè)數(shù)為MIMO雷達(dá)波形設(shè)計(jì)提供了更多自由度,設(shè)計(jì)波形的干擾抑制能力和正交性會(huì)在一定程度上得到改善。

    圖1 各陣元的功率譜密度比較曲線

    圖2 陣元1優(yōu)化波形在內(nèi)的功率分布情況比較

    圖3 優(yōu)化波形的ACF, CCF比較曲線

    表2自/互相關(guān)峰值旁瓣對(duì)比結(jié)果

    (dB)(dB) 陣元1陣元2陣元3陣元1陣元2陣元3 文獻(xiàn)[11]-16.92-17.26-17.14-17.73-18.12-18.05 文獻(xiàn)[12]-17.62-17.87-17.79-18.56-18.89-18.47 本文算法-19.83-20.41-20.33-21.42-21.30-20.96

    圖5 優(yōu)化波形的ACF, CCF比較曲線(單阻帶)

    圖6 權(quán)重因子對(duì)本文算法的影響

    5 結(jié)論

    本文提出了一種基于循環(huán)迭代的寬帶MIMO雷達(dá)正交稀疏頻譜波形設(shè)計(jì)方法,能夠有效解決寬帶MIMO雷達(dá)面臨的工作頻段擁塞和電磁干擾問(wèn)題。在設(shè)計(jì)過(guò)程中,利用了循環(huán)迭代的算法框架,通過(guò)FFT實(shí)現(xiàn)主要運(yùn)算,計(jì)算效率高。仿真結(jié)果表明,優(yōu)化后的發(fā)射波形功率譜密度曲線較為平坦,頻譜陷波性能更好;同時(shí)波形的自/互相關(guān)峰值旁瓣電平更低,能保證更好的脈沖壓縮性能。阻帶數(shù)目的減少為波形設(shè)計(jì)提供更多的自由度,能改善頻譜陷波性能和波形的正交性。權(quán)重因子可以根據(jù)實(shí)際電磁環(huán)境進(jìn)行合理設(shè)置,提高了波形設(shè)計(jì)方法的靈活性和適應(yīng)性。

    [1] 李春升, 王偉杰, 王鵬波, 等. 星載SAR技術(shù)的現(xiàn)狀與發(fā)展趨勢(shì)[J]. 電子與信息學(xué)報(bào), 2016, 38(1): 229-240. doi: 10.11999/ JEIT151116.

    LI Chunsheng, WANG Weijie, WANG Pengbo,. Current situation and development trends of spaceborne SAR technology[J].&, 2016, 38(1): 229-240. doi: 10.11999/JEIT151116.

    [2] ZHANG Yongshun, JIA Xin, and ZHU Weigang. Suppression of LFM interference in direct sequence spread spectrum communications based on compressive sensing[C]. IEEE International Conference on Signal Processing, Communications and Computing(ICSPCC), Ningbo, China, 2015: 1-5. doi: 10.1109/ICSPCC.2015.7338857.

    [3] 胡恒, 賀亞鵬, 莊珊娜, 等. 高頻地波雷達(dá)稀疏頻率波形優(yōu)化設(shè)計(jì)[J]. 電子與信息學(xué)報(bào), 2012, 34(6): 1291-1296. doi: 10. 3724/SP.J.1146.2011.00666.

    HU Heng, HE Yapeng, ZHUANG Shanna,. Sparse frequency waveform design for high frequency surface wave radar[J].&, 2012, 34(6): 1291-1296. doi: 10.3724/SP.J.1146.2011.00666.

    [4] WANG Guohua, MAI Chaoyun, SUN Jinping,. Sparse frequency waveform analysis and design based on ambiguity function theory[J].,&, 2016, 10(4): 707-717. doi: 10.1049/IET-RSN.2015.0270.

    [5] 申?yáng)|, 張林讓, 周宇, 等. 頻帶陷波優(yōu)化的低距離旁瓣波形設(shè)計(jì)方法[J]. 北京理工大學(xué)學(xué)報(bào), 2015, 35(3): 321-325. doi: 10.15918/J.TBIT1001-0645.2015.03029.

    SHEN Dong, ZHANG Linrang, ZHOU Yu,. Waveform design with stopband and low range sidelobes[J]., 2015, 35(3): 321-325. doi: 10.15918/J.TBIT1001-0645.2015.03029.

    [6] 吳灝, 范紅旗, 付強(qiáng). 低距離旁瓣雷達(dá)系統(tǒng)設(shè)計(jì)方法綜述[J]. 航空兵器, 2015, 53(4): 3-8. doi: 10.3969/J.ISSN.1673-5048. 2015.04.001.

    WU Hao, FAN Hongqi, and FU Qiang. Overview of design methods for radar system with low range sidelobes[J]., 2015, 53(4): 3-8. doi: 10.3969/J.ISSN.1673-5048. 2015.04.001.

    [7] ZHAO Dehua and WEI Yinsheng. Coherent process and optimal weighting for sparse frequency agility waveform[C]. IEEE Radar Conference, Arlington, 2015: 334-338. doi: 10.1109/RADAR.2015.7131020.

    [8] LINDENFELD M J. Sparse frequency transmit and receive waveform design[J]., 2004, 40(3): 851-861. doi: 10.1109/TAES. 2004.1337459.

    [9] HE Hao, STOICA P, and LI Jian. Waveform design with stopband and correlation constraints for cognitive radar[C]. International Workshop on Cognitive Information Processing (CIP), Elba, Italy, 2010: 344-349. doi: 10.1109/CIP.2010. 5604089.

    [10] AUBRY A, MAIO A De, PIEZZO M,. Radar waveform design in a spectrally crowded environment via nonconvex quadratic optimization[J]., 2014, 50(2): 1138-1152. doi: 10.1109/ TAES.2014.120731.

    [11] WANG G H and LU Y L. Designing single/multiple sparse frequency waveforms with sidelobe constraint[J].&, 2011, 5(1): 32-38. doi: 10.1049/IET- RSN.2009.0255.

    [12] 周宇, 張林讓, 趙珊珊. 組網(wǎng)雷達(dá)低自相關(guān)旁瓣和互相關(guān)干擾的稀疏頻譜波形設(shè)計(jì)方法[J]. 電子與信息學(xué)報(bào), 2014, 36(6): 1394-1399. doi: 10.3724/SP.J.1146.2013.00702.

    ZHOU Yu, ZHANG Linrang, and ZHAO Shanshan. Sparse frequency waveforms design with low correlation sidelobes for netted radar[J].&, 2014, 36(6): 1394-1399. doi: 10.3724/SP.J.1146. 2013.00702.

    [13] TROPP J A, DHILLON I S, HEATH R W,. Design structured tight frames via an alternating projection method[J]., 2005, 51(1): 188-209. doi: 10.1109/TIT.2004.839492.

    [14] HE Hao, STOICA P, and LI Jian. Designing unimodular sequence sets with good correlations-Including an application to MIMO radar[J]., 2009, 57(11): 4391-4405. doi: 10.1109/TSP.2009.2025108.

    [15] STOICA P, HE Hao, and LI Jian. New algorithms for designing unimodular sequences with good correlation properties[J]., 2009, 57(4): 1415-1425. doi: 10.1109/TSP.2009.2012562.

    Cyclic Iterative Method for Wideband MIMO Radar Orthogonal Sparse Frequency Waveform Design

    LI Yuxiang REN Xiukun SUN Yang ZHENG Nae

    (,,450001,)

    Wideband MIMO radar exhibits great potential in achieving the goals of high resolution imaging, but it also suffers from the electromagnetic signal congestion and interference problems, especially for radar that works in Very High Frequency (VHF) and Ultra High Frequency (UHF) band. To solve this problem, a cyclic iterative method for designing orthogonal sparse frequency waveforms is proposed. Firstly, the desired spectrum is used as an auxiliary variable, and a new objective function is constructed based on both the mean square error of the spectrum of transmitting waveform with the desired one and the integration side-lobe levels. The optimization model is established under the constraint that the waveform is constant envelope meanwhile the spectrum magnitude lies between the pre-established upper and lower bounds. Then, under the framework of cyclic iterative algorithm, fast Fourier transform and spectral decomposition techniques are used to improve computational efficiency. Simulation results show that the proposed method has good performance in designing orthogonal sparse frequency waveforms with low auto-correlation and cross-correlation side lobes.

    Wideband MIMO radar; Sparse frequency; Orthogonal waveform; Cyclic iterating

    TN958

    A

    1009-5896(2017)04-0953-07

    10.11999/JEIT160597

    2016-06-06;

    改回日期:2016-11-02;

    2016-12-02

    李玉翔 liyuxiangwork@163.com

    國(guó)家自然科學(xué)基金(41301481)

    The National Natural Science Foundation of China (41301481)

    李玉翔: 男,1987年生,博士生,研究方向?yàn)镸IMO雷達(dá)波形設(shè)計(jì)、陣列信號(hào)處理.

    任修坤: 男,1979年生,講師,研究方向?yàn)闊o(wú)線通信技術(shù)和空間信號(hào)處理.

    孫 揚(yáng): 男,1992年生,碩士生,研究方向?yàn)镸IMO雷達(dá)資源分配.

    鄭娜娥: 男,1984年生,講師,研究方向?yàn)镸IMO信號(hào)處理和無(wú)線資源分配.

    猜你喜歡
    阻帶陷波旁瓣
    一種新型可調(diào)雙阻帶濾波器設(shè)計(jì)
    基于圓柱陣通信系統(tǒng)的廣義旁瓣對(duì)消算法
    一種基于線性規(guī)劃的頻率編碼旁瓣抑制方法
    一種改進(jìn)的最大信雜比MTD濾波器設(shè)計(jì)算法
    基于加權(quán)積分旁瓣最小化的隨機(jī)多相碼設(shè)計(jì)
    基于四項(xiàng)最低旁瓣Nuttall窗的插值FFT諧波分析
    溫度對(duì)陷波網(wǎng)絡(luò)阻抗的影響
    遺傳算法在UWB陷波天線設(shè)計(jì)中的應(yīng)用
    一種基于互補(bǔ)環(huán)縫諧振器抑制SSN的新方法
    衛(wèi)星導(dǎo)航接收機(jī)基于IIR陷波器的單頻干擾抑制性能分析
    可以在线观看毛片的网站| 全区人妻精品视频| 亚洲在线观看片| 久久99精品国语久久久| 美女大奶头视频| 男女下面进入的视频免费午夜| 中文欧美无线码| 免费观看a级毛片全部| 国产成人免费观看mmmm| 天美传媒精品一区二区| 免费在线观看成人毛片| 国产综合懂色| 亚洲国产高清在线一区二区三| 亚洲美女视频黄频| 天美传媒精品一区二区| 亚洲av电影不卡..在线观看| 特级一级黄色大片| 国产 一区精品| 淫秽高清视频在线观看| 建设人人有责人人尽责人人享有的 | 精品午夜福利在线看| 18禁在线无遮挡免费观看视频| 国产精品蜜桃在线观看| www.av在线官网国产| 国产午夜精品久久久久久一区二区三区| 国产成人福利小说| 99热6这里只有精品| 中文字幕久久专区| 永久网站在线| 少妇的逼好多水| 亚洲精品国产av成人精品| 特大巨黑吊av在线直播| 国产色婷婷99| 麻豆av噜噜一区二区三区| 国产成人福利小说| 久久久久久久久大av| 高清日韩中文字幕在线| 国产一区二区亚洲精品在线观看| 国产亚洲最大av| 欧美zozozo另类| 亚洲av免费高清在线观看| 桃色一区二区三区在线观看| 久久精品国产鲁丝片午夜精品| 乱码一卡2卡4卡精品| 小蜜桃在线观看免费完整版高清| 天天一区二区日本电影三级| 色网站视频免费| 午夜福利网站1000一区二区三区| 91午夜精品亚洲一区二区三区| av在线天堂中文字幕| 一边摸一边抽搐一进一小说| 国产精品一二三区在线看| 人人妻人人看人人澡| 91午夜精品亚洲一区二区三区| 一边亲一边摸免费视频| 2021少妇久久久久久久久久久| 亚洲国产精品国产精品| 国产一级毛片七仙女欲春2| www.色视频.com| 毛片一级片免费看久久久久| 国产精品久久电影中文字幕| 看十八女毛片水多多多| 美女国产视频在线观看| 精品欧美国产一区二区三| 男女那种视频在线观看| 人妻夜夜爽99麻豆av| 长腿黑丝高跟| 人妻系列 视频| 欧美高清成人免费视频www| 波多野结衣高清无吗| 男女国产视频网站| 欧美丝袜亚洲另类| 国产免费男女视频| 久久婷婷人人爽人人干人人爱| 亚洲精品影视一区二区三区av| 高清毛片免费看| 亚洲av二区三区四区| 日韩精品有码人妻一区| 男插女下体视频免费在线播放| 中文在线观看免费www的网站| 精品国产一区二区三区久久久樱花 | 免费av毛片视频| 日韩av在线大香蕉| 韩国av在线不卡| 欧美人与善性xxx| 午夜亚洲福利在线播放| 99久久人妻综合| 国产伦理片在线播放av一区| 级片在线观看| av专区在线播放| 啦啦啦韩国在线观看视频| 国产伦精品一区二区三区四那| 中文字幕人妻熟人妻熟丝袜美| 国产亚洲91精品色在线| 有码 亚洲区| 免费大片18禁| 一级毛片久久久久久久久女| 国产亚洲精品久久久com| 国产精品国产三级国产专区5o | 麻豆久久精品国产亚洲av| 91aial.com中文字幕在线观看| 亚洲av男天堂| 欧美一区二区亚洲| 午夜免费激情av| 日本黄色片子视频| 亚洲在线观看片| 在线观看66精品国产| 深夜a级毛片| 亚洲熟妇中文字幕五十中出| 国产高清有码在线观看视频| 亚洲人与动物交配视频| 久久久久国产网址| 亚洲自偷自拍三级| 精品无人区乱码1区二区| 欧美人与善性xxx| 日韩精品有码人妻一区| 亚洲欧美中文字幕日韩二区| 日韩成人伦理影院| 菩萨蛮人人尽说江南好唐韦庄 | 少妇的逼水好多| 午夜福利在线在线| 亚洲aⅴ乱码一区二区在线播放| 久久精品国产亚洲av天美| 你懂的网址亚洲精品在线观看 | 亚洲,欧美,日韩| 精品久久久久久久久亚洲| av卡一久久| 日本一本二区三区精品| 中文欧美无线码| 一二三四中文在线观看免费高清| 久久久a久久爽久久v久久| 国产免费一级a男人的天堂| 国产亚洲午夜精品一区二区久久 | 最近最新中文字幕大全电影3| 国产91av在线免费观看| 亚洲精品乱码久久久v下载方式| 一级毛片aaaaaa免费看小| 三级男女做爰猛烈吃奶摸视频| 久久国产乱子免费精品| 九色成人免费人妻av| 午夜福利视频1000在线观看| 99久国产av精品| 成人一区二区视频在线观看| 亚洲经典国产精华液单| 亚洲精品日韩在线中文字幕| av免费在线看不卡| 日韩av在线大香蕉| 超碰97精品在线观看| 亚洲三级黄色毛片| 只有这里有精品99| 亚洲精品影视一区二区三区av| 十八禁国产超污无遮挡网站| 精品人妻偷拍中文字幕| 少妇裸体淫交视频免费看高清| 成年av动漫网址| 日本黄色视频三级网站网址| 18禁在线播放成人免费| 精品99又大又爽又粗少妇毛片| 夫妻性生交免费视频一级片| 亚洲精品国产成人久久av| 欧美成人免费av一区二区三区| 欧美最新免费一区二区三区| 男女国产视频网站| 精品久久久久久久久久久久久| 1000部很黄的大片| av视频在线观看入口| 亚洲天堂国产精品一区在线| 亚洲精品自拍成人| av国产久精品久网站免费入址| 中文字幕久久专区| 少妇人妻一区二区三区视频| 日本欧美国产在线视频| 伊人久久精品亚洲午夜| 99在线人妻在线中文字幕| 99热全是精品| 亚洲中文字幕日韩| 成人漫画全彩无遮挡| 内地一区二区视频在线| 午夜日本视频在线| 最后的刺客免费高清国语| 欧美潮喷喷水| 菩萨蛮人人尽说江南好唐韦庄 | 国产精品1区2区在线观看.| 国产精品一区二区三区四区久久| 欧美3d第一页| 人妻夜夜爽99麻豆av| 我要搜黄色片| 国产精品电影一区二区三区| 亚洲欧美成人综合另类久久久 | 日韩欧美精品v在线| 国模一区二区三区四区视频| 国产精品不卡视频一区二区| 久久婷婷人人爽人人干人人爱| 简卡轻食公司| av福利片在线观看| 18禁在线播放成人免费| 亚洲三级黄色毛片| 久久婷婷人人爽人人干人人爱| 国产片特级美女逼逼视频| 国产伦一二天堂av在线观看| 高清av免费在线| 色视频www国产| 国产高清不卡午夜福利| 欧美一区二区亚洲| 国产探花在线观看一区二区| 亚洲美女视频黄频| 国产精品一二三区在线看| 国产高潮美女av| 亚洲久久久久久中文字幕| 大又大粗又爽又黄少妇毛片口| 成人亚洲精品av一区二区| 有码 亚洲区| 男女边吃奶边做爰视频| 性色avwww在线观看| 亚洲精品乱码久久久v下载方式| 黄色配什么色好看| 美女大奶头视频| 免费观看性生交大片5| 成人二区视频| 亚洲性久久影院| 亚洲欧美精品专区久久| 精品久久久久久久久av| 日韩 亚洲 欧美在线| 国产在线男女| 久久精品人妻少妇| 国产一区二区亚洲精品在线观看| 51国产日韩欧美| kizo精华| 亚洲国产精品sss在线观看| 亚洲国产最新在线播放| 国产在视频线在精品| 日韩人妻高清精品专区| 2021少妇久久久久久久久久久| 搞女人的毛片| 国产亚洲最大av| 国产在视频线在精品| 亚洲乱码一区二区免费版| 2021少妇久久久久久久久久久| 内射极品少妇av片p| 国产精品一区二区三区四区久久| 直男gayav资源| 国产极品天堂在线| 99在线视频只有这里精品首页| 人妻制服诱惑在线中文字幕| 听说在线观看完整版免费高清| 乱码一卡2卡4卡精品| 国产精品无大码| 国产片特级美女逼逼视频| 99久久九九国产精品国产免费| 欧美性猛交黑人性爽| 久久精品国产鲁丝片午夜精品| 亚洲真实伦在线观看| 免费观看的影片在线观看| 久久鲁丝午夜福利片| 九色成人免费人妻av| av线在线观看网站| 五月伊人婷婷丁香| 欧美3d第一页| 国产精品国产三级专区第一集| 18禁裸乳无遮挡免费网站照片| 91精品一卡2卡3卡4卡| 色吧在线观看| 亚洲精品乱码久久久v下载方式| 日韩欧美 国产精品| 色播亚洲综合网| 久久久欧美国产精品| 国产精品,欧美在线| 久热久热在线精品观看| 午夜福利高清视频| 午夜精品在线福利| 日本免费在线观看一区| 97热精品久久久久久| 欧美人与善性xxx| 国产一区二区在线av高清观看| 尤物成人国产欧美一区二区三区| 久久精品国产自在天天线| 久久人人爽人人片av| 亚洲欧美精品自产自拍| 成人无遮挡网站| 久久精品国产亚洲av涩爱| 九色成人免费人妻av| 尤物成人国产欧美一区二区三区| 亚洲美女视频黄频| 婷婷色av中文字幕| 国产极品天堂在线| 免费黄色在线免费观看| 一区二区三区高清视频在线| 欧美色视频一区免费| 一级黄色大片毛片| 美女内射精品一级片tv| 国产精品国产高清国产av| 国产高清国产精品国产三级 | 深夜a级毛片| 亚洲人与动物交配视频| 免费观看在线日韩| 亚洲综合精品二区| av天堂中文字幕网| 美女内射精品一级片tv| 国语自产精品视频在线第100页| 欧美性猛交╳xxx乱大交人| 一级毛片aaaaaa免费看小| 免费电影在线观看免费观看| 精品久久久久久久久久久久久| 免费观看a级毛片全部| 中文字幕精品亚洲无线码一区| 日本免费在线观看一区| 一本久久精品| 你懂的网址亚洲精品在线观看 | 纵有疾风起免费观看全集完整版 | 少妇熟女欧美另类| 一区二区三区高清视频在线| 禁无遮挡网站| 国产午夜福利久久久久久| 国产成人a∨麻豆精品| 不卡视频在线观看欧美| 观看美女的网站| 国产又色又爽无遮挡免| 久久精品影院6| videos熟女内射| 在线观看一区二区三区| 成人亚洲精品av一区二区| 亚洲欧美日韩无卡精品| 欧美3d第一页| 久久国产乱子免费精品| 极品教师在线视频| av国产免费在线观看| 欧美丝袜亚洲另类| 日本与韩国留学比较| 少妇人妻精品综合一区二区| 亚洲精品日韩在线中文字幕| 少妇熟女欧美另类| 国产av在哪里看| 51国产日韩欧美| 国产高潮美女av| 一本一本综合久久| 1000部很黄的大片| 午夜爱爱视频在线播放| 22中文网久久字幕| 水蜜桃什么品种好| 亚洲激情五月婷婷啪啪| 国产精品美女特级片免费视频播放器| 日韩av不卡免费在线播放| 麻豆成人午夜福利视频| 22中文网久久字幕| 免费黄色在线免费观看| 国产极品天堂在线| 女人被狂操c到高潮| 波多野结衣高清无吗| 水蜜桃什么品种好| 午夜福利在线在线| 日本黄色片子视频| 日本色播在线视频| 精品99又大又爽又粗少妇毛片| 亚洲精品成人久久久久久| 亚洲一区高清亚洲精品| 亚洲成人中文字幕在线播放| 午夜精品国产一区二区电影 | 日韩国内少妇激情av| av在线天堂中文字幕| 97超碰精品成人国产| 久久久久久大精品| 中国美白少妇内射xxxbb| 狂野欧美激情性xxxx在线观看| 午夜视频国产福利| 噜噜噜噜噜久久久久久91| 波多野结衣巨乳人妻| 欧美一区二区国产精品久久精品| 小蜜桃在线观看免费完整版高清| 久久久精品94久久精品| 国产午夜精品一二区理论片| 男女国产视频网站| 亚洲性久久影院| 日本三级黄在线观看| 人妻夜夜爽99麻豆av| 国产亚洲最大av| 国产成人免费观看mmmm| 18禁在线播放成人免费| 中文字幕免费在线视频6| 久久精品国产亚洲网站| 亚洲人成网站在线播| 黄色一级大片看看| 中文精品一卡2卡3卡4更新| 特大巨黑吊av在线直播| 身体一侧抽搐| 日韩成人av中文字幕在线观看| 中文字幕人妻熟人妻熟丝袜美| 国产精品福利在线免费观看| 大香蕉97超碰在线| 欧美激情国产日韩精品一区| 国产精品一区二区三区四区久久| 老司机福利观看| 22中文网久久字幕| 97热精品久久久久久| 伦精品一区二区三区| 亚洲熟妇中文字幕五十中出| 成年av动漫网址| 亚洲精品亚洲一区二区| 国产亚洲精品久久久com| 国产白丝娇喘喷水9色精品| 成人特级av手机在线观看| 99热这里只有是精品50| 99久久中文字幕三级久久日本| 国产精品伦人一区二区| 亚洲精品亚洲一区二区| 超碰97精品在线观看| 少妇人妻一区二区三区视频| 亚洲精品乱码久久久v下载方式| 色综合亚洲欧美另类图片| 又爽又黄a免费视频| 日韩大片免费观看网站 | 黄色日韩在线| 插阴视频在线观看视频| 国产免费一级a男人的天堂| 久久久国产成人免费| 中文字幕av成人在线电影| 久久人人爽人人片av| 亚洲自拍偷在线| 如何舔出高潮| 亚洲综合色惰| 国产淫语在线视频| 最近视频中文字幕2019在线8| 尾随美女入室| 亚洲av熟女| 波野结衣二区三区在线| 国内精品宾馆在线| 久久欧美精品欧美久久欧美| 久久久国产成人精品二区| 日本色播在线视频| 又爽又黄a免费视频| av又黄又爽大尺度在线免费看 | 国产精品,欧美在线| 啦啦啦观看免费观看视频高清| 国产精品三级大全| 伦精品一区二区三区| 亚洲aⅴ乱码一区二区在线播放| 一级黄色大片毛片| av女优亚洲男人天堂| 久久久午夜欧美精品| 美女高潮的动态| 久久精品久久久久久噜噜老黄 | 国产亚洲最大av| 熟妇人妻久久中文字幕3abv| 性插视频无遮挡在线免费观看| 色网站视频免费| 国产精品嫩草影院av在线观看| 男女啪啪激烈高潮av片| 国产精品国产三级国产专区5o | 久久久久久久亚洲中文字幕| 国产亚洲午夜精品一区二区久久 | 国产精品电影一区二区三区| or卡值多少钱| 欧美+日韩+精品| 亚洲美女视频黄频| 欧美激情久久久久久爽电影| 国产精品三级大全| videos熟女内射| 久久99精品国语久久久| 最近2019中文字幕mv第一页| 色网站视频免费| 一夜夜www| 国产精品不卡视频一区二区| 一边亲一边摸免费视频| 十八禁国产超污无遮挡网站| 热99在线观看视频| 国产单亲对白刺激| 波多野结衣高清无吗| 精品一区二区三区视频在线| 插逼视频在线观看| 欧美激情在线99| 色5月婷婷丁香| 成年女人永久免费观看视频| 别揉我奶头 嗯啊视频| 日本午夜av视频| 成人鲁丝片一二三区免费| 久久久亚洲精品成人影院| 国产精华一区二区三区| 国产黄片美女视频| 国产高清有码在线观看视频| 亚洲性久久影院| 中文字幕熟女人妻在线| a级一级毛片免费在线观看| 国产精品一区二区三区四区久久| 高清在线视频一区二区三区 | 午夜亚洲福利在线播放| 三级男女做爰猛烈吃奶摸视频| 国产亚洲av嫩草精品影院| 亚洲av中文av极速乱| 亚洲国产成人一精品久久久| 3wmmmm亚洲av在线观看| 成年av动漫网址| 高清日韩中文字幕在线| 97超碰精品成人国产| 男女视频在线观看网站免费| ponron亚洲| 如何舔出高潮| 国产男人的电影天堂91| 亚洲在久久综合| 亚洲欧美清纯卡通| 日产精品乱码卡一卡2卡三| 精品99又大又爽又粗少妇毛片| 五月玫瑰六月丁香| 熟妇人妻久久中文字幕3abv| 禁无遮挡网站| 国产午夜精品久久久久久一区二区三区| 久久精品国产鲁丝片午夜精品| 麻豆国产97在线/欧美| 亚洲欧美日韩东京热| 亚洲最大成人手机在线| 亚洲欧美日韩无卡精品| 久久久久久久亚洲中文字幕| 国产精品女同一区二区软件| 欧美日韩综合久久久久久| 国产精品一区www在线观看| 婷婷色综合大香蕉| 亚洲三级黄色毛片| 免费观看性生交大片5| 永久网站在线| 嘟嘟电影网在线观看| 国产一区亚洲一区在线观看| 国产精品99久久久久久久久| 春色校园在线视频观看| 自拍偷自拍亚洲精品老妇| 天堂影院成人在线观看| 舔av片在线| 三级男女做爰猛烈吃奶摸视频| 狂野欧美白嫩少妇大欣赏| 国内揄拍国产精品人妻在线| 国产精品野战在线观看| 毛片一级片免费看久久久久| 日韩一本色道免费dvd| 亚洲美女搞黄在线观看| 成年女人永久免费观看视频| 国产黄片美女视频| 免费黄色在线免费观看| 国产精品久久久久久精品电影| 亚洲av中文字字幕乱码综合| 亚洲婷婷狠狠爱综合网| 2022亚洲国产成人精品| 乱码一卡2卡4卡精品| 欧美性猛交╳xxx乱大交人| 美女脱内裤让男人舔精品视频| 欧美人与善性xxx| 国产精品一区二区三区四区免费观看| 一个人免费在线观看电影| 久久久久免费精品人妻一区二区| 欧美一区二区亚洲| 内射极品少妇av片p| 亚洲国产欧美在线一区| 国产成人a∨麻豆精品| 国产精品乱码一区二三区的特点| av在线天堂中文字幕| 国产精品人妻久久久久久| 爱豆传媒免费全集在线观看| 国产在视频线在精品| 日韩大片免费观看网站 | 日韩,欧美,国产一区二区三区 | 亚洲欧美日韩东京热| 日韩av在线大香蕉| 小说图片视频综合网站| 亚洲精品国产av成人精品| 亚洲av日韩在线播放| 精品人妻熟女av久视频| 国产精品福利在线免费观看| 亚洲av不卡在线观看| 国内揄拍国产精品人妻在线| 在线播放国产精品三级| 国产成年人精品一区二区| 美女被艹到高潮喷水动态| 狂野欧美激情性xxxx在线观看| 97超碰精品成人国产| 日韩成人av中文字幕在线观看| 少妇高潮的动态图| 最近手机中文字幕大全| 久久精品人妻少妇| 国产精品蜜桃在线观看| 91aial.com中文字幕在线观看| 最近2019中文字幕mv第一页| 国产精品久久久久久久电影| 色综合站精品国产| 亚洲三级黄色毛片| 日韩大片免费观看网站 | 看免费成人av毛片| 性色avwww在线观看| 国产爱豆传媒在线观看| 一级av片app| 久久久成人免费电影| 联通29元200g的流量卡| 日韩强制内射视频| 少妇高潮的动态图| 午夜精品在线福利| 免费看a级黄色片| 亚洲精品乱码久久久v下载方式| 亚洲自偷自拍三级| kizo精华| 免费搜索国产男女视频| 午夜福利在线观看免费完整高清在| 国产麻豆成人av免费视频| 91av网一区二区| 国产 一区 欧美 日韩| 国产一区二区亚洲精品在线观看| 免费看光身美女| 久久久久九九精品影院| 久久精品人妻少妇| 亚洲国产成人一精品久久久| 免费不卡的大黄色大毛片视频在线观看 | 国产精品嫩草影院av在线观看| 久久久午夜欧美精品| 免费黄网站久久成人精品| 中文精品一卡2卡3卡4更新| 亚洲无线观看免费| 男女视频在线观看网站免费| 麻豆久久精品国产亚洲av| 麻豆av噜噜一区二区三区| 亚洲国产精品专区欧美| 久久精品熟女亚洲av麻豆精品 |