基于僅相位加權(quán)的寬零陷低副瓣波束賦形方法

周強(qiáng)鋒

(1.中國空空導(dǎo)彈研究院，河南 洛陽 471000;2.航空制導(dǎo)武器航空科技重點(diǎn)實(shí)驗(yàn)室，河南 洛陽 471000)

0 引言

陣列天線寬零陷波束賦形技術(shù)作為有效抑制或消除干擾的一種有效手段[1-2]，已廣泛應(yīng)用于雷達(dá)、通信、衛(wèi)星導(dǎo)航等領(lǐng)域。傳統(tǒng)寬零陷形成技術(shù)[3-6]一般是通過對(duì)陣列天線各單元同時(shí)施加幅度與相位分布激勵(lì)，或者通過對(duì)單元的幅度或位置進(jìn)行擾動(dòng)，實(shí)現(xiàn)寬零陷賦形的目的。由于陣列各單元的激勵(lì)幅度分布是通過饋電網(wǎng)絡(luò)實(shí)現(xiàn)的，且位置一般是固定的，實(shí)現(xiàn)傳統(tǒng)寬零陷形成技術(shù)將增加陣列天線系統(tǒng)的設(shè)計(jì)難度，如需要設(shè)計(jì)復(fù)雜饋電網(wǎng)絡(luò)實(shí)現(xiàn)幅度的改變、增加精確的機(jī)電傳動(dòng)裝置用于實(shí)現(xiàn)單位位置的改變等，缺乏靈活性，工程應(yīng)用較復(fù)雜。

對(duì)于相控陣天線的寬零陷波束賦形，解決以上問題的一種有效方法就是采用只改變饋電相位分布的僅相位加權(quán)方法，在不改變?cè)泄β史峙漯侂娋W(wǎng)絡(luò)和不增加新設(shè)備的情況下，通過相控陣天線移相器來改變各單元的相位分布，實(shí)現(xiàn)方向圖的寬零陷賦形。采用僅相位加權(quán)的零陷賦形方法已得到許多學(xué)者的關(guān)注。Hu[7]利用梯度搜索法實(shí)現(xiàn)了線陣的自適應(yīng)波束形成；Khzmalyan[8]利用一維求最小值法實(shí)現(xiàn)了波束零陷賦形；Gatti[9]利用擬牛頓法對(duì)描述相位加權(quán)分布的多項(xiàng)式的權(quán)值進(jìn)行優(yōu)化，實(shí)現(xiàn)了波束零陷賦形。但是這些方法缺乏對(duì)波束整體副瓣電平的約束且形成的零陷較窄、陡峭，當(dāng)對(duì)干擾方向估計(jì)產(chǎn)生誤差，可能會(huì)導(dǎo)致零陷位置偏出實(shí)際干擾方向，導(dǎo)致抑制干擾失效[10]。

為此，本文提出了一種新的僅相位加權(quán)的陣列天線寬零陷低副瓣波束賦形方法，通過建立對(duì)應(yīng)的優(yōu)化目標(biāo)模型，采用融合混沌搜索技術(shù)的粒子群算法(particle swarm optimization，PSO)進(jìn)行優(yōu)化。利用混沌具有隨機(jī)性、遍歷性等優(yōu)點(diǎn)進(jìn)行種群初始化，增強(qiáng)種群多樣性，并改進(jìn)粒子位置更新策略，對(duì)全局最優(yōu)位置更新采用變尺度混沌搜索技術(shù)，有效結(jié)合粒子群的全局搜索能力和混沌優(yōu)化局部搜索能力，提高了整個(gè)優(yōu)化過程的收斂速度和收斂精度。

1 數(shù)學(xué)模型

假定陣列天線單元方向圖為各向同性，對(duì)于一個(gè)N元等間距排列的線陣，其天線方向圖可表示為

(1)

式中:In為第n個(gè)天線單元的理想激勵(lì)幅度；λ為波長(zhǎng)；d為天線單元排列間距；θ為場(chǎng)空間角；θ0為波束掃描角。

采用僅相位加權(quán)方法方向圖賦形技術(shù)實(shí)質(zhì)是在保持原有指向下，增加一個(gè)用于實(shí)現(xiàn)形成預(yù)定波束形狀的陣列各單元的相位激勵(lì)分布，其基本原理可表示為

(2)

式中:φ=(φ1,φ2,…,φN)，φn為第n個(gè)單元的用于形成預(yù)定波束形狀的相位加權(quán)值。

