基于同心圓環的大間距子陣級陣列柵瓣抑制

李志剛 王鵬毅

(中國電子科技集團公司第五十四研究所,石家莊 050081)

基于同心圓環的大間距子陣級陣列柵瓣抑制

李志剛 王鵬毅

(中國電子科技集團公司第五十四研究所,石家莊 050081)

適當增加單元間距和應用子陣級陣列有利于降低相控陣的設計成本,但是陣列的稀疏會帶來不期望的柵瓣．為解決該問題,設計了一種無柵瓣陣列:子陣在平面內旋轉,圍繞陣列中心呈同心圓環分布．為保證子陣之間不會交疊,對圓環半徑范圍進行了分析．利用一種改進的窮舉算法,對每一層圓環上的子陣個數以及圓環半徑進行了計算機仿真優化．結果表明,以4×4子陣共1 024個單元為例,所提方法可抑制峰值旁瓣電平約-17.48 dB．

柵瓣抑制;大間距;子陣旋轉;同心圓環陣;稀疏陣

引 言

如何降低相控陣天線的設計成本一直是一個熱門話題,為此可以從兩方面考慮: ①在陣列孔徑尺寸一定條件下，合理減少單元數量,拉開單元間距,可減小互耦獲得高增益,但是會帶來柵瓣[1];②子陣級波束形成,與單元級波束形成相比子陣級波束形成極大地減小了通道數,降低了硬件成本和工程實現的難度[2].

在大間距陣列中,每一個子陣都有自己的柵瓣,整個陣列的柵瓣抑制可以通過以下幾種方法進行抑制:

1) 自適應波束形成[3]．該方法運用自適應加權使子陣方向圖在高旁瓣位置形成凹陷,但是需要實時調整單元權值,增加了計算負擔,甚至數據與權值有可能會失匹配．

2) 虛擬陣列[4]．該方法通過內插虛擬天線實現小間距組陣來抑制柵瓣,然而計算復雜度會隨之增加,且原來白噪聲也會受到污染變成色噪聲,波束形成性能受到影響．

3) 非規則子陣．采用多種子陣,不同子陣單元間距不同[5],或者單元間距相同、子陣規模不同[6]．該方法雖然可以抑制柵瓣,但是并不利于子陣的統一生產及降低制造成本．

4) 子陣間距優化[7-9]．采用同一子陣,對子陣中心位置進行優化．子陣個數越多優化的自由度越大[8],但計算復雜度會隨之成倍增加,并且子陣位置分布的規則性也會降低．

本文基于同一子陣,給出了一種具有對稱結構的無柵瓣布陣方法:子陣圍繞陣列中心等間距地呈同心圓環分布；然后設計了子陣的各項參數,并分析了圓環半徑的取值范圍；最后為最小化峰值旁瓣電平(Peak Side-lobe Level, PSL),對每層圓環的子陣個數和圓環半徑進行了仿真優化．

1 陣列布局

單元間距大于半波長時，陣列的方向圖會有柵瓣．不妨設單元間距d=λ,λ為波長，圖1顯示的是d=λ方形柵格陣列的方向圖．從圖1中可以看出,柵瓣增益與主瓣增益一樣高．

圖1 方形柵格陣列方向圖(d=λ)

文獻[8]提出把陣列分成Ns個相同的矩形柵格子陣,通過調整子陣位置,能夠達到最小柵瓣電平，可以由式(1)估計：

δ=-10lgNsdB,

(1)

子陣個數越多,公式(1)估計越精確．文獻[9]對104個方形柵格子陣進行平移優化,可抑制柵瓣-16.5dB,與公式(1)的理論值20.17dB還差約3.7dB．

除子陣平移外,子陣旋轉[8]可以分散柵瓣的能量,對柵瓣起到抑制效果．但是子陣無規律的旋轉可能會導致子陣之間重疊,因此子陣平移位置和旋轉角度都需要進行合理設計及優化．本文對子陣進行平面旋轉,并將子陣圍繞陣列中心呈同心圓環分布[10],效果如圖2所示．為便于陣列上下(左右)形成差波束,陣列分別關于X軸和Y軸對稱．

圖2 陣列布局效果示意圖 (1 024個單元分為64個子陣)

2 方向圖函數

子陣為Ne×Ne方形柵格結構,單元的水平和垂直間距均為d,那么子陣的邊長為L=Ned,單元的位置(xm,ym)可以求出．設各單元等幅激勵,子陣的方向圖可以表示為[11]

