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

    一種快速抗間歇采樣轉(zhuǎn)發(fā)干擾波形和濾波器聯(lián)合設計算法

    2022-04-30 02:00:34周凱何峰粟毅
    雷達學報 2022年2期
    關鍵詞:干擾信號間歇濾波器

    周凱 何峰 粟毅

    (國防科技大學電子科學學院 長沙 410073)

    1 引言

    隨著射頻存儲轉(zhuǎn)發(fā)技術(shù)的發(fā)展,基于采樣轉(zhuǎn)發(fā)的相參干擾能夠?qū)走_同時形成壓制和欺騙干擾效果[1–4]。間歇采樣轉(zhuǎn)發(fā)干擾是典型的基于存儲轉(zhuǎn)發(fā)的相干干擾,其給雷達系統(tǒng)帶來嚴峻挑戰(zhàn)[5–9]。為了提高雷達的抗干擾工作性能,抗間歇采樣轉(zhuǎn)發(fā)干擾技術(shù)已經(jīng)成為目前的研究熱點。

    目前抗間歇采樣轉(zhuǎn)發(fā)干擾技術(shù)包括接收端信號處理、發(fā)射端波形設計和聯(lián)合波形與濾波器設計方法。時頻分析是抑制間歇采樣轉(zhuǎn)發(fā)干擾的主要接收端信號處理方法[10–12]。文獻[10]根據(jù)干擾信號和目標回波之間的時頻差異實現(xiàn)了對3類典型間歇采樣轉(zhuǎn)發(fā)干擾的辨識。在此基礎上,其提出了一種基于時頻域濾波器設計的干擾抑制方法。文獻[11]則重點研究了間歇采樣轉(zhuǎn)發(fā)干擾的識別方法。文獻[13]估計了間歇采樣轉(zhuǎn)發(fā)干擾的關鍵參數(shù),然后通過干擾信號重構(gòu)有效抑制了間歇采樣轉(zhuǎn)發(fā)干擾。然而,現(xiàn)有基于時頻處理的方法要求干擾信號和雷達回波在時頻域可分辨[11]。當干擾信號產(chǎn)生的假目標和真實目標之間距離不可分時,此類方法的抗干擾性能將顯著下降。

    第2類方法主要通過設計稀疏多普勒特性波形或脈內(nèi)正交波形,從而達到抑制間歇采樣轉(zhuǎn)發(fā)干擾的目的[14–16]。然而,上述方法僅僅考慮了發(fā)射端的波形設計,而未在開展波形設計時考慮接收端的濾波器的設計問題。上述方法往往需要在接收端通過分段濾波和滑窗抽取檢測等信號處理方法進行干擾識別和抑制。因此,文獻[17]提出了一種聯(lián)合波形和濾波器設計抑制方法,能夠在保證脈沖壓縮低旁瓣性能的同時有效抑制間歇采樣轉(zhuǎn)發(fā)干擾。但是其優(yōu)化算法存在脈沖壓縮峰值增益不可控、計算復雜度高和運行時間長等問題。

    文獻[18]提出一種改進的抗間歇采樣轉(zhuǎn)發(fā)干擾模型,能有效控制脈沖壓縮峰值增益。然而,目標函數(shù)近似導致算法收斂速度較慢,且矩陣的乘法計算導致了計算復雜度較高。

    為解決上述問題,本文提出一種基于加速主分量最小化的抗間歇采樣轉(zhuǎn)發(fā)干擾聯(lián)合波形和濾波器設計算法。本文首先基于罰函數(shù)和帕累托最優(yōu)化原理建立了抗間歇采樣轉(zhuǎn)發(fā)干擾聯(lián)合設計的數(shù)學優(yōu)化模型。在此基礎上,基于主分量最小化原理推導了迭代過程中波形和濾波器的解析表達式,并通過推導求解過程中矩陣跡的解析表達式,有效避免了矩陣乘法運算,降低了算法計算復雜度。同時,提出了一種基于平方迭代加速方法,進一步加快了算法運行速度。最后,仿真對抗間歇采樣轉(zhuǎn)發(fā)干擾聯(lián)合設計數(shù)學模型關鍵參數(shù)進行了分析,驗證了本文所提算法的可行性。仿真結(jié)果表明,在設計波形和濾波器性能基本一致的情況下,本文所提算法比傳統(tǒng)方法具有更快的運行速度。并且,本文所提出的方法可解決傳統(tǒng)時頻處理方法干擾假目標和真實目標距離不可分辨時的抗干擾性能顯著惡化問題。

    2 聯(lián)合設計問題模型

    本節(jié)闡述抗間歇采樣轉(zhuǎn)發(fā)干擾聯(lián)合波形和非匹配濾波器設計的問題模型,然后推導得出聯(lián)合設計的數(shù)學模型。

    2.1 低脈沖壓縮積分旁瓣代價函數(shù)

    雷達發(fā)射波形為相位編碼信號,其數(shù)學表達式為x=[x1x2...xN]T,其中xn=,?n ∈(0,2π]。非匹配濾波器為,其中N表示雷達波形和非匹配濾波器長度。為建立波形和非匹配濾波卷積的矩陣表達式,定義矩陣[17–19]:

    則積分旁瓣能量為xHΛHW Λx或者hHΦHW Φh。其中,矩陣W表示維度為2N–1的對角方陣,除了第N個對角元素等于0之外,其余對角元素為1,且對角元素之外的元素均為0。此外,為約束脈沖壓縮峰值,引入目標函數(shù)其中bmax表示預設脈沖壓縮峰值。基于罰函數(shù)方法,抑制積分旁瓣的代價函數(shù)為

    2.2 抑制間歇采樣轉(zhuǎn)發(fā)干擾代價函數(shù)

    間歇采樣轉(zhuǎn)發(fā)干擾是利用數(shù)字射頻存儲技術(shù)和數(shù)字信號合成技術(shù)對雷達發(fā)射信號通過間斷采樣轉(zhuǎn)發(fā),從而對雷達實現(xiàn)干擾。對大時間帶寬積脈沖雷達而言,截獲到脈沖信號后,高保真度地采樣其中一小段信號進行轉(zhuǎn)發(fā),再采樣下一段并進行轉(zhuǎn)發(fā),如此收發(fā)分時、采樣和轉(zhuǎn)發(fā)交替工作直至大時寬信號結(jié)束。間歇采樣信號為矩形包絡脈沖串:

    其中,τ表示間歇采樣脈沖寬度,Ts表示間歇采樣重復周期。間歇采樣轉(zhuǎn)發(fā)干擾信號為采樣信號和雷達發(fā)射信號的乘積,即xjam(t)=x(t)p(t),其中x(t)表示雷達發(fā)射信號,xjam(t)表示雷達間歇采樣轉(zhuǎn)發(fā)干擾信號。間歇采樣轉(zhuǎn)發(fā)干擾則可以表示為xjam=x ⊙p。其中,p=[p1,p2,...,pN]T表示間歇采樣信號的離散序列,⊙表示哈達瑪積。同理,構(gòu)建矩陣

    其中,懲罰參數(shù)λ2的目的與式(3)類似,其目的是控制目標函數(shù)的數(shù)值,進而控制設計發(fā)射波形對應的間歇采樣轉(zhuǎn)發(fā)干擾信號非匹配濾波器輸出峰值(h ⊙p)Hx。因此,為抑制間歇采樣轉(zhuǎn)發(fā)干擾峰值,懲罰參數(shù)λ2應該取較大的數(shù)值[20,21]。利用帕累托優(yōu)化框架,同時優(yōu)化上述兩個目標函數(shù),則波形恒模和非匹配濾波器能量約束下的聯(lián)合設計優(yōu)化問題為

    其中,ε表示帕累托權(quán)值。由于式(8)約束了接收濾波器能量和波形能量,且最小化目標函數(shù)f1(x,h)可實現(xiàn)對波形脈沖壓縮峰值hHx的約束,所以本文可有效控制雷達波形的脈沖壓縮信號處理增益損耗。此外,占空比決定了算法中設計波形對應的干擾片段的長度,最小化目標函數(shù)意味著設計波形中與干擾信號對應的多段信號經(jīng)過非匹配濾波后將無法獲得信號處理增益,進而會導致信號非匹配處理增益損耗。因此,間歇采樣轉(zhuǎn)發(fā)干擾的占空比越大,則算法設計得到的信號處理增益損耗將會更大,反之亦然。因此,預設信號處理增益峰值bmax時,應充分考慮占空比數(shù)值。占空比越大,則預設的bmax應設置較小,以確保能夠獲得較好的脈沖壓縮和抗干擾性能。

    3 基于主分量最小化的交替聯(lián)合設計算法

    本節(jié)采用基于主分量最小化原理方法分別求取波形和濾波器的最優(yōu)解,并進一步提出一種加速的交替優(yōu)化設計方法。

    3.1 非匹配濾波器優(yōu)化

    問題(9)是二次約束二次規(guī)劃問題,由于約束條件是非凸集合,因此難以直接求解。為了求解該問題,先利用主分量最小化原理簡化問題。根據(jù)主分量最小化方法原理,假設z(l)是第l次迭代得到的最優(yōu)解,求目標函數(shù)的主分量函數(shù)u(z,z(l)),使其對任意z∈?滿足通過優(yōu)化主分量函數(shù),可得第l+1次的最優(yōu)解為z(l+1)=進而可得g(z(l+1))≤g(z(l))。相關證明過程可參考文獻[21]。

    引理1若對任意z ∈?滿足,則優(yōu)化問題

    可通過迭代優(yōu)化其主分量最小化問題獲得原非凸問題的最優(yōu)解:

    進一步,根據(jù)文獻[20]中引理1和文獻[21]引理1,對任意半正定矩陣L,均有L ?tr(L)IN。在給定z0情況下,zHLz的主分量最小化函數(shù)為

    由于文中R(l)滿足半正定特性,故hHR(l)h的主分量函數(shù)為

    因此,忽略無關常數(shù)項和系數(shù)后,問題(9)可簡化為

    由于式(15)為二次約束的線性規(guī)劃問題,易得非匹配濾波器表達式為

    其中

    3.2 波形優(yōu)化

    然后,給定非匹配濾波器h(l+1),將和f2的表達式代入式(8),波形設計優(yōu)化問題為

    其中,hjam=h ⊙p。因為矩陣T(l+1)為半正定矩陣,所以xHT(l+1)x的主分量函數(shù)為

    根據(jù)引理1,問題(18)可以簡化為

    式(21)為恒模約束線性規(guī)劃問題,很容易得出其解析表達式為

    其中

    3.3 算法收斂分析

    利用主分量最小化方法優(yōu)化非凸問題的收斂性分析已經(jīng)有較多的參考文獻[20,22–27]。歸納而言,為了證明算法的收斂性能,通常需要證明算法滿足以下3個條件:

    (1)算法單調(diào)遞減;

    (3) 目標函數(shù)的Kurdyka-Lojasiewicz特性[25]。

    上述3個條件中,前兩個條件是為了保證設計過程中得到的每個子序列均是目標函數(shù)的駐點。而第3個條件則是為了證明算法生成的序列是柯西序列,進而證明算法最終收斂。

    因此可以證明優(yōu)化問題中的目標函數(shù)是單調(diào)遞減的。第二,相對誤差條件的證明過程是一個標準推導過程,詳細的推導可以參考文獻[20]的附錄B。第三,由于文中的目標函數(shù)f(x,h)是可轉(zhuǎn)換為關于波形、濾波器實部和虛部的實函數(shù),因此,其相應的目標函數(shù)也可轉(zhuǎn)換為實函數(shù)。而另外,所有的實函數(shù)均是滿足Kurdyka-Lojasiewicz特性的[25]。且由于算法滿足單調(diào)遞減、相對誤差條件,故本文所提算法收斂性能得以證明。

    3.4 計算復雜度分析

    根據(jù)式(17)和式(23)可知,更新波形和濾波器的計算復雜度主要集中在其第1項矩陣向量乘積和第2項矩陣跡的計算表達式。其中,將矩陣R(l)表達式代入R(l)h(l)可知,本文式(17)的計算主要涉及矩陣和向量的乘法,因而可避免矩陣與矩陣之間的乘法運算。同理,式(23)的第1項計算復雜度同樣也僅來自矩陣和向量之間的乘法。因此,計(算本)文式(17)和式(23)第1項的計算復雜度均為O此外,若直接計算式(17)和式(23)第2項的跡,則需要完全計算出矩陣R(l)和T(l+1)的表達式,其中涉及矩陣之間的乘法,其計算復雜度為O(N3)。因此,若直接計算矩陣R(l)和T(l+1),算法的計算復雜度將為為了避免計算矩陣R(l)和T(l+1)的跡引入的復雜度,本文進一步推導上述矩陣R(l)和T(l+1)的解析表達式。首先提出下述引理。

    引理2由于波形和濾波器均滿足能量約束,因此矩陣R(l)和T(l+1)的跡滿足

    證明:首先,根據(jù)矩陣R(l)的定義和矩陣跡的計算引理tr(A+B)=tr(A)+tr(B)和tr(AB)=tr(BA),可得

    分別合并式(33)的前兩項和后兩項即可以得到式(26)所述的表達式。矩陣T(l+1)前3項的跡推導過程與矩陣R(l)的跡計算過程一樣,且計算結(jié)果相等,此處不再贅述。其最后一項中,由于下式滿足

    結(jié)合式(32)和上文推導矩陣R(l)的結(jié)果,則容易得到式(2 7)的表達式。定義兩個變量α=εN(N+λ1?1)+(1?ε)(τ/Ts)N(N+1)和β=εN·(N+λ1?1)+(1?ε)(τ/Ts)N+(1?ε)N2。由于R(l)和T(l+1)是半正定矩陣,故可得到下述關系式R(l)?αIN和T(l+1)?βIN。因此,根據(jù)主分量最小化原理,利用α和β分別替換矩陣R(l)和T(l+1)的跡仍可保證算法迭代過程中的目標函數(shù)是單調(diào)下降的。

    通過推導矩陣R(l)和T(l+1)的跡,有效避免了計算上述兩個矩陣引起的矩陣乘法運算。因此,本文所提算法的計算復雜度僅源于式(17)和式(23)第1項的矩陣和向量乘法運算。因此,算法的計算復雜度為O(N2)。

    3.5 算法加速

    在主分量最小化函數(shù)方法中,由于使用了主分量函數(shù)近似等價目標函數(shù),目標函數(shù)收斂速度將會大幅度下降。文獻[28]提出了平方迭代方法,其通過聯(lián)合利用最速下降法和兩點步長梯度方法可以實現(xiàn)對EM (Expectation-Maximization)算法加速。EM算法基本思路是先根據(jù)第l次迭代結(jié)果求取聯(lián)合條件概率期望的極大似然函數(shù),然后根據(jù)上述函數(shù)更新第l+1次的優(yōu)化解。另外,主分量最小化方法首先根據(jù)第l次迭代結(jié)果求目標函數(shù)的主分量最大化函數(shù),然后基于此求第l+1次迭代解析表達式[28]。因此,主分量最大化算法和EM算法工作流程基本是一致的。受此啟發(fā),本文將平方迭代方法應用于本文基于主分量最小化的聯(lián)合波形和濾波器設計方法。除此之外,平方迭代方法僅要求EM算法的更新策略,即波形和濾波器迭代更新表達式,故可直接在算法中進行應用。

    分別利用α和β替換式(16)和式(22)中的跡表達式,定義如下兩個映射函數(shù):

    利用平方迭代加速方法對算法進行加速,抗間歇采樣轉(zhuǎn)發(fā)干擾聯(lián)合優(yōu)化設計算法如表1所示。在本文所提方法中,分別對波形和濾波器迭代過程中應用平方迭代加速方法。其中,濾波器的更新加速對應表1的步驟4—步驟12,而波形的更新加速對應表1的步驟14—步驟22。其中,平方迭代加速方法的基本原理已在文獻[22,29]中詳細介紹,此處不再贅述。由于在加速過程并未引入額外的矩陣乘法等運算,故加速算法的計算復雜度仍為O(N2)。

    表1 快速抗間歇采樣轉(zhuǎn)發(fā)干擾波形和濾波器聯(lián)合設計算法Tab.1 Fast algorithm for joint waveform and filter design against interrupted sampling repeater jamming

    4 仿真

    為了驗證本文所提算法的有效性,本節(jié)開展7個仿真實驗。首先對聯(lián)合設計模型中關鍵參數(shù)進行仿真分析。通過對比傳統(tǒng)方法,分析了本文所提方法運行速度的性能。同時,為了評估聯(lián)合設計的抗干擾性能,將本文所提方法與傳統(tǒng)抗間歇采樣轉(zhuǎn)發(fā)干擾進行了對比分析。最后,對本文所提方法的干擾參數(shù)敏感性進行了分析。所有的仿真實驗均在同一個硬件和軟件平臺開展(PC with 4.20 GHz Intel Core i7-7700K CPU and 64GB RAM using Windows 7)。

    4.1 算法性能分析

    仿真1:假設波形和非匹配濾波器長度為512,波形非匹配濾波輸出歸一化峰值增益損耗為–1 dB,干擾信號非匹配濾波輸出歸一化峰值為–30 dB。算法截止運行閾值ζ=1e?7。間歇采樣信號占空比0.25,間歇采樣周期重復Ts=T/2。首先分析帕累托權(quán)值對算法性能的影響,帕累托權(quán)值ε=[0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9]。不同帕累托權(quán)值條件下的波形和非匹配濾波器初始序列分別為同一隨機相位編碼序列。

    圖1為波形非匹配脈沖壓縮積分旁瓣比和峰值損耗,以及干擾信號非匹配濾波輸出積分能量和峰值隨著帕累托權(quán)值的變化曲線圖。如圖1(a)和圖1(c)所示,隨著帕累托權(quán)值的增加,積分旁瓣比呈現(xiàn)下降趨勢,而峰值損耗值逐漸增加。由目標函數(shù)的定義可知,當帕累托權(quán)值增加時,目標函數(shù)中提高脈沖壓縮性能的權(quán)重增加。因此,在優(yōu)化過程中能夠使得目標函數(shù)f1(x,h)的數(shù)值更小,即積分旁瓣和峰值損耗目標函數(shù)的數(shù)值將會逐漸減小。其中,脈壓峰值目標函數(shù)減小會使得設計波形和濾波器的脈壓峰值和預期設定值之間的誤差值逐漸降低,所以設計波形和濾波器的脈壓峰值將逐漸靠近預期設定值,進而符合逐漸上升趨勢。此外,權(quán)值為0.3和0.8時積分旁瓣比水平值分別要小于權(quán)值為0.4和0.9時對應的數(shù)值,這可能與算法的迭代次數(shù)相關,即權(quán)值為0.3和0.8時迭代次數(shù)少于權(quán)值為0.4和0.9時對應的迭代次數(shù),進而導致其積分旁瓣值偏大。同時,圖1(b)和圖1(d)中干擾信號的積分能量比、干擾峰值隨著帕累托權(quán)值的增加,干擾信號非匹配濾波輸出峰值干擾增加。其原因是帕累托權(quán)值增加使得抑制干擾的代價函數(shù)權(quán)重降低,進而導致干擾信號非匹配濾波輸出峰值增加。

    由圖1仿真分析可知,隨著帕累托權(quán)值的增加,雷達波形脈沖壓縮性能與抗間歇采樣轉(zhuǎn)發(fā)干擾的性能變化趨勢是相反的。也就是說,很難同時獲得高脈沖壓縮性能和抗間歇采樣轉(zhuǎn)發(fā)干擾性能。實際上,參數(shù)的選擇可以根據(jù)場景中間歇采樣轉(zhuǎn)發(fā)干擾的功率強度決定。在干擾強度較大的情況下,可以設置較小的權(quán)值以使得聯(lián)合設計的波形和濾波器能夠獲得更優(yōu)的抗干擾性能。而當干擾強度較小時,則設置較大的權(quán)值,使得聯(lián)合設計的波形和濾波器能夠在維持更好的脈沖壓縮性能的同時獲得抗間歇采樣轉(zhuǎn)發(fā)干擾性能。為了平衡設計波形的抗干擾性能和脈沖壓縮性能,本文將帕累托權(quán)值設置為0.5。

    圖2為帕累托值為0.5時本文所提方法設計的波形實部和幅度圖,其滿足恒模約束。圖3為設計的波形及其對應的干擾信號經(jīng)過非匹配濾波輸出后的結(jié)果圖。圖3(a)為波形非匹配濾波脈沖壓縮結(jié)果圖,其歸一化峰值增益損耗為–1.03 dB,積分旁瓣比為–12.46 dB。圖3(b)為干擾信號非匹配濾波輸出結(jié)果圖,其歸一化峰值干擾為–28.43 dB,積分能量比為–14.22 dB。因此,其達到了同時抑制波形脈沖壓縮積分旁瓣和間歇采樣轉(zhuǎn)發(fā)干擾的目的。

    圖2 波形實部幅度圖Fig.2 The amplitude of the real part of waveform

    圖3 波形和干擾信號非匹配濾波輸出(N=512)Fig.3 Mismatch output of the waveform and jamming signal (N=512)

    仿真2:分析峰值增益損耗對波形非匹配濾波輸出積分旁瓣性能和干擾信號積分能量性能的影響。峰值增益損耗[–0.8,–1.0,–1.3,–1.6,–2.0],其余參數(shù)與上述一致。由于波形中間歇采樣部分的信號段被抑制,所以波形的非匹配濾波輸出存在峰值增益損耗。圖4為雷達波形非匹配濾波輸出積分旁瓣比和干擾信號積分能量隨著峰值增益損耗變化曲線,可以發(fā)現(xiàn),隨著峰值增益損耗絕對值的減小,積分旁瓣比和干擾信號積分能量逐漸增加。因此,為了有效抑制間歇采樣轉(zhuǎn)發(fā)干擾,并且提高波形脈沖壓縮積分旁瓣水平,峰值增益損耗不應太小。

    圖4 積分旁瓣比和積分能量比隨峰值增益損耗變化Fig.4 Integrated sidelobe levels and integrated levels versus loss of processing gain

    仿真3:對比分析本文所提算法和傳統(tǒng)方法的運行時間。假設波形和非匹配濾波器長度為128,干擾信號非匹配濾波輸出歸一化峰值為–30 dB。間歇采樣信號占空比0.25,間歇采樣周期重復Ts=Tp/4。算法截止運行閾值ζ=1e?7,帕累托權(quán)值ε=0.5。首先利用文獻[17]所提出的方法設計一組波形和非匹配濾波器,其性能如表2所示。為了公平比較,本文所提算法和文獻[18]方法的波形非匹配濾波輸出歸一化峰值增益損耗為LPG=–1.7 dB,權(quán)值λ1=λ2=256,且其余參數(shù)設置與文獻[17]算法一致。除此之外,所有對比方法的波形和非匹配濾波器初始序列均一致。圖5為兩種方法目標函數(shù)隨運行時間變化曲線,可以發(fā)現(xiàn)本文所提方法速度遠快于本文的對比方法。

    圖5 運行時間對比圖Fig.5 Comparison of running time

    表2 算法性能對比Tab.2 Performance comparison of tested algorithms

    表2為兩種算法設計波形和非匹配濾波器的性能對比,可以發(fā)現(xiàn)在相同峰值增益損耗情況下,本文所提方法的波形非匹配濾波脈沖壓縮積分旁瓣比和干擾信號非匹配濾波輸出積分能量值均低于文獻[17]方法。但是本文設計方法干擾信號非匹配濾波輸出峰值大于文獻[17]方法,這是因為本文設計方法的干擾峰值是通過最小化目標函數(shù)控制的,其值只能無限趨向于零,而難以精確控制為零。并且,可以發(fā)現(xiàn)在脈壓性能和抗干擾性能接近的情況下,本文所提算法的速度顯著提升。

    仿真4:為了評估本文分析提出的理論計算復雜度,開展了一組仿真。假設設計波形長度為[32,64,128,256,512,1024],對每個長度波形開展10組仿真,并記錄每次運行時間。其中算法的截止條件仍為ζ=1e?7。圖6為不同波形長度下的平均運行時間圖。不同信號長度波形設計的算法平均運行時間分別為2.23,10.11,43.57,125.91,506.61和2027.12 s??梢园l(fā)現(xiàn),信號長度每增加一倍,算法運行時間增加3倍。這說明,隨著信號長度增加,算法的運行時間按照平方次數(shù)增加,進而說明本文所提方法計算復雜度為O(N2)。

    圖6 不同波形長度下的運行時間圖Fig.6 Running time versus different code length

    4.2 抗干擾性能評估

    仿真5:評估本文所提方法抗干擾性能。仿真場景參數(shù)如表3所示。以SPL文獻中相關結(jié)果作為參考,設置本文方法相關參數(shù),其中峰值增益損耗為–1.6 dB。圖7為本文設計的波形和干擾信號非匹配濾波輸出結(jié)果。表4為兩種方法設計的仿真結(jié)果性能對比表??梢钥闯?,設計的波形非匹配濾波脈沖壓縮峰值增益損耗一樣,而本文所提方法的積分旁瓣比和積分能量比性能均優(yōu)于文獻[17]方法。

    表3 仿真參數(shù)Tab.3 Simulation parameters

    表4 兩組波形和非匹配濾波器性能對比Tab.4 The mismatch output performance of designed waveform via tested algorithms

    圖7 波形和干擾信號非匹配濾波輸出(N=800)Fig.7 Mismatch output of waveform and jamming signal (N=800)

    仿真6:為評估本文的抗干擾性能,將本文方法和文獻[10]、文獻[13]、文獻[17]得到的抗干擾結(jié)果進行了對比分析。仿真參數(shù)如表3所示,干擾機的轉(zhuǎn)發(fā)模式為重復轉(zhuǎn)發(fā)干擾,轉(zhuǎn)發(fā)次數(shù)為2次。信噪比為0 dB,信干比為–20 dB。圖8為抗干擾得到的仿真結(jié)果,其中所有結(jié)果均通過發(fā)射信號能量進行歸一化。圖8(a)—圖8(d)分別為文獻[13]、文獻[10]、文獻[17]和所提方法得到的仿真結(jié)果。表5為抗干擾仿真結(jié)果性能對比??梢园l(fā)現(xiàn),本文和文獻[17]提出的方法信號處理增益存在損耗,這是由非匹配濾波器引起的。相較于文獻[13]和文獻[17]提出的方法,本文所提出的方法的干擾峰值是最低的,并且干擾導致的假目標的數(shù)量也是最少的。但是,相較于文獻[10]提出的方法,本文的干擾峰值顯著偏高,并且其旁瓣的值也顯著低于其他3類對比方法。這是因為在干擾和目標距離可分辨情況下,通過干擾識別出目標所在距離單元后,文獻[10]提出的方法可在時頻域?qū)崿F(xiàn)對干擾能量的完全抑制。

    表5 抗干擾性能對比Tab.5 Performance comparison of jamming suppression

    圖8 抗間歇采樣轉(zhuǎn)發(fā)干擾評估Fig.8 Evaluation of interrupted sampling repeater jamming suppression

    其次,進一步分析干擾機位置變化對本文方法和文獻[10]抗干擾性能,以驗證本文所提方法的優(yōu)越性。假設干擾機的位置從2000 m到3000 m逐漸增加,增加的間距為20 m。其余的仿真參數(shù)與表2一致。同時每個位置開展100次蒙特卡羅仿真,并記錄其干擾峰值。圖9為本文方法和文獻[10]抗干擾性能隨著干擾機位置變化的結(jié)果圖。可以發(fā)現(xiàn),當干擾機和目標位置不可分離時,文獻[10]所提方法的抗干擾性能將顯著下降。這是因為,當干擾目標和實際目標處于臨近距離單元時,時頻域濾波將難以有效實現(xiàn)對干擾信號的抑制。這說明,當干擾目標和實際目標可分離時,文獻[10]提出的方法抗干擾性能優(yōu)于本文所提方法。但是,本文所提方法僅需要估計干擾機的占空比和周期參數(shù),設計發(fā)射波形和非匹配濾波器后即可實現(xiàn)干擾信號非匹配輸出峰值和積分能量抑制。相較于文獻[10]所提方法,本文所提方法在信號處理過程中僅需要利用非匹濾波進行脈壓即可得到抗干擾一維距離像,而無需根據(jù)脈壓后的時頻特性差異對干擾和目標回波進行辨識與抑制,故本文所提方法抗干擾性能不受干擾目標和真實目標之間的相對距離影響。當目標和干擾目標不可分離時,本文提出方法的性能顯著優(yōu)于文獻[10]提出的方法。因此,相較于文獻[10]提出的方法,本文所提出的方法抗干擾性能不受干擾目標和真實目標相對距離影響,普適性更強。

    圖9 干擾機位置對抗干擾性能的影響Fig.9 Jamming suppression performance versus location of jammer

    4.3 干擾參數(shù)敏感性分析

    仿真7:為了評估干擾采樣占空比和重復周期兩個關鍵參數(shù)對抗干擾性能的影響,開展了一組參數(shù)敏感性分析實驗。間歇采樣重復周期和占空比的估計誤差變化范圍為±10%,誤差間隔1%取一個周期和占空比數(shù)值。在每個誤差數(shù)值處分別開展100次蒙特卡羅仿真實驗,并記錄脈沖壓縮后的歸一化干擾峰值。本仿真中,設置信干比為–15 dB,信噪比為0 dB,干擾機位置為(0,2200,0) m,其余仿真參數(shù)如表2所示。由圖9可知,由于干擾機的位置為導致干擾目標和真實目標之間不可分離,因此文獻[10]所提方法此時無法有效抑制干擾。此外,我們與文獻[13]所示方法進行了進一步對比,且忽略其干擾信號重構(gòu)過程中的干擾信號幅度誤差。需要說明的是,由于假設其干擾重構(gòu)幅度是無誤差的,因此在間歇采樣重復周期和占空比精確也無誤差時,其抗干擾性能會比圖8(a)對應的結(jié)果更好。

    圖10為干擾峰值隨著重復周期和占空比的估計誤差變化圖。由圖10(a)可知,文獻[13]所提的方法對間歇采樣轉(zhuǎn)發(fā)重復周期十分敏感。這是因為,間歇采樣轉(zhuǎn)發(fā)重復周期的誤差均會導致重構(gòu)干擾信號與實際干擾信號之間存在顯著差異,進而導致抗干擾性能的下降。本文所提的抗干擾性能受間歇采樣轉(zhuǎn)發(fā)重復周期的影響較大,當相對誤差變化大于1.5%時,干擾輸出峰值大于–5 dB。由圖10(b)可知,當相對誤差為負數(shù),即實際干擾占空比小于聯(lián)合設計過程中設置的占空比,本文所提方法抗干擾性能幾乎不受影響,而文獻[13]所提方法的抗干擾性能則會逐漸降低。當實際干擾占空比大于參數(shù)估計得到占空比時,本文所提方法抗干擾性能惡化速度快于文獻[13]所提的方法。當相對誤差大于5%時,本文所提方法的干擾峰值大于–5 dB。因此,相較于文獻[13]所提的方法,本文所提方法與其間歇采樣周期敏感性基本一致。而兩種方法對占空比的敏感性則視估計誤差而定。如果估計的占空比值大于實際干擾機占空比,即估計誤差為負數(shù)時,本文所提方法對占空比的敏感性低于文獻[13]所提的方法;反之,則本文所提方法對占空比參數(shù)更敏感。

    圖10 干擾峰值隨著重復周期和占空比的估計誤差變化圖Fig.10 Jamming peak versus the relative error of pulse repetition period and duty ratio

    5 結(jié)論

    本文研究了抗間歇采樣轉(zhuǎn)發(fā)干擾的快速波形和濾波器聯(lián)合設計算法。相較于傳統(tǒng)方法,本文設計的波形和濾波器可以實現(xiàn)抑制間歇采樣轉(zhuǎn)發(fā)干擾。本文所提出的設計算法計算復雜度為O(N2),顯著優(yōu)于傳統(tǒng)方法。仿真結(jié)果驗證了本文所提加速算法的有效性,其算法運行速度顯著優(yōu)于傳統(tǒng)方法。此外,相較于傳統(tǒng)的時頻信號處理方法,本文所提方法能夠在干擾目標和真實目標距離不可分離時仍達到抑制干擾的目的。本文所提算法是在間歇采樣信號周期和占空比先驗知識抑制情況下開展的波形和濾波器聯(lián)合設計,并且本文所提方法對干擾采樣周期十分敏感。因此,干擾信號周期和占空比參數(shù)的精確估計是下一步研究的重點。

    猜你喜歡
    干擾信號間歇濾波器
    基于無擾濾波器和AED-ADT的無擾切換控制
    間歇供暖在散熱器供暖房間的應用
    煤氣與熱力(2022年4期)2022-05-23 12:44:46
    正弦采樣信號中單一脈沖干擾信號的快速剔除實踐方法
    基于粒子群算法的光纖通信干擾信號定位方法
    從濾波器理解卷積
    電子制作(2019年11期)2019-07-04 00:34:38
    開關電源EMI濾波器的應用方法探討
    電子制作(2018年16期)2018-09-26 03:26:50
    管群間歇散熱的土壤溫度響應與恢復特性
    淺析監(jiān)控干擾信號的優(yōu)化處置措施
    基于TMS320C6678的SAR方位向預濾波器的并行實現(xiàn)
    相參雷達典型干擾信號產(chǎn)生及關鍵技術(shù)
    久热这里只有精品99| 久久人人爽人人爽人人片va| 国产精品国产三级专区第一集| 肉色欧美久久久久久久蜜桃| 成人亚洲欧美一区二区av| 黑丝袜美女国产一区| 成年人免费黄色播放视频| 色5月婷婷丁香| 精品国产一区二区三区四区第35| 国产乱来视频区| 一二三四在线观看免费中文在 | 一级毛片我不卡| 在线观看免费高清a一片| 久久99热这里只频精品6学生| 亚洲情色 制服丝袜| 搡老乐熟女国产| 大片电影免费在线观看免费| av又黄又爽大尺度在线免费看| 桃花免费在线播放| 亚洲综合色惰| 制服诱惑二区| 欧美精品一区二区免费开放| 亚洲精品乱码久久久久久按摩| 一级片免费观看大全| 大陆偷拍与自拍| 国产av码专区亚洲av| 国产av精品麻豆| 宅男免费午夜| 欧美97在线视频| 国产免费又黄又爽又色| 亚洲图色成人| 伊人久久国产一区二区| av在线老鸭窝| 免费黄频网站在线观看国产| 高清在线视频一区二区三区| 亚洲精品一区蜜桃| 老女人水多毛片| 97在线人人人人妻| 狠狠精品人妻久久久久久综合| 亚洲国产成人一精品久久久| 韩国高清视频一区二区三区| 国产无遮挡羞羞视频在线观看| 亚洲成人一二三区av| 国产爽快片一区二区三区| 亚洲图色成人| 黄色 视频免费看| 久久99热6这里只有精品| 午夜激情av网站| 欧美成人午夜精品| 国产乱来视频区| 欧美3d第一页| 久久99热6这里只有精品| 国产精品国产三级国产专区5o| 久久精品久久久久久久性| 啦啦啦啦在线视频资源| 999精品在线视频| 成人二区视频| 国产爽快片一区二区三区| 国产精品久久久久久久久免| 婷婷成人精品国产| 中文字幕免费在线视频6| 老熟女久久久| 国产在线一区二区三区精| 亚洲国产毛片av蜜桃av| 日韩一区二区视频免费看| 赤兔流量卡办理| 久久久欧美国产精品| 欧美变态另类bdsm刘玥| 久久精品熟女亚洲av麻豆精品| 成人午夜精彩视频在线观看| 亚洲精华国产精华液的使用体验| 久久狼人影院| 免费少妇av软件| 久久精品熟女亚洲av麻豆精品| 国产精品秋霞免费鲁丝片| 亚洲美女黄色视频免费看| 欧美亚洲 丝袜 人妻 在线| 欧美亚洲 丝袜 人妻 在线| 超色免费av| 国产 一区精品| 国产极品粉嫩免费观看在线| 少妇人妻久久综合中文| 国产一区二区在线观看av| 搡女人真爽免费视频火全软件| 男人操女人黄网站| 午夜av观看不卡| 亚洲精品日韩在线中文字幕| 国产黄色免费在线视频| 亚洲精品乱码久久久久久按摩| 精品人妻偷拍中文字幕| 亚洲精品自拍成人| 亚洲精品日韩在线中文字幕| 欧美日韩成人在线一区二区| 午夜日本视频在线| 日产精品乱码卡一卡2卡三| 丝袜人妻中文字幕| 91国产中文字幕| 汤姆久久久久久久影院中文字幕| 天天操日日干夜夜撸| 亚洲欧美一区二区三区黑人 | 成人综合一区亚洲| 久久精品aⅴ一区二区三区四区 | 夜夜骑夜夜射夜夜干| 亚洲高清免费不卡视频| 国产免费视频播放在线视频| 王馨瑶露胸无遮挡在线观看| 久久97久久精品| 亚洲国产精品一区三区| 视频中文字幕在线观看| 亚洲 欧美一区二区三区| 日韩成人av中文字幕在线观看| 99久久综合免费| 热99国产精品久久久久久7| 美国免费a级毛片| av女优亚洲男人天堂| 午夜激情av网站| 侵犯人妻中文字幕一二三四区| 丰满饥渴人妻一区二区三| 久久久久久久国产电影| 精品99又大又爽又粗少妇毛片| 一个人免费看片子| 国产精品国产av在线观看| 亚洲,欧美精品.| 免费在线观看完整版高清| 亚洲性久久影院| 女性生殖器流出的白浆| 亚洲人成77777在线视频| 边亲边吃奶的免费视频| 天天躁夜夜躁狠狠久久av| 哪个播放器可以免费观看大片| 桃花免费在线播放| 国产精品一区二区在线不卡| av卡一久久| a级毛片黄视频| 国产亚洲一区二区精品| 纯流量卡能插随身wifi吗| av在线老鸭窝| 国产成人一区二区在线| 免费大片黄手机在线观看| 亚洲精品美女久久av网站| 精品少妇黑人巨大在线播放| 国产高清不卡午夜福利| 日本黄大片高清| 极品人妻少妇av视频| 国产激情久久老熟女| 黑人猛操日本美女一级片| 亚洲av电影在线观看一区二区三区| 高清不卡的av网站| 99九九在线精品视频| 高清欧美精品videossex| 在现免费观看毛片| 国产亚洲一区二区精品| 国产高清三级在线| 亚洲性久久影院| 插逼视频在线观看| 日韩av在线免费看完整版不卡| 国产精品成人在线| 国产精品国产三级国产av玫瑰| 狂野欧美激情性xxxx在线观看| 两个人看的免费小视频| 国产成人免费无遮挡视频| 亚洲精品日本国产第一区| 婷婷色av中文字幕| 国产色爽女视频免费观看| 少妇的逼水好多| 岛国毛片在线播放| 制服丝袜香蕉在线| 亚洲色图 男人天堂 中文字幕 | 久久久久久久精品精品| 欧美日韩国产mv在线观看视频| 国产精品99久久99久久久不卡 | 欧美日韩一区二区视频在线观看视频在线| 免费女性裸体啪啪无遮挡网站| 一本色道久久久久久精品综合| av片东京热男人的天堂| 国产精品国产三级国产av玫瑰| 免费不卡的大黄色大毛片视频在线观看| 免费av中文字幕在线| 久久人人爽av亚洲精品天堂| 天堂8中文在线网| videossex国产| 亚洲欧美中文字幕日韩二区| 丰满饥渴人妻一区二区三| 亚洲美女搞黄在线观看| 免费看不卡的av| 欧美激情极品国产一区二区三区 | 人人妻人人添人人爽欧美一区卜| 午夜91福利影院| 51国产日韩欧美| 日韩制服丝袜自拍偷拍| 亚洲欧美成人精品一区二区| 色网站视频免费| 国产精品国产三级专区第一集| 亚洲欧洲国产日韩| 午夜激情久久久久久久| xxxhd国产人妻xxx| 精品国产国语对白av| 国产无遮挡羞羞视频在线观看| 精品亚洲成国产av| 91精品三级在线观看| 亚洲国产精品一区二区三区在线| 国产熟女午夜一区二区三区| 男女边吃奶边做爰视频| 丰满少妇做爰视频| 水蜜桃什么品种好| 18禁观看日本| 亚洲丝袜综合中文字幕| 性高湖久久久久久久久免费观看| 男女边摸边吃奶| 七月丁香在线播放| 啦啦啦中文免费视频观看日本| 午夜老司机福利剧场| 秋霞在线观看毛片| 精品亚洲乱码少妇综合久久| 深夜精品福利| 人成视频在线观看免费观看| 日本色播在线视频| 免费久久久久久久精品成人欧美视频 | 国产xxxxx性猛交| 欧美+日韩+精品| 永久网站在线| 制服诱惑二区| 欧美精品av麻豆av| 一本大道久久a久久精品| 丝袜在线中文字幕| 成年女人在线观看亚洲视频| 国产一区二区激情短视频 | 国产成人a∨麻豆精品| 在现免费观看毛片| 久久久国产欧美日韩av| 深夜精品福利| 美女大奶头黄色视频| 国产男人的电影天堂91| 午夜91福利影院| 欧美日韩亚洲高清精品| 国产又色又爽无遮挡免| 亚洲欧美一区二区三区黑人 | 久久99蜜桃精品久久| 国产成人免费观看mmmm| 国产日韩一区二区三区精品不卡| 国产精品久久久久成人av| 亚洲美女搞黄在线观看| 精品人妻熟女毛片av久久网站| 蜜臀久久99精品久久宅男| 中文精品一卡2卡3卡4更新| 久久精品国产综合久久久 | 你懂的网址亚洲精品在线观看| 久久久久久伊人网av| 欧美人与性动交α欧美软件 | 在线看a的网站| 久久精品久久久久久噜噜老黄| 青春草亚洲视频在线观看| 午夜激情av网站| 97在线人人人人妻| 国产在线视频一区二区| 91午夜精品亚洲一区二区三区| 色吧在线观看| 9191精品国产免费久久| 国产69精品久久久久777片| 日韩视频在线欧美| 熟女电影av网| 99香蕉大伊视频| 国产成人av激情在线播放| 欧美成人午夜精品| 一个人免费看片子| 色94色欧美一区二区| 欧美日韩精品成人综合77777| 欧美国产精品一级二级三级| 又黄又粗又硬又大视频| 国产精品久久久久久精品古装| 欧美3d第一页| 久久婷婷青草| 国产在线一区二区三区精| 十八禁网站网址无遮挡| 黄片无遮挡物在线观看| 国产高清国产精品国产三级| www日本在线高清视频| 国产精品蜜桃在线观看| 亚洲,欧美精品.| 日韩av在线免费看完整版不卡| 少妇的丰满在线观看| 一级毛片我不卡| 亚洲图色成人| 黑人欧美特级aaaaaa片| 最近中文字幕高清免费大全6| 一边摸一边做爽爽视频免费| 亚洲第一av免费看| 天堂中文最新版在线下载| 亚洲成国产人片在线观看| 亚洲欧美日韩卡通动漫| 少妇 在线观看| 日本欧美视频一区| 精品熟女少妇av免费看| 国产成人91sexporn| 9191精品国产免费久久| 欧美 日韩 精品 国产| 人人妻人人爽人人添夜夜欢视频| 免费av不卡在线播放| 精品国产乱码久久久久久小说| 精品少妇久久久久久888优播| a级毛色黄片| 免费看不卡的av| 成人手机av| 成人漫画全彩无遮挡| 毛片一级片免费看久久久久| 日日摸夜夜添夜夜爱| 九九爱精品视频在线观看| 中文字幕人妻丝袜制服| 大香蕉久久成人网| 中文天堂在线官网| 桃花免费在线播放| 咕卡用的链子| 国产高清三级在线| 纵有疾风起免费观看全集完整版| 亚洲国产欧美在线一区| 在现免费观看毛片| 考比视频在线观看| 香蕉国产在线看| 综合色丁香网| 欧美成人精品欧美一级黄| 伦理电影大哥的女人| 亚洲国产av新网站| 丝袜脚勾引网站| 日韩精品免费视频一区二区三区 | 丝袜人妻中文字幕| 青春草亚洲视频在线观看| 欧美3d第一页| 五月开心婷婷网| 国产麻豆69| 性色avwww在线观看| 成人综合一区亚洲| 久久久亚洲精品成人影院| 国产精品欧美亚洲77777| 久久99热6这里只有精品| 国产av码专区亚洲av| xxx大片免费视频| 国产日韩欧美视频二区| 宅男免费午夜| 亚洲欧洲国产日韩| 久久国内精品自在自线图片| 久久婷婷青草| 亚洲精品自拍成人| 亚洲国产精品一区三区| 黄色 视频免费看| 国产成人精品在线电影| av又黄又爽大尺度在线免费看| 午夜激情av网站| 亚洲久久久国产精品| 婷婷色综合大香蕉| 一区二区三区乱码不卡18| 国产xxxxx性猛交| 日韩av在线免费看完整版不卡| 欧美少妇被猛烈插入视频| 热re99久久精品国产66热6| 久久女婷五月综合色啪小说| 中文精品一卡2卡3卡4更新| 久久国内精品自在自线图片| 精品国产一区二区三区四区第35| 国产免费一级a男人的天堂| 天堂8中文在线网| 亚洲精品一区蜜桃| 伊人亚洲综合成人网| 久久综合国产亚洲精品| 国产爽快片一区二区三区| 99视频精品全部免费 在线| 亚洲av福利一区| 国产精品久久久久久精品电影小说| 日韩在线高清观看一区二区三区| 久久精品国产自在天天线| 大片免费播放器 马上看| 成年女人在线观看亚洲视频| 免费高清在线观看视频在线观看| 黄色怎么调成土黄色| 91成人精品电影| 丝袜人妻中文字幕| 黑人欧美特级aaaaaa片| 韩国精品一区二区三区 | 午夜福利影视在线免费观看| 免费人成在线观看视频色| 亚洲精品,欧美精品| 美女内射精品一级片tv| 欧美变态另类bdsm刘玥| 五月伊人婷婷丁香| 国语对白做爰xxxⅹ性视频网站| 热99久久久久精品小说推荐| 成人黄色视频免费在线看| 久久国内精品自在自线图片| 18禁国产床啪视频网站| 大香蕉久久成人网| 91精品伊人久久大香线蕉| 美女中出高潮动态图| 街头女战士在线观看网站| 欧美日韩国产mv在线观看视频| 午夜视频国产福利| 国产一区二区三区av在线| 久久国内精品自在自线图片| 国产毛片在线视频| 日本wwww免费看| 日韩 亚洲 欧美在线| 亚洲性久久影院| 美女国产视频在线观看| 99热6这里只有精品| 成人免费观看视频高清| 五月玫瑰六月丁香| 欧美丝袜亚洲另类| 免费观看在线日韩| 青春草亚洲视频在线观看| 精品亚洲成国产av| 91国产中文字幕| 我要看黄色一级片免费的| 亚洲性久久影院| 成人国产av品久久久| 亚洲成人手机| 亚洲高清免费不卡视频| 美女中出高潮动态图| 亚洲国产av新网站| 久久精品久久久久久久性| 婷婷成人精品国产| 久久久久久久久久久免费av| 亚洲 欧美一区二区三区| 人人妻人人爽人人添夜夜欢视频| 午夜久久久在线观看| 欧美97在线视频| 极品人妻少妇av视频| 少妇熟女欧美另类| 亚洲av欧美aⅴ国产| 丝袜美足系列| 亚洲精华国产精华液的使用体验| 国产男人的电影天堂91| 熟妇人妻不卡中文字幕| 母亲3免费完整高清在线观看 | 国产免费又黄又爽又色| 欧美bdsm另类| 欧美丝袜亚洲另类| 亚洲精品aⅴ在线观看| 亚洲成人手机| 国产在线免费精品| 精品少妇久久久久久888优播| av视频免费观看在线观看| 国产日韩欧美在线精品| 亚洲av日韩在线播放| 人人妻人人添人人爽欧美一区卜| 日韩一区二区视频免费看| 人人妻人人澡人人看| 久久精品夜色国产| 欧美日韩一区二区视频在线观看视频在线| 综合色丁香网| 日韩精品免费视频一区二区三区 | 又大又黄又爽视频免费| 欧美日韩av久久| 99精国产麻豆久久婷婷| 久久99热这里只频精品6学生| 日韩一区二区视频免费看| 日本wwww免费看| 青春草国产在线视频| 99热这里只有是精品在线观看| 国产亚洲av片在线观看秒播厂| 中文欧美无线码| 亚洲国产精品成人久久小说| 免费高清在线观看日韩| 国国产精品蜜臀av免费| 午夜精品国产一区二区电影| 一级片免费观看大全| 99久久综合免费| 国产成人一区二区在线| 黑人高潮一二区| 18禁观看日本| 在线观看免费日韩欧美大片| 色婷婷久久久亚洲欧美| 全区人妻精品视频| 各种免费的搞黄视频| 人妻系列 视频| 日本-黄色视频高清免费观看| 免费观看性生交大片5| 日韩精品免费视频一区二区三区 | 亚洲欧美成人综合另类久久久| 麻豆精品久久久久久蜜桃| 国产精品久久久久久久电影| 最近中文字幕高清免费大全6| 久久久欧美国产精品| av又黄又爽大尺度在线免费看| 久热久热在线精品观看| 一级,二级,三级黄色视频| 毛片一级片免费看久久久久| 亚洲,欧美,日韩| 免费久久久久久久精品成人欧美视频 | 最近手机中文字幕大全| 色5月婷婷丁香| 一边亲一边摸免费视频| 国产色爽女视频免费观看| 久久久久精品久久久久真实原创| 亚洲精品第二区| 满18在线观看网站| 美女内射精品一级片tv| 亚洲在久久综合| 精品一区在线观看国产| 亚洲欧美色中文字幕在线| 国产国语露脸激情在线看| 国产亚洲一区二区精品| 亚洲激情五月婷婷啪啪| 精品一区二区三区视频在线| 久久国产精品大桥未久av| 天美传媒精品一区二区| 人妻 亚洲 视频| av又黄又爽大尺度在线免费看| 精品99又大又爽又粗少妇毛片| 欧美日韩视频精品一区| 99热6这里只有精品| 日韩一本色道免费dvd| 亚洲四区av| 99国产综合亚洲精品| 一二三四中文在线观看免费高清| 欧美激情 高清一区二区三区| 国产成人精品婷婷| 涩涩av久久男人的天堂| 国产麻豆69| h视频一区二区三区| 波多野结衣一区麻豆| 一边摸一边做爽爽视频免费| 午夜福利在线观看免费完整高清在| 18禁在线无遮挡免费观看视频| 日韩熟女老妇一区二区性免费视频| 午夜免费鲁丝| 人人妻人人爽人人添夜夜欢视频| 97在线视频观看| 午夜免费观看性视频| 中国三级夫妇交换| videos熟女内射| 久久久久精品久久久久真实原创| 国产亚洲一区二区精品| 亚洲激情五月婷婷啪啪| 亚洲国产精品成人久久小说| 亚洲国产精品一区二区三区在线| 丰满少妇做爰视频| 国产白丝娇喘喷水9色精品| 久久久久久久精品精品| 尾随美女入室| 久久ye,这里只有精品| 久久这里有精品视频免费| 亚洲精品色激情综合| 国产高清三级在线| 成人国语在线视频| 亚洲少妇的诱惑av| 搡女人真爽免费视频火全软件| 欧美日韩亚洲高清精品| 成人黄色视频免费在线看| 性色av一级| 如日韩欧美国产精品一区二区三区| 久久久精品区二区三区| 18在线观看网站| 少妇熟女欧美另类| 男女啪啪激烈高潮av片| 欧美日韩av久久| 国产精品 国内视频| 97人妻天天添夜夜摸| 交换朋友夫妻互换小说| 9热在线视频观看99| a级毛片在线看网站| 韩国av在线不卡| 国产黄色免费在线视频| 韩国av在线不卡| 国产精品国产三级专区第一集| 亚洲人与动物交配视频| 久久免费观看电影| 欧美日韩av久久| 久久久亚洲精品成人影院| 狂野欧美激情性bbbbbb| 天天躁夜夜躁狠狠久久av| 狂野欧美激情性bbbbbb| 亚洲一区二区三区欧美精品| a 毛片基地| 视频在线观看一区二区三区| 亚洲欧美一区二区三区国产| 亚洲av国产av综合av卡| 国产乱来视频区| 91国产中文字幕| 亚洲精品中文字幕在线视频| 国产欧美日韩综合在线一区二区| 久久久久久人妻| 亚洲伊人久久精品综合| 国产乱人偷精品视频| 最近最新中文字幕免费大全7| 老女人水多毛片| 国产激情久久老熟女| 少妇被粗大的猛进出69影院 | a级毛片黄视频| 久久久亚洲精品成人影院| 国产乱来视频区| 免费人妻精品一区二区三区视频| 婷婷色av中文字幕| 久久精品国产自在天天线| 九九爱精品视频在线观看| 99久久中文字幕三级久久日本| 国产日韩欧美在线精品| 国产精品免费大片| 国产精品国产av在线观看| 9色porny在线观看| 久久鲁丝午夜福利片| 亚洲婷婷狠狠爱综合网| 免费观看av网站的网址| 人妻少妇偷人精品九色| 免费高清在线观看视频在线观看| 亚洲精品中文字幕在线视频| 亚洲美女视频黄频| 国产亚洲av片在线观看秒播厂| 大码成人一级视频| 国产又色又爽无遮挡免| 日韩成人av中文字幕在线观看| 久久精品人人爽人人爽视色| 婷婷色综合www| 美女大奶头黄色视频| 建设人人有责人人尽责人人享有的| 婷婷色综合大香蕉| 亚洲国产精品国产精品| 亚洲精品乱久久久久久|