寬零陷低副瓣方向圖賦形的目的是在強(qiáng)干擾源方向變化的角度區(qū)域內(nèi)設(shè)置寬零陷，并保持整體波束的低副瓣來抑制其他方向的干擾。采用僅相位加權(quán)方法的寬零陷低副瓣波束形成技術(shù)的優(yōu)化目標(biāo)模型可描述如下：

|AF(θm)|,0)，

(3)

式中：M為計(jì)算中定義的角度采樣點(diǎn)數(shù)；θm為第m個(gè)采樣的角度；AF(θm)為理想設(shè)計(jì)的寬零陷方向圖。

由于式(3)的目標(biāo)函數(shù)與參數(shù)變量不構(gòu)成解析式，采用粒子群優(yōu)化算法進(jìn)行求解，獲得形成預(yù)定波束形狀的陣列天線各單元的相位加權(quán)值。

2 寬零陷低副瓣波束賦形的PSO算法

2.1 PSO算法

PSO算法是一種模擬鳥群捕食行為的智能全局優(yōu)化算法[11-12]。在PSO算法中，每個(gè)粒子(個(gè)體)利用自身找到的最優(yōu)解和整個(gè)種群的全局最優(yōu)解提供的信息，來更新自身的速度和位置，實(shí)現(xiàn)尋找最優(yōu)解的目的。

對(duì)于D維參數(shù)變量的優(yōu)化問題，假設(shè)第i個(gè)粒子所經(jīng)歷的最優(yōu)位置Pi=(pi1,pi2,…,piD)，所有粒子經(jīng)歷的全局最優(yōu)位置為Pg=(pg1,pg2,…,pgD)。在第k次循環(huán)過程中，標(biāo)準(zhǔn)PSO算法中第i個(gè)粒子的速度vi=(vi1,vi2,…,viD)和位置Xi=(xi1,xi2,…,xiD)的更新可以描述如下：

(4)

(5)

式中：ω為慣性權(quán)因子；c1和c2為加速度因子；r1和r2為[0,1]間的隨機(jī)數(shù)。慣性權(quán)因子ω常采用自適應(yīng)計(jì)算，即ω=ωmax-k(ωmax-ωmix)/kmax，其中k為當(dāng)前迭代次數(shù)，kmax為最大迭代次數(shù),ωmax和ωmix分別為最大和最小慣性常數(shù)。

2.2 混沌運(yùn)動(dòng)與PSO算法改進(jìn)策略

混沌運(yùn)動(dòng)是存在于非線性確定系統(tǒng)中的一種普遍現(xiàn)象，具有初值敏感性、隨機(jī)性、遍歷性以及自相似分形等非線性動(dòng)力學(xué)特性。利用混沌運(yùn)動(dòng)的這些性質(zhì)進(jìn)行優(yōu)化算法設(shè)計(jì)，是一種有潛力的新技術(shù)，并已經(jīng)用于多種優(yōu)化算法改進(jìn)[13-14]。

典型產(chǎn)生混沌序列的混沌模型有Logistic映射、Tent映射、Kent映射、Sinusoidal映射等，其中Sinusoidal映射已經(jīng)被證明相比于其他映射具有更好的混沌優(yōu)化效果[15]。

Sinusoidal混沌映射函數(shù)的定義如下：

(6)

式中：a=2.3。

為了克服基本PSO算法易出現(xiàn)早熟收斂、收斂緩慢等缺點(diǎn)，將混沌技術(shù)引人PSO算法優(yōu)化過程中，對(duì)種群初始化和粒子位置更新策略進(jìn)行改進(jìn)。

(1) 利用混沌技術(shù)初始化種群

利用混沌運(yùn)動(dòng)的隨機(jī)性和遍歷性優(yōu)點(diǎn)，通過Sinusoidal映射產(chǎn)生一組與優(yōu)化變量相同維數(shù)的混沌變量序列，即隨機(jī)產(chǎn)生一個(gè)D維且每個(gè)分量都在[0,1]間的隨機(jī)數(shù)向量Z1=(z11,z12,…,z1D)。根據(jù)Sinusoidal混沌映射函數(shù)式(6)得到M個(gè)混沌序列Z1，Z2，…,ZM，并將Zi的各個(gè)分量按式(7)映射到優(yōu)化變量的取值范圍[a,b]，利用式(3)計(jì)算每一個(gè)粒子的適應(yīng)度函數(shù)，并從中選取適應(yīng)度最好的Popsize(種群數(shù))個(gè)粒子組成初始種群。利用混沌技術(shù)初始化既保持初始化群時(shí)具有的隨機(jī)性本質(zhì)，又利用混沌技術(shù)提高了種群的多樣性和搜索的遍歷性。