(2)

式中:k=2π/λ，

(3)

表示觀察方向的方向余弦;u0和v0表示掃描方向的方向余弦;與陣面垂直方向俯仰角度φ為0°．

把子陣等間距地放置于同心圓環上,并圍繞陣列中心進行平面旋轉．陣列共有Ns個子陣被放置在Nr個同心圓環上,設第n(n=0,1,…,Nr)環上子陣的個數為Nn,圓環半徑為rn．特別的,n=0表示陣列的中心．那么第n環第m個子陣的中心坐標為

(xn,m,yn,m)=rn(cosαn,m,sinαn,m)

(4)

式中:αn,m=2π(m-1)/Nn+αn,0為子陣的旋轉角度,如圖2所示,αn,0為附加在第n環上的一個可以調節的旋轉角度．

把子陣看成一個單元,整個陣列的方向圖可以用子陣級方向圖表示為[11]

(5)

式中,Fn,m(u,v)表示第n環第m個子陣的方向圖,子陣方向圖可以通過對子陣旋轉和平移后求出,旋轉矩陣為

(6)

需要說明的是:

1) 對于線極化天線,若子陣不一致,旋轉陣列的波束形成會極為不方便;若子陣一致,不同子陣的極化平面可能會不同．

2) 對于圓極化天線,如果子陣被旋轉了,輻射波的相位會變化相同的角度．通過標準的校準技術進行陣列校準,子陣旋轉引起極化的幅相變化可以得到補償[8]．為便于分析,本文假設天線為理想的圓極化天線．

3 子陣位置優化

3.1 基本思路

影響方向圖函數的主要參數包括:子陣內參數,如單元間距d、子陣大小Ne×Ne、子陣個數Ns等;子陣間參數,如圓環上子陣個數Nn、圓環個數Nr、圓環半徑rn以及圓環的旋轉角度αn,0等．上述每一個參數的變化都會導致陣列方向圖的改變,需要優化的參數眾多,計算量龐大．我們優化陣列分布的基本思路是:固定子陣內的參數,優化子陣間的參數,使得峰值旁瓣電平達到最小．需要注意的是,在優化過程中還有一個限制條件,即子陣之間不能發生重疊．

3.2 子陣參數設計

3.3 子陣旋轉角度設計

圓環旋轉角度αn,0的設置．為了保持陣列的對稱結構,αn,0可以為0或者π/Nn．如果αn,0均設為0或者π/Nn,落在對稱軸上的子陣等效于只進行了平移,而沒有進行旋轉．子陣旋轉的目的就是為了把柵瓣能量打散．為了減少等效平移子陣的數量,可以讓某一層圓環的旋轉角度為0,讓相鄰的另一層圓環的旋轉角度為π/Nn,依法交替設置αn,0．

3.4 圓環半徑范圍

上述參數設定后,在一定孔徑范圍內對圓環半徑進行調節,以尋找最小的PSL．在調節過程中,子陣不應存在交疊,從而每一層圓環的半徑存在上邊界和下邊界．

對于最外層圓環,可以預先設定半徑的上邊界為一常數B;對于非最外層圓環,半徑rn的最大值受到外一層圓環半徑的限制,第n環子陣的外切圓半徑應小于第n+1環子陣的內切圓半徑,利用勾股定理有

(7)

化簡公式(7),得到rn的上邊界Bu(rn)為

(8)

當第n環的半徑rn減小到最小時,子陣的邊角緊緊相鄰,子陣朝向陣列中心的邊共同組成一個正多邊形,邊數為Nn．可以計算出該正多邊形內切圓半徑,再加上子陣邊長的一半,即得出rn的下邊界Bl(n)為

(9)

3.5 優化算法

下面優化子陣間參數:圓環半徑rn和每層圓環上子陣個數Nn．

1) 優化流程

針對陣列最佳分布綜合法的研究提出了很多種,如窮舉法、遺傳算法、模擬退火法等[12]．遺傳算法和模擬退火法容易陷入局部最優解,窮舉法可得到最優解,但是計算量龐大、難以實現．

本文對窮舉法進行了改進,對參數rn和Nn分別進行遍歷,流程如圖3所示.

圖3 一種改進的窮舉優化策略流程

步驟如下:

① 確定最外層圓環半徑B,固定子陣內參數,包括:單元間距d、子陣大小Ne×Ne、子陣個數Ns,固定圓環個數Nr;

