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

    抗間歇采樣轉(zhuǎn)發(fā)干擾的發(fā)射波形與失配濾波器聯(lián)合優(yōu)化算法

    2024-11-23 00:00:00劉茜戴奉周

    摘要: 間歇采樣轉(zhuǎn)發(fā)干擾(interrupted sampling repeater jamming, ISRJ)與發(fā)射信號之間具有強(qiáng)相關(guān)性,若在接收端采用匹配濾波處理則會產(chǎn)生多個數(shù)量可控的高逼真假目標(biāo),對雷達(dá)的檢測性能具有極大的干擾。針對上述問題,提出了一種抗ISRJ的發(fā)射波形和接收濾波器聯(lián)合優(yōu)化算法,將失配濾波體制下的信噪比(signal to noise ratio, SNR)損失、發(fā)射信號恒模約束以及濾波器能量約束考慮在內(nèi),以最小化干擾信號歸一化脈壓后峰值、失配濾波輸出信號積分副主比以及輸出干信比,并采用主分量最小化(majorization-minimization, MM)方法與平方迭代加速算法提高算法運(yùn)行速度。仿真結(jié)果表明,與其他同類方法相比,所提算法能在保證可靠抗干擾性能的同時,極大縮短運(yùn)行時間,具有較好的實(shí)時性。

    關(guān)鍵詞: 間歇采樣轉(zhuǎn)發(fā)干擾; 失配濾波; 主分量最小化; 平方迭代方法; 恒模約束

    中圖分類號: TN 974

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

    DOI:10.12305/j.issn.1001-506X.2024.06.04

    A joint optimization algorithm for transmitting waveform and mismatched filter against interrupted sampling repeater jamming

    LIU Qian, DAI Fengzhou*

    (School of Electronic Engineering, Xidian University, Xi’an 710071, China)

    Abstract: Interrupted sampling repeater jamming (ISRJ) has a strong correlation with the transmitting signal, if the matched filtering processing is still adopted at the receiving end, multiple high-fidelity 1 targets with controllable quantity will be generated, which seriously affects the detection performance of radar. To solve the above problem, a joint optimization algorithm for transmitting waveform and receiving filter against the ISRJ is proposed, which takes into account the signal to noise ratio (SNR) loss under the mismatched filtering system, the constant modulus constraint of transmitting signal and the energy constraint of filter, to minimize the normalized peak of the interference signal after pulse compression, the ratio of the energy of the integral side lobe to main lobe and the interference to signal ratio of the output of the mismatched filter, the majorization-minimization (MM) method and the square iterative acceleration algorithm are adopted to accelerate the running speed of the proposed algorithm. The simulation results show that compared with other similar methods, the proposed algorithm can greatly shorten the running time and has better real-time performance on the basis of ensuring reliable anti-interference performance.

    Keywords: interrupted sampling repeater jamming (ISRJ); mismatched filter; majorization-minimization (MM); squared iterative method; constant modulus constraint

    0 引 言

    間歇采樣轉(zhuǎn)發(fā)干擾(interrupted sampling repeater jamming, ISRJ)是一種基于數(shù)字射頻存儲技術(shù)的相干干擾[1-5,通過對雷達(dá)發(fā)射信號進(jìn)行間斷性采樣并轉(zhuǎn)發(fā),可以在接收端產(chǎn)生多個假目標(biāo)[6-8,對雷達(dá)檢測、跟蹤以及目標(biāo)識別具有極大的欺騙性。因此,抗ISRJ技術(shù)已成為當(dāng)前的研究熱點(diǎn)。

    當(dāng)前抗ISRJ的方法主要可分為兩大類,分別為在接收端通過信號處理濾除干擾的方法以及通過波形設(shè)計(jì)的抗干擾方法。對于第一類方法,文獻(xiàn)[1]指出ISRJ信號的時頻特性在脈沖持續(xù)時間內(nèi)是不連續(xù)的,并利用了ISRJ信號與目標(biāo)回波的這一時頻差異性,在接收端針對ISRJ信號構(gòu)造了多個帶通濾波器,以實(shí)現(xiàn)對ISRJ信號的抑制。文獻(xiàn)[9]基于雙向堆疊門控遞歸單元(stacked bidirectional gated recurrent unit, SBiGRU)網(wǎng)絡(luò)提出了抗ISRJ的方法:將信號提取轉(zhuǎn)化為時間分類問題,生成帶通濾波器,在抑制ISRJ信號的同時可保留真實(shí)目標(biāo)信號分量,從而保證可靠的抗干擾性能。然而,以上方案均為被動抗干擾策略,未考慮發(fā)射波形的優(yōu)化問題,使得算法可設(shè)計(jì)的自由度較低,且上述算法僅適用于信干比高的場景,而在信干比較低,甚至是信號被干擾信號淹沒的工作場景下則無法發(fā)揮理想的抗干擾效果。

    對于第二類抗ISRJ的方法,不同的發(fā)射信號波形無疑具有不同的抗干擾能力[10-18。當(dāng)前,為保障可靠的探測及抗干擾性能,在工程應(yīng)用中雷達(dá)廣泛采用脈內(nèi)調(diào)制信號,其主要分為脈內(nèi)調(diào)頻與脈內(nèi)調(diào)相兩種形式。較為典型的抗干擾信號波形即為宙斯盾系統(tǒng)[10中雷達(dá)發(fā)射的脈內(nèi)多載頻信號,其脈沖內(nèi)包含多個子脈沖,各子脈沖內(nèi)部進(jìn)行了相位編碼處理,且子脈沖之間伴有隨機(jī)的頻率編碼。研究表明,脈內(nèi)多載頻相位編碼信號具有比傳統(tǒng)的單載頻相位編碼信號更為可靠的抗干擾性能[10。張彥等12則對雙曲調(diào)頻(hyperbolic frequency modulated, HFM)波形的抗ISRJ性能進(jìn)行了分析,并指出了HFM波形相較于線性調(diào)頻信號具有更好的抗干擾性能。文獻(xiàn)[14]提出了一種用于抗ISRJ的自適應(yīng)相位編碼信號設(shè)計(jì)方案,利用遺傳算法優(yōu)化得出與干擾信號正交的發(fā)射信號,進(jìn)而可通過脈沖壓縮對ISRJ信號進(jìn)行有效抑制。文獻(xiàn)[15]同樣以設(shè)計(jì)相正交的工作波形與保護(hù)波形的思路來抑制ISRJ。但上述方法都僅僅對發(fā)射波形進(jìn)行優(yōu)化,而未考慮接收濾波器,算法可設(shè)計(jì)的自由度較低,為獲得更為可靠的抗干擾性能可對發(fā)射波形與接收濾波器進(jìn)行聯(lián)合優(yōu)化處理。因此,文獻(xiàn)[16]基于平方迭代加速方法,提出了一種發(fā)射波形和接收濾波器聯(lián)合設(shè)計(jì)的快速算法,在目標(biāo)與干擾機(jī)距離不可分時仍能保證可靠的抗干擾性能,但該算法對于間歇采樣重復(fù)周期的敏感度較高,因此對干擾信號采樣重復(fù)周期的估計(jì)精度提出了較高的要求。文獻(xiàn)[17]則以最小化失配濾波器輸出的旁瓣積分能量以及干擾積分能量為目標(biāo),提出了一種迭代算法,并分別使用了拉格朗日乘子法以及交替方向乘子法對發(fā)射波形和接收濾波器進(jìn)行優(yōu)化,可在抑制ISRJ信號的同時保持良好的脈沖壓縮性能,但該算法計(jì)算復(fù)雜度較高,算法收斂速度慢。

    針對上述問題,本文提出了一種抗ISRJ的發(fā)射波形和接收濾波器聯(lián)合設(shè)計(jì)算法。其基本思路是,以干擾信號歸一化脈壓峰值、失配濾波器輸出的積分副主比以及干信比最小為目標(biāo),同時考慮信號的恒模約束以及接收濾波器的能量約束等因素,并將失配濾波體制下的信噪比損失考慮在內(nèi),構(gòu)造雷達(dá)發(fā)射波形與接收濾波器聯(lián)合設(shè)計(jì)的數(shù)學(xué)模型。并在此基礎(chǔ)上,利用主分量最小化(majorization-minimization, MM)原理分別推導(dǎo)出了發(fā)射波形和失配濾波器在各次交替優(yōu)化迭代過程中的解析表達(dá)式,加快了算法運(yùn)行速度。同時,本文使用平方迭代加速方法,進(jìn)一步加快了算法運(yùn)行速度。最后,仿真結(jié)果表明,本文算法能夠以較小的信噪比損失(signal to noise ratio loss, SNRL)為代價換取可靠的抗干擾性能,且能夠有效地縮短算法運(yùn)行時間。

    1 發(fā)射波形與失配濾波器聯(lián)合優(yōu)化問題建模

    假設(shè)雷達(dá)發(fā)射信號x是恒模相位編碼信號,其長度N可表示為

    x=[x1,x2,…,xNT(1)

    式中:xn=ejφn,n=1,2,…,N,且其相位滿足φn∈[0,2π];(·)T為取轉(zhuǎn)置運(yùn)算。與之相應(yīng)的接收濾波器長度也為N,可表示為h=[h1,h2,…,hNT。進(jìn)一步可定義發(fā)射序列矩陣X與接收濾波器矩陣[16H:

    式中:(·)H表示取共軛轉(zhuǎn)置操作。信號經(jīng)失配濾波脈壓后的積分旁瓣能量可表示為hHXHWXh或xHHHWHx,主瓣能量則為hHXHXh-h(huán)HXHWXh=|hHx|2=hHxxHh,其中W為(2N-1)×(2N-1)的方陣,該矩陣中除第N個對角元素與所有非對角元素為0外,其余元素的為1。定義向量p=[p1,p2,…,pNT為用于間歇采樣的離散序列,則ISRJ可表示為xjam=x⊙p,其中⊙表示哈達(dá)瑪積。那么可進(jìn)一步定義ISRJ信號矩陣Xjam與相應(yīng)的失配濾波器矩陣16Hjam

    失配濾波器輸出的干擾信號積分能量可表示為hHXHjam·Xjamh或xHHHjamHjamx,期望失配濾波后的積分副主能量比以及干信比足夠小,以獲得可靠的目標(biāo)檢測性能,則可考慮如下約束函數(shù):

    g1(x,h)=hHXHWXh+hHXHjamXjamh|hHx|2(6)

    ISRJ是通過對發(fā)射信號進(jìn)行間歇性采樣所得,那么ISRJ信號經(jīng)過失配濾波后的脈壓峰值主要取決于零距離時延處脈壓值的大小,即干擾信號的脈壓峰值為hHxjam。為獲得可靠的目標(biāo)檢測性能,則應(yīng)保證干擾信號歸一化脈壓峰值足夠低,因此可考慮如下約束函數(shù):

    g2(x,h)=|hHxjam|2|hHx|2(7)

    此外,在發(fā)射波形與失配濾波器的聯(lián)合優(yōu)化問題中,為能夠?qū)SRJ起到有效的干擾抑制作用,會導(dǎo)致與ISRJ信號相對應(yīng)的發(fā)射信號片段經(jīng)過失配濾波器處理后只能獲得極小的信號處理增益,即與匹配濾波機(jī)制相比,會存在一定程度的SNRL,則在該聯(lián)合優(yōu)化問題中,需要將SNRL問題也考慮在內(nèi)。在失配濾波體制下,SNRL[16可表示為

    SNRL=10lgh2x2|hHx|2(8)

    式中:lg(·)表示取對數(shù)操作;·表示矢量的L2范數(shù)。那么在本文的聯(lián)合優(yōu)化問題中,SNRL可直接通過以下約束函數(shù)來控制:

    g3(x,h)=h2x2|hHx|2(9)

    式中:h=[h1,h2,…,hNT;x=[x1,x2,…,xNT。此外,假設(shè)接收濾波器滿足能量限制:hHh=N,且發(fā)射信號滿足恒模約束:|xn|=1, n=1,2,…,N。

    綜上所述,發(fā)射序列和接收濾波器的聯(lián)合優(yōu)化問題可以表示為多目標(biāo)優(yōu)化問題:

    minx,h α1g1(x,h)+α2g2(x,h)+α3g3(x,h)s.t. hHh=N

    |xn|=1, n=1,2,…,N(10)

    式中:α1、α2、α3均為可調(diào)節(jié)權(quán)重因子,可以通過增大某權(quán)重因子來增加其所對應(yīng)的約束函數(shù)在損失函數(shù)中所占的比重,從而獲得該方面性能較好的優(yōu)化結(jié)果。需要指出的是上述優(yōu)化問題無法使算法的抗干擾性能與信號處理增益性能同時達(dá)到最優(yōu)。為能夠?qū)SRJ起到有效的干擾抑制作用,會導(dǎo)致與ISRJ信號相對應(yīng)的信號片段經(jīng)過失配濾波器處理后只能獲得極小的信號處理增益,即存在一定的SNRL,而所允許的SNRL越小,優(yōu)化結(jié)果的抗干擾性能就會越差;相反,要想獲得較好的抗干擾性能,就必須犧牲一定的信號處理增益;那么,就需要在抗干擾性能與信號處理增益之間進(jìn)行折中。

    2 基于MM方法的模型求解

    觀察式(10)所表示的優(yōu)化問題,可知同時對發(fā)射信號和接收濾波器進(jìn)行優(yōu)化設(shè)計(jì)是十分困難的,因此本文采用交替迭代的方式分別對發(fā)射波形和接收濾波器進(jìn)行優(yōu)化,具體的優(yōu)化過程可以分為兩步:首先固定發(fā)射波形x(l),優(yōu)化失配濾波器h;然后再將優(yōu)化得到的失配濾波器h(l+1)固定,優(yōu)化發(fā)射波形x,重復(fù)進(jìn)行上述交替迭代過程直至算法收斂。那么,本文所提出的聯(lián)合優(yōu)化設(shè)計(jì)問題可表述為以下兩大步:

    h(l+1)=arg min g(x(l),h)(11)

    x(l+1)=arg min g(x,h(l+1))(12)

    式中:x(l)為第l次迭代所得的發(fā)射波形;h(l+1)和x(l+1)分別表示第l+1次迭代后所得到的接收濾波器和發(fā)射波形。

    2.1 非匹配濾波器的優(yōu)化問題

    固定第l次迭代所得到的發(fā)射波形x(l),優(yōu)化失配濾波器h,此時優(yōu)化問題式(10)可具體表示為

    minhα1hH(X(l)HWX(l)h+hH(X(l)jamHX(l)jamh|hHx(l)|22|hHx(l)jam|2|hHx(l)|23h2x(l)2|hHx(l)|2

    s.t. hHh=N(13)

    式中:x(l)為第l次迭代時所得到的發(fā)射波形,X(l)為第l次迭代所得到的發(fā)射序列矩陣,x(l)jam為與x(l)相對應(yīng)的間歇采樣干擾信號。而當(dāng)前的雷達(dá)發(fā)射信號x(l)已知且其滿足恒模約束,那么優(yōu)化問題式(13)可以寫為

    minhα1hH(X(l)HWX(l)h+hH(X(l)jamHX(l)jamhhHx(l)(x(l)Hh+α2hHx(l)jam(x(l)jamHhhHx(l)(x(l)Hh+α3NhHhhHx(l)(x(l)Hh

    s.t. hHh=N(14)

    為方便分析,定義兩個新函數(shù):

    g(x(l),h)=c(x(l),h)d(x(l),h)(15)

    e(x(l),h)=c(x(l),h)-g(x(l),h(l))d(x(l),h)(16)

    式中:

    c(x(l),h)=hH(α1((X(l)HWX(l)+(X(l)jamHX(l)jam)+

    α2x(l)jam(x(l)jamH3NIN)h

    d(x(l),h)=hHx(l)(x(l)Hh

    R(l)11((X(l)HWX(l)+(X(l)jamHX(l)jam)+α2x(l)jam(x(l)jamH3NIN

    R(l)2=g(x(l),h(l))x(l)(x(l)H

    根據(jù)式(16)對函數(shù)e(x(l),h)的定義知,當(dāng)h=h(l)時有:

    e(x(l),h(l))=c(x(l),h(l))-g(x(l),h(l))d(x(l),h(l))=0

    那么若第l+1次迭代所得到的接收濾波器h(l+1)能夠使函數(shù)e(x(l),h)值減小,則有:

    e(x(l),h(l+1))=c(x(l),h(l+1))-g(x(l),h(l))d(x(l),h(l+1))lt;

    e(x(l),h(l))=0(17)

    而其中d(x(l),h(l+1))是x(l)作為發(fā)射信號、h(l+1)作為接收濾波器時發(fā)射信號經(jīng)失配濾波器處理后的主瓣能量,其值始終大于0,則整理式(17)可得:

    c(x(l),h(l+1))d(x(l),h(l+1))=g(x(l),h(l+1))lt;g(x(l),h(l)

    即若e(x(l),h(l+1))lt;e(x(l),h(l))則一定滿足g(x(l),h(l+1))lt;g(x(l),h(l)),反之也成立。那么優(yōu)化問題式(14)可以轉(zhuǎn)化為

    minh hHR(l)1h-h(huán)HR(l)2h

    s.t. hHh=N(18)

    需要指出的是優(yōu)化問題式(18)的最優(yōu)解h*不一定是原優(yōu)化問題式(14)的最優(yōu)解,但一定可以保證在本次迭代中使原優(yōu)化問題式(14)中的目標(biāo)函數(shù)值減小,令R(l)=R(l)1-R(l)2。為保證R(l)正定,可使用對角加載方法令R~(l)=R(l)+λIN,其中λ=max{0,-λmin(R(l))}。顯然R~(l)為半正定矩陣,則根據(jù)文獻(xiàn)[19-20]中的引理1可知,可使用MM原理對優(yōu)化問題進(jìn)行簡化,利用MM法求解非凸優(yōu)化問題的收斂性[21-22,已有大量文獻(xiàn)進(jìn)行了分析。利用MM原理簡化并略去無關(guān)項(xiàng)與常數(shù)項(xiàng)后,最終可將優(yōu)化問題式(14)轉(zhuǎn)化為

    minh 2Re(h(R~(l)-tr(R~(l))IN)h(l)

    s.t. hHh=N(19)

    式中:Re(·)表示取實(shí)部運(yùn)算;tr(·)表示取矩陣的跡;IN為單位陣,其階次與接收濾波器長度一致,也為N。

    求解優(yōu)化問題式(19)可得,當(dāng)前條件下最優(yōu)接收濾波器的解析表達(dá)式為

    h(l+1)=-Nw(l)2w(l)(20)

    式中:N為接收濾波器的長度,且有:

    w(l)=(R~(l)-tr(R~(l))IN)h(l)(21)

    2.2 發(fā)射波形的優(yōu)化問題

    固定第(l+1)次迭代所得的接收濾波器h(l+1),優(yōu)化發(fā)射波形x,則優(yōu)化問題式(10)可具體表述為

    minx α1xH(H(l+1)HWH(l+1)x+xH(H(l+1)jamHH(l+1)jamx|xHh(l+1)|22|xHh(l+1)jam|2|xHh(l+1)|23h(l+1)2x2|xHh(l+1)|2

    s.t. |x(n)|=1, n=1,2,…,N(22)

    同樣地,利用接收濾波器的能量約束條件并使用MM原理23-26對優(yōu)化問題式(22)進(jìn)行簡化可得:

    式中:

    求解優(yōu)化問題式(23),當(dāng)前條件下最優(yōu)發(fā)射信號可表示為

    x(l+1)=-ej·arg(v(l+1))(25)

    式中:arg(·)表示取相位操作,且有:

    v(l+1)=(T~(l+1)-tr(T~(l+1))IN)x(l)(26)

    3 算法加速

    3.1 平方迭代算法加速優(yōu)化問題求解

    在上述推導(dǎo)中使用了MM原理提高了算法的運(yùn)行速度,在本節(jié)中將利用平方迭代加速算法[27-28來對聯(lián)合優(yōu)化問題進(jìn)一步加速。為方便表示,分別定義第l次迭代過程中求解失配濾波器與發(fā)射波形的函數(shù)為

    Fh(x(l),h(l))=-Nw(l)2w(l)(27)

    Fx(x(l),h(l+1))=-ej·arg(v(l+1))(28)

    式中:w(l)與v(l+1)的定義分別見式(21)、式(26)。

    使用MM原理對第2節(jié)中所構(gòu)造的優(yōu)化模型進(jìn)行求解的算法流程如算法1所示,但使用MM原理進(jìn)行優(yōu)化時算法的收斂速度與所定義的優(yōu)化函數(shù)性質(zhì)密切相關(guān)29,為進(jìn)一步提高算法的收斂速度,可用平方迭代加速方法來加速聯(lián)合優(yōu)化算法的收斂速度。根據(jù)上述分析,本文所提出的發(fā)射波形和接收濾波器聯(lián)合優(yōu)化抗ISRJ信號算法如算法2所示。

    3.2 算法復(fù)雜度分析

    結(jié)合分析式(20)、式(21)以及式(25)、式(26)可知,迭代更新發(fā)射波形x和接收濾波器h的計(jì)算復(fù)雜度主要是由矩陣R~(l)和T~(l+1)的計(jì)算更新引入的:而根據(jù)式(17)與式(24)可知要更新上述兩個矩陣會涉及矩陣之間的乘法運(yùn)算,其計(jì)算復(fù)雜度為O(N3);為減小算法計(jì)算復(fù)雜度,本文中直接計(jì)算R~(l)h(l)與T~(l+1)x(l),對發(fā)射波形x和接收濾波器h進(jìn)行更新,從而避免了由于直接進(jìn)行矩陣R~(l)和T~(l+1)更新引入的矩陣與矩陣之間的乘法運(yùn)算,此時其計(jì)算復(fù)雜度僅來自矩陣和向量之間的乘法運(yùn)算,計(jì)算復(fù)雜度為O(N2)。此外,本文根據(jù)發(fā)射波形的恒模約束以及接收濾波器的能量約束信息推出了矩陣R~(l)和T~(l+1)的跡滿足的關(guān)系式并利用MM原理對矩陣跡的求解進(jìn)行了簡化,從而避免了直接計(jì)算矩陣跡需要更新矩陣R~(l)和T~(l+1)而引入的矩陣之間的乘法運(yùn)算,因此本文算法的計(jì)算復(fù)雜度僅為O(N2)。本文推導(dǎo)得出矩陣R~(l)和T~(l+1)的跡滿足以下關(guān)系式:

    式中:D為間歇采樣轉(zhuǎn)發(fā)的占空比。

    下面給出證明:

    首先,利用矩陣跡[30的性質(zhì):tr(A+B)=tr(A)+tr(B),可得:

    tr(R~(l))=tr(R(l)1)-tr(R(l)2)+tr(λIN)(31)

    將R(l)1與R(l)2的表達(dá)式代入式(31),并利用矩陣跡[30的另一定理:tr(kA)=k·tr(A)可得:

    tr(R~(l))=α1(tr((X(l)HWX(l))+tr((X(l)jamHX(l)jam))+

    α2·tr(x(l)jam(x(l)jamH)+α3N·tr(IN)+λ·tr(IN)-

    g(x(l),h(l))tr(x(l)(x(l)H)(32)

    而本文中發(fā)射信號滿足恒模約束,則有:

    tr(x(l)(x(l)H)=N(33)

    tr(x(l)jam(x(l)jamH)=DN(34)

    此外,結(jié)合式(2)及矩陣W的定義得:

    tr((X(l)HWX(l))=N(N-1)(35)

    結(jié)合矩陣Xjam的定義式(4)分析可知矩陣(X(l)jamHX(l)jam的第n個對角元素為pnNi=1xix*i=Npn,那么其跡可以表示為

    tr((X(l)jamHX(l)jam)=N∑Nn=1pn(36)

    式中:pn為用于間歇采樣的向量p=[p1,p2,…,pNT的第n個元素,那么∑Nn=1pn可以用DN進(jìn)行表示,則有:

    tr((X(l)jamHX(l)jam)=DN2(37)

    至此,可以得出矩陣R~(l)的跡為

    tr(R~(l))=α1N((N-1)+DN)+α2DN+α3N2

    g(x(l),h(l))N+λN

    矩陣T~(l+1)跡的推導(dǎo)與矩陣R~(l)的推導(dǎo)類似,相同部分不再贅述,需要指出的是接收濾波器滿足能量約束hHh=N,則有:

    tr(h(l+1)jam(h(l+1)jamH)≤tr(h(l+1)(h(l+1)H)=N

    tr((H(l+1)jamHH(l+1)jam)≤tr((H(l+1)HH(l+1))=N2

    其中,當(dāng)且僅當(dāng)間歇采樣占空比D=100%時等號成立。那么,就可以得出:

    tr(T~(l+1))≤α1N((N-1)+N)+α2N+α3N2

    g(x(l),h(l+1))N+λN

    為方便分析,定義兩個變量:

    ε=α1N((N-1)+DN)+α2DN+α3N2

    g(x(l),h(l))N+λN(38)

    γ=α1N((N-1)+N)+α2N+α3N2

    g(x(l),h(l+1))N+λN(39)

    結(jié)合關(guān)系式(29)、式(30)以及矩陣R~(l)和T~(l+1)的半正定性,可知存在以下關(guān)系式:

    R~(l)pεIN

    T~(l+1)pγIN

    那么根據(jù)MM原理可得,利用ε和γ分別替換矩陣R~(l)和T~(l+1)跡的準(zhǔn)確值仍可保證算法迭代過程中目標(biāo)函數(shù)是單調(diào)下降的16,且可有效地避免由于直接計(jì)算矩陣R~(l)和T~(l+1)而引入的矩陣乘法運(yùn)算,從而減小算法的計(jì)算復(fù)雜度。

    4 仿真結(jié)果分析

    4.1 發(fā)射序列和失配濾波器性能仿真

    本節(jié)將對本文所提聯(lián)合設(shè)計(jì)方法的性能進(jìn)行驗(yàn)證。具體仿真參數(shù)的設(shè)置如表1所示,且算法的終止條件設(shè)為相鄰兩次迭代的相對誤差小于等于10-5,即當(dāng)

    g(x(l+1),h(l+1))-g(x(l),h(l))g(x(l),h(l)

    ≤10-5時算法終止迭代。為了量化算法的性能,定義第l次迭代時的歸一化損失函數(shù)29

    gNor(l)=10lg

    g(x(l),h(l))g(x(0),h(0)

    式中:x(0),h(0)分別為用隨機(jī)相位編碼序列初始化的發(fā)射波形和接收濾波器初值。

    圖1為歸一化損失函數(shù)在算法迭代過程中的變化曲線,從中可以看出,隨著迭代次數(shù)的增加,歸一化損失函數(shù)值逐漸減小且最終趨于穩(wěn)定,說明本文所提的聯(lián)合優(yōu)化算法可以有效地解決式(10)所示的優(yōu)化問題。圖2(a)、圖2(b)分別為所設(shè)計(jì)的發(fā)射信號與接收濾波器的零多普勒相關(guān)函數(shù)以及ISRJ與接收濾波器的零多普勒相關(guān)函數(shù)。從圖中可以看出,發(fā)射信號和失配濾波器脈壓輸出的峰值旁瓣比約為-29.6 dB,有效起到了抑制旁瓣能量的作用,且干擾信號和失配濾波器脈壓輸出的峰值可降低至-31.28 dB,說明本文算法所設(shè)計(jì)的發(fā)射序列與失配濾波器有較好的抗干擾性能。

    表2分別給出了不同信號長度下發(fā)射信號、ISRJ與接收濾波器相關(guān)函數(shù)的旁瓣峰值。從表中可以看出,隨著信號長度N的增大,旁瓣峰值和干擾脈壓峰值均逐漸降低,即所設(shè)計(jì)的發(fā)射序列和接收濾波器的相關(guān)性能與抗干擾性能逐漸提升,這是由于信號長度的增加,會使得模型優(yōu)化設(shè)計(jì)的自由度增大,從而能夠得到抗干擾性能更好的發(fā)射波形與接收濾波器。但需要指出的是,信號長度的增加會導(dǎo)致算法的運(yùn)算量增大,進(jìn)而導(dǎo)致算法的實(shí)時性變差(不同信號長度下該算法的具體運(yùn)行時間在表3中給出),因此在實(shí)際應(yīng)用中應(yīng)根據(jù)實(shí)際需求對算法實(shí)時性與抗干擾性能進(jìn)行綜合考慮,以選取合適的信號長度。

    4.2 算法實(shí)時性分析

    本節(jié)將對MM法、平方迭代加速M(fèi)M法以及文獻(xiàn)[16]算法的運(yùn)行時間進(jìn)行分析對比,驗(yàn)證本文算法的實(shí)時性。本節(jié)的仿真參數(shù)與表1中的設(shè)置一致。圖3中給出了MM法與平方迭代加速M(fèi)M法的損失函數(shù)收斂曲線,可以看出,隨著迭代次數(shù)的增加,兩種算法下的損失函數(shù)值均逐漸減小且最終均可收斂,但相比于MM法,加速M(fèi)M法的運(yùn)行時間明顯縮短,即其可以在不改變目標(biāo)函數(shù)收斂性的前提下加快算法收斂速度,同時也驗(yàn)證了本文算法2的可行性。

    為進(jìn)一步分析說明本文算法的實(shí)時性,本節(jié)在不同信號長度下對本文算法與文獻(xiàn)[16]算法的收斂時間進(jìn)行了對比。需要指出的是,為消除樣本偶然性帶來的影響,以下仿真中令兩種算法在各個信號長度下分別進(jìn)行50次蒙特卡羅實(shí)驗(yàn),并對各次運(yùn)行時間取平均,作為該參數(shù)體制下算法的運(yùn)行時間。表3記錄了在不同發(fā)射信號長度下,本文算法與文獻(xiàn)[16]算法的運(yùn)行時間??梢钥闯?,隨著信號與濾波器長度N的增加,兩種算法的運(yùn)行時間均明顯增加,這是由N增大會導(dǎo)致算法的計(jì)算量增大而造成的,但在長度N與迭代終止條件均相同的條件下本文所提算法的收斂速度明顯快于文獻(xiàn)[16]算法的收斂速度,即本文算法具有更好的實(shí)時性(本文中所有仿真均在同一平臺進(jìn)行)。

    4.3 算法抗干擾性能仿真

    本節(jié)將在單個目標(biāo)、單個干擾機(jī)的場景下對優(yōu)化得到的發(fā)射波形和接收濾波器進(jìn)行抗干擾性能分析。將雷達(dá)、目標(biāo)以及干擾機(jī)的位置分別設(shè)置為(0,0,0)m、(0,8 000,0)m、(0,6 000,0)m,假設(shè)該仿真場景下干擾機(jī)的工作模式為直接轉(zhuǎn)發(fā)干擾,并將信干比與信噪比分別設(shè)定為-10 dB、15 dB,其余仿真參數(shù)設(shè)置與表1一致。圖4給出了在發(fā)射序列和接收濾波器初始值相同的情況下分別使用兩種算法優(yōu)化所得的序列,并將其作為發(fā)射信號與接收濾波器后回波的失配濾波結(jié)果以及回波經(jīng)傳統(tǒng)匹配濾波后的結(jié)果,從中可以看出經(jīng)匹配濾波處理后干擾信號峰值高于目標(biāo)信號峰值,說明當(dāng)干擾信號較強(qiáng)時,傳統(tǒng)的匹配濾波處理已無法抑制ISRJ,而本文算法與文獻(xiàn)[16]方法下的輸出干擾峰值均在-10 dB之下,且本文算法具有更低的干擾信號脈壓峰值,抗干擾性能更為可靠。

    表4給出了在上述仿真場景下本文算法與文獻(xiàn)[16]方法的具體抗干擾性能對比,可以看出在SNRL基本相同的情況下,文獻(xiàn)[16]算法的輸出信干比略優(yōu)于本文算法,但本文算法的輸出干擾峰值以及輸出積分主副比性能明顯優(yōu)于文獻(xiàn)[16]算法,說明本文算法能夠以較小的SNRL為代價換取對ISRJ信號的有效抑制作用,進(jìn)而保證可靠的抗干擾性能。

    為進(jìn)一步分析本文算法的抗干擾性能,圖5給出了兩種算法下干擾機(jī)位于不同位置時輸出干擾峰值的變化。其中,雷達(dá)與目標(biāo)的位置設(shè)置不變,而干擾機(jī)的位置以20 m為間隔從(0,6 000,0)m變化至(0,7 980,0)m,在每個位置下均進(jìn)行100次蒙特卡羅實(shí)驗(yàn),并取均值作為該位置下的干擾峰值以減弱隨機(jī)噪聲的影響。從中可以看出,兩種算法在干擾機(jī)與目標(biāo)距離較近的情況下均可對干擾信號起到較好的抑制作用,且本文算法下的干擾峰值總體上低于文獻(xiàn)[16]方法。

    4.4 參數(shù)敏感性分析

    對于ISRJ而言,干擾機(jī)的輸出主要由間歇采樣重復(fù)周期、采樣占空比、采樣轉(zhuǎn)發(fā)起始位置以及轉(zhuǎn)發(fā)次數(shù)等關(guān)鍵參數(shù)所決定。在實(shí)際應(yīng)用中可基于雷達(dá)回波對上述干擾機(jī)參數(shù)進(jìn)行估計(jì):首先對雷達(dá)回波進(jìn)行均勻分段,對分段信號分別進(jìn)行最優(yōu)階次分?jǐn)?shù)階傅里葉變換后進(jìn)行拼接;然后改變信號分段長度,重復(fù)進(jìn)行上述操作,即可得到多組拼接信號;最后從多個拼接信號中尋找復(fù)雜度最低的信號,分析該拼接信號相鄰波峰的距離以及波峰分布的周期性,即可得到上述參數(shù)的估計(jì)值[31。

    從前文的分析中可知,在獲知干擾機(jī)參數(shù)的前提下,本文算法對ISRJ可以起到有效的抑制作用。但在實(shí)際應(yīng)用中上述關(guān)鍵參數(shù)的估計(jì)往往具有一定的誤差,為進(jìn)一步分析本文所提算法的性能,本節(jié)將對該算法的3個關(guān)鍵參數(shù):間歇采樣占空比、間歇采樣重復(fù)周期以及間歇采樣轉(zhuǎn)發(fā)起始位置進(jìn)行敏感性分析。

    4.4.1 間歇采樣占空比敏感性分析

    選取第4.3節(jié)中針對間歇采樣重復(fù)周期為10 μs、間歇采樣占空比D=0.25的情況下優(yōu)化得到的結(jié)果作為雷達(dá)發(fā)射序列和接收濾波器,具體仿真環(huán)境設(shè)置與第4.3節(jié)相同。圖6(a)~圖6(c)分別給出了間歇采樣重復(fù)周期固定為10 μs,改變實(shí)際間歇采樣占空比時本文算法與文獻(xiàn)[16]方法的輸出干擾峰值、積分主副比以及輸出信干比的變化曲線。

    從圖6(a)中可以看出,當(dāng)間歇采樣的實(shí)際占空比低于預(yù)設(shè)占空比時,兩種算法下的干擾峰值均在-10 dB之下,都能夠保證可靠的抗干擾性能,但此時本文算法對占空比有相對較低的敏感度;當(dāng)實(shí)際占空比高于預(yù)設(shè)值時,隨著偏離程度的增大,兩種算法下的干擾信號脈壓峰值均逐步增大,算法抗干擾性能變差,即當(dāng)實(shí)際占空比高于預(yù)設(shè)占空比時,兩種優(yōu)化算法對占空比都十分敏感,但可以看出,在相同的占空比偏離程度下本文算法有相對較低的干擾峰值。從圖6(b)中可以看出,隨著實(shí)際間歇采樣占空比與預(yù)設(shè)值之間偏離程度的增大,兩種算法下非匹配濾波輸出的積分主副比均有所減小,但在相同的占空比偏離程度下,使用本文算法優(yōu)化得到的序列可獲得較高的輸出積分主副比,對旁瓣的抑制作用優(yōu)于文獻(xiàn)[16]。從圖6(c)中可以看出,當(dāng)實(shí)際間歇采樣占空比低于預(yù)設(shè)值時,兩種算法下的輸出信干比基本一致,但當(dāng)實(shí)際采樣占空比高于預(yù)設(shè)值時,本文算法下的輸出信干比低于文獻(xiàn)[16]。綜合上述3項(xiàng)主要性能比較,可知本文算法對間歇采樣占空比的敏感度總體上低于文獻(xiàn)[16]方法。

    在實(shí)際應(yīng)用中,通常無法獲得間歇采樣占空比的準(zhǔn)確值,所估計(jì)出的占空比與實(shí)際占空比之間會存在一定程度的偏差,可以通過增大聯(lián)合優(yōu)化時的預(yù)設(shè)占空比來使得算法能夠在更大的實(shí)際間歇采樣占空比區(qū)間內(nèi)保證可靠的抗干擾性能,即獲得更大的對間歇采樣占空比不敏感的區(qū)間。但需要指出的是,預(yù)設(shè)間歇采樣占空比增大后,干擾信號與發(fā)射信號之間的相似性增加,會使得本文算法優(yōu)化設(shè)計(jì)的自由度減小,進(jìn)而導(dǎo)致算法的抗干擾性能有所減弱,且會導(dǎo)致SNRL增大。同時,從圖6(b)、圖6(c)中可以看出,隨著預(yù)設(shè)占空比正向偏離實(shí)際占空比程度的增大,輸出積分主副比與輸出信干比均逐漸減小,即會導(dǎo)致算法對干擾信號以及旁瓣的抑制作用變差,從而使得目標(biāo)檢測性能有一定程度的下降。因此,在實(shí)際應(yīng)用中應(yīng)根據(jù)間歇采樣占空比的估計(jì)值在一定范圍內(nèi)適當(dāng)增大占空比預(yù)設(shè)值。

    4.4.2 間歇采樣重復(fù)周期敏感性分析

    本節(jié)仍選取第4.3節(jié)仿真條件下優(yōu)化得到的結(jié)果分別作為發(fā)射波形和接收濾波器,具體仿真環(huán)境仍與第4.3節(jié)的設(shè)置相同。圖7(a)~圖7(c)分別給出了將實(shí)際間歇采樣占空比固定為D=0.25,改變實(shí)際間歇采樣重復(fù)周期時兩種算法下的輸出干擾峰值、積分主副比以及輸出信干比的變化曲線。從圖7(a)中可以看出,隨著實(shí)際間歇采樣重復(fù)周期與聯(lián)合優(yōu)化時預(yù)設(shè)采樣重復(fù)周期值之間偏離程度的增大,兩種算法下的干擾信號脈壓峰值均逐步上升,使算法抗干擾性能變差。當(dāng)偏離程度達(dá)到約±7%時,干擾信號峰值與目標(biāo)信號脈壓峰值基本相當(dāng),已無法滿足該場景下的目標(biāo)檢測需求,但在相同的間歇采樣重復(fù)周期偏離程度下本文算法有相對較低的干擾峰值,抗干擾性能更為可靠。從圖7(b)可以看出,隨著實(shí)際間歇采樣重復(fù)周期與預(yù)設(shè)值之間偏離程度的增大,兩種算法下的輸出積分主副比均顯著減小,但本文算法下的輸出積分主副比總體上高于文獻(xiàn)[16]方法,對旁瓣能量有較好的抑制作用。從圖7(c)中可以看出,兩種算法下的輸出信干比基本一致,且隨著實(shí)際間歇采樣周期與預(yù)設(shè)值之間偏離程度的增大,兩種算法下的輸出信干比均逐漸減小,當(dāng)偏離程度達(dá)到±10%時,輸出信干比接近-10 dB,基本無法發(fā)揮對干擾信號的抑制作用。綜合對以上3項(xiàng)重要性能指標(biāo)的討論,可知兩種算法對于間歇采樣重復(fù)周期的敏感度均較高。因此,在實(shí)際應(yīng)用中應(yīng)提高對間歇采樣重復(fù)周期的估計(jì)精度,以保證算法抗干擾性能得以發(fā)揮。

    4.4.3 間歇采樣起始位置敏感性分析

    本節(jié)仍選取第4.3節(jié)中的優(yōu)化結(jié)果作為雷達(dá)發(fā)射波形和接收濾波器,且具體仿真環(huán)境與第4.3節(jié)的設(shè)置相同。圖8(a)~圖8(c)分別給出了改變實(shí)際間歇采樣起始位置時本文算法與文獻(xiàn)[16]方法下的輸出干擾峰值、積分主副比以及輸出信干比的變化曲線。

    從圖8(a)中可以看出,當(dāng)間歇采樣轉(zhuǎn)發(fā)的實(shí)際起始位置與預(yù)估起始位置有相對±8%以內(nèi)的偏差時,兩種算法下的干擾峰值均在-10 dB之下,都能夠保證可靠的抗干擾性能,但可以看出,在相同的起始位置偏離程度下本文算法有相對較低的輸出干擾峰值。從圖8(b)中可以看出,隨著實(shí)際采樣轉(zhuǎn)發(fā)起始位置與預(yù)設(shè)值之間偏離程度的增大,兩種算法下非匹配濾波輸出的積分主副比均有所減小,但總體來看使用本文算法優(yōu)化得到的序列可獲得相對較高的輸出積分主副比,對旁瓣的抑制作用優(yōu)于文獻(xiàn)[16]方法。從圖8(c)中可以看出,隨著實(shí)際采樣轉(zhuǎn)發(fā)起始位置相對于預(yù)設(shè)起始位置偏離程度的增大,兩種算法下的輸出信干比均逐漸減小。綜合上述3項(xiàng)關(guān)鍵性能的比較,可以得出兩種算法對采樣轉(zhuǎn)發(fā)起始位置偏移量的敏感度均較高,這也對干擾機(jī)關(guān)鍵參數(shù)的估計(jì)精度提出了較高的要求。

    5 結(jié)束語

    本文提出了一種發(fā)射波形與接收濾波器聯(lián)合設(shè)計(jì)的抗ISRJ算法,仿真結(jié)果表明,本文方法能夠以較小的SNRL為代價換取可靠的抗干擾性能,同時有效地縮短了算法運(yùn)行時間,具有較好的實(shí)時性。需要指出的是,本文算法對干擾信號的間歇采樣重復(fù)周期有較高的敏感度,這就要求對間歇采樣重復(fù)周期的估計(jì)足夠精確,否則會導(dǎo)致本文算法的抗干擾性能下降。下一步的研究重點(diǎn)是尋找精準(zhǔn)估計(jì)間歇采樣重復(fù)周期的方法,且尋找更為高效的優(yōu)化算法。在保證可靠抗干擾性能的同時提升算法的實(shí)時性,也是未來研究的重點(diǎn)。

    參考文獻(xiàn)

    [1]GONG S X, WEI X Z, LI X. ECCM scheme against interrupted sampling repeater jammer based on time-frequency analysis[J]. Journal of Systems Engineering and Electronics, 2014, 25(6): 996-1003.

    [2]ZHANG Y, WEI Y S, YU L. Interrupted sampling repeater jamming recognition and suppression based on phase-coded signal processing[J]. Signal Processing, 2022, 198(14): 108596.

    [3]HE J Y, CHENG Z Y, HE Z S. Joint design of detecting and screening waveforms against interrupted sampling repeater jamming with inaccurate jamming knowledge[J]. Digital Signal Processing, 2022, 130(3): 103679.

    [4]LU L, GAO M G. A truncated matched filter method for interrupted sampling repeater jamming suppression based on jamming reconstruction[J]. Remote Sensing, 2021, 14(1): 97.

    [5]ZHOU C, LIU Q H, CHEN X L. Parameter estimation and suppression for DRFM-based interrupted sampling repeater jammer[J]. IET Radar, Sonar amp; Navigation, 2018, 12(1): 56-63.

    [6]LIU Y X, ZHANG Q, LIU Z D, et al. An anti-jmming method against interrupted sampling repeater jamming based on compressed sensing[J]. Sensors, 2022, 22(6): 2239.

    [7]FENG D J, XU L T, PAN X Y, et al. Jamming wideband radar using interrupted-sampling repeater[J].IEEE Trans.on Aerospace and Electronic Systems, 2017, 53(3): 1341-1354.

    [8]SUN J W, WANG C, SHI Q Z, et al. Interrupted sampling repeater signal based on phase modulation[C]∥Proc.of the International Conference on Microwave and Millimeter Wave Technology, 2021.

    [9]CHEN J, XU S Y, ZOU J W, et al. Interrupted-sampling repeater jamming suppression based on stacked bidirectional gated recurrent unit network and infinite training[J]. IEEE Access, 2019, 7: 107428-107437.

    [10]沈健君, 劉曉斌, 艾小鋒, 等. 宙斯盾雷達(dá)典型信號處理方法仿真分析[J]. 太赫茲科學(xué)與電子信息學(xué)報, 2022, 20(3): 228-234.

    SHEN J J, LIU X B, AI X F, et al. Simulation analysis of typical signal’s processing methods of Aegis radar[J]. Journal of Te-rahertz Science and Electronic Information Technology, 2022, 20(3): 228-234.

    [11]ZHONG S Z, HUANG X, WANG H J, et al. Anti-intermittent sampling repeater jamming waveform design based on immune genetics[C]∥Proc.of the IEEE International Conference on Power, Intelligent Computing and Systems, 2019: 553-559.

    [12]張彥, 葉春茂, 魯耀兵, 等. 雙曲調(diào)頻波形抑制間歇采樣轉(zhuǎn)發(fā)干擾效果分析[J]. 系統(tǒng)工程與電子技術(shù), 2021, 43(11): 3169-3176.

    ZHANG Y, YE C M, LU Y B, et al. Hyperbolic frequency modulated waveform for interrupted sampling repeater jamming suppression[J]. Systems Engineering and Electronics, 2021, 43(11): 3169-3176.

    [13]李思文, 王國宏, 張亮, 等. 一種強(qiáng)干擾背景下間歇采樣轉(zhuǎn)發(fā)干擾抑制算法[J]. 現(xiàn)代電子技術(shù), 2021, 44(15): 1-6.

    LI S W, WANG G H, ZHANG L, et al. Suppression algorithm for interrupted-sampling forwarding jamming under background of stronginte rference[J]. Modern Electronics Technique,2021, 44(15): 1-6.

    [14]ZHOU C, LIU F F, LIU Q H. An adaptive transmitting scheme for interrupted sampling repeater jamming suppression[J]. Sensors, 2017, 17(11): 2480.

    [15]何金陽, 程子揚(yáng), 李綻蕾, 等. 低多普勒敏感的抗間歇采樣干擾波形設(shè)計(jì)方法[J]. 系統(tǒng)工程與電子技術(shù), 2023, 45(5): 1333-1341.

    HE J Y, CHENG Z Y, LI Z L, et al. Low Doppler sensitive waveform design method against interrupted sampling repeater jamming[J]. Systems Engineering and Electronics, 2023, 45(5): 1333-1341.

    [16]周凱, 何峰, 粟毅. 一種快速抗間歇采樣轉(zhuǎn)發(fā)干擾波形和濾波器聯(lián)合設(shè)計(jì)算法[J]. 雷達(dá)學(xué)報, 2022, 11(2): 264-277.

    ZHOU K, HE F,SU Y. Fast algorithm for joint waveform and filter design against interrupted sampling repeater jamming[J]. Journal of Radar, 2022, 11(2): 264-277.

    [17]ZHOU K, LI D X, SU Y, et al. Joint design of transmit waveform and mismatch filter in the presence of interrupted sampling repeater jamming[J]. IEEE Signal Processing Letters, 2020, 27: 1610-1614.

    [18]ZHOU K, LI D X, QUAN S N, et al. SAR waveform and mismatched filter design for countering interrupted-sampling repeater jamming[J]. IEEE Tran.on Geoscience and Remote Sensing, 2022, 60: 5214514.

    [19]SONG J X, BABU P, PALOMAR D P. Optimization methods for designing sequences with low autocorrelation sidelobes[J]. IEEE Trans.on Signal Processing, 2015, 63(15): 3998-4009.

    [20]SONG J X, BABU P, PALOMAR D P. Sequence design to minimize the weighted integrated and peak sidelobe levels[J]. IEEE Trans.on Signal Processing, 2016, 64(8): 2051-2064.

    [21]GEIPING J, MOELLER M. Composite optimization by nonconvex majorization-minimization[J]. SIAM Journal on Imaging Sciences, 2018, 11(4): 2494-2528.

    [22]SUN Y, BABU P, PALOMAR D P. Majorization-minimization algorithms in signal processing, communications, and machine learning[J]. IEEE Trans.on Signal Processing, 2017, 65(3): 794-816.

    [23]HUNTER D R, LANGE K. A tutorial on MM algorithms[J]. The American Statistician, 2004, 58(1): 30-37.

    [24]何金陽, 程子揚(yáng), 何子述. 抗間歇采樣轉(zhuǎn)發(fā)干擾的認(rèn)知恒模波形設(shè)計(jì)方法[J]. 系統(tǒng)工程與電子技術(shù), 2021, 43(9): 2448-2456.

    HE J Y, CHENG Z Y, HE Z S. Cognitive constant modulus waveform design method against interrupted sampling repeater jamming[J]. Systems Engineering and Electronics, 2021, 43(9): 2448-2456.

    [25]WU L H, ZOU Y B. Majorization-minimization method for elliptic localization in the absence of transmitter position[J]. Sensors, 2022, 23(1). DOI:10.3390/s23010373.

    [26]GONG S Q, XING C W, LAU V K N, et al. Majorization-minimization aided hybrid transceivers for MIMO interference channels[J]. IEEE Trans.on Signal Processing, 2020, 68: 4903-4918.

    [27]ZHAO L C, SONG J X, BABU P, et al. A unified framework for low autocorrelation sequence design via majorization-minimization[J]. IEEE Trans.on Signal Processing, 2017, 65(2): 438-453.

    [28]VARADHAN R, ROLAND C.Simple and globally convergent methods for accelerating the convergence of any EM algorithm[J]. Scandinavian Journal of Statistics, 2008, 35(2): 335-353.

    [29]王福來, 龐晨, 殷加鵬, 等. 一種多普勒容忍的抗間歇采樣轉(zhuǎn)發(fā)干擾恒?;パa(bǔ)波形和接收濾波器聯(lián)合設(shè)計(jì)方法[J]. 雷達(dá)學(xué)報, 2022, 11(2): 278-288.

    WANG F L, PANG C, YIN J P, et al. Joint design of Doppler-tolerant complementary sequences and receiving filters against interrupted sampling repeater jamming[J]. Journal of Radars, 2022, 11(2): 278-288.

    [30]馬菊俠. 關(guān)于矩陣跡的運(yùn)算[J]. 陜西師范大學(xué)學(xué)報(自然科學(xué)版), 2003(S1): 39-41.

    MA J X. Operations about matrix trace[J]. Journal of Shaanxi Normal University (Natural Science Edition), 2003(S1): 39-41.

    [31]于雷, 位寅生, 孟蕓蕓, 等. 一種間歇采樣轉(zhuǎn)發(fā)干擾的參數(shù)辨識方法[P]. 中國: CN 110794374B, 2022-05-17.

    YU L, WEI Y S, MENG Y Y, et al. A parameter identification method for intermittent sampling and forwarding interference[P]. China: CN 110794374B, 2022-05-17.

    作者簡介

    劉 茜(1999—),女,博士研究生,主要研究方向?yàn)槟婧铣煽讖嚼走_(dá)成像、波形優(yōu)化及其抗干擾。

    戴奉周(1978—),男,教授,博士,主要研究方向?yàn)橹悄芨兄c信號處理、電磁超材料、毫米波雷達(dá)。

    欧美一区二区精品小视频在线| 99久久精品一区二区三区| 在线播放国产精品三级| 亚洲va日本ⅴa欧美va伊人久久| 特级一级黄色大片| 亚洲成av人片在线播放无| 脱女人内裤的视频| 欧美日本视频| 99国产精品一区二区三区| 精品久久久久久久久久久久久| 欧美成人免费av一区二区三区| 老司机午夜十八禁免费视频| 国产麻豆成人av免费视频| 一级作爱视频免费观看| 男插女下体视频免费在线播放| 日韩欧美精品免费久久 | 国产一区二区在线观看日韩 | 欧美在线一区亚洲| 国产伦人伦偷精品视频| 久久久久久大精品| 久久久国产精品麻豆| 久久精品影院6| 一二三四社区在线视频社区8| 国产免费av片在线观看野外av| 久久99热这里只有精品18| 男人舔女人下体高潮全视频| 亚洲成人精品中文字幕电影| 国产男靠女视频免费网站| 国产精品影院久久| 久久久久久九九精品二区国产| 在线天堂最新版资源| 搞女人的毛片| 国产亚洲精品综合一区在线观看| 欧美大码av| 老司机福利观看| 欧美黄色淫秽网站| 欧美色欧美亚洲另类二区| 美女 人体艺术 gogo| 国产麻豆成人av免费视频| 国产日本99.免费观看| 国产av一区在线观看免费| 免费人成在线观看视频色| a级一级毛片免费在线观看| 国产在线精品亚洲第一网站| 亚洲在线观看片| 级片在线观看| 午夜久久久久精精品| 国产精品久久电影中文字幕| 久久久国产成人免费| 成年人黄色毛片网站| 亚洲精品影视一区二区三区av| 丝袜美腿在线中文| 国产精品爽爽va在线观看网站| 欧美性猛交黑人性爽| 国产三级黄色录像| 国产蜜桃级精品一区二区三区| 欧美又色又爽又黄视频| 久久欧美精品欧美久久欧美| 精品国产亚洲在线| 中文字幕高清在线视频| 精品久久久久久久人妻蜜臀av| 国产激情偷乱视频一区二区| 啦啦啦韩国在线观看视频| 精品国内亚洲2022精品成人| 欧美最黄视频在线播放免费| 亚洲av中文字字幕乱码综合| 免费av观看视频| 午夜日韩欧美国产| 色综合婷婷激情| 在线观看舔阴道视频| 内射极品少妇av片p| 欧美+日韩+精品| 老司机深夜福利视频在线观看| 亚洲av五月六月丁香网| 女同久久另类99精品国产91| 日韩精品青青久久久久久| 婷婷亚洲欧美| 亚洲一区二区三区不卡视频| 日韩欧美在线乱码| 久久精品国产综合久久久| 一个人看视频在线观看www免费 | 久久精品影院6| 女生性感内裤真人,穿戴方法视频| 天堂√8在线中文| 乱人视频在线观看| 亚洲av日韩精品久久久久久密| 黄色视频,在线免费观看| av视频在线观看入口| 久久久久九九精品影院| 亚洲人成电影免费在线| 精品熟女少妇八av免费久了| xxxwww97欧美| 夜夜夜夜夜久久久久| 亚洲久久久久久中文字幕| 最近最新中文字幕大全电影3| 欧美最黄视频在线播放免费| 波多野结衣巨乳人妻| 757午夜福利合集在线观看| 国产一级毛片七仙女欲春2| 亚洲专区国产一区二区| 女同久久另类99精品国产91| 中文字幕av成人在线电影| 神马国产精品三级电影在线观看| www日本黄色视频网| 国产精品久久久久久精品电影| 国产av不卡久久| 一个人免费在线观看电影| 久久久久国内视频| 女警被强在线播放| 色av中文字幕| 午夜精品一区二区三区免费看| 亚洲av中文字字幕乱码综合| 99久久99久久久精品蜜桃| 又黄又爽又免费观看的视频| 国产亚洲精品久久久久久毛片| 一区二区三区高清视频在线| 男人舔奶头视频| 热99在线观看视频| 国产免费av片在线观看野外av| АⅤ资源中文在线天堂| av欧美777| 五月玫瑰六月丁香| 亚洲狠狠婷婷综合久久图片| 国产高清有码在线观看视频| 亚洲精品影视一区二区三区av| 亚洲欧美日韩高清专用| 午夜激情福利司机影院| av片东京热男人的天堂| 国产一级毛片七仙女欲春2| 手机成人av网站| 精品不卡国产一区二区三区| 在线观看免费视频日本深夜| 亚洲专区国产一区二区| 精品人妻偷拍中文字幕| 伊人久久精品亚洲午夜| 久久九九热精品免费| 午夜精品在线福利| 国产亚洲欧美98| 99久久成人亚洲精品观看| 99国产综合亚洲精品| 欧美+亚洲+日韩+国产| 亚洲色图av天堂| 人妻久久中文字幕网| 成年人黄色毛片网站| 欧美zozozo另类| 综合色av麻豆| 国产私拍福利视频在线观看| 在线观看一区二区三区| 欧美又色又爽又黄视频| 精品不卡国产一区二区三区| 亚洲人成网站在线播放欧美日韩| 久久久久久九九精品二区国产| 色噜噜av男人的天堂激情| 亚洲熟妇中文字幕五十中出| 欧美性猛交╳xxx乱大交人| 日本 欧美在线| 热99re8久久精品国产| 久久久久精品国产欧美久久久| 成年女人毛片免费观看观看9| 午夜激情欧美在线| 欧美日韩亚洲国产一区二区在线观看| 国产综合懂色| 亚洲中文日韩欧美视频| 少妇高潮的动态图| 不卡一级毛片| 亚洲真实伦在线观看| 黄色女人牲交| 深爱激情五月婷婷| 亚洲成av人片免费观看| 精品久久久久久成人av| 脱女人内裤的视频| 激情在线观看视频在线高清| 国内精品美女久久久久久| 欧美又色又爽又黄视频| 国产精品香港三级国产av潘金莲| 九色成人免费人妻av| 日本 欧美在线| 日韩亚洲欧美综合| 成年女人毛片免费观看观看9| 久久国产精品人妻蜜桃| 一本综合久久免费| 国产99白浆流出| 国产不卡一卡二| 久久精品综合一区二区三区| 老鸭窝网址在线观看| 欧美丝袜亚洲另类 | 亚洲专区国产一区二区| 国产免费一级a男人的天堂| www日本黄色视频网| 亚洲无线在线观看| 禁无遮挡网站| 亚洲男人的天堂狠狠| 免费一级毛片在线播放高清视频| 小说图片视频综合网站| 日本五十路高清| 十八禁网站免费在线| 亚洲黑人精品在线| 少妇人妻精品综合一区二区 | 少妇丰满av| 噜噜噜噜噜久久久久久91| 毛片女人毛片| 亚洲精品乱码久久久v下载方式 | 成人特级av手机在线观看| 男人和女人高潮做爰伦理| 亚洲av电影不卡..在线观看| 99久久精品一区二区三区| 悠悠久久av| 岛国视频午夜一区免费看| 日韩高清综合在线| 国产成人啪精品午夜网站| 日韩欧美在线乱码| 我的老师免费观看完整版| 国产欧美日韩一区二区三| 久久6这里有精品| 老司机午夜十八禁免费视频| 99久久精品一区二区三区| 日本黄大片高清| 久久久久久久久久黄片| 嫁个100分男人电影在线观看| 丰满人妻熟妇乱又伦精品不卡| 好看av亚洲va欧美ⅴa在| 亚洲专区中文字幕在线| 波多野结衣高清作品| 久久九九热精品免费| 两个人看的免费小视频| 变态另类成人亚洲欧美熟女| 国产伦精品一区二区三区视频9 | 国产真实伦视频高清在线观看 | 99久久99久久久精品蜜桃| 欧美三级亚洲精品| 村上凉子中文字幕在线| 国产免费av片在线观看野外av| 深爱激情五月婷婷| 午夜免费观看网址| 欧美成人性av电影在线观看| 成人无遮挡网站| 国产欧美日韩一区二区精品| 欧美黄色淫秽网站| 国产av一区在线观看免费| 国产av麻豆久久久久久久| 精品国产超薄肉色丝袜足j| 日韩精品青青久久久久久| 国产午夜福利久久久久久| 国产三级中文精品| 最新中文字幕久久久久| a级毛片a级免费在线| 91久久精品国产一区二区成人 | 天天躁日日操中文字幕| 又爽又黄无遮挡网站| 蜜桃久久精品国产亚洲av| 亚洲精品乱码久久久v下载方式 | 88av欧美| 亚洲狠狠婷婷综合久久图片| 热99在线观看视频| 国产蜜桃级精品一区二区三区| 久久香蕉国产精品| av视频在线观看入口| 亚洲最大成人手机在线| 国产综合懂色| 午夜日韩欧美国产| 少妇高潮的动态图| 国产真人三级小视频在线观看| 日韩欧美 国产精品| 日韩欧美一区二区三区在线观看| 女生性感内裤真人,穿戴方法视频| 91在线精品国自产拍蜜月 | 中文资源天堂在线| 日本熟妇午夜| 成人av一区二区三区在线看| 99精品欧美一区二区三区四区| 黄色视频,在线免费观看| 国产一区二区三区在线臀色熟女| 免费搜索国产男女视频| 国产精品三级大全| tocl精华| 久久精品国产99精品国产亚洲性色| 国产又黄又爽又无遮挡在线| 国产精品野战在线观看| 免费看光身美女| 色哟哟哟哟哟哟| 美女免费视频网站| 99久久成人亚洲精品观看| 51午夜福利影视在线观看| 一本一本综合久久| 中出人妻视频一区二区| 最近最新中文字幕大全电影3| 99久久99久久久精品蜜桃| bbb黄色大片| 人妻久久中文字幕网| 国产v大片淫在线免费观看| 亚洲最大成人中文| 色哟哟哟哟哟哟| 国产激情偷乱视频一区二区| 久久国产精品影院| 男女午夜视频在线观看| 日韩欧美三级三区| 母亲3免费完整高清在线观看| 尤物成人国产欧美一区二区三区| 欧美日本视频| 悠悠久久av| 国产三级黄色录像| 亚洲五月婷婷丁香| 最新在线观看一区二区三区| 人妻久久中文字幕网| 国产精品1区2区在线观看.| 国产99白浆流出| 免费电影在线观看免费观看| 久久性视频一级片| 蜜桃久久精品国产亚洲av| 一级毛片高清免费大全| 午夜福利在线在线| 亚洲18禁久久av| 男插女下体视频免费在线播放| 久久人妻av系列| 99久久九九国产精品国产免费| 日韩精品青青久久久久久| 久久婷婷人人爽人人干人人爱| 久久精品91无色码中文字幕| 很黄的视频免费| 制服人妻中文乱码| 三级国产精品欧美在线观看| 18禁黄网站禁片免费观看直播| 有码 亚洲区| 国产97色在线日韩免费| 国产亚洲精品综合一区在线观看| 亚洲国产精品sss在线观看| 亚洲av成人av| 高清毛片免费观看视频网站| 成人一区二区视频在线观看| 欧美区成人在线视频| 国产一区二区在线观看日韩 | 舔av片在线| 国产亚洲av嫩草精品影院| 免费观看精品视频网站| 香蕉久久夜色| 最好的美女福利视频网| 超碰av人人做人人爽久久 | 日本一二三区视频观看| 99精品在免费线老司机午夜| 欧美三级亚洲精品| 他把我摸到了高潮在线观看| 99久久无色码亚洲精品果冻| 国产亚洲欧美98| 看黄色毛片网站| 成人三级黄色视频| 综合色av麻豆| 国产精品亚洲av一区麻豆| 99精品在免费线老司机午夜| 99国产极品粉嫩在线观看| 十八禁网站免费在线| 亚洲最大成人手机在线| 亚洲国产欧美网| 波多野结衣高清无吗| 欧美在线一区亚洲| 婷婷丁香在线五月| 一本综合久久免费| 久久久成人免费电影| 亚洲五月天丁香| 脱女人内裤的视频| 国产亚洲精品av在线| 久久精品91蜜桃| 国产av不卡久久| 久久久国产成人精品二区| 人人妻人人看人人澡| 国内精品美女久久久久久| 国产精品女同一区二区软件 | 岛国在线观看网站| 12—13女人毛片做爰片一| 嫩草影院精品99| av女优亚洲男人天堂| 天堂网av新在线| 999久久久精品免费观看国产| 国产私拍福利视频在线观看| 亚洲精华国产精华精| 99热这里只有精品一区| 精品国产亚洲在线| 99在线人妻在线中文字幕| 成年免费大片在线观看| 色综合婷婷激情| 国产亚洲精品综合一区在线观看| 亚洲一区二区三区色噜噜| 久久中文看片网| 女人十人毛片免费观看3o分钟| 精品无人区乱码1区二区| 黑人欧美特级aaaaaa片| 久久精品国产综合久久久| 免费av观看视频| 好男人电影高清在线观看| eeuss影院久久| 久久亚洲真实| 精品久久久久久久久久免费视频| 欧美区成人在线视频| 日本精品一区二区三区蜜桃| 亚洲第一欧美日韩一区二区三区| 黄色日韩在线| 国产伦精品一区二区三区视频9 | 亚洲最大成人手机在线| 亚洲人成网站高清观看| 人妻夜夜爽99麻豆av| 国产精品精品国产色婷婷| 国产精品久久久久久久久免 | 岛国在线免费视频观看| 久久久久性生活片| 国产精品 欧美亚洲| 我的老师免费观看完整版| 天堂√8在线中文| 日本撒尿小便嘘嘘汇集6| 国产乱人伦免费视频| 国产欧美日韩一区二区三| 在线观看免费午夜福利视频| 日韩大尺度精品在线看网址| 国产乱人视频| 免费人成视频x8x8入口观看| 特大巨黑吊av在线直播| 少妇熟女aⅴ在线视频| 国产精品一区二区三区四区久久| 国产精华一区二区三区| 91九色精品人成在线观看| 美女大奶头视频| 国产精品嫩草影院av在线观看 | 一级黄片播放器| 波野结衣二区三区在线 | 国产男靠女视频免费网站| 日本 欧美在线| www.熟女人妻精品国产| 手机成人av网站| 欧美性猛交黑人性爽| 成人鲁丝片一二三区免费| 久久香蕉精品热| 波多野结衣高清作品| 狂野欧美激情性xxxx| av福利片在线观看| 乱人视频在线观看| 日本撒尿小便嘘嘘汇集6| 婷婷亚洲欧美| 岛国视频午夜一区免费看| 男人和女人高潮做爰伦理| 99精品在免费线老司机午夜| 国产av在哪里看| 高清毛片免费观看视频网站| 色综合婷婷激情| 国产91精品成人一区二区三区| 国产亚洲精品久久久com| 99热这里只有是精品50| 啦啦啦免费观看视频1| 一区二区三区免费毛片| 国产亚洲精品av在线| 18禁黄网站禁片免费观看直播| 日本 欧美在线| 一本综合久久免费| 国产成年人精品一区二区| 久9热在线精品视频| 国产成人影院久久av| www国产在线视频色| 亚洲成人久久性| 免费看美女性在线毛片视频| 脱女人内裤的视频| 久久国产乱子伦精品免费另类| 国产美女午夜福利| 特级一级黄色大片| 身体一侧抽搐| 成人精品一区二区免费| 他把我摸到了高潮在线观看| 欧美日韩中文字幕国产精品一区二区三区| 特大巨黑吊av在线直播| 日本免费一区二区三区高清不卡| 亚洲av一区综合| 黑人欧美特级aaaaaa片| 日韩欧美国产一区二区入口| 久久午夜亚洲精品久久| 国产亚洲精品综合一区在线观看| 亚洲国产欧洲综合997久久,| 中文字幕精品亚洲无线码一区| 黄色日韩在线| 深夜精品福利| 国产精品亚洲av一区麻豆| 亚洲七黄色美女视频| 国内精品久久久久精免费| 亚洲内射少妇av| 1024手机看黄色片| 成年女人看的毛片在线观看| 国产高清有码在线观看视频| 成年女人毛片免费观看观看9| 三级国产精品欧美在线观看| 欧美性感艳星| aaaaa片日本免费| 成年女人毛片免费观看观看9| 国产精品99久久99久久久不卡| 免费电影在线观看免费观看| 亚洲av电影不卡..在线观看| 国产男靠女视频免费网站| 国产精品久久视频播放| 欧美日韩国产亚洲二区| 搡女人真爽免费视频火全软件 | 国产在视频线在精品| av在线蜜桃| 一个人免费在线观看电影| 午夜激情福利司机影院| 午夜久久久久精精品| 男人舔奶头视频| 老鸭窝网址在线观看| 国产激情偷乱视频一区二区| 国产蜜桃级精品一区二区三区| 久99久视频精品免费| 国产精品98久久久久久宅男小说| 国产精品一区二区三区四区久久| bbb黄色大片| 人妻夜夜爽99麻豆av| 最近最新中文字幕大全免费视频| netflix在线观看网站| 久久99热这里只有精品18| 91麻豆av在线| 国产高清视频在线播放一区| 亚洲五月婷婷丁香| 波多野结衣巨乳人妻| 深夜精品福利| 欧美一级毛片孕妇| 午夜影院日韩av| 国产亚洲av嫩草精品影院| 亚洲专区国产一区二区| 日韩欧美三级三区| 女警被强在线播放| 久久久久免费精品人妻一区二区| 国产成人福利小说| 女人被狂操c到高潮| 丁香六月欧美| 精品久久久久久,| 久久精品91无色码中文字幕| 91九色精品人成在线观看| 国产精品影院久久| 久久久国产成人免费| 欧美不卡视频在线免费观看| 一个人免费在线观看电影| 啦啦啦韩国在线观看视频| 中文字幕人妻丝袜一区二区| 国产97色在线日韩免费| av天堂在线播放| 久久久久久久精品吃奶| 成人高潮视频无遮挡免费网站| 长腿黑丝高跟| 亚洲人成网站高清观看| 搡老岳熟女国产| 高清毛片免费观看视频网站| 亚洲精品一卡2卡三卡4卡5卡| 国产黄a三级三级三级人| 免费在线观看日本一区| 亚洲精品日韩av片在线观看 | 久久草成人影院| 最好的美女福利视频网| 黄片大片在线免费观看| 欧美区成人在线视频| 99在线人妻在线中文字幕| 免费在线观看日本一区| 国产主播在线观看一区二区| 国产成人啪精品午夜网站| 国产av在哪里看| 夜夜夜夜夜久久久久| 两个人的视频大全免费| 久久久久久大精品| 亚洲av成人精品一区久久| 99热只有精品国产| 成熟少妇高潮喷水视频| 少妇丰满av| 可以在线观看的亚洲视频| 欧美乱色亚洲激情| 制服丝袜大香蕉在线| 成人欧美大片| 少妇的逼好多水| 国产精品永久免费网站| 久久久久久久久大av| 国产真实乱freesex| 老熟妇乱子伦视频在线观看| 神马国产精品三级电影在线观看| aaaaa片日本免费| 国产一区二区三区视频了| 亚洲av免费高清在线观看| 手机成人av网站| 中亚洲国语对白在线视频| 欧美高清成人免费视频www| 久久久精品欧美日韩精品| 日韩大尺度精品在线看网址| 亚洲熟妇熟女久久| 国产又黄又爽又无遮挡在线| 久久精品国产清高在天天线| 成人鲁丝片一二三区免费| 91av网一区二区| 一级黄片播放器| 亚洲av成人不卡在线观看播放网| 亚洲最大成人手机在线| 国产免费一级a男人的天堂| 岛国在线免费视频观看| 欧美日韩精品网址| 国产精品爽爽va在线观看网站| 91在线观看av| 久久精品国产清高在天天线| 久久久久久久亚洲中文字幕 | 又紧又爽又黄一区二区| 嫁个100分男人电影在线观看| 免费看光身美女| 老鸭窝网址在线观看| 国产精品乱码一区二三区的特点| 亚洲黑人精品在线| 中文字幕av在线有码专区| 国产精品自产拍在线观看55亚洲| 啦啦啦观看免费观看视频高清| 99精品欧美一区二区三区四区| 欧美极品一区二区三区四区| 亚洲精品在线美女| 亚洲av免费高清在线观看| 亚洲va日本ⅴa欧美va伊人久久| 搡老妇女老女人老熟妇| 亚洲欧美日韩高清在线视频| 成年女人永久免费观看视频| 国产aⅴ精品一区二区三区波| 国产av不卡久久|