Xi=a+(b-a)Zi.

(7)

(2) 將混沌搜索嵌入粒子位置更新

根據(jù)基本PSO算法的粒子位置更新機(jī)制，對(duì)全局最優(yōu)的粒子位置更新缺乏有效引導(dǎo)，可能會(huì)影響全局的收斂精度與收斂速度。為進(jìn)一步挖掘PSO算法的局部搜索能力和提高優(yōu)化效率，全局最優(yōu)的粒子位置更新采用混沌搜索策略，即每次迭代過程中，當(dāng)粒子位置Xi屬于全局最優(yōu)的位置，對(duì)其位置更新時(shí)，采用如下更新策略：

(8)

式中：Δxij為利用Sinusoidal映射函數(shù)隨機(jī)產(chǎn)生的混沌序列變量；γ為搜索空間控制參數(shù)。

考慮到算法優(yōu)化迭代前期應(yīng)具有較強(qiáng)的全局搜索能力和跳出局部極值的能力，迭代后期算法應(yīng)具有較強(qiáng)的局部精細(xì)搜索能力，采用變尺度混沌優(yōu)化思想，即根據(jù)優(yōu)化進(jìn)程不斷縮小優(yōu)化變量的搜索空間，γ計(jì)算如下：

(9)

式中:α為搜索空間尺度參數(shù)；k為當(dāng)前迭代次數(shù)；kmax為最大迭代次數(shù)。

同時(shí)，為使全局最優(yōu)粒子位置始終保留在更新后的種群中，全局最優(yōu)的粒子位置更新引入競(jìng)爭(zhēng)機(jī)制，當(dāng)更新后的粒子位置適應(yīng)度更好，則保留更新，否則保留原有位置不變，具體公式為：

(10)

由于寬零陷低副瓣方向圖賦形的優(yōu)化變量是帶約束條件的，當(dāng)粒子位置更新完成后，可能會(huì)出現(xiàn)越界，需要強(qiáng)加一個(gè)簡(jiǎn)單的邊界約束。在完成粒子位置更新后，通過如下一個(gè)簡(jiǎn)單的步驟，把粒子位置映射到規(guī)定的區(qū)域內(nèi)：

(11)

2.3 融合混沌技術(shù)的改進(jìn)粒子群算法

對(duì)于采用僅相位加權(quán)方法的寬零陷低副瓣波束賦形優(yōu)化問題，融合混沌技術(shù)的改進(jìn)粒子群算法優(yōu)化步驟如下：

(1) 算法參數(shù)初始化，包括種群規(guī)模Popsize、搜索維數(shù)D、最大慣性常數(shù)ωmax、最小慣性常數(shù)ωmin、加速度因子c1和c2、最大迭代次數(shù)kmax、各粒子的初始速度向量為[-0.2,0.2]內(nèi)的隨機(jī)數(shù)。

(2) 利用Sinusoidal混沌映射初始化種群，計(jì)算每一個(gè)粒子的初始適應(yīng)度，并更新個(gè)體極值Pi和全體極值Pg。

(3) 執(zhí)行新的位置更新策略，利用式(4)和式(5)對(duì)非全局最優(yōu)粒子進(jìn)行位置更新，利用式(8)和式(10)對(duì)屬于全局最優(yōu)粒子的位置進(jìn)行位置更新。

(4) 執(zhí)行邊界約束處理，利用式(11)對(duì)粒子位置進(jìn)行映射到規(guī)定的位置范圍內(nèi)。

(5) 計(jì)算經(jīng)過位置更新后粒子所對(duì)應(yīng)的適應(yīng)度值，更新個(gè)體極值Pi和全體極值Pg。

(6) 判斷當(dāng)前迭代數(shù)k是否超出最大的迭代次數(shù)kmax(停止準(zhǔn)則)。若k

3 實(shí)驗(yàn)結(jié)果