② 暫時等比例設定rn,優化圓環上子陣個數Nn,為保持陣列軸對稱性,Nn設為偶數,Nn的步進長度為dΔ1=2,每次步進計算當前陣列的PSL;

③ 步驟② 循環結束后,尋找最小PSL對應的Nn;

④ 固定Nn,遍歷半徑rn,步長不妨設為dΔ2=λ/20,每次步進計算當前陣列的PSL;

⑤ 步驟④循環結束后,尋找最小PSL對應的rn,優化過程結束．

2) 運算量分析

(10)

如果采用窮舉法,rn和Nn聯合遍歷的次數由陣列孔徑的大小、子陣個數及層數共同決定．所提優化過程把rn和Nn的聯合遍歷分解成了兩個獨立的一維遍歷過程．

表1給出了一種運算量的示例．

表1 運算量示例(B=22λ,Ns=64,Nr=4)

每次計算PSL的時間由u,v的分辨率設置和計算機運行速度決定,可采用多核處理器并行計算．由表1看出,窮舉法對rn和Nn聯合遍歷,計算量龐大,往往難以實現;所提方法的計算量是兩維優化各自計算量的加和,以一種可行的計算量得到一組次優解．

4 仿真結果與分析

為得到最小的PSL,以MATLAB為仿真平臺,對圓環上的子陣個數Nn和圓環半徑rn進行了優化,其他的參數設置如表2所示．

表2 參數設置

首先,64個子陣分布在4層同心圓環上,圓環半徑rn暫時等比例固定,對每一層圓環上子陣的個數Nn進行優化,優化結果為分別4、12、20和28．對Nn優化后的方向圖如圖4所示,PSL為-14.54 dB．

圖4 rn等比例,對Nn優化后的方向圖

然后,每一層圓環上子陣的個數Nn固定后,對圓環半徑rn進行優化,經優化后結果分別為4.6λ、9.55λ、15.25λ和22λ．最終,優化結果如表3所示．

表3 1 024單元圓環半徑和圓環上

圖5顯示了rn和Nn優化后陣列未掃描時的方向圖．盡管單元之間的最小間距為d=λ,然而由圖4可以看出陣列并不存在柵瓣,PSL的最小值可以達到-17.48 dB,與公式(1)的理論估計結果-18.06 dB相差0.58 dB．

圖6是掃描角度為60°時的方向圖．由圖6看出： 當波束指向掃描時,PSL為-15.85 dB,與沒有掃描時相比,PSL抬高了1.63 dB．

圖5 rn和Nn優化后陣列的方向圖(未掃描時)

圖6 rn和Nn優化后陣列的方向圖(掃描60°時)

圖7 全方位掃描時歸一化峰值旁瓣電平的分布情況

文獻[9]把1 664個單元分成104個方形柵格子陣進行平移優化,未掃描和掃描10°時,方位向旁瓣電平均低于-16.5 dB．為了與文獻[9]進行對比,本文也采用相同的單元個數和子陣個數進行布陣,布陣優化結果如表4和圖8所示．由圖8看出,掃描10°時PSL可以達到-19.64 dB,優于文獻[9]的結果3.14 dB,與公式(1)的理論值相差約0.5 dB．

表4 參數優化結果(1 664個單元分為104個子陣)

圖8 1 664單元優化后方向圖(掃描10°時)

5 結 論

針對大間距、子陣級的無柵瓣布陣問題,給出了一種基于同心圓環的布陣方法,并對布陣參數進行了仿真優化．結果表明,以1 024單元為例,本文對包括柵瓣在內的旁瓣電平抑制結果與理論值相差約0.58 dB,且所設計的陣列呈軸對稱分布,易于形成差波束,為低成本、高增益相控陣天線的工程實現提供了一種理論參考．

[1] MAILOUX R J. Phased array antenna handbook [M]. Boston: Artech Hause, 2005: 80.

[2] KRIVOSHEEV Y V, SHISHLOV A V. Grating lobe suppression in phased arrays composed of identical or similar subarrays [C]//IEEE International Symposium on Phased Array Systems and Technology. Waltham, MA, USA, Oct. 2010: 724-730.

[3] 程乃平, 潘點飛. 大型陣列天線子陣劃分及柵瓣抑制方法[J]. 信號處理, 2014, 30(5): 535-543.

CHENG N P, PAN D F. Subarray partition method and grating lobe suppression for large array anenna [J]. Journal of signal processing, 2014, 30(5): 535-543. (in Chinese)

