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

    非合作無源探測中的參考信號提純新方法

    2022-03-29 07:49:50應(yīng)濤王雪寶田威周成侯小陽
    航空學(xué)報 2022年2期
    關(guān)鍵詞:輻射源無源傅里葉

    應(yīng)濤,王雪寶,田威,周成,侯小陽

    1.海軍工程大學(xué) 電子工程學(xué)院,武漢 430033

    2.海軍士官學(xué)校 電子對抗雷達(dá)聲納系,蚌埠 233012

    非合作無源探測系統(tǒng)是利用第三方非合作輻射源發(fā)射的信號探測目標(biāo)的雙/多基地?zé)o源雷達(dá)系統(tǒng)。該系統(tǒng)本身不發(fā)射能量,而是被動地接收目標(biāo)散射的非合作輻射源電磁信號,對目標(biāo)進(jìn)行探測、定位及跟蹤,具有戰(zhàn)場生存能力強(qiáng)、成本低、反隱身能力強(qiáng)等諸多優(yōu)點(diǎn),因此受到各國軍方重視。目前,非合作無源探測系統(tǒng)利用的第三方輻射源信號除FM 廣播信號、數(shù)字廣播信號、衛(wèi)星信 號、WIFI信號、UMTS 信號以及GSM 信號和LTE 信號等民用輻射源信號以外,還包括己方或敵方非合作雷達(dá)等軍用輻射源信號,英國雷卡公司防御系統(tǒng)分部利用倫敦蓋特威克機(jī)場的遠(yuǎn)程空管雷達(dá)開展了非合作無源探測系統(tǒng)的研究,國防科技大學(xué)選擇預(yù)警機(jī)雷達(dá)信號作為輻射源信號,開展了基于空基輻射源的非合作無源探測系統(tǒng)關(guān)鍵技術(shù)研究。戰(zhàn)時,由于可利用的民用輻射源匱乏,利用己方或敵方雷達(dá)作為輻射源進(jìn)行非合作無源探測將變得尤為重要。線性調(diào)頻(Linear Frequency Modulation,LFM)信號作為一種時寬-帶寬積較大的信號,具有低截獲概率特性和較好的抗干擾能力,廣泛應(yīng)用于雷達(dá)系統(tǒng)。因此,本文將以發(fā)射LFM 信號的雷達(dá)作為非合作輻射源,開展非合作無源探測系統(tǒng)中參考信號的提純研究。

    非合作無源探測系統(tǒng)有2個接收通道,即接收輻射源直達(dá)波作為參考信號的參考通道和接收目標(biāo)散射回波的監(jiān)測通道。參考通道中的直達(dá)波參考信號主要有2個方面的作用:①監(jiān)測通道接收到的目標(biāo)回波中,往往伴隨著大量的直達(dá)波干擾和多徑雜波干擾,且功率遠(yuǎn)大于目標(biāo)回波,這對目標(biāo)檢測造成很大干擾。這時,需要純凈的參考通道直達(dá)波信號作為參考信號,抑制監(jiān)測通道的直達(dá)波干擾和多徑雜波干擾;②參考通道直達(dá)波信號作為二維互相關(guān)檢測的參考信號,參考信號的純凈程度將直接決定系統(tǒng)的檢測性能。但在實(shí)際情況中,參考通道接收到的輻射源直達(dá)波信號不可避免地受到多徑信號的干擾,是輻射源直達(dá)波信號及多徑雜波信號的疊加,參考通道為監(jiān)測通道提供的參考信號不再“純凈”,對系統(tǒng)探測性能造成嚴(yán)重影響。

    目前,獲取“純凈”直達(dá)波參考信號主要有2條技術(shù)思路,即“重構(gòu)”和“提純”?!爸貥?gòu)”是指基于輻射源信號的特殊結(jié)構(gòu)特性,通過信道估計和同步解調(diào)來精確重構(gòu)參考信號,文獻(xiàn)[20]在壓縮感知的匹配追蹤算法基礎(chǔ)上,利用數(shù)字地面多媒體廣播信號的特性,提出了一種參考信號重構(gòu)新方法。“提純”則是利用時域和空域自適應(yīng)濾波等方法,通過抑制參考通道接收信號中的多徑干擾,獲取“純凈”的參考信號,但提取信號純度受限于天線方向性,且對參考通道接收天線提出了更高的要求。

    當(dāng)有用信號與干擾信號在時域或者頻域上重疊時,要想將有用信號從干擾信號中提取出來,通常從有用信號的特性出發(fā),考慮將混合信號變換到有用信號與干擾信號可分離的變換域上,然后采用濾波方法將有用信號提取出來。通常,在進(jìn)行非合作無源探測系統(tǒng)設(shè)計和部署時,會對輻射源信號進(jìn)行調(diào)研和分析,獲取有關(guān)先驗(yàn)知識,因此可對輻射源信號特性進(jìn)行研究,并充分利用輻射源信號的這些特性,將參考通道中直達(dá)波與多徑的混合信號變換到對信號時延敏感的變換域上,在此變換域上將直達(dá)波信號與多徑信號分離開,實(shí)現(xiàn)直達(dá)波信號的提取。從這個思路出發(fā),可以考慮將混合信號變換到對時延敏感的稀疏域上,由于信號在稀疏域上具有有效能量的數(shù)據(jù)點(diǎn)較少,直達(dá)波信號與多徑信號在稀疏域上具有低概率交疊特性,可以想象,最為理想的稀疏變換就是輻射源信號在此稀疏域上表示為一個點(diǎn)的沖激函數(shù),此時,只要該稀疏變換對時延敏感,則直達(dá)波與多徑的混合信號在此稀疏域上的表示將分別對應(yīng)于直達(dá)波和多徑且在稀疏域上相互分離的多個點(diǎn)的沖激函數(shù),則接下來就可以利用濾波方法或參數(shù)估計的方法來抑制多徑信號或恢復(fù)直達(dá)波信號,實(shí)現(xiàn)直達(dá)波信號的提取。本文將從雷達(dá)輻射源的信號特性出發(fā),充分利用線性調(diào)頻信號在分?jǐn)?shù)階傅里葉變換(Fractional Fourier Transformation,FRFT)域上的稀疏特性,根據(jù)變換關(guān)系,通過FRFT 域的尖峰參數(shù)推導(dǎo)出直達(dá)波信號的參數(shù),恢復(fù)直達(dá)波信號,從而實(shí)現(xiàn)直達(dá)波信號提純。

    1 問題描述

    假定系統(tǒng)是岸基固定的,且參考天線主瓣對準(zhǔn)輻射源,則進(jìn)入?yún)⒖纪ǖ赖哪繕?biāo)回波信號十分微弱可以忽略。實(shí)際情況中,參考通道接收到的輻射源直達(dá)波信號伴隨著多徑干擾信號,如圖1所示。假設(shè)在所接收到的參考通道中,存在個多徑信號,則

    圖1 參考信號受多徑干擾示意圖Fig.1 Reference signal interfered by multipath signals

    式中:()為參考通道接收信號;()為輻射源直達(dá)波信號;是參考通道中直達(dá)波信號的復(fù)包絡(luò);()為參考通道中的噪聲;cτ分別為第個多徑信號的復(fù)包絡(luò)和時延,不失一般性,令=1,ττ。

    為了進(jìn)一步分析參考通道中多徑信號對目標(biāo)檢測性能的影響,假設(shè)有2個強(qiáng)目標(biāo)和1個弱目標(biāo),在距離-多普勒曲面分別位于(19.5 km,50 Hz)、(60.0 km,-83 Hz)和(75.0 km,60 Hz),信噪比分別為-16.9、-18.3和-24.4 dB,監(jiān)測通道接收信號中除了3 個目標(biāo)回波外,還包含1個直達(dá)波干擾和10個多徑干擾。1個直達(dá)波干擾的干擾噪聲比是62.3 dB,10個多徑干擾的干擾噪聲比處于(4~34 d B)范圍內(nèi),仿真場景如圖2所示。

    圖2 仿真場景Fig.2 Simulation scenario

    當(dāng)參考信號不受多徑干擾時,采用文獻(xiàn)[22]中的方法,將監(jiān)測通道接收信號投影到多徑干擾的正交補(bǔ)子空間內(nèi),消除強(qiáng)直達(dá)波與強(qiáng)多徑干擾,并進(jìn)行二維互相關(guān)檢測,得到的結(jié)果如圖3所示。從圖3中可以看出,利用“純凈”的直達(dá)波信號作為參考信號進(jìn)行監(jiān)測通道的多徑雜波對消和二維互相關(guān)檢測,3個目標(biāo)可以很容易被檢測到。當(dāng)參考通道接收信號除了輻射源直達(dá)波信號外,還包含10個多徑干擾時,直達(dá)波與多徑雜波功率之比分 別 為2.1、3.2、4.8、6.5、9.2、12.3、11.6、14.9、17.2、19.8 dB,如果直接利用“不純凈”的參考信號進(jìn)行監(jiān)測通道的多徑雜波對消和二維互相關(guān)檢測,得到的結(jié)果如圖4 所示。與圖3 相比,2個強(qiáng)目標(biāo)雖然能夠被檢測到,但由于參考信號的“不純凈”,不能為多徑雜波對消和二維互相關(guān)檢測提供準(zhǔn)確樣本,導(dǎo)致強(qiáng)目標(biāo)產(chǎn)生模糊旁瓣,同時,還出現(xiàn)虛假目標(biāo),且弱目標(biāo)被抑制,很難被檢測到。因此,參考通道直達(dá)波信號提純對于系統(tǒng)目標(biāo)檢測性能十分重要。

    圖3 參考信號不含多徑雜波時的互模糊函數(shù)Fig.3 Cross correlation function for a reference signal without multipath

    圖4 參考信號含多徑信號時的互模糊函數(shù)Fig.4 Cross correlation function for a reference signal with multipath

    2 參考信號的稀疏性

    從數(shù)學(xué)角度,分?jǐn)?shù)階傅里葉變換可被理解為時頻平面內(nèi)坐標(biāo)軸繞原點(diǎn)逆時針旋轉(zhuǎn)任意角度后構(gòu)成的表示方法,是傳統(tǒng)傅里葉變換的廣義形式,可作為信號時頻分析的一種方法。

    從積分核的角度出發(fā),輸入函數(shù)()的階FRFT 定義為

    本文以發(fā)射LFM 脈沖信號的雷達(dá)作為非合作輻射源,非合作無源探測系統(tǒng)參考通道接收的直達(dá)波信號可表示為

    式中:為直達(dá)波信號包絡(luò);為起始頻率;=/為調(diào)頻斜率,為信號持續(xù)時間,為信號調(diào)頻帶寬。不失一般性,取=0,=1,則信號階FRFT 為

    將式(5)變量代換后簡化為

    由式(6)可看出,當(dāng)變換階次=arccot(-)/(π/2)時,直達(dá)波信號在FRFT 域上表現(xiàn)出了良好的能量聚集稀疏特性。也就是說,只要找到與直達(dá)波信號能量聚集特性對應(yīng)的最佳階次,與之對應(yīng)的FRFT域的正交基字典就可以被構(gòu)造出來,在該字典下,直達(dá)波信號將滿足稀疏性和正交性。FRFT的周期為4,∈(-2,2],對應(yīng)旋轉(zhuǎn)角度∈(-π,π],由FRFT 的時頻旋轉(zhuǎn)特性可知,當(dāng)變換階次從0變化到1,即可展示出信號從時域轉(zhuǎn)換到頻域的所有特征,根據(jù)對稱性,只要觀察∈(-1,1]就可找到符合條件的最佳變換階次。

    正交基字典矩陣是能否將信號進(jìn)行稀疏表示的關(guān)鍵,由于直達(dá)波信號是在FRFT 基上分解信號的,可考慮用離散分?jǐn)?shù)階傅里葉變換(Discrete Fractional Fourier Transformation,DFRFT)核矩陣來構(gòu)成正交字典K ,由上述分析可知,K 能很好地稀疏表示直達(dá)波信號,不同的變換階次對應(yīng)不同的核變換矩陣。下面將根據(jù)DFRFT 核構(gòu)造出正交基字典K ,并證明直達(dá)波信號在K 下的正交性和稀疏性。

    考慮直接采樣連續(xù)分?jǐn)?shù)階傅里葉變換核來獲得DFRFT核矩陣,進(jìn)一步由離散核矩陣來構(gòu)造FRFT正交基字典。這類方法中,Pei采樣型算法具有計算復(fù)雜度最低、運(yùn)算效率最高等特點(diǎn),得到了廣泛的使用。為了使離散分?jǐn)?shù)階傅里葉變換具備正交性和可逆性,其對連續(xù)分?jǐn)?shù)階傅里葉變換在分?jǐn)?shù)階傅里葉域和時域選擇合適的采樣間隔。本文基于分?jǐn)?shù)階傅里葉變換的正交基字典構(gòu)造方法步驟如下:

    對輸入函數(shù)()和輸出函數(shù)X ()分別以Δ和Δ為采樣間隔進(jìn)行采樣,即()=(Δ),Y ()=X (Δ),∈[-,],∈[-,]。

    由連續(xù)FRFT 定義式(2)可知

    由變換可逆性可知,式(9)對求和等于(-),即

    式中:是與2+1互質(zhì)的整數(shù)。

    將式(10)代入式(8)歸一化得

    令=sgn(sin)=±1,當(dāng)sin>0時,代入式(11)得

    采用正交基字典,令=,就可以得到所要構(gòu)造的DFRFT 正交基字典

    進(jìn)一步,可以得出利用FRFT 核矩陣構(gòu)成的字典K 也是參考通道接收信號()的正交基字典,證明詳見附錄A。

    3 參考信號提純

    由第2節(jié)的分析可知,通過對參考通道接收信號進(jìn)行FRFT,可將時域上混疊在一起的直達(dá)波信號和多徑信號在稀疏域上分離開,通常考慮采用濾波器技術(shù)就可將其他多徑信號抑制掉。但實(shí)際情況中,很難設(shè)計出一種理想的濾波器將多徑信號完全抑制的同時還能盡可能完整地保留直達(dá)波信號,尤其在多徑時延很小的情況。因此,本節(jié)考慮從FRFT 變換關(guān)系出發(fā),由參考通道接收信號在FRFT 域的尖峰參數(shù)推導(dǎo)出直達(dá)波的參數(shù),重構(gòu)提純直達(dá)波信號。

    由參考通道接收信號模型式(1)和式(3)可知,多徑信號需要估計的參數(shù)有:調(diào)頻斜率、多徑信號包絡(luò)c和時延τ,下面給出這些參數(shù)的估計。

    從分?jǐn)?shù)階傅里葉變換的定義出發(fā)進(jìn)行推導(dǎo)可以得出,直達(dá)波信號的包絡(luò)估計為

    得到直達(dá)波信號的參數(shù)估計后,就可以參考信號提純

    綜上所述,基于稀疏特性的直達(dá)波參考信號提取方法步驟如下:

    連續(xù)變化旋轉(zhuǎn)角度或者階次,對參考通道接收信號()進(jìn)行分?jǐn)?shù)階傅里葉變換,形成參考通道接收信號的能量在(,)平面上的三維分布圖。

    利用步驟3中估計出的直達(dá)波信號參數(shù),根據(jù)式(17)實(shí)現(xiàn)直達(dá)波信號提純。

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

    為了驗(yàn)證基于稀疏特性的直達(dá)波信號提純算法的有效性,進(jìn)行相應(yīng)的仿真實(shí)驗(yàn)。對雷達(dá)輻射源LFM 直達(dá)波信號仿真參數(shù)設(shè)置如下:信號帶寬=5 MHz,觀測時長=8μs,LFM 信號的調(diào)頻斜率=/=6.25×10Hz/s,信號采樣頻率=100 MHz,采樣間隔=1/=10s,信號采樣點(diǎn)數(shù)為800,信號歸一化幅度為1,信噪比SNR=8 dB的直達(dá)波信號如圖5所示。用不同的延時=1μs和=2μs以及對應(yīng)的幅度=05和=08來模擬2條多徑信號,參考通道中直達(dá)波與多徑信號的混合信號如圖6所示。

    圖5 直達(dá)波信號(SNR=8 dB)Fig.5 Direct signal(SNR=8 dB)

    圖6 參考通道接收的混合信號Fig.6 Mixture of direct signal and multipath signals received by reference channel

    圖7 參考通道接收信號的FRFT 三維圖Fig.7 3D plot of FRFT of the received signal in reference channel

    圖8 DFRFT 字典矩陣Fig.8 Dictionary matrix of DFRFT

    圖9 參考通道接收信號在最佳旋轉(zhuǎn)角度下的FRFTFig.9 FRFT of received signal in reference channel with the best rotation angle

    從分?jǐn)?shù)階傅里葉變換與輸入信號的對應(yīng)關(guān)系出發(fā),根據(jù)第3節(jié)中的參數(shù)估計方法消除多徑信號,實(shí)現(xiàn)參考信號提純后,參考通道接收信號的FRFT 三維圖和在最佳旋轉(zhuǎn)角度下的FRFT 如圖10和圖11所示。不難看出,參考通道接收信號中的多徑信號得到了很好地抑制,直達(dá)波信號得到了很好地保留,提取出的直達(dá)波波形與純凈的直達(dá)波波形十分接近,如圖12所示。

    圖10 參考通道接收信號的FRFT 三維圖(多徑信號消除后)Fig.10 3D plot of FRFT of received signal in reference channel after multipath removal

    圖11 參考通道接收信號在最佳旋轉(zhuǎn)角度下的FRFT(多徑信號消除后)Fig.11 FRFT of received signal in reference channel with the best rotation angle after multipath removal

    圖12 多徑信號消除后的直達(dá)波信號Fig.12 Direct signal after multipath removal

    為了進(jìn)一步定量評價本文方法參考信號提純的性能,定義波形失真率為

    式中:()為純凈的直達(dá)波信號;()為提取的直達(dá)波信號,波形失真率表示提取的直達(dá)波信號與純凈的直達(dá)波信號之間的相似程度。同時,與需要采用多通道接收機(jī)的空域自適應(yīng)濾波方法相比,本文算法可適用于單通道接收機(jī),無需增加接收機(jī)硬件復(fù)雜度,故將本文方法與時域自適應(yīng)濾波方法進(jìn)行仿真性能對比,圖13給出了不同信噪比下2 種方法提取的直達(dá)波信號波形失真率??梢钥闯?本文方法提純信號的波形失真率比時域自適應(yīng)濾波方法要小,直達(dá)波參考信號提純效果更好。隨著信噪比的增大,采用本文方法提取的直達(dá)波信號波形失真率在不斷減小,提取的直達(dá)波信號精度越高。這是由于隨著信噪比的增大,參數(shù)估計的精度會提高,字典的稀疏表達(dá)能力越強(qiáng),多徑信號的參數(shù)估計精度也越高,因此提取的直達(dá)波信號波形失真率減小,直達(dá)波信號提純效果越好。

    圖13 不同信噪比下2種方法的波形失真率Fig.13 Waveform Distortion rate of the two methods with different SNRs

    由于非合作無源探測中的目標(biāo)檢測是在參考信號與監(jiān)測信號進(jìn)行匹配濾波的基礎(chǔ)之上進(jìn)行的,波形失真率越低,參考信號提純效果越好,則匹配濾波獲得的目標(biāo)信號增益越大,越有利于目標(biāo)檢測。因此,匹配濾波后的目標(biāo)信號峰值信噪比可衡量目標(biāo)檢測性能的好壞。為便于分析,假定有一個目標(biāo)距離為15 km 的固定目標(biāo),圖14和圖15分別給出了波形失真率分別為62%和28%時,提純的參考信號與目標(biāo)信號進(jìn)行匹配濾波的結(jié)果。從圖中不難看出,當(dāng)波形失真率越低,參考信號提純效果越好,匹配濾波獲得的目標(biāo)信號增益越高,目標(biāo)信號信噪比就越高,目標(biāo)就越容易被檢測出來。

    圖14 波形失真率為62%時匹配濾波輸出Fig.14 Output of match filter of waveform distortion rate of 62%

    圖15 波形失真率為28%時匹配濾波輸出Fig.15 Output of match filter of waveform distortion rate of 28%

    5 總 結(jié)

    針對非合作無源探測中的參考信號提純問題,本文以發(fā)射線性調(diào)頻脈沖信號的雷達(dá)作為非合作輻射源,提出了一種基于稀疏特性的參考信號提純新方法。仿真實(shí)驗(yàn)結(jié)果表明,本文方法能夠充分利用LFM 信號在分?jǐn)?shù)階傅里葉變換域上的稀疏特性,將直達(dá)波信號與多徑信號在稀疏域上分離,通過參數(shù)估計實(shí)現(xiàn)參考信號有效提純,且與時域自適應(yīng)濾波方法相比,提純的直達(dá)波參考信號失真更小。雖然,本文方法主要適用于基于LFM 脈沖雷達(dá)輻射源的非合作無源探測系統(tǒng),但本文方法的技術(shù)思路仍可參考借鑒應(yīng)用于利用其他輻射源信號的系統(tǒng),關(guān)鍵是尋找到合適的稀疏域。

    利用FRFT 核矩陣構(gòu)成的字典K 是參考通道接收信號()的正交基字典。

    證明 本命題可分解為2個子命題:①證明字典K 是正交字典;②證明多徑信號在K 上是稀疏的。

    1)根據(jù)FRFT 定義式,K 是由一組Chirp信號展開成的滿足正交條件的完備正交基,在離散核矩陣來構(gòu)造FRFT 正交基字典過程中,各個基函數(shù)之間相互正交。

    2)令參考通道接收信號中直達(dá)波信號FRFT 為X (),當(dāng)變換階次滿足=arccot(-)/(π/2)時,大小為(),具有稀疏性。

    令=-,=-cos代入式(A1)得

    式(A2)化簡可得

    令=-cos,換元后可得

    可以看出,多徑信號的FRFT 與X ()只存在著不同的相移和時移因子,且同一階數(shù)下FRFT 將具有相同的模。

    選取最佳階次=arccot(-)/(π/2),=π/2,有X ()=(),則F [(-)]是一個經(jīng)過一定時移和相移后的沖激函數(shù),且域時移與域時延對應(yīng)關(guān)系為=cos,也具有稀疏性,證明完畢。

    猜你喜歡
    輻射源無源傅里葉
    基于博弈論的GRA-TOPSIS輻射源威脅評估方法
    一種三相無源逆變電源供電方案設(shè)計
    電子制作(2019年12期)2019-07-16 08:45:14
    雙線性傅里葉乘子算子的量化加權(quán)估計
    基于小波降噪的稀疏傅里葉變換時延估計
    數(shù)字電視外輻射源雷達(dá)多旋翼無人機(jī)微多普勒效應(yīng)實(shí)驗(yàn)研究
    外輻射源雷達(dá)直升機(jī)旋翼參數(shù)估計方法
    基于PCH模型的航天器姿態(tài)無源控制
    基于遷移成分分析的雷達(dá)輻射源識別方法研究
    無源互調(diào)干擾對TD-LTE系統(tǒng)的影響研究
    基于傅里葉變換的快速TAMVDR算法
    九色成人免费人妻av| 免费观看性生交大片5| 2022亚洲国产成人精品| 国产精品成人在线| 亚洲精品久久久久久婷婷小说| 最新的欧美精品一区二区| 国产午夜精品久久久久久一区二区三区| 国产熟女午夜一区二区三区 | 欧美精品国产亚洲| 精品一区二区免费观看| 国产在视频线精品| 中国国产av一级| 久久国产亚洲av麻豆专区| 国产 精品1| 欧美成人午夜免费资源| 丝袜喷水一区| 久久狼人影院| 亚洲经典国产精华液单| 男女啪啪激烈高潮av片| 欧美日韩精品成人综合77777| 成人国产av品久久久| 国产日韩欧美在线精品| 国产精品99久久久久久久久| 亚洲精品美女久久av网站| 大香蕉久久网| 亚洲欧美一区二区三区黑人 | 999精品在线视频| 亚洲av在线观看美女高潮| 极品少妇高潮喷水抽搐| 下体分泌物呈黄色| 秋霞在线观看毛片| 亚洲图色成人| 男人添女人高潮全过程视频| 一级毛片aaaaaa免费看小| 午夜免费观看性视频| 国产精品99久久久久久久久| 人人妻人人澡人人爽人人夜夜| 大话2 男鬼变身卡| 日本与韩国留学比较| 国产欧美日韩综合在线一区二区| 精品国产一区二区久久| 大话2 男鬼变身卡| 国产成人精品无人区| 制服诱惑二区| 免费不卡的大黄色大毛片视频在线观看| 久久人人爽人人爽人人片va| 青春草亚洲视频在线观看| 美女福利国产在线| a级毛片黄视频| 婷婷色av中文字幕| 亚洲av综合色区一区| 国产一区有黄有色的免费视频| 热99久久久久精品小说推荐| 性高湖久久久久久久久免费观看| 国产精品久久久久久久电影| 永久网站在线| 久久久午夜欧美精品| 欧美bdsm另类| 午夜激情久久久久久久| 国产黄片视频在线免费观看| 午夜老司机福利剧场| 熟女人妻精品中文字幕| 狂野欧美白嫩少妇大欣赏| 大又大粗又爽又黄少妇毛片口| 肉色欧美久久久久久久蜜桃| 各种免费的搞黄视频| 大片免费播放器 马上看| 亚洲国产成人一精品久久久| 高清视频免费观看一区二区| 观看av在线不卡| 综合色丁香网| 亚洲美女搞黄在线观看| 国产成人精品在线电影| 久久精品久久久久久噜噜老黄| 日本av手机在线免费观看| 建设人人有责人人尽责人人享有的| 少妇的逼水好多| 日本-黄色视频高清免费观看| 中国三级夫妇交换| 免费日韩欧美在线观看| 成人综合一区亚洲| 91久久精品电影网| 人人澡人人妻人| 亚洲欧美日韩另类电影网站| 麻豆成人av视频| 国产精品成人在线| 精品少妇久久久久久888优播| 人妻少妇偷人精品九色| 伦精品一区二区三区| 亚洲精品久久久久久婷婷小说| av一本久久久久| 少妇熟女欧美另类| 91精品三级在线观看| 街头女战士在线观看网站| 亚洲av电影在线观看一区二区三区| 99热6这里只有精品| av视频免费观看在线观看| 视频在线观看一区二区三区| 人人妻人人添人人爽欧美一区卜| 中文字幕免费在线视频6| 成年人免费黄色播放视频| 国产黄频视频在线观看| 成人影院久久| 卡戴珊不雅视频在线播放| 99热6这里只有精品| 精品亚洲成a人片在线观看| kizo精华| 亚洲欧美精品自产自拍| 男女高潮啪啪啪动态图| 少妇人妻精品综合一区二区| 免费黄色在线免费观看| 乱码一卡2卡4卡精品| 插逼视频在线观看| 午夜久久久在线观看| 久久精品夜色国产| 久久久久国产网址| 国产成人免费观看mmmm| 97精品久久久久久久久久精品| 另类亚洲欧美激情| 亚洲美女黄色视频免费看| 久久久久久久精品精品| 亚洲精品一区蜜桃| 久久免费观看电影| 国产成人av激情在线播放 | 色网站视频免费| av又黄又爽大尺度在线免费看| 亚洲欧美中文字幕日韩二区| 国产毛片在线视频| 国产一区二区在线观看日韩| 午夜91福利影院| 在线观看人妻少妇| 日产精品乱码卡一卡2卡三| 青青草视频在线视频观看| 婷婷成人精品国产| 国产av精品麻豆| 18禁裸乳无遮挡动漫免费视频| 女性被躁到高潮视频| 久久久精品免费免费高清| 18在线观看网站| 观看美女的网站| 18禁观看日本| 国产免费一级a男人的天堂| 午夜久久久在线观看| 99热这里只有是精品在线观看| 久久精品国产亚洲av涩爱| 少妇丰满av| 午夜老司机福利剧场| 黄色欧美视频在线观看| 国产精品无大码| 久久 成人 亚洲| 久久久久久久久久久丰满| 国产午夜精品一二区理论片| 2022亚洲国产成人精品| 久久精品久久精品一区二区三区| 久久久久久久亚洲中文字幕| 我的老师免费观看完整版| 日本色播在线视频| 99久国产av精品国产电影| 亚洲精品中文字幕在线视频| 色5月婷婷丁香| 最新中文字幕久久久久| 精品一区二区三区视频在线| 人妻制服诱惑在线中文字幕| 又粗又硬又长又爽又黄的视频| 狂野欧美白嫩少妇大欣赏| 精品少妇黑人巨大在线播放| 一本色道久久久久久精品综合| 亚洲av免费高清在线观看| 久久这里有精品视频免费| 中文欧美无线码| 大片免费播放器 马上看| 久久国产精品男人的天堂亚洲 | 伊人久久精品亚洲午夜| 精品久久蜜臀av无| 999精品在线视频| 国产高清不卡午夜福利| 国国产精品蜜臀av免费| 熟妇人妻不卡中文字幕| 亚洲一区二区三区欧美精品| 国产精品久久久久久av不卡| 精品视频人人做人人爽| 成人国产麻豆网| av不卡在线播放| 少妇被粗大的猛进出69影院 | 九九久久精品国产亚洲av麻豆| 最近手机中文字幕大全| 亚洲久久久国产精品| 精品国产一区二区久久| 三级国产精品片| av有码第一页| 97在线视频观看| 亚洲精品456在线播放app| 欧美日韩综合久久久久久| 色哟哟·www| 国产精品久久久久久精品古装| 久久久久久久久久人人人人人人| 一级二级三级毛片免费看| 精品99又大又爽又粗少妇毛片| 日本免费在线观看一区| 色视频在线一区二区三区| 亚洲精品一区蜜桃| 亚洲美女黄色视频免费看| 能在线免费看毛片的网站| 男女无遮挡免费网站观看| 欧美日韩在线观看h| videossex国产| h视频一区二区三区| 男人操女人黄网站| 国产黄色视频一区二区在线观看| 久久久久久久久久人人人人人人| 天堂8中文在线网| 伊人久久国产一区二区| 国产精品无大码| 一本久久精品| 日韩,欧美,国产一区二区三区| 中国三级夫妇交换| 中文字幕制服av| 永久网站在线| av卡一久久| 在线观看免费视频网站a站| 亚洲天堂av无毛| 视频中文字幕在线观看| 国产成人aa在线观看| 一本久久精品| 狂野欧美激情性xxxx在线观看| 老司机影院成人| 97在线人人人人妻| 色网站视频免费| xxx大片免费视频| a级毛片黄视频| h视频一区二区三区| 亚洲精品亚洲一区二区| 有码 亚洲区| 黄色视频在线播放观看不卡| 久久婷婷青草| 黄色欧美视频在线观看| 999精品在线视频| 久久国产精品大桥未久av| 国产极品粉嫩免费观看在线 | 国产精品久久久久久av不卡| 久久精品国产a三级三级三级| .国产精品久久| av国产精品久久久久影院| 在线观看人妻少妇| 制服丝袜香蕉在线| 夫妻性生交免费视频一级片| 97超碰精品成人国产| 丝袜喷水一区| 久久久a久久爽久久v久久| 中文字幕制服av| 999精品在线视频| 午夜精品国产一区二区电影| 一边摸一边做爽爽视频免费| 久久久久精品性色| 久久久久国产精品人妻一区二区| 久久久久久久精品精品| 亚洲精品久久午夜乱码| 午夜福利视频精品| 午夜福利影视在线免费观看| 大香蕉97超碰在线| 曰老女人黄片| 久久久久国产精品人妻一区二区| 国产视频内射| 999精品在线视频| 欧美日韩成人在线一区二区| 成人免费观看视频高清| 久久99蜜桃精品久久| 亚洲国产精品专区欧美| 一级毛片 在线播放| 男女边吃奶边做爰视频| 天美传媒精品一区二区| 国产在视频线精品| 美女国产高潮福利片在线看| 韩国高清视频一区二区三区| 九九爱精品视频在线观看| 天天影视国产精品| 欧美老熟妇乱子伦牲交| 精品酒店卫生间| 亚洲欧美一区二区三区黑人 | 极品人妻少妇av视频| 大片电影免费在线观看免费| 中国国产av一级| 日日摸夜夜添夜夜爱| 男人爽女人下面视频在线观看| 建设人人有责人人尽责人人享有的| 国产欧美日韩一区二区三区在线 | 汤姆久久久久久久影院中文字幕| 久久精品国产自在天天线| 一本—道久久a久久精品蜜桃钙片| 大片电影免费在线观看免费| 日日爽夜夜爽网站| 人妻人人澡人人爽人人| 亚洲精品成人av观看孕妇| 蜜桃久久精品国产亚洲av| 免费少妇av软件| 欧美+日韩+精品| 99国产综合亚洲精品| 国产一区二区三区av在线| 日韩大片免费观看网站| 欧美另类一区| 美女福利国产在线| 婷婷成人精品国产| 日韩在线高清观看一区二区三区| 啦啦啦视频在线资源免费观看| 一二三四中文在线观看免费高清| a级毛片黄视频| 亚洲国产欧美日韩在线播放| 亚洲丝袜综合中文字幕| 夫妻性生交免费视频一级片| 精品国产乱码久久久久久小说| 亚洲国产精品一区三区| 春色校园在线视频观看| 婷婷色av中文字幕| 国产精品国产av在线观看| 欧美日韩一区二区视频在线观看视频在线| 国产欧美日韩综合在线一区二区| a级毛片免费高清观看在线播放| av电影中文网址| 性色av一级| 日韩精品免费视频一区二区三区 | 大香蕉97超碰在线| 国产深夜福利视频在线观看| 国产精品国产av在线观看| 久久久久久久久久久久大奶| 中国三级夫妇交换| 亚洲,一卡二卡三卡| 啦啦啦啦在线视频资源| 午夜免费男女啪啪视频观看| 制服丝袜香蕉在线| 亚洲伊人久久精品综合| 视频中文字幕在线观看| 在线观看www视频免费| 亚洲精华国产精华液的使用体验| 亚洲精品aⅴ在线观看| 18禁裸乳无遮挡动漫免费视频| 91精品国产国语对白视频| 如何舔出高潮| 久久精品国产亚洲av涩爱| 这个男人来自地球电影免费观看 | 国产 一区精品| av一本久久久久| 最近2019中文字幕mv第一页| 美女中出高潮动态图| 人体艺术视频欧美日本| 久久99蜜桃精品久久| 视频区图区小说| 精品久久久久久电影网| 99热网站在线观看| 免费不卡的大黄色大毛片视频在线观看| 精品国产一区二区久久| freevideosex欧美| 99视频精品全部免费 在线| 国产精品不卡视频一区二区| 女人久久www免费人成看片| 欧美老熟妇乱子伦牲交| 精品少妇黑人巨大在线播放| 免费高清在线观看日韩| 国产免费现黄频在线看| 夫妻性生交免费视频一级片| 最近中文字幕2019免费版| 高清午夜精品一区二区三区| 少妇丰满av| 亚洲,欧美,日韩| 一级毛片黄色毛片免费观看视频| 精品久久久久久电影网| 欧美 亚洲 国产 日韩一| 大又大粗又爽又黄少妇毛片口| 免费高清在线观看视频在线观看| 色视频在线一区二区三区| 熟妇人妻不卡中文字幕| 亚洲国产日韩一区二区| 亚洲精品日韩av片在线观看| 国产成人a∨麻豆精品| 亚洲精品乱码久久久v下载方式| 亚洲欧美精品自产自拍| 国产片内射在线| 精品人妻熟女av久视频| 91精品一卡2卡3卡4卡| 一区在线观看完整版| 少妇猛男粗大的猛烈进出视频| 青春草国产在线视频| 女性生殖器流出的白浆| xxx大片免费视频| 国产乱来视频区| 黄色配什么色好看| 日韩中字成人| 一本大道久久a久久精品| 夜夜爽夜夜爽视频| 亚洲av国产av综合av卡| 美女xxoo啪啪120秒动态图| 久久久久久人妻| 人体艺术视频欧美日本| 2018国产大陆天天弄谢| 少妇熟女欧美另类| 国产精品一区二区在线观看99| 国产欧美另类精品又又久久亚洲欧美| a 毛片基地| 欧美+日韩+精品| 精品国产一区二区久久| 免费日韩欧美在线观看| 久久亚洲国产成人精品v| 亚洲欧洲国产日韩| 在线 av 中文字幕| 尾随美女入室| 美女国产视频在线观看| 精品一区二区三区视频在线| 中文字幕久久专区| 蜜臀久久99精品久久宅男| 久久 成人 亚洲| tube8黄色片| 久久免费观看电影| av在线播放精品| 一级爰片在线观看| 亚洲精品自拍成人| av国产精品久久久久影院| 久久久久久久亚洲中文字幕| 日韩人妻高清精品专区| 黄色怎么调成土黄色| 亚洲国产精品国产精品| 欧美精品人与动牲交sv欧美| 天堂8中文在线网| 自线自在国产av| 国产免费福利视频在线观看| 国产亚洲最大av| 蜜桃久久精品国产亚洲av| 成人毛片a级毛片在线播放| 国产成人av激情在线播放 | 精品一品国产午夜福利视频| 日产精品乱码卡一卡2卡三| 黄色毛片三级朝国网站| 精品亚洲成国产av| 婷婷色麻豆天堂久久| 亚洲成人手机| 大又大粗又爽又黄少妇毛片口| 成人手机av| 欧美日韩国产mv在线观看视频| av专区在线播放| 国产精品无大码| 国产精品久久久久成人av| 国产成人午夜福利电影在线观看| 视频区图区小说| 午夜激情福利司机影院| 十八禁网站网址无遮挡| 亚洲国产精品专区欧美| 男人爽女人下面视频在线观看| 国产在线视频一区二区| 亚洲综合色惰| 国产乱来视频区| 最近中文字幕2019免费版| 精品久久久久久久久av| 免费看av在线观看网站| 熟女电影av网| 中文字幕最新亚洲高清| av黄色大香蕉| 欧美老熟妇乱子伦牲交| 成年女人在线观看亚洲视频| 91国产中文字幕| 夜夜骑夜夜射夜夜干| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 精品久久国产蜜桃| 啦啦啦啦在线视频资源| 成年美女黄网站色视频大全免费 | 日本欧美视频一区| 国产精品一国产av| 五月天丁香电影| 国产探花极品一区二区| 黑人高潮一二区| 999精品在线视频| 男女免费视频国产| 亚洲精品av麻豆狂野| 赤兔流量卡办理| 少妇人妻久久综合中文| 午夜激情福利司机影院| 欧美日韩亚洲高清精品| 久久久国产精品麻豆| 国产精品三级大全| 久久久久久久久久久丰满| 一区二区三区精品91| 夜夜骑夜夜射夜夜干| 欧美激情极品国产一区二区三区 | 精品久久久久久久久亚洲| 国产精品嫩草影院av在线观看| 99国产精品免费福利视频| 成人国语在线视频| 晚上一个人看的免费电影| 满18在线观看网站| 男男h啪啪无遮挡| 欧美人与善性xxx| 五月玫瑰六月丁香| 亚洲av不卡在线观看| 日产精品乱码卡一卡2卡三| 在线 av 中文字幕| 亚洲国产日韩一区二区| 777米奇影视久久| 欧美成人午夜免费资源| 国产精品 国内视频| 一级毛片我不卡| 亚洲国产精品999| 久久人人爽人人爽人人片va| 国产深夜福利视频在线观看| 一区二区三区精品91| 婷婷色综合www| 色94色欧美一区二区| 国产综合精华液| 一级毛片aaaaaa免费看小| 搡女人真爽免费视频火全软件| 久久精品国产亚洲网站| 制服丝袜香蕉在线| 久久影院123| 少妇熟女欧美另类| 伦精品一区二区三区| 婷婷色综合大香蕉| 卡戴珊不雅视频在线播放| 日日啪夜夜爽| 日本黄大片高清| 欧美另类一区| 在线观看www视频免费| 麻豆成人av视频| 26uuu在线亚洲综合色| 午夜91福利影院| 又黄又爽又刺激的免费视频.| 欧美日韩一区二区视频在线观看视频在线| 中文精品一卡2卡3卡4更新| 97精品久久久久久久久久精品| 熟女av电影| 男女边摸边吃奶| 亚洲国产精品一区三区| 亚洲美女搞黄在线观看| 狠狠婷婷综合久久久久久88av| 午夜久久久在线观看| 亚洲国产欧美日韩在线播放| 国产精品不卡视频一区二区| 色哟哟·www| 99久久人妻综合| 少妇丰满av| 香蕉精品网在线| 街头女战士在线观看网站| 一级黄片播放器| 美女大奶头黄色视频| av女优亚洲男人天堂| 中文字幕人妻丝袜制服| 午夜老司机福利剧场| 99精国产麻豆久久婷婷| 王馨瑶露胸无遮挡在线观看| videosex国产| 卡戴珊不雅视频在线播放| 黄色配什么色好看| 久久久久久久精品精品| kizo精华| 激情五月婷婷亚洲| 赤兔流量卡办理| 日韩av免费高清视频| 亚洲第一av免费看| 在线观看三级黄色| 国产成人精品一,二区| av又黄又爽大尺度在线免费看| videos熟女内射| 我的女老师完整版在线观看| 丝袜喷水一区| 亚洲国产欧美在线一区| 亚洲欧美色中文字幕在线| 亚洲av成人精品一二三区| 久久精品国产鲁丝片午夜精品| 亚洲伊人久久精品综合| 国产亚洲精品第一综合不卡 | av福利片在线| 飞空精品影院首页| 五月伊人婷婷丁香| 国产欧美日韩一区二区三区在线 | 亚洲精品第二区| 日韩一本色道免费dvd| 九色亚洲精品在线播放| 中文天堂在线官网| 狂野欧美白嫩少妇大欣赏| 天天躁夜夜躁狠狠久久av| 国产黄色视频一区二区在线观看| 亚洲av中文av极速乱| 九色亚洲精品在线播放| 亚洲精品视频女| 妹子高潮喷水视频| 日韩中字成人| 自线自在国产av| 亚洲美女视频黄频| 伊人久久国产一区二区| 一本—道久久a久久精品蜜桃钙片| 婷婷色av中文字幕| av免费观看日本| 久久精品国产a三级三级三级| 我的女老师完整版在线观看| 午夜影院在线不卡| 一级片'在线观看视频| 久久精品熟女亚洲av麻豆精品| 一级片'在线观看视频| 成人午夜精彩视频在线观看| 97在线视频观看| 日日摸夜夜添夜夜添av毛片| 免费观看性生交大片5| 人人妻人人澡人人爽人人夜夜| 美女内射精品一级片tv| 美女脱内裤让男人舔精品视频| 美女内射精品一级片tv| 欧美三级亚洲精品| 人妻少妇偷人精品九色| 成人毛片a级毛片在线播放| 一区二区三区免费毛片| 日本黄色日本黄色录像| 卡戴珊不雅视频在线播放| 国产精品.久久久| 美女内射精品一级片tv| 大陆偷拍与自拍| 91精品国产国语对白视频| 中国美白少妇内射xxxbb| 一边亲一边摸免费视频| av有码第一页| 午夜91福利影院|