為了驗(yàn)證新方法的有效性，以半波長(zhǎng)等間距排列的60元線陣為例，采用本文改進(jìn)PSO算法(improved PSO,IPSO)和基本粒子群算法(PSO)分別進(jìn)行僅相位加權(quán)寬零陷方向圖賦形。仿真中，種群數(shù)M=60，搜索維數(shù)為D=60，最大迭代次數(shù)為kmax=3 000，最大最小慣性常數(shù)分別為ωmax=0.9和ωmin=0.4，加速度常數(shù)為c1=c2=2，相位加權(quán)值取值空間[a,b]=[0,2π]，搜索空間尺度參數(shù)α=6。

圖1是主波束指向0°時(shí)在(20°，25°) 角度內(nèi)形成寬5°、深-60 dB的零陷且副瓣電平不超過-30 dB的波束賦形仿真結(jié)果。圖中點(diǎn)線是未寬零陷波束賦形的原始方向圖，虛線是寬零陷低副瓣賦形的波形約束，虛點(diǎn)線和實(shí)線分別是采用基本PSO和本文IPSO形成的寬零陷方向圖。采用基本PSO算法未能優(yōu)化出符合要求的寬零陷方向圖，零陷深度和副瓣電平方面均沒有達(dá)到要求，而采用本文IPSO實(shí)現(xiàn)了在預(yù)定位置形成零陷深度-60 dB的寬零陷，且整體副瓣電平不超過-30 dB，其優(yōu)化獲得的激勵(lì)相位加權(quán)值見圖2。圖3是基本PSO和本文IPSO適應(yīng)度收斂曲線對(duì)比?；綪SO收斂速度非常慢，在迭代過程中陷入局部最優(yōu)，未能優(yōu)化出最優(yōu)解；而本文IPSO很好地跳出了局部最優(yōu)，并快速趨于收斂，獲得最優(yōu)解，大大加快了搜索速度。

圖4是主波束指向0°時(shí)在(30°，35°) 角度內(nèi)形成寬5°、深-60 dB的零陷且副瓣電平不超過-30 dB的波束賦形仿真結(jié)果。圖5是本文IPSO優(yōu)化獲得的激勵(lì)相位加權(quán)值。圖6是基本PSO和本文IPSO適應(yīng)度收斂曲線對(duì)比?；綪SO由于陷入局部最優(yōu)且收斂較慢，在(30°，35°) 角度內(nèi)形成的零陷深度約-53 dB；而采用本文IPSO很好地跳出了局部最優(yōu)，在迭代1 200次趨于收斂，實(shí)現(xiàn)了預(yù)定位置形成寬零陷的設(shè)計(jì)要求。仿真結(jié)果表明,采用本文IPSO算法有效地克服了基本粒子群算法易陷入局部最優(yōu)和收斂精度不高的缺點(diǎn)，能夠很好用于形成寬零陷低副瓣波束賦形。

為了進(jìn)一步分析改進(jìn)PSO算法進(jìn)行僅相位加權(quán)寬零陷方向圖賦形優(yōu)化的收斂精度和收斂速度，對(duì)上述給出的2個(gè)優(yōu)化問題各進(jìn)行30次獨(dú)立實(shí)驗(yàn)，結(jié)果如表1所示。

通過表1可以看出，本文改進(jìn)PSO算法用于僅相位加權(quán)寬零陷方向圖賦形優(yōu)化時(shí)，收斂精度高，能夠在3 000次迭代內(nèi)獲得符合要求的最優(yōu)解之一，收斂速度快，算法穩(wěn)健有效，便于工程應(yīng)用。

表1 本文改進(jìn)PSO算法的性能分析Table 1 Performance analysis of IPSO

4 結(jié)束語

針對(duì)寬零陷方向圖賦形問題，本文提出了一種僅改變相位分布的陣列天線寬零陷低副瓣波束賦形優(yōu)化方法；通過建立對(duì)應(yīng)的優(yōu)化目標(biāo)模型，采用融合混沌搜索技術(shù)的粒子群算法進(jìn)行優(yōu)化，獲得陣列天線各單元的饋電相位分布，在不改變各單元饋電幅度分布的前提下形成預(yù)定寬零陷波束形狀，并達(dá)到低副瓣設(shè)計(jì)要求。仿真結(jié)果驗(yàn)證了改進(jìn)的粒子群算法有效結(jié)合粒子群的全局搜索能力和混沌優(yōu)化局部搜索能力，能夠保持低副瓣下形成較深的寬零陷，具有快速收斂和較好穩(wěn)健性的優(yōu)點(diǎn)。