[4] 汪月清, 劉姜玲, 王小謨, 等. 虛擬陣元應用于平板端射天線的新型組陣方法研究[J]. 電波科學學報, 2014, 29(4): 663-667.

WANG Y Q, LIU J L, WANG X M, et al. Method of making up virtual plate end-fire antenna array [J]. Chinese jounal of radio science, 2014, 29(4): 663-667. (in Chinese)

[5] 朱瑞平, 王朝陽. 一種低柵瓣大瞬時帶寬相控陣天線的設計方法[J]. 電波科學學報, 2008, 23(4): 745-748.

ZHU R P, WANG C Y. Design method for wide band and low grating lobe phased array antennas [J]. Chinese jounal of radio science, 2008, 23(4): 745-748. (in Chinese)

[6] REN Y L, WANG J, LUO T G. Research on grating lobe reduction in limited scan phased array antenna [C]//Cross Strait Quad-Regional Radio Science and Wireless Technology Conference. Chengdu, China, 2013: 230-233.

[7] TOYAMA N. Aperiodic array consisting of subarrays for use in small mobile earth stations [J]. IEEE transactions on antennas and propagation, 2005, 53(6): 2004-2010.

[8] KRIVOSHEEV Y V, SHISHLOV A V, DENISENKO V V. Grating lobe suppression in aperiodic phased array antennas composed of periodic subarrays with large element spacing [J]. IEEE antennas and propagation magazine, 2015, 57(1): 76-85.

[9] 羅天光, 王建, 任永麗, 等. 采用DFP-BFGS優化大間距子陣級非周期八角陣[J]. 現代雷達, 2014, 36(3): 62-69.

LUO T G, WANG J, REN Y L, et al. An octagonal aperiodic array antenna optimizing large spacing’s subarray using DFP-BFGS algorithm [J]. Modern radar, 2014, 36(3): 62-69. (in Chinese)

[10] HAUPT R L. Optimized element spacing for low sidelobe concentric ring arrays [J]. IEEE transactions on antennas and propagation, 2008, 56(1): 266-268.

[11] 張光義. 相控陣雷達系統[M]. 北京: 國防工業出版社, 1994: 96-97.

[12] 劉恒, 劉波, 謝廣錢, 等. 迭代傅里葉算法用于六邊形稀疏陣列天線[J].電波科學學報, 2015, 30(6): 1235-1240.

LIU H, LIU B, XIE G Q, et al. Iterative Fourier technique applied for Hexagon thinned array [J]. Chinese journal of radio science, 2015, 30(6): 1235-1240. (in Chinese).

Grating lobe suppression in large spacing arrays composed of subarrays based on concentric ring

LI Zhigang WANG Pengyi

(The54thResearchInstituteofCETC,Shijiazhuang050081,China)

Increasing element spacing appropriately and using subarray-based arrays are favorable to reduce the production cost of phased arrays, however, the rarefied array has unexpected grating lobes. To solve the proplems, an array is designed with no grating lobe. The subarrays are rotated in plane, and are distributed in the concentric ring structure around the center of the array. The range of the ring radii are analyzed for ensuring that each subarray will not overlap. Computational simulation is carried out to optimize the number of the subarrays in each ring and the radii of the ring using an improved method of exhaustion. The results show that, taking 4×4 subarrays and 1 024 elements as example, the method provides suppression of peak side-lobe level of -17.46dB.

grating lobe suppression; large spacing; subarray rotation; concentric ring arrays; sparse arrays

10.13443/j.cjors.2016010101

2016-01-01

中國電子科技集團公司第五十四研究所發展基金(XX151060006)

TN957.2

A

1005-0388(2016)06-1132-06

李志剛 (1982-),男,河南人,博士,中國電子科技集團第五十四研究所高級工程師，主要研究方向為航天測控、陣列天線技術．

李志剛,王鵬毅. 基于同心圓環的大間距子陣級陣列柵瓣抑制[J]. 電波科學學報,2016,31(6):1132-1137.

LI Z G, WANG P Y. Grating lobe suppression in large spacing arrays composed of subarrays based on concentric ring[J]. Chinese journal of radio science,2016,31(6):1132-1137. (in Chinese). DOI: 10.13443/j.cjors.2016010101

聯系人： 李志剛 E-mail: lizhigang0372@126.com

DOI 10.13443/j.cjors.2016010101