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

    基于熵理論的間歇采樣轉(zhuǎn)發(fā)干擾識別方法*

    2016-07-21 06:52:47蔣瑩何明浩郁春來王冰切
    現(xiàn)代防御技術(shù) 2016年3期
    關(guān)鍵詞:支持向量機信息熵特征提取

    蔣瑩,何明浩,郁春來,王冰切

    (空軍預(yù)警學(xué)院,湖北 武漢 430019)

    ?

    探測跟蹤技術(shù)

    基于熵理論的間歇采樣轉(zhuǎn)發(fā)干擾識別方法*

    蔣瑩,何明浩,郁春來,王冰切

    (空軍預(yù)警學(xué)院,湖北 武漢430019)

    摘要:干擾識別歷來是抗干擾流程中的關(guān)鍵環(huán)節(jié),針對間歇采樣轉(zhuǎn)發(fā)干擾的識別問題,提出了一種基于熵理論的干擾識別方法,將熵的概念應(yīng)用于干擾信號的特征提取,并使用支持向量機進行分類識別。對干擾及回波信號進行FFT變換及預(yù)處理后,提取信號頻譜的3種熵特征,并進行仿真實驗檢驗熵特征的性能。仿真結(jié)果表明,熵特征受噪聲影響較小,將其作為特征參數(shù)進行干擾識別均能取得較為理想的識別結(jié)果,其中信息熵作為特征參數(shù)時的識別性能最優(yōu),指數(shù)熵次之,范數(shù)熵識別性能最差。

    關(guān)鍵詞:間歇采樣;干擾識別;特征提取;信息熵;指數(shù)熵;范數(shù)熵;支持向量機

    0引言

    隨著數(shù)字射頻存儲器(digitalradiofrequencymemory,DRFM)在干擾機中的應(yīng)用愈加廣泛,雷達(dá)干擾技術(shù)發(fā)展迅猛,新的干擾樣式不斷涌現(xiàn),間歇采樣轉(zhuǎn)發(fā)干擾尤為典型。工作在該模式下的干擾機,對接收到的雷達(dá)信號進行間接采樣,短暫存儲后直接轉(zhuǎn)發(fā)回雷達(dá),可以使雷達(dá)產(chǎn)生多個逼真假目標(biāo),從而喪失跟蹤檢測能力,其干擾性能優(yōu)越,相關(guān)針對性抗干擾技術(shù)的研究迫在眉睫[1-4]。

    眾所周知,干擾識別歷來是抗干擾技術(shù)中的關(guān)鍵環(huán)節(jié),正確地進行干擾類型的識別是采取針對性抗干擾措施的前提和基礎(chǔ)。文獻(xiàn)[5]提出了基于雙譜特征和模式識別技術(shù)的欺騙式干擾識別方法;文獻(xiàn)[6]以干擾誤差角為測度衡量了距離門內(nèi)有或無距離欺騙干擾2種情形下信號的相似程度;文獻(xiàn)[7]引入了短時傅里葉變換時頻分析方法以識別DRFM復(fù)制的干擾信號;文獻(xiàn)[8]提取了干擾信號的小波分解能量比作為干擾識別的特征參數(shù)。隨著現(xiàn)代信號處理方法在抗干擾領(lǐng)域的應(yīng)用,干擾識別方法不斷創(chuàng)新發(fā)展,但針對間歇采樣轉(zhuǎn)發(fā)干擾識別的研究目前仍較為缺乏。

    本文針對間歇采樣轉(zhuǎn)發(fā)干擾的識別問題展開研究,提出了基于熵理論的間歇采樣轉(zhuǎn)發(fā)干擾識別方法,利用干擾信號及目標(biāo)回波頻譜上的差異,提取信號頻譜的信息熵、指數(shù)熵以及范數(shù)熵特征,并使用支持向量機進行分類識別,仿真結(jié)果表明,3種熵特征作為干擾識別的特征參數(shù)均能取得較為理想的識別結(jié)果,其中信息熵作為特征參數(shù)時識別性能最優(yōu),指數(shù)熵次之,范數(shù)熵識別性能最差。

    1信號模型

    設(shè)雷達(dá)發(fā)射的信號為線性調(diào)頻信號

    (1)

    真實目標(biāo)回波信號可以表示為

    (2)

    式中:AR為回波信號幅度;ωd為目標(biāo)速度引起的多普勒頻移。本文的研究基于信號脈內(nèi)分析,因此未考慮目標(biāo)及干擾平臺與雷達(dá)之間距離引起的信號延遲。

    間歇采樣信號為矩形包絡(luò)脈沖串,其脈寬為τ,重復(fù)周期為Ts,包絡(luò)脈沖為

    (3)

    (4)

    其中,任意一個子脈沖信號可以表示為

    (5)

    式中:1≤n≤N,則采樣信號為

    xs(t)=

    (6)

    那么,間歇采樣直接轉(zhuǎn)發(fā)干擾信號可以表示為

    (7)

    間歇采樣重復(fù)轉(zhuǎn)發(fā)干擾信號可以表示為

    (8)

    式中:M為間歇采樣周期Ts內(nèi)可轉(zhuǎn)發(fā)當(dāng)前采樣的次數(shù),且M=?Ts/τ」-1,? 」表示向下取整。

    間歇采樣循環(huán)轉(zhuǎn)發(fā)干擾信號可以表示為

    (9)

    式中:N為脈沖持續(xù)時間t內(nèi)可進行間歇采樣的次數(shù),M為間歇采樣周期Ts內(nèi)可轉(zhuǎn)發(fā)采樣信號的次數(shù),且N=?T/Ts」+1,M=?Ts/τ」-1。

    2熵特征提取

    熵是熱力學(xué)中的一個重要參數(shù),其概念最早由克勞修斯(Clausius)在研究熱力學(xué)第二定律時提出,用來度量物質(zhì)系統(tǒng)中能的分布均勻性。此后,隨著在不同科學(xué)領(lǐng)域的應(yīng)用,熵的概念不斷得到派生和引申,主要用于對不同系統(tǒng)中雜亂無章、不平衡、不確定等無序狀態(tài)進行度量[9-12]。

    2.1信息熵

    1948年,香農(nóng)(Shannon)把熵的概念引入到信息論中,提出了信息熵(又名香農(nóng)熵,Shannonentropy,ShEn),以衡量信源的平均不確定度,解決對信息的量化度量問題。

    (10)

    信息熵的單位隨a的取值而變化,本文在研究過程中,取a=2,此時SH的單位為比特(bit)。

    信息熵應(yīng)用于信號處理,可以作為衡量信號狀態(tài)分布不確定性和信號復(fù)雜程度的特征參數(shù),信號不確定性越大、越復(fù)雜,信息熵就越大。

    2.2指數(shù)熵

    (11)

    2.3范數(shù)熵

    BoekeeD.E.在文獻(xiàn)[14]中提出所謂“R-范數(shù)信息度量”,前提條件與信息熵相同,具體定義為

    (12)

    張葛祥等人在此基礎(chǔ)上,提出了范數(shù)熵(NormEntropy,NoEn)的概念,并用其定量描述信號的能量分布情況。

    (13)

    式中:1

    2.4熵特征提取算法

    研究發(fā)現(xiàn),間歇采樣轉(zhuǎn)發(fā)干擾信號在頻譜上與真實目標(biāo)回波具有明顯差異。干擾信號頻譜的寬度和位置與真實目標(biāo)回波基本相同,不同的是干擾信號頻譜的帶內(nèi)、帶外均出現(xiàn)了許多起伏,帶內(nèi)起伏是各個子脈沖線性調(diào)頻脈沖頻譜的體現(xiàn),帶外起伏是各個子脈沖頻譜相互混疊所致。不同類型間歇采樣轉(zhuǎn)發(fā)干擾的頻譜之間也存在顯著的差異,干擾信號帶內(nèi)頻譜的分布主要受干擾信號參數(shù)設(shè)置影響[15]。

    本文利用干擾信號與真實目標(biāo)回波在頻譜上的差異,提取干擾及回波信號的信息熵、指數(shù)熵以及范數(shù)熵特征,對信號頻譜分布的不確定性進行定量的描述。

    提取熵特征的算法流程如下

    (1) 對3類間歇采樣轉(zhuǎn)發(fā)干擾與目標(biāo)回波信號進行FFT處理,將其從時域變換到頻域;

    (2) 求出信號頻譜的中心頻率和有效帶寬,并對帶寬進行歸一化處理,排除帶外噪聲及載頻的影響;

    (3) 對頻譜進行幅度歸一化處理,并按照范數(shù)熵的定義式,求出4種信號的頻譜范數(shù)熵SR;

    (5) 按照信息熵、指數(shù)熵的定義式,分別求出干擾及目標(biāo)回波的信息熵SH、指數(shù)熵SE。

    3仿真實驗

    為了檢驗3種熵特征的性能,設(shè)置仿真條件如下:雷達(dá)發(fā)射的信號為線性調(diào)頻信號,信號載頻30MBZ,帶寬5MBZ,脈寬100μs,采樣頻率120MBZ;目標(biāo)勻速運動,徑向速度380m/s;間歇采樣直接轉(zhuǎn)發(fā)干擾的間歇采樣周期為10μs,采樣占空比τ/Ts=0.5;重復(fù)轉(zhuǎn)發(fā)干擾的采樣周期為10μs,采樣時間1μs,單個周期轉(zhuǎn)發(fā)9次;循環(huán)轉(zhuǎn)發(fā)干擾的采樣周期為10μs,采樣時間1μs;干信比0dB,信噪比-5~25dB;信息熵參數(shù)取a=2,范數(shù)熵參數(shù)取R=1.5。在不同干信比條件下,對4種信號在每種信噪比下進行300次蒙特卡羅仿真,分別計算4種信號的信息熵、指數(shù)熵及范數(shù)熵特征,得到的特征參數(shù)均值隨信噪比變化情況如圖1~3所示。

    圖1 信息熵均值隨信噪比變化曲線Fig.1 Variation of mean value of ShEn with SNR

    圖2 指數(shù)熵均值隨信噪比變化曲線Fig.2 Variation of mean value of ExEn with SNR

    觀察圖1~3可以發(fā)現(xiàn),熵特征受信噪比影響較小。在信噪比大于0dB的時候,信息熵熵值曲線趨于平穩(wěn),穩(wěn)定后4種信號的熵值曲線沒有交疊,曲線間距離較遠(yuǎn);在信噪比大于0dB的時候,指數(shù)熵熵值曲線同樣達(dá)到穩(wěn)定狀態(tài),4種信號的熵值曲線基本沒有交疊,但曲線之間距離較近,真實目標(biāo)回波與間歇采樣直接轉(zhuǎn)發(fā)干擾的指數(shù)熵值尤為接近;而范數(shù)熵的熵值曲線同樣在信噪比大于0dB后趨于平穩(wěn),穩(wěn)定后4種信號的熵值曲線分離,由于范數(shù)熵值較大,曲線之間距離最遠(yuǎn)。

    圖3 范數(shù)熵均值隨信噪比變化曲線Fig.3 Variation of mean value of NoEn with SNR

    另外,對比圖1~3可以發(fā)現(xiàn),雖然ShEn,ExEn,NoEn的定義不同,但穩(wěn)定后4種信號的熵值大小次序一致,真實目標(biāo)回波的熵值最大,間歇采樣直接轉(zhuǎn)發(fā)干擾次之,間歇采樣重復(fù)轉(zhuǎn)發(fā)干擾的熵值最小,3種熵特征互相進行了印證。

    表1給出了仿真試驗中求取3種熵特征分別所用的時間,本文利用執(zhí)行300次蒙特卡羅實驗求取1 200個熵特征所耗費的時間來表征特征提取過程的時間復(fù)雜度。仿真基于Intel酷睿i3-2120處理器進行,仿真軟件Matlab版本為2014a。由表1可以看出,3種熵特征的時間復(fù)雜度較為接近。從確切的數(shù)值上來看,通過換算可得,計算單個特征的時間需要約70ms。在工程應(yīng)用中,干擾識別的實時性是十分重要的因素。而考慮到間歇采樣轉(zhuǎn)發(fā)干擾的主要應(yīng)用背景要求其干擾信號樣式預(yù)先設(shè)定,在干擾過程中不隨機發(fā)生變化,因此即使在雷達(dá)發(fā)射信號存在脈組跳頻、脈間跳頻等變化時,前期干擾識別的結(jié)果依然可以為雷達(dá)后續(xù)采取針對性抗干擾措施提供有用的參考信息。

    表1 時間復(fù)雜度

    為了考量3種熵特征的穩(wěn)定性,設(shè)定信噪比分別為-5dB和5dB,其余參數(shù)不變,進行300蒙特卡羅仿真,求得的特征樣本的方差如表2所示。從表2中可以看出,指數(shù)熵特征的方差最小,最為穩(wěn)定;范數(shù)熵方差最大,穩(wěn)定性最差。此外,信噪比越高,方差越小,熵特征穩(wěn)定性越好。

    表2 熵特征的方差

    為了驗證熵特征作為干擾識別特征參數(shù)的有效性,選用徑向基核函數(shù)設(shè)計的支持向量機對4種信號進行分類識別,分別選用ShEn,ExEn,NoEn作為特征參數(shù),在干信比JSR=0dB的條件下,對每種信號在信噪比為-5~25dB時分別進行300次蒙特卡羅仿真,得到300個特征樣本。選用其中100個樣本進行分類器訓(xùn)練,200個樣本用作信號識別測試,得到的識別準(zhǔn)確率如圖4~6所示。

    圖4 ShEn作為特征參數(shù)時的識別結(jié)果Fig.4 Recognition results while ShEn acts as parameter

    圖5 ExEn作為特征參數(shù)時的識別結(jié)果Fig.5 Recognition results while ExEn acts as parameter

    圖6 NoEn作為特征參數(shù)時的識別結(jié)果Fig.6 Recognition results while NoEn acts as parameter

    比較圖4~6可知,信息熵作為特征參數(shù)時識別效果最優(yōu),識別準(zhǔn)確率均在85%以上,信噪比高于0dB以后準(zhǔn)確率逼近100%;指數(shù)熵作為特征參數(shù)時識別結(jié)果稍次于信息熵,在低信噪比條件下,識別準(zhǔn)確率略低,隨著信噪比的增大,識別準(zhǔn)確率逐漸增大;范數(shù)熵作為特征參數(shù)時識別效果最差,此時低信噪比條件下識別準(zhǔn)確率與指數(shù)熵相近,隨著信噪比的增大,識別準(zhǔn)確率相應(yīng)有所提高,可達(dá)到95%以上,但在高信噪比時,識別準(zhǔn)確率仍存在明顯波動。

    綜合以上仿真結(jié)果,可以分析總結(jié)得到以下結(jié)論:①信息熵特征類間分離度好、穩(wěn)定性好、識別性能最優(yōu);②指數(shù)熵特征類間分離度稍差,但特征值穩(wěn)定,識別性能較優(yōu);③范數(shù)熵特征類間分離度好,但特征參數(shù)方差大、不穩(wěn)定,識別性能最差。

    4結(jié)束語

    干擾識別歷來是抗干擾的前提和基礎(chǔ),本文針對間歇采樣轉(zhuǎn)發(fā)干擾的識別問題展開研究,提出了基于熵特征的間歇采樣轉(zhuǎn)發(fā)干擾識別方法,將熵的概念引入到特征參數(shù)提取過程中,分別提取信號的信息熵、指數(shù)熵、范數(shù)熵特征,并結(jié)合仿真結(jié)果分析3種熵特征的性能。通過仿真實驗和數(shù)據(jù)分析,熵特征作為干擾識別的特征參數(shù)均能取得較為理想的識別結(jié)果,可以有效地區(qū)分不同干擾及真實目標(biāo)回波,其中信息熵作為特征參數(shù)時識別性能最優(yōu),范數(shù)熵作為特征參數(shù)時識別效果最差。

    參考文獻(xiàn):

    [1]趙鋒, 王雪松, 肖順平,等. 現(xiàn)代雷達(dá)電子戰(zhàn)系統(tǒng)建模與仿真[M]. 北京: 電子工業(yè)出版社, 2010.

    ZHAOFeng,WANGXue-song,XIAOShun-ping,etal.ModelingandSimulationofModernRadarandElectronicWarfareSystems[M].Beijing:PublishingHouseofElectronicIndustry, 2010.

    [2]劉忠, 王雪松, 劉建成,等. 基于數(shù)字射頻存儲器的間歇采樣重復(fù)轉(zhuǎn)發(fā)干擾[J]. 兵工學(xué)報, 2008, 29(4): 405-410.

    LIUZhong,WANGXue-song,LIUJian-cheng,etal.JammingTechniqueofInterrupted-SamplingandPeriodicRepeaterBasedonDigitalRadioFrequencyMemory[J].ActaArmamentarii, 2008, 29(4): 405-410.

    [3]房明星, 王杰貴. 卷積調(diào)制的SAR雷達(dá)二維間歇采樣轉(zhuǎn)發(fā)干擾技術(shù)[J]. 火力與指揮控制, 2014, 39(7): 59-62.

    FANGMing-xin,WANGJie-gui. 2-DIntermittentSamplingRepeaterJammingtoSARBasedonConvolutionModulation[J].FireControl&CommandControl, 2014, 39(7): 59-62.

    [4]潘小義, 王偉, 馮德軍,等. 對解線頻調(diào)ISAR的間歇采樣轉(zhuǎn)發(fā)干擾[J]. 宇航學(xué)報, 2013, 34(9): 1274-1280.

    PANXiao-yi,WANGWei,FENGDe-jun,etal.JammingDechirpingISARBasedonIntermittentSamplingRepeater[J].JournalofAstronautics, 2013, 34(9): 1274-1280.

    [5]李建勛, 唐斌, 呂強. 雙譜特征提取在欺騙式干擾方式識別中的應(yīng)用[J]. 電子科技大學(xué)學(xué)報, 2009, 38(3): 329-332.

    LIJian-xun,TANGBin,LüQiang.BispectrumFeatureExtractionUsedinDeceptiveJammingModesRecognition[J].JournalofUniversityofElectronicScienceandTechnologyofChina, 2009, 38(3): 329-332.

    [6]孫閩紅, 唐斌. 雷達(dá)DRFM欺騙干擾的檢測[J]. 信號處理, 2010, 26(5): 672-676.

    SUNMin-hong,TANGBin.DetectionofRadarDRFMDeceptionJamming[J].SignalProcessing, 2010, 26(5): 672-676.

    [7]蘇曉陽, 郝志梅, 曹蘭英. 基于時頻分析的數(shù)字射頻存儲干擾信號檢測[J]. 火力與指揮控制, 2011, 36(7): 195-196.

    SUXiao-yang,HAOZhi-mei,CAOLan-ying.DetectionofDRFMSignalsBasedonTime-FrequencyAnalysis[J].FireControl&CommandControl, 2011, 36(7): 195-196.

    [8]田曉, 唐斌. 基于歸一化小波分解能量比的雷達(dá)有源欺騙干擾識別[J]. 數(shù)據(jù)采集與處理, 2013, 28(4): 416-420.

    TIANXiao,TANGBin.ActiveDeceptionJammingRecognitionofRadarBasedonNormalizedWaveletDecompositionPowerRatio[J] .JournalofDataAcquisition&Processing, 2013, 28(4): 416-420.

    [9]李靖超, 李一兵, 林云. 熵值分析法在輻射源特征提取中的應(yīng)用[J]. 彈箭與制導(dǎo)學(xué)報, 2011, 31(5): 155-156.

    LIJing-chao,LIYi-bing,LINYun.TheApplicationofEntropyAnalysisinRadiationSourceFeatureExtraction[J].JournalofProjectiles,Rockets,MissilesandGuidance, 2011, 31(5): 155-156.

    [10]吳思東, 朱明, 付克昌. 基于多元集對分析的輻射源信號熵特征評價[J]. 電路與系統(tǒng)學(xué)報, 2013, 18(2): 298-304.

    WUSi-dong,ZHUMing,FUKe-chang.EntropyFeatureEvaluationofRadarEmitterSignalsBasedonSPA[J].JournalofCircuitsandSystems, 2013, 18(2): 298-304.

    [11]張葛祥, 胡來招, 金煒東. 基于熵特征的雷達(dá)輻射源信號識別[J]. 電波科學(xué)學(xué)報, 2005,20(4):440-445.

    ZHANGGe-xiang,HULai-zhao,JINWei-dong.RadarEmitterSignalRecognitionBasedonEntropyFeatures[J].ChineseJournalofRadioScience, 2005, 20(4):440-445.

    [12]李一兵, 常國彬, 葉方. 基于指數(shù)熵的認(rèn)知無線電頻譜感知算法[J]. 吉林大學(xué)學(xué)報:工學(xué)版, 2014, 44(5): 1506-1511.

    LIYi-bing,CHANGGuo-bin,YEFang.ExponentialEntropy-BasedSpectrumSensingAlgorithminCognitiveRadio[J].JournalofJilinUniversity:Eng.andTechnol.ed., 2014, 44(5): 1506-1511.

    [13]SankarKPal,NikhilRPal.Object-backgroundClassificationUsingANewDefinitionofEntropy[M].Proceedingsofthe1988IEEEIntermationalConferenceonSystems,man,andCybemetics,1988,2: 773-776.

    [14]BOEKEEDE,VANDERLubeJCA.TheR-normInformationMeasure[J].InformationandControl, 1980, 45: 136-155

    [15]劉忠. 基于DRFM的線性調(diào)頻脈沖壓縮雷達(dá)干擾新技術(shù)[D].長沙: 國防科學(xué)技術(shù)大學(xué), 2006.

    LIUZhong.JammingTechnologyforCounteringLFMPulseConpressionRaderBasedonDigitalRadioFrequencyMemory[D] .Changsha:NationalUniversityofDefenseTechnology, 2006.

    Novel Method of Interrupted-Sampling Repeater JammingRecognitionBasedonEntropyTheory

    JIANG Ying,HE Ming-hao, YU Chun-lai, WANG Bing-qie

    (Air Force Early Warning Academy,Hubei Wuhan 430019, China)

    Abstract:As jamming recognition is the key link of anti-jamming system, aiming at recognition of interrupted-sampling repeater jamming (ISRJ), the concept of entropy is applied to feature extraction and a jamming recognition scheme based on entropy theory is proposed. In addition, the support vector machine (SVM) is adopted to classify the target and jamming. After transforming signals through FFT transform, Shannon entropy, exponent entropy and norm entropy are extracted. Then, experiment is carried out to verify performance of entropy features. The simulation results show that entropy features are not sensitive to noise and desired recognition results can be obtained while entropy features are used as parameters to discriminate target and jamming. According to the performance of recognition, Shannon entropy is the best, exponent entropy is the next step and norm entropy is bad.

    Key words:interrupted-sampling; jamming recognition; feature extraction; shannon entropy; exponential entropy; norm entropy; support vector machine

    *收稿日期:2015-06-16;修回日期:2015-07-15

    作者簡介:蔣瑩(1991-),女,江蘇南京人。碩士生,研究方向為電子對抗信息處理。

    通信地址:430019湖北省武漢市江岸區(qū)黃浦大街288號研管21隊E-mail:jty614@163.com

    doi:10.3969/j.issn.1009-086x.2016.03.017

    中圖分類號:TN974;TP391.9

    文獻(xiàn)標(biāo)志碼:A

    文章編號:1009-086X(2016)-03-0104-06

    猜你喜歡
    支持向量機信息熵特征提取
    基于信息熵可信度的測試點選擇方法研究
    基于Daubechies(dbN)的飛行器音頻特征提取
    電子制作(2018年19期)2018-11-14 02:37:08
    基于信息熵的實驗教學(xué)量化研究
    電子測試(2017年12期)2017-12-18 06:35:48
    Bagging RCSP腦電特征提取算法
    一種基于信息熵的雷達(dá)動態(tài)自適應(yīng)選擇跟蹤方法
    動態(tài)場景中的視覺目標(biāo)識別方法分析
    論提高裝備故障預(yù)測準(zhǔn)確度的方法途徑
    價值工程(2016年32期)2016-12-20 20:36:43
    基于熵技術(shù)的公共事業(yè)費最優(yōu)組合預(yù)測
    價值工程(2016年29期)2016-11-14 00:13:35
    基于支持向量機的金融數(shù)據(jù)分析研究
    基于信息熵的IITFN多屬性決策方法
    午夜影院日韩av| 18禁裸乳无遮挡免费网站照片| 午夜免费成人在线视频| 欧美不卡视频在线免费观看| 啦啦啦观看免费观看视频高清| 成年女人毛片免费观看观看9| 我要看日韩黄色一级片| 他把我摸到了高潮在线观看| 久久国产乱子伦精品免费另类| 亚洲色图av天堂| 免费一级毛片在线播放高清视频| 久久99热6这里只有精品| 丝袜美腿在线中文| 嫩草影院入口| 亚洲欧美日韩卡通动漫| 成人午夜高清在线视频| 琪琪午夜伦伦电影理论片6080| 听说在线观看完整版免费高清| 精品久久国产蜜桃| 精品午夜福利在线看| 黄色女人牲交| 自拍偷自拍亚洲精品老妇| 波多野结衣巨乳人妻| 欧美极品一区二区三区四区| 成人性生交大片免费视频hd| 欧美成人性av电影在线观看| 亚洲av日韩精品久久久久久密| 国产欧美日韩一区二区三| 村上凉子中文字幕在线| 一进一出抽搐gif免费好疼| 欧美+亚洲+日韩+国产| 国产毛片a区久久久久| 日韩欧美在线乱码| 亚洲av五月六月丁香网| 老司机深夜福利视频在线观看| 久久久精品大字幕| 久久国产精品人妻蜜桃| 欧美成人免费av一区二区三区| x7x7x7水蜜桃| 99riav亚洲国产免费| 欧美日韩乱码在线| 亚洲性夜色夜夜综合| 一进一出好大好爽视频| 亚洲精品成人久久久久久| 神马国产精品三级电影在线观看| 国产伦精品一区二区三区四那| 久久性视频一级片| 最近中文字幕高清免费大全6 | 天堂网av新在线| 久久久国产成人精品二区| 免费看日本二区| 国内久久婷婷六月综合欲色啪| 天堂影院成人在线观看| 免费人成在线观看视频色| 天堂av国产一区二区熟女人妻| 欧美绝顶高潮抽搐喷水| 精品不卡国产一区二区三区| 91午夜精品亚洲一区二区三区 | 精品久久久久久久久亚洲 | 婷婷精品国产亚洲av在线| 亚洲人成网站在线播放欧美日韩| 欧美+日韩+精品| 国产精品电影一区二区三区| 久久精品国产亚洲av天美| xxxwww97欧美| 男人狂女人下面高潮的视频| 91麻豆av在线| 别揉我奶头 嗯啊视频| 高清在线国产一区| 小蜜桃在线观看免费完整版高清| 搞女人的毛片| 小说图片视频综合网站| 成人特级黄色片久久久久久久| 淫妇啪啪啪对白视频| 日韩亚洲欧美综合| 国产一区二区激情短视频| 村上凉子中文字幕在线| 人人妻人人澡欧美一区二区| 国产高潮美女av| 99国产极品粉嫩在线观看| 亚洲av不卡在线观看| 亚洲美女搞黄在线观看 | 国产主播在线观看一区二区| 我的女老师完整版在线观看| 亚洲国产精品sss在线观看| 日韩免费av在线播放| 波野结衣二区三区在线| 韩国av一区二区三区四区| 欧美精品啪啪一区二区三区| 欧美日韩瑟瑟在线播放| 少妇熟女aⅴ在线视频| 一本综合久久免费| 中文资源天堂在线| 91在线观看av| .国产精品久久| 亚洲国产高清在线一区二区三| ponron亚洲| 此物有八面人人有两片| 亚洲综合色惰| 国产三级在线视频| 亚洲乱码一区二区免费版| 色哟哟哟哟哟哟| 午夜福利18| 日韩免费av在线播放| 久久国产乱子伦精品免费另类| 国产主播在线观看一区二区| 国产国拍精品亚洲av在线观看| 久久久久久大精品| 一区二区三区四区激情视频 | 午夜激情福利司机影院| 色视频www国产| 国产精品伦人一区二区| 久久久久精品国产欧美久久久| 别揉我奶头 嗯啊视频| 亚洲第一欧美日韩一区二区三区| 91麻豆av在线| 国产淫片久久久久久久久 | 欧美黄色淫秽网站| 国产麻豆成人av免费视频| 一级av片app| 日本撒尿小便嘘嘘汇集6| bbb黄色大片| 桃红色精品国产亚洲av| 真人一进一出gif抽搐免费| 欧美色视频一区免费| 97超级碰碰碰精品色视频在线观看| 国产成人啪精品午夜网站| bbb黄色大片| 国产亚洲精品久久久com| 一个人观看的视频www高清免费观看| 日韩高清综合在线| 嫩草影院新地址| 女人十人毛片免费观看3o分钟| 深爱激情五月婷婷| 久久久久久久午夜电影| 国产一区二区三区在线臀色熟女| 亚洲av一区综合| 国产高清有码在线观看视频| 欧美日韩乱码在线| 制服丝袜大香蕉在线| 69av精品久久久久久| 毛片一级片免费看久久久久 | 亚洲国产高清在线一区二区三| 久久香蕉精品热| 日本黄色视频三级网站网址| 午夜福利在线观看免费完整高清在 | 精品免费久久久久久久清纯| 在线观看美女被高潮喷水网站 | 麻豆成人午夜福利视频| 久久久国产成人免费| bbb黄色大片| 99热6这里只有精品| 国产单亲对白刺激| 老司机午夜福利在线观看视频| 人妻久久中文字幕网| 亚洲欧美精品综合久久99| 亚洲国产精品久久男人天堂| 高清在线国产一区| 欧洲精品卡2卡3卡4卡5卡区| 两人在一起打扑克的视频| 成人鲁丝片一二三区免费| 1024手机看黄色片| 中文字幕精品亚洲无线码一区| 欧美中文日本在线观看视频| 国产单亲对白刺激| 一卡2卡三卡四卡精品乱码亚洲| 精品99又大又爽又粗少妇毛片 | 狂野欧美白嫩少妇大欣赏| 亚洲欧美日韩高清专用| 在线看三级毛片| 最新在线观看一区二区三区| 免费搜索国产男女视频| 国产精品精品国产色婷婷| 亚洲狠狠婷婷综合久久图片| 久久精品影院6| 我要搜黄色片| 在线观看美女被高潮喷水网站 | 91麻豆av在线| 又黄又爽又刺激的免费视频.| 欧美日本视频| 91狼人影院| 欧美精品啪啪一区二区三区| 在线a可以看的网站| 91久久精品电影网| av天堂中文字幕网| 亚洲 欧美 日韩 在线 免费| 国产三级中文精品| 3wmmmm亚洲av在线观看| 中文字幕熟女人妻在线| 亚洲五月婷婷丁香| 欧美丝袜亚洲另类 | 久久久色成人| 狠狠狠狠99中文字幕| 最近视频中文字幕2019在线8| 亚洲av一区综合| 夜夜爽天天搞| 一本久久中文字幕| 美女cb高潮喷水在线观看| 夜夜爽天天搞| 在线看三级毛片| 丰满人妻熟妇乱又伦精品不卡| 高潮久久久久久久久久久不卡| 欧美高清成人免费视频www| 人妻制服诱惑在线中文字幕| 人妻夜夜爽99麻豆av| 99国产极品粉嫩在线观看| 久久久成人免费电影| 国产一区二区激情短视频| 搡老妇女老女人老熟妇| 俺也久久电影网| 亚洲一区高清亚洲精品| 国产高清三级在线| 亚洲中文日韩欧美视频| 国产欧美日韩精品一区二区| 91狼人影院| 女同久久另类99精品国产91| 三级国产精品欧美在线观看| 成人三级黄色视频| 超碰av人人做人人爽久久| 18禁黄网站禁片午夜丰满| 丁香六月欧美| av在线观看视频网站免费| 少妇裸体淫交视频免费看高清| 中亚洲国语对白在线视频| 亚洲精品456在线播放app | 脱女人内裤的视频| 91在线观看av| 国产精品亚洲一级av第二区| 五月伊人婷婷丁香| 中文字幕熟女人妻在线| 国产真实乱freesex| 国产精品一区二区免费欧美| 亚洲经典国产精华液单 | 美女高潮的动态| 久久久久久久久中文| 成人毛片a级毛片在线播放| 欧美黑人欧美精品刺激| 高清毛片免费观看视频网站| 别揉我奶头 嗯啊视频| 欧美日韩瑟瑟在线播放| 日韩欧美精品v在线| 亚洲最大成人av| 亚洲成人久久爱视频| 精品国内亚洲2022精品成人| 欧美区成人在线视频| av欧美777| 69av精品久久久久久| 中出人妻视频一区二区| 欧美乱色亚洲激情| 欧美日本亚洲视频在线播放| 色哟哟·www| 国产淫片久久久久久久久 | 我的老师免费观看完整版| 五月玫瑰六月丁香| 男插女下体视频免费在线播放| 又爽又黄a免费视频| 婷婷精品国产亚洲av| 3wmmmm亚洲av在线观看| 搡老妇女老女人老熟妇| 又爽又黄无遮挡网站| 亚洲,欧美,日韩| 婷婷色综合大香蕉| 男人的好看免费观看在线视频| 欧美潮喷喷水| 夜夜看夜夜爽夜夜摸| 精品乱码久久久久久99久播| 最近最新免费中文字幕在线| 国产久久久一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 国产真实伦视频高清在线观看 | 日本a在线网址| 精品日产1卡2卡| 欧美激情久久久久久爽电影| 午夜福利欧美成人| 亚洲国产精品成人综合色| 蜜桃久久精品国产亚洲av| h日本视频在线播放| 精品国产三级普通话版| 中亚洲国语对白在线视频| 日本在线视频免费播放| 中文亚洲av片在线观看爽| 精品乱码久久久久久99久播| 毛片女人毛片| 亚洲av第一区精品v没综合| 国产亚洲精品久久久com| 日韩欧美国产在线观看| 精品不卡国产一区二区三区| 色吧在线观看| 精品欧美国产一区二区三| 精品国内亚洲2022精品成人| 日韩免费av在线播放| 亚洲最大成人av| 色视频www国产| 亚洲va日本ⅴa欧美va伊人久久| av国产免费在线观看| 99国产精品一区二区蜜桃av| 成人亚洲精品av一区二区| 国产黄片美女视频| 亚洲激情在线av| 亚洲成av人片免费观看| 可以在线观看毛片的网站| 久久久久亚洲av毛片大全| 欧美日韩乱码在线| 欧美日韩亚洲国产一区二区在线观看| 99久久九九国产精品国产免费| 午夜a级毛片| 毛片一级片免费看久久久久 | 午夜精品在线福利| av在线老鸭窝| 一区二区三区高清视频在线| 午夜福利免费观看在线| 欧美黄色淫秽网站| 成人国产综合亚洲| 国产高清激情床上av| 午夜免费激情av| 窝窝影院91人妻| 2021天堂中文幕一二区在线观| 天美传媒精品一区二区| 内地一区二区视频在线| 欧美日韩综合久久久久久 | 97碰自拍视频| 亚洲狠狠婷婷综合久久图片| 国产精品伦人一区二区| 黄色配什么色好看| 国产欧美日韩一区二区三| 精品欧美国产一区二区三| 国产在线男女| 日本与韩国留学比较| 欧美成人a在线观看| 久久精品国产99精品国产亚洲性色| 91在线观看av| 婷婷色综合大香蕉| 午夜老司机福利剧场| 欧美一级a爱片免费观看看| 日韩人妻高清精品专区| 一级毛片久久久久久久久女| 一卡2卡三卡四卡精品乱码亚洲| av天堂中文字幕网| 超碰av人人做人人爽久久| 久久久久精品国产欧美久久久| 国产精品久久久久久亚洲av鲁大| 99国产精品一区二区蜜桃av| 99热精品在线国产| 少妇被粗大猛烈的视频| 欧美午夜高清在线| 日本免费a在线| 日韩有码中文字幕| 乱码一卡2卡4卡精品| 久久久久国产精品人妻aⅴ院| 国产乱人视频| 97超视频在线观看视频| 欧美在线一区亚洲| 桃色一区二区三区在线观看| 脱女人内裤的视频| 九九在线视频观看精品| 色噜噜av男人的天堂激情| av在线天堂中文字幕| 麻豆成人午夜福利视频| 中文字幕人妻熟人妻熟丝袜美| 伊人久久精品亚洲午夜| 3wmmmm亚洲av在线观看| 99国产精品一区二区蜜桃av| 中文字幕免费在线视频6| 人妻制服诱惑在线中文字幕| 3wmmmm亚洲av在线观看| 一个人免费在线观看电影| 久久精品国产99精品国产亚洲性色| 国产精品99久久久久久久久| 少妇人妻一区二区三区视频| 一夜夜www| 成熟少妇高潮喷水视频| 黄色配什么色好看| 国产主播在线观看一区二区| 精品久久久久久久末码| 日韩欧美国产在线观看| 一区二区三区激情视频| 波多野结衣高清作品| 简卡轻食公司| 在线国产一区二区在线| av在线观看视频网站免费| 在线观看美女被高潮喷水网站 | 看免费av毛片| 大型黄色视频在线免费观看| 噜噜噜噜噜久久久久久91| 久久久国产成人精品二区| 人妻制服诱惑在线中文字幕| 亚洲美女视频黄频| 一本久久中文字幕| 亚洲三级黄色毛片| 日韩免费av在线播放| 舔av片在线| 如何舔出高潮| 亚洲专区国产一区二区| 97超视频在线观看视频| 一级作爱视频免费观看| 99久久99久久久精品蜜桃| 亚洲国产精品sss在线观看| 国产精品久久久久久人妻精品电影| 久久伊人香网站| 国产午夜福利久久久久久| 国产精品三级大全| 欧美黄色淫秽网站| 国产精品一区二区三区四区免费观看 | 欧美一区二区亚洲| 亚洲avbb在线观看| 精品不卡国产一区二区三区| АⅤ资源中文在线天堂| x7x7x7水蜜桃| 欧美日韩瑟瑟在线播放| 香蕉av资源在线| 亚洲人成电影免费在线| 91九色精品人成在线观看| 久久精品国产自在天天线| 久久久久九九精品影院| a级一级毛片免费在线观看| 看黄色毛片网站| 亚洲精品日韩av片在线观看| 亚洲五月天丁香| 人妻丰满熟妇av一区二区三区| 午夜福利在线在线| 丝袜美腿在线中文| 国产高清三级在线| 日韩精品青青久久久久久| 狠狠狠狠99中文字幕| 中文在线观看免费www的网站| 一区二区三区高清视频在线| 欧美一区二区亚洲| 国内久久婷婷六月综合欲色啪| av在线蜜桃| 偷拍熟女少妇极品色| 国内久久婷婷六月综合欲色啪| 午夜福利在线观看免费完整高清在 | 成人av一区二区三区在线看| 欧美在线一区亚洲| 欧美中文日本在线观看视频| 波野结衣二区三区在线| 久久久久久久午夜电影| 亚洲人成网站在线播| 91久久精品电影网| 99国产精品一区二区蜜桃av| 人妻丰满熟妇av一区二区三区| 日韩欧美精品v在线| 午夜福利视频1000在线观看| 好看av亚洲va欧美ⅴa在| 亚洲av.av天堂| 99热这里只有是精品在线观看 | 性插视频无遮挡在线免费观看| 成人欧美大片| 久久亚洲精品不卡| 狂野欧美白嫩少妇大欣赏| 一个人看的www免费观看视频| 国产一区二区亚洲精品在线观看| 精品99又大又爽又粗少妇毛片 | 男女下面进入的视频免费午夜| 一个人免费在线观看电影| 亚洲在线观看片| 亚洲精品成人久久久久久| 成年人黄色毛片网站| 国产免费男女视频| 亚洲国产欧洲综合997久久,| 欧美成人免费av一区二区三区| 亚洲成人久久性| 99久久无色码亚洲精品果冻| 国产欧美日韩精品一区二区| 婷婷六月久久综合丁香| 一级av片app| 久久国产乱子免费精品| 欧美中文日本在线观看视频| 日韩成人在线观看一区二区三区| 99热这里只有是精品在线观看 | 性色avwww在线观看| 午夜日韩欧美国产| 九九久久精品国产亚洲av麻豆| 欧美日韩福利视频一区二区| 亚洲最大成人中文| 一二三四社区在线视频社区8| 国产单亲对白刺激| 午夜福利在线观看免费完整高清在 | 此物有八面人人有两片| 毛片一级片免费看久久久久 | 简卡轻食公司| 国产精品一区二区三区四区免费观看 | 久久久久国内视频| 91在线观看av| 黄片小视频在线播放| 欧美黄色淫秽网站| 国产极品精品免费视频能看的| 免费大片18禁| 亚洲欧美日韩高清在线视频| 久久精品国产亚洲av香蕉五月| 又爽又黄无遮挡网站| 全区人妻精品视频| 听说在线观看完整版免费高清| 欧美国产日韩亚洲一区| 中文字幕av成人在线电影| 亚洲黑人精品在线| 欧美最黄视频在线播放免费| av国产免费在线观看| 男人和女人高潮做爰伦理| 国产精品99久久久久久久久| 亚洲,欧美精品.| 欧美色欧美亚洲另类二区| 欧美日韩综合久久久久久 | 国产一区二区三区视频了| 午夜福利在线观看吧| 久久久久国产精品人妻aⅴ院| 亚洲七黄色美女视频| 国产午夜精品论理片| 嫩草影院精品99| 国产精品不卡视频一区二区 | 欧美绝顶高潮抽搐喷水| 一级毛片久久久久久久久女| 欧美日韩国产亚洲二区| 久久久精品欧美日韩精品| 天堂网av新在线| 97人妻精品一区二区三区麻豆| 国产精品嫩草影院av在线观看 | 国产欧美日韩一区二区精品| 国产精品乱码一区二三区的特点| 久久精品影院6| 免费黄网站久久成人精品 | 99riav亚洲国产免费| 亚洲无线观看免费| 一二三四社区在线视频社区8| 人妻久久中文字幕网| 又粗又爽又猛毛片免费看| 欧美zozozo另类| 看片在线看免费视频| 亚洲av.av天堂| 女人被狂操c到高潮| 悠悠久久av| 国产精品久久视频播放| 韩国av一区二区三区四区| 亚洲自偷自拍三级| 亚洲男人的天堂狠狠| 精品人妻1区二区| 在现免费观看毛片| 欧美激情久久久久久爽电影| 亚洲 欧美 日韩 在线 免费| а√天堂www在线а√下载| 国产伦精品一区二区三区四那| 色噜噜av男人的天堂激情| 国产中年淑女户外野战色| 亚洲成人久久性| 色精品久久人妻99蜜桃| 国产精品国产高清国产av| 男女之事视频高清在线观看| 在线观看舔阴道视频| 90打野战视频偷拍视频| 国产精品野战在线观看| 国产欧美日韩精品亚洲av| 久久久久久久亚洲中文字幕 | 亚洲男人的天堂狠狠| 婷婷色综合大香蕉| 国产精品伦人一区二区| 亚洲乱码一区二区免费版| 无人区码免费观看不卡| 久久热精品热| 精品熟女少妇八av免费久了| 免费在线观看影片大全网站| 久久午夜亚洲精品久久| 成年女人看的毛片在线观看| 一本久久中文字幕| 中文亚洲av片在线观看爽| 性色avwww在线观看| 精品不卡国产一区二区三区| 久久天躁狠狠躁夜夜2o2o| 黄色日韩在线| 色播亚洲综合网| 黄色配什么色好看| 两个人视频免费观看高清| 在线十欧美十亚洲十日本专区| 一夜夜www| 人妻制服诱惑在线中文字幕| 国产高清有码在线观看视频| 人妻制服诱惑在线中文字幕| av视频在线观看入口| 成人国产综合亚洲| 免费人成视频x8x8入口观看| 伦理电影大哥的女人| 日本 欧美在线| 精品不卡国产一区二区三区| 性色avwww在线观看| 精品国产亚洲在线| 久久久色成人| 成人毛片a级毛片在线播放| 可以在线观看的亚洲视频| 大型黄色视频在线免费观看| 美女 人体艺术 gogo| av在线观看视频网站免费| 亚洲一区高清亚洲精品| 97人妻精品一区二区三区麻豆| 极品教师在线视频| 午夜a级毛片| 级片在线观看| 深夜精品福利| 成人一区二区视频在线观看| 国产精品日韩av在线免费观看| 色综合婷婷激情| 亚洲av一区综合| 国产精品国产高清国产av| 精品久久久久久久末码| 成人特级黄色片久久久久久久| 在现免费观看毛片| 欧美一区二区精品小视频在线| 午夜激情福利司机影院| 久久人人精品亚洲av| 亚洲avbb在线观看| 国产精品一区二区三区四区久久| 成人一区二区视频在线观看| 九九热线精品视视频播放| 好看av亚洲va欧美ⅴa在| 女人十人毛片免费观看3o分钟| 亚洲专区国产一区二区| 久久性视频一级片|