亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        連續(xù)凸逼近加權(quán)自相關(guān)波形優(yōu)化方法

        2023-09-28 07:12:48宋青青
        雷達(dá)與對(duì)抗 2023年3期
        關(guān)鍵詞:副瓣旁瓣波形

        張 勇,宋青青

        (1.南京六九零二科技有限公司,南京 210009;2.中國(guó)船舶集團(tuán)有限公司第八研究院,南京 211153)

        0 引 言

        波形設(shè)計(jì)是當(dāng)前最熱門的雷達(dá)系統(tǒng)技術(shù)研究領(lǐng)域之一。有源探測(cè)與通信系統(tǒng)均對(duì)具有良好自相關(guān)特性的波形序列有著十分迫切的需求,例如脈沖壓縮雷達(dá)系統(tǒng)收發(fā)信號(hào)、聲納探測(cè)器編碼序列與數(shù)字通信系統(tǒng)(如GPS接收器或CDMA蜂窩系統(tǒng))的同步、信道估計(jì)的導(dǎo)頻序列,甚至安全系統(tǒng)的密碼學(xué)[1-5]。在實(shí)際應(yīng)用中,由于序列生成硬件的限制(如信號(hào)最大幅度/功率起伏須保持在模數(shù)轉(zhuǎn)換器和功率放大器動(dòng)態(tài)范圍內(nèi)),通常要求發(fā)射波形具有恒模特性。

        近年來,恒模序列(也稱為多相序列)因具有強(qiáng)相關(guān)可塑性、高能量利用率等特點(diǎn)而備受關(guān)注[6-7]。相位編碼序列的獲取途徑主要包括固定計(jì)算(如Barker碼、Frank碼等)與迭代優(yōu)化(如遺傳算法、CAN算法等優(yōu)化架構(gòu))[8]。

        經(jīng)典的Barker碼序列屬于二進(jìn)制編碼序列,可由移位寄存器或綜合邏輯法得到,具有較好的自相關(guān)特性,自相關(guān)峰值旁瓣電平不高于1,但其序列編碼長(zhǎng)度較有限,現(xiàn)階段不超過13。1965年,Golomb和Scholtz開展了廣義Barker碼研究,仍遵循最小化峰值旁瓣比原則,擴(kuò)展了編碼元素集,尋求具有類似Barker序列自相關(guān)特性的相位編碼序列。以此為驅(qū)動(dòng),關(guān)于廣義Barker碼序列的研究不斷深化,序列長(zhǎng)度擴(kuò)展至77。是否存在更長(zhǎng)的廣義Barker序列當(dāng)前仍未得到證實(shí)。除了Barker序列,還存在多相碼簇可由閉式解析式得到,如Frank碼、Golomb序列。隨著序列長(zhǎng)度的增加,該類序列峰值旁瓣比增長(zhǎng)速度與序列長(zhǎng)度的平方根成正比,與Barker碼自相關(guān)特性差異化明顯。

        隨著數(shù)字技術(shù)的發(fā)展,波形碼元集不斷得到擴(kuò)充,碼型序列的復(fù)雜度呈指數(shù)上升態(tài)勢(shì)。面向復(fù)雜波形序列對(duì)強(qiáng)設(shè)計(jì)規(guī)則魯棒性、高算法通用性的需求,迭代優(yōu)化類算法應(yīng)運(yùn)而生,主要包括遺傳算法(Genetic Algorithm,GA)、模擬退火等全局尋優(yōu)、變換域優(yōu)化、乘子類等算法。GA最早是由美國(guó)的John Holland于20世紀(jì)70年代提出,是一種借鑒生物進(jìn)化規(guī)律的結(jié)構(gòu)化搜索全局最優(yōu)解的方法[9]。由于復(fù)雜波形全局尋優(yōu)的可行集巨大,該方法雖然不需要引入復(fù)雜的推導(dǎo)計(jì)算,但對(duì)于高維向量?jī)?yōu)化存在計(jì)算量大、耗時(shí)長(zhǎng)的缺點(diǎn),不利于波形集的快速更新。模擬退火算法(Simulated Annealing,SA)[10]思想最早由N Metropolis等人于1953年提出。在此基礎(chǔ)上,S Kirkpatrick等人于1983年將退火思想引入到組合優(yōu)化領(lǐng)域。它是基于Monte-Carlo迭代求解策略的一種隨機(jī)尋優(yōu)算法,其出發(fā)點(diǎn)是基于物理中固體物質(zhì)的退火過程與一般組合優(yōu)化問題之間的相似性,與遺傳算法一樣采用結(jié)構(gòu)化優(yōu)化架構(gòu),計(jì)算量巨大,亦不適合高維向量?jī)?yōu)化,對(duì)長(zhǎng)度大于等于103量級(jí)的序列優(yōu)化是不實(shí)際的。為了適應(yīng)高維度波形優(yōu)化,文獻(xiàn)[11]和[12]為了最小化自相關(guān)積分旁瓣(Integrated sidelobe level,ISL),采用其“幾乎等效”的頻域表達(dá)為目標(biāo)函數(shù),借助快速傅里葉變化提出了高效的CAN算法及加權(quán)循環(huán)算法(Weight cyclic algorithm new,WeCAN),可實(shí)現(xiàn)長(zhǎng)度不低于104的多相碼序列的快速優(yōu)化,遠(yuǎn)超F(xiàn)rank與Golomb序列長(zhǎng)度?;谏鲜鯥SL等效結(jié)果,以交替優(yōu)化為代表的乘子類算法[13]逐漸應(yīng)用在低ISL波形優(yōu)化中。峰值旁瓣(Peak sidelobe level,PSL)水平是另外一種較為常用的自相關(guān)副瓣考核指標(biāo),具有瞬變的特點(diǎn),相比于以ISL為目標(biāo)函數(shù)的波形設(shè)計(jì)問題,以PSL水平為優(yōu)化目標(biāo)的優(yōu)化問題更加難以駕馭,近端梯度下降法[14]解決了包含多個(gè)等式與不等式共存的非凸優(yōu)化問題,實(shí)現(xiàn)了低峰值旁瓣比的波形優(yōu)化設(shè)計(jì),但該方法對(duì)于子梯度選取具有不確定性,算法結(jié)果較難把握,而且對(duì)于自相關(guān)局部最小化問題無法求解。相比于PSL,以ISL為目標(biāo)函數(shù)更適合快速、可靠波形設(shè)計(jì)的應(yīng)用場(chǎng)合。

        文獻(xiàn)[11]使用了等效加權(quán)積分旁瓣為目標(biāo)函數(shù),無法使其與ISL概念完全一致,最優(yōu)解的ISL無法得到保證。本文以最小化自相關(guān)加權(quán)自相關(guān)副瓣積分為目標(biāo)函數(shù),對(duì)原非凸優(yōu)化問題目標(biāo)函數(shù)經(jīng)過連續(xù)兩次凸逼近,并將模量約束松弛為能量約束,使得原問題轉(zhuǎn)化為凸問題,利用成熟優(yōu)化技術(shù)對(duì)其求解。最后,將最優(yōu)解進(jìn)行幅度強(qiáng)制歸一化處理,獲得適用于飽和功放等器件使用的恒模波形序列,其中可通過調(diào)整權(quán)值分布對(duì)自相關(guān)形狀進(jìn)行有效控制。

        符號(hào)說明:本文采用粗體小寫字母表示矢量,用粗體大寫字母表示矩陣;(·)T、(·)H分別表示矩陣與向量的轉(zhuǎn)置和共軛轉(zhuǎn)置算子;| |、‖ ‖分別表示取輸入?yún)?shù)的絕對(duì)值和l2范數(shù);字母j為虛數(shù)單位;M和M分別表示M維實(shí)數(shù)向量和復(fù)數(shù)向量;M×N和M×N分別表示M×N維實(shí)數(shù)和復(fù)數(shù)矩陣;Tr(·)表示計(jì)算輸入矩陣的跡。

        1 信號(hào)模型與問題構(gòu)建

        (1)

        加權(quán)積分旁瓣電平γ為

        (2)

        式中,w為加權(quán)系數(shù)向量,其向量形式為

        w=[w1,w2,…,wN-1]T

        (3)

        當(dāng)w為全1矩陣時(shí),γ與常規(guī)積分旁瓣電平的含義相同。雷達(dá)、通信系統(tǒng)均需要低自相關(guān)旁瓣電平的波形序列。在實(shí)際電磁輻射系統(tǒng)中發(fā)射能量通常是一定,在此約束下,以最小化積分旁瓣電平為目標(biāo)函數(shù),構(gòu)建優(yōu)化問題為

        (4)

        若要實(shí)現(xiàn)自相關(guān)函數(shù)局部區(qū)域趨零或置零,則須適當(dāng)增加w中相關(guān)區(qū)域?qū)?yīng)的權(quán)值,權(quán)值增加,該區(qū)域?qū)δ繕?biāo)函數(shù)取值影響增大;反之易然。

        2 連續(xù)凸逼近算法

        由于目標(biāo)函數(shù)為非凸的,使用常規(guī)優(yōu)化算法難以解決。為了解決式(4),擬將尋找目標(biāo)函數(shù)凸上界函數(shù),并將其作為目標(biāo)函數(shù),構(gòu)建等效優(yōu)化問題進(jìn)行求解。

        為便于闡釋,采用矩陣表示式(4),則自相關(guān)函數(shù)矩陣形式為

        rk=Tr(ΕkxxH),k=0,…,N-1

        (5)

        式中,Εk∈N×N為第k對(duì)角線全為1的Toeplitz矩陣。

        則式(4)可進(jìn)一步表示為

        (6)

        |xn|=1,n=1,…,N

        由于自相關(guān)函數(shù)r的模值具有對(duì)稱性,因此式(6)可等價(jià)為

        subject toX=xxH

        (7)

        |xn|=1,n=1,…,N

        此時(shí),w維度擴(kuò)展為2N-1,完整表達(dá)式為

        w=[w1-N,…,0,…,wN-1]T

        (8)

        且w0=0,w-k=wk。

        將式(7)目標(biāo)函數(shù)中的矩陣向量化后,經(jīng)推導(dǎo),式(7)目標(biāo)函數(shù)可表示為

        (9)

        令xv=vec(X),ek=vec(Ek),則式(7)可重新整理為

        (10)

        假定

        (11)

        則可將式(10)中的目標(biāo)函數(shù)化簡(jiǎn)為二次型,即

        (12)

        為了最小化非凸目標(biāo)函數(shù),可通過最小化其上邊界凸函數(shù),達(dá)到原目標(biāo)函數(shù)最小化的目的。假定Q=λmax(L)I,因此(Q-L)為半正定矩陣,λmax(L)為L(zhǎng)的最大特征值。結(jié)合L的特殊結(jié)構(gòu),根據(jù)方陣特征值與特征向量的定義,可知

        因此,ek為屬于L的第k個(gè)特征向量,且wk(N-|k|)為特征向量ek對(duì)應(yīng)的特征值??赏ㄟ^比較所有wk(N-|k|),k=0,1,…,N-1值確定L的最大特征值。

        略去第一項(xiàng)與第三項(xiàng)常數(shù)部分,式(12)可重新構(gòu)建為

        (15)

        再將xv0、xv和L均還原為矩陣表達(dá)式,即

        (16)

        (17)

        為求解波形序列x,須將其顯性表示,故將X還原為向量形式,優(yōu)化式(17)可得

        (18)

        式中,

        (19)

        忽略常數(shù)項(xiàng)后,式(18)可轉(zhuǎn)化為

        subject to|xn|=1,n=1,…,N

        (21)

        式中,λg為Γ的最大特征值。

        優(yōu)化式(21)目標(biāo)為線性函數(shù)可示為凸函數(shù),然而由于恒模約束的存在,該問題仍為非凸優(yōu)化問題,難以直接求解,因此將優(yōu)化式(21)的約束進(jìn)行放松,可得

        (22)

        求解優(yōu)化式(22)后,可取其最優(yōu)解的相位向量,重新構(gòu)成恒模向量作為最終波形序列。具體算法流程可見表1。

        表1 連續(xù)凸逼近加權(quán)自相關(guān)波形優(yōu)化算法流程

        3 數(shù)值仿真

        本節(jié)將連續(xù)凸逼近加權(quán)自相關(guān)波形優(yōu)化算法與現(xiàn)有的WeCAN算法進(jìn)行比較,驗(yàn)證所提算法的有效性,開展了針對(duì)常規(guī)低副瓣波形和加權(quán)ISL波形的設(shè)計(jì)性能驗(yàn)證。仿真參數(shù)設(shè)置:非周期波形序列長(zhǎng)度為32、128;最大迭代次數(shù)15 000次;終止常數(shù)10-6;計(jì)算機(jī)配置為Intel(R)Core(TM)i5-7500 CPU @ 3.40 GHz,內(nèi)存4 G。

        本節(jié)從常規(guī)低副瓣波形設(shè)計(jì)入手,給出了恒模波形序列的相位、自相關(guān)函數(shù),并采用數(shù)值分析方法比較了WeCAN算法與連續(xù)凸逼近加權(quán)自相關(guān)波形優(yōu)化算法的收斂性能。進(jìn)一步,針對(duì)加權(quán)自相關(guān)波形優(yōu)化,驗(yàn)證所提算法的可行性,同樣從以上三方面進(jìn)行闡釋比較,證明所提算法的優(yōu)越性。

        (1)常規(guī)低副瓣波形設(shè)計(jì)

        常規(guī)低副瓣波形設(shè)計(jì)時(shí)權(quán)值設(shè)置為相同值。圖1(a)、(b)分別給出了長(zhǎng)度分別為32、128時(shí)WeCAN和連續(xù)凸逼近算法所得常規(guī)低副瓣波形的相位,二者存在較大差異,且類似隨機(jī)噪聲信號(hào)。圖2(a)、(b)分別給出了長(zhǎng)度32、128時(shí)兩算法所得的低副瓣波形自相關(guān)函數(shù),長(zhǎng)度相同時(shí)兩算法對(duì)應(yīng)最優(yōu)波形的自相關(guān)函數(shù)水平幾乎一致。圖3(a)、(b)分別給出了波形序列長(zhǎng)度為32和128時(shí),WeCAN與連續(xù)凸逼近加權(quán)自相關(guān)波形優(yōu)化算法的收斂曲線。從中可看出,WeCAN收斂速度優(yōu)于本文所提算法,此時(shí)頻域等效方法優(yōu)于時(shí)域直接優(yōu)化方法。

        (a)長(zhǎng)度32

        (a)長(zhǎng)度32

        (a)長(zhǎng)度32

        (2)加權(quán)副瓣波形設(shè)計(jì)

        加權(quán)副瓣波形設(shè)計(jì)時(shí),須先確定距離向干擾或雜波所處的距離段,再通過增加自相關(guān)函數(shù)中相應(yīng)延遲部分的權(quán)值,達(dá)到抑制干擾/雜波距離向處理增益的效果。參數(shù)設(shè)置如表2所示。

        表2 加權(quán)自相關(guān)波形優(yōu)化算法參數(shù)設(shè)置

        圖4(a)、(b)分別給出了長(zhǎng)度分別為32、128時(shí)WeCAN和連續(xù)凸逼近加權(quán)自相關(guān)波形優(yōu)化算法所得加權(quán)自相關(guān)副瓣波形的相位。圖5(a)、(b)給出了長(zhǎng)度分別為32、128時(shí)兩算法所得的加權(quán)副瓣波形自相關(guān)函數(shù),長(zhǎng)度相同時(shí)兩算法對(duì)應(yīng)最優(yōu)波形的自相關(guān)函數(shù)形狀一致,但凹口處連續(xù)凸逼近算法對(duì)應(yīng)的電平低于WeCAN,可以獲得更好的距離向干擾、雜波的抑制效果。圖6(a)、(b)分別給出了波形序列長(zhǎng)度為32、128時(shí),WeCAN與連續(xù)凸逼近算法的收斂曲線。從中可看出,本文所提連續(xù)凸逼近加權(quán)自相關(guān)波形優(yōu)化算法優(yōu)于WeCAN算法的收斂速度,且穩(wěn)定收斂時(shí)連續(xù)凸逼近算法目標(biāo)函數(shù)對(duì)應(yīng)的最優(yōu)值低于WeCAN算法,故而相比于WeCAN算法,本文所提算法獲得理想的凹口深度。

        (a)長(zhǎng)度32

        (a)長(zhǎng)度32

        (a)長(zhǎng)度32

        4 結(jié)束語

        本文研究了基于連續(xù)凸逼近加權(quán)自相關(guān)波形優(yōu)化算法,通過凸函數(shù)逼近將不定或非凸的優(yōu)化問題轉(zhuǎn)化為凸問題,并通過強(qiáng)制歸一化模值獲得恒模波形序列,并通過調(diào)整權(quán)值,控制自相關(guān)局部自相關(guān)電平,對(duì)自相關(guān)整體形狀進(jìn)行塑形。對(duì)于加權(quán)副瓣的波形設(shè)計(jì),與已有方法相比具有更快的收斂速度。因此,連續(xù)凸逼近算法可加快現(xiàn)代雷達(dá)波形集的構(gòu)建速度,增加復(fù)雜波形的多樣性,提升現(xiàn)代雷達(dá)距離向干擾與雜波的抑制能力。

        猜你喜歡
        副瓣旁瓣波形
        C波段高增益低副瓣微帶陣列天線設(shè)計(jì)
        基于圓柱陣通信系統(tǒng)的廣義旁瓣對(duì)消算法
        非均勻間距的低副瓣寬帶微帶陣列天線設(shè)計(jì)
        對(duì)《壓力容器波形膨脹節(jié)》2018版新標(biāo)準(zhǔn)的理解及分析
        一種基于線性規(guī)劃的頻率編碼旁瓣抑制方法
        基于LFM波形的靈巧干擾效能分析
        雷達(dá)副瓣跟蹤機(jī)理分析
        基于加權(quán)積分旁瓣最小化的隨機(jī)多相碼設(shè)計(jì)
        基于四項(xiàng)最低旁瓣Nuttall窗的插值FFT諧波分析
        基于ARM的任意波形電源設(shè)計(jì)
        久久久久久无码av成人影院| 韩国美女主播国产三级| 亚洲一区二区三区最新视频 | 成人欧美一区二区三区| 日本一区二区不卡视频| 亚洲中文无码精品久久不卡| 9l国产自产一区二区三区| 国产香蕉一区二区三区在线视频| 大地资源在线观看官网第三页| 成全视频高清免费| 少妇特殊按摩高潮惨叫无码| 日本免费大片一区二区三区| 亚洲爆乳无码精品aaa片蜜桃| 老熟妇仑乱视频一区二区 | 99热在线播放精品6| 日本女优中文字幕亚洲| 精品国产免费一区二区三区香蕉| 人与嘼交av免费| 99久久久精品免费| 亚洲精品一区二区网站| 亚洲欧美色一区二区三区| 女人与牲口性恔配视频免费| 国产一区二区三区av免费观看| 亚洲精品久久区二区三区蜜桃臀| 粗大的内捧猛烈进出视频| 美日韩毛片| 久久夜色精品国产三级| 精品精品国产高清a毛片| 失禁大喷潮在线播放| 国产经典免费视频在线观看| 男奸女永久免费视频网站| 国产一区二区女内射| 中文人成影院| 日本高清成人一区二区三区| 国产成a人亚洲精品无码樱花| 人妻熟妇乱又伦精品视频app| 亚洲黄片高清在线观看| 91精品国产色综合久久| 天天鲁一鲁摸一摸爽一爽| 日韩一二三四精品免费| 国产精品一二三区亚洲|